Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stability theory”

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 649 records · Page 36

Automatically adaptive stabilized finite elements and continuation analysis for compaction banding in geomaterials

Under compressive creep, viscoplastic solids experiencing internal mass transfer processes can accommodate singular cnoidal wave solutions as material instabilities at the stationary wave limit. These instabilities appear when the loading rate is significantly faster than the material's capacity to diffusive internal perturbations, leading to localized failure features (e.g., cracks and compaction bands). These cnoidal waves, generally found in fluids, have strong nonlinearities that produce periodic patterns. Due to the singular nature of the solutions, the applicability of the theory is currently limited. Additionally, practical simulation tools require proper regularization to overcome the challenges that singularity induces. We focus on the numerical treatment of the governing equation using a nonlinear approach building on a recent adaptive stabilized finite element method. This automatic refinement method provides an error estimate that drives mesh adaptivity, a crucial feature for the problem at hand. We compare the performance of this adaptive strategy against analytical and standard finite element solutions. We then investigate the sensitivity of the diffusivity ratio, the parameter controlling the process, and identify multiple possible solutions with several stress peaks. We also show the evolution of the spacing between peaks for all solutions as a function of that parameter.

42 ENGINEERING↗

Extended Lagrangian Born–Oppenheimer molecular dynamics for orbital-free density-functional theory and polarizable charge equilibration models

We report extended Lagrangian Born–Oppenheimer molecular dynamics (XL-BOMD) is formulated for orbital-free Hohenberg–Kohn density-functional theory and for charge equilibration and polarizable force-field models that can be derived from the same orbital-free framework. The purpose is to introduce the most recent features of orbital-based XL-BOMD to molecular dynamics simulations based on charge equilibration and polarizable force-field models. These features include a metric tensor generalization of the extended harmonic potential, preconditioners, and the ability to use only a single Coulomb summation to determine the fully equilibrated charges and the interatomic forces in each time step for the shadow Born–Oppenheimer potential energy surface. The orbital-free formulation has a charge-dependent, short-range energy term that is separate from long-range Coulomb interactions. This enables local parameterizations of the short-range energy term, while the long-range electrostatic interactions can be treated separately. The theory is illustrated for molecular dynamics simulations of an atomistic system described by a charge equilibration model with periodic boundary conditions. The system of linear equations that determines the equilibrated charges and the forces is diagonal, and only a single Ewald summation is needed in each time step. The simulations exhibit the same features in accuracy, convergence, and stability as are expected from orbital-based XL-BOMD.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Causal and stable first-order chiral hydrodynamics

We derive the set of inequalities that is necessary and sufficient for nonlinear causality and linear stability of first-order relativistic hydrodynamics with either a U ( 1 ) V conserved current or a U ( 1 ) A current with a chiral anomaly or both. Our results apply to generic hydrodynamic frames in which no relations among the transport parameters are imposed. Furthermore, our analysis yields, to the best of our knowledge, the first theory of viscous chiral hydrodynamics proven to be causal and stable. We find that causality demands the absence of vorticity-induced heat flux, forcing a departure from the thermodynamic frame in the chiral case. The inequalities for causality and stability define a hypervolume in the space of transport parameters, wherein each point corresponds to a consistent formulation. Notably, causality is determined by just three combinations of transport parameters. We present our results in a form amenable to numerical hydrodynamic simulations. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

Theoretical Investigation of the Oxygen Reduction Reaction over Platinum Catalysts Supported by Multi-Edged Vertically Aligned Carbon Nanofiber for Electrocatalyst Preparation

