The H$\alpha$ and H$\beta$ emission line luminosities measured in a single integrated spectrum are affected in non-trivial ways by point-to-point variations in dust attenuation in a galaxy. This work investigates the impact of this variation when estimating global H$\alpha$ luminosities corrected for the presence of dust by a global Balmer decrement. Analytical arguments show that the dust-corrected H$\alpha$ luminosity is always underestimated when using the global H$\alpha$/H$\beta$ flux ratio to correct for dust attenuation. We measure this effect on 156 face-on star-forming galaxies from the Mapping Nearby Galaxies at APO (MaNGA) survey. At 1-2 kpc spatial resolution, the effect is small but systematic, with the integrated dust-corrected H$\alpha$ luminosity underestimated by $2$-$4$ per cent (and typically not more than by $10$ per cent), and depends on the specific star formation rate of the galaxy. Given the spatial resolution of MaNGA, these are lower limits for the effect. From Multi Unit Spectroscopic Explorer (MUSE) observations of NGC 628 with a resolution of 36 pc we find the discrepancy between the globally and the point-by-point dust-corrected H$\alpha$ luminosity to be $14 \pm 1$ per cent, which may still underestimate the true effect. We use toy models and simulations to show that the true difference depends strongly on the spatial variance of the H$\alpha$/H$\beta$ flux ratio, and on the slope of the relation between H$\alpha$ luminosity and dust attenuation within a galaxy. Larger samples of higher spatial resolution observations are required to quantify the dependence of this effect as a function of galaxy properties.

We used the mark weighted correlation functions (MCFs), $W(s)$, to study the large scale structure of the Universe. We studied five types of MCFs with the weighting scheme $\rho^\alpha$, where $\rho$ is the local density, and $\alpha$ is taken as $-1,\ -0.5,\ 0,\ 0.5$, and 1. We found that different MCFs have very different amplitudes and scale-dependence. Some of the MCFs exhibit distinctive peaks and valleys that do not exist in the standard correlation functions. Their locations are robust against the redshifts and the background geometry, however it is unlikely that they can be used as ``standard rulers'' to probe the cosmic expansion history. Nonetheless we find that these features may be used to probe parameters related with the structure formation history, such as the values of $\sigma_8$ and the galaxy bias. Finally, after conducting a comprehensive analysis using the full shapes of the $W(s)$s and $W_{\Delta s}(\mu)$s, we found that, combining different types of MCFs can significantly improve the cosmological parameter constraints. Compared with using only the standard correlation function, the combinations of MCFs with $\alpha=0,\ 0.5,\ 1$ and $\alpha=0,\ -1,\ -0.5,\ 0.5,\ 1$ can improve the constraints on $\Omega_m$ and $w$ by $\approx30\%$ and $50\%$, respectively. We find highly significant evidence that MCFs can improve cosmological parameter constraints.

We developed a 3-D elastic Boundary Element Method (BEM) computer code, called AstroSeis, to model seismic wavefields in a body with an arbitrary shape, such as an asteroid. Besides the AstroSeis can handle arbitrary surface topography, it can deal with a liquid core in an asteroid model. Both the solid and liquid domains are homogenous in our current code. For seismic sources, we can use single forces or moment tensors. The AstroSeis is implemented in the frequency domain and the frequency-dependent Q can be readily incorporated. The code is in MATLAB and it is straightforward to set up the model to run our code. The frequency-domain calculation is advantageous to study the long-term elastic response of a celestial body due to a cyclic force such as the tidal force with no numerical dispersion issue suffered by many other methods requiring volume meshing. Our AstroSeis has been benchmarked with other methods such as normal modes summation and the direct solution method (DSM). This open-source AstroSeis will be a useful tool to study the interior and surface processes of asteroids.

The interstellar oxygen isotopic ratio of 18O/17O can reflect the relative amount of the secular enrichment by ejecta from high-mass versus intermediate-mass stars. Previous observations found a Galactic gradient of 18O/17O, i.e., low ratios in the Galactic center and large values in the Galactic disk, which supports the insideout formation scenario of our Galaxy. However, the observed objects are not many and, in particular, not so many at large galactocentric distances. Thus we started a systematic study on Galactic interstellar 18O/17O, through observations of C18O and C17O multi-transition lines toward a large sample of 286 sources (at least one order of magnitude larger than previous ones), from the Galactic center region to the far outer Galaxy (~22 kpc). Here we present our observations of J=1-0 lines of C18O and C17O, with the ARO12m and the IRAM 30m telescope. We detected successfully both C18O and C17O 1-0 lines for 34 sources among our IRAM30m sample of 50 targets and for 166 sources among our ARO12m sample of 260 targets. The C18O optical depth effect on our ratio results, evaluated by fitting results of C17O spectra with hyperfine components and our RADEX non-LTE model calculation for the strongest source, was found to be insignificant. Beam dilution does not seem to be a problem either, which was supported by the fact of no systematic variation between the isotopic ratio and the heliocentric distance, and consistent measured ratios from two telescopes for most of those detected sources. Our results, though there are still very few detections made for sources in the outer Galaxy, confirm the apparent 18O/17O gradient of 18O/17O = (0.10+-0.03)R_GC+(2.95+-0.30), with a Pearson's rank correlation coefficient R = 0.69. This is supported by the newest Galactic chemical evolution model including the impact of massive stellar rotators and novae.

The circumgalactic medium (CGM) of the Milky Way is mostly obscured by nearby gas in position-velocity space because we reside inside the Galaxy. Substantial biases exist in most studies on the Milky Way's CGM that focus on easier-to-detect high-velocity gas. With mock observations on a Milky-Way analog from the FOGGIE simulation, we investigate four observational biases related to the Milky Way's CGM. First, QSO absorption-line studies probe a limited amount of the CGM mass: only 35% of the mass is at high Galactic latitudes $|b|>20$ degrees, of which only half is moving at $|v_{\rm LSR}|\gtrsim100$ km s$^{-1}$. Second, the inflow rate ($\dot{M}$) of the cold gas observable in HI 21cm is reduced by a factor of $\sim10$ as we switch from the local standard of rest to the galaxy's rest frame; meanwhile $\dot{M}$ of the cool and warm gas does not change significantly. Third, OVI and NV are promising ions to probe the Milky Way's outer CGM ($r\gtrsim$15 kpc), but CIV may be less sensitive. Lastly, the scatter in ion column density is a factor of 2 higher if the CGM is observed from inside-out than from external views because of the gas radial density profile. Our work highlights that observations of the Milky Way's CGM, especially those using HI 21cm and QSO absorption lines, are highly biased. We demonstrate that these biases can be quantified and calibrated through synthetic observations with simulated Milky-Way analogs.

Recent integral field spectroscopic (IFS) surveys have revealed radial gradients in the optical spectral indices of post-starburst galaxies, which can be used to constrain their formation histories. We study the spectral indices of post-processed mock IFS datacubes of binary merger simulations, carefully matched to the properties of the MaNGA IFS survey, with a variety of black hole feedback models, progenitor galaxies, orbits and mass ratios. Based on our simulation sample, we find that only major mergers on prograde-prograde or retrograde-prograde orbits in combination with a mechanical black hole feedback model can form galaxies with weak enough ongoing star formation, and therefore absent H$\alpha$ emission, to be selected by traditional PSB selection methods. We find strong fluctuations in nebular emission line strengths, even within the PSB phase, suggesting that H$\alpha$ selected PSBs are only a subsample of the underlying population. The global PSB population can be more robustly identified using stellar continuum-based approaches. The difficulty in reproducing the very young PSBs in simulations potentially indicates that new sub-resolution star formation recipes are required to properly model the process of star formation quenching. In our simulations, we find that the starburst peaks at the same time at all radii, but is stronger and more prolonged in the inner regions. This results in a strong time evolution in the radial gradients of the spectral indices which can be used to estimate the age of the starburst without reliance on detailed star formation histories from spectral synthesis models.

The redshift-space distortion (RSD) in the observed distribution of galaxies is known as a powerful probe of cosmology. Observations of large-scale RSD have given tight constraints on the linear growth rate of the large-scale structures in the universe. On the other hand, the small-scale RSD, caused by galaxy random motions inside clusters, has not been much used in cosmology, but also has cosmological information because universes with different cosmological parameters have different halo mass functions and virialized velocities. We focus on the projected correlation function $w(r_p)$ and the multipole moments $\xi_l$ on small scales ($1.4$ to $30\ h^{-1}\rm{Mpc}$). Using simulated galaxy samples generated from a physically motivated most bound particle (MBP)-galaxy correspondence scheme in the Multiverse Simulation, we examine the dependence of the small-scale RSD on the cosmological matter density parameter $\Omega_m$, the satellite velocity bias with respect to MBPs, $b_v^s$, and the merger-time-scale parameter $\alpha$. We find that $\alpha=1.5$ gives an excellent fit to the $w(r_p)$ and $\xi_l$ measured from the SDSS-KIAS value added galaxy catalog. We also define the ``strength'' of Fingers-of-God as the ratio of the parallel and perpendicular size of the contour in the two-point correlation function set by a specific threshold value and show that the strength parameter helps constraining $(\Omega_m, b_v^s, \alpha)$ by breaking the degeneracy among them. The resulting parameter values from all measurements are $(\Omega_m,b_v^s)=(0.272\pm0.013,0.982\pm0.040)$, indicating a slight reduction of satellite galaxy velocity relative to the MBP. However, considering that the average MBP speed inside haloes is $0.94$ times the dark matter velocity dispersion, the main drivers behind the galaxy velocity bias are gravitational interactions, rather than baryonic effects.

Fast radio bursts (FRBs) are short pulses observed in radio band from cosmological distances, some of which emit repeating bursts. The physical origins of these mysterious events have been subject to wide speculations and heated debates. One class of models invoke soft gamma-ray repeaters (SGRs), or magnetars, as the sources of FRBs. Magnetars are rotating neutron stars with extremely strong magnetic field and can sporadically emit bursts from X-ray ($\sim$keV) to soft gamma-ray ($\sim$sub-MeV) with duration from $10^{-2}$~s to $10^2$~s. However, even though some bright radio bursts have been observed from some magnetars, no FRB-like events had been detected to be associated with any magnetar burst, including one giant flare, and no radio burst has been associated with any X-ray event from any magnetar. Therefore, there is still no observational evidence for magnetar-FRB association up to today. Recently, a pair of FRB-like bursts (FRB~200428 hereafter) separated by 30 milliseconds (ms) were detected from the general direction of the Galactic magnetar SGR~J1935+2154. Here we report the detection of a non-thermal X-ray burst in the 1-250 keV energy band with the Insight-HXMT satellite, which we identify as emitted from SGR~J1935+2154. The burst showed two hard peaks with a separation of $\sim30$ ms, consistent with the separation between the two bursts in FRB~200428. The delay time between the double radio and X-ray peaks is $\sim8.57$ s, fully consistent with the dispersion delay of FRB~200428. We thus identify the non-thermal X-ray burst is associated with FRB~200428 whose high energy counterpart is the two hard peaks in X-ray. Our results suggest that the non-thermal X-ray burst and FRB~200428 share the same physical origin in an explosive event from SGR~J1935+2154.

