M31N 2017-09a is a classical nova and was observed for some 160 days following its initial eruption, during which time it underwent a number of bright secondary outbursts. The light-curve is characterized by continual variation with excursions of at least 0.5 magnitudes on a daily time-scale. The lower envelope of the eruption suggests that a single power-law can describe the decline rate. The eruption is relatively long with $t_2 = 111$, and $t_3 = 153$ days.

Lithium depletion and enrichment in the cosmos is not yet well understood. To help tighten constraints on stellar and Galactic evolution models, we present the largest high-resolution analysis of Li abundances A(Li) to date, with results for over 100 000 GALAH field stars spanning effective temperatures $5900\,\mathrm{K} \lesssim \rm{T_{eff}} \lesssim7000\,\mathrm{K}$ and metallicities $-3 \lesssim \rm[Fe/H] \lesssim +0.5$. We separated these stars into two groups, on the warm and cool side of the so-called Li-dip, a localised region of the Kiel diagram wherein lithium is severely depleted. We discovered that stars in these two groups show similar trends in the A(Li)-[Fe/H] plane, but with a roughly constant offset in A(Li) of 0.4 dex, the warm group having higher Li abundances. At $\rm[Fe/H]\gtrsim-0.5$, a significant increasing in Li abundance with increasing metallicity is evident in both groups, signalling the onset of significant Galactic production. At lower metallicity, stars in the cool group sit on the Spite plateau, showing a reduced lithium of around 0.4 dex relative to the primordial value predicted from Big Bang nucleosynthesis (BBN). However, stars in the warm group between [Fe/H] = -1.0 and -0.5, form an elevated plateau that is largely consistent with the BBN prediction. This may indicate that these stars in fact preserve the primordial Li produced in the early Universe.

The first charge-coupled device photometric light curves in B, V, R${_c}$, and I${_c}$ bands of the short-period contact binary star OQ Cam are presented. Through analyzing the light curves with the Wilson-Devinney code, it is found that OQ Cam is a W-type shallow contact binary star with a contact degree of f =20.2%. Based on the O-C analyzing, the orbit period displays a long term increasing with a rate of $dP/dt=4.40\times10{^{-7}} day \cdot year{^{-1}}$. The increasing in orbit period can be interpreted by mass transfer from the less massive component to the more massive one. As the orbit period increasing, OQ Cam would evolve from the present shallow contact configuration to a none contact stage. So it may be a potential candidate to confirm the thermal relaxation oscillation theory.

In this paper, we present the first light curve synthesis and orbital period change analysis of nine contact binaries around the short period limit. It is found that all these systems are W-subtype contact binaries. One of them is a medium contact system while the others are shallow contact ones. Four of them manifest obvious O'Connell effect explained by a dark spot or hot spot on one of the component stars. Third light was detected in three systems. By investigating orbital period variations, we found that four of the targets display a secular period decrease while the others exhibit a long-term period increase. The secular period decrease is more likely caused by angular momentum loss while the long-term period increase is due to mass transfer from the less massive component to the more massive one. Based on the statistic of 19 ultrashort period contact binaries with known orbital period changes, we found that seven of them display long-term decrease (three of them also exhibit cyclic variations), ten of them manifest long-term increase while two of them only show cyclic variation and that most of them are shallow contact binaries supporting the long timescale angular momentum loss theory suggested by Stepien. For the three deep contact systems, we found that they are probably triple systems. The tertiary companion plays an essential role during their formation and evolution.

Sensitive radio continuum surveys of the Galactic plane are ideal for discovering new supernova remnants (SNRs). From the Sino-German {\lambda}6 cm polarisation survey of the Galactic plane, an extended shell-like structure has been found at l = 21.8 degree, b = -3.0 degree, which has a size of about 1 degree. New observations were made with the Effelsberg 100-m radio telescope at {\lambda}11 cm to estimate the source spectrum together with the Urumqi {\lambda}6 cm and the Effelsberg {\lambda}21 cm data. The spectral index of G21.8-3.0 was found to be {\alpha} = -0.72 {\pm} 0.16. Polarised emission was mostly detected in the eastern half of G21.8-3.0 at both {\lambda}6 cm and {\lambda}11 cm. These properties, together with the H{\alpha} filament along its northern periphery and the lack of infrared emission, indicate that the emission is non-thermal as is usual in shell-type SNRs.

In this paper we report the results of the first $\sim$four years of spectroscopic and photometric monitoring of the Type IIn supernova SN 2015da (also known as PSN J13522411+3941286, or iPTF16tu). The supernova exploded in the nearby spiral galaxy NGC 5337 in a relatively highly extinguished environment. The transient showed prominent narrow Balmer lines in emission at all times and a slow rise to maximum in all bands. In addition, early observations performed by amateur astronomers give a very well-constrained explosion epoch. The observables are consistent with continuous interaction between the supernova ejecta and a dense and extended H-rich circumstellar medium. The presence of such an extended and dense medium is difficult to reconcile with standard stellar evolution models, since the metallicity at the position of SN 2015da seems to be slightly subsolar. Interaction is likely the mechanism powering the light curve, as confirmed by the analysis of the pseudo bolometric light curve, which gives a total radiated energy $\gtrsim10^{51}\,\rm{erg}$. Modeling the light curve in the context of a supernova shock breakout through a dense circumstellar medium allowed us to infer the mass of the prexisting gas to be $\simeq8\,\rm{M_{\odot}}$, with an extreme mass-loss rate for the progenitor star $\simeq0.6\,\rm{M_{\odot}}\,\rm{yr^{-1}}$, suggesting that most of the circumstellar gas was produced during multiple eruptive events. Near- and mid-infrared observations reveal a flux excess in these domains, similar to those observed in SN 2010jl and other interacting transients, likely due to preexisting radiatively heated dust surrounding the supernova. By modeling the infrared excess, we infer a mass $\gtrsim0.4\times10^{-3}\,\rm{M_{\odot}}$ for the dust.

The expansion of hot electrons in flaring magnetic loops is crucial to understanding the dynamics of solar flares. In this paper we investigate, for the first time, the transport of hot electrons in a magnetic mirror field based on a 1-D particle-in-cell (PIC) simulation. The hot electrons with small pitch angle transport into the cold plasma, which leads to the generation of Langmuir waves in the cold plasma and ion acoustic waves in the hot plasma. The large pitch angle electrons can be confined by the magnetic mirror, resulting in the different evolution time scale between electron parallel and perpendicular temperature. This will cause the formation of electron temperature anisotropy, which can generate the whistler waves near the interface between hot electrons and cold electrons. The whistler waves can scatter the large pitch angle electrons to smaller value through the cyclotron resonance, leading to electrons escaping from the hot region. These results indicate that the whistler waves may play an important role in the transport of electrons in flaring magnetic loops. The findings from this study provide some new insights to understand the electron dynamics of solar flares.

We describe an ultra-wide-bandwidth, low-frequency receiver ("UWL") recently installed on the Parkes radio telescope. The receiver system provides continuous frequency coverage from 704 to 4032 MHz. For much of the band (~60%) the system temperature is approximately 22K and the receiver system remains in a linear regime even in the presence of strong mobile phone transmissions. We discuss the scientific and technical aspects of the new receiver including its astronomical objectives, as well as the feed, receiver, digitiser and signal-processor design. We describe the pipeline routines that form the archive-ready data products and how those data files can be accessed from the archives. The system performance is quantified including the system noise and linearity, beam shape, antenna efficiency, polarisation calibration and timing stability.

We investigate primordial perturbations and non-gaussianities in the Ho\v{r}ava-Lifshitz theory of gravitation. In the UV limit, the scalar perturbation in the Ho\v{r}ava theory is naturally scale-invariant, ignoring the details of the expansion of the Universe. One may thus relax the exponential inflation and the slow-roll conditions for the inflaton field. As a result, it is possible that the primordial non-gaussianities, which are ``slow-roll suppressed'' in the standard scenarios, become large. We calculate the non-gaussianities from the bispectrum of the perturbation and find that the equilateral-type non-gaussianity is of the order of unity, while the local-type non-gaussianity remains small, as in the usual single-field slow-roll inflation model in general relativity. Our result is a new constraint on the the Ho\v{r}ava-Lifshitz gravity.

We investigate the propagation of the gravitational waves in a cosmological background. Based on the framework of spatially covariant gravity, we derive the general quadratic action for the gravitational waves. The spatial derivatives of the extrinsic curvature and the parity-violating terms are systematically introduced. Special attention is paid to the propagation speed of the gravitational waves. We find that it is possible to make the two polarization modes propagate in the same speed, which may differ from that of the light, in the presence of parity-violating terms in the action. In particular, we identify a large class of spatially covariant gravity theories with parity violation, in which both the polarization modes propagate in the speed of light. Our results imply that there are more possibilities in the framework of spatially covariant gravity in light of the propagation speed of the gravitational waves.

By analyzing two sets of complete $BVR_cI_c$ light curves of V342 UMa and three sets of complete $BVR_cI_c$ light curves of V509 Cam, we determined that the two systems are both W-subtype contact binaries and that V342 UMa shows a shallow contact configuration, while V509 Cam exhibits a medium contact configuration. Since both of them are totally eclipsing binaries, the physical parameters derived only by the photometric light curves are reliable.Meanwhile, the period changes of the two targets were analyzed based on all available eclipsing times. We discovered that V342 UMa shows long-term period decrease with a rate of $-1.02(\pm0.54)\times10^{-7}$ days/year and that V509 Cam exhibits long-term period increase with a rate of $3.96(\pm0.90)\times10^{-8}$ days/year. Both the conservative mass transfer and AML via magnetic stellar winds can interpret the long-term period decrease of V342 UMa. The long-term period increase of V509 Cam can be explained by mass transfer from the less massive star to the more massive one. The absolute parameters of the two binaries were estimated according to the Gaia distances and our derived photometric solution results. This method can be extended to other contact binaries without radial velocities but with reliable photometric solutions. The evolutionary states of them were discussed, we found that they reveal identical properties of other W-subtype contact systems.

We report on a multi-observatory campaign to examine asteroid 3200 Phaethon during its December 2017 close approach to Earth, in order to improve our measurements of its fundamental parameters, and to search for surface variations, cometary activity and fragmentation. The mean colors of Phaethon are B-V = 0.702 +/- 0.004, V-R = 0.309 +/- 0.003, R-I = 0.266 +/- 0.004, neutral to slightly blue, consistent with previous classifications of Phaethon as a F-type or B-type asteroid. Variations in Phaethon's B-V colors (but not V-R or R-I) with observer sub-latitude are seen and may be associated with craters observed by the Arecibo radar. High cadence photometry over phases from 20 to 100 degrees allows a fit to the values of the HG photometric parameters; H = 14.57 +/- 0.02, 13.63 +/- 0.02, 13.28 +/- 0.02, 13.07 +/- 0.02; G = 0.00 +/- 0.01, -0.09 +/- 0.01, -0.10 +/- 0.01, -0.08 +/- 0.01 in the BVRI filters respectively; the negative G values are consistent with other observations of F type asteroids. Light curve variations were seen that are also consistent with concavities reported by Arecibo, indicative of large craters on Phaethon's surface whose ejecta may be the source of the Geminid meteoroid stream. A search for gas/dust production set an upper limit of 0.06 +/- 0.02 kg/s when Phaethon was 1.449 AU from the Sun, and 0.2 +/- 0.1 kg/s at 1.067 AU. A search for meter-class fragments accompanying Phaethon did not find any whose on-sky motion was not also consistent with background main belt asteroids.

