Recent Lyman-$\alpha$ forest tomography measurements of the intergalactic medium (IGM) have revealed a wealth of cosmic structures at high redshift ($z\sim 2.5$). In this work, we present the Tomographic Absorption Reconstruction and Density Inference Scheme (TARDIS), a new chrono-cosmographic analysis tool for understanding the formation and evolution of these observed structures. We use maximum likelihood techniques with a fast non-linear gravitational model to reconstruct the initial density field of the observed regions. We find that TARDIS allows accurate reconstruction of smaller scale structures than standard Wiener filtering techniques. Applying this technique to mock Lyman-$\alpha$ forest data sets that simulate ongoing and future surveys such as CLAMATO, Subaru-PFS or the ELTs, we are able to infer the underlying matter density field at observed redshift and classify the cosmic web structures. We find good agreement with the underlying truth both in the characteristic eigenvalues and eigenvectors of the pseudo-deformation tensor, with the eigenvalues inferred from 30m-class telescopes correlated at $r=0.95$ relative to the truth. As an output of this method, we are able to further evolve the inferred structures to late time ($z=0$), and also track the trajectories of coeval $z=2.5$ galaxies to their $z=0$ cosmic web environments.

We present an efficient implementation of Wiener filtering of real-space linear field and optimal quadratic estimator of its power spectrum Band-powers. We first recast the field reconstruction into an optimization problem, which we solve using quasi-Newton optimization. We then recast the power spectrum estimation into the field marginalization problem, from which we obtain an expression that depends on the field reconstruction solution and a determinant term. We develop a novel simulation based method for the latter. We extend the simulations formalism to provide the covariance matrix for the power spectrum. We develop a flexible framework that can be used on a variety of cosmological fields and present results for a variety of test cases, using simulated examples of projected density fields, projected shear maps from galaxy lensing, and observed Cosmic Microwave Background (CMB) temperature anisotropies, with a wide range of map incompleteness and variable noise. For smaller cases where direct numerical inversion is possible, we show that our solution matches that created by direct Wiener Filtering at a fraction of the overall computation cost. Even more significant reduction of computational is achieved by this implementation of optimal quadratic estimator due to the fast evaluation of the Hessian matrix. This technique allows for accurate map and power spectrum reconstruction with complex masks and nontrivial noise properties.

Cosmic Microwave Background (CMB) lensing is a powerful probe of the matter distribution in the Universe. The standard quadratic estimator, which is typically used to measure the lensing signal, is known to be suboptimal for low-noise polarization data from next-generation experiments. In this paper we explain why the quadratic estimator will also be suboptimal for measuring lensing on very small scales, even for measurements in temperature where this estimator typically performs well. Though maximum likelihood methods could be implemented to improve performance, we explore a much simpler solution, revisiting a previously proposed method to measure lensing which involves a direct inversion of the background gradient. An important application of this simple formalism is the measurement of cluster masses with CMB lensing. We find that directly applying a gradient inversion matched filter to simulated lensed images of the CMB can tighten constraints on cluster masses compared to the quadratic estimator. While the difference is not relevant for existing surveys, for future surveys it can translate to significant improvements in mass calibration for distant clusters, where galaxy lensing calibration is ineffective due to the lack of enough resolved background galaxies. Improvements can be as large as $\sim 50\%$ for a cluster at $z = 2$ and a next-generation CMB experiment with 1$\mu$K-arcmin noise, and over an order of magnitude for lower noise levels. For future surveys, this simple matched-filter or gradient inversion method approaches the performance of maximum likelihood methods, at a fraction of the computational cost.

