Source author record

George H. Rieke

George H. Rieke 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

60works
7topics
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

60 published item(s)

preprint2022arXiv

A Star-sized Impact-Produced Dust Clump in the Terrestrial Zone of the HD 166191 System

We report on five years of 3-5 micron photometry measurements obtained by warm Spitzer to track the dust debris emission in the terrestrial zone of HD 166191 in combination with simultaneous optical data. We show that the debris production in this young (~10 Myr) system increased significantly in early 2018 and reached a record high level (almost double by mid 2019) by the end of the Spitzer mission (early 2020), suggesting intense collisional activity in its terrestrial zone likely due to either initial assembling of terrestrial planets through giant impacts or dynamical shake-up from unseen planet-mass objects or recent planet migration. This intense activity is further highlighted by detecting a star-size dust clump, passing in front of the star, in the midst of its infrared brightening. We constrain the minimum size and mass of the clump using multiwavelength transit profiles and conclude that the dust clump is most likely created by a large impact involving objects of several hundred km in size with an apparent period of 142 days (i.e., 0.62 au assuming a circular orbit). The system's evolutionary state (right after the dispersal of its gas-rich disk) makes it extremely valuable to learn about the process of terrestrial planet formation and planetary architecture through future observations.

preprint2022arXiv

AGN Selection and Demographics in GOODS-S/HUDF from X-ray to Radio

We present a comprehensive census of the AGNs in the GOODS-S/HUDF region from the X-ray to the radio, covering both the obscured and unobscured populations. This work includes a robust analysis of the source optical-to-mid-IR SEDs featuring (semi-)empirical AGN and galaxy dust emission models and Baysian fitting techniques, ultra-deep VLA 3 and 6 GHz observations, and an integrated analysis of various AGN selection techniques, including X-ray properties, UV-to-MIR SED analysis, optical spectral features, mid-IR colors, radio loudness and spectral slope, and AGN variability. In total, we report $\sim$900 AGNs over the $\sim$170 arcmin$^2$ 3D-HST GOODS-S footprint, which has doubled the AGN number identified in the previous X-ray sample with $\sim$26\% of our sample undetected in the deepest Chandra image. With a summary of AGN demographics from different selection methods, we find that no one single band or technique comes close to selecting a complete AGN sample despite the great depth of the data in GOODS-S/HUDF. We estimate the yields of various approaches and explore the reasons for incompleteness. We characterize the statistical properties, such as source number density, obscuration fraction and luminosity function of the AGN sample in this field and discuss their immediate implications. We also provide some qualitative predictions of the AGN sample that might be discovered by the upcoming JWST surveys.

preprint2020arXiv

Dust emission from the first massive galaxies

We compare observed far-infrared (FIR) galaxy spectral energy distributions (SEDs) at $z\sim6$ and theoretical predictions for first massive population II (Pop II) galaxies. Observed FIR SEDs at $z\gtrsim5$ are broad and shifted to bluer wavelengths when compared to galaxies at $z\sim3$. By implementing an analytical model for dust emission from Pop II massive galaxies,we were able to reproduce the observed behaviour as a consequence of the high energy densities and silicate-rich dust composition of high-$z$ model galaxies. As noted in a previous work, the non-blackbody nature of galaxy SEDs at $z \sim 6$ should be taken into account when interpreting measurements of FIR luminosities to avoid underestimating star formation rates.

preprint2020arXiv

Mid-infrared Studies of HD 113766 and HD 172555: Assessing Variability in the Terrestrial Zone of Young Exoplanetary Systems

We present multi-epoch infrared photometry and spectroscopy obtained with warm Spitzer, Subaru and SOFIA to assess variability for the young ($\sim$20 Myr) and dusty debris systems around HD 172555 and HD 113766A. No variations (within 0.5%) were found for the former at either 3.6 or 4.5 $μ$m, while significant non-periodic variations (peak-to-peak of $\sim$10-15% relative to the primary star) were detected for the latter. Relative to the Spitzer IRS spectra taken in 2004, multi-epoch mid-infrared spectra reveal no change in either the shape of the prominent 10 $μ$m solid-state features or the overall flux levels (no more than 20%) for both systems, corroborating that the population of sub-$μ$m-sized grains that produce the pronounced solid-state features is stable over a decadal timescale. We suggest that these sub-$μ$m-sized grains were initially generated in an optically thick clump of debris of mm-sized vapor condensates resulting from a recent violent impact between large asteroidal or planetary bodies. Because of the shielding from the stellar photons provided by this clump, intense collisions led to an over-production of fine grains that would otherwise be ejected from the system by radiation pressure. As the clump is sheared by its orbital motion and becomes optically thin, a population of very fine grains could remain in stable orbits until Poynting-Robertson drag slowly spirals them into the star. We further suggest that the 3-5 $μ$m disk variation around HD 113766A is consistent with a clump/arc of such fine grains on a modestly eccentric orbit in its terrestrial zone.

preprint2020arXiv

New HST data and modeling reveal a massive planetesimal collision around Fomalhaut

The apparent detection of an exoplanet orbiting Fomalhaut was announced in 2008. However, subsequent observations of Fomalhaut b raised questions about its status: Unlike other exoplanets, it is bright in the optical and nondetected in the infrared, and its orbit appears to cross the debris ring around the star without the expected gravitational perturbations. We revisit previously published data and analyze additional Hubble Space Telescope (HST) data, finding that the source is likely on a radial trajectory and has faded and become extended. Dynamical and collisional modeling of a recently produced dust cloud yields results consistent with the observations. Fomalhaut b appears to be a directly imaged catastrophic collision between two large planetesimals in an extrasolar planetary system. Similar events should be very rare in quiescent planetary systems of the age of Fomalhaut, suggesting that we are possibly witnessing the effects of gravitational stirring due to the orbital evolution of hypothetical planet(s) around the star.

preprint2020arXiv

The quantum efficiency and diffractive image artifacts of Si:As IBC mid-IR detector arrays at 5 $-$ 10 $μ$m: Implications for the JWST/MIRI detectors

Arsenic doped back illuminated blocked impurity band (BIBIB) silicon detectors have advanced near and mid-IR astronomy for over thirty years; they have high quantum efficiency (QE), especially at wavelengths longer than 10 $μ$m, and a large spectral range. Their radiation hardness is also an asset for space based instruments. Three examples of Si:As BIBIB arrays are used in the Mid-InfraRed Instrument (MIRI) of the James Webb Space Telescope (JWST), observing between 5 and 28 $μ$m. In this paper, we analyze the parameters leading to high quantum efficiency (up to $\sim$ 60\%) for the MIRI devices between 5 and 10 $μ$m. We also model the cross-shaped artifact that was first noticed in the 5.7 and 7.8 $μ$m Spitzer/IRAC images and has since also been imaged at shorter wavelength ($\le 10~μ$m) laboratory tests of the MIRI detectors. The artifact is a result of internal reflective diffraction off the pixel-defining metallic contacts to the readout detector circuit. The low absorption in the arrays at the shorter wavelengths enables photons diffracted to wide angles to cross the detectors and substrates multiple times. This is related to similar behavior in other back illuminated solid-state detectors with poor absorption, such as conventional CCDs operating near 1 $μ$m. We investigate the properties of the artifact and its dependence on the detector architecture with a quantum-electrodynamic (QED) model of the probabilities of various photon paths. Knowledge of the artifact properties will be especially important for observations with the MIRI LRS and MRS spectroscopic modes.

preprint2016arXiv

A Comprehensive Dust Model Applied to the Resolved Beta Pictoris Debris Disk from Optical to Radio Wavelengths

We investigate whether varying the dust composition (described by the optical constants) can solve a persistent problem in debris disk modeling--the inability to fit the thermal emission without over-predicting the scattered light. We model five images of the beta Pictoris disk: two in scattered light from HST/STIS at 0.58 microns and HST/WFC3 at 1.16 microns, and three in thermal emission from Spitzer/MIPS at 24 microns, Herschel/PACS at 70 microns, and ALMA at 870 microns. The WFC3 and MIPS data are published here for the first time. We focus our modeling on the outer part of this disk, consisting of a parent body ring and a halo of small grains. First, we confirm that a model using astronomical silicates cannot simultaneously fit the thermal and scattered light data. Next, we use a simple, generic function for the optical constants to show that varying the dust composition can improve the fit substantially. Finally, we model the dust as a mixture of the most plausible debris constituents: astronomical silicates, water ice, organic refractory material, and vacuum. We achieve a good fit to all datasets with grains composed predominantly of silicates and organics, while ice and vacuum are, at most, present in small amounts. This composition is similar to one derived from previous work on the HR 4796A disk. Our model also fits the thermal SED, scattered light colors, and high-resolution mid-IR data from T-ReCS for this disk. Additionally, we show that sub-blowout grains are a necessary component of the halo.

preprint2016arXiv

A New Star-Formation Rate Calibration from Polycyclic Aromatic Hydrocarbon Emission Features and Application to High Redshift Galaxies

