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

The peculiar velocities of biased tracers of the cosmic density field contain important information about the growth of large scale structure and generate anisotropy in the observed clustering of galaxies. Using N-body data, we show that velocity expansions for halo redshift-space power spectra are converged at the percent-level at perturbative scales for most line-of-sight angles $\mu$ when the first three pairwise velocity moments are included, and that the third moment is well-approximated by a counterterm-like contribution. We compute these pairwise-velocity statistics in Fourier space using both Eulerian and Lagrangian one-loop perturbation theory using a cubic bias scheme and a complete set of counterterms and stochastic contributions. We compare the models and show that our models fit both real-space velocity statistics and redshift-space power spectra for both halos and a mock sample of galaxies at sub-percent level on perturbative scales using consistent sets of parameters, making them appealing choices for the upcoming era of spectroscopic, peculiar-velocity and kSZ surveys.

We study the ability of perturbative models with effective field theory contributions and infra-red resummation to model the redshift space clustering of biased tracers in models where the linear power spectrum has "features" -- either imprinted during inflation or induced by non-standard expansion histories. We show that both Eulerian and Lagrangian perturbation theory are capable of reproducing the Fourier space two-point functions of halos up to the non-linear scale from a suite of $4096^3$ particle N-body simulations. This is the first demonstration that perturbative models can accurately fit the redshift-space clustering of biased tracers in N-body simulations of such theories. By comparing different theoretical models and IR resummation schemes we assess the current theoretical uncertainty in predicting power spectra for models with features. Our results suggest that future surveys will be able to detect or tightly constrain features in the primordial spectrum below the one percent level across a wide range of scales.

We present the results of new radio observations carried out with the Karl G. Jansky Very Large Array C-configuration at 5.5 GHz for a sample of southern narrow-line Seyfert 1 galaxies (NLS1s). This work increases the number of known radio-detected NLS1s in the southern hemisphere, and confirms that the radio emission of NLS1s is mainly concentrated in a central region at kpc-scale and only a few sources show diffuse emission. In RQ NLS1s, the radio luminosity tends to be higher in steep-spectrum sources and be lower in flat-spectrum sources, which is opposite to RL NLS1s. This may be because the radio emission of steep NLS1s is dominated by misaligned jets, AGN-driven outflows, or star formation superposing on a compact core. Instead the radio emission of flat NLS1s may be produced by a central core which has not yet developed radio jets and outflows. We discover new NLS1s harboring kpc-scale radio jets and confirm that a powerful jet does not require a large-mass black hole to be generated. We also find sources dominated by star formation. These NLS1s could be new candidates in investigating the radio emission of different mechanisms.

Investigating in the laboratory the process of matter accretion onto forming stars through scaled experiments is important in order to better understand star and planetary systems formation and evolution. Such experiments can indeed complement observations by providing access to the processes with spatial and temporal resolution. A first step has been made in [G. Revet et al., Science Advances 3, e1700982 (2017), arXiv:1708.02528]. in allowing such investigations. It revealed the existence of a two components stream: a hot shell surrounding a cooler inner stream. The shell was formed by matter laterally ejected upon impact and refocused by the local magnetic field. That laboratory investigation was limited tonormal incidence impacts. However, in young stellar objects, complex structure of magnetic fields causes variability of the incidentangles of the accretion columns. This led us to undertake an investigation, using laboratory plasmas, of the consequence of having a slanted accretion impacting a young star. Here we use high power laser interactions and strong magnetic field generation in the laboratory, complemented by numerical simulations, to study the asymmetry induced upon accretion structures when columns of matter impact the surface of young stars with an oblique angle. Compared to the scenario where matter accretes normal to the star surface, we observe strongly asymmetric plasma structure, strong lateral ejecta of matter, poor confinement of the accreted material and reduced heating compared to the normal incidence case. Thus, slanted accretion is a configuration that seems to be capable of inducing perturbations of the chromosphere and hence possibly influence the level of activity of the corona.

The shaping of astrophysical outflows into bright, dense and collimated jets due to magnetic pressure is here investigated using laboratory experiments. We notably look at the impact on jet collimation of a misalignment between the outflow, as it stems from the source, and the magnetic field. For small misalignments, a magnetic nozzle forms and redirects the outflow in a collimated jet. For growing misalignments, this nozzle becomes increasingly asymmetric, disrupting jet formation. Our results thus suggest outflow/magnetic field misalignment to be a plausible key process regulating jet collimation in a variety of objects. Furthermore, they provide a possible interpretation for the observed structuring of astrophysical jets. Jet modulation could be interpreted as the signature of changes over time in the outflow/ambient field angle, and the change in the direction of the jet could be the signature of changes in the direction of the ambient field.

Transition Radiation (TR) plays an important role in particle identification in high-energy physics and its characteristics provide a feasible method of energy calibration in the energy range up to 10 TeV, which is of interest for dark matter searches in cosmic rays. In a Transition Radiation Detector (TRD), the TR signal is superimposed onto the ionization energy loss signal induced by incident charged particles. In order to make the TR signal stand out from the background of ionization energy loss in a significant way, we optimized both the radiators and the detector. We have designed a new prototype of regular radiator optimized for a maximal TR photon yield, combined with the Side-On TRD which is supposed to improve the detection efficiency of TR. We started a test beam experiment with the Side-On TRD at Conseil Europ\'{e}en pour la Recherche Nucl\'{e}aire (CERN), and found that the experimental data is consistent with the simulation results.

We identify a strong Ly-alpha damping wing profile in the spectrum of the quasar P183+05 at z=6.4386. Given the detection of several narrow metal absorption lines at z=6.40392, the most likely explanation for the absorption profile is that it is due to a damped Ly-alpha system. However, in order to match the data a contribution of an intergalactic medium 5-38% neutral or additional weaker absorbers near the quasar is also required. The absorption system presented here is the most distant damped Ly-alpha system currently known. We estimate an HI column density ($10^{20.68\pm0.25}\,$cm$^{-2}$), metallicity ([O/H]$=-2.92\pm 0.32$), and relative chemical abundances of a system consistent with a low-mass galaxy during the first Gyr of the universe. This object is among the most metal-poor damped Ly-alpha systems known and, even though it is observed only ~850 Myr after the big bang, its relative abundances do not show signatures of chemical enrichment by Population III stars.

We explore the deep ultraviolet (that is, short-distance) limit of the power spectrum (PS) and of the correlation function of a cold dark matter dominated Universe. While for large scales the PS can be written as a double series expansion, in powers of the linear PS and of the wavenumber $k$, we show that, in the opposite limit, it can be expressed via an expansion in powers of the form $1/k^{d+2n}$, where $d$ is the number of spatial dimensions, and $n$ is a non negative integer. The coefficients of the terms of the expansion are nonperturbative in the linear PS, and can be interpreted in terms of the probability density function for the displacement field, evaluated around specific configurations of the latter, that we identify. In the case of the Zel'dovich dynamics, these coefficients can be determined analytically, whereas for the exact dynamics they can be treated as fit, or nuisance, parameters. We confirm our findings with numerical simulations and discuss the necessary steps to match our results to those obtained for larger scales and to actual measurements.

Density-reconstruction sharpens the baryon acoustic oscillations signal by undoing some of the smoothing incurred by nonlinear structure formation. In this paper we present an analytical model for reconstruction based on the Zeldovich approximation, which for the first time includes a complete set of counterterms and bias terms up to quadratic order and can fit real and redshift-space data pre- and post-reconstruction data in both Fourier and configuration space over a wide range of scales. We compare our model to n-body data at $z = 0$ from the {\tt DarkSky} simulation \cite{skillman14}, finding sub-percent agreement in both real space and in the redshift-space power spectrum monopole out to $k = 0.4\ h$ Mpc$^{-1}$, and out to $k = 0.2\ h$ Mpc$^{-1}$ in the quadrupole, with comparable agreement in configuration space. We compare our model with several popular existing alternatives, updating existing theoretical results for exponential damping in wiggle/no-wiggle splits of the BAO signal and discuss the usually-ignored effect of higher bias contributions on the reconstructed signal. In the appendices, we re-derive the former within our formalism, present exploratory results on higher-order corrections due to nonlinearities inherent to reconstruction, and present numerical techniques with which to calculate the redshift-space power spectrum of biased tracers within the Zeldovich approximation.

We investigate the range of applicability of a model for the real-space power spectrum based on N-body dynamics and a (quadratic) Lagrangian bias expansion. This combination uses the highly accurate particle displacements that can be efficiently achieved by modern N-body methods with a symmetries-based bias expansion which describes the clustering of any tracer on large scales. We show that at low redshifts, and for moderately biased tracers, the substitution of N-body-determined dynamics improves over an equivalent model using perturbation theory by more than a factor of two in scale, while at high redshifts and for highly biased tracers the gains are more modest. This hybrid approach lends itself well to emulation. By removing the need to identify halos and subhalos, and by not requiring any galaxy-formation-related parameters to be included, the emulation task is significantly simplified at the cost of modeling a more limited range in scale.

We report on the observations of gamma-ray burst (GRB) 190114C by the Fermi Gamma-ray Space Telescope and the Neil Gehrels Swift Observatory. The early-time observations reveal multiple emission components that evolve independently, with a delayed power-law component that exhibits significant spectral attenuation above 40 MeV in the first few seconds of the burst. This power-law component transitions to a harder spectrum that is consistent with the afterglow emission observed at later times. This afterglow component is clearly identifiable in the GBM and BAT light curves as a slowly fading emission component on which the rest of the prompt emission is superimposed. As a result, we are able to constrain the transition from internal shock to external shock dominated emission. We find that the temporal and spectral evolution of the broadband afterglow emission can be well modeled as synchrotron emission from a forward shock propagating into a wind-like circumstellar environment and find that high-energy photons observed by Fermi LAT are in tension with the theoretical maximum energy that can be achieved through synchrotron emission from a shock. These violations of the maximum synchrotron energy are further compounded by the detection of very high energy (VHE) emission above 300 GeV by MAGIC concurrent with our observations. We conclude that the observations of VHE photons from GRB 190114C necessitates either an additional emission mechanism at very high energies that is hidden in the synchrotron component in the LAT energy range, an acceleration mechanism that imparts energy to the particles at a rate that is faster than the electron synchrotron energy loss rate, or revisions of the fundamental assumptions used in estimating the maximum photon energy attainable through the synchrotron process.

We report a search for astronomical neutrinos in the energy region from several GeV to TeV in the direction of the blazar TXS0506+056 using the Super-Kamiokande detector following the detection of a 100 TeV neutrino from the same location by the IceCube collaboration. Using Super-Kamiokande neutrino data across several data samples observed from April 1996 to February 2018 we have searched for both a total excess above known backgrounds across the entire period as well as localized excesses on smaller time scales in that interval. No significant excess nor significant variation in the observed event rate are found in the blazar direction. Upper limits are placed on the electron and muon neutrino fluxes at 90\% confidence level as $6.03 \times 10^{-7}$ and $4.52 \times 10^{-7}$ to $9.26 \times 10^{-10}$ [${\rm erg}/{\rm cm}^2/{\rm s}$], respectively.

Narrow-line Seyfert 1 galaxies (NLS1s) is one of the few classes of active galactic nuclei (AGN) harboring powerful relativistic jets and detected in $\gamma$ rays. NLS1s are well-known X-ray sources. While in non-jetted sources the origin of this X-ray emission may be a hot corona surrounding the accretion disk, in jetted objects, especially beamed ones, the contribution of corona and relativistic jet is difficult to disentangle without a proper sampling of the hard X-ray emission. For this reason, we observed with \textit{NuSTAR} the first four NLS1s detected at high energy $\gamma$ rays. These data, along with \textit{XMM-Newton} and \textit{Swift/XRT} observations, confirmed that X rays originate both in the jet and in the accretion disk corona. Time variability in hard X rays furthermore suggests that, as observed in flat-spectrum radio quasars, the dissipation region during flares could change its position from source to source, and it can be located both inside and outside the broad-line region. We find that jetted NLS1s, and other blazars as well, seem not to follow the classical fundamental plane of black hole activity, which therefore should be used as a black hole mass estimator in blazars with extreme care only. Our results strengthen the idea according to which $\gamma$-NLS1s are smaller and younger version of flat-spectrum radio quasars, in which both a Seyfert and a blazar component co-exist.

GRB 131108A is a bright long Gamma-Ray Burst (GRB) detected by the Large Area Telescope and the Gamma-ray Burst Monitor on board the \textit{Fermi Gamma-ray Space Telescope}. Dedicated temporal and spectral analyses reveal three $\gamma$-ray flares dominating above 100 MeV, which are not directly related to the prompt emission in the GBM band (10 keV--10 MeV). The high-energy light curve of GRB 131108A (100 MeV -- 10 GeV) shows an unusual evolution: a steep decay, followed by three flares with an underlying emission, and then a long-lasting decay phase. The detailed analysis of the $\gamma$-ray flares finds that the three flares are 6 -- 20 times brighter than the underlying emission and are similar to each other. The fluence of each flare, (1.6 $\sim$ 2.0) $\times$ 10$^{-6}$ erg cm$^{-2}$, is comparable to that of emission during the steep decay phase, 1.7 $\times$ 10$^{-6}$ erg cm$^{-2}$. The total fluence from three $\gamma$-ray flares is 5.3 $\times$ 10$^{-6}$ erg cm$^{-2}$. The three $\gamma$-ray flares show properties similar to the usual X-ray flares that are sharp flux increases, occurring in $\sim$ 50\% of afterglows, in some cases well after the prompt emission. Also, the temporal and spectral indices during the early steep decay phase and the decaying phase of each flare show the consistency with a relation of the curvature effect ($\hat{\alpha}$ = 2 + $\hat{\beta}$), which is the first observational evidence of the high-latitude emission in the GeV energy band.

The primary advantage of moderately superheated bubble chamber detectors is their simultaneous sensitivity to nuclear recoils from WIMP dark matter and insensitivity to electron recoil backgrounds. A comprehensive analysis of PICO gamma calibration data demonstrates for the first time that electron recoils in C$_3$F$_8$ scale in accordance with a new nucleation mechanism, rather than one driven by a hot-spike as previously supposed. Using this semi-empirical model, bubble chamber nucleation thresholds may be tuned to be sensitive to lower energy nuclear recoils while maintaining excellent electron recoil rejection. The PICO-40L detector will exploit this model to achieve thermodynamic thresholds as low as 2.8 keV while being dominated by single-scatter events from coherent elastic neutrino-nucleus scattering of solar neutrinos. In one year of operation, PICO-40L can improve existing leading limits from PICO on spin-dependent WIMP-proton coupling by nearly an order of magnitude for WIMP masses greater than 3 GeV c$^{-2}$ and will have the ability to surpass all existing non-xenon bounds on spin-independent WIMP-nucleon coupling for WIMP masses from 3 to 40 GeV c$^{-2}$.

Precise measurement about the cosmic-ray (CR) component knee is essential for revealing the mistery of CR's acceleration and propagation mechanism, as well as exploring the new physics. However, classification about the CR components is a tough task especially for the groups with the atomic number close to each other. Realizing that the deep learning has achieved a remarkable breakthrough in many fields, we seek for leveraging this technology to improve the classification performance about the CR Proton and Light groups on the LHAASO-KM2A experiment. In this work, we propose a fused Graph Neural Network model in combination of the KM2A arrays, in which the activated detectors are structured into graphs. We find that the signal and background can be effectively discriminated in this model, and its performance outperforms both the traditional physics-based method and the CNN-based model across the whole energy range.

We present a detailed study of 11 narrow-line Seyfert 1 galaxies (NLS1s) from the Six-degree Field Galaxy Survey (6dFGS) that both have optical and X-ray spectroscopic observations. There are five complex NLS1s (C-NLS1s) and six simple NLS1s (S-NLS1s). We propose a possible correlation between [O III] line asymmetry and X-ray complexity. The outflow or wind from the inner accretion disk is commonly present in NLS1s and mostly directed along the system axis. In C-NLS1s only weak wind effects are measured, the X-ray spectral complexity might be caused by the presence of ionized material in the wind. On the contrary, the wind in S-NLS1s is fast, the ionized material could be swept by such a strong wind, thus the complex feature is missing which results in a simple X-ray spectrum. Furthermore, this outflow scenario seems to be an inclination effect. Since the speed of the wind is higher in a small inclination while lower in a large inclination, S-NLS1s might be sources viewed at small angles while C-NLS1s might be sources viewed at large angles.

Supernova detection is a major objective of the Super-Kamiokande (SK) experiment. In the next stage of SK (SK-Gd), gadolinium (Gd) sulfate will be added to the detector, which will improve the ability of the detector to identify neutrons. A core-collapse supernova will be preceded by an increasing flux of neutrinos and anti-neutrinos, from thermal and weak nuclear processes in the star, over a timescale of hours; some of which may be detected at SK-Gd. This could provide an early warning of an imminent core-collapse supernova, hours earlier than the detection of the neutrinos from core collapse. Electron anti-neutrino detection will rely on inverse beta decay events below the usual analysis energy threshold of SK, so Gd loading is vital to reduce backgrounds while maximising detection efficiency. Assuming normal neutrino mass ordering, more than 200 events could be detected in the final 12 hours before core collapse for a 15-25 solar mass star at around 200 pc, which is representative of the nearest red supergiant to Earth, $\mathrm{\alpha}$Ori (Betelgeuse). At a statistical false alarm rate of 1 per century, detection could be up to 10 hours before core collapse, and a pre-supernova star could be detected by SK-Gd up to 600 pc away. A pre-supernova alert could be provided to the astrophysics community following gadolinium loading.

In this paper, we describe the International Pulsar Timing Array second data release, which includes recent pulsar timing data obtained by three regional consortia: the European Pulsar Timing Array, the North American Nanohertz Observatory for Gravitational Waves, and the Parkes Pulsar Timing Array. We analyse and where possible combine high-precision timing data for 65 millisecond pulsars which are regularly observed by these groups. A basic noise analysis, including the processes which are both correlated and uncorrelated in time, provides noise models and timing ephemerides for the pulsars. We find that the timing precisions of pulsars are generally improved compared to the previous data release, mainly due to the addition of new data in the combination. The main purpose of this work is to create the most up-to-date IPTA data release. These data are publicly available for searches for low-frequency gravitational waves and other pulsar science.

Magnetic reconnection occurs when two plasmas having co-planar but anti-parallel magnetic fields meet. At the contact point, the field is locally annihilated and the magnetic energy can be released into the surrounding plasma. Theory and numerical modelling still face many challenges in handling this complex process, the predictability of which remains elusive. Here we test, through a laboratory experiment conducted in a controlled geometry, the effect of changing the field topology from two-dimensional to three-dimensional. This is done by imposing an out-of-plane (guide) magnetic field of adjustable strength. A strong slowing down or even halting of symmetric reconnection is observed, even for a weak guide-field. Concomitantly, we observe a delayed heating of the plasma in the reconnection region and modified particle acceleration, with super-Alfvenic outflows ejected along the reconnection layer. These observations highlight the importance of taking into account three-dimensional effects in the many reconnection events taking place in natural and laboratory environments.

We present extensive datasets for a class of intermediate-luminosity optical transients known as "luminous red novae" (LRNe). They show double-peaked light curves, with an initial rapid luminosity rise to a blue peak (at -13 to -15 mag), which is followed by a longer-duration red peak that sometimes is attenuated, resembling a plateau. The progenitors of three of them (NGC4490-2011OT1, M101-2015OT1, and SNhunt248), likely relatively massive blue to yellow stars, were also observed in a pre-eruptive stage when their luminosity was slowly increasing. Early spectra obtained during the first peak show a blue continuum with superposed prominent narrow Balmer lines, with P Cygni profiles. Lines of Fe II are also clearly observed, mostly in emission. During the second peak, the spectral continuum becomes much redder, Halpha is barely detected, and a forest of narrow metal lines is observed in absorption. Very late-time spectra (~6 months after blue peak) show an extremely red spectral continuum, peaking in the infrared (IR) domain. Halpha is detected in pure emission at such late phases, along with broad absorption bands due to molecular overtones (such as TiO, VO). We discuss a few alternative scenarios for LRNe. Although major instabilities of single massive stars cannot be definitely ruled out, we favour a common envelope ejection in a close binary system, with possibly a final coalescence of the two stars. The similarity between LRNe and the outburst observed a few months before the explosion of the Type IIn SN 2011ht is also discussed.

