In Einstein's general relativity, gravity is mediated by a massless metric field. The extension of general relativity to consistently include a mass for the graviton has profound implications for gravitation and cosmology. Salient features of various massive gravity theories can be captured by Galileon models, the simplest of which is the cubic Galileon. The presence of the Galileon field leads to additional gravitational radiation in binary pulsars where the Vainshtein mechanism is less suppressed than its fifth-force counterpart, which deserves a detailed confrontation with observations. We prudently choose fourteen well-timed binary pulsars, and from their intrinsic orbital decay rates we put a new bound on the graviton mass, $m_g \lesssim 2 \times 10^{-28}\,{\rm eV}/c^2$ at the 95% confidence level, assuming a flat prior on $\ln m_g$. It is equivalent to a bound on the graviton Compton wavelength $\lambda_g \gtrsim 7 \times 10^{21}\,{\rm m}$. Furthermore, we extensively simulate times of arrival for pulsars in orbit around stellar-mass black holes and the supermassive black hole at the Galactic center, and investigate their prospects in probing the cubic Galileon theory in the near future.

In this paper, we propose a hexagonal description for the flavor composition of ultrahigh-energy (UHE) neutrinos and antineutrinos, which will hopefully be determined at the future large neutrino telescopes. With such a geometrical description, we are able to clearly separate the individual flavor composition of neutrinos from that of antineutrinos in one single regular hexagon, which can be regarded as a natural generalization of the widely-used ternary plot. For illustration, we consider the $pp$ or $p\gamma$ collisions as the dominant production mechanism for UHE neutrinos and antineutrinos in the cosmic accelerator, and investigate how neutrino oscillations in the standard picture and in the presence of Lindblad decoherence could change the flavor composition of neutrinos and antineutrinos at neutrino telescopes.

