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.

SourceResearch connectionFrom yearThrough year

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

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Investigation of the Effect of Framework Flexibility on Adsorption in SIFSIX-3-Cu using a Machine-Learned Force Field

Metal-organic frameworks (MOFs) are a promising class of adsorbents. The performance of MOF sorbents relies on high selectivity and low regeneration energy. This work focuses on the use of machine learned force fields (MLFFs) to model adsorption in a flexible MOF, SIFSIX-3-Cu. A DeePMD-based MLFF was trained to reproduce DFT (PBE+D3) energies, forces, and stresses, using an iterative sampling scheme combining sampling based on molecular dynamics, Monte Carlo, and geometry optimization to capture both attractive and repulsive regions of the potential energy surface. Flexibility of the MOF was explicitly included in this model. Hybrid Monte Carlo/molecular dynamics (MC/MD) simulations using the MLFF predicted adsorption isotherms in good agreement with experimental data for a range of pressures (40 Pa – 104 Pa) in contrast to rigid models, which overpredict CO2 adsorption at low pressures. The improvement was the result of a description of the variability of fluorine-fluorine diagonal distances at adsorption sites. This detailed description of flexibility afforded by the MLFF resulted in more accurate predictions adsorption isotherms when compared to the experimentally measured values. These results underscore the importance of including framework flexibility when modeling adsorption phenomena in MOFs, particularly for low pressure applications and provide a robust procedure for training MLFF models for MOFs.

Atomistic Simulation↗

DOE OSTI2023

Guiding Principles for the Rational Design of Hybrid Materials: Use of DFT Methodology for Evaluating Non‐Covalent Interactions in a Uranyl Tetrahalide Model System

Abstract Together with the synthesis and experimental characterization of 14 hybrid materials containing [UO 2 X 4 ] 2− (X=Cl − and Br − ) and organic cations, we report on novel methods for determining correlation trends in their formation enthalpy (ΔH f ) and observed vibrational signatures. Δ H f values were analyzed through isothermal acid calorimetry and a Density Functional Theory+Thermodynamics (DFT+T) approach with results showing good agreement between theory and experiment. Three factors (packing efficiency, cation protonation enthalpy, and hydrogen bonding energy [ ]) were assessed as descriptors for trends in Δ H f . Results demonstrated a strong correlation between and Δ H f , highlighting the importance of hydrogen bonding networks in determining the relative stability of solid‐state hybrid materials. Lastly, we investigate how hydrogen bonding networks affect the vibrational characteristics of uranyl solid‐state materials using experimental Raman and IR spectroscopy and theoretical bond orders and find that hydrogen bonding can red‐shift U≡O stretching modes. Overall, the tightly integrated experimental and theoretical studies presented here bridge the trends in macroscopic thermodynamic energies and spectroscopic features with molecular‐level details of the geometry and electronic structure. This modeling framework forms a basis for exploring 3D hydrogen bonding as a tunable design feature in the pursuit of supramolecular materials by rational design.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Guiding Principles for the Rational Design of Hybrid Materials: Use of DFT Methodology for Evaluating Non‐Covalent Interactions in a Uranyl Tetrahalide Model System

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Modeling Molecules under Pressure with Gaussian Potentials

The computational modeling of molecules under high pressure is a growing research area that augments experimental high-pressure chemistry. Here, a new electronic structure method for modeling atoms and molecules under pressure, Gaussians On Surface Tesserae Simulate HYdrostatic Pressure (GOSTSHYP) approach, is introduced. In this method, a set of Gaussian potentials is distributed evenly on the van der Waals surface of the investigated chemical system, leading to a compression of the electron density and the atomic scaffold. Since no parameters other than pressure need to be specified, GOSTSHYP allows straightforward geometry optimizations and ab initio molecular dynamics simulations of chemical systems under pressure for nonexpert users. In this work, calculated energies, bond lengths, and dipole moments under pressure fall within the range of established computational methods for high-pressure chemistry. A Diels-Alder reaction and the cyclotrimerization of acetylene showcase the ability of GOSTSHYP to model pressure-induced chemical reactions. The connection to mechanochemistry is pointed out.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Tris(carbene)borates; alternatives to cyclopentadienyls in organolanthanide chemistry