We study ram pressure stripping in simulated Milky Way-like halos at z>=2 from the Figuring Out Gas & Galaxies In Enzo (FOGGIE) project. These simulations reach exquisite resolution in their circumgalactic medium (CGM) gas owing to FOGGIE's novel refinement scheme. The CGM of each halo spans a wide dynamic range in density and velocity over its volume---roughly 6 dex and 1000 km/s, respectively---translating into a 5 dex range in ram pressure imparted to interacting satellites. The ram pressure profiles of the simulated CGM are highly stochastic, owing to kpc-scale variations of the density and velocity fields of the CGM gas. As a result, the efficacy of ram pressure stripping depends strongly on the specific path a satellite takes through the CGM. The ram-pressure history of a single satellite is generally unpredictable and not well correlated with its approach vector with respect to the host galaxy. The cumulative impact of ram pressure on the simulated satellites is dominated by only a few short strong impulses---on average, 90% of the total surface momentum gained through ram pressure is imparted in 20% or less of the total orbital time. These results reveal an erratic mode of ram pressure stripping in Milky-Way like halos at high redshift---one that is not captured by a smooth spherically-averaged model of the circumgalactic medium.

We use the In data collected during the 2019 outburst from X-ray pulsar 4U 1901+03 to complement the orbital parameters reported by Fermi/GBM. Using the Insight-HXMT, we examine the correlation between the derivative of the intrinsic spin frequency and bolometric flux based on accretion torque models. It was found that the pulse profiles significantly evolve during the outburst. The existence of two types of the profile's pattern discovered in the Insight-HXMT data indicates that this source experienced transition between a super-critical and a sub-critical accretion regime during its 2019 outburst. Based on the evolution of the pulse profiles and the torque model, we derive the distance to 4U 1901+03 as 12.4+-0.2 kpc.

We present a new class of non-singular bounce cosmology free from instabilities, using a single scalar field coupled to gravity within the framework of the Degenerate Higher-Order Scalar-Tensor (DHOST) theories. In this type of scenarios, the gradient instability that widely exists in non-singular bounce cosmologies in the framework of scalar-tensor and Horndeski/Galileon theories is removed by the effects of new operators introduced by the DHOST, due to the modification that they later bring about to the dispersion relation of perturbations. Hence, our results demonstrate that there is indeed a loophole for this type of bounce scenarios to be free from pathologies when primordial perturbations evolve through the bounce phase, and thus the theoretical {\it no-go} theorem for non-singular bounce cosmology of Horndeski/Galileon theories can be delicately evaded in DHOST extensions.

Taking a large amount of images, the Cassini Imaging Science Subsystem (ISS) has been routinely used in astrometry. In ISS images, disk-resolved objects often lead to false detection of stars that disturb the camera pointing correction. The aim of this study was to develop an automated processing method to remove the false image stars in disk-resolved objects in ISS images. The method included the following steps: extracting edges, segmenting boundary arcs, fitting circles and excluding false image stars. The proposed method was tested using 200 ISS images. Preliminary experimental results show that it can remove the false image stars in more than 95% of ISS images with disk-resolved objects in a fully automatic manner, i.e. outperforming the traditional circle detection based on Circular Hough Transform (CHT) by 17%. In addition, its speed is more than twice as fast as that of the CHT method. It is also more robust (no manual parameter tuning is needed) when compared with CHT. The proposed method was also applied to a set of ISS images of Rhea to eliminate the mismatch in pointing correction in automatic procedure. Experiment results showed that the precision of final astrometry results can be improve by roughly 2 times than that of automatic procedure without the method. It proved that the proposed method is helpful in the astrometry of ISS images in fully automatic manner.

We prove from the constrained modified gravity (MG) galaxy mock catalogs that the shape of the pairwise kinematic Sunyaev-Zeldovich (kSZ) power spectrum $P_{\rm kSZ}$ has a strong constraining power on discriminating different gravity theories on cosmological scales. By varying the effective optical depth $\tau_{\rm T}$ as a free parameter, we verify that the $\tau_{\rm T}$-$f$ (the linear growth rate) degeneracy in the linear theory of $P_{\rm kSZ}$ is broken down by the non-linear structure growth and the scale-dependence of $f$ in some MG theories. Equivalently speaking, the shape of $P_{\rm kSZ}$ alone could be used to tightly constrain the MG theories on cosmological scales. With a good knowledge of galaxy density biases, we verify that, a combination of the next generation galaxy spectroscopic redshift and CMB surveys, e.g. DESI+CMB-S4, could potentially discriminate $f(R)$ and nDGP models from the general relativity at $\sim5\sigma$ level using the shape of the galaxy pairwise kSZ dipole $P_{{\rm kSZ},\ell=1}$ alone, when $f_{R0}=10^{-5}$ and $H_0r_c=1.0$.

A new black hole X-ray binary (BHXRB) MAXI J1535-571 was discovered by MAXI during its outburst in 2017. Using observations taken by the first Chinese X-ray satellite, the Hard X-ray Modulation Telescope (dubbed as Insight-HXMT), we perform a joint spectral analysis (2-150 keV) in both energy and time domains. The energy spectra provide the essential input for probing the intrinsic Quasi-Periodic Oscillation (QPO) fractional rms spectra (FRS). Our results show that during the intermediate state, the energy spectra are in general consistent with those reported by Swift/XRT and NuSTAR. However, the QPO FRS become harder and the FRS residuals may suggest the presence of either an additional power-law component in the energy spectrum or a turn-over in the intrinsic QPO FRS at high energies.

