The "ATOMS," standing for {\it ALMA Three-millimeter Observations of Massive Star-forming regions}, survey has observed 146 active star forming regions with ALMA Band 3, aiming to systematically investigate the spatial distribution of various dense gas tracers in a large sample of Galactic massive clumps, to study the roles of stellar feedback in star formation, and to characterize filamentary structures inside massive clumps. In this work, the observations, data analysis, and example science of the "ATOMS" survey are presented, using a case study for the G9.62+0.19 complex. Toward this source, some transitions, commonly assumed to trace dense gas, including CS $J = 2-1$, HCO$^+$ $J = 1-0$ and HCN $J = 1-0$, are found to show extended gas emission in low density regions within the clump; less than 25\% of their emission is from dense cores. SO, CH$_3$OH, H$^{13}$CN and HC$_3$N show similar morphologies in their spatial distributions and reveal well the dense cores. Widespread narrow SiO emission is present (over $\sim$1 pc), which may be caused by slow shocks from large--scale colliding flows or H{\sc ii} regions. Stellar feedback from an expanding H{\sc ii} region has greatly reshaped the natal clump, significantly changed the spatial distribution of gas, and may also account for the sequential high-mass star formation in the G9.62+0.19 complex. The ATOMS survey data can be jointly analyzed with other survey data, e.g., "MALT90", "Orion B", "EMPIRE", "ALMA\_IMF", and "ALMAGAL", to deepen our understandings of "dense gas" star formation scaling relations and massive proto-cluster formation.

We propose a new model of Bayesian Neural Networks to not only detect the events of compact binary coalescence in the observational data of gravitational waves (GW) but also identify the time periods of the associated GW waveforms before the events. This is achieved by incorporating the Bayesian approach into the CLDNN classifier, which integrates together the Convolutional Neural Network (CNN) and the Long Short-Term Memory Recurrent Neural Network (LSTM). Our model successfully detect all seven BBH events in the LIGO Livingston O2 data, with the periods of their GW waveforms correctly labeled. The ability of a Bayesian approach for uncertainty estimation enables a newly defined `awareness' state for recognizing the possible presence of signals of unknown types, which is otherwise rejected in a non-Bayesian model. Such data chunks labeled with the awareness state can then be further investigated rather than overlooked. Performance tests show that our model recognizes 90% of the events when the optimal signal-to-noise ratio $\rho_\text{opt} >7$ (100% when $\rho_\text{opt} >8.5$) and successfully labels more than 95% of the waveform periods when $\rho_\text{opt} >8$. The latency between the arrival of peak signal and generating an alert with the associated waveform period labeled is only about 20 seconds for an unoptimized code on a moderate GPU-equipped personal computer. This makes our model possible for nearly real-time detection and for forecasting the coalescence events when assisted with deeper training on a larger dataset using the state-of-art HPCs.

We report studies of the relationships between the total bolometric luminosity ($L_{\rm bol}$ or $L_{\rm TIR}$) and the molecular line luminosities of $J=1-0$ transitions of H$^{13}$CN, H$^{13}$CO$^+$, HCN, and HCO$^+$ with data obtained from ACA observations in the "ATOMS" survey of 146 active Galactic star forming regions. The correlations between $L_{\rm bol}$ and molecular line luminosities $L'_{\rm mol}$ of the four transitions all appear to be approximately linear. Line emission of isotopologues shows as large scatters in $L_{\rm bol}$-$L'_{\rm mol}$ relations as their main line emission. The log($L_{\rm bol}$/$L'_{\rm mol}$) for different molecular line tracers have similar distributions. The $L_{\rm bol}$-to-$L'_{\rm mol}$ ratios do not change with galactocentric distances ($R_{\rm GC}$) and clump masses ($M_{\rm clump}$). The molecular line luminosity ratios (HCN-to-HCO$^+$, H$^{13}$CN-to-H$^{13}$CO$^+$, HCN-to-H$^{13}$CN and HCO$^+$-to-H$^{13}$CO$^+$) all appear constant against $L_{\rm bol}$, dust temperature ($T_{\rm d}$), $M_{\rm clump}$ and $R_{\rm GC}$. Our studies suggest that both the main lines and isotopologue lines are good tracers of the total masses of dense gas in Galactic molecular clumps. The large optical depths of main lines do not affect the interpretation of the slopes in star formation relations. We find that the mean star formation efficiency (SFE) of massive Galactic clumps in the "ATOMS" survey is reasonably consistent with other measures of the SFE for dense gas, even those using very different tracers or examining very different spatial scales.

We studied the filament structures and dense cores in OMC-2,3 region in Orion A North molecular cloud using the high-resolution N2H+ (1-0) spectral cube observed with the Atacama Large Millimeter/Submillimeter Array (ALMA). The filament network over a total length of 2 pc is found to contain 170 intersections and 128 candidate dense cores. The dense cores are all displaced from the infrared point sources (possible young stars), and the major fraction of cores (103) are located around the intersections. Towards the intersections, there is also an increasing trend for the total column density Ntot as well as the the power-law index of the column-density Probability Distribution Function (N-PDF), suggesting that the intersections would in general have more significant gas assembly than the other part of the filament paths. The virial analysis shows that the dense cores mostly have virial mass ratio of alpha_vir=M_vir/M_gas<1.0, suggesting them to be bounded by the self gravity. In the mean time, only about 23 percent of the cores have critical mass ratio of alpha_crit=M_crit/M_gas<1.0, suggesting them to be unstable against core collapse. Combining these results, it shows that the major fraction of the cold starless and possible prestellar cores in OMC-2,3 are being assembled around the intersections, and currently in a gravitationally bound state. But more extensive core collapse and star formation may still require continuous core-mass growth or other perturbatio

Radio-loud quasars (RLQs) are more X-ray luminous than predicted by the X-ray-optical/UV relation (i.e. $L_\mathrm{x}\propto L_\mathrm{uv}^\gamma$) for radio-quiet quasars (RQQs). The excess X-ray emission depends on the radio-loudness parameter ($R$) and radio spectral slope ($\alpha_\mathrm{r}$). We construct a uniform sample of 729 optically selected RLQs with high fractions of X-ray detections and $\alpha_\mathrm{r}$ measurements.We find that steep-spectrum radio quasars (SSRQs; $\alpha_\mathrm{r}\le-0.5$) follow a quantitatively similar $L_\mathrm{x}\propto L_\mathrm{uv}^\gamma$ relation as that for RQQs, suggesting a common coronal origin for the X-ray emission of both SSRQs and RQQs. However, the corresponding intercept of SSRQs is larger than that for RQQs and increases with $R$, suggesting a connection between the radio jets and the configuration of the accretion flow. Flat-spectrum radio quasars (FSRQs; $\alpha_\mathrm{r}>-0.5$) are generally more X-ray luminous than SSRQs at given $L_\mathrm{uv}$ and $R$, likely involving more physical processes. The emergent picture is different from that commonly assumed where the excess X-ray emission of RLQs is attributed to the jets. We thus perform model selection to comparecritically these different interpretations, which prefers the coronal scenario with a corona-jet connection. A distinct jet component is likely important for only a small portion of FSRQs.The corona-jet, disk-corona, and disk-jet connections of RLQs are likely driven by independent physical processes. Furthermore, the corona-jet connection implies that small-scale processesin the vicinity of SMBHs, probably associated with the magnetic flux/topology instead of black-hole spin, are controlling the radio-loudness of quasars.

We report on observations of redshifted CO(1-0) line emission and observed-frame $\rm\sim$ 30GHz radio continuum emission from five ultra-luminous, mid-IR selected hot, Dust-Obscured Galaxies (Hot DOGs) at $z\rm\gtrsim$ 3 using the Karl G. Jansky Very Large Array. We detect CO(1-0) line emission in all five Hot DOGs, with one of them at high signal to noise. We analyse FIR-radio spectral energy distributions, including dust, free-free and synchrotron emission for the galaxies. We find that most of the 115 GHz rest-frame continuum is mostly due to synchrotron or free-free emission, with only a potentially small contribution from thermal emission. We see a deficit in the rest-frame 115 GHz continuum emission compared to dusty star-forming galaxies (DSFGs) and sub-millimetre galaxies (SMGs) at high redshift, suggesting that Hot DOGs do not have similar cold gas reserves compared with star-forming galaxies. One target, W2305-0039, is detected in the FIRST 1.4 GHz survey, and is likely to possess compact radio jets. We compare to the FIR-radio correlation, and find that at least half of the Hot DOGs in our sample are radio-quiet with respect to normal galaxies. These findings suggest that Hot DOGs have comparably less cold molecular gas than star-forming galaxies at lower, $z\rm\sim$ 2 redshifts, and are dominated by powerful, yet radio-quiet AGN.

We present 11 detections of FRB 121102 in ~3 hours of observations during its 'active' period on the 10th of September 2019. The detections were made using the newly deployed MeerTRAP system and single pulse detection pipeline at the MeerKAT radio telescope in South Africa. Fortuitously, the Nancay radio telescope observations on this day overlapped with the last hour of MeerKAT observations and resulted in 4 simultaneous detections. The observations with MeerKAT's wide band receiver, which extends down to relatively low frequencies (900-1670 MHz usable L-band range), have allowed us to get a detailed look at the complex frequency structure, intensity variations and frequency-dependent sub-pulse drifting. The drift rates we measure for the full-band and sub-banded data are consistent with those published between 600-6500 MHz with a slope of -0.147 +/- 0.014 ms^-1. Two of the detected bursts exhibit fainter 'precursors' separated from the brighter main pulse by ~28 ms and ~34 ms. A follow-up multi-telescope campaign on the 6th and 8th October 2019 to better understand these frequency drifts and structures over a wide and continuous band was undertaken. No detections resulted, indicating that the source was 'inactive' over a broad frequency range during this time.

Hot Dust-Obscured Galaxies (Hot DOGs) are among the most luminous galaxies in the Universe. Powered by highly obscured, possibly Compton-thick, active galactic nuclei (AGNs), Hot DOGs are characterized by SEDs that are very red in the mid-IR yet dominated by the host galaxy stellar emission in the UV and optical. An earlier study identified a sub-sample of Hot DOGs with significantly enhanced UV emission. One target, W0204-0506, was studied in detail and, based on Chandra observations, it was concluded that the enhanced emission was most likely due to either extreme unobscured star-formation (${\rm SFR}>1000~M_{\odot}~\rm yr^{-1}$) or to light from the highly obscured AGN scattered by gas or dust into our line of sight. Here, we present a follow-up study of W0204-0506 as well as two more Hot DOGs with excess UV emission. For the two new objects we obtained Chandra/ACIS-S observations, and for all three targets we obtained HST/WFC3 F555W and F160W imaging. We conclude that the excess UV emission is primarily dominated by light from the central highly obscured, hyper-luminous AGN that has been scattered into our line of sight. We cannot rule out, however, that star-formation may significantly contribute to the UV excess of W0204-0506.

In this work, we present results that simultaneously constrain the dust size distribution, dust-to-gas ratio, fraction of dust re-deposition, and total mass production rates for comet 67P/Churyumov-Gerasimenko. We use a 3D Direct Simulation Monte Carlo (DSMC) gas dynamics code to simulate the inner gas coma of the comet for the duration of the Rosetta mission. The gas model is constrained by ROSINA/COPS data. Further, we simulate for different epochs the inner dust coma using a 3D dust dynamics code including gas drag and the nucleus' gravity. Using advanced dust scattering properties these results are used to produce synthetic images that can be compared to the OSIRIS data set. These simulations allow us to constrain the properties of the dust coma and the total gas and dust production rates. We determined a total volatile mass loss of $(6.1 \pm 1.5) \cdot 10^9$~kg during the 2015 apparition. Further, we found that power-laws with $q=3.7^{+0.57}_{-0.078}$ are consistent with the data. This results in a total of $5.1^{+6.0}_{-4.9}\cdot10^9$~kg of dust being ejected from the nucleus surface, of which $4.4^{+4.9}_{-4.2}\cdot10^9$~kg escape to space and $6.8^{+11}_{-6.8}\cdot10^8$~kg (or an equivalent of $14^{+22}_{-14}$~cm over the smooth regions) is re-deposited on the surface. This leads to a dust-to-gas ratio of $0.73^{+1.3}_{-0.70}$ for the escaping material and $0.84^{+1.6}_{-0.81}$ for the ejected material. We have further found that the smallest dust size must be strictly smaller than $\sim30\mu$m and nominally even smaller than $\sim12\mu$m.

The mass, origin and evolutionary stage of the binary system LB-1 has been the subject of intense debate, following the claim that it hosts an $\sim$70$M_{\odot}$ black hole, in stark contrast with the expectations for stellar remnants in the Milky Way. We conducted a high-resolution, phase-resolved spectroscopic study of the near-infrared Paschen lines in this system, using the 3.5-m telescope at Calar Alto Observatory. We find that Pa$\beta$ and Pa$\gamma$ (after proper subtraction of the stellar absorption component) are well fitted with a standard double-peaked model, typical of disk emission. We measured the velocity shifts of the red and blue peaks at 28 orbital phases: the line center has an orbital motion in perfect antiphase with the stellar motion, and the radial velocity amplitude ranges from 8 to 13 km/s for different choices of lines and profile modelling. We interpret this curve as proof that the disk is tracing the orbital motion of the primary, ruling out the circumbinary disk and the hierarchical triple scenarios. The phase-averaged peak-to-peak half-separation (proxy for the projected rotational velocity of the outer disk) is $\sim$70 km s$^{-1}$, larger than the stellar orbital velocity and also inconsistent with a circumbinary disk. From those results, we infer a primary mass 4--8 times higher than the secondary mass. Moreover, we show that the ratio of the blue and red peaks (V/R intensity ratio) has a sinusoidal behaviour in phase with the secondary star, which can be interpreted as the effect of external irradiation by the secondary star on the outer disk. Finally, we briefly discuss our findings in the context of alternative scenarios recently proposed for LB-1. Definitive tests between alternative solutions will require further astrometric data from $Gaia$.

To investigate the use of saliency-map analysis to aid in searches for transient signals, such as fast radio bursts and individual pulses from radio pulsars. We aim to demonstrate that saliency maps provide the means to understand predictions from machine learning algorithms and can be implemented in piplines used to search for transient events. We have implemented a new deep learning methodology to predict whether or not any segment of the data contains a transient event. The algorithm has been trained using real and simulated data sets. We demonstrate that the algorithm is able to identify such events. The output results are visually analysed via the use of saliency maps. We find that saliency maps can produce an enhanced image of any transient feature without the need for de-dispersion or removal of radio frequency interference. Such maps can be used to understand which features in the image were used in making the machine learning decision and to visualise the transient event. Even though the algorithm reported here was developed to demonstrate saliency-map analysis, we have detected, in archival data, a single burst event with dispersion measure of $41$\,cm$^{-3}$pc that is not associated with any currently known pulsar.

The European Space Agency's (ESA) Rosetta mission has returned a vast data set of measurements of the inner gas coma of comet 67P/Churyumov-Gerasimenko. These measurements have been used by different groups to determine the distribution of the gas sources at the nucleus surface. The solutions that have been found differ from each other substantially and illustrate the degeneracy of this issue. It is the aim of this work to explore the limitations that current gas models have in linking the coma measurements to the surface. In particular, we discuss the sensitivity of Rosetta's ROSINA/COPS, VIRTIS, and MIRO instruments to differentiate between vastly different spatial distributions of the gas emission from the surface. We have applied a state of the art 3D DSMC gas dynamics code to simulate the inner gas coma of different models that vary in the fraction of the surface that contains ice and in different sizes of active patches. These different distributions result in jet interactions that differ in their dynamical behaviour. We have found that ROSINA/COPS measurements by themselves cannot detect the differences in our models. While ROSINA/COPS measurements are important to constrain the regional inhomogeneities of the gas emission, they can by themselves not determine the surface emission distribution of the gas sources to a spatial accuracy of better than a few hundred metres (400 m). Any solutions fitting the ROSINA/COPS measurements is hence fundamentally degenerate, be it through a forward or inverse model. Only other instruments with complementary measurements can potentially lift this degeneracy as we show here for VIRTIS and MIRO. Finally, as a by-product, we have explored the effect of our activity distributions on lateral flow at the surface that may be responsible for some of the observed aeolian features.

We report results of a project to map HCN and HCO+ J = 1-0 emission toward a sample of molecular clouds in the inner Galaxy, all containing dense clumps that are actively engaged in star formation. We compare these two molecular line tracers with millimeter continuum emission and extinction, as inferred from 13CO, as tracers of dense gas in molecular clouds. The fraction of the line luminosity from each tracer that comes from the dense gas, as measured by AV > 8 mag, varies substantially from cloud to cloud. In all cases, a substantial fraction (in most cases, the majority) of the total luminosity arises in gas below the AV > 8 mag threshold and outside the region of strong mm continuum emission. Measurements of the luminosity of HCN toward other galaxies will likely be dominated by such gas at lower surface density. Substantial, even dominant, contributions to the total line luminosity can arise in gas with densities typical of the cloud as a whole (densities about 100 per cubic cm). Defining the dense clump from the HCN or HCO+ emission itself, similarly to previous studies, leads to a wide range of clump properties, with some being considerably larger and less dense than in previous studies. HCN and HCO+ have similar ability to trace dense gas for the clouds in this sample. For the two clouds with low virial parameters, the 13CO is definitely a worse tracer of the dense gas, but for the other four, it is equally good (or bad) at tracing dense gas.

We report the detection of diffuse Lyman-alpha halos (LAHs) around star-forming galaxies at $z \sim 5.7$ by stacking 310 spectroscopically confirmed Lyman-alpha (LAEs). The majority of the LAEs are identified from our spectroscopic survey of galaxies at $z>5.5$. They are all located in well-studied fields with deep narrowband and broadband imaging data. We combine the LAE sample and its subsamples in the narrowband NB816 (i.e., the Lyman-alpha band) and $z$ band (i.e., the continuum band). By comparing the stacked objects with their corresponding point spread functions, we clearly detect extended LAHs around these LAEs. We perform sophisticated simulations and analyses on statistical and systematic errors, and confirm that the detected halos are not caused by errors. The scale lengths of the LAHs, when described by a double-component model, range from 1.2 to 5.3 kpc with a typical value of $\sim$2 kpc. The LAH sizes from our sample are in agreement with those of individual LAEs at the similar redshift measured by VLT/MUSE, but are relatively smaller than those of photometrically selected LAEs in previous studies. We also find that LAEs with higher Lyman-alpha luminosities, higher UV continuum luminosities or smaller Lyman-alpha equivalent widths tend to have larger LAH sizes. Our results are consistent with a scenario that LAHs originate from the scattered light of the central galaxies by H I gas in the circumgalactic medium.

The Beijing-Arizona Sky Survey (BASS) is a wide and deep imaging survey to cover a 5400 deg$^2$ area in the Northern Galactic Cap with the 2.3m Bok telescope using two filters ($g$ and $r$ bands). The Mosaic $z$-band Legacy Survey (MzLS) covers the same area in $z$ band with the 4m Mayall telescope. These two surveys will be used for spectroscopic targeting of the Dark Energy Spectroscopic Instrument (DESI). The BASS survey observations were completed in 2019 March. This paper describes the third data release (DR3) of BASS, which contains the photometric data from all BASS and MzLS observations between 2015 January and 2019 March. The median astrometric precision relative to {\it Gaia} positions is about 17 mas and the median photometric offset relative to the PanSTARRS1 photometry is within 5 mmag. The median $5\sigma$ AB magnitude depths for point sources are 24.2, 23.6, and 23.0 mag for $g$, $r$, and $z$ bands, respectively. The photometric depth within the survey area is highly homogeneous, with the difference between the 20\% and 80\% depth less than 0.3 mag. The DR3 data, including raw data, calibrated single-epoch images, single-epoch photometric catalogs, stacked images, and co-added photometric catalogs, are publicly accessible at \url{http://batc.bao.ac.cn/BASS/doku.php?id=datarelease:home}.

We report the discovery of an extreme X-ray flux rise (by a factor of > 20) of the weak-line quasar SDSS J153913.47+395423.4 (hereafter SDSS J1539+3954) at z = 1.935. SDSS J1539+3954 is the most-luminous object among radio-quiet type 1 AGNs where such dramatic X-ray variability has been observed. Before the X-ray flux rise, SDSS J1539+3954 appeared X-ray weak compared with the expectation from its UV flux; after the rise, the ratio of its X-ray flux and UV flux is consistent with the majority of the AGN population. We also present a contemporaneous HET spectrum of SDSS J1539+3954, which demonstrates that its UV continuum level remains generally unchanged despite the dramatic increase in the X-ray flux, and its C iv emission line remains weak. The dramatic change only observed in the X-ray flux is consistent with a shielding model, where a thick inner accretion disk can block our line of sight to the central X-ray source. This thick inner accretion disk can also block the nuclear ionizing photons from reaching the high-ionization broad emission-line region, so that weak high-ionization emission lines are observed. Under this scenario, the extreme X-ray variability event may be caused by slight variations in the thickness of the disk. This event might also be explained by gravitational light-bending effects in a reflection model.

We present the results of processing an additional 44% of the High Time Resolution Universe South Low Latitude (HTRU-S LowLat) pulsar survey, the most sensitive blind pulsar survey of the southern Galactic plane to date. Our partially-coherent segmented acceleration search pipeline is designed to enable the discovery of pulsars in short, highly-accelerated orbits, while our 72-min integration lengths will allow us to discover pulsars at the lower end of the pulsar luminosity distribution. We report the discovery of 40 pulsars, including three millisecond pulsar-white dwarf binary systems (PSRs J1537-5312, J1547-5709 and J1618-4624), a black-widow binary system (PSR J1745-23) and a candidate black-widow binary system (PSR J1727-2951), a glitching pulsar (PSR J1706-4434), an eclipsing binary pulsar with a 1.5-yr orbital period (PSR J1653-45), and a pair of long spin-period binary pulsars which display either nulling or intermittent behaviour (PSRs J1812-15 and J1831-04). We show that the total population of 100 pulsars discovered in the HTRU-S LowLat survey to date represents both an older and lower-luminosity population, and indicates that we have yet to reach the bottom of the luminosity distribution function. We present evaluations of the performance of our search technique and of the overall yield of the survey, considering the 94% of the survey which we have processed to date. We show that our pulsar yield falls below earlier predictions by approximately 25% (especially in the case of millisecond pulsars), and discuss explanations for this discrepancy as well as future adaptations in RFI mitigation and searching techniques which may address these shortfalls.

A galaxy's morphological features encode details about its gas content, star formation history, and feedback processes which regulate its growth and evolution. We use deep convolutional neural networks (CNNs) to capture all of a galaxy's morphological information in order to estimate its neutral atomic hydrogen (HI) content directly from SDSS $gri$ image cutouts. We are able to predict a galaxy's HI mass fraction, $\mathcal M \equiv M_{\rm HI}/M_\star$, to within 0.25~dex accuracy using CNNs. The HI-morphology connection learned by the CNN appears to be constant in low- to intermediate-density galaxy environments, but it breaks down in the highest-density environments, i.e., for normalized overdensity parameter $\log(1+\delta_5) \gtrsim 0.5$ for the ALFALFA $\alpha.40$ sample, $\log(1+\delta_5) \gtrsim 0.1$ for the xGASS representative sample. This transition can be physically interpreted as the onset of ram pressure stripping, tidal effects, and other gas depletion processes in clustered environments. We also use a visualization algorithm, Gradient-weighted Class Activation Maps (Grad-CAM), to determine which morphological features are associated with low or high gas content. These results demonstrate that CNNs are powerful tools for understanding the connections between optical morphology and other properties, as well as for probing other latent variables, in a quantitative and interpretable manner.

We present VLT/XSHOOTER rest-frame UV-optical spectra of 10 Hot Dust-Obscured Galaxies (Hot DOGs) at $z\sim2$ to investigate AGN diagnostics and to assess the presence and effect of ionized gas outflows. Most Hot DOGs in this sample are narrow-line dominated AGN (type 1.8 or higher), and have higher Balmer decrements than typical type 2 quasars. Almost all (8/9) sources show evidence for ionized gas outflows in the form of broad and blueshifted [O III] profiles, and some sources have such profiles in H$\alpha$ (5/7) or [O II] (3/6). Combined with the literature, these results support additional sources of obscuration beyond the simple torus invoked by AGN unification models. Outflow rates derived from the broad [O III] line ($\rm \gtrsim10^{3}\,M_{\odot}\,yr^{-1}$) are greater than the black hole accretion and star formation rates, with feedback efficiencies ($\sim0.1-1\%$) consistent with negative feedback to the host galaxy's star formation in merger-driven quasar activity scenarios. We find the broad emission lines in luminous, obscured quasars are often better explained by outflows within the narrow line region, and caution that black hole mass estimates for such sources in the literature may have substantial uncertainty. Regardless, we find lower bounds on the Eddington ratio for Hot DOGs near unity.

All stellar mass black holes have hitherto been identified by X-rays emitted by gas that is accreting onto the black hole from a companion star. These systems are all binaries with black holes below 30 M$_{\odot}$$^{1-4}$. Theory predicts, however, that X-ray emitting systems form a minority of the total population of star-black hole binaries$^{5,6}$. When the black hole is not accreting gas, it can be found through radial velocity measurements of the motion of the companion star. Here we report radial velocity measurements of a Galactic star, LB-1, which is a B-type star, taken over two years. We find that the motion of the B-star and an accompanying H$\alpha$ emission line require the presence of a dark companion with a mass of $68^{+11}_{-13}$ M$_{\odot}$, which can only be a black hole. The long orbital period of 78.9 days shows that this is a wide binary system. The gravitational wave experiments have detected similarly massive black holes$^{7,8}$, but forming such massive ones in a high-metallicity environment would be extremely challenging to current stellar evolution theories$^{9-11}$.

We present VLT/SINFONI near-infrared (NIR) integral field spectroscopy of six $z \sim 0.2$ Lyman break galaxy "analogs" (LBAs), from which we detect HI, HeI, and [FeII] recombination lines, and multiple H$_2$ ro-vibrational lines in emission. Pa$\alpha$ kinematics reveal high velocity dispersions and low rotational velocities relative to random motions ($\langle v/\sigma \rangle = 1.2 \pm 0.8$). Matched-aperture comparisons of H$\beta$, H$\alpha$, and Pa$\alpha$ reveal that the nebular color excesses are lower relative to the continuum color excesses than is the case for typical local star-forming systems. We compare observed HeI/HI recombination line ratios to photoionization models to gauge the effective temperatures (T$_{\rm eff}$) of massive ionizing stars, finding the properties of at least one LBA are consistent with extra heating from an active galactic nucleus (AGN) and/or an overabundance of massive stars. We use H$_2$ 1-0 S($\cdot$) ro-vibrational spectra to determine rotational excitation temperature $T_{\rm ex} \sim 2000$ K for warm molecular gas, which we attribute to UV heating in dense photon-dominated regions. Spatially resolved NIR line ratios favor excitation by massive, young stars, rather than supernovae or AGN feedback. Our results suggest that the local analogs of Lyman break galaxies are primarily subject to strong feedback from recent star formation, with evidence for AGN and outflows in some cases.

The conversion between photons and axionlike particles (ALPs) in the Milky Way magnetic field could result in the detectable oscillation phenomena in $\gamma$-ray spectra of Galactic sources. In this work, the GeV (Fermi-LAT) and TeV (MAGIC/VERITAS/H.E.S.S.) data of three bright supernova remnants (SNRs, ie. IC443, W51C and W49B) have been adopted together to search such the oscillation effect. Different from our previous analysis of the sole Fermi-LAT data of IC443, we do not find any reliable signal for the photon-ALP oscillation in the joint broadband spectrum of each SNR. The reason for the inconsistence is that in this work we use the latest revision (P8R3) of Fermi-LAT data, updated diffuse emission templates and the new version of the source catalog (4FGL), which lead to some modification of the GeV spectrum of IC443. Then we set constraints on ALP parameters based on the combined analysis of all the three sources. Though these constraints are somewhat weaker than limits from the CAST experiment and globular clusters, they are supportive of and complementary to these other results.

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.

The Low-Frequency Array radio telescope discovered the $707$ Hz binary millisecond pulsar (MSP) J0952$-$0607 in a targeted radio pulsation search of an unidentified $\textit{Fermi}$ gamma-ray source. This source shows a weak energy flux of $F_\gamma = 2.6 \times 10^{-12}\,\text{erg}\,\text{cm}^{-2}\,\text{s}^{-1}$ in the energy range between $100\,\text{MeV}$ and $100\,\text{GeV}$. Here we report the detection of pulsed gamma-ray emission from PSR$\,$J0952$-$0607 in a very sensitive gamma-ray pulsation search. The pulsar's rotational, binary, and astrometric properties are measured over seven years of $\textit{Fermi}$-Large Area Telescope data. For this we take into account the uncertainty on the shape of the gamma-ray pulse profile. We present an updated radio-timing solution now spanning more than two years and show results from optical modeling of the black-widow-type companion based on new multi-band photometric data taken with HiPERCAM on the Gran Telescopio Canarias on La Palma and ULTRACAM on the New Technology Telescope at ESO La Silla. PSR$\,$J0952$-$0607 is now the fastest-spinning pulsar for which the intrinsic spin-down rate has been reliably constrained ($\dot{P}_\text{int} \lesssim 4.6 \times 10^{-21}\,\text{s}\,\text{s}^{-1}$). The inferred surface magnetic field strength of $B_\text{surf} \lesssim 8.2 \times 10^{7}\,\text{G}$ is among the ten lowest of all known pulsars. This discovery is another example of an extremely fast spinning black-widow pulsar hiding within an unidentified $\textit{Fermi} gamma-ray source. In the future such systems might help to pin down the maximum spin frequency and the minimum surface magnetic field strength of MSPs.

Most dynamically confirmed stellar-mass black holes and the candidates were originally selected from X-ray outbursts. In the present work, we search for black hole candidates in the LAMOST survey by using the spectra along with photometry from the ASAS-SN survey, where the orbital period of the binary may be revealed by the periodic light curve, such as the ellipsoidal modulation type. Our sample consists of 9 binaries, where each source contains a giant star with large radial velocity variation ($\Delta V_{\rm R} > 70~{\rm km~s^{-1}}$) and periods known from light curves. We focus on the 9 sources with long periods ($T_{\rm ph} > 5$ days) and evaluate the mass $M_2$ of the optically invisible companion. Since the observed $\Delta V_{\rm R}$ from only a few repeating spectroscopic observations is a lower limit of the real amplitude, the real mass $M_2$ can be significantly higher than the current evaluation. It is likely an efficient method to place constraints on $M_2$ by combining $\Delta V_{\rm R}$ from LAMOST and $T_{\rm ph}$ from ASAS-SN, particularly by the ongoing LAMOST Medium Resolution Survey.

We present our piggyback search for fast radio bursts using the Nanshan 26m Radio Telescope and the Kunming 40m Radio Telescope. The observations are performed in the L-band from 1380 MHz to 1700 MHz at Nanshan and S-band from 2170 MHz to 2310 MHz at Kunming. We built the \textsc{Roach2}-based FFT spectrometer and developed the real-time transient search software. We introduce a new radio interference mitigation technique named \emph{zero-DM matched filter} and give the formula of the signal-to-noise ratio loss in the transient search. Though we have no positive detection of bursts in about 1600 and 2400 hours data at Nanshan and Kunming respectively, an intriguing peryton was detected at Nanshan, from which hundreds of bursts were recorded. Perytons are terrestrial radio signals that mimic celestial fast radio bursts. They were first reported at Parkes and identified as microwave oven interferences later. The bursts detected at Nanshan show similar frequency swept emission and have double-peaked profiles. They appeared in different sky regions in about tens of minutes observations and the dispersion measure index is not exactly 2, which indicates the terrestrial origin. The peryton differs drastically from the known perytons detected at Parkes, because it appeared in a precise period of $p=1.71287\pm 0.00004$ s. Its origin remains unknown.

Due to ionosphere absorption and the interference by natural and artificial radio emissions, ground observation of the sky at the decameter or longer is very difficult. This unexplored part of electromagnetic spectrum has the potential of great discoveries, notably in the study of cosmic dark ages and dawn, but also in heliophysics and space weather, planets, cosmic ray and neutrinos, pulsar and interstellar medium, extragalactic radio sources, and even SETI. At a forum organized by the International Space Science Institute-Beijing (ISSI-BJ), we discussed the prospect of opening up this window for astronomical observations by using a constellation of small or micro-satellites. We discussed the past experiments and the current ones such as the low frequency payload on Chang'e-4 mission lander, relay satellite and the Longjiang satellite, and also the future DSL mission, which is a linear array on lunar orbit which can make synthesized map of the whole sky as well as measure the global spectrum. We also discuss the synergy with other experiments, including ground global experiments such as EDGES, SARAS, SCI-HI and High-z, PRIZM/Albatros, ground imaging facillities such as LOFAR and MWA, and space experiments such as SUNRISE, DARE/DAPPER and PRATUSH. We also discussed some technical aspects of the DSL concept.

The Gamma-Ray Integrated Detectors (GRID) is a space mission concept dedicated to monitoring the transient gamma-ray sky in the energy range from 10 keV to 2 MeV using scintillation detectors onboard CubeSats in low Earth orbits. The primary targets of GRID are the gamma-ray bursts (GRBs) in the local universe. The scientific goal of GRID is, in synergy with ground-based gravitational wave (GW) detectors such as LIGO and VIRGO, to accumulate a sample of GRBs associated with the merger of two compact stars and study jets and related physics of those objects. It also involves observing and studying other gamma-ray transients such as long GRBs, soft gamma-ray repeaters, terrestrial gamma-ray flashes, and solar flares. With multiple CubeSats in various orbits, GRID is unaffected by the Earth occultation and serves as a full-time and all-sky monitor. Assuming a horizon of 200 Mpc for ground-based GW detectors, we expect to see a few associated GW-GRB events per year. With about 10 CubeSats in operation, GRID is capable of localizing a faint GRB like 170817A with a 90% error radius of about 10 degrees, through triangulation and flux modulation. GRID is proposed and developed by students, with considerable contribution from undergraduate students, and will remain operated as a student project in the future. The current GRID collaboration involves more than 20 institutes and keeps growing. On August 29th, the first GRID detector onboard a CubeSat was launched into a Sun-synchronous orbit and is currently under test.

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.

We discuss the asymmetry of cosmic background evolution in time with respect to the quantum bounce in the Loop Quantum Cosmology (LQC), employing the value of scalar field at the bounce $\phi_{\rm B}$. We use the Chaotic and the $R^2$ potentials to demonstrate that a possible deflation before the bounce may counteract the inflation that is needed for resolving the cosmological conundrums, so a certain level of time asymmetry is required for the models in LQC. This $\phi_{\rm B}$ is model dependent and closely related to the amounts of deflation and inflation, so we may use observations to confine $\phi_{\rm B}$ and thus the model parameters. With further studies this formalism should be useful in providing an observational testbed for the LQC models.

With the observational advance in recent years, primordial gravitational waves (GWs), known as the tensor-mode cosmic perturbations, in the Loop Quantum Cosmology (LQC) are becoming testable and thus require better framework through which to bridge between the observations and the theories. In this work we present a new formalism that employs the transfer functions to bring the GWs from any epoch, even before the quantum bounce, to a later time, including the present. The evolutionary epochs considered here include the possible deflation, quantum bounce, and inflation. This formalism enables us to predict more accurately the GW power spectrum today. With the ADM formalism for the LQC background dynamics, our approach is equivalent to the commonly used Bogoliubov transformations for evolving the primordial GWs, but more transparent for discussions and easier to calculate due to its nature of being linear algebra dealing with linear perturbations. We utilize this advantage to have resolved the IR suppression problem. We also propose the field-free approximation for the effective mass in the quantum bounce epoch to largely improve the accuracy in the predicted GW power spectrum. Our transfer-function formalism is general in dealing with any linear problems, and thus expected to be equally useful under other context with linearity.

The Five-hundred-meter Aperture Spherical radio Telescope(FAST) is expected to complete its commissioning in 2019. FAST will soon begin the Commensal Radio Astronomy FasT Survey(CRAFTS), a novel and unprecedented commensal drift scan survey of the entire sky visible from FAST. The goal of CRAFTS is to cover more than 20000 $deg^{2}$ and reach redshift up to about 0.35. We provide empirical measurements of the beam size and sensitivity of FAST across the 1.05 to 1.45 GHz frequency range of the FAST L-band Array of 19-beams(FLAN). Using a simulated HI-galaxy catalogue based on the HI Mass Function(HIMF), we estimate the number of galaxies that CRAFTS may detect. At redshifts below 0.35, over $6\, \times \, 10^{5}$ HI galaxies may be detected. Below the redshift of 0.07, the CRAFTS HIMF will be complete above a mass threshold of $10^{9.5}\,M_{\odot}$. FAST will be able to investigate the environmental and redshift dependence of the HIMF to an unprecedented depth, shedding light onto the missing baryon and missing satellite problems.

Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions.

We propose a method to search for stellar-mass black hole (BH) candidates with giant companions from spectroscopic observations. Based on the stellar spectra of LAMOST Data Release 6, we obtain a sample of seven giants in binaries with large radial velocity variation $\Delta V_R > 80~{\rm km~s^{-1}}$. With the effective temperature, surface gravity, and metallicity provided by LAMOST, and the parallax given by {\it Gaia}, we can estimate the mass and radius of the giant, and therefore evaluate the possible mass of the optically invisible star in the binary. We show that the sources in our sample are potential BH candidates, and are worthy of dynamical measurement by further spectroscopic observations. Our method may be particularly valid for the selection of BH candidates in binaries with unknown orbital periods.

We train a deep residual convolutional neural network (CNN) to predict the gas-phase metallicity ($Z$) of galaxies derived from spectroscopic information ($Z \equiv 12 + \log(\rm O/H)$) using only three-band $gri$ images from the Sloan Digital Sky Survey. When trained and tested on $128 \times 128$-pixel images, the root mean squared error (RMSE) of $Z_{\rm pred} - Z_{\rm true}$ is only 0.085 dex, vastly outperforming a trained random forest algorithm on the same data set (RMSE $=0.130$ dex). The amount of scatter in $Z_{\rm pred} - Z_{\rm true}$ decreases with increasing image resolution in an intuitive manner. We are able to use CNN-predicted $Z_{\rm pred}$ and independently measured stellar masses to recover a mass-metallicity relation with $0.10$ dex scatter. Because our predicted MZR shows no more scatter than the empirical MZR, the difference between $Z_{\rm pred}$ and $Z_{\rm true}$ can not be due to purely random error. This suggests that the CNN has learned a representation of the gas-phase metallicity, from the optical imaging, beyond what is accessible with oxygen spectral lines.

We report initial results from a large Gemini program to observe z>~5.7 quasars with GNIRS near-IR spectroscopy. Our sample includes 50 quasars with simultaneous ~0.85-2.5 micron spectra covering the rest-frame ultraviolet and major broad emission lines from Ly-alpha to MgII. We present spectral measurements for these quasars and compare to their lower-redshift counterparts at z=1.5-2.3. We find that when quasar luminosity is matched, there are no significant differences between the rest-UV spectra of z>~5.7 quasars and the low-z comparison sample. High-z quasars have similar continuum and emission line properties and occupy the same region in the black hole mass and luminosity space as the comparison sample, accreting at an average Eddington ratio of ~0.3. There is no evidence for super-Eddington accretion or hypermassive (>10^10 Msun) black holes within our sample. We find a mild excess of quasars with weak CIV lines relative to the control sample. Our results, corroborating earlier studies but with better statistics, demonstrate that these high-z quasars are already mature systems of accreting supermassive black holes operating with the same physical mechanisms as those at lower redshifts.

The DArk Matter Particle Explorer (DAMPE) developed in China was designed to search for evidence of dark matter particles by observing primary cosmic rays and gamma rays in the energy range from 5 GeV to 10 TeV. Since its launch in December 2015, a large quantity of data has been recorded. With the data set acquired during more than a year of operation in space, a precise time-dependent calibration for the energy measured by the BGO ECAL has been developed. In this report, the instrumentation and development of the BGO Electromagnetic Calorimeter (BGO ECAL) are briefly described. The calibration on orbit, including that of the pedestal, attenuation length, minimum ionizing particle peak, and dynode ratio, is discussed, and additional details about the calibration methods and performance in space are presented.

The secondary-to-primary ratios are unique tools to investigate cosmic ray propagation mechanisms. In this work, we use the latest data of deuteron-to-helium~4 ratio and helium~3-to-helium~4 ratio measured by PAMELA combined with other Z$\leq$2 primary fluxes measured by PAMELA and Voyager-1, to constrain the cosmic ray acceleration and propagation models. The analysis is performed by interfacing statistical tools with the GALPROP propagation package. To better fit both the modulated and unmodulated low energy cosmic ray data, we find that a time-, charge- and rigidity-dependent solar modulation model is better than the force-field approximation. Among all the studied cosmic ray propagation models, the diffusion-reacceleration-convection model is strongly supported by the derived Bayesian evidence. The robustness of the estimated diffusion slope $\delta$ is cross-checked by another low-mass secondary-to-primary ratio, i.e. the antiproton-to-proton ratio. It is shown that the diffusion-reacceleration-convection model can reconcile well with the high energy antiproton-to-proton ratio. This indicates that the estimated value of $\delta$ is reliable. The well constraint $\delta$ from the `best' model is found to be close to 1/3, inferring a Kolmogorov-type interstellar magnetic turbulence.

