An implementation of radiative transfer in the cosmological simulation code gadget

We present a novel numerical implementation of radiative transfer in the cosmological smoothed particle hydrodynamics (SPH) simulation code gadget . It is based on a fast, robust and photon-conserving integration scheme where the radiation transport problem is approximated in terms of moments of the transfer equation and by using a variable Eddington tensor as a closure relation, following the Optically Thin Variable Eddington Tensor suggestion of Gnedin & Abel. We derive a suitable anisotropic diffusion operator for use in the SPH discretization of the local photon transport, and we combine this with an implicit solver that guarantees robustness and photon conservation. This entails a matrix inversion problem of a huge, sparsely populated matrix that is distributed in memory in our parallel code. We solve this task iteratively with a conjugate gradient scheme. Finally, to model photon sink processes we consider ionization and recombination processes of hydrogen, which is represented with a chemical network that is evolved with an implicit time integration scheme. We present several tests of our implementation, including single and multiple sources in static uniform density fields with and without temperature evolution, shadowing by a dense clump and multiple sources in a static cosmological density field. All tests agree quite well with analytical computations or with predictions from other radiative transfer codes, except for shadowing. However, unlike most other radiative transfer codes presently in use for studying re-ionization, our new method can be used on-the-fly during dynamical cosmological simulation, allowing simultaneous treatments of galaxy formation and the re-ionization process of the Universe.     

The escape of ionizing photons from supernova-dominated primordial galaxies

In order to assess the contribution of Lyman break galaxies (LBGs) and Lyman α emitters (LAEs) at redshifts  3 < z < 7  to the ionization of intergalactic medium (IGM), we investigate the escape fractions of ionizing photons from supernova-dominated primordial galaxies by solving the three-dimensional (3D) radiative transfer. The model galaxy is employed from an ultra-high-resolution chemodynamic simulation of a primordial galaxy by Mori & Umemura, which well reproduces the observed properties of LAEs and LBGs. The total mass of model galaxy is  10¹¹ M_⊙  . We solve not only photoionization but also collisional ionization by shocks. In addition, according to the chemical enrichment, we incorporate the effect of dust extinction, taking the size distributions of dust into account. As a result, we find that dust extinction reduces the escape fractions by a factor of 1.5–8.5 in the LAE phase and by a factor of 2.5–11 in the LBG phase, while the collisional ionization by shocks increases the escape fractions by a factor of  ≈2  . The resultant escape fractions are 0.07–0.47 in the LAE phase and 0.06–0.17 in the LBG phase. These results are well concordant with the recent estimations derived from the flux density ratio at 1500 to 900 Å of LAEs and LBGs. Combining the resultant escape fractions with the luminosity functions of LAEs and LBGs, we find that high- z LAEs and LBGs can ionize the IGM at   z = 3–5  . However, ionizing radiation from LAEs as well as LBGs falls short of ionizing the IGM at   z > 6  . That implies that additional ionization sources may be required at   z > 6  .     

Prompt optical emission and synchrotron self-absorption constraints on emission site of GRBs

We constrain the distance of the gamma-ray burst (GRB) prompt emission site from the explosion centre R , by determining the location of the electron's self-absorption frequency in the GRB prompt optical-to-X/γ-ray spectral energy distribution, assuming that the optical and the γ-ray emissions are among the same synchrotron radiation continuum of a group of hot electrons. All possible spectral regimes are considered in our analysis. The method has only two assumed parameters, namely the bulk Lorentz factor of the emitting source Γ and the magnetic field strength B in the emission region (with a weak dependence). We identify a small sample of four bursts that satisfy the following three criteria: (1) they all have simultaneous optical and γ-ray detections in multiple observational time intervals, (2) they all show temporal correlations between the optical and γ-ray light curves and (3) the optical emission is consistent with belonging to the same spectral component as the γ-ray emission. For all the time intervals of these four bursts, it is inferred that   R ≥ 10¹⁴  (Γ/300)^3/4 ( B /10⁵ G)^1/4  cm. For a small fraction of the sample, the constraint can be pinned down to   R ≈ 10¹⁴–10¹⁵ cm  for  Γ∼ 300  . For a second sample of bursts with prompt optical non-detections, only upper limits on R can be obtained. We find no inconsistency between the R -constraints for this non-detection sample and those for the detection sample.     

Lyman α radiative transfer in a multiphase medium

