Source author record

Li-Zhi Fang

Li-Zhi Fang appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

16works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

16 published item(s)

preprint2011arXiv

Effect of Dust on Lyman-alpha Photon Transfer in Optically Thick Halo

We investigate the effects of dust on Lyα photons emergent from an optically thick medium by solving the integro-differential equation of the radiative transfer of resonant photons. To solve the differential equations numerically we use the Weighted Essentially Non-oscillatory method (WENO). Although the effects of dust on radia-tive transfer is well known, the resonant scattering of Lyα photons makes the problem non-trivial. For instance, if the medium has the optical depth of dust absorption and scattering to be τa>>1, τ>>1, and τ>>τa, the effective absorption optical depth in a random walk scenario would be equal to \surd τa(τa+τ). We show, however, that for a resonant scattering at frequency ν0, the effective absorption optical depth would be even larger than τ(ν0). If the cross section of dust scattering and absorption is frequency-independent, the double-peaked structure of the frequency profile given by the resonant scattering is basically dust-independent. That is, dust causes neither narrowing nor widening of the width of the double peaked profile. One more result is that the time scales of the Lyα photon transfer in the optically thick halo are also basically independent of the dust scattering, even when the scattering is anisotropic. This is because those time scales are mainly determined by the transfer in the frequency space, while dust scattering, either isotropic or anisotropic, does not affect the behavior of the transfer in the frequency space when the cross section of scattering is wavelength-independent. This result does not support the speculation that dust will lead to the smoothing of the brightness distribution of Lyα photon source with optical thick halo.

preprint2011arXiv

Intermittence of the Map of Kinetic Sunyaev-Zel'dovich Effect and Turbulence of IGM

We investigate the possibility of detecting the turbulent state of the IGM with the kinetic Sunyaev-Zel'dovich (kSZ) effect. Being sensitive to the divergence-free component of the momentum field of the IGM, the kSZ effect might be used to probe the vorticity of the turbulent IGM. With cosmological hydrodynamical simulation in the concordance $Λ$CDM universe, we find that the structure functions of 2D kSZ maps show strong intermittence, and the intermittent exponents follow a law similar to the She-Leveque scaling formula of fully developed turbulence. We also find that the intermittence is weak in the maps of thermal Sunyaev-Zel'dovich (tSZ) effect. Nevertheless, the superposition of the kSZ and tSZ effects still contain significant intermittence. We conclude that the turbulent behavior of the IGM may be revealed by the observation of SZ effect on angular scales equal to or less than 0.5 arcminute, corresponding to the multipole parameter $l\geq 2 \times10^4$.

preprint2011arXiv

Time-Dependent Behavior of Lyman$α$ Photon Transfer in High Redshift Optically Thick Medium

With Monte Carlo simulation method, we investigate the time dependent behavior of Ly$α$ photon transfer in optically thick medium of the concordance $Λ$CDM universe. At high redshift, the Ly$α$ photon escaping from optically thick medium has a time scale as long as the age of the luminous object, or even comparable to the age of the universe. In this case, time-independent, or stationary solutions of the Ly$α$ photon transfer with resonant scattering will overlook important features of the escaped Ly$α$ photons in physical and frequency spaces. More seriously, the expansion of the universe leads to that the time-independent solutions of the Ly$α$ photon transfer may not exist. We show that time-dependent solutions sometimes are essential for understanding the Ly$α$ emission and absorption at high redshifts. For Lyαphotons from sources at redshift 1+z=10 and being surrounded by neutral hydrogen IGM of the $Λ$CDM universe, the escape coefficient is found to be always less, or much less than one, regardless of the age or life time of the sources. Under such environment, we also find that even when the Ly$α$ photon luminosity of the sources is stable, the mean surface brightness is gradually increasing in the first 10^6 years, and then decreasing with a power law of time, but never approaches a stable, time-independent state. That is, all 1+z=10 sources in a neutral Hubble expanding IGM with Ly$α$ luminosity L have their maximum of mean surface brightness ~ 10^{-21}(L/(10^{43}erg/s)) erg s^{-1} cm^{-2} arcsec^{-2} at the age of about 10^6 years. The time-dependent effects on the red damping wing profile are also addressed.

