Source author record

Alexandre Marcowith

Alexandre Marcowith 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

16works
5topics
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

16 published item(s)

preprint2026arXiv

Energetic particles accelerated via turbulent magnetic reconnection in protoplanetary discs -- I. Ionisation rates

Context. Ionisation controls the chemistry, thermal balance, and magnetic coupling in protoplanetary discs. However, standard ionisation vectors such as stellar UV, X-rays, Galactic cosmic rays (GCRs) might not be efficient enough, as UV/X-rays are attenuated rapidly with depth, while GCRs are modulated. Turbulence-induced magnetic reconnection in disc atmospheric layers offers a physically motivated, in situ source of energetic particles (EPs) that has never been considered. Aims. We quantify the ionisation and heating produced by EPs accelerated by turbulent reconnection, identify where they dominate over X-rays and GCRs, and determine energetic thresholds for their relevance. We provide scalable diagnostics tied to the local energy budget. Methods. We adopt a Fermi-like acceleration model with parameters linked to a turbulent reconnection geometry trigger by the magneto-rotational instability, yielding a steady-state energy distribution of the EP forming a power-law of index $p=2.5$. We propagate electrons and protons through the disc and compute primary and secondary ionisation and associated heating on a fiducial T Tauri disc model background. The non-thermal normalisation is set by the fraction of local viscous accretion energy dissipation channelled to EPs, parametrised by $κ$. Results. For $κ\gtrsim 0.4\%$, EPs ionisation overpass standard ionisation sources in the disc atmosphere and intermediate/deep layers out to radii of a few tens of astronomical units. Even at $κ\sim 0.025\%$, EPs contribute at the few-percent level, thus are chemically and dynamically relevant. These results identify EPs accelerated by turbulence-induced magnetic reconnection as a rather robust, disc-internal ionisation channel that should be included in thermo-chemical models of protoplanetary discs.

preprint2025arXiv

Pulsar+escape: the recipe for LHAASO sources ?

Observations of our Galaxy at very-high and ultra-high photon energies have revealed a rich population of sources, many of which have significant angular extension and/or are positionally coincident with energetic pulsars. We assessed whether the properties of these gamma-ray emitters are consistent with the predictions of a model for the non-thermal leptonic emission of a pulsar wind nebula expanding within a supernova remnant, including the possibility of particle escape across the different components of the system up to the ambient medium. We used a multi-zone model framework that describes the transport and radiation of a population of ultra-relativistic electron-positron pairs. The dynamics of the system is retrieved from numerical hydrodynamics simulations, which allows to follow the evolution through advanced stages. Model predictions for typical parameter setups are compared to source catalogs issued by the H.E.S.S. and LHAASO collaborations. If particles remain confined in the nebula, 1-100TeV emission is only partially consistent with the properties of observed sources in terms of surface brightness, angular extent, and photon index. In particular, the systems spend a long time in a bright and soft state that has no equivalent in observations. Conversely, including the possibility of energy-dependent particle escape across the object results in a much better agreement, both in the 1-10TeV range explored with H.E.S.S. or LHAASO-WCDA, and in the 20-100TeV range explored by LHAASO-KM2A. Emission components with intermediate to low brightnesses and large to very large sizes result from particles trapped in the remnant or escaped in the surrounding medium. Particle escape, clearly seen in the form of misaligned jets or halos around middle-aged pulsars, is a very important process in much younger systems with 1-100kyr ages, and shapes their appearance as TeV/PeV sources.

preprint2022arXiv

Diffusive shock acceleration at oblique high Mach number shocks

The current paradigm of cosmic ray (CR) origin states that the most part of galactic CRs is produced by supernova remnants. The interaction of supernova ejecta with the interstellar medium after supernova's explosions results in shocks responsible for CR acceleration via diffusive shock acceleration (DSA). We use particle-in-cell (PIC) simulations and a combined PIC-magnetohydrodynamic (PIC-MHD) technique to investigate whether DSA can occur in oblique high Mach number shocks. Using the PIC method, we follow the formation of the shock and determine the fraction of the particles that gets involved in DSA. Then, with this result, we use PIC-MHD simulations to model the large-scale structure of the plasma and the magnetic field surrounding the shock and find out whether or not the reflected particles can generate the upstream turbulence and trigger DSA. We find that the feasibility of this process in oblique shocks depends strongly on the Alfvenic Mach number, and the DSA process is more likely triggered at high Mach number shocks.

preprint2021arXiv

A diffusive origin for the cosmic-ray spectral hardening reveals signatures of a nearbysource in the leptons and protons data

