Jet Model Research Articles

Model of astrophysical jets in magnetic fields of laser relativistic plasmas

A brief review of the results of experimental modeling of cosmic jets in superstrong magnetic fields of laser relativistic plasma is given. It is noted that the development of cyclotron instability with the generation of cyclotron radiation plays a key role in a number of processes in a plasma with a magnetic field — self-localization of plasma in the form of solitons, conversion of rotational motion of plasma into translational motion, cyclotron acceleration of charged particles, separation (stratification) of a plasma jet into separate plasma formations.

Dynamics of the plasma ejection structure in laboratory modeling of young star jets at plasma focus facilities

The use of “plasma focus” type facilities, such as PF-3 (Kurchatov Institute), makes it possible to carry out well-controlled and diagnosable laboratory experiments to study laboratory jets with scale parameters close to the jets of young stars. In this paper, we present the results of numerical modeling of plasma outburst propagation in PF-3. A self-consistent configuration was chosen as the initial conditions, which correctly takes into account the internal structure of the jet. This allowed us to obtain a detailed structure of the interaction between the magnetized emission and the ambient gas. Due to the scalability of such a structure, one should expect such a structure from the head shock waves of jets of young stars.

High energy neutrino production in gamma-ray bursts: dependence of the neutrino signal on the jet composition

Heavy nuclei can be synthetized or entrained in gamma-ray bursts (GRBs) with implications on the high-energy neutrino emission. By means of a Monte-Carlo algorithm, we model nuclear cascades and investigate their impact on the neutrino production considering kinetic dominated jets (in the internal shock model, including a dissipative photosphere) as well as Poynting flux dominated jets (for a jet model invoking internal-collision-induced magnetic reconnection and turbulence, ICMART). We find that the ICMART model allows for efficient nuclear cascades leading to an overall larger neutrino fluence than in the other two jet models. The survival of nuclei and inefficient nuclear cascades lead to an overall reduction of the neutrino fluence up to one order of magnitude. However, if nuclei are disintegrated, the neutrino fluence may be comparable to the one emitted from a jet loaded with protons. Exploring the parameter space of jet properties, we conclude that the composition and the bulk Lorentz factor have significant impact on the efficiency of nuclear cascades as well as the spectral shape of the expected neutrino fluence. On the other hand, the neutrino spectral distribution is less sensitive to the power-law index of the accelerated population of protons or heavier nuclei. For what concerns the diffuse emission of neutrinos from GRBs, we find that the uncertainty due to the jet composition can be at most comparable to the one related to the GRB cosmological rate.

Afterglow Linear Polarization Signatures from Shallow GRB Jets: Implications for Energetic GRBs

Gamma-ray bursts (GRBs) are powered by ultrarelativistic jets. The launching sites of these jets are surrounded by dense media, which the jets must cross before they can accelerate and release high-energy emission. Interaction with the medium leads to the formation of a mildly relativistic sheath around the jet, resulting in an angular structure to the jet’s asymptotic Lorentz factor and energy per solid angle, which modifies the afterglow emission. We build a semi-analytical tool to analyze the afterglow light curve and polarization signatures of jets observed from a wide range of viewing angles, and focus on ones with slowly declining energy profiles known as shallow jets. We find overall lower polarization compared to the classical top hat jet model. We provide an analytical expression for the peak polarization degree as a function of the energy profile power-law index, magnetic field configuration, and viewing angle, and show that it occurs near the light-curve break time for all viewers. When applying our tool to GRB 221009A, suspected to originate from a shallow jet, we find that the suggested jet structures for this event agree with the upper limits placed on the afterglow polarization in the optical and X-ray bands. We also find that at early times the polarization levels may be significantly higher, allowing for a potential distinction between different jet structure models and possibly constraining the magnetization in both forward and reverse shocks at that stage.

Structured Jet Model for Multiwavelength Observations of the Jetted Tidal Disruption Event AT 2022cmc

AT 2022cmc is a recently documented tidal disruption event that exhibits a luminous jet, accompanied by fast-declining X-ray and long-lasting radio and millimeter emission. Motivated by the distinct spectral and temporal signatures between the X-ray and radio observations, we propose a multizone model involving relativistic jets with different Lorentz factors. We systematically study the evolution of faster and slower jets in an external density profile, considering the continuous energy injection rate associated with time-dependent accretion rates before and after the mass fallback time. We investigate time-dependent multiwavelength emission from both the forward shock (FS) and reverse shock (RS) regions of the fast and slow jets, in a self-consistent manner. Our analysis demonstrates that the energy injection rate can significantly impact the jet evolution and subsequently influence the lightcurves. We find that the X-ray spectra and lightcurves could be described by electron synchrotron emission from the RS of the faster jet, in which the late-time X-ray upper limits, extending to 400 days after the disruption, could be interpreted as a jet break. Meanwhile, the radio observations can be interpreted as a result of synchrotron emission from the FS region of the slower jet. We also discuss prospects for testing the model with current and future observations.

