{"id":"c0c53902-37aa-411d-b3ec-52c42900f7f6","arxiv_id":"2411.09912","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Using muonic 199Hg hyperfine data and atomic many-body calculations, the authors extract absolute Bohr-Weisskopf effects in Hg isotopes and find an empirical Moskowitz-Lombardi additive constant that is smaller than the Fujita-Arima value.","lead":"Mercury's nucleus is not a point magnet: how its magnetism is spread out changes the hyperfine structure of its atoms. The authors used muonic mercury data and atomic calculations to measure this Bohr-Weisskopf effect and found that the additive constant in a 50-year-old empirical rule is smaller than previously claimed.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Reported uncertainty on the ML additive constant omits the common-mode 199Hg anchor error; with it, c_exp is only ~1.5 sigma from the Fujita-Arima value, so the claimed significant deviation is not established.","rationale":"The reader's formal weakest assumption is Eq. (4), the neglect of the differential Breit-Rosenthal term in the neutral/ion screening ratio. That is a legitimate concern, and the numerical check should be shown. However, the most load-bearing issue for the paper's central claim is statistical: the fitted additive constant c is determined from an absolute anchor shared by every point in the fit. The quoted uncertainty 0.0008 is smaller than what an ordinary unweighted least-squares fit with the stated 0.44% errors would give, and far smaller than the common-mode anchor uncertainty, which cannot average down. This means the abstract's statement that the empirical value 'differs significantly' from Fujita-Arima is unsupported as written. The authors presumably used some weighting or omitted the common-mode term; either way, the fit covariance must be disclosed. The direct atomic extraction agrees with the muonic value but has twice the uncertainty and does not resolve the difference. The data may still support an empirical Moskowitz-Lombardi rule with c consistent with zero, but the claimed significant deviation from c_FA = -0.01 is not established. This is an addressable statistical issue, so I would keep the paper conditional rather than reject it outright, but the headline claim needs to be softened or the covariance shown to justify the reported precision.","tokens_in":11915,"tokens_out":12461,"duration_ms":134094,"concrete_test":"Refit Eq. (6) to the ten neutral-Hg epsilons of Table III using the full covariance matrix, with the 0.44% 199Hg-anchor uncertainty entered as a fully correlated common-mode component and the differential-anomaly errors from Ref. [60] entered as independent components. If the resulting sigma_c is about 0.0044 rather than the reported 0.0008, the significance claim must be withdrawn or reduced to 'consistent within about 1.5 sigma'; for robustness, also repeat the fit using the directly extracted ionic anchor -3.4(10)% instead of the muonic anchor.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing flaw is not Eq. (4) but the treatment of the fit that yields c_exp. The quoted c_exp = -3.6(8) x 10^-3 is used to claim a significant difference from c_FA = -0.01. All ten neutral-Hg epsilons in Table III share the same 199Hg anchor, epsilon_199 = -2.34(44)%. In the fit y_i = b + alpha/|mu_i| with y_i = -epsilon_i, a common shift of this anchor changes every y_i and therefore shifts b, and hence c = -b, by the full 0.0044; it cannot be reduced by the number of isotopes. Even a naive unweighted fit with the stated 0.44% errors gives sigma_c ~ 0.0021, and a correct common-mode treatment gives ~0.0044. With sigma_c = 0.0044, c_exp is only about 1.5 sigma from -0.01, so the headline result that the additive constant differs significantly from the Fujita-Arima value is not supported. The Eq. (4) screening-factor issue is secondary: a small residual there shifts c far less than this covariance effect.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper extracts the Bohr-Weisskopf (BW) effect for neutral and singly-ionized mercury isotopes from the measured hyperfine splitting of muonic 199Hg, using a two-step translation through H-like mercury with five nuclear wavefunction shapes, then applies electronic screening factors to obtain values for Hg+ and neutral Hg. The authors also perform an all-orders atomic many-body calculation for 199Hg+ and extract the BW effect directly from the experimental hyperfine constant. From the isotope-dependent BW values, they fit the Moskowitz-Lombardi rule epsilon = c - alpha/|mu| and report c_exp = -3.6(8) x 10^-3 and alpha_exp = 0.98(5) x 10^-2 muN for neutral Hg, concluding that the additive constant differs significantly from the Fujita-Arima value c_FA = -0.01.","tokens_in":12107,"tokens_out":7342,"duration_ms":73103,"significance":"If the central claim were fully supported, the paper would provide the first empirical determination of the absolute BW effect in mercury isotopes and a direct test of the additive constant in the Moskowitz-Lombardi rule, with implications for precision atomic clock measurements and nuclear structure studies. The extraction chain is explicit, the nuclear-model dependence is honestly assessed with five single-particle wavefunction shapes, and the independent all-orders atomic calculation serves as a valuable cross-check. However, the significance of the headline result is undermined by the common-mode uncertainty in the 199Hg anchor, as detailed below.","major_comments":[{"comment":"The quoted uncertainty on c_exp omits the fully correlated contribution from the common 199Hg anchor. All ten epsilons in Table III are derived by adding a differential anomaly to the 199Hg value epsilon_199 = -2.34(44)%, so a common shift of epsilon_199 shifts every point in the fit of Eq. (6) and therefore changes c by exactly the same amount. This 0.0044 absolute uncertainty is not reduced by the number of isotopes; the correct uncertainty on c is at least this large, and the quoted c_exp = -3.6(8) x 10^-3 uses an uncertainty of only 0.0008, which treats the ten values as independent. With the common-mode error included, the difference between c_exp and c_FA = -0.01 is only about 1.5 standard deviations, so the claim that the additive constant 'differs significantly' from the Fujita-Arima value is not supported. Please redo the fit with a covariance matrix that includes the fully correlated anchor error, and revise the abstract, title, and conclusions accordingly if the significance disappears.","section":"Table III, Eq. (10), surrounding text"},{"comment":"The assertion 'We have checked numerically that 199Delta201 is very small and may be omitted in the ratio (4)' is not backed by any numerical result or citation. Because Eq. (4) is the basis for translating the ion BW values to the neutral atom and therefore for all neutral-Hg entries in Table III, please provide the actual magnitude of 199Delta201 (or its ratio to the differential BW effect) so the reader can assess the size of the omitted Breit-Rosenthal term relative to the 0.4% uncertainty on the right-hand side of Eq. (4). Without this number, the systematic error on the neutral-atom BW values is not fully documented.","section":"Eq. (4) and following sentence"}],"minor_comments":[{"comment":"The two values listed for 199Hg+ (-2.59(49)% and -3.4(10)%) are visually confusing; consider using two explicitly labeled rows (e.g., 'muonic-deduced' and 'direct extraction') or a footnote that clearly separates the two determinations.","section":"Table II"},{"comment":"The convention for the isotope ordering in the differential anomaly 1Delta2 should be stated explicitly (i.e., whether it corresponds to epsilon(1) - epsilon(2) or the reverse), to avoid sign ambiguity when the definition is used with data from Ref. [60].","section":"Eq. (5)"},{"comment":"The statement 'The deviations arise from omitted many-body nuclear effects in the theory analysis of Ref. [32]' is a causal claim that goes beyond the data presented; consider softening it to 'may arise' unless a specific calculation supporting this attribution is included.","section":"Paragraph after Eq. (10)"}],"recommendation":"major_revision","confidential_remarks":"The paper contains a useful and honest extraction of absolute BW effects and a valuable atomic many-body cross-check, but the headline conclusion about the Moskowitz-Lombardi additive constant depends on a fit that ignores the common-mode error in the 199Hg anchor. The authors should redo the fit with a proper covariance treatment; if the statistical significance of the deviation from c_FA disappears, the title, abstract, and conclusions need to be revised accordingly. The data and method remain publishable even with a weaker claim."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThe thing to know about this paper: it gives the first empirical absolute Bohr-Weisskopf effects for electronic mercury, for ten isotopes and isomers, in both neutral and singly-ionized form, and fits the Moskowitz-Lombardi rule to them. The extraction chain from muonic 199Hg is explicit and careful, and the authors include an independent all-orders atomic calculation that agrees within its larger uncertainty. If the absolute values hold up, they are useful for hyperfine calculations in Hg clocks and EDM searches.\n\nThe soft spot is the fit that produces the headline claim. The BW effects in Table III all share the same anchor: the 199Hg muonic value, with an uncertainty of 0.44%. A shift in that anchor shifts every epsilon by the same amount, so the intercept of the ML fit, and therefore the additive constant c, carries at least that full 0.0044 uncertainty. The quoted c_exp = -3.6(8) x 10^-3 appears to treat the points as independent, giving an uncertainty that is far too small. With sigma_c = 0.0044, c_exp is only about 1.5 sigma from the Fujita-Arima value of -0.01. The claim that the additive constant 'differs significantly' is not supported. This is fixable: report the fit with a covariance matrix or a common-mode uncertainty, and temper the conclusion.\n\nA secondary weakness is Eq. (4), where the neutral/ion screening factor is taken from the ratio of measured differential anomalies and the differential Breit-Rosenthal contribution is dropped. The authors say they checked numerically that it is very small, but they do not show it. A referee should ask for that check.\n\nOverall, the paper is a serious contribution to a well-motivated problem. The nuclear model dependence is honestly treated with five wavefunction shapes, the uncertainty budget is mostly transparent, and the direct atomic extraction is a reasonable cross-check. I would send it to peer review, and I think it should be published after the statistics of the fit are corrected. The empirical BW values themselves are the real product, and they will be cited.","headline":"Valuable new absolute BW values for mercury, but the claimed significant ML-vs-FA deviation does not survive the common-mode muonic anchor error.","tokens_in":12696,"tokens_out":4104,"would_cite":true,"duration_ms":39724,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper determines the Bohr-Weisskopf hyperfine correction in mercury isotopes and shows that the additive constant in the Moskowitz-Lombardi rule is much smaller than the standard theory value.","keywords":["Bohr-Weisskopf effect","hyperfine anomaly","Moskowitz-Lombardi rule","mercury isotopes","muonic atoms","hyperfine structure","atomic many-body calculations","nuclear magnetic moments"],"falsifier":"Calculate the differential Breit-Rosenthal contribution $^{199}\\delta^{201}$ for the $6s6p\\,{}^3P_1$ state of neutral Hg with a relativistic many-body method; if it is comparable to the 0.4% uncertainty in the ratio 0.902(4), the neutral-atom BW values and the fitted constant $c_{\\mathrm{exp}}$ shift by more than their quoted errors. Alternatively, a new measurement of the muonic $^{199}$Hg 1s hyperfine splitting with uncertainty well below 0.12 keV would directly test the seed value $-68(8)\\%$.","tokens_in":11676,"feed_emoji":"⚛️","tokens_out":13065,"duration_ms":112235,"temperature":0.7,"pith_summary":"The paper establishes empirical values for the Bohr-Weisskopf (BW) effect, the correction to hyperfine structure from the finite radial distribution of the nuclear magnetic moment, in ten mercury isotopes and isomers. It tests the Moskowitz-Lombardi rule, which relates the relative BW effect $\\epsilon$ to $1/|\\mu|$ plus an additive constant; for mercury this constant had never been measured. From the measured hyperfine splitting in muonic $^{199}$Hg, translated to electronic systems with screening factors, the paper finds $\\epsilon = -2.34(44)\\%$ for neutral $^{199}$Hg and $-2.59(49)\\%$ for $^{199}$Hg$^{+}$. Fitting the rule to all isotopes gives $c_{\\mathrm{exp}} = -3.6(8)\\times10^{-3}$ and $\\alpha_{\\mathrm{exp}} = 0.98(5)\\times10^{-2}\\,\\mu_N$, so the additive constant is significantly smaller in magnitude than the widely used Fujita-Arima value $c=-0.01$. If the result holds, precision atomic calculations for mercury can correct for the finite nuclear magnetization with an empirical anchor instead of a model spread.","feed_headline":"Mercury's hyperfine anomaly measured: rule constant near zero","feed_subtitle":"Fifty years untested; the new fit puts the constant near zero, sharpening mercury clocks and EDM searches.","key_machinery":"The central identity is the split $A = A_0 + A_{\\mathrm{BW}} + A_{\\mathrm{QED}}$, where $A_0$ is the point-nucleus hyperfine constant, $A_{\\mathrm{BW}}$ is the Bohr-Weisskopf correction, and $A_{\\mathrm{QED}}$ is the QED correction; subtracting well-calculated $A_0$ and $A_{\\mathrm{QED}}$ from the measured $A$ isolates $A_{\\mathrm{BW}}$. The second load-bearing mechanism is the two-step translation of the muonic-atom BW effect to a H-like ion and then to many-electron mercury through electronic screening factors $x_{\\mathrm{scr}} = \\epsilon_{\\mathrm{atom}}/\\epsilon_{\\mathrm{H-like}}$, computed at the random-phase-approximation-with-exchange level. The third is the use of measured differential anomalies between isotopes, with the ratio in Eq. (4) relating neutral and singly-ionized mercury, which spreads the absolute $^{199}$Hg value to other isotopes and isomers. These pieces work together because only $s$ and $p_{1/2}$ orbitals penetrate the nucleus and they are proportional there, so the H-like result carries all the BW information.","core_discovery":"On its own terms, the paper claims that the absolute BW effect in $^{199}$Hg can be extracted from experiment rather than left to nuclear models. The extraction starts from the measured 1s hyperfine constant of muonic $^{199}$Hg, where the muon is essentially unscreened, giving $\\epsilon = -68(8)\\%$ in the muonic atom. This is translated to the H-like ion and then to many-electron Hg through electronic screening factors, yielding $-2.59(49)\\%$ for the $6s$ state of $^{199}$Hg$^{+}$ and $-2.34(44)\\%$ for the $6s6p\\,{}^3P_1$ state of neutral $^{199}$Hg. An independent direct extraction from the measured $^{199}$Hg$^{+}$ hyperfine constant using all-orders atomic many-body calculations gives $-3.4(10)\\%$, consistent with the muonic route. Combining the $^{199}$Hg anchor with measured differential anomalies gives BW effects for ten isotopes and isomers, and fitting $\\epsilon = c - \\alpha/|\\mu|$ to these data yields $c_{\\mathrm{exp}} = -3.6(8)\\times10^{-3}$ and $\\alpha_{\\mathrm{exp}} = 0.98(5)\\times10^{-2}\\,\\mu_N$ for neutral mercury. The paper concludes that the additive constant is much smaller than the Fujita-Arima value $c_{\\mathrm{FA}}=-0.01$, and that the discrepancy points to omitted many-body nuclear effects in the older theory.","pith_inferences":["Beyond the paper's claims, the agreement between the muonic and direct routes suggests that the same two-step translation could be applied to muonic-atom data for other heavy nuclei to test whether the Moskowitz-Lombardi rule's constants are universal or nucleus-specific.","A published numerical value for the differential Breit-Rosenthal factor $^{199}\\delta^{201}$ in the $6s6p\\,{}^3P_1$ state would close the only unshown check in the neutral-mercury chain; until then, the neutral screening factor carries a hidden systematic absent from the ion-only extraction.","If $c_{\\mathrm{exp}}$ is genuinely near zero, the Fujita-Arima factorization misses substantial many-body nuclear effects, which implies that absolute-anomaly measurements near other shell closures, not just isotopic differential anomalies, are the right way to map where the rule holds."],"forward_implications":["The absolute BW correction for $^{199}$Hg is now anchored empirically at the percent level, so hyperfine-structure calculations can use $-2.34(44)\\%$ (neutral) or $-2.59(49)\\%$ (ion) instead of a spread of nuclear-model predictions.","The fitted coefficients give a concrete empirical rule for neutral mercury, $\\epsilon(\\%) = -0.36 - 0.98/|\\mu(\\mu_N)|$, which predicts both absolute anomalies and isotopic trends.","Using the old Fujita-Arima constant $c=-0.01$ would make the predicted absolute anomaly for $^{199}$Hg about 45% larger in magnitude than the empirical value.","The ten isotope and isomer BW values provide direct benchmarks that microscopic nuclear-magnetization models must reproduce, which is relevant to Schiff-moment and electric-dipole-moment interpretations in mercury."],"supporting_citations":[{"why":"Defines the Moskowitz-Lombardi rule $\\epsilon=\\alpha_{\\mathrm{ML}}/|\\mu|$ that the paper tests.","marker":"[37]"},{"why":"Supplies the extended form $\\epsilon=c-\\alpha/|\\mu|$ and the Fujita-Arima constants $c=-0.01$, $\\alpha=1.2\\times10^{-2}\\,\\mu_N$ that the fit overturns.","marker":"[32]"},{"why":"Provides the measured muonic-199Hg 1s hyperfine constant used as the experimental seed.","marker":"[33]"},{"why":"Supplies the two-step method for translating the muonic-atom BW effect to H-like and many-electron systems.","marker":"[26]"},{"why":"Provides the single-particle model and five nuclear wave-function shapes used to propagate nuclear-model uncertainty from muonic to H-like Hg.","marker":"[56]"},{"why":"Supplies the electronic screening factors $x_{\\mathrm{scr}}$ that carry the H-like result into Hg$^{+}$ and neutral Hg.","marker":"[23]"},{"why":"Provides the measured differential anomalies $\\Delta_{199}$ for mercury isotopes used to extend the absolute 199Hg value to other isotopes.","marker":"[60]"},{"why":"Provides the measured 199Hg$^{+}$ hyperfine constant used for the independent direct extraction.","marker":"[43]"},{"why":"Supplies QED corrections for Hg$^{+}$ interpolated from Ba$^{+}$ and Ra$^{+}$, and benchmarking for the many-body calculation.","marker":"[15]"},{"why":"Provides the all-orders correlation-potential method used to compute the point-nucleus hyperfine constant for the direct check.","marker":"[64]"}],"fun_headline_variants":["Mercury hyperfine test: rule constant near zero after 50 years","Bohr-Weisskopf effect in mercury extracted from muonic data","Mercury's nuclear size effect measured: Moskowitz-Lombardi update","New constant for mercury hyperfine anomaly challenges old theory"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the measured ratio of differential anomalies between isotopes 199 and 201 in neutral mercury and in the ion equals the ratio of the Bohr-Weisskopf effects alone, with the small difference in how the finite nuclear charge distribution shifts the electron wave functions (the differential Breit-Rosenthal effect) neglected; the paper says it checked this numerically but does not show the check.","fun_headline_variants_meta":{"raw":{"variants":["Mercury hyperfine test: rule constant near zero after 50 years","Bohr-Weisskopf effect in mercury extracted from muonic data","Mercury's nuclear size effect measured: Moskowitz-Lombardi update","New constant for mercury hyperfine anomaly challenges old theory"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000233,"raw_usage":{"total_tokens":1551,"prompt_tokens":1064,"completion_tokens":487,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":680,"completion_tokens_details":{"reasoning_tokens":410}},"tokens_in":680,"tokens_out":487,"duration_ms":5038,"temperature":1.0,"reasoning_tokens":410,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T20:11:07.763839+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Calculate the differential Breit-Rosenthal contribution $^{199}\\delta^{201}$ for the $6s6p\\,{}^3P_1$ state of neutral Hg with a relativistic many-body method; if it is comparable to the 0.4% uncertainty in the ratio 0.902(4), the neutral-atom BW values and the fitted constant $c_{\\mathrm{exp}}$ shift by more than their quoted errors. Alternatively, a new measurement of the muonic $^{199}$Hg 1s hyperfine splitting with uncertainty well below 0.12 keV would directly test the seed value $-68(8)\\%$.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the Moskowitz-Lombardi rule $\\epsilon=\\alpha_{\\mathrm{ML}}/|\\mu|$ that the paper tests."},{"cited_title":"Fujita and A","cited_arxiv_id":null,"evidence_quote":"Supplies the extended form $\\epsilon=c-\\alpha/|\\mu|$ and the Fujita-Arima constants $c=-0.01$, $\\alpha=1.2\\times10^{-2}\\,\\mu_N$ that the fit overturns."},{"cited_title":"B¨ uttgenbach, Hyperfine Interact.20, 1 (1984)","cited_arxiv_id":null,"evidence_quote":"Provides the measured muonic-199Hg 1s hyperfine constant used as the experimental seed."},{"cited_title":"Sanamyan, B","cited_arxiv_id":null,"evidence_quote":"Supplies the two-step method for translating the muonic-atom BW effect to H-like and many-electron systems."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the single-particle model and five nuclear wave-function shapes used to propagate nuclear-model uncertainty from muonic to H-like Hg."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the electronic screening factors $x_{\\mathrm{scr}}$ that carry the H-like result into Hg$^{+}$ and neutral Hg."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the measured differential anomalies $\\Delta_{199}$ for mercury isotopes used to extend the absolute 199Hg value to other isotopes."},{"cited_title":"Engel, M","cited_arxiv_id":null,"evidence_quote":"Provides the measured 199Hg$^{+}$ hyperfine constant used for the independent direct extraction."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies QED corrections for Hg$^{+}$ interpolated from Ba$^{+}$ and Ra$^{+}$, and benchmarking for the many-body calculation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the all-orders correlation-potential method used to compute the point-nucleus hyperfine constant for the direct check."}],"review_version":1}