Search

We investigate the effects of stellar populations and sizes on Ly$\alpha$ escape in 27 spectroscopically confirmed and 35 photometric Lyman-Alpha Emitters (LAEs) at z $\approx$ 2.65 in seven fields of the Bo\"otes region of the NOAO Deep Wide-Field Survey. We use deep $HST$/WFC3 imaging to supplement ground-based observations and infer key galaxy properties. Compared to typical star-forming galaxies (SFGs) at similar redshifts, the LAEs are less massive ($M_{\star} \approx 10^{7} - 10^{9}~M_{\odot}$), younger (ages $\lesssim$ 1 Gyr), smaller ($r_{e} <$ 1 kpc), less dust-attenuated (E(B$-$V) $\le$ 0.26 mag), but have comparable star-formation-rates (SFRs $\approx 1 - 100~M_{\odot} {\rm yr^{-1}}$). Some of the LAEs in the sample may be very young galaxies having low nebular metallicities (${\rm Z_{neb} \lesssim 0.2 Z_{\odot}}$) and/or high ionization parameters ($\log{(\rm U)} \gtrsim -2.4$). Motivated by previous studies, we examine the effects of the concentration of star formation and gravitational potential on Ly$\alpha$ escape, by computing star-formation-rate surface density, $\Sigma_{\rm SFR}$ and specific star-formation-rate surface density, $\Sigma_{\rm sSFR}$. For a given $\Sigma_{\rm SFR}$, the Ly$\alpha$ escape fraction is higher for LAEs with lower stellar masses. LAEs have higher $\Sigma_{\rm sSFR}$ on average compared to SFGs. Our results suggest that compact star formation in a low gravitational potential yields conditions amenable to the escape of Ly$\alpha$ photons. These results have important implications for the physics of Ly$\alpha$ radiative transfer and for the type of galaxies that may contribute significantly to cosmic reionization., Comment: 36 pages, 15 figures; Accepted for publication in The Astronomical Journal

We study the process of panspermia in Milky Way-like galaxies by modeling the probability of successful travel of organic compounds between stars harboring potentially habitable planets. To this end, we apply the modified habitability recipe of Gobat & Hong (2016) to a model galaxy from the MUGS suite of zoom-in cosmological simulations. We find that, unlike habitability, which only occupies narrow dynamic range over the entire galaxy, the panspermia probability can vary be orders of magnitude between the inner ($R, b = 1-4 {\rm kpc}$) and outer disk. However, only a small fraction of star particles have very large values of panspermia probability and, consequently, the fraction of star particles where the panspermia process is more effective than prebiotic evolution is much lower than from na\"ive expectations based on the ratio between panspermia probability and natural habitability., Comment: 22 pages, 18 figures, ApJ accepted

By utilizing large-scale graph analytic tools implemented in the modern Big Data platform, Apache Spark, we investigate the topological structure of gravitational clustering in five different universes produced by cosmological $N$-body simulations with varying parameters: (1) a WMAP 5-year compatible $\Lambda$CDM cosmology, (2) two different dark energy equation of state variants, and (3) two different cosmic matter density variants. For the Big Data calculations, we use a custom build of stand-alone Spark/Hadoop cluster at Korea Institute for Advanced Study (KIAS) and Dataproc Compute Engine in Google Cloud Platform (GCP) with the sample size ranging from 7 millions to 200 millions. We find that among the many possible graph-topological measures, three simple ones: (1) the average of number of neighbors (the so-called average vertex degree) $\alpha$, (2) closed-to-connected triple fraction (the so-called transitivity) $\tau_\Delta$, and (3) the cumulative number density $n_{s\ge5}$ of subcomponents with connected component size $s \ge 5$, can effectively discriminate among the five model universes. Since these graph-topological measures are in direct relation with the usual $n$-points correlation functions of the cosmic density field, graph-topological statistics powered by Big Data computational infrastructure opens a new, intuitive, and computationally efficient window into the dark Universe., Comment: 16 pages, 11 figures, submitted to MNRAS