Hydrogen Lyman α (Lyα) is our primary emission-line window into high-redshift galaxies. Despite an extensive literature, Lyα radiative transfer in the most realistic case of a dusty, multiphase medium has received surprisingly little detailed theoretical attention. We investigate Lyα resonant scattering through an ensemble of dusty, moving, optically thick gas clumps. We treat each clump as a scattering particle and use Monte Carlo simulations of surface scattering to quantify continuum and Lyα surface scattering angles, absorption probabilities, and frequency redistribution, as a function of the gas dust content. This atomistic approach speeds up the simulations by many orders of magnitude, making possible calculations which are otherwise intractable. Our fitting formulae can be readily adapted for fast radiative transfer in numerical simulations. With these surface scattering results, we develop an analytic framework for estimating escape fractions and line widths as a function of gas geometry, motion, and dust content. Our simple analytic model shows good agreement with full Monte Carlo simulations. We show that the key geometric parameter is the average number of surface scatters for escape in the absence of absorption,     , and we provide fitting formulae for several geometries of astrophysical interest. We consider the following two interesting applications. (i) Equivalent widths ( EWs ). Lyα can preferentially escape from a dusty multiphase interstellar medium if most of the dust lies in cold neutral clouds, which Lyα photons cannot penetrate. This might explain the anomalously high EWs sometimes seen in high-redshift/submillimetre sources. (ii) Multiphase galactic outflows . We show the characteristic profile is asymmetric with a broad red tail, and relate the profile features to the outflow speed and gas geometry. Many future applications are envisaged.     

康普顿化的计算-推广的Kompaneets方程与Monte Carlo方法比较

在X射线天文学以及辐射物理学中,当硬X射线穿过" 冷"的等离子体时所发生的Compton软化是一种重要的辐射转移过程.简要介绍推广的Kompaneets方程,该方程在hv〈mec2及kTe〈mec2较宽松的条件下广泛成立,因此不仅能处理Compton硬化过程,而且也适用于Compton软化过程,后者在目前快速发展的X射线和γ射线天文学中十分重要.基于此方程,我们对天体物理中4种常见辐射谱(Gauss型发射谱线、黑体辐射谱、幂律辐射谱和热轫致辐射谱)在Compton软化情况下的谱演化进行了数值求解,并同Monte Carlo模拟结果做比较,证实了推广的Kompaneets方程的正确性和有效性.最后指出此方程在X射线天文学和γ射线天文学中的重要意义和潜在应用.     

The hydrogen emission spectrum of a shock wave in a stellar atmosphere

Based on a self-consistent solution of the equations of gas dynamics, kinetics of hydrogen atomic level populations, and radiative transfer, we analyze the structure of a shock wave that propagates in a partially ionized hydrogen gas. We consider the radiative transfer at the frequencies of spectral lines by taking into account the effects of a moving medium in the observer's frame of reference. The flux in Balmer lines is shown to be formed behind the shock discontinuity at the initial hydrogen recombination stage. The Doppler shift of the emission-line profile is approximately one and a half times smaller than the gas flow velocity in the Balmer emission region, because the radiation field of the shock wave is anisotropic. At Mach numbers M₁?10 and unperturbed gas densities σ₁=10^?10 g cm^?3, the Doppler shift is approximately one third of the shock velocity U₁. The FWHM of the emission-line profile δ _? is related to the shock velocity by δ _? ≈ k _? U₁, where k _? = 1, 0.6, and 0.65 for the Hα, Hβ, and Hγ lines, respectively.     

2D/3D Dust Continuum Radiative Transfer Codes to Analyze and Predict VLTI Observations

Radiative Transfer (RT) codes with image capability are a fundamental tool for preparing interferometric observations and for interpreting visibility data. In view of the upcoming VLTI facilities, we present the first comparison of images/visibilities coming from two 3D codes that use completely different techniques to solve the problem of self-consistent continuum RT. In addition, we focus on the astrophysical case of a disk distorted by tidal interaction with by-passing stars or internal planets and investigate for which parameters the distortion can be best detected in the mid-infrared using the mid-infrared interferometric device MIDI. This revised version was published online in July 2006 with corrections to the Cover Date.     

Solar Stokes Spectrum Observation at Yunnan Astronomical Observatory

1　ＳｉｇｎｉｆｉｃａｎｃｅｏｆＳｏｌａｒＳｔｏｋｅｓＳｐｅｃｔｒｕｍＯｂｓｅｒｖａｔｉｏｎＤｕｅｔｏｔｈｅｄｅｖｅｌｏｐｍｅｎｔｉｎｏｂｓｅｒｖｉｎｇｉｎｓｔｒｕｍｅｎｔ,ｔｈｅｃｏｎｃｅｐｔｏｆｓｏｌａｒｓｐｅｃｔｒｕｍｓｈｏｕｌｄｂｅｇｅｎｅｒａｌｉｚｅｄ .Ｉｎｔｈｅｐａｓｔ,ｉｔｉｓｏｎｌｙａｓｓｉｇｎｅｄｔｏｔｈｅｏｒｄｉｎａｒｙｏｒｕｎｐｏｌａｒｉｚｅｄｓｐｅｃｔｒｕｍｒｅｃｏｒｄｉｎｇｔｈｅｗａｖｅｌｅｎｇｔｈ ｄｉｓｔｒｉｂｕｔｅｄｉｎｔｅｎｓｉｔｙｏｆｔｈｅｃｏｎｔｉｎｕｕｍａｎｄ…