Nuclear Quantum Effects on the Thermodynamic Response Functions of a Polymorphic Waterlike Monatomic Liquid
Yizhi Liu,Gang Sun,Ali Eltareb,Gustavo E. Lopez,Nicolas Giovambattista,Limei Xu
DOI: https://doi.org/10.1103/physrevresearch.2.013153
2020-01-01
Physical Review Research
Abstract:Water and hydrogen are examples of substances proposed to exhibit a liquid-liquid critical point (LLCP) at conditions where nuclear quantum effects are relevant. The LLCP is usually accompanied by lines of maxima in density rho and thermodynamic response functions, such as isothermal compressibility kappa(T) and isobaric heat capacity C-P, in the supercritical region of the P-T plane. In the case of water, the rho- and kappa(T)-maxima lines can be accessed in experiments, while, instead, the LLCP has not been observed due to rapid crystallization. In this work, we study the nuclear quantum effects on a monatomic liquid that exhibits waterlike anomalous properties and a LLCP. By performing path-integral Monte Carlo simulations with different values of the Planck's constant h, we are able to explore how the location of the LLCP in the P-T plane and, in particular, the maxima lines in the supercritical region, shift as the system evolves from classical, h = 0, to quantum, h > 0. We find that as the quantum nature of the liquid (as quantified by h) increases and the atoms in the liquid become more delocalized, the LLCP shifts towards higher pressures and lower temperatures while the LLCP volume remains constant. Similar shifts (towards higher pressures and lower temperatures) are found in the case of the C-P- and kappa(T)-maxima lines. Instead, the.-maxima line extends towards higher temperatures and expands over a wider pressure interval as the liquid becomes more quantum. It follows that the nuclear quantum effects on the location of the LLCP may be estimated from the shift in C-P- and kappa(T)-maxima lines but not on measurements of the rho-maxima line. Interestingly, nuclear quantum effects considerably alter the slope of the liquid-liquid coexistence line and C-P-maxima line in the P-T plane while the slope of the kappa(T)-maxima line along the Widom line is barely affected. We discuss briefly the implications of our results to the case of H2O/D2O.