Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Molecular solids”

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

Topical Report on Water Balance Model Development and Initial Application

This topical report summarizes results from developing a generic dynamic water balance tool which can be used by coal-fired power plant operators to evaluate FGD water management strategies. EPRI is using a statistical software model created on the GoldSim platform and data from coal-fired energy plants to determine the fraction of power plant water needs that can be satisfied by electrodialysis reversal treatment and, correspondingly, the discharge volume and/or byproducts that still need to be managed. The generic model is designed to accept a number of easily obtainable process variables and generate probabilities of resultant water characteristics. The model is also constructed as a universal process flow; one that can be adjusted to accommodate most major components of an FGD blowdown handling system. The generic model is dynamic and utilizes user-inputted operational variation in water quality and quantity to best provide a probabilistic representation of potential performance outcomes. In order to account for complexities of water composition and treatment, a separate chemical modeling software is used to evaluate water chemistry and unit operation performance. The thermodynamic modeling software, OLI Studio, is employed to convert ion data from laboratory analyses into balanced molecular concentrations that are classified as dissolved and suspended solids. The dynamic water model is designed to import a total of 11 species commonly found in FGD blowdown water as molecular concentrations.Water management unit operations are used in the model and categorized by either physical unit operation blocks or chemical unit operations blocks. As the name implies, the physical water treatment blocks are unit operations that are exclusively physical changes to the system (e.g., separations and volume changes). The chemical water treatment blocks are unit operations where chemical changes occur and require the use of OLI Studio to characterize subsequent streams. Water and material balances from three unique coal-fired power plants were obtained and preliminary modeling on the dynamic water balance tool was performed using the extracted data from these power plants. These three power plants were chosen due to their unique configurations and water profiles that represented a variety of operational changes. Results from using the dynamic tool with the three plant configurations are presented and discussed in the context of the existing tool and planned future capabilities.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2022

Multiphase equation of state for Ta 2 O 5

A new equation of state for Ta 2 O 5 is presented. The EOS is constructed using the OpenSesame software and is referred to as SESAME 3530. The EOS uses a combination of density functional theory (DFT) calculations and experimental data. DFT calculations include cold curves and phonons of the solid phases, as well as DFT-based molecular dynamics simulations of the liquid phase. Experimental data includes isobaric, diamond anvil cell, and porous shock Hugoniot data. To fit the data, we create a multiphase EOS consisting of two solid phases and the liquid. Overall agreement with experimental data is shown, and we provide some suggestions for future experiments that could improve our knowledge of the phase diagram.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Experimental Constraints on Solid Nitride Phases in Rocky Mantles of Reduced Planets and Implications for Observable Atmosphere Compositions

Abstract Astronomical surveys have discovered thousands of transiting exoplanets, revealing that rocky planets are common in the galaxy. A planet's interior chemistry is frequently inferred by average density, described by mass‐radius (M‐R) relationships. However, M‐R relationships give rise to non‐unique interpretations of a planet's interior composition, an issue that limits our ability to characterize far‐away worlds. We present experimental and density functional theoretical results addressing the influence of an ultra‐reducing (oxygen‐poor) interior chemistry on rocky mantle phases and discuss the possible implications for atmospheric observables. We show that silicon carbide (SiC) and molecular nitrogen (N 2 ) react to form solid silicon nitride (γ‐Si 3 N 4 ) at high pressures and high temperatures in a laser‐heated diamond‐anvil cell, consistent with ab initio computations. Si 3 N 4 remains stable under extreme conditions and when quenched to ambient conditions. As SiC is a common compound found under very reducing conditions, these results indicate that nitrogen may form solid phases in an oxygen‐poor rocky planet. If, by sequestering nitrogen in a planet's mantle, the distribution of nitrogen between a planet's interior and atmosphere is altered (i.e., a nitrogen‐rich mantle and nitrogen‐poor atmosphere), these results indicate that there may be atmospheric observables connected to the mantle‐redox state of a rocky planet besides the oxygen‐containing phases ubiquitous in exoplanet literature.

