Source author record

Doris Folini

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

15works
12topics
4close collaborators

Actions

Connect this record

Log in to claim

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

15 published item(s)

preprint2022arXiv

The climate in climate economics

To analyze climate change mitigation strategies, economists rely on simplified climate models - climate emulators. We propose a generic and transparent calibration and evaluation strategy for these climate emulators that is based on Coupled Model Intercomparison Project, Phase 5 (CMIP5). We demonstrate that the appropriate choice of the free model parameters can be of key relevance for the predicted social cost of carbon. We propose to use four different test cases: two tests to separately calibrate and evaluate the carbon cycle and temperature response, a test to quantify the transient climate response, and a final test to evaluate the performance for scenarios close to those arising from economic models. We re-calibrate the climate part of the widely used DICE-2016: the multi-model mean as well as extreme, but still permissible climate sensitivities and carbon cycle responses. We demonstrate that the functional form of the climate emulator of the DICE-2016 model is fit for purpose, despite its simplicity, but its carbon cycle and temperature equations are miscalibrated. We examine the importance of the calibration for the social cost of carbon in the context of a partial equilibrium setting where interest rates are exogenous, as well as the simple general equilibrium setting from DICE-2016. We find that the model uncertainty from different consistent calibrations of the climate system can change the social cost of carbon by a factor of four if one assumes a quadratic damage function. When calibrated to the multi-model mean, our model predicts similar values for the social cost of carbon as the original DICE-2016, but with a strongly reduced sensitivity to the discount rate and about one degree less long-term warming. The social cost of carbon in DICE-2016 is oversensitive to the discount rate, leading to extreme comparative statics responses to changes in preferences.

preprint2020arXiv

Anomalous atmospheric circulation favored the spread of COVID-19 in Europe

The current pandemic caused by the coronavirus SARS-CoV-2 is having negative health, social and economic consequences worldwide. In Europe, the pandemic started to develop strongly at the end of February and beginning of March 2020. It has subsequently spread over the continent, with special virulence in northern Italy and inland Spain. In this study we show that an unusual persistent anticyclonic situation prevailing in southwestern Europe during February 2020 (i.e. anomalously strong positive phase of the North Atlantic and Arctic Oscillations) could have resulted in favorable conditions, in terms of air temperature and humidity, in Italy and Spain for a quicker spread of the virus compared with the rest of the European countries. It seems plausible that the strong atmospheric stability and associated dry conditions that dominated in these regions may have favored the virus's propagation, by short-range droplet transmission as well as likely by long-range aerosol (airborne) transmission.

preprint2020arXiv

Arctic Oscillation: possible trigger of COVID-19 outbreak

The current COVID-19 pandemic is having detrimental consequences worldwide. The pandemic started to develop strongly by the end of January and beginning of February 2020, first in China with subsequent rapid spread to other countries with new epicenters of the outbreaks concentrated mainly within the 30-50 degrees North latitudinal band (e.g., South Korea, Japan, Iran, Italy, Spain). Simultaneously, an unusual persistent anticyclonic situation prevailing at latitudes around 40 degrees North was observed on global scale, in line with an anomalously strong positive phase of the Arctic Oscillation. This atypical situation could have resulted in favorable meteorological conditions for a quicker spread of the virus over the latitude band detailed above. This possible connection needs further attention in order to understand the meteorological and climatological factors related to the COVID-19 outbreak, and for anticipating the spatio-temporal distribution of possible future pandemics.

preprint2020arXiv

Structure of a collisionless pair jet in a magnetized electron-proton plasma: Flow-aligned magnetic field

We present the results from a particle-in-cell (PIC) simulation that models the interaction between a spatially localized electron-positron cloud and an electron-ion plasma. The latter is permeated by a magnetic field that is initially spatially uniform and aligned with the mean velocity vector of the pair cloud. The pair cloud expels the magnetic field and piles it up into an electromagnetic piston. Its electromagnetic field is strong enough to separate the pair cloud from the ambient plasma in the direction that is perpendicular to the cloud propagation direction. The piston propagates away from the spine of the injected pair cloud and it accelerates the protons to a high nonrelativistic speed. The accelerated protons form an outer cocoon that will eventually become separated from the unperturbed ambient plasma by a fast magnetosonic shock. No electromagnetic piston forms at the front of the cloud and a shock is mediated here by the filamentation instability. The final plasma distribution resembles that of a hydrodynamic jet. Collisionless plasma jets may form in the coronal plasma of accreting black holes and the interaction between the strong magnetic field of the piston and the hot pair cloud may contribute to radio emissions by such objects.

preprint2016arXiv

Accuracy requirements to test the applicability of the random cascade model to supersonic turbulence

