Magneto-Centrifugal Launching of Jets from Accretion Disks. I: Cold Axisymmetric Flows

The magneto-centrifugal model for jet formation is studied by time-dependent simulations reaching steady state in a cold gas with negligible fluid pressure, in an axisymmetric geometry, using a modification of the Zeus3D code adapted to parallel computers. The number of boundary conditions imposed at the coronal base takes into account the existence of the fast and Alfvenic critical surfaces, avoiding over-determination of the flow. The size and shape of the computational box is chosen to include these critical surfaces, reducing the influence of the outer boundary conditions. As there is a region, near the origin, where the inclination of field lines to the axis is too small to drive a centrifugal wind, we inject a thin, axial jet, expected to form electromagnetically near black holes. Acceleration and collimation appear for wide generic conditions. A reference run is shown in detail, with a wind leaving the computational volume in the axial direction with a poloidal velocity equal to 4 times the poloidal Alfven speed, collimated inside 11 degrees. Finally, the critical surfaces, fieldlines, thrust, energy, torque and mass discharge of the outgoing wind are shown for simulations with various profiles of mass and magnetic flux at the base of the corona.Comment: 27 pages, including 10 figures and 2 tables. To appear in ApJ (Dec 1999). Revised version clarifies the abstract, section 3.2.4, conclusions and appendix, adds a simulation to section 4.2, and updates the reference

On the Link between Martian Total Ozone and Potential Vorticity

AbstractWe demonstrate for the first time that total ozone in the martian atmosphere is highly correlated with the dynamical tracer, potential vorticity, under certain conditions. The degree of correlation is investigated using a Mars global circulation model including a photochemical model. Potential vorticity is the quantity of choice to explore the dynamical nature of polar vortices because it contains information on winds and temperature in a single scalar variable. The correlation is found to display a distinct seasonal variation, with a strong positive correlation in both northern and southern winter at poleward latitudes in the northern and southern hemisphere respectively.The identified strong correlation implies variations in polar total ozone during winter are predominantly controlled by dynamical processes in these spatio-temporal regions. The weak correlation in northern and southern summer is due to the dominance of photochemical reactions resulting from extended exposure to sunlight. The total ozone/potential vorticity correlation is slightly weaker in southern winter due to topographical variations and the preference for ozone to accumulate in Hellas basin. In northern winter, total ozone can be used to track the polar vortex edge.The ozone/potential vorticity ratio is calculated for both northern and southern winter on Mars for the first time. Using the strong correlation in total ozone and potential vorticity in northern winter inside the polar vortex, it is shown that potential vorticity can be used as a proxy to deduce the distribution of total ozone where satellites cannot observe for the majority of northern winter. Where total ozone observations are available on the fringes of northern winter at poleward latitudes, the strong relationship of total ozone and potential vorticity implies that total ozone anomalies in the surf zone of the northern polar vortex can potentially be used to determine the origin of potential vorticity filaments

Magnetocentrifugal Winds in 3D: Nonaxisymmetric Steady State

Outflows can be loaded and accelerated to high speeds along rapidly rotating, open magnetic field lines by centrifugal forces. Whether such magnetocentrifugally driven winds are stable is a longstanding theoretical problem. As a step towards addressing this problem, we perform the first large-scale 3D MHD simulations that extend to a distance $\sim 10^2$ times beyond the launching region, starting from steady 2D (axisymmetric) solutions. In an attempt to drive the wind unstable, we increase the mass loading on one half of the launching surface by a factor of $\sqrt{10}$, and reduce it by the same factor on the other half. The evolution of the perturbed wind is followed numerically. We find no evidence for any rapidly growing instability that could disrupt the wind during the launching and initial phase of propagation, even when the magnetic field of the magnetocentrifugal wind is toroidally dominated all the way to the launching surface. The strongly perturbed wind settles into a new steady state, with a highly asymmetric mass distribution. The distribution of magnetic field strength is, in contrast, much more symmetric. We discuss possible reasons for the apparent stability, including stabilization by an axial poloidal magnetic field, which is required to bend field lines away from the vertical direction and produce a magnetocentrifugal wind in the first place.Comment: 10 pages, 2 figures, accepted for publication in ApJ

Spectral Analysis of the Chandra Comet Survey