We present timing solutions and analyses of 11 pulsars discovered by the Five-hundred-meter Aperture Spherical radio Telescope (FAST). These pulsars were discovered using an ultra-wide bandwidth receiver in drift-scan observations made during the commissioning phase of FAST, and were then confirmed and timed using the 64-m Parkes Radio Telescope. Each pulsar has been observed over a span of at least one year. Highlighted discoveries include PSR J0344-0901, which displays mode-changing behaviour and may belong to the class of so-called `swooshing' pulsars (alongside PSRs B0919+06 and B1859+07); PSR J0803-0942, whose emission is almost completely linearly polarised; and PSRs J1900-0134 and J1945+1211, whose well defined polarisation angle curves place stringent constraints on their emission geometry. We further discuss the detectability of these pulsars by earlier surveys, and highlight lessons learned from our work in carrying out confirmation and monitoring observations of pulsars discovered by a highly sensitive telescope, many of which may be applicable to next-generation pulsar surveys. This paper marks one of the first major releases of FAST-discovered pulsars, and paves the way for future discoveries anticipated from the Commensal Radio Astronomy FAST Survey (CRAFTS).

We investigate the stability and $e^+e^-$ pair creation of supercritically charged superheavy nuclei, $ud$QM nuggets, strangelets, and strangeon nuggets based on Thomas-Fermi approximation. The model parameters are fixed by reproducing their masses and charge properties reported in earlier publications. It is found that $ud$QM nuggets, strangelets, and strangeon nuggets may be more stable than ${}^{56}$Fe at $A\gtrsim 315$, $5\times10^4$, and $1.2\times10^8$, respectively. For those stable against neutron emission, the most massive superheavy element has a baryon number $\sim$965, while $ud$QM nuggets, strangelets, and strangeon nuggets need to have baryon numbers larger than $39$, 433, and $2.7\times10^5$. The $e^+e^-$ pair creation will inevitably start for superheavy nuclei with charge numbers $Z\geq177$, $ud$QM nuggets with $Z\geq163$, strangelets with $Z\geq 192$, and strangeon nuggets with $Z\geq 212$. A universal relation $Q/R_e = \left(m_e - \bar{\mu}_e\right)/\alpha$ is obtained at a given electron chemical potential $\bar{\mu}_e$, where $Q$ is the total charge and $R_e$ the radius of electron cloud. This predicts the maximum charge number by taking $\bar{\mu}_e=-m_e$. For supercritically charged objects with $\bar{\mu}_e<-m_e$, the decay rate for $e^+e^-$ pair production is estimated based on the JWKB approximation. It is found that most positrons are emitted at $t\lesssim 10^{-15}$ s, while a long lasting positron emission is observed for large objects with $R\gtrsim 1000$ fm. The emission and annihilation of positrons from supercritically charged objects may be partially responsible for the short $\gamma$-ray burst during the merger of binary compact stars, the 511 keV continuum emission, as well as the narrow faint emission lines in X-ray spectra from galaxies and galaxy clusters.

We develop a method to estimate the dust attenuation curve of galaxies from full spectral fitting of their optical spectra. Motivated from previous studies, we separate the small-scale features from the large-scale spectral shape, by performing a moving average method to both the observed spectrum and the simple stellar population model spectra. The intrinsic dust-free model spectrum is then derived by fitting the observed ratio of the small-scale to large-scale (S/L) components with the S/L ratios of the SSP models. The selective dust attenuation curve is then determined by comparing the observed spectrum with the dust-free model spectrum. One important advantage of this method is that the estimated dust attenuation curve is independent of the shape of theoretical dust attenuation curves. We have done a series of tests on a set of mock spectra covering wide ranges of stellar age and metallicity. We show that our method is able to recover the input dust attenuation curve accurately, although the accuracy depends slightly on signal-to-noise ratio of the spectra. We have applied our method to a number of edge-on galaxies with obvious dust lanes from the ongoing MaNGA survey, deriving their dust attenuation curves and $E(B-V)$ maps, as well as dust-free images in $g$, $r$, and $i$ bands. These galaxies show obvious dust lane features in their original images, which largely disappear after we have corrected the effect of dust attenuation. The vertical brightness profiles of these galaxies become axis-symmetric and can well be fitted by a simple model proposed for the disk vertical structure. Comparing the estimated dust attenuation curve with the three commonly-adopted model curves, we find that the Calzetti curve provides the best description of the estimated curves for the inner region of galaxies, while the Milky Way and SMC curves work better for the outer region.

Before massive stars heavier than $(8 \cdots 10)$ solar masses evolve to the phase of a gravitational core collapse, they will emit a huge number of MeV-energy neutrinos that are mainly produced in the thermal processes and nuclear weak interactions. The detection of such pre-supernova (pre-SN) neutrinos could provide an important and independent early warning for the optical observations of core-collapse SNe. In this paper, we investigate the prospects of future large liquid-scintillator detectors for the observation of pre-SN neutrinos in both $\bar{\nu}^{}_e + p \to e^+ + n$ and $\nu (\bar{\nu}) + e^- \to \nu (\bar{\nu}) + e^-$ reaction channels, where $\nu$ ($\bar{\nu}$) denotes neutrinos (antineutrinos) of all three flavors. We propose a quantitative assessment of the capability in terms of three working criteria, namely, how far the SN distance can be covered, how long the early warning before the core collapse can be sent out, and how well the direction pointing to the SN can be determined. The dependence of the final results on the different models of progenitor stars, neutrino flavor conversions and the relevant backgrounds is also discussed.

Mass-varying massive gravity allows the graviton mass to vary according to different environments. We investigate neutron star and white dwarf solutions in this theory and find that the graviton mass can become very large near the compact stars and settle down quickly to small cosmological values away the stars, similar to that of black holes in the theory. It is found that there exists a tower of compact star solutions where the graviton mass decreases radially to zero non-trivially. We compute the massive graviton effects on the mass-radius relations of the compact stars, and also compare the relative strengths between neutron stars and white dwarfs in constraining the parameter space of mass-varying massive gravity.

We measure the star formation histories (SFH) of a sample of low-mass galaxies with $M_\ast<10^9M_\odot$ from the SDSS-IV MaNGA survey. The large number of IFU spectra for each galaxy are either combined to reach a high signal to noise ratio or used to investigate spatial variations. We use Bayesian inferences based on full spectrum fitting. Our analysis based on Bayesian evidence ratio indicates a strong preference for a model that allows the presence of an old stellar population, and that an improper model for the SFH can significantly underestimate the old population in these galaxies. The addition of NIR photometry to the constraining data can further distinguish between different SFH model families and significantly tighten the constraints on the mass fraction in the old population. On average more than half of the stellar mass in present-day low-mass galaxies formed 8 Gyrs ago, while about 30\% within the past 4 Gyrs. Satellite galaxies on average have formed their stellar mass earlier than central galaxies, and stars in the outer regions of galaxy are younger than those in the central part. Our results suggest that most of the low-mass galaxies have an early episode of active star formation that produces a large fraction of their present stellar mass.

Within the so-called cosmological collider program, imprints of new particles on primordial non-Gaussianities have been studied intensively. In particular, their non-analytic features in the soft limit provide a smoking gun for new particles at the inflation scale. While this approach is very powerful to probe particles of the mass near the Hubble scale, the signal is exponentially suppressed for heavy particles. In this paper, to enlarge the scope of the cosmological collider, we explore a new approach to probing spins of heavy particles from signs of Wilson coefficients of the inflaton effective action and the corresponding primordial non-Gaussianities. As a first step, we focus on the regime where the de Sitter conformal symmetry is weakly broken. It is well known that the leading order effective operator $(\partial_\mu\phi\partial^\mu\phi)^2$ is universally positive as a consequence of unitarity. In contrast, we find that the sign of the six derivative operator $(\nabla_\mu\partial_\nu\phi)^2(\partial_\rho\phi)^2$ is positive for intermediate heavy scalars, whereas it is negative for intermediate heavy spinning states. Therefore, under the assumption of tree-level UV completion, the sign can be used to probe spins of heavy particles generating the effective interaction. We also study phenomenology of primordial non-Gaussianities thereof.

We discuss two possible scenarios, namely the curvaton mechanism and the dark matter density modulation, where non-Gaussianity signals of superheavy dark matter produced by gravity can be enhanced and observed. In both scenarios, superheavy dark matter couples to an additional light field as a mediator. In the case of derivative coupling, the resulting non-Gaussianities induced by the light field can be large, which can provide inflationary evidences for these superheavy dark matter scenarios.

Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) is acquiring integral-field spectroscopy for the largest sample of galaxies to date. By 2020, the MaNGA Survey --- one of three core programs in the fourth-generation Sloan Digital Sky Survey (SDSS-IV) --- will have observed a statistically representative sample of 10$^4$ galaxies in the local Universe ($z\lesssim0.15$). In addition to a robust data-reduction pipeline (DRP), MaNGA has developed a data-analysis pipeline (DAP) that provides higher-level data products. To accompany the first public release of its code base and data products, we provide an overview of the MaNGA DAP, including its software design, workflow, measurement procedures and algorithms, performance, and output data model. In conjunction with our companion paper Belfiore et al., we also assess the DAP output provided for 4718 observations of 4648 unique galaxies in the recent SDSS Data Release 15 (DR15). These analysis products focus on measurements that are close to the data and require minimal model-based assumptions. Namely, we provide stellar kinematics (velocity and velocity dispersion), emission-line properties (kinematics, fluxes, and equivalent widths), and spectral indices (e.g., D4000 and the Lick indices). We find that the DAP provides robust measurements and errors for the vast majority ($>$99%) of analyzed spectra. We summarize assessments of the precision and accuracy of our measurements as a function of signal-to-noise, and provide specific guidance to users regarding the limitations of the data. The MaNGA DAP software is publicly available and we encourage community involvement in its development.

We study the spectrum of semiclassical rotating strings on de Sitter space and its consistency. Even though a naive extrapolation of the linear Regge trajectory on flat space implies a violation of the Higuchi bound (a unitarity bound on the mass of higher-spin particles in de Sitter space), the curved space effects turn out to modify the trajectory to respect the bound. Interestingly, there exists a maximum spin for each Regge trajectory as a consequence of accelerated expansion, which is helpful to make the spectrum consistent with the Higuchi bound, but at the same time it could be an obstruction to stringy UV completion based on an infinite higher-spin tower. By pushing further this observation, we demonstrate that the vacuum energy $V$ inflating the universe has to be bounded by the string scale $M_s$ as $V\lesssim M_s^4$, if UV completion is achieved by the leading Regge trajectory with higher spin states up to the 4D Planck scale. Its application to inflation at the early universe implies an upper bound on the tensor-to-scalar ratio, $r\lesssim 0.01\times(M_s/10^{16} \text{GeV})^{4}$, which is within the scope of the near future CMB experiments. We also discuss another possibility that UV completion is achieved by multiple Regge trajectories.

We study the Schwinger effect during inflation and its imprints on the primordial power spectrum and bispectrum. The produced charged particles by Schwinger effect during inflation can leave a unique angular dependence on the primordial spectra.

We apply the recently developed positivity bounds for particles with spin, applied away from the forward limit, to the low energy effective theories of massive spin-1 and spin-2 theories. For spin-1 theories, we consider the generic Proca EFT which arises at low energies from a heavy Higgs mechanism, and the special case of a charged Galileon for which the EFT is reorganized by the Galileon symmetry. For spin-2, we consider generic $\Lambda_5$ massive gravity theories and the special `ghost-free' $\Lambda_3$ theories. Remarkably we find that at the level of 2-2 scattering, the positivity bounds applied to $\Lambda_5$ massive gravity theories impose the special tunings which generate the $\Lambda_3$ structure. For $\Lambda_3$ massive gravity theories, the island of positivity derived in the forward limit appears relatively stable against further bounds.

Gravitational wave echoes, if they exist, could encode important information of new physics from the strong gravity regime. Current echo searches usually assume constant interval echoes (CIEs) a priori, although unequal interval echoes (UIEs) are also possible. Despite of its simplicity, the using of CIE templates need to be properly justified, especially given the high sensitivity of future gravitational wave detectors. In this paper, we assess the necessity of UIE templates in echo searches. By reconstructing injected UIE signals with both CIE and UIE templates, we show that the CIE template may significantly misinterpret the echo signals if the variation of the interval is greater than the statistical errors of the interval, which is further confirmed by a Bayesian analysis on model stelection. We also forecast the constraints on the echo intervals given by future GW detectors such as Advanced LIGO and Einstein Telescope.

We study the gravitational production of super-Hubble-mass dark matter in the very early universe. We first review the simplest scenario where dark matter is produced mainly during slow roll inflation. Then we move on to consider the cases where dark matter is produced during the transition period between inflation and the subsequent cosmological evolution. The limits of smooth and sudden transitions are studied, respectively. The relic abundances and the cosmological collider signals are calculated.