In a new release of the fourth \emph{Fermi} Large Area Telescope source catalog (4FGL) in September (URL: https://fermi.gsfc.nasa.gov/ssc/data/access/lat/8yr_catalog/), 5065 $\gamma$-ray sources are reported, including 3207 active galactic nuclei (AGNs), 239 pulsars, 190 other sources and 1429 unassociated sources. We employ two different supervised machine learning classifiers, combined with the direct observation parameters given by the 4FGL fits (URL: https://fermi.gsfc.nasa.gov/ssc/data/access/lat/8yr_catalog/gll_psc_v20.fit), to search for sources potentially classified as AGNs and pulsars in the 1429 unassociated sources. In order to reduce the error caused by the large difference in the sizes of samples, we divide the classification process into two separate steps in order to identify the AGNs and the pulsars. First, we select the identified AGNs from all of the samples, and then select the identified pulsars from the remaining. Using the 4FGL sources associated or identified as AGNs, pulsars, and other sources with the features selected through the K-S test, we trained, optimized, and tested our classifier models. Then, the models are applied to classify the 1429 unassociated sources. According to the direct calculation results of the two classifiers, we show the sensitivity, specificity, accuracy in each step, and the class of unassociated sources given by each classifier. The accuracy obtained in the first step is approximately $95\%$; in the second step, the obtained overall accuracy is approximately $80\%$. Combining the results of the two classifiers, we predict that there are 674 AGN-type candidates, 86 pulsar-type candidates, 177 other types of $\gamma$-ray candidates, and 492 of uncertain type.

Compact group (CG) is a kind of special galaxy system where the galaxy members are separated at the distances of the order of galaxy size. The strong interaction between the galaxy members makes CGs ideal labs for studying the environmental effects on galaxy evolution. Traditional photometric selection algorithm biases against the CG candidates at low redshifts, while spectroscopic identification technique is affected by the spectroscopic incompleteness of sample galaxies and typically biases against the high redshift candidates. In this study, we combine these two methods and select CGs in the main galaxy sample of the Sloan Digital Sky Survey, where we also have taken the advantages of the complementary redshift measurements from the LAMOST spectral and GAMA surveys. We have obtained the largest and most complete CG samples to date. Our samples include 6,144 CGs and 8,022 CG candidates, which are unique in the studies of the nature of the CGs and the evolution of the galaxies inside.

We address the spatial scale, ionization structure, mass and metal content of gas at the Milky Way disk-halo interface detected as absorption in the foreground of seven closely-spaced, high-latitude halo blue horizontal branch stars (BHBs) with heights z = 3 - 14 kpc. We detect transitions that trace multiple ionization states (e.g. CaII, FeII, SiIV, CIV) with column densities that remain constant with height from the disk, indicating that the gas most likely lies within z < 3.4 kpc. The intermediate ionization state gas traced by CIV and SiIV is strongly correlated over the full range of transverse separations probed by our sightlines, indicating large, coherent structures greater than 1 kpc in size. The low ionization state material traced by CaII and FeII does not exhibit a correlation with either N$_{\rm HI}$ or transverse separation, implying cloudlets or clumpiness on scales less than 10 pc. We find that the observed ratio log(N_SiIV/ N_CIV), with a median value of -0.69+/-0.04, is sensitive to the total carbon content of the ionized gas under the assumption of either photoionization or collisional ionization. The only self-consistent solution for photoionized gas requires that Si be depleted onto dust by 0.35 dex relative to the solar Si/C ratio, similar to the level of Si depletion in DLAs and in the Milky Way ISM. The allowed range of values for the areal mass infall rate of warm, ionized gas at the disk-halo interface is 0.0003 < dM_gas / dtdA [M_sun kpc^-2 yr^-] < 0.006. Our data support a physical scenario in which the Milky Way is fed by complex, multiphase processes at its disk-halo interface that involve kpc-scale ionized envelopes or streams containing pc-scale, cool clumps.

Supermassive black holes (SMBHs) are thought to provide energy that prevents catastrophic cooling in the centers of massive galaxies and galaxy clusters. However, it remains unclear how this "feedback" process operates. We use high-resolution optical data to study the kinematics of multi-phase filamentary structures by measuring the velocity structure function (VSF) of the filaments over a wide range of scales in the centers of three nearby galaxy clusters: Perseus, Abell 2597 and Virgo. We find that the motions of the filaments are turbulent in all three clusters studied. There is a clear correlation between features of the VSFs and the sizes of bubbles inflated by SMBH driven jets. Our study demonstrates that SMBHs are the main driver of turbulent gas motions in the centers of galaxy clusters and suggests that this turbulence is an important channel for coupling feedback to the environment. Our measured amplitude of turbulence is in good agreement with Hitomi Doppler line broadening measurement and X-ray surface brightness fluctuation analysis, suggesting that the motion of the cold filaments is well-coupled to that of the hot gas. The smallest scales we probe are comparable to the mean free path in the intracluster medium (ICM). Our direct detection of turbulence on these scales provides the clearest evidence to date that isotropic viscosity is suppressed in the weakly-collisional, magnetized intracluster plasma.

The second observational campaign of gravitational waves organized by the LIGO/Virgo Collaborations has led to several breakthroughs such as the detection of gravitational wave signals from merger systems involving black holes or neutrons stars. During O2,14 gravitational wave alerts were sent to the astronomical community with sky regions covering mostly over hundreds of square degrees. Among them, 6 have been finally confirmed as real astrophysical events. Since 2013, a new set of ground-based robotic telescopes called GWAC and its pathfinder mini-GWAC have been developed to contribute to the various challenges of themulti-messenger and time domain astronomy. The GWAC system is built up in the framework of the ground-segment system of the SVOM mission that will be devoted to the study of the multi-wavelength transient sky in the next decade. During O2, only the mini-GWAC telescopenetwork was fully operational. Due to the wide field of view and fast automatic follow-up capabilities of the mini-GWAC telescopes, they were well adapted to efficiently cover the sky localization areas of the gravitational wave event candidates. In this paper, we present the mini-GWAC pipeline we have set up to respond to the GW alerts and we report our optical follow-up observations of 8 GW alerts detected during the O2 run. Our observations provided the largest coverage of the GW localization areas in a short latency made by any optical facility. We found tens of optical transient candidates in our images, but none of those could be securely associated with any confirmed black hole-black hole merger event. Based on this first experience and the near future technical improvements of our network system, we will be more competitive to detect the optical counterparts from some gravitational wave events that will be detected during the upcoming O3 run, especially those emerging from binary neutron star mergers.

We report a tentative detection of the circumgalactic medium (CGM) of WLM, an isolated, low-mass (log$M_*/M_\odot\approx7.6$), dwarf irregular galaxy in the Local Group (LG). We analyze an HST/COS archival spectrum of a quasar sightline (PHL2525), which is 45 kpc (0.5 virial radius) from WLM and close to the Magellanic Stream (MS). Along this sightline, two ion absorbers are detected in Si II, Si III, Si IV, C II, and C IV at velocities of $\sim$-220 km s$^{-1}$ (Component v-220) and $\sim$-150 km s$^{-1}$ (Component v-150). To identify their origins, we study the position-velocity alignment of the components with WLM and the nearby MS. Near the Magellanic longitude of PHL2525, the MS-related neutral and ionized gas moves at $\lesssim-190$ km s$^{-1}$, suggesting an MS origin for Component v-220, but not for Component v-150. Because PHL2525 passes near WLM and Component v-150 is close to WLM's systemic velocity ($\sim$-132 km s$^{-1}$), it is likely that Component v-150 arises from the galaxy's CGM. This results in a total Si mass in WLM's CGM of $M_{\rm Si}^{\rm CGM}\sim(0.2-1.0)\times10^5~M_\odot$ using assumption from other COS dwarf studies. Comparing $M_{\rm Si}^{\rm CGM}$ to the total Si mass synthesized in WLM over its lifetime ($\sim$1.3$\times10^5~M_\odot$), we find $\sim$3% is locked in stars, $\sim$6% in the ISM, $\sim$15%-77% in the CGM, and the rest ($\sim$14%-76%) is likely lost beyond the virial radius. Our finding resonates with other COS dwarf galaxy studies and theoretical predictions that low-mass galaxies can easily lose metals into their CGM due to stellar feedback and shallow gravitational potential.

We present a novel absorption line survey using 54 blue horizontal branch stars (BHBs) in the Milky Way halo as background sources for detecting gas flows at the disk-halo interface. Distance measurements to high-latitude ($b$ > 60{\deg}) background stars at 3.1-13.4 kpc, combined with unprecedented spatial sampling and spectral resolution, allow us to examine the 3-dimensional spatial distribution and kinematics of gas flows near the disk. We detect absorption signatures of extraplanar CaII and NaI in Keck HIRES spectra and find that their column densities exhibit no trend with distance to the background sources, indicating that these clouds lie within 3.1 kpc of the disk. We calculate covering fractions of $f_{CaII}$ = 63%, $f_{NaI}$ = 26%, and $f_{HI}$ = 52%, consistent with a picture of the CGM that includes multi-phase clouds containing small clumps of cool gas within hotter, more diffuse gas. Our measurements constrain the scale of any substructure within these cool clouds to <0.5 kpc. CaII and NaI absorption features exhibit an intermediate-velocity (IV) component inflowing at velocities of -75 km/s < v < -25 km/s relative to the local standard of rest, consistent with previously-studied HI structures in this region. We report the new detection of an inflow velocity gradient ${\Delta}v_z$ ~ 6-9 km/s/kpc across the Galactic plane. These findings place constraints on the physical and kinematic properties of CGM gas flows through the disk-halo interface, and support a galactic fountain model in which cold gas rains back onto the disk.

We develop an improved Alcock-Paczynski (AP) test method that uses the redshift-space two-point correlation function (2pCF) of galaxies. Cosmological constraints can be obtained by examining the redshift dependence of the normalized 2pCF, which should not change apart from the expected small non-linear evolution. An incorrect choice of cosmology used to convert redshift to comoving distance will manifest itself as redshift-dependent 2pCF. Our method decomposes the redshift difference of the two-dimensional correlation function into the Legendre polynomials whose amplitudes are modeled by radial fitting functions. Our likelihood analysis with this 2-D fitting scheme tightens the constraints on $\Omega_m$ and ${w}$ by $\sim 40\%$ compared to the method of Li et al. (2016, 2017, 2018) that uses one dimensional angular dependence only. We also find that the correction for the non-linear evolution in the 2pCF has a non-negligible cosmology dependence, which has been neglected in previous similar studies by Li et al.. With an accurate accounting for the non-linear systematics and use of full two-dimensional shape information of the 2pCF down to scales as small as $5~h^{-1}{\rm Mpc}$ it is expected that the AP test with redshift-space galaxy clustering anisotropy can be a powerful method to constrain the expansion history of the universe.

The mapping of galaxy clustering from real space to redshift space introduces the anisotropic property to the measured galaxy density power spectrum in redshift space, known as the redshift space distortion (RSD) effect. The mapping formula is intrinsically non-linear, which is complicated by the higher order polynomials due to indefinite orders of cross correlations between density and velocity fields, and the Finger--of--God (FoG) effect due to the randomness of the galaxy peculiar velocity field. In previous works, we have verified the robustness of advanced TNS mapping formula in our hybrid RSD model in dark matter case, where the halo bias models are not taken into account for the halo mapping formula in redshift space. Using 100 realizations of halo catalogs in N-body simulations, we find that our halo RSD model with the known halo bias model and the effective FoG function accurately predicts the halo power spectrum measurements, within 1$\sim$2% accuracy up to $k\sim 0.2h$/Mpc, depending on different halo masses and redshifts.

We formulate and solve the problem of spherically symmetric, steady state, adiabatic accretion onto a Schwarzschild-like black hole obtained recently. We derive the general analytic expressions for the critical points, the critical velocity, the critical speed of sound, and subsequently the mass accretion rate. The case for polytropic gas is discussed in detail. We find the parameter characterizing the breaking of Lorentz symmetry will slow down the mass accretion rate, while has no effect on the gas compression and the temperature profile below the critical radius and at the event horizon.

The expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $z < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. The large cosmological volume at $2 < z < 5$, together with the ability to efficiently target high-$z$ galaxies with known techniques, enables large gains in the study of Inflation and Dark Energy. A future spectroscopic survey can test the Gaussianity of the initial conditions up to a factor of ~50 better than our current bounds, crossing the crucial theoretical threshold of $\sigma(f_{NL}^{\rm local})$ of order unity that separates single field and multi-field models. Simultaneously, it can measure the fraction of Dark Energy at the percent level up to $z = 5$, thus serving as an unprecedented test of the standard model and opening up a tremendous discovery space.

The past decade has seen an explosion of discoveries and new insights into the diffuse gas within galaxies, galaxy clusters, and the filaments composing the Cosmic Web. A new decade will bring fresh opportunities to further this progress towards developing a comprehensive view of the composition, thermal state, and physical processes of diffuse gas in the Universe. Ultraviolet (UV) spectroscopy, probing diffuse 10^4-10^6 K gas at high spectral resolution, is uniquely poised to (1) witness environmental galaxy quenching processes in action, such as strangulation and tidal- and ram-pressure stripping, (2) directly account for the baryon content of galaxy clusters in the cold-warm (T<10^6 K) gas, (3) determine the phase structure and kinematics of gas participating in the equilibrium-regulating exchange of energy at the cores of galaxy clusters, and (4) map cold streams and filaments of the Cosmic Web that feed galaxies and clusters. With a substantial UV undertaking beyond the Hubble Space Telescope, all of the above would be achievable over the entire epoch of galaxy cluster formation. Such capabilities, coupled with already-planned advancements at other wavelengths, will transform extragalactic astronomy by revealing the dominant formation and growth mechanisms of gaseous halos over the mass spectrum, settling the debate between early- and late-time metal enrichment scenarios, and revealing how the ecosystems in which galaxies reside ultimately facilitate their demise.

Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with high precision. Nevertheless, a minimal deviation from Gaussianityis perhaps the most robust theoretical prediction of models that explain the observed Universe; itis necessarily present even in the simplest scenarios. In addition, most inflationary models produce far higher levels of non-Gaussianity. Since non-Gaussianity directly probes the dynamics in the early Universe, a detection would present a monumental discovery in cosmology, providing clues about physics at energy scales as high as the GUT scale.

Galaxies evolve under the influence of gas flows between their interstellar medium and their surrounding gaseous halos known as the circumgalactic medium (CGM). The CGM is a major reservoir of galactic baryons and metals, and plays a key role in the long cycles of accretion, feedback, and recycling of gas that drive star formation. In order to fully understand the physical processes at work within galaxies, it is therefore essential to have a firm understanding of the composition, structure, kinematics, thermodynamics, and evolution of the CGM. In this white paper we outline connections between the CGM and galactic star formation histories, internal kinematics, chemical evolution, quenching, satellite evolution, dark matter halo occupation, and the reionization of the larger-scale intergalactic medium in light of the advances that will be made on these topics in the 2020s. We argue that, in the next decade, fundamental progress on all of these major issues depends critically on improved empirical characterization and theoretical understanding of the CGM. In particular, we discuss how future advances in spatially-resolved CGM observations at high spectral resolution, broader characterization of the CGM across galaxy mass and redshift, and expected breakthroughs in cosmological hydrodynamic simulations will help resolve these major problems in galaxy evolution.

The discovery of Fast Radio Bursts (FRBs) at cosmological distances has opened a powerful window on otherwise unseen matter in the Universe. In the 2020s, observations of $>10^{4}$ FRBs will assess the baryon contents and physical conditions in the hot/diffuse circumgalactic, intracluster, and intergalactic medium, and test extant compact-object dark matter models.

The last decade has seen the rapid, concurrent development of new classes of energetic astrophysical transients, including Fast Radio Bursts, Superluminous Supernovae, and ultralong gamma-ray bursts. We consider the implications of an emerging consensus on the origin of these diverse phenomena as different signatures of magnetar birth. We describe how magnetar-powered radio transients from milliseconds to gigaseconds will be used to address diverse physical questions related to the stellar evolution and death, the circum- and intergalactic medium, and fundamental physics on scales ranging from the quantum to cosmological. We argue that high-sensitivity radio interferometer concepts in development will be critical to accessing this science. There is also a growing need for scientists with expertise in new statistical tools and software techniques.

There are still some important unanswered questions about the detailed particle acceleration and escape occurring during the quiescent epoches. As a result, the particle distribution that is adopted in the blazar quiescent spectral model have numerous unconstrained shapes. To help remedy this problem, we introduce a analytical particle transport model to reproduce quiescent broadband spectral energy distribution of blazar. In this model, the exact electron distribution is solved from a generalized transport equation that contains the terms describing first-order and secondary-order \emph{Fermi} acceleration, escape of particle due to both the advection and spatial diffusion, energy losses due to synchrotron emission and inverse-Compton scattering of an assumed soft photon field. We suggest that the advection is a significant escape mechanism in blazar jet. We find that in our model the advection process tends to harden the particle distribution, which enhances the high energy components of resulting synchrotron and synchrotron self-Comptom spectrum from jet. Our model is able to roughly reproduce the observed spectra of extreme BL Lac object 1ES 0414+009 with reasonable assumptions about the physical parameters.

Every quasar (QSO) spectrum contains absorption-line signatures from the interstellar medium, disk-halo interface, and circumgalactic medium (CGM) of the Milky Way (MW). We analyze Hubble Space Telescope/Cosmic Origins Spectrograph (COS) spectra of 132 QSOs to study the significance and origin of SiIV absorption at $|v_{\rm LSR}|\leq100$ km/s in the Galactic halo. The gas in the north predominantly falls in at $-50\lesssim v_{\rm LSR}\lesssim 0$ km/s, whereas in the south, no such pattern is observed. The SiIV column density has an average and a standard deviation of $\langle N_{\rm SiIV}\rangle=(3.8\pm1.4)\times10^{13}$ cm$^{-2}$. At $|b|\gtrsim 30$ degree, $N_{\rm SiIV}$ does not significantly correlate with $b$, which cannot be explained by a commonly adopted flat-slab geometry. We propose a two-component model to reconstruct the $N_{\rm SiIV}$-$b$ distribution: a plane-parallel component $N_{\rm DH}^{\perp}$ to account for the MW's disk-halo interface and a global component $N_{\rm G}$ to reproduce the weak dependence on $b$. We find $N_{\rm DH}^{\perp}=1.3^{+4.7}_{-0.7}\times10^{12}$ cm$^{-2}$ and $N_{\rm G}=(3.4\pm0.3)\times10^{13}$ cm$^{-2}$ on the basis of Bayesian analyses and block bootstrapping. The global component is most likely to have a Galactic origin, although its exact location is uncertain. If it were associated with the MW's CGM, we would find $M_{\rm gas, all}\gtrsim4.7\times10^9\ M_{\odot} (\frac{C_f}{1})(\frac{R}{75\ {\rm kpc}})^2 (\frac{f_{\rm SiIV}}{0.3})^{-1}(\frac{Z}{0.3\ Z_{\odot}})^{-1}$ for the cool gas at all velocities in the Galactic halo. Our analyses show that there is likely a considerable amount of gas at $|v_{\rm LSR}|\leq100$ km s$^{-1}$ hidden in the MW's CGM. Along with this work, we make our QSO dataset publicly available as the COS Quasar Database for Galactic Absorption Lines (COS-GAL).

Both multi-streaming (random motion) and bulk motion cause the Finger-of-God (FoG) effect in redshift space distortion (RSD). We apply a direct measurement of the multi-streaming effect in RSD from simulations, proving that it induces an additional, non-negligible FoG damping to the redshift space density power spectrum. We show that, including the multi-streaming effect, the RSD modelling is significantly improved. We also provide a theoretical explanation based on halo model for the measured effect, including a fitting formula with one to two free parameters. The improved understanding of FoG helps break the $f\sigma_8-\sigma_v$ degeneracy in RSD cosmology, and has the potential of significantly improving cosmological constraints.

The precise localization (<1'') of multiple fast radio bursts (FRBs) to z>0.1 galaxies has confirmed that the dispersion measures (DMs) of these enigmatic sources afford a new opportunity to probe the diffuse ionized gas around and in between galaxies. In this manuscript, we examine the signatures of gas in dark matter halos (aka halo gas) on DM observations in current and forthcoming FRB surveys. Combining constraints from observations of the high velocity clouds, OVII absorption, and the DM to the Large Magellanic Cloud with hydrostatic models of halo gas, we estimate that our Galactic halo will contribute ${\rm DM}_{\rm MW,halo} \approx 50-80 \rm pc/cm^{-3}$ from the Sun to 200 kpc independent of any contribution from the Galactic ISM. Extending analysis to the Local Group, we demonstrate that M31's halo will be easily detected by high-sample FRB surveys (e.g. CHIME) although signatures from a putative Local Group medium may compete. We then review current empirical constraints on halo gas in distant galaxies and discuss the implications for their DM contributions. We further examine the DM probability distribution function of a population of FRBs at z >> 0 using an updated halo mass function and new models for the halo density profile. Lastly, we illustrate the potential of FRB experiments for resolving the baryonic fraction of halos by analyzing simulated sightlines through the CASBaH survey. All of the code and data products of our analysis are available at https://github.com/FRBs.

Tidal force plays an important role in the evolution of the planet-moon system. The tidal force of a moon can excite seismic waves in the planet it is orbiting. A tidal-seismic resonance is expected when a tidal force frequency matches a free-oscillation frequency of the planet. Here we show that when the moon is close to the planet, the tidal-seismic resonance can cause large-amplitude seismic waves, which can change the shape of the planet and in turn exert a negative torque on the moon to cause it to fall rapidly toward the planet. We postulate that the tidal-seismic resonance may be an important mechanism which can accelerate planet accretion process. On the other hand, tidal-seismic resonance effect can also be used to interrogate planet interior by long term tracking of the orbital change of the moon.

We investigate whether the new horizon first law proposed recently still work in $f(R)$ theory. We identify the entropy and the energy of black hole as quantities proportional to the corresponding value of integration, supported by the fact that the new horizon first law holds true as a consequence of equations of motion in $f(R)$ theories. The formulas for the entropy and energy of black hole found here are in agreement with the results obtained in literatures. For applications, some nontrivial black hole solutions in $f(R)$ theories have been considered, the entropies and the energies of black holes in these models are firstly computed, which may be useful for future researches.

In this paper, we compile the very-high-energy and high-energy spectral indices of 43 BL Lac objects from the literature. Based on a simple math model, $\Delta\Gamma_{obs}=\alpha {\rm{z}}+\beta $, we present evidence for the origin of an observed spectral break that is denoted by the difference between the observed very-high-energy and high-energy spectral indices, $\Delta\Gamma_{obs}$. We find by linear regression analysis that $\alpha\ne 0$ and $\beta\ne 0$. These results suggest that the extragalactic background light attenuation and the intrinsic curvature dominate on the GeV-TeV $\gamma$-ray energy spectral break of BL Lac objects. We argue that the extragalactic background light attenuation is an exclusive explanation for the redshift evolution of the observed spectral break.

A long-standing issue in peculiar velocity cosmology is whether the halo/galaxy velocity bias $b_v=1$ at large scale. The resolution of this important issue must resort to high precision cosmological simulations. However, this is hampered by another long-standing `sampling artifact' problem in volume weighted velocity measurement. We circumvent this problem with a hybrid approach. We first measure statistics free of sampling artifact, then link them to volume weighted statistics in theory, finally solve for the velocity bias. $b_v$ determined by our method is not only free of sampling artifact, but also free of cosmic variance. We apply this method to a $\Lambda$CDM N-body simulation of $3072^3$ particles and $1200 Mpc/{\rm h}$ box size. For the first time, we determine the halo velocity bias to $0.1\%$-$1\%$ accuracy. Our major findings are as follows: (1) $b_v\neq 1$ at $k>0.1 h/{\rm Mpc}$. The deviation from unity ($|b_v-1|$) increases with $k$. Depending on halo mass and redshift, it may reach $\mathcal{O}(0.01)$ at $k=0.2 h/{\rm Mpc}$ and $\mathcal{O}(0.05)$ at $k\sim 0.3 h/{\rm Mpc}$. The discovered $b_v\neq 1$ has statistically significant impact on structure growth rate measurement by spectroscopic redshift surveys, including DESI, Euclid and SKA. (2) Both the sign and the amplitude of $b_v-1$ depend on mass and redshift. These results disagree with the peak model prediction in that $b_v$ has much weaker deviation from unity, varies with redshift, and can be bigger than unity. (3) Most of the mass and redshift dependences can be compressed into a single dependence on the halo density bias. Based on this finding, we provide an approximate two-parameter fitting formula.

We describe the theory of Fermi-type acceleration, including first-order Fermi acceleration at a parallel shock front and second-order Fermi acceleration in a test particle limit. Including the theory of the turbulent acceleration and the derivation of the general Fokker-Planck equation, we take into account the basic particle-transport equation and the Fokker-Planck equation with spatially homogeneous isotropic distribution. In the cases of some special physical processes, we construct the particle-transport equations, and compile the analytical or semi-analytical solutions in 11 different cases. Even though, traditionally, one of the electron energy distributions for these different cases can be used to reproduce the multi-wavelength emission of a blazar, due to very long integration times of high-energy observations, a time-integrated electron energy distribution should be introduced in approaches to modeling the multi-wavelength spectral energy distribution of blazars, especially in modeling the flaring state.

The present work considers a plane shock front propagating along a cylindrical jet. Electrons experience the diffusive shock acceleration around the shock front, and subsequently drift away into the downstream flow in which they emit most of their energy. Assuming a proper boundary condition at the interface between the shock zone and the downstream zone, we solve the transport equation for the electrons in the downstream flow zone, where the combined effects of escape, synchrotron and IC cooling in the Thomson regime are taken into account. Using the electron spectrum obtained in this manner we calculate the multi-wavelength spectral energy distribution of Mrk 501 in the synchrotron self-Compton scenario. We check numerically if the Klein-Nishina cross-section could be approximated to the Thomsom regime. We consider whether the model results yield physically reasonable parameters, and further discuss some of implications of the model results. It suggests that the process of diffusive shock acceleration operates in the outflow of Mrk 501.

Post-starburst galaxies can be identified via the presence of prominent Hydrogen Balmer absorption lines in their spectra. We present a comprehensive study of the origin of strong Balmer lines in a volume-limited sample of 189 galaxies with $0.01<z<0.05$, $\log(\mbox{M}_{\star}/\mbox{M}_{\odot})>9.5$ and projected axis ratio $b/a>0.32$. We explore their structural properties, environments, emission lines and star formation histories, and compare them to control samples of star-forming and quiescent galaxies, and simulated galaxy mergers. Excluding contaminants, in which the strong Balmer lines are most likely caused by dust-star geometry, we find evidence for three different pathways through the post-starburst phase, with most events occurring in intermediate-density environments: (1) a significant disruptive event, such as a gas-rich major merger, causing a starburst and growth of a spheroidal component, followed by quenching of the star formation (70% of post-starburst galaxies at $9.5<\log(\mbox{M}_{\star}/\mbox{M}_{\odot})<10.5$ and 60% at $\log(\mbox{M}_{\star}/\mbox{M}_{\odot})>10.5$); (2) at $9.5<\log(\mbox{M}_{\star}/\mbox{M}_{\odot})<10.5$, stochastic star formation in blue-sequence galaxies, causing a weak burst and subsequent return to the blue sequence (30%); (3) at $\log(\mbox{M}_{\star}/\mbox{M}_{\odot})>10.5$, cyclic evolution of quiescent galaxies which gradually move towards the high-mass end of the red sequence through weak starbursts, possibly as a result of a merger with a smaller gas-rich companion (40%). Our analysis suggests that AGN are `on' for $50%$ of the duration of the post-starburst phase, meaning that traditional samples of post-starburst galaxies with strict emission line cuts will be at least $50%$ incomplete due to the exclusion of narrow-line AGN.

The curvature of the $\gamma$-ray spectrum in blazars may reflect the intrinsic distribution of the emitting electron distribution, which will further give some information on the possible acceleration and cooling processes in the emitting region. The $\gamma$-ray spectra of Fermi blazars are normally fitted either by a single power-law (PL) or a log-normal (call Logarithmic Parabola, LP) form. The possible reason for this differnece is not unclear. We statistically explore this issue based on the different observational properties of 1419 Fermi blazars in the 3LAC Clean sample. We find that the $\gamma$-ray flux (100 MeV-100 GeV) and variability index follow bimodal distributions for PL and LP blazars, where $\gamma$-ray flux and variability index show {a positive correlation}. However, the distributions of the $\gamma$-ray luminosity and redshift follow a unimodal distribution. Our results suggest that the bimodal distribution of $\gamma$-ray flux for LP and PL blazars may be not intrinsic and all blazars may have an intrinsic curved $\gamma$-ray spectrum and the PL spectrum is just caused by the fitting effect due to the less photons.

In recent years there has been a lot of interest in discussing frame dependences/independences of the cosmological perturbations under the conformal transformations. This problem has previously been investigated in terms of the covariant approach for a single component universe, and it was found that the covariant approach is very powerful to pick out the perturbative variables which are both gauge and conformal invariant. In this work, we extend the covariant approach to a universe with multicomponent fluids. We find that similar results can be derived, as expected. In addition, some other interesting perturbations are also identified to be conformal invariant, such as entropy perturbation between two different components.

The observed power spectrum in redshift space appears distorted due to the peculiar motion of galaxies, known as redshift-space distortions (RSD). While all the effects in RSD are accounted for by the simple mapping formula from real to redshift spaces, accurately modeling redshift-space power spectrum is rather difficult due to the non-perturbative properties of the mapping. Still, however, a perturbative treatment may be applied to the power spectrum at large-scales, and on top of a careful modeling of the Finger-of-God effect caused by the small-scale random motion, the redshift-space power spectrum can be expressed as a series of expansion which contains the higher-order correlations of density and velocity fields. In our previous work [JCAP 8 (Aug., 2016) 050], we provide a perturbation-theory inspired model for power spectrum in which the higher-order correlations are evaluated directly from the cosmological $N$-body simulations. Adopting a simple Gaussian ansatz for Finger-of-God effect, the model is shown to quantitatively describe the simulation results. Here, we further push this approach, and present an accurate power spectrum template which can be used to estimate the growth of structure as a key to probe gravity on cosmological scales. Based on the simulations, we first calibrate the uncertainties and systematics in the pertrubation theory calculation in a fiducial cosmological model. Then, using the scaling relations, the calibrated power spectrum template is applied to a different cosmological model. We demonstrate that with our new template, the best-fitted growth functions are shown to reproduce the fiducial values in a good accuracy of 1 \% at $k<0.18 \hompc$ for cosmologies with different Hubble parameters.

It has recently been discovered that some, if not all, classical novae emit GeV gamma rays during outburst, but the mechanisms involved in the production of the gamma rays are still not well understood. We present here a comprehensive multi-wavelength dataset---from radio to X-rays---for the most gamma-ray luminous classical nova to-date, V1324 Sco. Using this dataset, we show that V1324 Sco is a canonical dusty Fe-II type nova, with a maximum ejecta velocity of 2600 km s$^{-1}$ and an ejecta mass of few $\times 10^{-5}$ M$_{\odot}$. There is also evidence for complex shock interactions, including a double-peaked radio light curve which shows high brightness temperatures at early times. To explore why V1324~Sco was so gamma-ray luminous, we present a model of the nova ejecta featuring strong internal shocks, and find that higher gamma-ray luminosities result from higher ejecta velocities and/or mass-loss rates. Comparison of V1324~Sco with other gamma-ray detected novae does not show clear signatures of either, and we conclude that a larger sample of similarly well-observed novae is needed to understand the origin and variation of gamma rays in novae.

Usually when applying the mimetic model to the early universe, higher derivative terms are needed to promote the mimetic field to be dynamical. However such models suffer from the ghost and/or the gradient instabilities and simple extensions cannot cure this pathology. We point out in this paper that it is possible to overcome this difficulty by considering the direct couplings of the higher derivatives of the mimetic field to the curvature of the spacetime.

Massive spectroscopic redshift surveys open a promising window to accurately measure peculiar velocity at cosmological distances through redshift space distortion (RSD). In paper I of this series of work we proposed to decompose peculiar velocity into three eigen-modes (v_\delta, v_S and v_B) in order to facilitate the RSD modeling and peculiar velocity reconstruction. In the current paper we measure the dark matter RSD related statistics of the velocity eigen-modes through a set of N-body simulations, including the velocity power spectra, correlation functions, one-point probability distribution functions, cumulants and the damping functions describing the Finger of God effect. (1) The power spectrum measurement shows that these velocity components have distinctly different spatial distribution and redshift evolution. In particular, we measure the window function \tilde{W}(k,z), which describes the impact of nonlinear evolution on the v_\delta-density relation. We confirm that it can induce a significant systematic error of O(10%) in RSD cosmology. We demonstrate that \tilde{W} can be accurately described by a simple fitting formula with one or two free parameters. (2) The correlation function measurement shows that the correlation length is O(100), O(10) and O(1) Mpc for v_\delta, v_S and v_B respectively. These correlation lengths determine where we can treat the velocity fields as spatially uncorrelated. (3) The velocity PDFs and cumulants quantify non-Gaussianities of the velocity fields. We confirm speculation in paper I that v_\delta is largely Gaussian, nevertheless with non-negligible non-Gaussianity, v_B is significantly non-Gaussian. We also measure the damping functions. Despite the observed non-Gaussianities, the damping functions and hence the FOG effect are all well approximated as Gaussian ones at scales of interest.

We report the detection of a largely ionized very-high velocity cloud (VHVC; $v_{\rm LSR}\sim-350$ km/s) toward M33 with the Hubble Space Telescope/Cosmic Origin Spectrograph. The VHVC is detected in OI, CII, SiII, and SiIII absorption along five sightlines separated by ~0.06-0.4 degree. On sub-degree scales, the velocities and ionic column densities of the VHVC remain relatively smooth with standard deviations of +/-14 km/s and +/-0.15 dex between the sightlines, respectively. The VHVC has a metallicity of [OI/HI]=-0.56+/-0.17 dex (Z=0.28+/-0.11 Z$_{\odot}$). Despite the position-velocity proximity of the VHVC to the ionized Magellanic Stream, the VHVC's higher metallicity makes it unlikely to be associated with the Stream, highlighting the complex velocity structure of this region of sky. We investigate the VHVC's possible origin by revisiting its surrounding HI environment. We find that the VHVC may be: (1) a MW CGM cloud, (2) related to a nearby HI VHVC -- Wright's Cloud, or (3) connected to M33's northern warp. Furthermore, the VHVC could be a bridge connecting Wright's Cloud and M33's northern warp, which would make it a Magellanic-like structure in the halo of M33.

Many technologies have been developed to help improve spatial resolution of observational images for ground-based solar telescopes, such as adaptive optics (AO) systems and post-processing reconstruction. As any AO system correction is only partial, it is indispensable to use post-processing reconstruction techniques. In the New Vacuum Solar Telescope (NVST), speckle masking method is used to achieve the diffraction limited resolution of the telescope. Although the method is very promising, the computation is quite intensive, and the amount of data is tremendous, requiring several months to reconstruct observational data of one day on a high-end computer. To accelerate image reconstruction, we parallelize the program package on a high performance cluster. We describe parallel implementation details for several reconstruction procedures. The code is written in C language using Message Passing Interface (MPI) and optimized for parallel processing in a multi-processor environment. We show the excellent performance of parallel implementation, and the whole data processing speed is about 71 times faster than before. Finally, we analyze the scalability of the code to find possible bottlenecks, and propose several ways to further improve the parallel performance. We conclude that the presented program is capable of executing in real-time reconstruction applications at NVST.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

A model-dependent method is proposed to determine the location of the $\gamma$-ray emitting region for a given flat spectrum radio quasar (FSRQ). In the model, the extra-relativistic electrons are injected at the base of the jet and non-thermal photons are produced by both synchrotron radiation and inverse-Comtpon (IC) scattering in the energy dissipation region. The target photons dominating inverse-Comtpon scattering originate from both synchrotron photons and external ambient photon fields, and the energy density of external radiation field is a function of the distance between the position of dissipation region and a central super-massive black hole, and their spectra are seen in the comoving frame. Moreover, the energy dissipation region could be determined by the model parameter through reproducing the $\gamma$-ray spectra. Such a model is applied to reproduce the quasi-simultaneous multi-wavelength observed data for 36 FSRQs. In order to define the width of the broad-line region shell and dusty molecular torus shell, a simple numerical constraint is used to determine the outer boundary of the broad-line region and dusty molecular torus. Our results show that 1) the $\gamma$-ray emitting regions are located at the range from 0.1 pc to 10 pc; 2) the $\gamma$-ray emitting regions are located outside the broad-line regions and within the dusty molecular tori; and 3) the $\gamma$-ray emitting region are located closer to the dusty molecular torus ranges than the broad-line regions. Therefore, it may concluded that a direct evidence for the \emph{far site} scenario could be obtained on the basis of the model results.

