Providing a possible connection between neutrino emission and gravitational-wave (GW) bursts is important to our understanding of the physical processes that occur when black holes or neutron stars merge. In the Daya Bay experiment, using data collected from December 2011 to August 2017, a search has been performed for electron-antineutrino signals coinciding with detected GW events, including GW150914, GW151012, GW151226, GW170104, GW170608, GW170814, and GW170817. We used three time windows of $\mathrm{\pm 10~s}$, $\mathrm{\pm 500~s}$, and $\mathrm{\pm 1000~s}$ relative to the occurrence of the GW events, and a neutrino energy range of 1.8 to 100 MeV to search for correlated neutrino candidates. The detected electron-antineutrino candidates are consistent with the expected background rates for all the three time windows. Assuming monochromatic spectra, we found upper limits (90% confidence level) on electron-antineutrino fluence of $(1.13~-~2.44) \times 10^{11}~\rm{cm^{-2}}$ at 5 MeV to $8.0 \times 10^{7}~\rm{cm^{-2}}$ at 100 MeV for the three time windows. Under the assumption of a Fermi-Dirac spectrum, the upper limits were found to be $(5.4~-~7.0)\times 10^{9}~\rm{cm^{-2}}$ for the three time windows.

The uncertainty of cosmological data generated from complex processes, such as observational Hubble parameter data (OHD) and the Type Ia supernova (SN Ia) data, cannot be accurately modeled by simple analytical probability distributions, e.g. Gaussian distribution. This necessitates the use of likelihood-free inference, which bypasses the direct calculation of likelihood. In this paper, we propose a new procedure to perform likelihood-free cosmological inference using two artificial neural networks (ANN), the Masked Autoregressive Flow (MAF) and denoising autoencoder (DAE). Our procedure is the first to use DAE to extract features from data in order to simplify the structure of MAF needed to estimate the posterior. Tested on simulated Hubble parameter data with a simple Gaussian likelihood, the procedure shows the capability of extracting feature from data and estimating posteriors without the need for tractable likelihood. We demonstrate that it can accurately approximate the real posterior, and achieves a performance comparable to the traditional MCMC method. We also discuss the application of the proposed procedure on OHD and Pantheon SN Ia data, and use them to constrain cosmological parameters from the non-flat $\Lambda$CDM model. With little prior knowledge added, we found constraints on $H_0, \Omega_m, \Omega_\Lambda$ similar to relevant work. In addition, this work is also the first to use Gaussian process in the procedure of OHD simulation.

We explore synergies between the space-based Wide-Field Infrared Survey Telescope (WFIRST) and the ground-based Rubin Observatory Legacy Survey of Space and Time (LSST). In particular, we consider a scenario where the currently envisioned survey strategy for WFIRST's High Latitude Survey (HLS), i.e., 2000 square degrees in four narrow photometric bands is altered in favor of a strategy that combines rapid coverage of the LSST area (to full LSST depth) in one band. We find that a 5-month WFIRST survey in the W-band can cover the full LSST survey area providing high-resolution imaging for >95% of the LSST Year 10 gold galaxy sample. We explore a second, more ambitious scenario where WFIRST spends 1.5 years covering the LSST area. For this second scenario we quantify the constraining power on dark energy equation of state parameters from a joint weak lensing and galaxy clustering analysis, and compare it to an LSST-only survey and to the Reference WFIRST HLS survey. Our survey simulations are based on the WFIRST exposure time calculator and redshift distributions from the CANDELS catalog. Our statistical uncertainties account for higher-order correlations of the density field, and we include a wide range of systematic effects, such as uncertainties in shape and redshift measurements, and modeling uncertainties of astrophysical systematics, such as galaxy bias, intrinsic galaxy alignment, and baryonic physics. Assuming the 5-month WFIRST wide scenario, we find a significant increase in constraining power for the joint LSST+WFIRST wide survey compared to LSST Y10 (FoM(Wwide)= 2.4 FoM(LSST)) and compared to LSST+WFIRST HLS (FoM(Wwide)= 5.5 FoM(HLS)).

We simulate the scientific performance of the Wide-Field Infrared Survey Telescope (WFIRST) High Latitude Survey (HLS) on dark energy and modified gravity. The 1.6 year HLS Reference survey is currently envisioned to image 2000 deg$^2$ in multiple bands to a depth of $\sim$26.5 in Y, J, H and to cover the same area with slit-less spectroscopy beyond z=3. The combination of deep, multi-band photometry and deep spectroscopy will allow scientists to measure the growth and geometry of the Universe through a variety of cosmological probes (e.g., weak lensing, galaxy clusters, galaxy clustering, BAO, Type Ia supernova) and, equally, it will allow an exquisite control of observational and astrophysical systematic effects. In this paper we explore multi-probe strategies that can be implemented given WFIRST's instrument capabilities. We model cosmological probes individually and jointly and account for correlated systematics and statistical uncertainties due to the higher order moments of the density field. We explore different levels of observational systematics for the WFIRST survey (photo-z and shear calibration) and ultimately run a joint likelihood analysis in N-dim parameter space. We find that the WFIRST reference survey alone (no external data sets) can achieve a standard dark energy FoM of >300 when including all probes. This assumes no information from external data sets and realistic assumptions for systematics. Our study of the HLS reference survey should be seen as part of a future community driven effort to simulate and optimize the science return of WFIRST.

