-
We present a multi-wavelength spectral and temporal analysis of the blazar OJ
287 during its recent activity between July 2016 -- July 2017. This period is
characterized by intense variability from optical to X-rays and frequent daily
detections at {\it Fermi}-LAT $\gamma$-ray energies but with no sign of
variations. The intense X-ray variation is accompanied by a change in spectral
index, implying an appearance of a new component affecting only the
X-ray-to-optical part of the spectral energy distribution (SED), without a
significant contribution to the \emph{Fermi}-LAT band. Cross-correlation
analyses during four continuous segments of the light curves, divided so as to
have similar base level fluxes before and after the rise in the
\emph{Swift}-UVOT bands, show simultaneous variations in optical-UV bands. The
X-ray and UV/optical are also simultaneous during the middle two segments. In
the first segment, the results suggest the X-rays may be lagging, while in the
last segment, the X-rays lead the optical/UV by $\sim$ 5-6 days and also show
some systematic trends with variations appearing first the higher energies
followed by lower energies. The LAT SED for the first half is similar to the
previous quiescent state SED while it becomes harder in the latter half,
coincident with the detections at VHEs by VERITAS with an extrapolated VHE SED
smoothly matching the LAT data. The overall broadband SEDs see to be a
combination of a typical OJ 287 SED and a high-energy peaked blazar (HBL) SED,
and can be explained through a two-zone leptonic model, with the second zone,
located at parsec scales beyond the BLR region, responsible for the HBL like
spectrum. The observed variations are accompanied by a significant change in
polarization degree and angle, with both systematic and chaotic changes
suggesting a dynamic role of the magnetic field.
-
We present a multi-wavelength spectral and temporal analysis of the blazar OJ
287 during its recent activity between December 2015 -- May 2016, showing
strong variability in the near-infrared (NIR) to X-ray energies with detection
at $\gamma$-ray energies as well. Most of the optical flux variations exhibit
strong changes in polarization angle and degree. All the inter-band time lags
are consistent with simultaneous emissions. Interestingly, on days with
excellent data coverage in the NIR--UV bands, the spectral energy distributions
(SEDs) show signatures of bumps in the visible--UV bands, never seen before in
this source. The optical bump can be explained as accretion-disk emission
associated with the primary black hole of mass $\sim \rm 1.8 \times10^{10}
M_{\odot}$ while the little bump feature in the optical-UV appears consistent
with line emission. Further, the broadband SEDs extracted during the first
flare and during a quiescent period during this span show very different
$\gamma$-ray spectra compared to previously observed flare or quiescent
spectra. The probable thermal bump in the visible seems to have been clearly
present since May 2013, as found by examining all available NIR-optical
observations, and favors the binary super-massive black hole model. The
simultaneous multi-wavelength variability and relatively weak $\gamma$-ray
emission that shows a shift in the SED peak is consistent with $\gamma$-ray
emission originating from inverse Compton scattering of photons from the line
emission that apparently contributes to the little blue bump.
-
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.
-
List of contributions from the Cherenkov Telescope Array (CTA) Consortium
presented at the 6th International Symposium on High-Energy Gamma-Ray Astronomy
(Gamma 2016), July 11-15, 2016, in Heidelberg, Germany.
-
Fast magnetic reconnection may occur in different astrophysical sources,
producing flare-like emission and particle acceleration. Currently, this
process is being studied as an efficient mechanism to accelerate particles via
a first-order Fermi process. In this work we analyse the acceleration rate and
the energy distribution of test particles injected in three-dimensional
magnetohydrodynamical (MHD) domains with large-scale current sheets where
reconnection is made fast by the presence of turbulence. We study the
dependence of the particle acceleration time with the relevant parameters of
the embedded turbulence, i.e., the Alfv\'en speed $V_{\rm A}$, the injection
power $P_{\rm inj}$ and scale $k_{\rm inj}$ ($k_{\rm inj} = 1/l_{\rm inj}$). We
find that the acceleration time follows a power-law dependence with the
particle kinetic energy: $t_{acc} \propto E^{\alpha}$, with $0.2 < \alpha <
0.6$ for a vast range of values of $c / V_{\rm A} \sim 20 - 1000$. The
acceleration time decreases with the Alfv\'en speed (and therefore with the
reconnection velocity) as expected, having an approximate dependence $t_{acc}
\propto (V_{\rm A} / c)^{-\kappa}$, with $\kappa \sim 2.1- 2.4$ for particles
reaching kinetic energies between $1 - 100 \, m_p c^2$, respectively.
Furthermore, we find that the acceleration time is only weakly dependent on the
$P_{\rm inj}$ and $l_{\rm inj}$ parameters of the turbulence. The particle
spectrum develops a high-energy tail which can be fitted by a hard power-law
already in the early times of the acceleration, in consistency with the results
of kinetic studies of particle acceleration by magnetic reconnection in
collisionless plasmas.
-
The phenomenon of solar "torsional oscillations" (TO) represents migratory
zonal flows associated with the solar cycle. These flows are observed on the
solar surface and, according to helioseismology, extend through the convection
zone. We study the origin of the TO using results from a global MHD simulation
of the solar interior that reproduces several of the observed characteristics
of the mean-flows and magnetic fields. Our results indicate that the magnetic
tension (MT) in the tachocline region is a key factor for the periodic changes
in the angular momentum transport that causes the TO. The torque induced by the
MT at the base of the convection zone is positive at the poles and negative at
the equator. A rising MT torque at higher latitudes causes the poles to
speed-up, whereas a declining negative MT torque at the lower latitudes causes
the equator to slow-down. These changes in the zonal flows propagate through
the convection zone up to the surface. Additionally, our results suggest that
it is the magnetic field at the tachocline that modulates the amplitude of the
surface meridional flow rather than the opposite as assumed by flux-transport
dynamo models of the solar cycle.
-
Turbulence in the weakly collisional intracluster medium of galaxies (ICM) is
able to generate strong thermal velocity anisotropies in the ions (with respect
to the local magnetic field direction), if the magnetic moment of the particles
is conserved in the absence of Coulomb collisions. In this scenario, the
anisotropic pressure magnetohydrodynamic (AMHD) turbulence shows a very
different statistical behaviour from the standard MHD one and is unable to
amplify seed magnetic fields, in disagreement with previous cosmological MHD
simulations which are successful to explain the observed magnetic fields in the
ICM. On the other hand, temperature anisotropies can also drive plasma
instabilities which can relax the anisotropy. This work aims to compare the
relaxation rate with the growth rate of the anisotropies driven by the
turbulence. We employ quasilinear theory to estimate the ions scattering rate
due to the parallel firehose, mirror, and ion-cyclotron instabilities, for a
set of plasma parameters resulting from AMHD simulations of the turbulent ICM.
We show that the ICM turbulence can sustain only anisotropy levels very close
to the instabilities thresholds. We argue that the AMHD model which bounds the
anisotropies at the marginal stability levels can describe the Alfvenic
turbulence cascade in the ICM.
-
Rotational shear layers at the boundary between radiative and convective
zones, tachoclines, play a key role in the process of magnetic field generation
in solar-like stars. We present two sets of global simulations of rotating
turbulent convection and dynamo. The first set considers a stellar convective
envelope only; the second one, aiming at the formation of a tachocline,
considers also the upper part of the radiative zone. Our results indicate that
the resulting mean-flows and dynamo properties like the growth rate, saturation
energy and mode depend on the Rossby (Ro) number. For the first set of models
either oscillatory (with ~2 yr period) or steady dynamo solutions are obtained.
The models in the second set naturally develop a tachocline which, in turn,
leads to the generation of strong mean magnetic field. Since the field is also
deposited into the stable deeper layer, its evolutionary time-scale is much
longer than in the models without a tachocline. Surprisingly, the magnetic
field in the upper turbulent convection zone evolves in the same time scale as
the deep field. These models result in either an oscillatory dynamo with ~30 yr
period or in a steady dynamo depending on Ro. In terms of the mean-field dynamo
coefficients computed using FOSA, the field evolution in the oscillatory models
without a tachocline seems to be consistent with dynamo waves propagating
according to the Parker-Yoshimura sign rule. In the models with tachoclines the
dynamics is more complex involving other transport mechanisms as well as
tachocline instabilities.
-
Realistic astrophysical environments are turbulent due to the extremely high
Reynolds numbers. Therefore, the theories of reconnection intended for
describing astrophysical reconnection should not ignore the effects of
turbulence on magnetic reconnection. Turbulence is known to change the nature
of many physical processes dramatically and in this review we claim that
magnetic reconnection is not an exception. We stress that not only
astrophysical turbulence is ubiquitous, but also magnetic reconnection itself
induces turbulence. Thus turbulence must be accounted for in any realistic
astrophysical reconnection setup. We argue that due to the similarities of MHD
turbulence in relativistic and non-relativistic cases the theory of magnetic
reconnection developed for the non-relativistic case can be extended to the
relativistic case and we provide numerical simulations that support this
conjecture. We also provide quantitative comparisons of the theoretical
predictions and results of numerical experiments, including the situations when
turbulent reconnection is self-driven, i.e. the turbulence in the system is
generated by the reconnection process itself. We show how turbulent
reconnection entails the violation of magnetic flux freezing, the conclusion
that has really far reaching consequences for many realistically turbulent
astrophysical environments. In addition, we consider observational testing of
turbulent reconnection as well as numerous implications of the theory. The
former includes the Sun and solar wind reconnection, while the latter include
the process of reconnection diffusion induced by turbulent reconnection, the
acceleration of energetic particles, bursts of turbulent reconnection related
to black hole sources as well as gamma ray bursts. Finally, we explain why
turbulent reconnection cannot be explained by turbulent resistivity or derived
through the mean field approach.
-
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.
-
List of contributions from the CTA Consortium presented at the 34th
International Cosmic Ray Conference, 30 July - 6 August 2015, The Hague, The
Netherlands.
-
Strong downstream magnetic fields of order of $\sim 1$G, with large
correlation lengths, are believed to cause the large synchrotron emission at
the afterglow phase of gamma ray bursts (GRBs). Despite of the recent
theoretical efforts, models have failed to fully explain the amplification of
the magnetic field, particularly in a matter dominated scenario. We revisit the
problem by considering the synchrotron emission to occur at the expanding shock
front of a weakly magnetized relativistic jet over a magnetized surrounding
medium. Analytical estimates and a number of high resolution 2D relativistic
magneto-hydrodynamical (RMHD) simulations are provided. Jet opening angles of
$\theta = 0^{\circ} - 20^{\circ}$, and ambient to jet density ratios of
$10^{-4} - 10^2$ were considered. We found that most of the amplification is
due to compression of the ambient magnetic field at the contact discontinuity
between the reverse and forward shocks at the jet head, with substantial
pile-up of the magnetic field lines as the jet propagates sweeping the ambient
field lines. The pile-up is maximum for $\theta \rightarrow 0$, decreasing with
$\theta$, but larger than in the spherical blast problem. Values obtained for
certain models are able to explain the observed intensities. The maximum
correlation lengths found for such strong fields is of $l_{\rm corr} \leq
10^{14}$ cm, $2 - 6$ orders of magnitude larger than the found in previous
works.
-
This paper presents the large-scale polarized sky as seen by Planck HFI at
353 GHz, which is the most sensitive Planck channel for dust polarization. We
construct and analyse large-scale maps of dust polarization fraction and
polarization direction, while taking account of noise bias and possible
systematic effects. We find that the maximum observed dust polarization
fraction is high (pmax > 18%), in particular in some of the intermediate dust
column density (AV < 1mag) regions. There is a systematic decrease in the dust
polarization fraction with increasing dust column density, and we interpret the
features of this correlation in light of both radiative grain alignment
predictions and fluctuations in the magnetic field orientation. We also
characterize the spatial structure of the polarization angle using the angle
dispersion function and find that, in nearby fields at intermediate latitudes,
the polarization angle is ordered over extended areas that are separated by
filamentary structures, which appear as interfaces where the magnetic field sky
projection rotates abruptly without apparent variations in the dust column
density. The polarization fraction is found to be anti-correlated with the
dispersion of the polarization angle, implying that the variations are likely
due to fluctuations in the 3D magnetic field orientation along the line of
sight sampling the diffuse interstellar medium.We also compare the dust
emission with the polarized synchrotron emission measured with the Planck LFI,
with low-frequency radio data, and with Faraday rotation measurements of
extragalactic sources. The two polarized components are globally similar in
structure along the plane and notably in the Fan and North Polar Spur regions.
A detailed comparison of these three tracers shows, however, that dust and
cosmic rays generally sample different parts of the line of sight and confirms
that much of the variation observed in the Planck data is due to the 3D
structure of the magnetic field.
-
For a molecular cloud clump to form stars some transport of magnetic flux is
required from the denser, inner regions to the outer regions of the cloud,
otherwise this can prevent the collapse. Fast magnetic reconnection which takes
place in the presence of turbulence can induce a process of reconnection
diffusion (RD). Extending earlier numerical studies of reconnection diffusion
in cylindrical clouds, we consider more realistic clouds with spherical
gravitational potentials and also account for the effects of the gas
self-gravity. We demonstrate that within our setup RD is efficient. We have
also identified the conditions under which RD becomes strong enough to make an
initially subcritical cloud clump supercritical and induce its collapse. Our
results indicate that the formation of a supercritical core is regulated by a
complex interplay between gravity, self-gravity, the magnetic field strength
and nearly transonic and trans-Alfv\'enic turbulence, confirming that RD is
able to remove magnetic flux from collapsing clumps, but only a few of them
become nearly critical or supercritical, sub-Alfv\'enic cores, which is
consistent with the observations. Besides, we have found that the supercritical
cores built up in our simulations develop a predominantly helical magnetic
field geometry which is also consistent with observations. Finally, we have
evaluated the effective values of the turbulent reconnection diffusion
coefficient and found that they are much larger than the numerical diffusion,
especially for initially trans-Alfv\'enic clouds, ensuring that the detected
magnetic flux removal is due to to the action of the RD rather than to
numerical diffusivity.
-
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.
-
In this study we present three-dimensional radiative cooling hydrodynamical
simulations of galactic winds generated particularly in M82-like starburst
galaxies. We have considered intermittent winds induced by SNe explosions
within super star clusters randomly distributed in the central region of the
galaxy and were able to reproduce the observed M82 wind conditions with its
complex morphological outflow structure. We have found that the environmental
conditions in the disk in nearly recent past are crucial to determine whether
the wind will develop a large scale rich filamentary structure, as in M82 wind,
or not. Also, the numerical evolution of the SN ejecta have allowed us to
obtain the abundance distribution over the first 3 kpc extension of the wind
and we have found that the SNe explosions change significantly the metallicity
only of the hot, low-density wind component. Moreover, we have found that the
SN-driven wind transports to outside the disk large amounts of energy, momentum
and gas, but the more massive high-density component reaches only intermediate
altitudes smaller than 1.5 kpc. Therefore, no significant amounts of gas mass
are lost to the IGM and the mass evolution of the galaxy is not much affected
by the starburst events occurring in the nuclear region.
-
Recent numerical analysis of Keplerian disk formation in turbulent,
magnetized cloud cores by Santos-Lima, de Gouveia Dal Pino, & Lazarian (2012)
demonstrated that reconnection diffusion is an efficient process to remove the
magnetic flux excess during the build up of a rotationally supported disk. This
process is induced by fast reconnection of the magnetic fields in a turbulent
flow. In a similar numerical study, Seifried et al. (2012) concluded that
reconnection diffusion or any other non-ideal MHD effects would not be
necessary and turbulence shear alone would provide a natural way to build up a
rotating disk without requiring magnetic flux loss. Their conclusion was based
on the fact that the mean mass-to-flux ratio ({\mu}) evaluated over a spherical
region with a radius much larger than the disk is nearly constant in their
models. In this letter we compare the two sets of simulations and show that
this averaging over large scales can mask significant real increases of {\mu}
in the inner regions where the disk is built up. We demonstrate that
turbulence-induced reconnection diffusion of the magnetic field happens in the
initial stages of the disk formation in the turbulent envelope material that is
accreting. Our analysis is suggestive that reconnection diffusion is present in
both sets of simulations and provides a simple solution for the "magnetic
braking catastrophe" which is discussed in the literature in relation to the
formation of protostellar accretion disks.
-
During structure formation, energetic events and random motions of the hot
gas residing inside galaxy clusters (the intracluster medium, ICM) generate
turbulent motions. Radio diffuse emission probes the presence of magnetic
fields and relativistic particles in the ICM, being a key ingredient for
understanding the physical processes at work in clusters of galaxies. In this
work, we present results from numerical simulations of magnetic field
turbulence in magnetohydrodynamics (MHD) and kinetic MHD (KMHD) frameworks. We
characterize magnetic field structures through the imprints left on
polarization maps.
-
In low-density plasma environments, such as the intracluster medium (ICM),
the Larmour frequency is much larger than the ion-ion collision frequency. In
such a case, the thermal pressure becomes anisotropic with respect to the
magnetic field orientation and the evolution of the turbulent gas is more
correctly described by a kinetic approach. A possible description of these
collisionless scenarios is given by the so-called kinetic magnetohydrodynamic
(KMHD) formalism, in which particles freely stream along the field lines, while
moving with the field lines in the perpendicular direction. In this way a
fluid-like behavior in the perpendicular plane is restored. In this work, we
study fast growing magnetic fluctuations in the smallest scales which operate
in the collisionless plasma that fills the ICM. In particular, we focus on the
impact of a particular evolution of the pressure anisotropy and its
implications for the turbulent dynamics of observables under the conditions
prevailing in the ICM. We present results from numerical simulations and
compare the results which those obtained using an MHD formalism.
-
The formation of protostellar disks out of molecular cloud cores is still not
fully understood. Under ideal MHD conditions, the removal of angular momentum
from the disk progenitor by the typically embedded magnetic field may prevent
the formation of a rotationally supported disk during the main protostellar
accretion phase of low mass stars. This has been known as the magnetic braking
problem and the most investigated mechanism to alleviate this problem and help
removing the excess of magnetic flux during the star formation process, the so
called ambipolar diffusion (AD), has been shown to be not sufficient to weaken
the magnetic braking at least at this stage of the disk formation. In this
work, motivated by recent progress in the understanding of magnetic
reconnection in turbulent environments, we appeal to the diffusion of magnetic
field mediated by magnetic reconnection as an alternative mechanism for
removing magnetic flux. We investigate numerically this mechanism during the
later phases of the protostellar disk formation and show its high efficiency.
By means 3D MHD simulations, we show that this mechanism is able to transport
magnetic flux to the outskirts of the disk progenitor at time scales compatible
with the collapse, allowing the formation of a rotationally supported disk
around the protostar of dimensions ~100 AU. Since MHD turbulence is expected to
be present in protostellar disks, this is a natural mechanism for removing
magnetic flux excess and allowing the formation of these disks. This mechanism
dismiss the necessity of postulating a hypothetical increase of the Ohmic
resistivity as discussed in the literature. Together with our earlier work
which showed that magnetic flux removal from molecular cloud cores is very
efficient, this work calls for reconsidering the relative role of AD for the
processes of star and planet formation.
-
The magnetic fields can change their topology through a process known as
magnetic reconnection. This process in not only important for understanding the
origin and evolution of the large-scale magnetic field, but is seen as a
possibly efficient particle accelerator producing cosmic rays mainly through
the first order Fermi process. In this work we study the properties of particle
acceleration in reconnection zones and show that the velocity component
parallel to the magnetic field of test particles inserted in nearly
non-resistive magnetohydrodynamic (MHD) domains of reconnection without
including kinetic effects, such as pressure anisotropy, the Hall term, or
anomalous effects, increases exponentially. Also, the acceleration of the
perpendicular component is always possible in such models. We have found that
within contracting magnetic islands or current sheets the particles accelerate
predominantly through the first order Fermi process, as previously described,
while outside the current sheets and islands the particles experience mostly
drift acceleration due to magnetic fields gradients. Considering two
dimensional MHD models without a guide field, we find that the parallel
acceleration stops at some level. This saturation effect is however removed in
the presence of an out-of-plane guide field or in three dimensional models.
Therefore, we stress the importance of the guide field and fully three
dimensional studies for a complete understanding of the process of particle
acceleration in astrophysical reconnection environments.
-
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.
-
Several galaxy clusters are known to present multiple and misaligned pairs of
cavities seen in X-rays, as well as twisted kiloparsec-scale jets at radio
wavelengths. It suggests that the AGN precessing jets play a role in the
formation of the misaligned bubbles. Also, X-ray spectra reveal that typically
these systems are also able to supress cooling flows, predicted theoretically.
The absence of cooling flows in galaxy clusters has been a mistery for many
years since numerical simulations and analytical studies suggest that AGN jets
are highly energetic, but are unable to redistribute it at all directions. We
performed 3D hydrodynamical simulations of the interaction between a precessing
AGN jet and the warm intracluster medium plasma, which dynamics is coupled to a
NFW dark matter gravitational potential. Radiative cooling has been taken into
account and the cooling flow problem was studied. We found that precession is
responsible for multiple pairs of bubbles, as observed. The misaligned bubbles
rise up to scales of tens of kiloparsecs, where the thermal energy released by
the jets are redistributed. After $\sim 150$ Myrs, the temperature of the gas
within the cavities is kept of order of $\sim 10^7$ K, while the denser plasma
of the intracluster medium at the central regions reaches $T \sim 10^5$ K. The
existence of multiple bubbles, at diferent directions, result in an integrated
temperature along the line of sight much larger than the simulations of
non-precessing jets. This result is in agreement with the observations. The
simulations reveal that the cooling flows cessed $\sim 50 - 70$ Myr after the
AGN jets are started.
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.
Fast magnetic reconnection may occur in different astrophysical sources, producing flare-like emission and particle acceleration. Currently, this process is being studied as an efficient mechanism to accelerate particles via a first-order Fermi process. In this work we analyse the acceleration rate and the energy distribution of test particles injected in three-dimensional magnetohydrodynamical (MHD) domains with large-scale current sheets where reconnection is made fast by the presence of turbulence. We study the dependence of the particle acceleration time with the relevant parameters of the embedded turbulence, i.e., the Alfv\'en speed $V_{\rm A}$, the injection power $P_{\rm inj}$ and scale $k_{\rm inj}$ ($k_{\rm inj} = 1/l_{\rm inj}$). We find that the acceleration time follows a power-law dependence with the particle kinetic energy: $t_{acc} \propto E^{\alpha}$, with $0.2 < \alpha < 0.6$ for a vast range of values of $c / V_{\rm A} \sim 20 - 1000$. The acceleration time decreases with the Alfv\'en speed (and therefore with the reconnection velocity) as expected, having an approximate dependence $t_{acc} \propto (V_{\rm A} / c)^{-\kappa}$, with $\kappa \sim 2.1- 2.4$ for particles reaching kinetic energies between $1 - 100 \, m_p c^2$, respectively. Furthermore, we find that the acceleration time is only weakly dependent on the $P_{\rm inj}$ and $l_{\rm inj}$ parameters of the turbulence. The particle spectrum develops a high-energy tail which can be fitted by a hard power-law already in the early times of the acceleration, in consistency with the results of kinetic studies of particle acceleration by magnetic reconnection in collisionless plasmas.
Turbulence in the weakly collisional intracluster medium of galaxies (ICM) is able to generate strong thermal velocity anisotropies in the ions (with respect to the local magnetic field direction), if the magnetic moment of the particles is conserved in the absence of Coulomb collisions. In this scenario, the anisotropic pressure magnetohydrodynamic (AMHD) turbulence shows a very different statistical behaviour from the standard MHD one and is unable to amplify seed magnetic fields, in disagreement with previous cosmological MHD simulations which are successful to explain the observed magnetic fields in the ICM. On the other hand, temperature anisotropies can also drive plasma instabilities which can relax the anisotropy. This work aims to compare the relaxation rate with the growth rate of the anisotropies driven by the turbulence. We employ quasilinear theory to estimate the ions scattering rate due to the parallel firehose, mirror, and ion-cyclotron instabilities, for a set of plasma parameters resulting from AMHD simulations of the turbulent ICM. We show that the ICM turbulence can sustain only anisotropy levels very close to the instabilities thresholds. We argue that the AMHD model which bounds the anisotropies at the marginal stability levels can describe the Alfvenic turbulence cascade in the ICM.
Rotational shear layers at the boundary between radiative and convective zones, tachoclines, play a key role in the process of magnetic field generation in solar-like stars. We present two sets of global simulations of rotating turbulent convection and dynamo. The first set considers a stellar convective envelope only; the second one, aiming at the formation of a tachocline, considers also the upper part of the radiative zone. Our results indicate that the resulting mean-flows and dynamo properties like the growth rate, saturation energy and mode depend on the Rossby (Ro) number. For the first set of models either oscillatory (with ~2 yr period) or steady dynamo solutions are obtained. The models in the second set naturally develop a tachocline which, in turn, leads to the generation of strong mean magnetic field. Since the field is also deposited into the stable deeper layer, its evolutionary time-scale is much longer than in the models without a tachocline. Surprisingly, the magnetic field in the upper turbulent convection zone evolves in the same time scale as the deep field. These models result in either an oscillatory dynamo with ~30 yr period or in a steady dynamo depending on Ro. In terms of the mean-field dynamo coefficients computed using FOSA, the field evolution in the oscillatory models without a tachocline seems to be consistent with dynamo waves propagating according to the Parker-Yoshimura sign rule. In the models with tachoclines the dynamics is more complex involving other transport mechanisms as well as tachocline instabilities.
Realistic astrophysical environments are turbulent due to the extremely high Reynolds numbers. Therefore, the theories of reconnection intended for describing astrophysical reconnection should not ignore the effects of turbulence on magnetic reconnection. Turbulence is known to change the nature of many physical processes dramatically and in this review we claim that magnetic reconnection is not an exception. We stress that not only astrophysical turbulence is ubiquitous, but also magnetic reconnection itself induces turbulence. Thus turbulence must be accounted for in any realistic astrophysical reconnection setup. We argue that due to the similarities of MHD turbulence in relativistic and non-relativistic cases the theory of magnetic reconnection developed for the non-relativistic case can be extended to the relativistic case and we provide numerical simulations that support this conjecture. We also provide quantitative comparisons of the theoretical predictions and results of numerical experiments, including the situations when turbulent reconnection is self-driven, i.e. the turbulence in the system is generated by the reconnection process itself. We show how turbulent reconnection entails the violation of magnetic flux freezing, the conclusion that has really far reaching consequences for many realistically turbulent astrophysical environments. In addition, we consider observational testing of turbulent reconnection as well as numerous implications of the theory. The former includes the Sun and solar wind reconnection, while the latter include the process of reconnection diffusion induced by turbulent reconnection, the acceleration of energetic particles, bursts of turbulent reconnection related to black hole sources as well as gamma ray bursts. Finally, we explain why turbulent reconnection cannot be explained by turbulent resistivity or derived through the mean field approach.
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.
Strong downstream magnetic fields of order of $\sim 1$G, with large correlation lengths, are believed to cause the large synchrotron emission at the afterglow phase of gamma ray bursts (GRBs). Despite of the recent theoretical efforts, models have failed to fully explain the amplification of the magnetic field, particularly in a matter dominated scenario. We revisit the problem by considering the synchrotron emission to occur at the expanding shock front of a weakly magnetized relativistic jet over a magnetized surrounding medium. Analytical estimates and a number of high resolution 2D relativistic magneto-hydrodynamical (RMHD) simulations are provided. Jet opening angles of $\theta = 0^{\circ} - 20^{\circ}$, and ambient to jet density ratios of $10^{-4} - 10^2$ were considered. We found that most of the amplification is due to compression of the ambient magnetic field at the contact discontinuity between the reverse and forward shocks at the jet head, with substantial pile-up of the magnetic field lines as the jet propagates sweeping the ambient field lines. The pile-up is maximum for $\theta \rightarrow 0$, decreasing with $\theta$, but larger than in the spherical blast problem. Values obtained for certain models are able to explain the observed intensities. The maximum correlation lengths found for such strong fields is of $l_{\rm corr} \leq 10^{14}$ cm, $2 - 6$ orders of magnitude larger than the found in previous works.
This paper presents the large-scale polarized sky as seen by Planck HFI at 353 GHz, which is the most sensitive Planck channel for dust polarization. We construct and analyse large-scale maps of dust polarization fraction and polarization direction, while taking account of noise bias and possible systematic effects. We find that the maximum observed dust polarization fraction is high (pmax > 18%), in particular in some of the intermediate dust column density (AV < 1mag) regions. There is a systematic decrease in the dust polarization fraction with increasing dust column density, and we interpret the features of this correlation in light of both radiative grain alignment predictions and fluctuations in the magnetic field orientation. We also characterize the spatial structure of the polarization angle using the angle dispersion function and find that, in nearby fields at intermediate latitudes, the polarization angle is ordered over extended areas that are separated by filamentary structures, which appear as interfaces where the magnetic field sky projection rotates abruptly without apparent variations in the dust column density. The polarization fraction is found to be anti-correlated with the dispersion of the polarization angle, implying that the variations are likely due to fluctuations in the 3D magnetic field orientation along the line of sight sampling the diffuse interstellar medium.We also compare the dust emission with the polarized synchrotron emission measured with the Planck LFI, with low-frequency radio data, and with Faraday rotation measurements of extragalactic sources. The two polarized components are globally similar in structure along the plane and notably in the Fan and North Polar Spur regions. A detailed comparison of these three tracers shows, however, that dust and cosmic rays generally sample different parts of the line of sight and confirms that much of the variation observed in the Planck data is due to the 3D structure of the magnetic field.
For a molecular cloud clump to form stars some transport of magnetic flux is required from the denser, inner regions to the outer regions of the cloud, otherwise this can prevent the collapse. Fast magnetic reconnection which takes place in the presence of turbulence can induce a process of reconnection diffusion (RD). Extending earlier numerical studies of reconnection diffusion in cylindrical clouds, we consider more realistic clouds with spherical gravitational potentials and also account for the effects of the gas self-gravity. We demonstrate that within our setup RD is efficient. We have also identified the conditions under which RD becomes strong enough to make an initially subcritical cloud clump supercritical and induce its collapse. Our results indicate that the formation of a supercritical core is regulated by a complex interplay between gravity, self-gravity, the magnetic field strength and nearly transonic and trans-Alfv\'enic turbulence, confirming that RD is able to remove magnetic flux from collapsing clumps, but only a few of them become nearly critical or supercritical, sub-Alfv\'enic cores, which is consistent with the observations. Besides, we have found that the supercritical cores built up in our simulations develop a predominantly helical magnetic field geometry which is also consistent with observations. Finally, we have evaluated the effective values of the turbulent reconnection diffusion coefficient and found that they are much larger than the numerical diffusion, especially for initially trans-Alfv\'enic clouds, ensuring that the detected magnetic flux removal is due to to the action of the RD rather than to numerical diffusivity.
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.
In this study we present three-dimensional radiative cooling hydrodynamical simulations of galactic winds generated particularly in M82-like starburst galaxies. We have considered intermittent winds induced by SNe explosions within super star clusters randomly distributed in the central region of the galaxy and were able to reproduce the observed M82 wind conditions with its complex morphological outflow structure. We have found that the environmental conditions in the disk in nearly recent past are crucial to determine whether the wind will develop a large scale rich filamentary structure, as in M82 wind, or not. Also, the numerical evolution of the SN ejecta have allowed us to obtain the abundance distribution over the first 3 kpc extension of the wind and we have found that the SNe explosions change significantly the metallicity only of the hot, low-density wind component. Moreover, we have found that the SN-driven wind transports to outside the disk large amounts of energy, momentum and gas, but the more massive high-density component reaches only intermediate altitudes smaller than 1.5 kpc. Therefore, no significant amounts of gas mass are lost to the IGM and the mass evolution of the galaxy is not much affected by the starburst events occurring in the nuclear region.
Recent numerical analysis of Keplerian disk formation in turbulent, magnetized cloud cores by Santos-Lima, de Gouveia Dal Pino, & Lazarian (2012) demonstrated that reconnection diffusion is an efficient process to remove the magnetic flux excess during the build up of a rotationally supported disk. This process is induced by fast reconnection of the magnetic fields in a turbulent flow. In a similar numerical study, Seifried et al. (2012) concluded that reconnection diffusion or any other non-ideal MHD effects would not be necessary and turbulence shear alone would provide a natural way to build up a rotating disk without requiring magnetic flux loss. Their conclusion was based on the fact that the mean mass-to-flux ratio ({\mu}) evaluated over a spherical region with a radius much larger than the disk is nearly constant in their models. In this letter we compare the two sets of simulations and show that this averaging over large scales can mask significant real increases of {\mu} in the inner regions where the disk is built up. We demonstrate that turbulence-induced reconnection diffusion of the magnetic field happens in the initial stages of the disk formation in the turbulent envelope material that is accreting. Our analysis is suggestive that reconnection diffusion is present in both sets of simulations and provides a simple solution for the "magnetic braking catastrophe" which is discussed in the literature in relation to the formation of protostellar accretion disks.
During structure formation, energetic events and random motions of the hot gas residing inside galaxy clusters (the intracluster medium, ICM) generate turbulent motions. Radio diffuse emission probes the presence of magnetic fields and relativistic particles in the ICM, being a key ingredient for understanding the physical processes at work in clusters of galaxies. In this work, we present results from numerical simulations of magnetic field turbulence in magnetohydrodynamics (MHD) and kinetic MHD (KMHD) frameworks. We characterize magnetic field structures through the imprints left on polarization maps.
In low-density plasma environments, such as the intracluster medium (ICM), the Larmour frequency is much larger than the ion-ion collision frequency. In such a case, the thermal pressure becomes anisotropic with respect to the magnetic field orientation and the evolution of the turbulent gas is more correctly described by a kinetic approach. A possible description of these collisionless scenarios is given by the so-called kinetic magnetohydrodynamic (KMHD) formalism, in which particles freely stream along the field lines, while moving with the field lines in the perpendicular direction. In this way a fluid-like behavior in the perpendicular plane is restored. In this work, we study fast growing magnetic fluctuations in the smallest scales which operate in the collisionless plasma that fills the ICM. In particular, we focus on the impact of a particular evolution of the pressure anisotropy and its implications for the turbulent dynamics of observables under the conditions prevailing in the ICM. We present results from numerical simulations and compare the results which those obtained using an MHD formalism.
The formation of protostellar disks out of molecular cloud cores is still not fully understood. Under ideal MHD conditions, the removal of angular momentum from the disk progenitor by the typically embedded magnetic field may prevent the formation of a rotationally supported disk during the main protostellar accretion phase of low mass stars. This has been known as the magnetic braking problem and the most investigated mechanism to alleviate this problem and help removing the excess of magnetic flux during the star formation process, the so called ambipolar diffusion (AD), has been shown to be not sufficient to weaken the magnetic braking at least at this stage of the disk formation. In this work, motivated by recent progress in the understanding of magnetic reconnection in turbulent environments, we appeal to the diffusion of magnetic field mediated by magnetic reconnection as an alternative mechanism for removing magnetic flux. We investigate numerically this mechanism during the later phases of the protostellar disk formation and show its high efficiency. By means 3D MHD simulations, we show that this mechanism is able to transport magnetic flux to the outskirts of the disk progenitor at time scales compatible with the collapse, allowing the formation of a rotationally supported disk around the protostar of dimensions ~100 AU. Since MHD turbulence is expected to be present in protostellar disks, this is a natural mechanism for removing magnetic flux excess and allowing the formation of these disks. This mechanism dismiss the necessity of postulating a hypothetical increase of the Ohmic resistivity as discussed in the literature. Together with our earlier work which showed that magnetic flux removal from molecular cloud cores is very efficient, this work calls for reconsidering the relative role of AD for the processes of star and planet formation.
The magnetic fields can change their topology through a process known as magnetic reconnection. This process in not only important for understanding the origin and evolution of the large-scale magnetic field, but is seen as a possibly efficient particle accelerator producing cosmic rays mainly through the first order Fermi process. In this work we study the properties of particle acceleration in reconnection zones and show that the velocity component parallel to the magnetic field of test particles inserted in nearly non-resistive magnetohydrodynamic (MHD) domains of reconnection without including kinetic effects, such as pressure anisotropy, the Hall term, or anomalous effects, increases exponentially. Also, the acceleration of the perpendicular component is always possible in such models. We have found that within contracting magnetic islands or current sheets the particles accelerate predominantly through the first order Fermi process, as previously described, while outside the current sheets and islands the particles experience mostly drift acceleration due to magnetic fields gradients. Considering two dimensional MHD models without a guide field, we find that the parallel acceleration stops at some level. This saturation effect is however removed in the presence of an out-of-plane guide field or in three dimensional models. Therefore, we stress the importance of the guide field and fully three dimensional studies for a complete understanding of the process of particle acceleration in astrophysical reconnection environments.
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.
Several galaxy clusters are known to present multiple and misaligned pairs of cavities seen in X-rays, as well as twisted kiloparsec-scale jets at radio wavelengths. It suggests that the AGN precessing jets play a role in the formation of the misaligned bubbles. Also, X-ray spectra reveal that typically these systems are also able to supress cooling flows, predicted theoretically. The absence of cooling flows in galaxy clusters has been a mistery for many years since numerical simulations and analytical studies suggest that AGN jets are highly energetic, but are unable to redistribute it at all directions. We performed 3D hydrodynamical simulations of the interaction between a precessing AGN jet and the warm intracluster medium plasma, which dynamics is coupled to a NFW dark matter gravitational potential. Radiative cooling has been taken into account and the cooling flow problem was studied. We found that precession is responsible for multiple pairs of bubbles, as observed. The misaligned bubbles rise up to scales of tens of kiloparsecs, where the thermal energy released by the jets are redistributed. After $\sim 150$ Myrs, the temperature of the gas within the cavities is kept of order of $\sim 10^7$ K, while the denser plasma of the intracluster medium at the central regions reaches $T \sim 10^5$ K. The existence of multiple bubbles, at diferent directions, result in an integrated temperature along the line of sight much larger than the simulations of non-precessing jets. This result is in agreement with the observations. The simulations reveal that the cooling flows cessed $\sim 50 - 70$ Myr after the AGN jets are started.