Search

Space weather predictions are necessary to avoid damage caused by intense geomagnetic storms. Such strong storms are usually caused by a co-rotating interaction region (CIR) passing at Earth or by the arrival of strong coronal mass ejections (CMEs). To mitigate the damage, the effect of propagating CMEs in the solar wind must be estimated accurately at Earth and other locations. Modelling solar wind accurately is crucial for space weather predictions, as it is the medium for CME propagation. The Icarus heliospheric modelling tool is upgraded to handle dynamic inner heliospheric driving instead of using steady boundary conditions. The ideal magnetohydrodynamic (MHD) solver and the automated grid-adaptivity are adjusted to the latest MPI-AMRVAC version. The inner boundary conditions, prescribed at 0.1 AU for the heliospheric model, are updated time-dependently throughout the simulation. The coronal model is computed repeatedly for selected magnetograms. The particle sampling within MPI-AMRVAC is extended to handle stretched spherical grid information. The solar wind obtained in the simulation is dynamic and shows significant variations throughout the evolution. When comparing the results with the observations, the dynamic solar wind results are more accurate than previous results obtained with purely steady boundary driving. The CMEs propagated through the dynamic solar wind background produce more similar signatures in the time-series data than in the steady solar wind. Dynamic boundary driving in Icarus results in a more self-consistent solar wind evolution in the inner heliosphere. The obtained space weather modelling tool for dynamic solar wind and CME simulations is better suited for space weather forecasting than a steady solar wind model., Comment: 10 pages, 8 figures

We investigate whether Alfv\'enic fluctuations (AFs) can affect the structure of magnetic ejecta (MEs) within interplanetary coronal mass ejections (ICMEs). We study an ICME observed on 2001 December 29 at 1 au by ACE and Wind, at a total angular separation of $\sim$0.8$^\circ$ ($\sim0.014$~au). We focus on the correlation and complexity of its magnetic structure measured between two spacecraft in association with large-amplitude AFs. The Alfv\'enicity of the ME is investigated in terms of the residual energy and cross helicity of fluctuations. We find that as for the event of interest, large-amplitude AFs occur in the rear region of the ME at both Wind and ACE with a duration of about six hours. We compare the correlation of the magnetic field strength and vector components measured between Wind and ACE, and investigate complexity in terms of the magnetic hodograms. The region showing AFs is found to be associated with a decreased correlation of the magnetic field components and an increased complexity of the ME magnetic configuration detected at ACE and Wind, which may be due to the fact that the two spacecraft crossing the same ME along different trajectories likely sampled AFs in different oscillation phases. Combining multi-point in-situ measurements and remote-sensing observations of the ICME source region, we further discuss different potential sources of the AFs., Comment: Accepted by ApJ

The microphysics of the intracluster medium (ICM) in galaxy clusters is still poorly understood. Observational evidence suggests that the effective viscosity is suppressed by plasma instabilities that reduce the mean free path of particles. Measuring the effective viscosity of the ICM is crucial to understanding the processes that govern its physics on small scales. The trails of ionized interstellar medium left behind by the so-called jellyfish galaxies can trace the turbulent motions of the surrounding ICM and constrain its local viscosity. We present the results of a systematic analysis of the velocity structure function (VSF) of the H$\alpha$ line for ten galaxies from the GASP sample. The VSFs show a sub-linear power law scaling below 10 kpc which may result from turbulent cascading and extends to 1 kpc, below the supposed ICM dissipation scales of tens of kpc expected in a fluid described by Coulomb collisions. Our result constrains the local ICM viscosity to be 0.3-25$\%$ of the expected Spitzer value. Our findings demonstrate that either the ICM particles have a smaller mean free path than expected in a regime defined by Coulomb collisions, or that we are probing effects due to collisionless physics in the ICM turbulence., Comment: 24 pages, 14 figures, accepted for publication on ApJ on 11/11/2024

We propose novel numerical schemes based on the Boris method in curved spacetime, incorporating both hadronic and radiative interactions for the first time. Once the proton has lost significant energy due to radiative and hadronic losses, and its gyroradius has decreased below typical scales on which the electromagnetic field varies, we apply a guiding center approximation (GCA). We fundamentally simulate collision processes either with a Monte-Carlo method or, where applicable, as a continuous energy loss, contingent on the local optical depth. To test our algorithm for the first time combining the effects of electromagnetic, gravitational, and radiation fields including hadronic interactions, we simulate highly relativistic protons traveling through various electromagnetic fields and proton backgrounds. We provide unit tests in various spatially dependent electromagnetic and gravitational fields and background photon and proton distributions, comparing the trajectory against analytic results. We propose that our method can be used to analyze hadronic interactions in black hole accretion disks, jets, and coronae to study the neutrino abundance from active galactic nuclei., Comment: 15 pages, 9 figures, submitted to Physical Review D

We present the largest 3D Particle-in-Cell shearing-box simulations of turbulence driven by the magnetorotational instability, for the first time employing the realistic proton-to-electron mass ratio. We investigate the energy partition between relativistically hot electrons and subrelativistic ions in turbulent accretion flows, a regime relevant to collisionless, radiatively inefficient accretion flows of supermassive black holes, such as those targeted by the Event Horizon Telescope. We provide a simple empirical formula to describe the measured heating ratio between ions and electrons, which can be used for more accurate global modeling of accretion flows with standard fluid approaches such as general-relativistic magnetohydrodynamics.