Multifrequency polarimetry of high-synchrotron peaked blazars probes the shape of their jets

Multifrequency polarimetry is emerging as a powerful probe of blazar jets, especially with the advent of the Imaging X-ray Polarimetry Explorer (IXPE) space observatory. We studied the polarization of high-synchrotron peaked (HSP) blazars, for which both optical and X-ray emission can be attributed to synchrotron radiation from a population of nonthermal electrons. We adopted an axisymmetric stationary force-free jet model in which the electromagnetic fields are determined by the jet shape. When the jet is nearly parabolic, the X-ray polarization degree is ΠX ∼ 15–50%, and the optical polarization degree is ΠO ∼ 5–25%. The polarization degree is strongly chromatic: ΠX/ΠO ∼ 2–9. This chromaticity is due to the softening of the electron distribution at high energies, and is much stronger than for a uniform magnetic field. The electric vector position angle (EVPA) is aligned with the projection of the jet axis on the plane of the sky. These results compare very well with multifrequency polarimetric observations of HSP blazars. When the jet is instead nearly cylindrical, the polarization degree is large and weakly chromatic (we find ΠX ∼ 70% and ΠO ∼ 60%, close to the expected values for a uniform magnetic field). The EVPA is perpendicular to the projection of the jet axis on the plane of the sky. A cylindrical geometry is therefore practically ruled out by current observations. The polarization degree and the EVPA may be less sensitive to the specific particle acceleration process (e.g., magnetic reconnection or shocks) than previously thought.

Dissipation of AGN Jets in a Clumpy Interstellar Medium

Accreting supermassive black holes frequently power jets that interact with the interstellar medium (ISM)/circumgalactic medium, 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 on the models of narrow conical jets interacting with a smooth ambient medium, including the effect of dense clouds, which 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 and cocoon speed, and jet radius. For the first time, we obtain the analytic condition on the ambient medium’s clumpiness 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 with more sophisticated numerical simulations and direct observations of jet–ISM interaction (e.g., quasar J1316+1753), and we 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.

The remarkable microquasar S26: A super-Eddington PeVatron

S26 is an extragalactic microquasar with the most powerful jets ever discovered. They have a kinetic luminosity of $L_ j $. This implies that the accretion power to the black hole should be super-Eddington, of the order of $L_ acc j $. However, the observed X-ray flux of this system indicates an apparent very sub-Eddington accretion luminosity of $L_ X We aim to characterize the nature of S26, explain the system emission, and study the feasibility of super-Eddington microquasars as potential PeVatron sources. We first analyze multi-epoch X-ray observations of S26 obtained with XMM-Newton and model the super-Eddington disk and its wind. We then develop a jet model and study the particle acceleration and radiative processes that occur in shocks generated near the base of the jet and in its terminal region. We find that the discrepancy between the jet and the apparent disk luminosities in S26 is caused by the complete absorption of the disk radiation by the wind ejected from the super-Eddington disk. The non-thermal X-rays are produced near the base of the jet, and the thermal X-rays are emitted in the terminal regions. The radio emission observed with the Australia Telescope Compact Array can be explained as synchrotron radiation produced at the reverse shock in the lobes. We also find that S26 can accelerate protons to PeV energies in both the inner jet and the lobes. The ultra-high energy protons accelerated in the lobes of S26 are injected into the interstellar medium with a total power of $ We conclude that S26 is a super-Eddington microquasar with a dense disk-driven wind that obscures the X-ray emission from the inner disk, and that the supercritical nature of the system allows the acceleration of cosmic rays to PeV energies.

Survey of non-thermal electrons around supermassive black holes through polarization flips

