Search

The relation between magnetic field strength B and gas density n in the interstellar medium is of fundamental importance to many areas of astrophysics, from protostellar disks to galaxy evolution. We present and compare Bayesian analyses of the B - n relation for a comprehensive observational data set, as well as a large body of numerical MHD simulations. We extend the original Zeeman relation of Crutcher et al. (2010) with a large body of magnetic data that includes 700 observations with the Davis-Chandrasekhar-Fermi method. By using a new multiparameter Bayesian analysis we present a new, more general, time-averaged observational relation: B \propto n^{0.27 \pm 0.017} for n \leq n_0 and B \propto n^{0.54 \pm 0.18} for n \geq n_0, with n_0 = 924^(+145-144) cm^-3. We perform a separate analysis on 19 numerical magnetohydrodynamics simulations that cover a wide range of scales, resolutions, initial conditions, and completed with a variety of codes: arepo, flash, pencil, and ramses. The power law exponents derived from the simulations depend on several physical factors including: dynamo effects, time scales, turbulence, and the initial seed field strength. In particular, early-time simulations where the density, velocity and magnetic fields are unevolved do not match the observational scalings. The simulations that trace the observed density range best, the evolved dwarf galaxy and Milky Way like galaxy simulations, settle into a near consistent exponent of = 0.5 in the dense gas, with variability in the diffuse gas exponent., Comment: 23 figures, 28 pages, submitted to MNRAS, Comments welcome

Submillimeter and millimeter wavelengths provide a unique view of the Universe, from the gas and dust that fills and surrounds galaxies to the chromosphere of our own Sun. Current single-dish facilities have presented a tantalising view of the brightest (sub-)mm sources, and interferometers have provided the exquisite resolution necessary to analyse the details in small fields, but there are still many open questions that cannot be answered with current facilities. In this report we summarise the science that is guiding the design of the Atacama Large Aperture Submillimeter Telescope (AtLAST). We demonstrate how tranformational advances in topics including star formation in high redshift galaxies, the diffuse circumgalactic medium, Galactic ecology, cometary compositions and solar flares motivate the need for a 50m, single-dish telescope with a 1-2 degree field of view and a new generation of highly multiplexed continuum and spectral cameras. AtLAST will have the resolution to drastically lower the confusion limit compared to current single-dish facilities, whilst also being able to rapidly map large areas of the sky and detect extended, diffuse structures. Its high sensitivity and large field of view will open up the field of submillimeter transient science by increasing the probability of serendipitous detections. Finally, the science cases listed here motivate the need for a highly flexible operations model capable of short observations of individual targets, large surveys, monitoring programmes, target of opportunity observations and coordinated observations with other observatories. AtLAST aims to be a sustainable, upgradeable, multipurpose facility that will deliver orders of magnitude increases in sensitivity and mapping speeds over current and planned submillimeter observatories., Comment: 47 pages, 12 figures. For further details on AtLAST see https://atlast.uio.no

Spectral energy distribution (SED) fitting is one of most commonly used techniques to study the dust properties in Active Galactic Nuclei (AGN). Works implementing this technique commonly use radiative transfer models that assume a variety of dust properties. Despite the key role of this aspect, limited effort has been put forward to explore the chemical composition, the role of different optical properties, and the grain size distribution of dust, all of which can have a substantial impact on the theoretical radiative transfer calculations. In this work, we explore the role of the dust chemical composition in the AGN dusty torus through SED fitting to \emph{Spitzer}/IRS spectra of a sample of 49 nearby AGN with silicate features in emission. We implement a mineralogy model including the popular astronomical silicates and a set of oxides and amorphous silicates with different grain sizes. We find that best fits use principally porous alumina, periclase, and olivine. In terms of mass fractions, $\sim99\%$ of the dust is composed of dust grains of size $\rm{0.1 \mu m}$, with a $<1\%$ contribution from $\rm{3 \mu m}$ grains. Moreover, the astronomical silicates have a very low occurrence in the best fits, suggesting that they are not the most suited dust species to reproduce the silicate features in our sample., Comment: 16 pages, 13 figures, Accepted by MNRAS

We present point-source photometry from the Spitzer Space Telescope's final survey of the Small Magellanic Cloud (SMC). We mapped 30 square degrees in two epochs in 2017, with the second extending to early 2018 at 3.6 and 4.5 microns using the Infrared Array Camera. This survey duplicates the footprint from the SAGE-SMC program in 2008. Together, these surveys cover a nearly 10 yr temporal baseline in the SMC. We performed aperture photometry on the mosaicked maps produced from the new data. We did not use any prior catalogs as inputs for the extractor in order to be sensitive to any moving objects (e.g., foreground brown dwarfs) and other transient phenomena (e.g., cataclysmic variables or FU Ori-type eruptions). We produced a point-source catalog with high-confidence sources for each epoch as well as combined-epoch catalog. For each epoch and the combined-epoch data, we also produced a more complete archive with lower-confidence sources. All of these data products will be available to the community at the Infrared Science Archive., Comment: 16 pages, 11 figures, 6 tables

