Search

Coronal-hole jets occur ubiquitously in the Sun's coronal holes, at EUV and X-ray bright points associated with intrusions of minority magnetic polarity. The embedded-bipole model for these jets posits that they are driven by explosive, fast reconnection between the stressed closed field of the embedded bipole and the open field of the surrounding coronal hole. Previous numerical studies in Cartesian geometry, assuming uniform ambient magnetic field and plasma while neglecting gravity and solar wind, demonstrated that the model is robust and can produce jet-like events in simple configurations. We have extended these investigations by including spherical geometry,gravity, and solar wind in a nonuniform, coronal hole-like ambient atmosphere. Our simulations confirm that the jet is initiated by the onset of a kink-like instability of the internal closed field, which induces a burst of reconnection between the closed and external open field, launching a helical jet. Our new results demonstrate that the jet propagation is sustained through the outer corona, in the form of a traveling nonlinear Alfven wave front trailed by slower-moving plasma density enhancements that are compressed and accelerated by the wave. This finding agrees well with observations of white-light coronal-hole jets, and can explain microstreams and torsional Alfven waves detected in situ in the solar wind. We also use our numerical results to deduce scaling relationships between properties of the coronal source region and the characteristics of the resulting jet, which can be tested against observations.

The mechanism that accelerates particles to the energies required to produce the observed high-energy impulsive emission in solar flares is not well understood. Drake et al. proposed a mechanism for accelerating electrons in contracting magnetic islands formed by kinetic reconnection in multi-layered current sheets (CSs). We apply these ideas to sunward-moving flux ropes (2.5D magnetic islands) formed during fast reconnection in a simulated eruptive flare. A simple analytic model is used to calculate the energy gain of particles orbiting the field lines of the contracting magnetic islands in our ultrahigh-resolution 2.5D numerical simulation. We find that the estimated energy gains in a single island range up to a factor of five. This is higher than that found by Drake et al. for islands in the terrestrial magnetosphere and at the heliopause, due to strong plasma compression that occurs at the flare CS. In order to increase their energy by two orders of magnitude and plausibly account for the observed high-energy flare emission, the electrons must visit multiple contracting islands. This mechanism should produce sporadic emission because island formation is intermittent. Moreover, a large number of particles could be accelerated in each magneto hydro dynamic-scale island, which may explain the inferred rates of energetic-electron production in flares. We conclude that island contraction in the flare CS is a promising candidate for electron acceleration in solar eruptions.

The mechanism that accelerates particles to the energies required to produce the observed high-energy emission in solar flares is not well understood. Drake et al. (2006) proposed a kinetic mechanism for accelerating electrons in contracting magnetic islands formed by reconnection. In this model, particles that gyrate around magnetic field lines transit from island to island, increasing their energy by Fermi acceleration in those islands that are contracting. Based on these ideas, we present an analytic model to estimate the energy gain of particles orbiting around field lines inside a flux rope (2.5D magnetic island). We calculate the change in the velocity of the particles as the flux rope evolves in time. The method assumes a simple profile for the magnetic field of the evolving island; it can be applied to any case where flux ropes are formed. In our case, the flux-rope evolution is obtained from our recent high-resolution, compressible 2.5D MHD simulations of breakout eruptive flares. The simulations allow us to resolve in detail the generation and evolution of large-scale flux ropes as a result of sporadic and patchy reconnection in the flare current sheet. Our results show that the initial energy of particles can be increased by 2-5 times in a typical contracting island, before the island reconnects with the underlying arcade. Therefore, particles need to transit only from 3-7 islands to increase their energies by two orders of magnitude. These macroscopic regions, filled with a large number of particles, may explain the large observed rates of energetic electron production in flares. We conclude that this mechanism is a promising candidate for electron acceleration in flares, but further research is needed to extend our results to 3D flare conditions.

