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 91 records · Page 5

r 2 SCAN-D4: Dispersion corrected meta-generalized gradient approximation for general chemical applications

We combine a regularized variant of the strongly constrained and appropriately normed semilocal density functional with the latest generation semi-classical London dispersion correction. The resulting density functional approximation r 2 SCAN-D4 has the speed of generalized gradient approximations while approaching the accuracy of hybrid functionals for general chemical applications. We demonstrate its numerical robustness in real-life settings and benchmark molecular geometries, general main group and organo-metallic thermochemistry, and non-covalent interactions in supramolecular complexes and molecular crystals. Main group and transition metal bond lengths have errors of just 0.8%, which is competitive with hybrid functionals for main group molecules and outperforms them for transition metal complexes. The weighted mean absolute deviation (WTMAD2) on the large GMTKN55 database of chemical properties is exceptionally small at 7.5 kcal/mol. This also holds for metal organic reactions with an MAD of 3.3 kcal/mol. The versatile applicability to organic and metal–organic systems transfers to condensed systems, where lattice energies of molecular crystals are within the chemical accuracy (errors <1 kcal/mol).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Subsurface spectroscopy of heterogeneous materials using optical wavefront shaping

Plastic-bonded explosives and propellants consist of energetic molecular crystals embedded in a polymer matrix that also includes additives, such as taggants, plasticizers, grit, antioxidants, etc. Furthermore, this heterogeneous composition renders these materials optically opaque, limiting optical characterization techniques (e.g., Raman spectroscopy) to probing chemical reactions at the sur-face. However, many reactions of interest are believed to occur inside the material where current optical techniques can not probe. To address this challenge, we have developed a new optical technique that utilizes feedback assisted wavefront shaping (using spatial light modulators) to focus probe light inside a heterogeneous material, such that we can perform subsurface spectroscopy. Recently, we demonstrated this technique by monitoring subsurface photodegradation and thermal degradation of Eu-doped molecular crystals using both fluorescence and Raman spectroscopy. Based off these successes we are now looking into improvements to this technique to allow for greater time resolution, including replacing our spatial light modulators with fast digital mirror devices, using ns optical phase conjugation, and self-healing optical beams.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Calculation and interpretation of classical turning surfaces in solids

Abstract Classical turning surfaces of Kohn–Sham potentials separate classically allowed regions (CARs) from classically forbidden regions (CFRs). They are useful for understanding many chemical properties of molecules but need not exist in solids, where the density never decays to zero. At equilibrium geometries, we find that CFRs are absent in perfect metals, rare in covalent semiconductors at equilibrium, but common in ionic and molecular crystals. In all materials, CFRs appear or grow as the internuclear distances are uniformly expanded. They can also appear at a monovacancy in a metal. Calculations with several approximate density functionals and codes confirm these behaviors. A classical picture of conduction suggests that CARs should be connected in metals, and disconnected in wide-gap insulators, and is confirmed in the limits of extreme compression and expansion. Surprisingly, many semiconductors have no CFR at equilibrium, a key finding for density functional construction. Nonetheless, a strong correlation with insulating behavior can still be inferred. Moreover, equilibrium bond lengths for all cases can be estimated from the bond type and the sum of the classical turning radii of the free atoms or ions.

Chemistry↗

Energy conversion and storage via photoinduced polarization change in non-ferroelectric molecular [CoGa] crystals

To alleviate the energy and environmental crisis, in the last decades, energy harvesting by utilizing optical control has emerged as a promising solution. Here we report a polar crystal that exhibits photoenergy conversion and energy storage upon light irradiation. The polar crystal consists of dinuclear [CoGa] molecules, which are oriented in a uniform direction inside the crystal lattice. Irradiation with green light induces a directional intramolecular electron transfer from the ligand to a low-spin Co III centre, and the resultant light-induced high-spin Co II excited state is trapped at low temperature, realizing energy storage. Additionally, electric current release is observed during relaxation from the trapped light-induced metastable state to the ground state, because the intramolecular electron transfer in the relaxation process is accompanied with macroscopic polarization switching at the single-crystal level. It demonstrates that energy storage and conversion to electrical energy is realized in the [CoGa] crystals, which is different from typical polar pyroelectric compounds that exhibit the conversion of thermal energy into electricity.

25 ENERGY STORAGE↗

End-to-end AI framework for interpretable prediction of molecular and crystal properties