The model-independent reconstruction of the energy spectra of $\overline{\nu}^{}_e$, $\nu^{}_e$ and $\nu^{}_x$ (i.e., $\nu^{}_\mu$, $\nu^{}_\tau$ and their antiparticles) from the future observation of a galactic core-collapse supernova (SN) is of crucial importance to understand the microscopic physics of SN explosions. To this end, we propose a practically useful method to combine the multi-channel detection of SN neutrinos in a large liquid-scintillator detector (e.g., JUNO), namely, the inverse beta decay $\overline{\nu}^{}_e + p \to e^+ + n$, the elastic neutrino-proton scattering $\nu + p \to \nu + p$ and the elastic neutrino-electron scattering $\nu + e^- \to \nu + e^-$, and reconstruct the energy spectra of $\overline{\nu}^{}_e$, $\nu^{}_e$ and $\nu^{}_x$ by making the best use of the observational data in those three channels. In addition, the neutrino energy spectra from the numerical simulations of the delayed neutrino-driven SN explosions are implemented to demonstrate the robustness of our method. Taking the ordinary matter effects into account, we also show how to extract the initial neutrino energy spectra in the presence of neutrino flavor conversions.

We analyze the stellar initial mass functions (IMF) of a large sample of early type galaxies (ETGs) provided by MaNGA. The large number of IFU spectra of individual galaxies provide high signal-to-noise composite spectra that are essential for constraining IMF and to investigate possible radial gradients of the IMF within individual galaxies. The large sample of ETGs also make it possible to study how the IMF shape depends on various properties of galaxies. We adopt a novel approach to study IMF variations in ETGs, use Bayesian inferences based on full spectrum fitting. The Bayesian method provides a statistically rigorous way to explore potential degeneracy in spectrum fitting, and to distinguish different IMF models with Bayesian evidence. We find that the IMF slope depends systematically on galaxy velocity dispersion, in that galaxies of higher velocity dispersion prefer a more bottom-heavy IMF, but the dependence is almost entirely due to the change of metallicity, $Z$, with velocity dispersion. The IMF shape also depends on stellar age, $A$, but the dependence is completely degenerate with that on metallicity through a combination $AZ^{-1.42}$. Using independent age and metallicity estimates we find that the IMF variation is produced by metallicity instead of age. The IMF near the centers of massive ETGs appears more bottom-heavy than that in the outer parts, while a weak opposite trend is seen for low-mass ETGs. Uncertainties produced by star formation history, dust extinction, $\alpha$-element abundance enhancement and noise in the spectra are tested.

The multi-ton-scale liquid xenon detectors, with an excellent energy resolution of a few keV, will be constructed to probe the dark-matter particles. In this paper, we show that precision measurements of the low-energy solar neutrinos via the elastic neutrino-electron scattering in this kind of detectors are able to improve the present limits on neutrino lifetimes and neutrino magnetic moments by about one order of magnitude. We carefully study the impact of the unknown neutrino mass spectrum on the ultimate limits in the case of non-radiative visible neutrino decays. In the case of invisible neutrino decays, the lower bounds $\tau^{}_1/m^{}_1 \gtrsim 3 \times 10^{-2}~{\rm s}/{\rm eV}$ and $\tau^{}_2/m^{}_2 \gtrsim 8\times 10^{-3}~{\rm s}/{\rm eV}$ at the $2\sigma$ level can be obtained for a total exposure of $70~{\rm ton}\cdot {\rm year}$. Furthermore, a restrictive constraint on the effective magnetic moment of neutrinos $\mu^{}_{\rm eff} \lesssim 2.6\times 10^{-12}~\mu^{}_{\rm B}$, with $\mu^{}_{\rm B}$ being the Bohr magneton, can be achieved. This is among the best results that will be available in the laboratory experiments in the near future.

We use the spinor helicity formalism to study KLT-like relations for the inflationary graviton four-point correlation function. New features are observed in this correlation function compared to the graviton scattering amplitude in flat spacetime. After obtaining the general momentum dependence, collinear, squeezed and collapsed limits are considered to further study the features of the correlation function, and the relation to the corresponding flat space scattering amplitude.

We develop the systematics for applying operators on Minkowski correlation functions to get the inflationary correlation functions. Simple structures and recursion relations are known for Minkowski correlation functions. Using the operator technique, various novel recursion relations for inflationary correlation functions are obtained.

Halo merger trees are constructed from ELUCID, a constrained $N$-body simulation in the Sloan Digital Sky Survey (SDSS) volume. These merger trees are used to populate dark matter halos with galaxies according to an empirical model of galaxy formation. Mock catalogs in the SDSS sky coverage are constructed, which can be used to study the spatial distribution of galaxies in the low-$z$ Universe. These mock catalogs are used to quantify the cosmic variance in the galaxy stellar mass function (GSMF) measured from the SDSS survey. The GSMF estimated from the SDSS magnitude-limited sample can be affected significantly by the presence of the under-dense region at $z<0.03$, so that the low-mass end of the function can be underestimated significantly. Several existing methods designed to deal with the effects of the cosmic variance in the estimate of GSMF are tested, and none is found to be able to fully account for the cosmic variance. We propose a method based on the conditional stellar mass functions in dark matter halos, which can provide an unbiased estimate of the global GSMF. The application of the method to the SDSS data shows that the GSMF has a significant upturn at $M_*< 10^{9.5} h^{-1}{\rm M}_\odot$, which has been missed in many earlier measurements of the local GSMF.

It has been recently suggested that oscillons produced in the early universe from certain asymmetric potentials continue to emit gravitational waves for a number of $e$-folds of expansion after their formation, leading to potentially detectable gravitational wave signals. We revisit this claim by conducting a convergence study using graphics processing unit (GPU)-accelerated lattice simulations and show that numerical errors accumulated with time are significant in low-resolution scenarios, or in scenarios where the run-time causes the resolution to drop below the relevant scales in the problem. Our study determines that the dominant, growing high frequency peak of the gravitational wave signals in the fiducial "hill-top model" in [arXiv:1607.01314] is a numerical artifact. This finding prompts the need for a more careful analysis of the numerical validity of other similar results related to gravitational waves from oscillon dynamics.

Data are gathered from the Parkes pulsar data archive of twelve young radio pulsars, with the intervals of data for each pulsar ranged between 2.8 years and 6.8 years. 31 glitches are identified by using phase connection from "pulsar timing" technology, ranging from $1.7\times10^{-9}$ to $8.5\times10^{-6}$ at the change in relative glitch sizes $\Delta\nu/\nu$, where $\nu=1/P$ is the pulse frequency. 8 post-glitch behaviours of 13 published glitches are updated with 4 exponential recoveries observed. Detecting 18 new glitches has great significance, but only one exponential relaxation was discovered in these new glitches. The bimodal distribution of $\Delta\nu/\nu$, as before, is validated with peaks at $\sim 10^{-6}$ and $\sim 10^{-9}$. Moreover, all exponential decays were observed in large glitches, and have very low $Q$. It takes short time toward the extrapolation of the pre-glitch pulse frequency with the timescale of 40 or 80 days. Besides, most glitches exhibit a linear decrease in slow-down rate $|\dot{\nu}|$ and a permanent change in $\ddot{\nu}$ after glitch. In special, PSR J1341-6220 was detected 4 new glitches to compose a total of 27 glitches, rising in third place of the most actively glitching pulsars known with a glitch rate of $1.17$ glitches per year. Unusual post-glitch behaviours demonstrate long-lasting increase of $\nu$ for hundreds of days in PSR J1112--6103 and PSR J1614--5048.

Gravitational wave (GW) ringdown waveforms may contain "echoes" that encode new physics in the strong gravity regime. It is commonly assumed that the new physics gives rise to the GW echoes whose intervals are constant. We point out that this assumption is not always applicable. In particular, if the post-merger object is initially a wormhole, which slowly pinches off and eventually collapses into a black hole, the late-time ringdown waveform exhibit a series of echoes whose intervals are increasing with time. We also assess how this affects the ability of Advanced LIGO/Virgo to detect these new signals.

