Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 1,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Learning electric vehicle driver range anxiety with an initial state of charge-oriented gradient boosting approach

This manuscript focuses on the modeling of electric vehicle (EV) driver’s range anxiety, a fear that a vehicle does not have sufficient range, or state of charge (SOC) of the battery pack, to reach its destination and would strand its occupants. Despite numerous research studies on the modeling of charging behaviors, modeling efforts to understand at what battery percentages do EV drivers charge their vehicles, and what are the associated contributing factors, are rather limited. To this end, an ensemble learning model based on gradient boosting is developed. The model sequentially fits new predictors to new residuals of the previous prediction and, then, minimizes the loss when adding the latest prediction. A total of 18 features are defined and extracted from the multisource data, which cover information on driver, vehicles, stations, traffic conditions, as well as spatial-temporal context information of the charging events. The analyzed dataset includes 4.5-year’s charging event log data from 3,096 users and 468 public charging stations in Kansas City Missouri, and the macroscopic travel demand model maintained by the metropolitan planning organization. Here, the result shows the proposed model achieved a satisfactory result with a R square value of 0.54 and root mean square error of 0.14, both better than multiple linear regression model and random forest model. To reduce range anxiety, it is suggested that the priorities of deploying new charging facilities should be given to the areas with higher daily traffic prediction, with more conservative EV users or that are further from residential areas.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2013

THE NATURE OF THE SECOND PARAMETER IN THE IRX-{beta} RELATION FOR LOCAL GALAXIES

We present an analysis of 98 galaxies of low-dust content, selected from the Spitzer Local Volume Legacy survey, aimed at examining the relation between the ultraviolet (UV) color and dust attenuation in normal star-forming galaxies. The IRX-{beta} diagram relates the total dust attenuation in a galaxy, traced by the far-IR (FIR) to UV ratio, to the observed UV color, indicated by {beta}. Previous research has indicated that while starburst galaxies exhibit a relatively tight IRX-{beta} relation, normal star-forming galaxies do not, and have a much larger spread in the total-IR to far-UV (FUV) luminosity for a fixed UV color. We examine the role that the age of the stellar population plays as the ''second parameter'' responsible for the observed deviation and spread of star-forming galaxies from the starburst relation. We model the FUV to FIR spectral energy distribution of each galaxy according to two broad bins of star formation history (SFH): constant and instantaneous burst. We find clear trends between stellar population mean age estimators (extinction-corrected FUV/NIR, U - B, and EW(H{alpha})) and the UV color {beta}; the trends are mostly driven by the galaxies best-described by instantaneous burst populations. We also find a significant correlation between {beta} and the mean age directly determined from the best-fit instantaneous models. As already indicated by other authors, the UV attenuation in star-forming galaxies may not be recovered with the UV color alone and is highly influenced by the stellar population's mean age and SFH. Overall, the scatter in the IRX-{beta} diagram is better correlated with {beta} than with the perpendicular distance, d{sub p}.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

A dendritic strontium river isoscape for fisheries applications in the Sacramento River basin, California, USA

