We present a Fisher information study of the statistical impact of galaxy bias and selection effects on the estimation of key cosmological parameters from galaxy redshift surveys; in particular, the angular diameter distance, Hubble parameter, and linear growth rate at a given redshift, the cold dark matter density, and the tilt and running of the primordial power spectrum. The line-of-sight-dependent selection contributions we include here are known to exist in real galaxy samples. We determine the maximum wavenumber included in the analysis by requiring that the next-order corrections to the galaxy power spectrum or bispectrum, treated here at next-to-leading and leading order, respectively, produce shifts of $\lesssim 0.25\sigma$ on each of the six cosmological parameters. With the galaxy power spectrum alone, selection effects can deteriorate the constraints severely, especially on the linear growth rate. Adding the galaxy bispectrum helps break parameter degeneracies significantly. We find that a joint power spectrum-bispectrum analysis of a Euclid-like survey can still measure the linear growth rate to 10% precision after complete marginalization over selection bias. We also discuss systematic parameter shifts arising from ignoring selection effects and/or other bias parameters, and emphasize that it is necessary to either control selection effects at the percent level or marginalize over them. We obtain similar results for the Roman Space Telescope and HETDEX.

Quasidilaton massive gravity offers a physically well-defined gravitational theory with non-zero graviton mass. We present the full set of dynamical equations governing the expansion history of the universe, valid during radiation domination, matter domination, and a late-time self-accelerating epoch related to the graviton mass. The existence of self-consistent solutions constrains the amplitude of the quasi-dilaton field and the graviton mass, as well as other model parameters. We point out that the effective mass of gravitational waves can be significantly larger than the graviton mass, opening the possibility that a single theory can explain both the late-time acceleration of the cosmic expansion and modifications of structure growth leading to the suppression of large-angle correlations observed in the cosmic microwave background.

We derive the evolution equation for the density matrix of a UV- and IR- limited band of comoving momentum modes of the canonically normalized scalar degree of freedom in two examples of nearly de Sitter universes. Including the effects of a cubic interaction term from the gravitational action and tracing out a set of longer wavelength modes, we find that the evolution of the system is non-Hamiltonian and non-Markovian. We find linear dissipation terms for a few modes with wavelength near the boundary between system and bath and nonlinear dissipation terms for all modes. The non-Hamiltonian terms persist to late times when the scalar field dynamics is such that the curvature perturbation continues to evolve on super-Hubble scales.

In this work, we report on the detection of enhanced TeV gamma- ray emission from the high synchrotron-peaked blazar Mrk 421 with the TACTIC telescope on the night of December 28, 2014 (MJD 57019). We use data from the TACTIC observations of Mrk 421 for one week during December 25-31, 2014 (MJD 57016-57022) in this study. The TACTIC observation on December 28, 2014 (MJD 57019) alone results in the detection of 86+-17 gamma- ray like events from Mrk 421 with a statistical significance of 5.17sigma in a livetime of ~ 2.2 hours above an energy threshold of 0.85 TeV. The high statistics (higher than three Crab Units) of TeV photons enables us to study the very high energy (VHE) gamma- ray emission from the source at shorter timescales. A minimum variability timescale of ~ 0.72 days is obtained for the TeV gamma- ray emission from Mrk 421 during the above flaring activity of the source. We have also analyzed the contemporaneous data from Fermi-LAT to study the high energy (HE) gamma- ray emission during the high activity state of the source. The HE gamma- ray emission is observed to be increasing after the TeV flaring activity detected with the TACTIC. We also use other near simultaneous archival data available from the Swift-BAT in hard X-rays and from SPOL at Steward Observatory in optical V and R bands to characterize the multi-wavelength emission of Mrk 421 during the high activity state observed at TeV energies. The TeV gamma- ray emission observed on December 28, 2014 (MJD 57019) is found to be more prominant than the emissions in lower energy bands during the same period. The TeV gamma- ray observation of Mrk 421 in high activity state with the TACTIC telescope is also used to understand the physical mechanism for blazar emission under the frame work of the leptonic single zone synchrotron self Compton process.

