To study the contribution of three-body dispersion to crystal lattice energies, we compute the three-body contributions to the lattice energies for crystalline benzene, carbon dioxide, and triazine using various computational methods. We show that these contributions converge quickly as the intermolecular distances between the monomers grow. In particular, the smallest value among the three pairwise intermonomer closest-contact distances, Rmin, shows a strong correlation with the three-body contribution to the lattice energy, and, here, the largest of the closest-contact distances, Rmax, serves as a cutoff criterion to limit the number of trimers to be considered. We considered all trimers up to Rmax=15Å. The trimers with Rmin<4Å contribute 90.4%, 90.6%, and 93.9% of the total three-body contributions for crystalline benzene, carbon dioxide, and triazine, respectively, for the coupled-cluster singles, doubles, and perturbative triples [CCSD(T)] method. For trimers with Rmin>4Å, the second-order Møller–Plesset perturbation theory (MP2) supplemented with the Axilrod–Teller–Muto (ATM) three-body dispersion correction reproduces the CCSD(T) values for the cumulative three-body contributions with errors of less than 0.1 kJ mol−1. Moreover, three-body contributions are converged within 0.15 kJ mol−1 by Rmax=10Å. From these results, it appears that in molecular crystals where dispersion dominates the three-body contribution to the lattice energy, the trimers with Rmin>4Å can be computed with the MP2+ATM method to reduce the computational cost, and those with Rmax>10Å appear to be basically negligible.
more »
« less
This content will become publicly available on June 21, 2024
Benchmark coupled-cluster lattice energy of crystalline benzene and assessment of multi-level approximations in the many-body expansion
The many-body expansion (MBE) is promising for the efficient, parallel computation of lattice energies in organic crystals. Very high accuracy should be achievable by employing coupled-cluster singles, doubles, and perturbative triples at the complete basis set limit [CCSD(T)/CBS] for the dimers, trimers, and potentially tetramers resulting from the MBE, but such a brute-force approach seems impractical for crystals of all but the smallest molecules. Here, we investigate hybrid or multi-level approaches that employ CCSD(T)/CBS only for the closest dimers and trimers and utilize much faster methods like Møller–Plesset perturbation theory (MP2) for more distant dimers and trimers. For trimers, MP2 is supplemented with the Axilrod–Teller–Muto (ATM) model of three-body dispersion. MP2(+ATM) is shown to be a very effective replacement for CCSD(T)/CBS for all but the closest dimers and trimers. A limited investigation of tetramers using CCSD(T)/CBS suggests that the four-body contribution is entirely negligible. The large set of CCSD(T)/CBS dimer and trimer data should be valuable in benchmarking approximate methods for molecular crystals and allows us to see that a literature estimate of the core-valence contribution of the closest dimers to the lattice energy using just MP2 was overbinding by 0.5 kJ mol−1, and an estimate of the three-body contribution from the closest trimers using the T0 approximation in local CCSD(T) was underbinding by 0.7 kJ mol−1. Our CCSD(T)/CBS best estimate of the 0 K lattice energy is −54.01 kJ mol−1, compared to an estimated experimental value of −55.3 ± 2.2 kJ mol−1.
more »
« less
- NSF-PAR ID:
- 10459879
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 158
- Issue:
- 23
- ISSN:
- 0021-9606
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Using the many-body expansion to predict crystal lattice energies (CLEs), a pleasantly parallel process, allows for flexibility in the choice of theoretical methods. Benchmark-level two-body contributions to CLEs of 23 molecular crystals have been computed using interaction energies of dimers with minimum inter-monomer separations (i.e., closest contact distances) up to 30 Å. In a search for ways to reduce the computational expense of calculating accurate CLEs, we have computed these two-body contributions with 15 different quantum chemical levels of theory and compared these energies to those computed with coupled-cluster in the complete basis set (CBS) limit. Interaction energies of the more distant dimers are easier to compute accurately and several of the methods tested are suitable as replacements for coupled-cluster through perturbative triples for all but the closest dimers. For our dataset, sub-kJ mol−1 accuracy can be obtained when calculating two-body interaction energies of dimers with separations shorter than 4 Å with coupled-cluster with single, double, and perturbative triple excitations/CBS and dimers with separations longer than 4 Å with MP2.5/aug-cc-pVDZ, among other schemes, reducing the number of dimers to be computed with coupled-cluster by as much as 98%.more » « less
-
more » « less
Abstract The global minima of urea and thiourea were characterized along with other low‐lying stationary points. Each structure was optimized with the CCSD(T) method and triple‐
ζ correlation consistent basis sets followed by harmonic vibrational frequency computations. Relative energies evaluated near the complete basis set limit with both canonical and explicitly correlated CCSD(T) techniques reveal several subtle but important details about both systems. These computations resolve a discrepancy by demonstrating that the electronic energy of the C2v second‐order saddle point of urea lies at least 1.5 kcal mol−1above the C2global minimum regardless of whether the structures were optimized with MP2, CCSD, or CCSD(T). Additionally, urea effectively has one minimum instead of two because the electronic barrier for inversion at one amino group in the Cslocal minimum vanishes at the CCSD(T) CBS limit. Characterization of both systems with the same ab initio methods and large basis sets conclusively establishes that the electronic barriers to inversion at one or both NH2groups in thiourea are appreciably smaller than in urea. CCSDT(Q)/cc‐pVTZ computations show higher‐order electron correlation effects have little impact on the relative energies and are consistently offset by core correlation effects of opposite sign and comparable magnitude. -
more » « less
Abstract The structures of zinc carbene ZnCH2and zinc carbyne HZnCH, and the conversion transition states between them are optimized at B3LYP/aug‐cc‐pVTZ, MP2/aug‐cc‐pVTZ, and CCSD/aug‐cc‐pVTZ levels of theory. The thermodynamic energies with CCSD(T) method are further extrapolated to basis set limit through a series of basis sets of aug‐cc‐pVXZ (X=D, T, Q, 5). The Zn−C bonding characteristics are interpreted by molecular plots, Laplacian of density plots, the integrated delocalization indices, net atomic charges, and derived atomic hardness. On the one hand, the studies demonstrated the efficiency of DFT method in structure optimizations and the accuracy of CBS method in obtaining thermodynamic energies; On the other hand, the density analysis of CCSD/aug‐cc‐pVDZ density demonstrates that both the sharing interaction and ionic interaction are important in ZnCH2ad HZnCH. The3B1state of ZnCH2is the global minimum and formed in visible light, but its small bond dissociation energy (47.0 kcal/mol) cannot keep the complex intact under UV light (79.4–102.1 kcal/mol). However, the3Σ−state of HZnCH can survive the UV light due to the greater Zn−C dissociation energy (100.7 kcal/mol). The delocalization indices of Zn…C in both3B1of ZnCH2(0.777) and the3Σ−state of HZnCH (0.785) are close to the delocalization index of the single C−C bond of ethane (0.841), i. e. the nomenclature of Zinc carbene and Zinc carbyne is incorrect. The stronger Zn−C bond in the3Σ−state of HZnCH than in the3B1state of ZnCH2can be attributed to the larger charge separation in the former. It was found that the Zn−C bonds in related Zinc organic compounds were also single bonds no matter whether the organic groups are CR, CR2, or CR3. The ionic interactions were discussed in terms of the atomic hardness that were in turn related to ionization energy and electron affinity. The unique combination of covalent and ionic characteristics in the Zn−C bonds of organic Zinc compounds could be the origin of many interesting applications of organic Zinc reagents.
-
more » « less
The optimized geometries, vibrational frequencies, and dissociation energies from MP2 and CCSD(T) computations with large correlation consistent basis sets are reported for (H2S)2and H2O/H2S. Anharmonic vibrational frequencies have also been computed with second‐order vibrational perturbation theory (VPT2). As such, the fundamental frequencies, overtones, and combination bands reported in this study should also provide a useful road map for future spectroscopic studies of the simple but important heterogeneous H2O/H2S dimer in which the hydrogen bond donor and acceptor can interchange, leading to two unique minima with very similar energies. Near the CCSD(T) complete basis set limit, the HOH⋯SH2configuration (H2O donor) lies only 0.2 kcal mol−1below the HSH⋯OH2structure (H2S donor). When the zero‐point vibrational energy is included, however, the latter configuration becomes slightly lower in energy than the former by <0.1 kcal mol−1. © 2018 Wiley Periodicals, Inc.
