Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular properties”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Chemical design of electronic and magnetic energy scales of tetravalent praseodymium materials

Lanthanides in the trivalent oxidation state are typically described using an ionic picture that leads to localized magnetic moments. The hierarchical energy scales associated with trivalent lanthanides produce desirable properties for e.g., molecular magnetism, quantum materials, and quantum transduction. Here, we show that this traditional ionic paradigm breaks down for praseodymium in the tetravalent oxidation state. Synthetic, spectroscopic, and theoretical tools deployed on several solid-state Pr 4+ -oxides uncover the unusual participation of 4f orbitals in bonding and the anomalous hybridization of the 4f 1 configuration with ligand valence electrons, analogous to transition metals. The competition between crystal-field and spin-orbit-coupling interactions fundamentally transforms the spin-orbital magnetism of Pr 4+ , which departs from the J eff = 1/2 limit and resembles that of high-valent actinides. Our results show that Pr 4+ ions are in a class on their own, where the hierarchy of single-ion energy scales can be tailored to explore new correlated phenomena in quantum materials.

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

DNA double-strand break movement in heterochromatin depends on the histone acetyltransferase dGcn5

Abstract Cells employ diverse strategies to repair double-strand breaks (DSBs), a dangerous form of DNA damage that threatens genome integrity. Eukaryotic nuclei consist of different chromatin environments, each displaying distinct molecular and biophysical properties that can significantly influence the DSB-repair process. DSBs arising in the compact and silenced heterochromatin domains have been found to move to the heterochromatin periphery in mouse and Drosophila to prevent aberrant recombination events. However, it is poorly understood how chromatin components, such as histone post-translational modifications, contribute to these DSB movements within heterochromatin. Using irradiation as well as locus-specific DSB induction in Drosophila tissues and cultured cells, we find enrichment of histone H3 lysine 9 acetylation (H3K9ac) at DSBs in heterochromatin but not euchromatin. We find this increase is mediated by the histone acetyltransferase dGcn5, which rapidly localizes to heterochromatic DSBs. Moreover, we demonstrate that in the absence of dGcn5, heterochromatic DSBs display impaired recruitment of the SUMO E3 ligase Nse2/Qjt and fail to relocate to the heterochromatin periphery to complete repair. In summary, our results reveal a previously unidentified role for dGcn5 and H3K9ac in heterochromatic DSB repair and underscore the importance of differential chromatin responses at heterochromatic and euchromatic DSBs to promote safe repair.

Biochemistry & Molecular Biology↗
DOE OSTI2023

Fractionated and purified hybrid poplar lignins as a polyol replacement in rigid polyurethane/polyisocyanurate foams

This study introduces fractionated lignin as an innovative component in the formulation of rigid polyurethane/polyisocyanurate (PUR/PIR) foams. Low-density PUR/PIR rigid foams were prepared by replacing 80% of the petrochemical-based polyol with a hybrid poplar (HP) lignin, recovered via alkaline pretreatment, and also with a fraction of this lignin, isolated via the Aqueous Lignin Purification with Hot Agents (ALPHA) process. The as-recovered HP and ALPHA-fractionated HP lignins were characterized to determine hydroxyl content, molecular weight (Mw) distribution, and pH. Both lignin-based foams met minimum standard requirements in terms of closed cell content, compression strength, and thermal conductivity. Notably, the foams made with ALPHA-based lignin outperformed the as-recovered lignin-based foams in all measured foam properties. Moreover, the ALPHA-based foam had comparable performance to the control foam (without lignin) except for density and even surpassed the control foam in closed cell content and compressive strength. This investigation of the molecular properties of lignin suggests that significant reductions in lignin average molecular weight and polydispersity (PDI) can positively impact the properties of lignin-based rigid foams.

09 BIOMASS FUELS↗
DOE OSTI2023

Properties of carbon up to 10 million kelvin from Kohn-Sham density functional theory molecular dynamics

