1
|
Kim S, Conrad JA, Tow GM, Maginn EJ, Boatz JA, Gordon MS. Intermolecular interactions in clusters of ethylammonium nitrate and 1-amino-1,2,3-triazole. Phys Chem Chem Phys 2023; 25:30428-30457. [PMID: 37917371 DOI: 10.1039/d3cp02407e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/04/2023]
Abstract
The intermolecular interaction energies, including hydrogen bonds (H-bonds), of clusters of the ionic liquid ethylammonium nitrate (EAN) and 1-amino-1,2,3-triazole (1-AT) based deep eutectic propellants (DeEP) are examined. 1-AT is introduced as a neutral hydrogen bond donor (HBD) to EAN in order to form a eutectic mixture. The effective fragment potential (EFP) is used to examine the bonding interactions in the DeEP clusters. The resolution of the Identity (RI) approximated second order Møller-Plesset perturbation theory (RI-MP2) and coupled cluster theory (RI-CCSD(T)) are used to validate the EFP results. The EFP method predicts that there are significant polarization and charge transfer effects in the EAN:1-AT complexes, along with Coulombic, dispersion and exchange repulsion interactions. The EFP interaction energies are in good agreement with the RI-MP2 and RI-CCSD(T) results. The quasi-atomic orbital (QUAO) bonding and kinetic bond order (KBO) analyses are additionally used to develop a conceptual and semi-quantitative understanding of the H-bonding interactions as a function of the size of the system. The QUAO and KBO analyses suggest that the H-bonds in the examined clusters follow the characteristic hydrogen bonding three-center four electron interactions. The strongest H-bonding interactions between the (EAN)1:(1-AT)n and (EAN)2:(1-AT)n (n = 1-5) complexes are observed internally within EAN; that is, between the ethylammonium cation [EA]+ and the nitrate anion ([NO3]-). The weakest H-bonding interactions occur between [NO3]- and 1-AT. Consequently, the average strengths of the H-bonds within a given (EAN)x:(1-AT)n complex decrease as more 1-AT molecules are introduced into the EAN monomer and EAN dimer. The QUAO bonding analysis suggests that 1-AT in (EAN)x:(1-AT)n can act as both a HBD and a hydrogen bond acceptor simultaneously. It is observed that two 1-AT molecules can form H-bonds to each other. Although the KBOs that correspond to H-bonding interactions in [EA]+:1-AT, [NO3]-:1-AT and between two 1-AT molecules are weaker than the H-bonds in EAN, those weak H-bond networks with 1-AT could be important to form a stable DeEP.
Collapse
Affiliation(s)
- Shinae Kim
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, IA 50011, USA.
- Combustion Research Facility, Sandia National Laboratories, Livermore, California 94550, USA
| | - Justin A Conrad
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, IA 50011, USA.
| | - Garrett M Tow
- Department of Chemical & Biomolecular Engineering, University of Notre Dame, Notre Dame, Indiana 46556, USA
| | - Edward J Maginn
- Department of Chemical & Biomolecular Engineering, University of Notre Dame, Notre Dame, Indiana 46556, USA
| | - Jerry A Boatz
- Aerospace Systems Directorate, Air Force Research Laboratory, Edwards Air Force Base, California 93524, USA
| | - Mark S Gordon
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, IA 50011, USA.
| |
Collapse
|
2
|
Oliveira RSS, Oliveira MAS, Borges I. The effect of the electronic structure method and basis set on the accuracy of the electric multipoles computed with the distributed multipole analysis (DMA). J Mol Model 2023; 29:357. [PMID: 37917318 DOI: 10.1007/s00894-023-05758-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2023] [Accepted: 10/18/2023] [Indexed: 11/04/2023]
Abstract
CONTEXT An accurate description of the molecular charge density is crucial for investigating intra- and inter-molecular properties. Among the different ways of describing and analyzing it, the widely used distributed multipole analysis (DMA) is an accurate method for decomposing the molecular charge density into atom-centered electric multipoles (monopole, dipole, quadrupole, and so on) that have a direct chemical interpretation. In this work, DMA was employed to decompose the molecular charge density of six chemically distinct molecules, namely, (2R)-2-amino-3-[(S)-prop-2-enylsulfinyl] propanoic acid (AAP), 4-amine-2-nitro-1,3,5 triazole (ANTA), (RS)-Propan-2-yl methylphosphonofluoridate (SARIN), chloromethane (CLMET), and 2-aminoacetic acid (GLY) into monopole, dipole, and quadrupole values. A hypothetical variation of ANTA built by exchanging all the nitrogen atoms with phosphorus that we named 4-phosphine-2-phosphite-1,3,5-phosphorine (ANTAP) was also studied. These molecules have different chemical structures bearing distinct carbon skeletons, electronegative atoms, and electron-withdrawing/donating groups. We found that although DFT multipole values can depend considerably on the exchange-correlation functional for specific atomic sites, the associated root-mean-square errors (RMSEs) compared to benchmark MP4 mainly were about [Formula: see text] The most significant variations were for monopoles and dipoles of sites highly polarized by adjacent atoms, and to a lesser degree, for the quadrupoles. The double hybrid B2PLYP and the hybrid meta M06-2X functionals, as expected in the framework of Jacob's ladder, overall give the most accurate results among the DFT methods. The MP2 DMA multipole values have an RMSE in relation to the MP4 benchmark mainly in the range [Formula: see text], thus representing a lower computational cost to obtain results with similar good accuracy without the ambiguity of choosing a DFT functional. The deviations of the HF multipoles from the benchmark in most cases were less than 20%, in agreement with the well-known fact that non-correlated charge densities have a slight dependence on the electronic correlation. We also confirmed that DMA values have a small dependence on the size of the basis set: deviations did not exceed 5% in most cases. However, the dependence of the DMA values on the size of the basis set increases with the rank of the electric multipole. To compute accurate values of DMA multipoles of an atom bonded to very electronegative atoms, especially dipoles and quadrupoles, a large basis set including diffuse functions is necessary. Despite that, for a given polarized basis set, the choice of the basis set to compute accurate DMA multipole values is not critical. METHODS The molecular charge densities were computed using the electronic structure methods Hartree-Fock (HF), MP2, MP4, DFT/PBE, DFT/B3LYP, DFT/B3PW91, DFT/M06-2X, and DFT/B2PLYP implemented in the Gaussian 09 package. MP4 was the benchmark method. The DMA multipoles were obtained with the GDMA program of Stone. The 6-311G + + (d,p) basis set was used for the production calculations, and the augmented correlation-consistent Dunning's hierarchy of basis sets was employed to evaluate the dependence of the DMA multipoles on the basis set size.
Collapse
Affiliation(s)
| | - Marco Aurélio Souza Oliveira
- Departamento de Química, Instituto Militar de Engenharia, Praça General Tibúrcio, 80, Rio de Janeiro, RJ, 22290-270, Brazil
| | - Itamar Borges
- Departamento de Química, Instituto Militar de Engenharia, Praça General Tibúrcio, 80, Rio de Janeiro, RJ, 22290-270, Brazil.
| |
Collapse
|
3
|
Zahariev F, Xu P, Westheimer BM, Webb S, Galvez Vallejo J, Tiwari A, Sundriyal V, Sosonkina M, Shen J, Schoendorff G, Schlinsog M, Sattasathuchana T, Ruedenberg K, Roskop LB, Rendell AP, Poole D, Piecuch P, Pham BQ, Mironov V, Mato J, Leonard S, Leang SS, Ivanic J, Hayes J, Harville T, Gururangan K, Guidez E, Gerasimov IS, Friedl C, Ferreras KN, Elliott G, Datta D, Cruz DDA, Carrington L, Bertoni C, Barca GMJ, Alkan M, Gordon MS. The General Atomic and Molecular Electronic Structure System (GAMESS): Novel Methods on Novel Architectures. J Chem Theory Comput 2023; 19:7031-7055. [PMID: 37793073 DOI: 10.1021/acs.jctc.3c00379] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/06/2023]
Abstract
The primary focus of GAMESS over the last 5 years has been the development of new high-performance codes that are able to take effective and efficient advantage of the most advanced computer architectures, both CPU and accelerators. These efforts include employing density fitting and fragmentation methods to reduce the high scaling of well-correlated (e.g., coupled-cluster) methods as well as developing novel codes that can take optimal advantage of graphical processing units and other modern accelerators. Because accurate wave functions can be very complex, an important new functionality in GAMESS is the quasi-atomic orbital analysis, an unbiased approach to the understanding of covalent bonds embedded in the wave function. Best practices for the maintenance and distribution of GAMESS are also discussed.
Collapse
Affiliation(s)
- Federico Zahariev
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50014, United States
| | - Peng Xu
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50014, United States
| | - Bryce M Westheimer
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50014, United States
| | - Simon Webb
- VeraChem LLC, 12850 Middlebrook Road, Suite 205, Germantown, Maryland 20874-5244, United States
| | - Jorge Galvez Vallejo
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50014, United States
- Research School of Computer Science, Australian National University, Canberra, ACT 2601, Australia
| | - Ananta Tiwari
- EP Analytics, Inc., 9909 Mira Mesa Boulevard, Suite 230, San Diego, California 92131, United States
| | - Vaibhav Sundriyal
- Department of Computational Modeling and Simulation Engineering, Old Dominion University, Norfolk, Virginia 23529, United States
| | - Masha Sosonkina
- Department of Computational Modeling and Simulation Engineering, Old Dominion University, Norfolk, Virginia 23529, United States
| | - Jun Shen
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
| | - George Schoendorff
- Propellants Branch, Rocket Propulsion Division, Aerospace Systems Directorate, Air Force Research Laboratory, AFRL/RQRP, Edwards Air Force Base, California 93524, United States
| | - Megan Schlinsog
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50014, United States
| | - Tosaporn Sattasathuchana
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50014, United States
| | - Klaus Ruedenberg
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50014, United States
| | - Luke B Roskop
- Hewlett-Packard Enterprise, 2131 Lindau Lane #1000, Bloomington, Minnesota 55425, United States
| | | | - David Poole
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50014, United States
- School of Chemistry & Biochemistry, Georgia Institute of Technology, Athens, Georgia 30332, United States
| | - Piotr Piecuch
- Department of Chemistry and Department of Physics and Astronomy, Michigan State University, East Lansing, Michigan 48824, United States
| | - Buu Q Pham
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50014, United States
| | - Vladimir Mironov
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Joani Mato
- Physical Sciences Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, MS K1-83, Richland, Washington 99352, United States
| | - Sam Leonard
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50014, United States
| | - Sarom S Leang
- EP Analytics, Inc., 9909 Mira Mesa Boulevard, Suite 230, San Diego, California 92131, United States
| | - Joe Ivanic
- Advanced Biomedical Computational Science, Frederick National Laboratory for Cancer Research, Frederick, Maryland 21702, United States
| | - Jackson Hayes
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50014, United States
| | - Taylor Harville
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50014, United States
| | - Karthik Gururangan
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
| | - Emilie Guidez
- Department of Chemistry, University of Colorado Denver, Denver, Colorado 80217, United States
| | - Igor S Gerasimov
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Christian Friedl
- Institut für Theoretische Physik, Johannes Kepler Universität Linz, Altenberger Str. 69, 4040 Linz, Austria
| | - Katherine N Ferreras
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50014, United States
| | - George Elliott
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50014, United States
| | - Dipayan Datta
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50014, United States
| | - Daniel Del Angel Cruz
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50014, United States
| | - Laura Carrington
- EP Analytics, Inc., 9909 Mira Mesa Boulevard, Suite 230, San Diego, California 92131, United States
| | - Colleen Bertoni
- Argonne Leadership Computing Facility, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Giuseppe M J Barca
- Research School of Computer Science, Australian National University, Canberra, ACT 2601, Australia
| | - Melisa Alkan
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50014, United States
| | - Mark S Gordon
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50014, United States
| |
Collapse
|
4
|
Li Y, Wang D, Fu F, Xia Q, Li W, Li S. Structures and properties of ionic crystals and condensed phase ionic liquids predicted with the generalized energy-based fragmentation method. J Comput Chem 2022; 43:704-716. [PMID: 35213748 DOI: 10.1002/jcc.26828] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2021] [Revised: 01/20/2022] [Accepted: 01/21/2022] [Indexed: 11/11/2022]
Abstract
The generalized energy-based fragmentation (GEBF) approach is extended to facilitate ab initio investigations of structures, lattice energies, vibrational spectra and 1 H NMR chemical shifts of ionic crystals and condensed-phase ionic liquids (ILs) with the periodic boundary conditions (PBC). For selected periodic systems, our results demonstrate that the so-called PBC-GEBF approach can provide satisfactory descriptions on ground-state energies, structures, and vibrational spectra of ionic crystals and IL crystals. The PBC-GEBF approach is then applied to three realistic condensed phase systems. For three ionic crystals (LiCl, NaCl, and KCl), we apply the PBC-GEBF approach with MP2 theory as well as some popular DFT methods to investigate their crystal structures and lattice energies. Our calculations indicate that the crystal structures obtained with PBC-GEBF-MP2/6-311 + G** are very close to the corresponding X-ray structures, while PBC-GEBF-ωB97X-D/6-311 + G** provides satisfactory prediction for crystal structures and lattice energies. For two polymorphs of [n-C4 mim][Cl] crystals, we find that the PBC-GEBF approach at the M06-2X/6-311 + G** level can give a satisfactory descriptions on structures and Raman spectra of these two crystals. Furthermore, for [C2 mim][BF4 ] ILs, we demonstrate that their 1 H NMR chemical shifts can be estimated from averaging over 5 typical snapshots (extracted from MD simulations) with the PBC-GEBF approach at the B97-2/pcSseg-2 level. The calculated results account for the observed experimental data quite well. Therefore, we expect that the PBC-GEBF approach, combined with various quantum chemistry methods, will become an effective tool in predicting structures and properties of ionic crystals and condensed-phase ILs.
Collapse
Affiliation(s)
- Yunzhi Li
- School of Chemistry and Chemical Engineering, Linyi University, Linyi, China.,School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing, China
| | - Dong Wang
- School of Chemistry and Chemical Engineering, Linyi University, Linyi, China
| | - Fangjia Fu
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing, China
| | - Qiying Xia
- School of Chemistry and Chemical Engineering, Linyi University, Linyi, China
| | - Wei Li
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing, China
| | - Shuhua Li
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing, China
| |
Collapse
|