The decay of massive particles during inflation generates characteristic signals in the squeezed limit of the primordial bispectrum. These signals are in particular distinctive in the regime of the quasi-single field inflation, where particles are oscillating with masses comparable to the Hubble scale. We apply the investigation to a class of scalar particles that experience a so-called waterfall phase transition in the isocurvature direction driven by the symmetry-breaking mechanism based on the hybrid inflation scenario. With a time-varying mass, a novel shape of oscillatory bispectrum is presented as the signature of a waterfall phase transition during inflation.

We show how cosmological correlation functions of massless fields can be rewritten in terms of Minkowski correlation functions, by extracting symmetry-breaking operators from the cosmological correlators. This technique simplifies some cosmological calculations. Also, known properties of Minkowski correlation functions can be translated to non-trivial properties of cosmological correlations. To illustrate this idea, inflation to Minkowski and matter bounce to Minkowski relations are presented for the interactions of general single field inflation. And a Minkowski recursion relation is translated to a novel relation for inflation.

In this paper, we show how to carry out a relatively more realistic and complete reconstruction of supernova neutrino spectra in the future large liquid-scintillator detectors, by implementing the method of singular value decomposition with a proper regularization. For a core-collapse supernova at a distance of $10~{\rm kpc}$ in the Milky Way, its $\overline{\nu}^{}_e$ spectrum can be precisely determined from the inverse beta-decay process $\overline{\nu}^{}_e + p \to e^+ + n$, for which a $20~{\rm kiloton}$ liquid-scintillator detector with the resolution similar to the Jiangmen Underground Neutrino Observatory (JUNO) may register more than 5000 events. We have to rely predominantly on the elastic neutrino-electron scattering $\nu + e^- \to \nu + e^-$ and the elastic neutrino-proton scattering $\nu + p \to \nu + p$ for the spectra of $\nu^{}_e$ and $\nu^{}_x$, where $\nu$ denotes collectively neutrinos and antineutrinos of all three flavors and $\nu^{}_x$ for $\nu^{}_\mu$ and $\nu^{}_\tau$ as well as their antiparticles. To demonstrate the validity of our approach, we also attempt to reconstruct the neutrino spectra by using the time-integrated neutrino data from the latest numerical simulations of delayed neutrino-driven supernova explosions.

We investigate possible interactions between neutrinos and massive scalar bosons via $g^{}_{\phi} \overline{\nu} \nu \phi$ (or massive vector bosons via $g^{}_V \overline{\nu} \gamma^\mu \nu V^{}_\mu$) and explore the allowed parameter space of the coupling constant $g^{}_{\phi}$ (or $g^{}_V$) and the scalar (or vector) boson mass $m^{}_\phi$ (or $m^{}_V$) by requiring that these secret neutrino interactions (SNIs) should not spoil the success of Big Bang nucleosynthesis (BBN). Incorporating the SNIs into the evolution of the early Universe in the BBN era, we numerically solve the Boltzmann equations and compare the predictions for the abundances of light elements with observations. It turns out that the constraint on $g^{}_{\phi}$ and $m^{}_\phi$ in the scalar-boson case is rather weak, due to a small number of degrees of freedom. However, in the vector-boson case, the most stringent bound on the coupling $g^{}_V \lesssim 6\times 10^{-10}$ at $95~\%$ confidence level is obtained for $m^{}_V \simeq 1~{\rm MeV}$, while the bound becomes much weaker $g^{}_V \lesssim 8\times 10^{-6}$ for smaller masses $m^{}_V \lesssim 10^{-4}~{\rm MeV}$. Moreover, we discuss in some detail how the SNIs affect the cosmological evolution and the abundances of the lightest elements.

We study quasi-single field inflation with a warm inflation background. The thermal effects at small scales can sufficiently enhance the magnitude of the primordial standard clock signal. This scenario offers us the possibility of probing the UV physics of the very early universe without the exponentially small Boltzmann factor when the mass of the isocurvaton is much heavier than Hubble. The thermal effects at small scales can be studied using the flat space thermal field theory, connected to an effective description using non-Bunch-Davies vacuum at large scales, with large clock signal.

Further studies on the corner singularity of GS reconstruction are are compiled in this paper. It's focused on solution of the Data Completion (DC) problem with the Extended Hilbert Transform (EHT) over plane rectangular region. Optimal selections of the regular parameters in{\it Tikhonov} solution of corresponding DC problem are developed in this study. The 6-parameter regular solutions and the {\it Jacobian} matrix and one {\it Hessian} tensor to the regular parameters are derived in this work. A concise formula for EHT in the near field of corners, which shows property for EHT near the corner, are also provided. It serves as the additional constraints for our parameter optimization problem (OP). Third, a nonlinear convex function defined by the regular solution and the corner constraints is introduced and is taken as the object function for the OP of the 6 regular parameters on half-space $(\mathbf{p}>0)$. Given an initial guess of $\mathbf{p}_0$, the optimal parameters are solved from the OP through a well known constrained nonlinear optimization method. Last, the benchmark tests to the proposed solution approach are carried out, and detailed results from totally 9 different bench-cases are tabulated. In contrast to solutions with given regular parameters, our bench results demonstrate that an objective way for selection of the optimal $\mathbf{p}$ is successfully laid out here. Robustness and efficiency of the suggested new approach are also highlighted in this study.

We study the evolution of cosmological perturbations in a non-singular bouncing cosmology with a bounce phase which has superimposed oscillations of the scale factor. We identify length scales for which the final spectrum of fluctuations obtains imprints of the non-trivial bounce dynamics. These imprints in the spectrum are manifested in the form of damped oscillation features at scales smaller than a characteristic value and an increased reddening of the spectrum at all the scales as the number of small bounces increases.

We study the power spectrum of quasi-single field inflation where strong coupling is considered. The contribution from the massive propagator can be divided into local and non-local contributions. The local one is the leading contribution and is power-law suppressed as a function of mass, while the non-local contribution is exponentially suppressed in the large mass limit. For the local contribution, it is possible to use the effective field theory approach to study the power spectrum in the strongly coupled region of the parameter space. For the non-local contribution, we developed a partial effective field theory method to simplify the calculation: When there are multiple massive propagators, one can fully compute it after integrating out all but one massive propagator by effective field theory. The result retains the "standard clock" signal, which is interesting for probing the expansion history of the primordial universe and the physics of a "cosmological collider". The error involved compared to the full calculation is power law suppressed by the effective mass of the heavy field.

The mass of the graviton, if nonzero, is usually considered to be very small, e.g. of the Hubble scale, from several observational constraints. In this paper, we propose a gravity model where the graviton mass is very small in the usual weak gravity environments, below all the current graviton mass bounds, but becomes much larger in the strong gravity regime such as a black hole's vicinity. For black holes in this model, significant deviations from general relativity emerge very close to the black hole horizon and alter the black hole quasi-normal modes, which can be extracted from the ringdown waveform of black hole binary mergers. Also, the enhancement of the graviton mass near the horizon can result in echoes in the late time ringdown, which can be verified in the upcoming gravitational wave observations of higher sensitivity.

We study the effective field theory (EFT) description of the virtual particle effects in quasi-single field inflation, which unifies the previous results on large mass and large mixing cases. By using a horizon crossing approximation and matching with known limits, approximate expressions for the power spectrum and the spectral index are obtained. The error of the approximate solution is within 10% in dominate parts of the parameter space, which corresponds to less-than-0.1% error in the $n_s$-$r$ diagram. The quasi-single field corrections on the $n_s$-$r$ diagram are plotted for a few inflation models. Especially, the quasi-single field correction drives $m^2\phi^2$ inflation to the best fit region on the $n_s$-$r$ diagram, with an amount of equilateral non-Gaussianity which can be tested in future experiments.

