BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Henry Krakauer

Henry Krakauer contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.mtrl-scicond-mat.str-elphysics.comp-phphysics.chem-phcond-mat.mes-hallcond-mat.othercond-mat.softcond-mat.supr-conmtrl-thphysics.geo-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
17works
0followers
10topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

17 published item(s)

preprint2016arXiv

Auxiliary-field quantum Monte Carlo calculations of the molybdenum dimer

Chemical accuracy is difficult to achieve for systems with transition metal atoms. Third row transition metal atoms are particularly challenging due to strong electron-electron correlation in localized d-orbitals. The Cr2 molecule is an outstanding example, which we previously treated with highly accurate auxiliary-field quantum Monte Carlo (AFQMC) calculations [Purwanto et al., J. Chem. Phys. 142, 064302 (2015)]. Somewhat surprisingly, computational description of the isoelectronic Mo2 dimer has also, to date, been scattered and less than satisfactory. We present high-level theoretical benchmarks of the Mo2 singlet ground state ($X ^1Σ_g^+$) and first triplet excited state ($a ^3Σ_u^+$), using the phaseless AFQMC calculations. Extrapolation to the complete basis set limit is performed. Excellent agreement with experimental spectroscopic constants is obtained. We also present a comparison of the correlation effects in Cr2 and Mo2.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Electron-phonon coupling and exchange-correlation effects in superconducting ${\rm H_3S}$ under high pressure

We investigate the ${\rm H_3S}$ phase of sulphur hydride under high pressure $\simeq$ 200 GPa by means of {\it ab-initio} calculations within the framework of the density-functional theory (DFT) with the PBE0 hybrid exchange-correlation ($E_{\rm xc}$) approximation. The choice of $E_{\rm xc}$ has the largest effect on the calculated electron-phonon coupling (EPC) matrix elements; the high pressure equation of state and phonon frequencies are only slightly modified. Mode-dependent EPC correction factors are determined from PBE0 using a frozen-phonon supercell approach, while standard density-functional perturbation theory is used to determine the EPC with PBE generalized-gradient approximation $E_{\rm xc}$. Our principle finding is that the calculated PBE0 $T_c$ is enhanced by 25\% compared to PBE. This is similar in magnitude, but in opposite direction, to the proposed suppression of $T_c$ by anharmonic effects [Errea {\em et al., Phys. Rev. Lett.} {\bf 114}, 157004 (2015)]. Our calculations demonstrate the importance of considering exchange-correlation approximations for calculations of superconducting properties for this class of materials.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

An Auxiliary-Field Quantum Monte Carlo Study of the Chromium Dimer

The chromium dimer (Cr2) presents an outstanding challenge for many-body electronic structure methods. Its complicated nature of binding, with a formal sextuple bond and an unusual potential energy curve, is emblematic of the competing tendencies and delicate balance found in many strongly correlated materials. We present a near-exact calculation of the potential energy curve (PEC) and ground state properties of Cr2, using the auxiliary-field quantum Monte Carlo (AFQMC) method. Unconstrained, exact AFQMC calculations are first carried out for a medium-sized but realistic basis set. Elimination of the remaining finite-basis errors and extrapolation to the complete basis set (CBS) limit is then achieved with a combination of phaseless and exact AFQMC calculations. Final results for the PEC and spectroscopic constants are in excellent agreement with experiment.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Quantum Monte Carlo calculations in solids with downfolded Hamiltonians

We present a systematic downfolding many-body approach for extended systems. Many-body calculations operate on a simpler Hamiltonian which retains material-specific properties. The Hamiltonian is systematically improvable and allows one to dial, in principle, between the simplest model and the original Hamiltonian. As a by-product, pseudopotential errors are essentially eliminated using a frozen-core treatment. The computational cost of the many-body calculation is dramatically reduced without sacrificing accuracy. We use the auxiliary-field quantum Monte Carlo (AFQMC) method to solve the downfolded Hamiltonian. Excellent accuracy is achieved for a range of solids, including semiconductors, ionic insulators, and metals. We further test the method by determining the spin gap in NiO, a challenging prototypical material with strong electron correlation effects. This approach greatly extends the reach of general, ab initio many-body calculations in materials.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Stability, Energetics, and Magnetic States of Cobalt Adatoms on Graphene

