Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular data”

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

Computing inelastic neutron scattering spectra from molecular dynamics trajectories

Abstract Inelastic neutron scattering (INS) provides a weighted density of phonon modes. Currently, INS spectra can only be interpreted for perfectly crystalline materials because of high computational cost for electronic simulations. INS has the potential to provide detailed morphological information if sufficiently large volumes and appropriate structural variety are simulated. Here, we propose a method that allows direct comparison between INS data with molecular dynamics simulations, a simulation method that is frequently used to simulate semicrystalline/amorphous materials. We illustrate the technique by analyzing spectra of a well-studied conjugated polymer, poly(3-hexylthiophene-2,5-diyl) (P3HT) and conclude that our technique provides improved volume and structural variety, but that the classical force field requires improvement before the morphology can be accurately interpreted.

36 MATERIALS SCIENCE↗
NASA NTRS1994

Far Outer Galaxy H II Regions

We have made a multifrequency (6, 3.6, and 2 cm), high-resolution (3"-6"), radio continuum survey of IRAS selected sources from the catalogue of Wouterloot & Brand (1989) to search for and study H II regions in the far outer Galaxy. We identified 31 sources in this catalog with well determined galactocentric distances, and with R approx.. greater than 15 kpc and L(sub FIR) approx.greater than 10(exp 4) solar luminosity, indicating the presence of high-mass star-formation. We have observed 11 of these sources with the Very Large Array (VLA). We observed the sources at 6 and 2 cm using "scaled arrays", making possible a direct and reliable comparison of the data at these two wavelengths for the determination of spectral indices. We detected a total of 12 radio sources, of which 10 have spectral indices consistent with optically-thin free-free emission from H II regions. Combined with previous VLA observations by other investigators, we have data on a total of 15 H II regions at galactocentric distances of 15 to 18.2kpc, among the most remote H II regions found in our Galaxy. The sizes of the H II regions range from approx. less than 0.10 to 2.3 pc. Using the measured fluxes and sizes, we determine the electron densities, emission measures, and excitation parameters of the H II regions, as well as the fluxes of Lyman continuum photons needed to keep the nebulae ionized. The sizes and electron densities are consistent with most of the sources detected in this survey being compact or ultracompact H II regions. Seven of the fifteen H II regions have sizes approx. less than 0.20 pc. Assuming simple pressure-driven expansion of the H II regions, these sizes indicate ages approx. less than 5 x 10(exp 4) yr, or only 1% of the lifetime of an O star, which implies an unlikely overabundance of O stars in the outer Galaxy. Thus, the large number of compact H II regions suggests that the time these regions spend in a compact phase must be much longer than their dynamical expansion times. Five of the fifteen H II regions have cometary shapes; the remainder are spherical or unresolved. Comparison of the radio continuum data with molecular line maps suggests that the cometary shape of the two H II regions in S 127 may be due to pressure confinement of the expanding ionized gas, as in the "blister" or "champagne flow" models of H II regions. Comparison of the radio continuum data with the IRAS far-infrared data indicates that the five most luminous H II regions are consistent with a single 0 or B star exciting a dust-free H II region. Subject headings: stars: formation - ISM: H II regions - ISM: individual objects: S 127 radio continuum: interstellar

Rudolph, A. L.↗
DOE OSTI2023

Rate coefficients for rotational state-to-state transitions in H 2 O + H 2 O collisions for cometary and planetary applications, as predicted by mixed quantum-classical theory

We present new calculations of collision cross sections for state-to-state transitions between the rotational states in an H 2 O + H 2 O system, which are used to generate a new database of collisional rate coefficients for cometary and planetary applications. Calculations were carried out using a mixed quantum-classical theory approach that is implemented in the code MQCT. The large basis set of rotational states used in these calculations permits us to predict thermally averaged cross sections for 441 transitions in para- and ortho-H 2 O in a broad range of temperatures. It is found that all state-to-state transitions in the H 2 O + H 2 O system split into two well-defined groups, one with higher cross-section values and lower energy transfer, which corresponds to the dipole-dipole driven processes. The other group has smaller cross sections and higher energy transfer, driven by higher-order interaction terms. We present a detailed analysis of the theoretical error bars, and we symmetrized the state-to-state transition matrixes to ensure that excitation and quenching processes for each transition satisfy the principle of microscopic reversibility. We also compare our results with other data available from the literature for H 2 O + H 2 O collisions.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2001

