{"id":"07540e29-bbd6-4057-8744-094232325400","arxiv_id":"2502.07099","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"First-order perturbative solutions of Maxwell's equations in elastically deformed fibers yield explicit phase-shift and birefringence formulas, with numerical estimates for the GRAVITES experiment.","lead":"This paper calculates how elastic deformations of optical fibers, caused by temperature, pressure, and gravity, change the phase and polarization of light traveling through them. The results provide quantitative estimates of environmental noise for the GRAVITES experiment, a planned kilometer-scale fiber interferometer measuring gravitational redshift on single photons.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The load-bearing concern is the homogeneous bare-cylinder elasticity input: real GRAVITES fibers have a doped core, polymer coating, and spool support, so the Table 3.3 numbers may shift unquantified.","rationale":"The paper's perturbative Maxwell machinery is internally coherent: the multiple-scales scheme, the cokernel condition, and the explicit source terms are laid out in detail, and the conclusion honestly notes that thermo-optic effects are not included. The central Jones-law form (3.51) is a symmetry result that does not depend on material homogeneity. However, the experimental predictions in Table 3.3 are the paper's main deliverable, and they are computed from the Michell displacement of a homogeneous, isotropic, bare cylinder. The reader's weakest-assumption identification is correct and load-bearing: real fibers have a doped core, a coating, and different boundary conditions, so the on-axis strain combinations that enter (3.51) will differ from the idealized values. The paper explicitly flags the homogeneous assumption but does not bound its error, and no code or error bars are provided. This does not invalidate the derivation, but it does mean the quantitative claims are conditional on an unvalidated idealization. A concrete composite-coating FEM check would settle whether the numbers are reliable; until then, the conditional verdict is the appropriate one.","tokens_in":21801,"tokens_out":9748,"duration_ms":94144,"concrete_test":"Implement a 2D plane-strain finite-element or composite-cylinder analytic solution with three layers (GeO2-doped core, silica cladding, acrylate coating) using published material constants, apply the same δT, δp, δg loads with both the rigid-plane support and a spool-wound boundary condition, extract (uxx+uyy)|_{r=0} and (uxx−uyy)|_{r=0}, and recompute Table 3.3 via (3.51). If the entries shift by more than roughly 10%, the quantitative central claim requires revision; if they stay within a few percent, the homogeneous approximation is adequate.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 3 (first paragraph) explicitly assumes 'elastic parameters of the core and cladding to be identical and constant throughout the medium', and Section 3.3 uses the Ref. [1] boundary condition of a bare waveguide resting on a rigid plane. The structural result (3.51) is a symmetry statement and is likely robust, but the quantitative entries in Table 3.3 inherit these idealizations. A real single-mode fiber has a GeO2-doped core with slightly different Young's modulus, Poisson ratio, thermal expansion coefficient, and photoelastic constants than the silica cladding, and it carries a polymer coating whose thermal expansion is roughly two orders of magnitude larger than glass; the GRAVITES fiber is wound on spools, not laid straight on a plane. Since the temperature phase shift (3190 rad for δT=0.01 K) is linear in (uxx+uyy)|_{r=0}, coating-induced thermal stress can add a contribution of the same order, and the gravity-gradient birefringence (−1.27e-6 rad) is proportional to b2, which is set by the support boundary conditions (3.12). These effects are not quantified in the paper, so the claimed experimental estimates rest on an unvalidated idealization.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a first-order perturbative treatment of Maxwell's equations in elastically deformed step-index optical fibers, using a multiple-scales scheme. It derives a Jones-vector propagation law in which the phase shift is controlled by the isotropic on-axis strain (uxx+uyy)|_{r=0} and the birefringence by the quadrupole strain (uxx−uyy)|_{r=0}, with coefficients ξp and ξb that depend only on the fiber's optical parameters. The elastic displacement field is imported from a previous companion paper, and the resulting formulas are evaluated numerically for a straight fiber with GRAVITES-like parameters, yielding concrete estimates in Table 3.3, e.g., 3190 rad photoelastic phase for δT=0.01 K, −53 rad for δp=10 Pa, and −1.27×10⁻⁶ rad birefringence for the gravity gradient.","tokens_in":22028,"tokens_out":10941,"duration_ms":98358,"significance":"If the result holds, it reduces the description of elastic environmental effects on fiber phase and polarization to two strain combinations, which is a clean and useful structural statement for fiber interferometry. The derivation is parameter-free in the sense that no experimental quantity is used to fit the model: all optical and elastic constants come from independent measurements or from the boundary-value problem of the companion paper. The paper also checks its confinement approximation against the full numerical evaluation (Table 3.3) and finds agreement to the displayed precision, which is a genuine internal consistency check. The numerical estimates for GRAVITES are directly usable for experimental design. The main weakness is that the quantitative predictions rest on a simplified elastic model (homogeneous fiber, bare cylinder on a rigid plane) whose error is not quantified.","major_comments":[{"comment":"The calculation assumes that the core and cladding have identical and constant elastic and photoelastic parameters, and that the fiber is a homogeneous cylinder resting on a rigid plane (the boundary condition from Ref. [1], Section 4.2.1). A real single-mode fiber has a GeO2-doped core with different Young's modulus, Poisson ratio, thermal expansion coefficient, and photoelastic constants than the silica cladding, and it carries a polymer coating whose thermal expansion is much larger than that of glass; moreover, the GRAVITES fiber is wound on spools, not laid straight on a plane. Since the Table 3.3 predictions—particularly the 3190 rad temperature phase and the −1.27×10⁻⁶ rad birefringence—are linear in the strain components that are set by this idealization, the claimed numerical results inherit an unquantified systematic uncertainty. The authors should estimate the size of these effects, e.g., by a two-layer elastic model for the core/cladding contrast and a simple model of coating-induced thermal stress, or at least state the expected range over which the numbers can shift.","section":"Section 3, first paragraph; Section 3.3, Table 3.3"},{"comment":"The lengthy source terms Σ(±,∆m) and ˜Γ(±,∆m) are presented as 'are given by' without derivation from Eqs. (3.40) and (3.9). These expressions are the central input to the Jones-matrix result (3.50), and the reader has no way to verify them without repeating a substantial algebra that the paper does not outline. The authors should provide at least a derivation sketch or an appendix showing how the general photoelastic source terms reduce to the displayed forms in terms of (uxx+uyy)|_{r=0} and (uxx−uyy)|_{r=0}.","section":"Section 3.2, Eqs. (3.46)–(3.49)"},{"comment":"There is a notational ambiguity in the ε-scaling of the perturbation. In Eq. (2.29) the perturbation is written as εΣ, and in Eq. (2.33) the bracket contains Σ without a prefactor ε, whereas Eqs. (3.36) and (3.46)–(3.47) define εΣ(±,∆m) as the full source term. It is not clear whether the matrix M in Eq. (2.44) is constructed from Σ or from εΣ; this changes the relation between the ζ-evolution in Eq. (2.45) and the physical phase accumulated over a fiber length L, and therefore affects the numerical values in Table 3.3. Please state the convention explicitly.","section":"Section 2.3 vs. Section 3"}],"minor_comments":[{"comment":"The phrase 'The first column differs due to a correction of the linear thermal expansion coefficient α' is ambiguous: it is not clear which column is meant (the table has a left 'Gravitational phase shift' column and then systematic-effect columns). Please specify that the temperature-phase entry differs from Ref. [1].","section":"Table 3.3 caption"},{"comment":"The displayed expressions contain stray trailing 'y' and 'z' characters at the end of some vector components (e.g., in (3.37) after the fifth component and in (3.49) after the fifth component). These appear to be typesetting artifacts that make the formulas hard to read; please check the source files.","section":"Eqs. (3.37) and (3.49)"},{"comment":"The numerical estimates are quoted without uncertainty propagation. Since the experiment aims to detect small phase shifts, an error budget based on the uncertainties of the material constants in Tables 3.1–3.2 would make the predictions more useful.","section":"Section 3.3"},{"comment":"The captions refer to 'Fiber 1' and 'Fiber 2' without repeating the definitions from the text; please add a sentence in each caption defining the two fibers, since the figures are stated to be the main result of the paper.","section":"Figures 1.2 and 1.3"}],"recommendation":"major_revision","confidential_remarks":"The paper is technically substantial and the structural result (3.51) is likely to be correct and useful. The main risk is not internal inconsistency but the gap between the idealized elastic model and the real experimental configuration; this is fixable by adding quantitative estimates or at least a careful discussion. The overlap with Ref. [1] is appropriate and not circular, since the elastic displacement is derived independently there. The ε-scaling ambiguity in Section 2.3 should be resolved before publication because it currently obstructs verification of the numerical values."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read this if you work on fiber interferometry or environmental noise in guided light. The paper gives the first explicit first-order formulas for elastically induced phase shift and birefringence in single-mode fibers, plus numbers for GRAVITES. The key result is (3.50)-(3.51): the Jones propagation matrix is controlled by two on-axis strain combinations, the trace (uxx+uyy) for phase and the quadrupole (uxx−uyy) for birefringence. That clean symmetry statement will survive refinement of the elastic model. The multiple-scales machinery comes from earlier work, but the combination, the explicit source terms, and the numerical estimates are new.\n\nWhat it does well: the confinement approximation reproduces the full numerical result in Table 3.3, which is good evidence the implementation is internally consistent. The paper also honestly says the thermo-optic effect is not included, and it compares photoelasticity with thermal expansion rather than overselling.\n\nThe soft spot is the one the stress-test flags. Section 3 assumes a homogeneous, bare cylinder with identical core and cladding elastic parameters. Real GRAVITES fibers have a doped core, a polymer coating, and sit on spools. The on-axis strains — and hence the numbers in Table 3.3 — will shift. The coating alone has a thermal expansion roughly two orders of magnitude above silica, so the 3190 rad temperature phase could move by a noticeable fraction. The birefringence entry depends on b2, which is set by the support boundary condition; change the support and you change the number. This is not quantified in the paper. That said, the structural result (3.51) does not depend on those details. The limitation is real but does not sink the derivation.\n\nTwo smaller gaps: the long source terms in Section 3.2 are stated without derivation, so a referee cannot easily verify the algebra; and there is no uncertainty propagation for material constants. Neither is fatal.\n\nWho is this for: the GRAVITES collaboration and anyone modeling environmental phase noise in long fiber interferometers. It deserves a serious referee. I would send it to review and ask the authors to quantify the homogeneous-cylinder assumption, ideally with a coated-fiber or core-cladding contrast check, before the absolute numbers are used in design.","headline":"First explicit formulas for elastic phase and birefringence noise in fiber interferometers; the structural result is solid, but the absolute numbers inherit an unquantified homogeneous-cylinder idealization.","tokens_in":22538,"tokens_out":2692,"would_cite":true,"duration_ms":25344,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["42.81.Gs","42.25.Lc"],"model":"deepseek-v4-flash","headline":"The paper establishes that, to first order, elastic phase shifts and birefringence in single-mode fibers are controlled by just two on-axis strain combinations.","keywords":["single-mode fiber","birefringence","photoelasticity","phase shift","elastic deformation","multiple-scales method","Jones vector","fiber interferometry"],"falsifier":"Measure the phase change of a well-characterized single-mode fiber, with known core and cladding composition, over roughly 100 km while stepping temperature by 0.01 K and pressure by 10 Pa; the observed slopes should match about +3190 rad and −53 rad. A deviation beyond the combined uncertainties of the elastic and photoelastic constants would indicate that the homogeneous-cylinder strain model is insufficient.","tokens_in":21598,"feed_emoji":"🔬","tokens_out":10943,"duration_ms":85499,"temperature":0.7,"pith_summary":"The paper aims to show that all first-order elastic effects of temperature, pressure, and gravity on light in a single-mode fiber reduce to two combinations of strain on the fiber axis: the isotropic sum $(u_{xx}+u_{yy})|_{r=0}$ controls the phase, and the quadrupole difference $(u_{xx}-u_{yy})|_{r=0}$ controls the birefringence. To establish this, it solves Maxwell's equations perturbatively in a deformed, photoelastically anisotropic fiber using a multiple-scales approximation, deriving a Jones-vector propagation law with coefficients that depend only on fiber geometry and refractive index. The concrete payoff is numerical: a $0.01$ K temperature change over $100$ km of fiber yields about $3190$ radians of photoelastic phase, a $10$ Pa pressure change yields about $-53$ radians, and the Earth's gravity gradient produces about $-1.27\\times 10^{-6}$ radians of birefringence. If correct, these formulas let long-baseline fiber interferometry experiments, including those searching for gravitationally induced phase shifts, subtract or cancel environmental elastic noise.","feed_headline":"Two on-axis strains set phase and birefringence in optical fibers","feed_subtitle":"A 0.01 K change over 100 km of fiber should shift phase by about 3190 radians; the paper shows why.","key_machinery":"The load-bearing device is the transport matrix $\\hat M$ of Eq. (2.45), which in all cases considered takes the form $\\hat M = \\xi_p(u_{xx}+u_{yy})|_{r=0}\\,\\sigma_0 + \\xi_b(u_{xx}-u_{yy})|_{r=0}\\,\\sigma_3$. It is assembled from four ingredients: the gauge-fixed Maxwell equations in the Gordon optical metric; the exact step-index fiber modes written with Bessel functions; a multiple-scales expansion along the fiber whose solvability conditions, enforced through the cokernel of the interface matrix, convert first-order perturbations into the Jones-vector law; and the Michell stress-function solution for a homogeneous isotropic elastic cylinder, whose only coupling coefficients are $d_0$ (isotropic strain) and $b_2$ (quadrupole strain). The key identities are $(u_{xx}+u_{yy})|_{r=0} = 2\\mu^{-1}(1-2\\nu)d_0 - 2\\nu\\kappa + 2(1+\\nu)\\alpha(T-T_0)$ and $(u_{xx}-u_{yy})|_{r=0} = -2\\mu^{-1} b_2$, which is why phase and birefringence are set entirely by axis strains.","core_discovery":"At first order in elastic perturbations, the propagation of the Jones vector (the two-component polarization state) along a single-mode fiber obeys $d/d\\zeta\\,(J_x,J_y)^T = i[\\xi_p(u_{xx}+u_{yy})|_{r=0}\\,\\sigma_0 + \\xi_b(u_{xx}-u_{yy})|_{r=0}\\,\\sigma_3](J_x,J_y)^T$. The isotropic on-axis strain $(u_{xx}+u_{yy})|_{r=0}$ multiplies the identity Pauli matrix and therefore advances both linear polarizations equally, which is phase shift. The quadrupole on-axis strain $(u_{xx}-u_{yy})|_{r=0}$ multiplies $\\sigma_3$ and moves the two polarizations oppositely, which is birefringence. The coefficients $\\xi_p$ and $\\xi_b$ depend only on the fiber core radius, the refractive indices, and the optical frequency, and are evaluated numerically from the unperturbed fiber modes. The paper applies this to a fused-silica fiber in a long-baseline interferometer, finding a photoelastic phase of $3190$ rad for $\\delta T = 0.01$ K, $-53$ rad for $\\delta p = 10$ Pa, and a gravity-gradient birefringence of $-1.27\\times 10^{-6}$ rad.","pith_inferences":["Not in the paper: replacing the homogeneous-cylinder assumption with a core-cladding elasticity solution should preserve the two-axis-strain structure but would shift $\\xi_p$ and $\\xi_b$, giving corrected predictions for doped fibers.","Not in the paper: the same perturbative machinery applies to time-dependent translation-invariant strains, such as acoustic waves, so it could model vibration-induced phase noise in fiber sensors and gyroscopes.","Not in the paper: because $\\xi_p$ and $\\xi_b$ depend only on core radius, index contrast, and frequency, the plotted curves could guide fiber design toward operating points with reduced photoelastic sensitivity.","Not in the paper: the close agreement between the confinement approximation and the full solution suggests that closed-form estimates may be accurate enough for engineering budgets in standard single-mode fibers."],"forward_implications":["Environmental phase and birefringence in a single-mode fiber can be budgeted from two on-axis scalar strain values rather than from a full two-dimensional integration over the fiber cross-section.","Photoelastic material response dominates the geometric core-cladding deformation by roughly four orders of magnitude, so interface-shape effects can be neglected in realistic models.","A 0.01 K temperature difference over 100 km of fiber contributes about 3190 radians of photoelastic phase, making temperature stability a primary systematic for long-baseline fiber interferometry.","A 10 Pa pressure difference contributes about $-53$ radians, while the gravity gradient contributes about $-1.1\\times 10^{-6}$ rad of phase and $-1.27\\times 10^{-6}$ rad of birefringence, quantifying the environmental stability needed for gravitational phase-shift experiments.","The confinement approximation, which expands the strain for $r\\ll a$, reproduces the full numerical values, so simplified formulas are available for standard fibers."],"supporting_citations":[{"why":"Supplies the elastic displacement and strain field of a homogeneous isotropic cylinder under gravity, pressure, and temperature, including the coefficients d0 and b2 used here.","marker":"[1]"},{"why":"Provides the original Michell stress-function solution for an elastic cylinder from which the displacement formulas are drawn.","marker":"[13]"},{"why":"Establishes the gauge-fixed Maxwell formulation and the unperturbed step-index fiber modes in terms of Bessel functions that serve as the zeroth-order solution.","marker":"[14]"},{"why":"Extends the perturbative fiber-optics scheme and supplies the Jones-vector transport framework used for polarization dynamics.","marker":"[16]"},{"why":"Provides the multiple-scales method and the cokernel solvability conditions that turn first-order field perturbations into a propagation law for the Jones vector.","marker":"[17]"},{"why":"Defines the experimental long-baseline interferometry parameters used for the numerical estimates in Table 3.3.","marker":"[9]"},{"why":"Gives the photoelastic constants of fused silica used to evaluate the material response numerically.","marker":"[21]"}],"fun_headline_variants":["On-axis fiber strain sets phase and birefringence","Elastic strain on fiber axis drives phase and birefringence","Phase shift and birefringence from on-axis strain","Strain on fiber axis controls light's phase and polarization","Tiny temperature change yields 3190 rad phase shift in fiber"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The computation assumes a homogeneous, isotropic elastic cylinder whose core and cladding share identical elastic and photoelastic constants, so the predicted on-axis strain in a real doped fiber could differ.","fun_headline_variants_meta":{"raw":{"variants":["On-axis fiber strain sets phase and birefringence","Elastic strain on fiber axis drives phase and birefringence","Phase shift and birefringence from on-axis strain","Strain on fiber axis controls light's phase and polarization","Tiny temperature change yields 3190 rad phase shift in fiber"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000267,"raw_usage":{"total_tokens":1614,"prompt_tokens":948,"completion_tokens":666,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":564,"completion_tokens_details":{"reasoning_tokens":578}},"tokens_in":564,"tokens_out":666,"duration_ms":6484,"temperature":1.0,"reasoning_tokens":578,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T13:47:40.191031+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the phase change of a well-characterized single-mode fiber, with known core and cladding composition, over roughly 100 km while stepping temperature by 0.01 K and pressure by 10 Pa; the observed slopes should match about +3190 rad and −53 rad. A deviation beyond the combined uncertainties of the elastic and photoelastic constants would indicate that the homogeneous-cylinder strain model is insufficient.","supporting_citations":[{"cited_title":"Barzegar, P","cited_arxiv_id":null,"evidence_quote":"Supplies the elastic displacement and strain field of a homogeneous isotropic cylinder under gravity, pressure, and temperature, including the coefficients d0 and b2 used here."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the original Michell stress-function solution for an elastic cylinder from which the displacement formulas are drawn."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the gauge-fixed Maxwell formulation and the unperturbed step-index fiber modes in terms of Bessel functions that serve as the zeroth-order solution."},{"cited_title":"Fiber optics in curved space-times","cited_arxiv_id":"2410.23048","evidence_quote":"Extends the perturbative fiber-optics scheme and supplies the Jones-vector transport framework used for polarization dynamics."},{"cited_title":"Hilweg, F","cited_arxiv_id":null,"evidence_quote":"Defines the experimental long-baseline interferometry parameters used for the numerical estimates in Table 3.3."},{"cited_title":"Primak and D","cited_arxiv_id":null,"evidence_quote":"Gives the photoelastic constants of fused silica used to evaluate the material response numerically."}],"review_version":1}