Path integral Monte Carlo study of 4He clusters doped with alkali and alkali-earth ions.

Path integral Monte Carlo calculations of (4)He nanodroplets doped with alkali (Na(+), K(+) and Cs(+)) and alkali-earth (Be(+) and Mg(+)) ions are presented. We study the system at T = 1 K and between 14 and 128 (4)He atoms. For all studied systems, we find that the ion is well localized at the center of the droplet with the formation of a "snowball" of well-defined shells of localized (4)He atoms forming solid-like order in at least the first surrounding shell. The number of surrounding helium shells (two or three) and the number of atoms per shell and the degree of localization of the helium atoms are sensitive to the type of ion. The number of (4)He atoms in the first shell varies from 12 for Na(+) to 18 for Mg(+) and depends weakly on the size of the droplet. The study of the density profile and of the angular correlations shows that the local solid-like order is more pronounced for the alkali ions with Na(+) giving a very stable icosahedral order extending up to three shells.

Liquid-solid transition in fully ionized hydrogen at ultra-high pressures.

We study the phase diagram of an effective ion model of fully ionized hydrogen at ultra-high pressure. We assume that the protons interact with a screened Coulomb potential derived from a static linear response theory. This model accurately reproduces the physical properties of hydrogen for densities greater than g/ρ(m)=10 cm(3) corresponding to the range of the coupling parameter r(s) ≲ 0.6. The pressure range, P ≳ 20 TPa, is well beyond present experimental limitations. Assuming classical protons, we find that the zero temperature enthalpy of the perfect bcc crystal is slightly lower than for other structures at g/ρ(m)=12.47 cm(3) while the fcc structure gains stability at higher density. Using Monte Carlo calculations, we compute the free energy of various phases and locate the melting transition versus density. We find that on melting, bcc is energetically favored with respect to fcc over the entire range investigated. In the solid phase the system undergoes a structural transition from bcc at higher temperature to fcc at lower temperature. The free energy difference between these two structures is very small so that obtaining a quantitative estimate of this second transition line requires accuracy beyond that provided by our method. We estimate the effect of proton zero point motion on the bcc melting line for hydrogen, deuterium, and tritium by a path integral Monte Carlo method. Although zero point effects on hydrogen are large, since the two competing phases (bcc and liquid) have locally similar environments, the effect on the melting line is small; the melting temperature for hydrogen is lowered by about 10% with respect to the classical value.

Evidence for a first-order liquid-liquid transition in high-pressure hydrogen from ab initio simulations.

Using quantum simulation techniques based on either density functional theory or quantum Monte Carlo, we find clear evidence of a first-order transition in liquid hydrogen, between a low conductivity molecular state and a high conductivity atomic state. Using the temperature dependence of the discontinuity in the electronic conductivity, we estimate the critical point of the transition at temperatures near 2,000 K and pressures near 120 GPa. Furthermore, we have determined the melting curve of molecular hydrogen up to pressures of 200 GPa, finding a reentrant melting line. The melting line crosses the metalization line at 700 K and 220 GPa using density functional energetics and at 550 K and 290 GPa using quantum Monte Carlo energetics.

Equation of state of metallic hydrogen from coupled electron-ion Monte Carlo simulations.

We present a study of hydrogen at pressures higher than molecular dissociation using the coupled electron-ion Monte Carlo method. These calculations use the accurate reptation quantum Monte Carlo method to estimate the electronic energy and pressure while doing a Monte Carlo simulation of the protons. In addition to presenting simulation results for the equation of state over a large region of the phase diagram, we report the free energy obtained by thermodynamic integration. We find very good agreement with density-functional theory based molecular-dynamics calculations for pressures beyond 600 GPa and densities above rho=1.4 g/cm(3) , both for thermodynamic and structural properties. This agreement provides a strong support to the different approximations employed in the density-functional treatment of the system, specifically the approximate exchange-correlation potential and the use of pseudopotentials for the range of densities considered. We find disagreement with chemical models, which suggests that a reinvestigation of planetary models--previously constructed using the Saumon-Chabrier-Van Horn equations of state--might be needed.

Hexatic and mesoscopic phases in a 2D quantum coulomb system.