The molecular responses of skeletal muscle satellite cells to continuous expression of IGF-1: implications for the rescue of induced muscular atrophy in aged rats

Approximately 50% of humans older than 85 years have physical frailty due to weak skeletal muscles. This indicates a need for determining mechanisms to combat this problem. A critical cellular factor for postnatal muscle growth is a population of myogenic precursor cells called satellite cells. Given the complex process of sarcopenia, it has been postulated that, at some point in this process, a limited satellite cell proliferation potential could become rate-limiting to the regrowth of old muscles. It is conceivable that if satellite cell proliferative capacity can be maintained or enhanced with advanced age, sarcopenia could potentially be delayed or prevented. Therefore, the purposes of this paper are to describe whether IGF-I can prevent muscular atrophy induced by repeated cycles of hindlimb immobilization, increase the in vitro proliferation in satellite cells from these muscles and, if so, the molecular mechanisms by which IGF-I mediates this increased proliferation. Our results provide evidence that IGF-I can enhance aged muscle regrowth possibly through increased satellite cell proliferation. The results also suggest that IGF-I enhances satellite cell proliferation by decreasing the cell cycle inhibitor, p27Kip1, through the PI3'-K/Akt pathway. These data provide molecular evidence for IGF-I's rescue effect upon aging-associated skeletal muscle atrophy.

NASA Discipline Musculoskeletal↗
NASA NTRS2022

Physical and Chemical Conditions in the N113 Star-Forming Region in the Low-Metallicity Large Magellanic Cloud

The Large Magellanic Cloud (LMC) is the nearest (~50 kpc) star-forming galaxy characterized by a low metallicity (Z~0.3-0.5 Z⨀) similar to galaxies during the early phases of their assembly. As a result, star formation studies in the LMC provide a stepping stone to understanding star formation at earlier epochs of the universe where these processes cannot be directly observed. N113 is one of the most prominent star-forming regions in the LMC hosting one of the most massive giant molecular clouds.N113 is small enough to be imaged in its entirety, but large enough to showcase many important phenomena such as multiple generations of stars, stellar feedback, and different environments. We present our findings from an investigation of the early stages of star formation in the N113 region using the Atacama Large Millimeter/submillimeter Array (ALMA)molecular line data probing a wide density range: 12CO, 13CO, and C18O (2-1), 13COand C18O (1-0), HCN (1-0), HCO+ (1-0), H13CN (1-0) and (3-2), H13CO+ (1-0) and (3-2), CS (2-1) and (5-4), as well as 1.3 mm and 3 mm continuum. We used the Python package quick clump to identify molecular clumps. We utilized the multiline non-LTE fitting tool based on models from RADEX developed by Finn et al. (2021, ApJ, 917, 106)to construct the CO, HCN, HCO+, and CS temperature and column density, and the H2density maps of N113. We constructed a catalog of molecular clumps including their physical properties, chemical abundances, sizes, velocities, and velocity dispersions. To establish the evolutionary status of the clumps, their positions were compared with previously identified young stellar objects (YSOs) from the Spitzer/SAGE and Herschel/HERITAGE surveys, as well as water and OH masers. We compared the properties of the clumps in N113 to those in the Galaxy and other regions in the LMC to assess the impact of the environment (e.g., metallicity, stellar feedback) on the star formation process.

Jonathon Noswitz↗
NASA NTRS2023

The First Extragalactic Detection of Higher-Order Hydrogen Recombination Lines

