In the outskirts of the intracluster medium (ICM) in galaxy clusters, the temperature decreases with radius. Due to the weakly collisional nature of the plasma, these regions are susceptible to the magneto-thermal instability (MTI), which can sustain turbulence and provide turbulent pressure support in the ICM. This instability arises due to heat conduction directed along the magnetic field, with a heat conductivity which is normally assumed to be given by the Spitzer value. Recent numerical studies of the ion mirror and the electron whistler instability using particle-in-cell codes have shown that microscale instabilities can lead to a reduced value for the heat conductivity in the ICM. This could in turn influence the efficiency with which the MTI drives turbulence. In this paper we investigate the influence of reduced heat transport on the nonlinear evolution of the MTI by employing a subgrid model that mimics the influence of the mirror instability on the heat conductivity. We use this subgrid model to assess the effect of microscales on the large scale dynamics of the ICM. We find that the nonlinear saturation of the MTI is surprisingly robust. Over a factor of $\sim 10^3$ the thermal-to-magnetic pressure ratio and collisionality we find at most modest changes to the saturation of the MTI with respect to reference simulations where heat transport is unsuppressed.

The assumption of a spatially uniform helium distribution in the intracluster medium can lead to biases in the estimates of key cluster parameters if composition gradients are present. The helium concentration profile in galaxy clusters is unfortunately not directly observable. Current models addressing the putative sedimentation are one-dimensional and parametrize the presence of magnetic fields in a crude way, ignoring the weakly-collisional, magnetized nature of the medium. When these effects are considered, a wide variety of instabilities can play an important role in the plasma dynamics. In a series of recent papers, we have developed the local, linear theory of these instabilities and addressed their non-linear development with a modified version of Athena. Here, we extend our study by developing a quasi-global approach that we use to simulate the mixing of helium as induced by generalizations of the heat-flux-driven buoyancy instability (HBI) and the magneto-thermal instability (MTI), which feed off thermal and composition gradients. In the inner region of the ICM, mixing can occur on few Gyrs, after which the average magnetic field inclination angle is $\sim 45^{\circ}$ resulting in an averaged Spitzer parameter higher by about 20 % than the value obtained in homogeneous simulations. In the cluster outskirts the instabilities are rather inefficient, due to the shallow gradients. This suggests that compositions gradients in cluster cores might be shallower than one-dimensional models predict. More quantitative statements demand more refined models that can incorporate the physics driving the sedimentation process and simultaneously account for the weakly-collisional nature of the plasma.

The distribution of Helium in the intracluster medium (ICM) permeating galaxy clusters is not well constrained due to the very high plasma temperature. Therefore, the plasma is often assumed to be homogeneous. A non-uniform Helium distribution can however lead to biases when measuring key cluster parameters. This has motivated one-dimensional models that evolve the ICM composition assuming that the effects of magnetic fields can be parameterized or ignored. Such models for non-isothermal clusters show that Helium can sediment in the cluster core leading to a peak in concentration offset from the cluster center. The resulting profiles have recently been shown to be linearly unstable when the weakly-collisional character of the magnetized plasma is considered. In this paper, we present a modified version of the MHD code Athena, which makes it possible to evolve a weakly-collisional plasma subject to a gravitational field and stratified in both temperature and composition. We thoroughly test our implementation and confirm excellent agreement against several analytical results. In order to isolate the effects of composition, in this initial study we focus our attention on isothermal plasmas. We show that plasma instabilities, feeding off gradients in composition, can induce turbulent mixing and saturate by re-arranging magnetic field lines and alleviating the composition gradient. Composition profiles that increase with radius lead to instabilities that saturate by driving the average magnetic field inclination to roughly $45^{\circ}$. We speculate that this effect may alleviate the core insulation observed in homogeneous settings, with potential consequences for the associated cooling flow problem.

We compare transport statistics of elongated incompressible shearing boxes for different Reynolds and magnetic Prandtl numbers, $Re$ and $Pm$, and aspect ratios, $L_z/L_x$. We find that at fixed aspect ratio $L_z/L_x=4$ and $Re = 10,000$, the turbulent stresses for $Pm \lesssim 1$ do not show considerable variation and follow a power law $\sim Pm^{3/2}$ for $Pm > 1$. This is qualitatively consistent with previous results based on net imposed flux and small box $L_z/L_x \sim 1$ simulations but the power law exponent is different. The saturated level of stresses, the ratio of Maxwell stress to the magnetic energy and Reynolds to Maxwell stress ratio are roughly invariant as $L_z/L_x$ is increased. For cases where the boxes are elongated in both the azimuth and vertical direction, the transport coefficient $\alpha \in [0.1,1.0]$ that is $10-100$ times larger than the case with $L_y/L_x = 2$ and large $L_z/L_x$. Overall, our results suggest that the saturated state of turbulence is sensitive to both dissipation coefficients and aspect ratio (both $L_z/L_x$, $L_y/L_x$) motivating further work on this problem.

Jupiter and Saturn formed early, before the gas disk dispersed. The presence of gap-opening planets affects the dynamics of the gas and embedded solids and halts the inward drift of grains above a certain size. A drift barrier can explain the absence of calcium aluminium rich inclusions (CAIs) in chondrites originating from parent bodies that accreted in the inner solar system. Employing an interdisciplinary approach, we use a $\mu$-X-Ray-fluorescence scanner to search for large CAIs and a scanning electron microscope to search for small CAIs in the ordinary chondrite NWA 5697. We carry out long-term, two-dimensional simulations including gas, dust, and planets to characterize the transport of grains within the viscous $\alpha$-disk framework exploring the scenarios of a stand-alone Jupiter, Jupiter and Saturn \textit{in situ}, or Jupiter and Saturn in a 3:2 resonance. In each case, we find a critical grain size above which drift is halted as a function of the physical conditions in the disk. From the laboratory search we find four CAIs with a largest size of $\approx$200$\,\mu$m. \Combining models and data, we provide an estimate for the upper limit of the $\alpha$-viscosity and the surface density at the location of Jupiter, using reasonable assumptions about the stellar accretion rate during inward transport of CAIs, and assuming angular momentum transport to happen exclusively through viscous effects. Moreover, we find that the compound gap structure in the presence of Saturn in a 3:2 resonance favors inward transport of grains larger than CAIs currently detected in ordinary chondrites.

