Ground-state energy and stability limit of small 3He drops

Small and stable drops of 3He atoms can only exist above a minimum number of particles, due to the combination of the 3He atom Fermi statistics and its light mass. An accurate estimation of this minimum number using microscopic theory has been difficult due to the inhomogeneous and fermionic nature of these systems. We present a diffusion Monte Carlo calculation of 3He drops with sizes near the minimum in order to determine the stability threshold. The results show that the minimum self-bound drop is formed by N=30 atoms with preferred orbitals for open shells corresponding to maximum value of the spin.Comment: 5 pages, 4 figure

Structure of metastable 2D liquid helium

We present diffusion Monte Carlo (DMC) results on a new metastable, superfluid phase above the crystal ground state in two-dimensional 4He at densities > 0.065 1/A^2. The state is anisotropic with hexatic orbital order. This implies that the liquid--solid phase transition has two stages: A second order phase transition from the isotropic superfluid to the hexatic superfluid, followed by a first order transition that localizes atoms into the triangular crystal order. This metastable hexatic phase has finite condensate fraction and it provides a natural explanation for the superflow in the supersolid grain boundaries

Layer- and bulk roton excitations of 4He in porous media

We examine the energetics of bulk and layer-roton excitations of 4He in various porous medial such as aerogel, Geltech, or Vycor, in order to find out what conclusions can be drawn from experiments on the energetics about the physisorption mechanism. The energy of the layer-roton minimum depends sensitively on the substrate strength, thus providing a mechanism for a direct measurement of this quantity. On the other hand, bulk-like roton excitations are largely independent of the interaction between the medium and the helium atoms, but the dependence of their energy on the degree of filling reflects the internal structure of the matrix and can reveal features of 4He at negative pressures. While bulk-like rotons are very similar to their true bulk counterparts, the layer modes are not in close relation to two-dimensional rotons and should be regarded as a third, completely independent kind of excitation

Boson localization and universality in YBa2Cu(3-x)M(x)O(7-delta)

We consider a two component mixture of charged fermions on neutralizing background with all sign combinations and arbitrarily small mass ratios. In the two impurity limit for the heavier component we show that the pair forms a bound state for all charge combinations. In the lowest order approximation we derive a closed form expression Veff(r) for the binding potential which has short-range repulsion followed by attraction. In the classical limit, when the mass of embedded particles is large m2 much greater than m, we can calculate from Veff(r) also the cohesive energy E and the bond length R of a metallic crystal such as lithium. The lowest order result is R = 3.1 A, E = -0.9 eV, not entirely different from the experimental result for lithium metal. The same interaction for two holes on a parabolic band with m2 greater than m gives the quantum mechanical bound state which one may interpret as a boson or local pair in the case of high-Te and heavy fermion superconductors. We also show that for compounds of the type YBa2Cu(3 - x)M(x)O(7 - delta) one can understand most of the experimental results for the superconducting and normal states with a single temperature dependent boson breaking function f(T) for each impurity content x governing the decay of bosons into pairing fermions. In the normal state f(T) turns out to be a linear, universal function, independent of the impurity content I and the oxygen content delta. We predict with universality a depression in Tc(x) with slight down bending in agreement with experiment. As a natural consequence of the model the bosons become localized slightly above Tc due to the Wigner crystallization, enhanced with lattice local field minima. The holes remain delocalized with a linearly increasing concentration in the normal state, thus explaining the rising Hall density. The boson localization temperature T(sub BL) shows up as a minimum in the Hall density R(sub ab)(exp -1). We also give explanation for very recently observed scaling of temperature dependent Hall effect in La(2 - x)Sr(x)CuO4

Heavy fermion behavior explained by bosons

Conventional heavy fermion (HF) theories require existence of massive fermions. We show that heavy fermion phenomena can also be simply explained by existence of bosons with moderate mass but temperature dependent concentration below the formation temperature T(sub B), which in turn is close to room temperature. The bosons B(++) are proposed to be in chemical equilibrium with a system of holes h(+): B(++) = h(+) + h(+). This equilibrium is governed by a boson breaking function f(T), which determines the decreasing boson density and the increasing fermion density with increasing temperature. Since HF-compounds are hybridized from minimum two elements, we assume in addition existence of another fermion component h(sub s)(+) with temperature independent density. This spectator component is thought to be the main agent in binding the bosons in analogy with electronic or muonic molecules. Using a linear boson breaking function we can explain temperature dependence of the giant linear specific heat coefficient gamma(T) coming essentially from bosons. The maxima in resistivity, Hall coefficient, and susceptibility are explained by boson localization effects due to the Wigner crystallization. The antiferromagnetic transitions in turn are explained by similar localization of the pairing fermion system when their density n(sub h)(T(sub FL)) becomes lower than n(sub WC), the critical density of Wigner crystallization. The model applies irrespective whether a compound is superconducting or not. The same model explains the occurrence of low temperature antiferromagnetism also in high-T(sub c) superconductors. The double transition in UPt3 is proposed to be due to the transition of the pairing fermion liquid from spin polarized to unpolarized state