Objective Understanding the origins and movements of fish is fundamental to effective conservation and fisheries management. Strontium isotope ratios ( 87 Sr/ 86 Sr) in otoliths provide a powerful tracer of natal origin and migratory pathways. However, existing 87 Sr/ 86 Sr isoscapes for the Sacramento River basin, an ecosystem that supports ecologically and economically important salmon populations, rely on discrete classification approaches that overlook unsampled habitats and do not incorporate spatial uncertainty. Our objective was to develop a continuous, network-explicit 87 Sr/ 86 Sr isoscape with quantified uncertainty to fill in data gaps and enable probabilistic assignments of fish origin and movement. Methods We used river water 87 Sr/ 86 Sr data from 106 sites (1997–2021) to develop spatial stream network models that use dendritic connectivity and watershed characteristics (lithology, bedrock age, and land cover) to predict river water 87 Sr/ 86 Sr throughout the basin. Models were fitted using maximum and restricted likelihood and were evaluated via Akaike’s information criterion and leave-one-out cross validation. We produced both historical (pre-dam) and present-day (below-dam) isoscapes, delineated uncertainty-informed isotopic ranges using k -means clustering, and applied a proof-of-concept Bayesian assignment to estimate natal origins and early rearing habitats for two endangered winter-run Chinook Salmon Oncorhynchus tshawytscha. Results Cross validation indicated strong performance of the 87 Sr/ 86 Sr model (leave-one-out cross validation: R 2 = 0.91; root mean square error = 0.0005). Uncertainty-informed clustering identified 19 isotopic “suites” (reaches with indistinguishable 87 Sr/ 86 Sr values) in present-day anadromous habitats and 25 suites in the historical network. Example natal and early rearing assignments included predictions that challenged expectations for juvenile salmon migration based on predicted river 87 Sr/ 86 Sr compositions. Conclusions This study developed a continuous, network-explicit 87 Sr/ 86 Sr isoscape that integrates existing river data to predict 87 Sr/ 86 Sr in unsampled reaches and the likely achievable range and resolution of otolith-based origin and life history inference. The resulting river isoscape provides a valuable tool to predict salmon movements and identify habitats supporting their survival and growth that otherwise might remain undetected. Coupling these predictions with complementary approaches that ground-truth juvenile presence (e.g., targeted fish surveys) represents an important step toward science-informed restoration and management of critical habitats throughout the Sacramento River basin.

Environmental sciences↗
DOE OSTI2023

A data-driven method for modelling dissipation rates in stratified turbulence

We present a deep probabilistic convolutional neural network (PCNN) model for predicting local values of small-scale mixing properties in stratified turbulent flows, namely the dissipation rates of turbulent kinetic energy and density variance, $\varepsilon$ and $\chi$ . Inputs to the PCNN are vertical columns of velocity and density gradients, motivated by data typically available from microstructure profilers in the ocean. The architecture is designed to enable the model to capture several characteristic features of stratified turbulence, in particular the dependence of small-scale isotropy on the buoyancy Reynolds number $Re_b:=\varepsilon /(\nu N^2)$ , where $\nu$ is the kinematic viscosity and $N$ is the background buoyancy frequency, the correlation between suitably locally averaged density gradients and turbulence intensity and the importance of capturing the tails of the probability distribution functions of values of dissipation. Empirically modified versions of commonly used isotropic models for $\varepsilon$ and $\chi$ that depend only on vertical derivatives of density and velocity are proposed based on the asymptotic regimes $Re_b\ll 1$ and $Re_b\gg 1$ , and serve as an instructive benchmark for comparison with the data-driven approach. When trained and tested on a simulation of stratified decaying turbulence which accesses a range of turbulent regimes (associated with differing values of $Re_b$ ), the PCNN outperforms assumptions of isotropy significantly as $Re_b$ decreases, and additionally demonstrates improvements over the fitted empirical models. A differential sensitivity analysis of the PCNN facilitates a comparison with the theoretical models and provides a physical interpretation of the features enabling it to make improved predictions.

42 ENGINEERING↗
DOE OSTI2017

Detection of Baryon Acoustic Oscillation features in the large-scale 3-point correlation function of SDSS BOSS DR12 CMASS galaxies

We present the large-scale 3-point correlation function (3PCF) of the SDSS DR12 CMASS sample of 777,202 Luminous Red Galaxies, the largest-ever sample used for a 3PCF or bispectrum measurement. We make the first high-significance (4.5σ) detection of Baryon Acoustic Oscillations (BAO) in the 3PCF. Using these acoustic features in the 3PCF as a standard ruler, we measure the distance to z=0.57 to 1.7% precision (statistical plus systematic). We find D V = 2024 ± 29Mpc (stat) ± 20Mpc(sys) for our fiducial cosmology (consistent with Planck 2015) and bias model. This measurement extends the use of the BAO technique from the 2-point correlation function (2PCF) and power spectrum to the 3PCF and opens an avenue for deriving additional cosmological distance information from future large-scale structure redshift surveys such as DESI. Our measured distance scale from the 3PCF is fairly independent from that derived from the pre-reconstruction 2PCF and is equivalent to increasing the length of BOSS by roughly 10%; reconstruction appears to lower the independence of the distance measurements. In conclusion, fitting a model including tidal tensor bias yields a moderate significance (2.6σ) detection of this bias with a value in agreement with the prediction from local Lagrangian biasing.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Exploring cryo-electron microscopy with molecular dynamics