We follow our first paper with an analysis of the ensemble of the extensive pre-explosion ground- and space-based infrared observations of the red supergiant (RSG) progenitor candidate for the nearby core-collapse supernova SN 2023ixf in Messier 101, together with optical data prior to explosion obtained with the Hubble Space Telescope (HST). We have confirmed the association of the progenitor candidate with the SN, as well as constrained the metallicity at the SN site, based on SN observations with instruments at Gemini-North. The internal host extinction to the SN has also been confirmed from a high-resolution Keck spectrum. We fit the observed spectral energy distribution (SED) for the star, accounting for its intrinsic variability, with dust radiative-transfer modeling, which assume a silicate-rich dust shell ahead of the underlying stellar photosphere. The star is heavily dust-obscured, likely the dustiest progenitor candidate yet encountered. We found median estimates of the star's effective temperature and luminosity of 2770 K and 9.0e4 L_Sun, with 68% credible intervals of 2340--3150 K and (7.5--10.9)e4 L_sun. The candidate may have a Galactic RSG analog, IRC -10414, with a strikingly similar SED and luminosity. Via comparison with single-star evolutionary models we have constrained the initial mass of the progenitor candidate from 12 M_sun to as high as 14 M_sun. We have had available to us an extraordinary view of the SN 2023ixf progenitor candidate, which should be further followed up in future years with HST and the James Webb Space Telescope., Comment: 40 pages, substantive modifications relative to the previous, although the overall conclusions remain the same; to appear in AAS Journals

Observational evidence points to a red supergiant (RSG) progenitor for SN 2023ixf. The progenitor candidate has been detected in archival images at wavelengths (>0.6 micron) where RSGs typically emit profusely. This object is distinctly variable in the infrared (IR). We characterize the variability using pre-explosion mid-IR (3.6 and 4.5 micron) Spitzer and ground-based near-IR (JHKs) archival data jointly covering 19 yr. The IR light curves exhibit significant variability with RMS amplitudes in the range of 0.2-0.4 mag, increasing with decreasing wavelength. From a robust period analysis of the more densely sampled Spitzer data, we measure a period of 1091+/-71 days. We demonstrate using Gaussian Process modeling that this periodicity is also present in the near-IR light curves, thus indicating a common physical origin, which is likely pulsational instability. We use a period-luminosity relation for RSGs to derive a value of M_K=-11.58+/-0.31 mag. Assuming a late M spectral type, this corresponds to log(L/L_sun)=5.27+/-0.12 at T_eff=3200 K and to log(L/L_sun)=5.37+/-0.12 at T_eff=3500 K. This gives an independent estimate of the progenitor's luminosity, unaffected by uncertainties in extinction and distance. Assuming the progenitor candidate underwent enhanced dust-driven mass-loss during the time of these archival observations, and using an empirical period-luminosity-based mass-loss prescription, we obtain a mass-loss rate of around (2-4)x10^-4 M_sun/yr. Comparing the above luminosity with stellar evolution models, we infer an initial mass for the progenitor candidate of 20+/-4 M_sun, making this one of the most massive progenitors for a Type II SN detected to-date., Comment: 20 pages, 2 tables, accepted to ApJ

We present a near-complete catalog of the metal-rich population of Thermally-Pulsing Asymptotic Giant Branch stars in the northwest quadrant of M31. This metal-rich sample complements the equally complete metal-poor Magellanic Cloud AGB catalogs produced by the SAGE program. Our catalog includes HST wide-band photometry from the Panchromatic Hubble Andromeda Treasury survey, HST medium-band photometry used to chemically classify a subset of the sample, and Spitzer mid- and far-IR photometry that we have used to isolate dust-producing AGB stars. We have detected 346,623 AGB stars; these include 4,802 AGB candidates producing considerable dust, and 1,356 AGB candidates that lie within clusters with measured ages, and in some cases metallicities. Using the Spitzer data and chemical classifications made with the medium-band data, we have identified both carbon- and oxygen-rich AGB candidates producing significant dust. We have applied color--mass-loss relations based on dusty AGB stars from the LMC to estimate the dust injection by AGB stars in the PHAT footprint. Applying our color relations to a subset of the chemically-classified stars producing the bulk of the dust, we find that ~97.8% of the dust is oxygen-rich. Using several scenarios for the dust lifetime, we have estimated the contribution of AGB stars to the global dust budget of M31 to be 0.9-35.5%, which is in line with previous estimates in the Magellanic Clouds. Follow-up observations of the M31 AGB candidates with the JWST will allow us to further constrain stellar and chemical evolutionary models, and the feedback and dust production of metal-rich evolved stars., Comment: Accepted for publication in ApJS. 35 pages, 25 Figures

We present the core mass function (CMF) of the massive star-forming clump G33.92+0.11 using 1.3 mm observations obtained with the Atacama Large Millimeter/submillimeter Array (ALMA). With a resolution of 1000 au, this is one of the highest resolution CMF measurements to date. The CMF is corrected by flux and number incompleteness to obtain a sample that is complete for gas masses $M\gtrsim2.0\ M_\odot$. The resulting CMF is well represented by a power-law function ($dN/d\log M\propto M^\Gamma$), whose slope is determined using two different approaches: $i)$ by least-squares fitting of power-law functions to the flux- and number-corrected CMF, and $ii)$ by comparing the observed CMF to simulated samples with similar incompleteness. We provide a prescription to quantify and correct a flattening bias affecting the slope fits in the first approach, which is caused by small-sample or edge effects when the data is represented by either classical histograms or a kernel density estimate, respectively. The resulting slopes from both approaches are in good agreement each other, with $\Gamma=-1.11_{-0.11}^{+0.12}$ being our adopted value. Although this slope appears to be slightly flatter than the Salpeter slope $\Gamma=-1.35$ for the stellar initial mass function (IMF), we find from Monte Carlo simulations that the CMF in G33.92+0.11 is statistically indistinguishable from the Salpeter representation of the stellar IMF. Our results are consistent with the idea that the form of the IMF is inherited from the CMF, at least at high masses and when the latter is observed at high-enough resolution., Comment: Accepted for publication in ApJ

