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

Carlo Pierleoni

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

ResearcherUnclaimed source record
cond-mat.softcond-mat.mtrl-sciphysics.comp-phcond-mat.stat-mechcond-mat.str-elphysics.chem-phBiological Physicscond-mat.mes-hallcond-mat.otherquant-ph

Catalog footprint

What is connected

27works
10topics
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

27 published item(s)

preprint2021arXiv

Multi-scale simulation of the adsorption of lithium ion on graphite surface: from Quantum Monte Carlo to Molecular Density Functional Theory

The structure of the double-layer formed at the surface of carbon electrodes is governed by the interactions between the electrode and the electrolyte species. However, carbon is notoriously difficult to simulate accurately, even with well-established methods such as electronic Density Functional Theory and Molecular Dynamics. Here we focus on the important case of a lithium ion in contact with the surface of graphite, and we perform a series of reference Quantum Monte Carlo calculations that allow us to benchmark various electronic Density Functional Theory functionals. We then fit an accurate carbon--lithium pair potential, which is used in molecular Density Functional Theory calculations to determine the free energy of the adsorption of the ion on the surface in the presence of water. The adsorption profile in solution differs markedly from the gas phase results, which emphasize the role of the solvent on the properties of the double-layer.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Electronic structure and optical properties of quantum crystals from first principles calculations in the Born-Oppenheimer approximation

We develop a formalism to accurately account for the renormalization of electronic structure due to quantum and thermal nuclear motions within the Born-Oppenheimer approximation. We focus on the fundamental energy gap obtained from electronic addition and removal energies from Quantum Monte Carlo calculations in either the canonical or grand canonical ensembles. The formalism applies as well to effective single electron theories such as those based on Density Functional Theory. We show that electronic (Bloch) crystal momentum can be restored by marginalizing the total electron-ion wave function with respect to the nuclear equilibrium distribution, and we describe an explicit procedure to establish the band structure of electronic excitations for quantum crystals within the Born-Oppenheimer approximation. Based on the Kubo-Greenwood equation, we discuss the effects of nuclear motion on optical conductivity. Our methodology applies to the low temperature regime where nuclear motion is quantized and in general differs from the semi-classical approximation. We apply our method to study the electronic structure of C2/c-24 crystalline hydrogen at 200K and 250 GPa and discuss the optical absorption profile of hydrogen crystal at 200K and carbon diamond at 297K.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Energy gap closure of crystalline molecular hydrogen with pressure

We study the gap closure with pressure of crystalline molecular hydrogen. The gaps are obtained from grand-canonical Quantum Monte Carlo methods properly extended to quantum and thermal crystals, simulated by Coupled Electron Ion Monte Carlo. Nuclear zero point effects cause a large reduction in the gap ($\sim 2eV$). \CP{Depending on the structure,} the fundamental indirect gap closes \CP{between 380GPa and} 530GPa for ideal crystals and 330-380GPa for quantum crystals. Beyond this pressure the system enters into a bad metal phase where the density of states at the Fermi level increases with pressure up to $\sim$450\CP{-500} GPa when the direct gap closes. Our work partially supports the interpretation of recent experiments in high pressure hydrogen.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Electronic band gaps from Quantum Monte Carlo methods

We develop a method for calculating the fundamental electronic gap of semiconductors and insulators using grand canonical Quantum Monte Carlo simulations. We discuss the origin of the bias introduced by supercell calculations of finite size and show how to correct the leading and subleading finite size errors either based on observables accessible in the finite-sized simulations or from DFT calculations. Our procedure is applied to solid molecular hydrogen and compared to experiment for carbon and silicon crystals.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Electron localization properties in high pressure hydrogen at the liquid-liquid phase transition by Coupled Electron-Ion Monte Carlo

We analyze in detail the electronic properties of high pressure hydrogen around the liquid-liquid phase transition based on Coupled Electron-Ion Monte Carlo calculations. Computing the off-diagonal single particle density matrix and the momentum distribution we discuss localization properties of the electrons. The abrupt changes of these distributions indicate a metal to insulator transition occurring together with the structural transition from the atomic to molecular fluid. We further discuss the electron-proton and electron-electron pair correlation functions, which also change abruptly at the transition.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the Properties of a Bundle of Flexible Actin Filaments in an Optical Trap

