共找到 20 条结果
Non-radiative shocks accelerate particles and heat astrophysical plasmas. While supernova remnants are the most well-studied example, neutron star (NS) bow shocks are also non-radiative and Balmer-dominated. NS bow shocks are likely ubiquitous in the interstellar medium due to their large speeds imparted at birth, and they are thought to be a discrete source population contributing to the Galactic cosmic ray spectrum. To date, nine NS bow shocks have been directly observed in H$α$ images. Most of these shocks have been characterized using narrowband H$α$ imaging and slit spectroscopy, which do not resolve the multi-component velocity structure of the shocks and their spatial geometry. Here we present integral field spectroscopy of three NS bow shocks: J0742$-$2822, J1741$-$2054, and J2225$+$6535 (the Guitar Nebula). We measure the shock properties simultaneously in four dimensions: the 2D projected shock morphology, the radial velocity structure, and the H$α$ flux. The broad-to-narrow line ratio ($I_{\rm b}/I_{\rm n}$) is inferred from radial velocity profiles, and for J1741$-$2054 the narrow line is detected in multiple regions of the shock. The inferred line ratios and widths sug
Cosmic rays are charged particles that are accelerated to relativistic speeds by astrophysical shocks. Numerical models have been successful in confirming the acceleration process for (quasi-)parallel shocks, which have the magnetic field aligned with the direction of the shock motion. However, the process is less clear when it comes to (quasi-)perpendicular shocks, where the field makes a large angle with the shock-normal. For such shocks, the angle between the magnetic field and flow ensures that only highly energetic particles can travel upstream at all, reducing the upstream current. This process is further inhibited for relativistic shocks, since the shock can become superluminal when the required particle velocity exceeds the speed of light, effectively inhibiting any upstream particle flow. In order to determine whether such shocks can accelerate particles, we use the particle-in-cell (PIC) method to determine what fraction of particles gets reflected initially at the shock. We then use this as input for a new simulation that combines the PIC method with grid-based magnetohydrodynamics to follow the acceleration (if any) of the particles over a larger time-period in a two-di
Residual energy quantifies the difference in energy between velocity and magnetic field fluctuations in a plasma. Recent observational evidence highlights that fast-mode interplanetary shock waves have positive residual energy, in sharp contrast to the negative residual energy of the turbulence and magnetic structures that constitute the vast majority of fluctuation power in the solar wind at magnetohydrodynamic (MHD) inertial scales. In this work, we apply the Rankine-Hugoniot conditions to derive an equation for the residual energy of an MHD shock jump as a function of the shock angle, density compression ratio and Alfvén Mach number upstream of the shock. An equation for the cross helicity is similarly derived. The residual energy equation gives only positive values for super-Alfvénic (i.e. fast-mode) shocks. The residual energy and cross helicity of slow-mode shocks and tangential, contact and rotational discontinuities are also determined. A simplified form of the residual energy equation applicable to perpendicular shocks has been verified against residual energy values directly estimated from observations of 141 interplanetary shocks; the equation is found to match well with
Supernova shocks into dense molecular cores in IC 443 (clumps B, C, and G) and 3C 391 were observed using the Stratospheric Observatory for Infrared Astronomy and complemented by archival data from the Herschel Space Observatory. The pure rotational transitions 0-0 S(1) and S(5) of H2, and the ground-state 110-101 transition of H2O, are all broadened, arising from molecules that survive the passage of the shock front. Theoretical models from the Paris-Durham shock code were analyzed to generate synthetic profiles that approximately match the observations. The observations can be fit with two shock conditions, which approximate the range of densities in the pre-shock molecular cloud. The width and brightness of the S(5) lines require shocks into gas with a density of order 2,000 cm-3, into which the IC 443 blast wave drives shocks with speed 60 km/s. The brightness and narrower width of the S(1) lines requires different shocks, into gas with density of order 10^5 cm-3, with shock speeds of 10 km/s. The H2O velocity distribution is also consistent with these shocks. The existence of shocks into dense gas shows that the bright shocked clumps in IC~443 were prestellar cores. It is unli
In the aftermath of the COVID-19 pandemic, empirical data have revealed that large-scale health crises not only cause immediate disruptions in mortality dynamics but also have persistent effects that may last for several years. Existing mortality models largely assume that mortality shocks are transitory and overlook how their effects can be long-lasting and heterogeneous across age groups and causes of death. In response to this limitation, we propose a novel stochastic mortality model that captures age- and cause-specific long-lasting effects of mortality jumps through a gamma-density-like decay function, estimated via a customized conditional maximum likelihood algorithm. Applying the model to recent U.S. mortality data, we reveal divergent persistence patterns across demographic groups and provide key insights into the tail risk profiles of life insurance and annuity products. Our scenario-based analyses further show that neglecting persistent shock effects can lead to suboptimal hedging, while the proposed model enables what-if testing to analyze such effects under potential future health crises.
We propose a method to learn the nonlinear impulse responses to structural shocks using neural networks, and apply it to uncover the effects of US financial shocks. The results reveal substantial asymmetries with respect to the sign of the shock. Adverse financial shocks have powerful effects on the US economy, while benign shocks trigger much smaller reactions. Instead, with respect to the size of the shocks, we find no discernible asymmetries.
Context: Cosmic ray acceleration in galaxy clusters is still an ongoing puzzle, with relativistic electrons forming radio relics at merger shocks and emitting synchrotron radiation. In the present work we perform hybrid-kinetic simulations in a range of various quasi-perpendicular foreshock conditions, including plasma beta, magnetic obliquity, and the shock Mach number. Aims: We study the ion kinetic physics, responsible for the shock structure and wave turbulence, that in turn affects the particle acceleration processes. Methods: We apply a recently developed generalized fluid-particle hybrid numerical code that can combine fluid modeling for both kinetic ions and fluid electrons. The model utilizes the exact form of the generalized Ohm's law, allowing for arbitrary choice of mass and energy densities, as well as the charge-to-mass ratio of the kinetic species. Results: We show that the properties of ion-driven multi-scale magnetic turbulence in merger shocks are in agreement with the ion structures observed in PIC simulations. In typical shocks with the sonic Mach number $M_s=3$, the magnetic structures and shock front density ripples grow and saturate at wavelengths reaching ap
Nonthermal particles can be produced due to incomplete thermalization at collisionless shocks and further accelerated to very high energies via diffusive shock acceleration. In a previous study we explored the cosmic ray (CR) acceleration at cosmic shocks through numerical simulations of CR modified, quasi-parallel shocks in 1D plane-parallel geometry with the physical parameters relevant for the shocks emerging in the large scale structure formation of the universe (Kang & Jones 2002). In order to consider the CR acceleration at shocks with a broader range of physical properties, in this contribution we present additional simulations. For a given Mach number, shocks with higher speeds accelerate CRs faster with a greater number of particles. However, two shocks with a same Mach number but with different shock speeds evolve qualitatively similarly when the results are presented in terms of diffusion length and time scales. Therefore, the time asymptotic value for the fraction of shock kinetic energy transferred to CRs is mainly controlled by shock Mach number rather than shock speed. We present the ``CR energy ratio'' for a wide range of shock parameters and for $ε=0.2-0.3$ at
We present results from a fully relativistic Monte Carlo simulation of diffusive shock acceleration (DSA) in unmodified shocks. The computer code uses a single algorithmic sequence to smoothly span the range from nonrelativistic speeds to fully relativistic shocks of arbitrary obliquity, providing a powerful consistency check. While known results are obtained for nonrelativistic and ultra-relativistic parallel shocks, new results are presented for the less explored trans- relativistic regime and for oblique, fully relativistic shocks. We find, for a wide trans-relativistic range extending to shock Lorentz factors >30, that the particle spectrum produced by DSA varies strongly from the canonical f(p) proportional to p^{-4.23} spectrum known to result in ultra-relativistic shocks. Trans- relativistic shocks may play an important role in gamma-ray bursts and other sources and most relativistic shocks will be highly oblique.
Extending our earlier work on Lax-type shocks of systems of conservation laws, we establish existence and stability of curved multidimensional shock fronts in the vanishing viscosity limit for general Lax- or undercompressive-type shock waves of nonconservative hyperbolic systems with parabolic regularization. The hyperbolic equations may be of variable multiplicity and the parabolic regularization may be of ``real'', or partially parabolic, type. We prove an existence result for inviscid nonconservative shocks that extends to multidimensional shocks a one-dimensional result of X. Lin proved by quite different methods. In addition, we construct families of smooth viscous shocks converging to a given inviscid shock as viscosity goes to zero, thereby justifying the small viscosity limit for multidimensional nonconservative shocks. In our previous work on shocks we made use of conservative form especially in parts of the low frequency analysis. Thus, most of the new analysis of this paper is concentrated in this area. By adopting the more general nonconservative viewpoint, we are able to shed new light on both the viscous and inviscid theory. For example, we can now provide a clearer
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.
The intracluster medium (ICM) is expected to experience on average about three passages of weak shocks with low sonic Mach numbers, $M\lesssim 3$, during the formation of galaxy clusters. Both protons and electrons could be accelerated to become high energy cosmic rays (CRs) at such ICM shocks via diffusive shock acceleration (DSA). We examine the effects of DSA by multiple shocks on the spectrum of accelerated CRs by including {\it in situ} injection/acceleration at each shock, followed by repeated re-acceleration at successive shocks in the test-particle regime. For simplicity, the accelerated particles are assumed to undergo adiabatic decompression without energy loss and escape from the system, before they encounter subsequent shocks. We show that in general the CR spectrum is flattened by multiple shock passages, compared to a single episode of DSA, and that the acceleration efficiency increases with successive shock passages. However, the decompression due to the expansion of shocks into the cluster outskirts may reduce the amplification and flattening of the CR spectrum by multiple shock passages. The final CR spectrum behind the last shock is determined by the accumulated e
We study the impact of fiscal revenue shocks on local fiscal policy. We focus on the very volatile revenues from the immovable property gains tax in the canton of Zurich, Switzerland, and analyze fiscal behavior following large and rare positive and negative revenue shocks. We apply causal machine learning strategies and implement the post-double-selection LASSO estimator to identify the causal effect of revenue shocks on public finances. We show that local policymakers overall predominantly smooth fiscal shocks. However, we also find some patterns consistent with fiscal conservatism, where positive shocks are smoothed, while negative ones are mitigated by spending cuts.
This is the second in a series of papers on the effects of dust on multifluid, MHD shock waves in weakly ionized molecular gas. We investigate the influence of dust on the critical shock speed, v_crit, above which C shocks cease to exist. Chernoff showed that v_crit cannot exceed the grain magnetosound speed, v_gms, if dust grains are dynamically well coupled to the magnetic field. We present numerical simulations of steady shocks where the grains may be well- or poorly coupled to the field. We use a time-dependent, multifluid MHD code that models the plasma as a system of interacting fluids: neutral particles, ions, electrons, and various ``dust fluids'' comprised of grains with different sizes and charges. Our simulations include grain inertia and grain charge fluctuations but to highlight the essential physics we assume adiabatic flow, single-size grains, and neglect the effects of chemistry. We show that the existence of a phase speed v_phi does not necessarily mean that C shocks will form for all shock speeds v_s less than v_phi. When the grains are weakly coupled to the field, steady, adiabatic shocks resemble shocks with no dust: the transition to J type flow occurs at v_cri
The interstellar medium (ISM) is constantly evolving due to unremitting injection of energy in various forms. Energetic radiation transfers energy to the ISM: from the UV photons, emitted by the massive stars, to X- and $γ$-ray ones. Cosmic rays are another source of energy. Finally, mechanical energy is injected through shocks or turbulence. Shocks are ubiquitous in the interstellar medium of galaxies. They are associated to star formation (through jets and bipolar outflows), life (via stellar winds), and death (in AGB stellar winds or supernovae explosion). The dynamical processes leading to the formation of molecular clouds also generate shocks where flows of interstellar matter collide. Because of their ubiquity, the study of interstellar shocks is also a useful probe to the other mechanisms of energy injection in the ISM. This study must be conducted in order to understand the evolution of the ISM as a whole, and to address various questions: what is the peculiar chemistry associated to shocks, and what is their contribution to the cycle of matter in galaxies ? What is the energetic impact of shocks on their surroundings on various scales, and hence what is the feedback of sta
We present the first results for the sputtering of grain mantles and cores obtained with self-consistent multifluid hydromagnetic models of C-type shocks propagating through dusty media. The threshold shock speed for mantle sputtering is about 10 km/s and is independent of density. The mantles are completely vapourised in shocks with speeds of 20-25 km/s. At such shock speeds core sputtering commences and gas-phase SiO forms. Core destruction is not total in any C-type shock because grains are not completely destroyed in shocks with speeds near the minimum speeds at which J-type shocks appear. Due to the density-dependence of the critical shock speed for this transition, higher gas-phase SiO fractional abundances are produced behind shocks propagating in lower density gas. For shock speeds near the threshold speeds for both core and mantle sputtering, sputtering is much greater for shock velocities at smaller angles relative to the upstream magnetic field. At higher shock speeds, the angular variation is still present but less pronounced.
We discuss 3D simulations of relativistic collisionless shocks in electron-positron pair plasmas using the particle-in-cell (PIC) method. The shock structure is mainly controlled by the shock's magnetization ("sigma" parameter). We demonstrate how the structure of the shock varies as a function of sigma for perpendicular shocks. At low magnetizations the shock is mediated mainly by the Weibel instability which generates transient magnetic fields that can exceed the initial field. At larger magnetizations the shock is dominated by magnetic reflections. We demonstrate where the transition occurs and argue that it is impossible to have very low magnetization collisionless shocks in nature (in more than one spatial dimension). We further discuss the acceleration properties of these shocks, and show that higher magnetization perpendicular shocks do not efficiently accelerate nonthermal particles in 3D. Among other astrophysical applications, this may pose a restriction on the structure and composition of gamma-ray bursts and pulsar wind outflows.
Collisionless shocks occur in various fields of physics. In the context of space and astrophysics they have been investigated for many decades. However, a thorough understanding of shock formation and particle acceleration is still missing. Collisionless shocks can be distinguished into electromagnetic and electrostatic shocks. Electromagnetic shocks are of importance mainly in astrophysical environments and they are mediated by the Weibel or filamentation instability. In such shocks, charged particles gain energy by diffusive shock acceleration. Electrostatic shocks are characterized by a strong electrostatic field, which leads to electron trapping. Ions are accelerated by reflection from the electrostatic potential. Shock formation and particle acceleration will be discussed in theory and simulations.
Perpendicular relativistic ($γ_0=10$) shocks in magnetized pair plasmas are investigated using two dimensional Particle-in-Cell simulations. A systematic survey, from unmagnetized to strongly magnetized shocks, is presented accurately capturing the transition from Weibel-mediated to magnetic-reflection-shaped shocks. This transition is found to occur for upstream flow magnetizations $10^{-3}<σ<10^{-2}$ at which a strong perpendicular net current is observed in the precursor, driving the so-called current-filamentation instability. The global structure of the shock and shock formation time are discussed. The MHD shock jump conditions are found in good agreement with the numerical results, except for $10^{-4} < σ< 10^{-2}$ where a deviation up to 10\% is observed. The particle precursor length converges toward the Larmor radius of particles injected in the upstream magnetic field at intermediate magnetizations. For $σ>10^{-2}$, it leaves place to a purely electromagnetic precursor following from the strong emission of electromagnetic waves at the shock front. Particle acceleration is found to be efficient in weakly magnetized perpendicular shocks in agreement with prev
We study the structure of shocks in clumpy media, using a multifluid formalism. As expected, shocks broaden as they weaken: for sufficiently weak shocks, no viscous subshock appears in the structure. This has significant implications for the survival of dense clouds in regions overrun by shocks in a wide range of astrophysical circumstances, from planetary nebulae to the nuclei of starburst galaxies.