Main-belt asteroid (6478) Gault was observed to show cometary features in early 2019. To investigate the cause, we conducted {\it BVR} observations at Xingming Observatory, China, from 2019 January to April. The two tails were formed around 2018 October 26--November 08, and 2018 December 29--2019 January 08, respectively, and consisted of dust grains of $\gtrsim$20 $\mu$m to 3 mm in radius ejected at a speed of $0.15 \pm 0.05$ m s$^{-1}$ and following a broken power-law size distribution bending at grain radius $\sim$70 $\mu$m (bulk density 1 g cm$^{-3}$ assumed). The total mass of dust within a $10^4$ km-radius aperture around Gault declined from $\sim$$9 \times 10^6$ kg since 2019 January at a rate of $2.28 \pm 0.07$ kg s$^{-1}$, but temporarily surged around 2019 March 25, because Earth thence crossed the orbital plane of Gault, within which the ejected dust was mainly distributed. No statistically significant colour or short-term lightcurve variation was seen. Nonetheless we argue that Gault is currently subjected to rotational instability. Using the available astrometry, we did not detect any nongravitational acceleration in the orbital motion of Gault.

The period distribution of contact binaries exhibits a very sharp short period cut-off at 0.22 days. In order to provide valuable information on this short period limit, we observed ten totally eclipsing contact binaries with orbital periods near this cut-off. By detailed analysis using W-D code, we determined that two of these systems are A-subtype contact binaries while the others are W-subtype contact binaries and all the targets show shallow contact configurations. Half of the targets exhibit stellar spot activity and four of them have third light. A statistical work on well studied USPCBs was carried out, and physical parameters of fifty-five USPCBs were obtained. Some common properties for these systems were derived. Due to the study of the period-color diagram of USPCBs, we found that the period-color relation of USPCBs is different from other W UMa type contact binaries and USPCBs are metal poor old stars. In addition, the evolutionary states of these systems were discussed by constructing the color-density diagram. We derived that the evolutionary states of the two components of USPCBs show identical characteristics of other contact binaries despite slower evolutionary status caused by smaller mass. We suggested that both the fully convective limit claimed by Rucinski (1992) and the AML theory recommended by Stepien (2006) can produce the short period cut-off and a tertiary companion is very important during the formation of the short period contact binaries by removing angular momentum from the host eclipsing pairs under certain circumstances.

We propose a new class of higher derivative scalar-tensor theories without the Ostrogradsky's ghost instabilities. The construction of our theory is originally motivated by a scalar field with spacelike gradient, which enables us to fix a gauge in which the scalar field appears to be non-dynamical. We dub such a gauge as the spatial gauge. Though the scalar field loses its dynamics, the spatial gauge fixing breaks the time diffeomorphism invariance and thus excites a scalar mode in the gravity sector. We generalize this idea and construct a general class of scalar-tensor theories through a non-dynamical scalar field, which preserves only spatial covariance. We perform a Hamiltonian analysis and confirm that there are at most three (two tensors and one scalar) dynamical degrees of freedom, which ensures the absence of a degree of freedom due to higher derivatives. Our construction opens a new branch of scalar-tensor theories with higher derivatives.

Large-scale radio continuum surveys provide data to get insights into the physical properties of radio sources. HII regions are prominent radio sources produced by thermal emission of ionised gas around young massive stars. We identify and analyse HII regions in the Sino-German 6cm polarisation survey of the Galactic plane. Objects with flat radio continuum spectra together with infrared and/or Halpha emission were identified as HII regions. For HII regions with small apparent sizes, we cross-matched the 6cm small-diameter source catalogue with the radio HII region catalogue compiled by Paladini and the infrared HII region catalogue based on the WISE data. Extended HII regions were identified by eye by overlaying the Paladini and the WISE HII regions onto the 6cm survey images for coincidences. The TT-plot method was employed for spectral index verification. A total of 401 HII regions were identified and their flux densities were determined with the Sino-German 6cm survey data. In the surveyed area, 76 pairs of sources are found to be duplicated in the Paladini HII region catalogue, mainly due to the non-distinction of previous observations with different angular resolutions, and 78 objects in their catalogue are misclassified as HII regions, being actually planetary nebulae, supernova remnants or extragalactic sources that have steep spectra. More than 30 HII regions and HII region candidates from our 6cm survey data, especially extended ones, do not have counterparts in the WISE HII region catalogue, of which 9 are identified for the first time. Based on the newly derived radio continuum spectra and the evidence of infrared emission, the previously identified SNRs G11.1-1.0, G20.4+0.1 and G16.4-0.5 are believed to be HII regions.

Open cluster members are coeval and share the same initial bulk chemical compositions. Consequently, differences in surface abundances between members of a cluster that are at different evolutionary stages can be used to study the effects of mixing and internal chemical processing. We carry out an abundance analysis of seven elements (Li, O, Na, Mg, Al, Si, Fe) in 66 stars belonging to the open cluster M67, based on high resolution GALAH spectra, 1D MARCS model atmospheres, and, for the first time for a large sample of M67 stars, non-local thermodynamic equilibrium (non-LTE) radiative transfer. From the non-LTE analysis, we find a typical star-to-star scatter in the abundance ratios of around 0.05 dex; this scatter is slightly but systematically larger when LTE is assumed instead. We find trends in the abundance ratios with effective temperature, indicating systematic differences in the surface abundances between turn-off and giant stars; these trends are more pronounced when LTE is assumed. However, in the non-LTE analysis, most of the element trends have been flattened. Two species are exceptions to this behaviour, namely Al and Si, which both clearly display remaining trends in the non-LTE analysis. We comment on the possible origin of these trends, by comparing them with recent stellar models that include atomic diffusion.

Context. Pulsars scintillate. Dynamic spectra show brightness variation of pulsars in the time and frequency domain. Secondary spectra demonstrate the distribution of fluctuation power in the dynamic spectra. Aims. Dynamic spectra strongly depend on observational frequencies, but were often observed at frequencies lower than 1.5 GHz. Scintillation observations at higher frequencies help to constrain the turbulence feature of the interstellar medium over a wide frequency range and can detect the scintillations of more distant pulsars. Methods. Ten pulsars were observed at 2250 MHz (S-band) with the Jiamusi 66 m telescope to study their scintillations. Their dynamic spectra were first obtained, from which the decorrelation bandwidths and time scales of diffractive scintillation were then derived by autocorrelation. Secondary spectra were calculated by forming the Fourier power spectra of the dynamic spectra. Results. Most of the newly obtained dynamic spectra are at the highest frequency or have the longest time span of any published data for these pulsars. For PSRs B0540+23, B2324+60 and B2351+61, these were the first dynamic spectra ever reported. The frequencydependence of scintillation parameters indicates that the intervening medium can rarely be ideally turbulent with a Kolmogorov spectrum. The thin screen model worked well at S-band for the scintillation of PSR B1933+16. Parabolic arcs were detected in the secondary spectra of three pulsars, PSRs B0355+54, B0540+23 and B2154+40, all of which were asymmetrically distributed. The inverted arclets of PSR B0355+54 were seen to evolve along the main parabola within a continuous observing session of 12 hours, from which the angular velocity of the pulsar was estimated that was consistent with the measurement by very long baseline interferometry (VLBI).

We investigate the inhomogeneous inflation, in which the space exponentially expands with inhomogeneities, and its cosmological perturbations. The inhomogeneous inflation is realized by introducing scalar fields with spacelike gradients that break the spatial symmetry. We find that the space can expand uniformly in different direction with the same rate. By using the perturbative method, we calculate the corrections to the power spectra of gravitational waves and curvature perturbation up to the linear order in the background inhomogeneities. Since the background is inhomogeneous, perturbations modes with different wave numbers get correlated. We show that generally the power spectra of perturbations depend on the ratio and the angle of wave numbers of the two correlated modes. In particular, the two circular polarization modes of the gravitational waves gain different powers when the background inhomogeneity is of vector or tensor type.

ERIS is an instrument that will both extend and enhance the fundamental diffraction limited imaging and spectroscopy capability for the VLT. It will replace two instruments that are now being maintained beyond their operational lifetimes, combine their functionality on a single focus, provide a new wavefront sensing module that makes use of the facility Adaptive Optics System, and considerably improve their performance. The instrument will be competitive with respect to JWST in several regimes, and has outstanding potential for studies of the Galactic Center, exoplanets, and high redshift galaxies. ERIS had its final design review in 2017, and is expected to be on sky in 2020. This contribution describes the instrument concept, outlines its expected performance, and highlights where it will most excel.

HiPERCAM is a quintuple-beam imager that saw first light on the 4.2m William Herschel Telescope (WHT) in October 2017 and on the 10.4m Gran Telescopio Canarias (GTC) in February 2018. The instrument uses re-imaging optics and 4 dichroic beamsplitters to record ugriz (300-1000nm) images simultaneously on its five CCD cameras. The detectors in HiPERCAM are frame-transfer devices cooled thermo-electrically to -90degC, thereby allowing both long-exposure, deep imaging of faint targets, as well as high-speed (over 1000 windowed frames per second) imaging of rapidly varying targets. In this paper, we report on the as-built design of HiPERCAM, its first-light performance on the GTC, and some of the planned future enhancements.

In this paper, with two-dimensional (2-D) hybrid simulations, we study the generation mechanism of filamentary structures downstream of a quasi-parallel shock. The results show that in the downstream both the amplitude of magnetic field and number density exhibit obvious filamentary structures, and the magnetic field and number density are anticorrelated. Detailed analysis find that these downstream compressive waves propagate almost perpendicular to the magnetic field, and the dominant wave number is around the inverse of ion kinetic scale. Their parallel and perpendicular components roughly satisfies(where and represent the parallel and in-plane perpendicular components of magnetic field, is the wave number in the perpendicular direction, and in the ion gyroradius), and their Alfven ratio also roughly agree with the analytical relation (where and indicate the Alfven ratio and plasma beta), while the corresponding cross helicity and compressibility show good agreement with previous theoretical calculations. All these properties are consistent with those of kinetic slow waves (KSWs). Therefore, we conclude that the filamentary structures in the downstream of a quasi-parallel shock are produced due to the excitation of KSWs.

The Galactic Archaeology with HERMES (GALAH) survey is a large-scale stellar spectroscopic survey of the Milky Way and designed to deliver chemical information complementary to a large number of stars covered by the $Gaia$ mission. We present the GALAH second public data release (GALAH DR2) containing 342,682 stars. For these stars, the GALAH collaboration provides stellar parameters and abundances for up to 23 elements to the community. Here we present the target selection, observation, data reduction and detailed explanation of how the spectra were analysed to estimate stellar parameters and element abundances. For the stellar analysis, we have used a multi-step approach. We use the physics-driven spectrum synthesis of Spectroscopy Made Easy (SME) to derive stellar labels ($T_\mathrm{eff}$, $\log g$, $\mathrm{[Fe/H]}$, $\mathrm{[X/Fe]}$, $v_\mathrm{mic}$, $v \sin i$, $A_{K_S}$) for a representative training set of stars. This information is then propagated to the whole survey with the data-driven method of $The~Cannon$. Special care has been exercised in the spectral synthesis to only consider spectral lines that have reliable atomic input data and are little affected by blending lines. Departures from local thermodynamic equilibrium (LTE) are considered for several key elements, including Li, O, Na, Mg, Al, Si, and Fe, using 1D MARCS stellar atmosphere models. Validation tests including repeat observations, Gaia benchmark stars, open and globular clusters, and K2 asteroseismic targets lend confidence in our methods and results. Combining the GALAH DR2 catalogue with the kinematic information from $Gaia$ will enable a wide range of Galactic Archaeology studies, with unprecedented detail, dimensionality, and scope.

