Search NASA⌕ Search

SEARCH · Search NASA

Results for “Non-equilibrium thermodynamics”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Elucidating key reducing species beyond ions in hydrogen plasma smelting reduction of iron ore

Hydrogen plasma smelting reduction (HPSR) of iron ore has attracted significant attention over the past decade due to its high-temperature operation, rapid plasma mediated reduction kinetics, and simpler density-based separation of molten iron product, compared to H2-based solid-state reduction. All of these attributes enable processing of low-grade ores for downstream use in electric-arc furnaces, as virgin iron with low gangue content is required for high quality steel and improved furnace operation. While positive ions exist within the plasma arc, this work demonstrates that near the anodic ore surface, hydrogen radicals and vibrationally excited hydrogen species dominate and their densities correlate well with observed reduction rates. Species concentrations in the transferred plasma arc and at the plasma-ore interface are evaluated using coupled thermal plasma and near-wall non-equilibrium plasma models. The thermal plasma model is validated against experimental voltage data and spectroscopic measurements of plasma temperature and density for varying current inputs. Modeling of the near surface thermochemical non-equilibrium and micrometer scale anode sheath layer reveals, in addition to the expected H + , significant concentrations of ArH + and H$^+_3$ ions, typically not observed in thermal plasmas under thermodynamic equilibrium. Our results show that the inverted sheath structure at the anodic ore surface strongly suppresses reactive positive ion fluxes, while non-equilibrium electron-impact processes generate abundant hydrogen radicals and vibrationally excited species. These findings highlight the critical role of non-equilibrium effects in hydrogen arc-driven iron ore reduction and advance understanding beyond prevailing hypotheses centered on hydrogen ion-driven mechanisms.

08 HYDROGEN↗
NASA NTRS2002

Laser Induced Aluminum Surface Breakdown Model

Laser powered propulsion systems involve complex fluid dynamics, thermodynamics and radiative transfer processes. Based on an unstructured grid, pressure-based computational aerothermodynamics; platform, several sub-models describing such underlying physics as laser ray tracing and focusing, thermal non-equilibrium, plasma radiation and air spark ignition have been developed. This proposed work shall extend the numerical platform and existing sub-models to include the aluminum wall surface Inverse Bremsstrahlung (IB) effect from which surface ablation and free-electron generation can be initiated without relying on the air spark ignition sub-model. The following tasks will be performed to accomplish the research objectives.

Chen, Yen-Sen↗
NASA NTRS2001

Laser Induced Aluminum Surface Breakdown Model

Laser powered propulsion systems involve complex fluid dynamics, thermodynamics and radiative transfer processes. Based on an unstructured grid., pressure-based computational aerothermodynamics, platform, several sub-nio"'dels describing such underlying physics as laser ray tracing and focusing, thermal non-equilibrium, plasma radiation and air spark ignition have been developed. This proposed work shall extend the numerical platform and existing sub-models to include the aluminum wall surface Inverse Bremsstrahlung (113) effect from which surface ablation and free-electron generation can be initiated without relying on the air spark ignition sub-model. The following tasks will be performed to accomplish the research objectives.

Chen, Yen-Sen↗
NASA NTRS1992

Innovative techniques for the production of energetic radicals for lunar processing including cold plasma processing of local planetary ores

Hydrogen reduction of ilmenite has been studied by a number of investigators as a potential means for recovery of oxygen from lunar soil. Interest in this process has always rested with the simplicity of the flow diagram and the utilization of established technology. Effective utilization of hydrogen in the reduction process at temperatures of 1200 C and below has always been disappointing and, as such, has led other investigators to focus attention on other systems. Effective utilization of hydrogen in the reduction of ilmenite can be significantly enhanced in the presence of a non-equilibrium hydrogen plasma. Ilmenite at solid specimen temperatures of 600 C to 970 C were reacted in a hydrogen plasma. Those experiments revealed that hydrogen utilization can be significantly enhanced. At a specimen temperature of 850 C the fraction of H2 reacted was 24 percent compared to the 7 percent theoretical limit calculated with thermodynamic theory for the same temperature. An added advantage for a hydrogen plasma involves further reduction of TiO2. Reduction of the iron oxide in ilmenite yields TiO2 and metallic iron as by products. Titanium forms a number of oxides including TiO, Ti2O3, Ti3O5 and the Magneli oxides (Ti4O7 to Ti50O99). In conventional processing of ilmenite with hydrogen it is possible to reduce TiO2 to Ti7O13 within approximately an hour, but with poor utilization of hydrogen on the order of one mole of H2 per thousand. In the cold or non-equilibrium plasma TiO2 can be rapidly reduced to Ti2O3 with hydrogen utilization exceeding 10 percent. Based on design considerations of the plasma reactor greater utilization of the hydrogen in the reduction of TiO2 is possible.

