Widom Line in Supercritical Water in Terms of Changes in Local Structure: Theoretical Perspective

Performing molecular dynamics simulations with the TIP4P/2005 water model along 9 isobars (from 175 to 375 bar) in the temperature range between 300 and 1100 K, we have found that the loci of the extrema in the rate of change of specific structural properties can be used to define purely structure-based Widom lines. We have examined several parameters that describe the local structure of water, such as the tetrahedral arrangement, nearest neighbor distance, local density around the molecules, and the size of the largest dense domain. The last two parameters were determined using the Voronoi polyhedral and density-based spatial clustering of applications with noise methods, respectively. By analyzing the moments of the associated distributions, we show that along a given isobar, the temperature at which we observe a maximum in the fluctuation, the rate of change of the average values, or in the skewness values unambiguously determines the Widom line that is in agreement with the experimentally detected, thermodynamic response function-based ones.