Vertically aligned carbon nanofibers (VACNFs) are promising supports for oxygen reduction reaction (ORR) electrocatalysts in fuel cells. Although experimentally these catalytic systems have shown great potential, there is lack of molecular understanding of the catalytic sites and reaction mechanisms. This work investigated the origin of the ORR reactivities of the platinum catalysts on multi-edged VACNFs (Pt/VACNF) using a multiscale modeling approach combining Density Functional Theory (DFT) and classical Molecular Dynamics (MD) simulations. Based on the ReaxFF potential, here all nanoscale Pt particles (Pt 55 , P 100 , and Pt 147 ) are stabilized by the open edges located axially along the VACNF walls. The calculated first-shell coordination numbers, CN, of surface Pt atoms are 6.63, 7.27, and 7.85, respectively, suggesting that the percentage of low coordination sites increases as the particle size decreases. The adsorption energies of OOH, O, and OH on Pt-55 were systematically probed using DFT calculations. These adsorption energies retain a linear correlation against the generalized coordination numbers ($\overline{CN}$). For Pt nanoparticles supported on VACNF, we found that the OOH and OH bind stronger than on Pt (111) by 0.14 and 0.17 eV, respectively, which can hinder the ORR activity with lower limiting potential than Pt (111). Our theoretical prediction is in good agreement with the linear sweeping voltammetry that revealed a left shift of the half-wave potential.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Cage Molecules Stabilize Lead Halide Perovskite Thin Films

The environmental stability of hybrid organic-inorganic perovskite (HOIP) materials needs to increase to enable their widespread adoption in thin-film solar and optoelectronic devices. Molecular additives have recently emerged as an effective strategy for regulating HOIP crystal growth and passivating defects. However, to date the choice of additives is largely limited to a dozen or so materials under the design philosophy that high crystallinity is a prerequisite for stable HOIP thin films. Here, in this study, we incorporate porous organic cages (POCs) as functional additives into perovskite thin films for the first time and investigate the HOIP-POC interaction via a combined experimental and computational approach. POCs are significantly larger than the small-molecule additives explored for HOIP synthesis to date but much smaller than polymeric sealants. Partially amorphized composites of MAPbI 3 (methylammonium lead iodide, HOIP) and RCC3 (an amine POC) form a network-like surface topography and lead to an increase in the optical bandgap from 1.60 to 1.63 eV. Further in situ optical imaging suggests that RCC3 can delay the MAPbI 3 film degradation onset up to 50x under heat and humidity stresses, showing promise for improving reliability in HOIP-based solar-cell and light-emitting applications. Furthermore, there is evidence of molecular interactions between RCC3 and MAPbI 3 , as fingerprinted by the suppressed N-H stretching mode in MA + from Fourier transform infrared (FTIR) spectra and density functional theory (DFT) simulations that suggest strong hydrogen bonding between MA + and RCC3. Given the diversity of POCs and HOIPs, our work opens a new avenue to stabilize HOIPs via tailored molecular interactions with functional organic materials.

36 MATERIALS SCIENCE↗

Accelerated Aging of a Polymer-Supported Molecular Amine Direct Air Capture Sorbent: Experiments and Theory

Amine-functionalized polymer sorbents are the preferred materials for CO2 extraction from air (DAC) due to their “high” CO2 adsorption capacity and selectivity at ultra-dilute conditions. The lifetime of adsorbents is a key property of equal importance as adsorption capacity, selectivity and kinetics with direct impact on the economics of commercial-scale operations. Therefore, a detailed investigation into the contributing factors (e.g. temperature, humidity) as well as the mechanism of amine degradation is vital for the development of amines with long-term stability. Herein, in situ attenuated total reflection (ATR) Fourier transform infrared spectroscopy (FTIR) is used to monitor the real time changes to amine-based sorbents under systematically varied conditions. In addition, IR spectra obtained from density functional theory molecular dynamics are used to obtain mechanistic insights pertaining to the degradation of amine groups in DAC sorbents.

Muldoon, Patrick↗

A Si-Substituted Spirobifluorene Hole-Transporting Material for Perovskite Solar Cells

