BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Prateek Sharma

Prateek Sharma contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.GAastro-ph.HEastro-ph.COphysics.flu-dynDistributed, Parallel, and Cluster Computingeess.SYgr-qcmath.NANumerical AnalysisPerformancephysics.ao-phSystems and Control

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
17works
0followers
12topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

17 published item(s)

preprint2023arXiv

Dissipation of AGN jets in a clumpy interstellar medium

Accreting supermassive black holes (SMBHs) frequently power jets that interact with the interstellar/circumgalactic medium (ISM/CGM), regulating star-formation in the galaxy. Highly supersonic jets launched by active galactic nuclei (AGN) power a cocoon that confines them and shocks the ambient medium. We build upon the models of narrow conical jets interacting with a smooth ambient medium, to include the effect of dense clouds that are an essential ingredient of a multiphase ISM. The key physical ingredient of this model is that the clouds along the supersonic jet-beam strongly decelerate the jet-head, but the subsonic cocoon easily moves around the clouds without much resistance. We propose scalings for important physical quantities -- cocoon pressure, head & cocoon speed, and jet radius. We obtain, for the first time, the analytic condition on clumpiness of the ambient medium for the jet to dissipate within the cocoon and verify it with numerical simulations of conical jets interacting with a uniform ISM with embedded spherical clouds. A jet is defined to be dissipated when the cocoon speed exceeds the speed of the jet-head. We compare our models to more sophisticated numerical simulations, direct observations of jet-ISM interaction (e.g., quasar J1316+1753), and discuss implications for the Fermi/eROSITA bubbles. Our work also motivates effective subgrid models for AGN jet feedback in a clumpy ISM unresolved by the present generation of cosmological galaxy formation simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Cosmic rays from massive star clusters : A close look at Westerlund 1

We study the effect of cosmic ray (CR) acceleration in the massive compact star cluster Westerlund 1 in light of its recent detection in $γ$-rays. Recent observations reveal a $1/r$ radial distribution of the CR energy density. Here we theoretically investigate whether or not this profile can help to distinguish between (1) continuous CR acceleration in the star cluster stellar wind-driven shocks and (2) discrete CR acceleration in multiple supernovae shocks -- which are often debated in the literature. Using idealized two-fluid simulations and exploring different acceleration sites and diffusion coefficients, we obtain the CR energy density profile and luminosity to find the best match for the $γ$-ray observations. We find that the inferred CR energy density profiles from observations of $γ$-ray luminosity and mass can be much different from the true radial profile. CR acceleration at either the cluster core region or the wind termination shock can explain the observations, if the diffusion coefficient is $κ_{\rm cr}\sim 10^{27}$ cm$^2$ s$^{-1}$ and a fraction of $\approx 10\%-20\%$ of the shock power/post-shock pressure is deposited into the CR component. We also study the possibility of discrete supernovae (SN) explosions being responsible for CR acceleration and find that with an injection rate of 1 SN in every $\sim 0.03$ Myr, one can explain the observed $γ$-ray profile. This multiple SN scenario is consistent with X-ray observations only if the thermal conductivity is close to the Spitzer value.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Cosmological evolution of gas and supermassive black holes in idealized isolated halos

We study the evolution of baryonic gas in cosmologically growing dark matter halos. To accurately model both the inner and outer regions of the halos, we use a dark matter density profile that transitions smoothly from the NFW profile within the virial radius to a more realistic flat profile far beyond the halo. We construct a dark matter gravitational potential consistent with this density profile, and we use a "cosmological" potential that accounts for gas evolution consistent with Hubble expansion at large radii. Gas is initialized with a density $\approx$ 0.2 times the dark matter density, consistent with the universal baryon fraction $ρ_{\rm g}/(ρ_{\rm g}+ρ_{\rm DM}) \approx 0.17$. We study the formation of the virial shock and evolution of the baryon fraction, including the effects of radiative cooling and AGN jet feedback. The feedback is powered by the accretion of cold gas onto a central supermassive black hole (SMBH). The cores of the halo exhibit heating and cooling cycles, whose strength and duration depend on the feedback efficiency and the halo mass. The central SMBH initially grows exponentially with time in the early quasar phase, but the growth slows down at later times. The baryon fraction in the core decreases with increasing feedback efficiency and decreasing halo mass. While the halo outskirts evolve self-similarly, the core density is non-evolving, in agreement with cluster observations. We analyze the correlations between the properties of the gas and the central SMBH, and explore the existence of a fundamental plane.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Deprojection of X-ray data in galaxy clusters: confronting simulations with observations