We present an asymptotically and unconditionally stable numerical method to account for the momentum transfer between multiple species. Momentum is conserved to machine precision. This implies that the asymptotic equilibrium corresponds to the velocity of the center of mass. Aimed at studying dust dynamics, we implement this numerical method in the publicly available code FARGO3D. To validate our implementation, we develop a test suite for an arbitrary number of species, based on analytical or exact solutions of problems related to perfect damping, damped sound waves, shocks, local and global gas-dust radial drift in a disk and, linear streaming instability. In particular, we obtain first-order, steady-state solutions for the radial drift of multiple dust species in protoplanetary disks, in which the pressure gradient is not necessarily small. We additionally present non-linear shearing-box simulations of the streaming instability and compare them with previous results obtained with Lagrangian particles. We successfully validate our implementation by recovering the solutions from the test suite to second- and first-order accuracy in space and time, respectively. From this, we conclude that our scheme is suitable, and very robust, to study the self-consistent dynamics of several fluids. In particular, it can be used for solving the collisions between gas and dust in protoplanetary disks, with any degree of coupling.

The streaming instability is thought to play a central role in the early stages of planet formation by enabling the efficient bypass of a number of barriers hindering the formation of planetesimals. We present the first study exploring the efficiency of the linear streaming instability when a particle-size distribution is considered. We find that, for a given dust-to-gas mass ratio, the multi-species streaming instability grows on timescales much longer than those expected when only one dust species is involved. In particular, distributions that contain close-to-order-unity dust-to-gas mass ratios lead to unstable modes that can grow on timescales comparable, or larger, with those of secular instabilities. We anticipate that processes leading to particle segregation and/or concentration can create favourable conditions for the instability to grow fast. Our findings may have important implications for a large number of processes in protoplanetary disks that rely on the streaming instability as usually envisioned for a unique dust species. Our results suggest that the growth rates of other resonant-drag-instabilities may also decrease considerably when multiple species are considered.

Imaging of the dust continuum emitted from disks around nearby protostars reveals diverse substructure. In recent years, theoretical efforts have been intensified to investigate how far the intrinsic dynamics of protoplanetary disks (PPDs) can lead to such features. Turbulence in the realm of non-ideal magnetohydrodynamics (MHD) is one candidate for explaining the generation of zonal flows which can lead to local dust enhancements. Adopting a radially varying cylindrical disk model, and considering combinations of vertical and azimuthal initial net flux, we perform 3D non-ideal MHD simulations aimed at studying self-organization induced by the Hall effect in turbulent PPDs. To this end, new modules have been incorporated into the NIRVANA-III and FARGO3D MHD codes. We moreover include dust grains, treated in the fluid approximation, in order to study their evolution subject to the emerging zonal flows. In the regime of a dominant Hall effect, we robustly obtain large-scale organized concentrations in the vertical magnetic field that remain stable for hundreds of orbits. For disks with vertical initial net flux alone, we confirm the presence of zonal flows and vortices that introduce regions of super-Keplerian gas flow. Including a moderately strong net-azimuthal magnetic flux can significantly alter the dynamics, partially preventing the self-organization of zonal flows. For plasma beta-parameters smaller than 50, large-scale, near-axisymmetric structures develop in the vertical magnetic flux. In all cases, we demonstrate that the emerging features are capable of accumulating dust grains for a range of Stokes numbers.

Fast inward migration of planetary cores is a common problem in the current planet formation paradigm. Even though dust is ubiquitous in protoplanetary disks, its dynamical role in the migration history of planetary embryos has not been assessed. In this Letter, we show that the scattered pebble-flow induced by a low-mass planetary embryo leads to an asymmetric dust-density distribution that is able to exert a net torque. By analyzing a large suite of multifluid hydrodynamical simulations addressing the interaction between the disk and a low-mass planet on a fixed circular orbit, and neglecting dust feedback onto the gas, we identify two different regimes, gas- and gravity-dominated, where the scattered pebble-flow results in almost all cases in positive torques. We collect our measurements in a first torque map for dusty disks, which will enable the incorporation of the effect of dust dynamics on migration into population synthesis models. Depending on the dust drift speed, the dust-to-gas mass ratio/distribution and the embryo mass, the dust-induced torque has the potential to halt inward migration or even induce fast outward migration of planetary cores. We thus anticipate that dust-driven migration could play a dominant role during the formation history of planets. Because dust torques scale with disk metallicity, we propose that dust-driven outward migration may enhance the occurrence of distant giant planets in higher-metallicity systems.

We present an analytical model to investigate the production of pebbles and their radial transport through a protoplanetary disk (PPD) with magnetically driven winds. While most of the previous analytical studies in this context assume that the radial turbulent coefficient is equal to the vertical dust diffusion coefficient, in the light of the results of recent numerical simulations, we relax this assumption by adopting effective parametrisations of the turbulent coefficients involved in terms of the strength of the magnetic fields driving the wind. Theoretical studies have already pointed out that even in the absence of winds, these coefficients are not necessarily equal, though its consequences regarding pebble production have not been explored. In this paper, we investigate the evolution of the pebble production line, the radial mass flux of the pebbles and their corresponding surface density as a function of the plasma parameter at the disk midplane. Our analysis explicitly demonstrates that the presence of magnetically-driven winds in a PPD leads to considerable reduction of the rate and duration of the pebble delivery. We show that when the wind is strong, the core growth in mass due to the pebble accretion is so slow that it is unlikely that a core could reach a pebble isolation mass during a PPD lifetime. When the mass of a core reaches this critical value, pebble accretion is halted due to core-driven perturbations in the gas. With decreasing wind strength, however, pebble accretion may, in a shorter time, increase the mass of a core to the pebble isolation mass.

