The overwhelming foreground causes severe contamination on the detection of 21-cm signal during the Epoch of Reionization (EoR). Among various foreground components, the Galactic free-free emission is less studied, so that its impact on the EoR observation remains unclear. To better constrain this emission, we perform the Monte Carlo simulation of H$\alpha$ emission, which comprises direct and scattered H$\alpha$ radiation from HII regions and warm ionized medium (WIM). The positions and radii of HII regions are quoted from the WISE HII catalog, and the WIM is described by an axisymmetric model. The scattering is off dust and free electrons that are realized by applying an exponential fitting to the HI4PI HI map and an exponential disk model, respectively. The simulated H$\alpha$ intensity, the Simfast21 software, and the latest SKA1-Low layout configuration are employed to simulate the SKA "observed" images of Galactic free-free emission and the EoR signal. By analyzing the one-dimensional power spectra, we find that the Galactic free-free emission can be about $10^{5.4}$-$10^{2.1}$, $10^{5.0}$-$10^{1.7}$, and $10^{4.3}$-$10^{1.1}$ times more luminous than the EoR signal on scales of $0.1~{\rm Mpc^{-1}} < k < 2~{\rm Mpc^{-1}}$ in the 116-124, 146-154, and 186-194 MHz frequency bands, respectively. We further calculate the two-dimensional power spectra inside the EoR window and show that the power leaked by Galactic free-free emission can still be significant, as the power ratios can reach about $110\%$-$8000\%$, $30\%$-$2400\%$, and $10\%$-$250\%$ on scales of $0.5~{\rm Mpc^{-1}} \lesssim k \lesssim 1~{\rm Mpc^{-1}}$ in three frequency bands. Therefore, we indicate that the Galactic free-free emission should be carefully treated in future EoR detections.

Based on the dynamics of single scalar field slow-roll inflation and the theory of reheating, we investigate the generalized natural inflationary (GNI) model. Concretely, we give constraints on the scalar spectral index $n_{s}$ and tensor-to scalar ratio $r$ for $\Lambda$CDM $+r$ model according to the latest data from Plack 2018 TT,TE,EE+lowE+lensing (P18) and BICEP2/Keck 2015 season (BK15), i.e., $n_{s}=0.9659\pm0.0044$ at $68\%$ confidence level (CL) and $r<0.0623$ at $95\%$CL. We find that the GNI model is favored by P18 plus BK15 in the ranges of $\log_{10}(f/M_{p})=0.62^{+0.17}_{-0.18}$ and $m=0.35^{+0.13}_{-0.23}$ at $68\%$CL. In addition, the corresponding predictions of the general and two-phase reheating are respectively discussed. It follows that the parameter $m$ has the significant effect on the model behaviors.

The star formation rate function (SFRF) and specific star formation rate function (sSFRF) from the observation are impacted by the Eddington bias, due to the uncertainties on the estimated SFR. We develop a novel method to correct the Eddington bias and obtained the intrinsic SFRF and sSFRF from the Sloan Digital Sky Survey Data Release 7. The intrinsic SFRF is in good agreement with measurements from previous data in the literature that relied on UV SFRs but its high star-forming end is slightly lower than those IR and radio tracers. We demonstrate that the intrinsic sSFRF from SDSS has a bi-modal form with the one peak found at ${\rm sSFR \sim 10^{-9.7} yr^{-1}}$ representing the star-forming objects while the other peak is found at ${\rm sSFR \sim 10^{-12} yr^{-1}}$ representing the quenched population. Furthermore, we compare our observations with the predictions from the IllustrisTNG and Illustris simulations and show that the ``TNG'' model performs much better than its predecessor. However, we show that the simulated SFRF and cosmic star formation density (CSFRD) of TNG simulations are highly dependent on resolution, reflecting the limitations of the model and today state-of-the-art simulations. We demonstrate that the bi-modal, two peaked sSFRF implied by the SDSS observations does not appear in TNG regardless of the adopted box-size or resolution. This tension reflects the need for inclusion of an additional efficient quenching mechanism to the TNG model.

We report new HI observations of four z$\sim$0.05 star-forming galaxies undertaken during the commissioning phase of the Five-hundred-meter Aperture Spherical Radio Telescope (FAST). FAST is the largest single-dish telescope with a 500 meter aperture and a 19-Beam receiver. Exploiting the unprecedented sensitivity provided by FAST, we aim to study the atomic gas, via the HI 21cm emission line, in low-$z$ star-forming galaxies taken from the Valpara\'iso ALMA/APEX Line Emission Survey (VALES) project. Together with previous ALMA CO($J=1-0$) observations, the HI data provides crucial information to measure the gas mass and dynamics. As a pilot HI survey, we targeted four local star-forming galaxies at $z\sim0.05$. In particular, one of them has already been detected in HI by the Arecibo Legacy Fast ALFA survey (ALFALFA), allowing a careful comparison. We use an ON-OFF observing approach that allowed us to reach an rms of 0.7mJy/beam at a 1.7km/s velocity resolution within only 20 minutes ON-target integration time. We demonstrate the great capabilities of the FAST 19-beam receiver for pushing the detectability of the HI emission line of extra-galactic sources. The HI emission line detected by FAST shows good consistency with the previous ALFALFA results. Our observations are put in context with previous multi-wavelength data to reveal the physical properties of these low-$z$ galaxies. We find that the CO($J=1-0$) and HI emission line profiles are similar. The dynamical mass estimated from the HI data is an order of magnitude higher than the baryon mass and the dynamical mass derived from the CO observations, implying that the mass probed by dynamics of HI is dominated by the dark matter halo. In one case, a target shows an excess of CO($J=1-0$) in the line centre, which can be explained by an enhanced CO($J=1-0$) emission induced by a nuclear starburst showing high velocity dispersion.

