Blue straggler stars (BSS) are well studied in globular clusters but their systematic study with secure membership determination is lacking in open clusters. We use Gaia DR2 data to determine accurate stellar membership for four intermediate-age open clusters, Melotte 66, NGC 2158, NGC 2506 and NGC 6819, and three old open clusters, Berkeley 39, NGC 188 and NGC 6791, to subsequently study their BSS populations. The BSS radial distributions of five clusters, Melotte 66, NGC 188, NGC 2158, NGC 2506, and NGC 6791, show bimodal distributions, placing them with Family II globular clusters which are of intermediate dynamical ages. The location of minima, $r_\mathrm{{min}}$, in the bimodal BSS radial distributions, varies from 1.5$r_c$ to 4.0$r_c$, where $r_c$ is the core radius of the clusters. We find a positive correlation between $r_\mathrm{{min}}$ and $N_{\mathrm{relax}}$, the ratio of cluster age to the current central relaxation time of the cluster. We further report that this correlation is consistent in its slope, within the errors, to the slope of the globular cluster correlation between the same quantities, but with a slightly higher intercept. This is the first example in open clusters that shows BSS radial distributions as efficient probes of dynamical age. The BSS radial distributions of the remaining two clusters, Berkeley 39 and NGC 6819, are flat. The estimated $N_{\mathrm{relax}}$ values of these two clusters, however, indicate that they are dynamically evolved. Berkeley 39 especially has its entire BSS population completely segregated to the inner regions of the cluster.

