Over the last few years it became clear that turbulent magnetic reconnection and magnetized turbulence are inseparable. It was not only shown that reconnection is responsible for violating the frozen-in condition in turbulence, but also that stochastic reconnection in 3D generates turbulence by itself. The actual mechanism responsible for this driving is still unknown. Processes such tearing mode or Kelvin-Helmholtz, among other plasma instabilities, could generate turbulence from irregular current sheets. We address the nature of driving mechanism for this process and consider a relative role of tearing and Kelvin-Helmholtz instabilities for the process of turbulence generation. In particular, we analyze the conditions for development of these two instabilities within three-dimensional reconnection regions. We show that both instabilities can excite turbulence fluctuations in reconnection regions. However, tearing mode has relatively slow growth rate, and at later times it becomes partially suppressed by transverse to the current sheet component of magnetic field, generated during the growth of turbulent fluctuations. On the contrary, the Kelvin-Helmholtz instability establishes quickly in the outflow region, and at later times it dominates the turbulence generation comparing to the contribution from tearing mode. Our results demonstrate that the tearing instability is subdominant compared to the the Kelvin-Helmholtz instability in terms of generation of turbulence in the 3D reconnection layers and therefore the self-driven reconnection is turbulent reconnection with tearing instability being important only at the initial stage of the reconnection.

The study of the acceleration of particles is an essential element of research in the heliospheric science. Here, we discuss the predisposition to the particle acceleration around coronal mass ejections (CMEs)-driven shocks with corrugated wave-like features. We adopt these attributes on shocks formed from disturbances due to the bimodal solar wind, CME deflection, irregular CME expansion, and the ubiquitous fluctuations in the solar corona. In order to understand the role of a wavy shock in particle acceleration, we define three initial smooth shock morphologies each one associated with a fast CME. Using polar Gaussian profiles we model these shocks in the low corona. We establish the corrugated appearance on smooth shock by using combinations of wave-like functions that represent the disturbances from medium and CME piston. For both shock types, smooth and corrugated, we calculate the shock normal angles between the shock normal and the radial upstream coronal magnetic field in order to classify the quasi-parallel and quasi-perpendicular regions. We consider that corrugated shocks are predisposed to different process of particle acceleration due to irregular distributions of shock normal angles around of the shock. We suggest that disturbances due to CME irregular expansion may be a decisive factor in origin of particle acceleration. Finally, we regard that accepting these features on shocks may be the start point for investigating some questions in the sheath and shock, like downstream-jets, instabilities, shock thermalization, shock stability, and injection particle process.

We investigate the general stability of 1D spherically symmetric ionized Bondi accretion onto a massive object in the specific context of accretion onto a young stellar object. We first derive a new analytic expression for a steady state two temperature solution that predicts the existence of compact and hypercompact HII regions. We then show that this solution is only marginally stable if ionization is treated self-consistently. This leads to a recurring collapse of the HII region over time. We derive a semi-analytic model to explain this instability, and test it using spatially converged 1D radiation hydrodynamical simulations. We discuss the implications of the 1D instability on 3D radiation hydrodynamics simulations of supersonic accreting flows.

This paper investigates the gravitational trapping of HII regions predicted by steady-state analysis using radiation hydrodynamical simulations. We present idealised spherically symmetric radiation hydrodynamical simulations of the early evolution of HII regions including the gravity of the central source. As with analytic steady state solutions of spherically symmetric ionised Bondi accretion flows, we find gravitationally trapped HII regions with accretion through the ionisation front onto the source. We found that, for a constant ionising luminosity, fluctuations in the ionisation front are unstable. This instability only occurs in this spherically symmetric accretion geometry. In the context of massive star formation, the ionising luminosity increases with time as the source accretes mass. The maximum radius of the recurring HII region increases on the accretion timescale until it reaches the sonic radius, where the infall velocity equals the sound speed of the ionised gas, after which it enters a pressure-driven expansion phase. This expansion prevents accretion of gas through the ionisation front, the accretion rate onto the star decreases to zero, and it stops growing from accretion. Because of the time required for any significant change in stellar mass and luminosity through accretion our simulations keep both mass and luminosity constant and follow the evolution from trapped to expanding in a piecewise manner. Implications of this evolution of HII regions include a continuation of accretion of material onto forming stars for a period after the star starts to emit ionising radiation, and an extension of the lifetime of ultracompact HII regions.

