Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular geometry”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Excited-state structural dynamics of nickel complexes probed by optical and X-ray transient absorption spectroscopies: insights and implications

Excited states of nickel complexes undergo a variety of photochemical processes, such as charge transfer, ligation/deligation, and redox reactions, relevant to solar energy conversion and photocatalysis. The efficiencies of the aforementioned processes are closely coupled to the molecular structures in the ground and excited states. The conventional optical transient absorption spectroscopy has revealed important excited-state pathways and kinetics, but information regarding the metal center, in particular transient structural and electronic properties, remains limited. These deficiencies are addressed by X-ray transient absorption (XTA) spectroscopy, a detailed probe of 3d orbital occupancy, oxidation state and coordination geometry. The examples of excited-state structural dynamics of nickel porphyrin and nickel phthalocyanine have been described from our previous studies with highlights on the unique structural information obtained by XTA spectroscopy. Furthermore, we close by surveying prospective applications of XTA spectroscopy to active areas of Ni-based photocatalysis based on the knowledge gained from our previous studies.

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

Comparing the structures and photophysical properties of two charge transfer co-crystals

Organic co-crystals have emerged as a promising class of semiconductors for next-generation optoelectronic devices due to their unique photophysical properties. This paper presents a joint experimental-theoretical study comparing the crystal structure, spectroscopy, and electronic structure of two charge transfer co-crystals. Reported herein is a novel co-crystal Npe:TCNQ, formed from 4-(1-naphthylvinyl)pyridine (Npe) and 7,7,8,8-tetracyanoquinodimethane (TCNQ) via molecular self-assembly. This work also presents a revised study of the co-crystal composed of Npe and 1,2,4,5-tetracyanobenzene (TCNB) molecules, Npe:TCNB, herein reported with a higher-symmetry (monoclinic) crystal structure than previously published. Npe:TCNB and Npe:TCNQ dimer clusters are used as theoretical model systems for the co-crystals; the geometries of the dimers are compared to geometries of the extended solids, which are computed with periodic boundary conditions density functional theory. UV-Vis absorption spectra of the dimers are computed with time-dependent density functional theory and compared to experimental UV-Vis diffuse reflectance spectra. Both Npe:TCNB and Npe:TCNQ are found to exhibit neutral character in the S 0 state and ionic character in the S 1 state. In conclusion, the high degree of charge transfer in the S 1 state of both Npe:TCNB and Npe:TCNQ is rationalized by analyzing the changes in orbital localization associated with the S 1 transitions.

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

Molecular concentration field design using closed-form steady-state solutions

Control over spatial concentration fields represents a fundamental challenge in designing synthetic biological systems and programmable soft materials. While nature creates morphogen gradients that orchestrate complex developmental processes, synthetic approaches have largely relied on empirical optimization and computationally intensive simulations. Here, we present an analytical framework for steady-state concentration fields generated by finite-sized localized sources in diffusion–degradation systems and derive closed-form solutions for one-, two-, and three-dimensional geometries. By expressing these solutions in dimensionless form, we show that gradient steepness and spatial structure are organized by the Thiele modulus, which captures the competition between diffusion and degradation length scales. The analysis reveals distinct design regimes: in degradation-dominated systems, gradient shape is governed by exponential decay and becomes dimension-independent, whereas in diffusion-dominated systems, gradient magnitude and extent follow dimension-dependent power-law scaling. Building on these results, we introduce a quantitative design strategy that uses threshold-based criteria to program concentration ranges by tuning physically accessible parameters, most directly the production rate, while holding transport and degradation properties fixed. Comparisons with numerical solutions and reported experimental systems demonstrate consistency with the predicted scaling behavior. Together, this work provides a generalizable and physically transparent framework for designing steady-state concentration fields in synthetic biological and soft matter systems, enabling predictive control of gradient-mediated organization without reliance on extensive numerical optimization.