Bullard, D.↗
NASA NTRS2001

Laser Powered Launch Vehicle Performance Analyses

The purpose of this study is to establish the technical ground for modeling the physics of laser powered pulse detonation phenomenon. Laser powered propulsion systems involve complex fluid dynamics, thermodynamics and radiative transfer processes. Successful predictions of the performance of laser powered launch vehicle concepts depend on the sophisticate models that reflects the underlying flow physics including the laser ray tracing the focusing, inverse Bremsstrahlung (IB) effects, finite-rate air chemistry, thermal non-equilibrium, plasma radiation and detonation wave propagation, etc. The proposed work will extend the base-line numerical model to an efficient design analysis tool. The proposed model is suitable for 3-D analysis using parallel computing methods.

Chen, Yen-Sen↗
NASA NTRS1991

Metastable phase formation in undercooled liquid lead alloys

During non-equilibrium processes metastable phases are formed instead of stable phases due to the operation of various kinetic or structural constraints. By removing the most effective nucleation sites for the stable phase using emulsified droplet samples, stable phase formation can be prohibited in a broad range of undercooling and the phase space available in the metastable regime can be mapped out. With this method the thermodynamic properties of the undercooled liquid and the metastable phase boundaries corresponding to reversible metastable (solid + liquid) phase equilibria are examined from experimental data. The analysis reveals important implications for the nucleation kinetics and the choice of the kinetically most favored solidification path.

Fecht, Hans J.↗
NASA NTRS2003

G-jitter Effects on Transport and Pattern Formation

The research performed under this grant has led to an number of new insights into two general categories of fluid flows in the presence of time-dependent acceleration, as outlined briefly below. These results have been widely communicated in the scientific community through seven presentations at international conferences (4 invited, 3 contributed), five published papers (4 journal articles and 1 conference proceeding), and images from the research featured on the cover of all 2003 editions of the research journal, Nonlinearity. The work performed under this proposal also contained a substantial educational component by contributed significantly to the scientific training of one postdoctoral associate, one Ph.D. student and five undergraduate researchers. One main area of focus in this research was convective flow with time-dependent acceleration. Convection is one class of behavior that can arise from g-jitter effects. Our research focused on studies of Rayleigh-Benard system, which is an important model for understanding thermal convection; studies of this problem in the presence of acceleration modulations provided insight into the nature of g-jitter induced flow and of the effects of modulation and noise on non-equilibrium pattern formation. Our experiments on vertically vibrated Rayleigh-Benard convection demonstrated the existence of two classes of pure flow patterns (synchronous & subharmonic) patterns) that had long been predicted by theory but never before observed experimentally. Detailed studies of ranges of parameters where both classes of patterns exist simultaneously led to the discovery of a new type of patterns (called superlattices) in systems driven out of thermodynamic equilibrium.

Schatz, Michael F.↗
NASA NTRS1993

Thermophysical properties of simple liquid metals: A brief review of theory