Accurately modeling dense plasmas over wide-ranging conditions of pressure and temperature is a grand challenge critically important to our understanding of stellar and planetary physics as well as inertial confinement fusion. In this work, we employ Kohn-Sham density functional theory (DFT) molecular dynamics (MD) to compute the properties of carbon at warm and hot dense matter conditions in the vicinity of the principal Hugoniot. In particular, we calculate the equation of state (EOS), Hugoniot, pair distribution functions, and diffusion coefficients for carbon at densities spanning 8 g/$\mathrm{cm^3}$ to 16 g/$\mathrm{cm^3}$ and temperatures ranging from 100 kK to 10 MK using the Spectral Quadrature method. Here, we find that the computed EOS and Hugoniot are in good agreement with path integral Monte Carlo results and the sesame database. Additionally, we calculate the ion-ion structure factor and viscosity for selected points. All results presented are at the level of full Kohn-Sham DFT-MD, free of empirical parameters, average-atom, and orbital-free approximations employed previously at such conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Toward Addressing the Challenge to Predict the Heat Capacities of RDX and HMX Energetic Materials

Availability of heat capacity as function of pressure and temperature is an essential prerequisite for development of a computational multiscale strategy capable to address the evolution of microstructure and energy release in advanced high energy density materials. In the case of 1,3,5-trinitro-1,3,5-triazinane (RDX) and 1,3,5,7-tetranitro-1,3,5,7-tetrazocane (HMX) systems as two of the most studied energetic materials, there are substantial gaps in experimental data, with available heat capacities values distributed only in a region close to standard ambient conditions. In this study we demonstrate how these major experimental limitations can be addressed in the case of the RDX and HMX systems based on the combined use of classical and quantum mechanical calculations. We show that by considering ideal gas properties evaluated using quantum mechanical methods, and residual properties obtained from molecular simulations using fully flexible atomistic force field models, excellent agreement can be obtained for the predicted heat capacities to the most recent experimental values. An important advantage of the current computational methodology is that it allows evaluation of both constant-volume and constant-pressure heat capacities for a broad interval of temperatures and pressures, which encompasses solid and liquid phases conditions. In the case of the solid α and γ phases of RDX and the β phase of HMX, the predicted results follow closely both the available experimental data at standard ambient conditions and the results obtained using density functional theory calculations at high pressures, a regime where experimental data are not available. A perspective to expand the current methodology is also discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Self-seeded growth of very large open-structured zeolite nanosheet assemblies with extraordinary micropore accessibility

Zeolite nanosheets (ZNs) offer improved micropore accessibilities and transport properties for enhanced molecular catalysis and separations. However, practical application of the ZN materials is hampered by the lack of efficient synthesis methods. Here, in this work, a ZN self-seeded method is demonstrated for single-step reproduction of flower-like assemblies of very large MFI ZN plates. The ZN plates are ~60 nm-thick stacks of 4-nm-thick single-crystal sheets. The ZN flower growth involves terrace nucleation on the seed surfaces and subsequent ZN epitaxial growth in [010] orientation directed by a diquaternary agent. The open architecture of the assemblies prevented collapse and agglomeration of ZNs during thermal activation that effectively preserved the interconnected intra-sheet and inter-sheet micropore system. Thus, the ZN assemblies exhibited markedly enhanced molecular adsorption capacity and transport diffusivity for the probing xylene molecules as compared to the conventional crystals. The ZN assembly and its harvestable very large-sized ZNs have the potential for developing high-performance ZN adsorbents, catalysts, and molecular-sieve membranes.

36 MATERIALS SCIENCE↗
DOE OSTI2025

The total neutron cross section of liquid and solid ammonia

Ammonia is a material of interest for future neutron moderators at high-power sources due to its high hydrogen density, low melting point, and resistance to polymerization in an intense radiation field. Its performance in such applications cannot currently be calculated due to the absence of suitable computer models for the interaction of neutrons with ammonia under relevant conditions. In an effort to develop suitable scattering kernels for computer simulations of moderator performance, we have conducted a series of Density Functional Theory and Molecular Dynamics calculations of the molecular-level thermal properties of ammonia at various temperatures within both the solid and liquid phases. In this paper, we compare computer calculations for the energy-dependent total neutron cross section of ammonia, based on these models, to experimental measurements of those cross sections at temperatures of 221 K, 180 K, and 35 K. The experimental data were collected over an energy range from 0.1 meV to 10 eV using time-of-flight techniques at the Low Energy Neutron Source (LENS) facility at Indiana University. This comparison provides a first validation in the development of thermal scattering libraries for Monte Carlo source design simulations based on liquid and solid ammonia. In conclusion, we also provide some insights into where additional development of tools for creating such models may be needed.

Ammonia↗
DOE OSTI2025

Understanding Differences in Water Adsorption Isotherms: Structural Variations, Force Fields, and Monte Carlo Simulation Approaches

