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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2012

PHOTOMETRICALLY DERIVED MASSES AND RADII OF THE PLANET AND STAR IN THE TrES-2 SYSTEM

We measure the mass and radius of the star and planet in the TrES-2 system using 2.7 years of observations by the Kepler spacecraft. The light curve shows evidence for ellipsoidal variations and Doppler beaming on a period consistent with the orbital period of the planet with amplitudes of 2.79{sup +0.44}{sub -0.62} and 3.44{sup +0.32}{sub -0.37} parts per million (ppm), respectively, and a difference between the dayside and the nightside planetary flux of 3.41{sup +0.55}{sub -0.82} ppm. We present an asteroseismic analysis of solar-like oscillations on TrES-2A which we use to calculate the stellar mass of 0.94 {+-} 0.05 M{sub Sun} and radius of 0.95 {+-} 0.02 R{sub Sun }. Using these stellar parameters, a transit model fit and the phase-curve variations, we determine the planetary radius of 1.162{sup +0.020}{sub -0.024} R{sub Jup} and derive a mass for TrES-2b from the photometry of 1.44 {+-} 0.21 M{sub Jup}. The ratio of the ellipsoidal variation to the Doppler beaming amplitudes agrees to better than 2{sigma} with theoretical predications, while our measured planet mass and radius agree within 2{sigma} of previously published values based on spectroscopic radial velocity measurements. We measure a geometric albedo of 0.0136{sup +0.0022}{sub -0.0033} and an occultation (secondary eclipse) depth of 6.5{sup +1.7}{sub -1.8} ppm which we combined with the day/night planetary flux ratio to model the atmosphere of TrES-2b. We find that an atmosphere model that contains a temperature inversion is strongly preferred. We hypothesize that the Kepler bandpass probes a significantly greater atmospheric depth on the night side relative to the day side.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Cosmological inference from an emulator based halo model. II. Joint analysis of galaxy-galaxy weak lensing and galaxy clustering from HSC-Y1 and SDSS

Here, we present high-fidelity cosmology results from a blinded joint analysis of galaxy-galaxy weak lensing (ΔΣ) and projected galaxy clustering (w p ) measured from the Hyper Suprime-Cam Year-1 (HSC-Y1) data and spectroscopic Sloan Digital Sky Survey (SDSS) galaxy catalogs in the redshift range 0.15 < z < 0.7. We define luminosity-limited samples of SDSS galaxies to serve as the tracers of w p in three spectroscopic redshift bins, and as the lens samples for ΔΣ. For the ΔΣ measurements, we select a single sample of 4×10 6 source galaxies over 140 deg 2 from HSC-Y1 with photometric redshifts (photo z) greater than 0.75, enabling a better handle of photo- z errors by comparing the ΔΣ amplitudes for the three lens redshift bins. The deep, high-quality HSC-Y1 data enable significant detections of the ΔΣ signals, with integrated signal-to-noise ratio S/N ~ 15 in the range 3 ≤ R/[h –1 Mpc] ≤ 30 for the three lens samples, despite the small area coverage. For cosmological parameter inference, we use an input galaxy-halo connection model built on the dark emulator package (which uses an ensemble set of high-resolution N-body simulations and enables fast, accurate computation of the clustering observables) with a halo occupation distribution that includes nuisance parameters to marginalize over modeling uncertainties. We model the ΔΣ and wp measurements on scales from R≃3 and 2h –1 Mpc , respectively, up to 30 h –1 Mpc (therefore excluding the baryon acoustic oscillations information) assuming a flat Λ CDM cosmology, marginalizing over about 20 nuisance parameters and demonstrating the robustness of our results to them. With various tests using mock catalogs described in Miyatake et al. [preceding paper, Phys. Rev. D 106, 083519 (2022)], we show that any bias in the clustering amplitude S 8 ≡ σ 8 (Ω m /0.3) 0.5 due to uncertainties in the galaxy-halo connection is less than ~ 50 % of the statistical uncertainty of S 8 , unless the assembly biaseffect is unexpectedly large. Our best-fit models have S 8 = 0.795$_{-0.042}^{+0.049}$ (mode and 68% credible interval) for the flat Λ CDM model; we find tighter constraints on the quantity S 8 (α = 0.17)≡σ 8 (Ω m /0.3) 0.17 = 0.745$_{-0.031}^{+0.039}$.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Multistage Core Formation in Planetesimals Revealed by Numerical Modeling and Hf-W Chronometry of Iron Meteorites

