BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Michele Casula

Michele Casula appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.str-elcond-mat.mtrl-scicond-mat.supr-conphysics.chem-phcond-mat.mes-hallcond-mat.othercond-mat.stat-mechphysics.comp-phquant-ph

Catalog footprint

What is connected

25works
9topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

25 published item(s)

preprint2023arXiv

Thermal dependence of the hydrated proton and optimal proton transfer

Water is a key ingredient for life and plays a central role as solvent in many biochemical reactions. However, the intrinsically quantum nature of the hydrogen nucleus, revealing itself in a large variety of physical manifestations, including proton transfer, gives rise to unexpected phenomena whose description is still elusive. Here we study, by an unprecedented combination of state-of-the-art quantum Monte Carlo methods and path-integral molecular dynamics, the structure and hydrogen-bond dynamics of the protonated water hexamer, the fundamental unit for the hydrated proton. We report a remarkably low thermal expansion of the hydrogen bond from zero temperature up to 300 K, owing to the presence of short-Zundel configurations, characterised by proton delocalisation and favoured by the synergy of nuclear quantum effects and thermal activation. The hydrogen bond strength progressively weakens above 300 K, when localised Eigen-like configurations become relevant. Our analysis, supported by the instanton statistics of shuttling protons, reveals that the near-room-temperature range from 250 K to 300 K is a ``sweet spot'' for proton transfer, and thus for many phenomena depending on it, including life.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Hydrogen phase-IV characterization by full account of quantum anharmonicity

We devise a framework to compute accurate phonons in molecular crystals even in case of strong quantum anharmonicity. Our approach is based on the calculation of the static limit of the phononic Matsubara Green's function from path integral molecular dynamics simulations. Our method enjoys a remarkably low variance, which allows one to compute accurate phonon frequencies after a few picoseconds of nuclear dynamics, and it is further stabilized by the use of appropriate constrained displacement operators. We applied it to solid hydrogen at high pressure. For phase III, our predicted infrared (IR) and Raman active vibrons agree very well with experiments. We then characterize the crystalline symmetry of phase IV by direct comparison with vibrational data and we determine the character of its Raman and IR vibron peaks.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Quantum Rényi entropy by optimal thermodynamic integration paths

Despite being a well-established operational approach to quantify entanglement, Rényi entropy calculations have been plagued by their computational complexity. We introduce here a theoretical framework based on an optimal thermodynamic integration scheme, where the Rényi entropy can be efficiently evaluated using regularizing paths. This approach avoids slowly convergent fluctuating contributions and leads to low-variance estimates. In this way, large system sizes and high levels of entanglement in model or first-principles Hamiltonians are within our reach. We demonstrate it in the one-dimensional quantum Ising model and perform the evaluation of entanglement entropy in the formic acid dimer, by discovering that its two shared protons are entangled even above room temperature.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Ground-state properties of the narrowest zigzag graphene nanoribbon from quantum Monte Carlo and comparison with density functional theory

By means of quantum Monte Carlo (QMC) calculations from first principles, we study the ground-state properties of the narrowest zigzag graphene nanoribbon, with an infinite linear acene structure. We show that this quasi-one-dimensional system is correlated and its ground state is made of localized $π$ electrons whose spins are antiferromagnetically (AFM) ordered. The AFM stablization energy (36(3) meV per carbon atom) and the absolute magnetization (1.13(1) $μ_\textrm{B}$ per unit cell) predicted by QMC are sizable, and they suggest the survival of antiferromagnetic correlations above room temperature. These values can be reproduced to some extent by density functional theory (DFT) only by assuming strong interactions, either within the DFT+U framework or using hybrid functionals. Based on our QMC results, we then provide the strength of Hubbard repulsion in DFT+U suitable for this class of systems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

TurboRVB: a many-body toolkit for {\it ab initio} electronic simulations by quantum Monte Carlo