Modifications of the matter power spectrum due to baryonic physics are one of the major theoretical uncertainties in cosmological weak lensing measurements. Developing robust mitigation schemes for this source of systematic uncertainty increases the robustness of cosmological constraints, and may increase their precision if they enable the use of information from smaller scales. Here we explore the performance of two mitigation schemes for baryonic effects in weak lensing cosmic shear: the PCA method and the halo-model approach in \textsc{HMcode}. We construct mock tomographic shear power spectra from four hydrodynamical simulations, and run simulated likelihood analyses with \textsc{CosmoLike} assuming LSST-like survey statistics. With an angular scale cut of $\ell_{\rm max}<2000$, both methods successfully remove the biases in cosmological parameters due to the various baryonic physics scenarios, with the PCA method causing less degradation in the parameter constraints than \textsc{HMcode}. For a more aggressive $\ell_{\rm max}=$5000, the PCA method performs well for all but one baryonic physics scenario, requiring additional training simulations to account for the extreme baryonic physics scenario of Illustris; \textsc{HMcode} exhibits tensions in the 2D posterior distributions of cosmological parameters due to lack of freedom in describing the power spectrum for $k > 10\ h^{-1}\mathrm{Mpc}$. We investigate variants of the PCA method and improve the bias mitigation through PCA by accounting for the noise properties in the data via Cholesky decomposition of the covariance matrix. Our improved PCA method allows us to retain more statistical constraining power while effectively mitigating baryonic uncertainties even for a broad range of baryonic physics scenarios.

Correlation between the variations of quasar absorption lines and the ionizing continuum have been recently confirmed in systematic studies. However, no convincing individual case is reported. We present a statistical analysis of the variable C iv broad absorption line (BAL) in the quasar SDSS J141007.74+541203.3, which have been observed with 44 epochs by the Sloan Digital Sky Survey Data Release 14. Grier et al. (2015) has recently concluded that the most likely cause of the variability of the BAL in SDSS J141007.74+541203.3 is a rapid response to changes in the incident ionizing continuum. In this paper, we confirm the anticorrelation between the equivalent width of BALs and the flux of the continuum based on the spectra of this quasar that show significant variations, which serve as another independent evidence for the view of Grier et al. (2015).

We investigate the electron-positron pair cascade taking place in the magnetosphere of a rapidly rotating black hole. Because of the spacetime frame dragging, the Goldreich-Julian charge density changes sign in the vicinity of the event horizon, which leads to an occurrence of a magnetic-field aligned electric field, in the same way as the pulsar outer-magnetospheric accelerator. In this lepton accelerator, electrons and positrons are accelerated in the opposite directions, to emit copious gamma-rays via the curvature and inverse-Compton processes. We examine a stationary pair cascade, and show that a stellar-mass black hole moving in a gaseous cloud can emit a detectable very-high-energy flux, provided that the black hole is extremely rotating and that the distance is less than about 1 kpc. We argue that the gamma-ray image will have a point-like morphology, and demonstrate that their gamma-ray spectra have a broad peak around 0.01-1 GeV and a sharp peak around 0.1 TeV, that the accelerators become most luminous when the mass accretion rate becomes about 0.01% of the Eddington rate, and that the predicted gamma-ray flux little changes in a wide range of magnetospheric currents. An implication of the stability of such a stationary gap is discussed.

Relativistic hydrodynamics is a powerful tool to simulate the evolution of the quark gluon plasma (QGP) in relativistic heavy ion collisions. Using 10000 initial and final profiles generated from 2+1-d relativistic hydrodynamics VISH2+1 with MC-Glauber initial conditions, we train a deep neural network based on stacked U-net, and use it to predict the final profiles associated with various initial conditions, including MC-Glauber, MC-KLN and AMPT and TRENTo. A comparison with the VISH2+1 results shows that the network predictions can nicely capture the magnitude and inhomogeneous structures of the final profiles, and nicely describe the related eccentricity distributions $P(\varepsilon_n)$ (n=2, 3, 4). These results indicate that deep learning technique can capture the main features of the non-linear evolution of hydrodynamics, showing its potential to largely accelerate the event-by-event simulations of relativistic hydrodynamics.