We study the Wigner crystal melting in a two-dimensional quantum system of distinguishable particles interacting via the 1/r Coulomb potential. We use quantum Monte Carlo methods to calculate its phase diagram, locate the Wigner crystal region, and analyze its instabilities towards the liquid phase. We discuss the role of quantum effects in the critical behavior of the system, and compare our numerical results with the classical theory of melting, and the microemulsion theory of frustrated Coulomb systems. We find a Pomeranchuk effect much larger then in solid helium. In addition, we find that the exponent for the algebraic decay of the hexatic phase differs significantly from the Kosterilitz-Thouless theory of melting. We search for the existence of mesoscopic phases and find evidence of metastable bubbles but no mesoscopic phase that is stable in equilibrium.

Electrical conductivity of high-pressure liquid hydrogen by quantum Monte Carlo methods.

We compute the electrical conductivity for liquid hydrogen at high pressure using Monte Carlo techniques. The method uses coupled electron-ion Monte Carlo simulations to generate configurations of liquid hydrogen. For each configuration, correlated sampling of electrons is performed in order to calculate a set of lowest many-body eigenstates and current-current correlation functions of the system, which are summed over in the many-body Kubo formula to give ac electrical conductivity. The extrapolated dc conductivity at 3000 K for several densities shows a liquid semiconductor to liquid-metal transition at high pressure. Our results are in good agreement with shock-wave data.

Hartree-Fock ground state of the three-dimensional electron gas.

In 1962, Overhauser showed that within Hartree-Fock (HF) the electron gas is unstable to a spin-density wave state. Determining the true HF ground state has remained a challenge. Using numerical calculations for finite systems and analytic techniques, we study the unrestricted HF ground state of the three-dimensional electron gas. At high density, we find broken spin symmetry states with a nearly constant charge density. Unlike previously discussed spin wave states, the observed wave vector of the spin-density wave is smaller than 2k(F). The broken-symmetry state originates from pairing instabilities at the Fermi surface, a model for which is proposed.

Interplay between magic number stabilities and superfluidity of small parahydrogen clusters.

The interplay between magic number stabilities and superfluidity of small parahydrogen clusters with sizes N=5 to 40 and temperatures 0.5 K

Off-diagonal long-range order in solid 4He.

Measurements of the moment of inertia by Kim and Chan have found that solid (4)He acts like a supersolid at low temperatures. To understand the order in solid 4(He), we have used path integral Monte Carlo simulations to calculate the off-diagonal long-range order (ODLRO) [equivalent to Bose-Einstein condensation (BEC)]. We do not find ODLRO in a defect-free hcp crystal of (4)He at the melting density. We discuss these results in relation to proposed quantum solid trial functions. We conclude that the solid (4)He wave function has correlations which suppress both vacancies and BEC.

Quantum Monte Carlo simulation of the high-pressure molecular-atomic crossover in fluid hydrogen.

A first-order liquid-liquid phase transition in high-pressure hydrogen between molecular and atomic fluid phases has been predicted in computer simulations using ab initio molecular dynamics approaches. However, experiments indicate that molecular dissociation may occur through a continuous crossover rather than a first-order transition. Here we study the nature of molecular dissociation in fluid hydrogen using an alternative simulation technique in which electronic correlation is computed within quantum Monte Carlo methods, the so-called coupled electron-ion Monte Carlo method. We find no evidence for a first-order liquid-liquid phase transition.

Superfluidity of dense 4He in vycor.

We calculate properties of a model of 4He in Vycor using the path integral Monte Carlo method. We find that 4He forms a distinct layered structure with a highly localized first layer, a disordered second layer with some atoms delocalized and able to give rise to the observed superfluid response, and higher layers of nearly perfect crystals. The addition of a single 3He atom was enough to bring down the total superfluidity by blocking the exchange in the second layer. Our results are consistent with the persistent liquid-layer model to explain the observations. Such a model may be relevant to the experiments on bulk solid 4He, if there is a fine network of grain boundaries in those systems.

Accurate, efficient, and simple forces computed with quantum Monte Carlo methods.

Computation of ionic forces using quantum Monte Carlo methods has long been a challenge. We introduce a simple procedure, based on known properties of physical electronic densities, to make the variance of the Hellmann-Feynman estimator finite. We obtain very accurate geometries for the molecules H(2), LiH, CH(4), NH(3), H(2)O, and HF, with a Slater-Jastrow trial wave function. Harmonic frequencies for diatomics are also in good agreement with experiment. An antithetical sampling method is also discussed for additional reduction of the variance.