Abstract Optically thick non-thermal synchrotron sources notably produce linear polarization vectors that are parallel to projected magnetic field lines on the observer’s screen, although they are perpendicular in well-known optically thin cases. To elucidate the complex relationship between the vectors and fields, and to investigate the energy and spatial distribution of non-thermal electrons through the images, we perform polarization radiative transfer calculations at submillimeter wavelengths. Here the calculations are based on semi-analytic force-free jet models with non-thermal electrons with a power-law energy distribution. In the calculated images, we find a $90{^\circ }$-flip of linear polarization vectors at the base of counter-side (receding) jet near a black hole, which occurs because of large optical depths for the synchrotron self-absorption effect. The $90{^\circ }$-flip of LP vectors is also seen on the photon ring at a high frequency, since the optical depth along the rays is large there due to the light bending effect. In addition, we see the flip of the sign of circular polarization components on the counter-jet and photon ring. Furthermore, we show that these polarization flips are synthesized with large values in the spectral index map, and also give rise to outstanding features in the Faraday Rotation Measure map. Since the conditions of flipping depend on the magnetic field strength and configuration and the energy distribution of electrons, we can expect that the polarization flips will provide us with observational evidence for the presence of non-thermal electrons around the black hole, and a clue to the magnetic driving mechanism of plasma jets.

Discovery of Limb Brightening in the Parsec-scale Jet of NGC 315 through Global Very Long Baseline Interferometry Observations and Its Implications for Jet Models

We report the first observation of the nearby giant radio galaxy NGC 315 using a global very long baseline interferometry (VLBI) array consisting of 22 radio antennas located across five continents, including high-sensitivity stations, at 22 GHz. Utilizing the extensive u v-coverage provided by the array, coupled with the application of a recently developed superresolution imaging technique based on the regularized maximum-likelihood method, we were able to transversely resolve the NGC 315 jet at parsec scales for the first time. Previously known for its central ridge-brightened morphology at similar scales in former VLBI studies, the jet now clearly exhibits a limb-brightened structure. This finding suggests an inherent limb brightening that was not observable before due to limited angular resolution. Considering that the jet is viewed at an angle of ∼50°, the observed limb brightening is challenging to reconcile with the magnetohydrodynamic models and simulations, which predict that the Doppler-boosted jet edges should dominate over the nonboosted central layer. The conventional jet model that proposes a fast spine and a slow sheath with uniform transverse emissivity may pertain to our observations. However, in this model, the relativistic spine would need to travel at speeds of Γ ≳ 6.0–12.9 along the deprojected jet distance of (2.3–10.8) × 103 gravitational radii from the black hole. We propose an alternative scenario that suggests higher emissivity at the jet boundary layer, resulting from more efficient particle acceleration or mass loading onto the jet edges, and consider prospects for future observations with even higher angular resolution.

Mapping the Distribution of the Magnetic Field Strength along the NGC 315 Jet

We study magnetic field strengths along the jet in NGC 315. First, we estimated the angular velocity of rotation in the jet magnetosphere by comparing the measured velocity profile of NGC 315 with the magnetohydrodynamic jet model proposed by Tomimatsu and Takahashi. Similar to the case of M87, we find that the model can reproduce the logarithmic feature of the velocity profile and suggest a slowly rotating black hole magnetosphere for NGC 315. By substituting the estimated Ω F into the jet power predicted by the Blandford–Znajek mechanism, we estimate the magnetic field strength near the event horizon of the central black hole as 5 × 103 G ≲ B H ≲ 2 × 104 G. We then estimate magnetic field strengths along the jet by comparing the spectral index distribution obtained from very long baseline interferometry (VLBI) observations with a synchrotron-emitting jet model. Then we constrain the magnetic field strength at a deprojected distance z from the black hole to be in the range 0.06 G ≲ B(z) ≲ 0.9 G for 5.2 × 103 r g ≲ z ≲ 4.9 × 104 r g , where r g represents the gravitational radius. By combining the obtained field strengths at the event horizon and the downstream section of the jet, we find that the accretion flow at the jet base is consistent with a magnetically arrested disk. We discuss a comparison of the jet power and the magnetic flux anchored to the event horizon in NGC 315 and M87.

Two-component off-axis jet model for radio flares of tidal disruption events

Two-component off-axis jet model for radio flares of tidal disruption events

The fundamental plane of black hole activity for low-luminosity radio active galactic nuclei across 0 < z < 4

