Numerical simulations are very helpful in understanding the physics of the formation of structure and galaxies. However, it is sometimes difficult to interpret model data with respect to observations, partly due to the difficulties and background noise inherent to observation. The goal, here, is to attempt to bridge this gap between simulation and observation by rendering the model output in image format which is then processed by tools commonly used in observational astronomy. Images are synthesized in various filters by folding the output of cosmological simulations of gasdynamics with star-formation and dark matter with the Bruzual- Charlot stellar population synthesis models. A variation of the Virgo-Gadget numerical simulation code is used with the hybrid gas and stellar formation models of Springel and Hernquist (2003). Outputs taken at various redshifts are stacked to create a synthetic view of the simulated star clusters. Source Extractor (SExtractor) is used to find groupings of stellar populations which are considered as galaxies or galaxy building blocks and photometry used to estimate the rest frame luminosities and distribution functions. With further refinements, this is expected to provide support for missions such as JWST, as well as to probe what additional physics are needed to model the data. The results show good agreement in many respects with observed properties of the galaxy luminosity function (LF) over a wide range of high redshifts. In particular, the slope (alpha) when fitted to the standard Schechter function shows excellent agreement both in value and evolution with redshift, when compared with observation. Discrepancies of other properties with observation are seen to be a result of limitations of the simulation and additional feedback mechanisms which are needed.

The first part of this dissertation presents the implementation of Bayesian statistics with galaxy surface luminosity (SL) prior probabilities to improve the ac- curacy of photometric redshifts. The addition of the SL prior probability helps break the degeneracy of spectro-photometric redshifts (SPZs) between low redshift 4000 A break galaxies and high redshift Lyman break galaxies which are mostly catas- trophic outliers. For a sample of 1138 galaxies with spectroscopic redshifts in the GOODS North and South fields at z < 1.6, the application of the surface luminosity prior reduces the fraction of galaxies with redshift deviation sigma(z) > 0.2 from 15.0% to 10.4%. The second part of this dissertation presents the study of the chemical evolution of the star-forming galaxies. The Hubble Space Telescope Probing Evolution and Reionization Spectroscopically (PEARS) grism Survey effectively selects emission line galaxies (ELGs) to mAB ~ 27. Follow-up Magellan LDSS3+IMACS spectroscopy of the HST/ACS PEARS ELGs confirms an accuracy of sigma_z = 0.006 for the HST/ACS PEARS grism redshifts. The luminosity-metallicity (L-Z) relation and the mass-metallicity (M-Z) relation of the PEARS ELGs at z ~ 0.6 are offset by ~ - 0.8 dex in metallicity for a given rest-frame B absolute magnitude and stellar mass relative to the local relations from SDSS galaxies. The offsets in both relations are ~ - 0.4 dex larger than that given by other samples at same redshifts, which are demonstrated to be due to the selection of different physical properties of the PEARS ELGs: low metallicities, very blue colors, small sizes, compact disturbed morphologies, high SSFR > 10^-9 yr^-1 , and high gas fraction. The downsizing effect, the tidal interacting induced inflow of metal-poor gas, and the SNe driven galactic winds outflows, may account for the significant offset of the PEARS galaxies in the L-Z and the M-Z relations relative to the local relations. The detection of the emission lines of ELGs down to m ~ 26 mag in the HST/ACS PEARS + HST/WCF3 ERS NIR composit grism spectra enables to extend the study of the evolution of the L-Z and M-Z relations to 0.6 < z < 2.4.

