EUSO-Balloon is a pathfinder mission for the Extreme Universe Space Observatory onboard the Japanese Experiment Module (JEM-EUSO). It was launched on the moonless night of the 25$^{th}$ of August 2014 from Timmins, Canada. The flight ended successfully after maintaining the target altitude of 38 km for five hours. One part of the mission was a 2.5 hour underflight using a helicopter equipped with three UV light sources (LED, xenon flasher and laser) to perform an inflight calibration and examine the detectors capability to measure tracks moving at the speed of light. We describe the helicopter laser system and details of the underflight as well as how the laser tracks were recorded and found in the data. These are the first recorded laser tracks measured from a fluorescence detector looking down on the atmosphere. Finally, we present a first reconstruction of the direction of the laser tracks relative to the detector.

We consider the problem of analyzing the structure of spectroscopic cubes using unsupervised machine learning techniques. We propose representing the target's signal as a homogeneous set of volumes through an iterative algorithm that separates the structured emission from the background while not overestimating the flux. Besides verifying some basic theoretical properties, the algorithm is designed to be tuned by domain experts, because its parameters have meaningful values in the astronomical context. Nevertheless, we propose a heuristic to automatically estimate the signal-to-noise ratio parameter of the algorithm directly from data. The resulting light-weighted set of samples ($\leq 1\%$ compared to the original data) offer several advantages. For instance, it is statistically correct and computationally inexpensive to apply well-established techniques of the pattern recognition and machine learning domains; such as clustering and dimensionality reduction algorithms. We use ALMA science verification data to validate our method, and present examples of the operations that can be performed by using the proposed representation. Even though this approach is focused on providing faster and better analysis tools for the end-user astronomer, it also opens the possibility of content-aware data discovery by applying our algorithm to big data.

Asteroseismology is witnessing a revolution thanks to high-precise asteroseismic space data (MOST, CoRoT, Kepler, BRITE), and their large ground-based follow-up programs. Those instruments have provided an unprecedented large amount of information, which allows us to scrutinize its statistical properties in the quest for hidden relations among pulsational and/or physical observables. This approach might be particularly useful for stars whose pulsation content is difficult to interpret. This is the case of intermediate-mass classical pulsating stars (i.e. gamma Dor, delta Scuti, hybrids) for which current theories do not properly predict the observed oscillation spectra. Here we establish a first step in finding such hidden relations from Data Mining techniques for these stars. We searched for those hidden relations in a sample of delta Scuti and hybrid stars observed by CoRoT and Kepler (74 and 153, respectively). No significant correlations between pairs of observables were found. However, two statistically significant correlations emerged from multivariable correlations in the observed seismic data, which describe the total number of observed frequencies and the largest one, respectively. Moreover, three different sets of stars were found to cluster according to their frequency density distribution. Such sets are in apparent agreement with the asteroseismic properties commonly accepted for A-F pulsating stars.

We study magnetic reconnection events in a turbulent plasma within the two-fluid theory. By identifying the diffusive regions, we measure the reconnection rates as function of the conductivity and current sheet thickness. We have found that the reconnection rate scales as the squared of the inverse of the current sheet's thickness and is independent of the aspect ratio of the diffusive region, in contrast to other analytical, e.g. the Sweet-Parker and Petscheck, and numerical models. Furthermore, while the reconnection rates are also proportional to the square inverse of the conductivity, the aspect ratios of the diffusive regions, which exhibit values in the range of $0.1-0.9$, are not correlated to the latter. Our findings suggest a new expression for the magnetic reconnection rate, which, after experimental verification, can provide a further understanding of the magnetic reconnection process.

Fluid dynamics in intrinsically curved geometries is encountered in many physical systems in nature, ranging from microscopic bio-membranes all the way up to general relativity at cosmological scales. Despite the diversity of applications, all of these systems share a common feature: the free motion of particles is affected by inertial forces originating from the curvature of the embedding space. Here we reveal a fundamental process underlying fluid dynamics in curved space: the free motion of fluids, in the complete absence of solid walls or obstacles, exhibits loss of energy due exclusively to the intrinsic curvature of space. We find that local sources of curvature generate viscous stresses as a result of the inertial forces. The curvature-induced viscous forces are shown to cause hitherto unnoticed and yet appreciable energy dissipation, which might play a significant role for a variety of physical systems involving fluid dynamics in curved spaces.

