Source author record

Jun Pan

Jun Pan appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.CO astro-ph astro-ph.GA chem-ph cond-mat physics.chem-ph

Catalog footprint

What is connected

16works

6topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2022arXiv

About One-point Statistics of the Ratio of Two Fourier-transformed Cosmic Fields and an Application

The Fourier transformation is an effective and efficient operation of Gaussianization at the one-point level. Using a set of N-body simulation data, we verified that the one-point distribution functions of the dark matter momentum divergence and density fields closely follow complex Gaussian distributions. The one-point distribution function of the quotient of two complex Gaussian variables is introduced and studied. Statistical theories are then applied to model one-point statistics about the growth of individual Fourier mode of the dark matter density field, which can be obtained by the ratio of two Fourier transformed cosmic fields. Our simulation results proved that the models based on the Gaussian approximation are impressively accurate, and our analysis revealed many interesting aspects about the growth of dark matter's density fluctuation in Fourier space.

preprint2021arXiv

Numerical investigation of non-Gaussianities in the phase and modulus of density Fourier modes

We numerically investigate non-Gaussianities in the late-time cosmological density field in Fourier space. We explore various statistics, including the two-point and three-point probability distribution function (PDF) of phase and modulus, and two \& three-point correlation function of of phase and modulus. We detect significant non-Gaussianity for certain configurations. We compare the simulation results with the theoretical expansion series of \citet{2007ApJS..170....1M}. We find that the $\mathcal{O}(V^{-1/2})$ order term alone is sufficiently accurate to describe all the measured non-Gaussianities in not only the PDFs, but also the correlations. We also numerically find that the phase-modulus cross-correlation contributes $\sim 50\%$ to the bispectrum, further verifying the accuracy of the $\mathcal{O}(V^{-1/2})$ order prediction. This work demonstrates that non-Gaussianity of the cosmic density field is simpler in Fourier space, and may facilitate the data analysis in the era of precision cosmology.

preprint2020arXiv

Simulating Kilonovae in the ΛCDM Universe

Kilonovae are optical flashes produced in the aftermath of neutron star-neutron star mergers (NNMs) or neutron star-black hole mergers (NBMs). In this work, we use the Millennium Simulation, combined with a semi-analytic galaxy formation model--GABE (Galaxy Assembly with Binary Evolution) which adopts binary stellar population synthesis models, to explore the cosmic event rate of kilonovae, and the properties of their host galaxies in a cosmological context. We find that model with supernova kick velocity of 0 km/s fits the observation best, in agreement with the exception of some formation channels of binary neutron star. This indicates that NNMs prefer to originate from binary systems with low kick velocities. With V$_{\rm kick}$=0 km/s, the cosmic event rate of NNMs and NBMs at z=0 are 283 Gpc$^{-3}$yr$^{-1}$ and 91 Gpc$^{-3}$yr$^{-1}$, respectively, marginally consistent with the constraint from LVC GWTC-1. For Milky Way-mass galaxies, we predict the NNM rate is $25.7^{+59.6}_{-7.1}$ Myr$^{-1}$, which is also in good agreement with the observed properties of binary neutron stars in the Milky Way. Taking all the NNMs into account in the history of Milky Way-mass galaxies, we find that the averaged r-process elements yield with A>79 in a NNM and NBM event should be 0.01 M$_{\odot}$ to be consistent with observation. We conclude that NGC 4993, the host galaxy of GW170817, is a typical host galaxy for NNMs. However, generally NNMs and NBMs tend to reside in young, blue, star-forming, late-type galaxies, with stellar mass and gaseous metallicity distribution peaking at $10^{10.65}$ M$_{\odot}$ and 8.72-8.85, respectively. By studying kilonovae host galaxies in the cosmological background, it is promising to constrain model details better when we have more events in the forthcoming future. (abridged)

preprint2013arXiv

An estimation of local bulk flow with the maximum-likelihood method