We investigate the spatial distribution of Lyman alpha emitting galaxies (LAEs) at $z \approx 2.67$, selected from the NOAO Deep Wide-Field Survey (NDWFS), using two-point statistics and topological diagnostics adopted from network science. We measure the clustering length, $r_0 \approx 4 h^{-1}$ Mpc, and the bias, $b_{LAE} = 2.2^{+0.2}_{-0.1}$. Fitting the clustering with halo occupation distribution (HOD) models results in two disparate possibilities: (1) where the fraction of central galaxies is $<$1% in halos of mass $>10^{12}$$M_\odot$; and (2) where the fraction is $\approx$20%. We refer to these two scenarios as the `Dusty Core Scenario' for Model#1 since most of central galaxies in massive halos are dead in Ly$\alpha$ emission, and the `Pristine Core Scenario' for Model#2 since the central galaxies are bright in Ly$\alpha$ emission. Traditional two-point statistics cannot distinguish between these disparate models given the current data sets. To overcome this degeneracy, we generate mock catalogs for each HOD model using a high resolution $N$-body simulation and adopt a network statistics approach, which provides excellent topological diagnostics for galaxy point distributions. We find three topological anomalies from the spatial distribution of observed LAEs, which are not reproduced by the HOD mocks. We find that Model#2 matches better all network statistics than Model#1, suggesting that the central galaxies in $> 10^{12} h^{-1} M_\odot$ halos at $z \approx 2.67$ need to be less dusty to be bright as LAEs, potentially implying some replenishing channels of pristine gas such as the cold mode accretion., Comment: 23 pages, 18 figures, accepted by MNRAS

We present a detailed census of galaxies in and around PC217.96+32.3, a spectroscopically confirmed Coma analog at z=3.78. Diverse galaxy types identified in the field include Lya emitters (LAEs), massive star-forming galaxies, and ultra-massive galaxies (log (Mstar/Msun)>= 11) which may have already halted their star formation. The sky distribution of the star-forming galaxies suggests the presence of a significant overdensity (delta_g=8+/-2), which is spatially offset from the previously confirmed members by 3-4 Mpc to the west. Candidate quiescent and post-starburst galaxies are also found in large excess (a factor of ~ 8-15 higher surface density than the field) although their redshifts are less certain. We estimate that the total enclosed mass traced by star-forming galaxy candidates is roughly comparable to that of PC217.96+32.3 traced by the LAEs. We speculate that the true extent of PC217.96+32.3 may be larger than previously known, a half of which is missed by our LAE selection. Alternatively, the newly discovered overdensity may belong to another Coma progenitor not associated with PC217.96+32.3. Expectations from theory suggest that both scenarios are equally unlikely (<1%), particularly in the cosmic volume probed in our survey. If confirmed as a single structure, its total mass will be well in excess of Coma, making this a singularly large cosmic structure rarely seen even in large cosmological simulations. Finally, we find that the protocluster galaxies follow the same SFR-M_star scaling relation as the field galaxies, suggesting that the environmental effect at z~4 is a subtle one at best for normal star-forming galaxies., Comment: 21 pages, 11 figures, accepted for publication in Astrophysical Journal

We report the discovery and spectroscopic confirmation of a very large star-forming Lyman Break galaxy, G6025, at z_spec=3.721+/-0.003. In the rest-frame ~2100A, G6025 subtends ~24 kpc in physical extent when measured from the 1.5-sigma isophote, in agreement with the parametric size measurements which yield the half-light radius of 4.9+/-0.5 kpc and the semi-major axis of 12.5+/-0.1 kpc. G6025 is also very UV-luminous (~5L*(z~4}) and young (~140+/-60 Myr). Despite its unusual size and luminosity, the stellar population parameters and dust reddening (M_star~M*(z~4)$, and E(B-V)=0.18+/-0.05) estimated from the integrated light, are similar to those of smaller galaxies at comparable redshifts. The ground-based morphology and spectroscopy show two dominant components, both located off-center, embedded in more diffuse emission. We speculate that G6025 may be a scaled-up version of chain galaxies seen in deep HST imaging, or alternatively, a nearly equal-mass merger involving two super-L* galaxies in its early stage. G6025 lies close to but not within a known massive protocluster at z=3.78. We find four companions within 6 Mpc from G6025, two of which lie within 1.6 Mpc. While the limited sensitivity of the existing spectroscopy does not allow us to robustly characterize the local environment of G6025, it likely resides in a locally overdense environment. The luminosity, size, and youth of G6025 make it uniquely suited to study the early formation of massive galaxies in the universe., Comment: accepted for publication