The Large Magellanic Cloud (LMC) is the nearest (~50 kpc) star-forming galaxy characterized by a low metallicity (Z~0.3-0.5 Z⨀) similar to galaxies during the early phases of their assembly. As a result, star formation studies in the LMC provide a stepping stone to understanding star formation at earlier epochs of the universe where these processes cannot be directly observed. N113 is one of the most prominent star-forming regions in the LMC hosting one of the most massive giant molecular clouds. N113 is small enough to be imaged in its entirety, but large enough to showcase many important phenomena such as multiple generations of stars, stellar feedback, and different environments. We present our findings from an investigation of the early stages of star formation in the N113 region using the Atacama Large Millimeter/submillimeter Array (ALMA) molecular line data probing a wide density range: 12CO, 13CO, and C18O (2-1), 13CO and C18O (1-0), HCN (1-0), HCO+ (1-0), H13CN (1-0) and (3-2), H13CO+ (1-0) and (3-2), CS (2-1) and (5-4), as well as 1.3 mm and 3 mm continuum. We used the Python package quickclump to identify molecular clumps. We utilized the multiline non-LTE fitting tool based on models from RADEX developed by Finn et al. (2021, ApJ, 917, 106) to construct the CO, HCN, HCO+, and CS temperature and column density, and the H2 density maps of N113. We constructed a catalog of molecular clumps including their physical properties, chemical abundances, sizes, velocities, and velocity dispersions. To establish the evolutionary status of the clumps, their positions were compared with previously identified young stellar objects (YSOs) from the Spitzer/SAGE and Herschel/HERITAGE surveys, as well as water and OH masers. We compared the properties of the clumps in N113 to those in the Galaxy and other regions in the LMC to assess the impact of the environment (e.g., metallicity, stellar feedback) on the star formation process.

Marta M Sewilo↗
NASA NTRSDate not supplied

Multiscale Modeling of Thermoplastics Using Atomistic-informed Micromechanics

A multiscale model was developed for predicting the thermoelastic behavior of semi-crystalline thermoplastic materials for composite aerospace applications. At the highest scale containing the semi-crystalline spherulite in an amorphous matrix, the generalized method of cells, or high fidelity method of cells, was used to perform the homogenization calculations to obtain the effective properties. Models were developed assuming a cubic, or spherical shape, for the spherulite to understand if the morphology of the spherulite affects the effective thermoelastic properties. The generalized method of cells was used to model at the repeating unit cells at the subscales of the microstructure including the lamellae stacks and granular crystal blocks. The scales are integrated using the multiscale micromechanics method in the NASA Multiscale Analysis Tool. Data from molecular dynamics simulations were used as inputs for the amorphous and crystalline constituents. Convergence studies were performed to determine the best level of discretization for the repeating unit cell at the highest scale. Effective Young’s modulus, shear modulus, Poisson’s ratio, coefficient of thermal expansion, and thermal conductivity were predicted for polyether ether ketone and polyether ketone ketone, and very good agreement between the model utilizing the cubic spherulite and the experimental data, where available, was observed for polyether ketone ketone. Normalization of the data for the bulk polyether ketone ketone, against amorphous data, improved the predictions as compared to experimental data. Overall, the high fidelity method of cells predicted a stiffer response then the generalized method of cells as the crystallinity was increased. The shape of the spherulite had a minimal effect on the predicted bulk properties of the polymers.

thermoplastics↗
NASA NTRS1992

Interstellar cloud shapes - A minimum-hypothesis account involving a purely gravitational effect

