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

Ali Alavi

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

ResearcherUnclaimed source record
physics.comp-phcond-mat.str-elphysics.chem-phcond-mat.mtrl-scicond-mat.quant-gasquant-phcond-mat.mes-hallcond-mat.stat-mechcond-mat.supr-con

Catalog footprint

What is connected

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

27 published item(s)

preprint2022arXiv

Ab initio wavefunction analysis of electron removal quasi-particle state of NdNiO$_2$ with fully correlated quantum chemical methods

The discovery of superconductivity in hole-doped infinite-layer NdNiO$_2$ -- a transition metal (TM) oxide that is both isostructural and isoelectronic to cuprate superconductors -- has lead to renewed enthusiasm in the hope of understanding the origin of unconventional superconductivity. Here, we investigate the electron-removal states in infinite-layered Ni$^{1+}$ oxide, NdNiO$_2$, which mimics hole-doping, with the state-of-the-art many-body multireference quantum chemistry methods. From the analysis of the many-body wavefunction, we find that the hole-doped $d^8$ ground state of NdNiO$_2$ is very different from the $d^8$ ground state in isostructural cuprate analog CaCuO$_2$, although the parent $d^9$ ground states are for the most part identical. We show that the doped hole in NdNiO$_2$ mainly localizes on the Ni $3d_{x^2-y^2}$ orbital to form a closed-shell singlet, and this singlet configuration contributes to $\sim$40% of the wavefunction. In contrast, in CaCuO$_2$ the Zhang-Rice singlet configurations contribute to $\sim$65% of the wavefunction. With the help of the quantum information concept of entanglement entropy, we quantify the different types of electronic correlations in the nickelate and cuprate compounds and find that the dynamic radial-type correlations within the Ni $d$ manifold are persistent in hole-doped NdNiO$_2$. As a result, the $d^8$ multiplet effects are stronger and the additional hole foot-print is more three-dimensional in NdNiO$_2$. Our analysis shows that the most commonly used three-band Hubbard model employed to express the doped scenario in cuprates represents $\sim$90% of the $d^8$ wavefunction for CaCuO$_2$, but such a model grossly approximates the $d^8$ wavefunction for NdNiO$_2$ as it only stands for $\sim$60% of the wavefunction.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

General analytical nuclear force and molecular potential energy surface from full configuration interaction quantum Monte Carlo

Full configuration interaction quantum Monte Carlo (FCIQMC) is a state-of-the-art stochastic electronic structure method, providing a methodology to compute FCI-level state energies of molecular systems within a quantum chemical basis. However, especially to probe {\em dynamics} at the FCIQMC level, it is necessary to devise more efficient schemes to produce nuclear forces and potential energy surfaces (PES) from FCIQMC. In this work, we derive the general formula for nuclear force from FCIQMC, and clarify different contributions of the total force. This method to obtain FCIQMC forces eliminates previous restrictions, and can be used with frozen core approximation and free selection of orbitals, making it promising for more efficient nuclear force calculations. After numerical check of this procedure on the binding curve of N$_2$ molecule, we use the FCIQMC energy and force to obtain the full-dimensional ground state PES of water molecule via Gaussian processes regression. The new water FCIQMC PES can be used as the basis for H$_2$O ground state nuclear dynamics, structure optimization, and rotation-vibrational spectrum calculation.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

General embedded cluster protocol for accurate modeling of oxygen vacancies in metal-oxides

The O vacancy (Ov) formation energy, $E_\textrm{Ov}$, is an important property of a metal-oxide, governing its performance in applications such as fuel cells or heterogeneous catalysis. These defects are routinely studied with density functional theory (DFT). However, it is well-recognized that standard DFT formulations (e.g. the generalized gradient approximation) are insufficient for modeling the Ov, requiring higher levels of theory. The embedded cluster method offers a promising approach to compute $E_\textrm{Ov}$ accurately, giving access to all electronic structure methods. Central to this approach is the construction of quantum(-mechanically treated) clusters placed within suitable embedding environments. Unfortunately, current approaches to constructing the quantum clusters either require large system sizes, preventing application of high-level methods, or require significant manual input, preventing investigations of multiple systems simultaneously. In this work, we present a systematic and general quantum cluster design protocol that can determine small converged quantum clusters for studying the Ov in metal-oxides with accurate methods such as local coupled cluster with single, double and perturbative triple excitations [CCSD(T)]. We apply this protocol to study the Ov in the bulk and surface planes of rutile TiO2 and rocksalt MgO, producing the first accurate and well-converged determinations of $E_\textrm{Ov}$ with this method. These reference values are used to benchmark exchange-correlation functionals in DFT and we find that all studied functionals underestimate $E_\textrm{Ov}$, with the average error decreasing along the rungs of Jacob's ladder. This protocol is automatable for high-throughput calculations and can be generalized to study other point defects or adsorbates.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Combined unitary and symmetric group approach applied to low-dimensional spin systems

