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 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

The Atacama Cosmology Telescope: likelihood for small-scale CMB data

The Atacama Cosmology Telescope has measured the angular power spectra of microwave fluctuations to arcminute scales at frequencies of 148 and 218 GHz, from three seasons of data. At small scales the fluctuations in the primordial Cosmic Microwave Background (CMB) become increasingly obscured by extragalactic foregounds and secondary CMB signals. We present results from a nine-parameter model describing these secondary effects, including the thermal and kinematic Sunyaev-Zel'dovich (tSZ and kSZ) power; the clustered and Poisson-like power from Cosmic Infrared Background (CIB) sources, and their frequency scaling; the tSZ-CIB correlation coefficient; the extragalactic radio source power; and thermal dust emission from Galactic cirrus in two different regions of the sky. In order to extract cosmological parameters, we describe a likelihood function for the ACT data, fitting this model to the multi-frequency spectra in the multipole range 500 < l < 10000. We extend the likelihood to include spectra from the South Pole Telescope at frequencies of 95, 150, and 220 GHz. Accounting for different radio source levels and Galactic cirrus emission, the same model provides an excellent fit to both datasets simultaneously, with χ{sup 2}/dof= 675/697 for ACT, and 96/107 for SPT. We then use the multi-frequency likelihood to estimate the CMB power spectrum from ACT in bandpowers, marginalizing over the secondary parameters. This provides a simplified 'CMB-only' likelihood in the range 500 < l < 3500 for use in cosmological parameter estimation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Does the first chaotic inflation model in supergravity provide the best fit to the Planck data?

I describe the first model of chaotic inflation in supergravity, which was proposed by Goncharov and the present author in 1983. The inflaton potential of this model has a plateau-type behavior V{sub 0}(1−(8/3) e{sup −√6|ϕ|}) at large values of the inflaton field. This model predicts n{sub s}=1−(2/N)≈0.967 and r=(4/(3N{sup 2}))≈4×10{sup −4}, in good agreement with the Planck data. I propose a slight generalization of this model, which allows to describe not only inflation but also dark energy and supersymmetry breaking.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

A viscous blast-wave model for high energy heavy-ion collisions

Employing a viscosity-based survival scale for initial geometrical perturbations formed in relativistic heavy-ion collisions, we model the radial flow velocity at freeze-out. Subsequently, we use the Cooper-Frye freeze-out prescription, with viscous corrections to the distribution function, to extract the transverse momentum dependence of particle yields and flow harmonics. We fit the model parameters for central collisions, by fitting the spectra of identified particles at the Large Hadron Collider (LHC), and estimate them for other centralities using simple hydrodynamic relations. We use the results of Monte Carlo Glauber model for initial eccentricities. We demonstrate that this improved viscous blast-wave model leads to good agreement with transverse momentum distribution of elliptic and triangular flow for all centralities and estimate the shear viscosity to entropy density ratio η/s ≃ 0.24 at the LHC.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2013

A CHARA ARRAY SURVEY OF CIRCUMSTELLAR DISKS AROUND NEARBY Be-TYPE STARS

We report on a high angular resolution survey of circumstellar disks around 24 northern sky Be stars. The K-band continuum survey was made using the CHARA Array long baseline interferometer (baselines of 30-331 m). The interferometric visibilities were corrected for the flux contribution of stellar companions in those cases where the Be star is a member of a known binary or multiple system. For those targets with good (u, v) coverage, we used a four-parameter Gaussian elliptical disk model to fit the visibilities and to determine the axial ratio, position angle, K-band photospheric flux contribution, and angular diameter of the disk's major axis. For the other targets with relatively limited (u, v) coverage, we constrained the axial ratio, inclination angle, and/or disk position angle where necessary in order to resolve the degeneracy between possible model solutions. We also made fits of the ultraviolet and infrared spectral energy distributions (SEDs) to estimate the stellar angular diameter and infrared flux excess of each target. The mean ratio of the disk diameter (measured in K-band emission) to stellar diameter (from SED modeling) is 4.4 among the 14 cases where we reliably resolved the disk emission, a value which is generally lower than the disk size ratio measured in the higher opacity H{alpha} emission line. We estimated the equatorial rotational velocity from the projected rotational velocity and disk inclination for 12 stars, and most of these stars rotate close to or at the critical rotational velocity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Full Modeling and parameter compression methods in configuration space for DESI 2024 and beyond

