An official website of the United States government
Quasi-harmonic approaches provide an economical route to modeling the temperature dependence of molecular crystal structures and properties. Several studies have demonstrated good performance of these models, at least for rigid molecules, when using fragment-based approaches with correlated wavefunction techniques. Many others have found success employing dispersion-corrected density functional theory (DFT). Here, a hierarchy of models in which the energies, geometries, and phonons are computed either with correlated methods or DFT are examined to identify which combinations produce useful predictions for properties such as the molar volume, enthalpy, and entropy as a function of temperature. The results demonstrate that refining DFT geometries and phonons with single-point energies based on dispersion-corrected second-order M{\o}ller-Plesset perturbation theory can provide clear improvements in the molar volumes and enthalpies compared to those obtained from DFT alone. Predicted entropies, which are governed by vibrational contributions, benefit less clearly from the hybrid schemes. Using these hybrid techniques, the room-temperature thermochemistry of acetaminophen (paracetamol) is predicted to address the discrepancy between two experimental sublimation enthalpy measurements.
more »
« less
Ab initio thermodynamic properties and their uncertainties for crystalline α-methanol
To investigate the performance of quasi-harmonic electronic structure methods for modeling molecular crystals at finite temperatures and pressures, thermodynamic properties are calculated for the low-temperature α polymorph of crystalline methanol. Both density functional theory (DFT) and ab initio wavefunction techniques up to coupled cluster theory with singles, doubles, and perturbative triples (CCSD(T)) are combined with the quasi-harmonic approximation to predict energies, structures, and properties. The accuracy, reliability, and uncertainties of the individual quantum-chemical methods are assessed via detailed comparison of calculated and experimental data on structural properties (density) and thermodynamic properties (isobaric heat capacity). Performance of individual methods is also studied in context of the hierarchy of the quantum-chemical methods. The results indicate that while some properties such as the sublimation enthalpy and thermal expansivity can be modeled fairly well, other properties such as the molar volume and isobaric heat capacities are harder to predict reliably. The errors among the energies, structures, and phonons are closely coupled, and most accurate predictions here appear to arise from fortuitous error compensation among the different contributions. This study highlights how sensitive molecular crystal property predictions can be to the underlying model approximations and the significant challenges inherent in first-principles predictions of solid state structures and thermochemistry.
more »
« less
- Award ID(s):
- 1665212
- PAR ID:
- 10062129
- Date Published:
- Journal Name:
- Physical Chemistry Chemical Physics
- Volume:
- 19
- Issue:
- 44
- ISSN:
- 1463-9076
- Page Range / eLocation ID:
- 29940 to 29953
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)The energy of the lowest-lying triplet state (T1) relative to the ground and first-excited singlet states (S0, S1) plays a critical role in optical multiexcitonic processes of organic chromophores. Focusing on triplet–triplet annihilation (TTA) upconversion, the S0 to T1 energy gap, known as the triplet energy, is difficult to measure experimentally for most molecules of interest. Ab initio predictions can provide a useful alternative, however low-scaling electronic structure methods such as the Kohn–Sham and time-dependent variants of Density Functional Theory (DFT) rely heavily on the fraction of exact exchange chosen for a given functional, and tend to be unreliable when strong electronic correlation is present. Here, we use auxiliary-field quantum Monte Carlo (AFQMC), a scalable electronic structure method capable of accurately describing even strongly correlated molecules, to predict the triplet energies for a series of candidate annihilators for TTA upconversion, including 9,10 substituted anthracenes and substituted benzothiadiazole (BTD) and benzoselenodiazole (BSeD) compounds. We compare our results to predictions from a number of commonly used DFT functionals, as well as DLPNO-CCSD(T 0 ), a localized approximation to coupled cluster with singles, doubles, and perturbative triples. Together with S1 estimates from absorption/emission spectra, which are well-reproduced by TD-DFT calculations employing the range-corrected hybrid functional CAM-B3LYP, we provide predictions regarding the thermodynamic feasibility of upconversion by requiring (a) the measured T1 of the sensitizer exceeds that of the calculated T1 of the candidate annihilator, and (b) twice the T1 of the annihilator exceeds its S1 energetic value. We demonstrate a successful example of in silico discovery of a novel annihilator, phenyl-substituted BTD, and present experimental validation via low temperature phosphorescence and the presence of upconverted blue light emission when coupled to a platinum octaethylporphyrin (PtOEP) sensitizer. The BTD framework thus represents a new class of annihilators for TTA upconversion. Its chemical functionalization, guided by the computational tools utilized herein, provides a promising route towards high energy (violet to near-UV) emission.more » « less
-
Ternary sulfides and selenides in the distorted-perovskite structure (“chalcogenide perovskites”) are predicted by theory to be semiconductors with a band gap in the visible-to-infrared and may be useful for optical, electronic, and energy conversion technologies. Here we use computational thermodynamics to predict the pressure–temperature phase diagrams for select chalcogenide perovskites. Our calculations incorporate formation energies calculated by density functional theory, and empirical estimates of heat capacities. We highlight the windows of thermodynamic equilibrium between solid chalcogenide perovskites and the vapor phase at high temperature and very low pressure. These results can guide the adsorption-limited growth of ternary chalcogenides by molecular beam epitaxy.more » « less
-
Bonding energies play an essential role in describing the relative stability of molecules in chemical space. Therefore, methods employed to search chemical space need to capture the bonding behavior for a wide range of molecules, including radicals. In this work, we investigate the ability of quantum alchemy to capture the bonding behavior of hypothetical chemical compounds, specifically diatomic molecules involving hydrogen with various electronic structures. We evaluate equilibrium bond lengths, ionization energies, and electron affinities of these fundamental systems. We compare and contrast how well manual quantum alchemy calculations, i.e., quantum mechanics calculations in which the nuclear charge is altered, and quantum alchemy approximations using a Taylor series expansion can predict these molecular properties. Our results suggest that while manual quantum alchemy calculations outperform Taylor series approximations, truncations of Taylor series approximations after the second order provide the most accurate Taylor series predictions. Furthermore, these results suggest that trends in quantum alchemy predictions are generally dependent on the predicted property (i.e., equilibrium bond length, ionization energy, or electron affinity). Taken together, this work provides insight into how quantum alchemy predictions using a Taylor series expansion may be applied to future studies of non-singlet systems as well as the challenges that remain open for predicting the bonding behavior of such systems.more » « less
-
The phenalene (triangulene) and olympicene molecules belong to the polycyclic aromatic hydrocarbons (PAHs) class, which attracted substantial technological interest due to their unique electronic properties. Electronic structure calculations serve as a valuable tool in investigating the stability and reactivity of these molecular systems. In the present work, the multireference calculations, namely the complete active space second-order perturbation theory (CASPT2) and multireference averaged quadratic coupled cluster (MR-AQCC), were employed to study the reactivity and stability of phenalene and olympicene isomers, as well as their modified structures where the sp3-carbon at the borders were removed. The harmonic oscillator model of aromaticity (HOMA) and the nucleus-independent chemical shift (NICS) as geometric and magnetic indexes calculated with density functional theory were utilized to assess the aromaticity of the studied molecules. These indexes were compared with properties such as the excitation energy and natural orbitals occupation. The reactivity analyzed using the HOMA index combined with MR-AQCC revealed the radical character of certain structures, as well as the weakening of their aromaticity. Moreover, the results suggest that the removal of sp3-carbon atoms and the addition of hydrogen atoms did not alter the π network and the excitation energies of the phenalene molecules.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1039/C7CP06605H