A maximum-likelihood method, tested as an unbiased estimator from numerical simulations, is used to estimate cosmic bulk flow from peculiar velocity surveys. The likelihood function is applied to four observational catalogues (ENEAR, SFI++, A1SN and SC) constructed from galaxy peculiar velocity surveys and Type-Ia supernovae data at low redshift ($z \leq 0.03$). We find that the Spiral Field {\it I}-band catalogue constrains the bulk flow to be $V=290 \pm 30 {\,{\rm km}\,{\rm s}^{-1}}$ towards $l=281^{\circ} \pm 7^{\circ}$, $b=8^{\circ +6^{\circ}}_{-5^{\circ}}$ on effective scales of $ 58 {\,{h^{-1} {\rm Mpc}}}$, which is the tightest constraints achievable at the present time. By comparing the amplitudes of our estimated bulk flows with theoretical prediction, we find excellent agreement between the two. In addition, directions of estimated bulk flows are also consistent with measurements in other studies.

preprint2013arXiv

Peculiar Velocity Decomposition, Redshift Space Distortion and Velocity Reconstruction in Redshift Surveys - I. The Methodology

Massive spectroscopic surveys will measure the redshift space distortion (RSD) induced by galaxy peculiar velocity to unprecedented accuracy. We develop a new method to improve the RSD modeling and to carry out robust reconstruction of the 3D large scale peculiar velocity through galaxy redshift surveys, in light of RSD. (1) We propose a mathematically unique and physically motivated decomposition of peculiar velocity into three eigen-components of different origins, different scale dependences and different impacts on RSD. (2) This decomposition has the potential to simplify and improve the RSD modeling. We derive a new formula for the redshift space power spectrum. Under the velocity decomposition scheme, all high order Gaussian corrections and non-Gaussian correction of order $δ^3$ can be taken into account without introducing extra model uncertainties. We find that the leading order contribution $\propto (1+f\tilde{W}(k)u^2)^2$ differs from the Kaiser formula by a function $\tilde{W}(k)\leq 1$. We hence identify a significant systematic error causing underestimation of the structure growth parameter $f$ by as much as $O(10%)$ even at relatively large scale $k=0.1h/$Mpc. The velocity decomposition also reveals the three origins of the finger of God (FOG) effect. (3) The velocity decomposition clarifies issues in peculiar velocity reconstruction through 3D galaxy distribution. We discuss two possible ways to carry out the velocity reconstruction. Both use the otherwise troublesome RSD in velocity reconstruction as a valuable source of information. Both have the advantage to render the reconstruction of a stochastic 3D field into the reconstruction of a deterministic window function $W^s(k,u)$ of limited degrees of freedom...[abridged]

preprint2012arXiv

Bulk flow of halos in ΛCDM simulation

Analysis of the Pangu N-body simulation validates that the bulk flow of halos follows a Maxwellian distribution which variance is consistent with the prediction of the linear theory of structure formation. We propose that the consistency between the observed bulk velocity and theories should be examined at the effective scale of the radius of a spherical top-hat window function yielding the same smoothed velocity variance in linear theory as the sample window function does. We compared some recently estimated bulk flows from observational samples with the prediction of the ΛCDM model we used; some results deviate from expectation at a level of ~ 3σbut the discrepancy is not as severe as previously claimed. We show that bulk flow is only weakly correlated with the dipole of the internal mass distribution, the alignment angle between the mass dipole and the bulk flow has a broad distribution peaked at ~ 30-50 deg., and also that the bulk flow shows little dependence on the mass of the halos used in the estimation. In a simulation of box size 1Gpc/h, for a cell of radius 100 Mpc/h the maximal bulk velocity is >500 km/s, dipoles of the environmental mass outside the cell are not tightly aligned with the bulk flow, but are rather located randomly around it with separation angles ~ 20-40 deg. In the fastest cell there is a slightly smaller number of low-mass halos; however halos inside are clustered more strongly at scales > ~ 20 Mpc/h, which might be a significant feature since the correlation between bulk flow and halo clustering actually increases in significance beyond such scales.

