Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Energy Density”

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Understanding the Electronic Behavior of f-Element Intermetallics

Due to their partially filled f-orbitals, lanthanide and actinide intermetallic materials exhibit intriguing electronic, structural, and magnetic properties. These materials provide a unique platform to explore exotic magnetic and quantum phenomena that arise from their strongly correlated f-electrons, offering insights into f-electron behavior. Incorporating light elements such as boron and carbon can significantly impact the electronic, chemical, and physical properties of these materials, with varying influence between lanthanides and actinides due to the more localized f-electrons in lanthanides. In this work, lanthanide-based intermetallic single crystals containing light atoms were synthesized using the molten metal flux growth method. In this method, one or more low-melting metals such as aluminum, gallium, tin, and bismuth are used in excess as the reaction medium. This solution-state approach allows for lower reaction temperatures compared to solid-state reactions and enables the isolation of kinetic products rather than thermodynamically stable compounds. The synthesized crystals were analyzed using an extensive array of analytical and computational characterization techniques, including single crystal X-ray diffraction, scanning electron microscopy, energy-dispersive spectroscopy, magnetization measurements, and density functional theory calculations. These analyses provided detailed insights into the structural and magnetic properties of the synthesized materials.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗
DOE OSTI2025

Computational methods based on density functional theory for reactions and processes involving electronic spin (Final Technical Report)

This award supports one post-doctoral researcher for 1.5 years. Publications that acknowledge this grant: Refs. 1–14. Refs. 3,10,12 assess current methodology for the evaluation of magnetic exchange couplings in transition metal complexes. In particular, Ref. 10 validates the use of an approximate (non-iterative) Green’s function approach for the calculation of magnetic exchange couplings and will be the foundation for Thrust 2 in this proposal. Refs. 3 and 12 focus on widely used density functional approaches based on the standard energy differences methodology for the particular case of oxo-bridged Fe(III) complexes. Refs. 2,4–7,11 apply current methodologies to problems of practical interest in molecular magnetism. Ref. 13 presents a methodology to explicitly simulate the dynamics of open quantum systems within density functional theory (DFT) calculations based on the Liouville-von Neumann equation of motion for quantum systems driven out-of-equilibrium. Ref. 8 uses non-collinear spin DFT to explain the mechanical behavior of magnetic mono-atomic Pt wires produced in break-junction experiments in the presence of a magnetic field.

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

Alkali triel chalcogenide nanocrystals: a molecular reactivity approach to ternary phase selectivity

Alkali-metal-based materials are promising building blocks for energy conversion and storage technologies. Here, we use a molecular reactivity-based solution-phase approach to selectively synthesize multiple phases and specific polymorphs of lithium- and sodium-containing triel chalcogenide nanocrystals, LiTrCh 2 , NaTrCh 2 , and NaTr 3 Ch 5 , where Tr = Ga and In and Ch = S, Se, and Te. Analogous to the case of binary II–VI and III–V tetrahedral semiconductors, where the two commonly isolated zinc blende and wurtzite polymorphs are separated by only 1–50 meV f.u. −1 , we find that LiTrCh 2 nanocrystals easily adopt tetragonal (chalcopyrite) and orthorhombic polymorphs separated by only 2.7–6.2 meV f.u. −1 Because of this small energy difference, soft colloidal synthesis succeeds in accessing either one of these polymorphs, depending on the specific dichalcogenide precursor used. Highly reactive diethyl diselenide favors the thermodynamically more stable tetragonal I$\bar{4}2$d phase, whereas mildly reactive diphenyl diselenide favors the kinetic, metastable orthorhombic Pna2 1 phase. Density functional theory calculations confirm the relative energies among multiple LiTrCh 2 polymorphs and also model the observed powder X-ray diffraction pattern of a new C2 NaIn 3 Te 5 phase. 7 Li, 69 Ga, and 77 Se solid-state NMR spectra are consistent with phase-pure ternary LiGaSe 2 nanocrystals. A majority of the nanocrystal compositions are visible-light emitters. This work opens the door to new Li/Na-based ternary triel chalcogenide nanostructures for energy storage and conversion applications.