Since its discovery in 2008, the Andromeda galaxy nova M31N 2008-12a has been observed in eruption every single year. This unprecedented frequency indicates an extreme object, with a massive white dwarf and a high accretion rate, which is the most promising candidate for the single-degenerate progenitor of a type-Ia supernova known to date. The previous three eruptions of M31N 2008-12a have displayed remarkably homogeneous multi-wavelength properties: (i) From a faint peak, the optical light curve declined rapidly by two magnitudes in less than two days; (ii) Early spectra showed initial high velocities that slowed down significantly within days and displayed clear He/N lines throughout; (iii) The supersoft X-ray source (SSS) phase of the nova began extremely early, six days after eruption, and only lasted for about two weeks. In contrast, the peculiar 2016 eruption was clearly different. Here we report (i) the considerable delay in the 2016 eruption date, (ii) the significantly shorter SSS phase, and (iii) the brighter optical peak magnitude (with a hitherto unobserved cusp shape). Early theoretical models suggest that these three different effects can be consistently understood as caused by a lower quiescence mass-accretion rate. The corresponding higher ignition mass caused a brighter peak in the free-free emission model. The less-massive accretion disk experienced greater disruption, consequently delaying re-establishment of effective accretion. Without the early refueling, the SSS phase was shortened. Observing the next few eruptions will determine whether the properties of the 2016 outburst make it a genuine outlier in the evolution of M31N 2008-12a.

The hybrid Tibet AS array was successfully constructed in 2014. It has 4500 m$^{2}$ underground water Cherenkov pools used as the muon detector (MD) and 789 scintillator detectors covering 36900 m$^{2}$ as the surface array. At 100 TeV, cosmic-ray background events can be rejected by approximately 99.99\%, according to the full Monte Carlo (MC) simulation for $\gamma$-ray observations. In order to use the muon detector efficiently, we propose to extend the surface array area to 72900 m$^{2}$ by adding 120 scintillator detectors around the current array to increase the effective detection area. A new prototype scintillator detector is developed via optimizing the detector geometry and its optical surface, by selecting the reflective material and adopting dynode readout. This detector can meet our physics requirements with a positional non-uniformity of the output charge within 10\% (with reference to the center of the scintillator), time resolution FWHM of $\sim$2.2 ns, and dynamic range from 1 to 500 minimum ionization particles.

Rapid satellite-to-site visibility determination is of great significance to coverage analysis of satellite constellations as well as onboard mission planning of autonomous spacecraft. This paper presents a novel self-adaptive Hermite interpolation technique for rapid satellite-to-site visibility determination. Piecewise cubic curves are utilized to approximate the waveform of the visibility function versus time. The fourth-order derivative is used to control the approximation error and to optimize the time step for interpolation. The rise and set times are analytically obtained from the roots of cubic polynomials. To further increase the computational speed, an interval shrinking strategy is adopted via investigating the geometric relationship between the ground viewing cone and the orbit trajectory. Simulation results show a 98% decrease in computation time over the brute force method. The method is suitable for all orbital types and analytical orbit propagators.

The Andromeda Galaxy recurrent nova M31N 2008-12a had been observed in eruption ten times, including yearly eruptions from 2008-2014. With a measured recurrence period of $P_\mathrm{rec}=351\pm13$ days (we believe the true value to be half of this) and a white dwarf very close to the Chandrasekhar limit, M31N 2008-12a has become the leading pre-explosion supernova type Ia progenitor candidate. Following multi-wavelength follow-up observations of the 2013 and 2014 eruptions, we initiated a campaign to ensure early detection of the predicted 2015 eruption, which triggered ambitious ground and space-based follow-up programs. In this paper we present the 2015 detection; visible to near-infrared photometry and visible spectroscopy; and ultraviolet and X-ray observations from the Swift observatory. The LCOGT 2m (Hawaii) discovered the 2015 eruption, estimated to have commenced at Aug. $28.28\pm0.12$ UT. The 2013-2015 eruptions are remarkably similar at all wavelengths. New early spectroscopic observations reveal short-lived emission from material with velocities $\sim13000$ km s$^{-1}$, possibly collimated outflows. Photometric and spectroscopic observations of the eruption provide strong evidence supporting a red giant donor. An apparently stochastic variability during the early super-soft X-ray phase was comparable in amplitude and duration to past eruptions, but the 2013 and 2015 eruptions show evidence of a brief flux dip during this phase. The multi-eruption Swift/XRT spectra show tentative evidence of high-ionization emission lines above a high-temperature continuum. Following Henze et al. (2015a), the updated recurrence period based on all known eruptions is $P_\mathrm{rec}=174\pm10$ d, and we expect the next eruption of M31N 2008-12a to occur around mid-Sep. 2016.

In the context of massive (bi-)gravity non-minimal matter couplings have been proposed. These couplings are special in the sense that they are free of the Boulware-Deser ghost below the strong coupling scale and can be used consistently as an effective field theory. Furthermore, they enrich the phenomenology of massive gravity. We consider these couplings in the framework of bimetric gravity and study the cosmological implications for background and linear tensor, vector, and scalar perturbations. Previous works have investigated special branch of solutions. Here we perform a complete perturbation analysis for the general background equations of motion completing previous analysis.

In this work we study the cosmological perturbations in massive bigravity in the presence of non-minimal derivative couplings. For this purpose we consider a specific subclass of Horndeski scalar-tensor interactions that live on the unique composite effective metric. For the viability of the model both metrics have to be dynamical. Nevertheless, the number of allowed kinetic terms is crucial. We adapt to the restriction of having one single kinetic term. After deriving the full set of equations of motion for flat Friedmann-Lemaitre-Robertson-Walker background, we study linear perturbations on top of it. We show explicitly that only four tensor, two vector and two scalar degrees of freedom propagate, one of which being the Horndeski scalar, while the Boulware-Deser ghost can be integrated out.

We explain some phenomena existing generally in the timing residuals: amplitude and sign of the second derivative of a pulsar's spin-frequency ($\ddot\nu$), some sophisticated residual patterns, which also change with the time span of data segments. The sample is taken from Hobbs et al.\,(2010), in which the pulsar's spin-frequency and its first derivative have been subtracted from the timing solution fitting. We first classify the timing residual patterns into different types based on the sign of $\ddot\nu$. Then we use the magnetic field oscillation model developed in our group \citep{zhang12a} to fit successfully the different kinds of timing residuals with the Markov Chain Monte Carlo method. Finally, we simulate the spin evolution over 20 years for a pulsar with typical parameters and analyze the data with the conventional timing solution fitting. By choosing different segments of the simulated data, we find that most of the observed residual patterns can be reproduced successfully. This is the first time that the observed residual patterns are fitted by a model and reproduced by simulations with very few parameters. From the distribution of the different residual patterns in the $P-\dot P$ diagram, we argue that (1) a single magnetic field oscillation mode exists commonly in all pulsars throughout their lifetimes; (2) there may be a transition period over the lifetimes of pulsars, in which multiple magnetic field oscillation modes exist.

In this paper, two-dimensional (2-D) hybrid simulations are performed to investigate ion dynamics at a rippled quasi-parallel shock. The results show that the ripples around the shock front are inherent structures of a quasi-parallel shock, and the reformation of the shock is not synchronous along the surface of the shock front. By following the trajectories of the upstream ions, we find that these ions behave differently when they interact with the shock front at different positions along the shock surface. The upstream particles are easier to transmit through the upper part of a ripple, and the bulk velocity in the corresponding downstream is larger, where a high-speed jet is formed. In the lower part of the ripple, the upstream particles tend to be reflected by the shock. For the reflected ions by the shock, they may suffer multiple stage acceleration when moving along the shock surface, or trapped between the upstream waves and the shock front. At last, these ions may escape to the further upstream or enter the downstream, therefore, the superthermal ions can be found in both the upstream and downstream.

We investigate the cosmological background evolution and perturbations in a general class of spatially covariant theories of gravity, which propagates two tensor modes and one scalar mode. We show that the structure of the theory is preserved under the disformal transformation. We also evaluate the primordial spectra for both the gravitational waves and the curvature perturbation, which are invariant under the disformal transformation. Due to the existence of higher spatial derivatives, the quadratic Lagrangian for the tensor modes itself cannot be transformed to the form in the Einstein frame. Nevertheless, there exists a one-parameter family of frames in which the spectrum of the gravitational waves takes the standard form in the Einstein frame.

XTE J1810-189 underwent an outburst in 2008, and was observed over $\sim 100$ d by RXTE. Performing a time-resolved spectral analysis on the photospheric radius expansion burst detected on 2008 May 4, we obtain the source distance in the range of 3.5--8.7 kpc for the first time. During its outburst, XTE J1810-189 did not enter into the high/soft state, and both the soft and hard colours decreased with decreasing flux. The fractional rms remained at high values ($\sim 30$ per cent). The RXTE/PCA spectra for 3-25 keV can be described by an absorbed power-law component with an additional Gaussian component, and the derived photon index $\Gamma$ increased from $1.84\pm0.01$ to $2.25\pm0.04$ when the unabsorbed X-ray luminosity in 3-25 keV dropped from $4\times10^{36}$ ergs s$^{-1}$ to $6\times10^{35}$ ergs s$^{-1}$. The relatively high flux, dense observations and broadband spectra allow us to provide strong evidence that the softening behaviour detected in the outburst of XTE J1810-189 originates from the evolution of non-thermal component rather than the thermal component (i.e. neutron star surface emission).

We present spectroscopic and photometric data of the Type Ibn supernova (SN) 2014av, discovered by the Xingming Observatory Sky Survey. Stringent pre-discovery detection limits indicate that the object was detected for the first time about 4 days after the explosion. A prompt follow-up campaign arranged by amateur astronomers allowed us to monitor the rising phase (lasting 10.6 days) and to accurately estimate the epoch of the maximum light, on 2014 April 23 (JD = 2456771.1 +/- 1.2). The absolute magnitude of the SN at the maximum light is M(R) = -19.76 +/- 0.16. The post-peak light curve shows an initial fast decline lasting about 3 weeks, and is followed by a slower decline in all bands until the end of the monitoring campaign. The spectra are initially characterized by a hot continuum. Later on, the temperature declines and a number of lines become prominent mostly in emission. In particular, later spectra are dominated by strong and narrow emission features of He I typical of Type Ibn supernovae (SNe), although there is a clear signature of lines from heavier elements (in particular O I, Mg II and Ca II). A forest of relatively narrow Fe II lines is also detected showing P-Cygni profiles, with the absorption component blue-shifted by about 1200 km/s. Another spectral feature often observed in interacting SNe, a strong blue pseudo-continuum, is seen in our latest spectra of SN 2014av. We discuss in this paper the physical parameters of SN 2014av in the context of the Type Ibn supernova variety.