preprint2012arXiv

Influence of baryons on spatial distribution of matter: higher order correlation functions

Baryonic physical processes could leave non-negligible imprint on cosmic matter distribution pattern. Series of high precision simulation data sets with identical initial condition are employed for count-in-cell (CIC) analysis, including one N-body dark matter run, one with adiabatic gas only and one with dissipative processes. Variances and higher order correlation functions of dark matter and gas are estimated. It is found that baryon physical processes mainly affected dark matter distribution at scales less than $1h^{-1}$Mpc. In comparison with the pure dark matter run, adiabatic process alone strengthens variance of dark matter by \sim 10% at scale $0.1h^{-1}$Mpc, while $S_n$s of dark matter deviate from pure dark matter case only mildly at a few percentages. Dissipative gas run does not differ much to the adiabatic run in dark matter variance, but renders significantly different $S_n$ parameters of dark matter, bringing about more than 10% enhancement to $S_3$ at $0.1h^{-1}$Mpc and $z=0$. Distribution patterns of gas in two hydrodynamical simulations are prominently different. Variance of gas at $z=0$ decreases by $\sim 30%$ in adiabatic simulation while by $\sim 60%$ in non-adiabatic simulation at $0.1h^{-1}$Mpc, the attenuation is weaker at larger scales but still obvious at $\sim 10h^{-1}$Mpc. $S_n$ parameters of gas are biased upward at scales $< \sim 4h^{-1}$Mpc, dissipative processes give $\sim 84%$ promotion at $z=0$ to $S_3$ at $0.1h^{-1}$Mpc against the moderate $\sim 7%$ in adiabatic simulation. The clustering segregation we observed between gas and dark matter could have intricate implication on modeling galaxy distribution and relevant cosmological application demanding fine details of matter distribution in strongly nonlinear regime.

preprint2011arXiv

On the fairness of the main galaxy sample of SDSS

Flux-limited and volume-limited galaxy samples are constructed from SDSS data releases DR4, DR6 and DR7 for statistical analysis. The two-point correlation functions $ξ(s)$, monopole of three-point correlation functions $ζ_0$, projected two-point correlation function $w_p$ and pairwise velocity dispersion $σ_{12}$ are measured to test if galaxy samples are fair for these statistics. We find that with increment of sky coverage of SDSS, $ξ(s)$ of flux-limited sample is extremely robust and insensitive to local structures at low redshift. But for volume-limited samples fainter than $L^*$ at large scales $s>\sim 10\hmpc$, deviation of $ξ(s)$ and $ζ_0$ of DR7 to those of DR4 and DR6 increases with larger absolute magnitude. In the weakly nonlinear regime, there is no agreement between $ζ_0$ of different data releases in all luminosity bins. Furthermore, $w_p$ of volume-limited samples of DR7 in luminosity bins fainter than $-M_{r,0.1}=[18.5,19.5]$ are significantly larger, and $σ_{12}$ of the two faintest volume-limited samples of DR7 display very different scale dependence than results of DR4 and DR6. Our findings call for cautions in understanding clustering analysis results of SDSS faint galaxy samples, and higher order statistics of SDSS volume-limited samples in the weakly nonlinear regime. The first zero-crossing points of $ξ(s)$ of volume-limited samples are also investigated and discussed.

preprint2011arXiv

The kinetic SZ tomography with spectroscopic redshift surveys