Iron meteorites provide some of the most direct insights into the processes and timescales of core formation in planetesimals. Of these, group IVB irons stand out by having one of the youngest 182 Hf- 182 W model ages for metal segregation (2.9 ± 0.6 Ma after solar system formation), as well as the lowest bulk sulfur content and hence highest liquidus temperature. Here in this paper, using a new model for the internal evolution of the IVB parent body, we show that a single stage of metal-silicate separation cannot account for the complete melting of pure Fe metal at the relatively late time given by the Hf-W model age. Instead, a complex metal-silicate separation scenario is required that includes migration of partial silicate melts, formation of a shallow magma ocean, and core formation in two distinct stages of metal segregation. In the first stage, a protocore formed at ≈1.5 Ma via settling of metal particles in a mantle magma ocean, followed by metal segregation from a shallow magma ocean at ≈5.4 Ma. As these stages of metal segregation occurred at different times, the two metal fractions had different 182 W compositions. Consequently, the final 182W composition of the IVB core does not correspond to a single differentiation event, but represents the average composition of early- and late-segregated core fractions. Our best fit model indicates an ≈100 km radius for the IVB parent body and provides an accretion age of ≈0.1–0.5 Ma after solar system formation. The computed solidification time is, furthermore, consistent with the Re-Os age for crystallization of the IVB core.

58 GEOSCIENCES↗
DOE OSTI2013

ORIGIN OF THE DENSE CORE MASS FUNCTION IN CONTRACTING FILAMENTS

Mass functions of starless dense cores (CMFs) may arise from contraction and dispersal of core-forming filaments. In an illustrative model, a filament contracts radially by self-gravity, increasing the mass of its cores. During this contraction, FUV photoevaporation and ablation by shocks and winds disperse filament gas and limit core growth. The stopping times of core growth are described by a waiting-time distribution. The initial filament column density profile and the resulting CMF each match recent Herschel observations in detail. Then low-mass cores have short growth ages and arise from the innermost filament gas, while massive cores have long growth ages and draw from more extended filament gas. The model fits the initial density profile and CMF best for mean core density 2 Multiplication-Sign 10{sup 4} cm{sup -3} and filament dispersal timescale 0.5 Myr. Then the typical core mass, radius, mean column density, and contraction speed are respectively 0.8 solar masses, 0.06 pc, 6 Multiplication-Sign 10{sup 21} cm{sup -2}, and 0.07 km s{sup -1}, also in accord with observed values.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

A Coupled Mathematical Model of the Intracellular Replication of Dengue Virus and the Host Cell Immune Response to Infection

Dengue virus (DV) is a positive-strand RNA virus of the Flavivirus genus. It is one of the most prevalent mosquito-borne viruses, infecting globally 390 million individuals per year. The clinical spectrum of DV infection ranges from an asymptomatic course to severe complications such as dengue hemorrhagic fever (DHF) and dengue shock syndrome (DSS), the latter because of severe plasma leakage. Given that the outcome of infection is likely determined by the kinetics of viral replication and the antiviral host cell immune response (HIR) it is of importance to understand the interaction between these two parameters. In this study, we use mathematical modeling to characterize and understand the complex interplay between intracellular DV replication and the host cells' defense mechanisms. We first measured viral RNA, viral protein, and virus particle production in Huh7 cells, which exhibit a notoriously weak intrinsic antiviral response. Based on these measurements, we developed a detailed intracellular DV replication model. We then measured replication in IFN competent A549 cells and used this data to couple the replication model with a model describing IFN activation and production of IFN stimulated genes (ISGs), as well as their interplay with DV replication. By comparing the cell line specific DV replication, we found that host factors involved in replication complex formation and virus particle production are crucial for replication efficiency. Regarding possible modes of action of the HIR, our model fits suggest that the HIR mainly affects DV RNA translation initiation, cytosolic DV RNA degradation, and naïve cell infection. We further analyzed the potential of direct acting antiviral drugs targeting different processes of the DV lifecycle in silico and found that targeting RNA synthesis and virus assembly and release are the most promising anti-DV drug targets.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Data-guided Multi-Map variables for ensemble refinement of molecular movies