Single particle analysis cryo-electron microscopy (EM) and molecular dynamics (MD) have been complimentary methods since cryo-EM was first applied to the field of structural biology. The relationship started by biasing structural models to fit low-resolution cryo-EM maps of large macromolecular complexes not amenable to crystallization. The connection between cryo-EM and MD evolved as cryo-EM maps improved in resolution, allowing advanced sampling algorithms to simultaneously refine backbone and sidechains. Moving beyond a single static snapshot, modern inferencing approaches integrate cryo-EM and MD to generate structural ensembles from cryo-EM map data or directly from the particle images themselves. We summarize the recent history of MD innovations in the area of cryo-EM modeling. The merits for the myriad of MD based cryo-EM modeling methods are discussed, as well as, the discoveries that were made possible by the integration of molecular modeling with cryo-EM. Lastly, current challenges and potential opportunities are reviewed.

Biochemistry & Molecular Biology↗
DOE OSTI2026

Linking Community‐Climate Disequilibrium to Ecosystem Function

Turnover in species composition often lags behind the pace of climate change, resulting in mismatches between climate and communities. However, the impact of these community‐climate disequilibria on ecosystem functions is rarely considered, and current methods for measuring disequilibria assume that species ranges were, until recently, in equilibrium with climate. Here, in this work, we develop a simple theoretical model to address both of these problems by linking community‐climate disequilibrium with ecosystem functioning. We show how disequilibrium can impair functioning in the near‐term even when climate change is expected to enhance functioning in the long‐term. Responses are most likely to change over time in communities where turnover is slow, the impact of disequilibrium counteracts the direct effects of climate on ecosystem function, and pre‐existing disequilibrium is large. These findings emphasise the importance of precise and unbiased estimates of community‐climate disequilibria for improving ecological forecasts. By fitting our model to time series of both climate and ecosystem function from a metacommunity simulation, we show the potential for community‐climate disequilibrium to be inferred without direct knowledge about species' distributions or climatic tolerances. We end by outlining a research agenda to apply dynamic disequilibrium concepts and test novel hypotheses across diverse ecosystems.

climate change↗
DOE OSTI2016

Probing primordial features with future galaxy surveys

We study the capability of future measurements of the galaxy clustering power spectrum to probe departures from a power-law spectrum for primordial fluctuations. On considering the information from the galaxy clustering power spectrum up to quasi-linear scales, i.e. k < 0.1 h Mpc{sup −1}, we present forecasts for DESI, Euclid and SPHEREx in combination with CMB measurements. As examples of departures in the primordial power spectrum from a simple power-law, we consider four Planck 2015 best-fits motivated by inflationary models with different breaking of the slow-roll approximation. At present, these four representative models provide an improved fit to CMB temperature anisotropies, although not at statistical significant level. As for other extensions in the matter content of the simplest ΛCDM model, the complementarity of the information in the resulting matter power spectrum expected from these galaxy surveys and in the primordial power spectrum from CMB anisotropies can be effective in constraining cosmological models. We find that the three galaxy surveys can add significant information to CMB to better constrain the extra parameters of the four models considered.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Volume-translated cubic EoS and PC-SAFT density models and a free volume-based viscosity model for hydrocarbons at extreme temperature and pressure conditions

