Source author record

J. M. Carpenter

J. M. Carpenter appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

28works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

28 published item(s)

preprint2026arXiv

The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) I: Motivation, sample, data reduction, and results overview

The outer regions of planetary systems host dusty debris discs analogous to the Kuiper belt (exoKuiper belts), which provide crucial constraints on planet formation and evolution processes. ALMA dust observations have revealed a great diversity, and that some belts contain CO gas, whose origin and implications are uncertain. Most of this progress, however, has been limited by low-resolution observations. We conducted the first ALMA large programme dedicated to debris discs: the ALMA survey to Resolve exoKuiper belt Substructures (ARKS). We selected the 24 most promising belts to constrain their detailed radial and vertical structure, and to characterise the gas content. We constrained the radial and vertical distribution of dust, as well as the presence of asymmetries. For a subset of six belts with CO gas, we constrained the gas distribution and kinematics. To interpret these observations, we used a wide range of dynamical models. The first ARKS results are presented as a series of ten papers. We discovered that up to 33% of our sample exhibits multiple dusty rings. For highly inclined belts, we found that non-Gaussian vertical distributions are common and are indicative of multiple dynamical populations. We also found that 10 of the 24 belts present asymmetries. We find that the CO gas is radially broader than the dust, but this could be an effect of optical depth. At least one system shows non-Keplerian kinematics due to strong pressure gradients, which may have triggered a vortex that trapped dust in an arc. Finally, we find evidence that the micron-sized grains may be affected by gas drag in gas rich systems. ARKS has revealed a great diversity of structures in exoKuiper belts that may arise when they are formed in protoplanetary discs or subsequently via interactions with planets and/or gas. We encourage the community to explore the reduced data and data products.

preprint2026arXiv

The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) IV: CO gas imaging and overview

CO gas is detected in a significant number of debris discs, but its origin and evolution remains unclear. Key constraints are its mass and spectro-spatial distribution, which are coupled through optical depth and have only been analysed at low to moderate resolution so far. The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) is the first ALMA large program to target debris discs at high spectro-spatial resolution. We used $^{12}$CO and $^{13}$CO J=3-2 line data of 18 ARKS debris belts, 5 of which were already known to host gas, to analyse the spectro-spatial distribution of CO, constrain the gas masses, and to search for gas in the remaining systems. We developed a line-imaging pipeline and produced line cubes for each disc, with a spatial resolution down to $\sim$70 mas and spectral resolution of 26 m s$^{-1}$. Using spectro-spatial shifting and stacking, we produced high signal-to-noise maps, and radial and spectral profiles that reveal the distribution and kinematics of gas in 5 gas-bearing discs. For these discs, we constrained the inner radius of the $^{12}$CO, and found the radial brightness profile of CO peaked interior to the dust ring, but that CO was more radially extended than the dust. We present the first radially resolved $^{12}$CO/$^{13}$CO isotopologue flux ratios in gas-bearing debris discs, which are constant with radius for the majority of systems, indicating $^{12}$CO and $^{13}$CO are both optically thick or thin throughout the discs. We report CO line fluxes/upper limits for all systems and optical depth dependant masses for the 5 gas-bearing systems. Finally, we analysed the $^{12}$CO J=3-2 line luminosities for the ARKS debris discs and discs from the literature. We confirm that gas is mostly detected in young systems. However, the high scatter seen in young/high fractional luminosity systems indicates no trend within the systems with detected gas.

preprint2026arXiv

The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) V: Comparison between scattered light and thermal emission

Debris discs are analogues to our own Kuiper belt around main-sequence stars and are therefore referred to as exoKuiper belts. They have been resolved at high angular resolution at wavelengths spanning the optical to the submillimetre-millimetre regime. Short wavelengths probe the light scattered by such discs, which is dominated by micron-sized dust particles, while millimetre wavelengths probe the thermal emission of millimetre-sized particles. Determining differences in the dust distribution between millimetre- and micron-sized dust is fundamental to revealing the dynamical processes affecting the dust in debris discs. We aim to compare the scattered light from the discs of the ALMA survey to Resolve exoKuiper belt Substructures (ARKS) with the thermal emission probed by ALMA. We focus on the radial distribution of the dust. We used high-contrast scattered light observations obtained with VLT/SPHERE, GPI, and the HST to uniformly study the dust distribution in those systems and compare it to the dust distribution extracted from the ALMA observations carried out in the course of the ARKS project. We also set constraints on the presence of planets by using these high-contrast images combined with exoplanet evolutionary models. 15 of the 24 discs comprising the ARKS sample are detected in scattered light, with TYC9340-437-1 being imaged for the first time at near-infrared wavelengths. For 6 of those 15 discs, the dust surface density seen in scattered light peaks farther out compared to that observed with ALMA. These 6 discs except one are known to also host cold CO gas. Conversely, the systems without significant offsets are not known to host gas, except one. This observational study suggests that the presence of gas in debris discs may affect the small and large grains differently, pushing the small dust to greater distances where the gas is less abundant.

preprint2026arXiv

The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) VII: Optically thick gas with broad CO gaussian local line profiles in the HD 121617 disc