Driving molecular dynamics simulations with data-guided collective variables offer a promising strategy to recover thermodynamic information from structure-centric experiments. In this study, the three-dimensional electron density of a protein, as it would be determined by cryo-EM or x-ray crystallography, is used to achieve simultaneously free-energy costs of conformational transitions and refined atomic structures. Unlike previous density-driven molecular dynamics methodologies that determine only the best map-model fits, our work employs the recently developed Multi-Map methodology to monitor concerted movements within equilibrium, non-equilibrium, and enhanced sampling simulations. Construction of all-atom ensembles along the chosen values of the Multi-Map variable enables simultaneous estimation of average properties, as well as real-space refinement of the structures contributing to such averages. Using three proteins of increasing size, we demonstrate that biased simulation along the reaction coordinates derived from electron densities can capture conformational transitions between known intermediates. The simulated pathways appear reversible with minimal hysteresis and require only low-resolution density information to guide the transition. The induced transitions also produce estimates for free energy differences that can be directly compared to experimental observables and population distributions. The refined model quality is superior compared to those found in the Protein Data Bank. We find that the best quantitative agreement with experimental free-energy differences is obtained using medium resolution density information coupled to comparatively large structural transitions. Practical considerations for probing the transitions between multiple intermediate density states are also discussed.

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

THE GALAXY-HALO/SUBHALO CONNECTION: MASS RELATIONS AND IMPLICATIONS FOR SOME SATELLITE OCCUPATIONAL DISTRIBUTIONS

We infer the local stellar-to-halo/subhalo mass relations (MRs) for central and satellite galaxies separately. Our statistical method is an extension of the abundance matching, halo occupation distribution, and conditional stellar mass function formalisms. We constrain the model using several combinations of observational data, consisting of the total galaxy stellar mass function (GSMF), its decomposition into centrals and satellites, and the projected two-point correlation functions (2PCFs) measured in different stellar mass (M{sub *}) bins. In addition, we use the {Lambda}CDM halo and subhalo mass functions. The differences among the resulting MRs are within the model-fit uncertainties (which are very small, smaller than the intrinsic scatter between galaxy and halo mass), no matter what combination of data are used. This shows that matching abundances or occupational numbers is equivalent, and that the GSMFs and 2PCFs are tightly connected. We also constrain the values of the intrinsic scatter around the central-halo (CH) and satellite-subhalo (SS) MRs assuming them to be constant: {sigma}{sub c} = 0.168 {+-} 0.051 dex and {sigma}{sub s} = 0.172 {+-} 0.057 dex, respectively. The CH and SS MRs are actually different, in particular when we take the subhalo mass at the present-day epoch instead of at their accretion time. When using the MRs for studying the satellite population (e.g., in the Milky Way, MW), the SS MR should be chosen instead of the average one. Our model allows one to calculate several population statistics. We find that the central galaxy M{sub *} is not on average within the mass distribution of the most massive satellite, even for cluster-sized halos, i.e., centrals are not a mere realization of the high end of the satellite mass function; however for >3 Multiplication-Sign 10{sup 13} M{sub Sun} halos, {approx}15% of centrals could be. We also find that the probabilities of MW-sized halos of having N Magellanic Cloud (MC) sized satellites agree well with observational measures; for a halo mass of 2 Multiplication-Sign 10{sup 12} M{sub Sun }, the probability to have two MCs is 5.4%, but if we exclude those systems with satellites larger than the MCs, then the probability decreases to <2.2%.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Constraining the properties of gaseous halos via cross-correlations of upcoming galaxy surveys and thermal Sunyaev-Zel’dovich maps

