Search for: All records

We revisit the double adiabatic evolution equations and extend them to the relativistic and ultrarelativistic regimes. We analytically solve the relativistic, time-dependent drift kinetic equation for a homogeneous, magnetized, collisionless plasma and obtain a solution explicitly dependent on the magnetic field and density variations. In the case of an initial relativistic Maxwellian distribution, a natural extension to an anisotropic Maxwell-Jüttner is obtained. We calculate the moments of this time-dependent solution and obtain analytical expressions for the evolution of the perpendicular and parallel pressures in the ultrarelativistic case. We numerically solve the moment equations in the relativistic case and obtain general expressions for the double-adiabatic equations in this regime. We confirm our results using fully kinetic particle-in-cell simulations of shearing and compressing boxes. Our findings can be readily applied to relativistic species including cosmic-rays and electron-positron pairs, present in astrophysical plasmas like pulsar wind nebulae, astrophysical jets, black hole accretion flows, and Van Allen radiation belts. 

Abstract Stars form in clusters from the gravitational collapse of giant molecular clouds, which is opposed by a variety of physical processes, including stellar feedback. The interplay between these processes determines the star formation rate of the clouds. To study how feedback controls star formation, we use a numerical framework that is optimized to simulate star cluster formation and evolution. This framework, calledTorch, combines the magnetohydrodynamical code FLASH with N-body and stellar evolution codes in the Astrophysical Multipurpose Software Environment (AMUSE).Torchincludes stellar feedback from ionizing and nonionizing radiation, stellar winds, and supernovae, but, until now, did not include protostellar jets. We present our implementation of protostellar jet feedback within theTorchframework and describe its free parameters. We then demonstrate our new module by comparing cluster formation simulations with and without jets. We find that the inclusion of protostellar jets slows star formation, even in clouds of up toM = 2 × 10⁴M_⊙. We also find that the star formation rate of our lower-mass clouds (M = 5 × 10³M_⊙) is strongly affected by both the inclusion of protostellar jets and the chosen jet parameters, including the jet lifetime and injection velocity. We follow the energy budget for each simulation and find that the inclusion of jets systematically increases the kinetic energy of the gas at early times. The implementation of protostellar jet feedback inTorchopens new areas of investigation regarding the role of feedback in star cluster formation and evolution. 

Abstract Collisionless low-Mach-number shocks are abundant in astrophysical and space plasma environments, exhibiting complex wave activity and wave–particle interactions. In this paper, we present 2D Particle-in-Cell (PIC) simulations of quasi-perpendicular nonrelativistic (v_sh≈ (5500–22000) km s⁻¹) low-Mach-number shocks, with a specific focus on studying electrostatic waves in the shock ramp and precursor regions. In these shocks, an ion-scale oblique whistler wave creates a configuration with two hot counterstreaming electron beams, which drive unstable electron acoustic waves (EAWs) that can turn into electrostatic solitary waves (ESWs) at the late stage of their evolution. By conducting simulations with periodic boundaries, we show that the EAW properties agree with linear dispersion analysis. The characteristics of ESWs in shock simulations, including their wavelength and amplitude, depend on the shock velocity. When extrapolated to shocks with realistic velocities (v_sh≈ 300 km s⁻¹), the ESW wavelength is reduced to one-tenth of the electron skin depth and the ESW amplitude is anticipated to surpass that of the quasi-static electric field by more than a factor of 100. These theoretical predictions may explain a discrepancy, between PIC and satellite measurements, in the relative amplitude of high- and low-frequency electric field fluctuations. 

