An efficient method for the Quantum Monte Carlo evaluation of the static density-response function of a many-electron system

In a recent Letter we introduced Hellmann-Feynman operator sampling in diffusion Monte Carlo calculations. Here we derive, by evaluating the second derivative of the total energy, an efficient method for the calculation of the static density-response function of a many-electron system. Our analysis of the effect of the nodes suggests that correlation is described correctly and we find that the effect of the nodes can be dealt with

The sonic analogue of black hole radiation

A microscopic description of Hawking radiation in sonic black holes has been recently presented (Giovanazzi S 2005 Phys. Rev. Lett. 94 061302). This exactly solvable model is formulated in terms of one-dimensional scattering of a Fermi gas. In this paper, the model is extended to account possible finite size effects of a realistic geometry. The flow of particles is maintained by a piston (i.e. an impenetrable barrier) moving slowly towards the sonic horizon. Using existing technologies the Hawking temperature can be of order of a few microkelvin in a realistic experiment.Comment: 14 pages, 7 figures, submitted to Journal of Physics B: Atomic, Molecular & Optical Physic

Spin fluctuations, susceptibility and the dipole oscillation of a nearly ferromagnetic Fermi gas

We discuss the spin fluctuations and the role played by the magnetic susceptibility in an atomic Fermi gas interacting with positive scattering length. Both thermal and quantum fluctuations are considered. Using a sum rule approach and recent {\it ab initio} Monte Carlo results for the magnetic susceptibility of uniform matter we provide explicit predictions for the frequency of the spin dipole oscillation of a gas trapped by a harmonic potential and discuss the deviations from the behaviour of an ideal gas when the system approaches the ferromagnetic transition. The role of the Landau's parameters in the characterization of the magnetic properties is also discussed.Comment: 5 pages, 2 figure

Isospin and density dependences of nuclear matter symmetry energy coefficients II

Symmetry energy coefficients of explicitly isospin asymmetric nuclear matter at variable densities (from .5$\rho_0$ up to 2 $\rho_0$) are studied as generalized screening functions. An extended stability condition for asymmetric nuclear matter is proposed. We find the possibility of obtaining stable asymmetric nuclear matter even in some cases for which the symmetric nuclear matter limit is unstable. Skyrme-type forces are extensively used in analytical expressions of the symmetry energy coefficients derived as generalized screening functions in the four channels of the particle hole interaction producing alternative behaviors at different $\rho$ and $b$ (respectively the density and the asymmetry coefficient). The spin and spin-isospin coefficients, with corrections to the usual Landau Migdal parameters, indicate the possibility of occurring instabilities with common features depending on the nuclear density and n-p asymmetry. Possible relevance for high energy heavy ions collisions and astrophysical objects is discussed.Comment: 16 pages (latex) plus twelve figures in four eps files, to be published in I.J.M.P.

The correlation energy functional within the GW-RPA approximation: exact forms, approximate forms and challenges

In principle, the Luttinger-Ward Green's function formalism allows one to compute simultaneously the total energy and the quasiparticle band structure of a many-body electronic system from first principles. We present approximate and exact expressions for the correlation energy within the GW-RPA approximation that are more amenable to computation and allow for developing efficient approximations to the self-energy operator and correlation energy. The exact form is a sum over differences between plasmon and interband energies. The approximate forms are based on summing over screened interband transitions. We also demonstrate that blind extremization of such functionals leads to unphysical results: imposing physical constraints on the allowed solutions (Green's functions) is necessary. Finally, we present some relevant numerical results for atomic systems.Comment: 3 figures and 3 tables, under review at Physical Review

Stability of the shell structure in 2D quantum dots

We study the effects of external impurities on the shell structure in semiconductor quantum dots by using a fast response-function method for solving the Kohn-Sham equations. We perform statistics of the addition energies up to 20 interacting electrons. The results show that the shell structure is generally preserved even if effects of high disorder are clear. The Coulomb interaction and the variation in ground-state spins have a strong effect on the addition-energy distributions, which in the noninteracting single-electron picture correspond to level statistics showing mixtures of Poisson and Wigner forms.Comment: 7 pages, 8 figures, submitted to Phys. Rev.

Nonlinear screening of charge impurities in graphene

It is shown that a ``vacuum polarization'' induced by Coulomb potential in graphene leads to a strong suppression of electric charges even for undoped case (no charge carriers). A standard linear response theory is therefore not applicable to describe the screening of charge impurities in graphene. In particular, it overestimates essentially the contributions of charge impurities into the resistivity of graphene.Comment: 3 pages, 1 figure; final version as published in the journa

Non-Adiabatic Spin Transfer Torque in Real Materials

The motion of simple domain walls and of more complex magnetic textures in the presence of a transport current is described by the Landau-Lifshitz-Slonczewski (LLS) equations. Predictions of the LLS equations depend sensitively on the ratio between the dimensionless material parameter $\beta$ which characterizes non-adiabatic spin-transfer torques and the Gilbert damping parameter $\alpha$. This ratio has been variously estimated to be close to 0, close to 1, and large compared to 1. By identifying $\beta$ as the influence of a transport current on $\alpha$, we derive a concise, explicit and relatively simple expression which relates $\beta$ to the band structure and Bloch state lifetimes of a magnetic metal. Using this expression we demonstrate that intrinsic spin-orbit interactions lead to intra-band contributions to $\beta$ which are often dominant and can be (i) estimated with some confidence and (ii) interpreted using the "breathing Fermi surface" model.Comment: 18 pages, 9 figures; submitted to Phys. Rev.

Comment on "Diffusion Monte Carlo study of jellium surfaces: Electronic densities and pair correlation functions"

In a fixed-node diffusion Monte Carlo calculation of the total energy of jellium slabs, Acioli and Ceperley [Phys. Rev. B {\bf 54}, 17199 (1996)] reported jellium surface energies that at low electron densities were significantly higher than those predicted in the local-density approximation (LDA) of density-functional theory. Assuming that the fixed-node error in the slab and the bulk calculations cancel out, we show that their data yield surface energies that are considerably closer to the LDA and in reasonable agreement with those obtained in the random-phase approximation.Comment: 3 pages, 2 figures, to appear in Phys. Rev.

Dynamic spin response of a strongly interacting Fermi gas

We present an experimental investigation of the dynamic spin response of a strongly interacting Fermi gas using Bragg spectroscopy. By varying the detuning of the Bragg lasers, we show that it is possible to measure the response in the spin and density channels separately. At low Bragg energies, the spin response is suppressed due to pairing, whereas the density response is enhanced. These experiments provide the first independent measurements of the spin-parallel and spin-antiparallel dynamic and static structure factors and open the way to a complete study of the structure factors at any momentum. At high momentum the spin-antiparallel dynamic structure factor displays a universal high frequency tail, proportional to $\omega^{-5/2}$, where $\hbar \omega$ is the probe energy.Comment: Replaced with final versio