We investigate the stability and electronic properties of single Co atoms on graphene with near-exact many-body calculations. A frozen-orbital embedding scheme was combined with auxiliary-field quantum Monte Carlo to increase the reach in system sizes. Several energy minima are found as a function of the distance $h$ between Co and graphene. Energetics only permit the Co atom to occupy the top site at $h = 2.2$ Å in a high-spin $3d^{8}4s^{1}$ state, and the van der Waals region at $h = 3.3$ Å in a high-spin $3d^{7}4s^{2}$ state. The findings provide an explanation for recent experimental results with Co on free-standing graphene.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Frozen-orbital and downfolding calculations with auxiliary-field quantum Monte Carlo

We describe the implementation of the frozen-orbital and downfolding approximations in the auxiliary-field quantum Monte Carlo (AFQMC) method. These approaches can provide significant computational savings compared to fully correlating all the electrons. While the many-body wave function is never explicit in AFQMC, its random walkers are Slater determinants, whose orbitals may be expressed in terms of any one-particle orbital basis. It is therefore straightforward to partition the full N-particle Hilbert space into active and inactive parts to implement the frozen-orbital method. In the frozen-core approximation, for example, the core electrons can be eliminated in the correlated part of the calculations, greatly increasing the computational efficiency, especially for heavy atoms. Scalar relativistic effects are easily included using the Douglas-Kroll-Hess theory. Using this method, we obtain a way to effectively eliminate the error due to single-projector, norm-conserving pseudopotentials in AFQMC. We also illustrate a generalization of the frozen-orbital approach that downfolds high-energy basis states to a physically relevant low-energy sector, which allows a systematic approach to produce realistic model Hamiltonians to further increase efficiency for extended systems.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Ab Initio Many-body Study of Cobalt Adatoms Adsorbed on Graphene

Many recent calculations have been performed to study a Co atom adsorbed on graphene, with significantly varying results on the nature of the bonding. We use auxiliary-field quantum Monte Carlo (AFQMC) and a size-correction embedding scheme to accurately calculate the binding energy of Co on graphene. We find that as a function of the distance h between the Co atom and the six-fold hollow site, there are three distinct ground states corresponding to three electronic configurations of the Co atom. Two of these states provide binding and exhibit a double-well feature with nearly equal binding energy of 0.4 eV at h = 1.51 and h = 1.65 angstroms, corresponding to low-spin $^2$Co (3d9 4s0) and high-spin $^4$Co (3d8 4s1), respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Excited state calculations in solids by auxiliary-field quantum Monte Carlo

We present an approach for ab initio many-body calculations of excited states in solids. Using auxiliary-field quantum Monte Carlo, we introduce an orthogonalization constraint to prevent collapse of the stochastic Slater determinants in the imaginary-time propagation. Trial wave functions from density-functional calculations are used for the constraints. Detailed band structures calculated for standard semiconductors are in good agreement with GW and experimental results. For the challenging ZnO wurtzite structure, we obtain a fundamental band gap of 3.26(16) eV, consistent with experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Assessing weak hydrogen binding on Ca+ centers: An accurate many-body study with large basis sets

Weak H2 physisorption energies present a significant challenge to even the best correlated theoretical many-body methods. We use the phaseless auxiliary-field quantum Monte Carlo (AFQMC) method to accurately predict the binding energy of Ca+ - 4H2. Attention has recently focused on this model chemistry to test the reliability of electronic structure methods for H2 binding on dispersed alkaline earth metal centers. A modified Cholesky decomposition is implemented to realize the Hubbard-Stratonovich transformation efficiently with large Gaussian basis sets. We employ the largest correlation-consistent Gaussian type basis sets available, up to cc-pCV5Z for Ca, to accurately extrapolate to the complete basis limit. The calculated potential energy curve exhibits binding with a double-well structure.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Finite-size correction in many-body electronic structure calculations of magnetic systems