CO gas has been detected in $\sim$20 debris discs. We present ALMA observations of the CO-rich HD 121617 debris disc from the ARKS survey. Using high-resolution Band 7 observations of $^{12}CO \ J=3-2$, we analyse local CO line profiles to investigate optical depth, CO mass, and temperature. Spectra are aligned and stacked in concentric annuli to produce local line profiles. The resulting profiles are Gaussian-shaped and broadened by Keplerian shear. The line profiles are modelled using both a simplified toy model and a RADMC-3D model including projection effects and Keplerian shear. Fitting the RADMC-3D model to the $^{13}$CO data, we find that an optically thick model with a temperature of 38 K and a CO mass of $2 \times 10^{-3}$ M$_{\oplus}$ reproduces the observations. The model reproduces the enhanced emission at orbital azimuths of $\sim \pm45^{\circ}$ and $\pm135^{\circ}$, forming an X-shaped structure in the velocity-integrated intensity map, as well as the broader $^{12}$CO linewidth relative to $^{13}$CO. Scaling the model by the ISM abundance ratio ($\sim$77) also reproduces the $^{12}$CO data, though high optical depths and model assumptions limit mass constraints. We find that azimuthally averaged local line profiles appear Gaussian regardless of optical depth, cautioning against their use for distinguishing optically thin and thick emission. We constrain the mean molecular weight to $12.6_{-1.1}^{+1.3}$, dependent on model assumptions. Our $^{13}$CO results suggest that C$^{18}$O may also be optically thick in CO-rich debris discs and that the mean molecular weight is significantly higher than if H$_2$ were the dominant gas species, suggesting a non-primordial composition.

preprint2026arXiv

The ALMA survey to Resolve exoKuiper belt Substructures (ARKS): VI. Asymmetries and offsets

Asymmetries in debris discs provide unique clues to understand the evolution and architecture of planetary systems.** The aim of the ALMA survey to Resolve exoKuiper belt Substructures (ARKS) is to expand our understanding of radial and vertical dust structures, as well as gas distributions and kinematics, in debris discs.** Here, in ARKS~VI, we present a systematic analysis of the asymmetries and stellocentric offsets present in the ALMA continuum data for the ARKS survey. Our aims are to identify asymmetries in debris disc dust distributions, quantify debris disc asymmetry properties, and discuss the potential origins of debris disc asymmetries.** We utilised empirical methods to identify emission asymmetries** and the presence of offset emission between disc centres and the locations of the host stars, via an analysis of their calibration procedures and disc properties. We associated observational asymmetry types** and plausible physical classes** associated with each source. We show that there are ten systems, almost half of the ARKS sample, that host either a continuum emission asymmetry or offset emission. Three systems host offsets (HD15115, HD32297, and HD109573 (HR4796)), four host azimuthal asymmetries (HD9672 (49Ceti), HD92945, HD107146, and HD121617), two host an asymmetry in their major axis (HD10647 (q$^1$ Eri), and HD39060 ($β$ Pic)), and one hosts an asymmetry in their minor axis (HD61005). We attribute the offset asymmetries to non-zero eccentricities, and three of the azimuthal asymmetries to arcs. The presence of an asymmetry or offset in the ARKS sample appears to be correlated with the fractional luminosity of cold dust.** Conclusions: This study demonstrates that debris disc asymmetries in the ARKS sample are common, and plausibly so in the wider population of debris discs at (sub)-millimetre wavelengths.** ** = ABRIDGED FOR ARXIV: FULL ABSTRACT IN PAPER

preprint2026arXiv

The ALMA survey to Resolve exoKuiper belt Substructures (ARKS). X. Interpreting the peculiar dust rings around HD 131835

Dusty discs detected around main-sequence stars are thought to be signs of planetesimal belts in which the dust distribution is shaped by collisional and dynamical processes, including interactions with gas if present. The debris disc around the young A-type star HD 131835 is composed of two dust rings at ~65 au and ~100 au, a third unconstrained innermost component, and a gaseous component centred at ~65 au. New ALMA observations show that the inner of the two dust rings is brighter than the outer one, in contrast with previous observations in scattered light. We explore two scenarios that could explain these observations: the two dust rings might represent distinct planetesimal belts with different collisional properties, or only the inner ring might contain planetesimals while the outer ring consists entirely of dust that has migrated outwards due to gas drag. To explore the first scenario, we employed a state-of-the-art collisional evolution code. To test the second scenario, we used a simple dynamical model of dust grain evolution in an optically thin gaseous disc. Collisional models of two planetesimal belts cannot fully reproduce the observations by only varying their dynamical excitation, and matching the data through a different material strength requires an extreme difference in dust composition. The gas-driven scenario can reproduce the location of the outer ring and the brightness ratio of the two rings from scattered light observations, but the resulting outer ring is too faint overall in both scattered light and sub-millimetre emission. The dust rings in HD 131835 could be produced from two planetesimal belts, although how these belts would attain the required extremely different properties needs to be explained. The dust-gas interaction is a plausible alternative explanation and deserves further study using a more comprehensive model.

preprint2022arXiv

Distribution of solids in the rings of the HD 163296 disk: a multiwavelength study

