1
|
Salem J, Tóbiás R, Császár AG, Mogren Al-Mogren M, Jaidane NE, Hochlaf M. Temperature-Dependent Line-Broadening Effects in CO 2 Caused by Ar. Chemphyschem 2024; 25:e202300467. [PMID: 37916391 DOI: 10.1002/cphc.202300467] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2023] [Revised: 10/28/2023] [Accepted: 10/30/2023] [Indexed: 11/03/2023]
Abstract
This computational study of line-broadening effects is based on an accurate, analytical representation of the intermonomer potential energy surface (PES) of the CO2 ⋅ Ar van der Waals (vdW) complex. The PES is employed to compute collisional broadening coefficients for rovibrational lines of CO2 perturbed by Ar. The semiclassical computations are performed using the modified Robert-Bonamy approach, including real and imaginary terms, and the exact trajectory model. The lines investigated are in the 10001←00011, 01101←00001, 00011←00001, and 00031←00001 vibrational bands and the computations are repeated at multiple temperatures. The computed results are in good agreement with the available experimental values, validating both the intermonomer PES developed and the methodology used. For lines in the 01101←00001 band of CO2 , temperature-dependent Ar-broadening coefficients are reported for the first time. The parameters presented should prove useful, among other applications, for the accurate experimental determination of CO2 and Ar abundances in planetary atmospheres.
Collapse
Affiliation(s)
- Jamel Salem
- Université de Gafsa, Laboratoire Technologie, Energie, Matériaux Innovants, Faculté des Sciences, Sidi Ahmed Zarroug -, 2112, Gafsa, Tunisia
| | - Roland Tóbiás
- Laboratory of Molecular Structure and Dynamics, Institute of Chemistry, ELTE Eötvös Loránd University, Pázmány Péter sétány 1/A, H-1117, Budapest, Hungary
| | - Attila G Császár
- Laboratory of Molecular Structure and Dynamics, Institute of Chemistry, ELTE Eötvös Loránd University, Pázmány Péter sétány 1/A, H-1117, Budapest, Hungary
- HUN-REN-ELTE Complex Chemical Systems Research Group, H-1518, Budapest 112, P.O. Box 32, Hungary
| | - Muneerah Mogren Al-Mogren
- Department of Chemistry, College of Sciences, King Saud University, P.O. Box 2455, Riyadh, 11451, Saudi Arabia
| | - Nejm-Edine Jaidane
- Laboratoire de Spectroscopie Atomique, Moléculaire et Applications LSAMA, Université de Tunis Al Manar, Tunis, Tunisia
| | - Majdi Hochlaf
- Université Gustave Eiffel, COSYS/IMSE, 5 Bd Descartes, 77454, Champs sur Marne, France
| |
Collapse
|
2
|
Tóbiás R, Simkó I, Császár AG. Unusual Dynamics and Vibrational Fingerprints of van der Waals Dimers Formed by Linear Molecules and Rare-Gas Atoms. J Chem Theory Comput 2023. [PMID: 38032107 DOI: 10.1021/acs.jctc.3c00914] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/01/2023]
Abstract
Detailed structural, dynamical, and vibrational analyses have been performed for systems composed of linear triatomic molecules solvated by a single rare-gas atom, He, Ne, or Ar. Among the chromophores of these van der Waals (vdW) dimers, there are four neutral molecules (CO2, CS2, N2O, and OCS) and six molecular cations (HHe2+, HNe2+, HAr2+, HHeNe+, HHeAr+, and HNeAr+), both of apolar and polar nature. Following the exploration of bonding preferences, high-level four-dimensional (4D) potential energy surfaces (PESs) have been developed for 24 vdW dimers, keeping the two intramonomer bond lengths fixed. For these 24 complexes, over 1500 bound vibrational states have been obtained via quasi-variational nuclear-motion computations, employing exact kinetic-energy operators together with the accurate 4D PESs and their 2D/3D cuts. The reduced-dimensional (2D to 4D) dimer models have been compared with full-dimensional (6D) ones in the cases of the neutral CO2·Ar and charged HHe2+·He dimers, corroborating the high accuracy of the 2D to 4D vibrational energies. The reduced-dimensional models suggest that (a) while the equilibrium structures are T-shaped and planar, the effective ground-state structures are nonplanar, (b) certain bound states belong to collinear molecular structures, even when they are not minima, (c) the vdW vibrations are heavily mixed and many states have amplitudes corresponding to both the T-shaped and collinear structures, (d) there are a few dimers, for which even some of the vdW fundamentals lie above the first dissociation limit, and (e) the vdW vibrations are almost fully decoupled from the intramonomer bending motion.
Collapse
Affiliation(s)
- Roland Tóbiás
- HUN-REN-ELTE Complex Chemical Systems Research Group, P.O. Box 32, H-1518 Budapest 112, Hungary
| | - Irén Simkó
- Hevesy György PhD School of Chemistry, ELTE Eötvös Loránd University, Pázmány Péter sétány 1/A, H-1117 Budapest, Hungary
- Institute of Chemistry, ELTE Eötvös Loránd University, Pázmány Péter sétány 1/A, H-1117 Budapest, Hungary
| | - Attila G Császár
- HUN-REN-ELTE Complex Chemical Systems Research Group, P.O. Box 32, H-1518 Budapest 112, Hungary
- Institute of Chemistry, ELTE Eötvös Loránd University, Pázmány Péter sétány 1/A, H-1117 Budapest, Hungary
| |
Collapse
|
3
|
Derbali E, Ajili Y, Mehnen B, Żuchowski PS, Kędziera D, Al-Mogren MM, Jaidane NE, Hochlaf M. Towards the generation of potential energy surfaces of weakly bound medium-sized molecular systems: the case of benzonitrile-He complex. Phys Chem Chem Phys 2023; 25:30198-30210. [PMID: 37807943 DOI: 10.1039/d3cp02720a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/10/2023]
Abstract
Currently, the explicitly correlated coupled cluster method is used routinely to generate the multi-dimensional potential energy surfaces (mD-PESs) of van der Waals complexes of small molecular systems relevant for atmospheric, astrophysical and industrial applications. Although very accurate, this method is computationally prohibitive for medium and large molecules containing clusters. For instance, the recent detections of complex organic molecules (COMs) in the interstellar medium, such as benzonitrile, revealed the need to establish an accurate enough electronic structure approach to map the mD-PESs of these species interacting with the surrounding gases. As a benchmark, we have treated the case of the polar molecule benzonitrile interacting with helium, where we use post-Hartree-Fock and symmetry-adapted perturbation theory (SAPT) techniques. Accordingly, we show that MP2 and distinguishable-cluster approximation (DCSD) cannot be used for this purpose, whereas accurate enough PESs may be obtained using the corresponding explicitly correlated versions (MP2-F12 or DCSD-F12) with a reduction in computational costs. Alternatively, computations revealed that SAPT(DFT) is as performant as CCSD(T)-F12/aug-cc-pVTZ, making it the method of choice for mapping the mD-PESs of COMs containing clusters. Therefore, we have used this approach to generate the 3D-PES of the benzonitrile-He complex along the intermonomer Jacobi coordinates. As an application, we have incorporated the analytic form of this PES into quantum dynamical computations to determine the cross sections of the rotational (de-)excitation of benzonitrile colliding with helium at a collision energy of 10 cm-1.
Collapse
Affiliation(s)
- Eya Derbali
- Laboratoire de Spectroscopie Atomique, Moléculaire et Applications LSAMA, Université de Tunis Al Manar, Tunis, Tunisia.
| | - Yosra Ajili
- Laboratoire de Spectroscopie Atomique, Moléculaire et Applications LSAMA, Université de Tunis Al Manar, Tunis, Tunisia.
| | - Bilel Mehnen
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University in Toruń, Grudziadz Street 5, 87-100 Toruń, Poland
| | - Piotr S Żuchowski
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University in Toruń, Grudziadz Street 5, 87-100 Toruń, Poland
| | - Dariusz Kędziera
- Faculty of Chemistry, Nicolaus Copernicus University in Toruń, ul. Gagarina 7, PL 87-100 Toruń, Poland.
| | - Muneerah Mogren Al-Mogren
- Department of Chemistry, College of Sciences, King Saud University, PO Box 2455, Riyadh 11451, Saudi Arabia
| | - Nejm-Edine Jaidane
- Laboratoire de Spectroscopie Atomique, Moléculaire et Applications LSAMA, Université de Tunis Al Manar, Tunis, Tunisia.
| | - Majdi Hochlaf
- Université Gustave Eiffel, COSYS/IMSE, 5 Bd Descartes 77454, Champs sur Marne, France.
| |
Collapse
|
4
|
Vindel-Zandbergen P, Kȩdziera D, Żółtowski M, Kłos J, Żuchowski P, Felker PM, Lique F, Bačić Z. H2O-HCN complex: A new potential energy surface and intermolecular rovibrational states from rigorous quantum calculations. J Chem Phys 2023; 159:174302. [PMID: 37909452 DOI: 10.1063/5.0173751] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2023] [Accepted: 10/13/2023] [Indexed: 11/03/2023] Open
Abstract
In this work the H2O-HCN complex is quantitatively characterized in two ways. First, we report a new rigid-monomer 5D intermolecular potential energy surface (PES) for this complex, calculated using the symmetry-adapted perturbation theory based on density functional theory method. The PES is based on 2833 ab initio points computed employing the aug-cc-pVQZ basis set, utilizing the autoPES code, which provides a site-site analytical fit with the long-range region given by perturbation theory. Next, we present the results of the quantum 5D calculations of the fully coupled intermolecular rovibrational states of the H2O-HCN complex for the total angular momentum J values of 0, 1, and 2, performed on the new PES. These calculations rely on the quantum bound-state methodology developed by us recently and applied to a variety of noncovalently bound binary molecular complexes. The vibrationally averaged ground-state geometry of H2O-HCN determined from the quantum 5D calculations agrees very well with that from the microwave spectroscopic measurements. In addition, the computed ground-state rotational transition frequencies, as well as the B and C rotational constants calculated for the ground state of the complex, are in excellent agreement with the experimental values. The assignment of the calculated intermolecular vibrational states of the H2O-HCN complex is surprisingly challenging. It turns out that only the excitations of the intermolecular stretch mode can be assigned with confidence. The coupling among the angular degrees of freedom (DOFs) of the complex is unusually strong, and as a result most of the excited intermolecular states are unassigned. On the other hand, the coupling of the radial, intermolecular stretch mode and the angular DOFs is weak, allowing straightforward assignment of the excitation of the former.
Collapse
Affiliation(s)
| | - Dariusz Kȩdziera
- Faculty of Chemistry, Nicolaus Copernicus University in Toruń, ul. Gagarina 7, 87-100 Toruń, Poland
| | - Michał Żółtowski
- University of Rennes, CNRS, IPR (Institut de Physique de Rennes) - UMR 6251, F-35000 Rennes, France
- LOMC - UMR 6294, CNRS-Université du Havre, 25 rue Philippe Lebon, BP1123, 76 063 Le Havre cedex, France
| | - Jacek Kłos
- Joint Quantum Institute, University of Maryland, College Park, Maryland 20742, USA
| | - Piotr Żuchowski
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University, ul. Grudziądzka 5, 87-100 Toruń, Poland
| | - Peter M Felker
- Department of Chemistry and Biochemistry, University of California, Los Angeles, California 90095-1569, USA
| | - François Lique
- University of Rennes, CNRS, IPR (Institut de Physique de Rennes) - UMR 6251, F-35000 Rennes, France
| | - Zlatko Bačić
- Department of Chemistry, New York University, New York, New York 10003, USA
- Simons Center for Computational Physical Chemistry at New York University, New York, New York 10003, USA
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshan Road North, Shanghai 200062, China
| |
Collapse
|
5
|
Ochieng SA, Patkowski K. Accurate three-body noncovalent interactions: the insights from energy decomposition. Phys Chem Chem Phys 2023; 25:28621-28637. [PMID: 37874287 DOI: 10.1039/d3cp03938b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2023]
Abstract
An impressive collection of accurate two-body interaction energies for small complexes has been assembled into benchmark databases and used to improve the performance of multiple density functional, semiempirical, and machine learning methods. Similar benchmark data on nonadditive three-body energies in molecular trimers are comparatively scarce, and the existing ones are practically limited to homotrimers. In this work, we present a benchmark dataset of 20 equilibrium noncovalent interaction energies for a small but diverse selection of 10 heteromolecular trimers. The new 3BHET dataset presents complexes that combine different interactions including π-π, anion-π, cation-π, and various motifs of hydrogen and halogen bonding in each trimer. A detailed symmetry-adapted perturbation theory (SAPT)-based energy decomposition of the two- and three-body interaction energies shows that 3BHET consists of electrostatics- and dispersion-dominated complexes. The nonadditive three-body contribution is dominated by induction, but its influence on the overall bonding type in the complex (as exemplified by its position on the ternary diagram) is quite small. We also tested the extended SAPT (XSAPT) approach which is capable of including some nonadditive interactions in clusters of any size. The resulting three-body dispersion term (obtained from the many-body dispersion formalism) is mostly in good agreement with the supermolecular CCSD(T)-MP2 values and the nonadditive induction term is similar to the three-body SAPT(DFT) data, but the overall three-body XSAPT energies are not very accurate as they are missing the first-order exchange terms.
Collapse
Affiliation(s)
- Sharon A Ochieng
- Department of Chemistry and Biochemistry, Auburn University, Auburn, Alabama 36849, USA.
| | - Konrad Patkowski
- Department of Chemistry and Biochemistry, Auburn University, Auburn, Alabama 36849, USA.
| |
Collapse
|
6
|
Hapka M, Krzemińska A, Modrzejewski M, Przybytek M, Pernal K. Efficient Calculation of the Dispersion Energy for Multireference Systems with Cholesky Decomposition: Application to Excited-State Interactions. J Phys Chem Lett 2023; 14:6895-6903. [PMID: 37494637 PMCID: PMC10405273 DOI: 10.1021/acs.jpclett.3c01568] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2023] [Accepted: 07/17/2023] [Indexed: 07/28/2023]
Abstract
Accurate and efficient prediction of dispersion interactions in excited-state complexes poses a challenge due to the complex nature of electron correlation effects that need to be simultaneously considered. We propose an algorithm for computing the dispersion energy in nondegenerate ground- or excited-state complexes with arbitrary spin. The algorithm scales with the fifth power of the system size due to employing Cholesky decomposition of Coulomb integrals and a recently developed recursive formula for density response functions of the monomers. As a numerical illustration, we apply the new algorithm in the framework of multiconfigurational symmetry adapted perturbation theory, SAPT(MC), to study interactions in dimers with localized excitons. The SAPT(MC) analysis reveals that the dispersion energy may be the main force stabilizing excited-state dimers.
Collapse
Affiliation(s)
- Michał Hapka
- Faculty
of Chemistry, University of Warsaw, ul. L. Pasteura 1, 02-093 Warsaw, Poland
| | - Agnieszka Krzemińska
- Institute
of Physics, Lodz University of Technology, ul. Wolczanska 217/221, 93-005 Lodz, Poland
| | - Marcin Modrzejewski
- Faculty
of Chemistry, University of Warsaw, ul. L. Pasteura 1, 02-093 Warsaw, Poland
| | - Michał Przybytek
- Faculty
of Chemistry, University of Warsaw, ul. L. Pasteura 1, 02-093 Warsaw, Poland
| | - Katarzyna Pernal
- Institute
of Physics, Lodz University of Technology, ul. Wolczanska 217/221, 93-005 Lodz, Poland
| |
Collapse
|
7
|
Tyagi R, Zen A, Voora VK. Quantifying the Impact of Halogenation on Intermolecular Interactions and Binding Modes of Aromatic Molecules. J Phys Chem A 2023. [PMID: 37406194 DOI: 10.1021/acs.jpca.3c02291] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/07/2023]
Abstract
Halogenation of aromatic molecules is frequently used to modulate intermolecular interactions with ramifications for optoelectronic and mechanical properties. In this work, we accurately quantify and understand the nature of intermolecular interactions in perhalogenated benzene (PHB) clusters. Using benchmark binding energies from the fixed-node diffusion Monte Carlo (FN-DMC) method, we show that generalized Kohn-Sham semicanonical projected random phase approximation (GKS-spRPA) plus approximate exchange kernel (AKX) provides reliable interaction energies with mean absolute error (MAE) of 0.23 kcal/mol. Using the GKS-spRPA+AXK method, we quantify the interaction energies of several binding modes of PHB clusters ((C6X6)n; X = F, Cl, Br, I; n = 2, 3). For a given binding mode, the interaction energies increase 3-4 times from X = F to X = I; the X-X binding modes have energies in the range of 2-4 kcal/mol, while the π-π binding mode has interaction energies in the range of 4-12 kcal/mol. SAPT-DFT-based energy decomposition analysis is then used to show that the equilibrium geometries are dictated primarily by the dispersion and exchange interactions. Finally, we test the accuracy of several dispersion-corrected density functional approximations and show that only the r2SCAN-D4 method has a low MAE and correct long-range behavior, which makes it suitable for large-scale simulations and for developing structure-function relationships of halogenated aromatic systems.
Collapse
Affiliation(s)
- Ritaj Tyagi
- Department of Chemical Sciences, Tata Institute of Fundamental Research, Homi Bhabha Road, Colaba, Mumbai 400005, India
| | - Andrea Zen
- Dipartimento di Fisica Ettore Pancini, Università di Napoli Federico II, Monte S. Angelo, I-80126 Napoli, Italy
| | - Vamsee K Voora
- Department of Chemical Sciences, Tata Institute of Fundamental Research, Homi Bhabha Road, Colaba, Mumbai 400005, India
| |
Collapse
|
8
|
Heindel JP, Herman KM, Xantheas SS. Many-Body Effects in Aqueous Systems: Synergies Between Interaction Analysis Techniques and Force Field Development. Annu Rev Phys Chem 2023; 74:337-360. [PMID: 37093659 DOI: 10.1146/annurev-physchem-062422-023532] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/25/2023]
Abstract
Interaction analysis techniques, including the many-body expansion (MBE), symmetry-adapted perturbation theory, and energy decomposition analysis, allow for an intuitive understanding of complex molecular interactions. We review these methods by first providing a historical context for the study of many-body interactions and discussing how nonadditivities emerge from Hamiltonians containing strictly pairwise-additive interactions. We then elaborate on the synergy between these interaction analysis techniques and the development of advanced force fields aimed at accurately reproducing the Born-Oppenheimer potential energy surface. In particular, we focus on ab initio-based force fields that aim to explicitly reproduce many-body terms and are fitted to high-level electronic structure results. These force fields generally incorporate many-body effects through (a) parameterization of distributed multipoles, (b) explicit fitting of the MBE, (c) inclusion of many-atom features in a neural network, and (d) coarse-graining of many-body terms into an effective two-body term. We also discuss the emerging use of the MBE to improve the accuracy and speed of ab initio molecular dynamics.
Collapse
Affiliation(s)
- Joseph P Heindel
- Department of Chemistry, University of Washington, Seattle, Washington, USA
| | - Kristina M Herman
- Department of Chemistry, University of Washington, Seattle, Washington, USA
| | - Sotiris S Xantheas
- Department of Chemistry, University of Washington, Seattle, Washington, USA
- Advanced Computing, Mathematics and Data Division, Pacific Northwest National Laboratory, Richland, Washington, USA; ,
| |
Collapse
|
9
|
Morera-Boado C, Bernal-Uruchurtu MI. Interaction energy of Cl 2 and Br 2 with H 2 O: Exchange, dispersion and density the crucial ingredients. J Comput Chem 2023; 44:1073-1087. [PMID: 36578228 DOI: 10.1002/jcc.27066] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2022] [Revised: 11/10/2022] [Accepted: 12/18/2022] [Indexed: 12/30/2022]
Abstract
Modern Density Functional Theory models are now suitable for many molecular and condensed phase studies. The study of noncovalent interactions, a well-known drawback, is no longer an insurmountable obstacle through design and empirical corrections. However, using empirical corrections as in the DFT-D methods might not be an all-in-one solution. This work uses a simple system, X2 -H2 O with X = Cl or Br, with two different interactions, halogen-bonded (XB) and hydrogen-halogen (HX), to investigate the capability of current density functional approximations (DFA) in predicting interaction energies with eight different exchange-correlation functionals. SAPT(DFT) provides, for all the studied cases, better predictions than the widely used supermolecular approach. In addition, the components of the interaction energy suggest where some of the shortcomings originate in each DFA. The analysis of the functionals used confirms that PBE0 and ω-B97X-D have a physically correct behavior. Using SAPT(DFT) and PBE0, and ω-B97X-D, we obtained the interaction energy of Cl2 and Br2 inside different clathrate cages and satisfactorily compared with wavefunction results; hence, the lower and upper limits of this value are defined: Cl2 @512 , -5.3 ± 0.3 kcal/mol; Cl2 @512 62 , -5.5 ± 0.1 kcal/mol; Br2 @512 62 , -7.6 ± 1.0 kcal/mol; Br2 @512 63 , -10.6 ± 1.0 kcal/mol; Br2 @512 64 , -10.9 ± 0.8 kcal/mol.
Collapse
Affiliation(s)
- Cercis Morera-Boado
- CONACYT - Centro de Investigaciones Químicas, Universidad Autónoma del Estado de Morelos, Cuernavaca, Morelos, Mexico.,Centro de Investigaciones Químicas, IICBA, Universidad Autónoma del Estado de Morelos, Cuernavaca, Mexico
| | | |
Collapse
|
10
|
Oliveira BGD. Why much of Chemistry may be indisputably non-bonded? Semina: Ciênc Ex Tech 2023. [DOI: 10.5433/1679-0375.2022v43n2p211] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023] Open
Abstract
In this compendium, the wide scope of all intermolecular interactions ever known has been revisited, in particular giving emphasis the capability of much of the elements of the periodic table to form non-covalent contacts. Either hydrogen bonds, dihydrogen bonds, halogen bonds, pnictogen bonds, chalcogen bonds, triel bonds, tetrel bonds, regium bonds, spodium bonds or even the aerogen bond interactions may be cited. Obviously that experimental techniques have been used in some works, but it was through the theoretical methods that these interactions were validate, wherein the QTAIM integrations and SAPT energy partitions have been useful in this regard. Therefore, the great goal concerns to elucidate the interaction strength and if the intermolecular system shall be total, partial or non-covalently bonded, wherein this last one encompasses the most majority of the intermolecular interactions what leading to affirm that chemistry is debatably non-bonded.
Collapse
|
11
|
Luu D, Patkowski K. Overcoming Artificial Multipoles in Intramolecular Symmetry-Adapted Perturbation Theory. J Phys Chem A 2023; 127:356-377. [PMID: 36563050 DOI: 10.1021/acs.jpca.2c06465] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Intramolecular symmetry-adapted perturbation theory (ISAPT) is a method to compute and decompose the noncovalent interaction energy between two molecular fragments A and B covalently connected via a linker C. However, the existing ISAPT algorithm displays several issues for many fragmentation patterns (that is, specific assignments of atoms to the A/B/C subsystems), including an artificially repulsive electrostatic energy (even when the fragments are hydrogen-bonded) and very large and mutually cancelling induction and exchange-induction terms. We attribute those issues to the presence of artificial dipole moments at the interfragment boundary, as the atoms of A and B directly connected to C are missing electrons on one of their hybrid orbitals. Therefore, we propose several new partitioning algorithms which reassign one electron, on a singly occupied link hybrid orbital, from C to each of A/B. Once the contributions from these link orbitals are added to fragment density matrices, the computation of ISAPT electrostatic, induction, and dispersion energies proceeds exactly as normal, and the exchange energy expressions need only minor modifications. Among the link partitioning algorithms introduced, the so-called ISAPT(SIAO1) approach (in which the link orbital is obtained by a projection onto the intrinsic atomic orbitals (IAOs) of a given fragment followed by orthogonalization to this fragment's occupied space) leads to reasonable values of all ISAPT corrections for all fragmentation patterns, and exhibits a fast and systematic basis set convergence. This improvement is made possible by a significant reduction in magnitude (even though not a complete elimination) of the unphysical dipole moments at the interfragment boundaries. We demonstrate the utility of the improved ISAPT partitioning by examining intramolecular interactions in several pentanediol isomers, examples of linear and branched alkanes, and the open and closed conformations of a family of N-arylimide molecular torsion balances.
Collapse
Affiliation(s)
- Du Luu
- Department of Chemistry and Biochemistry, Auburn University, Auburn, Alabama 36849, United States
| | - Konrad Patkowski
- Department of Chemistry and Biochemistry, Auburn University, Auburn, Alabama 36849, United States
| |
Collapse
|
12
|
Masoumifeshani E, Chojecki M, Rutkowska-Zbik D, Korona T. Association Complexes of Calix[6]arenes with Amino Acids Explained by Energy-Partitioning Methods. Molecules 2022; 27. [PMID: 36432040 DOI: 10.3390/molecules27227938] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 09/23/2022] [Revised: 11/09/2022] [Accepted: 11/09/2022] [Indexed: 11/18/2022]
Abstract
Intermolecular complexes with calixarenes are intriguing because of multiple possibilities of noncovalent binding for both polar and nonpolar molecules, including docking in the calixarene cavity. In this contribution calix[6]arenes interacting with amino acids are studied with an additional aim to show that tools such as symmetry-adapted perturbation theory (SAPT), functional-group SAPT (F-SAPT), and systematic molecular fragmentation (SMF) methods may provide explanations for different numbers of noncovalent bonds and of their varying strength for various calixarene conformers and guest molecules. The partitioning of the interaction energy provides an easy way to identify hydrogen bonds, including those with unconventional hydrogen acceptors, as well as other noncovalent bonds, and to find repulsive destabilizing interactions between functional groups. Various other features can be explained by energy partitioning, such as the red shift of an IR stretching frequency for some hydroxy groups, which arises from their attraction to the phenyl ring of calixarene. Pairs of hydrogen bonds and other noncovalent bonds of similar magnitude found by F-SAPT explain an increase in the stability of both inclusion and outer complexes.
Collapse
|
13
|
Sadia H, Ullah S, Ullah F, Jadoon T. DFT study about capturing of toxic sulfur gases over cyclic tetrapyrrole. COMPUT THEOR CHEM 2022. [DOI: 10.1016/j.comptc.2022.113966] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/02/2022]
|
14
|
Szalewicz K, Jeziorski B. Physical mechanisms of intermolecular interactions from symmetry-adapted perturbation theory. J Mol Model 2022; 28:273. [PMID: 36006512 DOI: 10.1007/s00894-022-05190-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2022] [Accepted: 05/12/2022] [Indexed: 10/15/2022]
Abstract
Symmetry-adapted perturbation theory (SAPT) is a method for computational studies of noncovalent interactions between molecules. This method will be discussed here from the perspective of establishing the paradigm for understanding mechanisms of intermolecular interactions. SAPT interaction energies are obtained as sums of several contributions. Each contribution possesses a clear physical interpretation as it results from some specific physical process. It also exhibits a specific dependence on the intermolecular separation R. The four major contributions are the electrostatic, induction, dispersion, and exchange energies, each due to a different mechanism, valid at any R. In addition, at large R, SAPT interaction energies are seamlessly connected with the corresponding terms in the asymptotic multipole expansion of interaction energy in inverse powers of R. Since such expansion explicitly depends on monomers' multipole moments and polarizabilities, this connection provides additional insights by rigorously relating interaction energies to monomers' properties.
Collapse
Affiliation(s)
- Krzysztof Szalewicz
- Department of Physics and Astronomy, University of Delaware, Newark, DE, 19716, USA.
| | - Bogumił Jeziorski
- Faculty of Chemistry, University of Warsaw, Pasteura 1, 02093, Warsaw, Poland
| |
Collapse
|
15
|
Abstract
About half a century after its little-known beginnings, the quantum topological approach called QTAIM has grown into a widespread, but still not mainstream, methodology of interpretational quantum chemistry. Although often confused in textbooks with yet another population analysis, be it perhaps an elegant but somewhat esoteric one, QTAIM has been enriched with about a dozen other research areas sharing its main mathematical language, such as Interacting Quantum Atoms (IQA) or Electron Localisation Function (ELF), to form an overarching approach called Quantum Chemical Topology (QCT). Instead of reviewing the latter’s role in understanding non-covalent interactions, we propose a number of ideas emerging from the full consequences of the space-filling nature of topological atoms, and discuss how they (will) impact on interatomic interactions, including non-covalent ones. The architecture of a force field called FFLUX, which is based on these ideas, is outlined. A new method called Relative Energy Gradient (REG) is put forward, which is able, by computation, to detect which fragments of a given molecular assembly govern the energetic behaviour of this whole assembly. This method can offer insight into the typical balance of competing atomic energies both in covalent and non-covalent case studies. A brief discussion on so-called bond critical points is given, highlighting concerns about their meaning, mainly in the arena of non-covalent interactions.
Collapse
|
16
|
Brinck T, Borrfors AN. The Importance of Electrostatics and Polarization for Noncovalent Interactions: Ionic Hydrogen Bonds vs Ionic Halogen Bonds. J Mol Model 2022; 28:275. [PMID: 36006525 DOI: 10.1007/s00894-022-05189-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2021] [Accepted: 01/19/2022] [Indexed: 11/05/2022]
Abstract
A series of 26 hydrogen-bonded complexes between Br− and halogen, oxygen and sulfur hydrogen-bond (HB) donors is investigated at the M06-2X/6–311 + G(2df,2p) level of theory. Analysis using a model in which Br− is replaced by a point charge shows that the interaction energy (\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\Delta E}_{Int}$$\end{document}ΔEInt) of the complexes is accurately reproduced by the scaled interaction energy with the point charge (\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\Delta E}_{Int}^{PC}$$\end{document}ΔEIntPC).This is demonstrated by \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\Delta E}_{Int}=0.86{\Delta E}_{Int}^{PC}$$\end{document}ΔEInt=0.86ΔEIntPC with a correlation coefficient, R2 =0.999. The only outlier is (Br-H-Br)−, which generally is classified as a strong charge-transfer complex with covalent character rather than a HB complex. \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\Delta E}_{Int}^{PC}$$\end{document}ΔEIntPC can be divided rigorously into an electrostatic contribution (\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\Delta E}_{ES}^{PC}$$\end{document}ΔEESPC) and a polarization contribution (\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\Delta E}_{Pol}^{PC}$$\end{document}ΔEPolPC).Within the set of HB complexes investigated, the former varies between -7.2 and -32.7 kcal mol−1, whereas the latter varies between -1.6 and -11.5 kcal mol−1. Compared to our previous study of halogen-bonded (XB) complexes between Br− and C–Br XB donors, the electrostatic contribution is generally stronger and the polarization contribution is generally weaker in the HB complexes. However, for both types of bonding, the variation in interaction strength can be reproduced accurately without invoking a charge-transfer term. For the Br−···HF complex, the importance of charge penetration on the variation of the interaction energy with intermolecular distance is investigated. It is shown that the repulsive character of \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\Delta E}_{Int}$$\end{document}ΔEInt at short distances in this complex to a large extent can be attributed to charge penetration.
Collapse
|
17
|
Abstract
An inexpensive and reliable method for molecular crystal structure predictions (CSPs) has been developed. The new CSP protocol starts from a two-dimensional graph of crystal’s monomer(s) and utilizes no experimental information. Using results of quantum mechanical calculations for molecular dimers, an accurate two-body, rigid-monomer ab initio-based force field (aiFF) for the crystal is developed. Since CSPs with aiFFs are essentially as expensive as with empirical FFs, tens of thousands of plausible polymorphs generated by the crystal packing procedures can be optimized. Here we show the robustness of this protocol which found the experimental crystal within the 20 most stable predicted polymorphs for each of the 15 investigated molecules. The ranking was further refined by performing periodic density-functional theory (DFT) plus dispersion correction (pDFT+D) calculations for these 20 top-ranked polymorphs, resulting in the experimental crystal ranked as number one for all the systems studied (and the second polymorph, if known, ranked in the top few). Alternatively, the polymorphs generated can be used to improve aiFFs, which also leads to rank one predictions. The proposed CSP protocol should result in aiFFs replacing empirical FFs in CSP research. Developing theoretical frameworks to predict new polymorphs is highly desirable. Here the authors present an ab initio based force-field approach for crystal structure prediction offering a dramatic computational speed-up over fully ab initio schemes.
Collapse
|
18
|
Barcza B, Szirmai ÁB, Szántó KJ, Tajti A, Szalay PG. Comparison of approximate intermolecular potentials for ab initio fragment calculations on medium sized N-heterocycles. J Comput Chem 2022; 43:1079-1093. [PMID: 35478353 PMCID: PMC9321956 DOI: 10.1002/jcc.26866] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2021] [Revised: 03/26/2022] [Accepted: 03/29/2022] [Indexed: 01/15/2023]
Abstract
The ground state intermolecular potential of bimolecular complexes of N‐heterocycles is analyzed for the impact of individual terms in the interaction energy as provided by various, conceptually different theories. Novel combinations with several formulations of the electrostatic, Pauli repulsion, and dispersion contributions are tested at both short‐ and long‐distance sides of the potential energy surface, for various alignments of the pyrrole dimer as well as the cytosine–uracil complex. The integration of a DFT/CCSD density embedding scheme, with dispersion terms from the effective fragment potential (EFP) method is found to provide good agreement with a reference CCSD(T) potential overall; simultaneously, a quantum mechanics/molecular mechanics approach using CHELPG atomic point charges for the electrostatic interaction, augmented by EFP dispersion and Pauli repulsion, comes also close to the reference result. Both schemes have the advantage of not relying on predefined force fields; rather, the interaction parameters can be determined for the system under study, thus being excellent candidates for ab initio modeling.
Collapse
Affiliation(s)
- Bónis Barcza
- Institute of Chemistry, Laboratory of Theoretical Chemistry, ELTE Eötvös Loránd University, Budapest, Hungary
| | - Ádám B Szirmai
- Institute of Chemistry, Laboratory of Theoretical Chemistry, ELTE Eötvös Loránd University, Budapest, Hungary
| | - Katalin J Szántó
- Institute of Chemistry, Laboratory of Theoretical Chemistry, ELTE Eötvös Loránd University, Budapest, Hungary
| | - Attila Tajti
- Institute of Chemistry, Laboratory of Theoretical Chemistry, ELTE Eötvös Loránd University, Budapest, Hungary
| | - Péter G Szalay
- Institute of Chemistry, Laboratory of Theoretical Chemistry, ELTE Eötvös Loránd University, Budapest, Hungary
| |
Collapse
|
19
|
Fedorov DG. Polarization energies in the fragment molecular orbital method. J Comput Chem 2022; 43:1094-1103. [PMID: 35446441 DOI: 10.1002/jcc.26869] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/04/2022] [Revised: 03/25/2022] [Accepted: 04/05/2022] [Indexed: 12/23/2022]
Abstract
Using isolated and polarized states of fragments, a method for computing the polarization energies in density functional theory (DFT) and density-functional tight-binding (DFTB) is developed in the framework of the fragment molecular orbital method. For DFTB, the method is extended into the use of periodic boundary conditions (PBC), for which a new component, a periodic self-polarization energy, is derived. The couplings of the polarization to other components in the pair interaction energy analysis (PIEDA) are derived for DFT and DFTB, and compared to Hartree-Fock and second-order Møller-Plesset perturbation theory (MP2). The effect of the self-consistent (DFT) and perturbative (MP2) treatment of the electron correlation on the polarization is discussed. The difference in the polarization in the bulk (PBC) and micro (cluster) solvation is elucidated.
Collapse
Affiliation(s)
- Dmitri G Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Tsukuba, Japan
| |
Collapse
|
20
|
Bamdad F, Farrokhpour H, Ashrafizaadeh M, Najafi B. Decomposition of the interaction energy of several flavonoids with Escherichia coli DNA Gyr using the SAPT (DFT) method: The relation between the interaction energy components, ligand structure, and biological activity. Biochim Biophys Acta Gen Subj 2022. [DOI: 10.1016/j.bbagen.2022.130111] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2021] [Revised: 01/19/2022] [Accepted: 02/07/2022] [Indexed: 12/28/2022]
|
21
|
Abstract
This paper accounts for a general procedure of bonding analysis that is, expectedly, adequate to describe any type of interaction involving the noble-gas (Ng) atoms. Building on our recently proposed classification of the Ng-X bonds (X = binding partner) [New J. Chem. 44, 15536 (2020)], these contacts are first distinguished into three types, namely, A, B, or C, based on the topology of the electron energy density H(r) and on the shape of its plotted form. Bonds of type B or C are, then, further assigned as B-loose (Bl) or B-tight (Bt) and C-loose (Cl) or C-tight (Ct) depending on the sign that H(r) takes along the Ng-X bond path located from the topological analysis of ρ(r), particularly at around the bond critical point (BCP). Any bond of type A, Bl/Bt, or Cl/Ct is, finally, assayed in terms of contribution of covalency. This is accomplished by studying the maximum, minimum, and average value of H(r) over the volume enclosed by the low-density reduced density gradient (RDG) isosurface associated with the bond (typically, the RDG isosurface including the BCP) and the average ρ(r) over the same volume. The bond assignment is also corroborated by calculating the values of quantitative indices specifically defined for the various types of interactions (A, B, or C). The generality of our taken approach should encourage its wide application to the study of Ng compounds.
Collapse
Affiliation(s)
- Stefano Borocci
- Dipartimento per la Innovazione nei sistemi Biologici, Agroalimentari e Forestali (DIBAF), Università della Tuscia, L.go dell'Università, s.n.c., 01100 Viterbo, Italy
| | - Felice Grandinetti
- Dipartimento per la Innovazione nei sistemi Biologici, Agroalimentari e Forestali (DIBAF), Università della Tuscia, L.go dell'Università, s.n.c., 01100 Viterbo, Italy
| | - Nico Sanna
- Dipartimento per la Innovazione nei sistemi Biologici, Agroalimentari e Forestali (DIBAF), Università della Tuscia, L.go dell'Università, s.n.c., 01100 Viterbo, Italy
| |
Collapse
|
22
|
Tang Z, Song Y, Zhang S, Wang W, Xu Y, Wu D, Wu W, Su P. XEDA, a fast and multipurpose energy decomposition analysis program. J Comput Chem 2021; 42:2341-2351. [PMID: 34626430 DOI: 10.1002/jcc.26765] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2021] [Revised: 09/16/2021] [Accepted: 09/22/2021] [Indexed: 01/06/2023]
Abstract
A fast and multipurpose energy decomposition analysis (EDA) program, called XEDA, is introduced for quantitative analysis of intermolecular interactions. This program contains a series of variational EDA methods, including LMO-EDA, GKS-EDA and their extensions, to analyze non-covalent interactions and strong chemical bonds in various environments. XEDA is highly efficient with a similar computational scaling of single point energy calculations. Its efficiency and universality are validated by a series of test examples including van der Waals interactions, hydrogen bonds, radical-radical interactions and strong covalent bonds.
Collapse
Affiliation(s)
- Zhen Tang
- The State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, China
| | - Yanlin Song
- The State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, China
| | - Shu Zhang
- The State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, China
| | - Wei Wang
- The State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, China
| | - Yuan Xu
- The State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, China
| | - Di Wu
- The State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, China
| | - Wei Wu
- The State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, China
| | - Peifeng Su
- The State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, China
| |
Collapse
|
23
|
Abstract
![]()
We present a formulation
of the multiconfigurational (MC) wave
function symmetry-adapted perturbation theory (SAPT). The method is
applicable to noncovalent interactions between monomers which require
a multiconfigurational description, in particular when the interacting
system is strongly correlated or in an electronically excited state.
SAPT(MC) is based on one- and two-particle reduced density matrices
of the monomers and assumes the single-exchange approximation for
the exchange energy contributions. Second-order terms are expressed
through response properties from extended random phase approximation
(ERPA). The dispersion components of SAPT(MC) have been introduced
in our previous works [HapkaM.2019, 15, 1016−102730525591; HapkaM.2019, 15, 6712–672331670950]. SAPT(MC) is applied either with generalized valence
bond perfect pairing (GVB) or with complete active space self-consistent
field (CASSCF) treatment of the monomers. We discuss two model multireference
systems: the H2 ··· H2 dimer
in out-of-equilibrium geometries and interaction between the argon
atom and excited state of ethylene. Using the C2H4* ··· Ar complex as an example, we examine second-order
terms arising from negative transitions in the linear response function
of an excited monomer. We demonstrate that the negative-transition
terms must be accounted for to ensure qualitative prediction of induction
and dispersion energies and develop a procedure allowing for their
computation. Factors limiting the accuracy of SAPT(MC) are discussed
in comparison with other second-order SAPT schemes on a data set of
small single-reference dimers.
Collapse
Affiliation(s)
- Michał Hapka
- Institute of Physics, Lodz University of Technology, ul. Wolczanska 219, 90-924 Lodz, Poland.,Faculty of Chemistry, University of Warsaw, ul. L. Pasteura 1, 02-093 Warsaw, Poland
| | - Michał Przybytek
- Faculty of Chemistry, University of Warsaw, ul. L. Pasteura 1, 02-093 Warsaw, Poland
| | - Katarzyna Pernal
- Institute of Physics, Lodz University of Technology, ul. Wolczanska 219, 90-924 Lodz, Poland
| |
Collapse
|
24
|
Cukras J, Sadlej J. Towards Quantum-Chemical Modeling of the Activity of Anesthetic Compounds. Int J Mol Sci 2021; 22:9272. [PMID: 34502179 DOI: 10.3390/ijms22179272] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2021] [Revised: 08/16/2021] [Accepted: 08/17/2021] [Indexed: 12/16/2022] Open
Abstract
The modeling of the activity of anesthetics is a real challenge because of their unique electronic and structural characteristics. Microscopic approaches relevant to the typical features of these systems have been developed based on the advancements in the theory of intermolecular interactions. By stressing the quantum chemical point of view, here, we review the advances in the field highlighting differences and similarities among the chemicals within this group. The binding of the anesthetics to their partners has been analyzed by Symmetry-Adapted Perturbation Theory to provide insight into the nature of the interaction and the modeling of the adducts/complexes allows us to rationalize their anesthetic properties. A new approach in the frame of microtubule concept and the importance of lipid rafts and channels in membranes is also discussed.
Collapse
|
25
|
Masoumifeshani E, Korona T. Symmetrized systematic molecular fragmentation model and its application for molecular properties. COMPUT THEOR CHEM 2021; 1202:113303. [DOI: 10.1016/j.comptc.2021.113303] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
|
26
|
Naseem-Khan S, Gresh N, Misquitta AJ, Piquemal JP. Assessment of SAPT and Supermolecular EDA Approaches for the Development of Separable and Polarizable Force Fields. J Chem Theory Comput 2021; 17:2759-2774. [DOI: 10.1021/acs.jctc.0c01337] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Sehr Naseem-Khan
- Laboratoire de Chimie Théorique, Sorbonne Université, UMR 7616 CNRS, 75005 Paris, France
| | - Nohad Gresh
- Laboratoire de Chimie Théorique, Sorbonne Université, UMR 7616 CNRS, 75005 Paris, France
| | - Alston J. Misquitta
- School of Physics and Astronomy and the Thomas Young Centre for Theory and Simulation of Materials at Queen Mary University of London, London E1 4NS, U.K
| | - Jean-Philip Piquemal
- Laboratoire de Chimie Théorique, Sorbonne Université, UMR 7616 CNRS, 75005 Paris, France
- Institut Universitaire de France, 75005 Paris, France
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| |
Collapse
|
27
|
Abstract
The Born-Oppenheimer potential energy surface (PES) has come a long way since its introduction in the 1920s, both conceptually and in predictive power for practical applications. Nevertheless, nearly 100 years later-despite astonishing advances in computational power-the state-of-the-art first-principles prediction of observables related to spectroscopy and scattering dynamics is surprisingly limited. For example, the water dimer, (H2O)2, with only six nuclei and 20 electrons, still presents a formidable challenge for full-dimensional variational calculations of bound states and is considered out of reach for rigorous scattering calculations. The extremely poor scaling of the most rigorous quantum methods is fundamental; however, recent progress in development of approximate methodologies has opened the door to fairly routine high-quality predictions, unthinkable 20 years ago. In this review, in relation to the workflow of spectroscopy and/or scattering studies, we summarize progress and challenges in the component areas of electronic structure calculations, PES fitting, and quantum dynamical calculations.
Collapse
Affiliation(s)
- Ernesto Quintas-Sánchez
- Department of Chemistry, Missouri University of Science and Technology, Rolla, Missouri 65409, USA;
| | - Richard Dawes
- Department of Chemistry, Missouri University of Science and Technology, Rolla, Missouri 65409, USA;
| |
Collapse
|
28
|
Li W, Chen J, Xu Y, Lu T, Gou Q, Feng G. Unveiling the structural and energetic properties of thiazole-water complex by microwave spectroscopy and theoretical calculations. Spectrochim Acta A Mol Biomol Spectrosc 2020; 242:118720. [PMID: 32736219 DOI: 10.1016/j.saa.2020.118720] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/26/2020] [Revised: 07/07/2020] [Accepted: 07/09/2020] [Indexed: 06/11/2023]
Abstract
The structure and non-covalent bonding features of the complex of thiazole and water were studied by using supersonic jet Fourier transform microwave spectroscopy and theoretical calculations. One isomer was observed which corresponds to the global minimum of the complex predicted theoretically. The rotational spectra of 9 additional isotopologues, including 5 mono-substituted heavy atoms of thiazole (34S, 13C and 15N), and 4 water isotopic species (H218O, DOH, HOD and D2O), were also measured and analyzed. The experimental spectroscopic parameters were used to determine the structural parameters of the observed isomer. Theoretical analyses based on quantum theory of atoms in molecules and natural bond orbital revealed that the two moieties are linked by a N···H-O hydrogen bond. The total interaction energy of the complex is calculated to be -25.1 kJmol-1 with electrostatics being the major term according to energy decomposition analysis.
Collapse
Affiliation(s)
- Wenqin Li
- School of Chemistry and Chemical Engineering, Chongqing University, Daxuecheng South Rd. 55, 401331 Chongqing, China
| | - Junhua Chen
- School of Chemistry and Chemical Engineering, Chongqing University, Daxuecheng South Rd. 55, 401331 Chongqing, China
| | - Yugao Xu
- School of Chemistry and Chemical Engineering, Chongqing University, Daxuecheng South Rd. 55, 401331 Chongqing, China
| | - Tao Lu
- School of Chemistry and Chemical Engineering, Chongqing University, Daxuecheng South Rd. 55, 401331 Chongqing, China
| | - Qian Gou
- School of Chemistry and Chemical Engineering, Chongqing University, Daxuecheng South Rd. 55, 401331 Chongqing, China
| | - Gang Feng
- School of Chemistry and Chemical Engineering, Chongqing University, Daxuecheng South Rd. 55, 401331 Chongqing, China.
| |
Collapse
|
29
|
Affiliation(s)
- C. David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia
30332-0400, USA
| | - David E. Manolopoulos
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, Oxford University, South Parks Road, Oxford OX1
3QZ, United Kingdom
| | - Todd J. Martínez
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305,
USA
| | - Angelos Michaelides
- Thomas Young Centre, London Centre for Nanotechnology and Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
| |
Collapse
|