The gas in the interstellar medium (ISM) of galaxies is supersonically turbulent. Measurements of turbulence typically rely on cold gas emission lines for low-z galaxies and warm ionized gas observations for z>0 galaxies. Studies of warm gas kinematics at z>0 conclude that the turbulence strongly evolves as a function of redshift, due to the increasing impact of gas accretion and mergers in the early Universe. However, recent findings suggest potential biases in turbulence measurements derived from ionized gas at high-z, impacting our understanding of turbulence origin, ISM physics and disk formation. We investigate the evolution of turbulence using velocity dispersion ($\sigma$) measurements from cold gas tracers (i.e., CO, [CI], [CII]) derived from a sample of 57 galaxy disks spanning the redshift range z=0-5. This sample consists of main-sequence and starburst galaxies with stellar masses $\gtrsim 10^{10} M_{\odot}$. The comparison with current H$\alpha$ kinematic observations and existing models demonstrates that the velocity dispersion inferred from cold gas tracers differ by a factor of $\approx 3$ from those obtained using emission lines tracing warm gas. We show that stellar feedback is the main driver of turbulence measured from cold gas tracers. This is fundamentally different from the conclusions of studies based on warm gas, which had to consider additional turbulence drivers to explain the high values of $\sigma$. We present a model predicting the redshift evolution of turbulence in galaxy disks, attributing the increase of $\sigma$ with redshift to the higher energy injected by supernovae due to the elevated star-formation rate in high-z galaxies. This supernova-driven model suggests that turbulence is lower in galaxies with lower stellar mass compared to those with higher stellar mass. Additionally, it forecasts the evolution of $\sigma$ in Milky-Way like progenitors., Comment: Accepted for publication in A&A. The abstract has been modified to comply with arXiv's character limit

The local gravitational instability of rotating discs is believed to be an important mechanism in different astrophysical processes, including the formation of gas and stellar clumps in galaxies. We aim to study in three dimensions the local gravitational instability of two-component thick discs. We take as starting point a recently proposed analytic three-dimensional (3D) instability criterion for discs with non-negligible thickness which takes the form $Q_{\rm 3D}<1$, where $Q_{\rm 3D}$ is a 3D version of the classical 2D Toomre $Q$ parameter for razor-thin discs. Here we extend the 3D stability analysis to two-component discs, considering first the influence on $Q_{\rm 3D}$ of a second unresponsive component, and then the case in which both components are responsive. We present the application to two-component discs with isothermal vertical distributions, which can represent, for instance, galactic discs with both stellar and gaseous components. Finally, we relax the assumption of vertical isothermal distribution, by studying one-component self-gravitating discs with polytropic vertical distributions for a range of values of the polytropic index corresponding to convectively stable configurations. We find that $Q_{\rm 3D}<1$, where $Q_{\rm 3D}$ can be computed from observationally inferred quantities, is a robust indicator of local gravitational instability, depending only weakly on the presence of a second component and on the vertical gradient of temperature or velocity dispersion. We derive a sufficient condition for local gravitational instability in the midplane of two-component discs, which can be employed when both components have $Q_{\rm 3D}>1$., Comment: 14 pages, 9 figures, accepted for publication in A&A

Turbulent states are ubiquitous in plasmas and the understanding of turbulence is fundamental in modern astrophysics. Numerical simulations, which are the state-of-the-art approach to the study of turbulence, require substantial computing resources. Recently, attention shifted to methods for generating synthetic turbulent magnetic fields, affordably creating fields with parameter-controlled characteristic features of turbulence. In this context, the BxC toolkit was developed and validated against direct numerical simulations (DNS) of isotropic turbulent magnetic fields. Here, we demonstrate novel extensions of BxC to generate realistic turbulent magnetic fields in a fast, controlled, geometric approach. First, we perform a parameter study to determine quantitative relations between the BxC input parameters and desired characteristic features of the turbulent power spectrum, such as the extent of the inertial range, its spectral slope, and the injection and dissipation scale. Second, we introduce in the model a set of structured background magnetic fields B0, as a natural and more realistic extension to the purely isotropic turbulent fields. Third, we extend the model to include anisotropic turbulence properties in the generated fields. With all these extensions combined, our tool can quickly generate any desired structured magnetic field with controlled, anisotropic turbulent fluctuations, faster by orders of magnitude with respect to DNSs. These can be used, e.g., to provide initial conditions for DNS simulations or easily generate synthetic data for many astrophysical settings, all at otherwise unaffordable resolutions.

The expanding solar wind plasma ubiquitously exhibits anisotropic non-thermal particle velocity distributions. Typically, proton Velocity Distribution Functions (VDFs) show the presence of a core and a field-aligned beam. Novel observations made by Parker Solar Probe (PSP) in the innermost heliosphere have revealed new complex features in the proton VDFs, namely anisotropic beams that sometimes experience perpendicular diffusion. In this study, we use a 2.5D fully kinetic simulation to investigate the stability of proton VDFs with anisotropic beams observed by PSP. Our setup consists of a core and an anisotropic beam populations that drift with respect to each other. This configuration triggers a proton-beam instability from which nearly parallel fast magnetosonic modes develop. Our results demonstrate that before this instability reaches saturation, the waves resonantly interact with the beam protons, causing perpendicular heating at the expense of the parallel temperature.

Local gravitational instability (LGI) is considered crucial for regulating star formation and gas turbulence in galaxy discs, especially at high redshift. Instability criteria usually assume infinitesimally thin discs or rely on approximations to include the stabilising effect of the gas disc thickness. We test a new 3D instability criterion for rotating gas discs that are vertically stratified in an external potential. This criterion reads $Q_{\rm3D}<1$, where $Q_{\rm3D}$ is the 3D analogue of the Toomre parameter $Q$. The advantage of $Q_{\rm3D}$ is that it allows us to study LGI in and above the galaxy midplane in a rigorous and self-consistent way. We apply the criterion to a sample of 44 star-forming galaxies at $0\lesssim\mathrm{z}\lesssim5$ hosting rotating discs of cold gas. The sample is representative of galaxies on the main sequence at $\mathrm{z}\approx 0$ and includes massive star-forming and starburst galaxies at $1\lesssim\mathrm{z}\lesssim5$. For each galaxy, we first apply the Toomre criterion for infinitesimally thin discs, finding 10 unstable systems. We then obtain maps of $Q_{\rm 3D}$ from a 3D model of the gas disc derived in the combined potential of dark matter, stars and the gas itself. According to the 3D criterion, two galaxies with $Q<1$ show no evidence of instability and the unstable regions that are 20% smaller than those where $Q<1$. No unstable disc is found at $0\lesssim\mathrm{z}\lesssim 1$, while $\approx 60$% of the systems at $2\lesssim\mathrm{z}\lesssim5$ are locally unstable. In these latter, a relatively small fraction of the total gas ($\approx 30$%) is potentially affected by the instability. Our results disfavour LGI as the main regulator of star formation and turbulence in moderately star-forming galaxies in the present-day Universe. LGI likely becomes important at high redshift, but the input by other mechanisms seems required [abridged], Comment: 19 pages (6 figures, 2 tables) + appendixes, accepted by Astronomy & Astrophysics

