skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Towards reliable ab initio sublimation pressures for organic molecular crystals – Are we there yet?
Knowledge of molecular crystal sublimation equilibrium data is vital in many industrial processes, but this data can be difficult to measure experimentally for low-volatility species. Theoretical prediction of sublimation pressures could provide a useful supplement to experiment, but the exponential temperature dependence of sublimation (or any saturated vapor) pressure curve makes this challenging. An uncertainty of only a few percent in the sublimation enthalpy or entropy can propagate to an error in the sublimation pressure exceeding several orders of magnitude for a given temperature interval. Despite this fundamental difficulty, this paper performs some of the first ab initio predictions of sublimation pressure curves. Four simple molecular crystals (ethane, methanol, benzene, and imidazole) have been selected for a case study showing the currently achievable accuracy of quantum chemistry calculations. Fragment-based ab initio techniques and the quasi-harmonic approximation are used for calculations of cohesive and phonon properties of the crystals, while the vapor phase is treated by the ideal gas model. Ab initio sublimation pressure curves for model compounds are compared against their experimental counterparts. The computational uncertainties are estimated, weak points of the computational methodology are identified, and further improvements are proposed.  more » « less
Award ID(s):
1665212
PAR ID:
10098715
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Physical Chemistry Chemical Physics
ISSN:
1463-9076
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Recent advances regarding the interplay between ab initio calculations and metrology are reviewed, with particular emphasis on gas-based techniques used for temperature and pressure measurements. Since roughly 2010, several thermophysical quantities – in particular, virial and transport coefficients – can be computed from first principles without uncontrolled approximations and with rigorously propagated uncertainties. In the case of helium, computational results have accuracies that exceed the best experimental data by at least one order of magnitude and are suitable to be used in primary metrology. The availability of ab initio virial and transport coefficients contributed to the recent SI definition of temperature by facilitating measurements of the Boltzmann constant with unprecedented accuracy. Presently, they enable the development of primary standards of thermodynamic temperature in the range 2.5–552 K and pressure up to 7 MPa using acoustic gas thermometry, dielectric constant gas thermometry, and refractive index gas thermometry. These approaches will be reviewed, highlighting the effect of first-principles data on their accuracy. The recent advances in electronic structure calculations that enabled highly accurate solutions for the many-body interaction potentials and polarizabilities of atoms – particularly helium – will be described, together with the subsequent computational methods, most often based on quantum statistical mechanics and its path-integral formulation, that provide thermophysical properties and their uncertainties. Similar approaches for molecular systems, and their applications, are briefly discussed. Current limitations and expected future lines of research are assessed. 
    more » « less
  2. 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
  3. This computational study characterises charge-transfer-to-solvent (CTTS) states of aqueous thiocyanate anion using equation-of-motion coupled-cluster methods combined with electrostatic embedding quantum mechanics/molecular mechanics (QM/MM) scheme. Equilibrium sampling was carried out using classical molecular dynamics (MD) with standard force-fields and QM/MM ab initio molecular dynamics (AIMD) using density functional theory. The two calculations yield significantly different local structure around solvated SCN− . Because of the diffuse character of CTTS states, they are very sensitive to the local structure of solvent around the solute and its dynamic fluctuations. Owing to this sensitivity, the spectra computed using MD and AIMD based snapshots differ considerably. This sensitivity suggests that the spectroscopy exploiting CTTS transitions can provide an experimental handle for assessing the quality of force-fields and density functionals. By combining CTTS-based spectroscopies with reliable theoretical modeling, detailed microscopic information of the solvent structure can be obtained. We present a robust computational protocol for modeling spectra of solvated anions and emphasise the use of an ab initio characterization of individual electronic transitions as CTTS or local excitations. 
    more » « less
  4. A molecular crystal structure prediction (CSP) protocol used in the seventh blind test is presented. The seventh blind test was divided into two stages and included seven targets, with crystals containing from one to three molecules in asymmetric units, monomers built of up to 100 atoms, and all targets containing monomers with flexible degrees of freedom. Some targets were cocrystals and one target was a salt. These diverse targets were treated using a CSP protocol starting from finding the global and local minima conformations of the target molecule. Subsequently, anab initiotwo-body rigid-monomer six-dimensional force field (aiFF) was developed for the global-minimum conformer. These aiFFs were then used in CSPs consisting of packing and lattice-energy minimization stages. Flexible-monomer CSPs were used for some targets. To describe the intramonomer FF, either generic empirical FFs or reparametrized FFs of this type were used, with some parameters fitted toab initioenergies of monomers in the latter case. A novel packing procedure was applied for two targets in stage 1. The success rate in the structure generation stage was 15% in submission phase and 54% in post-submission phase, while the corresponding values in the structure rating stage were 33% and 89%. We conclude that the inexpensive conformer-based approach with rigid-monomer CSPs can be recommended for investigations of crystals with flexible monomers. An advantage of this protocol is that it is fully based on first-principles quantum mechanics and generates tailor-made FFs suitable for use in subsequent molecular dynamics simulations investigating temperature-dependent effects. However, empirical intramonomer FFs reparametrized usingab initiodata are not yet adequate for CSPs. 
    more » « less
  5. Unlike naturally occurring oxide crystals such as ruby and gemstones, there are no naturally occurring nitride crystals because the triple bond of the nitrogen molecule is one of the strongest bonds in nature. Here, we report that when the transition metal scandium is subjected to molecular nitrogen, it self-catalyzes to break the nitrogen triple bond to form highly crystalline layers of ScN, a semiconductor. This reaction proceeds even at room temperature. Self-activated ScN films have a twin cubic crystal structure, atomic layering, and electronic and optical properties comparable to plasma-based methods. We extend our research to showcase Sc’s scavenging effect and demonstrate self-activated ScN growth under various growth conditions and on technologically significant substrates, such as 6H–SiC, AlN, and GaN. Ab initio calculations elucidate an energetically efficient pathway for the self-activated growth of crystalline ScN films from molecular N2. The findings open a new pathway to ultralow-energy synthesis of crystalline nitride semiconductor layers and beyond. 
    more » « less