Recently, aLIGO has announced the first direct detections of gravitational waves, a direct manifestation of the propagating degrees of freedom of gravity. The detected signals GW150914 and GW151226 have been used to examine the basic properties of these gravitational degrees of freedom, particularly setting an upper bound on their mass. It is timely to review what the mass of these gravitational degrees of freedom means from the theoretical point of view, particularly taking into account the recent developments in constructing consistent massive gravity theories. Apart from the GW150914 mass bound, a few other observational bounds have been established from the effects of the Yukawa potential, modified dispersion relation and fifth force that are all induced when the fundamental gravitational degrees of freedom are massive. We review these different mass bounds and examine how they stand in the wake of recent theoretical developments and how they compare to the bound from GW150914.

The EFT coefficients in any gapped, scalar, Lorentz invariant field theory must satisfy positivity requirements if there is to exist a local, analytic Wilsonian UV completion. We apply these bounds to the tree level scattering amplitudes for a massive Galileon. The addition of a mass term, which does not spoil the non-renormalization theorem of the Galileon and preserves the Galileon symmetry at loop level, is necessary to satisfy the lowest order positivity bound. We further show that a careful choice of successively higher derivative corrections are necessary to satisfy the higher order positivity bounds. There is then no obstruction to a local UV completion from considerations of tree level 2-to-2 scattering alone. To demonstrate this we give an explicit example of such a UV completion.

Assuming the existence of a local, analytic, unitary UV completion in a Poincar\'{e} invariant scalar field theory with a mass gap, we derive an infinite number of positivity requirements using the known properties of the amplitude at and away from the forward scattering limit. These take the form of bounds on combinations of the pole subtracted scattering amplitude and its derivatives. In turn, these positivity requirements act as constraints on the operator coefficients in the low energy effective theory. For certain theories these constraints can be used to place an upper bound on the mass of the next lightest state that must lie beyond the low energy effective theory if such a UV completion is to ever exist.

We present a comprehensive review of keV-scale sterile neutrino Dark Matter, collecting views and insights from all disciplines involved - cosmology, astrophysics, nuclear, and particle physics - in each case viewed from both theoretical and experimental/observational perspectives. After reviewing the role of active neutrinos in particle physics, astrophysics, and cosmology, we focus on sterile neutrinos in the context of the Dark Matter puzzle. Here, we first review the physics motivation for sterile neutrino Dark Matter, based on challenges and tensions in purely cold Dark Matter scenarios. We then round out the discussion by critically summarizing all known constraints on sterile neutrino Dark Matter arising from astrophysical observations, laboratory experiments, and theoretical considerations. In this context, we provide a balanced discourse on the possibly positive signal from X-ray observations. Another focus of the paper concerns the construction of particle physics models, aiming to explain how sterile neutrinos of keV-scale masses could arise in concrete settings beyond the Standard Model of elementary particle physics. The paper ends with an extensive review of current and future astrophysical and laboratory searches, highlighting new ideas and their experimental challenges, as well as future perspectives for the discovery of sterile neutrinos.

A characteristic signature of String Gas Cosmology is primordial power spectra for scalar and tensor modes which are almost scale-invariant but with a red tilt for scalar modes but a blue tilt for tensor modes. This feature, however, can also be realized in the so-called G-inflation model, in which Horndeski operators are introduced which leads to a blue tensor tilt by softly breaking the Null Energy Condition. In this article we search for potential observational differences between these two cosmologies by performing detailed perturbation analyses based on the Effective Field Theory approach. Our results show that, although both two models produce blue tilted tensor perturbations, they behave differently in three aspects. Firstly, String Gas Cosmology predicts a specific consistency relation between the index of the scalar modes ns and that of tensor ones nt, which is hard to be reproduced by G-inflation. Secondly, String Gas Cosmology typically predicts non-Gaussianities which are highly suppressed on observable scales, while G-inflation gives rise to observationally large non-Gaussianities because the kinetic terms in the action become important during inflation. However, after finely tuning the model parameters of G-inflation it is possible to obtain a blue tensor spectrum and negligible non-Gaussianities with a degeneracy between the two models. This degeneracy can be broken by a third observable, namely the scale dependence of the nonlinearity parameter, which vanishes for G-inflation but has a blue tilt in the case of String Gas Cosmology. Therefore, we conclude that String Gas Cosmology is in principle observationally distinguishable from the single field inflationary cosmology, even allowing for modifications such as G-inflation.

We use holographic techniques to compute inflationary non-Gaussianities for general single-field inflation, including models with a non-trivial sound speed. In this holographic approach, the inflationary dynamics is captured by a relevant deformation of the dual conformal field theory (CFT) in the UV, while the inflationary correlators are computed by conformal perturbation theory. In this paper, we discuss the effects of higher derivative operators, such as $(\partial_\mu\phi\partial^\mu\phi)^{m}$, which are known to induce a non-trivial sound speed and source potentially large non-Gaussianities. We compute the full inflationary bispectra from the deformed CFT correlators. We also discuss the squeezed limit of the bispectra from the viewpoint of operator product expansions. As is generic in the holographic description of inflation, our power spectrum is blue tilted in the UV region. We extend our bispectrum computation to the IR region by resumming the conformal perturbations to all orders. We provide a self-consistent setup which reproduces a red tilted power spectrum, as well as all possible bispectrum shapes in the slow-roll regime.

If massive neutrinos are Dirac particles, the proposed PTOLEMY experiment will hopefully be able to discover cosmic neutrino background via $\nu^{}_e + {^3}{\rm H} \to {^3}{\rm He} + e^-$ with a capture rate of $\Gamma^{}_{\rm D} \approx 4~{\rm yr}^{-1}$. Recently, it has been pointed out that right-handed components of Dirac neutrinos could also be copiously produced in the early Universe and become an extra thermal or nonthermal ingredient of cosmic relic neutrinos, enhancing the capture rate to $\Gamma^{}_{\rm D} \approx 5.1~{\rm yr}^{-1}$ or $\Gamma^{}_{\rm D} \approx 6.1~{\rm yr}^{-1}$. In this work, we investigate the possibility to distinguish between thermal and nonthermal spectra of cosmic relic neutrinos by measuring the annual modulation of the capture rate. For neutrino masses of $0.1~{\rm eV}$, we have found the amplitude of annual modulation in the standard case is ${\cal M} \approx 0.05\%$, which will be increased to $0.1\%$ and $0.15\%$ in the presence of additional thermal and nonthermal right-handed neutrinos, respectively. The future detection of such a modulation will be helpful in understanding the Majorana or Dirac nature of massive neutrinos.

Lorentz-invariant massive gravity is usually associated with a strong coupling scale $\Lambda_3$. By including non-trivial effects from the Stueckelberg modes, we show that about these vacua, one can push the strong coupling scale to higher values and evade the linear vDVZ-discontinuity. For generic parameters of the theory and generic vacua for the Stueckelberg fields, the $\Lambda_2$-decoupling limit of the theory is well-behaved and free of any ghost or gradient-like instabilities. We also discuss the implications for nonlinear sigma models with Lorentzian target spaces.

