Unstable magnetohydrodynamical continuous spectrum of accretion disks. A new route to magnetohydrodynamical turbulence in accretion disks

We present a detailed study of localised magnetohydrodynamical (MHD) instabilities occuring in two--dimensional magnetized accretion disks. We model axisymmetric MHD disk tori, and solve the equations governing a two--dimensional magnetized accretion disk equilibrium and linear wave modes about this equilibrium. We show the existence of novel MHD instabilities in these two--dimensional equilibria which do not occur in an accretion disk in the cylindrical limit. The disk equilibria are numerically computed by the FINESSE code. The stability of accretion disks is investigated analytically as well as numerically. We use the PHOENIX code to compute all the waves and instabilities accessible to the computed disk equilibrium. We concentrate on strongly magnetized disks and sub--Keplerian rotation in a large part of the disk. These disk equilibria show that the thermal pressure of the disk can only decrease outwards if there is a strong gravitational potential. Our theoretical stability analysis shows that convective continuum instabilities can only appear if the density contours coincide with the poloidal magnetic flux contours. Our numerical results confirm and complement this theoretical analysis. Furthermore, these results show that the influence of gravity can either be stabilizing or destabilizing on this new kind of MHD instability. In the likely case of a non--constant density, the height of the disk should exceed a threshold before this type of instability can play a role. This localised MHD instability provides an ideal, linear route to MHD turbulence in strongly magnetized accretion disk tori.Comment: 20 pages, 10 figures, accepted for publication in Astronomy & Astrophysic

Toward detailed prominence seismology - II. Charting the continuous magnetohydrodynamic spectrum

Starting from accurate MHD flux rope equilibria containing prominence condensations, we initiate a systematic survey of their linear eigenoscillations. To quantify the full spectrum of linear MHD eigenmodes, we require knowledge of all flux-surface localized modes, charting out the continuous parts of the MHD spectrum. We combine analytical and numerical findings for the continuous spectrum for realistic prominence configurations. The equations governing all eigenmodes for translationally symmetric, gravitating equilibria containing an axial shear flow, are analyzed, along with their flux-surface localized limit. The analysis is valid for general 2.5D equilibria, where either density, entropy, or temperature vary from one flux surface to another. We analyze the mode couplings caused by the poloidal variation in the flux rope equilibria, by performing a small gravity parameter expansion. We contrast the analytical results with continuous spectra obtained numerically. For equilibria where the density is a flux function, we show that continuum modes can be overstable, and we present the stability criterion for these convective continuum instabilities. Furthermore, for all equilibria, a four-mode coupling scheme between an Alfvenic mode of poloidal mode number m and three neighboring (m-1, m, m+1) slow modes is identified, occurring in the vicinity of rational flux surfaces. For realistically prominence equilibria, this coupling is shown to play an important role, from weak to stronger gravity parameter g values. The analytic predictions for small g are compared with numerical spectra, and progressive deviations for larger g are identified. The unstable continuum modes could be relevant for short-lived prominence configurations. The gaps created by poloidal mode coupling in the continuous spectrum need further analysis, as they form preferred frequency ranges for global eigenoscillations.Comment: Accepted by Astronmy & Astrophysics, 21 pages, 15 figure

Toward detailed prominence seismology - I. Computing accurate 2.5D magnetohydrodynamic equilibria

Context. Prominence seismology exploits our knowledge of the linear eigenoscillations for representative magnetohydro- dynamic models of filaments. To date, highly idealized models for prominences have been used, especially with respect to the overall magnetic configurations. Aims. We initiate a more systematic survey of filament wave modes, where we consider full multi-dimensional models with twisted magnetic fields representative of the surrounding magnetic flux rope. This requires the ability to compute accurate 2.5 dimensional magnetohydrodynamic equilibria that balance Lorentz forces, gravity, and pressure gradients, while containing density enhancements (static or in motion). Methods. The governing extended Grad-Shafranov equation is discussed, along with an analytic prediction for circular flux ropes for the Shafranov shift of the central magnetic axis due to gravity. Numerical equilibria are computed with a finite element-based code, demonstrating fourth order accuracy on an explicitly known, non-trivial test case. Results. The code is then used to construct more realistic prominence equilibria, for all three possible choices of a free flux-function. We quantify the influence of gravity, and generate cool condensations in hot cavities, as well as multi- layered prominences. Conclusions. The internal flux rope equilibria computed here have the prerequisite numerical accuracy to allow a yet more advanced analysis of the complete spectrum of linear magnetohydrodynamic perturbations, as will be demonstrated in the companion paper.Comment: Accepted by Astronomy & Astrophysics, 15 pages, 15 figure

Magneto-rotational overstability in accretion disks