We have observed Markarian 501 in Very High Energy (VHE) gamma-ray wavelength band for 70.6 hours from 15 April to 30 May, 2012 using TACTIC telescope. Detailed analysis of $\sim$66.3 hours of clean data revealed the presence of a TeV $\gamma$-ray signal (686$\pm$77 $\gamma$-ray events) from the source direction with a statistical significance of 8.89$\sigma$ above 850 GeV. Further, a total of 375 $\pm$ 47 $\gamma$-ray like events were detected in 25.2 hours of observation from 22 - 27 May, 2012 with a statistical significance of 8.05$\sigma$ indicating that the source has possibly switched over to a relatively high gamma-ray emission state. We have derived time-averaged differential energy spectrum of the state in the energy range 850 GeV - 17.24 TeV which fits well with a power law function of the form $dF/dE=f_{0}E^{-\Gamma}$ with $f_{0}= (2.27 \pm 0.38) \times 10 ^{-11} $ photons cm$^{-2}$ s$^{-1}$ TeV$^{-1}$ and $\Gamma=2.57 \pm 0.15$. In order to investigate the source state, we have also used almost simultaneous multiwavelength observations viz: high energy data collected by $\it{Fermi}$-LAT, X-ray data collected by $\it{Swift}$-XRT and MAXI, optical and UV data collected by $\it{Swift}$-UVOT, and radio data collected by OVRO, and reconstructed broad-band Spectral Energy Distribution (SED). The obtained SED supports leptonic model (homogeneous single zone) for VHE gamma-ray emission involving synchrotron and synchrotron self Compton (SSC) processes.

Strong X-ray and $\gamma$--ray flares have been detected in February 2010 from the high synchrotron peaked blazar Mrk 421 (z=0.031). With the motivation of understanding the physics involved in this flaring activity, we study the variability of the source in X-ray and $\gamma$--ray energy bands during the period February 10-23, 2010 (MJD 55237-55250). We use near simultaneous X-ray data collected by \emph{MAXI}, \emph{Swift}-XRT and $\gamma$--ray data collected by \emph{Fermi}-LAT and \emph{TACTIC} along with the optical V-band observations by \emph{SPOL} at Steward Observatory. We observe that the variation in the one day averaged flux from the source during the flare is characterized by fast rise and slow decay. Besides, the TeV $\gamma$--ray flux shows a strong correlation with the X-ray flux, suggesting the former to be an outcome of synchrotron self Compton emission process. To model the observed X-ray and $\gamma$--ray light curves, we numerically solve the kinetic equation describing the evolution of particle distribution in the emission region. The injection of particle distribution into the emission region, from the putative acceleration region, is assumed to be a time dependent power law. The synchrotron and synchrotron self Compton emission from the evolving particle distribution in the emission region are used to reproduce the X-ray and $\gamma$--ray flares successfully. Our study suggests that the flaring activity of Mrk 421 can be an outcome of an efficient acceleration process associated with the increase in underlying non-thermal particle distribution.

Non-Gaussianity in the distribution of inflationary perturbations, measurable in statistics of the cosmic microwave background (CMB) and large scale structure fluctuations, can be used to probe non-trivial initial quantum states for these perturbations. The bispectrum shapes predicted for generic non-Bunch-Davies initial states are non-factorizable ("non-separable") and are highly oscillatory functions of the three constituent wavenumbers. This can make the computation of CMB bispectra, in particular, computationally intractable. To efficiently compare with CMB data one needs to construct a separable template that has a significant similarity with the actual shape in momentum space. In this paper we consider a variety of inflationary scenarios, with different non-standard initial conditions, and how best to construct viable template matches. In addition to implementing commonly used separable polynomial and Fourier bases, we introduce a basis of localized piecewise spline functions. The spline basis is naturally nearly orthogonal, making it easy to implement and to extend to many modes. We show that, in comparison to existing techniques, the spline basis can provide better fits to the true bispectrum, as measured by the cosine between shapes, for sectors of the theory space of general initial states. As such, it offers a useful approach to investigate non-trivial features generated by fundamental properties of the inflationary Universe.

