Global Three-dimensional Magentohydrodynamic Simulations of Galactic Gaseous Disks. I. Amplification of Mean Magnetic Fields in an Axisymmetric Gravitational Potential

We carried out global three-dimensional resistive magnetohydrodynamic simulations of galactic gaseous disks to investigate how the galactic magnetic fields are amplified and maintained. We adopt a steady axisymmetric gravitational potential given by Miyamoto & Nagai and Miyamoto et al. As the initial condition, we assume a warm (T ~ 105 K) rotating gas torus centered at = 10 kpc threaded by weak azimuthal magnetic fields. Numerical results indicate that in differentially rotating galactic gaseous disks, magnetic fields are amplified due to magnetorotational instability and magnetic turbulence develops. After the amplification of magnetic energy saturates, the disk stays in a quasi-steady state. The mean azimuthal magnetic field increases with time and shows reversals with a period of 1 Gyr (2 Gyr for a full cycle). The amplitude of B near the equatorial plane is B ~ 1.5 μG at = 5 kpc. The magnetic fields show large fluctuations whose standard deviation is comparable to the mean field. The mean azimuthal magnetic field in the disk corona has a direction opposite to the mean magnetic field inside the disk. The mass accretion rate driven by the Maxwell stress is ~10-3 M☉ yr-1 at = 2.5 kpc when the mass of the initial torus is ~5 × 108 M☉. When we adopt an absorbing boundary condition at r = 0.8 kpc, the rotation curve obtained by numerical simulations almost coincides with the rotation curve of the stars and the dark matter. Thus, even when magnetic fields are not negligible for gas dynamics of a spiral galaxy, the galactic gravitational potential can be derived from observations of the rotation curve using the gas component of the disk.

[1]  E. Ostriker,et al.  Thermal and Magnetorotational Instability in the Interstellar Medium: Two-dimensional Numerical Simulations , 2003, astro-ph/0310510.

[2]  J. Sellwood,et al.  Differential Rotation and Turbulence in Extended H I Disks , 1998, astro-ph/9806307.

[3]  J. M. Stone,et al.  The Formation and Structure of a Strongly Magnetized Corona above a Weakly Magnetized Accretion Disk , 1999, astro-ph/9912135.

[4]  E. Parker.  The dynamical state of the interstellar gas and field. II. , 1966 .

[5]  C. Norman,et al.  Radio Jets and the Formation of Active Galaxies: Accretion Avalanches on the Torus by the Effect of a Large-Scale Magnetic Field , 1996 .

[6]  David Burstein,et al.  Rotation velocities of 16 SA galaxies and a comparison of Sa, SB and SC rotation properties. , 1985 .

[7]  Guohong Xu,et al.  Parasitic instabilities in magnetized, differentially rotating disks , 1994 .

[8]  K. Shibata,et al.  X-Ray Flares and Mass Outflows Driven by Magnetic Interaction between a Protostar and Its Surrounding Disk , 1996, astro-ph/9606157.

[9]  T. Stepinski,et al.  Generation of dynamo magnetic fields in protoplanetary and other astrophysical accretion disks , 1988 .

[10]  S. Inutsuka,et al.  Saturation and Thermalization of the Magnetorotational Instability: Recurrent Channel Flows and Reconnections , 2001, astro-ph/0110125.

[11]  Global Simulations of Differentially Rotating Magnetized Disks: Formation of Low-beta Filaments and Structured Coronae. , 1999, The Astrophysical journal.

[12]  P. Kronberg,et al.  New large-scale magnetic features of the Milky Way , 1979, Nature.

[13]  Constraining the Magnetic Effects on H I Rotation Curves and the Need for Dark Halos , 2004, astro-ph/0402022.

[14]  E. Parker.  The generation of magnetic fields in astrophysical bodies. II - The galactic field , 1971 .

[15]  S. Mineshige,et al.  ACCEPTED FOR PUBLICATION IN APJ Preprint typeset using LATEX style emulateapj v. 7/8/03 MAGNETOHYDRODYNAMICAL ACCRETION FLOWS: FORMATION OF MAGNETIC TOWER JET AND SUBSEQUENT QUASI-STEADY STATE , 2003 .