The effect produced by the anisotropic gravitational field on the shape of a nonspherical, self-gravitating interstellar cloud is examined. The projected axial ratios observed for molecular clouds can be accounted for if clouds are not generally in a state of strict dynamical equilibrium. For clouds having spectral line widths dominated by supersonic turbulent motions, this purely gravitational effect predicts a mean true axial ratio q about 0.3 for both oblate and prolate cloud geometries. Clouds having p less than about 0.3 are likely to be prolate and quite elongated. The observational data on molecular cloud cores, which have essentially thermal line widths, can be fitted equally well to either oblate or prolate spheroids having q between 0.1 and 0.4 with q about 0.2 providing the best match: highly elongated cores are again more likely to be prolate. While other processes certainly operate throughout the Galaxy to influence cloud shape, it appears plausible that at least some clouds may be shaped principally by their anisotropic gravitational field.

Fleck, Robert C., Jr.↗
DOE OSTI2024

Applying the Velocity Gradient Technique in NGC 1333: Comparison with Dust Polarization Observations

Magnetic fields (B-fields) are ubiquitous in the interstellar medium (ISM), and they play an essential role in the formation of molecular clouds and subsequent star formation. However, B-fields in interstellar environments remain challenging to measure, and their properties typically need to be inferred from dust polarization observations over multiple physical scales. In this work, we seek to use a recently proposed approach called the velocity gradient technique (VGT) to study B-fields in star-forming regions and compare the results with dust polarization observations in different wavelengths. The VGT is based on the anisotropic properties of eddies in magnetized turbulence to derive B-field properties in the ISM. We investigate that this technique is synergistic with dust polarimetry when applied to a turbulent diffused medium for the purpose of measuring its magnetization. Specifically, we use the VGT on molecular line data toward the NGC 1333 star-forming region ( 12 CO, 13 CO, C 18 O, and N 2 H + ), and we compare the derived B-field properties with those inferred from 214 and 850 μm dust polarization observations of the region using Stratospheric Observatory for Infrared Astronomy/High-Resolution Airborne Wide-band Camera Plus and James Clerk Maxwell Telescope/POL-2, respectively. We estimate both the inclination angle and the 3D Alfvénic Mach number M A from the molecular line gradients. Crucially, testing this technique on gravitationally bound, dynamic, and turbulent regions, and comparing the results with those obtained from polarization observations at different wavelengths, such as the plane-of-sky field orientation, is an important test on the applicability of the VGT in various density regimes of the ISM. We in general do not find a close correlation between the velocity gradient inferred orientations and the dust inferred magnetic field orientations.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2013

Rovibrational Spectroscopic Constants and Fundamental Vibrational Frequencies for Isotopologues of Cyclic and Bent Singlet HC2N isomers

Through established, highly-accurate ab initio quartic force fields (QFFs), a complete set of fundamental vibrational frequencies, rotational constants, and rovibrational coupling and centrifugal distortion constants have been determined for both the cyclic 1(sup 1) 1A' and bent 2(sup 1)A' DCCN, H(C13)CCN, HC(C-13)N, and HCC(N-15) isotopologues of HCCN. Spectroscopic constants are computed for all isotopologues using second-order vibrational perturbation theory (VPT2), and the fundamental vibrational frequencies are computed with VPT2 and vibrational configuration interaction (VCI) theory. Agreement between VPT2 and VCI results is quite good with the fundamental vibrational frequencies of the bent isomer isotopologues in accord to within a 0.1 to 3.2 / cm range. Similar accuracies are present for the cyclic isomer isotopologues. The data generated here serve as a reference for astronomical observations of these closed-shell, highly-dipolar molecules using new, high-resolution telescopes and as reference for laboratory studies where isotopic labeling may lead to elucidation of the formation mechanism for the known interstellar molecule: X 3A0 HCCN.

molecular data↗
DOE OSTI2022

ORNL_AISD_NiNb

