Search

This work presents a new detection of H$_2$ absorption arising in a high velocity cloud (HVC) associated with either the Milky Way or the Large Magellanic Cloud (LMC). The absorber was found in an archival Far Ultraviolet Spectroscopic Explorer spectrum of the LMC star Sk-70$^\circ$32. This is the fifth well-characterized H$_2$ absorber to be found in the Milky Way's halo and the second such absorber outside the Magellanic Stream and Bridge. The absorber has a local standard of rest central velocity of $+$140 km s$^{-1}$ and a H$_2$ column density of $10^{17.5}$ cm$^{-2}$. It is most likely part of a cool and relatively dense inclusion ($T\approx 75$ K, $n_{\rm H}\sim 100$ cm$^{-3}$) in a warmer and more diffuse halo cloud. This halo cloud may be part of a still-rising Milky Way Galactic fountain flow or an outflow from the Large Magellanic Cloud., Comment: Accepted to ApJ. 13 pages, 6 figures

This study addresses how the incidence rate of strong O VI absorbers in a galaxy's circumgalactic medium (CGM) depends on galaxy mass and, independently, on the amount of star formation in the galaxy. We use HST/COS absorption spectroscopy of quasars to measure O VI absorption within 400 projected kpc and 300 km s$^{-1}$ of 52 $M_{*}\sim 10^{10}$ $M_\odot$ galaxies. The galaxies have redshifts $0.12, Comment: 12 pages, 5 figures. Submitted to the AAS Journals

We combine 126 new galaxy-O VI absorber pairs from the CGM$^2$ survey with 123 pairs drawn from the literature to examine the simultaneous dependence of the column density of O VI absorbers ($N_{\rm O VI}$) on galaxy stellar mass, star formation rate, and impact parameter. The combined sample consists of 249 galaxy-O VI absorber pairs covering $z=0$-$0.6$, with host galaxy stellar masses $M^*=10^{7.8}$-$10^{11.2}$ $M_\odot$ and galaxy-absorber impact parameters $R_\perp=0$-$400$ proper kiloparsecs. In this work, we focus on the variation of $N_{\rm O VI}$ with galaxy mass and impact parameter among the star-forming galaxies in the sample. We find that the average $N_{\rm O VI}$ within one virial radius of a star-forming galaxy is greatest for star-forming galaxies with $M^*=10^{9.2}$-$10^{10}$ $M_\odot$. Star-forming galaxies with $M^*$ between $10^{8}$ and $10^{11.2}$ $M_\odot$ can explain most O VI systems with column densities greater than 10$^{13.5}$ cm$^{-2}$. 60% of the O VI mass associated with a star-forming galaxy is found within one virial radius and 35% is found between one and two virial radii. In general, we find that some departure from hydrostatic equilibrium in the CGM is necessary to reproduce the observed O VI amount, galaxy mass dependence, and extent. Our measurements serve as a test set for CGM models over a broad range of host galaxy masses., 25 pages, 9 figures. Accepted to ApJ