We present timing solutions and analyses of 11 pulsars discovered by the Five-hundred-meter Aperture Spherical radio Telescope (FAST). These pulsars were discovered using an ultra-wide bandwidth receiver in drift-scan observations made during the commissioning phase of FAST, and were then confirmed and timed using the 64-m Parkes Radio Telescope. Each pulsar has been observed over a span of at least one year. Highlighted discoveries include PSR J0344-0901, which displays mode-changing behaviour and may belong to the class of so-called `swooshing' pulsars (alongside PSRs B0919+06 and B1859+07); PSR J0803-0942, whose emission is almost completely linearly polarised; and PSRs J1900-0134 and J1945+1211, whose well defined polarisation angle curves place stringent constraints on their emission geometry. We further discuss the detectability of these pulsars by earlier surveys, and highlight lessons learned from our work in carrying out confirmation and monitoring observations of pulsars discovered by a highly sensitive telescope, many of which may be applicable to next-generation pulsar surveys. This paper marks one of the first major releases of FAST-discovered pulsars, and paves the way for future discoveries anticipated from the Commensal Radio Astronomy FAST Survey (CRAFTS).

Fast radio bursts (FRBs) are mysterious millisecond-duration radio transients observed from extragalactic distances. Magnetars have been long speculated as the possible engine to power repeating bursts from FRB sources, but no convincing evidence has been collected so far. Recently, a Galactic magnetar dubbed Soft Gamma-ray Repeater (SGR) J1935+2154 entered an active phase by emitting intense soft gamma-ray bursts. One fast radio burst with two peaks (hereafter FRB 200428) and a luminosity comparable to those of extragalactic FRBs was detected from the source in association with a soft $\gamma$-ray / hard X-ray flare. Here we report an eight-hour targeted observational campaign with the Five-hundred-meter Aperture Spherical radio Telescope (FAST) in four observing sessions, assisted by multi-wavelength (optical and hard X-rays) observations. During the third session, twenty-nine SGR bursts were detected by the {\it Fermi} Gamma-ray Burst Monitor (GBM). Throughout the observing period, especially during the scrutinized epochs when SGR bursts reached the Earth, no single dispersed pulsed emission was detected by FAST. This places a fluence upper limit eight orders of magnitude deeper than the fluence of FRB 200428. Our results suggest that FRB -- SGR burst associations are rather rare. FRBs may be highly relativistic and geometrically beamed or may carry narrow spectra with the characteristic frequencies of the most bursts far outside the FAST band. More intriguingly, it is also possible that the physical condition to achieve coherent radiation in SGR bursts is contrived. Only under extreme conditions could an FRB be made in association with an SGR burst.

To reach the full potential for the next generation of weak lensing surveys, it is necessary to mitigate the contamination of intrinsic alignments (IA) of galaxies in the observed cosmic shear signal. The self calibration (SC) of intrinsic alignments provides an independent method to measure the IA signal from the survey data and the photometric redshift information. It operates differently from the marginalization method based on IA modeling. In this work, we present the first application of SC to the KiDS450 data and the KV450 data, to split directly the intrinsic shape - galaxy density (Ig) correlation signal and the gravitational shear - galaxy density (Gg) correlation signal, using the information from photometric redshift (photo-z). We achieved a clear separation of the two signals and performed several validation tests. Our measured signals are found to be in general agreement with the KiDS450 cosmic shear best-fit cosmology, for both lensing and IA measurements. For KV450, we use partial (high-z) data, and our lensing measurements are also in good agreement with KV450 cosmic shear best-fit, however, our IA signal suggests a larger IA amplitude. We discussed the impact of photo-z quality on IA detection and several other potential systematic biases. Finally, we discuss the potential application of the information extracted for both the lensing signal and the IA signal in future surveys.

The overwhelming foreground contamination hinders the accurate detection of the 21-cm signal of neutral hydrogen during the Epoch of Reionization (EoR). Among various foreground components, the Galactic free-free emission is less studied, so that its impact on the EoR observations remains unclear. In this work, we employ the observed $\rm H\alpha$ intensity map with the correction of dust absorption and scattering, the Simfast21 software, and the latest SKA1-Low layout configuration to simulate the SKA "observed" images of Galactic free-free emission and the EoR signal. By calculating the one-dimensional power spectra from the simulated image cubes, we find that the Galactic free-free emission is about $10^{3.5}$-$10^{2.0}$, $10^{3.0}$-$10^{1.3}$, and $10^{2.5}$-$10^{1.0}$ times more luminous than the EoR signal on scales of $0.1~\rm Mpc^{-1} < k < 2~\rm Mpc^{-1}$ in the $116$-$124$, $146$-$154$, and $186$-$194$ ${\rm MHz}$ frequency bands. We further analyse the two-dimensional power spectra inside the properly defined EoR window and find that the leaked Galactic free-free emission can still cause non-negligible contamination, as the ratios of its power (amplitude squared) to the EoR signal power can reach about $200\%$, $60\%$, and $15\%$ on scales of $1.2~\rm Mpc^{-1}$ in three frequency bands, respectively. Therefore, we conclude that the Galactic free-free emission, as a severe contaminating foreground component, needs to be carefully treated in the forthcoming deep EoR observations.

Based mainly on X-ray observations, studies are made on interactions between the intra-cluster medium (ICM) in clusters of galaxies and their member galaxies. Through (magneto)hydrodynamic and gravitational channels, the moving galaxies are expected to drag the ICM around them, and transfer to the ICM some fraction of their dynamical energies on cosmological time scales. This hypothesis is in line with several observations, including the possible cosmological infall of galaxies towards the cluster center, found over redshifts of z~1 to z~0. Further assuming that the energy lost by the galaxies is first converted into ICM turbulence and then dissipated, this picture can explain the subsonic and uniform ICM turbulence, measured with Hitomi in the core region of the Perseus cluster. The scenario may also explain several other unanswered problems regarding clusters of galaxies, including what prevents the ICM from the expected radiative cooling, how the various mass components in nearby clusters have attained different radial distributions, and how a thermal stability is realized between hot and cool ICM components that co-exist around cD galaxies. This view is also considered to pertain to the general scenario of galaxy evolution, including their environmental effects.

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.

The weak-field approximation implying linear relationship between Stokes $V/I$ and longitudinal magnetic field, $B_{\Vert}$, often suffers from the saturation observed in strong magnetic field regions such as sunspot umbra. In this work we intend to improve the magnetic field observations carried out by the \textit{Solar Magnetic Field Telescope} (SMFT) at the Huairou Solar Observing Station, China. We propose to use non-linear relationship between Stokes $V/I$ and $B_{\Vert}$ to derive magnetic field. To determine the form of the relationship we perform cross-calibration of SMFT data and magnetograms provided by the \textit{Helioseismic and Magnetic Imager} on board the \textit{Solar Dynamics Observatory}. The algorithm of the magnetic field deriving is described in details. We show that using non-linear relationship between Stokes $V/I$ and $B_{\Vert}$ allows eliminating the magnetic field saturation inside sunspot umbra.

The emission from PSR J1107-5907 is erratic. Sometimes the radio pulse is undetectable, at other times the pulsed emission is weak, and for short durations the emission can be very bright. In order to improve our understanding of these state changes, we have identified archival data sets from the Parkes radio telescope in which the bright emission is present, and find that the emission never switches from the bright state to the weak state, but instead always transitions to the off state. Previous work had suggested the identification of the off state as an extreme manifestation of the weak state. However, the connection between the off and bright emission reported here suggests that the emission can be interpreted as undergoing only two emission states: a bursting state consisting of both bright pulses and nulls as well as the weak-emission state.

As China's first X-ray astronomical satellite, the Hard X-ray Modulation Telescope (HXMT), which was dubbed as Insight-HXMT after the launch on June 15, 2017, is a wide-band (1-250 keV) slat-collimator-based X-ray astronomy satellite with the capability of all-sky monitoring in 0.2-3 MeV. It was designed to perform pointing, scanning and gamma-ray burst (GRB) observations and, based on the Direct Demodulation Method (DDM), the image of the scanned sky region can be reconstructed. Here we give an overview of the mission and its progresses, including payload, core sciences, ground calibration/facility, ground segment, data archive, software, in-orbit performance, calibration, background model, observations and some preliminary results.

Corona cooling was detected previously from stacking a series of short type-I bursts occurred during the low/had state of atoll outburst. Type-I bursts are hence regarded as sharp probe to our better understanding on the basic property of the corona. The launch of the first Chinese X-ray satellite Insight-HXMT has large detection area at hard X-rays which provide almost unique chance to move further in this research field. We report the first detection of the corona cooling by Insight-HXMT from single short type-I burst showing up during {\bf flare} of 4U 1636-536. This type-I X-ray burst has a duration of $\sim$13 seconds and hard X-ray shortage is detected with significance 6.2~$\sigma$ in 40-70 keV. A cross-correlation analysis between the lightcurves of soft and hard X-ray band, shows that the corona shortage lag the burst emission by 1.6 $\pm$1.2~s. These results are consistent with those derived previously from stacking a large amount of bursts detected by RXTE/PCA within a series of {\bf flares} of 4U 1636-536. Moreover, the broad bandwidth of Insight-HXMT allows as well for the first time to infer the burst influence upon the continuum spectrum via performing the spectral fitting of the burst, which ends up with the finding that hard X-ray shortage appears at around 40 keV in the continuum spectrum. These results suggest that the evolution of the corona along with the outburst{\bf /flare} of NS XRB may be traced via looking into a series of embedded type-I bursts by using Insight-HXMT.

We present a study of timing properties of the accreting pulsar 2S 1417-624 observed during its 2018 outburst, based on Swift/BAT, Fermi/GBM, Insight-HXMT and NICER observations. We report a dramatic change of the pulse profiles with luminosity. The morphology of the profile in the range 0.2-10.0keV switches from double to triple peaks at $\sim2.5$ $\rm \times 10^{37}{\it D}_{10}^2\ erg\ s^{-1}$ and from triple to quadruple peaks at $\sim7$ $\rm \times 10^{37}{\it D}_{10}^2\ erg\ s^{-1}$. The profile at high energies (25-100keV) shows significant evolutions as well. We explain this phenomenon according to existing theoretical models. We argue that the first change is related to the transition from the sub to the super-critical accretion regime, while the second to the transition of the accretion disc from the gas-dominated to the radiation pressure-dominated state. Considering the spin-up as well due to the accretion torque, this interpretation allows to estimate the magnetic field self-consistently at $\sim7\times 10^{12}$G.

The precise measurement of the spectrum of protons, the most abundant component of the cosmic radiation, is necessary to understand the source and acceleration of cosmic rays in the Milky Way. This work reports the measurement of the cosmic ray proton fluxes with kinetic energies from 40 GeV to 100 TeV, with two and a half years of data recorded by the DArk Matter Particle Explorer (DAMPE). This is the first time an experiment directly measures the cosmic ray protons up to ~100 TeV with a high statistics. The measured spectrum confirms the spectral hardening found by previous experiments and reveals a softening at ~13.6 TeV, with the spectral index changing from ~2.60 to ~2.85. Our result suggests the existence of a new spectral feature of cosmic rays at energies lower than the so-called knee, and sheds new light on the origin of Galactic cosmic rays.

We report on analysis of observations of the bright transient X-ray pulsar \src obtained during its 2017-2018 giant outburst with Insight-HXMT, \emph{NuSTAR}, and \textit{Swift} observatories. We focus on the discovery of a sharp state transition of the timing and spectral properties of the source at super-Eddington accretion rates, which we associate with the transition of the accretion disk to a radiation pressure dominated (RPD) state, the first ever directly observed for magnetized neutron star. This transition occurs at slightly higher luminosity compared to already reported transition of the source from sub- to super-critical accretion regime associate with onset of an accretion column. We argue that this scenario can only be realized for comparatively weakly magnetized neutron star, not dissimilar to other ultra-luminous X-ray pulsars (ULPs), which accrete at similar rates. Further evidence for this conclusion is provided by the non-detection of the transition to the propeller state in quiescence which strongly implies compact magnetosphere and thus rules out magnetar-like fields.

When photons from distant galaxies and stars pass through our neighboring environment, the wavelengths of the photons would be shifted by our local gravitational potential. This local gravitational redshift effect can potentially have an impact on the measurement of cosmological distance-redshift relation. Using available supernovae data, Wojtak et al [1] found seemingly large biases of cosmological parameters for some extended models (non-flat $\Lambda$CDM, $w$CDM, etc.). Huang [2] pointed out that, however, the biases can be reduced to a negligible level if cosmic microwave background (CMB) data are added to break the strong degeneracy between parameters in the extended models. In this article we forecast the cosmological bias due to local gravitational redshifts for a future WFIRST-like supernovae survey. We find that the local gravitational redshift effect remains negligible, provided that CMB data or some future redshift survey data are added to break the degeneracy between parameters.

We present X-ray flux and spectral analyses of the three pointed Suzaku observations of the TeV high synchrotron peak blazar Mrk 421 taken throughout its complete operational duration. The observation taken on 5 May 2008 is, at 364.6 kiloseconds (i.e., 101.3 hours), the longest and most evenly sampled continuous observation of this source, or any blazar, in the X-ray energy 0.8 - 60 keV until now. We found large amplitude intra-day variability in all soft and hard bands in all the light curves. The discrete correction function analysis of the light curves in soft and hard bands peaks on zero lag, showing that the emission in hard and soft bands are cospatial and emitted from the same population of leptons. The hardness ratio plots imply that the source is more variable in the harder bands compared to the softer bands. The source is harder-when-brighter, following the general behavior of high synchrotron peak blazars. Power spectral densities of all three light curves are red noise dominated, with a range of power spectra slopes. If one assumes that the emission originates very close to the central super massive black hole, a crude estimate for its mass, of ~ 4 * 10^{8} M_{\odot}, can be made; but if the variability is due to perturbations arising there that are advected into the jet and are thus Doppler boosted, substantially higher masses are consistent with the quickest seen variations. We briefly discuss the possible physical mechanisms most likely responsible for the observed flux and spectral variability.

In this work, we report the in-orbit demonstration of X-ray pulsar navigation with Insight-Hard X-ray Modulation Telescope (Insight-HXMT), which was launched on Jun. 15th, 2017. The new pulsar navigation method 'Significance Enhancement of Pulse-profile with Orbit-dynamics' (SEPO) is adopted to determine the orbit with observations of only one pulsar. In this test, the Crab pulsar is chosen and observed by Insight-HXMT from Aug. 31th to Sept. 5th in 2017. Using the 5-day-long observation data, the orbit of Insight-HXMT is determined successfully with the three telescopes onboard - High Energy X-ray Telescope (HE), Medium Energy X-ray Telescope (ME) and Low Energy X-ray Telescope (LE) - respectively. Combining all the data, the position and velocity of the Insight-HXMT are pinpointed to within 10 km (3 sigma) and 10 m/s (3 sigma), respectively.

The recently discovered neutron star transient Swift J0243.6+6124 has been monitored by {\it the Hard X-ray Modulation Telescope} ({\it Insight-\rm HXMT). Based on the obtained data, we investigate the broadband spectrum of the source throughout the outburst. We estimate the broadband flux of the source and search for possible cyclotron line in the broadband spectrum. No evidence of line-like features is, however, found up to $\rm 150~keV$. In the absence of any cyclotron line in its energy spectrum, we estimate the magnetic field of the source based on the observed spin evolution of the neutron star by applying two accretion torque models. In both cases, we get consistent results with $B\rm \sim 10^{13}~G$, $D\rm \sim 6~kpc$ and peak luminosity of $\rm >10^{39}~erg~s^{-1}$ which makes the source the first Galactic ultraluminous X-ray source hosting a neutron star.

Cosmic backreaction as an additional source of the expansion of the universe has been a debate topic since the discovery of cosmic acceleration. The major concern is whether the self interaction of small-scale nonlinear structures would source gravity on very large scales. Gregory Ryskin argued against the additional inclusion of gravitational interaction energy of astronomical objects, whose masses are mostly inferred from gravitational effects and hence should already contain all sources with long-range gravity forces. Ryskin proposed that the backreaction contribution to the energy momentum tensor is instead from the rest of the universe beyond the observable patch. Ryskin's model elegantly solves the fine-tuning problem and is in good agreement with the Hubble diagram of Type Ia supernovae. In this article we revisit Ryskin's model and show that it is {\it inconsistent} with at least one of the following statements: (i) the universe is matter-dominated at low redshift ($z\lesssim 2$); (ii) the universe is radiation-dominated at sufficiently high redshift; (iii) matter density fluctuations are tiny ($\lesssim 10^{-4}$) at the recombination epoch.

The DArk Matter Particle Explorer (DAMPE), a satellite-based cosmic ray and gamma-ray detector, was launched on December 17, 2015, and began its on-orbit operation on December 24, 2015. In this work we document the on-orbit calibration procedures used by DAMPE and report the calibration results of the Plastic Scintillator strip Detector (PSD), the Silicon-Tungsten tracKer-converter (STK), the BGO imaging calorimeter (BGO), and the Neutron Detector (NUD). The results are obtained using Galactic cosmic rays, bright known GeV gamma-ray sources, and charge injection into the front-end electronics of each sub-detector. The determination of the boundary of the South Atlantic Anomaly (SAA), the measurement of the live time, and the alignments of the detectors are also introduced. The calibration results demonstrate the stability of the detectors in almost two years of the on-orbit operation.

The overwhelming foreground contamination is one of the primary impediments to probing the Epoch of Reionization (EoR) through measuring the redshifted 21 cm signal. Among various foreground components, radio halos are less studied and their impacts on the EoR observations are still poorly understood. In this work, we employ the Press-Schechter formalism, merger-induced turbulent re-acceleration model, and the latest SKA1-Low layout configuration to simulate the SKA "observed" images of radio halos. We calculate the one-dimensional power spectra from simulated images and find that radio halos can be about $10^4$, $10^3$ and $10^{2.5}$ times more luminous than the EoR signal on scales of $0.1\,\text{Mpc}^{-1} < k < 2\,\text{Mpc}^{-1}$ in the 120-128, 154-162, and 192-200 MHz bands, respectively. By examining the two-dimensional power spectra inside properly defined EoR windows, we find that the power leaked by radio halos can still be significant, as the power ratios of radio halos to the EoR signal on scales of $0.5\,\text{Mpc}^{-1} \lesssim k \lesssim 1\,\text{Mpc}^{-1}$ can be up to about 230-800%, 18-95%, and 7-40% in the three bands, when the 68% uncertainties caused by the variation of the number density of bright radio halos are considered. Furthermore, we find that radio halos located inside the far side-lobes of the station beam can also impose strong contamination within the EoR window. In conclusion, we argue that radio halos are severe foreground sources and need serious treatments in future EoR experiments.

We develop a novel method to explore the galaxy-halo connection using the galaxy imaging surveys by modeling the projected two-point correlation function measured from the galaxies with reasonable photometric redshift measurements. By assuming a Gaussian form of the photometric redshift errors, we are able to simultaneously constrain the halo occupation distribution (HOD) models and the effective photometric redshift uncertainties. Tests with mock galaxy catalogs demonstrate that this method can successfully recover (within $\sim 1\sigma$) the intrinsic large-scale galaxy bias, as well as the HOD models and the effective photometric redshift uncertainty. This method also works well even for galaxy samples with 10 per cent catastrophic photometric redshift errors.

The spatial-dependent propagation model has been successfully used to explain diverse observational phenomena, including the spectral hardening of cosmic-ray nuclei above $200$ GV, the large-scale dipole anisotropy and the diffusive gamma distribution. In this work, we further apply the spatial-dependent propagation model to both electrons and positrons. To account for the excess of positrons above $10$ GeV, an additional local source is introduced. And we also consider a more realistic spiral distribution of background sources. We find that due to the gradual hardening above $10$ GeV, the hardening of electron spectrum above tens of GeV can be explained in the SDP model and both positron and electron spectra less than TeV energies could be naturally described. The spatial-dependent propagation with spiral-distributed sources could conforms with the total electron spectrum in the whole energy. Meanwhile compared with the conventional model, the spatial-dependent propagation with spiral-distributed sources could produce larger background positron flux, so that the multiplier of background positron flux is $1.42$, which is much smaller than the required value by the conventional model. Thus the shortage of background positron flux could be solved. Furthermore we compute the anisotropy of electron under spatial-dependent propagation model, which is well below the observational limit of Fermi-LAT experiment.

When applying the foreground removal methods to uncover the faint cosmological signal from the epoch of reionization (EoR), the foreground spectra are assumed to be smooth. However, this assumption can be seriously violated in practice since the unresolved or mis-subtracted foreground sources, which are further complicated by the frequency-dependent beam effects of interferometers, will generate significant fluctuations along the frequency dimension. To address this issue, we propose a novel deep-learning-based method that uses a 9-layer convolutional denoising autoencoder (CDAE) to separate the EoR signal. After being trained on the SKA images simulated with realistic beam effects, the CDAE achieves excellent performance as the mean correlation coefficient ($\bar{\rho}$) between the reconstructed and input EoR signals reaches $0.929 \pm 0.045$. In comparison, the two representative traditional methods, namely the polynomial fitting method and the continuous wavelet transform method, both have difficulties in modelling and removing the foreground emission complicated with the beam effects, yielding only $\bar{\rho}_{\text{poly}} = 0.296 \pm 0.121$ and $\bar{\rho}_{\text{cwt}} = 0.198 \pm 0.160$, respectively. We conclude that, by hierarchically learning sophisticated features through multiple convolutional layers, the CDAE is a powerful tool that can be used to overcome the complicated beam effects and accurately separate the EoR signal. Our results also exhibit the great potential of deep-learning-based methods in future EoR experiments.

We present a morphological classification of 14,245 radio active galactic nuclei (AGNs) into six types, i.e., typical Fanaroff--Riley Class I / II (FRI/II), FRI/II-like bent-tailed, X-shaped radio galaxy, and ringlike radio galaxy, by designing a convolutional neural network (CNN) based autoencoder, namely MCRGNet, and applying it to a labeled radio galaxy (LRG) sample containing 1442 AGNs and an unlabeled radio galaxy (unLRG) sample containing 14,245 unlabeled AGNs selected from the Best--Heckman sample. We train MCRGNet and implement the classification task by a three-step strategy, i.e., pre-training, fine-tuning, and classification, which combines both unsupervised and supervised learnings. A four-layer dichotomous tree is designed to classify the radio AGNs, which leads to a significantly better performance than the direct six-type classification. On the LRG sample, our MCRGNet achieves a total precision of $\sim 93\%$ and an averaged sensitivity of $\sim 87\%$, which are better than those obtained in previous works. On the unLRG sample, whose labels have been human-inspected, the neural network achieves a total precision of $\sim 80\%$. Also, using the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) to calculate the $r$-band absolute magnitude ($M_\mathrm{opt}$) and using the flux densities to calculate the radio luminosity ($L_\mathrm{radio}$), we find that the distributions of the unLRG sources on the $L_\mathrm{radio}$--$M_\mathrm{opt}$ plane do not show an apparent redshift evolution and could confirm with a sufficiently large sample that there could not exist an abrupt separation between FRIs and FRIIs as reported in some previous works.

By analyzing the Chandra archival data of the nearby NGC 6338 galaxy group, we identify two X-ray bright clumps (N-clump and S-clump) within the central 100 kpc, and detect an arc-like X-ray brightness discontinuity at the south boundary of the N-clump, which is defined as a cold front with a gas flow Mach number of M<0.8. Furthermore, at the north-east boundary of the S-clump (dominated by galaxy NGC 6338) another X-ray edge is detected that corresponds to a weaker cold front. Therefore, the two clumps are approaching each other approximately from opposite directions, and the group is undergoing a head-on collision that is in a stage of pre-core passage. This merger scenario is also supported by the study of the line-of-sight velocity distribution of the group member galaxies. The merger mass ratio is about 1:1.8 as estimated from the central gas temperature of the two clumps, which suggests the merger is most likely to be a major merger. We also analyze the VLA 1.4 and 4.9 GHz radio data, but we do not detect any extended radio emission that is associated with the merger.

The origin of supermassive black holes (SMBHs) that inhabit the centers of massive galaxies is largely unconstrained. Remnants from supermassive stars (SMSs) with masses around 10,000 solar masses provide the ideal seed candidates, known as direct collapse black holes. However, their very existence and formation environment in the early Universe are still under debate, with their supposed rarity further exacerbating the problem of modeling their ab-initio formation. SMS models have shown that rapid collapse, with an infall rate above a critical value, in metal-free haloes is a requirement for the formation of a proto-stellar core which will then form an SMS. Using a radiation hydrodynamics simulation of early galaxy formation, we show the natural emergence of metal-free haloes both massive enough, and with sufficiently high infall rates, to form an SMS. We find that haloes that are exposed to both a Lyman-Werner intensity of J_LW ~ 3 J_21 and that undergo at least one period of rapid growth early in their evolution are ideal cradles for SMS formation. This rapid growth induces substantial dynamical heating, amplifying the existing Lyman-Werner suppression originating from a group of young galaxies 20 kiloparsecs away. Our results strongly indicate that structure formation dynamics, rather than a critical Lyman-Werner (LW) flux, may be the main driver of massive black hole formation in the early Universe. We find that massive black hole seeds may be much more common in overdense regions of the early Universe than previously considered with a comoving number density up to 10^-3 Mpc^-3.

The cosmic spatial curvature parameter $\Omega_k$ is constrained, primarily by cosmic microwave background (CMB) data, to be very small. Observations of the cosmic distance ladder and the large scale structure can provide independent checks of the cosmic flatness. Such late-universe constraints on $\Omega_k$, however, are sensitive to the assumptions of the nature of dark energy. For minimally coupled scalar-field models of dark energy, the equation of state $w$ has nontrivial dependence on the cosmic spatial curvature $\Omega_k$ (Miao and Huang, 2018). Such dependence has not been taken into account in previous studies of future observational projects. In this paper we use the $w$ parameterization proposed by Miao and Huang, where the dependence of $w$ on $\Omega_k$ is encoded, and perform Fisher forecast on mock data of three benchmark projects: a WFIRST-like Type Ia supernovae survey, an EUCLID-like spectroscopic redshift survey, and an LSST-like photometric redshift survey. We find that the correlation between $\Omega_k$ and $w$ is primarily determined by the data rather than by the theoretical prior. We thus validate the standard approaches of treating $\Omega_k$ and $w$ as independent quantities.

Forming macroscopic solid bodies in circumstellar discs requires local dust concentration levels significantly higher than the mean. Interactions of the dust particles with the gas must serve to augment local particle densities, and facilitate growth past barriers in the metre size range. Amongst a number of mechanisms that can amplify the local density of solids, aerodynamic streaming instability (SI) is one of the most promising. This work tests the physical assumptions of models that lead to SI in protoplanetary disks (PPDs). We conduct laboratory experiments in which we track the three-dimensional motion of spherical solid particles fluidized in a low-pressure, laminar, incompressible, gas stream. The particle sizes span the Stokes-Epstein drag regime transition and the overall dust-to-gas mass density ratio is close to unity. Lambrechts et al. (2016) established the similarity of the laboratory flow to a simplified PPD model flow. We obtain experimental results suggesting an instability due to particle-gas interaction: i) there exist variations in particle concentration in the direction of the mean drag forces; ii) the particles have a tendency to 'catch up' to one another when they are in proximity; iii) particle clumping occurs on very small scales, which implies local enhancements above the background dust-to-gas mass density ratio by factors of several tens; v) the presence of these density enhancements occurs for a mean dust-to-gas mass density ratio approaching or greater than 1; v) we find evidence for collective particle drag reduction when the local particle number density becomes high and when the background gas pressure is high so that the drag is in the continuum regime. The experiments presented here are precedent-setting for observing SI under controlled conditions and may lead to a deeper understanding of how it operates in nature.

We choose the bright compact group HCG 62, which was found to exhibit both excess X-ray emission and high Fe abundance to the southwest of its core, as an example to study the impact of mergers on chemical enrichment in the intragroup medium. We first reanalyze the high-quality Chandra and XMM-Newton archive data to search for the evidence for additional SN II yields, which is expected as a direct result of the possible merger-induced starburst. We reveal that, similar to the Fe abundance, the Mg abundance also shows a high value in both the innermost region and the southwest substructure, forming a high-abundance plateau, meanwhile all the SN Ia and SN II yields show rather flat distributions in $>0.1r_{200}$ in favor of an early enrichment. Then we carry out a series of idealized numerical simulations to model the collision of two initially isolated galaxy groups by using the TreePM-SPH GADGET-3 code. We find that the observed X-ray emission and metal distributions, as well as the relative positions of the two bright central galaxies with reference to the X-ray peak, can be well reproduced in a major merger with a mass ratio of 3 when the merger-induced starburst is assumed. The `best-match' snapshot is pinpointed after the third pericentric passage when the southwest substructure is formed due to gas sloshing. By following the evolution of the simulated merging system, we conclude that the effects of such a major merger on chemical enrichment are mostly restricted within the core region when the final relaxed state is reached.

We formulate a model of the conditional colour-magnitude distribution (CCMD) to describe the distribution of galaxy luminosity and colour as a function of halo mass. It consists of two populations of different colour distributions, dubbed pseudo-blue and pseudo-red, respectively, with each further separated into central and satellite galaxies. We define a global parameterization of these four colour-magnitude distributions and their dependence on halo mass, and we infer parameter values by simultaneously fitting the space densities and auto-correlation functions of 79 galaxy samples from the Sloan Digital Sky Survey defined by fine bins in the colour-magnitude diagram (CMD). The model deprojects the overall galaxy CMD, revealing its tomograph along the halo mass direction. The bimodality of the colour distribution is driven by central galaxies at most luminosities, though at low luminosities it is driven by the difference between blue centrals and red satellites. For central galaxies, the two pseudo-colour components are distinct and orthogonal to each other in the CCMD: at fixed halo mass, pseudo-blue galaxies have a narrow luminosity range and broad colour range, while pseudo-red galaxies have a narrow colour range and broad luminosity range. For pseudo-blue centrals, luminosity correlates tightly with halo mass, while for pseudo-red galaxies colour correlates more tightly (redder galaxies in more massive haloes). The satellite fraction is higher for redder and for fainter galaxies, with colour a stronger indicator than luminosity. We discuss the implications of the results and further applications of the CCMD model.