This dataset describes the nickel-niobium solid solution binary alloy, where the two constituent elements nickel (Ni) and niobium (Nb) are randomly placed on an underlying crystal lattice. This dataset for nickel-niobium (Ni-Nb) alloys available includes the formation energy and bulk modulus for each crystal structure. Each atomic sample has a disordered phase which is obtained starting from an initial regular crystal structure of type body-centered cubic (BCC), face-centered cubic (FCC), or hexagonal compact packed (HCP). The geometry optimization ensures that all the alloy samples reached the equilibrium with negative formation energy. We perform geometry optimizations using the LAMMPS simulation package [1], a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. We utilized the embedded atom model (EAM) potential for Ni and Nb developed in a previous study [2]. The potential could describe behaviors of the liquid and solid phases of Ni-Nb alloy. The structural factors and angular distributions of three atoms are well-matched with X-ray and ab initio-based molecular dynamics data. We prepared the three different crystals with different initial lattice parameters (3.52 Ã… for FCC, 3.32 Ã… for BCC, and 3.5 Ã… for HCP). We performed energy minimization in two steps. Firstly, we minimized the structures with an isotropic unit cell to minimize the side effects from our arbitrary lattice parameters for all other compositions. Then, we applied geometry optimization with a triclinic (non-orthogonal) unit cell to fully minimize the stress components to calculate the elastic constants. In this procedure, we chose 10,000 as the maximum number of allowable steps aimed at obtaining fully relaxed atomic geometries. The dataset consists of three sets of crystal structures. The first set contains 46,086 irregular crystal structures, each of them with 54 atoms, obtained through optimization starting from a regular BCC crystal structure. The second set contains 24,543 irregular crystal structures, each of them with 32 atoms, obtained through optimization starting from a regular FCC crystal structure. The third set contains 39,303 irregular crystal structures, each of them with 48 atoms, obtained through optimization starting from a regular HCP crystal structure. The atomic configurations within each set span the possible compositional range. The three sets have been unified in a global dataset, which is extremely heterogeneous in terms of crystal structures, lattice volumes, and atomic configurations. Organization of files inside the dataset: the dataset contains three subdirectories called • BCC_opt • FCC_opt • HCP_opt based on the type of initial regular structure used to start the geometry optimization. Inside each of these folders, every atomic structure is identified by a string “A_B_Câ€, where A denotes the number of Nb in the system, B denotes index of structure with a given Nb number, and C denotes the total number of structures generated with a given Nb number. For each optimized crystal structure identified by the unique string of characters “A_B_Câ€, three files are provided: • A_B_C_opt.xyz: The optimized geometries in xyz format • A_B_C_opt.cfg: The optimized geometries in cfg format. It includes cell information and atomic energy, and forces calculated from LAMMPS. • A_B_C.elastic: Raw data of 21 elastic constants from LAMMPS output. • A_B_C.bulk: Calculated upper and lower bounds of bulk modulus and averaged one based on Voigt-Reuss-Hill approach from *.elastic. References: [1] A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolintineanu, W. M. Brown, P. S. Crozier, P. J. in 't Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, R. Shan, M. J. Stevens, J. Tranchida, C. Trott, and S. J. Plimpton. LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. Comp. Phys. Comm., 271:108171, 2022. [2] Y Zhang, R Ashcraft, MI Mendelev, CZ Wang, and KF Kelton. Experimental and molecular dynamics simulation study of structure of liquid and amorphous ni62nb38 alloy. The Journal of chemical physics, 145(20):204505, 2016.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Elasticity of two-dimensional ferroelectrics across their paraelectric phase transformation

The mechanical behavior of two-dimensional (2D) materials across 2D phase changes is unknown, and the finite temperature (T) elasticity of paradigmatic SnSe monolayers—ferroelectric 2D materials turning paraelectric as their unit cell turns from a rectangle into a square—is described here in a progressive manner. To begin with, their zero–T total energy landscape gives way to (Boltzmann-like) averages from which the elastic behavior is determined. Furthermore, these estimates are complemented with results from the strain-fluctuation method, which employs the energy landscape or ab initio molecular dynamics data. All approaches capture the coalescence of elastic moduli < C 11 (T) > = < C 22 (T) > due to the structural transformation. The broad evolution and sudden changes of elastic parameters < C 11 (T) >, < C 22 (T) >, and < C 12 (T) > of these atomically thin phase-change membranes establishes a heretofore overlooked connection among 2D materials and soft matter.