Using large samples containing nearly $2300$ active galaxies of low radio luminosity (1.4 GHz luminosity between $2\times 10^{23}$ and $3\times 10^{25}$ W/Hz, essentially low-excitation radio galaxies) at $z\lesssim 0.3$, we present a self-contained analysis of the dependence of the nuclear radio activity on both intrinsic and extrinsic properties of galaxies, with the goal of identifying the best predictors of the nuclear radio activity. While confirming the established result that stellar mass must play a key role in the triggering of radio activities, we point out that for central, most massive galaxies, the radio activity also shows a strong dependence on halo mass, which is unlikely due to enhanced interaction rates in denser regions in massive, cluster-scale halos. We thus further investigate the effects of various properties of the intracluster medium (ICM) in massive clusters on the radio activities, employing two standard statistical tools, Principle Component Analysis and Logistic Regression. It is found that ICM entropy, local cooling time, and pressure are the most effective in predicting the radio activity, pointing to the accretion of gas cooling out of a hot atmosphere to be the likely origin in triggering such activities in galaxies residing in massive dark matter halos. Our analysis framework enables us to logically discern the mechanisms responsible for the radio activity separately for central and satellite galaxies.

We study the orientations of satellite galaxies in redMaPPer clusters constructed from the Sloan Digital Sky Survey at $0.1<z<0.35$ to determine whether there is any preferential tendency for satellites to point radially toward cluster centers. We analyze the satellite alignment (SA) signal based on three shape measurement methods (re-Gaussianization, de Vaucouleurs, and isophotal shapes), which trace galaxy light profiles at different radii. The measured SA signal depends on these shape measurement methods. We detect the strongest SA signal in isophotal shapes, followed by de Vaucouleurs shapes. While no net SA signal is detected using re-Gaussianization shapes across the entire sample, the observed SA signal reaches a statistically significant level when limiting to a subsample of higher luminosity satellites. We further investigate the impact of noise, systematics, and real physical isophotal twisting effects in the comparison between the SA signal detected via different shape measurement methods. Unlike previous studies, which only consider the dependence of SA on a few parameters, here we explore a total of 17 galaxy and cluster properties, using a statistical model averaging technique to naturally account for parameter correlations and identify significant SA predictors. We find that the measured SA signal is strongest for satellites with the following characteristics: higher luminosity, smaller distance to the cluster center, rounder in shape, higher bulge fraction, and distributed preferentially along the major axis directions of their centrals. Finally, we provide physical explanations for the identified dependences, and discuss the connection to theories of SA.

Seismic data will be a vital geophysical constraint on internal structure of Europa if we land instruments on the surface. Quantifying expected seismic activity on Europa both in terms of large, recognizable signals and ambient background noise is important for understanding dynamics of the moon, as well as interpretation of potential future data. Seismic energy sources will likely include cracking in the ice shell and turbulent motion in the oceans. We define a range of models of seismic activity in Europa's ice shell by assuming each model follows a Gutenberg-Richter relationship with varying parameters. A range of cumulative seismic moment release between $10^{16}$ and $10^{18}$ Nm/yr is defined by scaling tidal dissipation energy to tectonic events on the Earth's moon. Random catalogs are generated and used to create synthetic continuous noise records through numerical wave propagation in thermodynamically self-consistent models of the interior structure of Europa. Spectral characteristics of the noise are calculated by determining probabilistic power spectral densities of the synthetic records. While the range of seismicity models predicts noise levels that vary by 80 dB, we show that most noise estimates are below the self-noise floor of high-frequency geophones, but may be recorded by more sensitive instruments. The largest expected signals exceed background noise by $\sim$50 dB. Noise records may allow for constraints on interior structure through autocorrelation. Models of seismic noise generated by pressure variations at the base of the ice shell due to turbulent motions in the subsurface ocean may also generate observable seismic noise.

Ice-covered ocean worlds possess diverse energy sources and associated mechanisms that are capable of driving significant seismic activity, but to date no measurements of their seismic activity have been obtained. Such investigations could probe their transport properties and radial structures, with possibilities for locating and characterizing trapped liquids that may host life and yielding critical constraints on redox fluxes, and thus on habitability. Modeling efforts have examined seismic sources from tectonic fracturing and impacts. Here, we describe other possible seismic sources, their associations with science questions constraining habitability, and the feasibility of implementing such investigations. We argue, by analogy with the Moon, that detectable seismic activity on tidally flexed ocean worlds should occur frequently. Their ices fracture more easily than rocks, and dissipate more tidal energy than the <1 GW of the Moon and Mars. Icy ocean worlds also should create less thermal noise for a due to their greater distance and consequently smaller diurnal temperature variations. They also lack substantial atmospheres (except in the case of Titan) that would create additional noise. Thus, seismic experiments could be less complex and less susceptible to noise than prior or planned planetary seismology investigations of the Moon or Mars.