Ring exchanges and the supersolid phase of 4He.

Using path integral Monte Carlo simulations we calculate exchange frequencies in bulk hcp 4He as atoms undergo ring exchange. We fit the frequencies to a lattice model and examine whether such atoms could become a supersolid, that is, have a nonclassical rotational inertia. We find that the scaling with respect to the number of exchanging atoms is such that superfluid behavior will not be observed in a perfect 4He crystal.

Backflow correlations for the electron gas and metallic hydrogen.

We justify and evaluate backflow three-body wave functions for a two-component system of electrons and protons. Based on the generalized Feynman-Kacs formula, many-body perturbation theory, and band structure calculations, we analyze the use and the analytical form of the backflow function from different points of view. The resulting wave functions are used in variational and diffusion Monte Carlo calculations of the electron gas and of solid and liquid metallic hydrogen. For the electron gas, the purely analytic backflow and three-body form gives lower energies than those of previous calculations. For bcc hydrogen, analytical and optimized backflow-three-body wave functions lead to energies nearly as low as those from using local density approximation orbitals in the trial wave function. However, compared to wave functions constructed from density functional solutions, backflow wave functions have the advantage of only few parameters to estimate, the ability to include easily and accurately electron-electron correlations, and that they can be directly generalized from the crystal to a disordered liquid of protons.

Superfluidity in a doped helium droplet.

Path-integral Monte Carlo calculations of the superfluid density throughout 4He droplets doped with linear impurities are presented. After deriving a local estimator for the superfluid density distribution, we find a decreased superfluid response in the cylindrically symmetric region of the first solvation layer. The helium in this region has a superfluid transition temperature similar to that of a two-dimensional helium system and may be responsible for previously unexplained experimental Q-branch measurements.

Spin polarization of the low-density three-dimensional electron gas.

To determine the state of spin polarization of the three-dimensional electron gas at very low densities and zero temperature, we calculate the energy versus spin polarization using diffusion quantum Monte Carlo methods with backflow wave functions and twist averaged boundary conditions. We find a second-order phase transition to a partially polarized phase at r(s) approximately 50+/-2. The magnetic transition temperature is estimated using an effective mean-field method, the Stoner model.

Bose-Einstein condensation at a helium surface.

A path integral Monte Carlo method was used to calculate the Bose-Einstein condensate fraction at the surface of a helium film at T = 0.77 K, as a function of density. Moving from the center of the slab to the surface, the condensate fraction was found to initially increase with decreasing density to a maximum value of 0.9 before decreasing. Long wavelength density correlations were observed in the static structure factor at the surface of the slab. Finally, a surface dispersion relation was calculated from imaginary-time density-density correlations.

Twist-averaged boundary conditions in continuum quantum Monte Carlo algorithms.

We develop and test Quantum Monte Carlo algorithms that use a"twist" or a phase in the wave function for fermions in periodic boundary conditions. For metallic systems, averaging over the twist results in faster convergence to the thermodynamic limit than periodic boundary conditions for properties involving the kinetic energy and has the same computational complexity. We determine exponents for the rate of convergence to the thermodynamic limit for the components of the energy of coulomb systems. We show results with twist averaged variational Monte Carlo on free particles, the Stoner model and the electron gas using Hartree-Fock, Slater-Jastrow, and three-body and backflow wave function. We also discuss the use of twist averaging in the grand canonical ensemble, and numerical methods to accomplish the twist averaging.

Path integral Monte Carlo simulation of the low-density hydrogen plasma.

Restricted path integral Monte Carlo simulations are used to calculate the equilibrium properties of hydrogen in the density and temperature range of 9.83 x 10(-4)

Simulation of flux lines with columnar pins: bose glass and entangled liquids.

Using path integral Monte Carlo we simulate a 3D system of up to 1000 magnetic flux lines by mapping it onto interacting bosons in (2+1)D. With increasing temperatures we find first order melting from an ordered solid to an entangled liquid signaled by a finite entropy jump and sharp discontinuities of the defect density and the structure factor S(G). For a particular density of strong columnar pins the crystal is transformed into a Bose glass phase with patches of crystalline order disrupted by the trapped vortices at the pinning sites but with no overall positional or orientational order. This glassy phase melts into a defected entangled liquid through a continuous transition.