2-dimensional systems↗
DOE OSTI2022

Atomistic modeling of LiF microstructure ionic conductivity and its influence on nucleation and plating

We report the formation and degradation of the solid electrolyte interphase (SEI) and its underlying transport properties play an essential role in the overall performance of lithium-ion batteries. This paper presents classical molecular dynamics studies on polycrystalline inorganic lithium fluoride (LiF) layers to model and predict the SEI transport properties. The ionic conductivity is obtained from the lithium-ion diffusivity in polycrystalline structures of LiF using the Nernst-Einstein relation. The predicted molecular dynamics data are used in a continuum scale phase-field model to evaluate the plating kinetics under fast charging conditions. The analysis emphasizes that the SEI ionic conductivity properties impact the plating dynamics, where SEI's low ion conductivity value is prone to large plating and subsequent capacity degradation. The combination of atomic and continuum scale studies shown herein lays a foundation to tune in SEI transport properties to decrease the amount of lithium plating and improve the performance of fast-charging batteries.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Soil microbiome resilience to short-term (30 days, 90 days) and long-term (1000 days) drought

This dataset contains data used for the paper "Drought duration does not impact soil microbiome resilience". The Related References will be updated with a full citation when available. Increasing global droughts exert large but poorly understood effects on the microbial communities and ecology of soil. Microbial communities generally show resilience and return to pre-drought conditions when short-term droughted soils are rewet; soils exposed to long-term drought, however, often show a lag upon rewetting, after which microbial communities may or may not return to their pre-stressed conditions. Though short-term droughts have been widely studied, long-term drought manipulation experiments remain rare, especially those that compare microbial response to short-term and long-term drought in tandem. We conducted a 1000-day drought simulation in controlled laboratory conditions with soil cores collected from a tidal freshwater ecosystem in Washington state, USA, and subsequently exposed them to rewetting for two weeks. We also included short-term (30-day and 90-day) drought and rewet treatments to directly compare microbial community and organic matter responses across drought durations. We found distinct microbial taxa belonging to Firmicutes and Actinobacteria enriched after the 1000-day drought, but not after the short-term droughts. While we hypothesized that the microbial community would recover from a short-term drought after rewetting to resemble pre-drought conditions, our results revealed community dissimilarities between rewet and pre-drought conditions across all drought durations. These findings suggest unique microbial life history strategies within certain microbial phyla that make them successful colonizers during an extended drought period, and the influence of environmental and physiological context on microbial responses to rewetting. The 16SrRNA gene amplicon dataset contains processed DNA sequences in the form of an ASV table with raw unrarefied read counts and representative sequences in .fasta format as described in the ESS-DIVE amplicon sequence reporting format (https://ess-dive.gitbook.io/amplicon-sequencing-reporting-format/instructions). The Fourier Transform Ion Cyclotron Resonance Mass Spectrometry (FTICR-MS) dataset consists of processed files containing presence absence data of molecular formulae and molecular characterization of FTICR resolved peaks. The Nuclear Magnetic Resonance (NMR) dataset contains files relevant to NMR spectra and peaks. A sample key file and a sample metadata file is included for the FTICR/NMR and 16S dataset respectively.

1000-day drought↗
DOE OSTI2023

RNA target highlights in CASP15 : Evaluation of predicted models by structure providers

Abstract The first RNA category of the Critical Assessment of Techniques for Structure Prediction competition was only made possible because of the scientists who provided experimental structures to challenge the predictors. In this article, these scientists offer a unique and valuable analysis of both the successes and areas for improvement in the predicted models. All 10 RNA‐only targets yielded predictions topologically similar to experimentally determined structures. For one target, experimentalists were able to phase their x‐ray diffraction data by molecular replacement, showing a potential application of structure predictions for RNA structural biologists. Recommended areas for improvement include: enhancing the accuracy in local interaction predictions and increased consideration of the experimental conditions such as multimerization, structure determination method, and time along folding pathways. The prediction of RNA–protein complexes remains the most significant challenge. Finally, given the intrinsic flexibility of many RNAs, we propose the consideration of ensemble models.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Integration of a Wigner effect-based energy storage system with an advanced nuclear reactor

In this work, an innovative energy storage concept based on the purposeful creation of defects in crystalline material by neutron irradiation is presented. Lattice defects are generated when heavy particles collide with the atoms in a crystal structure, i.e., if the incoming particles have enough energy, recoil atoms are displaced from their initial lattice sites. Most of the displaced atoms will eventually combine with nearby vacancies, but some of them will come to rest in non-ideal locations. The energy held by displaced atoms is called Wigner energy. Lattice defects can migrate and form clusters, and the Wigner energy can be released from these groupings if sufficient activation energy is provided. In the nuclear industry, this effect is well-known since it represented an issue for graphite-moderated reactors. This work presents the conceptual design of an engineering system that exploits this physical process to store the energy of neutrons in advanced reactor concepts. In the first part of the paper, the theoretical performance of an energy storage system based on the Wigner effect is described. Given the lattice properties and the compatibility with the harsh reactor environment, graphite was selected as the candidate material for the irradiation targets. Both experimental data and molecular dynamics simulations confirmed that this system can achieve performance comparable with state-of-the-art batteries in terms of stored energy density. In the second part of the paper, the engineering challenges of this innovative technology and the proposed solutions are described. After defining the optimal irradiation conditions, the different steps of the operation of the proposed energy system (from energy storing to energy harvesting) were defined. Finally, the integration of this concept with advanced reactor designs, i.e., a Sodium-cooled Fast Reactor and a Molten Salt-cooled Reactor, was investigated and the corresponding performance was evaluated.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

The divergence of nearby trajectories in soft-sphere DEM

The n-body instability is investigated with the soft-sphere discrete element method. The divergence of nearby trajectories is quantified by the dynamical memory time. Using the inverse proportionality between the dynamical memory time and the largest Lyapunov exponent, the soft-sphere discrete element method results are compared to previous hard-sphere molecular dynamics data for the first time. Good agreement is observed at low concentrations and the degree of instability is shown to increase asymptotically with increasing spring stiffness. At particle concentrations above 30%, the soft-sphere Lyapunov exponents increase faster than the corresponding hard-sphere data. Finally, this paper concludes with a demonstration of how this case study may be used in conjunction with regression testing and code verification activities.

42 ENGINEERING↗
DOE OSTI2025

Role of Electron Correlation beyond the Active Space in Achieving Quantitative Predictions of Spin-Phonon Relaxation

Single-molecule magnets (SMMs) are promising candidates for molecular-scale data storage and processing due to their strong magnetic anisotropy and long spin relaxation times. However, as temperature rises, interactions between electronic states and lattice vibrations accelerate spin relaxation, significantly limiting their practical applications. Recently, ab initio simulations have made it possible to advance our understanding of phonon-induced magnetic relaxation, but significant deviations from experiments have often been observed. The description of molecules’ electronic structure has been mostly based on complete active space self-consistent field (CASSCF) calculations, and the impact of electron correlation beyond the active space remains largely unexplored. In this study, we provide the first systematic investigation of spin-phonon relaxation in SMMs with post-CASSCF multiconfigurational methods, specifically CAS followed by second-order perturbation theory and multiconfiguration pair-density functional theory. Taking Co(II)- and Dy(III)-based SMMs as case studies, we analyze how electron correlation influences spin-phonon relaxation rates across a range of temperatures, comparing theoretical predictions with experimental observations. Our findings demonstrate that post-CASSCF treatments make it possible to achieve quantitative predictions for Co(II)-based SMMs. For Dy(III)-based systems, however, accurate predictions require consideration of additional effects, underscoring the urgent necessity of further advancing the study of the effects of electronic correlation in these complex systems.

Energy↗