Tracy L. Beck

HERschel Inventory of The Agents of Galaxy Evolution (HERITAGE): The Large Magellanic Cloud dust

The HERschel Inventory of The Agents of Galaxy Evolution (HERITAGE) of the Magellanic Clouds will use dust emission to investigate the life cycle of matter in both the Large and Small Magellanic Clouds (LMC and SMC). Using the Herschel Space Observatory’s PACS and SPIRE photometry cameras, we imaged a 2° × 8° strip through the LMC, at a position angle of ~22.5° as part of the science demonstration phase of the Herschel mission. We present the data in all 5 Herschel bands: PACS 100 and 160 μm and SPIRE 250, 350 and 500 μm. We present two dust models that both adequately fit the spectral energy distribution for the entire strip and both reveal that the SPIRE 500 μm emission is in excess of the models by ~6 to 17%. The SPIRE emission follows the distribution of the dust mass, which is derived from the model. The PAH-to-dust mass (f_(PAH)) image of the strip reveals a possible enhancement in the LMC bar in agreement with previous work. We compare the gas mass distribution derived from the HI 21 cm and CO J = 1−0 line emission maps to the dust mass map from the models and derive gas-to-dust mass ratios (GDRs). The dust model, which uses the standard graphite and silicate optical properties for Galactic dust, has a very low GDR = 65^(+15) _(−18) making it an unrealistic dust model for the LMC. Our second dust model, which uses amorphous carbon instead of graphite, has a flatter emissivity index in the submillimeter and results in a GDR = 287^_(+25)_(−42) that is more consistent with a GDR inferred from extinction.

Feeding versus feedback in NGC 4151 probed with Gemini NIFS – II. Kinematics

We have used the Gemini Near-infrared Integral Field Spectrograph (NIFS) to map the gas kinematics of the inner ∼200 x 500 pc 2 of the Seyfert galaxy NGC 4151 in the Z, J, H and K bands at a resolving power of ≥5000 and spatial resolution of ∼8 pc. The ionized gas emission is most extended along the known ionization bi-cone at position angle PA = 60°-240°, but is observed also along its equatorial plane. This indicates that the active galactic nucleus (AGN) ionizes gas beyond the borders of the bi-cone, within a sphere with ≈1 arcsec radius around the nucleus. The ionized gas has three kinematic components: (1) one observed at the systemic velocity and interpreted as originating in the galactic disc; (2) one outflowing along the bi-cone, with line-of-sight velocities between -600 and 600 km s ―1 and strongest emission at ±(100―300) km s ―1 ; and (3) another component due to the interaction of the radio jet with ambient gas. The radio jet (at PA = 75°-255°) is not aligned with the narrow-line region (NLR) and produces flux enhancements mostly observed at the systemic velocity, suggesting that the jet is launched close to the plane of the galaxy (approximately plane of the sky). The mass outflow rate, estimated to be ≈ 1 M ⊙ yr ―1 along each cone, exceeds the inferred black hole accretion rate by a factor of ∼100. This can be understood if the NLR is formed mostly by entrained gas from the circumnuclear interstellar medium by an outflow probably originating in the accretion disc. This flow represents feedback from the AGN, estimated to release a kinetic power of Ė ≈ 2.4 x 10 41 erg s ―1 , which is only ∼0.3 per cent of the bolometric luminosity of the AGN. There is no evidence in our data for the gradual acceleration followed by gradual deceleration proposed by previous modelling of the [O III] emitting gas. Our data allow the possibility that the NLR clouds are accelerated close to the nucleus (within 0.1 arcsec, which corresponds to ≈7 pc at the galaxy) after which the flow moves at essentially constant velocity (≈600 km s ―1 ), being consistent with near-infrared emission arising predominantly from the interaction of the outflow with gas in the galactic disc. The molecular gas exhibits distinct kinematics relative to the ionized gas. Its emission arises in extended regions approximately perpendicular to the axis of the bi-cone and along the axis of the galaxys stellar bar, avoiding the innermost ionized regions. It does not show an outflowing component, being observed only at velocities very close to systemic, and is thus consistent with an origin in the galactic plane. This hot molecular gas may only be the tracer of a larger reservoir of colder gas which represents the AGN feeding.

