Exploring the Dust Content of Galactic Winds with Herschel. I. NGC 4631

We present a detailed analysis of deep far-infrared observations of the nearby edge-on star-forming galaxy NGC 4631 obtained with the Herschel Space Observatory. Our PACS images at 70 and 160 um show a rich complex of filaments and chimney-like features that extends up to a projected distance of 6 kpc above the plane of the galaxy. The PACS features often match extraplanar Halpha, radio-continuum, and soft X-ray features observed in this galaxy, pointing to a tight disk-halo connection regulated by star formation. On the other hand, the morphology of the colder dust component detected on larger scale in the SPIRE 250, 350, and 500 um data matches the extraplanar H~I streams previously reported in NGC 4631 and suggests a tidal origin. The PACS 70/160 ratios are elevated in the central ~3.0 kpc region above the nucleus of this galaxy (the "superbubble"). A pixel-by-pixel analysis shows that dust in this region has a higher temperature and/or an emissivity with a steeper spectral index (beta > 2) than the dust in the disk, possibly the result of the harsher environment in the superbubble. Star formation in the disk seems energetically insufficient to lift the material out of the disk, unless it was more active in the past or the dust-to-gas ratio in the superbubble region is higher than the Galactic value. Some of the dust in the halo may also have been tidally stripped from nearby companions or lifted from the disk by galaxy interactions.Comment: Accepted for publication in The Astrophysical Journa

Dust heating by the interstellar radiation field in models of turbulent molecular clouds

We have calculated the radiation field, dust grain temperatures, and far infrared emissivity of numerical models of turbulent molecular clouds. When compared to a uniform cloud of the same mean optical depth, most of the volume inside the turbulent cloud is brighter, but most of the mass is darker. There is little mean attenuation from center to edge, and clumping causes the radiation field to be somewhat bluer. There is also a large dispersion, typically by a few orders of magnitude, of all quantities relative to their means. However, despite the scatter, the 850 micron emission maps are well correlated with surface density. The fraction of mass as a function of intensity can be reproduced by a simple hierarchical model of density structure.Comment: 32 pages, 14 figures, submitted to Ap

Evolutionary multi-stage financial scenario tree generation

Multi-stage financial decision optimization under uncertainty depends on a careful numerical approximation of the underlying stochastic process, which describes the future returns of the selected assets or asset categories. Various approaches towards an optimal generation of discrete-time, discrete-state approximations (represented as scenario trees) have been suggested in the literature. In this paper, a new evolutionary algorithm to create scenario trees for multi-stage financial optimization models will be presented. Numerical results and implementation details conclude the paper

Physical Properties of Complex C Halo Clouds

Observations from the Galactic Arecibo L-Band Feed Array HI (GALFA-HI) Survey of the tail of Complex C are presented and the halo clouds associated with this complex cataloged. The properties of the Complex C clouds are compared to clouds cataloged at the tail of the Magellanic Stream to provide insight into the origin and destruction mechanism of Complex C. Magellanic Stream and Complex C clouds show similarities in their mass distributions (slope = -0.7 and -0.6, respectively) and have a common linewidth of 20 - 30 km/s (indicative of a warm component), which may indicate a common origin and/or physical process breaking down the clouds. The clouds cataloged at the tail of Complex C extend over a mass range of 10^1.1 to 10^4.8 solar masses, sizes of 10^1.2 to 10^2.6 pc, and have a median volume density of 0.065 cm^(-3) and median pressure of (P/k) = 580 K cm^{-3}. We do not see a prominent two-phase structure in Complex C, possibly due to its low metallicity and inefficient cooling compared to other halo clouds. From assuming the Complex C clouds are in pressure equilibrium with a hot halo medium, we find a median halo density of 5.8 x 10^(-4) cm^(-3), which given a constant distance of 10 kpc, is at a z-height of ~3 kpc. Using the same argument for the Stream results in a median halo density of 8.4 x 10^(-5) x (60kpc/d) cm^(-3). These densities are consistent with previous observational constraints and cosmological simulations. We also assess the derived cloud and halo properties with three dimensional grid simulations of halo HI clouds and find the temperature is generally consistent within a factor of 1.5 and the volume densities, pressures and halo densities are consistent within a factor of 3.Comment: Accepted for publication in AJ. 54 pages, including 6 tables and 16 figure

An ammonia spectral map of the L1495-B218 filaments in the Taurus molecular cloud. I. Physical properties of filaments and dense cores

We present deep NH3 observations of the L1495-B218 filaments in the Taurus molecular cloud covering over a 3° angular range using the K-band focal plane array on the 100 m Green Bank Telescope. The L1495-B218 filaments form an interconnected, nearby, large complex extending over 8 pc. We observed NH3 (1, 1) and (2, 2) with a spectral resolution of 0.038 km s−1 and a spatial resolution of 31''. Most of the ammonia peaks coincide with intensity peaks in dust continuum maps at 350 and 500 μm. We deduced physical properties by fitting a model to the observed spectra. We find gas kinetic temperatures of 8–15 K, velocity dispersions of 0.05–0.25 km s−1, and NH3 column densities of 5 × 1012 to 1 × 1014 cm−2. The CSAR algorithm, which is a hybrid of seeded-watershed and binary dendrogram algorithms, identifies a total of 55 NH3 structures, including 39 leaves and 16 branches. The masses of the NH3 sources range from 0.05 to 9.5 ${{M}_{\odot }}$. The masses of NH3 leaves are mostly smaller than their corresponding virial mass estimated from their internal and gravitational energies, which suggests that these leaves are gravitationally unbound structures. Nine out of 39 NH3 leaves are gravitationally bound, and seven out of nine gravitationally bound NH3 leaves are associated with star formation. We also found that 12 out of 30 gravitationally unbound leaves are pressure confined. Our data suggest that a dense core may form as a pressure-confined structure, evolve to a gravitationally bound core, and undergo collapse to form a protostar