In this paper we analyze new observations from ALMA and VLA, at a high angular resolution corresponding to 5 - 8 au, of the protoplanetary disk around HD 163296 to determine the dust spatial distribution and grain properties. We fit the spectral energy distribution as a function of the radius at five wavelengths from 0.9 to 9\,mm, using a simple power law and a physical model based on an analytic description of radiative transfer that includes isothermal scattering. We considered eight dust populations and compared the models' performance using Bayesian evidence. Our analysis shows that the moderately high optical depth ($τ$>1) at $λ\leq$ 1.3 mm in the dust rings artificially lower the millimeter spectral index, which should therefore not be considered as a reliable direct proxy of the dust properties and especially the grain size. We find that the outer disk is composed of small grains on the order of 200 $μ$m with no significant difference between rings at 66 and 100 au and the adjacent gaps, while in the innermost 30 au, larger grains ($\geq$mm) could be present. We show that the assumptions on the dust composition have a strong impact on the derived surface densities and grain size. In particular, increasing the porosity of the grains to 80\% results in a total dust mass about five times higher with respect to grains with 25\% porosity. Finally, we find that the derived opacities as a function of frequency deviate from a simple power law and that grains with a lower porosity seem to better reproduce the observations of HD163296. While we do not find evidence of differential trapping in the rings of HD163296, our overall results are consistent with the postulated presence of giant planets affecting the dust temperature structure and surface density, and possibly originating a second-generation dust population of small grains.

preprint2015arXiv

A multi-wavelength analysis for interferometric (sub-)mm observations of protoplanetary disks: radial constraints on the dust properties and the disk structure

Theoretical models of grain growth predict dust properties to change as a function of protoplanetary disk radius, mass, age and other physical conditions. We lay down the methodology for a multi-wavelength analysis of (sub-)mm and cm continuum interferometric observations to constrain self-consistently the disk structure and the radial variation of the dust properties. The computational architecture is massively parallel and highly modular. The analysis is based on the simultaneous fit in the uv-plane of observations at several wavelengths with a model for the disk thermal emission and for the dust opacity. The observed flux density at the different wavelengths is fitted by posing constraints on the disk structure and on the radial variation of the grain size distribution. We apply the analysis to observations of three protoplanetary disks (AS 209, FT Tau, DR Tau) for which a combination of spatially resolved observations in the range ~0.88mm to ~10mm is available (from SMA, CARMA, and VLA), finding evidence of a decreasing maximum dust grain size (a_max) with radius. We derive large a_max values up to 1 cm in the inner disk between 15 and 30 AU and smaller grains with a_max~1 mm in the outer disk (R > 80AU). In this paper we develop a multi-wavelength analysis that will allow this missing quantity to be constrained for statistically relevant samples of disks and to investigate possible correlations with disk or stellar parameters.

preprint2015arXiv

An ATCA survey of debris disks at 7 millimeters

We present ATCA continuum observations at a wavelength of 6.8 mm of five debris disks: $β$ Pictoris, q$^1$ Eridani, HD 107146, HD 181327, and HD 95086. These observations provide the detection at the longest wavelengths obtained to date for all these debris disks. By combining our 6.8 mm data with previous detections at shorter sub-millimeter/millimeter wavelengths we measure the long wavelength spectral index of these sources. We then use previous estimates for the temperature of the emitting dust to derive the spectral index of the dust emissivity. Under the assumption that all the detected flux comes from dust only, we constrain the slope of the solid size distribution, assumed to be a power-law. The values that we infer for the slope of the size distribution range between about 3.36 and 3.50. We compare our findings with the case of the Fomalhaut debris disk and use these results to test the predictions of collisional cascades of planetesimal belts.

preprint2014arXiv

ALMA observations of the debris disk around the young Solar Analog HD 107146

We present ALMA continuum observations at a wavelength of 1.25 mm of the debris disk surrounding the $\sim$ 100 Myr old solar analog HD 107146. The continuum emission extends from about 30 to 150 AU from the central star with a decrease in the surface brightness at intermediate radii. We analyze the ALMA interferometric visibilities using debris disk models with radial profiles for the dust surface density parametrized as i) a single power-law, ii) a single power-law with a gap, and iii) a double power-law. We find that models with a gap of radial width $\sim 8$ AU at a distance of $\sim 80$ AU from the central star, as well as double power-law models with a dip in the dust surface density at $\sim 70$ AU provide significantly better fits to the ALMA data than single power-law models. We discuss possible scenarios for the origin of the HD 107146 debris disk using models of planetesimal belts in which the formation of Pluto-sized objects trigger disruptive collisions of large bodies, as well as models which consider the interaction of a planetary system with a planetesimal belt and spatial variation of the dust opacity across the disk. If future observations with higher angular resolution and sensitivity confirm the fully-depleted gap structure discussed here, a planet with a mass of approximately a few Earth masses in a nearly circular orbit at $\sim 80$ AU from the central star would be a possible explanation for the presence of the gap.

preprint2014arXiv

On the structure of the transition disk around TW Hya

