Source author record

Robert A. Benjamin

Robert A. Benjamin 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

23works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

23 published item(s)

preprint2022arXiv

Evidence of a vertical kinematic oscillation beyond the Radcliffe Wave

The Radcliffe Wave (RW) is a recently discovered sinusoidal vertical feature of dense gas in the proximity of the Sun. In the disk plane, it is aligned with the Local Arm. However, the origin of its vertical undulation is still unknown. This study constrains the kinematics of the RW, using young stars and open clusters as tracers, and explores the possibility of this oscillation being part of a more extended vertical mode. We study the median vertical velocity trends of the young stars and clusters along with the RW and extend it further to the region beyond it. We discover a kinematic wave in the Galaxy, distinct from the warp, with the amplitude of oscillation depending on the age of the stellar population. We perform a similar analysis in the N-body simulation of a satellite as massive as the Sagittarius dwarf galaxy impacting the galactic disk. When projected in the plane, the spiral density wave induced by the satellite impact is aligned with the RW, suggesting that both may be the response of the disk to an external perturbation. However, the observed kinematic wave is misaligned. It appears as a kinematic wave travelling radially, winding up faster than the density wave matched by the RW, questioning its origin. If a satellite galaxy is responsible for this kinematic wave, we predict the existence of a vertical velocity dipole that should form across the disk and this may be measurable with the upcoming Gaia DR3 and DR4.

preprint2020arXiv

Discovery of diffuse optical emission lines from the inner Galaxy: Evidence for LI(N)ER-like gas

Optical emission lines are used to categorize galaxies into three groups according to their dominant central radiation source: active galactic nuclei, star formation, or low-ionization (nuclear) emission regions [LI(N)ERs] that may trace ionizing radiation from older stellar populations. Using the Wisconsin H-Alpha Mapper, we detect optical line emission in low-extinction windows within eight degrees of Galactic Center. The emission is associated with the 1.5-kiloparsec-radius "Tilted Disk" of neutral gas. We modify a model of this disk and find that the hydrogen gas observed is at least $48\%$ ionized. The ratio [NII] $λ$6584 $\overset{\lower.5em\circ}{\mathrm{A}}$/H$α$ $λ$6563 $\overset{\lower.5em\circ}{\mathrm{A}}$ increases from 0.3 to 2.5 with Galactocentric radius; [OIII] $λ$5007 $\overset{\lower.5em\circ}{\mathrm{A}}$ and H$β$ $λ$4861 $\overset{\lower.5em\circ}{\mathrm{A}}$ are also sometimes detected. The line ratios for most Tilted Disk sightlines are characteristic of LI(N)ER galaxies.

preprint2020arXiv

Discovery of High-Velocity H$α$ Emission in the Direction of the Fermi Bubble

Wisconsin H-Alpha Mapper (WHAM) observations reveal high-velocity and [NII]$\lambda6584$ emission lines in the same direction and velocity as ultraviolet absorption-line features that have been previously associated with the biconical gamma-ray lobes known as the Fermi Bubbles. We measure an extinction-corrected intensity of $I_{\textrm{H}α}=0.84^{+0.10}_{-0.09}$ Rayleigh for emission with line center $v_\textrm{LSR}=-221\pm3~\textrm{km}~\textrm{s}^{-1}$, corresponding to an emission measure of $EM = 2.00^{+0.64}_{-0.63}~\textrm{cm}^{-6}~\textrm{pc}$. This emission arises at the same velocity as Hubble Space Telescope/Cosmic Origins Spectrograph observations of ultraviolet absorption features detected in the PDS 456 quasar sight line that passes through the northern Bubble near $l = 10^\circ.4, b = +11^\circ.2$. We estimate the total column density of ionized gas in this velocity component to be $N(H^{+}) = \left(3.28 \pm 0.33\right) \times 10^{18}~\textrm{cm}^{-2}$. The comparison of ionized gas emission and absorption yields an estimate for the characteristic density of $n_{e,c} = 1.8 \pm 0.6~\textrm{cm}^{-3}$ and a characteristic length of $L_{c} =0.56 \pm 0.21~\textrm{pc}$ assuming $30\%$ solar metallicity. For a temperature of $T_{e}=8500^{+2700}_{-2600}$ K---consistent with the measured line widths and [NII]/H$α$ line ratio---the gas has a thermal pressure of $p/k = 32,000^{+15,000}_{-14,000}~\textrm{cm}^{-3}~\textrm{K}$. Assuming the gas is $\sim 6.5$ kpc distant, the derived density and pressure appear to be anomalously high for gas $\sim 1.3$ kpc above the Galactic midplane. The large thermal pressure is comparable to both a hot halo or Fermi Bubble model, but suggest that the H$α$ arises in an overpressurized zone.