We apply the picket fence treatment to model the effects brought about by spectral lines on the thermal structure of irradiated atmospheres. The lines may be due to purely absorption processes, purely coherent scattering processes or some combination of absorption and scattering. If the lines arise as a pure absorption process, the surface layers of the atmosphere are cooler whereas this surface cooling is completely absent if the lines are due to pure coherent isotropic scattering. The lines also lead to a warming of the deeper atmosphere. The warming of the deeper layers is, however, independent of the nature of line formation. Accounting for coherent isotropic scattering in the shortwave and longwave continuum results in anti-greenhouse cooling and greenhouse warming on an atmosphere-wide scale. The effects of coherent isotropic scattering in the line and continuum operate in tandem to determine the resulting thermal structure of the irradiated atmosphere.

Aerodynamic theory predicts that dust grains in protoplanetary disks will drift radially inward on comparatively short timescales. In this context, it has long been known that the presence of a gap opened by a planet can alter the dust dynamics significantly. In this paper, we carry out a systematic study employing long-term numerical simulations aimed at characterizing the critical particle-size for retention outside a gap as a function of particle size and for various key parameters defining the protoplanetary disk model. To this end, we perform multifluid hydrodynamical simulations in two dimensions, including different dust species, which we treat as pressureless fluids. We initialize the dust outside of the planet's orbit and study under which conditions dust grains are able to cross the gap carved by the planet. In agreement with previous work, we find that the permeability of the gap depends both on dust dynamical properties and the gas disk structure: while small dust follows the viscously accreting gas through the gap, dust grains approaching a critical size are progressively filtered out. Moreover, we introduce and compute a depletion factor that enables us to quantify the way in which higher viscosity, smaller planet mass, or a more massive disk can shift this critical size to larger values. Our results indicate that gap-opening planets may act to deplete the inner reaches of protoplanetary disks of large dust grains -- potentially limiting the accretion of solids onto forming terrestrial planets.

We investigate the stability of test-particle equilibrium orbits in axisymmetric, but otherwise arbitrary, gravitational and electromagnetic fields. We extend previous studies of this problem to include a toroidal magnetic field. We find that, even though the toroidal magnetic field does not alter the location of the circular orbits, it enters the problem as a gyroscopic force with the potential to provide gyroscopic stability. This is in essence similar to the situation encountered in the reduced three-body problem where rotation enables stability around the local maxima of the effective potential. Nevertheless, we show that gyroscopic stabilization by a toroidal magnetic field is impossible for axisymmetric force fields in source-free regions because in this case the effective potential does not possess any local maxima. As an example of an axisymmetric force field with sources, we consider the classical problem of a rotating, aligned magnetosphere. By analyzing the dynamics of halo and equatorial particle orbits we conclude that axisymmetric toroidal fields that are antisymmetric about the equator are unable to provide gyroscopic stabilization. On the other hand, a toroidal magnetic field that does not vanish at the equator can provide gyroscopic stabilization for positively charged particles in prograde equatorial orbits.

We present a new set of analytical solutions to model the steady state structure of a protoplanetary disc with a magnetically-driven wind. Our model implements a parametrization of the stresses involved and the wind launching mechanism in terms of the plasma parameter at the disc midplane, as suggested by the results of recent, local MHD simulations. When wind mass-loss is accounted for, we find that its rate significantly reduces the disc surface density, particularly in the inner disc region. We also find that models that include wind mass-loss lead to thinner dust layers. As an astrophysical application of our models, we address the case of HL Tau, whose disc exhibits a high accretion rate and efficient dust settling at its midplane. These two observational features are not easy to reconcile with conventional accretion disc theory, where the level of turbulence needed to explain the high accretion rate would prevent a thin dust layer. Our disc model that incorporates both mass-loss and angular momentum removal by a wind is able to account for HL Tau observational constraints concerning its high accretion rate and dust layer thinness.

The effect of magnetic field diffusion on the stability of accretion disks is a problem that has attracted considerable interest of late. In particular, the Hall effect has the potential to bring about remarkable changes in the dynamical behavior of disks that are without parallel. In this paper, we conduct a systematic examination of the linear eigenmodes in a weakly magnetized differentially rotating gas with special focus on Hall diffusion. We first develop a geometrical representation of the eigenmodes and provide a detailed quantitative description of the polarization properties of the oscillatory modes under the combined influence of the Coriolis and Hall effects. We also analyze the effects of magnetic diffusion on the structure of the unstable modes and derive analytical expressions for the kinetic and magnetic stresses and energy densities associated with the non-ideal MRI. Our analysis explicitly demonstrates that, if the dissipative effects are relatively weak, the kinetic stresses and energies make up the dominant contribution to the total stress and energy density when the equilibrium angular momentum and magnetic field vectors are anti-parallel. This is in sharp contrast to what is observed in the case of the ideal or dissipative MRI. We conduct shearing box simulations and find very good agreement with the results derived from linear analysis. As the modes in consideration are also exact solutions of the non-linear equations, the unconventional nature of the kinetic and magnetic stresses may have significant implications for the non-linear evolution in some regions of protoplanetary disks.

We show for the first time that sustained turbulence is possible at low magnetic Prandtl number for Keplerian flows with no mean magnetic flux. Our results indicate that increasing the vertical domain size is equivalent to increasing the dynamical range between the energy injection scale and the dissipative scale. This has important implications for a large variety of differentially rotating systems with low magnetic Prandtl number such as protostellar disks and laboratory experiments.

We use the machinery usually employed for studying the onset of Rayleigh--B\'enard convection in hydro- and magnetohydro-dynamic settings to address the onset of convection induced by the magnetothermal instability and the heat-flux-buoyancy-driven-instability in the weakly-collisional magnetized plasma permeating the intracluster medium. Since most of the related numerical simulations consider the plasma being bounded between two `plates' on which boundary conditions are specified, our strategy provides a framework that could enable a more direct connection between analytical and numerical studies. We derive the conditions for the onset of these instabilities considering the effects of induced magnetic tension resulting from a finite plasma beta. We provide expressions for the Rayleigh number in terms of the wave vector associated with a given mode, which allow us to characterize the modes that are first to become unstable. For both the heat-flux-buoyancy-driven-instability and the magnetothermal instability, oscillatory marginal stable states are possible.