The chemistry of the tris-carbene anion phenyltris(3-alkyl-imidazoline-2-yliden-1-yl)borate, [C3 Me ] - ligand, is initiated for f-block metal cations. Neutral, molecular complexes of the form Ln(C3) 2 I are formed for cerium(III), while a separated ion pair [Ln(C3) 2 ]I forms for ytterbium(III). DFT/QTAIM computational analyses of the complexes and related tridentate tris(pyrazolyl)borate (Tp) – supported analogs demonstrates the anticipated strength of the σ donation and confirms greater covalency in the metal–carbon bonds of the [C3 Me ] - complexes in comparison with those in the Tp Me,Me complexes. The DFT calculations demonstrate the crucial role of THF solvent in accurately reproducing the contrasting molecular and ion-pair geometries observed experimentally for the Ce and Yb complexes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Reduced scaling formulation of CASPT2 analytical gradients using the supporting subspace method

We present a reduced scaling and exact reformulation of state specific complete active space second-order perturbation (CASPT2) analytical gradients in terms of the MP2 and Fock derivatives using the supporting subspace method. This work follows naturally from the supporting subspace formulation of the CASPT2 energy in terms of the MP2 energy using dressed orbitals and Fock builds. For a given active space configuration, the terms corresponding to the MP2-gradient can be evaluated with O(N5) operations, while the rest of the calculations can be computed with O(N3) operations using Fock builds, Fock gradients, and linear algebra. When tensor-hyper-contraction is applied simultaneously, the computational cost can be further reduced to O(N4) for a fixed active space size. The new formulation enables efficient implementation of CASPT2 analytical gradients by leveraging the existing graphical processing unit (GPU)-based MP2 and Fock routines. We present benchmark results that demonstrate the accuracy and performance of the new method. Example applications of the new method in ab initio molecular dynamics simulation and constrained geometry optimization are given.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Geometry optimizations with spinor-based relativistic coupled-cluster theory

Development of analytic gradients for relativistic coupled-cluster singles and doubles augmented with a non-iterative triples [CCSD(T)] method using an all-electron exact two-component Hamiltonian with atomic mean-field spin–orbit integrals (X2CAMF) is reported. This enables efficient CC geometry optimizations with spin–orbit coupling included in orbitals. Here, the applicability of the implementation is demonstrated using benchmark X2CAMF-CCSD(T) calculations of equilibrium structures and harmonic vibrational frequencies for methyl halides, CH 3 X (X = Br, I, and At), as well as calculations of rotational constants and infrared spectrum for RaSH+, a radioactive molecular ion of interest to spectroscopic study.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

QMMM 2023: A program for combined quantum mechanical and molecular mechanical modeling and simulations

Combined quantum mechanical and molecular mechanical (QM/MM) methods play an important role in multiscale modeling and simulations. QMMM 2023 is a general-purpose program for single-point calculations, geometry optimizations, transition-state optimizations, and molecular dynamics (MD) at the QM/MM level. It calls a QM package and an MM package to perform the required single-level calculations and combines them into a QM/MM energy by a variety of schemes. QMMM 2023 supports GAMESS-US, Gaussian, and ORCA as QM packages and Tinker as the MM package. Four types of treatments are available for embedding the QM subsystem in the MM environment: mechanical embedding with gas-phase calculations of the QM region, electronic embedding that allows polarization of the QM region by the MM environment, polarizable embedding for mutual polarization of the QM and MM regions, and flexible embedding for both mutual polarization and partial charge transfer between the QM and MM regions. Boundaries between QM and MM regions that pass through covalent bonds can be treated by several methods, including the redistributed charge (RC) scheme, redistributed charge and dipole (RCD) scheme, balanced-RC scheme, balanced-RCD scheme, screened charge scheme that takes account of charge penetration effects, and smeared charge scheme that delocalizes the MM charges near the QM–MM boundary. Geometry optimization can be done using the optimizer implemented in QMMM 2023 or the Berny optimizer in Gaussian through external calls to Gaussian. Molecular dynamics simulations can be performed at the pure-MM level, pure-QM level, fixed-partitioning QM/MM level, and adaptive-partitioning QM/MM level. As a result, the adaptive-partitioning treatments permit on-the-fly relocation of the QM–MM boundary by dynamically reclassifying atoms or groups into the QM or MM subsystems.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Machine learning prediction of glass transition temperature of conjugated polymers from chemical structure