We report the detection of accreting ionized gas at the disk-halo interface of the nearby galaxy M33. We analyze HST/COS absorption-line spectra of seven ultraviolet-bright stars evenly distributed across the disk of M33. We find Si IV absorption components consistently redshifted relative to the bulk M33's ISM absorption along all the sightlines. The Si IV detection indicates an enriched, disk-wide, ionized gas inflow toward the disk. This inflow is most likely multi-phase as the redshifted components can also be observed in ions with lower ionization states (e.g., S II, P II, Fe II, Si II). Kinematic modeling of the inflow is consistent with an accreting layer at the disk-halo interface of M33, which has an accretion velocity of 110$^{+15}_{-20}$ km s$^{-1}$ at a distance of 1.5$^{+1.0}_{-1.0}$ kiloparsec above the disk. The modeling indicates a total mass of $\sim3.9\times10^7$ M$_{\odot}$ for the accreting material at the disk-halo interface on the near side of the M33 disk , with an accretion rate of $\sim2.9$ M$_{\odot}$ yr$^{-1}$. The high accretion rate and the level of metal-enrichment suggest the inflow is likely to be the fall back of M33 gas from a galactic fountain and/or the gas pulled loosed during a close interaction between M31 and M33. Our study of M33 is the first to unambiguously reveal the existence of a disk-wide, ionized gas inflow beyond the Milky Way, providing a better understanding of gas accretion in the vicinity of a galaxy disk.

