1
|
Gupta AK, Stulajter MM, Shaidu Y, Neaton JB, de Jong WA. Equivariant Neural Networks Utilizing Molecular Clusters for Accurate Molecular Crystal Lattice Energy Predictions. ACS OMEGA 2024; 9:40269-40282. [PMID: 39346862 PMCID: PMC11425815 DOI: 10.1021/acsomega.4c07434] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/13/2024] [Revised: 08/27/2024] [Accepted: 09/02/2024] [Indexed: 10/01/2024]
Abstract
Equivariant neural networks have emerged as prominent models in advancing the construction of interatomic potentials due to their remarkable data efficiency and generalization capabilities for out-of-distribution data. Here, we expand the utility of these networks to the prediction of crystal structures consisting of organic molecules. Traditional methods for computing crystal structure properties, such as plane-wave quantum chemical methods based on density functional theory (DFT), are prohibitively resource-intensive, often necessitating compromises in accuracy and the choice of exchange-correlation functional. We present an approach that leverages the efficiency, and transferability of equivariant neural networks, specifically Allegro, to predict molecular crystal structure energies at a reduced computational cost. Our neural network is trained on molecular clusters using a highly accurate Gaussian-type orbital (GTO)-based method as the target level of theory, eliminating the need for costly periodic DFT calculations, while providing access to all families of exchange-corelation functionals and post-Hartree-Fock methods. The trained model exhibits remarkable accuracy in predicting lattice energies, aligning closely with those computed by plane-wave based DFT methods, thus representing significant cost reductions. Furthermore, the Allegro network was seamlessly integrated with the USPEX framework, accelerating the discovery of low-energy crystal structures during crystal structure prediction.
Collapse
Affiliation(s)
- Ankur K Gupta
- Applied Mathematics and Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Miko M Stulajter
- Applied Mathematics and Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Yusuf Shaidu
- Department of Physics, University of California Berkeley, Berkeley, California 94720, United States
- Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Jeffrey B Neaton
- Department of Physics, University of California Berkeley, Berkeley, California 94720, United States
- Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
- Kavli Energy NanoSciences Institute at Berkeley, Berkeley, California 94720, United States
| | - Wibe A de Jong
- Applied Mathematics and Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
2
|
Carter-Fenk K, Head-Gordon M. Repartitioned Brillouin-Wigner perturbation theory with a size-consistent second-order correlation energy. J Chem Phys 2023; 158:234108. [PMID: 37338032 PMCID: PMC10284609 DOI: 10.1063/5.0150033] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2023] [Accepted: 06/05/2023] [Indexed: 06/21/2023] Open
Abstract
Second-order Møller-Plesset perturbation theory (MP2) often breaks down catastrophically in small-gap systems, leaving much to be desired in its performance for myriad chemical applications such as noncovalent interactions, thermochemistry, and dative bonding in transition metal complexes. This divergence problem has reignited interest in Brillouin-Wigner perturbation theory (BWPT), which is regular at all orders but lacks size consistency and extensivity, severely limiting its application to chemistry. In this work, we propose an alternative partitioning of the Hamiltonian that leads to a regular BWPT perturbation series that, through the second order, is size-extensive, size-consistent (provided its Hartree-Fock reference is also), and orbital invariant. Our second-order size-consistent Brillouin-Wigner (BW-s2) approach can describe the exact dissociation limit of H2 in a minimal basis set, regardless of the spin polarization of the reference orbitals. More broadly, we find that BW-s2 offers improvements relative to MP2 for covalent bond breaking, noncovalent interaction energies, and metal/organic reaction energies, although rivaling coupled-cluster with single and double substitutions for thermochemical properties.
Collapse
Affiliation(s)
- Kevin Carter-Fenk
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | | |
Collapse
|
3
|
Förster A. Assessment of the Second-Order Statically Screened Exchange Correction to the Random Phase Approximation for Correlation Energies. J Chem Theory Comput 2022; 18:5948-5965. [PMID: 36150190 PMCID: PMC9558381 DOI: 10.1021/acs.jctc.2c00366] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
With increasing interelectronic distance, the screening
of the
electron–electron interaction by the presence of other electrons
becomes the dominant source of electron correlation. This effect is
described by the random phase approximation (RPA) which is therefore
a promising method for the calculation of weak interactions. The success
of the RPA relies on the cancellation of errors, which can be traced
back to the violation of the crossing symmetry of the 4-point vertex,
leading to strongly overestimated total correlation energies. By the
addition of second-order screened exchange (SOSEX) to the correlation
energy, this issue is substantially reduced. In the adiabatic connection
(AC) SOSEX formalism, one of the two electron–electron interaction
lines in the second-order exchange term is dynamically screened (SOSEX(W, vc)). A
related SOSEX expression in which both electron–electron interaction
lines are statically screened (SOSEX(W(0), W(0))) is obtained from the G3W2 contribution to the electronic self-energy. In contrast to SOSEX(W, vc), the
evaluation of this correlation energy expression does not require
an expensive numerical frequency integration and is therefore advantageous
from a computational perspective. We compare the accuracy of the statically
screened variant to RPA and RPA+SOSEX(W, vc) for a wide range of chemical
reactions. While both methods fail for barrier heights, SOSEX(W(0), W(0)) agrees very well with SOSEX(W, vc) for
charged excitations and noncovalent interactions where they lead to
major improvements over RPA.
Collapse
Affiliation(s)
- Arno Förster
- Theoretical Chemistry, Vrije Universiteit, De Boelelaan 1083, NL-1081 HV, Amsterdam, The Netherlands
| |
Collapse
|
4
|
Shee J, Loipersberger M, Rettig A, Lee J, Head-Gordon M. Regularized Second-Order Møller-Plesset Theory: A More Accurate Alternative to Conventional MP2 for Noncovalent Interactions and Transition Metal Thermochemistry for the Same Computational Cost. J Phys Chem Lett 2021; 12:12084-12097. [PMID: 34910484 PMCID: PMC10037552 DOI: 10.1021/acs.jpclett.1c03468] [Citation(s) in RCA: 28] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/14/2023]
Abstract
Second-order Møller-Plesset theory (MP2) notoriously breaks down for π-driven dispersion interactions and dative bonds in transition metal complexes. Herein, we investigate three physically justified forms of single-parameter, energy-gap dependent regularization which can yield high and transferable accuracy for a variety of noncovalent interactions (including S22, S66, and L7 test sets) and (mostly closed shell) transition metal thermochemistry. Regularization serves to damp overestimated pairwise additive contributions, renormalizing first-order amplitudes such that the effects of higher-order correlations are incorporated. The optimal parameter values for the noncovalent and transition metal sets are 1.1, 0.7, and 0.4 for κ, σ, and σ2 regularizers, respectively. However, such regularization slightly degrades the accuracy of conventional MP2 for some small-molecule test sets, most of which have relatively large average frontier energy gaps. Our results suggest that appropriately regularized MP2 models may improve double hybrid density functionals, at no additional cost over conventional MP2.
Collapse
Affiliation(s)
- James Shee
- Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Matthias Loipersberger
- Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Adam Rettig
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Joonho Lee
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
5
|
Loipersberger M, Bertels LW, Lee J, Head-Gordon M. Exploring the Limits of Second- and Third-Order Møller-Plesset Perturbation Theories for Noncovalent Interactions: Revisiting MP2.5 and Assessing the Importance of Regularization and Reference Orbitals. J Chem Theory Comput 2021; 17:5582-5599. [PMID: 34382394 PMCID: PMC9948597 DOI: 10.1021/acs.jctc.1c00469] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
This work systematically assesses the influence of reference orbitals, regularization, and scaling on the performance of second- and third-order Møller-Plesset perturbation theory wave function methods for noncovalent interactions (NCIs). Testing on 19 data sets (A24, DS14, HB15, HSG, S22, X40, HW30, NC15, S66, AlkBind12, CO2Nitrogen16, HB49, Ionic43, TA13, XB18, Bauza30, CT20, XB51, and Orel26rad) covers a wide range of different NCIs including hydrogen bonding, dispersion, and halogen bonding. Inclusion of potential energy surfaces from different hydrogen bonds and dispersion-bound complexes gauges accuracy for nonequilibrium geometries. Fifteen methods are tested. In notation where nonstandard choices of orbitals are denoted as methods:orbitals, these are MP2, κ-MP2, SCS-MP2, OOMP2, κ-OOMP2, MP3, MP2.5, MP3:OOMP2, MP2.5:OOMP2, MP3:κ-OOMP2, MP2.5:κ-OOMP2, κ-MP3:κ-OOMP2, κ-MP2.5:κ-OOMP2, MP3:ωB97X-V, and MP2.5:ωB97X-V. Furthermore, we compare these methods to the ωB97M-V and B3LYP-D3 density functionals, as well as CCSD. We find that the κ-regularization (κ = 1.45 au was used throughout) improves the energetics in almost all data sets for both MP2 (in 17 out of 19 data sets) and OOMP2 (16 out of 19). The improvement is significant (e.g., the root-mean-square deviation (RMSD) for the S66 data set is 0.29 kcal/mol for κ-OOMP2 versus 0.67 kcal/mol for MP2) and for interactions between stable closed-shell molecules, not strongly dependent on the reference orbitals. Scaled MP3 (with a factor of 0.5) using κ-OOMP2 reference orbitals (MP2.5:κ-OOMP2) provides significantly more accurate results for NCIs across all data sets with noniterative O(N6) scaling (S66 data set RMSD: 0.10 kcal/mol). Across the entire data set of 356 points, the improvement over standard MP2.5 is approximately a factor of 2: RMSD for MP3:κ-OOMP2 is 0.25 vs 0.50 kcal/mol for MP2.5. The use of high-quality density functional reference orbitals (ωB97X-V) also significantly improves the results of MP2.5 for NCI over a Hartree-Fock orbital reference. All our assessments and conclusions are based on the use of the medium-sized aug-cc-pVTZ basis to yield results that are directly compared against complete basis set limit reference values.
Collapse
Affiliation(s)
| | - Luke W. Bertels
- Department of Chemistry, University of California, Berkeley, California 94720, USA,Present Address: Department of Chemistry, Virginia Tech, Blacksburg, VA 24061, USA
| | - Joonho Lee
- Department of Chemistry, University of California, Berkeley, California 94720, USA,Present Address: Department of Chemistry, Columbia University, NY
| | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, USA,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
6
|
Bertels LW, Lee J, Head-Gordon M. Third-Order Møller-Plesset Perturbation Theory Made Useful? Choice of Orbitals and Scaling Greatly Improves Accuracy for Thermochemistry, Kinetics, and Intermolecular Interactions. J Phys Chem Lett 2019; 10:4170-4176. [PMID: 31259560 DOI: 10.1021/acs.jpclett.9b01641] [Citation(s) in RCA: 29] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
We develop and test methods that include second- and third-order perturbation theory (MP3) using orbitals obtained from regularized orbital-optimized second-order perturbation theory, κ-OOMP2, denoted as MP3:κ-OOMP2. Testing MP3:κ-OOMP2 shows RMS errors that are 1.7-5 times smaller than those of MP3 across 7 data sets. To do still better, empirical training of the scaling factors for the second- and third-order correlation energies and the regularization parameter on one of those data sets led to an unregularized scaled (c2 = 1.0; c3 = 0.8) denoted as MP2.8:κ-OOMP2. MP2.8:κ-OOMP2 yields significant additional improvement over MP3:κ-OOMP2 in 4 of 6 test data sets on thermochemistry, kinetics, and noncovalent interactions. Remarkably, these two methods outperform coupled cluster with singles and doubles in 5 of the 7 data sets considered, at greatly reduced cost (no O(N6) iterations).
Collapse
Affiliation(s)
- Luke W Bertels
- Department of Chemistry , University of California , Berkeley , California 94720 , United States
| | - Joonho Lee
- Department of Chemistry , University of California , Berkeley , California 94720 , United States
| | - Martin Head-Gordon
- Department of Chemistry , University of California , Berkeley , California 94720 , United States
- Chemical Sciences Division , Lawrence Berkeley National Laboratory , Berkeley , California 94720 , United States
| |
Collapse
|
7
|
Cacelli I, Lipparini F, Greff da Silveira L, Jacobs M, Livotto PR, Prampolini G. Accurate interaction energies by spin component scaled Möller-Plesset second order perturbation theory calculations with optimized basis sets (SCS-MP2mod): Development and application to aromatic heterocycles. J Chem Phys 2019; 150:234113. [DOI: 10.1063/1.5094288] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Affiliation(s)
- Ivo Cacelli
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, I-56124 Pisa, Italy
- Istituto di Chimica dei Composti OrganoMetallici (ICCOM-CNR), Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, I-56124 Pisa, Italy
| | - Leandro Greff da Silveira
- Instituto de Química, Universidade Federal do Rio Grande do Sul, Avenida Bento Gonçalves 9500, CEP 91501-970 Porto Alegre, Brazil
| | - Matheus Jacobs
- Institut für Physik, Humboldt-Universität zu Berlin, Newtonstr. 15, 12489 Berlin, Germany
- IRIS Adelrshof, Humboldt-Universität zu Berlin, Zum Großen Windkanal 6, 12489 Berlin, Germany
| | - Paolo Roberto Livotto
- Instituto de Química, Universidade Federal do Rio Grande do Sul, Avenida Bento Gonçalves 9500, CEP 91501-970 Porto Alegre, Brazil
| | - Giacomo Prampolini
- Istituto di Chimica dei Composti OrganoMetallici (ICCOM-CNR), Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
| |
Collapse
|
8
|
Calabrese C, Li W, Prampolini G, Evangelisti L, Uriarte I, Cacelli I, Melandri S, Cocinero EJ. A General Treatment to Study Molecular Complexes Stabilized by Hydrogen‐, Halogen‐, and Carbon‐Bond Networks: Experiment and Theory of (CH
2
F
2
)
n
⋅⋅⋅(H
2
O)
m. Angew Chem Int Ed Engl 2019; 58:8437-8442. [DOI: 10.1002/anie.201902753] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2019] [Revised: 04/05/2019] [Indexed: 11/08/2022]
Affiliation(s)
- Camilla Calabrese
- Departament of Physical ChemistryUniversity of the Basque Country (UPV/EHU) Barrio Sarriena, S/N 48940 Leioa Spain
- Biofisika Institute, (CSIC, UPV/EHU) 48080 Bilbao Spain
| | - Weixing Li
- Dipartimento di Chimica “Giacomo Ciamician”Università degli Studi di Bologna via Selmi 2 I-40126 Bologna Italy
| | - Giacomo Prampolini
- Istituto di Chimica dei Composti OrganoMetallici (ICCOM-CNR)Area della Ricerca via G. Moruzzi 1 I-56124 Pisa Italy
| | - Luca Evangelisti
- Dipartimento di Chimica “Giacomo Ciamician”Università degli Studi di Bologna via Selmi 2 I-40126 Bologna Italy
- Department of ChemistryUniversity of Virginia McCormick Road VA 22903 Charlottesville USA
| | - Iciar Uriarte
- Departament of Physical ChemistryUniversity of the Basque Country (UPV/EHU) Barrio Sarriena, S/N 48940 Leioa Spain
- Biofisika Institute, (CSIC, UPV/EHU) 48080 Bilbao Spain
| | - Ivo Cacelli
- Istituto di Chimica dei Composti OrganoMetallici (ICCOM-CNR)Area della Ricerca via G. Moruzzi 1 I-56124 Pisa Italy
- Dipartimento di Chimica e Chimica IndustrialeUniversità di Pisa via Risorgimento 35 I-56126 Pisa Italy
| | - Sonia Melandri
- Dipartimento di Chimica “Giacomo Ciamician”Università degli Studi di Bologna via Selmi 2 I-40126 Bologna Italy
| | - Emilio J. Cocinero
- Departament of Physical ChemistryUniversity of the Basque Country (UPV/EHU) Barrio Sarriena, S/N 48940 Leioa Spain
- Biofisika Institute, (CSIC, UPV/EHU) 48080 Bilbao Spain
| |
Collapse
|
9
|
Calabrese C, Li W, Prampolini G, Evangelisti L, Uriarte I, Cacelli I, Melandri S, Cocinero EJ. A General Treatment to Study Molecular Complexes Stabilized by Hydrogen‐, Halogen‐, and Carbon‐Bond Networks: Experiment and Theory of (CH
2
F
2
)
n
⋅⋅⋅(H
2
O)
m
. Angew Chem Int Ed Engl 2019. [DOI: 10.1002/ange.201902753] [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)
- Camilla Calabrese
- Departament of Physical ChemistryUniversity of the Basque Country (UPV/EHU) Barrio Sarriena, S/N 48940 Leioa Spain
- Biofisika Institute, (CSIC, UPV/EHU) 48080 Bilbao Spain
| | - Weixing Li
- Dipartimento di Chimica “Giacomo Ciamician”Università degli Studi di Bologna via Selmi 2 I-40126 Bologna Italy
| | - Giacomo Prampolini
- Istituto di Chimica dei Composti OrganoMetallici (ICCOM-CNR)Area della Ricerca via G. Moruzzi 1 I-56124 Pisa Italy
| | - Luca Evangelisti
- Dipartimento di Chimica “Giacomo Ciamician”Università degli Studi di Bologna via Selmi 2 I-40126 Bologna Italy
- Department of ChemistryUniversity of Virginia McCormick Road VA 22903 Charlottesville USA
| | - Iciar Uriarte
- Departament of Physical ChemistryUniversity of the Basque Country (UPV/EHU) Barrio Sarriena, S/N 48940 Leioa Spain
- Biofisika Institute, (CSIC, UPV/EHU) 48080 Bilbao Spain
| | - Ivo Cacelli
- Istituto di Chimica dei Composti OrganoMetallici (ICCOM-CNR)Area della Ricerca via G. Moruzzi 1 I-56124 Pisa Italy
- Dipartimento di Chimica e Chimica IndustrialeUniversità di Pisa via Risorgimento 35 I-56126 Pisa Italy
| | - Sonia Melandri
- Dipartimento di Chimica “Giacomo Ciamician”Università degli Studi di Bologna via Selmi 2 I-40126 Bologna Italy
| | - Emilio J. Cocinero
- Departament of Physical ChemistryUniversity of the Basque Country (UPV/EHU) Barrio Sarriena, S/N 48940 Leioa Spain
- Biofisika Institute, (CSIC, UPV/EHU) 48080 Bilbao Spain
| |
Collapse
|
10
|
Lee J, Head-Gordon M. Regularized Orbital-Optimized Second-Order Møller–Plesset Perturbation Theory: A Reliable Fifth-Order-Scaling Electron Correlation Model with Orbital Energy Dependent Regularizers. J Chem Theory Comput 2018; 14:5203-5219. [DOI: 10.1021/acs.jctc.8b00731] [Citation(s) in RCA: 62] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Joonho Lee
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
11
|
Jacobs M, Greff Da Silveira L, Prampolini G, Livotto PR, Cacelli I. Interaction Energy Landscapes of Aromatic Heterocycles through a Reliable yet Affordable Computational Approach. J Chem Theory Comput 2018; 14:543-556. [DOI: 10.1021/acs.jctc.7b00602] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Affiliation(s)
- Matheus Jacobs
- Instituto
de Química, Universidade Federal do Rio Grande do Sul, Avenida Bento Gonçalves 9500, CEP 91501-970 Porto Alegre, Brazil
| | - Leandro Greff Da Silveira
- Instituto
de Química, Universidade Federal do Rio Grande do Sul, Avenida Bento Gonçalves 9500, CEP 91501-970 Porto Alegre, Brazil
- Departamento
de Ciências Exatas e da Terra, Universidade Regional Integrada do Alto Uruguay e da Missões (URI), Avenida Assis Brasil 709, CEP 98400-00 Frederico Westphalen, Brazil
| | - Giacomo Prampolini
- Istituto di Chimica
dei Composti OrganoMetallici (ICCOM-CNR), Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
| | - Paolo Roberto Livotto
- Instituto
de Química, Universidade Federal do Rio Grande do Sul, Avenida Bento Gonçalves 9500, CEP 91501-970 Porto Alegre, Brazil
| | - Ivo Cacelli
- Istituto di Chimica
dei Composti OrganoMetallici (ICCOM-CNR), Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
- Dipartimento
di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi
3, I-56124 Pisa, Italy
| |
Collapse
|
12
|
Barone V, Cacelli I, Ferretti A, Prampolini G. Noncovalent Interactions in the Catechol Dimer. Biomimetics (Basel) 2017; 2:E18. [PMID: 31105180 PMCID: PMC6352673 DOI: 10.3390/biomimetics2030018] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2017] [Revised: 09/04/2017] [Accepted: 09/05/2017] [Indexed: 12/02/2022] Open
Abstract
Noncovalent interactions play a significant role in a wide variety of biological processes and bio-inspired species. It is, therefore, important to have at hand suitable computational methods for their investigation. In this paper, we report on the contribution of dispersion and hydrogen bonds in both stacked and T-shaped catechol dimers, with the aim of delineating the respective role of these classes of interactions in determining the most stable structure. By using second-order Møller⁻Plesset (MP2) calculations with a small basis set, specifically optimized for these species, we have explored a number of significant sections of the interaction potential energy surface and found the most stable structures for the dimer, in good agreement with the highly accurate, but computationally more expensive coupled cluster single and double excitation and the perturbative triples (CCSD(T))/CBS) method.
Collapse
Affiliation(s)
- Vincenzo Barone
- Scuola Normale Superiore di Pisa, Piazza dei Cavalieri, I-56126 Pisa, Italy.
| | - Ivo Cacelli
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, I-56124 Pisa, Italy.
- Istituto di Chimica dei Composti OrganoMetallici (ICCOM-CNR), Area della Ricerca, Via G. Moruzzi 1, I-56124 Pisa, Italy.
| | - Alessandro Ferretti
- Istituto di Chimica dei Composti OrganoMetallici (ICCOM-CNR), Area della Ricerca, Via G. Moruzzi 1, I-56124 Pisa, Italy.
| | - Giacomo Prampolini
- Istituto di Chimica dei Composti OrganoMetallici (ICCOM-CNR), Area della Ricerca, Via G. Moruzzi 1, I-56124 Pisa, Italy.
| |
Collapse
|
13
|
Tan S, Barrera Acevedo S, Izgorodina EI. Generalized spin-ratio scaled MP2 method for accurate prediction of intermolecular interactions for neutral and ionic species. J Chem Phys 2017; 146:064108. [DOI: 10.1063/1.4975326] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Samuel Tan
- Monash Computational Chemistry Group, School of Chemistry, Monash University, 17 Rainforest Walk, Clayton VIC 3800, Australia
| | - Santiago Barrera Acevedo
- School of Mathematical Sciences, Monash University, 9 Rainforest Walk, Clayton VIC 3800, Australia
| | - Ekaterina I. Izgorodina
- Monash Computational Chemistry Group, School of Chemistry, Monash University, 17 Rainforest Walk, Clayton VIC 3800, Australia
| |
Collapse
|
14
|
Torubaev YV, Skabitskiy IV, Pavlova AV, Pasynskii AA. First structural evidence of a Se–Br–Br halogen-bonded molecular complex. NEW J CHEM 2017. [DOI: 10.1039/c6nj04096a] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Bromination of the bicyclic selenoether (2,6-dibromo-9-selenabicyclo[3.3.1]-nonane) affords the products ranging from the unique molecular complex C8H12Br2Se–Br2 to the ionic polybromide complex [C8H12Br2SeBr+]Br5.
Collapse
Affiliation(s)
- Yury V. Torubaev
- N.S. Kurnakov Institute of General and Inorganic Chemistry
- Russian Academy of Sciences 119991
- Moscow
- Russian Federation
| | - Ivan V. Skabitskiy
- N.S. Kurnakov Institute of General and Inorganic Chemistry
- Russian Academy of Sciences 119991
- Moscow
- Russian Federation
| | - Alina V. Pavlova
- N.S. Kurnakov Institute of General and Inorganic Chemistry
- Russian Academy of Sciences 119991
- Moscow
- Russian Federation
| | - Alexander A. Pasynskii
- N.S. Kurnakov Institute of General and Inorganic Chemistry
- Russian Academy of Sciences 119991
- Moscow
- Russian Federation
| |
Collapse
|
15
|
Grimme S, Hansen A, Brandenburg JG, Bannwarth C. Dispersion-Corrected Mean-Field Electronic Structure Methods. Chem Rev 2016; 116:5105-54. [DOI: 10.1021/acs.chemrev.5b00533] [Citation(s) in RCA: 799] [Impact Index Per Article: 99.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Affiliation(s)
- Stefan Grimme
- Mulliken Center for Theoretical
Chemistry, Universität Bonn, 53113 Bonn, Germany
| | - Andreas Hansen
- Mulliken Center for Theoretical
Chemistry, Universität Bonn, 53113 Bonn, Germany
| | | | - Christoph Bannwarth
- Mulliken Center for Theoretical
Chemistry, Universität Bonn, 53113 Bonn, Germany
| |
Collapse
|
16
|
Abstract
Interest in molecular crystals has grown thanks to their relevance to pharmaceuticals, organic semiconductor materials, foods, and many other applications. Electronic structure methods have become an increasingly important tool for modeling molecular crystals and polymorphism. This article reviews electronic structure techniques used to model molecular crystals, including periodic density functional theory, periodic second-order Møller-Plesset perturbation theory, fragment-based electronic structure methods, and diffusion Monte Carlo. It also discusses the use of these models for predicting a variety of crystal properties that are relevant to the study of polymorphism, including lattice energies, structures, crystal structure prediction, polymorphism, phase diagrams, vibrational spectroscopies, and nuclear magnetic resonance spectroscopy. Finally, tools for analyzing crystal structures and intermolecular interactions are briefly discussed.
Collapse
Affiliation(s)
- Gregory J O Beran
- Department of Chemistry, University of California , Riverside, California 92521, United States
| |
Collapse
|