In this Letter, we report the discovery of a strong correlation between the variability of narrow absorption lines (NALs) and the ionizing continuum from a two-epoch spectra sample of 40 quasars containing 52 variable C iv {\lambda}{\lambda}1548; 1551 absorption doublets. According to the concordance index, this sample is classified into two subsamples. Subsample I shows an anti-correlation between the variations of absorption lines and the continuum, while Subsample II exhibits a positive correlation. These results imply that these variable C iv {\lambda}{\lambda}1548; 1551 absorption doublets are intrinsic to the corresponding quasars and that their variations are caused primarily by the fluctuations of the ionizing continuum. Based on our analysis, we propose that there might be two kinds of absorption gas: one that is very sensitive to the continuum variations, the another that is not. In addition, we suggest that in many cases the emergence or disappearance of NALs is caused by fluctuations of the ionizing continuum.

In this Letter, we report the discovery of a strong correlation between the variability of narrow absorption lines (NALs) and the ionizing continuum from a two-epoch spectra sample of 40 quasars containing 52 variable C iv {\lambda}{\lambda}1548; 1551 absorption doublets. According to the concordance index, this sample is classified into two subsamples. Subsample I shows an anti-correlation between the variations of absorption lines and the continuum, while Subsample II exhibits a positive correlation. These results imply that these variable C iv {\lambda}{\lambda}1548; 1551 absorption doublets are intrinsic to the corresponding quasars and that their variations are caused primarily by the fluctuations of the ionizing continuum. Based on our analysis, we propose that there might be two kinds of absorption gas: one that is very sensitive to the continuum variations, the another that is not. In addition, we suggest that in many cases the emergence or disappearance of NALs is caused by fluctuations of the ionizing continuum.

We analyze sunspots rotation and magnetic transients in NOAA AR 11429 during two X-class (X5.4 and X1.3) flares using the data from the Helioseismic and Magnetic Imager on board the \emph{Solar Dynamics Observatory}. A large leading sunspot with positive magnetic polarity rotated counterclockwise. As expected, the rotation was significantly affected by the two flares. The magnetic transients induced by the flares were clearly evident in the sunspots with negative polarity. They were moving across the sunspots with speed of order $3-7\ \rm km \ s^{-1}$. Furthermore, the trend of magnetic flux evolution of these sunspots exhibited changes associated with the flares. These results may shed light on the understanding of the evolution of sunspots.

Geophysical measurements can reveal the structure of icy ocean worlds and cycling of volatiles. The associated density, temperature, sound speed, and electrical conductivity of such worlds thus characterizes their habitability. To explore the variability and correlation of these parameters, and to provide tools for planning and data analyses, we develop 1-D calculations of internal structure, which use available constraints on the thermodynamics of aqueous MgSO$_4$, NaCl (as seawater), and NH$_3$, water ices, and silicate content. Limits in available thermodynamic data narrow the parameter space that can be explored: insufficient coverage in pressure, temperature, and composition for end-member salinities of MgSO$_4$ and NaCl, and for relevant water ices; and a dearth of suitable data for aqueous mixtures of Na-Mg-Cl-SO$_4$-NH$_3$. For Europa, ocean compositions that are oxidized and dominated by MgSO$_4$, vs reduced (NaCl), illustrate these gaps, but also show the potential for diagnostic and measurable combinations of geophysical parameters. The low-density rocky core of Enceladus may comprise hydrated minerals, or anydrous minerals with high porosity comparable to Earth's upper mantle. Titan's ocean must be dense, but not necessarily saline, as previously noted, and may have little or no high-pressure ice at its base. Ganymede's silicious interior is deepest among all known ocean worlds, and may contain multiple phases of high-pressure ice, which will become buoyant if the ocean is sufficiently salty. Callisto's likely near-eutectic ocean cannot be adequately modeled using available data. Callisto may also lack high-pressure ices, but this cannot be confirmed due to uncertainty in its moment of inertia.

The shapes of cluster central galaxies are not randomly oriented, but rather exhibit coherent alignments with the shapes of their parent clusters as well as with the surrounding large-scale structures. In this work, we aim to identify the galaxy and cluster quantities that most strongly predict the central galaxy alignment phenomenon among a large parameter space with a sample of 8237 clusters and 94817 members within 0.1<z<0.35, based on the redMaPPer cluster catalog constructed from the Sloan Digital Sky Survey. We first quantify the alignment between the projected central galaxy shapes and the distribution of member satellites, to understand what central galaxy and cluster properties most strongly correlate with these alignments. Next, we investigate the angular segregation of satellites with respect to their central galaxy major axis directions, to identify the satellite properties that most strongly predict their angular segregation. We find that central galaxies are more aligned with their member galaxy distributions in clusters that are more elongated and have higher richness, and for central galaxies with larger physical size, higher luminosity and centering probability, and redder color. Satellites with redder color, higher luminosity, located closer to the central galaxy, and with smaller ellipticity show a stronger angular segregation toward their central galaxy major axes. Finally, we provide physical explanations for some of the identified correlations, and discuss the connection to theories of central galaxy alignments, the impact of primordial alignments with tidal fields, and the importance of anisotropic accretion.