Previously we identified a new class of early galaxy that we estimate contributes up to 30\% of the ionizing photons responsible for reionization. These are low mass halos in the range $M_h =10^{6.5}-10^{8} M_{\odot}$ that have been chemically enriched by supernova ejecta from prior Pop III star formation. Despite their low star formation rates, these Metal Cooling halos (MCs) are significant sources of ionizing radiation, especially at the onset of reionization, due to their high number density and ionizing escape fractions. Here we present a fully-coupled radiation hydrodynamic simulation of reionization that includes these MCs as well the more massive hydrogen atomic line cooling halos. Our method is novel: we perform halo finding inline with the radiation hydrodynamical simulation, and assign escaping ionizing fluxes to halos using a probability distribution function (PDF) measured from the galaxy-resolving Renaissance Simulations. The PDF captures the mass dependence of the ionizing escape fraction as well as the probability that a halo is actively forming stars. With MCs, reionization starts earlier than if only halos of $10^8 M_{\odot}$ and above are included, however the redshift when reionization completes is only marginally affected as this is driven by more massive galaxies. Because star formation is intermittent in MCs, the earliest phase of reionization exhibits a stochastic nature, with small H II regions forming and recombining. Only later, once halos of mass $\sim 10^9 M_{\odot}$ and above begin to dominate the ionizing emissivity, does reionization proceed smoothly in the usual manner deduced from previous studies. This occurs at $z\approx 10$ in our simulation.

We present the X-ray timing results of the new black hole candidate (BHC) MAXI J1535-571 during its 2017 outburst from Hard X-ray Modulation Telescope (\emph{Insight}-HXMT) observations taken from 2017 September 6 to 23. Following the definitions given by \citet{Belloni2010}, we find that the source exhibits state transitions from Low/Hard state (LHS) to Hard Intermediate state (HIMS) and eventually to Soft Intermediate state (SIMS). Quasi-periodic oscillations (QPOs) are found in the intermediate states, which suggest different types of QPOs. With the large effective area of \emph{Insight}-HXMT at high energies, we are able to present the energy dependence of the QPO amplitude and centroid frequency up to 100 keV which is rarely explored by previous satellites. We also find that the phase lag at the type-C QPOs centroid frequency is negative (soft lags) and strongly correlated with the centroid frequency. By assuming a geometrical origin of type-C QPOs, the source is consistent with being a high inclination system.

Turbulence dissipation is an important process affecting the energy balance in molecular clouds, the birth place of stars. Previously, the rate of turbulence dissipation is often estimated with semi-analytic formulae from simulation. Recently we developed a data analysis technique called core-velocity-dispersion (CVD), which, for the first time, provides direct measurements of the turbulence dissipation rate in Taurus, a star forming cloud. The thus measured dissipation rate of $(0.45\pm 0.05)\times 10^{33} {\rm erg\ s^{-1}}$ is similar to those from dimensional analysis and also consistent with the previous energy injection rate based on molecular outflows and bubbles.

We analyse the spatial distribution of current helicity in solar active regions. A comparison of current helicity maps derived from three different instruments (Helioseismic and Magnetc Imager on board the Solar Dynamics Observatory, SDO/HMI, Spectro-Polarimeter on board the Hinode, and Solar Magnetic Field Telescope at the Huairou Solar Observing Station, China, HSOS/SMFT) is carried out. The comparison showed an excellent correlation between the maps derived from the spaceborne instruments and moderate correlation between the maps derived from SDO/HMI and HSOS/SMFT vector magnetograms. The results suggest that the obtained maps characterize real spatial distribution of current helicity over an active region. To analyse intermittency of current helicity, we traditionally use the high-order structure functions and flatness function approach. The slope of a flatness function within some range of scales - the flatness exponent - is a measure of the degree of intermittency. SDO/HMI vector magnetograms for 3 ARs (NOAA 11158, 12494, and 12673) were used to calculate the flatness exponent time variations. All three ARs exhibited emergence of a new magnetic flux during the observational interval. The flatness exponent indicated the increase of intermittency 12-20 hours before the emergence of a new flux. We suppose that this behaviour can indicate subphotospheric fragmentation or distortion of the pre-existed electric current system by emerging magnetic flux.

Using isochronous mass spectrometry at the experimental storage ring CSRe in Lanzhou, the masses of $^{82}$Zr and $^{84}$Nb were measured for the first time with an uncertainty of $\sim 10$ keV, and the masses of $^{79}$Y, $^{81}$Zr, and $^{83}$Nb were re-determined with a higher precision. %The latter differ significantly from their literature values. The latter are significantly less bound than their literature values. Our new and accurate masses remove the irregularities of the mass surface in this region of the nuclear chart. Our results do not support the predicted island of pronounced low $\alpha$ separation energies for neutron-deficient Mo and Tc isotopes, making the formation of Zr-Nb cycle in the $rp$-process unlikely. The new proton separation energy of $^{83}$Nb was determined to be 490(400)~keV smaller than that in the Atomic Mass Evaluation 2012. This partly removes the overproduction of the $p$-nucleus $^{84}$Sr relative to the neutron-deficient molybdenum isotopes in the previous $\nu p$-process simulations.

This is the second installment for the Large Sky Area Multi-Object Fibre Spectroscopic Telescope (LAMOST) Quasar Survey, which includes quasars observed from September 2013 to June 2015. There are 9024 confirmed quasars in DR2 and 10911 in DR3. After cross-match with the SDSS quasar catalogs and NED, 12126 quasars are discovered independently. Among them 2225 quasars were released by SDSS DR12 QSO catalogue in 2014 after we finalised the survey candidates. 1801 sources were identified by SDSS DR14 as QSOs. The remaining 8100 quasars are considered as newly founded, and among them 6887 quasars can be given reliable emission line measurements and the estimated black hole masses. Quasars found in LAMOST are mostly located at low-to-moderate redshifts, with a mean value of 1.5. The highest redshift observed in DR2 and DR3 is 5. We applied emission line measurements to H$\alpha$, H$\beta$, Mg{\sc ii} and C{\sc iv}. We deduced the monochromatic continuum luminosities using photometry data, and estimated the virial black hole masses for the newly discovered quasars. Results are compiled into a quasar catalog, which will be available online.

In this paper, we will give a general introduction to the project of Ali CMB Polarization Telescope (AliCPT), which is a Sino-US joint project led by the Institute of High Energy Physics (IHEP) and has involved many different institutes in China. It is the first ground-based Cosmic Microwave Background (CMB) polarization experiment in China and an integral part of China's Gravitational Waves Program. The main scientific goal of AliCPT project is to probe the primordial gravitational waves (PGWs) originated from the very early Universe. The AliCPT project includes two stages. The first stage referred to as AliCPT-1, is to build a telescope in the Ali region of Tibet with an altitude of 5,250 meters. Once completed, it will be the worldwide highest ground-based CMB observatory and open a new window for probing PGWs in northern hemisphere. AliCPT-1 telescope is designed to have about 7,000 TES detectors at 90GHz and 150GHz. The second stage is to have a more sensitive telescope (AliCPT-2) with the number of detectors more than 20,000. Our simulations show that AliCPT will improve the current constraint on the tensor-to-scalar ratio $r$ by one order of magnitude with 3 years' observation. Besides the PGWs, the AliCPT will also enable a precise measurement on the CMB rotation angle and provide a precise test on the CPT symmetry. We show 3 years' observation will improve the current limit by two order of magnitude.

The effect of Quantum Gravity (QG) may bring a tiny light speed variation as $v(E)=c(1-E/E_{\rm LV})$, where $E$ is the photon energy and $E_{\rm LV}$ is a Lorentz violation scale. A remarkable regularity was suggested in previous studies to look for the light speed variation from high energy photon events of Gamma Ray Bursts (GRBs). We provide a general analysis on the data of 25 bright GRBs observed by the Fermi Gamma-ray Space Telescope (FGST). Such method allows a completed scan over all possibilities in a more clean and impartial way without any bias compared to previous intuitive analysis. The results show that with the increase in the intrinsic energies of photons, such regularity truly emerges and gradually becomes significant. For photons with intrinsic energies higher than 40~GeV, the regularity exists at a significance of 3-5~$\sigma$ with $E_{\rm LV}=3.6\times 10^{17}~\rm GeV$ determined by the GRB data.

We produce synthetic spectra and observations for metal-free stellar populations and high mass X-ray binaries in the Renaissance Simulations at a redshift of 15. We extend our methodology from the first paper in the series by modelling the production and extinction of emission lines throughout a dusty and metal-enriched interstellar and circum-galactic media extracted from the simulation, using a Monte Carlo calculation. To capture the impact of high-energy photons, we include all frequencies from hard X-ray to far infrared with enough frequency resolution to discern line emission and absorption profiles. The most common lines in our sample in order of their rate of occurrence are Ly$\alpha$, the C IV $\lambda\lambda1548,1551$ doublet, H-$\alpha$, and the Ca II $\lambda\lambda\lambda8498,8542,8662$ triplet. The best scenario for a direct observation of a metal-free stellar population is a merger between two Population III galaxies. In mergers between metal-enriched and metal-free stellar populations, some characteristics may be inferred indirectly. Single Population III galaxies are too dim to be observed photometrically at $z = 15$. Ly$\alpha$ emission is discernible by $JWST$ as an increase in $\rm{J_{200w} - J_{277w}}$ colour off the intrinsic stellar tracks. Observations of metal-free stars will be difficult, though not impossible, with the next generation of space telescopes.

High energy cosmic ray electrons plus positrons (CREs), which lose energy quickly during their propagation, provide an ideal probe of Galactic high-energy processes and may enable the observation of phenomena such as dark-matter particle annihilation or decay. The CRE spectrum has been directly measured up to $\sim 2$ TeV in previous balloon- or space-borne experiments, and indirectly up to $\sim 5$ TeV by ground-based Cherenkov $\gamma$-ray telescope arrays. Evidence for a spectral break in the TeV energy range has been provided by indirect measurements of H.E.S.S., although the results were qualified by sizeable systematic uncertainties. Here we report a direct measurement of CREs in the energy range $25~{\rm GeV}-4.6~{\rm TeV}$ by the DArk Matter Particle Explorer (DAMPE) with unprecedentedly high energy resolution and low background. The majority of the spectrum can be properly fitted by a smoothly broken power-law model rather than a single power-law model. The direct detection of a spectral break at $E \sim0.9$ TeV confirms the evidence found by H.E.S.S., clarifies the behavior of the CRE spectrum at energies above 1 TeV and sheds light on the physical origin of the sub-TeV CREs.

Finding the electromagnetic (EM) counterpart of binary compact star merger, especially the binary neutron star (BNS) merger, is critically important for gravitational wave (GW) astronomy, cosmology and fundamental physics. On Aug. 17, 2017, Advanced LIGO and \textit{Fermi}/GBM independently triggered the first BNS merger, GW170817, and its high energy EM counterpart, GRB 170817A, respectively, resulting in a global observation campaign covering gamma-ray, X-ray, UV, optical, IR, radio as well as neutrinos. The High Energy X-ray telescope (HE) onboard \textit{Insight}-HXMT (Hard X-ray Modulation Telescope) is the unique high-energy gamma-ray telescope that monitored the entire GW localization area and especially the optical counterpart (SSS17a/AT2017gfo) with very large collection area ($\sim$1000 cm$^2$) and microsecond time resolution in 0.2-5 MeV. In addition, \textit{Insight}-HXMT quickly implemented a Target of Opportunity (ToO) observation to scan the GW localization area for potential X-ray emission from the GW source. Although it did not detect any significant high energy (0.2-5 MeV) radiation from GW170817, its observation helped to confirm the unexpected weak and soft nature of GRB 170817A. Meanwhile, \textit{Insight}-HXMT/HE provides one of the most stringent constraints (~10$^{-7}$ to 10$^{-6}$ erg/cm$^2$/s) for both GRB170817A and any other possible precursor or extended emissions in 0.2-5 MeV, which help us to better understand the properties of EM radiation from this BNS merger. Therefore the observation of \textit{Insight}-HXMT constitutes an important chapter in the full context of multi-wavelength and multi-messenger observation of this historical GW event.

The DArk Matter Particle Explorer (DAMPE), one of the four scientific space science missions within the framework of the Strategic Pioneer Program on Space Science of the Chinese Academy of Sciences, is a general purpose high energy cosmic-ray and gamma-ray observatory, which was successfully launched on December 17th, 2015 from the Jiuquan Satellite Launch Center. The DAMPE scientific objectives include the study of galactic cosmic rays up to $\sim 10$ TeV and hundreds of TeV for electrons/gammas and nuclei respectively, and the search for dark matter signatures in their spectra. In this paper we illustrate the layout of the DAMPE instrument, and discuss the results of beam tests and calibrations performed on ground. Finally we present the expected performance in space and give an overview of the mission key scientific goals.

We study the dependence of galaxy clustering on atomic gas mass using a sample of $\sim$16,000 galaxies with redshift in the range of $0.0025<z<0.05$ and HI mass of $M_{\rm HI}>10^8M_{\odot}$, drawn from the 70% complete sample of the Arecibo Legacy Fast ALFA survey. We construct subsamples of galaxies with $M_{\rm HI}$ above different thresholds, and make volume-limited clustering measurements in terms of three statistics: the projected two-point correlation function, the projected cross-correlation function with respect to a reference sample selected from the Sloan Digital Sky Survey, and the redshift-space monopole moment. In contrast to previous studies, which found no/weak HI-mass dependence, we find both the clustering amplitude on scales above a few Mpc and the bias factors to increase significantly with increasing HI mass for subsamples with HI mass thresholds above $10^9M_{\odot}$. For HI mass thresholds below $10^9M_{\odot}$, while the measurements have large uncertainties caused by the limited survey volume and sample size, the inferred galaxy bias factors are systematically lower than the minimum halo bias factor from mass-selected halo samples. The simple halo model, in which galaxy content is only determined by halo mass, has difficulties in interpreting the clustering measurements of the HI-selected samples. We extend the simple model by including the halo formation time as an additional parameter. A model that puts HI-rich galaxies into halos that formed late can reproduce the clustering measurements reasonably well. We present the implications of our best-fitting model on the correlation of HI mass with halo mass and formation time, as well as the halo occupation distributions and HI mass functions for central and satellite galaxies. These results are compared with the predictions from semi-analytic galaxy formation models and hydrodynamic galaxy formation simulations.

We present synthetic observations for the first generations of galaxies in the Universe and make predictions for future deep field observations for redshifts greater than 6. Due to the strong impact of nebular emission lines and the relatively compact scale of HII regions, high resolution cosmological simulations and a robust suite of analysis tools are required to properly simulate spectra. We created a software pipeline consisting of FSPS, Hyperion, Cloudy and our own tools to generate synthetic IR observations from a fully three-dimensional arrangement of gas, dust, and stars. Our prescription allows us to include emission lines for a complete chemical network and tackle the effect of dust extinction and scattering in the various lines of sight. We provide spectra, 2-D binned photon imagery for both HST and JWST IR filters, luminosity relationships, and emission line strengths for a large sample of high redshift galaxies in the Renaissance Simulations. Our resulting synthetic spectra show high variability between galactic halos with a strong dependence on stellar mass, metallicity, gas mass fraction, and formation history. Halos with the lowest stellar mass have the greatest variability in [OIII]/H$\beta$, [OIII] and CIII] while halos with higher masses are seen to show consistency in their spectra and [OIII] equivalent widths (EW) between 1\AA\ and 10\AA. Viewing angle accounted for three-fold difference in flux due to the presence of ionized gas channels in a halo. Furthermore, JWST color plots show a discernible relationship between redshift, color, and mean stellar age.

We in this paper investigate the formation and evolution of primordial black holes (PBHs) in nonsingular bouncing cosmologies. We discuss the formation of PBH in the contracting phase and calculate the PBH abundance as a function of the sound speed and Hubble parameter. Afterwards, by taking into account the subsequent PBH evolution during the bouncing phase, we derive the density of PBHs and their Hawking radiation. Our analysis shows that nonsingular bounce models can be constrained from the backreaction of PBHs.

We present stellar properties (mass, age, radius, distances) of 57 stars from a seismic inference using full-length data sets from Kepler. These stars comprise active stars, planet-hosts, solar-analogs, and binary systems. We validate the distances derived from the astrometric Gaia-Tycho solution. Ensemble analysis of the stellar properties reveals a trend of mixing-length parameter with the surface gravity and effective temperature. We derive a linear relationship with the seismic quantity $\langle r_{02} \rangle$ to estimate the stellar age. Finally, we define the stellar regimes where the Kjeldsen et al (2008) empirical surface correction for 1D model frequencies is valid.

The Kepler space telescope yielded unprecedented data for the study of solar-like oscillations in other stars. The large samples of multi-year observations posed an enormous data analysis challenge that has only recently been surmounted. Asteroseismic modeling has become more sophisticated over time, with better methods gradually developing alongside the extended observations and improved data analysis techniques. We apply the latest version of the Asteroseismic Modeling Portal (AMP) to the full-length Kepler data sets for 57 stars and the Sun, comprising planetary hosts, binaries, solar-analogs, and active stars. From an analysis of the derived stellar properties for the full sample, we identify a variation of the mixing-length parameter with atmospheric properties. We also derive a linear relation between the stellar age and a characteristic frequency separation ratio. In addition, we find that the empirical correction for surface effects suggested by Kjeldsen and coworkers is adequate for solar-type stars that are not much hotter (Teff < 6200 K) or significantly more evolved (logg > 4.2, <Delta_nu> > 80 muHz) than the Sun. Precise parallaxes from the Gaia mission and future observations from TESS and PLATO promise to improve the reliability of stellar properties derived from asteroseismology.

Gamma-ray polarimetry is a new powerful tool to study the processes responsible for the emission from astrophysical sources and the environments in which this emission takes place. Few successful polarimetric measurements have however been performed thus far in the gamma-ray energy band due to the difficulties involved. POLAR is a dedicated polarimeter designed to perform high precision measurements of the polarization of the emission from gamma-ray burst in the 50-500 keV energy range. This new polarimeter is expected to detect approximately 50 gamma-ray bursts per year while performing high precision polarization measurements on approximately 10 bursts per year. The instrument was launched into lower earth orbit as part of the second Chinese space lab, the Tiangong-2, on September 15th 2016 and has been taking data successfully since being switched on one week after. The instrument uses a segmented scintillator array consisting of 1600 plastic scintillator bars, read out by 25 flat-panel multi-anode photomultipliers, to measure the Compton scattering angles of incoming photons. The small segmentation and relatively large uniform effective area allow the instrument to measure the polarization of a large number of transient events, such as gamma-ray bursts, with an unprecedented precision during its two year life-time. The final flight model underwent detailed calibration prior to launch as well as intensive space qualification tests, a summary of which will be presented in this paper. The instrument design will be discussed first followed by an overview of the on-ground tests, finally the in-orbit behavior as measured during the first weeks of the mission will be presented.