A model, which is widely used for inertial rang statistics of supersonic turbulence in the context of molecular clouds and star formation, expresses (measurable) relative scaling exponents Z_p of two-point velocity statistics as a function of two parameters, beta and Delta. The model relates them to the dimension D of the most dissipative structures, D=3-Delta/(1-beta). While this description has proved most successful for incompressible turbulence (beta=Delta=2/3, and D=1), its applicability in the highly compressible regime remains debated. For this regime, theoretical arguments suggest D=2 and Delta=2/3, or Delta=1. Best estimates based on 3D periodic box simulations of supersonic isothermal turbulence yield Delta=0.71 and D=1.9, with uncertainty ranges of Delta in [0.67, 0.78] and D in [2.04,1.60]. With these 5-10\% uncertainty ranges just marginally including the theoretical values of Delta=2/3 and D=2, doubts remain whether the model indeed applies and, if it applies, for what values of beta and Delta. We use a Monte Carlo approach to mimic actual simulation data and examine what factors are most relevant for the fit quality. We estimate that 0.1% (0.05%) accurate Z_p, with p=1...5, should allow for 2% (1%) accurate estimates of beta and Delta in the highly compressible regime, but not in the mildly compressible regime. We argue that simulation-based Z_p with such accuracy are within reach of today's computer resources. If this kind of data does not allow for the expected high quality fit of beta and Delta, then this may indicate the inapplicability of the model for the simulation data. In fact, other models than the one we examine here have been suggested.

preprint2015arXiv

A Jacobian-free Newton-Krylov method for time-implicit multidimensional hydrodynamics

This work is a continuation of our efforts to develop an efficient implicit solver for multidimensional hydrodynamics for the purpose of studying important physical processes in stellar interiors, such as turbulent convection and overshooting. We present an implicit solver that results from the combination of a Jacobian-Free Newton-Krylov method and a preconditioning technique tailored to the inviscid, compressible equations of stellar hydrodynamics. We assess the accuracy and performance of the solver for both 2D and 3D problems for Mach numbers down to $10^{-6}$. Although our applications concern flows in stellar interiors, the method can be applied to general advection and/or diffusion-dominated flows. The method presented in this paper opens up new avenues in 3D modeling of realistic stellar interiors allowing the study of important problems in stellar structure and evolution.

preprint2014arXiv

Relativistic magnetic reconnection in collisionless ion-electron plasmas explored with particle-in-cell simulations

Magnetic reconnection is a leading mechanism for magnetic energy conversion and high-energy non-thermal particle production in a variety of high-energy astrophysical objects, including ones with relativistic ion-electron plasmas (e.g., microquasars or AGNs) - a regime where first principle studies are scarce. We present 2D particle-in-cell (PIC) simulations of low $β$ ion-electron plasmas under relativistic conditions, i.e., with inflow magnetic energy exceeding the plasma rest-mass energy. We identify outstanding properties: (i) For relativistic inflow magnetizations (here $10 < σ_e < 360$), the reconnection outflows are dominated by thermal agitation instead of bulk kinetic energy. (ii) At large inflow electron magnetization ($σ_e > 80$), the reconnection electric field is sustained more by bulk inertia than by thermal inertia. It challenges the thermal-inertia-paradigm and its implications. (iii) The inflows feature sharp transitions at the entrance of the diffusion zones. These are not shocks but results from particle ballistic motions, all bouncing at the same location, provided that the thermal velocity in the inflow is far smaller than the inflow E cross B bulk velocity. (iv) Island centers are magnetically isolated from the rest of the flow, and can present a density depletion at their center. (v) The reconnection rates are slightly larger than in non-relativistic studies. They are best normalized by the inflow relativistic Alfvén speed projected in the outflow direction, which then leads to rates in a close range (0.14-0.25) thus allowing for an easy estimation of the reconnection electric field.

preprint2014arXiv

Simulation of microquasars -- the challenge of scales

We present first results of a long-term project which aims at multi-scale, multi-physics simulations of wind accretion in microquasars and high-mass X-ray binaries. The 3D hydrodynamical simulations cover all scales, from the circum-binary environment down to the immediate vicinity of the black hole. We first introduce the numerical method and parallelization strategy of the AMR A-MAZE code. We then discuss some preliminary results of how, and on what scales, an accretion disk is formed around the black hole. We finally present some characteristics of this disk, which is far from Keplerian. We emphasize that on all scales shocks play a decisive role for the accretion process and the process of structure formation -- for the formation of the large scale, nearly coherent structure of the disk, but also for the formation of turbulent fluctuations.

preprint2014arXiv

Supersonic turbulence in 3D isothermal flow collision