Screened modified gravity (SMG) is a kind of scalar-tensor theories with screening mechanisms, which can generate screening effect to suppress the fifth force in high density environments and pass the solar system tests. Meanwhile, the potential of scalar field in the theories can drive the acceleration of the late universe. In this paper, we calculate the parameterized post-Newtonian (PPN) parameters $\gamma$ and $\beta$, the effective gravitational constant $G_{\rm eff}$ and the effective cosmological constant $\Lambda$ for SMG with a general potential $V$ and coupling function $A$. The dependence of these parameters on the model parameters of SMG and/or the physical properties of the source object are clearly presented. As an application of these results, we focus on three specific theories of SMG (chameleon, symmetron and dilaton models). Using the formulae to calculate their PPN parameters and cosmological constant, we derive the constraints on the model parameters by combining the observations on solar system and cosmological scales.

The fact that the clustering of dark matter halos depends not only on their mass, but also the formation epoch, is a prominent, albeit subtle, feature of the cold dark matter structure formation theory, and is known as assembly bias. At low mass scales ($\sim 10^{12}\,h^{-1}M_\odot$), early-forming halos are predicted to be more strongly clustered than the late-forming ones. In this study we aim to robustly detect the signature of assembly bias observationally, making use of formation time indicators of central galaxies in low mass halos as a proxy for the halo formation history. Weak gravitational lensing is employed to ensure our early- and late-forming halo samples have similar masses, and are free of contamination of satellites from more massive halos. For the two formation time indicators used (resolved star formation history and current specific star formation rate), we do not find convincing evidence of assembly bias. For a pair of early- and late-forming galaxy samples with mean mass $M_{200c} \approx 9\times10^{11}\,h^{-1}M_\odot$, the relative bias is $1.00\pm 0.12$. We attribute the lack of detection to the possibilities that either the current measurements of these indicators are too noisy, or they do not correlate well with the halo formation history. Alternative proxies for the halo formation history that should perform better are suggested for future studies.

Understanding the formation history of brightest cluster galaxies is an important topic in galaxy formation. Utilizing the Planck Sunyaev-Zel'dovich cluster catalog, and applying the Ansatz that the most massive halos at one redshift remain among the most massive ones at a slightly later cosmic epoch, we have constructed cluster samples at redshift z~0.4 and z~0.2 that can be statistically regarded as progenitor-descendant pairs. This allows us to study the stellar mass assembly history of BCGs in these massive clusters at late times, finding the degree of growth between the two epochs is likely at only few percent level, which is far lower compared to the prediction from a state-of-the-art semi-analytic galaxy formation model.

The China Dark Matter Experiment reports results on light WIMP dark matter searches at the China Jinping Underground Laboratory with a germanium detector array with a total mass of 20 g. The physics threshold achieved is 177 eVee ("ee" represents electron equivalent energy) at 50% signal efficiency. With 0.784 kg-days of data, exclusion region on spin-independent coupling with the nucleon is derived, improving over our earlier bounds at WIMP mass less than 4.6 GeV.

On 13 December 2012, Chang'e-2 conducted a successful flyby of the near-Earth asteroid 4179 Toutatis at a closest distance of 770 $\pm$ 120 meters from the asteroid's surface. The highest-resolution image, with a resolution of better than 3 meters, reveals new discoveries on the asteroid, e.g., a giant basin at the big end, a sharply perpendicular silhouette near the neck region, and direct evidence of boulders and regolith, which suggests that Toutatis may bear a rubble-pile structure. Toutatis' maximum physical length and width are (4.75 $\times$ 1.95 km) $\pm$10$\%$, respectively, and the direction of the +$z$ axis is estimated to be (250$\pm$5$^\circ$, 63$\pm$5$^\circ$) with respect to the J2000 ecliptic coordinate system. The bifurcated configuration is indicative of a contact binary origin for Toutatis, which is composed of two lobes (head and body). Chang'e-2 observations have significantly improved our understanding of the characteristics, formation, and evolution of asteroids in general.