Numerical simulations with varying realism indicate an emergent principle -- multiphase condensation and large cavity power occur when the ratio of the cooling time to the free-fall time ($t_{\rm cool}/t_{\rm ff}$) falls below a threshold value close to 10. Observations indeed show cool-core signatures when this ratio falls below 20-30, but the prevalence of cores with $t_{\rm cool}/t_{\rm ff}$ ratio below 10 is rare as compared to simulations. In X-ray observations, we obtain projected spectra from which we have to infer radial gas density and temperature profiles. Using idealized models of X-ray cavities and multiphase gas in the core and 3-D hydro jet-ICM simulations, we quantify the biases introduced by deprojection based on the assumption of spherical symmetry in determining $t_{\rm cool}/t_{\rm ff}$. We show that while the used methods are able to recover the $t_{\rm cool}/t_{\rm ff}$ ratio for relaxed clusters, they have an uncertainty of a factor of $2-3$ in systems containing large cavities ($\gtrsim 20$ kpc). We also show that the mass estimates from these methods, in the absence of X-ray spectra close to the virial radius, suffer from a degeneracy between the virial mass ($M_{200}$) and the concentration parameter ($c$) in the form of $M_{200} c^2 \approx$ constant. Additionally, lack of soft-X-ray ($\lesssim 0.5$ keV) coverage and poor spatial resolution make us overestimate min($t_{\rm cool}/t_{\rm ff}$) by a factor of few in clusters with min($t_{\rm cool}/t_{\rm ff}$) $\lesssim 5$. This bias can largely explain the lack of cool-core clusters with min($t_{\rm cool}/t_{\rm ff}$) $\lesssim 5$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Multiphase turbulence in galactic halos: effect of the driving

Supernova explosions, active galactic nuclei jets, galaxy--galaxy interactions and cluster mergers can drive turbulence in the circumgalactic medium (CGM) and in the intracluster medium (ICM). However, the exact nature of turbulence forced by these sources and its impact on the different statistical properties of the CGM/ICM and their global thermodynamics is still unclear. To investigate the effects of different types of forcing, we conduct high resolution ($1008^3$ resolution elements) idealised hydrodynamic simulations with purely solenoidal (divergence-free) forcing, purely compressive (curl-free) forcing, and natural mixture forcing (equal fractions of the two components). The simulations also include radiative cooling. We study the impact of the three different forcing modes (sol, comp, mix) on the morphology of the gas, its temperature and density distributions, sources and sinks of enstrophy, i.e., solenoidal motions, as well as the kinematics of hot ($\sim10^7~\mathrm{K}$) X-ray emitting and cold ($\sim10^4~\mathrm{K}$) H$α$ emitting gas. We find that compressive forcing leads to stronger variations in density and temperature of the gas as compared to solenoidal forcing. The cold phase gas forms large-scale filamentary structures for compressive forcing and misty, small-scale clouds for solenoidal forcing. The cold phase gas has stronger large-scale velocities for compressive forcing. The natural mixture forcing shows kinematics and gas distributions intermediate between the two extremes, the cold-phase gas occurs as both large-scale filaments and small-scale misty clouds.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Shadow and weak gravitational lensing of a rotating regular black hole in a non-minimally coupled Einstein-Yang-Mills theory in the presence of plasma