Colliding supersonic bulk flows shape observable properties and internal physics of various astrophysical objects, like O-star winds, molecular clouds, galactic sheets, binaries, or gamma-ray bursts. Using numerical simulations, we show that the bulk flows leave a clear imprint on the collision zone, its mean properties and the turbulence it naturally develops. Our model setup consists of 3D head-on colliding isothermal hydrodynamical flows with Mach numbers between 2 and 43. Simulation results are in line with expectations from self-similarity: root mean square Mach numbers (Mrms) scale linearly with upstream Mach numbers, mean densities remain limited to a few times the upstream density. The density PDF is not log-normal. The turbulence is inhomogeneous: weaker in the zone center than close to the confining shocks. It is anisotropic: while Mrms is generally supersonic, Mrms transverse to the upstream flow is always subsonic. We argue that uniform, isothermal, head-on colliding flows generally disfavor isotropic, supersonic turbulence. The anisotropy carries over to other quantities like the density variance - Mach number relation. Structure functions differ depending on whether they are computed along a line-of-sight perpendicular or parallel to the upstream flow. We suggest that such line-of-sight effects should be kept in mind when interpreting turbulence characteristics derived from observations.

preprint2014arXiv

The energetics of relativistic magnetic reconnection: ion-electron repartition and particle distribution hardness

Collisionless magnetic reconnection is a prime candidate to account for flare-like or steady emission, outflow launching, or plasma heating, in a variety of high-energy astrophysical objects, including ones with relativistic ion-electron plasmas. But the fate of the initial magnetic energy in a reconnection event remains poorly known: what is the amount given to kinetic energy, the ion/electron repartition, and the hardness of the particle distributions? We explore these questions with 2D particle-in-cell simulations of ion-electron plasmas. We find that 45 to 75% of the total initial magnetic energy ends up in kinetic energy, this fraction increasing with the inflow magnetization. Depending on the guide field strength, ions get from 30 to 60% of the total kinetic energy. Particles can be separated into two populations that only weakly mix: (i) particles initially in the current sheet, heated by its initial tearing and subsequent contraction of the islands; and (ii) particles from the background plasma that primarily gain energy via the reconnection electric field when passing near the X-point. Particles (ii) tend to form a power-law with an index $p=-d\log n(γ)/d\logγ$, that depends mostly on the inflow Alfvén speed $V_A$ and magnetization $σ_s$ of species $s$, with for electrons $p=5$ to $1.2$ for increasing $σ_e$. The highest particle Lorentz factor, for ions or electrons, increases roughly linearly with time for all the relativistic simulations. This is faster, and the spectra can be harder, than for collisionless shock acceleration. We discuss applications to microquasar and AGN coronae, to extragalactic jets, and to radio lobes. We point out situations where effects such as Compton drag or pair creation are important.

preprint2013arXiv

Apar-T: code, validation, and physical interpretation of particle-in-cell results

We present the parallel particle-in-cell (PIC) code Apar-T and, more importantly, address the fundamental question of the relations between the PIC model, the Vlasov-Maxwell theory, and real plasmas. First, we present four validation tests: spectra from simulations of thermal plasmas, linear growth rates of the relativistic tearing instability and of the filamentation instability, and non-linear filamentation merging phase. For the filamentation instability we show that the effective growth rates measured on the total energy can differ by more than 50% from the linear cold predictions and from the fastest modes of the simulation. Second, we detail a new method for initial loading of Maxwell-Jüttner particle distributions with relativistic bulk velocity and relativistic temperature, and explain why the traditional method with individual particle boosting fails. Third, we scrutinize the question of what description of physical plasmas is obtained by PIC models. These models rely on two building blocks: coarse-graining, i.e., grouping of the order of p~10^10 real particles into a single computer superparticle, and field storage on a grid with its subsequent finite superparticle size. We introduce the notion of coarse-graining dependent quantities, i.e., quantities depending on p. They derive from the PIC plasma parameter Lambda^{PIC}, which we show to scale as 1/p. We explore two implications. One is that PIC collision- and fluctuation-induced thermalization times are expected to scale with the number of superparticles per grid cell, and thus to be a factor p~10^10 smaller than in real plasmas. The other is that the level of electric field fluctuations scales as 1/Lambda^{PIC} ~ p. We provide a corresponding exact expression. Fourth, we compare the Vlasov-Maxwell theory, which describes a phase-space fluid with infinite Lambda, to the PIC model and its relatively small Lambda.

preprint2013arXiv

Circum-stellar medium around rotating massive stars at solar metallicity