We present analytical and numerical studies of magnetorotational instabilities occuring in magnetized accretion disks. In these studies we make use of the linearised compressible MHD equations. These calculations are performed for general radially stratified disks in the cylindrical limit. In particular, we investigate the influence of nonvanishing toroidal magnetic field component on the growth rate and oscillation frequency of magnetorotational instabilities in Keplerian disks. We find the persistence of these instabilities in accretion disks close to equipartition. Our calculations show that these eigenmodes become overstable (complex eigenvalue), due to the presence of a toroidal magnetic field component, while their growth rate reduces slightly. Furthermore, we demonstrate the presence of magneto-rotational overstabilities in weakly magnetized sub-Keplerian rotating disks. We show that the growth rate scales with the rotation frequency of the disk. These eigenmodes also have a nonzero oscillation frequency, due to the presence of the dominant toroidal magnetic field component. The overstable character of the MRI increases as the rotation frequency of the disk decreases.Comment: 11 pager, 18 Postscript figures, accepted for publication in Astronomy & Astrophysic

Magnetohydrostatic solar prominences in near-potential coronal magnetic fields

We present numerical magnetohydrostatic solutions describing the gravitationally stratified, bulk equilibrium of cool, dense prominence plasma embedded in a near-potential coronal field. These solutions are calculated using the FINESSE magnetohydrodynamics equilibrium solver and describe the morphologies of magnetic field distributions in and around prominences and the cool prominence plasma that these fields support. The equilibrium condition for this class of problem is usually different in distinct subdomains, separated by free boundaries, across which solutions are matched by suitable continuity or jump conditions describing force balance. We employ our precise finite element elliptic solver to calculate solutions not accessible by previous analytical techniques with temperature or entropy prescribed as free functions of the magnetic flux function, including a range of values of the polytropic index, temperature variations mainly across magnetic field lines and photospheric field profiles sheared close to the polarity inversion line. Out of the many examples computed here, perhaps the most noteworthy is one which reproduces precisely the three-part structure often encountered in observations: a cool dense prominence within a cavity/flux rope embedded in a hot corona. The stability properties of these new equilibria, which may be relevant to solar eruptions, can be determined in the form of a full resistive MHD spectrum using a companion hyperbolic stability solver.Comment: To appear in ApJ August 200

Double-resonant fast particle-wave interaction

In future fusion devices fast particles must be well confined in order to transfer their energy to the background plasma. Magnetohydrodynamic instabilities like Toroidal Alfv\'en Eigenmodes or core-localized modes such as Beta Induced Alfv\'en Eigenmodes and Reversed Shear Alfv\'en Eigenmodes, both driven by fast particles, can lead to significant losses. This is observed in many ASDEX Upgrade discharges. The present study applies the drift-kinetic HAGIS code with the aim of understanding the underlying resonance mechanisms, especially in the presence of multiple modes with different frequencies. Of particular interest is the resonant interaction of particles simultaneously with two different modes, referred to as 'double-resonance'. Various mode overlapping scenarios with different q profiles are considered. It is found that, depending on the radial mode distance, double-resonance is able to enhance growth rates as well as mode amplitudes significantly. Surprisingly, no radial mode overlap is necessary for this effect. Quite the contrary is found: small radial mode distances can lead to strong nonlinear mode stabilization of a linearly dominant mode.Comment: 12 pages, 11 figures; Nuclear Fusion 52 (2012

The Role of Plasticity and Adaptation in the Incipient Speciation of a Fire Salamander Population

Phenotypic plasticity and local adaptation via genetic change are two major mechanisms of response to dynamic environmental conditions. These mechanisms are not mutually exclusive, since genetic change can establish similar phenotypes to plasticity. This connection between both mechanisms raises the question of how much of the variation observed between species or populations is plastic and how much of it is genetic. In this study, we used a structured population of fire salamanders (Salamandra salamandra), in which two subpopulations differ in terms of physiology, genetics, mate-, and habitat preferences. Our goal was to identify candidate genes for differential habitat adaptation in this system, and to explore the degree of plasticity compared to local adaptation. We therefore performed a reciprocal transfer experiment of stream- and pond-originated salamander larvae and analyzed changes in morphology and transcriptomic profile (using species-specific microarrays). We observed that stream- and pond-originated individuals diverge in morphology and gene expression. For instance, pond-originated larvae have larger gills, likely to cope with oxygen-poor ponds. When transferred to streams, pond-originated larvae showed a high degree of plasticity, resembling the morphology and gene expression of stream-originated larvae (reversion); however the same was not found for stream-originated larvae when transferred to ponds, where the expression of genes related to reduction-oxidation processes was increased, possibly to cope with environmental stress. The lack of symmetrical responses between transplanted animals highlights the fact that the adaptations are not fully plastic and that some level of local adaptation has already occurred in this population. This study illuminates the process by which phenotypic plasticity allows local adaptation to new environments and its potential role in the pathway of incipient speciation

Numerical simulations of kink instability in line-tied coronal loops

The results from numerical simulations carried out using a new shock-capturing, Lagrangian-remap, 3D MHD code, Lare3d are presented. We study the evolution of the m=1 kink mode instability in a photospherically line-tied coronal loop that has no net axial current. During the non-linear evolution of the kink instability, large current concentrations develop in the neighbourhood of the infinite length mode rational surface. We investigate whether this strong current saturates at a finite value or whether scaling indicates current sheet formation. In particular, we consider the effect of the shear, defined by where is the fieldline twist of the loop, on the current concentration. We also include a non-uniform resistivity in the simulations and observe the amount of free magnetic energy released by magnetic reconnection