For over a decade, the structure of the inner cavity in the transition disk of TW Hydrae has been a subject of debate. Modeling the disk with data obtained at different wavelengths has led to a variety of proposed disk structures. Rather than being inconsistent, the individual models might point to the different faces of physical processes going on in disks, such as dust growth and planet formation. Our aim is to investigate the structure of the transition disk again and to find to what extent we can reconcile apparent model differences. A large set of high-angular-resolution data was collected from near-infrared to centimeter wavelengths. We investigated the existing disk models and established a new self-consistent radiative-transfer model. A genetic fitting algorithm was used to automatize the parameter fitting. Simple disk models with a vertical inner rim and a radially homogeneous dust composition from small to large grains cannot reproduce the combined data set. Two modifications are applied to this simple disk model: (1) the inner rim is smoothed by exponentially decreasing the surface density in the inner ~3 AU, and (2) the largest grains (>100 um) are concentrated towards the inner disk region. Both properties can be linked to fundamental processes that determine the evolution of protoplanetary disks: the shaping by a possible companion and the different regimes of dust-grain growth, respectively. The full interferometric data set from near-infrared to centimeter wavelengths requires a revision of existing models for the TW Hya disk. We present a new model that incorporates the characteristic structures of previous models but deviates in two key aspects: it does not have a sharp edge at 4 AU, and the surface density of large grains differs from that of smaller grains. This is the first successful radiative-transfer-based model for a full set of interferometric data.

preprint2014arXiv

Young Stellar Object Variability (YSOVAR): Long Timescale Variations in the Mid-Infrared

The YSOVAR (Young Stellar Object VARiability) Spitzer Space Telescope observing program obtained the first extensive mid-infrared (3.6 & 4.5 um) time-series photometry of the Orion Nebula Cluster plus smaller footprints in eleven other star-forming cores (AFGL490, NGC1333, MonR2, GGD 12-15, NGC2264, L1688, Serpens Main, Serpens South, IRAS 20050+2720, IC1396A, and Ceph C). There are ~29,000 unique objects with light curves in either or both IRAC channels in the YSOVAR data set. We present the data collection and reduction for the Spitzer and ancillary data, and define the "standard sample" on which we calculate statistics, consisting of fast cadence data, with epochs about twice per day for ~40d. We also define a "standard sample of members", consisting of all the IR-selected members and X-ray selected members. We characterize the standard sample in terms of other properties, such as spectral energy distribution shape. We use three mechanisms to identify variables in the fast cadence data--the Stetson index, a chi^2 fit to a flat light curve, and significant periodicity. We also identified variables on the longest timescales possible of ~6 years, by comparing measurements taken early in the Spitzer mission with the mean from our YSOVAR campaign. The fraction of members in each cluster that are variable on these longest timescales is a function of the ratio of Class I/total members in each cluster, such that clusters with a higher fraction of Class I objects also have a higher fraction of long-term variables. For objects with a YSOVAR-determined period and a [3.6]-[8] color, we find that a star with a longer period is more likely than those with shorter periods to have an IR excess. We do not find any evidence for variability that causes [3.6]-[4.5] excesses to appear or vanish within our data; out of members and field objects combined, at most 0.02% may have transient IR excesses.

preprint2013arXiv

ALMA and CARMA observations of Brown Dwarfs disks: testing the models of dust evolution

The first steps toward planet formation involve the coagulation of small microscopic grains into larger and larger pebbles and rocks in gas-rich disks around young stars and brown dwarfs. Observations in the sub-millimeter can trace mm/cm-sized pebbles in the outer disks, and investigate the mechanisms of coagulation/fragmentation and radial migration of these solids. These represent key, yet not fully understood ingredients for our understanding of the formation of planetesimals, the building blocks of planets. Here we present the first results from an observational program using the ALMA and CARMA sub-mm/mm interferometers aimed at characterizing the dust properties and disk structure of young disks around brown dwarfs and very low mass stars. Given the physical conditions expected for these disks, they represent critical test beds for the models of the early stages of planet formation in proto-planetary disks.

preprint2013arXiv

CARMA interferometric observations of 2MASS J044427+2512: the first spatially resolved observations of thermal emission of a brown dwarf disk

We present CARMA 1.3 mm continuum data of the disk surrounding the young brown dwarf 2MASS J044427+2512 in the Taurus molecular cloud. The high angular resolution of the CARMA observations (0.16 arcsec) allows us to spatially resolve for the first time the thermal emission from dust around a brown dwarf. We analyze the interferometric visibilities and constrain the disk outer radius adopting disk models with power-law radial profiles of the dust surface density. In the case of a power-law index equal to or lower than 1, we obtain a disk radius in the range of about 15 - 30 AU, while larger disks are inferred for steeper radial profiles. By combining this information on the disk spatial extent with the sub-mm spectral index of this source we find conclusive evidence for mm-sized grains, or larger, in this brown dwarf disk. We discuss the implications of our results on the models of dust evolution in proto-planetary disks and brown dwarf formation.

preprint2013arXiv

The Ionized Circumstellar Envelopes of Orion Source I and the Becklin Neugebauer Object

