We present new observations of the active star-formation region NGC 1333 in the Perseus molecular cloud complex from the James Clerk Maxwell Telescope (JCMT) B-Fields In Star-forming Region Observations (BISTRO) survey with the POL-2 instrument. The BISTRO data cover the entire NGC 1333 complex (~ 1.5 pc x 2 pc) at 0.02 pc resolution and spatially resolve the polarized emission from individual filamentary structures for the first time. The inferred magnetic field structure is complex as a whole, with each individual filament aligned at different position angles relative to the local field orientation. We combine the BISTRO data with low- and high- resolution data derived from Planck and interferometers to study the multi-scale magnetic field structure in this region. The magnetic field morphology drastically changes below a scale of ~ 1 pc and remains continuous from the scales of filaments (~ 0.1 pc) to that of protostellar envelopes (~ 0.005 pc or ~ 1000 au). Finally, we construct simple models in which we assume that the magnetic field is always perpendicular to the long axis of the filaments. We demonstrate that the observed variation of the relative orientation between the filament axes and the magnetic field angles are well reproduced by this model, taking into account the projection effects of the magnetic field and filaments relative to the plane of the sky. These projection effects may explain the apparent complexity of the magnetic field structure observed at the resolution of BISTRO data towards the filament network.

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.

Uncertainty quantification has emerged as a rapidly growing field in nuclear science. Theoretical predictions of physical observables often involve extrapolations to regions that are poorly constrained by laboratory experiments and astrophysical observations. Without properly quantified theoretical errors, such model predictions are of limited value. Also, one often deals with theoretical constructs that involve fundamental quantities that are not accessible to experiment or observation. Particularly relevant in this context is the pressure of pure neutron matter. In this contribution we develop an analytic framework to answer the question of "How can new data reduce uncertainties of current theoretical models?" [P.-G. Reinhard and W. Nazarewicz, Phys. Rev. C81, 051303(R) (2010)]. Simple and insightful expressions are obtained to quantify the impact of one or two new observables on theoretical uncertainties in two critical quantities: the slope of the symmetry energy at saturation density and the pressure of pure neutron matter at twice nuclear matter saturation density.

Recently, repeating fast radio bursts (FRBs) with a period of $P_{\rm FRB}=16.35\pm0.18$ days from FRB 180916.J0158+65 had been reported. It still remains controversial how to give rise to such a periodicity of this FRB. In this Letter, based on an assumption of a young pulsar surrounding by a debris disk, we attempt to diagnose whether the Lense-Thirring precession of the disk on the emitter can produce the observed periodicity. Our calculations indicate that the Lense-Thirring effect of a tilted disk can result in a precession period of 16 days for a mass inflow rate of $0.5-1.5\times10^{18}~\rm g\,s^{-1}$, a spin period of 1-20 ms of the pulsar, and an extremely low viscous parameter $\alpha=10^{-8}$ in the disk. The disk mass and the magnetic field of the pulsar are also constrained to be $\sim10^{-3}~\rm M_{\odot}$ and $< 2.5\times 10^{13}~\rm G$. In our model, a new born pulsar with normal magnetic field and millisecond period would successively experience accretion phase, propeller phase, and is visible as a strong radio source in the current stage. The rotational energy of such a young NS can provide the observed radio bursting luminosity for $400$ years.

The scientific aim of the space gravitational wave (GW) detector Laser Interferometer Space Antenna (LISA) that was scheduled to launch in the early 2030s is to detect the low-frequency GW signals in the Galaxy. Its main candidate GW sources are compact binaries of white dwarfs and neutron stars. In this work, we examine whether compact intermediate-mass black hole (IMBH) X-ray binaries could be potential LISA sources. Our simulations indicate that IMBH binary with a 1000 $M_{\odot}$ IMBH and a 3 $M_{\odot}$ donor star in an orbit of initial orbital period near the so-called bifurcation period of 0.77 d could evolve into an ultra-compact X-ray binary, which will emit GW signals with maximum frequency of 2.5 mHz. According to the evolutionary tracks of characteristic strain, IMBH X-ray binaries with the initial donor-star masses of $1-3~M_{\odot}$ and the initial orbital periods slightly less than the bifurcation periods will be detectable by the LISA in a distance of 15 kpc. Assuming each of 60 Galactic globular clusters hosts a 1000 $M_{\odot}$ IMBH, the maximum number of compact IMBH X-ray binaries that LISA will detect in the Galaxy should be less than ten. Therefore, the detectability of compact IMBH X-ray binaries by the LISA is not optimistic.

Given its relative proximity (~430 pc), compact size (< 20'), young age (~ 0.5 Myr) and rich number of young stellar objects, the Serpens-South star forming region is a promising site for studying young sub-stellar objects, yet the low-mass members of this region remain largely undiscovered. In this paper we report on a deep photometric survey using a custom 1.45 um filter (W-band), as well as standard J and H near-IR filters, in order to identify candidate low-mass young brown-dwarfs in the Serpens-South region. We constructed a reddening-insensitive index (Q) by combining J, H and W-band photometry for survey objects, in order to identify candidate low-mass members of Serpens based on the strength of the water absorption feature at 1.45 um in the atmospheres of mid-M and later objects. We then conducted spectroscopic follow up to confirm youth and spectral type for our candidates. This is the first survey to identify the very low-mass and coolest members of Serpens-South. We identify 4 low-mass candidate Serpens members, which all display IR excess emission, indicating the likely presence of circumstellar disks around them. One of the four candidate low-mass members in our list, SERP182918-020245, exhibits Pa_beta and Br_gamma emission features, confirming its youth and ongoing magnetospheric accretion. Our new candidate members have spectral types >M4 and are the coolest and lowest mass candidate members yet identified in Serpens-South.

Small planets are common around late-M dwarfs and can be detected through highly precise photometry by the transit method. Planets orbiting nearby stars are particularly important as they are often the best-suited for future follow-up studies. We present observations of three nearby M-dwarfs referred to as EIC-1, EIC-2, and EIC-3, and use them to search for transits and set limits on the presence of planets. On most nights our observations are sensitive to Earth-sized transiting planets, and photometric precision is similar to or better than TESS for faint late-M dwarfs of the same magnitude (I=15 mag). We present our photometry and transit search pipeline, which utilizes simple median detrending in combination with transit least squares based transit detection (Hippke & Heller 2019).For these targets, and transiting planets between one and two Earth radii, we achieve an average transit detection probability of 60% between periods of 0.5 and 2 days, 30% between 2 and 5 days,and 10% between 5 and 10 days. These sensitivities are conservative compared to visual searches.

Two components of jets associated with the afterglow of the gamma-ray burst GRB 160623A were observed with multi-frequency observations including long-term monitoring in a sub-millimetre range (230 GHz) using the SMA. The observed light curves with temporal breaks suggests on the basis of the standard forward-shock synchrotron radiation model that the X-ray radiation is narrowly collimated with an opening angle $\theta_{n,j}<\sim6^{\circ}$ whereas the radio radiation originated from wider jets ($\sim27^{\circ}$). The temporal and spectral evolutions of the radio afterglow agree with those expected from a synchrotron radiation modelling with typical physical parameters except for the fact that the observed wide jet opening angle for the radio emission is significantly larger than the theoretical maximum opening angle. By contrast, the opening angle of the X-ray afterglow is consistent with the typical value of GRB jets. Since the theory of the relativistic cocoon afterglow emission is similar to that of a regular afterglow with an opening angle of $\sim30^{\circ}$, the observed radio emission can be interpreted as the shocked jet cocoon emission. This result therefore indicates that the two components of the jets observed in the GRB 160623A afterglow is caused by the jet and the shocked jet cocoon afterglows.

We report on a multi-band variability and correlation study of the TeV blazar Mrk 421 during an exceptional flaring activity observed from 2013 April 11 to 2013 April 19. The study uses, among others, data from GASP-WEBT, Swift, NuSTAR, Fermi-LAT, VERITAS, and MAGIC. The large blazar activity, and the 43 hours of simultaneous NuSTAR and MAGIC/VERITAS observations, permitted variability studies on 15 minute time bins, and over three X-ray bands (3-7 keV, 7-30 keV and 30-80 keV) and three very-high-energy (>0.1 TeV, hereafter VHE) gamma-ray bands (0.2-0.4 TeV, 0.4-0.8 TeV and >0.8 TeV). We detected substantial flux variations on multi-hour and sub-hour timescales in all the X-ray and VHE gamma-ray bands. The characteristics of the sub-hour flux variations are essentially energy-independent, while the multi-hour flux variations can have a strong dependence on the energy of the X-ray and the VHE gamma rays. The three VHE bands and the three X-ray bands are positively correlated with no time-lag, but the strength and the characteristics of the correlation changes substantially over time and across energy bands. Our findings favour multi-zone scenarios for explaining the achromatic/chromatic variability of the fast/slow components of the light curves, as well as the changes in the flux-flux correlation on day-long timescales. We interpret these results within a magnetic reconnection scenario, where the multi-hour flux variations are dominated by the combined emission from various plasmoids of different sizes and velocities, while the sub-hour flux variations are dominated by the emission from a single small plasmoid moving across the magnetic reconnection layer.

We report the results of decade-long (2008-2018) $\gamma$-ray to 1 GHz radio monitoring of the blazar 3C 279, including GASP/WEBT, $\it{Fermi}$ and $\it{Swift}$ data, as well as polarimetric and spectroscopic data. The X-ray and $\gamma$-ray light curves correlate well, with no delay > 3 hours, implying general co-spatiality of the emission regions. The $\gamma$-ray-optical flux-flux relation changes with activity state, ranging from a linear to a more complex dependence. The behaviour of the Stokes parameters at optical and radio wavelengths, including 43 GHz VLBA images, supports either a predominantly helical magnetic field or motion of the radiating plasma along a spiral path. Apparent speeds of emission knots range from 10 to 37c, with the highest values requiring bulk Lorentz factors close to those needed to explain $\gamma$-ray variability on very short time scales. The Mg II emission line flux in the `blue' and `red' wings correlates with the optical synchrotron continuum flux density, possibly providing a variable source of seed photons for inverse Compton scattering. In the radio bands we find progressive delays of the most prominent light curve maxima with decreasing frequency, as expected from the frequency dependence of the $\tau=1$ surface of synchrotron self-absorption. The global maximum in the 86 GHz light curve becomes less prominent at lower frequencies, while a local maximum, appearing in 2014, strengthens toward decreasing frequencies, becoming pronounced at $\sim5$ GHz. These tendencies suggest different Doppler boosting of stratified radio-emitting zones in the jet.

We present the stellar population, using {\it Gaia}\,DR2 parallax, kinematics, and photometry, of the young ($\sim 100$~Myr), nearby ($\sim 230$~pc) open cluster, Blanco1. A total of 644 member candidates are identified via the unsupervised machine learning method \textsc{StarGO} to find the clustering in the 5-dimensional position and proper motion parameter ($X$, $Y$, $Z$, $\mu_\alpha \cos\delta$, $\mu_\delta$) space. Within the tidal radius of $10.0 \pm 0.3$~pc, there are 488 member candidates, 3 times more than those outside. A leading tail and a trailing tail, each of 50--60~pc in the Galactic plane, are found for the first time for this cluster, with stars further from the cluster center streaming away faster, manifest stellar stripping. Blanco1 has a total detected mass of $285\pm32$~M$_\odot $ with a mass function consistent with a slope of $\alpha=1.35\pm0.2$ in the sense of $dN/dm \propto m^{-\alpha}$, in the mass range of 0.25--2.51~M$_\odot $, where $N$ is the number of members and $m$ is stellar mass. A Minimum Spanning Tree ($\Lambda_{\rm MSR}$) analysis shows the cluster to be moderately mass segregated among the most massive members ($\gtrsim 1.4$~M$_\odot$), suggesting an early stage of dynamical disintegration.

This white paper gives a brief summary of the time domain science that has been performed with the JCMT in recent years and highlights the opportunities for continuing work in this field over the next decade. The main focus of this document is the JCMT Transient Survey, a large program initiated in 2015 to measure the frequency and amplitude of variability events associated with protostars in nearby star-forming regions. After summarising the major accomplishments so far, an outline is given for extensions to the current survey, featuring a discussion on what will be possible with the new 850 micron camera that is expected to be installed in late 2022. We also discuss possible applications of submillimetre monitoring to active galactic nuclei, X-ray binaries, asymptotic giant branch stars, and flare stars.

The IC5146 cloud is a nearby star-forming region in Cygnus, consisting of molecular gas filaments in a variety of evolutionary stages. We used optical and near-infrared polarization data toward the IC5146 cloud, reported in the first paper of this series, to reveal the magnetic fields in this cloud. Using the newly released $Gaia$ data, we found that the IC5146 cloud may contain two separate clouds: a first cloud, including the densest main filament at a distance of $\sim$600 pc, and a second cloud, associated with the Cocoon Nebula at a distance of $\sim$800 pc. The spatially averaged H-band polarization map revealed a well-ordered magnetic field morphology, with the polarization segments perpendicular to the main filament but parallel to the nearby sub-filaments, consistent with models assuming that the magnetic field is regulating cloud evolution. We estimated the magnetic field strength using the Davis-Chandrasekhar-Fermi method, and found that the magnetic field strength scales with volume density with a power-law index of $\sim0.5$ in the density range from $N_{H_2}\sim$ 10 to 3000 cm$^{-3}$, which indicates an anisotropic cloud contraction with a preferred direction along the magnetic field. In addition, the mass-to-flux ratio of the cloud gradually changes from subcritical to supercritical from the cloud envelope to the deep regions. These features are consistent with strong magnetic field star-formation models and suggest that the magnetic field is important in regulating the evolution of the IC5146 cloud.

According to the recycling model, neutron stars in low-mass X-ray binaries were spun up to millisecond pulsars (MSPs), which indicates that all MSPs in the Galactic plane ought to be harbored in binaries. However, about $20\%$ Galactic field MSPs are found to be solitary. To interpret this problem, we assume that the accreting neutron star in binaries may collapse and become a strange star when it reaches some critical mass limit. Mass loss and a weak kick induced by asymmetric collapse during the phase transition (PT) from neutron star to strange star can result in isolated MSPs. In this work, we use a population-synthesis code to examine the PT model. The simulated results show that a kick velocity of $\sim60~{\rm km~s}^{-1}$ can produce $\sim6\times10^3$ isolated MSPs and birth rate of $\sim6.6\times10^{-7} {\rm ~yr}^{-1}$ in the Galaxy, which is approximately in agreement with predictions from observations. For the purpose of comparisons with future observation, we also give the mass distributions of radio and X-ray binary MSPs, along with the delay time distribution.

The All-sky Medium Energy Gamma-ray Observatory (AMEGO) is a probe class mission concept that will provide essential contributions to multimessenger astrophysics in the late 2020s and beyond. AMEGO combines high sensitivity in the 200 keV to 10 GeV energy range with a wide field of view, good spectral resolution, and polarization sensitivity. Therefore, AMEGO is key in the study of multimessenger astrophysical objects that have unique signatures in the gamma-ray regime, such as neutron star mergers, supernovae, and flaring active galactic nuclei. The order-of-magnitude improvement compared to previous MeV missions also enables discoveries of a wide range of phenomena whose energy output peaks in the relatively unexplored medium-energy gamma-ray band.

We present a multiwavelength study of the flat-spectrum radio quasar CTA 102 during 2013-2017. We use radio-to-optical data obtained by the Whole Earth Blazar Telescope, 15 GHz data from the Owens Valley Radio Observatory, 91 and 103 GHz data from the Atacama Large Millimeter Array, near-infrared data from the Rapid Eye Monitor telescope, as well as data from the Swift (optical-UV and X-rays) and Fermi ($\gamma$ rays) satellites to study flux and spectral variability and the correlation between flux changes at different wavelengths. Unprecedented $\gamma$-ray flaring activity was observed during 2016 November-2017 February, with four major outbursts. A peak flux of (2158 $\pm$ 63)$\times$10$^{-8}$ ph cm$^{-2}$ s$^{-1}$, corresponding to a luminosity of (2.2 $\pm$ 0.1)$\times$10$^{50}$ erg s$^{-1}$, was reached on 2016 December 28. These four $\gamma$-ray outbursts have corresponding events in the near-infrared, optical, and UV bands, with the peaks observed at the same time. A general agreement between X-ray and $\gamma$-ray activity is found. The $\gamma$-ray flux variations show a general, strong correlation with the optical ones with no time lag between the two bands and a comparable variability amplitude. This $\gamma$-ray/optical relationship is in agreement with the geometrical model that has successfully explained the low-energy flux and spectral behaviour, suggesting that the long-term flux variations are mainly due to changes in the Doppler factor produced by variations of the viewing angle of the emitting regions. The difference in behaviour between radio and higher energy emission would be ascribed to different viewing angles of the jet regions producing their emission.

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

The intra-night trajectories of asteroids can be approximated by straight lines, and so are their intra-night detections. Therefore, the Hough transform, a line detecting algorithm, can be used to connect the line-up detections to find asteroids. We applied this algorithm to a high-cadence Pan-STARRS 1 (PS1) observation, which was originally designed to collect asteroid light curves for rotation period measurements (Chang et al., 2019). The algorithm recovered most of the known asteroids in the observing fields and, moreover, discovered 3574 new asteroids with magnitude mainly of 21.5 < w_p1 < 22.5 mag. This magnitude range is equivalent to sub-kilometer main-belt asteroids (MBAs), which usually lack of rotation period measurements due to their faintness. Using the light curves of the 3574 new asteroids, we obtained 122 reliable rotation periods, of which 13 are super-fast rotators (SRFs; i.e., rotation period of < 2 hr). The required cohesion to survive these SFRs range from tens to thousands of Pa, a value consistent with the known SFRs and the regolith on the Moon and Mars. The higher chance of discovering SFRs here suggests that sub-kilometer MBAs probably harbor more SFRs.

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.

Recent timing observation provides an intriguing result for the braking index of the X-ray pulsar PSR J1640-4631, which has a measured braking index $n=3.15\pm0.03$. The decrease of the inclination angle between between the spin axis and the magnetic axis can be responsible for such a high braking index. However, the physical mechanisms causing the change of the magnetic inclination angle have not been fully understood. In this Letter, we apply a two-dipole model given by Hamil et al (2016) to explain the decrease of the magnetic inclination angle of PSR J1640-4631. The rotation effect of a charged sphere and the magnetization of ferromagnetically ordered material produce magnetic moments $M_{1}$ and $M_{2}$, respectively. There exist a minimum of the potential energy for the magnetic moment $M_{2}$ in the magnetic field of $M_{1}$, hence the $M_{2}$ will freely rotate around the minimum energy position (i. e equilibrium position), similar to a simple pendulum. Our calculation indicate the magnetic moment $M_{2}$ would evolve towards alignment with the spin axis for PSR J1640-4631, and cause the magnetic inclination angle to decrease. The single peak in the pulse profile favors a relatively low change rate of the magnetic inclination angle.

Geodetic Very Long Baseline Interferometry (VLBI) data are capable of measuring the light deflection caused by the gravitational field of the Sun and large planets with high accuracy. The parameter $\gamma$ of the parametrized Post-Newtonian (PPN) formalism estimated using observations of reference radio sources near the Sun should be equal to unity in the general relativity. We have run several VLBI experiments tracking reference radio sources from 1 to 3 degrees from the Sun. The best formal accuracy of the parameter $\gamma$ achieved in the single-session mode is less than 0.01 percent, or better than the formal accuracy obtained with a global solution included all available observations at arbitrary elongation from the Sun. We are planning more experiments starting from 2020 using better observing conditions near the minimum of the Solar activity cycle.

Matter evolved under influence of gravity from minuscule density fluctuations. Non-perturbative structure formed hierarchically over all scales, and developed non-Gaussian features in the Universe, known as the Cosmic Web. To fully understand the structure formation of the Universe is one of the holy grails of modern astrophysics. Astrophysicists survey large volumes of the Universe and employ a large ensemble of computer simulations to compare with the observed data in order to extract the full information of our own Universe. However, to evolve trillions of galaxies over billions of years even with the simplest physics is a daunting task. We build a deep neural network, the Deep Density Displacement Model (hereafter D$^3$M), to predict the non-linear structure formation of the Universe from simple linear perturbation theory. Our extensive analysis, demonstrates that D$^3$M outperforms the second order perturbation theory (hereafter 2LPT), the commonly used fast approximate simulation method, in point-wise comparison, 2-point correlation, and 3-point correlation. We also show that D$^3$M is able to accurately extrapolate far beyond its training data, and predict structure formation for significantly different cosmological parameters. Our study proves, for the first time, that deep learning is a practical and accurate alternative to approximate simulations of the gravitational structure formation of the Universe.

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 report on the highest energy photons from the Crab Nebula observed by the Tibet air shower array with the underground water-Cherenkov-type muon detector array. Based on the criterion of muon number measured in an air shower, we successfully suppress 99.92% of the cosmic-ray background events with energies $E>100$ TeV. As a result, we observed 24 photon-like events with $E>100$ TeV against 5.5 background events, which corresponds to 5.6$\sigma$ statistical significance. This is the first detection of photons with $E>100$ TeV from an astrophysical source.

We observed a field of $16'\times 16'$ in the star-forming region Pelican Nebula (IC 5070) at $BVRI$ wavelengths for 90 nights spread over one year in 2012-2013. More than 250 epochs in $VRI$-bands are used to identify and classify variables up to $V\sim 21$~mag. We present a catalogue of optical time-series photometry with periods, mean-magnitudes and classifications for 95 variable stars including 67 pre-main-sequence variables towards star-forming region IC 5070. The pre-main-sequence variables are further classified as candidate classical T Tauri and weak-line T Tauri stars based on their light curve variations and the locations on the color-color and color-magnitude diagrams using optical and infrared data together with Gaia DR2 astrometry. Classical T Tauri stars display variability amplitudes up to three times the maximum fluctuation in disk-free weak-line T Tauri stars, which show strong periodic variations. Short-term variability is missed in our photometry within single nights. Several classical T Tauri stars display long-lasting ($\geq 10$ days) single or multiple fading and brightening events up to a couple of magnitudes at optical wavelengths. The typical mass and age of the pre-main-sequence variables from the isochrone-fitting and spectral energy distributions are estimated to be $\le 1~M_\odot$ and $\sim 2$ Myr, respectively. We do not find any correlation between the optical amplitudes or periods with the physical parameters (mass and age) of pre-main-sequence stars.

We report the design concept and performance of a compact, light-weight, and economic imaging polarimeter, TRIPOL (the Triple Range Imager and POLarimeter), capable of simultaneous optical imagery and polarimetry. TRIPOL splits the beam from wavelength 400 to 830 nm into g'-, r'-, and i'-bands with two dichroic mirrors, and measures polarization with an achromatic half-waveplate and a wire-grid. The simultaneity makes TRIPOL a useful tool for small telescopes for photometry and polarimetry of time variable and wavelength dependent phenomena. TRIPOL is devised for a Cassegrain telescope of an aperture of ~1 m. This paper presents the engineering considerations of TRIPOL and compares the expected with the observed performance. Using the Lulin 1-m telescope and 100 seconds integration, the limiting magnitudes are g' ~ 19.0 mag, r' ~ 18.5 mag and i' ~ 18.0 mag with a signal-to-noise of 10, in agreement with design expectation. The instrumental polarization is measured to be ~ 0.3% at three bands. Two applications, one to the star-forming cloud IC 5146, and the other to the young variable GM Cep, are presented as demonstration.

We present the POL-2 850 $\mu$m linear polarization map of the Barnard 1 clump in the Perseus molecular cloud complex from the B-fields In STar-forming Region Observations (BISTRO) survey at the James Clerk Maxwell Telescope. We find a trend of decreasing polarization fraction as a function of total intensity, which we link to depolarization effects towards higher density regions of the cloud. We then use the polarization data at 850 $\mu$m to infer the plane-of-sky orientation of the large-scale magnetic field in Barnard 1. This magnetic field runs North-South across most of the cloud, with the exception of B1-c where it turns more East-West. From the dispersion of polarization angles, we calculate a turbulence correlation length of $5.0 \pm 2.5$ arcsec ($1500$ au), and a turbulent-to-total magnetic energy ratio of $0.5 \pm 0.3$ inside the cloud. We combine this turbulent-to-total magnetic energy ratio with observations of NH$_3$ molecular lines from the Green Bank Ammonia Survey (GAS) to estimate the strength of the plane-of-sky component of the magnetic field through the Davis-Chandrasekhar-Fermi method. With a plane-of-sky amplitude of $120 \pm 60$ $\mu$G and a criticality criterion $\lambda_c = 3.0 \pm 1.5$, we find that Barnard 1 is a supercritical molecular cloud with a magnetic field nearly dominated by its turbulent component.

Recently, the black hole X-ray binary (BHXB) Nova Muscae 1991 has been reported to be experiencing an extremely rapid orbital decay. So far, three BHXBs have anomalously high orbital period derivatives, which can not be interpreted by the standard stellar evolution theory. In this work, we investigate whether the resonant interaction between the binary and a surrounding circumbinary (CB) disk could produce the observed orbital period derivatives. Analytical calculations indicate that the observed orbital period derivatives of XTE J1118+480 and A0620-00 can originate from the tidal torque between the binary and a CB disk with a mass of $10^{-9}~\rm M_{\odot}$, which is approximately in agreement with the dust disk mass detected in these two sources. However, Nova Muscae 1991 was probably surrounded by a heavy CB disk with a mass of $10^{-7}~\rm M_{\odot}$. Based on the CB disk model and the anomalous magnetic braking theory, we simulate the evolution of the three BHXBs with intermediate-mass donor stars by using the MESA code. Our simulated results are approximately consistent with the observed donor star masses, orbital periods, and orbital-period derivatives. However, the calculated effective temperatures of the donor stars are higher than indicated by the observed spectral types of two sources.

We present the 850 $\mu$m polarization observations toward the IC5146 filamentary cloud taken using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) and its associated polarimeter (POL-2), mounted on the James Clerk Maxwell Telescope (JCMT), as part of the B-fields In STar forming Regions Observations (BISTRO). This work is aimed at revealing the magnetic field morphology within a core-scale ($\lesssim 1.0$ pc) hub-filament structure (HFS) located at the end of a parsec-scale filament. To investigate whether or not the observed polarization traces the magnetic field in the HFS, we analyze the dependence between the observed polarization fraction and total intensity using a Bayesian approach with the polarization fraction described by the Rice likelihood function, which can correctly describe the probability density function (PDF) of the observed polarization fraction for low signal-to-noise ratio (SNR) data. We find a power-law dependence between the polarization fraction and total intensity with an index of 0.56 in $A_V\sim$ 20--300 mag regions, suggesting that the dust grains in these dense regions can still be aligned with magnetic fields in the IC5146 regions. Our polarization maps reveal a curved magnetic field, possibly dragged by the contraction along the parsec-scale filament. We further obtain a magnetic field strength of 0.5$\pm$0.2 mG toward the central hub using the Davis-Chandrasekhar-Fermi method, corresponding to a mass-to-flux criticality of $\sim$ $1.3\pm0.4$ and an Alfv\'{e}nic Mach number of $<$0.6. These results suggest that gravity and magnetic field is currently of comparable importance in the HFS, and turbulence is less important.