In this work we aim at reproducing, simultaneously, the spectral feature at $\sim 10 \, \mathrm{TeV}$ in the cosmic-ray proton spectrum, recently reported by the DAMPE Collaboration, together with the spectral break at $\sim 1 \, \mathrm{TeV}$ measured by H.E.S.S. in the lepton spectrum. Those features are interpreted as signatures of one nearby hidden cosmic-ray accelerator. We show that this interpretation is consistent with the dipole-anisotropy data as long as the rigidity scaling of the diffusion coefficient features a hardening at $\sim 200 \, \mathrm{GV}$, as suggested by the light-nuclei data measured with high accuracy by the AMS-02 Collaboration. Such rigidity-dependent diffusion coefficient is applied consistently to the large-scale diffuse cosmic-ray sea as well as to the particles injected by the nearby source.

preprint2020arXiv

Impact of low-energy cosmic rays on star formation

In recent years, exciting developments have taken place in the identification of the role of cosmic rays in star-forming environments. Observations from radio to infrared wavelengths and theoretical modelling have shown that low-energy cosmic rays (<1 TeV) play a fundamental role in shaping the chemical richness of the interstellar medium, determining the dynamical evolution of molecular clouds. In this review we summarise in a coherent picture the main results obtained by observations and by theoretical models of propagation and generation of cosmic rays, from the smallest scales of protostars and circumstellar discs, to young stellar clusters, up to Galactic and extragalactic scales. We also discuss the new fields that will be explored in the near future thanks to new generation instruments, such as: CTA, for the $γ$-ray emission from high-mass protostars; SKA and precursors, for the synchrotron emission at different scales; and ELT/HIRES, JWST, and ARIEL, for the impact of cosmic rays on exoplanetary atmospheres and habitability.

preprint2020arXiv

Time-dependent high-energy gamma-ray signal from accelerated particles in core-collapse supernovae: the case of SN 1993J

Some core-collapse supernovae are likely to be efficient cosmic-ray accelerators up to the PeV range, and therefore, to potentially play an important role in the overall Galactic cosmic-ray population. The TeV gamma-ray domain can be used to study particle acceleration in the multi-TeV and PeV range. This motivates the study of the detectability of such supernovae by current and future gamma-ray facilities. The gamma-ray emission of core-collapse supernovae strongly depends on the level of the two-photon annihilation process: high-energy gamma-ray photons emitted at the expanding shock wave following the supernova explosion can interact with soft photons from the supernova photosphere through the pair production channel, thereby strongly suppressing the flux of gamma rays leaving the system. In the case of SN 1993J, whose photospheric and shock-related parameters are well measured, we calculate the temporal evolution of the expected gamma-ray attenuation by accounting for the temporal and geometrical effects. We find the attenuation to be of about $10$ orders of magnitude in the first few days after the SN explosion. The probability of detection of a supernova similar to SN 1993J with the Cherenkov Telescope Array is highest if observations are performed either earlier than 1 day, or later than 10 days after the explosion, when the gamma-ray attenuation decreases to about $2$ orders of magnitude.

preprint2016arXiv

Cosmic Ray propagation in sub-Alfvenic magnetohydrodynamic turbulence

This work has the main objective to provide a detailed investigation of cosmic ray propagation in magnetohydrodynamic turbulent fields generated by forcing the fluid velocity field at large scales. It provides a derivation of the particle mean free path dependences in terms of the turbulence level described by the Alfvénic Mach number and in terms of the particle rigidity. We use an upgrade version of the magnetohydrodynamic code {\tt RAMSES} which includes a forcing module and a kinetic module and solve the Lorentz equation for each particle. The simulations are performed using a 3 dimension periodical box in the test-particle and magnetostatic limits. The forcing module is implemented using an Ornstein-Uhlenbeck process. An ensemble average over a large number of particle trajectories is applied to reconstruct the particle mean free paths. We derive the cosmic ray mean free paths in terms of the Alfvénic Mach numbers and particle reduced rigidities in different turbulence forcing geometries. The reduced particle rigidity is $ρ=r_L/L$ where $r_L$ is the particle Larmor radius and $L$ is the simulation box length related to the turbulence coherence or injection scale $L_{inj}$ by $L \sim 5 L_{inj}$. We have investigated with a special attention compressible and solenoidal forcing geometries. We find that compressible forcing solutions are compatible with the quasi-linear theory or more advanced non-linear theories which predict a rigidity dependence as $ρ^{1/2}$ or $ρ^{1/3}$. Solenoidal forcing solutions at least at low or moderate Alfvénic numbers are not compatible with the above theoretical expectations and require more refined arguments to be interpreted. It appears especially for Alfvénic Mach numbers close to one that the wandering of field lines controls perpendicular mean free path solutions whatever the forcing geometry.

preprint2016arXiv