The 229 GHz (lambda 1.3mm) radio emission from Orion-KL was mapped with up to 0.14" angular resolution with CARMA, allowing measurements of the flux densities of Source I ('SrcI') and the Becklin-Neugebauer Object (BN), the 2 most massive stars in this region. We find integrated flux densities of 310 +/- 45 mJy for SrcI and 240 +/- 35 mJy for BN. SrcI is optically thick even at 229 GHz. No trace of the H30alpha recombination line is seen in its spectrum, although the v_2=1, 5(5,0)-6(4,3) transition of H2O, 3450 K above the ground state, is prominent. SrcI is elongated at position angle 140 degrees, as in 43 GHz images. These results are most easily reconciled with models in which the radio emission from SrcI arises via the H- free-free opacity in a T < 4500 K disk, as considered by Reid et al. (2007). By contrast, the radio spectrum of BN is consistent with p+/e- free-free emission from a dense (n_e ~ 5x10^7 cm^{-3}), but otherwise conventional, hypercompact HII region. The source is becoming optically thin at 229 GHz, and the H30alpha recombination line, at VLSR = 23.2 +/- 0.5 km/sec, is prominent in its spectrum. A Lyman continuum flux of 5x10^{45} photons/sec, consistent with that expected from a B star, is required to maintain the ionization. Supplementary 90 GHz observations were made to measure the H41alpha and H42alpha recombination lines toward BN. Published 43 and 86 GHz data suggest that SrcI brightened with respect to BN over the 15 year period from 1994 to 2009.

preprint2012arXiv

Confirming the Primarily Smooth Structure of the Vega Debris Disk at Millimeter Wavelengths

Clumpy structure in the debris disk around Vega has been previously reported at millimeter wavelengths and attributed to concentrations of dust grains trapped in resonances with an unseen planet. However, recent imaging at similar wavelengths with higher sensitivity has disputed the observed structure. We present three new millimeter-wavelength observations that help to resolve the puzzling and contradictory observations. We have observed the Vega system with the Submillimeter Array (SMA) at a wavelength of 880 um and angular resolution of 5"; with the Combined Array for Research in Millimeter-wave Astronomy (CARMA) at a wavelength of 1.3 mm and angular resolution of 5"; and with the Green Bank Telescope (GBT) at a wavelength of 3.3 mm and angular resolution of 10". Despite high sensitivity and short baselines, we do not detect the Vega debris disk in either of the interferometric data sets (SMA and CARMA), which should be sensitive at high significance to clumpy structure based on previously reported observations. We obtain a marginal (3-sigma) detection of disk emission in the GBT data; the spatial distribution of the emission is not well constrained. We analyze the observations in the context of several different models, demonstrating that the observations are consistent with a smooth, broad, axisymmetric disk with inner radius 20-100 AU and width >50 AU. The interferometric data require that at least half of the 860 um emission detected by previous single-dish observations with the James Clerk Maxwell Telescope be distributed axisymmetrically, ruling out strong contributions from flux concentrations on spatial scales of <100 AU. These observations support recent results from the Plateau de Bure Interferometer indicating that previous detections of clumpy structure in the Vega debris disk were spurious.

preprint2012arXiv

HerMES: Candidate Gravitationally Lensed Galaxies and Lensing Statistics at Submillimeter Wavelengths

We present a list of 13 candidate gravitationally lensed submillimeter galaxies (SMGs) from 95 square degrees of the Herschel Multi-tiered Extragalactic Survey, a surface density of 0.14\pm0.04deg^{-2}. The selected sources have 500um flux densities (S_500) greater than 100mJy. Gravitational lensing is confirmed by follow-up observations in 9 of the 13 systems (70%), and the lensing status of the four remaining sources is undetermined. We also present a supplementary sample of 29 (0.31\pm0.06deg^{-2}) gravitationally lensed SMG candidates with S_500=80--100mJy, which are expected to contain a higher fraction of interlopers than the primary candidates. The number counts of the candidate lensed galaxies are consistent with a simple statistical model of the lensing rate, which uses a foreground matter distribution, the intrinsic SMG number counts, and an assumed SMG redshift distribution. The model predicts that 32--74% of our S_500>100mJy candidates are strongly gravitationally lensed (mu>2), with the brightest sources being the most robust; this is consistent with the observational data. Our statistical model also predicts that, on average, lensed galaxies with S_500=100mJy are magnified by factors of ~9, with apparently brighter galaxies having progressively higher average magnification, due to the shape of the intrinsic number counts. 65% of the sources are expected to have intrinsic 500micron flux densities less than 30mJy. Thus, samples of strongly gravitationally lensed SMGs, such as those presented here, probe below the nominal Herschel detection limit at 500 micron. They are good targets for the detailed study of the physical conditions in distant dusty, star-forming galaxies, due to the lensing magnification, which can lead to spatial resolutions of ~0.01" in the source plane.

preprint2012arXiv

YSOVAR: Six pre-main-sequence eclipsing binaries in the Orion Nebula Cluster