Feeding versus feedback in NGC 4151 probed with Gemini NIFS – I. Excitation

We have used the Gemini Near-infrared Integral Field Spectrograph (NIFS) to map the emission-line intensity distributions and ratios in the narrow-line region (NLR) of the Seyfert galaxy NGC 4151 in the Z, J, H and K bands at a resolving power ≥5000, covering the inner ≈200 x 300 pc of the galaxy at a spatial resolution of ≈8 pc. We present intensity distributions in 14 emission lines, which show three distinct behaviours. (1) Most of the ionized gas intensity distributions are extended to ≈1100 pc from the nucleus along the region covered by the known biconical outflow (position angle, PA = 60/240°, NE-SW), consistent with an origin in the outflow; while the recombination lines show intensity profiles which decrease with distance r from the nucleus as I ∝ r -1 , most of the forbidden lines present a flat intensity profile (I ∝ r 0 ) or even increasing with distance from the nucleus towards the border of the NLR. (2) The H 2 emission lines show completely distinct intensity distributions, which avoid the region of the bicone, extending from ≈10 to ≈60 pc from the nucleus approximately along the large-scale bar, almost perpendicular to the bicone axis. This morphology supports an origin for the H 2 -emitting gas in the galaxy plane. (3) The coronal lines show a steep intensity profile, described by I ∝ r -2 ; the emission is clearly resolved only in the case of [Si VII], consistent with an origin in the inner NLR. Using the line-ratio maps [Fe II] 1.644/1.257 and Pa β/Br γ, we obtain an average reddening of E(B - V) ≈ 0.5 along the NLR and E(B - V) ≥ 1 at the nucleus. Our line-ratio map [Pe II] 1.257 μm/[P II] 1.189 μm of the NLR of NGC4151 is the first such map of an extragalactic source. Together with the [Fe II]/Paβ map, these line ratios correlate with the radio intensity distribution, mapping the effects of shocks produced by the radio jet on the NLR. These shocks probably release the Fe locked in grains and produce an enhancement of the [Fe II] emission at ≈1 arcsec from the nucleus. At these regions, we obtain electron densities N e ≈4000 cm -3 and temperatures T e ≈ 15 000 K for the [Fe II]-emitting gas. For the H 2 -emitting gas, we obtain much lower temperatures of T exc ≈ 2100 K and conclude that the gas is in thermal equilibrium. The heating necessary to excite the molecule may be due to X-rays escaping perpendicular to the cone (through the nuclear torus, if there is one) or to shocks probably produced by the accretion flow previously observed along the large-scale bar. The distinct intensity distributions and physical properties of the ionized and molecular gas, as well as their locations, the former along the outflowing cone, and the latter in the galaxy plane surrounding the nucleus, suggest that the H 2 -emitting gas traces the active galactic nuclei feeding, while the ionized gas traces its feedback.

Adaptive Optics Imaging of the Orion Trapezium Cluster

We have imaged an area ~5 arcmin2 at the center of the Trapezium cluster in Orion in the K band using the University of Hawaii (UH) adaptive optics system at the UH 2.2 m telescope. Our survey detects 292 stars brighter than K = 18.2 mag and resolves pairs to the 023 diffraction limit of the telescope. The binary fraction in the angular separation range 03–06, corresponding to 132–264 AU at 440 pc, is indistinguishable from that of the solar-like stars in the solar neighborhood. Proplyds are associated with both single stars and visual binaries. About half the stars in our sample have also been measured at V and I by Prosser et al.; most of these seem to be about 106 yr old, observed through moderate extinction, and having some excess emission at K. The luminosity function turns over before the stellar limit of ~13.5 mag, indicating that the cluster does not contain a large population of massive brown dwarfs.

The HERschel Inventory of the Agents of Galaxy Evolution in the Magellanic Clouds, a HERschel Open Time Key Program

