Source author record

J. D. Smith

J. D. Smith 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

59works
8topics
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

59 published item(s)

preprint2022arXiv

Search for Spatial Correlations of Neutrinos with Ultra-High-Energy Cosmic Rays

For several decades, the origin of ultra-high-energy cosmic rays (UHECRs) has been an unsolved question of high-energy astrophysics. One approach for solving this puzzle is to correlate UHECRs with high-energy neutrinos, since neutrinos are a direct probe of hadronic interactions of cosmic rays and are not deflected by magnetic fields. In this paper, we present three different approaches for correlating the arrival directions of neutrinos with the arrival directions of UHECRs. The neutrino data is provided by the IceCube Neutrino Observatory and ANTARES, while the UHECR data with energies above $\sim$50 EeV is provided by the Pierre Auger Observatory and the Telescope Array. All experiments provide increased statistics and improved reconstructions with respect to our previous results reported in 2015. The first analysis uses a high-statistics neutrino sample optimized for point-source searches to search for excesses of neutrinos clustering in the vicinity of UHECR directions. The second analysis searches for an excess of UHECRs in the direction of the highest-energy neutrinos. The third analysis searches for an excess of pairs of UHECRs and highest-energy neutrinos on different angular scales. None of the analyses has found a significant excess, and previously reported over-fluctuations are reduced in significance. Based on these results, we further constrain the neutrino flux spatially correlated with UHECRs.

preprint2022arXiv

The State of the Molecular Gas in Post-Starburst Galaxies

The molecular gas in galaxies traces both the fuel for star formation and the processes that can enhance or suppress star formation. Observations of the molecular gas state can thus point to when and why galaxies stop forming stars. In this study, we present ALMA observations of the molecular gas in galaxies evolving through the post-starburst phase. These galaxies have low current star formation rates, regardless of the SFR tracer used, with recent starbursts ending within the last 600 Myr. We present CO (3-2) observations for three post-starburst galaxies, and dense gas HCN/HCO+/HNC (1-0) observations for six (four new) post-starburst galaxies. The post-starbursts have low excitation traced by the CO spectral line energy distribution (SLED) up to CO (3-2), more similar to early-type than starburst galaxies. The low excitation indicates that lower density rather than high temperatures may suppress star formation during the post-starburst phase. One galaxy displays a blueshifted outflow traced by CO (3-2). MaNGA observations show that the ionized gas velocity is disturbed relative to the stellar velocity field, with a blueshifted component aligned with the molecular gas outflow, suggestive of a multiphase outflow. Low ratios of HCO+/CO, indicating low fractions of dense molecular gas relative to the total molecular gas, are seen throughout post-starburst phase, except for the youngest post-starburst galaxy considered here. These observations indicate that the impact of any feedback or quenching processes may be limited to low excitation and weak outflows in the cold molecular gas during the post-starburst phase.

preprint2021arXiv

Direct Far-Infrared Metal Abundances (FIRA) I: M101

Accurately determining gas-phase metal-abundances within galaxies is critical as metals strongly affect the physics of the interstellar medium (ISM). To date, the vast majority of widely-used gas-phase abundance-indicators rely on emission from bright optical-lines, whose emissivities are highly sensitive to the electron temperature. Alternatively, direct-abundance methods exist that measure the temperature of the emitting gas directly, though these methods usually require challenging observations of highly-excited auroral lines. Low-lying far-infrared (FIR) fine-structure lines are largely insensitive to electron temperature and thus provide an attractive alternative to optically-derived abundances. Here, we introduce the far-infrared abundances (FIRA) project, which employs these FIR transitions, together with both radio free-free emission and hydrogen recombination-lines, to derive direct, absolute gas-phase oxygen-abundances. Our first target is M101, a nearby spiral-galaxy with a relatively steep abundance gradient. Our results are consistent with the O$^{++}$ electron-temperatures and absolute oxygen-abundances derived using optical direct-abundance methods by the CHemical Abundance Of Spirals (CHAOS) program, with a small difference ($\sim$ 1.5$σ$) in the radial abundance-gradients derived by the FIR/free-free-normalized vs. CHAOS/direct-abundance techniques. This initial result demonstrates the validity of the FIRA methodology $-$ with the promise of determining absolute metal-abundances within dusty star-forming galaxies, both locally and at high redshift.

preprint2021arXiv

Molecular Gas Properties and CO-to-H2 Conversion Factors in the Central Kiloparsec of NGC 3351

The CO-to-H$_2$ conversion factor ($α_\rm{CO}$) is critical to studying molecular gas and star formation in galaxies. The value of $α_\rm{CO}$ has been found to vary within and between galaxies, but the specific environmental conditions that cause these variations are not fully understood. Previous observations on $\sim$kpc scales revealed low values of $α_\rm{CO}$ in the centers of some barred spiral galaxies, including NGC 3351. We present new ALMA Band 3, 6, and 7 observations of $^{12}$CO, $^{13}$CO, and C$^{18}$O lines on 100 pc scales in the inner $\sim$2 kpc of NGC 3351. Using multi-line radiative transfer modeling and a Bayesian likelihood analysis, we infer the H$_2$ density, kinetic temperature, CO column density per line width, and CO isotopologue abundances on a pixel-by-pixel basis. Our modeling implies the existence of a dominant gas component with a density of $2{-}3\times10^3$ $\rm{cm^{-3}}$ in the central ${\sim}$1 kpc and a high temperature of 30$-$60 K near the nucleus and near the contact points that connect to the bar-driven inflows. Assuming a CO/H$_2$ abundance of $3\times10^{-4}$, our analysis yields $α_\rm{CO}{\sim}0.5{-}2.0$ $\rm{M_\odot\,(K~km~s^{-1}~pc^2)^{-1}}$ with a decreasing trend with galactocentric radius in the central $\sim$1 kpc. The inflows show a substantially lower $α_\rm{CO} < 0.1$ $\rm{M_\odot\,(K~km~s^{-1}~pc^2)^{-1}}$, likely due to lower optical depths caused by turbulence or shear in the inflows. Over the whole region, this gives an intensity-weighted $α_\rm{CO}$ of ${\sim}1.5$ $\rm{M_\odot\,(K~km~s^{-1}~pc^2)^{-1}}$, which is similar to previous dust modeling based results at kpc scales. This suggests that low $α_\rm{CO}$ on kpc scales in the centers of some barred galaxies may be due to the contribution of low optical depth CO emission in bar-driven inflows.

preprint2021arXiv

Observation of Variations in Cosmic Ray Single Count Rates During Thunderstorms and Implications for Large-Scale Electric Field Changes

We present the first observation by the Telescope Array Surface Detector (TASD) of the effect of thunderstorms on the development of cosmic ray single count rate intensity over a 700 km$^{2}$ area. Observations of variations in the secondary low-energy cosmic ray counting rate, using the TASD, allow us to study the electric field inside thunderstorms, on a large scale, as it progresses on top of the 700 km$^{2}$ detector, without dealing with the limitation of narrow exposure in time and space using balloons and aircraft detectors. In this work, variations in the cosmic ray intensity (single count rate) using the TASD, were studied and found to be on average at the $\sim(0.5-1)\%$ and up to 2\% level. These observations were found to be both in excess and in deficit. They were also found to be correlated with lightning in addition to thunderstorms. These variations lasted for tens of minutes; their footprint on the ground ranged from 6 to 24 km in diameter and moved in the same direction as the thunderstorm. With the use of simple electric field models inside the cloud and between cloud to ground, the observed variations in the cosmic ray single count rate were recreated using CORSIKA simulations. Depending on the electric field model used and the direction of the electric field in that model, the electric field magnitude that reproduces the observed low-energy cosmic ray single count rate variations was found to be approximately between 0.2-0.4 GV. This in turn allows us to get a reasonable insight on the electric field and its effect on cosmic ray air showers inside thunderstorms.

preprint2020arXiv

Evidence for a Supergalactic Structure of Magnetic Deflection Multiplets of Ultra-High Energy Cosmic Rays