We establish the Statistical Mechanics framework for a bundle of Nf living and uncrosslinked actin filaments in a supercritical solution of free monomers pressing against a mobile wall. The filaments are anchored normally to a fixed planar surface at one of their ends and, because of their limited flexibility, they grow almost parallel to each other. Their growing ends hit a moving obstacle, depicted as a second planar wall, parallel to the previous one and subjected to a harmonic compressive force. The force constant is denoted as trap strength while the distance between the two walls as trap length to make contact with the experimental optical trap apparatus. For an ideal solution of reactive filaments and free monomers at fixed free monomers chemical potential, we obtain the general expression for the grand potential from which we derive averages and distributions of relevant physical quantities, namely the obstacle position, the bundle polymerization force and the number of filaments in direct contact with the wall. The grafted living filaments are modeled as discrete Wormlike chains, with Factin persistence length, subject to discrete contour length variations to model single monomer (de)polymerization steps. Rigid filaments, either isolated or in bundles, all provide average values of the stalling force in agreement with Hill's predictions, independent of the average trap length. Flexible filaments instead, for values of the trap strength suitable to prevent their lateral escape, provide an average bundle force and an average trap length slightly larger than the corresponding rigid cases (few percents). Still the stalling force remains nearly independent on the average trap length, but results from the product of two strongly L dependent contributions: the fraction of touching filaments and the single filament buckling force.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Theory of Finite Size Effects for Electronic Quantum Monte Carlo Calculations of Liquids and Solids

Concentrating on zero temperature Quantum Monte Carlo calculations of electronic systems, we give a general description of the theory of finite size extrapolations of energies to the thermodynamic limit based on one and two-body correlation functions. We introduce new effective procedures, such as using the potential and wavefunction split-up into long and short range functions to simplify the method and we discuss how to treat backflow wavefunctions. Then we explicitly test the accuracy of our method to correct finite size errors on example hydrogen and helium many-body systems and show that the finite size bias can be drastically reduced for even small systems.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A semi-flexible model prediction for the polymerization force exerted by a living F-actin filament on a fixed wall

We consider a single living semi-flexible filament with persistence length l_p in chemical equilibrium with a solution of free monomers at fixed monomer chemical potential mu_1 and fixed temperature T. While one end of the filament is chemically active with single monomer (de)polymerization steps, the other end is grafted normally to a rigid wall to mimick a rigid network from which the filament under consideration emerges. A second rigid wall, parallel to the grafting wall, is fixed at distance L<<l_p from the filament seed. In supercritical conditions the filament tends to grow and impinges onto the second surface which, in suitable conditions (non-escaping filament regime) stops the filament growth. We first establish the grand-potential and derive some general properties, in particular the filament size distribution and the force exerted by the living filament on the obstacle wall. We apply this formalism to the semi-flexible, living, discrete Wormlike chain (d-WLC) model with step size d and persistence length l_p, hitting a hard wall. By original Monte-Carlo calculations we justify the use of the weak bending universal expressions of Gholami et al. (Phys.Rev.E. 74,(2006), 041803) over the whole non escaping filament regime. Employing this universal form for living filaments, we find that the average force exerted by a living filament on a wall at distance L is in practice L independent and very close to the value predicted by Hill, his expression being strictly valid in the rigid filament limit. The average filament force results from the product of the cumulative size fraction x, where the filament is in contact with the wall, times the buckling force on a filament of size L_c ~ L. We discuss several consequences of the L independence of the stalling force for our specific filament model.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Coarse-graining polymer solutions: a critical appraisal of single- and multi-site models

We critically discuss and review the general ideas behind single- and multi-site coarse-grained (CG) models as applied to macromolecular solutions in the dilute and semi-dilute regime. We first consider single-site models with zero-density and density-dependent pair potentials. We highlight advantages and limitations of each option in reproducing the thermodynamic behavior and the large-scale structure of the underlying reference model. As a case study we consider solutions of linear homopolymers in a solvent of variable quality. Secondly, we extend the discussion to multi-component systems presenting, as a test case, results for mixtures of colloids and polymers. Specifically, we found the CG model with zero-density potentials to be unable to predict fluid-fluid demixing in a reasonable range of densities for mixtures of colloids and polymers of equal size. For larger colloids, the polymer volume fractions at which phase separation occurs are largely overestimated. CG models with density-dependent potentials are somewhat less accurate than models with zero-density potentials in reproducing the thermodynamics of the system and, although they presents a phase separation, they significantly underestimate the polymer volume fractions along the binodal. Finally, we discuss a general multi-site strategy, which is thermodynamically consistent and fully transferable with the number of sites, and that allows us to overcome most of the limitations discussed for single-site models.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Integral-equation analysis of single-site coarse-grained models for polymer-colloid mixtures

