Source author record

William T. Reach

William T. Reach 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

30works
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

30 published item(s)

preprint2022arXiv

Extragalactic magnetism with SOFIA (SALSA Legacy Program) -- III: First data release and on-the-fly polarization mapping characterization

We describe the data processing of the Survey on extragALactic magnetiSm with SOFIA (SALSA Legacy Program). This first data release presents 33% (51.34h out of 155.7h, including overheads) of the total awarded time taken from January 2020 to December 2021. Our observations were performed using the newly implemented on-the-fly mapping (OTFMAP) technique in the polarimetric mode. We present the pipeline steps to obtain homogeneously reduced high-level data products of polarimetric maps of galaxies for use in scientific analysis. Our approach has a general design and can be applied to sources smaller than the field-of-view of the HAWC+ array in any given band. We estimate that the OTFMAP polarimetric mode offers a reduction of observing overheads by a factor 2.34, and an improvement in sensitivity by a factor 1.80 when compared to previously obtained polarimetric observations using the chopping and nodding mode. The OTFMAP is a significant optimization of the polarimetric mode of HAWC+ as it ultimately reduces the cost of operations of SOFIA/HAWC+ by increasing the science collected per hour of observation up to an overall factor of 2.49. The OTFMAP polarimetric mode is the standard observing strategy of SALSA. The results and quantitative analysis of this first data release are presented in Papers IV and V of the series.

preprint2022arXiv

Extragalactic magnetism with SOFIA (SALSA Legacy Program) -- IV: Program overview and first results on the polarization fraction

We present the first data release of the Survey on extragALactic magnetiSm with SOFIA (SALSA Legacy Program) with a set of 14 nearby ($<20$ Mpc) galaxies with resolved imaging polarimetric observations using HAWC+ from $53$ to $214$ $μ$m at a resolution of $5-18$" ($90$ pc $-$ $1$ kpc). We introduce the definitions and background on extragalactic magnetism, and present the scientific motivation and sample selection of the program. Here, we focus on the general trends in the emissive polarization fraction. Far-infrared polarimetric observations trace the thermal polarized emission of magnetically aligned dust grains across the galaxy disks with polarization fractions of $P=0-15$% in the cold, $T_{\rm d} = [19,48]$ K, and dense, $\log_{10}(N_{\rm HI+H_{2}}) = [19.96,22.91]$, interstellar medium. The spiral galaxies show a median $\langle P_{154μm} \rangle = 3.3\pm0.9 $% across the disks. We report the first polarized spectrum of starburst galaxies showing a minimum within $89-154$ $μ$m. The falling $53-154$ $μ$m polarized spectrum may be due to a decrease in the dust grain alignment efficiency produced by variations in dust temperatures along the line-of-sight in the galactic outflow. We find that the starburst galaxies and the star-forming regions within normal galaxies have the lowest polarization fractions. We find that 50% (7 out of 14) of the galaxies require a broken power-law in the $P-N_{HI+H_{2}}$ and $P-T_{d}$ relations with three different trends. Group 1 has a relative increase of anisotropic random B-fields produced by compression or shear of B-fields in the galactic outflows, starburst rings, and inner-bar of galaxies; and Groups 2 and 3 have a relative increase of isotropic random B-fields driven by star-forming regions in the spiral arms, and/or an increase of dust grain alignment efficiency caused by shock-driven regions or evolutionary stages of a galaxy.

preprint2022arXiv

Extragalactic magnetism with SOFIA (SALSA Legacy Program). VI. The magnetic fields in the multi-phase interstellar medium of the Antennae galaxies

Mergers are thought to be a fundamental channel for galaxy growth, perturbing the gas dynamics and the magnetic fields (B-fields) in the interstellar medium (ISM). However, the mechanisms that amplify and dissipate B-fields during a merger remain unclear. We characterize the morphology of the ordered B-fields in the multi-phase ISM of the closest merger of two spiral galaxies, the Antennae galaxies. We compare the inferred B-fields using $154~μ$m thermal dust and $11$ cm radio synchrotron emission polarimetric observations. We find that the $154~μ$m B-fields are more ordered across the Antennae galaxies than the $11$ cm B-fields. The turbulent-to-ordered $154~μ$m B-field increases at the galaxy cores and star-forming regions. The relic spiral arm has an ordered spiral $154~μ$m B-field, while the $11$ cm B-field is radial. The $154~μ$m B-field may be dominated by turbulent dynamos with high $^{12}$CO(1-0) velocity dispersion driven by star-forming regions, while the $11$ cm B-field is cospatial with high HI velocity dispersion driven by galaxy interaction. This result shows the dissociation between the warm gas mainly disturbed by the merger, and the dense gas still following the dynamics of the relic spiral arm. We find a $\sim8.9$ kpc scale ordered B-field connecting the two galaxies. The base of the tidal tail is cospatial with the HI and $^{12}$CO(1-0) emission and has compressed and/or sheared $154~μ$m and $11$ cm B-fields driven by the merger. We suggest that amplify B-fields, with respect to the rest of the system and other spiral galaxies, may be supporting the gas flow between both galaxies and the tidal tail.