Within Type-I seesaw mechanism, Higgs mass can be dynamically generated via quantum effects of the right handed neutrinos assuming the potential is nearly conformal at the Ultra-Violet. The scenario, named as the "Neutrino Option" allows RH neutrino mass scale up to $M \lesssim$ $10^7$ GeV to be consistent with light neutrino masses, mixing and Higgs mass. Therefore, it is not consistent with standard hierarchical thermal leptogenesis. Parameter space for thermal resonant leptogenesis is highly constrained in this model. We point out that non-thermal pair production of RH neutrinos from inflaton decay allows non-resonant leptogenesis and allows RH mass scale to be smaller by an order of magnitude than what is obtained in the thermal case. Within the similar parameter space of thermal leptogenesis, RH neutrinos can also be produced from inflaton decay along with a Dark Matter having mass $M_{\rm DM}\lesssim$ 320 MeV. The main constraint in the latter scenario comes from the Ly$\alpha$ constraints on Dark Matter free streaming. We show in addition, that the Neutrino Option introduces a 'phantom window' for the RH mass scale, in which, contrary to the usual scenarios, CP asymmetry parameter for leptogenesis decreases with the increase of the RH mass scale and minimally fine-tuned seesaw models naturally exhibit this `phantom window'.

We perform a unique direct numerical simulation of Euler turbulence using delta-correlated velocity field as an initial condition, and report a full range of $k^2$ and $k$ energy spectra for 3D and 2D flows respectively, zero energy flux, and Maxwell-Boltzmann distribution for the velocity field. These are direct verification of the predictions of the absolute equilibrium theory of turbulence. For a coherent vortex as an initial condition, Euler turbulence transitions from a mixture of nonequilibrium-equilibrium state to a equilibrium state through a process called thermalization. In this letter, we present a model for thermalization in Euler turbulence.

A generic feature of models of inflation obtained from string compactifications is correlations between the model parameters and the post-inflationary evolution of the universe. Thus, the post-inflationary evolution depends on the inflationary model parameters and accurate inflationary predictions require that this be incorporated in the evolution of the primordial spectrum. The fibre inflation model is a promising model of inflation constructed in type IIB string theory. This model has two interesting features in its post-inflationary evolution. The reheating temperature of the model is directly correlated with the model parameters. The model also necessarily predicts some dark radiation, which can be sizable for certain choices of discrete parameters in the model. We analyse this model in detail using publicly available codes - ModeChord and CosmoMC with the latest Planck+BICEP2/Keck array data to constrain the model parameters and $N_{\rm pivot}$ (the number of e-foldings between horizon exit of CMB modes and the end of inflation). Our analysis sets up the basic methods necessary to extract precise inflationary prediction in string models incorporating correlations between model parameters and post-inflationary evolution.

The leading order effect of non-sphericity on the maximum size $R_{\rm TA,max}$ of realistic large scale bound cosmic structures is studied in the framework of $\Lambda{\rm CDM}$. Three simple models are used, two to describe intrinsically axisymmetric static mass distributions and the Kerr-de Sitter metric for stationary ones with non-vanishing angular momentum. In all cases the fractional change $\delta R_{\rm TA,max}(\theta)/R^{(0)}_{\rm TA,max}$ of $R_{\rm TA,max}$ of a given structure, compared to a spherical one with the same mass $M$, depends on the polar angle $\theta$ and is proportional to the product of the relevant eccentricity parameter, times the square of a small quantity, which in the static examples is the ratio of two characteristic length scales, while in the spinning case it is the ratio $v_{\rm out}/c$ of the azimuthal speed of the outmost members of the structure over the speed of light. Furthermore, the anglular average $\langle \delta R_{\rm TA,max}(\theta)/R^{(0)}_{\rm TA,max}\rangle$ is zero in the two static cases, while it is negative and proportional to ${\cal O}(v^2_{\rm out}/c^2)$ in the Kerr-de Sitter example. Thus, $\delta R_{\rm TA,max}(\theta)/R^{(0)}_{\rm TA,max}$ of axisymmetric structures is very small for practically any value of the eccentricity parameter. Based on the above, we argue that the maximum turnaround size of {\it realistic} superclusters has negligible dependence on their shapes.

For primordial black holes (PBHs) to form a considerable fraction of cold dark matter, the required amplitude of primordial scalar perturbations is quite large ($P_{\zeta}(k) \sim 10^{-2}$) if PBH is formed in radiation epoch. In alternate cosmological histories, where additional epochs of arbitrary equation of state precede radiation epoch, the dynamics of PBH formation and relevant mass ranges can be different leading to requirement of lower primordial power at smaller scales of inflation. Moreover, this alternate history can modify the predictions for the gravitational wave (GW) spectrum, which can be probed by upcoming GW observations. In this paper we show that an early kination epoch can lead to percent level abundance of PBH for a lower amplitude of $P_{\zeta}(k)$ as compared to PBH formation in a standard radiation epoch. Moreover, we calculate the effect of early kination epoch on the GW spectrum for first and second orders in perturbation which show enhancement in the amplitude of the GW spectrum in a kination epoch with respect to that in a standard radiation epoch.

Recent observations of the cosmic microwave background (CMB) indicate that a successful theory of cosmological inflation needs to have flat potential of the inflaton scalar field. Realizing the inflaton to be a pseudo-Nambu Goldstone boson (pNGB) could ensure the flatness and the sub-Planckian scales related to the dynamics of the paradigm. In this work, we have taken the most general form of such a scenario: Goldstone inflation, and studied the model in the noncanonical domain. Natural inflation is a limiting case of this model, which is also studied here in the noncanonical regime. Our result is compared with the recent release by Planck collaboration and it is shown that for some combination of the model parameters, a Goldstone inflationary model in the noncanonical realisation obeys the current observational bounds. Then, we studied the era of reheating after the end of inflation. For different choice of model parameters, constraints on the reheating temperature ($T_{\rm re}$) and number of e-folds during reheating($N_{\rm re}$) for the allowed inflationary observables (e.g. scalar spectral index($n_s$) and tensor to scalar ratio($r$)) are predicted for this model.

We analyse the epoch of reheating after an inflationary phase in the Randal Sundrum(RS) Type-$\rm II$ braneworld, where we did not consider any particular model of inflation, but rather reconstructed the inflationary potential solving the flow equations using Monte Carlo (MC) approach. It is shown numerically that a potential conceived through the MC reconstruction technique can be represented by an effective potential as a function of the number of e-foldings($N$). Then, the epoch of reheating is studied for this reconstructed potential. The relation between the reheating temperature ($T_{\rm reh}$) and the 5-dimensional Planck mass($M_5$) is established. Moreover, it is argued that there is a stringent bound on the critical reheating temperature that also translates to a tight bound on $M_5$ .

The age-velocity dispersion relation is an important tool to understand the evolution of the disc of the Andromeda galaxy (M31) in comparison with the Milky Way. We use Planetary Nebulae (PNe) to obtain the age-velocity dispersion relation in different radial bins of the M31 disc. We separate the observed PNe sample based on their extinction values into two distinct age populations. The observed velocities of our high- and low-extinction PNe, which correspond to higher and lower mass progenitors respectively, are fitted in de-projected elliptical bins to obtain their rotational velocities, $V_{\phi}$, and corresponding dispersions, $\rm\sigma_{\phi}$. We assign ages to the two PNe populations by comparing central-star properties of an archival sub-sample of PNe, having models fitted to their observed spectral features, to stellar evolution tracks. For the high- and low-extinction PNe, we find ages of $\sim2.5$ Gyr and $\sim4.5$ Gyr respectively, with distinct kinematics beyond a deprojected radius R$\rm_{GC}= 14$ kpc. At R$\rm_{GC}$=17--20 kpc, which is the equivalent distance in disc scale lengths of the Sun in the Milky Way disc, we obtain $\rm\sigma_{\phi,~2.5~Gyr}= 61\pm 14$ km s$^{-1}$ and $\rm\sigma_{\phi,~4.5~Gyr}= 101\pm 13$ km s$^{-1}$. The age-velocity dispersion relation for the M31 disc is obtained in two radial bins, R$\rm_{GC}$=14--17 and 17--20 kpc. The high- and low-extinction PNe are associated with the young thin and old thicker disc of M31 respectively, whose velocity dispersion values increase with age. These values are almost twice and thrice that of the Milky Way disc stellar population of corresponding ages. From comparison with simulations of merging galaxies, we find that the age-velocity dispersion relation in the M31 disc measured using PNe is indicative of a single major merger that occurred 2.5 -- 4.5 Gyr ago with an estimated merger mass ratio $\approx$ 1:5.

Hadronic inhomogeneities are formed after the quark hadron phase transition. The nature of the phase transition dictates the nature of the inhomogeneities formed. Recently some scenarios of inhomogeneities have been discussed where the strange quarks are in excess over the up and down quarks. The hadronization of these quarks will give rise to a large density of hyperons and kaons in addition to the protons and neutrons which are formed after the phase transition. These unstable hyperons decay into pions, muons and their respective neutrinos. Hence the plasma during this period consists of neutrons, protons, electrons, muons and neutrinos. Due to the decay of the hyperons, the muon component of the inhomogeneities will be very high. We study the diffusion of neutrons and protons in the presence of a large number of muons immediately after the quark hadron phase transition. We find that the presence of the muons enhances the diffusion coefficient of the neutrons/protons. As the diffusion coefficient is enhanced, the inhomogeneities will decay faster in the regions where the muon density is higher. Hence smaller muon rich inhomogeneities will be completely wiped out. The decay of the hyperons will also generate muon neutrinos. Since the big bang nucleosynthesis provides constraints on the neutrino degeneracies, we revisit the effect of non zero degeneracies on the primordial elements.

We propose a neutrino mass model with $\mu\tau$-flavored CP symmetry, where the effective light neutrino Lagrangian enjoys an additional invariance under a Friedberg-Lee (FL) transformation on the left-handed flavor neutrino fields, that leads to a highly predictive and testable scenario. While both types of the light neutrino mass ordering, i.e., Normal Ordering (NO) as well as the Inverted Ordering (IO) are allowed, the absolute scale of neutrino masses is fixed by the vanishing determinant of light Majorana neutrino mass matrix $M_\nu$. We show that for both types of mass ordering, whilst the atmospheric mixing angle $\theta_{23}$ is in general nonmaximal ($\theta_{23}\neq \pi/4$), the Dirac CP phase $\delta$ is exactly maximal ($\delta=\pi/2,3\pi/2$) for IO and nearly maximal for NO owing to $\cos\delta\propto \sin\theta_{13}$. For the NO, very tiny nonvanishing Majorana CP violation might appear through one of the Majorana phases $\beta$; otherwise the model predicts vanishing Majorana CP violation. Thus, despite the fact, that from the measurement of $\theta_{23}$, it is difficult to rule out the model, any large deviation of $\delta$ from its maximality, will surely falsify the scenario. For a comprehensive numerical analysis, beside fitting the neutrino oscillation global fit data, we also present a study on the $\nu_\mu\rightarrow \nu_e$ oscillation which is expected to show up Dirac CP violation in different long baseline experiments. Finally, assuming purely astrophysical sources, we calculate the Ultra High Energy (UHE) neutrino flavor flux ratios at neutrino telescopes, such as IceCube, from which statements on the octant of $\theta_{23}$ could be made in our model.

The Andromeda (M31) galaxy subtends nearly 100 sq. deg. on the sky, with severe contamination from the Milky Way halo stars whose surface density displays a steep gradient across the entire M31 field-of-view. Planetary Nebulae (PNe) are a population of stars firmly associated with M31, that are excellent tracers of light, chemistry and motion in galaxies. We present a 16 sq. deg. survey of the disk and inner halo of M31 with MegaCam@CFHT to identify PNe, characterize their luminosity-specific PN number and luminosity function (PNLF). PNe were identified based on their bright OIII 5007 $\unicode{x212B}$ emission and absence of a continuum. Subsamples of the faint PNe were independently confirmed by matching with resolved Hubble Space Telescope sources from the PHAT survey and spectroscopic follow-up observations with HectoSpec@MMT. The current survey reaches 2 mag fainter than the previous most-sensitive survey. We identify 4289 PNe, of which only 1099 were previously known. By comparing the PN number density with the surface brightness profile of M31 out to ~30 kpc along the minor-axis, we find that the stellar population in the inner halo has a 7 times larger luminosity-specific PN number value than that of the disk. It indicates that the stellar population at deprojected minor-axis radii larger than ~10 kpc is different from that in the M31 disk. We measure the PNLF and find a bright cut-off and a slope consistent with the previous determination by Ciardullo et al. (1989). Interestingly, it shows a significant rise at the faint end, present in all radial bins covered by the survey, much steeper than that observed for the Magellanic clouds and Milky Way bulge. M31 shows two major episodes of star formation and the rise in the faint end of the PNLF is possibly associated with the older stellar population. It may also be a result of varying opacity of the PNe.

Blue Straggler Stars (BSSs) are observed in Galactic globular clusters and old open clusters. The radial distribution of BSSs has been used to diagnose the dynamical evolution of globular clusters. For the first time, with a reliable sample of BSSs identified with Gaia DR2, we conduct such an analysis for an open cluster. We identify members, including BSSs, of the oldest known Galactic open cluster Berkeley 17 with the Gaia DR2 proper motions and parallaxes. We study the radial distribution of the BSS population to understand the dynamical evolution of the cluster. We select cluster members to populate the colour magnitude diagram in the Gaia filters. Cluster parameters are derived using the brightest members. The BSSs and giant branch stars are identified, and their radial distributions are compared. The segregation of BSSs is also evaluated with respect to the giant branch stars using the Minimum Spanning Tree analysis. We determine Berkeley 17 to be at $3138.6^{+285.5}_{-352.9}$ pc. We find 23 BSS cluster members, only two of which were previously identified. We find a bimodal radial distribution of BSSs supported by findings from the MST method. The bimodal radial distribution of BSSs in Berkeley 17 indicates that they have just started to sink towards the cluster center, placing Berkeley 17 with globular clusters of intermediate dynamical age. This is the first such determination for an open cluster.

In models of inflation obtained from string compactification, moduli vacuum misalignment leads to an epoch in the post-inflationary history of the universe when the energy density is dominated by cold moduli particles. This effect leads to a modification in the number of $e$-foldings ($N_{\rm pivot}$) between horizon exit of the CMB modes and the end of inflation. Taking K\"ahler moduli inflation as a prototype, the shift in $e$-foldings turns out to be a function of the model parameters which also determines the inflationary observables. We analyse this scenario numerically using publicly available {\sc{ModeChord}} and {\sc{CosmoMC}} with the latest \emph{Planck+BICEP2/Keck array} data to constrain the model parameters and $N_{\rm pivot}$. In light of the present and future precision data, the results show the importance of careful consideration of any post-inflationary non-standard epoch, as well as of the effects of reheating.

The scalar and vector cosmological perturbations at all length scales of our Universe are studied in the framework of the phantom braneworld model. The model is characterized by the parameter $\Omega_M\equiv M^3/2m^2H_0$, with $M$ and $m$ the 5- and 4-dimensional Planck scales, respectively, and $H_0$ the Hubble parameter today, while $\Omega_M\rightarrow 0$ recovers the $\Lambda\rm CDM$ model. Ignoring the backreaction due to the peculiar velocities and also the bulk cosmological constant, allows the explicit computation of the gravitational potentials, $\Phi$ and $\Psi$. They exhibit exponentially decreasing screening behaviour characterized by a screening length which is a function of the quasidensity parameter $ \Omega_M$.

In this work we study the cosmological attractor models of inflation in connection with certain scalar-tensor theories of gravity, e.g $f(R)$ gravity and Brans-Dicke theory. For some particular choices of the functional degrees of freedom in these theories, one obtains Starobinsky like predictions in the ($n_s$-$r$) observable plane. We have demonstrated that these choices in the Lagrangian density of certain $f(R)$ and Brans-Dicke theories fulfil the condition of the cosmological attractors. That explains why known predictions of $f(R)$ and Brans-Dicke theories in certain cases appear to be the predictions of the much discussed attractor theories. In addition, we did an analysis showing how the predictions of an attractor model is preserved with respect to the variation in the functional freedom of the theory.

According to the third law of Thermodynamics, it takes an infinite number of steps for any object, including black-holes, to reach zero temperature. For any physical system, the process of cooling to absolute zero corresponds to erasing information or generating pure states. In contrast with the ordinary matter, the black-hole temperature can be lowered only by adding matter-energy into it. However, it is impossible to remove the statistical fluctuations of the infalling matter-energy. The fluctuations lead to the fact the black-holes have a finite lower temperature and, hence, an upper bound on the horizon radius. We make an estimate of the upper bound for the horizon radius which is curiosly comparable to Hubble horizon. We compare this bound with known results and discuss its implications.

The formation and properties of star clusters formed at the edges of H II regions are poorly known. We study stellar content, physical conditions, and star formation processes around a relatively unknown young H II region IRAS 10427-6032, located in the southern outskirts of the Carina Nebula. We make use of near-IR data from VISTA, mid-IR from Spitzer and WISE, far-IR from Herschel, sub-mm from ATLASGAL, and 843 MHz radio-continuum data. Using multi-band photometry, we find a total of 5 Class I and 29 Class II young stellar object (YSO) candidates, most of which newly identified, in the 5'$\times$5' region centered on the IRAS source position. Modeling of the spectral energy distribution for selected YSO candidates using radiative transfer models shows that most of these candidates are intermediate mass YSOs in their early evolutionary stages. A majority of the YSO candidates are found to be coincident with the cold dense clump at the western rim of the H II region. Lyman continuum luminosity calculation using radio emission indicates the spectral type of the ionizing source to be earlier than B0.5-B1. We identified a candidate massive star possibly responsible for the H II region with an estimated spectral type B0-B0.5. The temperature and column density maps of the region constructed by performing pixel-wise modified blackbody fits to the thermal dust emission using the far-IR data show a high column density shell-like morphology around the H II region, and low column density (0.6 $\times$ 10$^{22}$ cm$^{-2}$) and high temperature ($\sim$21 K) matter within the H II region. Based on the morphology of the region in the ionized and the molecular gas, and the comparison between the estimated timescales of the H II region and the YSO candidates in the clump, we argue that the enhanced star-formation at the western rim of the H II region is likely due to compression by the ionized gas.

We confront the warm inflation observational predictions directly with the latest CMB data. We focus on a linear temperature (T) dissipative coefficient combined with the simplest model of inflation, a quartic chaotic potential. Although excluded in its standard cold inflation version, dissipation reduces the tensor-to-scalar ratio and brings the quartic chaotic model within the observable allowed range. We will use the CosmoMC package to derive constraints on the model parameters: the combination of coupling constants giving rise to dissipation, the effective number of relativistic degrees of freedom contributing to the thermal bath, and the quartic coupling in the inflaton potential. We do not assume a priori a power-law primordial spectrum, neither we fix the no. of e-folds at the horizon exit. The relation between the no. of e-folds and the comoving scale at horizon crossing is derived from the dynamics, depending on the parameters of the model, which allows us to obtain the $k$-dependent primordial power spectrum. We study the two possibilities considered in the literature for the spectrum, with the inflaton fluctuations having a thermal or a non-thermal origin, and discuss the ability of the data to constraint the model parameters.

Since the Schwarzschild-de Sitter spacetime is static inside the cosmological event horizon, if the dark energy state parameter is sufficiently close to $-1$, apparently one could still expect an effectively static geometry, in the attraction dominated region inside the maximum turn around radius, $R_{\rm TA, max}$, of a cosmic structure. We take the first order metric derived recently assuming a static and ideal dark energy fluid with equation of state $P(r)=\alpha\rho(r)$ as a source in Ref. [1], which reproduced the expression for $R_{\rm TA, max}$ found earlier in the cosmological McVittie spacetime. Here we show that the equality originates from the equivalence of geodesic motion in these two backgrounds, in the non-relativistic regime. We extend this metric up to the third order and compute the bending of light using the Rindler-Ishak method. For $ \alpha\neq -1$, a dark energy dependent term appears in the bending equation, unlike the case of the cosmological constant, $\alpha=-1$. Due to this new term in particular, existing data for the light bending at galactic scales yields, $(1+\alpha)\lesssim {\cal O}(10^{-14})$, thereby practically ruling out any such static and inhomogeneous dark energy fluid we started with. Implication of this result pertaining the uniqueness of the Schwarzschild-de Sitter spacetime in such inhomogeneous dark energy background is discussed.

We present predictions of the K\"ahler moduli inflation model for the spectral tilt by parametrising the reheating epoch by an effective equation-of-state parameter and the number of e-foldings of reheating; and taking into account the post-inflationary history of the model. This model has an epoch in which the energy density of the universe is dominated by cold moduli particles. We compare our results with data from the PLANCK mission and find that exotic reheating (with effective equation of state $w_{\rm re}$ greater than 1/3) is required to match the observations. For canonical reheating case with $w_{\rm re} = 0$, we deduce $\text{log}_{10}(T_{\rm re}/10^3~ \text{GeV}) \simeq 1190 (n_s - 0.956)$. We also analyse our results in the context of observations being planned for the future and their projected sensitivities.

Minimal scalar Higgs portal dark matter model is increasingly in tension with recent results form direct detection experiments like LUX and XENON. In this paper we make a systematic study of minimal extension of the $ \mathbb{Z}_2$ stabilised singlet scalar Higgs portal scenario in terms of their prospects at direct detection experiments. We consider both enlarging the stabilising symmetry to $\mathbb{Z}_3$ and incorporating multipartite features in the dark sector. We demonstrate that in these non-minimal models the interplay of annihilation, co-annihilation and semi-annihilation processes considerably relax constraints from present and proposed direct detection experiments while simultaneously saturating observed dark matter relic density. We explore in particular the resonant semi-annihilation channel within the multipartite $\mathbb{Z}_3$ framework which results in new unexplored regions of parameter space that would be difficult to constrain by direct detection experiments in the near future. The role of dark matter exchange processes within multi-component $\mathbb{Z}_3 \times \mathbb{Z}_3'$ framework is illustrated. We make quantitative estimates to elucidate the role of the various annihilation processes in the different allowed regions of parameter space of these models.

We consider the phantom braneworld model in the context of the maximum turn around radius, $R_{\rm TA,max}$, of a stable, spherical cosmic structure with a given mass. The maximum turn around radius is the point where the attraction due to the central inhomogeneity gets balanced with the repulsion of the ambient dark energy, beyond which a structure cannot hold any mass, thereby giving the maximum upper bound on the size of a stable structure. In this work we derive an analytical expression of $R_{\rm TA,max}$ for this model using cosmological scalar perturbation theory. Using this we numerically constrain the parameter space, including a bulk cosmological constant and the Weyl fluid, from the mass versus observed size data for some nearby, non-virial cosmic structures. We use different values of the matter density parameter $\Omega_m$, both larger and smaller than that of the $\Lambda{\rm CDM}$, as the input in our analysis. We show in particular, that a) with a vanishing bulk cosmological constant the predicted upper bound is always greater than what is actually observed; similar conclusion holds if the bulk cosmological constant is negative b) if it is positive, the predicted maximum size can go considerably below than what is actually observed and owing to the involved nature of the field equations, it leads to interesting constraints on not only the bulk cosmological constant itself but on the whole parameter space of the theory.

We combined Chandra ACIS observations of the globular cluster 47 Tucanae (hereafter, 47 Tuc) from 2000, 2002, and 2014-15 to create a deeper X-ray source list, and study some of the faint radio millisecond pulsars (MSPs) present in this cluster. We have detected 370 X-ray sources within the half-mass radius (2$'$.79) of the cluster, 81 of which are newly identified, by including new data and using improved source detection techniques. The majority of the newly identified sources are in the crowded core region, indicating cluster membership. We associate five of the new X-ray sources with chromospherically active BY Dra or W UMa variables identified by Albrow et al. (2001). We present alternative positions derived from two methods, centroiding and image reconstruction, for faint, crowded sources. We are able to extract X-ray spectra of the recently discovered MSPs 47 Tuc aa, 47 Tuc ab, the newly timed MSP 47 Tuc Z, and the newly resolved MSPs 47 Tuc S and 47 Tuc F. Generally, they are well fit by black body or neutron star atmosphere models, with temperatures, luminosities and emitting radii similar to those of other known MSPs in 47 Tuc, though 47 Tuc aa and 47 Tuc ab reach lower X-ray luminosities. We limit X-ray emission from the full surface of the rapidly spinning (542 Hz) MSP 47 Tuc aa, and use this limit to put an upper bound for amplitude of r-mode oscillations in this pulsar as $\alpha<2.5\times 10^{-9}$ and constrain the shape of the r-mode instability window.

We present the analysis of the morphological shape of Berkeley 17, the oldest known open cluster (~10 Gyr), using a probabilistic star counting of Pan-STARRS point sources, and confirm its core-tail shape, plus an antitail, previously detected with the 2MASS data. The stellar population, as diagnosed by the color-magnitude diagram and theoretical isochrones, shows many massive members in the cluster core, whereas there is a paucity of such members in both tails. This manifests mass segregation in this aged star cluster with the low-mass members being stripped away from the system. It has been claimed that Berkeley 17 is associated with an excessive number of blue straggler candidates. Comparison of nearby reference fields indicates that about half of these may be field contamination.

The maximum allowable size of a spherical cosmic structure as a function of its mass is determined by the maximum turn around radius $R_{\rm TA,max}$, the distance from its centre where the attraction on a radial test particle due to the spherical mass is balanced with the repulsion due to the ambient dark energy. In this work, we extend the existing results in several directions. (a) We first show that for $w\neq -1$, the expression for $R_{\rm TA, max}$ found earlier using the cosmological perturbation theory, can be derived using a static geometry as well. (b) In the generic dark energy model with arbitrary time dependent state parameter $w(t)$, taking into account the effect of inhomogeneities upon the dark energy as well, where it is shown that the data constrain $w(t={\rm today})>-2.3$, and (c) in the quintessence and the generalized Chaplygin gas models, both of which are shown to predict structure sizes consistent with observations.

The maximum size of a cosmic structure is given by the maximum turnaround radius -- the scale where the attraction due to its mass is balanced by the repulsion due to dark energy. We derive generic formulae for the estimation of the maximum turnaround radius in any theory of gravity obeying the Einstein equivalence principle, in two situations: on a spherically symmetric spacetime and on a perturbed Friedman-Robertson-Walker spacetime. We show that the two formulae agree. As an application of our formula, we calculate the maximum turnaround radius in the case of the Brans-Dicke theory of gravity. We find that for this theory, such maximum sizes always lie above the \LCDM value, by a factor $1 + \frac{1}{3\omega}$, where $\omega\gg 1$ is the Brans-Dicke parameter, implying consistency of the theory with current data.

Black-holes in asymptotically flat space-times have negative specific heat --- they get hotter as they loose energy. A clear statistical mechanical understanding of this has remained a challenge. In this work, we address this issue using fluid-gravity correspondence which aims to associate fluid degrees of freedom to the horizon. Using linear response theory and the teleological nature of event horizon, we show explicitly that the fluctuations of the horizon-fluid lead to negative specific heat for Schwarzschild black Hole. We also point out how the specific heat can be positive for Kerr-Newman or AdS black holes. Our approach constitutes an important advance as it allows us to apply canonical ensemble approach to study thermodynamics of asymptotically flat black-hole space-times.

In the background of a homogeneous and isotropic spacetime with zero spatial curvature, we consider interacting scenarios between two barotropic fluids, one is the pressureless dark matter (DM) and the other one is dark energy (DE), in which the equation of state (EoS) in DE is either constant or time dependent. In particular, for constant EoS in DE, we show that the evolution equations for both fluids can be analytically solved. For all these scenarios, the model parameters have been constrained using the current astronomical observations from Type Ia Supernovae, Hubble parameter measurements, and baryon acoustic oscillations distance measurements. Our analysis shows that both for constant and variable EoS in DE, a very small but nonzero interaction in the dark sector is favored while the EoS in DE can predict a slight phantom nature, i.e. the EoS in DE can cross the phantom divide line `$-1$'. On the other hand, although the models with variable EoS describe the observations better, but the Akaike Information Criterion supports models with minimal number of parameters. However, it is found that all the models are very close to the $\Lambda$CDM cosmology.