Eclipsing binaries (EBs) provide critical laboratories for empirically testing predictions of theoretical models of stellar structure and evolution. Pre-main-sequence (PMS) EBs are particularly valuable, both due to their rarity and the highly dynamic nature of PMS evolution, such that a dense grid of PMS EBs is required to properly calibrate theoretical PMS models. Analyzing multi-epoch, multi-color light curves for 2400 candidateOrion Nebula Cluster (ONC) members from our Warm Spitzer Exploration Science Program YSOVAR, we have identified 12 stars whose light curves show eclipse features. Four of these 12 EBs are previously known. Supplementing our light curves with follow-up optical and near-infrared spectroscopy, we establish two of the candidates as likely field EBs lying behind the ONC. We confirm the remaining six candidate systems, however, as newly identified ONC PMS EBs. These systems increase the number of known PMS EBs by over 50% and include the highest mass (Theta1 Ori E, for which we provide a complete set of well-determined parameters including component masses of 2.807 and 2.797 solar masses) and longest period (ISOY J053505.71-052354.1, P \sim 20 days) PMS EBs currently known. In two cases (Theta1 Ori E and ISOY J053526.88-044730.7), enough photometric and spectroscopic data exist to attempt an orbit solution and derive the system parameters. For the remaining systems, we combine our data with literature information to provide a preliminary characterization sufficient to guide follow-up investigations of these rare, benchmark systems.

preprint2011arXiv

Discovery of a Multiply-Lensed Submillimeter Galaxy in Early HerMES Herschel/SPIRE Data

We report the discovery of a bright ($f(250\mum) > 400$ mJy), multiply-lensed submillimeter galaxy \obj\ in {\it Herschel}/SPIRE Science Demonstration Phase data from the HerMES project. Interferometric 880\mum\ Submillimeter Array observations resolve at least four images with a large separation of $\sim 9\arcsec$. A high-resolution adaptive optics $K_p$ image with Keck/NIRC2 clearly shows strong lensing arcs. Follow-up spectroscopy gives a redshift of $z=2.9575$, and the lensing model gives a total magnification of $μ\sim 11 \pm 1$. The large image separation allows us to study the multi-wavelength spectral energy distribution (SED) of the lensed source unobscured by the central lensing mass. The far-IR/millimeter-wave SED is well described by a modified blackbody fit with an unusually warm dust temperature, $88 \pm 3$ K. We derive a lensing-corrected total IR luminosity of $(1.43 \pm 0.09) \times 10^{13}\, \mathrm{L}_{\odot}$, implying a star formation rate of $\sim 2500\, \mathrm{M}_{\odot}\, \mathrm{yr}^{-1}$. However, models primarily developed from brighter galaxies selected at longer wavelengths are a poor fit to the full optical-to-millimeter SED. A number of other strongly lensed systems have already been discovered in early {\it Herschel} data, and many more are expected as additional data are collected.

preprint2011arXiv

Dynamical Structure of the Molecular Interstellar Medium in an Extremely Bright, Multiply Lensed z~3 Submillimeter Galaxy Discovered with Herschel

We report the detection of CO 5-4, 3-2, and 1-0 emission in the strongly lensed, Herschel/SPIRE-selected submillimeter galaxy (SMG) HLSW-01 at z=2.9574+/-0.0001, using the Plateau de Bure Interferometer, the Combined Array for Research in Millimeter-wave Astronomy, and the Green Bank Telescope. The observations spatially resolve the molecular gas into four lensed images with a maximum separation of ~9", and reveal the internal gas dynamics in this system. We derive lensing-corrected CO line luminosities of L'(CO 1-0) = (4.17+/-0.41), L'(CO 3-2) = (3.96+/-0.20) and L'(CO 5-4) = (3.45+/-0.20) x 10^10 (mu_L/10.9)^-1 Kkm/s pc^2, corresponding to luminosity ratios of r_31 = 0.95+/-0.10, r_53 = 0.87+/-0.06, and r_51 = 0.83+/-0.09. This suggests a total molecular gas mass of Mgas = 3.3 x 10^10 (alpha_CO/0.8) (mu_L/10.9)^-1 Msun. The gas mass, gas mass fraction, gas depletion timescale, star formation efficiency, and specific star formation rate are typical for an SMG. The velocity structure of the gas reservoir suggests that the brightest two lensed images are dynamically resolved projections of the same dust-obscured region in the galaxy that are kinematically offset from the unresolved fainter images. The resolved kinematics appear consistent with the complex velocity structure observed in major, `wet' (i.e., gas-rich) mergers. Major mergers are commonly observed in SMGs, and are likely to be responsible for fueling their intense starbursts at high gas consumption rates. This study demonstrates the level of detail to which galaxies in the early universe can be studied by utilizing the increase in effective spatial resolution and sensitivity provided by gravitational lensing.

preprint2011arXiv

Redshift Determination and CO Line Excitation Modeling for the Multiply-Lensed Galaxy HLSW-01