A novel combined unitary and symmetric group approach is used to study the spin-$\frac{1}{2}$ Heisenberg model and related Fermionic systems in a spin-adapted representation, using a linearly-parameterised Ansatz for the many-body wave function. We show that a more compact ground state wave function representation is obtained when combining the symmetric group, $\mathcal{S}_n$, in the form of permutations of the underlying lattice site ordering, with the cumulative spin-coupling based on the unitary group, $\mathrm{U}(n)$. In one-dimensional systems the observed compression of the wave function is reminiscent of block-spin renormalization group approaches, and allows us to study larger lattices (here taken up to 80 sites) with the spin-adapted full configuration interaction quantum Monte Carlo method, which benefits from the sparsity of the Hamiltonian matrix and the corresponding sampled eigenstates that emerge from the reordering. We find that in an optimal lattice ordering the configuration state function with highest weight already captures with high accuracy the spin-spin correlation function of the exact ground state wave function. This feature is found for more general lattice models, such as the Hubbard model, and ab initio quantum chemical models, in this work exemplified by a one-dimensional hydrogen chain. We also provide numerical evidence that the optimal lattice ordering for the unitary group approach is not generally equivalent to the optimal ordering obtained for methods based on matrix-product states, such as the density-matrix renormalization group approach.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A comparative study using state-of-the-art electronic structure theories on solid hydrogen phases under high pressures

Identifying the atomic structure and properties of solid hydrogen under high pressures is a long-standing problem of high-pressure physics with far-reaching significance in planetary and materials science. Determining the pressure-temperature phase diagram of hydrogen is challenging for experiment and theory due to the extreme conditions and the required accuracy in the quantum mechanical treatment of the constituent electrons and nuclei, respectively. Here, we demonstrate explicitly that coupled cluster theory can serve as a computationally efficient theoretical tool to predict solid hydrogen phases with high accuracy. We present a first principles study of solid hydrogen phases at pressures ranging from 100 to 450 GPa. The computed static lattice enthalpies are compared to state-of-the-art diffusion Monte Carlo results and density functional theory calculations. Our coupled cluster theory results for the most stable phases including C2/c-24 and P21/c-24 are in good agreement with those obtained using diffusion Monte Carlo, with the exception of Cmca-4, which is predicted to be significantly less stable. We discuss the scope of the employed methods and how they can contribute as efficient and complementary theoretical tools to solve the long-standing puzzle of understanding solid hydrogen phases at high pressures.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

NECI: N-Electron Configuration Interaction with emphasis on state-of-the-art stochastic methods

We present NECI, a state-of-the-art implementation of the Full Configuration Interaction Quantum Monte Carlo algorithm, a method based on a stochastic application of the Hamiltonian matrix on a sparse sampling of the wave function. The program utilizes a very powerful parallelization and scales efficiently to more than 24000 CPU cores. In this paper, we describe the core functionalities of NECI and recent developments. This includes the capabilities to calculate ground and excited state energies, properties via the one- and two-body reduced density matrices, as well as spectral and Green's functions for ab initio and model systems. A number of enhancements of the bare FCIQMC algorithm are available within NECI, allowing to use a partially deterministic formulation of the algorithm, working in a spin-adapted basis or supporting transcorrelated Hamiltonians. NECI supports the FCIDUMP file format for integrals, supplying a convenient interface to numerous quantum chemistry programs and it is licensed under GPL-3.0.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Structural and electronic properties of solid molecular hydrogen from many-electron theories