preprint2021arXiv

Ionized carbon around IRC +10216

Asymptotic giant branch (AGB)stars create a rich inventory of molecules in their envelopes as they lose mass during later stages of their evolution. These molecules cannot survive the conditions in interstellar space, where they are exposed to ultraviolet photons of the interstellar radiation field. As a result, daughter molecules are the ones injected into space, and a halo of those molecules is predicted to exist around cool evolved stars. The most abundant molecule in the envelopes other than H2 is CO, which dissociates into C that is rapidly ionized in a halo around the star that is optically thin to the interstellar radiation field. We develop the specific predictions of the ionized carbon halo size and column density for the well-studied, nearby star IRC +10216. We compare those models to observations of the [C II] 157.7 micron far-infrared fine-structure line using SOFIA and Herschel. The combination of bright emission toward the star and upper limits to extended [C II] is inconsistent with any standard model. The presence of [C II] toward the star requires some dissociation and ionization in the inner part of the outflow, possibly due to a hot companion star. The lack of extended [C II] emission requires that daughter products from CO photodissociation in the outer envelope remain cold. The [C II] profile toward the star is asymmetric, with the blue-shifted absorption due to the cold outer envelope.

preprint2020arXiv

Observational evidence for rotational desorption of Complex Molecules by radiative torques from Orion BN/KL

Complex Organic Molecules (COMs) are believed to form in the ice mantle of dust grains and are released to the gas by thermal sublimation when grain mantles are heated to temperatures of $T_{\rm d}\gtrsim 100\,\rm K$. However, some COMs are detected in regions with temperatures below 100 K. Recently, a new mechanism of rotational desorption due to centrifugal stress induced by radiative torques (RATs) is proposed by Hoang & Tram 2020 that can desorb COMs at low temperatures. In this paper, we report observational evidence for rotational desorption of COMs toward the nearest massive star-forming region Orion BN/KL. We compare the abundance of three representative COMs which have very high binding energy computed by the rotational desorption mechanism with observations by ALMA, and demonstrate that the rotational desorption mechanism can explain the existence of such COMs. We also analyze the polarization data from SOFIA/HAWC+ and JCMT/SCUBA-2 and find that the polarization degree at far-infrared/submm decreases with increasing the grain temperature for $T_{\rm d}\gtrsim 71\,\rm K$. This is consistent with the theoretical prediction using the Radiative Torque (RAT) alignment theory and Radiative Torque Disruption (RATD) mechanism. Such an anti-correlation between dust polarization and dust temperature supports the rotational disruption as well as rotational desorption mechanism of COMs induced by RATs.

preprint2016arXiv

Mid-infrared spectra of comet nuclei

Jovian Trojan D-type asteroids have mid-infrared emissivity features strikingly similar to comet comae, suggesting that they have the same compositions and that the surfaces of the Trojans are highly porous. However, a direct comparison between a comet and asteroid surface has not been possible due to the paucity of spectra of comet nuclei at mid-infrared wavelengths. We present 5-35 μm thermal emission spectra of comets 10P/Tempel 2, and 49P/Arend-Rigaux observed with the Infrared Spectrograph on the Spitzer Space Telescope. Our analysis suggests the spectra are dominated by the comet nucleus. We fit each spectrum with the near-Earth asteroid thermal model (NEATM) and find sizes in agreement with previous values. However, the NEATM beaming parameters of the nuclei, 0.74 to 0.83, are systematically lower than the Jupiter-family comet population mean of 1.03+/-0.11, derived from 16- and 22-μm photometry. When the spectra are normalized by the NEATM model, a weak 10-μm silicate plateau is evident, with a shape similar to those seen in mid-infrared spectra of D-type asteroids. We compare, in detail, these comet nucleus emission features to those seen in spectra of the Jovian Trojan D-types (624) Hektor, (911) Agamemnon, and (1172) Aneas, as well as those seen in the spectra of seven comet comae. The comet comae present silicate features with two distinct shapes, either trapezoidal, or more rounded. The surfaces of Tempel 2, Arend-Rigaux, and Hektor best agree with the comae that present trapezoidal features. An emissivity minimum at 15 μm, present in the spectra of Tempel 2, Arend-Rigaux, Hektor, and Agamemnon, is also described, the origin of which remains unidentified. The compositional similarity between D-type asteroids and comets is discussed, and our data supports the hypothesis that they have similar origins in the early Solar System.

preprint2015arXiv

SOFIA Infrared Spectrophotometry of Comet C/2012 K1 (Pan-STARRS)

