{"id":"78a0c1f2-a6c6-4f1f-940d-830e7cfe80a9","arxiv_id":"2505.00170","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"A Monte Carlo model predicts that a three-grating neutron moiré interferometer with a one-ton lead mass can measure G to 150 ppm, and reanalysis suggests lunar tides may add tens of ppm to some prior torsion pendulum results.","lead":"This paper designs a neutron interferometer that could measure the gravitational constant G to about 150 parts per million, and it uses Monte Carlo simulations to estimate the signal and error sources. It matters because G is the least precisely known fundamental constant, and a measurement method with systematics independent of torsion balances could help resolve long-standing discrepancies between experiments.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. C6 inverts the contrast dependence of the phase uncertainty, understating the shot-noise limit; the 130 ppm statistical budget and the 150 ppm central claim are therefore unsupported.","rationale":"The central claim is the projected 150 ppm total uncertainty in G: 73 ppm systematic plus 130 ppm statistical. The statistical term is not a minor extrapolation; it is computed from Eq. C6, and that equation is internally inconsistent with Eq. C7 and with the preceding per-pixel derivation. A correct fringe-fitting variance scales as 2/(N C^2) for small C, so the required number of counts to reach 92.7 ppm of a 13.2 urad phase is about 2.7 x 10^19. The proposed beamline, even at 10^9 cm^-2 s^-1 through the stated slit, delivers about 3 x 10^15 counts in 120 d, roughly four orders of magnitude too few. This is not a question of achievable contrast; it is a shot-noise floor. The reader's concern about unvalidated 22.1% contrast is reasonable and remains, but the more decisive issue is that the paper's own statistical formula would need to be corrected before the projected uncertainty can be taken seriously. The model code and validation against two-grating PGMI data are welcome, but they do not address this analytic scaling error. For these reasons the current manuscript should not be accepted as a demonstration of 150 ppm feasibility; a revised version could resubmit with a corrected statistical analysis and, if the conclusion survives, a much longer beam-time estimate.","tokens_in":16951,"tokens_out":20405,"duration_ms":208588,"concrete_test":"Re-derive Eq. C6 from Eq. C3-C5 by summing the per-pixel inverse variances over M pixels: the correct result is (sigma_DeltaPhi)^2 = 1/[M N_bar (1 - sqrt(1-C^2))], not Eq. C6. Then evaluate it with N_total = 120 d x (10^9 cm^-2 s^-1 x 0.28 cm^2) and C = 0.221, computing sigma_Phi / 13.2 urad. If this ratio is about 9000 ppm, the statistical budget in Table II and the 150 ppm claim fail. As a cross-check, verify that Eq. C6 and Eq. C7 are mutually inconsistent for small C by a factor C^4/4.","verdict_should_be":"REJECT","load_bearing_attack":"Appendix C contains an internal inconsistency that is more load-bearing than the assumed contrast. Eq. C6 gives (sigma_Phi)^2 = [1 - sqrt(1-C^2)]/(N_dot T), which for C=0.221 is about 0.0247/(N_dot T). Summing the per-pixel variances from Eq. C3 over M pixels instead yields (sigma_Phi)^-2 = M N_bar [1 - sqrt(1-C^2)], so the factor [1 - sqrt(1-C^2)] belongs in the denominator, and the number of fitted pixels/total counts must appear. For small C this is the standard sigma_Phi approx sqrt(2/(N C^2)); it grows as contrast falls. Eq. C6, by contrast, shrinks as C -> 0, which is unphysical and also contradicts the paper's own Eq. C7, T approx 2/(N_dot (sigma_Phi C)^2), the standard result used for the C * sigma optimization. With the stated 10^9 cm^-2 s^-1 fluence through a 35 mm x 0.8 mm slit, 120 d provides at most about 3 x 10^15 detected counts. The corrected formula gives sigma_Phi approx 1.2 x 10^-7 rad, about 9000 ppm of the 13.2 urad gravitational phase, not 92.7 ppm. Even if the 22.1% contrast were achieved and all systematic terms were zero, the proposed experiment would not reach 150 ppm. The central feasibility claim therefore rests on a statistical error, not merely on an unvalidated contrast.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes a three-grating phase-grating moiré neutron interferometer (3-PGMI) with a 1 t lead source mass as a new way to measure the gravitational constant G. It introduces a Monte Carlo path-integral model that computes phases along all interfering neutron trajectories in arbitrary potentials, validates the model against the van der Zouw single-grating interferometer and against two-grating PGMI data, optimizes a specific interferometer geometry, and compiles an uncertainty budget. The central claim is that a 150 ppm measurement of G is achievable in the near term, with 73 ppm systematic and 130 ppm statistical contributions, along with a separate estimate that lunar tidal forces add 39.5 ppm to the uncertainty of the Luther–Towler torsion-pendulum result.","tokens_in":17260,"tokens_out":9051,"duration_ms":87772,"significance":"If the statistical and systematic estimates were correct, this would be a genuinely new technique for measuring G with systematic effects independent of pendulum and atom-interferometer experiments, and the Monte Carlo tool for nonlinear potentials would be a useful addition to the neutron-interferometry literature. The manuscript has real strengths: the model is validated against published experimental data, the code is made available in a public repository, the systematic catalog is fairly complete, and the lunar-tide calculation for the Luther–Towler experiment is a concrete and checkable result. However, the feasibility claim rests on a phase-fit statistics formula that is inconsistent with the paper's own derivation; the corrected formula changes the statistical uncertainty by roughly two orders of magnitude, so the central 150 ppm claim is not supported by the stated beam parameters.","major_comments":[{"comment":"Equation (C6) inverts the contrast dependence of the phase uncertainty. Equations (C3)–(C5) derive the sum of inverse per-pixel variances as \\bar{N}(1 - sqrt(1 - C^2)); by Eq. (C4) the phase variance is the reciprocal of this quantity, so the correct result is (\\sigma_\\Phi)^2 = 1/[\\dot{N} T (1 - sqrt(1 - C^2))], not (1/\\dot{N} T)(1 - sqrt(1 - C^2)) as printed. The printed version makes \\sigma_\\Phi shrink as C → 0, which is unphysical and contradicts the paper's own Eq. (C7). With N ≈ 2.9 × 10^15 counts in 120 d (from the 10^9 cm^-2 s^-1 fluence and 35 mm × 0.8 mm slit) and C = 0.221, the corrected formula gives \\sigma_\\Phi ≈ 1.2 × 10^-7 rad, about 9000 ppm of the 13.2 µrad gravitational phase, rather than 92.7 ppm; the combined null/source statistical uncertainty is then about 12,700 ppm, not 130 ppm. Even with perfect contrast, a 120 d run would give roughly 2000 ppm per phase measurement. The 150 ppm central claim is therefore not supported by the stated beam parameters.","section":"Appendix C, Eq. (C6)"},{"comment":"The source-mass density input in Table I is inconsistent with the budget entry in Table II. Table I lists a density uncertainty of 1.16 µg cm^-3 for lead (ρ = 11,342 kg/m^3), which is about 0.1 ppm and is the value quoted for PbWO4 optical characterization in Sec. IV. A 2.0 ppm contribution, as listed in Table II, would require a density uncertainty of about 23 µg cm^-3 for lead. The paper should either adopt a lead-specific density uncertainty consistent with 2 ppm or reduce the Table II entry to the 0.1 ppm level implied by Table I; as written, the density uncertainty is not traceable to the stated parameter.","section":"Table I, Sec. V.B, and Table II"},{"comment":"The assumed 22.1% contrast for the optimized 3-PGMI is load-bearing and is not validated by the experiments cited. Appendix A explicitly notes that three-grating PGMI experiments have so far achieved only low contrast [28], and the model validation in Figs. 5 and 6 uses a single-grating NI and two-grating PGMI data, not a 3-PGMI under gravitational loading. Because the statistical precision scales as 1/C in the corrected formula (and in the paper's own Eq. C7), an unvalidated factor of even 2–3 in contrast changes the feasibility conclusion. The authors should present a validated measurement of 3-PGMI contrast in the proposed geometry, or treat the contrast as an unknown and propagate its uncertainty into the final uncertainty budget.","section":"Appendix A and Sec. IV (contrast)"}],"minor_comments":[{"comment":"The phrase \"new meteorological devices\" should read \"new metrological devices.\"","section":"Section VI"},{"comment":"The dotted lines in Fig. 4 are labeled only as relative ppm values; the caption should state explicitly that the ppm is relative to the 13.2 µrad gravitational phase so that the reader can connect the figure to the uncertainty budget.","section":"Sec. IV, Fig. 4"},{"comment":"The Sagnac expression uses L/2 as the effective lever arm for the 3-PGMI area, but the derivation of this geometric factor for the three-grating arrangement is not given; a short derivation or reference would improve reproducibility.","section":"Eq. (15) and Sec. V.A.2"},{"comment":"Reference [51] is a NIST disclaimer rather than a specification for displacement sensors; please cite the specific commercial specifications or replace it with a metrology reference that supports the stated 3 µm and 0.5° uncertainties.","section":"Reference [51]"},{"comment":"The table lists \"Mass Surface Variation δLx,δLy,δLz\" with values of 5 µm but does not state the units in the table body; the caption should make the units explicit and consistent with the dimension uncertainties.","section":"Table I"}],"recommendation":"reject","confidential_remarks":"The Appendix C statistical error is not a presentation issue; it invalidates the paper's headline result. Even with the assumed contrast, the stated beam parameters give a statistical uncertainty roughly two orders of magnitude larger than claimed, so the proposed experiment cannot reach 150 ppm as designed. The Monte Carlo model and the lunar-tide calculation have value and could form the basis of a future manuscript if the feasibility analysis is redone with corrected statistics and a validated 3-PGMI contrast."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Kapahi et al. propose a three-grating phase-grating moiré neutron interferometer to measure G with a 1 t lead source mass, and they build a Monte Carlo path-integral model for nonlinear potentials. That model is the genuinely new piece: it is validated against the van der Zouw phase-grating NI and two-grating PGMI contrast data, and the code is public. The systematic uncertainty budget is thoughtful and mostly credible, with the atmospheric, Sagnac, alignment, and grating-period terms each worked out. Credit where due: this is a serious design study, not a crank proposal.\n\nThe problem is the statistical uncertainty, and it is load-bearing. Appendix C derives Eq. C6, (σΦ)^2 = [1−√(1−C^2)]/(ṆT). That has the contrast dependence inverted. Integrating the per-pixel variances in Eq. C3 gives (σΦ)^−2 = ṄT [1−√(1−C^2)], so the correct expression is (σΦ)^2 = 1/[ṆT(1−√(1−C^2))], which for small C is ≈ 2/(ṆT C^2). Their own Eq. C7 is the standard version and is consistent with the corrected result, not with Eq. C6. With the stated 10^9 cm^−2 s^−1 through the 35 mm × 0.8 mm slit, a 120 d run gives roughly 3×10^15 counts at best. The corrected formula yields σΦ ≈ 1.2×10^−7 rad, about 9000 ppm of the 13.2 µrad gravitational phase, not the 92.7 ppm they quote. Even summing null and source-mass runs, the statistical floor is ~1.3×10^4 ppm, two orders of magnitude above 150 ppm. The central feasibility claim therefore does not survive.\n\nOther issues are minor by comparison. The 22.1% contrast is an unvalidated assumption for the 3-PGMI, as the paper admits three-grating PGMI experiments have only shown low contrast. Table I lists a lead density uncertainty of 1.16 µg cm^−3, which they attribute to PbWO4 optical characterization, while Table II quotes 2.0 ppm from density; those numbers don't match for lead. The lunar-tide correction to Luther-Towler is stated with a sketch, not a full derivation. All of these are fixable.\n\nThe paper is worth a serious referee because the model and systematic treatment are reusable even if the headline sensitivity is wrong. But the conclusion needs a major revision: as written, the proposed experiment would require thousands of times more beam time or a much larger phase shift to reach 150 ppm. I'd recommend peer review with a request for a corrected statistical analysis, then re-evaluation.","headline":"The new PGMI model and systematic budget are real, but the central 150 ppm claim rests on an inverted statistics formula; the projected sensitivity is off by roughly two orders of magnitude.","tokens_in":17857,"tokens_out":6770,"would_cite":false,"duration_ms":58875,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A three-grating neutron moiré interferometer with a one-tonne lead source mass could measure the gravitational constant to 150 parts per million, using systematic errors independent of pendulum experiments.","keywords":["gravitational constant","neutron interferometry","phase-grating moiré interferometer","Monte Carlo simulation","WKB phase integral","lunar tidal forces","systematic uncertainty budget","cold neutrons"],"falsifier":"Set up the optimized 3-PGMI geometry (700 nm gratings, separations 80.0, 200.0, 200.9, and 399.1 cm, 4 Å polychromatic beam, 1 t lead source mass) and measure the fringe contrast and the source-mass phase shift. If the contrast falls far below 22.1% or the phase departs from the predicted $13.2\\,\\mu$rad by more than the model's error budget, the 150 ppm projection fails; matching both would directly support the design.","tokens_in":16706,"feed_emoji":"⚛️","tokens_out":14880,"duration_ms":137132,"temperature":0.7,"pith_summary":"The paper proposes a practical neutron interferometer that could measure the gravitational constant $G$ to about 150 parts per million, matching the precision of the 2022 recommended value but with a different set of systematic errors. The key idea is to use a three-grating phase-grating moiré interferometer (PGMI), which accepts a broad spectrum of neutron wavelengths and therefore far more flux than single-crystal neutron interferometers used in past gravity experiments. A new Monte Carlo model computes each neutron's phase along its actual trajectory through the gravitational field of a 1 t lead source mass, and the model is validated against earlier grating-interferometry measurements. The model yields an uncertainty budget dominated by grating period, neutron wavelength, and source-mass positioning, and it also shows that lunar tidal forces contribute roughly 40 ppm of extra uncertainty to a standard torsion-pendulum measurement of $G$.","feed_headline":"Neutron moiré interferometer aims to measure G to 150 ppm","feed_subtitle":"Three gratings and a 1-tonne lead mass would match the 2022 G precision with independent systematics.","key_machinery":"The central object is the three-grating phase-grating moiré interferometer (3-PGMI): three nanofabricated phase gratings with heights $\\pi/2$, $\\pi$, and $\\pi/2$ whose overlapping diffraction orders form a moiré fringe at the camera that is independent of incident neutron angle and wavelength. The carrying identity is the first-order WKB phase integral $\\gamma_{lmn} = \\int_{S_{lmn}} \\mathbf{k}(\\mathbf{r})\\cdot d\\mathbf{s}$, expanded as $k \\approx k_0\\left(1 + m_N U_{\\mathrm{tot}}/\\hbar^2 k_0^2\\right)$, which lets an arbitrary nonlinear potential—including the $1/r$ source-mass potential—be accumulated along each path. A Monte Carlo simulation varies neutron position, angle, and wavelength, sums interfering path pairs incoherently, and produces both the fringe pattern and the averaged expression for $G$ in which phase differences, corrections, and a path-integrated source-mass response factor are wavelength- and path-averaged.","core_discovery":"On its own terms, the claim is that a neutron 3-PGMI with a 1 t rectangular lead source mass can determine $G$ with a total relative uncertainty of 150 ppm (73 ppm systematic, 130 ppm statistical), using 240 days of beam time split equally between source-mass and null runs. The measurable is the phase shift of the interference fringes induced by the source-mass potential; the model projects this phase to $13.2\\,\\mu$rad at 22.1% fringe contrast. To make the extraction, the paper introduces a Monte Carlo path-integral model in which the WKB phase $\\gamma = \\int \\mathbf{k}(\\mathbf{r})\\cdot d\\mathbf{s}$ is accumulated along every neutron trajectory, with the neutron wavevector depending on the total potential including source mass, Earth rotation, air, and Moon. It uses this model to propagate experimental parameter uncertainties into $G$, and applies the lunar part of the model to a torsion-pendulum G measurement, finding that lunar gravitation adds about 39.5 ppm to that measurement's uncertainty unless corrected.","pith_inferences":["Not stated in the paper: the same Monte Carlo model could be inverted to calibrate the instrument in situ, using the fringe phase from a known calibration mass to recover the effective wavelength distribution and shrink the 40 ppm wavelength term.","Not stated in the paper: the lunar-tidal calculation generalizes to any G measurement whose test mass has finite extent, so recomputing historical torsion-pendulum results with the exact lunar ephemeris and laboratory coordinates could expose a correction-dependent shift between reported values.","Not stated in the paper: because the required beam time scales inversely with the square of the product of phase uncertainty and contrast, doubling the fringe contrast would cut a 240 d campaign to about one quarter of its length, making grating alignment a higher-leverage upgrade than a larger source mass."],"forward_implications":["A 240 d measurement, split between a source-mass run and a null run, would yield G with a total relative uncertainty of 150 ppm (73 ppm systematic, 130 ppm statistical), comparable to the 2022 recommended value.","The same geometry would need more than 800 d to reach 100 ppm, so near-term gains would have to come from higher flux, smaller grating periods, or a stronger source-mass potential.","Lunar tidal forces contribute about 39.5 ppm of additional uncertainty to a standard torsion-pendulum measurement, implying published G values may need date- and location-specific tidal corrections.","Monitoring grating alignment with x-rays between runs, independent of the gravitational signal, addresses the main experimental risk that the three-grating fringes have low contrast.","The broadband PGMI model also points to neutron measurements of magnetic structures and of the neutron electric dipole moment as further applications."],"supporting_citations":[{"why":"Defines the current recommended value of G and its roughly 22 ppm relative uncertainty, the precision target of the proposed measurement.","marker":"[1]"},{"why":"Supplies the WKB phase-integral treatment and the atmospheric phase-correction formula used to build the model.","marker":"[22]"},{"why":"Provides very-cold-neutron phase-grating interferometry data under gravity, used to validate the model against a gravitational potential.","marker":"[26]"},{"why":"Establishes the universal moiré effect that defines the interfering path pairs in the three-grating PGMI.","marker":"[27]"},{"why":"Demonstrates the three phase-grating moiré neutron interferometer that the proposal extends, and reports the low contrast achieved so far.","marker":"[28]"},{"why":"Supplies white-beam far-field PGMI experimental data used to validate the Monte Carlo model's fringe contrast.","marker":"[29]"},{"why":"Provides a momentum-space PGMI model and two-grating experimental parameters against which the new model is checked.","marker":"[31]"},{"why":"Gives the analytic gravitational potential of a right rectangular prism used for the source-mass phase integral.","marker":"[35]"},{"why":"The torsion-pendulum G measurement whose lunar tidal uncertainty is calculated as about 39.5 ppm.","marker":"[53]"},{"why":"Supplies the chopper-spectrometer wavelength-calibration method whose uncertainty sets the 40 ppm neutron-wavelength contribution.","marker":"[55]"}],"fun_headline_variants":["Neutron moiré interferometer targets G to 150 ppm","Moiré neutron interferometer to measure G at 150 ppm","New neutron interferometer design aims for G to 150 ppm","Moiré neutron interferometer could measure G to 150 ppm","Neutron moiré interferometer: 150 ppm G measurement"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The budget's load-bearing premise is that the three-grating interferometer with the 1 t lead mass in place will actually show the 22.1% fringe contrast the model predicts; three-grating neutron PGMI experiments to date have shown only low contrast, while the validation data come from two-grating and single-grating setups.","fun_headline_variants_meta":{"raw":{"variants":["Neutron moiré interferometer targets G to 150 ppm","Moiré neutron interferometer to measure G at 150 ppm","New neutron interferometer design aims for G to 150 ppm","Moiré neutron interferometer could measure G to 150 ppm","Neutron moiré interferometer: 150 ppm G measurement"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000799,"raw_usage":{"total_tokens":3566,"prompt_tokens":1049,"completion_tokens":2517,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":665,"completion_tokens_details":{"reasoning_tokens":2425}},"tokens_in":665,"tokens_out":2517,"duration_ms":18165,"temperature":1.0,"reasoning_tokens":2425,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T04:48:48.764051+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Set up the optimized 3-PGMI geometry (700 nm gratings, separations 80.0, 200.0, 200.9, and 399.1 cm, 4 Å polychromatic beam, 1 t lead source mass) and measure the fringe contrast and the source-mass phase shift. If the contrast falls far below 22.1% or the phase departs from the predicted $13.2\\,\\mu$rad by more than the model's error budget, the 150 ppm projection fails; matching both would directly support the design.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the current recommended value of G and its roughly 22 ppm relative uncertainty, the precision target of the proposed measurement."},{"cited_title":"Dynamic measure- ment of gravitational coupling between resonating beams in the hertz regime,","cited_arxiv_id":null,"evidence_quote":"Supplies the WKB phase-integral treatment and the atmospheric phase-correction formula used to build the model."},{"cited_title":"Test of a sin- gle crystal neutron interferometer,","cited_arxiv_id":null,"evidence_quote":"Provides very-cold-neutron phase-grating interferometry data under gravity, used to validate the model against a gravitational potential."},{"cited_title":"High-precision gravity measurements using atom interferometry,","cited_arxiv_id":null,"evidence_quote":"Establishes the universal moiré effect that defines the interfering path pairs in the three-grating PGMI."},{"cited_title":"Neutron Interferometry: Lessons in Experimental Quantum Mechanics,","cited_arxiv_id":null,"evidence_quote":"Supplies white-beam far-field PGMI experimental data used to validate the Monte Carlo model's fringe contrast."},{"cited_title":"Three Phase-Grating Moir´ e Neutron Interferometer for Large Interferometer Area Applications,","cited_arxiv_id":null,"evidence_quote":"Gives the analytic gravitational potential of a right rectangular prism used for the source-mass phase integral."},{"cited_title":"Redetermination of the newtonian gravitational constant G,","cited_arxiv_id":null,"evidence_quote":"The torsion-pendulum G measurement whose lunar tidal uncertainty is calculated as about 39.5 ppm."},{"cited_title":"The IAU 2009 system of astronomical constants: the report of the iau working group on numerical standards for fundamental astron- omy,","cited_arxiv_id":null,"evidence_quote":"Supplies the chopper-spectrometer wavelength-calibration method whose uncertainty sets the 40 ppm neutron-wavelength contribution."}],"review_version":1}