Jupiter-family comet 15P/Finlay has been reportedly quiet in activity for over a century but has harbored two outbursts during its 2014/2015 perihelion passage. Here we present an analysis of these two outbursts using a set of cometary observations. The outbursts took place between 2014 Dec. 15.4--16.0 UT and 2015 Jan. 15.5--16.0 UT as constrained by ground-based and spacecraft observations. We find a characteristic ejection speed of $V_0=300$ to $650 \mathrm{m \cdot s^{-1}}$ for the ejecta of the first outburst and $V_0=550$ to $750 \mathrm{m \cdot s^{-1}}$ for that of the second outburst using a Monte Carlo dust model. The mass of the ejecta is calculated to be $M_\mathrm{d}=2$ to $3\times10^5 \mathrm{kg}$ for the first outburst and $M_\mathrm{d}=4$ to $5\times10^5 \mathrm{kg}$ for the second outburst, corresponds to less than $10^{-7}$ of the nucleus mass. The specific energy of the two outbursts is found to be $0.3$ to $2\times10^5 \mathrm{J \cdot kg^{-1}}$. We also revisit the long-standing puzzle of the non-detection of the hypothetical Finlayid meteor shower by performing a cued search using the 13-year data from the Canadian Meteor Orbit Radar, which does not reveal any positives. The Earth will pass the 2014/2015 outburst ejecta around 2021 Oct. 6 at 22 h UT to Oct. 7 at 1 h UT, with a chance for some significant meteor activity in the radio range, which may provide further clues to the Finlayid puzzle. A southerly radiant in the constellation of Ara will favor the observers in the southern tip of Africa.

We present optical and near infrared observations of the type IIb supernova (SN) 2011fu from a few days to $\sim300$ d after explosion. The SN presents a double-peaked light curve (LC) similar to that of SN 1993J, although more luminous and with a longer cooling phase after the primary peak. The spectral evolution is also similar to SN 1993J's, with hydrogen dominating the spectra to $\sim40$ d, then helium gaining strength, and nebular emission lines appearing from $\sim60$ d post-explosion. The velocities derived from the P-Cygni absorptions are overall similar to those of other type IIb SNe. We have found a strong similarity between the oxygen and magnesium line profiles at late times, which suggests that these lines are forming at the same location within the ejecta. The hydrodynamical modelling of the pseudo-bolometric LC and the observed photospheric velocities suggest that SN 2011fu was the explosion of an extended star ($\rm R\sim450$ R$_\odot$), in which 1.3 $\times 10^{51}$ erg of kinetic energy were released and 0.15 M$_{\rm \odot}$ of $^{56}$Ni were synthesised. In addition, a better reproduction of the observed early pseudo-bolometric LC is achieved if a more massive H-rich-envelope than for other type IIb SNe is considered (0.3 M$_{\rm \odot}$). The hydrodynamical modelling of the LC and the comparison of our late-time spectra with nebular spectral models for type IIb SNe, point to a progenitor for SN 2011fu with a ZAMS mass of 13-18 M$_{\rm \odot}$.

It is speculated that some weakly active comets may be transitional objects between active and dormant comets. These objects are at a unique stage of the evolution of cometary nuclei, as they are still identifiable as active comets, in contrast to inactive comets that are observationally indistinguishable from low albedo asteroids. In this paper, we present a synthesis of comet and meteor observations of Jupiter-family comet 209P/LINEAR, one of the most weakly active comets recorded to-date. Images taken by the Xingming 0.35-m telescope and the Gemini Flamingo-2 camera are modeled by a Monte Carlo dust model, which yields a low dust ejection speed ($1/10$ of that of moderately active comets), dominance of large dust grains, and a low dust production of $0.4~\mathrm{kg \cdot s^{-1}}$ at 19~d after the 2014 perihelion passage. We also find a reddish nucleus of 209P/LINEAR that is similar to D-type asteroids and most Trojan asteroids. Meteor observations with the Canadian Meteor Orbit Radar (CMOR), coupled with meteoroid stream modeling, suggest a low dust production of the parent over the past few hundred orbits, although there are hints of a some temporary increase in activity in the 18th century. Dynamical simulations indicate 209P/LINEAR may have resided in a stable near-Earth orbit for $\sim 10^4$~yr, which is significantly longer than typical JFCs. All these lines of evidence imply that 209P/LINEAR as an aging comet quietly exhausting its remaining near surface volatiles. We also compare 209P/LINEAR to other low activity comets, where evidence for a diversity of the origin of low activity is seen.

We investigate the consequences of general curved trajectories in multi-field inflation. After setting up a completely general formalism using the mass basis, which naturally accommodates the notion of light and heavy modes, we study in detail the simple case of two successive turns in two-field system. We find the power spectrum of the curvature perturbation receives corrections that exhibit oscillatory features sinusoidal in the logarithm of the comoving wavenumber without slow-roll suppression. We show that this is because of the resonance of the heavy modes inside and outside the mass horizon.

The Andromeda Galaxy recurrent nova M31N 2008-12a has been caught in eruption nine times. Six observed eruptions in the seven years from 2008 to 2014 suggested a duty cycle of ~1 year, which makes this the most rapidly recurring system known and the leading single-degenerate Type Ia Supernova progenitor candidate; but no 2010 eruption has been found so far. Here we present evidence supporting the recovery of the 2010 eruption, based on archival images taken at and around the time. We detect the 2010 eruption in a pair of images at 2010 Nov 20.52 UT, with a magnitude of m_R = 17.84 +/- 0.19. The sequence of seven eruptions shows significant indications of a duty cycle slightly shorter than one year, which makes successive eruptions occur progressively earlier in the year. We compared three archival X-ray detections with the well observed multi-wavelength light curve of the 2014 eruption to accurately constrain the time of their optical peaks. The results imply that M31N 2008-12a might have in fact a recurrence period of ~6 months (175 +/- 11 days), making it even more exceptional. If this is the case, then we predict that soon two eruptions per year will be observable. Furthermore, we predict the next eruption will occur around late Sep 2015. We encourage additional observations.

Superbubbles and supershells are the channels for transferring mass and energy from the Galactic disk to the halo. Magnetic fields are believed to play a vital role in their evolution. We study the radio continuum and polarized emission properties of the W4 superbubble to determine its magnetic field strength. New sensitive radio continuum observations were made at 6 cm, 11 cm, and 21 cm. The total intensity measurements were used to derive the radio spectrum of the W4 superbubble. The linear polarization data were analysed to determine the magnetic field properties within the bubble shells. The observations show a multi-shell structure of the W4 superbubble. A flat radio continuum spectrum that stems from optically thin thermal emission is derived from 1.4 GHz to 4.8 GHz. By fitting a passive Faraday screen model and considering the filling factor fne , we obtain the thermal electron density ne = 1.0/\sqrt{fne} (\pm5%) cm^-3 and the strength of the line-of-sight component of the magnetic field B// = -5.0/\sqrt{fne} (\pm10%) {\mu}G (i.e. pointing away from us) within the western shell of the W4 superbubble. When the known tilted geometry of the W4 superbubble is considered, the total magnetic field Btot in its western shell is greater than 12 {\mu}G. The electron density and the magnetic field are lower and weaker in the high-latitude parts of the superbubble. The rotation measure is found to be positive in the eastern shell but negative in the western shell of the W4 superbubble, which is consistent with the case that the magnetic field in the Perseus arm is lifted up from the plane towards high latitudes. The magnetic field strength and the electron density we derived for the W4 superbubble are important parameters for evolution models of superbubbles breaking out of the Galactic plane.

The St\"uckelberg analysis of nonlinear massive gravity in the presence of a general fiducial metric is investigated. We develop a "covariant" formalism for the St\"uckelberg expansion by working with a local inertial frame, through which helicity modes can be characterized correctly. Within this covariant approach, an extended $\Lambda_3$ decoupling limit analysis can be consistently performed, which keeps $\bar{R}_{\mu\nu\rho\sigma}/m^2$ fixed with $\bar{R}_{\mu\nu\rho\sigma}$ the Riemann tensor of the fiducial metric. In this extended decoupling limit, the scalar mode $\pi$ acquires self-interactions due to the presence of the curvature of the fiducial metric. However, the equation of motion for $\pi$ remains of second order in derivatives, which extends the understanding of the absence of the Boulware Deser ghost in the case of a flat fiducial metric.

The simplest possible classical model leading to a cosmological bounce is examined in the light of the non-Gaussianities it can generate. Concentrating solely on the transition between contraction and expansion, and assuming initially purely Gaussian perturbations at the end of the contracting phase, we find that the bounce acts as a source such that the resulting value for the post-bounce $f_{\mathrm{NL}}$ may largely exceed all current limits, to the point of potentially casting doubts on the validity of the perturbative expansion. We conjecture that if one can assume that the non-Gaussianity production depends only on the bouncing behavior of the scale factor and not on the specifics of the model examined, then many realistic models in which a nonsingular classical bounce takes place could exhibit a generic non-Gaussianity excess problem that would need to be addressed for each case.

In the lights of current BICEP2 observations accompanied with the PLANCK satellite results, it has been observed that the simple single field chaotic inflationary models provide a good agreement with their spectral index n_s and large tensor-to-scalar ratio r (0.15 <r <0.26). To explore the other simple models, we consider the fractional-chaotic inflationary potentials of the form V_0 phi^(a/b) where a and b are relatively prime. We show that such kind of inflaton potentials can be realized elegantly in the supergravity framework with generalized shift symmetry and a nature bound a/b < 4 for consistency. Especially, for the number of e-folding from 50 to 60 and some a/b from 2 to 3, our predictions are nicely within at least 1 $\sigma$ region in the r-n_s plane. We also present a systematic investigation of such chaotic inflationary models with fractional exponents to explore the possibilities for the enhancement in the magnitude of running of spectral index (\alpha_{n_s}) beyond the simplistic models.

The origin of the knee in cosmic ray spectrum remains to be an unsolved fundamental problem. There are various kinds of models which predict different break positions and the compositions of the knee. In this work, we suggest to use diffuse $\gamma$-rays and neutrinos as probes to test these models. Based on several typical types of the composition models, the diffuse $\gamma$-ray and neutrino spectra are calculated, which show distinctive cutoff behaviours at energies from tens of TeV to multi-PeV. The expected flux will be observable by the newly upgraded Tibet-AS$\gamma$+MD (muon detector) experiment as well as more sensitive future projects, such as LHAASO and HiSCORE. By comparing the neutrino spectrum with the recent observations by IceCube experiment, we find that the diffuse neutrinos from interactions between the cosmic rays and the interstellar medium may not be responsible to the majority of the IceCube events. Future measurements of the neutrinos may be able to identify the Galactic diffuse component and further shed light on the problem of the knee of cosmic rays.

We successfully embed the Kim-Nilles-Peloso (KNP) alignment mechanism for enhancing the axion decay constant in the context of large volume type IIB orientifolds. The flat direction is generated in the plane of ($C_0$-$C_2$) axions corresponding to the involutively even universal axion $C_0$ and odd axion $C_2$, respectively. The moduli stabilization with large volume scheme has been established as well.

