Using Data-Consistent Inversion to Build Population-Informed Priors for Bayesian Inference
Explore the source record for details and available documents.
SEARCH · Search NASA
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.
Explore the source record for details and available documents.
ABSTRACT We present a summary of our project that studies galaxies hosting type Ia supernova (SN Ia) at different redshifts. We present Gran Telescopio de Canarias (GTC) optical spectroscopy of six SN Ia host galaxies at redshift z ∼ 0.4–0.5. They are joined to a set of SN Ia host galaxies at intermediate-high redshift, which include galaxies from surveys SDSS and COSMOS. The final sample, after a selection of galaxy spectra in terms of signal-to-noise and other characteristics, consists of 680 galaxies with redshift in the range 0.04 < z < 1. We perform an inverse stellar population synthesis with the code fado to estimate the star formation and enrichment histories of this set of galaxies, simultaneously obtaining their mean stellar age and metallicity and stellar mass. After analysing the correlations among these characteristics, we look for possible dependencies of the Hubble diagram residuals and supernova features (luminosity, colour and strength parameter) on these stellar parameters. We find that the Hubble residuals show a clear dependence on the stellar metallicity weighted by mass with a slope of −0.061 mag dex−1, when represented in logarithmic scale, log 〈ZM/Z⊙〉. This result supports our previous findings obtained from gas oxygen abundances for local and SDSS-survey galaxies. Comparing with other works from the literature that also use the stellar metallicity, we find a similar value, but with more precision and a better significance (2.08 versus ∼ 1.1), due to the higher number of objects and wider range of redshift of our sample.
Resonance frequency mixing in nonlinear interaction of two intense optical waves with three energy levels of atomic system
Population of atomic levels by cascade, dielectric and three-body electronic recombination, discussing spontaneous transition, electron impact and RF spectra
Excited state and steady state populations in high pressure plasma with heavy particle collisional ionization and electron-ion recombination
In order to quantitatively predict nano- as well as other particle-size distributions, one needs to have both a mathematical model and estimates of the parameters that appear in these models. Here, we show how one can use Bayesian inversion to obtain statistical estimates for the parameters that appear in recently derived mechanism-enabled population balance models (ME-PBM) of nanoparticle growth. The Bayesian approach addresses the question of “how well do we know our parameters, along with their uncertainties?.” The results reveal that Bayesian inversion statistical analysis on an example, prototype $\mathrm{lr(0)_n}$ nanoparticle formation system allows one to estimate not just the most likely rate constants and other parameter values, but also their SDs, confidence intervals, and other statistical information. Moreover, knowing the reliability of the mechanistic model's parameters in turn helps inform one about the reliability of the proposed mechanism, as well as the reliability of its predictions. Importantly, the paper can also be seen as a tutorial with the additional goal of achieving a “Gold Standard” Bayesian inversion ME-PBM benchmark that others can use as a control to check their own use of this methodology for other systems of interest throughout nature. Overall, the results provide strong support for the hypothesis that there is substantial value in using a Bayesian inversion methodology for parameter estimation in particle formation systems.
Nitrogen oxides (NO(x)) are key actors in air quality and climate change. Satellite remote sensing of tropospheric NO2 has developed rapidly with enhanced spatial and temporal resolution since initial observations in 1995. We have developed an improved algorithm and retrieved tropospheric NO2 columns from Ozone Monitoring Instrument. Column observations of tropospheric NO2 from the nadir-viewing satellite sensors contain large contributions from the boundary layer due to strong enhancement of NO2 in the boundary layer. We infer ground-level NO2 concentrations from the OMI satellite instrument which demonstrate significant agreement with in-situ surface measurements. We examine how NO2 columns measured by satellite, ground-level NO2 derived from satellite, and NO(x) emissions obtained from bottom-up inventories relate to world's urban population. We perform inverse modeling analysis of NO2 measurements from OMI to estimate "top-down" surface NO(x) emissions, which are used to evaluate and improve "bottom-up" emission inventories. We use NO2 column observations from OMI and the relationship between NO2 columns and NO(x) emissions from a GEOS-Chem model simulation to estimate the annual change in bottom-up NO(x) emissions. The emission updates offer an improved estimate of NO(x) that are critical to our understanding of air quality, acid deposition, and climate change.
Synchronised populations of mouse C3H/10T-1/2 cells were obtained by a stringent mitotic dislodgment procedure. Mitotic cells rapidly attach and progress sequentially through the cell cycle. Irradiation (3 Gy of X rays) was carried out at intervals from 0 to 18 h after initiating cell cycle progression of the mitotic cells. Oncogenic transformation was enhanced 10-fold over cells irradiated soon after replating (G1 and S phases) for cells in a near 2 h period corresponding to cells in G2 phase but not in mitosis. The cell surviving fraction had a 2-1/2-fold variation with resistant peaks corresponding to the late G1 and late S phases. These findings provide experimental support for the hypothesis initiated by Rossi and Kellerer and developed by Brenner and Hall to explain the LET dependent inverse dose rate effect for oncogenic transformation.
Modeling biogeochemical fluxes in the marine plankton requires the application of factors for extrapolation of biomass indicators measured in the field (chlorophyll a, adenosine triphosphate, bacterial counts) to biomass carbon or nitrogen. These are often inferred from culture studies and are poorly constrained for natural populations. At least squares inverse method with a simple linear model constrains the values of several common indicator ratios, giving self-consistent solutions that provide useful information about the structure of the microbial community at our North Pacific Ocean study site (Station ALOHA (A Long-term Oligotrophic Habitat Assessment)). These results indicate that the fraction of the microbial biomass that is autotrophic (pigmented) is greater in the mixed layer than at the deep chlorophyll maximum layer and that heterotrophic bacteria are a significant but not necessarily predominant component of the microbial community in the euphotic zone.
In recent years there has been significant experimental effort aimed at studying the impact of 16 O neutron poisoning on the weak 𝑠 process in rotating massive stars. Improving the understanding of energy levels in 21 Ne is crucial to reducing uncertainties in the rates of the 𝛼-induced reactions on 17 O that determine the overall efficiency of the weak 𝑠 process. This paper reports on one such experiment: a study of the 20 Ne (𝑑,𝑝) 21 Ne reaction in inverse kinematics. Deduced spin-parity assignments were made based on adiabatic distorted-wave approximation analysis and neutron partial widths were estimated by comparison with a previous experiment. Of particular significance for nuclear astrophysics is the resulting neutron partial width for the 7820-keV energy level, estimated to be 12 200 (2900) eV, and an upper limit of ≤ 9300 eV for that of the 7749-keV energy level. These results are in tension with previous studies and therefore this paper also discusses the current state of the research into the astrophysically relevant energy levels and highlights both areas of agreement and areas of disagreement between this and various other studies that have investigated this nucleus.
The first high-resolution in-beam γ -ray spectroscopy is reported for the neutron-rich nucleus 41 Si, a tenant of the N = 28 island of inversion. Excited states were populated in the direct one-proton removal reaction from 42 P projectiles and pn removal from 43 P. Seven γ-ray transitions were observed, only one of which had been reported previously in the literature. Furthermore, this makes 41 Si the most neutron-rich odd-even N = 27 isotone with high-resolution excited-state information. For the one-proton removal, the measured partial cross-section distribution to all observed bound final states is contrasted with results from direct one-proton removal calculations that combine eikonal reaction dynamics with SDPF-MU shell-model spectroscopic factors and assume various possible initial states for the poorly known 42 P projectile. Rather distinct calculated cross-section distributions emerge that, in comparison to the new data, imply that the initial state in 42 P is most likely 3 – or 2 – rather than 1 – or 0 – , the predicted shell-model ground state of 42 P. It is further shown that the level scheme from the novel VS-IMSRG calculation closely agrees with the one of SDPF-MU, the most successful phenomenological shell-model effective interaction in describing the much discussed neighboring isotope 42 Si, perhaps cross-validating these complementary approaches on the quest to model rapid shell evolution away from the valley of β stability.
ABSTRACT We investigate the kinematical and dynamical properties of cluster galaxy populations classified according to their dominant source of gas ionization, namely: star-forming (SF) galaxies, optical active galactic nuclei (AGNs), mixed SF plus AGN ionization (transition objects, T), and quiescent (Q) galaxies. We stack 8892 member galaxies from 336 relaxed galaxy clusters to build an ensemble cluster and estimate the observed projected profiles of numerical density and velocity dispersion, $\sigma _P(R)$, of each galaxy population. The MAMPOSSt code and the Jeans equations inversion technique are used to constrain the velocity anisotropy profiles of the galaxy populations in both parametric and non-parametric ways. We find that Q (SF) galaxies display the lowest (highest) typical cluster-centric distances and velocity dispersion values. Transition galaxies are more concentrated and tend to exhibit lower velocity dispersion values than SF galaxies. Galaxies that host an optical AGN are as concentrated as Q galaxies but display velocity dispersion values similar to those of the SF population. MAMPOSSt is able to find equilibrium solutions that successfully recover the observed $\sigma _P(R)$ profile only for the Q, T, and AGN populations. We find that the orbits of all populations are consistent with isotropy in the inner regions, becoming increasingly radial with the distance from the cluster centre. These results suggest that Q galaxies are in equilibrium within their clusters, while SF galaxies have more recently arrived in the cluster environment. Finally, the T and AGN populations appear to be in an intermediate dynamical state between those of the SF and Q populations.
This paper investigates the non-Boltzmann modeling of the radiating atomic and molecular electronic states present in lunar-return shock-layers. The Master Equation is derived for a general atom or molecule while accounting for a variety of excitation and de-excitation mechanisms. A new set of electronic-impact excitation rates is compiled for N, O, and N2+, which are the main radiating species for most lunar-return shock-layers. Based on these new rates, a novel approach of curve-fitting the non-Boltzmann populations of the radiating atomic and molecular states is developed. This new approach provides a simple and accurate method for calculating the atomic and molecular non-Boltzmann populations while avoiding the matrix inversion procedure required for the detailed solution of the Master Equation. The radiative flux values predicted by the present detailed non-Boltzmann model and the approximate curve-fitting approach are shown to agree within 5% for the Fire 1634 s case.
The density and configurational changes of crystal dislocations during plastic deformation influence the mechanical properties of materials. These influences have become clearest in nanoscale experiments, in terms of strength, hardness and work hardening size effects in small volumes. The mechanical characterization of a model crystal may be cast as an inverse problem of deducing the defect population characteristics (density, correlations) in small volumes from the mechanical behavior. In this work, we demonstrate how a deep residual network can be used to deduce the dislocation characteristics of a sample of interest using only its surface strain profiles at small deformations, and then statistically predict the mechanical response of size-affected samples at larger deformations. As a testbed of our approach, we utilize high-throughput discrete dislocation simulations for systems of widths that range from nano- to micro- meters. We show that the proposed deep learning model significantly outperforms a traditional machine learning model, as well as accurately produces statistical predictions of the size effects in samples of various widths. By visualizing the filters in convolutional layers and saliency maps, we find that the proposed model is able to learn the significant features of sample strain profiles.
Abstract The inverse problem of extracting the stellar population content of galaxy spectra is analysed here from a basic standpoint based on information theory. By interpreting spectra as probability distribution functions, we find that galaxy spectra have high entropy, thus leading to a rather low effective information content. The highest variation in entropy is unsurprisingly found in regions that have been well studied for decades with the conventional approach. We target a set of six spectral regions that show the highest variation in entropy – the 4000 Å break being the most informative one. As a test case with real data, we measure the entropy of a set of high-quality spectra from the Sloan Digital Sky Survey, and contrast entropy-based results with the traditional method based on line strengths. The data are classified into star-forming (SF), quiescent (Q), and active galactic nucleus (AGN) galaxies, and show – independently of any physical model – that AGN spectra can be interpreted as a transition between SF and Q galaxies, with SF galaxies featuring a more diverse variation in entropy. The high level of entanglement complicates the determination of population parameters in a robust, unbiased way, and affects traditional methods that compare models with observations, as well as machine learning (especially deep learning) algorithms that rely on the statistical properties of the data to assess the variations among spectra. Entropy provides a new avenue to improve population synthesis models so that they give a more faithful representation of real galaxy spectra.
The extended jets of the microquasar SS 433 have been observed in optical, radio, X-ray, and recently very-high-energy γ-rays by High Altitude Water Cerenkov (HAWC). The detection of HAWC γ-rays with energies as great as 25 TeV motivates searches for high-energy γ-ray counterparts in the Fermi-Large Area Telescope (LAT) data in the 100 MeV–300 GeV band. In this paper, we report on the first joint analysis of Fermi-LAT and HAWC observations to study the spectrum and location of γ-ray emission from SS 433. Our analysis finds common emission sites of GeV-to-TeV γ-rays inside the eastern and western lobes of SS 433. The total flux above 1 GeV is $\sim 1\times {10}^{-10}\,{\mathrm{cm}}^{-2}\,{{\rm{s}}}^{-1}$ in both lobes. The γ-ray spectrum in the eastern lobe is consistent with inverse Compton emission by an electron population that is accelerated by jets. To explain both the GeV and TeV flux, the electrons need to have a soft intrinsic energy spectrum, or undergo a quick cooling process due to synchrotron radiation in a magnetized environment.
The results are presented of a study to determine the performance of genetic direct search algorithms in solving function optimization problems arising in the optimal and adaptive control areas. The findings indicate that: (1) genetic algorithms can outperform standard algorithms in multimodal and/or noisy optimization situations, but suffer from lack of gradient exploitation facilities when gradient information can be utilized to guide the search. (2) For large populations, or low dimensional function spaces, mutation is a sufficient operator. However for small populations or high dimensional functions, crossover applied in about equal frequency with mutation is an optimum combination. (3) Complexity, in terms of storage space and running time, is significantly increased when population size is increased or the inversion operator, or the second level adaptation routine is added to the basic structure.
The high resolution spectrum of comet P/Halley has been measured in the region of the CS emission. CS in comet Halley is produced by photodissociation of CS2 and subsequently reexcited by broad band solar radiation. An inversion procedure is introduced to extract the rotational population distribution from the dispersed fluorescence spectrum. Intensity anomalies in the wavelength region corresponding to P and Q branch emission are considered.