We devise and explore an iterative optimization procedure for controlling particle populations in particle-in-cell (PIC) codes via merging and splitting of computational macro-particles. Our approach, is to compute an optimal representation of the global particle phase space structure while decreasing or increasing the entire particle population, based on k-means clustering of the data. In essence the procedure amounts to merging or splitting particles by statistical means, throughout the entire simulation volume in question, while minimizing a 6-dimensional total distance measure to preserve the physics. Particle merging is by far the most demanding procedure when considering conservation laws of physics; it amounts to lossy compression of particle phase space data. We demonstrate that our k-means approach conserves energy and momentum to high accuracy, even for high compression ratios, $\mathcal{R} \approx 3$ --- \emph{i.e.}, $N_{f} \lesssim 0.33N_{i}$. Interestingly, we find that an accurate particle splitting step can be performed using k-means as well; this from an argument of symmetry. The split solution, using k-means, places splitted particles optimally, to obtain maximal spanning on the phase space manifold. Implementation and testing is done using an electromagnetic PIC code, the \ppcode. Nonetheless, the k-means framework is general; it is not limited to Vlasov-Maxwell type PIC codes. We discuss advantages and drawbacks of this optimal phase space reconstruction.

Understanding whether Helium can sediment to the core of galaxy clusters is important for a number of problems in cosmology and astrophysics. All current models addressing this question are one-dimensional and do not account for the fact that magnetic fields can effectively channel ions and electrons, leading to anisotropic transport of momentum, heat, and particle diffusion in the weakly collisional intracluster medium (ICM). This anisotropy can lead to a wide variety of instabilities, which could be relevant for understanding the dynamics of heterogeneous media. In this paper, we consider the radial temperature and composition profiles as obtained from a state-of-the-art Helium sedimentation model and analyze its stability properties. We find that the associated radial profiles are unstable, to different kinds of instabilities depending on the magnetic field orientation, at all radii. The fastest growing modes are usually related to generalizations of the Magnetothermal Instability (MTI) and the Heat-flux-driven Buoyancy Instability (HBI) which operate in heterogeneous media. We find that the effect of sedimentation is to increase (decrease) the predicted growth rates in the inner (outer) cluster region. The unstable modes grow fast compared to the sedimentation timescale. This suggests that the composition gradients as inferred from sedimentation models, which do not fully account for the anisotropic character of the weakly collisional environment, might not be very robust. Our results emphasize the subtleties involved in understanding the gas dynamics of the ICM and argue for the need of a comprehensive approach to address the issue of Helium sedimentation beyond current models.

We conduct a comprehensive axisymmetric, local linear mode analysis of a stratified, differentially rotating disk permeated by a toroidal magnetic field which could provide significant pressure support. In the adiabatic limit, we derive a new stability criteria that differs from the one obtained for weak magnetic fields with a poloidal component and reduces continuously to the hydrodynamic Solberg-H{\o}iland criteria. Three fundamental unstable modes are found in the locally isothermal limit. They comprise of overstable: (i) acoustic oscillations, (ii) radial epicyclic (acoustic-inertial) oscillations and (iii) vertical epicyclic (or vertical shear) oscillations. All three modes are present for finite ranges of cooling times but they are each quickly quenched past respective cut-off times. The acoustic and acoustic-inertial overstable modes are driven by the background temperature gradient. When vertical structure is excluded, we find that the radial epicyclic modes appear as a nearly degenerate pair. One of these is the aforementioned acoustic-inertial mode and the other has been previously identified in a slightly different guise as the convective overstability. Inclusion of vertical structure leads to the development of overstable oscillations destabilized by vertical shear but also has the effect of suppressing the radial epicyclic modes. Although our study does not explicitly account for non-ideal effects, we argue that it may still shed light into the dynamics of protoplanetary disk regions where a strong toroidal field generates as a result of the Hall-shear instability.

Disks with a barotropic equilibrium structure, for which the pressure is only a function of the density, rotate on cylinders in the presence of a gravitational potential, so that the angular frequency of such a disk is independent of height. Such disks with barotropic equilibria can be approximately modeled using the shearing box framework, representing a small disk volume with height-independent angular frequency. If the disk is in baroclinic equilibrium, the angular frequency does generally depend on height, and it is thus necessary to go beyond the standard shearing box approach. In this paper, we show that given a global disk model, it is possible to develop approximate models that are local in horizontal planes without an expansion in height with shearing-periodic boundary conditions. We refer to the resulting framework as the vertically global shearing box (VGSB). These models can be non-axisymmetric for globally barotropic equilibria but should be axisymmetric for globally baroclinic equilibria. We provide explicit equations for this VGSB which can be implemented in standard magnetohydrodynamic codes by generalizing the shearing-periodic boundary conditions to allow for a height-dependent angular frequency and shear rate. We also discuss the limitations that result from the radial approximations that are needed in order to impose height-dependent shearing periodic boundary conditions. We illustrate the potential of this framework by studying a vertical shear instability and examining the modes associated with the magnetorotational instability.

