An unconventional accretion disc dynamo

In the light of recent results from numerical simulations of accretion disc MHD turbulence, we revisit the problem of the configuration of large-scale magnetic fields resulting from an α Ω dynamo operating in a thin accretion disc. In particular, we analyse the consequences of the peculiar sign of the α -effect suggested by numerical simulations . We determine the symmetry of the fastest-growing modes in the kinematic dynamo approximation and, in the framework of an ' α -quenched' dynamo model, study the evolution of the magnetic field. We find that the resulting field for this negative polarity of the α -effect generally has dipole symmetry with respect to the disc midplane, although the existence of an equilibrium configuration depends on the properties of the turbulence. The role of magnetic field dragging is discussed and, finally, the presence of an external uniform magnetic field is included to address the issue of magneto centrifugal wind launching from accretion discs.     

Mean-field effects in the Galloway–Proctor flow

In the framework of mean-field electrodynamics the coefficients defining the mean electromotive force in Galloway–Proctor flows are determined. These flows show a two-dimensional pattern and are helical. The pattern wobbles in its plane. Apart from one exception a circularly polarized Galloway–Proctor flow, i.e. a circular motion of the flow pattern is assumed. This corresponds to one of the cases considered recently by Courvoisier, Hughes & Tobias. An analytic theory of the α effect and related effects in this flow is developed within the second-order correlation approximation and a corresponding fourth-order approximation. In the validity range of these approximations there is an α effect but no γ effect, or pumping effect. Numerical results obtained with the test-field method, which are independent of these approximations, confirm the results for α and show that γ is in general non-zero. Both α and γ show a complex dependency on the magnetic Reynolds number and other parameters that define the flow, that is, amplitude and frequency of the circular motion. Some results for the magnetic diffusivity  η_t  and a related quantity are given, too. Finally, a result for α in the case of a randomly varying linearly polarized Galloway–Proctor flow, without the aforementioned circular motion, is presented. The flows investigated show quite interesting effects. There is, however, no straightforward way to relate these flows to turbulence and to use them for studying properties of the α effect and associated effects under realistic conditions.     

The role of the Yoshizawa effect in the Archontis dynamo

The generation of mean magnetic fields is studied for a simple non-helical flow where a net cross-helicity of either sign can emerge. This flow, which is also known as the Archontis flow, is a generalization of the Arnold–Beltrami–Childress flow, but with the cosine terms omitted. The presence of cross-helicity leads to a mean-field dynamo effect that is known as the Yoshizawa effect. Direct numerical simulations of such flows demonstrate the presence of magnetic fields on scales larger than the scale of the flow. Contrary to earlier expectations, the Yoshizawa effect is found to be proportional to the mean magnetic field and can therefore lead to its exponential instead of just linear amplification for magnetic Reynolds numbers that exceed a certain critical value. Unlike α effect dynamos, it is found that the Yoshizawa effect is not notably constrained by the presence of a conservation law. It is argued that this is due to the presence of a forcing term in the momentum equation, which leads to a non-zero correlation with the magnetic field. Finally, the application to energy convergence in solar wind turbulence is discussed.     

Mean-field dynamo in partially ionized plasmas – I

There are several astrophysical situations where one needs to study the dynamics of magnetic flux in partially ionized turbulent plasmas. In a partially ionized plasma, the magnetic induction is subjected to the ambipolar diffusion and the Hall effect in addition to the usual resistive dissipation. In this paper, we initiate the study of the kinematic dynamo in a partially ionized turbulent plasma. The Hall effect arises from the treatment of the electrons and the ions as two separate fluids and the ambipolar diffusion due to the inclusion of neutrals as the third fluid. It is shown that these non-ideal effects modify the so-called α effect and the turbulent diffusion coefficient β in a rather substantial way. The Hall effect may enhance or quench the dynamo action altogether. The ambipolar diffusion brings in an α which depends on the mean magnetic field. The new correlations embodying the coupling of the charged fluids and the neutral fluid appear in a decisive manner. The turbulence is necessarily magnetohydrodynamic with new spatial and time-scales. The nature of the new correlations is demonstrated by taking the Alfvénic turbulence as an example.     

The helicity constraint in turbulent dynamos with shear