Accurate prediction of water adsorption in micro- and mesoporous materials with hydrophobic pores is essential for the design and characterization of advanced adsorbent materials for separation and energy applications. Here, we assess the reproducibility and consistency of water adsorption isotherms in two microporous all-silica MFI zeolite structures (MFI-K and MFI-O) using two different zeolite force fields and three simulation approaches: grand canonical Monte Carlo (GCMC), Gibbs ensemble Monte Carlo (GEMC), and transition matrix Monte Carlo (TMMC). We demonstrate that consistent treatment of the bulk fluid phase in GCMC and TMMC simulations is critical for reconciling isotherms across methods, and we construct simulation-based equations of state for the TIP4P water model to enable rigorous fugacity-to-pressure conversions. Large shifts in the isotherms are observed for two zeolite force fields developed using different parametrization strategies, with the GCS force field representing implicitly a defect-containing all-silica zeolite, whereas the TraPPE-zeo force field accurately represents an essentially defect-free all-silica zeolite. While water in the van Koningsveld structure of MFI exhibits a first-order phase transition and condensation-like step for adsorption near room temperature, water in the Olson structure of MFI displays continuous adsorption, attributed to differences in the adsorption free energy landscapes. Structural analysis reveals that small geometric variations, particularly Si–O–Si bond angles near the strongest adsorption sites, lead to these substantial differences in adsorption behavior. Furthermore, our results highlight the sensitivity of simulated water adsorption isotherms in hydrophobic frameworks to seemingly small differences in the framework structures, force field parametrization, and simulation approaches.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A Meta-Generalized Gradient Approximation for the Cavity-Dependent Exchange-Correlation Interaction in Strongly Coupled Light–Matter Systems

Strong light–matter coupling in optical cavities enables the manipulation of chemical and physical properties without altering molecular composition. Theoretical modeling of such phenomena requires exchange-correlation (XC) functionals that account for both electron–electron and electron–photon (ep) interactions within quantum electrodynamical density functional theory (QEDFT). In this work, we develop a meta-generalized gradient approximation (meta-GGA) specifically targeting the cavity-dependent XC interaction in strongly coupled light–matter systems. This novel approximation is built upon a new semilocal polarizability approximation, which draws from the jellium-with-a-gap model, and can be extended to a “global hybrid” variant that goes beyond the isotropic model from previous approximations. The polarizability model yields significantly improved dispersion coefficients and benchmark calculations with the cavity-dependent XC functional demonstrate improved agreement with QED Hartree–Fock (QED-HF) reference energies. Application to the regioselectivity of brominated nitrobenzene intermediates reveals the functional’s capacity to capture cavity-induced energetic shifts. In conclusion, our results advance the Jacob’s ladder of functionals for QEDFT and provide a practical tool for modeling polaritonic chemistry.

Approximation↗
DOE OSTI2024

Emergence of Orchestrated and Dynamic Metabolism of Saccharomyces cerevisiae

Microbial metabolism is a fundamental cellular process that involves many biochemical events and is distinguished by its emergent properties. While the molecular details of individual reactions have been increasingly elucidated, it is not well understood how these reactions are quantitatively orchestrated to produce collective cellular behaviors. Here we developed a coarse-grained, systems, and dynamic mathematical framework, which integrates metabolic reactions with signal transduction and gene regulation to dissect the emergent metabolic traits of Saccharomyces cerevisiae. Our framework mechanistically captures a set of characteristic cellular behaviors, including the Crabtree effect, diauxic shift, diauxic lag time, and differential growth under nutrient-altered environments. It also allows modular expansion for zooming in on specific pathways for detailed metabolic profiles. Furthermore, this study provides a systems mathematical framework for yeast metabolic behaviors, providing insights into yeast physiology and metabolic engineering.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Quantum chemical calculations of lithium-ion battery electrolyte and interphase species

Abstract Lithium-ion batteries (LIBs) represent the state of the art in high-density energy storage. To further advance LIB technology, a fundamental understanding of the underlying chemical processes is required. In particular, the decomposition of electrolyte species and associated formation of the solid electrolyte interphase (SEI) is critical for LIB performance. However, SEI formation is poorly understood, in part due to insufficient exploration of the vast reactive space. The Lithium-Ion Battery Electrolyte (LIBE) dataset reported here aims to provide accurate first-principles data to improve the understanding of SEI species and associated reactions. The dataset was generated by fragmenting a set of principal molecules, including solvents, salts, and SEI products, and then selectively recombining a subset of the fragments. All candidate molecules were analyzed at the ω B97X-V/def2-TZVPPD/SMD level of theory at various charges and spin multiplicities. In total, LIBE contains structural, thermodynamic, and vibrational information on over 17,000 unique species. In addition to studies of reactivity in LIBs, this dataset may prove useful for machine learning of molecular and reaction properties.