The Sloan Digital Sky Survey has surveyed 14,555 square degrees of the sky, and delivered over a trillion pixels of imaging data. We present the large-scale clustering of 1.6 million quasars between z = 0.5 and z = 2.5 that have been classified from this imaging, representing the highest density of quasars ever studied for clustering measurements. This data set spans ~11,000 square degrees and probes a volume of 80(Gpc/h)^3. In principle, such a large volume and medium density of tracers should facilitate high-precision cosmological constraints. We measure the angular clustering of photometrically classified quasars using an optimal quadratic estimator in four redshift slices with an accuracy of ~25% over a bin width of l ~10 - 15 on scales corresponding to matter-radiation equality and larger (l ~ 2 - 30). Observational systematics can strongly bias clustering measurements on large scales, which can mimic cosmologically relevant signals such as deviations from Gaussianity in the spectrum of primordial perturbations. We account for systematics by employing a new method recently proposed by Agarwal et al. (2014) to the clustering of photometrically classified quasars. We carefully apply our methodology to mitigate known observational systematics and further remove angular bins that are contaminated by unknown systematics. Combining quasar data with the photometric luminous red galaxy (LRG) sample of Ross et al. (2011) and Ho et al. (2012), and marginalizing over all bias and shot noise-like parameters, we obtain a constraint on local primordial non-Gaussianity of fNL = -113+/-154 (1\sigma error). [Abridged]

In this paper, we present results of TeV $\gamma$--ray observations of the high synchrotron peaked BL Lac object 1ES 1218+304 (z=0.182) with the $TACTIC$ (TeV Atmospheric Cherenkov Telescope with Imaging Camera). The observations are primarily motivated by the unusually hard GeV-TeV spectrum of the source despite its relatively large redshift. The source is observed in the TeV energy range with the $TACTIC$ from March 1, 2013 to April 15, 2013 (MJD 56352--56397) for a total observation time of 39.62 h and no evidence of TeV $\gamma$--ray activity is found from the source. The corresponding 99$\%$ confidence level upper limit on the integral flux above a threshold energy of 1.1 TeV is estimated to be 3.41 $\times10^{-12}$ photons cm$^{-2}$ s$^{-1}$ (i.e $<23\%$ Crab Nebula flux) assuming a power law differential energy spectrum with photon index 3.0, as previously observed by the $MAGIC$ and $VERITAS$ telescopes. For the study of multi-wavelength emission from the source, we use nearly simultaneous optical, UV and and X--ray data collected by the UVOT and XRT instruments on board the \emph{Swift} satellite and high energy $\gamma$--ray data collected by the Large Area Telescope on board the \emph{Fermi} satellite. We also use radio data at 15 GHz from OVRO 40 m telescope in the same period. No significant increase of activity is detected from radio to TeV $\gamma$--rays from 1ES1218+304 during the period from March 1, 2013 to April 15, 2013.

The effective field theory (EFT) of inflation has become the preferred method for computing cosmological correlation functions of the curvature fluctuation, $\zeta$. It makes explicit use of the soft breaking of time diffeomorphisms by the inflationary background to organize the operators expansion in the action of the Goldstone mode $\pi$ associated with this breaking. Despite its ascendancy, there is another method for calculating $\zeta$ correlators, involving the direct calculation of the so-called Horndeski action order by order in powers of $\zeta$ and its derivatives. The question we address in this work is whether or not the $\zeta$ correlators calculated in these seemingly different ways are in fact the same. The answer is that the actions to cubic order in either set of variables do indeed give rise to the same $\zeta$ bispectra, but that to make this equivalence manifest requires a careful understanding of the non-linear transformations relating $\pi$ to $\zeta$ and how boundary terms in the actions are affected by imposing this relation. As a by product of our study we find that the calculations in the $\pi$ language can be simplified considerably in a way that allows us to use only the linear part of the $\pi-\zeta$ relation simply by changing the coefficients of some of the operators in the EFT. We also note that a proper accounting of the boundary terms will be of the greatest importance when computing the bispectrum for more general initial states than the Bunch-Davies one.