We investigate rapid TeV flaring in markarian 501 in the frame of a time-dependent one-zone synchrotron self-Compton (SSC) model. Inthis model, electrons are accelerated to extra-relativistic energy through the stochastic particle acceleration and evolve with the time, nonthermal photons are produced by both synchrotron and inverse Comtpon scattering off synchrotron photons. Moreover, nonthermal photons during a pre-flare are produced by the relativistic electrons in the steady state and those during a flare are produced by the electrons whose injection rate is changed at some time interval. We apply the model to the rapid flare of Markarian 501 on July 9, 2005 and obtain the multi-wavelength spectra during the pre-flare and during the flare. Our results show that the time-dependent properties of flares can be reproduced by varying the injection rate of electrons and a clear canonical anti-clockwise-loop can be given.

We develop a model for a possible origin of hard very high energy spectra from a distant blazar. In the model, both the primary photons produced in the source and secondary photons produced outside the source contribute to the observed high energy $\gamma$-rays emission. That is, the primary photons are produced in the source through the synchrotron self-Compton (SSC) process, and the secondary photons are produced outside the source through high energy protons interaction with the background photons along the line of sight. We apply the model to a characteristic case was the very high energy (VHE) $\gamma$-ray emissions in distant blazar 1ES 1101-232. Assuming a suitable electron and proton spectra, we obtain excellent fits to observed spectra of distant blazar 1ES 1101-232. This indicated that the surprisingly low attenuation of high energy $\gamma$-rays, especially for the shape of the very high energy $\gamma$-rays tail of the observed spectra, can be explained by secondary $\gamma$-rays produced in interactions of cosmic-ray protons with background photons in the intergalactic space.

