Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model fitting”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Analysis of co-spatial UV-optical HST/STIS spectra of planetary nebula NGC 3242

This project sought to consider two important aspects of the planetary nebula NGC 3242 using new long-slit HST/STIS spectra. First, we investigated whether this object is chemically homogeneous by spatially dividing the slit into different regions and calculating the abundances of each region. The major result is that the elements of He, C, O, and Ne are chemically homogeneous within uncertainties across the regions probed, implying that the stellar outflow was well-mixed. Second, we constrained the stellar properties using photoionization models computed by CLOUDY and tested the effects of three different density profiles on these parameters. The three profiles tested were a constant density profile, a Gaussian density profile, and a Gaussian with a power-law density profile. The temperature and luminosity were not affected significantly by the choice of density structure. The values for the stellar temperature and luminosity from our best-fit model are 89.7{sub −4.7}{sup +7.3} kK and log(L/L {sub ⊙}) = 3.36{sub −0.22}{sup +0.28}, respectively. Comparing to evolutionary models on an HR diagram, this corresponds to an initial and final mass of 0.95{sub −0.09}{sup +0.35}M{sub ⊙} and 0.56{sub −0.01}{sup +0.01}M{sub ⊙}, respectively.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

ON THE LAMPPOST MODEL OF ACCRETING BLACK HOLES

We study the lamppost model, in which the X-ray source in accreting black hole (BH) systems is located on the rotation axis close to the horizon. We point out a number of inconsistencies in the widely used lamppost model relxilllp, e.g., neglecting the redshift of the photons emitted by the lamppost that are directly observed. They appear to invalidate those model fitting results for which the source distances from the horizon are within several gravitational radii. Furthermore, if those results were correct, most of the photons produced in the lamppost would be trapped by the BH, and the luminosity generated in the source as measured at infinity would be much larger than that observed. This appears to be in conflict with the observed smooth state transitions between the hard and soft states of X-ray binaries. The required increase of the accretion rate and the associated efficiency reduction also present a problem for active galactic nuclei. Then, those models imply the luminosity measured in the local frame is much higher than that produced in the source and measured at infinity, due to the additional effects of time dilation and redshift, and the electron temperature is significantly higher than that observed. We show that these conditions imply that the fitted sources would be out of the e{sup ±} pair equilibrium. On the other hand, the above issues pose relatively minor problems for sources at large distances from the BH, where relxilllp can still be used.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

ON THE NATURE OF THE TRANSITION DISK AROUND LkCa 15

We present Combined Array for Research in Millimeter-wave Astronomy 1.3 mm continuum observations of the T Tauri star LkCa 15, which resolve the circumstellar dust continuum emission on angular scales between 0.''2 and 3'', corresponding to 28-420 AU at the distance of the star. The observations resolve the inner gap in the dust emission and reveal an asymmetric dust distribution in the outer disk. By comparing the observations with theoretical disk models, we calculate that 90% of the dust emission arises from an azimuthally symmetric ring that contains about 5 Multiplication-Sign 10{sup -4} M{sub Sun} of dust. A low surface-brightness tail that extends to the northwest out to a radius of about 300 AU contains the remaining 10% of the observed continuum emission. The ring is modeled with a rather flat surface density profile between 40 and 120 AU, while the inner cavity is consistent with either a sharp drop of the 1.3 mm dust optical depth at about 42 AU or a smooth inward decrease between 3 and 85 AU. We show that early science Atacama Large Millimeter Array observations will be able to disentangle these two scenarios. Within 40 AU, the observations constrain the amount of dust between 10{sup -6} and 7 Earth masses, where the minimum and maximum limits are set by the near-infrared spectral energy distribution modeling and by the millimeter-wave observations of the dust emission, respectively. In addition, we confirm the discrepancy in the outer disk radius inferred from the dust and gas, which corresponds to 150 AU and 900 AU, respectively. We cannot reconcile this difference by adopting an exponentially tapered surface density profile as suggested for other systems, but we instead suggest that the gas surface density in the outer disk decreases less steeply than that predicted by model fits to the dust continuum emission. The lack of continuum emission at radii larger than 120 AU suggests a drop of at least a factor of five in the dust-to-gas ratio or in the dust opacity. We show that a sharp dust opacity drop of this magnitude is consistent with a radial variation of the grain size distribution as predicted by existing grain growth models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