Photometric large scale structure (LSS) surveys probe the largest volumes in the Universe, but are inevitably limited by systematic uncertainties. Imperfect photometric calibration leads to biases in our measurements of the density fields of LSS tracers such as galaxies and quasars, and as a result in cosmological parameter estimation. Earlier studies have proposed using cross-correlations between different redshift slices or cross-correlations between different surveys to reduce the effects of such systematics. In this paper we develop a method to characterize unknown systematics. We demonstrate that while we do not have sufficient information to correct for unknown systematics in the data, we can obtain an estimate of their magnitude. We define a parameter to estimate contamination from unknown systematics using cross-correlations between different redshift slices and propose discarding bins in the angular power spectrum that lie outside a certain contamination tolerance level. We show that this method improves estimates of the bias using simulated data and further apply it to photometric luminous red galaxies in the Sloan Digital Sky Survey as a case study.

Primordial non-Gaussianity induces a scale-dependent bias in large scale structure (LSS) data, proportional to $f_{\rm NL}/k^2$ for the exact local ansatz. Recent work has shown that models of inflation that predict a large squeezed limit bispectrum, such as multi-field models and single field inflation with a modified initial state, typically give rise to a generalized local ansatz, with the scale-dependent bias now proportional to ${\cal A}_{\rm NL}/k^\alpha$. We use photometric measurements of the angular power spectrum of luminous red galaxies and quasars in the Sloan Digital Sky Survey Data Release Eight (SDSS DR8) with the above parameterization to constrain the amplitude ${\cal A}_{\rm NL}$ and scale-dependence $\alpha$. We find that the marginalized upper limit on $\alpha$ is 2.0 at the $95\%$ confidence level, consistent with the local ansatz. We also present Fisher forecasts for a survey of the same size as DR8 to assess the role of systematics in current photometric LSS data. Moreover, we present analytic results on the expected mass dependence of ${\cal A}_{\rm NL}$ for different inflationary models, which can be an important observable for future surveys, if primordial non-Gaussianity is non-zero.

We argue that a plausible operational definition for an initial state of the Universe is the initial quantum state of the curvature perturbations generated during inflation. We provide a parameterization of this state and generalize the standard in-in formalism to incorporate the structures in this state into the computation of correlators of the perturbations. Measurements of these correlators using both the CMB as well as large scale structure probe different structures in the initial state, as they give rise to bi- and tri-spectra peaked on different shapes of triangles and quadrilaterals in momentum space. In essence, the shapes implied by the correlators feed directly into information about the shape of the initial state and what physics could have preceded inflation to set this state up.

We study the effects of non-trivial initial quantum states for inflationary fluctuations within the context of the effective field theory for inflation constructed by Cheung et al. which allows us to discriminate between different initial states in a model-independent way. We develop a Green's function/path integral based formulation that incorporates initial state effects and use it to address questions such as how state-dependent is the consistency relation for the bispectrum, how many e-folds beyond the minimum required to solve the cosmological fine tunings of the big bang are we allowed so that some information from the initial state survives until late times, among others. We find that the so-called consistency condition relating the local limit of the bispectrum and the slow-roll parameter is a state-dependent statement that can be avoided for physically consistent initial states either with or without initial non-Gaussianities.

