1
|
Otis L, Neuscamman E. Optimization Stability in Excited-State-Specific Variational Monte Carlo. J Chem Theory Comput 2023; 19:767-782. [PMID: 36662538 DOI: 10.1021/acs.jctc.2c00642] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Abstract
We investigate the issue of optimization stability in variance-based state-specific variational Monte Carlo, discussing the roles of the objective function, the complexity of wave function ansatz, the amount of sampling effort, and the choice of minimization algorithm. Using a small cyanine dye molecule as a test case, we systematically perform minimizations using variants of the linear method as both a standalone algorithm and in a hybrid combination with accelerated descent. We demonstrate that adaptive step control is crucial for maintaining the linear method's stability when optimizing complicated wave functions and that the hybrid method enjoys both greater stability and minimization performance.
Collapse
Affiliation(s)
- Leon Otis
- Department of Physics, University of California Berkeley, Berkeley, California 94720, United States
| | - Eric Neuscamman
- Department of Chemistry, University of California Berkeley, Berkeley, California 94720, United States.,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
2
|
Otis L, Neuscamman E. A promising intersection of excited‐state‐specific methods from quantum chemistry and quantum Monte Carlo. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2023. [DOI: 10.1002/wcms.1659] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/11/2023]
Affiliation(s)
- Leon Otis
- Department of Physics University of California Berkeley Berkeley California USA
| | - Eric Neuscamman
- Department of Chemistry University of California Berkeley Berkeley California USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory Berkeley California USA
| |
Collapse
|
3
|
Chen S, Zhang S. A structural optimization algorithm with stochastic forces and stresses. NATURE COMPUTATIONAL SCIENCE 2022; 2:736-744. [PMID: 38177372 DOI: 10.1038/s43588-022-00350-w] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/07/2022] [Accepted: 10/07/2022] [Indexed: 01/06/2024]
Abstract
We propose an algorithm for optimizations in which the gradients contain stochastic noise. This arises, for example, in structural optimizations when computations of forces and stresses rely on methods involving Monte Carlo sampling, such as quantum Monte Carlo or neural network states, or are performed on quantum devices that have intrinsic noise. Our proposed algorithm is based on the combination of two ingredients: an update rule derived from the steepest-descent method, and a staged scheduling of the targeted statistical error and step size, with position averaging. We compare it with commonly applied algorithms, including some of the latest machine learning optimization methods, and show that the algorithm consistently performs efficiently and robustly under realistic conditions. Applying this algorithm, we achieve full-degree optimizations in solids using ab initio many-body computations, by auxiliary-field quantum Monte Carlo with plane waves and pseudopotentials. A potential metastable structure in Si is discovered using density-functional calculations with synthetic noisy forces.
Collapse
Affiliation(s)
- Siyuan Chen
- Department of Physics, College of William & Mary, Williamsburg, VA, USA.
| | - Shiwei Zhang
- Center for Computational Quantum Physics, Flatiron Institute, New York, NY, USA
| |
Collapse
|
4
|
Iyer GR, Rubenstein BM. Finite-Size Error Cancellation in Diffusion Monte Carlo Calculations of Surface Chemistry. J Phys Chem A 2022; 126:4636-4646. [PMID: 35820033 DOI: 10.1021/acs.jpca.2c01957] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The accurate prediction of reaction mechanisms in heterogeneous (surface) catalysis is one of the central challenges in computational chemistry. Quantum Monte Carlo methods─Diffusion Monte Carlo (DMC) in particular─are being recognized as higher-accuracy, albeit more computationally expensive, alternatives to Density Functional Theory (DFT) for energy predictions of catalytic systems. A major computational bottleneck in the broader adoption of DMC for catalysis is the need to perform finite-size extrapolations by simulating increasingly large periodic cells (supercells) to eliminate many-body finite-size effects and obtain energies in the thermodynamic limit. Here, we show that it is possible to significantly reduce this computational cost by leveraging the cancellation of many-body finite-size errors that accompanies the evaluation of energy differences when calculating quantities like adsorption (binding) energies and mapping potential energy surfaces. We analyze the cancellation and convergence of many-body finite-size errors in two well-known adsorbate/slab systems, H2O/LiH(001) and CO/Pt(111). Based on this analysis, we identify strategies for obtaining binding energies in the thermodynamic limit that optimally utilize error cancellation to balance accuracy and computational efficiency. Using one such strategy, we then predict the correct order of adsorption site preference on CO/Pt(111), a challenging problem for a wide range of density functionals. Our accurate and inexpensive DMC calculations are found to unambiguously recover the top > bridge > hollow site order, in agreement with experimental observations. We proceed to use this DMC method to map the potential energy surface of CO hopping between Pt(111) adsorption sites. This reveals the existence of an L-shaped top-bridge-hollow diffusion trajectory characterized by energy barriers that provide an additional kinetic justification for experimental observations of CO/Pt(111) adsorption. Overall, this work demonstrates that it is routinely possible to achieve order-of-magnitude speedups and memory savings in DMC calculations by taking advantage of error cancellation in the calculation of energy differences that are ubiquitous in heterogeneous catalysis and surface chemistry more broadly.
Collapse
Affiliation(s)
- Gopal R Iyer
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, United States
| | - Brenda M Rubenstein
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, United States
| |
Collapse
|
5
|
Smith JET, Lee J, Sharma S. Near-Exact Nuclear Gradients of Complete Active Space Self-Consistent Field Wave Functions. J Chem Phys 2022; 157:094104. [DOI: 10.1063/5.0085515] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
In this paper, we study the nuclear gradients of heat bath configuration interaction self-consistent field (HCISCF) wave functions and use them to optimize molecular geometries for various molecules.We show that the HCISCF nuclear gradients are fairly insensitive to the size of the "selected" variational space, which allows us to reduce the computational cost without introducing significant error.The ability of HCISCF to treat larger active spaces combined with the flexibility for users to control the computational cost makes the method very attractive for studying strongly correlated systems which require a larger active space than possible with complete active space self-consistent field (CASSCF).Finally, we study the realistic catalyst, Fe(PDI), and highlight some of the challenges this system poses for density functional theory (DFT).We demonstrate how HCISCF can clarify the energetic stability of geometries obtained from DFT when the results are strongly dependent on the functional.We also use the HCISCF gradients to optimize geometries for this species and study the adiabatic singlet-triplet gap. During geometry optimization, we find that multiple near-degenerate local minima exist on the triplet potential energy surface.
Collapse
Affiliation(s)
- James E. T. Smith
- Center for Computational Quantum Phyics, Flatiron Institute, United States of America
| | - Joonho Lee
- Chemistry, Columbia University Department of Chemistry, United States of America
| | - Sandeep Sharma
- University of Colorado at Boulder, United States of America
| |
Collapse
|
6
|
Dey G, Chakraborty A. Study of the conformations and tautomerisation pathway in (Z)-4-(hydroxypropyl) isochroman-1, 3‑dione: Analysis through energy, vibrational signatures and hardness profiles. J Mol Struct 2022. [DOI: 10.1016/j.molstruc.2021.132177] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
7
|
Tiihonen J, Kent PRC, Krogel JT. Surrogate Hessian accelerated structural optimization for stochastic electronic structure theories. J Chem Phys 2022; 156:054104. [DOI: 10.1063/5.0079046] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Juha Tiihonen
- Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - Paul R. C. Kent
- Computational Sciences and Engineering Division, Oak Ridge National Laboratory, Tennessee 37831, USA
| | - Jaron T. Krogel
- Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| |
Collapse
|
8
|
Chakraborty A, Dey G. Conformations and tautomerisation between (Z)-4-(hydroxyethyl) isochroman-1, 3-dione and and 4-acetyl-3-hydroxyisochroman-1-one: A computational study through Energy, electron Distribution, vibrational analysis and hardness profiles. COMPUT THEOR CHEM 2021. [DOI: 10.1016/j.comptc.2021.113486] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
|
9
|
Dey G, Chakraborty A. Conformational Landscape and Tautomerisation in (Z)-4-(hydroxymethylene) isochroman-1,3-dione: Analysis through Energy and Hardness profiles. J Mol Struct 2021. [DOI: 10.1016/j.molstruc.2021.130859] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
10
|
Wang J, Durbeej B. Thermal Fluctuations in Conjugation and their Effect on Calculated Excitation Energies: A Case Study on the Astaxanthin Carotenoid. CHEMPHOTOCHEM 2021. [DOI: 10.1002/cptc.202100178] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Jun Wang
- Jiangsu Key Laboratory for Chemistry of Low-Dimensional Materials Jiangsu Engineering Laboratory for Environment Functional Materials School of Chemistry and Chemical Engineering Huaiyin Normal University No. 111 West Changjiang Road 223300 Huaian Jiangsu Province China
| | - Bo Durbeej
- Division of Theoretical Chemistry IFM Linköping University 581 83 Linköping Sweden
| |
Collapse
|
11
|
Electronically excited state structures and stabilities of organic small molecules: A DFT study of triphenylamine derivatives. Chem Phys 2021. [DOI: 10.1016/j.chemphys.2021.111256] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
|
12
|
Chang JL, Jen CC, Huang YJ, Du ZK. Distinction of photoelectron spectroscopy of cis- and trans-acrolein explored by theoretical computation. Chem Phys 2021. [DOI: 10.1016/j.chemphys.2021.111166] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
13
|
Çaylak O, Baumeier B. Excited-State Geometry Optimization of Small Molecules with Many-Body Green's Functions Theory. J Chem Theory Comput 2021; 17:879-888. [PMID: 33399447 PMCID: PMC7876808 DOI: 10.1021/acs.jctc.0c01099] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
We present a benchmark study of gas
phase geometry optimizations
in the excited states of carbon monoxide, acetone, acrolein, and methylenecyclopropene
using many-body Green’s functions theory within the GW approximation and the Bethe–Salpeter equation
(BSE) employing numerical gradients. We scrutinize the influence of
several typical approximations in the GW-BSE framework;
we used one-shot G0W0 or eigenvalue self-consistent evGW, employing
a fully analytic approach or plasmon-pole model for the frequency
dependence of the electron self-energy, or performing the BSE step
within the Tamm–Dancoff approximation. The obtained geometries
are compared to reference results from multireference perturbation
theory (CASPT2), variational Monte Carlo (VMC) method, second-order
approximate coupled cluster (CC2) method, and time-dependent density-functional
theory (TDDFT). We find overall a good agreement of the structural
parameters optimized with the GW-BSE calculations
with CASPT2, with an average relative error of around 1% for the G0W0 and 1.5% for
the evGW variants based on a PBE0 ground state, respectively,
while the other approximations have negligible influence. The relative
errors are also smaller than those for CC2 and TDDFT with different
functionals and only larger than VMC, indicating that the GW-BSE method does not only yield excitation energies but
also geometries in good agreement with established higher-order wave
function methods.
Collapse
Affiliation(s)
- Onur Çaylak
- Department of Mathematics and Computer Science, Eindhoven University of Technology, P.O. Box 513, 5600MB Eindhoven, The Netherlands.,Institute for Complex Molecular Systems, Eindhoven University of Technology, P.O. Box 513, 5600MB Eindhoven, The Netherlands
| | - Björn Baumeier
- Department of Mathematics and Computer Science, Eindhoven University of Technology, P.O. Box 513, 5600MB Eindhoven, The Netherlands.,Institute for Complex Molecular Systems, Eindhoven University of Technology, P.O. Box 513, 5600MB Eindhoven, The Netherlands
| |
Collapse
|
14
|
Carvalho CMC, Gargano R, Martins JBL, Politi JRS. Accurate spectroscopic properties by diffusion quantum Monte Carlo calculations. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2020; 243:118707. [PMID: 32827906 DOI: 10.1016/j.saa.2020.118707] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/28/2020] [Revised: 07/05/2020] [Accepted: 07/13/2020] [Indexed: 06/11/2023]
Abstract
The capability of Diffusion Quantum Monte Carlo (DMC) to produce high quality potential energy curve (PEC) was evaluated. H2+, HeH+ and LiH PECs were built by all-electron fixed-node DMC calculations. Trial wave functions were obtained from Hartree-Fock (HF) (H2+), MCSCF and CI (HeH+ and LiH) calculations multiplied by Jastrow factor. The quality of these generated PECs was analyzed throughout equilibrium distance, dissociation energy, vibrational energies and rovibrational spectroscopic constants (ωe, ωexe, ωeye, αe, γe and Be). The Discrete Variable Representation (DVR) and the Dunham approaches were used to determine the rovibrational spectroscopic constants. The PECs and the aforementioned properties were also obtained by the following methods: MCSCF/aug-cc-pV5Z (LiH), CCSD(T)/aug-cc-pV5Z (HeH+ and LiH) and HF (H2+ and HeH+) levels. The results of these DMC computations, specially the DMC-DVR procedure, are the most accurate among others DMC calculations available in the literature for these systems. They suggest that DMC can be used to achieve accurate PECs to produce spectroscopic properties with the same level of accuracy of theoretical benchmarks and experimental data of the literature.
Collapse
Affiliation(s)
- Cassius M C Carvalho
- Institute of Chemistry, University of Brasília, Campus Darcy Ribeiro, Brasília, DF, Brazil
| | - Ricardo Gargano
- Institute of Physics, University of Brasília, Campus Darcy Ribeiro, Brasília, DF, Brazil
| | - João B L Martins
- Institute of Chemistry, University of Brasília, Campus Darcy Ribeiro, Brasília, DF, Brazil
| | - José Roberto S Politi
- Institute of Chemistry, University of Brasília, Campus Darcy Ribeiro, Brasília, DF, Brazil.
| |
Collapse
|
15
|
Chrayteh A, Blondel A, Loos PF, Jacquemin D. Mountaineering Strategy to Excited States: Highly Accurate Oscillator Strengths and Dipole Moments of Small Molecules. J Chem Theory Comput 2020; 17:416-438. [DOI: 10.1021/acs.jctc.0c01111] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Affiliation(s)
- Amara Chrayteh
- Université de Nantes, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| | - Aymeric Blondel
- Université de Nantes, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques, Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Denis Jacquemin
- Université de Nantes, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| |
Collapse
|
16
|
Hill NS, Coote ML. Rational design of photo-cleavable alkoxyamines for polymerization and synthesis. Phys Chem Chem Phys 2020; 22:19680-19686. [PMID: 32830210 DOI: 10.1039/d0cp02924f] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Theoretical calculations have been performed in order to investigate the impact of different substitution patterns on predicted photoreactivity of alkoxyamines fused to an anthraquinone chromophore. Amino and hydroxy groups (similar to those which have been previously synthesized) are introduced and their effect on excited state energies and charge transfer is assessed. Analogous to formally oxidized alkoxyamines, the charge-separated nNπ* state can undergo mesolytic cleavage or bimolecular or SN2 reactions with nucleophiles, according to the substitution patterns and other reagents present. While homolytic cleavage is in principle promoted by triplet ππ* states, the accessible ππ* triplet states in this system are centered on the chromophore and unreactive. We show that the reactive nNπ* state, which bears a negative charge, is stabilized by hydroxy substitution while amino substitution will destabilize it. After mesolysis to a carbon centred radical, the nitroxide radical re-forms; however, when carbocations are produced the remaining open-shell singlet is stable and unable to undergo coupling with the carbocation.
Collapse
Affiliation(s)
- Nicholas S Hill
- ARC Centre of Excellence for Electromaterials Science, Research School of Chemistry, Australian National University, Canberra, ACT 2601, Australia.
| | - Michelle L Coote
- ARC Centre of Excellence for Electromaterials Science, Research School of Chemistry, Australian National University, Canberra, ACT 2601, Australia.
| |
Collapse
|
17
|
Grotjahn R, Kaupp M. Validation of Local Hybrid Functionals for Excited States: Structures, Fluorescence, Phosphorescence, and Vibronic Spectra. J Chem Theory Comput 2020; 16:5821-5834. [DOI: 10.1021/acs.jctc.0c00520] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Robin Grotjahn
- Technische Universität Berlin, Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Straße des 17. Juni 135, D-10623 Berlin, Germany
| | - Martin Kaupp
- Technische Universität Berlin, Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Straße des 17. Juni 135, D-10623 Berlin, Germany
| |
Collapse
|
18
|
Wang J, Durbeej B. How accurate are TD‐DFT excited‐state geometries compared to DFT ground‐state geometries? J Comput Chem 2020; 41:1718-1729. [DOI: 10.1002/jcc.26213] [Citation(s) in RCA: 31] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2020] [Revised: 04/07/2020] [Accepted: 04/09/2020] [Indexed: 01/18/2023]
Affiliation(s)
- Jun Wang
- Division of Theoretical Chemistry, IFMLinköping University Linköping Sweden
- Institut de Química Computacional i Catàlisi, Facultat de CiènciesUniversitat de Girona Girona Spain
| | - Bo Durbeej
- Division of Theoretical Chemistry, IFMLinköping University Linköping Sweden
| |
Collapse
|
19
|
Loos PF, Scemama A, Jacquemin D. The Quest for Highly Accurate Excitation Energies: A Computational Perspective. J Phys Chem Lett 2020; 11:2374-2383. [PMID: 32125872 DOI: 10.1021/acs.jpclett.0c00014] [Citation(s) in RCA: 87] [Impact Index Per Article: 21.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/11/2023]
Abstract
We provide an overview of the successive steps that made it possible to obtain increasingly accurate excitation energies with computational chemistry tools, eventually leading to chemically accurate vertical transition energies for small- and medium-size molecules. First, we describe the evolution of ab initio methods employed to define benchmark values, with the original Roos CASPT2 method, then the CC3 method as in the renowned Thiel set, and more recently the resurgence of selected configuration interaction methods. The latter method has been able to deliver consistently, for both single and double excitations, highly accurate excitation energies for small molecules, as well as medium-size molecules with compact basis sets. Second, we describe how these high-level methods and the creation of representative benchmark sets of excitation energies have allowed the fair and accurate assessment of the performance of computationally lighter methods. We conclude by discussing possible future theoretical and technological developments in the field.
Collapse
Affiliation(s)
- Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques, Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Anthony Scemama
- Laboratoire de Chimie et Physique Quantiques, Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Denis Jacquemin
- Université de Nantes, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| |
Collapse
|
20
|
Loos PF, Lipparini F, Boggio-Pasqua M, Scemama A, Jacquemin D. A Mountaineering Strategy to Excited States: Highly Accurate Energies and Benchmarks for Medium Sized Molecules. J Chem Theory Comput 2020; 16:1711-1741. [PMID: 31986042 DOI: 10.1021/acs.jctc.9b01216] [Citation(s) in RCA: 107] [Impact Index Per Article: 26.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Following our previous work focusing on compounds containing up to 3 non-hydrogen atoms [J. Chem. Theory Comput. 2018, 14, 4360-4379], we present here highly accurate vertical transition energies obtained for 27 molecules encompassing 4, 5, and 6 non-hydrogen atoms: acetone, acrolein, benzene, butadiene, cyanoacetylene, cyanoformaldehyde, cyanogen, cyclopentadiene, cyclopropenone, cyclopropenethione, diacetylene, furan, glyoxal, imidazole, isobutene, methylenecyclopropene, propynal, pyrazine, pyridazine, pyridine, pyrimidine, pyrrole, tetrazine, thioacetone, thiophene, thiopropynal, and triazine. To obtain these energies, we use equation-of-motion/linear-response coupled cluster theory up to the highest technically possible excitation order for these systems (CC3, EOM-CCSDT, and EOM-CCSDTQ) and selected configuration interaction (SCI) calculations (with tens of millions of determinants in the reference space), as well as the multiconfigurational n-electron valence state perturbation theory (NEVPT2) method. All these approaches are applied in combination with diffuse-containing atomic basis sets. For all transitions, we report at least CC3/aug-cc-pVQZ vertical excitation energies as well as CC3/aug-cc-pVTZ oscillator strengths for each dipole-allowed transition. We show that CC3 almost systematically delivers transition energies in agreement with higher-level methods with a typical deviation of ±0.04 eV, except for transitions with a dominant double excitation character where the error is much larger. The present contribution gathers a large, diverse, and accurate set of more than 200 highly accurate transition energies for states of various natures (valence, Rydberg, singlet, triplet, n → π*, π → π*, ...). We use this series of theoretical best estimates to benchmark a series of popular methods for excited state calculations: CIS(D), ADC(2), CC2, STEOM-CCSD, EOM-CCSD, CCSDR(3), CCSDT-3, CC3, and NEVPT2. The results of these benchmarks are compared to the available literature data.
Collapse
Affiliation(s)
- Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques, CNRS et Université Toulouse III - Paul Sabatier, 118 route de Narbonne, 31062 Toulouse, France
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, Via Moruzzi 3, 56124 Pisa, Italy
| | - Martial Boggio-Pasqua
- Laboratoire de Chimie et Physique Quantiques, CNRS et Université Toulouse III - Paul Sabatier, 118 route de Narbonne, 31062 Toulouse, France
| | - Anthony Scemama
- Laboratoire de Chimie et Physique Quantiques, CNRS et Université Toulouse III - Paul Sabatier, 118 route de Narbonne, 31062 Toulouse, France
| | - Denis Jacquemin
- CEISAM Lab, UMR 6230, Université de Nantes, CNRS, F-44000 Nantes, France
| |
Collapse
|
21
|
Fihey A, Jacquemin D. Performances of Density Functional Tight-Binding Methods for Describing Ground and Excited State Geometries of Organic Molecules. J Chem Theory Comput 2019; 15:6267-6276. [DOI: 10.1021/acs.jctc.9b00688] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Affiliation(s)
- Arnaud Fihey
- Institut des Sciences Chimiques de Rennes, UMR 6226 CNRS, Université de Rennes 1, 263 Av. du Général Leclerc, 35042 Cedex Rennes, France
| | - Denis Jacquemin
- Laboratoire CEISAM - UMR CNRS 6230, Université de Nantes, 2 Rue de la Houssinière, BP 92208, 44322 Cedex 3 Nantes, France
| |
Collapse
|
22
|
Zhou X, Wang F. Singlet-triplet gaps in diradicals obtained with diffusion quantum Monte Carlo using a Slater-Jastrow trial wavefunction with a minimum number of determinants. Phys Chem Chem Phys 2019; 21:20422-20431. [PMID: 31501831 DOI: 10.1039/c9cp03045j] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Diradicals are essential species in a wide range of chemical processes, whereas the computational study of their electronic structure often remains a challenge due to near-degeneracy of the frontier molecular orbitals. The fixed-node diffusion quantum Monte Carlo (FN-DMC) method is employed to calculate adiabatic energy gaps of some typical diradicals with the Slater-Jastrow trial wavefunction. The antisymmetrized part of the trial wavefunction is taken to be a linear combination of a minimum number of determinants using RB3LYP orbitals from the closed-shell singlet state or ROB3LYP orbitals from the triplet state. Our results show that using the two-determinant-Jastrow trial wavefunction is necessary to achieve reliable energy differences between closed-shell singlet states. The energy of the triplet state with MS = 1 is calculated to be lower than that with MS = 0 with FN-DMC even using trial wavefunctions with spin-pure states as their antisymmetrized parts and this difference is reduced with better orbitals. This indicates that the fixed-node error is smaller for the triplet state with MS = 1. Adiabatic energy gaps obtained from the present FN-DMC calculations are in reasonable agreement with available experimental values. Compared with results of the high level EOM-SF-CC method, energy gaps of FN-DMC with RB3LYP orbitals are slightly better than those using ROB3LYP orbitals and results of EOM-SF-CCSD. The present FN-DMC calculations using the simplest ansatz for the trial wavefunction can achieve reasonable results for these diradicals and they can readily be applied to large diradicals.
Collapse
Affiliation(s)
- Xiaojun Zhou
- Institute of Atomic and Molecular Physics, Key Laboratory of High Energy Density Physics and Technology, Ministry of Education, Sichuan University, Chengdu, P. R. China.
| | | |
Collapse
|
23
|
Fang C, Durbeej B. Calculation of Free-Energy Barriers with TD-DFT: A Case Study on Excited-State Proton Transfer in Indigo. J Phys Chem A 2019; 123:8485-8495. [DOI: 10.1021/acs.jpca.9b05163] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Changfeng Fang
- Center for Optics Research and Engineering (CORE), Shandong University, Qingdao 266237, China
| | - Bo Durbeej
- Division of Theoretical Chemistry, IFM, Linköping University, SE-581 83 Linköping, Sweden
| |
Collapse
|
24
|
Dash M, Feldt J, Moroni S, Scemama A, Filippi C. Excited States with Selected Configuration Interaction-Quantum Monte Carlo: Chemically Accurate Excitation Energies and Geometries. J Chem Theory Comput 2019; 15:4896-4906. [PMID: 31348645 PMCID: PMC6740157 DOI: 10.1021/acs.jctc.9b00476] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We employ quantum Monte Carlo to obtain chemically accurate vertical and adiabatic excitation energies, and equilibrium excited-state structures for the small, yet challenging, formaldehyde and thioformaldehyde molecules. A key ingredient is a robust protocol to obtain balanced ground- and excited-state Jastrow-Slater wave functions at a given geometry, and to maintain such a balanced description as we relax the structure in the excited state. We use determinantal components generated via a selected configuration interaction scheme which targets the same second-order perturbation energy correction for all states of interest at different geometries, and fully optimize all variational parameters in the resultant Jastrow-Slater wave functions. Importantly, the excitation energies as well as the structural parameters in the ground and excited states are converged with very compact wave functions comprising few thousand determinants in a minimally augmented double-ζ basis set. These results are obtained already at the variational Monte Carlo level, the more accurate diffusion Monte Carlo method yielding only a small improvement in the adiabatic excitation energies. We find that matching Jastrow-Slater wave functions with similar variances can yield excitation energies compatible with our best estimates; however, the variance-matching procedure requires somewhat larger determinantal expansions to achieve the same accuracy, and it is less straightforward to adapt during structural optimization in the excited state.
Collapse
Affiliation(s)
- Monika Dash
- MESA+ Institute for Nanotechnology , University of Twente , P.O. Box 217, 7500 AE Enschede , The Netherlands
| | - Jonas Feldt
- MESA+ Institute for Nanotechnology , University of Twente , P.O. Box 217, 7500 AE Enschede , The Netherlands
| | - Saverio Moroni
- CNR-IOM DEMOCRITOS , Istituto Officina dei Materiali and SISSA Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , I-34136 Trieste , Italy
| | - Anthony Scemama
- Laboratoire de Chimie et Physique Quantiques , Université de Toulouse, CNRS, UPS , 31062 Toulouse , France
| | - Claudia Filippi
- MESA+ Institute for Nanotechnology , University of Twente , P.O. Box 217, 7500 AE Enschede , The Netherlands
| |
Collapse
|
25
|
Loos PF, Boggio-Pasqua M, Scemama A, Caffarel M, Jacquemin D. Reference Energies for Double Excitations. J Chem Theory Comput 2019; 15:1939-1956. [DOI: 10.1021/acs.jctc.8b01205] [Citation(s) in RCA: 79] [Impact Index Per Article: 15.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Martial Boggio-Pasqua
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Anthony Scemama
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Michel Caffarel
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Denis Jacquemin
- Laboratoire CEISAM (UMR 6230), CNRS, Université de Nantes, 44399 Cedex 3 Nantes, France
| |
Collapse
|
26
|
Influence of pseudopotentials on excitation energies from selected configuration interaction and diffusion Monte Carlo. RESULTS IN CHEMISTRY 2019. [DOI: 10.1016/j.rechem.2019.100002] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
|
27
|
Scemama A, Benali A, Jacquemin D, Caffarel M, Loos PF. Excitation energies from diffusion Monte Carlo using selected configuration interaction nodes. J Chem Phys 2018; 149:034108. [DOI: 10.1063/1.5041327] [Citation(s) in RCA: 44] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Anthony Scemama
- Laboratoire de Chimie et Physique Quantiques, Université de Toulouse, CNRS, UPS, Toulouse, France
| | - Anouar Benali
- Computational Science Division, Argonne National Laboratory, Argonne, Illinois 60439, USA
| | - Denis Jacquemin
- Laboratoire CEISAM—UMR CNRS 6230, Université de Nantes, 2 Rue de la Houssinière, BP 92208, 44322 Nantes Cedex 3, France
| | - Michel Caffarel
- Laboratoire de Chimie et Physique Quantiques, Université de Toulouse, CNRS, UPS, Toulouse, France
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques, Université de Toulouse, CNRS, UPS, Toulouse, France
| |
Collapse
|
28
|
Brémond E, Savarese M, Adamo C, Jacquemin D. Accuracy of TD-DFT Geometries: A Fresh Look. J Chem Theory Comput 2018; 14:3715-3727. [DOI: 10.1021/acs.jctc.8b00311] [Citation(s) in RCA: 47] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Affiliation(s)
- Eric Brémond
- Sorbonne Paris Cité, ITODYS, UMR CNRS 7086, Université Paris Diderot, 15 rue Jean-Antoine de Baı̈f, F-75013 Paris, France
| | - Marika Savarese
- D3 Compunet, Istituto Italiano di Tecnologia via Morego, 30 16163 Genova, Italy
| | - Carlo Adamo
- Chimie ParisTech, CNRS, Institut de Recherche de Chimie Paris, PSL Research University, 11 rue Pierre et Marie Curie, F-75005 Paris, France
| | - Denis Jacquemin
- Laboratoire CEISAM, UMR CNRS 6230, Université de Nantes, 2 rue de la Houssinière, BP 92208, 44322 Cedex 3 Nantes, France
| |
Collapse
|
29
|
Study of the low energy spectrum of titanium by using QMC methods. Chem Phys Lett 2018. [DOI: 10.1016/j.cplett.2018.01.004] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
30
|
Budzák Š, Scalmani G, Jacquemin D. Accurate Excited-State Geometries: A CASPT2 and Coupled-Cluster Reference Database for Small Molecules. J Chem Theory Comput 2017; 13:6237-6252. [PMID: 29140697 PMCID: PMC5729545 DOI: 10.1021/acs.jctc.7b00921] [Citation(s) in RCA: 44] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Abstract
![]()
We
present an investigation of the excited-state structural parameters
determined for a large set of small compounds with the dual goals
of defining reference values for further works and assessing the quality
of the geometries obtained with relatively cheap computational approaches.
In the first stage, we compare the excited-state geometries obtained
with ADC(2), CC2, CCSD, CCSDR(3), CC3, and CASPT2 and large atomic
basis sets. It is found that CASPT2 and CC3 results are generally
in very good agreement with one another (typical differences of ca.
3 × 10–3 Å) when all electrons are correlated
and when the aug-cc-pVTZ atomic basis set is employed with both methods.
In a second stage, a statistical analysis reveals that, on the one
hand, the excited-state (ES) bond lengths are much more sensitive
to the selected level of theory than their ground-state (GS) counterparts
and, on the other hand, that CCSDR(3) is probably the most cost-effective
method delivering accurate structures. Indeed, CCSD tends to provide
too compact multiple bond lengths on an almost systematic basis, whereas
both CC2 and ADC(2) tend to exaggerate these bond distances, with
more erratic error patterns, especially for the latter method. The
deviations are particularly marked for the polarized CO and CN bonds,
as well as for the puckering angle in formaldehyde homologues. In
the last part of this contribution, we provide a series of CCSDR(3)
GS and ES geometries of medium-sized molecules to be used as references
in further investigations.
Collapse
Affiliation(s)
- Šimon Budzák
- Department of Chemistry, Faculty of Natural Sciences, Matej Bel University , Tajovského 40, SK-97400 Banská Bystrica, Slovak Republic
| | - Giovanni Scalmani
- Gaussian Incorporated , 340 Quinnipiac Street, Building 40, Wallingford, Connecticut 06492 United States
| | - Denis Jacquemin
- Laboratoire CEISAM-UMR CNRS 6230, Université de Nantes , 2 Rue de la Houssiniére, BP 92208, 44322 Cedex 3 Nantes, France.,Institut Universitaire de France , 1 Rue Descartes, 75231 Cedex 5 Paris, France
| |
Collapse
|
31
|
McDaniel T, D’Azevedo EF, Li YW, Wong K, Kent PRC. Delayed Slater determinant update algorithms for high efficiency quantum Monte Carlo. J Chem Phys 2017; 147:174107. [DOI: 10.1063/1.4998616] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- T. McDaniel
- Department of Electrical Engineering and Computer Science, University of Tennessee, Knoxville, Tennessee 37996, USA
| | - E. F. D’Azevedo
- Computer Science and Mathematics Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - Y. W. Li
- National Center for Computational Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - K. Wong
- Joint Institute for Computational Sciences, University of Tennessee, Knoxville, Tennessee 37996, USA
| | - P. R. C. Kent
- Center for Nanophase Materials Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA and Computational Sciences and Engineering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| |
Collapse
|
32
|
Holmes AA, Umrigar CJ, Sharma S. Excited states using semistochastic heat-bath configuration interaction. J Chem Phys 2017; 147:164111. [DOI: 10.1063/1.4998614] [Citation(s) in RCA: 86] [Impact Index Per Article: 12.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Adam A. Holmes
- Department of Chemistry and Biochemistry, University of Colorado Boulder, Boulder, Colorado 80302, USA
- Laboratory of Atomic and Solid State Physics, Cornell University, Ithaca, New York 14853, USA
| | - C. J. Umrigar
- Laboratory of Atomic and Solid State Physics, Cornell University, Ithaca, New York 14853, USA
| | - Sandeep Sharma
- Department of Chemistry and Biochemistry, University of Colorado Boulder, Boulder, Colorado 80302, USA
| |
Collapse
|
33
|
Doblhoff-Dier K, Meyer J, Hoggan PE, Kroes GJ. Quantum Monte Carlo Calculations on a Benchmark Molecule-Metal Surface Reaction: H 2 + Cu(111). J Chem Theory Comput 2017; 13:3208-3219. [PMID: 28514594 PMCID: PMC5508338 DOI: 10.1021/acs.jctc.7b00344] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/02/2022]
Abstract
![]()
Accurate
modeling of heterogeneous catalysis requires the availability
of highly accurate potential energy surfaces. Within density functional
theory, these can—unfortunately—depend heavily on the
exchange-correlation functional. High-level ab initio calculations, on the other hand, are challenging due to the system
size and the metallic character of the metal slab. Here, we present
a quantum Monte Carlo (QMC) study for the benchmark system H2 + Cu(111), focusing on the dissociative chemisorption barrier height.
These computationally extremely challenging ab initio calculations agree to within 1.6 ± 1.0 kcal/mol with a chemically
accurate semiempirical value. Remaining errors, such as time-step
errors and locality errors, are analyzed in detail in order to assess
the reliability of the results. The benchmark studies presented here
are at the cutting edge of what is computationally feasible at the
present time. Illustrating not only the achievable accuracy but also
the challenges arising within QMC in such a calculation, our study
presents a clear picture of where we stand at the moment and which
approaches might allow for even more accurate results in the future.
Collapse
Affiliation(s)
- Katharina Doblhoff-Dier
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University , Post Office Box 9502, 2300 RA Leiden, The Netherlands
| | - Jörg Meyer
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University , Post Office Box 9502, 2300 RA Leiden, The Netherlands
| | - Philip E Hoggan
- Institute Pascal, UMR 6602 CNRS, University Blaise Pascal , 4 avenue Blaise Pascal, TSA 60026, CS 60026, 63178 Aubiere Cedex, France
| | - Geert-Jan Kroes
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University , Post Office Box 9502, 2300 RA Leiden, The Netherlands
| |
Collapse
|
34
|
Introduction to the Variational Monte Carlo Method in Quantum Chemistry and Physics. VARIATIONAL METHODS IN MOLECULAR MODELING 2017. [DOI: 10.1007/978-981-10-2502-0_10] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/10/2023]
|
35
|
Zhao L, Neuscamman E. Equation of Motion Theory for Excited States in Variational Monte Carlo and the Jastrow Antisymmetric Geminal Power in Hilbert Space. J Chem Theory Comput 2016; 12:3719-26. [DOI: 10.1021/acs.jctc.6b00480] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Luning Zhao
- Department
of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical
Science Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Eric Neuscamman
- Department
of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical
Science Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
36
|
Cleland DM, Per MC. Performance of quantum Monte Carlo for calculating molecular bond lengths. J Chem Phys 2016; 144:124108. [DOI: 10.1063/1.4944826] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Deidre M. Cleland
- CSIRO Virtual Nanoscience Laboratory, 343 Royal Parade, Parkville, Victoria 3052, Australia
| | - Manolo C. Per
- CSIRO Virtual Nanoscience Laboratory, 343 Royal Parade, Parkville, Victoria 3052, Australia
| |
Collapse
|
37
|
Kowalczyk T, Le K, Irle S. Self-Consistent Optimization of Excited States within Density-Functional Tight-Binding. J Chem Theory Comput 2015; 12:313-23. [DOI: 10.1021/acs.jctc.5b00734] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Tim Kowalczyk
- Institute
of Transformative Bio-Molecules (WPI-ITbM) and Department of Chemistry,
Graduate School of Science, Nagoya University, Nagoya 464-8602, Japan
- Department
of Chemistry, Advanced Materials Science and Engineering Center, and
Institute for Energy Studies, Western Washington University, Bellingham, Washington 98225, United States
| | - Khoa Le
- Department
of Chemistry, Advanced Materials Science and Engineering Center, and
Institute for Energy Studies, Western Washington University, Bellingham, Washington 98225, United States
| | - Stephan Irle
- Institute
of Transformative Bio-Molecules (WPI-ITbM) and Department of Chemistry,
Graduate School of Science, Nagoya University, Nagoya 464-8602, Japan
| |
Collapse
|
38
|
Dupuy N, Bouaouli S, Mauri F, Sorella S, Casula M. Vertical and adiabatic excitations in anthracene from quantum Monte Carlo: Constrained energy minimization for structural and electronic excited-state properties in the JAGP ansatz. J Chem Phys 2015; 142:214109. [DOI: 10.1063/1.4922048] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Affiliation(s)
- Nicolas Dupuy
- Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie, Université Pierre et Marie Curie, case 115, 4 place Jussieu, 75252 Paris Cedex 05, France
| | - Samira Bouaouli
- Laboratoire de Chimie Théorique, Université Pierre et Marie Curie, case 115, 4 place Jussieu, 75252 Paris Cedex 05, France
| | - Francesco Mauri
- CNRS and Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie, Université Pierre et Marie Curie, case 115, 4 place Jussieu, 75252 Paris Cedex 05, France
| | - Sandro Sorella
- International School for Advanced Studies (SISSA), Via Beirut 2-4, 34014 Trieste, Italy and INFM Democritos National Simulation Center, Trieste, Italy
| | - Michele Casula
- CNRS and Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie, Université Pierre et Marie Curie, case 115, 4 place Jussieu, 75252 Paris Cedex 05, France
| |
Collapse
|
39
|
Luo Y, Zen A, Sorella S. Ab initio molecular dynamics with noisy forces: Validating the quantum Monte Carlo approach with benchmark calculations of molecular vibrational properties. J Chem Phys 2014; 141:194112. [DOI: 10.1063/1.4901430] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
40
|
Guareschi R, Floris FM, Amovilli C, Filippi C. Solvent Effects on Excited-State Structures: A Quantum Monte Carlo and Density Functional Study. J Chem Theory Comput 2014; 10:5528-37. [DOI: 10.1021/ct500723s] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Affiliation(s)
- Riccardo Guareschi
- MESA+
Institute for Nanotechnology, University of Twente, P.O. Box 217, 7500 AE Enschede, The Netherlands
| | - Franca Maria Floris
- Dipartimento
di Chimica e Chimica Industriale, Università di Pisa, Via Giuseppe
Moruzzi 3, 56124 Pisa, Italy
| | - Claudio Amovilli
- Dipartimento
di Chimica e Chimica Industriale, Università di Pisa, Via Giuseppe
Moruzzi 3, 56124 Pisa, Italy
| | - Claudia Filippi
- MESA+
Institute for Nanotechnology, University of Twente, P.O. Box 217, 7500 AE Enschede, The Netherlands
| |
Collapse
|
41
|
Moroni S, Saccani S, Filippi C. Practical Schemes for Accurate Forces in Quantum Monte Carlo. J Chem Theory Comput 2014; 10:4823-9. [DOI: 10.1021/ct500780r] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Affiliation(s)
- S. Moroni
- CNR-IOM DEMOCRITOS,
Istituto Officina dei Materiali, and SISSA Scuola Internazionale Superiore
di Studi Avanzati, Via Bonomea 265, I-34136 Trieste, Italy
| | - S. Saccani
- SISSA Scuola Internazionale Superiore di Studi Avanzati, Via Bonomea 265, I-34136 Trieste, Italy
| | - C. Filippi
- MESA+
Institute for Nanotechnology, University of Twente, P.O. Box 217, 7500 AE Enschede, The Netherlands
| |
Collapse
|
42
|
Zen A, Trout BL, Guidoni L. Properties of reactive oxygen species by quantum Monte Carlo. J Chem Phys 2014; 141:014305. [DOI: 10.1063/1.4885144] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023] Open
Affiliation(s)
- Andrea Zen
- Dipartimento di Fisica, La Sapienza - Università di Roma, Piazzale Aldo Moro 2, 00185 Rome, Italy
| | - Bernhardt L. Trout
- Department of Chemical Engineering, Massachusetts Institute of Technology, 77 Massachusetts Ave, Cambridge, Massachusetts 02139, USA
| | - Leonardo Guidoni
- Dipartimento di Scienze Fisiche e Chimiche, Università degli studi de L'Aquila, Via Vetoio, 67100 Coppito, L'Aquila, Italy
| |
Collapse
|
43
|
Puzzarini C, Penocchio E, Biczysko M, Barone V. Molecular structure and spectroscopic signatures of acrolein: theory meets experiment. J Phys Chem A 2014; 118:6648-56. [PMID: 24842714 DOI: 10.1021/jp503672g] [Citation(s) in RCA: 32] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
A comprehensive study of the molecular structure and IR spectrum of cis and trans acrolein has been performed by an integrated computational approach coupling methods rooted in the coupled-cluster ansatz and the density functional theory. From the one side, DFT anharmonic force fields allow us to determine very reliable semiexperimental structures for both isomers, which are in remarkable agreement with the geometries issuing from CCSD(T) computations accounting for the extrapolation to the complete basis set and core correlation. The same kind of coupled-cluster computations provide dipole moment, relative energies, and interconversion barrier in remarkable agreement with experiments. Finally, harmonic CCSD(T) results coupled to DFT evaluation of mechanical and electrical anharmonicity allow us, in the framework of second-order perturbative vibrational theory, to confirm most of the experimental assignments of IR spectra, and to suggest some additional interpretations for congested regions including fundamental bands together with overtones and combination bands.
Collapse
Affiliation(s)
- Cristina Puzzarini
- Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna , Via Selmi 2, I-40126 Bologna, Italy
| | | | | | | |
Collapse
|