Hybrid organic-inorganic perovskite solar cells (PSCs) have shown promise for next-generation photovoltaics. Here, this study presents a simple approach for enhancing the performance and stability of PSCs by substituting the center carbon atom of the common hole transport material (HTM) Spiro-MeOTAD with a silicon atom. This modification, termed Si-Spiro, results in an increased hole mobility. A density functional theory simulation indicates that the enhanced hole mobility is due to the structural change of Si-Spiro. Electron beam-induced current microscopy measurements indicate improved charged extraction at the Si-Spiro/perovskite interface. A power conversion efficiency of 22.5% is achieved in Si-Spiro-based PSCs, outperforming standard Spiro-MeOTAD. Additionally, Si-Spiro-based PSCs demonstrate enhanced stability, maintaining over 90% of performance over 120 h of one-sun operation. Depth-profiling X-ray photoelectron spectroscopy revealed that Si-Spiro effectively blocks metal ion migration, which contributes to its enhanced stability. The findings suggest that Si-Spiro could be promising a HTM for high-performing, stable PSCs.

14 SOLAR ENERGY↗

Hybrid simulations of sub-cyclotron compressional and global Alfvén eigenmode stability in spherical tokamaks

A comprehensive numerical study has been conducted in order to investigate the stability of beam-driven, sub-cyclotron-frequency compressional Alfvén eigenmodes (CAEs) and global Alfvén eigenmodes (GAEs) in low-aspect-ratio plasmas for a wide range of beam parameters. The presence of CAEs and GAEs has previously been linked to anomalous electron temperature profile flattening at high beam powers in NSTX experiments, prompting a further examination of the conditions necessary for their excitation. Linear simulations have been performed with the hybrid MHD–kinetic initial value code HYM in order to capture the general Doppler-shifted cyclotron resonance that drives the modes. Additionally, three distinct types of modes were found in the simulations—co-CAEs, cntr-GAEs, and co-GAEs—with differing spectral and stability properties. The simulations revealed that unstable GAEs are more ubiquitous than unstable CAEs, which is consistent with experimental observations, as they are excited at lower beam energies and generally have larger growth rates. Local analytic theory is used to explain key features of the simulation results, including the preferential excitation of different modes based on beam injection geometry and the growth rate dependence on the beam injection velocity, critical velocity, and degree of velocity space anisotropy. The background damping rate is inferred from the simulations and analytically estimated for relevant sources absent from the simulation model, indicating that co-CAEs are closer to marginal stability than modes driven by the cyclotron resonances.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Anomalous stability of non-van der Waals bonded B4C nanosheets through surface reconstruction

Boron carbide (B4C) has been well studied both theoretically and experimentally in its bulk form due to its exceptional hardness and use as a high-temperature thermoelectric. However, the properties of its two-dimensional nanosheets are not well established. In this paper, using van der Waals-corrected density-functional theory simulations, we show that bulk B4C can be cleaved along different directions to form B4C nanosheets with low formation energies. We find that there is minimal dependence of formation energies on cleavage planes and surface terminations, even though the bulk is not van der Waals layered. This anomalous stability of B4C nanosheets is found to be a result of surface reconstructions that are unique to B-rich systems. While the density of states of the bulk B4C indicate that it is a semiconductor, the B4C nanosheets are found to be predominantly metallic. We attribute this metallic behavior to a redistribution of charges on the surface bonds of the films. The Seebeck coefficients of the B4C films remain comparable to those of the bulk and are nearly constant as a function of temperature. Our results provide guidance for experimental synthesis efforts and future application of B4C nanosheets in nanoelectronic and thermoelectric applications.

Physics↗

Decomposition pathways and mitigation strategies for highly-stable hydroxyphenazine flow battery anolytes

