A. de Vicente

AUTOMATIC UNSUPERVISED CLASSIFICATION OF ALL SLOAN DIGITAL SKY SURVEY DATA RELEASE 7 GALAXY SPECTRA

Using the k-means cluster analysis algorithm, we carry out an unsupervised classification of all galaxy spectra in the seventh and final Sloan Digital Sky Survey data release (SDSS/DR7). Except for the shift to restframe wavelengths, and the normalization to the g-band flux, no manipulation is applied to the original spectra. The algorithm guarantees that galaxies with similar spectra belong to the same class. We find that 99% of the galaxies can be assigned to only 17 major classes, with 11 additional minor classes including the remaining 1%. The classification is not unique since many galaxies appear in between classes, however, our rendering of the algorithm overcomes this weakness with a tool to identify borderline galaxies. Each class is characterized by a template spectrum, which is the average of all the spectra of the galaxies in the class. These low noise template spectra vary smoothly and continuously along a sequence labeled from 0 to 27, from the reddest class to the bluest class. Our Automatic Spectroscopic K-means-based (ASK) classification separates galaxies in colors, with classes characteristic of the red sequence, the blue cloud, as well as the green valley. When red sequence galaxies and green valley galaxies present emission lines, they are characteristic of AGN activity. Blue galaxy classes have emission lines corresponding to star formation regions. We find the expected correlation between spectroscopic class and Hubble type, but this relationship exhibits a high intrinsic scatter. Several potential uses of the ASK classification are identified and sketched, including fast determination of physical properties by interpolation, classes as templates in redshift determinations, and target selection in follow-up works (we find classes of Seyfert galaxies, green valley galaxies, as well as a significant number of outliers). The ASK classification is publicly accessible through various websites. Subject headings: catalogs – methods: statistical – galaxies: evolution – galaxies: fundamental parameters – galaxies: statistics

Rayleigh-Taylor instability in prominences from numerical simulations including partial ionization effects

We study the Rayleigh-Taylor instability (RTI) at a prominence-corona transition region in a non-linear regime. Our aim is to understand how the presence of neutral atoms in the prominence plasma influences the instability growth rate, and the evolution of velocity, magnetic field vector and thermodynamic parameters of turbulent drops. We perform 2.5D numerical simulations of the instability initiated by a multi-mode perturbation at the corona-prominence interface using a single-fluid MHD approach including a generalized Ohms law. The initial equilibrium configuration is purely hydrostatic and contains a homogeneous horizontal magnetic field forming an angle with the direction in which the plasma is perturbed. We analyze simulations with two different orientations of the magnetic field. For each field orientation we compare two simulations, one for the pure MHD case, and one including the ambipolar diffusion in the Ohms law (AD case). Other than that, both simulations for each field orientation are identical. The numerical results in the initial stage of the instability are compared with the analytical linear calculations. We find that the configuration is always unstable in the AD case. The growth rate of the small-scale modes in the non-linear regime is up to 50% larger in the AD case than in the purely MHD case and the average velocities of flows are a few percent larger. Significant drift momenta are found at the interface between the coronal and the prominence material at all stages of the instability, produced by the faster downward motion of the neutral component with respect to the ionized component. The differences in temperature of the bubbles between the ideal and non-ideal case are also significant, reaching 30%. There is an asymmetry between large rising bubbles and small-scale down flowing fingers, favoring the detection of upward velocities in observations.

ON THE ROBUSTNESS OF THE PENDULUM MODEL FOR LARGE-AMPLITUDE LONGITUDINAL OSCILLATIONS IN PROMINENCES

Large-amplitude longitudinal oscillations (LALOs) in prominences are spectacular manifestations of the solar activity. In such events nearby energetic disturbances induce periodic motions on filaments with displacements comparable to the size of the filaments themselves and with velocities larger than 20 km/s. The pendulum model, in which the gravity projected along a rigid magnetic field is the restoring force, was proposed to explain these events. However, it can be objected that in a realistic situation where the magnetic field reacts to the mass motion of the heavy prominence, the simplified pendulum model could be no longer valid. We have performed non-linear time-dependent numerical simulations of LALOs considering a dipped magnetic field line structure. In this work we demonstrate that for even relatively weak magnetic fields the pendulum model works very well. We therefore validate the pendulum model and show its robustness, with important implications for prominence seismology purposes. With this model it is possible to infer the geometry of the dipped field lines that support the prominence.

Simulated interaction of magnetohydrodynamic shock waves with a complex network-like region