Embedded star clusters as sources of high-energy cosmic rays: Modelling and constraints

Massive stars are mainly found in stellar associations. These massive star clusters occur in the heart of giant molecular clouds. The strong stellar wind activity in these objects generates large bubbles and induces collective effects that could accelerate particles up to high energy and produce gamma rays. The best way to input an acceleration origin to the stellar wind interaction in massive stellar cluster is to observe young massive star clusters in which no supernova explosion has occurred yet. This work aims to constrain the part of stellar wind mechanical energy that is converted into energetic particles using the sensitivity of the ongoing Fermi/LAT instrument. This work further provides detailed predictions of expected gamma-ray fluxes in the view of the on-set of the next generation of imaging atmospheric Cherenkov telescopes. A one-zone model where energetic particles are accelerated by repeated interactions with strong supersonic shocks occurring in massive star clusters was developed. The particle escape from the star cluster and subsequent interaction with the surrounding dense material and magnetic fields of the HII region was computed. We applied this model to a selection of eight embedded star clusters constricted by existing observations. We evaluated the gamma-ray signal from each object, combining both leptonic and hadronic contributions. We searched for these emissions in the Fermi/LAT observations in the energy range from 3 to 300 GeV and compared them to the sensitivity of the Cherenkov Telescope Array. No significant gamma-ray emission from these star clusters has been found. Less than 10% of stellar wind luminosities are supplied to the relativistic particles. Some clusters even show acceleration efficiency of less than 1%. The CTA would be able to detect gamma-ray emission from several clusters in the case of an acceleration efficiency of close to 1%.

preprint2016arXiv

Non-linear diffusion of cosmic rays escaping from supernova remnants - I. The effect of neutrals

Supernova remnants are believed to be the main sources of galactic Cosmic Rays (CR). Within this framework, particles are accelerated at supernova remnant shocks and then released in the interstellar medium. The mechanism through which CRs are released and the way in which they propagate still remain open issues. The main difficulty is the high non-linearity of the problem: CRs themselves excite the magnetic turbulence that confines them close to their sources. We solve numerically the coupled differential equations describing the evolution in space and time of the escaping particles and of the waves generated through the CR streaming instability. The warm ionized and warm neutral phases of the interstellar medium are considered. These phases occupy the largest fraction of the disc volume, where most supernovae explode, and are characterised by the significant presence of neutral particles. The friction between those neutrals and ions results in a very effective wave damping mechanism. It is found that streaming instability affects the propagation of CRs even in the presence of ion-neutral friction. The diffusion coefficient can be suppressed by more than a factor of $\sim 2$ over a region of few tens of pc around the remnant. The suppression increases for smaller distances. The propagation of $\approx 10$ GeV particles is affected for several tens of kiloyears after escape, while $\approx 1$ TeV particles are affected for few kiloyears. This might have a great impact on the interpretation of gamma-ray observations of molecular clouds located in the vicinity of supernova remnants.

preprint2015arXiv

Cosmic-ray acceleration in young protostars

The main signature of the interaction between cosmic rays and molecular clouds is the high ionisation degree. This decreases towards the densest parts of a cloud, where star formation is expected, because of energy losses and magnetic effects. However recent observations hint to high levels of ionisation in protostellar systems, therefore leading to an apparent contradiction that could be explained by the presence of energetic particles accelerated within young protostars. Our modelling consists of a set of conditions that has to be satisfied in order to have an efficient particle acceleration through the diffusive shock acceleration mechanism. We find that jet shocks can be strong accelerators of protons which can be boosted up to relativistic energies. Another possibly efficient acceleration site is located at protostellar surfaces, where shocks caused by impacting material during the collapse phase are strong enough to accelerate protons. Our results demonstrate the possibility of accelerating particles during the early phase of a proto-Solar-like system and can be used as an argument to support available observations. The existence of an internal source of energetic particles can have a strong and unforeseen impact on the star and planet formation process as well as on the formation of pre-biotic molecules.

preprint2015arXiv

Non-linear Cosmic Ray propagation close to the acceleration site

Recent advances on gamma-ray observations from SuperNova Remnants and Molecular Clouds offer the possibility to study in detail the properties of the propagation of escaping Cosmic Rays (CR). However, a complete theory for CR transport outside the acceleration site has not been developed yet. Two physical processes are thought to be relevant to regulate the transport: the growth of waves caused by streaming instability, and possible wave damping mechanisms that reduce the growth of the turbulence. Only a few attempts have been made so far to incorporate these mechanisms in the theory of CR diffusion. In this work we present recent advances in this subject. In particular, we show results obtained by solving the coupled equations for the diffusion of CRs and the evolution of Alfven waves. We discuss the importance of streaming instabilities and wave damping in different ISM phases.