preprint2010arXiv

Log-Poisson Non-Gaussianity of Ly$α$ Transmitted Flux Fluctuations at High Redshift

We investigate the non-Gaussian features of the IGM at redshift $z\sim 5 - 6$ using Ly$α$ transmitted flux of quasar absorption spectra and cosmological hydrodynamic simulation of the concordance $Λ$CDM universe. We show that the neutral hydrogen mass density field and Ly$α$ transmitted flux fluctuations possess all the non-Gaussian features predicted by the log-Poisson hierarchy, which depends only on two dimensionless parameters $β$ and $γ$, describing, respectively, the intermittence and singularity of the random fields. We find that the non-Gaussianity of the Ly$α$ transmitted flux of quasars from $z=4.9$ to $z=6.3$ can be well reconstructed by the hydrodynamical simulation samples. Although the Gunn-Peterson optical depth and its variance underwent a significant evolution in the redshift range of $5 - 6$, the intermittency measured by $β$ is almost redshift-independent in this range. More interesting, the intermittency of quasar's absorption spectra on physical scales $0.1-1$ h$^{-1}$Mpc in redshift $5 - 6$ are found to be about the same as that on physical scales $1-10$ h$^{-1}$Mpc at redshifts $2 - 4$. Considering the Jeans length is less than 0.1 h$^{-1}$Mpc at $z\sim 5$, and $1$ h$^{-1}$Mpc at $z\sim 2$, these results imply that the nonlinear evolution in high and low redshifts will lead the cosmic baryon fluid to a state similar to fully developed turbulence. The log-Poisson high order behavior of current high redshift data of quasar's spectrum can be explained by uniform UV background in the redshift range considered. We also studied the log-Poisson non-Gaussianity by considering inhomogeneous background. With several simplified models of inhomogeneous background, we found the effect of the inhomogeneous background on the log-Poisson non-Gaussianity is not larger than 1-sigma.

preprint2010arXiv

Resonant Scattering and Ly-alpha Radiation Emergent from Neutral Hydrogen Halos

With a state-of-the-art numerical method for solving the integral-differential equation of radiative transfer, we investigate the flux of the Ly$α$ photon $ν_0$ emergent from an optically thick halo containing a central light source. Our focus is on the time-dependent effects of the resonant scattering. We first show that the frequency distribution of photons in the halo are quickly approaching to a locally thermalized state around the resonant frequency, even when the mean intensity of the radiation is highly time-dependent. Since initial conditions are forgotten during the thermalization, some features of the flux, such as the two peak structure of its profile, actually are independent of the intrinsic width and time behavior of the central source, if the emergent photons are mainly from photons in the thermalized state. In this case, the difference $|ν_{\pm}-ν_0|$, where $ν_{\pm}$ are the frequencies of the two peaks of the flux, cannot be less than $2$ times of Doppler broadening. We then study the radiative transfer in the case where the light emitted from the central source is a flash. We calculate the light curves of the flux from the halo. It shows that the flux is still a flash. The time duration of the flash for the flux, however, is independent of the original time duration of the light source but depends on the optical depth of the halo. Therefore, the spatial transfer of resonant photons is a diffusion process, even though it is not a purely Brownian diffusion. This property enables an optically thick halo to trap and store thermalized photons around $ν_0$ for a long time after the cease of the central source emission. The photons trapped in the halo can yield delayed emission, of which the profile also shows typical two peak structure as that from locally thermalized photons. Possible applications of these results are addressed.

preprint2010arXiv

Statistical and dynamical decoupling of the IGM from Dark Matter