We report the detection of diffuse Ly$\alpha$ emission, or Ly$\alpha$ halos (LAHs), around star-forming galaxies at $z\approx3.78$ and $2.66$ in the NOAO Deep Wide-Field Survey Bo\"otes field. Our samples consist of a total of $\sim$1400 galaxies, within two separate regions containing spectroscopically confirmed galaxy overdensities. They provide a unique opportunity to investigate how the LAH characteristics vary with host galaxy large-scale environment and physical properties. We stack Ly$\alpha$ images of different samples defined by these properties and measure their median LAH sizes by decomposing the stacked Ly$\alpha$ radial profile into a compact galaxy-like and an extended halo-like component. We find that the exponential scale-length of LAHs depends on UV continuum and Ly$\alpha$ luminosities, but not on Ly$\alpha$ equivalent widths or galaxy overdensity parameters. The full samples, which are dominated by low UV-continuum luminosity Ly$\alpha$ emitters ($M_{\rm UV} \gtrsim -21$), exhibit LAH sizes of 5$\,-\,6\,$kpc. However, the most UV- or Ly$\alpha$-luminous galaxies have more extended halos with scale-lengths of 7$\,-\,9\,$kpc. The stacked Ly$\alpha$ radial profiles decline more steeply than recent theoretical predictions that include the contributions from gravitational cooling of infalling gas and from low-level star formation in satellites. On the other hand, the LAH extent matches what one would expect for photons produced in the galaxy and then resonantly scattered by gas in an outflowing envelope. The observed trends of LAH sizes with host galaxy properties suggest that the physical conditions of the circumgalactic medium (covering fraction, HI column density, and outflow velocity) change with halo mass and/or star-formation rates., Comment: published in ApJ, minor proof corrections applied

We present new observations of the field containing the z=3.786 protocluster, PC217.96+32.3. We confirm that it is one of the largest and most overdense high-redshift structures known. Such structures are rare even in the largest cosmological simulations. We used the Mayall/MOSAIC1.1 imaging camera to image a 1.2x0.6 deg area (~150x75 comoving Mpc) surrounding the protocluster's core and discovered 165 candidate Lyman Alpha emitting galaxies (LAEs) and 788 candidate Lyman Break galaxies (LBGs). There are at least 2 overdense regions traced by the LAEs, the largest of which shows an areal overdensity in its core (i.e., within a radius of 2.5 comoving Mpc) of 14+/-7 relative to the average LAE spatial density in the imaged field. Further, the average LAE spatial density in the imaged field is twice that derived by other field LAE surveys. Spectroscopy with Keck/DEIMOS yielded redshifts for 164 galaxies (79 LAEs and 85 LBGs); 65 lie at a redshift of 3.785+/-0.010. The velocity dispersion of galaxies near the core is 350+/-40 km/s, a value robust to selection effects. The overdensities are likely to collapse into systems with present-day masses of >10^{15} solar masses and >6x10^{14} solar masses. The low velocity dispersion may suggest a dynamically young protocluster. We find a weak trend between narrow-band (Lyman Alpha) luminosity and environmental density: the Lyman Alpha luminosity is enhanced on average by 1.35X within the protocluster core. There is no evidence that the Lyman Alpha equivalent width depends on environment. These suggest that star-formation and/or AGN activity is enhanced in the higher density regions of the structure. PC217.96+32.3 is a Coma cluster analog, witnessed in the process of formation., Comment: Accepted for publication in the Astrophysical Journal (March 27, 2016)

The concept of the cosmic web, viewing the Universe as a set of discrete galaxies held together by gravity, is deeply engrained in cosmology. Yet, little is known about the most effective construction and the characteristics of the underlying network. Here we explore seven network construction algorithms that use various galaxy properties, from their location, to their size and relative velocity, to assign a network to galaxy distributions provided by both simulations and observations. We find that a model relying only on spatial proximity offers the best correlations between the physical characteristics of the connected galaxies. We show that the properties of the networks generated from simulations and observations are identical, unveiling a deep universality of the cosmic web., Comment: 5 pages, 4 figures

