Curved boundary integral method for electromagnetic fields

The angular spectrum method is a rigorous method to synthesize near and far-field electromagnetic beams from planar field distributions. However, this limitation of planar surfaces has restricted its applicability to beams with simple focal planes. We propose a curved boundary integral method (CBIM) to synthesize electromagnetic beams from arbitrary surfaces to address this limitation and expand the method's scope to synthesize beams from and between shaped objects. This study presents a detailed theoretical framework behind the CBIM and validates its effectiveness and accuracy with a comprehensive set of simulations. Additionally, we present mathematical proof to support our proposal. The proposed method satisfies Maxwell's equations and significantly benefits optical systems and inverse beam design. It allows for analyzing electromagnetic forward/backward propagation between optical elements using a single method. It is also valuable for optical force beam design and analysis.