1
|
Wennrich JP, Holzenkamp C, Kolařík M, Maier W, Mándi A, Kurtán T, Ashrafi S, Ebada SS, Stadler M. Dactylfungins and Tetralones: Bioactive Metabolites from a Nematode-Associated Laburnicola nematophila. JOURNAL OF NATURAL PRODUCTS 2024; 87:1860-1871. [PMID: 39012621 PMCID: PMC11287750 DOI: 10.1021/acs.jnatprod.4c00623] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/27/2024] [Revised: 06/29/2024] [Accepted: 06/30/2024] [Indexed: 07/17/2024]
Abstract
A chemical investigation of Laburnicola nematophila, isolated from cysts of the plant parasitic nematode Heterodera filipjevi, affored three dactylfungin derivatives (1-3) and three tetralone congeners (4-6). Dactylfungin C (1), laburnicolin (4), and laburnicolenone (5) are previously undescribed natural products. Chemical structures of the isolated compounds were determined based on 1D and 2D NMR spectroscopic analyses together with HR-ESI-MS spectrometry and comparison with data reported in the literature. The relative configurations of compounds 1, 2, and 4-6 were determined based on their ROESY data and analysis of their coupling constants (J values). The absolute configurations of 4-6 were determined through the comparison of their measured and calculated TDDFT-ECD spectra. Compounds 1-3 were active against azole-resistant Aspergillus fumigatus.
Collapse
Affiliation(s)
- Jan-Peer Wennrich
- Department
of Microbial Drugs, Helmholtz Centre for
Infection Research (HZI) and German Centre for Infection Research, Inhoffenstrasse 7, 38124 Braunschweig, Germany
- Institute
of Microbiology, Technische Universität
Braunschweig, Spielmannstraße
7, 38106 Braunschweig, Germany
| | - Caren Holzenkamp
- Department
of Microbial Drugs, Helmholtz Centre for
Infection Research (HZI) and German Centre for Infection Research, Inhoffenstrasse 7, 38124 Braunschweig, Germany
- Institute
of Microbiology, Technische Universität
Braunschweig, Spielmannstraße
7, 38106 Braunschweig, Germany
| | - Miroslav Kolařík
- Institute
of Microbiology, Czech Academy of Science, Vídeňská 1083, 14220 Prague, Czech Republic
| | - Wolfgang Maier
- Institute
for Epidemiology and Pathogen Diagonstics, Julius Kühn Institut (JKI) - Federal Research Centre for Cultivated
Plants, Messeweg 11-12, 38104 Braunschweig, Germany
| | - Attila Mándi
- Department
of Organic Chemistry, University of Debrecen, P.O. Box 400, 4002 Debrecen, Hungary
| | - Tibor Kurtán
- Department
of Organic Chemistry, University of Debrecen, P.O. Box 400, 4002 Debrecen, Hungary
| | - Samad Ashrafi
- Institute
for Epidemiology and Pathogen Diagonstics, Julius Kühn Institut (JKI) - Federal Research Centre for Cultivated
Plants, Messeweg 11-12, 38104 Braunschweig, Germany
- Institute
for Crop and Soil Science, Julius Kühn
Institute (JKI) − Federal Research Centre for Cultivated Plants, Bundesallee 58, 38116 Braunschweig, Germany
| | - Sherif S. Ebada
- Department
of Microbial Drugs, Helmholtz Centre for
Infection Research (HZI) and German Centre for Infection Research, Inhoffenstrasse 7, 38124 Braunschweig, Germany
- Department
of Pharmacognosy, Faculty of Pharmacy, Ain
Shams University, Cairo 11566, Egypt
| | - Marc Stadler
- Department
of Microbial Drugs, Helmholtz Centre for
Infection Research (HZI) and German Centre for Infection Research, Inhoffenstrasse 7, 38124 Braunschweig, Germany
| |
Collapse
|
2
|
Miller ER, Hoehn SJ, Kumar A, Jiang D, Parker SM. Ultrafast photochemistry and electron diffraction for cyclobutanone in the S2 state: Surface hopping with time-dependent density functional theory. J Chem Phys 2024; 161:034105. [PMID: 39007373 DOI: 10.1063/5.0203679] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2024] [Accepted: 06/26/2024] [Indexed: 07/16/2024] Open
Abstract
We simulate the photodynamics of gas-phase cyclobutanone excited to the S2 state using fewest switches surface hopping (FSSH) dynamics powered by time-dependent density functional theory (TDDFT). We predict a total photoproduct yield of 8%, with a C3:C2 product ratio of 0 trajectories to 8 trajectories. One primary S2 → S1 conical intersection is identified involving the compression of an α-carbon-carbon-hydrogen bond angle. Excited state lifetimes computed with respect to electronic state populations were found to be 3.96 ps (S2 → S1) and 498 fs (S1 → S0). We also generate time-resolved difference pair distribution functions (ΔPDFs) from our TDDFT-FSSH dynamics results in order to generate direct comparisons with ultrafast electron diffraction experiment observables. Global and target analysis of time-resolved ΔPDFs produced a distinct set of lifetimes: (i) a 0.548 ps decay and (ii) a 1.69 ps decay, both resembling the S2 minimum, as well as (iii) a long decay that resembles the S1 minimum geometry and the fully separated C2 products. Finally, we contextualize our results by considering the impact of the most likely sources of significant errors.
Collapse
Affiliation(s)
- Ericka Roy Miller
- Department of Chemistry, Case Western Reserve University, 10900 Euclid Ave., Cleveland, Ohio 44106, USA
| | - Sean J Hoehn
- Department of Chemistry, Case Western Reserve University, 10900 Euclid Ave., Cleveland, Ohio 44106, USA
| | - Abhijith Kumar
- Department of Chemistry, Case Western Reserve University, 10900 Euclid Ave., Cleveland, Ohio 44106, USA
| | - Dehua Jiang
- Department of Chemistry, Case Western Reserve University, 10900 Euclid Ave., Cleveland, Ohio 44106, USA
| | - Shane M Parker
- Department of Chemistry, Case Western Reserve University, 10900 Euclid Ave., Cleveland, Ohio 44106, USA
| |
Collapse
|
3
|
Otlyotov AA, Moshchenkov AD, Minenkov Y. Ni, Cu, Zn, Pd, Ag and Cd Tetraphenylporphyrin Ab Initio Thermochemistry: Enthalpy of Formation of ZnTPP Revisited. Inorg Chem 2024; 63:10230-10239. [PMID: 38780084 DOI: 10.1021/acs.inorgchem.4c00662] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/25/2024]
Abstract
Groups 10-12 metalloporphyrins have been recognized for their numerous properties essential for the development of new sensing materials. In this work, accurate gas-phase enthalpies of formation, ΔfHm0(g,298.15), are predicted for the series of Ni, Cu, Zn, Pd, Ag, and Cd tetraphenylporphyrins (MTPPs) on the basis of the reaction-based Feller-Peterson-Dixon approach and high-level ab initio DLPNO-CCSD(T) calculations. Our recently developed automatic generator of the balanced chemical reactions was employed to reduce the bias of the theoretical ΔfHm0(g,298.15) toward a particular reaction. Theoretical ΔfHm0(g,298.15) for ZnTPP (227.0 ± 3.4 kcal mol-1) does not support the previously reported experimental value of 132 ± 2 kcal mol-1. The origin of the discrepancy probably lies in the experimental solid-state ΔfHm0(ZnTPP, cr,298.15) as it stems from our theoretical evaluations of the ΔfHm0(cr,298.15) values for the entire set of transition metal TPP complexes. The large discrepancy between experiment and theory also holds when different DFT functionals (ωB97M-V, PBE0-D4, and B3LYP-D4) paired with quadruple-ζ quality basis sets are used for the theoretical calculations. Experimental revisiting of the solid-state enthalpy of formation of ZnTPP and analogue measurements for other transition metal TPPs are needed to resolve the observed discrepancy. Based on the predicted enthalpies of formation of MTPPs, the relative energies of the metal-ligand bonding are evaluated and the trends are compared to those for the complexes of the unsubstituted porphyrin with the same set of metals derived in [Can. J. Chem., 2009, 87, 1063]. According to both studies, Pd complexes exhibit the strongest bonding, while the Cd species are the least stable metallocomplexes within the considered series.
Collapse
Affiliation(s)
- Arseniy A Otlyotov
- N. N. Semenov Federal Research Center for Chemical Physics RAS, Kosygina Street 4, 119991 Moscow, Russian Federation
| | - Andrey D Moshchenkov
- N. N. Semenov Federal Research Center for Chemical Physics RAS, Kosygina Street 4, 119991 Moscow, Russian Federation
| | - Yury Minenkov
- N. N. Semenov Federal Research Center for Chemical Physics RAS, Kosygina Street 4, 119991 Moscow, Russian Federation
| |
Collapse
|
4
|
Roy M, Dahmani R, Vallet V, Masella M. Early Steps in the O 2 Scavenger Process in the Aqueous Phase: Hydrazine vs DEHA. J Phys Chem A 2023; 127:10104-10117. [PMID: 37988629 DOI: 10.1021/acs.jpca.3c05383] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2023]
Abstract
We investigate the first direct proton abstraction reactions from reducing agents (RAHs) hydrazine and diethyl hydroxylamine (DEHA), toward dioxygen (O2) in the aqueous phase, spanning ambient to high-temperature conditions. Quantum chemistry methods and molecular dynamics simulations are employed in this study. Quantum chemistry methods are used to analyze the quasi-equilibrium between a reactive conformation and a transition state in the [RAH,O2] cluster. On the other hand, molecular dynamics simulations estimate the probability of observing a reactive conformation of the [RAH,O2] cluster in the solution. In this study, we assume that the energy barrier of the quasi-equilibrium is sufficiently high for the RAH/O2 association process to be at equilibrium. Our findings indicate that the first proton abstraction process from a reactive conformation cluster by DEHA is energetically favored compared to hydrazine. Conversely, the association process of hydrazine and O2 in solution is more favorable than that of DEHA. Consequently, the rate constant for the first proton abstraction process is similar for both hydrazine and DEHA, particularly at high temperatures, with activation energies of approximately 21.5 ± 1.5 kcal mol-1 for both compounds. These results align with recent experiments investigating the complete O2 scavenger process in liquid water with hydrazine and DEHA. Therefore, our findings support the assumption that first proton abstraction reactions are the rate-determining steps in O2 scavenger processes in the aqueous phase.
Collapse
Affiliation(s)
- Marion Roy
- Service de Physico-Chimie, Université Paris-Saclay, CEA, Gif-sur-Yvette 91190, France
| | - Rahma Dahmani
- Laboratoire de Biologie Structurale et Radiobiologie, Service de Bioénergétique, Biologie Structurale et Mécanismes, Institut de Biologie et de Technologies de Saclay, CEA Saclay, Gif sur Yvette Cedex F-91191, France
| | - Valérie Vallet
- Univ. Lille, CNRS, UMR 8523─PhLAM─Physique des Lasers Atomes et Molécules, Lille F-59000, France
| | - Michel Masella
- Laboratoire de Biologie Structurale et Radiobiologie, Service de Bioénergétique, Biologie Structurale et Mécanismes, Institut de Biologie et de Technologies de Saclay, CEA Saclay, Gif sur Yvette Cedex F-91191, France
| |
Collapse
|
5
|
Wang M, He X, Taylor M, Lorpaiboon W, Mun H, Ho J. Molecular Geometries and Vibrational Contributions to Reaction Thermochemistry Are Surprisingly Insensitive to the Choice of Basis Sets. J Chem Theory Comput 2023. [PMID: 37463146 DOI: 10.1021/acs.jctc.3c00388] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/20/2023]
Abstract
Calculation of molecular geometries and harmonic vibrational frequencies are pre-requisites for thermochemistry calculations. Contrary to conventional wisdom, this paper demonstrates that quantum chemical predictions of the thermochemistry of many gas and solution phase chemical reactions appear to be very insensitive to the choice of basis sets. For a large test set of 80 diverse organic and transition-metal-containing reactions, variations in reaction free energy based on geometries and frequencies calculated using a variety of double and triple-zeta basis sets from the Pople, Jensen, Ahlrichs, and Dunning families are typically less than 4 kJ mol-1, especially when the quasiharmonic oscillator correction is applied to mitigate the effects of low-frequency modes. Our analysis indicates that for many organic molecules and their transition states, high-level revDSD-PBEP86-D4 and DLPNO-CCSD(T)/(aug-)cc-pVTZ single-point energies usually vary by less than 2 kJ mol-1 on density functional theory geometries optimized using basis sets ranging from 6-31+G(d) to aug-pcseg-2 and aug-cc-pVTZ. In cases where these single-point energies vary significantly, indicating sensitivity of molecular geometries to the choice of basis set, there is often substantial cancellation of errors when the reaction energy or barrier is calculated. The study concludes that the choice of basis set for molecular geometry and frequencies, particularly those considered in this study, is not critical for the accuracy of thermochemistry calculations in the gas or solution phase.
Collapse
Affiliation(s)
- Minzhi Wang
- School of Chemistry, The University of New South Wales, Sydney, New South Wales 2052, Australia
| | - Xinlan He
- School of Chemistry, The University of New South Wales, Sydney, New South Wales 2052, Australia
| | - Mackenzie Taylor
- School of Chemistry, The University of New South Wales, Sydney, New South Wales 2052, Australia
| | - Wanutcha Lorpaiboon
- School of Chemistry, The University of New South Wales, Sydney, New South Wales 2052, Australia
| | - Haedam Mun
- School of Chemistry, The University of New South Wales, Sydney, New South Wales 2052, Australia
| | - Junming Ho
- School of Chemistry, The University of New South Wales, Sydney, New South Wales 2052, Australia
| |
Collapse
|
6
|
Minenkov Y. Solv: An Alternative Continuum Model Implementation Based on Fixed Atomic Charges, Scaled Particle Theory, and the Atom-Atom Potential Method. J Chem Theory Comput 2023. [PMID: 37390470 DOI: 10.1021/acs.jctc.3c00410] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/02/2023]
Abstract
An alternative continuum model implementation is reported. The electrostatic contribution to the solvation Gibbs free energy utilizes the noniterative conductor-like screening model of Vyboishchikov and Voityuk (DOI: 10.1002/jcc.26531) based on the fixed partial atomic charges. The nonelectrostatic solute-solvent dispersion-repulsion energy is calculated through the Caillet-Claverie atom-atom potential method employing the grid-based approach. The nonelectrostatic cavitation energy is computed within the scaled particle theory (SPT) formalism with the solute hard-sphere radius obtained via the Pierotti-Claverie (PC) scheme, from the solute molecular surface (SPT-S) or volume (SPT-V). The solvent hard-sphere radius is derived through the fitting to the experimental total solvation free energies of 2530 neutral species in 92 solvents. Application of the model to reproduce both the absolute and relative (reaction net) solvation free energies indicates that the SPT-V approach based on the CM5 charges is the best performer. The method is suggested for the solvation free energy calculation in the nonaqueous solvents.
Collapse
Affiliation(s)
- Yury Minenkov
- N. N. Semenov Federal Research Center for Chemical Physics RAS, Kosygina Street 4, 119991 Moscow, Russian Federation
- Joint Institute for High Temperatures, Russian Academy of Sciences, 13-2 Izhorskaya Street, Moscow 125412, Russian Federation
| |
Collapse
|
7
|
Charvati E, Sun H. Potential Energy Surfaces Sampled in Cremer-Pople Coordinates and Represented by Common Force Field Functionals for Small Cyclic Molecules. J Phys Chem A 2023; 127:2646-2663. [PMID: 36893434 DOI: 10.1021/acs.jpca.3c00095] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/11/2023]
Abstract
The complex conformations of the cyclic moieties impact the physical and chemical properties of molecules. In this work, we chose 22 molecules of four-, five-, and six-membered rings and performed a thorough conformational sampling using Cremer-Pople coordinates. With consideration of symmetries, we obtained a total of 1504 conformational structures for four-membered, 5576 for five-membered, and 13509 for six-membered rings. All well-known and many less well-known conformers for each molecule were identified. We represented the potential energy surfaces (PESs) by fitting the data to common analytical force field (FF) functional forms. We found that the general features of PESs can be described by the essential FF functional forms; however, the accuracy of representation can be improved remarkably by including the torsion-bond and torsion-angle coupling terms. The best fit yields R-squared (R2) values close to 1.0 and mean absolute errors in energy less than 0.3 kcal/mol.
Collapse
Affiliation(s)
- Evangelia Charvati
- School of Chemistry and Chemical Engineering, Materials Genome Initiative Center, and Key Laboratory of Scientific and Engineering Computing of Ministry of Education, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Huai Sun
- School of Chemistry and Chemical Engineering, Materials Genome Initiative Center, and Key Laboratory of Scientific and Engineering Computing of Ministry of Education, Shanghai Jiao Tong University, Shanghai 200240, China
| |
Collapse
|
8
|
Comparative analysis of the reactivity of anthocyanidins, leucoanthocyanidins, and flavonols using a quantum chemistry approach. J Mol Model 2023; 29:93. [PMID: 36905478 PMCID: PMC10006565 DOI: 10.1007/s00894-023-05468-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2022] [Accepted: 02/05/2023] [Indexed: 03/12/2023]
Abstract
Anthocyanidins, leucoanthocyanidins, and flavonols are natural compounds mainly known due to their reported biological activities, such as antiviral, antifungal, anti-inflammatory activities, and antioxidant activity. In the present study, we performed a comparative structural, conformational, electronic, and nuclear magnetic resonance analysis of the reactivity of the chemical structure of primary anthocyanidins, leucoanthocyanidins, and flavonoids. We focused our analysis on the following molecular questions: (i) differences in cyanidin catechols ( +)-catechin, leucocyanidin, and quercetin; (ii) the loss of hydroxyl presents in the R1 radical of leucoanthocyanidin in the functional groups linked to C4 (ring C); and (iii) the electron affinity of the 3-hydroxyl group (R7) in the flavonoids delphinidin, pelargonidin, cyanidin, quercetin, and kaempferol. We show unprecedented results for bond critical point (BCP) of leucopelargonidin and leucodelphirinidin. The BCP formed between hydroxyl hydrogen (R2) and ketone oxygen (R1) of kaempferol has the same degrees of covalence of quercetin. Kaempferol and quercetin exhibited localized electron densities between hydroxyl hydrogen (R2) and ketone oxygen (R1). Global molecular descriptors showed quercetin and leucocyanidin are the most reactive flavonoids in electrophilic reactions. Complementary, anthocyanidins are the most reactive in nucleophilic reactions, while the smallest gap occurs in delphinidin. Local descriptors indicate that anthocyanidins and flavonols are more prone to electrophilic attacks, while in leucoanthocyanidins, the most susceptible to attack are localized in the ring A. The ring C of anthocyanidins is more aromatic than the same found in flavonols and leucoanthocyanidins. METHODS: For the analysis of the molecular properties, we used the DFT to evaluate the formation of the covalent bonds and intermolecular forces. CAM-B3LYP functional with the def2TZV basis set was used for the geometry optimization. A broad analysis of quantum properties was performed using the assessment of the molecular electrostatic potential surface, electron localization function, Fukui functions, descriptors constructed from frontier orbitals, and nucleus independent chemical shift.
Collapse
|
9
|
Minenkov Y, Cavallo L, Peterson KA. Influence of the complete basis set approximation, tight weighted-core, and diffuse functions on the DLPNO-CCSD(T1) atomization energies of neutral H,C,O-compounds. J Comput Chem 2023; 44:687-696. [PMID: 36399072 DOI: 10.1002/jcc.27033] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2022] [Revised: 09/23/2022] [Accepted: 09/28/2022] [Indexed: 11/19/2022]
Abstract
The impact of complete basis set extrapolation schemes (CBS), diffuse functions, and tight weighted-core functions on enthalpies of formation predicted via the DLPNO-CCSD(T1) reduced Feller-Peterson-Dixon approach has been examined for neutral H,C,O-compounds. All tested three-point (TZ/QZ/5Z) extrapolation schemes result in mean unsigned deviation (MUD) below 2 kJ mol-1 relative to the experiment. The two-point QZ/5Z and TZ/QZ CBS 1 / l max 3 extrapolation schemes are inferior to their inverse power counterpart ( 1 / l max + 1 / 2 4 ) by 1.3 and 4.3 kJ mol-1 . The CBS extrapolated frozen core atomization energies are insensitive (within 1 kJ mol-1 ) to augmentation of the basis set with tight weighted core functions. The core-valence correlation effects converge already at triple-ζ, although double-ζ/triple-ζ CBS extrapolation performs better and is recommended. The effect of diffuse function augmentation converges slowly, and cannot be reproduced with double- ζ or triple- ζ calculations as these are plagued with basis set superposition and incompleteness errors.
Collapse
Affiliation(s)
- Yury Minenkov
- N.N. Semenov Federal Research Center for Chemical Physics RAS, Moscow, Russian Federation.,Joint Institute for High Temperatures, Russian Academy of Sciences, Moscow, Russian Federation
| | - Luigi Cavallo
- KAUST Catalysis Center (KCC), King Abdullah University of Science and Technology, Thuwal, Saudi Arabia
| | - Kirk A Peterson
- Department of Chemistry, Washington State University, Pullman, Washington, USA
| |
Collapse
|
10
|
Otlyotov AA, Itkis D, Yashina LV, Cavallo L, Minenkov Y. Physical and numerical aspects of sodium ion solvation free energies via the cluster-continuum model. Phys Chem Chem Phys 2022; 24:29927-29939. [PMID: 36468644 DOI: 10.1039/d2cp03583a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
Abstract
Sodium cation solvation Gibbs free energies (ΔGsolv(Na+)) have been obtained in water, dimethylformamide, dimethyl sulfoxide, ethanol, acetone, acetonitrile, and methanol through the "monomer cycle" cluster-continuum approach where a solvent reference state is described by infinitely separated molecules. The following steps are vital for obtaining reliable ΔGsolv(Na+) values: (a) a meticulous conformational search involving dispersion corrected density functional theory (DFT-D) and the continuum solvation model (CSM); (b) gas-phase DFT-D geometry optimization followed by single-point (SP) domain-based local pair natural orbital coupled clusters including single, double, and partly triple excitation (DLPNO-CCSD(T)) calculations in conjunction with the complete basis set extrapolation; (c) advanced statistical thermodynamic treatment of the low harmonic frequencies (<100 cm-1) to obtain the robust gas-phase Gibbs free energy correction; (d) gas-phase and dielectric continuum SP with non-electrostatic contributions included in the CSM; (e) an evaluation of the relative thermodynamic stability of the Na+(S)n clusters to identify the number of explicit solvent molecules n to be considered. Our refined computational protocol is promising with a Pearson correlation coefficient between the predicted and experimental data, ρ, of 0.82, and the mean signed and mean unsigned errors of 0.3 and 1.4 kcal mol-1, respectively.
Collapse
Affiliation(s)
- Arseniy A Otlyotov
- N.N. Semenov Federal Research Center for Chemical Physics RAS, Kosygina Street 4, 119991 Moscow, Russia.
| | - Daniil Itkis
- N.N. Semenov Federal Research Center for Chemical Physics RAS, Kosygina Street 4, 119991 Moscow, Russia. .,Lomonosov Moscow State University, Leninskie Gory 1, Bld. 3, 119991 Moscow, Russia
| | - Lada V Yashina
- N.N. Semenov Federal Research Center for Chemical Physics RAS, Kosygina Street 4, 119991 Moscow, Russia. .,Lomonosov Moscow State University, Leninskie Gory 1, Bld. 3, 119991 Moscow, Russia
| | - Luigi Cavallo
- KAUST Catalysis Center (KCC), King Abdullah University of Science and Technology, Thuwal-23955-6900, Saudi Arabia.
| | - Yury Minenkov
- N.N. Semenov Federal Research Center for Chemical Physics RAS, Kosygina Street 4, 119991 Moscow, Russia. .,Joint Institute for High Temperatures, Russian Academy of Sciences, 13-2 Izhorskaya Street, Moscow 125412, Russia
| |
Collapse
|
11
|
Otlyotov AA, Cavallo L, Minenkov Y. Cluster-Continuum Model as a Sanity Check of Sodium Ions’ Gibbs Free Energies of Transfer. Inorg Chem 2022; 61:18365-18379. [DOI: 10.1021/acs.inorgchem.2c02065] [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)
- Arseniy A. Otlyotov
- N. N. Semenov Federal Research Center for Chemical Physics RAS, Kosygina Street 4, Moscow 119991, Russian Federation
| | - Luigi Cavallo
- KAUST Catalysis Center (KCC), King Abdullah University of Science and Technology, Thuwal 23955-6900, Saudi Arabia
| | - Yury Minenkov
- N. N. Semenov Federal Research Center for Chemical Physics RAS, Kosygina Street 4, Moscow 119991, Russian Federation
- Joint Institute for High Temperatures, Russian Academy of Sciences, 13-2 Izhorskaya Street, Moscow 125412, Russian Federation
| |
Collapse
|
12
|
Abstract
![]()
In this paper, we extend the rank-reduced coupled-cluster
formalism
to the calculation of non-iterative energy corrections due to quadruple
excitations. There are two major components of the proposed formalism.
The first is an approximate compression of the quadruple excitation
amplitudes using the Tucker format. The second is a modified functional
used for the evaluation of the corrections which gives exactly the
same results for the exact amplitudes, but is less susceptible to
errors resulting from the aforementioned compression. We show, both
theoretically and numerically, that the computational cost of the
proposed method scales as the seventh power of the system size. Using
reference results for a set of small molecules, the method is calibrated
to deliver relative accuracy of a few percent in energy corrections.
To illustrate the potential of the theory, we calculate the isomerization
energy of ortho/meta benzyne (C6H4) and the barrier height for the Cope rearrangement
in bullvalene (C10H10). The method retains a
near-black-box nature of the conventional coupled-cluster formalism
and depends on only one additional parameter that controls the accuracy.
Collapse
Affiliation(s)
- Michał Lesiuk
- Faculty of Chemistry, University of Warsaw, Pasteura 1, Warsaw 02-093, Poland
| |
Collapse
|
13
|
Otlyotov AA, Minenkov Y, Zaitsau DH, Zherikova KV, Verevkin SP. "In Vitro" and "In Vivo" Diagnostic Check for the Thermochemistry of Metal-Organic Compounds. Inorg Chem 2022; 61:10743-10755. [PMID: 35797430 DOI: 10.1021/acs.inorgchem.2c00959] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Volatile metal β-diketonates are of interest from both practical and theoretical perspectives (manufacturing of film materials, catalysis, and the nature of metal-ligand bonding). Knowledge of their reliable thermochemical properties is essential for effective applications. However, there is an unacceptable scattering of the available data on the enthalpies of formation. In this work, we proposed "in vitro" and "in vivo" diagnostic tools to verify the available enthalpies of formation in both the crystalline and gaseous states for metal tris-β-diketonates. The "in vitro" procedure involved high-level quantum-chemical calculations and was applied to define a consistent data set on the enthalpies of formation for iron(III) β-diketonates. This data set has provided the basis for "in vivo" structure-property-based diagnostics to evaluate the robustness of the thermochemical data for β-diketonate tris-complexes with metals other than iron.
Collapse
Affiliation(s)
- Arseniy A Otlyotov
- N. N. Semenov Federal Research Center for Chemical Physics RAS, Moscow 119991, Russian Federation
| | - Yury Minenkov
- N. N. Semenov Federal Research Center for Chemical Physics RAS, Moscow 119991, Russian Federation.,Joint Institute for High Temperatures, Russian Academy of Sciences, Moscow 125412, Russian Federation
| | - Dzmitry H Zaitsau
- Department of Physical Chemistry and Faculty of Interdisciplinary Research, Competence Centre CALOR, University of Rostock, 18059 Rostock, Germany
| | - Kseniya V Zherikova
- Nikolaev Institute of Inorganic Chemistry of Siberian Branch of Russian Academy of Sciences, Novosibirsk 630090, Russian Federation
| | - Sergey P Verevkin
- Department of Physical Chemistry and Faculty of Interdisciplinary Research, Competence Centre CALOR, University of Rostock, 18059 Rostock, Germany.,Department of Physical Chemistry, Kazan Federal University, Kazan 420008, Russian Federation
| |
Collapse
|
14
|
Bakowies D. ATOMIC-2 Protocol for Thermochemistry. J Chem Theory Comput 2022; 18:4142-4163. [PMID: 35658473 DOI: 10.1021/acs.jctc.1c01272] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
ATOMIC is a midlevel thermochemistry protocol that uses Pople's concept of bond separation reactions (BSRs) as a theoretical framework to reduce computational demands in the evaluation of atomization energies and enthalpies of formation. Various composite models are available that approximate bond separation energies at the complete-basis-set limit of all-electron CCSD(T), each balancing computational cost with achievable accuracy. Evaluated energies are then combined with very high-level, precomputed atomization energies of all auxiliary molecules appearing in the BSR to obtain the atomization energy of the molecule under study. ATOMIC-2 is a new version of the protocol that retains the overall concept and all previously defined composite models but improves on ATOMIC-1 in various other ways: Geometry optimization and zero-point-energy evaluation are performed at the density functional level (PBE0-D3/6-311G(d)), which shows significant computational savings and better accuracy than the previously employed RI-MP2/cc-pVTZ. The BSR framework is improved, using more accurate complete-basis-set (CBS) extrapolations toward the Full CI limit for the atomization energies of all auxiliary molecules. Finally, and most importantly, an error and uncertainty model termed ATOMIC-2um is added that estimates average bias and uncertainty for each of the atomization energy contributions that arise from the simplified treatment of some contributions to bond separation energies (CCSD(T)) and the neglect of others (such as higher order, scalar relativistic, or diagonal Born-Oppenheimer corrections) or from residual error in the energies of auxiliary molecules. Large and diverse benchmarks including up to 1179 molecules are used to evaluate necessary reference data and to correlate the observed error for each of the contributions with appropriate proxies that are available without additional quantum-chemical calculations for a particular molecule and represent its size and type. The implementation of ATOMIC-2 considers neutral, closed-shell molecules containing H, C, N, O, and F atoms; compared to ATOMIC-1, the framework has been extended to cover a few challenging but rare bond topologies. In comparison to highly accurate reference data for 184 molecules taken from the ATcT database (V. 1.122r), regular ATOMIC-2 shows noticeable underbinding, but the bias-corrected protocol ATOMIC-2um is found to be more accurate than either ATOMIC-1 or standard Gaussian-4 theory, and the uncertainty model is consistent with statistics of actually observed errors. Problems arising from ambiguous or challenging Lewis-valence structures defining BSRs are discussed, and computational efficiency is demonstrated. Computer code is made available to perform ATOMIC-2um analyses.
Collapse
Affiliation(s)
- Dirk Bakowies
- Institute of Physical Chemistry, Department of Chemistry, University of Basel, Klingelbergstraße 80, CH 4056 Basel, Switzerland
| |
Collapse
|
15
|
Pernot P. The long road to calibrated prediction uncertainty in computational chemistry. J Chem Phys 2022; 156:114109. [DOI: 10.1063/5.0084302] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023] Open
Abstract
Uncertainty quantification (UQ) in computational chemistry (CC) is still in its infancy. Very few CC methods are designed to provide a confidence level on their predictions, and most users still rely improperly on the mean absolute error as an accuracy metric. The development of reliable UQ methods is essential, notably for CC to be used confidently in industrial processes. A review of the CC-UQ literature shows that there is no common standard procedure to report or validate prediction uncertainty. I consider here analysis tools using concepts (calibration and sharpness) developed in meteorology and machine learning for the validation of probabilistic forecasters. These tools are adapted to CC-UQ and applied to datasets of prediction uncertainties provided by composite methods, Bayesian ensembles methods, and machine learning and a posteriori statistical methods.
Collapse
Affiliation(s)
- Pascal Pernot
- Institut de Chimie Physique, UMR8000 CNRS, Université Paris-Saclay, 91405 Orsay, France
| |
Collapse
|
16
|
Vuckovic S. Quantification of Geometric Errors Made Simple: Application to Main-Group Molecular Structures. J Phys Chem A 2022; 126:1300-1311. [PMID: 35144382 DOI: 10.1021/acs.jpca.1c10688] [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/30/2022]
Abstract
Nearly all electronic structure simulations begin with obtaining approximate geometries, making a systematic quantification of errors in approximate molecular structures of key importance. Recently, the geometric energy offset (GEO) framework based on a single and natural measure for quantifying and analyzing these errors has been proposed ( J. Phys. Chem. Lett. 2020, 11, 99579964). An accurate and far less costly approximation to GEO is utilized here to readily quantify errors in main-group structures and analyze them in a chemically intuitive way. The use of semiexperimental geometries as a reference further simplifies the analysis. The analysis reveals new insights into the geometric performance of methods, their rankings, as well as patterns across different classes of methods and basis sets that arise from the analysis.
Collapse
Affiliation(s)
- Stefan Vuckovic
- Institute for Microelectronics and Microsystems (CNR-IMM), Via Monteroni, Campus Unisalento, 73100 Lecce, Italy.,Department of Chemistry & Pharmaceutical Sciences and Amsterdam Institute of Molecular and Life Sciences (AIMMS), Faculty of Science, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands
| |
Collapse
|
17
|
Schell J, Yang K, Glaser R. Computational Investigation of the Thermochemistry of the CO 2 Capture Reaction by Ethylamine, Propylamine, and Butylamine in Aqueous Solution Considering the Full Conformational Space via Boltzmann Statistics. J Phys Chem A 2021; 125:9578-9593. [PMID: 34714081 DOI: 10.1021/acs.jpca.1c06294] [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
Rubisco is the enzyme responsible for CO2 fixation in nature, and it is activated by CO2 addition to the amine group of its lysine 201 side chain. We are designing rubisco-based biomimetic systems for reversible CO2 capture from ambient air. The oligopeptide biomimetic capture systems are employed in aqueous solution. To provide a solid foundation for the experimental solution-phase studies of the CO2 capture reaction, we report here the results of computational studies of the thermodynamics of CO2 capture by small alkylamines in aqueous solution. We studied CO2 addition to methyl-, ethyl-, propyl-, and butylamine with the consideration of the full conformational space for the amine and the corresponding carbamic acids and with the application of an accurate solvation model for the potential energy surface analyses. The reaction energies of the carbamylation reactions were determined based on just the most stable structures (MSS) and based on the ensemble energies computed with the Boltzmann distribution (BD), and it is found that ΔGBD ≈ ΔGMSS. The effect of the proper accounting for the molecular translational entropies in solution with the Wertz approach are much more significant, and the free energy of the capture reactions ΔWGBD is more negative by 2.9 kcal/mol. Further accounting for volume effects in solution results in our best estimates for the reaction energies of the carbamylation reactions of ΔWABD = -5.4 kcal/mol. The overall difference is ΔGBD - ΔWABD = 2.4 kcal/mol for butylamine carbamylation. The full conformational space analyses inform about the conformational isomerizations of carbamic acids, and we determined the relevant rotational profiles and their transition-state structures. Our detailed studies emphasize that, more generally, solution-phase reaction energies should be evaluated with the Helmholtz free energy and can be affected substantially by solution effects on translational entropies.
Collapse
Affiliation(s)
- Joseph Schell
- Department of Chemistry, Missouri University of Science & Technology, Rolla, Missouri 65409, United States.,Department of Chemistry, University of Missouri, Columbia, Missouri 65211, United States
| | - Kaidi Yang
- Department of Chemistry, University of Missouri, Columbia, Missouri 65211, United States
| | - Rainer Glaser
- Department of Chemistry, Missouri University of Science & Technology, Rolla, Missouri 65409, United States
| |
Collapse
|