In the contemporary era of high-precision spectroscopic surveys, led by projects like DESI, there is an increasing demand for optimizing the extraction of cosmological information from clustering data. This work conducts a thorough comparison of various methodologies for modeling the full shape of the two-point statistics in configuration space. We investigate the performance of both direct fits (Full Modeling) and the parameter compression approaches (ShapeFit and Standard). We utilize the ABACUS-SUMMIT simulations, tailored to exceed DESI's precision requirements. Particularly, we fit the two-point statistics of three distinct tracers (LRG, ELG, and QSO), by employing a Gaussian Streaming Model in tandem with Convolution Lagrangian Perturbation Theory and Effective Field Theory. We explore methodological setup variations, including the range of scales, the set of galaxy bias parameters, the inclusion of the hexadecapole, as well as model extensions encompassing varying ns and allowing for w 0 w a CDM dark energy model. Throughout these varied explorations, while precision levels fluctuate and certain configurations exhibit tighter parameter constraints, our pipeline consistently recovers the parameter values of the mocks within 1σ in all cases for a 1-year DESI volume. Additionally, we compare the performance of configuration space analysis with its Fourier space counterpart using three models: PyBird, FOLPS and velocileptors, presented in companion papers. We find good agreement with the results from all these models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Color-magnitude distribution of face-on nearby galaxies in Sloan digital sky survey DR7

We have analyzed the distributions in the color-magnitude diagram (CMD) of a large sample of face-on galaxies to minimize the effect of dust extinctions on galaxy color. About 300,000 galaxies with log (a/b) < 0.2 and redshift z < 0.2 are selected from the Sloan Digital Sky Survey DR7 catalog. Two methods are employed to investigate the distributions of galaxies in the CMD, including one-dimensional (1D) Gaussian fitting to the distributions in individual magnitude bins and two-dimensional (2D) Gaussian mixture model (GMM) fitting to galaxies as a whole. We find that in the 1D fitting, two Gaussians are not enough to fit galaxies with the excess present between the blue cloud and the red sequence. The fitting to this excess defines the center of the green valley in the local universe to be (u – r){sub 0.1} = –0.121M {sub r,} 0{sub .1} – 0.061. The fraction of blue cloud and red sequence galaxies turns over around M {sub r,} {sub 0.1} ∼ –20.1 mag, corresponding to stellar mass of 3 × 10{sup 10} M {sub ☉}. For the 2D GMM fitting, a total of four Gaussians are required, one for the blue cloud, one for the red sequence, and the additional two for the green valley. The fact that two Gaussians are needed to describe the distributions of galaxies in the green valley is consistent with some models that argue for two different evolutionary paths from the blue cloud to the red sequence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

FLUORINE ABUNDANCES OF GALACTIC LOW-METALLICITY GIANTS

With abundances and 2{sigma} upper limits of fluorine (F) in seven metal-poor field giants, nucleosynthesis of stellar F at low metallicity is discussed. The measurements are derived from the HF(1-0) R9 line at 23358 A using near-infrared K-band high-resolution spectra obtained with CRIRES at the Very Large Telescope. The sample reaches lower metallicities than previous studies on F of field giants, ranging from [Fe/H] = -1.56 down to -2.13. Effects of three-dimensional model atmospheres on the derived F and O abundances are quantitatively estimated and shown to be insignificant for the program stars. The observed F yield in the form of [F/O] is compared with two sets of Galactic chemical evolution models, which quantitatively demonstrate the contribution of Type II supernova (SN II) {nu}-process and asymptotic giant branch/Wolf-Rayet stars. It is found that at this low-metallicity region, models cannot well predict the observed distribution of [F/O], while the observations are better fit by models considering an SN II {nu}-process with a neutrino energy of E {sub {nu}} = 3 Multiplication-Sign 10{sup 53} erg. Our sample contains HD 110281, a retrograde orbiting low-{alpha} halo star, showing a similar F evolution as globular clusters. This supports the theory that such halo stars are possibly accreted from dwarf galaxy progenitors of globular clusters in the halo.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A critical assessment of solutions to the galaxy diversity problem