The morphology and cluster membership of the Galactic open clusters - Czernik 20 and NGC 1857 were analyzed using two different clustering algorithms. We present the maiden use of density-based spatial clustering of applications with noise (DBSCAN) to determine open cluster morphology from spatial distribution. The region of analysis has also been spatially classified using a statistical membership determination algorithm. We utilized near infrared (NIR) data for a suitably large region around the clusters from the United Kingdom Infrared Deep Sky Survey Galactic Plane Survey star catalogue database, and also from the Two Micron All Sky Survey star catalogue database. The densest regions of the cluster morphologies (1 for Czernik 20 and 2 for NGC 1857) thus identified were analyzed with a K-band extinction map and color-magnitude diagrams (CMDs). To address significant discrepancy in known distance and reddening parameters, we carried out field decontamination of these CMDs and subsequent isochrone fitting of the cleaned CMDs to obtain reliable distance and reddening parameters for the clusters (Czernik 20: D = 2900 pc; E(J-K) = 0.33; NGC 1857: D = 2400 pc; E(J-K) = 0.18-0.19). The isochrones were also used to convert the luminosity functions for the densest regions of Czernik 20 and NGC 1857 into mass function, to derive their slopes. Additionally, a previously unknown over-density consistent with that of a star cluster is identified in the region of analysis.

