Likelihood informed dimension reduction for inverse problems in remote sensing of atmospheric constituent profiles

We use likelihood informed dimension reduction (LIS) (T. Cui et al. 2014) for inverting vertical profile information of atmospheric methane from ground based Fourier transform infrared (FTIR) measurements at Sodankyl\"a, Northern Finland. The measurements belong to the word wide TCCON network for greenhouse gas measurements and, in addition to providing accurate greenhouse gas measurements, they are important for validating satellite observations. LIS allows construction of an efficient Markov chain Monte Carlo sampling algorithm that explores only a reduced dimensional space but still produces a good approximation of the original full dimensional Bayesian posterior distribution. This in effect makes the statistical estimation problem independent of the discretization of the inverse problem. In addition, we compare LIS to a dimension reduction method based on prior covariance matrix truncation used earlier (S. Tukiainen et al. 2016)

Joint Bayesian Estimation of Quasar Continua and the Lyman-Alpha Forest Flux Probability Distribution Function

We present a new Bayesian algorithm making use of Markov Chain Monte Carlo sampling that allows us to simultaneously estimate the unknown continuum level of each quasar in an ensemble of high-resolution spectra, as well as their common probability distribution function (PDF) for the transmitted Ly$\alpha$ forest flux. This fully automated PDF regulated continuum fitting method models the unknown quasar continuum with a linear Principal Component Analysis (PCA) basis, with the PCA coefficients treated as nuisance parameters. The method allows one to estimate parameters governing the thermal state of the intergalactic medium (IGM), such as the slope of the temperature-density relation $\gamma-1$, while marginalizing out continuum uncertainties in a fully Bayesian way. Using realistic mock quasar spectra created from a simplified semi-numerical model of the IGM, we show that this method recovers the underlying quasar continua to a precision of $\simeq7\%$ and $\simeq10\%$ at $z=3$ and $z=5$, respectively. Given the number of principal component spectra, this is comparable to the underlying accuracy of the PCA model itself. Most importantly, we show that we can achieve a nearly unbiased estimate of the slope $\gamma-1$ of the IGM temperature-density relation with a precision of $\pm8.6\%$ at $z=3$, $\pm6.1\%$ at $z=5$, for an ensemble of ten mock high-resolution quasar spectra. Applying this method to real quasar spectra and comparing to a more realistic IGM model from hydrodynamical simulations would enable precise measurements of the thermal and cosmological parameters governing the IGM, albeit with somewhat larger uncertainties given the increased flexibility of the model.Comment: 21 pages (+ Appendix), accepted at Ap

Extraction of chemical information of suspensions using radiative transfer theory to remove multiple scattering effects : application to a model multicomponent system

The effectiveness of a scatter correction approach based on decoupling absorption and scattering effects through the use of the radiative transfer theory to invert a suitable set of measurements is studied by considering a model multicomponent suspension. The method was used in conjunction with partial least-squares regression to build calibration models for estimating the concentration of two types of analytes: an absorbing (nonscattering) species and a particulate (absorbing and scattering) species. The performances of the models built by this approach were compared with those obtained by applying empirical scatter correction approaches to diffuse reflectance, diffuse transmittance, and collimated transmittance measurements. It was found that the method provided appreciable improvement in model performance for the prediction of both types of analytes. The study indicates that, as long as the bulk absorption spectra are accurately extracted, no further empirical preprocessing to remove light scattering effects is required

Dynamical Masses of Young Star Clusters in NGC 4038/4039

In order to estimate the masses of the compact, young star clusters in the merging galaxy pair, NGC 4038/4039 (``the Antennae''), we have obtained medium and high resolution spectroscopy using ISAAC on VLT-UT1 and UVES on VLT-UT2 of five such clusters. The velocity dispersions were estimated using the stellar absorption features of CO at 2.29 microns and metal absorption lines at around 8500 \AA, including lines of the Calcium Triplet. The size scales and light profiles were measured from HST images. From these data and assuming Virial equilibrium, we estimated the masses of five clusters. The resulting masses range from 6.5 x 10^5 to 4.7 x 10^6 M_sun. These masses are large, factor of a few to more than 10 larger than the typical mass of a globular cluster in the Milky Way. The mass-to-light ratios for these clusters in the V- and K-bands in comparison with stellar synthesis models suggest that to first order the IMF slopes are approximately consistent with Salpeter for a mass range of 0.1 to 100 M_sun. However, the clusters show a significant range of possible IMF slopes or lower mass cut-offs and that these variations may correlate with the interstellar environment of the cluster. Comparison with the results of Fokker-Planck simulations of compact clusters with properties similar to the clusters studied here, suggest that they are likely to be long-lived and may lose a substantial fraction of their total mass. This mass loss would make the star clusters obtain masses which are comparable to the typical mass of a globular cluster.Comment: 16 pages, 12 figures, A&A accepte

Constraining mass ratio and extinction in the FU Orionis binary system with infrared integral field spectroscopy

We report low resolution near infrared spectroscopic observations of the eruptive star FU Orionis using the Integral Field Spectrograph Project 1640 installed at the Palomar Hale telescope. This work focuses on elucidating the nature of the faint source, located 0.5" south of FU Ori, and identified in 2003 as FU Ori S. We first use our observations in conjunction with published data to demonstrate that the two stars are indeed physically associated and form a true binary pair. We then proceed to extract J and H band spectro-photometry using the damped LOCI algorithm, a reduction method tailored for high contrast science with IFS. This is the first communication reporting the high accuracy of this technique, pioneered by the Project 1640 team, on a faint astronomical source. We use our low resolution near infrared spectrum in conjunction with 10.2 micron interferometric data to constrain the infrared excess of FU Ori S. We then focus on estimating the bulk physical properties of FU Ori S. Our models lead to estimates of an object heavily reddened, A_V =8-12, with an effective temperature of ~ 4000-6500 K . Finally we put these results in the context of the FU Ori N-S system and argue that our analysis provides evidence that FU Ori S might be the more massive component of this binary syste

The age-redshift relation for Luminous Red Galaxies in the Sloan Digital Sky Survey

We present a detailed analysis of 17,852 quiescent, Luminous Red Galaxies (LRGs) selected from Sloan Digital Sky Survey (SDSS) Data Release Seven (DR7) spanning a redshift range of 0.0 < z < 0.4. These galaxies are co-added into four equal bins of velocity dispersion and luminosity to produce high signal-to-noise spectra (>100A^{-1}), thus facilitating accurate measurements of the standard Lick absorption-line indices. In particular, we have carefully corrected and calibrated these indices onto the commonly used Lick/IDS system, thus allowing us to compare these data with other measurements in the literature, and derive realistic ages, metallicities ([Z/H]) and alpha-element abundance ratios ([alpha/Fe]) for these galaxies using Simple Stellar Population (SSP) models. We use these data to study the relationship of these galaxy parameters with redshift, and find little evidence for evolution in metallicity or alpha-elements (especially for our intermediate mass samples). This demonstrates that our subsamples are consistent with pure passive evolving (i.e. no chemical evolution) and represent a homogeneous population over this redshift range. We also present the age-redshift relation for these LRGs and clearly see a decrease in their age with redshift (5 Gyrs over the redshift range studied here) which is fully consistent with the cosmological lookback times in a concordance Lambda CDM universe. We also see that our most massive sample of LRGs is the youngest compared to the lower mass galaxies. We provide these data now to help future cosmological and galaxy evolution studies of LRGs, and provide in the appendices of this paper the required methodology and information to calibrate SDSS spectra onto the Lick/IDS system.Comment: 26 pages, with several appendices containing data. Accepted for publication in MNRA