The China Dark matter Experiment collaboration reports the first experimental limit on WIMP dark matter from 14.6 kg-day of data taken with a 994 g p-type point-contact germanium detector at the China Jinping underground Laboratory where the rock overburden is more than 2400 m. The energy threshold achieved was 400 eVee. According to the 14.6 kg-day live data, we placed the limit of N= 1.75 * 10^{-40} cm^{2} at 90% confidence level on the spin-independent cross-section at WIMP mass of 7 GeV before differentiating bulk signals from the surface backgrounds.

The transition between the proto-star, Class I, and the pre-main sequence star, Class II, phases is still one of the most uncertain, and important, stages in the knowledge of the process of formation of an individual star, because it is the stage that determines the final mass of the star. We observed the YSO "Holoea", associated with IRAS 05327+3404, which was classified as an object in transition between the Class I and Class II phases with several unusual properties, and appears to be surrounded by large amounts of circumstellar material. We used the SMA and BIMA telescopes at millimeter and sub-millimeter wavelengths to observe the dust continuum emission and the CO (1-0) and (2-1), HCO+ (1-0) and (3-2), and HCN (1-0) transitions in the region around IRAS 05327+3404. We detected two continuum emission peaks at 1.1-mm: SMM 1, the sub-mm counterpart of IRAS 05327+3404, and SMM 2, ~6 arcsec to the West. The emissions of the three molecules show marked differences. The CO emission near the systemic velocity is filtered out by the telescopes, and CO mostly traces the high-velocity gas. The HCO+ and HCN emissions are more centrally concentrated around the central parts of the region, and show several intensity peaks coincident with the sub-mm continuum peaks. We identify two main molecular outflows: a bipolar outflow in an E-W direction that would be powered by SMM 1 and another one in a NE direction, which we associate with SMM 2. We propose that the SMM sources are probably Class I objects, with SMM 1 in an earlier evolutionary stage.

The CDEX Collaboration has been established for direct detection of light dark matter particles, using ultra-low energy threshold p-type point-contact germanium detectors, in China JinPing underground Laboratory (CJPL). The first 1 kg point-contact germanium detector with a sub-keV energy threshold has been tested in a passive shielding system located in CJPL. The outputs from both the point-contact p+ electrode and the outside n+ electrode make it possible to scan the lower energy range of less than 1 keV and at the same time to detect the higher energy range up to 3 MeV. The outputs from both p+ and n+ electrode may also provide a more powerful method for signal discrimination for dark matter experiment. Some key parameters, including energy resolution, dead time, decay times of internal X-rays, and system stability, have been tested and measured. The results show that the 1 kg point-contact germanium detector, together with its shielding system and electronics, can run smoothly with good performances. This detector system will be deployed for dark matter search experiments.

Weakly Interacting Massive Particles (WIMPs) are the candidates of dark matter in our universe. Up to now any direct interaction of WIMP with nuclei has not been observed yet. The exclusion limits of the spin-independent cross section of WIMP-nucleon which have been experimentally obtained is about 10^{-7}pb at high mass region and only 10^{-5}pb} at low mass region. China Jin-Ping underground laboratory CJPL is the deepest underground lab in the world and provides a very promising environment for direct observation of dark matter. The China Dark Matter Experiment (CDEX) experiment is going to directly detect the WIMP flux with high sensitivity in the low mass region. Both CJPL and CDEX have achieved a remarkable progress in recent two years. The CDEX employs a point-contact germanium semi-conductor detector PCGe whose detection threshold is less than 300 eV. We report the measurement results of Muon flux, monitoring of radioactivity and Radon concentration carried out in CJPL, as well describe the structure and performance of the 1 kg PCGe detector CDEX-1 and 10kg detector array CDEX-10 including the detectors, electronics, shielding and cooling systems. Finally we discuss the physics goals of the CDEX-1, CDEX-10 and the future CDEX-1T detectors.

We collected 19.4 days of data from four 750 g TeO2 bolometers, and in three of them we were able to set the energy threshold around 3 keV using a new analysis technique. We found a background rate ranging from 25 cpd/keV/kg at 3 keV to 2 cpd/keV/kg at 25 keV, and a peak at 4.7 keV. The origin of this peak is presently unknown, but its presence is confirmed by a reanalysis of 62.7 kg.days of data from the finished CUORICINO experiment. Finally, we report the expected sensitivities of the CUORE0 (52 bolometers) and CUORE (988 bolometers) experiments to a WIMP annual modulation signal.

The extragalactic background light (EBL) contains important information about stellar and galaxy evolution. It leaves imprint on the very high energy $\gamma$-ray spectra from sources at cosmological distances due to the process of pair production. In this work we propose to {\em measure} the EBL directly by extracting the collective attenuation effects in a number of $\gamma$-ray sources at different redshifts. Using a Markov Chain Monte Carlo fitting method, the EBL intensities and the intrinsic spectral parameters of $\gamma$-ray sources are derived simultaneously. No prior shape of EBL is assumed in the fit. With this method, we can for the first time to derive the spectral shape of the EBL model-independently. Our result shows the expected features predicted by the present EBL models and thus support the understanding of the EBL origin.