Using $\delta N$ formalism, in the context of a generic multi-field inflation driven on a non-flat field space background, we revisit the analytic expressions of the various cosmological observables such as scalar/tensor power spectra, scalar/tensor spectral tilts, non-Gaussianity parameters, tensor-to-scalar ratio, and the various runnings of these observables. In our backward formalism approach, the subsequent expressions of observables automatically include terms beyond the leading order slow-roll expansion correcting many of the expression at subleading order. To connect our analysis properly with the earlier results, we rederive the (well) known (single field) expressions in the limiting cases of our generic formulae. Further, in the light of PLANCK results, we examine for the compatibility of the consistency relations within the slow-roll regime of a two-field roulette poly-instanton inflation realized in the context of large volume scenarios.

We compute the level of non-gaussianities produced by a cosmological bouncing phase in the minimal non-singular setup that lies within the context of General Relativity when the matter content consists of a simple scalar field with a standard kinetic term. Such a bouncing phase is obtained by requiring that the spatial sections of the background spacetime be positively curved. We restrict attention to the close vicinity of the bounce by Taylor expanding the scale factor, the scalar field and its potential in powers of the conformal time around the bounce. We find that possibly large non-gaussianities are generically produced at the bounce itself and also discuss which shapes of non-gaussianities are mostly likely to be produced.

Large-scale radio continuum surveys have good potential for discovering new Galactic supernova remnants (SNRs). Surveys of the Galactic plane are often limited in the Galactic latitude of |b| ~ 5 degree. SNRs at high latitudes, such as the Cygnus Loop or CTA~1, cannot be detected by surveys in such limited latitudes. Using the available Urumqi 6 cm Galactic plane survey data, together with the maps from the extended ongoing 6 cm medium latitude survey, we wish to discover new SNRs in a large sky area. We searched for shell-like structures and calculated radio spectra using the Urumqi 6 cm, Effelsberg 11 cm, and 21 cm survey data. Radio polarized emission and evidence in other wavelengths are also examined for the characteristics of SNRs. We discover an enclosed oval-shaped object G150.3+4.5 in the 6 cm survey map. It is about 2.5 degree wide and 3 degree high. Parts of the shell structures can be identified well in the 11 cm, 21 cm, and 73.5 cm observations. The Effelsberg 21 cm total intensity image resembles most of the structures of G150.3+4.5 seen at 6 cm, but the loop is not closed in the northwest. High resolution images at 21 cm and 73.5 cm from the Canadian Galactic Plane Survey confirm the extended emission from the eastern and western shells of G150.3+4.5. We calculated the radio continuum spectral indices of the eastern and western shells, which are $\beta \sim -2.4$ and $\beta \sim -2.7$ between 6 cm and 21 cm, respectively. The shell-like structures and their non-thermal nature strongly suggest that G150.3+4.5 is a shell-type SNR. For other objects in the field of view, G151.4+3.0 and G151.2+2.6, we confirm that the shell-like structure G151.4+3.0 very likely has a SNR origin, while the circular-shaped G151.2+2.6 is an HII region with a flat radio spectrum, associated with optical filamentary structure, H$\alpha$, and infrared emission.

We have observed 79 supernova remnants (SNRs) with the Urumqi 25m telescope at 6cm during the Sino-German 6cm polarization survey of the Galactic plane. We measured flux densities of SNRs at 6cm, some of which are the first ever measured or the measurements at the highest frequency, so that we can determine or improve spectra of SNRs. Our observations have ruled out spectral breaks or spectral flattening that were suggested for a few SNRs, and confirmed the spectral break of S147. By combining our 6cm maps with 11cm and 21cm maps from the Effelsberg 100 m telescope, we calculated the spectral index maps of several large SNRs. For many remnants we obtained for the first time polarization images, which show the intrinsic magnetic field structures at 6cm. We disapproved three objects as being SNRs, OA184, G192.8-1.1 and G16.8-1.1, which show a thermal spectrum and no polarization. We have discovered two large supernova remnants, G178.2-4.2 and G25.1-2.3., in the survey maps.

A radiation gene box (RGB) onboard the SJ-10 satellite is a device carrying mice and drosophila cells to determine the biological effects of space radiation environment. The shielded fluxes of different radioactive sources were calculated and the linear energy transfers of gamma-rays, electrons, protons and alpha-particles in tissue were acquired using A-150 tissue-equivalent plastic. Then, a conceptual model of a space radiation instrument employing three semiconductor sub-detectors for deriving the charged and uncharged radiation environment of the RGB was designed. The energy depositions in the three sub-detectors were classified into fifteen channels (bins) in an algorithm derived from the Monte Carlo method. The physical feasibility of the conceptual instrument was also verified by Monte Carlo simulations.

The expansion of interstellar bubbles is suggested to be an important mechanism of triggering material accumulation and star formation. In this work, we study the gaseous environment of a large sample of bubbles, aiming to explore the possible evidence of triggered gas accumulation and star formation in a statistical sense. By cross-matching 6124 Spitzer bubbles from the Milky Way Project (MWP) and more than 2500 Galactic HII regions, we obtain the velocity information for 818 MWP bubbles. To study the gaseous environment of the interstellar bubbles and get rid of the projection effect as much as possible, we constrain the velocity difference between the bubbles and the 13CO(1-0) emission extracted from the Galactic Ring Survey (GRS). Three methods: the mean azimuthally averaged radial profile of 13CO emission, the surface number density of molecular clumps, and the angular cross-correlation function of MWP bubbles and the GRS molecular clumps are adopted. Significant over density of molecular gas is found to be close to the bubble rims. 60 percent of the studied bubbles were found to have associated molecular clumps. By comparing the clump-associated and the clump-unassociated MWP bubbles, we reveal that the bubbles in associations tend to be larger and thicker in physical sizes. We speculate that some of the bubble-associated clumps result from the expansion of bubbles. The fraction of the molecular clumps associated with the MWP bubbles is estimated to be about 20 percent after considering the projection effect. For the bubble-clump complexes, we found that the bubbles in the complexes with associated massive young stellar object(s) (MYSO(s)) have larger physical sizes, hence the complexes tend to be older. We propose that an evolutionary sequence might exist between the relative younger MYSO-unassociated bubble-clump complexes and the MYSO-associated complexes.

Information of small-diameter sources is extracted from the Sino-German 6cm polarisation survey of the Galactic plane carried out with the Urumqi 25-m telescope. We performed two-dimensional elliptical Gaussian fits to the 6cm maps to obtain a list of sources with total-intensity and polarised flux densities. The source list contains 3832 sources with a fitted diameter smaller than 16 arcmin and a peak flux density exceeding 30 mJy, so about 5 times the rms noise, of the total-intensity data. The cumulative source count indicates completeness for flux densities exceeding about 60 mJy. We identify 125 linearly polarised sources at 6cm with a peak polarisation flux density greater than 10 mJy, so about 3 times the rms noise, of the polarised-intensity data. Despite lacking compact steep spectrum sources, the 6cm catalogue lists about 20 percent more sources than the Effelsberg 21cm source catalogue at the same angular resolution and for the same area. Most of the faint 6cm sources must have a flat spectrum and are either HII regions or extragalactic. When compared with the Green Bank 6cm (GB6) catalogue, we obtain higher flux densities for a number of extended sources with complex structures. Polarised 6cm sources density are uniformly distributed in Galactic latitude. Their number density decreases towards the inner Galaxy. More than 80 percent of the polarised sources are most likely extragalactic. With a few exceptions, the sources have a higher percentage polarisation at 6cm than at 21cm. Depolarisation seems to occur mostly within the sources with a minor contribution from the Galactic foreground emission.

Radio continuum observations detect non-thermal synchrotron and thermal bremsstrahlung radiation. Separation of the two different emission components is crucial to study the properties of diffuse interstellar medium. The Cygnus X region is one of the most complex areas in the radio sky which contains a number of massive stars and HII regions on the diffuse thermal and non-thermal background. More supernova remnants are expected to be discovered. We aim to develop a method which can properly separate the non-thermal and thermal radio continuum emission and apply it to the Cygnus X region. The result can be used to study the properties of different emission components and search for new supernova remnants in the complex. Multi-frequency radio continuum data from large-scale surveys are used to develop a new component separation method. Spectral analysis is done pixel by pixel for the non-thermal synchrotron emission with a realistic spectral index distribution and a fixed spectral index of beta = -2.1 for the thermal bremsstrahlung emission. With the new method, we separate the non-thermal and thermal components of the Cygnus X region at an angular resolution of 9.5arcmin. The thermal emission component is found to comprise 75% of the total continuum emission at 6cm. Thermal diffuse emission, rather than the discrete HII regions, is found to be the major contributor to the entire thermal budget. A smooth non-thermal emission background of 100 mK Tb is found. We successfully make the large-extent known supernova remnants and the HII regions embedded in the complex standing out, but no new large SNRs brighter than Sigma_1GHz = 3.7 x 10^-21 W m^-2 Hz^-1 sr^-1 are found.

We investigate non-singular bounce realizations in the framework of ghost-free generalized Galileon cosmology, which furthermore can be free of the anisotropy problem. Considering an Ekpyrotic-like potential we can obtain a total Equation-of-State (EoS) larger than one in the contracting phase, which is necessary for the evolution to be stable against small anisotropic fluctuations. Since such a large EoS forbids the Galileon field to generate the desired form of perturbations, we additionally introduce the curvaton field which can in general produce the observed nearly scale-invariant spectrum. In particular, we provide approximate analytical and exact semi-analytical expressions under which the bouncing scenario is consistent with observations. Finally, the combined Galileon-curvaton system is free of the Big-Rip after the bounce.

We investigate linear cosmological perturbations in multiple field inflationary models where some of the directions are light while others are heavy (with respect to the Hubble parameter). By integrating out the massive degrees of freedom, we determine the multi-dimensional effective theory for the light degrees of freedom and give explicitly the propagation matrix that replaces the effective sound speed of the one-dimensional case. We then examine in detail the consequences of a sudden turn along the inflationary trajectory, in particular the possible breakdown of the low energy effective theory in case the heavy modes are excited. Resorting to a new basis in field space, instead of the usual adiabatic/entropic basis, we study the evolution of the perturbations during the turn. In particular, we compute the power spectrum and compare with the result obtained from the low energy effective theory.

We investigate the coupling relations among perturbations in general multi-field models. We derived the equations of motion for both background and perturbations in a general basis. Within this formalism, we revisit the construction of kinematic orthogonal normal vectors using the successive time derivatives of the background field velocity. We show that the coupling relations among modes in this kinematic basis can be reduced, by employing the background equations of motion for the scalar fields and their high order time derivatives. There are two typical features in the field space: inflationary trajectory and geometry of the potential. Correspondingly, the couplings among modes fall into two categories: one is controlled only by the kinematic quantities, the other involves high order derivatives of the potential. Remarkably, the couplings of the first category, i.e. controlled by the kinematic quantities only, show a "chain" structure. That is, each mode is only coupled to its two neighbour modes.

In the context of two-field inflation characterized by a light direction and a heavy direction, we revisit the question of the impact of the massive modes on the power spectrum produced after a turn in the inflationary trajectory. We consider in particular the resonant effect due to the background oscillations following a sharp turn. Working in the mass basis, i.e. in the basis spanned by the eigenvectors of the effective mass matrix for the perturbations, we provide an analytical estimate of the resonant effect, using the in-in formalism. In comparison with earlier estimates, we find the same the spectral dependence but a smaller amplitude. We also compute, again via the in-in formalism, the effect of the direct coupling between the light and heavy modes at the instant of the turn and confirm our previous results obtained via a different method.