We investigate the stability of stationary integral solutions of an ideal irrotational fluid in a general static and spherically symmetric background, by studying the profile of the perturbation of the mass accretion rate. We consider low angular momentum axisymmetric accretion flows for three different accretion disk models and consider time dependent and radial linear perturbation of the mass accretion rate. First we show that the propagation of such perturbation can be determined by an effective $2\times2$ matrix, which has qualitatively similar acoustic causal properties as one obtains via the perturbation of the velocity potential. Next, using this matrix we analytically address the stability issues, for both standing and travelling wave configurations generated by the perturbation. Finally, based on this general formalism we briefly discuss the explicit example of the Schwarzschild spacetime and compare our results of stability with the existing literature, which instead address this problem via the perturbation of the velocity potential.

The maximum size of a bound cosmic structure is computed perturbatively as a function of its mass in the framework of the cubic galileon, proposed recently to model the dark energy of our Universe. Comparison of our results with observations constrains the matter-galileon coupling of the model to $0.033\lesssim \alpha \lesssim 0.17$, thus improving previous bounds based solely on solar system physics.

Ground and space-based sky surveys enable powerful cosmological probes based on measurements of galaxy properties and the distribution of galaxies in the Universe. These probes include weak lensing, baryon acoustic oscillations, abundance of galaxy clusters, and redshift space distortions; they are essential to improving our knowledge of the nature of dark energy. On the theory and modeling front, large-scale simulations of cosmic structure formation play an important role in interpreting the observations and in the challenging task of extracting cosmological physics at the needed precision. These simulations must cover a parameter range beyond the standard six cosmological parameters and need to be run at high mass and force resolution. One key simulation-based task is the generation of accurate theoretical predictions for observables, via the method of emulation. Using a new sampling technique, we explore an 8-dimensional parameter space including massive neutrinos and a variable dark energy equation of state. We construct trial emulators using two surrogate models (the linear power spectrum and an approximate halo mass function). The new sampling method allows us to build precision emulators from just 26 cosmological models and to increase the emulator accuracy by adding new sets of simulations in a prescribed way. This allows emulator fidelity to be systematically improved as new observational data becomes available and higher accuracy is required. Finally, using one LCDM cosmology as an example, we study the demands imposed on a simulation campaign to achieve the required statistics and accuracy when building emulators for dark energy investigations.

