{"id":"ebe24b7b-fcd4-4f8d-818d-fd4f78b96f15","arxiv_id":"2604.00794","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Apparent liquid-liquid transition signatures in supercooled water coincide with dynamical arrest, implying a glass transition near 189 K rather than a true second liquid.","lead":"Simulations suggest that what looks like a transition between two liquid forms of supercooled water is actually the liquid freezing into a glass. If correct, the long-sought second-liquid transition may instead be a glass transition near 189 K.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The inferred Tg=189±8 K depends on VFT extrapolation of relaxation times that may already include non-ergodic/trapping signatures; the data range is too short to discriminate functional forms, so the central quantitative claim is not robust.","rationale":"The paper presents compelling direct evidence—rate-dependent depletion in Fig. 2A/D, aging after rapid ramps, and the low-mobility plateau in Fig. 1—that the apparent low-density state is kinetically arrested rather than an equilibrium liquid on simulation timescales. This makes the qualitative reinterpretation plausible. However, the link to the quantitative Tg=189±8 K, which is crucial for the claim that the putative LLT coincides with a glass transition, depends on an extrapolation that is not tested. The VFT fit uses data in a window where the system may already be entering the non-ergodic regime, and alternative functional forms would alter the extrapolated Tg. The reader's weakest assumption is exactly this: the VFT extrapolation. I agree. The recommended verdict remains CONDITIONAL: accept the qualitative glassy-arrest scenario as a plausible reinterpretation, but with reservations about the specific Tg value and its coincidence with the LLT temperature. My proposed test would strengthen or weaken the quantitative claim.","tokens_in":15797,"tokens_out":11148,"duration_ms":113596,"concrete_test":"Re-fit Eq. (2) to the 1-bar relaxation times using only state points where the instantaneous MSD shows no arrest plateau (i.e., no trapping events longer than ~10% of trajectory). Then re-fit the same data with (i) the Bässler form tau=tau0 exp((T0/T)^2) and (ii) a VFT-plus-Arrhenius crossover model. If the resulting Tg at tau=100 s shifts by more than ±10 K from 189 K, the reported numeric claim is not robust; if it remains within ±8 K, the VFT extrapolation is validated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim—ambient-pressure Tg=189±8 K, close to the putative LLT temperature—rests on fitting the VFT equation (Eq. 2) to dielectric relaxation times that span only ~10^-9 to 10^-6 s at 1 bar (Fig. 3C) and then extrapolating eight orders of magnitude to tau=100 s. The paper's own observations make this extrapolation especially fragile. At the lowest simulated temperatures, the system is not a homogeneous equilibrium liquid: the trMSD traces (Fig. 1) show long-lived plateaus in the low-density state, and the two-state fluctuations are sensitive to barostat parameters (Fig. S1), indicating that the simulated dynamics include non-ergodic trapping. If any of the relaxation times used in the VFT fit originate from state points where the dipole autocorrelation has not fully decayed (or where the system intermittently visits an arrested state), the fitted VFT parameters are biased. Moreover, the reported ±8 K is the spread between three models, not an error bar on the extrapolation; alternative functional forms (e.g., Bässler, Arrhenius crossover) that fit the same short-time data equally well are common for glass formers and would shift Tg substantially. Because the entire argument that the LLT region coincides with vitrification depends on Tg landing between the simulated window and the experimental 100-s definition, the unvalidated VFT extrapolation is the most load-bearing assumption.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper argues that the two-state fluctuations widely interpreted as evidence of a liquid–liquid transition (LLT) in supercooled water are instead manifestations of dynamical arrest: the low-density state is a kinetically trapped glass, not an equilibrium second liquid. Evidence includes time-resolved mean-square displacement (trMSD) plateaus in three water models (DNN@SCAN, DNN@MB-pol, TIP4P/2005), rate-dependent vitrification in pressure- and temperature-ramp simulations, barostat sensitivity of the two-state fluctuations (SM Fig. S1), and the coincidence of kinetic arrest with the heat-capacity maximum. The central quantitative claim is an ambient-pressure glass-transition temperature of 189 ± 8 K, obtained by VFT extrapolation of simulated structural relaxation times to τ = 100 s, close to the temperature range where an LLT has been proposed. The authors also construct isochronal lines in the P–T plane for the DNN@SCAN model and compare simulated relaxation times with dielectric and NMR experiments at 1 bar and 5 kbar.","tokens_in":16194,"tokens_out":1925,"duration_ms":22042,"significance":"If correct, the reinterpretation is significant: it would unify several experimental and simulation observations under a glass-transition framework, reducing the need for an equilibrium LLT to explain water's anomalies. The paper has notable strengths: it uses three independent models, including two advanced machine-learning potentials with explicit melting-point corrections; it directly compares with experimental dielectric, NMR, calorimetric, and X-ray data; it provides a concrete, testable prediction (a high Tg near 189 K for bulk low-density water); and it includes supplementary mechanistic checks such as barostat sensitivity and aging behavior. The manuscript is also unusually candid about its own limitations, including the extrapolative nature of the VFT fit and the small-system dependence of the apparent phase fluctuations.","major_comments":[{"comment":"The headline Tg = 189 ± 8 K rests on VFT extrapolation of simulated relaxation times spanning roughly 10^-9 to 10^-6 s (Fig. 3C) to τ = 100 s, an extrapolation of eight orders of magnitude. The reported ±8 K is the spread among three models, not an uncertainty on the extrapolation itself. The paper acknowledges that the open symbol in Fig. 3B is a lower bound from a partially decayed correlation function, and that low-temperature simulations may be non-ergodic. This is load-bearing because the central claim — that vitrification coincides with the proposed LLT — depends on Tg landing near 189 K. A concrete robustness test would be to report VFT fit parameters with their confidence intervals, to show χ² residuals for alternative forms (e.g., Bässler, modified VFT, or Arrhenius crossover), and to test sensitivity of Tg to excluding the lowest-temperature points or to alternative definitions","section":"§4.4 / Fig. 3C, Eq. (2)"},{"comment":"The heat-capacity peak in Fig. 2E is obtained by differentiating a six-parameter empirical fit (Eq. S2) to enthalpy data from a single cooling trajectory at 2 K/ns. The fit parameters are smoothing parameters, not physical ones, and no uncertainty is provided. The claim that the Cp peak at ~243 K 'approaches' the experimental maximum at 228.9 K is presented without quantitative error bars or a rate-dependence scaling argument. This comparison is central to the argument that vitrification masks a thermodynamic divergence, so the lack of uncertainty propagation is a weakness. A simple bootstrap or multiple-trajectory estimate of the peak temperature would strengthen the claim.","section":"Fig. 2E and §3.2 (SM)"},{"comment":"The isochronal 'phase diagram' in Fig. 4 uses VFT extrapolations (open circles) to draw isochrones out to τ = 100 s, yet the VFT parameters at each pressure are not shown, nor are the number of state points per pressure or the quality of the fits. The claim that T_max∇Pτ = 193.5 K coincides with the proposed LLT is made from this extrapolated map; without fit-quality metrics or error bars on the isochrone positions, this coincidence is not quantitatively demonstrated. Providing the fit parameters and uncertainties, or at least a table of Tg(P) values, would make the map reproducible and the claim testable.","section":"Discussion, Fig. 4"},{"comment":"The barostat-sensitivity result (SM Fig. S1) is a strong point: it shows that the two-state fluctuations can be eliminated by changing the damping constants. However, the claim that the original-settings trajectory 'again exhibits trapping' in the low-energy state for ~0.5 μs is based on a single trajectory at 192 molecules. Given that the authors themselves argue that small boxes amplify pressure fluctuations, the generality of the trapping phenomenon should be quantified — e.g., how many independent trajectories show how long a plateau, or how the trMSD plateau duration varies with box size. If this is not presented, the claim that apparent LLT fluctuations are generally a glass-artifact remains suggestive rather than established.","section":"Fig. 1 / SM §2.1"}],"minor_comments":[{"comment":"The notation 'DNN@SCAN' and 'DNN@MB-pol' is used inconsistently: the models are also referred to as 'SCAN@DNN' in SM §2.1. Please standardize.","section":"Throughout"},{"comment":"The legend says 'experiment: NMR, 2.8nm; DS, 1.9nm; DS, LiCl-solution' but the caption does not specify which symbols correspond to which experimental dataset, nor the temperature ranges where each dataset is valid. Adding a legend entry with symbols and a sentence about the confinement sizes would improve readability.","section":"Fig. 3C legend"},{"comment":"The empirical pressure fit τ(P) = a1 exp(-P/P0) + a2 + a3 P^a4 uses five parameters for three temperatures; the exponent a4 is not reported for any temperature, and the fit is stated to be 'just a guide to the eye.' Since the open symbol at 2000 bar is obtained via frequency-pressure superposition of a partially decayed spectrum, the reader should be told explicitly whether the open symbol is included in the VFT fit at 1 bar or only in the pressure fit. This affects the interpretation of the low-pressure extrapolation.","section":"Eq. (S3)"},{"comment":"The inset of Fig. 2A shows physical aging after rapid depressurization, but the axis label 'Epot / kJ/mol' is ambiguous: is this per molecule or per mole? In the main panels it is stated to be per molecule; please clarify.","section":"Fig. 2A inset"},{"comment":"References 6, 7, 16, 36, 37, 46, 47 include years 2025–2026, which are plausible for a 2026 arXiv posting. However, for citations of unpublished results (e.g., Ref. 16 'Lunkenheimer, Reuter, Schulz, Wolf, Loidl, Phys. Rev. E 111(6), 065408 (2025)') the journal and year are given; please verify page/article numbers for all preprints, as some cited works appear to be very recent and not yet indexed.","section":"References"},{"comment":"The statement that 'many liquids follow the so-called two-thirds rule' is used to support Tg ≈ 182 K, but the rule is empirical and has exceptions; this is fine as a consistency check, but should be clearly labeled as 'not definitive' (which it already is, to the authors' credit). Perhaps add a sentence that the two-thirds rule for water's melting point at 273 K gives 182 K, which is within the stated uncertainty.","section":"SM §4.4"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is conceptually interesting and well-written, and I want to see it revised rather than rejected. The main blocker is the quantitative Tg extrapolation: the reported ±8 K is a model-spread, not an extrapolation error, and the VFT fit parameters and alternative functional-form tests are missing. This is fixable within the manuscript's scope — adding fit diagnostics, confidence intervals, and robustness checks does not require new physics. The second issue is the Cp comparison in Fig. 2E, which needs uncertainty quantification. The barostat-sensitivity result is a strength, but its single-trajectory basis should be acknowledged more explicitly. If the robustness checks are added, I would be comfortable with acceptance; if the authors cannot provide them, the headline Tg number should be downgraded to a qualitative range."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear Colleague,\n\nYou should know upfront: this is a paper that will get people talking. Pabst and Hassanali use three water models to argue that the supposed liquid-liquid transition in supercooled water is actually a glass transition—the low-density state is not a liquid at all but an amorphous solid that intermittently traps the system. The qualitative case is strong, and they support it with multiple independent probes: trMSD plateaus, pressure- and temperature-ramp vitrification, barostat sensitivity, and comparisons to X-ray scattering and heat capacity experiments. This is genuinely new as a systematic demonstration across ML models, and it reframes the long debate in a productive way.\n\nWhere I get uneasy is the headline number. The claim that ambient-pressure Tg = 189 ± 8 K comes from fitting VFT to dielectric relaxation times that span only ~10^-9 to 10^-6 s, then extrapolating eight orders of magnitude to tau = 100 s. The ±8 K is the spread between models, not a fit uncertainty, and no error bars are reported on the VFT parameters. At the lowest simulated temperatures the system is already close to arrest, so some of the points used in the fit may be biased by non-ergodicity. Alternative functional forms (Bässler, etc.) would give a different Tg. This is a real soft spot, and it's the load-bearing part of the quantitative conclusion. The paper would be much stronger if they released the tau(T) data and showed a stability analysis of Tg against fit choices and data range.\n\nThat said, the qualitative message does not depend on the exact Tg. The observations that two-state fluctuations coincide with kinetic arrest, that they vanish under different barostat settings, and that larger simulations don't show them are convincing. The authors are also candid about limits: they explicitly say they cannot exclude an LLT, and they note that their Tg estimate is model-dependent. That honesty helps.\n\nI would send this to peer review. A serious referee should push on the VFT extrapolation and ask for data/code, but the core reinterpretation deserves to be aired. I'd cite the qualitative result cautiously, not the specific Tg.","headline":"Strong qualitative case for glassy arrest behind water's apparent LLT; the 189 K Tg is a plausible but under-supported extrapolation.","tokens_in":16654,"tokens_out":4154,"would_cite":true,"duration_ms":41972,"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":"What is often interpreted as water's liquid-liquid transition is actually the onset of glassy arrest, with a glass-transition temperature near 189 K.","keywords":["supercooled water","glass transition","liquid-liquid transition","kinetic arrest","Vogel-Fulcher-Tammann equation","dielectric relaxation","two-state fluctuations","machine-learning potentials"],"falsifier":"Measure structural relaxation times of pure bulk water at ambient pressure down to ~185 K using an ultrafast probe that outruns crystallization; if the relaxation time at 189 K remains many orders of magnitude below the 100-s mark—i.e., the liquid is still mobile—the inferred Tg and the glass-LLT identification would be refuted.","tokens_in":15684,"feed_emoji":"💧","tokens_out":7248,"duration_ms":69714,"temperature":0.7,"pith_summary":"This paper argues that the long-sought liquid-liquid transition in supercooled water is not a transition between two equilibrium liquids. The features commonly read as evidence of a second low-density liquid—energy fluctuations, two-state switching, a low-q scattering response—coincide with periods of near-zero molecular mobility, meaning the low-density state is a kinetically arrested glass. By extrapolating well-equilibrated dielectric relaxation times with the Vogel-Fulcher-Tammann equation, the authors obtain an ambient-pressure glass-transition temperature of 189 ± 8 K, squarely in the range where a liquid-liquid transition has been proposed. If correct, water's celebrated two-state phenomenology would be the boundary between an ergodic liquid and a non-ergodic glass, not a hidden critical point.","feed_headline":"Water's second liquid is actually a glass, simulations show","feed_subtitle":"The purported second liquid is likely a glass; Tg of water sits near 189±8 K, close to the proposed LLT.","key_machinery":"The load-bearing tool is the Vogel-Fulcher-Tammann fit to dielectric relaxation times (obtained from the total dipole autocorrelation), which extrapolates simulated dynamics—covering roughly 10^-9 to 10^-6 s—to the conventional glass-transition definition tau = 100 s. Supporting measures include time-resolved mean-square displacement to identify dynamical arrest, pressure/cooling ramps that generate rate-dependent glassy states, and isochronal lines (equal-relaxation-time contours) that map the non-ergodic boundary in the pressure-temperature plane. The temperature at which relaxation time varies most steeply with pressure, near 193.5 K, coincides with the onset of the purported two-state fl","core_discovery":"The central claim is that the purported low-density liquid (LDL) phase of supercooled water is dynamically arrested rather than liquid: in trajectories where LDL-HDL fluctuations were previously reported, the low-energy state shows essentially flat time-resolved mean-square displacement, indicating vanishing diffusion over hundreds of nanoseconds. Pressure- and temperature-ramp simulations show classic glassy signatures—rate-dependent final energies, physical aging, a step-like then peak-like heat capacity—and the onset of kinetic arrest coincides with the heat-capacity maximum. VFT extrapolation of dielectric relaxation times to tau = 100 s yields Tg = 189 ± 8 K at ambient pressure across t","pith_inferences":["The maximal pressure sensitivity of relaxation near 193.5 K implies that modest pressure noise in experiments or simulations can flip small water samples across the ergodic boundary; this may explain why some ultrafast droplet experiments report two-liquid-like features while others do not.","If the glassy-arrest scenario holds, the singularity-free interpretation of water's anomalies becomes a natural framework, since sharp response-function maxima can emerge from kinetic slowing without an underlying critical point.","The supplemental barostat-dependence result suggests a concrete test: systematically varying thermostat/barostat coupling and system size should make two-state fluctuations vanish under softer pressure control, providing a falsifiable simulation benchmark.","The VFT parameters could be used to predict Tg for heavy water or for water under confinement, generating testable dielectric and NMR predictions beyond the paper's direct claims."],"forward_implications":["The low-density 'second liquid' would not exist as an equilibrium phase; simulated HDL-LDL coexistence reflects intermittent trapping in a glassy state.","Free-energy calculations that sample the low-density state would be biased by incomplete equilibration, so computed LLT free-energy differences need revisiting.","Experimental signatures—low-q scattering, sharp heat-capacity peaks, Widom-line-like response maxima—appear where dynamics turn non-ergodic and do not require a hidden critical point.","Water's ambient-pressure glass transition would be near 189 K rather than 136 K, reconciling confinement, electron-diffraction, and minimally perturbing salt-solution measurements.","The mechanism may generalize: first-order-like coexistence between liquid and arrested states could masquerade as liquid-liquid transitions in other polyamorphic or deeply supercooled systems."],"fun_headline_variants":["Water's second 'liquid' is really a glass, simulations show","Supercooled water's LDL phase is a glass, not a liquid","Water's fabled second liquid? It's a glass near 189 K","Water's second liquid is a glassy state, simulations say","Glassy signatures in water's elusive second liquid"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The central number rests on assuming the Vogel-Fulcher-Tammann equation fitted to relaxation times spanning only about a microsecond continues to hold across the roughly eight orders of magnitude needed to reach the 100-second definition of the glass transition, with no change in relaxation mechanism.","fun_headline_variants_meta":{"raw":{"variants":["Water's second 'liquid' is really a glass, simulations show","Supercooled water's LDL phase is a glass, not a liquid","Water's fabled second liquid? It's a glass near 189 K","Water's second liquid is a glassy state, simulations say","Glassy signatures in water's elusive second liquid"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00016,"raw_usage":{"total_tokens":1056,"prompt_tokens":716,"completion_tokens":340,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":460,"completion_tokens_details":{"reasoning_tokens":251}},"tokens_in":460,"tokens_out":340,"duration_ms":4040,"temperature":1.0,"reasoning_tokens":251,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T16:57:41.340111+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure structural relaxation times of pure bulk water at ambient pressure down to ~185 K using an ultrafast probe that outruns crystallization; if the relaxation time at 189 K remains many orders of magnitude below the 100-s mark—i.e., the liquid is still mobile—the inferred Tg and the glass-LLT identification would be refuted.","supporting_citations":[],"review_version":1}