We report the discovery of tidal structures around the intermediate-aged ($\sim$ 700--800~Myr), nearby ($\sim85$~pc) star cluster Coma Berenices. The spatial and kinematic grouping of stars is determined with the {\it Gaia} DR2 parallax and proper motion data, by a clustering analysis tool, \textsc{StarGO}, to map 5D parameters ($X, Y, Z$, $\mu_\alpha \cos\delta, \mu_\delta$) onto a 2D neural network. A leading and a trailing tails, each with an extension of $\sim50$~pc are revealed for the first time around this disrupting star cluster. The cluster members, totaling to $\sim115^{+5}_{-3}\,\rm {M_\odot}$, are clearly mass segregated, and exhibit a flat mass function with $\alpha \sim 0.79\pm0.16$, in the sense of $dN/dm \propto m^{-\alpha}$, where $N$ is the number of member stars and $m$ is stellar mass, in the mass range of $m=0.25$--$2.51~{\rm M_\odot}$. Within the tidal radius of $\sim$6.9~pc, there are 77 member candidates with an average position, i.e., as the cluster center, of R.A.= 186.8110~deg, and decl.= 25.8112~deg, and an average distance of 85.8~pc. Additional 120 member candidates reside in the tidal structures, i.e., outnumbering those in the cluster core. The expansion of escaping members lead to an anisotropy in the velocity field of the tidal tails. Our analysis also serendipitously uncovers an adjacent stellar group, part of which has been cataloged in the literature. We identify 218 member candidates, 10 times more than previously known. This star group is some 65~pc away from, and $\sim400$~Myr younger than, Coma Ber, but is already at the final stage of disruption.

We present the first release of the data and compact-source catalogue for the JCMT Large Program SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE). SCOPE consists of 850-um continuum observations of 1235 Planck Galactic Cold Clumps (PGCCs) made with the Submillimetre Common-User Bolometer Array 2 on the James Clerk Maxwell Telescope. These data are at an angular resolution of 14.4 arcsec, significantly improving upon the 353-GHz resolution of Planck at 5 arcmin, and allowing for a catalogue of 3528 compact sources in 558 PGCCs. We find that the detected PGCCs have significant sub-structure, with 61 per cent of detected PGCCs having 3 or more compact sources, with filamentary structure also prevalent within the sample. A detection rate of 45 per cent is found across the survey, which is 95 per cent complete to Planck column densities of $N_{H_{2}}$ $>$ 5 $\times$ 10$^{21}$ cm$^{-2}$. By positionally associating the SCOPE compact sources with YSOs, the star formation efficiency, as measured by the ratio of luminosity to mass, in nearby clouds is found to be similar to that in the more distant Galactic Plane, with the column density distributions also indistinguishable from each other.

We report 850~$\mu$m dust polarization observations of a low-mass ($\sim$12 $M_{\odot}$) starless core in the $\rho$ Ophiuchus cloud, Ophiuchus C, made with the POL-2 instrument on the James Clerk Maxwell Telescope (JCMT) as part of the JCMT B-fields In STar-forming Region Observations (BISTRO) survey. We detect an ordered magnetic field projected on the plane of sky in the starless core. The magnetic field across the $\sim$0.1~pc core shows a predominant northeast-southwest orientation centering between $\sim$40$^\circ$ to $\sim$100$^\circ$, indicating that the field in the core is well aligned with the magnetic field in lower-density regions of the cloud probed by near-infrared observations and also the cloud-scale magnetic field traced by Planck observations. The polarization percentage ($P$) decreases with an increasing total intensity ($I$) with a power-law index of $-$1.03 $\pm$ 0.05. We estimate the plane-of-sky field strength ($B_{\mathrm{pos}}$) using modified Davis-Chandrasekhar-Fermi (DCF) methods based on structure function (SF), auto-correlation (ACF), and unsharp masking (UM) analyses. We find that the estimates from the SF, ACF, and UM methods yield strengths of 103 $\pm$ 46 $\mu$G, 136 $\pm$ 69 $\mu$G, and 213 $\pm$ 115 $\mu$G, respectively. Our calculations suggest that the Ophiuchus C core is near magnetically critical or slightly magnetically supercritical (i.e. unstable to collapse). The total magnetic energy calculated from the SF method is comparable to the turbulent energy in Ophiuchus C, while the ACF method and the UM method only set upper limits for the total magnetic energy because of large uncertainties.

Blue Straggler Stars (BSSs) are observed in Galactic globular clusters and old open clusters. The radial distribution of BSSs has been used to diagnose the dynamical evolution of globular clusters. For the first time, with a reliable sample of BSSs identified with Gaia DR2, we conduct such an analysis for an open cluster. We identify members, including BSSs, of the oldest known Galactic open cluster Berkeley 17 with the Gaia DR2 proper motions and parallaxes. We study the radial distribution of the BSS population to understand the dynamical evolution of the cluster. We select cluster members to populate the colour magnitude diagram in the Gaia filters. Cluster parameters are derived using the brightest members. The BSSs and giant branch stars are identified, and their radial distributions are compared. The segregation of BSSs is also evaluated with respect to the giant branch stars using the Minimum Spanning Tree analysis. We determine Berkeley 17 to be at $3138.6^{+285.5}_{-352.9}$ pc. We find 23 BSS cluster members, only two of which were previously identified. We find a bimodal radial distribution of BSSs supported by findings from the MST method. The bimodal radial distribution of BSSs in Berkeley 17 indicates that they have just started to sink towards the cluster center, placing Berkeley 17 with globular clusters of intermediate dynamical age. This is the first such determination for an open cluster.