[16]  Samuel Z. Burstein,et al.  DIFFERENCE METHODS FOR THE INVISCID AND VISCOUS EQUATIONS OF A COMPRESSIBLE GAS. , 1967 .

[17]  Robert F. Stein,et al.  Dynamo-generated Turbulence and Large-Scale Magnetic Fields in a Keplerian Shear Flow , 1995 .

[18]  M. Merrifield.  The Rotation curve of the milky way to 2.5-R(0) from the thickness of the HI layer , 1992 .

[19]  Ryoji Matsumoto,et al.  Global Three-dimensional Magnetohydrodynamic Simulations of Black Hole Accretion Disks: X-Ray Flares in the Plunging Region , 2002 .

[20]  J. Hawley,et al.  A powerful local shear instability in weakly magnetized disks. IV: Nonaxisymmetric perturbations , 1992 .

[21]  Woong-Tae Kim,et al.  Magnetorotationally Driven Galactic Turbulence and the Formation of Giant Molecular Clouds , 2003, astro-ph/0309080.

[22]  C. Tout,et al.  Accretion disc viscosity: a simple model for a magnetic dynamo. , 1992 .

[23]  E. Battaner,et al.  Magnetic fields as an alternative explanation for the rotation curves of spiral galaxies , 1992, Nature.

[24]  T. Tajima,et al.  Two-dimensional magnetohydrodynamic model of emerging magnetic flux in the solar atmosphere , 1989 .

[25]  Charles F. Gammie,et al.  Local Three-dimensional Magnetohydrodynamic Simulations of Accretion Disks , 1995 .

[26]  A. G. Lyne,et al.  New rotation measures of distant pulsars in the inner Galaxy and magnetic field reversals , 1994 .

[27]  E. Blackman,et al.  Doubly Helical Coronal Ejections from Dynamos and Their Role in Sustaining the Solar Cycle , 2002, astro-ph/0212010.

[28]  E. Parker.  The generation of magnetic fields in astrophysical bodies. X - Magnetic buoyancy and the solar dynamo , 1975 .

[29]  Global Magnetohydrodynamical Simulations of Accretion Tori , 1999, astro-ph/9907385.

[30]  D. Sokoloff,et al.  Galactic Magnetism: Recent developments and perspectives , 1996 .

[31]  John F. Hawley,et al.  A powerful local shear instability in weakly magnetized disks. III - Long-term evolution in a shearing sheet. IV - Nonaxisymmetric perturbations , 1992 .

[32]  E. Spiegel,et al.  Magnetic activity and variations in solar luminosity , 1980, Nature.

[33]  Y. Sofue,et al.  GLOBAL STRUCTURE OF MAGNETIC FIELDS IN SPIRAL GALAXIES , 1986 .

[34]  L. Kitchatinov,et al.  Seed fields for galactic dynamos by the magnetorotational instability , 2004, astro-ph/0405228.

[35]  Ryoji Matsumoto,et al.  Magnetic viscosity by localized shear flow instability in magnetized accretion disks , 1995 .

[36]  Takaaki Yokoyama,et al.  What is the condition for fast magnetic reconnection , 1994 .

[37]  G. Rüdiger,et al.  Interstellar turbulence driven by the magnetorotational instability , 2004 .

[38]  S. Tremaine,et al.  The global properties of the Galaxy. II. The galactic rotation parameters from 21-cm H I observations. , 1979 .

[39]  C. Gammie,et al.  Local three-dimensional simulations of an accretion disk hydromagnetic dynamo , 1996 .

[40]  A. Lyne,et al.  Counterclockwise Magnetic Fields in the Norma Spiral Arm , 2002, astro-ph/0203517.

[41]  R. Rosner,et al.  Doubly diffusive magnetic buoyancy instability in the solar interior , 1983 .

[42]  J. Hawley,et al.  A powerful local shear instability in weakly magnetized disks. I - Linear analysis. II - Nonlinear evolution , 1990 .