We discuss the reliability of integral-equation methods based on several commonly used closure relations in determining the phase diagram of coarse-grained models of soft-matter systems characterized by mutually interacting soft and hard-core particles. Specifically, we consider a set of potentials appropriate to describe a system of hard-sphere colloids and linear homopolymers in good solvent, and investigate the behavior when the soft particles are smaller than the colloids, which is the regime of validity of the coarse-grained models. Using computer-simulation results as a benchmark, we find that the hypernetted-chain approximation provides accurate estimates of thermodynamics and structure in the colloid-gas phase in which the density of colloids is small. On the other hand, all closures considered appear to be unable to describe the behavior of the mixture in the colloid-liquid phase, as they cease to converge at polymer densities significantly smaller than those at the binodal. As a consequence, integral equations appear to be unable to predict a quantitatively correct phase diagram.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Phase diagram of mixtures of colloids and polymers in the thermal crossover from good to $θ$ solvent

We determine the phase diagram of mixtures of spherical colloids and neutral nonadsorbing polymers in the thermal crossover region between the $θ$ point and the good-solvent regime. We use the generalized free-volume theory (GFVT), which turns out to be quite accurate as long as $q = R_g/R_c\lesssim 1$ ($R_g$ is the radius of gyration of the polymer and $R_c$ is the colloid radius). Close to the $θ$ point the phase diagram is not very sensitive to solvent quality, while, close to the good-solvent region, changes of the solvent quality modify significantly the position of the critical point and of the binodals. We also analyze the phase behavior of aqueous solutions of charged colloids and polymers, using the extension of GFVT proposed by Fortini et al., J. Chem. Phys. 128, 024904 (2008).

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Accurate coarse-grained models for mixtures of colloids and linear polymers under good-solvent conditions

A coarse-graining strategy, previously developed for polymer solutions, is extended here to mixtures of linear polymers and hard-sphere colloids. In this approach groups of monomers are mapped onto a single pseudoatom (a blob) and the effective blob-blob interactions are obtained by requiring the model to reproduce some large-scale structural properties in the zero-density limit. We show that an accurate parametrization of the polymer-colloid interactions is obtained by simply introducing pair potentials between blobs and colloids. For the coarse-grained model in which polymers are modelled as four-blob chains (tetramers), the pair potentials are determined by means of the iterative Boltzmann inversion scheme, taking full-monomer pair correlation functions at zero-density as targets. For a larger number $n$ of blobs, pair potentials are determined by using a simple transferability assumption based on the polymer self-similarity. We validate the model by comparing its predictions with full-monomer results for the interfacial properties of polymer solutions in the presence of a single colloid and for thermodynamic and structural properties in the homogeneous phase at finite polymer and colloid density. The tetramer model is quite accurate for $q\lesssim 1$ ($q=\hat{R}_g/R_c$, where $\hat{R}_g$ is the zero-density polymer radius of gyration and $R_c$ is the colloid radius) and reasonably good also for $q=2$. For $q=2$ an accurate coarse-grained description is obtained by using the $n=10$ blob model. We also compare our results with those obtained by using single-blob models with state-dependent potentials.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Benchmark of Exchange-Correlation Functionals for High Pressure Hydrogen using Quantum Monte Carlo

The ab-initio phase diagram of dense hydrogen is very sensitive to errors in the treatment of electronic correlation. Recently, it has been shown that the choice of the density functional has a large effect on the predicted location of both the liquid-liquid phase transition and the solid insulator-to-metal transition in dense hydrogen. To identify the most accurate functional for dense hydrogen applications, we systematically benchmark some of the most commonly used functionals using Quantum Monte Carlo. By considering several measures of functional accuracy, we conclude that the van der Waals and hybrid functionals significantly out perform LDA and PBE. We support these conclusions by analyzing the impact of functional choice on structural optimization in the molecular solid, and on the location of the liquid-liquid phase transition.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Depletion effects in colloid-polymer solutions