We calibrate the integrated luminosity from the polycyclic aromatic hydrocarbon (PAH) features at 6.2\micron, 7.7\micron\ and 11.3\micron\ in galaxies as a measure of the star-formation rate (SFR). These features are strong (containing as much as 5-10\% of the total infrared luminosity) and suffer minimal extinction. Our calibration uses \spitzer\ Infrared Spectrograph (IRS) measurements of 105 galaxies at $0 < z < 0.4$, infrared (IR) luminosities of $10^9 - 10^{12} \lsol$, combined with other well-calibrated SFR indicators. The PAH luminosity correlates linearly with the SFR as measured by the extinction-corrected \ha\ luminosity over the range of luminosities in our calibration sample. The scatter is 0.14 dex comparable to that between SFRs derived from the \paa\ and extinction-corrected \ha\ emission lines, implying the PAH features may be as accurate a SFR indicator as hydrogen recombination lines. The PAH SFR relation depends on gas-phase metallicity, for which we supply an empirical correction for galaxies with $0.2 < \mathrm{Z} \lsim 0.7$~\zsol. We present a case study in advance of the \textit{James Webb Space Telescope} (\jwst), which will be capable of measuring SFRs from PAHs in distant galaxies at the peak of the SFR density in the universe ($z\sim2$) with SFRs as low as $\sim$~10~\sfrunits. We use \spitzer/IRS observations of the PAH features and \paa\ emission plus \ha\ measurements in lensed star-forming galaxies at $1 < z < 3$ to demonstrate the ability of the PAHs to derive accurate SFRs. We also demonstrate that because the PAH features dominate the mid-IR fluxes, broad-band mid-IR photometric measurements from \jwst\ will trace both the SFR and provide a way to exclude galaxies dominated by an AGN.

preprint2016arXiv

An Optical and Near-Infrared Study of the Type Ia/IIn Supernova PS15si

We present optical/near-infrared spectroscopy and photometry of the supernova (SN) PS15si. This object was originally identified as a Type IIn SN, but here we argue that it should be reclassified as a Type Ia SN with narrow hydrogen lines originating from interaction with circumstellar matter (CSM; i.e., SN Ia/IIn or SN Ia-CSM). Based on deep nondetections 27 days before discovery, we infer that this SN was discovered around or slightly before optical maximum, and we estimate the approximate time that it reached R-band maximum based on comparison with other SNe Ia/IIn. In terms of spectral morphology, we find that PS15si can be matched to a range of SN Ia spectral types, although SN 1991T-like SNe Ia provides the most self-consistent match. While this observation agrees with analysis of most other SNe Ia/IIn, we find that the implied CSM luminosity is too low to account for the overall luminosity of the SN at a time when the CSM should outshine the underlying SN by a few magnitudes. We infer that the similarity between PS15si and the hot, overluminous, high-ionisation spectrum of SN 1991T is a consequence of a spectrum that originates in ejecta layers that are heated by ultraviolet/X-ray radiation from CSM interaction. In addition, PS15si may have rebrightened over a short timescale in the B and V bands around 85 days after discovery, perhaps indicating that the SN ejecta are interacting with a local enhancement in CSM produced by clumps or a shell at large radii.

preprint2016arXiv

Distributions of quasar hosts on the galaxy main-sequence plane

The relation between star formation rates and stellar masses, i.e. the galaxy main sequence, is a useful diagnostic of galaxy evolution. We present the distributions relative to the main sequence of 55 optically-selected PG and 12 near-IR-selected 2MASS quasars at z <= 0.5. We estimate the quasar host stellar masses from Hubble Space Telescope or ground-based AO photometry, and the star formation rates through the mid-infrared aromatic features and far-IR photometry. We find that PG quasar hosts more or less follow the main sequence defined by normal star-forming galaxies while 2MASS quasar hosts lie systematically above the main sequence. PG and 2MASS quasars with higher nuclear luminosities seem to have higher specific SFRs (sSFRs), although there is a large scatter. No trends are seen between sSFRs and SMBH masses, Eddington ratios or even morphology types (ellipticals, spirals and mergers). Our results could be placed in an evolutionary scenario with quasars emerging during the transition from ULIRGs/mergers to ellipticals. However, combined with results at higher redshift, they suggest that quasars can be widely triggered in normal galaxies as long as they contain abundant gas and have ongoing star formation.

preprint2016arXiv

Photo-reverberation Mapping of a Protoplanetary Accretion Disk around a T Tauri Star

Theoretical models and spectroscopic observations of newborn stars suggest that protoplantary disks have an inner "wall" at a distance set by the disk interaction with the star. Around T Tauri stars, the size of this disk hole is expected to be on a 0.1-AU scale that is unresolved by current adaptive optics imaging, though some model-dependent constraints have been obtained by near-infrared interferometry. Here we report the first measurement of the inner disk wall around a solar-mass young stellar object, YLW 16B in the ρ Ophiuchi star-forming region, by detecting the light travel time of the variable radiation from the stellar surface to the disk. Consistent time lags were detected on two nights, when the time series in H (1.6 μm) and K (2.2 μm) bands were synchronized while the 4.5 μm emission lagged by 74.5 +/- 3.2 seconds. Considering the nearly edge-on geometry of the disk, the inner rim should be 0.084 AU from the protostar on average, with an error of order 0.01 AU. This size is likely larger than the range of magnetospheric truncations, and consistent with an optically and geometrically thick disk front at the dust sublimation radius at ~1500 K. The widths of the cross-correlation functions between the data in different wavebands place possible new constraints on the geometry of the disk.

preprint2016arXiv

The Correlation Between Metallicity and Debris Disk mass

We find that the initial dust masses in planetary debris disks are correlated with the metallicities of their central stars. We compiled a large sample of systems, including Spitzer, the Herschel DUNES and DEBRIS surveys, and WISE debris disk candidates. We also merged 33 metallicity catalogs to provide homogeneous [Fe/H] and $σ_{[Fe/H]}$ values. We analyzed this merged sample, including 222 detected disks (74 warm and 148 cold) around a total of 187 systems (some with multiple components) and 440 disks with only upper limits (125 warm and 315 cold), around a total of 360 systems. The disk dust masses at a common early evolutionary point in time were determined using our numerical disk evolutionary code, evolving a unique model for each of the 662 disks backward to an age of 1 Myr. We find that disk-bearing stars seldom have metallicities less than [Fe/H] = -0.2 and that the distribution of warm component masses lacks examples with large mass around stars of low metallicity ([Fe/H] < -0.085). Previous efforts to find a correlation have been largely unsuccessful; the primary improvements supporting our result are: 1.) basing the study on dust masses, not just infrared excess detections; 2.) including upper limits on dust mass in a quantitative way; 3.) accounting for the evolution of debris disk excesses as systems age; 4.) accounting fully for the range of uncertainties in metallicity measurements; and 5.) having a statistically large enough sample.

preprint2016arXiv

Variability in the Near-Infrared Synchrotron Emission From Cassiopeia A

We present multi-epoch Ks band imaging of the supernova remnant Cassiopeia A (Cas A). The morphology of the emission in this band is generally diffuse and filamentary, consistent with synchrotron radiation observed at radio wavelengths. However, in one region to the southwest of the remnant, compact knots of emission appear to be entrained in the ejecta and have the same proper motion as ejecta observed at similar projected radii. The presence of these knots suggests that material with high magnetic field strength contributes significantly to synchrotron emission at these wavelengths. We analyze these knots at J, H, and Ks bands as well as in 3.5-8 micron emission and at 6 cm where synchrotron emission is dominant and we find that the Ks band emission falls along the expected synchrotron spectrum. Using multi-epoch data, we calculate the magnetic field strength and electron density for a population of near-infrared synchrotron-emitting electrons. We find electron densities from 1,000-15,000 cm$^{-3}$ and magnetic field strengths from 1.3-5.8 mG. These magnetic field strengths are an order of magnitude higher than inferred from the much lower angular resolution gamma-ray observations toward Cas A. We conclude that dense knots of post-shock material behind the Cas A shock front are emitting synchrotron emission in a compressed and enhanced magnetic field.

preprint2015arXiv

A Systematic Survey for Broadened CO Emission Toward Galactic Supernova Remnants

We present molecular spectroscopy toward 50 Galactic supernova remnants (SNRs) taken at millimeter wavelengths in 12CO and 13CO J=2-1 with the Heinrich Hertz Submillimeter Telescope as part of a systematic survey for broad molecular line (BML) regions indicative of interactions with molecular clouds (MCs). These observations reveal BML regions toward nineteen SNRs, including nine newly identified BML regions associated with SNRs (G08.3-0.0, G09.9-0.8, G11.2-0.3, G12.2+0.3, G18.6-0.2, G23.6+0.3, 4C-04.71, G29.6+0.1, G32.4+0.1). The remaining ten SNRs with BML regions confirm previous evidence for MC interaction in most cases (G16.7+0.1, Kes 75, 3C 391, Kes 79, 3C 396, 3C 397, W49B, Cas A, IC 443), although we confirm that the BML region toward HB 3 is associated with the W3(OH) HII region, not the SNR. Based on the systemic velocity of each MC, molecular line diagnostics, and cloud morphology, we test whether these detections represent SNR-MC interactions. One of the targets (G54.1+0.3) had previous indications of a BML region, but we did not detect broadened emission toward it. Although broadened 12CO J=2-1 line emission should be detectable toward virtually all SNR-MC interactions we find relatively few examples; therefore, the number of interactions is low. This result favors mechanisms other than SN feedback as the basic trigger for star formation. In addition, we find no significant association between TeV gamma-ray sources and MC interactions, contrary to predictions that SNR-MC interfaces are the primary venues for cosmic ray acceleration.