We provide estimates of the wave energy reaching the solar chromosphere and corona in a network-like magnetic field topology, including a coronal null point. The waves are excited by an instantaneous strong subphotospheric source and propagate through the subphotosphere, photosphere, chromosphere, transition region, and corona with the plasma beta and other atmospheric parameters varying by several orders of magnitude. We compare two regimes of the wave propagation: a linear and nonlinear regime. While the amount of energy reaching the corona is similar in both regimes, this energy is transmitted at different frequencies. In both cases the dominant periods of waves at each height strongly depend on the local magnetic field topology, but this distribution is only in accordance with observations in the nonlinear case.

High-frequency waves in the corona due to null points

This work aims to understand the behavior of non-linear waves in the vicinity of a coronal null point. In previous works we have showed that high frequency waves are generated in such magnetic configuration. This paper studies those waves in detail in order to provide a plausible explanation of their generation. We demonstrate that slow magneto-acoustic shock waves generated in the chromosphere propagate through the null point and produce a train of secondary shocks that escape along the field lines. A particular combination of the shock wave speeds generates waves at a frequency of 80 mHz. We speculate that this frequency may be sensitive to the atmospheric parameters in the corona and therefore can be used to probe the structure of this solar layer.

Numerical simulations of quiet Sun magnetic fields seeded by the Biermann battery

The magnetic fields of the quiet Sun cover at any time more than 90% of its surface and their magnetic energy budget is crucial to explain the thermal structure of the solar atmosphere. One of the possible origins of these fields is the action of the local dynamo in the upper convection zone of the Sun. Existing simulations of the local solar dynamo require an initial seed field and sufficiently high spatial resolution in order to achieve the amplification of the seed field to the observed values in the quiet Sun. Here we report an alternative model of seeding based on the action of the Bierman battery effect. This effect generates a magnetic field due to the local imbalances in electron pressure in the partially ionized solar plasma. We show that the battery effect self-consistently creates from zero an initial seed field of a strength of the order of micro G, and together with dynamo amplification allows the generation of quiet Sun magnetic fields of a similar strength to those from solar observations.

MHDSTS: a new explicit numerical scheme for simulations of partially ionised solar plasma

The interaction of plasma with magnetic field in the partially ionised solar atmosphere is frequently modelled via a single-fluid approximation, which is valid for the case of a strongly coupled collisional media, such as solar photosphere and low chromosphere. Under the single-fluid formalism the main non-ideal effects are described by a series of extra terms in the generalised induction equation and in the energy conservation equation. These effects are: Ohmic diffusion, ambipolar diffusion, the Hall effect, and the Biermann battery effect. From the point of view of the numerical solution of the single-fluid equations, when ambipolar diffusion or Hall effects dominate can introduce severe restrictions on the integration time step and can compromise the stability of the numerical scheme. In this paper we introduce two numerical schemes to overcome those limitations. The first of them is known as Super Time-Stepping (STS) and it is designed to overcome the limitations imposed when the ambipolar diffusion term is dominant. The second scheme is called the Hall Diffusion Scheme (HDS) and it is used when the Hall term becomes dominant. These two numerical techniques can be used together by applying Strang operator splitting. This paper describes the implementation of the STS and HDS schemes in the single-fluid code Mancha3D. The validation for each of these schemes is provided by comparing the analytical solution with the numerical one for a suite of numerical tests.

Three-dimensional simulations of solar magneto-convection including effects of partial ionization

Over the last decades, realistic 3D radiative-MHD simulations have become the dominant theoretical tool for understanding the complex interactions between the plasma and the magnetic field on the Sun. Most of such simulations are based on approximations of magnetohydrodynamics, without directly considering the consequences of the very low degree of ionization of the solar plasma in the photosphere and bottom chromosphere. The presence of large amount of neutrals leads to a partial decoupling of the plasma and the magnetic field. As a consequence of that, a series of non-ideal effects (ambipolar diffusion, Hall effect and battery effect) arises. The ambipolar effect is the dominant one in the solar chromosphere. Here we report on the first three-dimensional realistic simulations of magneto-convection including ambipolar diffusion and battery effects. The simulations are done using the newly developed Mancha3D code. Our results reveal that ambipolar diffusion causes measurable effects on the amplitudes of waves excited by convection in the simulations, on the absorption of Poynting flux and heating and on the formation of chromospheric structures. We provide a low limit on the chromospheric temperature increase due to the ambipolar effect using the simulations with battery-excited dynamo fields.