The formation and evolution of a wide class of astrophysical objects is governed by turbulent, magnetized accretion disks. Understanding their secular dynamics is of primary importance. Apart from enabling mass accretion via the transport of angular momentum, the turbulence affects the long-term evolution of the embedded magnetic flux, which in turn regulates the efficiency of the transport. In this paper, we take a comprehensive next step towards an effective mean-field model for turbulent astrophysical disks by systematically studying the key properties of magnetorotational turbulence in vertically-stratified, isothermal shearing boxes. This allows us to infer emergent properties of the ensuing chaotic flow as a function of the shear parameter as well as the amount of net-vertical flux. Using the test-field method, we furthermore characterize the mean-field dynamo coefficients that describe the long-term evolution of large-scale fields. We simultaneously infer the vertical shape and the spectral scale dependence of these closure coefficients, with the latter describing non-local contributions to the turbulent electromotive force. Based on this, we infer a scale-separation ratio of about ten for the large-scale dynamo. We finally determine scaling properties of the mean-field dynamo coefficients. The relevant component of the dynamo {\alpha} effect is found to scale linearly with the shear rate, as is the corresponding turbulent diffusion, {\eta}. Together, these scalings allow us to predict, in a quantitative manner, the cycle period of the well-known butterfly diagram. This lends new support to the importance of the {\alpha}{\Omega} mechanism in determining the evolution of large-scale magnetic fields in turbulent accretion disks.

One of the most important predictions of any gap model for pulsar magnetospheres is the predicted $\gamma$-ray spectra. In the outer gap model, the properties of the synchro-curvature radiation are sensitive to many parameters, whose realistic ranges have been studied in detail in an accompanying paper. There we demonstrated that the uncertainty in the radius of curvature, the magnetic field geometry, and the X-ray surface flux may affect by orders of magnitude the predicted flux and spectral peak in the $\gamma$-ray regime. Here, we present a systematic, numerical study of the impact of the different parameters on the particle dynamics along the gap and calculate the emitted synchro-curvature radiation along the trajectory. By integrating the emitted radiation along the gap and convolving it with a parametrized particle distribution, we discuss how the comparison with the wealth of {\em Fermi}-LAT data can be used to constrain the applicability of the model. The resulting spectra show very different energy peaks, fluxes and shapes, qualitatively matching the great variety of the observed {\em Fermi}-LAT pulsars. In particular, if we see a large fraction of photons emitted from the initial part of the trajectory, we show that the spectra will be flatter at the low-energy {\it Fermi}-LAT regime (100 MeV -- 1 GeV). This provides a solution for such observed flat spectra, while still maintain synchro-curvature radiation as the origin of these photons.

We consider the fundamental problem of charged particles moving along and around a curved magnetic field line, revising the synchro-curvature radiation formulae introduced by Cheng and Zhang (1996). We provide more compact expressions to evaluate the spectrum emitted by a single particle, identifying the key parameter that controls the transition between the curvature-dominated and the synchrotron-dominated regime. This parameter depends on the local radius of curvature of the magnetic field line, the gyration radius, and the pitch angle. We numerically solve the equations of motion for the emitting particle by considering self-consistently the radiative losses, and provide the radiated spectrum produced by a particle when an electric acceleration is balanced by its radiative losses, as it is assumed to happen in the outer gaps of pulsar's magnetospheres. We compute the average spectrum radiated throughout the particle trajectory finding that the slope of the spectrum before the peak depends on the location and size of the emission region. We show how this effect could then lead to a variety of synchro-curvature spectra. Our results reinforce the idea that the purely synchrotron or curvature losses are, in general, inadequate to describe the radiative reaction on the particle motion, and the spectrum of emitted photons. Finally, we discuss the applicability of these calculations to different astrophysical scenarios.

The popular outer gap model of magnetospheric emission from pulsars has been widely applied to explain the properties observed in $\gamma$-rays. However, its quantitative predictions rely on a number of approximations and assumptions that are usually overlooked. Here we examine them, reviewing the main ingredients entering in the model, evaluating their range of uncertainties. Usually, in the quantitative applications of the model, key parameters like the radius of curvature and the energies of the interacting photons are taken to be a fixed, single value. Instead, here we explore their realistic ranges, and the impact of these on the consistency of the model itself. We conclude that the popular evaluation of the trans-field size of the gap as a function of period and period derivative, is unreliable and affected by a huge dispersion. Last, the exploration of the possible values for the radius of curvature, the local magnetic field and other quantities deserve more attention for quantitative applications of the outer gap model, like the calculation of $\gamma$-ray spectra, which is the subject of an accompanying paper.

The magnetorotational instability is thought to play an important role in enabling accretion in sufficiently ionized astrophysical disks. The rate at which MRI-driven turbulence transports angular momentum is related to both the strength of the amplitudes of the fluctuations on various scales and the degree of anisotropy of the underlying turbulence. This has motivated several studies of the distribution of turbulent power in spectral space. In this paper, we investigate the anisotropic nature of MRI-driven turbulence using a pseudo-spectral code and introduce novel ways to robustly characterize the underlying turbulence. We show that the general flow properties vary in a quasi-periodic way on timescales comparable to 10 inverse angular frequencies motivating the temporal analysis of its anisotropy. We introduce a 3D tensor invariant analysis to quantify and classify the evolution of the anisotropic turbulent flow. This analysis shows a continuous high level of anisotropy, with brief sporadic transitions towards two- and three-component isotropic turbulent flow. This temporal-dependent anisotropy renders standard shell-average, especially when used simultaneously with long temporal averages, inadequate for characterizing MRI-driven turbulence. We propose an alternative way to extract spectral information from the turbulent magnetized flow, whose anisotropic character depends strongly on time. This consists of stacking 1D Fourier spectra along three orthogonal directions that exhibit maximum anisotropy in Fourier space. The resulting averaged spectra show that the power along each of the three independent directions differs by several orders of magnitude over most scales, except the largest ones. Our results suggest that a first-principles theory to describe fully developed MRI-driven turbulence will likely have to consider the anisotropic nature of the flow at a fundamental level.