We present a study of energetic-electron trapping and acceleration in the Kelvin-Helmholtz-induced magnetohydrodynamic (MHD) turbulence of post-flare loops in the solar corona. Using the particle-tracing capabilities of MPI-AMRVAC 3.0, we evolve ensembles of test electrons (i.e. without feedback to the underlying MHD) inside the turbulent looptop, using the guiding-center approximation. With the MHD looptop model of Ruan et al. 2018, we investigate the relation between turbulence and particle trapping inside the looptop structure, showing that better-developed turbulent cascades result in more efficient trapping primarily due to mirror effects. We then quantify the electron acceleration in the time-evolving MHD turbulence, and find that ideal-MHD processes inside the looptop can produce nonthermal particle spectra from an initial Maxwellian distribution. Electrons in this turbulence are preferentially accelerated by mirror effects in the direction perpendicular to the local magnetic field while remaining confined within small regions of space between magnetic islands. Assuming dominance of Bremsstrahlung radiation mechanisms, we employ the resulting information from accelerated electrons (combined with the MHD background) to construct HXR spectra of the post-flare loop that include nonthermal-particle contributions. Our results pave the way to constructing more realistic simulations of radiative coronal structure for comparison with current and future observations., Comment: Accepted in MNRAS

Electromagnetic fundamental and harmonic emission is ubiquitously observed throughout the heliosphere, and in particular it is commonly associated with the occurrence of Type II and III solar radio bursts. Classical analytic calculations for the plasma-emission process, though useful, are limited to idealized situations; a conclusive numerical verification of this theory is still lacking, with earlier studies often providing contradicting results on e.g. the precise parameter space in which fundamental and harmonic emission can be produced. To accurately capture the chain of mechanisms underlying plasma emission - from precursor plasma processes to the generation of electromagnetic waves over long times - we perform large-scale, first-principles simulations of beam-plasma instabilities. By employing a very large number of computational particles we achieve very low numerical noise, and explore (with an array of simulations) a wide parameter space determined by the beam-plasma density ratio and the ion-to-electron temperature ratio. In particular, we observe direct evidence of both fundamental and harmonic plasma emission when the beam-to-background density ratio $\le$0.005 (with beam-to-background energy ratio ~0.5), tightly constraining this threshold. We observe that, asymptotically, in this regime ~0.1% of the initial beam energy is converted into harmonic emission, and ~0.001% into fundamental emission. In contrast with previous studies, we also find that this emission is independent of the ion-to-electron temperature ratio. In addition, we report the direct detection of third-harmonic emission in all of our simulations, at power levels compatible with observations. Our findings have important consequences for understanding the viable conditions leading to plasma emission in space systems, and for the interpretation of observed electromagnetic signals throughout the heliosphere., Comment: Accepted in MNRAS

We present 3D fully kinetic shearing-box simulations of pair-plasma magnetorotational turbulence with unprecedented macro-to-microscopic scale separation. While retrieving the expected fluid behavior of the plasma at large scales, we observe a steepening of turbulent spectra at kinetic scales and substantial angular-momentum transport linked with kinetic processes. For the first time, we provide a definitive demonstration of nonthermal particle acceleration in kinetic magnetorotational turbulence agnostically of initial conditions, by means of a novel strategy exploiting synchrotron cooling., Comment: Accepted in PRL

This study presents three-dimensional (3D) resistive Hall-magnetohydrodynamic simulations of the Kelvin-Helmholtz instability (KHI) dynamics at Earth's magnetospheric flanks during northward interplanetary magnetic field periods. By comparing two simulations with and without initial magnetic shear, we analyze the impact of distinct magnetic field orientations on plasma dynamics and magnetic reconnection events taking into account 3D mechanisms, such as KHI high latitude stabilization. The identical nature of the simulations, except for the presence/absence of an initial magnetic shear, enables, for the first time, a complete and coherent comparative analysis of the latitudinal distribution of KH vortices, current sheets, reconnection events, and the evolution of the mixing layer. In one configuration, a uniform magnetic field leads to double mid-latitude reconnection (MLR), while in the other, magnetic shear induces both type I vortex-induced reconnection (VIR) and MLR. Notably, the type I VIR observed in this second scenario results from the combined action of line advection and vortex-induced current sheet pinching (the classic mechanism driving two-dimensional type I VIR). Of particular importance is our quantification of newly closed field lines that experienced double reconnection, ultimately becoming embedded in solar wind plasma at low latitudes while remaining connected to magnetospheric plasma at high latitudes. The varying abundance of such lines in the two simulations holds implications for plasma transport at the magnetopause.

Using images collected with the WFC3 camera on board of the Hubble Space Telescope, we detect stellar clumps in continuum-subtracted $H\alpha$ and ultraviolet (F275W filter), such clumps are often embedded in larger regions (star-forming complexes) detected in the optical (F606W filter). We model the photometry of these objects using BAGPIPES to obtain their stellar population parameters. The median mass-weighted stellar ages are 27 Myr for $H\alpha$ clumps and 39 Myr for F275W clumps and star-forming complexes, the oldest stars in the complexes can be older than $\sim$300 Myr which indicates that star-formation is sustained for long periods of time. Stellar masses vary from 10$^{3.5}$ to 10$^{7.1}$ $M_\odot$, with star-forming complexes being more massive objects in the sample. Clumps and complexes found further away from the host galaxy are younger, less massive and less obscured by dust. We interpret these trends as due to the effect of ram-pressure in different phases of the interstellar medium. $H\alpha$ clumps form a well-defined sequence in the stellar mass--SFR plane with slope 0.73. Some F275W clumps and star-forming complexes follow the same sequence while others stray away from it and passively age. The difference in stellar age between a complex and its youngest embedded clump scales with the distance between the clump and the center of the complex, with the most displaced clumps being hosted by the most elongated complexes. This is consistent with a fireball-like morphology, where star-formation proceeds in a small portion of the complex while older stars are left behind producing a linear stellar population gradient. The stellar masses of star-forming complexes are consistent with the ones of globular clusters, but stellar mass surface densities are lower by 2 dex, and their properties are more consistent with the population of dwarf galaxies in clusters., Comment: Accepted for publication in A&A

