1
|
Gray M, Herbert JM. Assessing the domain-based local pair natural orbital (DLPNO) approximation for non-covalent interactions in sizable supramolecular complexes. J Chem Phys 2024; 161:054114. [PMID: 39105555 PMCID: PMC11305816 DOI: 10.1063/5.0206533] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2024] [Accepted: 07/18/2024] [Indexed: 08/07/2024] Open
Abstract
The titular domain-based local pair natural orbital (DLPNO) approximation is the most widely used method for extending correlated wave function models to large molecular systems, yet its fidelity for intermolecular interaction energies in large supramolecular complexes has not been thoroughly vetted. Non-covalent interactions are sensitive to tails of the electron density and involve nonlocal dispersion that is discarded or approximated if the screening of pair natural orbitals (PNOs) is too aggressive. Meanwhile, the accuracy of the DLPNO approximation is known to deteriorate as molecular size increases. Here, we test the DLPNO approximation at the level of second-order Møller-Plesset perturbation theory (MP2) and coupled-cluster theory with singles, doubles, and perturbative triples [CCSD(T)] for a variety of large supramolecular complexes. DLPNO-MP2 interaction energies are within 3% of canonical values for small dimers with ≲10 heavy atoms, but for larger systems, the DLPNO approximation is often quite poor unless the results are extrapolated to the canonical limit where the threshold for discarding PNOs is taken to zero. Counterpoise correction proves to be essential in reducing errors with respect to canonical results. For a sequence of nanoscale graphene dimers up to (C96H24)2, extrapolated DLPNO-MP2 interaction energies agree with canonical values to within 1%, independent of system size, provided that the basis set does not contain diffuse functions; these cause the DLPNO approximation to behave erratically, such that results cannot be extrapolated in a meaningful way. DLPNO-CCSD(T) calculations are typically performed using looser PNO thresholds as compared to DLPNO-MP2, but this significantly impacts accuracy for large supramolecular complexes. Standard DLPNO-CCSD(T) settings afford errors of 2-6 kcal/mol for dimers involving coronene (C24H12) and circumcoronene (C54H18), even at the DLPNO-CCSD(T1) level.
Collapse
Affiliation(s)
- Montgomery Gray
- Department of Chemistry & Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - John M. Herbert
- Department of Chemistry & Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| |
Collapse
|
2
|
Liu S, Yang Y, Song Q, Liu Z, Lu Y, Wang Z, Sivaguru P, Bi X. Tunable molecular editing of indoles with fluoroalkyl carbenes. Nat Chem 2024; 16:988-997. [PMID: 38443494 DOI: 10.1038/s41557-024-01468-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2023] [Accepted: 01/31/2024] [Indexed: 03/07/2024]
Abstract
Building molecular complexity from simple feedstocks through precise peripheral and skeletal modifications is central to modern organic synthesis. Nevertheless, a controllable strategy through which both the core skeleton and the periphery of an aromatic heterocycle can be modified with a common substrate remains elusive, despite its potential to maximize structural diversity and applications. Here we report a carbene-initiated chemodivergent molecular editing of indoles that allows both skeletal and peripheral editing by trapping an electrophilic fluoroalkyl carbene generated in situ from fluoroalkyl N-triftosylhydrazones. A variety of fluorine-containing N-heterocyclic scaffolds have been efficiently achieved through tunable chemoselective editing reactions at the skeleton or periphery of indoles, including one-carbon insertion, C3 gem-difluoroolefination, tandem cyclopropanation and N1 gem-difluoroolefination, and cyclopropanation. The power of this chemodivergent molecular editing strategy has been highlighted through the modification of the skeleton or periphery of natural products in a controllable and chemoselective manner. The reaction mechanism and origins of the chemo- and regioselectivity have been probed by both experimental and theoretical methods.
Collapse
Affiliation(s)
- Shaopeng Liu
- Department of Chemistry, Northeast Normal University, Changchun, China
| | - Yong Yang
- Department of Chemistry, Northeast Normal University, Changchun, China
| | - Qingmin Song
- Department of Chemistry, Northeast Normal University, Changchun, China
| | - Zhaohong Liu
- Department of Chemistry, Northeast Normal University, Changchun, China.
| | - Ying Lu
- Department of Chemistry, Northeast Normal University, Changchun, China
| | - Zhanjing Wang
- Department of Chemistry, Northeast Normal University, Changchun, China
| | | | - Xihe Bi
- Department of Chemistry, Northeast Normal University, Changchun, China.
- State Key Laboratory of Elemento-Organic Chemistry, Nankai University, Tianjin, China.
| |
Collapse
|
3
|
Ramek M, Sabolović J. Structure prediction of physiological bis(amino acidato)copper(II) species in aqueous solution: The copper(II) compounds with l-glutamine and l-histidine. J Inorg Biochem 2024; 251:112430. [PMID: 38006660 DOI: 10.1016/j.jinorgbio.2023.112430] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2023] [Revised: 10/27/2023] [Accepted: 11/13/2023] [Indexed: 11/27/2023]
Abstract
Neutral (l-histidinato)(l-glutaminato)copper(II) [Cu(His)(Gln)] has been established as the most abundant ternary copper(II) amino acid compound of the exchangeable copper(II) pool in blood plasma. The experimental studies of Cu(His)(Gln) and bis(glutaminato)copper(II) [Cu(Gln)2] in solutions did not specify their complete geometries. To determine the geometries, this paper investigates the conformers, energy landscapes, and a structure-magnetic parameters relation of Cu(Gln)2 and Cu(His)(Gln) by the density functional theory (DFT) calculations. We assume a glycine-like coordination of Gln (other coordination patterns are dismissed because of steric reasons), and three His in-plane copper(II) binding modes. The conformational analyses are performed in the gas phase and implicitly modeled aqueous solution. The reliability of the DFT relative electronic and Gibbs free energies of the Cu(His)(Gln) conformers is confirmed by benchmarking against the corresponding energies obtained by the domain-based local pair natural orbital coupled-cluster method with singles, doubles, and perturbative triples [DLPNO-CCSD(T)]. Several cis- and trans-Cu(His)(Gln) conformers with His in the histaminate-like and glycine-like modes have low Gibbs free energies, and the greatest estimated metal-binding affinities. The DFT-calculated magnetic parameters of the low-energy conformers reproduce best the experimental electron paramagnetic resonance parameters measured in aqueous solutions for trans- and cis-Cu(Gln)2 conformers having two oxygen atoms (either from Gln or water molecules) at the apical positions, and Cu(His)(Gln) conformers having His in the histaminate-like mode with an apically placed carboxylato oxygen atom. The predicted conformational flexibility of His‑copper(II)-amino acid compounds may be connected with their physiological abundance, and the role in copper(II) exchange reactions in blood plasma.
Collapse
Affiliation(s)
- Michael Ramek
- Institute of Physical and Theoretical Chemistry, Graz University of Technology, Stremayrgasse 9, A-8010 Graz, Austria
| | - Jasmina Sabolović
- Institute for Medical Research and Occupational Health, Ksaverska cesta 2, HR-10000 Zagreb, Croatia.
| |
Collapse
|
4
|
Wappett D, Goerigk L. Benchmarking Density Functional Theory Methods for Metalloenzyme Reactions: The Introduction of the MME55 Set. J Chem Theory Comput 2023; 19:8365-8383. [PMID: 37943578 PMCID: PMC10688432 DOI: 10.1021/acs.jctc.3c00558] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/25/2023] [Revised: 10/13/2023] [Accepted: 10/16/2023] [Indexed: 11/10/2023]
Abstract
We present a new benchmark set of metalloenzyme model reaction energies and barrier heights that we call MME55. The set contains 10 different enzymes, representing eight transition metals, both open and closed shell systems, and system sizes of up to 116 atoms. We use four DLPNO-CCSD(T)-based approaches to calculate reference values against which we then benchmark the performance of a range of density functional approximations with and without dispersion corrections. Dispersion corrections improve the results across the board, and triple-ζ basis sets provide the best balance of efficiency and accuracy. Jacob's ladder is reproduced for the whole set based on averaged mean absolute (percent) deviations, with the double hybrids SOS0-PBE0-2-D3(BJ) and revDOD-PBEP86-D4 standing out as the most accurate methods for the MME55 set. The range-separated hybrids ωB97M-V and ωB97X-V also perform well here and can be recommended as a reliable compromise between accuracy and efficiency; they have already been shown to be robust across many other types of chemical problems, as well. Despite the popularity of B3LYP in computational enzymology, it is not a strong performer on our benchmark set, and we discourage its use for enzyme energetics.
Collapse
Affiliation(s)
- Dominique
A. Wappett
- School of Chemistry, The University
of Melbourne, Melbourne, Victoria 3010, Australia
| | - Lars Goerigk
- School of Chemistry, The University
of Melbourne, Melbourne, Victoria 3010, Australia
| |
Collapse
|
5
|
Buvaylo EA, Nesterova OV, Goreshnik EA, Vyshniakova HV, Petrusenko SR, Nesterov DS. Supramolecular Diversity, Theoretical Investigation and Antibacterial Activity of Cu, Co and Cd Complexes Based on the Tridentate N,N,O-Schiff Base Ligand Formed In Situ. MOLECULES (BASEL, SWITZERLAND) 2022; 27:molecules27238233. [PMID: 36500325 PMCID: PMC9740120 DOI: 10.3390/molecules27238233] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/29/2022] [Revised: 11/18/2022] [Accepted: 11/22/2022] [Indexed: 11/29/2022]
Abstract
The four new complexes, [Cu(HL1)(L2)Cl] (1), [Cu(HL1)(L1)]∙Cl∙2H2O (2), [Co(L1)2]∙Cl (3) and [Cd(HL1)I2]∙dmso (4), have been prepared by one-pot reactions of the respective chloride or iodide metal salt with a non-aqueous solution of the polydentate Schiff base, HL1, resulted from in situ condensation of benzhydrazide and 2-pyridinecarboxaldehyde, while a ligand HL2, in case of 1, has been formed due to the oxidation of 2-pyridinecarboxaldehyde under reaction conditions. The crystallographic analysis revealed that the molecular building units in 1-4 are linked together into complex structures by hydrogen bonding, resulting in 1D, 2D and 3D supramolecular architectures for 1, 2 and 4, respectively, and the supramolecular trimer for 3. The electronic structures of 1-4 were investigated by the DFT theoretical calculations. The non-covalent interactions in the crystal structures of 1-4 were studied by means of the Hirshfeld surface analysis and the QTAIM theory with a special focus on the C-H⋯Cl bonding. From the DFT/DLPNO-CCSD(T) calculations, using a series of charged model {R3C-H}0⋯Cl- assemblies, we propose linear regressions for assessment of the interaction enthalpy (ΔH, kcal mol-1) and the binding energy (BE, kcal mol-1) between {R3C-H}0 and Cl- sites starting from the electron density at the bond critical point (ρ(rBCP), a.u.): ΔH = -678 × ρ(r) + 3 and BE = -726 × ρ(r) + 4. It was also has been found that compounds 1, 3 and 4 during in vitro screening showed an antibacterial activity toward the nine bacteria species, comprising both Gram-positive and Gram-negative, with MIC values ranging from 156.2 to 625 mg/L. The best results have been obtained against Acinetobacter baumannii MβL.
Collapse
Affiliation(s)
- Elena A. Buvaylo
- Department of Chemistry, Taras Shevchenko National University of Kyiv, Volodymyrska 64/13, 01601 Kyiv, Ukraine
| | - Oksana V. Nesterova
- Centro de Química Estrutural, Institute of Molecular Sciences, Instituto Superior Técnico, Universidade de Lisboa, Av. Rovisco Pais, 1049-001 Lisbon, Portugal
| | - Evgeny A. Goreshnik
- Department of Inorganic Chemistry and Technology, Jožef Stefan Institute, Jamova 39, 1000 Ljubljana, Slovenia
| | - Hanna V. Vyshniakova
- L.V. Gromashevsky Institute of Epidemiology and Infectious Diseases NAMS of Ukraine, M. Amosova 5, 03038 Kyiv, Ukraine
| | - Svitlana R. Petrusenko
- Department of Chemistry, Taras Shevchenko National University of Kyiv, Volodymyrska 64/13, 01601 Kyiv, Ukraine
| | - Dmytro S. Nesterov
- Centro de Química Estrutural, Institute of Molecular Sciences, Instituto Superior Técnico, Universidade de Lisboa, Av. Rovisco Pais, 1049-001 Lisbon, Portugal
- Correspondence:
| |
Collapse
|
6
|
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
|
7
|
Drosou M, Mitsopoulou CA, Pantazis DA. Reconciling Local Coupled Cluster with Multireference Approaches for Transition Metal Spin-State Energetics. J Chem Theory Comput 2022; 18:3538-3548. [PMID: 35582788 PMCID: PMC9202354 DOI: 10.1021/acs.jctc.2c00265] [Citation(s) in RCA: 21] [Impact Index Per Article: 10.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
![]()
Spin-state energetics
of transition metal complexes remain one
of the most challenging targets for electronic structure methods.
Among single-reference wave function approaches, local correlation
approximations to coupled cluster theory, most notably the domain-based
local pair natural orbital (DLPNO) approach, hold the promise of bringing
the accuracy of coupled cluster theory with single, double, and perturbative
triple excitations, CCSD(T), to molecular systems of realistic size
with acceptable computational cost. However, recent studies on spin-state
energetics of iron-containing systems raised doubts about the ability
of the DLPNO approach to adequately and systematically approximate
energetics obtained by the reference-quality complete active space
second-order perturbation theory with coupled-cluster semicore correlation,
CASPT2/CC. Here, we revisit this problem using a diverse set of iron
complexes and examine several aspects of the application of the DLPNO
approach. We show that DLPNO-CCSD(T) can accurately reproduce both
CASPT2/CC and canonical CCSD(T) results if two basic principles are
followed. These include the consistent use of the improved iterative
(T1) versus the semicanonical perturbative triple corrections
and, most importantly, a simple two-point extrapolation to the PNO
space limit. The latter practically eliminates errors arising from
the default truncation of electron-pair correlation spaces and should
be viewed as standard practice in applications of the method to transition
metal spin-state energetics. Our results show that reference-quality
results can be readily achieved with DLPNO-CCSD(T) if these principles
are followed. This is important also in view of the applicability
of the method to larger single-reference systems and multinuclear
clusters, whose treatment of dynamic correlation would be challenging
for multireference-based approaches.
Collapse
Affiliation(s)
- Maria Drosou
- Inorganic Chemistry Laboratory, National and Kapodistrian University of Athens, Panepistimiopolis, Zografou 15771, Greece
| | - Christiana A Mitsopoulou
- Inorganic Chemistry Laboratory, National and Kapodistrian University of Athens, Panepistimiopolis, Zografou 15771, Greece
| | - Dimitrios A Pantazis
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
8
|
Seeger ZL, Izgorodina EI. A DLPNO-CCSD(T) benchmarking study of intermolecular interactions of ionic liquids. J Comput Chem 2022; 43:106-120. [PMID: 34687062 DOI: 10.1002/jcc.26776] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2021] [Revised: 09/14/2021] [Accepted: 09/15/2021] [Indexed: 11/11/2022]
Abstract
The accuracy of correlation energy recovered by coupled cluster single-, double-, and perturbative triple-excitations, CCSD(T), has led to the method being considered the gold standard of computational chemistry. The application of CCSD(T) has been limited to medium-sized molecular systems due to its steep scaling with molecular size. The recent development of alternative domain-based local pair natural orbital coupled-cluster method, DLPNO-CCSD(T), has significantly broadened the range of chemical systems to which CCSD(T) level calculations can be applied. Condensed systems such as ionic liquids (ILs) have a large contribution from London dispersion forces of up to 150 kJ mol-1 in large-scale clusters. Ionic liquids show appreciable charge transfer effects that result in the increased valence orbital delocalization over the entire ionic network, raising the question whether the application of methods based on localized orbitals is reliable for these semi-Coulombic materials. Here the performance of DLPNO-CCSD(T) is validated for the prediction of correlation interaction energies of two data sets incorporating single-ion pairs of protic and aprotic ILs. DLPNO-CCSD(T) produced results within chemical accuracy with tight parameter settings and a non-iterative treatment of triple excitations. To achieve spectroscopic accuracy of 1 kJ mol-1 , especially for hydrogen-bonded ILs and those containing halides, the DLPNO settings had to be increased by two orders of magnitude and include the iterative treatment of triple excitations, resulting in a 2.5-fold increase in computational cost. Two new sets of parameters are put forward to produce the performance of DLPNO-CCSD(T) within chemical and spectroscopic accuracy.
Collapse
Affiliation(s)
- Zoe L Seeger
- School of Chemistry, Monash University, Clayton, Victoria, Australia
| | | |
Collapse
|
9
|
Amsler J, Bernart S, Plessow PN, Studt F. Theoretical investigation of the olefin cycle in H-SSZ-13 for the ethanol-to-olefins process using ab initio calculations and kinetic modeling. Catal Sci Technol 2022. [DOI: 10.1039/d1cy02289j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The formation of the hydrocarbon pool (HCP) in the ethanol-to-olefins (ETO) process catalyzed by H-SSZ-13 is studied in a kinetic model with ab initio computed reaction barriers.
Collapse
Affiliation(s)
- Jonas Amsler
- Institute of Catalysis Research and Technology, Karlsruhe Institute of Technology, Hermann-von-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany
| | - Sarah Bernart
- Institute of Catalysis Research and Technology, Karlsruhe Institute of Technology, Hermann-von-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany
| | - Philipp N. Plessow
- Institute of Catalysis Research and Technology, Karlsruhe Institute of Technology, Hermann-von-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany
| | - Felix Studt
- Institute of Catalysis Research and Technology, Karlsruhe Institute of Technology, Hermann-von-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany
- Institute for Chemical Technology and Polymer Chemistry, Karlsruhe Institute of Technology, Kaiserstr. 12, 76131 Karlsruhe, Germany
| |
Collapse
|
10
|
Kumar P, Escudero D. Computational Protocol to Calculate the Phosphorescence Energy of Pt(II) Complexes: Is the Lowest Triplet Excited State Always Involved in Emission? A Comprehensive Benchmark Study. Inorg Chem 2021; 60:17230-17240. [PMID: 34702026 DOI: 10.1021/acs.inorgchem.1c02562] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Abstract
The reliable calculation of phosphorescence energies of phosphor materials is at the core of designing efficient phosphorescent organic light-emitting diodes (PhOLEDs). Therefore, it is of paramount importance to have a robust computational protocol to perform those calculations in a black-box manner. In this work, we use Domain-Based Local Pair Natural Orbital Coupled Cluster theory with single, double, and perturbative triple excitation (DLPNO-CCSD(T)) calculations to attain the phosphorescence energies of a large pool of Pt(II) complexes. Several approaches to incorporate relativistic effects in our calculations were tested. In addition, we have used the DLPNO-CCSD(T) values (i.e., our best theoretical values) to assess the performance of different flavors of density functional theory including pure, hybrid, meta-hybrid, and range-separated functionals. Among the tested functionals, the M06HF functional provides the best values compared with the DLPNO-CCSD(T) ones, with a mean absolute deviation (MAD) value of 0.14 eV. In its turn, and thanks to the increased accuracy achieved in the calculation of phosphorescence energies, we also demonstrate that not all of the investigated complexes emit from their lowest-lying triplet state (T1). The outlier complexes include different complex photophysical scenarios and both Kasha and anti-Kasha types of complexes. Finally, we provide a general computational protocol to pre-screen whether T1 is actually the emissive state and to accurately calculate the phosphorescence energies of Pt(II) complexes.
Collapse
Affiliation(s)
- Prashant Kumar
- Quantum Chemistry and Physical Chemistry Section, Department of Chemistry, KU Leuven, Celestijnenlaan 200F, B-3001 Leuven, Belgium
| | - Daniel Escudero
- Quantum Chemistry and Physical Chemistry Section, Department of Chemistry, KU Leuven, Celestijnenlaan 200F, B-3001 Leuven, Belgium
| |
Collapse
|
11
|
Understanding water mediated proton migration in conversion of π-bond in olefinic carbon atoms into C–N bond to form β-amino adducts. Tetrahedron 2021. [DOI: 10.1016/j.tet.2021.132482] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
12
|
Maurer LR, Bursch M, Grimme S, Hansen A. Assessing Density Functional Theory for Chemically Relevant Open-Shell Transition Metal Reactions. J Chem Theory Comput 2021; 17:6134-6151. [PMID: 34546754 DOI: 10.1021/acs.jctc.1c00659] [Citation(s) in RCA: 55] [Impact Index Per Article: 18.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Due to the principle lack of systematic improvement possibilities of density functional theory, careful assessment of the performance of density functional approximations (DFAs) on well-designed benchmark sets, for example, for reaction energies and barrier heights, is crucial. While main-group chemistry is well covered by several available sets, benchmark data for transition metal chemistry is sparse. This is especially the case for larger, chemically relevant molecules. Addressing this issue, we recently introduced the MOR41 benchmark which covers chemically relevant reactions of closed-shell complexes. In this work, we extend these efforts to single-reference open-shell systems and introduce the "reactions of open-shell single-reference transition metal complexes" (ROST61) benchmark set. ROST61 includes accurate coupled-cluster reference values for 61 reaction energies with a mean reaction energy of -42.8 kcal mol-1. Complexes with 13-93 atoms covering 20 d-block elements are included, but due to the restriction to single-reference open-shell systems, important elements such as iron or platinum could not be taken into account, or only to a small extent. We assess the performance of 31 DFAs in combination with three London dispersion (LD) correction schemes. Further, DFT-based composite methods, MP2, and a few semiempirical quantum chemical methods are evaluated. Consistent with the results for the MOR41 closed-shell benchmark, we find that the ordering of DFAs according to Jacob's ladder is preserved and that adding an LD correction is crucial, clearly improving almost all tested methods. The recently introduced r2SCAN-3c composite method stands out with a remarkable mean absolute deviation (MAD) of only 2.9 kcal mol-1, which is surpassed only by hybrid DFAs with low amounts of Fock exchange (e.g., 2.3 kcal mol-1 for TPSS0-D4/def2-QZVPP) and double-hybrid (DH) DFAs but at a significantly higher computational cost. The lowest MAD of only 1.6 kcal mol-1 is obtained with the DH DFA PWPB95-D4 in the def2-QZVPP basis set approaching the estimated accuracy of the reference method. Overall, the ROST61 set adds important reference data to a sparsely sampled but practically relevant area of chemistry. At this point, it provides valuable orientation for the application and development of new DFAs and electronic structure methods in general.
Collapse
Affiliation(s)
- Leonard R Maurer
- Mulliken Center for Theoretical Chemistry, Institute for Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115 Bonn, Germany
| | - Markus Bursch
- Mulliken Center for Theoretical Chemistry, Institute for Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115 Bonn, Germany
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Institute for Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115 Bonn, Germany
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry, Institute for Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115 Bonn, Germany
| |
Collapse
|
13
|
Blaško M, Pašteka LF, Urban M. DFT Functionals for Modeling of Polyethylene Chains Cross-Linked by Metal Atoms. DLPNO-CCSD(T) Benchmark Calculations. J Phys Chem A 2021; 125:7382-7395. [PMID: 34428051 DOI: 10.1021/acs.jpca.1c04793] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Density functional theory (DFT) functionals for calculations of binding energies (BEs) of the polyethylene (PE) chains cross-linked by selected metal atoms (M) are benchmarked against DLPNO-CCSD(T) and DLPNO-CCSD(T1) data. PEX-M-PEX complexes as compared with plain parallel PEX···PEX chains with X = 3-9 carbon atoms are model species characterized by a cooperative effect of covalent C-M-C bonds and interchain dispersion interactions. The accuracy of DLPNO-CC methods was assessed by a comparison of BEs with the canonical CCSD(T) results for small PE3-M-PE3 complexes. Functionals for PEX···PEX and closed-shell PEX-M-PEX complexes (M = Be, Mg, Zn) were benchmarked against DLPNO-CCSD(T) BEs; open-shell complexes (M = Li, Ag, Au) were benchmarked against the DLPNO-CCSD(T1) method with iterative triples. Three dispersion corrections were combined with 25 DFT functionals for calculations of BEs with respect to PEX-M and PEX fragments employing def2-TZVPP and def2-QZVPP basis sets. Accuracy to within 5% for the closed-shell PEX-M-PEX complexes was achieved with five functionals. Less accurate are functionals for the open-shell PEX-M-PEX complexes; only two functionals deviate by less than 15% from DLPNO-CCSD(T1). Particularly problematic were PEX-Li-PEX complexes. A reasonable overall performance across all complexes in terms of the mean absolute percentage error is found for the range-separated hybrid functionals ωB97X-D3 and CAM-B3LYP/D3(BJ)-ABC.
Collapse
Affiliation(s)
- Martin Blaško
- FunGlass, A. Dubček University of Trenčín, Študentská 2, 911 50 Trenčín, Slovakia
| | - Lukáš F Pašteka
- Department of Physical and Theoretical Chemistry, Faculty of Natural Sciences, Comenius University in Bratislava, Ilkovičova 6, 84215 Bratislava, Slovakia
| | - Miroslav Urban
- Department of Physical and Theoretical Chemistry, Faculty of Natural Sciences, Comenius University in Bratislava, Ilkovičova 6, 84215 Bratislava, Slovakia
| |
Collapse
|
14
|
Zheng K, Li D, Jiang L, Li X, Xie C, Feng L, Qin J, Qian S, Pang Q. Revisiting stacking interactions in tetrathiafulvalene and selected derivatives using tight-binding quantum chemical calculations and local coupled-cluster method. ACTA CRYSTALLOGRAPHICA SECTION B, STRUCTURAL SCIENCE, CRYSTAL ENGINEERING AND MATERIALS 2021; 77:311-320. [PMID: 34096512 DOI: 10.1107/s2052520621003085] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/29/2020] [Accepted: 03/23/2021] [Indexed: 06/12/2023]
Abstract
The engineering of supramolecular architectures needs accurate descriptions of the intermolecular interactions in crystal structures. Tetrathiafulvalene (TTF) is an effective building block used in the construction of promising functional materials. The parallel packing of the neutral TTF-TTF system was studied previously using the high-level quantum chemical method, advancing it as a valuable model system. The recently developed tight-binding quantum chemical method GFN2-xTB and local coupled-cluster method DLPNO-CCSD(T) were used to investigate the stacking interactions of TTF and selected derivatives deposited in the Cambridge Structural Database. Using the interaction energy of the TTF-TTF dimer calculated at the CCSD(T)/CBS level as the reference, the accuracies of the two methods are investigated. The energy decomposition analysis within the DLPNO-CCSD(T) framework reveals the importance of dispersion interaction in the TTF-related stacking systems. The dispersion interaction density plot vividly shows the magnitude and distribution of the dispersion interaction, providing a revealing insight into the stacking interactions in crystal structures. The results show that the GFN2-xTB and DLPNO-CCSD(T) methods could achieve accuracy at an affordable computational cost, which would be valuable in understanding the nature of parallel stacking in supramolecular systems.
Collapse
Affiliation(s)
- Kang Zheng
- School of Life Sciences, Shandong University of Technology, 266 Xincun West Road, Zibo, Shandong 255049, People's Republic of China
| | - Danping Li
- School of Life Sciences, Shandong University of Technology, 266 Xincun West Road, Zibo, Shandong 255049, People's Republic of China
| | - Liu Jiang
- School of Life Sciences, Shandong University of Technology, 266 Xincun West Road, Zibo, Shandong 255049, People's Republic of China
| | - Xiaowei Li
- School of Materials Science and Engineering, Shandong University of Technology, 266 Xincun West Road, Zibo, Shandong 255049, People's Republic of China
| | - Changjian Xie
- School of Life Sciences, Shandong University of Technology, 266 Xincun West Road, Zibo, Shandong 255049, People's Republic of China
| | - Ling Feng
- School of Life Sciences, Shandong University of Technology, 266 Xincun West Road, Zibo, Shandong 255049, People's Republic of China
| | - Jie Qin
- School of Life Sciences, Shandong University of Technology, 266 Xincun West Road, Zibo, Shandong 255049, People's Republic of China
| | - Shaosong Qian
- School of Life Sciences, Shandong University of Technology, 266 Xincun West Road, Zibo, Shandong 255049, People's Republic of China
| | - Qiuxiang Pang
- School of Life Sciences, Shandong University of Technology, 266 Xincun West Road, Zibo, Shandong 255049, People's Republic of China
| |
Collapse
|
15
|
Zhang Z, Yang Z, Pu L, Chen X, Li Y, Wang J, Zhao L, King RB. Mechanism for the Reaction of White Phosphorus with Cp 2Cr 2(CO) 6 Leading Ultimately to the Triple-Decker Sandwich Cp 2Cr 2(μ-η 5,η 5-P 5): A Theoretical Study. Inorg Chem 2021; 60:5955-5968. [PMID: 33834774 DOI: 10.1021/acs.inorgchem.1c00382] [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
The experimentally known reaction of Cp2Cr2(CO)6 with white phosphorus (P4) to give CpCr(CO)2(η3-P3), Cp2Cr2(CO)4(μ-η,2η2-P2), and the triple-decker sandwich Cp2Cr2(μ-η,5η5-P5) is of interest since the P4 reactant having a tetrahedral cluster of four phosphorus atoms is converted to products having P2, P3, and P5 ligands. The mechanism of this obviously complicated reaction can be dissected into three stages using a coupled cluster theoretical method that has been benchmarked with the P2, Mn(CO)5, and CpCr(CO)3 dimerization processes. The first stage of the Cp2Cr2(CO)6/P4 reaction mechanism generates the unsaturated singlet intermediate Cp2Cr2(CO)5 that combines with the P4 reactant. Decarbonylation of the resulting Cp2Cr2(CO)5(P4) complex provides a singlet tetracarbonyl readily fragmenting into the stable triphosphacyclopropenyl complex CpCr(CO)2(η3-P3) and the chromium phosphide CpCr(CO)2(P). The isomeric triplet tetracarbonyl Cp2Cr2(CO)4(P4), readily fragments into CpCr(CO)2(η2-P2), which can generate the stable diphosphaacetylene complex Cp2Cr2(CO)4(η,2η2-P2) as well as the pentamer [CpCr(CO)2]5(P10). Combination of the coordinately unsaturated CpCr(CO)(η3-P3) with CpCr(CO)2(η2-P2) can lead to a ring expansion. This generates the P5 pentagonal ligand in a Cp2Cr2(CO)3(P5) precursor to the experimentally observed carbonyl-free triple-decker sandwich Cp2Cr2(μ-η,5η5-P5) after three successive decarbonylations.
Collapse
Affiliation(s)
- Zhong Zhang
- College of Chemistry & Pharmacy, Northwest A&F University, Yangling, Shaanxi 712100, P. R. China
| | - Zhipeng Yang
- College of Chemistry & Pharmacy, Northwest A&F University, Yangling, Shaanxi 712100, P. R. China
| | - Liang Pu
- College of Chemistry & Pharmacy, Northwest A&F University, Yangling, Shaanxi 712100, P. R. China
| | - Xian Chen
- College of Chemistry & Pharmacy, Northwest A&F University, Yangling, Shaanxi 712100, P. R. China
| | - Yun Li
- College of Chemistry & Pharmacy, Northwest A&F University, Yangling, Shaanxi 712100, P. R. China
| | - Jianping Wang
- Xi'an Key Laboratory of Advanced Photo-electronics Materials and Energy Conversion Device, Xijing University, Xi'an, Shaanxi 710123, P. R. China
| | - Lingzhi Zhao
- SCNU Qingyuan Institute of Science and Technology Innovation Co., Ltd., Qingyuan 511517, P. R. China
| | - R Bruce King
- Department of Chemistry and Center for Computational Chemistry, University of Georgia, Athens, Georgia 30602, United States
| |
Collapse
|
16
|
Khan A, Goepel M, Kubas A, Łomot D, Lisowski W, Lisovytskiy D, Nowicka A, Colmenares JC, Gläser R. Selective Oxidation of 5-Hydroxymethylfurfural to 2,5-Diformylfuran by Visible Light-Driven Photocatalysis over In Situ Substrate-Sensitized Titania. CHEMSUSCHEM 2021; 14:1351-1362. [PMID: 33453092 PMCID: PMC7986172 DOI: 10.1002/cssc.202002687] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/19/2020] [Revised: 12/30/2020] [Indexed: 06/12/2023]
Abstract
Solar energy-driven processes for biomass valorization are priority for the growing industrialized society. To address this challenge, efficient visible light-active photocatalyst for the selective oxidation of biomass-derived platform chemical is highly desirable. Herein, selective oxidation of 5-hydroxymethylfurfural (HMF) to 2,5-diformylfuran (DFF) was achieved by visible light-driven photocatalysis over titania. Pristine titania is photocatalytically inactive under visible light, so an unconventional approach was employed for the visible light (λ=515 nm) sensitization of titania via a formation of a visible light-absorbing complex of HMF (substrate) on the titania surface. Surface-complexation of HMF on titania mediated ligand-to-metal charge transfer (LMCT) under visible light, which efficiently catalyzed the oxidation of HMF to DFF. A high DFF selectivity of 87 % was achieved with 59 % HMF conversion after 4 h of illumination. The apparent quantum yield obtained for DFF production was calculated to be 6.3 %. It was proposed that the dissociative interaction of hydroxyl groups of HMF and the titania surface is responsible for the surface-complex formation. When the hydroxyl groups of titania were modified via surface-fluorination or calcination the oxidation of HMF was inhibited under visible light, signifying that hydroxyl groups are decisive for photocatalytic activity.
Collapse
Affiliation(s)
- Ayesha Khan
- Institute of Physical ChemistryPolish Academy of SciencesWarsaw01-224Poland
| | - Michael Goepel
- Institute of Chemical TechnologyLeipzig UniversityLeipzig04103Germany
| | - Adam Kubas
- Institute of Physical ChemistryPolish Academy of SciencesWarsaw01-224Poland
| | - Dariusz Łomot
- Institute of Physical ChemistryPolish Academy of SciencesWarsaw01-224Poland
| | - Wojciech Lisowski
- Institute of Physical ChemistryPolish Academy of SciencesWarsaw01-224Poland
| | - Dmytro Lisovytskiy
- Institute of Physical ChemistryPolish Academy of SciencesWarsaw01-224Poland
| | - Ariadna Nowicka
- Institute of Physical ChemistryPolish Academy of SciencesWarsaw01-224Poland
| | | | - Roger Gläser
- Institute of Chemical TechnologyLeipzig UniversityLeipzig04103Germany
| |
Collapse
|
17
|
Bursch M, Hansen A, Pracht P, Kohn JT, Grimme S. Theoretical study on conformational energies of transition metal complexes. Phys Chem Chem Phys 2021; 23:287-299. [PMID: 33336657 DOI: 10.1039/d0cp04696e] [Citation(s) in RCA: 43] [Impact Index Per Article: 14.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Abstract
Conformational energies are an important chemical property for which a performance assessment of theoretical methods is mandatory. Existing benchmark sets are often limited to biochemical or main group element containing molecules, while organometallic systems are generally less studied. A key problem herein is to routinely generate conformers for these molecules due to their complexity and manifold of possible coordination patterns. In this study we used our recently published CREST protocol [Pracht et al., Phys. Chem. Chem. Phys., 2020, 22, 7169-7192] to generate conformer ensembles for a variety of 40 challenging transition metal containing molecules, which were then used to form a comprehensive conformational energy benchmark set termed TMCONF40. Several low-cost semiempirical, density functional theory (DFT) and force-field methods were compared to high level DLPNO-CCSD(T1) and double-hybrid DFT reference values. Close attention was paid to the energetic ordering of the conformers in the statistical evaluation. With respect to the double-hybrid references, both tested low-cost composite DFT methods produce high Pearson correlation coefficients of rp,mean,B97-3c//B97-3c = 0.922 and rp,mean,PBEh-3c//B97-3c = 0.890, with mean absolute deviations close to or below 1 kcal mol-1. This good performance also holds for a comparison to DLPNO-CCSD(T1) reference energies for a smaller subset termed TMCONF5. Based on DFT geometries, the GFNn-xTB methods yield reasonable Pearson correlation coefficients of rp,mean,GFN1-xTB//B97-3c = 0.617 (MADmean = 2.15 kcal mol-1) and rp,mean,GFN2-xTB//B97-3c = 0.567 (MADmean = 2.68 kcal mol-1), outperforming the widely used PMx methods on the TMCONF40 test set. Employing the low-cost composite DFT method B97-3c on GFN2-xTB geometries yields an slightly improved correlation of rp,mean,B97-3c//GFN2-xTB = 0.632. Furthermore, for 68% of the investigated complexes at least one low-energy conformer was found that is more stable than the respective crystal structure conformation, which signals the importance of conformational studies. General recommendations for the application of the CREST protocol and DFT methods for transition metal conformational energies are given.
Collapse
Affiliation(s)
- Markus Bursch
- Mulliken Center for Theoretical Chemistry, Universität Bonn, Beringstr. 4, 53115 Bonn, Germany.
| | | | | | | | | |
Collapse
|
18
|
Fečík M, Plessow PN, Studt F. Theoretical investigation of the side-chain mechanism of the MTO process over H-SSZ-13 using DFT and ab initio calculations. Catal Sci Technol 2021. [DOI: 10.1039/d1cy00433f] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The side-chain mechanism of the methanol-to-olefins process over the H-SSZ-13 acidic zeolite was investigated using periodic density functional theory with corrections from highly accurate ab intio calculations on large cluster models.
Collapse
Affiliation(s)
- Michal Fečík
- Institute of Catalysis Research and Technology
- Karlsruhe Institute of Technology
- 76344 Eggenstein-Leopoldshafen
- Germany
| | - Philipp N. Plessow
- Institute of Catalysis Research and Technology
- Karlsruhe Institute of Technology
- 76344 Eggenstein-Leopoldshafen
- Germany
| | - Felix Studt
- Institute of Catalysis Research and Technology
- Karlsruhe Institute of Technology
- 76344 Eggenstein-Leopoldshafen
- Germany
- Institute for Chemical Technology and Polymer Chemistry
| |
Collapse
|
19
|
Itkis D, Cavallo L, Yashina LV, Minenkov Y. Ambiguities in solvation free energies from cluster-continuum quasichemical theory: lithium cation in protic and aprotic solvents. Phys Chem Chem Phys 2021; 23:16077-16088. [PMID: 34291782 DOI: 10.1039/d1cp01454d] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Gibbs free energies for Li+ solvation in water, methanol, acetonitrile, DMSO, dimethylacetamide, dimethoxyethane, dimethylformamide, gamma-butyrolactone, pyridine, and sulfolane have been calculated using the cluster-continuum quasichemical theory. With n independent solvent molecules S initial state forming the "monomer" thermodynamic cycle, Li+ solvation free energies are found to be on average 14 kcal mol-1 more positive compared to those from the "cluster" thermodynamic cycle where the initial state is the cluster Sn. We ascribe the inconsistency between the "monomer" and "cluster" cycles mainly to the incorrectly predicted solvation free energies of solvent clusters Sn from the SMD and CPCM continuum solvation models, which is in line with the earlier study of Bryantsev et al., J. Phys. Chem. B, 2008, 112, 9709-9719. When experimental-based solvation free energies of individual solvent molecules and solvent clusters are employed, the "monomer" and "cluster" cycles result in identical numbers. The best overall agreement with experimental-based "bulk" scale lithium cation solvation free energies was obtained for the "monomer" scale, and we recommend this set of values. We expect that further progress in the field is possible if (i) consensus on the accuracy of experimental reference values is achieved; (ii) the most recent continuum solvation models are properly parameterized for all solute-solvent combinations and become widely accessible for testing.
Collapse
Affiliation(s)
- Daniil Itkis
- N.N. Semenov Federal Research Center for Chemical Physics RAS, Kosygina Street 4, 119991 Moscow, Russia.
| | | | | | | |
Collapse
|
20
|
Sorathia K, Tew DP. Basis set extrapolation in pair natural orbital theories. J Chem Phys 2020; 153:174112. [PMID: 33167642 DOI: 10.1063/5.0022077] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present the results of a benchmark study of the effect of Pair Natural Orbital (PNO) truncation errors on the performance of basis set extrapolation. We find that reliable conclusions from the application of Helgaker's extrapolation method are only obtained when using tight PNO thresholds of at least 10-7. The use of looser thresholds introduces a significant risk of observing a false basis set convergence and underestimating the residual basis set errors. We propose an alternative extrapolation approach based on the PNO truncation level that only requires a single basis set and show that it is a viable alternative to hierarchical basis set extrapolation methods.
Collapse
Affiliation(s)
- Kesha Sorathia
- Max Planck Institute for Solid State Research, Heisenbergstrasse 1, 70569 Stuttgart, GermanyUniversity of Oxford, South Parks Road, Oxford OX1 3QZ, United Kingdom
| | - David P Tew
- Max Planck Institute for Solid State Research, Heisenbergstrasse 1, 70569 Stuttgart, Germany
| |
Collapse
|
21
|
Mallick S, Roy B, Kumar P. A comparison of DLPNO-CCSD(T) and CCSD(T) method for the determination of the energetics of hydrogen atom transfer reactions. COMPUT THEOR CHEM 2020. [DOI: 10.1016/j.comptc.2020.112934] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
|
22
|
Fečík M, Plessow PN, Studt F. A Systematic Study of Methylation from Benzene to Hexamethylbenzene in H-SSZ-13 Using Density Functional Theory and Ab Initio Calculations. ACS Catal 2020. [DOI: 10.1021/acscatal.0c02037] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Michal Fečík
- Institute of Catalysis Research and Technology, Karlsruhe Institute of Technology, Hermann-von-Helmholtz Platz 1, Eggenstein-Leopoldshafen 76344, Germany
| | - Philipp N. Plessow
- Institute of Catalysis Research and Technology, Karlsruhe Institute of Technology, Hermann-von-Helmholtz Platz 1, Eggenstein-Leopoldshafen 76344, Germany
| | - Felix Studt
- Institute of Catalysis Research and Technology, Karlsruhe Institute of Technology, Hermann-von-Helmholtz Platz 1, Eggenstein-Leopoldshafen 76344, Germany
- Institute for Chemical Technology and Polymer Chemistry, Karlsruhe Institute of Technology, Engesserstrasse 18, Karlsruhe 76131, Germany
| |
Collapse
|
23
|
Ye HZ, Tran HK, Van Voorhis T. Bootstrap Embedding For Large Molecular Systems. J Chem Theory Comput 2020; 16:5035-5046. [PMID: 32589842 DOI: 10.1021/acs.jctc.0c00438] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
Recent developments in quantum embedding theories have provided attractive approaches to correlated calculations for large systems. In this work, we extend our previous work [J. Chem. Theory Comput. 2019, 15, 4497-4506; J. Phys. Chem. Lett. 2019, 10, 6368-6374] on bootstrap embedding (BE) to enable correlated ab initio calculations at the coupled cluster with singles and doubles (CCSD) level for large molecules. We introduce several new algorithmic developments that significantly reduce the computational cost of BE, while maintaining its accuracy. The resulting implementation scales as O(N3) for the integral transform and O(N) for the CCSD calculation. Numerical results on a series of conjugated molecules suggest that BE with reasonably sized fragments can recover more than 99.5% of the total correlation energy of a full CCSD calculation, while the required computational resources (time and storage) compare favorably to one popular local correlation scheme: domain localized pair natural orbital (DLPNO). The largest BE calculation in this work involves ∼2900 basis functions and can be performed on a single node with 16 CPU cores and 64 GB of memory in a few days. We anticipate that these developments represent an important step toward the application of BE to solve practical problems.
Collapse
Affiliation(s)
- Hong-Zhou Ye
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Henry K Tran
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Troy Van Voorhis
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
24
|
Wu F, Deraedt C, Cornaton Y, Contreras-Garcia J, Boucher M, Karmazin L, Bailly C, Djukic JP. Making Base-Assisted C–H Bond Activation by Cp*Co(III) Effective: A Noncovalent Interaction-Inclusive Theoretical Insight and Experimental Validation. Organometallics 2020. [DOI: 10.1021/acs.organomet.0c00253] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Fule Wu
- Laboratoire de Chimie et Systémique Organométalliques, Institut de Chimie de Strasbourg UMR 7177 CNRS, Université de Strasbourg, 4 rue Blaise Pascal, 67000 Strasbourg, France
| | - Christophe Deraedt
- Laboratoire de Chimie et Systémique Organométalliques, Institut de Chimie de Strasbourg UMR 7177 CNRS, Université de Strasbourg, 4 rue Blaise Pascal, 67000 Strasbourg, France
| | - Yann Cornaton
- Laboratoire de Chimie et Systémique Organométalliques, Institut de Chimie de Strasbourg UMR 7177 CNRS, Université de Strasbourg, 4 rue Blaise Pascal, 67000 Strasbourg, France
| | - Julia Contreras-Garcia
- Laboratoire de Chimie Théorique UMR 7616 CNRS, Sorbonne Université, Site Jussieu, 4 place Jussieu, 75052 Paris cedex, France
| | - Mélanie Boucher
- Laboratoire de Chimie et Systémique Organométalliques, Institut de Chimie de Strasbourg UMR 7177 CNRS, Université de Strasbourg, 4 rue Blaise Pascal, 67000 Strasbourg, France
| | - Lydia Karmazin
- Service de Radiocristallographie, Fédération de Chimie Le Bel FR 2010, Université de Strasbourg, 1 rue Blaise Pascal, 67000 Strasbourg, France
| | - Corinne Bailly
- Service de Radiocristallographie, Fédération de Chimie Le Bel FR 2010, Université de Strasbourg, 1 rue Blaise Pascal, 67000 Strasbourg, France
| | - Jean-Pierre Djukic
- Laboratoire de Chimie et Systémique Organométalliques, Institut de Chimie de Strasbourg UMR 7177 CNRS, Université de Strasbourg, 4 rue Blaise Pascal, 67000 Strasbourg, France
| |
Collapse
|
25
|
Sylvetsky N. Toward Simple, Predictive Understanding of Protein-Ligand Interactions: Electronic Structure Calculations on Torpedo Californica Acetylcholinesterase Join Forces with the Chemist's Intuition. Sci Rep 2020; 10:9218. [PMID: 32513975 PMCID: PMC7280257 DOI: 10.1038/s41598-020-65984-0] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/19/2020] [Accepted: 05/13/2020] [Indexed: 11/09/2022] Open
Abstract
Contemporary efforts for empirically-unbiased modeling of protein-ligand interactions entail a painful tradeoff - as reliable information on both noncovalent binding factors and the dynamic behavior of a protein-ligand complex is often beyond practical limits. We demonstrate that information drawn exclusively from static molecular structures can be used for reproducing and predicting experimentally-measured binding affinities for protein-ligand complexes. In particular, inhibition constants (Ki) were calculated for seven different competitive inhibitors of Torpedo californica acetylcholinesterase using a multiple-linear-regression-based model. The latter, incorporating five independent variables - drawn from QM cluster, DLPNO-CCSD(T) calculations and LED analyses on the seven complexes, each containing active amino-acid residues found within interacting distance (3.5 Å) from the corresponding ligand - is shown to recover 99.9% of the sum of squares for measured Ki values, while having no statistically-significant residual errors. Despite being fitted to a small number of data points, leave-one-out cross-validation statistics suggest that it possesses surprising predictive value (Q2LOO=0.78, or 0.91 upon removal of a single outlier). This thus challenges ligand-invariant definitions of active sites, such as implied in the lock-key binding theory, as well as in alternatives highlighting shape-complementarity without taking electronic effects into account. Broader implications of the current work are discussed in dedicated appendices.
Collapse
Affiliation(s)
- Nitai Sylvetsky
- Department of Organic Chemistry, Weizmann Institute of Science, 7610001, Rehovot, Israel.
| |
Collapse
|
26
|
Lang J, Antalík A, Veis L, Brandejs J, Brabec J, Legeza Ö, Pittner J. Near-Linear Scaling in DMRG-Based Tailored Coupled Clusters: An Implementation of DLPNO-TCCSD and DLPNO-TCCSD(T). J Chem Theory Comput 2020; 16:3028-3040. [PMID: 32275424 DOI: 10.1021/acs.jctc.0c00065] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
We present a new implementation of density matrix renormalization group based tailored coupled clusters method (TCCSD), which employs the domain-based local pair natural orbital approach (DLPNO). Compared to the previous local pair natural orbital (LPNO) version of the method, the new implementation is more accurate, offers more favorable scaling, and provides more consistent behavior across the variety of systems. On top of the singles and doubles, we include the perturbative triples correction (T), which is able to retrieve even more dynamic correlation. The methods were tested on three systems: tetramethyleneethane, oxo-Mn(Salen), and iron(II)-porphyrin model. The first two were revisited to assess the performance with respect to LPNO-TCCSD. For oxo-Mn(Salen), we retrieved between 99.8 and 99.9% of the total canonical correlation energy which is an improvement of 0.2% over the LPNO version in less than 63% of the total LPNO runtime. Similar results were obtained for iron(II)-porphyrin. When the perturbative triples correction was employed, irrespective of the active space size or system, the obtained energy differences between two spin states were within the chemical accuracy of 1 kcal/mol using the default DLPNO settings.
Collapse
Affiliation(s)
- Jakub Lang
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic.,Faculty of Sciences, Charles University, Albertov 6, 128 00 Prague 2, Czech Republic
| | - Andrej Antalík
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic.,Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, 12116 Prague 2, Czech Republic
| | - Libor Veis
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Jan Brandejs
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic.,Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, 12116 Prague 2, Czech Republic
| | - Jiří Brabec
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Örs Legeza
- Strongly Correlated Systems "Lendület" Research group, Wigner Research Centre for Physics, H-1525 Budapest, Hungary
| | - Jiří Pittner
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| |
Collapse
|
27
|
Magalhães RP, Fernandes HS, Sousa SF. Modelling Enzymatic Mechanisms with QM/MM Approaches: Current Status and Future Challenges. Isr J Chem 2020. [DOI: 10.1002/ijch.202000014] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- Rita P. Magalhães
- UCIBIO@REQUIMTE, BioSIMDepartamento de Biomedicina, Faculdade de Medicina da Universidade do Porto Alameda Professor Hernâni Monteiro 4200-319 Porto Portugal
| | - Henriques S. Fernandes
- UCIBIO@REQUIMTE, BioSIMDepartamento de Biomedicina, Faculdade de Medicina da Universidade do Porto Alameda Professor Hernâni Monteiro 4200-319 Porto Portugal
| | - Sérgio F. Sousa
- UCIBIO@REQUIMTE, BioSIMDepartamento de Biomedicina, Faculdade de Medicina da Universidade do Porto Alameda Professor Hernâni Monteiro 4200-319 Porto Portugal
| |
Collapse
|
28
|
Fabrizio A, Petraglia R, Corminboeuf C. Balancing Density Functional Theory Interaction Energies in Charged Dimers Precursors to Organic Semiconductors. J Chem Theory Comput 2020; 16:3530-3542. [DOI: 10.1021/acs.jctc.9b01193] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Affiliation(s)
- Alberto Fabrizio
- Laboratory for Computational Molecular Design, Institute of Chemical Sciences and Engineering, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
- National Centre for Computational Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Riccardo Petraglia
- Laboratory for Computational Molecular Design, Institute of Chemical Sciences and Engineering, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Clemence Corminboeuf
- Laboratory for Computational Molecular Design, Institute of Chemical Sciences and Engineering, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
- National Centre for Computational Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| |
Collapse
|
29
|
Chen JL, Sun T, Wang YB, Wang W. Toward a less costly but accurate calculation of the CCSD(T)/CBS noncovalent interaction energy. J Comput Chem 2020; 41:1252-1260. [PMID: 32045021 DOI: 10.1002/jcc.26171] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2019] [Revised: 01/12/2020] [Accepted: 02/01/2020] [Indexed: 01/14/2023]
Abstract
The popular method of calculating the noncovalent interaction energies at the coupled-cluster single-, double-, and perturbative triple-excitations [CCSD(T)] theory level in the complete basis set (CBS) limit was to add a CCSD(T) correction term to the CBS second-order Møller-Plesset perturbation theory (MP2). The CCSD(T) correction term is the difference between the CCSD(T) and MP2 interaction energies evaluated in a medium basis set. However, the CCSD(T) calculations with the medium basis sets are still very expensive for systems with more than 30 atoms. Comparatively, the domain-based local pair natural orbital coupled-cluster method [DLPNO-CCSD(T)] can be applied to large systems with over 1,000 atoms. Considering both the computational accuracy and efficiency, in this work, we propose a new scheme to calculate the CCSD(T)/CBS interaction energies. In this scheme, the MP2/CBS term keeps intact and the CCSD(T) correction term is replaced by a DLPNO-CCSD(T) correction term which is the difference between the DLPNO-CCSD(T) and DLPNO-MP2 interaction energies evaluated in a medium basis set. The interaction energies of the noncovalent systems in the S22, HSG, HBC6, NBC10, and S66 databases were recalculated employing this new scheme. The consistent and tight settings of the truncation parameters for DLPNO-CCSD(T) and DLPNO-MP2 in this noncanonical CCSD(T)/CBS calculations lead to the maximum absolute deviation and root-mean-square deviation from the canonical CCSD(T)/CBS interaction energies of less than or equal to 0.28 kcal/mol and 0.09 kcal/mol, respectively. The high accuracy and low cost of this new computational scheme make it an excellent candidate for the study of large noncovalent systems.
Collapse
Affiliation(s)
- Jiu-Li Chen
- Department of Chemistry, and Key Laboratory of Guizhou High Performance Computational Chemistry, Guizhou University, Guiyang, China
| | - Tao Sun
- Department of Chemistry, and Key Laboratory of Guizhou High Performance Computational Chemistry, Guizhou University, Guiyang, China
| | - Yi-Bo Wang
- Department of Chemistry, and Key Laboratory of Guizhou High Performance Computational Chemistry, Guizhou University, Guiyang, China
| | - Weizhou Wang
- College of Chemistry and Chemical Engineering, and Henan Key Laboratory of Function-Oriented Porous Materials, Luoyang Normal University, Luoyang, China
| |
Collapse
|
30
|
Li J, Rogachev AY. Homolytic Versus Heterolytic Bond Breaking in Functionalized [R-C 20 H 10 ] + Systems. J Comput Chem 2020; 41:88-96. [PMID: 31495954 DOI: 10.1002/jcc.26065] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/10/2019] [Revised: 08/12/2019] [Accepted: 08/14/2019] [Indexed: 01/05/2023]
Abstract
The comprehensive theoretical investigation of stability of functionalized corannulene cations [R-C20 H10 ]+ with respect to two alternative bond-breaking mechanisms, namely, homolytic or radical ([R-C20 H10 ]+ → R• + C20 H10 +• ) and heterolytic or cationic ([R-C20 H10 ]+ → R+ + C20 H10 ), was accomplished. The special focus was on the influence of the nature of R-group on the energetics of the bond cleavage. Detailed study of energetics of both mechanisms has revealed that the systems with small alkyl groups such as methyl tend to undergo bond breaking in accordance with homolytic mechanism. Subsequent elongation of the chain of the R-group resulted in shifting the paradigm, making heterolytic path more energetically favorable. Subsequent analysis of different components of the bonding between R-group and corannulene polyaromatic core helped to shed light on trends observed. In both mechanisms, the covalent contribution was found to be dominating, whereas ionic part contributes ~25-27%. Two leading components of ΔEorb , C20 H10 → R and R → C20 H10 , were identified with NOCV-EDA approach. While the homolytic pathway is best described as R → C20 H10 process, the heterolytic mechanism shows domination of the C20 H10 → R term. Surprisingly, the preparation energy (ΔEprep ) was identified as a key player in stability tendencies found. In other words, the relative stability of corresponding molecular fragments (here R-groups as the corannulene fragment remains the same for all systems) in their cationic or radical forms determine the preference given to a specific bond breaking path and, as consequence, the total stability of target functionalized cations. These conclusions were further confirmed by extending a set of R-groups to conjugated (allyl, phenyl), bulky (iPr, tBu), β-silyl (CH2 SiH3 , CH2 SiMe3 ), and benzyl (CH2 Ph) groups. © 2019 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Jingbai Li
- Department of Chemistry, Illinois Institute of Technology, Chicago, Illinois, 60616
| | - Andrey Yu Rogachev
- Department of Chemistry, Illinois Institute of Technology, Chicago, Illinois, 60616
| |
Collapse
|
31
|
Obst MF, Gevorgyan A, Bayer A, Hopmann KH. Mechanistic Insights into Copper-Catalyzed Carboxylations. Organometallics 2020. [DOI: 10.1021/acs.organomet.9b00710] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Affiliation(s)
- Marc F. Obst
- Hylleraas Center for Quantum Molecular Sciences, Department of Chemistry, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Ashot Gevorgyan
- Department of Chemistry, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Annette Bayer
- Department of Chemistry, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Kathrin H. Hopmann
- Hylleraas Center for Quantum Molecular Sciences, Department of Chemistry, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| |
Collapse
|
32
|
Abstract
The computational efficiency of local correlation methods is strongly dependent on the size of the domain of functions used to expand local correlating orbitals such as orbital specific or pair natural orbitals. Here, we define a principal domain of order m as the subset of m one-particle functions that provides the best support for a given n-electron wave function by maximizing the partial trace of the one-body reduced density matrix. Principal domains maximize the overlap between the wave function and its approximant for two-electron systems and are the domain selection equivalent of Löwdin's natural orbitals. We present an efficient linear scaling greedy algorithm for obtaining principal domains of projected atomic orbitals and demonstrate its utility in the context of the pair natural orbital local correlation theory. We numerically determine thresholds such that the projected atomic orbital domain error is an order of magnitude smaller than the pair natural orbital truncation error.
Collapse
Affiliation(s)
- David P Tew
- Max Planck Institute for Solid State Research , Heisenbergstr. 1 , 70569 Stuttgart , Germany
| |
Collapse
|
33
|
Kulsha AV, Sharapa DI. Superhalogen and Superacid. J Comput Chem 2019; 40:2293-2300. [PMID: 31254480 DOI: 10.1002/jcc.26007] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2018] [Revised: 05/25/2019] [Accepted: 05/30/2019] [Indexed: 11/10/2022]
Abstract
A superhalogen F@C20 (CN)20 and a corresponding Brønsted superacid were designed and investigated on DFT and DLPNO-CCSD(T) levels of theory. Calculated compounds have outstanding electron affinity and deprotonation energy, respectively. We consider superacid H[F@C20 (CN)20 ] to be able to protonate molecular nitrogen. The stability of these structures is discussed, while some of the previous predictions concerning neutral Brønsted superacids of record strength are doubted. © 2019 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Andrey V Kulsha
- Lyceum of Belarusian State University, 8 Ulijanauskaja Str., Minsk, 220030, Belarus
| | - Dmitry I Sharapa
- Chair of Theoretical Chemistry and Interdisciplinary Center for Molecular Materials, Friedrich-Alexander-Universitat Erlangen-Nürnberg, Egerlandstraße 3, 91058, Erlangen, Germany.,Institute of Catalysis Research and Technology (IKFT), Hermann-von-Helmholtz-Platz 1, Eggenstein-Leopoldshafen, D-76344, Germany
| |
Collapse
|
34
|
Koszinowski K, Auth T. Modulation of Gas-Phase Lithium Cation Basicities by Microsolvation. JOURNAL OF THE AMERICAN SOCIETY FOR MASS SPECTROMETRY 2019; 30:1857-1866. [PMID: 31502224 DOI: 10.1007/s13361-019-02312-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/28/2019] [Revised: 08/12/2019] [Accepted: 08/12/2019] [Indexed: 06/10/2023]
Abstract
In contrast to the extensive knowledge of lithium cation affinities and basicities, the thermochemistry of microsolvated lithium cations is much less explored. Here, we determine the relative stabilities of Li(A,B)n+ complexes, n = 2 and 3, by monitoring their gas-phase reactions with A and B substrate molecules, A/B = Me2O, Et2O, tetrahydrofuran, and MeCN, in a three-dimensional quadrupole-ion trap mass spectrometer. Kinetic analysis of the observed ligand displacement reactions affords equilibrium constants, which are then converted into Gibbs reaction energies. In addition, we use high-level quantum chemical calculations to predict the structures and sequential ligand dissociation energies of the homoleptic Li(A)n+ complexes, n = 1-3. As expected, the ligands dissociate more easily from complexes in higher coordination states. However, the very nature of the ligand also matters. Ligands with different steric demands can, thus, invert their relative Li+ affinities depending on the coordination state of the metal center. This finding shows that microsolvation of Li+ can result in specific effects, which are not recognized if the analysis takes into account only simple lithium cation affinities and basicities.
Collapse
Affiliation(s)
- Konrad Koszinowski
- Institut für Organische und Biomolekulare Chemie, Universität Göttingen, Tammannstr. 2, 37077, Göttingen, Germany.
| | - Thomas Auth
- Institut für Organische und Biomolekulare Chemie, Universität Göttingen, Tammannstr. 2, 37077, Göttingen, Germany
| |
Collapse
|
35
|
Goncalves TJ, Plessow PN, Studt F. On the Accuracy of Density Functional Theory in Zeolite Catalysis. ChemCatChem 2019. [DOI: 10.1002/cctc.201900791] [Citation(s) in RCA: 33] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
- Tiago J. Goncalves
- Institute of Catalysis Research and TechnologyKarlsruhe Institute of Technology Hermann-von-Helmholtz Platz 1 76344 Eggenstein-Leopoldshafen Germany
| | - Philipp N. Plessow
- Institute of Catalysis Research and TechnologyKarlsruhe Institute of Technology Hermann-von-Helmholtz Platz 1 76344 Eggenstein-Leopoldshafen Germany
| | - Felix Studt
- Institute of Catalysis Research and TechnologyKarlsruhe Institute of Technology Hermann-von-Helmholtz Platz 1 76344 Eggenstein-Leopoldshafen Germany
- Institute for Chemical Technology and Polymer ChemistryKarlsruhe Institute of Technology Engesserstrasse 18 76131 Karlsruhe Germany
| |
Collapse
|
36
|
Kahremany S, Kubas A, Tochtrop GP, Palczewski K. Catalytic synthesis of 9-cis-retinoids: mechanistic insights. Dalton Trans 2019; 48:10581-10595. [PMID: 31218312 PMCID: PMC7004310 DOI: 10.1039/c9dt02189b] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
The regioselective Z-isomerization of thermodynamically stable all-trans retinoids remains challenging, and ultimately limits the availability of much needed therapeutics for the treatment of human diseases. We present here a novel, straightforward approach for the catalytic Z-isomerization of retinoids using conventional heat treatment or microwave irradiation. A screen of 20 transition metal-based catalysts identified an optimal approach for the regioselective production of Z-retinoids. The most effective catalytic system was comprised of a palladium complex with labile ligands. Several mechanistic studies, including isotopic H/D exchange and state-of-the-art quantum chemical calculations using coupled cluster methods indicate that the isomerization is initiated by catalyst dimerization followed by the formation of a cyclic, six-membered chloropalladate catalyst-substrate adduct, which eventually opens to produce the desired Z-isomer. The synthetic development described here, combined with thorough mechanistic analysis of the underlying chemistry, highlights the use of readily available transition metal-based catalysts in straightforward formats for gram-scale drug synthesis.
Collapse
Affiliation(s)
- Shirin Kahremany
- Gavin Herbert Eye Institute, Department of Ophthalmology, University of California, Irvine, CA 92697, USA. and Department of Pharmacology, Case Western Reserve University, Cleveland, OH 44106, USA
| | - Adam Kubas
- Institute of Physical Chemistry, Polish Academy of Sciences, Kasprzaka 44/52, 01-224 Warsaw, Poland.
| | - Gregory P Tochtrop
- Department of Chemistry, Case Western Reserve University, Cleveland, OH 44106, USA
| | - Krzysztof Palczewski
- Gavin Herbert Eye Institute, Department of Ophthalmology, University of California, Irvine, CA 92697, USA. and Department of Pharmacology, Case Western Reserve University, Cleveland, OH 44106, USA
| |
Collapse
|
37
|
Maity B, Minenkov Y, Cavallo L. Evaluation of experimental alkali metal ion–ligand noncovalent bond strengths with DLPNO-CCSD(T) method. J Chem Phys 2019; 151:014301. [DOI: 10.1063/1.5099580] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Affiliation(s)
- Bholanath Maity
- King Abdullah University of Science and Technology (KAUST), Physical Science and Engineering Division (PSE), KAUST Catalysis Center (KCC), 23955-6900 Thuwal, Saudi Arabia
| | - Yury Minenkov
- Moscow Institute of Physics and Technology, Institutskiy Pereulok 9, Dolgoprudny, Moscow Region 141700, Russian Federation
| | - Luigi Cavallo
- King Abdullah University of Science and Technology (KAUST), Physical Science and Engineering Division (PSE), KAUST Catalysis Center (KCC), 23955-6900 Thuwal, Saudi Arabia
| |
Collapse
|
38
|
Disale ST, Brahmmananda Rao CVS, Gopakumar G, Jayaram RV. Experimental and theoretical studies on actinide extraction: dibutyl phenyl phosphonate versus tri-n-butyl phosphate. J COORD CHEM 2019. [DOI: 10.1080/00958972.2019.1614175] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
- S. T. Disale
- Department of Chemistry, Kankavli College, Sindhudurg, India
| | - C. V. S. Brahmmananda Rao
- Homi Bhabha National Institute, Indira Gandhi Centre for Atomic Research, Kalpakkam, India
- Mass Spectrometry Studies Section, Fuel Chemistry Division, Materials Chemistry & Metal Fuel Cycle Group, Indira Gandhi Centre for Atomic Research, Kalpakkam, India
| | - Gopinadhanpillai Gopakumar
- Mass Spectrometry Studies Section, Fuel Chemistry Division, Materials Chemistry & Metal Fuel Cycle Group, Indira Gandhi Centre for Atomic Research, Kalpakkam, India
| | - R. V. Jayaram
- Department of Chemistry, Institute of Chemical Technology, Mumbai, India
| |
Collapse
|
39
|
Xu Y, Xi S, Wang F, Li X. Theoretical Study on Reactions of Alkylperoxy Radicals. J Phys Chem A 2019; 123:3949-3958. [DOI: 10.1021/acs.jpca.9b01496] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Yanjiao Xu
- Institute of Atomic and Molecular Physics, Key Laboratory of High Energy Density Physics and Technology, Ministry of Education, Sichuan University, Chengdu 610064, P. R. China
| | - Shuanghui Xi
- College of Chemical Engineering, Sichuan University, Chengdu 610064, P. R. China
| | - Fan Wang
- Institute of Atomic and Molecular Physics, Key Laboratory of High Energy Density Physics and Technology, Ministry of Education, Sichuan University, Chengdu 610064, P. R. China
| | - Xiangyuan Li
- College of Chemical Engineering, Sichuan University, Chengdu 610064, P. R. China
| |
Collapse
|
40
|
Lang J, Brabec J, Saitow M, Pittner J, Neese F, Demel O. Perturbative triples correction to domain-based local pair natural orbital variants of Mukherjee's state specific coupled cluster method. Phys Chem Chem Phys 2019; 21:5022-5038. [PMID: 30762044 DOI: 10.1039/c8cp03577f] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
In this article we report an implementation of the perturbative triples correction to Mukherjee's state-specific multireference coupled cluster method based on the domain-based pair natural orbital approach (DLPNO-MkCC). We tested the performance of DLPNO-MkCCSD(T) in calculations involving tetramethyleneethane and isomers of naphthynes. These tests show that more than 97% of triples energy was recovered with respect to the canonical MkCCSD(T) method, which together with the DLPNO-MkCCSD part accounts for about 99.70-99.85% of the total correlation energy. The applicability of the method was demonstrated on calculations of singlet-triplet gaps for several large systems: triangulene, dynemicin A, and a beryllium complex.
Collapse
Affiliation(s)
- Jakub Lang
- J. Heyrovský Institute of Physical Chemistry, v.v.i., Academy of Sciences of the Czech Republic, Dolejškova 3, 18223 Prague 8, Czech Republic.
| | | | | | | | | | | |
Collapse
|
41
|
Vogiatzis KD, Polynski MV, Kirkland JK, Townsend J, Hashemi A, Liu C, Pidko EA. Computational Approach to Molecular Catalysis by 3d Transition Metals: Challenges and Opportunities. Chem Rev 2019; 119:2453-2523. [PMID: 30376310 PMCID: PMC6396130 DOI: 10.1021/acs.chemrev.8b00361] [Citation(s) in RCA: 214] [Impact Index Per Article: 42.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2018] [Indexed: 12/28/2022]
Abstract
Computational chemistry provides a versatile toolbox for studying mechanistic details of catalytic reactions and holds promise to deliver practical strategies to enable the rational in silico catalyst design. The versatile reactivity and nontrivial electronic structure effects, common for systems based on 3d transition metals, introduce additional complexity that may represent a particular challenge to the standard computational strategies. In this review, we discuss the challenges and capabilities of modern electronic structure methods for studying the reaction mechanisms promoted by 3d transition metal molecular catalysts. Particular focus will be placed on the ways of addressing the multiconfigurational problem in electronic structure calculations and the role of expert bias in the practical utilization of the available methods. The development of density functionals designed to address transition metals is also discussed. Special emphasis is placed on the methods that account for solvation effects and the multicomponent nature of practical catalytic systems. This is followed by an overview of recent computational studies addressing the mechanistic complexity of catalytic processes by molecular catalysts based on 3d metals. Cases that involve noninnocent ligands, multicomponent reaction systems, metal-ligand and metal-metal cooperativity, as well as modeling complex catalytic systems such as metal-organic frameworks are presented. Conventionally, computational studies on catalytic mechanisms are heavily dependent on the chemical intuition and expert input of the researcher. Recent developments in advanced automated methods for reaction path analysis hold promise for eliminating such human-bias from computational catalysis studies. A brief overview of these approaches is presented in the final section of the review. The paper is closed with general concluding remarks.
Collapse
Affiliation(s)
| | | | - Justin K. Kirkland
- Department
of Chemistry, University of Tennessee, Knoxville, Tennessee 37996, United States
| | - Jacob Townsend
- Department
of Chemistry, University of Tennessee, Knoxville, Tennessee 37996, United States
| | - Ali Hashemi
- Inorganic
Systems Engineering group, Department of Chemical Engineering, Delft University of Technology, Van der Maasweg 9, 2629 HZ Delft, The Netherlands
| | - Chong Liu
- Inorganic
Systems Engineering group, Department of Chemical Engineering, Delft University of Technology, Van der Maasweg 9, 2629 HZ Delft, The Netherlands
| | - Evgeny A. Pidko
- TheoMAT
group, ITMO University, Lomonosova 9, St. Petersburg 191002, Russia
- Inorganic
Systems Engineering group, Department of Chemical Engineering, Delft University of Technology, Van der Maasweg 9, 2629 HZ Delft, The Netherlands
| |
Collapse
|
42
|
Pollice R, Chen P. Origin of the Immiscibility of Alkanes and Perfluoroalkanes. J Am Chem Soc 2019; 141:3489-3506. [DOI: 10.1021/jacs.8b10745] [Citation(s) in RCA: 28] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
|
43
|
Ma P, Chen H. Ligand-Dependent Multi-State Reactivity in Cobalt(III)-Catalyzed C–H Activations. ACS Catal 2019. [DOI: 10.1021/acscatal.8b04532] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- Pengchen Ma
- Beijing National Laboratory for Molecular Sciences (BNLMS), CAS Key Laboratory of Photochemistry, CAS Research/Education Center for Excellence in Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, P. R. China
- University of Chinese Academy of Sciences, Beijing 100049, P. R. China
| | - Hui Chen
- Beijing National Laboratory for Molecular Sciences (BNLMS), CAS Key Laboratory of Photochemistry, CAS Research/Education Center for Excellence in Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, P. R. China
| |
Collapse
|
44
|
Cornaton Y, Djukic JP. A noncovalent interaction insight onto the concerted metallation deprotonation mechanism. Phys Chem Chem Phys 2019; 21:20486-20498. [DOI: 10.1039/c9cp03650d] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The CMD/AMLA mechanisms of cyclopalladation and the parent fictitious cyclonickelation of N,N-dimethylbenzylamine have been investigated by joint DFT-D and DLPNO-CCSD(T) methods assisted by QTAIM.
Collapse
Affiliation(s)
- Yann Cornaton
- Laboratoire de Mathématiques et de Physique
- F-66860 Perpignan
- France
- Institut de Chimie de Strasbourg
- UMR 7177
| | | |
Collapse
|
45
|
Sharapa DI, Genaev A, Cavallo L, Minenkov Y. A Robust and Cost‐Efficient Scheme for Accurate Conformational Energies of Organic Molecules. Chemphyschem 2018; 20:92-102. [DOI: 10.1002/cphc.201801063] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/14/2018] [Indexed: 11/09/2022]
Affiliation(s)
- Dmitry I. Sharapa
- Institute of Catalysis Research and TechnologyKarlsruhe Institute of Technology (KIT) Hermann-von-Helmholtz Platz 1 Eggenstein-Leopoldshafen D-76344 Germany
| | - Alexander Genaev
- Vorozhtsov Novosibirsk Institute of Organic Chemistry Academician Lavrent'ev Ave., 9 Novosibirsk 630090 Russian Federation
| | - Luigi Cavallo
- KAUST Catalysis Center (KCC)King Abdullah University of Science and Technology Thuwal- 23955-6900 Saudi Arabia
| | - Yury Minenkov
- Moscow Institute of Physics and Technology Institutskiy Pereulok 9, Dolgoprudny Moscow Region 141700 Russia
| |
Collapse
|
46
|
Kumar S, Maji S, Gopakumar G, Joseph M, Sundararajan K, Sankaran K. Luminescent versus non-luminescent uranyl–picolinate complexes. J Radioanal Nucl Chem 2018. [DOI: 10.1007/s10967-018-6305-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
|
47
|
|
48
|
Li J, Tadakamalla D, Rogachev AY. Modulating stability of functionalized fullerene cations [R-C 60
] +
with the nature of R-group. J Comput Chem 2018; 39:2385-2396. [DOI: 10.1002/jcc.25579] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/19/2018] [Revised: 07/25/2018] [Accepted: 07/31/2018] [Indexed: 01/26/2023]
Affiliation(s)
- Jingbai Li
- Department of Chemistry; Illinois Institute of Technology; Chicago Illinois, 60616
| | - Divya Tadakamalla
- Department of Chemistry; Illinois Institute of Technology; Chicago Illinois, 60616
| | - Andrey Yu. Rogachev
- Department of Chemistry; Illinois Institute of Technology; Chicago Illinois, 60616
| |
Collapse
|
49
|
Gatineau D, Lesage D, Clavier H, Dossmann H, Chan CH, Milet A, Memboeuf A, Cole RB, Gimbert Y. Bond dissociation energies of carbonyl gold complexes: a new descriptor of ligand effects in gold(i) complexes? Dalton Trans 2018; 47:15497-15505. [PMID: 30338332 DOI: 10.1039/c8dt03721c] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Ligand electronic effects in gold(i) chemistry have been evaluated by means of the experimental determination of M-CO bond dissociation energies for 16 [L-Au-CO]+ complexes, bearing L ligands widely used in gold catalysis. Energy-resolved analyses have been made using tandem mass spectrometry with collision-induced dissociation. Coupled with DFT calculations, this approach enables the quantification of ligand effects based on the LAu-CO bond strength. A further energy decomposition analysis gives access to detailed insights into this bond's characteristics. Whereas small differences are observed between phosphine- and phosphite-containing gold complexes, carbene ligands are shown to stabilize the gold-carbonyl bond much more efficiently.
Collapse
Affiliation(s)
- David Gatineau
- Univ. Grenoble Alpes and CNRS, DCM (UMR 5250) BP 53, 38041 Grenoble Cedex9, France.
| | | | | | | | | | | | | | | | | |
Collapse
|
50
|
Bahri-Laleh N, Hanifpour A, Mirmohammadi SA, Poater A, Nekoomanesh-Haghighi M, Talarico G, Cavallo L. Computational modeling of heterogeneous Ziegler-Natta catalysts for olefins polymerization. Prog Polym Sci 2018. [DOI: 10.1016/j.progpolymsci.2018.06.005] [Citation(s) in RCA: 87] [Impact Index Per Article: 14.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|