Magnetic Energy of Force-Free Fields with Detached Field Lines

Using an axisymmetrical ideal MHD model in spherical coordinates, we present a numerical study of magnetic configurations characterized by a levitating flux rope embedded in a bipolar background field whose normal field at the solar surface is the same or very close to that of a central dipole. The characteristic plasma β (the ratio between gas pressure and magnetic pressure) is taken to be sosmall (β= 10^-4) that the magnetic field is close to being force-free. The system as a whole is then let evolve quasi-statically with a slow increase of either the annular magnetic flux or the axial magnetic flux of the rope, and the total magneticenergy of the system grows accordingly. It is found that there exists an energy threshold: the flux rope sticks to the solar surface in equilibrium if the magneticenergy of the system is below the threshold, whereas it loses equilibrium if the threshold is exceeded. The energy threshold is found to be larger than that of thecorresponding fully-open magnetic field by a factor of nearly 1.08 irrespective as towhether the background field is completely closed or partly open, or whether the magnetic energy is enhanced by an increase of annular or axial flux of the rope.This gives an example showing that a force-free magnetic field may have an energy larger than the corresponding open field energy if part of the field lines is allowed tobe detached from the solar surface. The implication of such a conclusion in coronal mass ejections is briefly discussed and some comments are made on the maximum energy of force-free magnetic fields.     

Energy Buildup,Flux Confinement and Helicity Accumulation in the Solar Corona

Starting from a dipole field and a given distribution of footpoint displace- ment of field lines on the photosphere,we find axisymmetric,force-free field solutions in spherical coordinates that have the same distribution of normal field on the photo- sphere and magnetic topology as the dipole field.A photospheric shear is introduced in the azimuthal direction in a region that strides across the equator and ends at latitude γ_s.The footpoint displacement has a sine distribution in latitude and a peak amplitude of (?)_m.The magnetic energy E,azimuthal flux F_(?),and magnetic helicity H_T in the solar corona are then calculated for each force-free field solution.It is found that for a given shear region range γ_s,all of the three quantities increase monotoni- cally with increasing (?)_m.In particular,both F_(?) and H_T have a linear dependence on (?)m.When (?)_m reaches a certain critical value (?)_(mc),the force-free field loses equilib- rium,leading to a partial opening of the field and the appearance of a current sheet in the equatorial plane.At this point,E,F_(?)and H_T reach their maximum values, E_c,F_((?)c) and H_(Tc).E_c increases,and F_((?)c) and H_(Tc)decrease with decreasing λ_s.It is found that E_c is always smaller than the open field energy,in agreement with the Aly conjecture.Of the three critical parameters,E_c has the weakest dependence on λ_s.Therefore,if one is interested in the transition of a magnetic configuration from a stable state to a dynamic one,the magnetic energy is probably the most appropriate marker of the transition.     

Evolution of filamentary molecular clouds in the presence of magnetic fields

The purpose of this paper is to explore the effect of magnetic fields on the dynamics of magnetized filamentary molecular clouds.We suppose there is a filament with cylindrical symmetry and two components of axial and toroidal magnetic fields.In comparison to previous works,the novelty in the present work involves a similarity solution that does not define a function of the magnetic fields or density.We consider the effect of the magnetic field on the collapse of the filament in both axial and toroidal directions and show that the magnetic field has a braking effect,which means that the increasing intensity of the magnetic field reduces the velocity of collapse.This is consistent with other studies.We find that the magnetic field in the central region tends to be aligned with the filament axis.Also,the magnitude and the direction of the magnetic field depend on the magnitude and direction of the initial magnetic field in the outer region.Moreover,we show that more energy dissipation from the filament causes a rise in the infall velocity.     

Research on the principle of space high-precision temperature control system of liquid crystals based Stokes polarimeter

The magnetic field is one of the most important parameters in solar physics,and a polarimeter is the key device to measure the solar magnetic field.Liquid crystals based Stokes polarimeter is a novel technology,and will be applied for magnetic field measurement in the first space-based solar observatory satellite developed by China,Advanced Space-based Solar Observatory.However,the liquid crystals based Stokes polarimeter in space is not a mature technology.Therefore,it is of great scientific significance to study the control method and characteristics of the device.The retardation produced by a liquid crystal variable retarder is sensitive to the temperature,and the retardation changes 0.09°per 0.10℃.The error in polarization measurement caused by this change is 0.016,which affects the accuracy of magnetic field measurement.In order to ensure the stability of its performance,this paper proposes a high-precision temperature control system for liquid crystals based Stokes polarimeter in space.In order to optimize the structure design and temperature control system,the temperature field of liquid crystals based Stokes polarimeter is analyzed by the finite element method,and the influence of light on the temperature field of the liquid crystal variable retarder is analyzed theoretically.By analyzing the principle of highprecision temperature measurement in space,a high-precision temperature measurement circuit based on integrated operational amplifier,programmable amplifier and 12 bit A/D is designed,and a high-precision space temperature control system is developed by applying the integral separation PI temperature control algorithm and PWM driving heating films.The experimental results show that the effect of temperature control is accurate and stable,whenever the liquid crystals based Stokes polarimeter is either in the air or vacuum.The temperature stability is within±0.0150℃,which demonstrates greatly improved stability for the liquid crystals based Stokes polarimeter.     