We present results from a non-cosmological, three-dimensional hydrodynamical simulation of the gas in the dwarf spheroidal galaxy Ursa Minor. Assuming an initial baryonic-to-dark-matter ratio derived from the cosmic microwave background radiation, we evolved the galactic gas distribution over 3 Gyr, taking into account the effects of the types Ia and II supernovae. For the first time, we used in our simulation the instantaneous supernovae rates derived from a chemical evolution model applied to spectroscopic observational data of Ursa Minor. We show that the amount of gas that is lost in this process is variable with time and radius, being the highest rates observed during the initial 600 Myr in our simulation. Our results indicate that types Ia and II supernovae must be essential drivers of the gas loss in Ursa Minor galaxy (and probably in other similar dwarf galaxies), but it is ultimately the combination of galactic winds powered by these supernovae and environmental effects (e.g., ram-pressure stripping) that results in the complete removal of the gas content.

Magnetic reconnection is a process that changes magnetic field topology in highly conducting fluids. Within the standard Sweet-Parker model, this process would be too slow to explain observations (e.g. solar flares). In reality, the process must be ubiquitous as astrophysical fluids are magnetized and motions of fluid elements necessarily entail crossing of magnetic frozen-in field lines and magnetic reconnection. In the presence of turbulence, the reconnection is independent of microscopic plasma properties, and may be much faster than previously thought, as proposed in Lazarian & Vishniac (1999) and tested in Kowal et al. (2009, 2012). However, the considered turbulence in the Lazarian-Vishniac model was imposed externally. In this work we consider reconnection-driven magnetized turbulence in realistic three-dimensional geometry initiated by stochastic noise. We demonstrate through numerical simulations that the stochastic reconnection is able to self-generate turbulence through interactions between its outflows. We analyze the statistical properties of velocity fluctuations using power spectra and anisotropy scaling, which demonstrates that the reconnection produces Kolmogorov-like turbulence, compatible with Goldreich-Sridhar (1995) model. Anisotropy statistics are, however, strongly affected by the dynamics of reconnection outflows. Once the broad turbulent region is formed, the typical anisotropy scaling $l_\parallel \propto l_\perp^{2/3}$ is formed, especially for high resolution models, were the broader range of scales is available. The decay of reconnection outflows to turbulent-like fluctuations, characterized by different anisotropy scalings, strongly depends on $\beta$ plasma. Moreover, the estimated reconnection rates are weakly dependent on the resolution, suggesting that no external processes are required to make reconnection fast.

Statistical analysis of Faraday Rotation Measure (RM) maps of the intracluster medium (ICM) of galaxy clusters provides a unique tool to evaluate some spatial features of the magnetic fields there. Its combination with numerical simulations of magnetohydrodynamic (MHD) turbulence allows the diagnosis of the ICM turbulence. Being the ICM plasma weakly collisional, the thermal velocity distribution of the particles naturally develops anisotropies as a consequence of the large scale motions and the conservation of the magnetic moment of the charged particles. A previous study (Paper I) analyzed the impact of large scale thermal anisotropy on the statistics of RM maps synthesized from simulations of turbulence; these simulations employed a collisionless MHD model which considered a tensor pressure with uniform anisotropy. In the present work, we extend that analysis to a collisionless MHD model in which the thermal anisotropy develops according to the conservation of the magnetic moment of the thermal particles. We also consider the effect of anisotropy relaxation caused by the micro-scale mirror and firehose instabilities. We show that if the relaxation rate is fast enough to keep the anisotropy limited by the threshold values of the instabilities, the dispersion and power spectrum of the RM maps are indistinguishable from those obtained from collisional MHD. Otherwise, there is a reduction in the dispersion and steepening of the power spectrum of the RM maps (compared to the collisional case). Considering the first scenario, the use of collisional MHD simulations for modeling the RM statistics in the ICM becomes better justified.

