Visualization of supercritical water pseudo-boiling at Widom line crossover

Substantial breakdown of the hydrogen-bonding network, local density inhomogeneities and fluid-liquid structural transitions in supercritical octanol-1: A molecular dynamics investigation

Molecular dynamics simulations have been employed to explore the hydrogen-bonding structure and dynamics in supercritical octanol-1 at a near-critical temperature and up to high densities and pressures. A substantial breakdown of the hydrogen-bonding network when going from ambient-liquid to supercritical conditions is revealed. The fraction of the non-hydrogen bonded molecules significantly increases in supercritical octanol-1, and a substantial decrease in the intermittent hydrogen-bond lifetime is observed. This behavior is also reflected on the maximum local density augmentation, which is comparable to the values obtained for non-polar and non-hydrogen bonded fluids. The existence of a structural transition from an inhomogeneous fluid phase to a soft-liquid one at densities higher than 2.0 ρc is also revealed. At higher densities, a significant change in the reorientational relaxation process is observed, reflected on the significant increase in the ratio of the Legendre reorientational times τ1R/τ2R. The latter becomes much higher than the value predicted by the Debye model of diffusive reorientation and the corresponding ratio for ambient liquid octanol-1. The non-polar tail of octanol-1 under supercritical conditions reorients more slowly in comparison with the polar tail. Interestingly, the opposite behavior is observed for the ambient liquid, further verifying the strong effect of the breakdown of the hydrogen bonding network on the properties of supercritical octanol-1. In accordance with the above-mentioned findings, the static dielectric constant of supercritical octanol-1 is very low even at high densities and pressures, comparable to the values obtained for non-polar and non-hydrogen bonded fluids.

Critical Threshold for Bubble-like Nucleation during Pseudoboiling at Supercritical Pressures

Supercritical pseudoboiling was proposed in the 1950s-1960s. Recently, evaporation-like and boiling-like heat transfer have been directly observed in macroscopic scales, and the contribution of pseudoboiling to the total heat transfer rate has been quantitatively characterized experimentally. Here, we explore the critical threshold to generate a bubble-like nucleus at supercritical pressure at the atomic scale, characterized by the total energy (Te = Ke + Pe, where Ke and Pe are kinetic energy and potential energy, respectively). Molecular dynamics simulations are performed, including an argon fluid box heated by a solid wall having its temperature above the fluid temperature. The fluid pressure is controlled by a movable piston wall opposite the heating wall. The effects of pressure, nonuniform heating, and surface wettability on pseudoboiling are investigated. It is found that the criterion Te > 0 should be satisfied for subcritical boiling, matching that reported previously. The criterion for supercritical pseudoboiling was newly obtained such that Te > 0.012 eV at 8 MPa for argon, but the threshold increases as pressure increases. Nonuniform heating and surface wettability do not affect the critical threshold of Te for bubble-like nucleation but affect the location of the initially generated bubble-like nucleus and the stabilized pseudofilm or pseudonucleate heat transfer modes, where the former is similar to (vapor) film boiling and the latter is similar to nucleate boiling at subcritical pressure. Because pseudoboiling occurs without surface tension at supercritical pressure, we observe that the bubble-like structure may not display a perfectly smooth gas-liquid interface but may display an irregular pattern instead. Our work explains pseudoboiling from the viewpoint of the competition between kinetic energy and potential energy and presents a link regarding boiling in the two domains of subcritical pressure and supercritical pressure.

Crossover from gas-like to liquid-like molecular diffusion in a simple supercritical fluid

According to textbooks, no physical observable can be discerned allowing to distinguish a liquid from a gas beyond the critical point. Yet, several proposals have been put forward challenging this view and various transition boundaries between a gas-like and a liquid-like behaviour, including the so-called Widom and Frenkel lines, and percolation line, have been suggested to delineate the supercritical state space. Here we report observation of a crossover from gas-like (Gaussian) to liquid-like (Lorentzian) self-dynamic structure factor by incoherent quasi-elastic neutron scattering measurements on supercritical fluid methane as a function of pressure, along the 200 K isotherm. The molecular self-diffusion coefficient was derived from the best Gaussian (at low pressures) or Lorentzian (at high pressures) fits to the neutron spectra. The Gaussian-to-Lorentzian crossover is progressive and takes place at about the Widom line intercept (59 bar). At considerably higher pressures, a liquid-like jump diffusion mechanism properly describes the supercritical fluid on both sides of the Frenkel line. The present observation of a gas-like to liquid-like crossover in the self dynamics of a simple supercritical fluid confirms emerging views on the unexpectedly complex physics of the supercritical state, and could have planet-wide implications and possible industrial applications in green chemistry.

