Triggered Star Formation in the Environment of Young Massive Stars

Recent observations with the Spitzer Space Telescope show clear evidence that star formation takes place in the surrounding of young massive O-type stars, which are shaping their environment due to their powerful radiation and stellar winds. In this work we investigate the effect of ionising radiation of massive stars on the ambient interstellar medium (ISM): In particular we want to examine whether the UV-radiation of O-type stars can lead to the observed pillar-like structures and can trigger star formation. We developed a new implementation, based on a parallel Smooth Particle Hydrodynamics code (called IVINE), that allows an efficient treatment of the effect of ionising radiation from massive stars on their turbulent gaseous environment. Here we present first results at very high resolution. We show that ionising radiation can trigger the collapse of an otherwise stable molecular cloud. The arising structures resemble observed structures (e.g. the pillars of creation in the Eagle Nebula (M16) or the Horsehead Nebula B33). Including the effect of gravitation we find small regions that can be identified as formation places of individual stars. We conclude that ionising radiation from massive stars alone can trigger substantial star formation in molecular clouds.Comment: 4 pages, 2 figures. To appear in: "Triggered Star Formation in a Turbulent ISM", IAU Symposium 237, Prague, Czech Republic, August 2006; eds. B.G.Elmegreen & J. Palou

Gravitational Collapse in Turbulent Molecular Clouds. II. Magnetohydrodynamical Turbulence

Hydrodynamic supersonic turbulence can only prevent local gravitational collapse if the turbulence is driven on scales smaller than the local Jeans lengths in the densest regions, a very severe requirement (Paper I). Magnetic fields have been suggested to support molecular clouds either magnetostatically or via magnetohydrodynamic (MHD) waves. Whereas the first mechanism would form sheet-like clouds, the second mechanism not only could exert a pressure onto the gas counteracting the gravitational forces, but could lead to a transfer of turbulent kinetic energy down to smaller spatial scales via MHD wave interactions. This turbulent magnetic cascade might provide sufficient energy at small scales to halt local collapse. We test this hypothesis with MHD simulations at resolutions up to 256^3 zones, done with ZEUS-3D. We first derive a resolution criterion for self-gravitating, magnetized gas: in order to prevent collapse of magnetostatically supported regions due to numerical diffusion, the minimum Jeans length must be resolved by four zones. Resolution of MHD waves increases this requirement to roughly six zones. We then find that magnetic fields cannot prevent local collapse unless they provide magnetostatic support. Weaker magnetic fields do somewhat delay collapse and cause it to occur more uniformly across the supported region in comparison to the hydrodynamical case. However, they still cannot prevent local collapse for much longer than a global free-fall time.Comment: 32 pages, 14 figures, accepted by Ap

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

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

Numerical Tests of Fast Reconnection in Weakly Stochastic Magnetic Fields

We study the effects of turbulence on magnetic reconnection using 3D numerical simulations. This is the first attempt to test a model of fast magnetic reconnection in the presence of weak turbulence proposed by Lazarian & Vishniac (1999). This model predicts that weak turbulence, generically present in most of astrophysical systems, enhances the rate of reconnection by reducing the transverse scale for reconnection events and by allowing many independent flux reconnection events to occur simultaneously. As a result the reconnection speed becomes independent of Ohmic resistivity and is determined by the magnetic field wandering induced by turbulence. To quantify the reconnection speed we use both an intuitive definition, i.e. the speed of the reconnected flux inflow, as well as a more sophisticated definition based on a formally derived analytical expression. Our results confirm the predictions of the Lazarian & Vishniac model. In particular, we find that Vrec Pinj^(1/2), as predicted by the model. The dependence on the injection scale for some of our models is a bit weaker than expected, i.e. l^(3/4), compared to the predicted linear dependence on the injection scale, which may require some refinement of the model or may be due to the effects like finite size of the excitation region. The reconnection speed was found to depend on the expected rate of magnetic field wandering and not on the magnitude of the guide field. In our models, we see no dependence on the guide field when its strength is comparable to the reconnected component. More importantly, while in the absence of turbulence we successfully reproduce the Sweet-Parker scaling of reconnection, in the presence of turbulence we do not observe any dependence on Ohmic resistivity, confirming that our reconnection is fast.Comment: 22 pages, 20 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

IVINE - Ionization in the parallel tree/sph code VINE: First results on the observed age-spread around O-stars

