{"id":"bee7bcfd-8576-43dc-801d-63f67e35c0ef","arxiv_id":"2607.24634","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Recommended B(T) for N2 (20-3000 K) and O2 (20-2000 K) are derived from ab initio pair potentials tuned to reanalyzed experimental data, with quantum effects included via path-integral Monte Carlo.","lead":"This paper presents new recommended values for the second virial coefficient of nitrogen and oxygen over wide temperature ranges, based on tuning quantum-chemical pair potentials to high-precision experimental data. A smart generalist should care because these coefficients are needed for precise metrology of air, pressure standards, and flow measurement.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Nitrogen B(T) may inherit bias from unpublished higher-virial inputs used in isotherm reanalysis.","rationale":"The reader's weakest assumption identified exactly the load-bearing concern: the nitrogen B values used for tuning are extracted from isotherms analyzed with higher virial coefficients from a preliminary, partly unpublished theory. The paper itself flags the three-body PES as available only on request (Sec. 10), and Sec. 2.1 reveals that the D, E, F values came from a version of the pair PES tuned to Egan and Yang (scaling 0.52) rather than the final retuned potential (scaling 0.472). Because the final tuning target is the B set produced under those constraints, any bias in D, E, or F is inherited by the recommended B(T). The 2000 cm9 mol-3 uncertainty assigned to D is an estimate without a documented derivation, and the paper's agreement checks are not independent. This concern is concrete and testable, and it is not mitigated by the paper's existing uncertainty propagation because the propagated component only covers the assumed D-uncertainty rather than a possible larger systematic error. The conditional verdict already captures this risk, since it asks for release of the three-body PES and further checks, so I recommend leaving the verdict unchanged. The acoustic virial discrepancies are a secondary but related indication that the uncertainty estimates may be optimistic; however, the D/E/F circularity is more central to the novelty of the nitrogen analysis.","tokens_in":17600,"tokens_out":6164,"duration_ms":54580,"concrete_test":"Compute D, E, and F for N2 using the final retuned pair PES and the nonadditive three-body PES at the Table 1 temperatures; reanalyze the Nowak and McLinden isotherms with these values instead of the preliminary ones. If the extracted B differs by more than the quoted U(B) of 0.05 cm3/mol over 150-340 K, the recommended B(T) carries a systematic bias not included in the uncertainty estimate. A minimal auxiliary check is to quantify |D_final - D_preliminary| against the assumed 2000 cm9 mol-3 uncertainty.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim for nitrogen rests on B values reanalyzed from density isotherms by fixing D, E, and F to values computed with a preliminary pair PES (scaling 0.52) and an unpublished nonadditive three-body PES (Sec. 2.1, Sec. 10). The pair potential is then retuned (scaling 0.472) to match these B values, so any systematic error in D, E, or F is absorbed into the extracted B and propagates into the final recommended B(T). The paper estimates the uncertainty of D as 2000 cm9 mol-3 and propagates it, but this is a heuristic: it is not derived from a comparison of D computed with the preliminary and final potentials, and the three-body PES is 'available upon reasonable request' rather than public. Because the final B(T) is tuned to data that were themselves corrected using this unpublished theory, the internal agreement shown in Fig. 1 cannot validate D, E, and F; agreement of Ccalc with Cexp is also not independent since Cexp comes from the same constrained fits. A bias in D larger than the estimate would shift the recommended B by roughly delta_D * rho_max^2, which at 300 K and 12 MPa gives about 0.05 cm3/mol for delta_D=2000, comparable to the quoted U(B); a larger systematic bias would exceed the quoted uncertainty and the claim of uncertainty similar to the best experiments would fail. The paper's own Sec. 4.1 notes unexplained deviations from acoustic virial data, an independent check that is not incorporated into the uncertainty estimate.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper recommends values for the second virial coefficient B(T) and its expanded uncertainty for molecular nitrogen (20-3000 K) and oxygen (20-2000 K). The recommendations are obtained by tuning previously published ab initio pair potentials to selected high-accuracy experimental data, with full quantum treatment via path-integral Monte Carlo. For nitrogen, the authors reanalyze density isotherms with the third through sixth virial coefficients constrained to values calculated using a preliminary pair PES and a new nonadditive three-body PES, and then retune the pair PES (scaling factor 0.472) to the extracted B values. For oxygen, the three singlet/triplet/quintet PES are retuned by shifting the isotropic C8 dispersion coefficient so that B at 293.15 K matches the refractometry data of Egan. The paper provides smooth correlations for B(T) and U(B), tables of parameters, PIMC data, and Fortran code for the tuned potentials.","tokens_in":18007,"tokens_out":4890,"duration_ms":42486,"significance":"If the recommended values are reliable, they are directly useful for pressure metrology, buoyancy corrections, humidity metrology, and future reference equations of state for air. The paper has clear strengths: the PIMC treatment of quantum effects is the appropriate rigorous method, the QFH semiclassical calculations provide a useful internal cross-check, the uncertainty estimation is transparent, and the authors provide the tuned potential routines and PIMC data as supplementary material. The central concern, discussed below, is that for nitrogen the B values used for tuning were extracted from isotherms constrained by higher virial coefficients from a preliminary and partly unpublished theory, so the final recommendation may inherit a systematic bias that is not fully captured by the quoted uncertainty.","major_comments":[{"comment":"The B values used to retune the nitrogen pair potential are extracted from isotherms with D, E, and F constrained to values calculated with a preliminary pair PES (scaling 0.52) and an unpublished nonadditive three-body PES. The final pair PES is then tuned to those extracted B values (scaling 0.472). Consequently, any systematic error in the constrained higher virial coefficients propagates directly into the extracted B values and into the final recommended B(T), and the internal agreement shown in Fig. 1 cannot validate D, E, and F. The paper's 2000 cm9 mol-3 uncertainty estimate for D is a heuristic: it is not derived from a comparison of D computed with the preliminary and final pair potentials, and no uncertainty is assigned to E and F. Table 1 reinforces this concern: at 130 K, Ccalc-Cexp = -48 cm6 mol-2 exceeds the reported U(Cexp) = 35 cm6 mol-2, indicating that the constrained fits may still contain systematic bias in the strongly correlated pair (B, C). I recommend that the authors publish the three-body PES, recompute D, E, and F with the final pair potential as well as the preliminary one, and perform sensitivity tests in which D (and, at least approximately, E and F) are varied within physically plausible bounds; the resulting effect on B should be incorporated into U(B).","section":"Section 2.1 and Table 1"},{"comment":"The acoustic second virial coefficient data of Ewing and Trusler and of Boyes deviate from the recommended correlation in a way that the authors do not explain quantitatively. Because these data provide an independent check that was not used in the tuning, the claim that U(B) = 0.05 cm3 mol-1 for nitrogen in the range 150-500 K is not supported by the full body of evidence. The authors should either reconcile the acoustic data with the correlation within a defensible uncertainty budget or quote a larger uncertainty that encompasses the discrepancy.","section":"Section 4.1 and Figure 2"},{"comment":"The low-temperature uncertainty estimate for nitrogen is obtained by rescaling the difference between the original and retuned potentials so that it equals 0.05 cm3 mol-1 at 150 K; this is a calibration to the tuning data rather than an independent bound on the potential error. The high-temperature estimate, a linear growth to 0.3 cm3 mol-1 at 3000 K, is justified by a statement about molecular flexibility effects in H2, which is a different molecule with different rotational constants. The authors should support these estimates with actual calculations of the flexibility correction for N2 and O2, or with a direct estimate of the rigid-rotor truncation error.","section":"Section 5"}],"minor_comments":[{"comment":"The procedure for assigning the plateau ranges in the virial fits is described as subjective; please provide the assigned minimum and maximum pressures for each isotherm in the supplementary material so that the analysis is reproducible.","section":"Section 2.1"},{"comment":"The statement that the PIMC calculations have an accuracy of 10 ppm would be more convincing if convergence data or a convergence plot were included in the supplementary material.","section":"Section 3.3"},{"comment":"Data availability via 'reasonable request' is not sufficient for a recommended reference dataset; the nonadditive three-body PES for N2 should be deposited in a public repository.","section":"Section 10"},{"comment":"The notation in Eq. (13) is potentially confusing because B appears both as the second density virial coefficient and as a refractivity-related coefficient; please add a clarifying phrase or a different symbol for the latter.","section":"Eq. (13)"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and the underlying approach is sound in principle. My recommendation of major revision is driven by the nitrogen analysis: the central B values used for tuning depend on higher virial coefficients from a preliminary and partly unpublished theory, so the reproducibility and the uncertainty claim both need strengthening. If the authors provide the three-body PES, quantify the sensitivity of B to D/E/F, and address the acoustic virial discrepancy, I would support acceptance."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This paper gives recommended B(T) correlations for nitrogen (20–3000 K) and oxygen (20–2000 K) with uncertainty estimates, built by retuning previously published ab initio pair potentials against selected high-accuracy experiments. That is genuinely useful: metrology applications (air density, refractometry-based pressure standards, flow meters) currently rely on decades-old correlations, and the oxygen case exploits new Egan refractometry data with 0.1 cm^3/mol uncertainty. The quantum treatment via PIMC is state of the art, and the authors are transparent about how uncertainties are constructed (heuristic in places, but clearly described).\n\nThe nitrogen part deserves scrutiny. The reanalysis of the Nowak and McLinden isotherms fixes D, E, F to values computed from a preliminary pair potential (scaling 0.52) and an unpublished three-body PES. Those constrained B values then tune the final pair potential (scaling 0.472). So the \"experimental\" B for nitrogen is partly theory; if D is biased, the bias leaks into the recommended B. The paper propagates an estimated uncertainty in D of 2000 cm^9/mol^3, which is a reasonable heuristic, but it's not derived from a comparison of D from the preliminary and final potentials, and the three-body PES is only available on request. The magnitude of the risk is real: at 300 K and 12 MPa, a 2000 cm^9/mol^3 shift in D moves B by about 0.05 cm^3/mol, which is comparable to the quoted U(B). So the uncertainty claim could be optimistic if the theory bias is larger than the estimate.\n\nThe acoustic virial data (Ewing & Trusler, Boyes) show deviations that the authors themselves call \"disappointing\" and can't explain. They don't fold this into the uncertainty. That is an honest statement, but it means the uncertainty budget is not conservative.\n\nThe oxygen part is cleaner: tuning against Egan's refractivity data, with a single parameter per surface and no entangled higher-virial extraction. I'd trust that one more.\n\nOverall: the paper is a solid contribution, clearly written, with reproducible potential codes and B(T) data in the supplementary material. The nitrogen circularity is a legitimate concern, but it's not disqualifying—the authors identify it and make a reasonable, if heuristic, uncertainty allowance. I'd like to see the three-body PES released and a more explicit treatment of the acoustic discrepancies, but the recommendations are probably better than anything currently available for these fluids.\n\nThis deserves a serious referee. The reader's conditional verdict is about right; I'd lean toward accepting after minor-to-moderate revision, with the nitrogen inputs and uncertainties revisited.","headline":"A careful, useful retuning of B(T) for N2 and O2, with a real reproducibility catch in the nitrogen higher-virial inputs.","tokens_in":18406,"tokens_out":1926,"would_cite":true,"duration_ms":16308,"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":"Retuned ab initio pair potentials, constrained by reanalyzed density and refractometry data, produce recommended second virial coefficients for nitrogen and oxygen with experimental-quality uncertainty across 20–3000 K.","keywords":["second virial coefficient","nitrogen","oxygen","ab initio pair potential","path-integral Monte Carlo","virial coefficient reanalysis","refractometry","thermophysical properties"],"falsifier":"Measure $B(T)$ for nitrogen between 130 and 340 K with a technique that does not rely on calculated higher virial coefficients, such as low-density speed-of-sound or dielectric-constant gas thermometry, and compare with the recommended curve; a deviation outside the claimed expanded uncertainty of 0.05 cm$^3$ mol$^{-1}$ over 150–500 K would show the higher-virial constraint is biased.","tokens_in":17337,"feed_emoji":"🌡️","tokens_out":14489,"duration_ms":112101,"temperature":0.7,"pith_summary":"The paper establishes recommended values for the second virial coefficient $B(T)$—the leading measure of how much a gas's density departs from ideal-gas behavior—for nitrogen over 20–3000 K and oxygen over 20–2000 K, each accompanied by an expanded-uncertainty curve. The route combines previously published ab initio pair potentials with selected high-accuracy experiments: the nitrogen potential is tuned to reanalyzed density data whose uncertainty is sharply reduced by constraining the fourth through sixth virial coefficients to calculated values, while the oxygen potentials are tuned to new refractometry data. Quantum effects are included exactly with path-integral Monte Carlo, so the curves are valid over a temperature range far wider than any single experiment covers. If the recommendations hold, they give metrology, atmospheric science, and reference equations of state a firmer basis for the nonideal behavior of air's two main components.","feed_headline":"Recommended virial coefficients now span 20–3000 K for N2 and O2","feed_subtitle":"Quantum-tuned potentials and reanalyzed density data give experimental-quality uncertainties over a much wider range.","key_machinery":"The load-bearing object is the two-body potential-energy surface (pair PES): an ab initio rigid-rotor interaction potential for N$_2$ and three spin-coupled surfaces (quintet, triplet, singlet) for O$_2$, each modified by one physical parameter so that the calculated $B(T)$ matches selected experiments. The $B(T)$ values are then computed with path-integral Monte Carlo (PIMC), which represents each quantum rotor pair as ring polymers and evaluates the Mayer-$f$ average exactly, including nuclear quantum effects; a quadratic Feynman–Hibbs (QFH) approximation supplies the higher virial coefficients $C$–$F$ used to constrain the nitrogen density isotherms and confirms the PIMC results at most temperatures. The recommended values are carried by the correlation of Eq. (15) with the parameters in Table 2, and the low-temperature uncertainties by Eq. (16).","core_discovery":"The central claim is that slightly retuning previously published ab initio pair potentials—scaling the difference between the CCSDT(Q) and CCSD(T) levels by 0.472 for nitrogen, and adjusting the isotropic $C_8$ dispersion coefficient of each of the three spin-coupled oxygen surfaces—makes the $B(T)$ calculated from the potentials agree with selected high-accuracy measurements while retaining a temperature range far wider than any experiment. For nitrogen, the most influential measurements are values of $B$ obtained by reanalyzing density isotherms with the fourth through sixth virial coefficients fixed to semiclassically calculated values; for oxygen, the tuning target is a single refractometry data set. The paper reports the recommended $B(T)$ as smooth correlations (Eq. (15) with Table 2) for 20–3000 K (N$_2$) and 20–2000 K (O$_2$), with expanded ($k=2$) uncertainties of 0.05 cm$^3$ mol$^{-1}$ for nitrogen over 150–500 K and 0.1 cm$^3$ mol$^{-1}$ for oxygen over 300–400 K, rising linearly outside those ranges. Quantum effects are included exactly with path-integral Monte Carlo, and the semiclassical quadratic Feynman–Hibbs results agree closely except below about 50–60 K.","pith_inferences":["An independent, openly available three-body surface for nitrogen would let the $D,E,F$ constraint be tested directly; a shift in $D$ beyond the assumed 2000 cm$^9$ mol$^{-3}$ would move the extracted $B$ values by a comparable amount.","The paper's explanation of the refractometry offset implies that a future ab initio calculation of the second refractivity virial coefficient $B_R$ could adjudicate between the recommended $B(T)$ and the values derived from refractivity.","No independent low-temperature oxygen data exist; the rescaled-tuning-difference method predicts larger low-temperature uncertainty for oxygen than nitrogen, so new measurements below 300 K would be the decisive test.","The per-spin-surface $B$ shifts in oxygen stay nearly equal up to 373 K, suggesting that the one-parameter tuning result does not depend on which detailed part of the potential is adjusted; the same strategy should transfer to other spin-coupled molecular pairs."],"forward_implications":["The recommended $B(T)$ curves provide low-density boundary conditions for updating the nitrogen reference equation of state (2000) and the oxygen reference equation of state (1985), which the paper identifies as outdated.","Applications in flow metering, refractometry-based pressure standards, and humidity or length metrology can use air's main components with uncertainties close to the best single measurements over a temperature span no experiment covers.","Viscosity and thermal conductivity values from the original potentials remain valid, because the tuning shifts them by far less than their uncertainties.","The tuned potentials form the foundation for a revised second virial coefficient of dry air, once the N$_2$–O$_2$ cross coefficient and the small argon contributions are added.","The unexplained discrepancies with existing acoustic virial data for nitrogen show that new acoustic measurements are needed to test and refine these recommendations."],"supporting_citations":[{"why":"Supplies the nitrogen density isotherms whose reanalysis, constrained by calculated higher virial coefficients, yields the low-uncertainty B values that dominate the tuning.","marker":"[19]"},{"why":"Supplies additional nitrogen density isotherms that validate the reanalysis at 293 K and 340 K and extend tuned B values to 480 K.","marker":"[21]"},{"why":"Supplies the refractometry-based nitrogen B(T) over 293–433 K used as the initial tuning target and as a consistency check.","marker":"[20]"},{"why":"Supplies the refractometry-based oxygen B(T) over 293–373 K with 0.1 cm^3 mol^-1 uncertainty that drives the retuning of all three oxygen surfaces.","marker":"[6]"},{"why":"Provides the ab initio nitrogen pair potential whose CCSDT(Q)-minus-CCSD(T) contribution is rescaled to tune B(T).","marker":"[18]"},{"why":"Provides the three ab initio oxygen pair potentials whose isotropic C8 dispersion coefficients are adjusted.","marker":"[23]"},{"why":"Supplies the Feynman–Hibbs path-integral formalism and the quadratic Feynman–Hibbs approximation used for the semiclassical higher virial coefficients.","marker":"[25]"},{"why":"Establishes the path-integral Monte Carlo route for computing molecular second virial coefficients from a six-dimensional potential.","marker":"[35]"}],"fun_headline_variants":["Tuned ab initio potentials yield virial coefficients for N2 (20–3000 K) and O2 (20–2000 K)","Quantum-tuned virial coefficients for N2 and O2: experimental accuracy from 20 K upward","Wide-range B(T) for N2 and O2: tuned potentials, matched to best experiments","Recommended virial data for N2 and O2: quantum effects, wide temperature span"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The nitrogen recommendations stand on the accuracy of the calculated fourth-through-sixth virial coefficients—theory values from a preliminary pair potential and a nonpublic three-body surface—used to constrain the reanalysis of density measurements; if those values are biased, the extracted $B$ values and the tuned potential inherit the bias.","fun_headline_variants_meta":{"raw":{"variants":["Tuned ab initio potentials yield virial coefficients for N2 (20–3000 K) and O2 (20–2000 K)","Quantum-tuned virial coefficients for N2 and O2: experimental accuracy from 20 K upward","Wide-range B(T) for N2 and O2: tuned potentials, matched to best experiments","Recommended virial data for N2 and O2: quantum effects, wide temperature span"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001246,"raw_usage":{"total_tokens":5121,"prompt_tokens":964,"completion_tokens":4157,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":580,"completion_tokens_details":{"reasoning_tokens":4048}},"tokens_in":580,"tokens_out":4157,"duration_ms":27729,"temperature":1.0,"reasoning_tokens":4048,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T15:25:54.184889+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure $B(T)$ for nitrogen between 130 and 340 K with a technique that does not rely on calculated higher virial coefficients, such as low-density speed-of-sound or dielectric-constant gas thermometry, and compare with the recommended curve; a deviation outside the claimed expanded uncertainty of 0.05 cm$^3$ mol$^{-1}$ over 150–500 K would show the higher-virial constraint is biased.","supporting_citations":[],"review_version":2}