We present pre-perihelion infrared 8 to 31 micron spectrophotometric and imaging observations of comet C/2012 K1 (Pan-STARRS), a dynamically new Oort Cloud comet, conducted with NASA's Stratospheric Observatory for Infrared Astronomy (SOFIA) facility (+FORCAST) in 2014 June. As a "new" comet (first inner solar system passage), the coma grain population may be extremely pristine, unencumbered by a rime and insufficiently irradiated by the Sun to carbonize its surface organics. The comet exhibited a weak 10 micron silicate feature ~1.18 +/- 0.03 above the underlying best-fit 215.32 +/- 0.95 K continuum blackbody. Thermal modeling of the observed spectral energy distribution indicates that the coma grains are fractally solid with a porosity factor D = 3 and the peak in the grain size distribution, a_peak = 0.6 micron, large. The sub-micron coma grains are dominated by amorphous carbon, with a silicate-to-carbon ratio of 0.80 (+0.25) (- 0.20). The silicate crystalline mass fraction is 0.20 (+0.30) (-0.10), similar to with other dynamically new comets exhibiting weak 10 micron silicate features. The bolometric dust albedo of the coma dust is 0.14 +/- 0.01 at a phase angle of 34.76 degrees, and the average dust production rate, corrected to zero phase, at the epoch of our observations was Afrho ~ 5340~cm.

preprint2015arXiv

Variations between Dust and Gas in the Diffuse Interstellar Medium

Using the Planck far-infrared and Arecibo GALFA 21-cm line surveys, we identified a set of isolated interstellar clouds (approximately degree-sized on the sky and comprising 100 solar masses) and assessed the ratio of gas mass to dust mass. Significant variations of the gas-to-dust ratio are found both from cloud to cloud and within regions of individual clouds; within the clouds, the atomic gas per unit dust decreases by more than a factor of 3 compared to the standard gas-to-dust ratio. Three hypotheses are considered. First, the apparently low gas-to-dust ratio could be due to molecular gas. Comparing to Planck CO maps, the brightest clouds have a H2/CO ratio comparable to galactic plane clouds, but a strong lower limit is placed on the ratio for other clouds, such that the required amount of molecular gas is far higher than would be expected based on the CO upper limits. Second, we consider self-absorbed 21-cm lines and find the optical depth must be approximately 3, significantly higher than found from surveys of radio sources. Third, grain properties may change within the clouds: they become more emissive when they are colder, while not utilizing heavy elements that already have their cosmic abundance fully locked into grains. It is possible all three processes are active, and follow-up studies will be required to disentangle them and measure the true total gas and dust content of interstellar clouds.

preprint2014arXiv

Molecular hydrogen emission in the interstellar medium of the Large Magellanic Cloud

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

preprint2013arXiv

Survey of cometary CO2, CO, and particulate emissions using the Spitzer Space Telescope: Smog check for comets

We surveyed 23 comets using the Infrared Array Camera on the Spitzer Space Telescope in wide filters centered at 3.6 and 4.5 microns. Emission in the 3.6 micron filter arises from sunlight scattered by dust grains; these images generally have a coma near the nucleus and a tail in the antisolar direction due to dust grains swept back by solar radiation pressure. The 4.5 micron filter contains the same dust grains, as well as strong emission lines from CO2 and CO gas; these show distinct morphologies, in which cases we infer they are dominated by gas. Based on the ratio of 4.5 to 3.6 micron brightness, we classify the survey comets as CO2+CO "rich" and "poor." This classification is correlated with previous classifications by A'Hearn based on carbon-chain molecule abundance, in the sense that comets classified as "depleted" in carbon-chain molecules are also "poor" in CO2+CO. The gas emission in the IRAC 4.5 micron images is characterized by a smooth morphology, typically a fan in the sunward hemisphere with a radial profile that varies approximately as the inverse of projected distance from the nucleus, as would apply for constant production and free expansion. There are very significant radial and azimuthal enhancements in many of the comets, and these are often distinct between the gas and dust, indicating that ejection of solid material may be driven either by H2O or CO2. Notable features in the images include the following. There is a prominent loop of gas emission from 103P/Hartley 2, possible due to an outburst of CO2 before the Spitzer image. Prominent, double jets are present in the image of 88P/Howell. A prominent single jet is evident for 3 comets. Spirals are apparent in 29P and C/2006 W3; we measure a rotation rate of 21 hr for the latter comet. Arcs (possibly parts of a spiral) are apparent in the images of 10P/Tempel 2, and 2P/Encke.

preprint2013arXiv

The Discovery of Cometary Activity in Near-Earth Asteroid (3552) Don Quixote