The thermal Sunyaev-Zel’dovich (tSZ) effect induces a Compton-$\mathit{y}$ distortion in cosmic microwave background (CMB) temperature maps that is sensitive to a line-of-sight integral of the ionized gas pressure. By correlating the positions of galaxies with maps of the Compton-$\mathit{y}$ distortion, one can probe baryonic feedback processes and study the thermodynamic properties of a significant fraction of the gas in the Universe. Using a model fitting approach, we forecast how well future galaxy and CMB surveys will be able to measure these correlations, and show that powerful constraints on halo pressure profiles can be obtained. In this work, our forecasts are focused on correlations between galaxies and halos identified by the upcoming Dark Energy Spectroscopic Instrument survey and tSZ maps from the Simons Observatory and CMB-S4 experiments, but have general applicability to other surveys, such as the Large Synoptic Survey Telescope. We include prescriptions for observational systematics, such as halo miscentering and halo mass bias, demonstrating several important degeneracies with pressure profile parameters. Assuming modest priors on these systematics, we find that measurements of halo-$\mathit{y}$ and galaxy-$\mathit{y}$ correlations with future surveys yield tight constraints on the pressure profiles of group-scale dark matter halos, and enable current feedback models to either be confirmed or ruled out.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A Qualitative Strategy for Fusion of Physics into Empirical Models for Process Anomaly Detection

To facilitate the automated online monitoring of power plants, a systematic and qualitative strategy for anomaly detection is presented. This strategy is essential to provide credible reasoning on why and when an empirical versus hybrid (i.e., physics-supported) approach should be used and to determine the ideal mix of these two approaches for a defined anomaly detection scope. Empirical methods are usually based on pattern, statistical, and causal inference. Hybrid methods include the use of physics models to train and test data methods, reduce data dimensionality, reduce data-model complexity, augment data, and reduce empirical uncertainty; hybrid methods also include the use of data to tune physics models. The presented strategy is driven by key decision points related to data relevance, simple modeling feasibility, data inference, physics-modeling value, data dimensionality, physics knowledge, method of validation, performance, data availability, and suitability for training and testing, cause-effect, entropy inference, and model fitting. The strategy is demonstrated through a pilot use case for the application of anomaly detection to capture a valve packing leak at the high-pressure coolant injection system of a nuclear power plant.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2012

FIRST SCIENCE RESULTS FROM SOFIA/FORCAST: THE MID-INFRARED VIEW OF THE COMPACT H II REGION W3A

The massive star-forming region W3 was observed with the faint object infrared camera for the SOFIA telescope as part of the Short Science program. The 6.4, 6.6, 7.7, 19.7, 24.2, 31.5, and 37.1 {mu}m bandpasses were used to observe the emission of polycyclic aromatic hydrocarbon (PAH) molecules, very small grains, and big grains. Optical depth and color temperature maps of W3A show that IRS2 has blown a bubble devoid of gas and dust of {approx}0.05 pc radius. It is embedded in a dusty shell of ionized gas that contributes 40% of the total 24 {mu}m emission of W3A. This dust component is mostly heated by far-ultraviolet, rather than trapped Ly{alpha} photons. This shell is itself surrounded by a thin ({approx}0.01 pc) photodissociation region where PAHs show intense emission. The infrared spectral energy distribution (SED) of three different zones located at 8'', 20'', and 25'' from IRS2 shows that the peak of the SED shifts toward longer wavelengths, when moving away from the star. Adopting the stellar radiation field for these three positions, DUSTEM model fits to these SEDs yield a dust-to-gas mass ratio in the ionized gas similar to that in the diffuse interstellar medium (ISM). However, the ratio of the IR-to-UV opacity of the dust in the ionized shell is increased by a factor of {approx_equal}3 compared to the diffuse ISM.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Characterization of Final State Interaction Strength in Plastic Scintillator by Muon-Neutrino Charged Current Charged Pion Production