preprint2015arXiv

Large Impacts around a Solar Analog Star in the Era of Terrestrial Planet Formation

The final assembly of terrestrial planets occurs via massive collisions, which can launch copious clouds of dust that are warmed by the star and glow in the infrared. We report the real-time detection of a debris-producing impact in the terrestrial planet zone around a 35-million year-old solar analog star. We observed a substantial brightening of the debris disk at 3-5 μm, followed by a decay over a year, with quasi-periodic modulations of the disk flux. The behavior is consistent with the occurrence of a violent impact that produced vapor out of which a thick cloud of silicate spherules condensed that were ground into dust by collisions. These results demonstrate how the time domain can become a new dimension for the study of terrestrial planet formation.

preprint2015arXiv

Magnetic grain trapping and the hot excesses around early-type stars

A significant fraction of main sequence stars observed interferometrically in the near infrared have slightly extended components that have been attributed to very hot dust. To match the spectrum appears to require the presence of large numbers of very small (< 200 nm in radius) dust grains. However, particularly for the hotter stars, it has been unclear how such grains can be retained close to the star against radiation pressure force. We find that the expected weak stellar magnetic fields are sufficient to trap nm-sized dust grains in epicyclic orbits for a few weeks or longer, sufficient to account for the hot excess emission. Our models provide a natural explanation for the requirement that the hot excess dust grains be smaller than 200 nm. They also suggest that magnetic trapping is more effective for rapidly rotating stars, consistent with the average vsini measurements of stars with hot excesses being larger (at about 2 sigma) than those for stars without such excesses.

preprint2015arXiv

Planetary Collisions outside the Solar System: Time Domain Characterization of Extreme Debris Disks

Luminous debris disks of warm dust in the terrestrial planet zones around solar-like stars are recently found to vary, indicative of ongoing large-scale collisions of rocky objects. We use Spitzer 3.6 and 4.5 μm time-series observations in 2012 and 2013 (extended to 2014 in one case) to monitor 5 more debris disks with unusually high fractional luminosities ("extreme debris disk"), including P1121 in the open cluster M47 (80 Myr), HD 15407A in the AB Dor moving group (80 Myr), HD 23514 in the Pleiades (120 Myr), HD 145263 in the Upper Sco Association (10 Myr), and the field star BD+20 307 (>1 Gyr). Together with the published results for ID8 in NGC 2547 (35 Myr), this makes the first systematic time-domain investigation of planetary impacts outside the solar system. Significant variations with timescales shorter than a year are detected in five out of the six extreme debris disks we have monitored. However, different systems show diverse sets of characteristics in the time domain, including long-term decay or growth, disk temperature variations, and possible periodicity.

preprint2015arXiv

Planetesimals in Debris Disks

Planetesimals form in gas-rich protoplanetary disks around young stars. However, protoplanetary disks fade in about 10 Myr. The planetesimals (and also many of the planets) left behind are too dim to study directly. Fortunately, collisions between planetesimals produce dusty debris disks. These debris disks trace the processes of terrestrial planet formation for 100 Myr and of exoplanetary system evolution out to 10 Gyr. This chapter begins with a summary of planetesimal formation as a prelude to the epoch of planetesimal destruction. Our review of debris disks covers the key issues, including dust production and dynamics, needed to understand the observations. Our discussion of extrasolar debris keeps an eye on similarities to and differences from Solar System dust.

preprint2015arXiv

The Inner Debris Structure in the Fomalhaut Planetary System

Fomalhaut plays an important role in the study of debris disks and small bodies in other planetary systems. The proximity and luminosity of the star make key features of its debris, like the water ice-line, accessible. Here we present ALMA cycle 1, 870 μm (345 GHz) observations targeted at the inner part of the Fomalhaut system with a synthesized beam of 0.45"x0.37" (~3 AU linear resolution at the distance of Fomalhaut) and a rms of 26 μJy/beam. The high angular resolution and sensitivity of the ALMA data enable us to place strong constraints on the nature of the warm excess revealed by Spitzer and Herschel observations. We detect a point source at the star position with a total flux consistent with thermal emission from the stellar photosphere. No structures that are brighter than 3σ are detected in the central 15 AU x 15 AU region. Modeling the spectral energy distribution using parameters expected for a dust-producing planetesimal belt indicates a radial location in the range ~8-15 AU. This is consistent with the location where ice sublimates in Fomalhaut, i.e., an asteroid-belt analog. The 3σ upper limit for such a belt is <1.3 mJy at 870 μm. We also interpret the 2 and 8-13 μm interferometric measurements to reveal the structure in the inner 10 AU region as dust naturally connected to this proposed asteroid belt by Poynting-Robertson drag, dust sublimation, and magnetically trapped nano grains.

preprint2015arXiv

The O- and B-Type Stellar Population in W3: Beyond the High-Density Layer

We present the first results from our survey of the star-forming complex W3, combining VRI photometry with multiobject spectroscopy to identify and characterize the high-mass stellar population across the region. With 79 new spectral classifications, we bring the total number of spectroscopically-confirmed O- and B-type stars in W3 to 105. We find that the high-mass slope of the mass function in W3 is consistent with a Salpeter IMF, and that the extinction toward the region is best characterized by an Rv of approximately 3.6. B-type stars are found to be more widely dispersed across the W3 giant molecular cloud (GMC) than previously realized: they are not confined to the high-density layer (HDL) created by the expansion of the neighboring W4 HII region into the GMC. This broader B-type population suggests that star formation in W3 began spontaneously up to 8--10 Myr ago, although at a lower level than the more recent star formation episodes in the HDL. In addition, we describe a method of optimizing sky subtraction for fiber spectra in regions of strong and spatially-variable nebular emission.

preprint2014arXiv

Debris Distribution in HD 95086 - A Young Analog of HR 8799

HD 95086 is a young early-type star that hosts (1) a 5 MJ planet at the projected distance of 56 AU revealed by direct imaging, and (2) a prominent debris disk. Here we report the detection of 69 um crystalline olivine feature from the disk using the Spitzer/MIPS-SED data covering 55-95 um. Due to the low resolution of MIPS-SED mode, this feature is not spectrally resolved, but is consistent with the emission from crystalline forsterite contributing 5% of the total dust mass. We also present detailed analysis of the disk SED and re-analysis of resolved images obtained by Herschel. Our results suggest that the debris structure around HD 95086 consists of a warm (175 K) belt, a cold (55 K) disk, and an extended disk halo (up to 800 AU), and is very similar to that of HR 8799. We compare the properties of the three debris components, and suggest that HD 95086 is a young analog of HR 8799. We further investigate and constrain single-planet, two-planet, three-planet and four-planet architectures that can account for the observed debris structure and are compatible with dynamical stability constraints. We find that equal-mass four-planet configurations of geometrically spaced orbits, with each planet of mass 5 MJ, could maintain the gap between the warm and cold debris belts, and also be just marginally stable for timescales comparable to the age of the system.

preprint2014arXiv

Exo-zodi modelling for the Large Binocular Telescope Interferometer

Habitable zone dust levels are a key unknown that must be understood to ensure the success of future space missions to image Earth analogues around nearby stars. Current detection limits are several orders of magnitude above the level of the Solar System's Zodiacal cloud, so characterisation of the brightness distribution of exo-zodi down to much fainter levels is needed. To this end, the large Binocular Telescope Interferometer (LBTI) will detect thermal emission from habitable zone exo-zodi a few times brighter than Solar System levels. Here we present a modelling framework for interpreting LBTI observations, which yields dust levels from detections and upper limits that are then converted into predictions and upper limits for the scattered light surface brightness. We apply this model to the HOSTS survey sample of nearby stars; assuming a null depth uncertainty of 10$^{-4}$ the LBTI will be sensitive to dust a few times above the Solar System level around Sun-like stars, and to even lower dust levels for more massive stars.

preprint2014arXiv

Interaction between Cassiopeia A and Nearby Molecular Clouds

We present spectroscopy of the supernova remnant Cassiopeia A (Cas A) observed at infrared wavelengths from 10 - 40 micron with the Spitzer Space Telescope and at millimeter wavelengths in 12CO and 13CO J = 2 - 1 (230 and 220 GHz) with the Heinrich Hertz Submillimeter Telescope. The IR spectra demonstrate high-velocity features toward a molecular cloud coincident with a region of bright radio continuum emission along the northern shock front of Cas A. The millimeter observations indicate that CO emission is broadened by a factor of two in some clouds toward Cas A, particularly to the south and west. We believe that these features trace interactions between the Cas A shock front and nearby molecular clouds. In addition, some of the molecular clouds that exhibit broadening in CO lie 1 - 2' away from the furthest extent of the supernova remnant shock front. We propose that this material may be accelerated by ejecta with velocity significantly larger than the observed free-expansion velocity of the Cas A shock front. These observations may trace cloud interactions with fast-moving outflows such as the bipolar outflow along the southwest to northeast axis of the Cas A supernova remnant, as well as fast-moving knots seen emerging in other directions.

preprint2014arXiv

