Anthony Scemama

Modeling Charge Resonance in Cationic Molecular Clusters: Combining DFT-Tight Binding with Configuration Interaction

In order to investigate charge resonance situations in molecular complexes, Wu et al. (J. Chem. Phys. 2007, 127, 164119) recently proposed a configuration interaction method with a valence bond-like multiconfigurational basis obtained from constrained DFT calculations. We adapt this method to the Self-Consistent Charge Density-Functional-based Tight Binding (SCC-DFTB) approach and provide expressions for the gradients of the energy with respect to the nuclear coordinates. It is shown that the method corrects the wrong SCC-DFTB behavior of the potential energy surface in the dissociation regions. This scheme is applied to determine the structural and stability properties of positively charged molecular dimers with full structural optimization, namely, the benzene dimer cation and the water dimer cation. The method yields binding energies in good agreement with experimental data and high-level reference calculations.

A theoretical study of linear beryllium chains: Full configuration interaction

We present a full configuration interaction study of Be(N) (N=2,3,4,5) linear chains. A comparative study of the basis-set effect on the reproduction of the energy profile has been reported. In particular, the 3s1p, 4s2p, 4s2p1d, 5s3p2d, and 5s3p2d1f bases were selected. For the smallest chains (i.e., Be(2) and Be(3)), smaller basis sets give dissociative energy profiles, so large basis set is demanded for the reproduction of equilibrium minima in the structures. For Be(4) and Be(5) linear chains, the energy profiles show a minimum also by using the smallest basis sets, but the largest ones give a much stronger stabilization energy. For all the structures, two spin states have been studied: the singlet and the triplet. It is shown that the energy separation of the two states, in the equilibrium region, is small and decays exponentially with respect to the number of atoms in the chain. Finally an interpolative technique allowing for the estimation of the long-chain parameters from shorter ones is presented.

Spin Density Distribution in Open-Shell Transition Metal Systems: A Comparative Post-Hartree-Fock, Density Functional Theory, and Quantum Monte Carlo Study of the CuCl2 Molecule.

We present a comparative study of the spatial distribution of the spin density of the ground state of CuCl2 using Density Functional Theory (DFT), quantum Monte Carlo (QMC), and post-Hartree-Fock wave function theory (WFT). A number of studies have shown that an accurate description of the electronic structure of the lowest-lying states of this molecule is particularly challenging due to the interplay between the strong dynamical correlation effects in the 3d shell and the delocalization of the 3d hole over the chlorine atoms. More generally, this problem is representative of the difficulties encountered when studying open-shell metal-containing molecular systems. Here, it is shown that qualitatively different results for the spin density distribution are obtained from the various quantum-mechanical approaches. At the DFT level, the spin density distribution is found to be very dependent on the functional employed. At the QMC level, Fixed-Node Diffusion Monte Carlo (FN-DMC) results are strongly dependent on the nodal structure of the trial wave function. Regarding wave function methods, most approaches not including a very high amount of dynamic correlation effects lead to a much too high localization of the spin density on the copper atom, in sharp contrast with DFT. To shed some light on these conflicting results Full CI-type (FCI) calculations using the 6-31G basis set and based on a selection process of the most important determinants, the so-called CIPSI approach (Configuration Interaction with Perturbative Selection done Iteratively) are performed. Quite remarkably, it is found that for this 63-electron molecule and a full CI space including about 10(18) determinants, the FCI limit can almost be reached. Putting all results together, a natural and coherent picture for the spin distribution is proposed.

Fixed-node diffusion Monte Carlo potential energy curve of the fluorine molecule F2 using selected configuration interaction trial wavefunctions