We analyse 2,015 mid-infrared (MIR) spectra of galaxies observed with Spitzer's Infrared Spectrograph, including objects with growing super-massive black holes and objects where most of the infrared emission originates from newly formed stars. We determine if and how accreting super-massive black holes at the centre of galaxies -- known as active galactic nuclei (AGN) -- heat and ionize their host galaxies' dust and molecular gas. We use four MIR diagnostics to estimate the contribution of the AGN to the total MIR emission. We refer to galaxies whose AGN contribute more than 50 per cent of the total MIR emission as AGN-dominated. We compare the relative strengths of PAH emission features and find that PAH grains in AGN-dominated sources have a wider range of sizes and fractional ionizations than PAH grains in non-AGN dominated sources. We measure rotational transitions of H_2 and estimate H_2 excitation temperatures and masses for individual targets, H_2 excitation temperatures for spectra stacked by their AGN contribution to the MIR, and the H_2 excitation temperature distributions via a hierarchical Bayesian model. We find an average 200 K difference between the excitation temperatures of the H_2 S(5) and H_2 S(7) pure rotational molecular hydrogen transition pair in AGN-dominated versus non-AGN dominated galaxies. Our findings suggest that AGN impact the interstellar medium of their host galaxies., Comment: Published in MNRAS, 20 pages, 16 figures, link to full tables and stacked spectra (http://pages.jh.edu/~elambri1/MIR_AGN_MolecularGas/)

Quasar outflows carry mass, momentum and energy into the surrounding environment, and have long been considered a potential key factor in regulating the growth of supermassive black holes and the evolution of their host galaxies. A crucial parameter for understanding the origin of these outflows and measuring their influence on their host galaxies is the distance (R) between the outflow gas and the galaxy center. While R has been measured in a number of individual galaxies, its distribution remains unknown. Here we report the distributions of R and the kinetic luminosities of quasars outflows, using the statistical properties of broad absorption line variability in a sample of 915 quasars from the Sloan Digital Sky Surveys. The mean and standard deviation of the distribution of R are 10^{1.4+/-0.5} parsecs. The typical outflow distance in this sample is tens of parsec, which is beyond the theoretically predicted location (0.01 ~ 0.1 parsecs) where the accretion disc line-driven wind is launched, but is smaller than the scales of most outflows that are derived using the excited state absorption lines. The typical value of the mass-flow rate is of tens to a hundred solar masses per year, or several times the accretion rate. The typical kinetic-to-bolometric luminosity ratio is a few per cent, indicating that outflows are energetic enough to influence the evolution of their host galaxies., Accepted by Nature Astronomy, 8 pages and 2 figures(main text), 8 supplementary figures

In this work we demonstrate that the Perseus Arm is not a continuous structure of molecular gas in the second quadrant. We first show that the observed, distanced-resolved velocity structure of the Galaxy in the outer disk is capable of creating illusory spiral arms, as was first proposed by Burton (1971). Second, we measure the distances to a collection of CO clouds at velocities consistent with the Perseus arm with $135^\circ < l < 160^\circ$. We find these distances using 3D dust maps from Green et al. (2019). We determine that these molecular cloud do not preferentially lie at the distance of a purported Perseus arm, but rather extend over 3 kpc in distance, with some evidence for a closer, high pitch angle structure between 1 and 1.5 kpc away. Finally, we demonstrate that velocity perturbations of the amplitude found near the Perseus arm can wreak havoc on our interpretation of the longitude-velocity diagram for more than half of the Milky Way disk., 14 pages, 7 figures. Accepted to the Astrophysical Journal

The metallicity and gas density dependence of interstellar depletions, the dust-to-gas (D/G), and dust-to-metal (D/M) ratios have important implications for how accurately we can trace the chemical enrichment of the universe; either by using FIR dust emission as a tracer of the ISM; or by using spectroscopy of damped Lyman-alpha systems (DLAs) to measure chemical abundances over a wide range of redshifts. We collect and compare large samples of depletion measurements in the Milky Way (MW), LMC (Z=0.5 Zsun), and SMC (Z=0.2 Zsun). The relation between the depletions of different elements do not strongly vary between the three galaxies, implying that abundance ratios should trace depletions accurately down to 20% solar metallicity. From the depletions, we derive D/G and D/M. The D/G increases with density, consistent with the more efficient accretion of gas-phase metals onto dust grains in the denser ISM. For log N(H) > 21 cm^-2, the depletion of metallicity tracers (S, Zn) exceeds -0.5 dex, even at 20% solar metallicity. The gas fraction of metals increases from the MW to the LMC (factor 3) and SMC (factor 6), compensating the reduction in total heavy element abundances and resulting in those three galaxies having the same neutral gas-phase metallicities. The D/G derived from depletions are a factor of 2 (LMC) and 5 (SMC) higher than the D/G derived from FIR, 21 cm, and CO emission, likely due to the combined uncertainties on the dust FIR opacity and on the depletion of carbon and oxygen., 26 pages, 12 figures

Photodissociation regions (PDRs) are parts of the ISM consisting of predominantly neutral gas, located at the interface between H II regions and molecular clouds. The physical conditions within these regions show variations on very short spatial scales, and therefore PDRs constitute ideal laboratories for investigating the properties and evolution of dust grains. We have mapped IC 63 at high resolution from the UV to the NIR (275 nm to 1.6 $\mu$m), using the Hubble Space Telescope WFC3. Using a Bayesian SED fitting tool, we simultaneously derive a set of stellar ($T_\text{eff}$, $\log(g)$, distance) and extinction ($A_V$, $R_V$) parameters for 520 background stars. We present maps of $A_V$ and $R_V$ with a resolution of 25 arcsec based on these results. The extinction properties vary across the PDR, with values for $A_V$ between 0.5 and 1.4 mag, and a decreasing trend in $R_V$, going from 3.7 at the front of the nebula to values as low as 2.5 further in. This provides evidence for evolution of the dust optical properties. We fit two modified blackbodies to the MIR and FIR SED, obtained by combining the $A_V$ map with data from Spitzer and Herschel. We derive effective temperatures (30 K and 227 K) and the ratio of opacities at 160 $\mu$m to V band $\kappa_{160} / \kappa_V$ ($7.0 \times 10^{-4}$ and $2.9 \times 10^{-9}$) for the two dust populations. Similar fits to individual pixels show spatial variations of $\kappa_{160} / \kappa_{V}$. The analysis of our HST data, combined with these Spitzer and Herschel data, provides the first panchromatic view of dust within a PDR., Comment: 28 pages, 20 figures, accepted for publication in The Astrophysical Journal

We present initial results from the \textit{COS and Gemini Mapping the Circumgalactic Medium} (\mbox{CGMCGM} $\equiv$ CGM$^{2}$) survey. The CGM$^{2}$ survey consists of 1689 galaxies, all with high-quality Gemini GMOS spectra, within 1 Mpc of twenty-two $z \lesssim 1$ quasars, all with S/N$\sim$10 {\emph{HST/COS}} G130M$+$G160M spectra. For 572 of these galaxies having stellar masses $10^{7} M_{\odot} < M_{\star} < 10^{11} M_{\odot}$ and $z \lesssim 0.5$, we show that the \ion{H}{1} covering fraction above a threshold of \NHI$>10^{14} $cm$^{-2}$ is $\gtrsim 0.5$ within 1.5 virial radii ($R_{\rm vir} \sim R_{200m}$). We examine the \ion{H}{1} kinematics and find that the majority of absorption lies within $\pm$ 250 km s$^{-1}$ of the galaxy systemic velocity. We examine \ion{H}{1} covering fractions over a range of impact parameters to infer a characteristic size of the CGM, $R^{14}_{\rm CGM}$, as a function of galaxy mass. $R^{14}_{\rm CGM}$ is the impact parameter at which the probability of observing an absorber with \NHI $>$ 10$^{14}$ cm$^{-2}$ is $>$ 50\%. In this framework, the radial extent of the CGM of $M_{\star} > 10^{9.9} M_{\odot}$ galaxies is $R^{14}_{\rm CGM} = 346^{+57}_{-53}$ kpc or $R^{14}_{\rm CGM} \simeq 1.2R_{\rm vir}$. Intermediate-mass galaxies with $10^{9.2} < M_{\star}/M_{\odot} < 10^{9.9}$ have an extent of $R^{14}_{\rm CGM} = 353^{+64}_{-50}$ kpc or $R^{14}_{\rm CGM} \simeq 2.4R_{\rm vir}$. Low-mass galaxies, $M_{\star} < 10^{9.2} M_{\odot}$, show a smaller physical scale $R^{14}_{\rm CGM} = 177_{-65}^{+70}$ kpc and extend to $R^{14}_{\rm CGM} \simeq 1.6R_{\rm vir}$. Our analysis suggests that using $R_{\rm vir}$ as a proxy for the characteristic radius of the CGM likely underestimates its extent., 30 pages, 14 figures. Submitted to ApJ, Comments welcome

For over 60 years, the scientific community has studied actively growing central super-massive black holes (active galactic nuclei -- AGN) but fundamental questions on their genesis remain unanswered. Numerical simulations and theoretical arguments show that black hole growth occurs during short-lived periods ($\sim$ 10$^{7}$ -10$^{8}$ yr) of powerful accretion. Major mergers are commonly invoked as the most likely dissipative process to trigger the rapid fueling of AGN. If the AGN-merger paradigm is true, we expect galaxy mergers to coincide with black hole accretion during a heavily obscured AGN phase (N$_H$ $ > 10^{23}$ cm$^{-2}$). Starting from one of the largest samples of obscured AGN at 0.5 $, Comment: 19 pages, 11 figures, accepted in ApJ

Significant galaxy mergers throughout cosmic time play a fundamental role in theories of galaxy evolution. The widespread usage of human classifiers to visually assess whether galaxies are in merging systems remains a fundamental component of many morphology studies. Studies that employ human classifiers usually construct a control sample, and rely on the assumption that the bias introduced by using humans will be evenly applied to all samples. In this work, we test this assumption and develop methods to correct for it. Using the standard binomial statistical methods employed in many morphology studies, we find that the merger fraction, error, and the significance of the difference between two samples are dependent on the intrinsic merger fraction of any given sample. We propose a method of quantifying merger biases of individual human classifiers and incorporate these biases into a full probabilistic model to determine the merger fraction and the probability of an individual galaxy being in a merger. Using 14 simulated human responses and accuracies, we are able to correctly label a galaxy as ''merger'' or ''isolated'' to within 1\% of the truth. Using 14 real human responses on a set of realistic mock galaxy simulation snapshots our model is able to recover the pre-coalesced merger fraction to within 10\%. Our method can not only increase the accuracy of studies probing the merger state of galaxies at cosmic noon, but also can be used to construct more accurate training sets in machine learning studies that use human classified data-sets., 23 pages, 16 figures, accepted for publication in ApJ, https://github.com/elambrid/merger_or_not

Absorption line spectroscopy is a powerful way of measuring properties of stars and the interstellar medium. Absorption spectra are often analyzed manually, an approach that limits reproducibility and which cannot practically be applied to modern datasets consisting of thousands or even millions of spectra. Simultaneous probabilistic modeling of absorption features and continuum shape is a promising approach for automating this analysis. Existing implementations of this approach use numerical methods such as Markov chain Monte Carlo (MCMC) to marginalize over the continuum parameters. Numerical marginalization over large numbers of continuum parameters is too slow to be convenient for exploratory analysis, can increase the dimensionality of an inference problem beyond the capacity of simple MCMC samplers, and is in general impractical for the analysis of large datasets. When continua are parameterized as linear functions such as polynomials or splines, it is possible to reduce continuum parameter marginalization to an integral over a multivariate normal distribution, which has a known closed form. In addition to speeding up probabilistic modeling, analytic marginalization makes it trivial to marginalize over continuum parameterizations and to combine continuum description marginalization with optimization for absorption line parameters. These new possibilities allow automatic, probabilistically justified continuum placement in analyses of large spectroscopic datasets. We compare the accuracy to within which absorption line parameters can be recovered using different continuum placement methods and find that marginalization is in many cases an improvement over other methods. We implement analytic marginalization over linear continuum parameters in the open-source package amlc., Submitted to AAS journals

We use deep $HST$ WFC3/IR imaging to study the Initial Mass Function (IMF) of the ultra faint dwarf galaxy Coma Berenices (Com Ber). Our observations reach the lowest stellar mass ever probed in a resolved galaxy, with 50\% completeness at $\sim 0.17$ M$_{\odot}$. Unresolved background galaxies however limit our purity below $\sim 0.23$ M$_{\odot}$. If modeled with a single power law, we find that the IMF slope is $-1.45^{+0.29}_{-0.3}$ (68\% credible intervals), compared to a Milky Way value of $-2.3$. For a broken power law, we obtain a low-mass slope of $-1.18_{-0.33}^{+0.49}$, a high-mass slope of $-1.88_{-0.49}^{+0.43}$ and a break mass of $0.57_{-0.08}^{+0.12}$ M$_{\odot}$, compared to $-1.3$, $-2.3$ and 0.5 M$_{\odot}$ for a Kroupa IMF. For a log-normal IMF model we obtain values of $0.33_{-0.16}^{+0.15}$ M$_{\odot}$ for the location parameter and of $0.68_{-0.12}^{+0.17}$ for $\sigma$ (0.22 M$_{\odot}$ and 0.57 for the Chabrier system IMF). All three parametrizations produce similar agreement with the data. Our results agree with previous analysis of shallower optical HST data. However analysis of similar optical data of other dwarfs finds IMFs significantly more bottom-light than in the Milky Way. These results suggest two, non mutually exclusive, possibilities: that the discrepancy of the dwarf galaxies IMF with respect to the Milky Way is, at least partly, an artifact of using a single power law model, and that there is real variance in the IMF at low masses between the currently studied nearby dwarfs, with Com Ber being similar to the Milky Way, but other dwarfs differing significantly., Comment: ApJ accepted

In this work, we derive a spatially resolved map of the line-of-sight velocity of the interstellar medium and use it, along with a second map of line-of-sight velocity from Paper I of this series, to determine the nature of gaseous spiral structure in the Milky Way. This map is derived from measurements of the 1.527 $\mu$m diffuse interstellar band (DIB) in stellar spectra from the APOGEE survey and covers the nearest 4-5 kpc of the Northern Galactic plane. We cross-check this new DIB-based line-of-sight velocity map with the map derived in Paper I and find that they agree. We then compare these maps with line-of-sight velocity maps derived from simulations of quasi-stationary density wave spiral structure and dynamic, or material, spiral structure in a Milky Way-like galaxy. While none of the maps derived from these simulations is an exact match to the measured velocity field of the Milky Way, the measurements are more consistent with simulations of dynamic spiral structure. In the dynamic spiral structure simulation that best matches the measurements, the Perseus spiral arm is being disrupted., Comment: Submitted to AAS journals

International audience; Not Available

The dust extinction curve is a critical component of many observational programs and an important diagnostic of the physics of the interstellar medium. Here we present new measurements of the dust extinction curve and its variation towards tens of thousands of stars, a hundred-fold larger sample than in existing detailed studies. We use data from the APOGEE spectroscopic survey in combination with ten-band photometry from Pan-STARRS1, 2MASS, and WISE. We find that the extinction curve in the optical through infrared is well characterized by a one-parameter family of curves described by R(V). The extinction curve is more uniform than suggested in past works, with sigma(R(V)) = 0.18, and with less than one percent of sight lines having R(V) > 4. Our data and analysis have revealed two new aspects of Galactic extinction: first, we find significant, wide-area variations in R(V) throughout the Galactic plane. These variations are on scales much larger than individual molecular clouds, indicating that R(V) variations must trace much more than just grain growth in dense molecular environments. Indeed, we find no correlation between R(V) and dust column density up to E(B-V) ~ 2. Second, we discover a strong relationship between R(V) and the far-infrared dust emissivity., 26 pages, 20 figures, comments welcome

Using deep observations obtained with the Advanced Camera for Surveys (ACS) on board the Hubble Space Telescope (HST), we demonstrate that the sub-solar stellar initial mass function (IMF) of six ultra-faint dwarf Milky Way satellites (UFDs) is more bottom light than the IMF of the Milky Way disk. Our data have a lower-mass limit of ~0.45 M_⊙, while the upper limit is ~0.8 M_⊙, set by the turnoff mass of these old, metal-poor systems. If formulated as a single power law, we obtain a shallower IMF slope than the Salpeter value of −2.3, ranging from −1.01 for Leo IV to −1.87 for Bootes I. The significance of these deviations depends on the galaxy and is typically 95% or more. When modeled as a log-normal, the IMF fit results in a higher peak mass than in the Milky Way disk, but a Milky Way disk value for the characteristic system mass (~0.22 M_⊙) is excluded at only 68% significance, and only for some UFDs in the sample. We find that the IMF slope correlates well with the galaxy mean metallicity, and to a lesser degree, with the velocity dispersion and the total mass. The strength of the observed correlations is limited by shot noise in the number of observed stars, but future space-based missions like the James Webb Space Telescope (JWST) and the Wide-Field Infrared Survey Telescope ( WFIRST) will enhance both the number of dwarf Milky Way satellites that can be studied in such detail and the observation depth for individual galaxies.

We study the stellar mass Tully-Fisher relation (TFR, stellar mass versus rotation velocity) for a morphologically blind selection of emission line galaxies in the field at redshifts 0.1 $, Comment: Accepted to MNRAS. 14 pages, 7 figures

We present a study of the composition of gas and dust in the Large and Small Magellanic Clouds (LMC and SMC, together -- the MCs) as measured by UV absorption spectroscopy. We have measured P II and Fe II along 85 sightlines toward the MCs using archival FUSE observations. For 16 of those sightlines, we have measured Si II, Cr II, and Zn II from new HST COS observations. We have combined these measurements with H I and H$_2$ column densities and reference stellar abundances from the literature to derive gas-phase abundances, depletions, and gas-to-dust ratios (GDRs). 80 of our 84 P measurements and 13 of our 16 Zn measurements are depleted by more than 0.1 decades, suggesting that P and Zn abundances are not accurate metallicity indicators at and above the metallicity of the SMC. The maximum P and Zn depletions are the same in the MW, LMC, and SMC. Si, Cr, and Fe are systematically less depleted in the SMC than in the MW or LMC. The minimum Si depletion in the SMC is consistent with zero. Our depletion-derived GDRs broadly agree with GDRs from the literature. The GDR varies from location to location within a galaxy by a factor of up to 2 in the LMC and up to 5 in the SMC. This variation is evidence of dust destruction and/or growth in the diffuse neutral phase of the interstellar medium., Comment: Submitted to ApJ

The nature, composition, abundance, and size distribution of dust in galaxies is determined by the rate at which it is created in the different stellar sources and destroyed by interstellar shocks. Because of their extensive wavelength coverage, proximity, and nearly face-on geometry, the Magellanic Clouds (MCs) provide a unique opportunity to study these processes in great detail. In this paper we use the complete sample of supernova remnants (SNRs) in the MCs to calculate the lifetime and destruction efficiencies of silicate and carbon dust in these galaxies. We find dust lifetimes of 22 +- 13 Myr (30 +- 17 Myr) for silicate (carbon) grains in the LMC, and 54 +- 32 Myr (72 +- 43 Myr) for silicate (carbon) grains in the SMC. The significantly shorter lifetimes in the MCs, as compared to the Milky Way, are explained as the combined effect of their lower total dust mass, and the fact that the dust-destroying isolated SNe in the MCs seem to be preferentially occurring in regions with higher than average dust-to-gas (D2G) mass ratios. We also calculate the supernova rate and the current star formation rate in the MCs, and use them to derive maximum dust injection rates by asymptotic giant branch stars and core collapse supernovae. We find that the injection rates are an order of magnitude lower than the dust destruction rates by the SNRs. This supports the conclusion that, unless the dust destruction rates have been considerably overestimated, most of the dust must be reconstituted from surviving grains in dense molecular clouds. More generally, we also discuss the dependence of the dust destruction rate on the local D2G mass ratio, the ambient gas density and metallicity, as well as the application of our results to other galaxies and dust evolution models., 15 pages, 8 figures, 5 tables, accepted to ApJ

The dust properties in the Large and Small Magellanic Clouds are studied using the HERITAGE Herschel Key Project photometric data in five bands from 100 to 500 micron. Three simple models of dust emission were fit to the observations: a single temperature blackbody modified by a power- law emissivity (SMBB), a single temperature blackbody modified by a broken power-law emissivity (BEMBB), and two blackbodies with different temperatures, both modified by the same power-law emissivity (TTMBB). Using these models we investigate the origin of the submm excess; defined as the submillimeter (submm) emission above that expected from SMBB models fit to observations < 200 micron. We find that the BEMBB model produces the lowest fit residuals with pixel-averaged 500 micron submm excesses of 27% and 43% for the LMC and SMC, respectively. Adopting gas masses from previous works, the gas-to-dust ratios calculated from our the fitting results shows that the TTMBB fits require significantly more dust than are available even if all the metals present in the interstellar medium (ISM) were condensed into dust. This indicates that the submm excess is more likely to be due to emissivity variations than a second population of colder dust. We derive integrated dust masses of (7.3 +/- 1.7) x 10^5 and (8.3 +/- 2.1) times 10^4 M(sun) for the LMC and SMC, respectively. We find significant correlations between the submm excess and other dust properties; further work is needed to determine the relative contributions of fitting noise and ISM physics to the correlations., 20 pages, 9 figures, accepted for publication in the ApJ

We have developed a method for deriving the distribution of the Milky Way's interstellar medium as a function of longitude, latitude, distance and line-of-sight velocity. This method takes as input maps of reddening as a function of longitude, latitude, and distance and maps of line emission as a function of longitude, latitude, and line-of-sight velocity. We have applied this method to datasets covering much of the Galactic plane. The output of this method correctly reproduces the line-of-sight velocities of high-mass star forming regions with known distances from Reid et al. (2014) and qualitatively agrees with results from the Milky Way kinematics literature. These maps will be useful for measuring flows of gas around the Milky Way's spiral arms and into and out of giant molecular clouds., Comment: 11 pages, 5 figures, accepted for publication in The Astronomical Journal

We report rigorous quantum five-dimensional (5D) calculations of the coupled translation-rotation (TR) energy levels and wave functions of an H2 molecule, in the ground (ν = 0) and vibrationally excited (ν = 1) states, confined inside the octahedral interstitial site of solid C60 with S6 symmetry. Translational and rotational excitations of H2 in this nanocavity have been measured by the inelastic neutron scattering (INS) and infrared (IR) spectroscopy, enabling direct comparison between theory and experiment. A pairwise additive 5D intermolecular potential energy surface (PES) was employed in the calculations. The quantum calculations cover the range of energies and types of translational and rotational excitations of the guest molecule which go substantially beyond those considered in the earlier theoretical investigations of this system, revealing new information about the TR energy level structure. The computed j = 1 and j = 2 rotational levels and their splittings, as well as the translational fundamental, are in semi-quantitative agreement with the available INS and IR data, indicating the need for a more accurate intermolecular PES. Our calculations reveal a strong dependence of the TR energy levels, in particular their splittings, on the setting angle which defines the orientation of the C60 molecules relative to their local threefold axes.

We present a new method for deriving the stellar birth function (SBF) of resolved stellar populations. The SBF (stars born per unit mass, time, and metallicity) is the combination of the initial mass function (IMF), the star-formation history (SFH), and the metallicity distribution function (MDF). The framework of our analysis is that of Poisson Point Processes (PPPs), a class of statistical models suitable when dealing with points (stars) in a multidimensional space (the measurement space of multiple photometric bands). The theory of PPPs easily accommodates the modeling of measurement errors as well as that of incompleteness. Compared to most of the tools used to study resolved stellar populations, our method avoids binning stars in the color-magnitude diagram and uses the entirety of the information (i.e., the whole likelihood function) for each data point; the proper combination of the individual likelihoods allows the computation of the posterior probability for the global population parameters. This includes unknowns such as the IMF slope and combination of SFH and MDF, which are rarely solved for simultaneously in the literature, however entangled and correlated they might be. Our method also allows proper inclusion of nuisance parameters, such as distance and extinction distributions. The aim of this paper, is to assess the validity of this new approach under a range of assumptions, using only simulated data. Forthcoming work will show applications to real data. Although it has a broad scope of possible applications, we have developed this method to study multi-band HST observations of the Milky Way Bulge. Therefore we will focus on simulations with characteristics similar to those of the Galactic Bulge., Comment: Accepted for publication in ApJ

We present the Bayesian Extinction And Stellar Tool (BEAST), a probabilistic approach to modeling the dust extinguished photometric spectral energy distribution of an individual star while accounting for observational uncertainties common to large resolved star surveys. Given a set of photometric measurements and an observational uncertainty model, the BEAST infers the physical properties of the stellar source using stellar evolution and atmosphere models and constrains the line of sight extinction using a newly developed mixture model that encompasses the full range of dust extinction curves seen in the Local Group. The BEAST is specifically formulated for use with large multi-band surveys of resolved stellar populations. Our approach accounts for measurement uncertainties and any covariance between them due to stellar crowding (both systematic biases and uncertainties in the bias) and absolute flux calibration, thereby incorporating the full information content of the measurement. We illustrate the accuracy and precision possible with the BEAST using data from the Panchromatic Hubble Andromeda Treasury. While the BEAST has been developed for this survey, it can be easily applied to similar existing and planned resolved star surveys., 20 pages, 19 figures, ApJ, in press