We present the results of three-dimensional simulations of quasar polarizations in the presence of pseudoscalar-photon mixing in the intergalactic medium. The intergalactic magnetic field is assumed to be uncorrelated in wave vector space but correlated in real space. Such a field may be obtained if its origin is primordial. Furthermore we assume that the quasars, located at cosmological distances, have negligible initial polarization. In the presence of pseudoscalar-photon mixing we show, through a direct comparison with observations, that this may explain the observed large scale alignments in quasar polarizations within the framework of big bang cosmology. We find that the simulation results give a reasonably good fit to the observed data.

We have observed the blazar Markarian 421 with the TACTIC $\gamma$-ray telescope at Mt. Abu, India, from 22 November 2009 to 16 May 2010 for 265 hours. Detailed analysis of the data so recorded revealed presence of a TeV $\gamma$-ray signal with a statistical significance of 12.12$\sigma$ at $E_{\gamma}\geq$ 1 TeV. We have estimated the time averaged differential energy spectrum of the source in the energy range 1.0 - 16.44 TeV. The spectrum fits well with the power law function of the form ($dF/dE=f_0 E^{-\Gamma}$) with $f_0=(1.39\pm0.239)\times 10^{-11}cm^{-2}s^{-1}TeV^{-1}$ and $\Gamma=2.31\pm0.14$.

We investigate the effects of pseudoscalar-photon mixing on electromagnetic radiation in the presence of correlated extragalactic magnetic fields. We model the Universe as a collection of magnetic domains and study the propagation of radiation through them. This leads to correlations between Stokes parameters over large scales and consistently explains the observed large-scale alignment of quasar polarizations at different redshifts within the framework of the big bang model.

We study the six-field dynamics of D3-brane inflation for a general scalar potential on the conifold, finding simple, universal behavior. We numerically evolve the equations of motion for an ensemble of more than 7 \times 10^7 realizations, drawing the coefficients in the scalar potential from statistical distributions whose detailed properties have demonstrably small effects on our results. When prolonged inflation occurs, it has a characteristic form: the D3-brane initially moves rapidly in the angular directions, spirals down to an inflection point in the potential, and settles into single-field inflation. The probability of N_{e} e-folds of inflation is a power law, P(N_{e}) \propto N_{e}^{-3}, and we derive the same exponent from a simple analytical model. The success of inflation is relatively insensitive to the initial conditions: we find attractor behavior in the angular directions, and the D3-brane can begin far above the inflection point without overshooting. In favorable regions of the parameter space, models yielding 60 e-folds of expansion arise approximately once in 10^3 trials. Realizations that are effectively single-field and give rise to a primordial spectrum of fluctuations consistent with WMAP, for which at least 120 e-folds are required, arise approximately once in 10^5 trials. The emergence of robust predictions from a six-field potential with hundreds of terms invites an analytic approach to multifield inflation.

We study a flat brane solution in an effective 5D action for cascading gravity and propose a mechanism to screen extrinsic curvature in the presence of a large tension on the brane. The screening mechanism leaves the bulk Riemann-flat, thus making it simpler to generalize large extra dimension dark energy models to higher codimensions. By studying an action with cubic interactions for the brane-bending scalar mode, we find that the perturbed action suffers from ghostlike instabilities for positive tension, whereas it can be made ghost-free for sufficiently small negative tension.

We develop a fully covariant, well-posed 5D effective action for the 6D cascading gravity brane-world model, and use this to study cosmological solutions. We obtain this effective action through the 6D decoupling limit, in which an additional scalar degree mode, \pi, called the brane-bending mode, determines the bulk-brane gravitational interaction. The 5D action obtained this way inherits from the sixth dimension an extra \pi self-interaction kinetic term. We compute appropriate boundary terms, to supplement the 5D action, and hence derive fully covariant junction conditions and the 5D Einstein field equations. Using these, we derive the cosmological evolution induced on a 3-brane moving in a static bulk. We study the strong- and weak-coupling regimes analytically in this static ansatz, and perform a complete numerical analysis of our solution. Although the cascading model can generate an accelerating solution in which the \pi field comes to dominate at late times, the presence of a critical singularity prevents the \pi field from dominating entirely. Our results open up the interesting possibility that a more general treatment of degravitation in a time-dependent bulk, or taking into account finite brane-thickness effects, may lead to an accelerating universe without a cosmological constant.