Probing the Terrestrial Regions of Planetary Systems: Warm Debris Disks with Emission Features

Observations of debris disks allow for the study of planetary systems, even where planets have not been detected. However, debris disks are often only characterized by unresolved infrared excesses that resemble featureless blackbodies, and the location of the emitting dust is uncertain due to a degeneracy with the dust grain properties. Here we characterize the Spitzer IRS spectra of 22 debris disks exhibiting 10 micron silicate emission features. Such features arise from small warm dust grains, and their presence can significantly constrain the orbital location of the emitting debris. We find that these features can be explained by the presence of an additional dust component in the terrestrial zones of the planetary systems, i.e. an exozodiacal belt. Aside from possessing exozodiacal dust, these debris disks are not particularly unique; their minimum grain sizes are consistent with the blowout sizes of their systems, and their brightnesses are comparable to those of featureless warm debris disks. These disks are in systems with a range of ages, although the older systems with features are found only around A-type stars. The features in young systems may be signatures of terrestrial planet formation. Analyzing the spectra of unresolved debris disks with emission features may be one of the simplest and most accessible ways to study the terrestrial regions of planetary systems.

preprint2014arXiv

The Herschel Cold Debris Disks: Confusion with the Extragalactic Background at 160 mu

The Herschel "DUst around NEarby Stars (DUNES)" survey has found a number of debris disk candidates that are apparently very cold, with temperatures near 22K. It has proven difficult to fit their spectral energy distributions with conventional models for debris disks. Given this issue we carefully examine the alternative explanation, that the detections arise from confusion with IR cirrus and/or background galaxies that are not physically associated with the foreground star. We find that such an explanation is consistent with all of these detections.

preprint2014arXiv

Tracing Ram-Pressure Stripping with Warm Molecular Hydrogen Emission

We use the Spitzer Infrared Spectrograph (IRS) to study four infalling cluster galaxies with signatures of on-going ram-pressure stripping. H$_2$ emission is detected in all four; two show extraplanar H$_2$ emission. The emission usually has a warm (T $\sim$ $115 - 160$K) and a hot (T $\sim$ 400 $-$ 600K) component that is approximately two orders of magnitude less massive than the warm one. The warm component column densities are typically $10^{19} - 10^{20}$ cm$^{-2}$ with masses of $10^6 - 10^8 M_\odot$. The warm H$_2$ is anomalously bright compared with normal star-forming galaxies and therefore may be excited by ram-pressure. In the case of CGCG 97-073, the H$_2$ is offset from the majority of star formation along the direction of the galaxy's motion in the cluster, suggesting it is forming in the ram-pressure wake of the galaxy. Another galaxy, NGC 4522, exhibits a warm H$_2$ tail approximately 4 kpc in length. These results support the hypothesis that H$_2$ within these galaxies is shock-heated from the interaction with the intracluster medium. Stripping of dust is also a common feature of the galaxies. For NGC 4522, where the distribution of dust at 8 $μ$m is well resolved, knots and ripples demonstrate the turbulent nature of the stripping process. The H$α$ and 24 $μ$m luminosities show that most of the galaxies have star formation rates comparable to similar mass counterparts in the field. Finally, we suggest a possible evolutionary sequence primarily related to the strength of ram-pressure a galaxy experiences to explain the varied results observed in our sample.

preprint2013arXiv

A Trend Between Cold Debris Disk Temperature and Stellar Type: Implications for the Formation and Evolution of Wide-Orbit Planets

Cold debris disks trace the limits of planet formation or migration in the outer regions of planetary systems, and thus have the potential to answer many of the outstanding questions in wide-orbit planet formation and evolution. We characterized the infrared excess spectral energy distributions of 174 cold debris disks around 546 main-sequence stars observed by both Spitzer IRS and MIPS. We found a trend between the temperature of the inner edges of cold debris disks and the stellar type of the stars they orbit. This argues against the importance of strictly temperature-dependent processes (e.g. non-water ice lines) in setting the dimensions of cold debris disks. Also, we found no evidence that delayed stirring causes the trend. The trend may result from outward planet migration that traces the extent of the primordial protoplanetary disk, or it may result from planet formation that halts at an orbital radius limited by the efficiency of core accretion.

preprint2013arXiv

Local Luminous Infrared Galaxies. III. Co-evolution of Black Hole Growth and Star Formation Activity?

Local luminous infrared (IR) galaxies (LIRGs) have both high star formation rates (SFR) and a high AGN (Seyfert and AGN/starburst composite) incidence. Therefore, they are ideal candidates to explore the co-evolution of black hole (BH) growth and star formation (SF) activity, not necessarily associated with major mergers. Here, we use Spitzer/IRS spectroscopy of a complete volume-limited sample of local LIRGs (distances of <78Mpc). We estimate typical BH masses of 3x10^7 M_sun using [NeIII]15.56micron and optical [OIII]5007A gas velocity dispersions and literature stellar velocity dispersions. We find that in a large fraction of local LIRGs the current SFR is taking place not only in the inner nuclear ~1.5kpc region, as estimated from the nuclear 11.3micron PAH luminosities, but also in the host galaxy. We next use the ratios between the SFRs and BH accretion rates (BHAR) to study whether the SF activity and BH growth are contemporaneous in local LIRGs. On average, local LIRGs have SFR to BHAR ratios higher than those of optically selected Seyferts of similar AGN luminosities. However, the majority of the IR-bright galaxies in the RSA Seyfert sample behave like local LIRGs. Moreover, the AGN incidence tends to be higher in local LIRGs with the lowest SFRs. All this suggests that in local LIRGs there is a distinct IR-bright star forming phase taking place prior to the bulk of the current BH growth (i.e., AGN phase). The latter is reflected first as a composite and then as a Seyfert, and later as a non-LIRG optically identified Seyfert nucleus with moderate SF in its host galaxy.

preprint2013arXiv

Signpost of Multiple Planets in Debris Disks

We review the nearby debris disk structures revealed by multi-wavelength images from Spitzer and Herschel, and complemented with detailed spectral energy distribution modeling. Similar to the definition of habitable zones around stars, debris disk structures should be identified and characterized in terms of dust temperatures rather than physical distances so that the heating power of different spectral type of stars is taken into account and common features in disks can be discussed and compared directly. Common features, such as warm (~150 K) dust belts near the water-ice line and cold (~50 K) Kuiper-belt analogs, give rise to our emerging understanding of the levels of order in debris disk structures and illuminate various processes about the formation and evolution of exoplanetary systems. In light of the disk structures in the debris disk twins (Vega and Fomalhaut), and the current limits on the masses of planetary objects, we suggest that the large gap between the warm and cold dust belts is the best signpost for multiple (low-mass) planets beyond the water-ice line.

preprint2013arXiv

Spitzer Spectroscopy of Infrared-Luminous Galaxies: Diagnostics of AGN and Star Formation and Contribution to Total Infrared Luminosity

We use Spitzer IRS spectroscopy to study the nature 65 IR-luminous galaxies at 0.02 < z < 0.6 with F(24micron) > 1.2mJy. The IRS spectra cover wavelengths spanning the PAH features and important atomic diagnostic lines. Our sample corresponds to L(IR) = L(8-1000micron) = 10^10-10^12 Lsun. We divide our galaxies into those with Spitzer IRAC colors indicative of warm dust heated by an AGN (IRAGN) and those whose colors indicate star-formation processes (non-IRAGN). Compared to the non-IRAGN, the IRAGN show smaller PAH emission EWs, which we attribute to an increase in mid-IR continuum from the AGN. We find that in both the IRAGN and non-IRAGN samples, the PAH luminosities correlate strongly with the [Ne II] emission, from which we conclude that the PAH luminosity directly traces the instantaneous SFR in both the IRAGN and non-IRAGN galaxies. We compare the ratio of PAH luminosity to the total IR luminosity and show that for most IRAGN star-formation accounts for 10-50% of the L(IR). We also find no measurable difference between the PAH luminosity ratios of L(11.3)/L(7.7) and L(6.2)/L(7.7) for the IRAGN and non-IRAGN, suggesting that AGN do not significantly excite or destroy PAH molecules on galaxy-wide scales. A small subset of galaxies show excess of [O IV] emission compared to their PAH emission, which indicates the presence of heavily-obscured AGN, including 3 galaxies that are not otherwise selected as IRAGN. The low PAH emission and low [Ne II] emission of the IRAGN and [O IV]-excess objects imply they have low SFRs and their IR luminosity is dominated by processes associated with the AGN.

preprint2013arXiv

The Collisional Evolution of Debris Disks

We explore the collisional decay of disk mass and infrared emission in debris disks. With models, we show that the rate of the decay varies throughout the evolution of the disks, increasing its rate up to a certain point, which is followed by a leveling off to a slower value. The total disk mass falls off ~ t^-0.35 at its fastest point (where t is time) for our reference model, while the dust mass and its proxy -- the infrared excess emission -- fades significantly faster (~ t^-0.8). These later level off to a decay rate of M_tot(t) ~ t^-0.08 and M_dust(t) or L_ir(t) ~ t^-0.6. This is slower than the ~ t^-1 decay given for all three system parameters by traditional analytic models. We also compile an extensive catalog of Spitzer and Herschel 24, 70, and 100 micron observations. Assuming a log-normal distribution of initial disk masses, we generate model population decay curves for the fraction of debris disk harboring stars observed at 24 micron and also model the distribution of measured excesses at the far-IR wavelengths (70-100 micron) at certain age regimes. We show general agreement at 24 micron between the decay of our numerical collisional population synthesis model and observations up to a Gyr. We associate offsets above a Gyr to stochastic events in a few select systems. We cannot fit the decay in the far infrared convincingly with grain strength properties appropriate for silicates, but those of water ice give fits more consistent with the observations.