Daviau, Kierstin↗
DOE OSTI2020

Intrinsic chemical reactivity of solid-electrolyte interphase components in silicon–lithium alloy anode batteries probed by FTIR spectroscopy

Herein we report the solid reaction products from the chemical reaction of aprotic battery electrolyte and three purported components of the Si-based anode SEI : SiO 2 nanoparticles (NPs), lithium silicate (Li x SiO y ) powders, and Si NPs. We use FTIR and classical molecular dynamics/density functional perturbation theory to assess the solid products remaining with these model materials after exposure to electrolyte. The absence of electrochemical bias provides a view of the chemical speciation resulting from early-stage chemical reactivity during battery assembly as well as under open circuit storage conditions. We believe these species represent the initial stages of SEI growth and predict they likely drive subsequent chemical and electrochemical reactions by controlling molecular interactons at the Si active material interface. We find that nominally equivalent materials react differently even before any electrochemistry is performed (e.g., acidic SiO 2 dissolves whereas alkaline SiO 2 is relatively robust), and derive new understanding of the chemical species that could and could not form stable SEI components in Si-based anodes. These results can be used to inform how to passivate Si anode surfaces and potentially generate an artificially engineered SEI that would be stable and enable next-generation battery anodes.

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

Orientational Order in Spin-Cast Lead-Iodide Perovskite Nanocrystal Solids

Combined synthetic control over size and composition renders colloidal lead-halide perovskite nanocrystals a tunable platform for high-efficiency optoelectronic applications. However, the properties and operational stability of devices based on nanocrystal solids are often dictated by the method of the evaporation-induced assembly. Ubiquitous slow evaporation techniques can produce highly ordered nanocrystal domains but limit the prospects for scalable fabrication of continuous device layers, calling for investigation of approaches to more rapidly form ordered perovskite nanocrystal solids. Here, we study orientationally ordered lead-iodide perovskite nanocrystal solids prepared by conventional spin coating with molecular additives (excess ligand) to enhance ordering within the arrays. In situ X-ray scattering measurements reveal that orientational ordering occurs rapidly upon solvent removal during spin coating and can be further enhanced by manipulating the spin speed. We vary the additive ligand length and explore trade-offs between ordering and layered perovskite impurity formation. Arrays treated with the intermediate-length octylamine ligand exhibit increased in-plane electronic conductivity, suggesting orientational ordering and internanocrystal electronic coupling can be enhanced by the treatment. In conclusion, these results highlight the prospects of establishing long-range order in lead-halide perovskite nanocrystal solids by using simple and fast coating methods.

14 SOLAR ENERGY↗
DOE OSTI2021

A Multi-modal Approach to Understanding Degradation of Organic Photovoltaic Materials

State-of-the-art organic photovoltaic (OPV) materials are composed of complex, chemically diverse polymeric and molecular structures that form highly intricate solid-state interactions, collectively yielding exceptional tunability in performance and aesthetics. These properties are especially attractive for semitransparent power-generating windows or shades in living environments, greenhouses, or other architectural integrations. However, before such a future is realized, a broader and deeper understanding of property stability must be acquired. Stability during operating and environmental conditions is critical, namely, material color steadfastness, optoelectronic performance retention, morphological rigidity, and chemical robustness. To date, no single investigation encompasses all four distinct, yet interconnected, metrics. Here, we present a multimodal strategy that captures a dynamic and interconnected evolution of each property during the course of an accelerated photobleaching experiment. We demonstrate this approach across relevant length scales (from molecular to visual macroscale) using X-ray photoelectron spectroscopy, grazing-incidence X-ray scattering, microwave conductivity, and time-dependent photobleaching spectroscopies for two high-performance semitransparent OPV blends—PDPP4T:PC 60 BM and PDPP4T:IEICO-4F, with comparisons to the stabilities of the individual components. We present direct evidence that specific molecular acceptor (fullerene vs nonfullerene) designs and the resulting donor–acceptor interactions lead to distinctly different mechanistic routes that ultimately arrive at what is termed “OPV degradation.” We directly observe a chemical oxidation of the cyano endcaps of the IEICO-4F that coincides with a morphological change and large loss in photoconductivity while the fullerene acceptor-containing blend demonstrates a significantly greater fraction of oxygen uptake but retains 55% of the photoconductivity. Finally, this experimental roadmap provides meaningful guidance for future high-throughput, multimodal studies, benchmarking the sensitivity of the different analytical techniques for assessing stability in printable active layers, independent of complete device architectures.

