MHD simulations of the long-term evolution of a dipolar magnetosphere surrounded by an accretion disk

The evolution of a stellar, initially dipole type magnetosphere interacting with an accretion disk is investigated using numerical ideal MHD simulations. The simulations follow several 1000 Keplerian periods of the inner disk (for animated movies see http://www.aip.de～cfendt).Our model prescribes a Keplerian disk around a rotating star as a fixed boundary condition. The initial magnetic field distribution remains frozen into the star and the disk. The mass flow rate into the corona is fixed for both components. The initial dipole type magnetic field develops into a spherically radial outflow pattern with two main components – a disk wind and a stellar wind – both evolving into a quasi-stationary final state. A neutral field line divides both components, along which small plasmoids are ejected in irregular time intervals. The half opening angle of the stellar wind cone varies from 30° to55° depending on the ratio of the mass flow rates of disk wind and stellar wind. The maximum speed of the outflow is about the Keplerian speed at the inner disk radius. An axial jet forms during the first decades of rotations. However, this feature does not survive on the very long time scale and a pressure driven low velocity flow along the axis evolves. Within a cone of 15° along the axis the formation of knots may be observed if the stellar wind is weak. With the chosen mass flow rates and field strength we see almost no indication for a flow self-collimation. This is due to the weak net poloidal electric current in the magnetosphere which is in difference to typical jet models.     

Numerical simulations of dissipationless disk accretion

Our goal is to study the regime of disk accretion in which almost all of the angular momentum and energy is carried away by the wind outflowing from the disk in numerical experiments. For this type of accretion the kinetic energy flux in the outflowing wind can exceed considerably the bolometric luminosity of the accretion disk, what is observed in the plasma flow from galactic nuclei in a number of cases. In this paper we consider the nonrelativistic case of an outflow from a cold Keplerian disk. All of the conclusions derived previously for such a system in the self-similar approximation are shown to be correct. The numerical results agree well with the analytical predictions. The inclination angle of the magnetic field lines in the disk is less than 60°, which ensures a free wind outflow from the disk, while the energy flux per wind particle is greater than the particle rotation energy in its Keplerian orbit by several orders of magnitude, provided that the ratio r _A/r ? 1, where r _A is the Alfvénic radius and r is the radius of the Keplerian orbit. In this case, the particle kinetic energy reaches half the maximum possible energy in the simulation region. The magnetic field collimates the outflowing wind near the rotation axis and decollimates appreciably the wind outflowing from the outer disk periphery.     

The equilibrium,stability and evolution of a rotating magnetized gaseous disk

The exact solutions for the equilibrium of rotating gaseous disk with poloidal magnetic field are obtained. The stability of the disk with respect to uniform expansion and contraction is investigated by means of the variational principle. It is shown that if the equilibrium is determined by gravitational and magnetic forces only, the disk is in neutral equilibrium with respect to perturbations of the form r=r. The instability to short-waves perturbations is studied by the quasi-classical method. The analysis shows that if the magnetic field isH>2G, where is the surface density, then these perturbations are stabilized. The configurations of the electrical field induced by the rotation of magnetized disk are found. In conclusion, the questions of the evolution of the disk are discussed in connection with the quasar model when pulsar-like radiation is taken into account.     

Three-Dimensional Simulations of Jets from Keplerian Disks: Stability Issues

Three-dimensional simulations of the time-dependent evolution of non-relativistic outflows from the surface of Keplerian accretion disks are presented. We investigate the outflow that arises from a magnetized accretion disk that is initially in hydrostatic balance with its surrounding cold corona. Our simulations show that jets maintain their long-term stability through a self-limiting process wherein the average Alfvénic Mach number within the jet is maintained to order unity. This is accomplished in at least two ways. First, poloidal magnetic field is concentrated along the central axis of the jet forming a `backbone' in which the Alfvén speed is sufficiently high to reduce the average jet Alfvénic Mach number to unity. Second, the onset of higher order Kelvin-Helmholtz `flute' modes (m ≥ 2) reduce the efficiency with which the jet material is accelerated, and transfer kinetic energy of the out flow into the stretched, poloidal field lines of the distorted jet. This too has the effect of increasing the Alfvén speed and thereby reducing the Alfvénic Mach number. The jet is able to survive the onset of the more destructive m=1 mode in this way.     

Average photospheric poloidal and toroidal magnetic field components near solar minimum

Average (over longitude and time) photospheric magnetic field components are derived from 3 Stanford magnetograms made near the solar minimum of cycle 21. The average magnetograph signal is found to behave as the projection of a vector for measurements made across the disk. The poloidal field exhibits the familiar dipolar structure near the poles, with a measured signal in the line Fe i 5250 Å of 1 G. At low latitudes the poloidal field has the polarity of the poles, but is of reduced magnitude ( 0.1 G). A net photospheric toroidal field with a broad latitudinal extent is found. The polarity of the toroidal field is opposite in the nothern and southern hemispheres and has the same sense as subsurface flux tubes giving rise to active regions of solar cycle 21.These observations are used to discusse large-scale electric currents crossing the photosphere and angular momentum loss to the solar wind.Now at Kitt Peak National Observatory, Tucson, Ariz. 85726, U.S.A.     

Cyclic Reversal of Magnetic Cloud Poloidal Field

Coronal mass ejections (CMEs) carry magnetic structure from the low corona into the heliosphere. The interplanetary CMEs (ICMEs) that exhibit the topology of helical magnetic fluxropes are traditionally called magnetic clouds (MCs). MC fluxropes with axis of low (high) inclination with respect to the ecliptic plane have been referred to as bipolar (unipolar) MCs. The poloidal field of bipolar MCs has a solar cycle dependence. We report a cyclic reversal of the poloidal field of low inclination MC fluxropes during 1976 to 2009. The MC poloidal field cyclic reversal on the same time scale of the solar magnetic cycle is evident over three sunspot cycles. Approximately 48% of ICMEs are MCs, and 40% of IMCs are bipolar MCs during solar cycle 23. The speed of the bipolar MCs has essentially the same distribution as all ICMEs, which implies that they are not from any special type of CMEs in terms of the solar origin. Although CME fluxropes may undergo a number of complications during the eruption and propagation, a significant group of MCs retains sufficient similarity to the source region magnetic field to posses the same cyclic periodicity in polarity reversal. The poloidal field of bipolar MCs gives the out-of-ecliptic-plane field or B _z component in the IMF time series. MCs with southward B _z field are particularly effective in causing geomagnetic disturbances. During the solar minima, the B _z field IMF sequence within MCs at the leading portion of a bipolar MC is the same with the solar global dipole field. Our finding shows that MCs preferentially remove the like polarity of the solar dipole field, and it supports the participation of CMEs in the solar magnetic cycle.     

On the MHD Acceleration of Astrophysical Jets

We present a 2.5D magnetohydrodynamic (MHD) simulation of the acceleration of a collimated jet from a magnetized accretion disk. We employ a MHD Adaptive Mesh Refinement (AMR) code (FLASH—University of Chicago). Thanks to this tool we can follow the evolution of the system for many dynamical timescales with a high-spatial resolution. Assuming an initial condition in which a Keplerian disk, thus with no accretion motions, is threaded by a uniform poloidal magnetic field, we show how both the accretion flow and the acceleration of the outflow occur, and we present in detail which are the forces responsible for the jet launching and collimation. Our simulation also shows how the collimating forces due to the self-generated toroidal magnetic field can produce some peculiar knotty features.     

The solar wind and the temperature-density structure of the solar corona

The influence of the solar wind on large-scale temperature and density distributions in the lower corona is studied. This influence is most profoundly felt through its effect upon the geometry of coronal magnetic fields since the presence of expansion divides the corona into magnetically open and closed regions. Each of these regions is governed by entirely different energy transport processes. This results in significant temperature differences since only the open field regions suffer outward conductive heat losses. Because the temperature influences the density in an exponential manner, large density inhomogeneities are to be expected.An approximate method for calculating the temperature and density distribution in a known magnetic field geometry is outlined and numerical estimates are carried out for representative coronal conditions. These estimates show that temperature differences of a factor of about two and density differences of ten can be expected in the lower corona even for uniform base conditions. As a result, we do not regard the so-called coronal holes necessairly as locations of reduced mechanical heating. Alternatively, we suggest that they are regions of open magnetic field lines being continuously drained of energy contert by the solar wind expansion and outward thermal conduction.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Phase relation between the poloidal and toroidal solar-cycle general magnetic fields and location of the origin of the surface magnetic fields

The phase relation of the poloidal and toroidal components of the solar-cycle general magnetic fields, which propagate along isorotation surfaces as dynamo waves, is investigated to infer the structure of the differential rotation and the direction of the regeneration action of the dynamo processes responsible for the solar cycle. It is shown that, from the phase relation alone, (i) the sign of the radial gradient of the differential rotation (/r) can be determined in the case that the radial gradient dominates the differential rotation, and (ii) the direction of the regeneration action can be determined in the case that the latitudinal gradient (/) dominates the differential rotation. Examining the observed poloidal and toroidal fields, it is concluded that (i) the / should dominate the differential rotation, and (ii) the determined sign of the regeneration factor (positive [negative] in the northern [southern] hemisphere) describing the direction of the regeneration action requires that the surface magnetic fields should originate from the upper part of the convection zone according to the model of the solar cycle driven by the dynamo action of the global convection.     

The geometrical height-scale and the pressure equilibrium in the sunspot umbra

Spectra of spots very near to the solar limb (limb distance 8) are used to determine the height difference between the levels of formation of the continuum and the line cores of 60 medium-strong Fraunhofer lines. For all lines (with Rowland Intensity < 10), this difference is < 1 (= 725 km) and well correlated with the Rowland intensity. The line absorption coefficient is calculated for some lines with known oscillator strength. This gives a possibility to deduce a value for the scale height of the umbra, which is found to be about 100 km, thus being equal to the photospheric scale height. Pure hydrostatic equilibrium exists, therefore, in the umbra, and vertical magnetic forces are negligible. Other methods for determining the scale height are discussed for comparison.The horizontal pressure equilibrium is discussed by taking into account the Wilson effect, and by neglecting dynamic terms (flow of matter). The magnetic field is confirmed to be force-free in higher layers (chromosphere). The pressure difference umbra-photosphere increases towards deeper layers, having a maximum at ^* - 1 which corresponds to about two times the magnetic pressure H ²/8. If rotational symmetry of the field is assumed, this can be explained by a minimum radius of curvature of the field lines of 1/4 spot radius.Mitteilungen aus dem Fraunhofer Institut Nr. 90.     

Stability of accretion discs threaded by a strong magnetic field

We study the stability of poloidal magnetic fields anchored in a thin accretion disc. The two-dimensional hydrodynamics in the disc plane is followed by a grid-based numerical simulation including the vertically integrated magnetic forces. The three-dimensional magnetic field outside the disc is calculated in a potential field approximation from the magnetic flux density distribution in the disc. For uniformly rotating discs we confirm numerically the existence of the interchange instability as predicted by Spruit, Stehle & Papaloizou . In agreement with predictions from the shearing sheet model, discs with Keplerian rotation are found to be stabilized by the shear, as long as the contribution of magnetic forces to support against gravity is small. When this support becomes significant, we find a global instability which transports angular momentum outwardly and allows mass to accrete inwardly. The instability takes the form of a m =1 rotating 'crescent', reminiscent of the purely hydrodynamic non-linear instability previously found in pressure-supported discs. A model where the initial surface mass density Σ( r ) and B _z ( r ) decrease with radius as power laws shows transient mass accretion during about six orbital periods, and settles into a state with surface density and field strength decreasing approximately exponentially with radius. We argue that this instability is likely to be the main angular momentum transport mechanism in discs with a poloidal magnetic field sufficiently strong to suppress magnetic turbulence. It may be especially relevant in jet-producing discs.     

A simulation study of the solar wind including the solar rotation effect

An axisymmetric solar wind structure including the solar rotation effect is studied by the method of MHD computer simulation. For the case of the radial magnetic field configuration, the simulation result is fairly well coincident with the steady-state solution. For the case of the dipole magnetic field configuration, the properties of the solution depend on the ratio of the gas pressure to the magnetic pressure-ratio) in the model. If the-ratio is small, a clearly defined stagnation region appears in the wind, in which the flow speed is very small and the azimuthal magnetic field is very weak because of the corotation of the plasma. If the-ratio is greater than 1, the plasma is not effectively trapped by the magnetic field so that the stagnation region is not clearly defined in the solution.     