Precise knowledge of neutrino-nucleus interactions is increasingly important as neutrino oscillation measurements transition into the systematics-limited era. In addition to modifying the initial interaction, the nuclear medium can scatter and absorb the interaction by-products through final state interactions, changing the types and kinematic distributions of particles seen by the detector. Recent neutrino pion production data from MiniBooNE is inconsistent with the final state interaction strength predicted by models and theoretical calculations, and some models fit best to the MiniBooNE data only after removing final state interactions entirely. This thesis presents a measurement of dσ/dT π and dσ/dθ π for muon-neutrino charged current charged pion production in the MINER A scintillator tracker. MINER A is a neutrino-nucleus scattering experiment installed in the few-GeV NuMI beam line at Fermilab. The analysis is limited to neutrino energies between 1.5-10 GeV. Dependence on invariant hadronic mass W is studied through two versions of the analysis that impose the limits W < 1.4 GeV and W < 1.8 GeV. The lower limit on W increases compatibility with the MiniBooNE pion data. The shapes of the differential cross sections, which depend strongly on the nature of final state interactions, are compared to Monte Carlo and theoretical predictions. It is shown that the measurements presented in this thesis favor models that contain final state interactions. Additionally, a variety of neutrino-nucleus interaction models are shown to successfully reproduce the thesis measurements, while simultaneously failing to describe the shape of the MiniBooNE data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Kinetic understanding of nitrogen supply condition on biosynthesis of polyhydroxyalkanoate from benzoate by Pseudomonas putida KT2440

Nitrogen supply is critical to the synthesis of intracellular PHA in various bacteria. However, the specific role of the nitrogen in synthesizing PHA from benzoate, a lignin model compound use for the study of bacteria catabolism of aromatics, is still not clear. Here, two culture conditions were maintained for Pseudomonas putida KT2440 to produce PHA using benzoate as a carbon source. Under nitrogen-limited and surplus conditions, the accumulation of PHA was to 37.3% and 0.25% of cell dry weight, respectively. A model fit to the kinetics of biomass growth and PHA accumulation showed good agreement with data. GC–MS and NMR showed that PHA contained six hydroxyl fatty acid monomers under nitrogen-limited conditions, while two monomers were identified under nitrogen surplus conditions. The average molecular weight of PHA increased after the nitrogen source was exhausted. Lastly, these results provide a promising strategy for optimization of lignin to PHA yields.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2014

Nature versus nurture: Luminous blue variable nebulae in and near massive stellar clusters at the galactic center

Three luminous blue variables (LBVs) are located in and near the Quintuplet Cluster at the Galactic center: the Pistol Star, G0.120-0.048, and qF362. We present imaging at 19, 25, 31, and 37 μm of the region containing these three LBVs, obtained with SOFIA using FORCAST. We argue that Pistol and G0.120-0.048 are identical 'twins' that exhibit contrasting nebulae due to the external influence of their different environments. Our images reveal the asymmetric, compressed shell of hot dust surrounding the Pistol Star and provide the first detection of the thermal emission from the symmetric, hot dust envelope surrounding G0.120-0.048. However, no detection of hot dust associated with qF362 is made. Dust and gas composing the Pistol nebula are primarily heated and ionized by the nearby Quintuplet Cluster stars. The northern region of the Pistol nebula is decelerated due to the interaction with the high-velocity (2000 km s{sup –1}) winds from adjacent Wolf-Rayet Carbon (WC) stars. From fits to the spectral energy distribution (SED) of the Pistol nebula with the DustEM code we determine that the Pistol nebula is composed of a distribution of very small, transiently heated grains (10 to ∼ 35 Å) having a total dust mass of 0.03 M {sub ☉}, and that it exhibits a gradient of decreasing grain size from south to north due to differential sputtering by the winds from the WC stars. The total IR luminosity of the Pistol nebula is 5.2 × 10{sup 5} L {sub ☉}. Dust in the G0.120-0.048 nebula is primarily heated by the central star; however, the nebular gas is ionized externally by the Arches Cluster. Unlike the Pistol nebula, the G0.120-0.048 nebula is freely expanding into the surrounding medium. A grain size distribution identical to that of the non-sputtered region of the Pistol nebula satisfies the constraints placed on the G0.120-0.048 nebula from DustEM model fits to its SED and implies a total dust mass of 0.021 M {sub ☉}. The total IR luminosity of the G0.120-0.048 nebula is ∼10{sup 5} L {sub ☉}. From Paschen-α and 6 cm observations we determine a total gas mass of 9.3 M {sub ☉} and 6.2 M {sub ☉} for the Pistol and G0.120-0.048 nebulae, respectively. Given the independent dust and gas mass estimates we find that the Pistol and G0.120-0.048 nebulae exhibit similar gas-to-dust mass ratios of 310{sub −52}{sup +77} and 293{sub −101}{sup +73}, respectively. Both nebulae share identical size scales (∼0.7 pc) which suggests that they have similar dynamical timescales of ∼10{sup 4} yr, assuming a shell expansion velocity of v {sub exp} = 60 km s{sup –1}.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Ordered magnetic fields around the 3C 84 central black hole

