Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Molecular crystals”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

How Accurate Are Simulations and Experiments for the Lattice Energies of Molecular Crystals?

Molecular crystals play a central role in a wide range of scientific fields, including pharmaceuticals and organic semiconductor devices. However, they are challenging systems to model accurately with computational approaches because of a delicate interplay of intermolecular interactions such as hydrogen bonding and Van der Waals dispersion forces. Here, by exploiting recent algorithmic developments, we report the first set of diffusion Monte Carlo lattice energies for all 23 molecular crystals in the popular and widely used X23 dataset. Comparisons with previous state-of-the-art lattice energy predictions (on a subset of the dataset) and a careful analysis of experimental sublimation enthalpies reveals that high-accuracy computational methods are now at least as reliable as (computationally derived) experiments for the lattice energies of molecular crystals. Overall, this work demonstrates the feasibility of high-level explicitly correlated electronic structure methods for broad benchmarking studies in complex condensed phase systems, and signposts a route towards closer agreement between experiment and simulation. Published by the American Physical Society 2024

Physics↗
DOE OSTI2023

How Accurate Can Crystal Structure Predictions Be for High-Energy Molecular Crystals?

Molecular crystals have shallow potential energy landscapes, with multiple local minima separated by very small differences in total energy. Predicting molecular packing and molecular conformation in the crystal generally requires ab initio methods of high accuracy, especially when polymorphs are involved. We used dispersion-corrected density functional theory (DFT-D) to assess the capabilities of an evolutionary algorithm (EA) for the crystal structure prediction (CSP) of well-known but challenging high-energy molecular crystals (HMX, RDX, CL-20, and FOX-7). While providing the EA with the experimental conformation of the molecule quickly re-discovers the experimental packing, it is more realistic to start instead from a naïve, flat, or neutral initial conformation, which reflects the limited experimental knowledge we generally have in the computational design of molecular crystals. By doing so, and using fully flexible molecules in fully variable unit cells, we show that the experimental structures can be predicted in fewer than 20 generations. Nonetheless, one must be aware that some molecular crystals have naturally hindered evolutions, requiring as many attempts as there are space groups of interest to predict their structures, and some may require the accuracy of all-electron calculations to discriminate between closely ranked structures. To save resources in this computationally demanding process, we showed that a hybrid xTB/DFT-D approach could be considered in a subsequent study to push the limits of CSP beyond 200+ atoms and for cocrystals.

42 ENGINEERING↗
DOE OSTI2026

Diabatization with Electrostatic Embedding for Studying Photophysics in Organic Molecular Crystals

Highly emissive organic molecular crystals find applications in several areas, such as organic electronics, solar cells, and sensors. Understanding the excited-state mechanisms underlying these applications is essential for optimizing and controlling them effectively. Exciton models coupled with nonadiabatic dynamics, particularly quantum dynamics, provide crucial insights into photochemical and photophysical processes in molecular crystals. Nevertheless, there remains a lack of general tools and automated workflows to facilitate such simulations. In this paper, we present a computational strategy to investigate the photoactivated dynamics of organic molecular crystals, bridging methodologies traditionally used for molecular systems and materials science, with a particular focus on the interplay between local excitations and charge transfer (CT) processes. We have implemented an interface between the fromage and Overdia programs, enabling the construction of vibronic Hamiltonians for molecular crystals within an excited-state ONIOM(QM:QM′) framework, incorporating long-range electrostatics through a RESP-based Ewald summation. Fragment-based diabatization provides a route to quantum dynamics simulations in weak-to-intermediate coupling regimes. The method was applied to the photophysics of dibenzo[g,p]chrysene (DBC) crystals using time-dependent DFT. The fromage/ Overdia interface was employed to compute the couplings of local excitations and CT states for 18 unique DBC dimers in the crystal and to quantify the influence of electrostatic embedding, which was found to be modest (10−20%). Simulations on π-stacked dimers reproduced the small red shift observed experimentally from solution to crystal, attributed to electronic interactions among fixed monomers rather than crystal electrostatics. Quantum dynamics simulations revealed ultrafast population transfer from bright local excitations to CT states. This approach establishes a robust framework linking molecular and solid-state excited-state dynamics, with potential applications for studying excitations, defects, and impurities in molecular crystals.

Crystals↗
DOE OSTI2026

Nanoconfined Grain Boundaries Increase the Conductivity of Polycrystalline Molecular Crystals