We study the extreme-ultraviolet (EUV) emissions modulated by leaky fast sausage modes (FSMs) in solar active region loops and examine their observational signatures via spectrometers like EIS. After computing fluid variables of leaky FSMs with MHD simulations, we forward-model the intensity and spectral properties of the Fe X 185~\AA~and Fe XII 195~\AA~lines by incorporating non-equilibrium ionization (NEI) in the computations of the relevant ionic fractions. The damping times derived from the intensity variations are then compared with the wave values, namely the damping times directly found from our MHD simulations. Our results show that in the equilibrium ionization cases, the density variations and the intensity variations can be either in phase or in anti-phase, depending on the loop temperature. NEI considerably impacts the intensity variations but has only marginal effects on the derived Doppler velocity or Doppler width. We find that the damping time derived from the intensity can largely reflect the wave damping time if the loop temperature is not drastically different from the nominal formation temperature of the corresponding emission line. These results are helpful for understanding the modulations to the EUV emissions by leaky FSMs and hence helpful for identifying FSMs in solar active region loops.

This white paper envisions a revolutionary post-DESI, post-LSST dark energy program based on intensity mapping of the redshifted 21cm emission line from neutral hydrogen at radio frequencies. The proposed intensity mapping survey has the unique capability to quadruple the volume of the Universe surveyed by optical programs, provide a percent-level measurement of the expansion history to $z \sim 6$, open a window to explore physics beyond the concordance $\Lambda$CDM model, and to significantly improve the precision on standard cosmological parameters. In addition, characterization of dark energy and new physics will be powerfully enhanced by cross-correlations with optical surveys and cosmic microwave background measurements. The rich dataset obtained by the proposed intensity mapping instrument will be simultaneously useful in exploring the time-domain physics of fast radio transients and pulsars, potentially in live "multi-messenger" coincidence with other observatories. The core dark energy/inflation science advances enabled by this program are the following: (i) Measure the expansion history of the universe over $z=0.3-6$ with a single instrument, extending the range deep into the pre-acceleration era, providing an unexplored window for new physics; (ii) Measure the growth rate of structure in the universe over the same redshift range; (iii) Observe, or constrain, the presence of inflationary relics in the primordial power spectrum, improving existing constraints by an order of magnitude; (iv) Observe, or constrain, primordial non-Gaussianity with unprecedented precision, improving constraints on several key numbers by an order of magnitude. Detailed mapping of the enormous, and still largely unexplored, volume of cosmic space will thus provide unprecedented information on fundamental questions of the vacuum energy and early-universe physics.

We present an analytic model to describe the supermassive black hole binary (SMBHB) merger rate in the Universe with astrophysical observables: galaxy stellar mass function, pair fraction, merger timescale and black hole - host galaxy relations. We construct observational priors and compute the allowed range of the characteristic spectrum $h_c$ of the gravitational wave background (GWB) to be $10^{-16}<h_c<10^{-15}$ at a frequency of $f=1/{\rm yr}$. We exploit our parametrization to tackle the problem of astrophysical inference from Pulsar Timing Array (PTA) observations. We simulate a series of upper limits and detections and use a nested sampling algorithm to explore the parameter space. Corroborating previous results, we find that the current PTA non-detection does not place significant constraints on any observables; however, either future upper limits or detections will significantly enhance our knowledge of the SMBHB population. If a GWB is not detected at a level of $h_c(f=1/{\rm yr})=10^{-17}$, our current understanding of galaxy and SMBHB mergers is disfavoured at a $5\sigma$ level, indicating a combination of severe binary stalling, over-estimating of the SMBH -- host galaxy relations, and extreme dynamical properties of merging SMBHBs. Conversely, future detections of a Square Kilometre Array (SKA)-type array will allow to constrain the normalization of the SMBHB merger rate in the Universe, the time between galaxy pairing and SMBHB merging, the normalization of the SMBH -- host galaxy relations and the dynamical binary properties, including their eccentricity and density of stellar environment.

We present high-cadence ultraviolet (UV), optical, and near-infrared (NIR) data on the luminous Type II-P supernova SN 2017gmr from hours after discovery through the first 180 days. SN 2017gmr does not show signs of narrow, high-ionization emission lines in the early optical spectra, yet the optical lightcurve evolution suggests that an extra energy source from circumstellar medium (CSM) interaction must be present for at least 2 days after explosion. Modeling of the early lightcurve indicates a ~500R$_{\odot}$ progenitor radius, consistent with a rather compact red supergiant, and late-time luminosities indicate up to 0.130 $\pm$ 0.026 M$_{\odot}$ of $^{56}$Ni are present, if the lightcurve is solely powered by radioactive decay, although the $^{56}$Ni mass may be lower if CSM interaction contributes to the post-plateau luminosity. Prominent multi-peaked emission lines of H$\alpha$ and [O I] emerge after day 154, as a result of either an asymmetric explosion or asymmetries in the CSM. The lack of narrow lines within the first two days of explosion in the likely presence of CSM interaction may be an example of close, dense, asymmetric CSM that is quickly enveloped by the spherical supernova ejecta.

We study synergies between three promising methods to measure $2<z<5$ large-scale structure in the next decade. Optical spectroscopic surveys are the most mature, but become increasingly difficult at $z>2$ and suffer from interloper problems even for spectroscopic surveys. Intensity mapping of the 21-cm signal can cover large volumes with exquisite fidelity, but is limited both by loss of information to foreground cleaning and by lack of knowledge of the mean signal. Cosmic microwave background (CMB) lensing is theoretically very clean, but ultimately measures just the projected variations in density. We find that cross-correlation between optical and radio can significantly improve the measurement of growth rate. Combining these with the CMB provides a promising avenue to detecting modified gravity at high redshifts, in particular by independently probing the Weyl and Newtonian potentials and by strengthening control of systematics. We find that cross-correlating a Stage {\sc ii} 21-cm survey with DESI quasars with a reasonable brightness temperature prior could enable measurements of the growth rate $f\sigma_8$ at sub 3\% and sub 8\% levels at $z = 3, 4$, representing a factor of 4 and 8 improvement over constraints obtainable from DESI quasars alone. Similarly, cross-correlating 21-cm data with a futuristic LBG survey to $m_{UV}<24.5$ over 1000 square degrees will make possible $f\sigma_8$ measurements at close to 1\% at $z = 3$ and 3\% at $z = 4$, and improve similar constraints at $z = 5$ by close to a factor of 3 to sub-10\% precision. Combining the above with CMB lensing from a Stage 4 CMB survey and LSST data can additionally constrain the gravitational slip $\gamma$ parameter to similar precision at these redshifts, enabling us to test the predictions of general relativity at large scales.

We explore Lagrangian perturbation theory (LPT) for biased tracers in the presence of two fluids, focusing on the case of cold dark matter (CDM) and baryons. The presence of two fluids induces corrections to the Lagrangian bias expansion and tracer advection, both of which we formulate as expansions in the three linear modes of the Lagrangian equations of motion. We compute the linear-order two-fluid corrections in the Zeldovich approximation, finding that modifications to the bias expansion and tracer advection both enter as percent-level corrections over a large range of wavenumbers at low redshift and draw parallels with the Eulerian formalism. We then discuss nonlinear corrections in the two-fluid picture, and calculate contributions from the relative velocity effect ($\propto v_r^2$) at one loop order. Finally, we conduct an exploratory Fisher analysis to assess the impact of two-fluid corrections on baryon acoustic oscillations (BAO) measurements, finding that while modest values of the relative bias parameters can introduce systematic biases in the measured BAO scale of up to $0.5\, \sigma$, fitting for these effects as additional parameters increases the error bar by less than $30\%$ across a wide range of bias values.

The Large High Altitude Air Shower Observatory (LHAASO) project is a new generation multi-component instrument, to be built at 4410 meters of altitude in the Sichuan province of China, with the aim to study with unprecedented sensitivity the spec trum, the composition and the anisotropy of cosmic rays in the energy range between 10$^{12}$ and 10$^{18}$ eV, as well as to act simultaneously as a wide aperture (one stereoradiant), continuously-operated gamma ray telescope in the energy range between 10$^{11}$ and $10^{15}$ eV. The experiment will be able of continuously surveying the TeV sky for steady and transient sources from 100 GeV to 1 PeV, t hus opening for the first time the 100-1000 TeV range to the direct observations of the high energy cosmic ray sources. In addition, the different observables (electronic, muonic and Cherenkov/fluorescence components) that will be measured in LHAASO will allow to investigate origin, acceleration and propagation of the radiation through a measurement of energy spec trum, elemental composition and anisotropy with unprecedented resolution. The remarkable sensitivity of LHAASO in cosmic rays physics and gamma astronomy would play a key-role in the comprehensive general program to explore the High Energy Universe. LHAASO will allow important studies of fundamental physics (such as indirect dark matter search, Lorentz invariance violation, quantum gravity) and solar and heliospheric physics. In this document we introduce the concept of LHAASO and the main science goals, providing an overview of the project.

The gamma-ray sky has been observed with unprecedented accuracy in the last decade by the Fermi large area telescope (LAT), allowing us to resolve and understand the high-energy Universe. The nature of the remaining unresolved emission (unresolved gamma-ray background, UGRB) below the LAT source detection threshold can be uncovered by characterizing the amplitude and angular scale of the UGRB fluctuation field. This work presents a measurement of the UGRB autocorrelation angular power spectrum based on eight years of Fermi LAT Pass 8 data products. The analysis is designed to be robust against contamination from resolved sources and noise systematics. The sensitivity to subthreshold sources is greatly enhanced with respect to previous measurements. We find evidence (with $\sim$3.7$\sigma$ significance) that the scenario in which two classes of sources contribute to the UGRB signal is favored over a single class. A double power law with exponential cutoff can explain the anisotropy energy spectrum well, with photon indices of the two populations being 2.55 $\pm$ 0.23 and 1.86 $\pm$ 0.15.

Magnetic field are transported and tangled by turbulence, even as they lose identity due to nonideal or resistive effects. On balance field lines undergo stretch-twist-fold processes. The curvature field, a scalar that measures the tangling of the magnetic field lines, is studied in detail here, in the context of magnetohydrodynamic turbulence. A central finding is that the magnitudes of the curvature and the magnetic field are anti-correlated. High curvature co-locates with low magnetic field, which gives rise to power-law tails of the probability density function of the curvature field. The curvature drift term that converts magnetic energy into flow and thermal energy, largely depends on the curvature field behavior, a relationship that helps to explain particle acceleration due to curvature drift. This adds as well to evidence that turbulent effects most likely play an essential role in particle energization since turbulence drives stronger tangled field configurations, and therefore curvature.

The expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $z < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. The large cosmological volume at $2 < z < 5$, together with the ability to efficiently target high-$z$ galaxies with known techniques, enables large gains in the study of Inflation and Dark Energy. A future spectroscopic survey can test the Gaussianity of the initial conditions up to a factor of ~50 better than our current bounds, crossing the crucial theoretical threshold of $\sigma(f_{NL}^{\rm local})$ of order unity that separates single field and multi-field models. Simultaneously, it can measure the fraction of Dark Energy at the percent level up to $z = 5$, thus serving as an unprecedented test of the standard model and opening up a tremendous discovery space.

Inspired by recent IRIS observations, we forward model the response of the Fe XXI 1354 \AA\ line to fundamental, standing, linear fast sausage modes (FSMs) in flare loops. Starting with the fluid parameters for an FSM in a straight tube with equilibrium parameters largely compatible with the IRIS measurements, we synthesize the line profiles by incorporating the non-Equilibrium Ionization (NEI) effect in the computation of the contribution function. We find that both the intensity and Doppler shift oscillate at the wave period ($P$). The phase difference between the two differs from the expected value ($90^\circ$) only slightly because NEI plays only a marginal role in determining the ionic fraction of Fe XXI in the examined dense loop. The Doppler width modulations, however, posses an asymmetry in the first and second halves of a wave period, leading to a secondary periodicity at $P/2$ in addition to the primary one at $P$. This behavior results from the competition between the broadening due to bulk flow and that due to temperature variations, with the latter being stronger but not overwhelmingly so. These expected signatures, with the exception of the Doppler width, are largely consistent with the IRIS measurements, thereby corroborating the reported detection of a fundamental FSM. The forward modeled signatures are useful for identifying fundamental FSMs in flare loops from measurements of the Fe \MyRoman{21} 1354 \AA\ line with instruments similar to IRIS, even though a much higher cadence is required for the expected behavior in the Doppler widths to be detected.

The hadronic interaction of cosmic rays with solar atmosphere can produce high energy gamma rays. The gamma-ray luminosity is correlated both with the flux of primary cosmic rays and the intensity of the solar magnetic field. The gamma rays below 200 GeV have been observed by $Fermi$ without any evident energy cutoff. The bright gamma-ray flux above 100 GeV has been detected only during solar minimum. The only available data in TeV range come from the HAWC observations, however outside the solar minimum. The ARGO-YBJ dataset has been used to search for sub-TeV/TeV gamma rays from the Sun during the solar minimum from 2008 to 2010, the same time period covered by the Fermi data. A suitable model containing the Sun shadow, solar disk emission and inverse-Compton emission has been developed, and the chi-square minimization method was used to quantitatively estimate the disk gamma-ray signal. The result shows that no significant gamma-ray signal is detected and upper limits to the gamma-ray flux at 0.3$-$7 TeV are set at 95\% confidence level. In the low energy range these limits are consistent with the extrapolation of the Fermi-LAT measurements taken during solar minimum and are compatible with a softening of the gamma-ray spectrum below 1 TeV. They provide also an experimental upper bound to any solar disk emission at TeV energies. Models of dark matter annihilation via long-lived mediators predicting gamma-ray fluxes > $10^{-7}$ GeV $cm^{-2}$ $s^{-1}$ below 1 TeV are ruled out by the ARGO-YBJ limits.

Young stars exhibit short-term photometric variability caused by mass accretion events from circumstellar disks, the presence of dusty warps within the inner disks, starspots that rotate across the stellar surfaces, and flares. Long-term variability also occurs owing to starspot longevity and cycles, and from changes in stellar angular momenta and activity as the stars age. We propose to observe the Carina star-forming region in different bands with a cadence of 30 minutes every night for one week per year to clarify the nature of both the short-term and long-term variability of the thousands of young stars in this region. By obtaining well-sampled multicolor lightcurves of this dense young cluster, LSST would acquire the first statistically significant data on how these objects vary on both short and long timescales. This information will allow us to relate the observed variability to stellar properties such as mass, age, binarity, and to environmental properties such as location within or exterior to the H II region, and to the presence or absence of a circumstellar disk.

Narrow-line Seyfert 1 galaxies (NLS1s) are a class of active galactic nuclei (AGN) which are known to be one of the few sources of $\gamma$ rays, which originate in a relativistic beamed jet. Becuase of their relatively large distance, a poorly investigated aspect of these jetted NLS1s is their environment, and in particular their host galaxy. In this work we present the results of a morphological analysis of the host galaxy of the jetted NLS1 IRAS 20181-2244 observed with the 6.5m Baade Telescope of the Las Campanas Observatory. The GALFIT analysis ran on the Ks image, along with additional spectroscopic observations performed with the Nordic Optical Telescope, clearly revealed the presence of an interacting system of two galaxies. The data suggest that this NLS1 is hosted by a late-type galaxy, although the result is not conclusive. This analysis, along with other results in the literature, might suggest that two populations of jetted NLS1 exist. Further morphological studies are needed to confirm or disprove this hypothesis.

Forward-modeling the emission properties in various passbands is important for confidently identifying magnetohydrodynamic waves in the structured solar corona. We examine how Non-equilibrium Ionization (NEI) affects the Extreme Ultraviolet (EUV) emissions modulated by standing fast sausage modes (FSMs) in coronal loops, taking the Fe IX 171 \AA\ and Fe XII 193 \AA\ emission lines as examples. Starting with the expressions for linear FSMs in straight cylinders, we synthesize the specific intensities and spectral profiles for the two spectral lines by incorporating the self-consistently derived ionic fractions in the relevant contribution functions. We find that relative to the case where Equilibrium Ionization (EI) is assumed, NEI considerably impacts the intensity modulations, but shows essentially no effect on the Doppler velocities or widths. Furthermore, NEI may affect the phase difference between intensity variations and those in Doppler widths for Fe XII 193 \AA\ when the line-of-sight is oblique to the loop axis. While this difference is $180^\circ$ when EI is assumed, it is $\sim 90^\circ$ when NEI is incorporated for the parameters we choose. We conclude that in addition to viewing angles and instrumental resolutions, NEI further complicates the detection of FSMs in spectroscopic measurements of coronal loops in the EUV passband.

We have studied numerically the shadows of a non-Kerr rotating compact object with quadrupole mass moment, which belongs to Manko-Novikov family. The non-integrable photon motion caused by quadrupole mass moment affects sharply the shadow of the compact object. As the deviation parameter related to quadrupole mass moment is negative, the shadow of compact object is prolate and there are two disconnected main shadows with eyebrows located symmetrically on both sides of the equatorial plane. As the deviation parameter is positive, the shadow becomes oblate and the main shadow is joined together in the equatorial plane. Moreover, in this positive cases, there is a disorder region in the left of shadow which increases with the quadrupole-deviation parameter. Interestingly, we also find that Einstein ring is broken as the deviation from Kerr metric is larger than a certain critical value. This critical value decreases with the rotation parameter of black hole. Especially, the observer on the direction of rotation axis will find some concentric bright rings in the black disc. Finally, supposing that the gravitational field of the supermassive central object of the galaxy described by this metric, we estimated the numerical values of the observables for the black hole shadow.

We have studied 4265 giant pulses (GPs) from the millisecond pulsar B1937+21; the largest-ever sample gathered for this pulsar, in observations made with the Large European Array for Pulsars. The pulse energy distribution of GPs associated with the interpulse are well-described by a power law, with index $\alpha = -3.99 \pm 0.04$, while those associated with the main pulse are best-described by a broken power law, with the break occurring at $\sim7$ Jy $\mu$s, with power law indices $\alpha_{\text{low}} = -3.48 \pm 0.04$ and $\alpha_{\text{high}} = -2.10 \pm 0.09$. The modulation indices of the GP emission are measured, which are found to vary by $\sim0.5$ at pulse phases close to the centre of the GP phase distributions. We find the frequency-resolved structure of GPs to vary significantly, and in a manner that cannot be attributed to the interstellar medium influence on the observed pulses. We examine the distribution of polarisation fractions of the GPs and find no correlation between GP emission phase and fractional polarisation. We use the GPs to time PSR B1937+21 and although the achievable time of arrival precision of the GPs is approximately a factor of two greater than that of the average pulse profile, there is a negligible difference in the precision of the overall timing solution when using the GPs.

We present a detailed overview of the cosmological surveys that will be carried out with Phase 1 of the Square Kilometre Array (SKA1), and the science that they will enable. We highlight three main surveys: a medium-deep continuum weak lensing and low-redshift spectroscopic HI galaxy survey over 5,000 sqdeg; a wide and deep continuum galaxy and HI intensity mapping survey over 20,000 sqdeg from z = 0.35 - 3; and a deep, high-redshift HI intensity mapping survey over 100 sqdeg from z = 3 - 6. Taken together, these surveys will achieve an array of important scientific goals: measuring the equation of state of dark energy out to z ~ 3 with percent-level precision measurements of the cosmic expansion rate; constraining possible deviations from General Relativity on cosmological scales by measuring the growth rate of structure through multiple independent methods; mapping the structure of the Universe on the largest accessible scales, thus constraining fundamental properties such as isotropy, homogeneity, and non-Gaussianity; and measuring the HI density and bias out to z = 6. These surveys will also provide highly complementary clustering and weak lensing measurements that have independent systematic uncertainties to those of optical surveys like LSST and Euclid, leading to a multitude of synergies that can improve constraints significantly beyond what optical or radio surveys can achieve on their own. This document, the 2018 Red Book, provides reference technical specifications, cosmological parameter forecasts, and an overview of relevant systematic effects for the three key surveys, and will be regularly updated by the Cosmology Science Working Group in the run up to start of operations and the Key Science Programme of SKA1.