Polarization Position Angle Swings caused by Relativistic Effects

Polarization position angle swings of - 180 ° observed in extragalactic radio sources are a regular behavior of variability in polarization. They should be due to some kind of physically regular process. We consider relativistic shocks which propagate through and 'illuminate' regular configurations of magnetic field, producing polarization angle swing events. Two magnetic field configurations (force-free field and homogeneous helical field) are considered to demonstrate the results. It is shown that the properties of polarization angle swings and the relationship between the swings and variations in total and polarized flux density are critically dependent on the configuration of magnetic field and the dynamical behavior of the shock. In particular, we find that in some cases polarization angle swings can occur when the total and polarized flux densities only vary by a very small amount. These results may be useful for understanding the polarization variability with both long and short timescales obser     

Testing a solar coronal magnetic field extrapolation code with the Titov–Dmoulin magnetic flux rope model

In the solar corona, the magnetic flux rope is believed to be a fundamental structure that accounts for magnetic free energy storage and solar eruptions. Up to the present, the extrapolation of the magnetic field from boundary data has been the primary way to obtain fully three-dimensional magnetic information about the corona. As a result, the ability to reliably recover the coronal magnetic flux rope is important for coronal field extrapolation. In this paper, our coronal field extrapolation code is examined with an analytical magnetic flux rope model proposed by Titov D′emoulin, which consists of a bipolar magnetic configuration holding a semi-circular line-tied flux rope in force-free equilibrium. By only using the vector field at the bottom boundary as input, we test our code with the model in a representative range of parameter space and find that the model field can be reconstructed with high accuracy. In particular, the magnetic topological interfaces formed between the flux rope and the surrounding arcade, i.e., the "hyperbolic flux tube" and "bald patch separatrix surface," are also reliably reproduced. By this test, we demonstrate that our CESE–MHD–NLFFF code can be applied to recovering the magnetic flux rope in the solar corona as long as the vector magnetogram satisfies the force-free constraints.     

Effects of dynamo magnetic fields on observational properties of accreting millisecond X-ray pulsars

In this paper,we have investigated accreting millisecond X-ray pulsars,which are rapidly rotating neutron stars in low-mass X-ray binaries.These systems exhibit coherent X-ray pulsations that arise when the accretion flow is magnetically channeled to the stellar surface.Here,we have developed the fundamental equations for an accretion disk around accreting millisecond X-ray pulsars in the presence of a dynamo generated magnetic field in the inner part of the disk.We have also formulated the numerical method for the structure equations in the inner region of the disk and the highest accretion rate is enough to form the inner region of the disk,which is overpowered by radiation pressure and electron scattering.Finally,we have examined our results with the effects of dynamo magnetic fields on accreting millisecond X-ray pulsars.     

A disk-corona model for the low/hard state of black hole X-ray binaries

A disk-corona model for fitting the low/hard（LH）state of the associated steady jet in black hole X-ray binaries（BHXBs）is proposed based on the large-scale magnetic field configuration that arises from the coexistence of the Blandford-Znajek（BZ）and Blandford-Payne（BP）processes,where the magnetic field configuration for the BP process is determined by the requirement of energy conversion from Poynting energy flux into kinetic energy flux in the jet.It is found that corona current is crucial to guarantee the consistency of the jet launching from the accretion disk.The relative importance of the BZ and BP processes in powering jets from black hole accretion disks is discussed,and the LH state of several BHXBs is fitted based on our model.In addition,we suggest that magnetic field configuration can be regarded as the second parameter for governing the state transition of BHXBs.     

A New Method of Identifying 3D Null Points in Solar Vector Magnetic Fields

Employing the Poincare index of isolated null-points in a vector field, we worked out a mathematical method of searching for 3D null-points in coronal magnetic fields. After introducing the relevant differential topology, we test the method by using the analytical model of Brown & Priest. The location of null-point identified by our method coincides precisely with the analytical solution. Finally we apply the method to the 3D coronal magnetic fields reconstructed from an observed MDI magnetogram of a super-active region (NOAA 10488). We find that the 3D null-point seems to be a key element in the magnetic topology associated with flare occurrence.     

Vector Magnetic Field Measurement of NOAA AR 10197

