Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:cond-mat/0702442 · Physical Review B

Quantum Monte Carlo simulation of spin-polarized H

Leandra Vranješ Markić, Jordi Boronat, J. Casulleras

Published Feb 15, 2007• 15 citations
arXivHTMLPDF

Abstract

The ground-state properties of spin polarized hydrogen $\mathrm{H}\ensuremath{\downarrow}$ are obtained by means of diffusion Monte Carlo calculations. Using the most accurate to date ab initio $\mathrm{H}\ensuremath{\downarrow}\text{\ensuremath{-}}\mathrm{H}\ensuremath{\downarrow}$ interatomic potential we have studied its gas phase, from the very dilute regime until densities above its freezing point. At very small densities, the equation of state of the gas is very well described in terms of the gas parameter $\ensuremath{\rho}{a}^{3}$, with $a$ the $s$-wave scattering length. The solid phase has also been studied up to high pressures. The gas-solid phase transition occurs at a pressure of $173\phantom{\rule{0.3em}{0ex}}\mathrm{bar}$, a much higher value than suggested by previous approximate descriptions.

DOI
10.1103/physrevb.75.064506

Related papers

Spin-polarized hydrogen adsorbed on the surface of superfluid 4He

J. M. Marín, Leandra Vranješ Markić, Jordi Boronat

The experimental realization of a thin layer of spin-polarized hydrogen H↓ adsorbed on top of the surface of superfluid (4)He provides one of the best examples of a stable, nearly two-dimensional (2D) quantum Bose gas. We report a theoretical study of this system using quantum Monte Carlo methods in the limit of zero temperature. Using the full Hamiltonian of the system, composed of a superfluid (4)He slab and the adsorbed H↓ layer, we calculate the main properties of its ground state using accurate models for the pair interatomic potentials. Comparing the results for the layer with the ones obtained for a strictly 2D setup, we analyze the departure from the 2D character when the density increases. Only when the coverage is rather small the use of a purely 2D model is justified. The condensate fraction of the layer is significantly larger than in 2D at the same surface density, being as large as 60% at the largest coverage studied.

arXiv:1312.4270•The Journal of Chemical Physics•Dec 13, 2013•1 citation

Monte Carlo simulations of the unitary Bose gas

M. Rossi, Luca Salasnich, Francesco Ancilotto +1 more

We investigate the zero-temperature properties of a diluted homogeneous Bose gas made of $N$ particles interacting via a two-body square-well potential by performing Monte Carlo simulations. We tune the interaction strength to achieve arbitrary positive values of the scattering length and compute by Monte Carlo quadrature the energy per particle $E/N$ and the condensate fraction ${N}_{0}/N$ of this system by using a Jastrow ansatz for the many-body wave function, which avoids the formation of the self-bound ground state and describes instead a (metastable) gaseous state with uniform density. In the unitarity limit, where the scattering length diverges while the range of the interatomic potential is much smaller than the average distance between atoms, we find a finite energy per particle [$E/N=0.70\phantom{\rule{4pt}{0ex}}{\ensuremath{\hbar}}^{2}{(6{\ensuremath{\pi}}^{2}n)}^{2/3}/2m$, with $n$ the number density] and a quite large condensate fraction (${N}_{0}/N=0.83$).

arXiv:1403.5145•Physical Review A•Apr 16, 2014•34 citations

Ground state of small mixed helium and spin-polarized tritium clusters: A quantum Monte Carlo study

Petar Stipanović, Leandra Vranješ Markić, Jordi Boronat +1 more

We report results for the ground-state energy and structural properties of small (4)He-T↓ clusters consisting of up to four T↓ and eight (4)He atoms. These results have been obtained using very well-known (4)He-(4)He and T↓- T↓ interaction potentials and several models for the (4)He- T↓ interatomic potential. All the calculations have been performed with variational and diffusion Monte Carlo methods. It takes at least three atoms to form a mixed bound state. In particular, for small clusters the binding energies are significantly affected by the precise form of the (4)He- T↓ interatomic potential but the stability limits remain unchanged. The only exception is the (4)He(2)T↓ trimer whose stability in the case of the weakest (4)He- T↓ interaction potential is uncertain while it seems stable for other potentials. The mixed trimer (4)He(T↓)(2), a candidate for the Borromean state, is not bound. All other studied clusters are stable. Some of the weakest bound clusters can be classified as quantum halo as a consequence of having high probability of being in a classically forbidden region.