Evidence for a large-scale supergalactic cosmic ray multiplet (arrival directions correlated with energy) structure is reported for ultra-high energy cosmic ray (UHECR) energies above 10$^{19}$ eV using seven years of data from the Telescope Array (TA) surface detector and updated to 10 years. Previous energy-position correlation studies have made assumptions regarding magnetic field shapes and strength, and UHECR composition. Here the assumption tested is that, since the supergalactic plane is a fit to the average matter density of the local Large Scale Structure (LSS), UHECR sources and intervening extragalactic magnetic fields are correlated with this plane. This supergalactic deflection hypothesis is tested by the entire field-of-view (FOV) behavior of the strength of intermediate-scale energy-angle correlations. These multiplets are measured in spherical cap section bins (wedges) of the FOV to account for coherent and random magnetic fields. The structure found is consistent with supergalactic deflection, the previously published energy spectrum anisotropy results of TA (the hotspot and coldspot), and toy-model simulations of a supergalactic magnetic sheet. The seven year data post-trial significance of this supergalactic structure of multiplets appearing by chance, on an isotropic sky, is found by Monte Carlo simulation to be 4.2$σ$. The ten years of data post-trial significance is 4.1$σ$. Furthermore, the starburst galaxy M82 is shown to be a possible source of the TA Hotspot, and an estimate of the supergalactic magnetic field using UHECR measurements is presented.

preprint2020arXiv

Measurement of the Proton-Air Cross Section with Telescope Array's Black Rock Mesa and Long Ridge Fluorescence Detectors, and Surface Array in Hybrid Mode

Ultra high energy cosmic rays provide the highest known energy source in the universe to measure proton cross sections. Though conditions for collecting such data are less controlled than an accelerator environment, current generation cosmic ray observatories have large enough exposures to collect significant statistics for a reliable measurement for energies above what can be attained in the lab. Cosmic ray measurements of cross section use atmospheric calorimetry to measure depth of air shower maximum ($X_{\mathrm{max}}$), which is related to the primary particle's energy and mass. The tail of the $X_{\mathrm{max}}$ distribution is assumed to be dominated by showers generated by protons, allowing measurement of the inelastic proton-air cross section. In this work the proton-air inelastic cross section measurement, $σ^{\mathrm{inel}}_{\mathrm{p-air}}$, using data observed by Telescope Array's Black Rock Mesa and Long Ridge fluorescence detectors and surface detector array in hybrid mode is presented. $σ^{\mathrm{inel}}_{\mathrm{p-air}}$ is observed to be $520.1 \pm 35.8$[Stat.] $^{+25.0}_{-40}$[Sys.]~mb at $\sqrt{s} = 73$ TeV. The total proton-proton cross section is subsequently inferred from Glauber formalism and is found to be $σ^{\mathrm{tot}}_{\mathrm{pp}} = 139.4 ^{+23.4}_{-21.3}$ [Stat.]$ ^{+15.0}_{-24.0}$[Sys.]~mb.

preprint2020arXiv

Observations of the Origin of Downward Terrestrial Gamma-Ray Flashes

In this paper we report the first close, high-resolution observations of downward-directed terrestrial gamma-ray flashes (TGFs) detected by the large-area Telescope Array cosmic ray observatory, obtained in conjunction with broadband VHF interferometer and fast electric field change measurements of the parent discharge. The results show that the TGFs occur during strong initial breakdown pulses (IBPs) in the first few milliseconds of negative cloud-to-ground and low-altitude intracloud flashes, and that the IBPs are produced by a newly-identified streamer-based discharge process called fast negative breakdown. The observations indicate the relativistic runaway electron avalanches (RREAs) responsible for producing the TGFs are initiated by embedded spark-like transient conducting events (TCEs) within the fast streamer system, and potentially also by individual fast streamers themselves. The TCEs are inferred to be the cause of impulsive sub-pulses that are characteristic features of classic IBP sferics. Additional development of the avalanches would be facilitated by the enhanced electric field ahead of the advancing front of the fast negative breakdown. In addition to showing the nature of IBPs and their enigmatic sub-pulses, the observations also provide a possible explanation for the unsolved question of how the streamer to leader transition occurs during the initial negative breakdown, namely as a result of strong currents flowing in the final stage of successive IBPs, extending backward through both the IBP itself and the negative streamer breakdown preceding the IBP.

preprint2020arXiv

Search for Large-scale Anisotropy on Arrival Directions of Ultra-high-energy Cosmic Rays Observed with the Telescope Array Experiment

Motivated by the detection of a significant dipole structure in the arrival directions of ultrahigh-energy cosmic rays above 8 EeV reported by the Pierre Auger Observatory (Auger), we search for a large-scale anisotropy using data collected with the surface detector array of the Telescope Array Experiment (TA). With 11 years of TA data, a dipole structure in a projection of the right ascension is fitted with an amplitude of 3.3+- 1.9% and a phase of 131 +- 33 degrees. The corresponding 99% confidence-level upper limit on the amplitude is 7.3%. At the current level of statistics, the fitted result is compatible with both an isotropic distribution and the dipole structure reported by Auger.

preprint2020arXiv

Search for point sources of ultra-high energy photons with the Telescope Array surface detector

The surface detector (SD) of the Telescope Array (TA) experiment allows one to indirectly detect photons with energies of order $10^{18}$ eV and higher and to separate photons from the cosmic-ray background. In this paper we present the results of a blind search for point sources of ultra-high energy (UHE) photons in the Northern sky using the TA SD data. The photon-induced extensive air showers (EAS) are separated from the hadron-induced EAS background by means of a multivariate classifier based upon 16 parameters that characterize the air shower events. No significant evidence for the photon point sources is found. The upper limits are set on the flux of photons from each particular direction in the sky within the TA field of view, according to the experiment's angular resolution for photons. Average 95% C.L. upper limits for the point-source flux of photons with energies greater than $10^{18}$, $10^{18.5}$, $10^{19}$, $10^{19.5}$ and $10^{20}$ eV are $0.094$, $0.029$, $0.010$, $0.0073$ and $0.0058$ km$^{-2}$yr$^{-1}$ respectively. For the energies higher than $10^{18.5}$ eV, the photon point-source limits are set for the first time. Numerical results for each given direction in each energy range are provided as a supplement to this paper.

preprint2016arXiv

First Upper Limits on the Radar Cross Section of Cosmic-Ray Induced Extensive Air Showers

TARA (Telescope Array Radar) is a cosmic ray radar detection experiment colocated with Telescope Array, the conventional surface scintillation detector (SD) and fluorescence telescope detector (FD) near Delta, Utah, U.S.A. The TARA detector combines a 40 kW, 54.1 MHz VHF transmitter and high-gain transmitting antenna which broadcasts the radar carrier over the SD array and within the FD field of view, towards a 250 MS/s DAQ receiver. TARA has been collecting data since 2013 with the primary goal of observing the radar signatures of extensive air showers (EAS). Simulations indicate that echoes are expected to be short in duration (~10 microseconds) and exhibit rapidly changing frequency, with rates on the order of 1 MHz/microsecond. The EAS radar cross-section (RCS) is currently unknown although it is the subject of over 70 years of speculation. A novel signal search technique is described in which the expected radar echo of a particular air shower is used as a matched filter template and compared to waveforms obtained by triggering the radar DAQ using the Telescope Array fluorescence detector. No evidence for the scattering of radio frequency radiation by EAS is obtained to date. We report the first quantitative RCS upper limits using EAS that triggered the Telescope Array Fluorescence Detector.

preprint2016arXiv

Ion clustering in aqueous salt solutions near the liquid/vapor interface

Molecular dynamics simulations of aqueous NaCl, KCl, NaI, and KI solutions are used to study the effects of salts on the properties of the liquid/vapor interface. The simulations use the models which include both charge transfer and polarization effects. Pairing and the formation of larger ion clusters occurs both in the bulk and surface region, with a decreased tendency to form larger clusters near the interface. An analysis of the roughness of the surface reveals that the chloride salts, which have less tendency to be near the surface, have a roughness that is less than pure water, while the iodide salts, which have a greater surface affinity, have a larger roughness. This suggests that ions away from the surface and ions near the surface affect the interface in opposite ways.

preprint2016arXiv

Search for correlations between the arrival directions of IceCube neutrino events and ultrahigh-energy cosmic rays detected by the Pierre Auger Observatory and the Telescope Array

This paper presents the results of different searches for correlations between very high-energy neutrino candidates detected by IceCube and the highest-energy cosmic rays measured by the Pierre Auger Observatory and the Telescope Array. We first consider samples of cascade neutrino events and of high-energy neutrino-induced muon tracks, which provided evidence for a neutrino flux of astrophysical origin, and study their cross-correlation with the ultrahigh-energy cosmic ray (UHECR) samples as a function of angular separation. We also study their possible directional correlations using a likelihood method stacking the neutrino arrival directions and adopting different assumptions on the size of the UHECR magnetic deflections. Finally, we perform another likelihood analysis stacking the UHECR directions and using a sample of through-going muon tracks optimized for neutrino point-source searches with sub-degree angular resolution. No indications of correlations at discovery level are obtained for any of the searches performed. The smallest of the p-values comes from the search for correlation between UHECRs with IceCube high-energy cascades, a result that should continue to be monitored.

