Efficient Monte Carlo Calculations of the One-Body Density

An alternative Monte Carlo estimator for the one-body density rho(r) is presented. This estimator has a simple form and can be readily used in any type of Monte Carlo simulation. Comparisons with the usual regularization of the delta-function on a grid show that the statistical errors are greatly reduced. Furthermore, our expression allows accurate calculations of the density at any point in space, even in the regions never visited during the Monte Carlo simulation. The method is illustrated with the computation of accurate Variational Monte Carlo electronic densities for the Helium atom (1D curve) and for the water dimer (3D grid containing up to 51x51x51=132651 points).Comment: 12 pages with 3 postscript figure

Optimizing the energy with quantum Monte Carlo: A lower numerical scaling for Jastrow-Slater expansions

We present an improved formalism for quantum Monte Carlo calculations of energy derivatives and properties (e.g. the interatomic forces), with a multideterminant Jastrow-Slater function. As a function of the number $N_e$ of Slater determinants, the numerical scaling of $O(N_e)$ per derivative we have recently reported is here lowered to $O(N_e)$ for the entire set of derivatives. As a function of the number of electrons $N$, the scaling to optimize the wave function and the geometry of a molecular system is lowered to $O(N^3)+O(N N_e)$, the same as computing the energy alone in the sampling process. The scaling is demonstrated on linear polyenes up to C$_{60}$H$_{62}$ and the efficiency of the method is illustrated with the structural optimization of butadiene and octatetraene with Jastrow-Slater wave functions comprising as many as 200000 determinants and 60000 parameters

Zero-Variance Zero-Bias Principle for Observables in quantum Monte Carlo: Application to Forces

A simple and stable method for computing accurate expectation values of observable with Variational Monte Carlo (VMC) or Diffusion Monte Carlo (DMC) algorithms is presented. The basic idea consists in replacing the usual ``bare'' estimator associated with the observable by an improved or ``renormalized'' estimator. Using this estimator more accurate averages are obtained: Not only the statistical fluctuations are reduced but also the systematic error (bias) associated with the approximate VMC or (fixed-node) DMC probability densities. It is shown that improved estimators obey a Zero-Variance Zero-Bias (ZVZB) property similar to the usual Zero-Variance Zero-Bias property of the energy with the local energy as improved estimator. Using this property improved estimators can be optimized and the resulting accuracy on expectation values may reach the remarkable accuracy obtained for total energies. As an important example, we present the application of our formalism to the computation of forces in molecular systems. Calculations of the entire force curve of the H$_2$,LiH, and Li$_2$ molecules are presented. Spectroscopic constants $R_e$ (equilibrium distance) and $\omega_e$ (harmonic frequency) are also computed. The equilibrium distances are obtained with a relative error smaller than 1%, while the harmonic frequencies are computed with an error of about 10%

Introduction to the variational and diffusion Monte Carlo methods

We provide a pedagogical introduction to the two main variants of real-space quantum Monte Carlo methods for electronic-structure calculations: variational Monte Carlo (VMC) and diffusion Monte Carlo (DMC). Assuming no prior knowledge on the subject, we review in depth the Metropolis-Hastings algorithm used in VMC for sampling the square of an approximate wave function, discussing details important for applications to electronic systems. We also review in detail the more sophisticated DMC algorithm within the fixed-node approximation, introduced to avoid the infamous Fermionic sign problem, which allows one to sample a more accurate approximation to the ground-state wave function. Throughout this review, we discuss the statistical methods used for evaluating expectation values and statistical uncertainties. In particular, we show how to estimate nonlinear functions of expectation values and their statistical uncertainties.Comment: Advances in Quantum Chemistry, 2015, Electron Correlation in Molecules -- ab initio Beyond Gaussian Quantum Chemistry, pp.000

Zero-variance zero-bias quantum Monte Carlo estimators of the spherically and system-averaged pair density

We construct improved quantum Monte Carlo estimators for the spherically- and system-averaged electron pair density (i.e. the probability density of finding two electrons separated by a relative distance u), also known as the spherically-averaged electron position intracule density I(u), using the general zero-variance zero-bias principle for observables, introduced by Assaraf and Caffarel. The calculation of I(u) is made vastly more efficient by replacing the average of the local delta-function operator by the average of a smooth non-local operator that has several orders of magnitude smaller variance. These new estimators also reduce the systematic error (or bias) of the intracule density due to the approximate trial wave function. Used in combination with the optimization of an increasing number of parameters in trial Jastrow-Slater wave functions, they allow one to obtain well converged correlated intracule densities for atoms and molecules. These ideas can be applied to calculating any pair-correlation function in classical or quantum Monte Carlo calculations.Comment: 13 pages, 9 figures, published versio

Simple formalism for efficient derivatives and multi-determinant expansions in quantum Monte Carlo