Interplanetary coronal mass ejections (ICMEs) are defined as ''coherent'' if they are capable of responding to external perturbations in a collective manner. This implies that information must be able to propagate across ICME structures, and if this is not the case, single-point in-situ measurements cannot be considered as indicative of global ICME properties. Here, we investigate the role of Alfvenic fluctuations (AFs) as mediators of ICME coherence. We consider multi-point magnetic field and plasma measurements of 10 ICMEs observed by the ACE and Wind spacecraft at 1 au at longitudinal separations of 0.5{\deg}-0.7{\deg}. For each event, we analyze the Alfvenicity in terms of the residual energy and cross helicity of fluctuations, and the coherence in terms of the magnetic correlation between Wind and ACE. We find that ~65% and 90% of ICME sheaths and magnetic ejecta (MEs), respectively, present extended AFs covering at least 20% of the structure. Cross helicity suggests AFs of solar and interplanetary origin may co-exist in the ICME population at 1 au. AFs are mainly concentrated downstream of shocks and in the back of MEs. The magnetic field is poorly correlated within sheaths, while the correlation decreases from the front to the back of the MEs for most magnetic field components. AFs are also associated with lower magnetic field correlations. This suggests either that ICME coherence is not mediated by Alfven waves, implying that the coherence scale may be smaller than previously predicted, or that the magnetic field correlation is not a measure of coherence.

Ram pressure stripping (RPS) is one of the most invoked mechanisms to explain the observed differences between cluster and field galaxies. In the local Universe, its effect on the galaxy star forming properties has been largely tackled and the general consensus is that this process first compresses the gas available in the galaxy disks, boosting the star formation for a limited amount of time, and then removes the remaining gas leading to quenching. Much less is known on the effect and preponderance of RPS at higher redshift, due to the lack of statistical samples. Exploiting VLT/MUSE observations of galaxies at 0.23sigma when considering the SFR occurring in both the tail and disk of galaxies. The enhancement is retrieved also on local scales: considering spatially resolved data, ram pressure stripped galaxies overall have large {\Sigma}SFR values, especially for Sigma_*>10^7.5M_sun kpc~2. RPS seems to leave the same imprint on the SFR-M* and Sigma_SFR-Sigma_* relations both in the Local Universe and at z~0.35., Comment: A&A in press

In this paper, we introduce and discuss an exactly energy-conserving Particle-in-Cell method for arbitrary curvilinear coordinates. The flexibility provided by curvilinear coordinates enables the study of plasmas in complex-shaped domains by aligning the grid to the given geometry, or by focusing grid resolution on regions of interest without overresolving the surrounding, potentially uninteresting domain. We have achieved this through the introduction of the metric tensor, the Jacobian matrix, and contravariant operators combined with an energy-conserving fully implicit solver. We demonstrate the method's capabilities using a Python implementation to study several one- and two-dimensional test cases: the electrostatic two-stream instability, the electromagnetic Weibel instability, and the geomagnetic environment modeling (GEM) reconnection challenge. The test results confirm the capability of our new method to reproduce theoretical expectations (e.g. instability growth rates) and the corresponding results obtained with a Cartesian uniform grid when using curvilinear grids. Simultaneously, we show that the method conserves energy to machine precision in all cases., Comment: 14 pages, 5 figures

We characterize the morphological properties of a statistically relevant sample of H$\alpha$ and UV young star-forming clumps and optical complexes, observed with the \textit{Hubble Space Telescope} in six galaxies of the GASP sample undergoing ram-pressure stripping. The catalogs comprise 2406 (323 in the tails) H$\alpha$ clumps, 3750 (899) UV clumps and 424 tail optical complexes. About 15-20\% of the clumps and 50\% of the complexes are resolved in size. We find that more than half of the complexes contain no H$\alpha$ clumps, while most of them contain at least one UV clump. The clump number and size increase with the complex size, while the median complex filling factor is larger for UV clumps ($0.27$) than for H$\alpha$ clumps ($0.10$) and does not correlate with almost any morphological property. This suggests that the clumps number and size grow with the complex keeping the filling factor constant. When studying the position of the clumps inside their complexes, H$\alpha$ clumps, and UV clumps to a lesser extent, show a displacement from the complex center of $0.1-1$ kpc and, in $\sim 60$\% of the cases, they are displaced away from the galactic disk. This is in accordance with the fireball configuration, already observed in the tails of stripped galaxies. Finally, the filling factor and the clump radius increase with the distance from the galactic disk, suggesting that the reciprocal displacement of the different stellar generations increases as a consequence of the velocity gradient caused by ram pressure., Comment: 20 pages, 14 figures, 2 tables. Accepted for publication in ApJ