THE EVOLUTION OF INNER DISK GAS IN TRANSITION DISKS

Investigating the molecular gas in the inner regions of protoplanetary disks (PPDs) provides insight into how the molecular disk environment changes during the transition from primordial to debris disk systems. We conduct a small survey of molecular hydrogen (H{sub 2}) fluorescent emission, using 14 well-studied Classical T Tauri stars at two distinct dust disk evolutionary stages, to explore how the structure of the inner molecular disk changes as the optically thick warm dust dissipates. We simulate the observed Hi-Lyman α-pumped H{sub 2} disk fluorescence by creating a 2D radiative transfer model that describes the radial distributions of H{sub 2} emission in the disk atmosphere and compare these to observations from the Hubble Space Telescope. We find the radial distributions that best describe the observed H{sub 2} FUV emission arising in primordial disk targets (full dust disk) are demonstrably different than those of transition disks (little-to-no warm dust observed). For each best-fit model, we estimate inner and outer disk emission boundaries (r{sub in} and r{sub out}), describing where the bulk of the observed H{sub 2} emission arises in each disk, and we examine correlations between these and several observational disk evolution indicators, such as n{sub 13–31}, r{sub in,} {sub CO}, and the mass accretion rate. We find strong, positive correlations between the H{sub 2} radial distributions and the slope of the dust spectral energy distribution, implying the behavior of the molecular disk atmosphere changes as the inner dust clears in evolving PPDs. Overall, we find that H{sub 2} inner radii are ∼4 times larger in transition systems, while the bulk of the H{sub 2} emission originates inside the dust gap radius for all transitional sources.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

THE IMPORTANCE OF PHYSICAL MODELS FOR DERIVING DUST MASSES AND GRAIN SIZE DISTRIBUTIONS IN SUPERNOVA EJECTA. I. RADIATIVELY HEATED DUST IN THE CRAB NEBULA

Recent far-infrared (IR) observations of supernova remnants (SNRs) have revealed significantly large amounts of newly condensed dust in their ejecta, comparable to the total mass of available refractory elements. The dust masses derived from these observations assume that all the grains of a given species radiate at the same temperature, regardless of the dust heating mechanism or grain radius. In this paper, we derive the dust mass in the ejecta of the Crab Nebula, using a physical model for the heating and radiation from the dust. We adopt a power-law distribution of grain sizes and two different dust compositions (silicates and amorphous carbon), and calculate the heating rate of each dust grain by the radiation from the pulsar wind nebula. We find that the grains attain a continuous range of temperatures, depending on their size and composition. The total mass derived from the best-fit models to the observed IR spectrum is 0.019-0.13 M{sub Sun }, depending on the assumed grain composition. We find that the power-law size distribution of dust grains is characterized by a power-law index of 3.5-4.0 and a maximum grain size larger than 0.1 {mu}m. The grain sizes and composition are consistent with what is expected for dust grains formed in a Type IIP supernova (SN). Our derived dust mass is at least a factor of two less than the mass reported in previous studies of the Crab Nebula that assumed more simplified two-temperature models. These models also require a larger mass of refractory elements to be locked up in dust than was likely available in the ejecta. The results of this study show that a physical model resulting in a realistic distribution of dust temperatures can constrain the dust properties and affect the derived dust masses. Our study may also have important implications for deriving grain properties and mass estimates in other SNRs and for the ultimate question of whether SNe are major sources of dust in the Galactic interstellar medium and in external galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The radio loudness of SDSS quasars from the LOFAR Two-metre Sky Survey: ubiquitous jet activity and constraints on star formation