We have observed the environments of a population of 33 heavily dust obscured, ultra-luminous, high-redshift galaxies, selected using WISE and NVSS at $z>$1.3 with the Infra-Red Array Camera on the $Spitzer$ Space Telescope over $\rm5.12\,'\times5.12\,'$ fields. Colour selections are used to quantify any potential overdensities of companion galaxies in these fields. We find no significant excess of galaxies with the standard colour selection for IRAC colours of $\rm[3.6]-[4.5]>-0.1$ consistent with galaxies at $z>$1.3 across the whole fields with respect to wide-area $Spitzer$ comparison fields, but there is a $\rm>2\sigma$ statistical excess within $\rm0.25\,'$ of the central radio-WISE galaxy. Using a colour selection of $\rm[3.6]-[4.5]>0.4$, 0.5 magnitudes redder than the standard method of selecting galaxies at $z>$1.3, we find a significant overdensity, in which $\rm76\%$ ($\rm33\%$) of the 33 fields have a surface density greater than the $\rm3\sigma$ ($\rm5\sigma$) level. There is a statistical excess of these redder galaxies within $\rm0.5\,'$, rising to a central peak $\rm\sim2$--4 times the average density. This implies that these galaxies are statistically linked to the radio-WISE selected galaxy, indicating similar structures to those traced by red galaxies around radio-loud AGN.

Galaxy mergers and gas accretion from the cosmic web drove the growth of galaxies and their central black holes at early epochs. We report spectroscopic imaging of a multiple merger event in the most luminous known galaxy, WISE J224607.56-052634.9 (W2246-0526), a dust-obscured quasar at redshift 4.6, 1.3 Gyr after the Big Bang. Far-infrared dust continuum observations show three galaxy companions around W2246-0526 with disturbed morphologies, connected by streams of dust likely produced by the dynamical interaction. The detection of tidal dusty bridges shows that W2246-0526 is accreting its neighbors, suggesting merger activity may be a dominant mechanism through which the most luminous galaxies simultaneously obscure and feed their central supermassive black holes.

One of the main purposes of the DArk Matter Particle Explorer (DAMPE) is to measure the cosmic ray nuclei up to several tens of TeV or beyond, whose origin and propagation remains a hot topic in astrophysics. The Plastic Scintillator Detector (PSD) on top of DAMPE is designed to measure the charges of cosmic ray nuclei from H to Fe and serves as a veto detector for discriminating gamma-rays from charged particles. We propose in this paper a charge reconstruction procedure to optimize the PSD performance in charge measurement. Essentials of our approach, including track finding, alignment of PSD, light attenuation correction, quenching and equalization correction are described detailedly in this paper after a brief description of the structure and operational principle of the PSD. Our results show that the PSD works very well and almost all the elements in cosmic rays from H to Fe are clearly identified in the charge spectrum.

We have investigated the jet-linked \mbox{X-ray} emission from highly radio-loud quasars (HRLQs; $\log R>2.5$) at high redshift. We studied the X-ray properties of 15 HRLQs at $z>4$, using new {\it Chandra} observations for six objects and archival {\it XMM-Newton} and {\it Swift} observations for the other nine. We focused on testing the apparent enhancement of jet-linked \mbox{X-ray} emission from HRLQs at $z>4$. Utilizing an enlarged (24 objects) optically flux-limited sample with complete X-ray coverage, we confirmed that HRLQs at $z>4$ have enhanced X-ray emission relative to that of HRLQs at $z\approx$ 1--2 with matched UV/optical and radio luminosity, at a \mbox{4.0--4.6}~$\sigma$ level; the X-ray enhancements are confirmed considering both two-point spectral indices and inspection of broad-band spectral energy distributions. The typical factor of enhancement is revised to $1.9^{+0.5}_{-0.4}$, which is smaller than but consistent with previous results. A fractional IC/CMB model can still explain our results at high redshift, which puts tighter constraints on the fraction of IC/CMB X-rays at lower redshifts, assuming the physical properties of quasar jets do not have a strong redshift dependence. A dominant IC/CMB model is inconsistent with our data.

Galaxy clusters trace the largest structures of the Universe and provide ideal laboratories for studying galaxy evolution and cosmology. Clusters with extended X-ray emission have been discovered at redshifts up to z ~ 2.5. Meanwhile, there has been growing interest in hunting for protoclusters, the progenitors of clusters, at higher redshifts. It is, however, very challenging to find the largest protoclusters at early times when they start to assemble. Here we report a giant protocluster of galaxies at redshift z = 5.7, when the Universe was only one billion years old. This protocluster occupies a volume of about 35x35x35 cubic co-moving megaparsecs. It is embedded in an even larger overdense region with at least 41 spectroscopically confirmed, luminous Lyman-alpha emitting galaxies (Lyman-alpha Emitters, or LAEs), including several previously reported LAEs. Its LAE density is 6.6 times the average density at z ~ 5.7. It is the only one of its kind in an LAE survey in four square degrees on the sky. Such a large structure is also rarely seen in current cosmological simulations. This protocluster will collapse into a galaxy cluster with a mass of (3.6+/-0.9) x 10^{15} solar masses, comparable to those of the most massive clusters or protoclusters known to date.

We use optical and near-infrared spectroscopy to observe rest-UV emission lines and estimate the black hole mass of WISEA J224607.56-052634.9 (W2246-0526) at z = 4.601, the most luminous hot dust-obscured galaxy yet discovered by WISE. From the broad component of the MgII-2799A emission line, we measure a black hole mass of log (M_BH/M_sun) = 9.6 +- 0.4. The broad CIV-1549A line is asymmetric and significantly blue-shifted. The derived M_BH from the blueshift-corrected broad CIV line width agrees with the MgII result. From direct measurement using a well-sampled SED, the bolometric luminosity is 3.6 * 10^14 L_sun. The corresponding Eddington ratio for W2246-0526 is lambda_Edd = L_AGN / L_Edd = 2.8. This high Eddington ratio may reach the level where the luminosity is saturating due to photon trapping in the accretion flow, and be insensitive to the mass accretion rate. In this case, the M_BH growth rate in W2246-0526 would exceed the apparent accretion rate derived from the observed luminosity.

The Plastic Scintillator Detector (PSD) of the DArk Matter Particle Explorer (DAMPE) is designed to measure cosmic ray charge (Z) and to act as a veto detector for gamma-ray identification. In order to fully exploit the charge identification potential of the PSD and to enhance its capability to identify the gamma ray events, we develop a PSD detector alignment method. The path length of a given track in the volume of a PSD bar is derived taking into account the shift and rotation alignment corrections. By examining energy spectra of corner-passing events and fully contained events, position shifts and rotations of all PSD bars are obtained, and found to be on average about 1mm and 0.0015 radian respectively. To validate the alignment method, we introduce the artificial shifts and rotations of PSD bars in the detector simulation. These shift and rotation parameters can be recovered successfully by the alignment procedure. As a result of the PSD alignment procedure, the charge resolution of the PSD is improved from $4\%$ to $8\%$ depending on the nuclei.

The Beijing-Arizona Sky Survey (BASS) is a new wide-field legacy imaging survey in the northern Galactic cap using the 2.3m Bok telescope. The survey will cover about 5400 deg$^2$ in the $g$ and $r$ bands, and the expected 5$\sigma$ depths (corrected for the Galactic extinction) in the two bands are 24.0 and 23.4 mag, respectively. BASS started observations in January 2015, and has completed about 41% of the whole area as of July 2016. The first data release contains both calibrated images and photometric catalogs obtained in 2015 and 2016. The depths of single-epoch images in the two bands are 23.4 and 22.9 mag, and the full depths of three epochs are about 24.1 and 23.5 mag, respectively.

We present XMM-Newton imaging spectroscopy of ten weak emission-line quasars (WLQs) at $0.928\leq z \leq 3.767$, six of which are radio quiet and four which are radio intermediate. The new X-ray data enabled us to measure the power-law photon index, at rest-frame energies $>2$ keV, in each source with relatively high accuracy. These measurements allowed us to confirm previous reports that WLQs have steeper X-ray spectra, suggesting higher accretion rates with respect to 'typical' quasars. A comparison between the photon indices of our radio-quiet WLQs and those of a control sample of 85 sources shows that the first are significantly higher, at the >~$3\sigma$ level. Collectively, the four radio-intermediate WLQs have lower photon indices with respect to the six radio-quiet WLQs, as may be expected if the spectra of the first group are contaminated by X-ray emission from a jet. Therefore, in the absence of significant jet emission along our line of sight, these results are in agreement with the idea that WLQs constitute the extreme high end of the accretion rate distribution in quasars. We detect soft excess emission in our lowest-redshift radio-quiet WLQ, in agreement with previous findings suggesting that the prominence of this feature is associated with a high accretion rate. We have not detected signatures of Compton reflection, Fe K$\alpha$ lines, or strong variability between two X-ray epochs in any of our WLQs, which can be attributed to their relatively high luminosity.

We present X-ray and multiwavelength analyses of 32 weak emission-line quasars (WLQs) selected in a consistent and unbiased manner. New $Chandra$ 3.1-4.8 ks observations were obtained for 14 of these WLQs with C IV rest-frame equivalent widths (REWs) of 5-15 \r{A}, and these serve as an X-ray observational "bridge" between previously studied WLQs with C IV REW $\lesssim$ 5 \r{A} and more-typical quasars with C IV REW $\approx$ 15-100 \r{A}. We have identified and quantified a strong dependence of the fraction of X-ray weak quasars upon C IV REW; this fraction declines by a factor of $\approx 13$ (from $\approx 44$% to $\approx 3$%) for C IV REW ranging from 4-50 \r{A}, and the rate of decline appears particularly strong in the 10-20 \r{A} range. The dependence broadly supports the proposed "shielding" model for WLQs, in which a geometrically and optically thick inner accretion disk, expected for a quasar accreting at a high Eddington ratio, both prevents ionizing EUV/X-ray photons from reaching the high-ionization broad emission-line region and also sometimes blocks the line-of-sight to the central X-ray emitting region. This model is also supported by the hard average spectral shape of X-ray weak WLQs (with a power-law effective photon index of $\Gamma_{\rm eff}=1.19^{+0.56}_{-0.45}$). Additionally, we have examined UV continuum/emission-line properties that might trace X-ray weakness among WLQs, confirming that red UV continuum color is the most-effective tracer.

We carried out multi-color optical monitoring of a sample of ten blazars from 2005 to 2011. The sample contains 3 LBLs, 2 IBLs, 4 HBLs, and 1 FSRQ. Our monitoring focused on the long-term variability and the sample included nine BL Lac objects and one flat-spectrum radio quasar. A total number of 14799 data points were collected. This is one of the largest optical database for a sample of ten blazars. All objects showed significant variability except OT 546. Because of the low sampling on each single night, only BL Lacertae was observed to have intra-day variability on 2006 November 6. Most BL Lac objects showed a bluer-when-brighter chromatism, while the flat-spectrum radio quasar, 3C 454.3, displayed a redder-when-brighter trend. The BWB color behaviors of most BL Lacs can be at least partly interpreted by the fact of increasing variation amplitude with increasing frequency observed in these objects. The average spectral index of LBLs is around 1.5, as expected from the model dominated by Synchrotron Self-Compton (SSC) loss. The optical emission of HBL is probably contaminated by the thermal emission from the host galaxies. Correlation analysis did not reveal any time delay between variations at different wavelengths.

We perform quasi-simultaneous optical multi-band monitoring of BL Lac object S5 0716+714 on seven nights from 2013 to 2016. Intra-day variability (IDV) is found on all seven nights. The source was faintest on JD 2456322 with 14.15 mags and brightest on JD 2457437 with 12.51 mags in the $R$ band. The maximum intra-day variation we observed is 0.15 mags in the $B$ band on JD 2456322. Both bluer-when-brighter and achromatic spectral behaviours were observed on the intra-day timescale. On the longer-term scale, the object exhibited a mild bluer-when-brighter behaviour between the $B$ and $R$ bands. We estimate the inter-band lags using two independent methods. The variation in the $B$ band was observed to lag that in the $I$ band by about 15 minutes on JD 2457315. We compare this lag with one reported previously and discussed the origin of these lags.

CREAM (Cosmic Ray Energetics And Mass) is a multi-flight balloon mission designed to collect direct data on the elemental composition and individual energy spectra of cosmic rays. Two instrument suites have been built to be flown alternately on a yearly base. The tungsten/Sci-Fi imaging calorimeter for the second flight, scheduled for December 2005, was calibrated with electron and proton beams at CERN. A calibration procedure based on the study of the longitudinal shower profile is described and preliminary results of the beam test are presented.

CREAM (Cosmic Ray Energetics And Mass) is a multi-flight balloon mission designed to collect direct data on the elemental composition and individual energy spectra of cosmic rays. Two instrument suites have been built to be flown alternately on a yearly base. The tungsten/Sci-Fi imaging calorimeter for the second flight, scheduled for December 2005, was calibrated with electron and proton beams at CERN. A calibration procedure based on the study of the longitudinal shower profile is described and preliminary results of the beam test are presented.

We consider the Palatini formulation of $f(R,T)$ gravity theory, in which a nonminimal coupling between the Ricci scalar and the trace of the energy-momentum tensor is introduced, by considering the metric and the affine connection as independent field variables. The field equations and the equations of motion for massive test particles are derived, and we show that the independent connection can be expressed as the Levi-Civita connection of an auxiliary, energy-momentum trace dependent metric, related to the physical metric by a conformal transformation. Similarly to the metric case, the field equations impose the non-conservation of the energy-momentum tensor. We obtain the explicit form of the equations of motion for massive test particles in the case of a perfect fluid, and the expression of the extra-force, which is identical to the one obtained in the metric case. The thermodynamic interpretation of the theory is also briefly discussed. We investigate in detail the cosmological implications of the theory, and we obtain the generalized Friedmann equations of the $f(R,T)$ gravity in the Palatini formulation. Cosmological models with Lagrangians of the type $f=R-\alpha ^2/R+g(T)$ and $f=R+\alpha ^2R^2+g(T)$ are investigated. These models lead to evolution equations whose solutions describe accelerating Universes at late times.

We performed a molecular line survey of 82 high-mass protostellar objects in search for inflow signatures associated with high-mass star formation. Using the H$^{13}$CO$^+$ (1$-$0) line as an optically thin tracer, we detected a statistically significant excess of blue asymmetric line profiles in the HCO$^+$ (1$-$0) transition, but nonsignificant excesses in the HCO$^+$ (3$-$2) and H$_2$CO (2$_{12}$$-$1$_{11}$) transitions. The negative blue excess for the HCN (3$-$2) transition suggests that the line profiles are affected by dynamics other than inflow motion. The HCO$^+$ (1$-$0) transition thus seems to be the suitable tracer of inflow motions in high-mass star-forming regions, as previously suggested. We found 27 inflow candidates that have at least one blue asymmetric profile and no red asymmetric profile, and derived the inflow velocities to be 0.23$-$2.00 km s$^{-1}$ for 20 of them using a simple two-layer radiative transfer model. Our sample is divided into two groups in different evolutionary stages. The blue excess of the group in relatively earlier evolutionary stages was estimated to be slightly higher than that of the other in the HCO$^+$ (1$-$0) transition.

The three Antarctic Survey Telescopes (AST3) aim to carry out time domain imaging survey at Dome A, Antarctica. The first of the three telescopes (AST3-1) was successfully deployed on January 2012. AST3-1 is a 500\,mm aperture modified Schmidt telescope with a 680\,mm diameter primary mirror. AST3-1 is equipped with a SDSS $i$ filter and a 10k $\times$ 10k frame transfer CCD camera, reduced to 5k $\times$ 10k by electronic shuttering, resulting in a 4.3 deg$^2$ field-of-view. To verify the capability of AST3-1 for a variety of science goals, extensive commissioning was carried out between March and May 2012. The commissioning included a survey covering 2000 deg$^2$ as well as the entire Large and Small Magellanic Clouds. Frequent repeated images were made of the center of the Large Magellanic Cloud, a selected exoplanet transit field, and fields including some Wolf-Rayet stars. Here we present the data reduction and photometric measurements of the point sources observed by AST3-1. We have achieved a survey depth of 19.3\,mag in 60 s exposures with 5\,mmag precision in the light curves of bright stars. The facility achieves sub-mmag photometric precision under stable survey conditions, approaching its photon noise limit. These results demonstrate that AST3-1 at Dome A is extraordinarily competitive in time-domain astronomy, including both quick searches for faint transients and the detection of tiny transit signals.

The 85-cm telescope at the Xinglong station is a well-operated prime focus system with high science outputs. The telescope has been upgraded since 2014 with new corrector, filters and camera, which are provided by Beijing Normal University (BNU). The filter set is Johnson-Cousins UBVRI system. We report the test results of the new system including the bias, dark current, linearity, gain and readout noise of the CCD camera . Then we derive accurate instrumental calibration coefficients in UBVRI bands with Landolt standard stars in the photometric nights. Finally, we give the limiting magnitudes with various exposure time and signal-to-noise ratio for observers as references.

Millisecond pulsars (MSPs) are old neutron stars that spin hundreds of times per second and appear to pulsate as their emission beams cross our line of sight. To date, radio pulsations have been detected from all rotation-powered MSPs. In an attempt to discover radio-quiet gamma-ray MSPs, we used the aggregated power from the computers of tens of thousands of volunteers participating in the Einstein@Home distributed computing project to search for pulsations from unidentified gamma-ray sources in Fermi Large Area Telescope data. This survey discovered two isolated MSPs, one of which is the only known rotation-powered MSP to remain undetected in radio observations. These gamma-ray MSPs were discovered in completely blind searches without prior constraints from other observations, raising hopes for detecting MSPs from a predicted Galactic bulge population.

Axionlike-particles (ALPs) are one promising type of dark matter candidate particle that may generate detectable effects on $\gamma$-ray spectra other than the canonical weakly interacting massive particles. In this work we search for such oscillation effects in the spectra of supernova remnants caused by the photon-ALP conversion, using the Fermi Large Area Telescope data. Three bright supernova remnants, IC443, W44, and W51C, are analyzed. The inclusion of photon-ALP oscillations yields an improved fit to the $\gamma$-ray spectrum of IC443, which gives a statistical significance of $4.2\sigma$ in favor of such spectral oscillation. However, the best-fit parameters of ALPs ($m_{a}=6.6\,{\rm neV}$, $g_{a\gamma}=13.4 \times 10^{-11}\,{\rm GeV}^{-1}$) are in tension with the upper bound ($g_{a\gamma}< 6.6 \times 10^{-11}\,{\rm GeV}^{-1}$) set by the CAST experiment. It is difficult to explain the results using the systematic uncertainties of the flux measurements. We speculate that the "irregularity" displayed in the spectrum of IC443 may be due to the superposition of the emission from different parts of the remnant.

We revisit the proposed extended Schmidt law (Shi et al. 2011) which points that the star formation efficiency in galaxies depends on the stellar mass surface density, by investigating spatially-resolved star formation rates (SFRs), gas masses and stellar masses of star formation regions in a vast range of galactic environments, from the outer disks of dwarf galaxies to spiral disks and to merging galaxies as well as individual molecular clouds in M33. We find that these regions are distributed in a tight power-law as Sigma_SFR ~(Sigma_star^0.5 Sigma_gas )^1.09, which is also valid for the integrated measurements of disk and merging galaxies at high-z. Interestingly, we show that star formation regions in the outer disks of dwarf galaxies with Sigma_SFR down to 10^(-5) Msun/yr/kpc^2, which are outliers of both Kennicutt-Schmidt and Silk-Elmegreen law, also follow the extended Schmidt law. Other outliers in the Kennicutt-Schmidt law, such as extremely-metal poor star-formation regions, also show significantly reduced deviations from the extended Schmidt law. These results suggest an important role for existing stars in helping to regulate star formation through the effect of their gravity on the mid-plane pressure in a wide range of galactic environments.