The near-Earth object (NEO) population, which mainly consists of fragments from collisions between asteroids in the main asteroid belt, is thought to include contributions from short-period comets as well. One of the most promising NEO candidates for a cometary origin is near-Earth asteroid (3552) Don Quixote, which has never been reported to show activity. Here we present the discovery of cometary activity in Don Quixote based on thermal-infrared observations made with the Spitzer Space Telescope in its 3.6 and 4.5 μm bands. Our observations clearly show the presence of a coma and a tail in the 4.5 μm but not in the 3.6 μm band, which is consistent with molecular band emission from CO2. Thermal modeling of the combined photometric data on Don Quixote reveals a diameter of 18.4 (-0.4/+0.3) km and an albedo of 0.03 (-0.01/+0.02), which confirms Don Quixote to be the third-largest known NEO. We derive an upper limit on the dust production rate of 1.9 kg s^-1 and derive a CO2 gas production rate of (1.1+-0.1)10^26 molecules s^-1. Spitzer IRS spectroscopic observations indicate the presence of fine-grained silicates, perhaps pyroxene rich, on the surface of Don Quixote. Our discovery suggests that CO2 can be present in near-Earth space over a long time. The presence of CO2 might also explain that Don Quixote's cometary nature remained hidden for nearly three decades.

preprint2013arXiv

The Persistent Activity of Jupiter-Family Comets at 3 to 7 AU

We present an analysis of comet activity based on the Spitzer Space Telescope component of the Survey of the Ensemble Physical Properties of Cometary Nuclei. We show that the survey is well suited to measuring the activity of Jupiter-family comets at 3-7 AU from the Sun. Dust was detected in 33 of 89 targets (37 +/- 6%), and we conclude that 21 comets (24 +/- 5%) have morphologies that suggest ongoing or recent cometary activity. Our dust detections are sensitivity limited, therefore our measured activity rate is necessarily a lower limit. All comets with small perihelion distances (q < 1.8 AU) are inactive in our survey, and the active comets in our sample are strongly biased to post-perihelion epochs. We introduce the quantity epsilon-f-rho, intended to be a thermal emission counterpart to the often reported A-f-rho, and find that the comets with large perihelion distances likely have greater dust production rates than other comets in our survey at 3-7 AU from the Sun, indicating a bias in the discovered Jupiter-family comet population. By examining the orbital history of our survey sample, we suggest that comets perturbed to smaller perihelion distances in the past 150 yr are more likely to be active, but more study on this effect is needed.

preprint2012arXiv

A Spitzer IRAC Measure of the Zodiacal Light

The dominant non-instrumental background source for space-based infrared observatories is the zo- diacal light. We present Spitzer Infrared Array Camera (IRAC) measurements of the zodiacal light at 3.6, 4.5, 5.8, and 8.0 μm, taken as part of the instrument calibrations. We measure the changing surface brightness levels in approximately weekly IRAC observations near the north ecliptic pole (NEP) over the period of roughly 8.5 years. This long time baseline is crucial for measuring the annual sinusoidal variation in the signal levels due to the tilt of the dust disk with respect to the ecliptic, which is the true signal of the zodiacal light. This is compared to both Cosmic Background Explorer Diffuse Infrared Background Experiment (COBE DIRBE) data and a zodiacal light model based thereon. Our data show a few percent discrepancy from the Kelsall et al. (1998) model including a potential warping of the interplanetary dust disk and a previously detected overdensity in the dust cloud directly behind the Earth in its orbit. Accurate knowledge of the zodiacal light is important for both extragalactic and Galactic astronomy including measurements of the cosmic infrared background, absolute measures of extended sources, and comparison to extrasolar interplanetary dust models. IRAC data can be used to further inform and test future zodiacal light models.

preprint2012arXiv

Nucleosynthetic Layers in the Shocked Ejecta of Cassiopeia A

We present a 3-dimensional analysis of the supernova remnant Cassiopeia A using high resolution spectra from the Spitzer Space Telescope. We observe supernova ejecta both immediately before and during the shock-ejecta interaction. We determine that the reverse shock of the remnant is spherical to within 7%, although the center of this sphere is offset from the geometric center of the remnant by 810 km/s. We determine that the velocity width of the nucleosynthetic layers is approximately 1000 km/s over 4000 square arcsecond regions, although the velocity width of a layer along any individual line of sight is <250 km/s. Si and O, which come from different nucleosynthetic layers in the progenitor star, are observed to be coincident in velocity space in some directions, but segregated by up to approximately 500 km/s in other directions. We compare these observations of the nucleosynthetic layers to predictions from supernova explosion models in an attempt to constrain such models. Finally, we observe small-scale, corrugated velocity structures that are likely caused during the supernova explosion itself, rather than hundreds of years later by dynamical instabilities at the remnant's reverse shock.

preprint2012arXiv

Trajectories and Distribution of Interstellar Dust Grains in the Heliosphere