Soft-solid molecular crystals consist of crystalline grains and fluid grain boundaries (GBs) that enhance the grain binding and transport of Li + ions between the grains. The total ionic conductivity consists of ion migration in both the grains and GBs. To unravel these contributions in adiponitrile (Adpn):LiPF 6 molecular crystals, the GB volume fraction was varied by changing the size of the crystals and the Adpn:LiPF 6 molar ratio. Molecular dynamics (MD) simulations indicate that ion motion was subdiffusive in the grains and “well-diffusive” in the GBs, with GBs characterized as disordered nanoconfined regions of higher charge carrier concentration (∼1 M) than in saturated Adpn:LiPF 6 solutions (0.04 M), and Li + ions predominantly solvated by cyano groups with few contact ion pairs. The diffusivity in the GBs is at least an order of magnitude higher than that in the crystalline grains. The emergent picture is the grains as a reservoir of ions that migrate to faster-conducting GBs.

Crystals↗
DOE OSTI2020

Genarris 2.0: A Random Structure Generator for Molecular Crystals

Genarris is an open source Python package for generating random molecular crystal structures with physical constraints for seeding crystal structure prediction algorithms and training machine learning models. Here we present a new version of the code, containing several major improvements. A MPI-based parallelization scheme has been implemented, which facilitates the seamless sequential execution of user-defined workflows. A new method for estimating the unit cell volume based on the single molecule structure has been developed using a machine-learned model trained on experimental structures. A new algorithm has been implemented for generating crystal structures with molecules occupying special Wyckoff positions. A new hierarchical structure check procedure has been developed to detect unphysical close contacts efficiently and accurately. New intermolecular distance settings have been implemented for strong hydrogen bonds. To demonstrate these new features, we study two specific cases: benzene and glycine. Genarris finds the experimental structures of the two polymorphs of benzene and the three polymorphs of glycine. Program summary Program Title: Genarris 2.0 Program Files doi: http://dx.doi.org/10.17632/grx6mz4pjn.1 Licensing provisions: BSD-3 Clause Programming language: Python, C External routines/libraries: Spglib, ASE, pymatgen, SciPy, mpi4py, scikit-learn, PyTorch, FHI-aims. Nature of problem: Molecular crystal structure prediction. Solution method: Genarris 2.0 generates molecular crystal structures over the 230 space groups, on general and special Wyckoff positions, using physical constraints. Down-sampling of the generated structures may be performed subsequently, based on molecular crystal packing descriptors and an unsupervised machine learning algorithm. Lastly, ab initio structure relaxation may be performed for the final pool. Depending on the user-defined workflow implemented, Genarris may be used to generate diverse molecular crystal datasets to seed evolutionary algorithms or to train machine learning algorithms or as a standalone crystal structure prediction method. Restrictions: For crystal structure generation, the molecule of interest must be semi-rigid with no bond rotational degrees of freedom. Unusual features: Genarris 2.0 is a highly distributed program, making use of MPI for Python parallelization. The user has the ability to design and implement workflows by executing a user-defined list of procedures. Genarris 2.0 offers new features including a machine learning model for estimating the molecular volume in the solid state from the single molecule structure, structure generation in special Wyckoff positions of space groups, hierarchical structure checks including rigorous treatment of non-orthogonal structures, and clustering and down-selection workflows combining first principles simulations with machine learning. (C) 2020 Elsevier B.V. All rights reserved.

Crystal structure prediction↗
DOE OSTI2025

PAH101: A GW+BSE Dataset of 101 Polycyclic Aromatic Hydrocarbon (PAH) Molecular Crystals

Abstract The excited-state properties of molecular crystals are important for applications in organic electronic devices. TheGWapproximation and Bethe-Salpeter equation (GW+BSE) is the state-of-the-art method for calculating the excited-state properties of crystalline solids with periodic boundary conditions. We present the PAH101 dataset ofGW+BSE calculations for 101 molecular crystals of polycyclic aromatic hydrocarbons (PAHs) with up to ~500 atoms in the unit cell. To the best of our knowledge, this is the firstGW+BSE dataset for molecular crystals. The data records include theGWquasiparticle band structure, the fundamental band gap, the static dielectric constant, the first singlet exciton energy (optical gap), the first triplet exciton energy, the dielectric function, and optical absorption spectra for light polarized along the three lattice vectors. The dataset can be used to (i) discover materials with desired electronic/optical properties, (ii) identify correlations between DFT andGW+BSE quantities, and (iii) train machine learned models to help in materials discovery efforts.

Science & Technology - Other Topics↗
DOE OSTI2022

Can a shock-induced phonon up-pumping model relate to impact sensitivity of molecular crystals, polymorphs and cocrystals?

