Acoustic/shock wave heating in the gravitationally stratified partially ionized plasmas: the two-fluid effects

<p> The chromospheric heating problem is a long-standing intriguing topic of solar physics, and the acoustic wave/shock wave heating in the chromospheric plasma has been investigated in the last several decades. It has been confirmed that acoustic waves, and the shock waves induced by the steepening acoustic waves in the gravitationally stratified chromospheric plasma, are able to transport energy to the chromosphere, but the energy supplied in this way is not necessarily sufficient for heating the chromosphere. Here, we further investigate the acoustic/shock wave heating process while taking into account the two-fluid effects.</p><p> As the plasma in the chromosphere is weakly or partially ionized,  neutrals play an important role in wave propagation in this region. Therefore,  a two-fluid computational model treating neutrals and charged particles (electrons and ions) as two separate fluids is used for modelling the acoustic/shock wave propagation in idealised partially ionized plasmas, while taking into account the ion-neutral collisions, ionization and recombination. We have thus investigated  the collisional and reactive interactions between separated ions and neutrals, as well as the resulting effects in the acoustic/shock wave propagation and damping. In the numerical simulations, both the initial hydrostatic equilibrium and chemical equilibrium are taken into account to provide different density profiles for comparison.</p><p>We have found that the shock heating in the partially ionized plasmas strongly depends on the ionization fraction. In particular, the relatively smaller ionization fraction resulting from the initial chemical equilibrium significantly enhances the shock wave heating, which dominates the overall heating effect according to an approximated quantitative comparison. Moreover, the decoupling between ions and neutrals is also enhanced while considering ionization and recombination, resulting in stronger collisional heating.</p>

Revealing the cosmic reionization history with fast radio bursts in the era of Square Kilometre Array

ABSTRACT Revealing the cosmic reionization history is at the frontier of extragalactic astronomy. The power spectrum of the cosmic microwave background (CMB) polarization can be used to constrain the reionization history. Here, we propose a CMB-independent method using fast radio bursts (FRBs) to directly measure the ionization fraction of the intergalactic medium (IGM) as a function of redshift. FRBs are new astronomical transients with millisecond time-scales. Their dispersion measure (DMIGM) is an indicator of the amount of ionized material in the IGM. Since the differential of DMIGM against redshift is proportional to the ionization fraction, our method allows us to directly measure the reionization history without any assumption on its functional shape. As a proof of concept, we constructed mock non-repeating FRB sources to be detected with the Square Kilometre Array, assuming three different reionization histories with the same optical depth of Thomson scattering. We considered three cases of redshift measurements: (A) spectroscopic redshift for all mock data, (B) spectroscopic redshift for 10 per cent of mock data, and (C) redshift estimated from an empirical relation of FRBs between their time-integrated luminosity and rest-frame intrinsic duration. In all cases, the reionization histories are consistently reconstructed from the mock FRB data using our method. Our results demonstrate the capability of future FRBs in constraining the reionization history.

Tracers of the ionization fraction in dense and translucent gas

Context. The ionization fraction in the neutral interstellar medium (ISM) plays a key role in the physics and chemistry of the ISM, from controlling the coupling of the gas to the magnetic field to allowing fast ion-neutral reactions that drive interstellar chemistry. Most estimations of the ionization fraction have relied on deuterated species such as DCO+, whose detection is limited to dense cores representing an extremely small fraction of the volume of the giant molecular clouds that they are part of. As large field-of-view hyperspectral maps become available, new tracers may be found. The growth of observational datasets is paralleled by the growth of massive modeling datasets and new methods need to be devised to exploit the wealth of information they contain. Aims. We search for the best observable tracers of the ionization fraction based on a grid of astrochemical models, with the broader aim of finding a general automated method applicable to searching for tracers of any unobservable quantity based on grids of models. Methods. We built grids of models that randomly sample a large range of physical conditions (unobservable quantities such as gas density, temperature, elemental abundances, etc.) and computed the corresponding observables (line intensities, column densities) and the ionization fraction. We estimated the predictive power of each potential tracer by training a random forest model to predict the ionization fraction from that tracer, based on these model grids. Results. In both translucent medium and cold dense medium conditions, we found several observable tracers with very good predictive power for the ionization fraction. Many tracers in cold dense medium conditions are found to be better and more widely applicable than the traditional DCO+/HCO+ ratio. We also provide simpler analytical fits for estimating the ionization fraction from the best tracers, and for estimating the associated uncertainties. We discuss the limitations of the present study and select a few recommended tracers in both types of conditions. Conclusions. The method presented here is very general and can be applied to the measurement of any other quantity of interest (cosmic ray flux, elemental abundances, etc.) from any type of model (PDR models, time-dependent chemical models, etc.).

