The Effect of the Hall Term on the Nonlinear Evolution of the Magnetorotational Instability: II. Saturation Level and Critical Magnetic Reynolds Number

The nonlinear evolution of the magnetorotational instability (MRI) in weakly ionized accretion disks, including the effect of the Hall term and ohmic dissipation, is investigated using local three-dimensional MHD simulations and various initial magnetic field geometries. When the magnetic Reynolds number, Re_M \equiv v_A^2 / \eta \Omega (where v_A is the Alfven speed, \eta the magnetic diffusivity, and \Omega the angular frequency), is initially larger than a critical value Re_{M, crit}, the MRI evolves into MHD turbulence in which angular momentum is transported efficiently by the Maxwell stress. If Re_M < Re_{M, crit}, however, ohmic dissipation suppresses the MRI, and the stress is reduced by several orders of magnitude. The critical value is in the range of 1 - 30 depending on the initial field configuration. The Hall effect does not modify the critical magnetic Reynolds number by much, but enhances the saturation level of the Maxwell stress by a factor of a few. We show that the saturation level of the MRI is characterized by v_{Az}^2 / \eta \Omega, where v_{Az} is the Alfven speed in the nonlinear regime along the vertical component of the field. The condition for turbulence and significant transport is given by v_{Az}^2 / \eta \Omega \gtrsim 1, and this critical value is independent of the strength and geometry of the magnetic field or the size of the Hall term. If the magnetic field strength in an accretion disk can be estimated observationally, and the magnetic Reynolds number v_A^2 / \eta \Omega is larger than about 30, this would imply the MRI is operating in the disk.Comment: 43 pages, 8 tables, 20 figures, accepted for publication in ApJ, postscript version also available from http://www.astro.umd.edu/~sano/publications

Axisymmetric Magnetorotational Instability in Viscous Accretion Disks