In this paper, we review the current theory of the thermophysical properties of simple liquid metals. The emphasis is on thermodynamic properties, but we also briefly discuss the nonequilibrium properties of liquid metals. We begin by defining a 'simple liquid metal' as one in which the valence electrons interact only weakly with the ionic cores, so that the interaction can be treated by perturbation theory. We then write down the equilibrium Hamiltonian of a liquid metal as a sum of five terms: the bare ion-ion interaction, the electron-electron interaction, the bare electron-ion interaction, and the kinetic energies of electrons and ions. Since the electron-ion interaction can be treated by perturbation, the electronic part contributes in two ways to the Helmholtz free energy: it gives a density-dependent term which is independent of the arrangement of ions, and it acts to screen the ion-ion interaction, giving rise to effective ion-ion pair potentials which are density-dependent, in general. After sketching the form of a typical pair potential, we briefly enumerate some methods for calculating the ionic distribution function and hence the Helmholtz free energy of the liquid: monte Carlo simulations, molecular dynamics simulations, and thermodynamic perturbation theory. The final result is a general expression for the Helmholtz free energy of the liquid metal. It can be used to calculate a wide range of thermodynamic properties of simple metal liquids, which we enumerate. They include not only a range of thermodynamic coefficients of both metals and alloys, but also many aspects of the phase diagram, including freezing curves of pure elements and phase diagrams of liquid alloys (including liquidus and solidus curves). We briefly mention some key discoveries resulting from previous applications of this method, and point out that the same methods work for other materials not normally considered to be liquid metals (such as colloidal suspensions, in which the suspended microspheres behave like ions screened by the salt solution in which they are suspended). We conclude with a brief discussion of some non-equilibrium (i.e., transport) properties which can be treated by an extension of these methods. These include electrical resistivity, thermal conductivity, viscosity, atomic self-diffusion coefficients, concentration diffusion coefficients in alloys, surface tension and thermal emissivity. Finally, we briefly mention two methods by which the theory might be extended to non-simple liquid metals: these are empirical techniques (i.e., empirical two- and three-body potentials), and numerical many-body approaches. Both may be potentially applicable to extremely complex systems, such as nonstoichiometric liquid semiconductor alloys.

Stroud, David↗
DOE OSTI2025

Potentials of mean force fail to describe chemical bond-breaking in solution

Many liquid phase studies assume that the potential energy surfaces of reacting molecules are the same as in the gas phase, neglecting complex solvent dynamics that can completely alter the nature of chemical reactivity. Even studies that include solvent effects typically only consider them in an average, equilibrium way as part of a potential of mean force (PMF). In this work, we use mixed quantum/classical simulations to compare how equilibrium and non-equilibrium solvent motions affect the photodissociation of a simple diatomic molecule, NaK + , in liquid tetrahydrofuran. A PMF analysis shows that as the excited-state molecule dissociates with the solvent at equilibrium, the bonding electron remains associated with K + at short bond distances but eventually localizes on Na + at the end of dissociation. When we examine non-equilibrium dynamical photodissociation trajectories, however, we find that they fall into three distinct categories: about a quarter of them have the bonding electron mainly associated with Na + , another quarter stay mainly associated with K + , and about half have the bonding electron shared roughly equally between the two ions. The results show that equilibrium PMFs cannot accurately describe the dynamics of bond-breaking chemical reactions in solution because there is insufficient time for the solvent to reach equilibrium on the time scale over which bond dissociation occurs. Furthermore, our analysis shows that the solvent coupling between the electronic energy surfaces is similar at and away from equilibrium, suggesting that other factors, such as solute velocity-driven solvent memory effects, play a more important role in explaining the failure of the equilibrium PMF to predict the non-equilibrium dynamics.

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

Incorporating the molecular-scale into a hydrodynamic description of confined aqueous systems

Hydrodynamics provides a continuum-level description of fluid motion, but its applicability at the nanoscale becomes uncertain due to the emerging importance of molecular-level effects such as spatial heterogeneity. Hydrodynamic boundary conditions that incorporate molecular details allow us to partition the system into a near-wall region and a bulk fluid region. We identify a hydrodynamic wall located inside the fluid that determines where slip begins. By extending the hydrodynamic wall with the slip length, the position of the extrapolated wall is established. This offers a unified description of both slip and stagnant flow behaviors, with wall hydrophobicity characterized by the relative location of the extrapolated wall with respect to the physical wall. Employing this concept in analyses of equilibrium molecular dynamics (MD) and non-equilibrium MD simulations of Couette and Poiseuille flows, our results demonstrate consistency between equilibrium and non-equilibrium approaches across different flow types and confinement levels. This demonstrates the robust nature of linear response theory. We then explore the effects of fluid-wall and bulk fluid interactions on the hydrodynamic properties. These findings enhance the effectiveness of molecular-based simulations for investigating complex confined systems in nanofluidics, biology, and colloidal science, offering a complementary molecular-scale perspective to traditional continuum approaches.

flow simulations↗
DOE OSTI2026

A modified scheil approach for nucleation-dependent solidification pathways