Context. The fundamental plane of black hole activity describes the correlation between radio luminosity (LR), X-ray luminosity (LX), and black hole mass (MBH). It reflects a connection between the accretion disc and the jet. However, the dependence of the fundamental plane on various physical properties of active galactic nuclei (AGNs) and host galaxies remains unclear, especially for low-luminosity AGNs, which is important for understanding the accretion physics in AGNs. Aims. Here, we explore the dependence of the fundamental plane on the radio loudness, Eddington-ratio (λEdd), redshift, and galaxy star formation properties (star-forming galaxies and quiescent galaxies) across 0.1 < z ≤ 4 for radio AGNs. Based on current deep and large surveys, our studies can extend to lower luminosities and higher redshifts. Methods. From the deep and large multi-wavelength surveys in the GOODS-N, GOODS-S, and COSMOS/UltraVISTA fields, we constructed a large and homogeneous radio AGN sample consisting of 208 objects with available estimates for LR and LX. Then we divided the radio AGN sample into 141 radio-quiet AGNs and 67 radio-loud AGNs according to the radio loudness defined by the ratio of LR to LX, and explored the dependence of the fundamental plane on different physical properties of the two populations, separately. Results. The ratio of LR to LX shows a bimodal distribution that is well described by two single Gaussian models. The cross point between these two Gaussian components corresponds to a radio-loudness threshold of log(LR/LX) = − 2.73. The radio-quiet AGNs have a significantly larger Eddington ratio than the radio-loud AGNs. Our radio-quiet and radio-loud AGNs show a significantly different fundamental plane, which indicates a significant dependence of the fundamental plane on the radio loudness. For both radio-quiet and radio-loud AGNs, the fundamental plane shows a significant dependence on λEdd, but no dependence on redshift. The fundamental plane shows a significant dependence on the galaxy star formation properties for radio-quiet AGNs, while for radio-loud AGNs this dependence disappears. Conclusions. The fundamental plane sheds important light on the accretion physics and X-ray emission origins of central engines. X-ray emission of radio-quiet AGNs at 0.01 < λEdd < 0.1 are produced by a combination of advection-dominated accretion flow (ADAF) and synchrotron radiation from the jet, while at 0.1 < λEdd < 1 they mainly follow the synchrotron jet model. The origins of X-ray emission of radio-loud AGNs are consistent with a combination of ADAF and the synchrotron jet model at λEdd < 0.01, agree with the synchrotron jet model at 0.01 < λEdd < 0.1, and follow a combination of the standard thin disc and a jet model at λEdd > 0.1.

Non-modal analysis of transient growth in a liquid annular jet surrounded by gas flow

Transient energy growth is a common mathematical concept in many fluid flow systems, and it has been widely investigated in recent years using non-modal analysis. Non-modal analysis can characterize the short-term energy amplification of perturbations, which is influenced by the Reynolds number, the Weber number, and the initial conditions such as the wavenumber. In gas–liquid coaxial nozzles, annular jets are often produced, and their breakup process is influenced by transient energy growth. However, research in this area has been limited so far. This paper for the first time investigates the transient energy growth of an annular liquid jet in static gas and validates it using a modified annular jet model. In the derivation process, the gas–liquid interfaces inside and outside the annular liquid film are taken into account. It has been found that there exists an optimal initial condition for a certain Reynolds number and a Weber number. The increase in the Reynolds number and ratio of inner and outer radius of the annular jet can maximize the transient growth under a specific initial wavenumber, while the increase in gas/liquid density ratio and the Weber number will minimize the transient growth. It is also found that transient energy growth is caused by the displacement of the free boundary.

Weak long‐range diffusion

AbstractWe study the spreading solutions of the nonlinear diffusion equation when the far‐field diffusivity is small. The method of strained coordinates is used to construct a uniform asymptotic correction to the similarity solution of the unperturbed problem. The equation provides a possible analogue to similar models of fluid jets and plumes.

Integrated Study of X-Ray Spectrum and Time Lags for HBL Mrk 421 within the Framework of the Multiple-zone Leptonic Model

We present the timing analysis of 10 archived XMM-Newton observations with an exposure of >40 ks of Markarian 421. Mrk 421 is the brightest high-frequency-peaked BL Lac object emitting in X-rays produced by electrons accelerated in the innermost regions of a relativistic jet pointing toward us. For each observation, we construct averaged X-ray spectra in 0.5–10 keV band, as well as 100 s binned light curves (LCs) in various subbands. During these observations, the source exhibited various intensity states differing by close to an order of magnitude in flux, with the fractional variability amplitude increasing with energy through the X-ray band. Bayesian power spectral density analysis reveals that the X-ray variability can be characterized by a colored noise, with an index ranging from ∼ −1.9 to −3.0. Moreover, both the standard cross-correlation function and cross-spectral methods indicate that the amount of time lags increases with the energy difference between two compared LCs. A time-dependent two-zone jet model is developed to extract physical information from the X-ray emission of Mrk 421. In the model, we assume that the jet emission mostly comprises a quasi-stationary component and a highly variable one. Our results show that the two-zone model can simultaneously provide a satisfactory description for both the X-ray spectra and time lags observed in different epochs, with the model parameters constrained in a fully acceptable interval. We suggest that shocks within the jets may be the primary energy dissipation process responsible for triggering the rapid variability, although magnetic reconnection cannot be excluded.