preprint2020arXiv

Magnetic Field Morphology in Interstellar Clouds with the Velocity Gradient Technique

Magnetic fields, while ubiquitous in many astrophysical environments, are challenging to measure observationally. Based on the properties of anisotropy of eddies in magnetized turbulence, the Velocity Gradient Technique is a method synergistic to dust polarimetry that is capable of tracing plane-of-the-sky magnetic field, measuring the magnetization of interstellar media and estimating the fraction of gravitational collapsing gas in molecular clouds using spectral line observations. In this paper, we apply this technique to five low-mass star-forming molecular clouds in the Gould Belt and compare the results to the magnetic-field orientation obtained from polarized dust emission. We find the estimates of magnetic field orientations and magnetization for both methods are statistically similar. We estimate the fraction of collapsing gas in the selected clouds. By means of the Velocity Gradient Technique, we also present the plane-of-the-sky magnetic field orientation and magnetization of the Smith cloud, for which dust polarimetry data are unavailable.

preprint2020arXiv

The Effect of Bars on the Ionized ISM: Optical Emission Lines from Milky Way Analogs

Gas interior to the bar of the Milky Way has recently been shown as the closest example of a Low Ionization (Nuclear) Emission Region--LI(N)ER--in the universe. To better understand the nature of this gas, a sample of face-on galaxies with integral field spectroscopy are used to study the ionized gas conditions of 240 barred and 250 nonbarred galaxies, focusing on those that are most similar to the Milky Way. Strong optical line emission of $[NII]$ $λ6584$, H$α$, $[OIII]$ $λ5007$, and H$β$ are used to diagnose the dominant ionization mechanisms of gas across galaxies and the Galaxy via Baldwin-Phillips-Terlevich (BPT) Diagrams. Barred galaxies show a strong suppression of star formation and an increase in composite and LI(N)ER like spectra in their inner regions when compared with similar nonbarred counterparts. This effect is lessened in galaxies of very low ($\log_{10}(M_\star/M_\odot) \lesssim 10.4$) or very high ($\log_{10}(M_\star/M_\odot) \gtrsim 11.1$) total stellar mass. Bar masks from Galaxy Zoo:3D show the bar's non-axisymmetric effect on the ionized gas and help predict the face-on distribution of ionized gas conditions near the bar of the Milky Way.

preprint2020arXiv

VVV WIN 1733$-$3349: a low extinction window to probe the far side of the Milky Way bulge

Windows of low extinction in the Milky Way (MW) have been used along the past decades for the study of the Galactic structure and the stellar population across the inner bulge and disk. Here we report the analysis of another low extinction near-IR window discovered by the VISTA Variables in the Vía Láctea Survey. VVV WIN 1733$-$3349 is about half a degree in size and is conveniently located right in the MW plane, at Galactic coordinates $(l, b) = (-5.2, -0.3)$. The mean extinction of VVV WIN 1733$-$3349 is $A_{Ks} = 0.61\pm0.08$ mag, which is much smaller than the extinction in the surrounding area. The excess in the star counts is consistent with the reduced extinction, and complemented by studying the distribution of red clump (RC) stars. Thanks to the strategic low-latitude location of VVV WIN 1733$-$3349, we are able to interpret their RC density fluctuations with the expected overdensities due to the presence of the spiral arms beyond the bulge. In addition, we find a clear excess in the number of microlensing events within the window, which corroborates our interpretation that VVV WIN 1733$-$3349 is revealing the far side of the MW bulge.