The solar wind carves a bubble in the surrounding interstellar medium (ISM), known as the heliosphere. Charged interstellar dust grains (ISDG) encountering the heliosphere may be diverted around the heliopause or penetrate it depending on their charge-to-mass ratio. We present new calculations of trajectories of ISDG in the heliosphere, and the dust density distributions that result. We include up-to-date grain charging calculations using a realistic UV radiation field and full 3-D magnetohydrodynamic fluid + kinetic models for the heliosphere. Models with two different (constant) polarities for the solar wind magnetic field (SWMF) are used, with the grain trajectory calculations done separately for each polarity. Small grains a_gr ~ 0.01 micron are completely excluded from the inner heliosphere. Large grains, a_gr ~ 1.0 micron pass into the inner solar system and are concentrated near the Sun by its gravity. Trajectories of intermediate size grains depend strongly on the SWMF polarity. When the field has magnetic north pointing to ecliptic north, the field de-focuses the grains resulting in low densities in the inner heliosphere, while for the opposite polarity the dust is focused near the Sun. The ISDG density outside the heliosphere inferred from applying the model results to in situ dust measurements is inconsistent with local ISM depletion data for both SWMF polarities, but is bracketed by them. This result points to the need to include the time variation in the SWMF polarity during grain propagation. Our results provide valuable insights for interpretation of the in situ dust observations from Ulysses.

preprint2011arXiv

Surveying the Agents of Galaxy Evolution in the Tidally-Stripped, Low Metallicity Small Magellanic Cloud (SAGE-SMC). I. Overview

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

preprint2010arXiv

Explosion of Comet 17P/Holmes as revealed by the Spitzer Space Telescope

An explosion on comet 17P/Holmes occurred on 2007 Oct 23, projecting particulate debris of a wide range of sizes into the interplanetary medium. We observed the comet using the Spitzer spectrograph on 2007 Nov 10 and 2008 Feb 27, and the photometer, on 2008 Mar 13. The fresh ejecta have detailed mineralogical features from small crystalline silicate grains. The 2008 Feb 27 spectra, and the central core of the 2007 Nov 10 spectral map, reveal nearly featureless spectra, due to much larger grains that were ejected from the nucleus more slowly. We break the infrared image into three components (size, speed) that also explain the temporal evolution of the mm-wave flux. Optical images were obtained on multiple dates spanning 2007 Oct 27 to 2008 Mar 10 at the Holloway Comet Observatory and 1.5-m telescope at Palomar Observatory. The orientation of the leading edge of the ejecta shell and the ejecta blob, relative to the nucleus, do not change as the orientation of the Sun changes; instead, the configuration was imprinted by the orientation of the initial explosion. The kinetic energy of the ejecta >1e21 erg is greater than the gravitational binding energy of the nucleus. We model the explosion as being due to crystallization and release of volatiles from interior amorphous ice within a subsurface cavity; once the pressure in the cavity exceeded the surface strength, the material above the cavity was propelled from the comet. The size of the cavity and the tensile strength of the upper layer of the nucleus are constrained by the observed properties of the ejecta; tensile strengths on >10 m scale must be greater than 10 kPa. The appearance of the 2007 outburst is similar to that witnessed in 1892, but the 1892 explosion was less energetic by a factor of about 20.

preprint2010arXiv

The Dust Trail of Comet 67P/Churyumov-Gerasimenko between 2004 and 2006

We report on observations of the dust trail of comet 67P/Churyumov-Gerasimenko (CG) in visible light with the Wide Field Imager at the ESO/MPG 2.2m telescope at 4.7 AU before aphelion, and at 24 micron with the MIPS instrument on board the Spitzer Space Telescope at 5.7 AU both before and after aphelion. The comet did not appear to be active during our observations. Our images probe large dust grains emitted from the comet that have a radiation pressure parameter beta<0.01. We compare our observations with simulated images generated with a dynamical model of the cometary dust and constrain the emission speeds, size distribution, production rate and geometric albedo of the dust. We achieve the best fit to our data with a differential size distribution exponent of -4.1, and emission speeds for a beta=0.01 particle of 25 m/s at perihelion and 2 m/s at 3 AU. The dust production rate in our model is on the order of 1000 kg/s at perihelion and 1 kg/s at 3 AU, and we require a dust geometric albedo between 0.022 and 0.044. The production rates of large (>10 micron) particles required to reproduce the brightness of the trail are sufficient to also account for the coma brightness observed while the comet was inside 3 AU, and we infer that the cross-section in the coma of CG may be dominated by grains of the order of 60-600 micron.

preprint2010arXiv

The Three-Dimensional Structure of Cassiopeia A