We present a simple and general formalism to compute efficiently the derivatives of a multi-determinant Jastrow-Slater wave function, the local energy, the interatomic forces, and similar quantities needed in quantum Monte Carlo. Through a straightforward manipulation of matrices evaluated on the occupied and virtual orbitals, we obtain an efficiency equivalent to algorithmic differentiation in the computation of the interatomic forces and the optimization of the orbital paramaters. Furthermore, for a large multi-determinant expansion, the significant computational gain recently reported for the calculation of the wave function is here improved and extended to all local properties in both all-electron and pseudopotential calculations.Comment: 15 pages, 3 figure

Ab initio lifetime correction to scattering states for time-dependent electronic-structure calculations with incomplete basis sets

We propose a method for obtaining effective lifetimes of scattering electronic states for avoiding the artificially confinement of the wave function due to the use of incomplete basis sets in time-dependent electronic-structure calculations of atoms and molecules. In this method, using a fitting procedure, the lifetimes are extracted from the spatial asymptotic decay of the approximate scattering wave functions obtained with a given basis set. The method is based on a rigorous analysis of the complex-energy solutions of the Schr{\"o}dinger equation. It gives lifetimes adapted to any given basis set without using any empirical parameters. The method can be considered as an ab initio version of the heuristic lifetime model of Klinkusch et al. [J. Chem. Phys. 131, 114304 (2009)]. The method is validated on the H and He atoms using Gaussian-type basis sets for calculation of high-harmonic-generation spectra

Not only P-glycoprotein: amplification of the ABCB1-containing chromosome region 7q21 confers multidrug resistance upon cancer cells by coordinated overexpression of an assortment of resistance-related proteins

The development of drug resistance continues to be a dominant hindrance toward curative cancer treatment. Overexpression of a wide-spectrum of ATP-dependent efflux pumps, and in particular of ABCB1 (P-glycoprotein or MDR1) is a well-known resistance mechanism for a plethora of cancer chemotherapeutics including for example taxenes, anthracyclines, Vinca alkaloids, and epipodopyllotoxins, demonstrated by a large array of published papers, both in tumor cell lines and in a variety of tumors, including various solid tumors and hematological malignancies. Upon repeated or even single dose treatment of cultured tumor cells or tumors in vivo with anti-tumor agents such as paclitaxel and doxorubicin, increased ABCB1 copy number has been demonstrated, resulting from chromosomal amplification events at 7q11.2-21 locus, leading to marked P-glycoprotein overexpression, and multidrug resistance (MDR). Clearly however, additional mechanisms such as single nucleotide polymorphisms (SNPs) and epigenetic modifications have shown a role in the overexpression of ABCB1 and of other MDR efflux pumps. However, notwithstanding the design of 4 generations of ABCB1 inhibitors and the wealth of information on the biochemistry and substrate specificity of ABC transporters, translation of this vast knowledge from the bench to the bedside has proven to be unexpectedly difficult. Many studies show that upon repeated treatment schedules of cell cultures or tumors with taxenes and anthracyclines as well as other chemotherapeutic drugs, amplification, and/or overexpression of a series of genes genomically surrounding the ABCB1 locus, is observed. Consequently, altered levels of other proteins may contribute to the establishment of the MDR phenotype, and lead to poor clinical outcome. Thus, the genes contained in this ABCB1 amplicon including ABCB4, SRI, DBF4, TMEM243, and RUNDC3B are overexpressed in many cancers, and especially in MDR tumors, while TP53TG1 and DMTF1 are bona fide tumor suppressors. This review describes the role of these genes in cancer and especially in the acquisition of MDR, elucidates possible connections in transcriptional regulation (co-amplification/repression) of genes belonging to the same ABCB1 amplicon region, and delineates their novel emerging contributions to tumor biology and possible strategies to overcome cancer MDR

Curing basis-set convergence of wave-function theory using density-functional theory: a systematically improvable approach

The present work proposes to use density-functional theory (DFT) to correct for the basis-set error of wave-function theory (WFT). One of the key ideas developed here is to define a range-separation parameter which automatically adapts to a given basis set. The derivation of the exact equations are based on the Levy-Lieb formulation of DFT, which helps us to define a complementary functional which corrects uniquely for the basis-set error of WFT. The coupling of DFT and WFT is done through the definition of a real-space representation of the electron-electron Coulomb operator projected in a one-particle basis set. Such an effective interaction has the particularity to coincide with the exact electron-electron interaction in the limit of a complete basis set, and to be finite at the electron-electron coalescence point when the basis set is incomplete. The non-diverging character of the effective interaction allows one to define a mapping with the long-range interaction used in the context of range-separated DFT and to design practical approximations for the unknown complementary functional. Here, a local-density approximation is proposed for both full-configuration-interaction (FCI) and selected configuration-interaction approaches. Our theory is numerically tested to compute total energies and ionization potentials for a series of atomic systems. The results clearly show that the DFT correction drastically improves the basis-set convergence of both the total energies and the energy differences. For instance, a sub kcal/mol accuracy is obtained from the aug-cc-pVTZ basis set with the method proposed here when an aug-cc-pV5Z basis set barely reaches such a level of accuracy at the near FCI level