Individual highly magnified stars have been recently discovered at lookback times of more than half the age of the Universe, in lensed galaxies that straddle the critical curves of massive galaxy clusters. Having confirmed their detectability, it is now important to carry out systematic searches for them in order to establish their frequency, and in turn learn about the statistical properties of high-redshift stars and of the granularity of matter in the foreground deflector. Here we report the discovery of a highly magnified star at redshift $z = 0.94$ in a strongly lensed arc behind a Hubble Frontier Field galaxy cluster, MACS J0416.1-2403, discovered as part of a systematic archival search. The bright transient (dubbed "Warhol") was discovered in Hubble Space Telescope data taken on 2014 September 15 and 16. This single image faded over a period of two weeks, and observations taken on 2014 September 1 show that the duration of the microlensing event was at most four weeks in total. The light curve may also exhibit slow changes over a period of years consistent with the level of microlensing expected from stars responsible for the intracluster light (ICL) of the cluster. Optical and infrared observations taken near peak brightness can be fit by a stellar spectrum with moderate host-galaxy extinction. A blue supergiant matches the measured spectral energy distribution near peak, implying a temporary magnification of at least several thousand. While the spectrum of an O-type star would also fit the transient's spectral energy distribution, extremely luminous O-type stars are much less common than blue supergiants. The short timescale of the event and the estimated effective temperature indicate that the lensed source is an extremely magnified star.

UX Orionis stars (UXors) are Herbig Ae/Be or T Tauri stars exhibiting sporadic occultation of stellar light by circumstellar dust. GM\,Cephei is such a UXor in the young ($\sim4$~Myr) open cluster Trumpler\,37, showing prominent infrared excess, emission-line spectra, and flare activity. Our photometric monitoring (2008--2018) detects (1)~an $\sim$3.43~day period, likely arising from rotational modulation by surface starspots, (2)~sporadic brightening on time scales of days due to accretion, (3)~irregular minor flux drops due to circumstellar dust extinction, and (4)~major flux drops, each lasting for a couple of months with a recurrence time, though not exactly periodic, of about two years. The star experiences normal reddening by large grains, i.e., redder when dimmer, but exhibits an unusual "blueing" phenomenon in that the star turns blue near brightness minima. The maximum extinction during relatively short (lasting $\leq 50$~days) events, is proportional to the duration, a consequence of varying clump sizes. For longer events, the extinction is independent of duration, suggestive of a transverse string distribution of clumps. Polarization monitoring indicates an optical polarization varying $\sim3\%$--8$\%$, with the level anticorrelated with the slow brightness change. Temporal variation of the unpolarized and polarized light sets constraints on the size and orbital distance of the circumstellar clumps in the interplay with the young star and scattering envelope. These transiting clumps are edge-on manifestations of the ring- or spiral-like structures found recently in young stars with imaging in infrared of scattered light, or in submillimeter of thermalized dust emission.

Pan-STARRS1 has carried out a set of distinct synoptic imaging sky surveys including the $3\pi$ Steradian Survey and the Medium Deep Survey in 5 bands ($grizy_{P1}$). The mean 5$\sigma$ point source limiting sensitivities in the stacked 3$\pi$ Steradian Survey in $grizy_{P1}$ are (23.3, 23.2, 23.1, 22.3, 21.4) respectively. The upper bound on the systematic uncertainty in the photometric calibration across the sky is 7-12 millimag depending on the bandpass. The systematic uncertainty of the astrometric calibration using the Gaia frame comes from a comparison of the results with Gaia: the standard deviation of the mean and median residuals ($ \Delta ra, \Delta dec $) are (2.3, 1.7) milliarcsec, and (3.1, 4.8) milliarcsec respectively. The Pan-STARRS system and the design of the PS1 surveys are described and an overview of the resulting image and catalog data products and their basic characteristics are described together with a summary of important results. The images, reduced data products, and derived data products from the Pan-STARRS1 surveys are available to the community from the Mikulski Archive for Space Telescopes (MAST) at STScI.

A class of asteroids, called large super-fast rotators (large SFRs), have rotation periods shorter than 2 hours and diameters larger than ~ 0.3 km. They pose challenges to the usual interior rubble-pile structure unless a relatively high bulk density is assumed. So far, only six large SFRs have been found. Therefore, we present a survey of asteroid rotation period using Pan-STARRS 1 telescopes during 2016 October 26 to 31 to search more large SFRs and study their properties. A total of 876 reliable rotation periods are measured, among which seven are large SFRs, thereby increasing the inventory of known large SFRs. These seven newly discovered large SFRs have diverse colors and locations in the main asteroid belt, suggesting that the taxonomic tendency and the location preference in the inner main belt of the six perviously known large SFRs could be a bias due to various observational limits. Interestingly, five out of the seven newly discovered large SFRs are mid main-belt asteroids. Considering the rare discovery rates of large SFR in the previously similar surveys (Chang et al., 2015, 2016) and the survey condition in this work, the chance of detecting a large SFR in the inner main belt seems to be relatively low. This probably suggests that the inner main belt harbors less large SFRs than the mid main belt. From our survey, we also found the drop in number appearing at f > 5 rev/day on the spin-rate distribution for the outer main-belt asteroids of D < 3 km, which was reported for the inner and mid main belt by (Chang et al., 2015, 2016).

Luminous high-redshift quasars can be used to probe of the intergalactic medium (IGM) in the early universe because their UV light is absorbed by the neutral hydrogen along the line of sight. They help us to measure the neutral hydrogen fraction of the high-z universe, shedding light on the end of reionization epoch. In this paper, we present a discovery of a new quasar (PSO J006.1240+39.2219) at redshift $z=6.61\pm0.02$ from Panoramic Survey Telescope & Rapid Response System 1. Including this quasar, there are nine quasars above $z>6.5$ up to date. The estimated continuum brightness is $M_\text{1450}$=$-25.96\pm0.08$. PSO J006.1240+39.2219 has a strong Ly~$\alpha$ emission compared with typical low-redshift quasars, but the measured near-zone region size is $R_\text{NZ}=3.2\pm1.1$ proper megaparsecs, which is consistent with other quasars at z$\sim$6.

The binary T Tauri system JW 566 in the Orion Molecular Cloud underwent an energetic, short-lived flare observed at submillimetre wavelengths by the SCUBA-2 instrument on 26 November 2016 (UT). The emission faded by nearly 50% during the 31 minute integration. The simultaneous source fluxes averaged over the observation are 500 +/- 107 mJy/beam at 450 microns and 466 +/- 47 mJy/beam at 850 microns. The 850 micron flux corresponds to a radio luminosity of $L_{\nu}=8\times10^{19}$ erg/s/Hz, approximately one order of magnitude brighter (in terms of $\nu L_{\nu}$) than that of a flare of the young star GMR-A, detected in Orion in 2003 at 3mm. The event may be the most luminous known flare associated with a young stellar object and is also the first coronal flare discovered at sub-mm wavelengths. The spectral index between 450 microns and 850 microns of $\alpha = 0.11$ is broadly consistent with non-thermal emission. The brightness temperature was in excess of $6\times10^{4}$ K. We interpret this event to be a magnetic reconnection that energised charged particles to emit gyrosynchrotron/synchrotron radiation.

TeV emission from blazars can be used to probe the intergalactic magnetic fields and measure their intensity, coherence length, and helicity. Intergalactic magnetic fields deflect the electron-positron pairs produced by very-high-energy gamma-rays from blazars, resulting in broadened beams of cascade gamma-rays developing along the projected direction of the blazar jet. We present an analysis that uses for the first time the jet orientation of 12 high-synchrotron-peaked (HSP) BL~Lac-type blazars from VLBA radio observations to search for signatures of pair haloes in the Fermi-LAT data. Our search improves the sensitivity of previous studies by taking the asymmetry of the pair haloes into account, increasing the signal to noise and reducing the possibility of systematics in determining the point spread function. Although there is no significant detection, a hint for an offset halo with a global significance of 2$\sigma$ is found in the $30-300$ GeV energy range, corresponding to an intergalactic magnetic field with $\sim 10^{-15}$ Gauss, consistent with the inferred field from a prior study using stacked HSP BL Lac objects.

A multiwavelength campaign was organized to take place between March and July of 2012. Excellent temporal coverage was obtained with more than 25 instruments, including the MAGIC, FACT and VERITAS Cherenkov telescopes, the instruments on board the Swift and Fermi spacecraft, and the telescopes operated by the GASP-WEBT collaboration. Mrk 501 showed a very high energy (VHE) gamma-ray flux above 0.2 TeV of $\sim$0.5 times the Crab Nebula flux (CU) for most of the campaign. The highest activity occurred on 2012 June 9, when the VHE flux was $\sim$3 CU, and the peak of the high-energy spectral component was found to be at $\sim$2 TeV. This study reports very hard X-ray spectra, and the hardest VHE spectra measured to date for Mrk 501. The fractional variability was found to increase with energy, with the highest variability occurring at VHE, and a significant correlation between the X-ray and VHE bands. The unprecedentedly hard X-ray and VHE spectra measured imply that their low- and high-energy components peaked above 5 keV and 0.5 TeV, respectively, during a large fraction of the observing campaign, and hence that Mrk 501 behaved like an extreme high-frequency- peaked blazar (EHBL) throughout the 2012 observing season. This suggests that being an EHBL may not be a permanent characteristic of a blazar, but rather a state which may change over time. The one-zone synchrotron self-Compton (SSC) scenario can successfully describe the segments of the SED where most energy is emitted, with a significant correlation between the electron energy density and the VHE gamma-ray activity, suggesting that most of the variability may be explained by the injection of high-energy electrons. The one-zone SSC scenario used reproduces the behaviour seen between the measured X-ray and VHE gamma-ray fluxes, and predicts that the correlation becomes stronger with increasing energy of the X-rays.

We have identified stellar and substellar members in the nearby star cluster Coma Berenices, using photometry, proper motions, and distances of a combination of 2MASS, UKIDSS, URAT1, and {\it Gaia}/DR2 data. Those with {\it Gaia}/DR2 parallax measurements provide the most reliable sample to constrain the distance, averaging 86.7~pc with a dispersion 7.1~pc, and the age $\sim800$~Myr, of the cluster. This age is older than the 400--600~Myr commonly adopted in the literature. Our analysis, complete within 5deg of the cluster radius, leads to identification of 192 candidates, among which, after field contamination is considered, about 148 are true members. The members have $J\sim3$~mag to $\sim17.5$~mag, corresponding to stellar masses 2.3--0.06~$M_\odot$. The mass function of the cluster peaks around 0.3~$M_\odot$ and, in the sense of $dN/dm = m^{-\alpha}$, where $N$ is the number of members and $m$ is stellar mass, has a slope $\alpha\approx 0.49\pm0.03$ in the mass range 0.3--2.3~$M_\odot$. This is much shallower than that of the field population in the solar neighborhood. The slope $\alpha=-1.69\pm0.14$ from 0.3~$M_\odot$ to 0.06~$M_\odot$, the lowest mass in our sample. The cluster is mass segregated and has a shape elongated toward the Galactic plane. Our list contains nine substellar members, including three new discoveries of an M8, an L1 and an L4 brown dwarfs, extending from the previously known coolest members of late-M types to even cooler types.

We examine the possible influence of Earth-directed coronal mass ejections (ECMEs) on the Sun's shadow in the 3~TeV cosmic-ray intensity observed by the Tibet-III air shower (AS) array. We confirm a clear solar-cycle variation of the intensity deficit in the Sun's shadow during ten years between 2000 and 2009. This solar-cycle variation is overall reproduced by our Monte Carlo (MC) simulations of the Sun's shadow based on the potential field model of the solar magnetic field averaged over each solar rotation period. We find, however, that the magnitude of the observed intensity deficit in the Sun's shadow is significantly less than that predicted by MC simulations, particularly during the period around solar maximum when a significant number of ECMEs is recorded. The $\chi^2$ tests of the agreement between the observations and the MC simulations show that the difference is larger during the periods when the ECMEs occur, and the difference is reduced if the periods of ECMEs are excluded from the analysis. This suggests the first experimental evidence of the ECMEs affecting the Sun's shadow observed in the 3~TeV cosmic-ray intensity.

We present the results of dust emission polarization measurements of Ophiuchus-B (Oph-B) carried out using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) camera with its associated polarimeter (POL-2) on the James Clerk Maxwell Telescope (JCMT) in Hawaii. This work is part of the B-fields In Star-forming Region Observations (BISTRO) survey initiated to understand the role of magnetic fields in star formation for nearby star-forming molecular clouds. We present a first look at the geometry and strength of magnetic fields in Oph-B. The field geometry is traced over $\sim$0.2 pc, with clear detection of both of the sub-clumps of Oph-B. The field pattern appears significantly disordered in sub-clump Oph-B1. The field geometry in Oph-B2 is more ordered, with a tendency to be along the major axis of the clump, parallel to the filamentary structure within which it lies. The degree of polarization decreases systematically towards the dense core material in the two sub-clumps. The field lines in the lower density material along the periphery are smoothly joined to the large scale magnetic fields probed by NIR polarization observations. We estimated a magnetic field strength of 630$\pm$410 $\mu$G in the Oph-B2 sub-clump using a Davis-Chandeasekhar-Fermi analysis. With this magnetic field strength, we find a mass-to-flux ratio $\lambda$= 1.6$\pm$1.1, which suggests that the Oph-B2 clump is slightly magnetically supercritical.

The Hubble Space Telescope observation of the double white dwarf (WD) binary SDSS J125733.63$+$542850.5 reveal that the massive WD has a surface gravity log$g_1\sim8.7$ (which implies a mass of $M_1\sim1.06~{\rm M_\odot}$) and an effective temperature $T_1\sim13000$ K, while the effective temperature of the low-mass WD ($M_2<0.24$ {\rm M$_\odot$}) is $T_2\sim6400K$. Therefore, the massive and the low-mass WDs have a cooling age $\tau_1\sim1$ {\rm Gyr} and $\tau_2\geq5$ {\rm Gyr}, respectively. This is in contradiction with traditional binary evolution theory. In this Letter, we propose a strange dwarf scenario to explain the formation of this double WD binary. We assume that the massive WD is a strange dwarf originating from a phase transition in a $\sim1.1$ M$_\odot$ WD, which has experienced accretion and spin-down processes. Its high effective temperature could arise from the heating process during the phase transition. Our simulations suggest that the progenitor of SDSS J125733.63$+$542850.5 can be a binary system consisting of a $0.65~\rm M_{\odot}$ WD and a $1.5~\rm M_{\odot}$ main sequence star in a 1.492 day orbit. Especially, the secondary star (i.e., the progenitor of the low mass WD) is likely to have an ultra-low metallicity of $Z=0.0001$.

We present 850 $\mu$m imaging polarimetry data of the $\rho$ Oph-A core taken with the Submillimeter Common-User Bolometer Array-2 (SCUBA-2) and its polarimeter (POL-2), as part of our ongoing survey project, BISTRO (B-fields In STar forming RegiOns). The polarization vectors are used to identify the orientation of the magnetic field projected on the plane of the sky at a resolution of 0.01 pc. We identify 10 subregions with distinct polarization fractions and angles in the 0.2 pc $\rho$ Oph A core; some of them can be part of a coherent magnetic field structure in the $\rho$ Oph region. The results are consistent with previous observations of the brightest regions of $\rho$ Oph-A, where the degrees of polarization are at a level of a few percents, but our data reveal for the first time the magnetic field structures in the fainter regions surrounding the core where the degree of polarization is much higher ($> 5 \%$). A comparison with previous near-infrared polarimetric data shows that there are several magnetic field components which are consistent at near-infrared and submillimeter wavelengths. Using the Davis-Chandrasekhar-Fermi method, we also derive magnetic field strengths in several sub-core regions, which range from approximately 0.2 to 5 mG. We also find a correlation between the magnetic field orientations projected on the sky with the core centroid velocity components.

Jets and outflows are the first signposts of stellar birth. Emission in the H2 1-0S(1) line at 2.122 micron is a powerful tracer of shock excitation in these objects. Here we present the analysis of 2.0 x 0.8 square degrees data from the UK Widefield Infrared Survey for H2 (UWISH2) in the 1-0S(1) line to identify and characterize the outflows of the M17 complex. We uncover 48 probable outflows, of which, 93 per cent are new discoveries. We identified driving source candidates for 60 per cent of the outflows. Among the driving source candidate YSOs: 90 per cent are protostars and the reminder 10 per cent are Class II YSOs. Comparing with results from other surveys, we suggest that H2 emission fades very quickly as the objects evolve from protostars to pre-main-sequence stars. We fit SED models to 14 candidate outflow driving sources and conclude that the outflows of our sample are mostly driven by moderate-mass YSOs that are still actively accreting from their protoplanetary disc. We examined the spatial distribution of the outflows with the gas and dust distribution of the complex, and observed that the filamentary dark-cloud "M17SWex" located at the south-western side of the complex, is associated with a greater number of outflows. We find our results corroborate previous suggestions, that in the M17 complex, M17SWex is the most active site of star formation. Several of our newly identified outflow candidates are excellent targets for follow up studies to better understand very early phase of protostellar evolution.

