The main focus of this thesis is the investigation of an interesting new population of massive passively evolving galaxies found at high redshifts. We use a sample of these galaxies at redshifts 1<z<2 drawn from the Gemini Deep Deep Survey sample and measure their structural parameters based on the Hubble Space Telescope imaging in the rest frame visible and ultraviolet wavelength range. Our analysis shows that a fraction of these systems are very compact, with effective radii of R_e<1 kpc, even when observed in rest frame visible light. The average size of these objects is 2-5 times smaller than the typical size of an early-type galaxy of similar mass found locally. We combine the results from our morphological analysis with data from published spectroscopic samples of quiescent systems with known structural parameters. Analysis of these data for galaxies over the redshift range 0<z<2.5 shows that passively evolving galaxies are continuously and gradually growing in size. We also find smooth evolution of the stellar mass density within the central kiloparsec of these systems. The stellar mass density grows by a factor of 3 from z=0 to z~2.5. None of the models proposed to drive the structural evolution of early-type objects can explain all the observed aspects of this process. Because these massive compact galaxies have such small angular sizes, future studies of these systems will benefit from adaptive optics. In order to compile a large statistical sample of these objects suitable for adaptive optics follow-up, we first need to find a large number of targets with nearby bright natural guide stars. As a first step in this process, we describe the properties of a set of one square degree regions of the sky we have located that have a rare combination of high stellar surface density and low levels of extinction. We demonstrate that the adaptive optics-related properties of these fields are in some cases orders of magnitude better than those of existing deep fields.

Massive young star clusters and OB associations (M > 10 000 Msun) dominate the energetic feedback from stars into the interstellar medium. They contain the most massive and luminous stars in the Galaxy, which shape their environments through winds, ionizing flux, radiation pressure, and eventually supernovae, destroying their natal molecular clouds and inflating superbubbles. Few such clusters have been identified in our Galaxy. We systematically investigate the most luminous H II regions, which we identify using the WMAP foreground maps. We find that the 13 most luminous sources produce one-third of the Galaxy’s total ionizing luminosity, all with expected powering populations of M > 40 000 Msun. These populations are grouped in small numbers of clusters or associations for each WMAP source. The emission from these regions is dominated by the diffuse component at large radii (∼10-70 pc) indicating a high leaking fraction of ionizing photons. Using 8 μm maps from Spitzer GLIMPSE and published radio recombination line observations, we resolve the large (> 1◦) WMAP sources into 40 star forming complexes (SFCs) exhibiting shell morphology with evidence of expansion due to a central powering source. We develop a method, based on differential extinction of the galactic disk, to identify the SFC’s powering cluster candidates with 2MASS. We identify 22 cluster candidates within the 40 SFCs having extinctions consistent with their distances. With near-infrared spectroscopy from the New Technology Telescope, we have confirmed the existence of the most massive of these associations, the Dragonfish Association, with M = 100 000 Msun. Of the 50 sampled stars, we identify 2 Luminous Blue Variable candidates, a Wolf-Rayet, and 15 O-type stars, consistent with the yield expected from the candidate contamination rate, verifying the candidate cluster identification method. This investigation doubles the number of massive young star clusters and OB associations known and produces the most complete picture of the upper end of the Galaxy’s cluster mass function to date.

The main focus of this thesis is the investigation of an interesting new population of massive passively evolving galaxies found at high redshifts. We use a sample of these galaxies at redshifts 1<z<2 drawn from the Gemini Deep Deep Survey sample and measure their structural parameters based on the Hubble Space Telescope imaging in the rest frame visible and ultraviolet wavelength range. Our analysis shows that a fraction of these systems are very compact, with effective radii of R_e<1 kpc, even when observed in rest frame visible light. The average size of these objects is 2-5 times smaller than the typical size of an early-type galaxy of similar mass found locally. We combine the results from our morphological analysis with data from published spectroscopic samples of quiescent systems with known structural parameters. Analysis of these data for galaxies over the redshift range 0<z<2.5 shows that passively evolving galaxies are continuously and gradually growing in size. We also find smooth evolution of the stellar mass density within the central kiloparsec of these systems. The stellar mass density grows by a factor of 3 from z=0 to z~2.5. None of the models proposed to drive the structural evolution of early-type objects can explain all the observed aspects of this process. Because these massive compact galaxies have such small angular sizes, future studies of these systems will benefit from adaptive optics. In order to compile a large statistical sample of these objects suitable for adaptive optics follow-up, we first need to find a large number of targets with nearby bright natural guide stars. As a first step in this process, we describe the properties of a set of one square degree regions of the sky we have located that have a rare combination of high stellar surface density and low levels of extinction. We demonstrate that the adaptive optics-related properties of these fields are in some cases orders of magnitude better than those of existing deep fields.