In this work we present an alternative derivation of the general relativistic acoustic analogue geometry by perturbing the mass accretion rate or flux of an ideal fluid flowing radially in a general static and spherically symmetric spacetime. To the best of our knowledge, this has so far been done in non-relativistic scenario. The resulting causal structure of the two dimensional acoustic geometry is qualitatively similar to that one derives via the perturbation of the velocity potential. Using this, we then briefly discuss the stability issues by studying the wave configurations generated by the perturbation of the mass accretion rate, and formally demonstrate the stability of the accretion process. This is in qualitative agreement with earlier results on stability, established via study of wave configurations generated by the perturbation of velocity potential, by using the acoustic geometry associated with it. We further discuss explicit examples of the Schwarzschild and Rindler spacetimes.

We introduce a novel formalism to investigate the role of the spin angular momentum of astrophysical black holes in influencing the behaviour of low angular momentum general relativistic accretion. We propose a metric independent analysis of axisymmetric general relativistic flow, and consequently formulate the space and time dependent equations describing the general relativistic hydrodynamic accretion flow in the Kerr metric. The associated stationary critical solutions for such flow equations are provided and the stability of the stationary transonic configuration is examined using an elegant linear perturbation technique. We examine the properties of infalling material for both prograde and retrograde accretion as a function of the Kerr parameter at extremely close proximity to the event horizon. Our formalism can be used to identify a new spectral signature of black hole spin, and has the potential of performing the black hole shadow imaging corresponding to the low angular momentum accretion flow.

Modeling large-scale sky survey observations is a key driver for the continuing development of high resolution, large-volume, cosmological simulations. We report the first results from the 'Q Continuum' cosmological N-body simulation run carried out on the GPU-accelerated supercomputer Titan. The simulation encompasses a volume of (1300 Mpc)^3 and evolves more than half a trillion particles, leading to a particle mass resolution of ~1.5 X 10^8 M_sun. At this mass resolution, the Q Continuum run is currently the largest cosmology simulation available. It enables the construction of detailed synthetic sky catalogs, encompassing different modeling methodologies, including semi-analytic modeling and sub-halo abundance matching in a large, cosmological volume. Here we describe the simulation and outputs in detail and present first results for a range of cosmological statistics, such as mass power spectra, halo mass functions, and halo mass-concentration relations for different epochs. We also provide details on challenges connected to running a simulation on almost 90% of Titan, one of the fastest supercomputers in the world, including our usage of Titan's GPU accelerators.

We propose a novel linear perturbation scheme to study the stability properties of the stationary transonic integral solutions for axisymmetric matter flow around astrophysical black holes for the Schwarzschild as well as for rotating Rindler space time. We discuss the emergence of the relativistic acoustic geometry as a consequence of such stability analysis. Our work thus makes a crucial connection between two apparently non-overlapping fields of research - the accretion astrophysics and the analogue gravity phenomena.

We propose a method to estimate the location of the Sun in the disk of the Milky Way using a method based on the Hellinger distance and construct confidence sets on our estimate of the unknown location using a bootstrap based method. Assuming the Galactic disk to be two-dimensional, the sought solar location then reduces to the radial distance separating the Sun from the Galactic center and the angular separation of the Galactic center to Sun line, from a pre-fixed line on the disk. On astronomical scales, the unknown solar location is equivalent to the location of us earthlings who observe the velocities of a sample of stars in the neighborhood of the Sun. This unknown location is estimated by undertaking pairwise comparisons of the estimated density of the observed set of velocities of the sampled stars, with densities estimated using synthetic stellar velocity data sets generated at chosen locations in the Milky Way disk according to four base astrophysical models. The "match" between the pair of estimated densities is parameterized by the affinity measure based on the familiar Hellinger distance. We perform a novel cross-validation procedure to establish a desirable "consistency" property of the proposed method.

We present a measurement of the angular bispectrum of the millimeter-wave sky in observing bands centered at roughly 95, 150, and 220 GHz, on angular scales of $1^\prime \lesssim \theta \lesssim 10^\prime$ (multipole number $1000 \lesssim l \lesssim 10000$). At these frequencies and angular scales, the main contributions to the bispectrum are expected to be the thermal Sunyaev-Zel'dovich (tSZ) effect and emission from extragalactic sources, predominantly dusty, star-forming galaxies (DSFGs) and active galactic nuclei. We measure the bispectrum in 800 $\mathrm{deg}^2$ of three-band South Pole Telescope data, and we use a multi-frequency fitting procedure to separate the bispectrum of the tSZ effect from the extragalactic source contribution. We simultaneously detect the bispectrum of the tSZ effect at $>$10$\sigma$, the unclustered component of the extragalactic source bispectrum at $>$5$\sigma$ in each frequency band, and the bispectrum due to the clustering of DSFGs---i.e., the clustered cosmic infrared background (CIB) bispectrum---at $>$5$\sigma$. This is the first reported detection of the clustered CIB bispectrum. We use the measured tSZ bispectrum amplitude, compared to model predictions, to constrain the normalization of the matter power spectrum to be $\sigma_8 = 0.787 \pm 0.031$ and to predict the amplitude of the tSZ power spectrum at $l = 3000$. This prediction improves our ability to separate the thermal and kinematic contributions to the total SZ power spectrum. The addition of bispectrum data improves our constraint on the tSZ power spectrum amplitude by a factor of two compared to power spectrum measurements alone and demonstrates a preference for a nonzero kinematic SZ (kSZ) power spectrum, with a derived constraint on the kSZ amplitude at $l=3000$ of A_kSZ $ = 2.9 \pm 1.6 \ \mu$K$^2$, or A_kSZ $ = 2.6 \pm 1.8 \ \mu$K$^2$ if the default A_kSZ > 0 prior is removed.

For Rindler observers accelerating close to the horizon in local patches around a spacetime point, the matter-energy passing through the horizon increases the entropy and heat energy. Jacobson has showed that the Einstein equation can be derived from the consideration of this thermodynamic process. This, however, works only if the acceleration $a$ is much larger than the scale set by the curvature of the spacetime. It is explored here whether an extension is possible to the case with no lower bound on $a$. We show that this is possible if one assumes that in a locally accelerating frame, the matter-energy passing through null hypersurfaces could result in an increase in the heat energy and the entropy. Such a generalisation extends the thermodynamic derivation of gravity to include any non-freely falling observer. A new method of determining the temperature detected by such locally accelerating observers is also presented. By considering only the quantisation of sufficiently localised wave modes of a field, it is shown that the observer finds himself in a thermal environment.

The lensing power spectrum from cosmic microwave background (CMB) temperature maps will be measured with unprecedented precision with upcoming experiments, including upgrades to ACT and SPT. Achieving significant improvements in cosmological parameter constraints, such as percent level errors on sigma_8 and an uncertainty on the total neutrino mass of approximately 50 meV, requires percent level measurements of the CMB lensing power. This necessitates tight control of systematic biases. We study several types of biases to the temperature-based lensing reconstruction signal from foreground sources such as radio and infrared galaxies and the thermal Sunyaev-Zel'dovich effect from galaxy clusters. These foregrounds bias the CMB lensing signal due to their non-Gaussian nature. Using simulations as well as some analytical models we find that these sources can substantially impact the measured signal if left untreated. However, these biases can be brought to the percent level if one masks galaxies with fluxes at 150 GHz above 1 mJy and galaxy clusters with masses above M_vir = 10^14 M_sun. To achieve such percent level bias, we find that only modes up to a maximum multipole of l_max ~ 2500 should be included in the lensing reconstruction. We also discuss ways to minimize additional bias induced by such aggressive foreground masking by, for example, exploring a two-step masking and in-painting algorithm.