The Sino-German 6cm polarization survey has mapped in total intensity and polarization intensity over an area of approximately 2200 square degrees in the Galactic disk. This survey provides an opportunity to search for Galactic supernova remnants (SNRs) that were previously unknown. We discovered the new SNRs G178.2-4.2 and G25.1-2.3 which have non-thermal spectra, using the 6cm data together with the observations with the Effelsberg telescope at 11 cm and 21 cm. Both G178.2-4.2 and G25.1-2.3 are faint and have an apparent diameter greater than 1deg. G178.2-4.2 shows a polarized shell. HI data suggest that G25.1-2.3 might have a distance of about 3 kpc. The 6cm survey data were also very important to identify two other new SNRs, G152.4-2.1 and G190.9-2.2.

In the context of Type IIB LARGE volume orientifold setup equipped with poly-instanton corrections, the standard single-field poly-instanton inflation driven by a 'Wilson' divisor volume modulus is generalized by the inclusion of respective axion modulus. This two-field dynamics results in a "Roulette" type inflation with the presence of several inflationary trajectories which could produce 50 (or more) e-foldings. The evolution of various trajectories along with physical observables are studied. The possibility of generating primordial non-Gaussianities in the slow-roll as well as in the beyond slow-roll region is investigated. We find that although the non- linearity parameters are quite small during the slow-roll regime, the same are significantly enhanced in the beyond slow-roll regime investigated up to the end of inflation.

We compute the full bispectra, namely both auto- and cross- bispectra, of primordial curvature and tensor perturbations in the most general single-field inflation model whose scalar and gravitational equations of motion are of second order. The formulae in the limits of k-inflation and potential-driven inflation are also given. These expressions are useful for estimating the full bispectra of temperature and polarization anisotropies of the cosmic microwave background radiation.

SCUBA-2 is an innovative 10000 pixel bolometer camera operating at submillimetre wavelengths on the James Clerk Maxwell Telescope (JCMT). The camera has the capability to carry out wide-field surveys to unprecedented depths, addressing key questions relating to the origins of galaxies, stars and planets. With two imaging arrays working simultaneously in the atmospheric windows at 450 and 850 microns, the vast increase in pixel count means that SCUBA-2 maps the sky 100-150 times faster than the previous SCUBA instrument. In this paper we present an overview of the instrument, discuss the physical characteristics of the superconducting detector arrays, outline the observing modes and data acquisition, and present the early performance figures on the telescope. We also showcase the capabilities of the instrument via some early examples of the science SCUBA-2 has already undertaken. In February 2012, SCUBA-2 began a series of unique legacy surveys for the JCMT community. These surveys will take 2.5 years and the results are already providing complementary data to the shorter wavelength, shallower, larger-area surveys from Herschel. The SCUBA-2 surveys will also provide a wealth of information for further study with new facilities such as ALMA, and future telescopes such as CCAT and SPICA.

The Origem Loop in the Galactic anticentre was discovered in 1970s and suggested to be a large supernova remnant. It was argued later to be a chance superposition of unrelated radio sources. We attempt to understand the properties of the Origem Loop. Available multi-frequency radio data were used for the determination of radio spectra of different parts of the Origem Loop and the polarization properties of the loop. Newly available sensitive observations show that the Origem Loop is a loop of more than 6 deg in diameter. It consists of a large non-thermal arc in the north, which we call the Origem Arc, and several known thermal H II regions in the south. Polarized radio emission associated with the arc was detected at 6 cm, revealing tangential magnetic fields. The arc has a brightness temperature spectral index of \beta = -2.70, indicating its non-thermal nature as a supernova remnant. We estimate the distance to the Origem Arc to be about 1.7 kpc, similar to those of some H II regions in the southern part of the loop. The Origem Loop is a visible loop in the sky, which consists of a supernova remnant arc in the north and H II regions in the south.

Light dark matter (DM) with a large DM-nucleon spin-independent cross section and furthermore proper isospin violation (ISV) $f_n/f_p\approx-0.7$ may provide a way to understand the confusing DM direct detection results. Combing with the stringent astrophysical and collider constraints, we systematically investigate the origin of ISV first via general operator analyses and further via specifying three kinds of (single) mediators: A light $Z'$ from chiral $U(1)_X$, an approximate spectator Higgs doublet (It can explain the $W+jj$ anomaly simultaneously) and color triplets. In addition, although $Z'$ from an exotic $U(1)_X$ mixing with $U(1)_Y$ generating $f_n=0$, we can combine it with the conventional Higgs to achieve proper ISV. As a concrete example, we propose the $U(1)_X$ model where the $U(1)_X$ charged light sneutrino is the inelastic DM, which dominantly annihilates to light dark states such as $Z'$ with sub-GeV mass. This model can address the recent GoGeNT annual modulation consistent with other DM direct detection results and free of exclusions.

In this paper we study possible signatures of binary planets or exomoons on the Rossiter-McLaughlin (R-M) effect. Our analyses show that the R-M effect for a binary planet or exomoon during its complete transit phase can be divided into two parts. The first is the conventional one similar to the R-M effect from the transit of a single planet, of which the mass and the projected area are the combinations of the binary components; and the second is caused by the orbital rotation of the binary components, which may add a sine- or linear-mode deviation to the stellar radial velocity curve. We find that the latter effect can be up to several or several ten m/s. By doing numerical simulations as well as analytical analyses, we illustrate that the distribution and dispersion of the latter effects obtained from multiple transit events can be used to constrain the dynamical configuration of the binary planet, such as, how the inner orbit of the binary planet is inclined to its orbit rotating around the central star. We find that the signatures caused by the orbital rotation of the binary components are more likely to be revealed if the two components of binary planet have different masses and mass densities, especially if the heavy one has a high mass density and the light one has a low density. Similar signature on the R-M effect may also be revealed in a hierarchical triple star system containing a dark compact binary and a tertiary star.

We have finished the 6cm polarization survey of the Galactic plane using the Urumqi 25m radio telescope. It covers 10deg<l<230deg in Galactic longitude and |b| <5deg in Galactic latitude. The new polarization maps not only reveal new properties of the diffuse magnetized interstellar medium, but also are very useful for studying individual objects such as Hii regions, which may act as Faraday screens with strong regular magnetic fields inside, and supernova remnants for their polarization properties and spectra. The high sensitivity of the survey enables us to discover two new SNRs G178.2-4.2 and G25.3-2.1 and a number of Hii regions.

We study leading order perturbative corrections to the two point correlation function of the scalar field describing the curvature perturbation in a slow-roll inflationary background, paying particular attention to the contribution of entropy mode loops. We find that the infrared divergences are worse than in pure de Sitter space: they are power law rather than logarithmic. The validity of perturbation theory and thus of the effective field theory of cosmological perturbations leads to stringent constraints on the coupling constants describing the interactions, in our model the quartic self-interaction coupling constant of the entropy field. If the self coupling constant is larger than some critical value which depends in particular on the duration of the inflationary phase, then perturbation theory breaks down. Our analysis may have implications for the stability of de Sitter space: the quantum effects which lead to an instability of de Sitter space will be larger in magnitude in the presence of entropy fluctuations.

We set up cosmological perturbation theory and study the cosmological implications of the so-called ``generalized Galileon'' developed in \cite{Deffayet:2011gz,horndeski}. This is the most general scalar field theory whose Lagrangian contains derivatives up to second order while keeping second order equations of motion, and contains as sub-cases $k$-inflation, $G$-inflation and many other models. We calculate the power spectrum of the primordial curvature perturbation, finding a modification of the usual consistency relation of the tensor-to-scalar ratio in $k$-inflation or perfect fluid models. Finally we also calculate the bispectrum, which contains no new shapes beyond those of $k$-inflation.

We prove the existence of a fully nonlinear conserved curvature perturbation on large scales in Galileon-type scalar field models in two approaches. The first approach is based on the conservation of energy-momentum tensor of the Galileon field, which is also the familiar approach in understanding the conservation in $k$-essence or perfect fluid models. We show that the fluid corresponding to the Galileon field becomes perfect and barotropic on large scales, which is responsible to the conservation. The difference from $k$-essence model is that, besides the energy-momentum conservation, the Einstein equation must be employed in order to complete the proof of barotropy. In the second approach, we derive the fully non-perturbative action for the curvature perturbation $\zeta$ in Galileon model on large scales, and argue that $\zeta=const$ is indeed an exact solution on large scales. This conservation of curvature perturbation is important since it relates the later and the primordial universe.

We study the spectral and polarization properties of supernova remnants (SNRs) based on our 6cm survey data. The observations were taken from the Sino-German 6cm polarization survey of the Galactic plane. By using the integrated flux densities at 6cm together with measurements at other wavelengths from the literature we derive the global spectra of 50 SNRs. In addition, we use the observations at 6cm to present the polarization images of 24 SNRs. We derived integrated flux densities at 6cm for 51 small SNRs with angular sizes less than 1 degree. Global radio spectral indices were obtained in all the cases except for Cas A. For SNRs G15.1-1.6, G16.2-2.7, G16.4-0.5, G17.4-2.3, G17.8-2.6, G20.4+0.1, G36.6+2.6, G43.9+1.6, G53.6-2.2, G55.7+3.4, G59.8+1.2, G68.6-1.2, and G113.0+0.2, the spectra have been significantly improved. From our analysis we argue that the object G16.8-1.1 is probably an HII region instead of a SNR. Cas A shows a secular decrease in total intensity, and we measured a flux density of 688+/-35 Jy at 6cm between 2004 and 2008. Polarized emission from 25 SNRs were detected. For G16.2-2.7, G69.7+1.0, G84.2-0.8 and G85.9-0.6, the polarized emission is detected for the first time confirming them as SNRs. High frequency observations of SNRs are rare but important to establish their spectra and trace them in polarization in particular towards the inner Galaxy where Faraday effects are important.

We completely clarify the feature of primordial non-Gaussianities of tensor perturbations in generalized G-inflation, i.e., the most general single-field inflation model with second order field equations. It is shown that the most general cubic action for the tensor perturbation (gravitational wave) $h_{ij}$ is composed only of two contributions, one with two spacial derivatives and the other with one time derivative on each $h_{ij}$. The former is essentially identical to the cubic term that appears in Einstein gravity and predicts a squeezed shape, while the latter newly appears in the presence of the kinetic coupling to the Einstein tensor and predicts an equilateral shape. Thus, only two shapes appear in the graviton bispectrum of the most general single-field inflation model, which could open a new clue to the identification of inflationary gravitational waves in observations of cosmic microwave background anisotropies as well as direct gravitational wave detection experiments.

We show that modified gravity presents distinctive nonlinear features on the Cosmic Microwave Background (CMB) anisotropies comparing with General Relativity (GR). We calculate the contribution to the CMB non-Gaussianity from nonlinear Sachs-Wolfe effect in $f(R)$ gravity and show that, contrary to GR's contribution which is typically $\lesssim \mathcal{O}(1)$, the contribution in $f(R)$ gravity is sensitive to the nonlinear structure of $f(R)$ and can be large in principle. Optimistically, this gives an alternative origin for the possibly observed large CMB non-Gaussianities besides the primordial ones. On the other hand, such nonlinear features can be employed to provide a new cosmological test of $f(R)$ or other modified theories of gravitation, which is unique and independent of previously known tests.