25 ENERGY STORAGE↗
DOE OSTI2022

Ethylene polymerization with a crystallographically well-defined metal–organic framework supported catalyst

The inherent crystallinity of metal–organic framework (MOF) catalysts offers the possibility to understand the structure of the active site at the molecular level. This property is often lacking in traditional amorphous supports. Cr-SIM-NU-1000, a MOF-supported Cr 3+ heterogeneous catalyst, is shown to be competent for ethylene polymerization after activation with AlEt 2 Cl (DEAC), producing crystalline linear polyethylene (PE). The polymer produced has a low polydispersity (D = 2.0), in marked contrast to the Phillips supported chromium catalyst, Cr@SiO 2 (D ~ 8–65). Cr-SIM-NU-1000 achieves a turnover frequency of 2.6 × 10 3 h –1 under 40 bar ethylene pressure at room temperature, with corresponding PE productivity of 1.3 × 10 5 g PE mol –1 Cr per h. Single crystal X-ray diffraction (SC-XRD) of the pre-catalyst was conducted by collecting a structure of alkyl aluminum (DEAC) co-catalyst treated Cr-SIM-NU-1000. This crystal structure provides insight into the interactions between DEAC co-catalyst and Cr active site, revealing a Cr–C bond after treatment with DEAC. Furthermore, DEAC is crystallographically resolved at the terminal oxy-ligands of the node and likely also exists within the window pores of the framework between nodes, based on electron density mapping. Cr-SIM-NU-1000 offers the opportunity to study a structurally well-defined olefin polymerization system, with atomically precise characterization of the pre-catalyst structure. Furthermore, this allows the proposal of a mechanism and feeds into future development of next-generation heterogeneous catalyst systems.

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

Highly-destabilized ligand field excited states of iron carbene complexes and their relation to charge transfer state lifetimes

Lifetimes of photoexcited charge transfer (CT) states in transition metal chromophores are influenced by low-lying ligand field (LF) excited states, especially for 3d metal complexes. To manipulate interactions between LF and CT states, it is important to be able to control LF excited state energies using tunable synthetic variables. In this report, we use Fe 2p3d L 3 -edge resonant inelastic X-ray scattering (RIXS) to measure LF excited state energies of three homoleptic iron chromophores coordinated by strong-field N-heterocyclic carbenes (NHCs). We investigate the effect of oxidation state and ligand scaffold on LF energies and covalency parameters. A cyclometalated bis(NHC) ligand affords both high LF excited state energies (and thus high 10 Dq) as well as high metal–ligand covalency compared to other iron complexes with very strong-field ligands. However, for the set of complexes investigated, we do not observe meaningful correlation between the LF excited state energies and the CT excited state lifetimes. These results illustrate that targeting long-lived CT excited states necessitates control of multiple molecular excited state properties, with destabilization of the LF excited state energies proving necessary, but insufficient, to control the CT excited state lifetime in Fe carbene complexes.

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

Resistive switching properties of monolayer h-BN atomristors with different electrodes

Resistive switching properties based on molecular beam epitaxy-grown monolayer hexagonal boron nitride (h-BN) atomristors are studied by using metal insulator metal configurations with different electrode materials. Au/monolayer h-BN/Ni devices demonstrate a forming-free bipolar resistive switching (BRS) behavior, a good endurance with up to 97 cycles at a high compliance current of 100 mA, an average on/off ratio of 10 3 , and a low set/reset voltage variability. Metal/monolayer h-BN/graphite/Co devices exhibit self-compliant current BRS characteristics. Both metal/h-BN/Ni and metal/h-BN/graphite/Co devices show the coexistence of BRS, unipolar resistive switching (URS), and nonvolatile threshold switching (TH) modes. The formation of conductive filaments is attributed to the diffusion and trapping of metal ions on the defect sites driven by the electric field, while the rupture is driven by the electric field in BRS and by Joule heating in URS and TH modes.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Fast uncertainty estimates in deep learning interatomic potentials