Optimal growth of small perturbations in thin gaseous disks

A thin gaseous disk with a nearly Keplerian rotation profile and free boundaries in the external gravitational field of a point gravitating object does not generate any growing perturbation eigenmodes. In spite of this, a significant transient growth of linear perturbations measured by the evolution of their total acoustic energy is possible in such a disk. This is shown within the framework of the simplest model of an inviscid polytropic thin disk with a finite radial extent in which small adiabatic perturbations that are a linear combination of neutral eigenmodes with a corotation radius beyond the outer flow boundary are considered.     

Solar wind speed and its acceleration inferred using the interplanetary scintillation method in carrington rotation 1753

Solar wind speeds (SWSs) estimated by interplanetary scintillation (IPS) observations during Carrington rotation 1753 are projected onto the so-called source-surface of 2.5 solar radii along magnetic field lines in interplanetary space. The following two working hypotheses are examined from different points of view: (1) The SWS is a weighted mean along the line of sight to a radio source; the weight for the SWS depends on the distance from theP-point, the closest approach to the Sun on the line of sight. (2) The weighting function has a very sharp peak at theP-point, so that the SWS shows a real solar wind speed at theP-point. In both the two cases, the SWSs projected onto the source surface are further projected onto the photosphere along magnetic field lines in the corona. Footpoints of these field lines are inferred as photospheric source regions of the solar wind. The intensity of the Hei (1083 nm) absorption line (HEI) in the chromosphere corresponding to these photospheric sources is interpolated from observational data. The weighted mean of the HEI is calculated in case 1. The HEI corresponding to theP-point is used in case 2. The SWS is compared with the HEI in both the two cases. It is found that the correlation between the SWS and the HEI is better in case 2 than in case 1. It is further inferred by correlation analysis between the SWS and the HEI that the solar wind is accelerated within 27 solar radii on average. Although the data examined in this paper were limited to just one solar rotation, these results suggest that the SWS estimated by the IPS technique corresponds to the solar wind speed near theP-point and the weighting function along the line of sight may have a very sharp peak near theP-point.     