We introduce an end-to-end computational framework that allows for hyperparameter optimization using the DeepHyper library, accelerated model training, and interpretable AI inference. The framework is based on state-of-the-art AI models including CGCNN, PhysNet, SchNet, MPNN, MPNN-transformer, and TorchMD-NET. We employ these AI models along with the benchmark QM9, hMOF, and MD17 datasets to showcase how the models can predict user-specified material properties within modern computing environments. We demonstrate transferable applications in the modeling of small molecules, inorganic crystals and nanoporous metal organic frameworks with a unified, standalone framework. We have deployed and tested this framework in the ThetaGPU supercomputer at the Argonne Leadership Computing Facility, and in the Delta supercomputer at the National Center for Supercomputing Applications to provide researchers with modern tools to conduct accelerated AI-driven discovery in leadership-class computing environments. We release these digital assets as open source scientific software in GitLab, and ready-to-use Jupyter notebooks in Google Colab.

36 MATERIALS SCIENCE↗

Crystal microbalance measures condensable molecular fluxes

Quartz crystal quantitatively measures molecular fluxes emanating from and condensing on spacecraft surfaces. Vibrating in a thickness shear mode, the crystal is frequency sensitive to changes in mass on its surface and can measure a fractional monolayer of a condensate.

Stephens, J. B.↗

Effect of melt shearing on D -mannitol crystal twisting in the presence of small molecule and macromolecular additives

Here, many molecular crystals grow as twisted lamellae from the melt when the driving force is high, but the mechanisms governing the spontaneous formation of helicoidal crystallites are varied and those operative in particular cases are difficult to identify. Polyvinylpyrrolidone (PVP) induces the crystallization of D-mannitol as banded spherulites comprised of twisted fibrils emanating radially from spherulite nucleation centers. The molecular weight of the PVP phase affects the twisting pitch, P, (rotation of the crystallite by 180°), in D-mannitol thin films crystallized from the melt. At a relatively high crystallization temperature (T c =130°C), P was sensitive to the PVP molecular weight, ranging from 390±20 μm for D-mannitol films incorporating 15 wt% 10 kDa PVP to 20±3 μm for the same weight of 1300 kDa PVP. Magnitudes of complex viscosities, |η*|, of D-mannitol/PVP melts measured via small-amplitude oscillatory shearing were strongly dependent on the PVP molecular weight, but were not correlated to P(T). Instead, P was sensitive to the dynamic PVP chain conformation during D-mannitol crystallization. Under steady torsional shear, P decreased from~30 μm to~8μm with increasing shear rates from 0.01 to 100 s -1 for D-mannitol films crystallized at 130°C in the presence of 15 wt% 10 or 1300 kDa PVP. Shear forces decrease the entanglement density of polymer chains while orienting the chains along the stream lines of the viscometric flow, indicating that the conformations of macromolecular additives can affect the pitch of banded spherulites. By contrast, P in D-mannitol twisted by D-sorbitol was independent of shear rate.

36 MATERIALS SCIENCE↗

Suppression of Disorder in Benzamide and Thiobenzamide Crystals by Fluorine Substitution

Disorder is a common feature of molecular crystals that complicates determination of structures and can potentially affect electric and mechanical properties. Suppression of disorder is observed in otherwise severely disordered benzamide and thiobenzamide crystals by substituting hydrogen with fluorine in the ortho-position of the phenyl ring. Fluorine occupancies of 20–30% are sufficient to suppress disorder without changing the packing motif. Crystal structure prediction calculations reveal a much denser lattice energy landscape for benzamide compared to 2-fluorobenzamide, suggesting that fluorine substitution makes disorder less likely.

36 MATERIALS SCIENCE↗

On-lattice kinetic Monte Carlo approaches for modeling molecular anisotropy in resveratrol crystallization

Stilbenes are a class of organic compounds with broad-ranging pharmaceutical and agricultural applications, which are typically isolated and purified through recrystallization. We are motivated by reducing experimental waste and optimizing yield via developing predictive simulations for processing-dependent crystal morphologies. Using resveratrol as a model stilbene system, we have developed an approach for simulating crystallization with molecular resolution using on-lattice kinetic Monte Carlo. In this work, we highlight modifications to the Stochastic Parallel PARticle Kinetic Simulator (SPPARKS) software package, which were essential to this application. Key enhancements include the incorporation of non-orthogonal cell shapes and monomer anisotropy approximations using bound hard spheres. This new SPPARKS application has been applied to resveratrol with attachment energy libraries obtained from density functional theory, resulting in excellent agreement with experimental morphology prediction.

crystallization↗

The mechanism driving a solid–solid phase transition in a biomacromolecular crystal

