1,720,981 research outputs found
Ground State Geometries of Polyacetylene Chains from Many-Particle Quantum Mechanics
Full text linkDue to the crucial role played by electron correlation, the accurate determination of ground state geometries of π-conjugated molecules is still a challenge for many quantum chemistry methods. Because of the high parallelism of the algorithms and their explicit treatment of electron correlation effects, Quantum Monte Carlo calculations can offer an accurate and reliable description of the electronic states and of the geometries of such systems, competing with traditional quantum chemistry approaches. Here, we report the structural properties of polyacetylene chains H-(C2H2)N-H up to N = 12 acetylene units, by means of Variational Monte Carlo (VMC) calculations based on the multi-determinant Jastrow Antisymmetrized Geminal Power (JAGP) wave function. This compact ansatz can provide for such systems an accurate description of the dynamical electronic correlation as recently detailed for the 1,3-butadiene molecule [J. Chem. Theory Comput. 2015 11 (2), 508-517]. The calculated Bond Length Alternation (BLA), namely the difference between the single and double carbon bonds, extrapolates, for N → ∞, to a value of 0.0910(7) Å, compatible with the experimental data. An accurate analysis was able to distinguish between the influence of the multi-determinantal AGP expansion and of the Jastrow factor on the geometrical properties of the fragments. Our size-extensive and self-interaction-free results provide new and accurate ab initio references for the structures of the ground state of polyenes. © 2015 American Chemical Society
π-Conjugation in trans -1,3-Butadiene: Static and dynamical electronic correlations described through quantum monte carlo
No full textWe investigate the effects of the static and dynamical electronic correlations on the level of conjugation of the trans-1,3-butadiene molecule through Quantum Monte Carlo methods applied to an Antisymmetrized Geminal Power (AGP) wave function, with a Jastrow factor similar to the Gutzwiller ansatz. The degree of conjugation is measured through the convergence of the structural properties of 1,3-butadiene and in particular of the Bond Length Alternation (BLA), that is the difference between the lengths of the single and double carbon bonds. After verifying the different roles of the Fermionic AGP part of our wave function and of the Jastrow factor in recovering electronic correlation, we study the effects of a constrained Active Space AGP (AGPAS), similar to that used in the Complete Active Space (CAS) representation. Through this AGPAS, we are able to identify the effect of the limited active space on the degree of conjugation, showing that in the limit of infinite active space the structural properties converge exactly to those of the atomic AGP, giving a BLA for 1,3-butadiene around 0.1244(5) Å. © 2015 American Chemical Society
Geometries of low spin states of multi-centre transition metal complexes through extended broken symmetry variational Monte Carlo
No full textThe correct description of the ground state electronic and geometrical properties of multi-centre transition metal complexes necessitates of a high-level description of both dynamical and static correlation effects. In di-metallic complexes, the ground state low spin properties can be computed starting from single-determinants High-Spin (HS) and Broken Symmetry (BS) states by reconstructing an approximated low spin potential energy surface through the extended broken symmetry approach, based on the Heisenberg Hamiltonian. In the present work, we first apply this approach within the variational Monte Carlo method to tackle the geometry optimization of a Fe2S2(SH)4(2-) model complex. To describe the HS and BS wavefunctions, we use a fully optimized unrestricted single determinant with a correlated Jastrow factor able to recover a large amount of dynamical correlation. We compared our results with those obtained by density functional theory and other multiconfigurational approaches, discussing the role of the nodal surface on the structural parameters
Kohn-Sham orbitals and potentials from quantum Monte Carlo molecular densities
No full textIn this work we show the possibility to extract Kohn-Sham orbitals, orbital energies, and exchange correlation potentials from accurate Quantum Monte Carlo (QMC) densities for atoms (He, Be, Ne) and molecules (H-2, Be-2, H2O, and C2H4). The Variational Monte Carlo (VMC) densities based on accurate Jastrow Antisymmetrised Geminal Power wave functions are calculated through different estimators. Using these reference densities, we extract the Kohn-Sham quantities with the method developed by Zhao, Morrison, and Parr (ZMP) [Phys. Rev. A 50, 2138 (1994)]. We compare these extracted quantities with those obtained form CISD densities and with other data reported in the literature, finding a good agreement between VMC and other high-level quantum chemistry methods. Our results demonstrate the applicability of the ZMP procedure to QMC molecular densities, that can be used for the testing and development of improved functionals and for the implementation of embedding schemes based on QMC and Density Functional Theory. (C) 2014 AIP Publishing LLC
Structural optimization by quantum Monte Carlo: Investigating the low-lying excited states of ethylene
No full textWe present full structural optimizations of the ground state and of the low lying triplet state of the ethylene molecule by means of Quantum Monte Carlo methods. Using the efficient structural optimization method based on renormalization techniques and on adjoint differentiation algorithms recently proposed [Sore S.; Capriotti, L. J. Chem. Phys. 2010, 133, 234111], we present the variational convergence of both wave function parameters and atomic positions. All of the calculations were done using an accurate and compact wave function based on Pauling's resonating valence bond representation: the Jastrow Antisymmetrized Geminal Power (JAGP). All structural and wave function parameters are optimized, including coefficients and exponents of the Gaussian primitives of the AGP and the Jastrow atomic orbitals. Bond lengths and bond angles are calculated with a statistical error of about 0.1% and are in good agreement with the available experimental data. The Variational and Diffusion Monte Carlo calculations estimate vertical and adiabatic excitation energies in the ranges 4.623(10)-4.688(5) eV and 3.001(5)-3.091(5) eV, respectively. The adiabatic gap, which is in line with other correlated quantum chemistry methods, is slightly higher than the value estimated by recent photodissociation experiments. Our results demonstrate how Quantum Monte Carlo calculations have become a promising and computationally affordable tool for the structural optimization of correlated molecular systems
Role of Electron Correlation along the Water Splitting Reaction
Full text linkElectron correlation plays a crucial role in the energetics of reactions catalyzed by transition metal complexes, such as water splitting. In the present work we exploit the performance of various methods to describe the thermodynamics of a simple but representative model of water splitting reaction, based on a single cobalt ion as catalyst. Density Functional Theory (DFT) calculations show a significant dependence on the adopted functional, and not negligible differences with respect to CCSD(T) findings are found along the reaction cycle. We performed quantum Monte Carlo calculations using an unrestricted single Slater determinant wave function multiplied by a Jastrow factor using both DFT and fully optimized orbitals. Variational and Lattice Regularized Diffusion Monte Carlo results are in overall agreement with the CCSD(T) free-energy profile, even though differences in the description of the thermodynamics of the reaction cycle are found. NEVPT2 calculations reveal that the role of the static correlation of the different reaction steps is not large, and it is limited to only a few intermediate structures. Finally, the free-energy difference of the overall water splitting reaction computed at the quantum Monte Carlo level shows an excellent match with the experimental value of 4.92 eV, underlining the capability of these techniques to properly describe the dynamical correlation of such reactions
Molecular Electrical Properties from Quantum Monte Carlo Calculations: Application to Ethyne
No full text"We used Quantum Monte Carlo (QMC) methods to study the polarizability and the quadrupole moment of the ethyne molecule using the Jastrow-Antisymmetrised Geminal Power (JAGP) wave function, a compact and strongly correlated variational ansatz. The compactness of the functional form and the full optimization of all its variational parameters, including linear and exponential coefficients in atomic orbitals, allow us to observe a fast convergence of the electrical properties with the size of the atomic and Jastrow basis sets. Both variational results on isotropic polarizability and quadrupole moment based on Gaussian type and Slater type basis sets are very close to the Lattice Regularized Diffusion Monte Carlo values and in very good agreement with experimental data and with other quantum chemistry calculations. We also study the electronic density along the C C and C-H bonds by introducing a generalization for molecular systems of the small-variance improved estimator of the electronic density proposed by Assaraf et al. (Assaraf, R; Caffarel, M.; Scemama, A. Phys. Rev. E, 2007, 75, 035701).
Angle-resolved photoemission spectroscopy from first-principles quantum Monte Carlo
Full text linkAngle-resolved photoemission spectroscopy allows one to visualize in momentum space the probability weight maps of electrons subtracted from molecules deposited on a substrate. The interpretation of these maps usually relies on the plane wave approximation through the Fourier transform of single particle orbitals obtained from density functional theory. Here we propose a first-principle many-body approach based on quantum Monte Carlo (QMC) to directly calculate the quasi-particle wave functions (also known as Dyson orbitals) of molecules in momentum space. The comparison between these correlated QMC images and their single particle counterpart highlights features that arise from many-body effects. We test the QMC approach on the linear C2H2, CO2, and N2 molecules, for which only small amplitude remodulations are visible. Then, we consider the case of the pentacene molecule, focusing on the relationship between the momentum space features and the real space quasi-particle orbital. Eventually, we verify the correlation effects present in the metal CuCl42- planar complex
Correlation Effects in Scanning Tunneling Microscopy Images of Molecules Revealed by Quantum Monte Carlo
No full textScanning tunneling microscopy (STM) and spectroscopy probe the local density of states of single molecules electrically insulated from the substrate. The experimental images, although usually interpreted in terms of single-particle molecular orbitals, are associated with quasiparticle wave functions dressed by the whole electron-electron interaction. Here we propose an ab initio approach based on quantum Monte Carlo to calculate the quasiparticle wave functions of molecules. Through the comparison between Monte Carlo wave functions and their uncorrelated Hartree-Fock counterparts we visualize the electronic correlation embedded in the simulated STM images, highlighting the many-body features that might be observed
Role of Electron Correlation along the Water Splitting Reaction
No full textElectron correlation
plays a crucial role in the energetics of
reactions catalyzed by transition metal complexes, such as water splitting.
In the present work we exploit the performance of various methods
to describe the thermodynamics of a simple but representative model
of water splitting reaction, based on a single cobalt ion as catalyst.
Density Functional Theory (DFT) calculations show a significant dependence
on the adopted functional, and not negligible differences with respect
to CCSD(T) findings are found along the reaction cycle. We performed
quantum Monte Carlo calculations using an unrestricted single Slater
determinant wave function multiplied by a Jastrow factor using both
DFT and fully optimized orbitals. Variational and Lattice Regularized
Diffusion Monte Carlo results are in overall agreement with the CCSD(T)
free-energy profile, even though differences in the description of
the thermodynamics of the reaction cycle are found. NEVPT2 calculations
reveal that the role of the static correlation of the different reaction
steps is not large, and it is limited to only a few intermediate structures.
Finally, the free-energy difference of the overall water splitting
reaction computed at the quantum Monte Carlo level shows an excellent
match with the experimental value of 4.92 eV, underlining the capability
of these techniques to properly describe the dynamical correlation
of such reactions
- …