Data-driven surrogate modeling and optimization of supercritical jet into supersonic crossflow

For the design and optimization of advanced aero-engines, the prohibitively computational resources required for numerical simulations pose a significant challenge, due to the extensive exploration of design parameters across a vast design space. Surrogate modeling techniques offer a viable alternative for efficiently emulating numerical results within a notably compressed timeframe. This study introduces parametric Reduced-Order Models (ROMs) based on Convolutional AutoEncoders (CAE), Fully Connected AutoEncoders (FCAE), and Proper Orthogonal Decomposition (POD) to fast emulate spatial distributions of physical variables for a supercritical jet into a supersonic crossflow under different operating conditions. To further accelerate the decision-making process, an optimization model is developed to enhance fuel-oxidizer mixing efficiency while minimizing total pressure loss. Results indicate that CAE-based ROMs exhibit superior prediction accuracy while FCAE-based ROMs show inferior predictive accuracy but minimal uncertainty. The latter may be ascribed to the markedly greater number of hyperparameters. POD-based ROMs underperform in regions of strong nonlinear flow dynamics, coupled with higher overall prediction uncertainties. Both AE- and POD-based ROMs achieve online predictions approximately 9 orders of magnitude faster than conventional simulations. The established optimization model enables the attainment of Pareto-optimal frontiers for spatial mixing deficiencies and total pressure recovery coefficient.

Numerical simulations of downburst wind fields: A comparative analysis of stationary and moving storms using the impinging jet model

Numerical simulations of downburst wind fields: A comparative analysis of stationary and moving storms using the impinging jet model

A model for general relativistic magnetohydrodynamic spine jets

Context. We study jets using a semi-analytical model of the general relativistic magnetohydrodynamic (GRMHD) equations in the Kerr metric that describes them near the rotation axis, assuming a steady state, and axisymmetry. Aims. The goal is to model the inner spine of a relativistic jet in order to solve for the bulk acceleration and the shape of the jet and understand how these quantities depend on the enthalpy and the magnetic field. Methods. The model is constructed by expanding the rotating black hole metric and forces with respect to the polar angle about the rotation axis. This results in a system of ordinary differential equations that determine the dependence on the radial distance. The difference with previous semi-analytical models that expand the metric around the rotation axis is that the flow is governed by a polytropic equation of state. Results. The solutions in this work start from a stagnation surface very close to the event horizon and become highly relativistic, achieving large Lorentz factors at large distances.

The variable radio jet of the accreting neutron star the Rapid Burster

ABSTRACT The Rapid Burster is a unique neutron star low-mass X-ray binary system, showing both thermonuclear v-I and accretion-driven Type-II X-ray bursts. Recent studies have demonstrated how coordinated observations of X-ray and radio variability can constrain jet properties of accreting neutron stars – particularly when the X-ray variability is dominated by discrete changes. We present a simultaneous very large array, Swift, and INTErnational Gamma-Ray Astrophysics Laboratory observing campaign of the Rapid Burster to investigate whether its jet responds to Type-II bursts. We observe the radio counterpart of the X-ray binary at its faintest-detected radio luminosity, while the X-ray observations reveal prolific, fast X-ray bursting. A time-resolved analysis reveals that the radio counterpart varies significantly between observing scans, displaying a fractional variability of $38 \pm 5$ per cent. The radio faintness of the system prevents the robust identification of a causal relation between individual Type-II bursts and the evolution of the radio jet. However, based on a comparison of its low-radio luminosity with archival Rapid Burster observations and other accreting neutron stars, and on a qualitative assessment of the X-ray and radio light curves, we explore the presence of a tentative connection between bursts and jet: i.e. the Type-II bursts may weaken or strengthen the jet. The former of those two scenarios would fit with magnetorotational jet models; we discuss three lines of future research to establish this potential relation between Type-II bursts and jets more confidently.