Excitations in confined helium

We design models for helium in matrices like aerogel, Vycor or Geltech from a manifestly microscopic point of view. For that purpose, we calculate the dynamic structure function of 4He on Si substrates and between two Si walls as a function of energy, momentum transfer, and the scattering angle. The angle--averaged results are in good agreement with the neutron scattering data; the remaining differences can be attributed to the simplified model used here for the complex pore structure of the materials. A focus of the present work is the detailed identification of coexisting layer modes and bulk--like excitations, and, in the case of thick films, ripplon excitations. Involving essentially two--dimensional motion of atoms, the layer modes are sensitive to the scattering angle.Comment: Phys. Rev. B (2003, in press

Boson and fermion dynamics in quasi-one-dimensional flat band lattices

The difference between boson and fermion dynamics in quasi-one-dimensional lattices is studied with exact simulations of particle motion and by calculating the persistent current in small quantum rings. We consider three different lattices which in the tight binding model exhibit flat bands. The physical realization is considered to be an optical lattice with bosonic or fermionic atoms. The atoms are assumed to interact with a repulsive short range interaction. The different statistics of bosons and fermions causes different dynamics. Spinless fermions are easily trapped in the flat band states due to the Pauli exclusion principle, which prevents them from interacting, while boson are able to push each other out from the flat band states

Pair Excitations and Vertex Corrections in Fermi Fluids

Based on an equations--of--motion approach for time--dependent pair correlations in strongly interacting Fermi liquids, we have developed a theory for describing the excitation spectrum of these systems. Compared to the known ``correlated'' random--phase approximation (CRPA), our approach has the following properties: i) The CRPA is reproduced when pair fluctuations are neglected. ii) The first two energy--weighted sumrules are fulfilled implying a correct static structure. iii) No ad--hoc assumptions for the effective mass are needed to reproduce the experimental dispersion of the roton in 3He. iv) The density response function displays a novel form, arising from vertex corrections in the proper polarisation. Our theory is presented here with special emphasis on this latter point. We have also extended the approach to the single particle self-energy and included pair fluctuations in the same way. The theory provides a diagrammatic superset of the familiar GW approximation. It aims at a consistent calculation of single particle excitations with an accuracy that has previously only been achieved for impurities in Bose liquids.Comment: to be published in: JLTP (2007) Proc. Int. Symp. QFS2006, 1-6 Aug. 2006, Kyoto, Japa

Many-body aspects of positron annihilation in the electron gas

We investigate positron annihilation in electron liquid as a case study for many-body theory, in particular the optimized Fermi Hypernetted Chain (FHNC-EL) method. We examine several approximation schemes and show that one has to go up to the most sophisticated implementation of the theory available at the moment in order to get annihilation rates that agree reasonably well with experimental data. Even though there is basically just one number to look at, the electron-positron pair distribution function at zero distance, it is exactly this number that dictates how the full pair distribution behaves: In most cases, it falls off monotonously towards unity as the distance increases. Cases where the electron-positron pair distribution exhibits a dip are precursors to the formation of bound electron--positron pairs. The formation of electron-positron pairs is indicated by a divergence of the FHNC-EL equations, from this we can estimate the density regime where positrons must be localized. This occurs in our calculations in the range 9.4 <= r_s <=10, where r_s is the dimensionless density parameter of the electron liquid.Comment: To appear in Phys. Rev. B (2003

Mott Transition and Spin Structures of Spin-1 Bosons in Two-Dimensional Optical Lattice at Unit Filling

We study the ground state properties of spin-1 bosons in a two-dimensional optical lattice, by applying a variational Monte Carlo method to the S=1 Bose-Hubbard model on a square lattice at unit filling. A doublon-holon binding factor introduced in the trial state provides a noticeable improvement in the variational energy over the conventional Gutzwiller wave function and allows us to deal effectively with the inter-site correlations of particle densities and spins. We systematically show how spin-dependent interactions modify the superfluid-Mott insulator transitions in the S=1 Bose-Hubbard model due to the interplay between the density and spin fluctuations of bosons. Furthermore, regarding the magnetic phases in the Mott region, the calculated spin structure factor elucidates the emergence of nematic and ferromagnetic spin orders for antiferromagnetic ($U_2>0$) and ferromagnetic ($U_2<0$) couplings, respectively.Comment: 5 pages, 5 figures, to appear in Journal of the Physical Society of Japa