Radiative Transfer Calculations Research Articles

Multiline study of the radial extent of SiO, CS, and SiS in asymptotic giant branch envelopes

Circumstellar envelopes around asymptotic giant branch (AGB) stars contain a rich diversity of molecules, whose spatial distribution is regulated by different chemical and physical processes. In the outer circumstellar layers, all molecules are efficiently destroyed due to interactions with interstellar ultraviolet photons. Here we aim to carry out a coherent and uniform characterization of the radial extent of three molecules (SiO, CS, and SiS) in envelopes around O- and C-rich AGB stars, and to study their dependence on mass-loss rate. To this end, we observed a reduced sample of seven M-type and seven C-type AGB envelopes in multiple lines of SiO, CS, and SiS with the Yebes\,40m and IRAM\,30m telescopes. The selected sources cover a wide range of mass-loss rates, from sim $ $M_ odot $ yr$^ $ to a few times 10$^ $ $M_ odot $ yr$^ $, and the observed lines cover a wide range of upper-level energies, from 2 K to 130 K. We carried out excitation and radiative transfer calculations over a wide parameter space in order to characterize the abundance and radial extent of each molecule. A $ analysis indicates that the abundance is usually well constrained while the radial extent is in some cases more difficult to constrain. Our results indicate that the radial extent of the molecules considered here increases with increasing envelope density, in agreement with previous observational findings. At high envelope densities of $ M /v_ infty $ $M_ odot $ yr$^ $ km$^ $ s, SiO, CS, and SiS show a similar radial extent, while at low envelope densities of $ M /v_ infty $ $M_ odot $ yr$^ $ km$^ $ s, differences in radial extent appear among the three molecules, in agreement with theoretical expectations based on destruction due to photodissociation. At low envelope densities, we find a sequence of increasing radial extent, SiS\,rightarrow \,CS\,rightarrow \,SiO. We also find a tentative dependence of the radial extent on the chemical type (O- or C-rich) of the star for SiO and CS. Interferometric observations and further investigation of the photodissociation of SiO, CS, and SiS should provide clarification of the situation in regards to the relative photodissociation radii of SiO, CS, and SiS in AGB envelopes and their dependence on envelope density and C/O ratio.

Correction of cloud optical thickness retrievals from nadir reflectances in the presence of 3D radiative effects. Part I: concept and tests on 3D RT simulations

3D effects cause substantial underestimation of cloud optical thickness (COT) in airborne and satellite retrievals based on 1D radiative transfer computations (such as in the case of widely used bispectral technique). For a single-layer isolated cloud we propose a simple linear correction of the retrieved COT with the renormalization factor dependent on the cloud’s aspect ratio (the ratio between vertical and horizontal dimensions of the cloud). This is an empirical assumption which we successfully test using synthetic 3D RT data. We introduce a heuristic “block model” of 3D radiative effects and show that the functional form of the renormalization factor is consistent with the process of radiation escape from cloud sides in an essentially 3D geometry. We also extend the block model to the case of single-layer broken cloud field with radiative interaction between the neighboring clouds. In this case the renormalization factor depends also on the distance between clouds.

Dual Role of Accretion Disk Winds as X-Ray Obscurers and UV Line Absorbers in AGN

X-ray obscuration of active galactic nuclei (AGNs) is considered in the context of ionized winds of stratified structure launched from accretion disks. We argue that a Compton-thick layer of a large-scale disk wind can obscure continuum X-rays and also lead to broad UV absorption, such as in the blue wing of C iv; the former originates from the inner wind and the latter from the outer wind, as a dual role. Motivated by a number of lines of observational evidence showing strong AGN obscuration phenomena in Seyfert 1 AGNs such as NGC 5548, we demonstrate in this work, by utilizing a physically motivated wind model coupled to post-process radiative transfer calculations, that an extended disk wind under certain physical conditions (e.g., morphology and density) could naturally cause a sufficient obscuration qualitatively consistent with UV/X-ray observations. Predicted UV/X-ray correlation is also presented as a consequence of variable spatial size of the wind in this scenario.

Satellite-Based Estimation of the Role of Cloud-Radiative Interaction in Accelerating Tropical Cyclone Development