The target space of a nonlinear sigma model is usually required to be positive definite to avoid ghosts. We introduce a unique class of nonlinear sigma models where the target space metric has a Lorentzian signature, thus the associated group being non-compact. We show that the would-be ghost associated with the negative direction is fully projected out by 2 second-class constraints, and there exist stable solutions in this class of models. This result also has important implications for Lorentz--invariant massive gravity: There exist stable nontrivial vacua in massive gravity that are free from any linear vDVZ--discontinuity and a $\Lambda_2$ decoupling limit can be defined on these vacua.

It remains to be determined experimentally if massive neutrinos are Majorana or Dirac particles. In this connection, it has been recently suggested that the detection of cosmic neutrino background of left-handed neutrinos $\nu^{}_{\rm L}$ and right-handed antineutrinos $\overline{\nu}^{}_{\rm R}$ in future experiments of neutrino capture on beta-decaying nuclei (e.g., $\nu^{}_e + {^3{\rm H}} \to {^3}{\rm He} + e^-$ for the PTOLEMY experiment) is likely to distinguish between Majorana and Dirac neutrinos, since the capture rate is twice larger in the former case. In this paper, we investigate the possible impact of right-handed neutrinos on the capture rate, assuming that massive neutrinos are Dirac particles and both right-handed neutrinos $\nu^{}_{\rm R}$ and left-handed antineutrinos $\overline{\nu}^{}_{\rm L}$ can be efficiently produced in the early Universe. It turns out that the capture rate can be enhanced at most by $28\%$ due to the presence of relic $\nu^{}_{\rm R}$ and $\overline{\nu}^{}_{\rm L}$ with a total number density of $95~{\rm cm}^{-3}$, which should be compared to the number density $336~{\rm cm}^{-3}$ of cosmic neutrino background. The enhancement has actually been limited by the latest cosmological and astrophysical bounds on the effective number of neutrino generations $N^{}_{\rm eff} = 3.14^{+0.44}_{-0.43}$ at the $95\%$ confidence level. For illustration, two possible scenarios have been proposed for thermal production of right-handed neutrinos in the early Universe.

In de Rham-Gabadadze-Tolley (dRGT) massive gravity and bigravity, a non-minimal matter coupling involving both metrics generically re-introduces the Boulware--Deser (BD) ghost. A non-minimal matter coupling via a simple, yet specific composite metric has been proposed, which eliminates the BD ghost below the strong coupling scale. Working explicitly in the metric formulation and for arbitrary spacetime dimensions, we show that this composite metric is the unique consistent non-minimal matter coupling below the strong coupling scale, which emerges out of two diagnostics, namely, absence of Ostrogradski ghosts in the decoupling limit and absence of the BD ghost from matter quantum loop corrections.

We compute the stochastic gravitational wave production from Affleck-Dine condensate fragmentation in the early universe, focusing on an effective potential with a logarithmic mass correction that typically arises in gravity mediated supersymmetry breaking scenarios. We find that a significant gravitational wave background can be generated when Q-balls are being formed out of the condensate fragmentation. This gravitational wave background has a distinct multi-peak power spectrum where the trough is closely linked to the supersymmetry breaking scale and whose frequencies are peaked around kHz for TeV supersymmetry breaking.

Motivated by the string landscape, inflation may happen on a high dimensional complicated potential. We propose a new way to construct some high dimensional random potentials, and study inflation on top of that, for up to 50-dimensions in field space. Especially, random bifurcations of classical inflationary trajectory are studied. It is shown that the bifurcation probability increases as a function of number of dimensions. Those random bifurcations are not consistent with observations, and dramatically limit the parameter space of inflation on a complicated landscape. For example, in 10 dimensions, only $10^{-3} \sim 10^{-6}$ of the parameter space volume leads to unique classical trajectories. The rest is ruled out by random bifurcations.

We investigate the coupling between the inflaton and massive vector fields. All renormalizable couplings with shift symmetry of the inflaton are considered. The massive vector can be decomposed into a scalar mode and a divergence-free vector mode. We show that the former naturally interacts with the inflaton and the latter decouples at tree level. The model in general predicts $f_{NL}^\mathrm{equil} = \mathcal{O}(1)$, while in some regions of the parameter space large non-Gaussianity can arise.

In a recent Letter we have shown that in shift-symmetric Horndeski theory the scalar field is forced to obtain a nontrivial configuration in black hole spacetimes, unless a linear coupling with the Gauss-Bonnet invariant is tuned away. As a result, black holes generically have hair in this theory. In this companion paper, we first review our argument and discuss it in more detail. We then present actual black hole solutions in the simplest case of a theory with the linear scalar-Gauss-Bonnet coupling. We generate exact solutions numerically for a wide range of values of the coupling and also construct analytic solutions perturbatively in the small coupling limit. Comparison of the two types of solutions indicates that non-linear effects that are not captured by the perturbative solution lead to a finite area, as opposed to a central, singularity. Remarkably, black holes have a minimum size, controlled by the length scale associated with the scalar-Gauss-Bonnet coupling. We also compute some phenomenological observables for the numerical solution for a wide range of values of the scalar-Gauss-Bonnet coupling. Deviations from the Schwarzschild geometry are generically very small.

Supersonic radiation diffusion approximation is a useful way to study the radiation transportation. Considering the bent Marshak wave theory in 2-dimensions, and an invariable source temperature, we get the supersonic radiation diffusion conditions which are about the Mach number $M>8(1+\sqrt{\ep})/3$, and the optical depth $\tau>1$. A large Mach number requires a high temperature, while a large optical depth requires a low temperature. Only when the source temperature is in a proper region these conditions can be satisfied. Assuming the material opacity and the specific internal energy depend on the temperature and the density as a form of power law, for a given density, these conditions correspond to a region about source temperature and the length of the sample. This supersonic diffusion region involves both lower and upper limit of source temperature, while that in 1-dimension only gives a lower limit. Taking $\rm SiO_2$ and the Au for example, we show the supersonic region numerically.

Q-balls are non-topological solitonic solutions to a wide class of field theories that possess global symmetries. Here we show that in these same theories there also exists a tower of novel composite Q-ball solutions where, within one composite Q-ball, positive and negative charges co-exist and swap at a frequency lower than the natural frequency of an individual Q-ball. These charge-swapping Q-balls are constructed by assembling Q-balls and anti-Q-balls tightly such that their nonlinear cores overlap. We explain why charge-swapping Q-balls can form and why they swap charges.

In this work we show that the general singlet extension of the MSSM can naturally provide a self-interacting singlino dark matter to solve the small cosmological scale anomalies (a large Sommerfeld enhancement factor can also be obtained). However, we find that the NMSSM (the singlet extension of the MSSM with Z3 symmetry) cannot achieve this due to the restricted parameter space. In our analysis we introduce the concept of symmetric and antisymmetric viscosity cross sections to deal with the non-relativistic Majorana-fermion dark matter scattering.

Light-like galileon solutions have been used to investigate the chronology problem in galileon-like theories, and in some cases may also be considered as solitons, evading a non-existence constraint from a zero-mode argument. Their stabilities have been analyzed via "local" approximation, which appears to suggest that all these light-like solutions are stable. We re-analyze the stability problem by solving the linear perturbation equation \emph{exactly}, and point out that the finite energy condition is essential for the light-like solitons to be stable. We also clarify potential ghost instabilities and why the zero-mode argument can not be naively generalized to include the light-like solitons.