We extend the post-processing finite-size (FS) correction method, developed by Kwee, Zhang, and Krakauer [Phys. Rev. Lett. 100, 126404 (2008)], to spin polarized systems. The method estimates the FS effects in many-body electronic structure calculations of extended systems by a modified density functional theory (DFT) calculation, without having to repeat expensive many-body simulations. We construct a unified FS DFT exchange-correlation functional for spin unpolarized and fully spin polarized systems, under the local density approximation. The results are then interpolated to arbitrary spin polarizations. Generalization to other functional forms in DFT are discussed. The application of this FS correction method to several typical magnetic systems with varying supercell sizes demonstrates that it consistently removes most of the FS errors, leading to rapid convergence of the many-body results to the infinite size limit.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

First-principles calculations of $^{17}$O NMR chemical shielding in Pb(Zr$_{1/2}$Ti$_{1/2}$)O$_3$ and Pb(Mg$_{1/3}$Nb$_{2/3}$)O$_3$: linear dependence on transition-metal/oxygen bond lengths

First-principles density functional theory (DFT) oxygen chemical shift tensors were calculated for A(B,B$^{\prime}$)O$_3$ perovskite alloys Pb(Zr$_{1/2}$Ti$_{1/2}$)O$_3$ (PZT) and Pb(Mg$_{1/3}$Nb$_{2/3}$)O$_3$ (PMN). Quantum chemistry methods for embedded clusters and the GIPAW method [C.\ J.\ Pickard and F.\ Mauri, {\it Phys. Rev. B} {\bf 63} 245101 (2001)] for periodic boundary conditions were used. Results from both methods are in good agreement for PZT and prototypical perovskites. PMN results were obtained using only GIPAW. Both isotropic $δ_\mathrm{iso}$ and axial $δ_\mathrm{ax}$ chemical shifts were found to vary approximately linearly as a function of the nearest-distance transition-metal/oxygen bond length,$r_\mathrm{s}$. Using these results, we argue against Ti clustering in PZT, as conjectured from recent $^{17}$O NMR magic-angle-spinning measurements. Our findings indicate that $^{17}$O NMR measurements, coupled with first-principles calculations, can be an important probe of local structure in complex perovskite solid solutions.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

High sensitivity of 17O NMR to p-d hybridization in transition metal perovskites: first principles calculations of large anisotropic chemical shielding

A first principles embedded cluster approach is used to calculate O chemical shielding tensors, sigma, in prototypical transition metal oxide ABO_3 perovskite crystals. Our principal findings are 1) a large anisotropy of sigma between deshielded sigma_x ~ sigma_y and shielded sigma_z components (z along the Ti-O bond); 2) a nearly linear variation, across all the systems studied, of the isotropic sigma_iso and uniaxial sigma_ax components, as a function of the B-O-B bond asymmetry. We show that the anisotropy and linear variation arise from large paramagnetic contributions to sigma_x and sigma_y due to virtual transitions between O(2p) and unoccupied B(nd) states. The calculated isotropic delta_iso and uniaxial delta_ax chemical shifts are in good agreement with recent BaTiO_3 and SrTiO_3 single crystal 17O NMR measurements. In PbTiO_3 and PbZrO_3, calculated delta_iso are also in good agreement with NMR powder spectrum measurements. In PbZrO_3, delta_iso calculations of the five chemically distinct sites indicate a correction of the experimental assignments. The strong dependence of sigma on covalent O(2p)-B(nd) interactions seen in our calculations indicates that 17O NMR spectroscopy, coupled with first principles calculations, can be an especially useful tool to study the local structure in complex perovskite alloys.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Excited state calculations using phaseless auxiliary-field quantum Monte Carlo: potential energy curves of low lying C2 singlet states

We show that the recently developed phaseless auxiliary-field quantum Monte Carlo (AFQMC) method can be used to study excited states, providing an alternative to standard quantum chemistry methods. The phaseless AFQMC approach, whose computational cost scales as M^3-M^4 with system size M, has been shown to be among the most accurate many-body methods in ground state calculations. For excited states, prevention of collapse into the ground state and control of the Fermion sign/phase problem are accomplished by the approximate phaseless constraint with a trial wave function. Using the challenging C2 molecule as a test case, we calculate the potential energy curves of the ground and two low-lying singlet excited states. The trial wave function is obtained by truncating complete active space wave functions, with no further optimization. The phaseless AFQMC results using a small basis set are in good agreement with exact full configuration interaction calculations, while those using large basis sets are in good agreement with experimental spectroscopic constants.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Kinetic Monte Carlo Simulations of Crystal Growth in Ferroelectric Alloys