Cluster galaxies are affected by the surrounding environment, which influences, in particular, their gas, stellar content and morphology. In particular, the ram-pressure exerted by the intracluster medium promotes the formation of multi-phase tails of stripped gas detectable both at optical wavelengths and in the sub-mm and radio regimes, tracing the cold molecular and atomic gas components, respectively. In this work we analyze a sample of sixteen galaxies belonging to clusters at redshift $\sim 0.05$ showing evidence of an asymmetric HI morphology (based on MeerKAT observations) with and without a star forming tail. To this sample we add three galaxies with evidence of a star forming tail and no HI detection. Here we present the galaxies $\rm H_{2}$ gas content from APEX observations of the CO(2-1) emission. We find that in most galaxies with a star forming tail the $\rm H_{2}$ global content is enhanced with respect to undisturbed field galaxies with similar stellar masses, suggesting an evolutionary path driven by the ram-pressure stripping. As galaxies enter into the clusters their HI is displaced but also partially converted into $\rm H_{2}$, so that they are $\rm H_{2}$ enriched when they pass close to the pericenter, i. e. when they develop the star forming tails that are visible in UV/B broad bands and in H$\alpha$ emission. An inspection of the phase-space diagram for our sample suggests an anticorrelation between the HI and $\rm H_{2}$ gas phases as galaxies fall into the cluster potential. This peculiar behaviour is a key signature of the ram-pressure stripping in action., Comment: 19 pages, 8 figures, accepted for publication in ApJ

The magnetorotational instability (MRI) plays a crucial role in regulating the accretion efficiency in astrophysical accretion disks. In low-luminosity disks around black holes, such as Sgr A* and M87, Coulomb collisions are infrequent, making the MRI physics effectively collisionless. The collisionless MRI gives rise to kinetic plasma effects that can potentially affect its dynamic and thermodynamic properties. We present 2D and 3D particle-in-cell (PIC) plasma simulations of the collisionless MRI in stratified disks using shearing boxes with net vertical field. We use pair plasmas, with initial $\beta=100$ and concentrate on sub-relativistic plasma temperatures ($k_BT \lesssim mc^2$). Our 2D and 3D runs show disk expansion, particle and magnetic field outflows, and a dynamo-like process. They also produce magnetic pressure dominated disks with (Maxwell stress dominated) viscosity parameter $\alpha \sim 0.5-1$. By the end of the simulations, the dynamo-like magnetic field tends to dominate the magnetic energy and the viscosity in the disks. Our 2D and 3D runs produce fairly similar results, and are also consistent with previous 3D MHD simulations. Our simulations also show nonthermal particle acceleration, approximately characterized by power-law tails with temperature dependent spectral indices $-p$. For temperatures $k_BT \sim 0.05-0.3\, mc^2$, we find $p\approx 2.2-1.9$. The maximum accelerated particle energy depends on the scale separation between MHD and Larmor-scale plasma phenomena in a way consistent with previous PIC results of magnetic reconnection-driven acceleration. Our study constitutes a first step towards modeling from first principles potentially observable stratified MRI effects in low-luminosity accretion disks around black holes.

The main component of the NextGeneration EU (NGEU) program is the Recovery and Resilience Facility (RRF), spanning an implementation period between 2021 and 2026. The RRF also includes a monitoring system: every six months, each country is required to send an update on the progress of the plan against 14 common indicators, measured on specific quantitative scales. The aim of this paper is to present the first empirical evidence on this system, while, at the same time, emphasizing the potential of its integration with the sustainable development framework (SDGs). We propose to develop a first linkage between the 14 common indicators and the SDGs which allows us to produce a composite index (SDGs-RRF) for France, Germany, Italy, and Spain for the period 2014-2021. Over this time, widespread improvements in the composite index across the four countries led to a partial reduction of the divergence. The proposed approach represents a first step towards a wider use of the SDGs for the assessment of the RRF, in line with their use in the European Semester documents prepared by the European Commission., Comment: 30 pages, in Italian, 10 figures, 13 tables

We present a novel Relativistic Semi-Implicit Method (RelSIM) for particle-in-cell (PIC) simulations of astrophysical plasmas, implemented in a code framework ready for production runs. While explicit PIC methods have gained widespread recognition in the astrophysical community as a reliable tool to simulate plasma phenomena, implicit methods have been seldom explored. This is partly due to the lack of a reliable relativistic implicit PIC formulation that is applicable to state-of-the-art simulations. We propose the RelSIM to fill this gap: our new method is relatively simple, being free of nonlinear iterations and only requiring a global linear solve of the field equations. With a set of one- and two-dimensional tests, we demonstrate that the RelSIM produces more accurate results with much smaller numerical errors in the total energy than standard explicit PIC, particularly when characteristic plasma scales (skin depth and plasma frequency) are heavily underresolved on the numerical grid. By construction, the RelSIM also performs much better than the Relativistic Implicit-Moment Method (RelIMM), originally proposed for semi-implicit PIC simulations in the relativistic regime. Our results are promising to conduct large-scale (in terms of duration and domain size) PIC simulations of astrophysical plasmas, potentially reaching physical regimes inaccessible by standard explicit PIC codes., Comment: Accepted in ApJS

Computational astrophysics routinely combines grid-adaptive capabilities with modern shock-capturing, high resolution spatio-temporal schemes on multi-dimensional hydro- and magnetohydrodynamics. We provide an update on developments within the open-source MPI-AMRVAC code. With online documentation, the MPI-AMRVAC 3.0 release includes several added equation sets, and many options to explore and quantify the influence of implementation details. Showcasing this on a variety of hydro and MHD tests, we document new modules of interest for state-of-the-art solar applications. Test cases address how higher order reconstructions impact long term simulations of shear layers, with and without gas-dust coupling, how runaway radiative losses transit to intricate multi-temperature, multi-phase dynamics, and how different flavors of spatio-temporal schemes and magnetic monopole control produce consistent MHD results in combination with adaptive meshes. We demonstrate Super-Time-Stepping strategies for specific parabolic terms and give details on all implemented Implicit-Explicit integrators. A new magnetofrictional module can be used for computing force-free magnetic fields or for data-driven time-dependent evolutions, while the Regularized-Biot-Savart-Law approach can insert fluxropes in 3D domains. Synthetic observations of 3D MHD simulations can be rendered on-the-fly, or in post-processing, in many spectral wavebands. A particle module and a generic fieldline tracing, compatible with the hierarchical meshes, can be used to sample information at prescribed locations, to follow dynamics of charged particles, or realize two-way coupled simulations between MHD setups and field-aligned non-thermal processes. Highlighting the latest additions and various technical aspects, our open-source strategy welcomes any further code usage, contribution, or spin-off development., Comment: Accepted for A&A