Using high-resolution three-dimensional adaptive mesh refinement simulations I study the interaction between primordial minihalo, a clump of baryonic and dark matter with a virial temperature below the atomic cooling limit, and a galaxy outflow. In Chapter 2 I concentrate on the formation of molecular coolants and their effect on the evolution of the minihalo gas. Molecular coolants are important since they allow gas to cool below 10000 K. Therefore, I implement a primordial chemistry and cooling network that tracks the evolution and cooling from these species. I show that the shock from the galaxy outflow produces an abundance of coolants in the primordial gas which allows the gas to cool to below 10000 K. I also show that this interaction produces compact stellar clusters that are ejected from their parent dark matter halos. In Chapter 3 I look at the turbulent mixing of metals that occur between the minihalo and outflow. To do this, I develop a sub-grid model for turbulence that reproduces three primary fluid instabilities. I find that the metals from the outflow are well mixed throughout the minihalo gas. In addition, the metal abundance found roughly corresponds to the observed abundances in halo globular clusters. In Chapter 4, I conduct a suite of simulations that follow this interaction over a wide range of parameters. In almost all cases, the shocked minihalos form molecules and cool rapidly to become compact, chemically homogenous stellar clusters. Furthermore, I show that the unique properties of these clusters make them a prime observational target for study with the next generation of telescopes. Given the unique properties of these clusters there are reasons to suspect that their low-redshift counterparts are halo globular clusters. I outline this comparison in Chapter 5 and give my conclusions in Chapter 6. Finally, I summarize my current work in Chapter 7 and future extensions in Chapter 8. By the end, I hope to convince you that the interaction between a galaxy outflow and a primordial minihalo provides a formation pathway for present day halo globular clusters.

Most stars form in groups, and these clusters are themselves nestled within larger associations and stellar complexes. It is not yet clear, however, whether stars cluster on preferred size scales within galaxies, or if stellar groupings have a continuous size distribution. I have developed two methods to select stellar groupings across a wide range of size-scales in order to assess trends in the size distribution and other basic properties of stellar groupings. The first method uses visual inspection of color-magnitude and color-color diagrams of clustered stars to assess whether the compact sources within the potential association are coeval, and thus likely to be born from the same parentmolecular cloud. This method was developed using the stellar associations in the M51/NGC 5195 interacting galaxy system. This process is highly effective at selecting single-aged stellar associations, but in order to assess properties of stellar clustering in a larger sample of nearby galaxies, an automated method for selecting stellar groupings is needed. I have developed an automated stellar grouping selection method that is sensitive to stellar clustering on all size scales. Using the Source Extractor software package on Gaussian-blurred images of NGC 4214, and the annular surface brightness to determine the characteristic size of each cluster/association, I eliminate much of the size and density biases intrinsic to other methods. This automated method was tested in the nearby dwarf irregular galaxy NGC 4214, and can detect stellar groupings with sizes ranging from compact clusters to stellar complexes. In future work, the automatic selection method developed in this dissertation will be used to identify stellar groupings in a set of nearby galaxies to determine if the size scales for stellar clustering are uniform in the nearby universe or if it is dependent on local galactic environment. Once the stellar clusters and associations have been identified and age-dated, this information can be used to deduce disruption times from the age distribution as a function of the position of the stellar grouping within the galaxy, the size of the cluster or association, and the morphological type of the galaxy. The implications of these results for galaxy formation and evolution are discussed.

There are many lines of evidence for anisotropy at all scales in the explosions of core collapse supernovae, e.g. visual inspection of the images of resolved supernova remnants, polarization measurements, velocity profiles, "natal kicks" of neutron stars, or spectroscopic observations of different regions of remnants. Theoretical stability considerations and detailed numerical simulations have shown that Rayleigh-Taylor (RT) instabilities arise in the star after the explosion, which leads to the early fragmentation of parts of the ejecta. The clumps thus created are of interest to a variety of topics, one of them being the formation environment of the solar system. There is a high probability that the solar system formed in the vicinity of a massive star that, shortly after its formation, exploded as a core collapse supernova. As argued in this thesis as well as other works, a core collapse supernova generally is a good candidate for chemically enriching the forming solar system with material. As forming proto--planetary systems in general have a high probability of being contaminated with supernova material, a method was developed for detecting tracer elements indicative supernova contamination in proto--planetary systems.The degree of the anisotropy of the supernova explosion can have dramatic effects on the mode of delivery of that material to the solar system, or proto--planetary systems in general. Thus it is of particular interest to be able to predict the structure of the supernova ejecta. Numerical simulations of the explosions of core collapse supernovae were done in 3 dimensions in order to study the formation of structure. It is found that RT instabilities result in clumps in the He- and C+O rich regions in the exploding star that are overdense by 1-2 orders of magnitude. These clumps are potential candidates for enriching the solar system with material. In the course of the further evolution of the supernova remnant, these RT clumps are likely to evolve into ejecta knots of the type observed in the Cassiopeia A supernova remnant.