Galactic rotation curves exhibit a diverse range of inner slopes. Observational data indicates that explaining this diversity may require a mechanism that correlates a galaxy's surface brightness with the central-most region of its dark matter halo. Here, in this work, we compare several concrete models that capture the relevant physics required to explain the galaxy diversity problem. We focus specifically on a Self-Interacting Dark Matter (SIDM) model with an isothermal core and two Cold Dark Matter (CDM) models with/without baryonic feedback. In contrast to the CDM case, the SIDM model can lead to the formation of an isothermal core in the halo, and is also mostly insensitive to baryonic feedback processes, which act on longer time-scales. Using rotation curves from 90 galaxies in the Spitzer Photometry & Accurate Rotation Curves (SPARC) catalog, we perform a comprehensive model comparison that addresses issues of statistical methodology from prior works. The best-fit halo models that we recover are consistent with standard CDM concentration-mass and abundance matching relations. We find that both the SIDM and feedback-affected CDM models are better than a CDM model with no feedback in explaining the rotation curves of low and high surface brightness galaxies in the sample. However, when compared to each other, there is no strong statistical preference for either the SIDM or the feedback-affected CDM halo model as the source of galaxy diversity in the SPARC catalog.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Robustness of cosmic birefringence measurement against Galactic foreground emission and instrumental systematics

The polarization of the cosmic microwave background (CMB) can be used to search for parity-violating processes like that predicted by a Chern-Simons coupling to a light pseudoscalar field. Such an interaction rotates E modes into E modes in the observed CMB signal through an effect known as cosmic birefringence. Even though isotropic birefringence can be confused with the rotation produced by a miscalibration of the detectors' polarization angles, the degeneracy between both effects is broken when Galactic foreground emission is used as a calibrator. Here, in this work, we use realistic simulations of the High-Frequency Instrument of the Planck mission to test the impact that Galactic foreground emission and instrumental systematics have on the recent birefringence measurements obtained through this technique. Our results demonstrate the robustness of the methodology against the miscalibration of polarization angles and other systematic effects, like intensity-to-polarization leakage, beam leakage, or cross-polarization effects. However, our estimator is sensitive to the EB correlation of polarized foreground emission. Here we propose to correct the bias induced by dust EB by modeling the foreground signal with templates produced in Bayesian component-separation analyses that fit parametric models to CMB data. Acknowledging the limitations of currently available dust templates like that of the Commander sky model, high-precision CMB data and a characterization of dust beyond the modified blackbody paradigm are needed to obtain a definitive measurement of cosmic birefringence in the future.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Three-band tight-binding model for monolayers of group-VIB transition metal dichalcogenides

We present a three-band tight-binding (TB) model for describing the low-energy physics in monolayers of group-VIB transition metal dichalcogenides MX2 (M = Mo, W; X = S, Se, Te). As the conduction- and valence-band edges are predominantly contributed by the dz2 , dxy, and dx2 y2 orbitals of M atoms, the TB model is constructed using these three orbitals based on the symmetries of the monolayers. Parameters of the TB model are fitted from the first-principles energy bands for all MX2 monolayers. The TB model involving only the nearest-neighbor M-M hoppings is sufficient to capture the band-edge properties in the K valleys, including the energy dispersions as well as the Berry curvatures. The TB model involving up to the third-nearest-neighbor M-M hoppings can well reproduce the energy bands in the entire Brillouin zone. Spin-orbit coupling in valence bands is well accounted for by including the on-site spin-orbit interactions ofM atoms. The conduction band also exhibits a small valley-dependent spin splitting which has an overall sign difference between MoX2 and WX2. We discuss the origins of these corrections to the three-band model. The three-band TB model developed here is efficient to account for low-energy physics in MX2 monolayers, and its simplicity can be particularly useful in the study of many-body physics and physics of edge states.

Liu, G. B.↗
DOE OSTI2024

Reduced-Order Modeling of Hydrogen Releases from Vent Stacks and with Wind

This report describes computational modeling in the HyRAM+ software for study of hydrogen behavior in two common scenarios. First, models of unignited plumes exiting a vent stack were considered. It was shown that entrainment and vent backpressure were major factors in plume physics. Second, HyRAM+ was extended to include wind effects on plumes by modifying the plume momentum balance and entrainment modeling. Use of the model showed that plume shape and length changed with wind speed and direction; in all cases, wind causes a shortening of the plume along the streamline. While the no-wind case in HyRAM+ has been validated and the newly developed wind model was fitted to very limited experimental data, more controlled experimental configurations would help validate the models and ensure accurate simulation of hydrogen plume behavior for vent stack releases or in wind.