arXiv:1012.3262•The Journal of Chemical Physics•Feb 2, 2011•16 citations

EFFECTIVE RANGE THEORY AND MANY-BODY PERTURBATION THEORY APPLIED TO ELECTRON SPIN-POLARIZED ATOMIC HYDROGEN (H+)

Yea-Hwang Uang, W. C. Stwalley

Accurate scattering lengths and effective ranges for hydrogen atomic collisions are determined from effective range theory. The precise ab initio b3Σ+u potential energy curve values of H2 calculated by Kolos and Wolniewicz and C6, C8, C10 long-range expansion coefficients are used in obtaining accurate representation of the two-body interaction potential. The s-wave scattering length obtained for this potential (A =1/33455 ao) is used to calculate the zero temperature energy as a function of density of electron spin-polarized ("spin-aligned") atomic hydrogen (H+) using previous many-body perturbation theory results for bosons. Neglecting small potential shape dependent terms, one obtains E/N (energy per particle ) as a function of density, which may be compared with similar results obtained from Monte-Carlo or variational calculations. The scattering lengths and effective ranges for the other hydrogen isotopic systems are also presented. It is shown that the first bound state for the b3Σ+u potential (within the Born-Oppenheimer approximation) only exists for a hypothetical isotopic species with mass larger than a value m = 3.258176 g/mole (where the scattering length becomes infinite). In terms of the quantum parameter, η ≡ [MATH]2/ (2µεσ2), two-body binding occurs in the b3Σ+u potential only for η < 0.17094. The second order liquid-to-gas transition (corresponding to the onset of infinite-body binding) occurs at the accurately determined critical value of the quantum parameter ηCB = 0.43891. This corresponds to zero scattering length.

Le Journal de Physique Colloques•Jul 1, 1980•0 citations

Spin-polarized hydrogen and its isotopes: A rich class of quantum phases (Review Article)

I. Bešlić, Leandra Vranješ Markić, Jordi Boronat

We review the recent activity in the theoretical description of spin-polarized atomic hydrogen and its isotopes at very low temperatures. Spin-polarized hydrogen is the only system in nature that remains stable in the gas phase even in the zero temperature limit due to its small mass and weak interatomic interaction. Hydrogen and its heavier isotope tritium are bosons, the heavier mass of tritium producing a self-bound (liquid) system at zero temperature. The other isotope, deuterium, is a fermion with nuclear spin one making possible the study of three different quantum systems depending on the population of the three degenerate spin states. From the theoretical point of view, spin-polarized hydrogen is specially appealing because its interatomic potential is very accurately known making possible its precise quantum many-body study. The experimental study of atomic hydrogen has been very difficult due to its high recombination rate, but it finally led to its Bose-Einstein condensate state in 1998. Degeneracy has also been observed in thin films of hydrogen adsorbed on the 4He surface allowing for the possibility of observing the Berezinskii–Kosterlitz–Thouless superfluid transition.

Low Temperature Physics•Oct 1, 2013•9 citations

Path-integral ground state and superfluid hydrodynamics of a bosonic gas of hard spheres

M. Rossi, Luca Salasnich