Predicting the glass transition temperature (T g ) is of critical importance as it governs the thermomechanical performance of conjugated polymers (CPs). Here, we report a predictive modeling framework to predict T g of CPs through the integration of machine learning (ML), molecular dynamics (MD) simulations, and experiments. With 154 T g data collected, an ML model is developed by taking simplified “geometry” of six chemical building blocks as molecular features, where side-chain fraction, isolated rings, fused rings, and bridged rings features are identified as the dominant ones for T g . MD simulations further unravel the fundamental roles of those chemical building blocks in dynamical heterogeneity and local mobility of CPs at a molecular level. The developed ML model is demonstrated for its capability of predicting T g of several new high-performance solar cell materials to a good approximation. The established predictive framework facilitates the design and prediction of T g of complex CPs, paving the way for addressing device stability issues that have hampered the field from developing stable organic electronics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Modeling Pb(II) Adsorption on Mineral Surfaces: Bridging Density Functional Theory and Experiment with Thermodynamic Insights

Despite decades of work on aqueous lead (Pb) adsorption on a-Fe2O3 (hematite) and a-Al2O3 (alumina), gaps between measurements and modeling obscure molecular-level understanding. Achieving well-matched geometries between theory and experimental for mineral-water interfaces is a hurdle, as surface functional group type and distribution must be accounted for in determining mechanisms. Additionally, computational methods that can describe the substrate are often not appropriate to capture aqueous effects. Progress requires focusing on well-studied and relevant systems, such as key facets (001),(012), and (110) of hematite and alumina, and ubiquitous contaminants such as aqueous Pb. In the past, bulk-parametrized bond-valence principles were used to rationalize Pb(II) adsorption trends. These approaches can break down at surfaces, where flexible bonding environments and adsorption-induced surface relaxations play a critical role. Here, we adapt and apply a density functional theory (DFT) and thermodynamics framework, integrating DFT-calculated energies with experimental data and electrochemical principles, to predict Pb(II) adsorption. Our model results capture trends across the full set of surfaces and predict that inner-sphere Pb(II) sorption on (001) alumina varies from unfavorable to weakly favorable across a range of pH conditions. This aligns with experiment insights that Pb(II) interacts at that surface through outer-sphere interactions. Extending to Fe(II) adsorption, we demonstrate a coverage-dependent site preference, potentially explaining disorder in overlayers grown by the oxidative adsorption of Fe(II) on hematite (001).

lead contamination↗

DOE OSTI2024

Comparison of Machine Learning Approaches for Prediction of the Equivalent Alkane Carbon Number for Microemulsions Based on Molecular Properties

The chemical properties of oils are vital in the design of microemulsion systems. The hydrophilic–lipophilic difference equation used to predict microemulsions’ phase behavior expresses the oils’ physiochemical properties as the equivalent alkane carbon number (EACN). The experimental determination of EACN requires knowledge of the temperature dependence of the microemulsion system and the effects of different surfactant concentrations. Thus, the experimental determination is time-intensive and tedious, requiring days to months for proper separations. Furthermore, the experiments require high purity of chemicals because microemulsions are sensitive to impurities. Our work focuses on the quick and reliable predictions of the EACN with machine learning (ML) models. Due to the immaturity of ML chemical predictions, we compare three graph neural networks (GNNs) and a gradient-boosted tree algorithm, known as XGBoost. The GNNs use the molecular structures represented as simplified molecular-input line-entry system (SMILES) codes for the initial input, which allows us to assess whether geometry optimization is necessary for reliable results. The XGBoost model also begins with the SMILES representations of the molecules but uses molecular descriptors instead of geometry optimizations. As a result, the best model tested (crystal graph convolutional neural network with Merck molecular force field-94) has an error of 1.15 EACN units of the true EACN for unknown data with the errors skewed toward zero and an R² score of 0.9

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Orientation-Dependent Nanomechanical Behavior of Pentaerythritol Tetranitrate as Probed by Multiple Nanoindentation Tip Geometries