Context. Jets are dynamic, impulsive, well collimated plasma events developing at many different scales and in differentlayers of the solar atmosphere.Aims. Jets are believed to be induced by magnetic reconnection a process central to many astrophysical phenomena.Studying their dynamics can help us to better understand the processes acting in larger eruptive events (e.g. flares andcoronal mass ejection) as well as mass, magnetic helicity and energy transfer at all scale in the solar atmosphere. Therelative simplicity of their magnetic geometry and topology, compared with larger solar active events, makes jets idealcandidate for studying the fundamental role of reconnection in energetic events. Methods. In this study, using our state-of-the-art numerical solver ARMS, we present several parametric studies of anumerical three-dimensional magneto-hydrodynamic model of solar jet-like events. We study the impact of the magneticfield inclination and photospheric field distribution on the generation and properties of two types of observed solar jets:the so-called standard and blowout jets.Results. The present parametric studies validate our model of jets for different geometric properties of the magneticconfiguration. On one hand, a blowout jets is eventually always trigger for the range of parameters that we tested. The3D magnetic null-point configuration is a very robust structure for the energy storage and impulsive release characteristicof blowout jets. On the other hand, the existence of the standard jet depends on the magnetic geometry. We show thatthe reconnection occurring during the standard jet phase influences the properties of the trigger of the blowout jet.Conclusions. Our results allow us to better understand the energization, triggering and driving processes of standardand blowout jets. Our model allows us to predict the impulsiveness and energetics of jets in terms of the surroundingmagnetic field configuration. Finally we discuss the interpretation of standard jets and blowout jets and the physicalfactors that determine which type of jets will occur.

We investigate the process of formation and subsequent evolution of prominence plasma in a filament channel and its overlying arcade. We construct a three-dimensional time-dependent model of a filament-channel prominence suitable to be compared with observations. We combine this magnetic field structure with one-dimensional independent simulations of many flux tubes. The magnetic structure is a three-dimensional sheared double arcade, and the thermal non-equilibrium process governs the plasma evolution. We have found that the condensations in the corona can be divided into two populations: threads and blobs. Threads are massive condensations that linger in the field line dips. Blobs are ubiquitous small condensations that are produced throughout the filament and overlying arcade magnetic structure, and rapidly fall to the chromosphere. The total prominence mass is in agreement with observations. The threads are the principal contributors to the total mass, whereas the blob contribution is small. The motion of the threads is basically horizontal, while blobs move in all directions along the field. The peak velocities for both populations are comparable, but there is a weak tendency for the velocity to increase with the inclination, and the blobs with motion near vertical have the largest values of the velocity. We have generated synthetic images of the whole structure in an H proxy and in two EUV channels of the AIA instrument aboard SDO. These images show the plasma at cool, warm and hot temperatures. The theoretical differential emission measure of our system agrees very well with observations in the temperature range log T = 4.6-5.7. We conclude that the sheared-arcade magnetic structure and plasma dynamics fit well the abundant observational evidence.

The largest and most efficient particle accelerators in the solar system are the giant events consisting of a fast coronal mass ejection (CME) and an intense X-class solar flare. Both flares and CMEs can produce l0(exp 32) ergs or more in nonthermal particles. Two general processes are believed to be responsible: particle acceleration at the strong shock ahead of the CME, and reconnection-driven acceleration in the flare current sheet. Although shock acceleration is relatively well understood, the mechanism by which flare reconnection produces nonthermal particles is still an issue of great debate. We address the question of CME/flare particle acceleration in the context of the breakout model using 2.5D MHD simulations with adaptive mesh refinement (AMR). The AMR capability allows us to achieve ultra-high numerical resolution and, thereby, determine the detailed structure and dynamics of the flare reconnection region. Furthermore, we employ newly developed numerical analysis tools for identifying and characterizing magnetic nulls, so that we can quantify accurately the number and location of magnetic islands during reconnection. Our calculations show that flare reconnection is dominated by the formation of magnetic islands. In agreement with many other studies, we find that the number of islands scales with the effective Lundquist number. This result supports the recent work by Drake and co-workers that postulates particle acceleration by magnetic islands. On the other hand, our calculations also show that the flare reconnection region is populated by numerous shocks and other indicators of strong turbulence, which can also accelerate particles. We discuss the implications of our calculations for the flare particle acceleration mechanism and for observational tests of the models.

The integration of kinetic effects into macroscopic numerical models is currently of great interest to the heliophysics community, particularly in the context of magnetic reconnection. Reconnection governs the large-scale energy release and topological rearrangement of magnetic fields in a wide variety of laboratory, heliophysical, and astrophysical systems. We are examining the formation and reconnection of current sheets in a simple, two-dimensional X-line configuration using high-resolution particle-in-cell (PIC) simulations. The initial minimum-energy, potential magnetic field is perturbed by excess thermal pressure introduced into the particle distribution function far from the X line. Subsequently, the relaxation of this added stress leads self-consistently to the development of a current sheet that reconnects for imposed stress of sufficient strength. We compare the time-dependent evolution and final state of our PIC simulations with macroscopic magnetohydrodynamic simulations assuming both uniform and localized electrical resistivities (C. R. DeVore et al., this meeting), as well as with force-free magnetic-field equilibria in which the amount of reconnection across the X line can be constrained to be zero (ideal evolution) or optimal (minimum final magnetic energy). We will discuss implications of our results for understanding magnetic-reconnection onset and cessation at kinetic scales in dynamically formed current sheets, such as those occurring in the solar corona and terrestrial magnetotail.