Cosmic reionization is thought to be primarily fueled by the first generations of galaxies. We examine their stellar and gaseous properties, focusing on the star formation rates and the escape of ionizing photons, as a function of halo mass, redshift, and environment using the full suite of the {\it Renaissance Simulations} with an eye to provide better inputs to global reionization simulations. This suite, carried out with the adaptive mesh refinement code Enzo, is unprecedented in terms of their size and physical ingredients. The simulations probe overdense, average, and underdense regions of the universe of several hundred comoving Mpc$^3$, each yielding a sample of over 3,000 halos in the mass range $10^7 - 10^{9.5}~\Ms$ at their final redshifts of 15, 12.5, and 8, respectively. In the process, we simulate the effects of radiative and supernova feedback from 5,000 to 10,000 metal-free (Population III) stars in each simulation. We find that halos as small as $10^7~\Ms$ are able to form stars due to metal-line cooling from earlier enrichment by massive Population III stars. However, we find such halos do not form stars continuously. Using our large sample, we find that the galaxy-halo occupation fraction drops from unity at virial masses above $10^{8.5}~\Ms$ to $\sim$50\% at $10^8 ~\Ms$ and $\sim$10\% at $10^7~\Ms$, quite independent of redshift and region. Their average ionizing escape fraction is $\sim$5\% in the mass range $10^8 - 10^9~\Ms$ and increases with decreasing halo mass below this range, reaching 40--60\% at $10^7~\Ms$. Interestingly, we find that the escape fraction varies between 10--20\% in halos with virial masses $\sim 3 \times 10^9~\Ms$. Taken together, our results confirm the importance of the smallest galaxies as sources of ionizing radiation contributing to the reionization of the universe.

Due to their long mean free paths, X-rays are expected to have many significant impacts globally on the properties of the intergalactic medium (IGM) by their heating and ionizing processes on large scales. At high redshifts, X-rays from Population (Pop) III binaries might have important effects on cosmic reionization and the Lyman alpha forest. As a continuation of our previous work on Pop III binary X-rays (Xu et al. 2014), we use the Pop III distribution and evolution from the Renaissance Simulations, a suite of self-consistent cosmological radiation hydrodynamics simulations of the formation of the first galaxies, to calculate the X-ray luminosity density and background over the redshift range 20 > z > 7.6. As we find that Pop III star formation continues at a low, nearly constant rate to the end of reionization, X-rays are being continuously produced at significant rates compared to other possible X-ray sources, such as AGNs and normal X-ray binaries during the same period of time. We estimate that Pop III binaries produce approximately 6 eV of energy in the X-rays per hydrogen atom. We calculate the X-ray background for different monochromatic photon energies. KeV X-rays redshift and accumulate to produce a strong X-ray background spectrum extending to roughly 500 eV. The X-ray background is strong enough to heat the IGM to ~ 1000 K and to ionize a few percent of the neutral hydrogen. These effects are important for an understanding of the neutral hydrogen hyperfine transition 21-cm line signatures, the Ly alpha forest, and optical depth of the CMB to Thomson scattering.

We present a catalog of 624 radio sources detected around the North Celestial Pole (NCP) with the 21 Centimeter Array (21CMA), a radio interferometer dedicated to the statistical measurement of the epoch of reionization (EoR). The data are taken from a 12 h observation made on 2013 April 13, with a frequency coverage from 75 to 175 MHz and an angular resolution of ~ 4 arcmin. The catalog includes flux densities at eight sub-bands across the 21CMA bandwidth and provides the in-band spectral indices for the detected sources. To reduce the complexity of interferometric imaging from the so-called "w" term and ionospheric effects, the present analysis are restricted to the east-west baselines within 1500 m only. The 624 radio sources are found within 5 degrees around the NCP down to ~ 0.1 Jy. Our source counts are compared, and also exhibit a good agreement, with deep low-frequency observations made recently with the GMRT and MWA. In particular, for fainter radio sources below ~ 1 Jy, we find a flattening trend of source counts towards lower frequencies. While the thermal noise (~0.4 mJy) is well controlled to below the confusion limit, the dynamical range (~10^4) and sensitivity of current 21CMA imaging is largely limited by calibration and deconvolution errors, especially the grating lobes of very bright sources, such as 3C061.1, in the NCP field which result from the regular spacings of the 21CMA. We note that particular attention should be paid to the extended sources, and their modeling and removals may constitute a large technical challenge for current EoR experiments. Our analysis may serve as a useful guide to design of next generation low-frequency interferometers like the Square Kilometre Array.

The effect of quantum gravity can bring a tiny light speed variation which is detectable through energetic photons propagating from gamma ray bursts (GRBs) to an observer such as the space observatory. Through an analysis of the energetic photon data of the GRBs observed by the Fermi Gamma-ray Space Telescope (FGST), we reveal a surprising regularity of the observed time lags between photons of different energies with respect to the Lorentz violation factor due to the light speed energy dependence. Such regularity suggests a linear form correction of the light speed $v(E)=c(1-E/E_{\rm LV})$, where $E$ is the photon energy and $E_{\rm LV}=(3.60 \pm 0.26) \times 10^{17}~ \rm GeV$ is the Lorentz violation scale measured by the energetic photon data of GRBs. The results support an energy dependence of the light speed in cosmological space.

Additional physics beyond standard hydrodynamics is needed to fully model the intracluster medium (ICM); however, as we move to more sophisticated models, it is important to consider the role of magnetic fields and the way the fluid approximation breaks down. This paper represents a first step towards developing a self-consistent model of the ICM by characterizing the statistical properties of magnetic fields in cosmological simulations of galaxy clusters. We find that plasma conditions are largely homogeneous across a range of cluster masses and relaxation states. We also find that the magnetic field length scales are resolution dependent and not based on any particular physical process. Energy transfer mechanisms and scales are also identified, and imply the existence of small scale dynamo action. The scales of the small scale dynamo are resolution limited and driven by numerical resistivity and viscosity.

It is postulated in Einstein's relativity that the speed of light in vacuum is a constant for all observers. However, the effect of quantum gravity could bring an energy dependence of light speed. Even a tiny speed variation, when amplified by the cosmological distance, may be revealed by the observed time lags between photons with different energies from astrophysical sources. From the newly detected long gamma ray burst GRB~160509A, we find evidence to support the prediction for a linear form modification of light speed in cosmological space.

Faint red galaxies in the Sloan Digital Sky Survey show a puzzling clustering pattern in previous measurements. In the two-point correlation function (2PCF), they appear to be strongly clustered on small-scales, indicating a tendency to reside in massive haloes as satellite galaxies. However, their weak clustering on large scales suggests that they are more likely to be found in low mass haloes. The interpretation of the clustering pattern suffers from the large sample variance in the 2PCF measurements, given the small volume of the volume-limited sample of such faint galaxies. We introduce a method to improve the clustering measurements of faint galaxies by making a full use of a flux-limited sample to obtain volume-limited measurements with an increased effective volume. In the improved 2PCF measurements, the fractional uncertainties on large-scales drop by more than 40 per cent, and the strong contrast between small-scale and large-scale clustering amplitudes seen in previous work is no longer prominent. From halo occupation distribution modelling of the measurements, we find that a considerable fraction of faint red galaxies to be satellites in massive haloes, a senario supported by the strong covariance of small-scale 2PCF measurements and the relative spatial distribution of faint red galaxies and luminous galaxies. However, the satellite fraction is found to be degenerate with the slope of the distribution profile of satellites in inner haloes. We compare the modelling results with semi-analytic model predictions and discuss the implications.

In this work we propose a new diagnostic to segregate cool core (CC) clusters from non-cool core (NCC) clusters by studying the two-dimensional power spectra of the X-ray images observed with the Chandra X-ray observatory. Our sample contains 41 members ($z=0.01\sim 0.54$), which are selected from the Chandra archive when a high photon count, an adequate angular resolution, a relatively complete detector coverage, and coincident CC-NCC classifications derived with three traditional diagnostics are simultaneously guaranteed. We find that in the log-log space the derived image power spectra can be well represented by a constant model component at large wavenumbers, while at small wavenumbers a power excess beyond the constant component appears in all clusters, with a clear tendency that the excess is stronger in CC clusters. By introducing a new CC diagnostic parameter, i.e., the power excess index (PEI), we classify the clusters in our sample and compare the results with those obtained with three traditional CC diagnostics. We find that the results agree with each other very well. By calculating the PEI values of the simulated clusters, we find that the new diagnostic works well at redshifts up to 0.5 for intermediately sized and massive clusters with a typical Chandra or XMM pointing observation. The new CC diagnostic has several advantages over its counterparts, e.g., it is free of the effects of the commonly seen centroid shift of the X-ray halo caused by merger event, and the corresponding calculation is straightforward, almost irrelevant to the complicated spectral analysis.

We present results on the formation of Pop III stars at redshift 7.6 from the Renaissance Simulations, a suite of extremely high-resolution and physics-rich radiation transport hydrodynamics cosmological adaptive-mesh refinement simulations of high redshift galaxy formation performed on the Blue Waters supercomputer. In a survey volume of about 220 comoving Mpc$^3$, we found 14 Pop III galaxies with recent star formation. The surprisingly late formation of Pop III stars is possible due to two factors: (i) the metal enrichment process is local and slow, leaving plenty of pristine gas to exist in the vast volume; and (ii) strong Lyman-Werner radiation from vigorous metal-enriched star formation in early galaxies suppresses Pop III formation in ("not so") small primordial halos with mass less than $\sim$ 3 $\times$ 10$^7$ M$_\odot$. We quantify the properties of these Pop III galaxies and their Pop III star formation environments. We look for analogues to the recently discovered luminous Ly $\alpha$ emitter CR7 (Sobral et al. 2015), which has been interpreted as a Pop III star cluster within or near a metal-enriched star forming galaxy. We find and discuss a system similar to this in some respects, however the Pop III star cluster is far less massive and luminous than CR7 is inferred to be.

Utilizing an archived Suzaku data acquired on 2007 December 25 for 46 ks, X-ray spectroscopic properties of the dipping and eclipsing low-mass X-ray binary EXO 0748$-$676 were studied. At an assumed distance of 7.1 kpc, the data gave a persistent unabsorbed luminosity of $3.4\times10^{36}$ erg cm$^{-2}$ s$^{-1}$ in 0.6 $-$ 55 keV. The source was in a relatively bright low/hard state, wherein the 0.6 $-$ 55 keV spectrum can be successfully explained by a "double-seed" Comptonization model, incorporating a common corona with an electron temperature of $\sim13$ keV. The seed photons are thought to be supplied from both the neutron star surface, and a cooler truncated disk. Compared to a sample of non-dipping low-mass X-ray binaries in the low/hard state, the spectrum is subject to stronger Comptonization, with a relatively larger Comptonizing $y$-parameter of $\sim1.4$ and a larger coronal optical depth of $\sim5$. This result, when attributed to the high inclination of EXO 0748$-$676, suggests that the Comptonizing corona may elongate along the disk plane, and give a longer path for the seed photons when viewed from edge-on inclinations.

Detection and mitigation of radio frequency interference (RFI) is the first and also the key step for data processing in radio observations, especially for ongoing low frequency radio experiments towards the detection of the cosmic dawn and epoch of reionization (EoR). In this paper we demonstrate the technique and efficiency of RFI identification and mitigation for the 21 Centimeter Array (21CMA), a radio interferometer dedicated to the statistical measurement of EoR. For terrestrial, man-made RFI, we concentrate mainly on a statistical approach by identifying and then excising non-Gaussian signatures, in the sense that the extremely weak cosmic signal is actually buried under thermal and therefore Gaussian noise. We also introduce the so-called visibility correlation coefficient instead of conventional visibility, which allows a further suppression of rapidly time-varying RFI. Finally, we briefly discuss removals of the sky RFI, the leakage of sidelobes from off-field strong radio sources with time-invariant power and a featureless spectrum. It turns out that state of the art technique should allow us to detect and mitigate RFI to a satisfactory level in present low frequency interferometer observations such as those acquired with the 21CMA, and the accuracy and efficiency can be greatly improved with the employment of low-cost, high-speed computing facilities for data acquisition and processing.

To study systematically the evolution on the angular extents of the galaxy, ICM, and dark matter components in galaxy clusters, we compiled the optical and X-ray properties of a sample of 340 clusters with redshifts $<0.5$, based on all the available data with the Sloan Digital Sky Survey (SDSS) and {\it Chandra}/{\it XMM-Newton}. For each cluster, the member galaxies were determined primarily with photometric redshift measurements. The radial ICM mass distribution, as well as the total gravitational mass distribution, were derived from a spatially-resolved spectral analysis of the X-ray data. When normalizing the radial profile of galaxy number to that of the ICM mass, the relative curve was found to depend significantly on the cluster redshift; it drops more steeply towards outside in lower redshift subsamples. The same evolution is found in the galaxy-to-total mass profile, while the ICM-to-total mass profile varies in an opposite way. We interpret that the galaxies, the ICM, and the dark matter components had similar angular distributions when a cluster was formed, while the galaxies travelling interior of the cluster have continuously fallen towards the center relative to the other components, and the ICM has slightly expanded relative to the dark matter although it suffers strong radiative loss. This cosmological galaxy infall, accompanied by an ICM expansion, can be explained by considering that the galaxies interact strongly with the ICM while they are moving through it. The interaction is considered to create a large energy flow of $10^{44-45}$ erg $\rm s^{-1}$ per cluster from the member galaxies to their environment, which is expected to continue over cosmological time scales.

We apply a number of statistical and machine learning techniques to classify and rank gamma-ray sources from the Third Fermi Large Area Telescope (LAT) Source Catalog (3FGL), according to their likelihood of falling into the two major classes of gamma-ray emitters: pulsars (PSR) or Active Galactic Nuclei (AGN). Using 1904 3FGL sources that have been identified/associated with AGN (1738) and PSR (166), we train (using 70% of our sample) and test (using 30%) our algorithms and find that the best overall accuracy (>96%) is obtained with the Random Forest (RF) technique, while using a logistic regression (LR) algorithm results in only marginally lower accuracy. We apply the same techniques on a sub-sample of 142 known gamma-ray pulsars to classify them into two major subcategories: young (YNG) and millisecond pulsars (MSP). Once more, the RF algorithm has the best overall accuracy (~90%), while a boosted LR analysis comes a close second. We apply our two best models (RF and LR) to the entire 3FGL catalog, providing predictions on the likely nature of {\it unassociated} sources, including the likely type of pulsar (YNG or MSP). We also use our predictions to shed light on the possible nature of some gamma-ray sources with known associations (e.g. binaries, SNR/PWN). Finally, we provide a list of plausible X-ray counterparts for some pulsar candidates, obtained using Swift, Chandra, and XMM. The results of our study will be of interest for both in-depth follow-up searches (e.g. pulsar) at various wavelengths, as well as for broader population studies.

The extent of nucleosynthesis in models of type I X-ray bursts and the associated impact on the energy released in these explosive events are sensitive to nuclear masses and reaction rates around the $^{64}$Ge waiting point. Using the well known mass of $^{64}$Ge, the recently measured $^{65}$As mass, and large-scale shell model calculations, we have determined new thermonuclear rates of the $^{64}$Ge($p$,$\gamma$)$^{65}$As and $^{65}$As($p$,$\gamma$)$^{66}$Se reactions with reliable uncertainties. The new reaction rates differ significantly from previously published rates. Using the new data we analyze the impact of the new rates and the remaining nuclear physics uncertainties on the $^{64}$Ge waiting point in a number of representative one-zone X-ray burst models. We find that in contrast to previous work, when all relevant uncertainties are considered, a strong $^{64}$Ge $rp$-process waiting point cannot be ruled out. The nuclear physics uncertainties strongly affect X-ray burst model predictions of the synthesis of $^{64}$Zn, the synthesis of nuclei beyond $A=64$, energy generation, and burst light curve. We also identify key nuclear uncertainties that need to be addressed to determine the role of the $^{64}$Ge waiting point in X-ray bursts. These include the remaining uncertainty in the $^{65}$As mass, the uncertainty of the $^{66}$Se mass, and the remaining uncertainty in the $^{65}$As($p$,$\gamma$)$^{66}$Se reaction rate, which mainly originates from uncertain resonance energies.

Gamma Ray Bursts (GRBs) are the strongest explosions in the universe which might be associated with creation of black holes. Magnetic field structure and burst dynamics may influence polarization of the emitted gamma-rays. Precise polarization detection can be an ultimate tool to unveil the true GRB mechanism. POLAR is a space-borne Compton scattering detector for precise measurements of the GRB polarization. It consists of a 40$\times$40 array of plastic scintillator bars read out by 25 multi-anode PMTs (MaPMTs). It is scheduled to be launched into space in 2016 onboard of the Chinese space laboratory TG2. We present a dedicated methodology for POLAR calibration and some calibration results based on the combined use of the laboratory radioactive sources and polarized X-ray beams from the European Synchrotron Radiation Facility. They include calibration of the energy response, computation of the energy conversion factor vs. high voltage as well as determination of the threshold values, crosstalk contributions and polarization modulation factors.

In order to investigate the spatial distribution of the ICM temperature in galaxy clusters in a quantitative way and probe the physics behind, we analyze the X-ray spectra of a sample of 50 galaxy clusters, which were observed with the Chandra ACIS instrument in the past 15 years, and measure the radial temperature profiles out to $0.45r_{500}$. We construct a physical model that takes into account the effects of gravitational heating, thermal history (such as radiative cooling, AGN feedback, and thermal conduction) and work done via gas compression, and use it to fit the observed temperature profiles by running Bayesian regressions. The results show that in all cases our model provides an acceptable fit at the 68% confidence level. To further validate this model we select nine clusters that have been observed with both Chandra (out to $\gtrsim 0.3r_{500}$) and Suzaku (out to $\gtrsim 1.5r_{500}$), fit their Chandra spectra with our model, and compare the extrapolation of the best-fits with the Suzaku measurements. We find that the model profiles agree with the Suzaku results very well in seven clusters. In the rest two clusters the difference between the model and observation is possibly caused by local thermal substructures. Our study also implies that for most of the clusters the assumption of hydrostatic equilibrium is safe out to at least $0.5r_{500}$, and the non-gravitational interactions between dark matter and its luminous counterpart is consistent with zero.