Who made the noise? Systematic approach for the assessment of neutron imaging scintillators

We propose a method to analyze the characteristics of scintillator screens for neutron imaging applications. Using calculations based on the theory of cascaded linear steps as well as experimental measurements, we compared the characteristics of different lithium- and gadolinium-based scintillator screens. Our results show that, despite their much lower light output, gadolinium-based scintillators outperform lithium-based scintillators in terms of noise characteristics for a variety of imaging setups. However, the relative performance of scintillator screens is highly dependent on the other setup characteristics such as the beam spectrum, field of view, used optical lens and size of the camera sensor. Consequently, the selection of the best scintillator screen - as well as the scintillator characteristics assessment in new developments - requires a systematic consideration of all these elements, as enabled by the framework presented here.

Supercritical fluids behave as complex networks

Supercritical fluids play a key role in environmental, geological, and celestial processes, and are of great importance to many scientific and engineering applications. They exhibit strong variations in thermodynamic response functions, which has been hypothesized to stem from the microstructural behavior. However, a direct connection between thermodynamic conditions and the microstructural behavior, as described by molecular clusters, remains an outstanding issue. By utilizing a first-principles-based criterion and self-similarity analysis, we identify energetically localized molecular clusters whose size distribution and connectivity exhibit self-similarity in the extended supercritical phase space. We show that the structural response of these clusters follows a complex network behavior whose dynamics arises from the energetics of isotropic molecular interactions. Furthermore, we demonstrate that a hidden variable network model can accurately describe the structural and dynamical response of supercritical fluids. These results highlight the need for constitutive models and provide a basis to relate the fluid microstructure to thermodynamic response functions.

Fingerprints of the Crossing of the Frenkel and Melting Line on the Properties of High-Pressure Supercritical Water

Using molecular dynamics simulations in combination with the two-phase thermodynamic model, we reveal novel characteristic fingerprints of the crossing of the Frenkel and melting line on the properties of high-pressure water at a near-critical temperature (1.03T_c). The crossing of the Frenkel line at about 1.17 GPa is characterized by a crossover in the rotational and translational entropy ratio S^rot/S^trans, indicating a change in the coupling between translational and rotational motions which is also reflected in the shape of the rotational density of states. The observed isosbestic points in the translational and rotational density of states are also blue-shifted at density and pressure conditions higher than the ones corresponding to the Frenkel line. The first-order phase transition from a rigid liquid to a face-centered cubic plastic crystal phase at about 8.5 GPa is reflected in the discontinuous changes in the translational and rotational entropy, particularly in the significant increase of the ratio S^rot/S^trans. A noticeable discontinuous increase of the dielectric constant has also been revealed when crossing this melting line, which is attributed to the different arrangement of the water molecules in the plastic crystal phase. The reorientational dynamics in the plastic crystal phase is faster in comparison with the "rigid" liquid-like phase, but it remains unchanged upon a further pressure increase in the range of 8.5-11 GPa.

Gas–liquid crossover in the Lennard-Jones system

It is demonstrated that the crossover between gas- and liquid-like regions on the phase diagram of the Lennard-Jones system occurs at a fixed value of the density divided by its value at the freezing point, ρ/ ρfr ≃ 0.35. This definition is consistent with other definitions proposed recently. As a result, a very simple practical expression for the gas-to-liquid crossover line emerges.

Numerical Investigation of Heat Transfer Characteristics of scCO2 Flowing in a Vertically-Upward Tube with High Mass Flux

In this work, the heat transfer characteristics of supercritical pressure CO₂ in vertical heating tube with 10 mm inner diameter under high mass flux were investigated by using an SST k-ω turbulent model. The influences of inlet temperature, heat flux, mass flux, buoyancy and flow acceleration on the heat transfer of supercritical pressure CO₂ were discussed. Our results show that the buoyancy and flow acceleration effect based on single phase fluid assumption fail to explain the current simulation results. Here, supercritical pseudo-boiling theory is introduced to deal with heat transfer of scCO₂. scCO₂ is treated to have a heterogeneous structure consisting of vapor-like fluid and liquid-like fluid. A physical model of scCO₂ heat transfer in vertical heating tube was established containing a gas-like layer near the wall and a liquid-like fluid layer. Detailed distribution of thermophysical properties and turbulence in radial direction show that scCO₂ heat transfer is greatly affected by the thickness of gas-like film, thermal properties of gas-like film and turbulent kinetic energy in the near-wall region. Buoyancy parameters Bu < 10⁻⁵, Bu* < 5.6 × 10⁻⁷ and flow acceleration parameter K_v < 3 × 10⁻⁶ in this paper, which indicate that buoyancy effect and flow acceleration effect has no influence on heat transfer of scCO₂ under high mass fluxes. This work successfully explains the heat transfer mechanism of supercritical fluid under high mass flux.