ABSTRACT We examine the distribution of radio emission from ∼42 000 quasars from the Sloan Digital Sky Survey, as measured in the LOFAR Two-metre Sky Survey (LoTSS). We present a model of the radio luminosity distribution of the quasars that assumes that every quasar displays a superposition of two sources of radio emission: active galactic nuclei (jets) and star formation. Our two-component model provides an excellent match to the observed radio flux density distributions across a wide range of redshifts and quasar optical luminosities; this suggests that the jet-launching mechanism operates in all quasars but with different powering efficiency. The wide distribution of jet powers allows for a smooth transition between the ‘radio-quiet’ and ‘radio-loud’ quasar regimes, without need for any explicit bimodality. The best-fitting model parameters indicate that the star formation rate of quasar host galaxies correlates strongly with quasar luminosity and also increases with redshift at least out to z ∼ 2. For a model where star formation rate scales as $L_{\rm bol}^{\alpha } (1+z)^{\beta }$, we find α = 0.47 ± 0.01 and β = 1.61 ± 0.05, in agreement with far-infrared studies. Quasars contribute ≈0.15 per cent of the cosmic star formation rate density at z = 0.5, rising to 0.4 per cent by z ∼ 2. The typical radio jet power is seen to increase with both increasing optical luminosity and black hole mass independently, but does not vary with redshift, suggesting intrinsic properties govern the production of the radio jets. We discuss the implications of these results for the triggering of quasar activity and the launching of jets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

The epsilon Eridani system resolved by millimeter interferometry