The existence of Gravitational Waves from binary black holes is one of the most interesting predictions of General Relativity. These ripples in space-time should be visible to ground-based gravitational wave detectors worldwide in the next few years. One such detector, the Laser Interferometer Gravitational-wave Observatory (LIGO) is in the process of being upgraded to its Advanced sensitivity which should make gravitational wave detections routine. Even so, the signals that LIGO will detect will be faint compared to the detector noise, and so accurate waveform templates are crucial. In this thesis, we present a detailed analysis of the accuracy of hybrid gravitational waveforms. Hybrids are created by stitching a long post-Newtonian inspiral to the late inspiral, merger, and ringdown produced by numerical relativity simulations. We begin our investigation with a study of the systematic errors in the numerical waveform, and errors due to hybridization and choice of detector noise. For current NR waveforms, the largest source of error comes from the unknown high-order terms in the post-Newtonian waveform, which we first explore for equal-mass, non-spinning binaries, and also for unequal-mass, non-spinning binaries. We then consider the potential reduction in hybrid errors if these higher-order terms were known. Finally, we investigate the possibility of using hybrid waveforms as a detection template bank and integrating NR+PN hybrids into the LIGO detection pipeline.

The existence of Gravitational Waves from binary black holes is one of the most interesting predictions of General Relativity. These ripples in space-time should be visible to ground-based gravitational wave detectors worldwide in the next few years. One such detector, the Laser Interferometer Gravitational-wave Observatory (LIGO) is in the process of being upgraded to its Advanced sensitivity which should make gravitational wave detections routine. Even so, the signals that LIGO will detect will be faint compared to the detector noise, and so accurate waveform templates are crucial. In this thesis, we present a detailed analysis of the accuracy of hybrid gravitational waveforms. Hybrids are created by stitching a long post-Newtonian inspiral to the late inspiral, merger, and ringdown produced by numerical relativity simulations. We begin our investigation with a study of the systematic errors in the numerical waveform, and errors due to hybridization and choice of detector noise. For current NR waveforms, the largest source of error comes from the unknown high-order terms in the post-Newtonian waveform, which we first explore for equal-mass, non-spinning binaries, and also for unequal-mass, non-spinning binaries. We then consider the potential reduction in hybrid errors if these higher-order terms were known. Finally, we investigate the possibility of using hybrid waveforms as a detection template bank and integrating NR+PN hybrids into the LIGO detection pipeline.