In general relativity gravitational waves propagate at the speed of light, however in alternative theories of gravity that might not be the case. We study the effects of a modified speed of gravity, $c_T^2$, on the B-modes of the Cosmic Microwave Background (CMB) anisotropy in polarisation. We find that a departure from the light speed value would leave a characteristic imprint on the BB spectrum part induced by tensors, manifesting as a shift in the angular scale of its peaks. We derive constraints by using the available {\it Planck} and BICEP2 datasets showing how $c_T^2$ can be measured, albeit obtaining weak constraints due to the overall poor accuracy of the current BB power spectrum measurements. The present constraint corresponds to $c_T^2 = 1.30 \pm 0.79$ and $c_T^2< 2.85$ at $95%$ C.L. by assuming a power law primordial tensor power spectrum and $c_T^2<2.33$ at $95%$ C.L. if the running of the spectral index is allowed. We derive forecasts for the next generation CMB satellites, which we find capable of tightly constraining $c_T^2$ at percent level, comparable with bounds from binary pulsar measurements, largely due to the absence of degeneracy with other cosmological parameters.

It has recently been shown that the graviton can consistently gain a constant mass without introducing the Boulware-Deser ghost. We propose a gravity model where the graviton mass is set by a scalar field and prove that this model is free of the Boulware-Deser ghost by analyzing its constraint system and showing that two constraints arise. We also initiate the study of the model's cosmic background evolution and tentatively discuss possible cosmological implications of this model. In particular, we consider a simple scenario where the scalar field setting the graviton mass is identified with the inflaton and the graviton mass evolves from a high to a low energy scale, giving rise to the current cosmic acceleration.

Oscillons are long-lived, localized excitations of nonlinear scalar fields which may be copiously produced during preheating after inflation, leading to a possible oscillon-dominated phase in the early Universe. For example, this can happen after axion monodromy inflation, on which we run our simulations. We investigate the stochastic gravitational wave background associated with an oscillon-dominated phase. An isolated oscillon is spherically symmetric and does not radiate gravitational waves, and we show that the flux of gravitational radiation generated between oscillons is also small. However, a significant stochastic gravitational wave background may be generated during preheating itself (i.e, when oscillons are forming), and in this case the characteristic size of the oscillons is imprinted on the gravitational wave power spectrum, which has multiple, distinct peaks.

In Phys. Rev. Lett. 109 (2012) 231301, van den Aarssen et al. attempt to solve all the small-scale structure problems in the cold-dark-matter scenario by introducing interaction between the dark-matter particle and an MeV-mass vector boson. In addition, the latter is coupled to ordinary neutrinos, and the neutrino-vector boson coupling constant 10^{-5} < g_nu < 0.1 is required for that purpose. We show that this dark-matter scenario is subject to restrictive constraints from Big Bang Nucleosynthesis.

We extend the four-dimensional de Rham-Gabadadze-Tolley (dRGT) massive gravity model to a general scalar massive-tensor theory in arbitrary dimensions, coupling a dRGT massive graviton to multiple scalars and allowing for generic kinetic and mass matrix mixing between the massive graviton and the scalars, and derive its Hamiltonian formulation and associated constraint system. When passing to the Hamiltonian formulation, two different sectors arise: a general sector and a special sector. Although obtained via different ways, there are two second class constraints in either of the two sectors, eliminating the BD ghost. However, for the special sector, there are still ghost instabilities except for the case of two dimensions. In particular, for the special sector with one scalar, there is a "second BD ghost".

We develop the gradient expansion formalism for shift-symmetric Galileon-type actions. We focus on backgrounds that undergo inflation, work in the synchronous gauge, and obtain a general solution up to second order without imposing extra conditions at first order. The solution simplifies during the late stages of inflation. We also define a curvature perturbation conserved up to first order.

We continue to introduce bi-galileon theory, the generalisation of the single galileon model introduced by Nicolis et al. The theory contains two coupled scalar fields and is described by a Lagrangian that is invariant under Galilean shifts in those fields. This paper is the second of two, and focuses on the phenomenology of the theory. We are particularly interesting in models that admit solutions that are asymptotically self accelerating or asymptotically self tuning. In contrast to the single galileon theories, we find examples of self accelerating models that are simultaneously free from ghosts, tachyons and tadpoles, able to pass solar system constraints through Vainshtein screening, and do not suffer from problems with superluminality, Cerenkov emission or strong coupling. We also find self tuning models and discuss how Weinberg's no go theorem is evaded by breaking Poincar\'e invariance in the scalar sector. Whereas the galileon description is valid all the way down to solar system scales for the self-accelerating models, unfortunately the same cannot be said for self tuning models owing to the scalars backreacting strongly on to the geometry.

Galileon models arise in certain braneworld scenarios as modifications to General Relativity, and are also interesting as field theories in their own right. We show how the galileon model can be naturally generalized to include local gauge symmetries, by allowing for couplings to Yang-Mills fields. The resulting theories have at most second order spacetime derivatives in any representation of the gauge group, thereby avoiding Ostrogradski ghosts. We also extend the models to include curved space, and show how in that case we need to include non-minimal couplings between the galileons and the curvature tensors for the theory to retain its second order nature.

We re-examine the interesting possibility of utilizing the Glashow resonance (GR) channel nu_ebar + e^- to W^- to anything to discriminate between the UHE cosmic neutrinos originating from p-gamma and pp collisions in an optically thin source of cosmic rays. We propose a general parametrization of the initial neutrino flavor composition by allowing the ratios Phi^{p gamma}_{pi^-}/Phi^{p gamma}_{pi^+} and Phi^{pp}_{pi^-}/Phi^{pp}_{pi^+} to slightly deviate from their conventional values. A relationship between the typical source parameter kappa = (Phi^{p gamma}_{pi^+} + Phi^{p gamma}_{pi^-})/(Phi^{pp}_{pi^+} + Phi^{pp}_{pi^-} + Phi^{p gamma}_{pi^+} + Phi^{p gamma}_{pi^-}) and the working observable of the GR R_0 = Phi^T_{nu_ebar}/ (Phi^T_{nu_mu} + Phi^T_{nu_mu}) at a neutrino telescope is derived, and the numerical dependence of R_0 on kappa is illustrated by taking account of the latest experimental data on three neutrino mixing angles. It is shown that a measurement of R_0 is in principle possible to identify the pure p-gamma interaction (kappa =1), the pure pp interaction (kappa =0) or a mixture of both of them (0 < kappa < 1) at a given source of UHE cosmic neutrinos. The event rate of the GR signal against the background is also estimated.

Active-sterile neutrino mixing is strongly constrained for m_s > 100 keV to avoid excessive energy losses from supernova cores. For smaller m_s, matter effects suppress the effective mixing angle except for a resonant range of energies where it is enhanced. We study the case of \nu_tau-\nu_s-mixing where a \nu_tau-\bar\nu_tau asymmetry builds up due to the strong excess of \nu_s over \bar\nu_s emission or vice versa, reducing the overall emission rate. In the warm dark matter range m_s < 10 keV the mixing angle is essentially unconstrained.

The galileon model was recently proposed to locally describe a class of modified gravity theories, including the braneworld DGP model. We discuss spontaneous symmetry breaking of the self-accelerating branch in a multi-galileon theory with internal global symmetries. We show a modified version of Goldstone's theorem is applicable to the symmetry breaking pattern and discuss its implications. We also derive the Hamiltonian of a general multi-galileon theory and discuss its implications.

