{"id":"91a7722a-53f5-417b-a696-5daeb515634f","arxiv_id":"2411.14048","paper_version":1,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Focused Gaussian beams acquire spin and orbital angular momentum from wavefront curvature, even when the input beam has zero angular momentum.","lead":"A linearly polarized Gaussian beam develops spin and orbital angular momentum near the focal plane because the wavefront is curved, even though the beam starts with no angular momentum. The result suggests wavefront curvature can be used as a new control knob for optical tweezing and light-matter interactions at the nanoscale.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No significant objection identified: the main analytical claims survive scrutiny.","rationale":"The reader identified Eq. (5) as the weakest assumption, and I agree that this is the natural place to attack the derivation. However, closer order-by-order analysis shows the neglected z-variations do not feed into the specific second-order quantities that carry the central claims: the longitudinal electric spin s_z and the longitudinal orbital density l_z for the linearly polarized Gaussian. The slow envelope derivative is a higher-order effect for s_z and has vanishing angular derivative for l_z. The vectorial diffraction modelling provides independent, if qualitative, confirmation. No internal inconsistency, same-order omission, or circular step was found. The paper's self-declared limitation in the Discussion (neglecting z-variations of the amplitude) is real, but it is consistent with an O(ε³) error that does not alter the O(ε²) results highlighted here. Therefore the ACCEPT verdict stands. I propose a concrete symbolic/quantitative check as a prudent verification step, not because the concern is currently persuasive, but because the lengthy algebra is otherwise unverified and the numerical comparison is only qualitative.","tokens_in":14824,"tokens_out":31180,"duration_ms":307907,"concrete_test":"Recompute Eq. (6) with a symbolic Maxwell solver that retains ∂u/∂z = (i/2k)∇⊥²u in the Faraday and Ampere iteration steps, and verify that the resulting s_z and l_z for a linearly polarized fundamental Gaussian (α=1, ℓ=p=0) coincide with Eqs. (13) and (20). In addition, evaluate the vectorial diffraction code at a single point, e.g., r = w0/√2, φ = π/4, z = z_R, with w0 = λ and a stated NA (say NA = 0.5), and compare the value of s_z with Eq. (13); agreement within about 10% would close the residual concern about the truncated expansion.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most plausible concern is the reader's weakest assumption: Eq. (5) replaces the z-derivative of the paraxial envelope by ik, neglecting envelope and Gouy-phase variations. If this omission entered at the same order as the claimed second-order spin/orbital densities, the central attribution to wavefront curvature would be compromised. Order counting shows it does not. The neglected envelope derivative is ∂u/∂z = (i/2k)∇⊥²u (from the paraxial equation), which is O(ε²) relative to iku, where ε = 1/(kw0). Thus corrections to L1 are O(ε³), and corrections to T2 are O(ε²) but only in certain components. For the linearly polarized Gaussian case (α=1, ℓ=p=0): s_z (Eq. 13) is built from Im(u* T2_y); T2_y is obtained from ∂_x B_z, with B_z = (1/iω)(-∂_y E_x), and is free of slow-z corrections because B_x=0 enters the Ampere step. l_z (Eq. 20) is built from Im(u* ∂_ϕ T2_x); the slow correction to T2_x is proportional to ∇⊥²u, which is ϕ-independent for a fundamental Gaussian, so its ∂_ϕ vanishes. Hence the reported second-order densities are unaffected at the quoted order. The vectorial diffraction comparison (Figs. 1 vs 4, S4-S5) independently supports the qualitative picture. The only residual weakness is that the comparison is visual rather than quantitative, but that does not rise to a load-bearing objection to the central claim.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript constructs a systematic nonparaxial description of focused Laguerre-Gaussian beams by iteratively applying Maxwell's equations, starting from the paraxial field T0 and generating first-order longitudinal (L1) and second-order transverse (T2) corrections. The authors then evaluate electric spin and orbital angular momentum densities to second order in the paraxial parameter. The central claims are (i) a 2D linearly polarized Gaussian beam acquires a longitudinal electric spin density s_z = -2ε0 r^2/(k w^2 R[z]) sin 2φ |u|^2 (Eq. 13) and an equal-and-opposite orbital density l_z (Eq. 20) around the focal plane, so that the total longitudinal angular momentum density vanishes; and (ii) circularly polarized beams acquire a helicity-dependent transverse spin component from the wavefront-curvature term kr/R[z] (Eq. 14), in contrast to the helicity-independent transverse spin of evanescent waves. The analytical results are compared with vectorial diffraction theory, with qualitative agreement shown in Figs. 4, S4, and S5. The paper also discusses optical helicity and dual-symmetric spin, noting that the longitudinal spin for the linearly polarized case is electric-only and carries zero helicity.","tokens_in":15110,"tokens_out":9239,"duration_ms":83202,"significance":"If correct, the result is significant because it identifies wavefront curvature as a physical control knob for local optical angular momentum, shows that focusing a simple linearly polarized Gaussian beam produces local spin and orbital densities, and predicts a helicity-dependent transverse spin in focused circularly polarized beams. The derivation is first-principles, contains no fitted parameters, and is independently validated by vectorial diffraction. The explicit analytical formulas (Eqs. 12-16 and 20) are useful and falsifiable, and the open-source numerical implementation is a further strength. The authors are transparent about the dual-symmetric nature of the spin and about the vanishing optical helicity in the linearly polarized case. The main weaknesses are presentation-level: the numerical validation is visual rather than quantitative, and some terminology around 'gradient of wavefront curvature' is imprecise, but neither undermines the central derivation.","major_comments":[],"minor_comments":[{"comment":"The validation against vectorial diffraction theory is presented only as color maps; because Eq. (5) is the key approximation that underpins the analytical field, a quantitative comparison (line profiles or normalized errors) for the spin densities, and ideally also for the orbital density of Eq. (20), would make the claimed 'extremely good match' explicit and reproducible.","section":"Section VI, Figs. 4, S4, S5"},{"comment":"The phrase 'gradient of the wavefront curvature' is imprecise: the derived expressions are proportional to 1/R[z], i.e., to the wavefront curvature itself, not to its derivative along z. Please define the intended meaning, for example 'the transverse gradient of the wavefront phase due to curvature.'","section":"Section II, Eq. (7) and Eqs. (12)-(14)"},{"comment":"The statement that R[z] 'diverges to +∞ at both the focal plane and infinity' is not accurate for negative z, where R[z] = (z^2 + z_R^2)/z is negative; rephrase as 'diverges at the focal plane and grows in magnitude toward infinity.'","section":"Section III, text near Eq. (13)"},{"comment":"The typeset Eq. (11) contains a stray 'n' and some unbalanced braces in the displayed expression, which makes the formula difficult to parse; please correct the equation formatting.","section":"Eq. (11)"},{"comment":"There are typos in the supplementary figure captions ('propgation', 'Corrsponding', 'distsnce', 'ta') that should be corrected before publication.","section":"Supplementary figures S1 and S3"},{"comment":"The statement 'No data were generated or analyzed in the presented research' sits oddly beside the vectorial diffraction calculations used to produce Figs. 4, S4, and S5; please clarify what code or parameters are available for reproducing those simulations.","section":"Data availability statement"}],"recommendation":"minor_revision","confidential_remarks":"For the editor: this is a solid, original contribution that is squarely within the scope of the journal. The analytical construction is transparent and the stress-test checks confirm that the approximation behind Eq. (5) does not enter at the claimed order. The only substantive caveat, which the authors themselves discuss, is that the longitudinal spin in the linearly polarized case is electric-only, with zero helicity and zero dual-symmetric spin; I would encourage the authors to reflect this qualification more prominently in the abstract. The remaining issues are local and should be resolved with minor revisions."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The headline: I think this paper is right, and the central result is more interesting than the modest packaging suggests. For a linearly polarized Gaussian beam, the authors derive second-order nonparaxial corrections and show the wavefront curvature gradient produces longitudinal spin and orbital densities with s_z = -l_z, so total j_z vanishes. That is a clean result. They also find a helicity-dependent transverse spin for circularly polarized beams, which is a nice contrast to evanescent-wave transverse spin.\n\nWhat is actually new: prior work had longitudinal spin in focused linearly polarized vortex beams; this paper shows the same mechanism operates for an ordinary Gaussian with no OAM, and traces it to the gradient of R[z]. To my knowledge that identification is new. The derivation is systematic: start with paraxial LG, impose Gauss's law to get E_z, use Faraday and Ampère to get B and E to second order, then evaluate spin/orbital densities. No parameters are fitted. The j_z = 0 cancellation is a good consistency check.\n\nSoft spots: the main approximation is Eq. (5), replacing the z-derivative by ik. The stress-test note argues by order counting that the neglected slow-z terms don't affect the second-order densities at the quoted order. I checked the logic and it seems right: for the fundamental Gaussian the dangerous correction is φ-independent and drops out of l_z. That said, the validation against vectorial diffraction is visual, not quantitative; I would like to see a numerical error plot. The algebra is long and I haven't verified every component. No experiment, of course, but the paper is upfront about that, and the analytical derivation plus the diffraction comparison is enough for a theory paper.\n\nThe paper also honestly discusses its approximation in the Discussion, which is good. The novelty claim is a bit stronger than the literature allows in places, but they cite the relevant work on vortex longitudinal spin and situate themselves correctly.\n\nWho is this for: anyone working on focused beams, transverse spin, optical tweezing, or chiral light-matter interactions. It deserves a serious referee. I'd send it to review and expect a useful exchange about the approximations and maybe a request for quantitative validation.","headline":"Focusing a linearly polarized Gaussian beam really does generate local spin and orbital angular momentum from wavefront curvature, and this paper's analytical derivation is the cleanest treatment I've seen.","tokens_in":15623,"tokens_out":2336,"would_cite":true,"duration_ms":23520,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["42.25.Bs","42.25.Ja"],"model":"deepseek-v4-flash","headline":"This paper claims that wavefront curvature, present in strongly focused beams, is a physical source of optical angular momentum: a simply polarized Gaussian beam acquires longitudinal spin and orbital angular momentum near focus, and…","keywords":["optical angular momentum","wavefront curvature","focused Gaussian beams","nonparaxial optics","spin angular momentum","orbital angular momentum","vectorial diffraction","optical helicity"],"falsifier":"Measure the local spin torque on a small absorbing probe placed just off the focal plane of a tightly focused linearly polarized Gaussian beam: Eq. (13) predicts a torque density proportional to $\\sin 2\\phi/(w^2 R[z])|u_{0,0}|^2$, vanishing at $z=0$ and changing sign across the focus. Alternatively, an exact numerical solution that retains all $z$-gradients of the field, rather than truncating as in Eq. (5), should reproduce Eq. (13) in the strong-focusing limit; disagreement would falsify the identification of wavefront curvature as the source.","tokens_in":14633,"feed_emoji":"🌀","tokens_out":6741,"duration_ms":64349,"temperature":0.7,"pith_summary":"Strongly focusing a beam of light makes it nonparaxial, so the field develops a component along the propagation direction and the usual textbook picture of a fixed polarization and zero angular momentum breaks down. The paper shows that for a linearly polarized Gaussian beam, this breakdown produces nonzero longitudinal spin and orbital angular momentum densities near the focal plane, even though the input beam carried neither. For circularly polarized beams, it produces a transverse spin whose sign follows the handedness of the input polarization. The paper traces both effects to the gradient of the wavefront curvature and supports the analytics with vectorial diffraction calculations. If correct, this means an ordinary focused laser beam can transfer angular momentum to small particles without any engineered vortex structure.","feed_headline":"Wavefront curvature creates spin and orbital angular momentum","feed_subtitle":"Focusing a linearly polarized Gaussian beam creates spin and orbital angular momentum near the focal plane.","key_machinery":"The engine is an iterative Maxwell correction to a paraxial Laguerre-Gaussian field, ordered by the paraxial parameter $1/kw$. The first correction is the longitudinal field $E_z^{L1} \\approx (i/k)\\nabla_\\perp\\cdot\\mathbf{E}_{T0}$, which converts the transverse divergence of the input field into a $z$-component; the second-order transverse field $\\mathbf{E}_{T2}$ is then obtained from the Maxwell-Ampere law. In the spin density $\\mathbf{s}_E = (\\epsilon_0/2)\\operatorname{Im}(\\mathbf{E}^*\\times\\mathbf{E})$ and orbital density $\\mathbf{l}_E = (\\epsilon_0/2)\\operatorname{Im}[\\mathbf{r}\\times(\\mathbf{E}^*\\cdot\\nabla)\\mathbf{E}]$, the curvature term enters through $\\operatorname{Im}\\gamma = kr/R[z]$, where $R[z] = (z^2+z_R^2)/z$ is the wavefront curvature. That term is what makes a spinless linearly polarized Gaussian beam acquire longitudinal spin and orbital densities, and a circularly polarized beam acquire helicity-dependent transverse spin, in Eqs. (13), (14), and (20).","core_discovery":"The central claim is that the gradient of the wavefront curvature of a strongly focused Gaussian beam is a genuine source of local optical angular momentum. Working to second order in the paraxial parameter $1/kw$, the authors derive the nonparaxial electromagnetic field of a focused Laguerre-Gaussian beam and compute the electric spin and orbital angular momentum densities. For an input beam with linear polarization, the longitudinal spin density is $s_z = -2\\epsilon_0 r^2/(k w^2 R[z]) \\sin 2\\phi |u_{0,0}|^2$, which is nonzero at all $z$ except the focal plane; the longitudinal orbital density is exactly opposite, so the total longitudinal angular momentum density vanishes. For circular polarization, the transverse spin density acquires a term proportional to $\\sigma k r/R[z]$ that is helicity-dependent, in contrast to the helicity-independent transverse spin of evanescent waves. These are second-order nonparaxial effects, invisible in a paraxial description, and the authors verify them against vectorial diffraction theory.","pith_inferences":["The gradient-of-curvature mechanism suggests that any focused beam with a curved wavefront, not only Gaussian or Laguerre-Gaussian modes, should show analogous angular momentum densities wherever $R[z]$ varies; astigmatic or aberrated wavefronts could produce structured spin patterns.","The helicity-dependent transverse spin identified here points toward a possible route to helicity-sensitive lateral forces on chiral or magnetic nanoparticles, complementing the zero-helicity longitudinal spin.","A position-resolved torque measurement on a small particle near the focus of a linearly polarized Gaussian beam could map the predicted $\\sin 2\\phi$ spatial pattern and give a direct experimental test that does not require a vortex mask."],"forward_implications":["A simple linearly polarized Gaussian beam focused by a high-aperture lens can exert spin torques on absorbing particles near focus, even though the beam carries zero spin angular momentum in the far field.","Focused circularly polarized beams gain a transverse spin whose sign follows the input helicity, offering a way to reverse the direction of transverse spin forces by flipping the handedness of the input light.","The longitudinal spin generated from a linearly polarized Gaussian beam carries zero optical helicity, so it can torque achiral absorbers but cannot drive chiral differential absorption or chiral radiation pressure.","For the linearly polarized Gaussian case the total longitudinal angular momentum density remains zero, meaning spin and orbital parts appear as a local redistribution rather than a net creation of longitudinal angular momentum.","The magnitude and sign of all these angular momentum densities are tunable through the wavefront curvature, i.e. through focusing strength and distance from the focal plane."],"supporting_citations":[{"why":"Supplies the approximation $E_z \\approx (i/k)\\nabla_\\perp\\cdot\\mathbf{E}_{T0}$ from which the first-order longitudinal field and all wavefront-curvature terms follow.","marker":"[24]"},{"why":"Provides the electric spin and orbital angular momentum density definitions used to compute Eqs. (8) and (17).","marker":"[25]"},{"why":"Gives the standard treatment of nonparaxial focused fields with three-dimensional field components, framing the physics the paper builds on.","marker":"[20]"},{"why":"Supplies the Richards-Wolf angular spectrum integration method used for the semi-analytical vectorial diffraction validation.","marker":"[32]"},{"why":"Provides the general vectorial diffraction implementation through which the focused vortex field is propagated and compared with the analytical model.","marker":"[34]"},{"why":"Establishes the transverse and longitudinal angular momenta of light, including the helicity-independent transverse spin of evanescent waves that the paper contrasts with its new helicity-dependent result.","marker":"[9]"},{"why":"Underpins the zero optical helicity result for linearly polarized vortex beams, which the paper extends to its linearly polarized Gaussian spin densities.","marker":"[13]"},{"why":"Supplies the recent result on orbit-induced spin density and zero helicity for tightly focused linearly polarized vortices, used in Section V to connect the new spin density to zero optical helicity.","marker":"[14]"}],"fun_headline_variants":["Curved wavefronts turn simple Gaussian beams into angular momentum sources","Focusing alone can impart spin and orbital momentum to light","Wavefront curvature is a hidden knob for optical angular momentum","Nonparaxial focusing: wavefront curvature generates spin and orbital momentum"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The derivation hinges on Eq. (5), which assumes the field varies along $z$ mainly through the factor $e^{ikz}$, so the longitudinal component is obtained from the transverse divergence of the input field; if the variation of the amplitude or Gouy phase along $z$ is significant in the strong-focusing regime, the wavefront-curvature terms identified here would be modified.","fun_headline_variants_meta":{"raw":{"variants":["Curved wavefronts turn simple Gaussian beams into angular momentum sources","Focusing alone can impart spin and orbital momentum to light","Wavefront curvature is a hidden knob for optical angular momentum","Nonparaxial focusing: wavefront curvature generates spin and orbital momentum"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001005,"raw_usage":{"total_tokens":4239,"prompt_tokens":923,"completion_tokens":3316,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":539,"completion_tokens_details":{"reasoning_tokens":3245}},"tokens_in":539,"tokens_out":3316,"duration_ms":25181,"temperature":1.0,"reasoning_tokens":3245,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T15:35:36.488520+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the local spin torque on a small absorbing probe placed just off the focal plane of a tightly focused linearly polarized Gaussian beam: Eq. (13) predicts a torque density proportional to $\\sin 2\\phi/(w^2 R[z])|u_{0,0}|^2$, vanishing at $z=0$ and changing sign across the focus. Alternatively, an exact numerical solution that retains all $z$-gradients of the field, rather than truncating as in Eq. (5), should reproduce Eq. (13) in the strong-focusing limit; disagreement would falsify the identification of wavefront curvature as the source.","supporting_citations":[{"cited_title":"From maxwell to paraxial wave optics,","cited_arxiv_id":null,"evidence_quote":"Supplies the approximation $E_z \\approx (i/k)\\nabla_\\perp\\cdot\\mathbf{E}_{T0}$ from which the first-order longitudinal field and all wavefront-curvature terms follow."},{"cited_title":"Vector fields in a tight laser focus: com- parison of models,","cited_arxiv_id":null,"evidence_quote":"Provides the electric spin and orbital angular momentum density definitions used to compute Eqs. (8) and (17)."},{"cited_title":"Optical chirality of vortex beams at the nanoscale,","cited_arxiv_id":null,"evidence_quote":"Gives the standard treatment of nonparaxial focused fields with three-dimensional field components, framing the physics the paper builds on."},{"cited_title":"Conservation of the spin and orbital angular momenta in electromagnetism,","cited_arxiv_id":null,"evidence_quote":"Supplies the Richards-Wolf angular spectrum integration method used for the semi-analytical vectorial diffraction validation."},{"cited_title":"Electromagnetic diffraction in optical systems, II. Structure of the image field in an aplanatic system,","cited_arxiv_id":null,"evidence_quote":"Provides the general vectorial diffraction implementation through which the focused vortex field is propagated and compared with the analytical model."},{"cited_title":"parallel","cited_arxiv_id":null,"evidence_quote":"Establishes the transverse and longitudinal angular momenta of light, including the helicity-independent transverse spin of evanescent waves that the paper contrasts with its new helicity-dependent result."},{"cited_title":"Advances in light transverse momenta and opti- cal lateral forces,","cited_arxiv_id":null,"evidence_quote":"Underpins the zero optical helicity result for linearly polarized vortex beams, which the paper extends to its linearly polarized Gaussian spin densities."},{"cited_title":"Orbit-induced localized spin angular momentum in the tight focusing of linearly polarized vortex beams,","cited_arxiv_id":null,"evidence_quote":"Supplies the recent result on orbit-induced spin density and zero helicity for tightly focused linearly polarized vortices, used in Section V to connect the new spin density to zero optical helicity."}],"review_version":1}