We report on the redshift measurement and CO line excitation of HERMES J105751.1+573027 (HLSW-01), a strongly lensed submillimeter galaxy discovered in Herschel/SPIRE observations as part of the Herschel Multi-tiered Extragalactic Survey (HerMES). HLSW-01 is an ultra-luminous galaxy with an intrinsic far-infrared luminosity of 1.4x10^(13) solar luminosities, and is lensed by a massive group of galaxies into at least four images with a total magnification of 10.9+/-0.7. With the 100 GHz instantaneous bandwidth of the Z-Spec instrument on the Caltech Submillimeter Observatory, we robustly identify a redshift of z=2.958+/-0.007 for this source, using the simultaneous detection of four CO emission lines (J = 7-6, J = 8-7, J = 9-8, and J = 10-9). Combining the measured line fluxes for these high-J transitions with the J = 1-0, J = 3-2 and J = 5-4 line fluxes measured with the Green Bank Telescope, the Combined Array for Research in Millimeter Astronomy, and the Plateau de Bure Interferometer, respectively, we model the physical properties of the molecular gas in this galaxy. We find that the full CO spectral line energy distribution is well described by warm, moderate-density gas with Tkin = 86-235 K and n(H2) = (1.1-3.5)x10^3 cm^(-3). However, it is possible that the highest-J transitions are tracing a small fraction of very dense gas in molecular cloud cores, and two-component models that include a warm/dense molecular gas phase with Tkin ~ 200 K, n(H2) ~ 10^5 cm^(-3) are also consistent with these data. Higher signal-to-noise measurements of the J(upper) > 7 transitions with high spectral resolution, combined with high spatial resolution CO maps, are needed to improve our understanding of the gas excitation, morphology, and dynamics of this interesting high-redshift galaxy.

preprint2011arXiv

The Birth of a Relativistic Outflow in the Unusual γ-ray Transient Swift J164449.3+573451

Active galactic nuclei (AGN), powered by long-term accretion onto central supermassive black holes, produce relativistic jets with lifetimes of greater than one million yr that preclude observations at birth. Transient accretion onto a supermassive black hole, for example through the tidal disruption of a stray star, may therefore offer a unique opportunity to observe and study the birth of a relativistic jet. On 2011 March 25, the Swift γ-ray satellite discovered an unusual transient source (Swift J164449.3+573451) potentially representing such an event. Here we present the discovery of a luminous radio transient associated with Swift J164449.3+573451, and an extensive set of observations spanning centimeter to millimeter wavelengths and covering the first month of evolution. These observations lead to a positional coincidence with the nucleus of an inactive galaxy, and provide direct evidence for a newly-formed relativistic outflow, launched by transient accretion onto a million solar mass black hole. While a relativistic outflow was not predicted in this scenario, we show that the tidal disruption of a star naturally explains the high-energy properties, radio luminosity, and the inferred rate of such events. The weaker beaming in the radio compared to γ-rays/X-rays, suggests that radio searches may uncover similar events out to redshifts of z ~ 6.

preprint2011arXiv

The North American and Pelican Nebulae II. MIPS Observations and Analysis

We present observations of ~7 square degrees of the North American and Pelican Nebulae region at 24, 70, and 160 microns with the Spitzer Space Telescope Multiband Imaging Photometer for Spitzer (MIPS). We incorporate the MIPS observations with earlier Spitzer Infrared Array Camera (IRAC) observations, as well as archival near-infrared (IR) and optical data. We use the MIPS data to identify 1286 young stellar object (YSO) candidates. IRAC data alone can identify 806 more YSO candidates, for a total of 2076 YSO candidates. Prior to the Spitzer observations, there were only ~200 YSOs known in this region. Three subregions within the complex are highlighted as clusters: the Gulf of Mexico, the Pelican, and the Pelican's Hat. The Gulf of Mexico cluster is subject to the highest extinction (Av at least ~30) and has the widest range of infrared colors of the three clusters, including the largest excesses and by far the most point-source detections at 70 microns. Just 3% of the cluster members were previously identified; we have redefined this cluster as about 10-100 times larger (in projected area) than was previously realized.

preprint2011arXiv

YSOVAR: the first sensitive, wide-area, mid-IR photometric monitoring of the ONC

We present initial results from time series imaging at infrared wavelengths of 0.9 sq. degrees in the Orion Nebula Cluster (ONC). During Fall 2009 we obtained 81 epochs of Spitzer 3.6 and 4.5 micron data over 40 consecutive days. We extracted light curves with ~3% photometric accuracy for ~2000 ONC members ranging from several solar masses down to well below the hydrogen burning mass limit. For many of the stars, we also have time-series photometry obtained at optical (Ic) and/or near-infrared (JKs) wavelengths. Our data set can be mined to determine stellar rotation periods, identify new pre-main-sequence (PMS) eclipsing binaries, search for new substellar Orion members, and help better determine the frequency of circumstellar disks as a function of stellar mass in the ONC. Our primary focus is the unique ability of 3.6 & 4.5 micron variability information to improve our understanding of inner disk processes and structure in the Class I and II young stellar objects (YSOs). In this paper, we provide a brief overview of the YSOVAR Orion data obtained in Fall 2009, and we highlight our light curves for AA-Tau analogs - YSOs with narrow dips in flux, most probably due to disk density structures passing through our line of sight. Detailed follow-up observations are needed in order to better quantify the nature of the obscuring bodies and what this implies for the structure of the inner disks of YSOs.

preprint2010arXiv

Recurring millimeter flares as evidence for star-star magnetic reconnection events in the DQ Tauri PMS binary system