We study the structural and electronic properties of phase III of solid hydrogen using accurate many-electron theories and compare to state-of-the-art experimental findings. The atomic structures of phase III modelled by C2/c-24 crystals are fully optimized on the level of second-order perturbation theory, demonstrating that previously employed structures optimized on the level of approximate density functionals exhibit errors in the H$_2$ bond lengths that cause significant discrepancies in the computed quasi particle band gaps and vibrational frequencies compared to experiment. Using the newly optimized atomic structures, we study the band gap closure and change in vibrational frequencies as a function of pressure. Our findings are in good agreement with recent experimental observations and may prove useful in resolving long-standing discrepancies between experimental estimates of metallization pressures possibly caused by disagreeing pressure calibrations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Adaptive Shift Method in Full Configuration Interaction Quantum Monte Carlo: Development and Applications

In a recent paper, we proposed the adaptive shift method for correcting the undersampling bias of the initiator-FCIQMC. The method allows faster convergence with the number of walkers to the FCI limit than the normal initiator method, particularly for large systems. In its application to strongly correlated molecules, however, the method is prone to overshooting the FCI energy at intermediate walker numbers, with convergence to the FCI limit from below. In this paper, we present a solution to the overshooting problem in strongly correlated molecules, as well as further accelerating convergence to the FCI energy. This is achieved by offsetting the reference energy to a value typically below the Hartree-Fock energy but above the exact energy. This offsetting procedure does not change the exactness property of the algorithm, namely convergence to the exact FCI solution in the large-walker limit, but at its optimal value greatly accelerates convergence. There is no overhead cost associated with this offsetting procedure, and is therefore a pure and substantial computational gain. We illustrate the behavior of this offset adaptive shift method by applying it to the N$_2$ molecule, the ozone molecule at three different geometries (equilibrium open minimum, a hypothetical ring minimum, and a transition state) in three basis sets (cc-pV$X$Z, $X$=D,T,Q), and the chromium dimer in cc-pVDZ basis set, correlating 28 electrons in 76 orbitals. We show that in most cases the offset adaptive shift method converges much faster than both the normal initiator method and the original adaptive shift method.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The color center singlet state of oxygen vacancies in TiO$_2$

Oxygen vacancies are ubiquitous in TiO$_2$ and play key roles in catalysis and magnetism applications.Despite being extensively investigated, the electronic structure of oxygen vacancies in TiO$_2$ remains controversial both experimentally and theoretically.Here we report a study of a neutral oxygen vacancy in TiO$_2$ using state-of-the-art quantum chemical electronic structure methods.We find that the ground state is a color center singlet state in both the rutile and the anatase phase of TiO$_2$. Specifically, embedded CCSD(T) calculations find, for an oxygen vacancy in rutile, that the lowest triplet state energy is 0.6 eV above the singlet state, and in anatase the triplet state energy is higher by 1.4 eV. Our study provides fresh insights on the electronic structure of the oxygen vacancy in TiO$_2$, clarifying earlier controversies and potentially inspiring future studies of defects with correlated wave function theories.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Small polarons and the Janus nature of $\text{TiO}_\text{2}(110)$

Polarons are ubiquitous in many semiconductors and have been linked with conductivity and optical response of materials for photovoltaics and heterogeneous catalysis, yet how surface polarons influence adsorption remains unclear. Here, by modelling the surface of rutile titania using density functional theory, we reveal the effect of small surface polarons on water adsorption, dissociation, and hydrogen bonding. On the one hand the presence of such polarons significantly suppresses dissociation of water molecules that are bonded directly to polaronic sites. On the other hand, polarons facilitate water dissociation at certain non-polaronic sites. Furthermore, polarons strengthen hydrogen bonds, which in turn affects water dissociation in hydrogen bonded overlayer structures. This study reveals that polarons at the rutile surface have complex, multi-faceted, effects on water adsorption, dissociation and hydrogen bonding, highlighting the importance of polarons on water structure and dynamics on such surfaces. We expect that many of the physical properties of surface polarons identified here will apply more generally to surfaces and interfaces that can host small polarons, beyond titania.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Assessment of Multireference Approaches to Explicitly Correlated Full Configuration Interaction Quantum Monte Carlo