Deciphering the molecular mechanism of water boiling at heterogeneous interfaces

Water boiling control evolution of natural geothermal systems is widely exploited in industrial processes due to the unique non-linear thermophysical behavior. Even though the properties of water both in the liquid and gas state have been extensively studied experimentally and by numerical simulations, there is still a fundamental knowledge gap in understanding the mechanism of the heterogeneous nucleate boiling controlling evaporation and condensation. In this study, the molecular mechanism of bubble nucleation at the hydrophilic and hydrophobic solid-water interface was determined by performing unbiased molecular dynamics simulations using the transition path sampling scheme. Analyzing the liquid to vapor transition path, the initiation of small void cavities (vapor bubbles nuclei) and their subsequent merging mechanism, leading to successively growing vacuum domains (vapor phase), has been elucidated. The molecular mechanism and the boiling nucleation sites' location are strongly dependent on the solid surface hydrophobicity and hydrophilicity. Then simulations reveal the impact of the surface functionality on the adsorbed thin water molecules film structuring and the location of high probability nucleation sites. Our findings provide molecular-scale insights into the computational aided design of new novel materials for more efficient heat removal and rationalizing the damage mechanisms.

From Critical Point to Critical Point: The Two-States Model Describes Liquid Water Self-Diffusion from 623 to 126 K

Liquid’s behaviour, when close to critical points, is of extreme importance both for fundamental research and industrial applications. A detailed knowledge of the structural–dynamical correlations in their proximity is still today a target to reach. Liquid water anomalies are ascribed to the presence of a second liquid–liquid critical point, which seems to be located in the very deep supercooled regime, even below 200 K and at pressure around 2 kbar. In this work, the thermal behaviour of the self-diffusion coefficient for liquid water is analyzed, in terms of a two-states model, for the first time in a very wide thermal region (126 K < T < 623 K), including those of the two critical points. Further, the corresponding configurational entropy and isobaric-specific heat have been evaluated within the same interval. The two liquid states correspond to high and low-density water local structures that play a primary role on water dynamical behavior over 500 K.

On the Different Faces of the Supercritical Phase of Water at a Near-Critical Temperature: Pressure-Induced Structural Transitions Ranging from a Gaslike Fluid to a Plastic Crystal Polymorph

The present study reports a systematic analysis of a wide variety of structural, thermodynamic, and dynamic properties of supercritical water along the near-critical isotherm of T = 1.03T_c and up to extreme pressures, using molecular dynamics and Monte Carlo simulations. The methodology employed provides solid evidence about the existence of a structural transition from a liquidlike fluid to a compressed, tightly packed liquid, in the density and pressure region around 3.4ρ_c and 1.17 GPa, introducing an alternative approach to locate the crossing of the Frenkel line. Around 8.5 GPa another transition to a face-centered-cubic plastic crystal polymorph with density 5.178ρ_c is also observed, further confirmed by Gibbs free energy calculations using the two-phase thermodynamic model. The isobaric heat capacity maximum, closely related to the crossing of the Widom line, has also been observed around 0.8ρ_c, where the local density augmentation is also maximized. Another structural transition has been observed at 0.2ρ_c, related to the transformation of the fluid to a dilute gas at lower densities. These findings indicate that a near-critical isotherm can be divided into different domains where supercritical water exhibits distinct behavior, ranging from a gaslike one to a plastic crystal one.

Phase distribution including a bubblelike region in supercritical fluid