Aims. Observations show nebulae around some massive stars but not around others. If observed, their chemical composition is far from homogeneous. Our goal is to put these observational features into the context of the evolution of massive stars and their circumstellar medium (CSM) and, more generally, to quantify the role of massive stars for the chemical and dynamical evolution of the ISM. Methods. Using the A-MAZE code, we perform 2d-axisymmetric hydrodynamical simulations of the evolution of the CSM, shaped by stellar winds, for a whole grid of massive stellar models from 15 to 120 Msun and following the stellar evolution from the zero-age main-sequence to the time of supernova explosion. In addition to the usual quantities, we also follow five chemical species: H, He, C, N, and O. Results. We show how various quantities evolve as a function of time: size of the bubble, position of the wind termination shock, chemical composition of the bubble, etc. The chemical composition of the bubble changes considerably compared to the initial composition, particularly during the red-supergiant (RSG) and Wolf-Rayet (WR) phases. In some extreme cases, the inner region of the bubble can be completely depleted in hydrogen and nitrogen, and is mainly composed of carbon, helium and oxygen. We argue why the bubble typically expands at a lower rate than predicted by self-similarity theory. In particular, the size of the bubble is very sensitive to the density of the ISM, decreasing by a factor of around 2.5 for each additional dex in ISM density. The bubble size also decreases with the metallicity of the central star, as low-metallicity stars have weaker winds. Our models qualitatively fit the observations of WR ejecta nebulae.

preprint2013arXiv

Differences between real and particle-in-cell plasmas: effects of coarse-graining

The PIC model relies on two building blocks. The first stems from the capability of computers to handle only up to $\sim10^{10}$ particles, while real plasmas contain from $10^4$ to $10^{20}$ particles per Debye sphere: a coarse-graining step must be used, whereby of the order of $p\sim10^{10}$ real particles are represented by a single computer superparticle. The second is field storage on a grid with its subsequent finite superparticle size. We introduce the notion of coarse-graining dependent quantities, i.e. physical quantities depending on the number $p$. They all derive from the plasma parameter $Λ$, which we show to be proportional to $1/p$. We explore three examples: the rapid collision- and fluctuation-induced thermalization of plasmas with different temperatures, that scale with the number of superparticles per grid cell and are a factor $p\sim10^{10}$ faster than in real plasmas; the high level of electrostatic fluctuations in a thermal plasma, with corrections due to the finite superparticle sizes; and the blurring of the linear spectrum of the filamentation instability, where the fastest growing modes do not dominate the total energy because of a high level of fluctuations. We stress that the enhanced collisions and correlations of PIC plasmas must be kept negligible toward kinetic physics.

preprint2011arXiv

Magnetic fields in massive stars, their winds, and their nebulae

Massive stars are crucial building blocks of galaxies and the universe, as production sites of heavy elements and as stirring agents and energy providers through stellar winds and supernovae. The field of magnetic massive stars has seen tremendous progress in recent years. Different perspectives -- ranging from direct field measurements over dynamo theory and stellar evolution to colliding winds and the stellar environment -- fruitfully combine into a most interesting and still evolving overall picture, which we attempt to review here. Zeeman signatures leave no doubt that at least some O- and early B-type stars have a surface magnetic field. Indirect evidence, especially non-thermal radio emission from colliding winds, suggests many more. The emerging picture for massive stars shows similarities with results from intermediate mass stars, for which much more data are available. Observations are often compatible with a dipole or low order multi-pole field of about 1 kG (O-stars) or 300 G to 30 kG (Ap / Bp stars). Weak and unordered fields have been detected in the O-star zeta-Ori A and in Vega, the first normal A-type star with a magnetic field. Theory offers essentially two explanations for the origin of the observed surface fields: fossil fields, particularly for strong and ordered fields, or different dynamo mechanisms, preferentially for less ordered fields. Numerical simulations yield the first concrete stable (fossil) field configuration, but give contradictory results as to whether dynamo action in the radiative envelope of massive main sequence stars is possible. Internal magnetic fields, which may not even show up at the stellar surface, affect stellar evolution as they lead to a more uniform rotation, with more slowly rotating cores and faster surface rotation. Surface metallicities may become enhanced, thus affecting the mass-loss rates.

preprint2011arXiv

Simulating the circum-stellar environment of supernova and GRB progenitors by combining stellar evolution models and hydrodynamical code

The medium around massive stars is strongly shaped by the stellar winds. Those winds depend on various stellar parameters (effective temperature, luminosity, chemical composition, rotation, ...), which are varying as a function of the time. Using the wind properties obtained by classical stellar evolution code allows for the multi- D hydrodynamical simulation of the circum-stellar medium accounting for the time variations of the wind. We present here the preliminary results of the simulation of the medium around a fast rotating star. Comparing the results of such simulations with the observed properties of the circum-stellar medium will allow to better understand the interactions of the stellar winds with the interstellar medium, and could give hints on the past history of the mass loss around massive stars, providing constraints for stellar models.