Spectral energy distribution (SED) in protostellar disks is determined by the disks'internal dissipation and reprocessing of irradiation from their host stars. Around T Tauri stars, most mid-infrared (MIR) radiation (in a few to a few ten {\mu}m wavelength range) emerge from regions around a fraction to a few AU's. This region is interesting because it contains both the habitable zone and the snow line. Recent observations reveal that SED variations, in the MIR wavelength range. These variations are puzzling because they occur on time scale (a few days) which is much shorter than the dynamical (months to years) time scale at 1AU to a few AU's. They are probably caused by shadows casted by inner onto outer disk regions. Interaction between disks and their misaligned magnetized host stars can lead to warped structure and periodic SED modulations. Rapid aperiodic SED variations may also be induced by observed X-ray flares from T Tauri stars. These flares can significantly modulate the ionization fraction of the gas and the net charge carried by the grains near the surface of the inner disk. The newly charged grains may be accelerated by the stellar or disk magnetic field and adjust their distances from the midplane. Shadows casted by these grains attenuates the flux of stellar photons irradiated onto regions at several AU's from the central stars. We use this model to account for the observed rapid aperiodic SED variabilities. We suggest regular monitoring of SED variations will not only provide valuable information on the distribution of the disk aspect ratio near the habitable zone but also provide a probe on the interaction between the inner regions of the disk with the magnetosphere of their host stars.

We have systematically studied the X-ray emission properties of globular cluster millisecond pulsars in order to evaluate their spectral properties and luminosities in a uniform way. Cross-correlating the radio timing positions of the cluster pulsars with the high resolution Chandra images revealed 31 X-ray counterparts identified in nine different globular cluster systems, including those in 47 Tuc. Timing analysis has been performed for all sources corresponding to the temporal resolution available in the archival Chandra data. Making use of unpublished data on M28, M4 and NGC 6752 allowed us to obtain further constraints for the millisecond pulsar counterparts located in these clusters. Counting rate and energy flux upper limits were computed for those 36 pulsars for which no X-ray counterparts could be detected. Comparing the X-ray and radio pulse profiles of PSR J1821-2452 in M28 and the 47 Tuc pulsars PSR J0024-7204D,O,R indicated some correspondence between both wavebands. The X-ray efficiency of the globular cluster millisecond pulsars was found to be in good agreement with the efficiency Lx ~ 10^-3 Edot observed in Galactic field rotation-powered pulsars. Millisecond pulsars in the galactic plane and in globular clusters appear to show no distinct differences in their X-ray emission properties.

We report on XMM-Newton observations of the "Black Widow pulsar", PSR B1957+20. The pulsar's X-ray emission is non-thermal and best modeled with a single powerlaw spectrum of photon index 2.03(+0.51/-0.36). No coherent X-ray pulsations at the pulsar's spin-period could be detected, though a strong binary-phase dependence of the X-ray flux is observed for the first time. The data suggest that the majority of the pulsar's X-radiation is emitted from a small part of the binary orbit only. We identified this part as being near to where the radio eclipse takes place. This could mean that the X-rays from PSR B1957+20 are mostly due to intra-shock emission which is strongest when the pulsar wind interacts with the ablated material from the companion star.

We report on XMM-Newton observations of the Black Widow pulsar, PSR B1957+20. The pulsar's X-ray emission is non-thermal and best modeled with a single powerlaw spectrum of photon index 2.03^{+0.51}_{-0.36}. No coherent X-ray pulsations at the pulsar's spin-period could be detected, though a strong binary-phase dependence of the X-ray flux is observed for the first time. The data suggest that the majority of the pulsar's X-radiation is emitted from a small part of the binary orbit only. We identified this part as being near to where the radio eclipse takes place. This could mean that the X-rays from PSR B1957+20 are mostly due to intra-shock emission which is strongest when the pulsar wind interacts with the ablated material from the companion star.

We have mapped the SiO J=5-4 line at 217GHz from the HH211 molecular outflow with the Submillimeter Array (SMA). The high resolution map (1.6''x0.9'') shows that the SiO J=5-4 emission comes from the central narrow jet along the outflow axis with a width of ~0.8'' (~250 AU) FWHM. The SiO jet consists of a chain of knots separated by 3-4'' (~1000 AU) and most of the SiO knots have counterparts in shocked H_2 emission seen in a new, deep VLT near-infrared image of the outflow. A new, innermost pair of knots are discovered at just +/-2'' from the central star. The line ratio between the SiO J=5-4 data and upper limits from the SiO J=1-0 data of Chandler & Richer (2001) suggests that these knots have a temperature in excess of 300-500 K and a density of (0.5-1) x10^7 cm^{-3}. The radial velocity measured for these knots is ~30 km/s, comparable to the maximum velocity seen in the entire jet. The high temperature, high density, and velocity structure observed in this pair of SiO knots suggest that they are closely related to the primary jet launched close to the protostar.