We report the AGILE detection and the results of the multifrequency follow-up observations of a bright $\gamma$-ray flare of the blazar 3C 279 in June 2015. We use AGILE-GRID and Fermi-LAT $\gamma$-ray data, together with Swift-XRT, Swift-UVOT, and ground-based GASP-WEBT optical observations, including polarization information, to study the source variability and the overall spectral energy distribution during the $\gamma$-ray flare. The $\gamma$-ray flaring data, compared with as yet unpublished simultaneous optical data which allow to set constraints on the big blue bump disk luminosity, show very high Compton dominance values of $\sim 100$, with a ratio of $\gamma$-ray to optical emission rising by a factor of three in a few hours. The multi-wavelength behavior of the source during the flare challenges one-zone leptonic theoretical models. The new observations during the June 2015 flare are also compared with already published data and non-simultaneous historical 3C 279 archival data.

The formation and properties of star clusters formed at the edges of H II regions are poorly known. We study stellar content, physical conditions, and star formation processes around a relatively unknown young H II region IRAS 10427-6032, located in the southern outskirts of the Carina Nebula. We make use of near-IR data from VISTA, mid-IR from Spitzer and WISE, far-IR from Herschel, sub-mm from ATLASGAL, and 843 MHz radio-continuum data. Using multi-band photometry, we find a total of 5 Class I and 29 Class II young stellar object (YSO) candidates, most of which newly identified, in the 5'$\times$5' region centered on the IRAS source position. Modeling of the spectral energy distribution for selected YSO candidates using radiative transfer models shows that most of these candidates are intermediate mass YSOs in their early evolutionary stages. A majority of the YSO candidates are found to be coincident with the cold dense clump at the western rim of the H II region. Lyman continuum luminosity calculation using radio emission indicates the spectral type of the ionizing source to be earlier than B0.5-B1. We identified a candidate massive star possibly responsible for the H II region with an estimated spectral type B0-B0.5. The temperature and column density maps of the region constructed by performing pixel-wise modified blackbody fits to the thermal dust emission using the far-IR data show a high column density shell-like morphology around the H II region, and low column density (0.6 $\times$ 10$^{22}$ cm$^{-2}$) and high temperature ($\sim$21 K) matter within the H II region. Based on the morphology of the region in the ionized and the molecular gas, and the comparison between the estimated timescales of the H II region and the YSO candidates in the clump, we argue that the enhanced star-formation at the western rim of the H II region is likely due to compression by the ionized gas.

In this paper, we report our investigation of pulsar scintillation phenomena by monitoring PSR B0355$+$54 at 2.25 GHz for three successive months using \emph{Kunming 40-m radio telescope}. We have measured the dynamic spectrum, the two-dimensional correlation function, and the secondary spectrum. In those observations with high signal-to-noise ratio ($S/N\ge100$), we have detected the scintillation arcs, which are rarely observable using such a small telescope. The sub-microsecond scale width of the scintillation arc indicates that the transverse scale of structures on scattering screen is as compact as AU size. Our monitoring has also shown that both the scintillation bandwidth, timescale, and arc curvature of PSR B0355$+$54 were varying temporally. The plausible explanation would need to invoke multiple-scattering-screen or multiple-scattering-structure scenario that different screens or ray paths dominate the scintillation process at different epochs.

It has been suggested that accretion-induced collapse (AIC) is a non-negligible path for the formation of the observed neutron stars (NSs). An ONe white dwarf (WD) that accretes material from a He star may experience AIC process and eventually produce intermediate-mass binary pulsars (IMBPs), named as the ONe WD+He star scenario. Note that previous studies can only account for part of the observed IMBPs with short orbital periods. In this work, we investigate the evolution of about 900 ONe WD+He star binaries to explore the distribution of IMBPs. We found that the ONe WD+He star scenario could form IMBPs including pulsars with 5-340 ms spin periods and 0.75-1.38 Msun WD companions, in which the orbital periods range from 0.04 to 900 d. Compared with the 20 observed IMBPs, this scenario can cover the parameters of 13 sources in the final orbital period-WD mass plane and the Corbet diagram, most of which has short orbital periods. We found that the ONe WD+He star scenario can explain almost all the observed IMBPs with short orbital periods. This work can well match the observed parameters of PSR J1802-2124 (one of the two precisely observed IMBPs), providing a possible evolutional path for its formation. We also speculate that the compact companion of HD 49798 (a hydrogen depleted sdO6 star) may be not a NS based on the present work.

We analyze the Sun's shadow observed with the Tibet-III air shower array and find that the shadow's center deviates northward (southward) from the optical solar disc center in the "Away" ("Toward") IMF sector. By comparing with numerical simulations based on the solar magnetic field model, we find that the average IMF strength in the "Away" ("Toward") sector is $1.54 \pm 0.21_{\rm stat} \pm 0.20_{\rm syst}$ ($1.62 \pm 0.15_{\rm stat} \pm 0.22_{\rm syst}$) times larger than the model prediction. These demonstrate that the observed Sun's shadow is a useful tool for the quantitative evaluation of the average solar magnetic field.

We analyze results from the first eighteen months of monthly sub-mm monitoring of eight star-forming regions in the JCMT Transient Survey. In our search for stochastic variability in 1643 bright peaks, only the previously identified source, EC53, shows behavior well above the expected measurement uncertainty. Another four sources, two disks and two protostars, show moderately-enhanced standard deviations in brightness, as expected for stochastic variables. For the two protostars, this apparent variability is the result of single epochs that are much brighter than the mean. In our search for secular brightness variations that are linear in time, we measure the fractional brightness change per year for 150 bright peaks, fifty of which are protostellar. The ensemble distribution of slopes is well fit by a normal distribution with sigma ~ 0.023. Most sources are not rapidly brightening or fading in the sub-mm. Comparison against time-randomized realizations shows that the width of the distribution is dominated by the uncertainty in the individual brightness measurements of the sources. A toy model for secular variability reveals that an underlying Gaussian distribution of linear fractional brightness change sigma = 0.005 would be unobservable in the present sample, whereas an underlying distribution with sigma = 0.02 is ruled out. Five protostellar sources, 10% of the protostellar sample, are found to have robust secular measures deviating from a constant flux. The sensitivity to secular brightness variations will improve significantly with a larger time sample, with a factor of two improvement expected by the conclusion of our 36-month survey.

GRB 170817A is a weak short gamma-ray burst (GRB) accompanied by the gravitational wave (GW) event GW170817. It is believed, that an off beaming relativistic jet, produces this weak GRB. Here we use the $E_{\rm p,i}- E_{\rm iso}$ and $\Gamma - E_{\rm iso}$ relations to determine the Lorentz factor $\Gamma$ and the viewing angle from the edge of the jet $\theta'_{\rm obs}$ of GRB 170817A. Our result shows, $\Gamma = 13.4 ^{+9.8}_{-5.5}$ and $\theta'_{\rm obs} = 4.3^{+1.8}_{-1.5} \,^{\circ}$. This corresponds to the on-axis $E_{\rm p,i} = 415^{+361}_{-167}$ keV and $E_{\rm iso} = (2.4^{+1.6}_{-1.9}) \times 10^{47}$ ergs of the GRB, which is an intrinsically weak short GRB. Interestingly, the Doppler factor and the luminosity follow a universal relation from GRBs and blazars, which indicates they may share similar radiation mechanism.

Mergers of neutron star-neutron star (NS-NS) or neutron star-black hole (NS-BH) binaries are candidate sources of gravitational waves (GWs). At least a fraction of the merger remnant should be a stellar mass BH with a sub-relativistic ejecta. A collimated jet is launched via Blandford-Znajek mechanism from the central BH to trigger a short gamma-ray burst (SGRB). At the same time, a near-isotropic wind can be driven by the Blandford-Payne mechanism (BP). In previous work, additional energy injection to the ejecta from the BP mechanism was ignored, and radioactive decay has long been thought as the main source of energy for "kilonova" in the BH scenario. In this letter, we propose that the wind driven by BP mechanism from the new-born BH-disk can heat up and push the ejecta during the prompt emission phase or even at late time when there is mass fall-back. Such a BP-powered merger-nova could be bright at optical band even for a low-luminosity SGRB. The detection of a GW event with a merger product of BH, and accompanied by a bright merger-nova, would be a robust test of our model.