Aqueous organic redox flow batteries are a promising technology for large-scale energy storage. The stability of the redox active organic molecules is increasingly being recognized as one of the major hurdles. Upon extended flow battery cycling, 7,8-dihydroxyphenazine-2 sulfonic acid (DHPS) undergoes desulfonation and reduction of a phenolic C–O bond to yield a mixture of 7/8-hydroxyphenazine-2-sulfonic acid, as well as hydrogenation of the aromatic ring system. Density functional theory (DFT) analysis of the charged DHPS, its ring-hydrogenated products, and variably substituted hydroxy phenazines has led to the development of a series of dihydroxylated phenazine isomers which provide insight into the effects of substitution pattern on solubility and stability. Seven dihydroxyphenazine (DHP) isomers were synthesized and their solubilities, electrochemical properties, flow battery cycling performance, and degradation pathways have been investigated. Based on theoretical and experimental results, hydroxyl substitution at the 1, 4, 6 and 9 positions yields highly stable derivatives, while substitution at the 2, 3, 7, and 8 positions results in unstable derivatives. Flow cells of 1,4- and 1,6-DHPs coupled with ferro/ferricyanide achieved high stabilities, with temporal capacity loss of 0.029 and 0.031% per day, respectively. Decomposition of 1,8- and 2,7-DHPs were found to arise from irreversible hydrogen rearrangement (tautomerization), yielding redox-inactive species. These results provide a detailed understanding of decomposition pathways and mitigation strategies for phenazine-based anolytes, and can provide general design guidelines for the development of stable redox-active organics.

25 ENERGY STORAGE↗

Metal oxide candidates for thermochemical water splitting obtained with a generative diffusion model

Generative diffusion models (DMs) for inorganic crystalline materials are being actively investigated for their potential to expand the chemical and structural design spaces for known functional materials. Generative candidates are particularly useful for applications where few functional, let alone commercially viable, materials currently exist, such as metal oxides for thermochemical water-splitting, which have strict requirements for defect thermodynamics and host stability. Here, we critically examine generated metal oxides from the M ATTER G EN DM conditioned on select chemical systems for thermochemical water splitting applications. Perhaps most notably, we find that M ATTER G EN predicts a novel, thermodynamically stable, quinary metal oxide, Ba 2 SrInFeO 6 , although this compound represents an ordered and layered substitution within the same A 3 B 2 O 6 structural prototype as its two ternary end members. Detailed density functional theory calculations and spin configuration sampling for this material and its possible decomposition products—beyond what existed in M ATTER G EN training data—are required to quantitatively validate hull energy predictions and conclusions of stability. Furthermore, the material exhibits oxygen defect formation energies appropriate for thermochemical water splitting, warranting targeted investigation in an experimental validation campaign, along with other future M ATTER G EN candidates in this application space.

36 MATERIALS SCIENCE↗

High Capacity Adsorption—Dominated Potassium and Sodium Ion Storage in Activated Crumpled Graphene

In this study, structurally and chemically defective activated-crumbled graphene (A-CG) is employed to achieve unique synergy of large reversible potassium (K) and sodium (Na) ion storage capacity with fast charging and extended cyclability. A-CG synthesis consists of low temperature spraying of graphene oxide slurry, followed by partial reduction annealing and air activation. For K storage, the reversible capacities are 340 mAh g -1 at 0.04 A g -1 , 261 mAh g -1 at 0.5 A g -1 , and 210 mAh g -1 at 2 A g -1 . For Na storage, the reversible capacities are 280 mAh g -1 at 0.04 A g -1 , 191 mAh g -1 at 0.5 A g -1 , and 151 mAh g -1 at 2 A g -1 . A-CG shows a stable intermediate rate (0.5 Ag -1 ) cycling with both K and Na, with minimal fade after 2800 and 8000 cycles. These are among the most favorable capacity—rate capability—cyclability combinations recorded for potassium-ion battery and sodium-ion battery carbons. Electroanalytical studies (cyclic voltammetry, galvanostatic intermittent titration technique, b -value) and density functional theory (DFT) reveal that enhanced electrochemical performance originates from ion adsorption at various defects, such as Stone–Wales defects. Moreover, DFT highlights enhanced thermodynamic stability of A-CG with adsorbed K versus with adsorbed Na, explaining the unexpected higher reversible capacity with the former.