08 HYDROGEN↗
DOE OSTI2017

Limb Darkening and Planetary Transits: Testing Center-to-limb Intensity Variations and Limb-darkening Directly from Model Stellar Atmospheres

The transit method, employed by Microvariability and Oscillation of Stars ( MOST ), Kepler , and various ground-based surveys has enabled the characterization of extrasolar planets to unprecedented precision. These results are precise enough to begin to measure planet atmosphere composition, planetary oblateness, starspots, and other phenomena at the level of a few hundred parts per million. However, these results depend on our understanding of stellar limb darkening, that is, the intensity distribution across the stellar disk that is sequentially blocked as the planet transits. Typically, stellar limb darkening is assumed to be a simple parameterization with two coefficients that are derived from stellar atmosphere models or fit directly. In this work, we revisit this assumption and compute synthetic planetary-transit light curves directly from model stellar atmosphere center-to-limb intensity variations (CLIVs) using the plane-parallel Atlas and spherically symmetric SAtlas codes. We compare these light curves to those constructed using best-fit limb-darkening parameterizations. We find that adopting parametric stellar limb-darkening laws leads to systematic differences from the more geometrically realistic model stellar atmosphere CLIV of about 50–100 ppm at the transit center and up to 300 ppm at ingress/egress. While these errors are small, they are systematic, and they appear to limit the precision necessary to measure secondary effects. Our results may also have a significant impact on transit spectra.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

A small-box approach to the local crystal structure of Y 3 NbO 7

The local crystal structure of Y 3 NbO 7 was probed using neutron total scattering. Five different structural models were fit to the experimental pair distribution function of this material up to 15 Å. Fits of defect fluorite (Fm$\overline{3}$m), pyrochlore (Fd$\overline{3}$m), and orthorhombic C weberite models (C222 1 ) were attempted based on prior structural investigations of Y 3 NbO 7 and related fluorite-based materials. Orthorhombic P weberite (P222 1 ) and monoclinic weberite models (P112 1 ) were derived via symmetry reduction of the orthorhombic C weberite structure. All models except monoclinic weberite failed to reproduce the experimental pair distribution function. Misfits were clearly visible for peaks arising from atom pairs within the first coordination shell (defect fluorite and pyrochlore) and beyond (orthorhombic C and P weberite). By contrast, monoclinic weberite (γ ≈ 91.8°) provided an excellent fit, with no major misfits observed. This unit cell features distorted YO 8 dodecahedra, YO 7 monocapped octahedra, and NbO 6 octahedra as building blocks. Comparison between the topologies of the cationic and anionic substructures in monoclinic weberite and its parent orthorhombic C structure revealed displacements of the metal atoms along the [100] direction and a less compact spatial distribution of oxygen atoms along the a axis. In conclusion, the monoclinic weberite model presented herein may serve as a starting point to revisit structural modeling of fluorite-based materials with challenging local structures.

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

Parameterized yields of semivolatile products from isoprene oxidation under different NOx levels: impacts of chemical aging and wall-loss of reactive gases

We developed a new parameterization for the yields of semi-volatile/intermediate volatility products from isoprene oxidation that took into account the chemical aging under different NOx levels, including functionalization and fragmentation reactions within the gas phase, oligomerization and photolysis within the condensed phase, as well as the loss of particles and gases to chamber walls. We used a box model to fit the dynamic evolution of measured secondary organic aerosol (SOA) concentrations from chamber experiments. The explicit inclusion of multigenerational gas phase chemistry and SOA oligomerization/photolysis during fitting itself had large impacts on the fitted yields that relate to the initial volatility distribution of isoprene oxidation products. During further multigenerational aging of these initial organics, we varied the branching ratio between functionalization and fragmentation. The best fits suggested increasing amounts of fragmentation with increasing initial NOx levels for mixed experiments. Using a box model to simulate evolution of SOA in the atmosphere at constant dilution rate (in absence of wall losses), we found that the yields determined by including different chemical scenarios during fitting changed predicted SOA mass concentrations substantially over aging timescale of 1-2 days. Our results demonstrated the need for better constraining smog-chamber-derived volatility distributions of both gas- and particle-phase SOA constituents with measurements.