We investigate whether photon ring observations in black hole imaging are able to distinguish between the Kerr black hole in general relativity and alternative black holes that deviate from Kerr. Certain aspects of photon rings have been argued to be robust observables in Very-Long-Baseline Interferometry (VLBI) black hole observations which carry imprints of the underlying spacetime. The photon ring shape, as well as its Lyapunov exponent (which encodes the narrowing of successive photon subrings), are detailed probes of the underlying geometry; measurements thereof have been argued to provide a strong null test of general relativity and the Kerr metric. However, a more complicated question is whether such observations of the photon ring properties can distinguish between Kerr and alternative black holes. We provide a first answer to this question by calculating photon rings of the Johannsen, Rasheed-Larsen, and Manko-Novikov black holes. We find that large deviations from Kerr and large observer inclinations are needed to obtain measurable differences in the photon ring shape. In other words, the Kerr photon ring shape appears to be the universal shape even for deviating black holes at low inclinations. On the other hand, the Lyapunov exponent shows more marked variations for deviations from the Kerr metric. Our analysis lays out the groundwork to determine deviations from the Kerr spacetime in photon rings that are potentially detectable by future observing missions., Comment: 31 pages + appendices; 20 figures

Cluster galaxies are subject to the ram pressure exerted by the intracluster medium, which can perturb or even strip away their gas while leaving the stars undisturbed. We model the distribution and kinematics of the stars and the molecular gas in four late-type cluster galaxies (JO201, JO204, JO206, and JW100), which show tails of atomic and ionized gas indicative of ongoing ram pressure stripping. We analyze MUSE@VLT data and CO data from ALMA searching for signatures of radial gas flows, ram pressure stripping, and other perturbations. We find that all galaxies, with the possible exception of JW100, host stellar bars. Signatures of ram pressure are found in JO201 and JO206, which also shows clear indications of ongoing stripping in the molecular disk outskirts. The stripping affects the whole molecular gas disk of JW100. The molecular gas kinematics in JO204 is instead dominated by rotation rather than ram pressure. We also find indications of enhanced turbulence of the molecular gas compared to field galaxies. Large-scale radial flows of molecular gas are present in JO204 and JW100, but more uncertain in JO201 and JO206. We show that our sample follows the molecular gas mass-size relation, confirming that it is essentially independent of environment even for the most extreme cases of stripping. Our findings are consistent with the molecular gas being affected by the ram pressure on different timescales and less severely than the atomic and ionized gas phases, likely because the molecular gas is denser and more gravitationally bound to the galaxy., Comment: 29 pages, 14 figures, accepted by ApJ

We study the accretion of collisionless plasma onto a rotating black hole from first principles using axisymmetric general-relativistic particle-in-cell simulations. We carry out a side-by-side comparison of these results to analogous general-relativistic magnetohydrodynamic simulations. Although there are many similarities in the overall flow dynamics, three key differences between the kinetic and fluid simulations are identified. Magnetic reconnection is more efficient, and rapidly accelerates a nonthermal particle population, in our kinetic approach. In addition, the plasma in the kinetic simulations develops significant departures from thermal equilibrium, including pressure anisotropy that excites kinetic-scale instabilities, and a large field-aligned heat flux near the horizon that approaches the free-streaming value. We discuss the implications of our results for modeling event-horizon scale observations of Sgr A* and M87 by GRAVITY and the Event Horizon Telescope., Comment: 10 pages, 8 figures, submitted

HI and CO observations indicate that the cold gas in galaxies is very turbulent. However, the turbulent energy is expected to be quickly dissipated, implying that some energy source is needed to explain the observations. The nature of such turbulence was long unclear, as even the main candidate, supernova (SN) feedback, seemed insufficient. Other mechanisms have been proposed, but without reaching a general consensus. The key novelty of our work is considering that the gas disc thickness and flaring increase the dissipation timescale of turbulence, thus reducing the energy injection rate required to sustain it. In excellent agreement with the theoretical expectations, we found that the fraction of the SN energy (a.k.a. SN coupling efficiency) needed to maintain the cold gas turbulence is $\sim 1$%, solving a long-standing conundrum., Comment: 3 pages, 1 figure, Proceedings IAU Symposium 373, accepted

Star formation laws are empirical relations between the cold gas (HI+H$_2$) content of a galaxy and its star formation rate (SFR), being crucial for any model of galaxy formation and evolution. A well known example of such laws is the Schmidt-Kennicutt law, which is based on the projected surface densities. However, it has been long unclear whether a more fundamental relation exists between the intrinsic volume densities. By assuming the vertical hydrostatic equilibrium, we infer radial profiles for the thickness of gaseous discs in a sample of 23 local galaxies, and use these measurements to convert the observed surface densities of the gas and the SFR into the de-projected volume densities. We find a tight correlation linking these quantities, that we call the volumetric star formation law. This relation and its properties have crucial implications for our understanding of the physics of star formation., Comment: 3 pages, 1 figure, Proceedings IAU Symposium 373, accepted

Magnetic reconnection can power spectacular high-energy astrophysical phenomena by producing non-thermal energy distributions in highly magnetized regions around compact objects. By means of two-dimensional fully kinetic particle-in-cell (PIC) simulations we investigate relativistic collisionless plasmoid-mediated reconnection in magnetically dominated pair plasmas with and without guide field. In X-points, where diverging flows result in a non-diagonal thermal pressure tensor, a finite residence time for particles gives rise to a localized collisionless effective resistivity. Here, for the first time for relativistic reconnection in a fully developed plasmoid chain we identify the mechanisms driving the non-ideal electric field using a full Ohm's law by means of a statistical analysis based on our PIC simulations. We show that the non-ideal electric field is predominantly driven by gradients of nongyrotropic thermal pressures. We propose a kinetic physics motivated non-uniform effective resistivity model, which is negligible on global scales and becomes significant only locally in X-points, that captures the properties of collisionless reconnection with the aim of mimicking its essentials in non-ideal magnetohydrodynamic descriptions. This effective resistivity model provides a viable opportunity to design physically grounded global models for reconnection-powered high-energy emission., Comment: 12 pages, 4 figures

