We present the results of the search for ultra-high-energy photons with nine years of data from the Telescope Array surface detector. A multivariate classifier is built upon 16 reconstructed parameters of the extensive air shower. These parameters are related to the curvature and the width of the shower front, the steepness of the lateral distribution function, and the timing parameters of the waveforms sensitive to the shower muon content. A total number of two photon candidates found in the search is fully compatible with the expected background. The $95\%\,$CL limits on the diffuse flux of the photons with energies greater than $10^{18.0}$, $10^{18.5}$, $10^{19.0}$, $10^{19.5}$ and $10^{20.0}$ eV are set at the level of $0.067$, $0.012$, $0.0036$, $0.0013$, $0.0013~\mbox{km}^{-2}\mbox{yr}^{-1}\mbox{sr}^{-1}$ correspondingly.

The results on ultra-high-energy cosmic rays (UHECR) mass composition obtained with the Telescope Array surface detector are presented. The analysis employs the boosted decision tree (BDT) multivariate analysis built upon 14 observables related to both the properties of the shower front and the lateral distribution function. The multivariate classifier is trained with Monte-Carlo sets of events induced by the primary protons and iron. An average atomic mass of UHECR is presented for energies $10^{18.0}-10^{20.0}\ \mbox{eV}$. The average atomic mass of primary particles shows no significant energy dependence and corresponds to $\langle \ln A \rangle = 2.0 \pm 0.1 (stat.) \pm 0.44 (syst.)$. The result is compared to the mass composition obtained by the Telescope Array with $\mbox{X}_{\mbox{max}}$ technique along with the results of other experiments. Possible systematic errors of the method are discussed.

An analysis is given of inflation based on a supersymmetric Dirac-Born-Infeld (DBI) action in an axionic landscape. The DBI model we discuss involves a landscape of chiral superfields with one $U(1)$ shift symmetry which is broken by instanton type non-perturbative terms in the superpotential. Breaking of the shift symmetry leads to one pseudo-Nambu-Goldstone-boson which acts as the inflaton while the remaining normalized phases of the chiral fields generically labeled axions are invariant under the $U(1)$ shift symmetry. The analysis is carried out in the vacuum with stabilized saxions, which are the magnitudes of the chiral fields. Regions of the parameter space where slow-roll inflation occurs are exhibited and the spectral indices as well as the ratio of the tensor to the scalar power spectrum are computed. An interesting aspect of supersymmetric DBI models analyzed is that in most of the parameter space tensor to scalar ratio and scalar spectral index are consistent with Planck data if slow roll occurs and is not eternal. Also interesting is that the ratio of the tensor to the scalar power spectrum can be large and can lie close to the experimental upper limit and thus testable in improved experiment. Non-Gaussianity in this class of models is explored.

We discuss inflation models within supersymmetry and supergravity frameworks with a landscape of chiral superfields and one $U(1)$ shift symmetry which is broken by non-perturbative symmetry breaking terms in the superpotential. We label the pseudo scalar component of the chiral fields axions and their real parts saxions. Thus in the models only one combination of axions will be a pseudo-Nambu-Goldstone-boson which will act as the inflaton. The proposed models constitute consistent inflation for the following reasons: The inflation potential arises dynamically with stabilized saxions, the axion decay constant can lie in the sub-Planckian region, and consistency with the Planck data is achieved. The axion landscape consisting of $m$ axion pairs is assumed with the axions in each pair having opposite charges. A fast roll--slow roll splitting mechanism for the axion potential is proposed which is realized with a special choice of the axion basis. In this basis the $2m$ coupled equations split into $2m-1$ equations which enter in the fast roll and there is one unique linear combination of the $2m$ fields which controls the slow roll and thus the power spectrum of curvature and tensor perturbations. It is shown that a significant part of the parameter space exists where inflation is successful, i.e., $N_{\rm pivot} = {[50, 60]}$, the spectral index $n_s$ of curvature perturbations, and the ratio $r$ of the power spectrum of tensor perturbations and curvature perturbations, lie in the experimentally allowed regions given by the Planck experiment. Further, it is shown that the model allows for a significant region of the parameter space where the effective axion decay constant can lie in the sub-Planckian domain.

The luminosity of fading type Ia supernovae is governed by radioactive decays of 56Ni and 56Co. The decay rates are proportional to the Fermi coupling constant G_F and, therefore, are determined by the vacuum expectation value v of the Brout-Englert-Higgs field. We use publicly available sets of light curves of type Ia supernova at various redshifts to constrain possible spacetime variations of the 56Ni decay rate. The resulting constraint is not very tight; however, it is the only direct bound on the variation of the decay rate for redshifts up to z~1. We discuss potential applications of the result to searches for non-constancy of G_F and v.

It is known that the high energy $(> 100\,\text{MeV})$ emission of gamma-ray bursts is delayed with respect to the low energy emission. However, the dependence of light curves on energy has not been studied for the high energy bands. In this paper we consider the bursts observed by Fermi LAT from 2008 August 4 to 2011 August 1, for which at least $10$ photons were observed with the energy greater than $1\,\text{GeV}$. These include $4$ bursts: GRB 080916C, GRB 090510, GRB 090902B, and GRB 090926A. We use the Kolmogorov-Smirnov test to compare the light curves in the two bands, $100\,\text{MeV} < E < 1\,\text{GeV}$ and $E > 1\,\text{GeV}$. For GRB 080916C and GRB 090510 the light curves in the two bands are statistically compatible. However, for GRB 090926A, the higher-energy light curve is stretched compared to the lower-energy one with a statistical significance of $3.3 \sigma$ and, for GRB 090902B, on the contrary, the lower-energy curve is stretched with $2.3 \sigma$ significance. We argue that the observed diversity of stretching factors may be explained in a simple geometrical model. The model assumes that the jet opening angle depends on the emission energy in a way that the most energetic photons are radiated near the axis of the jet. All the bursts are considered equivalent in their rest frames and the observed light curves differ only due to different redshifts and view directions. The model conforms to the total burst energy constraint and matches the Fermi LAT observations of the fraction of GRBs visible in $100\,\text{MeV} < E < 1\,\text{GeV}$ band, which may be observed at higher energies. The model predicts the distribution of observable stretching factors, which may be tested in the future data. Finally, we propose a method to estimate observer's off-axis angle based on the stretching factor and the fraction of the high-energy photons.