Search NASA⌕ Search

SEARCH · Search NASA

Results for “Markov Chain Monte Carlo Solar”

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.

Recognizing Chromospheric Objects via Markov Chain Monte Carlo

The solar chromosphere consists of three classes which contribute differentially to ultraviolet radiation reaching the earth. We describe a data set of solar images, means of segmenting the images into the constituent classes, and a novel high-level representation for compact objects based on a triangulated spatial membership function.

Markov Chain Monte Carlo Solar↗

Metallicities of Emission-Line Galaxies from HST ACS PEARS and HST WFC3 ERS Grism Spectroscopy at 0.6 is less than z is less than 2.4

Galaxies selected on the basis of their emission line strength. show low metallicities, regardless of their redshifts. We conclude this from a sample of faint galaxies at redshifts between 0.6 < z < 2.4, selected by their prominent emission lines in low resolution grism spectra in the optiCa.i with the Advanced Camera for Surveys (ACS) on the Hubble Space Telescope (HST) and in the near-infrared using Wide-Field Camera 3 (WFC3). Using a sample of 11 emission line galaxies (ELGs) at 0.6 < z < 2.4 with luminosities of -22 approx < MB approx -19 which have [OII], H-Beta, and [OIII] line flux measurements from the combination of two grism spectral surveys, we use the R23 method to derive the gas-phase oxygen abundances: 7.5 <12+log(0/H)<8.5. The galaxy stellar masses are derived using Bayesian based Markov Chain Monte Carlo (pi MC(exp 2)) fitting of their Spectral Energy Distribution (SED), and span the mass range 8.1 < log(M(stellar)/M(solar)) < 10.1. These galaxies show a mass-metal1icity (M-L) and Luminosity-Metallicity (LZ) relation, which is offset by -<J.6 dex in meta1licity at given absolute magnitude and stellar mass relative to the local SDSS galaxies, as well as continuum selected DEEP2 samples at similar redshifts. The emission-line selected galaxies most resemble the local "green peas" galaxies and Lyman-alpha galaxies at z approx = 0.3 and z approx = 2.3 in the M-Z and L-Z relations and their morphologies. The G - M(sub 20) morphology analysis shows that 10 out of 11 show disturbed morphology, even as the star-forming regions are compact. These galaxies may be intrinsically metal poor, being at early stages of formation, or the low metallicities may be due to gas infall and accretion due to mergers.

Xia, Lifang↗

Representing Solar Active Regions with Triangulations

The solar chromosphere consists of three classes which contribute differently to ultraviolet radiation reaching the earth. We describe a data set of solar images, means of segmenting the images into the constituent classes, and novel high-level representation for compact objects based on a triangulation spatial 'membership function'.

Triangulation chromosphere Markov chain Monte Carl↗

Alma Images the Eccentric HD 53143 Debris Disk

We present ALMA 1.3 mm observations of the HD 53143 debris disk—the first infrared or millimeter image produced of this ∼1 Gyr old solar analog. Previous HST STIS coronagraphic imaging did not detect flux along the minor axis of the disk, which could suggest a face-on geometry with two clumps of dust. These ALMA observations reveal a disk with a strikingly different structure. In order to fit models to the millimeter visibilities and constrain the uncertainties on the disk parameters, we adopt a Markov Chain Monte Carlo approach. This is the most eccentric debris disk observed to date with a forced eccentricity of 0.21 ± 0.02, nearly twice that of the Fomalhaut debris disk, and also displays an apocenter glow. Although this eccentric model fits the outer debris disk well, significant interior residuals remain, which may suggest a possible edge-on inner disk, which remains unresolved in these observations. Combined with the observed structure difference between HST and ALMA, these results suggest a potential previous scattering event or dynamical instability in this system. We also note that the stellar flux changes considerably over the course of our observations, suggesting flaring at millimeter wavelengths. Using simultaneous TESS observations, we determine the stellar rotation period to be 9.6 ± 0.1 days.

Radio astronomy↗

A Space-Based Observational Strategy for Characterizing the First Stars and Galaxies Using the Redshifted 21cm Global Spectrum

The redshifted 21 cm monopole is expected to be a powerful probe of the epoch of the first stars and galaxies(10 less than z less than 35). The global 21 cm signal is sensitive to the thermal and ionization state of hydrogen gas and thusprovides a tracer of sources of energetic photonsprimarily hot stars and accreting black holeswhich ionize andheat the high redshift intergalactic medium (IGM). This paper presents a strategy for observations of the globalspectrum with a realizable instrument placed in a low-altitude lunar orbit, performing night-time 40120 MHzspectral observations, while on the farside to avoid terrestrial radio frequency interference, ionospheric corruption,and solar radio emissions. The frequency structure, uniformity over large scales, and unpolarized state of theredshifted 21 cm spectrum are distinct from the spectrally featureless, spatially varying, and polarized emissionfrom the bright foregrounds. This allows a clean separation between the primordial signal and foregrounds. Forsignal extraction, we model the foreground, instrument, and 21 cm spectrum with eigenmodes calculated viaSingular Value Decomposition analyses. Using a Markov Chain Monte Carlo algorithm to explore the parameterspace defined by the coefficients associated with these modes, we illustrate how the spectrum can be measured andhow astrophysical parameters (e.g., IGM properties, first star characteristics) can be constrained in the presence offoregrounds using the Dark Ages Radio Explorer (DARE).

Burns, Jack O.↗

Estimating the Properties of Hard X-Ray Solar Flares by Constraining Model Parameters