Deep learning has emerged as a promising paradigm to give access to highly accurate predictions of molecular and material properties. A common short-coming shared by current approaches, however, is that neural networks only give point estimates of their predictions and do not come with predictive uncertainties associated with these estimates. Existing uncertainty quantification efforts have primarily leveraged the standard deviation of predictions across an ensemble of independently trained neural networks. This incurs a large computational overhead in both training and prediction, resulting in order-of-magnitude more expensive predictions. Here, we propose a method to estimate the predictive uncertainty based on a single neural network without the need for an ensemble. This allows us to obtain uncertainty estimates with virtually no additional computational overhead over standard training and inference. We demonstrate that the quality of the uncertainty estimates matches those obtained from deep ensembles. We further examine the uncertainty estimates of our methods and deep ensembles across the configuration space of our test system and compare the uncertainties to the potential energy surface. Finally, we study the efficacy of the method in an active learning setting and find the results to match an ensemble-based strategy at order-of-magnitude reduced computational cost.

Chemistry↗
DOE OSTI2023

Roadmap on electronic structure codes in the exascale era

Electronic structure calculations have been instrumental in providing many important insights into a range of physical and chemical properties of various molecular and solid-state systems. Their importance to various fields, including materials science, chemical sciences, computational chemistry and device physics, is underscored by the large fraction of available public supercomputing resources devoted to these calculations. As we enter the exascale era, exciting new opportunities to increase simulation numbers, sizes, and accuracies present themselves. In order to realize these promises, the community of electronic structure software developers will however first have to tackle a number of challenges pertaining to the efficient use of new architectures that will rely heavily on massive parallelism and hardware accelerators. This roadmap provides a broad overview of the state-of-the-art in electronic structure calculations and of the various new directions being pursued by the community. It covers 14 electronic structure codes, presenting their current status, their development priorities over the next five years, and their plans towards tackling the challenges and leveraging the opportunities presented by the advent of exascale computing.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Analysis of diagonal G and subspace W approximations within fully self-consistent GW calculations for bulk semiconducting systems

Fully self-consistent GW (sc-GW) methods are now available to evaluate quasiparticle and spectral properties of various molecular and bulk systems. However, such techniques based on the full matrix of G and W are computationally demanding. Additionally, the routinely used single-shot GW approximation (G 0 W 0 ) has an undesirable dependency on the choice of initial exchange-correlation functional. In the literature, many so-called self-consistent GW methods are based on diagonal approximation of G and low-ranking approximation of W. It is thus worth checking how good such approximations are in comparison with the full matrix method. In this work, we consider AlAs, AlP, GaP, and ZnS as the prototype systems to perform sc-GW calculations by expressing the full G matrix using a plane-wave basis set. We compared our sc-GW results with the diagonal G and subspace W approximated sc-GW results (sc-GW-diagG and sc-GW-subW methods). In the sc-GW-diagG method, interacting G is expanded in the eigenvectors of noninteracting G such that only diagonal elements are retained, whereas the number of eigenmodes is truncated in sc-GW-subW calculations. A systematic analysis of the results obtained from the above techniques is presented. The differences in the quasiparticle band gap between the approximated and the full matrix sc-GW approaches are mostly less than 1.7%, which validates such widely adopted approximations, and also shows how such low-ranking approximation can be used to include higher-order terms such as the vertex correction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Atoms to fibers: Identifying novel processing methods in the synthesis of pitch-based carbon fibers

Understanding and optimizing the key mechanisms used in the synthesis of pitch-based carbon fibers (CFs) are challenging, because unlike polyacrylonitrile-based CFs, the feedstock for pitch-based CFs is chemically hetero- geneous, resulting in complex fabrication leading to inconsistency in the final properties. In this work, we use molecular dynamics simulations to explore the processing and chemical phase space through a framework of CF models to identify their effects on elastic performance. The results are in excellent agreement with experiments. We find that density, followed by alignment, and functionality of the molecular constituents dictate the CF mechanical properties more strongly than their size and shape. Last, we propose a previously unexplored fabrication route for high-modulus CFs. Unlike graphitization, this results in increased sp 3 fraction, achieved via generating high-density CFs. In addition, the high sp 3 fraction leads to the fabrication of CFs with isometric compressive and tensile moduli, enabling their potential applications for compressive loading.

01 COAL, LIGNITE, AND PEAT↗