Observations of the intracluster medium (ICM) in galaxy clusters suggest for the presence of turbulence and the magnetic fields existence has been proved through observations of Faraday Rotation and synchrotron emission. The ICM is also known to be filled by a rarefied weakly collisional plasma. In this work we study the possible signatures left on Faraday Rotation maps by collisionless instabilities. For this purpose we use a numerical approach to investigate the dynamics of the turbulence in collisionless plasmas based on an magnetohydrodynamical (MHD) formalism taking into account different levels of pressure anisotropy. We consider models covering the sub/super-Alfv\'enic and trans/supersonic regimes, one of them representing the fiducial conditions corresponding to the ICM. From the simulated models we compute Faraday Rotation maps and analyze several statistical indicators in order to characterize the magnetic field structure and compare the results obtained with the collisionless model to those obtained using standard collisional MHD framework. We find that important imprints of the pressure anisotropy prevails in the magnetic field and also manifest in the associated Faraday Rotation maps which evidence smaller correlation lengths in the collisionless MHD case. These points are remarkably noticeable for the case mimicking the conditions prevailing in ICM. Nevertheless, in this study we have neglected the decrease of pressure anisotropy due to the feedback of the instabilities that naturally arise in collisionless plasmas at small scales. This decrease may not affect the statistical imprint differences described above, but should be examined elsewhere.

Within ten nearby (d < 450 pc) Gould Belt molecular clouds we evaluate statistically the relative orientation between the magnetic field projected on the plane of sky, inferred from the polarized thermal emission of Galactic dust observed by Planck at 353 GHz, and the gas column density structures, quantified by the gradient of the column density, $N_H$. The selected regions, covering several degrees in size, are analyzed at an effective angular resolution of 10' FWHM, thus sampling physical scales from 0.4 to 40 pc in the nearest cloud. The column densities in the selected regions range from $N_H \approx 10^{21}$ to $10^{23}$ cm$^{-2}$, and hence they correspond to the bulk of the molecular clouds. The relative orientation is evaluated pixel by pixel and analyzed in bins of column density using the novel statistical tool called "Histogram of Relative Orientations". Throughout this study, we assume that the polarized emission observed by Planck at 353 GHz is representative of the projected morphology of the magnetic field in each region, i.e., we assume a constant dust grain alignment efficiency, independent of the local environment. Within most clouds we find that the relative orientation changes progressively with increasing $N_H$, from preferentially parallel or having no preferred orientation to preferentially perpendicular. In simulations of magnetohydrodynamic turbulence in molecular clouds this trend in relative orientation is a signature of Alfv\'enic or sub-Alfv\'enic turbulence, implying that the magnetic field is significant for the gas dynamics at the scales probed by Planck. We compare the deduced magnetic field strength with estimates we obtain from other methods and discuss the implications of the Planck observations for the general picture of molecular cloud formation and evolution.

We report observations of $\eta$ Carinae obtained with ALMA in the continuum of 100, 230, 280 and 660 GHz in 2012 November, with a resolution that varied from 2.88" to 0.45" for the lower and higher frequencies respectively. The source is not resolved, even at the highest frequency; its spectrum is characteristic of thermal bremsstrahlung of a compact source, but different from the spectrum of optically thin wind. The recombination lines H42$\alpha$, He42$\alpha$, H40$\alpha$, He40$\alpha$, H50$\beta$, H28$\alpha$, He28$\alpha$, H21$\alpha$ and He21$\alpha$ were also detected and their intensities reveal non local thermodynamic equilibrium (NLTE) effects. We found that the line profiles could only be fit by an expanding shell of dense and ionized gas, which produces a slow shock in the surroundings of $\eta$ Carinae. Combined with fittings to the continuum, we were able to constrain the shell size, radius, density, temperature and velocity. The detection of the He recombination lines is compatible with the high temperature gas and requires a high energy ionizing photon flux, which must be provided by the companion star. The mass loss rate and wind velocity, necessary to explain the formation of the shell, are compatible with a LBV eruption. The position, velocity and physical parameters of the shell coincide with those of the Weigelt blobs. The dynamics found for the expanding shell corresponds to matter ejected by $\eta$ Carinae in 1941, in an event similar to that which formed the Little Homunculus ; for that reason we called the new ejecta the "Baby Homunculus".