A solid-solid phase transition (SSPT) occurs between distinguishable crystalline forms. Because of its importance in application and theory in material science and condensed matter physics, SSPT has been studied most extensively in metallic alloys, inorganic salt or small organic molecular crystals, but much less so in biomacromolecular crystals. In general, the mechanism of SSPT at atomic and molecular levels is not well understood. Here, we describe the ordered molecular rearrangements in biomacromolecular crystals of the adenine riboswitch (riboA) aptamer using real-time serial crystallography and solution atomic force microscopy (AFM). The large, ligand-induced conformational changes drive the initial phase transition from the apo unit cell (AUC) to the trans unit cell 1 (TUC1). During this transition, coaxial stacking of P1 duplexes becomes the dominant packing interface, whereas P2-P2 interactions are almost completely disrupted, resulting in “floating” layers of molecules. The coupling points in TUC1 and their local conformational flexibility allow the molecules to reorganize to achieve the more densely packed and energetically favorable bound unit cell (BUC). Our study thus reveals the interplay between the conformational changes and the crystal phases—the underlying mechanism that drives the phase transition. Using polarized video microscopy (PVM) to monitor the SSPT in small crystals at high ligand concentration, we have identified the time window during which the major conformational changes take place, and simulated the in crystallo kinetics. Together, these results provide the spatiotemporal information necessary for informing time-resolved crystallography (TRX) experiments. Moreover, this study illustrates a practical approach to characterize SSPT in transparent crystals.

36 MATERIALS SCIENCE↗

Protein crystal growth in microgravity: Temperature induced large scale crystallization of insulin

One of the major stumbling blocks that prevents rapid structure determination using x-ray crystallography is macro-molecular crystal growth. There are many examples where crystallization takes longer than structure determination. In some cases, it is impossible to grow useful crystals on earth. Recent experiments conducted in conjuction with NASA on various Space Shuttle missions have demonstrated that protein crystals often grow larger and display better internal molecular order than their earth-grown counterparts. This paper reports results from three Shuttle flights using the Protein Crystallization Facility (PCF). The PCF hardware produced large, high-quality insulin crystals by using a temperature change as the sole means to affect protein solubility and thus, crystallization. The facility consists of cylinders/containers with volumes of 500, 200, 100, and 50 ml. Data from the three Shuttle flights demonstrated that larger, higher resolution crystals (as evidenced by x-ray diffraction data) were obtained from the microgravity experiments when compared to earth-grown crystals.

Long, Marianna M.↗

The Potential Energy Hotspot: Effects of Impact Velocity, Defect Geometry, and Crystallographic Orientation

In energetic materials, the localization of energy into “hotspots” is known to dictate the initiation of chemical reactions and detonation. Recent all-atom simulations have shown that more energy is localized as internal potential energy (PE) than can be inferred from the kinetic energy (KE) alone. The mechanisms associated with pore collapse and hotspot formation are known to depend on pore geometry and dynamic material response such as plasticity. Therefore, we use molecular dynamics (MD) simulations to characterize shock-induced pore collapse and the subsequent formation of hotspots in 1,3,5-triamino-2,4,6-trinitrobenzene (TATB), a highly anisotropic molecular crystal, for various defect shapes, shock strengths, and crystallographic orientations. We find that the localization of energy as PE is consistently larger than the KE in cases with significant plastic deformation. Additionally, an analysis of MD trajectories reveals the underlying molecular- and crystal-level processes that govern the effect of orientation and pore shape on PE localization. We find that the regions of highest PE relate to the areas of maximum plastic deformation, while KE is maximized at the point of impact. Comparisons against octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) reveal less energy localization in TATB, which could be a contributing factor to the latter’s insensitivity.

36 MATERIALS SCIENCE↗

Intermediate polaronic charge transport in organic crystals from a many-body first-principles approach

Abstract Charge transport in organic molecular crystals (OMCs) is conventionally categorized into two limiting regimes − band transport, characterized by weak electron-phonon (e-ph) interactions, and charge hopping due to localized polarons formed by strong e-ph interactions. However, between these two limiting cases there is a less well understood intermediate regime where polarons are present but transport does not occur via hopping. Here we show a many-body first-principles approach that can accurately predict the carrier mobility in this intermediate regime and shed light on its microscopic origin. Our approach combines a finite-temperature cumulant method to describe strong e-ph interactions with Green-Kubo transport calculations. We apply this parameter-free framework to naphthalene crystal, demonstrating electron mobility predictions within a factor of 1.5−2 of experiment between 100 and 300 K. Our analysis reveals the formation of a broad polaron satellite peak in the electron spectral function and the failure of the Boltzmann equation in the intermediate regime.

36 MATERIALS SCIENCE↗