25 ENERGY STORAGE↗

Precisely doping the surface of tin-based electrocatalysts for improved CO 2 conversion to liquid chemicals

Doping tin catalysts with sulfur can improve the electrochemical CO 2 conversion into formate/formic acid, but the lack of composition-dependent activity trends hinders further catalyst development. Here, we precisely controlled the composition of sulfur-doped Sn catalysts to show that sulfur doping only improves CO 2 conversion over a very narrow composition range, achieving maximum activity at 1.4 at% S. In situ Raman spectroscopy indicted working catalysts were in a primarily metallic state (e.g. S-Sn), and we achieved some of the highest reported partial current densities in both H-cell and full-cell electrolyzer configurations. Density Functional Theory calculations predicted S atoms preferentially occupied the catalyst surface and improved CO 2 reduction by localizing charge density at the catalyst/intermediate interface, which stabilized the *OCOH intermediate and lowered the CO 2 conversion thermodynamic barrier. Our work quantifies the composition-dependent influence of S dopants on Sn-based CO 2 reduction catalysts and provides a pathway for maximizing their CO 2 conversion activity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Surface structures of magnetostrictive D0 3 -Fe 3 Ga(0 0 1)

First-principles total energy calculations and scanning tunneling microscopy experiments were performed to study the surface reconstruction of the magnetostrictive Fe 3 Ga alloy. The inverse magnetostrictive behavior was evaluated in the bulk by compressing and stretching its lattice parameter, showing an increase in magnetic moments as strain increases. Surface analysis demonstrates two thermodynamically stable surfaces, the (1 × 1) and (3 × 1). The (1 × 1) is an ideal FeGa terminated surface, whereas the (3 × 1) is also FeGa terminated but it has a first-layer Fe atom substituted by a Ga atom every three unit-cells, forming a row-like surface structure. Tersoff–Hamann scanning tunneling microscopy simulations were obtained and compared with experimental results. We found good agreement between theory and experiment, in which the distance between rows is ~12.3 Å. Theoretical findings suggest that the substrate-induced strain may increase the stability of the (3 × 1) reconstruction. Analysis of the magnetic moments in the reconstructions showed that their behavior is affected by a surface effect, as well as by the inverse magnetostriction of the structure. A good understanding of the FeGa surface reconstructions is an important step towards further improvements in magnetic storage devices and sensors.

36 MATERIALS SCIENCE↗

First-principles investigation of lanthanides diffusion in HCP zirconium via vacancy-mediated transport

The diffusion of lanthanide fission products plays an important role in the growth of the fuel-cladding chemical interaction (FCCI) region in metallic fuels. The use of a Zr interdiffusion barrier may mitigate the transport of lanthanides from the fuel to the cladding, but the efficacy of such a liner is not yet known. In this paper, the stability and vacancy-mediated diffusion of La, Ce, Pr, and Nd in hexagonal close-packed (HCP) Zr is investigated via density functional theory (DFT) calculations and self-consistent mean field (SCMF) analysis. DFT is used to calculate the formation, binding, and migration energies of vacancies and vacancy-solute pairs. The DFT energetics are used in the KineCluE code to calculate the Onsager transport coefficients. La is found to be the fastest diffusing species in HCP Zr and experiences an almost isotropic diffusion behavior. The other three species (Ce, Pr, and Nd) demonstrate anisotropic diffusion where the diffusion in the basal planes is significantly faster than that along the c-axis. Further, the calculated lanthanide diffusivities in HCP Zr are fitted to an Arrhenius relation and the activation energies and prefactors are reported for the first time. Furthermore, the vacancy drag and the segregation tendencies were analyzed using the calculated off-diagonal transport coefficients. According to our vacancy-mediated diffusion model, lanthanides are expected to be enriched at vacancy sinks at low temperatures, while at high temperatures, lanthanides are depleted at sinks and will preferably diffuse into the bulk. The enrichment/depletion transition temperature depends on the diffusion direction (basal or axial) and hence will be controlled by the grain texture and orientation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Structure, Bonding, and Vibrational Dynamics of a Triamine High Energy Density Material under Pressure