14 SOLAR ENERGY↗
DOE OSTI2023

CHGNet as a pretrained universal neural network potential for charge-informed atomistic modelling

Large-scale simulations with complex electron interactions remain one of the greatest challenges for atomistic modelling. Although classical force fields often fail to describe the coupling between electronic states and ionic rearrangements, the more accurate ab initio molecular dynamics suffers from computational complexity that prevents long-time and large-scale simulations, which are essential to study technologically relevant phenomena. Here we present the Crystal Hamiltonian Graph Neural Network (CHGNet), a graph neural network-based machine-learning interatomic potential (MLIP) that models the universal potential energy surface. CHGNet is pretrained on the energies, forces, stresses and magnetic moments from the Materials Project Trajectory Dataset, which consists of over 10 years of density functional theory calculations of more than 1.5 million inorganic structures. The explicit inclusion of magnetic moments enables CHGNet to learn and accurately represent the orbital occupancy of electrons, enhancing its capability to describe both atomic and electronic degrees of freedom. We demonstrate several applications of CHGNet in solid-state materials, including charge-informed molecular dynamics in Li x MnO 2 , the finite temperature phase diagram for Li x FePO 4 and Li diffusion in garnet conductors. We highlight the significance of charge information for capturing appropriate chemistry and provide insights into ionic systems with additional electronic degrees of freedom that cannot be observed by previous MLIPs.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Multistage nucleation pathway in LiF molten salt mirrors the crystal–melt interface structure

Despite over a century of studies, fundamental questions remain about the processes governing crystal nucleation from melts or solutions. Research over the past three decades has presented mounting evidence for kinetic pathways of crystal nucleation that are more complex than envisioned by the simplest forms of classical theory. Such observations have been presented for colloidal and elemental systems with covalent and metallic bonding. Despite the technological and geochemical importance of molten salts, similar studies for these ionically bonded systems are currently lacking. Here we develop a machine learning interatomic potential for a model ionic system: LiF. The potential features quantum-level accuracy for both liquid and multiple solid polymorphs over wide temperature and pressure ranges and accurately reproduces experimentally measured properties. Thanks to the efficiency of the potential, which enables microsecond-scale molecular dynamics simulations, induction times for nucleation of LiF solids from their melts are computed over a range of undercoolings. With the aid of a set of robust local order parameters established here, the simulations reveal that homogeneous crystal nucleation in undercooled melts preferentially initiates from liquid regions showing slow dynamics and high bond orientational order simultaneously, and the second-shell order of both precritical nuclei and the surface of postcritical nuclei is dominated by hexagonal close packing and body-centered cubic local structure, even though the nucleus core is dominated by face-centered cubic structure corresponding to the stable rocksalt crystal structure. Finally, we establish a connection between the crystallization pathway and the equilibrium crystal-melt interface structure.

Applied Physical Sciences↗
DOE OSTI2023

Using limited neural networks to assess relative mechanistic influence on shock heating in granular solids