The mean mass densities of cosmic dark matter is larger than that of baryonic matter by a factor of about 5 in the $Λ$CDM universe. Therefore, the gravity on large scales should be dominant by the distribution of dark matter in the universe. However, a series of observations incontrovertibly show that the velocity and density fields of baryonic matter are decoupling from underlying dark matter field. This paper shows our attemps to unveil the physics behind this puzzle. In linear approximation, the dynamics of the baryon fluid is completely governed by the gravity of the dark matter. Consequently, the mass density field of baryon matter $ρ_b({\bf r},t)$ will be proportional to that of dark matter $ρ_{\rm dm}({\bf r},t)$, even though they are different from each other initially. In weak and moderate nonlinear regime, the dynamics of the baryon fluid can be sketched by Burgers equation. A basic feature of the Burgers dynamics is to yield shocks. When the Reynolds number is large, the Burgers fluid will be in the state of Burgers turbulence, which consists of shocks and complex structures. On the other hand, the collisionless dark matter may not show such shock, but a multivalued velocity field. Therefore, the weak and moderate nonlinear evolution leads to the IGM-dark matter deviation. Yet, the velocity field of Burgers fluid is still irrotational, as gravity is curl-free. In fully nonlinear regime, the vorticity of velocity field developed, and the cosmic baryonic fluid will no longer be potential, as the dynamics of vorticity is independent of gravity and can be self maintained by the nonlinearity of hydrodynamics. In this case, the cosmic baryon fluid is in the state of fully developed turbulence, which is statistically and dynamically decoupling from dark matter. This scenario provides a mechanism of cohenent explanation of observations.

preprint2010arXiv

Vorticity of IGM Velocity Field on Large Scales

We investigate the vorticity of the IGM velocity field on large scales with cosmological hydrodynamic simulation of the concordance model of LCDM. We show that the vorticity field is significantly increasing with time as it can effectively be generated by shocks and complex structures in the IGM. Therefore, the vorticity field is an effective tool to reveal the nonlinear behavior of the IGM, especially the formation and evolution of turbulence in the IGM. We find that the vorticity field does not follow the filaments and sheets structures of underlying dark matter density field and shows highly non- Gaussian and intermittent features. The power spectrum of the vorticity field is used to measure the development of turbulence in Fourier space. We show that the relation between the power spectra of vorticity and velocity fields is perfectly in agreement with the prediction of a fully developed homogeneous and isotropic turbulence from 0.2 to 3 h^{-1} Mpc at z~0. This indicates that cosmic baryonic field is in the state of fully developed turbulence on scales less than about 3 h^{-1} Mpc. The random field of the turbulent fluid yields turbulent pressure to prevent the gravitational collapsing of the IGM. The vorticity and turbulent pressure are strong inside and even outside of high density regions. In IGM regions with 10 times mean overdensity, the turbulent pressure can be on an average equivalent to the thermal pressure of the baryonic gas with a temperature of 10^5 K. The fully developed turbulence would prevent the baryons in the IGM from falling into the gravitational well of dark matter halos. Moreover, turbulent pressure is dynamical and non-thermal, which makes it different from pre-heating mechanism as it does not affect the thermal state and ionizing process of hydrogen in the IGM.

preprint2009arXiv

A wavelet-Galerkin algorithm of the E/B decomposition of CMB polarization maps

We develop an algorithm of separating the $E$ and $B$ modes of the CMB polarization from the noisy and discretized maps of Stokes parameter $Q$ and $U$ in a finite area. A key step of the algorithm is to take a wavelet-Galerkin discretization of the differential relation between the $E$, $B$ and $Q$, $U$ fields. This discretization allows derivative operator to be represented by a matrix, which is exactly diagonal in scale space, and narrowly banded in spatial space. We show that the effect of boundary can be eliminated by dropping a few DWT modes located on or nearby the boundary. This method reveals that the derivative operators will cause large errors in the $E$ and $B$ power spectra on small scales if the $Q$ and $U$ maps contain Gaussian noise. It also reveals that if the $Q$ and $U$ maps are random, these fields lead to the mixing of the $E$ and $B$ modes. Consequently, the $B$ mode will be contaminated if the powers of $E$ modes are much larger than that of $B$ modes. Nevertheless, numerical tests show that the power spectra of both $E$ and $B$ on scales larger than the finest scale by a factor of 4 and higher can reasonably be recovered, even when the power ratio of $E$- to $B$-modes is as large as about 10$^2$, and the signal-to-noise ratio is equal to 10 and higher. This is because the Galerkin discretization is free of false correlations, and keeps the contamination under control. As wavelet variables contain information of both spatial and scale spaces, the developed method is also effective to recover the spatial structures of the $E$ and $B$ mode fields.