The mapping of dark matter clustering from real space to redshift space introduces the anisotropic property to the measured density power spectrum in redshift space, known as the redshift space distortion effect. The mapping formula is intrinsically non-linear, which is complicated by the higher order polynomials due to indefinite cross correlations between the density and velocity fields, and the Finger-of-God effect due to the randomness of the peculiar velocity field. Whilst the full higher order polynomials remain unknown, the other systematics can be controlled consistently within the same order truncation in the expansion of the mapping formula, as shown in this paper. The systematic due to the unknown non-linear density and velocity fields is removed by separately measuring all terms in the expansion directly using simulations. The uncertainty caused by the velocity randomness is controlled by splitting the FoG term into two pieces, 1) the "one-point" FoG term being independent of the separation vector between two different points, and 2) the "correlated" FoG term appearing as an indefinite polynomials which is expanded in the same order as all other perturbative polynomials. Using 100 realizations of simulations, we find that the Gaussian FoG function with only one scale-independent free parameter works quite well, and that our new mapping formulation accurately reproduces the observed 2-dimensional density power spectrum in redshift space at the smallest scales by far, up to $k\sim 0.2h$Mpc, considering the resolution of future experiments.

We investigate MeV-GeV $\gamma$-ray outbursts in 3C 273 in the frame of a time-dependent one-zone synchrotron self-Compton (SSC) model. In this model, electrons are accelerated to extra-relativistic energy through the stochastic particle acceleration and evolve with the time, nonthermal photons are produced by both synchrotron and inverse Compton scattering of synchrotron photons. Moreover, nonthermal photons during a quiescent are produced by the relativistic electrons in the steady state and those during a outburst are produced by the electrons whose injection rate is changed at some time interval. We apply the model to two exceptionally luminous $\gamma$-ray outbursts observed by the Fermi-LAT from 3C 273 in September, 2009 and obtain the multi-wavelength spectra during the quiescent and during the outburst states, respectively. Our results show that the time-dependent properties of outbursts can be reproduced by adopting the appropriate injection rate function of the electron population.

We present an ordinary case in the momentum diffusion equation for the electron spectrum evolution and investigate the energy spectra and time dependent properties of flare in Mrk 421 in the frame of a time-dependent one-zone synchrotron self-Compton model. In this model, electrons are accelerated to extra-relativistic energy through the strong magnetic turbulence and evolve with time, and non-thermal photons are produced by both synchrotron and inverse Comtpon scattering off synchrotron photons. Moreover, non-thermal photons during a pre-flare are produced by the relativistic electrons in the steady state and those during a flare are produced by the electrons whose injection rate is changed during some time interval. We apply the model to the energy spectra and time-dependent properties of flares in Mrk 421 by reproducing the pre-burst spectrum of the source and varying the injection rate in the Bohm diffusion (q=1) and the hard sphere approximation (q=2) case, respectively. Our results show that Bohm diffusion case leads to hard photon spectra, and hard sphere case approximation seems to reproduce the energy spectra and the time dependent properties of flare still better.

We present medium-resolution, near-ultraviolet VLT/FLAMES observations of the star USNO-A0600-15865535. We adapt a standard method of stellar typing to our measurement of the shape of the Balmer epsilon absorption line to demonstrates that USNO-A0600-15865535 is a blue horizontal branch star, residing in the lower stellar halo at a distance of 4.4 kpc from the Sun. We measure the H & K lines of singly-ionized calcium and find two isolated velocity components, one originating in the disk, and one associated with high-velocity cloud complex WD. This detection demonstrated that complex WD is closer than ~4.4 kpc and is the first distance constraint on the +100 km/s Galactic complex of clouds. We find that Complex WD is not in corotation with the Galactic disk as has been assumed for decades. We examine a number of scenarios, and find that the most likely is that Complex WD was ejected from the solar neighborhood and is only a few kpc from the Sun.

We employ a length-dependent conical jet model for the jet structure and emission properties of flat spectrum radio quasar 3C 279 in the steady state. In the model, ultra-relativistic leptons are injected at the base of jet and propagate along the jet structure, non-thermal photons are produced by both synchrotron emission and inverse Comtpon scattering off synchrotron photons and external soft photons at each segment of jet. We derive the total energy spectra contribution through integrating every segment. We apply the model to the quasi-simultaneous multi-wavelength observed data of two quiescent epoches. Using the observed radio data of source, we determine the length of jet $L\sim 100$ pc, and magnetic field $B_{0}\sim 0.1-1$ G at the base of jet. Assuming a steady geometry of the jet structure and suitable physical parameters, we reproduce the multi-wavelength spectra during two quiescent observed epoches, respectively. Our results show that the initial $\gamma$-ray emission site takes place at $\sim 0.5$ pc from black hole.

Context. Very high-energy (VHE) $\gamma$-ray measurements of distant TeV blazars can be nicely explained by TeV spectra induced by ultra high-energy cosmic rays. Aims. We develop a model for a plausible origin of hard spectra in distant TeV blazars. Methods. In the model, the TeV emission in distant TeV blazars is dominated by two mixed components. The first is the internal component with the photon energy around 1 TeV produced by inverse Compton scattering of the relativistic electrons on the synchrotron photons (SSC) with a correction for extragalactic background light absorbtion and the other is the external component with the photon energy more than 1 TeV produced by the cascade emission from high-energy protons propagating through intergalactic space. Results. Assuming suitable model parameters, we apply the model to observed spectra of distant TeV blazars of 1ES 0229+200. Our results show that 1) the observed spectrum properties of 1ES 0229+200, especially the TeV $\gamma$-ray tail of the observed spectra, could be reproduced in our model and 2) an expected TeV $\gamma$-ray spectrum with photon energy $>$1 TeV of 1ES 0229+200 should be comparable with the 50-hour sensitivity goal of the Cherenkov Telescope Array (CTA) and the differential sensitivity curve for the one-year observation with the Large High Altitude Air Shower Observatory (LHAASO). Conclusions. We argue that strong evidence for the Bethe-Heitler cascades along the line of sight as a plausible origin of hard spectra in distant TeV blazars could be obtained from VHE observations with CTA, LHAASO, HAWC, and HiSCORE.

PKS 1424+240 is a distant very high energy gamma-ray BL Lac object with redshift $z=0.601$. It was found that pure synchrotron self-Compton (SSC) process normally need extreme input parameters (e.g., very low magnetic field intensity and extraordinarily large Doppler factor) to explain its multi-wavelength spectral energy distributions (SEDs). To avoid the extreme model parameters, different models have been proposed (e.g., two-zone SSC model or lepto-hadronic model). In this work, we employ the traditional one-zone leptonic model after including a weak external Compton component to re-explore the simultaneous multi-wavelength SEDs of PKS 1424+240 in both high (2009) and low (2013) states. We find that the input parameters of magnetic field and Doppler factor are roughly consistent with those of other BL Lacs if a weak external photon field from either broad line region (BLR) or the dust torus. However, the required energy density of seed photons from BLR or torus is about 3 orders of magnitude less than that constrained in luminous quasars (e.g., flat-spectrum radio quasars, FSRQs). This result suggests that the BLR/torus in BL Lacs is much weaker than that of luminous FSRQs (but not fully disappear), and the inverse-Compton of external photons from BLR/torus may still play a role even in high synchrotron peaked blazars.

The importance of shocks in nova explosions has been highlighted by Fermi's discovery of \gamma-ray producing novae. Over three years of multi-band VLA radio observations of the 2010 nova V1723 Aql show that shocks between fast and slow flows within the ejecta led to the acceleration of particles and the production of synchrotron radiation. Soon after the start of the eruption, shocks in the ejecta produced an unexpected radio flare, resulting in a multi-peaked radio light curve. The emission eventually became consistent with an expanding thermal remnant with mass $2 \times 10^{-4} M_\odot$ and temperature $10^4$ K. However, during the first two months, the $\gtrsim 10^6$ K brightness temperature at low frequencies was too high to be due to thermal emission from the small amount of X-ray producing shock-heated gas. Radio imaging showed structures with velocities of 400 km s$^{-1}$ (d/6 kpc) in the plane of the sky, perpendicular to a more elongated 1500 km s$^{-1}$ (d/6 kpc) flow. The morpho-kinematic structure of the ejecta from V1723 Aql appears similar to nova V959 Mon, where collisions between a slow torus and a faster flow collimated the fast flow and gave rise to \gamma -ray producing shocks. Optical spectroscopy and X-ray observations of V1723 Aql during the radio flare are consistent with this picture. Our observations support the idea that shocks in novae occur when a fast flow collides with a slow collimating torus. Such shocks could be responsible for hard X-ray emission, \gamma -ray production, and double-peaked radio light curves from some classical novae.

We have monitored the BL Lac object Mrk 501 in optical $V$, $R$ and $I$ bands from 2010 to 2015. For Mrk 501, the presence of strong host galaxy component can affect the results of photometry. After subtracting the host galaxy contributions, the source shows intraday and long-term variabilities for optical flux and color indices. The average variability amplitudes of $V$, $R$ and $I$ bands are $22.05\%, 22.25\%, 23.82\%$ respectively, and the value of duty cycle 14.87 per cent. A minimal variability timescale of 106 minutes is detected. No significant time lag between $V$ and $I$ bands is found on one night. The bluer-when-brighter (BWB) trend is dominant for Mrk 501 on intermediate, short and intraday timescales which supports the shock-in-jet model. For the long timescale, Mrk 501 in different state can have different BWB trend. The corresponding results of non-correcting host galaxy contributions are also presented.