A set of two-dimensional Stokes spectral data of NOAA AR 10197 obtained by the Solar Stokes Spectral Telescope (S3T) at the Yunnan Observatory are qualitatively analyzed. The three components of the vector magnetic field, the strength H, inclinationγand azimuth X, are derived. Based on the three components, we contour the distributions of the longitudinal magnetic field and transverse magnetic field. The active region (AR) has two different magnetic polarities apparent in the longitudinal magnetic map due to projection effect. There is a basic agreement on the longitudinal magnetic fields between the S3T and SOHO/MDI magnetograms, with a correlation coefficientρBl = 0.911. The transverse magnetic field of the AR has a radial distribution from a center located in the southwest of the AR. It is also found that the transverse magnetic fields obtained by Huairou Solar Observing Station (HRSOS) have a similar radial distribution. The distributions of transverse magnetic field obtained by S3T and HRSOS have correlation coefficients,ρAzimu = 0.86 andρBt = 0.883, in regard to the azimuthal angle and intensity.     

Nonlinear Force-Free Magnetic Field Extrapolations: Comparison of the Grad Rubin and Wheatland Sturrock Roumeliotis Algorithm

We compare the performance of two alternative algorithms which aim to construct a force-free magnetic field given suitable boundary conditions. For this comparison, we have implemented both algorithms on the same finite element grid which uses Whitney forms to describe the fields within the grid cells. The additional use of conjugate gradient and multigrid iterations result in quite effective codes. The Grad Rubin and Wheatland Sturrock Roumeliotis algorithms both perform well for the reconstruction of a known analytic force-free field. For more arbitrary boundary conditions the Wheatland Sturrock Roumeliotis approach has some difficulties because it requires overdetermined boundary information which may include inconsistencies. The Grad Rubin code on the other hand loses convergence for strong current densities. For the example we have investigated, however, the maximum possible current density seems to be not far from the limit beyond which a force-free field cannot exist anymore for a given normal magnetic field intensity on the boundary.     

Calculating and Testing Nonlinear Force-Free Fields

Improvements to an existing method for calculating nonlinear force-free magnetic fields (Wheatland, Solar Phys. 238, 29, 2006) are described. In particular a solution of the 3-D Poisson equation using 2-D Fourier transforms is presented. The improved nonlinear force-free method is demonstrated in application to linear force-free test cases with localized nonzero values of the normal component of the field in the boundary. These fields provide suitable test cases for nonlinear force-free calculations because the boundary conditions involve localized nonzero values of the normal components of the field and of the current density, and because (being linear force-free fields) they have more direct numerical solutions. Despite their simplicity, fields of this kind have not been recognized as test cases for nonlinear methods before. The examples illustrate the treatment of the boundary conditions on current in the nonlinear force-free method, and in particular the limitations imposed by field lines that connect outside of the boundary region.     

A class of force-free magnetic fields for modeling pre-flare coronal magnetic configurations

Three-dimensional, quasi-static evolutions of coronal magnetic fields driven by photospheric flux emergence are modeled by a class of analytic force-free magnetic fields. Our models relate commonly observed photospheric magnetic phenomena, such as the formation and growth of sunspots, the emergence of an X-type separator, and the collision and merging of sunspots, to the three-dimensional magnetic fields in the corona above. By tracking the evolution in terms of a continuous sequence of force-free states, we show that flux emergence and submergence along magnetic neutral lines in the photosphere are essential processes in all these photospheric phenomena. The analytic solutions we present have a parametric regime within which the magnetic energy attained by an evolving force-free field may be of the order of 10³⁰ ergs to several 10³¹ ergs, depending on the magnetic environment into which an emerging flux intrudes. The commonly used indicators of magnetic shear in magnetogram interpretation are discussed in terms of field connectivity in our models. It is demonstrated that the crossing angle of the photospheric transverse magnetic field with the neutral line may not be a reliable indicator of the magnetic shear in the coronal field above, due to the complexity of three-dimensionality. The poorly understood constraint of magnetic-helicity conservation on the availability of magnetic free energy for a flare is briefly discussed.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Force-free magnetic fields in the solar atmosphere

Reliable measurements of the solar magnetic field are restricted to the level of the photosphere. For about half a century attempts have been made to calculate the field in the layers above the photosphere, i.e. in the chromosphere and in the corona, from the measured photospheric field. The procedure is known as magnetic field extrapolation. In the superphotospheric parts of active regions the magnetic field is approximately force-free, i.e. electric currents are aligned with the magnetic field. The practical application to solar active regions has been largely confined to constant-α or linear force-free fields, with a spatially constant ratio, α, between the electric current and the magnetic field. We review results obtained from extrapolations with constant-α force-free fields, in particular on magnetic topologies favourable for flares and on magnetic and current helicities. Presently, different methods are being developed to calculate non-constant-α or nonlinear force-free fields from photospheric vector magnetograms. We also briefly discuss these methods and present a comparison of a linear and a nonlinear force-free magnetic field extrapolation applied to the same photospheric boundary data. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