The thermal emission of dust is one of the most important tracers of the interstellar medium: multi-wavelength photometry in the far-infrared (FIR) and submillimeter (submm) can be fitted with a model, providing estimates of the dust mass. The fit results depend on the assumed value for FIR/submm opacity, which in most models - due to the scarcity, until recently, of experimental measurements - is extrapolated from shorter wavelengths. Lab measurements of dust analogues, however, show that FIR opacities are usually higher than the values used in models and depend on temperature, which suggests that dust mass estimates may be biased. To test the extent of this bias, we create multi-wavelength synthetic photometry for dusty galaxies at different temperatures and redshifts, using experimental results for FIR/submm dust opacity, then we fit the synthetic data using standard dust models. We find that the dust masses recovered by typical models are overestimated by a factor 2 to 20, depending on how the experimental opacities are treated. If the experimental dust samples are accurate analogues of interstellar dust, therefore, current dust masses are overestimated by up to a factor 20. The implications for our understanding of dust, both Galactic and at high redshift, are discussed., Comment: 26 pages, 23 figures, to be published in MNRAS

The nearby ($\sim$500 kpc) metal-poor ([Fe/H] $\approx$ -1.2; $Z$ $\approx$ 30% $Z_{\odot}$) star-forming galaxy NGC 6822 has a metallicity similar to systems at the epoch of peak star formation. Through identification and study of dusty and dust-producing stars, it is therefore a useful laboratory to shed light on the dust life cycle in the early Universe. We present a catalog of sources combining near- and mid-IR photometry from the United Kingdom Infrared Telescope (UKIRT; $J$, $H$, and $K$) and the $Spitzer$ $Space$ $Telescope$ (IRAC 3.6, 4.5, 5.8, and 8.0 $\mu$m and MIPS 24 $\mu$m). This catalog is employed to identify dusty and evolved stars in NGC 6822 utilizing three color-magnitude diagrams (CMDs). With diagnostic CMDs covering a wavelength range spanning the near- and mid-IR, we develop color cuts using kernel density estimate (KDE) techniques to identify dust-producing evolved stars, including red supergiant (RSG) and thermally-pulsing asymptotic giant branch (TP-AGB) star candidates. In total, we report 1,292 RSG candidates, 1,050 oxygen-rich AGB star candidates, and 560 carbon-rich AGB star candidates with high confidence in NGC 6822. Our analysis of the AGB stars suggests a robust population inhabiting the central stellar bar of the galaxy, with a measured global stellar metallicity of [Fe/H] = -1.286 $\pm$ 0.095, consistent with previous studies. In addition, we identify 277 young stellar object (YSO) candidates. The detection of a large number of YSO candidates within a centrally-located, compact cluster reveals the existence of an embedded, high-mass star-formation region that has eluded previous detailed study. Spitzer I appears to be younger and more active than the other prominent star-forming regions in the galaxy., Comment: Accepted for publication in The Astrophysical Journal; 27 pages, 13 figures, 12 tables (master catalog will be available from publisher)

This white paper discusses recent progress in the field of evolved stars, primarily highlighting the contributions of the James Clerk Maxwell Telescope. It discusses the ongoing project, the \emph{Nearby Evolved Stars Survey} (NESS), and the potential to build upon NESS in the next decade. It then outlines a number of science cases which may become feasible with the proposed 850\,$\mu$m camera which is due to become available at the JCMT in late 2022. These include mapping the extended envelopes of evolved stars, including in polarisation, and time-domain monitoring of their variations. The improved sensitivity of the proposed instrument will facilitate statistical studies that put the morphology, polarisation properties and sub-mm variability in perspective with a relatively modest commitment of time that would be impossible with current instrumentation. We also consider the role that could be played by other continuum wavelengths, heterodyne instruments or other facilities in contributing towards these objectives., Comment: 11 pages, 2 figures; Part of the 2019 EAO Submillimetre Futures Paper Series, https://www.eaobservatory.org//jcmt/wp-content/uploads/sites/2/2019/11/EAOSubmmFuturesWPs.pdf

We present a comprehensive study of massive young stellar objects (YSOs) in the metal-poor galaxy NGC 6822 using IRAC and MIPS data obtained from the {\em Spitzer Space Telescope}. We find over 500 new YSO candidates in seven massive star-formation regions; these sources were selected using six colour-magnitude cuts. Via spectral energy distribution fitting to the data with YSO radiative transfer models we refine this list, identifying 105 high-confidence and 88 medium-confidence YSO candidates. For these sources we constrain their evolutionary state and estimate their physical properties. The majority of our YSO candidates are massive protostars with an accreting envelope in the initial stages of formation. We fit the mass distribution of the Stage I YSOs with a Kroupa initial mass function and determine a global star-formation rate of 0.039 $M_{\odot} yr^{-1}$. This is higher than star-formation rate estimates based on integrated UV fluxes. The new YSO candidates are preferentially located in clusters which correspond to seven active high-mass star-formation regions which are strongly correlated with the 8 and 24 $\mu$m emission from PAHs and warm dust. This analysis reveals an embedded high-mass star-formation region, Spitzer I, which hosts the highest number of massive YSO candidates in NGC 6822. The properties of Spitzer I suggest it is younger and more active than the other prominent H\,{\sc ii} and star-formation regions in the galaxy., Comment: 16 pages, 10 figures, 9 tables; accepted for publication in MNRAS