We present an overview of the HERschel Inventory of The Agents of Galaxy Evolution (HERITAGE) in the Magellanic Clouds project, which is a Herschel Space Observatory open time key program. We mapped the Large Magellanic Cloud (LMC) and Small Magellanic Cloud (SMC) at 100, 160, 250, 350, and 500 μm with the Spectral and Photometric Imaging Receiver (SPIRE) and Photodetector Array Camera and Spectrometer (PACS) instruments on board Herschel using the SPIRE/PACS parallel mode. The overriding science goal of HERITAGE is to study the life cycle of matter as traced by dust in the LMC and SMC. The far-infrared and submillimeter emission is an effective tracer of the interstellar medium (ISM) dust, the most deeply embedded young stellar objects (YSOs), and the dust ejected by the most massive stars. We describe in detail the data processing, particularly for the PACS data, which required some custom steps because of the large angular extent of a single observational unit and overall the large amount of data to be processed as an ensemble. We report total global fluxes for the LMC and SMC and demonstrate their agreement with measurements by prior missions. The HERITAGE maps of the LMC and SMC are dominated by the ISM dust emission and bear most resemblance to the tracers of ISM gas rather than the stellar content of the galaxies. We describe the point source extraction processing and the criteria used to establish a catalog for each waveband for the HERITAGE program. The 250 μm band is the most sensitive and the source catalogs for this band have ~25,000 objects for the LMC and ~5500 objects for the SMC. These data enable studies of ISM dust properties, submillimeter excess dust emission, dust-to-gas ratio, Class 0 YSO candidates, dusty massive evolved stars, supernova remnants (including SN1987A), H II regions, and dust evolution in the LMC and SMC. All images and catalogs are delivered to the Herschel Science Center as part of the community support aspects of the project. These HERITAGE images and catalogs provide an excellent basis for future research and follow up with other facilities.

The Infrared Light Curve of SN 2011fe in M101 and the Distance to M101

We present near-infrared light curves of supernova (SN) 2011fe in M101, including 34 epochs in H band starting 14 days before maximum brightness in the B band. The light curve data were obtained with the WIYN High-Resolution Infrared Camera. When the data are calibrated using templates of other Type Ia SNe, we derive an apparent H-band magnitude at the epoch of B-band maximum of 10.85 ± 0.04. This implies a distance modulus for M101 that ranges from 28.86 to 29.17 mag, depending on which absolute calibration for Type Ia SNe is used.

V1647 ORIONIS: ONE YEAR INTO QUIESCENCE

We present new optical, near-IR, and mid-IR observations of the young eruptive variable star V1647 Orionis that went into outburst in late 2004 for approximately two years. Our observations, taken one year after the star had faded to its pre-outburst optical brightness, show that V1647 Ori is still actively accreting circumstellar material. We compare and contrast these data with existing observations of the source from both pre-outburst and outburst phases. From near-IR spectroscopy we identify photospheric absorption features for the first time that allow us to constrain the classification of the young star itself. Our best-fit spectral type is M0 ± 2 sub-classes with a visual extinction of 19 ± 2 magnitudes and a K-band veiling of rK ~ 1.5 ± 0.2. We estimate that V1647 Ori has a quiescent bolometric luminosity of ~ 9.5 L☉ and a mass accretion rate of ~ 1 × 10−6 M☉ yr−1. Our derived mass and age, from comparison with evolutionary models, are 0.8 ± 0.2 M☉ and 0.5 Myr, respectively. The presence toward the star of shock-excited optical [S II] and [Fe II] emission as well as near-IR H2 and [Fe II] emission perhaps suggests that a new Herbig-Haro flow is becoming visible close to the star.

Possible planet formation in the young, low-mass, multiple stellar system GG Tau A.