We investigate the effects of the presence of relic classical superhorizon inhomogeneities during inflation. This superhorizon inhomogeneity appears as a gradient locally and picks out a preferred direction. Quantum fluctuations on this slightly inhomogeneous background are generally statistical anisotropic. We find a quadrupole modification to the ordinary isotropic spectrum. Moreover, this deviation from statistical isotropy is scale-dependent, with a $\sim -1/k^2$ factor. This result implies that the statistical anisotropy mainly appears on large scales, while the spectrum on small scales remains highly isotropic. Moreover, due to this $-1/k^2$ factor, the power on large scales is suppressed. Thus, our model can simultaneously explain the observed anisotropic alignments of the low-$\ell$ multipoles and their low power.

Supernova remnants (SNRs) were often discovered in radio surveys of the Galactic plane. Because of the surface-brightness limit of previous surveys, more faint or confused SNRs await discovery. The Sino-German $\lambda$6\ cm Galactic plane survey is a sensitive survey with the potential to detect new low surface-brightness SNRs. We want to identify new SNRs from the $\lambda$6\ cm survey map of the Galactic plane. We searched for new shell-like objects in the $\lambda$6\ cm survey maps, and studied their radio emission, polarization, and spectra using the $\lambda$6\ cm maps together with the $\lambda$11\ cm and $\lambda$21\ cm Effelsberg observations. Extended polarized objects with non-thermal spectra were identified as SNRs. We have discovered two new, large, faint SNRs, G178.2-4.2 and G25.1-2.3, both of which show shell structure. G178.2-4.2 has a size of 72 arcmin x 62 arcmin with strongly polarized emission being detected along its northern shell. The spectrum of G178.2-4.2 is non-thermal, with an integrated spectral index of $\alpha = -0.48\pm0.13$. Its surface brightness is $\Sigma_{1 GHz} = 7.2 x 10^{-23}{Wm^{-2} Hz^{-1} sr^{-1}}$, which makes G178.2-4.2 the second faintest known Galactic SNR. G25.1-2.3 is revealed by its strong southern shell which has a size of 80 arcmin x 30\arcmin. It has a non-thermal radio spectrum with a spectral index of $\alpha = -0.49\pm0.13$. Two new large shell-type SNRs have been detected at $\lambda$6\ cm in an area of 2200 deg^2 along the the Galactic plane. This demonstrates that more large and faint SNRs exist, but are very difficult to detect.

The external Compton (EC) model is used to study the high energy emission of some blazars, in which the external photon field is considered to dominate inverse Compton radiation. We explore the properties of external photon field through analyzing the FERMI LAT bright AGN sample within three months detection. In the sample, assuming the high energy radiation of low synchrotron peaked blazars from the EC process, we find that the external photon parameter Uext/\nuext may not be a constant. Calculating synchrotron and inverse Compton (IC) luminosity from the quasi-simultaneous broadband spectrum energy distributions (SEDs), we find that they have an approximately linear relation. This indicates that the ratio of external photon and magnetic energy density is a constant in the comoving frame, implying that the Lorentz factor of the emitting blob depends on external photon field and magnetic field. The result gives a strong constraint on the jet dynamic model.

Observations of large supernova remnants (SNRs) at high frequencies are rare, but provide valuable information about their physical properties. The total intensity and polarization properties of 16 large SNRs in the Galactic plane were investigated based on observations of the Urumqi $\lambda$6\ cm polarization survey of the Galactic plane with an angular resolution of 9$\farcm$5. We extracted total intensity and linear polarization maps of large SNRs from the Urumqi $\lambda$6\ cm survey, obtained their integrated flux densities, and derived the radio spectra in context with previously published flux densities at various frequencies. In particular, Effelsberg $\lambda$11\ cm and $\lambda$21\ cm survey data were used for calculating integrated flux densities. The $\lambda$6\ cm polarization data also delineate the magnetic field structures of the SNRs. We present the first total intensity maps at $\lambda$6\ cm for SNRs G106.3+2.7, G114.3+0.3, G116.5+1.1, G166.0+4.3 (VRO 42.05.01), G205.5+0.5 (Monoceros Nebula) and G206.9+2.3 (PKS 0646+06) and the first polarization measurements at $\lambda$6\ cm for SNRs G82.2+5.3 (W63), G106.3+2.7, G114.3+0.3, G116.5+1.1, G166.0+4.3 (VRO 42.05.01), G205.5+0.5 (Monoceros Nebula) and G206.9+2.3 (PKS 0646+06). Most of the newly derived integrated radio spectra are consistent with previous results. The new flux densities obtained from the Urumqi $\lambda$6\ cm, Effelsberg $\lambda$11\ cm and $\lambda$21\ cm surveys are crucial to determine the spectra of SNR G65.1+0.6, G69.0+2.7 (CTB 80), G93.7-0.2 and G114.3+0.3. We find that G192.8$-$1.1 (PKS 0607+17) consists of background sources, \ion{H}{II} regions and the extended diffuse emission of thermal nature, and conclude that G192.8$-$1.1 is not a SNR.

We perform a detailed phase-space analysis of a class of k-essence cosmology. We find the critical points can be divided into three classes: points unstable but the model stable, both points and the model stable, points stable but the model unstable. Like the case of points unstable but the model stable, the case of points stable but the model unstable is not relevant from a cosmological point of view, though they can be late-time attractors for the universe. So in order to study the possible final state of the universe, it is important to investigate not only the stability of the critical points but also the stability of the model. The case of both points and the model stable can further be divided into two classes: points only presenting decelerating phases and points at which all decelerating, constant-speed, and accelerating phases can appear; the final state of the universe dependents on the potential.

In this letter, we introduce a new method of image stacking to directly study the undetected but possible gamma-ray point sources. Applying the method to the Australia Telescope 20 GHz Survey (AT20G) sources which have not been detected by LAT on Fermi, we find that the sources contribute (10.5+/-1.1)% and (4.3+/-0.9)% of the extragalactic gamma-ray background (EGB) and have a very soft spectrum with the photon indexes of 3.09+/-0.23 and 2.61+/-0.26, in the 1-3 and 3-300GeV energy ranges. In the 0.1-1GeV range, they probably contribute more large faction to the EGB, but it is not quite sure. It maybe not appropriate to assume that the undetected sources have the similar property to the detected sources.

We make an complementary investigation of the primordial trispectrum from exchanging intermediate scalar modes in multi-field inflation models with generalized kinetic terms. Together with the calculation of irreducible contributions to the primordial trispectrum in Ref.[103], we give the full leading-order primordial trispectrum in generalized multi-field models.

Non-linear CMB temperature anisotropies up to the third-order on large scales are calculated. On large scales and in the Sachs-Wolfe limit, we give the explicit expression for the observed temperature anisotropy in terms of the primordial curvature perturbation up to the third-order. We derived the final bispectrum and trispectrum of anisotropies and the corresponding non-linear parameters, in which the contributions to the observed non-Gaussianity from primordial perturbations and from the non-linear mapping from primordial curvature perturbation to the temperature anisotropy are transparently separated.

Linearly polarized Galactic synchrotron emission provides valuable information about the properties of the Galactic magnetic field and the interstellar magneto-ionic medium, when Faraday rotation along the line of sight is properly taken into account. We aim to survey the Galactic plane at 6 cm including linear polarization. At such a short wavelength Faraday rotation effects are in general small and the Galactic magnetic field properties can be probed to larger distances than at long wavelengths. The Urumqi 25-m telescope is used for a sensitive 6 cm survey in total and polarized intensities. WMAP K-band (22.8 GHz) polarization data are used to restore the absolute zero-level of the Urumqi U and Q maps by extrapolation. Total intensity and polarization maps are presented for a Galactic plane region of 129 degree < l < 230 degree and |b| < 5 degree in the anti-centre with an angular resolution of 9'5 and an average sensitivity of 0.6 mK and 0.4 mK Tb in total and polarized intensity, respectively. We briefly discuss the properties of some extended Faraday Screens detected in the 6 cm polarization maps. The Sino-German 6 cm polarization survey provides new information about the properties of the magnetic ISM. The survey also adds valuable information for discrete Galactic objects and is in particular suited to detect extended Faraday Screens with large rotation measures hosting strong regular magnetic fields.

We investigate the effects of couplings between curvature and isocurvature perturbations before and around horizon-crossings during cosmological inflation. We consider a generalized two-field inflation model, in which the non-canonical kinetic term allows us arbitrary sound speeds of curvature and isocurvature perturbations. By using the field-theoretical perturbative analysis, we calculate the cross-spectrum between curvature and isocurvature perturbations and the corrections to curvature and isocurvature power spectra due to the presence of couplings between them. Our analysis confirms previous results that the cross-correlations are generated and amplified when perturbations cross the horizons. Moreover, we find the cross-correlation, which was previously shown to be first-order in slow-roll parameter, can be enhanced when the sound speed of isocurvature perturbation is much smaller than that of the curvature perturbation. This is because in this case the isocurvature perturbation exits its horizon much earlier than the curvature perturbation and acts as a nearly constant source on the curvature perturbation.

We analyze the large-scale two-dimensional sidereal anisotropy of multi-TeV cosmic rays by Tibet Air Shower Array, with the data taken from 1999 November to 2008 December. To explore temporal variations of the anisotropy, the data set is divided into nine intervals, each in a time span of about one year. The sidereal anisotropy of magnitude about 0.1% appears fairly stable from year to year over the entire observation period of nine years. This indicates that the anisotropy of TeV Galactic cosmic rays remains insensitive to solar activities since the observation period covers more than a half of the 23rd solar cycle.

Using the Tibet-III air shower array, we search for TeV gamma-rays from 27 potential Galactic sources in the early list of bright sources obtained by the Fermi Large Area Telescope at energies above 100 MeV. Among them, we observe 7 sources instead of the expected 0.61 sources at a significance of 2 sigma or more excess. The chance probability from Poisson statistics would be estimated to be 3.8 x 10^-6. If the excess distribution observed by the Tibet-III array has a density gradient toward the Galactic plane, the expected number of sources may be enhanced in chance association. Then, the chance probability rises slightly, to 1.2 x 10^-5, based on a simple Monte Carlo simulation. These low chance probabilities clearly show that the Fermi bright Galactic sources have statistically significant correlations with TeV gamma-ray excesses. We also find that all 7 sources are associated with pulsars, and 6 of them are coincident with sources detected by the Milagro experiment at a significance of 3 sigma or more at the representative energy of 35 TeV. The significance maps observed by the Tibet-III air shower array around the Fermi sources, which are coincident with the Milagro >=3sigma sources, are consistent with the Milagro observations. This is the first result of the northern sky survey of the Fermi bright Galactic sources in the TeV region.

Ho\v{r}ava-Lifshitz gravity is a potentially UV complete theory with important implications for the very early universe. In particular, in the presence of spatial curvature it is possible to obtain a non-singular bouncing cosmology. The bounce is realized as a consequence of higher order spatial curvature terms in the gravitational action. Here, we extend the study of linear cosmological perturbations in Ho\v{r}ava-Lifshitz gravity coupled to matter in the case when spatial curvature is present. As in the case without spatial curvature, we find that there is no extra dynamical degree of freedom for scalar metric perturbations. We study the evolution of fluctuations through the bounce and show that the solutions remain non-singular throughout. If we start with quantum vacuum fluctuations on sub-Hubble scales in the contracting phase, and if the contracting phase is dominated by pressure-less matter, then for $\lambda = 1$ and in the infrared limit the perturbations at late times are scale invariant. Thus, Ho\v{r}ava-Lifshitz gravity can provide a realization of the ``matter bounce'' scenario of structure formation.