Context. 3C 84 is a nearby radio source with a complex total intensity structure, showing linear polarisation and spectral patterns. A detailed investigation of the central engine region necessitates the use of very-long-baseline interferometry (VLBI) above the hitherto available maximum frequency of 86 GHz.Aims. Using ultrahigh resolution VLBI observations at the currently highest available frequency of 228 GHz, we aim to perform a direct detection of compact structures and understand the physical conditions in the compact region of 3C 84.Methods. We used Event Horizon Telescope (EHT) 228 GHz observations and, given the limited (u, v)-coverage, applied geometric model fitting to the data. Furthermore, we employed quasi-simultaneously observed, ancillary multi-frequency VLBI data for the source in order to carry out a comprehensive analysis of the core structure.Results. We report the detection of a highly ordered, strong magnetic field around the central, supermassive black hole of 3C 84. The brightness temperature analysis suggests that the system is in equipartition. We also determined a turnover frequency of νm = (113 ± 4) GHz, a corresponding synchrotron self-absorbed magnetic field of BSSA = (2.9 ± 1.6) G, and an equipartition magnetic field of Beq = (5.2 ± 0.6) G. Three components are resolved with the highest fractional polarisation detected for this object (mnet = (17.0 ± 3.9)%). The positions of the components are compatible with those seen in low-frequency VLBI observations since 2017–2018. We report a steeply negative slope of the spectrum at 228 GHz. We used these findings to test existing models of jet formation, propagation, and Faraday rotation in 3C 84.Conclusions. The findings of our investigation into different flow geometries and black hole spins support an advection-dominated accretion flow in a magnetically arrested state around a rapidly rotating supermassive black hole as a model of the jet-launching system in the core of 3C 84. However, systematic uncertainties due to the limited (u, v)-coverage, however, cannot be ignored. Our upcoming work using new EHT data, which offer full imaging capabilities, will shed more light on the compact region of 3C 84.Key words: techniques: high angular resolution / techniques: interferometric / galaxies: active / galaxies: individual: NGC 1275 / galaxies: jets

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Investigating Hα, UV, and IR star-formation rate diagnostics for a large sample of z ∼ 2 galaxies

We use a sample of 262 spectroscopically confirmed star-forming galaxies at redshifts 2.08⩽z⩽2.51 to compare Hα, ultraviolet (UV), and IR star formation rate (SFR) diagnostics and to investigate the dust properties of the galaxies. At these redshifts, the Hα line shifts to the K{sub s} band. By comparing K{sub s}-band photometry to underlying stellar population model fits to other UV, optical, and near-infrared data, we infer the Hα flux for each galaxy. We obtain the best agreement between Hα- and UV-based SFRs if we assume that the ionized gas and stellar continuum are reddened by the same value and that the Calzetti attenuation curve is applied to both. Aided with MIPS 24 μm data, we find that an attenuation curve steeper than the Calzetti curve is needed to reproduce the observed IR/UV ratios of galaxies younger than 100 Myr. Furthermore, using the bolometric SFR inferred from the UV and mid-IR data (SFR{sub IR}+SFR{sub UV}), we calculated the conversion between the Hα luminosity and SFR to be (7.5±1.3)×10{sup −42} for a Salpeter initial mass function, which is consistent with the Kennicutt conversion. The derived conversion factor is independent of any assumption of the dust correction and is robust to stellar population model uncertainties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Blanco DECam Bulge Survey (BDBS)