An Improved Approach to Non-Force-Free Coronal Magnetic Field Extrapolation

We develop an approach to deriving the three-dimensional non-force-free coronal magnetic field from vector magnetograms. Based on the principle of minimum dissipation rate, a general non-force-free magnetic field is expressed as the superposition of one potential field and two constant-α (linear) force-free fields. Each is extrapolated from its bottom boundary data, providing the normal component only. The constant-α parameters are distinct and determined by minimizing the deviations between the numerically computed and measured transverse magnetic field at the bottom boundary. The boundary conditions required are at least two layers of vector magnetograms, one at the photospheric level and the other at the chromospheric level, presumably. We apply our approach to a few analytic test cases, especially to two nonlinear force-free cases examined by Schrijver et al. (Solar Phys. 235, 161, 2006). We find that for one case with small α parameters, the quantitative measures of the quality of our result are better than the median values of those from a set of nonlinear force-free methods. The reconstructed magnetic-field configuration is valid up to a vertical height of the transverse scale. For the other cases, the results remain valid to a lower vertical height owing to the limitations of the linear force-free-field solver. Because our method is based on the fast-Fourier-transform algorithm, it is much faster and easy to implement. We discuss the potential usefulness of our method and its limitations.     

Force-free magnetic arcades relevant to two-ribbon solar flares

Simple analytic models for the passive evolution of arcade-like magnetic fields through a series of force-free equilibria are presented. At the photospheric boundary, the normal magnetic field component is prescribed together with either the longitudinal field component or the photospheric shear. Analytic progress is made by considering either cylindrically symmetric solutions or using the separation of variables technique. Two distinct cylindrically symmetric force-free fields are obtained that possess the same normal field component and photospheric shear. The scond field contains a magnetic bubble. As the shear increases beyond a critical value, so the magnetic energy of the first configuration exceeds that of the second. The possibility is therefore suggested of an eruption of the first field outwards towards the second. Such an eruptive instability is proposed as the origin of a two-ribbon solar flare.A new analytic solution to the force-free field equations, of separable form, is discovered and it is pointed out that the existence of shear in a magnetic field does not preclude it from being potential.Now at AWRE, Aldermaston, Reading, Berkshire.     

Magnetic field configurations associated with polarity intrusion in a solar active region

This paper presents a new class of exact solutions describing the non-linear force-free field above a spatially localized photospheric bipolar magnetic region. An essential feature is the variation in all three Cartesian directions and this could not be modelled adequately with previously known symmetric force-free fields. Sequences of force-free fields are constructed and analyzed to simulate the slow growth of a pair of spots on the photosphere. The axis connecting the spots executes rotational motion, distorting the photospheric neutral line separating fluxes of opposite signs. We show directly from the analytic solutions that the resulting reversal of the positions of the spots relative to the background field is associated with (i) the creation of magnetic free energy, (ii) the severe shearing of localized low-lying loops in the vicinity where the photospheric transverse field aligns with the photospheric neutral line, and (iii) the emergence and disappearance of flux from the photosphere at these highly stressed regions. The model relates theoretically for the first time these different magnetic field features that have been suggested by observation and theoretical considerations to be flare precursors. A general formula, based on the virial theorem, is also given for the free energy of a force-free field, strictly in terms of the field value at the photosphere. This formula has obvious practical application.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Non-ideal evolution of non-axisymmetric, force-free magnetic fields in a magnetar

Recent numerical magnetohydrodynamic calculations by Braithwaite and collaborators support the 'fossil field' hypothesis regarding the origin of magnetic fields in compact stars and suggest that the resistive evolution of the fossil field can explain the reorganization and decay of magnetar magnetic fields. Here, these findings are modelled analytically by allowing the stellar magnetic field to relax through a quasi-static sequence of non-axisymmetric, force-free states, by analogy with spheromak relaxation experiments, starting from a random field. Under the hypothesis that the force-free modes approach energy equipartition in the absence of resistivity, the output of the numerical calculations is semiquantitatively recovered: the field settles down to a linked poloidal–toroidal configuration, which inflates and becomes more toroidal as time passes. A qualitatively similar (but not identical) end state is reached if the magnetic field evolves by exchanging helicity between small and large scales according to an α-dynamo-like, mean-field mechanism, arising from the fluctuating electromotive force produced by the initial random field. The impossibility of matching a force-free internal field to a potential exterior field is discussed in the magnetar context.     

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.