Non-perturbative many-body treatment of molecular magnets.

Molecular magnets have received significant attention because of their potential applications in quantum information and quantum computing. A delicate balance of electron correlation, spin-orbit coupling (SOC), ligand field splitting, and other effects produces a persistent magnetic moment within each molecular magnet unit. The discovery and design of molecular magnets with improved functionalities would be greatly aided by accurate computations. However, the competition among the different effects poses a challenge for theoretical treatments. Electron correlation plays a central role since d- or f-element ions, which provide the magnetic states in molecular magnets, often require explicit many-body treatments. SOC, which expands the dimensionality of the Hilbert space, can also lead to non-perturbative effects in the presence of strong interaction. Furthermore, molecular magnets are large, with tens of atoms in even the smallest systems. We show how an ab initio treatment of molecular magnets can be achieved with auxiliary-field quantum Monte Carlo, in which electron correlation, SOC, and material specificity are included accurately and on an equal footing. The approach is demonstrated by an application to compute the zero-field splitting of a locally linear Co2+ complex.

Ab initio calculations in atoms, molecules, and solids, treating spin-orbit coupling and electron interaction on an equal footing.

We incorporate explicit, non-perturbative treatment of spin-orbit coupling into ab initio auxiliary-field quantum Monte Carlo (AFQMC) calculations. The approach allows a general computational framework for molecular and bulk systems in which material specificity, electron correlation, and spin-orbit coupling effects can be captured accurately and on an equal footing, with favorable computational scaling vs system size. We adopt relativistic effective-core potentials that have been obtained by fitting to fully relativistic data and that have demonstrated a high degree of reliability and transferability in molecular systems. This results in a two-component spin-coupled Hamiltonian, which is then treated by generalizing the ab initio AFQMC approach. We demonstrate the method by computing the electron affinity in Pb, the bond dissociation energy in Br2 and I2, and solid Bi.

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 [W. 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.

Quantum Monte Carlo Calculations in Solids with Downfolded Hamiltonians.

We present a combination of a downfolding many-body approach with auxiliary-field quantum Monte Carlo (AFQMC) calculations 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 frozen orbitals constructed adaptively from the solid environment. The computational cost of the many-body calculation is dramatically reduced without sacrificing accuracy. Excellent accuracy is achieved for a range of solids, including semiconductors, ionic insulators, and metals. We apply the method to calculate the equation of state of cubic BN under ultrahigh pressure, and determine the spin gap in NiO, a challenging prototypical material with strong electron correlation effects.

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 (PEC), is emblematic of the competing tendencies and delicate balance found in many strongly correlated materials. We present an accurate calculation of the 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 limit are 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.

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 calculations 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.

First-principles calculations of 17O nuclear magnetic resonance chemical shielding in Pb(Zr(1/2)Ti(1/2))O3 and Pb(Mg(1/3)Nb(2/3))O3: linear dependence on transition-metal/oxygen bond lengths.

First-principles density functional theory oxygen chemical shift tensors were calculated for A(B,B')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 gauge including projector augmented waves (GIPAW) method [C. J. Pickard and F. Mauri, Phys. Rev. B 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 δ(iso) and axial δ(ax) chemical shifts were found to vary approximately linearly as a function of the nearest-distance transition-metal/oxygen bond length, r(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.

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

Weak H(2) physisorption energies present a significant challenge to even the best correlated theoretical many-body methods. We use the phaseless auxiliary-field quantum Monte Carlo method to accurately predict the binding energy of Ca(+)-4H(2). Attention has recently focused on this model chemistry to test the reliability of electronic structure methods for H(2) 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.

Excited state calculations using phaseless auxiliary-field quantum Monte Carlo: Potential energy curves of low-lying C(2) 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 C(2) 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.

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) approximately 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 (17)O 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 (17)O NMR spectroscopy, coupled with first principles calculations, can be an especially useful tool to study the local structure in complex perovskite alloys.

Local Embedding and Effective Downfolding in the Auxiliary-Field Quantum Monte Carlo Method.

A local embedding and effective downfolding scheme has been developed and implemented in the auxiliary-field quantum Monte Carlo (AFQMC) method. A local cluster in which electrons are fully correlated is defined, and the frozen orbital method is used on the remainder of the system to construct an effective Hamiltonian, which operates within the local cluster. Local embedding, which involves only the occupied sector, has previously been employed in the context of Co/graphene. Here, the methodology is extended in order to allow for effective downfolding of the virtual sector, thus allowing for significant reduction in the computational effort required for AFQMC calculations. The system size, which can be feasibly treated with AFQMC, is therefore greatly extended as only a single local cluster is explicitly correlated at the AFQMC level of theory. The approximation is controlled by the separate choice of the spatial size of the active occupied region ( Ro) and of the active virtual region ( Rv). The systematic dependence of the AFQMC energy on Ro and Rv is investigated, and it is found that relative AFQMC energies of physical and chemical interest converge rapidly to the full AFQMC treatment (i.e., using no embedding or downfolding).

Projected survival benefit as criterion for listing and organ allocation in heart transplantation.

BACKGROUND: Current policies for the selection of candidates and the allocation of hearts for transplantation give priority to patients at greatest risk if not transplanted. However, to achieve best use of the donated organs, it is necessary to estimate the net benefit associated with transplantation. METHODS: The survival benefit associated with being listed or not, with being transplanted or left on the waiting list, or with being transplanted or being denied the opportunity for a transplant can be estimated by means of time-to-event modeling of competing risks with intervening states. The data were obtained from the Organ Procurement and Transplantation Network and describe the outcomes of listings from 1997 to June 1999. Our analyses assessed 9,059 heart transplantation candidates, who were followed for at least 1 and up to 2 years after listing. RESULTS: The probability of receiving a heart transplant does not increase with the probability of death while awaiting the transplant. It is comparable in the second and tenth deciles of risk (measured by the probability of death while awaiting a transplant) at 1 month after listing (15% vs 18%), but is considerably higher in the second decile at 6 months (53% vs 38%) and increasingly more so thereafter. The estimate of survival benefit stabilizes within 1 year of follow-up. Through the fourth decile of risk, the benefit of being placed on the waiting list is negligible at best, but becomes substantial (10%) for patients in the highest 2 deciles. Heart transplantation may reduce survival in the least ill patients but is clearly strongly beneficial for severely ill patients, offering reductions of 20 to 35 percentage points in probability of death when compared with remaining on the waiting list or not receiving a transplant at all. CONCLUSIONS: Our analyses indicate that criteria other than severity of illness as measured by the probability of death are, in practice, dominant in the allocation of donated hearts for transplantation. High percentages of patients listed as well as those transplanted are not expected to undergo a substantial increase in probability of survival, and some are likely to be harmed. A survival benefit is anticipated only for severely ill patients. Estimation of the projected survival benefit of listing and of transplantation is feasible and may be used to prioritize patients and lead to the best use of donated organs.

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 of 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.

Finite-size correction in many-body electronic structure calculations.

Finite-size (FS) effects are a major source of error in many-body (MB) electronic structure calculations of extended systems. A method is presented to correct for such errors. We show that MB FS effects can be effectively included in a modified local density approximation calculation. A parametrization for the FS exchange-correlation functional is obtained. The method is simple and gives post-processing corrections that can be applied to any MB results. Applications to a model insulator (P2 in a supercell), to semiconducting Si, and to metallic Na show that the method delivers greatly improved FS corrections.

Eliminating spin contamination in auxiliary-field quantum Monte Carlo: realistic potential energy curve of F(2).

The use of an approximate reference state wave function mid R:Phi(r) in electronic many-body methods can break the spin symmetry of Born-Oppenheimer spin-independent Hamiltonians. This can result in significant errors, especially when bonds are stretched or broken. A simple spin-projection method is introduced for auxiliary-field quantum Monte Carlo (AFQMC) calculations, which yields spin-contamination-free results, even with a spin-contaminated mid R:Phi(r). The method is applied to the difficult F(2) molecule, which is unbound within unrestricted Hartree-Fock (UHF). With a UHF mid R:Phi(r), spin contamination causes large systematic errors and long equilibration times in AFQMC in the intermediate, bond-breaking region. The spin-projection method eliminates these problems and delivers an accurate potential energy curve from equilibrium to the dissociation limit using the UHF mid R:Phi(r). Realistic potential energy curves are obtained with a cc-pVQZ basis. The calculated spectroscopic constants are in excellent agreement with experiment.

A study of H+H2 and several H-bonded molecules by phaseless auxiliary-field quantum Monte Carlo with plane wave and Gaussian basis sets.

The authors present phaseless auxiliary-field (AF) quantum Monte Carlo (QMC) calculations of the ground states of some hydrogen-bonded systems. These systems were selected to test and benchmark different aspects of the new phaseless AF QMC method. They include the transition state of H+H(2) near the equilibrium geometry and in the van der Walls limit, as well as the H(2)O, OH, and H(2)O(2) molecules. Most of these systems present significant challenges for traditional independent-particle electronic structure approaches, and many also have exact results available. The phaseless AF QMC method is used either with a plane wave basis with pseudopotentials or with all-electron Gaussian basis sets. For some systems, calculations are done with both to compare and characterize the performance of AF QMC under different basis sets and different Hubbard-Stratonovich decompositions. Excellent results are obtained using as input single Slater determinant wave functions taken from independent-particle calculations. Comparisons of the Gaussian based AF QMC results with exact full configuration interaction show that the errors from controlling the phase problem with the phaseless approximation are small. At the large basis-size limit, the AF QMC results using both types of basis sets are in good agreement with each other and with experimental values.

Bond breaking with auxiliary-field quantum Monte Carlo.

Bond stretching mimics different levels of electron correlation and provides a challenging test bed for approximate many-body computational methods. Using the recently developed phaseless auxiliary-field quantum Monte Carlo (AF QMC) method, we examine bond stretching in the well-studied molecules BH and N(2) and in the H(50) chain. To control the sign/phase problem, the phaseless AF QMC method constrains the paths in the auxiliary-field path integrals with an approximate phase condition that depends on a trial wave function. With single Slater determinants from unrestricted Hartree-Fock as trial wave function, the phaseless AF QMC method generally gives better overall accuracy and a more uniform behavior than the coupled cluster CCSD(T) method in mapping the potential-energy curve. In both BH and N(2), we also study the use of multiple-determinant trial wave functions from multiconfiguration self-consistent-field calculations. The increase in computational cost versus the gain in statistical and systematic accuracy are examined. With such trial wave functions, excellent results are obtained across the entire region between equilibrium and the dissociation limit.

Auxiliary-field quantum Monte Carlo calculations of molecular systems with a Gaussian basis.

We extend the recently introduced phaseless auxiliary-field quantum Monte Carlo (QMC) approach to any single-particle basis and apply it to molecular systems with Gaussian basis sets. QMC methods in general scale favorably with the system size as a low power. A QMC approach with auxiliary fields, in principle, allows an exact solution of the Schrodinger equation in the chosen basis. However, the well-known sign/phase problem causes the statistical noise to increase exponentially. The phaseless method controls this problem by constraining the paths in the auxiliary-field path integrals with an approximate phase condition that depends on a trial wave function. In the present calculations, the trial wave function is a single Slater determinant from a Hartree-Fock calculation. The calculated all-electron total energies show typical systematic errors of no more than a few millihartrees compared to exact results. At equilibrium geometries in the molecules we studied, this accuracy is roughly comparable to that of coupled cluster with single and double excitations and with noniterative triples [CCSD(T)]. For stretched bonds in H(2)O, our method exhibits a better overall accuracy and a more uniform behavior than CCSD(T).

Auxiliary-field quantum Monte Carlo study of first- and second-row post-d elements.

A series of calculations for the first- and second-row post-d elements (Ga-Br and In-I) are presented using the phaseless auxiliary-field quantum Monte Carlo (AF QMC) method. This method is formulated in a Hilbert space defined by any chosen one-particle basis and maps the many-body problem into a linear combination of independent-particle solutions with external auxiliary fields. The phase/sign problem is handled approximately by the phaseless formalism using a trial wave function, which in our calculations was chosen to be the Hartree-Fock solution. We used the consistent correlated basis sets of Peterson et al. [J. Chem. Phys. 119, 11099 (2003); 119, 11113 (2003)], which employ a small-core relativistic pseudopotential. The AF QMC results are compared with experiment and with those from density functional (generalized gradient approximation and B3LYP) and CCSD(T) calculations. The AF QMC total energies agree with CCSD(T) to within a few millihartrees across the systems and over several basis sets. The calculated atomic electron affinities, ionization energies, and spectroscopic properties of dimers are, at large basis sets, in excellent agreement with experiment.