The rate of magnetic field diffusion plays an essential role in several astrophysical plasma processes. It has been demonstrated that the omnipresent turbulence in astrophysical media induces fast magnetic reconnection, which consequently leads to la
rge-scale magnetic flux diffusion at a rate independent of the plasma microphysics. This process is called ``reconnection diffusion (RD) and allows for the diffusion of fields which are dynamically important. The current theory describing RD is based on incompressible magnetohydrodynamic (MHD) turbulence. In this work, we have tested quantitatively the predictions of the RD theory when magnetic forces are dominant in the turbulence dynamics (Alfv{e}nic Mach number $M_A < 1$). We employed the textsc{Pencil Code} to perform numerical simulations of forced MHD turbulence, extracting the values of the diffusion coefficient $eta_{RD}$ using the Test-Field method. Our results are consistent with the RD theory ($eta_{RD} sim M_A^{3}$ for $M_A < 1$) when turbulence approaches the incompressible limit (sonic Mach number $M_S lesssim 0.02$), while for larger $M_S$ the diffusion is faster ($eta_{RD} sim M_A^{2}$). This work shows for the first time simulations of compressible MHD turbulence with the suppression of the cascade in the direction parallel to the mean magnetic field, which is consistent with incompressible weak turbulence theory. We also verified that in our simulations the energy cascading time does not follow the scaling with $M_A$ predicted for the weak regime, in contradiction with the RD theory assumption. Our results generally support and expand the RD theory predictions.
A super-massive black hole (SMBH) binary in the core of the blazar OJ 287 has been invoked in previous works to explain its observed optical flare quasi-periodicity. Following this picture, we investigate a hadronic origin for the X-ray and $gamma$-r
ay counterparts of the November 2015 major optical flare of this source. An impact outflow must result after the lighter SMBH (the secondary) crosses the accretion disc of the heavier one (the primary). We then consider acceleration of cosmic-ray (CR) protons in the shock driven by the impact outflow as it expands and collides with the active galactic nucleus (AGN) wind of the primary SMBH. We show that the emission of these CRs can reproduce the X-ray and $gamma$-ray flare data self-consistently with the optical component of the November 2015 major flare. The derived emission models are consistent with a magnetic field $B sim 5$ G in the emission region and a power-law index of $qsim2.2$ for the energy distribution of the emitting CRs. The mechanical luminosity of the AGN wind represents $lesssim 50%$ of the mass accretion power of the primary SMBH in all the derived emission profiles.
The latest flare of the regular $sim$ 12 years quasi-periodic optical outbursts in the binary SMBH candidate system OJ 287 occurred in December 2015. Following this, the source has exhibited enhanced multi-wavelength (MW) variability in spectral, tem
poral and polarization domains with new features never seen before. Our MW investigation show that the overall MW variability can be divided into two-phase, (i) November 2015 -- May 2016 with variability from near-infrared (NIR) to Fermi-LAT $rm gamma$-ray energies (0.1 -- 300 GeV), and (ii) September 2016 -- July 2017 with intense NIR to X-ray variability but without any activity in the Fermi-LAT band, and the very first detection at very high energies (VHEs, E $>$ 100 GeV) by VERITAS. The broadband SEDs during the first phase show a thermal bump in the NIR-optical region and a hardening in the $rm gamma$-ray spectra with a shift in its peak. The thermal bump like feature is consistent with the description of the standard accretion-disk associated with the primary SMBH of mass $sim 1.8times10^{10} M_odot$ while the $rm gamma$-ray emission can be naturally reproduced by inverse Compton scattering of photons from the broad line region which has been seen during the close encounter duration of the binary SMBHs, thereby suggesting a sub-parsec scale origin. The SEDs during the second phase (VHE detection) is a mixture of typical OJ 287 SED with hardened $rm gamma$-ray spectra and an HBL SED and can be explained in a two-zone model, one located at sub-parsec scales and other at parsec scales. During both the phases, the MW variability is simultaneous and almost always accompanied by changes in the polarization properties, exhibiting random and systematic variations, suggesting a strong role of magnetic field and turbulence.
We present a multi-wavelength (MW) spectral and temporal study of the recent activity of the claimed super-massive binary black hole system OJ 287 since December 2015. The overall MW activity can be divided into two durations: December 2015 - April 2
016 (MJD: 57360 - 57500), showing strong activity from near-infrared (NIR) to $gamma$-rays and September 2016 - June 2017 (MJD: 57650 - 57930), showing intense NIR to X-ray variability concurrent with detection at very high energies (VHE) by VERITAS, but without any signatures of variability in the Fermi-LAT band. In the first duration, the variations are almost simultaneous and the SEDs show new components in NIR-optical and optical-UV region. The NIR-optical bump is consistent with standard accretion-disk (AD) description while the optical-UV appears consistent with contributions from the broad-line region. The extracted broadband SEDs also show a clear shift in gamma-ray SED peak and can be explained with inverse Compton scattering of photons from broad line region. In the second period, the variations are also simultaneous except for one duration during which X-ray leads the optical/UV by $sim$ 5-6 day. The broadband SEDs, on the other hand, show mixture of a typical OJ 287 SED and an HBL SED, consistent with an origin from two different zones, one located at sub-parsec scales and other at parsec scales.
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 or
igin 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.
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 a
nd 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.
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 cluster
s 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.
NGC1275, the central galaxy in the Perseus cluster, is the host of gigantic hot bipolar bubbles inflated by AGN jets observed in the radio as Perseus A. It presents a spectacular $H{alpha}$-emitting nebulosity surrounding NGC1275, with loops and fila
ments of gas extending to over 50 kpc. The origin of the filaments is still unknown, but probably correlates with the mechanism responsible for the giant buoyant bubbles. We present 2.5 and 3-dimensional MHD simulations of the central region of the cluster in which turbulent energy, possibly triggered by star formation and supernovae (SNe) explosions is introduced. The simulations reveal that the turbulence injected by massive stars could be responsible for the nearly isotropic distribution of filaments and loops that drag magnetic fields upward as indicated by recent observations. Weak shell-like shock fronts propagating into the ICM with velocities of 100-500 km/s are found, also resembling the observations. The isotropic outflow momentum of the turbulence slows the infall of the intracluster medium, thus limiting further starburst activity in NGC1275. As the turbulence is subsonic over most of the simulated volume, the turbulent kinetic energy is not efficiently converted into heat and additional heating is required to suppress the cooling flow at the core of the cluster. Simulations combining the MHD turbulence with the AGN outflow can reproduce the temperature radial profile observed around NGC1275. While the AGN mechanism is the main heating source, the supernovae are crucial to isotropize the energy distribution.
In this work, considering the impact of a SNR with a neutral magnetized cloud we derived analytically a set of conditions which are favorable for driving gravitational instability in the cloud and thus star formation. We have built diagrams of the SN
R radius, versus the cloud density, that constrain a domain in the parameter space where star formation is allowed. The diagrams are also tested with fully 3-D MHD simulations involving a SNR and a self-gravitating cloud and we find that the numerical analysis is consistent with the results predicted by the diagrams. While the inclusion of a homogeneous magnetic field approximately perpendicular to the impact velocity of the SNR with an intensity ~1 $ mu$G results only a small shrinking of the star formation triggering zone in the diagrams, a larger magnetic field (~10 $ mu$G) causes a significant shrinking, as expected. Applications of the diagrams to a few regions of our own galaxy have revealed that star formation in those sites could have been triggered by shock waves from SNRs. Finally, we have evaluated the effective star formation efficiency for this sort of interaction and found that it is smaller than the observed values in our own Galaxy (sfe ~0.01-0.3). This result is consistent with previous work in the literature and also suggests that the mechanism presently investigated, though very powerful to drive structure formation, supersonic turbulence and eventually, local star formation, does not seem to be sufficient to drive global star formation in normal star forming galaxies, not even when the magnetic field in the neutral clouds is neglected. (abridged)
The sunspot solar cycle has been usually explained as the result of a dynamo process operating in the sun. This is a classical problem in Astrophysics that until the present is not fully solved. Here we discuss current problems and limitations with t
he solar dynamo modeling and their possible solutions using a kinematic dynamo model within the Babcock-Leighton approximation as a tool. In particular, we discuss the importance of the turbulent magnetic pumping versus the meridional flow circulation in the dynamo operation.