preprint2012arXiv

Modeling Collisional Cascades In Debris Disks: Steep Dust-Size Distributions

We explore the evolution of the mass distribution of dust in collision-dominated debris disks, using the collisional code introduced in our previous paper. We analyze the equilibrium distribution and its dependence on model parameters by evolving over 100 models to 10 Gyr. With our numerical models, we confirm that systems reach collisional equilibrium with a mass distribution that is steeper than the traditional solution by Dohnanyi (1969). Our model yields a quasi steady-state slope of n(m) ~ m^{-1.88} [n(a) ~ a^{-3.65}] as a robust solution for a wide range of possible model parameters. We also show that a simple power-law function can be an appropriate approximation for the mass distribution of particles in certain regimes. The steeper solution has observable effects in the submillimeter and millimeter wavelength regimes of the electromagnetic spectrum. We assemble data for nine debris disks that have been observed at these wavelengths and, using a simplified absorption efficiency model, show that the predicted slope of the particle mass distribution generates SEDs that are in agreement with the observed ones.

preprint2012arXiv

Modeling Collisional Cascades In Debris Disks: The Numerical Method

We develop a new numerical algorithm to model collisional cascades in debris disks. Because of the large dynamical range in particle masses, we solve the integro-differential equations describing erosive and catastrophic collisions in a particle-in-a-box approach, while treating the orbital dynamics of the particles in an approximate fashion. We employ a new scheme for describing erosive (cratering) collisions that yields a continuous set of outcomes as a function of colliding masses. We demonstrate the stability and convergence characteristics of our algorithm and compare it with other treatments. We show that incorporating the effects of erosive collisions results in a decay of the particle distribution that is significantly faster than with purely catastrophic collisions.

preprint2012arXiv

Resolved Debris Discs Around A Stars in the Herschel DEBRIS Survey

The majority of debris discs discovered so far have only been detected through infrared excess emission above stellar photospheres. While disc properties can be inferred from unresolved photometry alone under various assumptions for the physical properties of dust grains, there is a degeneracy between disc radius and dust temperature that depends on the grain size distribution and optical properties. By resolving the disc we can measure the actual location of the dust. The launch of Herschel, with an angular resolution superior to previous far-infrared telescopes, allows us to spatially resolve more discs and locate the dust directly. Here we present the nine resolved discs around A stars between 20 and 40 pc observed by the DEBRIS survey. We use these data to investigate the disc radii by fitting narrow ring models to images at 70, 100 and 160 μm and by fitting blackbodies to full spectral energy distributions. We do this with the aim of finding an improved way of estimating disc radii for unresolved systems. The ratio between the resolved and blackbody radii varies between 1 and 2.5. This ratio is inversely correlated with luminosity and any remaining discrepancies are most likely explained by differences to the minimum size of grain in the size distribution or differences in composition. We find that three of the systems are well fit by a narrow ring, two systems are borderline cases and the other four likely require wider or multiple rings to fully explain the observations, reflecting the diversity of planetary systems.

preprint2012arXiv

Variability of the Infrared Excess of Extreme Debris Disks

Debris disks with extremely large infrared excesses (fractional luminosities $> 10^{-2}$) are rare. Those with ages between 30 and 130 Myr are of interest because their evolution has progressed well beyond that of protoplanetary disks (which dissipate with a timescale of order 3 Myr), yet they represent a period when dynamical models suggest that terrestrial planet building may still be progressing through large, violent collisions that could yield large amounts of debris and large infrared excesses. For example, our Moon was formed through a violent collision of two large proto-planets during this age range. We report two disks around the solar-like stars ID8 and HD 23514 in this age range where the 24 μm infrared excesses vary on timescales of a few years, even though the stars are not variable in the optical. Variations this rapid are difficult to understand if the debris is produced by collisional cascades, as it is for most debris disks. It is possible that the debris in these two systems arises in part from condensates from silicate-rich vapor produced in a series of violent collisions among relatively large bodies. If their evolution is rapid, the rate of detection of extreme excesses would indicate that major collisions may be relatively common in this age range.

preprint2011arXiv

Local Luminous Infrared Galaxies. II. AGN Activity from Spitzer/IRS spectra

We quantify the active galactic nucleus (AGN) contribution to the mid-infrared (mid-IR) and the total infrared (IR, 8-1000micron) emission in a complete volume-limited sample of 53 local luminous infrared galaxies (LIRGs). We decompose the Spitzer Infrared Spectrograph (IRS) low-resolution 5-38micron spectra of the LIRGs into AGN and starburst components using clumpy torus models and star-forming galaxy templates, respectively. We find that 50% (25/50) of local LIRGs have an AGN component detected with this method. There is good agreement between these AGN detections through mid-IR spectral decomposition and other AGN indicators, such as the optical spectral class, mid-IR spectral features and X-ray properties. Taking all the AGN indicators together, the AGN detection rate in the individual nuclei of LIRGs is ~62%. The derived AGN bolometric luminosities are in the range L_bol(AGN)=0.4 -50x10^{43} erg/s. The AGN bolometric contribution to the IR luminosities of the galaxies is generally small, with 70% of LIRGs having L_bol(AGN)/L_IR<0.05. Only ~8% of local LIRGs have a significant AGN bolometric contribution L_bol(AGN)/L_IR > 0.25. From the comparison of our results with literature results of ultraluminous infrared galaxies, we confirm that in the local universe the AGN bolometric contribution to the IR luminosity increases with the IR luminosity of the galaxy/system. If we add up the AGN bolometric luminosities we find that AGNs only account for 5%^{+8%}_{-3%} of the total IR luminosity produced by local LIRGs (with and without AGN detections). This proves that the bulk of the IR luminosity of local LIRGs is due to star formation activity. Taking the newly determined IR luminosity density of LIRGs in the local universe, we then estimate an AGN IR luminosity density of Omega_IR(AGN) = 3x10^5 L_sun Mpc^{-3}$ in LIRGs.

preprint2011arXiv

SDSS-III: Massive Spectroscopic Surveys of the Distant Universe, the Milky Way Galaxy, and Extra-Solar Planetary Systems

Building on the legacy of the Sloan Digital Sky Survey (SDSS-I and II), SDSS-III is a program of four spectroscopic surveys on three scientific themes: dark energy and cosmological parameters, the history and structure of the Milky Way, and the population of giant planets around other stars. In keeping with SDSS tradition, SDSS-III will provide regular public releases of all its data, beginning with SDSS DR8 (which occurred in Jan 2011). This paper presents an overview of the four SDSS-III surveys. BOSS will measure redshifts of 1.5 million massive galaxies and Lya forest spectra of 150,000 quasars, using the BAO feature of large scale structure to obtain percent-level determinations of the distance scale and Hubble expansion rate at z<0.7 and at z~2.5. SEGUE-2, which is now completed, measured medium-resolution (R=1800) optical spectra of 118,000 stars in a variety of target categories, probing chemical evolution, stellar kinematics and substructure, and the mass profile of the dark matter halo from the solar neighborhood to distances of 100 kpc. APOGEE will obtain high-resolution (R~30,000), high signal-to-noise (S/N>100 per resolution element), H-band (1.51-1.70 micron) spectra of 10^5 evolved, late-type stars, measuring separate abundances for ~15 elements per star and creating the first high-precision spectroscopic survey of all Galactic stellar populations (bulge, bar, disks, halo) with a uniform set of stellar tracers and spectral diagnostics. MARVELS will monitor radial velocities of more than 8000 FGK stars with the sensitivity and cadence (10-40 m/s, ~24 visits per star) needed to detect giant planets with periods up to two years, providing an unprecedented data set for understanding the formation and dynamical evolution of giant planet systems. (Abridged)

preprint2011arXiv

Spitzer IRS Spectroscopy of the 10 Myr-old EF Cha Debris Disk: Evidence for Phyllosilicate-Rich Dust in the Terrestrial Zone

We describe Spitzer IRS spectroscopic observations of the 10 Myr-old star, EF Cha. Compositional modeling of the spectra from 5 μm to 35 μm confirms that it is surrounded by a luminous debris disk with LD/L\star ~ 10-3, containing dust with temperatures between 225 K and 430 K characteristic of the terrestrial zone. The EF Cha spectrum shows evidence for many solid-state features, unlike most cold, low-luminosity debris disks but like some other 10-20 Myr-old luminous, warm debris disks (e.g. HD 113766A). The EF Cha debris disk is unusually rich in a species or combination of species whose emissivities resemble that of finely powdered, laboratory-measured phyllosilicate species (talc, saponite, and smectite), which are likely produced by aqueous alteration of primordial anhydrous rocky materials. The dust and, by inference, the parent bodies of the debris also contain abundant amorphous silicates and metal sulfides, and possibly water ice. The dust's total olivine to pyroxene ratio of ~ 2 also provides evidence of aqueous alteration. The large mass volume of grains with sizes comparable to or below the radiation blow-out limit implies that planetesimals may be colliding at a rate high enough to yield the emitting dust but not so high as to devolatize the planetesimals via impact processing. Because phyllosilicates are produced by the interactions between anhydrous rock and warm, reactive water, EF Cha's disk is a likely signpost for water delivery to the terrestrial zone of a young planetary system.