Inflation is now an accepted paradigm in standard cosmology, with its predictions consistent with observations of the cosmic microwave background. It lacks, however, a firm physical theory, with many possible theoretical origins beyond the simplest, canonical, slow-roll inflation, including Dirac-Born-Infeld inflation and k-inflation. We discuss how a hierarchy of Hubble flow parameters, extended to include the evolution of the inflationary sound speed, can be applied to compare a general, single field inflationary action with cosmological observational data. We show that it is important to calculate the precise scalar and tensor primordial power spectra by integrating the full flow and perturbation equations, since values of observables can deviate appreciably from those obtained using typical second-order Taylor expanded approximations in flow parameters. As part of this, we find that a commonly applied approximation for the tensor to scalar ratio, r \approx 16 c_s\epsilon, becomes poor (deviating by as much as 50%) as c_s deviates from 1 and hence the Taylor expansion including next-to-leading order contribution terms involving c_s is required. By integrating the full flow equations, we use a Monte-Carlo-Markov-Chain approach to impose constraints on the parameter space of general single field inflation, and reconstruct the properties of such an underlying theory in light of recent cosmic microwave background and large-scale structure observations.

We model the effect of photon and ultra-light pseudoscalar mixing on the propagation of electromagnetic radiation through the extragalactic medium. The medium is modelled as a large number of magnetic domains, uncorrelated with one another. We obtain an analytic expression for the different Stokes parameters in the limit of small mixing angle. The different Stokes parameters are found to increase linearly with the number of domains. We also verify this result by direct numerical simulations. We use this formalism to estimate the effect of pseudoscalar-photon mixing on the Cosmic Microwave Background (CMB) polarization. We impose limits on the model parameters by the CMB observations. We find that the currently allowed parameter range admits a CMB circular polarization up to order $10^{-7}$.

We establish the dynamical attractor behavior in scalar-tensor theories of dark energy, providing a powerful framework to analyze classes of theories, predicting common evolutionary characteristics that can be compared against cosmological constraints. In the Jordan frame the theories are viewed as a coupling between a scalar field, \Phi, and the Ricci scalar, R, F(\Phi)R. The Jordan frame evolution is described in terms of dynamical variables m \equiv d\ln F/d\ln \Phi and r \equiv -\Phi F/f, where F(\Phi) = d f(\Phi)/d\Phi. The evolution can be alternatively viewed in the Einstein frame as a general coupling between scalar dark energy and matter, \beta. We present a complete, consistent picture of evolution in the Einstein and Jordan frames and consider the conditions on the form of the coupling F and \beta required to give the observed cold dark matter (CDM) dominated era that transitions into a late time accelerative phase, including transitory accelerative eras that have not previously been investigated. We find five classes of evolutionary behavior of which four are qualitatively similar to those for f(R) theories (which have \beta=1/2). The fifth class exists only for |\beta| < \sqrt{3}/4, i.e. not for f(R) theories. In models giving transitory late time acceleration, we find a viable accelerative region of the (r,m) plane accessible to scalar-tensor theories with any coupling, \beta (at least in the range |\beta| \leq 1/2, which we study in detail), and an additional region open only to theories with |\beta| < \sqrt{3}/4.

The results of our observational campaigns on the two extragalactic sources Mkn 501 and Mkn 421, carried out with the Imaging Element of the TACTIC array, during March-May, 1998 and April-May, 1999, are presented. The results indicate that the two BL Lac objects (Mkn 501 and Mkn 421) were in a `low' gamma-ray emission state during both epochs of our observations.