This research focuses on providing the petroleum reservoir engineering community with robust models of hydrocarbon density and viscosity at the extreme temperature and pressure conditions (up to 533 K and 276 MPa, respectively) characteristic of ultra-deep reservoirs, such as those associated with the deepwater wells in the Gulf of Mexico. Our strategy is to base the volume-translated (VT) Peng–Robinson (PR) and Soave–Redlich–Kwong (SRK) cubic equations of state (EoSs) and perturbed-chain, statistical associating fluid theory (PC-SAFT) on an extensive data base of high temperature (278–533 K), high pressure (6.9–276 MPa) density rather than fitting the models to low pressure saturated liquid density data. This high-temperature, high-pressure (HTHP) data base consists of literature data for hydrocarbons ranging from methane to C{sub 40}. The three new models developed in this work, HTHP VT-PR EoS, HTHP VT-SRK EoS, and hybrid PC-SAFT, yield mean absolute percent deviation values (MAPD) for HTHP hydrocarbon density of ~2.0%, ~1.5%, and <1.0%, respectively. An effort was also made to provide accurate hydrocarbon viscosity models based on literature data. Viscosity values are estimated with the frictional theory (f-theory) and free volume (FV) theory of viscosity. The best results were obtained when the PC-SAFT equation was used to obtain both the attractive and repulsive pressure inputs to f-theory, and the density input to FV theory. Both viscosity models provide accurate results at pressures to 100 MPa but experimental and model results can deviate by more than 25% at pressures above 200 MPa.

02 PETROLEUM↗
DOE OSTI2021

Mechanistic basis of post-treatment control of SIV after anti-α4β7 antibody therapy

Treating macaques with an anti-α4β7 antibody under the umbrella of combination antiretroviral therapy (cART) during early SIV infection can lead to viral remission, with viral loads maintained at < 50 SIV RNA copies/ml after removal of all treatment in a subset of animals. Depletion of CD8+ lymphocytes in controllers resulted in transient recrudescence of viremia, suggesting that the combination of cART and anti-α4β7 antibody treatment led to a state where ongoing immune responses kept the virus undetectable in the absence of treatment. A previous mathematical model of HIV infection and cART incorporates immune effector cell responses and exhibits the property of two different viral load set-points. While the lower set-point could correspond to the attainment of long-term viral remission, attaining the higher set-point may be the result of viral rebound. Here we expand that model to include possible mechanisms of action of an anti-α4β7 antibody operating in these treated animals. We show that the model can fit the longitudinal viral load data from both IgG control and anti-α4β7 antibody treated macaques, suggesting explanations for the viral control associated with cART and an anti-α4β7 antibody treatment. This effective perturbation to the virus-host interaction can also explain observations in other nonhuman primate experiments in which cART and immunotherapy have led to post-treatment control or resetting of the viral load set-point. Interestingly, because the viral kinetics in the various treated animals differed—some animals exhibited large fluctuations in viral load after cART cessation—the model suggests that anti-α4β7 treatment could act by different primary mechanisms in different animals and still lead to post-treatment viral control. This outcome is nonetheless in accordance with a model with two stable viral load set-points, in which therapy can perturb the system from one set-point to a lower one through different biological mechanisms.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

A Multi-frequency analysis of dark matter annihilation interpretations of recent anti-particle and γ-ray excesses in cosmic structures