This paper describes the Third Public Data Release (DR3) of the Calar Alto Legacy Integral Field Area (CALIFA) survey. Science-grade quality data for 667 galaxies are made public, including the 200 galaxies of the Second Public Data Release (DR2). Data were obtained with the integral-field spectrograph PMAS/PPak mounted on the 3.5m telescope at the Calar Alto Observatory. Three different spectral setups are available, i) a low-resolution V500 setup covering the wavelength range 3749-7500 AA (4240-7140 AA unvignetted) with a spectral resolution of 6.0 AA (FWHM), for 646 galaxies, ii) a medium-resolution V1200 setup covering the wavelength range 3650-4840 AA (3650-4620 AA unvignetted) with a spectral resolution of 2.3 AA (FWHM), for 484 galaxies, and iii) the combination of the cubes from both setups (called COMBO), with a spectral resolution of 6.0 AA and a wavelength range between 3700-7500 AA (3700-7140 AA unvignetted), for 446 galaxies. The Main Sample, selected and observed according to the CALIFA survey strategy covers a redshift range between 0.005 and 0.03, spans the color-magnitude diagram and probes a wide range of stellar mass, ionization conditions, and morphological types. The Extension Sample covers several types of galaxies that are rare in the overall galaxy population and therefore not numerous or absent in the CALIFA Main Sample. All the cubes in the data release were processed using the latest pipeline, which includes improved versions of the calibration frames and an even further improved im- age reconstruction quality. In total, the third data release contains 1576 datacubes, including ~1.5 million independent spectra. It is available at http://califa.caha.es/DR3.

The search of extragalactic regions with conspicuous presence of Wolf-Rayet (WR) stars outside the Local Group is challenging task due to the difficulties in detecting their faint spectral features. In this exploratory work, we develop a methodology to perform an automated search of WR signatures through a pixel-by-pixel analysis of integral field spectroscopy (IFS) data belonging to the Calar Alto Legacy Integral Field Area survey, CALIFA. This technique allowed us to build the first catalogue of Wolf-Rayet rich regions with spatially-resolved information, allowing to study the properties of these complexes in a 2D context. The detection technique is based on the identification of the blue WR bump (around He II 4686 {\AA}, mainly associated to nitrogen-rich WR stars, WN) and the red WR bump (around C IV 5808 {\AA} and associated to carbon-rich WR stars, WC) using a pixel-by-pixel analysis. We identified 44 WR-rich regions with blue bumps distributed in 25 galaxies of a total of 558. The red WR bump was identified only in 5 of those regions. We found that the majority of the galaxies hosting WR populations in our sample are involved in some kind of interaction process. Half of the host galaxies share some properties with gamma-ray burst (GRB) hosts where WR stars, as potential candidates to being the progenitors of GRBs, are found. We also compared the WR properties derived from the CALIFA data with stellar population synthesis models, and confirm that simple star models are generally not able to reproduce the observations. We conclude that other effects, such as the binary star channel (which could extend the WR phase up to 10 Myr), fast rotation or other physical processes that causes the loss of observed Lyman continuum photons, are very likely affecting the derived WR properties, and hence should be considered when modelling the evolution of massive stars.

This paper aims at providing aperture corrections for emission lines in a sample of spiral galaxies from the Calar Alto Legacy Integral Field Area Survey (CALIFA) database. In particular, we explore the behavior of the log([OIII]5007/Hbeta)/([NII]6583/Halpha) (O3N2) and log[NII]6583/Halpha (N2) flux ratios since they are closely connected to different empirical calibrations of the oxygen abundances in star forming galaxies. We compute median growth curves of Halpha, Halpha/Hbeta, O3N2 and N2 up to 2.5R_50 and 1.5 disk R_eff. The growth curves simulate the effect of observing galaxies through apertures of varying radii. The median growth curve of the Halpha/Hbeta ratio monotonically decreases from the center towards larger radii, showing for small apertures a maximum value of ~10% larger than the integrated one. The median growth curve of N2 shows a similar behavior, decreasing from the center towards larger radii. No strong dependence is seen with the inclination, morphological type and stellar mass for these growth curves. Finally, the median growth curve of O3N2 increases monotonically with radius. However, at small radii it shows systematically higher values for galaxies of earlier morphological types and for high stellar mass galaxies. Applying our aperture corrections to a sample of galaxies from the SDSS survey at 0.02<=z<=0.3 shows that the average difference between fiber-based and aperture corrected oxygen abundances, for different galaxy stellar mass and redshift ranges, reaches typically to ~11%, depending on the abundance calibration used. This average difference is found to be systematically biased, though still within the typical uncertainties of oxygen abundances derived from empirical calibrations. Caution must be exercised when using observations of galaxies for small radii (e.g. below 0.5R_eff) given the high dispersion shown around the median growth curves.