Abstract Recent modeling studies have suggested a potentially important role of cloud-radiative interactions in accelerating tropical cyclone (TC) development, but there has been only limited investigation of this in observations. Here, we investigate this by performing radiative transfer calculations based on cloud property retrievals from the CloudSat Tropical Cyclone (CSTC) dataset. We examine the radius–height structure of radiative heating anomalies, compute the resulting radiatively driven circulations, and use the moist static energy variance budget to compute radiative feedbacks. We find that inner-core midlevel ice water content and anomalous specific humidity increase with TC intensification rate, resulting in enhanced inner-core deep-layer longwave warming anomalies and shortwave cooling anomalies in rapidly intensifying TCs. This leads to a stronger radiatively driven deep in-up-and-out overturning circulation and inner-core radiative feedback in rapidly intensifying TCs. The longwave-driven circulation provides radially inward momentum fluxes and upward moisture fluxes, which benefit TC development, while the shortwave-driven circulation suppresses TC development. The longwave anomalies, which dominate the inner-core positive radiative feedback, are mainly generated from cloud-radiative interactions, with ice particles dominating the deep-layer circulation and liquid droplets and water vapor contributing to the shallow circulation. Moreover, the variability in ice water content, as opposed to the variability in liquid water content and the effective radii of ice particles and liquid droplets, dominates the uncertainty in TC-radiative interaction. These results provide observational evidence for the importance of cloud-radiative interactions in TC development and suggest that the amount and spatial structure of ice water content are critical for determining the strength of this interaction. Significance Statement The limited investigation of tropical cyclone (TC)-radiative interaction in observations impedes our understanding of TC development. This study aims to quantitatively show the spatial variation in radiation in TCs and their effect on TC development by using a set of satellite-based observations. We relate TC-radiative interaction to TC intensification and emphasize the inner-core features. Moreover, we quantitatively demonstrate the relative contribution from clouds, liquid droplets, ice particles, and water vapor to TC-radiative interaction as well as the source of the variation in radiative properties. These results provide an additional observational foundation for the importance of cloud-radiative interactions in TC development and support a quantitative validation for numerical modeling.

Dynamics and emission properties of flux ropes from two-temperature GRMHD simulations with multiple magnetic loops

Flux ropes erupting in the vicinity of a black hole are thought to be a potential model for the flares observed in Sagittarius\,A$^*$. In this study, we examine the radiative properties of flux ropes that emerged in the vicinity of a black hole. We performed three-dimensional two-temperature general relativistic magnetohydrodynamic (GRMHD) simulations of magnetized accretion flows with alternating multiple magnetic loops and general relativistic radiation transfer (GRRT) calculations. In the GRMHD simulations, we implemented two different sizes of initial magnetic loops. In the small loop case, magnetic dissipation leads to a weaker excitement of magneto-rotational instability inside the torus, which generates a lower accretion rate compared to the large loop case. However, the small loop case generates more flux ropes due to frequent reconnection by magnetic loops with different polarities. By calculating the thermal synchrotron emission, we found that the variability of light curves and the emitting region are tightly related. At $230\ GHz$ and higher, the emission from the flux ropes is relatively stronger compared to the background, which is responsible for the filamentary structure in the images. At lower frequencies (e.g. $43\ GHz$), emission comes from more extended regions, which have a less filamentary structure in the image. Our study shows that self-consistent electron temperature models are essential for the calculation of thermal synchrotron radiation and the morphology of the GRRT images. Flux ropes contribute considerable emission at frequencies $ GHz$.

H2CO and CS in diffuse clouds: Excitation and abundance