Pseudoboiling in supercritical fluid (SF) has been paid great attention in recent years. Available works mainly focus on thermodynamics analysis. Fewer studies were reported on the spatial time phase distribution. Here, SF is investigated in a multiphase fluid framework using molecular dynamics (MD) simulations. A simulation box contains 10 976 argon atoms, with periodic boundary conditions applied on all the box surfaces. Pressure and temperature are well controlled. Based on MD simulation results, an onset pseudoboiling temperature T^{-} and a termination pseudoboiling temperature T^{+} are defined using the neighboring molecules method, the radial distribution function method, and the two-body excess entropy method. The two transition temperatures divide the whole phase diagram into three regimes of liquidlike, two-phase-like (TPL), and gaslike, and the MD determined T^{-} and T^{+} well matched the thermodynamics-determined values. In the TPL regime, nanovoids are observed to have two distinct characteristics: (1) Particles are sparsely distributed to have gas density inside the void, but are densely populated to have liquid density outside the void. (2) Voids have a curved interface. These characteristics are very similar to bubble characteristics in subcritical pressure. Hence, voids in the supercritical state are called "bubblelike" in this paper. Nonlinear dynamics demonstrates chaotic behavior in the TPL regime, similar to the two-phase regime in the subcritical domain. The above findings give strong evidence that SF in the TPL regime consists of a mixture of bubblelike voids and surrounding liquids. Our work highlights the multiphase feature of a SF, hence, the well-established multiphase theory in subcritical pressures can be introduced to handle the complex SF.

Manipulation of epsilon-near-zero wavelength for the optimization of linear and nonlinear absorption by supercritical fluid

We introduce supercritical fluid (SCF) technology to epsilon-near-zero (ENZ) photonics for the first time and experimentally demonstrate the manipulation of the ENZ wavelength for the enhancement of linear and nonlinear optical absorption in ENZ indium tin oxide (ITO) nanolayer. Inspired by the SCF’s applications in repairing defects, reconnecting bonds, introducing dopants, and boosting the performance of microelectronic devices, here, this technique is used to exploit the influence of the electronic properties on optical characteristics. By reducing oxygen vacancies and electron scattering in the SCF oxidation process, the ENZ wavelength is shifted by 23.25 nm, the intrinsic loss is reduced by 20%, and the saturable absorption modulation depth is enhanced by > 30%. The proposed technique offers a time-saving low-temperature technique to optimize the linear and nonlinear absorption performance of plasmonics-based ENZ nanophotonic devices.

Quasi-equilibrium phase coexistence in single component supercritical fluids

In their supercritical state simple fluids are generally thought to assume a homogeneous phase throughout all combinations of pressures and temperatures, although various response functions or transport properties may exhibit anomalous behavior, characterizing a state point as either more gas-like or liquid-like, respectively. While a large body of results has been compiled in the last two decades regarding the details of the supercritical phase in thermodynamic equilibrium, far less studies have been dedicated to out-of-equilibrium situations that nevertheless occur along with the handling of substances such as carbon dioxide or Argon. Here we consider successive compression-expansion cycles of equal amounts of Argon injected into a high-pressure chamber, traversing the critical pressure at two times the critical temperature. Due to expansion cooling, the fluid temporarily becomes sub-critical, and light scattering experiments show the formation of sub-micron-sized droplets and nanometer-scale clusters, both of which are distinct from spontaneous density fluctuations of the supercritical background and persist for a surprisingly long time. A kinetic rate model of the exchange of liquid droplets with the smaller clusters can explain this behavior. Our results indicate non-equilibrium aspects of supercritical fluids that may prove important for their processing in industrial applications.

In their supercritical state simple fluids are generally thought to assume a homogeneous phase throughout. Lee et al. find that liquid droplets temporarily formed in a supercritical background after sub-critical injection can survive for a surprisingly long time.

Dynamic Crossover in Fluids: From Hard Spheres to Molecules

We propose a simple and generic definition of a demarcation reconciling structural and dynamic frameworks when combined with the entropy scaling framework. This crossover line between gas- and liquid-like behaviors is defined as the curve for which an individual property, the contribution to viscosity due to molecules' translation, is exactly equal to a collective property, the contribution to viscosity due to molecular interactions. Such a definition is shown to be consistent with the one based on the minima of the kinematic viscosity. For the hard sphere, this is shown to be an exact solution. For Lennard-Jones spheres and dimers and for some simple real fluids, this relation holds very well. This crossover line passes nearby the critical point, and for all studied fluids, it is well captured by the critical excess entropy curve for atomic fluids, emphasizing the link between transport properties and local structure.

A quantum magnetic analogue to the critical point of water