The growth rates and chemical ordering of ferroelectric alloys are studied with kinetic Monte Carlo (KMC) simulations using an electrostatic model with long-range Coulomb interactions, as a function of temperature, chemical composition, and substrate orientation. Crystal growth is characterized by thermodynamic processes involving adsorption and evaporation, with solid-on-solid restrictions and excluding diffusion. A KMC algorithm is formulated to simulate this model efficiently in the presence of long-range interactions. Simulations were carried out on Ba(Mg_{1/3}Nb_{2/3})O_3 (BMN) type materials. Compared to the simple rocksalt ordered structures, ordered BMN grows only at very low temperatures and only under finely tuned conditions. For materials with tetravalent compositions, such as (1-x)Ba(Mg_{1/3}Nb_{2/3})O_3 + xBaZrO_3 (BMN-BZ), the model does not incorporate tetravalent ions at low-temperature, exhibiting a phase-separated ground state instead. At higher temperatures, tetravalent ions can be incorporated, but the resulting crystals show no chemical ordering in the absence of diffusive mechanisms.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Polarization Dependence of Born Effective Charge and Dielectric Constant in KNbO$_3$

The Born effective charge Z^{*} and dielectric tensor ε_{\infty} of KNbO_3 are found to be very sensitive to the atomic geometry, changing by as much as 27% between the paraelectric cubic and ferroelectric tetragonal and rhombohedral phases. Subtracting the bare ionic contribution reveals changes of the dynamic component of Z^{*} as large as 50%, for atomic displacements that are typically only a few percent of the lattice constant. Z^{*}, ε_{\infty} and all phonon frequencies at the Brillouin zone center were calculated using the {\it ab initio} linearized augmented plane-wave linear response method with respect to the reference cubic, experimental tetragonal, and theoretically determined rhombohedral ground state structures. The ground state rhombohedral structure of KNbO_3 was determined by minimizing the forces on the relaxed atoms. By contrast with the cubic structure, all zone center phonon modes of the rhombohedral structure are stable and their frequencies are in good agreement with experiment. In the tetragonal phase, one of the soft zone center modes in the cubic phase is stablized. In view of the small atomic displacements involved in the ferroelectric transitions, it is evident that not only the soft mode frequencies but also the Born effective charge and dielectric constants are very sensitive to the atomic geometry.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Prediction of Phase Transition in CaSiO$_3$ Perovskite and Implications for Lower Mantle Structure

First principles linear response calculations are used to investigate the lattice dynamics of what is thought to be the third most abundant phase in the lower mantle, CaSiO_3 perovskite. The commonly assumed cubic structure (Pm3m) is found to be dynamically unstable at all pressures, exhibiting unstable modes along the Brillouin zone edge from the M-point to the R-point. Based on these results, we predict that the ground state structure of CaSiO_3 perovskite is a distorted phase with lower than cubic symmetry. We predict that a phase transition occurs in CaSiO_3 perovskite within the earth's lower mantle from the low temperature distorted phase to the cubic phase at high temperature. The predicted phase transition provides a possible explanation of some of the seismological observations of reflective features within the lower mantle.

0 citations0 reviewsNo code linkedOpen access
preprint1994arXiv

First-Principles Determination of Chain-Structure Instability in KNbO$_3$

A complete mapping in the Brillouin zone of the structural instability associated with the ferroelectric phase transitions of KNbO$_3$ has been obtained by first-principles calculations using an LAPW linear response approach. The wavevector dependence of the instability reveals pronounced two-dimensional character, which corresponds to chains oriented along $\left<100\right>$\ directions of displaced Nb atoms. The results are discussed in relation to models of the ferroelectric phase transitions.

0 citations0 reviewsNo code linkedOpen access