SEARCH · Search NASA

Results for “Entropy Production”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Entanglement and correlations in fast collective neutrino flavor oscillations

Collective neutrino oscillations play a crucial role in transporting lepton flavor in astrophysical settings like supernovae and neutron star binary merger remnants, which are characterized by large neutrino densities. In these settings, simulations in the mean-field approximation show that neutrino-neutrino interactions can overtake vacuum oscillations and give rise to fast collective flavor evolution on time-scales t ∝ μ –1 , with μ proportional to the local neutrino density. In this work, we study the full out-of-equilibrium flavor dynamics in simple multi-angle geometries displaying fast oscillations in the mean field linear stability analysis. Focusing on simple initial conditions, we analyze the production of pair correlations and entanglement in the complete many-body-dynamics as a function of the number N of neutrinos in the system, for up to thousands of neutrinos. Similarly to simpler geometries with only two neutrino beams, we identify three regimes: stable configurations with vanishing flavor oscillations, marginally unstable configurations with evolution occurring on long time scales τ ≈ µ –1 $\sqrt{N}$, and unstable configurations showing flavor evolution on short time scales τ ≈ μ –1 log(N). Here, we present evidence that these fast collective modes are generated by the same dynamical phase transition which leads to the slow bipolar oscillations, establishing a connection between these two phenomena and explaining the difference in their time scales. We conclude by discussing a semi-classical approximation which reproduces the entanglement entropy at short to medium time scales and can be potentially useful in situations with more complicated geometries where classical simulation methods starts to become inefficient.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Finite-temperature many-body perturbation theory for anharmonic vibrations: Recursions, algebraic reduction, second-quantized reduction, diagrammatic rules, linked-diagram theorem, finite-temperature self-consistent field, and general-order algorithm

A unified theory is presented for finite-temperature many-body perturbation expansions of the anharmonic vibrational contributions to thermodynamic functions, i.e., the free energy, internal energy, and entropy. The theory is diagrammatically size-consistent at any order, as ensured by the linked-diagram theorem proved in this study, and, thus, applicable to molecular gases and solids on an equal footing. It is also a basis-set-free formalism, just like its underlying Bose–Einstein theory, capable of summing anharmonic effects over an infinite number of states analytically. It is formulated by the Rayleigh–Schrödinger-style recursions, generating sum-over-states formulas for the perturbation series, which unambiguously converges at the finite-temperature vibrational full-configuration-interaction limits. Two strategies are introduced to reduce these sum-over-states formulas into compact sum-over-modes analytical formulas. One is a purely algebraic method that factorizes each many-mode thermal average into a product of one-mode thermal averages, which are then evaluated by the thermal Born–Huang rules. Canonical forms of these rules are proposed, dramatically expediting the reduction process. The other is finite-temperature normal-ordered second quantization, which is fully developed in this study, including a proof of thermal Wick’s theorem and the derivation of a normal-ordered vibrational Hamiltonian at finite temperature. The latter naturally defines a finite-temperature extension of size-extensive vibrational self-consistent field theory. These reduced formulas can be represented graphically as Feynman diagrams with resolvent lines, which include anomalous and renormalization diagrams. Two order-by-order and one general-order algorithms of computing these perturbation corrections are implemented and applied up to the eighth order. The results show no signs of Kohn–Luttinger-type nonconvergence.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2023

Effect of a micro-scale dislocation pileup on the atomic-scale multi-variant phase transformation and twinning

