Beyond the locality approximation in the standard diffusion Monte Carlo method

We present a way to include non local potentials in the standard Diffusion Monte Carlo method without using the locality approximation. We define a stochastic projection based on a fixed node effective Hamiltonian, whose lowest energy is an upper bound of the true ground state energy, even in the presence of non local operators in the Hamiltonian. The variational property of the resulting algorithm provides a stable diffusion process, even in the case of divergent non local potentials, like the hard-core pseudopotentials. It turns out that the modification required to improve the standard Diffusion Monte Carlo algorithm is simple.Comment: 4 pages, 3 figures, to appear in Physical Review

Charge and spin correlations of a one dimensional electron gas on the continuum

We present a variational Monte Carlo study of a model one dimensional electron gas on the continuum, with long-range interaction (1/r decay). At low density the reduced dimensionality brings about pseudonodes of the many-body wavefunction, yielding non-ergodic behavior of naive Monte Carlo sampling, which affects the evaluation of pair correlations and the related structure factors. The problem is however easily solved and we are able to carefully analyze the structure factors obtained from an optimal trial function, finding good agreement with the exact predictions for a Luttinger-like hamiltonian with an interaction similar to the one used in the present study.Comment: 4 pages, 3 figure

Dynamical screening effects in correlated materials: plasmon satellites and spectral weight transfers from a Green's function ansatz to extended dynamical mean field theory

Dynamical screening of the Coulomb interactions in correlated electron systems results in a low-energy effective problem with a dynamical Hubbard interaction U(omega). We propose a Green's function ansatz for the Anderson impurity problem with retarded interactions, in which the Green's function factorizes into a contribution stemming from an effective static-U problem and a bosonic high-energy part introducing collective plasmon excitations. Our approach relies on the scale separation of the low-energy properties, related to the instantaneous static U, from the intermediate to high energy features originating from the retarded part of the interaction. We argue that for correlated materials where retarded interactions arise from downfolding higher-energy degrees of freedom, the characteristic frequencies are typically in the antiadiabatic regime. In this case, accurate approximations to the bosonic factor are relatively easy to construct, with the most simple being the boson factor of the dynamical atomic limit problem. We benchmark the quality of our method against numerically exact continuous time quantum Monte Carlo results for the Anderson-Holstein model both, at half- and quarter-filling. Furthermore we study the Mott transition within the Hubbard-Holstein model within extended dynamical mean field theory. Finally, we apply our technique to a realistic three-band Hamiltonian for SrVO3. We show that our approach reproduces both, the effective mass renormalization and the position of the lower Hubbard band by means of a dynamically screened U, previously determined ab initio within the constrained random phase approximation. Our approach could also be used within schemes beyond dynamical mean field theory, opening a quite general way of describing satellites and plasmon excitations in correlated materials.Comment: 13 pages, 11 figure

Local and non-local electron-phonon couplings in K3Picene and the effect of metallic screening

We analyze the properties of electron-phonon couplings in K3Picene by exploiting a molecular orbital representation derived in the maximally localized Wannier function formalism. This allows us to go beyond the analysis done in Phys. Rev. Lett. 107, 137006 (2011), and separate not only the intra- and intermolecular phonon contributions but also the local and non-local electronic states in the electron-phonon matrix elements. Despite the molecular nature of the crystal, we find that the purely molecular contributions (Holstein-like couplings where the local deformation potential is coupled to intramolecular phonons) account for only 20% of the total electron-phonon interaction lambda. In particular, the Holstein-like contributions to lambda in K3Picene are four times smaller than those computed for an isolated neutral molecule, as they are strongly screened by the metallic bands of the doped crystal. Our findings invalidate the use of molecular electron-phonon calculations to estimate the total electron-phonon coupling in metallic picene, and possibly in other doped metallic molecular crystals. The major contribution (80%) to lambda in K3Picene comes from non-local couplings due to phonon modulated hoppings. We show that the crystal geometry together with the molecular picene structure leads to a strong 1D spatial anisotropy of the non-local couplings. Finally, based on the parameters derived from our density functional theory calculations, we propose a lattice modelization of the electron-phonon couplings in K3Picene which gives 90% of ab-initio lambda.Comment: 13 pages, 8 figures, 3 table

Ground state properties of the one dimensional Coulomb gas