`Internal plateau' followed by a sharp decay is commonly seen in short gamma-ray burst (GRB) light curves. The plateau component is usually interpreted as the dipole emission from a supra-massive magnetar, and the sharp decay may imply the collapse of the magnetar to a black hole (BH). Fall-back accretion onto the new-born BH could produce long-lasting activities via the Blandford-Znajek (BZ) process. The magnetic flux accumulated near the BH would be confined by the accretion disks for a period of time. As the accretion rate decreases, the magnetic flux is strong enough to obstruct gas infall, leading to a magnetically-arrested disk (MAD). Within this scenario, we show that the BZ process could produce two types of typical X-ray light curves: type I exhibits a long-lasting plateau, followed by a power-law decay with slopes ranging from 5/3 to 40/9; type II shows roughly a single power-law decay with slope of 5/3. The former requires low magnetic filed strength, while the latter corresponds to relatively high values. We search for such signatures of the new-born BH from a sample of short GRBs with an internal plateau, and find two candidates: GRB 101219A and GRB 160821B, corresponding to type II and type I light curve, respectively. It is shown that our model can explain the data very well.

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.

Berkeley 59 is a nearby ($\sim$1 kpc) young cluster associated with the Sh2-171 H{\sc ii} region. We present deep optical observations of the central $\sim$2.5$\times$2.5 pc$^2$ area of the cluster, obtained with the 3.58-m Telescopio Nazionale Galileo. The $V$/($V$-$I$) color-magnitude diagram manifests a clear pre-main-sequence (PMS) population down to $\sim$ 0.2 M$_\odot$. Using the near-infrared and optical colors of the low-mass PMS members we derive a global extinction of A$_V$= 4 mag and a mean age of $\sim$ 1.8 Myr, respectively, for the cluster. We constructed the initial mass function and found that its global slopes in the mass ranges of 0.2 - 28 M$_\odot$ and 0.2 - 1.5 M$_\odot$ are -1.33 and -1.23, respectively, in good agreement with the Salpeter value in the solar neighborhood. We looked for the radial variation of the mass function and found that the slope is flatter in the inner region than in the outer region, indicating mass segregation. The dynamical status of the cluster suggests that the mass segregation is likely primordial. The age distribution of the PMS sources reveals that the younger sources appear to concentrate close to the inner region compared to the outer region of the cluster, a phenomenon possibly linked to the time evolution of star-forming clouds is discussed. Within the observed area, we derive a total mass of $\sim$ 10$^3$ M$_\odot$ for the cluster. Comparing the properties of Berkeley 59 with other young clusters, we suggest it resembles more to the Trapezium cluster.

The variability of young stellar objects is mostly driven by star-disk interactions. In long-term photometric monitoring of the accreting T Tauri star GI Tau, we detect extinction events with typical depths of $\Delta V \sim 2.5$ mag that last for days-to-months and often appear to occur stochastically. In 2014 - 2015, extinctions that repeated with a quasi-period of 21 days over several months is the first empirical evidence of slow warps predicted from MHD simulations to form at a few stellar radii away from the central star. The reddening is consistent with $R_V=3.85\pm0.5$ and, along with an absence of diffuse interstellar bands, indicates that some dust processing has occurred in the disk. The 2015 -- 2016 multi-band lightcurve includes variations in spot coverage, extinction, and accretion, each of which results in different traces in color-magnitude diagrams. This lightcurve is initially dominated by a month-long extinction event and return to the unocculted brightness. The subsequent light-curve then features spot modulation with a 7.03 day period, punctuated by brief, randomly-spaced extinction events. The accretion rate measured from $U$-band photometry ranges from $1.3\times10^{-8}$ to $1.1\times10^{-10}$ M$_\odot$ yr$^{-1}$ (excluding the highest and lowest 5% of high- and low- accretion rate outliers), with an average of $4.7 \times 10^{-9}$ M$_\odot$ yr$^{-1}$. A total of 50% of the mass is accreted during bursts of $>12.8\times10^{-9}$ M$_\odot$ yr${^{-1}}$, which indicates limitations on analyses of disk evolution using single-epoch accretion rates.

We present optical and near-infrared stellar polarization observations toward the dark filamentary clouds associated with IC5146. The data allow us to investigate the dust properties (this paper) and the magnetic field structure (Paper II). A total of 2022 background stars were detected in $R_{c}$-, $i'$-, $H$-, and/or $K$-bands to $A_V \lesssim 25$ mag. The ratio of the polarization percentage at different wavelengths provides an estimate of $\lambda_{max}$, the wavelength of peak polarization, which is an indicator of the small-size cutoff of the grain size distribution. The grain size distribution seems to significantly change at $A_V \sim$ 3 mag, where both the average and dispersion of $P_{R_c}/P_{H}$ decrease. In addition, we found $\lambda_{max}$ $\sim$ 0.6-0.9 $\mu$m for $A_V>2.5$ mag, which is larger than the $\sim$ 0.55 $\mu$m in the general ISM, suggesting that grain growth has already started in low $A_V$ regions. Our data also reveal that polarization efficiency (PE $\equiv P_{\lambda}/A_V$) decreases with $A_V$ as a power-law in $R_c$-, $i'$-, and $K$-bands with indices of -0.71$\pm$0.10, -1.23$\pm$0.10 and -0.53$\pm$0.09. However, $H$-band data show a power index change; the PE varies with $A_V$ steeply (index of -0.95$\pm$0.30) when $A_V < 2.88\pm0.67$ mag but softly (index of -0.25$\pm$0.06) for greater $A_V$ values. The soft decay of PE in high $A_V$ regions is consistent with the Radiative Aligned Torque model, suggesting that our data trace the magnetic field to $A_V \sim 20$ mag. Furthermore, the breakpoint found in $H$-band is similar to the $A_V$ where we found the $P_{R_c}/P_{H}$ dispersion significantly decreased. Therefore, the flat PE-$A_V$ in high $A_V$ regions implies that the power index changes result from additional grain growth.

Investigating variability at the earliest stages of low-mass star formation is fundamental in understanding how a protostar assembles mass. While many simulations of protostellar disks predict non-steady accretion onto protostars, deeper investigation requires robust observational constraints on the frequency and amplitude of variability events characterised across the observable SED. In this study, we develop methods to robustly analyse repeated observations of an area of the sky for submillimetre variability in order to determine constraints on the magnitude and frequency of deeply embedded protostars. We compare \mbox{850 $\mu$m} JCMT Transient Survey data with archival JCMT Gould Belt Survey data to investigate variability over 2-4 year timescales. Out of 175 bright, independent emission sources identified in the overlapping fields, we find 7 variable candidates, 5 of which we classify as \textit{Strong} and the remaining 2 as \textit{Extended} to indicate the latter are associated with larger-scale structure. For the \textit{Strong} variable candidates, we find an average fractional peak brightness change per year of |4.0|\% yr$^{-1}$ with a standard deviation of $2.7\%\mathrm{\:yr}^{-1}$. In total, 7\% of the protostars associated with \mbox{850 $\mu$m} emission in our sample show signs of variability. Four of the five \textit{Strong} sources are associated with a known protostar. The remaining source is a good follow-up target for an object that is anticipated to contain an enshrouded, deeply embedded protostar. In addition, we estimate the \mbox{850 $\mu$m} periodicity of the submillimetre variable source, EC 53, to be \mbox{567 $\pm$ 32 days} based on the archival Gould Belt Survey data.

The influence of magnetic fields (B-fields) in the formation and evolution of bipolar bubbles, due to the expanding ionization fronts (I-fronts) driven by the Hii regions that are formed and embedded in filamentary molecular clouds, has not been well-studied yet. In addition to the anisotropic expansion of I-fronts into a filament, B-fields are expected to introduce an additional anisotropic pressure which might favor expansion and propagation of I-fronts to form a bipolar bubble. We present results based on near-infrared polarimetric observations towards the central $\sim$8'$\times$8' area of the star-forming region RCW57A which hosts an Hii region, a filament, and a bipolar bubble. Polarization measurements of 178 reddened background stars, out of the 919 detected sources in the JHKs-bands, reveal B-fields that thread perpendicular to the filament long axis. The B-fields exhibit an hour-glass morphology that closely follows the structure of the bipolar bubble. The mean B-field strength, estimated using the Chandrasekhar-Fermi method, is 91$\pm$8 {\mu}G. B-field pressure dominates over turbulent and thermal pressures. Thermal pressure might act in the same orientation as those of B-fields to accelerate the expansion of those I-fronts. The observed morphological correspondence among the B-fields, filament, and bipolar bubble demonstrate that the B-fields are important to the cloud contraction that formed the filament, gravitational collapse and star formation in it, and in feedback processes. The latter include the formation and evolution of mid-infrared bubbles by means of B-field supported propagation and expansion of I-fronts. These may shed light on preexisting conditions favoring the formation of the massive stellar cluster in RCW57A.

We report here the results of searching for inelastic scattering of dark matter (initial and final state dark matter particles differ by a small mass splitting) with nucleon with the first 79.6-day of PandaX-II data (Run 9). We set the upper limits for the spin independent WIMP-nucleon scattering cross section up to a mass splitting of 300 keV/c$^2$ at two benchmark dark matter masses of 1 and 10 TeV/c$^2$.

The X-ray plateau followed by a steep decay ("internal plateau") has been observed in both long and short gamma-ray burst (GRBs), implying a millisecond magnetar operating in some GRBs. The sharp decay at the end of plateau, marking the abrupt cessation of the magnetar central engine, has been considered as the collapse of a supra-massive magnetar to a black hole (BH) when it spins down. If "internal plateau" is indeed the evidence of a magnetar central engine, a natural expectation is a signature from the new-born BH in some candidates. In this work, we find that GRB 070110 is a particular case, which shows a small X-ray bump following its "internal plateau". We interpret the plateau with a spin-down supra-massive magnetar and the X-ray bump with a fall-back BH accretion. This indicates that the new-born BH is likely active in some GRBs. Therefore, GRB 070110-like events may provide a further support to the magnetar central engine model and enable us to investigate the properties of the magnetar as well as the new-born BH.

We report a new search of weakly interacting massive particles (WIMPs) using the combined low background data sets in 2016 and 2017 from the PandaX-II experiment in China. The latest data set contains a new exposure of 77.1 live day, with the background reduced to a level of 0.8$\times10^{-3}$ evt/kg/day, improved by a factor of 2.5 in comparison to the previous run in 2016. No excess events were found above the expected background. With a total exposure of 5.4$\times10^4$ kg day, the most stringent upper limit on spin-independent WIMP-nucleon cross section was set for a WIMP with mass larger than 100 GeV/c$^2$, with the lowest exclusion at 8.6$\times10^{-47}$ cm$^2$ at 40 GeV/c$^2$.

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.

Most protostars have luminosities that are fainter than expected from steady accretion over the protostellar lifetime. The solution to this problem may lie in episodic mass accretion -- prolonged periods of very low accretion punctuated by short bursts of rapid accretion. However, the timescale and amplitude for variability at the protostellar phase is almost entirely unconstrained. In "A JCMT/SCUBA-2 Transient Survey of Protostars in Nearby Star Forming Regions", we are monitoring monthly with SCUBA-2 the sub-mm emission in eight fields within nearby (<500 pc) star forming regions to measure the accretion variability of protostars. The total survey area of ~1.6 sq.deg. includes ~105 peaks with peaks brighter than 0.5 Jy/beam (43 associated with embedded protostars or disks) and 237 peaks of 0.125-0.5 Jy/beam (50 with embedded protostars or disks). Each field has enough bright peaks for flux calibration relative to other peaks in the same field, which improves upon the nominal flux calibration uncertainties of sub-mm observations to reach a precision of ~2-3% rms, and also provides quantified confidence in any measured variability. The timescales and amplitudes of any sub-mm variation will then be converted into variations in accretion rate and subsequently used to infer the physical causes of the variability. This survey is the first dedicated survey for sub-mm variability and complements other transient surveys at optical and near-IR wavelengths, which are not sensitive to accretion variability of deeply embedded protostars.

We present the analysis of the morphological shape of Berkeley 17, the oldest known open cluster (~10 Gyr), using a probabilistic star counting of Pan-STARRS point sources, and confirm its core-tail shape, plus an antitail, previously detected with the 2MASS data. The stellar population, as diagnosed by the color-magnitude diagram and theoretical isochrones, shows many massive members in the cluster core, whereas there is a paucity of such members in both tails. This manifests mass segregation in this aged star cluster with the low-mass members being stripped away from the system. It has been claimed that Berkeley 17 is associated with an excessive number of blue straggler candidates. Comparison of nearby reference fields indicates that about half of these may be field contamination.

Sco X-1 is the brightest persistent X-ray in the sky. It is generally believed that Sco X-1 is a low-mass X-ray binary containing a neutron star accreting from a low-mass donor star where mass transfer is driven by the magnetic braking. However, the mass transfer rate predicted by the standard magnetic braking model is at least one order of magnitude lower than the one inferred by X-ray luminosity. In this work, we investigate whether this source could evolved from an intermediate-mass X-ray binary including Ap/Bp stars with a slightly strong magnetic field of 300 - 1000 G. The coupling between the magnetic field and an irradiation-driven wind induced by the X-ray flux from the accretor can yield a strong magnetic braking, which could give rise to a relatively high mass transfer rate. According to the observed orbital period, the mass transfer rate, the mass ratio, and the donor star spectral type, the progenitor of Sco X-1 should be an intermediate-mass X-ray binary including a 1.6 $-$ 1.8 $\rm M_{\odot}$ Ap/Bp donor star in a 1.3 $-$ 1.5 day orbit. Therefore, we propose that anomalous magnetic braking of Ap/Bp stars provides an alternative evolutionary channel to a part of luminous X-ray sources.

We make use of a catalog of 1600 Pan-STARRS1 groups produced by the probability friends-of-friends algorithm to explore how the galaxy properties, i.e. the specific star formation rate (SSFR) and quiescent fraction, depend on stellar mass and group-centric radius. The work is the extension of Lin et al. (2014). In this work, powered by a stacking technique plus a background subtraction for contamination removal, a finer correction and more precise results are obtained than in our previous work. We find that while the quiescent fraction increases with decreasing group-centric radius the median SSFRs of star-forming galaxies in groups at fixed stellar mass drop slightly from the field toward the group center. This suggests that the major quenching process in groups is likely a fast mechanism. On the other hand, a reduction in SSFRs by ~0.2 dex is seen inside clusters as opposed to the field galaxies. If the reduction is attributed to the slow quenching effect, the slow quenching process acts dominantly in clusters. In addition, we also examine the density-color relation, where the density is defined by using a sixth-nearest neighbor approach. Comparing the quiescent fractions contributed from the density and radial effect, we find that the density effect dominates over the massive group or cluster galaxies, and the radial effect becomes more effective in less massive galaxies. The results support mergers and/or starvation as the main quenching mechanisms in the group environment, while harassment and/or starvation dominate in clusters.

Distant luminous quasars provide important information on the growth of the first supermassive black holes, their host galaxies and the epoch of reionization. The identification of quasars is usually performed through detection of their Lyman-$\alpha$ line redshifted to $\sim$ 0.9 microns at z>6.5. Here, we report the discovery of a very Lyman-$\alpha$ luminous quasar, PSO J006.1240+39.2219 at redshift z=6.618, selected based on its red colour and multi-epoch detection of the Lyman-$\alpha$ emission in a single near-infrared band. The Lyman-$\alpha$-line luminosity of PSO J006.1240+39.2219 is unusually high and estimated to be 0.8$\times$10$^{12}$ Solar luminosities (about 3% of the total quasar luminosity). The Lyman-$\alpha$ emission of PSO J006.1240+39.2219 shows fast variability on timescales of days in the quasar rest frame, which has never been detected in any of the known high-redshift quasars. The high luminosity of the Lyman-$\alpha$ line, its narrow width and fast variability resemble properties of local Narrow-Line Seyfert 1 galaxies which suggests that the quasar is likely at the active phase of the black hole growth accreting close or even beyond the Eddington limit.

We present the first results from the B-fields In STar-forming Region Observations (BISTRO) survey, using the Sub-millimetre Common-User Bolometer Array 2 (SCUBA-2) camera, with its associated polarimeter (POL-2), on the James Clerk Maxwell Telescope (JCMT) in Hawaii. We discuss the survey's aims and objectives. We describe the rationale behind the survey, and the questions which the survey will aim to answer. The most important of these is the role of magnetic fields in the star formation process on the scale of individual filaments and cores in dense regions. We describe the data acquisition and reduction processes for POL-2, demonstrating both repeatability and consistency with previous data. We present a first-look analysis of the first results from the BISTRO survey in the OMC 1 region. We see that the magnetic field lies approximately perpendicular to the famous 'integral filament' in the densest regions of that filament. Furthermore, we see an 'hour-glass' magnetic field morphology extending beyond the densest region of the integral filament into the less-dense surrounding material, and discuss possible causes for this. We also discuss the more complex morphology seen along the Orion Bar region. We examine the morphology of the field along the lower-density north-eastern filament. We find consistency with previous theoretical models that predict magnetic fields lying parallel to low-density, non-self-gravitating filaments, and perpendicular to higher-density, self-gravitating filaments.

W4 is a giant H{\sc ii} region ionized by the OB stars of the cluster IC~1805. The H{\sc ii} region/cluster complex has been a subject of numerous investigations as it is an excellent laboratory for studying the feedback effect of massive stars on the surrounding region. However, the low-mass stellar content of the cluster IC~1805 remains poorly studied till now. With the aim to unravel the low-mass stellar population of the cluster, we present the results of a multiwavelength study based on deep optical data obtained with the Canada-France-Hawaii Telescope, infrared data from 2MASS, $Spitzer$ Space Telescope and X-ray data from $Chandra$ Space Telescope. The present optical dataset is complete enough to detect stars down to 0.2~M$_\odot$, which is the deepest optical observations so far for the cluster. We identified 384 candidate young stellar objects (YSOs; 101 Class I/II and 283 Class III) within the cluster using various colour-colour and colour-magnitude diagrams. We inferred the mean age of the identified YSOs to be $\sim$ 2.5 Myr and mass in the range 0.3 - 2.5 M$_\odot$. The mass function of our YSO sample has a power law index of -1.23 $\pm$ 0.23, close to the Salpeter value (-1.35), and consistent with those of other star-forming complexes. We explored the disk evolution of the cluster members and found that the diskless sources are relatively older compared to the disk bearing YSO candidates. We examined the effect of high-mass stars on the circumstellar disks and found that within uncertainties, the influence of massive stars on the disk fraction seems to be insignificant. We also studied the spatial correlation of the YSOs with the distribution of gas and dust of the complex to conclude that IC 1805 would have formed in a large filamentary cloud.

Apollo-type NEA (3200) Phaethon, classified at the B/F-type taxonomy, probably the main mass of the Phaethon-Geminid stream complex (PGC), can be the most metamorphic C-complex asteroid in our solar system, since it is heated up to ~1000 K by the solar radiation around its perihelion passages. Hence, its surface material may be easily decomposed in near-sun environment. Phaethon's spectrum exhibits extremely blue-slope in the VIS-NIR region (so-called Phaethon Blue). Another candidate large member of the PGC, (155140) 2005 UD, shows a B/F-type color, however with a C-type-like red color over its ~1/4 rotational part, which implies an exposition of less metamorphosed primordial internal structure of the PGC precursor by a splitting or breakup event long ago. If so, some rotational part of Phaethon should show the C-type color as well as 2005 UD. Hence, we carried out the time-series VIS-spectroscopic observations of Phaethon using 1-m telescope in order to detect such a signature. Also, R-band photometries were simultaneously performed in order to complement our spectroscopy. Consequently, we obtained a total of 68 VIS-spectrophotometric data, 78% of which show the B-type blue-color, as against the rest of 22% showing the C-type red-color. We successfully acquired rotationally time-resolved spectroscopic data, of which particular rotational phase shows a red-spectral slope as the C-type color, as 2005 UD does, suggesting longitudinal inhomogeneity on Phaethon's surface. We constrained this C-type red-colored area in the mid-latitude in Phaethon's southern hemisphere based on the rotationally time-resolved spectroscopy along with Phaethon's axial rotation state, of which size suggests the impact-induced origin of the PGC. We also surveyed the meteoritic analog of Phaethon's surface blue-color, and found thermally metamorphosed CI/CM chondrites as likely candidates.

Some short-period, detached binary systems have recently been reported to experience very rapid orbital decay, much faster than is expected from the angular-momentum loss caused by gravitational radiation alone. As these systems contain fully convective stars, magnetic braking is not believed to be operative, making the large orbital-period derivative puzzling. Here we explore whether a resonant interaction between the binary and a surrounding circumbinary (CB) disk could account for the observed orbital decay. Our calculations indicate that the observed orbital-period derivatives in seven detached binaries can be produced by the resonant interaction between the binary and a CB disk if the latter has a mass in the range of $10^{-4}-10^{-2}~ M_{\odot}$, which is of the same order as the inferred disk mass ($\sim2.4\times 10^{-4}~ M_{\odot}$) in the post-common-envelope binary NN Ser.

It is generally thought that an ultracompact X-ray binary is composed with a neutron star and a helium white dwarf donor star. As one of the most compact binaries, 4U 1820-303 in globular cluster NGC 6624 was predicted an orbital-period derivative at a rate of $\dot{P}/P\sim1.1\times10^{-7}$ ${\rm yr^{-1}}$ if the mass transfer is fully driven by gravitational radiation. However, the recent analysis of the 16 ${\rm yr}$ data from \textit{Rossi X-ray Timing Explorer} and other historical records yielded a negative orbital-period derivative in the past $35$ yr. In this work, we propose an evolutionary circumbinary (CB) disk model to account for this anomalous orbital-period derivative. 4U 1820-30 was known to undergo superbust events caused by runaway thermal nuclear burning on the neutron star. We assume that for a small fraction of the superbursts, part of the ejected material may form a CB disk around the binary. If the recurrence time of such superbursts is $\sim10,000$ yr and $\sim10$% of the ejected mass feeds a CB disk, the abrupt angular-momentum loss causes a temporary orbital shrink, and the donor's radius and its Roche-lobe radius do not keep in step. Driven by mass transfer and angular-momentum loss, the binary would adjust its orbital parameters to recover a new stable stage. Based on the theoretical analysis and numerical simulation, we find that the required feed mass at the CB disk is approximately $\sim 10^{-8}$ ${\rm M_{\odot}}$.

We present an extensive study of the BL Lac object Mrk 501 based on a data set collected during the multi-instrument campaign spanning from 2009 March 15 to 2009 August 1 which includes, among other instruments, MAGIC, VERITAS, Whipple 10-m, Fermi-LAT, RXTE, Swift, GASP-WEBT and VLBA. We find an increase in the fractional variability with energy, while no significant interband correlations of flux changes are found in the acquired data set. The higher variability in the very high energy (>100 GeV, VHE) gamma-ray emission and the lack of correlation with the X-ray emission indicate that the highest-energy electrons that are responsible for the VHE gamma-rays do not make a dominant contribution to the ~1 keV emission. Alternatively, there could be a very variable component contributing to the VHE gamma-ray emission in addition to that coming from the synchrotron self-Compton (SSC) scenarios. The space of SSC model parameters is probed following a dedicated grid-scan strategy, allowing for a wide range of models to be tested and offering a study of the degeneracy of model-to-data agreement in the individual model parameters. We find that there is some degeneracy in both the one-zone and the two-zone SSC scenarios that were probed, with several combinations of model parameters yielding a similar model-to-data agreement, and some parameters better constrained than others. The SSC model grid-scan shows that the flaring activity around 2009 May 22 cannot be modeled adequately with a one-zone SSC scenario, while it can be suitably described within a two-independent-zone SSC scenario. The observation of an electric vector polarization angle rotation coincident with the gamma-ray flare from 2009 May 1 resembles those reported previously for low frequency peaked blazars, hence suggesting that there are many similarities in the flaring mechanisms of blazars with different jet properties.

Recently, ultraluminous X-ray source (ULX) M82 X-2 has been identified to be an accreting neutron star, which has a $P=1.37$ s spin period, and is spinning up at a rate $\dot{P}=-2.0\times 10^{-10}~\rm s\,s^{-1}$. Interestingly, its isotropic X-ray luminosity $L_{\rm iso}=1.8\times 10^{40}~\rm erg\,s^{-1}$ during outbursts is 100 times the Eddington limit for a $1.4~\rm M_{\odot}$ neutron star. In this Letter, based on the standard accretion model we attempt to constrain the dipolar magnetic field of the pulsar in ULX M82 X-2. Our calculations indicate that the accretion rate at the magnetospheric radius must be super-Eddington during outbursts. To support such a super-Eddington accretion, a relatively high multipole field ($\ga 10^{13}$ G) near the surface of the accretor is invoked to produce an accreting gas column. However, our constraint shows that the surface dipolar magnetic field of the pulsar should be in the range of $1.0-3.5\times 10^{12}$ G. Therefore, our model supports that the neutron star in ULX M82 X-2 could be a low magnetic field magnetar (proposed by Tong) with a normal dipolar field ($\sim 10^{12}$ G) and relatively strong multipole field. For the large luminosity variations of this source, our scenario can also present a self-consistency interpretation.

Recent timing analysis reveals that the orbital period of the first discovered accreting millisecond pulsar SAX J1808.4-3658 is increasing at a rate $\dot{P}_{\rm orb}=(3.89\pm0.15)\times 10^{-12}~\rm s\,s^{-1}$, which is at least one order of magnitude higher than the value arising from the conservative mass transfer. An ejection of mass loss rate of $10^{-9}~\rm M_{\odot}{\rm yr}^{-1}$ from the donor star at the inner Lagrangian point during the quiescence state could interpret the observed orbital period derivative. However, it is unknown whether this source can offer such a high mass loss rate. In this work, we attempt to investigate an evolutionary channel towards SAX J1808.4-3658. Once the accretion disk becomes thermally and viscously unstable, the spin-down luminosity of the millisecond pulsar and the X-ray luminosity during outbursts are assumed to evaporate the donor star, and the resulting winds carry away the specific orbital angular momentum at the inner Lagrangian point. Our scenario could yield the observed orbital period, the orbital period derivative, the peak X-ray luminosity during outbursts. Low-mass X-ray binaries with a $1.0~\rm M_{\odot}$ donor star, and an orbital period in the range of 0.8 - 1.5 d, may be the progenitor of SAX J1808.4-3658. Our numerical calculations propose that the current donor star mass is $0.044 ~\rm M_{\odot}$, which is in approximately agreement with the minimum mass of the donor star. In addition, our scenario can also account for the formation of black widows or the diamond planet like PSR J1719-1438.

We observed RZ LMi, which is renowned for the extremely (~19d) short supercycle and is a member of a small, unusual class of cataclysmic variables called ER UMa-type dwarf novae, in 2013 and 2016. In 2016, the supercycles of this object substantially lengthened in comparison to the previous measurements to 35, 32, 60d for three consecutive superoutbursts. We consider that the object virtually experienced a transition to the novalike state (permanent superhumper). This observed behavior extremely well reproduced the prediction of the thermal-tidal instability model. We detected a precursor in the 2016 superoutburst and detected growing (stage A) superhumps with a mean period of 0.0602(1)d in 2016 and in 2013. Combined with the period of superhumps immediately after the superoutburst, the mass ratio is not as small as in WZ Sge-type dwarf novae, having orbital periods similar to RZ LMi. By using least absolute shrinkage and selection operator (Lasso) two-dimensional power spectra, we detected possible negative superhumps with a period of 0.05710(1)d. We estimated the orbital period of 0.05792d, which suggests a mass ratio of 0.105(5). This relatively large mass ratio is even above ordinary SU UMa-type dwarf novae, and it is also possible that the exceptionally high mass-transfer rate in RZ LMi may be a result of a stripped core evolved secondary which are evolving toward an AM CVn-type object.

In this work we report the detection of seven Neptune Trojans (NTs) in the Pan-STARRS 1 (PS1) survey. Five of these are new discoveries, consisting of four L4 Trojans and one L5 Trojan. Our orbital simulations show that the L5 Trojan stably librates for only several million years. This suggests that the L5 Trojan must be of recent capture origin. On the other hand, all four new L4 Trojans stably occupy the 1:1 resonance with Neptune for more than 1 Gyr. They can, therefore, be of primordial origin. Our survey simulation results show that the inclination width of the Neptune Trojan population should be between $7^{\circ}$ and $27^{\circ}$ at $>$ 95% confidence, and most likely $\sim 11^{\circ}$. In this paper, we describe the PS1 survey, the Outer Solar System pipeline, the confirming observations, and the orbital/physical properties of the new Neptune Trojans.

It is generally believed that Ultracompact X-ray binaries (UCXBs) evolved from binaries consisting of a neutron star accreting from a low-mass white dwarf or helium star where mass transfer is driven by gravitational radiation. However, the standard white-dwarf evolutionary channel cannot produce the relatively long-period ($40 - 60$\,min) UCXBs with high time-averaged mass-transfer rate. In this work, we explore an alternative evolutionary route toward UCXBs where the companions evolve from intermediate-mass Ap/Bp stars with an anomalously strong magnetic field ($100 - 10000$\,G). Including the magnetic braking caused by the coupling between the magnetic field and an irradiation-driven wind induced by the X-ray flux from the accreting component, we show that intermediate-mass X-ray binaries (IMXBs) can evolve into UCXBs. Using the \emph{MESA} code, we have calculated evolutionary sequences for a large number of IMXBs. The simulated results indicate that, for a small wind-driving efficiency $f=10^{-5}$, the anomalous magnetic braking can drive IMXBs to an ultra-short period of 11 min. Comparing our simulated results with the observed parameters of fifteen identified UCXBs, the anomalous magnetic braking evolutionary channel can account for the formation of seven and eight sources with $f=10^{-3}$, and $10^{-5}$, respectively. In particular, a relatively large value of $f$ can fit three of the long-period, persistent sources with high mass-transfer rate. Though the proportion of Ap/Bp stars in intermediate-mass stars is only 5\%, the lifetime of the UCXB phase is $\gtrsim$ 2 Gyr, producing a relatively high number of observable systems, making this an alternative evolutionary channel for the formation of UCXBs.

Although the majority of Centaurs are thought to have originated in the scattered disk, with the high-inclination members coming from the Oort cloud, the origin of the high inclination component of trans-Neptunian objects (TNOs) remains uncertain. We report the discovery of a retrograde TNO, which we nickname "Niku", detected by the Pan-STARRS 1 Outer Solar System Survey. Our numerical integrations show that the orbital dynamics of Niku are very similar to that of 2008 KV$_{42}$ (Drac), with a half-life of $\sim 500$ Myr. Comparing similar high inclination TNOs and Centaurs ($q > 10$ AU, $a < 100$ AU and $i > 60^\circ$), we find that these objects exhibit a surprising clustering of ascending node, and occupy a common orbital plane. This orbital configuration has high statistical significance: 3.8-$\sigma$. An unknown mechanism is required to explain the observed clustering. This discovery may provide a pathway to investigate a possible reservoir of high-inclination objects.

Recent timing observation constrained the braking index of the X-ray pulsar PSR J1640-4631 to be $n=3.15\pm0.03$, which is the highest value of all pulsars with measured braking indices so far. In this Letter, we investigate whether pulsar braking by combined between the magnetic dipole emission and the gravitational radiation might have a braking index greater than three. For conventional neutron star and low mass quark star candidates, the inferred ellipticities derived by the observed braking index are obviously much larger than the theoretical estimated maximum value. If PSR J1640-4631 is a low-mass neutron star with a mass of $0.1~ \rm M_{\odot}$, the inferred ellipticity can be approximately equal to the theoretical estimated maximum value. Because of the radio-quiet nature of this source, we employ the vacuum gap model developed by Ruderman and Sutherland to constrain the inclination angle to be $87.2-90^{\circ}$. Based on this, we propose that a low-mass neutron star with a large inclination angle can interpret the high braking index and the radio-quiet nature of this source. Future observations such as gravitational wave detection and long-term timing for this source are required to confirm or confute our scenario.

How black holes accrete surrounding matter is a fundamental, yet unsolved question in astrophysics. It is generally believed that matter is absorbed into black holes via accretion disks, the state of which depends primarily on the mass-accretion rate. When this rate approaches the critical rate (the Eddington limit), thermal instability is supposed to occur in the inner disc, causing repetitive patterns of large-amplitude X-ray variability (oscillations) on timescales of minutes to hours. In fact, such oscillations have been observed only in sources with a high mass accretion rate, such as GRS 1915+105. These large-amplitude, relatively slow timescale, phenomena are thought to have physical origins distinct from X-ray or optical variations with small amplitudes and fast ($\lesssim$10 sec) timescales often observed in other black hole binaries (e.g., XTE J1118+480 and GX 339-4). Here we report an extensive multi-colour optical photometric data set of V404 Cygni, an X-ray transient source containing a black hole of nine solar masses (and a conpanion star) at a distance of 2.4 kiloparsecs. Our data show that optical oscillations on timescales of 100 seconds to 2.5 hours can occur at mass-accretion rates more than ten times lower than previously thought. This suggests that the accretion rate is not the critical parameter for inducing inner-disc instabilities. Instead, we propose that a long orbital period is a key condition for these large-amplitude oscillations, because the outer part of the large disc in binaries with long orbital periods will have surface densities too low to maintain sustained mass accretion to the inner part of the disc. The lack of sustained accretion -- not the actual rate -- would then be the critical factor causing large-amplitude oscillations in long-period systems.

We present long-term photometric observations of the young open cluster IC 348 with a baseline time-scale of 2.4 yr. Our study was conducted with several telescopes from the Young Exoplanet Transit Initiative (YETI) network in the Bessel $R$ band to find periodic variability of young stars. We identified 87 stars in IC 348 to be periodically variable; 33 of them were unreported before. Additionally, we detected 61 periodic non-members of which 41 are new discoveries. Our wide field of view was the key to those numerous newly found variable stars. The distribution of rotation periods in IC 348 has always been of special interest. We investigate it further with our newly detected periods but we cannot find a statistically significant bimodality. We also report the detection of a close eclipsing binary in IC 348 composed of a low-mass stellar component ($M \gtrsim 0.09\,\mathrm{M}_{\odot}$) and a K0 pre-main sequence star ($M \approx 2.7\,\mathrm{M}_{\odot}$). Furthermore, we discovered three detached binaries among the background stars in our field of view and confirmed the period of a fourth one.

After several years of quiescence, the blazar CTA 102 underwent an exceptional outburst in 2012 September-October. The flare was tracked from gamma-ray to near-infrared frequencies, including Fermi and Swift data as well as photometric and polarimetric data from several observatories. An intensive GASP-WEBT collaboration campaign in optical and NIR bands, with an addition of previously unpublished archival data and extension through fall 2015, allows comparison of this outburst with the previous activity period of this blazar in 2004-2005. We find remarkable similarity between the optical and gamma-ray behaviour of CTA 102 during the outburst, with a time lag between the two light curves of ~1 hour, indicative of co-spatiality of the optical and gamma-ray emission regions. The relation between the gamma-ray and optical fluxes is consistent with the SSC mechanism, with a quadratic dependence of the SSC gamma-ray flux on the synchrotron optical flux evident in the post-outburst stage. However, the gamma-ray/optical relationship is linear during the outburst; we attribute this to changes in the Doppler factor. A strong harder-when-brighter spectral dependence is seen both the in gamma-ray and optical non-thermal emission. This hardening can be explained by convexity of the UV-NIR spectrum that moves to higher frequencies owing to an increased Doppler shift as the viewing angle decreases during the outburst stage. The overall pattern of Stokes parameter variations agrees with a model of a radiating blob or shock wave that moves along a helical path down the jet.

Aims. We study a sample of eight B[e] stars with uncertain evolutionary status to shed light on the origin of their circumstellar dust. Methods. We performed a diagnostic analysis on the spectral energy distribution beyond infrared wavelengths, and conducted a census of neighboring region of each target to ascertain its evolutionary status. Results. In comparison to pre-main sequence Herbig stars, these B[e] stars show equally substantial excess emission in the near-infrared, indicative of existence of warm dust, but much reduced excess at longer wavelengths, so the dusty envelopes should be compact in size. Isolation from star-forming regions excludes the possibility of their pre-main sequence status. Six of our targets, including HD 50138, HD 45677, CD-245721, CD-49 3441, MWC 623, and HD 85567, have been previously considered as FS CMa stars, whereas HD 181615/6 and HD 98922 are added to the sample by this work. We argue that the circumstellar grains of these isolated B[e] stars, already evolved beyond the pre-main sequence phase, should be formed in situ. This is in contrast to Herbig stars, which inherit large grains from parental molecular clouds. It has been thought that HD 98922, in particular, is a Herbig star because of its large infrared excess, but we propose it being in a more evolved stage. Because dust condenses out of stellar mass loss in an inside-out manner, the dusty envelope is spatially confined, and anisotropic mass flows, or anomalous optical properties of tiny grains, lead to the generally low line-of-sight extinction toward these stars.

We use high time cadence, high spectral resolution optical observations to detect excess H-alpha emission from the 2 - 3 Myr old weak lined T Tauri star PTFO8-8695. This excess emission appears to move in velocity as expected if it were produced by the suspected planetary companion to this young star. The excess emission is not always present, but when it is, the predicted velocity motion is often observed. We have considered the possibility that the observed excess emission is produced by stellar activity (flares), accretion from a disk, or a planetary companion; we find the planetary companion to be the most likely explanation. If this is the case, the strength of the H-alpha line indicates that the emission comes from an extended volume around the planet, likely fed by mass loss from the planet which is expected to be overflowing its Roche lobe.

We study the multi-band variability and correlations of the TeV blazar Mrk 421 on year time scales, which can bring additional insight on the processes responsible for its broadband emission. We observed Mrk 421 in the very high energy (VHE) gamma-ray range with the Cherenkov telescope MAGIC-I from March 2007 to June 2009 for a total of 96 hours of effective time after quality cuts. The VHE flux variability is quantified with several methods, including the Bayesian Block algorithm, which is applied to data from Cherenkov telescopes for the first time. The 2.3 year long MAGIC light curve is complemented with data from the Swift/BAT and RXTE/ASM satellites and the KVA, GASP-WEBT, OVRO, and Mets\"ahovi telescopes from February 2007 to July 2009, allowing for an excellent characterisation of the multi-band variability and correlations over year time scales. Mrk 421 was found in different gamma-ray emission states during the 2.3 year long observation period. Flares and different levels of variability in the gamma-ray light curve could be identified with the Bayesian Block algorithm. The same behaviour of a quiet and active emission was found in the X-ray light curves measured by Swift/BAT and the RXTE/ASM, with a direct correlation in time. The behaviour of the optical light curve of GASP-WEBT and the radio light curves by OVRO and Mets\"ahovi are different as they show no coincident features with the higher energetic light curves and a less variable emission. The fractional variability is overall increasing with energy. The comparable variability in the X-ray and VHE bands and their direct correlation during both high- and low-activity periods spanning many months show that the electron populations radiating the X-ray and gamma-ray photons are either the same, as expected in the Synchrotron-Self-Compton mechanism, or at least strongly correlated, as expected in electromagnetic cascades.

The ~2 Myr old classical T Tauri star CI Tau shows periodic variability in its radial velocity (RV) variations measured at infrared (IR) and optical wavelengths. We find that these observations are consistent with a massive planet in a ~9-day period orbit. These results are based on 71 IR RV measurements of this system obtained over 5 years, and on 26 optical RV measurements obtained over 9 years. CI Tau was also observed photometrically in the optical on 34 nights over ~one month in 2012. The optical RV data alone are inadequate to identify an orbital period, likely the result of star spot and activity induced noise for this relatively small dataset. The infrared RV measurements reveal significant periodicity at ~9 days. In addition, the full set of optical and IR RV measurements taken together phase coherently and with equal amplitudes to the ~9 day period. Periodic radial velocity signals can in principle be produced by cool spots, hot spots, and reflection of the stellar spectrum off the inner disk, in addition to resulting from a planetary companion. We have considered each of these and find the planet hypothesis most consistent with the data. The radial velocity amplitude yields an Msin(i) of ~8.1 M_Jup; in conjunction with a 1.3 mm continuum emission measurement of the circumstellar disk inclination from the literature, we find a planet mass of ~11.3 M_Jup, assuming alignment of the planetary orbit with the disk.

A variety of methods have been proposed to define and to quantify galaxy environments. While these techniques work well in general with spectroscopic redshift samples, their application to photometric redshift surveys remains uncertain. To investigate whether galaxy environments can be robustly measured with photo-z samples, we quantify how the density measured with the nearest neighbor approach is affected by photo-z uncertainties by using the Durham mock galaxy catalogs in which the 3D real-space environments and the properties of galaxies are exactly known. Furthermore, we present an optimization scheme in the choice of parameters used in the 2D projected measurements which yield the tightest correlation with respect to the 3D real-space environments. By adopting the optimized parameters in the density measurements, we show that the correlation between the 2D projected optimized density and real-space density can still be revealed, and the color-density relation is also visible out to $z \sim 0.8$ even for a photo-z uncertainty ($\sigma_{\Delta_{z}/(1+z)}$) up to 0.06. We find that at the redshift $0.3 < z < 0.5$ a deep ($i \sim 25$) photometric redshift survey with $\sigma_{\Delta_{z}/(1+z)} = 0.02$ yields a comparable performance of small-scale density measurement to a shallower $i \sim$ 22.5 spectroscopic sample with $\sim$ 10% sampling rate. Finally, we discuss the application of the local density measurements to the Pan-STARRS1 Medium Deep survey, one of the largest deep optical imaging surveys. Using data from $\sim5$ square degrees of survey area, our results show that it is possible to measure local density and to probe the color-density relation with 3$\sigma$ confidence level out to $z \sim 0.8$ in the PS-MDS. The color-density relation, however, quickly degrades for data covering smaller areas.

We report the detection of orbital modulation, a model solution, and X-ray properties of a newly discovered contact binary, 2MASS J11201034$-$2201340. We serendipitously found this X-ray point source outside the error ellipse when searching for possible X-ray counterparts of $\gamma$-ray millisecond pulsars among the unidentified objects detected by the {\it Fermi Gamma-ray Space Telescope}. The optical counterpart of the X-ray source (unrelated to the $\gamma$-ray source) was then identified using archival databases. The long-term CRTS survey detected a precise signal with a period of $P=0.28876208(56)$ days. A follow-up observation made by the SLT telescope of Lulin Observatory revealed the binary nature of the object. Utilizing archived photometric data of multi-band surveys, we construct the spectral energy distribution, which is well fitted by a K2V spectral template. The fitting result of the orbital profile using the Wilson--Devinney code suggests that 2MASS J11201034-2201340 is a short-period A-type contact binary and the more massive component has a cool spot. The X-ray emission was first noted in observations made by \swift\, then further confirmed and characterized by an \xmm\ observation. The X-ray spectrum can be described by a power law or thermal Bremsstrahlung. Unfortunately, we could not observe significant X-ray orbital modulation. Finally, according to the spectral energy distribution, this system is estimated to be 690 pc from Earth with a calculated X-ray intensity of $(0.7-1.5)\times 10^{30}$ erg s$^{-1}$, which is in the expected range of an X-ray emitting contact binary.

We present monitoring observations of the active T Tauri star RW Aur, from 2010 October to 2015 January, using optical high-resolution (R>10000) spectroscopy with CFHT-ESPaDOnS. Optical photometry in the literature shows bright, stable fluxes over most of this period, with lower fluxes (by 2-3 mag.) in 2010 and 2014. In the bright period our spectra show clear photospheric absorption, complicated variation in the Ca II 8542 A emission}profile shapes, and a large variation in redshifted absorption in the O I 7772 and 8446 A and He I 5876 A lines, suggesting unstable mass accretion during this period. In contrast, these line profiles are relatively uniform during the faint periods, suggesting stable mass accretion. During the faint periods the photospheric absorption lines are absent or marginal, and the averaged Li I profile shows redshifted absorption due to an inflow. We discuss (1) occultation by circumstellar material or a companion and (2) changes in the activity of mass accretion to explain the above results, together with near-infrared and X-ray observations from 2011-2015. Neither scenario can simply explain all the observed trends, and more theoretical work is needed to further investigate their feasibilities.

We report the discovery of 2 new Be stars, and re-identify one known Be star in the open cluster NGC 6830. Eleven H-alpha emitters were discovered using the H-alpha imaging photometry of the Palomar Transient Factory Survey. Stellar membership of the candidates was verified with photometric and kinematic information using 2MASS data and proper motions. The spectroscopic confirmation was carried out by using the Shane 3-m telescope at Lick observatory. Based on their spectral types, three H-alpha emitters were confirmed as Be stars with H-alpha equivalent widths > -10 Angstrom. Two objects were also observed by the new spectrograph SED-Machine on the Palomar 60 inch Telescope. The SED-Machine results show strong H-alpha emission lines, which are consistent with the results of the Lick observations. The high efficiency of the SED-Machine can provide rapid observations for Be stars in a comprehensive survey in the future.

The Karin family is a very young asteroid family created by an asteroid breakup 5.8 Myr ago. Since the members of this family probably have not experienced significant orbital or collisional evolution yet, it is possible that they still preserve properties of the original family-forming event in terms of their spin state. As we carried out a series of photometric observations of the Karin family asteroids, here we report an analysis result of lightcurves including the rotation period of eleven members. The mean rotation rate of the Karin family members turned out to be much lower than those of NEAs or smaller MBAs (diameter D<12 km), and even lower than that of larger MBAs (D>130 km). We investigated a correlation between the peak-to-peak variation magnitude reduced to zero solar phase angle and the rotation period of the eleven Karin family asteroids, and found a possible trend that elongated members have lower spin rate, and less elongated members have higher spin rate. However, this has to be confirmed by another series of future observations.

We present coordinated multiwavelength observations of the bright, nearby BL Lac object Mrk 421 taken in 2013 January-March, involving GASP-WEBT, Swift, NuSTAR, Fermi-LAT, MAGIC, VERITAS, and other collaborations and instruments, providing data from radio to very-high-energy (VHE) gamma-ray bands. NuSTAR yielded previously unattainable sensitivity in the 3-79 keV range, revealing that the spectrum softens when the source is dimmer until the X-ray spectral shape saturates into a steep power law with a photon index of approximately 3, with no evidence for an exponential cutoff or additional hard components up to about 80 keV. For the first time, we observed both the synchrotron and the inverse-Compton peaks of the spectral energy distribution (SED) simultaneously shifted to frequencies below the typical quiescent state by an order of magnitude. The fractional variability as a function of photon energy shows a double-bump structure which relates to the two bumps of the broadband SED. In each bump, the variability increases with energy which, in the framework of the synchrotron self-Compton model, implies that the electrons with higher energies are more variable. The measured multi-band variability, the significant X-ray-to-VHE correlation down to some of the lowest fluxes ever observed in both bands, the lack of correlation between optical/UV and X-ray flux, the low degree of polarization and its significant (random) variations, the short estimated electron cooling time, and the significantly longer variability timescale observed in the NuSTAR light curves point toward in-situ electron acceleration, and suggest that there are multiple compact regions contributing to the broadband emission of Mrk 421 during low-activity states.

We present time series photometry of 57 variable stars in the cluster region NGC 7380. The association of these variable stars to the cluster NGC 7380 has been established on the basis of two colour diagrams and colour-magnitude diagrams. Seventeen stars are found to be main-sequence variables, which are mainly B type stars and are classified as slowly pulsating B stars, $\beta$ Cep or $\delta$ Scuti stars. Some of them may belong to new class variables as discussed by Mowlavi et al. (2013) and Lata et al. (2014). Present sample also contains 14 pre-main-sequence stars, whose ages and masses are found to be mostly $\lesssim$ 5 Myr and range 0.60 $\lesssim M/M_{\odot} \lesssim$ 2.30 and hence should be T-Tauri stars. About half of the weak line T-Tauri stars are found to be fast rotators with a period of $\lesssim$ 2 days as compared to the classical T-Tauri stars. Some of the variables belong to the field star population.

Recent timing observation reported that the radio pulsar PSR J1734 - 3333 with a rotating period $P=1.17~\rm s$ is slowing down with a period derivative $\dot{P}=2.28\times 10^{-12}\rm s\,s^{-1}$. Its derived braking index $n=0.9 \pm 0.2$ is the lowest value among young radio pulsars with the measured braking indices. In this Letter, we attempt to investigate the influence of the braking torque caused by the interaction between the fall-back disk and the strong magnetic field of the pulsar on the spin evolution of PSR J1734 - 3333. Analytical result show that this braking torque is obviously far more than that by magnetic dipole radiation for pulsars with spin period of $> 0.1$ s, and play an important role during the spin-down of the pulsars. Our simulated results indicate that, for some typical neutron star parameters, the braking index and the period derivative approximately in agreement with the measured value of PSR J1734 - 3333 if the material inflow rate in the fallback disk is $2 \times 10^{17} \rm g\,s^{-1}$. In addition, our scenario can account for the measured braking indices of four young pulsars. However, our predicted X-ray luminosity are 1 -2 order of magnitude higher than the observation. We proposed that this discrepancy may originate from the instability of fall-back disk.

Recently, compact black hole X-ray binaries XTE J 1118+480 and A0620-00 have been reported to be experiencing a fast orbital period decay, which is two orders of magnitude higher than expected with gravitational wave radiation. Magnetic braking of an Ap/Bp star has been suggested to account for the period change when the surface magnetic field of the companion star $B_{\rm s}\ga 10^{4}$ G. However, our calculation indicates that anomalous magnetic braking cannot significantly contribute to the large orbital period decay rates observed in these two sources even if $B_{\rm s}\ga 10^{4}$ G. Observations have provided evidence that circumbinary disks around two compact black hole X-ray binaries may exist. Our analysis shows that, for some reasonable parameters, tidal torque between the circumbinary disk and the binary can efficiently extract the orbital angular momentum from the binary, and result in a large orbital period change rate. Based on the circumbinary disk model, we simulate the evolution of XTE J 1118+480 via a stellar evolution code. Our computations are approximatively in agreement with the observed data (the masses of two components, donor star radius, orbital period, and orbital period derivative). The mass transfer rate and circumbinary disk mass are obviously far greater than the inferred values from observations. Therefore, it seems that the circumbinary disk is unlikely to be the main cause of the rapid orbital decay observed in some compact black hole X-ray binaries.

We report on simultaneous broadband observations of the TeV-emitting blazar Markarian 501 between 1 April and 10 August 2013, including the first detailed characterization of the synchrotron peak with Swift and NuSTAR. During the campaign, the nearby BL Lac object was observed in both a quiescent and an elevated state. The broadband campaign includes observations with NuSTAR, MAGIC, VERITAS, the Fermi Large Area Telescope (LAT), Swift X-ray Telescope and UV Optical Telescope, various ground-based optical instruments, including the GASP-WEBT program, as well as radio observations by OVRO, Mets\"ahovi and the F-Gamma consortium. Some of the MAGIC observations were affected by a sand layer from the Saharan desert, and had to be corrected using event-by-event corrections derived with a LIDAR (LIght Detection And Ranging) facility. This is the first time that LIDAR information is used to produce a physics result with Cherenkov Telescope data taken during adverse atmospheric conditions, and hence sets a precedent for the current and future ground-based gamma-ray instruments. The NuSTAR instrument provides unprecedented sensitivity in hard X-rays, showing the source to display a spectral energy distribution between 3 and 79 keV consistent with a log-parabolic spectrum and hard X-ray variability on hour timescales. None (of the four extended NuSTAR observations) shows evidence of the onset of inverse-Compton emission at hard X-ray energies. We apply a single-zone equilibrium synchrotron self-Compton model to five simultaneous broadband spectral energy distributions. We find that the synchrotron self-Compton model can reproduce the observed broadband states through a decrease in the magnetic field strength coinciding with an increase in the luminosity and hardness of the relativistic leptons responsible for the high-energy emission.

We report the discovery of two methane dwarfs in the dark cloud L 1688 of the Rho Oph star-forming region. The two objects were among the T dwarf candidates with possible methane absorption and cool atmospheres, as diagnosed by infrared colors using deep WIRCam/CFHT HK plus CH4ON images, and IRAC/Spitzer c2d data. Follow-up spectroscopic observations with the FLAMINGOS-2/Gemini South confirmed the methane absorption at 1.6 micron. Compared with spectral templates of known T dwarfs in the field, i.e., of the old populations, Oph J162738-245240 (Oph-T3) is a T0/T1 type, whereas Oph J162645-241949 (Oph-T17) is consistent with a T3/T4 type in the H band but an L8/T1 in the K band. Compared with the BT-Settl model, both Oph-T3 and Oph-T17 are consistent with being cool, ~ 1000 K and ~ 900 K, respectively, and of low surface gravity, log(g) = 3.5. With an age no more than a couple Myr, these two methane dwarfs thereby represent the youngest T dwarfs ever confirmed. A young late L dwarf, OphJ162651-242110, was found serendipitously in our spectroscopic observations.

As part of our ongoing effort to investigate transit timing variations (TTVs) of known exoplanets, we monitored transits of the four exoplanets HAT-P-18b, HAT-P-19b, HAT-P-27b/WASP-40b and WASP-21b. All of them are suspected to show TTVs due to the known properties of their host systems based on the respective discovery papers. During the past three years 46 transit observations were carried out, mostly using telescopes of the Young Exoplanet Transit Initiative. The analyses are used to refine the systems orbital parameters. In all cases we found no hints for significant TTVs, or changes in the system parameters inclination, fractional stellar radius and planet to star radius ratio. However, comparing our results with those available in the literature shows that we can confirm the already published values.

We report on the results of a search for $\gamma$-ray pair halos with a stacking analysis of low-redshift blazars using data from the Fermi Large Area Telescope. For this analysis we used a number of a-priori selection criteria, including the spatial and spectral properties of the Fermi sources. The angular distribution of $\sim$ 1GeV photons around 24 stacked isolated high-synchrotron-peaked BL Lacs with redshift $z<0.5$ shows an excess over that of point-like sources. A statistical analysis yields a Bayes factor of $\mathrm{log}_{10}B_{10}>2$, providing evidence in favor of extended emission against the point-source hypothesis, consistent with expectations for pair halos produced in the IGMF with strength $B_{\mathrm{IGMF}}\sim 10^{-17}-10^{-15}\mathrm{G}$.

A flare from the TeV blazar Mrk 421, occurring in March 2010, was observed for 13 consecutive days from radio to very high energy (VHE, E > 100 GeV) gamma-rays with MAGIC, VERITAS, Whipple, FermiLAT, MAXI, RXTE, Swift, GASP-WEBT, and several optical and radio telescopes. We model the day-scale SEDs with one-zone and two-zone synchrotron self-Compton (SSC) models, investigate the physical parameters, and evaluate whether the observed broadband SED variability can be associated to variations in the relativistic particle population. Flux variability was remarkable in the X-ray and VHE bands while it was minor or not significant in the other bands. The one-zone SSC model can describe reasonably well the SED of each day for the 13 consecutive days. This flaring activity is also very well described by a two-zone SSC model, where one zone is responsible for the quiescent emission while the other smaller zone, which is spatially separated from the first one, contributes to the daily-variable emission occurring in X-rays and VHE gamma-rays. Both the one-zone SSC and the two-zone SSC models can describe the daily SEDs via the variation of only four or five model parameters, under the hypothesis that the variability is associated mostly to the underlying particle population. This shows that the particle acceleration and cooling mechanism producing the radiating particles could be the main one responsible for the broadband SED variations during the flaring episodes in blazars. The two-zone SSC model provides a better agreement to the observed SED at the narrow peaks of the low- and high-energy bumps during the highest activity, although the reported one-zone SSC model could be further improved by the variation of the parameters related to the emitting region itself ($\delta$, $B$ and $R$), in addition to the parameters related to the particle population.

Recent progress in the determination of both masses and radii of neutron stars are starting to place stringent constraints on the dense matter equation of state. In particular, new theoretical developments together with improved statistical tools seem to favor stellar radii that are significantly smaller than those predicted by models using purely nucleonic equations of state. Given that the underlying equation of state must also account for the observation of $2M_{\odot}$ neutron stars, theoretical approaches to the study of the dense matter equation of state are facing serious challenges. In response to this challenge, we compute in a model-independent way the underlying equation of state associated with an assumed mass-radius template similar to the "common radius" assumption used in recent studies. Once such a mass-radius template is adopted, the equation of state follows directly from the implementation of Lindblom's algorithm; assumptions on the nature or composition of the dense stellar core are not required. By analyzing mass-radius profiles with a maximum mass consistent with observation and common radii in the 8 to 11 km range, a lower limit on the stellar radius of a $1.4M_{\odot}$ neutron star of $R_{\rm NS}\!\gtrsim\!10.7$ km is obtained in order to prevent the equation of state from violating causality.

We performed a 4.5-month multi-instrument campaign (from radio to VHE gamma rays) on Mrk421 between January 2009 and June 2009, which included VLBA, F-GAMMA, GASP-WEBT, Swift, RXTE, Fermi-LAT, MAGIC, and Whipple, among other instruments and collaborations. Mrk421 was found in its typical (non-flaring) activity state, with a VHE flux of about half that of the Crab Nebula, yet the light curves show significant variability at all wavelengths, the highest variability being in the X-rays. We determined the power spectral densities (PSD) at most wavelengths and found that all PSDs can be described by power-laws without a break, and with indices consistent with pink/red-noise behavior. We observed a harder-when-brighter behavior in the X-ray spectra and measured a positive correlation between VHE and X-ray fluxes with zero time lag. Such characteristics have been reported many times during flaring activity, but here they are reported for the first time in the non-flaring state. We also observed an overall anti-correlation between optical/UV and X-rays extending over the duration of the campaign. The harder-when-brighter behavior in the X-ray spectra and the measured positive X-ray/VHE correlation during the 2009 multi-wavelength campaign suggests that the physical processes dominating the emission during non-flaring states have similarities with those occurring during flaring activity. In particular, this observation supports leptonic scenarios as being responsible for the emission of Mrk421 during non-flaring activity. Such a temporally extended X-ray/VHE correlation is not driven by any single flaring event, and hence is difficult to explain within the standard hadronic scenarios. The highest variability is observed in the X-ray band, which, within the one-zone synchrotron self-Compton scenario, indicates that the electron energy distribution is most variable at the highest energies.

Non-vanishing parity-odd correlators of gamma ray arrival directions observed by Fermi-LAT indicate the existence of a helical intergalactic magnetic field (Tashiro et al.2013). We successfully test this hypothesis using more stringent cuts of the data, Monte Carlo simulations with Fermi-LAT time exposure information, separate analyses for the northern and southern galactic hemispheres, and confirm predictions made in Tashiro & Vachaspati (2014). With some further technical assumptions, we show how to reconstruct the magnetic helicity spectrum from the parity-odd correlators.

We present Be star candidates in the open cluster NGC\,663, identified by H$\alpha$ imaging photometry with the Palomar Transient Factory Survey, as a pilot program to investigate how the Be star phenomena, the emission spectra, extended circumstellar envelopes, and fast rotation, correlate with massive stellar evolution. Stellar membership of the candidates was verified by 2MASS magnitudes and colors, and by proper motions. We discover 4 new Be stars and exclude one known Be star from being a member due to its inconsistent proper motions. The fraction of Be stars to member stars [N(Be)/N(members)] in NGC\,663 is 3.5\%. The spectral type of the 34 Be stars in NGC\,663 shows bimodal peaks at B0--B2 and B5--B7, which is consistent with the statistics in most star clusters. Additionally, we also discover 23 emission-line stars of different types, including non-member Be stars, dwarfs, and giants.

Mrk 501 is one of the brightest blazars at TeV energies and has been extensively studied since its first VHE detection in 1996. Our goal is to characterize in detail the source gamma-ray emission, together with the radio-to-X-ray emission, during the non-flaring (low) activity, which is less often studied than the occasional flaring (high) activity. We organized a multiwavelength (MW) campaign on Mrk 501 between March and May 2008. This multi-instrument effort included the most sensitive VHE gamma-ray instruments in the northern hemisphere, namely the imaging atmospheric Cherenkov telescopes MAGIC and VERITAS, as well as Swift, RXTE, the F-GAMMA, GASP-WEBT, and other collaborations and instruments. Mrk 501 was found to be in a low state of activity during the campaign, with a VHE flux in the range of 10%-20% of the Crab nebula flux. Nevertheless, significant flux variations were detected with various instruments, with a trend of increasing variability with energy. The broadband spectral energy distribution during the two different emission states of the campaign can be adequately described within the homogeneous one-zone synchrotron self-Compton model, with the (slightly) higher state described by an increase in the electron number density. This agrees with previous studies of the broadband emission of this source during flaring and non-flaring states. We report for the first time a tentative X-ray-to-VHE correlation during a low VHE activity. Although marginally significant, this positive correlation between X-ray and VHE, which has been reported many times during flaring activity, suggests that the mechanisms that dominate the X-ray/VHE emission during non-flaring-activity are not substantially different from those that are responsible for the emission during flaring activity.

The existence of a cosmological magnetic field could be revealed by the effects of non-trivial helicity on large scales. We evaluate a CP odd statistic, $Q$, using gamma ray data obtained from Fermi satellite observations at high galactic latitudes to search for such a signature. Observed values of $Q$ are found to be non-zero; the probability of a similar signal in Monte Carlo simulations is $\sim 0.2\%$. Contamination from the Milky Way does not seem to be responsible for the signal since it is present even for data at very high galactic latitudes. Assuming that the signal is indeed due to a helical cosmological magnetic field, our results indicate left-handed magnetic helicity and field strength $\sim 10^{-14}~{\rm G}$ on $\sim 10~{\rm Mpc}$ scales.

Theoretical predictions of physical observables often involve extrapolations to regions that are poorly constrained by laboratory experiments and astrophysical observations. Without properly quantified theoretical errors, such model predictions are of very limited utility. In this contribution we use maximum-likelihood estimation to compute theoretical errors and assess correlations between physical observables. We illustrate the power and elegance of these methods using examples of both pedagogical and realistic interest. In particular, we implement a gaussian approximation to the likelihood function to develop a new relativistic effective interaction constrained by ground-state properties of finite nuclei, their monopole response, and masses of neutron stars.

Among more than fifty blazars detected in very high energy (VHE, E>100GeV) gamma-rays, only three belong to the subclass of Flat Spectrum Radio Quasars (FSRQs). MAGIC observed FSRQ PKS 1510-089 in February-April 2012 during a high activity state in the high energy (HE, E>100 MeV) gamma-ray band observed by AGILE and Fermi. MAGIC observations result in the detection of a source with significance of 6.0 sigma. In agreement with the previous VHE observations of the source, we find no statistically significant variability during the MAGIC observations in daily, weekly or monthly time scales. The other two known VHE FSRQs have shown daily scale to sub-hour variability. We study the multifrequency behaviour of the source at the epoch of MAGIC observation, collecting quasi-simultaneous data at radio and optical (GASP-WEBT and F-Gamma collaborations, REM, Steward, Perkins, Liverpool, OVRO and VLBA telescopes), X-ray (Swift satellite) and HE gamma-ray frequencies. The gamma-ray SED combining AGILE, Fermi and MAGIC data joins smoothly and shows no hint of a break. The multifrequency light curves suggest a common origin for the millimeter radio and HE gamma-ray emission and the HE gamma-ray flaring starts when the new component is ejected from the 43GHz VLBA core. The quasi-simultaneous multifrequency SED is modelled with a one-zone inverse Compton model. We study two different origins of the seed photons for the inverse Compton scattering, namely the infra-red torus and a slow sheath surrounding the jet around the VLBA core. Both models fit the data well. However, the fast HE gamma-ray variability requires that within the modelled large emitting region, there must exist more compact regions. We suggest that these observed signatures would be most naturally explained by a turbulent plasma flowing at a relativistic speed down the jet and crossing a standing conical shock.

Bright-rimmed clouds (BRCs), illuminated and shaped by nearby OB stars, are potential sites of recent/ongoing star formation. Here we present an optical and infrared photometric study of three BRCs: BRC 5, BRC 7 and BRC 39 to obtain a census of the young stellar population, thereby inferring the star formation scenario, in these regions. In each BRC, the Class I sources are found to be located mostly near the bright rim or inside the cloud, whereas the Class II sources are preferentially outside, with younger sources closer to the rim. This provides strong support to sequential star formation triggered by radiation driven implosion due to the UV radiation. Moreover, each BRC contains a small group of young stars being revealed at its head, as the next-generation stars. In particular, the young stars at the heads of BRC 5 and BRC 7 are found to be intermediate/high mass stars, which, under proper conditions, may themselves trigger further star birth, thereby propagating star formation out to long distances.

It is difficult for the intermediate-mass X-ray binaries (IMXBs) evolutionary channel to form intermediate-mass binary pulsars (IMBPs) with a short orbital period (less than 3 d) via stable mass transfer. The main reason is that the magnetic braking mechanisms are generally thought not to work for donor stars with a mass of greater than 1.5 $\rm M_{\odot}$ in the canonical model. However, some intermediate-mass stars have anomalously strong magnetic fields (about 100 -- 10000 G), i. e. so-called Ap or Bp stars. With the coupling between the magnetic field and the irradiation-driven wind from the surface of Ap/Bp stars, a plausible magnetic braking mechanism should be expected. In this work, we attempt to investigate if IMXBs with Ap/Bp stars can produce IMBPs with a short orbital period (less than 3 d) by such an anomalous magnetic braking mechanism. Using a stellar evolution code, we have simulated the evolution of a large number of IMXBs consisting of a NS and an Ap/Bp star. For the spin evolution of the NS, we consider the accretion torque, the propeller torque, and the spin-down torque caused by the interaction between the magnetic field and the accretion disc. The calculated results show that, employing anomalous magnetic braking of Ap/Bp stars, IMXBs can evolve into compact IMBPs with short orbital periods of less than 3 d. However, there exists significant discrepancy between the spin periods of IMBPs in our simulated results and those observed.

Aims: We performed a detailed photometric analysis of the lensed system UM673 (Q0142-100) and an analysis of the tentative lens models. Methods: High-resolution adaptive optics images of UM673 taken with the Subaru telescope in the H band were examined. We also analysed the J, H and K-band observational data of UM673 obtained with the 1.3m telescope at the CTIO observatory. Results: We present photometry of quasar components A and B of UM673, the lens, and the nearby bright galaxy using H-band observational data obtained with the Subaru telescope. Based on the CTIO observations of UM673, we also present J- and H-band photometry and estimates of the J, H and K-band flux ratios between the two UM673 components in recent epochs. The near-infrared fluxes of the A and B components of UM673 and their published optical fluxes are analysed to measure extinction properties of the lensing galaxy. We estimate the extinction-corrected flux ratio between components A and B to be about 2.14 mag. We discuss lens models for the UM673 system constrained with the positions of the UM673 components, their flux ratio, and the previously measured time delay

GRB 131231A was detected by the Large Area Telescope onboard Fermi Space Gamma-ray Telescope. The high energy gamma-ray ($> 100$ MeV) afterglow emission spectrum is $F_\nu \propto \nu^{-0.54\pm0.15}$ in the first $\sim 1300$ s after the trigger and the most energetic photon has an energy $\sim 62$ GeV arriving at $t\sim 520$ s. With reasonable parameters of the GRB outflow as well as the density of the circum-burst medium, the synchrotron radiation of electrons or protons accelerated at an external forward shock have difficulty accounting for the data. The synchrotron self-Compton radiation of the forward shock-accelerated electrons, instead, can account for both the spectrum and temporal behavior of the GeV afterglow emission. We also show that the prospect for detecting GRB 131231A$-$like GRBs with Cherenkov Telescope Array (CTA) is promising.

We report our investigation of the first transiting planet candidate from the YETI project in the young (~4 Myr old) open cluster Trumpler 37. The transit-like signal detected in the lightcurve of the F8V star 2M21385603+5711345 repeats every 1.364894+/-0.000015 days, and has a depth of 54.5+/-0.8 mmag in R. Membership to the cluster is supported by its mean radial velocity and location in the color-magnitude diagram, while the Li diagnostic and proper motion are inconclusive in this regard. Follow-up photometric monitoring and adaptive optics imaging allow us to rule out many possible blend scenarios, but our radial-velocity measurements show it to be an eclipsing single-lined spectroscopic binary with a late-type (mid-M) stellar companion, rather than one of planetary nature. The estimated mass of the companion is 0.15-0.44 solar masses. The search for planets around very young stars such as those targeted by the YETI survey remains of critical importance to understand the early stages of planet formation and evolution.

P/2011 S1 (Gibbs) is an outer solar system comet or active Centaur with a similar orbit to that of the famous 29P/Schwassmann-Wachmann 1. P/2011 S1 (Gibbs) has been observed by the Pan-STARRS 1 (PS1) sky survey from 2010 to 2012. The resulting data allow us to perform multi-color studies of the nucleus and coma of the comet. Analysis of PS1 images reveals that P/2011 S1 (Gibbs) has a small nucleus $< 4$ km radius, with colors $g_{P1}-r_{P1} = 0.5 \pm 0.02$, $r_{P1}-i_{P1} = 0.12 \pm 0.02$ and $i_{P1}-z_{P1} = 0.46 \pm 0.03$. The comet remained active from 2010 to 2012, with a model-dependent mass-loss rate of $\sim100$ kg s$^{-1}$. The mass-loss rate per unit surface area of P/2011 S1 (Gibbs) is as high as that of 29P/Schwassmann-Wachmann 1, making it one of the most active Centaurs. The mass-loss rate also varies with time from $\sim 40$ kg s$^{-1}$ to 150 kg s$^{-1}$. Due to its rather circular orbit, we propose that P/2011 S1 (Gibbs) has 29P/Schwassmann-Wachmann 1-like outbursts that control the outgassing rate. The results indicate that it may have a similar surface composition to that of 29P/Schwassmann-Wachmann 1. Our numerical simulations show that the future orbital evolution of P/2011 S1 (Gibbs) is more similar to that of the main population of Centaurs than to that of 29P/Schwassmann-Wachmann 1. The results also demonstrate that P/2011 S1 (Gibbs) is dynamically unstable and can only remain near its current orbit for roughly a thousand years.

Membership identification is the first step to determine the properties of a star cluster. Low-mass members in particular could be used to trace the dynamical history, such as mass segregation, stellar evaporation, or tidal stripping, of a star cluster in its Galactic environment. We identified member candidates with stellar masses $\sim$0.11--2.4 Msun of the intermediate-age Praesepe cluster (M44), by using Pan-STARRS and 2MASS photometry, and PPMXL proper motions. Within a sky area of 3 deg radius, 1040 candidates are identified, of which 96 are new inclusions. Using the same set of selection criteria on field stars led to an estimate of a false positive rate 16%, suggesting 872 of the candidates being true members. This most complete and reliable membership list allows us to favor the BT-Settl model in comparison with other stellar models. The cluster shows a distinct binary track above the main sequence, with a binary frequency of 20--40%, and a high occurrence rate of similar mass pairs. The mass function is consistent with that of the disk population but shows a deficit of members below 0.3 solar masses. A clear mass segregation is evidenced, with the lowest-mass members in our sample being evaporated from this disintegrating cluster.

The XENON100 dark matter experiment uses liquid xenon in a time projection chamber (TPC) to measure xenon nuclear recoils resulting from the scattering of dark matter Weakly Interacting Massive Particles (WIMPs). In this paper, we report the observation of single-electron charge signals which are not related to WIMP interactions. These signals, which show the excellent sensitivity of the detector to small charge signals, are explained as being due to the photoionization of impurities in the liquid xenon and of the metal components inside the TPC. They are used as a unique calibration source to characterize the detector. We explain how we can infer crucial parameters for the XENON100 experiment: the secondary-scintillation gain, the extraction yield from the liquid to the gas phase and the electron drift velocity.

Using a large sample of field and group galaxies drawn from the Pan-STARRS1 Medium-Deep Survey, we present the specific star formation rate (SSFR) - stellar mass (M*) relation, as well as the quiescent fraction versus M* relation in different environments. We confirm that the fraction of quiescent galaxies is strongly dependent on environment at a fixed stellar mass, but that the amplitude and the slope of the star-forming sequence is similar between the field and groups: in other words, the SSFR-density relation at a fixed stellar mass is primarily driven by the change in the star-forming and quiescent fractions between different environments rather than a global suppression in the star formation rate for the star-forming population. However, when we restrict our sample to the cluster-scale environments ($M>10^{14}M_{solar}$), we find a global reduction in the SSFR of the star forming sequence of $17\%$ at 4$\sigma$ confidence as opposed to its field counterpart. After removing the stellar mass dependence of the quiescent fraction seen in field galaxies, the excess in the quiescent fraction due to the environment quenching in groups and clusters is found to increase with stellar mass. We argue that these results are in favor of galaxy mergers to be the primary environment quenching mechanism operating in galaxy groups whereas strangulation is able to reproduce the observed trend in the environment quenching efficiency and stellar mass relation seen in clusters. Our results also suggest that the relative importance between mass quenching and environment quenching depends on stellar mass -- the mass quenching plays a dominant role in producing quiescent galaxies for more massive galaxies, while less massive galaxies are quenched mostly through the environmental effect, with the transition mass around $1-2\times10^{10}M_{solar}$ in the group/cluster environment. (abridged)

Continuing the attempt to understand the properties of the stellar content in the young cluster NGC 1893 we have carried out a comprehensive multi-wavelength study of the region. The present study focuses on the X-ray properties of T-Tauri Stars (TTSs) in the NGC 1893 region. We found a correlation between the X-ray luminosity, $L_X$, and the stellar mass (in the range 0.2$-$2.0 \msun) of TTSs in the NGC 1893 region, similar to those reported in some other young clusters, however the value of the power-law slope obtained in the present study ($\sim$ 0.9) for NGC 1893 is smaller than those ($\sim$1.4 - 3.6) reported in the case of TMC, ONC, IC 348 and Chameleon star forming regions. However, the slope in the case of Class III sources (Weak line TTSs) is found to be comparable to that reported in the case of NGC 6611 ($\sim$ 1.1). It is found that the presence of circumstellar disks has no influence on the X-ray emission. The X-ray luminosity for both CTTSs and WTTSs is found to decrease systematically with age (in the range $\sim $ 0.4 Myr - 5 Myr). The decrease of the X-ray luminosity of TTSs (slope $\sim$ -0.6) in the case of NGC 1893 seems to be faster than observed in the case of other star-forming regions (slope -0.2 to -0.5). There is indication that the sources having relatively large NIR excess have relatively lower $L_X$ values. TTSs in NGC 1893 do not follow the well established X-ray activity - rotation relation as in the case of main-sequence stars.