In this paper, we perform concurrent atomistic–continuum (CAC) simulations to assess the contribution of the internal stress induced by the microscale dislocation pileup at an atomically structured interface to the atomic-scale phase transformations (PTs), reverse PTs, and twinning. The main novelty of this work is to unify the atomistic description of the interface and the coarse-grained (CG) description of the lagging dislocations away from the interface within one single framework. Our major findings are: (a) the interface dynamically responds to a pileup by forming steps/ledges, the height of which is proportional to the number of dislocations arriving at the interface; (b) the pileup-induced internal stress concentration profile follows neither the classical Eshelby model nor the super-dislocation model alone, but a combination of them; (c) when the pre-sheared sample is compressed, a direct square-to-hexagonal PT occurs ahead of the pileup tip and eventually grows into a wedge shape. The two variants of the hexagonal phases form a twin with respect to each other; (d) upon a further increase of the loading, part of the newly formed hexagonal phase transforms back to the square phase. The square product phase resulting from this reverse PT forms a twin with respect to the initial square phase. All phase boundaries (PBs) and twin boundaries (TBs) are stationary and correspond to zero thermodynamic Eshelby driving forces; and (e) the microscale dislocation pileup-induced internal shear stress and the structural change at the atomic-scale interface reduces the stress required for initiating a PT by a factor of 5.5, comparing with that in the sample containing no dislocations. This work is the first characterization of the behavior of PTs/twinning resulting from the reaction between a microscale dislocation slip and an atomically structured interface. In conclusion, the gained knowledge will advance our understanding of how the multi-phase material behaves in many complex physical processes, such as the synthesis of multi-phase high-entropy alloys or superhard ceramics under high-pressure torsion, deep mantle earthquakes in geophysics, and so on, which all involve dislocation slip, PTs, twinning, and their interactions across from the atomistic to the microscale and beyond.

36 MATERIALS SCIENCE↗

DOE OSTI2026

A Robust Methodology to Elucidate Kinetics of Room Temperature Electrochemical Propane Adsorption on Platinum

Electrocatalytic activation of alkanes can further decarbonize chemical manufacturing by leveraging affordable renewable electricity and readily available shale gas reserves in the United States. Earlier works have identified the unique role of Pt in adsorbing and activating alkanes, like propane, at room temperature in acidic, aqueous electrolytes, revealing spontaneous formation of deeply dehydrogenated propane-derived surface species with an intact C 3 - backbone. Although an adsorption mechanism was hypothesized, it has not been explicitly investigated to date, preventing the quantification of kinetic rate parameters. A robust methodology to investigate and benchmark propane adsorption kinetics on Pt is critical for the rational design of electrocatalysts that exhibit higher selectivity toward desired partially oxidized products. Herein, we analyze an oxidative current transience that appears during the adsorption of propane on Pt in aqueous electrochemical conditions and develop a methodology that elucidates the adsorption mechanism and enables quantification of rate parameters such as order dependences and apparent activation barriers. This method yields an expected first-order dependence with respect to propane concentration at low coverage and reveals a second-order dependence with respect to the concentration of surface active sites. Additionally, the apparent activation barrier for propane adsorption was calculated using an Arrhenius analysis of the current transience under temperature control. The experimentally measured activation barrier of 35 kJ mol –1 is in excellent agreement with the theoretical barrier calculated by density functional theory (DFT). The kinetic analysis was extended, via the use of transition state theory, to extract entropy and enthalpy of activation, yielding consistent results with the proposed two-step adsorption mechanism and DFT calculations. These results demonstrate reliable quantification of kinetic parameters for electrocatalytic activation of C–H bonds in alkanes that can be employed for rational catalyst development for a versatile range of electrocatalytic conditions.

alkane activation↗

DOE OSTI2021

Protective effects of non–catalytic proteins on endoglucanase activity at air and lignin interfaces