The rapid compaction of granular media results in localized heating that can induce chemical reactions, phase transformations, and melting. However, there are numerous mechanisms in play that can be dependent on a variety of microstructural features. Machine learning techniques such as neural networks offer a ubiquitous method to develop models for physical processes. Limiting what kind of microstructural information is used as an input and assessing normalized changes in network error, the relative importance of different mechanisms can be inferred. Here we utilize binned, initial density information as network inputs to predict local shock heating in a granular high explosive trained from large-scale molecular dynamics simulations. In this study, the spatial extent of the density field used in the network is altered to assess the importance and relevant length scales of the physical mechanisms in play, where different microstructural features result in different predictive capabilities.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Noncovalent Immobilization of Pentamethylcyclopentadienyl Iridium Complexes on Ordered Mesoporous Carbon for Electrocatalytic Water Oxidation

The attachment of molecular catalysts to conductive supports for the preparation of solid-state anodes is important for the development of devices for electrocatalytic water oxidation. The preparation and characterization of three molecular cyclopentadienyl iridium(III) complexes, Cp*Ir(1-pyrenyl(2-pyridyl)ethanolate-κO,κN)Cl (1) (Cp* = pentamethylcyclopentadienyl), Cp*Ir(diphenyl(2-pyridyl)methanolate-κO,κN)Cl (2), and [Cp*Ir(4-(1-pyrenyl)-2,2'-bipyridine)Cl]Cl (3), as precursors for electrochemical water oxidation catalysts, are reported. These complexes contain aromatic groups that can be attached via noncovalent π-stacking to ordered mesoporous carbon (OMC). The resulting iridium-based OMC materials (Ir-1, Ir-2, and Ir-3) were tested for electrocatalytic water oxidation leading to turnover frequencies (TOFs) of 0.9–1.6 s -1 at an overpotential of 300 mV under acidic conditions. The stability of the materials is demonstrated by electrochemical cycling and X-ray absorption spectroscopy analysis before and after catalysis. Theoretical studies on the interactions between the molecular complexes and the OMC support provide insight onto the noncovalent binding and are in agreement with the experimental loadings.

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

Rational design of heterogeneous single-site catalysts via surface organometallic chemistry

Single-site heterogeneous catalysts offer an attractive route to unite the molecular precision of homogeneous catalysis with the durability and practical advantages of solids. Surface organometallic chemistry (SOMC) provides a particularly powerful strategy for this purpose by grafting molecular precursors onto tailored surfaces and converting support functionalities into ligand environments for isolated metal centers. As a result, SOMC brings the language and logic of coordination chemistry to heterogeneous catalysis, where the support becomes an integral part of the active site coordination sphere. This Review surveys recent progress in the rational design of SOMC-derived single-site catalysts, with emphasis on synthetic routes, post synthetic transformations, and the deliberate tuning of catalytic behavior through metal-support interactions. Discussions are made on how support identity, hydroxyl topology, acidity, and redox activity shape the geometry, electronic structure, and oxidation state of supported metal sites, as well as how these factors determine activity, selectivity, and stability. We also examine a central limitation of these systems: despite their molecularly informed design, supported single sites often exist as structurally distributed ensembles rather than uniform species, particularly on amorphous supports. This site heterogeneity, along with catalyst dynamics under operating conditions, remains a major barrier to definitive structure-activity relationships. Therefore, emerging approaches that combine advanced characterization, first-principles modeling, ensemble kinetics, and machine learning to resolve active-site structure and guide catalyst development are highlighted. Together, these advances position SOMC as a versatile coordination chemistry framework for the predictive design of heterogeneous catalysts with well-defined molecularly tailored active sites.

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

Elucidating Lithium Transport Mechanisms in Disordered LiF from Machine-Learning Molecular Dynamics Simulations

Lithium fluoride (LiF) is a ubiquitous component of solid- and cathode–electrolyte interphases, yet its functional role remains unclear under the structural and chemical heterogeneity typical of cycling batteries. Here, we systematically quantify how structural disorder, off-stoichiometry, and strain govern Li-ion transport in LiF. Using a machine-learning potential to enable extensive molecular-dynamics sampling, we compare crystalline and amorphous LiF, Li 0.95 F, and LiF 0.95 , and evaluate the impact of small homogeneous deformations. Defect-free crystalline LiF is effectively ion-blocking at 300−500 K, whereas amorphization generates free-volume–assisted percolation pathways that facilitates Li-ion diffusion. At elevated temperatures, thermodynamically driven crystallization disrupts these pathways, leading to non-Arrhenius behavior. In crystalline phases, Li deficiency activates vacancy-mediated diffusion, while in amorphous LiF, transport is governed primarily by network connectivity. Strain is found to have only a marginal effect on Li mobility in both crystalline and amorphous structures.

Batteries↗
DOE OSTI2020

Mean force kinetic theory applied to self-diffusion in supercritical Lennard-Jones fluids

A new kinetic-theory-based calculation of the self-diffusion coefficient for dense supercritical Lennard-Jones fluids is presented. The mean force kinetic theory, which was recently developed for transport in dense plasmas, is applied for the calculation of diffusion in dense neutral fluids. The calculation only requires the pair distribution function, a quantity that is readily calculable from equilibrium statistical mechanics for many systems, including the Lennard-Jones fluid. The self-diffusion coefficients are compared with calculations from molecular dynamics simulations, and good agreement at high density is demonstrated, even in the vicinity of the solid–fluid coexistence line. A comparison of different kinetic models with molecular dynamics simulations demonstrates that the transport coefficients have important contributions due to particle interaction via a potential of mean force and local correlations, which increase the collision rate. The new calculations compare well to those from free-volume theory and overcome a limitation of this theory that prevents its use in systems that interact via long range monotonic potentials. It is expected that this approach will also apply to other systems, including neutral-plasma and neutral-electrolyte mixtures.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Influence of Charge Block Length on Conformation and Solution Behavior of Polyampholytes

In this paper, we investigate the effect of charge block length on polyampholyte chain conformation and phase behavior using small-angle X-ray scattering (SAXS) and implicit-solvent molecular simulations. To this end, we use solid phase peptide synthesis to precision-tailor a series of polyampholytes consisting of l-glutamic acid (E) and l-lysine (K) monomers arranged in alternating blocks from 2 to 16 monomers. We observe that the polyampholytes tend to phase separate as block size increases. With addition of NaCl, phase separated polyampholytes exhibit a salting-in effect dependent on charge block length. Fourier-transform infrared (FTIR) spectroscopy reveals the presence of intramolecular hydrogen bonds that are disrupted upon the addition of NaCl, implicating both electrostatic interactions and hydrogen bonding in the phase behavior. SAXS spectra at no-added salt conditions show minimal dependence of charge block length on the radius of gyration (R g ) for soluble polyampholytes, but local chain stiffening is found to be dependent on charge block length. With increasing NaCl, consistent with electrostatic screening, all polyampholytes expand and behave as neutral or swollen chains in good solvent conditions. Molecular simulations are qualitatively consistent with experiments. Implications for understanding intracellular condensates and material design are noted.

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

Computational Modeling of Battery Materials

This chapter provides a broad overview of the various computational modeling techniques used to gain fundamental insights into coupled electrochemical processes that occur in battery materials at electronic-to-mesoscopic scales. Furthermore, representative successes of these techniques in modeling electrodes, electrolytes, and electrode-electrolyte interfaces are highlighted to establish the current state-of-the-art in the field.

25 ENERGY STORAGE↗
DOE OSTI2021

Dynamics in Flexible Pillar[ n ]arenes Probed by Solid-State NMR

Pillar[n]arenes are supramolecular assemblies that can perform a range of technologically important molecular separations which are enabled by their molecular flexibility. Here, we probe dynamical behavior by performing a range of variable-temperature solid-state NMR experiments on microcrystalline perethylated pillar[n]arene (n = 5, 6) and the corresponding three pillar[6]arene xylene adducts in the 100–350 K range. This was achieved either by measuring site-selective motional averaged 13 C 1 H heteronuclear dipolar couplings and subsequently accessing order parameters or by determining 1 H and 13 C spin–lattice relaxation times and extracting correlation times based on dipolar and/or chemical shift anisotropy relaxation mechanisms. We demonstrate fast motional regimes at room temperature and highlight a significant difference in dynamics between the core of the pillar[n]arenes, the protruding flexible ethoxy groups, and the adsorbed xylene guest. Additionally, unexpected and sizable 13 C 1 H heteronuclear dipolar couplings for a quaternary carbon were observed for p-xylene adsorbed in pillar[6]arene only, indicating a strong host–guest interaction and establishing the p-xylene location inside the host, confirming structural refinements.

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

Tuning the melting point of selected ionic liquids through adjustment of the cation’s dipole moment.

In prior research with thermally robust salts [Cassity et al., PCCP, 2017, 19, 31560] it was noted that an increase in the dipole moment of the cation generally led to a decrease in the melting point. Molecular dynamics simulations of the liquid state revealed that an increased dipole moment reduces cation-cation repulsions through dipole-dipole alignment. This was believed to reduce the liquid phase enthalpy, which would tend to lower the melting point of the IL. Here, we also test this principle by replacing hydrogen atoms with fluorine atoms at selected positions within the cation. This allows us to alter the electrostatics of the cation without substantially affecting the sterics. Furthermore, the strength of the dipole moment can be controlled by choosing different positions within the cation for replacement. We studied variants of four different parent cations paired with bistriflimide and determined their melting points, and enthalpies and entropies of fusion through DSC experiments. The decreases in the melting point were determined to be enthalpically driven. We found that the dipole moment of the cation, as determined by quantum chemical calculations, is inversely correlated with the melting point of the given compound. Molecular dynamics simulations of the crystalline and solid states of two isomers showed differences in their enthalpies of fusion that closely matched those seen experimentally. Moreover, this reduction in the enthalpy of fusion was determined to be caused by an increase in the enthalpy of the crystalline state. We provide evidence that dipole-dipole interactions between cations leads to the formation of cationic domains in the crystalline state. These cationic associations partially block favourable cation–anion interactions, which are recovered upon melting. If, however, the dipole- dipole interactions between cations is too strong they have a tendency to form glasses. This study provides a design rule for lowering the melting point of structurally similar ILs by altering their dipole moment.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Solid Electrolyte Interphase Formation at the Ionic Liquid Electrolyte—Lithium-Metal Interface Using an Ab Initio Molecular Dynamics Approach

An interfacial study is performed using ab initio molecular dynamics (AIMD) simulations to elucidate the electrochemical interfacial stability of an ionic liquid (IL) as component of an electrolyte in contact with a Li-metal anode when an additive, vinylene carbonate (VC), is added to the electrolyte. The IL electrolyte is composed of 1-butyl-1-methylpyrrolidinium bis(fluoro-sulfonyl)-imide (PYR14 + FSI − ) with lithium bis(trifluoromethane-sulfonyl)imide salt (Li + TFSI − ) in a 8:2 concentration ratio. Reactions observed at the interface electrolyte/Li-metal during the AIMD simulations are verified by electronic structure calculations using density functional theory calculations considering a continuum solvent model. We found that the addition of VC to the IL prevents further reduction of the counter-ion TFSI − and retards reduction of FSI − in comparison to the case when the additive VC is absent. The PYR14 + reduction is never observed with or without the additive. The 1st, 2nd, and 3rd electron affinities of the PYR14 + are −0.08, −0.36 and −3.65 eV, indicating little to large repulsive behavior against electrons. The opposite is observed for the two negative ions, FSI − and TFSI − , which yield 1st electron affinities of 1.73 and 1.83 eV, respectively. Reactions obtained under open circuit conditions. Dissociation mechanisms are presented for all IL components reduced at the IL/Li-metal interface.

Electrochemistry↗