There has been considerable interest in sausage modes in photospheric waveguides like pores and sunspots, and slow surface sausage modes (SSSMs) have been suggested to damp ufficiently rapidly to account for chromospheric heating. Working in the framework of linear resistive magnetohydrodynamics, we examine how efficient electric resistivity and resonant absorption in the cusp continuum can be for damping SSSMs in a photospheric waveguide with equilibrium parameters compatible with recent measurements of a photospheric pore. For SSSMs with the measured wavelength, we find that the damping rate due to the cusp resonance is substantially less strong than theoretically expected with the thin-boundary approximation. The damping-time-to-period ratio ($\tau/P$) we derive for standing modes, equivalent to the damping-length-to-wavelength ratio for propagating modes given the extremely weak dispersion, can reach only $\sim 180$. However, the accepted values for electric resistivity ($\eta$) correspond to a regime where both the cusp resonance and resistivity play a role. The values for $\tau/P$ attained at the largest allowed $\eta$ may reach $\sim 30$. We conclude that electric resistivity can be considerably more efficient than the cusp resonance for damping SSSMs in the pore in question, and it needs to be incorporated into future studies on the damping of SSSMs in photospheric waveguides in general.

Narrow-line Seyfert 1 galaxies (NLS1s) are a subclass of active galactic nuclei (AGN). It is often believed that these AGN have small black hole mass, which is responsible for the narrowness of the permitted lines. They are also characterised by a high accretion rate, typically closer to the Eddington limit. Nevertheless, narrow permitted lines might also be caused by a disk-like broad-line region (BLR) viewed pole-on. This class of objects presents strong X-ray emission, which is characterised by a very steep spectral index described by a single power law. In particular, some of them exhibit particular features around the iron K-shell energy at 6-8 keV. Recently, this different spectral behaviour was attributed to inclination. In this work we are going to analyse optical spectra to measure in different ways the width of H$\beta$, which is another potential inclination indicator. Our aim is to search for a correlation between the high-energy spectral complexity and FWHM of H$\beta$, in order to verify whether or not the broad-line region could be flattened.

We present a new accurate sample of narrow-line Seyfert 1 galaxies (NLS1s) in the southern hemisphere from the Six-degree Field Galaxy Survey (6dFGS). Based on the optical spectral features, 167 sources were classified as NLS1s. We derived flux-calibrated spectra in the sample for the first time. Strong luminosity correlations between the continuum and the emission lines were found. We estimated their central black hole masses and Eddington ratios, which are lying in a typical range of NLS1s. In the sample, 23 NLS1s were detected at radio frequencies and 12 of them are radio-loud. We publish the X-ray data analysis of a campaign of observations carried out by the Swift X-ray Telescope.

In the context of space and astrophysical plasma turbulence and particle heating, several vocabularies emerge for estimating turbulent energy dissipation rate, including Kolmogorov-Yaglom third-order law and, in its various forms, $\boldsymbol{j}\cdot\boldsymbol{E}$ (work done by the electromagnetic field on particles), and $-\left( \boldsymbol{P} \cdot \nabla \right) \cdot \boldsymbol{u}$ (pressure-strain interaction), to name a couple. It is now understood that these energy transfer channels, to some extent, are correlated with coherent structures. In particular, we find that different energy dissipation proxies, although not point-wise correlated, are concentrated in proximity to each other, for which they decorrelate in a few $d_i$(s). However, the energy dissipation proxies dominate at different scales. For example, there is an inertial range over which the third-order law is meaningful. Contributions from scale bands stemming from scale-dependent spatial filtering show that, the energy exchange through $\boldsymbol{j}\cdot\boldsymbol{E}$ mainly results from large scales, while the energy conversion from fluid flow to internal through $-\left( \boldsymbol{P} \cdot \nabla \right) \cdot \boldsymbol{u}$ dominates at small scales.

We present a new accurate catalog of narrow-line Seyfert 1 galaxies (NLS1s) in the southern hemisphere from the Six-degree Field Galaxy Survey (6dFGS) final data release, which is currently the most extensive spectroscopic survey available in the southern sky whose database has not yet been systematically explored. We classified 167 sources as NLS1s based on their optical spectral properties. We derived flux-calibrated spectra for the first time that the 6dFGS indeed does not provide. By analyzing these spectra, we obtained strong correlations between the monochromatic luminosity at 5100 A and the luminosity of Hbeta and [O III]5007 lines. The estimated central black hole mass and Eddington ratio have an average value of 8.6 x 10^6 M_Sun and 0.96 L_Edd respectively, which is a typical value for NLS1s. In the sample, 23 (13.8%) NLS1s were detected at radio frequencies, and 12 (7.0%) of them are radio-loud. Our results confirmed that radio-loud sources tend to have a higher redshift, more massive black hole, and higher radio and optical luminosity than radio-quiet sources.

Using parallaxes from Gaia DR2, we estimate the distance to the globular clusters 47 Tuc and NGC 362, taking advantage of the background stars in the Small Magellanic Cloud and quasars to account for various parallax systematics. We found the parallax to be dependent on the Gaia DR2 G apparent magnitude for stars with $13<\textrm{G}<18$ , where brighter stars have a lower parallax zero point than fainter stars. The distance to 47 Tuc was found to be $4.45\pm0.01\pm0.12$ kpc, and for NGC 362 $8.54\pm0.20\pm0.44$ kpc with random and systematic errors listed respectively. This is the first time a precise distance measurement directly using parallaxes has been done for either of these two globular clusters.

Thanks to the execution of extensive spectroscopic surveys that have covered large fractions of the sky down to magnitudes as faint as $i \approx 19$, it has been possible to identify several narrow-line Seyfert 1 galaxies (NLS1s) and to investigate their properties over a large range of the electro-magnetic spectrum. The interpretation of their nature, however, is still hampered by the statistical uncertainties related to the way in which NLS1 candidates are selected. In this contribution, we present a study to detect and to model emission lines in optical spectra extracted from the Sloan Digital Sky Survey (SDSS), adopting the most proper strategy to identify the source of line excitation and to produce a detailed model with measurements of several emission-line parameters. We demonstrate the application of this technique to explore fundamental questions, such as the presence of gas and dust around the core of AGNs and the spectral energy distribution of their ionizing radiation. We compare the spectral properties of NLS1s with those of other type 1 active galaxies and we summarize the potential of this approach to identify NLS1s in present day and future spectroscopic surveys. We finally consider the implications of multi-frequency data analysis in the debate concerning the intrinsic nature of NLS1s.

The cosmic-ray (CR) knee and the compositions contain abundant information for probing the CR's origin, acceleration and propagation mechanisms, as well as the frontier of the fundamental physics. Realizing that major proposals toward the knee's shape can be divided into two categories: the rigidity-dependent (also Z-dependent) knee and the mass-dependent (also A-dependent) knee, where the former one relates to the acceleration or the propagation mechanisms, and the other one is often associated with the new physics, it is essential to precisely measure the individual compositions. Benefit from the high altitude and hybrid detection methods, the LHAASO experiment has the ability in determining the individual component and brings us an opportunity in discriminating these two models. We test this expected ability of LHAASO from 100 TeV to 10 PeV with 3-year observation. And find the dominant component at the knee is essential to this issue, while much heavier nuclei occupying the knee leads to higher significance. In the analysis, the He-dominant knee under the A-dependent case can be recognized at the significance about 6.6 $\sigma$, while the P-dominant knee under the Z-dependent case will be classified with 2 $\sigma$ significance.

The Gemini Infrared Multi-Object Spectrograph (GIRMOS) is a powerful new instrument being built to facility-class standards for the Gemini telescope. It takes advantage of the latest developments in adaptive optics and integral field spectrographs. GIRMOS will carry out simultaneous high-angular-resolution, spatially-resolved infrared ($1-2.4$ $\mu$m) spectroscopy of four objects within a two-arcminute field-of-regard by taking advantage of multi-object adaptive optics. This capability does not currently exist anywhere in the world and therefore offers significant scientific gains over a very broad range of topics in astronomical research. For example, current programs for high redshift galaxies are pushing the limits of what is possible with infrared spectroscopy at $8-10$-meter class facilities by requiring up to several nights of observing time per target. Therefore, the observation of multiple objects simultaneously with adaptive optics is absolutely necessary to make effective use of telescope time and obtain statistically significant samples for high redshift science. With an expected commissioning date of 2023, GIRMOS's capabilities will also make it a key followup instrument for the James Webb Space Telescope when it is launched in 2021, as well as a true scientific and technical pathfinder for future Thirty Meter Telescope (TMT) multi-object spectroscopic instrumentation. In this paper, we will present an overview of this instrument's capabilities and overall architecture. We also highlight how this instrument lays the ground work for a future TMT early-light instrument.

As we demonstrated in Paper I, the quenched fractions of central and satellite galaxies as function of halo mass are extremely similar, as long as one controls for stellar mass. The same holds for the quenched fractions as a function of central velocity dispersion, which is tightly correlated with black hole mass, as long as one controls for both stellar and halo mass. Here we use mock galaxy catalogs constructed from the latest semi-analytic model, L-GALAXIES, and the state-of-the-art hydrodynamical simulation, EAGLE, to investigate whether these models can reproduce the trends seen in the data. We also check how the group finder used to identify centrals and satellites impacts our results. We find that L-GALAXIES fails to reproduce the trends. The predicted quenched fraction of central galaxies increases sharply with halo mass around $10^{12.5}h^{-1}M_{\odot}$ and with black hole mass around $\sim10^{6.5}M_{\odot}$, while the predicted quenched fraction of satellites increases with both halo and black hole masses gradually. In contrast, centrals and satellites in EAGLE follow almost the same trend as seen in the data. We discuss the implications of our results for how feedback processes regulate galaxy quenching.

This paper reports on the observation of the sidereal large-scale anisotropy of cosmic rays using data collected by the ARGO-YBJ experiment over 5 years (2008$-$2012). This analysis extends previous work limited to the period from 2008 January to 2009 December,near the minimum of solar activity between cycles 23 and 24.With the new data sample the period of solar cycle 24 from near minimum to maximum is investigated. A new method is used to improve the energy reconstruction, allowing us to cover a much wider energy range, from 4 to 520 TeV. Below 100 TeV, the anisotropy is dominated by two wide regions, the so-called "tail-in" and "loss-cone" features. At higher energies, a dramatic change of the morphology is confirmed. The yearly time dependence of the anisotropy is investigated. Finally, no noticeable variation of cosmic-ray anisotropy with solar activity is observed for a median energy of 7 TeV.

We investigate the quenching properties of central and satellite galaxies, utilizing the halo masses and central-satellite identifications from the SDSS galaxy group catalog of Yang et al. We find that the quenched fractions of centrals and satellites of similar stellar masses have similar dependence on host halo mass. The similarity of the two populations is also found in terms of specific star formation rate and 4000 \AA\ break. The quenched fractions of centrals and satellites of similar masses show similar dependencies on bulge-to-total light ratio, central velocity dispersion and halo-centric distance in halos of given halo masses. The prevalence of optical/radio-loud AGNs is found to be similar for centrals and satellites at given stellar masses. All these findings strongly suggest that centrals and satellites of similar masses experience similar quenching processes in their host halos. We discuss implications of our results for the understanding of galaxy quenching.

We report the results of a neutrino search in Super-Kamiokande for coincident signals with the first detected gravitational wave produced by a binary neutron star merger, GW170817, which was followed by a short gamma-ray burst, GRB170817A, and a kilonova/macronova. We searched for coincident neutrino events in the range from 3.5 MeV to $\sim$100 PeV, in a time window $\pm$500 seconds around the gravitational wave detection time, as well as during a 14-day period after the detection. No significant neutrino signal was observed for either time window. We calculated 90% confidence level upper limits on the neutrino fluence for GW170817. From the upward-going-muon events in the energy region above 1.6 GeV, the neutrino fluence limit is $16.0^{+0.7}_{-0.6}$ ($21.3^{+1.1}_{-0.8}$) cm$^{-2}$ for muon neutrinos (muon antineutrinos), with an error range of $\pm5^{\circ}$ around the zenith angle of NGC4993, and the energy spectrum is under the assumption of an index of $-2$. The fluence limit for neutrino energies less than 100 MeV, for which the emission mechanism would be different than for higher-energy neutrinos, is also calculated. It is $6.6 \times 10^7$ cm$^{-2}$ for anti-electron neutrinos under the assumption of a Fermi-Dirac spectrum with average energy of 20 MeV.

We report the first results of a survey on 74 narrow-line Seyfert 1 galaxies (NLS1s) carried out in 2015 with the Karl G. Jansky Very Large Array (JVLA) at 5 GHz in A-configuration. So far, this is the largest survey aimed to image the radio continuum of NLS1s. We produced radio maps in order to compare the general properties of three different samples of objects: radio-quiet NLS1s (RQNLS1s), steep-spectrum radio-loud NLS1s (S-NLS1s), and flat-spectrum radio-loud NLS1s (F-NLS1s). We find that the three classes correspond to different radio morphologies, with F-NLS1s being more compact, and RQNLS1s often showing diffuse emission on kpc scales. We also find that F-NLS1s might be low-luminosity and possibly young blazars, and that S-NLS1s are part of the parent population of F-NLS1s. Dedicated studies to RQNLS1s are needed in order to fully understand their role in the unification pictures.

A correlation between the secondary cosmic ray flux and the near-earth electric field intensity, measured during thunderstorms, has been found by analyzing the data of the ARGO-YBJ experiment, a full coverage air shower array located at the Yangbajing Cosmic Ray Laboratory (4300 m a. s. l., Tibet, China). The counting rates of showers with different particle multiplicities, have been found to be strongly dependent upon the intensity and polarity of the electric field measured during the course of 15 thunderstorms. In negative electric fields (i.e. accelerating negative charges downwards), the counting rates increase with increasing electric field strength. In positive fields, the rates decrease with field intensity until a certain value of the field EFmin (whose value depends on the event multiplicity), above which the rates begin increasing. By using Monte Carlo simulations, we found that this peculiar behavior can be well described by the presence of an electric field in a layer of thickness of a few hundred meters in the atmosphere above the detector, which accelerates/decelerates the secondary shower particles of opposite charge, modifying the number of particles with energy exceeding the detector threshold. These results, for the first time, give a consistent explanation for the origin of the variation of the electron/positron flux observed for decades by high altitude cosmic ray detectors during thunderstorms.

Impulsively generated sausage wave trains in coronal structures are important for interpreting a substantial number of observations of quasi-periodic signals with quasi-periods of order seconds. We have previously shown that the Morlet spectra of these wave trains in coronal tubes depend crucially on the dispersive properties of trapped sausage waves, the existence of cutoff axial wavenumbers and the monotonicity of the dependence of the axial group speed on the axial wavenumber in particular. This study examines the difference a slab geometry may introduce, for which purpose we conduct a comprehensive eigenmode analysis, both analytically and numerically, on trapped sausage modes in coronal slabs with a considerable number of density profiles. For the profile descriptions examined, coronal slabs can trap sausage waves with longer axial wavelengths, and the group speed approaches the internal Alfv\'en speed more rapidly at large wavenumbers in the cylindrical case. However, common to both geometries, cutoff wavenumbers exist only when the density profile falls sufficiently rapidly at distances far from coronal structures. Likewise, the monotonicity of the group speed curves depends critically on the profile steepness right at the structure axis. Furthermore, the Morlet spectra of the wave trains are shaped by the group speed curves for coronal slabs and tubes alike. Consequently, we conclude that these spectra have the potential for telling the sub-resolution density structuring inside coronal structures, although their detection requires an instrumental cadence of better than $\sim 1$ second.

We examine the dispersive properties of linear fast standing modes in transversely nonuniform solar coronal slabs with finite gas pressure, or, equivalently, finite plasma beta. We derive a generic dispersion relation governing fast waves in coronal slabs for which the continuous transverse distributions of the physical parameters comprise a uniform core, a uniform external medium, and a transition layer (TL) in between. The profiles in the TL are allowed to be essentially arbitrary. Restricting ourselves to the first several branches of fast modes, which are of most interest from the observational standpoint, we find that a finite plasma beta plays an at most marginal role in influencing the periods ($P$), damping times ($\tau$), and critical longitudinal wavenumbers ($k_{\rm c}$), when both $P$ and $\tau$ are measured in units of the transverse fast time. However, these parameters are in general significantly affected by how the TL profiles are described. We conclude that, for typical coronal structures, the dispersive properties of the first several branches of fast standing modes can be evaluated with the much simpler theory for cold slabs provided that the transverse profiles are properly addressed and the transverse Alfv\'en time in cold MHD is replaced with the transverse fast time.

In this paper, I investigated the Doppler magnification effect in the flux-limited galaxy number counts with finite redshift bin width. In contradiction with our intuition, the number counts correction formula for different redshift bin are different. For the redshift window as delta function and constant function, the correction formulas have been derived analytically in this paper. These two windows correspond to two extreme cases(i.e. extremely narrow bin and extremely broad bin). An simulation have been implemented to test these formulas. The simulation results indicate the perturbation changes smoothly from one extreme case to another while extending the redshift bin width from half bin width $0.05$ to $0.4$. As a result, Doppler magnification caused number counts perturbation of finite redshift bin width can NOT be compute via redshift integration of the perturbation formula derived with delta function redshift window. It is only a good approximation when the bin width is small. These results are important for galaxy redshift survey kinematic dipole estimation, and meaningful for the relativistic galaxy number counts angular power spectrum estimation, such as CLASSgal and CAMB source.

We present the results of a search in the Super-Kamiokande (SK) detector for excesses of neutrinos with energies above a few GeV that are in the direction of the track events reported in IceCube. Data from all SK phases (SK-I through SK-IV) were used, spanning a period from April 1996 to April 2016 and corresponding to an exposure of 225 kilotonne-years . We considered the 14 IceCube track events from a data set with 1347 livetime days taken from 2010 to 2014. We use Poisson counting to determine if there is an excess of neutrinos detected in SK in a 10 degree search cone (5 degrees for the highest energy data set) around the reconstructed direction of the IceCube event. No significant excess was found in any of the search directions we examined. We also looked for coincidences with a recently reported IceCube multiplet event. No events were detected within a $\pm$ 500 s time window around the first detected event, and no significant excess was seen from that direction over the lifetime of SK.

In the concordance model of the Universe, the matter distribution - as observed in galaxy number counts or the intensity of line emission (such as the 21cm line of neutral hydrogen) - should have a kinematic dipole due to the Sun's motion relative to the CMB rest-frame. This dipole should be aligned with the kinematic dipole in the CMB temperature. Accurate measurement of the direction of the matter dipole will become possible with future galaxy surveys, and this will be a critical test of the foundations of the concordance model. The amplitude of the matter dipole is also a potential cosmological probe. We derive formulas for the amplitude of the kinematic dipole in galaxy redshift and intensity mapping surveys, taking into account the Doppler, aberration and other relativistic effects. The amplitude of the matter dipole can be significantly larger than that of the CMB dipole. Its redshift dependence encodes information on the evolution of the Universe and on the tracers, and we discuss possible ways to determine the amplitude.

The Large European Array for Pulsars (LEAP) is a European Pulsar Timing Array project that combines the Lovell, Effelsberg, Nan\c{c}ay, Sardinia, and Westerbork radio telescopes into a single tied-array, and makes monthly observations of a set of millisecond pulsars (MSPs). The overview of our experiment is presented in Bassa et al. (2016). Baseband data are recorded at a central frequency of 1396 MHz and a bandwidth of 128 MHz at each telescope, and are correlated offline on a cluster at Jodrell Bank Observatory using a purpose-built correlator, detailed in Smits et al. (2017). LEAP offers a substantial increase in sensitivity over that of the individual telescopes, and can operate in timing and imaging modes (notably in observations of the galactic centre radio magnetar; Wucknitz 2015). To date, 4 years of observations have been reduced. Here, we report on the scientific projects which have made use of LEAP data.

