At a fixed camera-axis direction, the grating equation gives , so adjacent central wavelengths differ by . This differs from the exact wavelength cell where an order is nearest the optical camera axis. For and , equal distances of orders and mean opposite diffraction angles; adding their diffraction grating equations gives . The full order- cell is bounded by and , hence has width for and accessible neighbouring orders. Both widths approach at high order.
For a circular illuminated aperture of diameter , the characteristic angular diffraction scale is . The first zero of the ideal intensity pattern lies at approximately .
Geometrical optics Created 2026-10-05 Updated 2026-10-06
Geometrical optics describes propagation by rays and Snell's law, neglecting diffraction on scales much larger than the wavelength.
Put
In free space, the parabolic wave equation is . Under the Fourier transform convention , it becomes
The Gaussian integral is legitimate because . Invert the transform after multiplying by . A second Gaussian integral, or equivalently one-dimensional transverse Fresnel propagation, gives
Choose the square-root branch continuously from ; for real there is no zero of . This gives the correct incident field at , and direct differentiation verifies the free parabolic wave equation.
For clarity, the squared envelope magnitude is
The Gaussian beam with one transverse coordinate remains Gaussian, with one-transverse-coordinate amplitude factor , not the of a beam with two transverse coordinates. The negative initial quadratic phase produces focusing for ; diffraction prevents a singularity at . These expressions describe the paraxial approximation to free propagation, rather than an exact unrestricted Helmholtz equation beam.
Write and . Substituting into part (a), without prematurely replacing by its linear part, gives
Taking the expectation yields the exact mean equation within the paraxial approximation:
The correlation term cannot be discarded just because : the field depends on the same random medium as . Nor can be factored exactly from the given variance alone. The printed statistics do not specify the autocorrelation function of a random field, so they do not determine a unique local closed equation for .
A useful covariance-dependent equation follows to second order in . Define
and let
be the Fresnel propagator kernel for one transverse coordinate. The Duhamel principle applied to the linear fluctuation term gives
where . Multiply by and average. Also . Since , replacing by in terms already multiplied by preserves second-order accuracy. We obtain the weak-fluctuation memory equation for the coherent field
This perturbation statement is for fixed propagation intervals with sufficient covariance regularity and moment bounds. The kernel is understood as an oscillatory integral. If the field is jointly Gaussian, odd perturbative contributions vanish and the remainder can be improved under the corresponding expansion hypotheses. Joint Gaussianity itself still leaves the covariance unspecified. The term is the quadratic refractive-index shift in the coherent field; it comes from the mean of .
For the usual local answer, impose the additional Markov approximation for a random medium. Its longitudinal correlation length must be much shorter than the distance on which the envelope varies, and diffraction across that correlation distance must not resolve substantial transverse covariance variation. Away from the entry layer, replace the slowly varying factors in the memory integral by their local values and extend the longitudinal integral to infinity. Define the integrated covariance
For an integrable real stationary covariance, is even, and is its zero-longitudinal-frequency spectral density. The conventional leading weak-index model keeps only the linear random potential ; in that model the local first-moment equation is
The factor can be checked independently in the white-noise model. If is a Brownian field with , write
Converting the Stratonovich integral to an Itô integral adds . The remaining stochastic integral has zero mean, giving the boxed equation and
If the finite-correlation model retains consistently through order , the local second-order expansion also has , adding the phase factor . This phase is computed before the white-noise idealization: squaring ideal Gaussian white noise is not the finite-variance operation .
The need for a correlation assumption can be seen without any closure argument. The stationary Gaussian random field , with a single standard normal variable , has all the printed one-point statistics. In the linear weak-index model it gives
Its attenuation is quadratic in propagation distance, whereas the Markov model gives linear-distance exponential attenuation. The local attenuation equation requires the additional correlation/Markov assumption; it does not follow from statistical homogeneity and unit variance alone.
The optical telescope focuses the sky onto an optical slit; the collimator turns the selected light into a beam incident on the reflection diffraction grating; the optical camera images each diffracted direction onto the photodetector. The drawing also shows the monochromatic optical slit image for uniform optical slit illumination in the geometric, slit-limited approximation.
Figure 1.
Reflection-grating spectrograph and slit-limited line profile
.
Near the chosen wavelength , a small wavelength increment moves the image by , where is the magnitude of the linear spectral dispersion. An optical slit image of physical width therefore has an apparent spectral width . With one optical slit width as the adopted separation criterion,
Uniform illumination gives a top-hat monochromatic profile of width ; an unresolved stellar image need not illuminate the optical slit uniformly, so its profile follows the actual illumination convolved with the instrument response. Negligible optical slit diffraction and optical camera optical aberrations do not remove the finite diffraction grating diffraction width. The formula is the slit-limited result when that width is small compared with , rather than a universal exact line-profile criterion.
Substitution of gives
The incident and outgoing optical path length contributions between the two ends of the illuminated diffraction grating add. Their total is , so is the edge-to-edge optical path difference in the selected diffracted direction. This is a length, not a dimensionless resolving power. With illuminated grooves, and the ideal diffraction grating diffraction-limited spectral resolving power is . A wide optical slit generally gives a smaller slit-limited ; this distinction limits how far the slit-width scaling may be extrapolated. All of these formulae assume the stated reflection-grating sign convention.
The fundamental angular floor is the diffraction limit of a telescope, of order . The ideal diffraction grating has , whereas its slit-limited resolving power has . Requiring gives
The coefficient depends on the adopted line-resolution and aperture conventions. Below this scale, a narrower spaxel oversamples the same spatial mode rather than creating a new independently resolved element. Equivalently a spatial optical mode occupies etendue of order .
There is also an engineering floor before arbitrary shrinkage: the optical camera focal length, photodetector size and support clearances do not shrink, while the number of photodetectors increases as . Thus the shrinking-beam model's constant combined optical camera volume is eventually overtaken by fixed per-camera overheads. For fixed brightness, smaller spaxels receive fewer photons, so read noise and photon shot noise can constrain useful sampling before the formal optical floor. Diffraction limits independent spatial information; fixed photodetector and mechanical dimensions limit practical volume reduction.
Use the time-harmonic wave convention and write , , with a deterministic incident wave envelope. The printed speed ratio is inconsistent with : the refractive index used below is . Substituting in the Helmholtz equation gives
The paraxial approximation discards . The usual weak-fluctuation model also linearizes , giving
Dropping the quadratic contrast is an additional weak-fluctuation assumption, not a consequence of small propagation angles. Although the linearized random potential generates attenuation of order , it does not retain every effect of that order in the literal finite-correlation index: the discarded quadratic contrast can also produce a mean wave phase shift. This linearization must precede a Gaussian white noise limit: the square of ideal white noise has no ordinary pointwise meaning.
The split-step Fourier method alternates free-space diffraction, , with a random phase screen,
For jointly Gaussian random fields, Gaussian phase averaging gives . More generally, a product of fields and conjugate fields picks up , with signs or . Its screen average is determined entirely by the wave phase covariance matrix. The deterministic diffraction step acts on each coordinate, with opposite signs on factors formed by complex conjugation. This is the basis of the field-moment equations.
There is an important closure qualification. A stationary Gaussian random field need not have independent longitudinal increments. For finite-correlation fluctuations, the unlinearized parabolic wave equation instead gives
In the weak-fluctuation model, the exact first-moment equation is
The last term cannot in general be replaced by a constant times . An exact generic solution of the linearized model is
where denotes time ordering of the propagation operators. For the unlinearized model, add inside the propagation generator. The covariance function of the medium is needed to evaluate this expression; One-point Gaussian distributions alone would not even determine the joint wave phase statistics.
For example, omit diffraction and take a longitudinal autocovariance function . Direct Gaussian phase averaging gives
Even here, unit-variance stationary Gaussian random fields with different give different answers. In the general problem the diffraction and multiplication operators do not commute, so this scalar attenuation cannot simply be multiplied by without an additional approximation.
The standard closed answer uses the Markov approximation for a random medium, made explicit in part ii. Let the longitudinally integrated covariance kernel be
Replace the medium by longitudinal Gaussian white noise with this strength. A screen of thickness then has and is independent of the incoming field. For the moment , expanding both steps to order gives
In particular, the coherent attenuation in a white-noise random medium and its Fresnel propagator solution are
Equivalently, write , where is Brownian motion in with transverse covariance kernel . The Stratonovich integral formulation is . Its Itô integral form is
The mean of the Itô integral vanishes, independently confirming the attenuation drift. For in two transverse dimensions,
Thus a unit plane wave has . For a general incident wave envelope, the Fresnel propagator supplies its spreading. The attenuation is redistribution between the coherent and diffuse wave fields, rather than wave absorption: for a field and its conjugate at the same point, the screen contribution in the second-moment equation cancels.
In a classical Cassegrain reflector, a concave parabolic primary sends light toward an intermediate focus, but a convex hyperbolic secondary intercepts it before it reaches that focus. The secondary returns the beam through a central hole in the primary to a focal plane behind the primary. The two relevant foci of the secondary’s hyperbola are the primary’s would-be focus and the final focus.
The secondary magnifies the effective focal length, giving a compact tube and convenient rear-mounted instruments. The classical conics correct on-axis spherical aberration, but not the off-axis coma, astigmatism or field curvature. There is secondary obscuration, diffraction from its supports, and sensitivity to mirror alignment. A long effective focal length is useful for a small angular image scale per detector pixel, but yields a small field for a fixed detector size.
Parabolic primary → convex hyperbolic secondary before prime focus → rear focus.
A Newtonian reflector uses a concave primary generated by a parabola, followed by a flat secondary inclined at to the optical axis. Parallel marginal rays converge toward the primary focus; the secondary intercepts that converging beam and folds it sideways to an accessible focal plane.
The Newtonian reflector needs only one powered mirror, so it is relatively simple to fabricate and inexpensive. A parabolic primary has no on-axis spherical aberration. Its principal wide-field limitation is coma, accompanied by field curvature and off-axis astigmatism. The diagonal and its supports obstruct the entrance pupil and introduce diffraction; the tube is long compared with a folded two-powered-mirror design, and heavy instruments at the side focus can be awkward to support. These comments also answer the unheaded design-comparison clause on the next PDF page.
Parabolic primary → flat diagonal → side focus.
The telescope forms the sky image at the entrance slit. A collimator makes the transmitted beam parallel, a reflection diffraction grating disperses it, and a camera’s optical lens focuses each wavelength to a different detector position.
Figure 1.
Reflection-grating spectrograph and slit-limited line profile
.
With a uniformly illuminated slit, ideal imaging and negligible slit diffraction, each monochromatic slit image has an approximately top-hat intensity profile of physical width . A small wavelength separation displaces two such images by , using the positive magnitude of grating dispersion. In the geometrical, slit-limited convention, resolution occurs when the displacement is of order the slit-image width. Hence
Here is the slit’s apparent wavelength extent. A finite diffraction grating broadens the sharp edges; the top-hat sketch and the following identities presume slit-limited resolving power of a grating, rather than the regime where the grating’s own diffraction determines the line width. A precise resolution criterion for nonideal profiles can change order-unity factors.
Substitute the grating equation into the preceding result:
Moving across the illuminated grating by one groove spacing changes the incident-plus-outgoing optical path length by . Across length , the total path difference is therefore . Thus is the optical path difference between contributions from the two ends of the illuminated grating.
The number of coherent phase cycles across it is , which is the intrinsic diffraction-limited spectral resolving power of the diffraction grating under the usual first-minimum criterion. In the slit-limited regime, also states how much sky angle can be accepted at a given resolution and aperture. Holding slit angle and resolution fixed while increasing telescope diameter requires a larger optical path span. The geometry gives ; large incidence and diffraction angles increase resolution per unit grating length, though grazing beams become impractical.
The slit equations do not imply unlimited resolution when tends to zero. Once is comparable to , finite-aperture diffraction matters and the actual resolving power is bounded by about .
Let the Bravais lattice have primitive vectors , and put
The reciprocal lattice has the basis
The scalar triple product and cross product identities give
Consequently every satisfies for every lattice point , and conversely these three conditions force the coefficients of in the reciprocal basis to be integers. This proves that the displayed vectors generate .
In the Born approximation, suppose that the crystal potential is the sum of translates of one atomic potential over the finite set of lattice sites. Its Fourier transform at the momentum transfer factors as
Thus the single-atom scattering amplitude is multiplied by the crystal lattice structure factor . Writing
separates the sum into three finite geometric series. For this gives
When the are large, a factor has a sharp maximum when . All three factors are therefore simultaneously large precisely when . These are the reciprocal-lattice peaks of a finite crystal.
For the stated Body-centered cubic lattice basis,
the reciprocal-basis formula gives
Equivalently,
The shortest nonzero reciprocal vectors have squared Euclidean norm , so
In elastic scattering, . If and is the angle between and , Euclidean geometry gives
The first possible diffraction peak therefore occurs at
as , by the small-angle approximation. This is the elastic Bragg scattering condition for the shortest reciprocal-lattice vector.