Kim, Dong Woo [Johns Hopkins University, Baltimore↗
DOE OSTI2024

High-throughput dataset of impurity adsorption on common catalysts in biomass upgrading applications

Abstract An extensive dataset consisting of adsorption energies of pernicious impurities present in biomass upgrading processes on common catalysts and support materials has been generated. This work aims to inform catalyst and process development for the conversion of biomass-derived feedstocks to fuels and chemicals. A high-throughput workflow was developed to execute density functional theory calculations for a diverse set of atomic (Al, B, Ca, Cl, Fe, K, Mg, Mn, N, Na, P, S, Si, Zn) and molecular (COS, H 2 S, HCl, HCN, K 2 O, KCl, NH 3 ) species on 35 unique surfaces for transition-metal (Ag, Au, Co, Cu, Fe, Ir, Ni, Pd, Pt, Re, Rh, Ru) and metal-oxide (Al 2 O 3 , MgO, anatase-TiO 2 , rutile-TiO 2 , ZnO, ZrO 2 ) catalysts and supports. Approximately 3,000 unique adsorption geometries and corresponding adsorption energies were obtained.

09 BIOMASS FUELS↗
DOE OSTI2026

Solution structures in alkali nitrates and nitrites at high concentrations

In highly concentrated electrolyte solutions, where classical models often fail, specific ion–solvent interactions dictate bulk properties. Here, in this study, we combine small-angle X-ray scattering (SAXS) and Raman spectroscopy to link molecular coordination to mesoscale structure across alkali nitrate and nitrite solutions. Our results reveal a structural hierarchy driven by cation identity in that smaller cations (Li + , Na + ) form discrete contact ion pairs, while larger cations (K + , Rb + ) promote increasingly disordered local coordination and extended solute–solvent networks. Using nitrite salts as a structural control, we confirm this organizing principle across different anion geometries. Cs + undergoes a concentration-induced transition to a uniquely ordered state, resulting in a highly structured solution. This work provides a direct, mechanistic explanation for how cation choice dictates the ‘rules’ of solution architecture, offering a predictive basis for understanding phase behavior in complex industrial and environmental systems.

Mergelsberg, Sebastian T. [Pacific Northwest Natio↗
DOE OSTI2025

Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials: Part II Resveratrol Exemplar

This SAND report summarizes work supported by an Engineering Sciences Research Foundation (ESRF) Lab Directed Research and Development (LDRD) project entitled “Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials.” This SAND report is written in two parts with Part 1 discusses recrystallization of our explosive exemplar and Part 2 summarizing our work with recrystallization of resveratrol. We have studied resveratrol recrystallization with a multiscale approach combining experiments, modeling and simulation. At the single crystal scale, microscopy experiments illuminate crystal time-dependent growth rates using advanced image analysis. Bench scale experiments were carried out to look at growth of multiple particles in a small reactor creating thousands of particles and analyzing the results with microscopy and μCT. For the modeling we combine kinetic Monte Carlo (kMC) models with subscale information from density functional theory (DFT) or molecular dynamics. This work is discussed in Part 1 and can also be found in a paper from the project discussing a coarse-grained kMC model specifically developed for resveratrol. For well-mixed systems, we have population balance equations (PBE) linked with species mass conservation forming a set of ordinary differential equations that can be solved quickly. For more complicated geometries, such as the vat crystallization used throughout the complex, a coupled computational fluid dynamic (CFD)/PBE method was developed to account for gradients in temperature and concentration and differences in crystallization rates throughout the domain. These simulations are more complex and require high performance computing. We present results for two cases: 5% seed fast cool with parameters fit to the well-mixed case and 5% seed slow cool using the same parameters. We show reasonable agreement with experiments though are particles are significantly larger than the experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2023

An Ab Initio -Derived Force Field for Amorphous Silica Interfaces for Use in Molecular Dynamics Simulations

Here, we present a classical interatomic force field, silica-DDEC, to describe the interactions of amorphous and crystalline silica surfaces, parametrized using density functional theory-based charges. Charge schemes for silica surfaces were developed using the density-derived electrostatic and chemical (DDEC) method, which reproduces atomic charges of the periodic models as well as the electrostatic potential away from the atom sites. Lennard–Jones parameters were determined by requiring the correct description of (i) the amorphous silica density, coordination defects, and local coordination geometry, relative to experimental measurements, and (ii) water-silica interatomic distances compared with ab initio results. Deprotonated surface silanol sites are also described within the model based on DDEC charges. The result is a general electronic structure-derived model for describing fully flexible amorphous and crystalline silica surfaces and interactions of liquids with silica surfaces of varying structure and protonation state.

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

Tailoring Polymer Colloids Derived Porous Carbon Spheres Based on Specific Chemical Reactions

Porous carbon spheres derived from polymer colloids with regular geometry, monodispersed morphology, well-controlled contents and structures play important roles in many areas of application, such as energy storage/conversion, gas adsorption/separation, catalysis, and chemo-photothermal therapy. Suitable polymerization reaction and synthetic strategy are both critical for the obtainment of stable polymer colloids as carbon precursors. Basic polymerization reactions are the cornerstones of synthetic strategies, which directly provides the direct molecular-based design of functionalized polymer/ carbon spheres. Thus, this progress report mainly focuses on the summary of suitable polymerization reactions for colloidal polymer derived porous carbon spheres. Recent advances in the synthetic strategies and applications are also discussed, including their corresponding polymerization reactions. Finally, the perspectives for the development of polymer derived porous carbon spheres are provided based on the controlled synthesis of polymer colloids and optimization over the carbonization process to achieve highly functionalized carbon spheres for practical applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Ultrafast Correlation Energy Estimator

A virtually no-cost method is proposed that can compute the correlation energies of general, covalently bonded, organic, and inorganic molecules (including conjugated π-electron systems) with a well-defined dominant Lewis structure at the accuracy of 99.5% of the near-exact values determined by the coupled-cluster singles, doubles, and perturbative triples [CCSD(T)] in the complete-basis-set (CBS) limit. This Correlation Energy Per Bond (CEPB) method assigns a partial correlation energy to each bond type (characterized by the identities of the two atoms forming the bond and its integer bond order) and to a lone pair, regardless of the bond length, bond angle, sp-hybridization, π-electron conjugation, ionicity, noncovalent interactions, etc. At its current stage, the method is mainly suitable for near-equilibrium geometries. The correlation energies per bond are determined by a fit to the CCSD(T)/CBS benchmarks. It can neither improve the equilibrium structures nor discern conformers or positional isomers, yet its accuracy for reaction energies rivals that of the second-order Møller–Plesset perturbation theory, which is far more expensive. Its promising performance underscores the possibility that surprisingly compact, chemically intuitive molecular fragments exist into which correlation energies can be partitioned, leading to various ultrafast correlation-energy estimators tailored to different purposes.

Binding energy↗
DOE OSTI2021

Self-similar Reynolds-averaged mechanical–scalar turbulence models for Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing in the small Atwood number limit

Analytical self-similar solutions to two-, three-, and four-equation Reynolds-averaged mechanical–scalar turbulence models describing incompressible turbulent Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing in planar geometry are derived in the small Atwood number (Boussinesq) limit. The models are based on the turbulent kinetic energy K and its dissipation rate ε, together with the scalar (heavy-fluid mass fraction) variance S and its dissipation rate Χ modeled either differentially or algebraically. The models allow for a simultaneous description of mechanical and scalar mixing, i.e., mixing layer growth and molecular mixing, respectively. Mixing layer growth parameters and other physical observables relevant to each instability are obtained explicitly as functions of the model coefficients. The turbulent fields are also expressed in terms of the model coefficients, with their temporal power-law scalings obtained by requiring that the self-similar equations are explicitly time-independent. The model calibration methodology is described and discussed. Expressions for a subset of the various physical observables are used to calibrate each of the two-, three-, and four-equation models, such that the self-similar solutions are consistent with experimental and numerical simulation data corresponding to these values of the observables and to specific canonical Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz turbulent flows. A calibrated four-equation model is then used to reconstruct the mean and turbulent fields, and late-time turbulent equation budgets for each instability-induced flow across the mixing layer. The reference solutions derived here can provide systematic calibrations and better understanding of mechanical–scalar turbulence models and their predictions for instability-induced turbulent mixing in the very large Reynolds number limit.

42 ENGINEERING↗
DOE OSTI2022

Closing the loop between microstructure and charge transport in conjugated polymers by combining microscopy and simulation

A grand challenge in materials science is to identify the impact of molecular composition and structure across a range of length scales on macroscopic properties. We demonstrate a unified experimental–theoretical framework that coordinates experimental measurements of mesoscale structure with molecular-level physical modeling to bridge multiple scales of physical behavior. Here we apply this framework to understand charge transport in a semiconducting polymer. Spatially-resolved nanodiffraction in a transmission electron microscope is combined with a self-consistent framework of the polymer chain statistics to yield a detailed picture of the polymer microstructure ranging from the molecular to device relevant scale. Using these data as inputs for charge transport calculations, the combined multiscale approach highlights the underrepresented role of defects in existing transport models. Short-range transport is shown to be more chaotic than is often pictured, with the drift velocity accounting for a small portion of overall charge motion. Local transport is sensitive to the alignment and geometry of polymer chains. At longer length scales, large domains and gradual grain boundaries funnel charges preferentially to certain regions, creating inhomogeneous charge distributions. While alignment generally improves mobility, these funneling effects negatively impact mobility. The microstructure is modified in silico to explore possible design rules, showing chain stiffness and alignment to be beneficial while local homogeneity has no positive effect. This combined approach creates a flexible and extensible pipeline for analyzing multiscale functional properties and a general strategy for extending the accesible length scales of experimental and theoretical probes by harnessing their combined strengths.

42 ENGINEERING↗
DOE OSTI2020

Improved chemistry restraints for crystallographic refinement by integrating the Amber force field into Phenix

The refinement of biomolecular crystallographic models relies on geometric restraints to help to address the paucity of experimental data typical in these experiments. Limitations in these restraints can degrade the quality of the resulting atomic models. Here, an integration of the full all-atom Amber molecular-dynamics force field into Phenix crystallographic refinement is presented, which enables more complete modeling of biomolecular chemistry. The advantages of the force field include a carefully derived set of torsion-angle potentials, an extensive and flexible set of atom types, Lennard–Jones treatment of nonbonded interactions and a full treatment of crystalline electrostatics. The new combined method was tested against conventional geometry restraints for over 22 000 protein structures. Structures refined with the new method show substantially improved model quality. On average, Ramachandran and rotamer scores are somewhat better, clashscores and MolProbity scores are significantly improved, and the modeling of electrostatics leads to structures that exhibit more, and more correct, hydrogen bonds than those refined using traditional geometry restraints. In general it is found that model improvements are greatest at lower resolutions, prompting plans to add the Amber target function to real-space refinement for use in electron cryo-microscopy. This work opens the door to the future development of more advanced applications such as Amber -based ensemble refinement, quantum-mechanical representation of active sites and improved geometric restraints for simulated annealing.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Carbon Kagome nanotubes—quasi-one-dimensional nanostructures with flat bands

In recent years, a number of bulk materials and heterostructures have been explored due their connections with exotic materials phenomena emanating from flat band physics and strong electronic correlation. The possibility of realizing such fascinating material properties in simple realistic nanostructures is particularly exciting, especially as the investigation of exotic states of electronic matter in wire-like geometries is relatively unexplored in the literature. Motivated by these considerations, we introduce in this work carbon Kagome nanotubes (CKNTs)—a new allotrope of carbon formed by rolling up Kagome graphene, and investigate this material using specialized first principles calculations. We identify two principal varieties of CKNTs—armchair and zigzag, and find both varieties to be stable at room temperature, based on ab initio molecular dynamics simulations. CKNTs are metallic and feature dispersionless states (i.e., flat bands) near the Fermi level throughout their Brillouin zone, along with an associated singular peak in the electronic density of states. We calculate the mechanical and electronic response of CKNTs to torsional and axial strains, and show that CKNTs appear to be more mechanically compliant than conventional carbon nanotubes (CNTs). Additionally, we find that the electronic properties of CKNTs undergo significant electronic transitions—with emergent partial flat bands and tilted Dirac points—when twisted. We develop a relatively simple tight-binding model that can explain many of these electronic features. We also discuss possible routes for the synthesis of CKNTs. Overall, CKNTs appear to be unique and striking examples of realistic elemental quasi-one-dimensional materials that may display fascinating material properties due to strong electronic correlation. Distorted CKNTs may provide an interesting nanomaterial platform where flat band physics and chirality induced anomalous transport effects may be studied together.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Role of mechanical stress localizations on the radiation hardness of AlGaN/GaN high electron mobility transistors

Multi-material, multi-layered systems such as AlGaN/GaN high electron mobility transistors (HEMTs) contain residual mechanical stresses that arise from sharp contrasts in device geometry and materials parameters. These stresses, which can be either tensile or compressive, are difficult to detect and eliminate because of their highly localized nature. We propose that their high-stored internal energy makes potential sites for defect nucleation sites under radiation, particularly if their locations coincide with the electrically sensitive regions of a transistor. In this study, we validate this hypothesis with molecular dynamic simulation and experiments exposing both pristine and annealed HEMTS to 2.8 MeV Au +3 irradiation. Our unique annealing process uses mechanical momentum of electrons, also known as the electron wind force (EWF) to mitigate the residual stress at room temperature. High-resolution transmission electron microscopy and cathodoluminescence spectra reveal the reduction of point defects and dislocations near the two-dimensional electron gas region of EWF-treated devices compared to pristine devices. The EWF-treated HEMTs showed relatively higher resilience with approximately 10% less degradation of drain saturation current and ON-resistance and 5% less degradation of peak transconductance. Both mobility and carrier concentration of the EWF-treated devices were less impacted compared to the pristine devices. Our results suggest that the lower density of nanoscale stress localization contributed to the improved radiation tolerance of the EWF-treated devices. Intriguingly, the EWF is found to modulate the defect distribution by moving the defects to electrically less sensitive regions in the form of dislocation networks, which act as sinks for the radiation induced defects and this assisted faster dynamic annealing.

AlGaN/GaN HEMTs↗
DOE OSTI2021

Radial Two-Dimensional Ion Crystals in a Linear Paul Trap

We experimentally study two-dimensional (2D) Coulomb crystals in the “radial-2D” phase of a linear Paul trap. This phase is identified by a 2D ion lattice aligned entirely with the radial plane and is created by imposing a large ratio of axial to radial trapping potentials. Using arrays of up to 19 171 Yb + ions, we demonstrate that the structural phase boundaries of such crystals are well described by the pseudopotential approximation, despite the time-dependent ion positions driven by intrinsic micromotion. We further observe that micromotion-induced heating of the radial-2D crystal is confined to the radial plane. Finally, we verify that the transverse motional modes, which are used in most ion-trap quantum simulation schemes, are well-predictable numerically and remain decoupled and cold in this geometry. Our results establish radial-2D ion crystals as a robust experimental platform for realizing a variety of theoretical proposals in quantum simulation and computation.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Modulating Transition Metal Reactivity with Force

The reactivity and selectivity of a transition metal catalyst is intimately related to its ligand-sphere geometry, and, in many cases, the ideal ligand geometry for one step of a catalytic cycle is poorly matched to the ideal ligand geometry for another. For this reason, methods for reversibly modulating ligand geometry on the time scale of catalytic turnover or monomer enchainment are highly desirable. Mechanical force represents a heretofore untapped approach to modulate catalyst geometry and/or reactivity, with the potential to do so on the timescale of catalytic turnover or monomer enchainment. Macroscopic mechanical forces are large, directional and localized to an extent that differentiates them from other forms of energy input such as heat or light. Here, in this Concept, we describe our efforts to address the fundamental challenges associated with force-modulated transition metal catalysis by employing molecular force probe ligands comprising a stiff stilbene photoswitch tethered to rotationally flexible biaryl bisphosphine ligand. Our efforts to date include the modulation of catalytic activity through force-mediated ligand perturbations, quantification of the force-coupled ligand effects on the energetics of elementary organometallic transformations, and evaluation of the mechanisms of force transduction in these systems.

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

All-electron molecular tunnel ionization based on the weak-field asymptotic theory in the integral representation

Tunnel ionization (TI) underlies many important ultrafast processes, such as high-harmonic generation andstrong-field ionization. Among the existing theories for TI, many-electron weak-field asymptotic theory (ME-WFAT) is by design capable of accurately treating many-electron effects in TI. An earlier version of ME-WFATrelied on an accurate representation of the asymptotic tail of the orbitals, which hindered its implementation inGaussian-basis-set-based quantum chemistry programs. In this work, we reformulate ME-WFAT in the integralrepresentation, which makes the quality of the asymptotic tail much less critical, hence greatly facilitating itsimplementation in standard quantum chemistry packages. The integral reformulation introduced here is thereforemuch more robust when applied to molecules with arbitrary geometry. Here, we present several case studies, amongwhich is the CO molecule where some earlier theories disagree with experiments. Here we find that ME-WFATproduces the largest ionization probability when the field points from C to O, as experiments suggest. Anattractive feature of ME-WFAT is that it can be used with various types of multielectron methods whether ofdensity functional or multiconfiguration types, this inturn facilitates tunnel ionization calculation in systems exhibiting a strong multireference character.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Distinct Mechanisms of Triplet Pair Decay in Amorphous and Crystalline Heteroacene Thin Films

Triplet pairs (TT) in crystalline molecular semiconductors have unique spin properties of interest for quantum information or enhancing solar photoconversion. The population and diffusion dynamics of TT have been the subject of recent studies, both in covalent dimers and in crystalline systems. Here, in this study, we monitor the triplet population in neat polycrystalline and amorphous films of a heteroacene with known TT spectral properties and tunable spin polarization depending on the intermolecular geometry. Transient measurements reveal an anomalous power dependence in polycrystalline films that we attribute to the fast diffusion and interaction of dissociated triplet pairs confined to one-dimensional stacks of strongly coupled molecules. The nongeminate triplet interaction after dephasing facilitates conversion to the triplet 3 TT and eventually T 1 +S 0 . Amorphous films have no power dependence and proceed directly from 1 TT to 3 TT and subsequently T 1 +S 0 via state mixing facilitated by nonparallel geometries and weak exchange coupling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