We have conducted a highly sensitive census of the evolved-star population in the metal-poor dwarf galaxy Leo P and detected four asymptotic giant branch (AGB) star candidates. Leo P is one of the best examples of a nearby analog of high-redshift galaxies because of its primitive metal content (2% of the solar value), proximity, and isolated nature, ensuring a less complicated history. Using medium-band optical photometry from the Hubble Space Telescope (HST), we have classified the AGB candidates by their chemical type. We have identified one oxygen-rich source which appears to be dusty in both the HST and Spitzer observations. Its brightness, however, suggests it may be a planetary nebula or post-AGB object. We have also identified three carbon-rich candidates, one of which may be dusty. Follow-up observations are needed to confirm the nature of these sources and to study the composition of any dust that they produce. If dust is confirmed, these stars would likely be among the most metal-poor examples of dust-producing stars known and will provide valuable insight into our understanding of dust formation at high redshift., Comment: Accepted for publication in ApJ

We present the highest resolution single-dish submillimetre observations of the detached shell source U Antliae to date. The observations were obtained at $450~\micron$ and $850~\micron$ with SCUBA-2 instrument on the James Clerk Maxwell Telescope as part of the Nearby Evolved Stars Survey. The emission at $850~\micron$ peaks at $40\arcsec$ with hints of a second peak seen at $\sim 20\arcsec$. The emission can be traced out to a radius of $56\arcsec$ at a $3\sigma$ level. The outer peak observed at $850~\micron$ aligns well with the peak observed at Herschel/PACS wavelengths. With the help of spectral energy distribution fitting and radiative transfer calculations of multiple-shell models for the circumstellar envelope, we explore the various shell structures and the variation of grain sizes along the in the circumstellar envelope. We determine a total shell dust mass of $(2.0 \pm 0.3) \times 10^{-5}$ M$_{\odot}$ and established that the thermal pulse which gave rise to the detached shell occurred 3500 $\pm$ 500 years ago., Comment: Accepted for publication in MNRAS

We present extensive optical photometric and spectroscopic observations, from 4 to 482 days after explosion, of the Type II-plateau (II-P) supernova (SN) 2017eaw in NGC 6946. SN 2017eaw is a normal SN II-P intermediate in properties between, for example, SN 1999em and SN 2012aw and the more luminous SN 2004et, also in NGC 6946. We have determined that the extinction to SN 2017eaw is primarily due to the Galactic foreground and that the SN site metallicity is likely subsolar. We have also independently confirmed a tip-of-the-red-giant-branch (TRGB) distance to NGC 6946 of 7.73+/-0.78 Mpc. The distances to the SN that we have also estimated via both the standardized candle method and expanding photosphere method corroborate the TRGB distance. We confirm the SN progenitor identity in pre-explosion archival Hubble Space Telescope (HST) and Spitzer Space Telescope images, via imaging of the SN through our HST Target of Opportunity program. Detailed modeling of the progenitor's spectral energy distribution indicates that the star was a dusty, luminous red supergiant consistent with an initial mass of ~15 Msuns., Comment: Appearing in ApJ

We present JCMT SCUBA-2 $450\mu$m and $850\mu$m observations of 14 Asymptotic Giant Branch (AGB) stars (9 O--rich, 4 C-rich and 1 S--type) and one Red Supergiant (RSG) in the Solar Neighbourhood. We combine these observations with \emph{Herschel}/PACS observations at $70\mu$m and $160\mu$m and obtain azimuthally-averaged surface-brightness profiles and their PSF subtracted residuals. The extent of the SCUBA-2 850 $\mu$m emission ranges from 0.01 to 0.16 pc with an average of $\sim40\%$ of the total flux being emitted from the extended component. By fitting a modified black-body to the four-point SED at each point along the radial profile we derive the temperature ($T$), spectral index of dust emissivity ($\beta$) and dust column density ($\Sigma$) as a function of radius. For all the sources, the density profile deviates significantly from what is expected for a constant mass-loss rate, showing that all the sources have undergone variations in mass-loss during this evolutionary phase. In combination with results from CO line emission, we determined the dust-to-gas mass ratio for all the sources in our sample. We find that, when sources are grouped according to their chemistry, the resulting average dust-to-gas ratios are consistent with the respective canonical values. However we see a range of values with significant scatter which indicate the importance of including spatial information when deriving these numbers., Comment: 18 pages plus 13 pages of supplementary figures. Accepted for publication on MNRAS on 24/05/2018. Updated 05/07/2018: Fixed secondary axes (pc axes) on Radial Profiles in all figures

The tori around active galactic nuclei (AGN) are potential formation sites for large amounts of dust, and they may help resolve the so-called dust budget crisis at high redshift. We investigate the dust composition in 53 of the 87 Palomar Green (PG) quasars showing the 9.7 micron silicate feature in emission. By simultaneously fitting the mid-infrared spectroscopic features and the underlying continuum, we estimate the mass fraction in various amorphous and crystalline dust species. We find that the dust consists predominantly of alumina and amorphous silicates, with a small fraction in crystalline form. The mean crystallinity is 8 +/- 6%, with more than half of the crystallinities greater than 5%, well above the upper limit determined for the Galaxy. Higher values of crystallinity are found for higher oxide fractions and for more luminous sources., Comment: 24 pages, 9 figures, accepted for publication in Planetary and Space Science

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