TurboRVB is a computational package for {\it ab initio} Quantum Monte Carlo (QMC) simulations of both molecular and bulk electronic systems. The code implements two types of well established QMC algorithms: Variational Monte Carlo (VMC), and Diffusion Monte Carlo in its robust and efficient lattice regularized variant. A key feature of the code is the possibility of using strongly correlated many-body wave functions. The electronic wave function (WF) is obtained by applying a Jastrow factor, which takes into account dynamical correlations, to the most general mean-field ground state, written either as an antisymmetrized geminal product with spin-singlet pairing, or as a Pfaffian, including both singlet and triplet correlations. This wave function can be viewed as an efficient implementation of the so-called resonating valence bond (RVB) ansatz, first proposed by L. Pauling and P. W. Anderson in quantum chemistry and condensed matter physics, respectively. The RVB ansatz implemented in TurboRVB has a large variational freedom, including the Jastrow correlated Slater determinant as its simplest, but nontrivial case. Moreover, it has the remarkable advantage of remaining with an affordable computational cost, proportional to the one spent for the evaluation of a single Slater determinant. The code implements the adjoint algorithmic differentiation that enables a very efficient evaluation of energy derivatives, comprising the ionic forces. Thus, one can perform structural optimizations and molecular dynamics in the canonical NVT ensemble at the VMC level. For the electronic part, a full WF optimization is made possible thanks to state-of-the-art stochastic algorithms for energy minimization. The code has been efficiently parallelized by using a hybrid MPI-OpenMP protocol, that is also an ideal environment for exploiting the computational power of modern GPU accelerators.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Photoinduced renormalization of Dirac states in BaNiS$_2$

By means of pump-probe time- and angle-resolved photoelectron spectroscopy, we provide evidence of a sizeable reduction of the Fermi velocity of out-of-equilibrium Dirac bands in the quasi-two-dimensional semimetal BaNiS$_2$. First-principle calculations indicate that this band renormalization is ascribed to a change in non-local electron correlations driven by a photo-induced enhancement of screening properties. This effect is accompanied by a slowing down of the Dirac fermions and by a non-rigid shift of the bands at the center of the Brillouin zone. This result suggests that other similar electronic structure renormalizations may be photoinduced in other materials in presence of strong non-local correlations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Anomalous metallic state in quasi-two-dimensional BaNiS$_{2}$

We report on a systematic study of the thermodynamic, electronic and charge transport properties of high-quality single crystals of BaNiS$_2$, the metallic end-member of the quasi-twodimensional BaCo$_{1-x}$Ni$_x$S$_2$ system characterized by a metal-insulator transition at $x_{cr}=0.22$. Our analysis of magnetoresistivity and specific heat data consistently suggests a picture of compensated semimetal with two hole- and one electron-bands, where electron-phonon scattering dominates charge transport and the minority holes exhibit, below $\sim$100 K, a very large mobility, $μ_h\sim$ 15000 cm$^2$V$^{-1}$s$^{-1}$, which is explained by a Dirac-like band. Evidence of unconventional metallic properties is given by an intriguing crossover of the resistivity from a Bloch-Grüneisen regime to a linear$-T$ regime occurring at 2 K and by a strong linear term in the paramagnetic susceptibility above 100 K. We discuss the possibility that these anomalies reflect a departure from conventional Fermi-liquid properties in presence of short-range AF fluctuations and of a large Hund coupling.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Competing collinear magnetic structures in superconducting FeSe by first principles quantum Monte Carlo calculations

Resolving the interplay between magnetic interactions and structural properties in strongly correlated materials through a quantitatively accurate approach has been a major challenge in condensed matter physics. Here we apply highly accurate first principles quantum Monte Carlo (QMC) techniques to obtain structural and magnetic properties of the iron selenide (FeSe) superconductor under pressure. Where comparable, the computed properties are very close to the experimental values. Of potential ordered magnetic configurations, collinear spin configurations are the most energetically favorable over the explored pressure range. They become nearly degenerate in energy with bicollinear spin orderings at around 7 GPa, when the experimental critical temperature $T_c$ is the highest. On the other hand, ferromagnetic, checkerboard, and staggered dimer configurations become relatively higher in energy as the pressure increases. The behavior under pressure is explained by an accurate analysis of the charge compressibility and the orbital occupation as described by the QMC many-body wave function, which reveals how spin, charge and orbital degrees of freedom are strongly coupled in this compound. This remarkable pressure evolution suggests that stripe-like magnetic fluctuations may be responsible for the enhanced $T_c$ in FeSe and that higher T$_c$ is associated with nearness to a crossover between collinear and bicollinear ordering.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dynamical screening in correlated electron systems - from lattice models to realistic materials