We investigate Poiseuille channel flow through intrinsically curved media, equipped with localized metric perturbations. To this end, we study the flux of a fluid driven through the curved channel in dependence of the spatial deformation, characterized by the parameters of the metric perturbations (amplitude, range and density). We find that the flux depends only on a specific combination of parameters, which we identify as the average metric perturbation, and derive a universal flux law for the Poiseuille flow. For the purpose of this study, we have improved and validated our recently developed lattice Boltzmann model in curved space by considerably reducing discrete lattice effects.

In this paper, we develop a lattice Boltzmann model for relativistic magnetohydrodynamics (MHD). Even though the model is derived for resistive MHD, it is shown that it is numerically robust even in the high conductivity (ideal MHD) limit. In order to validate the numerical method, test simulations are carried out for both ideal and resistive limits, namely the propagation of Alfv\'en waves in the ideal MHD and the evolution of current sheets in the resistive regime, where very good agreement is observed comparing to the analytical results. Additionally, two-dimensional magnetic reconnection driven by Kelvin-Helmholtz instability is studied and the effects of different parameters on the reconnection rate are investigated. It is shown that the density ratio has negligible effect on the magnetic reconnection rate, while an increase in shear velocity decreases the reconnection rate. Additionally, it is found that the reconnection rate is proportional to $\sigma^{-\frac{1}{2}}$, $\sigma$ being the conductivity, which is in agreement with the scaling law of the Sweet-Parker model. Finally, the numerical model is used to study the magnetic reconnection in a stellar flare. Three-dimensional simulation suggests that the reconnection between the background and flux rope magnetic lines in a stellar flare can take place as a result of a shear velocity in the photosphere.

We have found that the relation between the flow through campylotic (generically curved) media, consisting of randomly located curvature perturbations, and the average Ricci scalar of the system exhibits two distinct functional expressions (hysteresis), depending on whether the typical spatial extent of the curvature perturbation lies above or below the critical value maximizing the overall Ricci curvature. Furthermore, the flow through such systems as a function of the number of curvature perturbations presents a sublinear behavior for large concentrations due to the interference between curvature perturbations that, consequently, produces a less curved space. For the purpose of this study, we have developed and validated a lattice kinetic model capable of describing fluid flow in arbitrarily curved manifolds, which allows to deal with highly complex spaces in a very compact and efficient way.

In the Bona-Masso formulation, Einstein equations are written as a set of flux conservative first order hyperbolic equations that resemble fluid dynamics equations. Based on this formulation, we construct a lattice Boltzmann model for Numerical Relativity. Our model is validated with well-established tests, showing good agreement with analytical solutions. Furthermore, we show that by increasing the relaxation time, we gain stability at the cost of losing accuracy, and by decreasing the lattice spacings while keeping a constant numerical diffusivity, the accuracy and stability of our simulations improves. Finally, in order to show the potential of our approach a linear scaling law for parallelisation with respect to number of CPU cores is demonstrated. Our model represents the first step in using lattice kinetic theory to solve gravitational problems.

We provide numerical evidence that the Richtmyer-Meshkov (RM) instability contributes to the cooling of a relativistic fluid. Due to the presence of jet particles traveling throughout the medium, shock waves are generated in the form of Mach cones. The interaction of multiple shock waves can trigger the RM instability, and we have found that this process leads to a down-cooling of the relativistic fluid. To confirm the cooling effect of the instability, shock tube Richtmyer-Meshkov instability simulations are performed. Additionally, in order to provide an experimental observable of the RM instability resulting from the Mach cone interaction, we measure the two particle correlation function and highlight the effects of the interaction. The simulations have been performed with an improved version of the relativistic lattice Boltzmann model, including general equations of state and external forces.