The surface tension, the adsorption, and the depletion thickness of polymers close to a single nonadsorbing colloidal sphere are computed by means of Monte Carlo simulations. We consider polymers under good-solvent conditions and in the thermal crossover region between good-solvent and $θ$ behavior. In the dilute regime we consider a wide range of values of $q$, from $q = 0$ (planar surface) up to $q\approx 30$-50, while in the semidilute regime, for $ρ_p/ρ_p^*\le 4$ ($ρ_p$ is the polymer concentration and $ρ_p^*$ is its value at overlap), we only consider $q = 0,0.5,1$ and 2. The results are compared with the available theoretical predictions, verifying the existing scaling arguments. Field-theoretical results, both in the dilute and in the semidilute regime, are in good agreement with the numerical estimates for polymers under good-solvent conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Molecular-Atomic Transition in the Deuterium Hugoniot with Coupled Electron Ion Monte Carlo

We have performed accurate simulations of the Deuterium Hugoniot using Coupled Electron Ion Monte Carlo (CEIMC). Using highly accurate quantum Monte Carlo methods for the electrons, we study the region of maximum compression along the principal Hugoniot, where the system undergoes a continuous transition from a molecular fluid to a monatomic fluid. We include all relevant physical corrections so that a direct comparison to experiment can be made. Around 50 GPa we found a maximum compression of 4.85, roughly 10% larger than previous theoretical predictions and experimental data but still compatible with the latter because of their large uncertainty.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Consistent coarse-graining strategy for polymer solutions in the thermal crossover from Good to Theta solvent