We present observations of ϵ Eridani from the Submillimeter Array (SMA) at 1.3 mm and from the Australia Telescope Compact Array at 7 mm that reach an angular resolution of ∼4″ (13 AU). These first millimeter interferometer observations of ϵ Eridani, which hosts the closest debris disk to the Sun, reveal two distinct emission components: (1) the well-known outer dust belt, which, although patchy, is clearly resolved in the radial direction, and (2) an unresolved source coincident with the position of the star. We use direct model-fitting of the millimeter visibilities to constrain the basic properties of these two components. A simple Gaussian shape for the outer belt fit to the SMA data results in a radial location of 64.4{sub −3.0}{sup +2.4} AU and FWHM of 20.2{sub −8.2}{sup +6.0} AU (fractional width ΔR/R = 0.3). Similar results are obtained taking a power law radial emission profile for the belt, though the power law index cannot be usefully constrained. Within the noise obtained (0.2 mJy beam{sup −1}), these data are consistent with an axisymmetric belt model and show no significant azimuthal structure that might be introduced by unseen planets in the system. These data also limit any stellocentric offset of the belt to <9 AU, which disfavors the presence of giant planets on highly eccentric (>0.1) and wide (10's of AU) orbits. The flux density of the unresolved central component exceeds predictions for the stellar photosphere at these long wavelengths, by a marginally significant amount at 1.3 mm but by a factor of a few at 7 mm (with brightness temperature 13000 ± 1600 K for a source size of the optical stellar radius). We attribute this excess emission to ionized plasma from a stellar corona or chromosphere.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Systematic Matter and Binding-Energy Distributions from a Dispersive Optical Model Analysis

We present the first systematic dispersive-optical-model analysis of 16,18 O, 40,48 Ca, 58,64 Ni, 112,124 Sn, and 208 Pb. In all systems, roughly half the total nuclear binding energy is associated with the most-bound 10% of the total nucleon density. The extracted neutron skins reveal interplay of asymmetry, Coulomb, and shell effects on the skin thickness. Here, our results indicate that optical-model fits on isotopic pairs are effective for constraining nuclear structural quantities otherwise difficult to observe experimentally.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

Fault geometry inversion and slip distribution of the 2010 M w 7.2 El Mayor-Cucapah earthquake from geodetic data

The 4 April 2010 M w 7.2 El Mayor-Cucapah (EMC) earthquake in Baja, California, and Sonora, Mexico, had primarily right-lateral strike-slip motion and a minor normal-slip component. The surface rupture extended about 120 km in a NW-SE direction, west of the Cerro Prieto fault. We use geodetic measurements including near- to far-field GPS, interferometric synthetic aperture radar (InSAR), and subpixel offset measurements of radar and optical images to characterize the fault slip during the EMC event. We use dislocation inversion methods and determine an optimal nine-segment fault geometry, as well as a subfault slip distribution from the geodetic measurements. With systematic perturbation of the fault dip angles, randomly removing one geodetic data constraint, or different data combinations, we are able to explore the robustness of the inferred slip distribution along fault strike and depth. The model fitting residuals imply contributions of early postseismic deformation to the InSAR measurements as well as lateral heterogeneity in the crustal elastic structure between the Peninsular Ranges and the Salton Trough. We also find that with incorporation of near-field geodetic data and finer fault patch size, the shallow slip deficit is reduced in the EMC event by reductions in the level of smoothing. Finally, these results show that the outcomes of coseismic inversions can vary greatly depending on model parameterization and methodology.

58 GEOSCIENCES↗
DOE OSTI2012

CORRELATED SPECTRAL AND TEMPORAL BEHAVIOR OF LATE-TIME AFTERGLOWS OF GAMMA-RAY BURSTS

The cannonball (CB) model of gamma-ray bursts (GRBs) predicts that the asymptotic behavior of the spectral energy density of GRB afterglows is a power law in time and in frequency, and the difference between the temporal and spectral power-law indices, {alpha}{sub X} - {beta}{sub X}, is restricted to the values 0, 1/2, and 1. Here we report the distributions of the values {alpha}{sub X} and {beta}{sub X}, and their difference for a sample of 315 Swift GRBs. This sample includes all Swift GRBs that were detected before 2012 August 1, whose X-ray afterglow extended well beyond 1 day and the estimated error in {alpha}{sub X} - {beta}{sub X} was {<=}0.25. The values of {alpha}{sub X} were extracted from the CB-model fits to the entire light curves of their X-ray afterglow while the spectral index was extracted by the Swift team from the time-integrated X-ray afterglow of these GRBs. We found that the distribution of the difference {alpha}{sub X} - {beta}{sub X} for these 315 Swift GRBs has three narrow peaks around 0, 1/2, and 1 whose widths are consistent with being due to the measurement errors, in agreement with the CB-model prediction.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Evaporation and scattering of momentum- and velocity-dependent dark matter in the Sun

Dark matter with momentum- or velocity-dependent interactions with nuclei has shown significant promise for explaining the so-called Solar Abundance Problem, a longstanding discrepancy between solar spectroscopy and helioseismology. The best-fit models are all rather light, typically with masses in the range of 3–5 GeV. This is exactly the mass range where dark matter evaporation from the Sun can be important, but to date no detailed calculation of the evaporation of such models has been performed. Here we carry out this calculation, for the first time including arbitrary velocity- and momentum-dependent interactions, thermal effects, and a completely general treatment valid from the optically thin limit all the way through to the optically thick regime. We find that depending on the dark matter mass, interaction strength and type, the mass below which evaporation is relevant can vary from 1 to 4 GeV. This has the effect of weakening some of the better-fitting solutions to the Solar Abundance Problem, but also improving a number of others. As a by-product, we also provide an improved derivation of the capture rate that takes into account thermal and optical depth effects, allowing the standard result to be smoothly matched to the well-known saturation limit.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

Laboratory and Lysimeter Experimentation and Transport Modeling of Neptunium and Strontium in Savannah River Site Sediments

The Savannah River Site (SRS) conducts performance assessment (PA) calculations to determine the appropriate amount of low-level radiological waste that can be safely disposed on site. Parameters are included in these calculations that account for the interaction between the immobile solid phase and the mobile aqueous phase. These parameters are either the distribution coefficient (K$_{d}$ value) or the apparent solubility value (K$_{sp}$). These parameters are readily found in the literature and are used throughout the DOE complex. One shortcoming of K$_{d}$ values is that they are only applicable to a given set of solid and aqueous phase conditions. Therefore, a given radionuclide may have several K$_{d}$ values as it moves between formations and comes into contact with different solids and different aqueous phases. It is expected that the K$_{d}$ construct will be appropriate to use for a majority of the PA and for a majority of the radionuclides. However, semi-mechanistic models would be more representative in isolated cases where the chemistry is especially transitory or the radionuclide chemistry is especially complex, bringing to bear multiple species of varying sorption tendencies to the sediment. Semi-mechanistic models explicitly accommodate the dependency of K$_{d}$ values, or other sorption parameters, on contaminant concentration, competing ion concentrations, pH-dependent surface charge on the adsorbent, and solute species distribution. Incorporating semi-mechanistic concepts into geochemical models is desirable to make the models more robust and technically defensible. Furthermore, these alternative models could be used to augment or validate a K$_{d}$-based DOE Order 435.1 Performance Assessment. The objectives of this study were to: 1) develop a quantitative thermodynamically-based model for neptunium sorption to SRS sediments, and 2) determine a sorption constant from an SRS 11-year lysimeter study. The modeling studies were conducted with existing data sets. The first data set used laboratory generated Np sorption data as a function of concentration (three orders of magnitude) and as a function of pH (four orders of magnitude of proton concentration). In this modeling exercise, a very simple solution was identified by assuming that all sorption occurred only to the iron oxides in the sediment and that all the added NpO$_{4}$$^{-}$ remained in the oxidized state and was not reduced to the Np(IV) state (as occurs rapidly with Pu(V)). With rather limited input data, very good agreement between experimental and modeling results was observed. This modeling approach would be easy to add to the PA with little additional data requirements. This model would be useful in a system where pH is expected to change greatly, such as directly beneath a grout or concrete structure. The second model discussed in the report was to derive strontium K$_{d}$ values from data collected in an 11-year-old field transport study. In this controlled lysimeter study, a sensitivity analysis was conducted of hydrological and chemical processes that influence contaminant transport, including diffusion coefficients, seepage velocity, and K$_{d}$ value. The best overall K$_{d}$ derived from the model fit to the data was 32 L kg$^{-1}$, which was the same value that was previously measured in traditional laboratory batch sorption studies. This was an unexpected result given the differences in experimental conditions between the batch test and the lysimeter flow through test, in particular the differences between strontium adsorption and desorption processes occurring in the latter test and not in the former. There were some trends in the lysimeter strontium data that were not predicted by the K$_{d}$ model, which suggest that other geochemical processes are likely also controlling strontium transport. Strontium release and cation exchange are being evaluated. These results suggest that future modeling efforts (e.g., PAs) could be improved by employing a more robust semi-empirical modeling approach to transient or complex conditions.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Bay Area Regional Energy: Network Integrated Commercial Retrofits (BRICR) Project. Final Report