We measure the cross-power spectrum of the projected mass density as traced by the convergence of the cosmic microwave background lensing field from the South Pole Telescope (SPT) and a sample of Type 1 and 2 (unobscured and obscured) quasars at z~1 selected with the Wide-field Infrared Survey Explorer, over 2500 deg^2. The cross-power spectrum is detected at ~7-sigma, and we measure a linear bias b=1.67+/-0.24, consistent with clustering analyses. Using an independent lensing map, derived from Planck observations, to measure the cross-spectrum, we find excellent agreement with the SPT analysis. The bias of the combined sample of Type 1 and 2 quasars determined in this work is similar to that previously determined for Type 1 quasars alone; we conclude that that obscured and unobscured quasars must trace the matter field in a similar way. This result has implications for our understanding of quasar unification and evolution schemes.

We use a temperature map of the cosmic microwave background (CMB) obtained using the South Pole Telescope at 150 GHz to construct a map of the gravitational convergence to z ~ 1100, revealing the fluctuations in the projected mass density. This map shows individual features that are significant at the ~ 4 sigma level, providing the first image of CMB lensing convergence. We cross-correlate this map with Herschel/SPIRE maps covering 90 square degrees at wavelengths of 500, 350, and 250 microns. We show that these submillimeter-wavelength (submm) maps are strongly correlated with the lensing convergence map, with detection significances in each of the three submm bands ranging from 6.7 to 8.8 sigma. We fit the measurement of the cross power spectrum assuming a simple constant bias model and infer bias factors of b=1.3-1.8, with a statistical uncertainty of 15%, depending on the assumed model for the redshift distribution of the dusty galaxies that are contributing to the Herschel/SPIRE maps.

We perform a detailed study of the bispectrum of the Sunyaev-Zel'dovich effect. Using an analytical model for the pressure profiles of the intracluster medium, we demonstrate the SZ bispectrum to be a sensitive probe of the amplitude of the matter power spectrum parameter sigma_8. We find that the bispectrum amplitude scales as B_SZ ~ sigma_8^{11-12}, compared to that of the power spectrum, which scales as A_tSZ ~ sigma_8^{7-9}. We show that the SZ bispectrum is principally sourced by massive clusters at redshifts around z~0.4, which have been well-studied observationally. This is in contrast to the SZ power spectrum, which receives a significant contribution from less-well understood low-mass and high-redshift groups and clusters. Therefore, the amplitude of the bispectrum at l~3000 is less sensitive to astrophysical uncertainties than the SZ power spectrum. We show that current high resolution CMB experiments should be able to detect the SZ bispectrum amplitude with high significance, in part due to the low contamination from extra-galactic foregrounds. A combination of the SZ bispectrum and the power spectrum can sharpen the measurements of thermal and kinetic SZ components and help distinguish cosmological and astrophysical information from high-resolution CMB maps.

The concentration-mass relation for dark matter-dominated halos is one of the essential results expected from a theory of structure formation. We present a simple prediction scheme, a cosmic emulator, for the c-M relation as a function of cosmological parameters for wCDM models. The emulator is constructed from 37 individual models, with three nested N-body gravity-only simulations carried out for each model. The mass range covered by the emulator is 2 x 10^{12} M_sun < M <10^{15} M_sun with a corresponding redshift range of z=0 -1. Over this range of mass and redshift, as well as the variation of cosmological parameters studied, the mean halo concentration varies from c ~ 2 to c ~ 8. The distribution of the concentration at fixed mass is Gaussian with a standard deviation of one-third of the mean value, almost independent of cosmology, mass, and redshift over the ranges probed by the simulations. We compare results from the emulator with previously derived heuristic analytic fits for the c-M relation, finding that they underestimate the halo concentration at high masses. Using the emulator to investigate the cosmology dependence of the c-M relation over the currently allowable range of values, we find -- not surprisingly -- that \sigma_8 and \omega_m influence it considerably, but also that the dark energy equation of state parameter, w, has a substantial effect. In general, the concentration of lower-mass halos is more sensitive to changes in cosmological parameters as compared to cluster mass halos.

Dark matter-dominated cluster-scale halos act as an important cosmological probe and provide a key testing ground for structure formation theory. Focusing on their mass profiles, we have carried out (gravity-only) simulations of the concordance LCDM cosmology, covering a mass range of 2.10^{12}-2.10^{15} solar mass/h and a redshift range of z=0-2, while satisfying the associated requirements of resolution and statistical control. When fitting to the Navarro-Frenk-White profile, our concentration-mass (c-M) relation differs in normalization and shape in comparison to previous studies that have limited statistics in the upper end of the mass range. We show that the flattening of the c-M relation with redshift is naturally expressed if c is viewed as a function of the peak height parameter, \nu. Unlike the c-M relation, the slope of the c-\nu relation is effectively constant over the redshift range z=0-2, while the amplitude varies by ~30% for massive clusters. This relation is, however, not universal: Using a simulation suite covering the allowed wCDM parameter space, we show that the c-\nu relation varies by about +/- 20% as cosmological parameters are varied. At fixed mass, the c(M) distribution is well-fit by a Gaussian with \sigma_c/c = 0.33, independent of the radius at which the concentration is defined, the halo dynamical state, and the underlying cosmology. We compare the LCDM predictions with observations of halo concentrations from strong lensing, weak lensing, galaxy kinematics, and X-ray data, finding good agreement for massive clusters (M > 4.10^{14} solar mass/h), but with some disagreements at lower masses. Because of uncertainty in observational systematics and modeling of baryonic physics, the significance of these discrepancies remains unclear.

An excited state of $^{12}$C having excitation energy E$_x \sim$ 9.65 $\pm$ 0.02 MeV and width (FWHM) $\sim607\pm$ 55 keV, which decays to three $ \alpha $-particles via Hoyle state ($E_x \sim$ 7.65 MeV), has been directly identified for the first time in the exclusive inelastic scattering of 60 MeV $^{4}$He on $^{12}$C, measured in coincidence with the recoiling $^{12}$C$ ^* $ Hoyle state (decaying mostly as $^{12}$C$ ^* $ $\rightarrow \ ^{8} $Be + $ \alpha $ $\rightarrow \ \alpha + \alpha + \alpha$) by event-by-event kinematic reconstruction of the completely detected (4$ \alpha $) events. This state is likely to be a candidate for 2$_2^+$ first excited Hoyle state, the existence of which has recently been indirectly evidenced from the recent inclusive inelastic scattering studies.

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$.

In this work, we have studied the thermodynamic quantities like temperature of the universe, heat capacity and squared speed of sound in generalized gravity theories like Brans-Dicke, Ho$\check{\text r}$ava-Lifshitz and $f(R)$ gravities. We have considered the universe filled with dark matter and dark energy. Also we have considered the equation of state parameters for open, closed and flat models. We have observed that in all cases the equation of state behaves like quintessence. The temperature and heat capacity of the universe are found to decrease with the expansion of the universe in all cases. In Brans-Dicke and $f(R)$ gravity theories the squared speed of sound is found to exhibit increasing behavior for open, closed and flat models and in Ho$\check{\text r}$ava-Lifshitz gravity theory it is found to exhibit decreasing behavior for open and closed models with the evolution of the universe. However, for flat universe, the squared speed of sound remains constant in Ho$\check{\text r}$ava-Lifshitz gravity.

We study the collapse of a massless scalar field coupled to gravity. A class of blackhole solutions are identified. We also report on a class of solutions where collapse starts from a regular spacelike surface but then the collapsing scalar field freezes. As a result, in these solutions, a black hole does not form, neither is there any singularity in the future.