We investigate the primordial non-Gaussianities from the trispectra in multi-field inflation models, which can be seen as generalization of multi-field $k$-inflation and multi-DBI inflation. We derive the full fourth-order perturbation action for the inflaton fields and evaluate the four-point correlation functions for the perturbations in the limit $\ca \ll 1$ and $\ce \ll1$. There are three types of momentum-dependent shape functions which arise from three types of four-point interaction vertices. The final trispectrum of the curvature perturbation can be expressed in terms of the deformations and permutations of these three shape functions, and is determined by $\ca$, $\ce$, $\lambda$, $\Pi$ which depend on the non-linear structure of the model and also the transfer function $T_{\Rc\Sc}$. We also discuss the parameter space for the trispectrum and plot the shape diagrams for the trispectrum both for visualization and for distinguishing different shapes from each other.

We develop a model anisotropy best-fitting to the two-dimensional sky-map of multi-TeV galactic cosmic ray (GCR) intensity observed with the Tibet III air shower (AS) array. By incorporating a pair of intensity excesses in the hydrogen deflection plane (HDP) suggested by Gurnett et al., together with the uni-directional and bi-directional flows for reproducing the observed global feature, this model successfully reproduces the observed sky-map including the "skewed" feature of the excess intensity from the heliotail direction, whose physical origin has long remained unknown. These additional excesses are modeled by a pair of the northern and southern Gaussian distributions, each placed ~50 degree away from the heliotail direction. The amplitude of the southern excess is as large as ~0.2 %, more than twice the amplitude of the northern excess. This implies that the Tibet AS experiment discovered for the first time a clear evidence of the significant modulation of GCR intensity in the heliotail and the asymmetric heliosphere.

We investigate the primordial trispectra of the general multifield DBI inflationary model. In contrast with the single field model, the entropic modes can source the curvature perturbations on the super horizon scales, so we calculate the contributions from the interaction of four entropic modes mediating one adiabatic mode to the trispectra, at the large transfer limit ($T_{RS}\gg1$). We obtained the general form of the 4-point correlation functions, plotted the shape diagrams in two specific momenta configurations, "equilateral configuration" and "specialized configuration". Our figures showed that we can easily distinguish the two different momenta configurations.

We investigate one-loop quantum corrections to the power spectrum of adiabatic perturbation from entropy modes/adiabatic mode cross-interactions in multiple DBI inflationary models. We find that due to the non-canonical kinetic term in DBI models, the loop corrections are enhanced by slow-varying parameter $\epsilon$ and small sound speed $c_s$. Thus, in general the loop-corrections in multi-DBI models can be large. Moreover, we find that the loop-corrections from adiabatic/entropy cross-interaction vertices are IR finite.

The spectral energy distribution (SED) of TeV blazars has a double-humped shape that is usually interpreted as Synchrotron Self Compton (SSC) model. The one zone SSC model is used broadly but cannot fit the high energy tail of SED very well. It need bulk Lorentz factor which is conflict with the observation. Furthermore one zone SSC model can not explain the entire spectrum. In the paper, we propose a new model that the high energy emission is produced by the accelerated protons in the blob with a small size and high magnetic field, the low energy radiation comes from the electrons in the expanded blob. Because the high and low energy photons are not produced at the same time, the requirement of large Doppler factor from pair production is relaxed. We present the fitting results of the SEDs for Mkn 501 during April 1997 and Mkn 421 during March 2001 respectively.

The synchrotron self-Compton (SSC) models and External Compton (EC) models of AGN jets with continually longitudinal and transverse bulk velocity structures are constructed. The observed spectra show complex and interesting patterns in different velocity structures and viewing angles. These models are used to calculate the synchrotron and inverse Compton spectra of two typical BL Lac objects (BLO) (Mrk 421 and 0716+714) and one Flat Spectrum Radio Quasars (FSRQs) (3c 279), and to discuss the implications of jet bulk velocity structures in unification of the BLO and FR I radio galaxies (FRI). By calculating the synchrotron spectra and SSC spectra of BL Lac object jets with continually bulk velocity structures, we find that the spectra are much different from ones in jets with uniform velocity structure under the increase of viewing angles. The unification of BLO and FRI is less constrained by viewing angles and would be imprinted by velocity structures intrinsic to the jet themselves. By considering the jets with bulk velocity structures constrained by apparent speed, we discuss the velocity structures imprinted on the observed spectra for different viewing angles. We find that the spectra are greatly impacted by longitudinal velocity structures, becasue the volume elements are compressed or expanded. Finally, we present the EC spectra of FSRQs and FR II radio galaxies (FRII) and find that they are weakly affected by velocity structures compared to synchrotron and SSC spectra.

We study cosmological perturbations in Ho\v{r}ava-Lifshitz Gravity. We consider scalar metric fluctuations about a homogeneous and isotropic space-time. Starting from the most general metric, we work out the complete second order action for the perturbations. We then make use of the residual gauge invariance and of the constraint equations to reduce the number of dynamical degrees of freedom. After introducing the Sasaki-Mukhanov variable, the combination of spatial metric fluctuation and matter inhomogeneity for which the action in General Relativity has canonical form, we find that this variable has the standard time derivative term in the second order action, and that the extra degree of freedom is non-dynamical. The limit $\lambda \to 1$ is well-behaved, unlike what is obtained when performing incomplete analyses of cosmological fluctuations. Thus, there is no strong coupling problem for Ho\v{r}ava-Lifshitz gravity when considering cosmological solutions. We also compute the spectrum of cosmological perturbations. If the potential in the action is taken to be of "detailed balance" form, we find a cancelation of the highest derivative terms in the action for the curvature fluctuations. As a consequence, the initial spectrum of perturbations will not be scale-invariant in a general spacetime background. As an application, we consider fluctuations in an inflationary background and draw connections with the "trans-Planckian problem" for cosmological perturbations. In the special case in which the potential term in the action is of detailed balance form and in which $\lambda = 1$, the equation of motion for cosmological perturbations in the far UV takes the same form as in GR. However, in general the equation of motion is characterized by a modified dispersion relation.

We perform a general study of the primordial scalar non-Gaussianities in multi-field inflationary models in Einstein gravity. We consider models governed by a Lagrangian which is a general function of the scalar fields and their first spacetime derivatives. We use $\delta N$ formalism to relate scalar fields and curvature perturbations. We calculate the explicit cubic order perturbation action and the three-point function of curvature perturbation evaluated at horizon-crossing. Under reasonable assumptions, in the limit of small slow-varying parameters and a sound speed $c_s$ close to one, we find that the non-Gaussianity is completely determined by these slow-varying parameters and some other parameters determined by the structure of the inflationary models. Our work generalizes previous results, and would be useful to study non-Gaussianity in multi-field inflationary models that will be constructed in the future.

Evidence for a spectral flattening of the supernova remnant (SNR) HB3 (G132.7+1.3) was recently claimed in literature based on previously published total flux density data, and the flattening was further interpreted as the discovery of thermal bremsstrahlung emission in the shell of HB3. A spectral flattening has never been observed from any SNR before. Reliable observations of HB3 at frequencies above 3000 MHz are crucial to confirm such a spectral behaviour. We extracted 4800 MHz total intensity and polarisation data of HB3 from the Sino-German 6 cm polarisation survey of the Galactic plane made with the Urumqi 25 m telescope, and analysed the spectrum of HB3, together with Effelsberg data at 1408 MHz and 2675 MHz. We found an overall spectral index of HB3 of alpha=-0.61+-0.06 between 1408 MHz and 4800 MHz, similar to the index at lower frequencies. There is no spectral flattening at high frequencies. We detected strong polarised emission from HB3 at 4800 MHz. Our 4800 MHz data show a tangential field orientation in the HB3 shell.

Using the data taken from Tibet II High Density (HD) Array (1997 February-1999 September) and Tibet-III array (1999 November-2005 November), our previous northern sky survey for TeV $\gamma-$ray point sources has now been updated by a factor of 2.8 improved statistics. From $0.0^{\circ}$ to $60.0^{\circ}$ in declination (Dec) range, no new TeV $\gamma-$ray point sources with sufficiently high significance were identified while the well-known Crab Nebula and Mrk421 remain to be the brightest TeV $\gamma-$ray sources within the field of view of the Tibet air shower array. Based on the currently available data and at the 90% confidence level (C.L.), the flux upper limits for different power law index assumption are re-derived, which are approximately improved by 1.7 times as compared with our previous reported limits.

The amplitude of the Compton-Getting (CG) anisotropy contains the power-law index of the cosmic-ray energy spectrum. Based on this relation and using the Tibet air-shower array data, we measure the cosmic-ray spectral index to be $-3.03 \pm 0.55_{stat} \pm < 0.62_{syst}$ between 6 TeV and 40 TeV, consistent with $-$2.7 from direct energy spectrum measurements. Potentially, this CG anisotropy analysis can be utilized to confirm the astrophysical origin of the ``knee'' against models for non-standard hadronic interactions in the atmosphere.

The Tibet air shower array, which has an effective area of 37,000 square meters and is located at 4300 m in altitude, has been observing air showers induced by cosmic rays with energies above a few TeV. We are planning to add a large muon detector array to it for the purpose of increasing its sensitivity to cosmic gamma rays in the 100 TeV (10 - 1000 TeV) energy region by discriminating them from cosmic-ray hadrons. We report on the possibility of detection of gamma rays in the 100 TeV energy region in our field of view, based on the improved sensitivity of our air shower array deduced from the full Monte Carlo simulation.

The Tibet air shower array, which has an effective area of 37,000 square meters and is located at 4300 m in altitude, has been observing air showers induced by cosmic rays with energies above a few TeV. We have a plan to add a large muon detector array to it for the purpose of increasing its sensitivity to cosmic gamma rays in the 100 TeV energy region by discriminating them from cosmic-ray hadrons. We have deduced the attainable sensitivity of the muon detector array using our Monte Carlo simulation. We report here on the detailed procedure of our Monte Carlo simulation.

We report the discovery of five SiC grains and one graphite grain isolated from the Murchison carbonaceous meteorite whose major-element isotopic compositions indicate an origin in nova explosions. The grains are characterized by low 12C/13C (4-9) and 14N/15N (5-20) ratios, large excesses in 30Si (30Si/28Si ratios range to 2.1 times solar) and high 26Al/27Al ratios. These isotopic signatures are theoretically predicted for the ejecta from ONe novae and cannot be matched by any other stellar sources. Previous studies of presolar grains from primitive meteorites have shown that the vast majority formed in red giant outflows and supernova ejecta. Although a classical nova origin was suggested for a few presolar graphite grains on the basis of 22Ne enrichments, this identification is somewhat ambiguous since it is based only on one trace element. Our present study presents the first evidence for nova grains on the basis of major element isotopic compositions of single grains. We also present the results of nucleosynthetic calculations of classical nova models and compare the predicted isotopic ratios with those of the grains. The comparison points toward massive ONe novae if the ejecta are mixed with material of close-to-solar composition.