preprint2008arXiv

Log-Poisson Hierarchical Clustering of Cosmic Neutral Hydrogen and Ly-alpha Transmitted Flux of QSO Absorption Spectrum

we study, in this paper, the non-Gaussian features of the mass density field of neutral hydrogen fluid and the Ly-alpha transmitted flux of QSO absorption spectrum from the point-of-view of self-similar log-Poisson hierarchy. It has been shown recently that, in the scale range from the onset of nonlinear evolution to dissipation, the velocity and mass density fields of cosmic baryon fluid are extremely well described by the She-Leveque's scaling formula, which is due to the log-Poisson hierarchical cascade. Since the mass density ratio between ionized hydrogen to total hydrogen is not uniform in space, the mass density field of neutral hydrogen component is not given by a similar mapping of total baryon fluid. Nevertheless, we show, with hydrodynamic simulation samples of the concordance $Λ$CDM universe, that the mass density field of neutral hydrogen, is also well described by the log-Poisson hierarchy. We then investigate the field of Ly$α$ transmitted flux of QSO absorption spectrum. Due to redshift distortion, Ly$α$ transmitted flux fluctuations are no longer to show all features of the log-Poisson hierarchy. However, some non-Gaussian features predicted by the log-Poisson hierarchy are not affected by the redshift distortion. We test these predictions with the high resolution and high S/N data of quasars Ly$α$ absorption spectra. All results given by real data, including $β$-hierarchy, high order moments and scale-scale correlation, are found to be well consistent with the log-Poisson hierarchy. We compare the log-Poisson hierarchy with the popular log-normal model of the Ly$α$ transmitted flux. The later is found to yield too strong non-Gaussianity at high orders, while the log-Poisson hierarchy is in agreement with observed data.

preprint2007arXiv

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 $(r=ct)$ 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.

preprint2005arXiv

The velocity field of baryonic gas in the universe

The dynamic evolution of the baryonic intergalactic medium (IGM) caused by the underlying dark matter gravity is governed by the Navier-Stokes equations in which many cooling and heating processes are involved. However, it has long been recognized that the growth mode dynamics of cosmic matter clustering can be sketched by a random force driven Burgers' equation if cooling and heating are ignored. Just how well the dynamics of the IGM can be described as a Burgers fluid has not been fully investigated probably because cooling and heating are essential for a detailed understanding of the IGM. Using IGM samples produced by a cosmological hydrodynamic simulation in which heating and cooling processes are properly accounted for, we show that the IGM velocity field in the nonlinear regime shows the features of a Burgers fluid, that is, when the Reynolds number is high, the velocity field consists of an ensemble of shocks. Consequently, (1) the IGM velocity $v$ is generally smaller than that of dark matter; (2) for the smoothed field, the IGM velocity shows tight correlation with dark matter given by $v \simeq s v_{dm}$, with $s<1$, such that the lower the redshift, the smaller $s$; (3) the velocity PDFs are asymmetric between acceleration and deceleration events; (4) the PDF of velocity difference $Δv=v(x+r)-v(x)$ satisfies the scaling relation for a Burgers fluid, i.e., $P(Δv)=(1 r^y)F(Δv/r^y)$. We find the scaling function and parameters for the IGM which are applicable to the entire scale range of the samples (0.26 - 8 h$^{-1}$ Mpc). These properties show that the similarity mapping between the IGM and dark matter is violated on scales much larger than the Jeans length of the IGM.

preprint2001arXiv

The large-scale rms bulk velocity estimated from QSOs' Ly$α$ forests