At the liquid-gas phase transition in water, the density has a discontinuity at atmospheric pressure; however, the line of these first-order transitions defined by increasing the applied pressure terminates at the critical point¹, a concept ubiquitous in statistical thermodynamics². In correlated quantum materials, it was predicted³ and then confirmed experimentally^4,5 that a critical point terminates the line of Mott metal-insulator transitions, which are also first-order with a discontinuous charge carrier density. In quantum spin systems, continuous quantum phase transitions⁶ have been controlled by pressure^7,8, applied magnetic field^9,10 and disorder¹¹, but discontinuous quantum phase transitions have received less attention. The geometrically frustrated quantum antiferromagnet SrCu₂(BO₃)₂ constitutes a near-exact realization of the paradigmatic Shastry-Sutherland model^12-14 and displays exotic phenomena including magnetization plateaus¹⁵, low-lying bound-state excitations¹⁶, anomalous thermodynamics¹⁷ and discontinuous quantum phase transitions^18,19. Here we control both the pressure and the magnetic field applied to SrCu₂(BO₃)₂ to provide evidence of critical-point physics in a pure spin system. We use high-precision specific-heat measurements to demonstrate that, as in water, the pressure-temperature phase diagram has a first-order transition line that separates phases with different local magnetic energy densities, and that terminates at an Ising critical point. We provide a quantitative explanation of our data using recently developed finite-temperature tensor-network methods^17,20-22. These results further our understanding of first-order quantum phase transitions in quantum magnetism, with potential applications in materials where anisotropic spin interactions produce the topological properties^23,24 that are useful for spintronic applications.

Thermodynamics and Dynamics of Supercritical Water Pseudo-Boiling

Supercritical fluid pseudo-boiling (PB), recently brought to the attention of the scientific community, is the phenomenon occurring when fluid changes its structure from liquid-like (LL) to gas-like (GL) states across the Widom line. This work provides the first quantitative analysis on the thermodynamics and the dynamics of water's PB, since the understanding of this phase transition is mandatory for the successful implementation of technologies using supercritical water (scH2O) for environmental, energy, and nanomaterial applications. The study combines computational techniques with in situ neutron imaging measurements. The results demonstrate that, during isobaric heating close to the critical point, while water density drops by a factor of three in the PB transitional region, the system needs >16 times less energy to increase its temperature by 1 K than to change its structure from LL to GL phase. Above the PB-Widom line, the structure of LL water consists mainly of tetramers and trimers, while below the line mostly dimers and monomers form in the GL phase. At atomic level, the PB dynamics are similar to those of the subcritical water vaporization. This fundamental knowledge has great impact on water science, as it helps to establish the structure-properties relationship of scH2O.

Universality, Scaling, and Collapse in Supercritical Fluids

Supercritical fluid (SCF) is known to exhibit salient dynamic and thermodynamic crossovers and an inhomogeneous molecular distribution. However, the question as to what basic physics underlies these microscopic and macroscopic anomalies remains open. Here, using an order parameter extracted by machine learning, the fraction of gas-like (or liquid-like) molecules, we find simplicity and universality in SCF: First, all isotherms of a given fluid collapse onto a single master curve described by a scaling relation. The observed power law holds from the high-temperature and -pressure regime down to the critical point where it diverges. Second, phase diagrams of different compounds collapse onto their master curves by the same scaling exponent, thereby demonstrating a putative law of corresponding supercritical states in simple fluids. The reported results support a model of the SCF as a mixture of two interchangeable microstates, whose spatiotemporal dynamics gives rise to unique macroscopic properties.

The use of process simulation in supercritical fluids applications

Modelling and simulation from micro- to macro-scale are needed to attain a broader commercialization of supercritical technologies.

Supercritical water anomalies in the vicinity of the Widom line

Supercritical water is used in a variety of chemical and industrial applications. As a consequence, a detailed knowledge of the structure-properties correlations is of uttermost importance. Although supercritical water was considered as a homogeneous fluid, recent studies revealed an anomalous behaviour due to nanoscale density fluctuations (inhomogeneity). The inhomogeneity is clearly demarked through the Widom line (maxima in response factions) and drastically affect the properties. In the current study the physical properties of supercritical water have been determined by classical molecular dynamics simulations using a variety of polarized and polarizable interatomic potentials. Their validity which was not available at supercritical conditions has been assessed based on the ability to reproduce experimental data. Overall, the polarized TIP4P/2005 model accurately predicted the properties of water in both liquid-like and gas-like regions. All interatomic potentials captured the anomalous behaviour providing a direct evidence of molecular-scale inhomogeneity.