BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Luca Del Zanna

Luca Del Zanna appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.HEhep-phastro-ph.SRgr-qcnucl-thnucl-exphysics.comp-phphysics.flu-dynphysics.plasm-ph

Catalog footprint

What is connected

11works
9topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2024arXiv

A GPU-Accelerated Modern Fortran Version of the ECHO Code for Relativistic Magnetohydrodynamics

The numerical study of relativistic magnetohydrodynamics (MHD) plays a crucial role in high-energy astrophysics, but unfortunately is computationally demanding, given the complex physics involved (high Lorentz factor flows, extreme magnetization, curved spacetimes near compact objects) and the large variety of spatial scales needed to resolve turbulent motions. A great benefit comes from the porting of existing codes running on standard processors to GPU-based platforms. However, this usually requires a drastic rewriting of the original code, the use of specific languages like CUDA, and a complex analysis of data management and optimization of parallel processes. Here we describe the porting of the ECHO code for special and general relativistic MHD to accelerated devices, simply based on native Fortran language built-in constructs, especially 'do concurrent' loops, few OpenACC directives, and the straightforward data management provided by the Unified Memory option of NVIDIA compilers.Thanks to these very minor modifications to the original code, the new version of ECHO runs at least 16 times faster on GPU platforms compared to CPU-based ones. The chosen benchmark is the 3D propagation of a relativistic MHD Alfvén wave, for which strong and weak scaling tests performed on the LEONARDO pre-exascale supercomputer at CINECA are provided (using up to 256 nodes corresponding to 1024 GPUs, and over 14 billion cells). Finally, an example of high-resolution relativistic MHD Alfvénic turbulence simulation is shown, demonstrating the potential for astrophysical plasmas of the new GPU-based version of ECHO.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Magnetic fields in heavy ion collisions: flow and charge transport

At the earliest times after a heavy-ion collision, the magnetic field created by the spectator nucleons will generate an extremely strong, albeit rapidly decreasing in time, magnetic field. The impact of this magnetic field may have detectable consequences, and is believed to drive anomalous transport effects like the Chiral Magnetic Effect (CME). We detail an exploratory study on the effects of a dynamical magnetic field on the hydrodynamic medium created in the collisions of two ultrarelativistic heavy-ions, using the framework of numerical ideal MagnetoHydroDynamics (MHD) with the ECHO-QGP code. In this study, we consider a magnetic field captured in a conducting medium, where the conductivity can receive contributions from the electromagnetic conductivity $σ$ and the chiral magnetic conductivity $σ_χ$. We first study the elliptic flow of pions, which we show is relatively unchanged by the introduction of a magnetic field. However, by increasing the magnitude of the magnetic field, we find evidence for an enhancement of the elliptic flow in peripheral collisions. Next, we explore the impact of the chiral magnetic conductivity on electric charges produced at the edges of the fireball. This initial $σ_χ$ can be understood as a long-wavelength effective description of chiral fermion production. We then demonstrate that this chiral charge, when transported by the MHD medium, produces a charge dipole perpendicular to the reaction plane which extends a few units in rapidity. Assuming charge conservation at the freeze-out surface, we show that the produced charge imbalance can have measurable effects on some experimental observables, like $v_1$ or $\langle \sin ϕ\rangle$. This demonstrates the ability of a MHD fluid to transport the signature of the initial chiral magnetic fields to late times.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Constraints on particle acceleration sites in the Crab Nebula from relativistic MHD simulations

The Crab Nebula is one of the most efficient accelerators in the Galaxy and the only galactic source showing direct evidence of PeV particles. In spite of this, the physical process behind such effective acceleration is still a deep mystery. While particle acceleration, at least at the highest energies, is commonly thought to occur at the pulsar wind termination shock, the properties of the upstream flow are thought to be non-uniform along the shock surface, and important constraints on the mechanism at work come from exact knowledge of where along this surface particles are being accelerated. Here we use axisymmetric relativistic MHD simulations to obtain constraints on the acceleration site(s) of particles of different energies in the Crab Nebula. Various scenarios are considered for the injection of particles responsible for synchrotron radiation in the different frequency bands, radio, optical and X-rays. The resulting emission properties are compared with available data on the multi wavelength time variability of the inner nebula. Our main result is that the X-ray emitting particles are accelerated in the equatorial region of the pulsar wind. Possible implications on the nature of the acceleration mechanism are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