We present results of the analysis of cometary X-ray spectra with an extended version of our charge exchange emission model (Bodewits et al. 2006). We have applied this model to the sample of 8 comets thus far observed with the Chandra X-ray observatory and ACIS spectrometer in the 300-1000 eV range. The surveyed comets are C/1999 S4 (LINEAR), C/1999 T1 (McNaught-Hartley), C/2000 WM1 (LINEAR), 153P/2002 (Ikeya-Zhang), 2P/2003 (Encke), C/2001 Q4 (NEAT), 9P/2005 (Tempel 1) and 73P/2006-B (Schwassmann-Wachmann 3) and the observations include a broad variety of comets, solar wind environments and observational conditions. The interaction model is based on state selective, velocity dependent charge exchange cross sections and is used to explore how cometary X-ray emission depend on cometary, observational and solar wind characteristics. It is further demonstrated that cometary X-ray spectra mainly reflect the state of the local solar wind. The current sample of Chandra observations was fit using the constrains of the charge exchange model, and relative solar wind abundances were derived from the X-ray spectra. Our analysis showed that spectral differences can be ascribed to different solar wind states, as such identifying comets interacting with (I) fast, cold wind, (II), slow, warm wind and (III) disturbed, fast, hot winds associated with interplanetary coronal mass ejections. We furthermore predict the existence of a fourth spectral class, associated with the cool, fast high latitude wind.Comment: 16 pages, 16 figures, and 7 Tables; accepted A&A (Due to space limits, this version has lower resolution jpeg images.

A reanalysis of ozone on Mars from assimilation of SPICAM observations

We have assimilated for the first time SPICAM retrievals of total ozone into a Martian global circulation model to provide a global reanalysis of the ozone cycle. Disagreement in total ozone between model prediction and assimilation is observed between 45°S–10°S from LS=135–180° and at northern polar (60°N–90°N) latitudes during northern fall (LS=150–195°). Large percentage differences in total ozone at northern fall polar latitudes identified through the assimilation process are linked with excessive northward transport of water vapour west of Tharsis and over Arabia Terra. Modelling biases in water vapour can also explain the underestimation of total ozone between 45°S–10°S from LS=135–180°. Heterogeneous uptake of odd hydrogen radicals are unable to explain the outstanding underestimation of northern polar total ozone in late northern fall. Assimilation of total ozone retrievals results in alterations of the modelled spatial distribution of ozone in the southern polar winter high altitude ozone layer. This illustrates the potential use of assimilation methods in constraining total ozone where SPICAM cannot observe, in a region where total ozone is especially important for potential investigations of the polar dynamics

Disc formation in turbulent massive cores: Circumventing the magnetic braking catastrophe

We present collapse simulations of 100 M_{\sun}, turbulent cloud cores threaded by a strong magnetic field. During the initial collapse phase filaments are generated which fragment quickly and form several protostars. Around these protostars Keplerian discs with typical sizes of up to 100 AU build up in contrast to previous simulations neglecting turbulence. We examine three mechanisms potentially responsible for lowering the magnetic braking efficiency and therefore allowing for the formation of Keplerian discs. Analysing the condensations in which the discs form, we show that the build-up of Keplerian discs is neither caused by magnetic flux loss due to turbulent reconnection nor by the misalignment of the magnetic field and the angular momentum. It is rather a consequence of the turbulent surroundings of the disc which exhibit no coherent rotation structure while strong local shear flows carry large amounts of angular momentum. We suggest that the "magnetic braking catastrophe", i.e. the formation of sub-Keplerian discs only, is an artefact of the idealised non-turbulent initial conditions and that turbulence provides a natural mechanism to circumvent this problem.Comment: 6 pages, 5 figures, accepted by MNRAS Letters, updated to final versio

Magneto-Centrifugal Launching of Jets from Accretion Disks. II: Inner Disk-Driven Winds

We follow numerically the time evolution of axisymmetric outflows driven magneto-centrifugally from the inner portion of accretion disks, from their launching surface to large, observable distances. Special attention is paid to the collimation of part of the outflow into a dense, narrow jet around the rotation axis, after a steady state has been reached. For parameters typical of T Tauri stars, we define a fiducial ``jet'' as outlined by the contour of constant density at 10^4 cm^{-3}. We find that the jet, so defined, appears nearly cylindrical well above the disk, in agreement with previous asymptotic analyses. Closer to the equatorial plane, the density contour can either bulge outwards or pinch inwards, depending on the conditions at the launching surface, particularly the mass flux distribution. We find that even though a dense, jet-like feature is always formed around the axis, there is no guarantee that the high-density axial jet would dominate the more tenuous, wide-angle part of the wind. Specifically, on the 100 AU scale, resolvable by HST and ground-based adaptive optics for nearby T Tauri winds, the fraction of the wind mass flux enclosed by the fiducial jet can vary substantially, again depending on the launching conditions. We show two examples in which the fraction is ~20% and ~45%. These dependences may provide a way to constrain the conditions at the launching surface, which are poorly known at present.Comment: 11 pages, 6 figures. Accepted for publication in ApJ, scheduled for vol. 595, October 1, 200

Are Magnetic Wind-Driving Disks Inherently Unstable?

There have been claims in the literature that accretion disks in which a centrifugally driven wind is the dominant mode of angular momentum transport are inherently unstable. This issue is considered here by applying an equilibrium-curve analysis to the wind-driving, ambipolar diffusion-dominated, magnetic disk model of Wardle & Konigl (1993). The equilibrium solution curves for this class of models typically exhibit two distinct branches. It is argued that only one of these branches represents unstable equilibria and that a real disk/wind system likely corresponds to a stable solution.Comment: 5 pages, 2 figures, to be published in ApJ, vol. 617 (2004 Dec 20). Uses emulateapj.cl

Incorporating Ambipolar and Ohmic Diffusion in the AMR MHD code RAMSES

We have implemented non-ideal Magneto-Hydrodynamics (MHD) effects in the Adaptive Mesh Refinement (AMR) code RAMSES, namely ambipolar diffusion and Ohmic dissipation, as additional source terms in the ideal MHD equations. We describe in details how we have discretized these terms using the adaptive Cartesian mesh, and how the time step is diminished with respect to the ideal case, in order to perform a stable time integration. We have performed a large suite of test runs, featuring the Barenblatt diffusion test, the Ohmic diffusion test, the C-shock test and the Alfven wave test. For the latter, we have performed a careful truncation error analysis to estimate the magnitude of the numerical diffusion induced by our Godunov scheme, allowing us to estimate the spatial resolution that is required to address non-ideal MHD effects reliably. We show that our scheme is second-order accurate, and is therefore ideally suited to study non-ideal MHD effects in the context of star formation and molecular cloud dynamics