Freezing of 4He and its liquid-solid interface from density functional theory

We show that, at high densities, fully variational solutions of solidlike types can be obtained from a density functional formalism originally designed for liquid 4He . Motivated by this finding, we propose an extension of the method that accurately describes the solid phase and the freezing transition of liquid 4He at zero temperature. The density profile of the interface between liquid and the (0001) surface of the 4He crystal is also investigated, and its surface energy evaluated. The interfacial tension is found to be in semiquantitative agreement with experiments and with other microscopic calculations. This opens the possibility to use unbiased density functional (DF) methods to study highly nonhomogeneous systems, like 4He interacting with strongly attractive impurities and/or substrates, or the nucleation of the solid phase in the metastable liquid.

Metastable liquid-liquid phase transition in a single-component system with only one crystal phase and no density anomaly

We investigate the phase behavior of a single-component system in three dimensions with spherically-symmetric, pairwise-additive, soft-core interactions with an attractive well at a long distance, a repulsive soft-core shoulder at an intermediate distance, and a hard-core repulsion at a short distance, similar to potentials used to describe liquid systems such as colloids, protein solutions, or liquid metals. We showed [Nature (London) 409, 692 (2001)] that, even with no evidence of the density anomaly, the phase diagram has two first-order fluid-fluid phase transitions, one ending in a gas¿low-density-liquid (LDL) critical point, and the other in a gas¿high-density-liquid (HDL) critical point, with a LDL-HDL phase transition at low temperatures. Here we use integral equation calculations to explore the three-parameter space of the soft-core potential and perform molecular dynamics simulations in the interesting region of parameters. For the equilibrium phase diagram, we analyze the structure of the crystal phase and find that, within the considered range of densities, the structure is independent of the density. Then, we analyze in detail the fluid metastable phases and, by explicit thermodynamic calculation in the supercooled phase, we show the absence of the density anomaly. We suggest that this absence is related to the presence of only one stable crystal structure.

Thermal nucleation of cavities in liquid helium at negative pressures

We have investigated the nucleation rate at which cavities are formed in 4He and 3He at negative pressures due to thermal fluctuations. To this end, we have used a density functional that reproduces the He liquid-gas interface along the coexistence line. The inclusion of thermal effects in the calculation of the barrier against nucleation results in a sizable decrease of the absolute value of the tensile strength above 1.5 K.

Two-dimensional clusters of liquid 4He

The binding energies of two-dimensional clusters (puddles) of¿4He are calculated in the framework of the diffusion Monte Carlo method. The results are well fitted by a mass formula in powers of x=N-1/2, where N is the number of particles. The analysis of the mass formula allows for the extraction of the line tension, which turns out to be 0.121 K/Å. Sizes and density profiles of the puddles are also reported.

Cavitation in 3He-4He liquid mixtures at low temperatures

Bubble formation in solutions of 3He and 4He is studied within a density-functional approach. In particular, the temperature dependence of the cavitation pressure for different 3He concentrations is calculated at low temperatures and compared to that of pure 4He. The presence of Andreev states lowers the surface tension and, consequently, nucleation barriers are drastically reduced. This fact means that even at low 3He concentrations the cavitation process takes place at higher pressures than the spinodal pressure, which is not the case for pure 4He.

High-Momentum Response of Liquid 3He

A final-state-effects formalism suitable to analyze the high-momentum response of Fermi liquids is presented and used to study the dynamic structure function of liquid 3He. The theory, developed as a natural extension of the Gersch-Rodriguez formalism, incorporates the Fermi statistics explicitly through a new additive term which depends on the semidiagonal two-body density matrix. The use of a realistic momentum distribution, calculated using the diffusion Monte Carlo method, and the inclusion of this additive correction allows for good agreement with available deep-inelastic neutron scattering data.

Model-independent bounds for the potential and kinetic energy of liquid 4He at zero temperature