We wish to better constrain the properties of solar flares by exploring how parameterized models of solar flares interact with uncertainty estimation methods. We compare four different methods of calculating uncertainty estimates in fitting parameterized models to Ramaty High Energy Solar Spectroscopic Imager X-ray spectra, considering only statistical sources of error. Three of the four methods are based on estimating the scale-size of the minimum in a hypersurface formed by the weighted sum of the squares of the differences between the model fit and the data as a function of the fit parameters, and are implemented as commonly practiced. The fourth method is also based on the difference between the data and the model, but instead uses Bayesian data analysis and Markov chain Monte Carlo (MCMC) techniques to calculate an uncertainty estimate. Two flare spectra are modeled: one from the Geostationary Operational Environmental Satellite X1.3 class flare of 2005 January 19, and the other from the X4.8 flare of 2002 July 23.We find that the four methods give approximately the same uncertainty estimates for the 2005 January 19 spectral fit parameters, but lead to very different uncertainty estimates for the 2002 July 23 spectral fit. This is because each method implements different analyses of the hypersurface, yielding method-dependent results that can differ greatly depending on the shape of the hypersurface. The hypersurface arising from the 2005 January 19 analysis is consistent with a normal distribution; therefore, the assumptions behind the three non- Bayesian uncertainty estimation methods are satisfied and similar estimates are found. The 2002 July 23 analysis shows that the hypersurface is not consistent with a normal distribution, indicating that the assumptions behind the three non-Bayesian uncertainty estimation methods are not satisfied, leading to differing estimates of the uncertainty. We find that the shape of the hypersurface is crucial in understanding the output from each uncertainty estimation technique, and that a crucial factor determining the shape of hypersurface is the location of the low-energy cutoff relative to energies where the thermal emission dominates. The Bayesian/MCMC approach also allows us to provide detailed information on probable values of the low-energy cutoff, Ec, a crucial parameter in defining the energy content of the flare-accelerated electrons. We show that for the 2002 July 23 flare data, there is a 95% probability that Ec lies below approximately 40 keV, and a 68% probability that it lies in the range 7-36 keV. Further, the low-energy cutoff is more likely to be in the range 25-35 keV than in any other 10 keV wide energy range. The low-energy cutoff for the 2005 January 19 flare is more tightly constrained to 107 +/- 4 keV with 68% probability.

X-rays↗

Exocometary Gas in the HD 181327 Debris Ring

An increasing number of observations have shown that gaseous debris discs are not an exception. However, until now, we only knew of cases around A stars. Here we present the first detection of 12CO (2-1) disc emission around an F star, HD 181327, obtained with the Atacama Large Millimeter/submillimeter Array (ALMA) observations at 1.3 mm. The continuum and CO emission are resolved into an axisymmetric disc with ring-like morphology. Using a Markov chain Monte Carlo method coupled with radiative transfer calculations, we study the dust and CO mass distribution. We find the dust is distributed in a ring with a radius of 86.0 +/- 0.4 au and a radial width of 23.2 +/- 1.0 au. At this frequency, the ring radius is smaller than in the optical, revealing grain size segregation expected due to radiation pressure. We also report on the detection of low-level continuum emission beyond the main ring out to approximately 200 au. We model the CO emission in the non-local thermodynamic equilibrium regime and we find that the CO is co-located with the dust, with a total CO gas mass ranging between 1.2 x 10(exp -6) solar mass and 2.9 x 10(exp -6) solar mass, depending on the gas kinetic temperature and collisional partners densities. The CO densities and location suggest a secondary origin, i.e. released from icy planetesimals in the ring. We derive a CO+CO2 cometary composition that is consistent with Solar system comets. Due to the low gas densities, it is unlikely that the gas is shaping the dust distribution.

exocometary gas↗

STIS Coronagraphic Imaging of Fomalhaut: Main Belt Structure and the Orbit of Fomalhaut b

We present new optical coronagraphic data of Fomalhaut obtained with HST/STIS in 2010 and 2012. Fomalhaut b is recovered at both epochs to high significance. The observations include the discoveries of tenuous nebulosity beyond the main dust belt detected to at least 209AU projected radius, and a approx. 50AU wide azimuthal gap in the belt northward of Fomalhaut b. The two epochs of Space Telescope Imaging Spectrograph (STIS) photometry exclude optical variability greater than 35%. A Markov chain Monte Carlo analysis demonstrates that the orbit of Fomalhaut b is highly eccentric, with e = 0.8 +/- 0.1, a = 177 +/- 68AU, and q = 32 +/- 24AU. Fomalhaut b is apsidally aligned with the belt and 90% of allowed orbits have mutual inclination <=36 deg. Fomalhaut b's orbit is belt crossing in the sky plane projection, but only 12% of possible orbits have ascending or descending nodes within a 25AU wide belt annulus. The high eccentricity invokes a dynamical history where Fomalhaut b may have experienced a significant dynamical interaction with a hypothetical planet Fomalhaut c, and the current orbital configuration may be relatively short-lived. The Tisserand parameter with respect to a hypothetical Fomalhaut planet at 30AU or 120AU lies in the range 2-3, similar to highly eccentric dwarf planets in our solar system. We argue that Fomalhaut b's minimum mass is that of a dwarf planet in order for a circumplanetary satellite system to remain bound to a sufficient radius from the planet to be consistent with the dust scattered light hypothesis. In the coplanar case, Fomalhaut b will collide with the main belt around 2032, and the subsequent emergent phenomena may help determine its physical nature.

Tisserand parameter↗