The potential energy curve of the F2 molecule is calculated with Fixed-Node Diffusion Monte Carlo (FN-DMC) using Configuration Interaction (CI)-type trial wavefunctions. To keep the number of determinants reasonable and thus make FN-DMC calculations feasible in practice, the CI expansion is restricted to those determinants that contribute the most to the total energy. The selection of the determinants is made using the CIPSI approach (Configuration Interaction using a Perturbative Selection made Iteratively). The trial wavefunction used in FN-DMC is directly issued from the deterministic CI program; no Jastrow factor is used and no preliminary multi-parameter stochastic optimization of the trial wavefunction is performed. The nodes of CIPSI wavefunctions are found to reduce significantly the fixed-node error and to be systematically improved upon increasing the number of selected determinants. To reduce the non-parallelism error of the potential energy curve, a scheme based on the use of a R-dependent number of determinants is introduced. Using Dunnings cc-pVDZ basis set, the FN-DMC energy curve of F2 is found to be of a quality similar to that obtained with full configuration interaction/cc-pVQZ.

Electron Pair Localization Function (EPLF) for Density Functional Theory and ab Initio Wave Function-Based Methods: A New Tool for Chemical Interpretation

We present a modified definition of the Electron Pair Localization Function (EPLF), initially defined within the framework of quantum Monte Carlo approaches [ Scemama , A. ; Caffarel , M. ; Chaquin , P. J. Chem. Phys. 2004 , 121 , 1725 ] to be used in Density Functional Theories (DFT) and ab initio wave-function-based methods. This modified version of the EPLF-while keeping the same physical and chemical contents-is built to be analytically computable with standard wave functions or Kohn-Sham representations. It is illustrated that the EPLF defines a simple and powerful tool for chemical interpretation via selected applications including atomic and molecular closed-shell systems, σ and π bonds, radical and singlet open-shell systems, and molecules having a strong multiconfigurational character. Some applications of the EPLF are presented at various levels of theory and compared to Becke and Edgecombes Electron Localization Function (ELF). Our open-source parallel software implementation of the EPLF opens the possibility of its use by a large community of chemists interested in the chemical interpretation of complex electronic structures.

An efficient sampling algorithm for variational Monte Carlo

We propose a new algorithm for sampling the N-body density mid R:Psi(R)mid R:(2)R(3N)mid R:Psimid R:(2) in the variational Monte Carlo framework. This algorithm is based upon a modified Ricci-Ciccotti discretization of the Langevin dynamics in the phase space (R,P) improved by a Metropolis-Hastings accept/reject step. We show through some representative numerical examples (lithium, fluorine, and copper atoms and phenol molecule) that this algorithm is superior to the standard sampling algorithm based on the biased random walk (importance sampling).

Accurate nonrelativistic ground-state energies of 3d transition metal atoms

We present accurate nonrelativistic ground-state energies of the transition metal atoms of the 3d series calculated with Fixed-Node Diffusion Monte Carlo (FN-DMC). Selected multi-determinantal expansions obtained with the CIPSI (Configuration Interaction using a Perturbative Selection made Iteratively) method and including the most prominent determinants of the full configuration interaction expansion are used as trial wavefunctions. Using a maximum of a few tens of thousands determinants, fixed-node errors on total DMC energies are found to be greatly reduced for some atoms with respect to those obtained with Hartree-Fock nodes. To the best of our knowledge, the FN-DMC/(CIPSI nodes) ground-state energies presented here are the lowest variational total energies reported so far. They differ from the recently recommended non-variational values of McCarthy and Thakkar [J. Chem. Phys. 136, 054107 (2012)] only by a few percents of the correlation energy. Thanks to the variational property of FN-DMC total energies, our results provide exact lower bounds for the absolute value of all-electron correlation energies, |Ec|.

Simple and efficient approach to the optimization of correlated wave functions

We present a simple and efficient method to optimize within energy minimization the determinantal component of the many-body wave functions commonly used in quantum Monte Carlo calculations. The approach obtains the optimal wave function as an approximate perturbative solution of an effective Hamiltonian iteratively constructed via Monte Carlo sampling. The effectiveness of the method, as well as its ability to substantially improve the accuracy of quantum Monte Carlo calculations, are demonstrated by optimizing a large number of parameters for the ground state of acetone and the difficult case of the

Quantum Monte Carlo for large chemical systems: Implementing efficient strategies for petascale platforms and beyond

1^{1}B_{u}

Code interoperability and standard data formats in quantum chemistry and quantum dynamics: The Q5/D5Cost data model.

state of hexatriene.