The kinetic Sunyaev Zel'dovich effect (kSZ) effect is a potentially powerful probe to the missing baryons. However, the kSZ signal is overwhelmed by various contaminations and the cosmological application is hampered by loss of redshift information due to the projection effect. We propose a kSZ tomography method to alleviate these problems, with the aid of galaxy spectroscopic redshift surveys. We propose to estimate the large scale peculiar velocity through the 3D galaxy distribution, weigh it by the 3D galaxy density and adopt the product projected along the line of sight with a proper weighting as an estimator of the true kSZ temperature fluctuation $Θ$. We thus propose to measure the kSZ signal through the $\HatΘ$-$Θ$ cross correlation. This approach has a number of advantages (see details in the abstract of the paper). We test the proposed kSZ tomography against non-adiabatic and adiabatic hydrodynamical simulations. We confirm that $\hatΘ$ is indeed tightly correlated with $Θ$ at $k\la 1h/$Mpc, although nonlinearities in the density and velocity fields and nonlinear redshift distortion do weaken the tightness of the $\hatΘ$-$Θ$ correlation. We further quantify the reconstruction noise in $\HatΘ$ from galaxy distribution shot noise. Based on these results, we quantify the applicability of the proposed kSZ tomography for future surveys. We find that, in combination with the BigBOSS-N spectroscopic redshift survey, the PLANCK CMB experiment will be able to detect the kSZ with an overall significance of $\sim 50σ$ and further measure its redshift distribution at many redshift bins over $0<z<2$.

preprint2010arXiv

A redshift distortion free correlation function at third order in the nonlinear regime

The zeroth-order component of the cosine expansion of the projected three-point correlation function is proposed for clustering analysis of cosmic large scale structure. These functions are third order statistics but can be measured similarly to the projected two-point correlations. Numerical experiments with N-body simulations indicate that the advocated statistics are redshift distortion free within 10% in the non-linear regime on scales ~0.2-10Mpc/h. Halo model prediction of the zeroth-order component of the projected three-point correlation function agrees with simulations within ~10%. This lays the ground work for using these functions to perform joint analyses with the projected two-point correlation functions, exploring galaxy clustering properties in the framework of the halo model and relevant extensions.

preprint2010arXiv

The ergodicity bias in the observed galaxy distribution

The spatial distribution of galaxies we observed is subject to the given condition that we, human beings are sitting right in a galaxy -- the Milky Way. Thus the ergodicity assumption is questionable in interpretation of the observed galaxy distribution. The resultant difference between observed statistics (volume average) and the true cosmic value (ensemble average) is termed as the ergodicity bias. We perform explicit numerical investigation of the effect for a set of galaxy survey depths and near-end distance cuts. It is found that the ergodicity bias in observed two- and three-point correlation functions in most cases is insignificant for modern analysis of samples from galaxy surveys and thus close a loophole in precision cosmology. However, it may become non-negligible in certain circumstances, such as those applications involving three-point correlation function at large scales of local galaxy samples. Thus one is reminded to take extra care in galaxy sample construction and interpretation of the statistics of the sample, especially when the characteristic redshift is low.

preprint2009arXiv

Clustering of K-band selected local galaxies

We present detailed clustering analysis of a large K-band selected local galaxy sample, which is constructed from the 2MASS and the SDSS and consists of $82,486$ galaxies with $10 < K < 13.5$ and $0.01 < z < 0.1$. The two-point correlation function of the magnitude-limited sample in real space at small scales is well described by a power law $ξ(r)=(r/6.44\pm0.23)^{-1.81\pm0.02}$. The pairwise velocity dispersion is derived from the anisotropic two-point correlation function and find the dispersion $σ_{12}=685\pm 17\kms$ if its scale invariance is assumed, which is larger than values measured in optical bands selected galaxy samples. We further investigate the dependence of the two-point correlation function and the $σ_{12}$ on the $g-r$ color and the $K$-band luminosity, obtain similar results to previous works in optical bands. Comparing a mock galaxy sample with our real data indicates that the semi-analytical model can not mimic the $σ_{12}$ in observation albeit it can approximate the two-point correlation function within measurement uncertainties.

preprint2006arXiv

Note on Redshift Distortion in Fourier Space