We present MeerKAT HI observations of six jellyfish candidate galaxies (JFCGs) in the galaxy cluster, A2626. Two of the six galaxies JW100 and JW103, that were identified as JFCGs from B-band images, are confirmed as jellyfish galaxies (JFGs). Both of the JFGs have low HI content, reside in the cluster core, and move at very high velocities ($\sim$ 3$\sigma_{cl}$). The other JFCGs, identified as non-jellyfish galaxies, are HI rich, with HI morphologies revealing warps, asymmetries, and possible tidal interactions. Both the A2626 JFGs and three other confirmed JFGs from the GASP sample show that these galaxies are HI stripped but not yet quenched. We detect HI, Halpha, and CO tails of similar extent ($\sim$ 50 kpc) in JW100. Comparing the multi-phase velocity channels, we do not detect any HI or CO emission in the northern section of the tail where Halpha emission is present, possibly due to prolonged interaction between the stripped gas and the ICM. We also observe an anti-correlation between HI and CO, which hints at an efficient conversion of HI to H2 in the southern part of the tail. We find that both RPS and HI-to-H2 conversion are significant depletion channels for atomic gas. HI-to-H2 conversion is more efficient in the disc than in the tail., Comment: 15 pages, 10 figures, accepted for publication in MNRAS

We present a study of the orbits, environments and morphologies of 13 ram-pressure stripped galaxies in the massive, intermediate redshift (z$\sim0.3-0.4$) galaxy clusters A2744 and A370, using MUSE integral-field spectroscopy and HST imaging from the Frontier Fields Program. We compare different measures of the locations and morphologies of the stripped sample with a sample of 6 poststarburst galaxies identified within the same clusters, as well as the general cluster population. We calculate the phase space locations of all cluster galaxies and carry out a substructure analysis, finding that the ram-pressure stripped galaxies in A370 are not associated with any substructures, but are likely isolated infalling galaxies. In contrast, the ram-pressure stripped galaxies in A2744 are strictly located within a high-velocity substructure, moving through a region of dense X-ray emitting gas. We conclude that their ram-pressure interactions are likely to be the direct result of the merger between two components of the cluster. Finally, we study the morphologies of the stripped and poststarburst galaxies, using numerical measures to quantify the level of visual disturbances. We explore any morphological deviations of these galaxies from the cluster population, particularly the weaker cases which have been confirmed via the presence of ionised gas tails to be undergoing ram-pressure stripping, but are not strongly visually disturbed in the broad-band data. We find that the stripped sample galaxies are generally divergent from the general cluster sample, with poststarburst galaxies being intermediary in morphology between stripped galaxies and red passive cluster members., Comment: 27 pages, 17 figures

In 1980 Kowal and Drake found that in December 1612 and January 1613 Galileo observed the planet Neptune. At that time, according to these authors, Galileo was able to measure angular separations with an accuracy of about 10 seconds of arc. However, as noticed by Kowal and Drake, the position of Neptune reported by Galileo is wrong with respect to the position computed with the modern ephemeris of about 1 minute of arc. This led Kowal and Drake to speculate on the possible errors of modern ephemeris of Neptune and sparked some debate about Neptune's ephemeris and/or possible errors in Galileo's measures. Until today this anomaly has remained without a conclusive answer. Here we show that, in addition to the random errors, there are other significant measurement errors present in Galileo's observations. These errors may help clarify the origin of the alleged anomalies in the position of Neptune.

The magnetorotational instability (MRI) is a fundamental mechanism determining the macroscopic dynamics of astrophysical accretion disks. In collisionless accretion flows around supermassive black holes, MRI-driven plasma turbulence cascading to microscopic (i.e. kinetic) scales can result in enhanced angular-momentum transport and redistribution, nonthermal particle acceleration, and a two-temperature state where electrons and ions are heated unequally. However, this microscopic physics cannot be captured with standard magnetohydrodynamic (MHD) approaches typically employed to study the MRI. In this work, we explore the nonlinear development of MRI turbulence in a pair plasma, employing fully kinetic Particle-in-Cell (PIC) simulations in two and three dimensions. First, we thoroughly study the axisymmetric MRI with 2D simulations, explaining how and why the 2D geometry produces results that differ substantially from MHD expectations. We then perform the largest (to date) 3D simulations, for which we employ a novel shearing-box approach, demonstrating that 3D PIC models can reproduce the mesoscale (i.e. MHD) MRI dynamics in sufficiently large runs. With our fully kinetic simulations, we are able to describe the nonthermal particle acceleration and angular-momentum transport driven by the collisionless MRI. Since these microscopic processes ultimately lead to the emission of potentially measurable radiation in accreting plasmas, our work is of prime importance to understand current and future observations from first principles, beyond the limitations imposed by fluid (MHD) models. While in this first study we focus on pair plasmas for simplicity, our results represent an essential step toward designing more realistic electron-ion simulations, on which we will focus in future work., Comment: Submitted to ApJ

