Computational materials discovery efforts are enabled by large databases of properties derived from high-throughput density functional theory (DFT), which now contain millions of calculations at the generalized gradient approximation (GGA) level of theory. It is now feasible to carry out high-throughput calculations using more accurate methods, such as meta-GGA DFT; however recomputing an entire database with a higher-fidelity method would not effectively leverage the enormous investment of computational resources embodied in existing (GGA) calculations. Instead, we propose here a general procedure by which higher-fidelity, low-coverage calculations (e.g., meta-GGA calculations for selected chemical systems) can be combined with lower-fidelity, high-coverage calculations (e.g., an existing database of GGA calculations) in a robust and scalable manner. We then use legacy PBE(+
- Award ID(s):
- 1665212
- NSF-PAR ID:
- 10098716
- Date Published:
- Journal Name:
- Physical Chemistry Chemical Physics
- ISSN:
- 1463-9076
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract U ) GGA calculations and new r2SCAN meta-GGA calculations from the Materials Project database to demonstrate that our scheme improves solid and aqueous phase stability predictions, and discuss practical considerations for its implementation. -
Ices in the interstellar medium largely exist as amorphous solids composed of small molecules including ammonia, water, and carbon dioxide. Describing gas-phase molecules can be readily accomplished with current high-level quantum chemical calculations with the description of crystalline solids becoming more readily accomplished. Differently, amorphous solids require more novel approaches. The present work describes a method for generating amorphous structures and constructing electronic spectra through a combination of quantum chemical calculations and statistical mechanics. The structures are generated through a random positioning program and DFT methods, such as ωB97-XD and CAM-B3LYP. A Boltzmann distribution weights the excitations to compile a final spectrum from a sampling of molecular clusters. Three ice analogs are presented herein consisting of ammonia, carbon dioxide, and water. Ammonia and carbon dioxide provide semi-quantitative agreement with experiment for CAM-B3LYP/6-311++G(2d,2p) from 30 clusters of 8 molecules. Meanwhile, the amorphous water description improves when the sample size is increased in cluster size and count to as many as 105 clusters of 32 water molecules. The described methodology can produce highly comparative descriptions of electronic spectra for ice analogs and can be used to predict electronic spectra for other ice analogs.more » « less
-
Terahertz vibrational spectroscopy has emerged as a powerful spectroscopic technique, providing valuable information regarding long-range interactions – and associated collective dynamics – occurring in solids. However, the terahertz sciences are relatively nascent, and there have been significant advances over the last several decades that have profoundly influenced the interpretation and assignment of experimental terahertz spectra. Specifically, because there do not exist any functional group or material-specific terahertz transitions, it is not possible to interpret experimental spectra without additional analysis, specifically, computational simulations. Over the years simulations utilizing periodic boundary conditions have proven to be most successful for reproducing experimental terahertz dynamics, due to the ability of the calculations to accurately take long-range forces into account. On the other hand, there are numerous reports in the literature that utilize gas phase cluster geometries, to varying levels of apparent success. This perspective will provide a concise introduction into the terahertz sciences, specifically terahertz spectroscopy, followed by an evaluation of gas phase and periodic simulations for the assignment of crystalline terahertz spectra, highlighting potential pitfalls and good practice for future endeavors.more » « less
-
more » « less
Understanding, predicting, and ultimately controlling exciton band structure and exciton dynamics are central to diverse chemical and materials problems. Here, we have developed a first-principles method to determine exciton dispersion and exciton–phonon interaction in semiconducting and insulating solids based on time-dependent density functional theory. The first-principles method is formulated in planewave bases and pseudopotentials and can be used to compute exciton band structures, exciton charge density, ionic forces, the non-adiabatic coupling matrix between excitonic states, and the exciton–phonon coupling matrix. Based on the spinor formulation, the method enables self-consistent noncollinear calculations to capture spin-orbital coupling. Hybrid exchange-correlation functionals are incorporated to deal with long-range electron–hole interactions in solids. A sub-Hilbert space approximation is introduced to reduce the computational cost without loss of accuracy. For validations, we have applied the method to compute the exciton band structure and exciton–phonon coupling strength in transition metal dichalcogenide monolayers; both agree very well with the previous GW-Bethe–Salpeter equation and experimental results. This development paves the way for accurate determinations of exciton dynamics in a wide range of solid-state materials.
-
more » « less
Abstract Developing an atomistic understanding of ionizing radiation induced changes to organic materials is necessary for intentional design of greener and more sustainable materials for radiation shielding and detection. Cocrystals are promising for these purposes, but a detailed understanding of how the specific intermolecular interactions within the lattice upon exposure to radiation affect the structural stability of the organic crystalline material is unknown. This study evaluates atomistic‐level effects of γ radiation on both single‐ and multicomponent organic crystalline materials and how specific noncovalent interactions and packing within the crystalline lattice enhance structural stability. Dose studies were performed on all crystalline systems and evaluated via experimental and computational methods. Changes in crystallinity were evaluated by p‐XRD and free radical formation was analyzed via EPR spectroscopy. Type of intermolecular interactions and packing within the crystal lattice was delineated and related to the specific free radical species formed and the structural integrity of each material. Periodic DFT and HOMO‐LUMO surface mapping calculations provided atomistic‐level identifications of the most probable sites for the radicals formed upon exposure to γ radiation and relate intermolecular interactions and molecular packing within the crystalline lattice to experimental results.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1039/C9CP01666J