Asymptotic giant branch (AGB) stars are significant contributors to the chemical enrichment of the interstellar medium (ISM) of galaxies. It is therefore essential to constrain the AGB contribution to the dust budget in galaxies. Recent estimates of the total dust injection rate to the Large and Small Magellanic Clouds (LMC and SMC; Riebel et al. 2012, Boyer et al. 2012, Srinivasan et al. in prep) have used data from the Spitzer Space Telescope SAGE (Surveying the Agents of Galaxy Evolution; Meixner et al. 2006) and SAGE-SMC (Gordon et al. 2011) surveys. When sorted by dust chemistry, the data allow for a comparison of O-rich and carbonaceous dust-production rates. In the LMC, for instance, the rate of dust production from carbon stars is about two and a half times that from oxygen-rich AGBs. A reliable determination of the fractional contributions of the two types of dust would serve as input to models of chemical evolution. However, the Spitzer IRAC photometric bands do not sufficiently probe the characteristic mid-infrared spectral features that can distinguish O-rich AGBs from carbon stars - namely, the 9.7 $\mu$m silicate feature and the 11.3 $\mu$m silicon carbide feature. With the continuous spectral coverage in the 4-30 $\mu$m range, SPICA has the potential to distinguish these two types of chemistries. In this contribution, synthetic photometry from the model grid of AGB stars, GRAMS (Sargent et al. 2011; Srinivasan et al. 2011) will be used to discuss the science possibilities that SPICA might offer this study., Comment: Accepted for publication in the proceedings of the "SPICA Science Conference 2013, From Exoplanets to Distant Galaxies: SPICA's New Window on the Cool Universe", to be published in the ASP Conference Series

The lifecycle of dust in the interstellar medium (ISM) is heavily influenced by outflows from asymptotic giant branch (AGB) and red supergiant (RSG) stars, a large fraction of which is contributed by a few very dusty sources. We compute the dust input to the Small Magellanic Cloud (SMC) by fitting the multi-epoch mid-infrared spectral energy distributions (SEDs) of AGB/RSG candidates with models from the {\em G}rid of {\em R}SG and {\em A}GB {\em M}odel{\em S} (GRAMS) grid, allowing us to estimate the luminosities and dust-production rates (DPRs) of the entire population. By removing contaminants, we guarantee a high-quality dataset with reliable DPRs and a complete inventory of the dustiest sources. We find a global AGB/RSG dust-injection rate of $(1.3\pm 0.1)\times 10^{-6}$ \msunperyr, in agreement with estimates derived from mid-infrared colours and excess fluxes. As in the LMC, a majority (66\%) of the dust arises from the extreme AGB stars, which comprise only $\approx$7\% of our sample. A handful of far-infrared sources, whose 24 \mic\ fluxes exceed their 8 \mic\ fluxes, dominate the dust input. Their inclusion boosts the global DPR by $\approx$1.5$\times$, making it necessary to determine whether they are AGB stars. Model assumptions, rather than missing data, are the major sources of uncertainty; depending on the choice of dust shell expansion speed and dust optical constants, the global DPR can be up to $\approx$10 times higher. Our results suggest a non-stellar origin for the SMC dust, barring as yet undiscovered evolved stars with very high DPRs., Comment: 28 pages, 22 figures, accepted for publication in MNRAS

Ground-based observations at terahertz (THz) frequencies are a newly explorable area of astronomy for the next ten years. We discuss science cases for a first-generation 10-m class THz telescope, focusing on the Greenland Telescope as an example of such a facility. We propose science cases and provide quantitative estimates for each case. The largest advantage of ground-based THz telescopes is their higher angular resolution (~ 4 arcsec for a 10-m dish), as compared to space or airborne THz telescopes. Thus, high-resolution mapping is an important scientific argument. In particular, we can isolate zones of interest for Galactic and extragalactic star-forming regions. The THz windows are suitable for observations of high-excitation CO lines and [N II] 205 um lines, which are scientifically relevant tracers of star formation and stellar feedback. Those lines are the brightest lines in the THz windows, so that they are suitable for the initiation of ground-based THz observations. THz polarization of star-forming regions can also be explored since it traces the dust population contributing to the THz spectral peak. For survey-type observations, we focus on ``sub-THz'' extragalactic surveys, whose uniqueness is to detect galaxies at redshifts z ~ 1--2, where the dust emission per comoving volume is the largest in the history of the Universe. Finally we explore possibilities of flexible time scheduling, which enables us to monitor active galactic nuclei, and to target gamma-ray burst afterglows. For these objects, THz and submillimeter wavelength ranges have not yet been explored., Comment: 39 pages, 17 figures, accepted for publication in PASJ (as a review paper)

We reconsider the origin and processing of dust in elliptical galaxies. We theoretically formulate the evolution of grain size distribution, taking into account dust supply from asymptotic giant branch (AGB) stars and dust destruction by sputtering in the hot interstellar medium (ISM), whose temperature evolution is treated by including two cooling paths: gas emission and dust emission (i.e. gas cooling and dust cooling). With our new full treatment of grain size distribution, we confirm that dust destruction by sputtering is too efficient to explain the observed dust abundance even if AGB stars continue to supply dust grains, and that, except for the case where the initial dust-to-gas ratio in the hot gas is as high as $\sim 0.01$, dust cooling is negligible compared with gas cooling. However, we show that, contrary to previous expectations, cooling does not help to protect the dust; rather, the sputtering efficiency is raised by the gas compression as a result of cooling. We additionally consider grain growth after the gas cools down. Dust growth by the accretion of gas-phase metals in the cold medium increase the dust-to-gas ratio up to $\sim 10^{-3}$ if this process lasts >10/(n_H/10^3 cm^{-3}) Myr, where $n_\mathrm{H}$ is the number density of hydrogen nuclei. We show that the accretion of gas-phase metals is a viable mechanism of increasing the dust abundance in elliptical galaxies to a level consistent with observations, and that the steepness of observed extinction curves is better explained with grain growth by accretion., Comment: 15 pages, 11 figures, accepted for publication in MNRAS