Xing, Li↗
DOE OSTI2017

Estimating the reproductive number, total outbreak size, and reporting rates for Zika epidemics in South and Central America

As South and Central American countries prepare for increased birth defects from Zika virus outbreaks and plan for mitigation strategies to minimize ongoing and future outbreaks, understanding important characteristics of Zika outbreaks and how they vary across regions is a challenging and important problem. We developed a mathematical model for the 2015/2016 Zika virus outbreak dynamics in Colombia, El Salvador, and Suriname. We fit the model to publicly available data provided by the Pan American Health Organization, using Approximate Bayesian Computation to estimate parameter distributions and provide uncertainty quantification. The model indicated that a country-level analysis was not appropriate for Colombia. We then estimated the basic reproduction number to range between 4 and 6 for El Salvador and Suriname with a median of 4.3 and 5.3, respectively. We estimated the reporting rate to be around 16% in El Salvador and 18% in Suriname with estimated total outbreak sizes of 73,395 and 21,647 people, respectively. The uncertainty in parameter estimates highlights a need for research and data collection that will better constrain parameter ranges.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2012

WARM SPITZER OBSERVATIONS OF THREE HOT EXOPLANETS: XO-4b, HAT-P-6b, AND HAT-P-8b

We analyze Warm Spitzer/Infrared Array Camera observations of the secondary eclipses of three planets, XO-4b, HAT-P-6b, and HAT-P-8b. We measure secondary eclipse amplitudes at 3.6 {mu}m and 4.5 {mu}m for each target. XO-4b exhibits a stronger eclipse depth at 4.5 {mu}m than at 3.6 {mu}m, which is consistent with the presence of a temperature inversion. HAT-P-8b shows a stronger eclipse amplitude at 3.6 {mu}m and is best described by models without a temperature inversion. The eclipse depths of HAT-P-6b can be fitted with models with a small or no temperature inversion. We consider our results in the context of a postulated relationship between stellar activity and temperature inversion and a relationship between irradiation level and planet dayside temperature, as discussed by Knutson et al. and Cowan and Agol, respectively. Our results are consistent with these hypotheses, but do not significantly strengthen them. To measure accurate secondary eclipse central phases, we require accurate ephemerides. We obtain primary transit observations and supplement them with publicly available observations to update the orbital ephemerides of the three planets. Based on the secondary eclipse timing, we set upper boundaries for ecos ({omega}) for HAT-P-6b, HAT-P-8b, and XO-4b and find that the values are consistent with circular orbits.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Tape-in-tape-out (TITO): a new approach for in-situ contact resistance measurements in high temperature superconducting CORC® cables

One of the ongoing development challenges with ReBCO high-temperature superconducting (HTS) cables is normal zone initiation and local heating, which is associated with over-critical current flowing through and around local performance reductions in individual tapes. Although inter-tape contact resistances are well-reported for individual tapes and HTS cables, these measurements are a challenge in CORC® cables as current percolates through much of the cable. In this work, developments in tape-in-tape-out (TITO) automated individual tape powering experiments are presented, and a simplified modeling approach for current percolation in CORC® cables is developed. An optimization is formulated to fit the model parameters to a large set of TITO experiments on a single cable, allowing the layer-dependent inter-tape contact resistance to be extracted. Measurements are presented and discussed for a straight CORC® cable and a cable bent to a 152 mm and 76 mm radius. The approach provides both quantitative insight and facilitates better understanding of current sharing phenomena in HTS cables.

Teyber, Reed↗
DOE OSTI2012

Cosmological solutions in bimetric gravity and their observational tests

We obtain the general cosmological evolution equations for a classically consistent theory of bimetric gravity. Their analytic solutions are demonstrated to generically allow for a cosmic evolution starting out from a matter dominated FLRW universe and relaxing towards a de Sitter (anti-de Sitter) phase at late cosmic time. In particular, we examine a subclass of models which contain solutions that are able to reproduce the expansion history of the cosmic concordance model inspite of the nonlinear couplings of the two metrics. This is demonstrated explicitly by fitting these models to observational data from Type Ia supernovae, Cosmic Microwave Background and Baryon Acoustic Oscillations. In the appendix we comment on the relation to massive gravity.

79 ASTRONOMY AND ASTROPHYSICS↗