Recent progress in treating the dynamically screened nature of the Coulomb interaction in strongly correlated lattice models and materials is reviewed with a focus on computational schemes based on the dynamical mean field approximation. We discuss approximate and exact methods for the solution of impurity models with retarded interactions, and explain how these models appear as auxiliary problems in various extensions of the dynamical mean field formalism. The current state of the field is illustrated with results from recent applications of these schemes to U-V Hubbard models and correlated materials.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Geminal embedding scheme for optimal atomic basis set construction in correlated calculations

We introduce an efficient method to construct optimal and system adaptive basis sets for use in electronic structure and quantum Monte Carlo calculations. The method is based on an embedding scheme in which a reference atom is singled out from its environment, while the entire system (atom and environment) is described by a Slater determinant or its antisymmetrized geminal power (AGP) extension. The embedding procedure described here allows for the systematic and consistent contraction of the primitive basis set into geminal embedded orbitals (GEOs), with a dramatic reduction of the number of variational parameters necessary to represent the many-body wave function, for a chosen target accuracy. Within the variational Monte Carlo method, the Slater or AGP part is determined by a variational minimization of the energy of the whole system in presence of a flexible and accurate Jastrow factor, representing most of the dynamical electronic correlation. The resulting GEO basis set opens the way for a fully controlled optimization of many-body wave functions in electronic structure calculation of bulk materials, namely, containing a large number of electrons and atoms. We present applications on the water molecule, the volume collapse transition in cerium, and the high-pressure liquid hydrogen.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Downfolding electron-phonon Hamiltonians from ab-initio calculations: application to K$_3$Picene

We propose an electron-phonon parameterization which reliably reproduces the geometry and harmonic frequencies of a real system. With respect to standard electron-phonon models, it adds a "double-counting" correction, which takes into account the lattice deformation as the system is dressed by low-energy electron-phonon processes. We show the importance of this correction by studying potassium-doped picene (K$_3$Picene), recently claimed to be a superconductor with a $T_c$ of up to 18 K. The Hamiltonian parameters are derived from ab-initio density functional theory, and the lattice model is solved by dynamical mean-field theory. Our calculations include the effects of electron-electron interactions and local electron-phonon couplings. Even with the inclusion of a strongly coupled molecular phonon, the Hubbard repulsion prevails and the system is an insulator with a small Mott gap of $\approx$ 0.2 eV.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Dynamical correlations and screened exchange on the experimental bench: spectral properties of the cobalt pnictide BaCo2As2

Understanding the Fermi surface and low-energy excitations of iron or cobalt pnictides is crucial for assessing electronic instabilities such as magnetic or superconducting states. Here, we propose and implement a new approach to compute the low-energy properties of correlated electron materials, taking into account both screened exchange beyond the local density approximation and local dynamical correlations. The scheme allows us to resolve the puzzle of BaCo2As2, for which standard electronic structure techniques predict a ferromagnetic instability not observed in nature.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Asymmetric band widening by screened exchange competing with local correlations in SrVO3: new surprises on an old compound from combined GW and dynamical mean field theory GW+DMFT