Nanoindentation can be leveraged to aid in the high fidelity modeling of dislocation mediated plasticity in pentaerythritol tetranitrate (PETN), an anisotropic energetic molecular crystal. Moreover, nanoindentation tip parameters such as tip geometry, size, and degree of acuity can be utilized to target anisotropic behavior. In this work, nanoindentation was conducted across a range of orientations on the (110) face of PETN to characterize resultant yield behavior, mechanical property measurements, and resultant slip behavior and fracture initiation. Three different indentation tips were utilized: a 3-sided pyramidal Berkovich tip, a 4-sided high aspect ratio Knoop tip, and a 90° conical tip. Ultimately, indenter tip radius was documented to impact yield behavior, whereas tip geometry affected larger scale processes such as slip, and tip acuity was the dominating factor that led to fracture. The axisymmetric conical tip, serving as a baseline, showed the least amount of variation in mechanical property measurements but also the largest distribution of maximum shear stress at which initial yielding occurred. Its high degree of acuity, however, was more prone to induce fracture at higher loads. The Knoop tip was shown to be suitable for average measurements, but also for elucidation of certain anisotropic features. A distinctly higher perceived hardness at 45° was measured with the Knoop tip, indicating less dislocation motion in that direction also observed in this work via scanning probe microscopy. Lastly, the commonly used Berkovich tip was a good compromise whereby it provided a representative volume element describing the average behavior of the material. These results can be utilized to target desired anisotropic behavior in a wider range of molecular crystals, as well as to inform theoretical considerations for dislocation mediated plasticity in PETN.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Nonadditive Interactions Unlock Small-Particle Mobility in Binary Colloidal Monolayers

We examine the organization and dynamics of binary colloidal monolayers composed of micron-scale silica particles interspersed with smaller-diameter silica particles that serve as minority component impurities. These binary monolayers are prepared at the surface of ionic liquid droplets over a range of size ratios (σ = 0.16–0.66) and are studied with low-dose minimally perturbative scanning electron microscopy (SEM). The high resolution of SEM imaging provides direct tracking of all particle coordinates over time, enabling a complete description of the microscopic state. In these bidisperse size mixtures, particle interactions are nonadditive because interfacial pinning to the droplet surface causes the equators of differently sized particles to lie in separate planes. By varying the size ratio, we control the extent of nonadditivity in order to achieve phase behavior inaccessible to additive 2D systems. Across the range of size ratios, we tune the system from a mobile small-particle phase (σ < 0.24) to an interstitial solid (0.24 < σ < 0.33) and furthermore to a disordered glass (σ > 0.33). These distinct phase regimes are classified through measurements of hexagonal ordering of the large-particle host lattice and the lattice’s capacity for small-particle transport. Altogether, we explain these structural and dynamic trends by considering the combined influence of interparticle interactions and the colloidal packing geometry. Our measurements are reproduced in molecular dynamics simulations of 2D nonadditive disks, suggesting an efficient method for describing confined systems with reduced dimensionality representations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Low Energy Implantation into Transition-Metal Dichalcogenide Monolayers to Form Janus Structures

Atomically-thin two-dimensional (2D) materials face significant energy barriers for synthesis and processing into functional metastable phases such as Janus structures. Here, in this study, the controllable implantation of hyperthermal species from pulsed laser deposition (PLD) plasmas is introduced as a top-down method to compositionally engineer 2D monolayers. The kinetic energies of Se clusters impinging on suspended monolayer WS 2 crystals were controlled in the <10 eV/atom range with in situ plasma diagnostics to determine the thresholds for selective top layer replacement of sulfur by selenium for the formation of high quality WSSe Janus monolayers at low (300 °C) temperatures, and bottom layer replacement for complete conversion to WSe 2 . Atomic-resolution electron microscopy and spectroscopy in tilted geometry confirm the WSSe Janus monolayer. Molecular dynamics simulations reveal that Se clusters implant to form disordered metastable alloy regions, which then recrystallize to form highly ordered structures, opening the door for low-energy implantation by PLD as a novel method to explore the synthesis of 2D Janus layers and alloys of variable composition.