The BRICR project applied large-scale building energy modeling concepts with the aim of reducing the cost of energy efficiency targeting, design, and project development, and measurement of energy savings for energy efficiency programs implemented by local governments that serve small and medium commercial buildings (SMB). The project leveraged the services and resources of existing local government energy programs serving disadvantaged and hard-to-reach SMB customers. In contrast to programs run by utilities, local government programs generally do not have direct access to energy billing records for an entire class of customers in a geographic area, which prior research demonstrated useful for large-scale building energy model baseline development and calibration. , However, local governments are rich in public records that offer important clues about physical attributes and uses that, along with behavior, determine energy use. Relying only on public records, BRICR demonstrated development of credible baseline energy models for 3,792 office, retail, and hotel buildings. Publicly disclosed annual energy use data from a local energy benchmarking program and anonymized data from the Building Performance Database, the nation’s largest dataset about energy-related characteristics of buildings, were utilized to validate and calibrate energy models via an innovative method comparing distributions of energy intensity by fuel type for portfolios of buildings of similar size, vintage, and use. Portfolio calibration does not provide certainty that an energy model fits an individual building; the method is useful when billing data is not accessible – a common situation for researchers, energy service providers and ESCOs, local governments, and any party other than a utility. A software component was developed, the BRICR gem, which automates simulation when relevant data is added or edited by the user to a file saved in the standardized BuildingSync XML schema for energy audit data. The component was demonstrated as a simplified means to generate a mass of energy models corresponding to public records containing basic attributes such as building scale, location, use, year built, and aspect ratio in combination with building energy code prototype data corresponding to use and vintage. The component was also demonstrated as a simplified means to automate energy simulation when attributes are revised; the intention was to enable iterative improvement of the baseline model and energy savings estimates for common energy conservation measures as users revise relevant attributes based on their observations. In the context of institutional change and uncertainty for the participating local government energy programs, 13 whole building retrofits were completed. Impacts were measured by applying the CalTRACK2.0 methods to standardize measurement of normalized metered energy consumption. The GRIDMeter methods of stratified sampling and individual load shape analysis were applied to adjust for impacts of the effect of COVID-19 on retrofitted buildings in the context of all local buildings of similar size and use. Excluding impacts of the pandemic, retrofitted buildings demonstrated between 1.6% and 25.1% reduction in energy use. The project contributed use cases and feedback that helped inform evolution of the software tools and data formats that were combined for the first time in the BRICR project.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