Accretion bursts in magnetized gas-dust protoplanetary disks

Aims. Accretion bursts triggered by the magnetorotational instability (MRI) in the innermost disk regions were studied for protoplanetary gas-dust disks that formed from prestellar cores of a various mass Mcore and mass-to-magnetic flux ratio λ. Methods. Numerical magnetohydrodynamics simulations in the thin-disk limit were employed to study the long-term (~1.0 Myr) evolution of protoplanetary disks with an adaptive turbulent α-parameter, which explicitly depends on the strength of the magnetic field and ionization fraction in the disk. The numerical models also feature the co-evolution of gas and dust, including the back-reaction of dust on gas and dust growth. Results. A dead zone with a low ionization fraction of x≲10−13 and temperature on the order of several hundred Kelvin forms in the inner disk soon after its formation, extending from several to several tens of astronomical units depending on the model. The dead zone features pronounced dust rings that are formed due to the concentration of grown dust particles in the local pressure maxima. Thermal ionization of alkaline metals in the dead zone trigger the MRI and associated accretion burst, which is characterized by a sharp rise, small-scale variability in the active phase, and fast decline once the inner MRI-active region is depleted of matter. The burst occurrence frequency is highest in the initial stages of disk formation and is driven by gravitational instability (GI), but it declines with diminishing disk mass-loading from the infalling envelope. There is a causal link between the initial burst activity and the strength of GI in the disk fueled by mass infall from the envelope. We find that the MRI-driven burst phenomenon occurs for λ = 2–10, but diminishes in models with Mcore ≲ M⊙, suggesting a lower limit on the stellar mass for which the MRI-triggered burst can occur. Conclusions. The MRI-triggered bursts occur for a wide range of mass-to-magnetic flux ratios and initial cloud core masses. The burst occurrence frequency is highest in the initial disk formation stage and reduces as the disk evolves from a gravitationally unstable to a viscous-dominated state. The MRI-triggered bursts are intrinsically connected with the dust rings in the inner disk regions, and both can be a manifestation of the same phenomenon, that is to say the formation of a dead zone.

Implications of baryon–dark matter interaction on IGM temperature and tSZ effect with magnetic field

ABSTRACT We show that the combined effect of cosmic magnetic field and a possible non-standard interaction between baryons and dark matter (DM) has interesting consequences on the thermal Sunyaev−Zel’dovich (tSZ) effect depending on the temperature and the ionization state of the intergalactic medium. The drag force between the baryons and DM due to the relative velocity between them, and their temperature difference results in heat transfer between these two species. At the same time, the ambipolar diffusion and the decaying magnetic turbulence tends to heat up the baryons. This interplay of these two processes give rise to different evolution histories of the thermal and ionization state of the universe and hence influences the cosmic microwave background (CMB) spectrum at small scales through the tSZ effect. In this work, we have computed the evolution of the temperature, ionization fraction, and the y-parameter of the CMB for different strengths of the magnetic field and the interaction cross-section. We note that the y-parameter can be significantly enhanced with the inclusion of magnetic field and baryon–DM interaction as compared to the case when these are absent. The enhancement depends on the strength of the magnetic field.

Analysing the Epoch of Reionization with three-point correlation functions and machine learning techniques