(abridged) The large-scale distribution of galaxies is generally analyzed using the two-point correlation function. However, this statistic does not capture the topology of the distribution, and it is necessary to resort to higher order correlations to break degeneracies. We demonstrate that an alternate approach using network analysis can discriminate between topologically different distributions that have similar two-point correlations. We investigate two galaxy point distributions, one produced by a cosmological simulation and the other by a L\'evy walk. For the cosmological simulation, we adopt the redshift $z = 0.58$ slice from Illustris (Vogelsberger et al. 2014A) and select galaxies with stellar masses greater than $10^8$$M_\odot$. The two point correlation function of these simulated galaxies follows a single power-law, $\xi(r) \sim r^{-1.5}$. Then, we generate L\'evy walks matching the correlation function and abundance with the simulated galaxies. We find that, while the two simulated galaxy point distributions have the same abundance and two point correlation function, their spatial distributions are very different; most prominently, \emph{filamentary structures}, absent in L\'evy fractals. To quantify these missing topologies, we adopt network analysis tools and measure diameter, giant component, and transitivity from networks built by a conventional friends-of-friends recipe with various linking lengths. Unlike the abundance and two point correlation function, these network quantities reveal a clear separation between the two simulated distributions; therefore, the galaxy distribution simulated by Illustris is not a L\'evy fractal quantitatively. We find that the described network quantities offer an efficient tool for discriminating topologies and for comparing observed and theoretical distributions., Comment: 11 pages, 5 figures, submitted to MNRAS on 12/15/2015, now fully reviewed for publication; more information about our network analyses can be found at https://sites.google.com/site/shongscience/research

We apply simple analyses techniques developed for the study of complex networks to the study of the cosmic web, the large scale galaxy distribution. In this paper, we measure three network centralities (ranks of topological importance), Degree Centrality (DC), Closeness Centrality (CL), and Betweenness Centrality (BC) from a network built from the Cosmological Evolution Survey (COSMOS) catalog. We define 8 galaxy populations according to the centrality measures; Void, Wall, and Cluster by DC, Main Branch and Dangling Leaf by BC, and Kernel, Backbone, and Fracture by CL. We also define three populations by voronoi tessellation density to compare these with the DC selection. We apply the topological selections to galaxies in the (photometric) redshift range $0.91

We report the discovery of a large-scale structure containing multiple protoclusters at z=3.78 in the Bo\"otes field. The spectroscopic discovery of five galaxies at z=3.783+/-0.002 lying within 1 Mpc of one another led us to undertake a deep narrow- and broad-band imaging survey of the surrounding field. Within a comoving volume of 72x72x25 Mpc^3, we have identified 65 Lyman alpha emitter (LAE) candidates at z=3.795+/-0.015, and four additional galaxies at z_spec=3.730,3.753,3.780,3.835. The galaxy distribution within the field is highly non-uniform, exhibiting three large (~3-5x) overdensities separated by 8-14 Mpc (physical) and possibly connected by filamentary structures traced by LAEs. The observed number of LAEs in the entire field is nearly twice the average expected in field environments, based on estimates of the Lya luminosity function at these redshifts. We estimate that by z=0 the largest overdensity will grow into a cluster of mass 10^15 Msun; the two smaller overdensities will grow into clusters of mass (2-6)x10^14 Msun. The highest concentration of galaxies is located at the southern end of the image, suggesting that the current imaging may not map the true extent of the large scale structure. Finding three large protocluster candidates within a single 0.3 deg^2 field is highly unusual; expectations from theory suggest that such alignments should occur less than 2% of the time. Searching for and characterizing such structures and accurately measuring their volume space density can therefore place constraints on the theory of structure formation. Such regions can also serve as laboratories for the study of galaxy formation in dense environments., Comment: 10 pages, 4 figures, submitted to Astrophysical Journal