preprint2016arXiv

Science Impacts of the SPHEREx All-Sky Optical to Near-Infrared Spectral Survey: Report of a Community Workshop Examining Extragalactic, Galactic, Stellar and Planetary Science

SPHEREx is a proposed SMEX mission selected for Phase A. SPHEREx will carry out the first all-sky spectral survey and provide for every 6.2" pixel a spectra between 0.75 and 4.18 $μ$m [with R$\sim$41.4] and 4.18 and 5.00 $μ$m [with R$\sim$135]. The SPHEREx team has proposed three specific science investigations to be carried out with this unique data set: cosmic inflation, interstellar and circumstellar ices, and the extra-galactic background light. It is readily apparent, however, that many other questions in astrophysics and planetary sciences could be addressed with the SPHEREx data. The SPHEREx team convened a community workshop in February 2016, with the intent of enlisting the aid of a larger group of scientists in defining these questions. This paper summarizes the rich and varied menu of investigations that was laid out. It includes studies of the composition of main belt and Trojan/Greek asteroids; mapping the zodiacal light with unprecedented spatial and spectral resolution; identifying and studying very low-metallicity stars; improving stellar parameters in order to better characterize transiting exoplanets; studying aliphatic and aromatic carbon-bearing molecules in the interstellar medium; mapping star formation rates in nearby galaxies; determining the redshift of clusters of galaxies; identifying high redshift quasars over the full sky; and providing a NIR spectrum for most eROSITA X-ray sources. All of these investigations, and others not listed here, can be carried out with the nominal all-sky spectra to be produced by SPHEREx. In addition, the workshop defined enhanced data products and user tools which would facilitate some of these scientific studies. Finally, the workshop noted the high degrees of synergy between SPHEREx and a number of other current or forthcoming programs, including JWST, WFIRST, Euclid, GAIA, K2/Kepler, TESS, eROSITA and LSST.

preprint2014arXiv

Halpha and [S II] emission from warm ionized gas in the Scutum-Centaurus Arm

We present Wisconsin H-Alpha Mapper [S II] λ6716 and Hα spectroscopic maps of the warm ionized medium (WIM) in the Scutum-Centaurus Arm at Galactic longitudes 310° < l < 345°. Using extinction-corrected Hα intensities (IHαc), we measure an exponential scale height of electron density-squared in the arm of H_ne^2 = 0.30 kpc (assuming a distance of 3.5 kpc), intermediate between that observed in the inner Galaxy and in the Perseus Arm. The [S II]/Hα line ratio is enhanced at large |z| and in sightlines with faint IHαc. We find that the [S II]/Hα line ratio has a power law relationship with IHαc from a value of ~=1.0 at IHαc < 0.2 R (Rayleighs) to a value of ~=0.08 at IHαc >= 100 R. The line ratio is better correlated with Hα intensity than with height above the plane, indicating that the physical conditions within the WIM vary systematically with electron density. We argue that the variation of the line ratio with height is a consequence of the decrease of electron density with height. Our results reinforce the well-established picture in which the diffuse Hα emission is due primarily to emission from in situ photoionized gas, with scattered light only a minor contributor.

preprint2014arXiv

Probing the Fermi Bubbles in Ultraviolet Absorption: A Spectroscopic Signature of the Milky Way's Biconical Nuclear Outflow

Giant lobes of plasma extend 55 degrees above and below the Galactic Center, glowing in emission from gamma rays (the Fermi Bubbles) to microwaves (the WMAP haze) and polarized radio waves. We use ultraviolet absorption-line spectra from the Hubble Space Telescope to constrain the velocity of the outflowing gas within these regions, targeting the quasar PDS 456 (Galactic coordinates l,b=10.4, +11.2 degrees). This sightline passes through a clear biconical structure seen in hard X-ray and gamma-ray emission near the base of the northern Fermi Bubble. We report two high-velocity metal absorption components, at v_LSR=-235 and +250 km/s, which cannot be explained by co-rotating gas in the Galactic disk or halo. Their velocities are suggestive of an origin on the front and back side of an expanding biconical outflow emanating from the Galactic Center. We develop simple kinematic biconical outflow models that can explain these observed profiles with an outflow velocity of ~900 km/s and a full opening angle of ~110 degrees (matching the X-ray bicone). This indicates Galactic Center activity over the last ~2.5-4.0 Myr, in line with age estimates of the Fermi Bubbles. The observations illustrate the use of UV absorption-line spectroscopy to probe the properties of swept-up gas venting into the Fermi Bubbles.