Pavel, Md Riad Sarkar [Iowa State Univ., Ames, IA ↗
DOE OSTI2026

An Exploration of the Equation of State Dependence of Core-collapse Supernova Explosion Outcomes and Signatures

We explore, using a state-of-the-art simulation code in 3D and late-enough times to witness final observables, the dependence of core-collapse supernova explosions on the nuclear equation of state (EOS). Going beyond questions of explodability, we compare final explosion energies, nucleosynthetic yields, recoil kicks, and gravitational-wave and neutrino signatures using the SFHo and DD2 nuclear EOSs for a 9 M ⊙ /solar-metallicity progenitor star. The DD2 EOS is stiffer and has a lower effective nucleon mass. The result is a more extended protoneutron star (PNS) and lower central densities. As a consequence, the mean neutrino energies, final explosion energy, and recoil kick speed are lower. Moreover, the evolution of PNS convection differs between the two EOS models in significant ways. This translates in part into interestingly altered neutrino “light” curves and noticeably altered gravitational-wave signal strengths and frequency characteristics that may be diagnostic. The faster exploding model (SFHo) yields slightly more neutron-rich ejecta and more species with atomic weights between 60 and 90 and a weak r-process. However, this is merely a preliminary study. The next step is a more comprehensive and multiprogenitor set of 3D supernova simulations for various EOSs to late times when the observables have asymptoted. Such a future investigation will have a direct bearing on the neutron star and black hole birth mass functions and the quest toward a fully quantitative theory of supernova observables.

Rusakov, Aleksandr [Princeton University, NJ (Unit↗
DOE OSTI2025

Defect-Mediated Diffusion Pathways in Spodumene Accelerate Lithium Transport

Lithium extraction from naturally occurring α-spodumene is hindered by poor lithium diffusivity, necessitating high-temperature phase transformation to a low-density β polymorph. Although β spodumene exhibits up to 5 orders of magnitude higher lithium-ion diffusivity, both phases have diffusion activation energies between 0.8 and 1 eV, indicating that polymorph density is not the controlling factor over diffusivity. We show that aluminum vacancies facilitate lithium-ion diffusion in α-spodumene by reducing the migration barrier from 2.4 to 0.9 eV. Bond valence site energy and nudged elastic band calculations show a new lithium local minimum site which promotes a one-dimensional percolation network by reducing the lithium intersite distance from 4.5 Å to 2.9 Å. However, aluminum vacancies are energetically unfavorable to percolate through the whole structure, resulting in very low net lithium diffusivity and highlighting the critical role of nonstoichiometric defects in facilitating lithium transport in rigid aluminosilicate structures.

Chemical structure↗
DOE OSTI2024

Medium-density amorphous ice unveils shear rate as a new dimension in water’s phase diagram

Recent experiments revealed a new amorphous ice phase, medium-density amorphous ice (MDA), formed by ball-milling ice I h at 77 K [Rosu-Finsen et al., Science 379, 474–478 (2023)]. MDA has density between that of low-density amorphous (LDA) and high-density amorphous (HDA) ices, adding to the complexity of water’s phase diagram, known for its glass polyamorphism and two-state thermodynamics. The nature of MDA and its relation to other amorphous ices and liquid water remain unsolved. Here, we use molecular simulations under controlled pressure and shear rate at 77 K to produce and investigate MDA. Here. we find that MDA formed at constant shear rate is a steady-state nonequilibrium shear-driven amorphous ice (SDA), that can be produced by shearing ice I h , LDA, or HDA. Our results suggest that MDA could be obtained by ball-milling water glasses without crystallization interference. Increasing the shear rate at ambient pressure produces SDAs with densities ranging from LDA to HDA, revealing shear rate as a new thermodynamic variable in the nonequilibrium phase diagram of water. Indeed, shearing provides access to amorphous states inaccessible by controlling pressure and temperature alone. SDAs produced with shearing rates as high as 10 6 s −1 sample the same region of the potential energy landscape than hyperquenched glasses with identical density, pressure, and temperature. Intriguingly, SDAs obtained by shearing at ~10 8 s −1 have density, enthalpy, and structure indistinguishable from those of water “instantaneously” quenched from room temperature to 77 K over 10 ps, making them good approximants for the “true glass” of ambient liquid water.

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

pH-Dependent Vibrational Dynamics Drives Excited-State Quenching in the Phycobiliprotein Complex PC645

Phycocyanin 645 (PC645) is a closed-form lightharvesting complex found in the lumen of the photosynthetic membrane of cryptophyte algae. These peripheral antenna complexes contain bilin chromophores that absorb sunlight and transfer excitation energy to the core antenna complexes embedded in the thylakoid membrane. The location of cryptophyte antenna complex on the luminal side of the membrane is unusual. During photosynthetic activity, the pH of the lumen drops, by up to two pH units. There is little known about how this pH-change affects the light-harvesting complexes. In this study, we report multiscale simulations using a computationally efficient density functional tight-binding framework to investigate the spectroscopy and excitation energy transfer in the PC645 complex. Complementary experiments were conducted using both steady-state and time-resolved spectroscopic measurements at low, neutral, and high pH values. Our study shows that (de)protonation of specific bilin pigments, namely, the mesobiliverdins (MBVs), modulates the excitation energies, excitonic couplings, and spectral densities. These changes cause excitation transfer rates to increase by up to a factor of two to three, leading to pH-dependent energy transfer pathways in the complex. Using this model, we calculated the pH-dependent fluorescence quantum yield of the system, obtaining quantitative agreement with the experimental results. These computational simulations, supported by experiments, identify MBVs as a more prominent excitation sink than previously realized, and that this role is tuned by pH.

Maity, Sayan [Constructor Univ., Bremen (Germany);↗
DOE OSTI2025

Urban heat islands can influence the wind energy resource during heatwaves

Urban wind energy is critical for sustainable electricity generation in cities. However, little research has explored how the urban heat island (UHI) effect influences wind energy, particularly in heatwaves when energy demand surges. In this study, we examine wind energy distribution in the Boston–Providence metropolitan area during heatwaves, using Weather Research and Forecasting (WRF) model integrated with Building Energy Parameterization/Building Energy Model (BEP/BEM). Two scenarios, a realistic case and a hypothetical case without urban warmth, were compared to isolate UHI impacts. Results reveal that UHI induces a "wind energy loss zone" in this urban area, reducing wind power density (WPD) by 20–30 W/m 2 at 50–100 m, while suburban/rural areas exhibit a "wind energy gain zone," with WPD increases up to 40 W/m 2 at 150–200 m. These losses diminish with distance from urban centers and become negligible beyond main urban and suburban sprawl. Heatwave expands the urban "loss zone", while amplifying wind energy gains in suburban/rural areas, driven by stronger thermal gradients and weakened background winds that intensify air convergence in urban and urban-rural circulations, thereby exacerbating urban wind energy losses by 15–20 %. An analysis of 235 wind farms using turbine power curves reveals that built areas dependent on stand-alone or off-grid turbines face significant energy deficits during a heatwave. Wind energy drops by up to 25 %, while cooling-related building energy demand rises 30–40 % during a heatwave. These findings underscore the need for strategic urban wind energy planning to ensure reliable power during extreme heat.

Energy - Wind↗
DOE OSTI2025

Likelihood-Based Particle Identification in the Short-Baseline Near Detector

Accurate particle identification is crucial in any high-energy physics experiment, allowing scientists to understand the unique interactions and mechanisms at play in a detector. In this project, I develop and study a new particle identification (PID) algorithm for the Short-Baseline Near Detector, a likelihood-based approach, different from out current $\chi^2$ method. A likelihood estimation offers a more physically motivated strategy for PID. The distribution random energy losses of charged particles traveling through a medium are described by the Vavilov probability density function. By using this model, we can account for random energy losses and construct likelihood functions specific to each particle type, potentially enabling a more accurate method for PID.

Vanderwaal, Sophia [U. Alabama, Huntsville] (ORCID↗
DOE OSTI2024

One-dimensional mapping of femtosecond laser filaments using coherent microwave scattering

Here, this paper reports on the use of coherent microwave scattering (CMS) for spatially resolved electron number density measurements of elongated plasma structures induced at mid-IR femtosecond filamentation in air. The presented studies comprise one-dimensional mapping of laser filaments induced via 3.9 µm, 127.3 fs laser pulses at output energies up to 15 mJ. The axial electron number density was measured to be invariant (about 2 × 10 15 c⁢m −3 ) along the entire filament length and for all tested laser pulse energies 5–15 mJ, and the corresponding laser intensity in the middle portion of the filament was estimated to be nearly constant for 5–15 mJ pulse energies (about 30−40 TW/c⁢m 2 ). These findings support that intensity clamping conditions were achieved in the experiments. The proposed approach enables capabilities that are currently unavailable to perform absolute and longitudinally resolved measurements of electron number density in laser filaments and to precisely characterize conditions associated with self-focusing and intensity clamping.

plasma diagnostic techniques↗
DOE OSTI2026

Custom-trained Machine-learning Interatomic Potentials: ZnCl2 Aqueous Solution

This dataset was generated using an iterative active-learning strategy implemented in the ArcaNN software package (https://github.com/arcann-chem/arcann_training) to train machine-learning interatomic potentials for aqueous ZnCl2 solutions. Each active-learning cycle consisted of three stages: training, exploration, and labeling. The initial training set combined configurations generated in this work from enhanced-sampling ab initio molecular dynamics simulations with configurations from a previously reported neural-network-potential study of aqueous ZnCl2. The enhanced-sampling ab initio molecular dynamics simulations involved Zn–Cl separation and the chloride coordination number around Zn²? as collective variables. These configurations served as the seed dataset. Subsequent active-learning cycles expanded the training set by identifying and labeling configurations that were poorly represented by the current models, thereby improving coverage of ion-association states and changes in local coordination and charge-state environments relevant to the solution free-energy landscape. For all selected configurations, single-point calculations of the total energies and atomic forces were performed within density functional theory using the CP2K Quickstep module. Reference calculations employed the revPBE-D3 and r2SCAN exchange-correlation functionals. Motivated by recent work on aqueous Zn²?, the main revPBE calculations omitted D3 dispersion contributions involving Zn²?, while retaining the D3 correction for water and chloride. For comparison, fully dispersion-corrected revPBE-D3 reference calculations were also performed, with D3 applied to all species, including Zn²?. Valence electrons were treated explicitly, while core electrons were represented using norm-conserving Goedecker–Teter–Hutter pseudopotentials. The wave functions were expanded using the mixed Gaussian-and-plane-wave scheme with TZV2P-MOLOPT basis sets for all elements and a 600 Ry auxiliary plane-wave cutoff for the electron density. Self-consistent-field convergence was accelerated using the orbital-transformation and Direct Inversion in the Iterative Subspace algorithms, with a convergence threshold of 10?6. All single-point calculations were performed in periodic orthorhombic cells. The CELL_REF keyword in CP2K was used to define a fixed reference cell with a box length of 25 Å. This treatment ensured a consistent reference for configurations extracted from NpT trajectories with fluctuating cell dimensions. The resulting DFT energies and atomic forces constitute the ground-truth labels used to train the MLIPs. The resulting MLIP was trained for aqueous ZnCl2 solutions spanning concentrations from 0 to 30 molal and a broad pH range, from strongly acidic to strongly basic conditions. Representative examples of configurations included in the MLIP training dataset are provided below. These include 1) Representative configurations from the dataset labeled at the revPBE-D3 level, with D3 dispersion interactions involving Zn2+ excluded (revPBE-wo-D3). 2) Representative configurations from the dataset labeled at the fully dispersion-corrected revPBE-D3 level, with D3 interactions applied to all species, including Zn2+ (revPBE-D3). 3) Representative configurations from the dataset labeled at the r2SCAN level of theory (r2SCAN).

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗
DOE OSTI2024

Computational descriptor for electrochemical currents of carbon dioxide reduction on Cu facets

Computation screening is crucial for designing efficient electrochemical catalysts for carbon dioxide (CO 2 R) reduction that produce valuable hydrocarbons and oxygenates. In this work, leveraging density functional theory calculations for the CO adsorption energy ΔE CO on seventeen Cu terminations, we discover a strong linear correlation between ΔE CO and the experimentally measured CO 2 R electrochemical currents (ACS Catal. 2022, 12, 11, 6578–6588). Examining ab initio thermodynamics of early critical intermediates CO*, COH*, and CHO*, we find that CO* → CHO* is the thermodynamically controlling step. Beyond the general CO adsorption energy that only shows a linear trend with CO 2 R activity, we show that the reaction free energy of CO* → CHO* is the descriptor for the overall CO 2 R activity for Cu facets, as it displays a volcano relationship with the experimental current. Importantly, we show that high step and kink density of the Cu terminations not only enhances CO adsorption strength but also modulates the CO* → CHO* pathway, as respectively exemplified in the (941) and (741) facets. In addition, we explain that the high activity of (741) is due to its relatively low hydrogen evolution reaction activity compared with the other Cu surfaces.

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

Jordan–Wigner Transformation for the Description of Strong Correlation in Fermionic Systems

Seniority is a useful way of organizing Hilbert space for strongly correlated systems. The exact zero-seniority wave function, doubly occupied configuration interaction (DOCI), provides accurate results (given the right orbitals) for many strongly correlated electronic systems but has a combinatorial computational cost. In many cases, pair coupled cluster doubles provide a polynomial-cost approximation that closely reproduces the energies of DOCI, but it breaks down in some cases and, as shown herein, it does not provide particularly good density matrices. In this article, we demonstrate that by using the Jordan–Wigner transformation to turn the seniority zero problem back into a Fermionic one, we can provide mean-field variational results of DOCI quality for the Hubbard model and a few small molecular dissociation examples, with polynomial cost, both for the energies and for density matrices, all while being protected from collapse. This success is rooted in the proof we provide, showing that the Hartree–Fock wave function on the Jordan–Wigner-transformed Hamiltonian transforms back to variational coupled cluster doubles in the seniority zero representation, but restricted to have determinant rather than permanent amplitude coefficients, without compromising its overall accuracy.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Energy Barriers for H 2 Uptake and Release by Defects in Boron Nitride

Density functional theory simulations (DFT) have been used in combination with the climbing image nudged elastic band method (CI-NEB) to determine the energy barriers required to dissociate H 2 at four types of h-BN defects and evaluate the potential for reacted H species to hop to adjacent available sites. The defects were selected based on previous work in which ab initio thermodynamics was used to identify specific defects with optimal thermodynamics for reversible H 2 storage, that is, the defects would dissociate molecular H 2 while not binding H species too strongly so that hydrogen release is favorable. The four defects investigated are the N and B monovacancies, which only offer either B or N defective sites respectively, as well as larger cluster vacancy defects named 3V(1B2N) and hexagonal 6V(3B3N), which present a mixture of B and N defect sites available for H binding. While homogeneous H 2 dissociation by B and N monovacancies have energy barriers larger than 0.4 eV, we found that mixed BN terminated defects have energy barrier lower than 0.4 eV for H 2 dissociation. In particular, H 2 dissociation by the hexagonal 6V(3B3N) defect, exclusively made of Frustrated Lewis Pairs (FLPs), is almost barrierless. The energy barrier for subsequent H hopping to available adjacent B or N sites of the defect were also calculated. It was found that the hopping of H species generally involves high energy barriers ranging from 0.26 eV to 2.57 eV. In conclusion, this work completes previous thermodynamic theoretical investigations and provides a more comprehensive picture of H 2 reaction with defective h-BN for potential H 2 storage and recovery.

Sassi, Michel J. [Pacific Northwest National Labor↗
DOE OSTI2026

Effects of heteroatom doping on hydrogen uptake in tungsten oxide

Redox-active transition metal oxides (TMOs) that can undergo proton-insertion coupled electron transfer (PICET) are promising candidates for catalyzing molecular conversion reactions which require the transfer of hydrogen atoms (or the thermochemical equivalent, H + , e − ). Herein, we studied the effects of isovalent (Mo 6+ ) and aliovalent (V 5+ and Nb 5+ ) heteroatom doping on the electrochemical PICET behavior of monoclinic tungsten oxide (WO 3 ). Cyclic voltammetry in aqueous acidic electrolytes shows that the addition of redox-active heteroatoms (Mo 6+ and V 5+ ) leads to systematic shifts in redox couple half-wave potentials (E 1/2 ), broadening, and an overall decrease in the current response. Conversely, the non-redox active heteroatom (Nb 5+ ) only reduces the current response with no observable peak-current broadening. This broadening is attributed to changes in the proton binding affinities of oxygen in different chemical environments, i.e., bridging different pairs of redox-active transition metal cations. We determined the hydrogen bond dissociation free energy (H BDFE) values to elucidate the thermodynamic effect of heteroatom substitution. Density functional theory calculations reveal a differentiation in the hydrogen binding and oxygen vacancy formation energies between heteroatom doped structures. The PICET-induced structural phase transitions of the pristine and doped samples were further probed with operando electrochemical X-ray diffraction (EC-XRD) and with ex situ chemical reduction. The broadening of the potential-dependent current response with increased heteroatom doping manifests in the operando EC-XRD results as prolonged structural regions where multiple hydrogen bronze phases exist and the appearance of cubic bronze phases at lower degrees of reduction compared to pristine WO 3 .

Holzapfel, Noah P. [North Carolina State Universit↗
DOE OSTI2025

Predicting interface structure using the minima hopping method

Here, we adapt the minima hopping method (MHM) to the problem of interfacial structure prediction and apply it to study a canonical problem, the tilt grain boundaries in SrTiO 3 . Our method employs a hybrid approach by first exploring the potential energy surface (PES) of different grain boundary samplings with an empirical force field, among which the fifteen candidates with lower energies are then refined using ab initio density functional theory (DFT) calculations. During the exploratory stage, we bias the search using a local order parameter to primarily sample various reconstructions in the vicinity of the interface, while preserving the crystallinity of the bulk regions. We further enhance the search by incorporating initial structures with rigid body displacements to account for translational variations between bulk phases, enabling the MHM to effectively generate both stoichiometric and nonstoichiometric SrTiO 3 Σ⁢3(111)[110] and Σ⁢3(112)[110] grain boundaries. From an algorithmic standpoint, MHM outperforms earlier studies based on genetic algorithms (GA) by identifying more stable interfacial structures of several SrTiO 3 grain boundaries. The performance of the present implementation of the MHM approach is primarily limited by exploring an approximate description of the PES with a rather simple Buckingham potential. This limitation leads to variations in performance when compared to approaches utilizing more advanced surrogate PES models, such as direct DFT-PES sampling or GA with the embedded atom method (EAM). Despite the present limitations, the MHM approach is able to yield interfacial structures with comparable or lower interfacial energies in specific cases, such as Σ⁢3(111)[110] Γ=1, ±0.5 and Σ⁢3(112)[110] Γ= ±1, −2, underscoring the robustness of the MHM approach even with a simple approximation of the DFT PES. The MHM interfacial structure prediction method thus offers an efficient approach to understanding the grain boundaries and heterointerfaces at the atomic scale, providing an important prerequisite for effective materials design.

density functional theory↗
DOE OSTI2025

Be₂C as a neutron moderator for molten salt reactors: A first-principles study of structural, electronic, and defect properties

Beryllium carbide (Be 2 C), valued for its high neutron moderation efficiency and low absorption cross section, is a promising high-temperature neutron moderator for molten salt reactors. Its practical adoption, however, demands significant technological maturation, requiring comprehensive theoretical and experimental studies of its response to different conditions, including high temperature and irradiation. Here, we report initial results on the fundamental properties and radiation-induced defects of Be 2 C, focusing on antisites, vacancies, interstitial atoms, and Frenkel pairs in the Be and C sublattices. Using density functional theory (DFT) and ab initio molecular dynamics (AIMD), we calculate the defects formation and binding energies, evaluating their dependence on the supercell size, charge states, and chemical environment. In general, carbon defects exhibit higher formation energies, greater sensitivity to cell size, and stronger impacts on the density of states compared to beryllium defects, with charged state the effects being more pronounced. Static DFT reveals multiple metastable interstitial configurations, while AIMD identifies ground states as C-C <100> dumbbells and octahedral Be interstitials. In conclusion, the diversity of metastable configurations and defect states complicates the diffusion mechanisms, requiring further molecular dynamics analysis to elucidate the mechanisms and rates of radiation-induced atomic transport, as well as the structural stability of Be2C.

Ab initio modeling↗
DOE OSTI2025

Real-Space Local Dynamics in 1,2,3-Triazole Using Inelastic Neutron Scattering

Enhancing proton transport in polymer electrolytes is crucial for advancing next-generation solid-state batteries, yet our understanding of proton conductivity in nonaqueous environments remains limited due to a lack of atomic-scale insights. Here, in this study, we investigated the atomic-scale dynamics of 1,2,3-triazole, a small molecule capable of dynamic hydrogen bonding, as a model system for proton hopping in nonaqueous environments. Using the real-space correlation function determined by the double Fourier transformation of inelastic neutron scattering spectra, we identified that the self-motion of protons and intermolecular dynamics occur on comparable time scales. Furthermore, we observed that the activation energy associated with the intermolecular dynamics matches the energy barrier for molecular rotations determined through Density Functional Theory calculations. These findings underscore the importance of controlling molecular dynamics at the atomic scale to control proton transport. Additionally, we demonstrated that intermolecular dynamics in systems involving protons can be studied using inelastic neutron scattering even without deuteration, thereby providing a broader avenue for studying atomic-scale dynamics in soft matter systems.

Shinohara, Yuya [Oak Ridge National Laboratory (OR↗