Context. Diffuse interstellar clouds present an active chemistry despite their relatively low density and the ubiquitous presence of far-UV radiation. Aims. To provide constraints on the chemical processes responsible for the observed columns of organic species, we used the NOEMA interferometer to observe the sight line toward NRAO150 (B0355+508) in the 2 mm spectral window. Methods. We targeted the low excitation lines of ortho H2 CO (21,1–11,0) and para H2 CO (20,2–10,1) as well as the nearby transitions of CS (3–2) and c-C3H2 (31,2–22,1), (41,4–30,3), and (22,0–11,1). We combined these data with previous observations of the same sight line to determine the excitation conditions, column densities, and abundances relative to H2 in the different velocity components. We performed non-LTE radiative transfer calculations including collision cross sections with ortho and para H2 and with electrons. New collision cross sections with electrons were computed for ortho and para formaldehyde. Results. All targeted lines were detected. The c-C3H2 line profiles are very similar to those of HCO+ and CCH, while the CS absorption features are narrower and mostly concentrated in two main velocity components at VLSR = −17.2 and −10.4 km s−1. H2CO absorption lines present an intermediate pattern with absorption in all velocity components but larger opacities in the two main velocity components. The ortho-to-para ratios of H2CO and c-C3H2 are consistent with the statistical value of three. While the excitation temperature of all c-C3H2 velocity components is consistent with the Cosmic Microwave Background (CMB), the two strong components detected in CS show a clear excess over the CMB indicating that CS resides at higher densities than other species along this particular sightline, n(H2) ~ 2500 cm−3 while n(H2) < 500 cm−3 for the other velocity components. We detected faint absorption from o-H213CO and C34S allowing us to derive isotopic ratios: o-H2CO/o-H213CO = 61 ± 12 and C32S/C34S = 24 ± 6. The excitation of the 11,0−11,1 line of formaldehyde at 4.8 GHz is sensitive to the electron fraction and its excitation temperature is predicted to be lower than the CMB at low and moderate electron fractions (x(e) < 6 × 10−5), and to rise above the CMB at high electron fractions (x(e) > 10−4).

Global climate modeling of the Jupiter troposphere and effect of dry and moist convection on jets

The atmosphere of Jupiter is characterized by banded jets, including an equatorial super-rotating jet, by an intense moist convective activity, and by perturbations exerted by vortices, waves, and turbulence. Even after space exploration missions to Jupiter and detailed numerical modeling of Jupiter, questions remain about the mechanisms underlying the banded jets and the role played by dry and moist convection in maintaining these jets. We report three-dimensional simulations of the Jupiter weather layer using a global climate model (GCM) called Jupiter-DYNAMICO, which couples hydrodynamical integrations on an icosahedral grid with detailed radiative transfer computations. We added a thermal plume model for Jupiter that emulates the effect of mixing of heat, momentum, and tracers by dry and moist convective plumes that are left unresolved in the GCM mesh spacing with a physics-based approach. Our Jupiter-DYNAMICO global climate simulations show that the large-scale Jovian flow, in particular the jet structure, could be highly sensitive to the water abundance in the troposphere and that an abundance threshold exists at which equatorial super-rotation develops. In contrast to our dry (or weakly moist) simulations, simulations that include the observed amount of tropospheric water exhibit a clear-cut super-rotating eastward jet at the equator and a dozen eastward mid-latitude jets that do not migrate poleward. The magnitudes agree with the observations. The convective activity simulated by our thermal plume model is weaker in the equatorial regions than in mid to high latitudes, as indicated by lightning observations. Regardless of whether they are dry or moist, our simulations exhibit the observed inverse energy cascade from small (eddies) to large scales (jets) in a zonostrophic regime.

Radiative transfer of 21-cm line through ionised cavities in an expanding universe

Abstract The optical depth parameterisation is typically used to study the 21-cm signals associated with the properties of the neutral hydrogen (H i) gas and the ionisation morphology during the Epoch of Reionisation (EoR), without solving the radiative transfer equation. To assess the uncertainties resulting from this simplification, we conduct explicit radiative transfer calculations using the cosmological 21-cm radiative transfer (C21LRT) code and examine the imprints of ionisation structures on the 21-cm spectrum. We consider a globally averaged reionisation history and implement fully ionised cavities (H ii bubbles) of diameters d ranging from 0.01 Mpc to 10 Mpc at epochs within the emission and the absorption regimes of the 21-cm global signal. The single-ray C21LRT calculations show that the shape of the imprinted spectral features are primarily determined by d and the 21-cm line profile, which is parametrised by the turbulent velocity of the H i gas. It reveals the spectral features tied to the transition from ionised to neutral regions that calculations based on the optical depth parametrisation were unable to capture. We also present analytical approximations of the calculated spectral features of the H ii bubbles. The multiple-ray calculations show that the apparent shape of a H ii bubble (of d = 5 Mpc at z = 8), because of the finite speed of light, differs depending on whether the bubble’s ionisation front is stationary or expanding. Our study shows the necessity of properly accounting for the effects of line-continuum interaction, line broadening and cosmological expansion to correctly predict the EoR 21-cm signals.

A systematic analysis of energy-saving potential of spectrum-selective radiative cooling materials in buildings