The Star Formation Rate (SFR) is one of the main parameters used to analyze the evolution of galaxies through time. The need for recovering the light reprocessed by dust commonly requires the use of low spatial resolution far-infrared data. Recombination-line luminosities provide an alternative, although uncertain dust-extinction corrections based on narrow-band imaging or long-slit spectroscopy have traditionally posed a limit to their applicability. Integral Field Spectroscopy (IFS) is clearly the way to overcome such limitation. We obtain integrated H{\alpha}, ultraviolet (UV) and infrared (IR)-based SFR measurements for 272 galaxies from the CALIFA survey at 0.005 < z < 0.03 using single-band and hybrid tracers. We provide updated calibrations, both global and split by properties (including stellar mass and morphological type), referred to H{\alpha}. The extinction-corrected H{\alpha} luminosity agrees with the updated hybrid SFR estimators based on either UV or H{\alpha} plus IR luminosity over the full range of SFRs (0.03-20 M$_{\odot}$ yr$^{-1}$). The coefficient that weights the amount of energy produced by newly-born stars that is reprocessed by dust on the hybrid tracers, a$_{IR}$, shows a large dispersion. However, it does not became increasingly small at high attenuations, as expected if significant highly-obscured H$\alpha$ emission would be missed. Lenticulars, early-type spirals and type-2 AGN host galaxies show smaller coefficients due to the contribution of optical photons and AGN to dust heating. In the Local Universe the H{\alpha} luminosity derived from IFS observations can be used to measure SFR, at least in statistically-significant, optically-selected galaxy samples. The analysis of the SFR calibrations by galaxies properties could be potentially used by other works to study the impact of different selection criteria in the SFR values derived.

Variations in the stellar initial mass function (IMF) have been invoked to explain the spectroscopic and dynamical properties of early-type galaxies. However, no observations have yet been able to disentangle the physical driver. We analyse here a sample of 24 early-type galaxies drawn from the CALIFA survey, deriving in a homogeneous way their stellar population and kinematic properties. We find that the local IMF is tightly related to the local metallicity, becoming more bottom-heavy towards metal-rich populations. Our result, combined with the galaxy mass-metallicity relation, naturally explains previous claims of a galaxy mass-IMF relation, derived from non-IFU spectra. If we assume that - within the star formation environment of early-type galaxies - metallicity is the main driver of IMF variations, a significant revision of the interpretation of galaxy evolution observables is necessary.

We present a new approach to find accurate solutions to the Poisson equation, as obtained from the steady-state limit of a diffusion equation with strong source terms. For this purpose, we start from Boltzmann's kinetic theory and investigate the influence of higher order terms on the resulting macroscopic equations. By performing an appropriate expansion of the equilibrium distribution, we provide a method to remove the unnecessary terms up to a desired order and show that it is possible to find, with high level of accuracy, the steady-state solution of the diffusion equation for sizeable Knudsen numbers. In order to test our kinetic approach, we discretise the Boltzmann equation and solve the Poisson equation, spending up to six order of magnitude less computational time for a given precision than standard lattice Boltzmann methods.

While studies of gas-phase metallicity gradients in disc galaxies are common, very little has been done in the acquisition of stellar abundance gradients in the same regions. We present here a comparative study of the stellar metallicity and age distributions in a sample of 62 nearly face-on, spiral galaxies with and without bars, using data from the CALIFA survey. We measure the slopes of the gradients and study their relation with other properties of the galaxies. We find that the mean stellar age and metallicity gradients in the disc are shallow and negative. Furthermore, when normalized to the effective radius of the disc, the slope of the stellar population gradients does not correlate with the mass or with the morphological type of the galaxies. Contrary to this, the values of both age and metallicity at $\sim$2.5 scale-lengths correlate with the central velocity dispersion in a similar manner to the central values of the bulges, although bulges show, on average, older ages and higher metallicities than the discs. One of the goals of the present paper is to test the theoretical prediction that non-linear coupling between the bar and the spiral arms is an efficient mechanism for producing radial migrations across significant distances within discs. The process of radial migration should flatten the stellar metallicity gradient with time and, therefore, we would expect flatter stellar metallicity gradients in barred galaxies. However, we do not find any difference in the metallicity or age gradients in galaxies with without bars. We discuss possible scenarios that can lead to this absence of difference.