The evolution of magnetic fields is studied using simulations of forced helical turbulence with strong imposed shear. After some initial exponential growth, the magnetic field develops a large-scale travelling wave pattern. The resulting field structure possesses magnetic helicity, which is conserved in a periodic box by the ideal magnetohydrodynamics equations and can hence only change on a resistive time-scale. This strongly constrains the growth time of the large-scale magnetic field, but less strongly constrains the length of the cycle period. Comparing this with the case without shear, the time-scale for large-scale field amplification is shortened by a factor Q , which depends on the relative importance of shear and helical turbulence, and which also controls the ratio of toroidal to poloidal field. The results of the simulations can be reproduced qualitatively and quantitatively with a mean-field α Ω-dynamo model with alpha-effect and turbulent magnetic diffusivity coefficients that are less strongly quenched than in the corresponding α ²-dynamo.     

Simple tests for the ideal MHD stability of line-tied coronal magnetic fields

Methods for investigating the stability of line-tied, cylindrically-symmetric magnetic fields are presented. The energy method is used and the perturbed potential energy integral is manipulated to produce simple tests that predict either stability to general coronal disturbances or instability to localized modes, both satisfying photospheric line-tying. Using these tests the maximum amount of magnetic energy, that can be stored in the coronal magnetic field prior to an instability, can be estimated. The tests are applied to four different classes of equilibria and results are obtained for both arcade and loop geometries.     

Viscous-resistive ADAF with a general large-scale magnetic field

We have studied the structure of hot accretion flow bathed in a general large-scale magnetic field. We have considered magnetic parameters , where are the Alfvén sound speeds in three direction of cylindrical coordinate (r,φ,z). The dominant mechanism of energy dissipation is assumed to be the magnetic diffusivity due to turbulence and viscosity in the accretion flow. Also, we adopt a more realistic model for kinematic viscosity (ν=αc _s H), with both c _s and H as a function of magnetic field. As a result in our model, the kinematic viscosity and magnetic diffusivity (η=η ₀ c _s H) are not constant. In order to solve the integrated equations that govern the behavior of the accretion flow, a self-similar method is used. It is found that the existence of magnetic resistivity will increase the radial infall velocity as well as sound speed and vertical thickness of the disk. However the rotational velocity of the disk decreases by the increase of magnetic resistivity. Moreover, we study the effect of three components of global magnetic field on the structure of the disk. We found out that the radial velocity and sound speed are Sub-Keplerian for all values of magnetic field parameters, but the rotational velocity can be Super-Keplerian by the increase of toroidal magnetic field. Also, Our numerical results show that all components of magnetic field can be important and have a considerable effect on velocities and vertical thickness of the disk.     

Computation of inner coronal magnetic fields from longitudinal field components on a spherical photosphere

A method is developed which for a certain day permits the approximate calculation of closed small and large scale magnetic field lines. From the photospheric longitudinal components of the magnetic field measured at this day normal components are derived taking into account the curvature of the solar surface. The magnetic fields are assumed to be potential or force-free fields.The method is applied to observations of September 5 and September 7, 1973. The projected magnetic field lines are compared with the loop structures which are visible in XUV pictures taken on these days. In the cases where no good agreement could be obtained for potential fields, force-free fields are calculated and fitted to the observed structures.     

Turbulent dynamo in asymptotic giant branch stars

Using recent results on the operation of turbulent dynamos, we show that a turbulent dynamo may amplify a large-scale magnetic field in the envelopes of asymptotic giant branch (AGB) stars. We propose that a slow rotation of the AGB envelope can fix the symmetry axis, leading to the formation of an axisymmetric magnetic field structure. Unlike solar-type αω dynamos, the rotation has only a small role in amplifying the toroidal component of the magnetic field; instead of an αω dynamo we propose an α ² ω . The magnetic field may reach a value of     , where B _e is the equipartition (between the turbulent and magnetic energy densities) magnetic field. The large-scale magnetic field is strong enough for the formation of magnetic cool spots on the AGB stellar surface. The spots may regulate dust formation, and hence the mass-loss rate, leading to axisymmetric mass loss and the formation of elliptical planetary nebulae (PNe). Despite its role in forming cool spots, the large-scale magnetic field is too weak to play a dynamic role and directly influence the wind from the AGB star, as required by some models. We discuss other possible problems in models where the magnetic field plays a dynamic role in shaping the AGB winds, and argue that they cannot explain the formation of non-spherical PNe.     

Resolving Azimuth Ambiguity Using Vertical Nature of Solar Quiet-Sun Magnetic Fields