General relativistic radiation hydrodynamics of accretion flows. I: Bondi-Hoyle accretion

We present a new code for performing general-relativistic radiation-hydrodynamics simulations of accretion flows onto black holes. The radiation field is treated in the optically-thick approximation, with the opacity contributed by Thomson scattering and thermal bremsstrahlung. Our analysis is concentrated on a detailed numerical investigation of hot two-dimensional, Bondi-Hoyle accretion flows with various Mach numbers. We find significant differences with respect to purely hydrodynamical evolutions. In particular, once the system relaxes to a radiation-pressure dominated regime, the accretion rates become about two orders of magnitude smaller than in the purely hydrodynamical case, remaining however super-Eddington as are the luminosities. Furthermore, when increasing the Mach number of the inflowing gas, the accretion rates become smaller because of the smaller cross section of the black hole, but the luminosities increase as a result a stronger emission in the shocked regions. Overall, our approach provides the first self-consistent calculation of the Bondi-Hoyle luminosity, most of which is emitted within r~100 M from the black hole, with typical values L/L_Edd ~ 1-7, and corresponding energy efficiencies eta_BH ~ 0.09-0.5. The possibility of computing luminosities self-consistently has also allowed us to compare with the bremsstrahlung luminosity often used in modelling the electromagnetic counterparts to supermassive black-hole binaries, to find that in the optically-thick regime these more crude estimates are about 20 times larger than our radiation-hydrodynamics results.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Resistive magnetohydrodynamic simulations of the ideal tearing mode

We study the linear and nonlinear evolution of the tearing instability on thin current sheets by means of two-dimensional numerical simulations, within the framework of compressible, resistive magnetohydrodynamics. In particular we analyze the behavior of current sheets whose inverse aspect ratio scales with the Lundquist number $S$ as $S^{-1/3}$. This scaling has been recently recognized to yield the threshold separating fast, ideal reconnection, with an evolution and growth which are independent of $S$ provided this is high enough, as it should be natural having the ideal case as a limit for $S\to\infty$. Our simulations confirm that the tearing instability growth rate can be as fast as $γ\approx 0.6\,{τ_A}^{-1}$, where $τ_A$ is the ideal Alfvénic time set by the macroscopic scales, for our least diffusive case with $S=10^7$. The expected instability dispersion relation and eigenmodes are also retrieved in the linear regime, for the values of $S$ explored here. Moreover, in the nonlinear stage of the simulations we observe secondary events obeying the same critical scaling with $S$, here calculated on the \emph{local}, much smaller lengths, leading to increasingly faster reconnection. These findings strongly support the idea that in a fully dynamic regime, as soon as current sheets develop, thin and reach this critical threshold in their aspect ratio, the tearing mode is able to trigger plasmoid formation and reconnection on the local (ideal) Alfvénic timescales, as required to explain the explosive flaring activity often observed in solar and astrophysical plasmas.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A perturbative approach to the hydrodynamics of heavy ion collisions

Initial fluctuations in hydrodynamic fields such as energy density or flow velocity give access to understanding initial state and equilibration physics as well as thermodynamic and transport properties. We provide evidence that the fluid dynamic propagation of fluctuations of realistic size can be based on a background-fluctuation splitting and a systematic perturbative expansion in the fluctuating fields. Initial conditions are characterized by a Bessel-Fourier expansion for single events, event-by-event correlations and probability distributions. The evolution equations can be solved order-by-order in the expansion which allows to study the fluid dynamical propagation of single modes, the study of interaction effects between modes, the determination of the associated particle spectra and the generalization of the whole program to event-by-event correlations and distributions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Heavy Ion Collision evolution modeling with ECHO-QGP