The mass distribution of halos, as specified by the halo mass function, is a key input for several cosmological probes. The sizes of $N$-body simulations are now such that, for the most part, results need no longer be statistics-limited, but are still subject to various systematic uncertainties. We investigate and discuss some of the reasons for these differences. Quantifying error sources and compensating for them as appropriate, we carry out a high-statistics study of dark matter halos from 67 $N$-body simulations to investigate the mass function and its evolution for a reference $\Lambda$CDM cosmology and for a set of $w$CDM cosmologies. For the reference $\Lambda$CDM cosmology (close to WMAP5), we quantify the breaking of universality in the form of the mass function as a function of redshift, finding an evolution of as much as 10% away from the universal form between redshifts $z=0$ and $z=2$. For cosmologies very close to this reference we provide a fitting formula to our results for the (evolving) $\Lambda$CDM mass function over a mass range of $6\cdot 10^{11}-3\cdot 10^{15}$ M$_{\odot}$ to an estimated accuracy of about 2%. The set of $w$CDM cosmologies is taken from the Coyote Universe simulation suite. The mass functions from this suite (which includes a $\Lambda$CDM cosmology and others with $w\simeq-1$) are described by the fitting formula for the reference $\Lambda$CDM case at an accuracy level of 10%, but with clear systematic deviations. We argue that, as a consequence, fitting formulae based on a universal form for the mass function may have limited utility in high precision cosmological applications.

Upcoming imaging surveys such as the Large Synoptic Survey Telescope will repeatedly scan large areas of sky and have the potential to yield million-supernova catalogs. Type Ia supernovae are excellent standard candles and will provide distance measures that suffice to detect mean pairwise velocities of their host galaxies. We show that when combining these distance measures with photometric redshifts for either the supernovae or their host galaxies, the mean pairwise velocities of the host galaxies will provide a dark energy probe which is competitive with other widely discussed methods. Adding information from this test to type Ia supernova photometric luminosity distances from the same experiment, plus the cosmic microwave background power spectrum from the Planck satellite, improves the Dark Energy Task Force Figure of Merit by a factor of 1.8. Pairwise velocity measurements require no additional observational effort beyond that required to perform the traditional supernova luminosity distance test, but may provide complementary constraints on dark energy parameters and the nature of gravity. Incorporating additional spectroscopic redshift follow-up observations could provide important dark energy constraints from pairwise velocities alone. Mean pairwise velocities are much less sensitive to systematic redshift errors than the luminosity distance test or weak lensing techniques, and also are only mildly affected by systematic evolution of supernova luminosity.

We study a class of early dark energy (EDE) models, in which, unlike in standard dark energy models, a substantial amount of dark energy exists in the matter-dominated era. We self-consistently include dark energy perturbations, and show that these models may be successfully constrained using future observations of galaxy clusters, in particular the redshift abundance, and the Sunyaev-Zel'dovich (SZ) power spectrum. We make predictions for EDE models, as well as LCDM for incoming X-ray (eROSITA) and microwave (South Pole Telescope) observations. We show that galaxy clusters' mass function and the SZ power spectrum will put strong constraints both on the equation of state of \de today and the redshift at which EDE transits to present-day LCDM like behavior for these models, thus providing complementary information to the geometric probes of dark energy. Not including perturbations in EDE models leads to those models being practically indistinguishable from LCDM. An MCMC analysis of future galaxy cluster surveys provides constraints for EDE parameters that are competitive with and complementary to background expansion observations such as supernovae.

We use an analytic model to investigate the theoretical uncertainty on the thermal Sunyaev-Zel'dovich (SZ) power spectrum due to astrophysical uncertainties in the thermal structure of the intracluster medium. Our model accounts for star formation and energy feedback (from supernovae and active galactic nuclei) as well as radially dependent non-thermal pressure support due to random gas motions, the latter calibrated by recent hydrodynamical simulations. We compare the model against X-ray observations of low redshift clusters, finding excellent agreement with observed pressure profiles. Varying the levels of feedback and non-thermal pressure support can significantly change both the amplitude and shape of the thermal SZ power spectrum. Increasing the feedback suppresses power at small angular scales, shifting the peak of the power spectrum to lower ell. On the other hand, increasing the non-thermal pressure support has the opposite effect, significantly reducing power at large angular scales. In general, including non-thermal pressure at the level measured in simulations has a large effect on the power spectrum, reducing the amplitude by 50% at angular scales of a few arcminutes compared to a model without a non-thermal component. Our results demonstrate that measurements of the shape of the power spectrum can reveal useful information on important physical processes in groups and clusters, especially at high-redshift where there exists little observational data. Comparing with the recent South Pole Telescope measurements of the small-scale cosmic microwave background power spectrum, we find our model reduces the tension between the values of sigma_8 measured from the SZ power spectrum and from cluster abundances.

Cosmological constraints from cluster surveys rely on accurate mass estimates from the mass-observable relations. In order to avoid systematic biases and reduce uncertainties, we study the form and physical origin of the intrinsic scatter about the mean Sunyaev-Zel'dovich (SZ) flux-mass relation using a hydrodynamical simulation of galaxy cluster formation. We examine the assumption of lognormal scatter and detect non-negligible positive skewness and kurtosis (> 0.5) for a wide range of limiting masses and redshifts. These higher-order moments should be included in the parametrization of scatter in order not to bias cosmological constraints. We investigate the sources of the scatter by correlating it with measures of cluster morphology, halo concentration, and dynamical state, and we quantify the individual contribution from each source. We find that statistically the impact of dynamical state is weak, so the selection bias due to mergers is negligible. On the other hand, there is a strong correlation between the scatter and halo concentration, which can be used to reduce the scatter significantly (from 12.07% to 7.34% or by ~40% for clusters at z = 0). We also show that a cross-calibration by combining information from X-ray followups can be used to reduce the scatter in the flux-mass relation and also identify outliers in both X-ray and SZ cluster surveys.

In this work, we have described the Brans-Dicke theory of gravity and given a particular solution by choosing a power law form of scalar field $\phi$ and constant $\omega$. If we assume first law and entropy formula on apparent horizon then we recover Friedmann equations. Next, assuming first law of thermodynamics, the validity conditions of GSL on event horizon are presented. Also without use first law, if we impose the entropy relation on the horizon, then we also obtain the condition of validity of GSL on event horizon. The validity of GSL completely depends on the model of BD scalar field solutions. We have justified that on the apparent horizon the two process are equivalent, but on the event horizon they are not equivalent. If first law is valid on the event horizon then GSL may be satisfied in BD solution, but if first law is not satisfied then GSL is not satisfied in BD solution. So first law always favours GSL on event horizon. In our effective approach, the first law and GSL is always satisfied in apparent horizon, which do not depend on BD theory of gravity.

In this work, we have investigated the validity of GSL of thermodynamics in a universe (open, closed and flat) governed by Ho$\check{\text r}$ava-Lifshitz gravity. If the universe contains barotropic fluid the corresponding solutions have been obtained. The validity of GSL have been examined by two approaches: (i) robust approach and (ii) effective approach. In robust approach, we have considered the universe contains only matter fluid and the effect of the gravitational sector of HL gravity was incorporated through the modified black hole entropy on the horizon. Effective approach is that all extra information of HL gravity into an effective dark energy fluid and so we consider the universe contains matter fluid plus this effective fluid. This approach is essentially same as the Einstein's gravity theory. The general prescription for validity of GSL have been discussed. Graphically we have shown that the GSL may be satisfied for open, closed and flat universe on the different horizons with different conditions.

We perform a multilepton channel analysis in the context of the Large Hadron Collider (LHC) for Wilkinson Microwave Anisotropy Probe (WMAP) compatible points in a model with non-universal scalar masses, which admits a Higgs funnel region of supersymmetry dark matter even for a small $\tan\beta$. In addition to two and three-lepton final states, four-lepton events, too, are shown to be useful for this purpose. We also compare the collider signatures in similar channels for WMAP compatible points in the minimal supergravity (mSUGRA) framework with similar gluino masses. Some definite features of such non-universal scenario emerge from the analysis.

We present redshifts and optical richness properties of 21 galaxy clusters uniformly selected by their Sunyaev-Zel'dovich signature. These clusters, plus an additional, unconfirmed candidate, were detected in a 178 square-degree area surveyed by the South Pole Telescope in 2008. Using griz imaging from the Blanco Cosmology Survey and from pointed Magellan telescope observations, as well as spectroscopy using Magellan facilities, we confirm the existence of clustered red-sequence galaxies, report red-sequence photometric redshifts, present spectroscopic redshifts for a subsample, and derive R_200 radii and M_200 masses from optical richness. The clusters span redshifts from 0.15 to greater than 1, with a median redshift of 0.74; three clusters are estimated to be at z > 1. Redshifts inferred from mean red-sequence colors exhibit 2% RMS scatter in sigma_z/(1+z) with respect to the spectroscopic subsample for z < 1. We show that M_200 cluster masses derived from optical richness correlate with masses derived from South Pole Telescope data and agree with previously derived scaling relations to within the uncertainties. Optical and infrared imaging is an efficient means of cluster identification and redshift estimation in large Sunyaev-Zel'dovich surveys, and exploiting the same data for richness measurements, as we have done, will be useful for constraining cluster masses and radii for large samples in cosmological analysis.