The integration of kinetic effects into macroscopic numerical models is currently of great interest to the plasma physics community, particularly in the context of magnetic reconnection. We are examining the formation and reconnection of current sheets in a simple, two-dimensional X-line configuration using high resolution particle-in-cell (PIC) simulations. The initial potential magnetic field is perturbed by thermal pressure introduced into the particle distribution far from the X line. The relaxation of this added stress leads to the development of a current sheet, which reconnects for imposed stress of sufficient strength. We compare the evolution and final state of our PIC simulations with magnetohydrodynamic simulations assuming both uniform and localized resistivities, and with force-free magnetic-field equilibria in which the amount of reconnect ion across the X line can be constrained to be zero (ideal evolution) or optimal (minimum final magnetic energy). We will discuss implications of our results for reconnection onset and cessation at kinetic scales in dynamically formed current sheets, such as those occurring in the terrestrial magnetotail and solar corona.

For understanding and eventually predicting coronal mass ejections/eruptive flares, two critical questions must be answered: What is the mechanism for eruption onset, and what is the mechanism for the rapid acceleration? We address these questions in the context of the breakout model using 2.5D MHD simulations with adaptive mesh refinement (AMR). The AMR capability allowed us to achieve ultra-high numerical resolution and, thereby, determine the influence of the effective Lundquist number on the eruption. Our calculations show that, at least, for the breakout model, the onset of reconnection external to the highly sheared filament channel is the onset mechanism. Once this reconnection turns on, eruption is inevitable. However, as long as this is the only reconnection in the system, the eruption remains slow. We find that the eruption undergoes an abrupt "take-off" when the flare reconnection below the erupting plasmoid develops significant reconnection jets. We conclude that in fast CMEs, flare reconnection is the primary mechanism responsible for both flare heating and CME acceleration. We discuss the implications of these results for SDO observations and describe possible tests of the model.

A brief review is presented of work done over the last eight years investigating the fundamental physics of plasmas and magnetic fields under conditions similar to those that are thought to be present in the outer layers of the solar atmosphere, including the transition region and the corona. The models used to study the coronal structures and the thermal instability in the solar atmosphere are discussed. The results of studies of magnetic energy release in the corona and MHD turbulence in the solar wind are examined.

Spectral transfer has been proposed as the primary mechanism for generating outward-propagating Alfven waves in the solar wind. This process has been investigated extensively for imcompressible magnetofluids, but the issue of whether it occurs in compressible magnetofluids such as the solar wind remains unresolved. The results of direct numerical simulations of nonisentropic-compressible two-dimensional MHD turbulence indicate that, for systems with finite initial cross helicity, the correlation between the fluctuating velocity field and the fluctuating magnetic field grows as a function of time. This growth of correlation can be interpreted as a turbulent process, as shown by examination of modal wavenumber spectra.

The nonlinear evolution of a radiation-driven thermally unstable planar fluid is simulated numerically using a semiimplicit finite-difference algorithm. When the equilibrium state of the fluid is perturbed by random initial excitation of the velocity field, dense, cool, two-dimensional structures are found to form in a rarer, warmer surrounding medium. The nonlinear phase of evolution is characterized by the turbulent contraction of the condensed region, accompanied by a significant increase in the amount of energy radiated. It is found that, if the random velocity perturbation has a sufficiently large amplitude, the fluid will not form condensed structures. Finally, the relationship of these results to observations of the solar chromosphere, transition region, and corona is discussed.

A characterization is given of the preflare magnetic field, using theoretical models of force free fields together with observed field structure to determine the general morphology. Direct observational evidence for sheared magnetic fields is presented. The role of this magnetic shear in the flare process is considered within the context of a MHD model that describes the buildup of magnetic energy, and the concept of a critical value of shear is explored. The related subject of electric currents in the preflare state is discussed next, with emphasis on new insights provided by direct calculations of the vertical electric current density from vector magnetograph data and on the role of these currents in producing preflare brightenings. Results from investigations concerning velocity fields in flaring active regions, describing observations and analyses of preflare ejecta, sheared velocities, and vortical motions near flaring sites are given. This is followed by a critical review of prevalent concepts concerning the association of flux emergence with flares