The Fermi-LAT observation of a γ-ray excess from the galactic-centre, as well as the PAMELA, AMS, and AMS-2 anti-particle excesses, and the recent indications of a Fermi-LAT γ-ray excess in the Reticulum II dwarf galaxy have all been variously put forward as possible indirect signatures of supersymmetric neutralino dark matter. These are of particular interest as the neutralino annihilation models which fit these observations must have observable consequences across the frequency spectrum, from radio to γ-ray emission. Moreover, since dark matter is expected to be a major constituent of cosmic structure, these multi-frequency consequences should be common to such structures across the mass spectrum, from dwarf galaxies to galaxy clusters. Thus, in this work we make predictions for the multi-frequency spectra of three well-known sources dominated by dark matter on cluster, galaxy and dwarf galaxy scales, e.g. the Coma cluster, the galaxy M81, and the Draco dwarf galaxy, using models favoured by dark matter interpretations of the aforementioned observations. We pay special attention to the consequences for these models when their cross-sections are renormalised to reproduce the recent γ-ray excess observed in the Reticulum II dwarf galaxy, as well as using cross-sections from the Fermi-LAT dwarf galaxy limits, which throw a dark matter interpretation of this excess into doubt. We find that the multi-frequency data of Coma and Draco are in conflict with the dark matter interpretation of the AMS, PAMELA and Fermi positron excess. Additionally, models derived from Fermi-LAT galactic centre observations, and AMS-2 re-analysis, present similar but less extensive conflicts. Using the sensitivity projections for the Square Kilometre Array, the Cherenkov Telescope Array, as well as the ASTROGAM and ASTRO-H satellites, we determine the detection prospects for a subset of neutralino models that remain consistent with Planck cosmological constraints. Although the SKA has the greatest sensitivity to dark matter models, we demonstrate that ASTRO-H is well positioned to probe the inverse-Compton scattering emissions from neutralino annihilation and identify characteristics of the spectra which contain information about the neutralino mass and annihilation channel. This means that, given environments with favourable X-ray backgrounds, multi-frequency observation with the next generation of experiments will allow for unprecedented sensitivity to the neutralino parameter space as well as offsetting the individual weaknesses of each observation mode. Finally we show that all of the studied models can be better tested with the SKA phase 1.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

bmdrc: Python package for quantifying phenotypes from chemical exposures with benchmark dose modeling