preprint2015arXiv

A Cryogenic Space Telescope for Far-Infrared Astrophysics: A Vision for NASA in the 2020 Decade

Many of the transformative processes in the Universe have taken place in regions obscured by dust, and are best studied with far-IR spectroscopy. We present the Cryogenic-Aperture Large Infrared-Submillimeter Telescope Observatory (CALISTO), a 5-meter class, space-borne telescope actively cooled to 4 K, emphasizing moderate-resolution spectroscopy in the crucial 35 to 600 micron band. CALISTO will enable NASA and the world to study the rise of heavy elements in the Universe's first billion years, chart star formation and black hole growth in dust-obscured galaxies through cosmic time, and conduct a census of forming planetary systems in our region of the Galaxy. CALISTO will capitalize on rapid progress in both format and sensitivity of far-IR detectors. Arrays with a total count of a few 100,000 detector pixels will form the heart of a suite of imaging spectrometers in which each detector reaches the photon background limit. This document contains a large overview paper on CALISTO, as well as six 2-3 page scientific white papers, all prepared in response to NASA's Cosmic Origins Program Analysis Group (COPAG's) request for input on future mission concepts. The Far-IR Science Interest Group will meet from 3-5 June 2015 with the intention of reaching consensus on the architecture for the Far-IR Surveyor mission. This white paper describes one of the architectures to be considered by the community. One or more companion papers will describe alternative architectures.

preprint2015arXiv

A Northern Sky Survey for Point-Like Sources of EeV Neutral Particles with the Telescope Array Experiment

We report on the search for steady point-like sources of neutral particles around 10$^{18}$ eV between 2008 May and 2013 May with the scintillator surface detector of the Telescope Array experiment. We found overall no significant point-like excess above 0.5 EeV in the northern sky. Subsequently, we also searched for coincidence with the Fermi bright Galactic sources. No significant coincidence was found within the statistical uncertainty. Hence, we set an upper limit on the neutron flux that corresponds to an averaged flux of 0.07 km$^{-2}$ yr$^{-1}$ for $E>1$ EeV in the northern sky at the 95% confidence level. This is the most stringent flux upper limit in a northern sky survey assuming point-like sources. The upper limit at the 95% confidence level on the neutron flux from Cygnus X-3 is also set to 0.2 km$^{-2}$ yr$^{-1}$ for $E>0.5$ EeV. This is an order of magnitude lower than previous flux measurements.

preprint2015arXiv

Heating and cooling of the neutral ISM in the NGC4736 circumnuclear ring

The manner in which gas accretes and orbits within circumnuclear rings has direct implications for the star formation process. In particular, gas may be compressed and shocked at the inflow points, resulting in bursts of star formation at these locations. Afterwards the gas and young stars move together through the ring. In addition, star formation may occur throughout the ring, if and when the gas reaches sufficient density to collapse under gravity. These two scenarios for star formation in rings are often referred to as the `pearls on a string' and `popcorn' paradigms. In this paper, we use new Herschel PACS observations, obtained as part of the KINGFISH Open Time Key Program, along with archival Spitzer and ground-based observations from the SINGS Legacy project, to investigate the heating and cooling of the interstellar medium in the nearby star-forming ring galaxy, NGC4736. By comparing spatially resolved estimates of the stellar FUV flux available for heating, with the gas and dust cooling derived from the FIR continuum and line emission, we show that while star formation is indeed dominant at the inflow points in NGC 4736, additional star formation is needed to balance the gas heating and cooling throughout the ring. This additional component most likely arises from the general increase in gas density in the ring over its lifetime. Our data provide strong evidence, therefore, for a combination of the two paradigms for star formation in the ring in NGC4736.

preprint2015arXiv

Measurement of the Proton-Air Cross Section with Telescope Array's Middle Drum Detector and Surface Array in Hybrid Mode

In this work we are reporting on the measurement of the proton-air inelastic cross section $σ^{\rm inel}_{\rm p-air}$ using the Telescope Array (TA) detector. Based on the measurement of the $σ^{\rm inel}_{\rm p-air}$ the proton-proton cross section $σ_{\rm p-p}$ value is also determined at $\sqrt{s} = 95_{-8}^{+5}$ TeV. Detecting cosmic ray events at ultra high energies with Telescope Array enables us to study this fundamental parameter that we are otherwise unable to access with particle accelerators. The data used in this report is the hybrid events observed by the Middle Drum fluorescence detector together with the surface array detector collected over five years. The value of the $σ^{\rm inel}_{\rm p-air}$ is found to be equal to $567.0 \pm 70.5 [{\rm Stat.}] ^{+29}_{-25} [{\rm Sys.}]$ mb. The total proton-proton cross section is subsequently inferred from Glauber Formalism and Block, Halzen and Stanev QCD inspired fit and is found to be equal to $170_{-44}^{+48} [{\rm Stat.}] _{-17}^{+19} [{\rm Sys.}] $mb.

preprint2015arXiv

Spatially Resolved Spitzer-IRS Spectral Maps of the Superwind in M82

We have mapped the superwind/halo region of the nearby starburst galaxy M82 in the mid-infrared with $Spitzer-IRS$. The spectral regions covered include the H$_2 S(1)-S(3)$, [NeII], [NeIII] emission lines and PAH features. We estimate the total warm H$_2$ mass and the kinetic energy of the outflowing warm molecular gas to be between $M_{warm}\sim5-17\times10^6$ M$_{\odot}$ and $E_{K}\sim6-20\times10^{53}$ erg. Using the ratios of the 6.2, 7.7 and 11.3 micron PAH features in the IRS spectra, we are able to estimate the average size and ionization state of the small grains in the superwind. There are large variations in the PAH flux ratios throughout the outflow. The 11.3/7.7 and the 6.2/7.7 PAH ratios both vary by more than a factor of five across the wind region. The Northern part of the wind has a significant population of PAH's with smaller 6.2/7.7 ratios than either the starburst disk or the Southern wind, indicating that on average, PAH emitters are larger and more ionized. The warm molecular gas to PAH flux ratios (H$_2/PAH$) are enhanced in the outflow by factors of 10-100 as compared to the starburst disk. This enhancement in the H$_2/PAH$ ratio does not seem to follow the ionization of the atomic gas (as measured with the [NeIII]/[NeII] line flux ratio) in the outflow. This suggests that much of the warm H$_2$ in the outflow is excited by shocks. The observed H$_2$ line intensities can be reproduced with low velocity shocks ($v < 40$ km s$^{-1}$) driven into moderately dense molecular gas ($10^2 <n_H < 10^4$ cm$^{-3}$) entrained in the outflow.

preprint2015arXiv

The IceCube Neutrino Observatory, the Pierre Auger Observatory and the Telescope Array: Joint Contribution to the 34th International Cosmic Ray Conference (ICRC 2015)

We have conducted three searches for correlations between ultra-high energy cosmic rays detected by the Telescope Array and the Pierre Auger Observatory, and high-energy neutrino candidate events from IceCube. Two cross-correlation analyses with UHECRs are done: one with 39 cascades from the IceCube `high-energy starting events' sample and the other with 16 high-energy `track events'. The angular separation between the arrival directions of neutrinos and UHECRs is scanned over. The same events are also used in a separate search using a maximum likelihood approach, after the neutrino arrival directions are stacked. To estimate the significance we assume UHECR magnetic deflections to be inversely proportional to their energy, with values $3^\circ$, $6^\circ$ and $9^\circ$ at 100 EeV to allow for the uncertainties on the magnetic field strength and UHECR charge. A similar analysis is performed on stacked UHECR arrival directions and the IceCube sample of through-going muon track events which were optimized for neutrino point-source searches.

preprint2014arXiv

[CII] 158 $μ$m Emission as a Star Formation Tracer

The [CII] 157.74 $μ$m transition is the dominant coolant of the neutral interstellar gas, and has great potential as a star formation rate (SFR) tracer. Using the Herschel KINGFISH sample of 46 nearby galaxies, we investigate the relation of [CII] surface brightness and luminosity with SFR. We conclude that [CII] can be used for measurements of SFR on both global and kiloparsec scales in normal star-forming galaxies in the absence of strong active galactic nuclei (AGN). The uncertainty of the $Σ_{\rm [CII]}-Σ_{\rm SFR}$ calibration is $\pm$0.21 dex. The main source of scatter in the correlation is associated with regions that exhibit warm IR colors, and we provide an adjustment based on IR color that reduces the scatter. We show that the color-adjusted $Σ_{\rm[CII]}-Σ_{\rm SFR}$ correlation is valid over almost 5 orders of magnitude in $Σ_{\rm SFR}$, holding for both normal star-forming galaxies and non-AGN luminous infrared galaxies. Using [CII] luminosity instead of surface brightness to estimate SFR suffers from worse systematics, frequently underpredicting SFR in luminous infrared galaxies even after IR color adjustment (although this depends on the SFR measure employed). We suspect that surface brightness relations are better behaved than the luminosity relations because the former are more closely related to the local far-UV field strength, most likely the main parameter controlling the efficiency of the conversion of far-UV radiation into gas heating. A simple model based on Starburst99 population-synthesis code to connect SFR to [CII] finds that heating efficiencies are $1\%-3\%$ in normal galaxies.

preprint2014arXiv

Indications of Intermediate-Scale Anisotropy of Cosmic Rays with Energy Greater Than 57 EeV in the Northern Sky Measured with the Surface Detector of the Telescope Array Experiment

We have searched for intermediate-scale anisotropy in the arrival directions of ultrahigh-energy cosmic rays with energies above 57~EeV in the northern sky using data collected over a 5 year period by the surface detector of the Telescope Array experiment. We report on a cluster of events that we call the hotspot, found by oversampling using 20$^\circ$-radius circles. The hotspot has a Li-Ma statistical significance of 5.1$σ$, and is centered at R.A.=146.7$^{\circ}$, Dec.=43.2$^{\circ}$. The position of the hotspot is about 19$^{\circ}$ off of the supergalactic plane. The probability of a cluster of events of 5.1$σ$ significance, appearing by chance in an isotropic cosmic-ray sky, is estimated to be 3.7$\times$10$^{-4}$ (3.4$σ$).

preprint2014arXiv

Searches for Large-Scale Anisotropy in the Arrival Directions of Cosmic Rays Detected above Energy of $10^{19}$ eV at the Pierre Auger Observatory and the Telescope Array

Spherical harmonic moments are well-suited for capturing anisotropy at any scale in the flux of cosmic rays. An unambiguous measurement of the full set of spherical harmonic coefficients requires full-sky coverage. This can be achieved by combining data from observatories located in both the northern and southern hemispheres. To this end, a joint analysis using data recorded at the Telescope Array and the Pierre Auger Observatory above $10^{19}$ eV is presented in this work. The resulting multipolar expansion of the flux of cosmic rays allows us to perform a series of anisotropy searches, and in particular to report on the angular power spectrum of cosmic rays above $10^{19}$ eV. No significant deviation from isotropic expectations is found throughout the analyses performed. Upper limits on the amplitudes of the dipole and quadrupole moments are derived as a function of the direction in the sky, varying between 7% and 13% for the dipole and between 7% and 10% for a symmetric quadrupole.

preprint2014arXiv

Study of Ultra-High Energy Cosmic Ray Composition Using Telescope Array's Middle Drum Detector and Surface Array in Hybrid Mode

Previous measurements of the composition of Ultra-High Energy Cosmic Rays(UHECRs) made by the High Resolution Fly's Eye(HiRes) and Pierre Auger Observatory(PAO) are seemingly contradictory, but utilize different detection methods, as HiRes was a stereo detector and PAO is a hybrid detector. The five year Telescope Array(TA) Middle Drum hybrid composition measurement is similar in some, but not all, respects in methodology to PAO, and good agreement is evident between data and a light, largely protonic, composition when comparing the measurements to predictions obtained with the QGSJetII-03 and QGSJet-01c models. These models are also in agreement with previous HiRes stereo measurements, confirming the equivalence of the stereo and hybrid methods. The data is incompatible with a pure iron composition, for all models examined, over the available range of energies. The elongation rate and mean values of Xmax are in good agreement with Pierre Auger Observatory data. This analysis is presented using two methods: data cuts using simple geometrical variables and a new pattern recognition technique.

preprint2013arXiv

Energy Spectrum of Ultra-High Energy Cosmic Rays Observed with the Telescope Array Using a Hybrid Technique

We measure the spectrum of cosmic rays with energies greater than $10^{18.2}$ eV with the Fluorescence Detectors (FDs) and the Surface Detectors (SDs) of the Telescope Array Experiment using the data taken in our first 2.3-year observation from May 27 2008 to September 7 2010. A hybrid air shower reconstruction technique is employed to improve accuracies in determination of arrival directions and primary energies of cosmic rays using both FD and SD data. The energy spectrum presented here is in agreement with our previously published spectra and the HiRes results.

preprint2013arXiv

New Radial Abundance Gradients for NGC 628 and NGC 2403

Motived by recent ISM studies, we present high quality MMT and Gemini spectroscopic observations of H II regions in the nearby spiral galaxies NGC 628 and NGC 2403 in order to measure their chemical abundance gradients. Using long-slit and multi-object mask optical spectroscopy, we obtained measurements of the temperature sensitive auroral lines [O III] λ4363 and/or [N II] λ5755 at a strength of 4σ or greater in 11 H II regions in NGC 628 and 7 regions in NGC 2403. These observations allow us, for the first time, to derive an oxygen abundance gradient in NGC 628 based solely on "direct" oxygen abundances of H II regions: 12 + log(O/H) = (8.43+/-0.03) + (-0.017+/-0.002) x Rg (dex/kpc), with a dispersion in log(O/H) of σ = 0.10 dex, from 14 regions with a radial coverage of ~2-19 kpc. This is a significantly shallower slope than found by previous "strong-line" abundance studies. In NGC 2403, we derive an oxygen abundance gradient of 12 + log(O/H) = (8.48+/-0.04) + (-0.032+/-0.007) x Rg (dex/kpc), with a dispersion in log(O/H) of σ = 0.07 dex, from 7 H II with a radial coverage of ~1-10 kpc. Additionally, we measure the N, S, Ne, and Ar abundances. We find the N/O ratio decreases with increasing radius for the inner disk, but reaches a plateau past R25 in NGC 628. NGC 2403 also has a negative N/O gradient with radius, but we do not sample the outer disk of the galaxy past R25 and so do not see evidence for a plateau. This bi-modal pattern measured for NGC 628 indicates dominant contributions from secondary nitrogen inside of the R25 transition and dominantly primary nitrogen farther out. As expected for α-process elements, S/O, Ne/O, and Ar/O are consistent with constant values over a range in oxygen abundance.

preprint2013arXiv

Search for correlations of the arrival directions of ultra-high energy cosmic ray with extragalactic objects as observed by the telescope array experiment

We search for correlations between positions of extragalactic objects and arrival directions of Ultra-High Energy Cosmic Rays (UHECRs) with primary energy $E \ge 40$ EeV as observed by the surface detector array of the Telescope Array (TA) experiment during the first 40 months of operation. We examined several public astronomical object catalogs, including the Veron-Cetty and Veron catalog of active galactic nuclei. We counted the number of TA events correlated with objects in each catalog as a function of three parameters: the maximum angular separation between a TA event and an object, the minimum energy of the events, and the maximum redshift of the objects. We determine combinations of these parameters which maximize the correlations, and calculate the chance probabilities of having the same levels of correlations from an isotropic distribution of UHECR arrival directions. No statistically significant correlations are found when penalties for scanning over the above parameters and for searching in several catalogs are taken into account.

preprint2013arXiv

Shock Excited Molecules in NGC 1266: ULIRG conditions at the center of a Bulge Dominated Galaxy

We investigate the far infrared spectrum of NGC 1266, a S0 galaxy that contains a massive reservoir of highly excited molecular gas. Using the SPIRE-FTS, we detect the $^{12}$CO ladder up to J=(13-12), [C I] and [N II] lines, and also strong water lines more characteristic of UltraLuminous IR Galaxies (ULIRGs). The 12CO line emission is modeled with a combination of a low-velocity C-shock and a PDR. Shocks are required to produce the H2O and most of the high-J 12CO emission. Despite having an infrared luminosity thirty times less than a typical ULIRG, the spectral characteristics and physical conditions of the ISM of NGC 1266 closely resemble those of ULIRGs, which often harbor strong shocks and large-scale outflows.

preprint2013arXiv

The Cosmic Ray Energy Spectrum Observed with the Surface Detector of the Telescope Array Experiment

The Telescope Array (TA) collaboration has measured the energy spectrum of ultra-high energy cosmic rays with primary energies above 1.6 x 10^(18) eV. This measurement is based upon four years of observation by the surface detector component of TA. The spectrum shows a dip at an energy of 4.6 x 10^(18) eV and a steepening at 5.4 x 10^(19) eV which is consistent with the expectation from the GZK cutoff. We present the results of a technique, new to the analysis of ultra-high energy cosmic ray surface detector data, that involves generating a complete simulation of ultra-high energy cosmic rays striking the TA surface detector. The procedure starts with shower simulations using the CORSIKA Monte Carlo program where we have solved the problems caused by use of the "thinning" approximation. This simulation method allows us to make an accurate calculation of the acceptance of the detector for the energies concerned.

preprint2013arXiv

The Energy Spectrum of Ultra-High-Energy Cosmic Rays Measured by the Telescope Array FADC Fluorescence Detectors in Monocular Mode

We present a measurement of the energy spectrum of ultra-high-energy cosmic rays performed by the Telescope Array experiment using monocular observations from its two new FADC-based fluorescence detectors. After a short description of the experiment, we describe the data analysis and event reconstruction procedures. Since the aperture of the experiment must be calculated by Monte Carlo simulation, we describe this calculation and the comparisons of simulated and real data used to verify the validity of the aperture calculation. Finally, we present the energy spectrum calculated from the merged monocular data sets of the two FADC-based detectors, and also the combination of this merged spectrum with an independent, previously published monocular spectrum measurement performed by Telescope Array's third fluorescence detector (Abu-Zayyad {\it et al.}, {Astropart. Phys.} 39 (2012), 109). This combined spectrum corroborates the recently published Telescope Array surface detector spectrum (Abu-Zayyad {\it et al.}, {Astrophys. Journ.} 768 (2013), L1) with independent systematic uncertainties.

preprint2013arXiv

Towards A Removal of Temperature Dependencies from Abundance Determinations

The metal content of a galaxy is a key property for distinguishing between viable galaxy evolutionary scenarios, and it strongly influences many of the physical processes in the interstellar medium. An absolute and robust determination of extragalactic metallicities is essential in constraining models of chemical enrichment and chemical evolution, however, current gas phase abundance determinations from optical fine-structure lines are uncertain to 0.8 dex as conversion of these optical line fluxes to abundances is strongly dependent on the electron temperature of the ionized gas. In contrast, the far-IR emission lines can be used to derive an O++ abundance that is relatively insensitive to temperature, while the ratio of the optical to far-IR lines provides a consistent temperature to be used in the derivation of an O$^+$ abundance. We present observations of the [O III] 88 \micron\ fine-structure line in NGC 628 that were obtained as part of the KINGFISH program. These data are combined with optical IFU data to derive oxygen abundances for seven HII regions. We find the abundance of these regions to all lie between the high and low values of strong line calibrations and in agreement with estimates that assume temperature fluctuations are present in the HII regions.

preprint2013arXiv

Upper limit on the flux of photons with energies above 10^19 eV using the Telescope Array surface detector

We search for ultra-high energy photons by analyzing geometrical properties of shower fronts of events registered by the Telescope Array surface detector. By making use of an event-by-event statistical method, we derive upper limits on the absolute flux of primary photons with energies above 10^19, 10^19.5 and 10^20 eV based on the first three years of data taken.

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

Resolving the Far-IR Line Deficit: Photoelectric Heating and Far-IR Line Cooling in NGC 1097 and NGC 4559

The physical state of interstellar gas and dust is dependent on the processes which heat and cool this medium. To probe heating and cooling of the ISM over a large range of infrared surface brightness, on sub-kiloparsec scales, we employ line maps of [C \ii] 158 $μ$m, [O \one] 63 $μ$m, and [N \ii] 122 $μ$m in NGC 1097 and NGC 4559, obtained with the PACS spectrometer onboard {\it Herschel}. We matched new observations to existing Spitzer-IRS data that trace the total emission of polycyclic aromatic hydrocarbons (PAHs). We confirm at small scales in these galaxies that the canonical measure of photoelectric heating efficiency, ([C \ii] + [O \one])/TIR, decreases as the far-infrared color, $νf_ν$(70 $μ$m)/$νf_ν$(100 $μ$m), increases. In contrast, the ratio of far-infrared (far-IR) cooling to total PAH emission, ([C \ii] + [O \one])/PAH, is a near constant $\sim$6% over a wide range of far-infrared color, 0.5 \textless\ $νf_ν$(70 $μ$m)/$νf_ν$(100 $μ$m) $\lesssim$ 0.95. In the warmest regions, where $νf_ν$(70 $μ$m)/$νf_ν$(100 $μ$m) $\gtrsim$ 0.95, the ratio ([C \ii] + [O \one])/PAH drops rapidly to 4%. We derived representative values of the local UV radiation density, $G_0$, and the gas density, $n_H$, by comparing our observations to models of photodissociation regions. The ratio $G_0/n_H$, derived from fine-structure lines, is found to correlate with the mean dust-weighted starlight intensity, $<U>$ derived from models of the IR SED. Emission from regions that exhibit a line deficit is characterized by an intense radiation field, indicating that small grains are susceptible to ionization effects. We note that there is a shift in the 7.7 / 11.3 $μ$m PAH ratio in regions that exhibit a deficit in ([C \ii] + [O \one])/PAH, suggesting that small grains are ionized in these environments.

preprint2012arXiv

Search for Anisotropy of Ultra-High Energy Cosmic Rays with the Telescope Array Experiment

We study the anisotropy of Ultra-High Energy Cosmic Ray (UHECR) events collected by the Telescope Array (TA) detector in the first 40 months of operation. Following earlier studies, we examine event sets with energy thresholds of 10 EeV, 40 EeV, and 57 EeV. We find that the distributions of the events in right ascension and declination are compatible with an isotropic distribution in all three sets. We then compare with previously reported clustering of the UHECR events at small angular scales. No significant clustering is found in the TA data. We then check the events with E>57 EeV for correlations with nearby active galactic nuclei. No significant correlation is found. Finally, we examine all three sets for correlations with the large-scale structure of the Universe. We find that the two higher-energy sets are compatible with both an isotropic distribution and the hypothesis that UHECR sources follow the matter distribution of the Universe (the LSS hypothesis), while the event set with E>10 EeV is compatible with isotropy and is not compatible with the LSS hypothesis at 95% CL unless large deflection angles are also assumed. We show that accounting for UHECR deflections in a realistic model of the Galactic magnetic field can make this set compatible with the LSS hypothesis.

preprint2012arXiv

The Energy Spectrum of Telescope Array's Middle Drum Detector and the Direct Comparison to the High Resolution Fly's Eye Experiment

The Telescope Array's Middle Drum fluorescence detector was instrumented with telescopes refurbished from the High Resolution Fly's Eye's HiRes-1 site. The data observed by Middle Drum in monocular mode was analyzed via the HiRes-1 profile-constrained geometry reconstruction technique and utilized the same calibration techniques enabling a direct comparison of the energy spectra and energy scales between the two experiments. The spectrum measured using the Middle Drum telescopes is based on a three-year exposure collected between December 16, 2007 and December 16, 2010. The calculated difference between the spectrum of the Middle Drum observations and the published spectrum obtained by the data collected by the HiRes-1 site allows the HiRes-1 energy scale to be transferred to Middle Drum. The HiRes energy scale is applied to the entire Telescope Array by making a comparison between Middle Drum monocular events and hybrid events that triggered both Middle Drum and the Telescope Array's scintillator Ground Array.

preprint2012arXiv

The surface detector array of the Telescope Array experiment

The Telescope Array (TA) experiment, located in the western desert of Utah,USA, is designed for observation of extensive air showers from extremely high energy cosmic rays. The experiment has a surface detector array surrounded by three fluorescence detectors to enable simultaneous detection of shower particles at ground level and fluorescence photons along the shower track. The TA surface detectors and fluorescence detectors started full hybrid observation in March, 2008. In this article we describe the design and technical features of the TA surface detector.

preprint2011arXiv

New air fluorescence detectors employed in the Telescope Array experiment

Since 2007, the Telescope Array (TA) experiment, based in Utah, USA, has been observing ultra high energy cosmic rays to understand their origins. The experiment involves a surface detector (SD) array and three fluorescence detector (FD) stations. FD stations, installed surrounding the SD array, measure the air fluorescence light emitted from extensive air showers (EASs) for precise determination of their energies and species. The detectors employed at one of the three FD stations were relocated from the High Resolution Fly's Eye experiment. At the other two stations, newly designed detectors were constructed for the TA experiment. An FD consists of a primary mirror and a camera equipped with photomultiplier tubes. To obtain the EAS parameters with high accuracies, understanding the FD optical characteristics is important. In this paper, we report the characteristics and installation of new FDs and the performances of the FD components. The results of the monitored mirror reflectance during the observation time are also described in this report.

preprint2011arXiv

Spitzer IRS Spectral Mapping of the Toomre Sequence: Spatial Variations of PAH, Gas, and Dust Properties in Nearby Major Mergers

We have mapped the key mid-IR diagnostics in eight major merger systems of the Toomre Sequence (NGC4676, NGC7592, NGC6621, NGC2623, NGC6240, NGC520, NGC3921, and NGC7252) using the Spitzer Infrared Spectrograph (IRS). With these maps, we explore the variation of the ionized-gas, PAH, and warm-gas (H_2) properties across the sequence and within the galaxies. While the global PAH interband strength and ionized gas flux ratios ([Ne III]/[Ne II]) are similar to those of normal star forming galaxies, the distribution of the spatially resolved PAH and fine structure line flux ratios is significant different from one system to the other. Rather than a constant H_2/PAH flux ratio, we find that the relation between the H_2 and PAH fluxes is characterized by a power law with a roughly constant exponent (0.61+/-0.05) over all merger components and spatial scales. While following the same power law on local scales, three galaxies have a factor of ten larger integrated (i.e. global) H_2/PAH flux ratio than the rest of the sample, even larger than what it is in most nearby AGNs. These findings suggest a common dominant excitation mechanism for H_2 emission over a large range of global H_2/PAH flux ratios in major mergers. Early merger systems show a different distribution between the cold (CO J=1-0) and warm (H_2) molecular gas component, which is likely due to the merger interaction. Strong evidence for buried star formation in the overlap region of the merging galaxies is found in two merger systems (NGC6621 and NGC7592) as seen in the PAH, [Ne II], [Ne III], and warm gas line emission, but with no apparent corresponding CO (J=1-0) emission. Our findings also demonstrate that the variations of the physical conditions within a merger are much larger than any systematic trends along the Toomre Sequence.

preprint2011arXiv

The deeply obscured AGN of NGC4945 I. Spitzer-IRS maps of [Ne V], [Ne II], H2 0-0 S(1), S(2), and other tracers

The nearly edge-on galaxy NGC4945 is one of the closest galaxies where an AGN and starburst coexist, and is one of the brightest sources at 100 keV. Near and mid-infrared spectroscopy have shown very strong obscuration of its central region, rivaled only in strength by some of the most deeply obscured ULIRGs. We aim to determine the spatial distribution of ISM features in the central 426x426 pc^2 of NGC4945. We map the central region of NGC4945 in three Spitzer-IRS modules (SH, SL and LL). We produce maps of the flux distribution of the starburst tracers [Ne II], [Ne III], [S III] and [S IV] at 12.81, 15.56, 18.71 and 10.51 mum, respectively, and a map of the AGN narrow-line region tracer [Ne V] at 14.32 mum. We also mapped the S(1), S(2) and S(3) pure rotational lines of H2, which trace the distribution of warm molecular hydrogen. We obtained an extinction map (A_V) based on the apparent strength of the 9.7 mum silicate absorption feature. Our A_V map traces the contours of the starburst ring but the highest extinction (A_V(9.85 mum)~60) is found at the nucleus. Within the uncertainty of the astrometry all emission lines are found to peak on the nucleus, except for the warm molecular hydrogen emission which shows a maximum 60-100 pc NW of the nucleus. We favour a scenario in which the lower H2 0-0 S(1) and S(2) rotational lines originate mainly from an unobscured extra-nuclear component associated with the super-wind cone observed in the HST NICMOS map of the H2 1-0 S(1) vibrational line. For the [Ne V] emission we infer an attenuation of a factor 12-160 (A_V=55-112) based on a comparison of the ratio of our [Ne V] flux and the absorption-corrected 14-195 keV Swift-BAT flux to the average [Ne V]/BAT ratio for Seyfert 1 nuclei. The high attenuation indicates that [Ne V] and [O IV] cannot be used as extinction-free tracers of AGN power in galaxies with deeply buried nuclei.

preprint2011arXiv

The Displaced Dusty ISM of NGC 3077: Tidal Stripping in the M 81 Triplet

We present the detection of extended (~30 kpc^2) dust emission in the tidal \hi\ arm near NGC 3077 (member of the M\,81 triplet) using SPIRE on board Herschel. Dust emission in the tidal arm is typically detected where the \hi\ column densities are >10^21 cm^-2. The SPIRE band ratios show that the dust in the tidal arm is significantly colder (~13 K) than in NGC 3077 itself (~31 K), consistent with the lower radiation field in the tidal arm. The total dust mass in the tidal arm is ~1.8 x 10^6 M_sun (assuming beta=2), i.e. substantially larger than the dust mass associated with NGC 3077 (~2 x 10^5 M_sun). Where dust is detected, the dust-to-gas ratio is 6+/-3 x 10^-3, consistent within the uncertainties with what is found in NGC 3077 and nearby spiral galaxies with Galactic metallicities. The faint HII regions in the tidal arm can not be responsible for the detected enriched material and are not the main source of the dust heating in the tidal arm. We conclude that the interstellar medium (atomic HI, molecules and dust) in this tidal feature was pre-enriched and stripped off NGC 3077 during its recent interaction (~3 x 10^8 yr ago) with M 82 and M 81. This implies that interaction can efficiently remove heavy elements and enriched material (dust, molecular gas) from galaxies. As interactions were more frequent at large lookback times, it is conceivable that they could substantially contribute (along with galactic outflows) to the enrichment of the intergalactic medium.

preprint2011arXiv

The Emission by Dust and Stars of Nearby Galaxies in the Herschel KINGFISH Survey

Using new far-infrared imaging from the Herschel Space Observatory with ancillary data from ultraviolet to submillimeter wavelengths, we estimate the total emission from dust and stars of 62 nearby galaxies in the KINGFISH survey in a way that is as empirical and model-independent as possible. We collect and exploit these data in order to measure from the spectral energy distributions (SEDs) precisely how much stellar radiation is intercepted and re-radiated by dust, and how this quantity varies with galaxy properties. By including SPIRE data, we are more sensitive to emission from cold dust grains than previous analyses at shorter wavelengths, allowing for more accurate estimates of dust temperatures and masses. The dust/stellar flux ratio, which we measure by integrating the SEDs, has a range of nearly three decades. The inclusion of SPIRE data shows that estimates based on data not reaching these far-IR wavelengths are biased low. We find that the dust/stellar flux ratio varies with morphology and total IR luminosity. We also find that dust/stellar flux ratios are related to gas-phase metallicity, while the dust/stellar mass ratios are less so. The substantial scatter between dust/stellar flux and dust/stellar mass indicates that the former is a poor proxy of the latter. Comparing the dust/stellar flux ratios and dust temperatures, we show that early-types tend to have slightly warmer temperatures than spiral galaxies, which may be due to more intense interstellar radiation fields, or to different dust grain compositions. Finally, we show that early-types and early-type spirals have a strong correlation between the dust/stellar flux ratio and specific star formation rate, which suggests that the relatively bright far-IR emission of some of these galaxies is due to ongoing star formation and the radiation field from older stars.

preprint2010arXiv

Analysis of large-scale anisotropy of ultra-high energy cosmic rays in HiRes data

Stereo data collected by the HiRes experiment over a six year period are examined for large-scale anisotropy related to the inhomogeneous distribution of matter in the nearby Universe. We consider the generic case of small cosmic-ray deflections and a large number of sources tracing the matter distribution. In this matter tracer model the expected cosmic ray flux depends essentially on a single free parameter, the typical deflection angle theta. We find that the HiRes data with threshold energies of 40 EeV and 57 EeV are incompatible with the matter tracer model at a 95% confidence level unless theta is larger than 10 degrees and are compatible with an isotropic flux. The data set above 10 EeV is compatible with both the matter tracer model and an isotropic flux.

preprint2010arXiv

Indications of Proton-Dominated Cosmic Ray Composition above 1.6 EeV

We report studies of ultra-high energy cosmic ray composition via analysis of depth of airshower maximum (Xmax), for airshower events collected by the High Resolution Fly's Eye (HiRes) observatory. The HiRes data are consistent with a constant elongation rate d<Xmax>/d(log(E)) of 47.9 +- 6.0 (stat.) +- 3.2 (syst.)g/cm^2/decade for energies between 1.6 EeV and 63 EeV, and are consistent with a predominantly protonic composition of cosmic rays when interpreted via the QGSJET01 and QGSJET-II high-energy hadronic interaction models. These measurements constrain models in which the galactic-to-extragalactic transition is the cause of the energy spectrum "ankle' at 4 EeV.

preprint2010arXiv

Infrared Luminosities and Aromatic Features in the 24um Flux Limited Sample of 5MUSES

We study a 24\,$μ$m selected sample of 330 galaxies observed with the Infrared Spectrograph for the 5\,mJy Unbiased Spitzer Extragalactic Survey. We estimate accurate total infrared luminosities by combining mid-IR spectroscopy and mid-to-far infrared photometry, and by utilizing new empirical spectral templates from {\em Spitzer} data. The infrared luminosities of this sample range mostly from 10$^9$L$_\odot$ to $10^{13.5}$L$_\odot$, with 83% in the range 10$^{10}$L$_\odot$$<$L$_{\rm IR}$$<10^{12}$L$_\odot$. The redshifts range from 0.008 to 4.27, with a median of 0.144. The equivalent widths of the 6.2\,$μ$m aromatic feature have a bimodal distribution. We use the 6.2\,$μ$m PAH EW to classify our objects as SB-dominated (44%), SB-AGN composite (22%), and AGN-dominated (34%). The high EW objects (SB-dominated) tend to have steeper mid-IR to far-IR spectral slopes and lower L$_{\rm IR}$ and redshifts. The low EW objects (AGN-dominated) tend to have less steep spectral slopes and higher L$_{\rm IR}$ and redshifts. This dichotomy leads to a gross correlation between EW and slope, which does not hold within either group. AGN dominated sources tend to have lower log(L$_{\rm PAH 7.7μm}$/L$_{\rm PAH 11.3μm}$) ratios than star-forming galaxies, possibly due to preferential destruction of the smaller aromatics by the AGN. The log(L$_{\rm PAH 7.7μm}$/L$_{\rm PAH 11.3μm}$) ratios for star-forming galaxies are lower in our sample than the ratios measured from the nuclear spectra of nearby normal galaxies, most probably indicating a difference in the ionization state or grain size distribution between the nuclear regions and the entire galaxy. Finally, we provide a calibration relating the monochromatic 5.8, 8, 14 and 24um continuum or Aromatic Feature luminosity to L$_{\rm IR}$ for different types of objects.

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

Physical Conditions in the Ionized Gas of 30 Doradus

We present a mid-infrared spectroscopic data cube of the central part of 30 Doradus, observed with Spitzer's IRS and MIPS/SED mode. Aromatic dust emission features and emission lines from molecular and atomic hydrogen are detected but not particularly strong. The dominant spectral features are emission lines from moderately ionized species of argon, neon, and sulphur, which are used to determine the physical conditions in the ionized gas. The ionized gas excitation shows strong variations on parsec scales, some of which can plausibly be associated with individual hot stars. We fit the ionic line strengths with photoionization and shock models, and find that photoionization dominates in the region. The ionization parameter U traces the rim of the central bubble, as well as highlighting isolated sources of ionization, and at least one quiescent clump. The hardness of the ionizing radiation field T_rad reveals several "hot spots" that are either the result of individual very hot stars or trace the propagation of the diffuse ionizing field through the surrounding neutral cloud. Consistent with other measurements of giant molecular hydrogen regions, log(U) ranges between -3 and -0.75, and T_rad between 30000 and 85000K.

preprint2009arXiv

Powerful H$_2$ Line-cooling in Stephan's Quintet : I - Mapping the Significant Cooling Pathways in Group-wide Shocks

We present results from the mid-infrared spectral mapping of Stephan's Quintet using the Spitzer Space Telescope. A 1000 km/s collision has produced a group-wide shock and for the first time the large-scale distribution of warm molecular hydrogen emission is revealed, as well as its close association with known shock structures. In the main shock region alone we find 5.0 $\times10^{8}$ M$_{\odot}$ of warm H$_2$ spread over $\sim$ 480 kpc$^2$ and additionally report the discovery of a second major shock-excited H$_2$ feature. This brings the total H$_2$ line luminosity of the group in excess of 10$^42$ erg/s. In the main shock, the H$_2$ line luminosity exceeds, by a factor of three, the X-ray luminosity from the hot shocked gas, confirming that the H$_2$-cooling pathway dominates over the X-ray. [Si II]34.82$μ$m emission, detected at a luminosity of 1/10th of that of the H$_2$, appears to trace the group-wide shock closely and in addition, we detect weak [FeII]25.99$μ$m emission from the most X-ray luminous part of the shock. Comparison with shock models reveals that this emission is consistent with regions of fast shocks (100 < $V_{s}$ < 300 km/s) experiencing depletion of iron and silicon onto dust grains. Star formation in the shock (as traced via ionic lines, PAH and dust emission) appears in the intruder galaxy, but most strikingly at either end of the radio shock. The shock ridge itself shows little star formation, consistent with a model in which the tremendous H$_{2}$ power is driven by turbulent energy transfer from motions in a post-shocked layer. The significance of the molecular hydrogen lines over other measured sources of cooling in fast galaxy-scale shocks may have crucial implications for the cooling of gas in the assembly of the first galaxies.

preprint2009arXiv

Spitzer-IRS Study of the Antennae Galaxies NGC 4038/39

Using the Infrared Spectrograph on the Spitzer Space Telescope, we observed the Antennae galaxies obtaining spectral maps of the entire central region and high signal-to-noise 5-38 um spectra of the two galactic nuclei and six infrared-luminous regions. The total infrared luminosity of our six IR peaks plus the two nuclei is L_IR = 3.8x10^10 L_o, with their derived star formation rates ranging between 0.2 and 2 M_o/yr, with a total of 6.6 M_o/yr. The hardest and most luminous radiation originates from two compact clusters in the southern part of the overlap region, which also have the highest dust temperatures. PAH emission and other tracers of softer radiation are spatially extended throughout and beyond the overlap region, but regions with harder and intenser radiation field show a reduced PAH strength. The strong H_2 emission is rather confined around the nucleus of NGC 4039, where shocks appear to be the dominant excitation mechanism, and the southern part of the overlap region, where it traces the most recent starburst activity. The luminosity ratio between the warm molecular gas (traced by the H_2 lines) and the total far-IR emission is ~1.6x10^-4, similar to that found in many starburst and ULIRGs. The total mass of warm H_2 in the Antennae is 2.5x10^7 M_o, with a fraction of warm to total H_2 gas mass of about 0.35%. The average warm H_2 temperature is 302+/-26 K and appears anti-correlated with the radiation field hardness, possibly due to an evolution of the PDR morphology. The previously reported tight correlation between the H_2 and PAH emission was not found but higher total PAH emission to continuum ratios were found in PDRs with warmer gas.

preprint2008arXiv

New Young Brown Dwarfs in the Orion Molecular Cloud 2/3 Region

Forty new low mass members with spectral types ranging from M4-M9 have been confirmed in the Orion Molecular Cloud 2/3 region. Through deep, I, z', J, H, K photometry of a 20' x 20' field in OMC 2/3, we selected brown dwarf candidates for follow-up spectroscopy. Low resolution far-red and near-infrared spectra were obtained for the candidates, and 19 young brown dwarfs in the OMC 2/3 region are confirmed. They exhibit spectral types of M6.5-M9, corresponding to approximate masses of 0.075-0.015 M_solar using the evolutionary models of Baraffe et al. (1998). At least one of these bona fide young brown dwarfs has strong Halpha emission, indicating that it is actively accreting. In addition, we confirm 21 new low mass members with spectral types of M4-M6, corresponding to approximate masses of 0.35-0.10 M_solar in OMC 2/3. By comparing pre-main sequence tracks to the positions of the members in the H-R diagram, we find that most of the brown dwarfs are less than 1 Myr, but find a number of low mass stars with inferred ages greater than 3 Myr. The discrepancy in the stellar and substellar ages is due to our selection of only low luminosity sources; however, the presence of such objects implies the presence of an age spread in the OMC 2/3 region. We discuss possible reasons for this apparent age spread.

preprint2007arXiv

Analysis of the emission of very small dust particles from Spitzer spectro-imagery data using blind signal separation methods

This work was conducted as part of the SPECPDR program, dedicated to the study of very small particles and astrochemistry, in Photo-Dissociation Regions (PDRs). We present the analysis of the mid-IR spectro-imagery observations of Ced 201, NCG 7023 East and North-West and rho-Ophiuchi West filament. Using the data from all four modules of the InfraRed Spectrograph onboard the Spitzer Space Telescope, we produced a spectral cube ranging from 5 to 35 um for each one of the observed PDRs. The resulting cubes were analysed using Blind Signal Separation methods (NMF and FastICA). For Ced 201, rho-Ophiuchi West filament and NGC 7023 East, we find that two signals can be extracted from the original data cubes, which are 5 to 35 um spectra. The main features of the first spectrum are a strong continuum emission at long wavelengths, and a broad 7.8 um band. On the contrary, the second spectrum exhibits the classical Aromatic Infrared Bands (AIBs) and no continuum. The reconstructed spatial distribution maps show that the latter spectrum is mainly present at the cloud surface, close to the star whereas the first one is located slightly deeper inside the PDR. The study of the spectral energy distribution of Ced 201 up to 100 um suggests that, in cool PDRs, the 5-25 um continuum is carried by Very Small Grains (VSGs). The AIB spectra in the observed objects can be interpreted as the contribution of neutral and positively-charged Polycyclic Aromatic Hydrocarbons (PAHs). We extracted the 5 to 25 um emission spectrum of VSGs in cool PDRs, these grains being most likely carbonaceous. We show that the variations of the mid-IR (5-35 um) spectra of PDRs can be explained by the photo-chemical processing of these VSGs and PAHs, VSGs being the progenitors of free PAH.

preprint2007arXiv

Dust and Atomic Gas in Dwarf Irregular Galaxies of the M81 Group: The SINGS and THINGS view

[abridged] We present observations of the dust and atomic gas phase in seven dwarf irregular galaxies of the M81 group from the SINGS and THINGS surveys. The Spitzer observations provide a first glimpse of the nature of the non-atomic ISM in these metal-poor (Z~0.1 Z_sun), quiescent (SFR~0.001-0.1 M_sun/yr) dwarf galaxies. Dust emission is detected in five out of the seven targets. Most detected dust emission is restricted to HI column densities >1x10^21 cm^-2. Spitzer spectroscopy of two regions in the brightest galaxies (IC 2574 and Holmberg II) show distinctly different spectral shapes. The spectrum of IC 2574 shows aromatic features that are less luminous (relative to the FIR luminosity) compared to an average SINGS spiral galaxy by a factor of \~7 . The aromatic features in Holmberg~II (which has only a slightly lower gas-phase metallicity) are fainter than in IC 2574 by an order of magnitude. This result emphazises that the strength of the aromatic features is not a simple linear function of metallicity. We estimate dust masses of ~10^4-10^6 M_sun for the M81 dwarf galaxies, resulting in an average dust--to--gas ratio (M_dust/M_HI) of ~3x10^-4 (1.5x10^-3 if only the HI that is associated with dust emission is considered); this is an order of magnitude lower than the typical value derived for the SINGS spirals. The dwarf galaxies are underluminous per unit star formation rate at 70um as compared to the more massive galaxies in SINGS by a factor of ~2. However, the average 70um/160um ratio in the sample dwarf galaxies is higher than what is found in the other galaxies of the SINGS sample. This can be explained by a combination of a lower dust content in conjunction with a higher dust temperature in the dwarfs.

preprint2007arXiv

Spitzer IRS mapping of the central kpc of Centaurus A

We report on the results of infrared spectroscopic mapping observations carried out in the nuclear region of Centaurus A (NGC5128). The 500 pc bipolar dust shell discovered by Quillen et al.(2006) is even more clearly seen in the 11.3 micron dust emission feature than previous broad band imaging. The pure rotational lines of molecular hydrogen other than the S(0) line are detected above the dusty disk and associated with the oval dust shell. The molecular hydrogen transitions indicate the presence of warm gas at temperatures 250--720K. The ratio of the dust emission features at 7.7 and 11.3 micron and the ratio of the [NeII](12.8) and 11.3 dust emission feature are lower in the 500 pc dust shell than in the star forming disk. The clearer shell morphology at 11.3 micron, warm molecular hydrogen emission in the shell, and variation in line ratios in the shell compared to those in the disk, confirm spectroscopically that this shell is a separate coherent entity and is unlikely to be a chance superposition of dust filaments. The physical conditions in the shell are most similar to Galactic supernova remnants where blast waves encounter molecular clouds. The lines requiring the highest level of ionization, [NeV] and [OIV], are detected 20--25 arcsec north-east and south-west of the nucleus and at position angles near the radio jet axis. Fine structure line ratios and limits from this region suggest that the medium is low density and illuminated by a hard radiation field at low ionization parameter. These higher S molecular hydrogen pure rotational transitions are also particularly bright in the same region as the [OIV] and [NeV] emission. This suggests that the gas associated with the dust shell has been excited near the jet axis and is part of an ionization cone.

preprint2007arXiv

Warm molecular hydrogen in the Spitzer SINGS galaxy sample

(simplified) Results on the properties of warm H2 in 57 normal galaxies are derived from H2 rotational transitions, obtained as part of SINGS. This study extends previous extragalactic surveys of H2, the most abundant constituent of the molecular ISM, to more common systems (L_FIR = e7 to 6e10 L_sun) of all morphological and nuclear types. The S(1) transition is securely detected in the nuclear regions of 86% of SINGS galaxies with stellar masses above 10^9.5 M_sun. The derived column densities of warm H2 (T > ~100 K), even though averaged over kiloparsec-scale areas, are commensurate with those of resolved PDRs; the median of the sample is 3e20 cm-2. They amount to between 1% and >30% of the total H2. The power emitted in the sum of the S(0) to S(2) transitions is on average 30% of the [SiII] line power, and ~4e-4 of the total infrared power (TIR) within the same area for star-forming galaxies, which is consistent with excitation in PDRs. The fact that H2 emission scales tightly with PAH emission, even though the average radiation field intensity varies by a factor ten, can also be understood if both tracers originate predominantly in PDRs, either dense or diffuse. A large fraction of the 25 LINER/Sy targets, however, strongly depart from the rest of the sample, in having warmer H2 in the excited states, and an excess of H2 emission with respect to PAHs, TIR and [SiII]. We propose a threshold in H2 to PAH power ratios, allowing the identification of low-luminosity AGNs by an excess H2 excitation. A dominant contribution from shock heating is favored in these objects. Finally, we detect, in nearly half the star-forming targets, non-equilibrium ortho to para ratios, consistent with FUV pumping combined with incomplete ortho-para thermalization by collisions, or possibly non-equilibrium PDR fronts advancing into cold gas.

preprint2006arXiv

Discovery of a 500 pc shell in the nucleus of Centaurus A

Spitzer Space Telescope mid-infrared images of the radio galaxy Centaurus A reveal a shell-like, bipolar, structure 500 pc to the north and south of the nucleus. This shell is seen in 5.8, 8.0 and 24 micron broad-band images. Such a remarkable shell has not been previously detected in a radio galaxy and is the first extragalactic nuclear shell detected at mid-infrared wavelengths. We estimate that the shell is a few million years old and has a mass of order million solar masses. A conservative estimate for the mechanical energy in the wind driven bubble is 10^53 erg. The shell could have created by a small few thousand solar mass nuclear burst of star formation. Alternatively, the bolometric luminosity of the active nucleus is sufficiently large that it could power the shell. Constraints on the shell's velocity are lacking. However, if the shell is moving at 1000 km/s then the required mechanical energy would be 100 times larger.

preprint2000arXiv

Discovery of a Bright Field Methane (T-type) Brown Dwarf by 2MASS

We report the discovery of a bright (J = 13.83$\pm$0.03) methane brown dwarf, or T dwarf, by the Two Micron All Sky Survey. This object, 2MASSI J0559191-140448, is the first brown dwarf identified by the newly commissioned CorMASS instrument mounted on the Palomar 60-inch Telescope. Near-infrared spectra from 0.9 - 2.35 $\micron$ show characteristic CH$_4$ bands at 1.1, 1.3, 1.6, and 2.2 $\micron$, which are significantly shallower than those seen in other T dwarfs discovered to date. Coupled with the detection of an FeH band at 0.9896 $\micron$ and two sets of K I doublets at J-band, we propose that 2MASS J0559-14 is a warm T dwarf, close to the transition between L and T spectral classes. The brightness of this object makes it a good candidate for detailed investigation over a broad wavelength regime and at higher resolution.

preprint1999arXiv

Mid-Infrared Spectra of Be Stars

We present the first medium-resolution ($R\sim 600$) mid-infrared (8-13.3\micron) spectra of 11 Be stars. A large number of lines are observed and identified in these spectra, including, as an example, 39 hydrogen recombination lines in the spectrum of $γ$ Cas. In the majority of our spectra, all of the observed lines are attributable to hydrogen recombination. Two of the sources, $β$ Lyr and MWC 349 also show emission from other species. Both of these objects show evidence of [Ne II] emission, and $β$ Lyr also shows evidence of He I emission. We tabulate the effective line strength and line widths for the observed lines, and briefly discuss the physical implications of the observed line series. We also use a simple model of free-free emission to characterize the disks around these sources.