We present new Herschel observations of four massive, Sunyaev-Zel'dovich Effect (SZE)-selected clusters at $0.3 \leq z \leq 1.1$, two of which have also been observed with ALMA. We detect 19 Herschel/PACS counterparts to spectroscopically confirmed cluster members, five of which have redshifts determined via CO($4-3$) and [CI](${}^3P_1 - {}^3P_0$) lines. The mean [CI]/CO line ratio is $0.19 \pm 0.07$ in brightness temperature units, consistent with previous results for field samples. We do not detect significant stacked ALMA dust continuum or spectral line emission, implying upper limits on mean interstellar medium (H$_2$ + HI) and molecular gas masses. An apparent anticorrelation of $L_{IR}$ with clustercentric radius is driven by the tight relation between star formation rate and stellar mass. We find average specific star formation rate log(sSFR/yr$^{-1}$) = -10.36, which is below the SFR$-M_*$ correlation measured for field galaxies at similar redshifts. The fraction of infrared-bright galaxies (IRBGs; $\log (L_{IR}/L_\odot) > 10.6$) per cluster and average sSFR rise significantly with redshift. For CO detections, we find $f_{gas} \sim 0.2$, comparable to those of field galaxies, and gas depletion timescales of about 2 Gyr. We use radio observations to distinguish active galactic nuclei (AGNs) from star-forming galaxies. At least four of our 19 Herschel cluster members have $q_{IR} < 1.8$, implying an AGN fraction $f_{AGN} \gtrsim 0.2$ for our PACS-selected sample.

We report on the analysis of 13 gamma-ray pulsars discovered in the Einstein@Home blind search survey using Fermi Large Area Telescope (LAT) Pass 8 data. The 13 new gamma-ray pulsars were discovered by searching 118 unassociated LAT sources from the third LAT source catalog (3FGL), selected using the Gaussian Mixture Model (GMM) machine learning algorithm on the basis of their gamma-ray emission properties being suggestive of pulsar magnetospheric emission. The new gamma-ray pulsars have pulse profiles and spectral properties similar to those of previously-detected young gamma-ray pulsars. Follow-up radio observations have revealed faint radio pulsations from two of the newly-discovered pulsars, and enabled us to derive upper limits on the radio emission from the others, demonstrating that they are likely radio-quiet gamma-ray pulsars. We also present results from modeling the gamma-ray pulse profiles and radio profiles, if available, using different geometric emission models of pulsars. The high discovery rate of this survey, despite the increasing difficulty of blind pulsar searches in gamma rays, suggests that new systematic surveys such as presented in this article should be continued when new LAT source catalogs become available.

Hot, Dust-Obscured Galaxies, or "Hot DOGs", are a rare, dusty, hyperluminous galaxy population discovered by the WISE mission. Predominantly at redshifts 2-3, they include the most luminous known galaxies in the universe. Their high luminosities likely come from accretion onto highly obscured super massive black holes (SMBHs). We have conducted a pilot survey to measure the SMBH masses of five z~2 Hot DOGs via broad H_alpha emission lines, using Keck/MOSFIRE and Gemini/FLAMINGOS-2. We detect broad H_alpha emission in all five Hot DOGs. We find substantial corresponding SMBH masses for these Hot DOGs (~ 10^{9} M_sun), and their derived Eddington ratios are close to unity. These z~2 Hot DOGs are the most luminous AGNs at given BH masses, suggesting they are accreting at the maximum rates for their BHs. A similar property is found for known z~6 quasars. Our results are consistent with scenarios in which Hot DOGs represent a transitional, high-accretion phase between obscured and unobscured quasars. Hot DOGs may mark a special evolutionary stage before the red quasar and optical quasar phases, and they may be present at other cosmic epochs.

Studies of the most luminous quasars at high redshift directly probe the evolution of the most massive black holes in the early Universe and their connection to massive galaxy formation. However, extremely luminous quasars at high redshift are very rare objects. Only wide area surveys have a chance to constrain their population. The Sloan Digital Sky Survey (SDSS) has so far provided the most widely adopted measurements of the quasar luminosity function (QLF) at $z>3$. However, a careful re-examination of the SDSS quasar sample revealed that the SDSS quasar selection is in fact missing a significant fraction of $z\gtrsim3$ quasars at the brightest end. We have identified the purely optical color selection of SDSS, where quasars at these redshifts are strongly contaminated by late-type dwarfs, and the spectroscopic incompleteness of the SDSS footprint as the main reasons. Therefore we have designed the Extremely Luminous Quasar Survey (ELQS), based on a novel near-infrared JKW2 color cut using WISE AllWISE and 2MASS all-sky photometry, to yield high completeness for very bright ($m_{\rm{i}} < 18.0$) quasars in the redshift range of $3.0\leq z\leq5.0$. It effectively uses random forest machine-learning algorithms on SDSS and WISE photometry for quasar-star classification and photometric redshift estimation. The ELQS will spectroscopically follow-up $\sim 230$ new quasar candidates in an area of $\sim12000\,\rm{deg}^2$ in the SDSS footprint, to obtain a well-defined and complete quasars sample for an accurate measurement of the bright-end quasar luminosity function at $3.0\leq z\leq5.0$. In this paper we present the quasar selection algorithm and the quasar candidate catalog.

Hot Dust-Obscured Galaxies (Hot DOGs) are hyperluminous ($L_{\mathrm{8-1000\,\mu m}}>10^{13}\,\mathrm{L_\odot}$) infrared galaxies with extremely high (up to hundreds of K) dust temperatures. The sources powering both their extremely high luminosities and dust temperatures are thought to be deeply buried and rapidly accreting supermassive black holes (SMBHs). Hot DOGs could therefore represent a key evolutionary phase in which the SMBH growth peaks. X-ray observations can be used to study their obscuration levels and luminosities. In this work, we present the X-ray properties of the 20 most-luminous ($L_{\mathrm{bol}}\gtrsim10^{14}\, L_\odot$) known Hot DOGs at $z=2-4.6$. Five of them are covered by long-exposure ($10-70$ ks) Chandra and XMM-Newton observations, with three being X-ray detected, and we study their individual properties. One of these sources (W0116$-$0505) is a Compton-thick candidate, with column density $N_H=(1.0-1.5)\times10^{24}\,\mathrm{cm^{-2}}$ derived from X-ray spectral fitting. The remaining 15 Hot DOGs have been targeted by a Chandra snapshot (3.1 ks) survey. None of these 15 is individually detected; therefore we applied a stacking analysis to investigate their average emission. From hardness-ratio analysis, we constrained the average obscuring column density and intrinsic luminosity to be log$N_H\,\mathrm{[cm^{-2}]}>23.5$ and $L_X\gtrsim10^{44}\,\mathrm{erg\,cm^{-2}\,s^{-1}}$, which are consistent with results for individually detected sources. We also investigated the $L_X-L_{6\mu\mathrm{m}}$ and $L_X-L_{bol}$ relations, finding hints that Hot DOGs are typically X-ray weaker than expected, although larger samples of luminous obscured QSOs are needed to derive solid conclusions.

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.

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.

In this paper, we presented a study of temperature effects on BGO calorimeters using proton MIP's collected in the first year operation of DAMPE. By directly comparing MIP calibration constants used by DAMPE data production pipe line, we found an experimental relation between temperature and signal amplitudes of each BGO bar: a general deviation of -1.162%/$^{\circ}$C,and -0.47%/$^{\circ}$C to -1.60%/$^{\circ}$C statistically for each detector element. During 2016, DAMPE's temperature changed by about 7 degrees due to solar elevation angle and the corresponding energy scale bias is about 8%. By frequent MIP calibration operation, this kind of bias is eliminated to an acceptable value.

We present a spectroscopic survey of high-redshift, luminous galaxies over four square degrees on the sky, aiming to build a large and homogeneous sample of Ly$\alpha$ emitters (LAEs) at $z\approx5.7$ and 6.5, and Lyman-break galaxies (LBGs) at $5.5<z<6.8$. The fields that we choose to observe are well-studied, such as SXDS and COSMOS. They have deep optical imaging data in a series of broad and narrow bands, allowing efficient selection of galaxy candidates. Spectroscopic observations are being carried out using the multi-object spectrograph M2FS on the Magellan Clay telescope. M2FS is efficient to identify high-redshift galaxies, owing to its 256 optical fibers deployed over a circular field-of-view 30 arcmin in diameter. We have observed $\sim2.5$ square degrees. When the program is completed, we expect to identify more than 400 bright LAEs at $z\approx5.7$ and 6.5, and a substantial number of LBGs at $z\ge6$. This unique sample will be used to study a variety of galaxy properties and to search for large protoclusters. Furthermore, the statistical properties of these galaxies will be used to probe cosmic reionization. We describe the motivation, program design, target selection, and M2FS observations. We also outline our science goals, and present a sample of the brightest LAEs at $z\approx5.7$ and 6.5. This sample contains 32 LAEs with Ly$\alpha$ luminosities higher than 10$^{43}$ erg s$^{-1}$. A few of them reach $\ge3\times10^{43}$ erg s$^{-1}$, comparable to the two most luminous LAEs known at $z\ge6$, `CR7' and `COLA1'. These LAEs provide ideal targets to study extreme galaxies in the distant universe.

We monitored BL Lacertae simultaneously in the optical B, V, R and I bands for 13 nights during the period 2012-2016. The variations were well correlated in all bands and the source showed significant intraday variability (IDV). We also studied its optical flux and colour behaviour, and searched for inter-band time lags. A strong bluer-when-brighter chromatism was found on the intra-night time-scale. The spectral changes are not sensitive to the host galaxy contribution. Cross-correlation analysis revealed possible time delay of about 10 min between variations in the V and R bands. We interpreted the observed flares in terms of the model consisting of individual synchrotron pulses.

We present HST WFC3 F160W imaging and infrared spectral energy distributions for twelve extremely luminous, obscured AGN at $1.8<z<2.7$, selected via "Hot, Dust Obscured" mid-infrared colors. Their infrared luminosities span $2-15\times10^{13}$L$_{\odot}$, making them among the most luminous objects in the Universe at $z\sim2$. In all cases the infrared emission is consistent with arising at least in most part from AGN activity. The AGN fractional luminosities are higher than those in either sub-millimeter galaxies, or AGN selected via other mid-infrared criteria. Adopting the $G$, M$_{20}$ and $A$ morphological parameters, together with traditional classification boundaries, infers that three quarters of the sample as mergers. Our sample do not, however, show any correlation between the considered morphological parameters and either infrared luminosity or AGN fractional luminosity. Moreover, their asymmetries and effective radii are distributed identically to those of massive galaxies at $z\sim2$. We conclude that our sample is not preferentially associated with mergers, though a significant merger fraction is still plausible. Instead, we propose that our sample are examples of the massive galaxy population at $z\sim2$ that harbor a briefly luminous, "flickering" AGN, and in which the $G$ and M$_{20}$ values have been perturbed, due to either the AGN, and/or the earliest formation stages of a bulge in an inside-out manner. Furthermore, we find that the mass assembly of the central black holes in our sample leads the mass assembly of any bulge component. Finally, we speculate that our sample represent a small fraction of the immediate antecedents of compact star-forming galaxies at $z\sim2$.

A software system has been developed for the DArk Matter Particle Explorer (DAMPE) mission, a satellite-based experiment. The DAMPE software is mainly written in C++ and steered using Python script. This article presents an overview of the DAMPE offline software, including the major architecture design and specific implementation for simulation, calibration and reconstruction. The whole system has been successfully applied to DAMPE data analysis, based on which some results from simulation and beam test experiments are obtained and presented.

GC-1 and GC-2 are two globular clusters (GCs) in the remote halo of M81 and M82 in the M81 group discovered by Jang et al. using the {\it Hubble Space Telescope} ({\it HST}) images. These two GCs were observed as part of the Beijing--Arizona--Taiwan--Connecticut (BATC) Multicolor Sky Survey, using 14 intermediate-band filters covering a wavelength range of 4000--10000 \AA. We accurately determine these two clusters' ages and masses by comparing their spectral energy distributions (from 2267 to 20000~{\AA}, comprising photometric data in the near-ultraviolet of the {\it Galaxy Evolution Explorer}, 14 BATC intermediate-band, and Two Micron All Sky Survey near-infrared $JHK_{\rm s}$ filters) with theoretical stellar population-synthesis models, resulting in ages of $15.50\pm3.20$ for GC-1 and $15.10\pm2.70$ Gyr for GC-2. The masses of GC-1 and GC-2 obtained here are $1.77-2.04\times 10^6$ and $5.20-7.11\times 10^6 \rm~M_\odot$, respectively. In addition, the deep observations with the Advanced Camera for Surveys and Wide Field Camera 3 on the {\it HST} are used to provide the surface brightness profiles of GC-1 and GC-2. The structural and dynamical parameters are derived from fitting the profiles to three different models; in particular, the internal velocity dispersions of GC-1 and GC-2 are derived, which can be compared with ones obtained based on spectral observations in the future. For the first time, in this paper, the $r_h$ versus $M_V$ diagram shows that GC-2 is an ultra-compact dwarf in the M81 group.

We investigate extremely luminous dusty galaxies in the environments around WISE-selected hot dust obscured galaxies (Hot DOGs) and WISE/radio-selected active galactic nuclei (AGNs) at average redshifts of z = 2.7 and z = 1.7, respectively. Previous observations have detected overdensities of companion submillimetre-selected sources around 10 Hot DOGs and 30 WISE/radio AGNs, with overdensities of ~ 2 - 3 and ~ 5 - 6 , respectively. We find that the space densities in both samples to be overdense compared to normal star-forming galaxies and submillimetre galaxies (SMGs) in the SCUBA-2 Cosmology Legacy Survey (S2CLS). Both samples of companion sources have consistent mid-IR colours and mid-IR to submm ratios as SMGs. The brighter population around WISE/radio AGNs could be responsible for the higher overdensity reported. We also find the star formation rate density (SFRDs) are higher than the field, but consistent with clusters of dusty galaxies. WISE-selected AGNs appear to be good signposts for protoclusters at high redshift on arcmin scales. The results reported here provide an upper limit to the strength of angular clustering using the two-point correlation function. Monte Carlo simulations show no angular correlation, which could indicate protoclusters on scales larger than the SCUBA-2 1.5arcmin scale maps.

Milky Way-like galaxies are predicted to host a very large number of dark matter subhalos. Some massive and nearby subhalos could generate detectable gamma-rays, appearing as unidentified, spatially-extended and stable gamma-ray sources. We search for such sources in the third Fermi Large Area Telescope source List (3FGL) and report the identification of a new candidate, 3FGL J1924.8-1034. With the Fermi-LAT Pass 8 data, we find that 3FGL J1924.8-1034 is spatially-extended at a high confidence level of $5.4\sigma$, with a best-fit extension radius of $\sim0.15^{\circ}$. No significant variability has been found and its gamma-ray spectrum is well fitted by the dark matter annihilation into $b\bar{b}$ with a mass of $\sim 43$ GeV. All these facts make 3FGL J1924.8-1034 a possible dark matter subhalo candidate. However, due to the limited angular resolution, the possibility of that the spatial extension of 3FGL J1924.8-1034 is caused by the contamination from the other un-resolved point source can not be ruled out.

Primary cosmic-ray elemental spectra have been measured with the balloon-borne Cosmic Ray Energetics And Mass (CREAM) experiment since 2004. The third CREAM payload (CREAM-III) flew for 29 days during the 2007-2008 Antarctic season. Energies of incident particles above 1 TeV are measured with a calorimeter. Individual elements are clearly separated with a charge resolution of ~0.12 e (in charge units) and ~0.14 e for protons and helium nuclei, respectively, using two layers of silicon charge detectors. The measured proton and helium energy spectra at the top of the atmosphere are harder than other existing measurements at a few tens of GeV. The relative abundance of protons to helium nuclei is 9.53+-0.03 for the range of 1 TeV/n to 63 TeV/n. The ratio is considerably smaller than other measurements at a few tens of GeV/n. The spectra become softer above ~20 TeV. However, our statistical uncertainties are large at these energies and more data are needed.

In an effort to understand the cause of the apparent depletion in the number density of radio-loud AGNs at $z>3$, this work investigates the viability of the so-called Cosmic Microwave Background (CMB) quenching mechanism of intrinsically jetted, high-z AGNs, whereby Inverse Compton scattering of CMB photons off electrons within the extended lobes results in a substantial dimming of the lobe synchrotron emission at GHz frequencies, while simultaneously boosting their diffuse X-ray signal. We focus on five $z>3.5$ radio galaxies that have sufficiently deep Chandra exposure (> 50 ks) to warrant a meaningful investigation of any extended X-ray emission. For those objects with evidence for statistically significant extended X-ray lobes (4C 41.17 and 4C 03.24), we combine the Chandra measurements with literature data at lower frequencies to assemble the systems' Spectral Energy Distributions (SEDs), and utilize state-of-the-art SED modelling (Ghisellini et al. 2015) -- including emission from the disk, torus, jet, hotspots, and lobes -- to infer their physical parameters. For both radio galaxies, the magnetic energy density in the hotspots is found to exceed the energy density in CMB photons, wheres the opposite is true for the lobes. This implies that any extended synchrotron emission likely originates from the hotspots themselves, rather than the lobes. Conversely, Inverse Compton scattering of CMB photons dominates the extended X-ray emission from the lobes, which are effectively "radio-quenched". As a result, CMB quenching is effective in these systems in spite of the fact that the observed X-ray to radio luminosity ratio does not bear the signature $(1+z)^4$ dependence of the CMB energy density.

The volume of data generated by modern astronomical telescopes is extremely large and rapidly growing. However, current high-performance data processing architectures/frameworks are not well suited for astronomers because of their limitations and programming difficulties. In this paper, we therefore present OpenCluster, an open-source distributed computing framework to support rapidly developing high-performance processing pipelines of astronomical big data. We first detail the OpenCluster design principles and implementations and present the APIs facilitated by the framework. We then demonstrate a case in which OpenCluster is used to resolve complex data processing problems for developing a pipeline for the Mingantu Ultrawide Spectral Radioheliograph. Finally, we present our OpenCluster performance evaluation. Overall, OpenCluster provides not only high fault tolerance and simple programming interfaces, but also a flexible means of scaling up the number of interacting entities. OpenCluster thereby provides an easily integrated distributed computing framework for quickly developing a high-performance data processing system of astronomical telescopes and for significantly reducing software development expenses.

We report on the results of a recent blind search survey for gamma-ray pulsars in Fermi Large Area Telescope (LAT) data being carried out on the distributed volunteer computing system, Einstein@Home. The survey has searched for pulsations in 118 unidentified pulsar-like sources, requiring about 10,000 years of CPU core time. In total, this survey has resulted in the discovery of 17 new gamma-ray pulsars, of which 13 are newly reported in this work, and an accompanying paper. These pulsars are all young, isolated pulsars with characteristic ages between 12 kyr and 2 Myr, and spin-down powers between $10^{34}$ and $4\times10^{36}$ erg/s. Two of these are the slowest spinning gamma-ray pulsars yet known. One pulsar experienced a very large glitch $\Delta f/f \approx 3.5\times10^{-6}$ during the Fermi mission. In this, the first of two associated papers, we describe the search scheme used in this survey, and estimate the sensitivity of our search to pulsations in unidentified Fermi-LAT sources. One such estimate results in an upper limit of 57% for the fraction of pulsed emission from the gamma-ray source associated with the Cas A supernova remnant, constraining the pulsed gamma-ray photon flux that can be produced by the neutron star at its center. We also present the results of precise timing analyses for each of the newly detected pulsars.

We present the results of a photometric and spectroscopic analysis of the Galactic Bulge Survey X-ray source CXOGBS J174954.5-294335 (hereafter, referred to as CX19). CX19 is a long period, eclipsing intermediate polar type cataclysmic variable with broad, single-peaked Balmer and Paschen emission lines along with HeII $\lambda4686$ and Bowen blend emission features. With coverage of one full and two partial eclipses and archival photometry, we determine the ephemeris for CX19 to be HJD(eclipse) = 2455691.8581(5) + 0.358704(2)$\times$N. We also recovered the white dwarf spin period of P$_{\rm spin}$ = 503.32(3) seconds which gives a P$_{\rm spin}$/P$_{\rm orb}$ = 0.016(6), comparable to several confirmed, long period intermediate polars. CX19 also shows a clear X-ray eclipse in the 0.3-8.0 keV range observed with Chandra. Two optical outbursts were observed lasting between 6-8 hours (lower limits) reaching $\sim$1.3 mags in amplitude. The outbursts, both in duration and magnitude, the accretion disc dominated spectra and hard X-ray emission are reminiscent of the intermediate polar V1223 Sgr sharing many of the same characteristics. If we assume a main sequence companion, we estimate the donor to be an early G-type star and find a minimum distance of $d \approx$ 2.1 kpc and a 0.5-10.0 keV X-ray luminosity upper limit of 2.0 $\times$ 10$^{33}$ erg s$^{-1}$. Such an X-ray luminosity is consistent with a white dwarf accretor in a magnetic cataclysmic variable system. To date, CX19 is only the second deeply-eclipsing intermediate polar with X-ray eclipses and the first which is optically accessible.

Hot, Dust-Obscured Galaxies (Hot DOGs), selected from the WISE all sky infrared survey, host some of the most powerful Active Galactic Nuclei (AGN) known, and might represent an important stage in the evolution of galaxies. Most known Hot DOGs are at $z> 1.5$, due in part to a strong bias against identifying them at lower redshift related to the selection criteria. We present a new selection method that identifies 153 Hot DOG candidates at $z\sim 1$, where they are significantly brighter and easier to study. We validate this approach by measuring a redshift $z=1.009$, and an SED similar to higher redshift Hot DOGs for one of these objects, WISE J1036+0449 ($L_{\rm\,Bol}\simeq 8\times 10^{46}\rm\,erg\,s^{-1}$), using data from Keck/LRIS and NIRSPEC, SDSS, and CSO. We find evidence of a broadened component in MgII, which, if due to the gravitational potential of the supermassive black hole, would imply a black hole mass of $M_{\rm\,BH}\simeq 2 \times 10^8 M_{\odot}$, and an Eddington ratio of $\lambda_{\rm\,Edd}\simeq 2.7$. WISE J1036+0449 is the first Hot DOG detected by NuSTAR, and the observations show that the source is heavily obscured, with a column density of $N_{\rm\,H}\simeq(2-15)\times10^{23}\rm\,cm^{-2}$. The source has an intrinsic 2-10 keV luminosity of $\sim 6\times 10^{44}\rm\,erg\,s^{-1}$, a value significantly lower than that expected from the mid-infrared/X-ray correlation. We also find that the other Hot DOGs observed by X-ray facilities show a similar deficiency of X-ray flux. We discuss the origin of the X-ray weakness and the absorption properties of Hot DOGs. Hot DOGs at $z\lesssim1$ could be excellent laboratories to probe the characteristics of the accretion flow and of the X-ray emitting plasma at extreme values of the Eddington ratio.

We report the discovery and timing measurements of PSR J1208-6238, a young and highly magnetized gamma-ray pulsar, with a spin period of 440 ms. The pulsar was discovered in gamma-ray photon data from the Fermi Large Area Telescope (LAT) during a blind-search survey of unidentified LAT sources, running on the distributed volunteer computing system Einstein@Home. No radio pulsations were detected in dedicated follow-up searches with the Parkes radio telescope, with a flux density upper limit at 1369 MHz of 30 $\mu$Jy. By timing this pulsar's gamma-ray pulsations, we measure its braking index over five years of LAT observations to be $n = 2.598 \pm 0.001 \pm 0.1$, where the first uncertainty is statistical and the second estimates the bias due to timing noise. Assuming its braking index has been similar since birth, the pulsar has an estimated age of around 2,700 yr, making it the youngest pulsar to be found in a blind search of gamma-ray data and the youngest known radio-quiet gamma-ray pulsar. Despite its young age the pulsar is not associated with any known supernova remnant or pulsar wind nebula. The pulsar's inferred dipolar surface magnetic field strength is $3.8 \times 10^{13}$ G, almost 90% of the quantum-critical level. We investigate some potential physical causes of the braking index deviating from the simple dipole model but find that LAT data covering a longer time interval will be necessary to distinguish between these.

We present the discovery of nine quasars at $z\sim6$ identified in the Sloan Digital Sky Survey (SDSS) imaging data. This completes our survey of $z\sim6$ quasars in the SDSS footprint. Our final sample consists of 52 quasars at $5.7<z\le6.4$, including 29 quasars with $z_{\rm AB}\le20$ mag selected from 11,240 deg$^2$ of the SDSS single-epoch imaging survey (the main survey), 10 quasars with $20\le z_{\rm AB}\le20.5$ selected from 4223 deg$^2$ of the SDSS overlap regions (regions with two or more imaging scans), and 13 quasars down to $z_{\rm AB}\approx22$ mag from the 277 deg$^2$ in Stripe 82. They span a wide luminosity range of $-29.0\le M_{1450}\le-24.5$. This well-defined sample is used to derive the quasar luminosity function (QLF) at $z\sim6$. After combining our SDSS sample with two faint ($M_{1450}\ge-23$ mag) quasars from the literature, we obtain the parameters for a double power-law fit to the QLF. The bright-end slope $\beta$ of the QLF is well constrained to be $\beta=-2.8\pm0.2$. Due to the small number of low-luminosity quasars, the faint-end slope $\alpha$ and the characteristic magnitude $M_{1450}^{\ast}$ are less well constrained, with $\alpha=-1.90_{-0.44}^{+0.58}$ and $M^{\ast}=-25.2_{-3.8}^{+1.2}$ mag. The spatial density of luminous quasars, parametrized as $\rho(M_{1450}<-26,z)=\rho(z=6)\,10^{k(z-6)}$, drops rapidly from $z\sim5$ to 6, with $k=-0.72\pm0.11$. Based on our fitted QLF and assuming an IGM clumping factor of $C=3$, we find that the observed quasar population cannot provide enough photons to ionize the $z\sim6$ IGM at $\sim90$\% confidence. Quasars may still provide a significant fraction of the required photons, although much larger samples of faint quasars are needed for more stringent constraints on the quasar contribution to reionization.

The Yuan-Tseh Lee Array for Microwave Background Anisotropy (AMiBA) is a co-planar interferometer array operating at a wavelength of 3mm to measure the Sunyaev-Zeldovich effect (SZE) of galaxy clusters. In the first phase of operation -- with a compact 7-element array with 0.6m antennas (AMiBA-7) -- we observed six clusters at angular scales from 5\arcmin to 23\arcmin. Here, we describe the expansion of AMiBA to a 13-element array with 1.2m antennas (AMiBA-13), its subsequent commissioning, and our cluster SZE observing program. The most important changes compared to AMiBA-7 are (1) array re-configuration with baselines ranging from 1.4m to 4.8m covering angular scales from 2\arcmin to 11.5\arcmin, (2) thirteen new lightweight carbon-fiber-reinforced plastic (CFRP) 1.2m reflectors, and (3) additional correlators and six new receivers. From the AMiBA-13 SZE observing program, we present here maps of a subset of twelve clusters. In highlights, we combine AMiBA-7 and AMiBA-13 observations of Abell 1689 and perform a joint fitting assuming a generalized NFW pressure profile. Our cylindrically integrated Compton-y values for this cluster are consistent with the BIMA/OVRA, SZA, and Planck results. We report the first targeted SZE detection towards the optically selected galaxy cluster RCS J1447+0828, and we demonstrate the ability of AMiBA SZE data to serve as a proxy for the total cluster mass. Finally, we show that our AMiBA-SZE derived cluster masses are consistent with recent lensing mass measurements in the literature.

The Wide Field Camera 3 (WFC3) on the Hubble Space Telescope (HST) enabled the search for the first galaxies observed at z ~ 8 - 11 (500 - 700 Myr after the Big Bang). To continue quantifying the number density of the most luminous galaxies (M_AB ~ -22.0) at the earliest epoch observable with HST, we search for z ~ 10 galaxies (F125W-dropouts) in archival data from the Brightest of Reionizing Galaxies (BoRG[z8]) survey, originally designed for detection of z ~ 8 galaxies (F098M-dropouts). By focusing on the deepest 293 arcmin^2 of the data along 62 independent lines of sight, we identify six z ~ 10 candidates satisfying the color selection criteria, detected at S/N > 8 in F160W with M_AB = -22.8 to -21.1 if at z = 10. Three of the six sources, including the two brightest, are in a single WFC3 pointing (~ 4 arcmin^2), suggestive of significant clustering, which is expected from bright galaxies at z ~ 10. However, the two brightest galaxies are too extended to be likely at z ~ 10, and one additional source is unresolved and possibly a brown dwarf. The remaining three candidates have m_AB ~ 26, and given the area and completeness of our search, our best estimate is a number density of sources that is marginally higher but consistent at 2{\sigma} with searches in legacy fields. Our study highlights that z ~ 10 searches can yield a small number of candidates, making tailored follow-ups of HST pure-parallel observations viable and effective.

As part of the Chandra Galactic Bulge Survey (GBS), we present a catalogue of optical sources in the GBS footprint. This consists of two regions centered at Galactic latitude b = 1.5 degrees above and below the Galactic Centre, spanning (l x b) = (6x1) degrees. The catalogue consists of 2 or more epochs of observations for each line of sight in r', i' and H{\alpha} filters. It is complete down to r' = 20.2 and i' = 19.2 mag; the mean 5{\sigma} depth is r' = 22.5 and i' = 21.1 mag. The mean root-mean-square residuals of the astrometric solutions is 0.04 arcsec. We cross-correlate this optical catalogue with the 1640 unique X-ray sources detected in Chandra observations of the GBS area, and find candidate optical counterparts to 1480 X-ray sources. We use a false alarm probability analysis to estimate the contamination by interlopers, and expect ~ 10 per cent of optical counterparts to be chance alignments. To determine the most likely counterpart for each X-ray source, we compute the likelihood ratio for all optical sources within the 4{\sigma} X-ray error circle. This analysis yields 1480 potential counterparts (~ 90 per cent of the sample). 584 counterparts have saturated photometry (r'<17, i'<16), indicating these objects are likely foreground sources and the real counterparts. 171 candidate counterparts are detected only in the i'-band. These sources are good qLMXB and CV candidates as they are X-ray bright and likely located in the Bulge.