The electric current helicity density $\displaystyle \chi=\langle\epsilon_{ijk}b_i\frac{\partial b_k}{\partial x_j}\rangle$ contains six terms, where $b_i$ are components of the magnetic field. Due to the observational limitations, only four of the above six terms can be inferred from solar photospheric vector magnetograms. By comparing the results for simulation we distinguished the statistical difference of above six terms for isotropic and anisotropic cases. We estimated the relative degree of anisotropy for three typical active regions and found that it is of order 0.8 which means the assumption of local isotropy for the observable current helicity density terms is generally not satisfied for solar active regions. Upon studies of the statistical properties of the anisotropy of magnetic field of solar active regions with latitudes and with evolution in the solar cycle, we conclude that the consistency of that assumption of local homogeneity and isotropy requires further analysis in the light of our findings.

We use imaging from the Next Generation Virgo cluster Survey (NGVS) to present a comparative study of ultra-compact dwarf (UCD) galaxies associated with three prominent Virgo sub-clusters: those centered on the massive, red-sequence galaxies M87, M49 and M60. We show how UCDs can be selected with high completeness using a combination of half-light radius and location in color-color diagrams ($u^*iK_s$ or $u^*gz$). Although the central galaxies in each of these sub-clusters have nearly identical luminosities and stellar masses, we find large differences in the sizes of their UCD populations, with M87 containing ~3.5 and 7.8 times more UCDs than M49 and M60, respectively. The relative abundance of UCDs in the three regions scales in proportion to sub-cluster mass, as traced by X-ray gas mass, total gravitating mass, number of globular clusters, and number of nearby galaxies. We find that the UCDs are predominantly blue in color, with ~85% of the UCDs having colors similar to blue GCs and stellar nuclei of dwarf galaxies. We present evidence that UCDs surrounding M87 and M49 may follow a morphological sequence ordered by the prominence of their outer, low surface brightness envelope, ultimately merging with the sequence of nucleated low-mass galaxies, and that envelope prominence correlates with distance from either galaxy. Our analysis provides evidence that tidal stripping of nucleated galaxies is an important process in the formation of UCDs.

In this paper, we present the first results from the Renaissance Simulations, a suite of extremely high-resolution and physics-rich AMR calculations of high redshift galaxy formation performed on the Blue Waters supercomputer. These simulations contain hundreds of well-resolved galaxies at $z \sim 25-8$, and make several novel, testable predictions. Most critically, we show that the ultraviolet luminosity function of our simulated galaxies is consistent with observations of high-z galaxy populations at the bright end of the luminosity function (M$_{1600} \leq -17$), but at lower luminosities is essentially flat rather than rising steeply, as has been inferred by Schechter function fits to high-z observations, and has a clearly-defined lower limit in UV luminosity. This behavior of the luminosity function is due to two factors: (i) the strong dependence of the star formation rate on halo virial mass in our simulated galaxy population, with lower-mass halos having systematically lower star formation rates and thus lower UV luminosities; and (ii) the fact that halos with virial masses below $\simeq 2 \times 10^8$ M$_\odot$ do not universally contain stars, with the fraction of halos containing stars dropping to zero at $\simeq 7 \times 10^6$ M$_\odot$. Finally, we show that the brightest of our simulated galaxies may be visible to current and future ultra-deep space-based surveys, particularly if lensed regions are chosen for observation.