Galaxies with strong Lyman-alpha (Lya) emission line (also called Lya galaxies or emitters) offer an unique probe of the epoch of reionization - one of the important phases when most of the neutral hydrogen in the universe was ionized. In addition, Lya galaxies at high redshifts are a powerful tool to study low-mass galaxy formation. Since current observations suggest that the reionization is complete by redshift z~ 6, it is therefore necessary to discover galaxies at z > 6, to use their luminosity function (LF) as a probe of reionization. I found five z = 7.7 candidate Lya galaxies with line fluxes > 7x10-18 erg/s/cm/2 , from three different deep near-infrared (IR) narrowband (NB) imaging surveys in a volume > 4x104Mpc3. From the spectroscopic followup of four candidate galaxies, and with the current spectroscopic sensitivity, the detection of only the brightest candidate galaxy can be ruled out at 5 sigma level. Moreover, these observations successfully demonstrate that the sensitivity necessary for both, the NB imaging as well as the spectroscopic followup of z~ 8 Lya galaxies can be reached with the current instrumentation. While future, more sensitive spectroscopic observations are necessary, the observed Lya LF at z = 7.7 is consistent with z = 6.6 LF, suggesting that the intergalactic medium (IGM) is relatively ionized even at z = 7.7, with neutral fraction xHI≤ 30%. On the theoretical front, while several models of Lya emitters have been developed, the physical nature of Lya emitters is not yet completely known. Moreover, multi-parameter models and their complexities necessitates a simpler model. I have developed a simple, single-parameter model to populate dark mater halos with Lya emitters. The central tenet of this model, different from many of the earlier models, is that the star-formation rate (SFR), and hence the Lya luminosity, is proportional to the mass accretion rate rather than the total halo mass. This simple model is successful in reproducing many observable including LFs, stellar masses, SFRs, and clustering of Lya emitters from z~ 3 to z~ 7. Finally, using this model, I find that the mass accretion, and hence the star-formation in > 30% of Lya emitters at z~ 3 occur through major mergers, and this fraction increases to ~ 50% at z~7.

Type Ia supernovae are important, but mysterious cosmological tools. Their standard brightnesses have enabled cosmologists to measure extreme distances and to discover dark energy. However, the nature of their progenitor mechanisms remains elusive, with many competing models offering only partial clues to their origins. Here, type Ia supernova delay times are explored using analytical models. Combined with a new observation technique, this model places new constraints on the characteristic time delay between the formation of stars and the first type Ia supernovae. This derived delay time (500 million years) implies low-mass companions for single degenerate progenitor scenarios. In the latter portions of this dissertation, two progenitor mechanisms are simulated in detail; white dwarf collisions and mergers. From the first of these simulations, it is evident that white dwarf collisions offer a viable and unique pathway to producing type Ia supernovae. Many of the combinations of masses simulated produce sufficient quantities of 56Ni (up to 0.51 solar masses) to masquerade as normal type Ia supernovae. Other combinations of masses produce 56Ni yields that span the entire range of supernova brightnesses, from the very dim and underluminous, with 0.14 solar masses, to the over-bright and superluminous, with up to 1.71 solar masses. The 56Ni yield in the collision simulations depends non-linearly on total system mass, mass ratio, and impact parameter. Using the same numerical tools as in the collisions examination, white dwarf mergers are studied in detail. Nearly all of the simulations produce merger remnants consisting of a cold, degenerate core surrounded by a hot accretion disk. The properties of these disks have strong implications for various viscosity treatments that have attempted to pin down the accretion times. Some mass combinations produce super-Chandrasekhar cores on shorter time scales than viscosity driven accretion. A handful of simulations also exhibit helium detonations on the surface of the primary that bear a resemblance to helium novae. Finally, some of the preliminary groundwork that has been laid for constructing a new numerical tool is discussed. This new tool advances the merger simulations further than any research group has done before, and has the potential to answer some of the lingering questions that the merger study has uncovered. The results of thermal diffusion tests using this tool have a remarkable correspondence to analytical predictions.