We used the Spitzer Space Telescope's Infrared Spectrograph to map nearly the entire extent of Cassiopeia A between 5-40 micron. Using infrared and Chandra X-ray Doppler velocity measurements, along with the locations of optical ejecta beyond the forward shock, we constructed a 3-D model of the remnant. The structure of Cas A can be characterized into a spherical component, a tilted thick disk, and multiple ejecta jets/pistons and optical fast-moving knots all populating the thick disk plane. The Bright Ring in Cas A identifies the intersection between the thick plane/pistons and a roughly spherical reverse shock. The ejecta pistons indicate a radial velocity gradient in the explosion. Some ejecta pistons are bipolar with oppositely-directed flows about the expansion center while some ejecta pistons show no such symmetry. Some ejecta pistons appear to maintain the integrity of the nuclear burning layers while others appear to have punched through the outer layers. The ejecta pistons indicate a radial velocity gradient in the explosion. In 3-D, the Fe jet in the southeast occupies a "hole" in the Si-group emission and does not represent "overturning", as previously thought. Although interaction with the circumstellar medium affects the detailed appearance of the remnant and may affect the visibility of the southeast Fe jet, the bulk of the symmetries and asymmetries in Cas A are intrinsic to the explosion.

preprint2010arXiv

The Three-Dimensional Structure of Interior Ejecta in Cassiopeia A at High Spectral Resolution

We used the Spitzer Space Telescope's Infrared Spectrograph to create a high resolution spectral map of the central region of the Cassiopeia A supernova remnant, allowing us to make a Doppler reconstruction of its 3D structure. The ejecta responsible for this emission have not yet encountered the remnant's reverse shock or the circumstellar medium, making it an ideal laboratory for exploring the dynamics of the supernova explosion itself. We observe that the O, Si, and S ejecta can form both sheet-like structures as well as filaments. Si and O, which come from different nucleosynthetic layers of the star, are observed to be coincident in velocity space in some regions, and separated by 500 km/s or more in others. Ejecta traveling toward us are, on average, ~900 km/s slower than the material traveling away from us. We compare our observations to recent supernova explosion models and find that no single model can simultaneously reproduce all the observed features. However, models of different supernova explosions can collectively produce the observed geometries and structures of the interior emission. We use the results from the models to address the conditions during the supernova explosion, concentrating on asymmetries in the shock structure. We also predict that the back surface of Cassiopeia A will begin brightening in ~30 years, and the front surface in ~100 years.

preprint2009arXiv

Distribution and properties of fragments and debris from the split comet 73P/Schwassmann-Wachmann 3 as revealed by Spitzer Space Telescope

During 2006 Mar - 2007 Jan, we used the IRAC and MIPS instruments on the Spitzer Space Telescope to study the infrared emission from the ensemble of fragments, meteoroids, and dust tails in the more than 3 degree wide 73P/Schwassmann-Wachmann 3 debris field. We also investigated contemporaneous ground based and HST observations. In 2006 May, 55 fragments were detected in the Spitzer image. The wide spread of fragments along the comet's orbit indicates they were formed from the 1995 splitting event. While the number of major fragments in the Spitzer image is similar to that seen from the ground by optical observers, the correspondence between the fragments with optical astrometry and those seen in the Spitzer images cannot be readily established, due either to strong non-gravitational terms, astrometric uncertainties, or transience of the fragments outgassing. The Spitzer data resolve the structure of the dust comae at a resolution of 1000 km, and they reveal the infrared emission due to large (mm to cm size) particles in a continuous dust trail that closely follows the projected orbit. We detect fluorescence from outflowing CO2 gas from the largest fragments (B and C), and we measure the CO2:H2O proportion (1:10 and 1:20, respectively). Three dimensionless parameters to explain dynamics of the solid particles: alpha (sublimation reaction), beta (radiation pressure), and nu (ejection velocity). The major fragments have nu>alpha>beta and are dominated by the kinetic energy imparted to them by the fragmentation process. The small, ephemeral fragments seen by HST in the tails of the major fragments have alpha>nu>beta dominated by rocket forces. The meteoroids along the projected orbit have beta~nu>>alpha. Dust in the fragments' tails has beta>>(nu+alpha) and is dominated by radiation pressure.

preprint2009arXiv

Exploring the Optical Transient Sky with the Palomar Transient Factory

The Palomar Transient Factory (PTF) is a wide-field experiment designed to investigate the optical transient and variable sky on time scales from minutes to years. PTF uses the CFH12k mosaic camera, with a field of view of 7.9 deg^2 and a plate scale of 1 asec/pixel, mounted on the the Palomar Observatory 48-inch Samuel Oschin Telescope. The PTF operation strategy is devised to probe the existing gaps in the transient phase space and to search for theoretically predicted, but not yet detected, phenomena, such as fallback supernovae, macronovae, .Ia supernovae and the orphan afterglows of gamma-ray bursts. PTF will also discover many new members of known source classes, from cataclysmic variables in their various avatars to supernovae and active galactic nuclei, and will provide important insights into understanding galactic dynamics (through RR Lyrae stars) and the Solar system (asteroids and near-Earth objects). The lessons that can be learned from PTF will be essential for the preparation of future large synoptic sky surveys like the Large Synoptic Survey Telescope. In this paper we present the scientific motivation for PTF and describe in detail the goals and expectations for this experiment.

preprint2009arXiv