We present new Chandra X-ray observations of seven low-mass black holes (~1e6 Msun) accreting at low Eddington ratios between -2.0<log L/Ledd<-1.5. We compare the X-ray properties of these seven low-mass active galactic nuclei (AGN) to a total of 73 other low-mass AGN in the literature with published Chandra observations (with Eddington ratios extending from -2.0<log L/Ledd<-0.1). We do not find any statistical differences between low- and high-Eddington ratio low-mass AGN in the distributions of their X-ray to ultraviolet luminosity ratios (aox), or in their X-ray spectral shapes. Furthermore, the aox distribution of low-L/Ledd AGN displays an X-ray weak tail that is also observed within high-L/Ledd objects. Our results indicate that between -2<log L/Ledd<-0.1, there is no systematic change in the structure of the accretion flow for active galaxies hosting 1e6 Msun black holes. We examine the accuracy of current bolometric luminosity estimates for our low-L/Ledd objects with new Chandra observations, and it is plausible that their Eddington ratios could be underestimated by up to an order of magnitude. If so, then in analogy with weak emission line quasars, we suggest that accretion from a geometrically thick, radiatively inefficient `slim disk' could explain their diverse properties in aox. Alternatively, if current Eddington ratios are in fact correct (or overestimated), then the X-ray weak tail would imply that there is diversity in disk/corona couplings among individual low-mass objects. Finally, we conclude by noting that the aox distribution for low-mass black holes may have favorable consequences for the epoch of cosmic reionization being driven by AGN.

We present a catalog of low-mass dense cores observed with the SHARC-II instrument at 350 microns. Our observations have an effective angular resolution of 10", approximately 2.5 times higher than observations at the same wavelength obtained with the Herschel Space Observatory, albeit with lower sensitivity, especially to extended emission. The catalog includes 81 maps covering a total of 164 detected sources. For each detected source, we tabulate basic source properties including position, peak intensity, flux density in fixed apertures, and radius. We examine the uncertainties in the pointing model applied to all SHARC-II data and conservatively find that the model corrections are good to within ~3", approximately 1/3 of the SHARC-II beam. We examine the differences between two array scan modes and find that the instrument calibration, beam size, and beam shape are similar between the two modes. We also show that the same flux densities are measured when sources are observed in the two different modes, indicating that there are no systematic effects introduced into our catalog by utilizing two different scan patterns during the course of taking observations. We find a detection rate of 95% for protostellar cores but only 45% for starless cores, and demonstrate the existence of a SHARC-II detection bias against all but the most massive and compact starless cores. Finally, we discuss the improvements in protostellar classification enabled by these 350 micron observations.

We analyze the physical properties and dynamical origin of a curved jet of comet 67P/Churyumov-Gerasimenko that was observed repeatedly in several nucleus rotations starting on May 30 and persisting until early August, 2015. We simulated the motion of dust grains ejected from the nucleus surface under the influence of the gravity and viscous drag effect of the expanding gas flow from the rotating nucleus. The formation of the curved jet is a combination of the size of the dust particles (~0.1-1 mm) and the location of the source region near the nucleus equator. This enhances the spiral feature of the collimated dust stream after the dust is accelerated to a terminal speed on the order of m/s.

Hot Dust-Obscured Galaxies (Hot DOGs) are a population of hyper-luminous infrared galaxies identified by the WISE mission from their very red mid-IR colors, and characterized by hot dust temperatures ($T>60~\rm K$). Several studies have shown clear evidence that the IR emission in these objects is powered by a highly dust-obscured AGN that shows close to Compton-thick absorption at X-ray wavelengths. Thanks to the high AGN obscuration, the host galaxy is easily observable, and has UV/optical colors usually consistent with those of a normal galaxy. Here we discuss a sub-population of 8 Hot DOGs that show enhanced rest-frame UV/optical emission. We discuss three scenarios that might explain the excess UV emission: (i) unobscured light leaked from the AGN by reflection over the dust or by partial coverage of the accretion disk; (ii) a second unobscured AGN in the system; or (iii) a luminous young starburst. X-ray observations can help discriminate between these scenarios. We study in detail the blue excess Hot DOG WISE J020446.13-050640.8, which was serendipitously observed by Chandra/ACIS-I for 174.5 ks. The X-ray spectrum is consistent with a single, hyper-luminous, highly absorbed AGN, and is strongly inconsistent with the presence of a secondary unobscured AGN. Based on this, we argue that the excess blue emission in this object is most likely either due to reflection or a co-eval starburst. We favor the reflection scenario as the unobscured star-formation rate needed to power the UV/optical emission would be $\gtrsim 1000~\rm M_{\odot}~\rm yr^{-1}$. Deep polarimetry observations could confirm the reflection hypothesis.

We present a dynamical study of the Galactic black hole binary system Nova Muscae 1991 (GS/GRS 1124-683). We utilize 72 high resolution Magellan Echellette (MagE) spectra and 72 strictly simultaneous V-band photometric observations; the simultaneity is a unique and crucial feature of this dynamical study. The data were taken on two consecutive nights and cover the full 10.4-hour orbital cycle. The radial velocities of the secondary star are determined by cross-correlating the object spectra with the best-match template spectrum obtained using the same instrument configuration. Based on our independent analysis of five orders of the echellette spectrum, the semi-amplitude of the radial velocity of the secondary is measured to be K_2 = 406.8+/-2.7 km/s, which is consistent with previous work, while the uncertainty is reduced by a factor of 3. The corresponding mass function is f(M) = 3.02+/-0.06 M_\odot. We have also obtained an accurate measurement of the rotational broadening of the stellar absorption lines (v sin i = 85.0+/-2.6 km/s) and hence the mass ratio of the system q = 0.079+/-0.007. Finally, we have measured the spectrum of the non-stellar component of emission that veils the spectrum of the secondary. In a future paper, we will use our veiling-corrected spectrum of the secondary and accurate values of K_2 and q to model multi-color light curves and determine the systemic inclination and the mass of the black hole.

The bright soft X-ray transient Nova Muscae 1991 was intensively observed during its entire 8-month outburst using the Large Area Counter (LAC) onboard the Ginga satellite. Recently, we obtained accurate estimates of the mass of the black hole primary, the orbital inclination angle of the system, and the distance. Using these crucial input data and Ginga X-ray spectra, we have measured the spin of the black hole using the continuum-fitting method. For four X-ray spectra of extraordinary quality we have determined the dimensionless spin parameter of the black hole to be a/M = 0.63 (-0.19, +0.16) (1 sigma confidence level), a result that we confirm using eleven additional spectra of lower quality. Our spin estimate challenges two published results: It is somewhat higher than the value predicted by a proposed relationship between jet power and spin; and we find that the spin of the black hole is decidedly prograde, not retrograde as has been claimed.

The optical counterpart of the black-hole soft X-ray transient Nova Muscae 1991 has brightened by $\Delta{V}\approx0.8$ mag since its return to quiescence 23 years ago. We present the first clear evidence that the brightening of soft X-ray transients in quiescence occurs at a nearly linear rate. This discovery, and our precise determination of the disk component of emission obtained using our $simultaneous$ photometric and spectroscopic data, have allowed us to identify and accurately model archival ellipsoidal light curves of the highest quality. The simultaneity, and the strong constraint it provides on the component of disk emission, is a key element of our work. Based on our analysis of the light curves, and our earlier measurements of the mass function and mass ratio, we have obtained for Nova Muscae 1991 the first accurate estimates of its systemic inclination $i=43.2^{+2.1}_{-2.7}$ deg, and black hole mass $M=11.0^{+2.1}_{-1.4}\ M_\odot$. Based on our determination of the radius of the secondary, we estimate the distance to be $D=4.95^{+0.69}_{-0.65}$ kpc. We discuss the implications of our work for future dynamical studies of black-hole soft X-ray transients.

We have monitored the BL Lacertae object S5 0716+714 simultaneously in the B, R and I bands on three nights in November 2014. The average time resolution is quite high (73s, 34s, 58s for the filters B, R and I), which can help us trace the profile of the variation and search for the short inter-band time delay. Intra-day variability was about 0.1 mag on the first two nights and more than 0.3 mag on the third. A bluer-when-brighter color behavior was found. An clear loop path can be seen on the color-magnitude diagram of the third night, revealing possible time delays between variations at high and low energies. It is the first time that the intra-day spectral hysteresis loop has been found so obviously in the optical band. We used the interpolated cross-correlation function method to further confirm the time delay and calculated the values of lag between light curves at different wavelengths on each night. On the third night, variations in the R and B bands is approximately 1.5 minutes lagging behind the I band. Such optical time delay is probably due to the interplay of different processes of electrons in the jet of the blazar.

We present Submillimeter Array (SMA) observations of the HCN\,(3--2) and HCO$^+$\,(3--2) molecular lines toward the W3(H$_2$O) and W3(OH) star-forming complexes. Infall and outflow motions in the W3(H$_2$O) have been characterized by observing HCN and HCO$^+$ transitions. High-velocity blue/red-shifted emission, tracing the outflow, show multiple knots, which might originate in episodic and precessing outflows. `Blue-peaked' line profiles indicate that gas is infalling onto the W3(H$_2$O) dust core. The measured large mass accretion rate, 2.3$\times$10$^{-3}$~M$_{\odot}$~yr$^{-1}$, together with the small free-fall time scale, 5$\times$10$^{3}$~yr, suggest W3(H$_2$O) is in an early evolutionary stage of the process of formation of high-mass stars. For the W3(OH), a two-layer model fit to the HCN and HCO$^+$ spectral lines and Spizter/IRAC images support that the W3(OH) H{\sc ii} region is expanding and interacting with the ambient gas, with the shocked neutral gas being expanding with an expansion timescale of 6.4$\times$10$^{3}$~yr. The observations suggest different kinematical timescales and dynamical states for the W3(H$_2$O) and W3(OH).

Observed at z = 4.601 and with L_bol = 3.5 x 10^14 Lsun, W2246-0526 is the most luminous galaxy known in the Universe, and hosts a deeply-buried active galactic nucleus (AGN)/super-massive black hole (SMBH). Discovered using the Wide-field Infrared Survey Explorer (WISE), W2246-0526 is classified as a Hot Dust Obscured Galaxy (Hot DOG), based on its luminosity and dust temperature. Here we present spatially resolved ALMA [CII]157.7um observations of W2246-0526, providing unique insight into the kinematics of its interstellar medium (ISM). The measured [CII]-to-far-infrared ratio is ~2 x 10^-4, implying ISM conditions that compare only with the most obscured, compact starbursts and AGN in the local Universe today. The spatially resolved [CII] line is strikingly uniform and very broad, 500-600 km/s wide, extending throughout the entire galaxy over about 2.5 kpc, with modest shear. Such a large, homogeneous velocity dispersion indicates a highly turbulent medium. W2246-0526 is unstable in terms of the energy and momentum that are being injected into the ISM, strongly suggesting that the gas is being blown away from the system isotropically, likely reflecting a cathartic state on its road to becoming an un-obscured quasar. W2246-0526 provides an extraordinary laboratory to study and model the properties and kinematics of gas in an extreme environment under strong feedback, at a time when the Universe was 1/10 of its current age: a system pushing the limits that can be reached during galaxy formation.

We present ALMA 870 micron (345 GHz) data for 49 high redshift (0.47<z<2.85), luminous (11.7 < log L(bol) (Lsun) < 14.2) radio-powerful AGN, obtained to constrain cool dust emission from starbursts concurrent with highly obscured radiative-mode black hole (BH) accretion in massive galaxies which possess a small radio jet. The sample was selected from WISE with extremely steep (red) mid-infrared (MIR) colors and with compact radio emission from NVSS/FIRST. Twenty-six sources are detected at 870 microns, and we find that the sample has large mid- to far-infrared luminosity ratios consistent with a dominant and highly obscured quasar. The rest-frame 3 GHz radio powers are 24.7 < log P3.0 GHz (W/Hz) < 27.3, and all sources are radio-intermediate or radio-loud. BH mass estimates are 7.7 < log M(BH) (Msun) < 10.2. The rest frame 1-5 um SEDs are very similar to the "Hot DOGs" (Hot Dust Obscured Galaxies), and steeper (redder) than almost any other known extragalactic sources. ISM masses estimated for the ALMA detected sources are 9.9 < log M(ISM) (Msun) < 11.75 assuming a dust temperature of 30K. The cool dust emission is consistent with star formation rates (SFRs) reaching several thousand Msun/yr, depending on the assumed dust temperature, however we cannot rule out the alternative that the AGN powers all the emission in some cases. Our best constrained source has radiative transfer solutions with ~ equal contributions from an obscured AGN and a young (10-15 Myr) compact starburst.

We report the discovery of PSR J1906+0722, a gamma-ray pulsar detected as part of a blind survey of unidentified Fermi Large Area Telescope (LAT) sources being carried out on the volunteer distributed computing system, Einstein@Home. This newly discovered pulsar previously appeared as the most significant remaining unidentified gamma-ray source without a known association in the second Fermi-LAT source catalog (2FGL) and was among the top ten most significant unassociated sources in the recent third catalog (3FGL). PSR J1906+0722 is a young, energetic, isolated pulsar, with a spin frequency of $8.9$ Hz, a characteristic age of $49$ kyr, and spin-down power $1.0 \times 10^{36}$ erg s$^{-1}$. In 2009 August it suffered one of the largest glitches detected from a gamma-ray pulsar ($\Delta f / f \approx 4.5\times10^{-6}$). Remaining undetected in dedicated radio follow-up observations, the pulsar is likely radio-quiet. An off-pulse analysis of the gamma-ray flux from the location of PSR J1906+0722 revealed the presence of an additional nearby source, which may be emission from the interaction between a neighboring supernova remnant and a molecular cloud. We discuss possible effects which may have hindered the detection of PSR J1906+0722 in previous searches and describe the methods by which these effects were mitigated in this survey. We also demonstrate the use of advanced timing methods for estimating the positional, spin and glitch parameters of difficult-to-time pulsars such as this.

We are conducting an archival Swift program to measure multiwavelength variability in active galactic nuclei (AGN). This variability information will provide constraints on the geometry, physical conditions and processes of the structures around the central black holes that emit and reprocess the observed flux. Among our goals are: (1) to produce a catalog of type 1 AGN with time-resolved multi-wavelength data; (2) to characterize variability in the optical, UV and X-ay bands as well as changes in spectral slope; (3) to quantify the impact of variability on multi-wavelength properties; and (4) to measure correlated variability between bands. Our initial efforts have revealed a UVOT calibration issue that can cause a few percent of measured UV fluxes to be anomalously low, by up to 30%.

During the "WISE at 5: Legacy and Prospects" conference in Pasadena, CA -- which ran from February 10 - 12, 2015 -- attendees were invited to engage in an interactive session exploring the future uses of the Wide-field Infrared Survey Explorer (WISE) data. The 65 participants -- many of whom are extensive users of the data -- brainstormed the top questions still to be answered by the mission, as well as the complementary current and future datasets and additional processing of WISE/NEOWISE data that would aid in addressing these most important scientific questions. The results were mainly bifurcated between topics related to extragalactic studies (e.g. AGN, QSOs) and substellar mass objects. In summary, participants found that complementing WISE/NEOWISE data with cross-correlated multiwavelength surveys (e.g. SDSS, Pan-STARRS, LSST, Gaia, Euclid, etc.) would be highly beneficial for all future mission goals. Moreover, developing or implementing machine-learning tools to comb through and understand cross-correlated data was often mentioned for future uses. Finally, attendees agreed that additional processing of the data such as co-adding WISE and NEOWISE and extracting a multi-epoch photometric database and parallax and proper motion catalog would greatly improve the scientific results of the most important projects identified. In that respect, a project such as MaxWISE which would execute the most important additional processing and extraction as well as make the data and catalogs easily accessible via a public portal was deemed extremely important.

We present 20 WISE-selected galaxies with bolometric luminosities L_bol > 10^14 L_sun, including five with infrared luminosities L_IR = L(rest 8-1000 micron) > 10^14 L_sun. These "extremely luminous infrared galaxies," or ELIRGs, were discovered using the "W1W2-dropout" selection criteria which requires marginal or non-detections at 3.4 and 4.6 micron (W1 and W2, respectively) but strong detections at 12 and 22 micron in the WISE survey. Their spectral energy distributions are dominated by emission at rest-frame 4-10 micron, suggesting that hot dust with T_d ~ 450K is responsible for the high luminosities. These galaxies are likely powered by highly obscured AGNs, and there is no evidence suggesting these systems are beamed or lensed. We compare this WISE-selected sample with 116 optically selected quasars that reach the same L_bol level, corresponding to the most luminous unobscured quasars in the literature. We find that the rest-frame 5.8 and 7.8 micron luminosities of the WISE-selected ELIRGs can be 30-80% higher than that of the unobscured quasars. The existence of AGNs with L_bol > 10^14 L_sun at z > 3 suggests that these supermassive black holes are born with large mass, or have very rapid mass assembly. For black hole seed masses ~ 10^3 M_sun, either sustained super-Eddington accretion is needed, or the radiative efficiency must be <15%, implying a black hole with slow spin, possibly due to chaotic accretion.

We present time-resolved optical spectroscopy of the counterpart to the high-inclination black hole low-mass X-ray binary Swift J1357.2-0933 in quiescence. Absorption features from the mass donor star were not detected. Instead the spectra display prominent broad double-peaked Halpha emission and weaker HeI emission lines. From the Halpha peak-to-peak separation we constrain the radial velocity semi-amplitude of the donor star to > 789 km/s. Further analysis through radial velocity and equivalent width measurements indicates that the Halpha line is free of variability due to S-wave components or disc eclipses. From our data and previous observations during outburst, we conclude that long-term radial velocity changes ascribed to a precessing disc were of low amplitude or not present. This implies that the centroid position of the line should closely represent the systemic radial velocity. Using the derived systemic velocity of -150 km/s and the best available limits on the source distance, we infer that the black hole is moving towards the Plane in its current Galactic orbit unless the proper motion is substantial. Finally, the depth of the central absorption in the double peaked profiles adds support for Swift J1357.2-0933 as a high-inclination system. On the other hand, we argue that the low hydrogen column density inferred from X-ray fitting suggests that the system is not seen edge-on.

We present Gemini spectroscopy for 21 candidate optical counterparts to X-ray sources discovered in the Galactic Bulge Survey (GBS). For the majority of the 21 sources, the optical spectroscopy establishes that they are indeed the likely counterparts. One of the criteria we used for the identification was the presence of an Ha emission line. The spectra of several sources revealed an Ha emission line only after careful subtraction of the F or G stellar spectral absorption lines. In a sub-class of three of these sources the residual Halpha emission line is broad (> 400 km/s) which suggests that it is formed in an accretion disk, whereas in other cases the line width is such that we currently cannot determine whether the line emission is formed in an active star/binary or in an accretion disk. GBS source CX377 shows this hidden accretion behaviour most dramatically. The previously-identified broad Ha emission of this source is not present in its Gemini spectra taken about 1 year later. However, broad emission is revealed after subtracting an F6 template star spectrum. The Gemini spectra of three sources (CX446, CX1004, and CXB2) as well as the presence of possible eclipses in light curves of these sources suggest that these sources are accreting binaries viewed under a high inclination.

Over the past 15 years, examples of exotic radio-quiet quasars with intrinsically weak or absent broad emission line regions (BELRs) have emerged from large-scale spectroscopic sky surveys. Here, we present spectroscopy of seven such weak emission line quasars (WLQs) at moderate redshifts (z=1.4-1.7) using the X-shooter spectrograph, which provides simultaneous optical and near-infrared spectroscopy covering the rest-frame ultraviolet through optical. These new observations effectively double the number of WLQs with spectroscopy in the optical rest-frame, and they allow us to compare the strengths of (weak) high-ionization emission lines (e.g., CIV) to low-ionization lines (e.g., MgII, Hb, Ha) in individual objects. We detect broad Hb and Ha emission in all objects, and these lines are generally toward the weaker end of the distribution expected for typical quasars (e.g., Hb has rest-frame equivalent widths ranging from 15-40 Ang.). However, these low-ionization lines are not exceptionally weak, as is the case for high-ionization lines in WLQs. The X-shooter spectra also display relatively strong optical FeII emission, Hb FWHM <4000 km/s, and significant CIV blueshifts (1000-5500 km/s) relative to the systemic redshift; two spectra also show elevated ultraviolet FeII emission, and an outflowing component to their (weak) MgII emission lines. These properties suggest that WLQs are exotic versions of "wind-dominated" quasars. Their BELRs either have unusual high-ionization components, or their BELRs are in an atypical photoionization state because of an unusually soft continuum.

We present JCMT SCUBA-2 850microns submillimetre (submm) observations of 30 mid-infrared (mid-IR) luminous AGN, detected jointly by the WISE all-sky IR survey and the NVSS/FIRST radio survey. These rare sources are selected by their extremely red mid-infrared spectral energy distributions (SEDs) and compact radio counterparts. Further investigations show that they are highly obscured, have abundant warm AGN-heated dust and are thought to be experiencing intense AGN feedback. These galaxies appear to be consistent with an AGN-dominated galaxy, and could be a transient phase of merging galaxies. When comparing the number of submm galaxies (SMGs) detected serendipitously in the surrounding 1.5-arcmin to those in blank-field submm surveys, there is a very significant overdensity, of order 5, but no sign of radial clustering centred at our primary objects. The WISE/radio-selected AGN thus reside in 10-Mpc-scale overdense environments, that could be forming in pre-viralised clusters of galaxies. WISE/radio-selected AGNs appear to be the strongest signposts of high-density regions of active, luminous and dusty galaxies. SCUBA-2 850microns observations indicate that their submm fluxes are low compared to many popular AGN SED templates, hence the WISE/radio-selected AGNs have either less cold and/or more warm dust emission than normally assumed for typical AGN. Most of the targets are not detected, only four targets are detected at SCUBA-2 850microns, and have total IR luminosities >= 10^13 L_solar, if their redshifts are consistent with the subset of the 10 SCUBA-2 undetected targets with known redshifts, z ~ 0.44 - 2.86.

We present an X-ray and multiwavelength study of 33 weak emission-line quasars (WLQs) and 18 quasars that are analogs of the extreme WLQ, PHL 1811, at z ~ 0.5-2.9. New Chandra 1.5-9.5 ks exploratory observations were obtained for 32 objects while the others have archival X-ray observations. Significant fractions of these luminous type 1 quasars are distinctly X-ray weak compared to typical quasars, including 16 (48%) of the WLQs and 17 (94%) of the PHL 1811 analogs with average X-ray weakness factors of 17 and 39, respectively. We measure a relatively hard ($\Gamma=1.16_{-0.32}^{+0.37}$) effective power-law photon index for a stack of the X-ray weak subsample, suggesting X-ray absorption, and spectral analysis of one PHL 1811 analog, J1521+5202, also indicates significant intrinsic X-ray absorption. We compare composite SDSS spectra for the X-ray weak and X-ray normal populations and find several optical-UV tracers of X-ray weakness; e.g., Fe II rest-frame equivalent width and relative color. We describe how orientation effects under our previously proposed "shielding-gas" scenario can likely unify the X-ray weak and X-ray normal populations. We suggest that the shielding gas may naturally be understood as a geometrically thick inner accretion disk that shields the broad line region from the ionizing continuum. If WLQs and PHL 1811 analogs have very high Eddington ratios, the inner disk could be significantly puffed up (e.g., a slim disk). Shielding of the broad emission-line region by a geometrically thick disk may have a significant role in setting the broad distributions of C IV rest-frame equivalent width and blueshift for quasars more generally.

We present a spectroscopic analysis of the black hole binary Nova Muscae 1991 in quiescence using data obtained in 2009 with MagE on the Magellan Clay telescope and in 2010 with IMACS on the Magellan Baade telescope at the Las Campanas Observatory. Emission from the disc is observed in H alpha, H beta and Ca II (8662 A). A prominent hotspot is observed in the Doppler maps of all three emission lines. The existence of this spot establishes ongoing mass transfer from the donor star in 2009-2010 and, given its absence in the 1993-1995 observations, demonstrates the presence of a variable hotspot in the system. We find the radial distance to the hotspot from the black hole to be consistent with the circularization radius. Our tomograms are suggestive of stream-disc overflow in the system. We also detect possible Ca II (8662 A) absorption from the donor star.

The WISE mission has unveiled a rare population of high-redshift ($z=1-4.6$), dusty, hyper-luminous galaxies, with infrared luminosities $L_{\rm IR} > 10^{13}~L_{\odot}$, and sometimes exceeding $10^{14}~L_{\odot}$. Previous work has shown that their dust temperatures and overall far-IR SEDs are significantly hotter than expected for star-formation. We present here an analysis of the rest-frame optical through mid-IR SEDs for a large sample of these so-called "Hot, Dust-Obscured Galaxies" (Hot DOGs). We find that the SEDs of Hot DOGs are generally well modeled by the combination of a luminous, yet obscured AGN that dominates the rest-frame emission at $\lambda > 1\mu\rm m$ and the bolometric luminosity output, and a less luminous host galaxy that is responsible for the bulk of the rest optical/UV emission. Even though the stellar mass of the host galaxies may be as large as $10^{11}-10^{12}~M_{\odot}$, the AGN emission, with luminosities comparable to those of the most luminous QSOs known, require that either Hot DOGs have black hole masses significantly in excess of the local relations, or that they radiate significantly above the Eddington limit. We show that, while rare, the number density of Hot DOGs is comparable to that of equally luminous but unobscured (i.e., Type 1) QSOs. This is inconsistent with the trend of a diminishing fraction of obscured objects with increasing luminosity found for less luminous QSOs, possibly indicating a reversal in this relation at high luminosity, and that Hot DOGs are not the torus-obscured counterparts of the known optically selected, largely unobscured Hyper-Luminous QSOs. Hot DOGs may represent a different type of galaxy and thus a new component of the galaxy evolution paradigm. Finally, we discuss the environments of Hot DOGs and show that these objects are in regions as dense as those of known high-redshift proto-clusters.(Abridged)

We report on the Submillimeter Array (SMA) observations of molecular lines at 270 GHz toward W3(OH) and W3(H$_2$O) complex. Although previous observations already resolved the W3(H$_2$O) into two or three sub-components, the physical and chemical properties of the two sources are not well constrained. Our SMA observations clearly resolved W3(OH) and W3(H$_2$O) continuum cores. Taking the advantage of the line fitting tool XCLASS, we identified and modeled a rich molecular spectrum in this complex, including multiple CH$_3$CN and CH$_3$OH transitions in both cores. HDO, C$_2$H$_5$CN, O$^{13}$CS, and vibrationally excited lines of HCN, CH$_3$CN, and CH$_3$OCHO were only detected in W3(H$_2$O). We calculate gas temperatures and column densities for both cores. The results show that W3(H$_{2}$O) has higher gas temperatures, and larger column densities than W3(OH) as previously observed, suggesting physical and chemical differences between the two cores. We compare the molecular abundances in W3(H$_2$O) to those in the Sgr B2(N) hot core, the Orion KL hot core and the Orion Compact Ridge, and discuss the chemical origin of specific species. An east-west velocity gradient is seen in W3(H$_2$O), and the extension is consistent with the bipolar outflow orientation traced by water masers and radio jets. A north-south velocity gradient across W3(OH) is also observed. However, with current observations we can not assure if the velocity gradients are caused by rotation, outflow or radial velocity differences of the sub-components in W3(OH).

This paper presents CCD multicolor photometry for 304 old star clusters in the nearby spiral galaxy M31. Of which photometry of 55 star clusters is first obtained. The observations were carried out as a part of the Beijing--Arizona--Taiwan--Connecticut (BATC) Multicolor Sky Survey from 1995 February to 2008 March, using 15 intermediate-band filters covering 3000--10000 \AA. Detailed comparisons show that our photometry is in agreement with previous measurements. Based on the ages and metallicities from Caldwell et al. and the photometric measurements here, we estimated the clusters' masses by comparing their multicolor photometry with stellar population synthesis models. The results show that the sample clusters have masses between $\sim 3\times10^4 M_\odot$ and $\sim 10^7 M_\odot$ with the peak of $\sim 4\times10^5 M_\odot$. The masses here are in good agreement with those in previous studies. Combined with the masses of young star clusters of M31 from Wang et al., we find that the peak of mass of old clusters is ten times that of young clusters.

We have investigated the field around the radio-quiet $\gamma$-ray pulsar, PSR J2021+4026, with a ~140 ks XMM-Newton observation and a ~56 ks archival Chandra data. Through analyzing the pulsed spectrum, we show that the X-ray pulsation is purely thermal in nature which suggests the pulsation is originated from a hot polar cap with $T\sim3\times10^{6}$ K on the surface of a rotating neutron star. On the other hand, the power-law component that dominates the pulsar emission in the hard band is originated from off-pulse phases, which possibly comes from a pulsar wind nebula. In re-analyzing the Chandra data, we have confirmed the presence of bow-shock nebula which extends from the pulsar to west by ~10 arcsec. The orientation of this nebular feature suggests that the pulsar is probably moving eastward which is consistent with the speculated proper motion by extrapolating from the nominal geometrical center of the supernova remnant (SNR) G78.2+2.1 to the current pulsar position. For G78.2+2.1, our deep XMM-Newton observation also enables a study of the central region and part of the southeastern region with superior photon statistics. The column absorption derived for the SNR is comparable with that for PSR J2021+4026, which supports their association. The remnant emission in both examined regions are in an non-equilibrium ionization state. Also, the elapsed time of both regions after shock-heating is apparently shorter than the Sedov age of G78.2+2.1. This might suggest the reverse shock has reached the center not long ago. Apart from PSR J2021+4026 and G78.2+2.1, we have also serendipitously detected an X-ray flash-like event XMM J202154.7+402855 from this XMM-Newton observation.

We present the results of the first systematic long-term multi-color optical monitoring of the BL Lacertae object 1ES 0806+52.4. The monitoring was performed in multiple passbands with a 60/90 cm Schmidt telescope from December 2005 to February 2011. The overall brightness of this object decreased from 2005 December to 2008 December, and regained after that. A sharp outburst probably occurred around the end of our monitoring program. Overlapped on the long-term trend are some short-term small-amplitude oscillations. No intra-night variability was found in the object, which is in accord with the historical observations before 2005. By investigating the color behavior, we found strong bluer-when-brighter chromatism for the long-term variability of 1ES 0806+52.4. The total amplitudes at the c, i and o bands are 1.18, 1.12, and 1.02 mags, respectively. The amplitudes tend to increase toward shorter wavelength, which may be the major cause of bluer-when-brighter. Such bluer-when-brighter is also found in other blazars like S5 0716+714, OJ 287, etc. The hard X-ray data collected from the Swift/BAT archive was correlated with our optical data. No positive result was found, the reason of which may be that the hard X-ray flux is a combination of the synchrotron and inverse Compton emission but with different timescales and cadences under the leptonic Synchrotron-Self-Compton (SSC) model.