A moving Type IV radio burst and its relation to the coronal magnetic field

A moving Type IV burst, observed with the Culgoora radioheliograph on 1970 April 29, moved out to about 3 R and attained high circular polarization before fading. The appearance of the moving Type IV source suggests an isolated, self-contained, synchrotron emitting plasmoid. Magnetic field maps of the corona derived from photospheric observations indicate that the plasmoid moved almost radially outward from the flare region along open field lines. To explain the observed source structure and high unipolar polarization, we suggest that a ring of electric current was ejected from the low corona and guided by coronal magnetic field lines; the radio emission was synchrotron radiation generated by mildly-relativistic electrons trapped in the poloidal magnetic field of the ring current.Part of the research reported here was carried out while the author was at the Division of Radiophysics, C.S.I.R.O., Sydney, Australia.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Hydromagnetic oscillations of a self-gravitating fluid spheroid

The oscillations of a homogeneous, compressible, self gravitating fluid spheroid in static equilibrium with a poloidal magnetic field inside and a dipole type field outside are studied using the second order tensor virial equations. It is found that for small values of the eccentricity, the equilibrium model is dynamically stable provided the usual criterion of pulsative stability in the absence of a magnetic field (>4/3) is satisfied. The magnetic field removes the accidental degeneracy of the radial and the non-radial modes of oscillation which exists for =1.6 in the absence of a magnetic field.     