Abstract General relativistic magnetohydrodynamic (GRMHD) simulations of black hole tilted disks—where the angular momentum of the accretion flow at large distances is misaligned with respect to the black hole spin—commonly display standing shocks within a few to tens of gravitational radii from the black hole. In GRMHD simulations of geometrically thick, optically thin accretion flows, applicable to low-luminosity sources like Sgr A* and M87*, the shocks have transrelativistic speed, moderate plasma beta (the ratio of ion thermal pressure to magnetic pressure isβ_pi1∼ 1–8), and low sonic Mach number (the ratio of shock speed to sound speed isM_s∼ 1–6). We study such shocks with 2D particle-in-cell simulations, and we quantify the efficiency and mechanisms of electron heating for the special case of preshock magnetic fields perpendicular to the shock direction of propagation. We find that the postshock electron temperatureT_e2exceeds the adiabatic expectationT_e2,adby an amount $T_{\mathrm{e}2}/T_{\mathrm{e}2,\mathrm{ad}}-1\simeq 0.0016M_s^{3.6}$ , nearly independent of the plasma beta and of the preshock electron-to-ion temperature ratioT_e1/T_i1, which we vary from 0.1 to unity. We investigate the heating physics forM_s∼ 5–6 and find that electron superadiabatic heating is governed by magnetic pumping atT_e1/T_i1= 1, whereas heating byB-parallel electric fields (i.e., parallel to the local magnetic field) dominates atT_e1/T_i1= 0.1. Our results provide physically motivated subgrid prescriptions for electron heating at the collisionless shocks seen in GRMHD simulations of black hole accretion flows. 

The mode of star formation that results in the formation of globular clusters and young massive clusters is difficult to constrain through observations. We present models of massive star cluster formation using the TORCHframework, which uses the Astrophysical MUltipurpose Software Environment (AMUSE) to couple distinct multi-physics codes that handle star formation, stellar evolution and dynamics, radiative transfer, and magnetohydrodynamics. We upgraded TORCHby implementing the N-body code PETAR, thereby enabling TORCHto handle massive clusters forming from 10⁶M_⊙clouds with ≥10⁵individual stars. We present results from TORCHsimulations of star clusters forming from 10⁴,  10⁵, and 10⁶M_⊙turbulent spherical gas clouds (named M4, M5, M6) of radiusR= 11.7 pc. We find that star formation is highly efficient and becomes more so at a higher cloud mass and surface density. For M4, M5, and M6 with initial surface densities 2.325 × 10^1,2,3M_⊙pc⁻², after a free-fall time oft_ff= 6.7,2.1,0.67 Myr, we find that ∼30%, 40%, and 60% of the cloud mass has formed into stars, respectively. The end of simulation-integrated star formation efficiencies for M4, M5, and M6 areϵ_⋆ = M_⋆/M_cloud = 36%, 65%, and 85%. Observations of nearby clusters similar in mass and size to M4 have instantaneous star formation efficiencies ofϵ_inst ≤ 30%, which is slightly lower than the integrated star formation efficiency of M4. The M5 and M6 models represent a different regime of cluster formation that is more appropriate for the conditions in starburst galaxies and gas-rich galaxies at high redshift, and that leads to a significantly higher efficiency of star formation. We argue that young massive clusters build up through short efficient bursts of star formation in regions that are sufficiently dense (Σ ≥ 10²M_⊙pc⁻²) and massive (M_cloud≥ 10⁵M_⊙). In such environments, stellar feedback from winds and radiation is not strong enough to counteract the gravity from gas and stars until a majority of the gas has formed into stars. 

Abstract In galaxy clusters, the intracluster medium (ICM) is expected to host a diffuse, long-lived, and invisible population of “fossil” cosmic-ray electrons (CRe) with 1–100 MeV energies. These CRe, if reaccelerated by 100× in energy, can contribute synchrotron luminosity to cluster radio halos, relics, and phoenices. Reacceleration may be aided by CRe scattering upon the ion-Larmor-scale waves that spawn when ICM is compressed, dilated, or sheared. We study CRe scattering and energy gain due to ion cyclotron (IC) waves generated by continuously driven compression in 1D fully kinetic particle-in-cell simulations. We find that pitch-angle scattering of CRe by IC waves induces energy gain via magnetic pumping. In an optimal range of IC-resonant momenta, CRe may gain up to ∼10%–30% of their initial energy in one compression/dilation cycle with magnetic field amplification ∼3–6×, assuming adiabatic decompression without further scattering and averaging over initial pitch angle. 