As-solidified microstructures of near-eutectic alloys often contain multiple primary phases that are not expected from equilibrium phase diagrams. Such microstructures are caused by cooling-rate-dependent solidification pathways, a factor not captured by the Scheil–Gulliver model or variations thereof. Here, we present a model and algorithm that incorporate the critical nucleation undercooling for each solid phase into the Scheil–Gulliver model. We hypothesize that the non-equilibrium microstructure formation is primarily governed by a nucleation-competition mechanism. This mechanism accounts for both stable/metastable phase selection and primary-phase formation within eutectic regions driven by asymmetric nucleation barriers. The model is validated against a hypereutectic Al-Fe alloy, where it successfully reproduces the observed microstructural constituents, revealing the key dependencies of solidification microstructure on nucleation kinetics. Applicability to multicomponent systems is demonstrated through a hypereutectic Al–Fe–Si ternary alloy, where the model successfully predicts divorced eutectic microstructures and the associated oscillatory solidification pathways along univariant lines. As a result, the proposed framework establishes a nucleation-dependent computational approach for interpreting and predicting solidification microstructures.

Alloy design↗
DOE OSTI2024

Flame kinetics at scramjet-engine-relevant conditions: Role of prompt dissociation of weakly-bound radicals

Combustion in high-speed ram-based propulsion engines occurs under distinct thermodynamic conditions of high reactant temperatures (greater than 1000 K) and relatively low pressures (<5 atm). There is a lack of fundamental flame measurements at such conditions that result in adiabatic flame temperatures (T ad ) exceeding 2500 K. In this work, we have measured laminar flame speeds of oxygen-enriched CH 4 /oxidizer mixtures at sub-atmospheric conditions to probe kinetics at high T ad using the isobaric spherically expanding flame approach. Simulations with recent kinetic models revealed increasing differences between data and model predictions with increasing T ad , reaching up to 25 %. Kinetic analyses reveal that at the thermodynamic conditions in these O 2 -enriched flames, i.e., lower pressures and higher T ad , the effects of HCO prompt dissociation are accentuated. In addition to HCO, the prompt dissociations of CH 2 OH and C 2 H 5 are also considered. Here, the prompt dissociations of all three radicals were evaluated and their effects considered in flame speed simulations. Reaction path analysis for the present flames revealed that approximately half of the reaction flux for HCO formation undergoes prompt dissociation to H + CO. Furthermore, these analyses also revealed that the pathways and sensitive reactions are similar between oxygen-enriched fuel/oxidizer mixtures and preheated fuel/air mixtures, if both have similar T ad . Thus, flames of oxygen-enriched mixtures could be a surrogate to probe the flame chemistry of highly preheated mixtures at relatively low pressures that are often encountered in ram-based propulsion engine combustors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2007

Large Solar-Rejection Filter