Modified gravity theories with the Gauss-Bonnet term $G=R^2-4R^{\mu\nu}R_{\mu\nu}+R^{\mu\nu\rho\sigma}R_{\mu\nu\rho\sigma}$ have recently gained a lot of attention as a possible explanation of dark energy. We perform a thorough phase space analysis on the so-called $f(G)$ models, where $f(G)$ is some general function of the Gauss-Bonnet term, and derive conditions for the cosmological viability of $f(G)$ dark energy models. Following the $f(R)$ case, we show that these conditions can be nicely presented as geometrical constraints on the derivatives of $f(G)$. We find that for general $f(G)$ models there are two kinds of stable accelerated solutions, a de Sitter solution and a phantom-like solution. They co-exist with each other and which solution the universe evolves to depends on the initial conditions. Finally, several toy models of $f(G)$ dark energy are explored. Cosmologically viable trajectories that mimic the $\Lambda$CDM model in the radiation and matter dominated periods, but have distinctive signatures at late times, are obtained.

A measurement of the ultrahigh-energy (UHE) cosmic neutrinos at a km^3-size neutrino telescope will open a new window to constrain the 3 x 3 neutrino mixing matrix V and probe possible new physics. We point out that it is in principle possible to examine the non-unitarity of V, which is naturally expected in a class of seesaw models with one or more TeV-scale Majorana neutrinos, by using neutrino telescopes. Considering the UHE neutrinos produced from the decays of charged pions arising from proton-proton and (or) proton-gamma collisions at a distant astrophysical source, we show that their flavor ratios at a terrestrial neutrino telescope may deviate from the democratic flavor distribution \phi^T_e : \phi^T_\mu : \phi^T_\tau = 1 : 1 : 1 due to the seesaw-induced unitarity violation of V. Its effect can be as large as several percent and can serve for an illustration of how sensitive a neutrino telescope should be to this kind of new physics.

Reliable values of quark and lepton masses are important for model building at a fundamental energy scale, such as the Fermi scale M_Z \approx 91.2 GeV and the would-be GUT scale \Lambda_GUT \sim 2 \times 10^16 GeV. Using the latest data given by the Particle Data Group, we update the running quark and charged-lepton masses at a number of interesting energy scales below and above M_Z. In particular, we take into account the possible new physics scale (\mu \sim 1 TeV) to be explored by the LHC and the typical seesaw scales (\mu \sim 10^9 GeV and \mu \sim 10^12 GeV) which might be relevant to the generation of neutrino masses. For illustration, the running masses of three light Majorana neutrinos are also calculated. Our up-to-date table of running fermion masses are expected to be very useful for the study of flavor dynamics at various energy scales.

We take a new approach to construct Quintessential models. With this approach, we first easily obtain a tracker solution that is different from those discovered before and straightforwardly find a solution of multiple attractors, i.e., a solution with more than one attractor for a given set of parameters. Then we propose a scenario of Quintessence where the field jumps out of the scaling attractor to the de-Sitter-like attractor, by introducing a field whose value changes a certain amount in a short time, leading to the current acceleration. We also calculate the change the field needs for a successful jump and suggest a possible mechanism that involves spontaneous symmetry breaking to realize the sudden change of the field value.

Following Georgi's unparticle scheme, we examine the effective couplings between neutrinos and unparticle operators. As an immediate consequence, neutrinos become unstable and can decay into the unparticle stuff. Assuming the dimension transmutation scale is around $\Lambda^{}_{\cal U} \sim 1 ~{\rm TeV}$, we implement the cosmological limit on the neutrino lifetime to constrain the neutrino-unparticle couplings for different scaling dimensions $d$. In addition, provided that the electron-unparticle coupling is restricted due to the precise measurement of the anomalous magnetic moment of electron, we calculate the unparticle contribution to the neutrino-electron elastic scattering. It is more important to jointly deal with the couplings of the unparticle to the standard model particles rather than separately. Taking into account both electron- and neutrino-unparticle couplings, we find that the scaling dimension of the scalar unparticle should lie in the narrow range $1 < d < 2$ by requiring the observables to be physically meaningful. However, there is no consistent range of $d$ for the vector unparticle operator.

We present a review of neutrino phenomenology in the minimal seesaw model (MSM), an economical and intriguing extension of the Standard Model with only two heavy right-handed Majorana neutrinos. Given current neutrino oscillation data, the MSM can predict the neutrino mass spectrum and constrain the effective masses of the tritium beta decay and the neutrinoless double-beta decay. We outline five distinct schemes to parameterize the neutrino Yukawa-coupling matrix of the MSM. The lepton flavor mixing and baryogenesis via leptogenesis are investigated in some detail by taking account of possible texture zeros of the Dirac neutrino mass matrix. We derive an upper bound on the CP-violating asymmetry in the decay of the lighter right-handed Majorana neutrino. The effects of the renormalization-group evolution on the neutrino mixing parameters are analyzed, and the correlation between the CP-violating phenomena at low and high energies is highlighted. We show that the observed matter-antimatter asymmetry of the Universe can naturally be interpreted through the resonant leptogenesis mechanism at the TeV scale. The lepton-flavor-violating rare decays, such as $\mu \to e + \gamma$, are also discussed in the supersymmetric extension of the MSM.

We propose a simple but useful parametrization of the flavor composition of ultrahigh-energy neutrino fluxes produced from distant astrophysical sources: $\phi^{}_e : \phi^{}_\mu : \phi^{}_\tau = \sin^2 \xi \cos^2 \zeta : \cos^2 \xi \cos^2 \zeta : \sin^2 \zeta$. We show that it is possible to determine or constrain $\xi$ and $\zeta$ by observing two independent neutrino flux ratios at the second-generation neutrino telescopes, provided three neutrino mixing angles and the Dirac CP-violating phase have been well measured in neutrino oscillations. Any deviation of $\zeta$ from zero will signify the existence of cosmic $\nu^{}_\tau$ and $\bar{\nu}^{}_\tau$ neutrinos at the source, and an accurate value of $\xi$ can be used to test both the conventional mechanism and the postulated scenarios for cosmic neutrino production.

We have observed 23 Class 0 sources in the HCO+ J=4-3 and 3-2 lines. The mean bolometric temperature of the 16 sources with well-determined values is 44 K and the mean luminosity is 5.7 L_sun, excluding two sources of considerably higher luminosity. Nine sources, including three sources previously suggested to be collapsing, have the correct (blue) spectral line asymmetry for infall in both lines. Three sources have the opposite (red) asymmetry in both lines, and one source, L1157, has a red asymmetry in HCO+ J=4-3 and a blue asymmetry in J=3-2. The rest have no significant or consistent asymmetry. The H13CO+ J=4-3 and 3-2 lines were also observed to find the velocity of the ambient gas, and sources with an interesting line asymmetry were mapped. A Monte Carlo code was used to produce an evolutionary sequence of collapsing cloud models of the HCO+ J=4-3 and 3-2 lines and to compare various diagnostics of the resulting line profiles. The same code was used to compare infall models to the observations in one source, L1527. The results were consistent with previous collapse models. Based on integrated intensity maps of the line peaks and wings, as well as the velocity of the H13CO+ line, we select six of the nine sources with a blue line asymmetry as good candidates for protostellar collapse. Further evidence is needed to establish that infall is taking place. The HCO+ spectra are not conclusive because bipolar outflows produce strong emission, which can confuse the issue in any individual source. However, the predominance of blue asymmetries over red asymmetries is not naturally explained in outflow models, whereas it is expected in collapse models.

We have observed rovibrational transitions of acetylene and HCN near 13 microns in absorption toward GL2591. We also observed rotational lines of CS, HCN, H2CO, and HCO+. The combined data are analyzed in terms of models with a cloud envelope with density gradients and discrete regions of hot, dense gas, probably near the infrared source. The abundance of HCN is enhanced by a factor of 400 in the gas producing the infrared absorption, in agreement with chemical models which involve depletion of molecules onto grains and subsequent sublimation when temperatures are raised.