Though chemical exposures are known to potentially have negative impacts on health, including contributing to chronic diseases such as cancer, the quantitative contribution of risk is not fully understood for every chemical. A commonly used approach to quantify levels of risk is to measure the proportion of organisms (such as a total number of zebrafish on a plate or mice in a cage) with abnormal behavioral responses or morphology at increasing concentrations of chemical exposure. A particular challenge with processing the proportional data from these assays is the appropriate estimation of chemical concentration levels that result in malformations or acute toxicity, as these values typically vary between experimental measurements. The recommended approach by the Environmental Protection Agency (EPA) is to fit benchmark dose curves with specific filters and model fitting steps, which are crucial to properly processing the proportional data. Several tools exist for the fitting of benchmark dose response curves, but none are standalone Python libraries built to process both morphological and behavioral data as proportions with all the EPA recommended filters, filter parameters, models, and model parameters. Thus, here we present the benchmark dose response curve (bmdrc) Python library, which was built to closely follow these EPA guidelines with helpful visualizations of filters and fitted model curves, and reports for reproducibility purposes. bmdrc is open-source and has demonstrated utility as a support package to an existing web portal for information on chemicals (https://srp.pnnl.gov). Our package will support any toxicology analysis where the response is a proportional value at increasing levels of a concentration of a chemical or chemical mixture.

Superfund↗
DOE OSTI2020

Probing Transit Timing Variation and Its Possible Origin with 12 New Transits of TrES-3b

We present 12 new transit light curves of the hot-Jupiter TrES-3b observed during 2012−2018 to probe the transit timing variation (TTV). By combining the midtransit times determined from these 12 transit data with those reestimated through uniform procedure from 71 transit data available in the literature, we derive new linear ephemeris and obtain the timing residuals that suggest the possibility of TTV in the TrES-3 system. However, the frequency analysis shows that the possible TTV is unlikely to be periodic, indicating the absence of an additional body in this system. To explore the other possible origins of TTV, the orbital decay and apsidal precession ephemeris models are fitted to the transit time data. We find the decay rate of TrES-3b to be P-dot {sub q}=−4.1±3.1 ms yr{sup −1}, and the corresponding estimated modified stellar tidal quality factor of Q{sub ∗}{sup ′}∼1.11×10{sup 5} is consistent with the theoretically predicted values for the stars hosting the hot-Jupiters. The shift in the transit arrival time of TrES-3b after 11 years is expected to be T {sub shift} ∼ 69.55 s, which is consistent with the rms of the timing residuals. Besides, we find that the apsidal precession ephemeris model is statistically less probable than the other considered ephemeris models. It is also discussed that despite the fact that the linear ephemeris model appears to be the most plausible model to represent the transit time data, the possibility of the orbital decay cannot be completely ruled out in the TrES-3 system. To confirm this, further high-precision and high-cadence follow-up observation of transits of TrES-3b would be important.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Mathematical modeling and parameter estimation of MFI membranes for para/ortho‐xylene separation

Abstract A mathematical model of a membrane with a thin, oriented, and selective MFI layer, which also includes contributions from defects, pore blockages, support layer, and external mass transfer, was formulated based on the Maxwell‐Stefan equations. It was validated using reported (Kim et al., Angew Chemie Int Ed, 2018, 57:480–485; Jeon et al., Nature, 2017, 543:690–694) para/ortho‐xylene separation data from five MFI membranes. The diffusivities of the xylenes were considered to be the same for all membranes, while the thickness and the density of defects and pore blockage were treated as unique properties of each membrane, and their contributions were estimated by fitting the model to the corresponding separation data. The effects of these properties and the role of permeate pressure on the separation performance were subsequently analyzed. The proposed modeling, parameter estimation, and analysis framework allow one to quantitatively interpret the variation of separation performance, to understand separation bottlenecks, and to provide guidance for designing membranes with desired performance.

Liu, Jingjun↗
DOE OSTI2012

CHANDRA AND HST IMAGING OF THE QUASARS PKS B0106+013 AND 3C 345: INVERSE COMPTON X-RAYS AND MAGNETIZED JETS

We present results from deep ({approx}70 ks) Chandra/ACIS observations and Hubble Space Telescope (HST) Advanced Camera for Surveys F475W observations of two highly optically polarized quasars belonging to the MOJAVE blazar sample, viz., PKS B0106+013 and 1641+399 (3C 345). These observations reveal X-ray and optical emissions from the jets in both sources. X-ray emission is detected from the entire length of the 0106+013 radio jet, which shows clear bends or wiggles-the X-ray emission is brightest at the first prominent kiloparsec jet bend. A picture of a helical kiloparsec jet with the first kiloparsec-scale bend representing a jet segment moving close(r) to our line of sight, and getting Doppler boosted at both radio and X-ray frequencies, is consistent with these observations. The X-ray emission from the jet end, however, peaks at about 0.''4 ({approx}3.4 kpc) upstream of the radio hot spot. Optical emission is detected both at the X-ray jet termination peak and at the radio hot spot. The X-ray jet termination peak is found upstream of the radio hot spot by around 0.''2 ({approx}1.3 kpc) in the short projected jet of 3C 345. HST optical emission is seen in an arc-like structure coincident with the bright radio hot spot, which we propose is a sharp (apparent) jet bend instead of a terminal point, that crosses our line of sight and consequently has a higher Doppler beaming factor. A weak radio hot spot is indeed observed less than 1'' downstream of the bright radio hot spot, but has no optical or X-ray counterpart. By making use of the parsec-scale radio and the kiloparsec-scale radio/X-ray data, we derive constraints on the jet Lorentz factors ({Gamma}{sub jet}) and inclination angles ({theta}): for a constant jet speed from parsec to kiloparsec scales, we obtain a {Gamma}{sub jet} of {approx}70 for 0106+013 and {approx}40 for 3C 345. On relaxing this assumption, we derive a {Gamma}{sub jet} of {approx}2.5 for both the sources. Upper limits on {theta} of {approx}13 Degree-Sign are obtained for the two quasars. Broadband (radio-optical-X-ray) spectral energy distribution (SED) modeling of individual jet components in both quasars suggests that the optical emission is from the synchrotron mechanism, while the X-rays are produced via the inverse Compton mechanism from relativistically boosted cosmic microwave background seed photons. The locations of the upstream X-ray termination peaks strongly suggest that the sites of bulk jet deceleration lie upstream (by a few kiloparsecs) of the radio hot spots in these quasars. These regions are also the sites of shocks or magnetic field dissipation, which reaccelerate charged particles and produce high-energy optical and X-ray photons. This is consistent with the best-fit SED modeling parameters of magnetic field strength and electron power-law indices being higher in the jet termination regions compared to the cores. The shocked jet regions upstream of the radio hot spots, the kiloparsec-scale jet wiggles and a 'nose cone'-like jet structure in 0106+013, and the V-shaped radio structure in 3C 345, are all broadly consistent with instabilities associated with Poynting-flux-dominated jets. A greater theoretical understanding and more sensitive numerical simulations of jets spanning parsec to kiloparsec scales are needed, however, to make direct quantitative comparisons.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Metal-poor dwarf galaxies in the SIGRID galaxy sample. II. The electron temperature-abundance calibration and the parameters that affect it

In this paper, we use the Mappings photoionization code to explore the physical parameters that impact on the measurement of electron temperature and abundance in H II regions. In our previous paper, we presented observations and measurements of physical properties from the spectra of 17 H II regions in 14 isolated dwarf irregular galaxies from the SIGRID sample. Here, we analyze these observations further, together with three additional published data sets. We explore the effects of optical thickness, electron density, ionization parameter, ionization source, and non-equilibrium effects on the relation between electron temperature and metallicity. We present a standard model that fits the observed data remarkably well at metallicities between one-tenth and 1 solar. We investigate the effects of optically thin H II regions, and show that they can have a considerable effect on the measured electron temperature, and that there is evidence that some of the observed objects are optically thin. We look at the role of the ionization parameter and find that lower ionization parameter values give better fits at higher oxygen abundance. We show that higher pressures combined with low optical depth, and also κ electron energy distributions at low κ values, can generate the apparent high electron temperatures in low-metallicity H II regions, and that the former provides the better fit to observations. We examine the effects of these parameters on the strong line diagnostic methods. We extend this to three-dimensional diagnostic grids to confirm how well the observations are described by the grids.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

3D modeling of electric fields in the LUX detector

This work details the development of a three-dimensional (3D) electric field model for the LUX detector. The detector took data to search for weakly interacting massive particles (WIMPs) during two periods. After the first period completed, a time-varying non-uniform negative charge developed in the polytetrafluoroethylene (PTFE) panels that define the radial boundary of the detector's active volume. This caused electric field variations in the detector in time, depth and azimuth, generating an electrostatic radially-inward force on electrons on their way upward to the liquid surface. To map this behavior, 3D electric field maps of the detector's active volume were generated on a monthly basis. This was done by fitting a model built in COMSOL Multiphysics to the uniformly distributed calibration data that were collected on a regular basis. The modeled average PTFE charge density increased over the course of the exposure from -3.6 to -5.5 μC/m 2 . Here, from our studies, we deduce that the electric field magnitude varied locally while the mean value of the field of ~200 V/cm remained constant throughout the exposure. As a result of this work the varying electric fields and their impact on event reconstruction and discrimination were successfully modeled.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2019

Modeling the Connection between Subhalos and Satellites in Milky Way–like Systems

We develop a comprehensive and flexible model for the connection between satellite galaxies and dark matter subhalos in dark matter-only zoom-in simulations of Milky Way (MW)–mass host halos. We systematically identify the physical and numerical uncertainties in the galaxy–halo connection and simulations underlying our method, including (i) the influence of host halo properties; (ii) the relationship between satellite luminosities and subhalo properties, including the effects of reionization; (iii) the relationship between satellite and subhalo locations; (iv) the relationship between satellite sizes and subhalo properties, including the effects of tidal stripping; (v) satellite and subhalo disruption due to baryonic effects; and (vi) artificial subhalo disruption and orphan satellites. To illustrate our approach, we fit this model to the luminosity distribution of both classical MW satellites and those discovered in the Sloan Digital Sky Survey by performing realistic mock observations that depend on the luminosity, size, and distance of our predicted satellites, and we infer the total satellite population that will be probed by upcoming surveys. As a result, we argue that galaxy size and surface brightness modeling will play a key role in interpreting current and future observations, as the expected number of observable satellites depends sensitively on their surface brightness distribution.

79 ASTRONOMY AND ASTROPHYSICS↗