We study the ground state properties of a quasi one dimensional electron gas, interacting via an effective potential with a harmonic transversal confinement and long range Coulomb tail. The exact correlation energy has been calculated for a wide range of electron densities by using the lattice regularized diffusion Monte Carlo method, which is a recent development of the standard projection Monte Carlo technique. In this case it is particularly useful as it allows to sample the exact ground state of the system, even in the low density regime when the exchange between electrons is extremely small. For different values of the width parameter b (0.1 a*_0 <= b <= 4 a*_0), we give a simple parametrization of the correlation energy, which provides an accurate local density energy functional for quasi one dimensional systems. Moreover we show that static correlations are in qualitative agreement with those obtained for the Luttinger liquid model with long range interactions.Comment: 15 pages, 13 figures, to appear in Phys. Rev.

A consistent description of the iron dimer spectrum with a correlated single-determinant wave function

We study the iron dimer by using an accurate ansatz for quantum chemical calculations based on a simple variational wave function, defined by a single geminal expanded in molecular orbitals and combined with a real space correlation factor. By means of this approach we predict that, contrary to previous expectations, the neutral ground state is $^7 \Delta$ while the ground state of the anion is $^8 \Sigma_g^-$, hence explaining in a simple way a long standing controversy in the interpretation of the experiments. Moreover, we characterize consistently the states seen in the photoemission spectroscopy by Leopold \emph{et al.}. It is shown that the non-dynamical correlations included in the geminal expansion are relevant to correctly reproduce the energy ordering of the low-lying spin states.Comment: 5 pages, submitted to the Chemical Physics Letter

Size-consistent variational approaches to non-local pseudopotentials: standard and lattice regularized diffusion Monte Carlo methods revisited

We propose improved versions of the standard diffusion Monte Carlo (DMC) and the lattice regularized diffusion Monte Carlo (LRDMC) algorithms. For the DMC method, we refine a scheme recently devised to treat non-local pseudopotential in a variational way. We show that such scheme --when applied to large enough systems-- maintains its effectiveness only at correspondingly small enough time-steps, and we present two simple upgrades of the method which guarantee the variational property in a size-consistent manner. For the LRDMC method, which is size-consistent and variational by construction, we enhance the computational efficiency by introducing (i) an improved definition of the effective lattice Hamiltonian which remains size-consistent and entails a small lattice-space error with a known leading term, and (ii) a new randomization method for the positions of the lattice knots which requires a single lattice-space.Comment: 10 pages, 4 figures, submitted to the Journal of Chemical Physic

Computation of the Modes of Elliptic Waveguides with a Curvilinear 2D Frequency-Domain Finite-Difference Approach

A scalar Frequency-Domain Finite-Difference approach to the mode computation of elliptic waveguides is presented. The use of an elliptic cylindrical grid allows us to take exactly into account the curved boundary of the structure and a single mesh has been used both for TE and TM modes. As a consequence, a high accuracy is obtained with a reduced computational burden, since the resulting matrix is highly sparse

Resonating Valence Bond wave function: from lattice models to realistic systems

Although mean field theories have been very successful to predict a wide range of properties for solids, the discovery of high temperature superconductivity in cuprates supported the idea that strongly correlated materials cannot be qualitatively described by a mean field approach. After the original proposal by Anderson, there is now a large amount of numerical evidence that the simple but general resonating valence bond (RVB) wave function contains just those ingredients missing in uncorrelated theories, so that the main features of electron correlation can be captured by the variational RVB approach. Strongly correlated antiferromagnetic (AFM) systems, like Cs2CuCl4, displaying unconventional features of spin fractionalization, are also understood within this variational scheme. From the computational point of view the remarkable feature of this approach is that several resonating valence bonds can be dealt simultaneously with a single determinant, at a computational cost growing with the number of electrons similarly to more conventional methods, such as Hartree-Fock or Density Functional Theory. Recently several molecules have been studied by using the RVB wave function; we have always obtained total energies, bonding lengths and binding energies comparable with more demanding multi configurational methods, and in some cases much better than single determinantal schemes. Here we present the paradigmatic case of benzene.Comment: 14 pages, 4 figures. Proceedings of the Conference on Computational Physics CCP2004. To appear in Computer Physics Communication