WISE has discovered an extraordinary population of hyper-luminous dusty galaxies which are faint in the two bluer passbands ($3.4\, \mu$m and $4.6\, \mu$m) but are bright in the two redder passbands of WISE ($12\, \mu$m and $22\, \mu$m). We report on initial follow-up observations of three of these hot, dust-obscured galaxies, or Hot DOGs, using the CARMA and SMA interferometer arrays at submm/mm wavelengths. We report continuum detections at $\sim$ 1.3 mm of two sources (WISE J014946.17+235014.5 and WISE J223810.20+265319.7, hereafter W0149+2350 and W2238+2653, respectively), and upper limits to CO line emission at 3 mm in the observed frame for two sources (W0149+2350 and WISE J181417.29+341224.8, hereafter W1814+3412). The 1.3 mm continuum images have a resolution of 1-2 arcsec and are consistent with single point sources. We estimate the masses of cold dust are 2.0$\times 10^{8} M_{\odot}$ for W0149+2350 and 3.9$\times 10^{8} M_{\odot}$ for W2238+2653, comparable to cold dust masses of luminous quasars. We obtain 2$\sigma$ upper limits to the molecular gas masses traced by CO, which are 3.3$\times 10^{10} M_{\odot}$ and 2.3$\times 10^{10} M_{\odot}$ for W0149+2350 and W1814+3412, respectively. We also present high-resolution, near-IR imaging with WFC3 on the Hubble Space Telescope for W0149+2653 and with NIRC2 on Keck for W2238+2653. The near-IR images show morphological structure dominated by a single, centrally condensed source with effective radius less than 4 kpc. No signs of gravitational lensing are evident.

We present SCUBA-2 850um submillimetre (submm) observations of the fields of 10 dusty, luminous galaxies at z ~ 1.7 - 4.6, detected at 12um and/or 22um by the WISE all-sky survey, but faint or undetected at 3.4um and 4.6um; dubbed hot, dust-obscured galaxies (Hot DOGs). The six detected targets all have total infrared luminosities greater than 10^13 L_sun, with one greater than 10^14 L_sun. Their spectral energy distributions (SEDs) are very blue from mid-infrared to submm wavelengths and not well fitted by standard AGN SED templates, without adding extra dust extinction to fit the WISE 3.4um and 4.6um data. The SCUBA-2 850um observations confirm that the Hot DOGs have less cold and/or more warm dust emission than standard AGN templates, and limit an underlying extended spiral or ULIRG-type galaxy to contribute less than about 2% or 55% of the typical total Hot DOG IR luminosity, respectively. The two most distant and luminous targets have similar observed submm to mid-infrared ratios to the rest, and thus appear to have even hotter SEDs. The number of serendipitous submm galaxies (SMGs) detected in the 1.5-arcmin-radius SCUBA-2 850um maps indicates there is a significant over-density of serendipitous sources around Hot DOGs. These submm observations confirm that the WISE-selected ultra-luminous galaxies have very blue mid-infrared to submm SEDs, suggesting that they contain very powerful AGN, and are apparently located in unusual arcmin-scale overdensities of very luminous dusty galaxies.

We report on a NuSTAR and XMM-Newton program that has observed a sample of three extremely luminous, heavily obscured WISE-selected AGN at z~2 in a broad X-ray band (0.1 - 79 keV). The parent sample, selected to be faint or undetected in the WISE 3.4um (W1) and 4.6um (W2) bands but bright at 12um (W3) and 22um (W4), are extremely rare, with only ~1000 so-called W1W2-dropouts across the extragalactic sky. Optical spectroscopy reveals typical redshifts of z~2 for this population, implying rest-frame mid-IR luminosities of L(6um)~6e46 erg/s and bolometric luminosities that can exceed L(bol)~1e14 L(sun). The corresponding intrinsic, unobscured hard X-ray luminosities are L(2-10)~4e45 erg/s for typical quasar templates. These are amongst the most luminous AGN known, though the optical spectra rarely show evidence of a broad-line region and the selection criteria imply heavy obscuration even at rest-frame 1.5um. We designed our X-ray observations to obtain robust detections for gas column densities N(H)<1e24 /cm2. In fact, the sources prove to be fainter than these predictions. Two of the sources were observed by both NuSTAR and XMM-Newton, with neither being detected by NuSTAR and one being faintly detected by XMM-Newton. A third source was observed only with XMM-Newton, yielding a faint detection. The X-ray data require gas column densities N(H)>1e24 /cm2, implying the sources are extremely obscured, consistent with Compton-thick, luminous quasars. The discovery of a significant population of heavily obscured, extremely luminous AGN does not conform to the standard paradigm of a receding torus, in which more luminous quasars are less likely to be obscured. If a larger sample conforms with this finding, then this suggests an additional source of obscuration for these extreme sources.

We report on XMM-Newton and Chandra observations of the Galactic supernova remnant candidate G38.7-1.4, together with complementary radio, infrared, and gamma-ray data. An approximately elliptical X-ray structure is found to be well correlated with radio shell as seen by the Very Large Array. The X-ray spectrum of G38.7-1.4 can be well-described by an absorbed collisional ionization equilibrium plasma model, which suggests the plasma is shock heated. Based on the morphology and the spectral behaviour, we suggest that G38.7-1.4 is indeed a supernova remnant belongs to a mix-morphology category.

Fitting the continuum component of a quasar spectrum in UV/optical band is challenging due to contamination of numerous emission lines. Traditional fitting algorithms such as the least-square fitting and the Levenberg-Marquardt algorithm (LMA) are fast but are sensitive to initial values of fitting parameters. They cannot guarantee to find global optimum solutions when the object functions have multiple minima. In this work, we attempt to fit a typical quasar spectrum using the Covariance Matrix Adaptation Evolution Strategy (CMA-ES). The spectrum is generated by composing a number of real quasar spectra from the Sloan Digital Sky Survey (SDSS) quasar catalog data release 3 (DR3) so it has a higher signal-to-noise ratio. The CMA-ES algorithm is an evolutionary algorithm that is designed to find the global rather than the local minima. The algorithm we implemented achieves an improved fitting result than the LMA and unlike the LMA, it is independent of initial parameter values. We are looking forward to implementing this algorithm to real quasar spectra in UV/optical band.

We have investigated the X-ray and $\gamma$-ray properties of the redback millisecond pulsar PSR J1723-2837 with XMM-Newton, Chandra and Fermi. We have discovered the X-ray orbital modulation of this binary system with the minimum that coincides with the phases of radio eclipse. The X-ray emission is clearly non-thermal in nature which can be well described by a simple power-law with a photon index of $\sim1.2$. The phase-averaged luminosity is $\sim9\times10^{31}$ erg/s in 0.3-10 keV which consumes $\sim0.2\%$ of the spin-down power. We have detected the $\gamma-$ray emission in $0.1-300$ GeV from this system at a significance of $\sim6\sigma$ for the first time. The $\gamma-$rays in this energy range consumes $\sim2\%$ of the spin-down power and can be modeled by a power-law with a photon index of $\sim2.6$. We discuss the high energy properties of the new redback in the context of a intrabinary shock model.

We provide the Chandra source list for the last ~quarter of the area covered by the Galactic Bulge Survey (GBS). The GBS targets two strips of 6\degr x 1\degr (12 square degrees in total), one above (1\degr<b<2\degr) and one below (-2\degr<b<-1\degr) the Galactic plane in the direction of the Galactic Center at X-ray, optical and near-infrared wavelengths. For the X-ray part of the survey we use 2 ks per Chandra pointing. We find 424 X-ray sources in the 63 Chandra observations we report on here. These sources are in addition to the 1216 X-ray sources discovered in the first part of the GBS survey described before. We discuss the characteristics and the X-ray variability of the brightest of the sources as well as the radio properties from existing radio surveys. We point out an interesting asymmetry in the number of X-ray sources as a function of their Galactic l and b coordinates which is probably caused by differences in average extinction towards the different parts of the GBS survey area.

The GeV$-$TeV $\gamma-$ray line signal is the smoking gun signature of the dark matter annihilation or decay. The detection of such a signal is one of the main targets of some space-based telescopes, including Fermi-LAT and the upcoming CALET, DAMPE and Gamma-400. An important feature of the dark-matter-annihilation originated $\gamma-$ray line photons is their concentration at the center of the Galaxy. So far no reliable $\gamma-$ray line has been detected by Fermi-LAT and the upper limits on the cross section of annihilation into $\gamma-$rays have been reported. We use these upper limits to estimate the "maximal" number of $\gamma-$ray line photons detectable for Fermi-LAT, DAMPE and Gamma-400 and then investigate the spatial distribution of these photons. We show that usually the center of the distribution region will be offset from the Galactic centre (Sgr A$^{\star}$) due to the limited statistics. Such a result is almost independent of the dark matter distribution models and renders the reconstruction of the dark matter distribution with the $\gamma-$ray line signal very challenging for the foreseeable space-based detectors.

On behalf of the WISE Science team, we present the discovery of a class of distant dust-enshrouded galaxies with extremely high luminosity. These galaxies are selected to have extreme red colors in the mid-IR using NASA's Wide-field Infrared Survey Explorer (WISE). They are faint in the optical and near-IR, predominantly at z=2-4, and with IR luminosity > $10^{13}\, L_{Sun}$, making them Hyper-Luminous Infrared Galaxies (HyLIRGs). SEDs incorporating the WISE, Spitzer, and Herschel PACS and SPIRE photometry indicate hot dust dominates the bolometric luminosity, presumably powered by AGN. Preliminary multi-wavelength follow-up suggests that they are different from normal populations in the local M-sigma relation. Their low source density implies that these objects are either intrinsically rare, or a short-lived phase in a more numerous population. If the latter is the case, these hot, dust-enshrouded galaxies may be an early stage in the interplay between AGN and galaxies.

Context. The international whole earth blazar telescope (WEBT) consortium planned and carried out three days of intensive micro-variability observations of S5 0716+714 from February 22, 2009 to February 25, 2009. This object was chosen due to its bright apparent magnitude range, its high declination, and its very large duty cycle for micro-variations. Aims. We report here on the long continuous optical micro-variability light curve of 0716+714 obtained during the multi-site observing campaign during which the Blazar showed almost constant variability over a 0.5 magnitude range. The resulting light curve is presented here for the first time. Observations from participating observatories were corrected for instrumental differences and combined to construct the overall smoothed light curve. Methods. Thirty-six observatories in sixteen countries participated in this continuous monitoring program and twenty of them submitted data for compilation into a continuous light curve. The light curve was analyzed using several techniques including Fourier transform, Wavelet and noise analysis techniques. Those results led us to model the light curve by attributing the variations to a series of synchrotron pulses. Results. We have interpreted the observed microvariations in this extended light curve in terms of a new model consisting of individual stochastic pulses due to cells in a turbulent jet which are energized by a passing shock and cool by means of synchrotron emission. We obtained an excellent fit to the 72-hour light curve with the synchrotron pulse model.

We present radio continuum mapping, optical imaging and spectroscopy of the newly discovered double-peaked broad-lined AGN WISE J233237.05-505643.5 at redshift z = 0.3447. This source exhibits an FR-I and FR-II hybrid-morphology, characterized by bright core, jet, and Doppler-boosted lobe structures in ATCA continuum maps at 1.5, 5.6, and 9 GHz. Unlike most FR-II objects, W2332-5056 is hosted by a disk-like galaxy. The core has a projected 5" linear radio feature that is perpendicular to the curved primary jet, hinting at unusual and complex activity within the inner 25 kpc. The multi-epoch optical-near-IR photometric measurements indicate significant variability over a 3-20 year baseline from the AGN component. Gemini-South optical data shows an unusual double-peaked emission-line features: the centroids of the broad-lined components of H-alpha and H-beta are blueshifted with respect to the narrow lines and host galaxy by ~ 3800 km/s. We examine possible cases which involve single or double supermassive black holes in the system, and discuss required future investigations to disentangle the mystery nature of this system.

We present WISE All-Sky mid-infrared (IR) survey detections of 55% (17/31) of the known QSOs at z>6 from a range of surveys: the SDSS, the CFHT-LS, FIRST, Spitzer and UKIDSS. The WISE catalog thus provides a substantial increase in the quantity of IR data available for these sources: 17 are detected in the WISE W1 (3.4-micron) band, 16 in W2 (4.6-micron), 3 in W3 (12-micron) and 0 in W4 (22-micron). This is particularly important with Spitzer in its warm-mission phase and no faint follow-up capability at wavelengths longwards of 5 microns until the launch of JWST. WISE thus provides a useful tool for understanding QSOs found in forthcoming large-area optical/IR sky surveys, using PanSTARRS, SkyMapper, VISTA, DES and LSST. The rest-UV properties of the WISE-detected and the WISE-non-detected samples differ: the detections have brighter i/z-band magnitudes and redder rest-UV colors. This suggests that a more aggressive hunt for very-high-redshift QSOs, by combining WISE W1 and W2 data with red observed optical colors could be effective at least for a subset of dusty candidate QSOs. Stacking the WISE images of the WISE-non-detected QSOs indicates that they are on average significantly fainter than the WISE-detected examples, and are thus not narrowly missing detection in the WISE catalog. The WISE-catalog detection of three of our sample in the W3 band indicates that their mid-IR flux can be detected individually, although there is no stacked W3 detection of sources detected in W1 but not W3. Stacking analyses of WISE data for large AGN samples will be a useful tool, and high-redshift QSOs of all types will be easy targets for JWST.

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

Stern et al.(2012) presented a study of WISE selection of AGN in the 2 deg^2 COSMOS field, finding that a simple criterion W1-W2>=0.8 provides a highly reliable and complete AGN sample for W2<15.05, where the W1 and W2 passbands are centered at 3.4 and 4.6 microns, respectively. Here we extend this study using the larger 9 deg^2 NOAO Deep Wide-Field Survey Bootes field which also has considerably deeper WISE observations than the COSMOS field, and find that this simple color-cut significantly loses reliability at fainter fluxes. We define a modified selection criterion combining the W1-W2 color and the W2 magnitude to provide highly reliable or highly complete AGN samples for fainter WISE sources. In particular, we define a color-magnitude cut that finds 130+/-4 deg^-2 AGN candidates for W2<17.11 with 90% reliability. Using the extensive UV through mid-IR broad-band photometry available in this field, we study the spectral energy distributions of WISE AGN candidates. As expected, the WISE AGN selection is biased towards objects where the AGN dominates the bolometric luminosity output, and that it can identify highly obscured AGN. We study the distribution of reddening in the AGN sample and discuss a formalism to account for sample incompleteness based on the step-wise maximum-likelihood method of Efstathiou et al.(1988). The resulting dust obscuration distributions depend strongly on AGN luminosity, consistent with the trend expected for a Simpson (2005) receding torus. At L_AGN~3x10^44 erg/s, 29+/-7% of AGN are observed as Type 1, while at ~4x10^45 erg/s the fraction is 64+/-13%. The distribution of obscuration values suggests that dust in the torus is present as both a diffuse medium and in optically thick clouds.

We present a new way to solve the platform deformation problem of co-planar interferometers. The platform of a co-planar interferometer can be deformed due to driving forces and gravity. A deformed platform will induce extra components into the geometric delay of each baseline, and change the phases of observed visibilities. The reconstructed images will also be diluted due to the errors of the phases. The platform deformations of The Yuan-Tseh Lee Array for Microwave Background Anisotropy (AMiBA) were modelled based on photogrammetry data with about 20 mount pointing positions. We then used the differential optical pointing error between two optical telescopes to fit the model parameters in the entire horizontal coordinate space. With the platform deformation model, we can predict the errors of the geometric phase delays due to platform deformation with given azimuth and elevation of the targets and calibrators. After correcting the phases of the radio point sources in the AMiBA interferometric data, we recover 50% - 70% flux loss due to phase errors. This allows us to restore more than 90% of a source flux. The method outlined in this work is not only applicable to the correction of deformation for other co-planar telescopes but also to single dish telescopes with deformation problems. This work also forms the basis of the upcoming science results of AMiBA-13.

We present the results of high spatial resolution observations at 1.1 mm waveband, with the Submillimetre Array (SMA), towards the protocluster G10.6-0.4. The 1.1 mm continuum emission reveals seven dense cores, towards which infall motions are all detected with the red-shifted absorption dips in HCN (3--2) line. This is the first time that infall is seen towards multiple sources in a protocluster. We also identified four infrared point sources in this region, which are most likely Class 0/I protostars. Two jet-like structures are also identified from Spitzer/IRAC image. The dense core located in the centre has much larger mass than the off-centre cores. The clump is in overall collapse and the infall motion is supersonic. The standard deviation of core velocities and the velocity differences between the cores and the cloud/clump are all larger than the thermal velocity dispersion. The picture of G10.6-0.4 seems to favor the "competitive accretion" model but needs to be tested by further observations.

We report the discovery of X-ray periodicity of ~265.3 ms from a deep XMM-Newton observation of the radio-quiet gamma-ray pulsar, PSR J2021+4026, located at the edge of the supernova remnant G78.2+2.1 (gamma-Cygni).The detected frequency is consistent with the gamma-ray pulsation determined by the observation of Fermi Gamma-ray Space Telescope at the same epoch. The X-ray pulse profile resembles the modulation of hot spot on the surface of the neutron star. The phase-averaged spectral analysis also suggests that the majority of the observed X-rays have a thermal origin. This is the third member in the class of radio-quiet pulsars with the significant pulsations detected from both X-rays and gamma-ray regimes.

The satellite-borne experiment PAMELA has been used to make new measurements of cosmic ray H and He isotopes. The isotopic composition was measured between 100 and 600 MeV/n for hydrogen and between 100 and 900 MeV/n for helium isotopes over the 23rd solar minimum from July 2006 to December 2007. The energy spectrum of these components carries fundamental information regarding the propagation of cosmic rays in the galaxy which are competitive with those obtained from other secondary to primary measurements such as B/C.

The ratios of fluxes of heavy nuclei from sulfur (Z=16) to chromium (Z=24) to the flux of iron were measured by the ATIC-2 experiment. The ratios are decreasing functions of energy from 5 GeV/n to approximately 80 GeV/n, as expected. However, an unexpected sharp upturn in the ratios are observed for energies above 100 GeV/n for all elements from Z=16 to Z=24. Similar upturn but with lower amplitude was also discovered in the ATIC-2 data for the ratio of fluxes of abundant even nuclei (C, O, Ne, Mg, Si) to the flux of iron. Therefore the spectrum of iron is significantly different from the spectra of other abundant even nuclei.

Fe exists abundantly in the universe. In particular, the dynamical structures and transport properties of warm dense Fe are crucial for understanding the evolution and structures of giant planets. In this article, we present the ionic structures, equation of states, diffusion and viscosity of Fe at two typical densities of 33.385 g/cm$^3$ and 45 g/cm$^3$ in the temperature range of 1 eV and 10 eV, giving the data by the first principles calculations using quantum Langevin molecular dynamics (QLMD). Furthermore, the validation of Stokes-Einstein (SE) relation in this regime is discussed, showing the importance of choosing the effective atomic diameter. The results remind us of the careful usage of the SE relation under extreme conditions.

We study the perturbation of the effective Hamiltonian constraint with holonomy correction from Euclidean loop quantum gravity. The Poisson bracket between the corrected Hamiltonian constraint and the diffeomorphism constraint is derived for vector modes. Some specific form of the holonomy correction function $f^{i}_{cd}$ is found, which satisfies that the constraint algebra is anomaly-free. This result confirms the possibility of non-trivial holonomy corrections from full theory while preserving anomaly-free constraint algebra in the perturbation framework. It also gives valuable hints on the possible form of holonomy corrections in effective loop quantum gravity.

Using the data from the Large Area Telescope on board the Fermi Gamma-ray Space Telescope, we have searched for the gamma-ray pulsations from the direction of globular cluster M28 (NGC 6626). We report the discovery of a signal with the frequency consistent with that of the energetic millisecond pulsar (MSP) PSR B1821-24 in M28. A weighted H-test test statisic (TS) of 28.8 is attained which corresponds to a chance probability of ~1e-5 (4.3-sigma detection). With a phase-resolved analysis, the pulsed component is found to contribute ~25% of the total observed gamma-ray emission from the cluster. On the other hand, the unpulsed level provides a constraint for the underlying MSP population and the fundamental plane relations for the scenario of inverse Compton scattering. Follow-up timing observations in radio/X-ray are encouraged for further investigating this periodic signal candidate.

The energy spectra of galactic cosmic rays carry fundamental information regarding their origin and propagation. These spectra, when measured near Earth, are significantly affected by the solar magnetic field. A comprehensive description of the cosmic radiation must therefore include the transport and modulation of cosmic rays inside the heliosphere. During the end of the last decade the Sun underwent a peculiarly long quiet phase well suited to study modulation processes. In this paper we present proton spectra measured from July 2006 to December 2009 by PAMELA. The large collected statistics of protons allowed the time variation to be followed on a nearly monthly basis down to 400 MV. Data are compared with a state-of-the-art three-dimensional model of solar modulation.

(Abridged) We present a systematic study of the X-ray and multiwavelength properties of a sample of 17 highly radio-loud quasars (HRLQs) at z > 4 with sensitive X-ray coverage from new Chandra and archival Chandra, XMM-Newton, and Swift observations. Eight of the new and archival observations are reported in this work for the first time. New Chandra observations of two moderately radio-loud and highly optically luminous quasars at z > 4 are also reported. Our HRLQ sample represents the top ~5% of radio-loud quasars in terms of radio loudness. We found that our HRLQs have an X-ray emission enhancement over HRLQs at lower redshifts (by a typical factor of ~3), and this effect, after controlling for several factors which may introduce biases, has been solidly estimated to be significant at the 3-4 sigma level. HRLQs at z=3-4 are also found to have a similar X-ray emission enhancement over z < 3 HRLQs, which supports further the robustness of our results. We discuss models for the X-ray enhancement's origin including a fractional contribution from inverse Compton scattering of cosmic microwave background photons. No strong correlations are found between the relative X-ray brightness and optical/UV emission-line rest-frame equivalent widths (REWs) for radio-loud quasars. However, the line REWs are positively correlated with radio loudness, which suggests that relativistic jets make a negligible contribution to the optical/UV continua of these HRLQs (contrary to the case where the emission lines are diluted by the relativistically boosted continuum). Our HRLQs are generally consistent with the known anti-correlation between radio loudness and X-ray power-law photon index. We also found that the two moderately radio-loud quasars appear to have the hardest X-ray spectra among our objects, suggesting that intrinsic X-ray absorption (N_H~10^23 cm^-2) may be present.

We have monitored the BL Lac object S5 0716+714 in five intermediate optical wavebands from 2004 September to 2011 April. Here we present the data that include 8661 measurements which represents one of the largest databases obtained for an object at optical domain. A simple analysis of the data indicates that the object was active in most time, and intraday variability was frequently observed. In total, the object varied by 2.614 magnitudes in the $i$ band. Strong bluer-when-brighter chromatism was observed on long, intermediate, and short timescales.

We present X-ray and multiwavelength studies of a sample of eight high-luminosity active galactic nuclei (AGNs) with disc-like H\beta emission-line profiles selected from the Sloan Digital Sky Survey Data Release 7. These sources have higher redshift (z~0.6) than the majority of the known disc-like emitters, and they occupy a largely unexplored space in the luminosity-redshift plane. Seven sources have typical AGN X-ray spectra with power-law photon indices of \Gamma~1.4-2.0; two of them show some X-ray absorption (column density N_H~10^{21}-10^{22} cm^{-2}$ for neutral gas). The other source, J0850+4451, has only three hard X-ray photons detected and is probably heavily obscured (N_H>3x10^{23} cm^{-2}). This object is also identified as a low-ionization broad absorption line (BAL) quasar based on Mg II \lambda2799 absorption; it is the first disc-like emitter reported that is also a BAL quasar. The IR-to-UV spectral energy distributions (SEDs) of these eight sources are similar to the mean SEDs of typical quasars with a UV "bump", suggestive of standard accretion discs radiating with high efficiency, which differs from low-luminosity disc-like emitters. Studies of the X-ray-to-optical power-law slope parameters (\alpha_{OX}) indicate that there is no significant excess X-ray emission in these high-luminosity disc-like emitters. Energy budget analysis suggests that for disc-like emitters in general, the inner disc must illuminate and ionize the outer disc efficiently (~15% of the nuclear ionizing radiation is required on average) via direct illumination and/or scattering. Warped accretion discs are probably needed for direct illumination to work in high-luminosity objects, as their geometrically thin inner discs decrease the amount of direct illumination possible for a flat disc.

We present X-ray and multiwavelength studies of a sample of eight high-luminosity active galactic nuclei (AGNs) with disc-like H\beta emission-line profiles selected from the Sloan Digital Sky Survey Data Release 7. These sources have higher redshift (z~0.6) than the majority of the known disc-like emitters, and they occupy a largely unexplored space in the luminosity-redshift plane. Seven sources have typical AGN X-ray spectra with power-law photon indices of \Gamma~1.4-2.0; two of them show some X-ray absorption (column density N_H~10^{21}-10^{22} cm^{-2}$ for neutral gas). The other source, J0850+4451, has only three hard X-ray photons detected and is probably heavily obscured (N_H>3x10^{23} cm^{-2}). This object is also identified as a low-ionization broad absorption line (BAL) quasar based on Mg II \lambda2799 absorption; it is the first disc-like emitter reported that is also a BAL quasar. The IR-to-UV spectral energy distributions (SEDs) of these eight sources are similar to the mean SEDs of typical quasars with a UV "bump", suggestive of standard accretion discs radiating with high efficiency, which differs from low-luminosity disc-like emitters. Studies of the X-ray-to-optical power-law slope parameters (\alpha_{OX}) indicate that there is no significant excess X-ray emission in these high-luminosity disc-like emitters. Energy budget analysis suggests that for disc-like emitters in general, the inner disc must illuminate and ionize the outer disc efficiently (~15% of the nuclear ionizing radiation is required on average) via direct illumination and/or scattering. Warped accretion discs are probably needed for direct illumination to work in high-luminosity objects, as their geometrically thin inner discs decrease the amount of direct illumination possible for a flat disc.

In this paper, we presented metal abundance properties of 144 M81 globular clusters. These globulars consist of the largest globular cluster sample in M81 till now. Our main results are: the distribution of metallicities are bimodal, with metallicity peaks at [Fe/H]\sim-1.51 and -0.58, and the metal-poor globular clusters tend to be less spatially concentrated than the metal-rich ones; the metal-rich globular clusters in M81 do not demonstrate a centrally concentrated spatial distribution as the metal-rich ones in M31 do; like our Galaxy and M31, the globular clusters in M81 have a small radial metallicity gradient. These results are consistent with those obtained based on a small sample of M81 globular clusters. In addition, this paper showed that there is evidence that a strong rotation of the M81 globular cluster system around the minor axis exists, and that rotation is present in the metal-rich globular cluster subsample, while the metal-poor globular cluster subsample shows no evidence for rotation. The most significant difference between the rotation of the metal-rich and metal-poor globular clusters occurs at intermediate projected galactocentric radii. The results of this paper confirm the conclusion of Schroder et al. that M81's metal-rich globular clusters at intermediate projected radii were associated with a thick disk of M81.

We report on evidence for orbital phase-dependence of the gamma-ray emission from PSR B1957+20 black widow system by using the data of the Fermi Large Area Telescope. We divide an orbital cycle into two regions: a region containing the inferior conjunction, and the other region containing rest of the orbital cycle. We show that the observed spectra for the different orbital regions are fitted by different functional forms. The spectrum of the orbital region containing inferior conjunction can be described by a power-law with an exponential cutoff (PLE) model, which gives the best-fit model for the orbital phase that does not contain the inferior conjunction, plus an extra component above ~2.7 GeV. The emission above 3 GeV in this region is detected with a ~7-sigma confidence level. The gamma-ray data above ~2.7 GeV are observed to be modulated at the orbital period at the ~2.3-sigma level. We anticipate that the PLE component dominating below ~2.7 GeV originates from the pulsar magnetosphere. We also show that the inverse-Compton scattering of the thermal radiation of the companion star off a "cold" ultra-relativistic pulsar wind can explain the extra component above ~2.7 GeV. The black widow pulsar PSR B1957+20 may be the member of a new class of object, in the sense that the system is showing gamma-ray emission with both magnetospheric and pulsar wind origins.

Anomalous x-ray pulsars (AXPs) are thought to be magnetars which are young isolated neutron stars with extremely strong magnetic fields of > 10^14Gauss. Their tremendous magnetic fields inferred from the spin parameters provide a huge energy reservoir to power the observed x-ray emission. High-energy emission above 0.3 MeV has never been detected despite intensive search. Here, we present the possible Fermi Large Area Telescope (LAT) detection of {\gamma}-ray pulsations above 200 MeV from the AXP, 1E 2259+586, which puts the current theoretical models of {\gamma}-ray emission mechanisms of magnetars into challenge. We speculate that the high-energy {\gamma}-rays originate from the outer magnetosphere of the magnetar.

We present Spitzer 3.6 and 4.5 $\mu$m photometry and positions for a sample of 1510 brown dwarf candidates identified by the WISE all-sky survey. Of these, 166 have been spectroscopically classified as objects with spectral types M(1), L(7), T(146), and Y(12); Sixteen other objects are non-(sub)stellar in nature. The remainder are most likely distant L and T dwarfs lacking spectroscopic verification, other Y dwarf candidates still awaiting follow-up, and assorted other objects whose Spitzer photometry reveals them to be background sources. We present a catalog of Spitzer photometry for all astrophysical sources identified in these fields and use this catalog to identify 7 fainter (4.5 $\mu$m $\sim$ 17.0 mag) brown dwarf candidates, which are possibly wide-field companions to the original WISE sources. To test this hypothesis, we use a sample of 919 Spitzer observations around WISE-selected high-redshift hyper-luminous infrared galaxy (HyLIRG) candidates. For this control sample we find another 6 brown dwarf candidates, suggesting that the 7 companion candidates are not physically associated. In fact, only one of these 7 Spitzer brown dwarf candidates has a photometric distance estimate consistent with being a companion to the WISE brown dwarf candidate. Other than this there is no evidence for any widely separated ($>$ 20 AU) ultra-cool binaries. As an adjunct to this paper, we make available a source catalog of $\sim$ 7.33 $\times 10^5$ objects detected in all of these Spitzer follow-up fields for use by the astronomical community. The complete catalog includes the Spitzer 3.6 and 4.5 $\mu$m photometry, along with positionally matched $B$ and $R$ photometry from USNO-B; $J$, $H$, and $K_s$ photometry from 2MASS; and $W1$, $W2$, $W3$, and $W4$ photometry from the WISE all-sky catalog.