Evaluation of PV Module Packaging Strategies of Monofacial and Bifacial PERC Using Degradation Pathway Network Modeling

As the PV industry is rapidly expanding, it is important to thoroughly investigate the long-term impact of packaging strategies on the performance of PV modules. In this study, the variants in sets differ on the basis of manufacturer (A/B), encapsulant (EVA/POE), rear encapsulant (UV-cutoff/opaque/transparent), module architecture (GB/DG) and cell type (monofacial/bifacial). The minimodules were exposed for 2520 hours in modified damp heat, with or without full spectrum light. Every 504 hours, stepwise electrical characterization techniques were employed to track changes in minimodules. Degradation pathway modeling using network structural equation modeling was employed to study pairwise relationships between variables and service lifetime prediction in minimodules. Through this study, differences in quality control are identified in minimodules made by different manufacturers. Minimodules with UV-cutoff rear encapsulant show relatively better stability, whereas the ones with opaque rear encapsulant show greater power loss. In addition, GB having UV-cutoff rear encapsulation and GB with POE having opaque rear encapsulation were identified to be stable as they lack a best model fit. Here, the primary power loss mechanism in degrading variants is interconnect corrosion.

14 SOLAR ENERGY↗
DOE OSTI2017

Scattering Properties of Large Irregular Cosmic Dust Particles at Visible Wavelengths

The effect of internal inhomogeneities and surface roughness on the scattering behavior of large cosmic dust particles is studied by comparing model simulations with laboratory measurements. The present work shows the results of an attempt to model a dust sample measured in the laboratory with simulations performed by a ray-optics model code. We consider this dust sample as a good analogue for interplanetary and interstellar dust as it shares its refractive index with known materials in these media. Several sensitivity tests have been performed for both structural cases (internal inclusions and surface roughness). Three different samples have been selected to mimic inclusion/coating inhomogeneities: two measured scattering matrices of hematite and white clay, and a simulated matrix for water ice. These three matrices are selected to cover a wide range of imaginary refractive indices. The selection of these materials also seeks to study astrophysical environments of interest such as Mars, where hematite and clays have been detected, and comets. Based on the results of the sensitivity tests shown in this work, we perform calculations for a size distribution of a silicate-type host particle model with inclusions and surface roughness to reproduce the experimental measurements of a dust sample. The model fits the measurements quite well, proving that surface roughness and internal structure play a role in the scattering pattern of irregular cosmic dust particles.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Kinetics of D/H isotope fractionation between molecular hydrogen and water

