The Effects of Rotation on the Evolution of Rising Omega-loops in a Stratified Model Convection Zone

We present three-dimensional MHD simulations of buoyant magnetic flux tubes that rise through a stratified model convection zone in the presence of solar rotation. The equations of MHD are solved in the anelastic approximation, and the results are used to determine the effects of solar rotation on the dynamic evolution an Omega-loop. We find that the Coriolis force significantly suppresses the degree of fragmentation at the apex of the loop during its ascent toward the photosphere. If the initial axial field strength of the tube is reduced, then, in the absence of forces due to convective motions, the degree of apex fragmentation is also reduced. We show that the Coriolis force slows the rise of the tube, and induces a retrograde flow in both the magnetized and unmagnetized plasma of an emerging active region. Observationally, we predict that this flow will appear to originate at the leading polarity, and will terminate at the trailing polarity.Comment: 25 pages, 8 figures, ApJ in pres

Estimating Electric Fields from Vector Magnetogram Sequences

Determining the electric field (E-field) distribution on the Sun's photosphere is essential for quantitative studies of how energy flows from the Sun's photosphere, through the corona, and into the heliosphere. This E-field also provides valuable input for data-driven models of the solar atmosphere and the Sun-Earth system. We show how Faraday's Law can be used with observed vector magnetogram time series to estimate the photospheric E-field, an ill-posed inversion problem. Our method uses a "poloidal-toroidal decomposition" (PTD) of the time derivative of the vector magnetic field. The PTD solutions are not unique; the gradient of a scalar potential can be added to the PTD E-field without affecting consistency with Faraday's Law. We present an iterative technique to determine a potential function consistent with ideal MHD evolution; but this E-field is also not a unique solution to Faraday's Law. Finally, we explore a variational approach that minimizes an energy functional to determine a unique E-field, similar to Longcope's "Minimum Energy Fit". The PTD technique, the iterative technique, and the variational technique are used to estimate E-fields from a pair of synthetic vector magnetograms taken from an MHD simulation; and these E-fields are compared with the simulation's known electric fields. These three techniques are then applied to a pair of vector magnetograms of solar active region NOAA AR8210, to demonstrate the methods with real data.Comment: 41 pages, 10 figure

Flux-loss of buoyant ropes interacting with convective flows

We present 3-d numerical magneto-hydrodynamic simulations of a buoyant, twisted magnetic flux rope embedded in a stratified, solar-like model convection zone. The flux rope is given an initial twist such that it neither kinks nor fragments during its ascent. Moreover, its magnetic energy content with respect to convection is chosen so that the flux rope retains its basic geometry while being deflected from a purely vertical ascent by convective flows. The simulations show that magnetic flux is advected away from the core of the flux rope as it interacts with the convection. The results thus support the idea that the amount of toroidal flux stored at or near the bottom of the solar convection zone may currently be underestimated.Comment: 5 pages, 3 figures. Accepted for publication in Astronomy & Astrophysic

The Three-dimensional Evolution of Rising, Twisted Magnetic Flux Tubes in a Gravitationally Stratified Model Convection Zone

We present three-dimensional numerical simulations of the rise and fragmentation of twisted, initially horizontal magnetic flux tubes which evolve into emerging Omega-loops. The flux tubes rise buoyantly through an adiabatically stratified plasma that represents the solar convection zone. The MHD equations are solved in the anelastic approximation, and the results are compared with studies of flux tube fragmentation in two dimensions. We find that if the initial amount of field line twist is below a critical value, the degree of fragmentation at the apex of a rising Omega-loop depends on its three-dimensional geometry: the greater the apex curvature of a given Omega-loop, the lesser the degree of fragmentation of the loop as it approaches the photosphere. Thus, the amount of initial twist necessary for the loop to retain its cohesion can be reduced substantially from the two-dimensional limit. The simulations also suggest that as a fragmented flux tube emerges through a relatively quiet portion of the solar disk, extended crescent-shaped magnetic features of opposite polarity should form and steadily recede from one another. These features eventually coalesce after the fragmented portion of the Omega-loop emerges through the photosphere.Comment: 17 pages, 17 figures, uses AAS LaTeX macros v5.0. ApJ, in pres

Radiative Cooling in MHD Models of the Quiet Sun Convection Zone and Corona

We present a series of numerical simulations of the quiet Sun plasma threaded by magnetic fields that extend from the upper convection zone into the low corona. We discuss an efficient, simplified approximation to the physics of optically thick radiative transport through the surface layers, and investigate the effects of convective turbulence on the magnetic structure of the Sun's atmosphere in an initially unipolar (open field) region. We find that the net Poynting flux below the surface is on average directed toward the interior, while in the photosphere and chromosphere the net flow of electromagnetic energy is outward into the solar corona. Overturning convective motions between these layers driven by rapid radiative cooling appears to be the source of energy for the oppositely directed fluxes of electromagnetic energy.Comment: 20 pages, 5 figures, Solar Physics, in pres

Convective Dynamos and the Minimum X-ray Flux in Main Sequence Stars

