Stellar streams that emerge from globular clusters (GCs) are thin stellar structures spread along the orbits of progenitor clusters. Numerical modeling of these streams is essential for understanding their interaction with the host galaxy's mass distribution. Traditional methods are either computationally expensive or oversimplified, motivating us to develop a fast and accurate approach using a particle spray algorithm. By conducting a series of N-body simulations of GCs orbiting a host galaxy, we find that the position and velocity distributions of newly-escaped stream particles are consistent across various GC masses and orbital parameters. Based on these distributions, we develop a new algorithm that avoids computing the detailed internal cluster dynamics by directly drawing tracer particles from these distributions. This algorithm correctly reproduces the action space distribution of stream particles and achieves a 10% accuracy in stream morphology and velocities compared to N-body simulations.
Dark matter (DM) bars that shadow stellar bars have been previously shown to form in idealized simulations of isolated disk galaxies. Here, we show that DM bars commonly occur in barred disk galaxies in the TNG50 cosmological simulation suite, but do not appear in unbarred disk galaxies. Consistent with isolated simulations, DM bars are typically shorter than their stellar counterparts and are $75\%$ weaker as measured by the Fourier $A_2$ moment. DM bars dominate the shape of the inner halo potential and are easily identified in the time series of quadrupolar coefficients. We present two novel methods for measuring the bar pattern speed using these coefficients, and use them to make a measurement of the pattern speed and rotation axis orientation for one sample galaxy located in one of the high-time resolution subboxes of TNG50. The stellar and dark bars in this galaxy remain co-aligned throughout the last 8 Gyr and have identical pattern speeds throughout. Both the pattern speed and rotation axis orientation of the bars evolve considerably over the last 8 Gyr, consistent with torques on the bars due to dynamical friction and gas accretion. While the bar pattern speed generally decreases over the time course, it is seen to increase after mergers. Rather than remaining static in time, the rotation axis displays both precession and nutation possibly caused by torques outside the plane of rotation. We find that the shape of the stellar and DM mass distributions are tightly correlated with the bar pattern speed.
In this paper, we calculate the energy, signal-to-noise ratio, detection range, and angular anisotropy of the matter, matter memory, and neutrino memory gravitational wave (GW) signatures of 21 three-dimensional initially non-rotating core-collapse supernova (CCSN) models carried to late times. We find that inferred energy, signal-to-noise ratio, and detection range are angle-dependent quantities, and that the spread of possible energy, signal-to-noise, and detection ranges across all viewing angles generally increases with progenitor mass. When examining the low-frequency matter memory and neutrino memory components of the signal, we find that the neutrino memory is the most detectable component of a CCSN GW signal, and that DECIGO is best-equipped to detect both matter memory and neutrino memory. Moreover, we find that the polarization angle between the $h_+$ and $h_{\times}$ strains serves as a unique identifier of matter and neutrino memory. Finally, we develop a galactic density- and stellar mass-weighted formalism to calculate the rate at which we can expect to detect CCSN GW signals with Advanced LIGO. When considering only the matter component of the signal, the aLIGO detection rate is around 65$\%$ of the total galactic supernova rate, but increases to 90$\%$ when incorporating the neutrino memory component. We find that all future detectors (ET, CE, DECIGO) will be able to detect CCSN GW signals from the entire galaxy, and for the higher-mass progenitors even into the local group of galaxies.
Short Large Amplitude Magnetic Structures (SLAMS) are frequently detected during spacecraft crossings over the Earth bow shock. We investigate the existence of such structures at astrophysical shocks, where they could result from the steepening of cosmic-ray (CR) driven waves. Using kinetic particle-in-cell simulations, we study the growth of SLAMS and the appearance of associated transient shocks in the upstream region of quasi-parallel, non-relativistic, high-Mach number collisionless shocks. We find that high-energy CRs significantly enhance the transverse magnetic field within SLAMS, producing highly inclined field lines. As SLAMS are advected towards the shock, these fields lines form an intermittent superluminal configuration which traps magnetized electrons at fast shocks. Due to their oscillatory nature, SLAMS are periodically separated by subluminal gaps with lower transverse magnetic field strength. In these regions, electrons diffuse and accelerate by bouncing between the shock and the approaching SLAMS region through a mechanism that we call quasi-periodic shock acceleration (QSA). We analytically derive the distribution of electrons accelerated via QSA, $f(p)\sim p^{[-4.7,-5.7]}$, which agrees well with the simulation spectra. We find that the electron power law remains steep until the end of our longest runs, providing a possible explanation for the steep electron spectra observed at least up to GeV energies in young and fast supernova remnants.
arXiv:1904.02180 by other authors