High energy density materials have complex intermolecular interactions which influence their stability and performance. We used a combination of synchrotron X-ray diffraction, synchrotron infrared spectroscopy, and Raman vibrational spectroscopy, supplemented by density functional theory calculations, to probe pressure-induced changes in structure and intermolecular interactions of 1H,4'H-[3,3'-bis(1,2,4-triazole)]-4',5,5'-triamine as a model high energy density material up to 40 GPa. We find that compression of the triamine is accompanied by increased intermolecular interactions that give rise to an interesting evolution of the structure and bonding with pressure. Analysis of the equation of state determined from the X-ray diffraction indicates a change in compression mechanism near 19 GPa consistent with changes in vibrational spectra that provide evidence for a structural rearrangement associated with changes in hydrogen bonding near that pressure. As a result, the overall compressional behavior calculated theoretically agrees with that observed experimentally though differences are found that indicate the need for improved treatment of the intermolecular interactions including hydrogen bonding under pressure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Experimental and Computational Evaluation of Nicotinamide Cofactor Biomimetics

Oxidoreductase enzymes are widely used biocatalysts due to their high enantioselectivity and broad substrate compatibility in useful transformations. Many oxidoreductases require nicotinamide cofactors (i.e., NAD(P)H). To replace this costly natural cofactor, synthetic nicotinamide cofactor biomimetics (NCBs) offer different shapes, binding affinities, and reducing potentials that exceed the capabilities of wild-type NAD(P)H. However, the ill-defined structure–activity relationships (SARs) of various NCBs slow rationally guided innovation, such as customized reducing potentials. Here, we dissect two essential elements of NCB design, holding the nicotinamide invariant. First, the linker length between the nicotinamide and an unconjugated aromatic ring uncovered unexpected benefits to redox activity for two or three carbon linkers. Second, substitution on this unconjugated aryl group (Ring 2) might not be expected to affect activity. However, SAR trends demonstrate substantial benefits to reductive potential conferred by electron-donating functionalities on Ring 2. Furthermore, catalysis by two enzymes demonstrates enzyme-dependent tolerance or sensitivity to the NCB structures. Density functional theory (DFT) and computational modeling provide a theoretical framework to understand and build upon these observations. Ring 2 reaches up to the nicotinamide to stabilize its positive charge after oxidation through π–π stacking and charge transfer. Thus, the systematic examination of NCB’s stability, electrochemical redox potentials, and kinetics uncovers trends for the improved design of NCBs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electromagnetic signatures of a chiral quantum spin liquid

Quantum spin liquids (QSL) have emerged as a captivating subject within interacting spin systems that exhibit no magnetic ordering even at the lowest temperature accessible experimentally. However, definitive experimental evidence remains elusive. In light of the recent surge in theoretical and experimental interest in the half-filled Hubbard model on a triangular lattice, which offers the potential for stabilizing a chiral QSL, we investigate the electromagnetic signatures of this phase to facilitate experimental detection. Utilizing a combination of parton mean-field theory and unbiased density-matrix renormalization group calculations, we systematically examine the electrical charge and orbital electrical current associated with a spinon excitation in the chiral QSL. Additionally, we calculate the longitudinal and transverse optical conductivities below the Mott gap. Furthermore, employing quantum field theory analysis, we unravel the connection between spinon excitations and emergent as well as physical gauge fields. Our results demonstrate that the chiral QSL phase exhibits a distinct electromagnetic response, even within a Mott insulator regime. This finding holds great potential for enabling the experimental detection of this long-sought-after phase.

36 MATERIALS SCIENCE↗