preprint2011arXiv

The Relationship Between Black Hole Growth and Star Formation in Seyfert Galaxies

We present estimates of black hole accretion rates and nuclear, extended, and total star-formation rates for a complete sample of Seyfert galaxies. Using data from the Spitzer Space Telescope, we measure the active galactic nucleus (AGN) luminosity using the [O IV] 25.89 micron emission line and the star-forming luminosity using the 11.3 micron aromatic feature and extended 24 micron continuum emission. We find that black hole growth is strongly correlated with nuclear (r<1 kpc) star formation, but only weakly correlated with extended (r>1 kpc) star formation in the host galaxy. In particular, the nuclear star-formation rate (SFR) traced by the 11.3 micron aromatic feature follows a relationship with the black hole accretion rate (BHAR) of the form SFR\proptoBHAR^0.8, with an observed scatter of 0.5 dex. This SFR-BHAR relationship persists when additional star formation in physically matched r=1 kpc apertures is included, taking the form SFR\proptoBHAR^0.6. However, the relationship becomes almost indiscernible when total SFRs are considered. This suggests a physical connection between the gas on sub-kpc and sub-pc scales in local Seyfert galaxies that is not related to external processes in the host galaxy. It also suggests that the observed scaling between star formation and black hole growth for samples of AGNs will depend on whether the star formation is dominated by a nuclear or extended component. We estimate the integrated black hole and bulge growth that occurs in these galaxies and find that an AGN duty cycle of 5-10% would maintain the ratio between black hole and bulge masses seen in the local universe.

preprint2011arXiv

The X-ray emission of local luminous infrared galaxies

We study the X-ray emission of a representative sample of 27 local luminous infrared galaxies (LIRGs). The median IR luminosity of our sample is log L_IR/L_sun = 11.2, thus the low-luminosity end of the LIRG class is well represented. We used new XMM-Newton data as well as Chandra and XMM-Newton archive data. The soft X-ray (0.5-2 keV) emission of most of the galaxies (>80%), including LIRGs hosting a Seyfert 2 nucleus, is dominated by star-formation related processes. These LIRGs follow the star-formation rate (SFR) versus soft X-ray luminosity correlation observed in local starbursts. We find that ~15% of the non-Seyfert LIRGs (3 out of 20) have an excess hard X-ray emission relative to that expected from star-formation that might indicate the presence of an obscured AGN. The rest of the non-Seyfert LIRGs follow the SFR versus hard X-ray (2-10 keV) luminosity correlation of local starbursts. The non-detection of the 6.4 keV Fe K alpha emission line in the non-Seyfert LIRGs allows us to put an upper limit to the bolometric luminosity of an obscured AGN, L_bol <1043 erg s-1 . That is, in these galaxies, if they hosted a low luminosity AGN, its contribution to total luminosity would be less than 10%. Finally we estimate that the AGN contribution to the total luminosity for our sample of local LIRGs is between 7% and 10%.

preprint2010arXiv

A Warm Molecular Hydrogen Tail Due to Ram Pressure Stripping of a Cluster Galaxy

We have discovered a remarkable warm (130-160 K) molecular hydrogen tail with a H_2 mass of approximately 4*10^7 solar masses extending 20 kpc from a cluster spiral galaxy, ESO 137-001, in Abell 3627. At least half of this gas is lost permanently to the intracluster medium, as the tail extends beyond the tidal radius of the galaxy. We also detect a hot (400-550 K) component in the tail that is approximately 1% of the mass. The large H_2 line to IR continuum luminosity ratio in the tail indicates that star formation is not a major excitation source and that the gas is possibly shock-heated. This discovery confirms that the galaxy is currently undergoing ram-pressure stripping, as also indicated by its X-ray and Halpha tails found previously. We estimate the galaxy is losing its warm H_2 gas at a rate of ~ 2-3 solar masses per year. The true mass loss rate is likely higher if we account for cold molecular gas and atomic gas. We predict that the galaxy will lose most of its gas in a single pass through the core and place a strong upper limit on the ram-pressure timescale of 1 Gyr. We also study the star-forming properties of the galaxy and its tail. We identify most of the previously discovered external Halpha sources within the tail in our 8 um data but not in our 3.6 um data; IRS spectroscopy of the region containing these Halpha sources also reveals aromatic features typically associated with star formation. From the positions of these HII regions, it appears that star formation is not occurring throughout the molecular hydrogen tail but only immediately downstream of the galaxy. Some of these HII regions lie outside the tidal radius of the galaxy, indicating that ram-pressure stripping can be a source of intracluster stars.

preprint2010arXiv

An X-ray Study of Local Infrared Bright Galaxies

We are carrying out detailed study of the X-ray and infrared (IR) properties of a sample of local (d < 70 Mpc) luminous infrared galaxies (LIRGs) using XMM-Newton and Spitzer (imaging and spectroscopy). The main goal is to study the extreme processes of star formation and/or active galactic nuclei (AGN) taking place in this cosmologically important class of galaxies. In this proceedings we present the preliminary results obtained from the analysis of the XMM-Newton X-ray images and the X-ray spectral modeling.

preprint2010arXiv

Could the planets around HR 8799 be brown dwarfs?

We consider the limiting case for orbital stability of the companions to HR 8799. This case is only consistent with ages for the system of ~100 Myr, not with the 1 Gyr age proposed from astroseismology. The discrepancy probably arises because the inclination of the star is smaller than assumed in analyzing the astroseismology data. Given this young age, the best estimates of the companion masses place them by a small margin on the planet side of the division between planets and brown dwarfs.

preprint2010arXiv

Local Luminous Infrared Galaxies. I. Spatially resolved observations with Spitzer/IRS

We present results from the Spitzer/IRS spectral mapping observations of 15 local luminous infrared galaxies (LIRGs). In this paper we investigate the spatial variations of the mid-IR emission which includes: fine structure lines, molecular hydrogen lines, polycyclic aromatic features (PAHs), continuum emission and the 9.7um silicate feature. We also compare the nuclear and integrated spectra. We find that the star formation takes place in extended regions (several kpc) as probed by the PAH emission as well as the [NeII] and [NeIII] emissions. The behavior of the integrated PAH emission and 9.7um silicate feature is similar to that of local starburst galaxies. We also find that the minima of the [NeIII]/[NeII] ratio tends to be located at the nuclei and its value is lower than that of HII regions in our LIRGs and nearby galaxies. It is likely that increased densities in the nuclei of LIRGs are responsible for the smaller nuclear [NeIII]/[NeII] ratios. This includes the possibility that some of the most massive stars in the nuclei are still embedded in ultracompact HII regions. In a large fraction of our sample the 11.3um PAH emission appears more extended than the dust 5.5um continuum emission. We find a dependency of the 11.3um PAH/7.7 um PAH and [NeII]/11.3um PAH ratios with the age of the stellar populations. Smaller and larger ratios respectively indicate recent star formation. The estimated warm (300 K < T < 1000 K) molecular hydrogen masses are of the order of 10^8 M_Sun, which are similar to those found in ULIRGs, local starbursts and Seyfert galaxies. Finally we find that the [NeII] velocity fields for most of the LIRGs in our sample are compatible with a rotating disk at ~kpc scales, and they are in a good agreement with H-alpha velocity fields.

preprint2010arXiv

Locating the planetesimals belts in the multiple-planet systems HD 128311, HD 202206, HD 82943 and HR 8799

In addition to the Sun, six other stars are known to harbor multiple planets and debris disks: HD 69830, HD 38529, HD 128311, HD 202206, HD 82943 and HR 8799. In this paper we set constraints on the location of the dust-producing planetesimals around the latter four systems. We use a radiative transfer model to analyze the spectral energy distributions of the dust disks (including two new Spitzer IRS spectra presented in this paper), and a dynamical model to assess the long-term stability of the planetesimals' orbits. As members of a small group of stars that show evidence of harboring a multiple planets and planetesimals, their study can help us learn about the diversity of planetary systems.

preprint2010arXiv

Morphology and Size Differences between Local & High Redshift Luminous Infrared Galaxies

We show that the star-forming regions in high-redshift luminous and ultraluminous infrared galaxies (LIRGs and ULIRGs) and submillimeter galaxies (SMGs) have similar physical scales to those in local normal star-forming galaxies. To first order, their higher infrared (IR) luminosities result from higher luminosity surface density. We also find a good correlation between the IR luminosity and IR luminosity surface density in starburst galaxies across over five orders of magnitude of IR luminosity from local normal galaxies to z ~ 2 SMGs. The intensely star-forming regions of local ULIRGs are significantly smaller than those in their high-redshift counterparts and hence diverge significantly from this correlation, indicating that the ULIRGs found locally are a different population from the high-redshift ULIRGs and SMGs. Based on this relationship, we suggest that luminosity surface density should serve as a more accurate indicator for the IR emitting environment, and hence the observable properties, of star-forming galaxies than their IR luminosity. We demonstrate this approach by showing that ULIRGs at z ~ 1 and a lensed galaxy at z ~ 2.5 exhibit aromatic features agreeing with local LIRGs that are an order of magnitude less luminous, but have similar IR luminosity surface density. A consequence of this relationship is that the aromatic emission strength in star-forming galaxies will appear to increase at z > 1 for a given IR luminosity compared to their local counterparts.