Solar flares stem from the reconnection of twisted magnetic field lines in the solar photosphere. The energy and waiting time distributions of these events follow complex patterns that have been carefully considered in the past and that bear some resemblance with earthquakes and stockmarkets. Here we explore in detail the tangling motion of interacting flux tubes anchored in the plasma and the energy ejections resulting when they recombine. The mechanism for energy accumulation and release in the flow is reminiscent of self-organized criticality. From this model we suggest the origin for two important and widely studied properties of solar flare statistics, including the time-energy correlations. We first propose that the scale-free energy distribution of solar flares is largely due to the twist exerted by the vorticity of the turbulent photosphere. Second, the long-range temporal and time-energy correlations appear to arise from the tube-tube interactions. The agreement with satellite measurements is encouraging.

H ii regions in galaxies are the sites of star formation and thus particular places to understand the build-up of stellar mass in the universe. The line ratios of this ionized gas are frequently used to characterize the ionization conditions. We use the Hii regions catalogue from the CALIFA survey (~5000 H ii regions), to explore their distribution across the classical [OIII]/Hbeta vs. [NII]/Halpha diagnostic diagram, and how it depends on the oxygen abundance, ionization parameter, electron density, and dust attenuation. We compared the line ratios with predictions from photoionization models. Finally, we explore the dependences on the properties of the host galaxies, the location within those galaxies and the properties of the underlying stellar population. We found that the location within the BPT diagrams is not totally predicted by photoionization models. Indeed, it depends on the properties of the host galaxies, their galactocentric distances and the properties of the underlying stellar population. These results indicate that although H ii regions are short lived events, they are affected by the total underlying stellar population. One may say that H ii regions keep a memory of the stellar evolution and chemical enrichment that have left an imprint on the both the ionizing stellar population and the ionized gas

This paper describes the Second Public Data Release (DR2) of the Calar Alto Legacy Integral Field Area (CALIFA) survey. The data for 200 objects are made public, including the 100 galaxies of the First Public Data Release (DR1). Data were obtained with the integral-field spectrograph PMAS/PPak mounted on the 3.5 m telescope at the Calar Alto observatory. Two different spectral setups are available for each galaxy, (i) a low-resolution V500 setup covering the wavelength range 3745-7500 \AA\ with a spectral resolution of 6.0 \AA\ (FWHM), and (ii) a medium-resolution V1200 setup covering the wavelength range 3650-4840 \AA\ with a spectral resolution of 2.3 \AA\ (FWHM). The sample covers a redshift range between 0.005 and 0.03, with a wide range of properties in the Color-Magnitude diagram, stellar mass, ionization conditions, and morphological types. All released cubes were reduced with the latest pipeline, including improved spectrophotometric calibration, spatial registration and spatial resolution. The spectrophotometric calibration is better than 6% and the median spatial resolution is 2.5". Altogether the second data release contains over 1.5 million spectra. It is available at http://califa.caha.es/DR2.

[Abridged] We use optical IFS of nearby SN host galaxies provided by the CALIFA Survey with the goal of finding correlations in the environmental parameters at the location of different SN types. We recover the sequence in association of different SN types to the star-forming regions by using several indicators of the ongoing and recent SF related to both the ionized gas and the stellar populations. While the total ongoing SF is on average the same for the three SN types, SNe Ibc/IIb tend to happen closer to star-forming regions and occur in higher SF density locations compared to SNe II and SNe~Ia, the latter showing the weakest correlation. SNe~Ia host galaxies have on average masses that are $\sim$0.3-0.8~dex higher than CC SNe hosts due to a larger fraction of old stellar populations in the SNe~Ia hosts. Using the recent SN~Ia delay-time distribution and the SFHs of the galaxies, we show that the SN~Ia hosts in our sample should presently produce a factor 2 more SNe~Ia than the CC~SN hosts. Since both types are in hosts with similar SF rate and hence similar CC~SN rate, this can explain the mass difference between the SN~Ia and CC~SN hosts, and reinforce the finding that at least part of SNe~Ia should originate from very old progenitors. Comparing the mean SFH of the eight least massive galaxies to that of the massive SF SN~Ia hosts we find that the low-mass galaxies formed their stars over more extended time than the massive SN~Ia hosts. We estimate that the low-mass galaxies should produce by a factor of 10 less SNe~Ia, and a factor of 3 less CC~SNe than the high-mass group. Therefore the ratio between the number of CC~SNe and SNe~Ia is expected to increase with decreasing the galaxy mass. CC~SNe tend to explode at positions with younger stellar populations than the galaxy average, but the galaxy properties at SNe~Ia locations are one average the same as the global ones.