ABSTRACT Three-point and high-order clustering statistics of the high-redshift 21 cm signal contain valuable information about the Epoch of Reionization (EoR). We present 3PCF-fast, an optimized code for estimating the three-point correlation function (3PCF) of 3D pixelized data such as the outputs from numerical and seminumerical simulations. After testing 3PCF-fast on data with known analytical 3PCF, we use machine learning techniques to recover the mean bubble size and global ionization fraction from correlations in the outputs of the publicly available 21cmfast code. We assume that foregrounds have been perfectly removed and negligible instrumental noise. Using ionization fraction data, our best multilayer perceptron (MLP) model recovers the mean bubble size with a median prediction error of around $10 {{\ \rm per\ cent}}$, or from the 21 cm differential brightness temperature with median prediction error of around $14 {{\ \rm per\ cent}}$. A further two MLP models recover the global ionization fraction with median prediction errors of around $4 {{\ \rm per\ cent}}$ (using ionization fraction data) or around $16 {{\ \rm per\ cent}}$ (using brightness temperature). Our results indicate that clustering in both the ionization fraction field and the brightness temperature field encode useful information about the progress of the EoR in a complementary way to other summary statistics. Using clustering would be particularly useful in regimes where high signal-to-noise ratio prevents direct measurement of bubble size statistics. We compare the quality of MLP models using the power spectrum, and find that using the 3PCF outperforms the power spectrum at predicting both global ionization fraction and mean bubble size.

Probing the size and charge of polycyclic aromatic hydrocarbons

ABSTRACT We present a new method to accurately describe the ionization fraction and the size distribution of polycyclic aromatic hydrocarbons (PAHs) within astrophysical sources. To this purpose, we have computed the mid-infrared emission spectra of 308 PAH molecules of varying sizes, symmetries, and compactness, generated in a range of radiation fields. We show that the intensity ratio of the solo CH out-of-plane bending mode in PAH cations and anions (referred to as the ‘11.0’ μm band, falling in the 11.0–11.3 μm region for cations and anions) to their 3.3 μm emission scales with PAH size, similarly to the scaling of the 11.2/3.3 ratio with the number of carbon atoms (NC) for neutral molecules. Among the different PAH emission bands, it is the 3.3 μm band intensity that has the strongest correlation with NC, and drives the reported PAH intensity ratio correlations with NC for both neutral and ionized PAHs. The 6.2/7.7 intensity ratio, previously adopted to track PAH size, shows no evident scaling with NC in our large sample. We define a new diagnostic grid space to probe PAH charge and size, using the (11.2 + 11.0)/7.7 and (11.2 + 11.0)/3.3 PAH intensity ratios, respectively. We demonstrate the application of the (11.2 + 11.0)/7.7–(11.2 + 11.0)/3.3 diagnostic grid for galaxies M82 and NGC 253, for the planetary nebula NGC 7027, and the reflection nebulae NGC 2023 and NGC 7023. Finally, we provide quantitative relations for PAH size determination depending on the ionization fraction of the PAHs and the radiation field they are exposed to.

Photophoresis in the circumjovian disk and its impact on the orbital configuration of the Galilean satellites

Jupiter has four large regular satellites called the Galilean satellites: Io, Europa, Ganymede, and Callisto. The inner three of the Galilean satellites orbit in a 4:2:1 mean motion resonance; therefore their orbital configuration may originate from the stopping of the migration of Io near the bump in the surface density distribution and following resonant trapping of Europa and Ganymede. The formation mechanism of the bump near the orbit of the innermost satellite, Io, is not yet understood, however. Here, we show that photophoresis in the circumjovian disk could be the cause of the bump using analytic calculations of steady-state accretion disks. We propose that photophoresis in the circumjovian disk could stop the inward migration of dust particles near the orbit of Io. The resulting dust-depleted inner region would have a higher ionization fraction, and thus admit increased magnetorotational-instability-driven accretion stress in comparison to the outer region. The increase of the accretion stress at the photophoretic dust barrier would form a bump in the surface density distribution, halting the migration of Io.