Gas discs of late-type galaxies are flared, with scale heights increasing with the distance from the galaxy centres and often reaching kpc scales. We study the effects of gas disc flaring on the recovered dark matter halo parameters from rotation curve decomposition. For this, we carefully select a sample of 32 dwarf and spiral galaxies with high-quality neutral gas, molecular gas, and stellar mass profiles, robust H\,{\sc i} rotation curves obtained via 3D kinematic modelling, and reliable bulge-disc decomposition. By assuming vertical hydrostatic equilibrium, we derive the scale heights of the atomic and molecular gas discs and fit dark matter haloes to the rotation curves self-consistently. We find that the effect of the gas flaring in the rotation curve decomposition can play an important role only for the smallest, gas-dominated dwarfs, while for most of the galaxies the effect is minor and can be ignored. We revisit the stellar- and baryon-to-halo mass relations ($M_\ast-M_{200}$ and $M_{\rm bar}-M_{200}$). Both relations increase smoothly up to $M_{200} \approx 10^{12}~\rm{ M_\odot}$, with galaxies at this end having high $M_\ast/M_{200}$ and $M_{\rm bar}/M_{200}$ ratios approaching the cosmological baryon fraction. At higher $M_{200}$ the relations show a larger scatter. Most haloes of our galaxy sample closely follow the concentration-mass ($c_{200}-M_{\rm 200}$) relation resulting from N-body cosmological simulations. Interestingly, the galaxies deviating above and below the relation have the highest and lowest stellar and baryon factions, respectively, which suggests that the departures from the $c_{200}-M_{\rm 200}$ law are regulated by adiabatic contraction and an increasing importance of feedback., Comment: accepted for publication in MNRAS. Supplementary material will be available from MNRAS and meantime at this link: https://www.dropbox.com/sh/768phag6ttgym0j/AABZwerBgFofkkZ5npHW17YBa

We present results from MUSE spatially-resolved spectroscopy of 21 post-starburst galaxies in the centers of 8 clusters from $z\sim0.3$ to $z\sim0.4$. We measure spatially resolved star-formation histories (SFHs), the time since quenching ($t_Q$) and the fraction of stellar mass assembled in the past 1.5 Gyr ($\mu_{1.5}$). The SFHs display a clear enhancement of star-formation prior to quenching for 16 out of 21 objects, with at least 10% (and up to $>50$%) of the stellar mass being assembled in the past 1.5 Gyr and $t_Q$ ranging from less than 100 Myrs to $\sim800$ Myrs. By mapping $t_Q$ and $\mu_{1.5}$, we analyze the quenching patterns of the galaxies. Most galaxies in our sample have quenched their star-formation from the outside-in or show a side-to-side/irregular pattern, both consistent with quenching by ram-pressure stripping. Only three objects show an inside-out quenching pattern, all of which are at the high-mass end of our sample. At least two of them currently host an active galactic nucleus. In two post-starbursts, we identify tails of ionized gas indicating that these objects had their gas stripped by ram pressure very recently. Post-starburst features are also found in the stripped regions of galaxies undergoing ram-pressure stripping in the same clusters, confirming the link between these classes of objects. Our results point to ram-pressure stripping as the main driver of fast quenching in these environments, with active galactic nuclei playing a role at high stellar masses., Comment: Accepted for publication in ApJ

Using Large Binocular Telescope deep imaging data from the Smallest Scale of Hierarchy Survey (SSH) and archival Hubble Space Telescope data, we reveal the presence of two elongated stellar features contiguous to a bar-like stellar structure in the inner regions of the dwarf irregular galaxy NGC 3741. These structures are dominated by stars younger than a few hundred Myr and collectively are about twice as extended as the old stellar component. These properties are very unusual for dwarf galaxies in the nearby Universe and difficult to explain by hydro-dynamical simulations. From the analysis of archival 21-cm observations, we find that the young stellar "bar" coincides with an HI high-density region proposed by previous studies to be a purely gaseous bar; we furthermore confirm radial motions of a few km/s, compatible with an inflow/outflow, and derive a steeply-rising rotation curve and high HI surface density at the center, indicating a very concentrated mass distribution. We propose that the peculiar properties of the stellar and gaseous components of NGC 3741 may be explained by a recent merger or ongoing gas accretion from the intergalactic medium, which caused gas inflows towards the galaxy center and triggered star formation a few hundred Myr ago. This event may explain the young and extended stellar features, the bar-like structure, the very extended HI disc and the central HI spiral arms. The high central HI density and the steeply rising rotation curve suggest that NGC 3741 may be the progenitor or the descendant of a starburst dwarf., Comment: Accepted for publication in MNRAS. Corrected typo in the list of authors

We present the results of a VLA HI imaging survey aimed at understanding why some galaxies develop long extraplanar H$\alpha$ tails, becoming extreme jellyfish galaxies. The observations are centered on five extreme jellyfish galaxies, optically selected from the WINGS and OmegaWINGS surveys and confirmed to have long H$\alpha$ tails through MUSE observations. Each galaxy is located in a different cluster. In the observations there are in total 88 other spiral galaxies within the field of view (40'x40') and observed bandwidth (6500 km s$^{-1}$). We detect 13 of these 88 spirals, plus one uncatalogued spiral, with HI masses ranging from 1 to 7 $\times$ 10${^9}$ M$_{\odot}$. Many of these detections have extended HI disks, two show direct evidence for ram pressure stripping, while others are possibly affected by tidal forces and/or ram-pressure stripping. We stack the 75 non-detected spiral galaxies and find an average HI mass of 1.9 $\times$ 10$^{8}$ M$_{\odot}$, which given their average stellar mass, implies they are very HI deficient. Comparing the extreme jellyfish galaxies to the other disk galaxies, we find that they have a larger stellar mass than almost all disk galaxies and than all HI detected galaxies, they are at smaller projected distance from the cluster center and at higher relative velocity to the cluster mean than all HI detections and most non-detections. We conclude that the high stellar mass allows extreme jellyfish galaxies to fall deeply into the cluster before being stripped and the surrounding ICM pressure gives rise to their spectacular star-forming tails., Comment: 25 pages, 13 figures, accepted to ApJ

With first-principles kinetic simulations, we show that a large-scale Alfv\'en wave (AW) propagating in an inhomogeneous background decays into kinetic Alfv\'en waves (KAWs), triggering ion and electron energization. We demonstrate that the two species can access unequal amounts of the initial AW energy, experiencing differential heating. During the decay process, the electric field carried by KAWs produces non-Maxwellian features in the particle VDFs, in accordance with space observations. The process we present solely requires the interaction of a large-scale AW with a magnetic shear and may be relevant for several astrophysical and laboratory plasmas., Comment: Accepted in PRL