Abstract Turbulence driven by active galactic nuclei activity, cluster mergers, and galaxy motion constitutes an attractive energy source for heating the intracluster medium (ICM). How this energy dissipates into the ICM plasma remains unclear, given its low collisionality and high magnetization (precluding viscous heating by Coulomb processes). Kunz et al. proposed a viable heating mechanism based on the anisotropy of the plasma pressure under ICM conditions. The present paper builds upon that work and shows that particles can be heated by large-scale turbulent fluctuations via magnetic pumping. We study how the anisotropy evolves under a range of forcing frequencies, what waves and instabilities are generated, and demonstrate that the particle distribution function acquires a high-energy tail. For this, we perform particle-in-cell simulations where we periodically vary the mean magnetic field B ( t ). When B ( t ) grows (dwindles), a pressure anisotropy P ⊥ > P ∥ ( P ⊥ < P ∥ ) builds up ( P ⊥ and P ∥ are, respectively, the pressures perpendicular and parallel to B ( t )). These pressure anisotropies excite mirror ( P ⊥ > P ∥ ) and oblique firehose ( P ∥ > P ⊥ ) instabilities, which trap and scatter the particles, limiting the anisotropy, and providing a channel to heat the plasma. The efficiency of this mechanism depends on the frequency of the large-scale turbulent fluctuations and the efficiency of the scattering the instabilities provide in their nonlinear stage. We provide a simplified analytical heating model that captures the phenomenology involved. Our results show that this process can be relevant in dissipating and distributing turbulent energy at kinetic scales in the ICM. 

ABSTRACT Protoplanetary discs spend their lives in the dense environment of a star-forming region. While there, they can be affected by nearby stars through external photoevaporation and dynamic truncations. We present simulations that use the amuse framework to couple the torch model for star cluster formation from a molecular cloud with a model for the evolution of protoplanetary discs under these two environmental processes. We compare simulations with and without extinction of photoevaporation-driving radiation. We find that the majority of discs in our simulations are considerably shielded from photoevaporation-driving radiation for at least 0.5 Myr after the formation of the first massive stars. Radiation shielding increases disc lifetimes by an order of magnitude and can let a disc retain more solid material for planet formation. The reduction in external photoevaporation leaves discs larger and more easily dynamically truncated, although external photoevaporation remains the dominant mass-loss process. Finally, we find that the correlation between disc mass and projected distance to the most massive nearby star (often interpreted as a sign of external photoevaporation) can be erased by the presence of less massive stars that dominate their local radiation field. Overall, we find that the presence and dynamics of gas in embedded clusters with massive stars is important for the evolution of protoplanetary discs. 

ABSTRACT We perform simulations of star cluster formation to investigate the morphological evolution of embedded star clusters in the earliest stages of their evolution. We conduct our simulations with Torch, which uses the Amuse framework to couple state-of-the-art stellar dynamics to star formation, radiation, stellar winds, and hydrodynamics in Flash. We simulate a suite of 104 M⊙ clouds at 0.0683 pc resolution for ∼2 Myr after the onset of star formation, with virial parameters α = 0.8, 2.0, 4.0 and different random samplings of the stellar initial mass function and prescriptions for primordial binaries. Our simulations result in a population of embedded clusters with realistic morphologies (sizes, densities, and ellipticities) that reproduce the known trend of clouds with higher initial α having lower star formation efficiencies. Our key results are as follows: (1) Cluster mass growth is not monotonic, and clusters can lose up to half of their mass while they are embedded. (2) Cluster morphology is not correlated with cluster mass and changes over ∼0.01 Myr time-scales. (3) The morphology of an embedded cluster is not indicative of its long-term evolution but only of its recent history: radius and ellipticity increase sharply when a cluster accretes stars. (4) The dynamical evolution of very young embedded clusters with masses ≲1000 M⊙ is dominated by the overall gravitational potential of the star-forming region rather than by internal dynamical processes such as two- or few-body relaxation.