Emerging magnetic flux plays an important role in the development of active regions on the sun and, perhaps, in the subsequent activation of flares. However, the energy input that produces preflare brightenings and flares probably does not come from the flux emergence itself but from one or more associated energy-releasing processes - likely candidates include magnetic reconnection and various current-driven plasma micro-instabilities. Here the interplay between the changing physical characteristics of an emerging magnetic-flux tube and the onset and evolution of a representative 'bump-on-tail' plasma current-driven instability is investigated. The microinstability heats the ambient material, thus changing the macroscopic characteristics of the plasma in which the model and current-driven instability occurs.

The large flare of November 7, 1980 provides a unique opportunity to investigate the upward-moving plasma seen during the early stages of many flares. Soft X-ray spectroscopic data obtained by the Solar Flare X-ray (SOLFLEX) instruments on board the Air Force P78-1 satellite have been used to determine the spatial extent, turbulent velocity, temperature, and emission measure of the blueshifted and stationary plasmas, as well as the upward velocity of the blueshifted component alone. Two geometries are considered in calculating the resultant mass and energy balance. In addition, coincident hard X-ray data was acquired from the HXRBS instrument on board the SMM satellite to determine the relative timing and enertics of the hard and soft X-ray flare plasmas. These results are compared with the predictions of the chromospheric evaporation hypothesis. It is concluded that electron-induced evaporation plays a minor role in this flare, and that another mechanism must account for the observed blueshifted emission.

Coronal bullets are small ejecta of cool, dense plasma observed to accelerate through the solar atmosphere from 20 to 450 km/s. The NRL Dynamic Flux Tube Model has been used to simulate the evolving physical properties of these dynamic events. The present calculations utilize an adaptive-gridding technique to resolve the fine structure within and around the bullets. In this work, an identification was made of a component of shocked plasma which piles up ahead of the bullet and eventually dominates both the dynamics and heating of the original bullet mass. The observational consequences of this shocked component are discussed in terms of the available HRTS EUV data, and suggestions are made for optimizing future observations of this phenomenon. An investigation has also been conducted of the structure of the bullet material visible in EUV spectral lines and the observable characteristics of the EUV-emitting plasma. Finally, the most likely mechanisms for accelerating the bullets, as well as favorable sites of origin are evaluated.

A numerical model is defined for asymmetric full solar flare loop heating and comparisons are made with observational data. The Dynamic Flux Tube Model is used to describe the heating process in terms of one-dimensional, two fluid conservation equations of mass, energy and momentum. An adaptive grid allows for the downward movement of the transition region caused by an advancing conduction front. A loop 20,000 km long is considered, along with a flare heating system and the hydrodynamic evolution of the loop. The model was applied to generating line profiles and spatial X-ray and UV line distributions, which were compared with SMM, P78-1 and Hintori data for Fe, Ca and Mg spectra. Little agreement was obtained, and it is suggested that flares be treated as multi-loop phenomena. Finally, it is concluded that chromospheric evaporation is not an effective mechanism for generating the soft X-ray bursts associated with flares.

Results are presented of computer simulations of the jets that accelerate through the corona at velocities of 50 to 400 km/s. Particular emphasis is placed on the sensitivity of the induced acceleration to the form in which energy is put into the system. A comparison is made between the observed and predicted physical characteristics of the high-velocity bullets; the potential contribution of the bullets to the mass and energy balance of the solar corona is considered. It is found that the velocity and temperature evolution of the bullets can be modeled successfully by assuming energy input in the form of an external force, pushing continuously on the ejected material. From the physical characteristics of the model bullets and the energy input required to reproduce the observations, it is concluded that the bullets may constitute a significant fraction of the coronal mass flux but only a negligible component of the coronal energy budget.

An explanation in terms of betatron acceleration is attempted for two prominent features of solar hard X-ray emission profiles: the multiplicity manifested by complex impulsive events, and the intense, long-enduring radiation typical of two-stage or extended events. The model employs the basic characteristics of a simple magnetic trap. The observed spectral evolution can be represented effectively by a set of characteristics derived in a thick-target representation. The application to multiple-spike and two-stage bursts is studied, and it is found that impulsive events show no signs of betatron acceleration. However, during the first few minutes of the second stage, the nonimpulsive X-ray emission appears to be consistent with the predictions of the model. Several intense extended events are analyzed and the betatron process is established as a viable candidate for producing the early second-stage emission.