The null geodesics of the regular and rotating magnetically charged black hole in a non-minimally coupled Einstein-Yang-Mills theory surrounded by a plasma medium is studied. The effect of magnetic charge and Yang-Mills parameter on the effective potential and radius of photon orbits has investigated. We then study the shadow of a regular and rotating magnetically charged black hole along with the observables in presence of the plasma medium. The presence of plasma medium affects the apparent size of the shadow of a regular rotating black hole in comparison to vacuum case. Variation of shadow radius and deformation parameter with Yang-Mills and plasma parameter has examined. Furthermore, the deflection angle of the massless test particles in weak field approximation around this black hole spacetime in presence of homogeneous plasma medium is also investigated. Finally, we have compared the obtained results with Kerr-Newman and Schwarzschild black hole solutions in general relativity (GR).

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A method for determining the parameters in a rheological model for viscoelastic materials by minimizing Tikhonov functionals

Mathematical models describing the behavior of viscoelastic materials are often based on evolution equations that measure the change in stress depending on its material parameters such as stiffness, viscosity or relaxation time. In this article, we introduce a Maxwell-based rheological model, define the associated forward operator and the inverse problem in order to determine the number of Maxwell elements and the material parameters of the underlying viscoelastic material. We perform a relaxation experiment by applying a strain to the material and measure the generated stress. Since the measured data varies with the number of Maxwell elements, the forward operator of the underlying inverse problem depends on parts of the solution. By introducing assumptions on the relaxation times, we propose a clustering algorithm to resolve this problem. We provide the calculations that are necessary for the minimization process and conclude with numerical results by investigating unperturbed as well as noisy data. We present different reconstruction approaches based on minimizing a least squares functional. Furthermore, we look at individual stress components to analyze different displacement rates. Finally, we study reconstructions with shortened data sets to obtain assertions on how long experiments have to be performed to identify conclusive material parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A numerical approach to the non-uniqueness problem of cosmic ray two-fluid equations at shocks

Cosmic rays (CRs) are frequently modeled as an additional fluid in hydrodynamic (HD) and magnetohydrodynamic (MHD) simulations of astrophysical flows. The standard CR two-fluid model is described in terms of three conservation laws (expressing conservation of mass, momentum and total energy) and one additional equation (for the CR pressure) that cannot be cast in a satisfactory conservative form. The presence of non-conservative terms with spatial derivatives in the model equations prevents a unique weak solution behind a shock. We investigate a number of methods for the numerical solution of the two-fluid equations and find that, in the presence of shock waves, the results generally depend on the numerical details (spatial reconstruction, time stepping, the CFL number, and the adopted discretization). All methods converge to a unique result if the energy partition between the thermal and non-thermal fluids at the shock is prescribed using a subgrid prescription. This highlights the non-uniqueness problem of the two-fluid equations at shocks. From our numerical investigations, we report a robust method for which the solutions are insensitive to the numerical details even in absence of a subgrid prescription, although we recommend a subgrid closure at shocks using results from kinetic theory. The subgrid closure is crucial for a reliable post-shock solution and also its impact on large scale flows because the shock microphysics that determines CR acceleration is not accurately captured in a fluid approximation. Critical test problems, limitations of fluid modeling, and future directions are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Characterising the Dynamo in a Radiatively Inefficient Accretion Flow

