Angular spectrum method 2026-11-05
The angular spectrum representation expresses a monochromatic wave as a superposition of plane waves. In a homogeneous, lossless, source-free half-space with real wavenumber , transverse Fourier transform components propagate independently, providing a useful solution of the Helmholtz equation.[ref]
Use angular frequency and time dependence , and define . Then
The outgoing branch of the angular spectrum has nonnegative real and imaginary parts of for propagation toward . A propagating angular-spectrum component has real ; an evanescent angular-spectrum component has and decays as . The source-plane amplitudes and their complex arguments matter, not just their squared magnitudes.[ref]
For an original example, a transverse period gives , and . After distance , its amplitude is multiplied by . Such fine spatial detail does not propagate with unchanged amplitude to the far field.
The paraxial angular-spectrum propagator uses . For the source-plane Gaussian beam , its transverse spectrum is proportional to . The on-axis Gaussian integral therefore gives
The denominator produces the axial Gouy phase . This connects the phase shift with interference of transverse spectral components, rather than an abrupt phase jump at the focus.[ref]
The scalar propagation formula does not authorize arbitrary independent vector components. A physical electromagnetic field must satisfy the transverse polarization constraint in an angular spectrum; in a uniform source-free medium,. Boundary conditions and source consistency accompany a vector reconstruction.
Aurèle J.L. Adam Created 2026-10-18 Updated 2026-11-05
An optical physicist who coauthored an interactive teaching text on wave optics. Its polarization and coherence chapters relate field components and correlations to measurable intensity, helping explain why combining beams does not always produce stable interference fringes.[ref]
An evanescent angular-spectrum component has transverse spatial frequency greater than the medium wavenumber. In a homogeneous lossless medium, with , its forward longitudinal wavenumber is imaginary:
It oscillates transversely while decaying away from the source plane. The bounded sign is fixed by the outgoing branch of the angular spectrum.[ref]
For an original transverse period , take , . Then and the amplitude-decay distance is . At , the amplitude is about 0.0658 of its source value, while its squared magnitude is about 0.00433. Amplitude and squared-amplitude attenuation must not be interchanged.
This explains a basic limitation of far-field reconstruction: fine transverse structure can be present at a boundary yet be strongly suppressed before reaching a distant detector. Inverting the factor also amplifies errors. A tiny measured high-frequency component may therefore be an unstable basis for reconstructing fine source detail.
The electromagnetic field is still constrained by Maxwell equations. In particular, remains valid for complex . A single ideal evanescent wave in a lossless medium does not carry time-averaged energy flux normally away from its interface in the same way as a propagating wave. Near-field coupling or a superposition involving another boundary can nevertheless transfer energy, so decay is not a universal prohibition on energy exchange.
An illustrative experiment compares the same patterned source with detectors at increasing distances and accounts for the detector's own spatial response. The decay predicted here concerns a resolved spectral component under the stated homogeneous-medium model, not a guarantee that every observed loss of contrast is evanescence. Absorption, finite aperture and instrument resolution can also suppress structure.
The evanescent field in total internal reflection oscillates and advances in wave phase along the boundary while its wave amplitude decays away from it. With the coordinates above and time dependence , a TE polarization has complex amplitude
The physical electric field is . The decaying sign excludes a field that grows without bound at . Dielectric-interface boundary conditions determine ; for nongrazing incidence it is generally nonzero.[ref]
Zero normal energy flux follows directly from Maxwell equations. The wavevector is ; the relation gives
Consequently, the cycle-averaged Poynting vector is
The electric field and are a quarter-cycle apart, so their product has zero cycle average. Tangential energy flux survives. This is a local calculation for the single decaying evanescent wave in the infinite second half-space.
To check reflected power independently, put . The existing transverse-electric Fresnel reflection coefficient becomes . Its numerator and denominator are complex conjugates, so , while supplies the surface electric field relative to the incident wave amplitude. Thus a unit optical reflectance does not require a zero transmitted field. The same normal power balance holds for TM polarization; its reflected wave phase generally differs from the TE polarization phase. A mixture can therefore change its polarization of an electromagnetic wave on reflection even though all normal incident power returns.
H. Paul Urbach Created 2026-10-18 Updated 2026-11-05
An optical physicist who coauthored a teaching text connecting electromagnetic wave descriptions with optical measurement. Its accounts of polarization and interference make the coherence assumptions behind an intensity prediction explicit, an essential step when applying ideal wave models to real sources.[ref]
Chapter of Interactive Optics by Sander Konijnenberg, Aurèle J.L. Adam and H. Paul Urbach. Whole-work author credits apply; no separate chapter author byline is established. No chapter publication or revision date is established; the URL path does not establish a date or edition.
Cited passages: Scalar Diffraction Optics: Fresnel Approximation, equations (381)–(385), following paraxial equation remark; Fresnel and Fraunhofer Approximations, equations (386)–(392): omitted aperture quadratic phase, full-diameter Fresnel number; About the authors: biographies and final whole-book author credit.
The outgoing branch of the angular spectrum selects the longitudinal wavenumber that represents forward radiation and bounded evanescent decay in a specified source-free half-space. For real , time dependence and propagation into ,
The first root is positive real and the second positive imaginary. They accompany the factor .[ref]
The scalar Helmholtz equation permits both signs. For , and describe opposite longitudinal phase propagation under the stated time convention. For , the corresponding factors are decaying and growing exponentials. A source-free upper half-space with no incident field from infinity uses the forward, bounded choice.
For an original transverse spatial frequency , the decaying root is . At , amplitude changes by . Choosing the other sign would instead multiply it by approximately 111.3. A sign error can therefore become an enormous near-field propagation error without changing the input data.
This choice is a statement about the physical region and boundary conditions, not an instruction to delete reflected waves from every problem. An interface or cavity can require both propagation directions in a finite region. Likewise, propagation towards negative must use a convention consistent with that region. At , the root is zero and the component is grazing; formulas that divide by need separate attention.
Backward reconstruction of a decayed evanescent angular-spectrum component multiplies measurement error by an increasing exponential. That is an inverse problem's conditioning issue, rather than evidence that the forward outgoing branch should grow. Stating the time convention, physical half-space and radiation condition prevents the square-root symbol from hiding these choices.
The paraxial angular-spectrum propagator is the quadratic-phase approximation to forward free propagation of transverse Fourier transform components concentrated near the optical axis. For real medium wavenumber , write . With time convention , the exact factor is . Expanding at gives
after removing the common carrier .[ref]
The approximation preserves the magnitude of each retained spectral component and changes its wave phase. Recombining those components produces diffraction; it need not preserve the shape or intensity at an individual spatial point. In coordinate space the envelope satisfies the parabolic wave equation, .
Small propagation angle and acceptable accumulated phase error are separate conditions. The first omitted term produces a phase discrepancy of order . For an original component with , the exact difference between approximate and exact phases is
It is about 0.00789 radians at , but about 0.789 radians at . The direction is still only about off axis. Propagation distance matters when phase accuracy is the observable.
For a source-plane Gaussian beam , the spectral width is of order . Thus concentrates the important components in the paraxial region. That provides a reason for the familiar wide-waist condition rather than treating all Gaussian formulas as exact electromagnetic solutions.
A calculation containing substantial evanescent angular-spectrum components or large-angle components requires the exact branch and, when appropriate, the vector transverse polarization constraint in an angular spectrum. The paraxial transfer factor has unit magnitude even at large , so extending it blindly into the evanescent region replaces physical decay by an artificial phase rotation.
Fresnel diffraction is interference during propagation when curvature of the optical path across an aperture matters. It provides a continuous account of fringes near geometric edges and their evolution into the far-field pattern. A screen that is many aperture widths away is not automatically in the Fraunhofer diffraction limit; wavelength also enters the comparison.
For plane-wave illumination of a transmission at , the Huygens-Fresnel principle in the scalar paraxial approximation gives
It follows by expanding the propagation distance as while treating the amplitude prefactor as slowly varying. Small propagation angles justify dropping higher path terms. Finite coherent illumination, vector edge fields and very small apertures can require a more complete theory. The Fresnel number , with a clearly declared aperture scale , measures whether the aperture's quadratic phase is significant. When it is negligible, the remaining linear phase produces the aperture Fourier transform.[ref]
Define the optical Fresnel integrals
Their limiting values are , with odd continuation to negative arguments. The Cornu spiral plots : the field contributed by a finite interval is proportional to the chord between its endpoints. Intensities require squared chord lengths; summing zone intensities would lose the cancellation mechanism.
For an opaque straight edge transmitting and all , let , so is the geometrically illuminated side. Integrating in the unrestricted transverse direction and using the half-line in the other gives
hence
Far into the shadow, ; far into the illuminated region, . At the geometric shadow boundary, , not half the incident intensity: the field is halved there. On the illuminated side, alternating constructive and destructive contributions give bright and dark fringes whose contrast declines away from the edge. Their positions scale with and are not equally spaced. A finite detector averages fine fringes and tends to recover the geometric-optics step. In a controlled experiment, translate a straight edge through a broad collimated beam, scan a small photodetector across the screen and normalize against an unobstructed reference; vary to test the square-root scaling. Source size, spectral bandwidth, edge roughness and detector aperture wash out contrast and must be distinguished from a failure of the propagation law.
For a rectangular aperture of full widths , define
The two transverse integrals separate:
and
A long slit leaves only one varying Fresnel factor. For a square centred on axis, and give
Opposite edge contributions overlap; changing aperture size or propagation distance changes their relative phase. Consequently a geometrically illuminated centre can brighten or darken, rather than remaining a fixed top-hat plateau. The original rectangular-aperture map shows fringes parallel to both edges and their crossing structure. Outside the geometric rectangle the wave also extends into the shadow.
At large Fresnel number an interior region approaches geometric illumination while edge fringes remain on a scale . At intermediate values the opposite-edge fringes interact across the whole aperture image. In the far field, the formula tends to the product of sinc-squared factors,
The far-field lobe scale grows linearly with , unlike the local near-edge scale. The transition is controlled by the neglected phase, not an abrupt physical boundary.
To test the main effect, record normalized profiles at several distances, calculate the Fresnel integral with the measured aperture dimensions and compare fringe positions before fitting a free intensity scale. Use the same transverse coordinate origin and check that the incident beam is sufficiently uniform. Babinet's principle relates complementary aperture and obstruction fields through ; their near-field intensities are not generally equal, because the unobstructed field remains present and interferes with the diffracted field. This is why a far-field obstruction argument cannot simply be carried over as equality of Fresnel intensities.
Figure 1.
Original Fresnel-integral calculations for a straight edge, slits, rectangular aperture and Cornu spiral
. The plots use the stated scalar paraxial model and chosen dimensionless aperture parameters. The straight-edge boundary value is one quarter. The aperture panels show quadratic-phase propagation and the Cornu panel identifies field amplitudes as chords rather than summed intensities. Details.
A propagating angular-spectrum component is a transverse Fourier transform component whose longitudinal wavenumber is real in a homogeneous lossless medium. With and the forward outgoing branch of the angular spectrum,
Its propagation factor has unit magnitude and a distance-dependent wave phase.[ref]
For an original component with and , and . It repeats its phase along over distance , even though its wavelength measured along its propagation direction remains . The axial phase period is not an energy-transit time.
Preserving each component's magnitude does not preserve the reconstructed spatial intensity. Different values of accumulate different phases; their coherent superposition changes the interference pattern and can spread a beam. This is the angular-spectrum description of diffraction.
Energy crossing a reference plane also depends on direction. For a vacuum plane electromagnetic wave with electric-field magnitude , the time-averaged Poynting vector has magnitude along the propagation direction. Its normal component through a plane of constant is that magnitude times , where is the vacuum impedance. Equal electric-field magnitudes at angles zero and therefore give normal fluxes in the ratio .
The boundary case is grazing, with and no forward normal phase advance. For , the root becomes imaginary and the component is instead an evanescent angular-spectrum component. These classifications apply to the stated real-wavenumber, lossless model; material absorption requires a consistent complex-wavevector treatment.
Sander Konijnenberg Created 2026-10-18 Updated 2026-11-05
A physicist who coauthored Interactive Optics by Sander Konijnenberg, Aurèle J.L. Adam and H. Paul Urbach. Its treatment connects wave polarization, interference and coherence to optical experiments. Writing those topics in a common wave framework helps readers move from a field equation to a prediction for the intensity measured by an instrument.[ref]
Total internal reflection Created 2026-09-29 Updated 2026-10-18
Total internal reflection occurs when a wave travels from higher to lower refractive index above the critical angle for total internal reflection, leaving no propagating transmitted plane electromagnetic wave normal to the interface.
At an infinite planar interface, total internal reflection returns all incident normal power in the reflected plane electromagnetic wave. Assume homogeneous dielectrics with isotropy, linearity, no dielectric loss, magnetic permeability , and positive real refractive indices . Above the critical angle for total internal reflection, both TE polarization and TM polarization have unit optical reflectance despite the external evanescent wave.[ref]
This makes total internal reflection useful for guiding light in an optical fiber: suitable internal incidence angles keep energy inside successive boundaries. Its more surprising feature is that perfect optical reflectance coexists with an external electric field. The evanescent field in total internal reflection can interact with a nearby object, as the frustrated total internal reflection experiment demonstrates.
For an original idealized example, choose , at vacuum wavelength . At angle of incidence , Snell law for electromagnetic waves gives an angle of refraction , so a transmitted light ray still emerges. At , its direction becomes tangent to the boundary. At , the proposed sine of a real angle of refraction would be : the transmitted wavevector must instead have an imaginary normal component. These are predictions for the stipulated media, not measured properties of a particular block.
Figure 1.
Constructed ray geometry just below, at, and above the critical angle for refractive indices 1.5 and 1; arrows indicate directions rather than relative powers
. Details.
A real semicircular-block experiment isolates this threshold. The University of Maryland laboratory instructions place a laser, motor-driven rotation stage and photodetector together; an aperture limits the detector voltage. The semicircular block is centered on the rotation axis, with the photodetector aligned along its flat face. The instructions identify a strong grazing-beam signal and a weaker signal just before the critical angle for total internal reflection. The readout is detector voltage versus stage angle. Retracing the incident laser beam identifies normal incidence for alignment.[ref]
The semicircular shape supplies a geometric control: a light ray passing radially through the curved face enters at normal incidence, avoiding additional entrance refraction. The weaker signal near the flat face locates a limiting angle of refraction; a vanishing remote-detector signal alone would not establish that the external electric field vanishes. Beam width and alignment limit angular resolution; no measured threshold or decay depth is supplied.[ref]
The transverse polarization constraint in an angular spectrum requires each source-free electromagnetic plane wave component to have electric field orthogonal to its wavevector in the algebraic sense:
For an isotropic homogeneous medium with constant permittivity, this follows by inserting a Fourier component into . It restricts the vector fields that can accompany a scalar angular spectrum method calculation.[ref]
With time dependence and , the relation becomes
When , the longitudinal field is therefore . It is not an arbitrary third scalar image that can be propagated independently of the first two.
For an original propagating angular-spectrum component with , choosing and forces . Polarization along instead gives a purely transverse-to-axis field for this propagation plane. Thus electromagnetic transversality relative to does not require zero field along the laboratory axis.
For an evanescent angular-spectrum component with , the same choice gives . The dot product is not a Hermitian inner product: there is no complex conjugation of in the Fourier-transformed divergence equation. The imaginary longitudinal amplitude represents a phase relationship, not permission to discard the constraint.
At grazing propagation, , division by fails and the original equation must be used. Boundary conditions, material response and the corresponding magnetic field remain necessary to reconstruct a physical solution. In the strongly paraxial approximation, longitudinal components may be small, but that is an approximation derived from the allowed spectrum rather than a replacement for the vector condition.