We propose a method for estimating the large-scale rms bulk velocity of the cosmic mass field from the transmitted fluxes of Ly$α$ forests. It is based on two linear relationships on large scales: 1) the relation between the fluctuations of the transmission and the underlying density field, and 2) the relation between the density fluctuations and the peculiar velocity field. We show that, with a multiscale decomposition, the two relations can be effectively employed for predicting the rms bulk velocity. Since QSO's Ly$α$ forest is due to the absorptions of diffusely distributed and photoionized IGM, this method provides an independent estimate of the rms bulk velocity at high redshifts, on large scales, and free from the bias of galaxies. Using the transmitted flux of 60 moderate-resolution QSO spectra, the rms bulk velocity is found to be 230$\pm$50 km s$^{-1}$ around redshift $z=2.25$ on scale 23 h$^{-1}$Mpc, and down to 110$\pm$45 km s$^{-1}$ around $z=3.25$ on scale 92 h$^{-1}$Mpc for an LCDM universe ($Ω=0.3$ and $Λ=0.7$). The results are basically consistent with the linear evolution theory.

preprint2000arXiv

Intermittent behavior of cosmic mass field revealed by QSO's Ly_alpha forests

The intermittent behavior of the space-scale distribution of Ly$α$ transmitted flux of QSO HS1700+64 has been analyzed via a discrete wavelet transform. We found that there are strong indications of intermittency on scales down to about 10 $h^{-1}$ kpc. These are: 1.) the probability distribution function of the local fluctuations of the flux is significantly long-tailed on small scales, and 2.) the local power spectrum of the flux shows prominent spiky structures on small scales. Moreover, the local power spectrum averaged on regions with different sizes shows similar spiky structures. Therefore, the random mass density field traced by the Ly$α$ forests is rougher on smaller scales, consistent with singular clustering.

preprint1997arXiv

Large-scale Structures revealed by Wavelet Decomposition

We present a detailed review of large-scale structure (LSS) study using the discrete wavelet transform (DWT). After describing how one constructs a wavelet decomposition we show how this bases can be used as a complete statistical discription of LSS. Among the topics studied are the the DWT estimation of the probability distribution function; the reconstruction of the power spectrum; the regularization of complex geometry in observational samples; cluster identification; extraction and identification of coherent structures; scale-decomposition of non-Gaussianity, such as spectra of skewnes and kurtosis and scale-scale correlations. These methods are applied to both observational and simulated samples of the QSO Lyman-alpha forests. It is clearly demonstrated that the statistical measures developed using the DWT are needed to distinguish between competing models of structure formation. The DWT also reveals physical features in these distributions not detected before. We conclude with a look towards the future of the use of the DWT in LSS.

preprint1995arXiv

Thermally Induced Density Perturbations in the Inflation Era

The possibility of thermally induced initial density perturbations in inflationary cosmology is examined. The fluctuation dynamics of a scalar field plus thermal bath system during slow roll is described by a Langevin-like equation. Fluctuation-dissipation arguments show that for a wide parameter range within the standard inflation model, the thermal fluctuations of the scalar field can dominate its quantum fluctuations. The initial amplitude of density perturbations is found to lie in a range which is consistent with the recent observations of cosmic temperature fluctuations.

preprint1994arXiv

On the measurement of the Hubble constant in a local low-density universe

Astrophysical observations indicate that the ``Local Universe" has a relatively lower matter density ($Ω_0$) than the predictions of the standard inflation cosmology and the large-scale motions of galaxies which provide a mean mass density to be very close to unity. In such a local underdense region the Hubble expansion may not be representative of the global behaviour. Utilizing an underdense sphere embedded in a flat universe as the model of our ``Local Universe", we show that the local Hubble constant would be 1.2 -- 1.4 times larger than the global value on scale of $\sim80$ Mpc, depending on the variation of $Ω_0$. This may account for the recent measurements of the unpleasantly large Hubble constant of $\sim$80 km/s/Mpc using the Cepheid variables in the Virgo cluster and the relative distance between Virgo and Coma cluster and removes the resulted apparent paradox of the age of our universe.