Using the recently-commissioned Submillimeter Array (SMA), we have detected several complex organic molecules, including (CH3)2O, C2H5OH, C2H5CN, and tentatively CH2CDCN, toward the protostellar hot cores of IRAS 16293-2422 at arcsecond-resolution (~< 400 AU in radius). Vibrationally excited transitions of SO, SO2 and HCN with energy levels up to 1800 K were also observed. In addition to the other organic molecules (HC3N, CH2CO, CH3OH, CH2CHCN and HCOOCH3) previously reported by us (Kuan et al. 2004), these results clearly indicate the existence of a rich organic chemistry in low-mass `hot corinos'. From the observation of optically thin HC15N emission, we conclude I16293A is a rotating circumstellar disk lying along the north-south direction ~10 deg to the east and with an inclination ~30 deg to the sky. We suggest that the observed vibrational SO and SO2 emission may originate from shock waves near or in the circumstellar disks. Between the two cores, we find a strong anticorrelation in emission from C2H5OH and C2H5CN. The relative contribution of gas phase and grain-surface chemistries to the production of the observed complex molecules is discussed. We point out the shortcomings underlying recent claims that all the O-bearing organics are formed on grains. The presence of so many well-known interstellar molecules in solar-type hot corinos strengthens the link between molecular cloud chemistry, the starting materials of protoplanetary disks such as the protosolar nebula, and the composition of comets. Establishing the fine details of this connection is crucial in answering fundamental questions concerning the importance of galactic astrochemistry for astrobiology.

Arcsecond-resolution spectral observations toward the protobinary system IRAS 16293-2422 at 344 and 354 GHz were conducted using the Submillimeter Array. Complex organic molecules such as CH3OH and HCOOCH3 were detected. Together with the rich organic inventory revealed, it clearly indicates the existence of two, rather than one, compact hot molecular cores (smaller than or equal to 400 AU in radius) associated with each of the protobinary components identified by their dust continuum emission in the inner star-forming core.

We present results from an observing campaign of the molecular content of the coma of comet C/1999 S4 (LINEAR) carried out jointly with the millimeter-arrays of the Berkeley-Illinois-Maryland Association (BIMA) and the Owens Valley Radio Observatory (OVRO). Using the BIMA array in autocorrelation (`single-dish') mode, we detected weak HCN J=1-0 emission from comet C/1999 S4 (LINEAR) at 14 +- 4 mK km/s averaged over the 143" beam. The three days over which emission was detected, 2000 July 21.9-24.2, immediately precede the reported full breakup of the nucleus of this comet. During this same period, we find an upper limit for HCN 1-0 of 144 mJy/beam km/s (203 mK km/s) in the 9"x12" synthesized beam of combined observations of BIMA and OVRO in cross-correlation (`imaging') mode. Together with reported values of HCN 1-0 emission in the 28" IRAM 30-meter beam, our data probe the spatial distribution of the HCN emission from radii of 1300 to 19,000 km. Using literature results of HCN excitation in cometary comae, we find that the relative line fluxes in the 12"x9", 28" and 143" beams are consistent with expectations for a nuclear source of HCN and expansion of the volatile gases and evaporating icy grains following a Haser model.

We have searched three hot molecular cores for submillimeter emission from the nucleic acid building-block pyrimidine. We obtain upper limits to the total pyrimidine (beam-averaged) column densities towards Sgr B2(N), Orion KL and W51 e1/e2 of 1.7E+14 cm^{-2}, 2.4E+14 cm^{-2} and 3.4E+14 cm^{-2}, respectively. The associated upper limits to the pyrimidine fractional abundances lie in the range (0.3-3)E-10. Implications of this result for interstellar organic chemistry, and for the prospects of detecting nitrogen heterocycles in general, are briefly discussed.

We review recent studies of organic molecule formation in dense molecular clouds and in comets. We summarise the known organic inventories of molecular clouds and recent comets, particularly Hale-Bopp. The principal chemical formation pathways involving gas phase reactions, as well as formation by catalytic reactions on grain surfaces or through dust fragmentation, are identified for both dense clouds and cometary comae. The processes leading to organic molecules with known biological function, carbon chains, deuterium fractionation, HNC and S-bearing compounds are described. Observational searches for new interstellar organics are outlined and the connection between observed interstellar organics and those detected in comets Hale-Bopp and Hyakutake are discussed.