The Blanco DECam Bulge Survey (BDBS) has imaged more than 200 square degrees of the southern Galactic bulge, providing photometry in the ugrizy filters for ~250 million unique stars. The presence of a strong foreground disk population, along with complex reddening and extreme image crowding, has made it difficult to constrain the presence of young and intermediate age stars in the bulge population. We employed an accurate cross-match of BDBS with the latest data release (EDR3) from the Gaia mission, matching more than 140 million sources with BDBS photometry and Gaia EDR3 photometry and astrometry. We relied on Gaia EDR3 astrometry, without any photometric selection, to produce clean BDBS bulge colour-magnitude diagrams (CMDs). Gaia parallaxes were used to filter out bright foreground sources, and a Gaussian mixture model fit to Galactic proper motions could identify stars kinematically consistent with bulge membership. We applied this method to 127 different bulge fields of 1 deg 2 each, with |$\ell$| ≤ 9.5 ° and –9.5 ° ≤ b ≤ –2.5 °. The astrometric cleaning procedure removes the majority of blue stars in each field, especially near the Galactic plane, where the ratio of blue to red stars is .10%, increasing to values ~20% at higher Galactic latitudes. We rule out the presence of a widespread population of stars younger than 2 Gyr. The vast majority of blue stars brighter than the turnoff belong to the foreground population, according to their measured astrometry. We introduce the distance between the observed red giant branch bump and the red clump as a simple age proxy for the dominant population in the field, and we confirm the picture of a predominantly old bulge. Further work is needed to apply the method to estimate ages to fields at higher latitudes, and to model the complex morphology of the Galactic bulge. We also produce transverse kinematic maps, recovering expected patterns related to the presence of the bar and of the X-shaped nature of the bulge.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

A scripting interface for GSAS-II

The GSAS-II suite is a complete crystallographic analysis package for all types of X-ray and neutron diffraction data, suited for fitting models ranging from simple materials to proteins, and designed around an integrated graphical user interface. This article describes a new module within GSAS-II , GSASIIscript­able , which provides access to the GSAS-II data structures and an extensible mechanism to run a significant fraction of GSAS-II functionality, particularly for powder diffraction applications, from within Python scripts or directly from commands in a shell/batch script. In conclusion, this allows parallelization in a high-performance computing environment; near ideal speedup is seen with up to 240 simultaneous processes.

47 OTHER INSTRUMENTATION↗
DOE OSTI2016

TU-FG-BRB-03: Basis Vector Model Based Method for Proton Stopping Power Estimation From Experimental Dual Energy CT Data

Purpose: This work aims at reducing the uncertainty in proton stopping power (SP) estimation by a novel combination of a linear, separable basis vector model (BVM) for stopping power calculation (Med Phys 43:600) and a statistical, model-based dual-energy CT (DECT) image reconstruction algorithm (TMI 35:685). The method was applied to experimental data. Methods: BVM assumes the photon attenuation coefficients, electron densities, and mean excitation energies (I-values) of unknown materials can be approximated by a combination of the corresponding quantities of two reference materials. The DECT projection data for a phantom with 5 different known materials was collected on a Philips Brilliance scanner using two scans at 90 kVp and 140 kVp. The line integral alternating minimization (LIAM) algorithm was used to recover the two BVM coefficient images using the measured source spectra. The proton stopping powers are then estimated from the Bethe-Bloch equation using electron densities and I-values derived from the BVM coefficients. The proton stopping powers and proton ranges for the phantom materials estimated via our BVM based DECT method are compared to ICRU reference values and a post-processing DECT analysis (Yang PMB 55:1343) applied to vendorreconstructed images using the Torikoshi parametric fit model (tPFM). Results: For the phantom materials, the average stopping power estimations for 175 MeV protons derived from our method are within 1% of the ICRU reference values (except for Teflon with a 1.48% error), with an average standard deviation of 0.46% over pixels. The resultant proton ranges agree with the reference values within 2 mm. Conclusion: Our principled DECT iterative reconstruction algorithm, incorporating optimal beam hardening and scatter corrections, in conjunction with a simple linear BVM model, achieves more accurate and robust proton stopping power maps than the post-processing, nonlinear tPFM based DECT analysis applied to conventional reconstructions of low and high energy scans. Funding Support: NIH R01CA 75371; NCI grant R01 CA 149305.

60 APPLIED LIFE SCIENCES↗