We present the largest and most homogeneous catalog of HII regions and associations compiled so far. The catalog comprises more than 7000 ionized regions, extracted from 306 galaxies observed by the CALIFA survey. We describe the procedures used to detect, select, and analyse the spectroscopic properties of these ionized regions. In the current study we focus on the characterization of the radial gradient of the oxygen abundance in the ionized gas, based on the study of the deprojected distribution of HII regions. We found that all galaxies without clear evidence of an interaction present a common gradient in the oxygen abundance, with a characteristic slope of alpha = -0.1 dex/re between 0.3 and 2 disk effective radii, and a scatter compatible with random fluctuations around this value, when the gradient is normalized to the disk effective radius. The slope is independent of morphology, incidence of bars, absolute magnitude or mass. Only those galaxies with evidence of interactions and/or clear merging systems present a significant shallower gradient, consistent with previous results. The majority of the 94 galaxies with H ii regions detected beyond 2 disk effective radii present a flattening in the oxygen abundance. The flattening is statistically significant. We cannot provide with a conclusive answer regarding the origin of this flattening. However, our results indicate that its origin is most probably related to the secular evolution of galaxies. Finally, we find a drop/truncation of the oxygen abundance in the inner regions for 26 of the galaxies. All of them are non-interacting, mostly unbarred, Sb/Sbc galaxies. This feature is associated with a central star-forming ring, which suggests that both features are produced by radial gas flows induced by resonance processes.

Theoretical and numerical analysis of the relativistic effects on the Richtmyer-Meshkov (RM) instability reveals new and potentially very useful effects. We find that, in contrast with the non- relativistic case, the growth rate of the RM instability depends strongly on the equation of state of the fluid, opening up the possibility to infer equations of state from experimental observations of the RM instability. As opposed to the non-relativistic case, we also discover that, above a critical value of the fluid velocity, the growth rate of the instability counter-intuitively decreases due to the Lorentz's factor, and vanishes in the ultrarelativistic limit, as the speed of the particles approaches the speed of light. Both effects might prove very useful for leading-edge applications, such as the study of the equation of state of quark-gluon matter, and the design of fast ignition inertial confinement fusion (ICF) schemes. We perform a linear stability analysis to characterize the instability, for an arbitrary equation of state, and implement numerical simulations to study the instability in the non-linear regime, using the equation of state of an ideal gas. Furthermore, based on the numerical results, we propose a general expression that characterizes the long term evolution of the instability.

We develop a relativistic lattice Boltzmann model capable of describing relativistic fluid dynamics at ultra-high velocities, with Lorentz factors up to $\gamma \sim 10$. To this purpose, we first build a new lattice kinetic scheme by expanding the Maxwell-J\"uttner distribution function in an orthogonal basis of polynomials and applying an appropriate quadrature, providing the discrete versions of the relativistic Boltzmann equation and the equilibrium distribution. To achieve ultra-high velocities, we include a flux limiter scheme, and introduce the bulk viscosity by a suitable extension of the discrete relativistic Boltzmann equation. The model is validated by performing simulations of shock waves in viscous quark-gluon plasmas and comparing with existing models, finding very good agreement. To the best of our knowledge, we for the first time successfully simulate viscous shock waves in the highly relativistic regime. Moreover, we show that our model can also be used for near-inviscid flows even at very high velocities. Finally, as an astrophysical application, we simulate a relativistic shock wave, generated by, say, a supernova explosion, colliding with a massive interstellar cloud, e.g. molecular gas.

A detailed derivation of the Lattice Boltzmann (LB) scheme for relativistic fluids recently proposed in Ref. [1], is presented. The method is numerically validated and applied to the case of two quite different relativistic fluid dynamic problems, namely shock-wave propagation in quark-gluon plasmas and the impact of a supernova blast-wave on massive interstellar clouds. Close to second order convergence with the grid resolution, as well as linear dependence of computational time on the number of grid points and time-steps, are reported.