We explore the MRI driven dynamo in a radiatively inefficient accretion flow (RIAF) using the mean field dynamo paradigm. Using singular value decomposition (SVD) we obtain the least squares fitting dynamo coefficients $α$ and $γ$ by comparing the time series of the turbulent electromotive force and the mean magnetic field. Our study is the first one to show the poloidal distribution of these dynamo coefficients in global accretion flow simulations. Surprisingly, we obtain a high value of the turbulent pumping coefficient $γ$ which transports the mean magnetic flux radially outward. This would have implications for the launching of magnetised jets which are produced efficiently in presence a large-scale poloidal magnetic field close to the compact object. We present a scenario of a truncated disc beyond the RIAF where a large scale dynamo-generated poloidal magnetic field can aid jet-launching close to the black hole. Magnitude of all the calculated coefficients decreases with radius. Meridional variations of $α_{ϕϕ}$, responsible for toroidal to poloidal field conversion, is very similar to that found in shearing box simulations using the `test field' (TF) method. By estimating the relative importance of $α$-effect and shear, we conclude that the MRI driven large-scale dynamo, which operates at high latitudes beyond a disc scale height, is essentially of the $α-Ω$ type.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Cloud-scale VM Deflation for Running Interactive Applications On Transient Servers

Transient computing has become popular in public cloud environments for running delay-insensitive batch and data processing applications at low cost. Since transient cloud servers can be revoked at any time by the cloud provider, they are considered unsuitable for running interactive application such as web services. In this paper, we present VM deflation as an alternative mechanism to server preemption for reclaiming resources from transient cloud servers under resource pressure. Using real traces from top-tier cloud providers, we show the feasibility of using VM deflation as a resource reclamation mechanism for interactive applications in public clouds. We show how current hypervisor mechanisms can be used to implement VM deflation and present cluster deflation policies for resource management of transient and on-demand cloud VMs. Experimental evaluation of our deflation system on a Linux cluster shows that microservice-based applications can be deflated by up to 50\% with negligible performance overhead. Our cluster-level deflation policies allow overcommitment levels as high as 50\%, with less than a 1\% decrease in application throughput, and can enable cloud platforms to increase revenue by 30\%.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Cool-Core Cycles and Phoenix

Recent observations show that the star formation rate (SFR) in the {\it Phoenix} cluster's central galaxy is $\sim 500$ M$_\odot$ yr$^{-1}$. Even though {\it Phoenix} is a massive cluster ($M_{200} \approx 2.0\times 10^{15}$ M$_\odot$; $z\approx 0.6$) such a high central SFR is not expected in a scenario in which feedback from an active galactic nucleus (AGN) maintains the intracluster medium (ICM) in a state of rough thermal balance. It has been argued that either AGN feedback saturates in very massive clusters or the central supermassive black hole (SMBH) is too small to produce enough kinetic feedback and hence is unable to quench the catastrophic cooling. In this work, we present an alternate scenario wherein intense short-lived cooling and star formation phases followed by strong AGN outbursts are part of the AGN feedback loop. Using results from a 3D hydrodynamic simulation of a standard cool-core cluster ($M_{200}\sim 7\times10^{14}$ M$_\odot$; $z=0$), scaled to account for differences in mass and redshift, we argue that {\it Phoenix} is at the end of a cooling phase in which an AGN outburst has begun but has not yet arrested core cooling. This state of high cooling rate and star formation is expected to last for $\lesssim$ 100 Myr in {\it Phoenix}.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Modeling The Temporally Constrained Preemptions of Transient Cloud VMs

Transient cloud servers such as Amazon Spot instances, Google Preemptible VMs, and Azure Low-priority batch VMs, can reduce cloud computing costs by as much as $10\times$, but can be unilaterally preempted by the cloud provider. Understanding preemption characteristics (such as frequency) is a key first step in minimizing the effect of preemptions on application performance, availability, and cost. However, little is understood about temporally constrained preemptions---wherein preemptions must occur in a given time window. We study temporally constrained preemptions by conducting a large scale empirical study of Google's Preemptible VMs (that have a maximum lifetime of 24 hours), develop a new preemption probability model, new model-driven resource management policies, and implement them in a batch computing service for scientific computing workloads. Our statistical and experimental analysis indicates that temporally constrained preemptions are not uniformly distributed, but are time-dependent and have a bathtub shape. We find that existing memoryless models and policies are not suitable for temporally constrained preemptions. We develop a new probability model for bathtub preemptions, and analyze it through the lens of reliability theory. To highlight the effectiveness of our model, we develop optimized policies for job scheduling and checkpointing. Compared to existing techniques, our model-based policies can reduce the probability of job failure by more than $2\times$. We also implement our policies as part of a batch computing service for scientific computing applications, which reduces cost by $5\times$ compared to conventional cloud deployments and keeps performance overheads under $3\%$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On modelling CC85 wind in an expanding local box

The Chevalier & Clegg (1985) spherical wind is widely used to model galactic outflows. Efforts are underway to understand the production of multiphase gas in such outflows. Since the important scales for the multiphase gas are very small compared to the size of the wind, it is very useful to have a local model of the cold gas interacting with a volume-filling expanding wind. In this note we explicitly derive the governing hydrodynamic equations in such a local box for the first time. We highlight that the expansion in a wind occurs only in the transverse direction, and this results in anisotropic source terms in the momentum equation. Our equations will be very useful in realistic modelling of the small scale structure of multiphase galactic outflows. [original note does not include an abstract]

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Resolving small-scale cold circumgalactic gas in TNG50

We use the high-resolution TNG50 cosmological magnetohydrodynamical simulation to explore the properties and origin of cold circumgalactic medium (CGM) gas around massive galaxies (M* > 10^11 Msun) at intermediate redshift (z~0.5). We discover a significant abundance of small-scale, cold gas structure in the CGM of &#39;red and dead&#39; elliptical systems, as traced by neutral HI and MgII. Halos can host tens of thousands of discrete absorbing cloudlets, with sizes of order a kpc or smaller. With a Lagrangian tracer analysis, we show that cold clouds form due to strong drho/rho >> 1 gas density perturbations which stimulate thermal instability. These local overdensities trigger rapid cooling from the hot virialized background medium at ~10^7 K to radiatively inefficient ~10^4 K clouds, which act as cosmologically long-lived, &#39;stimulated cooling&#39; seeds in a regime where the global halo does not satisfy the classic tcool/tff < 10 criterion. Furthermore, these small clouds are dominated by magnetic rather than thermal pressure, with plasma beta << 1, suggesting that magnetic fields may play an important role. The number and total mass of cold clouds both increase with resolution, and the ~8x10^4 Msun cell mass of TNG50 enables the ~few hundred pc, small-scale CGM structure we observe to form. Finally, we make a preliminary comparison against observations from the COS-LRG, LRG-RDR, COS-Halos, and SDSS LRG surveys. We broadly find that our recent, high-resolution cosmological simulations produce sufficiently high covering fractions of extended, cold gas as observed to surround massive galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Turbulence in stratified atmospheres: implications for the intracluster medium

The gas motions in the intracluster medium (ICM) are governed by stratified turbulence. Stratified turbulence is fundamentally different from Kolmogorov (isotropic, homogeneous) turbulence; kinetic energy not only cascades from large to small scales, but it is also converted into buoyancy potential energy. To understand the density and velocity fluctuations in the ICM, we conduct high-resolution ($1024^2\times 1536$ grid points) hydrodynamical simulations of subsonic turbulence (with rms Mach number $\mathcal{M}\approx 0.25$) and different levels of stratification, quantified by the Richardson number $\mathrm{Ri}$, from $\mathrm{Ri}=0$ (no stratification) to $\mathrm{Ri}=13$ (strong stratification). We quantify the density, pressure and velocity fields for varying stratification because observational studies often use surface brightness fluctuations to infer the turbulent gas velocities of the ICM. We find that the standard deviation of the logarithmic density fluctuations ($σ_s$), where $s=\ln(ρ/\left<ρ(z)\right>)$, increases with $\mathrm{Ri}$. For weakly stratified subsonic turbulence ($\mathrm{Ri}\lesssim10$, $\mathcal{M}<1$), we derive a new $σ_s$--$\mathcal{M}$--$\mathrm{Ri}$ relation, $σ_s^2=\ln(1+b^2\mathcal{M}^4+0.09\mathcal{M}^2\mathrm{Ri}H_P/H_S)$, where $b=1/3$--$1$ is the turbulence driving parameter, and $H_P$ and $H_S$ are the pressure and entropy scale heights respectively. We further find that the power spectrum of density fluctuations, $P(ρ_k/\left<ρ\right>)$, increases in magnitude with increasing $\mathrm{Ri}$, whereas the velocity power spectrum is invariant. Thus, the ratio between density and velocity power spectra strongly depends on $\mathrm{Ri}$. Pressure fluctuations, on the other hand, are independent of stratification and only depend on $\mathcal{M}$.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Radio halos of star forming galaxies

We study the synchrotron radio emission from extra-planar regions of star forming galaxies. We use ideal magneto-hydrodynamical (MHD) simulations of a rotating Milky Way-type disk galaxy with distributed star formation sites for three star formation rates (SFRs) (0.3, 3, 30 M$_{\odot}$ yr$^{-1}$). From our simulations, we see emergence of galactic-scale magnetised outflows, carrying gas from the disk. We compare the morphology of the outflowing gas with hydrodynamic (HD) simulations. We look at the spatial distribution of magnetic field in the outflows. Assuming that a certain fraction of gas energy density is converted into cosmic ray energy density, and using information about the magnetic field, we obtain synchrotron emissivity throughout the simulation domain. We generate the surface brightness maps at a frequency of 1.4 GHz. The outflows are more extended in the vertical direction than radial and hence have an oblate shape. We further find that the matter right behind the outer shock, shines brighter in these maps than that above or below. To understand whether this feature can be observed, we produce vertical intensity profiles. We convolve the vertical intensity profile with the typical beam sizes of radio telescopes, for a galaxy located at 10 Mpc (similar to NGC 891) in order to estimate the radio scale height to compare with observations. We find that for our SFRs this feature will lie below the RMS noise limit of instruments. The radio scale height is found to be $\sim 300-1200$ pc , depending on the resolution of the telescope. We relate the advection speed of the outer shock with the surface density of star formation as $\rm{v}_{\rm adv} \propto Σ_{\rm SFR}^{0.3}$ which is consistent with earlier observations and analytical estimates.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Realistic modeling of wind and supernovae shocks in star clusters: addressing ${\rm ^{22}Ne/^{20}Ne}$ and other problems in Galactic cosmic rays

Cosmic ray (CR) sources leave signatures in the isotopic abundances of CRs. Current models of Galactic CRs that consider supernovae (SNe) shocks as the main sites of particle acceleration cannot satisfactorily explain the higher ${\rm ^{22}Ne/^{20}Ne}$ ratio in CRs compared to the interstellar medium. Although stellar winds from massive stars have been invoked, their contribution relative to SNe ejecta has been taken as a free parameter. Here we present a theoretical calculation of the relative contributions of wind termination shocks (WTSs) and SNe shocks in superbubbles, based on the hydrodynamics of winds in clusters, the standard stellar mass function, and stellar evolution theory. We find that the contribution of WTSs towards the total CR production is at least $25\%$, which rises to $\gtrsim 50\%$ for young ($\lesssim 10$ Myr) clusters, and explains the observed $^{22}{\rm Ne}/^{20} {\rm Ne}$ ratio. We argue that since the progenitors of apparently isolated supernovae remnants (SNRs) are born in massive star clusters, both WTS and SNe shocks can be integrated into a combined scenario of CRs being accelerated in massive clusters. This scenario is consistent with the observed ratio of SNRs to $γ$-ray bright ($L_γ\gtrsim 10^{35}$ erg s$^{-1}$) star clusters, as predicted by star cluster mass function. Moreover, WTSs can accelerate CRs to PeV energies, and solve other longstanding problems of the standard supernova paradigm of CR acceleration.

0 citations0 reviewsNo code linkedOpen access