The Millisecond Pulsars in NGC 6760

We present the results of recent Arecibo and Green Bank observations of the globular cluster NGC 6760. Using Arecibo, a phase-coherent timing solution has been obtained for the previously known binary pulsar in this cluster, PSR J1911+0102A. We have also discovered a new millisecond pulsar in NGC 6760, PSR J1911+0101B, an isolated object with a rotational period of 5.38 ms and a dispersion measure DM = 196.7 cm-3 pc. Both pulsars are located within 1.3 core radii of the cluster center and have negative period derivatives. The resulting lower limits for the accelerations of the pulsars are within the range expected given a simple model of the cluster. A search for eclipses in the PSR J1911+0102A binary system using both telescopes yielded negative results. The corresponding limits on the extra gas column density at superior conjunction are consistent with the hypothesis that the observational properties of ultra-low-mass binary pulsars like PSR J1911+0102A are strongly affected by the inclination of the orbital plane of the system. Among globular cluster pulsar populations, that of NGC 6760 exhibits one of the largest known spreads in DM. This quantity seems to be roughly proportional to a cluster's central DM; this suggests that the observed spread is caused by a turbulent interstellar medium at spatial scales of 1 pc.Comment: 10 pages in referee format, 4 figures, one table, re-submitted to the Astrophysical Journa

The Frequency of Mid-Infrared Excess Sources in Galactic Surveys

We have identified 230 Tycho-2 Spectral Catalog stars that exhibit 8 micron mid-infrared extraphotospheric excesses in the MidCourse Space Experiment (MSX) and Spitzer Space Telescope Galactic Legacy MidPlane Survey Extraordinaire (GLIMPSE) surveys. Of these, 183 are either OB stars earlier than B8 in which the excess plausibly arises from a thermal bremsstrahlung component or evolved stars in which the excess may be explained by an atmospheric dust component. The remaining 47 stars have spectral classifications B8 or later and appear to be main sequence or late pre-main-sequence objects harboring circumstellar disks. Six of the 47 stars exhibit multiple signatures characteristic of pre-main-sequence circumstellar disks, including emission lines, near-infrared K-band excesses, and X-ray emission. Approximately one-third of the remaining 41 sources have emission lines suggesting relative youth. Of the 25 GLIMPSE stars with SST data at >24 microns, 20 also show an excess at 24 microns. Three additional objects have 24 micron upper limits consistent with possible excesses, and two objects have photospheric measurements at 24 microns. Six MSX sources had a measurement at wavelengths >8 microns. We modeled the excesses in 26 stars having two or more measurements in excess of the expected photospheres as single-component blackbodies. We determine probable disk temperatures and fractional infrared luminosities in the range 191 < T < 787 and 3.9x10^-4 < L_IR/L_* < 2.7x10^-1. We estimate a lower limit on the fraction of Tycho-2 Spectral Catalog main-sequence stars having mid-IR, but not near-IR, excesses to be 1.0+-0.3%.Comment: Accepted to Ap

Scenario trees and policy selection for multistage stochastic programming using machine learning

We propose a hybrid algorithmic strategy for complex stochastic optimization problems, which combines the use of scenario trees from multistage stochastic programming with machine learning techniques for learning a policy in the form of a statistical model, in the context of constrained vector-valued decisions. Such a policy allows one to run out-of-sample simulations over a large number of independent scenarios, and obtain a signal on the quality of the approximation scheme used to solve the multistage stochastic program. We propose to apply this fast simulation technique to choose the best tree from a set of scenario trees. A solution scheme is introduced, where several scenario trees with random branching structure are solved in parallel, and where the tree from which the best policy for the true problem could be learned is ultimately retained. Numerical tests show that excellent trade-offs can be achieved between run times and solution quality

The GALFA-HI Compact Cloud Catalog

We present a catalog of 1964 isolated, compact neutral hydrogen clouds from the Galactic Arecibo L-Band Feed Array Survey Data Release One (GALFA-HI DR1). The clouds were identified by a custom machine-vision algorithm utilizing Difference of Gaussian kernels to search for clouds smaller than 20'. The clouds have velocities typically between |VLSR| = 20-400 km/s, linewidths of 2.5-35 km/s, and column densities ranging from 1 - 35 x 10^18 cm^-2. The distances to the clouds in this catalog may cover several orders of magnitude, so the masses may range from less than a Solar mass for clouds within the Galactic disc, to greater than 10^4 Solar Masses for HVCs at the tip of the Magellanic Stream. To search for trends, we separate the catalog into five populations based on position, velocity, and linewidth: high velocity clouds (HVCs); galaxy candidates; cold low velocity clouds (LVCs); warm, low positive-velocity clouds in the third Galactic Quadrant; and the remaining warm LVCs. The observed HVCs are found to be associated with previously-identified HVC complexes. We do not observe a large population of isolated clouds at high velocities as some models predict. We see evidence for distinct histories at low velocities in detecting populations of clouds corotating with the Galactic disc and a set of clouds that is not corotating.Comment: 34 Pages, 9 Figures, published in ApJ (2012, ApJ, 758, 44), this version has the corrected fluxes and corresponding flux histogram and masse