The Full Configuration Interaction Quantum Monte Carlo (FCIQMC) method has proved able to provide near-exact solutions to the electronic Schrödinger equation within a finite orbital basis set, without relying on an expansion about a reference state. However, a drawback to the approach is that being based on an expansion of Slater determinants, the FCIQMC method suffers from a basis set incompleteness error that decays very slowly with the size of the employed single particle basis. The FCIQMC results obtained in a small basis set can be improved significantly with explicitly correlated techniques. Here, we present a study that assesses and compares two contrasting `universal' explicitly correlated approaches that fit into the FCIQMC framework; the $[2]_{R12}$ method of Valeev {\em et al.}, and the explicitly correlated canonical transcorrelation approach of Yanai {\em et al}. The former is an {\em a posteriori} internally-contracted perturbative approach, while the latter transforms the Hamiltonian prior to the FCIQMC simulation. These comparisons are made across the 55 molecules of the G1 standard set. We found that both methods consistently reduce the basis set incompleteness, for accurate atomization energies in small basis sets, reducing the error from 28~mE$_{\text h}$ to 3-4~mE$_{\text h}$. While many of the conclusions hold in general for any combination of multireference approaches with these methodologies, we also consider FCIQMC-specific advantages of each approach.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Combining internally contracted states and matrix product states to perform multireference perturbation theory

We present two efficient and intruder-free methods for treating dynamic correlation on top of general multi-configuration reference wave functions---including such as obtained by the density matrix renormalization group (DMRG) with large active spaces. The new methods are the second order variant of the recently proposed multi-reference linearized coupled cluster method (MRLCC) [S. Sharma, A. Alavi, J. Chem. Phys. 143, 102815 (2015)], and of N-electron valence perturbation theory (NEVPT2), with expected accuracies similar to MRCI+Q and (at least) CASPT2, respectively. Great efficiency gains are realized by representing the first-order wave function with a combination of internal contraction (IC) and matrix product state perturbation theory (MPSPT). With this combination, only third order reduced density matrices (RDMs) are required. Thus, we obviate the need for calculating (or estimating) RDMs of fourth or higher order; these had so far posed a severe bottleneck for dynamic correlation treatments involving the large active spaces accessible to DMRG. Using several benchmark systems, including first and second row containing small molecules, Cr$_2$, pentacene and oxo-Mn(Salen), we shown that active spaces containing at least 30 orbitals can be treated using this method. On a single node, MRLCC2 and NEVPT2 calculations can be performed with over 550 and 1100 virtual orbitals, respectively. We also critically examine the errors incurred due to the three sources of errors introduced in the present implementation - calculating second order instead of third order energy corrections, use of internal contraction and approximations made in the reference wavefunction due to DMRG.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

An excited-state approach within full configuration interaction quantum Monte Carlo

We present a new approach to calculate excited states with the full configuration interaction quantum Monte Carlo (FCIQMC) method. The approach uses a Gram-Schmidt procedure, instantaneously applied to the stochastically evolving distributions of walkers, to orthogonalize higher energy states against lower energy ones. It can thus be used to study several of the lowest-energy states of a system within the same symmetry. This additional step is particularly simple and computationally inexpensive, requiring only a small change to the underlying FCIQMC algorithm. No trial wave functions or partitioning of the space is needed. The approach should allow excited states to be studied for systems similar to those accessible to the ground-state method, due to a comparable computational cost. As a first application we consider the carbon dimer in basis sets up to quadruple-zeta quality, and compare to existing results where available.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Analytic nuclear forces and molecular properties from full configuration interaction quantum Monte Carlo

Unbiased stochastic sampling of the one- and two-body reduced density matrices is achieved in full configuration interaction quantum Monte Carlo with the introduction of a second, "replica" ensemble of walkers, whose population evolves in imaginary time independently from the first, and which entails only modest additional computational overheads. The matrices obtained from this approach are shown to be representative of full configuration-interaction quality, and hence provide a realistic opportunity to achieve high-quality results for a range of properties whose operators do not necessarily commute with the hamiltonian. A density-matrix formulated quasi-variational energy estimator having been already proposed and investigated, the present work extends the scope of the theory to take in studies of analytic nuclear forces, molecular dipole moments and polarisabilities, with extensive comparison to exact results where possible. These new results confirm the suitability of the sampling technique and, where sufficiently large basis sets are available, achieve close agreement with experimental values, expanding the scope of the method to new areas of investigation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Krylov-projected quantum Monte Carlo

We present an approach to the calculation of arbitrary spectral, thermal and excited state properties within the full configuration interaction quantum Monte Carlo framework. This is achieved via an unbiased projection of the Hamiltonian eigenvalue problem into a space of stochastically sampled Krylov vectors, thus enabling the calculation of real-frequency spectral and thermal properties and avoiding explicit analytic continuation. We use this approach to calculate temperature-dependent properties and one- and two-body spectral functions for various Hubbard models, as well as isolated excited states in ab initio systems.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Multireference linearized Coupled Cluster theory for strongly correlated systems using Matrix Product States

We propose a multireference linearized coupled cluster theory using matrix product states (MPS-LCC) which provides remarkably accurate ground-state energies, at a computational cost that has the same scaling as multireference configuration interaction singles and doubles (MRCISD), for a wide variety of electronic Hamiltonians. These range from first-row dimers at equilibrium and stretched geometries, to highly multireference systems such as the chromium dimer and lattice models such as periodic two-dimensional 1-band and 3-band Hubbard models. The MPS-LCC theory shows a speed up of several orders of magnitude over the usual DMRG algorithm while delivering energies in excellent agreement with converged DMRG calculations. Also, in all the benchmark calculations presented here MPS-LCC outperformed the commonly used multi-reference quantum chemistry methods in some cases giving energies in excess of an order of magnitude more accurate. As a size-extensive method that can treat large active spaces, MPS-LCC opens up the use of multireference quantum chemical techniques in strongly-correlated \emph{ab-initio} Hamiltonians, including two and three-dimensional solids.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Semi-stochastic full configuration interaction quantum Monte Carlo: developments and application

We expand upon the recent semi-stochastic adaptation to full configuration interaction quantum Monte Carlo (FCIQMC). We present an alternate method for generating the deterministic space without a priori knowledge of the wave function and present stochastic efficiencies for a variety of both molecular and lattice systems. The algorithmic details of an efficient semi-stochastic implementation are presented, with particular consideration given to the effect that the adaptation has on parallel performance in FCIQMC. We further demonstrate the benefit for calculation of reduced density matrices in FCIQMC through replica sampling, where the semi-stochastic adaptation seems to have even larger efficiency gains. We then combine these ideas to produce explicitly correlated corrected FCIQMC energies for the Beryllium dimer, for which stochastic errors on the order of wavenumber accuracy are achievable.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Stochastic multi-configurational self-consistent field theory

The multi-configurational self-consistent field theory is considered the standard starting point for almost all multireference approaches required for strongly-correlated molecular problems. The limitation of the approach is generally given by the number of strongly-correlated orbitals in the molecule, as its cost will grow exponentially with this number. We present a new multi-configurational self-consistent field approach, wherein linear determinant coefficients of a multi-configurational wavefunction are optimized via the stochastic full configuration interaction quantum Monte Carlo technique at greatly reduced computational cost, with non-linear orbital rotation parameters updated variationally based on this sampled wavefunction. This extends this approach to strongly-correlated systems with far larger active spaces than it is possible to treat by conventional means. By comparison with this traditional approach, we demonstrate that the introduction of stochastic noise in both the determinant amplitudes and the gradient and Hessian of the orbital rotations does not preclude robust and reliable convergence of the orbital optimization. It can even improve the ability to avoid convergence to local minima in the orbital space, and therefore aid in finding variationally lower-energy solutions. We consider the effect on the convergence of the orbitals as the number of walkers and the sampling time within the active space increases, as well as the effect on the final energy and error. The scope of the new protocol is demonstrated with a study of the increasingly strongly correlated electronic structure in a series of polycyclic aromatic hydrocarbons, up to the large coronene molecule in a complete active space of 24 pi electrons in 24 orbitals, requiring only modest computational resources.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Unbiased Reduced Density Matrices and Electronic Properties from Full Configuration Interaction Quantum Monte Carlo

Properties that are necessarily formulated within pure (symmetric) expectation values are difficult to calculate for projector quantum Monte Carlo approaches, but are critical in order to compute many of the important observable properties of electronic systems. Here, we investigate an approach for the sampling of unbiased reduced density matrices within the Full Configuration Interaction Quantum Monte Carlo dynamic, which requires only small computational overheads. This is achieved via an independent replica population of walkers in the dynamic, sampled alongside the original population. The resulting reduced density matrices are free from systematic error (beyond those present via constraints on the dynamic itself), and can be used to compute a variety of expectation values and properties, with rapid convergence to an exact limit. A quasi-variational energy estimate derived from these density matrices is proposed as an accurate alternative to the projected estimator for multiconfigurational wavefunctions, while its variational property could potentially lend itself to accurate extrapolation approaches in larger systems.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Explicitly correlated plane waves: Accelerating convergence in periodic wavefunction expansions

We present an investigation into the use of an explicitly correlated plane wave basis for periodic wavefunction expansions at the level of second-order Møller-Plesset perturbation theory (MP2). The convergence of the electronic correlation energy with respect to the one-electron basis set is investigated and compared to conventional MP2 theory in a finite homogeneous electron gas model. In addition to the widely used Slater-type geminal correlation factor, we also derive and investigate a novel correlation factor that we term Yukawa-Coulomb. The Yukawa-Coulomb correlation factor is motivated by analytic results for two electrons in a box and allows for a further improved convergence of the correlation energies with respect to the employed basis set. We find the combination of the infinitely delocalized plane waves and local short-ranged geminals provides a complementary, and rapidly convergent basis for the description of periodic wavefunctions. We hope that this approach will expand the scope of discrete wavefunction expansions in periodic systems.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Linear-scaling and parallelizable algorithms for stochastic quantum chemistry

For many decades, quantum chemical method development has been dominated by algorithms which involve increasingly complex series of tensor contractions over one-electron orbital spaces. Procedures for their derivation and implementation have evolved to require the minimum amount of logic and rely heavily on computationally efficient library-based matrix algebra and optimized paging schemes. In this regard, the recent development of exact stochastic quantum chemical algorithms to reduce computational scaling and memory overhead requires a contrasting algorithmic philosophy, but one which when implemented efficiently can often achieve higher accuracy/cost ratios with small random errors. Additionally, they can exploit the continuing trend for massive parallelization which hinders the progress of deterministic high-level quantum chemical algorithms. In the Quantum Monte Carlo community, stochastic algorithms are ubiquitous but the discrete Fock space of quantum chemical methods is often unfamiliar, and the methods introduce new concepts required for algorithmic efficiency. In this paper, we explore these concepts and detail an algorithm used for Full Configuration Interaction Quantum Monte Carlo (FCIQMC), which is implemented and available in MOLPRO and as a standalone code, and is designed for high-level parallelism and linear-scaling with walker number. Many of the algorithms are also in use in, or can be transferred to, other stochastic quantum chemical methods and implementations. We apply these algorithms to the strongly correlated Chromium dimer, to demonstrate their efficiency and parallelism.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

An explicitly correlated approach to basis set incompleteness in Full Configuration Interaction Quantum Monte Carlo

By performing a stochastic dynamic in a space of Slater determinants, the Full Configuration Interaction Quantum Monte Carlo (FCIQMC) method has been able to obtain energies which are essentially free from systematic error to the basis set correlation energy, within small and systematically improvable errorbars. However, the weakly exponential scaling with basis size makes converging the energy with respect to basis set costly and in larger systems, impossible. To ameliorate these basis set issues, here we use perturbation theory to couple the FCIQMC wave function to an explicitly correlated strongly orthogonal basis of geminals, following the [2]_{\textrm{R12}} approach of Valeev {\em et al.}. The required one- and two-particle density matrices are computed on-the-fly during the FCIQMC dynamic, using a sampling procedure which incurs relatively little additional computation expense. The F12 energy corrections are shown to converge rapidly as a function of sampling, both in imaginary time, and number of walkers. Our pilot calculations on the binding curve for carbon dimer, which exhibits strong correlation effects as well as substantial basis set dependence, demonstrate that the accuracy of the FCIQMC-F12 method surpasses that of all previous FCIQMC calculations, and that the F12 correction improves accuracy equivalent to increasing the quality of the one-electron basis by two cardinal numbers.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Convergence of many-body wavefunction expansions using a plane wave basis: from the homogeneous electron gas to the solid state

Using the finite simulation-cell homogeneous electron gas (HEG) as a model, we investigate the convergence of the correlation energy to the complete basis set (CBS) limit in methods utilising plane-wave wavefunction expansions. Simple analytic and numerical results from second-order Møller-Plesset theory (MP2) suggest a 1/M decay of the basis-set incompleteness error where M is the number of plane waves used in the calculation, allowing for straightforward extrapolation to the CBS limit. As we shall show, the choice of basis set truncation when constructing many-electron wavefunctions is far from obvious, and here we propose several alternatives based on the momentum transfer vector, which greatly improve the rate of convergence. This is demonstrated for a variety of wavefunction methods, from MP2 to coupled-cluster doubles theory (CCD) and the random-phase approximation plus second-order screened exchange (RPA+SOSEX). Finite basis-set energies are presented for these methods and compared with exact benchmarks. A transformation can map the orbitals of a general solid state system onto the HEG plane wave basis and thereby allow application of these methods to more realistic physical problems.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Emergence of Critical Phenomena in Full Configuration Interaction Quantum Monte Carlo

There has been recent literature discussion on the origin and severity of the `sign problem' in full configuration interaction quantum Monte Carlo (FCIQMC) and its `initiator' adaptation (i-FCIQMC), methods of interest and potential because they allow for exact (FCI) ground-state solutions to be obtained often at a much reduced computational cost. In this study we aim to use a simple order parameter, describing the `sign structure' of the stochastic wavefunction representation, to empirically characterise the fundamentally different collective behaviour of the walker population in both methods.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Investigation of the Full Configuration Interaction Quantum Monte Carlo Method Using Homogeneous Electron Gas Models

Using the homogeneous electron gas (HEG) as a model, we investigate the sources of error in the `initiator' adaptation to Full Configuration Interaction Quantum Monte Carlo (i-FCIQMC), with a view to accelerating convergence. In particular we find that the fixed shift phase, where the walker number is allowed to grow slowly, can be used to effectively assess stochastic and initiator error. Using this approach we provide simple explanations for the internal parameters of an i-FCIQMC simulation. We exploit the consistent basis sets and adjustable correlation strength of the HEG to analyze properties of the algorithm, and present finite basis benchmark energies for N=14 over a range of densities $0.5 \leq r_s \leq 5.0$ a.u. A \emph{single-point extrapolation} scheme is introduced to produce complete basis energies for 14, 38 and 54 electrons. It is empirically found that, in the weakly correlated regime, the computational cost scales linearly with the plane wave basis set size, which is justifiable on physical grounds. We expect the fixed shift strategy to reduce the computational cost of many \iFCIQMC calculations of weakly correlated systems. In addition, we provide benchmarks for the electron gas, to be used by other quantum chemical methods in exploring periodic solid state systems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Full Configuration Interaction Perspective on the Homogeneous Electron Gas