We report the discovery by the Wide-field Infrared Survey Explorer of the z = 2.452 source WISE J181417.29+341224.9, the first hyperluminous source found in the WISE survey. WISE 1814+3412 is also the prototype for an all-sky sample of ~1000 extremely luminous "W1W2-dropouts" (sources faint or undetected by WISE at 3.4 and 4.6 microns and well detected at 12 or 22 microns). The WISE data and a 350 micron detection give a minimum bolometric luminosity of 3.7 x 10^13 Lsun, with ~10^14 Lsun plausible. Followup images reveal four nearby sources: a QSO and two Lyman Break Galaxies (LBGs) at z = 2.45, and an M dwarf star. The brighter LBG dominates the bolometric emission. Gravitational lensing is unlikely given the source locations and their different spectra and colors. The dominant LBG spectrum indicates a star formation rate ~300 Msun/yr, accounting for < 10% of the bolometric luminosity. Strong 22 micron emission relative to 350 microns implies that warm dust contributes significantly to the luminosity, while cooler dust normally associated with starbursts is constrained by an upper limit at 1.1 mm. Radio emission is ~10x above the far-infrared/radio correlation, indicating an active galactic nucleus is present. An obscured AGN combined with starburst and evolved stellar components can account for the observations. If the black hole mass follows the local M_BH-bulge mass relation, the implied Eddington ratio is >~4. WISE 1814+3412 may be a heavily obscured object where the peak AGN activity occurred prior to the peak era of star formation.

We have used the Caltech Submillimeter Observatory (CSO) to follow-up a sample of WISE-selected, hyperluminous galaxies, so called W1W2-dropout galaxies. This is a rare (~ 1000 all-sky) population of galaxies at high redshift (peaks at z=2-3), that are faint or undetected by WISE at 3.4 and 4.6 um, yet are clearly detected at 12 and 22 um. The optical spectra of most of these galaxies show significant AGN activity. We observed 14 high-redshift (z > 1.7) W1W2-dropout galaxies with SHARC-II at 350 to 850 um, with 9 detections; and observed 18 with Bolocam at 1.1 mm, with five detections. Warm Spitzer follow-up of 25 targets at 3.6 and 4.5 um, as well as optical spectra of 12 targets are also presented in the paper. Combining WISE data with observations from warm Spitzer and CSO, we constructed their mid-IR to millimeter spectral energy distributions (SEDs). These SEDs have a consistent shape, showing significantly higher mid-IR to submm ratios than other galaxy templates, suggesting a hotter dust temperature. We estimate their dust temperatures to be 60-120 K using a single-temperature model. Their infrared luminosities are well over 10^{13} Lsun. These SEDs are not well fitted with existing galaxy templates, suggesting they are a new population with very high luminosity and hot dust. They are likely among the most luminous galaxies in the Universe. We argue that they are extreme cases of luminous, hot dust-obscured galaxies (DOGs), possibly representing a short evolutionary phase during galaxy merging and evolution. A better understanding of their long-wavelength properties needs ALMA as well as Herschel data.

We present simultaneous Planck, Swift, Fermi, and ground-based data for 105 blazars belonging to three samples with flux limits in the soft X-ray, hard X-ray, and gamma-ray bands. Our unique data set has allowed us to demonstrate that the selection method strongly influences the results, producing biases that cannot be ignored. Almost all the BL Lac objects have been detected by Fermi-LAT, whereas ~40% of the flat-spectrum radio quasars (FSRQs) in the radio, soft X-ray, and hard X-ray selected samples are still below the gamma-ray detection limit even after integrating 27 months of Fermi-LAT data. The radio to sub-mm spectral slope of blazars is quite flat up to ~70GHz, above which it steepens to <\alpha>~-0.65. BL Lacs have significantly flatter spectra than FSRQs at higher frequencies. The distribution of the rest-frame synchrotron peak frequency (\nupS) in the SED of FSRQs is the same in all the blazar samples with <\nupS>=10^13.1 Hz, while the mean inverse-Compton peak frequency, <\nupIC>, ranges from 10^21 to 10^22 Hz. The distributions of \nupS and of \nupIC of BL Lacs are much broader and are shifted to higher energies than those of FSRQs and strongly depend on the selection method. The Compton dominance of blazars ranges from ~0.2 to ~100, with only FSRQs reaching values >3. Its distribution is broad and depends strongly on the selection method, with gamma-ray selected blazars peaking at ~7 or more, and radio-selected blazars at values ~1, thus implying that the assumption that the blazar power is dominated by high-energy emission is a selection effect. Simple SSC models cannot explain the SEDs of most of the gamma-ray detected blazars in all samples. The SED of the blazars that were not detected by Fermi-LAT may instead be consistent with SSC emission. Our data challenge the correlation between bolometric luminosity and \nupS predicted by the blazar sequence.

Studying the acceleration and propagation mechanisms of Galactic cosmic rays can provide information regarding astrophysical sources, the properties of our Galaxy, and possible exotic sources such as dark matter. To understand cosmic ray acceleration and propagation mechanisms, accurate measurements of different cosmic ray elements over a wide energy range are needed. The PAMELA experiment is a satellite-borne apparatus which allows different cosmic ray species to be identified over background. Measurements of the cosmic ray antiproton flux and the antiproton-to-proton flux ratio from 1.5 GeV to 180 GeV are presented in this thesis. Compared to previous experiments, PAMELA extends the energy range of antiproton measurements and provides significantly higher statistics. The derived antiproton flux and antiproton-to-proton flux ratio are consistent with previous measurements and generally considered to be produced as secondary products when cosmic ray protons and helium nuclei interact with the interstellar medium. To constrain cosmic ray acceleration and propagation models, the antiproton data measured by PAMELA were further used together with the proton spectrum reported by PAMELA, as well as the B/C data provided by other experiments. Statistical tools were interfaced with the cosmic ray propagation package GALPROP to perform the constraining analyses. Diffusion models with a linear diffusion coefficient and modified diffusion models with a low energy dependence of the diffusion coefficient were studied in the $\chi^{2}$ study. Uncertainties on the parameters and the goodness of fit of each model were given. Some models are further studied using the Bayesian inference. Posterior means and errors of the parameters base on our prior knowledge on them were obtained in the Bayesian framework. This method also allowed us to understand the correlation between parameters and compare models.

We report a new technique to select 1.6<z<4.6 dusty Lyman-alpha emitters (LAEs), over a third of which are `blobs' (LABs) with emission extended on scales of 30-100kpc. Combining data from the NASA Wide-field Infrared Survey Explorer (WISE) mission with optical spectroscopy from the W.M. Keck telescope, we present a color criteria that yields a 78% success rate in identifying rare, dusty LAEs of which at least 37% are LABs. The objects have a surface density of only ~0.1 per square degree, making them rare enough that they have been largely missed in narrow surveys. We measured spectroscopic redshifts for 92 of these WISE-selected, typically radio-quiet galaxies and find that the LAEs (LABs) have a median redshift of 2.3 (2.5). The WISE photometry coupled with data from Herschel reveals that these galaxies have extreme far-infrared luminosities (L_IR>10^{13-14}L_sun) and warm colors, typically larger than submillimeter-selected galaxies (SMGs) and dust-obscured galaxies (DOGs). These traits are commonly associated with the dust being energized by intense AGN activity. We hypothesize that the combination of spatially extended Lyman-alpha, large amounts of warm IR-luminous dust, and rarity (implying a short-lived phase) can be explained if the galaxies are undergoing strong `feedback' transforming them from an extreme dusty starburst to a QSO.

We present a method for attaining sub-arcsecond pointing stability during sub- orbital balloon flights, as designed for in the High Altitude Lensing Observatory (HALO) concept. The pointing method presented here has the potential to perform near-space quality optical astronomical imaging at 1-2% of the cost of space-based missions. We also discuss an architecture that can achieve sufficient thermomechanical stability to match the pointing stability. This concept is motivated by advances in the development and testing of Ultra Long Duration Balloon (ULDB) flights which promise to allow observation campaigns lasting more than three months. The design incorporates a multi-stage pointing architecture comprising: a gondola coarse azimuth control system, a multi-axis nested gimbal frame structure with arcsecond stability, a telescope de-rotator to eliminate field rotation, and a fine guidance stage consisting of both a telescope mounted angular rate sensor and guide CCDs in the focal plane to drive a fast-steering mirror. We discuss the results of pointing tests together with a preliminary thermo-mechanical analysis required for sub-arcsecond pointing at high altitude. Possible future applications in the areas of wide-field surveys and exoplanet searches are also discussed.

We have compiled a catalog of optically-selected quasars with simultaneous observations in UV/optical and X-ray bands by the Swift Gamma Ray Burst Explorer. Objects in this catalog are identified by matching the Swift pointings with the Sloan Digital Sky Survey Data Release 5 quasar catalog. The final catalog contains 843 objects, among which 637 have both UVOT and XRT observations and 354 of which are detected by both instruments. The overall X-ray detection rate is ~60% which rises to ~85% among sources with at least 10 ks of XRT exposure time. We construct the time-averaged spectral energy distribution for each of the 354 quasars using UVOT photometric measurements and XRT spectra. From model fits to these SEDs, we find that the big blue bump contributes about 0.3 dex to the quasar luminosity. We re-visit the alpha_ox-L_uv relation by selecting a clean sample with only type 1 radio-quiet quasars; the dispersion of this relation is reduced by at least 15% compared to studies that use non-simultaneous UV/optical and X-ray data. We only found a weak correlation between L/L_Edd and alpha_uv. We do not find significant correlations between alpha_x and alpha_ox, alpha_ox and alpha_uv, and alpha_x and Log L(0.3-10 keV). The correlations between alpha_uv and alpha_x, alpha_ox and alpha_x, alpha_ox and alpha_uv, L/L_Edd and alpha_x, and L/L_Edd and alpha_ox are stronger amongst low-redshift quasars, indicating that these correlations are likely driven by the changes of SED shape with accretion state.

This paper presents multicolor optical photometry of the nearby galaxy cluster Abell 119 (z = 0:0442) with the Beijing-Arizona-Taiwan-Connecticut (BATC) system of 15 intermediate bands. Within the BATC viewing field of 58'* 58', there are 368 galaxies with known spectroscopic redshifts, including 238 member galaxies (called sample I). Based on the spectral energy distributions (SEDs) of 1376 galaxies brighter than iBATC = 19:5, photometric redshift technique and the color-magnitude relation of earlytype galaxies are applied to select faint member galaxies. As a result, 117 faint galaxies were selected as new member galaxies. Combined with sample I, an enlarged sample (called sample II) of 355 member galaxies is obtained. Spatial distribution and localized velocity structure for two samples demonstrate that A119 is a dynamically complex cluster with at least three prominent substructures in the central region within 1 Mpc. A large velocity dispersion for the central clump indicates a merging along the line of sight. No significant evidences for morphology and luminosity segregations are found in both samples. With the evolutionary synthesis model PEGASE, environmental effect on the star formation properties is confirmed. Faint galaxies in low-density region tend to have longer time scales of star formation, smaller mean stellar ages, and lower metallicities of interstellar medium, which is in agreement with the context of hierarchical cosmological scenario.

We describe the infrared properties of sources detected over ~36 deg^2 of sky in the GAMA 15-hr equatorial field, using data from both the Herschel Astrophysical Terahertz Large-Area Survey (H-ATLAS) and Wide-field Infrared Survey (WISE). With 5-sigma point-source depths of 34 and 0.048 mJy at 250 micron and 3.4 micron, respectively, we are able to identify 50.6% of the H-ATLAS sources in the WISE survey, corresponding to a surface density of ~630 deg^{-2}. Approximately two-thirds of these sources have measured spectroscopic or optical/near-IR photometric redshifts of z<1. For sources with spectroscopic redshifts at z<0.3, we find a linear correlation between the infrared luminosity at 3.4 micron and that at 250 micron, with +-50% scatter over ~1.5 orders of magnitude in luminosity, ~10^9 - 10^{10.5} L_sun. By contrast, the matched sources without previously measured redshifts (r>~20.5) have 250-350 micron flux density ratios that suggest either high-redshift galaxies (z>~1.5) or optically faint low-redshift galaxies with unusually low temperatures (T<~20). Their small 3.4-250 micron flux ratios favor a high-redshift galaxy population, as only the most actively star-forming galaxies at low redshift (e.g., Arp 220) exhibit comparable flux density ratios. Furthermore, we find a relatively large AGN fraction (~30%) in a 12 micron flux-limited subsample of H-ATLAS sources, also consistent with there being a significant population of high-redshift sources in the no-redshift sample.

We present the results of our optical monitoring of the BL Lac object S5 0716+714 on seven nights in 2006 December. The monitoring was carried out simultaneously at three optical wavelengths with a novel photometric system. The object did not show large-amplitude internight variations during this period. Intranight variations were observed on four nights and probably on one more. Strong bluer-when-brighter chromatism was detected on both intranight and internight timescales. The intranight variation amplitude decreases in the wavelength sequence of B', R', and V'. Cross correlation analyses revealed that the variability at the $B'$ and $V'$ bands lead that at the $R'$ band by about 30 minutes on one night.

We report the discovery of pulsed {\gamma}-ray emission and X-ray emission from the black widow millisecond pulsar PSR J2051-0827 by using the data from the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope and the Advanced CCD Imaging Spectrometer array (ACIS-S) on the Chandra X-ray Observatory. Using 3 years of LAT data, PSR J2051-0827 is clearly detected in {\gamma}-ray with a signicance of \sim 8{\sigma} in the 0.2 - 20 GeV band. The 200 MeV - 20 GeV {\gamma}-ray spectrum of PSR J2051-0827 can be modeled by a simple power- law with a photon index of 2.46 \pm 0.15. Significant (\sim 5{\sigma}) {\gamma}-ray pulsations at the radio period were detected. PSR J2051-0827 was also detected in soft (0.3-7 keV) X-ray with Chandra. By comparing the observed {\gamma}-rays and X-rays with theoretical models, we suggest that the {\gamma}-ray emission is from the outer gap while the X-rays can be from intra-binary shock and pulsar magnetospheric synchrotron emissions.

We report on a population of X-ray weak quasars with similar UV emission-line properties to those of the remarkable quasar PHL 1811. All radio-quiet PHL 1811 analogs are notably X-ray weak by a mean factor of ~13, with hints of heavy X-ray absorption. Correlations between the X-ray weakness and UV emission-line properties suggest that PHL 1811 analogs may have extreme wind-dominated broad emission-line regions (BELRs). We propose an AGN geometry that can potentially unify the PHL 1811 analogs and the general population of weak-line quasars.

X-ray loud M dwarfs are a major source of by-product (contamination) in the X-ray band of the multiwavelength quasar survey (MWQS). As a by-product, the low dispersion spectra of 22 M dwarfs are obtained in which the spectra of 16 sources are taken for the first time. The spectral types and distance of the sample are given based on spectral indices CaH2, CaH3, and TiO5. The parameter {\zeta}TiO/CaH is calculated to make the metallicity class separation among dwarfs, subdwarfs and extreme subdwarfs. We also discuss the distributions in the diagrams of Log(Lx/Lbol) versus spectral type and infrared colors.

We use data from the Wide-Field Infrared Survey Explorer (WISE) to perform a statistical study on the mid-infrared (IR) properties of a large number ($\sim10^2$) of BL Lac objects --- low-luminosity Active Galactic Nuclei (AGN) with a jet beamed toward the Earth. As expected, many BL Lac objects are so highly beamed that their jet synchrotron emission dominates their IR spectral energy distributions. In other BL Lac objects, however, the jet is not strong enough to completely dilute the rest of the AGN emission. We do not see observational signatures of the dusty torus from these weakly beamed BL Lac objects. The lack of observable torus emission is consistent with suggestions that BL Lac objects are fed by radiatively inefficient accretion disks. Implications for the "nature vs. nurture" debate for FR I and FR II radio galaxies are briefly discussed. Our study supports the notion that, beyond orientation, accretion rate plays an important role in AGN unification.

Spectral energy distributions (SEDs) and radio continuum spectra are presented for a northern sample of 104 extragalactic radio sources, based on the Planck Early Release Compact Source Catalogue (ERCSC) and simultaneous multifrequency data. The nine Planck frequencies, from 30 to 857 GHz, are complemented by a set of simultaneous observations ranging from radio to gamma-rays. This is the first extensive frequency coverage in the radio and millimetre domains for an essentially complete sample of extragalactic radio sources, and it shows how the individual shocks, each in their own phase of development, shape the radio spectra as they move in the relativistic jet. The SEDs presented in this paper were fitted with second and third degree polynomials to estimate the frequencies of the synchrotron and inverse Compton (IC) peaks, and the spectral indices of low and high frequency radio data, including the Planck ERCSC data, were calculated. SED modelling methods are discussed, with an emphasis on proper, physical modelling of the synchrotron bump using multiplecomponents. Planck ERCSC data also suggest that the original accelerated electron energy spectrum could be much harder than commonly thought, with power-law index around 1.5 instead of the canonical 2.5. The implications of this are discussed for the acceleration mechanisms effective in blazar shock. Furthermore in many cases the Planck data indicate that gamma-ray emission must originate in the same shocks that produce the radio emission.

(Abridged) We report on the X-ray and multiwavelength properties of 11 radio-quiet quasars with weak or no emission lines identified by the Sloan Digital Sky Survey (SDSS) with redshift z=0.4-2.5. The distribution of relative X-ray brightness for our low-redshift weak-line quasar (WLQ) candidates is significantly different from that of typical radio-quiet quasars, having an excess of X-ray weak sources, but it is consistent with that of high-redshift WLQs. The X-ray weak sources generally show similar UV emission-line properties to those of the X-ray weak quasar PHL 1811; they may belong to the notable class of PHL 1811 analogs. The average X-ray spectrum of these sources is somewhat harder than that of typical radio-quiet quasars. Several other low-redshift WLQ candidates have normal ratios of X-ray-to-optical/UV flux, and their average X-ray spectral properties are also similar to those of typical radio-quiet quasars. The X-ray weak and X-ray normal WLQ candidates may belong to the same subset of quasars having high-ionization "shielding gas" covering most of the wind-dominated broad emission-line region, but be viewed at different inclinations. The mid-infrared-to-X-ray spectral energy distributions (SEDs) of these sources are generally consistent with those of typical SDSS quasars, showing that they are not likely to be BL Lac objects with relativistically boosted continua and diluted emission lines. However, one source in our X-ray observed sample is remarkably strong in X-rays, indicating that a small fraction of low-redshift WLQ candidates may actually be BL Lacs residing in the radio-faint tail of the BL Lac population. We also investigate universal selection criteria for WLQs over a wide range of redshift, finding that it is not possible to select WLQ candidates in a fully consistent way using different prominent emission lines as a function of redshift.

B514 is a remote M31 globular cluster which locating at a projected distance of R_p~55 kpc. Deep observations with the Advanced Camera for Surveys (ACS) on the Hubble Space Telescope (HST) are used to provide the accurate integrated light and star counts of B514. By coupling analysis of the distribution of the integrated light with star counts, we are able to reliably follow the profile of the cluster out to ~40". Based on the combined profile, we study in detail its surface brightness distribution in F606W and F814W filters, and determine its structural parameters by fitting a single-mass isotropic King model. The results showed that, the surface brightness distribution departs from the best-fit King model for r>10". B514 is quite flatted in the inner region, and has a larger half-light radius than majority of normal globular clusters of the same luminosity. It is interesting that, in the M_V versus log R_h plane, B514 lies nearly on the threshold for ordinary globular clusters as defined by Mackey & van den Bergh. In addition, B514 was observed as part of the Beijing-Arizona-Taiwan-Connecticut (BATC) Multicolor Sky Survey, using 13 intermediate-band filters covering a wavelength range of 3000--8500 \AA. Based on aperture photometry, we obtain its SEDs as defined by the 13 BATC filters. We determine the cluster's age and mass by comparing its SEDs (from 2267 to 20000{\AA}, comprising photometric data in the near-ultraviolet of GALEX, 5 SDSS bands, 13 BATC intermediate-band, and 2MASS near-infrared JHKs} filters) with theoretical stellar population synthesis models, resulting in age of $11.5\pm3.5$ Gyr. This age confirms the previous suggestion that B514 is an old GC in M31. B514 has a mass of $0.96-1.08 \times 10^6 \rm M_sun$, and is a medium-mass globular cluster in M31.

The accretion efficiency for individual black holes is very difficult to determine accurately. There are many factors that can influence each step of the calculation, such as the dust and host galaxy contribution to the observed luminosity, the black hole mass and more importantly, the uncertainties on the bolometric luminosity measurement. Ideally, we would measure the AGN emission at every wavelength, remove the host galaxy and dust, reconstruct the AGN spectral energy distribution and integrate to determine the intrinsic emission and the accretion rate. In reality, this is not possible due to observational limitations and our own galaxy line of sight obscuration. We have then to infer the bolometric luminosity from spectral measurements made in discontinuous wavebands and at different epochs. In this paper we tackle this issue by exploring different methods to determine the bolometric luminosity. We first explore the trend of accretion efficiency with black hole mass (efficiency proportional to M^{\sim 0.5}) found in recent work by Davis & Laor and discuss why this is most likely an artefact of the parameter space covered by their PG quasar sample. We then target small samples of AGN at different redshifts, luminosities and black hole masses to investigate the possible methods to calculate the accretion efficiency. For these sources we are able to determine the mass accretion rate and, with some assumptions, the accretion efficiency distributions. Even though we select the sources for which we are able to determine the parameters more accurately, there are still factors affecting the measurements that are hard to constrain. We suggest methods to overcome these problems based on contemporaneous multi-wavelength data measurements and specifically targeted observations for AGN in different black hole mass ranges.

We report the discovery of GeV emission at the position of supernova remnant Kes 17 by using the data from the Large Area Telescope on board the Fermi Gamma-ray Space Telescope. Kes 17 can be clearly detected with a significance of ~12 sigma in the 1 - 20 GeV range. Moreover, a number of gamma-ray sources were detected in its vicinity. The gamma-ray spectrum of Kes 17 can be well described by a simple power-law with a photon index of ~ 2.4. Together with the multi-wavelength evidence for its interactions with the nearby molecular cloud, the gamma-ray detection suggests that Kes 17 is a candidate acceleration site for cosmic-rays.

The colour behaviour of blazars is a subject of much debate. One argument is that the BL Lac objects show bluer-when-brighter chromatism while the flat-spectrum radio quasars (FSRQs) display redder-when-brighter trend. Base on a 3.5-year three-colour monitoring programme, we studied the optical variability and colour behaviour of one FSRQ, 3C 345. There is at least one outburst in this period. The overall variation amplitude is 2.640 mags in the $i$ band. Intra-night variability was observed on two nights. The bluer-when-brighter and redder-when-brighter chromatisms were simultaneously observed in this object when using different pairs of passbands to compute the colours. The bluer-when-brighter chromatism is a shared property with the BL Lacs, while the redder-when-brighter trend is likely due to two less variable emission features, the Mg\,{\sc ii} line and the blue bump, at short wavelengths. With numerical simulations, we show that some other strong but less variable emission lines in the spectrum of FSRQs may also significantly alter their colour behaviour. Then the colour behaviour of an FSRQ is linked not only to the emission process in the relativistic jet, but also to the redshift, the passbands used for computing the colour and the strengths of the less variable emission features relative to the strength of the non-thermal continuum.

We analyze the HCO+ 3-2 and H13CO+ 3-2 line profiles of 27 high-mass star-forming regions to identify asymmetries that are suggestive of mass inflow. Three quantitative measures of line asymmetry are used to indicate whether a line profile is blue, red or neither - the ratio of the temperature of the blue and red peaks, the line skew and the dimensionless parameter delta_v. We find nine HCO+ 3-2 line profiles with a significant blue asymmetry and four with significant red asymmetric profiles. Comparing our HCO+ 3-2 results to HCN 3-2 observations from Wu et al. (2003, 2010), we find that eight of the blue and three of red have profiles with the same asymmetry in HCN. The eight sources with blue asymmetries in both tracers are considered strong candidates for inflow. Quantitative measures of the asymmetry (e.g. delta_v) tend to be larger for HCN. This, combined with possible HCO+ abundance enhancements in outflows, suggests that HCN may be a better tracer of inflow. Understanding the behavior of common molecular tracers like HCO+ in clumps of different masses is important for properly analyzing the line profiles seen in a sample of sources representing a broad range of clump masses. Such studies will soon be possible with the large number of sources with possible self-absorption seen in spectroscopic follow-up observations of clumps identified in the Bolocam Galactic Plane Survey.

The existence of a significant flux of antiprotons confined to Earth's magnetosphere has been considered in several theoretical works. These antiparticles are produced in nuclear interactions of energetic cosmic rays with the terrestrial atmosphere and accumulate in the geomagnetic field at altitudes of several hundred kilometers. A contribution from the decay of albedo antineutrons has been hypothesized in analogy to proton production by neutron decay, which constitutes the main source of trapped protons at energies above some tens of MeV. This Letter reports the discovery of an antiproton radiation belt around the Earth. The trapped antiproton energy spectrum in the South Atlantic Anomaly (SAA) region has been measured by the PAMELA experiment for the kinetic energy range 60--750 MeV. A measurement of the atmospheric sub-cutoff antiproton spectrum outside the radiation belts is also reported. PAMELA data show that the magnetospheric antiproton flux in the SAA exceeds the cosmic-ray antiproton flux by three orders of magnitude at the present solar minimum, and exceeds the sub-cutoff antiproton flux outside radiation belts by four orders of magnitude, constituting the most abundant source of antiprotons near the Earth.

We present the space spectrometer PAMELA observations of proton and helium fluxes during the December 13 and 14, 2006 solar particle events. This is the first direct measurement of the solar energetic particles in space with a single instrument in the energy range from $\sim$ 80 MeV/n up to $\sim$ 3 GeV/n. In the event of December 13 measured energy spectra of solar protons and helium were compared with results obtained by neutron monitors and other detectors. Our measurements show a spectral behaviour different from those derived from the neutron monitor network. No satisfactory analytical fitting was found for the energy spectra. During the first hours of the December 13 event solar energetic particles spectra were close to the exponential form demonstrating rather significant temporal evolution. Solar He with energy up to ~1 GeV/n was recorded on December 13. In the event of December 14 energy of solar protons reached ~600 MeV whereas maximum energy of He was below 100 MeV/n. The spectra were slightly bended in the lower energy range and preserved their form during the second event. Difference in the particle flux appearance and temporal evolution in these two events may argue for a special conditions leading to acceleration of solar particles up to relativistic energies.

We report two new low metallicity blue compact dwarf galaxies (BCDs), WISEP J080103.93+264053.9 (hereafter W0801+26) and WISEP J170233.53+180306.4 (hereafter W1702+18), discovered using the Wide-field Infrared Survey Explorer (WISE). We identified these two BCDs from their extremely red colors at mid-infrared wavelengths, and obtained follow-up optical spectroscopy using the Low Resolution Imaging Spectrometer on Keck I. The mid-infrared properties of these two sources are similar to the well studied, extremely low metallicity galaxy SBS 0335-052E. We determine metallicities of 12 + log(O/H) = 7.75 and 7.63 for W0801+26 and W1702+18, respectively, placing them amongst a very small group of very metal deficient galaxies (Z < 1/10 Zsun). Their > 300 Angstrom Hbeta equivalent widths, similar to SBS 0335-052E, imply the existence of young (< 5 Myr) star forming regions. We measure star formation rates of 2.6 and 10.9 Msun/yr for W0801+26 and W1702+18, respectively. These BCDs, showing recent star formation activity in extremely low metallicity environments, provide new laboratories for studying star formation in extreme conditions and are low-redshift analogs of the first generation of galaxies to form in the universe. Using the all-sky WISE survey, we discuss a new method to identify similar star forming, low metallicity BCDs.