The measurement of solar magnetic fields using the Zeeman effect diagnostics has a fundamental 180° ambiguity in the determination of the azimuth angle of the transverse field component. There are several methods that are used in the community and each one has its merits and demerits. Here we present a disambiguation idea that is based on the assumption that most of the magnetic field on the sun is predominantly vertical. While the method is not applicable to penumbra or other features harboring predominantly horizontal fields like the sheared neutral lines, it is useful for regions where fields are predominantly vertical like network and plage areas. The method is tested with the full-disk solar vector magnetograms observed by the SOLIS/VSM instrument. We find that statistically about 60?–?85 % of the pixels in a typical full-disk magnetogram has a field inclination in the range of 0?–?30° with respect to the local solar normal, and thus can be successfully disambiguated by the proposed method. Due to its non-iterative nature, the present method is extremely fast and therefore can be used as a good initial guess for iterative schemes like the non-potential field computation (NPFC) method. Furthermore, the method is insensitive to noisy pixels as it does not depend upon the neighboring pixels or derivatives.     

Magnetic helicity in stellar dynamos: new numerical experiments

The theory of large scale dynamos is reviewed with particular emphasis on the magnetic helicity constraint in the presence of closed and open boundaries. In the presence of closed or periodic boundaries, helical dynamos respond to the helicity constraint by developing small scale separation in the kinematic regime, and by showing long time scales in the nonlinear regime where the scale separation has grown to the maximum possible value. A resistively limited evolution towards saturation is also found at intermediate scales before the largest scale of the system is reached. Larger aspect ratios can give rise to different structures of the mean field which are obtained at early times, but the final saturation field strength is still decreasing with decreasing resistivity. In the presence of shear, cyclic magnetic fields are found whose period is increasing with decreasing resistivity, but the saturation energy of the mean field is in strong super‐equipartition with the turbulent energy. It is shown that artificially induced losses of small scale field of opposite sign of magnetic helicity as the large scale field can, at least in principle, accelerate the production of large scale (poloidal) field. Based on mean field models with an outer potential field boundary condition in spherical geometry, we verify that the sign of the magnetic helicity flux from the large scale field agrees with the sign of α. For solar parameters, typical magnetic helicity fluxes lie around 10⁴⁷ Mx² per cycle.     

Magnetic field dragging and the vertical structure of thin accretion discs

The presence of an imposed vertical magnetic field may drastically influence the structure of thin accretion discs. If the field is sufficiently strong, the rotation law can depart from the Keplerian one. We consider the structure of a disc for a given eddy magnetic diffusivity but neglect details of the energy transport. The magnetic field is assumed to be in balance with the internal energy of the accretion flow. The thickness of the disc as well as the turbulent magnetic Prandtl number and the viscosity, α , are the key parameters of our model. The calculations show that the radial velocity can reach the sound speed for a magnetic disc if the thickness is comparable to that of a non-magnetic one. This leads to a strong amplification of the accretion rate for a given surface density. The inclination angle of the magnetic field lines can exceed the critical value 30° (required to launch cold jets) even for a relatively small magnetic Prandtl number of order unity. The toroidal magnetic fields induced at the disc surface are smaller than predicted in previous studies.     

Topology of Magnetic Field and Coronal Heating in Solar Active Regions

Force-free magnetic fields can be computed by making use of a new numerical technique, in which the fields are represented by a boundary integral equation based on a specific Green's function. Vector magnetic fields observed on the photospheric surface can be taken as the boundary conditions of this equation. In this numerical computation, the following two points are emphasized: (1) A new method for data reduction is proposed, for removing uncertainties in boundary data and determining the parameter in this Green's function, which is important for solving the boundary integral equation. In this method, the transverse components of the observed boundary field are calibrated with a linear force-free field model without changing their azimuth. (2) The computed 3-D fields satisfy the divergence-free and force-free conditions with high precision. The alignment of these field lines is mostly in agreement with structures in Hα and Yohkoh soft X-ray images. Since the boundary data are calibrated with a linear force-free field model, the computed 3-D magnetic field can be regarded as a quasi-linear force-free field approximation. The reconstruction of 3-D magnetic field in active region NOAA 7321 was taken as an example to quantitatively exhibit the capability of our new numerical technique.     

Planetary dynamos: effects of electrically conducting flows overlying turbulent regions of magnetic field generation

A fully three-dimensional, nonlinear, time-dependent, multi-layered spherical kinematic dynamo model is used to study the effect on the observable external magnetic field of flow in an electrically conducting layer above a spherical turbulent dynamo region in which the α effect generates the magnetic field. It is shown that the amplitude and structure of an observable planetary magnetic field are largely determined by the magnitude and structure of the flow in the overlying layer. It is also shown that a strong-field planetary dynamo can be readily produced by the effect of an electrically conducting flow layer at the top of a convective core. The overlying layer and the underlying convective region constitute a magnetically strongly coupled system. Such overlying layers might exist at the top of the Earth's core due to chemical or thermal causes, in the cores of other terrestrial planets for similar reasons, and in Saturn due to the differentiation of helium from hydrogen. An electrically conducting and differentially rotating layer could exist above the metallic hydrogen region in Jupiter and affect the jovian magnetic field similar to the overlying layers in other planets. Lateral temperature gradients resulting in thermal winds drive the flow in the overlying layers. All planetary magnetic fields could be maintained by similar turbulent convective dynamos in the field-generation regions of planets with the differences among observable magnetic fields due to different circulations in the overlying electrically conducting layers.     