The amplification and maintenance of the observed magnetic fields in the ICM are usually attributed to the turbulent dynamo action. This is generally derived employing a collisional MHD model. However, in the ICM the ion mean free path between collisions is of the order of the dynamical scales, thus requiring a collisionless MHD description. Unlike collisional MHD simulations, our study uses an anisotropic plasma pressure with respect to the direction of the local magnetic field, which brings the plasma within a parameter space where collisionless instabilities should take place. Within the adopted model these instabilities are contained at bay through the relaxation term of the pressure anisotropy which simulates the feedback of the mirror and firehose instabilities. This relaxation acts to get the plasma distribution function consistent with the empirical studies of collisionless plasmas. Our 3D numerical simulations of forced transonic turbulence motivated by modeling of the turbulent ICM are performed for different initial values of the magnetic field intensity, and different relaxation rates of the pressure anisotropy. We found that in the high beta plasma regime corresponding to the ICM conditions, a fast anisotropy relaxation rate gives results which are similar to the collisional-MHD model as far as the statistical properties of the turbulence are concerned. Also, the amplification of seed magnetic fields due to the turbulent dynamo action is similar to the collisional-MHD model. Our simulations that do not employ the anisotropy relaxation deviate significantly from the collisional-MHD results, and show more power at the small-scale fluctuations of both density and velocity representing the results of the instabilities. For these simulations the large scale fluctuations in the magnetic field are mostly suppressed and the turbulent dynamo fails in amplifying seed magnetic fields.

The role of turbulence in astrophysical environments and its interplay with magnetic fields is still highly debated. In this lecture, we will discuss this issue in the framework of dynamo processes. We will first present a very brief summary of turbulent dynamo theories, then will focus on small scale turbulent dynamos and their particular relevance on the origin and maintenance of magnetic fields in the intra-cluster media (ICM) of galaxies. In these environments, the very low density of the flow requires a collisionless-MHD treatment. We will show the implications of this approach in the turbulent amplification of the magnetic fields in these environments. To finalize, we will also briefly address the connection between MHD turbulence and fast magnetic reconnection and its possible implications in the diffusion of magnetic flux in the dynamo process.

Theoretical $\Lambda$CDM cosmological models predict a much larger number of low mass dark matter haloes than has been observed in the Local Group of galaxies. One possible explanation is the increased difficulty of detecting these haloes if most of the visible matter is lost at early evolutionary phases through galactic winds. In this work we study the current models of triggering galactic winds in dwarf spheroidal galaxies (dSph) from supernovae, and study, based on 3D hydrodynamic numerical simulations, the correlation of the mass loss rates and important physical parameters as the dark matter halo mass and its radial profile, and the star formation rate. We find that the existence of winds is ubiquitous, independent on the gravitational potential. Our simulations revealed that the Rayleigh-Taylor Instability (RTI) may play a major role on pushing matter out of these systems, even for very massive haloes. The instability is responsible for 5 - 40% of the mass loss during the early evolution of the galaxy, being less relevant at $t > 200$Myrs. There is no significant difference in the mass loss rates obtained for the different dark matter profiles studied (NFW and logarithmic). We have also found a correlation between the mass loss rate and both the halo mass and the rate of supernovae, as already reported in previous works. Besides, the epoch in which most of the baryon galactic matter is removed from the galaxy varies depending on the SN rate and gravitational potential. The later, combined to the importance of the RTI in each model, may change our understanding about the chemical evolution of dwarf galaxies, as well as in the heavy element contamination of the intergalactic medium at high redshifts.

We report a study of the stellar content of the Near-infrared cluster [DBS2003]\,157 embedded in the extended H\,{\sc ii} region GAL\,331.31-00.34, which is associated with the IRAS source 16085-5138. $JHK$ photometry was carried out in order to identify potential ionizing candidates, and the follow-up NIR spectroscopy allowed the spectral classification of some sources, including two O-type stars. A combination of NIR photometry and spectroscopy data was used to obtain the distance of these two stars, with the method of spectroscopic parallax: IRS\,298 (O6\,{\sc V}, $3.35 \pm 0.61$\,kpc) and IRS\,339 (O9\,{\sc V}, $3.24 \pm 0.56$\,kpc). Adopting the average distance of $3.29 \pm 0.58$\,kpc and comparing the Lyman continuum luminosity of these stars with that required to account for the radio continuum flux of the H\,{\sc ii} region, we conclude that these two stars are the ionizing sources of GAL\,331.31-00.34. Young stellar objects (YSOs) were searched by using our NIR photometry and MIR data from the GLIMPSE survey. The analysis of NIR and MIR colour-colour diagrams resulted in 47 YSO candidates. The GLIMPSE counterpart of IRAS\,16085-5138, which presents IRAS colour indices compatible with an ultra-compact H\,{\sc ii} region, has been identified. The analysis of its spectral energy distribution between 2 and $100\,\mu$m revealed that this source shows a spectral index $\alpha = 3.6$ between 2 and $25\,\mu$m, which is typical of a YSO immersed in a protostellar envelope. Lower limits to the bolometric luminosity and the mass of the embedded protostar have been estimated as $L=7.7\times10^3L_{\sun}$ and $M = 10\,M_{\sun}$, respectively, which corresponds to a B0--B1\,{\sc V} ZAMS star.