We have identified a new class of Asymptotic Giant Branch (AGB) stars in the Small and Large Magellanic Clouds (SMC/LMC) using optical to infrared photometry, light curves, and optical spectroscopy. The strong dust production and long-period pulsations of these stars indicate that they are at the very end of their AGB evolution. Period-mass-radius relations for the fundamental-mode pulsators give median current stellar masses of 1.14 M_sun in the LMC and 0.94 M_sun in the SMC (with dispersions of 0.21 and 0.18 M_sun, respectively), and models suggest initial masses of <1.5 M_sun and <1.25 M_sun, respectively. This new class of stars includes both O-rich and C-rich chemistries, placing the limit where dredge-up allows carbon star production below these masses. A high fraction of the brightest among them should show S star characteristics indicative of atmospheric C/O ~ 1, and many will form O-rich dust prior to their C-rich phase. These stars can be separated from their less-evolved counterparts by their characteristically red J-[8] colors., Comment: 16 pages, 18 figures, accepted for publication in ApJ

We present the first detection of 24 {\mu}m variability in 24 sources in the Local Group galaxy M33. These results are based on 4 epochs of MIPS observations, which are irregularly spaced over ~750 days. We find that these sources are constrained exclusively to the Holmberg radius of the galaxy, which increases their chances of being members of M33. We have constructed spectral energy distributions (SEDs) ranging from the optical to the sub-mm to investigate the nature of these objects. We find that 23 of our objects are most likely heavily self-obscured, evolved stars; while the remaining source is the Giant HII region, NGC 604. We believe that the observed variability is the intrinsic variability of the central star reprocessed through their circumstellar dust shells. Radiative transfer modeling was carried out to determine their likely chemical composition, luminosity, and dust production rate (DPR). As a sample, our modeling has determined an average luminosity of (3.8 $\pm$ 0.9) x 10$^4$ L$_\odot$ and a total DPR of (2.3 $\pm$ 0.1) x 10$^{-5}$ M$_\odot$ yr$^{-1}$. Most of the sources, given the high DPRs and short wavelength obscuration, are likely "extreme" AGB (XAGB) stars. Five of the sources are found to have luminosities above the classical AGB limit (M$_{\rm bol}$ < -7.1, L > 54,000 L$_\odot$), which classifies them as probably red supergiants (RSGs). Almost all of the sources are classified as oxygen rich. As also seen in the LMC, a significant fraction of the dust in M33 is produced by a handful of XAGB and RSG stars., Comment: 36 pages, 14 figures, 4 tables, Accepted for publication in AJ

We present the detection and analysis of molecular hydrogen emission toward ten interstellar regions in the Large Magellanic Cloud. We examined low-resolution infrared spectral maps of twelve regions obtained with the Spitzer infrared spectrograph (IRS). The pure rotational 0--0 transitions of H$_2$ at 28.2 and 17.1${\,\rm \mu m}$ are detected in the IRS spectra for ten regions. The higher level transitions are mostly upper limit measurements except for three regions, where a 3$\sigma$ detection threshold is achieved for lines at 12.2 and 8.6${\,\rm \mu m}$. The excitation diagrams of the detected H$_2$ transitions are used to determine the warm H$_2$ gas column density and temperature. The single-temperature fits through the lower transition lines give temperatures in the range $86-137\,{\rm K}$. The bulk of the excited H$_2$ gas is found at these temperatures and contributes $\sim$5-17% to the total gas mass. We find a tight correlation of the H$_2$ surface brightness with polycyclic aromatic hydrocarbon and total infrared emission, which is a clear indication of photo-electric heating in photodissociation regions. We find the excitation of H$_2$ by this process is equally efficient in both atomic and molecular dominated regions. We also present the correlation of the warm H$_2$ physical conditions with dust properties. The warm H$_2$ mass fraction and excitation temperature show positive correlations with the average starlight intensity, again supporting H$_2$ excitation in photodissociation regions., Comment: Accepted for publication in MNRAS

We show that Montecarlo simulations of the TP-AGB stellar population in the LMC and SMC galaxies using the CB* models produce LF and color distributions that are in closer agreement with observations than those obtained with the BC03 and CB07 models. This is a progress report of work that will be published elsewhere., Comment: 5 pages, 4 figures, in proceedings of IAU Symposium 295, "The intriguing life of massive galaxies", eds. D. Thomas, A. Pasquali and I. Ferreras

We study the emission by dust and stars in the Large and Small Magellanic Clouds, a pair of low-metallicity nearby galaxies, as traced by their spatially resolved spectral energy distributions (SEDs). This project combines Herschel Space Observatory PACS and SPIRE far-infrared photometry with other data at infrared and optical wavelengths. We build maps of dust and stellar luminosity and mass of both Magellanic Clouds, and analyze the spatial distribution of dust/stellar luminosity and mass ratios. These ratios vary considerably throughout the galaxies, generally between the range $0.01\leq L_{\rm dust}/L_\ast\leq 0.6$ and $10^{-4}\leq M_{\rm dust}/M_\ast\leq 4\times10^{-3}$. We observe that the dust/stellar ratios depend on the interstellar medium (ISM) environment, such as the distance from currently or previously star-forming regions, and on the intensity of the interstellar radiation field (ISRF). In addition, we construct star formation rate (SFR) maps, and find that the SFR is correlated with the dust/stellar luminosity and dust temperature in both galaxies, demonstrating the relation between star formation, dust emission and heating, though these correlations exhibit substantial scatter., Comment: 15 pages, 18 figures; ApJ, in press; version published in the journal will have higher-resolution figures