We present a numerical code modeling the evolution of the medium formed in relativistic heavy ion collisions, ECHO-QGP. The code solves relativistic hydrodynamics in $(3+1)-$D, with dissipative terms included within the framework of Israel-Stewart theory; it can work both in Minkowskian and in Bjorken coordinates. Initial conditions are provided through an implementation of the Glauber model (both Optical and Monte Carlo), while freezeout and particle generation are based on the Cooper-Frye prescription. The code is validated against several test problems and shows remarkable stability and accuracy with the combination of a conservative (shock-capturing) approach and the high-order methods employed. In particular it beautifully agrees with the semi-analytic solution known as Gubser flow, both in the ideal and in the viscous Israel-Stewart case, up to very large times and without any ad hoc tuning of the algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Numerical simulations of stellar jets and comparison between synthetic and observed maps: clues to the launch mechanism

High angular resolution spectra obtained with the Hubble Space Telescope Imaging Spectrograph (HST/STIS) provide rich morphological and kinematical information about the stellar jet phenomenon, which allows us to test theoretical models efficiently. In this work, numerical simulations of stellar jets in the propagation region are executed with the PLUTO code, by adopting inflow conditions that arise from former numerical simulations of magnetized outflows, accelerated by the disk-wind mechanism in the launching region. By matching the two regions, information about the magneto-centrifugal accelerating mechanism underlying a given astrophysical object can be extrapolated by comparing synthetic and observed position-velocity diagrams (PVDs). We show that quite different jets, like those from the young T Tauri stars DG-Tau and RW-Aur, may originate from the same disk-wind model for different configurations of the magnetic field at the disk surface. This result supports the idea that all the observed jets may be generated by the same mechanism.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

How (non-) linear is the hydrodynamics of heavy ion collisions?

We provide evidence from full numerical solutions that the hydrodynamical evolution of initial density fluctuations in heavy ion collisions can be understood order-by-order in a perturbative series in deviations from a smooth and azimuthally symmetric background solution. To leading linear order, modes with different azimuthal wave numbers do not mix. Quadratic and higher order corrections are small and can be understood as overtones with corresponding wave numbers.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On the MHD modeling of the Crab Nebula radio emission

In recent years, it has become a well-established paradigm that many aspects of the physics of Pulsar Wind Nebulae (PWNe) can be fully accounted for within a relativistic MHD description. Numerical simulations have proven extremely successful in reproducing the X-ray morphology of the Crab Nebula, down to very fine detail. Radio emission, instead, is currently one of the most obscure aspects of the physics of these objects, and one that holds important information about pulsar properties and their role as antimatter factories. Here we address the question of radio emission morphology and integrated spectrum from the Crab Nebula, by using for the first time an axisymmetric dynamical model with parameters chosen to best reproduce its X-ray morphology. Based on our findings we discuss constraints on the origin of the radio emitting particles.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

EM counterparts of recoiling black holes: general relativistic simulations of non-Keplerian discs

We investigate the dynamics of a circumbinary disc that responds to the loss of mass and to the recoil velocity of the black hole produced by the merger of a binary system of supermassive black holes. We perform the first two-dimensional general relativistic hydrodynamics simulations of \textit{extended} non-Keplerian discs and employ a new technique to construct a "shock detector", thus determining the precise location of the shocks produced in the accreting disc by the recoiling black hole. In this way we can study how the properties of the system, such as the spin, mass and recoil velocity of the black hole, affect the mass accretion rate and are imprinted on the electromagnetic emission from these sources. We argue that the estimates of the bremsstrahlung luminosity computed without properly taking into account the radiation transfer yield cooling times that are unrealistically short. At the same time we show, through an approximation based on the relativistic isothermal evolution, that the luminosity produced can reach a peak value above $L \simeq 10^{43} \ {\rm erg/s} $ at about $\sim 30\,{\rm d}$ after the merger of a binary with total mass $M\simeq 10^6 M_\odot$ and persist for several days at values which are a factor of a few smaller. If confirmed by more sophisticated calculations such a signal could indeed lead to an electromagnetic counterpart of the merger of binary black-hole system.

0 citations0 reviewsNo code linkedOpen access