Spatial variations of dust abundances across the Large Magellanic Cloud

Using the data obtained with the Spitzer Space telescope as part of the Surveying the Agents of a Galaxy's Evolution (SAGE) legacy survey, we have studied the variations of the dust composition and abundance across the Large Magellanic Cloud (LMC). Such variations are expected, as the explosive events which have lead to the formation of the many HI shells observed should have affected the dust properties. Using a model and comparing with a reference spectral energy distribution from our Galaxy, we deduce the relative abundance variations of small dust grains across the LMC. We examined the infrared color ratios as well as the relative abundances of very small grains (VSGs) and polycyclic aromatic hydrocarbons (PAHs) relative to the big grain (BG) abundance. Results show that each dust component could have different origins or evolution in the interstellar medium (ISM). The VSG abundance traces the star formation activity and could result from shattering of larger grains, whereas the PAH abundance increases around molecular clouds as well as in the stellar bar, where they could have been injected into the ISM during mass loss from old stars.

preprint2009arXiv

Spitzer Observations of Molecular Hydrogen in Interacting Supernova Remnants

With Spitzer IRS we have obtained sensitive low-resolution spectroscopy from 5 to 35 microns for six supernova remnants (SNRs) that show evidence of shocked molecular gas: Kes 69, 3C 396, Kes 17, G346.6-0.2, G348.5-0.0 and G349.7+0.2. Bright, pure-rotational lines of molecular hydrogen are detected at the shock front in all remnants, indicative of radiative cooling from shocks interacting with dense clouds. We find the excitation of H2 S(0)-S(7) lines in these SNRs requires two non-dissociative shock components: a slow, 10 km/s C- shock through clumps of density 10^6 cm^-3, and a faster, 40-70 km/s C- shock through a medium of density 10^4 cm^-3. The ortho-to-para ratio for molecular hydrogen in the warm shocked gas is typically found to be much less than the LTE value, suggesting that these SNRs are propagating into cold quiescent clouds. Additionally a total of thirteen atomic fine-structure transitions of Ar+, Ar++, Fe+, Ne+, Ne++, S++, and Si+ are detected. The ionic emitting regions are spatially segregated from the molecular emitting regions within the IRS slits. The presence of ionic lines with high appearance potential requires the presence of much faster, dissociative shocks through a lower density medium. The IRS slits are sufficiently wide to include regions outside the SNR which permits emission from diffuse gas around the remnants to be separated from the shocked emission. We find the diffuse molecular hydrogen gas projected outside the SNR is excited to a temperature of 100 to 300 K with a warm gas fraction of 0.5 to 15 percent along the line of sight.

preprint2009arXiv

Spitzer Observations of the Oldest White Dwarfs in the Solar Neighborhood

We present Spitzer 5-15 micron spectroscopy of one cool white dwarf and 3.6-8 micron photometry of 51 cool white dwarfs with T_eff < 6000 K. The majority of our targets have accurate BVRIJHK photometry and trigonometric parallax measurements available, which enables us to perform a detailed model atmosphere analysis using their optical, near- and mid-infrared photometry with state- of-the-art model atmospheres. We demonstrate that the optical and infrared spectral energy distributions of cool white dwarfs are well reproduced by our grid of models. Our best fit models are consistent with the observations within 5% in all filters except the IRAC 8 micron band, which has the lowest signal- to-noise ratio photometry. Excluding the ultracool white dwarfs, none of the stars in our sample show significant mid-infrared flux deficits or excesses. The non-detection of mid-infrared excess flux around our 2-9 Gyr old targets constrain the fraction of cool white dwarfs with warm debris disks to 0.8% (+1.5% -0.8%).

preprint2008arXiv

Spitzer Planet Limits around the Pulsating White Dwarf GD66

We present infrared observations in search of a planet around the white dwarf, GD66. Time-series photometry of GD66 shows a variation in the arrival time of stellar pulsations consistent with the presence of a planet with mass > 2.4Mj. Any such planet is too close to the star to be resolved, but the planet's light can be directly detected as an excess flux at 4.5um. We observed GD66 with the two shorter wavelength channels of IRAC on Spitzer but did not find strong evidence of a companion, placing an upper limit of 5--7Mj on the mass of the companion, assuming an age of 1.2--1.7Gyr.

preprint2008arXiv

The Dust cloud around the White Dwarf G 29-38. 2. Spectrum from 5-40 microns and mid-infrared variability