Over the last decade, substantial efforts have been devoted to understanding the stability properties, transport phenomena, and long-term evolution of weakly-collisional, magnetized plasmas which are stratified in temperature. These studies have improved our understanding of the physics governing the intra-cluster medium (ICM), but assumed that ICM is a homogeneous. This, however, might not be a good approximation if heavy elements sediment in the inner region of the galaxy cluster. In this paper, we analyze the stability of a weakly-collisional, magnetized plane-parallel atmosphere which is stratified in both temperature and composition. This allows us to discuss for the first time the dynamics of weakly-collisional environments where heat conduction, momentum transport, and ion-diffusion are anisotropic with respect to the direction of the magnetic field. We show that, depending on the relative signs and magnitudes of the gradients in the temperature and the mean molecular weight, the plasma can be subject to a wide variety of unstable modes which include modifications to the magnetothermal instability (MTI), the heat-flux-driven buoyancy instability (HBI), and overstable gravity modes previously studied in homogeneous media. We discuss the astrophysical implications of our findings for a representative galaxy cluster where helium has sedimented. Our findings suggest that the core insulation that results from the magnetic field configurations that arise as a natural consequence of the HBI, which would be MTI stable in a homogeneous medium, could be alleviated if the mean molecular weight gradient is steep enough, i.e., $(\nabla \mu)/\mu > (\nabla T)/T$. This study constitutes a first step toward understanding the interaction between magnetic turbulence and the diffusion of heavy elements in the ICM. (abridged)

The physical modeling of the accretion disk boundary layer, the region where the disk meets the surface of the accreting star, usually relies on the assumption that angular momentum transport is opposite to the radial angular frequency gradient of the disk. The standard model for turbulent shear viscosity, widely adopted in astrophysics, satisfies this assumption by construction. However, this behavior is not supported by numerical simulations of turbulent magnetohydrodynamic (MHD) accretion disks, which show that angular momentum transport driven by the magnetorotational instability is inefficient in this inner disk region. I will discuss the results of a recent study on the generation of hydromagnetic stresses and energy density in the boundary layer around a weakly magnetized star. Our findings suggest that although magnetic energy density can be significantly amplified in this region, angular momentum transport is rather inefficient. This seems consistent with the results obtained in numerical simulations and suggests that the detailed structure of turbulent MHD boundary layers could differ appreciably from those derived within the standard framework of turbulent shear viscosity.

Detailed calculations of the physical structure of accretion disk boundary layers, and thus their inferred observational properties, rely on the assumption that angular momentum transport is opposite to the radial angular frequency gradient of the disk. The standard model for turbulent shear viscosity satisfies this assumption by construction. However, this behavior is not supported by numerical simulations of turbulent magnetohydrodynamic (MHD) accretion disks, which show that angular momentum transport driven by the magnetorotational instability (MRI) is inefficient in disk regions where, as expected in boundary layers, the angular frequency increases with radius. In order to shed light into physically viable mechanisms for angular momentum transport in this inner disk region, we examine the generation of hydromagnetic stresses and energy density in differentially rotating backgrounds with angular frequencies that increase outward in the shearing-sheet framework. We isolate the modes that are unrelated to the standard MRI and provide analytic solutions for the long-term evolution of the resulting shearing MHD waves. We show that, although the energy density of these waves can be amplified significantly, their associated stresses oscillate around zero, rendering them an inefficient mechanism to transport significant angular momentum (inward). These findings are consistent with the results obtained in numerical simulations of MHD accretion disk boundary layers and challenge the standard assumption of efficient angular momentum transport in the inner disk regions. This suggests that the detailed structure of turbulent MHD accretion disk boundary layers could differ appreciably from those derived within the standard framework of turbulent shear viscosity.

The magnetorotational instability (MRI) is considered a key process for driving efficient angular momentum transport in astrophysical disks. Understanding its non-linear saturation constitutes a fundamental problem in modern accretion disk theory. The large dynamical range in physical conditions in accretion disks makes it challenging to address this problem only with numerical simulations. We analyze the concept that (secondary) parasitic instabilities are responsible for the saturation of the MRI. Our approach enables us to explore dissipative regimes that are relevant to astrophysical and laboratory conditions that lie beyond the regime accessible to current numerical simulations. We calculate the spectrum and physical structure of parasitic modes that feed off the fastest, exact (primary) MRI mode when its amplitude is such that the fastest parasitic mode grows as fast as the MRI. We argue that this "saturation" amplitude provides an estimate of the magnetic field that can be generated by the MRI before the secondary instabilities suppress its growth significantly. Recent works suggest that the saturation amplitude of the MRI depends mainly on the magnetic Prandtl number. Our results suggest that, as long as viscous effects do not dominate the fluid dynamics, the saturation level of the MRI depends only on the Elsasser number $\Lambda_\eta$. We calculate the ratio between the stress and the magnetic energy density, $\alpha_{\rm sat}\beta_{\rm sat}$, associated with the primary MRI mode. We find that for $\Lambda_\eta >1$ Kelvin-Helmholtz modes are responsible for saturation and $\alpha_{\rm sat}\beta_{\rm sat} = 0.4$, while for $\Lambda_\eta < 1$ tearing modes prevail and $\alpha_{\rm sat}\beta_{\rm sat} \simeq 0.5 \, \Lambda_\eta$. Several features of MRI simulations in accretion disks surrounding young stars and compact objects can be interpreted in terms of our findings.

We analyze the stability of a dilute plasma with thermal and composition gradients in the limit where conduction is slow compared to the dynamical timescale. We find necessary and sufficient conditions for stability when the background magnetic field is either parallel or perpendicular to the thermal and composition gradients that are parallel to the gravitational field. We provide approximate solutions for all the relevant modes involved, which are driven by gravity, conduction, and diffusion. We discuss the astrophysical implications of our findings for a representative galaxy cluster where helium has sedimented.

We examine the effects of density stratification on magnetohydrodynamic turbulence driven by the magnetorotational instability in local simulations that adopt the shearing box approximation. Our primary result is that, even in the absence of explicit dissipation, the addition of vertical gravity leads to convergence in the turbulent energy densities and stresses as the resolution increases, contrary to results for zero net flux, unstratified boxes. The ratio of total stress to midplane pressure has a mean of ~0.01, although there can be significant fluctuations on long (>~50 orbit) timescales. We find that the time averaged stresses are largely insensitive to both the radial or vertical aspect ratio of our simulation domain. For simulations with explicit dissipation, we find that stratification extends the range of Reynolds and magnetic Prandtl numbers for which turbulence is sustained. Confirming the results of previous studies, we find oscillations in the large scale toroidal field with periods of ~10 orbits and describe the dynamo process that underlies these cycles.

