gaussian_beam — OPTICS wave op

• Data kinds: कोई नहीं → table (केवल arguments से तय होने वाला op — यह कोई image या data input नहीं लेता)

• Call: import fullseye as fs; fs.ledger.gaussian_beam(waist_um=100.0, wavelength_um=1.064, distance_mm=0.0, n_medium=1.0) (implementation को सीधे call करना हो तो import optics; optics.gaussian_beam(waist_um=100.0, wavelength_um=1.064, distance_mm=0.0, n_medium=1.0), registry से लेना हो तो opsoptics.get("gaussian_beam"))

उपयोग

Gaussian-beam propagation: spot size, wavefront curvature और Gouy phase।

> नीचे का विस्तृत विवरण मूल text में है — सारांश और headings का अनुवाद हो चुका है।

The q-parameter solution of the paraxial wave equation, evaluated at *distance_mm* from the waist. With zR = pi*w0^2*n/lambda:

w(z) = w0*sqrt(1 + (z/zR)^2) · R(z) = z*(1 + (zR/z)^2) · psi(z) = arctan(z/zR) (Gouy) · far-field half-divergence theta = lambda/(pi*n*w0).

Returns a dict: radius_um the 1/e^2 field radius w(z) · rayleigh_mm zR · wavefront_radius_mm R(z) · curvature_per_mm 1/R(z) · gouy_deg · divergence_mrad the far-field half angle · waist_um and distance_mm echoed back.

wavefront_radius_mm is inf exactly at the waist, by contract:

the wavefront there is planar and a plane has infinite radius. That is why curvature_per_mm is also returned — it is 0 at the waist and finite everywhere, so downstream arithmetic never has to divide by an infinity. A negative R (before the waist, z < 0) means a converging wavefront.

Ground truth it reproduces exactly (verified at machine precision): at z = zR the spot is sqrt(2)*w0, the wavefront radius is 2*zR (its minimum over all z) and the Gouy phase is exactly 45 degrees; the beam-parameter product w0*theta = lambda/(pi*n) holds for every input.

Raises ValueError: non-positive or non-finite *waist_um*, *wavelength_um*, *n_medium*; non-finite *distance_mm*.

Paraxial and fundamental-mode (TEM00) only: a real multimode beam needs an M^2 factor, which this op deliberately does not fake by accepting one silently.

पूरी family का input contract (fail-closed)

optics के सभी ops गणना से पहले input की जाँच करते हैं (कुछ भी चुपचाप पास नहीं होता):

• Unit argument के नाम में ही लिखी होती है — _mm / _um / _deg / _mrad। mm और µm की गड़बड़ी crash नहीं करती, बल्कि "सही दिखने वाला पर गलत उत्तर" देती है, इसलिए नाम से इसे रोका जाता है। मान के आकार से unit का अनुमान बिल्कुल नहीं लगाया जाता।

• String से ValueError — float('50') सफल हो जाता है, इसलिए बिना parse किया config value length के रूप में निकल जाता (मापा गया: thin_lens('50', '200') सही दिखने वाला 66.667 mm लौटा रहा था)। bool को भी True == 1 वाले implicit promotion के रूप में अस्वीकार किया जाता है।

• complex / masked array से ValueError (केवल real-valued slots; imaginary part को चुपचाप छोड़ना या mask हटाना अस्वीकार)। हर input पर NaN/Inf से ValueError।

• शून्य से भाग और उससे जुड़े मामलों को नाम लेकर अस्वीकार किया जाता है: focal length 0, radius of curvature 0, refractive index <= 0, पूरी तरह opaque aperture (सब 0, इसलिए normalisation 0/0), ऐसा PSF जिसका योग <= 0 हो, S0 = 0 वाला Stokes vector, और front focal point पर रखा object (image अनंत पर)।

• केवल 2 ops non-finite value लौटाते हैं, और वह भी contract के रूप में लिखा हुआ: depth_of_field hyperfocal distance से आगे far_mm = inf (hyperfocal distance की परिभाषा यही है) और gaussian_beam waist पर wavefront_radius_mm = inf (plane wavefront की radius of curvature)। दोनों साथ में एक finite साथी (far_is_infinite / curvature_per_mm) भी लौटाते हैं। इसके अलावा कोई भी चुपचाप आया NaN/Inf अंदर ही पकड़ा जाता है और ValueError देता है — "float64 overflow हुआ" और "उत्तर अनंत है" अलग-अलग दावे हैं, इसलिए पहले को दूसरे का चेहरा पहनाकर नहीं लौटाया जाता।

• Size की सीमा: generate होने वाले grid के लिए optics.MAX_GRID (4096), दिए गए field/PSF/aperture के लिए optics.MAX_FIELD_ELEMENTS (2^24), ABCD element chain के लिए optics.MAX_SYSTEM_ELEMENTS (1024), Zernike के लिए MAX_ZERNIKE_TERMS (512) / MAX_ZERNIKE_ORDER (40) / MAX_ZERNIKE_BASIS (2^25)। छोटे argument से बहुत बड़ा internal allocation होने वाले रास्तों (मापा गया: n_max=40 × 4096² पर 108 GB) को fail-closed ढंग से बंद किया गया है।

• भौतिक रूप से असंभव states भी अस्वीकार: degree of polarisation > 1 वाला Stokes vector, ऋणात्मक transmittance, ऋणात्मक intensity, और n-|m| odd जैसे अमान्य Zernike indices।

उपयोग की विस्तृत guide

• optics_imaging family guide

संदर्भ (sample data, साहित्य)

• Sample data catalog (DL URL / license) — 2-D के लिए skimage.data (BSD/public) + synthetic, 3-D के लिए real data sources (Stanford/PDS आदि) के DL URL।

• Operators का स्रोत और संदर्भ — इस op family के मूल research/methods के स्रोत।

• Algorithm का मानक स्रोत (लेखक, वर्ष) और उपयोग ऊपर की family usage guide में दिए गए हैं।

चलाए जा सकने वाले उदाहरण (verified samples जो सच में यह op call करते हैं)

• optics_imaging — py -3.11 examples/optics_imaging.py

वे अगले ops जिनसे type जुड़ता है (table को input के रूप में लेते हैं)

abcd_matrix · wavefront_stats · sfr_from_edge · paraxial_trace · seidel_coefficients · spot_stats · tolerance_analysis · wavefront_from_opd

इसी category के (wave)

airy_pattern · angular_spectrum_propagate · fraunhofer_pattern · fourier_plane_filter · four_f_filter · defocus_from_shift · pupil_psf · pupil_blur


*Provenance: optics.py — OPTICS operator registry. यह per-op note tools/opdocs.py md से अपने-आप generate होता है (हाथ से edit न करें)।*

© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.