Passive cooling of buildings through radiative cooling (RC) is considered an effective method for reducing the energy consumption of air conditioning. In recent years, spectrum-selective RC materials (RCMs) have received extensive attention owing to their low solar absorptivity and high thermal emissivity in the atmospheric window. The energy-saving performance of RCMs is related to their own properties and local meteorological conditions, building type, building form, and energy use patterns. Although several studies have explored the performance of RCMs, most existing studies have considered a single building type and meteorological condition, and simple building forms that differ significantly from real buildings. Hence, the application analyses of RCMs to buildings are not sufficiently comprehensive. This study focused on the energy-saving potential of RCMs for building applications, and presented a systematic framework to analyze the energy-saving potential of RCMs under different conditions. DeST, which is embedded with a long-wave radiative heat transfer calculation module based on the equivalent sky temperature, was adopted as the simulation software. It can accurately calculate the radiative heat transfer in different wavelength bands, thus making the simulation results more accurate. Based on the calculation results of large-scale cases considering four major factors, namely climate, building type, construction year, and installation location of RCMs, the energy-saving performance of RCMs in different cases was analyzed, and the return on investment is performed, and more comprehensive suggestions were provided for applying RCMs in future buildings.

Towards an observationally motivated AGN dusty torus model – I. Dust chemical composition from the modelling of Spitzer spectra

ABSTRACT Spectral energy distribution (SED) fitting is one of the most commonly used techniques to study the dust properties in active galactic nuclei (AGN). Works implementing this technique commonly use radiative transfer models that assume a variety of dust properties. Despite the key role of this aspect, limited effort has been put forward to explore the chemical composition, the role of different optical properties, and the grain size distribution of dust, all of which can have a substantial impact on the theoretical radiative transfer calculations. In this work, we explore the role of the dust chemical composition in the AGN dusty torus through SED fitting to Spitzer/IRS spectra of a sample of 49 nearby AGN with silicate features in emission. We implement a mineralogy model including the popular astronomical silicates and a set of oxides and amorphous silicates with different grain sizes. Wefind that best fits use principally porous alumina, periclase, and olivine. In terms of mass fractions, $\sim 99~{{\ \rm per\ cent}}$ of the dust is composed of dust grains of size $\rm {0.1 \, \mu m}$, with a $\lt 1~{{\ \rm per\ cent}}$ contribution from $\rm {3 \, \mu m}$ grains. Moreover, the astronomical silicates have a very low occurrence in the best fits, suggesting that they are not the most suited dust species to reproduce the silicate features in our sample.

CAI formation in the early Solar System

Calcium-aluminium-rich inclusions (CAIs) are the oldest dated solid materials in the Solar System, and are found as light-coloured crystalline ingredients in carbonaceous chondrite meteorites. Their formation time is commonly associated with age zero of the Solar System. Nevertheless, the physical and chemical processes that once led to the formation of these submillimetre- to centimetre-sized mineral particles in the early solar nebula are still a matter of debate. In this paper, we propose a pathway to form such inclusions during the earliest phases of disc evolution. We combine 1D viscous disc evolutionary models with 2D radiative transfer, equilibrium condensation, and new dust opacity calculations. We show that the viscous heating associated with the high accretion rates in the earliest evolutionary phases causes the midplane inside of about 0.5\,au to heat up to limiting temperatures of about 1500-1700\,K, but no further. These high temperatures force all refractory material components of the inherited interstellar dust grains to sublimate --- except for a few Al-Ca-Ti oxides, such as Al$_2$O$_3$, Ca$_2$Al$_2$SiO$_7$, and CaTiO$_3$. This is a recurring and very stable result in all our simulations, because these minerals form a natural thermostat. Once the Mg-Fe silicates are gone, the dust becomes more transparent and the heat is more efficiently transported to the disc surface, which prevents further warming. This thermostat mechanism keeps these minerals above their annealing temperature for hundreds of thousands of years, allowing them to form large pure crystalline particles. These particles are dragged out by the viscously spreading disc, and once they reach a distance of about 0.5\,au, the silicates recondense on the surface of the Ca-Al-rich particles, adding an amorphous silicate matrix. We estimate that this mechanism of CAI production works during the first 50\,000 years of disc evolution. These particles then continue to move outward and populate the entire disc up to radii of about 50\,au, before the accretion rate eventually subsides, the disc cools, and the particles start to drift inwards.

Determination of interparticle radiative heat transfer and radiation absorption for ellipsoidal particle beds using Monte Carlo ray tracing