We present final statistics from a survey for intervening MgII absorption towards 100 quasars with emission redshifts between $z=3.55$ and $z=7.08$. Using infrared spectra from Magellan/FIRE, we detect 279 cosmological MgII absorbers, and confirm that the incidence rate of $W_r>0.3 \AA$ MgII absorption per comoving path length does not evolve measurably between $z=0.25$ and $z=7$. This is consistent with our detection of seven new MgII systems at $z>6$, a redshift range that was not covered in prior searches. Restricting to relatively strong MgII systems ($W_r>1$\AA), there is significant evidence for redshift evolution. These systems roughly double in number density between $z=0$ and $z=2$-$3$, but decline by an order of magnitude from this peak by $z\sim 6$. This evolution mirrors that of the global star formation rate density, which could reflect a connection between star formation feedback and strong MgII absorbers. We compared our results to the Illustris cosmological simulation at $z=2$-$4$ by assigning absorption to catalogued dark-matter halos and by direct extraction of spectra from the simulation volume. To reproduce our results using the halo catalogs, we require circumgalactic (CGM) MgII envelopes within halos of progressively smaller mass at earlier times. This occurs naturally if we define the lower integration cutoff using SFR rather than mass. MgII profiles calculated directly from the Illustris volume yield far too few strong absorbers. We argue that this arises from unresolved phase space structure of CGM gas, particularly from turbulent velocities on sub-mesh scales. The presence of CGM MgII at $z>6$-- just $\sim 250$ Myr after the reionization redshift implied by Planck--suggests that enrichment of intra-halo gas may have begun before the presumed host galaxies' stellar populations were mature and dynamically relaxed. [abridged]

The absence of any confirmative signals from extensive DM searching motivates us to go beyond the conventional WIMPs scenario. The feebly interacting massive particles (FIMPs) paradigm provides a good alternative which, despite of its feebly interaction with the thermal particles, still could correctly produce relic abundance without conventional DM signals. The Infrared-FIMP based on the renormalizable operators is usually suffering the very tiny coupling drawback, which can be overcome in the UltraViolet-FIMP scenario based on high dimensional effective operators. However, it is sensitive to the history of the very early Universe. The previous works terminates this sensitivity at the reheating temperature $T_{RH}$. We, motivated by its UV-sensitivity, investigate the effects from the even earlier Universe, reheating era. We find that in the usual case with $T_{RH}\gg m_{\rm DM}$, the production rate during reheating is very small as long as the effective operators dimension $d \leq 8$. Besides, we consider the contribution from the mediator, which may be produced during reheating. Moreover, we study the situation when $T_{RH}$ is even lower than $m_{\rm DM}$ and DM can be directly produced during reheating if its mass does not exceed $T_{MAX}$.

We examine the quenched fraction of central and satellite galaxies as a function of galaxy stellar mass, halo mass, and the matter density of their large scale environment. Matter densities are inferred from our ELUCID simulation, a constrained simulation of local Universe sampled by SDSS, while halo masses and central/satellite classification are taken from the galaxy group catalog of Yang et al. The quenched fraction for the total population increases systematically with the three quantities. We find that the `environmental quenching efficiency', which quantifies the quenched fraction as function of halo mass, is independent of stellar mass. And this independence is the origin of the stellar mass-independence of density-based quenching efficiency, found in previous studies. Considering centrals and satellites separately, we find that the two populations follow similar correlations of quenching efficiency with halo mass and stellar mass, suggesting that they have experienced similar quenching processes in their host halo. We demonstrate that satellite quenching alone cannot account for the environmental quenching efficiency of the total galaxy population and the difference between the two populations found previously mainly arises from the fact that centrals and satellites of the same stellar mass reside, on average, in halos of different mass. After removing these halo-mass and stellar-mass effects, there remains a weak, but significant, residual dependence on environmental density, which is eliminated when halo assembly bias is taken into account. Our results therefore indicate that halo mass is the prime environmental parameter that regulates the quenching of both centrals and satellites.

The spectra of Active Galactic Nuclei (AGNs) are often characterized by a wealth of emission lines with different profiles and intensity ratios that led to a complicated classification. Their electro-magnetic radiation spans more than 10 orders of magnitude in frequency. In spite of the differences between various classes, the origin of their activity is attributed to a combination of emitting components, surrounding an accreting Super Massive Black Hole, in the so called Unified Model. Currently, the execution of sky surveys, with instruments operating at various frequencies, provides the possibility to detect and to investigate the properties of AGNs on very large statistical samples. Thanks to the spectroscopic surveys that allow investigation of many objects, we have the opportunity to place new constraints on the nature and evolution of AGNs. In this contribution we present the results obtained by working on multi-frequency data and we discuss their relations with the available optical spectra. We compare our findings with the AGN Unified Model predictions, and we present a revised technique to select AGNs of different types from other line emitting objects. We discuss the multi-frequency properties in terms of the innermost structures of the sources.

Active Galactic Nuclei (AGNs) are characterized by emission of radiation over more than 10 orders of magnitude in frequency. Therefore, the execution of extensive surveys of the sky, with different types of detectors, is providing the attractive possibility to identify and to investigate the properties of AGNs on very large statistical samples. Thanks to the large spectroscopic surveys that allow detailed investigation of many of these sources, we have the opportunity to place new constraints on the nature and evolution of AGNs and to investigate their relations with the host systems. In this contribution we present the results that can be obtained by using a new interactive catalogue that we developed to investigate the range of AGN spectral energy distributions (SEDs). We present simple SED models based on data collected in the catalogue and discuss their relations with optical spectra obtained by follow up observations. We compare our findings with the expectations based on the AGN Unification Model, and we discuss the perspectives of multi-wavelength approaches to address AGN related processes such as black hole accretion and acceleration of relativistic jets.

Multi-wavelength astronomical studies brings a wealth of science within reach. One way to achieve a cross-wavelength analysis is via `stacking', i.e. combining precise positional information from an image at one wavelength with data from one at another wavelength in order to extract source-flux distributions and other derived quantities. For the first time we extend stacking to include the effects of confusion. We develop our algorithm in a fully Bayesian framework and apply it to the Square Kilometre Array Design Study (SKADS) simulation in order to extract galaxy number counts. Previous studies have shown that recovered source counts are highly biased high when confusion is non-negligible. With this new method, source counts are returned correctly. We also describe a novel estimator for quantifying the impact of confusion on stacking analyses. This method is an essential step in exploiting scientific return for upcoming deep radio surveys, e.g. MIGHTEE on MeerKAT.

Accretion dynamics in the forming of young stars is still object of debate because of limitations in observations and modelling. Through scaled laboratory experiments of collimated plasma accretion onto a solid in the presence of a magnetic field, we open first window on this phenomenon by tracking, with spatial and temporal resolution, the dynamics of the system and simultaneously measuring multiband emissions. We observe in these experiments that matter, upon impact, is laterally ejected from the solid surface, then refocused by the magnetic field toward the incoming stream. Such ejected matter forms a plasma shell that envelops the shocked core, reducing escaped X-ray emission. This demonstrates one possible structure reconciling current discrepancies between mass accretion rates derived from X-ray and optical observations.

Pulsar timing arrays (PTAs) are presently the only means to search for the gravitational wave stochastic background from supermassive black hole binary populations, considered to be within the grasp of current or near future observations. However, the stringent upperlimit set by the Parkes PTA (Shannon et al. 2013, 2015) has been interpreted as excluding at $> 90\%$ confidence the current paradigm of binary assembly through galaxy mergers and hardening via stellar interactions, suggesting evolution is accelerated (by stars and/or gas) or stalled. Using Bayesian hierarchical modelling, we consider implications of this upperlimit for a comprehensive range of astrophysical scenarios, without invoking stalling nor more exotic physical processes. We find they are fully consistent with the upperlimit, but (weak) bounds on population parameters can be inferred. Bayes factors between models vary between $\approx 1.03$ -- $5.81$ and Kullback-Leibler divergences between characteristic amplitude prior and posterior lie between $0.37$ -- $0.85$. Considering prior astrophysical information on galaxy merger rates, recent upwards revisions of the black hole-galaxy bulge mass relation (Kormendy & Ho 2013) are disfavoured at $1.6\sigma$ against lighter models (e.g. Shankar et al. 2016). We also show, if no detection is achieved once sensitivity improves by an order of magnitude, the most optimistic scenario is disfavoured at $3.9\sigma$.

We present a fast and versatile method to calculate the characteristic spectrum $h_c$ of the gravitational wave background (GWB) emitted by a population of eccentric massive black hole binaries (MBHBs). We fit the spectrum of a reference MBHB with a simple analytic function and show that the spectrum of any other MBHB can be derived from this reference spectrum via simple scalings of mass, redshift and frequency. We then apply our calculation to a realistic population of MBHBs evolving via 3-body scattering of stars in galactic nuclei. We demonstrate that our analytic prescription satisfactorily describes the signal in the frequency band relevant to pulsar timing array (PTA) observations. Finally we model the high frequency steepening of the GWB to provide a complete description of the features characterizing the spectrum. For typical stellar distributions observed in massive galaxies, our calculation shows that 3-body scattering alone is unlikely to affect the GWB in the PTA band and a low frequency turnover in the spectrum is caused primarily by high eccentricities.

We consider the inverse problem in pulsar timing array (PTA) analysis, investigating what astrophysical information about the underlying massive black hole binary (MBHB) population can be recovered from the detection of a stochastic gravitational wave background (GWB). We employ a physically motivated model that connects the GWB spectrum to a series of parameters describing the underlying redshift evolution of the MBHB mass function and to the typical eccentricity they acquire while interacting with the dense environment of post merger galactic nuclei. This allows the folding in of information about the spectral shape of the GWB into the analysis. The priors on the model parameters are assumed to be uninformative and consistent with the current lack of secure observations of sub-parsec MBHBs. We explore the implications of current upper limits as well as of future detections with a variety of PTA configurations. We confirm our previous finding that current upper limits can only place an upper bound on the overall MBHB merger rate. Depending on the properties of the array, future detections can also constrain several MBHB population models at different degrees of fidelity. In particular, a simultaneous detection of a steepening of the spectrum at high frequency and a bending at low frequency will place strong constraints on both the MBHB mass function and on the typical eccentricity of inspiralling MBHBs, providing insights on MBHB astrophysics unlikely to be achievable by any other means.

The ARGO-YBJ experiment, a full coverage extensive air shower (EAS) detector located at high altitude (4300 m a.s.l.) in Tibet, China, has smoothly taken data, with very high stability, since November 2007 to the beginning of 2013. The array consisted of a carpet of about 7000 m$^2$ Resistive Plate Chambers (RPCs) operated in streamer mode and equipped with both digital and analog readout, providing the measurement of particle densities up to few particles per cm$^2$. The unique detector features (full coverage, readout granularity, wide dynamic range, etc) and location (very high altitude) allowed a detailed study of the lateral density profile of charged particles at ground very close to the shower axis and its description by a proper lateral distribution function (LDF). In particular, the information collected in the first 10 m from the shower axis have been shown to provide a very effective tool for the determination of the shower development stage ("age") in the energy range 50 TeV - 10 PeV. The sensitivity of the age parameter to the mass composition of primary Cosmic Rays is also discussed.

Kinetic plasma turbulence cascade spans multiple scales ranging from macroscopic fluid flow to sub-electron scales. Mechanisms that dissipate large scale energy, terminate the inertial range cascade and convert kinetic energy into heat are hotly debated. Here we revisit these puzzles using fully kinetic simulation. By performing scale-dependent spatial filtering on the Vlasov equation, we extract information at prescribed scales and introduce several energy transfer functions. This approach allows highly inhomogeneous energy cascade to be quantified as it proceeds down to kinetic scales. The pressure work, $-\left( \boldsymbol{P} \cdot \nabla \right) \cdot \boldsymbol{u}$, can trigger a channel of the energy conversion between fluid flow and random motions, which is a collision-free generalization of the viscous dissipation in collisional fluid. Both the energy transfer and the pressure work are strongly correlated with velocity gradients.

The ARGO-YBJ detector, located at the Yangbajing Cosmic Ray Laboratory (4300 m a. s. l., Tibet, China), was a full coverage air shower array dedicated to gamma ray astronomy and cosmic ray studies. The wide field of view (~ 2 sr) and high duty cycle (> 86%), made ARGO-YBJ suitable to search for short and unexpected gamma ray emissions like gamma ray bursts (GRBs). Between 2007 November 6 and 2013 February 7, 156 satellite-triggered GRBs (24 of them with known redshift) occurred within the ARGO-YBJ field of view. A search for possible emission associated to these GRBs has been made in the two energy ranges 10-100 GeV and 10-1000 GeV. No significant excess has been found in time coincidence with the satellite detections nor in a time window of one hour after the bursts. Taking into account the EBL absorption, upper limits to the energy fluence at 99% of confidence level have been evaluated,with values ranging from ~ 10-5 erg cm-2 to ~10-1 erg cm-2.

The behavior of the axial group speeds of trapped sausage modes plays an important role in determining impulsively generated wave trains, which have often been invoked to account for quasi-periodic signals with quasi-periods of order seconds in a considerable number of coronal structures. We conduct a comprehensive eigenmode analysis, both analytically and numerically, on the dispersive properties of sausage modes in pressureless tubes with three families of continuous radial density profiles. We find a rich variety of the dependence on the axial wavenumber $k$ of the axial group speed $v_{\rm gr}$. Depending on the density contrast and profile steepness as well as on the detailed profile description, the $v_{\rm gr}-k$ curves either possess or do not possess cutoff wavenumbers, and they can behave in either a monotonical or non-monotonical manner. With time-dependent simulations, we further show that this rich variety of the group speed characteristics heavily influences the temporal evolution and Morlet spectra of impulsively generated wave trains. In particular, the Morlet spectra can look substantially different from "crazy tadpoles" found for the much-studied discontinuous density profiles. We conclude that it is necessary to re-examine available high-cadence data to look for the rich set of temporal and spectral features, which can be employed to discriminate between the unknown forms of the density distributions transverse to coronal structures.

Although still under debate, observations generally suggest that dwarf spheroidal (dSph) galaxies exhibit large constant-density cores in the centers, which can hardly be explained by dissipationless cold dark matter simulations without baryonic feedback. Wave dark matter (${\psi {\rm DM}}$), characterized by a single parameter, the dark matter particle mass $m_{\psi}$, predicts a central soliton core in every galaxy arising from quantum pressure against gravity. Here we apply Jeans analysis assuming a soliton core profile to the kinematic data of eight classical dSphs so as to constrain $m_{\psi}$, and obtain $m_{\psi}=1.18_{-0.24}^{+0.28}\times10^{-22}{\,\rm eV}$ and $m_{\psi}=1.79_{-0.33}^{+0.35}\times10^{-22}{\,\rm eV}~(2\sigma)$ using the observational data sets of Walker et al. (2007) and Walker et al. (2009b), respectively. We show that the estimate of $m_{\psi}$ is sensitive to the dSphs kinematic data sets and is robust to various models of stellar density profile. We also consider multiple stellar subpopulations in dSphs and find consistent results. This mass range of $m_{\psi}$ is in good agreement with other independent estimates, such as the high-redshift luminosity functions, the reionization history, and the Thomson optical depth to the cosmic microwave background.

The frequency dependence of the longitudinal group speeds of trapped sausage waves plays an important role in determining impulsively generated wave trains, which have often been invoked to account for quasi-periodic signals in coronal loops. We examine how the group speeds ($v_{\rm gr}$) depend on angular frequency ($\omega$) for sausage modes in pressureless coronal tubes with continuous transverse density distributions by solving the dispersion relation pertinent to the case where the density inhomogeneity of arbitrary form takes place in a transition layer of arbitrary thickness. We find that in addition to the transverse lengthscale $l$ and density contrast $\rho_{\rm i}/\rho_{\rm e}$, the group speed behavior depends also on the detailed form of the density inhomogeneity. For parabolic profiles, $v_{\rm gr}$ always decreases with $\omega$ first before increasing again, as happens for the much studied top-hat profiles. For linear profiles, however, the behavior of the $\omega-v_{\rm gr}$ curves is more complex. When $\rho_{\rm i}/\rho_{\rm e} \lesssim 6$, the curves become monotonical for large values of $l$. On the other hand, for higher density contrasts, a local maximum $v_{\rm gr}^{\rm max}$ exists in addition to a local minimum $v_{\rm gr}^{\rm min}$ when coronal tubes are diffuse. With time-dependent computations, we show that the different behavior of group speed curves, the characteristic speeds $v_{\rm gr}^{\rm min}$ and $v_{\rm gr}^{\rm max}$ in particular, is reflected in the temporal evolution and Morlet spectra of impulsively generated wave trains. We conclude that the observed quasi-periodic wave trains not only can be employed to probe such key parameters as density contrasts and profile steepness, but also have the potential to discriminate between the unknown forms of the transverse density distribution.

While standing fast sausage modes in flare loops are often invoked to interpret quasi-periodic pulsations (QPPs) in solar flares, it is unclear as to how they are influenced by the combined effects of a continuous transverse structuring and a finite internal plasma beta ($\beta_{\rm i}$). We derive a generic dispersion relation (DR) governing linear sausage waves in straight magnetic tubes for which plasma pressure is not negligible and the density and temperature inhomogeneities of essentially arbitrary form take place in a layer of arbitrary width. Focusing on fast modes, we find that $\beta_{\rm i}$ only weakly influences $k_{\rm c}$, the critical longitudinal wavenumber separating the leaky from trapped modes. Likewise, for both trapped and leaky modes, the periods $P$ in units of the transverse fast time depend only weakly on $\beta_{\rm i}$, which is compatible with the fact that the effective wavevectors of fast sausage modes are largely perpendicular to the background magnetic field. However, a weak $\beta_{\rm i}$ dependence of the damping times $\tau$ is seen only when the length-to-radius ratio $L/R$ is $\sim 50\%$ larger than some critical value $\pi/(k_{\rm c} R)$, which itself rather sensitively depends on the density contrast, profile steepness as well as on how the transverse structuring is described. In the context of QPPs, we conclude that the much simpler zero-beta theory can be employed for trapped modes, as long as one sees the deduced internal Alfv\'en speed as actually being the fast speed. In contrast, effects due to a finite beta in flare loops should be considered when leaky modes are exploited.

The inverse seesaw mechanism provides an attractive approach to generate small neutrino mass, which origins from a tiny $U(1)_L$ breaking. In this paper, we work in the supersymmetric version of this mechanism, where the singlet-like sneutrino could be an asymmetric dark matter (ADM) candidate in the maximally $U(1)_{L}$ symmetric limit. However, even a tiny $\delta m$, the mass splitting between sneutrino and anti-sneutrino as a result of the tiny $U(1)_{L}$ breaking effect, could lead to fast oscillation between sneutrino and anti-sneutrino and thus spoils the ADM scenario. We study the evolution of this oscillation and find that a weak scale sneutrino, which tolerates a relatively larger $\delta m\sim 10^{-5}$ eV, is strongly favored. We also investigate possible natural ways to realize that small $\delta m$ in the model.

A comprehensive study on the atmospheric neutrino flux in the energy region from sub-GeV up to several TeV using the Super-Kamiokande water Cherenkov detector is presented in this paper. The energy and azimuthal spectra of the atmospheric ${\nu}_e+{\bar{\nu}}_e$ and ${\nu}_{\mu}+{\bar{\nu}}_{\mu}$ fluxes are measured. The energy spectra are obtained using an iterative unfolding method by combining various event topologies with differing energy responses. The azimuthal spectra depending on energy and zenith angle, and their modulation by geomagnetic effects, are also studied. A predicted east-west asymmetry is observed in both the ${\nu}_e$ and ${\nu}_{\mu}$ samples at 8.0 {\sigma} and 6.0 {\sigma} significance, respectively, and an indication that the asymmetry dipole angle changes depending on the zenith angle was seen at the 2.2 {\sigma} level. The measured energy and azimuthal spectra are consistent with the current flux models within the estimated systematic uncertainties. A study of the long-term correlation between the atmospheric neutrino flux and the solar magnetic activity cycle is also performed, and a weak indication of a correlation was seen at the 1.1 {\sigma} level, using SK I-IV data spanning a 20 year period. For particularly strong solar activity periods known as Forbush decreases, no theoretical prediction is available, but a deviation below the typical neutrino event rate is seen at the 2.4 {\sigma} level.

We report the results from a search in Super-Kamiokande for neutrino signals coincident with the first detected gravitational wave events, GW150914 and GW151226, using a neutrino energy range from 3.5 MeV to 100 PeV. We searched for coincident neutrino events within a time window of $\pm$500 seconds around the gravitational wave detection time. Four neutrino candidates are found for GW150914 and no candidates are found for GW151226. The remaining neutrino candidates are consistent with the expected background events. We calculated the 90\% confidence level upper limits on the combined neutrino fluence for both gravitational wave events, which depends on event energy and topologies. Considering the upward going muon data set (1.6 GeV - 100 PeV) the neutrino fluence limit for each gravitational wave event is 14 - 37 (19 - 50) cm$^{-2}$ for muon neutrinos (muon antineutrinos), depending on the zenith angle of the event. In the other data sets, the combined fluence limits for both gravitational wave events range from 2.4$\times 10^{4}$ to 7.0$\times 10^{9}$ cm$^{-2}$.