2D materials↗

DOE OSTI2024

Driving macro-scale transformations in three-dimensional-printed biopolymers through controlled induction of molecular anisotropy at the nanoscale

Motivated by the need to harness the properties of renewable and biodegradable polymers for the design and manufacturing of multi-scale structures with complex geometries, we have employed our additive manufacturing platform that leverages molecular self-assembly for the production of metre-scale structures characterized by complex geometries and heterogeneous material composition. As a precursor material, we used chitosan, a chemically modified form of chitin, an abundant and sustainable structural polysaccharide. We demonstrate the ability to control concentration-dependent crystallization as well as the induction of the preferred orientation of the polymer chains through the combination of extrusion-based robotic fabrication and directional toolpathing. Anisotropy is demonstrated and assessed through high-resolution micro-X-ray diffraction in conjunction with finite element simulations. Using this approach, we can leverage controlled and user-defined small-scale propagation of residual stresses to induce large-scale folding of the resulting structures.

36 MATERIALS SCIENCE↗

DOE OSTI2022

The Molecular Ordering and Double-Channel Carrier Generation of Nonfullerene Photovoltaics within Multi-Length-Scale Morphology

The success of nonfullerene acceptor (NFA) solar cells lies in their unique physical properties beyond the extended absorption and suitable energy levels. The current study investigates the morphology and photophysical behavior of PBDB-T donor blending with ITIC, 4TIC, and 6TIC acceptors. Single-crystal study shows that the π–π stacking and side-chain interaction dictate molecular assembly, which can be carried to blended films, forming a multi-length-scale morphology. Spontaneous carrier generation is seen in ITIC, 4TIC, and 6TIC neat films and their blended thin films using the PBDB-T donor, providing a new avenue of zero-energy-loss carrier formation. The molecular packing associated with specific contacts and geometry is key in influencing the photophysics, as demonstrated by the charge transfer and carrier lifetime results. The 2D layer of 6TIC facilitates the exciton-to-polaron conversion, and the largest photogenerated polaron yield is obtained. Finally, the new mechanism, together with the highly efficient blending region carrier generation, has the prospect of the fundamental advantage for NFA solar cells, from molecular assembly to thin-film morphology.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Bond-centric modular design of protein assemblies

Directional interactions that generate regular coordination geometries are a powerful means of guiding molecular and colloidal self-assembly, but implementing such high-level interactions with proteins remains challenging due to their complex shapes and intricate interface properties. Here we describe a modular approach to protein nanomaterial design inspired by the rich chemical diversity that can be generated from the small number of atomic valencies. We design protein building blocks using deep learning-based generative tools, incorporating regular coordination geometries and tailorable bonding interactions that enable the assembly of diverse closed and open architectures guided by simple geometric principles. Experimental characterization confirms the successful formation of more than 20 multicomponent polyhedral protein cages, two-dimensional arrays and three-dimensional protein lattices, with a high (10%–50%) success rate and electron microscopy data closely matching the corresponding design models. Due to modularity, individual building blocks can assemble with different partners to generate distinct regular assemblies, resulting in an economy of parts and enabling the construction of reconfigurable networks for designer nanomaterials.

Biomaterials – proteins↗

DOE OSTI2020

Molecular-Level Exploration of the Structure-Function Relations Underlying Interfacial Charge Transfer in the Subphthalocyanine/C 60 Organic Photovoltaic System

The arrangement of organic molecules at the donor-acceptor interface in an organic photovoltaic (OPV) cell can have a strong effect on the generation of charge carriers and thereby cell performance. In this paper, we report the molecular-level exploration of the ensemble of interfacial donor-acceptor pair geometries and the charge-transfer (CT) rates to which they give rise. Our approach combines molecular-dynamics simulations, electronic structure calculations, machine learning, and rate theory. This approach is applied to the boron subphthalocyanine chloride (donor) and C 60 (acceptor) OPV system. In this work we find that the interface is dominated by a previously unreported donor-acceptor pair edge geometry, which contributes significantly to device performance in a manner that depends on the initial conditions. Quantitative relations between the morphology and CT rates are established, which can be used to advance the design of more efficient OPV devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