A Non-Linear Force-Free Field Model for the Evolving Magnetic Structure of Solar Filaments

In this paper the effect of a small magnetic element approaching the main body of a solar filament is considered through non-linear force-free field modeling. The filament is represented by a series of magnetic dips. Once the dips are calculated, a simple hydrostatic atmosphere model is applied to determine which structures have sufficient column mass depth to be visible in Hα. Two orientations of the bipole are considered, either parallel or anti-parallel to the overlying arcade. The magnetic polarity that lies closest to the filament is then advected towards the filament. Initially for both the dominant and minority polarity advected elements, right/left bearing barbs are produced for dextral/sinsitral filaments. The production of barbs due to dominant polarity elements is a new feature. In later stages the filament breaks into two dipped sections and takes a highly irregular, non-symmetrical form with multiple pillars. The two sections are connected by field lines with double dips even though the twist of the field is less than one turn. Reconnection is not found to play a key role in the break up of the filament. The non-linear force-free fields produce very different results to extrapolated linear-force free fields. For the cases considered here the linear force-free field does not produce the break up of the filament nor the production of barbs as a result of dominant polarity elements.     

Topology of Magnetic Field and Coronal Heating in Solar Active Regions – II. The Role of Quasi-Separatrix Layers

The photospheric vector magnetic fields, H and soft X-ray images of AR 7321 were simultaneously observed with the Solar Flare Telescope at Mitaka and the Soft X-ray Telescope of Yohkoh on October 26, 1992, when there was no important activity in this region. Taking the observed photospheric vector magnetic fields as the boundary condition, 3D magnetic fields above the photosphere were computed with a new numerical technique. Then quasi-separatrix layers (QSLs), i.e., regions where 3D magnetic reconnection takes place, were determined in the computed 3D magnetic fields. Since Yohkoh data and Mitaka data were obtained in well-arranged time sequences during the day, the evolution of 3D fields, H features and soft X-ray features in this region can be studied in detail. Through a comparison among the 3D magnetic fields, H features and soft X-ray features, the following results have been obtained: (a) H plages are associated with the portions of QSLs in the chromosphere; (b) diffuse coronal features (DCFs) and bright coronal features (BCFs) are morphologically confined by the coronal linkage of the field lines related to the QSLs; (c) BCFs are associated with a part of the magnetic field lines related to the QSLs. These results suggest that as the likely places where energy release may occur by 3D magnetic reconnection, QSLs play an important role in the chromospheric and coronal heating in this active region.     

Large-Scale Magnetic Field and Sunspot Cycles

Hα magnetic synoptic charts of the Sun are processed for 1915–1999 and the spherical harmonics are calculated. It is shown that the polarity distribution of the magnetic field on Hα charts is similar to the polarity distribution of the Stanford magnetic field observations during 1975–1999. The index of activity of the large-scale magnetic field A(t), representing the sum of the intensities of dipole and octupole components, is introduced. It is shown that the cycle of the large-scale magnetic field of the Sun precedes on the average by 5.5 years the sunspot activity cycle, W(t). This means that the weak large-scale magnetic fields of the Sun do not result from decay and diffusion of strong fields from active regions as it is supposed in all modern theories of the solar cycle. On the basis of the new data the intensity of the current solar cycle 23 is predicted and some aspects of the theory of the solar cycle are discussed.     

Including stereoscopic information in the reconstruction of coronal magnetic fields

We present a method to include stereoscopic information about the three-dimensional structure of flux tubes into the reconstruction of the coronal magnetic field. Due to the low plasma beta in the corona we can assume a force-free magnetic field, with the current density parallel to the magnetic field lines. Here we use linear force-free fields for simplicity. The method uses the line-of-sight magnetic field on the photosphere as observational input. The value of is determined iteratively by comparing the reconstructed magnetic field with the observed structures. The final configuration is the optimal linear force-free solution constrained by both the photospheric magnetogram and the observed plasma structures. As an example we apply our method to SOHO MDI/EIT data of an active region. In the future it is planned to apply the method to analyse data from the SECCHI instrument aboard the STEREO mission.