The DUNE/LBNF program aims to address key questions in neutrino physics and astroparticle physics. Realizing DUNE's potential to reconstruct low-energy particles in the 10-100 MeV energy range will bring significant benefits for all DUNE's science goals. In neutrino physics, low-energy sensitivity will improve neutrino energy reconstruction in the GeV range relevant for the kinematics of DUNE's long-baseline oscillation program. In astroparticle physics, low-energy capabilities will make DUNE's far detectors the world's best apparatus for studying the electron-neutrino flux from a supernova. This will open a new window to unrivaled studies of the dynamics and neutronization of a star's central core in real time, the potential discovery of the neutrino mass hierarchy, provide new sensitivity to physics beyond the Standard Model, and evidence of neutrino quantum-coherence effects. The same capabilities will also provide new sensitivity to `boosted dark matter' models that are not observable in traditional direct dark matter detectors.

Nested multilayer mirrors are commonly used in X-ray telescope structures to increase the collecting area. To balance the difficulty and cost of producing these mirrors, the classical Wolter-I structure has previously been replaced with a conical Wolter-I structure, but it can lead to significantly poorer angular resolution. In this paper, we consider changing one of the mirror shapes (paraboloid or hyperboloid) of the Wolter-I structure to a conical mirror shape, while the other mirror shape remains a quadric surface-type structure, which could thus ensure the imaging quality. Based on a theoretical analysis and the results of simulations, the cone-hyperboloid and paraboloid-cone structures are found to provide the best angular resolutions, and the cone-hyperboloid structure is shorter than the paraboloid-cone structure at the focal length. The cone-hyperboloid structure is then nested to obtain the best on-axis angular resolution and off-axis images.

The Rapid Transient Surveyor (RTS) is a proposed rapid-response, high-cadence adaptive optics (AO) facility for the UH 2.2-m telescope on Maunakea. RTS will uniquely address the need for high-acuity and sensitive near-infrared spectral follow-up observations of tens of thousands of objects in mere months by combining an excellent observing site, unmatched robotic observational efficiency, and an AO system that significantly increases both sensitivity and spatial resolving power. We will initially use RTS to obtain the infrared spectra of ~4,000 Type Ia supernovae identified by the Asteroid Terrestrial-Impact Last Alert System over a two year period that will be crucial to precisely measuring distances and mapping the distribution of dark matter in the z < 0.1 universe. RTS will comprise an upgraded version of the Robo-AO laser AO system and will respond quickly to target-of-opportunity events, minimizing the time between discovery and characterization. RTS will acquire simultaneous-multicolor images with an acuity of 0.07-0.10" across the entire visible spectrum (20% i'-band Strehl in median conditions) and <0.16" in the near infrared, and will detect companions at 0.5" at contrast ratio of ~500. The system will include a high-efficiency prism integral field unit spectrograph: R = 70-140 over a total bandpass of 840-1830 nm with an 8.7" by 6.0" field of view (0.15" spaxels). The AO correction boosts the infrared point-source sensitivity of the spectrograph against the sky background by a factor of seven for faint targets, giving the UH 2.2-m the H-band sensitivity of a 5.7-m telescope without AO.

We measure the angular two-point correlation and angular power spectrum from the NRAO VLA Sky Survey (NVSS) of radio galaxies. They are found to be consistent with the best-fit cosmological model from the Planck analysis, and with the redshift distribution obtained from the Combined EIS-NVSS Survey Of Radio Sources (CENSORS). Our analysis is based on an optimal estimation of the two-point correlation function and makes use of a new mask, that takes into account direction dependent effects of the observations, sidelobe effects of bright sources and galactic foreground. We also set a flux threshold and take the cosmic radio dipole into account. The latter turns out to be an essential step in the analysis. This improved cosmological analysis of the NVSS emphasizes the importance of a flux calibration that is robust and stable on large angular scales for future radio continuum surveys.

Supernova remnants (SNRs) are believed to be the most important acceleration sites for cosmic rays (CRs) below $\sim10^{15}$ eV in the Galaxy. High energy photons, either directly from the shocks of the SNRs or indirectly from the interaction between SNRs and the nearby clouds, are crucial probes for the CR acceleration. Big progresses on observations of SNRs have been achieved by space- and ground-based $\gamma$-ray facilities. However, whether $\gamma$-rays come from accelerated hadrons or not, as well as their connection with the CRs observed at Earth, remains in debate. Large High Altitude Air Shower Observatory (LHAASO), the next generation experiment, is designed to survey the northern part of the very high energy $\gamma$-ray sky from $\sim 0.3$ TeV to PeV with the sensitivity of $\lesssim1\%$ of the Crab nebula flux. In this paper, we indicate that LHAASO will be dedicated to enlarging the $\gamma$-ray SNR samples and improving the spectral and morphological measurements. These measurements, especially at energies above 30 TeV, will be important for us to finally understand the CR acceleration in SNRs.

We report on a measurement of thermal neutrons, generated by the hadronic component of extensive air showers (EAS), by means of a small array of EN-detectors developed for the PRISMA project (PRImary Spectrum Measurement Array), novel devices based on a compound alloy of ZnS(Ag) and $^{6}$LiF. This array has been operated within the ARGO-YBJ experiment at the high altitude Cosmic Ray Observatory in Yangbajing (Tibet, 4300 m a.s.l.). Due to the tight correlation between the air shower hadrons and thermal neutrons, this technique can be envisaged as a simple way to estimate the number of high energy hadrons in EAS. Coincident events generated by primary cosmic rays of energies greater than 100 TeV have been selected and analyzed. The EN-detectors have been used to record simultaneously thermal neutrons and the air shower electromagnetic component. The density distributions of both components and the total number of thermal neutrons have been measured. The correlation of these data with the measurements carried out by ARGO-YBJ confirms the excellent performance of the EN-detector.

Large scale tidal field estimated directly from the distribution of dark matter halos is used to investigate how halo shapes and spin vectors are aligned with the cosmic web. The major, intermediate and minor axes of halos are aligned with the corresponding tidal axes, and halo spin axes tend to be parallel with the intermediate axes and perpendicular to the major axes of tidal field. The strengths of these alignments generally increase with halo mass and redshift, but the dependencies are only through the peak height, {\nu}. The scaling relations of the alignment strengths with the value of {\nu} indicate that the alignment strengths remain roughly constant when the structures within which the halos reside are still in quasi-linear regime, but decreases as nonlinear evolution becomes more important. We also calculate the alignments in projection so that our results can be compared directly with observations. Finally, we investigate the alignments of tidal tensors on large scales, and use the results to understand alignments of halo pairs separated at various distances. Our results suggest coherent structure of the tidal field is the underlying reason for the alignments of halos and galaxies seen in numerical simulations and in observations.

We present a real-time supernova neutrino burst monitor at Super-Kamiokande (SK). Detecting supernova explosions by neutrinos in real time is crucial for giving a clear picture of the explosion mechanism. Since the neutrinos are expected to come earlier than light, a fast broadcasting of the detection may give astronomers a chance to make electromagnetic radiation observations of the explosions right at the onset. The role of the monitor includes a fast announcement of the neutrino burst detection to the world and a determination of the supernova direction. We present the online neutrino burst detection system and studies of the direction determination accuracy based on simulations at SK.

Standing oscillations with multiple periods were found in a number of atmospheric structures on the Sun. The ratio of the period of the fundamental to twice the one of its first overtone, $P_1/2P_2$, is important in applications of solar magneto-seismology. We examine how field-aligned flows impact $P_1/2P_2$ of standing modes in solar magnetic cylinders. For coronal loops, the flow effects are significant for both fast kink and sausage modes. For kink ones, they reduce $P_1/2P_2$ by up to 17\% relative to the static case even when the density contrast between the loop and its surroundings approaches infinity. For sausage modes, the reduction in $P_1/2P_2$ due to flow is typically $\lesssim 5.5\%$ compared with the static case. However, the threshold aspect ratio, only above which can trapped sausage modes be supported, may increase dramatically with the flow magnitude. For photospheric tubes, the flow effect on $P_1/2P_2$ is not as strong. However, when applied to sausage modes, introducing field-aligned flows offers more possibilities in interpreting the multiple periods recently measured. We conclude that field-aligned flows should be taken into account to help better understand what causes the departure of $P_1/2P_2$ from unity.

Standing fast sausage modes in flare loops were suggested to account for a considerable number of quasi-periodic pulsations (QPPs) in the light curves of solar flares. This study continues our investigation into the possibility to invert the measured periods $P$ and damping times $\tau$ of sausage modes to deduce the transverse Alfv\'en time $R/v_{\rm Ai}$, density contrast $\rho_{\rm i}/\rho_{\rm e}$, and the steepness of the density distribution transverse to flare loops. A generic dispersion relation (DR) governing linear sausage modes is derived for pressureless cylinders where density inhomogeneity of arbitrary form takes place within the cylinder. We show that in general the inversion problem is under-determined for QPP events where only a single sausage mode exists, be the measurements spatially resolved or unresolved. While $R/v_{\rm Ai}$ can be inferred to some extent, the range of possible steepness parameters may be too broad to be useful. However, for spatially resolved measurements where an additional mode is present, it is possible to deduce self-consistently $\rho_{\rm i}/\rho_{\rm e}$, the profile steepness, and the internal Alfv\'en speed $v_{\rm Ai}$. We show that at least for a recent QPP event that involves a fundamental kink mode in addition to a sausage one, flare loop parameters are well constrained, even if the specific form of the transverse density distribution remains unknown. We conclude that spatially resolved, multi-mode QPP measurements need to be pursued for inferring flare loop parameters.

We report on the extensive multi-wavelength observations of the blazar Markarian 421 (Mrk 421) covering radio to gamma-rays, during the 4.5 year period of ARGO-YBJ and Fermi common operation time, from August 2008 to February 2013. In particular, thanks to the ARGO-YBJ and Fermi data, the whole energy range from 100 MeV to 10 TeV is covered without any gap. In the observation period, Mrk 421 showed both low and high activity states at all wavebands. The correlations among flux variations in different wavebands were analyzed. Seven large flares, including five X-ray flares and two GeV gamma-ray flares with variable durations (3-58 days), and one X-ray outburst phase were identified and used to investigate the variation of the spectral energy distribution with respect to a relative quiescent phase. During the outburst phase and the seven flaring episodes, the peak energy in X-rays is observed to increase from sub-keV to few keV. The TeV gamma-ray flux increases up to 0.9-7.2 times the flux of the Crab Nebula. The behavior of GeV gamma-rays is found to vary depending on the flare, a feature that leads us to classify flares into three groups according to the GeV flux variation. Finally, the one-zone synchrotron self-Compton model was adopted to describe the emission spectra. Two out of three groups can be satisfactorily described using injected electrons with a power-law spectral index around 2.2, as expected from relativistic diffuse shock acceleration, whereas the remaining group requires a harder injected spectrum. The underlying physical mechanisms responsible for different groups may be related to the acceleration process or to the environment properties.

We examine the influence of a continuous density structuring transverse to coronal slabs on the dispersive properties of fundamental standing kink and sausage modes supported therein. We derive generic dispersion relations (DRs) governing linear fast waves in pressureless straight slabs with general transverse density distributions, and focus on the cases where the density inhomogeneity takes place in a layer of arbitrary width and in arbitrary form. The physical relevance of the solutions to the DRs is demonstrated by the corresponding time-dependent computations. For all profiles examined, the lowest-order kink modes are trapped regardless of longitudinal wavenumber $k$. A continuous density distribution introduces a difference to their periods of $\lesssim 13\%$ when $k$ is the observed range, relative to the case where the density profile takes a step-function form. Sausage modes and other branches of kink modes are leaky at small $k$, and their periods and damping times are heavily influenced by how the transverse density profile is prescribed, the lengthscale in particular. These modes have sufficiently high quality to be observable only for physical parameters representative of flare loops. We conclude that while the simpler DR pertinent to a step-function profile can be used for the lowest-order kink modes, the detailed information on the transverse density structuring needs to be incorporated into studies of sausage modes and higher-order kink modes.

As one of the most active blazars, Mrk421 is an excellent candidate for the study of the physical processes within the jets of AGN. Here we report on the extensive multi-wavelength observations of Mrk 421 over 4.5 years, from 2008 August to 2013 February. This source was simultaneously monitored by several experiments at different wavelengths: ARGO-YBJ in TeV $\gamma$-rays, $Fermi$-LAT in GeV $\gamma$-rays, $Swift$-BAT in hard X-rays, $RXTE$-ASM, $MAXI$ and $Swift$-XRT in soft X-rays, $Swift$-UVOT in ultraviolet, and OVRO in radio frequencies. In particular, thanks to the ARGO-YBJ and $Fermi$ data, the whole energy range from 100 MeV to 10 TeV is covered without any gap. According to the observed light curves, ten states (including seven large flares, two quiescent phases and one outburst) were selected. For the first time, the multi-wavelength spectral evolutions of Mrk 421 during different states were systematically analyzed. During the outburst phase and the seven flaring episodes, the peak energy in X-rays is observed to increase from sub-keV to few keV. The TeV $\gamma$-ray flux increases up to 0.9$-$7.2 times the flux of the Crab Nebula. The behavior of GeV $\gamma$-rays is found to vary depending on the flare, a feature that leads us to classify flares into three groups according to the GeV flux variation. Finally, the observed radiation spectra above 0.3 keV of different states can be reasonably described by a simple one-zone synchrotron self-Compton model. The underlying physical mechanisms responsible for different states may be related to the acceleration process or to variations of the ambient medium.

In this work a series of methods are developed for understanding the Friedmann equation when it is beyond the reach of the Chebyshev theorem. First it will be demonstrated that every solution of the Friedmann equation admits a representation as a roulette such that information on the latter may be used to obtain that for the former. Next the Friedmann equation is integrated for a quadratic equation of state and for the Randall--Sundrum II universe, leading to a harvest of a rich collection of new interesting phenomena. Finally an analytic method is used to isolate the asymptotic behavior of the solutions of the Friedmann equation, when the equation of state is of an extended form which renders the integration impossible, and to establish a universal exponential growth law.

Standing sausage modes in flare loops are important for interpreting quasi-periodic pulsations (QPPs) in solar flare lightcurves. We propose an inversion scheme that consistently uses their periods $P$ and damping times $\tau$ to diagnose flare loop parameters. We derive a generic dispersion relation governing linear sausage waves in pressure-less straight tubes, for which the transverse density inhomogeneity takes place in a layer of arbitrary width $l$ and is of arbitrary form. We find that $P$ and $\tau$ depend on the combination of $[R/v_{\rm Ai}, L/R, l/R, \rho_{\rm i}/\rho_{\rm e}]$, where $R$ is the loop radius, $L$ is the looplength, $v_{\rm Ai}$ is the internal Alfv\'en speed, and $\rho_{\rm i}/\rho_{\rm e}$ is the density contrast. For all the density profiles examined, $P$ and $\tau$ experience saturation when $L/R \gg 1$, yielding an inversion curve in the $[R/v_{\rm Ai}, l/R, \rho_{\rm i}/\rho_{\rm e}]$ space with a specific density profile when $L/R$ is sufficiently large. When applied to a spatially unresolved QPP event, the scheme yields that $R/v_{\rm Ai}$ is the best constrained, whereas $l/R$ corresponds to the other extreme. For spatially resolved QPPs, while $L/R \gg 1$ cannot be assumed beforehand, an inversion curve remains possible due to additional geometrical constraints. When a spatially resolved QPP event involves another mode, as is the case for a recent event, the full set of $[v_{\rm Ai}, l, \rho_{\rm i}/\rho_{\rm e}]$ can be inferred. We conclude that the proposed scheme provides a useful tool for magneto-seismologically exploiting QPPs.

Sausage modes are important in coronal seismology. Spatially damped propagating sausage waves were recently observed in the solar atmosphere. We examine how wave leakage influences the spatial damping of sausage waves propagating along coronal structures modeled by a cylindrical density enhancement embedded in a uniform magnetic field. Working in the framework of cold magnetohydrodynamics, we solve the dispersion relation (DR) governing sausage waves for complex-valued longitudinal wavenumber $k$ at given real angular frequencies $\omega$. For validation purposes, we also provide analytical approximations to the DR in the low-frequency limit and in the vicinity of $\omega_{\rm c}$, the critical angular frequency separating trapped from leaky waves. In contrast to the standing case, propagating sausage waves are allowed for $\omega$ much lower than $\omega_{\rm c}$. However, while able to direct their energy upwards, these low-frequency waves are subject to substantial spatial attenuation. The spatial damping length shows little dependence on the density contrast between the cylinder and its surroundings, and depends only weakly on frequency. This spatial damping length is of the order of the cylinder radius for $\omega \lesssim 1.5 v_{\rm Ai}/a$, where $a$ and $v_{\rm Ai}$ are the cylinder radius and the Alfv\'en speed in the cylinder, respectively. We conclude that if a coronal cylinder is perturbed by symmetric boundary drivers (e.g., granular motions) with a broadband spectrum, wave leakage efficiently filters out the low-frequency components.

The events recorded by ARGO-YBJ in more than five years of data collection have been analyzed to determine the diffuse gamma-ray emission in the Galactic plane at Galactic longitudes 25{\deg} < l < 100{\deg} and Galactic latitudes . The energy range covered by this analysis, from ~350 GeV to ~2 TeV, allows the connection of the region explored by Fermi with the multi-TeV measurements carried out by Milagro. Our analysis has been focused on two selected regions of the Galactic plane, i.e., 40{\deg} < l < 100{\deg} and 65{\deg} < l < 85{\deg} (the Cygnus region), where Milagro observed an excess with respect to the predictions of current models. Great care has been taken in order to mask the most intense gamma-ray sources, including the TeV counterpart of the Cygnus cocoon recently identified by ARGO-YBJ, and to remove residual contributions. The ARGO-YBJ results do not show any excess at sub-TeV energies corresponding to the excess found by Milagro, and are consistent with the predictions of the Fermi model for the diffuse Galactic emission. From the measured energy distribution we derive spectral indices and the differential flux at 1 TeV of the diffuse gamma-ray emission in the sky regions investigated.

Standing sausage modes are important in interpreting quasi-periodic pulsations in the lightcurves of solar flares. Their periods and damping times play an important role in seismologically diagnosing key parameters like the magnetic field strength in regions where flare energy is released. Usually such applications are based on theoretical results neglecting unresolved fine structures in magnetized loops. However, the existence of fine structuring is suggested on both theoretical and observational grounds. Adopting the framework of cold magnetohydrodynamics (MHD), we model coronal loops as magnetized cylinders with a transverse equilibrium density profile comprising a monolithic part and a modulation due to fine structuring in the form of concentric shells. The equation governing the transverse velocity perturbation is solved with an initial-value-problem approach, and the effects of fine structuring on the periods $P$ and damping times $\tau$ of global, leaky, standing sausage modes are examined. A parameter study shows that fine structuring, be it periodically or randomly distributed, brings changes of only a few percent to $P$ and $\tau$ when there are more than about ten shells. The monolithic part, its steepness in particular, plays a far more important role in determining $P$ and $\tau$. We conclude that when measured values of $P$ and $\tau$ of sausage modes are used for seismological purposes, it is justified to use theoretical results where the effects due to fine structuring are neglected.

Providing an early warning of galactic supernova explosions from neutrino signals is important in studying supernova dynamics and neutrino physics. A dedicated supernova trigger system has been designed and installed in the data acquisition system at Daya Bay and integrated into the worldwide Supernova Early Warning System (SNEWS). Daya Bay's unique feature of eight identically-designed detectors deployed in three separate experimental halls makes the trigger system naturally robust against cosmogenic backgrounds, enabling a prompt analysis of online triggers and a tight control of the false-alert rate. The trigger system is estimated to be fully sensitive to 1987A-type supernova bursts throughout most of the Milky Way. The significant gain in sensitivity of the eight-detector configuration over a mass-equivalent single detector is also estimated. The experience of this online trigger system is applicable to future projects with spatially distributed detectors.

The ARGO-YBJ experiment has been in stable data taking from November 2007 till February 2013 at the YangBaJing Cosmic Ray Observatory (4300 m a.s.l.). The detector consists of a single layer of Resistive Plate Chambers (RPCs) ( about 6700 m^2}) operated in streamer mode. The signal pick-up is obtained by means of strips facing one side of the gas volume. The digital readout of the signals, while allows a high space-time resolution in the shower front reconstruction, limits the measurable energy to a few hundred TeV. In order to fully investigate the 1-10 PeV region, an analog readout has been implemented by instrumenting each RPC with two large size electrodes facing the other side of the gas volume. Since December 2009 the RPC charge readout has been in operation on the entire central carpet (about 5800 m^2). In this configuration the detector is able to measure the particle density at the core position where it ranges from tens to many thousands of particles per m^2. Thus ARGO-YBJ provides a highly detailed image of the charge component at the core of air showers. In this paper we describe the analog readout of RPCs in ARGO-YBJ and discuss both the performance of the system and the physical impact on the EAS measurements.

The ARGO-YBJ experiment is a full-coverage air shower detector located at the Yangbajing Cosmic Ray Observatory (Tibet, People's Republic of China, 4300 m a.s.l.). The high altitude, combined with the full-coverage technique, allows the detection of extensive air showers in a wide energy range and offer the possibility of measuring the cosmic ray proton plus helium spectrum down to the TeV region, where direct balloon/space-borne measurements are available. The detector has been in stable data taking in its full configuration from November 2007 to February 2013. In this paper the measurement of the cosmic ray proton plus helium energy spectrum is presented in the region 3-300 TeV by analyzing the full collected data sample. The resulting spectral index is $\gamma = -2.64 \pm 0.01$. These results demonstrate the possibility of performing an accurate measurement of the spectrum of light elements with a ground based air shower detector.

This paper is a continuation of our earlier study on the integrability of the Friedmann equations in the light of the Chebyshev theorem. Our main focus will be on a series of important, yet not previously touched, problems when the equation of state for the perfect-fluid universe is nonlinear. These include the generalized Chaplygin gas, two-term energy density, trinomial Friedmann, Born--Infeld, and two-fluid models. We show that some of these may be integrated using Chebyshev's result while other are out of reach by the theorem but may be integrated explicitly by other methods. With the explicit integration, we are able to understand exactly the roles of the physical parameters in various models play in the cosmological evolution. For example, in the Chaplygin gas universe, it is seen that, as far as there is a tiny presence of nonlinear matter, linear matter makes contribution to the dark matter, which becomes significant near the phantom divide line. The Friedmann equations also arise in areas of physics not directly related to cosmology. We provide some examples ranging from geometric optics and central orbits to soap films and the shape of glaciated valleys to which our results may be applied.

The search for Gamma Ray Burst (GRB) emission in the energy range 1-100 GeV in coincidence with the satellite detection has been carried out using the Astrophysical Radiation with Ground-based Observatory at YangBaJing (ARGO-YBJ) experiment. The high altitude location (4300 m a.s.l.), the large active surface ($\sim$ 6700 m$^2$ of Resistive Plate Chambers), the wide field of view ($\sim 2~$sr, limited only by the atmospheric absorption) and the high duty cycle ($>$ 86 %) make the ARGO-YBJ experiment particularly suitable to detect short and unexpected events like GRBs. With the scaler mode technique, i.e., counting all the particles hitting the detector with no measurement of the primary energy and arrival direction, the minimum threshold of $\sim$ 1 GeV can be reached, overlapping the direct measurements carried out by satellites. During the experiment lifetime, from December 17, 2004 to February 7, 2013, a total of 206 GRBs occurring within the ARGO-YBJ field of view (zenith angle $\theta$ $\le$ 45$^{\circ}$) have been analyzed. This is the largest sample of GRBs investigated with a ground-based detector. Two lightcurve models have been assumed and since in both cases no significant excess has been found, the corresponding fluence upper limits in the 1-100 GeV energy region have been derived, with values as low as 10$^{-5}~$erg cm$^{-2}$. The analysis of a subset of 24 GRBs with known redshift has been used to constrain the fluence extrapolation to the GeV region together with possible cutoffs under different assumptions on the spectrum.

The measurement of cosmic ray energy spectra, in particular for individual species, is an essential approach in finding their origin. Locating the "knees" of the spectra is an important part of the approach and has yet to be achieved. Here we report a measurement of the mixed Hydrogen and Helium spectrum using the combination of the ARGO-YBJ experiment and of a prototype Cherenkov telescope for the LHAASO experiment. A knee feature at 640+/-87 TeV, with a clear steepening of the spectrum, is observed. This gives fundamental inputs to galactic cosmic ray acceleration models.

We investigate the differential number counts of sources in radio continuum surveys, including all terms at linear order in cosmological perturbations. Our framework does not assume a specific gauge condition. This general approach allows us to recover gauge invariance explicitly. With the complete derivations of the covariant volume integral on the past light cone, we have identified several contributions in the number counts. To clarify their underlying physics, we present each contributions in terms of scalar, vector and tensor modes. This theoretical framework promises to be widely applicable to continuum radio galaxy surveys to model the expected angular power spectrum and two-point correlation.

Luminous distant quasars are unique probes of the high redshift intergalactic medium (IGM) and of the growth of massive galaxies and black holes in the early universe. Absorption due to neutral Hydrogen in the IGM makes quasars beyond a redshift of z~6.5 very faint in the optical $z$-band, thus locating quasars at higher redshifts require large surveys that are sensitive above 1 micron. We report the discovery of three new z>6.5 quasars, corresponding to an age of the universe of <850 Myr, selected as z-band dropouts in the Pan-STARRS1 survey. This increases the number of known z>6.5 quasars from 4 to 7. The quasars have redshifts of z=6.50, 6.52, and 6.66, and include the brightest z-dropout quasar reported to date, PSO J036.5078+03.0498 with M_1450=-27.4. We obtained near-infrared spectroscopy for the quasars and from the MgII line we estimate that the central black holes have masses between 5x10^8 and 4x10^9 M_sun, and are accreting close to the Eddington limit (L_Bol/L_Edd=0.13-1.2). We investigate the ionized regions around the quasars and find near zone radii of R_NZ=1.5-5.2 proper Mpc, confirming the trend of decreasing near zone sizes with increasing redshift found for quasars at 5.7<z<6.4. By combining R_NZ of the PS1 quasars with those of 5.7<z<7.1 quasars in the literature, we derive a luminosity corrected redshift evolution of R_NZ,corrected=(7.2+/-0.2)-(6.1+/-0.7)x(z-6) Mpc. However, the large spread in R_NZ in the new quasars implies a wide range in quasar ages and/or a large variation in the neutral Hydrogen fraction along different lines of sight.

Continuum and HI surveys with the Square Kilometre Array (SKA) will allow us to probe some of the most fundamental assumptions of modern cosmology, including the Cosmological Principle. SKA all-sky surveys will map an enormous slice of space-time and reveal cosmology at superhorizon scales and redshifts of order unity. We illustrate the potential of these surveys and discuss the prospects to measure the cosmic radio dipole at high fidelity. We outline several potentially transformational tests of cosmology to be carried out by means of SKA all-sky surveys.

PSR B1821$-$24 is a solitary millisecond pulsar (MSP) which radiates multi-wavelength pulsed photons. It has complex radio, X-ray and $\gamma$-ray pulse profiles with distinct peak phase-separations that challenge the traditional caustic emission models. Using the single-pole annular gap model with suitable magnetic inclination angle ($\alpha=40^\circ$) and viewing angle ($\zeta=75^\circ$), we managed to reproduce its pulse profiles of three wavebands. It is found that the middle radio peak is originated from the core gap region at high altitudes, and the other two radio peaks are originated from the annular gap region at relatively low altitudes. Two peaks of both X-ray and $\gamma$-ray wavebands are fundamentally originated from annular gap region, while the $\gamma$-ray emission generated from the core gap region contributes somewhat to the first $\gamma$-ray peak. Precisely reproducing the multi-wavelength pulse profiles of PSR B1821$-$24 enables us to understand emission regions of distinct wavebands and justify pulsar emission models.

This short but systematic work demonstrates a link between Chebyshev's theorem and the explicit integration in cosmological time $t$ and conformal time $\eta$ of the Friedmann equations in all dimensions and with an arbitrary cosmological constant $\Lambda$. More precisely, it is shown that for spatially flat universes an explicit integration in $t$ may always be carried out, and that, in the non-flat situation and when $\Lambda$ is zero and the ratio $w$ of the pressure and energy density in the barotropic equation of state of the perfect-fluid universe is rational, an explicit integration may be carried out if and only if the dimension $n$ of space and $w$ obey some specific relations among an infinite family. The situation for explicit integration in $\eta$ is complementary to that in $t$. More precisely, it is shown in the flat-universe case with $\Lambda\neq0$ that an explicit integration in $\eta$ can be carried out if and only if $w$ and $n$ obey similar relations among a well-defined family which we specify, and that, when $\Lambda=0$, an explicit integration can always be carried out whether the space is flat, closed, or open. We also show that our method may be used to study more realistic cosmological situations when the equation of state is nonlinear.

We present an overview of the design of IRIS, an infrared (0.84 - 2.4 micron) integral field spectrograph and imaging camera for the Thirty Meter Telescope (TMT). With extremely low wavefront error (<30 nm) and on-board wavefront sensors, IRIS will take advantage of the high angular resolution of the narrow field infrared adaptive optics system (NFIRAOS) to dissect the sky at the diffraction limit of the 30-meter aperture. With a primary spectral resolution of 4000 and spatial sampling starting at 4 milliarcseconds, the instrument will create an unparalleled ability to explore high redshift galaxies, the Galactic center, star forming regions and virtually any astrophysical object. This paper summarizes the entire design and basic capabilities. Among the design innovations is the combination of lenslet and slicer integral field units, new 4Kx4k detectors, extremely precise atmospheric dispersion correction, infrared wavefront sensors, and a very large vacuum cryogenic system.

Maximizing the grating efficiency is a key goal for the first light instrument IRIS (Infrared Imaging Spectrograph) currently being designed to sample the diffraction limit of the TMT (Thirty Meter Telescope). Volume Phase Holographic (VPH) gratings have been shown to offer extremely high efficiencies that approach 100% for high line frequencies (i.e., 600 to 6000l/mm), which has been applicable for astronomical optical spectrographs. However, VPH gratings have been less exploited in the near-infrared, particularly for gratings that have lower line frequencies. Given their potential to offer high throughputs and low scattered light, VPH gratings are being explored for IRIS as a potential dispersing element in the spectrograph. Our team has procured near-infrared gratings from two separate vendors. We have two gratings with the specifications needed for IRIS current design: 1.51-1.82{\mu}m (H-band) to produce a spectral resolution of 4000 and 1.19- 1.37 {\mu}m (J-band) to produce a spectral resolution of 8000. The center wavelengths for each grating are 1.629{\mu}m and 1.27{\mu}m, and the groove densities are 177l/mm and 440l/mm for H-band R=4000 and J-band R=8000, respectively. We directly measure the efficiencies in the lab and find that the peak efficiencies of these two types of gratings are quite good with a peak efficiency of ~88% at the Bragg angle in both TM and TE modes at H-band, and 90.23% in TM mode, 79.91% in TE mode at J-band for the best vendor. We determine the drop in efficiency off the Bragg angle, with a 20-23% decrease in efficiency at H-band when 2.5 degree deviation from the Bragg angle, and 25%-28% decrease at J-band when 5{\deg} deviation from the Bragg angle.

We present the efficiency of near-infrared reflective ruled diffraction gratings designed for the InfraRed Imaging Spectrograph (IRIS). IRIS is a first light, integral field spectrograph and imager for the Thirty Meter Telescope (TMT) and narrow field infrared adaptive optics system (NFIRAOS). We present our experimental setup and analysis of the efficiency of selected reflective diffraction gratings. These measurements are used as a comparison sample against selected candidate Volume Phase Holographic (VPH) gratings (see Chen et al., this conference). We investigate the efficiencies of five ruled gratings designed for IRIS from two separate vendors. Three of the gratings accept a bandpass of 1.19-1.37 {\mu}m (J band) with ideal spectral resolutions of R=4000 and R=8000, groove densities of 249 and 516 lines/mm, and blaze angles of 9.86 and 20.54 degrees, respectively. The other two gratings accept a bandpass of 1.51-1.82 {\mu}m (H Band) with an ideal spectral resolution of R=4000, groove density of 141 lines/mm, and blaze angle of 9.86{\deg}. We measure the efficiencies off blaze angle for all gratings and the efficiencies between the polarization transverse magnetic (TM) and transverse electric (TE) states. The peak reflective efficiencies are 98.90 +/- 3.36% (TM) and 84.99 +/- 2.74% (TM) for the H-band R=4000 and J-band R=4000 respectively. The peak reflective efficiency for the J-band R=8000 grating is 78.78 +/- 2.54% (TE). We find that these ruled gratings do not exhibit a wide dependency on incident angle within +/-3{\deg}. Our best-manufactured gratings were found to exhibit a dependency on the polarization state of the incident beam with a ~10-20% deviation, consistent with the theoretical efficiency predictions.

The extended TeV gamma-ray source ARGO J2031+4157 (or MGRO J2031+41) is positionally consistent with the Cygnus Cocoon discovered by $Fermi$-LAT at GeV energies in the Cygnus superbubble. Reanalyzing the ARGO-YBJ data collected from November 2007 to January 2013, the angular extension and energy spectrum of ARGO J2031+4157 are evaluated. After subtracting the contribution of the overlapping TeV sources, the ARGO-YBJ excess map is fitted with a two-dimensional Gaussian function in a square region of $10^{\circ}\times 10^{\circ}$, finding a source extension $\sigma_{ext}$= 1$^{\circ}$.8$\pm$0$^{\circ}$.5. The observed differential energy spectrum is $dN/dE =(2.5\pm0.4) \times 10^{-11}(E/1 TeV)^{-2.6\pm0.3}$ photons cm$^{-2}$ s$^{-1}$ TeV$^{-1}$, in the energy range 0.2-10 TeV. The angular extension is consistent with that of the Cygnus Cocoon as measured by $Fermi$-LAT, and the spectrum also shows a good connection with the one measured in the 1-100 GeV energy range. These features suggest to identify ARGO J2031+4157 as the counterpart of the Cygnus Cocoon at TeV energies. The Cygnus Cocoon, located in the star-forming region of Cygnus X, is interpreted as a cocoon of freshly accelerated cosmic rays related to the Cygnus superbubble. The spectral similarity with Supernova Remnants indicates that the particle acceleration inside a superbubble is similar to that in a SNR. The spectral measurements from 1 GeV to 10 TeV allows for the first time to determine the possible spectrum slope of the underlying particle distribution. A hadronic model is adopted to explain the spectral energy distribution.

Magnetohydrodynamic waves are important for diagnosing the physical parameters of coronal plasmas. Field-aligned flows appear frequently in coronal loops.We examine the effects of transverse density and plasma flow structuring on standing sausage modes trapped in coronal loops, and examine their observational implications. We model coronal loops as straight cold cylinders with plasma flow embedded in a static corona. An eigen-value problem governing propagating sausage waves is formulated, its solutions used to construct standing modes. Two transverse profiles are distinguished, one being the generalized Epstein distribution (profile E) and the other (N) proposed recently in Nakariakov et al.(2012). A parameter study is performed on the dependence of the maximum period $P_\mathrm{max}$ and cutoff length-to-radius ratio $(L/a)_{\mathrm{cutoff}}$ in the trapped regime on the density parameters ($\rho_0/\rho_\infty$ and profile steepness $p$) and flow parameters (magnitude $U_0$ and profile steepness $u$). For either profile, introducing a flow reduces $P_\mathrm{max}$ relative to the static case. $P_\mathrm{max}$ depends sensitively on $p$ for profile N but is insensitive to $p$ for profile E. By far the most important effect a flow introduces is to reduce the capability for loops to trap standing sausage modes: $(L/a)_{\mathrm{cutoff}}$ may be substantially reduced in the case with flow relative to the static one. If the density distribution can be described by profile N, then measuring the sausage mode period can help deduce the density profile steepness. However, this practice is not feasible if profile E better describes the density distribution. Furthermore, even field-aligned flows with magnitudes substantially smaller than the ambient Alfv\'en speed can make coronal loops considerably less likely to support trapped standing sausage modes.

MGRO J2019+37 within the Cygnus region is a bright and extended source revealed by Milagro at 12-35 TeV. This source is almost as bright as Crab Nebula in northern sky, while it is not confirmed by ARGO-YBJ around TeV. Up to now, no obvious counterpart at low energy wavelengths has been found. Hence, MGRO J2019+37 becomes mysterious and its VHE gamma-ray emission mechanism is attractive. In this paper, a brief summary of the multi-wavelength observations from Radio to gamma-ray is presented. All the available data from the XMM-Newton and INTEGRAL at X-ray, and the Fermi-LAT at gamma-ray bands were used to get constraint on its emission flux at low energy wavelengths. Then, its possible counterparts and the VHE emission mechanism are discussed.

The $^{18}$Ne($\alpha$,$p$)$^{21}$Na reaction is thought to be one of the key breakout reactions from the hot CNO cycles to the rp-process in type I x-ray bursts. In this work, the resonant properties of the compound nucleus $^{22}$Mg have been investigated by measuring the resonant elastic scattering of $^{21}$Na+$p$. An 89 MeV $^{21}$Na radioactive beam delivered from the CNS Radioactive Ion Beam Separator bombarded an 8.8 mg/cm$^2$ thick polyethylene (CH$_{2}$)$_{n}$ target. The $^{21}$Na beam intensity was about 2$\times$10$^{5}$ pps, with a purity of about 70% on target. The recoiled protons were measured at the center-of-mass scattering angles of $\theta_{c.m.}$$\approx$175.2${^\circ}$, 152.2${^\circ}$, and 150.5${^\circ}$ by three sets of $\Delta E$-$E$ telescopes, respectively. The excitation function was obtained with the thick-target method over energies $E_x$($^{22}$Mg)=5.5--9.2 MeV. In total, 23 states above the proton-threshold in $^{22}$Mg were observed, and their resonant parameters were determined via an $R$-matrix analysis of the excitation functions. We have made several new $J^{\pi}$ assignments and confirmed some tentative assignments made in previous work. The thermonuclear $^{18}$Ne($\alpha$,$p$)$^{21}$Na rate has been recalculated based on our recommended spin-parity assignments. The astrophysical impact of our new rate has been investigated through one-zone postprocessing x-ray burst calculations. We find that the $^{18}$Ne($\alpha$,$p$)$^{21}$Na rate significantly affects the peak nuclear energy generation rate, reaction fluxes, as well as the onset temperature of this breakout reaction in these astrophysical phenomena.

The energy spectrum of cosmic Hydrogen and Helium nuclei has been measured, below the so-called "knee", by using a hybrid experiment with a wide field-of-view Cherenkov telescope and the Resistive Plate Chamber (RPC) array of the ARGO-YBJ experiment at 4300 m above sea level. The Hydrogen and Helium nuclei have been well separated from other cosmic ray components by using a multi-parameter technique. A highly uniform energy resolution of about 25% is achieved throughout the whole energy range (100 TeV - 700 TeV). The observed energy spectrum is compatible with a single power law with index gamma=-2.63+/-0.06.

Procedures and results on hardware level detector calibration in Super-Kamiokande (SK) are presented in this paper. In particular, we report improvements made in our calibration methods for the experimental phase IV in which new readout electronics have been operating since 2008. The topics are separated into two parts. The first part describes the determination of constants needed to interpret the digitized output of our electronics so that we can obtain physical numbers such as photon counts and their arrival times for each photomultiplier tube (PMT). In this context, we developed an in-situ procedure to determine high-voltage settings for PMTs in large detectors like SK, as well as a new method for measuring PMT quantum efficiency and gain in such a detector. The second part describes the modeling of the detector in our Monte Carlo simulation, including in particular the optical properties of its water target and their variability over time. Detailed studies on the water quality are also presented. As a result of this work, we achieved a precision sufficient for physics analysis over a wide energy range (from a few MeV to above a TeV). For example, the charge determination was understood at the 1% level, and the timing resolution was 2.1 nsec at the one-photoelectron charge level and 0.5 nsec at the 100-photoelectron charge level.

The ARGO-YBJ detector is an extensive air shower array that has been used to monitor the northern $\gamma$-ray sky at energies above 0.3 TeV from 2007 November to 2013 January. In this paper, we present the results of a sky survey in the declination band from $-10^{\circ}$ to $70^{\circ}$, using data recorded over the past five years. With an integrated sensitivity ranging from 0.24 to $\sim$1 Crab units depending on the declination, six sources have been detected with a statistical significance greater than 5 standard deviations. Several excesses are also reported as potential $\gamma$-ray emitters. The features of each source are presented and discussed. Additionally, $95\%$ confidence level upper limits of the flux from the investigated sky region are shown. Specific upper limits for 663 GeV $\gamma$-ray AGNs inside the ARGO-YBJ field of view are reported. The effect of the absorption of $\gamma$-rays due to the interaction with extragalactic background light is estimated.

{\it Fermi}-LAT has detected pulsed gamma-ray emissions with high confidences from more than 40 millisecond pulsars (MSPs). Here we study the phase-averaged gamma-ray properties of MSPs by using revised version of a self-consistent outer gap model. In this model, a strong multipole magnetic field near the stellar surface for a MSP is assumed and such a field will be close to the surface magnetic fields ($\sim 10^{11}- 10^{12}$ G) of young ulsars; the outer gap of a MSP is controlled by photon-photon pair production process, where the effects of magnetic inclination angle ($\alpha$) and magnetic geometry have been taken into account, therefore the fractional size of the outer gap is a function of not only pulsar's period and magnetic field strength but also magnetic inclination angle and radial distance to the neutron star, the inner boundary of the outer gap can be estimated by the pair production process of the gamma-ray photons which are produced by the back-flowing particles through the null charge surface; inside the outer gap, a Gaussian distribution of the parallel electric field along the trans-field thickness is assumed, and the gamma-ray emission is represented by the emission from the average radial distance along the central field lines of the outer gap. Using this model, the phase-averaged gamma-ray spectra are calculated and compared with the observed spectra of 37 MSPs given by the second {\it Fermi}-LAT catalog of gamma-ray pulsars, our results show that the {\it Fermi}-LAT results can be well explained by this model. The thermal X-ray emission properties from MSPs are also investigated.

Fundamental standing modes and their overtones play an important role in coronal seismology. We examine how a significant field-aligned flow affects standing modes supported by coronal loops, modeled here as cold magnetic slabs. Of particular interest are the period ratios of the fundamental to its $(n-1)$-th overtone ($P_1/nP_n$) for both kink and sausage modes, and the threshold half-width-to-length ratio for sausage modes. For standing kink modes, the flow significantly reduces $P_1/nP_n$ in general, the effect being particularly strong for larger $n$ and when the density contrast $\rho_0/\rho_e$ between loops and their surroundings is weak. That said, even when $\rho_0/\rho_e$ approaches infinity, this effect is still substantial, reducing the minimal $P_1/nP_n$ by up to 13.7% (24.5%) for $n=2$ ($n=4$) relative to the static case, when the Alfv\'en Mach number $M_A$ reaches 0.8 where $M_A$ measures the loop flow speed in units of the internal Alfv\'en speed. For standing sausage modes, though not negligible, the flow effect in reducing $P_1/nP_n$ is not as strong. However, the threshold half-width-to-length ratio is considerably larger in the flowing case than its static counterpart. For $\rho_0/\rho_e$ in the range $[9, 1024]$ and $M_A$ in the range $[0, 0.5]$, an exhaustive parameter study yields that this threshold is well fitted by our Equation (23) which involves the two parameters in a simple way. This allows one to analytically constrain the combination $(\rho_0/\rho_e, M_A)$ for a loop with known width-to-length ratio when a standing sausage oscillation is identified therein. It also allows one to further examine the idea of partial sausage modes, and the flow is found to reduce significantly the spatial extent where partial modes are allowed.

The $^{18}$Ne($\alpha$,$p$)$^{21}$Na reaction provides a pathway for breakout from the hot CNO cycles to the $rp$-process in type I x-ray bursts. To better determine this astrophysical reaction rate, the resonance parameters of the compound nucleus $^{22}$Mg have been investigated by measuring the resonant elastic scattering of $^{21}$Na+$p$. An 89 MeV $^{21}$Na radioactive ion beam was produced at the CNS Radioactive Ion Beam Separator and bombarded an 8.8 mg/cm$^2$ thick polyethylene target. The recoiled protons were measured at scattering angles of $\theta_{c.m.}$$\approx 175 {^\circ}$ and 152${^\circ}$ by three $\Delta E$-$E$ silicon telescopes. The excitation function was obtained with a thick-target method over energies $E_x$($^{22}$Mg)=5.5--9.2 MeV. The resonance parameters have been determined through an $R$-matrix analysis. For the first time, the $J^{\pi}$ values for ten states above the alpha threshold in $^{22}$Mg have been experimentally determined in a single consistent measurement. We have made three new $J^{\pi}$ assignments and confirmed seven of the ten tentative assignments in the previous work. The $^{18}$Ne($\alpha$,$p$)$^{21}$Na reaction rate has been recalculated, and the astrophysical impact of our new rate has been investigated through one-zone postprocessing x-ray burst calculations. We find that the $^{18}$Ne($\alpha$,$p$)$^{21}$Na rate significantly affects the peak nuclear energy generation rate and the onset temperature of this breakout reaction in these phenomena.

The ground-based EAS array is usually operated with a high duty cycle (> 90%) and a large field of view (~2 sr), which can continuously monitor the sky. It is essential and irreplaceable to understand the gamma-ray emission mechanism and intrinsic physics progress of the variable source AGN. The EAS arrays, AS-r experiment (since 1990) and ARGO-YBJ experiment (since 2007), have continuously monitored the northern sky at energies above 3 TeV and 0.3 TeV, respectively. They have made substantial contributions for long-term monitoring of Mrk 421 and Mrk 501. In this paper, we will review the results obtained by the EAS arrays. The next generation of EAS array, LHAASO project, will boost the sensitivity of current EAS array at least up to 30 times with a much wider energy range from 40 GeV to 1 PeV. Beside increasing the number of VHE gamma-ray sources, it will guide us look sight into the properties of jet, and throw light on the determining of the EBL, intergalactic magnetic fields, and the validity of the Lorentz Invariance.

We have tested the effect of spatial gradients in stellar mass-to-light ratio (Y) on measurements of black hole masses (MBH) derived from stellar orbit superposition models. Such models construct a static gravitational potential for a galaxy and its central black hole, but typically assume spatially uniform Y. We have modeled three giant elliptical galaxies with gradients alpha = d(log Y)/d(log r) from -0.2 to +0.1. Color and line strength gradients suggest mildly negative alpha in these galaxies. Introducing a negative (positive) gradient in Y increases (decreases) the enclosed stellar mass near the center of the galaxy and leads to systematically smaller (larger) MBH measurements. For models with alpha = -0.2, the best-fit values of MBH are 28%, 27%, and 17% lower than the constant-Y case, in NGC 3842, NGC 6086, and NGC 7768, respectively. For alpha = +0.1, MBH are 14%, 22%, and 17% higher than the constant-Y case for the three respective galaxies. For NGC 3842 and NGC 6086, this bias is comparable to the statistical errors from individual modeling trials. At larger radii, negative (positive) gradients in Y cause the total stellar mass to decrease (increase) and the dark matter fraction within one effective radius to increase (decrease).

We report the observation of a very high energy \gamma-ray source, whose position is coincident with HESS J1841-055. This source has been observed for 4.5 years by the ARGO-YBJ experiment from November 2007 to July 2012. Its emission is detected with a statistical significance of 5.3 standard deviations. Parameterizing the source shape with a two-dimensional Gaussian function we estimate an extension \sigma=(0.40(+0.32,-0.22}) degree, consistent with the HESS measurement. The observed energy spectrum is dN/dE =(9.0-+1.6) x 10^{-13}(E/5 TeV)^{-2.32-+0.23} photons cm^{-2} s^{-1} TeV^{-1}, in the energy range 0.9-50 TeV. The integral \gamma-ray flux above 1 TeV is 1.3-+0.4 Crab units, which is 3.2-+1.0 times the flux derived by HESS. The differences in the flux determination between HESS and ARGO-YBJ, and possible counterparts at other wavelengths are discussed.

We present the discovery and timing solutions of five new pulsars by students involved in the Pulsar Search Collaboratory (PSC), a NSF-funded joint program between the National Radio Astronomy Observatory and West Virginia University designed to excite and engage high-school students in Science, Technology, Engineering, and Mathematics (STEM) and related fields. We encourage students to pursue STEM fields by apprenticing them within a professional scientific community doing cutting edge research, specifically by teaching them to search for pulsars. The students are analyzing 300 hours of drift-scan survey data taken with the Green Bank Telescope at 350 MHz. These data cover 2876 square degrees of the sky. Over the course of five years, more than 700 students have inspected diagnostic plots through a web-based graphical interface designed for this project. The five pulsars discovered in the data have spin periods ranging from 3.1 ms to 4.8 s. Among the new discoveries are - PSR J1926-1314, a long period, nulling pulsar; PSR J1821+0155, an isolated, partially recycled 33-ms pulsar; and PSR J1400-1438, a millisecond pulsar in a 9.5-day orbit whose companion is likely a white dwarf star.

As one of the brightest active blazars in both X-ray and very high energy $\gamma$-ray bands, Mrk 501 is very useful for physics associated with jets from AGNs. The ARGO-YBJ experiment is monitoring it for $\gamma$-rays above 0.3 TeV since November 2007. Starting from October 2011 the largest flare since 2005 is observed, which lasts to about April 2012. In this paper, a detailed analysis is reported. During the brightest $\gamma$-rays flaring episodes from October 17 to November 22, 2011, an excess of the event rate over 6 $\sigma$ is detected by ARGO-YBJ in the direction of Mrk 501, corresponding to an increase of the $\gamma$-ray flux above 1 TeV by a factor of 6.6$\pm$2.2 from its steady emission. In particular, the $\gamma$-ray flux above 8 TeV is detected with a significance better than 4 $\sigma$. Based on time-dependent synchrotron self-Compton (SSC) processes, the broad-band energy spectrum is interpreted as the emission from an electron energy distribution parameterized with a single power-law function with an exponential cutoff at its high energy end. The average spectral energy distribution for the steady emission is well described by this simple one-zone SSC model. However, the detection of $\gamma$-rays above 8 TeV during the flare challenges this model due to the hardness of the spectra. Correlations between X-rays and $\gamma$-rays are also investigated.

We study the strong gravitational lensing on the equatorial plane of a quasi-Kerr compact object with arbitrary quadrupole moments which can be used to model the super-massive central object of the galaxy. We find that, when the quadrupolar correction parameter $\xi$ takes the positive (negative) value, the photon-sphere radius $r_{ps}$, the minimum impact parameter $u_{ps}$, the coefficient $\bar{b}$, the relative magnitudes $r_m$ and the angular position of the relativistic images $\theta_{\infty}$ are larger (smaller) than the results obtained in the Kerr black hole, but the coefficient $\bar{a}$, the deflection angle $\alpha(\theta)$ and the angular separation $s$ are smaller (larger) than that in the Kerr black hole. These features may offer a way to probe special properties for some rotating compact objects by the astronomical instruments in the future.

We present the result of an indirect search for high energy neutrinos from WIMP annihilation in the Sun using upward-going muon (upmu) events at Super-Kamiokande. Datasets from SKI-SKIII (3109.6 days) were used for the analysis. We looked for an excess of neutrino signal from the Sun as compared with the expected atmospheric neutrino background in three upmu categories: stopping, non-showering, and showering. No significant excess was observed. The 90% C.L. upper limits of upward-going muon flux induced by WIMPs of 100 GeV/c$^2$ were 6.4$\times10^{-15}$ cm$^{-2}$ sec$^{-1}$ and 4.0$\times10^{-15}$ cm$^{-2}$ sec$^{-1}$ for the soft and hard annihilation channels, respectively. These limits correspond to upper limits of 4.5$\times10^{-39}$ cm$^{-2}$ and 2.7$\times10^{-40}$ cm$^{-2}$ for spin-dependent WIMP-nucleon scattering cross sections in the soft and hard annihilation channels, respectively.

We study a model of a scalar field minimally coupled to gravity, with a specific potential energy for the scalar field, and include curvature and radiation as two additional parameters. Our goal is to obtain analytically the complete set of configurations of a homogeneous and isotropic universe as a function of time. This leads to a geodesically complete description of the universe, including the passage through the cosmological singularities, at the classical level. We give all the solutions analytically without any restrictions on the parameter space of the model or initial values of the fields. We find that for generic solutions the universe goes through a singular (zero-size) bounce by entering a period of antigravity at each big crunch and exiting from it at the following big bang. This happens cyclically again and again without violating the null energy condition. There is a special subset of geodesically complete non-generic solutions which perform zero-size bounces without ever entering the antigravity regime in all cycles. For these, initial values of the fields are synchronized and quantized but the parameters of the model are not restricted. There is also a subset of spatial curvature-induced solutions that have finite-size bounces in the gravity regime and never enter the antigravity phase. These exist only within a small continuous domain of parameter space without fine tuning initial conditions. To obtain these results, we identified 25 regions of a 6-parameter space in which the complete set of analytic solutions are explicitly obtained.

The properties of the ergosphere and energy extraction by Penrose process in a rotating non-Kerr black hole are investigated. It is shown that the ergosphere is sensitive to the deformation parameter $\epsilon$ and the shape of the ergosphere becomes thick with increase of the parameter $\epsilon$. It is of interest to note that, comparing with the Kerr black hole, the deformation parameter $\epsilon$ can enhance the maximum efficiency of the energy extraction process greatly. Especially, for the case of $a>M$, the non-Kerr metric describes a superspinning compact object and the maximum efficiency can exceed 60%, while it is only 20.7% for the extremal Kerr black hole.

We study the strong gravitational lensing in the background of a rotating non-Kerr compact object with a deformed parameter $\epsilon$ and an unbound rotation parameter $a$. We find that the photon sphere radius and the deflection angle depend sharply on the parameters $\epsilon$ and $a$. For the case in which the black hole is more prolate than a Kerr black hole, the photon sphere exists only in the regime $\epsilon\leq\epsilon_{max}$ for prograde photon. The upper limit $\epsilon_{max}$ is a function of the rotation parameter $a$. As $\epsilon>\epsilon_{max}$, the deflection angle of the light ray closing very to the naked singularity is a positive finite value, which is different from those in both the usual Kerr black hole spacetime and in the rotating naked singularity described by Janis-Newman-Winicour metric. For the oblate black hole and the retrograde photon, there does not exist such a threshold value. Modelling the supermassive central object of the Galaxy as a rotating non-Kerr compact object, we estimated the numerical values of the coefficients and observables for gravitational lensing in the strong field limit.

We study the accretion process in the thin disk around a rotating non-Kerr black hole with a deformed parameter and an unbound rotation parameter. Our results show that the presence of the deformed parameter $\epsilon$ modifies the standard properties of the disk. For the case in which the black hole is more oblate than a Kerr black hole, the larger deviation leads to the smaller energy flux, the lower radiation temperature and the fainter spectra luminosity in the disk. For the black hole with positive deformed parameter, we find that the effect of the deformed parameter on the disk becomes more complicated. It depends not only on the rotation direction of the black hole and the orbit particles, but also on the sign of the difference between the deformed parameter $\epsilon$ and a certain critical value $\epsilon_{c}$. These significant features in the mass accretion process may provide a possibility to test the no-hair theorem in the strong field regime in future astronomical observations.

GUT monopoles captured by the Sun's gravitation are expected to catalyze proton decays via the Callan-Rubakov process. In this scenario, protons, which initially decay into pions, will ultimately produce \nu_{e}, \nu_{\mu} and \bar{\nu}_{\mu}. After undergoing neutrino oscillation, all neutrino species appear when they arrive at the Earth, and can be detected by a 50,000 metric ton water Cherenkov detector, Super-Kamiokande (SK). A search for low energy neutrinos in the electron total energy range from 19 to 55 MeV was carried out with SK and gives a monopole flux limit of F_M(\sigma_0/1 mb) < 6.3 \times 10^{-24} (\beta_M/10^{-3})^2 cm^{-2} s^{-1} sr^{-1} at 90% C.L., where \beta_M is the monopole velocity in units of the speed of light and \sigma_0 is the catalysis cross section at \beta_M=1. The obtained limit is more than eight orders of magnitude more stringent than the current best cosmic-ray supermassive monopole flux limit, F_M < 1 \times 10^{-15} cm^{-2} s^{-1} sr^{-1} for \beta_M < 10^{-3} and also two orders of magnitude lower than the result of the Kamiokande experiment, which used a similar detection method.

Parity violation found in the Cosmic Microwave Background (CMB) radiation is a crucial clue for the non-standard cosmological model or the possible contamination of various foreground residuals and/or calibration of the CMB data sets. In this paper, we study the directional properties of the CMB parity asymmetry by excluding the $m=0$ modes in the definition of parity parameters. We find that the preferred directions of the parity parameters coincide with the CMB kinematic dipole, which implies that the CMB parity asymmetry may be connected with the possible contamination of the residual dipole component. We also find that such tendency is not only localized at $l=2,3$, but in the extended multipole ranges up to $l\sim 22$.

We report the observation of TeV gamma-rays from the Cygnus region using the ARGO-YBJ data collected from 2007 November to 2011 August. Several TeV sources are located in this region including the two bright extended MGRO J2019+37 and MGRO J2031+41. According to the Milagro data set, at 20 TeV MGRO J2019+37 is the most significant source apart from the Crab Nebula. No signal from MGRO J2019+37 is detected by the ARGO-YBJ experiment, and the derived flux upper limits at 90% confidence level for all the events above 600 GeV with medium energy of 3 TeV are lower than the Milagro flux, implying that the source might be variable and hard to be identified as a pulsar wind nebula. The only statistically significant (6.4 standard deviations) gamma-ray signal is found from MGRO J2031+41, with a flux consistent with the measurement by Milagro.

A wide field of view Cherenkov/fluorescence telescope array is one of the main components of the Large High Altitude Air Shower Observatory project. To serve as Cherenkov and fluorescence detectors, a flexible and mobile design is adopted for easy reconfiguring of the telescope array. Two prototype telescopes have been constructed and successfully run at the site of the ARGO-YBJ experiment in Tibet. The features and performance of the telescopes are presented.

The ARGO-YBJ experiment is an extensive air shower array with full coverage RPC detectors located at Yangbajing (4300 m asl, Tibet, China). It is operated with high duty cycle (>86%) and a large field of view ($\sim$ 2sr). It continuously monitors the entire overhead sky at $\gamma$-ray energies above 0.1 TeV. In the talk, we will present the result of the northern sky survey (between declinations of -10$^{\circ}$ and 70$^{\circ}$) from an analysis of ~4 years of the ARGO-YBJ data (between July 2006 and February 2011). There are four known TeV sources observed with significance greater than 5 S.D.. The significance from Crab Nebula is more than 16 S.D.. 90% confidence level upper limits to the flux from all directions in the sky are also presented, which vary from 0.09 to 0.53 Crab unit for Crab-like point sources.

We show that the f(T) gravitational paradigm, in which gravity is described by an arbitrary function of the torsion scalar, can provide a mechanism for realizing bouncing cosmologies, thereby avoiding the Big Bang singularity. After constructing the simplest version of an f(T) matter bounce, we investigate the scalar and tensor modes of cosmological perturbations. Our results show that metric perturbations in the scalar sector lead to a background-dependent sound speed, which is a distinguishable feature from Einstein gravity. Additionally, we obtain a scale-invariant primordial power spectrum, which is consistent with cosmological observations, but suffers from the problem of a large tensor-to-scalar ratio. However, this can be avoided by introducing extra fields, such as a matter bounce curvaton.

We have studied the quasinormal modes of a massive scalar field coupling to Einstein's tensor in the spacelike stretched $AdS_3$ black hole spacetime. We find that both the right-moving and left-moving quasinormal frequencies depend not only on the warped parameter $v$ of black hole, but also on the coupling between the scalar field and Einstein's tensor. Moreover, we also discuss the warped $AdS/CFT$ correspondence from the quasinormal modes and probe the effects of the coupling on the left and right conformal weights $h_L$ and $h_R$ of the operators dual to the scalar field in the boundary.

We study the accretion process in the thin disk around a squashed Kaluza-Klein black hole and probe the effects of the extra dimensional scale $\rho_0$ on the physical properties of the disk. Our results show that with the increase of the parameter $\rho_0$, the energy flux, the conversion efficiency, the radiation temperature, the spectra luminosity and the spectra cut-off frequency of the thin accretion disk decrease, but the inner border of the disk increases. This implies that the extra dimension scale imprints in the mass accretion process in the disk.

We study the collision of two particles with the different rest masses moving in the equatorial plane of a Kerr-Taub-NUT spacetime and get the center-of-mass (CM) energy for the particles. We find that the CM energy depends not only on the rotation parameter, $a$, but also on the NUT charge of the Kerr-Taub-NUT spacetime, $n$. Especially, for the extremal Kerr-Taub-NUT spacetime, an unlimited CM energy can be approached if the parameter $a$ is in the range $[1,\sqrt{2}]$, which is different from that of the Kerr and Kerr-Newman black holes.

We explore the influence of new physics on extremely high energy cosmic ray (EXCER) particles. In particular, we devote our mind to one example of new physics, unparticle stuff, on one specific process that EXECRs participate in, photopion production of the EXECR nucleon with 2.7K cosmic microwave background radiation (CMBR). Through computing the differential cross section of virtual exchange of unparticle in the $p+\gamma\to p+ \pi^0$ process, we acquire the general consequence of new physics on the EXECR propagation. It is astonishing but reasonable that due to the lowness of interaction energy $\sqrt{s}$, the new physics will play a nearly negligible role in EXECR interaction.

ARGO-YBJ is an air shower detector array with a fully covered layer of resistive plate chambers. It is operated with a high duty cycle and a large field of view. It continuously monitors the northern sky at energies above 0.3 TeV. In this paper, we report a long-term monitoring of Mrk 421 over the period from 2007 November to 2010 February. This source was observed by the satellite-borne experiments Rossi X-ray Timing Explorer and Swift in the X-ray band. Mrk 421 was especially active in the first half of 2008. Many flares are observed in both X-ray and gamma-ray bands simultaneously. The gamma-ray flux observed by ARGO-YBJ has a clear correlation with the X-ray flux. No lag between the X-ray and gamma-ray photons longer than 1 day is found. The evolution of the spectral energy distribution is investigated by measuring spectral indices at four different flux levels. Hardening of the spectra is observed in both X-ray and gamma-ray bands. The gamma-ray flux increases quadratically with the simultaneously measured X-ray flux. All these observational results strongly favor the synchrotron self-Compton process as the underlying radiative mechanism.

We investigate the strong gravitational lensing in a squashed Kaluza-Klein black hole immersed in the G\"{o}del universe with global rotation. Our result show that the strong gravitational lensing in the squashed Kaluza-Klein G\"{o}del black hole spacetime has some distinct behaviors from that in the Kerr case. In the squashed Kaluza-Klein G\"{o}del black hole spacetime, the photon sphere radius, the minimum impact parameter, the coefficient $\bar{a}$, $\bar{b}$ and the deflection angle $\alpha(\theta)$ in the $\phi$ direction are independent of whether the photon goes with or against the global rotation of the G\"{o}del Universe. While in the Kerr black hole, the values of these quantities for the prograde photons are different from those for the retrograde photons. Moreover, the coefficient of $\bar{b}$ increases with $j$ in the squashed Kaluza-Klein G\"{o}del black hole, but decreases with $a$ in the Kerr case. We also probe the influence of the squashed effect on the strong gravitational lensing in this black hole and find that in the extremely squashed case $\rho_0=0$, the coefficient $\bar{a}$ is a constant 1 and is independent of the global rotation of the G\"{o}del Universe. Furthermore, we assume that the gravitational field of the supermassive central object of the Galaxy can be described by this metric and estimate the numerical values of the coefficients and the main observables in the strong gravitational lensing.

We present analytic solutions to a class of cosmological models described by a canonical scalar field minimally coupled to gravity and experiencing self interactions through a hyperbolic potential. Using models and methods inspired by 2T-physics, we show how analytic solutions can be obtained in flat/open/closed Friedmann-Robertson-Walker universes. Among the analytic solutions, there are many interesting geodesically complete cyclic solutions in which the universe bounces at either zero or finite sizes. When geodesic completeness is imposed, it restricts models and their parameters to a certain parameter subspace, including some quantization conditions on initial conditions in the case of zero-size bounces, but no conditions on initial conditions for the case of finite-size bounces. We will explain the theoretical origin of our model from the point of view of 2T-gravity as well as from the point of view of the colliding branes scenario in the context of M-theory. We will indicate how to associate solutions of the quantum Wheeler-deWitt equation with our classical analytic solutions, mention some physical aspects of the cyclic solutions, and outline future directions.

We study collision of two general geodesic particles around the Kerr-Newman (KN) black hole and get the center-of-mass (CM) energy of the non-marginally and marginally bound critical particles in the direct collision and LSO collision scenarios. We find the constraint conditions that arbitrarily high CM energy can be obtained for the near-horizon collision of two general geodesic particles in the extremal KN black hole, and note that the charge decreases the value of the latitude in which arbitrarily high CM energy can occurs. We also interpret why the high-velocity collision belt centers to the equator with the increase of the charge.

Noncommutative geometry may be a starting point to a quantum gravity. We study the influence of the spacetime noncommutative parameter on the strong field gravitational lensing in the noncommutative Schwarzschild black-hole spacetime and obtain the angular position and magnification of the relativistic images. Supposing that the gravitational field of the supermassive central object of the galaxy described by this metric, we estimate the numerical values of the coefficients and observables for strong gravitational lensing. Comparing to the Reissner-Norstr\"{om} black hole, we find that the influences of the spacetime noncommutative parameter is similar to those of the charge, just these influences are much smaller. This may offer a way to distinguish a noncommutative black hole from a Reissner-Norstr\"{om} black hole, and may probe the spacetime noncommutative constant $\vartheta$ \cite{foot} by the astronomical instruments in the future.

The sun blocks cosmic ray particles from outside the solar system, forming a detectable shadow in the sky map of cosmic rays detected by the ARGO-YBJ experiment in Tibet. Because the cosmic ray particles are positive charged, the magnetic field between the sun and the earth deflects them from straight trajectories and results in a shift of the shadow from the true location of the sun. Here we show that the shift measures the intensity of the field which is transported by the solar wind from the sun to the earth.

A new approach is given for the implementation of boundary conditions used in solving the Mukhanov-Sasaki equation in the context of inflation. The familiar quantization procedure is reviewed, along with a discussion of where one might expect deviations from the standard approach to arise. The proposed method introduces a (model dependent) fitting function for the z"/z and a"/a terms in the Mukhanov-Sasaki equation for scalar and tensor modes, as well as imposes the boundary conditions at a finite conformal time. As an example, we employ a fitting function, and compute the spectral index, along with its running, for a specific inflationary model which possesses background equations that are analytically solvable. The observational upper bound on the tensor to scalar ratio is used to constrain the parameters of the boundary conditions in the tensor sector as well. An overview on the generalization of this method is also discussed.

We study the second-order phase transition (SOPT) for the spherically symmetric Kehagias-Sfetsos (KS) black hole in the deformed H\v{o}rava-Lifshitz gravity by applying the methods of equilibrium and non-equilibrium fluctuations. We find that, although the KS black hole has only one mass parameter as the usual Schwarzschild ones, the SOPT will take place if the mass of the KS black hole changes across the critical point $\frac{\sqrt{5+\sqrt{33}}(\sqrt{33}-1)} {16\sqrt{\omega}} $. The result show us that there is difference between the H\v{o}rava-Lifshitz gravity and the Einstein's gravity theory.

We investigate the cosmological perturbations in f(T) gravity. Examining the pure gravitational perturbations in the scalar sector using a diagonal vierbien, we extract the corresponding dispersion relation, which provides a constraint on the f(T) ansatzes that lead to a theory free of instabilities. Additionally, upon inclusion of the matter perturbations, we derive the fully perturbed equations of motion, and we study the growth of matter overdensities. We show that f(T) gravity with f(T) constant coincides with General Relativity, both at the background as well as at the first-order perturbation level. Applying our formalism to the power-law model we find that on large subhorizon scales (O(100 Mpc) or larger), the evolution of matter overdensity will differ from LCDM cosmology. Finally, examining the linear perturbations of the vector and tensor sectors, we find that (for the standard choice of vierbein) f(T) gravity is free of massive gravitons.

The Large High Altitude Air Shower Observatory project is proposed to study high energy gamma ray astronomy ( 40 GeV-1 PeV ) and cosmic ray physics ( 20 TeV-1 EeV ). The wide field of view Cherenkov telescope array, as a component of the LHAASO project, will be used to study energy spectrum and compositions of cosmic ray by measuring the total Cherenkov light generated by air showers and shower maximum depth. Two prototype telescopes have been in operation since 2008. The pointing accuracy of each telescope is crucial to the direction reconstruction of the primary particles. On the other hand the primary energy reconstruction relies on the shape of the Cherenkov image on the camera and the unrecorded photons due to the imperfect connections between photomultiplier tubes. UV bright stars are used as point-like objects to calibrate the pointing and to study the optical properties of the camera, the spot size and the fractions of unrecorded photons in the insensitive areas of the camera.

In a recent paper, four of the present authors proposed a class of dark matter models where generalized parity symmetry leads to equality of dark matter abundance with baryon asymmetry of the Universe and predicts dark matter mass to be around 5 GeV. In this note we explore how this model can be tested in direct search experiments. In particular, we point out that if the dark matter happens to be the mirror neutron, the direct detection cross section has the unique feature that it increases at low recoil energy unlike the case of conventional WIMPs. It is also interesting to note that the predicted spin-dependent scattering could make significant contribution to the total direct detection rate, especially for light nucleus. With this scenario, one could explain recent DAMA and CoGeNT results.

The Spectral and Photometric Imaging Receiver (SPIRE), is the Herschel Space Observatory`s submillimetre camera and spectrometer. It contains a three-band imaging photometer operating at 250, 350 and 500 microns, and an imaging Fourier Transform Spectrometer (FTS) which covers simultaneously its whole operating range of 194-671 microns (447-1550 GHz). The SPIRE detectors are arrays of feedhorn-coupled bolometers cooled to 0.3 K. The photometer has a field of view of 4' x 8', observed simultaneously in the three spectral bands. Its main operating mode is scan-mapping, whereby the field of view is scanned across the sky to achieve full spatial sampling and to cover large areas if desired. The spectrometer has an approximately circular field of view with a diameter of 2.6'. The spectral resolution can be adjusted between 1.2 and 25 GHz by changing the stroke length of the FTS scan mirror. Its main operating mode involves a fixed telescope pointing with multiple scans of the FTS mirror to acquire spectral data. For extended source measurements, multiple position offsets are implemented by means of an internal beam steering mirror to achieve the desired spatial sampling and by rastering of the telescope pointing to map areas larger than the field of view. The SPIRE instrument consists of a cold focal plane unit located inside the Herschel cryostat and warm electronics units, located on the spacecraft Service Module, for instrument control and data handling. Science data are transmitted to Earth with no on-board data compression, and processed by automatic pipelines to produce calibrated science products. The in-flight performance of the instrument matches or exceeds predictions based on pre-launch testing and modelling: the photometer sensitivity is comparable to or slightly better than estimated pre-launch, and the spectrometer sensitivity is also better by a factor of 1.5-2.

