Path Integral Calculations of exchange in solid 4He

Recently there have been experimental indications that solid 4He might be a supersolid. We discuss the relation of supersolid behavior to ring exchange. The tunnelling frequencies for ring exchanges in quantum solids are calculated using Path Integral Monte Carlo by finding the free energy for making a path that begins with the atoms in one configuration and ends with a permutation of those positions. We find that the exchange frequencies in solid 4He are described by a simple lattice model which does not show supersolid behavior. Thus, the PIMC calculations constrain the mechanism for the supersolid behavior. We also look at the characteristics of very long exchanges needed for macroscopic mass transport

Critical temperature of the superfluid transition in bose liquids

A phenomenological criterion for the superfluid transition is proposed, which is similar to the Lindemann criterion for the crystal melting. Then we derive a new formula for the critical temperature, relating $T_{\lambda}$ to the mean kinetic energy per particle above the transition. The suppression of the critical temperature in a sufficiently dense liquid is described as a result of the quantum decoherence phenomenon. The theory can account for the observed dependence of $T_{\lambda}$ on density in liquid helium and results in an estimate $T_{\lambda} \sim 1.1$ K for molecular hydrogen.Comment: 4 pages, 1 fi

Proof for an upper bound in fixed-node Monte Carlo for lattice fermions

We justify a recently proposed prescription for performing Green Function Monte Carlo calculations on systems of lattice fermions, by which one is able to avoid the sign problem. We generalize the prescription such that it can also be used for problems with hopping terms of different signs. We prove that the effective Hamiltonian, used in this method, leads to an upper bound for the ground-state energy of the real Hamiltonian, and we illustrate the effectiveness of the method on small systems.Comment: 14 pages in revtex v3.0, no figure

Excitation spectra of a 3He impurity on 4He clusters

The diffusion Monte Carlo technique is used to calculate and analyze the excitation spectrum of a single 3He atom bound to a cluster with N 4He atoms, with the aim of establishing the most adequate filling ordering of single-fermion orbits to the mixed clusters with a large number of 3He atoms. The resulting ordering looks like the rotational spectrum of a diatomic molecule, being classified only by the angular momentum of the level, although vibrational-like excitations appear at higher energies for sufficiently large N

Condensate fraction in liquid 4He at zero temperature

We present results of the one-body density matrix (OBDM) and the condensate fraction n_0 of liquid 4He calculated at zero temperature by means of the Path Integral Ground State Monte Carlo method. This technique allows to generate a highly accurate approximation for the ground state wave function Psi_0 in a totally model-independent way, that depends only on the Hamiltonian of the system and on the symmetry properties of Psi_0. With this unbiased estimation of the OBDM, we obtain precise results for the condensate fraction n_0 and the kinetic energy K of the system. The dependence of n_0 with the pressure shows an excellent agreement of our results with recent experimental measurements. Above the melting pressure, overpressurized liquid 4He shows a small condensate fraction that has dropped to 0.8% at the highest pressure of p = 87 bar.Comment: 12 pages. 4 figures. Accepted for publication on "Journal of Low Temperature Physics

The Debye-Waller Factor in solid 3He and 4He

The Debye-Waller factor and the mean-squared displacement from lattice sites for solid 3He and 4He were calculated with Path Integral Monte Carlo at temperatures between 5 K and 35 K, and densities between 38 nm^(-3) and 67 nm^(-3). It was found that the mean-squared displacement exhibits finite-size scaling consistent with a crossover between the quantum and classical limits of N^(-2/3) and N^(-1/3), respectively. The temperature dependence appears to be T^3, different than expected from harmonic theory. An anisotropic k^4 term was also observed in the Debye-Waller factor, indicating the presence of non-Gaussian corrections to the density distribution around lattice sites. Our results, extrapolated to the thermodynamic limit, agree well with recent values from scattering experiments.Comment: 5 figure

Hydrogen-Helium Mixtures at High Pressure

The properties of hydrogen-helium mixtures at high pressure are crucial to address important questions about the interior of Giant planets e.g. whether Jupiter has a rocky core and did it emerge via core accretion? Using path integral Monte Carlo simulations, we study the properties of these mixtures as a function of temperature, density and composition. The equation of state is calculated and compared to chemical models. We probe the accuracy of the ideal mixing approximation commonly used in such models. Finally, we discuss the structure of the liquid in terms of pair correlation functions.Comment: Proceedings article of the 5th Conference on Cryocrystals and Quantum Crystals in Wroclaw, Poland, submitted to J. Low. Temp. Phys. (2004

Improved tensor-product expansions for the two-particle density matrix

We present a new density-matrix functional within the recently introduced framework for tensor-product expansions of the two-particle density matrix. It performs well both for the homogeneous electron gas as well as atoms. For the homogeneous electron gas, it performs significantly better than all previous density-matrix functionals, becoming very accurate for high densities and outperforming Hartree-Fock at metallic valence electron densities. For isolated atoms and ions, it is on a par with previous density-matrix functionals and generalized gradient approximations to density-functional theory. We also present analytic results for the correlation energy in the low density limit of the free electron gas for a broad class of such functionals.Comment: 4 pages, 2 figure

Incorporation of Density Matrix Wavefunctions in Monte Carlo Simulations: Application to the Frustrated Heisenberg Model

We combine the Density Matrix Technique (DMRG) with Green Function Monte Carlo (GFMC) simulations. The DMRG is most successful in 1-dimensional systems and can only be extended to 2-dimensional systems for strips of limited width. GFMC is not restricted to low dimensions but is limited by the efficiency of the sampling. This limitation is crucial when the system exhibits a so-called sign problem, which on the other hand is not a particular obstacle for the DMRG. We show how to combine the virtues of both methods by using a DMRG wavefunction as guiding wave function for the GFMC. This requires a special representation of the DMRG wavefunction to make the simulations possible within reasonable computational time. As a test case we apply the method to the 2-dimensional frustrated Heisenberg antiferromagnet. By supplementing the branching in GFMC with Stochastic Reconfiguration (SR) we get a stable simulation with a small variance also in the region where the fluctuations due to minus sign problem are maximal. The sensitivity of the results to the choice of the guiding wavefunction is extensively investigated. We analyse the model as a function of the ratio of the next-nearest to nearest neighbor coupling strength. We observe in the frustrated regime a pattern of the spin correlations which is in-between dimerlike and plaquette type ordering, states that have recently been suggested. It is a state with strong dimerization in one direction and weaker dimerization in the perpendicular direction.Comment: slightly revised version with added reference

Annealing Effect for Supersolid Fraction in 4^4He

We report on experimental confirmation of the non-classical rotational inertia (NCRI) in solid helium samples originally reported by Kim and Chan. The onset of NCRI was observed at temperatures below ~400 mK. The ac velocity for initiation of the NCRI suppression is estimated to be ~10 $\mu$m/sec. After an additional annealing of the sample at $T= 1.8$ K for 12 hours, ~ 10% relative increase of NCRI fraction was observed. Then after repeated annealing with the same conditions, the NCRI fraction was saturated. It differs from Reppy's observation on a low pressure solid sample.Comment: to be published in J. of Low Temp. Phys. (QFS2006 proceedings