We present a three-dimensional, fully parallelized, efficient implementation of ionizing ultraviolet (UV) radiation for smoothed particle hydrodynamics (sph) including self-gravity. Our method is based on the sph/tree code vine. We therefore call it iVINE (for Ionization + VINE). This approach allows detailed high-resolution studies of the effects of ionizing radiation from, for example, young massive stars on their turbulent parental molecular clouds. In this paper, we describe the concept and the numerical implementation of the radiative transfer for a plane-parallel geometry and we discuss several test cases demonstrating the efficiency and accuracy of the new method. As a first application, we study the radiatively driven implosion of marginally stable molecular clouds at various distances of a strong UV source and show that they are driven into gravitational collapse. The resulting cores are very compact and dense exactly as it is observed in clustered environments. Our simulations indicate that the time of triggered collapse depends on the distance of the core from the UV source. Clouds closer to the source collapse several 105 yr earlier than more distant clouds. This effect can explain the observed age spread in OB associations where stars closer to the source are found to be younger. We discuss possible uncertainties in the observational derivation of shock front velocities due to early stripping of protostellar envelopes by ionizing radiation

Turbulent Control of the Star Formation Efficiency

Supersonic turbulence plays a dual role in molecular clouds: On one hand, it contributes to the global support of the clouds, while on the other it promotes the formation of small-scale density fluctuations, identifiable with clumps and cores. Within these, the local Jeans length \Ljc is reduced, and collapse ensues if \Ljc becomes smaller than the clump size and the magnetic support is insufficient (i.e., the core is ``magnetically supercritical''); otherwise, the clumps do not collapse and are expected to re-expand and disperse on a few free-fall times. This case may correspond to a fraction of the observed starless cores. The star formation efficiency (SFE, the fraction of the cloud's mass that ends up in collapsed objects) is smaller than unity because the mass contained in collapsing clumps is smaller than the total cloud mass. However, in non-magnetic numerical simulations with realistic Mach numbers and turbulence driving scales, the SFE is still larger than observational estimates. The presence of a magnetic field, even if magnetically supercritical, appears to further reduce the SFE, but by reducing the probability of core formation rather than by delaying the collapse of individual cores, as was formerly thought. Precise quantification of these effects as a function of global cloud parameters is still needed.Comment: Invited review for the conference "IMF@50: the Initial Mass Function 50 Years Later", to be published by Kluwer Academic Publishers, eds. E. Corbelli, F. Palla, and H. Zinnecke

Physics of the Galactic Center Cloud G2, on its Way towards the Super-Massive Black Hole

The origin, structure and evolution of the small gas cloud, G2, is investigated, that is on an orbit almost straight into the Galactic central supermassive black hole (SMBH). G2 is a sensitive probe of the hot accretion zone of Sgr A*, requiring gas temperatures and densities that agree well with models of captured shock-heated stellar winds. Its mass is equal to the critical mass below which cold clumps would be destroyed quickly by evaporation. Its mass is also constrained by the fact that at apocenter its sound crossing timescale was equal to its orbital timescale. Our numerical simulations show that the observed structure and evolution of G2 can be well reproduced if it formed in pressure equilibrium with the surrounding in 1995 at a distance from the SMBH of 7.6e16 cm. If the cloud would have formed at apocenter in the 'clockwise' stellar disk as expected from its orbit, it would be torn into a very elongated spaghetti-like filament by 2011 which is not observed. This problem can be solved if G2 is the head of a larger, shell-like structure that formed at apocenter. Our numerical simulations show that this scenario explains not only G2's observed kinematical and geometrical properties but also the Br_gamma observations of a low surface brightness gas tail that trails the cloud. In 2013, while passing the SMBH G2 will break up into a string of droplets that within the next 30 years mix with the surrounding hot gas and trigger cycles of AGN activity.Comment: 22 pages, 13 figures, submitted to Ap

Fast Reconnection in a Two-Stage Process

Magnetic reconnection plays an essential role in the generation and evolution of astrophysical magnetic fields. The best tested and most robust reconnection theory is that of Parker and Sweet. According to this theory, the reconnection rate scales with magnetic diffusivity lambda as lambda^0.5. In the interstellar medium, the Parker-Sweet reconnection rate is far too slow to be of interest. Thus, a mechanism for fast reconnection seems to be required. We have studied the magnetic merging of two oppositely directed flux systems in weakly ionized, but highly conducting, compressible gas. In such systems, ambipolar diffusion steepens the magnetic profile, leading to a thin current sheet. If the ion pressure is small enough, and the recombination of ions is fast enough, the resulting rate of magnetic merging is fast, and independent of lambda. Slow recombination or sufficiently large ion pressure leads to slower merging which scales with lambda as lambda^0.5. We derive a criterion for distinguishing these two regimes, and discuss applications to the weakly ionized ISM and to protoplanetary accretion disks.Comment: 21 pages, 13 figures, submitted to Ap