Homogenization of radiation transport calculations for multi-physics simulations
Explore the source record for details and available documents.
Abrantes, Richard June Espino [Sandia National Lab↗
SEARCH · Search NASA
Results for “density functional calculations”
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.
Cosmic neutrino decoupling and its observable imprints: insights from entropic-dual transport
Abstract Very different processes characterize the decoupling of neutrinos to form the cosmic neutrino background (CνB) and the much later decoupling of photons from thermal equilibrium to form the cosmic microwave background (CMB). The CνB emerges from the fuzzy, energy-dependent neutrinosphere and encodes the physics operating in the early universe in the temperature rangeT∼ 10 MeV toT∼ 10 keV. This is the epoch where beyond Standard Model (BSM) physics, especially in the neutrino sector, may be influential in setting the light element abundances, the necessarily distorted fossil neutrino energy spectra, and other light particle energy density contributions. Here we use techniques honed in extensive CMB studies to analyze the CνB as calculated in detailed neutrino energy transport and nuclear reaction simulations of the protracted weak decoupling and primordial nucleosynthesis epochs. Our moment method, relative entropy, and differential visibility approach can leverage future high precision CMB and light element primordial abundance measurements to provide new insights into the CνB and any BSM physics it encodes. We demonstrate that the evolution of the energy spectrum of the CνB throughout the weak decoupling epoch is accurately captured in the Standard Model by only three parameters per species, a non-trivial conclusion given the deviation from thermal equilibrium and the impact of the decrease of electron-positron pairs. Furthermore, we can interpret each of the three parameters as physical characteristics of a non-equilibrium system. Though the treatment presented here makes some simplifying assumptions including ignoring neutrino flavor oscillations, the success of our compact description within the Standard Model motivates its use also in BSM scenarios. We further demonstrate how observations of primordial light element abundances can be used to place constraints on the CνB energy spectrum, deriving response functions that can be applied for general deviations from a thermal spectrum. Combined with the description of those deviations that we develop here, our methods provide a convenient and powerful framework to constrain the impact of BSM physics on the CνB.
Astronomy & Astrophysics↗
RMCProfile7 : reverse Monte Carlo for multiphase systems
This work introduces a completely rewritten version of the programRMCProfile(version 7), big-box, reverse Monte Carlo modelling software for analysis of total scattering data. The major new feature ofRMCProfile7is the ability to refine multiple phases simultaneously, which is relevant for many current research areas such as energy materials, catalysis and engineering. Other new features include improved support for molecular potentials and rigid-body refinements, as well as multiple different data sets. An empirical resolution correction and calculation of the pair distribution function as a back-Fourier transform are now also available.RMCProfile7is freely available for download at https://rmcprofile.ornl.gov/.
Chemistry↗
Defect-induced displacement of topological surface state in quantum magnet MnBi2Te4
The topological magnet MnBi2Te4 (MBT), with gapped topological surface state, is an attractive platform for realizing quantum anomalous Hall and axion insulator states. However, the experimentally observed surface state gaps fail to meet theoretical predictions, although the exact mechanism behind the gap suppression has been debated. Recent theoretical studies suggest that intrinsic antisite defects push the topological surface state away from the MBT surface, closing its gap and making it less accessible to scanning probe experiments. Here, we report on the local effect of defects on the MBT surface states and demonstrate that high defect concentrations lead to a displacement of the surface states well into the MBT crystal, validating the theorized mechanism. The local and global influence of antisite defects on the topological surface states are studied with samples of varying defect densities by combining scanning tunneling microscopy, angle-resolved photoemission spectroscopy, and density functional theory. Our findings identify a combination of increased defect density and reduced defect spacing as the primary factors underlying the displacement of the surface states and suppression of surface gap, guiding further development of topological quantum materials.
Lupke, Felix [Carnegie Mellon University (CMU)]↗
R&D GREET Battery Carbon Footprint Calculator
The Battery Carbon Footprint (CF) Calculator was developed to help U.S. battery manufacturers meet the carbon footprint reporting requirements of the EU Battery Regulation (EU) 2023/1542. The calculator incorporates several major battery carbon footprint frameworks, including the Joint Research Centre's Rules for the Calculation of the Carbon Footprint of Electric Vehicle Batteries (CFB-EV), RECHARGE's Product Environmental Footprint Category Rules for High Specific Energy Rechargeable Batteries for Mobile Applications (PEFCR), the Catena-X Product Carbon Footprint Rulebook (CX-PCF Rules), Battery Pass's Battery Carbon Footprint: Rules for Calculating the Carbon Footprint of the "Distribution" and "End-of-Life and Recycling" Life Cycle Stages, the Global Battery Alliance's Greenhouse Gas Rulebook: Generic Rules, Version 2.1, and the Ministry of Economy, Trade and Industry's draft Carbon Footprint Calculation Method for Automotive Batteries. The tool pairs these frameworks with foreground data from Argonne's R&D GREET models and integrates user-supplied background data covering battery manufacturing and supply chain activities. By bringing multiple international methodologies together in a single platform, the calculator enables manufacturers to evaluate product carbon footprints, improve data consistency, and prepare for evolving regulatory compliance and global market reporting requirements.
Zhang, Jingyi↗
Experimental and analytical ion thruster research
The results of further spectroscopic studies on the plume from a 3 cm ion source operated on an argon propellant is reported on. In particular, it is shown that it should be possible to use the spectroscopic technique to measure the plasma density of the ion plume close to the grids, where it is difficult to use electrical probe measurements. How the technique, along with electrical probe measurements in the far downstream region of the plume, can be used to characterize the operation of a three-grid, 15 cm diameter thruster from NASA JPL is outlined. Pumping speed measurements on the Vacuum Research Facility have shown that this facility should be adequate for testing the JPL thruster at pressures in the low 10(exp -5) Torr range. Finally, we describe a simple analytical model which can be used to calculate the grid impingement current which results from charge-exchange collisions in the ion plume.
Wilhelmus M. Ruyten↗
Phonon-mediated ultrafast dynamics in self-assembled monolayers of 4-mercaptobenzoic acid on gold
Non-equilibrium interactions between plasmonic metals and adsorbed molecules lie at the heart of emerging applications such as plasmonic photocatalysis and sensing, though the ultrafast charge and energy transfer mechanisms arising from these interactions are not well understood. Herein, we investigate the ultrafast dynamics of Au nano-islands tethered with a self-assembled monolayer (SAM) of electron-withdrawing 4-mercaptobenzoic acid (4MBA) molecules. Ultrafast UV-visible transient absorption spectroscopy following excitation of the interband transition in Au reveals three well-known, characteristic time constants that quantify electron–electron (el–el), electron–phonon (el–ph) and phonon–phonon (ph–ph) scattering lifetimes. When comparing the dynamics of bare Au and 4MBA-Au, we find that the el–ph and ph–ph scattering lifetimes are notably longer in 4MBA-Au. Density functional perturbation theory calculations ascribe the elongation in el–ph lifetimes in 4MBA-Au to the significant coupling of acoustic phonon modes of Au with certain molecular vibrations of 4MBA, leading to decreased spatial overlap between carrier electronic states and the acoustic modes. We speculate that the elongation of ph–ph scattering lifetimes in 4MBA-Au arises due to poor thermal conductivity of the SAM which disrupts efficient energy dissipation from Au to the environment, thus slowing down the thermalization of phonons. This work provides a glimpse into how molecular adsorbates modify the charge carrier and phonon dynamics of Au and sets the stage for further systematic exploration of plasmonic metal–molecule interactions.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Recommended Practices in Laser-Induced Fluorescence (LIF) Diagnostics for Electric Propulsion
Laser-induced fluorescence (LIF) spectroscopy uses Doppler-shifted laser photons to excite a bound electron transition in an atom or ion, with subsequent fluorescence emission detected in order to measure the local velocity distribution function (VDF) and/or particle density. Due to its non-invasive nature, high spatial resolution, and acceptable difficulty of implementation, LIF has become a favored diagnostic for electric propulsion (EP) plasmas, particularly in the study of Hall thrusters. This paper presents recommended best practices for LIF measurements in electric thrusters, drawn from over three decades of implementation heritage in the EP community. Focusing on both single-photon LIF (typically used to measure VDFs) and two-photon LIF (TALIF, typically used to determine densities), the paper covers selection of atomic transitions to target, experimental setups, and interpretation of data. A number of special topics and advanced applications are discussed, including uncertainty analysis, techniques for high-speed measurements, and applications to a variety of propellants.
plasma plume↗
Eddy Covariance Theory: A Review
Eddy covariance (EC), the gold standard for measuring ecosystem scale gas and heat exchanges, has transformed our understanding of the breathing of the biosphere, and thus global change biology. Despite numerous methodological improvements and insights gained from the technique, the community faces persistent challenges that have been present since the first EC measurements. Here, we review the theoretical developments underpinning EC. We present theoretical developments in four important areas that have relevance to EC measurements of the net ecosystem exchanges (NEE) of gases and heat from a single tower: (i) measuring the total vertical flux density, (ii) flux attenuation, (iii) coordinate rotations, and (iv) energy balance closure. Persistent problems with EC measurements, such as the inability to close the energy budget, led us to identify two priorities for revisiting the theory underlying: (i) sensible heat flux calculations, and (ii) constraining the mean vertical wind velocity. We present a framework for improved calculation of sensible heat flux derived from first principles of fluid mechanics and thermodynamics that considers coupled heat and mass transfer so that conservation of both is obeyed. These refinements are motivated by the need for unbiased measurements of energy and mass transfer between the land surface and atmosphere for ecosystem research and to validate satellite observations and land surface models.
ecosystem fluxes↗
Comparative Assessment of Battery Carbon Footprint Calculation Frameworks to Support U.S. Battery Manufacturing
A “battery passport” is a digital record that provides comprehensive information about an individual battery across its life cycle. This concept has been introduced by Battery Regulation (EU) 2023/1542 [1], which applies to batteries for electric vehicles (EVs), industrial batteries with a capacities greater than 2 kWh, and batteries for light means of transport (LMT) greater than 2kWh sold in the European Union (EU) market – impacting manufacturers and exporters across multiple jurisdictions globally. Among its reporting mandates, a key feature of battery passports is the requirement for a carbon footprint (CF) calculation methodology supported by enhanced data granularity to enable traceable and verifiable CF results. While no other jurisdictions have yet adopted formal battery passport requirements like the EU’s, some are developing CF calculation methods to comply with the EU Battery Regulation or are creating CF-related regulations that could evolve in a similar direction, reflecting the growing importance of battery CF guidelines for manufacturers seeking to remain competitive in global markets.
Zhang, Jingyi [Argonne National Laboratory (ANL), ↗
Dimensional Evolution Guides Property Control in the A n Cu 4– n TiS 4 Semiconductor Series
Through progressive reduction of the three-dimensional (3D) covalent network of Cu 4 TiS 4 , we isolate seven new members of the A n Cu 4–n TiS 4 family (A = alkali metal; n = 0–4), spanning 3D, 2D, 1D, and 0D structural fragments. The dimensional reduction is rational, as it preserves the edge-sharing connectivity between [CuS 4 ] 7– and [TiS 4 ] 4– tetrahedra across the series. This structural evolution is driven by the stepwise substitution of Cu with alkali metals, guiding the formation of fragments with reduced dimensionality. The effects of “n” and “A” on the crystal structures, stabilities, electronic structures, and optoelectronic properties are profound, demonstrating that the manipulation of alkali metal size and A n Cu 4–n TiS 4 stoichiometry enables predictable variations in structure and properties. For example, the n = 0 and n = 4 end members of the A n Cu 4–n TiS 4 family set the range of achievable band gaps with 2.00 eV for Cu 4 TiS 4 , 2.60 eV for Na 4 TiS 4 , and intermediate values for the n = 1–3 members. Notably, CsCu 3 TiS 4 exhibits exceptional air stability and congruent melting, with density functional theory (DFT) calculating moderate hole and electron effective masses in specific crystallographic directions (mh = 1.24m 0 , me = 0.87m 0 ). Additionally, A 3 CuTiS 4 (A = Na, K, Rb) displays direct band gap behavior and long photoluminescence lifetimes of 2.3–8.6 μs, and K 3 CuTiS 4 has a PLQY of 5.19%. These findings underscore the potential of the A n Cu 4–n TiS 4 family for applications in optoelectronics and demonstrate widely applicable design concepts that unveil rational stoichiometries within a given composition space to generate a series of crystal structures related through an evolving covalent dimensionality that corresponds to a predictable electronic structure and property progression.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Calculating a minimum overlap period for successful intercalibration of soil moisture sensors
Abstract Long‐term in situ soil moisture monitoring inevitably requires sensors to be replaced. Ensuing discontinuities in the data record can be mitigated by intercalibration, however it is unclear how long the existing sensor needs to remain alongside the newly installed before there is enough overlapping data to generate a robust intercalibration. We used 154 pairs of established and newly installed sensors within the Marena, Oklahoma, In Situ Sensor Testbed to determine if there is a minimum overlap time that should be considered when planning upcoming replacements. Hourly observations of the existing sensor were linearly calibrated to those of the newly installed sensor with coefficients determined from overlap periods incremented by 30 days until a reference period of 2 years was reached. The resulting bias, root‐mean‐square error, and correlation coefficient for sensor pairs indicate that a minimum of 6 to 9 months of overlapping data are required to generate a successful intercalibration. Extending that to a full year before decommissioning the old sensor results in a stable intercalibration with higher confidence.
Agriculture↗
CASM Monte Carlo: Calculations of the thermodynamic and kinetic properties of complex multicomponent crystals
Monte Carlo techniques play a central role in statistical mechanics approaches that connect macroscopic thermodynamic and kinetic properties to the electronic structure of a material. This paper describes the implementation of Monte Carlo techniques for the study of multicomponent crystalline materials within the Clusters Approach to Statistical Mechanics (CASM) software suite, and demonstrates their use in model systems to calculate free energies and kinetic coefficients, study phase transitions, and construct phase diagrams from first principles. Many crystal structures are complex, with multiple sublattices occupied by differing sets of chemical species, along with the presence of vacancies or interstitial species. This imposes constraints on concentration variables, the form of thermodynamic potentials, and the values of kinetic transport coefficients. The framework used by CASM to formulate thermodynamic potentials and kinetic transport coefficients accounting for arbitrarily complex crystal structures is presented and demonstrated with examples of increasing complexity. Additionally, an overview of the capabilities of the CASM software specific to Monte Carlo methods is given, and a new CASM software package is introduced, casm-flow, which helps automate the setup, submission, management, and analysis of Monte Carlo simulations.
Cluster expansion↗
The CP-PAW Code Package for First-Principles Calculations from a User’s Perspective
CP-PAW is a combined electronic structure and ab initio molecular dynamics code to perform mixed quantum and classical simulations of atomistic condensed phase systems, such as solids, liquids, and molecular systems. As the name suggests, the CP-PAW code unifies the all-electron projector augmented-wave (PAW) method with the Car–Parrinello (CP) approach to determine not only the electronic and nuclear ground states of condensed matter but also to study their properties and dynamics. In addition to briefly outlining the underlying theory, the focus will be on the unique aspects of CP-PAW and how to correctly employ them as a user. How to install CP-PAW using the new build system will also be briefly mentioned.
Blöchl, Peter E [Institute for Theoretical Physic↗
Reducing experimental uncertainty in the calculation of cavity swelling in a transmission electron microscope through crystallographically aided void volume tracking
Accurate measurement of cavity swelling in a transmission electron microscope is essential to define material performance under irradiation, and the conventionally applied spherical assumption for the calculation of cavity volumes in irradiated materials can result in errors ranging between a 25% underestimation and 72% overestimation of volume purely based on the assumption of shape. This assumption is undeniably expedient for calculation but does not sufficiently account for the 3D nature of polyhedral cavities and their shape projection in the transmission electron microscope, and therefore presents too large of an associated uncertainty in swelling determination for faceted cavities. This uncertainty has been defined for many common cavity shapes in FCC and BCC crystal systems, and has been tabulated across the cubic fundamental region. A revised methodology for crystallographically aided void volume tracking, or CAVV-T, is presented and demonstrated on a specimen of neutron irradiated Ni. In-depth discussion on the application of this technique is provided along with resources to allow for conversion between the spherical assumption and this revised method. This work seeks to increase experimental confidence in the characterization and quantification of critical aspects of irradiation damage in materials by applying a crystallographically-resolved approach for cavity swelling calculation.
Haag IV, James V. [Pacific Northwest National Labo↗
Charmonium Transport in Heavy-Ion Collisions at the LHC
We provide an update on our semi-classical transport approach for quarkonium production in high-energy heavy-ion collisions, focusing on 𝐽/𝜓 and 𝜓(2𝑆) mesons in 5.02 TeV Pb-Pb collisions at the Large Hadron Collider (LHC) at both forward and mid-rapidity. In particular, we employ the most recent charm-production cross sections reported in pp collisions, which are pivotal for the magnitude of the regeneration contribution, and their modifications due to cold-nuclear-matter (CNM) effects. Multi-differential observables are calculated in terms of nuclear modification factors as a function of centrality, transverse momentum, and rapidity, including the contributions from feeddown from bottom hadron decays. For our predictions for 𝜓(2𝑆) production, the mechanism of sequential regeneration relative to the more strongly bound 𝐽/𝜓 meson plays an important role in interpreting recent ALICE data.
Astronomy & Astrophysics↗
Calculations of the Moon's Thermal History at Different Concentrations of Radioactive Elements, Taking Into Account Differentiation on Melting
Calculations of the thermal history of the moon were carried out by solving the thermal conductivity equation for the case in which the heat sources are the long-lived radioactive elements Th, U, and K 40 . The concentrations of these elements were adjusted to give four variations of the heat flow: 1.35 x 10 -8 and 0.91 x 10 -6 cal cm -2 s -1 (I and I**), 0.61 x 10 -6 (variant II, the terrestrial mixture of Lyubimova and Starkova), and 0.236 x 10 -6 (variant III, the chondrite model of Urey and MacDonald). In the same calculations, we considered layering of the differentiated material with transport to the surface of the radioactive elements after the temperature of the layer rose to 200 K above the melting temperature, which is considered in five variants that differ in the amount of transported radioactive elements: 100 percent ( n = 1), 80 percent n = 0.8), 40 percent, 20 percent ( n = 0.6, 0.4, 0.2). During fusion the heat capacity and heat conductivity were changed. We considered two variants of an initially cold (273 K) and hot (900 K) Moon. Calculations show that the interior of the Moon was heated to melting during the first 0.7 to 2.3 x 10 9 years. The maximum fusion involved practically the entire Moon to a distance from 15 to 45 km beneath the surface, and started 3.5 to 4.0 x 10 9 years ago (I, I**), or 2.5 to 3.0 x 10 9 years ago (II, III) and continued for 1 to 2 x 10 9 years. Today the Moon is cooling. The current thickness of the solid crust is from 150 to 200 km and the heat flow exceeds the stationary value 1.5 fold. Apparently the most realistic variant is II (terrestrial mixture) for an initially hot Moon, and gives, regardless of the moderate concentration of radioactive elements, a heat flow of 0.9 to 0.95 x 10 -6 cal cm -2 s -1 , which agrees with the radioastronomical measurements of Troitsky and Krotikov and with the Apollo 15 data.
O I Ornatskaya↗
Thermodynamically informed priors for uncertainty propagation in first-principles statistical mechanics
Here, this work demonstrates how first-principles statistical mechanics approaches within a Bayesian framework can quantify and propagate uncertainties to downstream thermodynamic calculations. To address the issue of Bayesian prior selection, knowledge of 0 K ground states in the material system of interest is incorporated into the prior. The effectiveness of this framework is shown by creating a phase diagram for the fcc zirconium nitride system, including confidence intervals on order-disorder transition temperatures.
Bayesian methods↗