We study a bosonic gas of hard spheres by using the exact zero-temperature path-integral ground-state (PIGS) Monte Carlo method and the equations of superfluid hydrodynamics. The PIGS method is implemented to calculate for the bulk system the energy per particle and the condensate fraction through a large range of the gas parameter $n{a}^{3}$ (with $n$ the number density and $a$ the $s$-wave scattering length), going from the dilute gas into the solid phase. The Maxwell construction is then adopted to determine the freezing at $n{a}^{3}=0.264\ifmmode\pm\else\textpm\fi{}0.003$ and the melting at $n{a}^{3}=0.290\ifmmode\pm\else\textpm\fi{}0.003$. In the liquid phase, where the condensate fraction is finite, the equations of superfluid hydrodynamics, based on the PIGS equation of state, are used to find other relevant quantities as a function of the gas parameter: the chemical potential, the pressure, and the sound velocity. In addition, within Feynman's approximation, from the PIGS static structure factor we determine the full excitation spectrum, which displays a maxon-roton behavior when the gas parameter is close to the freezing value. Finally, the equations of superfluid hydrodynamics with the PIGS equation of state are solved for the bosonic system under axially symmetric harmonic confinement obtaining its collective breathing modes.

arXiv:1309.7326•Physical Review A•Nov 14, 2013•33 citations
The ground-state properties of spin polarized hydrogen H\ensuremath↓\mathrm{H}\ensuremath{\downarrow}H\ensuremath↓ are obtained by means of diffusion Monte Carlo calculations. Using the most accurate to date ab initio H\ensuremath↓\ensuremath-H\ensuremath↓\mathrm{H}\ensuremath{\downarrow}\text{\ensuremath{-}}\mathrm{H}\ensuremath{\downarrow}H\ensuremath↓\ensuremath-H\ensuremath↓ interatomic potential we have studied its gas phase, from the very dilute regime until densities above its freezing point. At very small densities, the equation of state of the gas is very well described in terms of the gas parameter \ensuremathρa3\ensuremath{\rho}{a}^{3}\ensuremathρa3, with aaa the sss-wave scattering length. The solid phase has also been studied up to high pressures. The gas-solid phase transition occurs at a pressure of 173bar173\phantom{\rule{0.3em}{0ex}}\mathrm{bar}173bar, a much higher value than suggested by previous approximate descriptions.
The experimental realization of a thin layer of spin-polarized hydrogen H↓ adsorbed on top of the surface of superfluid (4)He provides one of the best examples of a stable, nearly two-dimensional (2D) quantum Bose gas. We report a theoretical study of this system using quantum Monte Carlo methods in the limit of zero temperature. Using the full Hamiltonian of the system, composed of a superfluid (4)He slab and the adsorbed H↓ layer, we calculate the main properties of its ground state using accurate models for the pair interatomic potentials. Comparing the results for the layer with the ones obtained for a strictly 2D setup, we analyze the departure from the 2D character when the density increases. Only when the coverage is rather small the use of a purely 2D model is justified. The condensate fraction of the layer is significantly larger than in 2D at the same surface density, being as large as 60% at the largest coverage studied.
We investigate the zero-temperature properties of a diluted homogeneous Bose gas made of NNN particles interacting via a two-body square-well potential by performing Monte Carlo simulations. We tune the interaction strength to achieve arbitrary positive values of the scattering length and compute by Monte Carlo quadrature the energy per particle E/NE/NE/N and the condensate fraction N0/N{N}_{0}/NN0​/N of this system by using a Jastrow ansatz for the many-body wave function, which avoids the formation of the self-bound ground state and describes instead a (metastable) gaseous state with uniform density. In the unitarity limit, where the scattering length diverges while the range of the interatomic potential is much smaller than the average distance between atoms, we find a finite energy per particle [$E/N=0.70\phantom{\rule{4pt}{0ex}}{\ensuremath{\hbar}}^{2}{(6{\ensuremath{\pi}}^{2}n)}^{2/3}/2m$, with nnn the number density] and a quite large condensate fraction (N0/N=0.83{N}_{0}/N=0.83N0​/N=0.83).
We report results for the ground-state energy and structural properties of small (4)He-T↓ clusters consisting of up to four T↓ and eight (4)He atoms. These results have been obtained using very well-known (4)He-(4)He and T↓- T↓ interaction potentials and several models for the (4)He- T↓ interatomic potential. All the calculations have been performed with variational and diffusion Monte Carlo methods. It takes at least three atoms to form a mixed bound state. In particular, for small clusters the binding energies are significantly affected by the precise form of the (4)He- T↓ interatomic potential but the stability limits remain unchanged. The only exception is the (4)He(2)T↓ trimer whose stability in the case of the weakest (4)He- T↓ interaction potential is uncertain while it seems stable for other potentials. The mixed trimer (4)He(T↓)(2), a candidate for the Borromean state, is not bound. All other studied clusters are stable. Some of the weakest bound clusters can be classified as quantum halo as a consequence of having high probability of being in a classically forbidden region.
Accurate scattering lengths and effective ranges for hydrogen atomic collisions are determined from effective range theory. The precise ab initio b3Σ+u potential energy curve values of H2 calculated by Kolos and Wolniewicz and C6, C8, C10 long-range expansion coefficients are used in obtaining accurate representation of the two-body interaction potential. The s-wave scattering length obtained for this potential (A =1/33455 ao) is used to calculate the zero temperature energy as a function of density of electron spin-polarized ("spin-aligned") atomic hydrogen (H+) using previous many-body perturbation theory results for bosons. Neglecting small potential shape dependent terms, one obtains E/N (energy per particle ) as a function of density, which may be compared with similar results obtained from Monte-Carlo or variational calculations. The scattering lengths and effective ranges for the other hydrogen isotopic systems are also presented. It is shown that the first bound state for the b3Σ+u potential (within the Born-Oppenheimer approximation) only exists for a hypothetical isotopic species with mass larger than a value m = 3.258176 g/mole (where the scattering length becomes infinite). In terms of the quantum parameter, η ≡ [MATH]2/ (2µεσ2), two-body binding occurs in the b3Σ+u potential only for η < 0.17094. The second order liquid-to-gas transition (corresponding to the onset of infinite-body binding) occurs at the accurately determined critical value of the quantum parameter ηCB = 0.43891. This corresponds to zero scattering length.
We review the recent activity in the theoretical description of spin-polarized atomic hydrogen and its isotopes at very low temperatures. Spin-polarized hydrogen is the only system in nature that remains stable in the gas phase even in the zero temperature limit due to its small mass and weak interatomic interaction. Hydrogen and its heavier isotope tritium are bosons, the heavier mass of tritium producing a self-bound (liquid) system at zero temperature. The other isotope, deuterium, is a fermion with nuclear spin one making possible the study of three different quantum systems depending on the population of the three degenerate spin states. From the theoretical point of view, spin-polarized hydrogen is specially appealing because its interatomic potential is very accurately known making possible its precise quantum many-body study. The experimental study of atomic hydrogen has been very difficult due to its high recombination rate, but it finally led to its Bose-Einstein condensate state in 1998. Degeneracy has also been observed in thin films of hydrogen adsorbed on the 4He surface allowing for the possibility of observing the Berezinskii–Kosterlitz–Thouless superfluid transition.
We study a bosonic gas of hard spheres by using the exact zero-temperature path-integral ground-state (PIGS) Monte Carlo method and the equations of superfluid hydrodynamics. The PIGS method is implemented to calculate for the bulk system the energy per particle and the condensate fraction through a large range of the gas parameter na3n{a}^{3}na3 (with nnn the number density and aaa the sss-wave scattering length), going from the dilute gas into the solid phase. The Maxwell construction is then adopted to determine the freezing at $n{a}^{3}=0.264\ifmmode\pm\else\textpm\fi{}0.003$ and the melting at $n{a}^{3}=0.290\ifmmode\pm\else\textpm\fi{}0.003$. In the liquid phase, where the condensate fraction is finite, the equations of superfluid hydrodynamics, based on the PIGS equation of state, are used to find other relevant quantities as a function of the gas parameter: the chemical potential, the pressure, and the sound velocity. In addition, within Feynman's approximation, from the PIGS static structure factor we determine the full excitation spectrum, which displays a maxon-roton behavior when the gas parameter is close to the freezing value. Finally, the equations of superfluid hydrodynamics with the PIGS equation of state are solved for the bosonic system under axially symmetric harmonic confinement obtaining its collective breathing modes.