Search NASA⌕ Search

Engineering topics

Scalmani, Giovanni

Publications and source records attributed to Scalmani, Giovanni.

8 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Modeling Multi-Step Organic Reactions: Can Density Functional Theory Deliver Misleading Chemistry?

Many organic reactions are characterized by a complex mechanism with a variety of transition states and intermediates of different chemical natures. Their correct and accurate theoretical characterization critically depends on the accuracy of the computational method used. In this work, we study a complex ambimodal cycloaddition with five transition states, two intermediates, and three products, and we ask whether density functional theory (DFT) can provide a correct description of this type of complex and multifaceted reaction. Our work fills a gap in that most systematic benchmarks of DFT for chemical reactions have considered much simpler reactions. Our results show that many density functionals not only lead to seriously large errors but also differ from one another in predicting whether the reaction is ambimodal. Only a few of the available functionals provide a balanced description of the complex and multifaceted reactions. The parameters varied in the tested functionals are the ingredients, the treatment of medium-range and nonlocal correlation energy, and the inclusion of Hartree–Fock exchange. Furthermore, these results show a clear need for more benchmarks on the mechanisms of large molecules in complex reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Barrier Heights for Diels–Alder Transition States Leading to Pentacyclic Adducts: A Benchmark Study of Crowded, Strained Transition States of Large Molecules

Theoretical characterization of reactions of complex molecules depends on providing consistent accuracy for the relative energies of intermediates and transition states. Here we employ the DLPNO-CCSD(T) method with core–valence correlation, large basis sets, and extrapolation to the CBS limit to provide benchmark values for Diels–Alder transition states leading to competitive strained pentacyclic adducts. We then used those benchmarks to test a diverse set of wave function and density functional methods for the absolute and relative barrier heights of these transition states. Our results show that only a few of the tested density functionals can predict the absolute barrier heights satisfactorily, although relative barrier heights are more accurate. The most accurate functionals tested are ωB97M-V, M11plus, ωB97X-V, PBE-D3(0), M11, and MN15 with MUDs from best estimates less than 3.0 kcal. Furthermore, these findings can guide selection of density functionals for future studies of crowded, strained transition states of large molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

M11plus: A Range-Separated Hybrid Meta Functional with Both Local and Rung-3.5 Correlation Terms and High Across-the-Board Accuracy for Chemical Applications

The way to improve Kohn–Sham density functional theory is to improve the exchange–correlation functionals, and functionals have been successively improved by adding new ingredients, especially local spin density gradients, nonlocal Hartree–Fock exchange, and local meta terms based on kinetic energy density. Here, we present a new kind of functional obtained by adding rung-3.5 terms to a functional including local gradients, local meta terms, and range-separated Hartree–Fock exchange. A rung-3.5 term has short-range nonlocality designed to account for nondynamic correlation; we add two kinds of rung-3.5 terms, one kind modeled on position-dependent Hartree–Fock exchange and another modeled on the spin density at a point interacting with the opposite-spin exchange hole at the same point. Here, optimization of the functional yields broad accuracy for both ground states and excited states with especially significant improvement for systems with strong correlation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

Effective Inclusion of Mechanical and Electrical Anharmonicity in Excited Electronic States: VPT2-TDDFT Route

Here, we present a reliable and cost-effective procedure for the inclusion of anharmonic effects in excited-state energies and spectroscopic intensities by means of second-order vibrational perturbation theory. This development is made possible thanks to a recent efficient implementation of excited-state analytic Hessians and properties within the time-dependent density functional theory framework. As illustrated in this work, by taking advantage of such algorithmic developments, it is possible to perform calculations of excited-state infrared spectra of medium-large isolated molecular systems, with anharmonicity effects included in both the energy and property surfaces. We also explore the use of this procedure for the inclusion of anharmonic effects in the simulation of vibronic bandshapes of electronic spectra, and compare the results with previous, more approximate models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

Quantum, classical, and hybrid QM/MM calculations in solution: General implementation of the ddCOSMO linear scaling strategy

We present the general theory and implementation of the Conductor-like Screening Model according to the recently developed ddCOSMO paradigm. The various quantities needed to apply ddCOSMO at different levels of theory, including quantum mechanical descriptions, are discussed in detail, with a particular focus on how to compute the integrals needed to evaluate the ddCOSMO solvation energy and its derivatives. The overall computational cost of a ddCOSMO computation is then analyzed and decomposed in the various steps: the different relative weights of such contributions are then discussed for both ddCOSMO and the fastest available alternative discretization to the COSMO equations. Finally, the scaling of the cost of the various steps with respect to the size of the solute is analyzed and discussed, showing how ddCOSMO opens significantly new possibilities when cheap or hybrid molecular mechanics/quantum mechanics methods are used to describe the solute.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

How far do electrons delocalize?

Electron delocalization is central to chemical bonding, but it is also a fundamentally nonclassical and nonintuitive quantum mechanical phenomenon. Tools to quantify and visualize electron delocalization help to understand, teach, and predict chemical reactivity. We develop a new approach to quantify and visualize electron delocalization in real space. Our electron delocalization range function EDR(r{sup -vector};u) quantifies the degree to which electrons at point r{sup -vector} in a calculated wavefunction delocalize over length scale u. Its predictions are physically reasonable. For example, EDR(r{sup -vector};u=0.25 bohr) is close to one at points r{sup -vector} in the cores of first-row atoms, consistent with the localization of core electrons to ∼0.25 bohr. EDR(r{sup -vector};u=1 bohr) is close to one at points r{sup -vector} in typical covalent bonds, consistent with electrons delocalizing over the length of the bond. Our approach provides a rich representation of atomic shell structure; covalent and ionic bonding; the delocalization of excited states, defects, and solvated electrons; metallic and insulating systems; and bond stretching and strong correlation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

Practical auxiliary basis implementation of Rung 3.5 functionals

Approximate exchange-correlation functionals for Kohn-Sham density functional theory often benefit from incorporating exact exchange. Exact exchange is constructed from the noninteracting reference system's nonlocal one-particle density matrix γ(r{sup -vector},r{sup -vector}′). Rung 3.5 functionals attempt to balance the strengths and limitations of exact exchange using a new ingredient, a projection of γ(r{sup -vector},r{sup -vector} ′) onto a semilocal model density matrix γ{sub SL}(ρ(r{sup -vector}),∇ρ(r{sup -vector}),r{sup -vector}−r{sup -vector} ′). γ{sub SL} depends on the electron density ρ(r{sup -vector}) at reference point r{sup -vector}, and is closely related to semilocal model exchange holes. We present a practical implementation of Rung 3.5 functionals, expanding the r{sup -vector}−r{sup -vector} ′ dependence of γ{sub SL} in an auxiliary basis set. Energies and energy derivatives are obtained from 3D numerical integration as in standard semilocal functionals. We also present numerical tests of a range of properties, including molecular thermochemistry and kinetics, geometries and vibrational frequencies, and bandgaps and excitation energies. Rung 3.5 functionals typically provide accuracy intermediate between semilocal and hybrid approximations. Nonlocal potential contributions from γ{sub SL} yield interesting successes and failures for band structures and excitation energies. The results enable and motivate continued exploration of Rung 3.5 functional forms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