No Arabic abstract
We report on the calculation of the ground-state atomic kinetic energy, $E_{k}$, and momentum distribution of solid Ne by means of the diffusion Monte Carlo method and Aziz HFD-B pair potential. This approach is shown to perform notably for this crystal since we obtain very good agreement with respect to experimental thermodynamic data. Additionally, we study the structural properties of solid Ne at densities near the equilibrium by estimating the radial pair-distribution function, Lindemanns ratio and atomic density profile around the positions of the perfect crystalline lattice. Our value for $E_{k}$ at the equilibrium density is $41.51(6)$ K, which agrees perfectly with the recent prediction made by Timms {it et al.}, $41(2)$ K, based on their deep-inelastic neutron scattering experiments carried out over the temperature range $4 - 20$ K, and also with previous path integral Monte Carlo results obtained with the Lennard-Jones and Aziz HFD-C2 atomic pairwise interactions. The one-body density function of solid Ne is calculated accurately and found to fit perfectly, within statistical uncertainty, to a Gaussian curve. Furthermore, we analyze the degree of anharmonicity of solid Ne by calculating some of its microscopic ground-state properties within traditional harmonic approaches. We provide insightful comparison to solid $^4$He in terms of the Debye model, in order to size the relevance of anharmonic effects in Ne.
We study the zero-temperature equation of state (EOS) of solid 4He in the hexagonal closed packet (hcp) phase over the 0-57 GPa pressure range by means of the Diffusion Monte Carlo (DMC) method and the semi-empirical Aziz pair potential HFD-B(HE). In the low pressure regime (P ~ 0-1 GPa) we assess excellent agreement with experiments and we give an accurate description of the atomic kinetic energy, Lindemann ratio and Debye temperature over a wide range of molar volumes (22-6 cm^{3}/mol). However, on moving to higher pressures our calculated P-V curve presents an increasingly steeper slope which ultimately provides differences within ~40 % with respect to measurements. In order to account for many-body interactions arising in the crystal with compression which are not reproduced by our model, we perform additional electronic density-functional theory (DFT) calculations for correcting the computed DMC energies in a perturbative way. We explore both generalized gradient and local density approximations (GGA and LDA, respectively) for the electronic exchange-correlation potential. By proceeding in this manner, we show that discrepancies with respect to high pressure data are reduced to 5-10 % with few computational extra cost. Further comparison between our calculated EOSs and ab initio curves deduced for the perfect crystal and corrected for the zero-point motion of the atoms enforces the reliability of our approach.
The momentum distribution function for the two-component 1D gases of bosons and fermions is studied in the limit of strong interatomic repulsion. A pronounced reconstruction of the distribution is found at a temperature much smaller than the Fermi temperature. This new temperature scale, which equals the Fermi temperature divided by the dimensionless coupling strength, is a feature of the two-component model and does not exist in the one-component case. We estimate the parameters relevant for the experimental observation of the crossover effect.
We show that thermalization of the motion of atoms at negative temperature is possible in an optical lattice, for conditions that are feasible in current experiments. We present a method for reversibly inverting the temperature of a trapped gas. Moreover, a negative-temperature ensemble can be cooled, reducing abs(T), by evaporation of the lowest-energy particles. This enables the attainment of the Bose-Einstein condensation phase transition at negative temperature.
Light nuclei at room temperature and below exhibit a kinetic energy which significantly deviates from the predictions of classical statistical mechanics. This quantum kinetic energy is responsible for a wide variety of isotope effects of interest in fields ranging from chemistry to climatology. It also furnishes the second moment of the nuclear momentum distribution, which contains subtle information about the chemical environment and has recently become accessible to deep inelastic neutron scattering experiments. Here we show how, by combining imaginary time path integral dynamics with a carefully designed generalized Langevin equation, it is possible to dramatically reduce the expense of computing the quantum kinetic energy. We also introduce a transient anisotropic Gaussian approximation to the nuclear momentum distribution which can be calculated with negligible additional effort. As an example, we evaluate the structural properties, the quantum kinetic energy, and the nuclear momentum distribution for a first-principles simulation of liquid water.
We study molecular para-hydrogen (p-${rm H_{2}}$) and ortho-deuterium (o-${rm D_{2}}$) in two dimensions and in the limit of zero temperature by means of the diffusion Monte Carlo method. We report energetic and structural properties of both systems like the total and kinetic energy per particle, radial pair distribution function, and Lindemanns ratio in the low pressure regime. By comparing the total energy per particle as a function of the density in liquid and solid p-${rm H_{2}}$, we show that molecular para-hydrogen, and also ortho-deuterium, remain solid at zero temperature. Interestingly, we assess the quality of three different symmetrized trial wave functions, based on the Nosanow-Jastrow model, in the p-${rm H_{2}}$ solid film at the variational level. In particular, we analyze a new type of symmetrized trial wave function which has been used very recently to describe solid $^{4}$He and found that also characterizes hydrogen satisfactorily. With this wave function, we show that the one-body density matrix $varrho_{1} (r)$ of solid p-${rm H_{2}}$ possesses off-diagonal long range order, with a condensate fraction that increases sizably in the negative pressure regime.