The calculation of radiative heat transfer between particles in particle-based solar systems presents significant challenges in terms of computational cost and accuracy. This study addresses these challenges by employing the discrete element method to generate random particle packings and the Monte Carlo ray tracing method to investigate the influence of particle shape and volume fraction (ϕ) on interparticle radiation transfer and the absorption of incident radiation. The results reveal that, at low volume fractions, the shape and orientation of ellipsoidal particles have a pronounced impact on direct radiative transfer. Oblate particles exhibit the highest average radiation distribution factor (RDF), and the dispersion of RDF data is more prominent for oblate and elongated prolate particles compared to spheres. As ϕ increases to 0.35, the contribution of indirect transfer from neighboring particles becomes significant for oblate particles, resulting in a reduction in the dispersion of RDF data. Furthermore, when ϕ reaches 0.55, the decisive factor in determining RDF becomes the shielding effect of particles positioned between the emitting and receiving particles. The average RDF tends to converge among particles of different shapes, although oblate and elongated prolate particles still exhibit considerably larger RDF dispersion compared to spherical and near-spherical particles. Additionally, the effective absorptivity of the particle bed for incident radiation decreases as the aspect ratio of particles increases at all volume fractions. Moreover, the transmitted fraction demonstrates an exponential relationship with the depth, and beds composed of spherical particles exhibit better radiation transmittance compared to beds containing non-spherical particles.

Enhanced Blandford Znajek jet in loop quantum black hole

The Blandford-Znajek (BZ) process powers energetic jets by extracting the rotating energy of a Kerr black hole. It is important to understand this process in non-Kerr black hole spacetimes. In this study, we conduct two-dimensional and three-dimensional two-temperature General Relativistic Magnetohydrodynamic (GRMHD) simulations of magnetized accretion flows onto a rotating Loop-Quantum black hole (LQBH). Our investigation focuses on the accretion flow structure and jet launching dynamics from our simulations.We observe that the loop quantum effects increase the black hole angular frequency for spinning black holes.This phenomenon intensifies the frame-dragging effect, leading to an amplification of the toroidal magnetic field within the funnel region and enhancement of the launching jet power.It is possible to fit the jet power following a similar fitting formula of the black hole angular frequency as seen in the Kerr black hole.Based on the General Relativistic Radiation Transfer (GRRT) calculation, we find that the jet image from LQBH has a wider opening angle and an extended structure than the Kerr BH.

TeV pion bumps in the gamma-ray spectra of flaring blazars

Context. Very high-energy (VHE, E > 100 GeV) observations of the blazar Mrk 501 with MAGIC in 2014 provided evidence for an unusual narrow spectral feature at about 3 TeV during an extreme X-ray flaring activity. The one-zone synchrotron-self Compton scenario, widely used in blazar broadband spectral modeling, fails to explain the narrow TeV component. Aims. Motivated by this rare observation, we propose an alternative model for the production of narrow features in the VHE spectra of flaring blazars. These spectral features may result from the decay of neutral pions (π0 bumps) that are in turn produced via interactions of protons (of tens of TeV energy) with energetic photons, whose density increases during hard X-ray flares. Methods. We explored the conditions needed for the emergence of narrow π0 bumps in VHE blazar spectra during X-ray flares reaching synchrotron energies ∼100 keV using time-dependent radiative transfer calculations. We focused on high-synchrotron peaked (HSP) blazars, which comprise the majority of VHE-detected extragalactic sources. Results. We find that synchrotron-dominated flares with peak energies ≳100 keV can be ideal periods for the search of π0 bumps in the VHE spectra of HSP blazars. The flaring region is optically thin to photopion production, its energy content is dominated by the relativistic proton population, and the inferred jet power is highly super-Eddington. Application of the model to the spectral energy distribution of Mrk 501 on MJD 56857.98 shows that the VHE spectrum of the flare is described well by the sum of a synchrotron self-Compton (SSC) component and a distinct π0 bump centered at 3 TeV. Spectral fitting of simulated SSC+π0 spectra for the Cherenkov Telescope Array (CTA) show that a π0 bump could be detected at a 5σ significance level with a 30-min exposure. Conclusions. A harder VHE γ-ray spectrum than the usual SSC prediction or, more occasionally, a distinct narrow bump at VHE energies during hard X-ray flares, can be suggestive of a relativistic hadronic component in blazar jets that otherwise would remain hidden. The production of narrow features or spectral hardenings due to π0 decay in the VHE spectra of blazars is testable with the advent of CTA.

