No Arabic abstract
We introduce a new general-purpose time-dependent ionisation network (IN) and a radiation transport (RT) module in the magneto-hydrodynamic (MHD) code PLUTO. Our ionisation network is reliable for temperatures ranging from 5e3 to 3e8 K, and includes all ionisation states of H, He, C, N, O, Ne, Mg, Si, S and Fe making it suitable for studying a variety of astrophysical scenarios. Radiation loss for each ion-electron pair is calculated using CLOUDY-17 data on-the-fly. Photo-ionisation and charge exchange are the chemical heating mechanisms. The IN is fully coupled to the radiation transport module over a very large range of opacities at different frequencies. The RT module employs a method of short characteristics assuming spherical symmetry. The radiation module requires the assumption of spherical symmetry, while the IN is compatible with full 3D. We also include a simple prescription for dust opacity, grain destruction, and the dust contribution to radiation pressure. We present numerical tests to show the reliability and limitations of the new modules. We also present a post-processing tool to calculate projected column densities and emission spectra.
We present numerical simulations of the evolution of a supernova (SN) remnant expanding into a uniform background medium with density $n_H = 1.0$ cm$^{-3}$ and temperature of $10^4$ K. We include a dynamically evolving non-equilibrium ionisation (NEI) network (consisting of all the ions of H, He, C, N, O, Ne, Mg, Si, S, Fe), frequency dependent radiation transfer (RT), thermal conduction, and a simple dust evolution model, all intra-coupled to each other and to the hydrodynamics. We assume spherical symmetry. Photo-ionisation, radiation losses, photo-heating, charge-exchange heating/cooling and radiation pressure are calculated on-the-fly depending on the local radiation field and ion fractions. We find that the dynamics and energetics (but not the emission spectra) of the SN remnants can be well modelled by collisional equilibrium cooling curves even in the absence of non-equilibrium cooling and radiative transport. We find that the effect of precursor ionising radiation at different stages of SN remnant are dominated by rapid cooling of the shock and differ from steady state shocks. The predicted column densities of different ions such as N II, C IV and N V, can be higher by up to several orders of magnitude compared to steady state shocks. We also present some high esolution emission spectra that can be compared with the observed remnants to obtain important information about the physical and chemical states of the remnant, as well as constrain the background ISM.
Persephone is a NASA concept mission study that addresses key questions raised by New Horizons encounters with Kuiper Belt objects (KBOs), with arguably the most important being Does Pluto have a subsurface ocean?. More broadly, Persephone would answer four significant science questions: (1) What are the internal structures of Pluto and Charon? (2) How have the surfaces and atmospheres in the Pluto system evolved? (3) How has the KBO population evolved? (4) What are the particles and magnetic field environments of the Kuiper Belt? To answer these questions, Persephone has a comprehensive payload, and would both orbit within the Pluto system and encounter other KBOs. The nominal mission is 30.7 years long, with launch in 2031 on a Space Launch System (SLS) Block 2 rocket with a Centaur kick stage, followed by a 27.6 year cruise powered by existing radioisotope electric propulsion (REP) and a Jupiter gravity assist to reach Pluto in 2058. En route to Pluto, Persephone would have one 50- to 100-km-class KBO encounter before starting a 3.1 Earth-year orbital campaign of the Pluto system. The mission also includes the potential for an 8-year extended mission, which would enable the exploration of another KBO in the 100- to 150-km-size class. The mission payload includes 11 instruments: Panchromatic and Color High-Resolution Imager; Low-Light Camera; Ultra-Violet Spectrometer; Near-Infrared (IR) Spectrometer; Thermal IR Camera; Radio Frequency Spectrometer; Mass Spectrometer; Altimeter; Sounding Radar; Magnetometer; and Plasma Spectrometer. The nominal cost of this mission is $3.0B, making it a large strategic science mission.
Resonant lines are powerful probes of the interstellar and circumgalactic medium of galaxies. Their transfer in gas being a complex process, the interpretation of their observational signatures, either in absorption or in emission, is often not straightforward. Numerical radiative transfer simulations are needed to accurately describe the travel of resonant line photons in real and in frequency space, and to produce realistic mock observations. This paper introduces RASCAS, a new public 3D radiative transfer code developed to perform the propagation of any resonant line in numerical simulations of astrophysical objects. RASCAS was designed to be easily customisable and to process simulations of arbitrarily large sizes on large supercomputers. RASCAS performs radiative transfer on an adaptive mesh with an octree structure using the Monte Carlo technique. RASCAS features full MPI parallelisation, domain decomposition, adaptive load-balancing, and a standard peeling algorithm to construct mock observations. The radiative transport of resonant line photons through different mixes of species (e.g. ion{H}{i}, ion{Si}{ii}, ion{Mg}{ii}, ion{Fe}{ii}), including their interaction with dust, is implemented in a modular fashion to allow new transitions to be easily added to the code. RASCAS is very accurate and efficient. It shows perfect scaling up to a minimum of a thousand cores. It has been fully tested against radiative transfer problems with analytic solutions and against various test cases proposed in the literature. Although it was designed to describe accurately the many scatterings of line photons, RASCAS may also be used to propagate photons at any wavelength (e.g. stellar continuum or fluorescent lines), or to cast millions of rays to integrate the optical depths of ionising photons, making it highly versatile.
We examine morphology-separated color-mass diagrams to study the quenching of star formation in $sim 100,000$ ($zsim0$) Sloan Digital Sky Survey (SDSS) and $sim 20,000$ ($zsim1$) Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (CANDELS) galaxies. To classify galaxies morphologically, we developed Galaxy Morphology Network (GaMorNet), a convolutional neural network that classifies galaxies according to their bulge-to-total light ratio. GaMorNet does not need a large training set of real data and can be applied to data sets with a range of signal-to-noise ratios and spatial resolutions. GaMorNets source code as well as the trained models are made public as part of this work ( http://www.astro.yale.edu/aghosh/gamornet.html ). We first trained GaMorNet on simulations of galaxies with a bulge and a disk component and then transfer learned using $sim25%$ of each data set to achieve misclassification rates of $lesssim5%$. The misclassified sample of galaxies is dominated by small galaxies with low signal-to-noise ratios. Using the GaMorNet classifications, we find that bulge- and disk-dominated galaxies have distinct color-mass diagrams, in agreement with previous studies. For both SDSS and CANDELS galaxies, disk-dominated galaxies peak in the blue cloud, across a broad range of masses, consistent with the slow exhaustion of star-forming gas with no rapid quenching. A small population of red disks is found at high mass ($sim14%$ of disks at $zsim0$ and $2%$ of disks at $z sim 1$). In contrast, bulge-dominated galaxies are mostly red, with much smaller numbers down toward the blue cloud, suggesting rapid quenching and fast evolution across the green valley. This inferred difference in quenching mechanism is in agreement with previous studies that used other morphology classification techniques on much smaller samples at $zsim0$ and $zsim1$.
In most of Seyfert-1 active galactic nucei (AGN) the optical linear continuum polarization degree is usually small (less than 1%) and the polarization position angle is nearly parallel to the AGN radio-axis. However, there are many types-1 AGNs with unexplained intermediate values for both positional angles and polarization degrees. Our explanation of polarization degree and positional angle of Seyfert-1 AGNs focuses on the reflection of non-polarized radiation from sub-parsec jets in optically thick accretion discs. The presence of a magnetic field surrounding the scattering media will induce Faraday rotation of the polarization plane that may explain the intermediate values of positional angles if there is a magnetic field component normal to the accretion disc. The Faraday rotation depolarization effect in disc diminishes the competition between polarization of the reflected radiation with the parallel component of polarization and the perpendicular polarization from internal radiation of disc (the Milne problem) in favor of polarization of reflected radiation. This effect allows us to explain the observed polarization of Seyfert-1 AGN radiation even though the jet optical luminosity is much lower than the luminosity of disc. We present the calculation of polarization degrees for a number of Seyfert-1 AGNs.