The Small Magellanic Cloud (SMC) provides a unique laboratory for the study of the lifecycle of dust given its low metallicity (~1/5 solar) and relative proximity (~60 kpc). This motivated the SAGE-SMC (Surveying the Agents of Galaxy Evolution in the Tidally-Stripped, Low Metallicity Small Magellanic Cloud) Spitzer Legacy program with the specific goals of studying the amount and type of dust in the present interstellar medium, the sources of dust in the winds of evolved stars, and how much dust is consumed in star formation. This program mapped the full SMC (30 sq. deg.) including the Body, Wing, and Tail in 7 bands from 3.6 to 160 micron using the IRAC and MIPS instruments on the Spitzer Space Telescope. The data were reduced, mosaicked, and the point sources measured using customized routines specific for large surveys. We have made the resulting mosaics and point source catalogs available to the community. The infrared colors of the SMC are compared to those of other nearby galaxies and the 8 micron/24 micron ratio is somewhat lower and the 70 micron/160 micron ratio is somewhat higher than the average. The global infrared spectral energy distribution shows that the SMC has ~3X lower aromatic emission/PAH (polycyclic aromatic hydrocarbon) abundances compared to most nearby galaxies. Infrared color-magnitude diagrams are given illustrating the distribution of different asymptotic giant branch stars and the locations of young stellar objects. Finally, the average spectral energy distribution (SED) of HII/star formation regions is compared to the equivalent Large Magellanic Cloud average HII/star formation region SED. These preliminary results are expanded in detail in companion papers., Comment: 16 pages, 11 figures, AJ, in press

We investigate the infrared (IR) properties of cool, evolved stars in the Small Magellanic Cloud (SMC), including the red giant branch (RGB) stars and the dust-producing red supergiant (RSG) and asymptotic giant branch (AGB) stars using observations from the Spitzer Space Telescope Legacy program entitled: "Surveying the Agents of Galaxy Evolution in the Tidally-stripped, Low Metallicity SMC", or SAGE-SMC. The survey includes, for the first time, full spatial coverage of the SMC bar, wing, and tail regions at infrared (IR) wavelengths (3.6 - 160 microns). We identify evolved stars using a combination of near-IR and mid-IR photometry and point out a new feature in the mid-IR color-magnitude diagram that may be due to particularly dusty O-rich AGB stars. We find that the RSG and AGB stars each contribute ~20% of the global SMC flux (extended + point-source) at 3.6 microns, which emphasizes the importance of both stellar types to the integrated flux of distant metal-poor galaxies. The equivalent SAGE survey of the higher-metallicity Large Magellanic Cloud (SAGE-LMC) allows us to explore the influence of metallicity on dust production. We find that the SMC RSG stars are less likely to produce a large amount of dust (as indicated by the [3.6]-[8] color). There is a higher fraction of carbon-rich stars in the SMC, and these stars appear to able to reach colors as red as their LMC counterparts, indicating that C-rich dust forms efficiently in both galaxies. A preliminary estimate of the dust production in AGB and RSG stars reveals that the extreme C-rich AGB stars dominate the dust input in both galaxies, and that the O-rich stars may play a larger role in the LMC than in the SMC., Comment: Accepted for publication in AJ. 25 pages, 36 figures, Table 4 will be available electronically from AJ

The total dust return rate from AGB and RSG star outflows is an important constraint to galactic chemical evolution models. However, this requires detailed radiative transfer (RT) modeling of individual stars, which becomes impractical for large data sets. Another approach is to select the best-fit spectral energy distribution (SED) from a grid of dust shell models, allowing for a faster determination of the luminosities and mass-loss rates for entire samples. We have developed the Grid of RSG and AGB ModelS (GRAMS) to measure the mass-loss return from evolved stars. The models span the range of stellar, dust shell and grain properties relevant to evolved stars. In this paper we present the carbon-star grid and compare our results with data of Large Magellanic Cloud (LMC) carbon stars from the SAGE and SAGE-Spec programs. We generate spherically symmetric dust shell models using the 2Dust code, with hydrostatic models for the central stars. We explore five values of the inner radius R_in of the dust shell (1.5, 3, 4.5, 7 and 12 R_star). We use amorphous carbon dust mixed with 10% silicon carbide by mass. The grain sizes follows a KMH distribution. The models span 26 values of 11.3 um optical depth, ranging from 0.001 to 4. For each model, 2Dust calculates the output SED from 0.2 to 200 um. Over 12,000 models have dust temperatures below 1800 K. The GRAMS synthetic photometry is in good agreement with SAGE photometry for LMC carbon-rich and extreme AGB star candidates, as well as spectroscopically confirmed carbon stars from the SAGE-Spec study. Our models reproduce the IRAC colors of most of the extreme AGB star candidates, consistent with the expectation that a majority of these enshrouded stars have carbon-rich dust. Finally, we fit the SEDs of some well-studied carbon stars and compare the resulting luminosities and mass-loss rates with those from previous studies., Comment: 16 pages, 13 figures, accepted for publication in Astronomy & Astrophysics

We present a 2Dust model for the dust shell around a LMC long-period variable (LPV) previously studied as part of the OGLE survey. OGLE LMC LPV 28579 (SAGE J051306.40-690946.3) is a carbon-rich asymptotic giant branch (AGB) star for which we have photometry and spectra from the Spitzer SAGE and SAGE-Spec programs along with UBVIJHK_s photometry. By modeling this source, we obtain a baseline set of dust properties to be used in the construction of a grid of models for carbon stars. We reproduce its spectral energy distribution using a mixture of AmC and SiC (15% by mass). The grain sizes are distributed according to the KMH model. The best-fit model has an optical depth of 0.28 for the shell at the peak of the SiC feature, with R_in~1430 R_sun or 4.4 R_star. The temperature at this inner radius is 1310 K. Assuming an expansion velocity of 10 km s^-1, we obtain a dust mass-loss rate of 2.5x10^-9 M_sun yr-1. We calculate a 15% variation in this rate by testing the fit sensitivity against variation in input parameters. We also present a simple model for the molecular gas in the extended atmosphere that could give rise to the 13.7 \mu m feature seen in the spectrum. We find that a combination of CO and C_2H_2 gas at an excitation temperature of about 1000 K and column densities of 3x10^21 cm^-2 and 10^19 cm^-2 respectively are able to reproduce the observations. Given that the excitation temperature is close to T_dust(R_in), most of the molecular contribution probably arises from the inner shell region. The luminosity corresponding to the first epoch of SAGE observations is 6580 L_sun. For an effective temperature of about 3000 K, this implies a stellar mass of 1.5-2 M_sun and an age of 1-2.5 Gyr. For a gas:dust ratio of 200, we obtain a gas mass-loss rate of 5.0x10^-7 M_sun yr^-1, consistent with the gas mass-loss rates estimated from the period, color and 8 \mu m flux of the source., Comment: 14 pages, 5 figures, accepted for publication in Astronomy & Astrophysics