On 7 January 2014 an X1.2 flare and CME with a radial speed $\approx$2500 km s$^{-1}$ was observed from near an active region close to disk center. This led many forecasters to estimate a rapid arrival at Earth ($\approx$36 hours) and predict a strong geomagnetic storm. However, only a glancing CME arrival was observed at Earth with a transit time of $\approx$49 hours and a $K_{\rm P}$ geomagnetic index of only $3-$. We study the interplanetary propagation of this CME using the ensemble Wang-Sheeley-Arge (WSA)-ENLIL+Cone model, that allows a sampling of CME parameter uncertainties. We explore a series of simulations to isolate the effects of the background solar wind solution, CME shape, tilt, location, size, and speed, and the results are compared with observed in-situ arrivals at Venus, Earth, and Mars. Our results show that a tilted ellipsoid CME shape improves the initial real-time prediction to better reflect the observed in-situ signatures and the geomagnetic storm strength. CME parameters from the Graduated Cylindrical Shell model used as input to WSA--ENLIL+Cone, along with a tilted ellipsoid cloud shape, improve the arrival-time error by 14.5, 18.7, 23.4 hours for Venus, Earth, and Mars respectively. These results highlight that CME orientation and directionality with respect to observatories play an important role in understanding the propagation of this CME, and for forecasting other glancing CME arrivals. This study also demonstrates the importance of three-dimensional CME fitting made possible by multiple viewpoint imaging.

A profound assumption in peculiar velocity cosmology is $b_v=1$ at sufficiently large scales, where $b_v$ is the volume weighted halo(galaxy) velocity bias with respect to the matter velocity field. However, this fundamental assumption has not been robustly verified in numerical simulations. Furthermore, it is challenged by structure formation theory (BBKS, 1986, ApJ; Desjacques and Sheth, 2010, PRD), which predicts the existence of velocity bias (at least for proto-halos) due to the fact that halos reside in special regions (local density peaks). The major obstacle to measure the volume weighted velocity from N-body simulations is an unphysical sampling artifact. It is entangled in the measured velocity statistics and becomes significant for sparse populations. With recently improved understanding of the sampling artifact (Zhang, Zheng and Jing, 2015, PRD; Zheng, Zhang and Jing, 2015, PRD), for the first time we are able to {\it appropriately correct this sampling artifact and then robustly measure the volume weighted halo velocity bias}. (1) We verify $b_v=1$ within $2\%$ model uncertainty at $k\lesssim 0.1h/$Mpc and $z=0$-$2$ for halos of mass $\sim 10^{12}$-$10^{13} h^{-1} M_\odot$, and, therefore, consolidates a foundation of the peculiar velocity cosmology. (2) We also find statistically significant signs of $b_v\neq 1$ at $k\gtrsim 0.1h/$Mpc. Unfortunately, whether this is real or caused by residual sampling artifact requires further investigation. Nevertheless, cosmology based on $k\gtrsim 0.1h/$Mpc velocity data shall be careful this potential velocity bias.

Ensemble modeling of CMEs provides a probabilistic forecast of CME arrival time which includes an estimation of arrival time uncertainty from the spread and distribution of predictions and forecast confidence in the likelihood of CME arrival. The real-time ensemble modeling of CME propagation uses the WSA-ENLIL+Cone model installed at the CCMC and executed in real-time. The current implementation evaluates the sensitivity of WSA-ENLIL+Cone model simulations of CME propagation to initial CME parameters. We discuss the results of real-time ensemble simulations for a total of 35 CME events between January 2013 - July 2014. For the 17 events where the CME was predicted to arrive at Earth, the mean absolute arrival time prediction error was 12.3 hours, which is comparable to the errors reported in other studies. For predictions of CME arrival at Earth the correct rejection rate is 62% and the false-alarm rate is 38%. The arrival time was within the range of the ensemble arrival predictions for 8 out of 17 events. The Brier Score for CME arrival predictions is 0.15 (where 1 is a perfect forecast), indicating that on average, the predicted likelihood of CME arrival is fairly accurate. The reliability of ensemble CME arrival predictions is heavily dependent on the initial distribution of CME input parameters, particularly the median and spread. Preliminary analysis of the probabilistic forecasts suggests undervariability, indicating that these ensembles do not sample a wide enough spread in CME input parameters. Prediction errors can also arise from ambient model parameters, the accuracy of the solar wind background derived from coronal maps, or other model limitations. Finally, predictions of the Kp geomagnetic index differ from observed values by less than one for 11 out of 17 of the ensembles and Kp prediction errors computed from the mean predicted Kp show a mean absolute error of 1.3.

We assess the fraction of the Milky Way's circumgalactic medium (CGM) eluding detection due to its velocity being similar to gas in the disk. This is achieved using synthetic observations of the CGM in a simulated MW-mass galaxy that shows similar CGM kinematics to the MW and external L$\sim$L$_*$ galaxies. As viewed by a mock observer at a location similar to the Sun, only 50$\%$ (by mass) of the gas moves at high velocity ($|v_{\rm LSR}|\geq$100 km s$^{-1}$ or $|v_{\rm DEV}|\geq$50 km s$^{-1}$) in the simulated CGM and would be observable. The low velocity gas is thermodynamically similar to the high velocity gas, indicating the 50$\%$ observable fraction is applicable to spectral lines from the radio to the ultraviolet. We apply the observable mass fraction (50$\%$) to current estimates of the MW's CGM, and find a corrected total mass of 2.8$\times$10$^{8} M_{\odot}$ for gas below 10$^6$K within $\sim15$ kpc (this excludes the Magellanic System). This is less than the total mass of the CGM extending out to $\sim$150 kpc in other L$\sim$L$_*$ galaxies. However, we find similar OVI column densities when the discrepancy in path length between the MW and external galaxies is considered. The coherent spatial and kinematic distribution of low velocity gas in the simulated CGM suggests that current HI observations of the MW's CGM may miss large low velocity HI complexes. In addition, current mass estimates of the MW's CGM based on high-velocity line observations with distance constraints may miss a non-negligible fraction of gas in the outer halo which can be obscured if it moves at a velocity similar to the gas in the lower halo.

Cosmology based on large scale peculiar velocity preferes volume weighted velocity statistics. However, measuring the volume weighted velocity statistics from inhomogeneously distributed galaxies (simulation particles/halos) suffer from an inevitable and significant sampling artifact. We study this sampling artifact in the velocity power spectrum measured by the nearest-particle (NP) velocity assignment method(Zheng et al. 2013, PRD). We derive the analytical expression of leading and higher order terms. We find that the sampling artifact suppresses the $z=0$ E-mode velocity power spectrum by $\sim 10\%$ at $k=0.1h/$Mpc , for samples with number density $10^{-3}({\rm Mpc}/h)^{-3}$. This suppression becomes larger for larger $k$ and for sparser samples. We argue that, this source of systematic errors in peculiar velocity cosmology, albeit severe, can be self-calibrated in the framework of our theoretical modelling. We also work out the sampling artifact in the density-velocity cross power spectrum measurement. More robust evaluation of related statistics through simulations will be presented in a companion paper (Zheng, Zhang & Jing, 2015, PRD). We also argue that similar sampling artifact exists in other velocity assignment methods and hence must be carefully corrected to avoid systematic bias in peculiar velocity cosmology.

Measuring the volume weighted velocity power spectrum suffers from a severe systematic error, due to imperfect sampling of the velocity field from inhomogeneous distribution of dark matter particles/halos in simulations or galaxies with velocity measurement. This "sampling artifact" depends on both the mean particle number density $\bar{n}_P$ and the intrinsic large scale structure (LSS) fluctuation in the particle distribution. (1) We report robust detection of this sampling artifact in N-body simulations. It causes $\sim 12$% underestimation of the velocity power spectrum at $k=0.1$h/Mpc for samples with $\bar{n}_P=6\times10^{-3}$ (Mpc/h)$^{-3}$. This systematic underestimation increases with decreasing $\bar{n}_P$ and increasing $k$. Its dependence on the intrinsic LSS fluctuations is also robustly detected. (2) All these findings are expected by our theoretical modelling in paper I \cite{Zhang14}. In particular, the leading order theoretical approximation agrees quantitatively well with simulation result for $\bar{n}_P\gtrsim6\times 10^{-4}$(Mpc/h)$^{-3}$. Furthermore, we provide an ansatz to take high order terms into account. It improves the model accuracy to $\lesssim1$% at $k\lesssim0.1$h/Mpc over 3 orders of magnitude in $\bar{n}_P$ and over typical LSS clustering from $z=0$ to $z=2$. (3) The sampling artifact is determined by the deflection ${\bf D}$ field, which is straightforwardly available in both simulations and data of galaxy velocity. Hence the sampling artifact in the velocity power spectrum measurement can be self-calibrated within our framework. By applying such self-calibration in simulations, it becomes promising to determine the {\it real} large scale velocity bias of $10^{13}M_\odot$ halos with $\sim 1$% accuracy, and that of lower mass halos by better accuracy. ...[abridged]

We have developed and implemented an iterative algorithm of flux calibration for the LAMOST Spectroscopic Survey of the Galactic anti-center (LSS-GAC). For a given LSS-GAC plate, the spectra are first processed with a set of nominal spectral response curves (SRCs) and used to derive initial stellar atmospheric parameters (effective temperature $T_{\rm eff}$, surface gravity log\,$g$ and metallicity [Fe/H]) as well as dust reddening $E(B-V)$ of all targeted stars. For each of the sixteen spectrographs, several F-type stars of good signal-to-noise ratios (SNRs) are then selected as flux standard stars for further, iterative spectral flux calibration. Comparison of spectrophotometric colours, deduced from the flux-calibrated spectra, with the photometric measurements yields average differences of 0.02$\pm$0.07 and $-$0.04$\pm$0.09\,mag for the $(g-r)$ and $(g-i)$, respectively. The relatively large negative offset in $(g-i)$ is due to the fact that we have opted not to correct for the telluric bands, most notably the atmospheric A-band in the wavelength range of $i$-band. Comparison of LSS-GAC multi-epoch observations of duplicate targets indicates that the algorithm has achieved an accuracy of about 10 per cent in relative flux calibration for the wavelength range 4000 -- 9000\,\AA. The shapes of SRC deduced for the individual LAMOST spectrographs are found to vary by up to 30 per cent for a given night, and larger for different nights, indicating that the derivation of SRCs for the individual plates is essential in order to achieve accurate flux calibration for the LAMOST spectra.

It is suggested by Chen {\it et al.} that the Earth's surface Ultraviolet irradiance ($280-400$ nm) could influence the Earth's surface temperature variation by "Highly Excited Water Vapor" (HEWV) effect. In this manuscript, we reconstruct the developing history of the ozone layer depth variation from 1860 to 2011 based on the HEWV effect. It is shown that the reconstructed ozone layer depth variation correlates with the observational variation from 1958 to 2005 very well ($R=0.8422$, $P>99.9\%$). From this reconstruction, we may limit the spectra band of the surface Ultraviolet irradiance referred in HEWV effect to Ultraviolet B ($280-320$ nm).