preprint2014arXiv

The Bones of the Milky Way

The very long and thin infrared dark cloud "Nessie" is even longer than had been previously claimed, and an analysis of its Galactic location suggests that it lies directly in the Milky Way's mid-plane, tracing out a highly elongated bone-like feature within the prominent Scutum-Centaurus spiral arm. Re-analysis of mid-infrared imagery from the Spitzer Space Telescope shows that this IRDC is at least 2, and possibly as many as 5 times longer than had originally been claimed by Nessie's discoverers (Jackson et al. 2010); its aspect ratio is therefore at least 300:1, and possibly as large as 800:1. A careful accounting for both the Sun's offset from the Galactic plane ($\sim 25$ pc) and the Galactic center's offset from the $(l^{II},b^{II})=(0,0)$ position shows that the latitude of the true Galactic mid-plane at the 3.1 kpc distance to the Scutum-Centaurus Arm is not $b=0$, but instead closer to $b=-0.4$, which is the latitude of Nessie to within a few pc. An analysis of the radial velocities of low-density (CO) and high-density (${\rm NH}_3$) gas associated with the Nessie dust feature suggests that Nessie runs along the Scutum-Centaurus Arm in position-position-velocity space, which means it likely forms a dense `spine' of the arm in real space as well. The Scutum-Centaurus arm is the closest major spiral arm to the Sun toward the inner Galaxy, and, at the longitude of Nessie, it is almost perpendicular to our line of sight, making Nessie the easiest feature to see as a shadow elongated along the Galactic Plane from our location. Future high-resolution dust mapping and molecular line observations of the harder-to-find Galactic "bones" should allow us to exploit the Sun's position above the plane to gain a (very foreshortened) view "from above" of the Milky Way's structure.

preprint2013arXiv

A Low Metallicity Molecular Cloud In The Lower Galactic Halo

We find evidence for the impact of infalling, low-metallicity gas on the Galactic disk. This is based on FUV absorption line spectra, 21-cm emission line spectra, and FIR mapping to estimate the abundance and physical properties of IV21 (IVC135+54-45), a galactic intermediate-velocity molecular cloud (IVMC) that lies ~300 pc above the disk. The metallicity of IV21 was estimated using observations toward the sdB star PG1144+615, located at a projected distance of 16 pc from the cloud's densest core, by measuring ion and HI column densities for comparison with known solar abundances. Despite the cloud's bright FIR emission and large column densities of molecular gas as traced by CO, we find that it has a sub-solar metallicity of log(Z/Z_Sun)=-0.43 +/- 0.12dex. IV21 is thus the first known sub-solar metallicity cloud in the solar neighborhood. In contrast, most intermediate-velocity clouds (IVC) have near-solar metallicities and are believed to originate in the Galactic Fountain. The cloud's low metallicity is also atypical for Galactic molecular clouds, especially in the light of the bright FIR emission which suggest a substantial dust content. The measured I_100mu/N(HI) ratio is a factor of three below the average found in high latitude \HI clouds within the solar neighborhood. We argue that IV21 represents the impact of an infalling, low-metallicity high-velocity cloud (HVC) that is mixing with disk gas in the lower Galactic halo.

preprint2013arXiv

Magnetized Gas in the Smith High Velocity Cloud

We report the first detection of magnetic fields associated with the Smith High Velocity Cloud. We use a catalog of Faraday rotation measures towards extragalactic radio sources behind the Smith Cloud, new HI observations from the Green Bank Telescope, and a spectroscopic map of Hα from the Wisconsin H-Alpha Mapper Northern Sky Survey. There are enhancements in rotation measure of approximately 100 rad m^(-2) which are generally well correlated with decelerated Hα emission. We estimate a lower limit on the line-of-sight component of the field of approximately 8 μG along a decelerated filament; this is a lower limit due to our assumptions about the geometry. No RM excess is evident in sightlines dominated by HI or Hα at the velocity of the Smith Cloud. The smooth Hα morphology of the emission at the Smith Cloud velocity suggests photoionization by the Galactic ionizing radiation field as the dominant ionization mechanism, while the filamentary morphology and high (approximately 1 Rayleigh) Hα intensity of the lower-velocity magnetized ionized gas suggests an ionization process associated with shocks due to interaction with the Galactic interstellar medium. The presence of the magnetic field may contribute to the survival of high velocity clouds like the Smith Cloud as they move from the Galactic halo to the disk. We expect these data to provide a test for magnetohydrodynamic simulations of infalling gas.

preprint2012arXiv

A self-consistent model of Galactic stellar and dust infrared emission and the abundance of polycyclic aromatic hydrocarbons

We present a self-consistent three-dimensional Monte-Carlo radiative transfer model of the stellar and dust emission in the Milky-Way, and have computed synthetic observations of the 3.6 to 100 microns emission in the Galactic mid-plane. In order to compare the model to observations, we use the GLIMPSE, MIPSGAL, and IRAS surveys to construct total emission spectra, as well as longitude and latitude profiles for the emission. The distribution of stars and dust is taken from the SKY model, and the dust emissivities includes an approximation of the emission from polycyclic aromatic hydrocarbons in addition to thermal emission. The model emission is in broad agreement with the observations, but a few modifications are needed to obtain a good fit. Firstly, by adjusting the model to include two major and two minor spiral arms rather than four equal spiral arms, the fit to the longitude profiles for |l|>30 degrees can be improved. Secondly, introducing a deficit in the dust distribution in the inner Galaxy results in a better fit to the shape of the IRAS longitude profiles at 60 and 100 microns. With these modifications, the model fits the observed profiles well, although it systematically under-estimates the 5.8 and 8.0 microns fluxes. One way to resolve this discrepancy is to increase the abundance of PAH molecules by 50% compared to the original model, although we note that changes to the dust distribution or radiation field may provide alternative solutions. Finally, we use the model to quantify which stellar populations contribute the most to the heating of different dust types, and which stellar populations and dust types contribute the most to the emission at different wavelengths.

preprint2012arXiv

Characterizing Transition Temperature Gas in the Galactic Corona

We present a study of the properties of the transition temperature (T~10^5 K) gas in the Milky Way corona, based on measurements of OVI, NV, CIV, SiIV and FeIII absorption lines seen in the far ultraviolet spectra of 58 sightlines to extragalactic targets, obtained with Far-Ultraviolet Spectroscopic Explorer (FUSE) and Space Telescope Imaging Spectrograph. In many sightlines the Galactic absorption profiles show multiple components, which are analyzed separately. We find that the highly-ionized atoms are distributed irregularly in a layer with a scaleheight of about 3 kpc, which rotates along with the gas in the disk, without an obvious gradient in the rotation velocity away from the Galactic plane. Within this layer the gas has randomly oriented velocities with a dispersion of 40-60 km/s. On average the integrated column densities are log N(OVI)=14.3, log N(NV)=13.5, log N(CIV)=14.2, log N(SiIV)=13.6 and log N(FeIII)=14.2, with a dispersion of just 0.2 dex in each case. In sightlines around the Galactic Center and Galactic North Pole all column densities are enhanced by a factor ~2, while at intermediate latitudes in the southern sky there is a deficit in N(OVI) of about a factor 2, but no deficit for the other ions. We compare the column densities and ionic ratios to a series of theoretical predictions: collisional ionization equilibrium, shock ionization, conductive interfaces, turbulent mixing, thick disk supernovae, static non-equilibrium ionization (NIE) radiative cooling and an NIE radiative cooling model in which the gas flows through the cooling zone. None of these models can fully reproduce the data, but it is clear that non-equilibrium ionization radiative cooling is important in generating the transition temperature gas.

preprint2012arXiv

Distances to Dark Clouds: Comparing Extinction Distances to Maser Parallax Distances

We test two different methods of using near-infrared extinction to estimate distances to dark clouds in the first quadrant of the Galaxy using large near infrared (2MASS and UKIDSS) surveys. VLBI parallax measurements of masers around massive young stars provide the most direct and bias-free measurement of the distance to these dark clouds. We compare the extinction distance estimates to these maser parallax distances. We also compare these distances to kinematic distances, including recent re-calibrations of the Galactic rotation curve. The extinction distance methods agree with the maser parallax distances (within the errors) between 66% and 100% of the time (depending on method and input survey) and between 85% and 100% of the time outside of the crowded Galactic center. Although the sample size is small, extinction distance methods reproduce maser parallax distances better than kinematic distances; furthermore, extinction distance methods do not suffer from the kinematic distance ambiguity. This validation gives us confidence that these extinction methods may be extended to additional dark clouds where maser parallaxes are not available.

preprint2012arXiv

Multi-wavelength Observations of the Radio Magnetar PSR J1622-4950 and Discovery of its Possibly Associated Supernova Remnant

We present multi-wavelength observations of the radio magnetar PSR J1622-4950 and its environment. Observations of PSR J1622-4950 with Chandra (in 2007 and 2009) and XMM (in 2011) show that the X-ray flux of PSR J1622-4950 has decreased by a factor of ~50 over 3.7 years, decaying exponentially with a characteristic time of 360 +/- 11 days. This behavior identifies PSR J1622-4950 as a possible addition to the small class of transient magnetars. The X-ray decay likely indicates that PSR J1622-4950 is recovering from an X-ray outburst that occurred earlier in 2007, before the 2007 Chandra observations. Observations with the Australia Telescope Compact Array show strong radio variability, including a possible radio flaring event at least one and a half years after the 2007 X-ray outburst that may be a direct result of this X-ray event. Radio observations with the Molonglo Observatory Synthesis Telescope reveal that PSR J1622-4950 is 8' southeast of a diffuse radio arc, G333.9+0.0, which appears non-thermal in nature and which could possibly be a previously undiscovered supernova remnant. If G333.9+0.0 is a supernova remnant then the estimates of its size and age, combined with the close proximity and reasonable implied velocity of PSR J1622-4950, suggests that these two objects could be physically associated.

preprint2012arXiv

Vertical structure of a supernova-driven turbulent magnetized ISM

Stellar feedback drives the circulation of matter from the disk to the halo of galaxies. We perform three-dimensional magnetohydrodynamic simulations of a vertical column of the interstellar medium with initial conditions typical of the solar circle in which supernovae drive turbulence and determine the vertical stratification of the medium. The simulations were run using a stable, positivity-preserving scheme for ideal MHD implemented in the FLASH code. We find that the majority (\approx 90 %) of the mass is contained in thermally-stable temperature regimes of cold molecular and atomic gas at T < 200 K or warm atomic and ionized gas at 5000 K < T < 10^{4.2} K, with strong peaks in probability distribution functions of temperature in both the cold and warm regimes. The 200 - 10^{4.2} K gas fills 50-60 % of the volume near the plane, with hotter gas associated with supernova remnants (30-40 %) and cold clouds (< 10 %) embedded within. At |z| ~ 1-2 kpc, transition-temperature (10^5 K) gas accounts for most of the mass and volume, while hot gas dominates at |z| > 3 kpc. The magnetic field in our models has no significant impact on the scale heights of gas in each temperature regime; the magnetic tension force is approximately equal to and opposite the magnetic pressure, so the addition of the field does not significantly affect the vertical support of the gas. The addition of a magnetic field does reduce the fraction of gas in the cold (< 200 K) regime with a corresponding increase in the fraction of warm (~ 10^4 K) gas. However, our models lack rotational shear and thus have no large-scale dynamo, which reduces the role of the field in the models compared to reality. The supernovae drive oscillations in the vertical distribution of halo gas, with the period of the oscillations ranging from ~ 30 Myr in the T < 200 K gas to ~ 100 Myr in the 10^6 K gas, in line with predictions by Walters & Cox.

preprint2011arXiv

Infrared Spectroscopy of Halos of Edge-on Galaxies

We present a study of ionized gas, PAHs, and molecular hydrogen emission in the halos of three edge-on galaxies, NGC 891, NGC 5775 and NGC 3044, based on 10-20 micron Spitzer Space Telescope spectra. The [Ne III]/[Ne II] ratio, an excellent measure of radiation hardness, rises with z in the halo of NGC 891. It is also higher in the halo of NGC 5775 than in the disk. NGC 3044 presents a more confusing situation. To explain the [Ne III]/[Ne II] as well as optical line ratio behavior in NGC 891, we carry out a simple exploration of parameter space with CLOUDY, which indicates a large increase in radiation temperature with height. Illustrative examples of physical models using a Monte Carlo radiative transfer code show that the rising neon ratio may be explained by adding a vertically extended, hot stellar source to a thin disk of massive stars. However, several other sources of hard spectra may be relevant. PAH features have scale heights of 430--530 pc in NGC 891 and 720--1080 pc in NGC 5775, suggesting they can be transported by disk-halo flows. Within NGC 891 and NGC 5775, scale heights are similar for all PAHs. For NGC 891, the scale heights exceed that of 8 micron emission, indicating a transition from more ionized to more neutral PAHs with height. Most PAH equivalent widths are higher in the halos. Molecular hydrogen 17.03 micron emission with scale heights of 550-580 pc in NGC 891 and 850 pc in NGC 5775 suggests a molecular component in a surprisingly thick layer.

preprint2011arXiv

MHD simulations of a supernova-driven ISM and the warm ionized medium using a positivity preserving ideal MHD scheme

We present new 3D magnetohydrodynamic (MHD) simulations of a supernova-driven, stratified interstellar medium. These simulations were run using the Waagan (2009) positivity preserving scheme for ideal MHD implemented in the Flash code. The scheme is stable even for the Mach numbers approaching 100 found in this problem. We have previously shown that the density distribution arising from hydrodynamical versions of these simulations creates low-density pathways through which Lyman continuum photons can travel to heights |z| > 1 kpc. This naturally produces the warm ionized medium through photoionization due primarily to O stars near the plane. However, our earlier models reproduce the peak but not the width of the observed emission measure distribution. Here, we examine whether inclusion of magnetic fields and a greater vertical extent to the simulation domain produce a gas distribution that better matches the observations. We further study the change of magnetic energy over time in our models, showing that it appears to reach a steady state after a few hundred megayears, presumably supported by a turbulent dynamo driven by the supernova explosions.

preprint2010arXiv

Photoionization of High Altitude Gas in a Supernova-Driven Turbulent Interstellar Medium

We investigate models for the photoionization of the widespread diffuse ionized gas in galaxies. In particular we address the long standing question of the penetration of Lyman continuum photons from sources close to the galactic midplane to large heights in the galactic halo. We find that recent hydrodynamical simulations of a supernova-driven interstellar medium have low density paths and voids that allow for ionizing photons from midplane OB stars to reach and ionize gas many kiloparsecs above the midplane. We find ionizing fluxes throughout our simulation grids are larger than predicted by one dimensional slab models, thus allowing for photoionization by O stars of low altitude neutral clouds in the Galaxy that are also detected in Halpha. In previous studies of such clouds the photoionization scenario had been rejected and the Halpha had been attributed to enhanced cosmic ray ionization or scattered light from midplane H II regions. We do find that the emission measure distributions in our simulations are wider than those derived from Halpha observations in the Milky Way. In addition, the horizontally averaged height dependence of the gas density in the hydrodynamical models is lower than inferred in the Galaxy. These discrepancies are likely due to the absence of magnetic fields in the hydrodynamic simulations and we discuss how magnetohydrodynamic effects may reconcile models and observations. Nevertheless, we anticipate that the inclusion of magnetic fields in the dynamical simulations will not alter our primary finding that midplane OB stars are capable of producing high altitude diffuse ionized gas in a realistic three-dimensional interstellar medium.

preprint2007arXiv

Integral Field Unit Observations of NGC 4302: Kinematics of the Diffuse Ionized Gas Halo

We present moderate resolution spectroscopy of extraplanar diffuse ionized gas (EDIG) emission in the edge-on spiral galaxy NGC 4302. The spectra were obtained with the SparsePak integral field unit (IFU) at the WIYN Observatory. The spectra are used to construct position-velocity (PV) diagrams at several ranges of heights above the midplane. Azimuthal velocities are directly extracted from the PV diagrams using the envelope tracing method, and indicate an extremely steep dropoff in rotational velocity with increasing height, with magnitude ~30 km/s/kpc. We find evidence for a radial variation in the velocity gradient on the receding side. We have also performed artificial observations of galaxy models in an attempt to match the PV diagrams. The results of a statistical analysis also favor a gradient of ~30 km/s/kpc. We compare these results with an entirely ballistic model of disk-halo flow, and find a strong dichotomy between the observed kinematics and those predicted by the model. The disagreement is worse than we have found for other galaxies in previous studies. The conclusions of this paper are compared to results from two other galaxies, NGC 5775 and NGC 891. We find that the vertical gradient in rotation speed, per unit EDIG scale height, for all three galaxies is consistent with a constant magnitude (within the errors) of approximately 15-25 km/s/scaleheight, independent of radius. This relationship is also true within the galaxy NGC 4302. We also discuss how the gradient depends on the distribution and morphology of the EDIG and the star formation rates of the galaxies, and consequences for the origin of the gas.

preprint2007arXiv

The Milky Way's Kiloparsec Scale Wind: A Hybrid Cosmic-Ray and Thermally Driven Outflow

We apply a wind model, driven by combined cosmic-ray and thermal-gas pressure, to the Milky Way, and show that the observed Galactic diffuse soft X-ray emission can be better explained by a wind than by previous static gas models. We find that cosmic-ray pressure is essential to driving the observed wind. Having thus defined a "best-fit" model for a Galactic wind, we explore variations in the base parameters and show how the wind's properties vary with changes in gas pressure, cosmic-ray pressure and density. We demonstrate the importance of cosmic rays in launching winds, and the effect cosmic rays have on wind dynamics. In addition, this model adds support to the hypothesis of Breitschwerdt and collaborators that such a wind may help explain the relatively small gradient observed in gamma-ray emission as a function of galactocentric radius.

preprint2006arXiv

Integral Field Unit Observations of NGC 891: Kinematics of the Diffuse Ionized Gas Halo

We present high and moderate spectral resolution spectroscopy of diffuse ionized gas (DIG) emission in the halo of NGC 891. The data were obtained with the SparsePak integral field unit at the WIYN Observatory. The wavelength coverage includes the [NII]6548,6583, Halpha, and [SII]6716,6731 emission lines. Position-velocity (PV) diagrams, constructed using spectra extracted from four SparsePak pointings in the halo, are used to examine the kinematics of the DIG. Using two independent methods, a vertical gradient in azimuthal velocity is found to be present in the northeast quadrant of the halo, with magnitude approximately 15-18 km/s/kpc, in agreement with results from HI observations. The kinematics of the DIG suggest that this gradient begins at approximately 1 kpc above the midplane. In another part of the halo, the southeast quadrant, the kinematics are markedly different, and suggest rotation at about 175 km/s, much slower than the disk but with no vertical gradient. We utilize an entirely ballistic model of disk-halo flow in an attempt to reproduce the kinematics observed in the northeast quadrant. Analysis shows that the velocity gradient predicted by the ballistic model is far too shallow. Based on intensity cuts made parallel to the major axis in the ballistic model and an Halpha image of NGC 891 from the literature, we conclude that the DIG halo is much more centrally concentrated than the model, suggesting that hydrodynamics dominate over ballistic motion in shaping the density structure of the halo. Velocity dispersion measurements along the minor axis of NGC 891 seem to indicate a lack of radial motions in the halo, but the uncertainties do not allow us to set firm limits.