We analyze large scale mapping observations of the molecular lines in the 12CO (J=2-1), 12CO (J=3-2), 13CO (J=2-1), and 13CO (J=3-2) transition emissions toward the Cepheus B molecular cloud with the KOSMA 3m-telescope. The integrated intensity map of the 12CO (J=2-1) transition has shown a structure with a compact core and a compact ridge extended in the north-west of the core. The cloud is surrounded by an optically bright rim, where the radiation-driven implosion (RDI) may greatly change the gas properties. The intensities of the CO (J=3-2) transition are higher than those of the CO (J=2-1) transition along the rim area.We find characteristic RDI structure in positionvelocity diagrams. Non-LTE Large velocity gradient (LVG) model analysis shows that the density and temperature at the edge are higher than that in the center. Our results provide evidences that RDI is taking place in Cepheus B molecular cloud.

We propose a model of asymmetric dark matter (DM) where the dark sector is an identical copy of both forces and matter of the standard model (SM) as in the mirror universe models discussed in literature. In addition to being connected by gravity, the SM and DM sectors are also connected at high temperature by a common set of heavy right-handed Majorana neutrinos via their Yukawa couplings to leptons and Higgs bosons. The lightest nucleon in the dark (mirror) sector is a candidate for dark matter. The out of equilibrium decay of right-handed neutrino produces equal lepton asymmetry in both sectors via resonant leptogenesis which then get converted to baryonic and dark baryonic matter. The dark baryon asymmetry due to higher dark nucleon masses leads to higher dark matter density compared to the familiar baryon density that is observed. The standard model neutrinos in this case acquire masses from the inverse seesaw mechanism. A kinetic mixing between the U(1) gauge fields of the two sectors is introduced to guarantee the success of Big-Bang Nucleosynthesis.