Axisymmetric magnetorotational instability (MRI) in viscous accretion disks is investigated by linear analysis and two-dimensional nonlinear simulations. The linear growth of the viscous MRI is characterized by the Reynolds number defined as $R_{\rm MRI} \equiv v_A^2/\nu\Omega$, where $v_A$ is the Alfv{\'e}n velocity, $\nu$ is the kinematic viscosity, and $\Omega$ is the angular velocity of the disk. Although the linear growth rate is suppressed considerably as the Reynolds number decreases, the nonlinear behavior is found to be almost independent of $R_{\rm MRI}$. At the nonlinear evolutionary stage, a two-channel flow continues growing and the Maxwell stress increases until the end of calculations even though the Reynolds number is much smaller than unity. A large portion of the injected energy to the system is converted to the magnetic energy. The gain rate of the thermal energy, on the other hand, is found to be much larger than the viscous heating rate. Nonlinear behavior of the MRI in the viscous regime and its difference from that in the highly resistive regime can be explained schematically by using the characteristics of the linear dispersion relation. Applying our results to the case with both the viscosity and resistivity, it is anticipated that the critical value of the Lundquist number $S_{\rm MRI} \equiv v_A^2/\eta\Omega$ for active turbulence depends on the magnetic Prandtl number $S_{{\rm MRI},c} \propto Pm^{1/2}$ in the regime of $Pm \gg 1$ and remains constant when $Pm \ll 1$, where $Pm \equiv S_{\rm MRI}/R_{\rm MRI} = \nu/\eta$ and $\eta$ is the magnetic diffusivity.Comment: Accepted for publication in ApJ -- 18 pages, 9 figures, 1 tabl

A Local One-Zone Model of MHD Turbulence in Dwarf Nova Disks

The evolution of the magnetorotational instability (MRI) during the transition from outburst to quiescence in a dwarf nova disk is investigated using three-dimensional MHD simulations. The shearing box approximation is adopted for the analysis, so that the efficiency of angular momentum transport is studied in a small local patch of the disk: this is usually referred as to a one-zone model. To take account of the low ionization fraction of the disk, the induction equation includes both ohmic dissipation and the Hall effect. We induce a transition from outburst to quiescence by an instantaneous decrease of the temperature. The evolution of the MRI during the transition is found to be very sensitive to the temperature of the quiescent disk. As long as the temperature is higher than a critical value of about 2000 K, MHD turbulence and angular momentum transport is sustained by the MRI. However, MHD turbulence dies away within an orbital time if the temperature falls below this critical value. In this case, the stress drops off by more than 2 orders of magnitude, and is dominated by the Reynolds stress associated with the remnant motions from the outburst. The critical temperature depends slightly on the distance from the central star and the local density of the disk.Comment: 20 pages, 2 tables, 6 figures, accepted for publication in Ap

Electronic States and Superconducting Transition Temperature based on the Tomonaga-Luttinger liquid in Pr2_{2}Ba4_{4}Cu7_{7}O15−δ_{15-\delta}

An NQR experiment revealed superconductivity of Pr$_2$Ba$_4$Cu$_7$O$_{15-\delta}$ (Pr247) to be realized on CuO double chain layers and suggests possibility of novel one-dimensional(1D) superconductivity. To clarify the nature of the 1D superconductivity, we calculate the band dispersions of Pr247 by using the generalized gradient approximation(GGA). It indicates that Fermi surface of CuO double chains is well described to the electronic structure of a quasi-1D system. Assuming the zigzag Hubbard chain model to be an effective model of the system, we derive tight binding parameters of the model from a fit to the result of GGA. Based on the Tomonaga-Luttinger liquid theory, we estimate transition temperature ($T_c$) of the quasi-1D zigzag Hubbard model from the calculated value of the Luttinger liquid parameter $K_{\rho}$. The result of $T_c$ is consistent with that of experiments in Pr247 and it suggests that the mechanism of the superconductivity is well understood within the concept of the Tomonaga-Luttinger liquid.Comment: 4 pages, 5 figure

The Effect of the Hall Term on the Nonlinear Evolution of the Magnetorotational Instability: I. Local Axisymmetric Simulations

The effect of the Hall term on the evolution of the magnetorotational instability (MRI) in weakly ionized accretion disks is investigated using local axisymmetric simulations. First, we show that the Hall term has important effects on the MRI when the temperature and density in the disk is below a few thousand K and between 10^13 and 10^18 cm^{-3} respectively. Such conditions can occur in the quiescent phase of dwarf nova disks, or in the inner part (inside 10 - 100 AU) of protoplanetary disks. When the Hall term is important, the properties of the MRI are dependent on the direction of the magnetic field with respect to the angular velocity vector \Omega. If the disk is threaded by a uniform vertical field oriented in the same sense as \Omega, the axisymmetric evolution of the MRI is an exponentially growing two-channel flow without saturation. When the field is oppositely directed to \Omega, however, small scale fluctuations prevent the nonlinear growth of the channel flow and the MRI evolves into MHD turbulence. These results are anticipated from the characteristics of the linear dispersion relation. In axisymmetry on a field with zero-net flux, the evolution of the MRI is independent of the size of the Hall term relative to the inductive term. The evolution in this case is determined mostly by the effect of ohmic dissipation.Comment: 31 pages, 3 tables, 12 figures, accepted for publication in ApJ, postscript version also available from http://www.astro.umd.edu/~sano/publications

Ferromagnetism and Superconductivity in the multi-orbital Hubbard Model: Hund's Rule Coupling versus Crystal-Field Splitting

The multi-orbital Hubbard model in one dimension is studied using the numerical diagonalization method. Due to the effect of the crystal-field splitting $\Delta$, the fully polarized ferromagnetism which is observed in the strong coupling regime becomes unstable against the partially polarized ferromagnetism when the Hund's rule coupling $J$ is smaller than a certain critical value of order of $\Delta$. In the vicinity of the partially polarized ferromagnetism, the orbital fluctuation develops due to the competition between the Hund's rule coupling and the crystal-field splitting. The superconducting phase with the Luttinger liquid parameter $K_{\rho}>1$ is observed for the singlet ground state in this region.Comment: 4 pages,5 figures,submitted to J.Phys.Soc.Jp

Angular Momentum Transport by MHD Turbulence in Accretion Disks: Gas Pressure Dependence of the Saturation Level of the Magnetorotational Instability

The saturation level of the magnetorotational instability (MRI) is investigated using three-dimensional MHD simulations. The shearing box approximation is adopted and the vertical component of gravity is ignored, so that the evolution of the MRI is followed in a small local part of the disk. We focus on the dependence of the saturation level of the stress on the gas pressure, which is a key assumption in the standard alpha disk model. From our numerical experiments it is found that there is a weak power-law relation between the saturation level of the Maxwell stress and the gas pressure in the nonlinear regime; the higher the gas pressure, the larger the stress. Although the power-law index depends slightly on the initial field geometry, the relationship between stress and gas pressure is independent of the initial field strength, and is unaffected by Ohmic dissipation if the magnetic Reynolds number is at least 10. The relationship is the same in adiabatic calculations, where pressure increases over time, and nearly-isothermal calculations, where pressure varies little with time. Our numerical results are qualitatively consistent with an idea that the saturation level of the MRI is determined by a balance between the growth of the MRI and the dissipation of the field through reconnection. The quantitative interpretation of the pressure-stress relation, however, may require advances in the theoretical understanding of non-steady magnetic reconnection.Comment: 45 pages, 5 tables, 17 figures, accepted for publication in Ap

Superconductivity in a Two-Orbital Hubbard Model with Electron and Hole Fermi Pockets: Application in Iron Oxypnictide Superconductors

We investigate the electronic states of a one-dimensional two-orbital Hubbard model with band splitting by the exact diagonalization method. The Luttinger liquid parameter $K_{\rho}$ is calculated to obtain superconducting (SC) phase diagram as a function of on-site interactions: the intra- and inter-orbital Coulomb $U$ and $U'$, the Hund coupling $J$, and the pair transfer $J'$. In this model, electron and hole Fermi pockets are produced when the Fermi level crosses both the upper and lower orbital bands. We find that the system shows two types of SC phases, the SC \Roman{u'-large} for $U>U'$ and the SC \Roman{u-large} for $U<U'$, in the wide parameter region including both weak and strong correlation regimes. Pairing correlation functions indicate that the most dominant pairing for the SC \Roman{u'-large} (SC \Roman{u-large}) is the intersite (on-site) intraorbital spin-singlet with (without) sign reversal of the order parameters between two Fermi pockets. The result of the SC \Roman{u'-large} is consistent with the sign-reversing s-wave pairing that has recently been proposed for iron oxypnictide superconductors.Comment: 5 pages, 8 figures, to appear in J. Phys. Soc. Jpn., Vol.78, No.12, p.12470

Modeling the gamma-ray emission produced by runaway cosmic rays in the environment of RX J1713.7-3946

Diffusive shock acceleration in supernova remnants is the most widely invoked paradigm to explain the Galactic cosmic ray spectrum. Cosmic rays escaping supernova remnants diffuse in the interstellar medium and collide with the ambient atomic and molecular gas. From such collisions gamma-rays are created, which can possibly provide the first evidence of a parent population of runaway cosmic rays. We present model predictions for the GeV to TeV gamma-ray emission produced by the collisions of runaway cosmic rays with the gas in the environment surrounding the shell-type supernova remnant RX J1713.7-3946. The spectral and spatial distributions of the emission, which depend upon the source age, the source injection history, the diffusion regime and the distribution of the ambient gas, as mapped by the LAB and NANTEN surveys, are studied in detail. In particular, we find for the region surrounding RX J1713-3946, that depending on the energy one is observing at, one may observe startlingly different spectra or may not detect any enhanced emission with respect to the diffuse emission contributed by background cosmic rays. This result has important implications for current and future gamma-ray experiments.Comment: version published on PAS

Molecular Clouds as Cosmic-Ray Barometers

The advent of high sensitivity, high resolution gamma-ray detectors, together with a knowledge of the distribution of the atomic hydrogen and especially of the molecular hydrogen in the Galaxy on sub-degree scales creates a unique opportunity to explore the flux of cosmic rays in the Galaxy. We here present the new data on the distribution of the molecular hydrogen from a large region of the inner Galaxy obtained by the NANTEN Collaboration. We then introduce a methodology which aims to provide a test bed for current and future gamma-ray observatories to explore the cosmic ray flux at various positions in our Galaxy. In particular, for a distribution of molecular clouds, as provided by the NANTEN survey, and local cosmic ray density as measured at the Earth, we estimate the expected GeV to TeV gamma-ray signal, which can then be compared with observations and use to test the cosmic ray flux.Comment: PASJ (in press