We explore features of redshift distortion in Fourier analysis of N-body simulations. The phases of the Fourier modes of the dark matter density fluctuation are generally shifted by the peculiar motion along the line of sight, the induced phase shift is stochastic and has probability distribution function (PDF) symmetric to the peak at zero shift while the exact shape depends on the wave vector, except on very large scales where phases are invariant by linear perturbation theory. Analysis of the phase shifts motivates our phenomenological models for the bispectrum in redshift space. Comparison with simulations shows that our toy models are very successful in modeling bispectrum of equilateral and isosceles triangles at large scales. In the second part we compare the monopole of the power spectrum and bispectrum in the radial and plane-parallel distortion to test the plane-parallel approximation. We confirm the results of Scoccimarro (2000) that difference of power spectrum is at the level of 10%, in the reduced bispectrum such difference is as small as a few percents. However, on the plane perpendicular to the line of sight of k_z=0, the difference in power spectrum between the radial and plane-parallel approximation can be more than 10%, and even worse on very small scales. Such difference is prominent for bispectrum, especially for those configurations of tilted triangles. The non-Gaussian signals under radial distortion on small scales are systematically biased downside than that in plane-parallel approximation, while amplitudes of differences depend on the opening angle of the sample to the observer. The observation gives warning to the practice of using the power spectrum and bispectrum measured on the k_z=0 plane as estimation of the real space statistics.

preprint2004arXiv

Three-point Correlation Functions of SDSS Galaxies in Redshift Space: Morphology, Color, and Luminosity Dependence

We present measurements of the redshift--space three-point correlation function of galaxies in the Sloan Digital Sky Survey (SDSS). For the first time, we analyze the dependence of this statistic on galaxy morphology, color and luminosity. In order to control systematics due to selection effects, we used $r$--band, volume-limited samples of galaxies, constructed from the magnitude-limited SDSS data ($14.5<r<17.5$), and further divided the samples into two morphological types (early and late) or two color populations (red and blue). The three-point correlation function of SDSS galaxies follow the hierarchical relation well and the reduced three-point amplitudes in redshift--space are almost scale-independent ($Q_z=0.5\sim1.0$). In addition, their dependence on the morphology, color and luminosity is not statistically significant. Given the robust morphological, color and luminosity dependences of the two-point correlation function, this implies that galaxy biasing is complex on weakly non-linear to non-linear scales. We show that simple deterministic linear relation with the underlying mass could not explain our measurements on these scales.

preprint1995arXiv

Chemistry of Nanoscale Semiconductor Clusters

The ground state structures of small silicon clusters are determined through exhaustive tight-binding molecular dynamics simulation studies. These simulations revealed that \Si{11} is an icosahedron with one missing cap, \Si{12} is a complete icosahedron, \Si{13} is a surface capped icosahedron, \Si{14} is a 4-4-4 layer structure with two caps, \Si{15} is a 1-5-3-5-1 layer structure, and \Si{16} is a partially closed cage consisting of five-membered rings. The characteristic feature of these clusters is that they are all surface. Smalley and co-workers discovered that chemisorption reactivities of silicon clusters vary over three orders of magnitude as a function of cluster size. In particular, they found that \Si{33}, \Si{39}, and \Si{45} clusters are least reactive towards various reagents compared to their immediate neighbors in size. We provide insights into this observed reactivity pattern through our stuffed fullerene model. This structural model consists of bulk-like core of five atoms surrounded by fullerene-like surface. Reconstruction of the ideal fullerene geometry gives rise to four-fold coordinated crown atoms and $π$-bonded dimer pairs. This model yields unique structures for \Si{33}, \Si{39}, and \Si{45} clusters without any dangling bonds and thus explains their lowest reactivity towards chemisorption of closed shell reagents. We also explain why a) these clusters are substantially unreactive compared to bulk surfaces and b) dissociative chemisorption occurs on bulk surfaces while molecular chemisorption occurs on cluster surfaces. Finally, experiments on Si$_x$X$_y$ (X = B, Al, Ga, P, As, AlP, GaAs) are suggested as a means of verifying the proposed model.