Panchromatic spectral energy distribution (SED) from the ultraviolet (UV), optical to infrared (IR) photometry of NGC 628, combined with the evolutionary stellar population synthesis, is used to derive the spatially resolved age, metallicity and reddening maps. These parameter distributions show that the bulge of this galaxy is a disk-like pseudobulge, which has the S{\'e}rsic index close to the exponential law, rich gas, and a young circumnuclear ring structure. We also discover the disk has two distinct regions with different radial age and metallicity gradients. The inner region is older and has a much steeper age gradient than the outer region of the disk. Both these two regions and the central young structure can be seen in the radial profile of the optical color. Based on the age and reddening distributions, we consider that the pseudobulge and disk are likely to have grown via the secular evolution, which is the redistribution of mass and energy through the angular momentum transport caused by the non-axisymmetric potential of the spirals. However, possible gas accretion events could affect the outer region of the disk, due to abundant H{\sc i} gas accumulating in the outer disk.

We present a detailed analysis of the gamma-ray emission from HESS J1745-303 with the data obtained by the Fermi Gamma-ray Space Telescope in the first ~29 months observation.The source can be clearly detected at the level of ~18-sigma and ~6-sigma in 1-20 GeV and 10-20 GeV respectively. Different from the results obtained by the Compton Gamma-ray Observatory, we do not find any evidence of variability. Most of emission in 10-20 GeV is found to coincide with the region C of HESS J1745-303. A simple power-law is sufficient to describe the GeV spectrum with a photon index of ~2.6. The power-law spectrum inferred in the GeV regime can be connected to that of a particular spatial component of HESS J1745-303 in 1-10 TeV without any spectral break. These properties impose independent constraints for understanding the nature of this "dark particle accelerator".

We present observations of HCO+ and H^13CO+, N2H+, HCS+, HCN and HN^13C, SO and ^34SO, CCH, SO_2, and CH_3OH-E towards a sample of 27 high-mass clumps coincident with water maser emission. All transitions are observed with or convolved to nearly identical resolution (30"), allowing for inter-comparison of the clump properties derived from the mapped transitions. We find N2H+ emission is spatially differentiated compared to the dust and the other molecules towards a few very luminous cores (10 of 27) and the N2H+ integrated intensity does not correlate well with dust continuum flux. We calculate the effective excitation density, n_eff, the density required to excite a 1 K line in T_kin=20 K gas for each molecular tracer. The intensity of molecular tracers with larger effective excitation densities (n_eff > 10^5 cm^-3) appear to correlate more strongly with the submillimeter dust continuum intensity. The median sizes of the clumps are anti-correlated with the n_eff of the tracers (which span more than three orders of magnitude). Virial mass is not correlated with n_eff, especially where the lines are optically thick as the linewidths may be broadened significantly by non-virial motions. The median mass surface density and median volume density of the clumps is correlated with n_eff indicating the importance of understanding the excitation conditions of the molecular tracer when deriving the average properties of an ensemble of cores.

We report the results from Chandra and XMM-Newton observations of a sample of 10 type 1 quasars selected to have unusual UV emission-line properties (weak and blueshifted high-ionization lines; strong UV Fe emission) similar to those of PHL 1811, a confirmed intrinsically X-ray weak quasar. These quasars were identified by the Sloan Digital Sky Survey at high redshift (z~2.2); eight are radio quiet while two are radio intermediate. All of the radio-quiet PHL 1811 analogs are notably X-ray weak by a mean factor of ~13. These sources lack broad absorption lines and have blue UV/optical continua, suggesting they are intrinsically X-ray weak. However, their average X-ray spectrum appears to be harder than those of typical quasars, which may indicate the presence of heavy intrinsic X-ray absorption. Our radio-quiet PHL 1811 analogs support a connection between an X-ray weak spectral energy distribution and PHL 1811-like UV emission lines; this connection provides an economical way to identify X-ray weak type 1 quasars. The fraction of radio-quiet PHL 1811 analogs in the radio-quiet quasar population is estimated to be < 1.2%. We have investigated correlations between relative X-ray brightness and UV emission-line properties for a sample combining radio-quiet PHL 1811 analogs, PHL 1811, and typical type 1 quasars. These correlation analyses suggest that PHL 1811 analogs may have extreme wind-dominated broad emission-line regions. Observationally, radio-quiet PHL 1811 analogs appear to be a subset (~30%) of radio-quiet weak-line quasars. The existence of a subset of quasars in which high-ionization "shielding gas" covers most of the BELR, but little more than the BELR, could potentially unify the PHL 1811 analogs and WLQs. The two radio-intermediate PHL 1811 analogs are X-ray bright. One of them appears to have jet-dominated X-ray emission, while the nature of the other remains unclear.

A strong excess in a form of a wide peak in the energy range of 300-800 GeV was discovered in the first measurements of the electron spectrum in the energy range from 20 GeV to 3 TeV by the balloon-borne experiment ATIC (J. Chang et al. Nature, 2008). The experimental data processing and analysis of the electron spectrum with different criteria for selection of electrons, completely independent of the results reported in (J. Chang et al. Nature, 2008) is employed in the present paper. The new independent analysis generally confirms the results of (J. Chang et al. Nature, 2008), but shows that the spectrum in the region of the excess is represented by a number of narrow peaks. The measured spectrum is compared to the spectrum of (J. Chang et al. Nature, 2008) and to the spectrum of the Fermi/LAT experiment.

Protons and helium nuclei are the most abundant components of the cosmic radiation. Precise measurements of their fluxes are needed to understand the acceleration and subsequent propagation of cosmic rays in the Galaxy. We report precision measurements of the proton and helium spectra in the rigidity range 1 GV-1.2 TV performed by the satellite-borne experiment PAMELA. We find that the spectral shapes of these two species are different and cannot be well described by a single power law. These data challenge the current paradigm of cosmic-ray acceleration in supernova remnants followed by diffusive propagation in the Galaxy. More complex processes of acceleration and propagation of cosmic rays are required to explain the spectral structures observed in our data.

Precision measurements of the electron component in the cosmic radiation provide important information about the origin and propagation of cosmic rays in the Galaxy. Here we present new results regarding negatively charged electrons between 1 and 625 GeV performed by the satellite-borne experiment PAMELA. This is the first time that cosmic-ray electrons have been identified above 50 GeV. The electron spectrum can be described with a single power law energy dependence with spectral index -3.18 +- 0.05 above the energy region influenced by the solar wind (> 30 GeV). No significant spectral features are observed and the data can be interpreted in terms of conventional diffusive propagation models. However, the data are also consistent with models including new cosmic-ray sources that could explain the rise in the positron fraction.

The delensing procedure is an effective tool for removing lensing-induced $B$-mode polarization in the Cosmic Microwave Background to allow for deep searches of primordial $B$-modes. However, the delensing algorithm existing in the literature breaks down if the target $B$-mode signals overlap significantly with the lensing $B$-mode ($\sim 300<l<2000$) in multipole-$l$ space. In this paper, we identify the cause of the breakdown to be correlations between the input $B$-map and the deflection field estimator ($EB$). The amplitude of this bias is quantified by numerical simulations and compared to the analytically derived functional form. We also propose a revised delensing procedure that circumvents the bias. While the newly identified bias does not affect the search of degree scale $B$-mode generated by tensor perturbations, the modified delensing algorithm makes it possible to perform deep searches of high-$l$ $B$-modes such as those generated by patchy reionization, cosmic strings, or rotation of polarization angles. Finally, we estimate how well future polarization experiments can do in detecting tensor- and cosmic string- generated $B$-mode after delensing and comment on different survey strategies.

Cosmic-ray proton and helium spectra have been measured with the balloon-borne Cosmic Ray Energetics And Mass experiment flown for 42 days in Antarctica in the 2004-2005 austral summer season. High-energy cosmic-ray data were collected at an average altitude of ~38.5 km with an average atmospheric overburden of ~3.9 g cm$^{-2}$. Individual elements are clearly separated with a charge resolution of ~0.15 e (in charge units) and ~0.2 e for protons and helium nuclei, respectively. The measured spectra at the top of the atmosphere are represented by power laws with a spectral index of -2.66 $\pm$ 0.02 for protons from 2.5 TeV to 250 TeV and -2.58 $\pm$ 0.02 for helium nuclei from 630 GeV/nucleon to 63 TeV/nucleon. They are harder than previous measurements at a few tens of GeV/nucleon. The helium flux is higher than that expected from the extrapolation of the power law fitted to the lower-energy data. The relative abundance of protons to helium nuclei is 9.1 $\pm$ 0.5 for the range from 2.5 TeV/nucleon to 63 TeV/nucleon. This ratio is considerably smaller than the previous measurements at a few tens of GeV/nucleon.

In this paper we present a multicolor photometry for A2589 ($z=0.0414$) with 15 intermediate bands in the Beijing-Arizona-Taiwan-Connecticut (BATC) system which covers an optical wavelength range from 3000 \AA\ to 10000 \AA. The spectral energy distributions (SEDs) for more than 5000 sources are achieved down to {\it V} $\sim$ 20 mag in about 1 deg$^{2}$ field. A2589 has been also covered by the Sloan Digital Sky Survey (SDSS) in photometric mode only. A cross-identification of the BATC-detected galaxies with the SDSS photometric catalog achieves 1199 galaxies brighter than $i=19.5$ mag, among which 68 member galaxies with known spectroscopic redshifts are found. After combining the SDSS five-band photometric data and the BATC SEDs, the technique of photometric redshift is applied to these galaxies for selecting faint member galaxies. The color-magnitude relation is taken as a further restriction of early-type cluster galaxies. As a result, 106 galaxies are newly selected as member galaxies. Spatial distribution of member galaxies shows a north-south elongation which agrees with the X-ray brightness profile and the orientation of central cD galaxy, NGC 7647. No substructures are detected on the basis of positions and radial velocities of cluster galaxies, indicating that A2589 is a well-relaxed system. The luminosity function of A2589 exhibits a peak at $M_{R} \sim -20$ mag and a dip at $M_{R} \sim -19 $ mag. The low-density outer regions are the preferred habitat of faint galaxies. With the evolutionary population synthesis model, PEGASE, the environmental effect on the star formation properties for 68 spectroscopically confirmed member galaxies is studied. The outlier faint galaxies tend to have longer time scales of star formation, shorter mean stellar ages, and lower metallicities of interstellar medium, which can be interpreted in the context of hierarchical cosmological scenario.

The final results of processing the data from the balloon-born experiment ATIC-2 (Antarctica, 2002-2003) for the energy spectra of protons and He, C, O, Ne, Mg, Si, and Fe nuclei, the spectrum of all particles, and the mean logarithm of atomic weight of primary cosmic rays as a function of energy are presented. The final results are based on improvement of the methods used earlier, in particular, considerably increased resolution of the charge spectrum. The preliminary conclusions on the significant difference in the spectra of protons and helium nuclei (the proton spectrum is steeper) and the non-power character of the spectra of protons and heavier nuclei (flattening of carbon spectrum at energies above 10 TeV) are confirmed. A complex structure of the energy dependence of the mean logarithm of atomic weight is found.

The i-band observing conditions at Dome A on the Antarctic plateau have been investigated using data acquired during 2008 with the Chinese Small Telescope ARray. The sky brightness, variations in atmospheric transparency, cloud cover, and the presence of aurorae are obtained from these images. The median sky brightness of moonless clear nights is 20.5 mag arcsec^{-2} in the SDSS $i$ band at the South Celestial Pole (which includes a contribution of about 0.06 mag from diffuse Galactic light). The median over all Moon phases in the Antarctic winter is about 19.8 mag arcsec^{-2}. There were no thick clouds in 2008. We model contributions of the Sun and the Moon to the sky background to obtain the relationship between the sky brightness and transparency. Aurorae are identified by comparing the observed sky brightness to the sky brightness expected from this model. About 2% of the images are affected by relatively strong aurorae.

VDB0-B195D is a massive, blue star cluster in M31. It was observed as part of the Beijing-Arizona-Taiwan-Connecticut (BATC) Multicolor Sky Survey using 15 intermediate-band filters covering a wavelength range of 3000--10,000 \AA. Based on aperture photometry, we obtain its spectral-energy distribution (SED) as defined by the 15 BATC filters. We apply previously established relations between the BATC intermediate-band and the Johnson-Cousins $UBVRI$ broad-band systems to convert our BATC photometry to the standard system. A detailed comparison shows that our newly derived $VRI$ magnitudes are fully consistent with previous results, while our new $B$ magnitude agrees to within $2\sigma$. In addition, we determine the cluster's age and mass by comparing its SED (from 3000 to 20,000{\AA}, comprising photometric data in the 15 BATC intermediate bands, optical broad-band $BVRI$, and 2MASS near-infrared $JHK_{\rm s}$ data) with theoretical stellar population synthesis models, resulting in age and mass determinations of $60.0\pm 8.0$~Myr and $(1.1-1.6) \times 10^5 M_\odot$, respectively. This age and mass confirms previous suggestions that VDB0-B195D is a young massive cluster in M31.

Interferometric millimeter observations of the cosmic microwave background and clusters of galaxies with arcmin resolutions require antenna arrays with short spacings. Having all antennas co-mounted on a single steerable platform sets limits to the overall weight. A 25 kg lightweight novel carbon-fiber design for a 1.2 m diameter Cassegrain antenna is presented. The finite element analysis predicts excellent structural behavior under gravity, wind and thermal load. The primary and secondary mirror surfaces are aluminum coated with a thin TiO$_2$ top layer for protection. A low beam sidelobe level is achieved with a Gaussian feed illumination pattern with edge taper, designed based on feedhorn antenna simulations and verified in a far field beam pattern measurement. A shielding baffle reduces inter-antenna coupling to below $\sim$ -135 dB. The overall antenna efficiency, including a series of efficiency factors, is estimated to be around 60%, with major losses coming from the feed spillover and secondary blocking. With this new antenna, a detection rate of about 50 clusters per year is anticipated in a 13-element array operation.

This paper aims to study the color behavior of the BL Lac object OJ 287 during optical outburst. According to the revisit of the data from the OJ-94 monitoring project and the analysis the data obtained with the 60/90 cm Schmidt Telescope of NAOC, we found a bluer-when-brighter chromatism in this object. The amplitude of variation tends to decrease with the decrease of frequency. These results are consistent with the shock-in-jet model. We made some simulations and confirmed that both amplitude difference and time delay between variations at different wavelengths can result in the phenomenon of bluer-when-brighter. Our observations confirmed that OJ 287 underwent a double-peaked outburst after about 12 years from 1996, which provides further evidence for the binary black hole model in this object.

We provide a semi-analytic study of the small scale aspects of the power spectra of warm dark matter (WDM) candidates that decoupled while relativistic with arbitrary distribution functions. These are characterized by two widely different scales $k_{eq} \sim 0.01\,(\mathrm{Mpc})^{-1}$ and $k_{fs}= \sqrt{3}\,k_{eq}/2\,< V^2_{eq} >^{1/2} $ with $< V^2_{eq} >^{1/2} \ll 1 $ the velocity dispersion at matter radiation equality. Density perturbations evolve through three stages: radiation domination when the particle is relativistic and non-relativistic and matter domination. An early ISW effect during the first stage leads to an enhancement of density perturbations and a plateau in the transfer function for $k \lesssim k_{fs}$. An effective fluid description emerges at small scales which includes the effects of free streaming in initial conditions and inhomogeneities. The transfer function features \emph{WDM-acoustic oscillations} at scales $k \gtrsim 2 \,k_{fs}$. We study the power spectra for two models of sterile neutrinos with $m \sim \,\mathrm{keV}$ produced non-resonantly, at the QCD and EW scales respectively. The latter case yields acoustic oscillations on mass scales $\sim 10^{8}\,M_{\odot}$. Our results reveal a \emph{quasi-degeneracy} between the mass, distribution function and decoupling temperature suggesting caveats on the constraints on the mass of a sterile neutrino from current WDM N-body simulations and Lyman-$\alpha$ forest data. A simple analytic interpolation of the power spectra between large and small scales and its numerical implementation is given.

We present the results of an investigation into the X-ray properties of radio-intermediate and radio-loud quasars (RIQs and RLQs, respectively). We combine large, modern optical (e.g., SDSS) and radio (e.g., FIRST) surveys with archival X-ray data from Chandra, XMM-Newton, and ROSAT to generate an optically selected sample that includes 188 RIQs and 603 RLQs. This sample is constructed independently of X-ray properties but has a high X-ray detection rate (85%); it provides broad and dense coverage of the l-z plane, including at high redshifts (22% of objects have z=2-5), and it extends to high radio-loudness values (33% of objects have R*=3-5, using logarithmic units). We measure the "excess" X-ray luminosity of RIQs and RLQs relative to radio-quiet quasars (RQQs) as a function of radio loudness and luminosity, and parameterize the X-ray luminosity of RIQs and RLQs both as a function of optical/UV luminosity and also as a joint function of optical/UV and radio luminosity. RIQs are only modestly X-ray bright relative to RQQs; it is only at high values of radio-loudness (R*>3.5) and radio luminosity that RLQs become strongly X-ray bright. We find no evidence for evolution in the X-ray properties of RIQs and RLQs with redshift (implying jet-linked IC/CMB emission does not contribute substantially to the nuclear X-ray continuum). Finally, we consider a model in which the nuclear X-ray emission contains both disk/corona-linked and jet-linked components and demonstrate that the X-ray jet-linked emission is likely beamed but to a lesser degree than applies to the radio jet. This model is used to investigate the increasing dominance of jet-linked X-ray emission at low inclinations.

This is the first paper of our series of high resolution (1") studies of the massive star forming region G10.6--0.4. We present the emission line observations of the hot core type tracers (O$^{13}$CS, OCS, SO$_{2}$) with $\sim$0$"$.5 resolution. By comparing the results to the high--resolution NH$_{3}$ absorption line observation, we confirm for the first time the rotationally flattened hot toroid in the central $<$0.1 pc region, which has a rotational axis perpendicular to its geometrical major axis. In addition, we present the observations of NH$_{3}$, $^{13}$CS, and CH$_{3}$CN with $\sim$1$"$ resolution, and follow the dynamics of the molecular accretion flow from the 0.3 pc radius to the inner 0.03 pc radius. With reference to the rotational axis of the hot toroid, we measure the rotational velocity from the molecular emission in the region. The results are consistent with an envelope with a rapid decrease of the specific angular momentum from the outer to the inner region. These new results improve the current understanding of the molecular accretion flow in an ultracompact (UC) H\textsc{ii} region created by the embedded O-type cluster.

We investigate the contamination of the Sunyaev--Zel'dovich (SZ) effect for six galaxy clusters, A1689, A1995, A2142, A2163, A2261, and A2390, observed by the Y. T. Lee Array for Microwave Background Anisotropy during 2007. With the range of baselines used, we find that the largest effect (of order 13%-50% of the central SZ flux density) comes from primary anisotropies in the cosmic microwave background and exceeds the thermal noise in all six cases. Contamination from discrete radio sources is estimated to be at a level of (3%-60%) of the central SZ flux density. We use the statistics of these contaminating sources to estimate and correct the errors in the measured SZ effects of these clusters.

Clusters of galaxies have been used extensively to determine cosmological parameters. A major difficulty in making best use of Sunyaev-Zel'dovich (SZ) and X-ray observations of clusters for cosmology is that using X-ray observations it is difficult to measure the temperature distribution and therefore determine the density distribution in individual clusters of galaxies out to the virial radius. Observations with the new generation of SZ instruments are a promising alternative approach. We use clusters of galaxies drawn from high-resolution adaptive mesh refinement (AMR) cosmological simulations to study how well we should be able to constrain the large-scale distribution of the intra-cluster gas (ICG) in individual massive relaxed clusters using AMiBA in its configuration with 13 1.2-m diameter dishes (AMiBA13) along with X-ray observations. We show that non-isothermal beta models provide a good description of the ICG in our simulated relaxed clusters. We use simulated X-ray observations to estimate the quality of constraints on the distribution of gas density, and simulated SZ visibilities (AMiBA13 observations) for constraints on the large-scale temperature distribution of the ICG. We find that AMiBA13 visibilities should constrain the scale radius of the temperature distribution to about 50% accuracy. We conclude that the upgraded AMiBA, AMiBA13, should be a powerful instrument to constrain the large-scale distribution of the ICG.

We have compiled a sample of 14 of the optically brightest radio-quiet quasars ($m_{i}$~$\le$~17.5 and $z$~$\ge$~1.9) in the Sloan Digital Sky Survey Data Release 5 quasar catalog that have C IV mini-BALs present in their spectra. X-ray data for 12 of the objects were obtained via a Chandra snapshot survey using ACIS-S, while data for the other two quasars were obtained from archival XMM-Newton observations. Joint X-ray spectral analysis shows the mini-BAL quasars have a similar average power-law photon index ($\Gamma\approx1.9$) and level of intrinsic absorption ($N_H \lesssim 8\times 10^{21} \ {\rm cm}^{-2}$) as non-BMB (neither BAL nor mini-BAL) quasars. Mini-BAL quasars are more similar to non-BMB quasars than to BAL quasars in their distribution of relative X-ray brightness (assessed with $\Delta\alpha_{\rm ox}$). Relative colors indicate mild dust reddening in the optical spectra of mini-BAL quasars. Significant correlations between $\Delta\alpha_{\rm ox}$ and UV absorption properties are confirmed for a sample of 56 sources combining mini-BAL and BAL quasars with high signal-to-noise ratio rest-frame UV spectra, which generally supports models in which X-ray absorption is important in enabling driving of the UV absorption-line wind. We also propose alternative parametrizations of the UV absorption properties of mini-BAL and BAL quasars, which may better describe the broad absorption troughs in some respects.

M31 globular cluster B379 is the first extragalactic cluster, the age of which was determined by main-sequence photometry. In this method, the age of a cluster is obtained by fitting its CMD with stellar evolutionary models. However, different stellar evolutionary models use different parameters of stellar evolution, such as range of stellar masses, different opacities and equations of state, and different recipes, and so on. So, it is interesting to check whether different stellar evolutionary models can give consistent results for the same cluster. Brown et al. (2004a) constrained the age of B379 by comparing its CMD with isochrones of the 2006 VandenBerg models. Using SSP models of BC03 and its multi-photometry, Ma et al. (2007) independently determined the age of B379, which is in good agreement with the determination of Brown et al. (2004a). The BC03 models are calculated based on the Padova evolutionary tracks. It is necessary to check whether the age of B379 which, being determined based on the Padova evolutionary tracks, is in agreement with the determination of Brown et al. (2004a). So, in this paper, we re-determine its age using isochrones of the Padova stellar evolutionary models. In addition, the metal abundance, the distance modulus, and the reddening value for B379 are also determined in this paper. The results obtained in this paper are consistent with the previous determinations, which including the age obtained by Brown et al. (2004a). So, this paper confirms the consistence of the age scale of B379 between the Padova isochrones and the 2006 VandenBerg isochrones, i.e. the results' comparison between Brown et al. (2004a) and Ma et al. (2007) is meaningful. The results obtained in this paper are: the metallicity [M/H]=-0.325, the age $\tau=11.0\pm1.5$ Gyr, the reddening value E(B-V)=0.08, and the distance modulus $(m-M)_{0}=24.44\pm0.10$.

We employ detailed photoionization models to infer the physical conditions of intrinsic narrow absorption line systems found in high resolution spectra of three quasars at z=2.6-3.0. We focus on a family of intrinsic absorbers characterized by N V lines that are strong relative to the Ly-alpha lines. The inferred physical conditions are similar for the three intrinsic N V absorbers, with metallicities greater than 10 times the solar value (assuming a solar abundance pattern), and with high ionization parameters (log U ~ 0). Thus, we conclude that the unusual strength of the N V lines results from a combination of partial coverage, a high ionization state, and high metallicity. We consider whether dilution of the absorption lines by flux from the broad-emission line region can lead us to overestimate the metallicities and we find that this is an unlikely possibility. The high abundances that we infer are not surprising in the context of scenarios in which metal enrichment takes place very early on in massive galaxies. We estimate that the mass outflow rate in the absorbing gas (which is likely to have a filamentary structure) is less than a few solar masses per year under the most optimistic assumptions, although it may be embedded in a much hotter, more massive outflow.

The satellite-borne experiment PAMELA has been used to make a new measurement of the cosmic-ray antiproton flux and the antiproton-to-proton flux ratio which extends previously published measurements down to 60 MeV and up to 180 GeV in kinetic energy. During 850 days of data acquisition approximately 1500 antiprotons were observed. The measurements are consistent with purely secondary production of antiprotons in the galaxy. More precise secondary production models are required for a complete interpretation of the results.

In this paper we investigate the scalar mode of first-order metric perturbations over spatially flat FRW spacetime when the holonomy correction is taken into account in the semi-classical framework of loop quantum cosmology. By means of the Hamiltonian derivation, the cosmological perturbation equations is obtained in longitudinal gauge. It turns out that in the presence of metric perturbation the holonomy effects influence both background and perturbations, and contribute the non-trivial terms $S_{h1}$ and $S_{h2}$ in the cosmological perturbation equations.

We have mapped over 50 massive, dense clumps with four dense gas tracers: HCN J=1-0 and 3-2; and CS J=2-1 and 7-6 transitions. Spectral lines of optically thin H^{13}CN 3-2 and C^{34}S 5-4 were also obtained towards the map centers. These maps usually demonstrate single well-peaked distributions at our resolution, even with higher J transitions. The size, virial mass, surface density, and mean volume density within a well-defined angular size (FWHM) were calculated from the contour maps for each transition. We found that transitions with higher effective density usually trace the more compact, inner part of the clumps but have larger linewidths, leading to an inverse linewidth-size relation using different tracers. The mean surface densities are 0.29, 0.33, 0.78, 1.09 g cm^{-2} within FWHM contours of CS 2-1, HCN 1-0, HCN 3-2 and CS 7-6, respectively. We find no correlation of L_{IR} with surface density and a possible inverse correlation with mean volume density, contrary to some theoretical expectations. We see no evidence in the data for the relation between L'_{mol} and mean density posited by modelers. The correlation between L'_{mol} and the virial mass is roughly linear for each dense gas tracer. A nearly linear correlation was found between the infrared luminosity and the line luminosity of all dense gas tracers for these massive, dense clumps, with a lower cutoff in luminosity at L_{IR}=10^{4.5} Lsun. The L_{IR}-L'_{HCN1-0} correlation agrees well with the one found in galaxies. These correlations indicate a constant star formation rate per unit mass from the scale of dense clumps to that of distant galaxies when the mass is measured for dense gas. These results support the suggestion that starburst galaxies may be understood as having a large fraction of gas in dense clumps.

We identified a large sample of radio quasars, including those with complex radio morphology, from the Sloan Digital Sky Survey (SDSS) and the Faint Images of Radio Sky at Twenty-cm (FIRST). Using this sample, we inspect previous radio quasar samples for selection effects resulting from complex radio morphologies and adopting positional coincidence between radio and optical sources alone. We find that 13.0% and 8.1% radio quasars do not show a radio core within 1.2 and 2 arcsecs of their optical position, and thus are missed in such samples. Radio flux is under-estimated by a factor of more than 2 for an additional 8.7% radio quasars. These missing radio extended quasars are more radio loud with a typical radio-to-optical flux ratio namely radio loudness RL >100, and radio power P >10^{25} W/Hz. They account for more than one third of all quasars with RL>100. The color of radio extended quasars tends to be bluer than the radio compact quasars. This suggests that radio extended quasars are more radio powerful sources, e.g., Fanaroff-Riley type 2 (FR-II) sources, rather than the compact ones viewed at larger inclination angles. By comparison with the radio data from the NRAO VLA Sky Survey (NVSS), we find that for sources with total radio flux less than 3 mJy, low surface brightness components tend to be underestimated by FIRST, indicating that lobes in these faint radio sources are still missed.

The balloon-borne Cosmic Ray Energetics And Mass (CREAM) experiment launched five times from Antarctica has achieved a cumulative flight duration of about 156 days above 99.5% of the atmosphere. The instrument is configured with complementary and redundant particle detectors designed to extend direct measurements of cosmic-ray composition to the highest energies practical with balloon flights. All elements from protons to iron nuclei are separated with excellent charge resolution. Here we report results from the first two flights of ~70 days, which indicate hardening of the elemental spectra above ~200 GeV/nucleon and a spectral difference between the two most abundant species, protons and helium nuclei. These results challenge the view that cosmic-ray spectra are simple power laws below the so-called knee at ~1015 eV. This discrepant hardening may result from a relatively nearby source, or it could represent spectral concavity caused by interactions of cosmic rays with the accelerating shock. Other possible explanations should also be investigated.

During its second Antarctic flight, the CREAM (Cosmic Ray Energetics And Mass) balloon experiment collected data for 28 days, measuring the charge and the energy of cosmic rays (CR) with a redundant system of particle identification and an imaging thin ionization calorimeter. Preliminary direct measurements of the absolute intensities of individual CR nuclei are reported in the elemental range from carbon to iron at very high energy.

We present new measurements of the energy spectra of cosmic-ray (CR) nuclei from the second flight of the balloon-borne experiment CREAM (Cosmic Ray Energetics And Mass). The instrument (CREAM-II) was comprised of detectors based on different techniques (Cherenkov light, specific ionization in scintillators and silicon sensors) to provide a redundant charge identification and a thin ionization calorimeter capable of measuring the energy of cosmic rays up to several hundreds of TeV. The data analysis is described and the individual energy spectra of C, O, Ne, Mg, Si and Fe are reported up to ~ 10^14 eV. The spectral shape looks nearly the same for all the primary elements and can be expressed as a power law in energy E^{-2.66+/-0.04}. The nitrogen absolute intensity in the energy range 100-800 GeV/n is also measured.

The Cosmic Ray Energetics And Mass (CREAM) is a balloon-borne experiment designed to measure the composition and energy spectra of cosmic rays of charge Z = 1 to 26 up to an energy of ~ 10^15 eV. CREAM had two successful flights on long-duration balloons (LDB) launched from Mc- Murdo Station, Antarctica, in December 2004 and December 2005. CREAM-I achieves a substantial measurement redundancy by employing multiple detector systems, namely a Timing Charge Detector and a Silicon Charge Detector (SCD) for particle identification, and a Transition Radiation Detector and a sampling tungsten/scintillating-fiber ionization calorimeter (CAL) for energy measurement. In this paper, preliminary energy spectra of various elements measured with CAL/SCD during the first 42-day flight are presented.

The Cosmic Ray Energetics And Mass (CREAM) balloon experiment had two successful flights in 2004/05 and 2005/06. It was designed to perform energy measurements from a few GeV up to 1000 TeV, taking advantage of different detection techniques. The first instrument, CREAM-1, combined a transition radiation detector with a calorimeter to provide independent energy measurements of cosmicraynuclei. Each detector was calibrated with particle beams in a limited range of energies. In order to assess the absolute energy scale of the instrument and to investigate the systematic effects of each technique, a cross-calibration was performed by comparing the two independent energy estimates on selected samples of oxygen and carbon nuclei.

The Cosmic Ray Energetics And Mass (CREAM) calorimeter is designed to measure the spectra of cosmic-ray particles over the energy range from ~10^11 eV to ~10^15 eV. Its first flight as part of the CREAM-I balloon-borne payload in Antarctica during the 2004/05 season resulted in a recordbreaking 42 days of exposure. Calorimeter calibration using various beam test data will be discussed in an attempt to assess the uncertainties of the energy measurements.

The CREAM calorimeter, designed to measure the spectra of cosmic-ray nuclei from under 1 TeV to 1000 TeV, is a 20 radiation length (X0) deep sampling calorimeter. The calorimeter is comprised of 20 layers of tungsten interleaved with 20 layers of scintillating fiber ribbons, and is preceded by a pair of graphite interaction targets providing about 0.42 proton interaction lengths (\lambda int). The calorimeter was placed in one of CERN's SPS accelerator beams for calibration and testing. Beams of 150 GeV electrons were used for calibration, and a variety of electron, proton, and nuclear fragment beams were used to test the simulation model of the detector. In this paper we discuss the performance of the calorimeter in the electron beam and compare electron beam data with simulation results.

A wideband analog correlator has been constructed for the Yuan-Tseh Lee Array for Microwave Background Anisotropy. Lag correlators using analog multipliers provide large bandwidth and moderate frequency resolution. Broadband IF distribution, backend signal processing and control are described. Operating conditions for optimum sensitivity and linearity are discussed. From observations, a large effective bandwidth of around 10 GHz has been shown to provide sufficient sensitivity for detecting cosmic microwave background variations.

The dynamics of a tachyon field plus a barotropic fluid is investigated in spatially curved FRW universe. We perform a phase-plane analysis and obtain scaling solutions accompanying with a discussion on their stability. Furthermore, we construct the form of scalar potential which may give rise to stable solutions for spatially open and closed universe separately.

The Sunyaev-Zel'dovich Effect (SZE) has been observed toward six massive galaxy clusters, at redshifts 0.091 \leq z \leq 0.322 in the 86-102 GHz band with the Y. T. Lee Array for Microwave Background Anisotropy (AMiBA). We modify an iterative method, based on the isothermal \beta-models, to derive the electron temperature T_e, total mass M_t, gas mass M_g, and integrated Compton Y within r_2500, from the AMiBA SZE data. Non-isothermal universal temperature profile (UTP) \beta models are also considered in this paper. These results are in good agreement with those deduced from other observations. We also investigate the embedded scaling relations, due to the assumptions that have been made in the method we adopted, between these purely SZE-deduced T_e, M_t, M_g and Y. Our results suggest that cluster properties may be measurable with SZE observations alone. However, the assumptions built into the pure-SZE method bias the results of scaling relation estimations and need further study.

The thermal non-Gaussianity in Near-Milne universe is investigated in this letter. Through classifying thermal fluctuations into two types, one characterized by a phase transition and the other without phase transition, we show that for fluctuations undergoing a phase transition, the non-Gaussianity depends linearly on $k$, but in non-phase transition case it is proportional to $k^{-3\gamma},\gamma>0$. Moreover, if the phase transition scenario is similar to that in holographic cosmology, we find that the non-Gaussianity $f_{NL}^{equil}$ can reach ${\cal O}(1)$ or larger than ${\cal O}(1)$ by fine tuning of the equation of state $\omega_{C}$. Especially at the limit of $\omega_{C}\to -5/3$, the non-Gaussianity can be very large. On the other hand, for non-phase transition case, the non-Gaussianity estimator $f_{NL}^{equil}$ is approximately $f_{NL}^{equil}\simeq 5/108$ when the energy is sub-extensive, especially when $\gamma \to 0^{+}$. However, when $\gamma \to 4$, large non-Gaussianity can be obtained.

The Chinese Small Telescope ARray (hereinafter CSTAR) is the first Chinese astronomical instrument on the Antarctic ice cap. The low temperature and low pressure testing of the data acquisition system was carried out in a laboratory refrigerator and on the 4500m Pamirs high plateau, respectively. The results from the final four nights of test observations demonstrated that CSTAR was ready for operation at Dome A, Antarctica. In this paper we present a description of CSTAR and the performance derived from the test observations.

In 2008 January the 24th Chinese expedition team successfully deployed the Chinese Small Telescope ARray (CSTAR) to DomeA, the highest point on the Antarctic plateau. CSTAR consists of four 14.5cm optical telescopes, each with a different filter (g, r, i and open) and has a 4.5degree x 4.5degree field of view (FOV). It operates robotically as part of the Plateau Observatory, PLATO, with each telescope taking an image every 30 seconds throughout the year whenever it is dark. During 2008, CSTAR #1 performed almost flawlessly, acquiring more than 0.3 million i-band images for a total integration time of 1728 hours during 158 days of observations. For each image taken under good sky conditions, more than 10,000 sources down to 16 mag could be detected. We performed aperture photometry on all the sources in the field to create the catalog described herein. Since CSTAR has a fixed pointing centered on the South Celestial Pole (Dec =-90 degree), all the sources within the FOV of CSTAR were monitored continuously for several months. The photometric catalog can be used for studying any variability in these sources, and for the discovery of transient sources such as supernovae, gamma-ray bursts and minor planets.

The PAMELA satellite experiment has measured the cosmic-ray positron fraction between 1.5 GeV and 100 GeV. The need to reliably discriminate between the positron signal and proton background has required the development of an ad hoc analysis procedure. In this paper, a method for positron identification is described and its stability and capability to yield a correct background estimate is shown. The analysis includes new experimental data, the application of three different fitting techniques for the background sample and an estimate of systematic uncertainties due to possible inaccuracies in the background selection. The new experimental results confirm both solar modulation effects on cosmic-rays with low rigidities and an anomalous positron abundance above 10 GeV.

(Abridged) OJ287 is a BL Lac object that has shown double-peaked bursts at regular intervals of ~12 yr during the last ~40 yr. We analyse optical photopolarimetric monitoring data from 2005-2009, during which the latest double-peaked outburst occurred. The aim of this study is twofold: firstly, we aim to analyse variability patterns and statistical properties of the optical polarization light-curve. We find a strong preferred position angle in optical polarization. The preferred position angle can be explained by separating the jet emission into two components: an optical polarization core and chaotic jet emission. The optical polarization core is stable on time scales of years and can be explained as emission from an underlying quiescent jet component. The chaotic jet emission sometimes exhibits a circular movement in the Stokes plane. We interpret these events as a shock front moving forwards and backwards in the jet, swiping through a helical magnetic field. Secondly, we use our data to assess different binary black hole models proposed to explain the regularly appearing double-peaked bursts in OJ287. We compose a list of requirements a model has to fulfil. The list includes not only characteristics of the light-curve but also other properties of OJ287, such as the black hole mass and restrictions on accretion flow properties. We rate all existing models using this list and conclude that none of the models is able to explain all observations. We discuss possible new explanations and propose a new approach to understanding OJ287. We suggest that both the double-peaked bursts and the evolution of the optical polarization position angle could be explained as a sign of resonant accretion of magnetic field lines, a 'magnetic breathing' of the disc.

We study the production of sterile neutrinos in the region $T\sim M_W$ in an extension beyond the standard model with the see-saw mass matrix originating in Yukawa couplings to Higgs-like scalars with masses and vev's of the order of the electroweak scale. Sterile neutrinos are produced by the decay of scalars and standard model vector bosons. We obtain the index of refraction, dispersion relations, mixing angles in the medium and production rates including those for right-handed sterile neutrinos, from the standard model and beyond the standard model self-energies. For $1 \lesssim M_W/T \lesssim 3$ we find narrow MSW resonances with $k \lesssim T$ for both left and right handed neutrinos even in absence of a lepton asymmetry in the (active) neutrino sector, as well as very low energy ($k/T \ll |\xi|$) narrow MSW resonances in the presence of a lepton asymmetry consistent with the bounds from WMAP and BBN. For small vacuum mixing angle, consistent with observational bounds, the absorptive part of the self-energies lead to a strong damping regime very near the resonances resulting in the \emph{exact} degeneracy of the propagating modes with a concomitant breakdown of adiabaticity. We argue that cosmological expansion sweeps through the resonances, \emph{resonant and non-resonant} sterile neutrino production results in a highly \emph{non-thermal} distribution function enhanced at small momentum $k < T$, with potentially important consequences for their free streaming length and transfer function at small scales.

We investigate the scaling relations between the X-ray and the thermal Sunyaev-Zel'dovich Effect (SZE) properties of clusters of galaxies, using data taken during 2007 by the Y.T. Lee Array for Microwave Background Anisotropy (AMiBA) at 94 GHz for the six clusters A1689, A1995, A2142, A2163, A2261, and A2390. The scaling relations relate the integrated Compton-y parameter Y_{2500} to the X-ray derived gas temperature T_{e}, total mass M_{2500}, and bolometric luminosity L_X within r_{2500}. Our results for the power-law index and normalization are both consistent with the self-similar model and other studies in the literature except for the Y_{2500}-L_X relation, for which a physical explanation is given though further investigation may be still needed. Our results not only provide confidence for the AMiBA project but also support our understanding of galaxy clusters.

We present new measurements of the energy spectra of cosmic-ray (CR) nuclei from the second flight of the balloon-borne experiment Cosmic Ray Energetics And Mass (CREAM). The instrument included different particle detectors to provide redundant charge identification and measure the energy of CRs up to several hundred TeV. The measured individual energy spectra of C, O, Ne, Mg, Si, and Fe are presented up to $\sim 10^{14}$ eV. The spectral shape looks nearly the same for these primary elements and it can be fitted to an $E^{-2.66 \pm 0.04}$ power law in energy. Moreover, a new measurement of the absolute intensity of nitrogen in the 100-800 GeV/$n$ energy range with smaller errors than previous observations, clearly indicates a hardening of the spectrum at high energy. The relative abundance of N/O at the top of the atmosphere is measured to be $0.080 \pm 0.025 $(stat.)$ \pm 0.025 $(sys.) at $\sim $800 GeV/$n$, in good agreement with a recent result from the first CREAM flight.

We report observations of 15 high redshift (z = 1-5) galaxies at 350 microns using the Caltech Submillimeter Observatory and SHARC-II array detector. Emission was detected from eight galaxies, for which far-infrared luminosities, star formation rates, total dust masses, and minimum source size estimates are derived. These galaxies have star formation rates and star formation efficiencies comparable to other high redshift molecular emission line galaxies. The results are used to test the idea that star formation in these galaxies occurs in a large number of basic units, the units being similar to star-forming clumps in the Milky Way. The luminosity of these extreme galaxies can be reproduced in a simple model with (0.9-30) *10^6 dense clumps, each with a luminosity of 5 *10^5 Lsun, the mean value for such clumps in the Milky Way. Radiative transfer models of such clumps can provide reasonable matches to the overall SEDs of the galaxies. They indicate that the individual clumps are quite opaque in the far-infrared. Luminosity to mass ratios vary over two orders of magnitude, correlating strongly with the dust temperature derived from simple fits to the SED. The gas masses derived from the dust modeling are in remarkable agreement with those from CO luminosities, suggesting that the assumptions going into both calculations are reasonable.

The Array for Microwave Background Anisotropy (AMiBA) is a radio interferometer for research in cosmology, currently operating 7 0.6m diameter antennas co-mounted on a 6m diameter platform driven by a hexapod mount. AMiBA is currently the largest hexapod telescope. We briefly summarize the hexapod operation with the current pointing error model. We then focus on the upcoming 13-element expansion with its potential difficulties and solutions. Photogrammetry measurements of the platform reveal deformations at a level which can affect the optical pointing and the receiver radio phase. In order to prepare for the 13-element upgrade, two optical telescopes are installed on the platform to correlate optical pointing tests. Being mounted on different locations, the residuals of the two sets of pointing errors show a characteristic phase and amplitude difference as a function of the platform deformation pattern. These results depend on the telescope's azimuth, elevation and polarization position. An analytical model for the deformation is derived in order to separate the local deformation induced error from the real hexapod pointing error. Similarly, we demonstrate that the deformation induced radio phase error can be reliably modeled and calibrated, which allows us to recover the ideal synthesized beam in amplitude and shape of up to 90% or more. The resulting array efficiency and its limits are discussed based on the derived errors.

We present observations at 1.2 mm with MAMBO-II of a sample of z>~2 radio-intermediate obscured quasars, as well as CO observations of two sources with the Plateau de Bure Interferometer. Five out of 21 sources (24%) are detected at a significance of >=3sigma. Stacking all sources leads to a statistical detection of <S_1.2mm>= 0.96+-0.11 mJy and stacking only the non-detections also yields a statistical detection, with <S_1.2mm>= 0.51+-0.13 mJy. This corresponds to a typical far-infrared luminosity L_FIR~4x10^12 Lsol. If the far-infrared luminosity is powered entirely by star-formation, and not by AGN-heated dust, then the characteristic inferred star-formation rate is ~700 Msol yr-1. This far-infrared luminosity implies a dust mass of M_dust~3x10^8 Msol. We estimate that such large dust masses on kpc scales can plausibly cause the obscuration of the quasars. We present dust SEDs for our sample and derive a mean SED for our sample. This mean SED is not well fitted by clumpy torus models, unless additional extinction and far-infrared re-emission due to cool dust are included. There is a hint that the host galaxies of obscured quasars must have higher far-infrared luminosities and cool-dust masses and are therefore often found at an earlier evolutionary phase than those of unobscured quasars. For one source at z=2.767, we detect the CO(3-2) transition, with S_CO Delta nu=630+-50 mJy km s-1, corresponding to L_CO(3-2)= 3.2x10^7 Lsol, or L'_CO(3-2)=2.4x10^10 K km s-1 pc2. For another source at z=4.17, the lack of detection of the CO(4-3) line yields a limit of L'_CO(4-3)<1x10^10 K km s-1 pc2. Molecular gas masses, gas depletion timescales and gas-to-dust ratios are estimated (Abridged).

An optical photometric observation with the Beijing-Arizona-Taiwan-Connecticut (BATC) multicolor system is carried out for A98 (z=0.104), a galaxy cluster with two large enhancements in X-ray surface brightness. The spectral energy distributions (SEDs) covering 15 intermediate bands are obtained for all sources detected down to V ~ 20 mag in a field of $58' \times 58'$. After the star-galaxy separation by the color-color diagrams, a photometric redshift technique is applied to the galaxy sample for further membership determination. The color-magnitude relation is taken as a further restriction of the early-type cluster galaxies. As a result, a list of 198 faint member galaxies is achieved. Based on newly generated sample of member galaxies, the dynamical substructures, A98N, A98S, and A98W, are investigated in detail. A separate galaxy group, A98X, is also found to the south of main concentration of A98, which is gravitationally unbound to A98. For 74 spectroscopically confirmed member galaxies, the environmental effect on the star formation histories is found. The bright galaxies in the core region are found to have shorter time scales of star formation, longer mean stellar ages, and higher metallicities of interstellar medium, which can be interpreted in the context of hierarchical cosmological scenario.

We use UV/optical and X-ray observations of 272 radio-quiet Type 1 AGNs and quasars to investigate the CIV Baldwin Effect (BEff). The UV/optical spectra are drawn from the Hubble Space Telescope, International Ultraviolet Explorer and Sloan Digital Sky Survey archives. The X-ray spectra are from the Chandra and XMM-Newton archives. We apply correlation and partial-correlation analyses to the equivalent widths, continuum monochromatic luminosities, and alpha_ox, which characterizes the relative X-ray to UV brightness. The equivalent width of the CIV 1549 emission line is correlated with both alpha_ox and luminosity. We find that by regressing l_UV with EW(CIV) and alpha_ox, we can obtain tighter correlations than by regressing l_UV with only EW(CIV). Both correlation and regression analyses imply that l_UV is not the only factor controlling the changes of EW(CIV); alpha_ox (or, equivalently, the soft X-ray emission) plays a fundamental role in the formation and variation of CIV. Variability contributes at least 60% of the scatter of the EW(CIV)-l_UV relation and at least 75% of the scatter of the of the EW(CIV)-alpha_ox relation. In our sample, narrow Fe Kalpha 6.4 keV emission lines are detected in 50 objects. Although narrow Fe Kalpha exhibits a BEff similar to that of CIV, its equivalent width has almost no dependence on either alpha_ox or EW(CIV). This suggests that the majority of narrow Fe Kalpha emission is unlikely to be produced in the broad emission-line region. We do find suggestive correlations between the emission-line luminosities of CIV and Fe Kalpha, which could be potentially used to estimate the detectability of the Fe Kalpha line of quasars from rest-frame UV spectroscopic observations.

We use UV/optical and X-ray observations of 272 radio-quiet Type 1 AGNs and quasars to investigate the CIV Baldwin Effect (BEff). The UV/optical spectra are drawn from the Hubble Space Telescope, International Ultraviolet Explorer and Sloan Digital Sky Survey archives. The X-ray spectra are from the Chandra and XMM-Newton archives. We apply correlation and partial-correlation analyses to the equivalent widths, continuum monochromatic luminosities, and alpha_ox, which characterizes the relative X-ray to UV brightness. The equivalent width of the CIV 1549 emission line is correlated with both alpha_ox and luminosity. We find that by regressing l_UV with EW(CIV) and alpha_ox, we can obtain tighter correlations than by regressing l_UV with only EW(CIV). Both correlation and regression analyses imply that l_UV is not the only factor controlling the changes of EW(CIV); alpha_ox (or, equivalently, the soft X-ray emission) plays a fundamental role in the formation and variation of CIV. Variability contributes at least 60% of the scatter of the EW(CIV)-l_UV relation and at least 75% of the scatter of the of the EW(CIV)-alpha_ox relation. In our sample, narrow Fe Kalpha 6.4 keV emission lines are detected in 50 objects. Although narrow Fe Kalpha exhibits a BEff similar to that of CIV, its equivalent width has almost no dependence on either alpha_ox or EW(CIV). This suggests that the majority of narrow Fe Kalpha emission is unlikely to be produced in the broad emission-line region. We do find suggestive correlations between the emission-line luminosities of CIV and Fe Kalpha, which could be potentially used to estimate the detectability of the Fe Kalpha line of quasars from rest-frame UV spectroscopic observations.

We report the results of a time-series CCD photometric survey of variable stars in the field of the open cluster NGC 2126. In about one square degree field covering the cluster, a total number of 21 variable candidates are detected during this survey, of which 16 are newly found. The periods, classifications and spectral types of 14 newly discovered variables are discussed, which consist of six eclipsing binaries systems, three pulsating variable stars, three long period variables, one RS CVn star, and one W UMa or $\delta$ Scuti star. In addition, there are two variable candidates, the properties of which cannot be determined in this paper. By a method based on fitting spectral energy distributions(SEDs) of stars with theoretical ones, the membership probabilities and the fundamental parameters of this cluster are determined. As a result, five variables are probably members of NGC 2126. The fundamental parameters of the this cluster are determined as: the metallicity to be 0.008 $Z_\odot$, the age $\log(t)$= 8.95, the distance modulus $(m-M)_{0}$= 10.34 and the reddening value $E(B-V)$=0.55 mag.

The 15th of June 2006, the PAMELA satellite-borne experiment was launched from the Baikonur cosmodrome and it has been collecting data since July 2006. The apparatus comprises a time-of-flight system, a silicon-microstrip magnetic spectrometer, a silicon-tungsten electromagnetic calorimeter, an anticoincidence system, a shower tail counter scintillator and a neutron detector. The combination of these devices allows precision studies of the charged cosmic radiation to be conducted over a wide energy range (100 MeV -- 100's GeV) with high statistics. The primary scientific goal is the measurement of the antiproton and positron energy spectrum in order to search for exotic sources, such as dark matter particle annihilations. PAMELA is also searching for primordial antinuclei (anti-helium) and testing cosmic-ray propagation models through precise measurements of the anti-particle energy spectrum and precision studies of light nuclei and their isotopes. Moreover, PAMELA is investigating phenomena connected with solar and earth physics.

Titanium is a rare, secondary nucleus among Galactic cosmic rays. Using the Silicon matrix in the ATIC experiment, Titanium has been separated. The energy dependence of the Ti to Fe flux ratio in the energy region from 5 GeV per nucleon to about 500 GeV per nucleon is presented.

This paper supplements Jiang et al. (2003), who studied 172 M31 globular clusters (GCs) and globular cluster candidates from Battistini et al. (1987) on the basis of integrated photometric measurements in the Beijing-Arizona-Taiwan-Connecticut (BATC) photometric system. Here, we present multicolor photometric CCD data (in the BATC system) for the remaining 39 M31 GCs and candidates. In addition, the ages of 35 GCs are constrained by comparing our accurate photometry with updated theoretical stellar synthesis models. We use photometric measurements from GALEX in the far- and near-ultraviolet and 2MASS infrared $JHK_s$ data, in combination with optical photometry. Except for two clusters, the ages of the other sample GCs are all older than 1 Gyr. Their age distribution shows that most sample clusters are younger than 6 Gyr, with a peak at ~3 Gyr, although the `usual' complement of well-known old GCs (i.e., GCs of similar age as the majority of the Galactic GCs) is present as well.

The Y. T. Lee Array for Microwave Background (AMiBA) has reported the first science results on the detection of galaxy clusters via the Sunyaev Zel'dovich effect. The science objectives required small reflectors in order to sample large scale structures (20') while interferometry provided modest resolutions (2'). With these constraints, we designed for the best sensitivity by utilizing the maximum possible continuum bandwidth matched to the atmospheric window at 86-102GHz, with dual polarizations. A novel wide-band analog correlator was designed that is easily expandable for more interferometer elements. MMIC technology was used throughout as much as possible in order to miniaturize the components and to enhance mass production. These designs will find application in other upcoming astronomy projects. AMiBA is now in operations since 2006, and we are in the process to expand the array from 7 to 13 elements.

The Y.T. Lee Array for Microwave Background Anisotropy (AMiBA) started scientific operation in early 2007. This work describes the optimization of the system performance for the measurements of the Sunyaev-Zel'dovich effect for six massive galaxy clusters at redshifts $0.09 - 0.32$. We achieved a point source sensitivity of $63\pm 7$ mJy with the seven 0.6m dishes in 1 hour of on-source integration in 2-patch differencing observations. We measured and compensated for the delays between the antennas of our platform-mounted interferometer. Beam switching was used to cancel instrumental instabilities and ground pick up. Total power and phase stability were good on time scales of hours, and the system was shown to integrate down on equivalent timescales of 300 hours per baseline/correlation, or about 10 hours for the entire array. While the broadband correlator leads to good sensitivity, the small number of lags in the correlator resulted in poorly measured bandpass response. We corrected for this by using external calibrators (Jupiter and Saturn). Using Jupiter as the flux standard, we measured the disk brightness temperature of Saturn to be $149^{+5}_{-12}$ K.

AMiBA is the largest hexapod astronomical telescope in current operation. We present a description of this novel hexapod mount with its main mechanical components -- the support cone, universal joints, jack screws, and platform -- and outline the control system with the pointing model and the operating modes that are supported. The AMiBA hexapod mount performance is verified based on optical pointing tests and platform photogrammetry measurements. The photogrammetry results show that the deformations in the inner part of the platform are less than 120 micron rms. This is negligible for optical pointing corrections, radio alignment and radio phase errors for the currently operational 7-element compact configuration. The optical pointing error in azimuth and elevation is successively reduced by a series of corrections to about 0.4 arcmin rms which meets our goal for the 7-element target specifications.

We present Submillimeter Array observations toward the 10^{4.7} Lsun star-forming region G240.31+0.07, in the J=2-1 transition of 12CO and 13CO and at 1.3 mm continuum, as well as the 12CO and 13CO observations from the Caltech Submillimeter Observatory to recover the extended emission filtered out by the interferometer. Maps of the 12CO and 13CO emission show a bipolar, wide-angle, quasi-parabolic molecular outflow, roughly coincident with an IR nebula revealed by the Spitzer 3.6 and 4.5 micron emission. The outflow has ~98 Msun molecular gas, making it one of the most massive molecular outflows known, and resulting in a very high mass-loss rate of 4.1 by 10^{-3} Msun yr^{-1} over a dynamical timescale of 2.4 by 10^4 yr. The 1.3 mm continuum observations with a 4" by 3" beam reveal a flattened dusty envelope of ~150 Msun, which is further resolved with a 1.2" by 1" beam into three dense cores with a total mass of ~40 Msun. The central mm core, showing evidence of active star formation, approximately coincides with the geometric center of the bipolar outflow thus most likely harbors the powering source of the outflow. Overall our observations provide the best case to date of a well-defined wide-angle molecular outflow in a >10^4 Lsun star-forming region. The outflow is morphologically and kinematically similar to low-mass protostellar outflows but has two to three orders of magnitude greater mass, momentum, and energy, and is apparently driven by an underlying wide-angle wind, hence further supports that high-mass stars up to late-O types, even in a crowded clustering environment, can form as a scaled-up version of low-mass star formation.

We present observations, analysis and results for the first-year operation of AMiBA, an interferometric experiment designed to study cosmology via the measurement of Cosmic Microwave Background (CMB). AMiBA is the first CMB interferometer operating at 3 mm to have reported successful results, currently with seven close-packed antennas of 60-cm diameter giving a synthesized resolution of around 6'. During 2007 AMiBA detected the Sunyaev-Zel'dovich effects (SZE) of six galaxy clusters at redshift 0.091 <= z <= 0.322. An observing strategy with on-off-source switching is used to minimize the effects from electronic offset and ground pickup. Planets were used to test the observational capability of AMiBA and to calibrate the conversion from correlator time-lag data to visibilities. The detailed formalism for data analysis is given. We summarize our early tests including observations of planets and quasars, and present images, visibility profiles, the estimated central coordinates, sizes, and SZE amplitudes of the galaxy clusters. Science results are summarized. We also discuss possible systematic effects in the results.

We present results from and the analysis of data from MAXIPOL, a balloon-borne experiment designed to measure the polarization in the Cosmic Microwave Background (CMB). MAXIPOL is the first CMB experiment to obtain results using a rotating half-wave plate as a rapid polarization modulator. We report results from observations of a sky area of 8 deg^2 with 10-arcmin resolution, providing information up to l~700. We use a maximum-likelihood method to estimate maps of the Q and U Stokes parameters from the demodulated time streams, and then both Bayesian and frequentist approaches to compute the EE, EB, and BB power spectra. Detailed formalisms of the analyses are given. A variety of tests show no evidence for systematic errors. The Bayesian analysis gives weak evidence for an EE signal. The EE power is 55^{+51}_{-45} \mu K^2 at the 68% confidence level for l=151-693. Its likelihood function is asymmetric and skewed positive such that with a uniform prior the probability of a positive EE power is 96%. The powers of EB and BB signals at the 68% confidence level are 18^{+27}_{-34} \mu K^2 and -31^{+31}_{-19} \mu K^2 respectively and thus consistent with zero. The upper limit of the BB-mode at the 95% confidence level is 9.5 \mu K. Results from the frequentist approach are in agreement within statistical errors. These results are consistent with the current concordance LCDM model.

The source J 1128+5925 was found recently to show strong intraday variability at radio wavelengths and its radio variability may come from interstellar scintillation. In optical, the object was quiet in our 2007 monitoring session. Here we report the results of our new optical monitoring of this source in 2008. In addition to confirm our 2007 results, that the object did not display any clear variation on timescales from hour--day to month, we provide evidence that the object does not vary on timescale of one year, and it is probably intrinsically quiet in optical domain. Its very different behaviors in optical and radio regimes can be naturally explained if its strong radio variability comes from interstellar scintillation.

The Yuan-Tseh Lee Array for Microwave Background Anisotropy (AMiBA) is the first interferometer dedicated to studying the cosmic microwave background (CMB) radiation at 3mm wavelength. The choice of 3mm was made to minimize the contributions from foreground synchrotron radiation and Galactic dust emission. The initial configuration of seven 0.6m telescopes mounted on a 6-m hexapod platform was dedicated in October 2006 on Mauna Loa, Hawaii. Scientific operations began with the detection of a number of clusters of galaxies via the thermal Sunyaev-Zel'dovich effect. We compare our data with Subaru weak lensing data in order to study the structure of dark matter. We also compare our data with X-ray data in order to derive the Hubble constant.

We present a multiwavelength analysis of a sample of four hot (T_X>8keV) X-ray galaxy clusters (A1689, A2261, A2142, and A2390) using joint AMiBA Sunyaev-Zel'dovich effect (SZE) and Subaru weak lensing observations, combined with published X-ray temperatures, to examine the distribution of mass and the intracluster medium (ICM) in massive cluster environments. Our observations show that A2261 is very similar to A1689 in terms of lensing properties. Many tangential arcs are visible around A2261, with an effective Einstein radius \sim 40 arcsec (at z \sim 1.5), which when combined with our weak lensing measurements implies a mass profile well fitted by an NFW model with a high concentration c_{vir} \sim 10, similar to A1689 and to other massive clusters. The cluster A2142 shows complex mass substructure, and displays a shallower profile (c_{vir} \sim 5), consistent with detailed X-ray observations which imply recent interaction. The AMiBA map of A2142 exhibits an SZE feature associated with mass substructure lying ahead of the sharp north-west edge of the X-ray core suggesting a pressure increase in the ICM. For A2390 we obtain highly elliptical mass and ICM distributions at all radii, consistent with other X-ray and strong lensing work. Our cluster gas fraction measurements, free from the hydrostatic equilibrium assumption, are overall in good agreement with published X-ray and SZE observations, with the sample-averaged gas fraction of <f_{gas}> = 0.133 \pm 0.027, for our sample <M_{vir}> = (1.2 \pm 0.1) \times 10^{15} M_{sun} h^{-1}. When compared to the cosmic baryon fraction f_b = \Omega_b/\Omega_m constrained by the WMAP 5-year data, this indicates <f_{gas}>/f_b = 0.78 \pm 0.16, i.e., (22 \pm 16)% of the baryons are missing from the hot phase of clusters.

ASTROD I is a planned interplanetary space mission with multiple goals. The primary aims are: to test General Relativity with an improvement in sensitivity of over 3 orders of magnitude, improving our understanding of gravity and aiding the development of a new quantum gravity theory; to measure key solar system parameters with increased accuracy, advancing solar physics and our knowledge of the solar system and to measure the time rate of change of the gravitational constant with an order of magnitude improvement and the anomalous Pioneer acceleration, thereby probing dark matter and dark energy gravitationally. It is an international project, with major contributions from Europe and China and is envisaged as the first in a series of ASTROD missions. ASTROD I will consist of one spacecraft carrying a telescope, four lasers, two event timers and a clock. Two-way, two-wavelength laser pulse ranging will be used between the spacecraft in a solar orbit and deep space laser stations on Earth, to achieve the ASTROD I goals. A second mission, ASTROD II is envisaged as a three-spacecraft mission which would test General Relativity to one part per billion, enable detection of solar g-modes, measure the solar Lense-Thirring effect to 10 parts per million, and probe gravitational waves at frequencies below the LISA bandwidth. In the third phase (ASTROD III or Super-ASTROD), larger orbits could be implemented to map the outer solar system and to probe primordial gravitational-waves at frequencies below the ASTROD II bandwidth.

This paper reports the discovery of a Very Low Luminosity Object (VeLLO) in the "starless" dense core L328, using the Spitzer Space Telescope and ground based observations from near-infrared to millimeter wavelengths. The Spitzer 8 micron image indicates that L328 consists of three subcores of which the smallest one may harbor a source, L328-IRS while two other subcores remain starless. L328-IRS is a Class 0 protostar according to its bolometric temperature (44 K) and the high fraction ~72 % of its luminosity emitted at sub-millimeter wavelengths. Its inferred "internal luminosity" (0.04 - 0.06 Lsun) using a radiative transfer model under the most plausible assumption of its distance as 200 pc is much fainter than for a typical protostar, and even fainter than other VeLLOs studied previously. Note, however, that its inferred luminosity may be uncertain by a factor of 2-3 if we consider two extreme values of the distance of L328-IRS (125 or 310 pc). Low angular resolution observations of CO do not show any clear evidence of a molecular outflow activity. But broad line widths toward L328, and Spitzer and near-infrared images showing nebulosity possibly tracing an outflow cavity, strongly suggest the existence of outflow activity. Provided that an envelope of at most ~0.1 Msunis the only mass accretion reservoir for L328-IRS, and the star formation efficiency is close to the canonical value ~30%, L328-IRS has not yet accreted more than 0.05 Msun. At the assumed distance of 200 pc, L328-IRS is destined to be a brown dwarf.