Highly accurate results for the homogeneous electron gas (HEG) have only been achieved to date within a diffusion Monte Carlo (DMC) framework. Here, we introduce a newly developed stochastic technique, Full Configuration Interaction Quantum Monte Carlo (FCIQMC), which samples the exact wavefunction expanded in plane wave Slater determinants. Despite the introduction of a basis set incompleteness error, we obtain a finite-basis energy which is significantly, and variationally lower than any previously published work for the 54-electron HEG at $r_s$ = 0.5 a.u., in a Hilbert space of $10^{108}$ Slater determinants. At this value of $r_s$, as well as of 1.0 a.u., we remove the remaining basis set incompleteness error by extrapolation, yielding results comparable or better than state-of-the-art DMC backflow energies. In doing so, we demonstrate that it is possible to yield highly accurate results with the FCIQMC method in sizable periodic systems.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Dispersion interactions between semiconducting wires

The dispersion energy between extended molecular chains (or equivalently infinite wires) with non-zero band gaps is generally assumed to be expressible as a pair-wise sum of atom-atom terms which decay as $R^{-6}$. Using a model system of two parallel wires with a variable band gap, we show that this is not the case. The dispersion interaction scales as $z^{-5}$ for large interwire separations $z$, as expected for an insulator, but as the band gap decreases the interaction is greatly enhanced; while at shorter (but non-overlapping) separations it approaches a power-law scaling given by $z^{-2}$, \emph{i.e.} the dispersion interaction expected between \emph{metallic} wires. We demonstrate that these effects can be understood from the increasing length scale of the plasmon modes (charge fluctuations), and their increasing contribution to the molecular dipole polarizability and the dispersion interaction, as the band gaps are reduced. This result calls into question methods which invoke locality assumptions in deriving dispersion interactions between extended small-gap systems.

0 citations0 reviewsNo code linkedOpen access