NASA utilized Image Intensified Video Cameras for ATV data acquisition from a jet flying at 12.8 km. Afterwards the video was digitized and then analyzed with a modified commercial software package, Image Systems Trackeye. Astrometric results were limited by saturation, plate scale, and imposed linear plate solution based on field reference stars. Time-dependent fragment angular trajectories, velocities, accelerations, and luminosities were derived in each video segment. It was evident that individual fragments behave differently. Photometric accuracy was insufficient to confidently assess correlations between luminosity and fragment spatial behavior (velocity, deceleration). Use of high resolution digital video cameras in future should remedy this shortcoming. A scenically accurate description of matter interpreted as a substance made up of corpuscular constituents was established during the course of the 19th century. In this description, atoms--the building blocks of the matter--form molecules. The properties of the molecules were described by chemistry or thermodynamics depending on what characteristics of the matter were investigated. In both theories, the molecules can dissociate to atoms when the kinetic energies of the atoms exceed the strength of the chemical bonds. The number of atoms is always preserved in a closed system. This is not true, however, when the matter takes up much higher energies at relativistic scales. New particles can be produced at the expense of the kinetic energy. The number of particles is no longer preserved. There are other conserved quantities, however, these quantities, the charge, baryon number, lepton number, are associated with particles that are considered elementary today. The properties and behavior of these elementary particles is the subject of Particle Physics or High Energy Physics. Practice Page: A crewmember touching a positively charged surface was thought to be galvanically isolated from the vehicle ground analogous to a bird on a high voltage power wire. Recent analysis confirms that positive floating potentials, ionospheric currents to the EVA suit, can be hazardous. The analysis is wrong in that the ionospheric plasma itself can close the circuit. Parametric analysis of very low voltage exposures (2 to 15 volts) could cause pain and/or involuntary muscle tetani or spinal cord shock. NASA worked with the Naval Health Research Center Detachment Directed Energy Bioeffects Laboratory to examine the affects electrical hazards could have on extravehicular activity using two models. The results of the two computational models were combined to predict areas of the body in which neurons of different diameters would be excited. They predicted that physiologically active current could be conducted across the crew member causing catastrophic hazards. Future work to analyze additional current paths was proposed. The FUSE spectrum of BB Dor, observed in a high state, is modeled with an accretion disk with a very low inclination (possibly lower than 10 degrees). Assuming an average WD mass of 0.8 solar mass leads to a distance of the order of approximately 650pc, consistent with the extremely low galactic reddening in its direction, and a mass accretion rate of 10 (exp -9) solar mass a year. The spectrum presents some broad and deep silicon and sulfur absorption lines, indicating that these elements are over-abundant: silicon is 3 times solar, and sulfur is 20 times solar. The FUSE spectrum of BB Dor, observed in a high state, is modeled with an accretion disk with a very low inclination (possibly lower than 10 degrees). Assuming an average WD mass of 0.8 solar mass leads to a distance of the order of approximately 650pc, consistent with the extremely low galactic reddening in its direction, and a mass accretion rate of 10 (exp -9) solar mass a year. The spectrum presents some broad and deep silicon and sulfur absorption lines, indicating that these elements are over-abundant: silicon is times solar, and sulfur is 20 times solar. From its previously observed optical emission lines, P831-57 (WD 0334 6400 or Ret 1 in A Catalog and Atlas of Cataclysmic Variables: Living Edition) has been suspected to contain an accretion disk associated with a companion star in orbit around a subdwarf star with a temperature T is greater than 21,000K. P831-57 has therefore been classified as a nova-like. However, our present observations show it to be a DA + dMe binary. The analysis of its Far Ultraviolet Spectroscopic Explorer (FUSE) spectrum (continuum and lines) reveal an average mass white dwarf (Log(g) approximately equals 7.8 plus or minus 0.1)with a temperature T approximately equals 37,000 plus or minus 500K, an extremely low projected rotational velocity, and a distance of about 115 plus or minus 5pc. The photosphere contains C, N, Si, and S (at about 1% of solar abundances). The dMe star is seen as a flux excess in near-infrared photometry and appears to show occasional flaring of about one magnitude as seen in the Harvard plates. There is no evidence of periodic variability in the spectroscopic or photometric data. We find no evidence of a an accretion disk, instead we find evidence of wind accretion as the stellar carbon abundance (N(C)/N(H) = 2.5 x 10 (exp -6) is about ten times larger than predicted by radiative levitation for such a gravity and temperature. The power needs and solutions for the space exploration and lunar mobility program are discussed. Long term missions in space and on the lunar surface require high energy batteries. Rechargeable batteries for mobility systems and portable utility pallet are needed for successful exploration missions. Nanomaterial usage increases the energy density of the cells apart from increasing the power density. The symptoms and threats from acute mountain sickness (AMS) are discussed. The underlying assumptions concerning spacecraft atmosphere mean there is a potential risk to astronauts. The baseline worst case scenario is about 25% based on direct ascent to 8.0 psia with 32% O2. The practical mitigation approach is a staged depressurization scheme. Due to the uncertainty about potential AMS risk the flight surgeon should prepare, focused research should proceed, and current analytical efforts should continue. This slide presentation reviews the threat astronauts face from acute mountain sickness (AMS). The incidence of AMS is highly variable and the risk to astronauts is unclear. However, the potential risk is higher than current EAAs suggest and precautions should be implemented. Due to the uncertain nature of the risk flight surgeons should prepare, and focused research should proceed. Quasiparticle tunneling can be understood at low temperatures using a non-equilibrium kinetic theory. Low-temperature tunneling can be suppressed using existing techniques: island-lead gap engineering, lead cooling with SIN junctions, and quasiparticle traps. The slide presentation reviews the science goals and mission planning for the proposed Jupiter System Observer (JSO). The spacecraft would accommodate 9 instruments making the instrument payload more capable than any previously flown to Jupiter. The would-be mission includes study programs for Jupiter's satellites, interior, magnetosphere and atmosphere. The power needs and solutions for the space exploration and lunar mobility program are discussed. Long term missions in space and on the lunar surface require high energy batteries. Rechargeable batteries for mobility systems and portable utility pallet are needed for successful exploration missions. Nanomaterial usage increases the energy density of the cells apart from increasing the power density. This slide presentation reviews the power requirements for the space exploration and the lunar surface mobility programs. It includes information about the specifications for high energy batteries and the power requirements for lunar rovers,

Roberts, William↗
DOE OSTI2025

Design of novel refractory equiatomic multi-principal elemental alloys based on Mo-Nb-Ti system for Gen IV reactor applications

Excellent irradiation damage resistance demonstrated by multi-principal elemental alloys (MPEAs) has sparked significant interest among researchers, prompting exploration into their vast compositional space, to validate their suitability for nuclear applications. A combined approach of thermodynamic and empirical parameters calculations alongside CALPHAD (CALculation of PHAse Diagrams) for phase formation predictions enable high-throughput material selection for sophisticated applications like nuclear, overcoming laborious and time-consuming experiments. Key thermodynamic and empirical parameters for eight novel equiatomic MPEAs, based on seven low thermal neutron cross section refractory elements, for predicting phase formation were calculated, and equilibrium and non-equilibrium simulations in CALPHAD were employed to comprehensively model the systems. Pseudo binary phase diagram simulations showed that Zr, V or equiatomic CrV additions to the base MoNbTi alloy (MoNbTi-Zr, MoNbTi-V and MoNbTi-CrV alloys) favor the formation of isomorphous body-centered cubic (BCC) phase at high temperatures, while Cr, Al, equiatomic ZrV, or equiatomic CrAl additions (MoNbTi-Cr, MoNbTi-Al, MoNbTi-ZrV or MoNbTi-CrAl alloys) limit the solubility of them. Equilibrium CALPHAD simulations at 750 oC were consistent with XRD results on MoNbTi, MoNbTiZr and MoNbTiCr alloys, and partially for others. Notably, elemental segregation observed in the backscattered electron (BSE) scanning electron microscopy (SEM) images of the alloys was accurately simulated through non-equilibrium Scheil solidification calculations in CALPHAD, further verified by experiments. The precipitation of TiCr2 Laves phase in Cr containing MoNbTiCr and MoNbTiCrAl was accurately predicted while discrepancies were noted in MoNbTiCrV. The equilibrium simulations also provided insights into phase compositions at specific temperatures offering a pathway for tailoring the desired microstructure and properties of these systems. Empirical parameters calculations successfully predicted random solid solution in the base MoNbTi alloy, and with an exception in MoNbTiV and MoNbTiAl, predicted intermetallic precipitation in the rest, especially, Laves phase precipitation in Cr containing alloys.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

Pseudo-equilibrium theory for extrinsic doping control of the topological semimetal Cd 3 As 2

The standard approach for predicting defect equilibria from first principles assumes that the solid-state system is initially in a thermodynamic equilibrium with the external atomic reservoirs. This “growth step” is then often followed by a temperature quench in a “pseudo-equilibrium” in which some or all defect concentrations are frozen in until only the Fermi level E F remains to be equilibrated. However, this protocol does not account for the possibility of site exchanges which can create important defect redistributions as long as short-range defect migration is kinetically permissible. To model this redistribution, we developed an approach to solve for the non-equilibrium chemical potentials as a function of temperature while maintaining the overall defect stoichiometry. We then apply this approach to the Dirac semimetal Cd 3 As 2 to model extrinsic doping with group 1/11 and 14 elements. Undoped Cd 3 As 2 exhibits an undesirable mismatch between E F and the Dirac point. This unintentional electron doping originates from intrinsic defects and is difficult to overcome through adjustment of synthesis conditions alone. Employing our pseudo-equilibrium modeling, we identify extrinsic doping strategies for realizing doping-balanced Cd 3 As 2 at the relatively low temperatures accessible in thin-film growth of this material.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Thermodynamic consistency and fluctuations in mesoscopic stochastic simulations of reactive gas mixtures

It is essential that mesoscopic simulations of reactive systems reproduce the correct statistical distributions at thermodynamic equilibrium. By considering a compressible fluctuating hydrodynamics (FHD) simulation method of ideal gas mixtures undergoing reversible reactions described by the chemical Langevin equations, we show that thermodynamic consistency in reaction rates and the use of instantaneous temperatures for the evaluation of reaction rates is required for fluctuations for the overall system to be correct. We then formulate the required properties of a thermodynamically consistent reaction (TCR) model. As noted in the literature, while reactions are often discussed in terms of forward and reverse rates, these rates should not be modeled independently because they must be compatible with thermodynamic equilibrium for the system. Using a simple TCR model where each chemical species has constant heat capacity, we derive the explicit condition that the forward and reverse reaction rate constants must satisfy in order for the system to be thermodynamically consistent. We perform equilibrium and non-equilibrium simulations of ideal gas mixtures undergoing a reversible dimerization reaction to measure the fluctuational behavior of the system numerically. We confirm that FHD simulations with the TCR model give the correct static structure factor of equilibrium fluctuations. For the statistically steady simulation of a gas mixture between two isothermal walls with different temperatures, we show using the TCR model that the temperature variance agrees with the corresponding thermodynamic-equilibrium temperature variance in the interior of the system, whereas noticeable deviations are present in regions near walls, where chemistry is far from equilibrium.

Polimeno, Matteo [University of California, Merced↗
DOE OSTI2026

Spectral deconvolution without the deconvolution: Extracting temperature from x-ray Thomson scattering spectra without the source-and-instrument function

X-ray Thomson scattering (XRTS) probes the dynamic structure factor of the system, but the measured spectrum is broadened by the combined source-and-instrument function (SIF) of the setup. In order to extract properties such as temperature from an XRTS spectrum, the broadening by the SIF needs to be removed. Recent work [Dornheim et al. Nat. Commun. 13 , 7911 (2022)] has suggested that the SIF may be deconvolved using the two-sided Laplace transform. However, the extracted information can depend strongly on the shape of the input SIF, and the SIF is in practice challenging to measure accurately. Here, we propose an alternative approach: we demonstrate that considering ratios of Laplace-transformed XRTS spectra collected at different scattering angles is equivalent to performing the deconvolution, but without the need for explicit knowledge of the SIF. From these ratios, it is possible to directly extract the temperature from the scattering spectra, when the system is in thermal equilibrium. We find the method to be generally robust to spectral noise and physical differences between the spectrometers, and we explore situations in which the method breaks down. Furthermore, the fact that consistent temperatures can be extracted for systems in thermal equilibrium indicates that non-equilibrium effects could be identified by inconsistent temperatures of a few eV between the ratios of three or more scattering angles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1996

Nonlinear Dynamics and Nucleation Kinetics in Near-Critical Liquids

The objective of our study is to model the nonlinear behavior of a near-critical liquid following a rapid change of the temperature and/or other thermodynamic parameters (pressure, external electric or gravitational field). The thermodynamic critical point is manifested by large, strongly correlated fluctuations of the order parameter (particle density in liquid-gas systems, concentration in binary solutions) in the critical range of scales. The largest critical length scale is the correlation radius r(sub c). According to the scaling theory, r(sub c) increases as r(sub c) = r(sub 0)epsilon(exp -alpha) when the nondimensional distance epsilon = (T - T(sub c))/T(sub c) to the critical point decreases. The normal gravity alters the nature of correlated long-range fluctuations when one reaches epsilon approximately equal to 10(exp -5), and correspondingly the relaxation time, tau(r(sub c)), is approximately equal to 10(exp -3) seconds; this time is short when compared to the typical experimental time. Close to the critical point, a rapid, relatively small temperature change may perturb the thermodynamic equilibrium on many scales. The critical fluctuations have a hierarchical structure, and the relaxation involves many length and time scales. Above the critical point, in the one-phase region, we consider the relaxation of the liquid following a sudden temperature change that simultaneously violates the equilibrium on many scales. Below T(sub c), a non-equilibrium state may include a distribution of small scale phase droplets; we consider the relaxation of such a droplet following a temperature change that has made the phase of the matrix stable.

Patashinski, Alexander Z.↗