Impact sensitivity engineering of high-energy molecular crystals requires accurate predictive models. For this purpose, the promising multi-phonon based approach is selected, assessing a bit more its strengths and weaknesses. Presently used with high-quality phonon calculations of 22 molecular crystals, using a physics-based criterion to determine the phonon bath extent, the resulting intrinsic shock sensitivity index (SSI) is compared to the most common marker of impact sensitivity, h 50 , as determined from drop-weight impact tests. Selecting a data subset from experiments performed under very similar conditions (2.5 kg hammer with grit and 30–40 mg samples), the model can predict h 50 values for mono-molecular crystals with very good accuracy, including the ability to discriminate the polymorphs of HMX and CL20. This very good agreement validates an initial indirect up-pumping mechanism occurring under these conditions, where the doorway modes also interact with the phonon bath. However, the phonon bath criterion for mono-molecular crystals does not transfer well to cocrystals. Owing to the vibrational coupling of the co-molecules, it seems a broader phonon bath should be considered. Additionally recalling experimental uncertainty and various experimental factors affecting h 50 values for a given compounds, we recommend that the density of the sample, granularity and morphology be systematically considered and reported along with measurements, which will in turn allow for more systematic data and predictive capabilities for sensitivity models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Quantum Vibronic Effects on the Electronic Properties of Molecular Crystals

We present a study of molecular crystals, focused on the effect of nuclear quantum motion and anharmonicity on their electronic properties. We consider a system composed of relatively rigid molecules, a diamondoid crystal, and one composed of floppier molecules, NAI-DMAC, a thermally activated delayed fluorescence compound. We compute fundamental electronic gaps at the density functional theory (DFT) level of theory, with the Perdew–Burke–Erzenhof (PBE) and strongly constrained and approximately normed (SCAN) functionals, by coupling first-principles molecular dynamics with a nuclear quantum thermostat. We find a sizable zero-point renormalization (ZPR) of the band gaps, which is much larger in the case of diamondoids (0.6 eV) than for NAI-DMAC (0.22 eV). We show that the frozen phonon (FP) approximation, which neglects intermolecular anharmonic effects, leads to a large error (~50%) in the calculation of the band gap ZPR. Instead, when using a stochastic method, we obtain results in good agreement with those of our quantum simulations for the diamondoid crystal. However, the agreement is worse for NAI-DMAC where intramolecular anharmonicities contribute to the ZPR. Our results highlight the importance of accurately including nuclear and anharmonic quantum effects to predict the electronic properties of molecular crystals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Relative molecular orientation can impact the onset of plasticity in molecular crystals

Abstract Creating or moving dislocations is the first step to dissipating mechanical energy via plastic deformation under contact loading. In molecular crystals there is both a lattice that defines crystal orientation and a relative orientation of the basis of the molecules. We define a normalization parameter which relates strain at yield, the hardness of the bulk crystal, and a distance parameter analogous to a Burgers vector that nominally predicts the relative ease of initiating plasticity in this broad class of materials. Analyzing the yield behavior of 10 different molecular crystals of varying space groups shows the inter-molecular orientation predicts the experimentally observed applied stress needed to nucleate dislocations. When molecules are oriented ‘parallel’ relative to one another the normalized maximum shear stress at the onset of plasticity is on the order of 3–5 times lower than when molecules within the crystal are ‘anti-parallel’, and molecules with a more equiaxed shape fall in between these bounds. This provides an initial indication of a structural feature which predicts the relative ease of initiating plasticity during contact loading in molecular crystals.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Quantification of Dynamic Scattering Effects in Molecular Crystals using Large Angle Rocking Beam Electron Diffraction

Electron crystallography provides a pathway to solve structure of small crystals (< 1um) in size, and thus overcomes difficult synthesis constraints involved in growing large crystals. Generally, electron diffraction data is collected in the form of integrated intensity using continuous rotation or by precession. The measured intensities in 3D are utilized for structure solution. Using this approach, structure solution of difficult crystals, such as small crystals of zeolites, metal-organic frameworks, molecular crystals, and proteins , can be solved by electron diffraction. However, electron structure solutions are regularly reported with higher Rvalues than x-ray or neutron diffraction. While similar structures are found despite the high R-values, the differences in the measure intensity and theory calculated intensity limit information that can be extracted by electron diffraction. Previous work demonstrated that including multiple-scattering effects significantly reduce the R-values . Thus, it is critical to be able to quantify the dynamical diffraction effects in molecular crystals.

Crystallography↗
DOE OSTI2026

Environment-imposed selection rules for nuclear-spin conversion of H 2 in molecular crystals