We model the mineralogy and distribution of dust around the white dwarf G29-39 using the infrared spectrum from 1-35 microns and combining a wide range of materials based on spectral studies of comets and debris disks. In order of their contribution to the mid-infrared emission, the most abundant minerals for G29-38 are amorphous carbon, amorphous and crystalline silicates, water ice, and metal sulfides. The amorphous C can be equivalently replaced by other materials (like metallic Fe) with featureless infrared spectra. The best-fitting crystalline silicate is Fe-rich pyroxene. In order to absorb enough starlight to power the observed emission, the disk must either be much thinner than the stellar radius (so that it can be heated from above and below) or it must have an opening angle wider than 2 degrees. A `moderately optically thick' torus model with mass 2x10^19 g fits the spectrum well. A physically thin (less than the white dwarf radius) and optically thick disk can contribute to the near-infrared continuum only; such a disk cannot explain the longer-wavelength continuum or strong silicate features. The silicate composition contains minerals found from cometary spectra and meteorites, but Fe-rich pyroxene is more abundant than enstatite (Mg-rich pyroxene) or forsterite (Mg-rich olivine) in G29-38, in contrast to what is found in most comet or meteorite mineralogies. Enstatite meteorites may be the most similar solar system materials to the G29-38 dust. Finally, we suggest the surviving core of a `hot jupiter' as an alternative (neither cometary nor asteroidal) origin for the debris, though further theoretical work is needed to determine if this hypothesis is viable.

preprint2007arXiv

The New Class of Dusty DAZ White Dwarfs

Our mid-infrared survey of 124 white dwarfs with the Spitzer Space Telescope and the IRAC imager has revealed an infrared excess associated with the white dwarf WD 2115-560 naturally explained by circumstellar dust. This object is the fourth white dwarf observed to have circumstellar dust. All four are DAZ white dwarfs, i.e. they have both photospheric Balmer lines and photospheric metal lines. We discuss these four objects as a class, which we abbreviate "DAZd", where the "d" stands for "dust". Using an optically-thick, geometrically-thin disk model analogous to Saturn's rings, we find that the inner disk edges are at >~0.1 to 0.2 Ro and that the outer disk edges are ~0.3 to 0.6 Ro. This model naturally explains the accretion rates and lifetimes of the detected WD disks and the accretion rates inferred from photospheric metal abundances.

preprint2004arXiv

Shocked Molecular Gas in the Supernova Remnants W 28 and W 44: Near-infrared and millimeter-wave observations

High resolution millimeter-wave and near-infrared observations of the supernova remnants W28 and W44 reveal extensive shocked molecular gas where supernova blast waves are propagating into giant molecular clouds. New CO observations were carried out with the IRAM 30-m and ARO 12-m telescopes, and the near-infrared observations were with Palomar 200-inch telescope. The near-infrared observations reveal shocked H2 emission from both supernova remnants, showing intricate networks of filaments on arcsec scales, following the bright ridges of the radio shells. The CO and CS linewidths, indicative of the shock speed, are 20-30 km/s. Both the near-infrared and millimeter-wave emission are attributed to shocks into gas with density >1e3 cm-3. Individual shock structures are resolved in the H2 emission, with inferred edge-on shock thickness ~1e17 cm, consistent with non-dissociative shocks into gas densities of 1e3-1e4 cm-3. Bright 1720 MHz OH masers are located within the shocked H2 gas complexes and highlight only localized areas where the conditions for masing are optimal. The Halpha and X-ray emission, have morphologies very different from the radio. We find a detailed correlation of the radio and H2 emission for some long filaments, indicating cosmic ray acceleration or re-acceleration due to the shocks into moderately dense gas. The different morphologies of these two remnants at different wavelengths is explained by a highly nonuniform structure for giant molecular clouds.

preprint2000arXiv

The formation of Encke meteoroids and dust trail

We observed comet 2P/Encke with the Infrared Space Observatory ISOCAM on July 14, 1997 from a particularly favorable viewing geometry above the comet's orbital plane and at a distance of 0.25 AU. A structured coma was observed, along with a long, straight dust trail. For the first time, we are able to observe the path of particles as they evolve from the nucleus to the trail. The particles that produce the infrared coma are large, with a radiation to gravitational force ratio beta<0.001 (corresponding to >mm-sized particles). The dust trail follows the orbit of the comet across our image, with a central core that is 20,000 km wide, composed of particles with beta<1e-5 (size $\sim 5$ cm) from previous apparitions. The abundant large particles near the comet pose a significant hazard to spacecraft. There is no evidence of a classical cometary dust tail due to small particles with beta>0.001, in marked contrast to other comets like P/Halley or C/Hale-Bopp. The structure of the coma requires anisotropic emission and that the spin axis of the nucleus to be nearly parallel to the orbital plane, resulting in strong seasonal variations of the particle emission. While most of the infrared coma emission is due to dust produced during the 1997 apparition, the core of the dust trail requires emissions from previous apparitions. The total mass lost during the 1997 apparition is estimated to be 2-6e13 g. Comparing to the gas mass loss from ultraviolet observations, the dust-to-gas mass ratio is 10-30, much higher than has ever been suggested from visual light observations. Using the recently-measured nuclear diameter, we find that Encke can only last 3000-10,000 rhoN yr (where rhoN is the nuclear density in g/cc) at its present mass loss rate.