In the past few decades detailed observations of radio and X-rays emission from massive binary systems revealed a whole new physics present in such systems. Both thermal and non-thermal components of this emission indicate that most of the radiation at these bands originates in shocks. OB and WR stars present supersonic and massive winds that, when colliding, emit largely due to the free-free radiation. The non-thermal radio and X-ray emissions are due to synchrotron and inverse compton processes, respectively. In this case, magnetic fields are expected to play an important role on the emission distribution. In the past few years the modeling of the free-free and synchrotron emissions from massive binary systems have been based on purely hydrodynamical simulations, and ad hoc assumptions regarding the distribution of magnetic energy and the field geometry. In this work we provide the first full MHD numerical simulations of wind-wind collision in massive binary systems. We study the free-free emission characterizing its dependence on the stellar and orbital parameters. We also study self-consistently the evolution of the magnetic field at the shock region, obtaining also the synchrotron energy distribution integrated along different lines of sight. We show that the magnetic field in the shocks is larger than that obtained when the proportionality between $B$ and the plasma density is assumed. Also, we show that the role of the synchrotron emission relative to the total radio emission has been underestimated.

Continuing our series of papers on the 3-D structure and accurate distances of Planetary Nebulae (PNe), we present here the results obtained for the planetary nebula NGC\,40. Using data from different sources and wavelengths, we construct 3-D photoionization models and derive the physical quantitities of the ionizing source and nebular gas. The procedure, discussed in detail in the previous papers, consists of the use of 3-D photoionization codes constrained by observational data to derive the three-dimensional nebular structure, physical and chemical characteristics and ionizing star parameters of the objects by simultaneously fitting the integrated line intensities, the density map, the temperature map, and the observed morphologies in different emission lines. For this particular case we combined hydrodynamical simulations with the photoionization scheme in order to obtain self-consistent distributions of density and velocity of the nebular material. Combining the velocity field with the emission line cubes we also obtained the synthetic position-velocity plots that are compared to the observations. Finally, using theoretical evolutionary tracks of intermediate and low mass stars, we derive the mass and age of the central star of NGC\,40 as $(0.567 \pm 0.06)$M$_{\odot}$ and $(5810 \pm 600)$yrs, respectively. The distance obtained from the fitting procedure was $(1150 \pm 120)$pc.

The standard magnetohydrodynamic (MHD) description of the plasma in the hot, magnetized gas of the intra-cluster (ICM) medium is not adequate because it is weakly collisional. In such collisionless magnetized gas, the microscopic velocity distribution of the particles is not isotropic, giving rise to kinetic effects on the dynamical scales. These kinetic effects could be important in understanding the turbulence, as so as the amplification and maintenance of the magnetic fields in the ICM. It is possible to formulate fluid models for collisonless or weakly collisional gas by introducing modifications in the MHD equations. These models are often referred as kinetic MHD (KMHD). Using a KMHD model based on the CGL-closure, which allows the adiabatic evolution of the two components of the pressure tensor (the parallel and perpendicular components with respect to the local magnetic field), we performed 3D numerical simulations of forced turbulence in order to study the amplification of an initially weak seed magnetic field. We found that the growth rate of the magnetic energy is comparable to that of the ordinary MHD turbulent dynamo, but the magnetic energy saturates in a level smaller than of the MHD case. We also found that a necessary condition for the dynamo works is to impose limits to the anisotropy of the pressure.

The transport of magnetic flux to outside of collapsing molecular clouds is a required step to allow the formation of stars. Although ambipolar diffusion is often regarded as a key mechanism for that, it has been recently argued that it may not be efficient enough. In this review, we discuss the role that MHD turbulence plays in the transport of magnetic flux in star forming flows. In particular, based on recent advances in the theory of fast magnetic reconnection in turbulent flows, we will show results of three-dimensional numerical simulations that indicate that the diffusion of magnetic field induced by turbulent reconnection can be a very efficient mechanism, especially in the early stages of cloud collapse and star formation. To conclude, we will also briefly discuss the turbulence-star formation connection and feedback in different astrophysical environments: from galactic to cluster of galaxy scales.