The gradual soft X-ray emissions associated with a homogeneous set of solar flares have been investigated in the context of a thermal model proposed to explain the impulsive components. The parametric techniques which successfully characterized the hard X-ray and microwave observations are employed in an event-by-event analysis to test for quantitative and correlative relationships between the impulsive and gradual emissions. The results of this investigation are consistent with the hypothesis that the hard X-ray and microwave emissions are produced by bulk heating of a common thermal source. The quantitative relationships require an additional source to explain the soft X-ray observations, consistent with previous results. Correlations between the energetics of the impulsive and gradual emissions, identified in the present work, provide the first clear evidence that their energizing mechanisms are related.

A set of solar hard X-ray bursts observed with the hard X-ray burst spectrometer on board the OSO-5 satellite was analyzed. The multiply-impulsive two stage events were selected on the basis of both morphological characteristics and association with appropriate phenomena at other wavelengths. Coincident radio, soft X-ray, H-alpha interplanetary particle, and magnetographic data were obtained from several observatories, to aid in developing a comprehensive picture of the physical processes underlying these complex bursts. Two classes of multiply impulsive bursts were identified: events whose components spikes apparently originate in one location, and events in which groups of spikes appear to come from separate regions which flare sequentially. The origin of multiplicity in the case of a single source region remains unidentified. Purely impulsive emissions show no sign of betatron acceleration, thus eliminating this mechanisn as a candidate for inducing multiply spiked structure. The majority of the two stage bursts, however, exhibited spectral behavior consistent with the betatron model, for the first few minutes of the second stage. Betatron acceleration thus has been identified as a common second stage phenomenon.

It is proposed that the source for solar hard X-ray bursts may consist of a distribution of many small impulsively-heated kernels, each cooled by anomalous conduction, with lifetimes shorter than current burst data temporal resolution. The dynamic spectra of bursts are governed by the dynamic evolution of the kernel production process, such as magnetic-field dissipation in the tearing mode. An integral equation is formulated, the solution of which yields information on this kernel production process from dynamic burst spectra, for any kernel model.

Some results from the analysis of a set of multiply impulsive hard X-ray and microwave solar bursts are presented, showing that some bursts can exhibit widely different magnetic-field strengths at different times. Two categories of microwave spectral behavior are identified: those events during which the microwave turnover frequency (MTF) and spectral shape (SS) remain the same from peak to peak, and those during which the MTF and SS change significantly. These categories correspond to two classes of multiply impulsive bursts: those for which the emission can be characterized by a constant magnetic field and therefore a single source region, and those in which groups of component spikes appear to originate in regions of different magnetic-field strengths, corresponding to separate source regions which flare sequentially. With regard to the latter type, examples are presented, the discrete flaring regions are examined, and their spatial separations are estimated.

The possible production of Li-7 by nuclear reactions in UV Ceti flares has been considered. By utilizing solar observations and theory, a relationship is derived between flare energy and production rates for Li-7; approximately 100 erg of total flare energy is found to denote the formation of a Li-7 atom. Based on this value and best estimates of UV Ceti-type flare rates, it is concluded that less than 10% of the Li-7 observed in the intestellar medium may have been produced by this mechanism. Formation of significant amounts of interstellar deuterium by this method is ruled out.

Hard X-ray and microwave observations of multiply-impulsive solar bursts, identified in the OSO-5 data were analyzed. Spectra in both frequency ranges were used to determine whether or not the source properties change from peak to peak within individual bursts. Two categories of microwave spectral behavior were identified: those events during which the microwave turnover frequency and spectral shape remain the same from peak to peak, and those during which the turnover frequency and spectral shape change significantly. These categories correspond to two classes of multiply-impulsive bursts: those for which the emission can be characterized by a constant magnetic field and therefore a single source region, in which case the multiplicity may be due to modulation of the emission process; and those in which groups of component spikes appear to originate in regions of different magnetic-field strengths, corresponding to separate source regions which flare sequentially. Examples of the latter type of events are presented. The discrete flaring regions are analyzed and their spatial separations estimated.

Coordinated X ray, optical, and radio observations of the flare star YZ CMi are reported. Twenty-two minor optical flares and twelve radio events were recorded. No major optical flares, greater than 3 magnitudes, were observed. Although no flare related X ray emission was observed, the measured upper limits in this band enable meaningful comparisons with published flare star models. Three of the five models predicting the relative X ray to optical or radio flare luminosities are in serious disagreement with the observations. For the largest optical flare with coincident X ray coverage, the 3 sigma upper limit on X ray emission in the 0.15 to 0.8 keV band is 8.7 x 10 to the 28th power erg/s, corresponding to a ratio of X ray to B-band luminosity of less than 0.3. Based on the present results, the contribution of the flares UV Ceti flare stars to the galactic component of the diffuse soft X ray background is less than 0.2 percent.