{"id":"704a5f57-ddeb-4348-8f1f-91804a07fd73","arxiv_id":"2412.12451","paper_version":3,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Quantum anharmonic effects, computed with neural canonical transformations, lower the bcc-fcc transition temperature in lithium, and hybrid functional corrections stabilize the high-pressure oC88 phase.","lead":"A neural-network quantum method was applied to lithium crystals, showing that nuclear quantum effects lower the predicted bcc-fcc transition temperature and that the high-pressure oC88 phase is stabilized by electronic structure accuracy rather than thermal motion. The work demonstrates a practical route to include quantum anharmonicity in first-principles phase diagrams.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"oC88 stabilization claim rests on a single-point HSE06 correction at 70 GPa applied as a constant shift; HSE relaxation and pressure transferability are untested.","rationale":"The reader identified essentially the same weak point: the oC88 stabilization conclusion depends on single-point HSE corrections without re-relaxation or validation against a higher-level method, applied across a pressure range from one pressure point. I agree that this is the most load-bearing concern. The NCT machinery itself appears credible: the method is benchmarked on a strongly anharmonic double-well potential, the code is open source, and the ABACUS and FHI-aims HSE results agree with each other. However, those two HSE calculations share the same functional and the same PBE/NCT geometries, so they do not test pressure transferability or whether structures shift under HSE. The HSE correction is large enough (about 6–8 meV/atom) to reverse the PBE ordering, but it is applied as a rigid shift, and a small pressure dependence or structural relaxation effect could change or eliminate the predicted oC88 stability window. For the bcc-fcc transition, the unpropagated DP-model systematic errors are also a real issue, but I do not elevate that to the primary concern because both NCT and classical MD use the same DP surface, so the qualitative direction of the quantum correction is less sensitive to the model error; the absolute transition temperatures remain less certain. Thus the reader's conditional verdict is appropriate, and no change to that verdict is needed.","tokens_in":18243,"tokens_out":5631,"duration_ms":55536,"concrete_test":"Run a full HSE06 structural relaxation of both cI16 and oC88 at 70 GPa, relaxing supercell volume and internal coordinates starting from the NCT/PBE geometries, and recompute the enthalpy difference; then repeat the HSE06 single-point energy difference at 60 and 80 GPa on the HSE-relaxed geometries. If the relaxed oC88−cI16 HSE enthalpy gap is not negative, or if the HSE correction changes by more than 1 meV/atom between 60 and 70 GPa, the predicted 62 GPa phase transition is not robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's most consequential claim is that oC88 is stabilized not by nuclear quantum effects but by the electronic-structure functional. The evidence is SM Table S2: single-point HSE06 on NCT/PBE geometries at 70 GPa changes the oC88−cI16 potential-energy difference by −6.17 meV/atom (ABACUS) or −7.67 meV/atom (FHI-aims), and this shift is drawn in Fig. 3(c) as a constant line over 60–90 GPa. Load-bearing assumptions are: (1) the HSE correction is pressure-independent or only weakly so; (2) PBE-optimized volumes and internal coordinates remain near HSE equilibrium for both the 88- and 432-atom cells; and (3) the PBE-based NCT free-energy contributions (anharmonic ZPE, entropy) are unchanged on the HSE surface. The correction is comparable to the nuclear quantum effects it dominates (about 4.9 meV ZPE difference), so a 1–2 meV/atom error or drift across 60–70 GPa could shift the predicted boundary by several GPa or remove it. Since HSE06 is not benchmarked against a higher-level electronic-structure method for this poor metal, and no re-relaxation is performed, the central high-pressure conclusion is not yet settled.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents neural canonical transformations (NCT), a variational density-matrix method that combines normalizing-flow phonon wave functions with a product ansatz for phonon occupation probabilities, and applies this method to solid lithium. Using a deep-potential surrogate for the PBE Born-Oppenheimer energy surface, the authors compute anharmonic free energies for bcc/fcc at 0-2 GPa and for cI16/oC88 at high pressure. They report that quantum anharmonicity lowers the bcc-fcc transition temperature relative to classical molecular dynamics (84, 142, and 196 K versus 144, 185, and 218 K at 0, 1, and 2 GPa), that the predicted cI16 fractional coordinates agree with experiment, and that a single-point HSE06 electronic-structure correction makes oC88 more stable than cI16 below roughly 62 GPa at 100 K, leading them to conclude that the experimentally observed oC88 phase is stabilized by electronic-structure effects rather than by nuclear quantum or thermal effects.","tokens_in":18468,"tokens_out":8102,"duration_ms":72322,"significance":"If fully supported, these results would be significant for both methodology and lithium phase diagram physics. The paper introduces an open-source, variational framework that goes beyond Gaussian phonon ansätze, and the low-pressure conclusion that quantum anharmonicity lowers the bcc-fcc transition temperature is a concrete, falsifiable prediction. The high-pressure claim—that the poor-metal oC88 structure is stabilized by the electronic-structure functional—offers a new resolution to a known discrepancy. The strengths of the manuscript include the principled variational formalism, the reproducible open-source code, the use of two independent electronic-structure codes for the HSE correction, and the explicit reporting of statistical errors. However, the high-pressure conclusion rests on untested assumptions about the pressure transferability and structural transferability of a single-point HSE correction, and the bcc-fcc transition temperatures are not error-propagated with respect to the quoted deep-potential systematic errors.","major_comments":[{"comment":"The central high-pressure conclusion is based on a single-point HSE06 correction evaluated at 70 GPa and applied in Fig. 3(c) as a pressure-independent shift over the 60-90 GPa range. The predicted cI16-oC88 transition occurs near 62 GPa, 8 GPa away from the pressure at which the correction is computed, and no HSE calculation at other pressures or HSE-level volume relaxation is reported. Because the correction (6.17-7.67 meV/atom) is larger than the PBE-level free-energy difference it is correcting (approximately 3.96 meV/atom at 70 GPa), a modest pressure dependence of the HSE shift could move or eliminate the predicted transition. The authors should assess pressure transferability, for example by computing HSE at 62 and 80 GPa or by re-relaxing both structures under HSE.","section":"Fig. 3(c) and SM §SII.D, Table S2"},{"comment":"The HSE stabilization is computed as a potential-energy difference at fixed PBE/NCT-optimized structures, while the NCT free-energy contributions (zero-point energy, anharmonic entropy, thermal terms) are taken unchanged from the PBE surface. This assumes that HSE and PBE have essentially the same equilibrium volumes, internal coordinates, and phonon curvatures for both the 432-atom cI16 and 352-atom oC88 cells. The claim that oC88 is stabilized by the potential energy surface 'rather than thermal or quantum nuclear effects' is therefore only as strong as this assumption. A concrete test would be to re-relax both structures with HSE, or at least to recompute the electronic energy along NCT-sampled configurations, and to verify that the relative free energy remains negative when the HSE surface is used.","section":"SM §SII.D, Table S2; main-text section on high-pressure structural stability"},{"comment":"The bcc-fcc transition temperatures are extracted from Gibbs free-energy differences of order 0.1-0.6 meV/atom, whereas the quoted systematic error of the low-pressure deep-potential model is 0.2-0.5 meV/atom. Since the manuscript itself notes that an error of 1 meV can shift the transition temperature by more than 100 K, the reported transition temperatures of 84, 142, and 196 K carry unquantified systematic uncertainty. The comparison with classical MD on the same energy surface supports the direction of the quantum-anharmonicity correction, but the precise transition temperatures, and the claimed agreement with experiment (77-88 K at 0 GPa), need either an error analysis with multiple DP models or direct DFT validation at selected thermodynamic points.","section":"SM §SII.B and Fig. 2"},{"comment":"The claim that the poor-metal oC88 phase is stabilized by HSE is not benchmarked against a higher-level electronic-structure method. For a system whose relative stability changes by about 6-8 meV/atom when switching from PBE to HSE, the possibility of overcorrection is a genuine correctness risk. Before asserting the mechanism as established, the authors should validate the HSE energy ordering with an independent higher-level method, such as quantum Monte Carlo or RPA, or at least with a different hybrid-functional parameterization.","section":"Main text 'High-pressure structural stability' and SM §SII.D"}],"minor_comments":[{"comment":"The text states that the anharmonic zero-point energy of oC88 decreases by 1.67 meV/atom, but the harmonic and anharmonic values in the figure (89.43 and 88.10 meV/atom) give a decrease of 1.33 meV/atom; this inconsistency should be corrected.","section":"Fig. 3(b) and main text"},{"comment":"The caption and legend contain typographical errors: 'Neural cannonical transformation' should read 'Neural canonical transformation', and the space-group symbol 'c2mb' in the oC88 panel should be 'C2mb'.","section":"Fig. 3(a)"},{"comment":"The caption reads 'Ironic entropy' and should read 'Ionic entropy'.","section":"SM Fig. S4(b)"},{"comment":"Several axis labels and legends contain garbled or missing symbols, such as 'opt  only' and 'T emperature'; these should be cleaned before publication.","section":"Figs. 1 and 2"},{"comment":"The energy-level cutoff K=20 for the product spectrum ansatz is stated, but no convergence test with respect to K is reported; a brief convergence statement would strengthen the variational free-energy values.","section":"SM §SII.B"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of the journal and the NCT methodology is interesting. The most novel claim—that oC88 is stabilized by the electronic-structure functional—is not yet settled because it relies on a single-point HSE correction at one pressure applied as a constant shift. I recommend requesting the additional HSE calculations and error-propagation analysis described in the major comments before accepting the paper."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"What you should know: this is a careful application of the neural canonical transformation method (prior work from the same group) to lithium, and it produces two genuinely new results. First, quantum anharmonicity lowers the bcc-fcc transition temperature compared to classical MD on the same deep-potential surface, bringing predictions closer to experiment. Second, the paper argues that the high-pressure oC88 phase is stabilized by electronic-structure accuracy (HSE06) rather than thermal or nuclear quantum effects—a direct reversal of earlier attributions. The cI16 fractional coordinates also match experiment, which is a nice validation of the technique.\n\nThe paper does several things well. The variational framework is principled, statistical errors are tiny, and they check their anharmonic machinery against a double-well benchmark. The HSE correction is calculated with two independent electronic-structure codes (ABACUS and FHI-aims) and the code is open-sourced. The comparison to classical MD on the same potential is the right way to isolate quantum effects.\n\nThe soft spots are real but not disqualifying. The high-pressure conclusion rests on a single-point HSE06 correction at 70 GPa, applied as a constant shift over 60–90 GPa. There is no re-relaxation under HSE and no test of pressure transferability. The correction is 6–8 meV/atom, which is the same order as the nuclear quantum effects it overrides, so a modest error could shift the boundary by several GPa. The stress-test note is right to flag this. The bcc-fcc transition temperatures are also more uncertain than the paper suggests: the free-energy differences are sub-meV while the deep-potential systematic error is quoted as 0.2–0.5 meV/atom. They do not propagate this error, though their MD comparison on the same potential does partially address it. The paper is transparent about these limitations in the SM, which helps.\n\nWho is this for? People working on lithium phase diagrams, anharmonic quantum solids, or nuclear quantum effect methods. The bcc-fcc result is a solid data point; the oC88 claim is a well-posed hypothesis that deserves further testing, ideally with HSE structure relaxation or a higher-level electronic-structure benchmark.\n\nRecommendation: this deserves a serious referee. I would send it to peer review and ask the authors to address the pressure dependence of the HSE correction and quantify the DP-model error more explicitly.","headline":"Solid application of a known method with a plausible but not yet airtight claim that oC88 stability comes from electronic structure, not nuclear quantum effects.","tokens_in":19010,"tokens_out":2126,"would_cite":true,"duration_ms":20044,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["64.70.kd","63.20.-e","71.15.Mb"],"model":"deepseek-v4-flash","headline":"Neural canonical transformations reveal that quantum anharmonicity lowers lithium's bcc-fcc transition temperature, and that the high-pressure oC88 phase is stabilized by the electronic potential energy surface rather than by nuclear…","keywords":["neural canonical transformation","quantum anharmonicity","lithium phase diagram","bcc-fcc transition","oC88 structure","hybrid functional correction","normalizing flow","variational free energy"],"falsifier":"Perform full structure relaxations of cI16 and oC88 with the HSE06 hybrid functional at 60, 70, and 80 GPa, and check whether oC88's enthalpy falls below cI16 around 62 GPa; if it does not, the paper's stabilization mechanism fails.","tokens_in":18026,"feed_emoji":"⚛️","tokens_out":8913,"duration_ms":73570,"temperature":0.7,"pith_summary":"This paper asks how nuclear quantum motion and anharmonicity affect the phase stability of solid lithium, and whether the puzzling high-pressure oC88 phase can be explained at all. The authors use neural canonical transformations, a variational density-matrix method that combines normalizing-flow phonon wave functions with learned occupation probabilities, to compute free energies of bcc, fcc, cI16, and oC88 lithium. They find that quantum anharmonicity lowers the bcc-fcc transition temperature to 84, 142, and 196 K at 0, 1, and 2 GPa, closer to experiment than classical molecular dynamics. They also find that anharmonic effects actually widen the free-energy gap between oC88 and cI16, and that it is only after applying an HSE06 hybrid-functional energy correction that oC88 becomes stable, predicting a cI16-oC88 transition near 62 GPa at 100 K. The message is that for lithium, the electronic potential energy surface, not nuclear quantum or thermal effects, is what stabilizes oC88.","feed_headline":"Electrons, not vibrations, stabilize lithium's oC88 phase","feed_subtitle":"Hybrid-functional corrections put the cI16-oC88 transition near 62 GPa at 100 K, matching experiments.","key_machinery":"The load-bearing object is the neural canonical transformation ansatz for the variational density matrix, $\\rho = \\sum_n p_n |\\Psi_n\\rangle\\langle\\Psi_n|$, with occupation probabilities $p_n$ modeled by a product spectrum ansatz and phonon wave functions $|\\Psi_n\\rangle = U_\\theta |\\Phi_n\\rangle$ built by a normalizing flow that maps phonon coordinates $q$ to quasi-phonon coordinates $\\xi = f_\\theta(q)$ and enters the wave function through a Jacobian determinant. This gives orthogonal, non-Gaussian wave functions for excited vibrational states of hundreds-atom supercells, so the anharmonic potential energy surface is treated exactly rather than through a Taylor expansion or a Gaussian variational density matrix. The free energy $F = \\mathbb{E}[k_B T \\ln p_n + \\text{local energy}]$ is minimized jointly over the probability and wave-function parameters, and constant-pressure calculations are closed by computing the stress tensor; the anharmonic phonon frequencies and zero-point energies appearing in the phase-stability comparison are read off from single-phonon excitation energies of the optimized states. The Born-Oppenheimer surface is represented by a machine-learned potential trained on PBE density functional theory, which supplies energies and forces for the large supercells.","core_discovery":"On the paper's own terms, the central discovery is that quantum anharmonicity and electronic-structure accuracy play opposite roles in lithium's high-pressure phase diagram. Neural canonical transformations show that anharmonicity lowers the zero-point energy of cI16 by about 4.92 meV/atom but oC88 by only about 1.67 meV/atom, so nuclear quantum effects actually worsen oC88's stability relative to cI16. A single-point HSE06 calculation on the NCT-optimized structures reverses the picture: it lowers oC88 relative to cI16 by 6.17 meV/atom compared with PBE, enough to place the cI16-oC88 transition at roughly 62 GPa and 100 K, in line with experiments that put oC88 between about 62 and 70 GPa. In the same framework, the bcc-fcc transition temperature comes out at 84, 142, and 196 K at 0, 1, and 2 GPa, consistently below classical MD values of 144, 185, and 218 K and closer to the experimental low-pressure boundary. The paper also reports that NCT-optimized fractional coordinates of the cI16 Wyckoff position reproduce the experimental values reported for the 16c site.","pith_inferences":["Because the HSE06 correction in the paper is a single-point calculation at 70 GPa on geometries relaxed with a PBE-trained potential, repeating the correction at 60 and 80 GPa on re-relaxed structures would directly test whether the predicted 62 GPa transition is an artifact of the single-point approximation.","The same variational machinery could be applied to other poor metals where PBE over-stabilizes metallic states; hybrid-functional corrections might expose narrow stability windows that classical and PBE-based structure searches miss.","The reported bcc-fcc transition temperatures carry an unquantified systematic error from the machine-learned potential (0.2-0.5 meV/atom), so a higher-accuracy potential could move the 84 K boundary by tens of kelvin; this is worth testing with path-integral calculations on a denser DFT dataset.","If anharmonicity consistently softens the more symmetric phase (cI16 here) while the electronic correction favors the less symmetric one (oC88), then the two effects oppose each other; in other light-element solids the same competition could produce phase sequences that depend sensitively on both functional accuracy and nuclear mass."],"forward_implications":["The bcc-fcc phase boundary in lithium is set roughly 40-60 K lower than classical molecular dynamics predicts, so quantum anharmonicity must be included when comparing computed with measured transition temperatures.","oC88's stability window between cI16 and oC40 is controlled by the electronic potential energy surface, meaning future crystal-structure searches for lithium should not rely on PBE alone for poor metallic phases.","Anharmonicity actually increases the free-energy difference between oC88 and cI16, so nuclear quantum effects are not the missing ingredient that earlier studies invoked.","Neural canonical transformations deliver anharmonic phonon spectra and zero-point energies for supercells of several hundred atoms, making such calculations feasible for other quantum solids like hydrogen, helium, and hydrides.","The cI16-oC88 transition near 62 GPa and 100 K predicted after the HSE06 correction matches the experimentally observed narrow stability range of oC88."],"supporting_citations":[{"why":"Experimental observation of the oC88 phase between 62 and 70 GPa; the transition pressure the paper aims to reproduce.","marker":"[17]"},{"why":"Experimental and DFT fcc-bcc phase boundary; the low-pressure transition line the NCT results are compared against.","marker":"[22]"},{"why":"Earlier zero-temperature calculation concluding oC88 is second-most stable and blaming ZPE/thermal effects; the hypothesis this paper overturns.","marker":"[27]"},{"why":"Earlier PBE-based claim that oC88 becomes stable with harmonic ZPE at 200 K; the result the paper shows is not reproduced by anharmonic NCT.","marker":"[28]"},{"why":"Source of the machine-learned potential and the recent conclusion that anharmonicity does not stabilize oC88; provides the base model and the competing result.","marker":"[29]"},{"why":"The specific neural canonical transformation variational density-matrix formulation that the paper extends to lithium.","marker":"[39]"},{"why":"PBE functional whose energies define the base Born-Oppenheimer surface used in the machine-learned potential.","marker":"[49]"},{"why":"Experimental fractional coordinates for cI16 that the NCT-optimized structure reproduces.","marker":"[63]"},{"why":"Defines the HSE06 hybrid functional used for the single-point electronic-structure correction.","marker":"[80]"}],"fun_headline_variants":["Electrons, not vibrations, settle lithium's oC88 phase","Hybrid functional, not nuclear motion, stabilizes lithium's oC88","Electronic accuracy dethrones quantum motion for lithium's oC88","Neural canonical transforms: electrons, not zero-point motion, pick oC88","Hybrid functional flips lithium's oC88 stability away from vibrations"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The high-pressure conclusion assumes that a single HSE06 energy correction, computed once at 70 GPa on structures relaxed with a PBE-trained machine-learned potential, gives the correct relative stability of cI16 and oC88 across the whole 60-80 GPa window without re-relaxing the structures under HSE.","fun_headline_variants_meta":{"raw":{"variants":["Electrons, not vibrations, settle lithium's oC88 phase","Hybrid functional, not nuclear motion, stabilizes lithium's oC88","Electronic accuracy dethrones quantum motion for lithium's oC88","Neural canonical transforms: electrons, not zero-point motion, pick oC88","Hybrid functional flips lithium's oC88 stability away from vibrations"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00086,"raw_usage":{"total_tokens":3759,"prompt_tokens":999,"completion_tokens":2760,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":615,"completion_tokens_details":{"reasoning_tokens":2673}},"tokens_in":615,"tokens_out":2760,"duration_ms":17810,"temperature":1.0,"reasoning_tokens":2673,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T14:03:26.276157+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Perform full structure relaxations of cI16 and oC88 with the HSE06 hybrid functional at 60, 70, and 80 GPa, and check whether oC88's enthalpy falls below cI16 around 62 GPa; if it does not, the paper's stabilization mechanism fails.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Experimental observation of the oC88 phase between 62 and 70 GPa; the transition pressure the paper aims to reproduce."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Experimental and DFT fcc-bcc phase boundary; the low-pressure transition line the NCT results are compared against."},{"cited_title":"Marqu ´es, M","cited_arxiv_id":null,"evidence_quote":"Earlier zero-temperature calculation concluding oC88 is second-most stable and blaming ZPE/thermal effects; the hypothesis this paper overturns."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Earlier PBE-based claim that oC88 becomes stable with harmonic ZPE at 200 K; the result the paper shows is not reproduced by anharmonic NCT."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Source of the machine-learned potential and the recent conclusion that anharmonicity does not stabilize oC88; provides the base model and the competing result."},{"cited_title":"Zhang, R.-S","cited_arxiv_id":null,"evidence_quote":"The specific neural canonical transformation variational density-matrix formulation that the paper extends to lithium."},{"cited_title":"Hanfland, K","cited_arxiv_id":null,"evidence_quote":"Experimental fractional coordinates for cI16 that the NCT-optimized structure reproduces."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the HSE06 hybrid functional used for the single-point electronic-structure correction."}],"review_version":1}