Non-thermal X-ray emission in compact accretion engines can be interpreted to result from magnetic dissipation in an optically thin magnetized corona above an optically thick accretion disk. If coronal magnetic field originates in the disk and the disk is turbulent, then only magnetic structures large enough for their turbulent shredding time to exceed their buoyant rise time survive the journey to the corona. We use this concept and a physical model to constrain the minimum fraction of magnetic energy above the critical scale for buoyancy as a function of the observed coronal to bolometric emission. Our results suggest that a significant fraction of the magnetic energy in accretion disks resides in large scale fields, which in turn provides circumstantial evidence for significant non-local transport phenomena and the need for large scale magnetic field generation. For the example of Seyfert AGN, for which of order 30 per cent of the bolometric flux is in the X-ray band, we find that more than 20 per cent of the magnetic energy must be of large enough scale to rise and dissipate in the corona.

We investigate the stability of incompressible, exact, non-ideal magnetorotational (MRI) modes against parasitic instabilities. Both Kelvin-Helmholtz and tearing-mode parasitic instabilities may occur in the dissipative regimes accessible to current numerical simulations. We suppose that a primary MRI mode saturates at an amplitude such that its fastest parasite has a growth rate comparable to its own. The predicted alpha parameter then depends critically on whether the fastest primary and parasitic modes fit within the computational domain and whether non-axisymmetric parasitic modes are allowed. Hence even simulations that resolve viscous and resistive scales may not saturate properly unless the numerical domain is large enough to allow the free evolution of both MRI and parasitic modes. To minimally satisfy these requirements in simulations with vertical background fields, the vertical extent of the domain should accommodate the fastest growing MRI mode while the radial and azimuthal extents must be twice as large. The fastest parasites have horizontal wavelengths roughly twice as long as the vertical wavelength of the primary.

Numerical simulations of turbulent, magnetized, differentially rotating flows driven by the magnetorotational instability are often used to calculate the effective values of alpha viscosity that is invoked in analytical models of accretion discs. In this paper we use various dynamical models of turbulent magnetohydrodynamic stresses, as well as numerical simulations of shearing boxes, to show that angular momentum transport in MRI-driven accretion discs cannot be described by the standard model for shear viscosity. In particular, we demonstrate that turbulent magnetorotational stresses are not linearly proportional to the local shear and vanish identically for angular velocity profiles that increase outwards.

We carry out a comprehensive analysis of the behavior of the magnetorotational instability (MRI) in viscous, resistive plasmas. We find exact, non-linear solutions of the non-ideal magnetohydrodynamic (MHD) equations describing the local dynamics of an incompressible, differentially rotating background threaded by a vertical magnetic field when disturbances with wavenumbers perpendicular to the shear are considered. We provide a geometrical description of these viscous, resistive MRI modes and show how their physical structure is modified as a function of the Reynolds and magnetic Reynolds numbers. We demonstrate that when finite dissipative effects are considered, velocity and magnetic field disturbances are no longer orthogonal (as it is the case in the ideal MHD limit) unless the magnetic Prandtl number is unity. We generalize previous results found in the ideal limit and show that a series of key properties of the mean Reynolds and Maxwell stresses also hold for the viscous, resistive MRI. In particular, we show that the Reynolds stress is always positive and the Maxwell stress is always negative. Therefore, even in the presence of viscosity and resistivity, the total mean angular momentum transport is always directed outwards. We also find that, for any combination of the Reynolds and magnetic Reynolds numbers, magnetic disturbances dominate both the energetics and the transport of angular momentum and that the total mean energy density is an upper bound for the total mean stress responsible for angular momentum transport. The ratios between the Maxwell and Reynolds stresses and between magnetic and kinetic energy densities increase with decreasing Reynolds numbers for any magnetic Reynolds number; the lowest limit of both ratios is reached in the ideal MHD regime.

We present a scaling law that predicts the values of the stresses obtained in numerical simulations of saturated MRI-driven turbulence in non-stratified shearing boxes. It relates the turbulent stresses to the strength of the vertical magnetic field, the sound speed, the vertical size of the box, and the numerical resolution and predicts accurately the results of 35 numerical simulations performed for a wide variety of physical conditions. We use our result to show that the saturated stresses in simulations with zero net magnetic flux depend linearly on the numerical resolution and would become negligible if the resolution were set equal to the natural dissipation scale in astrophysical disks. We conclude that, in order for MRI-driven turbulent angular momentum transport to be able to account for the large value of the effective alpha viscosity inferred observationally, the disk must be threaded by a significant vertical magnetic field and the turbulent magnetic energy must be in near equipartition with the thermal energy. This result has important implications for the spectra of accretion disks and their stability.

We develop a local model for the exponential growth and saturation of the Reynolds and Maxwell stresses in turbulent flows driven by the magnetorotational instability. We first derive equations that describe the effects of the instability on the growth and pumping of the stresses. We highlight the relevance of a new type of correlations that couples the dynamical evolution of the Reynolds and Maxwell stresses and plays a key role in developing and sustaining the magnetorotational turbulence. We then supplement these equations with a phenomenological description of the triple correlations that lead to a saturated turbulent state. We show that the steady-state limit of the model describes successfully the correlations among stresses found in numerical simulations of shearing boxes.