We present our prediction for the local 21 cm differential brightness temperature ($\delta T_{b}$) from a set of strongly clustered sources of Population III (Pop III) and II (Pop II) objects in the early Universe, by a numerical simulation of their formation and radiative feedback. These objects are located inside a highly biased environment, which is a rare, high-density peak ("Rarepeak") extending to $\sim7$ comoving Mpc. We study the impact of ultraviolet (UV) and X-ray photons on the intergalactic medium (IGM) and the resulting $\delta T_{b}$, when Pop III stars are assumed to emit X-ray photons by forming X-ray binaries very efficiently. We parameterize the rest-frame spectral energy distribution (SED) of X-ray photons, which regulates X-ray photon-trapping, IGM-heating, secondary Lyman-alpha pumping and the resulting morphology of $\delta T_{b}$. A combination of emission ($\delta T_{b}>0$) and absorption ($\delta T_{b}<0$) regions appears in varying amplitudes and angular scales. The boost of the signal by the high-density environment ($\delta\sim0.64$) and on a relatively large scale combine to make Rarepeak a discernible, spatially-extended ($\theta\sim10'$) object for 21 cm observation at $13\lesssim z\lesssim17$, which is found to be detectable as a single object by SKA with integration time of $\sim1000$ hours. Power spectrum analysis by some of the SKA precursors (LOFAR, MWA, PAPER) of such rare peaks is found difficult due to the rarity of these peaks, and the contribution only by these rare peaks to the total power spectrum remains subdominant compared to that by all astrophysical sources.

Galaxy clusters are unique laboratories to investigate turbulent fluid motions and large scale magnetic fields. Synchrotron radio halos at the center of merging galaxy clusters provide the most spectacular and direct evidence of the presence of relativistic particles and magnetic fields associated with the intracluster medium. The study of polarized emission from radio halos is extremely important to constrain the properties of intracluster magnetic fields and the physics of the acceleration and transport of the relativistic particles. However, detecting this polarized signal is a very hard task with the current radio facilities.We use cosmological magneto-hydrodynamical simulations to predict the expected polarized surface brightness of radio halos at 1.4 GHz. We compare these expectations with the sensitivity and the resolution reachable with the SKA1. This allows us to evaluate the potential for studying intracluster magnetic fields in the surveys planned for SKA1.

The turbulent motion within molecular clouds is a key factor controlling star formation. Turbulence supports molecular cloud cores from evolving to gravitational collapse and hence sets a lower bound on the size of molecular cloud cores in which star formation can occur. On the other hand, without a continuous external energy source maintaining the turbulence, such as in molecular clouds, the turbulence decays with an energy dissipation time comparable to the dynamic timescale of clouds, which could change the size limits obtained from Jean's criterion by assuming constant turbulence intensities. Here we adopt scaling relations of physical variables in decaying turbulence to analyze its specific effects on the formation of stars. We find that the decay of turbulence provides an additional approach for Jeans' criterion to be achieved, after which gravitational infall governs the motion of the cloud core. This epoch of turbulence decay is defined as cloud core relaxation. The existence of cloud core relaxation provides a more complete understanding in the competition between turbulence and gravity on the dynamics of molecular cloud cores and star formation.

Recently the number of main-sequence and subgiant stars exhibiting solar-like oscillations that are resolved into individual mode frequencies has increased dramatically. While only a few such data sets were available for detailed modeling just a decade ago, the Kepler mission has produced suitable observations for hundreds of new targets. This rapid expansion in observational capacity has been accompanied by a shift in analysis and modeling strategies to yield uniform sets of derived stellar properties more quickly and easily. We use previously published asteroseismic and spectroscopic data sets to provide a uniform analysis of 42 solar-type Kepler targets from the Asteroseismic Modeling Portal (AMP). We find that fitting the individual frequencies typically doubles the precision of the asteroseismic radius, mass and age compared to grid-based modeling of the global oscillation properties, and improves the precision of the radius and mass by about a factor of three over empirical scaling relations. We demonstrate the utility of the derived properties with several applications.

The first galaxies in the Universe are the building blocks of all observed galaxies. We present scaling relations for galaxies forming at redshifts $z \ge 15$ when reionization is just beginning. We utilize the ``Rarepeak' cosmological radiation hydrodynamics simulation that captures the complete star formation history in over 3,300 galaxies, starting with massive Population III stars that form in dark matter halos as small as ~$10^6 M_\odot$. We make various correlations between the bulk halo quantities, such as virial, gas, and stellar masses and metallicities and their respective accretion rates, quantifying a variety of properties of the first galaxies up to halo masses of $10^9 M_\odot$. Galaxy formation is not solely relegated to atomic cooling halos with virial temperatures greater than $10^4$ K, where we find a dichotomy in galaxy properties between halos above and below this critical mass scale. Halos below the atomic cooling limit have a stellar mass -- halo mass relationship $\log M_\star \simeq 3.5 + 1.3\log(M_{\rm vir} / 10^7 M_\odot)$. We find a non-monotonic relationship between metallicity and halo mass for the smallest galaxies. Their initial star formation events enrich the interstellar medium and subsequent star formation to a median of $10^{-2} Z_\odot$ and $10^{-1.5} Z_\odot$, respectively, in halos of total mass $10^7 M_\odot$ that is then diluted by metal-poor inflows, well beyond Population III pre-enrichment levels of $10^{-3.5} Z_\odot$. The scaling relations presented here can be employed in models of reionization, galaxy formation and chemical evolution in order to consider these galaxies forming prior to reionization.

Due to their long mean free path, X-rays are expected to have an important impact on cosmic reionization by heating and ionizing the IGM on large scales, especially after simulations have suggested that Population III stars may form in pairs at redshifts as high as 20 - 30. We use the Pop III distribution and evolution from a self-consistent cosmological radiation hydrodynamics simulation of the formation of the first galaxies and a simple Pop III X-ray binaries model to estimate their X-ray output in a high density region larger than 100 comoving (Mpc)$^3$. We then combine three different methods --- ray tracing, a one-zone model, and X-ray background modeling --- to investigate the X-ray propagation, intensity distribution, and long term effects on the IGM thermal and ionization state. The efficiency and morphology of photo-heating and photo-ionization are dependent on the photon energies. The sub-keV X-rays only impact the IGM near the sources while the keV photons contribute significantly to the X-ray background and heat and ionize the IGM smoothly. The X-rays just below 1 keV are most efficient in heating and ionizing the IGM. We find that the IGM might be heated to over 100 K by $z=10$ and the high density source region might reach 10$^4$ K, limited by atomic hydrogen cooling. This may be important for predicting the 21-cm neutral hydrogen signals. But, on the other hand, the free electrons from X-ray ionizations are not enough to contribute significantly to the optical depth of CMB to the Thomson scattering.

We investigate the luminosity and colour dependence of clustering of CMASS galaxies in the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey Tenth Data Release. The halo occupation distribution framework is adopted to model the projected two-point correlation function measurements on small and intermediate scales (from $0.02$ to $60\,h^{-1}{\rm {Mpc}}$) and to interpret the observed trends and infer the connection of galaxies to dark matter halos. We find that luminous red galaxies reside in massive halos of mass $M{\sim}10^{13}$--$10^{14}\,h^{-1}{\rm M_\odot}$ and more luminous galaxies are more clustered and hosted by more massive halos. The strong small-scale clustering requires a fraction of these galaxies to be satellites in massive halos, with the fraction at the level of 5--8 per cent and decreasing with luminosity. The characteristic mass of a halo hosting on average one satellite galaxy above a luminosity threshold is about a factor $8.7$ larger than that of a halo hosting a central galaxy above the same threshold. At a fixed luminosity, progressively redder galaxies are more strongly clustered on small scales, which can be explained by having a larger fraction of these galaxies in the form of satellites in massive halos. Our clustering measurements on scales below $0.4\,h^{-1}{\rm {Mpc}}$ allow us to study the small-scale spatial distribution of satellites inside halos. While the clustering of luminosity-threshold samples can be well described by a Navarro-Frenk-White (NFW) profile, that of the reddest galaxies prefers a steeper or more concentrated profile. Finally, we also use galaxy samples of constant number density at different redshifts to study the evolution of luminous galaxies, and find the clustering to be consistent with passive evolution in the redshift range of $0.5 \lesssim z \lesssim 0.6$.

We present the observation of a major solar eruption that is associated with fast sunspot rotation. The event includes a sigmoidal filament eruption, a coronal mass ejection, and a GOES X2.1 flare from NOAA active region 11283. The filament and some overlying arcades were partially rooted in a sunspot. The sunspot rotated at $\sim$10$^\circ$ per hour rate during a period of 6 hours prior to the eruption. In this period, the filament was found to rise gradually along with the sunspot rotation. Based on the HMI observation, for an area along the polarity inversion line underneath the filament, we found gradual pre-eruption decreases of both the mean strength of the photospheric horizontal field ($B_h$) and the mean inclination angle between the vector magnetic field and the local radial (or vertical) direction. These observations are consistent with the pre-eruption gradual rising of the filament-associated magnetic structure. In addition, according to the Non-Linear Force-Free-Field reconstruction of the coronal magnetic field, a pre-eruption magnetic flux rope structure is found to be in alignment with the filament, and a considerable amount of magnetic energy was transported to the corona during the period of sunspot rotation. Our study provides evidences that in this event sunspot rotation plays an important role in twisting, energizing, and destabilizing the coronal filament-flux rope system, and led to the eruption. We also propose that the pre-event evolution of $B_h$ may be used to discern the driving mechanism of eruptions.

The existence of supermassive black holes as early as z ~ 7 is one of the great unsolved problems in cosmological structure formation. One leading theory argues that they are born during catastrophic baryon collapse in z ~ 15 protogalaxies in strong Lyman-Werner UV backgrounds. Atomic line cooling in such galaxies fragments baryons into massive clumps that are thought to directly collapse to 10^4 - 10^5 solar-mass black holes. We have now discovered that some of these fragments can instead become supermassive stars that eventually explode as pair-instability supernovae with energies of ~ 10^55 erg, the most energetic explosions in the universe. We have calculated light curves and spectra for supermassive Pop III PI SNe with the Los Alamos RAGE and SPECTRUM codes. We find that they will be visible in NIR all-sky surveys by Euclid out to z ~ 10 - 15 and by WFIRST and WISH out to z ~ 15 - 20, perhaps revealing the birthplaces of the first quasars.

This paper describes the open-source code Enzo, which uses block-structured adaptive mesh refinement to provide high spatial and temporal resolution for modeling astrophysical fluid flows. The code is Cartesian, can be run in 1, 2, and 3 dimensions, and supports a wide variety of physics including hydrodynamics, ideal and non-ideal magnetohydrodynamics, N-body dynamics (and, more broadly, self-gravity of fluids and particles), primordial gas chemistry, optically-thin radiative cooling of primordial and metal-enriched plasmas (as well as some optically-thick cooling models), radiation transport, cosmological expansion, and models for star formation and feedback in a cosmological context. In addition to explaining the algorithms implemented, we present solutions for a wide range of test problems, demonstrate the code's parallel performance, and discuss the Enzo collaboration's code development methodology.

Recent simulations of Population III star formation have suggested that some fraction form in binary systems, in addition to having a characteristic mass of tens of solar masses. The deaths of metal-free stars result in the initial chemical enrichment of the universe and the production of the first stellar-mass black holes. Here we present a cosmological adaptive mesh refinement simulation of an overdense region that forms a few 10^9 Msun dark matter halos and over 13,000 Population III stars by redshift 15. We find that most halos do not form Population III stars until they reach Mvir ~ 10^7 Msun because this biased region is quickly enriched from both Population III and galaxies, which also produce high levels of ultraviolet radiation that suppress H2 formation. Nevertheless, Population III stars continue to form, albeit in more massive halos, at a rate of ~ 10^{-4} Msun yr^{-1} Mpc^{-3} at redshift 15. The most massive starless halo has a mass of 7 x 10^7 Msun, which could host massive black hole formation through the direct gaseous collapse scenario. We show that the multiplicity of the Population III remnants grows with halo mass above 10^8 Msun, culminating in 50 remnants located in 10^9 Msun halos on average. This has implications that high mass X-ray binaries and intermediate mass black holes that originate from metal-free stars may be abundant in high-redshift galaxies.

We propose a continuous wavelet transform based non-parametric foreground subtraction method for the detection of redshifted 21 cm signal from the epoch of reionization. This method works based on the assumption that the foreground spectra are smooth in frequency domain, while the 21 cm signal spectrum is full of saw-tooth-like structures, thus their characteristic scales are significantly different. We can distinguish them in the wavelet coefficient space easily and perform the foreground subtraction. Compared with the traditional spectral fitting based method, our method is more tolerant to complex foregrounds. Furthermore, we also find that when the instrument has uncorrected response error, our method can also work significantly better than the spectral fitting based method. Our method can obtain similar results with the Wp smoothing method, which is also a non-parametric method, but our method consumes much less computing time.

In this work, we present a high spatial resolution study of a fossil cluster, SDSS J0150-1005 ($z \simeq 0.364$), based on the imaging spectroscopic analysis of \chandra\ observation. The \chandra\ X-ray image shows a relax and symmetric morphology, which indicates SDSS J0150-1005 is a well-developed galaxy cluster with no sign of recent merger. Its global gas temperature is 5.73 $\pm$ 0.80 keV, and the virial mass is 6.23 $\pm$ 1.34 $\times 10^{14} $ M$_{\odot}$ according to the isothermal model. Compared with polytropic temperature model, the mass calculated based on isothermal model overestimates 49 $\pm$ 11%. The central gas entropy, $S_{0.1r_{200}} = 143.9 \pm 18.3$ keV cm$^2$, is significantly lower than average value of the normal galaxy clusters with similar temperatures. Our results indicate the early formation epoch of SDSS J0150-1005.

Synchrotron radio halos at the center of merging galaxy clusters provide the most spectacular and direct evidence of the presence of relativistic particles and magnetic fields associated with the intracluster medium. The study of polarized emission from radio halos has been shown to be extremely important to constrain the properties of intracluster magnetic fields. However, detecting this polarized signal is a very hard task with the current radio facilities.We investigate whether future radio observatories, such as the Square Kilometer Array (SKA) and its precursors and pathfinders, will be able to detect the polarized emission of radio halos in galaxy clusters.On the basis of cosmological magnetohydrodynamical simulations with initial magnetic fields injected by active galactic nuclei, we predict the expected radio halo polarized signal at 1.4 GHz. We compare these expectations with the limits of current radio facilities and explore the potential of the forthcoming radio interferometers to investigate intracluster magnetic fields through the detection of polarized emission from radio halos.The resolution and sensitivity values that are expected to be obtained in future sky surveys performed at 1.4 GHz using the SKA precursors and pathfinders (like APERTIF and ASKAP) are very promising for the detection of the polarized emission of the most powerful (L1.4GHz>10e25 Watt/Hz) radio halos. Furthermore, the JVLA have the potential to already detect polarized emission from strong radio halos, at a relatively low resolution.However, the possibility of detecting the polarized signal in fainter radio halos (L1.4GHz~10e24 Watt/Hz) at high resolution requires a sensitivity reachable only with SKA.

Based on optical and X-ray data for a sample of 34 relaxed rich clusters of galaxies with redshifts of 0.1-0.9, we studied relative spatial distributions of the two major baryon contents, the cluster galaxies and the hot plasmas. Using multi-band photometric data taken with the UH88 telescope, we determined the integrated (two dimensional) radial light profiles of member galaxies in each cluster using two independent approaches, i.e., the background subtraction and the color-magnitude filtering. The ICM mass profile of each cluster in our sample, also integrated in two dimensions, was derived from a spatially-resolved spectral analysis using XMM-Newton and Chandra data. Then, the radially-integrated light profile of each cluster was divided by its ICM mass profile, to obtain a profile of "galaxy light vs. ICM mass ratio". The ratio profiles over the central 0.65 R500 regions were found to steepen from the higher- to lower- redshift subsamples, meaning that the galaxies become more concentrated in the ICM sphere towards lower redshifts. The evolution is also seen in galaxy number vs. ICM mass ratio profiles. A range of systematic uncertainties in the galaxy light measurements, as well as many radius-/redshift- dependent biases to the galaxy vs. ICM profiles have been assessed, but none of them is significant against the observed evolution. Besides, the galaxy light vs. total mass ratio profiles also exhibit gradual concentration towards lower redshift. We interpret in the context that the galaxies, the ICM, and the dark matter components followed a similar spatial distribution in the early phase (z>0.5), while the galaxies have fallen towards the center relative to the others at a later phase.

Cosmological MHD simulations of galaxy cluster formation show a significant amplification of seed magnetic fields. We developed a novel method to decompose cluster magnetized turbulence into modes and showed that the fraction of the fast mode is fairly large, around 1/4 in terms of energy. This is larger than that was estimated before, which implies that cluster turbulence interacts with cosmic rays rather efficiently. We propose a framework to deal with electron and proton reacceleration in galaxy clusters that includes feedback on turbulence. In particular, we establish a new upper limit on proton and electron fluxes based on turbulence intensity. These findings, along with detailed modeling of reacceleration, will help to reconcile the observed giant radio haloes and the unobserved diffuse gamma-ray emission from these clusters.

We investigate to what extent the current helicity distribution observed in solar active regions is compatible with solar dynamo models. We use an advanced 2D mean-field dynamo model with dynamo action largely concentrated near the bottom of the convective zone, and dynamo saturation based on the evolution of the magnetic helicity and algebraic quenching. For comparison, we also studied a more basic 2D mean-field dynamo model with simple algebraic alpha quenching only. Using these numerical models we obtain butterfly diagrams for both the small-scale current helicity and the large-scale magnetic helicity, and compare them with the butterfly diagram for the current helicity in active regions obtained from observations. This comparison shows that the current helicity of active regions, as estimated by $-A \cdot B$ evaluated at the depth from which the active region arises, resembles the observational data much better than the small-scale current helicity calculated directly from the helicity evolution equation. Here $B$ and $A$ are respectively the dynamo generated mean magnetic field and its vector potential.

Aiming to correctly restore the redshifted 21 cm signals emitted by the neutral hydrogen during the cosmic reionization processes, we re-examine the separation approaches based on the quadratic polynomial fitting technique in frequency space to investigate whether they works satisfactorily with complex foreground, by quantitatively evaluate the quality of restored 21 cm signals in terms of sample statistics. We construct the foreground model to characterize both spatial and spectral substructures of the real sky, and use it to simulate the observed radio spectra. By comparing between different separation approaches through statistical analysis of restored 21 cm spectra and corresponding power spectra, as well as their constraints on the mean halo bias $b$ and average ionization fraction $x_e$ of the reionization processes, at $z=8$ and the noise level of 60 mK we find that, although the complex foreground can be well approximated with quadratic polynomial expansion, a significant part of Mpc-scale components of the 21 cm signals (75% for $\gtrsim 6h^{-1}$ Mpc scales and 34% for $\gtrsim 1h^{-1}$ Mpc scales) is lost because it tends to be mis-identified as part of the foreground when single-narrow-segment separation approach is applied. The best restoration of the 21 cm signals and the tightest determination of $b$ and $x_e$ can be obtained with the three-narrow-segment fitting technique as proposed in this paper. Similar results can be obtained at other redshifts.

Aims. A long timeline kinematic study of the archetypal CSO OQ 208 sheds light on the physical properties of the most compact radio sources. Methods. Archival data from the VLBA at 15 GHz over a time span of 13.6 yr are used to investigate the kinematics of the radio source. The flux density monitoring data obtained at the Michigan 26-meter radio telescope are also used as supplementary information. Results. At 8.4 and 15 GHz, the two lobes are resolved into two sub-components, identified as hotspots. A knotty jet is linked with the NE hotspot and traces back toward the geometric center. The core is too weak to be detected. Significant flux density variation is found in the primary hotspots with the maximum level of 62% (NE1) and 19% (SW1). The peak in the flux density of NE1 leads that of SW1 by approximately 5.00 yr, suggesting that the northeast lobe is advancing and the southwest lobe is receding. This light travel difference indicates a radial distance difference between the two hotspots of 1.53 pc, which indicates an inclination angle of about 80.8 degree between the radio jet and the line of sight. The angular separation rate between NE1 and SW1 is 0.027 mas/yr (or 0.133 c). The inner jet knot moves at 0.047 mas/yr (or 0.230 c), about 3.5 times the hotspot advancing speed. Conclusions. The large viewing angle and the modest jet speed suggest a mildly relativistic jet. The jet axis is close to the plane of the sky. The separation rate and the distance between the two primary hotspots result in a kinematic age of 255$\pm$17 yr, confirming that OQ 208 is indeed a young radio source. In addition to the hotspot advancing motions, sideways motions provide evidence that the lobes are obstructed by the external interstellar medium.

Non-thermal radio emission from cosmic ray electrons in the vicinity of merging galaxy clusters is an important tracer of cluster merger activity, and is the result of complex physical processes that involve magnetic fields, particle acceleration, gas dynamics, and radiation. In particular, objects known as radio relics are thought to be the result of shock-accelerated electrons that, when embedded in a magnetic field, emit synchrotron radiation in the radio wavelengths. In order to properly model this emission, we utilize the adaptive mesh refinement simulation of the magnetohydrodynamic evolution of a galaxy cluster from cosmological initial conditions. We locate shock fronts and apply models of cosmic ray electron acceleration that are then input into radio emission models. We have determined the thermodynamic properties of this radio-emitting plasma and constructed synthetic radio observations to compare to observed galaxy clusters. We find a significant dependence of the observed morphology and radio relic properties on the viewing angle of the cluster, raising concerns regarding the interpretation of observed radio features in clusters. We also find that a given shock should not be characterized by a single Mach number. We find that the bulk of the radio emission comes from gas with T>5x10^7, \rho~10^(-28)-10^(-27) g/cm^3, with magnetic field strengths of 0.1-1.0 \mu G and shock Mach numbers of M~3-6. We present an analysis of the radio spectral index which suggests that the spatial variation of the spectral index can mimic synchrotron aging. Finally, we examine the polarization fraction and position angle of the simulated radio features, and compare to observations.

Based on a 21.5 ks \chandra\ observation of A2556, we identify an edge on the surface brightness profile (SBP) at about 160$h_{71}^{-1}$ kpc northeast of the cluster center, and it corresponds to a shock front whose Mach number $\mathcal{M}$ is calculated to be $1.25_{-0.03}^{+0.02}$. No prominent substructure, such as sub-cluster, is found in either optical or X-ray band that can be associated with the edge, suggesting that the conventional super-sonic motion mechanism may not work in this case. As an alternative solution, we propose that the nonlinear steepening of acoustic wave, which is induced by the turbulence of the ICM at the core of the cluster, can be used to explain the origin of the shock front. Although nonlinear steepening weak shock is expected to occur frequently in clusters, why it is rarely observed still remains a question that requires further investigation, including both deeper X-ray observation and extensive theoretical studies.

The redshifted 21cm signal of neutral hydrogen from the epoch of reionization (EoR) is extremely weak and its first detection is therefore expected to be statistical with first-generation low-frequency radio interferometers. In this letter we propose a method to extract the angular power spectrum of EoR from the visibility correlation coefficients p_{ij}(u,v), instead of the visibilities V_{ij}(u,v) measured directly by radio interferometers in conventional algorithm. The visibility correlation coefficients are defined as p_{ij}(u,v)=V_{ij}(u,v)/\sqrt{|V_{ii}||V_{jj}|} by introducing the auto-correlation terms V_{ii} and V_{jj} such that the angular power spectrum C_{\ell} can be obtained through C_{\ell}=4pi^2T_0^2<|p_{ij}(u,v)|^2>, independently of the primary beams of antennas. This also removes partially the influence of receiver gains in the measurement of C_{\ell} because the amplitudes of the gains cancel each other out in the statistical average operation of <|p_{ij}(u,v)|^2>.We use the average system temperature T_0 as a calibrator of C_{\ell}, which is dominated by the Milky Way and extragalactic sources in our interested frequency range below 200 MHz. Finally we demonstrate the feasibility of the novel method using the simulated sky maps as targets and the 21 CentiMeter Array (21CMA) as interferometer.

The Swift/XRT light curve of the ultraluminous X-ray (ULX) source NGC 5408 X-1 was re-analyzed with two new numerical approaches, Weighted Wavelet $Z$-transform (WWZ) and CLEANest, that are different from previous studies. Both techniques detected a prominent periodicity with a time scale of $115.5\pm1.5$ days, in excellent agreement with the detection of the same periodicity first reported by Strohmayer (2009). Monte Carlo simulation was employed to test the statisiticak confidence of the 115.5-day periodicity, yielding a statistical significance of $> 99.98%$ (or $>3.8\sigma$). The robust detection of the 115.5-day quasi-periodic oscillations (QPOs), if it is due to the orbital motion of the binary, would infer a mass of a few thousand $M_\odot$ for the central black hole, implying an intermediate-mass black hole in NGC 5408 X-1.

Radio observations of galaxy clusters show that there are $\mu$G magnetic fields permeating the intra-cluster medium (ICM), but it is hard to accurately constrain the strength and structure of the magnetic fields without the help of advanced computer simulations. We present qualitative comparisons of synthetic VLA observations of simulated galaxy clusters to radio observations of Faraday Rotation Measure (RM) and radio halos. The cluster formation is modeled using adaptive mesh refinement (AMR) magneto-hydrodynamic (MHD) simulations with the assumption that the initial magnetic fields are injected into the ICM by active galactic nuclei (AGNs) at high redshift. In addition to simulated clusters in Xu et al. (2010, 2011), we present a new simulation with magnetic field injections from multiple AGNs. We find that the cluster with multiple injection sources is magnetized to a similar level as in previous simulations with a single AGN. The RM profiles from simulated clusters, both $|RM|$ and the dispersion of RM ($\sigma_{RM}$), are consistent at a first-order with the radial distribution from observations. The correlations between the $\sigma_{RM}$ and X-ray surface brightness from simulations are in a broad agreement with the observations, although there is an indication that the simulated clusters could be slightly over-dense and less magnetized with respect to those in the observed sample. In addition, the simulated radio halos agree with the observed correlations between the radio power versus the cluster X-ray luminosity and between the radio power versus the radio halo size. These studies show that the cluster wide magnetic fields that originate from AGNs and are then amplified by the ICM turbulence (Xu et al. 2010) match observations of magnetic fields in galaxy clusters.

Presence of spectral curvature due to synchrotron self-absorption (SSA) of extragalactic radio sources may break down the spectral smoothness feature-the premise that bright radio foreground can be successfully removed in the 21cm experiments of searching for the epoch of reionization (EOR).We present a quantitative estimate of the effect on the measurement of the angular power spectrum of the low-frequency sky,incorporating a phenomenological model,characterized by the fraction of radio sources with turnover frequencies in 100-1000MHz range and a broken power law for the spectral transition around turnover frequencies nu_m,into the simulated radio sources over a small sky area of 10x10 deg^2.We compare statistically the changes in their residual maps with/without inclusion of the SSA after the bright sources of S_150MHz>=100mJy are excised and the best-fitted polynomials in frequency domain on each pixel are further subtracted.It has been shown that the effect of SSA on the detection of EOR depends sensitively on the spectral profiles of radio sources around the turnover frequencies: A hard transition model described by the broken power law with the turnover of spectral index at nu_m would leave pronounced imprints on the residual background and cause serious confusion with the EOR signal.However,the spectral signatures on the angular power spectrum of extragalactic foreground generated by a soft transition model,in which the rise and fall power laws of spectral distribution around nu_m are connected through a smooth transition spanning >=200 MHz in characteristic width,can be fitted and consequently subtracted by employment of polynomials to an acceptable degree(delta T<1mK).As this latter scenario seems to be favored by both theoretical expectation and radio spectral observations,we conclude that the influence of SSA on the 21cm experiments is probably very minor.

We demonstrate that the current helicity observed in solar active regions traces the magnetic helicity of the large-scale dynamo generated field. We use an advanced 2D mean-field dynamo model with dynamo saturation based on the evolution of the magnetic helicity and algebraic quenching. For comparison, we also studied a more basic 2D mean-field dynamo model with simple algebraic alpha quenching only. Using these numerical models we obtained butterfly diagrams both for the small-scale current helicity and also for the large-scale magnetic helicity, and compared them with the butterfly diagram for the current helicity in active regions obtained from observations. This comparison shows that the current helicity of active regions, as estimated by $-{\bf A \cdot B}$ evaluated at the depth from which the active region arises, resembles the observational data much better than the small-scale current helicity calculated directly from the helicity evolution equation. Here ${\bf B}$ and ${\bf A}$ are respectively the dynamo generated mean magnetic field and its vector potential. A theoretical interpretation of these results is given.

Based on a detailed analysis of the high-quality Chandra, XMM-Newton, and Suzaku data of the X-ray bright cluster of galaxies Abell 1795, we report clear evidence for a two-phase intracluster medium (ICM) structure, which consists of a cool (with a temperature T = 2.0-2.2 keV) and a hot (T = 5.0-5.7 keV) component that coexist and dominate the X-ray emission at least in the central 80 kpc. A third weak emission component (T = 0.8 keV) is also detected within the innermost 144 kpc and is ascribed to a portion of inter-stellar medium (ISM) of the cD galaxy. Deprojected spectral analysis reveals flat radial temperature distributions for both the hot phase and cool phase components. These results are consistent with the ASCA measurements reported in Xu et al. (1998), and resemble the previous findings for the Centaurus cluster (e.g., Takahashi et al. 2009). By analyzing the emission measure ratio and gas metal abundance maps created from the Chandra data, we find that the cool phase component is more metal-enriched than the hot phase one in 50-100 kpc region, which agrees with that found in M87 (Simionescu et al. 2008). The coexistence of the cool phase and hot phase ICM cannot be realized by bubble uplifting from active galactic nuclei (AGN) alone. Instead, the two-phase ICM properties are better reconciled with a cD corona model (Makishima et al. 2001). (Abridged)

We examine the effects of self-gravity and magnetic fields on supersonic turbulence in isothermal molecular clouds with high resolution simulations and adaptive mesh refinement. These simulations use large root grids (512^3) to capture turbulence and four levels of refinement to capture high density, for an effective resolution of 8,196^3. Three Mach 9 simulations are performed, two super-Alfv\'enic and one trans-Alfv\'enic. We find that gravity splits the clouds into two populations, one low density turbulent state and one high density collapsing state. The low density state exhibits properties similar to non-self-gravitating in this regime, and we examine the effects of varied magnetic field strength on statistical properties: the density probability distribution function is approximately lognormal; velocity power spectral slopes decrease with field strength; alignment between velocity and magnetic field increases with field; the magnetic field probability distribution can be fit to a stretched exponential. The high density state is characterized by self-similar spheres; the density PDF is a power-law; collapse rate decreases with increasing mean field; density power spectra have positive slopes, P({\rho},k) \propto k; thermal-to-magnetic pressure ratios are unity for all simulations; dynamic-to-magnetic pressure ratios are larger than unity for all simulations; magnetic field distribution is a power-law. The high Alfv\'en Mach numbers in collapsing regions explain recent observations of magnetic influence decreasing with density. We also find that the high density state is found in filaments formed by converging flows, consistent with recent Herschel observations. Possible modifications to existing star formation theories are explored.

Theory and simulations suggest that magnetic fields from radio jets and lobes powered by their central super massive black holes can be an important source of magnetic fields in the galaxy clusters. This is paper II in a series of studies where we present self-consistent high-resolution adaptive mesh refinement cosmological magnetohydrodynamic (MHD) simulations that simultaneously follow the formation of a galaxy cluster and evolution of magnetic fields ejected by an active galactic nucleus (AGN). We studied 12 different galaxy clusters with virial masses ranging from 1 $\times$ 10$^{14}$ to 2 $\times$ 10$^{15}$ M$_{\odot}$. In this work we examine the effects of the mass and merger history on the final magnetic properties. We find that the evolution of magnetic fields is qualitatively similar to those of previous studies. In most clusters, the injected magnetic fields can be transported throughout the cluster and be further amplified by the intra-cluster medium (ICM) turbulence during the cluster formation process with hierarchical mergers, while the amplification history and the magnetic field distribution depend on the cluster formation and magnetism history. This can be very different for different clusters. The total magnetic energies in these clusters are between 4 $\times$ 10$^{57}$ and $10^{61}$ erg, which is mainly decided by the cluster mass, scaling approximately with the square of the total mass. Dynamically older relaxed clusters usually have more magnetic fields in their ICM. The dynamically very young clusters may be magnetized weakly since there is not enough time for magnetic fields to be amplified.

We employ simulations of supersonic super-Alfvenic turbulence decay as a benchmark test problem to assess and compare the performance of nine astrophysical MHD methods actively used to model star formation. The set of nine codes includes: ENZO, FLASH, KT-MHD, LL-MHD, PLUTO, PPML, RAMSES, STAGGER, and ZEUS. We present a comprehensive set of statistical measures designed to quantify the effects of numerical dissipation in these MHD solvers. We compare power spectra for basic fields to determine the effective spectral bandwidth of the methods and rank them based on their relative effective Reynolds numbers. We also compare numerical dissipation for solenoidal and dilatational velocity components to check for possible impacts of the numerics on small-scale density statistics. Finally, we discuss convergence of various characteristics for the turbulence decay test and impacts of various components of numerical schemes on the accuracy of solutions. We show that the best performing codes employ a consistently high order of accuracy for spatial reconstruction of the evolved fields, transverse gradient interpolation, conservation law update step, and Lorentz force computation. The best results are achieved with divergence-free evolution of the magnetic field using the constrained transport method, and using little to no explicit artificial viscosity. Codes which fall short in one or more of these areas are still useful, but they must compensate higher numerical dissipation with higher numerical resolution. This paper is the largest, most comprehensive MHD code comparison on an application-like test problem to date. We hope this work will help developers improve their numerical algorithms while helping users to make informed choices in picking optimal applications for their specific astrophysical problems.

We present a series of cosmological magnetohydrodynamic (MHD) simulations that simultaneously follow the formation of a galaxy cluster and evolution of magnetic fields ejected by an Active Galactic Nucleus (AGN). Specifically, we investigate the influence of both the epoch of AGN (z $\sim$ 3-0.5) and the AGN energy ($\sim$ 3 $\times$ 10$^{57}$ - 2 $\times$ 10$^{60}$ ergs)on the final magnetic field distribution in a relatively massive cluster (M$_{vir}$ $\sim$10$^{15}$ M$_\odot$). We find that as long as the AGN magnetic fields are ejected before the major mergers in the cluster formation history, magnetic fields can be transported throughout the cluster and can be further amplified by the intra-cluster medium (ICM) turbulence cause by hierarchical mergers during the cluster formation process. The total magnetic energy in the cluster can reach $\sim$ $10^{61}$ ergs, with micro Gauss fields distributed over $\sim$ Mpc scale. The amplification of the total magnetic energy by the ICM turbulence can be significant, up to $\sim$1000 times in some cases. Therefore even weak magnetic fields from AGNs can be used to magnetize the cluster to the observed level. The final magnetic energy in the ICM is determined by the ICM turbulent energy, with a weak dependence on the AGN injection energy. We discuss the properties of magnetic fields throughout the cluster and the synthetic Faraday rotation measure maps they produce. We also show that high spatial resolution over most of the magnetic regions of the cluster is very important to capture the small scale dynamo process and maintain the magnetic field structure in our simulations.

Properties of proton resonances in $^{18}$Ne have been investigated efficiently by utilizing a technique of proton resonant elastic scattering with a $^{17}$F radioactive ion (RI) beam and a thick proton target. A 4.22~MeV/nucleon $^{17}$F RI beam was produced via a projectile-fragmentation reaction, and subsequently separated by a Radioactive Ion Beam Line in Lanzhou ({\tt RIBLL}). Energy spectra of the recoiled protons were measured by two sets of $\Delta$E-E silicon telescope at center-of-mass scattering angles of $\theta_{c.m.}$$\approx$175${^\circ}$$\pm$5${^\circ}$, $\theta_{c.m.}$$\approx$152${^\circ}$$\pm$8${^\circ}$, respectively. Several proton resonances in $^{18}$Ne were observed, and their resonant parameters have been determined by an $R$-matrix analysis of the differential cross sections in combination with the previous results. The resonant parameters are related to the reaction-rate calculation of the stellar $^{14}$O($\alpha$,$p$)$^{17}$F reaction, which was thought to be the breakout reaction from the hot CNO cycles into the $rp$-process in x-ray bursters. Here, $J^\pi$=(3$^-$, 2$^-$) are tentatively assigned to the 6.15-MeV state which was thought the key 1$^-$ state previously. In addition, a doublet structure at 7.05 MeV are tentatively identified, and its contribution to the resonant reaction rate of $^{14}$O($\alpha$,$p$)$^{17}$F could be enhanced by at least factors of about 4$\sim$6 in comparison with the previous estimation involving only a singlet. The present calculated resonant rates are much larger than those previous values, and it may imply that this breakout reaction could play a crucial role under x-ray bursters conditions.

In this work, we present the mass and magnetic distributions found in a recent Adaptive Mesh Refinement (AMR) MHD simulation of supersonic, \sa, self gravitating turbulence. Powerlaw tails are found in both volume density and magnetic field probability density functions, with $P(\rho) \propto \rho^{-1.67}$ and $P(B)\propto B^{-2.74}$. A power law is also found between magnetic field strength and density, with $B\propto \rho^{0.48}$, throughout the collapsing gas. The mass distribution of gravitationally bound cores is shown to be in excellent agreement with recent observation of prestellar cores. The mass to flux distribution of cores is also found to be in excellent agreement with recent Zeeman splitting measurements.

In this work we simulate the $50-200$ MHz radio sky that is constrained in the field of view ($5^{\circ}$ radius) of the 21 Centimeter Array (21CMA), by carrying out Monte-Carlo simulations to model redshifted cosmological reionization signals and strong contaminating foregrounds, including emissions from our Galaxy, galaxy clusters, and extragalactic point sources. As an improvement of previous works, we consider in detail not only random variations of morphological and spectroscopic parameters within the ranges allowed by multi-band observations, but also evolution of radio halos in galaxy clusters, assuming that relativistic electrons are re-accelerated in the ICM in merger events and lose energy via both synchrotron emission and inverse Compton scattering with CMB photons. By introducing a new approach designed on the basis of independent component analysis (ICA) and wavelet detection algorithm, we prove that, with a cumulative observation of one month with the 21CMA array, about $80\%$ of galaxy clusters with central brightness temperatures of $> 10~{\rm K}$ at 65 MHz can be safely identified and separated from the overwhelmingly bright foreground. We find that the morphological and spectroscopic distortions are extremely small as compared to the input simulated clusters, and the reduced $\chi^2$ of brightness temperature profiles and spectra are controlled to be $\lesssim 0.5$ and $\lesssim 1.3$, respectively. These results robustly indicate that in the near future a sample of dozens of bright galaxy clusters will be disentangled from the foreground in 21CMA observations, the study of which will greatly improve our knowledge about cluster merger rates, electron acceleration mechanisms in cluster radio halos, and magnetic field in the ICM.

By analyzing Chandra X-ray data of a sample of 21 galaxy groups and 19 galaxy clusters, we find that in 31 sample systems there exists a significant central ($R^{<}_{\sim} 10h_{71}^{-1}$ kpc) gas entropy excess ($\Delta K_{0}$), which corresponds to $\simeq 0.1-0.5$ keV per gas particle, beyond the power-law model that best fits the radial entropy profile of outer regions. We also find a distinct correlation between the central entropy excess $\Delta K_{0}$ and $K$-band luminosity $L_{K}$ of the central dominating galaxies (CDGs), which is scaled as $\Delta K_{0} \propto L_{K}^{1.6\pm0.4}$, where $L_{K}$ is tightly associated with the mass of the supermassive black hole hosted in the CDG. In fact, if an effective mass-to-energy conversion-efficiency of 0.02 is assumed for the accretion process, the cumulative AGN feedback $E^{\rm AGN}_{\rm feedback} \simeq \eta M_{\rm BH}c^{2}$ yields an extra heating of $\simeq 0.5-17.0$ keV per particle, which is sufficient to explain the central entropy excess. In most cases the AGN contribution can compensate the radiative loss of the X-ray gas within the cooling radius ($\simeq 0.002-2.2$ keV per particle), and apparently exceeds the energy required to deviate the scaling relations from the self-similar predictions ($\simeq 0.2-1.0$ keV per particle). In contrast to the AGN feedback, the extra heating provided by supernova explosions accounts for $\simeq 0.01-0.08$ keV per particle in groups and is almost negligible in clusters. Therefore, the observed correlation between $\Delta K_{0}$ and $L_{K}$ can be considered as a direct evidence for AGN feedback in galaxy groups and clusters.

By analysing the Chandra and XMM-Newton archived data of the nearby galaxy cluster Abell 3158, which was reported to possess a relatively regular, relaxed morphology in the X-ray band in previous works, we identify a bow edge-shaped discontinuity in the X-ray surface brightness distribution at about $120h_{71}^{-1}$ kpc west of the X-ray peak. This feature is found to be associated with a massive, off-centre cool gas clump, and actually forms the west boundary of the cool clump. We find that the cool gas clump is moving at a subsonic velocity of ~700 km/s toward west on the sky plane. We exclude the possibility that this cool clump was formed by local inhomogeneous radiative cooling in the intra-cluster medium, due to the effectiveness of the thermal conduction on the time-scale of $\sim 0.3$ Gyr. Since no evidence for central AGN activity has been found in Abell 3158, and this cool clump bears many similarities to the off-centre cool gas clumps detected in other merging clusters in terms of their mass, size, location, and thermal properties (e.g. lower temperature and higher abundance as compared with the environment), we speculate that the cool clump in Abell 3158 was caused by a merger event, and is the remnant of the original central cool-core of the main cluster or the infalling sub-cluster. This idea is supported not only by the study of line-of-sight velocity distribution of the cluster member galaxies, but also by the study of gas entropy-temperature correlation. This example shows that the appearance of such massive, off-centre cool gas clumps can be used to diagnose the dynamical state of a cluster, especially when prominent shocks and cold fronts are absent.

By creating and analyzing the two dimensional gas temperature and abundance maps of the RGH 80 compact galaxy group with the high-quality Chandra data, we detect a high-abundance ($\simeq 0.7$ $Z_\odot$) arc, where the metal abundance is significantly higher than the surrounding regions by $\simeq 0.3$ Z_\odot$. This structure shows tight spatial correlations with the member galaxy PGC 046529, as well as with the arm-like feature identified on the X-ray image in the previous work of Randall et al. (2009). Since no apparent signature of AGN activity is found associated with PGC 046529 in multi-band observations, and the gas temperature, metallicity, and mass of the high-abundance arc resemble those of the ISM of typical early-type galaxies, we conclude that this high-abundance structure is the remnant of the ISM of PGC 046529, which was stripped out of the galaxy by ram pressure stripping due to the motion of PGC 046529 in RGH 80. This novel case shows that ram pressure stripping can work efficiently in the metal enrichment process in galaxy groups, as it can in galaxy clusters.

Previously unobservable mirror asymmetry of the solar magnetic field -- a key ingredient of the dynamo mechanism which is believed to drive the 11-year activity cycle -- has now been measured. This was achieved through systematic monitoring of solar active regions carried out for more than 20 years at observatories in Mees, Huairou, and Mitaka. In this paper we report on detailed analysis of vector magnetic field data, obtained at Huairou Solar Observing Station in China. Electric current helicity (the product of current and magnetic field component in the same direction) was estimated from the data and a latitude-time plot of solar helicity during the last two solar cycles has been produced. We find that like sunspots helicity patterns propagate equatorwards but unlike sunspot polarity helicity in each solar hemisphere does not change sign from cycle to cycle - confirming the theory. There are, however, two significant time-latitudinal domains in each cycle when the sign does briefly invert. Our findings shed new light on stellar and planetary dynamos and has yet to be included in the theory.

By analyzing the gas temperature maps created from the Chandra archive data, we reveal the prevailing existence of temperature substructures on ~100 kpc scales in the central regions of nine intermediate-redshift (z~0.1) galaxy clusters, which resemble those found in the Virgo and Coma Clusters. Each substructure contains a clump of hot plasma whose temperature is about 2-3 keV higher than the environment, corresponding to an excess thermal energy of ~1E58-1E60 erg per clump. Since if there were no significant non-gravitational heating sources, these substructures would have perished in 1E8-1E9 yrs due to thermal conduction and turbulent flows, whose velocity is found to range from about 200 to 400 km/s, we conclude that the substructures cannot be created and sustained by inhomogeneous radiative cooling. We also eliminate the possibilities that the temperature substructures are caused by supernova explosions, or by the non-thermal X-ray emission due to the inverse-Comptonization of the CMB photons. By calculating the rising time of AGN-induced buoyant bubbles, we speculate that the intermittent AGN outbursts (~ 1E60 erg per burst) may have played a crucial role in the forming of the high temperature substructures. Our results are supported by recent study of McNamara & Nulsen (2007), posing a tight observational constraint on future theoretical and numerical studies.

We present self-consistent cosmological magnetohydrodynamic (MHD) simulations that simultaneously follow the formation of a galaxy cluster and the magnetic field ejection by an active galactic nucleus (AGN). We find that the magnetic fields ejected by the AGNs, though initially distributed in relatively small volumes, can be transported throughout the cluster and be further amplified by the intra-cluster medium (ICM) turbulence during the cluster formation process. The ICM turbulence is shown to be generated and sustained by the frequent mergers of smaller halos. Furthermore, a cluster-wide dynamo process is shown to exist in the ICM and amplify the magnetic field energy and flux. The total magnetic energy in the cluster can reach $\sim$ $10^{61}$ ergs while micro Gauss ($\mu$G) fields can distribute over $\sim$ Mpc scales throughout the whole cluster. This finding shows that magnetic fields from AGNs, being further amplified by the ICM turbulence through small-scale dynamo processes, can be the origin of cluster-wide magnetic fields.

We present an observational search for the possible correlation between cosmic microwave background (CMB) polarization map and soft X-ray background (SXRB) based on the ROSAT All-sky Survey (RASS) archive and WMAP five-year observations. This is motivated by the fact that some of the CMB polarization may arise from the scattering of CMB photons due to the free electrons generated by X-ray heating sources in the epoch of reionization. Detection of such a correlation allows one to study the role of X-ray heating in the process of reionization. However, the cross angular power spectrum of the CMB polarization and SXRB maps constructed from ROSAT RASS and WMAP five-year maps is consistent with no correlation. We attribute this negative detection to both the extremely weak signals and large instrumental noises. While a future search is needed with high sensitivity instruments for both CMB polarization and soft X-ray photons, our current results can still be used as a useful constraint on the effect of X-ray heating in the epoch of reionization.

In this work, we present MHDEnzo, the extension of the cosmological code Enzo to include the effects magnetic fields through the ideal MHD approximation. We use a higher order Godunov Riemann solver for the computation of interface fluxes. We use two constrained transport methods to compute the electric field from those interface fluxes, which simultaneously advances the induction equation and maintains the divergence of the magnetic field. A third order divergence free reconstruction technique is used to interpolate the magnetic fields in the block structured AMR framework already extant in Enzo. This reconstruction also preserves the divergence of the magnetic field to machine precision. We use operator splitting to include gravity and cosmological expansion. We then present a series of cosmological and non cosmological tests problems to demonstrate the quality of solution resulting from this combination of solvers.

We present cosmological magnetohydrodynamic simulations of the formation of a galaxy cluster with magnetic energy feedback from an active galactic nuclei (AGN). We demonstrate that X-ray cavities can be produced by the magnetically dominated jet-lobe system that is supported by a central axial current. The cavities are magnetically dominated and their morphology is determined jointedly by the magnetic fields and the background cluster pressure profile. The expansion and motion of the cavities are driven initially by the Lorentz force of the magnetic fields, and the cavities only become buoyant at late stages ($> 500$ Myr). We find that up to $80%-90%$ of the injected magnetic energy goes into doing work against the hot cluster medium, heating it, and lifting it in the cluster potential.

We report the results of the first self-consistent three-dimensional adaptive mesh refinement magnetohydrodynamical simulations of Population III star formation including the Biermann Battery effect. We find that the Population III stars formed including this effect are both qualitatively and quantitatively similar to those from hydrodynamics-only (non-MHD) cosmological simulations. We observe peak magnetic fields of ~10^-9 G in the center of our star-forming halo at z ~ 17.55. The magnetic fields created by the Biermann Battery effect are predominantly formed early in the evolution of the primordial halo at low density and large spatial scales, and then grow through compression and by shear flows. The fields seen in this calculation are never large enough to be dynamically important (with beta >= 10^{15} at all times), and should be considered the minimum possible fields in existence during Population III star formation, and may be seed fields for the stellar dynamo or the magnetorotational instability at higher densities and smaller spatial scales.

Magnetic fields play an important role in almost all astrophysical phenomena including star formation. But due to the difficulty in analytic modeling and observation, magnetic fields are still poorly studied and numerical simulation has become a major tool. We have implemented a cosmological magnetohydrodynamics package for Enzo which is an AMR hydrodynamics code designed to simulate structure formation. We use the TVD solver developed by S. Li as the base solver. In addition, we employ the constrained transport(CT) algorithm as described by D. Balsara. For interpolation magnetic fields to fine grids we used a divergence free quadratic reconstruction, also described by Balsara. We present results from several test problems including MHD caustics, MHD pancake and galaxy cluster formation with magnetic fields. We also discuss possible applications of our AMR MHD code to first star research.

By analyzing the Chandra ACIS S3 data we studied the X-ray properties of the low-mass and high-mass X-ray binary populations in the nearby spiral galaxy NGC 5055. A total of 43 X-ray point sources were detected within the 2 effective radii, with 31 sources located on the disk and the rest 12 sources in the bulge. The resolved point sources dominate the total X-ray emission of the galaxy by accounting for about 80% of the total counts in 0.3--10 keV. By carrying out the spectral fittings we calculated the 0.3--10.0 keV luminosities of all the detected X-ray point sources and found that they span a wide range from a few 10^{37} erg s^{-1} to over 10^{39} erg s^{-1}. After compensating for the incompleteness at the low luminosity end, we find that the corrected XLF of the bulge population is well fitted with a broken power-law model with a break at 1.57^{+0.21}_{-0.20}\times 10^{38} erg s^{-1}, while the profile of the disk population's XLF agrees with a single power-law distribution with a slope of 0.93^{+0.07}_{-0.06}. The disk population is significantly richer at {}^{>}_{\sim}2\times10^{38} erg s^{-1} than the bulge population, inferring that the disk may have undergone relatively recent, strong starbursts that significantly increased the HMXB population, although ongoing starbursts are also observed in the nuclear region. Similar XLF profiles of the bulge and disk populations were found in M81. However, in most other spiral galaxies different patterns of the spatial variation of the XLF profiles from the bulge to the disk have been observed, indicating that the star formation and evolution history may be more complex than we have expected.

We present a Chandra study of the metal distribution in the X-ray bright compact group of galaxies HCG 62. We find that the diffuse X-ray emission is peaked at the core of the central galaxy NGC 4778, and is dominated by the contribution of the hot gas. The diffuse emission is roughly symmetric within \simeq0.25^{\prime}, which is straddled by double-sided X-ray cavities aligned in the northeast-southwest direction. By mapping the emission hardness ratio distributions and by performing the 2-dimensional spectral analysis, we identify a remarkable high-abundance arc region at about 2^{\prime} (33.6h_{70}^{-1} kpc) from the X-ray peak that spans over a vast region from south to northwest, a part of which roughly coinciding with the outer edge of the southwest X-ray cavity. The measured average abundance in this arc is higher than that in its neighboring regions by a factor of about 2, and the abundance ratios therein are nicely consistent with the dominance of the SN Ia yields. We estimate that the mass of iron contained in the arc is >3\times 10^{6}h_{70}^{-2.5} solar mass, which accounts for >3% of the iron synthesized in the galaxy. The high-abundance arc could have been formed by the AGN activities. However, it is also possible that the arc was formed in a recent merger as is implied by the recent optical kinematic study (Spavone et al. 2006), which implies that mergers may be as important as AGN activities in metal redistributions in early-type galaxies and their associated groups or clusters.

We present a study of the mass distributions in the bright E0 galaxy NGC 1407 and its associated group by analyzing the high quality Chandra and ROSAT X-ray spectroscopic data. In order to probe the stellar mass distribution we calculated the B-band mass-to-light ratio profile by comparing the observed line-strength indices and multi-color photometric data with different stellar synthesis model predictions. We find that the gas is single-phase with a temperature of ~0.7 keV within 1Re. Outside 1Re the gas temperature increases quickly outwards to >1 keV, indicating its group origin. We reveal that the X-ray surface brightness profile shows a central excess in the innermost region, and on both the total mass and dark matter profiles there is a flattened feature at about <1Re, which coincides with the gas temperature transition from the galactic level to the group level. The total mass and dark matter distributions within 0.85Re are cuspy and can be approximated by power-law profiles with indices of ~2, which are marginally consistent with the generalized NFW profiles with zeta=2. The mass in outer regions can be well fitted by a single NFW profile, with the derived concentration parameter c=18.6\pm1.5. We find that the NGC 1407 group has a baryon-dominated core, while the mass in the >1Re is dominated by dark matter. At the virial radius r200=572\pm118 kpc, the inferred mass and mass-to-light ratio are M200=2.20\pm0.42E13 solarM and Mvir/LB=311\pm60 solarM/solarL,respectively, showing that the NGC 1407 group is an extremely dark system even comparable to many clusters of galaxies. Since the obtained total mass is lower than those given in the earlier galaxy kinematic works, we speculate that NGC 1400 is not a virialized member in the group's gravitational potential well.

We present a Chandra ACIS study of the early-type galaxy NGC 4552. We detect 47 X-ray point sources, most of which are likely LMXBs, within 4 R_e. The brightest X-ray source coincides with the optical/UV/radio center of the galaxy, and shows variability on >1 hr timescales, indicating the possible existence of a LLAGN. The 46 off-center sources and the unresolved point sources contribute about 29% and 20% to the total luminosity of the galaxy, respectively. The corrected cumulative XLF of the off-center sources is best fit by a broken power-law model with a break at L_b=4.4+2.0-1.4 e38 ergs/s. We identified 210 GCs candidates in a HST WFPC2 optical image of the galaxy's central region. Of the 25 off-center LMXBs that fall within the WFPC2 FOV, 10 sources are coincident with a GC. Thus the fraction of the GCs hosting bright LMXBs and the fraction of the LMXBs associated with GCs are 4.8% and 40%, respectively. In the V and I bands, the GCs hosting bright LMXBs are typically 1-2 magnitudes brighter than the GCs with no detected LMXBs. There are about 1.9+-0.4 times as many LMXBs in the red, metal-rich GCs as there are in the blue, metal-poor ones. We find no obvious difference between the luminosity distributions of LMXBs in GCs and in field, but the cumulative spectrum of the LMXBs in GCs tends to be softer than that of the LMXBs in field. We detected 3 X-ray sources that have isotropic luminosities larger than 1e39 erg/s. The one located in the joint Chandra-HST field is found to be associated with a GC. By studying the ACIS spectra we infer that the this may be a candidate black hole system with a mass of 15-135 M_sun. One of the other sources with a luminosity brighter than 1e39 ergs/s reveals temporal variations in brightness on timescales greater than an hour.

We present the Chandra study of the intermediately distant cluster of galaxies MS 0839.9+2938. By performing both the projected and deprojected spectral analyses, we find that the gas temperature is approximately constant at about 4 keV in 130-444h_70^-1 kpc. In the inner regions, the gas temperature descends towards the center, reaching <~ 3 keV in the central 37h_70^-1 kpc. This infers that the lower and upper limits of the mass deposit rate are 9-34 M_sun yr^-1 and 96-126 M_sun yr^-1, respectively within 74h_70^-1 kpc where the gas is significantly colder. Along with the temperature drop, we detect a significant inward iron abundance increase from about 0.4 solar in the outer regions to about 1 solar within the central 37h_70^-1 kpc. Thus MS 0839.9+2938 is the cluster showing the most significant central iron excess at z>~ 0.2. We argue that most of the excess iron should have been contributed by SNe Ia. By utilizing the observed SN Ia rate and stellar mass loss rate, we estimate that the time needed to enrich the central region with excess iron is 6.4-7.9 Gyr, which is similar to those found for the nearby clusters. Coinciding with the optical extension of the cD galaxy (up to about 30h_70^-1 kpc), the observed X-ray surface brightness profile exhibits an excess beyond the distribution expected by either the beta model or the NFW model, and can be well fitted with an empirical two-beta model that leads to a relatively flatter mass profile in the innermost region.

We present joint X-ray and optical observations of the high redshift (z~0.83) lensing cluster CLJ0152.7-1357 made with the Chandra X-ray Observatory and the Keck telescope. We confirm the existence of significant substructure at both X-ray and optical wavelengths in the form of two distinct clumps, whose temperatures are 6.6(-1.5,+2.4) keV and 5.7(-1.6,+2.9) keV, respectively. The X-ray surface brightness profiles of the two clumps can be fitted by either a single beta-model or an NFW-like profile; the latter giving better fits to the central regions. We find that the X-ray derived mass of this cluster is in good agreement with independent lensing measurements. While its appearance indicates that the cluster has not reached a dynamical equilibrium state, its X-ray luminosity L_x, temperature T and dynamical mass M are consistent with the well-defined L_x-T and M-T relations for low-redshift galaxy clusters, which suggests that the dynamical properties of the clusters have remained almost unchanged since z~0.8.

In many clusters of galaxies, the cooling time at the core of the intracluster medium is much less than the age of the system, suggesting that the the gas should continually lose energy by radiation. Simple thermodynamic arguments then require that the expected "cooling flow" should exhibit a specific spectroscopic signature, characterized by a differential emission measure distribution that is inversely proportional to the cooling function of the plasma. That prediction can be quantitatively tested for the first time by the Reflection Grating Spectrometer (RGS) experiment on XMM-Newton, which provides high resolution X-ray spectra, even for moderately extended sources like clusters. We present RGS data on 14 separate cooling flow clusters, sampling a very wide range in mass deposition rate. Surprisingly, in all cases we find a systematic deficit of low temperature emission relative to the predictions of the cooling flow models. However, we do see evidence for cooling flow gas at temperatures just below the cluster background temperature, $T_0$, roughly down to $T_0/2$. These results are difficult to reconcile with most of the possible explanations for the cooling flow problem that have been proposed to date. We also present RGS data on the massive elliptical galaxy NGC 4636. In this case, we detect evidence for resonance emission line scattering of high oscillator strength Fe L-shell emission lines within the gaseous halo of the galaxy. The detection of that effect leads to very tight constraints on physical conditions within the halo. However, here again, the expected signature of a cooling flow is not detected, perhaps suggesting some fundamental uncertainty in our understanding of radiative cooling in low density cosmic plasmas.

We present the first high spectral resolution X-ray observation of the giant elliptical galaxy NGC 4636, obtained with the Reflection Grating Spectrometer on-board the XMM-Newton Observatory. The resulting spectrum contains a wealth of emission lines from various charge states of oxygen, neon, magnesium, and iron. Examination of the cross-dispersion profiles of several of these lines provides clear, unambiguous evidence of resonance scattering by the highest oscillator strength lines, as well as a weak temperature gradient in the inner regions of the interstellar medium. We invoke a sophisticated new Monte Carlo technique which allows us to properly account for these effects in performing quantitative fits to the spectrum. Our spectral fits are not subject to many of the systematics that have plagued earlier investigations. The derived metal abundances are higher than have been inferred from prior, lower spectral resolution observations of this source (Awaki et al. 1994), but are still incompatible with conventional chemical enrichment models of elliptical galaxies. In addition, our data are incompatible with standard cooling flow models for this system - our derived upper limit to the mass deposition rate is below the predicted value by a factor of 3--5.

Results of ASCA X-ray study of central regions of medium-richness clusters of galaxies are summarized, emphasizing differences between cD and non-cD clusters. The intra-cluster medium (ICM) is likely to consist of two (hot and cool) phases within $\sim 100$ kpc of a cD galaxy, where the ICM metallicity is also enhanced. In contrast, the ICM in non-cD clusters appears to be isothermal with little metallicity gradient right to the center. The gravitational potential exhibits a hierarchical nesting around cD galaxies, while a total mass-density profile with a central cusp is indicated for a non-cD cluster Abell~1060. The iron-mass-to-light ratio of the ICM decreases toward the center in both types of clusters, although it is radially constant in peripheral regions. The silicon-to-iron abundance ratio in the ICM increases with the cluster richness, but remains close to the solar ratio around cD galaxies. These overall results are interpreted without appealing to the popular cooling-flow hypothesis. Instead, an emphasis is put on the halo-in-halo structure formed around cD galaxies.

We present a new investigation of the mass-temperature (\MTrelation) relation of 22 nearby clusters based on the analysis of their ROSAT X-ray surface brightness profiles ($S_{X}(r)$) and their ASCA emission weighted temperatures. Two methods of the cluster mass estimations are employed and their results are compared: (1) the conventional $\beta$ model for gas distribution along with the isothermal and hydrostatic equilibrium assumptions, and (2) the NFW profile for dark matter distribution whose characteristic density and length are determined by the observed $S_{X}(r)$. These two models yield essentially the same goodness of fits for $S_{X}(r)$ and the similar \MTrelation relations, with the latter demonstrating a significant departure from the simple gravitational scaling of $M_{tot}\propto T_X^{3/2}$. It is also shown that the best-fit \MTrelation relations could be reconciled with the theoretical expectation if the low-temperature clusters ($T_X<3.5$ keV) are excluded from the list, which lends support to the scenario that the intra-cluster medium is preheated in the early phase of cluster formation. Together with the entropy-temperature distribution, the existence of a similarity break at $T_X=3-4$ keV in the dynamical scaling relations for galaxy clusters has been confirmed.

X-ray observations of three Shakhbazyan's Compact Groups of Galaxies, SCGG 202, SCGG 205, and SCGG 223, are presented for the first time. Extended X-ray emission was detected from SCGG 202 and SCGG 223 with 0.5-2 keV luminosity of 1E42 erg/s and 3E42 erg/s (for H_0 = 75 km/s/Mpc), respectively, while no significant emission was detected from SCGG 205 with an upper limit of 0.7E42 erg/s. The X-ray spectra of SCGG 202 and SCGG 223 can be described with thin-thermal plasma emission with temperature about 1 keV. X-ray properties of SCGG 202 and SCGG 223 are in good agreement with those of other known groups of galaxies, proving the physical nature of their grouping.

Modeling of the x-ray spectra of the Galactic superluminal jet sources GRS 1915+105 and GRO J1655-40 reveal a three-layered atmospheric structure in the inner region of their accretion disks. Above the cold and optically thick disk of a temperature 0.2-0.5 keV, there is a warm layer with a temperature of 1.0-1.5 keV and an optical depth around 10. Sometimes there is also a much hotter, optically thin corona above the warm layer, with a temperature of 100 keV or higher and an optical depth around unity. The structural similarity between the accretion disks and the solar atmosphere suggest that similar physical processes may be operating in these different systems.

X-ray spectral and imaging data from ASCA and ROSAT were used to measure the total mass profile in the central region of Abell 1060, a nearby and relatively poor cluster of galaxies. The ASCA X-ray spectra, after correcting for the spatial response of the X-ray telescope, show an isothermal distribution of the intra-cluster medium (ICM) within at least $\sim$ 12' (or $160h_{70}^{-1}$ kpc; $H_0 = 70 h_{70}$ km s$^{-1}$Mpc$^{-1}$) in radius of the cluster center. The azimuthally averaged surface brightness profile from the ROSAT PSPC exhibits a central excess above an isothermal $\beta$ model. The ring-sorted ASCA GIS spectra and the radial surface brightness distribution from the ROSAT PSPC were simultaneously utilized to constrain the gravitational potential profile. Some analytic models of the total mass density profile were examined. The ICM density profile was also specified by analytic forms. The ICM temperature distribution was constrained to satisfy the hydrostatic equilibrium, and to be consistent with the data. Then, the total mass distribution was found to be described better by the universal dark halo profile proposed by Navarro, Frenk, and White (1996;1997) than by a King-type model with a flat density core. A profile with a central cusp together with a logarithmic radial slope of $\sim 1.5$ was also consistent with the data. Discussions are made concerning the estimated dark matter distribution around the cluster center.

ASCA and ROSAT data of the Centaurus cluster were analyzed. A central excess in the radial brightness profile is found in the hard energy band up to 10keV. This requires a deeper central potential than a King-type one. A double-beta brightness distribution gives a good account of the data. A deprojected energy spectrum within a spherical region of radius 30 kpc at the cluster center indicates that the ICM cannot be isothermal there. Simultaneous fitting of the three (PSPC, GIS, and SIS) energy spectra extracted from the central 3' region showed that a two-temperature model that has temperatures of 1.4keV and 4.4keV and a common metallicity of 1.0 solar with the Galactic absorption gives the best description to the spectra. A cooling-flow spectral model did not give satisfactory fit. A 3-dimensional cluster model incorporating the two-phase picture, the double-beta brightness distribution, and the central metallicity increase reproduced the overall ASCA and ROSAT data successfully. The derived spatial distribution of the dark matter deviates from a King-type profile exhibiting a central excess, which is consistent with the Navvaro-Frenk-White model. The iron is revealed to be more widely spread than stars, which may have been caused by energetic supernovae and the motion of the cD galaxy. Since the derived characteristics of the cool phase including the temperature, angular extent, gas mass, and metallicity, are on a smooth extension of those of inter-stellar medium (ISM) of elliptical galaxies, the cool phase can be regarded as the ISM associated with the cD galaxy. The cool-phase X-ray emission is presumably sustained by energies dissipated by infalling mass to the cD galaxy in the course of cluster evolution.

We present a detailed analysis on the variability of the Fe K emission line in NGC 4051 using ASCA data. Through simple Gaussian line fits, we find not only obvious Fe K line variability with no significant difference in the X-ray continuum flux between two ASCA observations which were separated by $\sim$ 440 days, but also rapid variability of Fe K line on time scales $\sim 10^4$ s within the second observation. During the second observation, the line is strong (EW = 733$^{+206}_{-219}$ eV) and broad ($\sigma = 0.96^{+0.49}_{-0.35}$ keV) when the source is brightest, and become weaker (EW = 165$^{+87}_{-86}$ eV) and narrower ($\sigma<0.09$ keV) whilst the source is weakest. The equivalent width of Fe K line correlates positively with the continuum flux, which shows an opposite trend with another Seyfert 1 galaxy MCG --6-30-15.

Using the X-ray data from \ASCA, spectral and spatial properties of the intra-cluster medium (ICM) of the cD cluster Abell 1795 are studied, up to a radial distance of $\sim 12'$ ($\sim 1.3$ $h_{50}^{-1}$ kpc). The ICM temperature and abundance are spatially rather constant, although the cool emission component is reconfirmed in the central region. The azimuthally- averaged radial X-ray surface brightness profiles are very similar between soft (0.7--3 keV) and hard (3--10 keV) energy bands, and neither can be fitted with a single-$\beta$ model due to a strong data excess within $\sim5'$ of the cluster center. In contrast, double-$\beta$ models can successfully reproduce the overall brightness profiles both in the soft and hard energy bands, as well as that derived with the \ROSAT PSPC. Properties of the central excess brightness are very similar over the 0.2--10 keV energy range spanned by \ROSAT and \ASCA. Thus, the excess X-ray emission from the core region of this cluster is confirmed for the first time in hard X-rays above 3 keV. This indicates that the shape of the gravitational potential becomes deeper than the King-type one towards the cluster center. Radial profiles of the total gravitating matter, calculated using the double-$\beta$ model, reveal an excess mass of $\sim 3 \times 10^{13}~ M_{\odot}$ within $\sim 150 h^{-1}_{50}$ kpc of the cluster center. This suggests a hierarchy in the gravitational potential corresponding to the cD galaxy and the entire cluster.

We analyzed the ASCA X-ray data of 40 nearby clusters of galaxies, whose intracluster-medium temperature distributes in the range of 0.9--10 keV. We measured the Si and Fe abundances of the intracluster medium, spatially averaging over each cluster, but excluding the central $\sim 0.15 h_{50}^{-1}$ Mpc region in order to avoid any possible abundance gradients and complex temperature structures. The Fe abundances of these clusters are 0.2--0.3 solar, with only weak dependence on the temperature of the intracluster medium, hence on the cluster richness. In contrast, the Si abundance is observed to increase from 0.3 to 0.6--0.7 solar from the poorer to richer clusters. These results suggest that the supernovae of both type-Ia and type-II significantly contribute to the metal enrichment of the intracluster medium, with the relative contribution of type-II supernovae increasing towards richer clusters. We suggest a possibility that a considerable fraction of type-II supernova products escaped from poorer systems.

A large-scale gradient in the metal abundance has been detected with ASCA from an X-ray bright cluster of galaxies AWM7. The metal abundance shows a peak of 0.5 solar at the center and smoothly declines to <~ 0.2 solar at a radius of 500 kpc. The gas temperature is found to be constant at 3.8 keV. The radial distribution of iron can be fit with a beta-model with beta ~ 0.8 assuming the same core radius (115 kpc) as that of the intracluster medium. The metal distribution in AWM7 suggests that the gas injected from galaxies is not efficiently mixed in the cluster space and traces the distribution of galaxies.