We combine variability information from the MAssive Compact Halo Objects (MACHO) survey of the Large Magellanic Cloud (LMC) with infrared photometry from the Spitzer Space Telescope Surveying the Agents of a Galaxy's Evolution (SAGE) survey to create a dataset of ~30 000 variable red sources. We photometrically classify these sources as being on the first ascent of the Red Giant Branch (RGB), or as being in one of three stages along the Asymptotic Giant Branch (AGB): oxygen-rich, carbon-rich, or highly reddened with indeterminate chemistry ("extreme" AGB candidates). We present linear period-luminosity relationships for these sources using 8 separate infrared bands (J, H, K, 3.6, 4.5, 5.8, 8.0, and 24 micron) as proxies for the luminosity. We find that the wavelength dependence of the slope of the period-luminosity relationship is different for different photometrically determined classes of AGB stars. Stars photometrically classified as O-rich show the least variation of slope with wavelength, while dust enshrouded extreme AGB stars show a pronounced trend toward steeper slopes with increasing wavelength. We find that O-rich AGB stars pulsating in the fundamental mode obey a period-magnitude relation with a slope of -3.41 +/- 0.04 when magnitude is measured in the 3.6 micron band, in contrast to C-rich AGB stars, which obey a relation of slope -3.77 +/- 0.05.

The asymptotic giant branch (AGB) phase is the penultimate stage of evolution for low- and intermediate-mass stars. AGB star outflows inject a significant amount of material into the interstellar medium (ISM), seeding new star formation. AGB mass loss is thus a crucial component of galactic chemical evolution. The Large Magellanic Cloud (LMC) is an excellent site for AGB studies. Over 40,000 AGB candidates have been identified using photometric data from the Spitzer Space Telescope Surveying The Agents of a Galaxy's Evolution (SAGE) mid-infrared (MIR) survey, including about 35,000 oxygen-rich, 7000 carbon-rich and 1400 "extreme" sources. For the first time, SAGE photometry reveals two distinct populations of O-rich sources in the LMC: a faint population that gradually evolves into C-rich stars and a bright, massive population that circumvents this evolution, remaining O-rich. This work aims to quantify the mass-loss return from AGB stars to the LMC, a rough estimate for which is derived from the amount of MIR dust emission in excess of that from starlight. I show that this excess flux is a good proxy for the mass-loss rate, and I calculate the total AGB injection rate to be (5.9-13) x 10{-3} Msun/yr. A more accurate determination requires detailed dust radiative transfer (RT) modeling. For this purpose, I present a grid of C-rich AGB models generated by the RT code 2DUST, spanning a range of effective temperatures, gravities, dust shell radii and optical depths as well as a baseline set of dust properties obtained by modeling a carbon star, data for which was acquired as part of the spectroscopic follow-up to SAGE. AGB stars are the best laboratories for dust studies, and the development of a model grid will reinforce future research in this field., Comment: Ph.D. Thesis, 101 pages, 7 tables, 31 figures

We present empirical relations describing excess emission from evolved stars in the Large Magellanic Cloud (LMC) using data from the SAGE (Surveying the Agents of a Galaxy's Evolution) survey which includes the IRAC 3.6, 4.5, 5.8 and 8.0 \mu m and MIPS 24, 70 and 160 \mu m bands. We combine the SAGE data with the Two Micron All Sky Survey (2MASS; J, H and Ks) and the optical Magellanic Cloud Photometric Survey (MCPS; U, B, V and I) point source catalogs to create complete spectral energy distributions (SEDs) of the asymptotic giant branch (AGB) star candidates in the LMC. AGB star outflows are among the main producers of dust in a galaxy, and this mass loss results in an excess in the fluxes observed in the 8 and 24 \mic m bands. We identify oxygen-rich, carbon-rich and extreme AGB star populations in our sample based on their 2MASS and IRAC colors. We calculate excesses fluxes in the mid-IR bands by comparison of the SEDs with model photospheres. We find about 16,000 O-rich, 6300 C-rich and 1000 extreme sources with reliable 8 \mu m excesses, and about 4500 O-rich, 5300 C-rich and 960 extreme sources with reliable 24 \mic m excesses. The excesses are in the range 0.1 mJy--5 Jy. The 8 and 24 \mic m excesses for all three types of AGB candidates show a general increasing trend with luminosity. The dust color temperature derived from the ratio of the 8 and 24 \mic m excesses decreases with an increase in excess, while the 24 \mic m optical depth increases with excess. The extreme AGB candidates are the major contributors to the mass loss, and we estimate the total AGB mass-loss return to the LMC to be (5.9--13)\times 10^{-3} Msun yr^{-1}., Comment: 41 pages, 18 figures, accepted for publication in the Astronomical Journal