The objective of this paper is to investigate whether a convective dynamo can account quantitatively for the observed lower limit of X-ray surface flux in solar-type main sequence stars. Our approach is to use 3D numerical simulations of a turbulent dynamo driven by convection to characterize the dynamic behavior, magnetic field strengths, and filling factors in a non-rotating stratified medium, and to predict these magnetic properties at the surface of cool stars. We use simple applications of stellar structure theory for the convective envelopes of main-sequence stars to scale our simulations to the outer layers of stars in the F0--M0 spectral range, which allows us to estimate the unsigned magnetic flux on the surface of non-rotating reference stars. With these estimates we use the recent results of \citet{Pevtsov03} to predict the level of X-ray emission from such a turbulent dynamo, and find that our results compare well with observed lower limits of surface X-ray flux. If we scale our predicted X-ray fluxes to \ion{Mg}{2} fluxes we also find good agreement with the observed lower limit of chromospheric emission in K dwarfs. This suggests that dynamo action from a convecting, non-rotating plasma is a viable alternative to acoustic heating models as an explanation for the basal emission level seen in chromospheric, transition region, and coronal diagnostics from late-type stars.Comment: ApJ, accepted, 30 pages with 7 figure

Is null-point reconnection important for solar flux emergence?

The role of null-point reconnection in a 3D numerical MHD model of solar emerging flux is investigated. The model consists of a twisted magnetic flux tube rising through a stratified convection zone and atmosphere to interact and reconnect with a horizontal overlying magnetic field in the atmosphere. Null points appear as the reconnection begins and persist throughout the rest of the emergence, where they can be found mostly in the model photosphere and transition region, forming two loose clusters on either side of the emerging flux tube. Up to 26 nulls are present at any one time, and tracking in time shows that there is a total of 305 overall, despite the initial simplicity of the magnetic field configuration. We find evidence for the reality of the nulls in terms of their methods of creation and destruction, their balance of signs, their long lifetimes, and their geometrical stability. We then show that due to the low parallel electric fields associated with the nulls, null-point reconnection is not the main type of magnetic reconnection involved in the interaction of the newly emerged flux with the overlying field. However, the large number of nulls implies that the topological structure of the magnetic field must be very complex and the importance of reconnection along separators or separatrix surfaces for flux emergence cannot be ruled out.Comment: 26 pages, 12 figures. Added one referenc

Can we Determine Electric Fields and Poynting Fluxes from Vector Magnetograms and Doppler Measurements?

The availability of vector magnetogram sequences with sufficient accuracy and cadence to estimate the time derivative of the magnetic field allows us to use Faraday's law to find an approximate solution for the electric field in the photosphere, using a Poloidal-Toroidal Decomposition (PTD) of the magnetic field and its partial time derivative. Without additional information, however, the electric field found from this technique is under-determined -- Faraday's law provides no information about the electric field that can be derived the gradient of a scalar potential. Here, we show how additional information in the form of line-of-sight Doppler flow measurements, and motions transverse to the line-of-sight determined with ad-hoc methods such as local correlation tracking, can be combined with the PTD solutions to provide much more accurate solutions for the solar electric field, and therefore the Poynting flux of electromagnetic energy in the solar photosphere. Reliable, accurate maps of the Poynting flux are essential for quantitative studies of the buildup of magnetic energy before flares and coronal mass ejections.Comment: Solar Physics, in press. 14 pages, 3 figure

New Insights into White-Light Flare Emission from Radiative-Hydrodynamic Modeling of a Chromospheric Condensation

(abridged) The heating mechanism at high densities during M dwarf flares is poorly understood. Spectra of M dwarf flares in the optical and near-ultraviolet wavelength regimes have revealed three continuum components during the impulsive phase: 1) an energetically dominant blackbody component with a color temperature of T $\sim$ 10,000 K in the blue-optical, 2) a smaller amount of Balmer continuum emission in the near-ultraviolet at lambda $<$ 3646 Angstroms and 3) an apparent pseudo-continuum of blended high-order Balmer lines. These properties are not reproduced by models that employ a typical "solar-type" flare heating level in nonthermal electrons, and therefore our understanding of these spectra is limited to a phenomenological interpretation. We present a new 1D radiative-hydrodynamic model of an M dwarf flare from precipitating nonthermal electrons with a large energy flux of $10^{13}$ erg cm$^{-2}$ s$^{-1}$. The simulation produces bright continuum emission from a dense, hot chromospheric condensation. For the first time, the observed color temperature and Balmer jump ratio are produced self-consistently in a radiative-hydrodynamic flare model. We find that a T $\sim$ 10,000 K blackbody-like continuum component and a small Balmer jump ratio result from optically thick Balmer and Paschen recombination radiation, and thus the properties of the flux spectrum are caused by blue light escaping over a larger physical depth range compared to red and near-ultraviolet light. To model the near-ultraviolet pseudo-continuum previously attributed to overlapping Balmer lines, we include the extra Balmer continuum opacity from Landau-Zener transitions that result from merged, high order energy levels of hydrogen in a dense, partially ionized atmosphere. This reveals a new diagnostic of ambient charge density in the densest regions of the atmosphere that are heated during dMe and solar flares.Comment: 50 pages, 2 tables, 13 figures. Accepted for publication in the Solar Physics Topical Issue, "Solar and Stellar Flares". Version 2 (June 22, 2015): updated to include comments by Guest Editor. The final publication is available at Springer via http://dx.doi.org/10.1007/s11207-015-0708-