Observations of the T Tauri spectroscopic binary DQ Tau in April 2008 captured an unusual flare at 3 mm, which peaked at an observed max flux of 0.5 Jy (about 27x the quiescent value). Here we present follow-up mm observations that demonstrate a periodicity to the phenomenon. While monitoring 3 new periastron encounters, we detect flares within 17.5 hrs (or 4.6%) of the orbital phase of the first reported flare, and we constrain the main emitting region to a stellar height of 3.7-6.8 Rstar. The recorded activity is consistent with the proposed picture for synchrotron emission initiated by a magnetic reconnection event when the two stellar magnetospheres of the highly eccentric (e=0.556) binary are believed to collide near periastron as the stars approach a minimum separation of 8 Rstar (~13 Rsolar). The similar light curve decay profiles allow us to estimate an average flare duration of 30 hrs. Assuming one mm flare per orbit, DQ Tau could spend approximately 8% of its 15.8-d orbital period in an elevated flux state. Our analysis of the mm emission provides an upper limit of 5% on the linear polarization. We discuss the extent to which a severely entangled magnetic field structure and Faraday rotation effects are likely to reduce the observed polarization fraction. We also predict that, for the current picture, the stellar magnetospheres must be misaligned at a significant angle or, alternatively, that the topologies of the outer magnetospheres are poorly described by a well-ordered dipole inside a radius of 7 Rstar. Finally, to investigate whether reorganization of the magnetic field during the interaction affects mass accretion, we also present simultaneous optical (VRI) monitoring, as an established tracer of accretion activity in this system. We find that an accretion event can occur coincident in both time and duration with the synchrotron fallout of a magnetic reconnection event.

preprint2009arXiv

The North American and Pelican Nebulae I. IRAC Observations

We present a 9 deg^2 map of the North American and Pelican Nebulae regions obtained in all four IRAC channels with the Spitzer Space Telescope. The resulting photometry is merged with that at JHKs from 2MASS and a more spatially limited $BVI$ survey from previous ground-based work. We use a mixture of color- color diagrams to select a minimally contaminated set of more than 1600 objects that we claim are young stellar objects (YSOs) associated with the star forming region. Because our selection technique uses IR excess as a requirement, our sample is strongly biased against inclusion of Class III YSOs. The distribution of IRAC spectral slopes for our YSOs indicates that most of these objects are Class II, with a peak towards steeper spectral slopes but a substantial contribution from a tail of flat spectrum and Class I type objects. By studying the small fraction of the sample that is optically visible, we infer a typical age of a few Myr for the low mass population. The young stars are clustered, with about a third of them located in eight clusters that are located within or near the LDN 935 dark cloud. Half of the YSOs are located in regions with surface densities higher than 1000 YSOs / deg^2. The Class I objects are more clustered than the Class II stars.

preprint2008arXiv

2MASS J01542930+0053266 : A New Eclipsing M-dwarf Binary System

We report on 2MASS J01542930+0053266, a faint eclipsing system composed of two M dwarfs. The variability of this system was originally discovered during a pilot study of the 2MASS Calibration Point Source Working Database. Additional photometry from the Sloan Digital Sky Survey yields an 8-passband lightcurve, from which we derive an orbital period of 2.6390157 +/- 0.0000016 days. Spectroscopic followup confirms our photometric classification of the system, which is likely composed of M0 and M1 dwarfs. Radial velocity measurements allow us to derive the masses (M_1 = 0.66 +/- 0.03 M_sun; M_2 = 0.62 +/- 0.03 M_sun) and radii (R_1 = 0.64 +/- 0.08 R_sun; R_2 = 0.61 +/- 0.09 R_sun) of the components, which are consistent with empirical mass-radius relationships for low-mass stars in binary systems. We perform Monte Carlo simulations of the lightcurves which allow us to uncover complicated degeneracies between the system parameters. Both stars show evidence of H-alpha emission, something not common in early-type M dwarfs. This suggests that binarity may influence the magnetic activity properties of low-mass stars; activity in the binary may persist long after the dynamos in their isolated counterparts have decayed, yielding a new potential foreground of flaring activity for next generation variability surveys.

preprint2007arXiv

Detection of [Ne II] Emission from Young Circumstellar Disks

We report the detection of [Ne II] emission at 12.81 micron in four out of the six optically thick dust disks observed as part of the FEPS Spitzer Legacy program. In addition, we detect a H I(7-6) emission line at 12.37 micron from the source RXJ1852.3-3700. Detections of [Ne II] lines are favored by low mid-infrared excess emission. Both stellar X-rays and extreme UV (EUV) photons can sufficiently ionize the disk surface to reproduce the observed line fluxes, suggesting that emission from Ne+ originates in the hot disk atmosphere. On the other hand, the H I(7-6) line is not associated with the gas in the disk surface and magnetospheric accretion flows can account only for at most ~30% of the observed flux. We conclude that accretion shock regions and/or the stellar corona could contribute to most of the H I(7-6)emission. Finally, we discuss the observations necessary to identify whether stellar X-rays or EUV photons are the dominant ionization mechanism for Ne atoms. Because the observed [Ne II] emission probes very small amounts of gas in the disk surface (~10^{-6} Jupiter masses) we suggest using this gas line to determine the presence or absence of gas in more evolved circumstellar disks.