The pure self-gravitating system in this paper refers to a multi-body gaseous system where the self-gravity plays a dominant role and the intermolecular interactions can be neglected. Therefore its total mass must be much more than a limit mass, the minimum mass of the system exhibiting long-range nature. Thee method to estimate the limit mass is then proposed. The nonequlibrium stationary state in the system is identical to the Tsallis equilibrium state, at which the Tsallis entropy approaches to its maximum. On basis of this idea, we introduce a new concept of the temperature whose expression includes the gravitational potential and therefore we call it gravitational temperature. Accordingly, the gravitational thermal capacity is also introduced and it can be used to verify the thermodynamic stability of the astrophysical systems.

Whether natural factors could interpret the rise of the Earth's surface temperature is still controversial. Though numerous recent researches have reported apparent correlations between solar activity and the Earth's climate, solar activity has encountered a big problem when describing the rapid global warming after 1970s. Our investigation shows the good positive correlations between the Earth's surface Ultraviolet irradiance (280-400 nm) and the Earth's surface temperature both in temporal and spatial variations by analyzing the global surface Ultraviolet irradiance (280-400 nm) and global surface temperature data from 1980-1999. The rise of CO$_2$ cannot interpret the good positive correlations, and we could even get an opposite result to the good correlations when employing the rise of CO$_2$ to describe the relation between them. Based on the good positive correlations, we suggest a new effect, named "Highly Excited Water Vapor" (HEWV) effect, which can interpret how the Sun influences the Earth's surface temperature reasonably, including the rapid warming after 1970s.

Classical novae are the most common astrophysical thermonuclear explosions, occurring on the surfaces of white dwarf stars accreting gas from companions in binary star systems. Novae typically expel ~10^(-4) solar masses of material at velocities exceeding 1,000 kilometres per second. However, the mechanism of mass ejection in novae is poorly understood, and could be dominated by the impulsive flash of thermonuclear energy, prolonged optically thick winds, or binary interaction with the nova envelope. Classical novae are now routinely detected in gigaelectronvolt gamma-ray wavelengths, suggesting that relativistic particles are accelerated by strong shocks in the ejecta. Here we report high-resolution radio imaging of the gamma-ray-emitting nova V959 Mon. We find that its ejecta were shaped by the motion of the binary system: some gas was expelled rapidly along the poles as a wind from the white dwarf, while denser material drifted out along the equatorial plane, propelled by orbital motion. At the interface between the equatorial and polar regions, we observe synchrotron emission indicative of shocks and relativistic particle acceleration, thereby pinpointing the location of gamma-ray production. Binary shaping of the nova ejecta and associated internal shocks are expected to be widespread among novae, explaining why many novae are gamma-ray emitters.

We present two possible physical explanations of the nonextensive parameter in a self-gravitating system. One is related to the detailed balance in such system. With the detailed balance, the statistical balance of molecular energy is reached, under which the radiation energy "stored" in the molecules becomes important. Then the relation between the nonextensive parameter and the storage coefficient which measures the ability of molecules to store the radiation energy is constructed. The other explanation is from one dimensionless quantity called (logarithmic) temperature gradient defined in the stellar physics. In view of this, we find that in order to keep the nuclear reactions alive inside the Sun, the ratio of radiation pressure to the total pressure in the core must be more than 0.2833.

Ultracompact dark matter minihalos (UCMHs) would be formed during the earlier universe if there were large density perturbations. If the dark matter can decay into the standard model particles, such as neutrinos, these objects would become the potential astrophysical sources and could be detected by the related instruments, such as IceCube. In this paper, we investigate the neutrino signals from the nearby UCMHs due to the gravitino dark matter decay and compare these signals with the background neutrino flux which is mainly from the atmosphere to get the constraints on the abundance of UCMHs.

With the assumption of local Tsallis equilibrium, the newly defined gravitational temperature is calculated in the solar interior, whose distribution curve can be divided into three parts, the solar core region, radiation region and convection region, in excellent agreement with the solar hierarchical structure. By generalizing the Fourier law, one new mechanism of heat conduction, based on the gradient of the gravitational temperature, is introduced into the astrophysical system. This mechanism is related to the self-gravity of such self-gravitating system whose characteristic scale is large enough. It perhaps plays an important role in the astrophysical system which, in the solar interior, leads to the heat accumulation at the bottom of the convection layer and then motivates the convection motion.

3C 66A was monitored by the BATC(Beijing-Arizona-Taipei-Connecticut) telescope from 2005 to 2008,1994 observations were obtained on 89 nights. Detailed research and analysis was performed on these observations in this paper. A long term burst occurred in the whole light curve. No intra-day variability was claimed in our campaign by intra-night light curve analysis. Time lag of shorter wavelenth preceding longer wavelength was shown by correlation analysis. The results showed that the optical spectral shape turned flatter when the source brightened, and the spectral variability indicator was bigger on shorter time-scale as determined by the color indices variation analysis.

The radio light curves of novae rise and fall over the course of months to years, allowing for detailed observations of the evolution of the nova shell. However, the main parameter determined by radio models of nova explosions - the mass of the ejecta - often seems to exceed theoretical expectations by an order of magnitude. With the recent technological improvements on the Karl G. Jansky Very Large Array (VLA), new observations can test the assumptions upon which ejecta mass estimates are based. Early observations of the classical nova V1723 Aql showed an unexpectedly rapid rise in radio flux density and a distinct bump in the radio light curve on the rise to radio maximum, which is inconsistent with the simple model of spherical ejecta expelled in a single discrete event. This initial bump appears to indicate the presence of shocked material in the outer region of the ejected shell, with the emission from the shocks fading over time. We explore possible origins for this emission and its relation to the mass loss history of the nova. The evolution of the radio spectrum also reveals the density profile, the mass of the ejected shell, and other properties of the ejecta. These observations comprise one of the most complete, longterm set of multi-wavelength radio observations for any classical nova to date.

Quasars have reigned as the most luminous and distant objects in the Universe. These concepts were based on the hypothesis that the redshifts of quasars are obeying the Hubble low. This hypothesis has little serious competition today, but has met some difficulties when using it to determine the distances of quasars. However, the majority of the redshifts of quasars can be explained by the soft photon process proposed by Zheng (2013), then quasars are not located at the distances as that calculated from observed redshifts according to the Hubble low. Therefore distances of quasars can be reduced to a reasonable range, and the difficulties have met by quasars disappeared.

In 1929 Zwicky proposed a tired light theory to interpret the Hubble law (Hubble 1929). The key of the tired light theory is to interpret the observed redshift of galaxy as the non-Doppler effect. But the derivation of the dark matter in the Coma cluster proposed by Zwicky in 1933 was based on the interpretation that redshifts of galaxies were the Doppler effect, and the non-Doppler effect was not considered at all. However, if there is a reasonable non-Doppler effect and the great majority of the observed redshifts of galaxies in the Coma cluster can be interpreted by the non-Doppler effect, then it's not needed to introduce the dark matter in the Coma cluster to keep the stability of the cluster. In this paper it is shown that the great majority of the redshifts of galaxies in the Coma cluster are caused by the non-Doppler redshift proposed by Zheng (2013), and a numerical estimation of the redshifts of galaxies in the Coma cluster is presented.

Anomalous redshifts of some galactic objects such as binary stars, early-type stars in the solar neighborhood, and O stars in a star clusters are discussed. It is shown that all these phenomena have a common characteristic, that is, the redshifts of stars increase as the temperature rises. This characteristic cannot be explained by means of the Doppler Effect but can by means of the soft-photon process proposed by Yijia Zheng (arXiv:1305.0427 [astro-ph.HE]).

The quiet solar coronal heating problem and the observed center-to-limb wavelength variations of the solar lines (limb effect) can be explained. In this paper the quantitative calculations for these two phenomena are presented.

The dilute ionized gas is very popular in the Universe. Usually only the Compton interactions, the "Sunyaev-Zel'dovich" effect, were considered while photons propagated in this medium. In this paper the "soft-photon process" is considered. Due to the soft photons emitted during the propagation of a photon in the dilute ionized gas, the main photon (propagating in the original direction) will be redshifted. The formula to calculate this redshift is derived.

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 $\delta^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]

Upon their formation, dynamically cool (collapsing) star clusters will, within only a few million years, achieve stellar mass segregation for stars down to a few solar masses, simply because of gravitational two-body encounters. Since binary systems are, on average, more massive than single stars, one would expect them to also rapidly mass segregate dynamically. Contrary to these expectations and based on high-resolution Hubble Space Telescope observations, we show that the compact, 15-30 Myr-old Large Magellanic Cloud cluster NGC 1818 exhibits tantalizing hints at the >= 2 sigma level of significance (> 3 sigma if we assume a power-law secondary-to-primary mass-ratio distribution) of an increasing fraction of F-star binary systems (with combined masses of 1.3-1.6 Msun) with increasing distance from the cluster center, specifically between the inner 10 to 20" (approximately equivalent to the cluster's core and half-mass radii) and the outer 60 to 80". If confirmed, this will offer support of the theoretically predicted but thus far unobserved dynamical disruption processes of the significant population of 'soft' binary systems---with relatively low binding energies compared to the kinetic energy of their stellar members---in star clusters, which we have access to here by virtue of the cluster's unique combination of youth and high stellar density.

An antenna array devoted to the autonomous radio-detection of high energy cosmic rays is being deployed on the site of the 21 cm array radio telescope in XinJiang, China. Thanks in particular to the very good electromagnetic environment of this remote experimental site, self-triggering on extensive air showers induced by cosmic rays has been achieved with a small scale prototype of the foreseen antenna array. We give here a detailed description of the detector and present the first detection of extensive air showers with this prototype.

Many quasi-simultaneous optical observations of 17 blazars are obtained from previous papers published over the last 19 years in order to investigate the spectral slope variability and understand the radiation mechanism of blazars. The long-period dereddened optical spectral slopes are calculated. We analyse the average spectral slope distribution, which suggests that the spectra of flat spectrum radio quasars (FSRQs) and high energy peaked BL Lac objects (HBLs) are probably deformed by other emission components. The average spectral slopes of low energy peaked BL Lac objects(LBLs), which scatter around 1.5, show a good accordance with the synchrotron self-Compton (SSC) loss-dominated model. We present and discuss the variability between the spectral slope and optical luminosity. The spectra of all HBLs and LBLs get flatter when they turn brighter, while for FSRQs this trend does not exist or may even be reversed. This phenomenon may imply that there is a thermal contribution to the optical spectrum for FSRQs. For the FSRQ 1156+295, there is a hint that the slope gets flatter at both the brightest and faintest states. Our result shows that three subclasses locate in different regions in the pattern of slope variability indicator versus average spectral slope. The relativistic jet mechanism is supported by the significant correlation between the optical Doppler factor and the average spectral slope.

The neutral hydrogen 21-cm spectral line (1420.4 MHz) and the four 18-cm lines of the hydroxyl molecule (1612-1720 MHz) are observable at redshifts which put their measured line frequencies well below their protected frequency bands. Part of the redshift ranges (z = 0.171-0.477 for HI and z = 0.37-0.73 for OH) fall in the 960 to 1215 MHz band that is allocated to aircraft navigation. Most of the signals in this band are pulsed emissions of low duty cycle so much of the time between pulses is interference free. This paper outlines the structure and measured properties of signals in this band and demonstrates a signal processing strategy that is effective at removing the pulsed signals from spectra at sensitivities produced by several hours of integration.