preprint2015arXiv

Shape and evolution of wind-blown bubbles of massive stars: on the effect of the interstellar magnetic field

The winds of massive stars create large (>10 pc) bubbles around their progenitors. As these bubbles expand they encounter the interstellar coherent magnetic field which, depending on its strength, can influence the shape of the bubble. We wish to investigate if, and how much, the interstellar magnetic field can contribute to the shape of an expanding circumstellar bubble around a massive star. We use the MPI-AMRVAC code to make magneto-hydrodynamical simulations of bubbles, using a single star model, combined with several different field strengths: B = 5, 10, and 20 muG for the interstellar magnetic field. This covers the typical field strengths of the interstellar magnetic fields found in the galactic disk and bulge. Furthermore, we present two simulations that include both a 5 muG interstellar magnetic field and a 10,000 K interstellar medium and two different ISM densities to demonstrate how the magnetic field can combine with other external factors to influence the morphology of the circumstellar bubbles. Our results show that low magnetic fields, as found in the galactic disk, inhibit the growth of the circumstellar bubbles in the direction perpendicular to the field. As a result, the bubbles become ovoid, rather than spherical. Strong interstellar fields, such as observed for the galactic bulge, can completely stop the expansion of the bubble in the direction perpendicular to the field, leading to the formation of a tube-like bubble. When combined with a warm, high-density ISM the bubble is greatly reduced in size, causing a dramatic change in the evolution of temporary features inside the bubble. The magnetic field of the interstellar medium can affect the shape of circumstellar bubbles. This effect may have consequences for the shape and evolution of circumstellar nebulae and supernova remnants, which are formed within the main wind-blown bubble.

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

Stochastic acceleration in accreting black holes

We present an extension of the BELM code (Belmont et al 2008) to investigate the microphysics of particle acceleration in black holes accretion disc corona. The updated version of the code accounts for the dynamics of resonant slab waves as well as their interaction with both leptons or protons. It is found that the proton temperature is an important regulating effect of the stochastic particle acceleration process in accretion disk corona. We present a preliminary fit of the high soft spectral state of Cygnus X-1.

preprint2012arXiv

A hydrodynamical model of the circumstellar bubble created by two massive stars

Numerical models of the wind-blown bubble of massive stars usually only account for the wind of a single star. However, since massive stars are usually formed in clusters, it would be more realistic to follow the evolution of a bubble created by several stars. We develope a two-dimensional (2D) model of the circumstellar bubble created by two massive stars, a 40 solar mass star and a 25 solar mass star, and follow its evolution. The stars are separated by approximately 16 pc and surrounded by a cold medium with a density of 20 particles per cubic cm. We use the MPI-AMRVAC hydrodynamics code to solve the conservation equations of hydrodynamics on a 2D cylindrical grid using time-dependent models for the wind parameters of the two stars. At the end of the stellar evolution (4.5 and 7.0 million years for the 40 and 25 solar mass stars, respectively), we simulate the supernova explosion of each star. Each star initially creates its own bubble. However, as the bubbles expand they merge, creating a combined, aspherical bubble. The combined bubble evolves over time, influenced by the stellar winds and supernova explosions. The evolution of a wind-blown bubble created by two stars deviates from that of the bubbles around single stars. In particular, once one of the stars has exploded, the bubble is too large for the wind of the remaining star to maintain and the outer shell starts to disintegrate. The lack of thermal pressure inside the bubble also changes the behavior of circumstellar features close to the remaining star. The supernovae are contained inside the bubble, which reflects part of the energy back into the circumstellar medium.

preprint2009arXiv

On the shape of the spectrum of cosmic-rays accelerated inside superbubbles

Supernova remnants are believed to be a major source of cosmic-rays in the Galaxy. As their progenitors are commonly found clustered in OB associations, one has to consider the possibility of collective effects in the acceleration process. In this work we investigate the shape of the spectrum of protons accelerated inside the superbubbles blown around clusters of massive stars. To do so we embed semi-analytical models of particle acceleration and transport inside Monte-Carlo simulations of OB associations timelines. We consider regular acceleration at the shock front of supernova remnants, as well as stochastic re-acceleration and escape occurring between the shocks. We observe that particle spectra, although strongly intermittent, get a distinctive shape resulting from a competition between acceleration and escape: they are harder at the lowest energies and softer at the highest energies. The momentum at which this spectral break occurs depends on a single dimensionless parameter, which we evaluate for a selection of objects. The behaviour of a superbubble regarding acceleration depends on the magnetic turbulence: if B is low the superbubble is simply the host of a collection of individual supernovae shocks, but if B is high enough it acts as a global accelerator, producing distinctive spectra ? which has important implications on the high-energy emission from these objects.