Dark matter is a focused issue in galactic evolution and cosmology. Axion is a viable particle candidate for dark matter. Its interaction with photon is an effective way to detect it, e.g., pseudoscalar-photon interaction will generate vacuum dichroism in a magnetic field. Motivated to measure the QED vacuum birefringence and to detect pseudoscalar-photon interaction, we started to build up the Q & A experiment (QED [Quantum Electrodynamics] and Axion experiment) in 1994. In this talk, we first give a brief historical account of planet hunting and dark matter evidence. We then review our 3.5 m Fabry-Perot interferometer together with our results of measuring vacuum dichroism and gaseous Cotton-Mouton effects. Our first results give (-0.2 $\pm$ 2.8) $\times$ 10$^{-13}$ rad/pass, at 2.3 T with 18,700 passes through a 0.6 m long magnet, for vacuum dichroism measurement. We are upgrading our interferometer to 7 m arm-length with a new 1.8 m 2.3 T permanent magnet capable of rotation up to 13 cycles per second. We aim at 10 nrad/Hz$^{1/2}$ optical sensitivity with 532 nm cavity finesse around 100,000. When achieved, vacuum dichroism would be measured to 8.6 $\times$ 10$^{-17}$ rad/pass in about 50 days, and QED birefringence would be measured to 28 %.