Using the experimental values of the chemical potentials of liquid 4He and of a 3He impurity in liquid 4He, we derive a model-independent lower (upper) bound to the kinetic (potential) energy per particle at zero temperature. The values of the bounds at the experimental saturation density are 13.42 K for the kinetic energy and -20.59 K for the potential energy. All the theoretical calculations based on the Lennard-Jones potential violate the upper-bound condition for the potential energy.

Intramolecular coupling as a mechanism for a liquid-liquid phase transition

We study a model for water with a tunable intramolecular interaction Js, using mean-field theory and off-lattice Monte Carlo simulations. For all Js>~0, the model displays a temperature of maximum density. For a finite intramolecular interaction Js>0, our calculations support the presence of a liquid-liquid phase transition with a possible liquid-liquid critical point for water, likely preempted by inevitable freezing. For J=0, the liquid-liquid critical point disappears at T=0.

Static and dynamic structure of liquid metals: Role of the different parts of the interaction potential

The influence of different parts of the interaction potential on the microscopic behavior of simple liquid metals is investigated by molecular dynamics simulation. The role of the soft-core repulsive, short-range attractive, and long-range oscillatory forces on the properties of liquid lithium close to the triple point is analyzed by comparing the results from simulations of identical systems but truncating the potential at different distances. Special attention is paid to dynamic collective properties such as the dynamic structure factors, transverse current correlation functions, and transport coefficients. It is observed that, in general, the effects of the short-range attractive forces are important. On the contrary, the influence of the oscillatory long-range interactions is considerably less, being the most pronounced for the dynamic structure factor at long wavelengths. The results of this work suggest that the influence of the attractive forces becomes less significant when temperature and density increase.

Finite-size scaling investigation of the liquid-liquid critical point in ST2 water and its stability with respect to crystallization.

The liquid-liquid critical point scenario of water hypothesizes the existence of two metastable liq- uid phases low-density liquid (LDL) and high-density liquid (HDL) deep within the supercooled region. The hypothesis originates from computer simulations of the ST2 water model, but the stabil- ity of the LDL phase with respect to the crystal is still being debated. We simulate supercooled ST2 water at constant pressure, constant temperature, and constant number of molecules N for N ≤ 729 and times up to 1 μs. We observe clear differences between the two liquids, both structural and dynamical. Using several methods, including finite-size scaling, we confirm the presence of a liquid-liquid phase transition ending in a critical point. We find that the LDL is stable with respect to the crystal in 98% of our runs (we perform 372 runs for LDL or LDL-like states), and in 100% of our runs for the two largest system sizes (N = 512 and 729, for which we perform 136 runs for LDL or LDL-like states). In all these runs, tiny crystallites grow and then melt within 1 μs. Only for N ≤ 343 we observe six events (over 236 runs for LDL or LDL-like states) of spontaneous crystal- lization after crystallites reach an estimated critical size of about 70 ± 10 molecules.

Urban population density functions : the case of the Barcelona region

L'anàlisi de la densitat urbana és utilitzada per examinar la distribució espacial de la població dins de les àrees urbanes, i és força útil per planificar els serveis públics. En aquest article, s'estudien setze formes funcionals clàssiques de la relació existent entre la densitat i la distancia en la regió metropolitana de Barcelona i els seus onze subcentres.

Cubic spline population density functions and subcentre delimitation. The case of Barcelona

The presence of subcentres cannot be captured by an exponential function. Cubic spline functions seem more appropriate to depict the polycentricity pattern of modern urban systems. Using data from Barcelona Metropolitan Region, two possible population subcentre delimitation procedures are discussed. One, taking an estimated derivative equal to zero, the other, a density gradient equal to zero. It is argued that, in using a cubic spline function, a delimitation strategy based on derivatives is more appropriate than one based on gradients because the estimated density can be negative in sections with very low densities and few observations, leading to sudden changes in estimated gradients. It is also argued that using as a criteria for subcentre delimitation a second derivative with value zero allow us to capture a more restricted subcentre area than using as a criteria a first derivative zero. This methodology can also be used for intermediate ring delimitation.