1,633 research outputs found
A WENO Algorithm of the Temperature and Ionization Profiles around a Point Source
We develop a numerical solver for radiative transfer problems based on the
weighted essentially nonoscillatory (WENO) scheme modified with anti-diffusive
flux corrections, in order to solve the temperature and ionization profiles
around a point source of photons in the reionization epoch. Algorithms for such
simulation must be able to handle the following two features: 1. the sharp
profiles of ionization and temperature at the ionizing front (I-front) and the
heating front (T-front), and 2. the fraction of neutral hydrogen within the
ionized sphere is extremely small due to the stiffness of the rate equations of
atom processes. The WENO scheme can properly handle these two features, as it
has been shown to have high order of accuracy and good convergence in capturing
discontinuities and complicated structures in fluid as well as to be
significantly superior over piecewise smooth solutions containing
discontinuities. With this algorithm, we show the time-dependence of the
preheated shell around a UV photon source. In the first stage the I-front and
T-front are coincident, and propagate with almost the speed of light. In later
stage, when the frequency spectrum of UV photons is hardened, the speeds of
propagation of the ionizing and heating fronts are both significantly less than
the speed of light, and the heating front is always beyond the ionizing front.
In the spherical shell between the I- and T-fronts, the IGM is heated, while
atoms keep almost neutral. The time scale of the preheated shell evolution is
dependent on the intensity of the photon source. We also find that the details
of the pre-heated shell and the distribution of neutral hydrogen remained in
the ionized sphere are actually sensitive to the parameters used. The WENO
algorithm can provide stable and robust solutions to study these details.Comment: 24 pages, 7 figures, accepted in New Astronom
A WENO Algorithm for the Radiative Transfer and Ionized Sphere at Reionization
We show that the algorithm based on the weighted essentially nonoscillatory
(WENO) scheme with anti-diffusive flux corrections can be used as a solver of
the radiative transfer equations. This algorithm is highly stable and robust
for solving problems with both discontinuities and smooth solution structures.
We test this code with the ionized sphere around point sources. It shows that
the WENO scheme can reveal the discontinuity of the radiative or ionizing
fronts as well as the evolution of photon frequency spectrum with high accuracy
on coarse meshes and for a very wide parameter space. This method would be
useful to study the details of the ionized patch given by individual source in
the epoch of reionization. We demonstrate this method by calculating the
evolution of the ionized sphere around point sources in physical and frequency
spaces. It shows that the profile of the fraction of neutral hydrogen and the
ionized radius are sensitively dependent on the intensity of the source.Comment: Elsart Latex file, 20 pages, 8 figures included, accepted for
publication in New Astronom
21 cm Signals from Early Ionizing Sources
We investigate the 21 cm signals from the UV ionizing sources in the
reionization epoch. The formation and evolution of 21 cm emission and
absorption regions depend essentially on the kinetics of photons in the
physical and frequency spaces. To solve the radiative transfer equation, we use
the WENO algorithm, which is effective to capture the sharp ionization profile
and the cut-off at the front of light and to handle the small fraction
of neutral hydrogen and helium in the ionized sphere. We show that a spherical
shell of 21 cm emission and absorption will develop around a point source once
the speed of the ionization front (I-front) is significantly lower than the
speed of light. The 21 cm shell extends from the I-front to the front of light;
its inner part is the emission region and its outer part is the absorption
region. The 21 cm emission region depends strongly on the intensity,
frequency-spectrum and life-time of the UV ionizing source. For a source of
short life-time, no 21 cm emission region can be formed if the source dies out
before the I-front speed is significantly lower than the speed of light. Yet, a
21 cm absorption region can form and develop even after the emission of the
source ceases.Comment: 25 pages, 9 figures, accepted by Ap
Facile hydrothermal synthesis and optical limiting properties of TiO 2 -reduced graphene oxide nanocomposites
TiO2/reduced graphene oxide (RGO) nanocomposites Gx (RGO titania nanocomposite, x grams tetrabutyl titanate per 0.03 g RGO, x = 0.25, 0.50, 1.00) were prepared by a hydrothermal method: graphene oxide was reduced to RGO in a 2:1 water:ethanol mixture in the presence of varying quantities of tetrabutyl titanate, which deposited as TiO2 on the RGO sheets. The nanocomposites were characterized by a combination of Fourier transform infrared spectroscopy, diffuse reflectance ultraviolet-visible spectroscopy, photoluminescence spectroscopy, Raman spectroscopy, X-ray powder diffraction, X-ray photoelectron spectroscopy and transmission electron microscopy studies. The nanocomposite G0.25 exhibits enhanced nonlinear optical properties compared to its individual components, which is ascribed to a combination of mechanisms. The role of defects and electron/energy transfer in the optical limiting performance of G0.25 was clarified with the help of Raman and photoluminescence spectroscopies. Intensity-dependent switching between reverse saturable absorption and saturable absorption behavior was observed with the G0.50 nanocomposite
PHoToNs–A parallel heterogeneous and threads oriented code for cosmological N-body simulation
We introduce a new code for cosmological simulations, PHoToNs, which incorporates features for performing massive cosmological simulations on heterogeneous high performance computer (HPC) systems and threads oriented programming. PHoToNs adopts a hybrid scheme to compute gravitational force, with the conventional Particle-Mesh (PM) algorithm to compute the long-range force, the Tree algorithm to compute the short range force and the direct summation Particle-Particle (PP) algorithm to compute gravity from very close particles. A self-similar space filling a Peano-Hilbert curve is used to decompose the computing domain. Threads programming is advantageously used to more flexibly manage the domain communication, PM calculation and synchronization, as well as Dual Tree Traversal on the CPU+MIC platform. PHoToNs scales well and efficiency of the PP kernel achieves 68.6% of peak performance on MIC and 74.4% on CPU platforms. We also test the accuracy of the code against the much used Gadget-2 in the community and found excellent agreement
- …