The very first dynamical implementation of the combined GW and dynamical mean field scheme "GW+DMFT" for a real material was achieved recently [J.M. Tomczak et al., Europhys. Lett. 100 67001 (2012)], and applied to the ternary transition metal oxide SrVO3. Here, we review and extend that work, giving not only a detailed account of full GW+DMFT calculations, but also discussing and testing simplified approximate schemes. We give insights into the nature of exchange and correlation effects: Dynamical renormalizations in the Fermi liquid regime of SrVO3 are essentially local, and nonlocal correlations mainly act to screen the Fock exchange term. The latter substantially widens the quasi-particle band structure, while the band narrowing induced by the former is accompanied by a spectral weight transfer to higher energies. Most interestingly, the exchange broadening is much more pronounced in the unoccupied part of spectrum. As a result, the GW+DMFT electronic structure of SrVO3 resembles the conventional density functional based dynamical mean field (DFT+DMFT) description for occupied states, but is profoundly modified in the empty part. Our work leads to a reinterpretation of inverse photoemission spectroscopy (IPES) data. Indeed, we assign a prominent peak at about 2.7 eV dominantly to eg states, rather than to an upper Hubbard band of t2g character. Similar surprises can be expected for other transition metal oxides, calling for more detailed investigations of the conduction band states.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Improper s-wave symmetry for the electronic pairing in iron-based superconductors by first-principles calculation

By means of space-group symmetry arguments, we argue that the electronic pairing in iron-based high temperature superconductors shows a structure which is a linear combination of planar s-wave and d-wave symmetry channels, both preserving the 3-dimensional A_1g irreducible representation of the corresponding crystal point-group. We demonstrate that the s- and d-wave channels are determined by the parity under reflection of the electronic orbitals through the iron planes, and by improper rotations around the iron sites. We provide evidence of these general properties by performing accurate quantum Monte Carlo ab-initio calculations of the pairing function, for a FeSe lattice with tetragonal experimental geometry at ambient pressure. We find that this picture survives even in the FeSe under pressure and at low temperatures, when the tetragonal point-group symmetry is slightly broken. In order to achieve a higher resolution in momentum space we introduce a BCS model that faithfully describes our QMC variational pairing function on the simulated 4x4 FeSe unit cell. This allows us to provide a k-resolved image of the pairing function, and show that non-isotropic contributions in the BCS gap function are related to the improper s-wave symmetry. Our theory can rationalize and explain a series of contradictory experimental findings, such as the observation of twofold symmetry in the FeSe superconducting phase, the anomalous drop of T_c with Co-impurity in LaFeAsO_(1-x)F_x, the s-to-d-wave gap transition in BaFe_2As_2 under K doping, and the nodes appearing in the LiFeAs superconducting gap upon P isovalent substitution.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Quantum Monte Carlo study of the protonated water dimer

We report an extensive theoretical study of the protonated water dimer (Zundel ion) by means of the highly correlated variational Monte Carlo and lattice regularized Monte Carlo approaches. This system represents the simplest model for proton transfer (PT) and a correct description of its properties is essential in order to understand the PT mechanism in more complex acqueous systems. Our Jastrow correlated AGP wave function ensures an accurate treatment of electron correlations. Exploiting the advantages of contracting the primitive basis set over atomic hybrid orbitals, we are able to limit dramatically the number of variational parameters with a systematic control on the numerical precision, crucial in order to simulate larger systems. We investigate energetics and geometrical properties of the Zundel ion as a function of the oxygen-oxygen distance, taken as reaction coordinate. In both cases, our QMC results are found in excellent agreement with coupled cluster CCSD(T) technique, the quantum chemistry "gold standard". Calculations on proton transfer static barriers and dissociation energies display the same agreement. A comparison with density functional theory (DFT) results in the PBE approximation points out the crucial role of electron correlations for a correct description of the PT in the dimer. Our method is able to resolve the tiny energy differences (~ 0.1 Kcal/mol) and the corresponding structural variations at which the proton transfer takes place; it combines these features with a N^3-N^4 scaling with number of particles (favorable with respect to other post-DFT methods), as proven by a benchmark calculation on a larger protonated water cluster. The QMC approach used here represents a promising candidate to provide the first high-level ab initio description of PT in water.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Combined GW and dynamical mean field theory: Dynamical screening effects in transition metal oxides

We present the first dynamical implementation of the combined GW and dynamical mean field scheme ("GW+DMFT") for first principles calculations of the electronic properties of correlated materials. The application to the ternary transition metal oxide SrVO3 demonstrates that this schemes inherits the virtues of its two parent theories: a good description of the local low energy correlation physics encoded in a renormalized quasi-particle band structure, spectral weight transfer to Hubbard bands, and the physics of screening driven by long-range Coulomb interactions. Our data is in good agreement with available photoemission and inverse photoemission spectra; our analysis leads to a reinterpretation of the commonly accepted "three-peak structure" as originating from orbital effects rather than from the electron addition peak within the t2g manifold.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Intercalant and intermolecular phonon assisted superconductivity in K-doped picene