The method of baryon acoustic oscillation (BAO) is among the best probes of the dark energy equation of state, and worldwide efforts are being invested in order to perform measurements that are accurate at the percent level. In current data analyses, however, estimates of the error about the BAO are based on the assumption that the density field can be treated as Gaussian, an assumption that becomes less accurate as smaller scales are included in the measurement. It was recently shown from large samples of N-body simulations that the error bars about the BAO obtained this way are in fact up to 15-20 per cent too small. This important bias has shaken the confidence in the way error bars are calculated, and is motivating developments of analyses pipelines that include non-Gaussian features in the matter density fields. In this thesis, we propose general strategies to incorporate non-Gaussian effects in the context of a survey. After describing the high performance N-body code that we used, we present novel properties of the non-Gaussian uncertainty about the matter power spectrum, and explain how these combine with a general survey selection function. Assuming that the non-Gaussian features that are observed in the simulations correspond to those of Nature, this approach is the first unbiased measurement of the error bar about the power spectrum, which simultaneously removes the undesired bias on the BAO error. We then relax this assumption about the similitude of the non-Gaussian natures in simulations and data, and develop tools that aim at measuring the non-Gaussian error bars exclusively from the data. It is possible to improve the constraining power of non-Gaussian analyses with `Gaussianizations' techniques, which map the observed fields into something more Gaussian. We show that two of such techniques maximally recover degrees of freedom that were lost in the gravitational collapse. Finally, from a large sample of high resolution N-body realizations, we construct a series of weak gravitational lensing distortion maps and provide high resolution halo catalogues that are used by the CFTHLenS community to calibrate their estimators and study many secondary effects with unprecedented accuracy.

In this thesis, we examine the role of accretion disks in the formation of stellar systems, focusing on young massive disks which regulate the flow of material from the parent molecular core down to the star. We study the evolution of disks with high infall rates that develop strong gravitational instabilities. We begin in chapter 1 with a review of the observations and theory which underpin models for the earliest phases of star formation and provide a brief review of basic accretion disk physics, and the numerical methods which we employ. In chapter 2 we outline the current models of binary and multiple star formation, and review their successes and shortcomings from a theoretical and observational perspective. In chapter 3 we begin with a relatively simple analytic model for disks around young, very massive stars, showing that instability in these disks may be responsible for the higher multiplicity fraction of massive stars, and perhaps the upper mass to which they grow. We extend these models in chapter 4 to explore the properties of disks and the formation of binary companions across a broad range of stellar masses. In particular, we model the role of global and local mechanisms for angular momentum transport in regulating the relative masses of disks and stars. We follow the evolution of these disks throughout the main accretion phase of the system, and predict the trajectory of disks through parameter space. We follow up on the predictions made in our analytic models with a series of high resolution, global numerical experiments in chapter 5. Here we propose and test a new parameterization for describing rapidly accreting, gravitationally unstable disks. We find that disk properties and system multiplicity can be mapped out well in this parameter space. Finally, in chapter 6, we address whether our studies of unstable disks are relevant to recently detected massive planets on wide orbits around their central stars.

Massive young star clusters and OB associations (M > 10 000 Msun) dominate the energetic feedback from stars into the interstellar medium. They contain the most massive and luminous stars in the Galaxy, which shape their environments through winds, ionizing flux, radiation pressure, and eventually supernovae, destroying their natal molecular clouds and inflating superbubbles. Few such clusters have been identified in our Galaxy. We systematically investigate the most luminous H II regions, which we identify using the WMAP foreground maps. We find that the 13 most luminous sources produce one-third of the Galaxy’s total ionizing luminosity, all with expected powering populations of M > 40 000 Msun. These populations are grouped in small numbers of clusters or associations for each WMAP source. The emission from these regions is dominated by the diffuse component at large radii (∼10-70 pc) indicating a high leaking fraction of ionizing photons. Using 8 μm maps from Spitzer GLIMPSE and published radio recombination line observations, we resolve the large (> 1◦) WMAP sources into 40 star forming complexes (SFCs) exhibiting shell morphology with evidence of expansion due to a central powering source. We develop a method, based on differential extinction of the galactic disk, to identify the SFC’s powering cluster candidates with 2MASS. We identify 22 cluster candidates within the 40 SFCs having extinctions consistent with their distances. With near-infrared spectroscopy from the New Technology Telescope, we have confirmed the existence of the most massive of these associations, the Dragonfish Association, with M = 100 000 Msun. Of the 50 sampled stars, we identify 2 Luminous Blue Variable candidates, a Wolf-Rayet, and 15 O-type stars, consistent with the yield expected from the candidate contamination rate, verifying the candidate cluster identification method. This investigation doubles the number of massive young star clusters and OB associations known and produces the most complete picture of the upper end of the Galaxy’s cluster mass function to date.

The method of baryon acoustic oscillation (BAO) is among the best probes of the dark energy equation of state, and worldwide efforts are being invested in order to perform measurements that are accurate at the percent level. In current data analyses, however, estimates of the error about the BAO are based on the assumption that the density field can be treated as Gaussian, an assumption that becomes less accurate as smaller scales are included in the measurement. It was recently shown from large samples of N-body simulations that the error bars about the BAO obtained this way are in fact up to 15-20 per cent too small. This important bias has shaken the confidence in the way error bars are calculated, and is motivating developments of analyses pipelines that include non-Gaussian features in the matter density fields. In this thesis, we propose general strategies to incorporate non-Gaussian effects in the context of a survey. After describing the high performance N-body code that we used, we present novel properties of the non-Gaussian uncertainty about the matter power spectrum, and explain how these combine with a general survey selection function. Assuming that the non-Gaussian features that are observed in the simulations correspond to those of Nature, this approach is the first unbiased measurement of the error bar about the power spectrum, which simultaneously removes the undesired bias on the BAO error. We then relax this assumption about the similitude of the non-Gaussian natures in simulations and data, and develop tools that aim at measuring the non-Gaussian error bars exclusively from the data. It is possible to improve the constraining power of non-Gaussian analyses with `Gaussianizations' techniques, which map the observed fields into something more Gaussian. We show that two of such techniques maximally recover degrees of freedom that were lost in the gravitational collapse. Finally, from a large sample of high resolution N-body realizations, we construct a series of weak gravitational lensing distortion maps and provide high resolution halo catalogues that are used by the CFTHLenS community to calibrate their estimators and study many secondary effects with unprecedented accuracy.

Master of Science / Department of Physics / Bharat Ratra / In this report, we discuss the four standard energy conditions of General Relativity (null, weak, dominant, and strong) and investigate their cosmological consequences. We note that these energy conditions can be compatible with cosmic acceleration provided that a repulsive cosmological constant exists and the acceleration stays within certain bounds. Scalar fields and dark energy, and their relationships to the energy conditions, are also discussed. Special attention is paid to the 1988 Ratra-Peebles scalar field model, which is notable in that it provides a physical self-consistent framework for the phenomenology of dark energy. Appendix B, which is part of joint-research with Anatoly Pavlov, Khaled Saaidi, and Bharat Ratra, reports on the existence of the Ratra-Peebles scalar field tracker solution in a curvature-dominated universe, and discusses the problem of investigating the evolution of long-wavelength inhomogeneities in this solution while taking into account the gravitational back-reaction (in the linear perturbative approximation).