Recent X-ray variability studies suggest that the log of the square of the fractional rms variability amplitude, rms^2, seems to correlate with the log of the AGN black-hole mass, M_BH, with larger black holes being less variable for a fixed time interval. This has motivated the theoretical modeling of the rms^2-M_BH correlation with the aim of constraining AGN masses based on X-ray variability. A viable approach to addressing this problem is to assume an underlying power spectral density with a suitable mass dependence, derive the functional form of the rms^2-M_BH correlation for a given sampling pattern, and investigate whether the result is consistent with the observations. For simplicity, previous studies, inspired by the similarities shared by the timing properties of AGN and X-ray binaries, have explored model power spectral densities characterized by broken power laws. and ignored, in general, the distorting effects that the particular sampling pattern imprints in the observed power spectral density. Motivated by the latest timing results from X-ray binaries, obtained with RXTE, we propose that AGN broad-band noise spectra consist of a small number of Lorentzian components. This assumption allows, for the first time, to fully account for sampling effects in theoretical models of X-ray variability in an analytic manner. We show that, neglecting sampling effects when deriving the fractional rms from the model power spectral density can lead to underestimating it by a factor of up to 80% with respect to its true value for the typical sampling patterns used to monitor AGN. We discuss the implications of our results for the derivation of AGN masses using theoretical models of the rms^2-M_BH correlation. (Abridged)

The magnetorotational instability is thought to be responsible for the generation of magnetohydrodynamic turbulence that leads to enhanced outward angular momentum transport in accretion discs. Here, we present the first formal analytical proof showing that, during the exponential growth of the instability, the mean (averaged over the disc scale-height) Reynolds stress is always positive, the mean Maxwell stress is always negative, and hence the mean total stress is positive and leads to a net outward flux of angular momentum. More importantly, we show that the ratio of the Maxwell to the Reynolds stresses during the late times of the exponential growth of the instability is determined only by the local shear and does not depend on the initial spectrum of perturbations or the strength of the seed magnetic. Even though we derived these properties of the stress tensors for the exponential growth of the instability in incompressible flows, numerical simulations of shearing boxes show that this characteristic is qualitatively preserved under more general conditions, even during the saturated turbulent state generated by the instability.

During the last decade it has become evident that the magnetorotational instability is at the heart of the enhanced angular momentum transport in weakly magnetized accretion disks around neutron stars and black holes. In this paper, we investigate the local linear stability of differentially rotating, magnetized flows and the evolution of the magnetorotational instability beyond the weak-field limit. We show that, when superthermal toroidal fields are considered, the effects of both compressibility and magnetic tension forces, which are related to the curvature of toroidal field lines, should be taken fully into account. We demonstrate that the presence of a strong toroidal component in the magnetic field plays a non-trivial role. When strong fields are considered, the strength of the toroidal magnetic field not only modifies the growth rates of the unstable modes but also determines which modes are subject to instabilities. We find that, for rotating configurations with Keplerian laws, the magnetorotational instability is stabilized at low wavenumbers for toroidal Alfven speeds exceeding the geometric mean of the sound speed and the rotational speed. We discuss the significance of our findings for the stability of cold, magnetically dominated, rotating fluids and argue that, for these systems, the curvature of toroidal field lines cannot be neglected even when short wavelength perturbations are considered. We also comment on the implications of our results for the validity of shearing box simulations in which superthermal toroidal fields are generated.

We present a detailed study of the effects of gravitational microlensing on compact and distant $\gamma$-ray blazars. These objects have $\gamma$-ray emitting regions which are small enough as to be affected by microlensing effects produced by stars lying in intermediate galaxies. We analyze the temporal evolution of the gamma-ray magnification for sources moving in a caustic pattern field, where the combined effects of thousands of stars are taken into account using a numerical technique. We propose that some of the unidentified $\gamma$-ray sources (particularly some of those lying at high galactic latitude whose gamma-ray statistical properties are very similar to detected $\gamma$-ray blazars) are indeed the result of gravitational lensing magnification of background undetected Active Galactic Nuclei (AGNs).

We present a Monte Carlo analysis of the recently introduced variability indices $\tau$ (Tompkins 1999) and $I$ (Zhang et al. 2000 & Torres et al. 2001) for $\gamma$-ray sources. We explore different variability criteria and prove that these two indices, despite the very different approaches used to compute them, are statistically correlated (5 to 7$\sigma$). This conclusion is maintained also for the subset of AGNs and high latitude ($|b|>10$ deg) sources, whereas the correlation is lowered for the low latitude ones, where the influence of the diffuse galactic emission background is strong.

We analyze how the thermal history of the universe is influenced by the statistical description, assuming a deviation from the usual Bose-Einstein, Fermi-Dirac and Boltzmann-Gibbs distribution functions. These deviations represent the possible appearance of non-extensive effects related with the existence of long range forces, memory effects, or evolution in fractal or multi-fractal space. In the early universe, it is usually assumed that the distribution functions are the standard ones. Then, considering the evolution in a larger theoretical framework will allow to test this assumption and to place limits to the range of its validity. The corrections obtained will change with temperature, and consequently, the bounds on the possible amount of non-extensivity will also change with time. We generalize results which can be used in other contexts as well, as the Boltzmann equation and the Saha law, and provide an estimate on how known cosmological bounds on the masses of neutrinos are modified by a change in the statistics. We particularly analyze here the recombination epoch, making explicit use of the chemical potentials involved in order to attain the necessary corrections. All these results constitute the basic tools needed for placing bounds on the amount of non-extensivity that could be present at different eras and will be later used to study primordial nucleosynthesis.

If the universe is slightly non-extensive, and the distribution functions are not exactly given by those of Boltzmann-Gibbs, the primordial production of light elements will be non-trivially modified. In particular, the principle of detailed balance (PDB), of fundamental importance in the standard analytical analysis, is no longer valid, and a non-extensive correction appears. This correction is computed and its influence is studied and compared with previous works, where, even when the universe was considered as an slightly non-extensive system, the PDB was assumed valid. We analytically track the formation of Helium and Deuterium, and study the kind of deviation one could expect from the standard regime. The correction to the capture time, the moment in which Deuterium can no longer be substantially photo-disintegrated, is also presented. This allows us to take into account the process of the free decay of neutrons, which was absent in all previous treatments of the topic. We show that even when considering a first (linear) order correction in the quantum distribution functions, the final output on the primordial nucleosynthesis yields cannot be reduced to a linear correction in the abundances. We finally obtain new bounds upon the non-extensive parameter, both comparing the range of physical viability of the theory, and using the latest observational data.