We study the effects of unparticle physics on the cosmic ray photon and $e^\pm$, including on the pair production (PP) and elastic scattering (ES) of cosmic ray photon off various background radiations, and on the inverse Compton scattering of cosmic ray $e^\pm$ with cosmic radiations. We compute the spin-averaged amplitudes squared of three processes and find that the advent of unparticle will never significantly change the interactions of cosmic ray photon and $e^\pm$ with various background radiations, although the available papers show that ES which occurs in the tree-level through unparticle exchanges will easily surpass PP in the approximate parameter regions.

We propose a new class of R-parity violating extension of MSSM with type II seesaw mechanism for neutrino masses where an unstable gravitino is the dark matter of the Universe. It decays predominantly into three leptons final states, thereby providing a natural explanation of the positron excess but no antiproton excess in the PAMELA experiment. The model can explain neutrino masses without invoking any high scale physics while keeping the pre-existing baryon asymmetry of the universe in tact.

A new method is proposed to compute the bulk viscosity in strange quark matter at high densities. Using the method it is straightforward to prove that the bulk viscosity is positive definite, which is not so easy to accomplish in other approaches especially for multi-component fluids like strange quark matter with light up and down quarks and massive strange quarks.

Unparticle $\U$ with scaling dimension $d_\U$ has peculiar thermal properties due to its unique phase space structure. We find that the equation of state parameter $\omega_\U$, the ratio of pressure to energy density, is given by $1/(2d_\U +1)$ providing a new form of energy in our universe. In an expanding universe, the unparticle energy density $\rho_\U(T)$ evolves dramatically differently from that for photons. For $d_\U >1$, even if $\rho_\U(T_D)$ at a high decoupling temperature $T_D$ is very small, it is possible to have a large relic density $\rho_\U(T^0_\gamma)$ at present photon temperature $T^0_\gamma$, large enough to play the role of dark matter. We calculate $T_D$ and $\rho_\U(T^0_\gamma)$ using photon-unparticle interactions for illustration.

We present the results from arcsecond resolution observations of various line transitions at 1.3 mm toward hypercompact HII region G28.20-0.04N. With the SMA data, we have detected and mapped the transitions in the CH$_{3}$CN, CO, $^{13}$CO, SO$_{2}$, OCS, and CH$_{3}$OH molecular lines as well as the radio recombination line H30$\alpha$. The observations and analysis indicate a hot core associated with G28.20-0.04N. The outflow and possible rotation are detected in this region.

There are obvious discrepancies among various experimental constraints on the variation of the fine-structure constant, $\alpha$. We attempt to discuss the issue in the framework of de Sitter invariant Special Relativity (${\cal SR}_{c,R}$) and to present a possible solution to the disagreement. In addition, on the basis of the observational data and the discussions presented in this Letter, we derive a rough theoretical estimate of the radius of the Universe.

A number of experiments are underway to detect vacuum birefringence and dichroism -- PVLAS, Q & A, and BMV. Recently, PVLAS experiment has observed optical rotation in vacuum by a magnetic field (vacuum dichroism). Theoretical interpretations of this result include a possible pseudoscalar-photon interaction and the existence of millicharged fermions. Here, we report the progress and first results of Q & A (QED [quantum electrodynamics] and Axion) experiment proposed and started in 1994. A 3.5-m high-finesse (around 30,000) Fabry-Perot prototype detector extendable to 7-m has been built and tested. We use X-pendulums and automatic control schemes developed by the gravitational-wave detection community for mirror suspension and cavity control. To polarize the vacuum, we use a 2.3-T dipole permanent magnet, with 27-mm-diameter clear borehole and 0.6-m field length,. In the experiment, the magnet is rotated at 5-10 rev/s to generate time-dependent polarization signal with twice the rotation frequency. Our ellipsometer/polarization-rotation-detection-system is formed by a pair of Glan-Taylor type polarizing prisms with extinction ratio lower than 10-8 together with a polarization modulating Faraday Cell with/without a quarter wave plate. We made an independent calibration of our apparatus by performing a measurement of gaseous Cotton-Mouton effect of nitrogen. We present our first experimental results and give a brief discussion of our experimental limit on pseudo-scalar-photon interaction and millicharged fermions.

Dust in the host galaxies of gamma-ray bursts (GRBs) dims and reddens their afterglow spectra. Knowledge of the nature of this dust is crucial for correcting for extinction, providing clues to the nature of GRB progenitors, and probing the interstellar medium of high-redshift galaxies as well as the nature of cosmic dust when the universe was much younger and galaxies were much less evolved. The dust and extinction properties of GRB host galaxies are still poorly known. Unlike previous work, we derive in this Letter the extinction curves for 10 GRB host galaxies without a priori assumption of any specific extinction types (such as that of the Milky Way, or the Small/Large Magellanic Clouds). It is found that there appears to exist two different types of extinction curves: one is relatively flat and gray, the other displays a steeper dependence on inverse wavelength, closely resembling that of the Milky Way but with the 2175$\Angstrom$ feature removed.

We discuss the threshold anomaly in ultra-high energy cosmic ray physics by assuming that the matter world just be perturbation of a de Sitter universe, which is consistent with the recent astronomical observations: about two thirds of the whole energy in the universe is contributed by a small positive cosmological constant. One-particle states are presented explicitly. It is noticed that the dispersion relation of free particles is dependent on the degrees of freedom of angular momentum. This fact can be regarded as the effects of the cosmological constant on kinematics of particles.

We show that a tiny but non-zero positive cosmological constant, which has been strongly suggested by the recent astronomical observations on supernovae and CMBR, may change notably the behaviors of the ultrahigh energy cosmic ray interacting with soft photons. The threshold anomaly of the ultrahigh energy cosmic ray disappears naturally after the effects of the cosmological constant are taken into account.

We present a covariant framework of kinematics in the dS spacetime, which is a natural postulate of recent astronomical observations ($\Lambda>0)$. One-particle states are presented explicitly. It is noticed that the dispersion relation of free particles is dependent on the degrees of freedom of angular momentum and spin. This fact can be referred to as the effects of the cosmological constant on kinematics of particles. The kinematics in dS spacetime is used to investigate the phenomenon of ultra high energy cosmic rays. We emphasize the possibility of solving the threshold anomalies of the interactions between ultra high energy cosmic rays and soft photons in the covariant framework of kinematics in the dS spacetime.