Revisiting the rovibrational (de-)excitation of molecular hydrogen by helium

Context. The collisional (de-)excitation of H2 by He plays an important role in the thermal balance and chemistry of various astro-physical environments, making accurate rate coefficients essential for interpreting observations of the interstellar medium. Aims. Our goal is to utilize a state-of-the-art potential energy surface (PES) to provide comprehensive state-to-state rate coefficients for He-induced transitions among rovibrational levels of H2. Methods. We performed quantum scattering calculations for the H2-He system. Thus, we were able to provide state-to-state rate coefficients for 1059 transitions between rovibrational levels of H2, with internal energies up to ≃15 000 cm−1, for temperatures ranging from 20 to 8000 K. Results. Our results demonstrate a good agreement with previous calculations for pure rotational transitions between low-lying rotational levels. However, we do find significant discrepancies for rovibrational processes involving highly-excited rotational and vibrational states. We attribute these differences to two key factors: (1) the broader range of intramolecular distances covered by ab initio points and (2) the superior accuracy of the PES, resulting from the utilization of the state-of-the-art quantum chemistry methods, compared to previous lower-level calculations. Conclusions. Radiative transfer calculations performed with the new collisional data indicate that the population of rotational levels in excited vibrational states experiences significant modifications, highlighting the critical need for this updated dataset in models of high-temperature astrophysical environments.

The effects of changes in HITRAN and the water vapor continuum model on infrared radiative transfer calculations and remote sensing applications

This study examines the implications of recent updates to the High Resolution Transmission (HITRAN) database and the MT_CKD water vapor continuum model on infrared radiative transfer calculations and satellite remote sensing applications. Specifically, it assesses the impact of the latest versions of the HITRAN database (i.e. HITRAN2016 and HITRAN2020) and the MT_CKD WV continuum model (i.e. versions 3.2 and 4.1.1) on clear-sky infrared flux calculations and the simulation of satellite infrared channels, specifically those of the Geostationary Operational Environmental Satellite (GOES) Advanced Baseline Imager (ABI). The results indicate that updates to the MT_CKD model tend to reduce atmospheric opacity, while updates to the HITRAN database have the opposite effect. The most significant differences noted are attributed to updates in the MT_CKD, leading to an increase in upward flux of 0.13 Wm−2 for the TRO and a decrease in downward flux of 0.28 Wm−2 for the MLS profile. However, these impacts are not uniform across the spectrum, atmospheric height, or among different atmospheric profiles. Notably, changes in the WV continuum model have significant effects in the far-infrared region. The updates in HITRAN exhibit pronounced differences associated with changes in the spectroscopic parameters of H2O, CO2 and O3. The results of GOES’ ABI channels indicate that changes in HITRAN lead to an average brightness temperature decrease of 0.3 K for channel 12, which exceeds the channel’s sensitivity. This research underscores the importance of continuous updates and evaluations of spectroscopic databases and radiative transfer models to improve the accuracy of atmospheric remote sensing data.

Impacts of nonthermal emission on the images of a black hole shadow and extended jets in two-temperature GRMHD simulations

The recent 230\,GHz observations from the Event Horizon Telescope (EHT) collaboration are able to image the innermost structure of the M87 galaxy showing the shadow of the black hole, a photon ring, and a ring-like structure that agrees with thermal synchrotron emission from the accretion disc. However, at lower frequencies, M87 is characterized by a large-scale jet with clear signatures of nonthermal emission. It is necessary to explore the impacts of nonthermal emission on black hole shadow images and extended jets, especially at lower frequencies. In this study, we aim to compare models with different electron heating prescriptions to one another and to investigate how these prescriptions and nonthermal electron distributions may affect black hole shadow images and the broadband spectrum energy distribution (SED) function. We performed general relativistic radiative transfer (GRRT) calculations in various two-temperature general relativistic magnetohydrodynamic (GRMHD) models utilizing different black hole spins and different electron heating prescriptions coupled with different electron distribution functions (eDFs). Through a comparison with GRRT images and SEDs, we found that when considering a variable kappa eDF, the parameterized prescription of the $R- electron temperature model with h $ = 1 is similar to the model with electron heating in the morphology of images, and the SEDs at a high frequency. This is consistent with previous studies using the thermal eDF. However, the nuance between them could be differentiated through the diffuse extended structure seen in GRRT images, especially at a lower frequency, and the behavior of SEDs at low frequency. The emission from the nearside jet region is enhanced in the case of electron heating provided by magnetic reconnection and it will increase if the contribution from the regions with stronger magnetization is included or if the magnetic energy contribution to kappa eDF mainly in the magnetized regions is considered. Compared with the thermal eDF, the peaks of the SEDs shift to a lower frequency when we consider nonthermal eDF.