Thermal Sunyaev-Zeldovich (tSZ) power spectrum is one of the most sensitive methods to constrain cosmological parameters, scaling as the amplitude $\sigma_8^8$. It is determined by the integral over the halo mass function multiplied by the total pressure content of clusters, and further convolved by the cluster gas pressure profile. It has been shown that various feedback effects can change significantly the pressure profile, possibly even pushing the gas out to the virial radius and beyond, strongly affecting the tSZ power spectrum at high $l$. Energetics arguments and SZ-halo mass scaling relations suggest feedback is unlikely to significantly change the total pressure content, making low $l$ tSZ power spectrum more robust against feedback effects. Furthermore, the separation between the cosmic infrared background (CIB) and tSZ is more reliable at low $l$. Low $l$ modes are however probing very small volumes, giving rise to very large non-gaussian sampling variance errors. By computing the trispectrum contribution we identify $90<l<350$ as the minimum variance scale where the combined error is minimized. Using the measurement at this $l$ provides constraints that are almost as strong as using the entire power spectrum. We extend the work by the Planck collaboration by including a full trispectrum and feedback effects in the analysis. We perform a Markov-Chain Monte Carlo analysis over the two dimensional parameter space and find constraints on $\sigma_8$ by marginalizing over the feedback nuisance parameter.We obtain $\sigma_8 =0.820^{+0.021}_{-0.009}\left(\Omega_m/0.31\right)^{0.4}$, when fixing other parameters to Planck cosmology values. Our results suggest it is possible to obtain competitive cosmological constraints from tSZ without cluster redshift information, and that the current tSZ power spectrum shows no evidence for a low amplitude of $\sigma_8$.

Much attention has been drawn to the recent discoveries by the Advanced Laser Interferometer Gravitational-Wave Observatory (LIGO) of merging intermediate mass black holes. Of particular interest is the possibility that the merger events detected could be evidence of dark matter in the form of primordial black holes (PBHs). It has been argued that the presence of many black holes would effect the thermal and ionization history of the universe via their accretion of matter which would have strong signatures in the Cosmic Microwave Background's (CMB) power spectra evident in the damping of anisotropies and change in low-$l$ polarization power. In general the accretion is quite sensitive to the specific physics involved and the conditions of the early universe. In this work, we take a minimal approach and find constraints on PBHs not including the model dependent effects of nonlinear structure of formation or transition between different accretion models which would work to increase the effect. In addition, we include the relative velocity between dark matter and baryonic matter including the effects of supersonic streaming at high redshift which work to significantly reduce the constraining power. We also examine the constraints on more astrophysically-motivated extended black hole mass functions and discuss how mergers might effect this distribution. We find constraints on PBHs in the range $ \approx 30 M_\odot$, finding that they could not compose more than $10\%$ of the total dark matter content.

A search for RR Lyrae stars (RRLS) in 840 sq.deg. of the sky in right ascension 150 - 210 deg and declination -10 - +10 deg yielded 1013 type ab and 359 type c RRLS. This sample is used to study the density profile of the Galactic halo, halo substructures, and the Oosterhoff type of the halo over distances from 5 to 80 kpc. The halo is flattened towards the Galactic plane, and its density profile steepens in slope at galactocentric distances greater than 25 kpc. The RRLS in the stellar stream from the Sagittarius dwarf spheroidal (dSph) galaxy match well the model of Law and Majewski for the stars that were stripped 1.3 to 3.2 Gyr ago, but not for the ones stripped 3.2 to 5.0 Gyr ago. Over densities are found at the locations of the Virgo Overdensity and the Virgo Stellar Stream. Within 1 deg of 1220-1, which Jerjen et al. identify as a halo substructure at a distance of 24 kpc, there are 4 RRLS that are possibly members. Away from substructures, the RRLS are a mixture of Oosterhoff types I and II, but mostly OoI (73%). The accretion of galaxies resembling in RRLS content the most massive Milky Way satellites (LMC, SMC, For,Sgr) may explain this preponderance of OoI. Six new RRLS and 3 new anomalous Cepheids were found in the Sextans dSph galaxy.

We describe the La Silla-QUEST (LSQ) Variability Survey. LSQ is a dedicated wide-field synoptic survey in the Southern Hemisphere, focussing on the discovery and study of transients ranging from low redshift (z < 0.1) SN Ia, Tidal Disruption events, RR Lyr{\ae} variables, CVs, Quasars, TNOs and others. The survey utilizes the 1.0-m Schmidt Telescope of the European Southern Observatory at La Silla, Chile, with the large-area QUEST camera, a mosaic of 112 CCDs with field of view of 9.6 square degrees. The LSQ Survey was commissioned in 2009, and is now regularly covering ~1000 square deg per night with a repeat cadence of hours to days. The data are currently processed on a daily basis. We present here a first look at the photometric capabilities of LSQ and we discuss some of the most interesting recent transient detections.