The accelerating expansion of the Universe at recent epochs has called into question the validity of general relativity on cosmological scales. One probe of gravity is a comparison of expansion history of the Universe with the history of structure growth via gravitational instability: general relativity predicts a specific relation between these two observables. Here we show that the mean pairwise streaming velocity of galaxy clusters provides a useful method of constraining this relation. Galaxy cluster velocities can be measured via the kinetic Sunyaev-Zeldovich distortion of the cosmic microwave background radiation; future surveys can provide large enough catalogs of cluster velocities to discriminate between general relativity and other proposed gravitational theories.

Recent technological advances have led to a flood of new data on cosmology rich in information about the formation and evolution of the universe, e.g., the data collected in Sloan Digital Sky Survey (SDSS) for more than 200 million objects. The analyses of such data demand cutting edge statistical technologies. Here, we have used the concept of mixture model within Bayesian semiparametric methodology to fit the regression curve with the bivariate data for the apparent magnitude and redshift for Quasars in SDSS (2007) catalogue. Associated with the mixture modeling is a highly efficient curve-fitting procedure, which is central to the application considered in this paper. Moreover, we adopt a new method for analysing the posterior distribution of clusterings, also generated as a by-product of our methodology. The results of our analysis of the cosmological data clearly indicate the existence of four change points on the regression curve andposssibiltiy of clustering of quasars specially at high redshift. This sheds new light not only on the issue of evolution, existence of acceleration or decceleration and environment around quasars at high redshift but also help us to estimate thecosmological parameters related to acceleration or decceleration.

We study the utility of a large sample of type Ia supernovae that might be observed in an imaging survey that rapidly scans a large fraction of the sky for constraining dark energy. We consider information from the traditional luminosity distance test as well as the spread in SNeIa fluxes at fixed redshift induced by gravitational lensing. We include a treatment of photometric redshift uncertainties in our analysis. Our primary result is that the information contained in the mean distance moduli of SNeIa and the dispersion among SNeIa distance moduli complement each other, breaking a degeneracy between the present dark energy equation of state and its time variation without the need for a high-redshift supernova sample. To address photometric redshift uncertainties, we present dark energy constraints as a function of the size of an external set of spectroscopically-observed SNeIa that may be used for redshift calibration, nspec. We find that an imaging survey can constrain the dark energy equation of state at the epoch where it is best constrained with a 1-sigma error of sigma(wpiv)~0.03-0.09$, depending upon various assumptions. In addition, the marginal improvement in the error sigma(wpiv) from an increase in the spectroscopic calibration sample drops once nspec ~ 10^3. This result is important because it is of the order of the size of calibration samples likely to be compiled in the coming decade and because, for samples of this size, the spectroscopic and imaging surveys individually place comparable constraints on the dark energy equation of state. In all cases, it is best to calibrate photometric redshifts with a set of spectroscopically-observed SNeIa with relatively more objects at high redshift than the parent sample of imaging SNeIa.

Galaxy cluster surveys compiled via the Sunyaev-Zeldovich Effect have the potential to place strong constraints on cosmology, and in particular the nature of dark energy. Here we consider cluster velocity surveys using kinetic Sunyaev-Zeldovich measurements. Cluster velocities closely trace the large-scale velocity field independent of cluster mass; we demonstrate that two useful cluster velocity statistics are nearly independent of cluster mass, in marked contrast to cluster number count statistics. On the other hand, cluster velocity determinations from three-band observations of Sunyaev-Zeldovich distortions can require additional cluster data or assumptions, and are complicated by microwave emission from dusty galaxies and radio sources, which may be correlated with clusters. Systematic errors in velocity due to these factors can give substantial biases in determination of dark energy parameters, although large cluster velocity surveys will contain enough information that the errors can be modeled using the data itself, with little degradation in cosmological constraints.

We study the effect of quasar feedback on distributions of baryons and properties of intracluster medium in galaxy groups using high-resolution numerical simulations. We use the entropy-conserving Gadget code that includes gas cooling and star formation, modified to include a physically-based model of quasar feedback. For a sample of ten galaxy group-sized dark matter halos with masses in the range of 1 to $5\times 10^{13} M_{\odot}/h$, star formation is suppressed by more than 50% in the inner regions due to the additional pressure support by quasar feedback, while gas is driven from the inner region towards the outer region of the halos. As a result, the average gas density is 50% lower in the inner region and 10% higher in the outer region in the simulation, compared to a similar simulation with no quasar feedback. Gas pressure is lowered by about 40% in the inner region and higher in the outer region, while temperature and entropy are enhanced in the inner region by about 20-40%. The total group gas fraction in the two simulations generally differs by less than 10%. We also find a small change of the total thermal Sunyaev-Zeldovich distortion, leading to 10% changes in the microwave angular power spectrum at angular scales below two arcminutes.

Future multifrequency microwave background experiments with arcminute resolution and micro-Kelvin temperature sensitivity will be able to detect the kinetic Sunyaev-Zeldovich (kSZ) effect, providing a way to measure radial peculiar velocities of massive galaxy clusters. We show that cluster peculiar velocities have the potential to constrain several dark energy parameters. We compare three velocity statistics (the distribution of radial velocities, the mean pairwise streaming velocity, and the velocity correlation function) and analyze the relative merits of these statistics in constraining dark energy parameters. Of the three statistics, mean pairwise streaming velocity provides constraints that are least sensitive to velocity errors: the constraints on parameters degrades only by a factor of two when the random error is increased from 100 to 500 km/s. We also compare cluster velocities with other dark energy probes proposed in the Dark Energy Task Force report. For cluster velocity measurements with realistic priors, the eventual constraints on the dark energy density, the dark energy equation of state and its evolution are comparable to constraints from supernovae measurements, and better than cluster counts and baryon acoustic oscillations; adding velocity to other dark energy probes improves constraints on the figure of merit by more than a factor of two. For upcoming Sunyaev-Zeldovich galaxy cluster surveys, even velocity measurements with errors as large as 1000 km/s will substantially improve the cosmological constraints compared to using the cluster number density alone.

Future microwave sky surveys will have the sensitivity to detect the kinematic Sunyaev-Zeldovich signal from moving galaxy clusters, thus providing a direct measurement of their line-of-sight peculiar velocity. We show that cluster peculiar velocity statistics applied to foreseeable surveys will put significant constraints on fundamental cosmological parameters. We consider three statistical quantities that can be constructed from a cluster peculiar velocity catalog: the probability density function, the mean pairwise streaming velocity, and the pairwise velocity dispersion. These quantities are applied to an envisioned data set which measures line-of-sight cluster velocities with normal errors of 100 km/s for all clusters with masses larger than $10^{14}$ solar masses over a sky area of up to 5000 square degrees. A simple Fisher matrix analysis of this survey shows that the normalization of the matter power spectrum and the dark energy equation of state can be constrained to better than 10 percent, and the Hubble constant and the primordial power spectrum index can be constrained to a few percent, independent of any other cosmological observations. We also find that the current constraint on the power spectrum normalization can be improved by more than a factor of two using data from a 400 square degree survey and WMAP third-year priors. We also show how the constraints on cosmological parameters changes if cluster velocities are measured with normal errors of 300 km/s.

Angular correlations of ultra-high energy cosmic rays with cosmologically distant sources may provide clues to these mysterious events. We compare cosmic ray tracks with energies above $ 10^{20}eV$ to a compilation of radio-loud compact QSO positions. The statistical method emphasizes invariant quantities and a test of {\it statistical independence} of track and source distributions. Statistical independence is ruled out by several independent statistics at confidence levels of less that $10^{-3}$ ($99.9% $.), indicating that track directions and QSO source positions are correlated at a highly significant level.