Path integral Monte Carlo study of quantum-hard sphere solids

Real Space Triplets in Quantum Condensed Matter: Numerical Experiments Using Path Integrals, Closures, and Hard Spheres

Path integral Monte Carlo and closure computations are utilized to study real space triplet correlations in the quantum hard-sphere system. The conditions cover from the normal fluid phase to the solid phases face-centered cubic (FCC) and cI16 (de Broglie wavelengths , densities ). The focus is on the equilateral and isosceles features of the path-integral centroid and instantaneous structures. Complementary calculations of the associated pair structures are also carried out to strengthen structural identifications and facilitate closure evaluations. The three closures employed are Kirkwood superposition, Jackson-Feenberg convolution, and their average (AV3). A large quantity of new data are reported, and conclusions are drawn regarding (i) the remarkable performance of AV3 for the centroid and instantaneous correlations, (ii) the correspondences between the fluid and FCC salient features on the coexistence line, and (iii) the most conspicuous differences between FCC and cI16 at the pair and the triplet levels at moderately high densities (. This research is expected to provide low-temperature insights useful for the future related studies of properties of real systems (e.g., helium, alkali metals, and general colloidal systems).

Computation of static quantum triplet structure factors of liquid para-hydrogen

The instantaneous and centroid triplet structure factors, S ( 3 ) ( k 1 , k 2 ) , of liquid (one-center) para-hydrogen are computed on the crystallization line for temperatures T/K ≤ 33. The focus is on salient equilateral and isosceles features, and the methods utilized are path integral Monte Carlo (PIMC) simulations and Ornstein-Zernike (OZ) integral equations, which involve Jackson-Feenberg convolution (JF3) and other distinct closures. Long path integral simulation runs are carried out in the canonical ensemble, so as to obtain sufficiently accurate direct PI triplet results. Conclusions are drawn regarding general triplet structure features and the role and usefulness of the OZ closures. The equilateral features are studied in more detail, and one finds that (a) PIMC results point to the existence of regularity in the centroid main peak amplitudes; (b) some of the studied closures give qualitative descriptions for wave numbers below k ≈ 1 Å^-1, but they all fail to describe the main peak amplitude regions (1.75 < k/Å^-1 < 2.5); and (c) JF3 plays the role of a limit closure that is valid for increasing wave numbers (k ≥ 2.6 Å^-1). In addition, representative isosceles PI features turn out to be reasonably bounded (within Δk = 0.1 Å^-1) by those of some closures.

Path integral Monte Carlo on a lattice. II. Bound states

The equilibrium properties of a single quantum particle (qp) interacting with a classical gas for a wide range of temperatures that explore the system's behavior in the classical as well as in the quantum regime is investigated. Both the qp and the atoms are restricted to sites on a one-dimensional lattice. A path integral formalism developed within the context of the canonical ensemble is utilized, where the qp is represented by a closed, variable-step random walk on the lattice. Monte Carlo methods are employed to determine the system's properties. To test the usefulness of the path integral formalism, the Metropolis algorithm is employed to determine the equilibrium properties of the qp in the context of a square well potential, forcing the qp to occupy bound states. We consider a one-dimensional square well potential where all atoms on the lattice are occupied with one atom with an on-site potential except for a contiguous set of sites of various lengths centered at the middle of the lattice. Comparison of the potential energy, the energy fluctuations, and the correlation function are made between the results of the Monte Carlo simulations and the numerical calculations.

Path-integral and Ornstein-Zernike study of quantum fluid structures on the crystallization line

Liquid neon, liquid para-hydrogen, and the quantum hard-sphere fluid are studied with path integral Monte Carlo simulations and the Ornstein-Zernike pair equation on their respective crystallization lines. The results cover the whole sets of structures in the r-space and the k-space and, for completeness, the internal energies, pressures and isothermal compressibilities. Comparison with experiment is made wherever possible, and the possibilities of establishing k-space criteria for quantum crystallization based on the path-integral centroids are discussed. In this regard, the results show that the centroid structure factor contains two significant parameters related to its main peak features (amplitude and shape) that can be useful to characterize freezing.

Path-integral simulation of solids

The path-integral formulation of the statistical mechanics of quantum many-body systems is described, with the purpose of introducing practical techniques for the simulation of solids. Monte Carlo and molecular dynamics methods for distinguishable quantum particles are presented, with particular attention to the isothermal-isobaric ensemble. Applications of these computational techniques to different types of solids are reviewed, including noble-gas solids (helium and heavier elements), group-IV materials (diamond and elemental semiconductors), and molecular solids (with emphasis on hydrogen and ice). Structural, vibrational, and thermodynamic properties of these materials are discussed. Applications also include point defects in solids (structure and diffusion), as well as nuclear quantum effects in solid surfaces and adsorbates. Different phenomena are discussed, as solid-to-solid and orientational phase transitions, rates of quantum processes, classical-to-quantum crossover, and various finite-temperature anharmonic effects (thermal expansion, isotopic effects, electron-phonon interactions). Nuclear quantum effects are most remarkable in the presence of light atoms, so that especial emphasis is laid on solids containing hydrogen as a constituent element or as an impurity.