Two-dimensional mapping of the Sun with the RATAN-600

We present two-dimensional solar maps at 2.7, 3.2, 4, and 8.2 cm computed from one-dimensional observations with the RATAN-600, using Earth rotation aperture synthesis techniques. Before the calculation of maps, the position of each scan was corrected with respect to the center of the solar disk and the scans were calibrated. The circular polarization scans were corrected for polarization cross-talk between the I and V channels. Subsequently, the quiet-Sun background emission was subtracted. After all corrections, a dirty map was computed by combining the scans at different position angles. The last step of the processing was an attempt to free the dirty map of the sidelobes, using the standard CLEAN procedure. The resolution of the clean maps at 2.7 cm was 0.5 by 6. Both active regions which were present on the solar disk were mapped. We studied the flux spectra of different types of sources: one was associated with a sunspot, the second was located over the neutral line of an active region, and the other was associated with the plage. The emission mechanism of the former was attributed to the gyroresonance process, while the short wavelength emission of the others was attributed to the free-free process. For the sunspot-associated source we estimated the magnetic field strengths at the base of the transition region and in the low corona.     

The circumstellar gas of sigma orionis E

In order to explain the variable H emission and the eclipse-like light variation of Ori E, we investigated the circumstellar gas trapped by the stellar magnetic field and corotating with the star. By considering the potential along the magnetic field line, we found that the gas concentrates to a potential minimum. The circumstellar gas forms either two condensations or a disk, depending on the inclination of the magnetic dipole to the stellar rotation axis. The geometrical thickness of the circumstellar disk, of about 0.2 stellar radii, and the distance from the center of the star to the inner edge of the disk, of about 3 stellar radii, were obtained. The H emission line profile at its maximum phase and the amplitude of light variation were calculated by assuming the isothermal gas in LTE with the maximum gas density which the magnetic field can hold. The model gives good agreement with observation in the low obliquity case, and also explains the phase correlation among the H emission maximum, the light minimum, and the magnetic extreme. The model, however, failed to explain the large IR excess in theM band.     

A Rossby-wave dynamo for the Sun,II

To make the analysis more tractable, we simplify the equations of Part I to apply to two superposed layers of fluid, with horizontal variations in the motion and magnetic field represented by a small number of Fourier harmonics. The resulting set of eighteen ordinary nonlinear differential equations in time for the Fourier amplitudes is integrated numerically. We analyze in detail the dynamo action from a typical Rossby wave motion and compare it with the solar cycle.The field reversal process is similar in some respects to that put forth by Babcock. Toroidal fields are dragged up by vertical motions in the Rossby waves to form large-scale vertical fields, whose polarities alternate with longitude roughly like bipolar magnetic regions. Vertical fields of preferentially one polarity are carried toward the pole by the meridional motion in the wave to form an axisymmetric poloidal field. This poloidal field is then stretched out by the differential rotation into a new toroidal field of the opposite sign from the original. The poloidal field changes sign when the toroidal and bipolar region like fields are maximum, and vice versa.For the case studied, the reversal period is too short ( 2 years) and the poloidal fields too large ( 40 G) for the sun. Improvements for the model are discussed.Part I has been published in Solar Phys. 8, 316.