K$_3$-picene is a superconducting molecular crystal with critical temperature T$_c=7$K or 18K, depending on preparation conditions. Using density functional theory we show that electron-phonon interaction accounts for T$_c ~ 3-8$ K. The average electron-phonon coupling, calculated by including the phonon energy scale in the electron-phonon scattering, is $λ=0.73$ and $ω_{\rm log}=18.0\,$ meV. Intercalant and intermolecular phonon-modes contribute substantially (40%) to $λ$ as also shown by the isotope exponents of potassium (0.19) and carbon (0.31). The relevance of these modes makes superconductivity in K-doped picene peculiar and different from that of fullerenes.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Satellites and large doping- and temperature-dependence of electronic properties in hole-doped BaFe2As2

Over the last years, superconductivity has been discovered in several families of iron-based compounds. Despite intense research, even basic electronic properties of these materials, such as Fermi surfaces, effective electron masses, or orbital characters are still subject to debate. Here, we address an issue that has not been considered before, namely the consequences of dynamical screening of the Coulomb interactions among Fe-d electrons. We demonstrate its importance not only for correlation satellites seen in photoemission spectroscopy, but also for the low-energy electronic structure. From our analysis of the normal phase of BaFe2As2 emerges the picture of a strongly correlated compound with strongly doping- and temperature-dependent properties. In the hole overdoped regime, an incoherent metal is found, while Fermi-liquid behavior is recovered in the undoped compound. At optimal doping, the self-energy exhibits an unusual square-root energy dependence which leads to strong band renormalizations near the Fermi level.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Ab-initio calculations for the beta-tin diamond transition in Silicon: comparing theories with experiments

We investigate the pressure-induced metal-insulator transition from diamond to beta-tin in bulk Silicon, using quantum Monte Carlo (QMC) and density functional theory (DFT) approaches. We show that it is possible to efficiently describe many-body effects, using a variational wave function with an optimized Jastrow factor and a Slater determinant. Variational results are obtained with a small computational cost and are further improved by performing diffusion Monte Carlo calculations and an explicit optimization of molecular orbitals in the determinant. Finite temperature corrections and zero point motion effects are included by calculating phonon dispersions in both phases at the DFT level. Our results indicate that the theoretical QMC (DFT) transition pressure is significantly larger (smaller) than the accepted experimental value. We discuss the limitation of DFT approaches due to the choice of the exchange and correlation functionals and the difficulty to determine consistent pseudopotentials within the QMC framework, a limitation that may significantly affect the accuracy of the technique.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Hexatic and mesoscopic phases in the 2D quantum Coulomb system

We study the Wigner crystal melting in a two dimensional quantum system of particles interacting via the 1/r Coulomb potential. We use quantum Monte Carlo methods to calculate its phase diagram, locate the Wigner crystal region, and analyze its instabilities towards the liquid phase. We discuss the role of quantum effects in the critical behavior of the system, and compare our numerical results with the classical theory of melting, and the microemulsion theory of frustrated Coulomb systems. We find a Pomeranchuk effect much larger then in solid helium. In addition, we find that the exponent for the algebraic decay of the hexatic phase differs significantly from the Kosterilitz-Thouless theory of melting. We search for the existence of mesoscopic phases and find evidence of metastable bubbles but no mesoscopic phase that is stable in equilibrium.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Spin resolved energy parametrization of a quasi-one-dimensional electron gas

By carrying out extensive lattice regularized diffusion Monte Carlo calculations, we study the spin and density dependence of the ground state energy for a quasi-one-dimensional electron gas, with harmonic transverse confinement and long-range $1/r$ interactions. We present a parametrization of the exchange-correlation energy suitable for spin density functional calculations, which fulfills exact low and high density limits.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

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.

0 citations0 reviewsNo code linkedOpen access