preprint2010arXiv

The Effect of Active Galactic Nuclei on the Mid-Infrared Aromatic Features

We present Spitzer measurements of the aromatic (also known as PAH) features for 35 Seyfert galaxies from the revised Shapley-Ames sample and find that the relative strengths of the features differ significantly from those observed in star-forming galaxies. Specifically, the features at 6.2, 7.7, and 8.6 micron are suppressed relative to the 11.3 micron feature in Seyferts. Furthermore, we find an anti-correlation between the L(7.7 micron)/L(11.3 micron) ratio and the strength of the rotational H2 (molecular hydrogen) emission, which traces shocked gas. This suggests that shocks suppress the short-wavelength features by modifying the structure of the aromatic molecules or destroying the smallest grains. Most Seyfert nuclei fall on the relationship between aromatic emission and [Ne II] emission for star-forming galaxies, indicating that aromatic-based estimates of the star-formation rate in AGN host galaxies are generally reasonable. For the outliers from this relationship, which have small L(7.7 micron)/L(11.3 micron) ratios and strong H2 emission, the 11.3 micron feature still provides a valid measure of the star-formation rate.

preprint2010arXiv

The Evolution of the Star Formation Rate of Galaxies at 0.0 < z < 1.2

We present the 24 micron rest-frame luminosity function (LF) of star-forming galaxies in the redshift range 0.0 < z < 0.6 constructed from 4047 spectroscopic redshifts from the AGN and Galaxy Evolution Survey of 24 micron selected sources in the Bootes field of the NOAO Deep Wide-Field Survey. This sample provides the best available combination of large area (9 deg^2), depth, and statistically complete spectroscopic observations, allowing us to probe the evolution of the 24 micron LF of galaxies at low and intermediate redshifts while minimizing the effects of cosmic variance. In order to use the observed 24 micron luminosity as a tracer for star formation, active galactic nuclei (AGNs) that could contribute significantly at 24 micron are identified and excluded from our star-forming galaxy sample based on their mid-IR spectral energy distributions or the detection of X-ray emission. The evolution of the 24 micron LF of star-forming galaxies for redshifts of z < 0.65 is consistent with a pure luminosity evolution where the characteristic 24 micron luminosity evolves as (1+z)^(3.8+/-0.3). We extend our evolutionary study to encompass 0.0 < z < 1.2 by combining our data with that of the Far-Infrared Deep Extragalactic Legacy Survey. Over this entire redshift range the evolution of the characteristic 24 micron luminosity is described by a slightly shallower power law of (1+z)^(3.4+/-0.2). We find a local star formation rate density of (1.09+/-0.21) x 10^-2 Msun/yr/Mpc^-3, and that it evolves as (1+z)^(3.5+/-0.2) over 0.0 < z < 1.2. These estimates are in good agreement with the rates using optical and UV fluxes corrected for the effects of intrinsic extinction in the observed sources. This agreement confirms that star formation at z <~ 1.2 is robustly traced by 24 micron observations and that it largely occurs in obscured regions of galaxies. (Abridged)

preprint2010arXiv

The Mid-Infrared High-Ionization Lines from Active Galactic Nuclei and Star-Forming Galaxies

We used Spitzer/IRS spectroscopic data on 426 galaxies including quasars, Seyferts, LINER and HII galaxies to investigate the relationship among the mid-IR emission lines. There is a tight linear correlation between the [Ne V]14.3 um and 24.3 um (97.1 eV) and the [O IV]25.9 um (54.9 eV) high-ionization emission lines. The correlation also holds for these high-ionization emission lines and the [Ne III]15.56 um (41 eV) emission line, although only for active galaxies. We used these correlations to calculate the [Ne III] excess due to star formation in Seyfert galaxies. We also estimated the [O IV] luminosity due to star formation in active galaxies and determined that it dominates the [O IV] emission only if the contribution of the active nucleus to the total luminosity is below 5%. We find that the AGN dominates the [O IV] emission in most Seyfert galaxies, whereas star-formation adequately explains the observed [O IV] emission in optically classified HII galaxies. Finally we computed photoionization models to determine the physical conditions of the narrow line region where these high-ionization lines originate. The estimated ionization parameter range is -2.8 < log U < -2.5 and the total hydrogen column density range is 20 < log nH (cm-2) < 21.

preprint2010arXiv

The Structure of the β Leonis Debris Disk

We combine nulling interferometry at 10 μm using the MMT and Keck Telescopes with spectroscopy, imaging, and photometry from 3 to 100 μm using Spitzer to study the debris disk around β Leo over a broad range of spatial scales, corresponding to radii of 0.1 to ~100 AU. We have also measured the close binary star o Leo with both Keck and MMT interferometers to verify our procedures with these instruments. The β Leo debris system has a complex structure: 1.) relatively little material within 1 AU; 2.) an inner component with a color temperature of ~600 K, fitted by a dusty ring from about 2 to 3 AU; and 3.) a second component with a color temperature of ~120 K fitted by a broad dusty emission zone extending from about ~5 AU to ~55 AU. Unlike many other A-type stars with debris disks, β Leo lacks a dominant outer belt near 100 AU.

preprint2010arXiv

Unobscured Type 2 AGNs

Type 2 AGNs with intrinsically weak broad emission lines (BELs) would be exceptions to the unified model. After examining a number of proposed candidates critically, we find that the sample is contaminated significantly by objects with BELs of strengths indicating that they actually contain intermediate-type AGNs, plus a few Compton-thick sources as revealed by extremely low ratios of X-ray to nuclear IR luminosities. We develop quantitative metrics that show two (NGC 3147 and NGC 4594) of the remaining candidates to have BELs 2-3 orders of magnitude weaker than those of typical type-1 AGNs. Several more galaxies remain as candidates to have anomalously weak BELs, but this status cannot be confirmed with the existing information. Although the parent sample is poorly defined, the two confirmed objects are well under 1% of its total number of members, showing that the absence of a BEL is possible, but very uncommon in AGN. We evaluate these two objects in detail using multi-wavelength measurements. They have little X-ray extinction with N_H < 10^21 cm^{-2}. Their IR spectra show strong silicate emission (NGC 4594) or weak aromatic features on a generally power law continuum with a suggestion of silicates in emission (NGC 3147). No polarized BEL is detected in NGC 3147. These results indicate that the two unobscured type-2 objects have circumnuclear tori that are approximately face-on. Combined with their X-ray and optical/UV properties, this behavior implies that we have an unobscured view of the nuclei and thus that they have intrinsically weak BELs. We compare their properties with those of the other less-extreme candidates. We then compare the distributions of bolometric luminosities and accretion rates of these objects with theoretical models that predict weak BELs.

preprint2009arXiv

Cosmic Evolution of Star Formation In SDSS Quasar Hosts Since z=1

We present Spitzer IRS observations of a complete sample of 57 SDSS type-1 quasars at z~1. Aromatic features at 6.2 and/or 7.7 um are detected in about half of the sample and show profiles similar to those seen in normal galaxies at both low- and high-redshift, indicating a star-formation origin for the features. Based on the ratio of aromatic to star-formation IR (SFIR) luminosities for normal star-forming galaxies at z~1, we have constructed the SFIR luminosity function (LF) of z~1 quasars. As we found earlier for low-redshift PG quasars, these z~1 quasars show a flatter SFIR LF than do z~1 field galaxies, implying the quasar host galaxy population has on average a higher SFR than the field galaxies do. As measured from their SFIR LF, individual quasar hosts have on average LIRG-level SFRs, which mainly arise in the circumnuclear regions. By comparing with similar measurements of low-redshift PG quasars, we find that the comoving SFIR luminosity density in quasar hosts shows a much larger increase with redshift than that in field galaxies. The behavior is consistent with pure density evolution since the average SFR and the average SFR/BH-accretion-rate in quasar hosts show little evolution with redshift. For individual quasars, we have found a correlation between the aromatic-based SFR and the luminosity of the nuclear radiation, consistent with predictions of some theoretical models. We propose that type 1 quasars reside in a distinct galaxy population that shows elliptical morphology but that harbors a significant fraction of intermediate-age stars and is experiencing intense circumnuclear star formation.

preprint2009arXiv

Evidence for Dynamical Changes in a Transitional Protoplanetary Disk with Mid-infrared Variability

We present multi-epoch Spitzer Space Telescope observations of the transitional disk LRLL 31 in the 2-3 Myr-old star forming region IC 348. Our measurements show remarkable mid-infrared variability on timescales as short as one week. The infrared continuum emission exhibits systematic wavelength-dependent changes that suggest corresponding dynamical changes in the inner disk structure and variable shadowing of outer disk material. We propose several possible sources for the structural changes, including a variable accretion rate or a stellar or planetary companion embedded in the disk. Our results indicate that variability studies in the infrared can provide important new constraints on protoplanetary disk behavior.