We present an objective method to remove the stellar continuum emission from narrow-band images to derive emission-line images. The method is based on the skewness of the pixel histogram of the residual images. Specifically, we exploit a transition in the skewness of the signal in the continuum-subtracted image, which appears when the image changes from being under-subtracted to over-subtracted. Tests on one-dimensional artificial images demonstrate that the transition identifies the optimal scaling factor {\mu} to be used on the broad-band image IB in order to produce the optimal line-emission image IE, i.e., IE =IN - {\mu} IB, with IN the original (un-subtracted) narrow-band image. The advantage of this method is that it uses all information-bearing pixels in the final image, and not just a sub-set of those pixels (the latter being common in many traditional approaches to stellar continuum removal from narrow-band images). We apply our method to actual images, both from ground-based and space facilities, in particular to WFPC2 and ACS images from the Hubble Space Telescope, and we show that it is successful irrespective of the nature of the sources (point-like or extended). We also discuss the impact on the accuracy of the method of non-optimal images, such as those containing saturated sources or non-uniform background, and present `workarounds' for those problems., Comment: 22 pages, 19 figures, accepted for publication in PASP; The program and tutorial video about this method are available at https://sites.google.com/site/shongcontisub/home

Modern spectroscopic surveys produce large spectroscopic databases, generally with sizes well beyond the scope of manual investigation. The need arises, therefore, for an automated line detection method with objective indicators for detection significance. In this paper, we present an automated and objective method for emission line detection in spectroscopic surveys and apply this technique to 1574 spectra, obtained with the Hectospec spectrograph on the MMT Observatory (MMTO), to detect Lyman alpha emitters near z ~ 2.7. The basic idea is to generate on-source (signal plus noise) and off-source (noise only) mock observations using Monte Carlo simulations, and calculate completeness and reliability values, (C, R), for each simulated signal. By comparing the detections from real data with the Monte Carlo results, we assign the completeness and reliability values to each real detection. From 1574 spectra, we obtain 881 raw detections and, by removing low reliability detections, we finalize 649 detections from an automated pipeline. Most of high completeness and reliability detections, (C, R) ~ (1.0, 1.0), are robust detections when visually inspected; the low C and R detections are also marginal on visual inspection. This method at detecting faint sources is dependent on the accuracy of the sky subtraction., Comment: 19 pages, 16 figures, submitted to MNRAS

We present HST/WFC3 narrow-band imaging of the starburst galaxy M83 targeting the hydrogen recombination lines (H$\beta$, H$\alpha$ and Pa$\beta$), which we use to investigate the dust extinction in the HII regions. We derive extinction maps with 6 parsec spatial resolution from two combinations of hydrogen lines (H$\alpha$/H$\beta$ and H$\alpha$/Pa$\beta$), and show that the longer wavelengths probe larger optical depths, with $A_V$ values larger by $\gtrsim$1 mag than those derived from the shorter wavelengths. This difference leads to a factor $\gtrsim$2 discrepancy in the extinction-corrected H$\alpha$ luminosity, a significant effect when studying extragalactic HII regions. By comparing these observations to a series of simple models, we conclude that a large diversity of absorber/emitter geometric configurations can account for the data, implying a more complex physical structure than the classical foreground "dust screen" assumption. However, most data points are bracketed by the foreground screen and a model where dust and emitters are uniformly mixed. When averaged over large ($\gtrsim$100--200 pc) scales, the extinction becomes consistent with a "dust screen", suggesting that other geometries tend to be restricted to more local scales. Moreover, the extinction in any region can be described by a combination of the foreground screen and the uniform mixture model with weights of 1/3 and 2/3 in the center ($\lesssim$2 kpc), respectively, and 2/3 and 1/3 for the rest of the disk. This simple prescription significantly improves the accuracy of the dust extinction corrections and can be especially useful for pixel-based analyses of galaxies similar to M83., Comment: 7 pages, 4 figures, 1 table. Accepted for publication in ApJ Letters

(abridged) We investigate the properties of feedback-driven shocks in 8 nearby starburst galaxies using narrow-band imaging data from the Hubble Space Telescope (HST). We identify the shock--ionized component via the line diagnostic diagram \oiii/\hb vs. \sii (or \nii)/\ha, applied to resolved regions 3--15 pc in size. We divide our sample into three sub-samples: sub-solar (Holmberg II, NGC 1569, NGC 4214, NGC 4449, and NGC 5253), solar (He 2-10, NGC 3077) and super-solar (NGC 5236) for consistent shock measurements. For the sub-solar sub-sample, we derive three scaling relations: (1) $L_{shock} \propto {SFR}^{~0.62}$, (2) $L_{shock} \propto {\Sigma_{SFR,HL}}^{~0.92}$, and (3) $L_{shock}/L_{tot} \propto {(L_H/L_{\odot,H})}^{-0.65}$, where $L_{shock}$ is the \ha luminosity from shock--ionized gas, ${\Sigma_{SFR,HL}}$ the SFR per unit half-light area, $L_{tot}$ the total \ha luminosity, and $L_H/L_{\odot,H}$ the absolute H-band luminosity from 2MASS normalized to solar luminosity. The other two sub--samples do not have enough number statistics, but appear to follow the first scaling relation. The energy recovered indicates that the shocks from stellar feedback in our sample galaxies are fully radiative. If the scaling relations are applicable in general to stellar feedback, our results are similar to those by Hopkins et al. (2012) for galactic super winds. This similarity should, however, be taken with caution at this point, as the underlying physics that enables the transition from radiative shocks to gas outflows in galaxies is still poorly understood., Comment: 29 pages, 14 figures, accepted for publication in the ApJ

We investigate the ionization structure of the nebular gas in M83 using the line diagnostic diagram, [O III](5007 \degA)/H{\beta} vs. [S II](6716 \deg A+6731 \deg A)/H{\alpha} with the newly available narrowband images from the Wide Field Camera 3 (WFC3) of the Hubble Space Telescope (HST). We produce the diagnostic diagram on a pixel-by-pixel (0.2" x 0.2") basis and compare it with several photo- and shock-ionization models. For the photo-ionized gas, we observe a gradual increase of the log([O III]/H{\beta}) ratios from the center to the spiral arm, consistent with the metallicity gradient, as the H II regions go from super solar abundance to roughly solar abundance from the center out. Using the diagnostic diagram, we separate the photo-ionized from the shock-ionized component of the gas. We find that the shock-ionized H{\alpha} emission ranges from ~2% to about 15-33% of the total, depending on the separation criteria used. An interesting feature in the diagnostic diagram is an horizontal distribution around log([O III]/H{\beta}) ~ 0. This feature is well fit by a shock-ionization model with 2.0 Z\odot metallicity and shock velocities in the range of 250 km/s to 350 km/s. A low velocity shock component, < 200 km/s, is also detected, and is spatially located at the boundary between the outer ring and the spiral arm. The low velocity shock component can be due to : 1) supernova remnants located nearby, 2) dynamical interaction between the outer ring and the spiral arm, 3) abnormal line ratios from extreme local dust extinction. The current data do not enable us to distinguish among those three possible interpretations. Our main conclusion is that, even at the HST resolution, the shocked gas represents a small fraction of the total ionized gas emission at less than 33% of the total. However, it accounts for virtually all of the mechanical energy produced by the central starburst in M83., Comment: Accepted for publication in ApJ. aastex preprint 12pt, 21 pages, 13 figures

We use Spitzer 24 $\mu$m, 70 $\mu$m and ground based H$\alpha$ data for a sample of 40 SINGS galaxies to establish a star formation rate (SFR) indicator using 70 $\mu$m emission for sub--galactic ($\sim0.05-2\ \rm{kpc}$) line-emitting regions and to investigate limits in application. A linear correlation between 70 $\mu$m and SFR is found and a star formation indicator SFR(70) is proposed for line-emitting sub-galactic regions as $\rm \Sigma(SFR)\ ({M_{\odot}\cdot yr^{-1}\cdot kpc^{-2}})=9.4\times10^{-44}\ \Sigma(70)\ \rm{(ergs\cdot s^{-1}\cdot kpc^{-2})}$, for regions with $12+\rm{log(O/H)}\gtrsim8.4$ and $\rm \Sigma(SFR)\gtrsim10^{-3}\ (M_{\odot}\cdot yr^{-1}\cdot kpc^{-2})$, with a 1-$\sigma$ dispersion around the calibration of $\sim0.16$ dex. We also discuss the influence of metallicity on the scatter of the data. Comparing with the SFR indicator at 70 $\mu$m for integrated light from galaxies, we find that there is $\sim40%$ excess 70 $\mu$m emission in galaxies, which can be attributed to stellar populations not involved in the current star formation activity., Comment: 36 pages, 1 table, 18 figures, accepted by ApJ

The kinematics of damped Lyman alpha absorbers (DLAs) are difficult to reproduce in hierarchical galaxy formation models, particularly the preponderance of wide systems. We investigate DLA kinematics at z=3 using high-resolution cosmological hydrodynamical simulations that include a heuristic model for galactic outflows. Without outflows, our simulations fail to yield enough wide DLAs, as in previous studies. With outflows, predicted DLA kinematics are in much better agreement with observations. Comparing two outflow models, we find that a model based on momentum-driven wind scalings provides the best match to the observed DLA kinematic statistics of Prochaska & Wolfe. In this model, DLAs typically arise a few kpc away from galaxies that would be identified in emission. Narrow DLAs can arise from any halo and galaxy mass, but wide ones only arise in halos with mass >10^11 Mo, from either large central or small satellite galaxies. This implies that the success of this outflow model originates from being most efficient at pushing gas out from small satellite galaxies living in larger halos. This increases the cross-section for large halos relative to smaller ones, thereby yielding wider kinematics. Our simulations do not include radiative transfer effects or detailed metal tracking, and outflows are modeled heuristically, but they strongly suggest that galactic outflows are central to understanding DLA kinematics. An interesting consequence is that DLA kinematics may place constraints on the nature and efficiency of gas ejection from high-z galaxies., Comment: submitted to MNRAS