Supernova remnants of red supergiants: From barrels to loops

Core-collapse (CC) supernova remnants (SNRs) are the nebular leftovers of defunct massive stars that died during a supernova explosion, mostly while undergoing the red supergiant phase of their evolution. The morphology and emission properties of those remnants are a function of the distribution of circumstellar material at the moment of the supernova, as well as the intrinsic properties of the explosion and those of the ambient medium. By means of 2.5-dimensional (2.5D) numerical magneto-hydrodynamic (MHD) simulations, we modelled the long-term evolution of SNRs generated by runaway rotating massive stars moving into a magnetised interstellar medium (ISM). Radiative transfer calculations reveal that the projected non-thermal emission of SNRs decreases over time, namely: older remnants are fainter than younger ones. Older ($80\ kyr$) SNRs whose progenitors were moving with a space velocity corresponding to a Mach number of $M=1$ ($v_ km\ $) in the Galactic plane of the interstellar medium ($n_ ISM $) are brighter in synchrotron than when moving with a Mach number of $M=2$ ($v_ km\ $). We show that runaway red supergiant progenitors first induce an asymmetric non-thermal $1.4\ GHz$ barrel-like synchrotron SNRs (at the age of about $8\ kyr$), before further evolving to adopt a Cygnus-loop-like shape (at about $80\ kyr$). It is conjectured that a significative fraction of SNRs are currently in this bilateral-to-Cygnus loop evolutionary sequence. Therefore, this population should be taken into account with repect to interpreting the data as part of the forthcoming Cherenkov Telescope Array (CTA) observatory.

Uncertainties in cloud-radiative heating within an idealized extratropical cyclone

Abstract. Cloud-radiative heating (CRH) within the atmosphere affects the dynamics and predictability of extratropical cyclones. However, CRH is uncertain in numerical weather prediction and climate models, and this could affect model predictions of extratropical cyclones. In this paper, we present a systematic quantification of CRH uncertainties. To this end, we study an idealized extratropical cyclone simulated at a convection-permitting resolution of 2.5 km and combine large-eddy-model simulations at a 300 m resolution with offline radiative transfer calculations. We quantify four factors contributing to the CRH uncertainty in different regions of the cyclone: 3D cloud-radiative effects, parameterization of ice optical properties, cloud horizontal heterogeneity, and cloud vertical overlap. The last two factors can be considered essentially resolved at 300 m but need to be parameterized at a 2.5 km resolution. Our results indicate that parameterization of ice optical properties and cloud horizontal heterogeneity are the two factors contributing most to the mean uncertainty in CRH at larger spatial scales and can be more relevant for the large-scale dynamics of the cyclone. On the other hand, 3D cloud-radiative effects are much smaller on average, especially for stratiform clouds within the warm conveyor belt of the cyclone. Our analysis in particular highlights the potential to improve the simulation of CRH by better representing ice optical properties. Future work should, in particular, address how uncertainty in ice optical properties affects the dynamics and predictability of extratropical cyclones.

Radiative Effect of Two Contrail Cirrus Outbreaks Over Western Europe Estimated Using Geostationary Satellite Observations and Radiative Transfer Calculations

AbstractEstimation of the perturbation to the Earth's energy budget by contrail outbreaks is required for estimating the climate impact of aviation and verifying the climate benefits of proposed contrail avoidance strategies such as aircraft rerouting. Here we identified two successive large‐scale contrail outbreaks developing in clear‐sky conditions in geostationary and polar‐orbiting satellite infrared images of Western Europe lasting from 22–23 June 2020. Their hourly cloud radiative effect, obtained using geostationary satellite cloud retrievals and radiative transfer calculations, is negative or weakly positive during daytime and positive during nighttime. The cumulative energy forcing of the two outbreaks is 7 PJ and −8.5 PJ, with uncertainties of 3 PJ, stemming each from approximately 15–20 flights over periods of 19 and 7 hr, respectively. This study suggests that an automated quantification of contrail outbreak radiative effect is possible, at least for contrails forming in clear sky conditions.