preprint2009arXiv

Isotropic Luminosity Indicators in a Complete AGN Sample

The [O IV] 25.89 micron line has been shown to be an accurate indicator of active galactic nucleus (AGN) intrinsic luminosity in that it correlates well with hard (10-200 keV) X-ray emission. We present measurements of [O IV] for 89 Seyfert galaxies from the unbiased Revised Shapley-Ames (RSA) sample. The [O IV] luminosity distributions of obscured and unobscured Seyferts are indistinguishable, indicating that their intrinsic AGN luminosities are quite similar and that the RSA sample is well suited for tests of the unified model. In addition, we analyze several commonly used proxies for AGN luminosity, including [O III] 5007 A, 6 cm radio, and 2-10 keV X-ray emission. We find that the radio luminosity distributions of obscured and unobscured AGNs show no significant difference, indicating that radio luminosity is a useful isotropic luminosity indicator. However, the observed [O III] and 2-10 keV luminosities are systematically smaller for obscured Seyferts, indicating that they are not emitted isotropically.

preprint2009arXiv

Local Luminous Infrared Galaxies: Spatially resolved mid-infrared observations with Spitzer/IRS

Luminous Infrared (IR) Galaxies (LIRGs) are an important cosmological class of galaxies as they are the main contributors to the co-moving star formation rate density of the universe at z=1. In this paper we present a GTO Spitzer IRS program aimed to obtain spectral mapping of a sample of 14 local (d<76Mpc) LIRGs. The data cubes map, at least, the central 20arcsec x 20arcsec to 30arcsec x 30arcsec regions of the galaxies, and use all four IRS modules covering the full 5-38micron spectral range. The final goal of this project is to characterize fully the mid-IR properties of local LIRGs as a first step to understanding their more distant counterparts. In this paper we present the first results of this GTO program. The IRS spectral mapping data allow us to build spectral maps of the bright mid-IR emission lines (e.g., [NeII], [NeIII], [SIII], H_2), continuum, the 6.2 and 11.3micron PAH features, and the 9.7micron silicate feature, as well as to extract 1D spectra for regions of interest in each galaxy. The IRS data are used to obtain spatially resolved measurements of the extinction using the 9.7micron silicate feature, and to trace star forming regions using the neon lines and the PAH features. We also investigate a number of AGN indicators, including the presence of high excitation emission lines and a strong dust continuum emission at around 6micron. We finally use the integrated Spitzer/IRS spectra as templates of local LIRGs. We discuss several possible uses for these templates, including the calibration of the star formation rate of IR-bright galaxies at high redshift. We also predict the intensities of the brightest mid-IR emission lines for LIRGs as a function of redshift, and compare them with the expected sensitivities of future space IR missions.

preprint2009arXiv

Spitzer and HHT observations of starless cores: masses and environments

We present Spitzer observations of a sample of 12 starless cores selected to have prominent 24 micron shadows. The Spitzer images show 8 and 24 micron shadows and in some cases 70 micron shadows; these spatially resolved absorption features trace the densest regions of the cores. We have carried out a 12CO (2-1) and 13CO (2-1) mapping survey of these cores with the Heinrich Hertz Telescope (HHT). We use the shadow features to derive optical depth maps. We derive molecular masses for the cores and the surrounding environment; we find that the 24 micron shadow masses are always greater than or equal to the molecular masses derived in the same region, a discrepancy likely caused by CO freeze--out onto dust grains. We combine this sample with two additional cores that we studied previously to bring the total sample to 14 cores. Using a simple Jeans mass criterion we find that ~ 2/3 of the cores selected to have prominent 24 micron shadows are collapsing or near collapse, a result that is supported by millimeter line observations. Of this subset at least half have indications of 70 micron shadows. All cores observed to produce absorption features at 70 micron are close to collapse. We conclude that 24 micron shadows, and even more so the 70 micron ones, are useful markers of cloud cores that are approaching collapse.

preprint2009arXiv

The Extreme Star Formation Activity of Arp299 Revealed by Spitzer IRS Spectral Mapping

We present Spitzer/IRS spectral mapping observations of the luminous infrared galaxy (LIRG) Arp299 (IC694 + NGC3690) covering the central 45arcsec ~ 9kpc. The integrated mid-IR spectrum of Arp299 is similar to that of local starbursts despite its strongly interacting nature and high infrared luminosity, L_IR ~ 6x10^11 Lsun. This is explained because the star formation (probed by e.g. high [NeIII]15.56micron/[NeII]micron line ratios) is spread across at least 6-8kpc. Moreover, a large fraction of this star formation is taking place in young regions of moderate mid-IR optical depths such as the C+C' complex in the overlap region between the two galaxies and in HII regions in the disks of the galaxies. It is only source A, the nuclear region of IC694, that shows the typical mid-IR characteristics of ultraluminous infrared galaxies (ULIRGs, L_IR > 10^12 Lsun), that is, very compact (less than 1kpc) and dust-enshrouded star formation resulting in a deep silicate feature and moderate equivalent widths of the PAHs. The nuclear region of NGC3690, known as source B1, hosts a low-luminosity AGN and is surrounded by regions of star formation. Although the high excitation [NeV]14.32micron line typical of AGN is not detected in B1, its upper limit is consistent with the value expected from the X-ray luminosity. The AGN emission is detected in the form of a strong hot dust component that accounts for 80-90% of the 6micron luminosity of B1. The similarity between the Arp299 integrated mid-IR spectrum and those of high-z ULIRGs suggests that Arp299 may represent a local example, albeit with lower IR luminosity and possibly higher metallicity, of the star-formation processes occurring at high-z.

preprint2006arXiv

Warm Dust and Spatially Variable PAH Emission in the Dwarf Starburst Galaxy NGC 1705

We present Spitzer observations of the dwarf starburst galaxy NGC 1705 obtained as part of SINGS. The galaxy morphology is very different shortward and longward of ~5 microns: short-wavelength imaging shows an underlying red stellar population, with the central super star cluster (SSC) dominating the luminosity; longer-wavelength data reveals warm dust emission arising from two off-nuclear regions offset by ~250 pc from the SSC. These regions show little extinction at optical wavelengths. The galaxy has a relatively low global dust mass (~2E5 solar masses, implying a global dust-to-gas mass ratio ~2--4 times lower than the Milky Way average). The off-nuclear dust emission appears to be powered by photons from the same stellar population responsible for the excitation of the observed H Alpha emission; these photons are unassociated with the SSC (though a contribution from embedded sources to the IR luminosity of the off-nuclear regions cannot be ruled out). Low-resolution IRS spectroscopy shows moderate-strength PAH emission in the 11.3 micron band in the eastern peak; no PAH emission is detected in the SSC or the western dust emission complex. There is significant diffuse 8 micron emission after scaling and subtracting shorter wavelength data; the spatially variable PAH emission strengths revealed by the IRS data suggest caution in the interpretation of diffuse 8 micron emission as arising from PAH carriers alone. The metallicity of NGC 1705 falls at the transition level of 35% solar found by Engelbracht and collaborators; the fact that a system at this metallicity shows spatially variable PAH emission demonstrates the complexity of interpreting diffuse 8 micron emission. A radio continuum non-detection, NGC 1705 deviates significantly from the canonical far-IR vs. radio correlation. (Abridged)

preprint2005arXiv

Understanding Radio-Selected Thermal Sources in M 33: Ultraviolet, Optical, Near-Infrared, Spitzer Mid-Infrared, and Radio Observations

We present ultraviolet, optical, near-infrared, Spitzer mid-infrared, and radio images of 14 radio-selected objects in M 33. These objects are thought to represent the youngest phase of star cluster formation. We have detected the majority of cluster candidates in M 33 at all wavelengths. From the near-IR images, we derived ages 2-10 Myr, K_S-band extinctions (A_K_S) of 0-1 mag, and stellar masses of 10^3-10^4 M_solar. We have generated spectral energy distributions (SEDs) of each cluster from 0.1 micron to 160 microns. From these SEDs, we have modeled the dust emission around these star clusters to determine the dust masses (1-10^3 M_solar) and temperatures (40-90 K) of the clusters' local interstellar medium. Extinctions derived from the JHK_S, Halpha, and UV images are similar to within a factor of 2 or 3. These results suggest that eleven of the fourteen radio-selected objects are optically-visible young star clusters with a surrounding H II region, that two are background objects, possibly AGN, and that one is a Wolf-Rayet star with a surrounding H II region.

preprint2003arXiv

A Medium-Resolution Near-Infrared Spectral Library of Late Type Stars: I

We present an empirical infrared spectral library of medium resolution (R~2000-3000) H (1.6 micron) and K (2.2 micron) band spectra of 218 red stars, spanning a range of [Fe/H] from ~-2.2 to ~+0.3. The sample includes Galactic disk stars, bulge stars from Baade's window, and red giants from Galactic globular clusters. We report the values of 19 indices covering 12 spectral features measured from the spectra in the library. Finally, we derive calibrations to estimate the effective temperature, and diagnostic relationships to determine the luminosity classes of individual stars from near-infrared spectra. This paper is part of a larger effort aimed at building a near-IR spectral library to be incorporated in population synthesis models, as well as, at testing synthetic stellar spectra.