Nuclear-spin conversion in molecular hydrogen is governed by strict symmetry rules that typically require magnetic fields or catalytic surfaces to break. Here we demonstrate that the intrinsic tensor composition of a non-magnetic molecular crystal field can impose and relax these rules without external fields. High-resolution infrared spectra of H 2 in crystalline CO 2 reveal large rank-2 (quadrupolar) crystal-field splittings of the 𝑚 sublevels, while nuclear-spin conversion occurs only through 𝛥𝑚 = 0 channels. Replacing CO 2 with polar N 2 O introduces rank-1 (dipole) components that partially open 𝛥⁢𝑚 ≠ 0 pathways, while incorporation of paramagnetic NO 2 fully lifts the restriction. Furthermore, these results establish a direct correspondence between crystal-field tensor rank and nuclear-spin dynamics, introducing a general symmetry-based framework for designing and controlling spin-isomer populations and quantum-state connectivity in molecular solids.

McLane, Nathan [University of Maryland, College Pa↗
DOE OSTI2020

Automated Identification of Molecular Crystals’ Packing Motifs

Packing motifs—patterns in how molecules orient relative to one another in a crystal structure—are an important concept in many subdisciplines of materials science because of correlations observed between specific packing motifs and properties of interest. That said, packing motif data sets have remained small and noisy due to intensive manual labeling processes and insufficient labeling schemes. The most prominent labeling algorithms calculate relative interplanar angles of nearest neighbor molecules to determine the packing motif of a molecular crystal, but this simple approach can fail when neighbors are naively sampled isotropically around the crystal structure. To remedy this issue, here we propose an optimization algorithm, which rotates the molecular crystal structure to find representative molecules that inform the packing motif. We package this algorithm into an automated framework—Autopack—which both optimally rotates the crystal structure and labels the packing motif based on the appropriate neighboring molecules. In this work, we detail the Autopack framework and its performance, which shows improvements compared to previous state-of-the-art labeling methods, providing the first quantitative point of comparison for packing motif labeling algorithms. Furthermore, using Autopack (available at https://ipo.llnl.gov/technologies/software/autopack), we perform the first large-scale study of potential relationships between chemicals’ compositions and packing motifs, which shows that these relationships are more complex than previously hypothesized from studies that used only tens of polycyclic aromatic hydrocarbon molecules. Autopack’s capabilities help pose next steps for crystal engineering research focusing not only on a molecule’s adoption of a specific packing motif but also on new structure–property relationships.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

The mechanical response of micron-sized molecular crystals

Abstract Microstructures and corresponding properties of compacted powders ultimately depend on the mechanical response of individual particles. In principle, computational simulations can predict the results of powder compaction processes, but the selection of appropriate models for both particle–particle interactions and particle deformations across all relevant length scales remain nontrivial tasks, especially in material systems lacking detailed mechanical property information. The work presented here addresses these issues by conducting uniaxial compressions in situ inside of a scanning electron microscope to characterize the mechanical response of individual micron-sized particles of a molecular crystal, hexanitrohexaazaisowurtzitane (CL-20). This experimental approach enabled the collection of quantitative force and displacement data alongside simultaneous imaging to capture morphology changes. The results reveal information about elastic deformation, yield, plastic deformation, creep, and fracture phenomena. Accordingly, this work demonstrates a generalizable approach for assessing the mechanical response of individual micron-sized molecular crystal particles and utilizing those responses in particle-level models. Graphic abstract

36 MATERIALS SCIENCE↗
DOE OSTI2020

Uncovering the Functional Role of Coherent Phonons during the Photoinduced Phase Transition in a Molecular Crystal

The atomic motions that make up phonons and molecular vibrations in molecular crystals influence their photophysical and electronic properties, including polaron formation, carrier mobility, and phase transitions. Discriminating between spectator and driving motions is a significant challenge hindering optimization. Unlocking this information and developing fine-tuned controls over actively participating phonon modes would not only lead to a stronger understanding of photochemistry but also provide a significant new tool in controlling solid state chemistry. We present a strategy using rationally designed double pulses to unveil the unique function of specific excited state phonon modes. Using ultrafast spectroscopy, we identified 50 and 90 cm –1 phonons involved in modulating the photoinduced spin-Peierls melting of potassium tetracyanoquinodimethane crystals. Here, we show that the 50 cm –1 phonon specifically corresponds to the coherent nuclear wavepacket involved in the charge transfer component of the overall spin-Peierls phase melting process, while the 90 cm –1 phonon facilitates the phase transition component.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Determination of Vibrational Motions Driving Photoinduced Electron Transfer Reactions in Molecular Crystals and Organic Thin Films

The objective of this research is to develop an understanding of how nuclear motions can be used to drive and control photo-induced charge transfer and singlet fission reactions in molecular crystals and thin films. By following the structural evolution of reacting molecules and nanostructures on the timescale of their nuclear motion, we will determine what interplay of nuclear coordinates is most efficient in driving electron transfers, a necessary first step in most photovoltaic and photocatalytic processes. We will identify which nuclear coordinates effectively promote desirable electron transfer processes, or ineffectively dissipate photon energy into unwanted thermal energy, thus guiding the design and implementation of molecular systems tailored for high efficiency processes. Specifically, we aim to identify which nuclear coordinates are responsible for driving electron transfer and singlet fission reactions at different points along the reaction coordinate, particularly in regimes which deviate from Marcus theory predictions. We will do so by investigating systems with electron transfer and singlet fission occurring from non-thermalized states due to strong donor-acceptor electronic coupling. First, we will prove that specific vibrational modes are responsible for charge transfer processes in molecular crystals. Secondly, we will investigate how vibrational coherences and specific nuclear motions can be used to drive high yield singlet fission processes. Finally, we will determine how complex film morphology affects long-range charge transport, a key limiting factor in our ability to make gains in efficiency.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Phonon-Induced Localization of Excitons in Molecular Crystals from First Principles

The spatial extent of excitons in molecular systems underpins their photophysics and utility for optoelectronic applications. Phonons are reported to lead to both exciton localization and delocalization. However, a microscopic understanding of phonon-induced (de)localization is lacking, in particular, how localized states form, the role of specific vibrations, and the relative importance of quantum and thermal nuclear fluctuations. Here, we present a first-principles study of these phenomena in solid pentacene, a prototypical molecular crystal, capturing the formation of bound excitons, exciton-phonon coupling to all orders, and phonon anharmonicity, using density functional theory, the ab initio GW-Bethe-Salpeter equation approach, finite-difference, and path integral techniques. In this study, we find that for pentacene zero-point nuclear motion causes uniformly strong localization, with thermal motion providing additional localization only for Wannier-Mott-like excitons. Anharmonic effects drive temperature-dependent localization, and, while such effects prevent the emergence of highly delocalized excitons, we explore the conditions under which these might be realized

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Spectra-Orthogonal Optical Anisotropy in Wafer-Scale Molecular Crystal Monolayers

Controlling the spectral and polarization response of two-dimensional (2D) crystals is vital for developing ultrathin platforms for compact optoelectronic devices. However, independently tuning optical anisotropy and spectral response remains challenging in conventional semiconductors due to the intertwined nature of their lattice and electronic structures. Here, we report spectra-orthogonal optical anisotropy─where polarization anisotropy is tuned independently of spectral response─in wafer-scale, one-atom-thick 2D molecular crystal (2DMC) monolayers synthesized on monolayer transition-metal dichalcogenide (TMD) crystals. Utilizing the concomitant spectral consistency and structural tunability of perylene derivatives, we demonstrate tunable optical polarization anisotropy in 2DMCs with similar spectral profiles, as confirmed by room-temperature scanning tunneling microscopy and cross-polarized reflectance microscopy. Additional angle-dependent analysis of the single-crystal and polycrystalline molecular domains reveals an epitaxial relationship between the 2DMC and TMD. In conclusion, our results establish a scalable, molecule-based 2D crystalline platform for unique and tunable functionalities unattainable in covalent 2D solids.

2D materials↗
DOE OSTI2025

Second‐ and Third‐Order Elastic Constants of Inert and Energetic Molecular Crystals From Density Functional Theory

Complete tensors of the second- and third-order elastic constants of the organic molecular crystals acetaminophen, pentaerythritol tetranitrate (PETN), cyclotrimethylene trinitramine (RDX), cyclotetramethylene tetranitramine (HMX), 1,1-diamino-2,2-dinitroethylene (FOX-7), hexanitrohexaazaisowurtzitane (CL-20), and erythritol tetranitrate (ETN) have been calculated using dispersion-corrected density functional theory. The sets of second- and third-order elastic constants are expected to provide a more accurate and reliable description of the behavior of these materials under nonhydrostatic loads than pressure- and volume-dependent second-order elastic constants. The tensors of second-order constants have been compared with experimental data and/or other calculations when possible, and with the exception of results for CL-20 from Brillouin scattering experiments, we find good agreement. The calculated third-order elastic constants of PETN are in very good agreement with the subset of third-order constants derived from experimental wave speed measurements. The elastic anisotropies of the crystals have been estimated using the universal elastic anisotropy index, which shows that the crystals fall into three groups with low anisotropy (PETN, RDX, and CL-20), moderate anisotropy (acetaminophen, HMX, and ETN), and high elastic anisotropy (FOX-7).

36 MATERIALS SCIENCE↗