The manner in which added non-catalytic proteins during enzymatic hydrolysis of lignocellulosic substrates enhances hydrolysis mechanisms is not completely understood. Prior research has indicated that a reduction in the non-specific adsorption of enzymes on lignin, and deactivation of enzymes exposed to air–liquid interface provide rationale. This work investigated root causes including effects of the air–liquid interface on non-catalytic proteins, and effects of lignin on endoglucanase. Three different experimental designs and three variables (air–liquid interfacial area, the types of lignin (acid or enzymatic lignin), and the presence of non-enzymatic protein (bovine serum albumin [BSA] or soy proteins ) were used. The results showed that acid isolated lignin adsorbed almost all endoglucanase activity initially present in supernatant, independent of air interface conditions (25 or 250 ml flasks) with the presence of BSA preventing this effect. Endoglucanase lost 30%–50% of its activity due to an air–liquid interface in the presence of lignin while addition of non-enzymatic protein helped to preserve this enzyme's activity. Langmuir and Freundlich models applied to experimental data indicated that the adsorption increases with increasing temperature for both endoglucanase and BSA. Adsorption of the enzyme and protein were endothermic with an increase in entropy. Finally, these results, combined, show that hydrophobicity plays a strong role in the adsorption of both endoglucanase and BSA on lignin.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS2010

Context Modeler for Wavelet Compression of Spectral Hyperspectral Images

A context-modeling sub-algorithm has been developed as part of an algorithm that effects three-dimensional (3D) wavelet-based compression of hyperspectral image data. The context-modeling subalgorithm, hereafter denoted the context modeler, provides estimates of probability distributions of wavelet-transformed data being encoded. These estimates are utilized by an entropy coding subalgorithm that is another major component of the compression algorithm. The estimates make it possible to compress the image data more effectively than would otherwise be possible. The following background discussion is prerequisite to a meaningful summary of the context modeler. This discussion is presented relative to ICER-3D, which is the name attached to a particular compression algorithm and the software that implements it. The ICER-3D software is summarized briefly in the preceding article, ICER-3D Hyperspectral Image Compression Software (NPO-43238). Some aspects of this algorithm were previously described, in a slightly more general context than the ICER-3D software, in "Improving 3D Wavelet-Based Compression of Hyperspectral Images" (NPO-41381), NASA Tech Briefs, Vol. 33, No. 3 (March 2009), page 7a. In turn, ICER-3D is a product of generalization of ICER, another previously reported algorithm and computer program that can perform both lossless and lossy wavelet-based compression and decompression of gray-scale-image data. In ICER-3D, hyperspectral image data are decomposed using a 3D discrete wavelet transform (DWT). Following wavelet decomposition, mean values are subtracted from spatial planes of spatially low-pass subbands prior to encoding. The resulting data are converted to sign-magnitude form and compressed. In ICER-3D, compression is progressive, in that compressed information is ordered so that as more of the compressed data stream is received, successive reconstructions of the hyperspectral image data are of successively higher overall fidelity.

Kiely, Aaron↗

NASA NTRS1976

Noise produced by fluid inhomogeneities

In an actual engine system, one mechanism for production of hot spots is the burning of various fluid elements at various different mixture ratios. Variable mixture ratio means variable temperature, and this is the effect that has been studied insofar as a noise source is concerned. However, variable mixture ratio also implies variable molecular weight and heat capacity. The paper investigates whether either of these last two variations may be responsible for a significant noise source. The analysis is made within the context of one-dimensional unsteady flow as in the work of Candel (1972). A mixture of thermally perfect gases is assumed as the working fluid, and the fluid composition consists of species 1 and a small and variable mole fraction of species 2 which has different molecular weight and specific heats as compared with species 1. In the absence of changes in the ratio of specific heat (gamma), the entropy variations due to temperature and molecular weight variations are equivalent as a sound source. The portion of sound called 'gamma prime' noise is discussed.

Strahle, W. C.↗

DOE OSTI2023

Molecular-Weight-Dependent Infiltration and Adsorption of Polymers into Nanochannels

The mesoporous silica shell coating hydrogenolysis nano-catalysts alters the molecular weight distributions of cleaved polymer chains compared to catalysts without a shell. The shell, composed of radially aligned narrow cylindrical nanopores, reduces the formation of low-valued gaseous products and increases the median molecular weight of the product, thus enhancing the value of the products for polymer upcycling. Here, to understand the role of the mesoporous shell, we have studied the spatial distribution of polystyrene chains, used as a model polymer, in the nanochannels in both the melt phase and solution phase. In the melt, we observed from small-angle X-ray scattering experiments that the infiltration rate of the polymer into the nanochannels is inversely proportional to the molecular weight, which is consistent with theory. In theta solution experiments using UV–vis spectroscopy, we found that the shell significantly enhances polymer adsorption compared to nanoparticles without pores. In addition, the degree of polymer adsorption is not a monotonic function of molecular weight but initially increases with the molecular weight before eventually decreasing. The molecular weight for the peak adsorption increases with the pore diameter. This adsorption behavior is rationalized as resulting from a balance between the mixing entropy gain by surface adsorption and the conformational entropy penalty incurred by chains confined in the nanochannels. The spatial distribution of polymer chains in the nanochannels is visualized by energy-dispersive X-ray spectroscopy (EDX), and inverse Abel-transformed data reveals a less uniform polymer distribution along the primary pore axis for longer chains.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2021

Maximum entropy modeling to identify physical drivers of shallow snowpack heterogeneity using unpiloted aerial system (UAS) lidar

Understanding the spatial variability of the snowpack is valuable for hydrologists and ecologists seeking to predict hydrological processes in a cold region. Snow distribution is a function of interactions among static variables, such as terrain, vegetation, and soil properties, and dynamic meteorological variables, such as solar radiation, wind speed and direction, and soil moisture. However, identifying the dominant physical drivers responsible for spatial patterns of the snowpack, particularly for ephemeral, shallow snowpacks, has been challenging due to the lack of the high-resolution snowpack and physical variables with high vertical accuracy as well as inherent limitations in traditional approaches. This study uses an Unpiloted Aerial System (UAS) lidar-based snow depth and static variables (1-m spatial resolution) to analyze field-scale spatial structures of snow depth and apply the Maximum Entropy (MaxEnt) model to identify primary controls over open terrain and forests at the University of New Hampshire Thompson Farm Research Observatory, New Hampshire, United States. We found that, among nine topographic and soil variables, plant functional type and terrain roughness contribute up to 80% and 76% of relative importance in the MaxEnt framework to predict locations of deeper or shallower snowpacks, respectively, across a mixed temperate forested and field landscape. Soil variables, such as organic matter and saturated hydraulic conductivity, were also important controls (up to 70% and 81%) on snow depth spatial variations for both open and forested landscapes suggesting spatial variations in soil variables under snow can control thermal transfer among soil, snowpack, and surface-atmosphere. This work contributes to improving land surface and snow models by informing parameterization of the sub-grid scale snow depths, down-scaling remotely sensed snow products, and understanding field scale snow states.

Eunsang Cho↗

NASA NTRS2022

Bryce Canyon Water Resources: Monitoring Vegetation Health and Water Availability in Bryce Canyon National Park for Drought Stress Mitigation Planning

Bryce Canyon National Park is home to groundwater-dependent ecosystems (GDEs) that are threatened by a multidecadal drought and increased groundwater extraction due to a spike in tourism. These ecosystems contain unique species that are only found in areas where near-surface groundwater is present, such as aspen groves and fens. These species contribute to the high biodiversity found in Bryce Canyon, which boosts an ecosystem’s productivity and the services it provides to the park. Unfortunately, many of these GDEs are too small to identify with traditional Earth observation platforms and are difficult to physically reach for monitoring purposes. This project partnered with the National Park Service to identify springs and seeps as a proxy for GDEs within Bryce Canyon from 2013–2022. Furthermore, this project tested the feasibility of various methods to detect and monitor springs and seeps and therefore facilitate the partner’s efforts to conserve these ecologically valuable GDEs in Bryce Canyon. The team mapped groundwater discharge with high resolution National Agriculture Imagery Program (NAIP) and assessed park vegetation trends with Landsat 8 Operational Land Imager (OLI) and PlanetScope imagery. In-situ precipitation data and the Western Land Data Assimilation System (WLDAS) were used to produce time series of climatic variables. Seeps and spring locations were predicted using random forest classification and maximum entropy machine learning models.

Groundwater dependent ecosystems↗

DOE OSTI2022

Reversible Ligand Protonation in Noninnocent Constrained-Geometry-Like Group 4 Complexes

The proligands C 5 H 5 CMe 2 CHPhOx R (R = Me 2 , CHMe 2 , and CMe 3 ) react with M(NMe 2 ) 4 (M = Ti, Zr, and Hf) to give monodeprotonated (Ox R CHPhCMe 2 C 5 H 4 )M(NMe 2 ) 3 , doubly deprotonated (C 5 H 4 CMe 2 CPhOx R )M(NMe 2 ) 2 , or a mixture of both. The observed products depend on reaction conditions, the oxazoline substituent, and the metal center, with 4-t-butyl-oxazoline or Ti giving only the monodeprotonated (Ox R CHPhCMe 2 C 5 H 4 )M(NMe 2 ) 3 . Amine elimination from (Ox R CHPhCMe 2 C 5 H 4 )M(NMe 2 ) 3 to give (C 5 H 4 CMe 2 CPhOx R )M(NMe 2 ) 2 is reversible for 4,4-dimethyl- and 4-isopropyl-oxazoline-based ligands in Zr or Hf complexes. Temperature-dependent kinetic studies of the equilibration of (C 5 H 4 CMe 2 CPhOx Me2 )Zr(NMe 2 ) 2 , HNMe 2 and (Ox Me2 CHPhCMe 2 C 5 H 4 )Zr(NMe 2 ) 3 provide the experimental thermodynamic parameters ΔS° = -17.4 ± 2.6 cal·mol –1 K –1 and ΔH° = -6.8 ± 0.8 kcal·mol –1 . An Eyring plot of the rate constants, determined from the system as it approaches equilibrium, gives the activation entropy and activation enthalpy for the addition of (C 5 H 4 CMe 2 CPhOx Me2 )Zr(NMe 2 ) 2 and HNMe 2 as -36 ± 4 cal·mol –1 K –1 and 8.9 ± 1.3 kcal·mol –1 , respectively; the elimination of HNMe 2 from (Ox Me2 CHPhCMe 2 C 5 H 4 )Zr(NMe 2 ) 3 is characterized by ΔS ‡ = -19 ± 5 cal·mol –1 K –1 and ΔH ‡ = 15.7 ± 1.5 kcal·mol –1 . DFT computational models indicate a single-step, nonlinear transfer of the H between the benzylic position of the noninnocent, oxazoline-coordinated ligand and NMe 2 . Computations also confirm the negative activation entropy and the trends in the barriers support the experimental results. Together, these studies indicate the importance of steric effects from the oxazoline ligand, metal center, ancillary ligands, and leaving group on the shuttling of the proton between HNMe 2 and the noninnocent ligand. Finally, these effects suggest that coordination of oxazoline to the metal center is a key part of the benzylic deprotonation and noninnocent behavior of the cyclopentadienyl-oxazoline ligand.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Spinel-Structured High Entropy Oxides: Low Temperature Synthesis, Characterization, and Potential Applications

High entropy oxides (HEOs) have recently attracted increasing attention due to their remarkable properties and relatively low cost. Herein we report a simple, highly scalable, and low temperature method for synthesizing spinel (FeNiCoCuZn) 3 O 4 (HEO). We heated an aqueous solution containing divalent cations in high alkali environments to temperatures of 25 - 95 °C for 24 h under atmospheric pressure. The HEO, synthesized at 95 °C for 24 h in 1 M KOH, was paramagnetic at room temperature, with a magnetic mass susceptibility of χ = (7.5 ± 0.2) × 10 –7 m 3 ·kg –1 . It demonstrated stable electrochemical lithium storage performance, with a gravimetric capacity of ∼300 mAh·g –1 at 100 mA·g –1 . It was also active in the electrocatalytic oxygen evolution reaction with an overpotential of 460 mV in alkaline media. The band gap energies were in the range of 2.4 eV. Our advancement in the synthesis and processing of transition metal-based HEOs will undoubtedly render them a viable solution for next generation materials for energy production and storage.

diffraction↗

NASA NTRS2005

Control of Combustion-Instabilities Through Various Passive Devices

It is well known that under some operating conditions, rocket engines (using solid or liquid fuels) exhibit unstable modes of operation that can lead to engine malfunction and shutdown. The sources of these instabilities are diverse and are dependent on fuel, chamber geometry and various upstream sources such as pumps, valves and injection mechanism. It is believed that combustion-acoustic instabilities occur when the acoustic energy increase due to the unsteady heat release of the flame is greater than the losses of acoustic energy from the system [1, 2]. Giammar and Putnam [3] performed a comprehensive study of noise generated by gasfired industrial burners and made several key observations; flow noise was sometimes more intense than combustion roar, which tended to have a characteristic frequency spectrum. Turbulence was amplified by the flame. The noise power varied directly with combustion intensity and also with the product of pressure drop and heat release rate. Karchmer [4] correlated the noise emitted from a turbofan jet engine with that in the combustion chamber. This is important, since it quantified how much of the noise from an engine originates in the combustor. A physical interpretation of the interchange of energy between sound waves and unsteady heat release rates was given by Rayleigh [5] for inviscid, linear perturbations. Bloxidge et al [6] extended Rayleigh s criterion to describe the interaction of unsteady combustion with one-dimensional acoustic waves in a duct. Solutions to the mass, momentum and energy conservation equations in the pre- and post-flame zones were matched by making several assumptions about the combustion process. They concluded that changes in boundary conditions affect the energy balance of acoustic waves in the combustor. Abouseif et al [7] also solved the one-dimensional flow equations, but they used a onestep reaction to evaluate the unsteady heat release rate by relating it to temperature and velocity perturbations. Their analysis showed that oscillations arise from coupling between entropy waves produced at the flame and pressure waves originating from the nozzle. Yang and Culick [8] assumed a thin flame sheet, which is distorted by velocity and pressure oscillations. Conservation equations were expressed in integral form and solutions for the acoustic wave equations and complex frequencies were obtained. The imaginary part of the frequency indicated stability regions of the flame. Activation energy asymptotics together with a one-step reaction were used by McIntosh [9] to study the effects of acoustic forcing and feedback on unsteady, one-dimensional flames. He found that the flame stability was altered by the upstream acoustic feedback. Shyy et al [10] used a high-accuracy TVD scheme to simulate unsteady, one-dimensional longitudinal, combustion instabilities. However, numerical diffusion was not completely eliminated. Recently, Prasad [11] investigated numerically the interactions of pressure perturbations with premixed flames. He used complex chemistry to study responses of pressure perturbations in one-dimensional combustors. His results indicated that reflected and transmitted waves differed significantly from incident waves.

Frendi, Kader↗

DOE OSTI2021

Transformative Materials for High-Efficiency Thermochemical Production of Solar Fuels (Final Report)

Metal oxide-based two-step solar thermochemical (STC) H 2 O and CO 2 -splitting cycles are a promising route to convert solar thermal energy into fuels. The metal oxide materials are reduced at high temperatures (Step 1), and then at low (but still elevated) temperatures, the reduced oxide is used to split H 2 O or CO 2 (Step 2). However, current applications of these cycles are limited by the efficiency of the metal oxide materials. A lower temperature for reduction is desirable, but that brings a concomitant reduction in the driving force for gas splitting. So, designing novel, high-efficiency materials is challenging. Here, we devised and performed a joint computational-experimental project, combined with materials design strategies and high-throughput approaches with the goal to quickly discover and demonstrate novel thermochemical materials with superior properties. Our approach involves an active feedback loop between experimental and computational research. We have use our previously developed materials design map, to efficiently screen results from high-throughput first- principles computation to predict new compositions and subsequently to experimentally study the properties of novel, predicted materials. Experimental measurements of the thermodynamic properties of selected materials serve as a critical benchmark, not only in evaluating the STCH properties of synthesized compounds, but also serving as a validation dataset for computational approaches. We experimentally explored a set of predicted ABO3 perovskites and focus on obtaining high quality thermodynamic properties for validation of computational prediction of enthalpy and entropy of reduction. These thermodynamic quantities play a major role in designing materials with reduced temperatures of reduction but sufficient gas-splitting rates. Finding good validation between experimental and computational data for enthalpies of reduction, we turn to promising ABO3 materials modified by A and B site substitutions, opening an enormous combinatorial space of materials. We use our high-throughput approach to “tune in” the desired solar thermochemical (STC) properties for STC. This vast composition space can only be reasonably explored using the high-throughput approaches, both computational and experimental, of this proposal. Several promising novel oxide materials have been predicted, synthesized, and verified using our materials design approach.

08 HYDROGEN↗

DOE OSTI2022

Impact of multiplicity fluctuations on entropy scaling across system size

The initial state is one of the greatest uncertainties in heavy-ion collisions. A model-agnostic approach is taken in the phenomenological Trento framework which constrains parameters using Bayesian analysis. However, the color-glass condensate (CGC) effective theory predicts initial energy densities that lie outside the recent Bayesian analyses due, in part, to the assumption in Trento of event-by-event multiplicity fluctuations following a Γ distribution. We compare the Trento-preferred √T A T B scaling to CGC-like T A T B scaling coupled with log-normal fluctuations in AuAu and dAu collisions and find there is a significant impact on the multiplicity distributions and on the eccentricities, which may affect the extraction of viscosity in small systems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Evolution of oxygen and stratification and their relationship in the North Pacific Ocean in CMIP6 Earth system models

Abstract. This study examines the linkages between the upper-ocean (0–200 m) oxygen (O2) content and stratification in the North Pacific Ocean using four Earth system models (ESMs), an ocean hindcast simulation, and an ocean reanalysis. The trends and variability in oceanic O2 content are driven by the imbalance between physical supply and biological demand. Physical supply is primarily controlled by ocean ventilation, which is responsible for the transport of O2-rich surface waters to the subsurface. Isopycnic potential vorticity (IPV), a quasi-conservative tracer proportional to density stratification that can be evaluated from temperature and salinity measurements, is used herein as a dynamical proxy for ocean ventilation. The predictability potential of the IPV field is evaluated through its information entropy. The results highlight a strong O2–IPV connection and somewhat higher (as compared to the rest of the basin) predictability potential for IPV across the tropical Pacific, where the El Niño–Southern Oscillation occurs. This pattern of higher predictability and strong anticorrelation between O2 and stratification is robust across multiple models and datasets. In contrast, IPV at mid-latitudes has low predictability potential and its center of action differs from that of O2. In addition, the locations of extreme events or hotspots may or may not differ between the two fields, with a strong model dependency, which persists in future projections. On the one hand, these results suggest that it may be possible to monitor ocean O2 in the tropical Pacific based on a few observational sites co-located with the more abundant IPV measurements; on the other, they lead us to question the robustness of the IPV–O2 relationship in the extratropics. The proposed framework helps to characterize and interpret O2 variability in relation to physical variability and may be especially useful in the analysis of new observation-based data products derived from the BGC-Argo float array in combination with the traditional but far more abundant Argo data.

Novi, Lyuba↗

DOE OSTI2022

Replica symmetry breaking for the integrable two-site Sachdev–Ye–Kitaev model

We analyze a two-body non-Hermitian two-site Sachdev–Ye–Kitaev (SYK) model with the couplings of one site complex conjugated to the other site. This model, with no explicit coupling between the sites, shows an infinite number of second-order phase transitions, which is a consequence of the factorization of the partition function into a product over Matsubara frequencies. We calculate the quenched free energy in two different ways: first in terms of the single-particle energies and second by solving the Schwinger–Dyson equations of the two-site model. The first calculation can be done entirely in terms of a one-site model. The conjugate replica enters due to non-analyticities when Matsubara frequencies enter the spectral support of the coupling matrix. The second calculation is based on the replica trick of the two-site partition function. Both methods give the same result. The free-fermion partition function can also be rephrased as a matrix model for the coupling matrix. Up to minor details, this model is the random matrix model that describes the chiral phase transition of QCD, and the order parameter of the two-body model corresponds to the chiral condensate of QCD. Comparing to the corresponding four-body model, we are able to determine which features of the free energy are due to the chaotic nature of the four-body model. The high-temperature phase of both models is entropy dominated, and in both cases, the free energy is determined by the spectral density. The chaotic four-body SYK model has a low-temperature phase whose free energy is almost temperature-independent, signaling an effective gap of the theory even though the actual spectrum does not exhibit a gap. On the other hand, the low-temperature free energy of the two-body SYK model is not flat; in fact, it oscillates to arbitrarily low temperature. This indicates a less desirable feature that the entropy of the two-body model is not always positive in the low-temperature phase, which most likely is a consequence of the non-hermiticity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS2002

Stability of Pseudobrookite-Type Titanium Oxides

Orthorhombic, (Bbmm), (Al, Fe, Cr, Ti)(sub 2) TiO5-(Mg, Fe)Ti2O5 solid solutions (pseudobrookites, s.l.) are found either as an oxidation product of ilmenite and/or spinel or a primary crystallizing phase in igneous and metamorphic rocks on Earth (e.g., basalt flows, crustal and mantle xenoliths, hornfels), and basaltic rocks on the Moon. Moreover, orthorhombic oxides are often part of the crystalline matrix in glass/ceramics with useful applications, and play a major role in the industrial production of TiO2. To fully exploit the potential of these compounds as petrogenetic indicators and/or useful materials we need to quantitatively understand the factors controlling their properties and stability, and thus, to extrapolate beyond the calibrating experiments. For that purpose, we need to combine thermochemistry, phase equilibrium, and in situ P-V-T-cation disorder experimental data that presently either are incomplete or lacking. Perhaps, the most complete data set is that for MgTi2O5 (karrooite) which allows the calibration of models for the Gibbs free energy of the MgTi2O5 as a function of pressure, temperature, and the Mg2+-Ti4+ distribution between the two nonequivalent octahedral sites. Consequently, the effect of cation disorder on MgTi2O5 stability, and the phase relations among MgTi2O5, other titanium oxides, and silicate minerals can be examined. Calculated phase relations in the Mg-Ti-Si-O system and phase equilibrium experiments in Fe-bearing compositions suggest that pseudobrookite-type oxides may be a more common in rocks than previously realized. However, homogeneous and heterogeneous equilibria, and crystallization paths likely affect their stability. For example, isobaric increases in temperature favor disordering and thus entropy-stabilization, in contrast, isothermal increases in pressure have the opposite effect. Although, currently, the potential effect of composition to cation disorder cannot be fully explored, it appears that enrichment in trivalent cations probably enhances entropy-stabilization and thus may increase the stability of (Al, Fe, Cr, Ti)-rich pseudobrookites relative to that of (Mg, Fe)-rich ones. In addition, high-temperature, nearly isothermal, decompression paths of olivine+orthopyroxene+oxide assemblages may favor pseudobrookites (s.l.) over rutile and/or ilmenite, in contrast, cooling at low pressures seems to favor ilmenite and/or rutile. In the case of crustal and mantle xenoliths, the presence or absence of orthorhombic oxides is probably controlled by reactions with olivine, orthopyroxene, ilmenite, and rutile. In oceanic mantle xenoliths such reactions may also involve a TiO2-enriched but not SiO2-enriched melt/fluid, because pseudobrookites (s.l.) would react with the SiO2-enriched melt/fluid to form orthopyroxene and rutile. Parenthetically, experiments and model calculations in the Mg-Ti-Si-O system suggest that low degree partial melting of low-TiO2 bulk compositions may produce Ti-enriched liquids in equilibrium with olivine, orthopyroxen ad=nd MgTi2O5, rutile or ilmenite.

Xirouchakis, Dimistrios↗