At equilibrium, the D/H isotope fractionation factor between H 2 and H 2 O (α H 2 O-H 2 (eq) ) is a sensitive indicator of temperature, and has been used as a geothermometer for natural springs and gas discharges. However, δD H 2 measured in spring waters may underestimate subsurface temperatures of origin due to partial isotopic re-equilibration during ascent and cooling. In this work, we present new experimental data on the kinetics of D–H exchange for H 2 dissolved in liquid water at temperatures below 100 °C. Comparing these results with published exchange rates obtained from gas phase experiments (100–400 °C), we derive a consistent activation energy of 52 kJ/mol, and the following rate expressions; ln k =9.186 - 6298/ T and k 1 = 9764.61[H 2 O] e-6298/ T where T is absolute temperature (K), k is the universal rate constant ([L/mol]/hr), and k 1 is a pseudo-first-order constant (hr -1 ) applicable to water-dominated terrestrial systems by constraining [H 2 O] as the density of H 2 O (in mol/L) at the P-T of interest. The density-dependent rate constant accounts for the kinetic disparity of D–H exchange with H 2 when dissolved in liquid H 2 O relative to a gas/steam phase, exemplifed by 1/ k 1 at 100 °C of ~2 days in liquid, versus ~7 yrs in saturated steam. This difference may explain the high variability of δDH 2 observed in fumarolic gases. Fluids convecting in the crust frequently reach T > 225 °C, where isotopic equilibrium is rapidly attained (<1 hr). We compare fractionation factors measured in natural fluids (α OBS ) with values expected for equilibrium at the T of acquisition. Where these values differ, we use kinetic models to estimate cooling rates during upward advection that account for the observed disequilibrium. Models fit to fluids from Yellowstone Park and the Lost City (deep-sea) vent field, both recovered at ~90 °C, require respective transit times of ~7 hrs and ~11 days between higher temperature reaction zones and the surface. Using estimates of subsurface depths of origin, however, suggests similar mean fluid flow rates (10 s of meters/hr). Additional complications must be considered when interpreting the δDH 2 of lower-temperature effluent. When applied to data from deep-sea hydrothermal systems, our kinetic models indicate microbial catalysis accelerates D–H exchange once fluids cool below ~60 °C. The H 2 measured in both continental alkaline springs and fracture fluids from Precambrian shield rock is likely produced at T < 100 °C, through processes such as serpentinization. In these settings, δDH 2 values appear closer to equilibrium with H 2 O than those from geothermal systems. Considering kinetic isotope effects may yield H 2 that is out of equilibrium when generated at lower temperatures, we calculate maximum (isothermal) times to apparent isotopic equilibrium of 1.3 yrs at 50 °C, 9 yrs at 25 °C, and 35 yrs at 10 °C. A similar calculation applied to Antarctic brines (-13 °C), where measured δDH 2 is far from equilibrium, yields ~350 yrs. This time is shorter than the fluids have been isolated (2.8 ka), suggesting kinetic isotope effects associated with H 2 destruction or loss via diffusion may also be possible.

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

Thermal-model-based characterization of heavy-ion-collision systems at chemical freeze-out

We investigate the chemical freeze-out in heavy-ion collisions (HICs) and the impact of the hadronic spectrum on thermal model analyses [1, 2]. Detailed knowledge of the hadronic spectrum is still an open question, which has phenomenological consequences on the study of HICs. By varying the number of resonances included in Hadron Resonance Gas (HRG) Model calculations, we can shed light on which particles may be produced. Furthermore, we study the influence of the number of states on the so-called two flavor freezeout scenario, in which strange and light particles can freeze-out separately. We consider results for the chemical freeze-out parameters obtained from thermal model fits and from calculating net-particle fluctuations. We will show the effect of using one global temperature to fit all particles and alternatively, allowing particles with and without strange quarks to freeze-out separately.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

The Three-dimensional Expansion of the Ejecta from Tycho's Supernova Remnant

We present the first 3D measurements of the velocity of various ejecta knots in Tycho’s supernova remnant, known to result from a Type Ia explosion. Chandra X-ray observations over a 12 yr baseline from 2003 to 2015 allow us to measure the proper motion of nearly 60 “tufts” of Si-rich ejecta, giving us the velocity in the plane of the sky. For the line-of-sight velocity, we use two different methods: a nonequilibrium ionization model fit to the strong Si and S lines in the 1.2–2.8 keV regime, and a fit consisting of a series of Gaussian lines. These methods give consistent results, allowing us to determine the redshift or blueshift of each of the knots. Assuming a distance of 3.5 kpc, we find total velocities that range from 2400 to 6600 km s{sup −1}, with a mean of 4430 km s{sup −1}. We find several regions where the ejecta knots have overtaken the forward shock. These regions have proper motions in excess of 6000 km s{sup −1}. Some SN Ia explosion models predict a velocity asymmetry in the ejecta. We find no such velocity asymmetries in Tycho, and we discuss our findings in light of various explosion models, favoring those delayed-detonation models with relatively vigorous and symmetrical deflagrations. Finally, we compare measurements with models of the remnant’s evolution that include both smooth and clumpy ejecta profiles, finding that both ejecta profiles can be accommodated by the observations.

79 ASTRONOMY AND ASTROPHYSICS↗