1
|
Boninsegni M. The Solid Phase of 4He: A Monte Carlo Simulation Study. ENTROPY (BASEL, SWITZERLAND) 2023; 25:1114. [PMID: 37628144 PMCID: PMC10453513 DOI: 10.3390/e25081114] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/26/2023] [Revised: 07/21/2023] [Accepted: 07/21/2023] [Indexed: 08/27/2023]
Abstract
The thermodynamics of solid (hcp) 4He is studied theoretically by means of unbiased Monte Carlo simulations at finite temperature, in a wide range of density. This study complements and extends previous theoretical work, mainly by obtaining results at significantly lower temperatures (down to 60 mK) and for systems of greater size, by including in full the effect of quantum statistics, and by comparing estimates yielded by different pair potentials. All the main thermodynamic properties of the crystal, e.g., the kinetic energy per atom, are predicted to be essentially independent of temperature below ∼ 1 K. Quantum-mechanical exchanges are virtually non-existent in this system, even at the lowest temperature considered. However, effects of quantum statistics are detectable in the momentum distribution. Comparison with available measurements shows general agreement within the experimental uncertainties.
Collapse
Affiliation(s)
- Massimo Boninsegni
- Department of Physics, University of Alberta, Edmonton, AB T6G 2H5, Canada
| |
Collapse
|
2
|
Real Space Triplets in Quantum Condensed Matter: Numerical Experiments Using Path Integrals, Closures, and Hard Spheres. ENTROPY 2020; 22:e22121338. [PMID: 33266522 PMCID: PMC7759805 DOI: 10.3390/e22121338] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/28/2020] [Revised: 11/20/2020] [Accepted: 11/21/2020] [Indexed: 11/17/2022]
Abstract
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).
Collapse
|
3
|
Herrero CP, Ramírez R. Path-integral simulation of solids. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2014; 26:233201. [PMID: 24810944 DOI: 10.1088/0953-8984/26/23/233201] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/16/2023]
Abstract
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.
Collapse
Affiliation(s)
- C P Herrero
- Departamento de Teoria y Simulation de Materiales Instituto de Ciencia de Materiales de Madrid, Consejo Superior de Investigaciones Científicas (CSIC), Campus de Cantoblanco, 28049 Madrid, Spain
| | | |
Collapse
|
4
|
Cencek W, Patkowski K, Szalewicz K. Full-configuration-interaction calculation of three-body nonadditive contribution to helium interaction potential. J Chem Phys 2009; 131:064105. [DOI: 10.1063/1.3204319] [Citation(s) in RCA: 41] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
5
|
Ramírez R, Herrero CP. Quantum path-integral study of the phase diagram and isotope effects of neon. J Chem Phys 2008; 129:204502. [DOI: 10.1063/1.3023036] [Citation(s) in RCA: 26] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
6
|
Sesé LM. Computational study of the structures of gaseous helium-3 at low temperature. J Phys Chem B 2008; 112:10241-54. [PMID: 18651763 DOI: 10.1021/jp802542u] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
A computational study of gaseous helium-3 at T=5.23 K, for number densities rho N<0.0021 A(-3), analyzing the different pair and triplet structures in both the r and the k spaces, is presented. Structures in r space (i.e., instantaneous, total continuous linear response, and centroids) are determined via path-integral Monte Carlo simulations in the canonical ensemble by utilizing the Aziz-Slaman and the SAPT2 interatomic potentials. Additional results obtained with the application of two closures for triplets in r space, the Kirkwood superposition approximation and the Jackson-Feenberg convolution, are also reported. Besides, an analysis of the nonsuitability of quantum hard spheres for describing this system is included. The pair structures in k space are fixed via Ornstein-Zernike schemes appropriate for dealing with quantum diffraction effects in fluids. The effect on the quality of the computed isothermal compressibilities brought about by increasing the sample size in the simulations and by the subsequent application of a grand ensemble correction to the asymptotic behavior of the canonical pair radial correlation functions is also investigated. Furthermore, it is demonstrated analytically that the methods of classical statistical mechanics for dealing with the higher-order direct correlation functions remain fully valid for studying the higher-order correlations of path-integral centroids. By taking advantage of this result, the triplet structure factors for the centroid (also for the instantaneous) correlations are computed by following a number of distinct closures and methods that involve triplet direct correlation functions. The latter computations are intended to explore an alternative scheme to the very expensive fixing of triplet structure factors through direct path integral simulations, an alternative which is expected to yield the main features of these triplet quantities for this gas. Comparison with experiment is made wherever possible, and the results presented allow one to explain the substantial structural features existing in gaseous helium-3.
Collapse
Affiliation(s)
- Luis M Sesé
- Departamento de Ciencias y Técnicas Fisicoquímicas, Facultad de Ciencias, Universidad Nacional de Educación a Distancia, Madrid, Spain.
| |
Collapse
|
7
|
Ramírez R, Herrero CP, Antonelli A, Hernández ER. Path integral calculation of free energies: Quantum effects on the melting temperature of neon. J Chem Phys 2008; 129:064110. [DOI: 10.1063/1.2966006] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
8
|
Sesé LM. Feynman-Hibbs potentials and path integrals for quantum Lennard-Jones systems: Theory and Monte Carlo simulations. Mol Phys 2008. [DOI: 10.1080/00268979500101571] [Citation(s) in RCA: 81] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
9
|
Sesé LM. Computational study of the melting-freezing transition in the quantum hard-sphere system for intermediate densities. I. Thermodynamic results. J Chem Phys 2007; 126:164508. [PMID: 17477615 DOI: 10.1063/1.2718523] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The points where the fluid-solid (face-centered-cubic) transition takes place in the quantum hard-sphere system, for reduced densities 0.85>rhoN*>0.5 (reduced de Broglie wavelengths lambdaB*<or=0.8), have been determined via calculations of Helmholtz free energies. A number of complementary methods have been utilized, namely, path-integral Monte Carlo simulations for fixing the basic thermodynamic and structural quantities, Ornstein-Zernike computations of the fluid isothermal compressibilities using the centroid correlations, and applications of the Einstein crystal technique. Attention is paid to the evaluation of the statistical uncertainties in the isothermal compressibilities and also to the quantum implementation of the Einstein crystal technique by including explicitly the constraint of fixed center of mass. The equation of state along the fluid lambdaB* branches studied has been determined with two methods, one based on the isothermal compressibilities and the other on the usual virial estimator. Along the solid lambdaB* branches the equation of state has been fixed with the virial estimator. The results indicate that the phase transition investigated is governed by entropic effects and that the fluid-solid coexistence densities are arranged along a straight line rhoFCC*=rho(rhoF*), a behavior which at least holds even for lambdaB*<2, as revealed by completing the present analysis with data available in the literature.
Collapse
Affiliation(s)
- Luis M Sesé
- Departamento de Ciencias y Técnicas Fisicoquímicas, Facultad de Ciencias, Universidad Nacional de Educación a Distancia, Paseo Senda del Rey 9, 28040 Madrid, Spain
| |
Collapse
|
10
|
Ujevic S, Vitiello SA. Analysis of the contributions of three-body potentials in the equation of state of 4He. J Chem Phys 2003. [DOI: 10.1063/1.1611872] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022] Open
|
11
|
|
12
|
Chang SY, Boninsegni M. Ab initio potentials and the equation of state of condensed helium at high pressure. J Chem Phys 2001. [DOI: 10.1063/1.1386657] [Citation(s) in RCA: 30] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
13
|
Shinoda K, Miura S, Okazaki S. A molecular approach to quantum fluids based on a generalized Ornstein–Zernike integral equation. J Chem Phys 2001. [DOI: 10.1063/1.1359480] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
14
|
Larrimore L, McFarland RN, Sterne PA, Bug ALR. A two-chain path integral model of positronium. J Chem Phys 2000. [DOI: 10.1063/1.1323979] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
15
|
Chakravarty C, Lynden-Bell RM. Landau free energy curves for melting of quantum solids. J Chem Phys 2000. [DOI: 10.1063/1.1316105] [Citation(s) in RCA: 18] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
16
|
Moroni S, Pederiva F, Fantoni S, Boninsegni M. Equation of state of solid 3He. PHYSICAL REVIEW LETTERS 2000; 84:2650-2653. [PMID: 11017291 DOI: 10.1103/physrevlett.84.2650] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/21/1999] [Indexed: 05/23/2023]
Abstract
We present results of diffusion Monte Carlo calculations for the bcc and hcp phases of solid 3He, using a recent ab initio interatomic potential, including two- and three-body terms. This potential is found to yield an equation of state for condensed 4He in excellent agreement with experiment, in a wide density range. For 3He, we find a systematic discrepancy, worth 0.7 K, between our computed equation of state and a commonly accepted experimental one. We attribute such a discrepancy to an improper choice of reference energy in the determination of the experimental equation of state.
Collapse
Affiliation(s)
- S Moroni
- INFM, Universita di Roma La Sapienza, Roma, Italy
| | | | | | | |
Collapse
|
17
|
Bafile U, Zoppi M, Barocchi F, Magli R, Mayers J. Density evolution of the kinetic energy of fluid and solid 4He at T=6.1 K. PHYSICAL REVIEW. B, CONDENSED MATTER 1996; 54:11969-11972. [PMID: 9985044 DOI: 10.1103/physrevb.54.11969] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/12/2023]
|
18
|
Ceperley DM, Simmons RO, Blasdell RC. Kinetic Energy of Liquid and Solid 4He. PHYSICAL REVIEW LETTERS 1996; 77:115-118. [PMID: 10061785 DOI: 10.1103/physrevlett.77.115] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/23/2023]
|
19
|
Müser MH. The Path-Integral Monte Carlo of Rigid Linear Molecules in Three Dimensions. MOLECULAR SIMULATION 1996. [DOI: 10.1080/08927029608024102] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
20
|
Schell N, Simmons RO, Kaprolat A, Schülke W, Burkel E. Electronic excitations in hcp 4He at 61.5 MPa and 4.3 K studied by inelastic x-ray scattering spectroscopy. PHYSICAL REVIEW LETTERS 1995; 74:2535-2538. [PMID: 10057952 DOI: 10.1103/physrevlett.74.2535] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/23/2023]
|
21
|
Müser MH, Nielaba P, Binder K. Path-integral Monte Carlo study of crystalline Lennard-Jones systems. PHYSICAL REVIEW. B, CONDENSED MATTER 1995; 51:2723-2731. [PMID: 9979043 DOI: 10.1103/physrevb.51.2723] [Citation(s) in RCA: 30] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/12/2023]
|
22
|
Pierleoni C, Ceperley DM, Bernu B, Magro WR. Equation of state of the hydrogen plasma by path integral Monte Carlo simulation. PHYSICAL REVIEW LETTERS 1994; 73:2145-2149. [PMID: 10056984 DOI: 10.1103/physrevlett.73.2145] [Citation(s) in RCA: 28] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/23/2023]
|