We extend our previously developed coarse-graining strategy for linear polymers with a tunable number n of effective atoms (blobs) per chain [D'Adamo et al., J. Chem. Phys. 137, 4901 (2012)] to polymer systems in thermal crossover between the good-solvent and the Theta regimes. We consider the thermal crossover in the region in which tricritical effects can be neglected, i.e. not too close to the Theta point, for a wide range of chain volume fractions Phi=c/c* (c* is the overlap concentration), up to Phi=30. Scaling crossover functions for global properties of the solution are obtained by Monte-Carlo simulations of the Domb-Joyce model. They provide the input data to develop a minimal coarse-grained model with four blobs per chain. As in the good-solvent case, the coarse-grained model potentials are derived at zero density, thus avoiding the inconsistencies related to the use of state-dependent potentials. We find that the coarse-grained model reproduces the properties of the underlying system up to some reduced density which increases when lowering the temperature towards the Theta state. Close to the lower-temperature crossover boundary, the tetramer model is accurate at least up to Phi<10, while near the good-solvent regime reasonably accurate results are obtained up to Phi<2. The density region in which the coarse-grained model is predictive can be enlarged by developing coarse-grained models with more blobs per chain. We extend the strategy used in the good-solvent case to the crossover regime. This requires a proper treatment of the length rescalings as before, but also a proper temperature redefinition as the number of blobs is increased. The case n=10 is investigated. Comparison with full-monomer results shows that the density region in which accurate predictions can be obtained is significantly wider than that corresponding to the n=4 case.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Predicting the thermodynamics by using state-dependent interactions

We reconsider the structure-based route to coarse graining in which the coarse-grained model is defined in such a way to reproduce some distributions functions of the original system as accurately as possible. We consider standard expressions for pressure and chemical potential applied to this family of coarse-grained models with density-dependent interactions and show that they only provide approximations to the pressure and chemical potential of the underlying original system. These approximations are then carefully compared in two cases: we consider a generic microscopic system in the low-density regime and polymer solutions under good-solvent conditions. Moreover, we show that the state-dependent potentials depend on the ensemble in which they have been derived. Therefore, care must be used in applying canonical state-dependent potentials to predict phase lines, which is typically performed in other ensembles.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Coarse-graining strategies in polymer solutions

We review a coarse-graining strategy (multiblob approach) for polymer solutions in which groups of monomers are mapped onto a single atom (a blob) and effective blob-blob interactions are obtained by requiring the coarse-grained model to reproduce some coarse-grained features of the zero-density isolated-chain structure. By tuning the level of coarse graining, i.e. the number of monomers to be mapped onto a single blob, the model should be adequate to explore the semidilute regime above the collapse transition, since in this case the monomer density is very small if chains are long enough. The implementation of these ideas has been previously based on a transferability hypothesis, which was not completely tested against full-monomer results (Pierleoni et al., J. Chem. Phys, 127, 171102 (2007)). We study different models proposed in the past and we compare their predictions to full-monomer results for the chain structure and the thermodynamics in the range of polymer volume fractions Φbetween 0 and 8. We find that the transferability assumption has a limited predictive power if a thermodynamically consistent model is required. We introduce a new tetramer model parametrized in such a way to reproduce not only zero-density intramolecular and intermolecular two-body probabilities, but also some intramolecular three-body and four-body distributions. We find that such a model correctly predicts three-chain effects, the structure and the thermodynamics up to Φ~ 2, a range considerably larger than that obtained with previous simpler models using zero-density potentials. Our results show the correctness of the ideas behind the multiblob approach but also that more work is needed to understand how to develop models with more effective monomers which would allow us to explore the semidilute regime at larger chain volume fractions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Consistent and transferrable coarse-grained model for semidilute polymer solutions in good solvent

We present a coarse-grained model for linear polymers with a tunable number of effective atoms (blobs) per chain interacting by intra- and inter-molecular potentials obtained at zero density. We show how this model is able to accurately reproduce the universal properties of the underlying solution of athermal linear chains at various levels of coarse-graining and in a range of chain densities which can be widened by increasing the spatial resolution of the multiblob representation, i.e., the number of blobs per chain. The present model is unique in its ability to quantitatively predict thermodynamic and large scale structural properties of polymer solutions deep in the semidilute regime with a very limited computational effort, overcoming most of the problems related to the simulations of semidilute polymer solutions in good solvent conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Crystalline free energies of micelles of diblock copolymer solutions

We report a characterization of the relative stability and structural behavior of various micellar crystals of an athermal model of AB-diblock copolymers in solution. We adopt a previously devel- oped coarse-graining representation of the chains which maps each copolymer on a soft dumbbell. Thanks to this strong reduction of degrees of freedom, we are able to investigate large aggregated systems, and for a specific length ratio of the blocks f = MA/(MA + MB) = 0.6, to locate the order-disorder transition of the system of micelles. Above the transition, mechanical and thermal properties are found to depend on the number of particles per lattice site in the simulation box, and the application of a recent methodology for multiple occupancy crystals (B.M. Mladek et al., Phys. Rev. Lett. 99, 235702 (2007)) is necessary to correctly define the equilibrium state. Within this scheme we have performed free energy calculations at two reduced density ρ/ρ\ast = 4,5 and for several cubic structures as FCC,BCC,A15. At both densities, the BCC symmetry is found to correspond to the minimum of the unconstrained free energy, that is to the stable symmetry among the few considered, while the A15 structure is almost degenerate, indicating that the present sys- tem prefers to crystallize in less packed structures. At ρ/ρ\ast = 4 close to melting, the Lindemann ratio is fairly high (~ 0.29) and the concentration of vacancies is roughly 6%. At ρ/ρ\ast = 5 the mechanical stability of the stable BCC structure increases and the concentration of vacancies ac- cordingly decreases. The ratio of the corona layer thickness to the core radius is found to be in good agreement with experimental data for poly(styrene-b-isoprene)(22-12) in isoprene selective solvent which is also reported to crystallize in the BCC structure.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Polymers as compressible soft spheres

We consider a coarse-grained model in which polymers under good-solvent conditions are represented by soft spheres whose radii, which should be identified with the polymer radii of gyrations, are allowed to fluctuate. The corresponding pair potential depends on the sphere radii. This model is a single-sphere version of the one proposed in Vettorel et al., Soft Matter 6, 2282 (2010), and it is sufficiently simple to allow us to determine all potentials accurately from full-monomer simulations of two isolated polymers (zero-density potentials). We find that in the dilute regime (which is the expected validity range of single-sphere coarse-grained models based on zero-density potentials) this model correctly reproduces the density dependence of the radius of gyration. However, for the thermodynamics and the intermolecular structure, the model is largely equivalent to the simpler one in which the sphere radii are fixed to the average value of the radius of gyration and radiiindependent potentials are used: for the thermodynamics there is no advantage in considering a fluctuating sphere size.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Free energy methods in Coupled Electron Ion Monte Carlo

Recent progress in simulation methodologies and in computer power allow first principle simulations of condensed systems with Born-Oppenheimer electronic energies obtained by Quantum Monte Carlo methods. Computing free energies and therefore getting a quantitative determination of phase diagrams is one step more demanding in terms of computer resources. In this paper we derive a general relation to compute the free energy of an ab-initio model with Reptation Quantum Monte Carlo (RQMC) energies from the knowledge of the free energy of the same ab-initio model in which the electronic energies are computed by the less demanding but less accurate Variational Monte Carlo (VMC) method. Moreover we devise a procedure to correct transition lines based on the use of the new relation. In order to illustrate the procedure, we consider the liquid-liquid phase transition in hydrogen, a first order transition between a lower pressure, molecular and insulating phase and a higher pressure, partially dissociated and conducting phase. We provide new results along the T = 600K isotherm across the phase transition and find good agreement between the transition pressure and specific volumes at coexistence for the model with RQMC accuracy between the prediction of our procedure and the values that can be directly inferred from the observed plateau in the pressure-volume curve along the isotherm. This work paves the way for future use of VMC in first principle simulations of high pressure hydrogen, an essential simplification when considering larger system sizes or quantum proton effects by Path Integral Monte Carlo methods.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The momentum distribution of the homogeneous electron gas

We calculate the off-diagonal density matrix of the homogeneous electron gas at zero temperature using unbiased Reptation Monte Carlo for various densities and extrapolate the momentum distribution, and the kinetic and potential energies to the thermodynamic limit. Our results on the renormalization factor allows us to validate approximate G_0W_0 calculations concerning quasiparticle properties over a broad density region (1 <= r_s <= 10) and show that near the Fermi surface, vertex corrections and self-consistency aspects almost cancel each other out.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Electrical conductivity of high-pressure liquid hydrogen by quantum Monte Carlo methods

We compute the electrical conductivity for liquid hydrogen at high pressure using quantum Monte Carlo. The method uses Coupled Electron-Ion Monte Carlo to generate configurations of liquid hydrogen. For each configuration correlated sampling of electrons is performed in order to calculate a set of lowest many-body eigenstates and current-current correlation functions of the system, which are summed over in the many-body Kubo formula to give AC electrical conductivity directly. The extrapolated DC conductivity at 3000 K for several densities shows a liquid semiconductor to liquid-metal transition at high pressure. Our results are in good agreement with shock-wave data.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Equation of state of metallic hydrogen from Coupled Electron-Ion Monte Carlo simulations

We present a study of hydrogen at pressures higher than molecular dissociation using the Coupled Electron-Ion Monte Carlo method. These calculations use the accurate Reptation Quantum Monte Carlo method to estimate the electronic energy and pressure while doing a Monte Carlo simulation of the protons. In addition to presenting simulation results for the equation of state over a large region of phase space, we report the free energy obtained by thermodynamic integration. We find very good agreement with DFT calculations for pressures beyond 600 GPa and densities above $ρ=1.4 g/cm^3$. Both thermodynamic as well as structural properties are accurately reproduced by DFT calculations. This agreement gives a strong support to the different approximations employed in DFT, specifically the approximate exchange-correlation potential and the use of pseudopotentials for the range of densities considered. We find disagreement with chemical models, which suggests a reinvestigation of planetary models, previously constructed using the Saumon-Chabrier-Van Horn equations of state.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

The induced charge in a Frohlich polaron: Sum rule and spatial extent

Within the path integral formalism, we derive exact expressions for correlation functions measuring the lattice charge induced by an electron and associated polarization in Frohlich polaron problem. We prove that a sum rule for the total induced charge, already obtained within approximated approaches is indeed exact. As a consequence the total induced charge is shown rigorously to be temperature independent. In addition we perform path integral Monte Carlo calculations of the correlation functions and we compare with variational results based on Feynman method. As the temperature increases the polaron radius decreases. On the other hand at high temperatures the electron motion is not hindered by the lattice. These apparently contradictory results are discussed.

0 citations0 reviewsNo code linkedOpen access