The formation of planets around binary stars may be more difficult than around single stars. In a close binary star (with a separation of less than a hundred astronomical units), theory predicts the presence of circumstellar disks around each star, and an outer circumbinary disk surrounding a gravitationally cleared inner cavity around the stars. Given that the inner disks are depleted by accretion onto the stars on timescales of a few thousand years, any replenishing material must be transferred from the outer reservoir to fuel planet formation (which occurs on timescales of about one million years). Gas flowing through disk cavities has been detected in single star systems. A circumbinary disk was discovered around the young low-mass binary system GG Tau A (ref. 7), which has recently been shown to be a hierarchical triple system. It has one large inner disk around the single star, GG Tau Aa, and shows small amounts of shocked hydrogen gas residing within the central cavity, but other than a single weak detection, the distribution of cold gas in this cavity or in any other binary or multiple star system has not hitherto been determined. Here we report imaging of gas fragments emitting radiation characteristic of carbon monoxide within the GG Tau A cavity. From the kinematics we conclude that the flow appears capable of sustaining the inner disk (around GG Tau Aa) beyond the accretion lifetime, leaving time for planet formation to occur there. These results show the complexity of planet formation around multiple stars and confirm the general picture predicted by numerical simulations.

Dynamical Mass Estimates for Incomplete Orbits: Young Multiple Stars in Taurus and Ophiuchus

We present recent measurements of the orbital motion in the binaries DF Tau and ZZ Tau and in the triples Elias 12, T Tau, and V853 Oph. We observed these systems with the Fine Guidance Sensors on the Hubble Space Telescope and with adaptive optics imaging at the W. M. Keck and Gemini North Observatories. Based on our measurements and those presented in the literature we perform preliminary orbital analyses for DF Tau, ZZ Tau, Elias 12 Na-Nb, and T Tau Sa-Sb. Because the orbital coverage in most of these systems does not yet span a sufficient portion of the orbit, we are not able to find definitive orbit solutions. By using a Monte Carlo search technique we explored the orbital parameter space allowed by the current set of data available for each binary. We constructed weighted distributions for the total mass of the binaries derived from a large sample of possible orbits that fit the data. These mass distributions show that the total mass is already well defined. We compute total mass estimates of 0.78, 0.66, 1.13, and 4.13 M⊙ for DF Tau, ZZ Tau, Elias 12 Na-Nb, and T Tau Sa-Sb, respectively, using a distance of 140 pc. For Elias 12 Na-Nb, where the orbital coverage spans ~164°, we compute a preliminary orbit solution with a period of ~9-12 yr. By including an earlier lunar occultation measurement, we also find a likely orbit solution for ZZ Tau with a period of ~32 yr. With additional measurements to continue mapping the orbits the derived dynamical masses will be useful in constraining the theoretical tracks of pre-main-sequence evolution.

Kinematics and excitation of the molecular hydrogen accretion disc in NGC 1275

We report the results of high spatial and spectral resolution integral-field spectroscopy of the central ~3 x 3 arcsec^2 of the active galaxy NGC 1275 (Perseus A), based on observations with the Near-infrared Integral Field Spectrograph (NIFS) and the ALTAIR adaptive-optics system on the Gemini North telescope. The circum-nuclear disc in the inner R~50 pc of NGC 1275 is seen in both the H2 and [FeII] lines. The disc is interpreted as the outer part of a collisionally-excited turbulent accretion disc. The kinematic major axis of the disc at a position angle of 68 deg is oriented perpendicular to the radio jet. A streamer-like feature to the south-west of the disc, detected in H2 but not in [FeII], is discussed as one of possibly several molecular streamers, presumably falling into the nuclear region. Indications of an ionization structure within the disc are deduced from the HeI and Br gamma emission lines, which may partially originate from the inner portions of the accretion disc. The kinematics of these two lines agrees with the signature of the circum-nuclear disc, but both lines display a larger central velocity dispersion than the H2 line. The rovibrational H2 transitions from the core of NGC 1275 are indicative of thermal excitation caused by shocks and agree with excitation temperatures of ~1360 and ~4290 K for the lower- and higher-energy H2 transitions, respectively. The data suggest X-ray heating as the dominant excitation mechanism of [FeII] emission in the core, while fast shocks are a possible alternative. The [FeII] lines indicate an electron density of ~4000 cm^{-3}. The H2 disc is modelled using simulated NIFS data cubes of H2 emission from inclined discs in Keplerian rotation around a central mass. Assuming a disc inclination of 45 deg +/- 10 deg, the best-fitting models imply a central mass of (8^{+7}_{-2}) x 10^8 Msun. (abridged)