Looking for bright galaxies born in the early universe is fundamental to investigating the Epoch of Reionization, the era when the first stars and galaxies ionized the intergalactic medium. We utilize Hubble Space Telescope pure parallel imaging to select galaxy candidates at a time 500 to 650 million years after the Big Bang, which corresponds to redshifts z ~ 8-10. These data come from the Brightest of Reionizing Galaxies Survey (BoRG) Cycle 22 dataset, which consists of pure-parallel imaging in ~ 90 different lines of sight that sum up to an area of ~ 420 arcmin^2. This survey uses five filters and has the advantage (compared to the Cycle 21 BoRG program) of including imaging in the JH140 band, covering continuous wavelengths from the visible to near-infrared (lambda = 0.35um - 1.7um). This allows us to perform reliable selection of galaxies at z>8 using the photometric redshift technique. We use these galaxy candidates to constrain the bright end of the rest-frame ultraviolet luminosity function in this epoch. These candidates are excellent targets for follow-up observations, particularly with the James Webb Space Telescope.

We present a catalog of 208 $0.3 < z < 2.1$ Emission Line Galaxies (ELG) selected from 1D slitless spectroscopy obtained using Hubble's WFC3 G102 grism, as part of the Faint Infrared Grism Survey (FIGS). We identify ELG candidates by searching for significant peaks in all continuum-subtracted G102 spectra, and, where possible, confirm candidates by identifying consistent emission lines in other available spectra or with published spectroscopic redshifts. We provide derived emission line fluxes and errors, redshifts, and equivalent widths (EW) for H$\alpha$ $\lambda6563$, [OIII]$\lambda\lambda4959,5007$, and [OII]$\lambda\lambda3727$ emission lines, for emission line galaxies down to AB(F105W) $ > 28$ and $> 10^{-17}$ erg cm$^{-2}$ s$^{-1}$ line flux. We use the resulting line catalog to investigate a possible relationship between line emission and a galaxy's environment. We use 7th-nearest-neighbor distances to investigate the typical surroundings of ELGs compared to non-ELGs, and we find that [OIII] emitters are preferentially found at intermediate galaxy densities near galaxy groups. We characterize these ELGs in terms of the galaxy specific star formation rate (SSFR) versus stellar mass, and find no significant influence of environment on that relation. We calculate star formation rates (SFR), and find no dependence of SFR on local galaxy surface density for $0.3 < z < 0.8$ H$\alpha$ emitters and for $0.8<z<1.3$ [OIII] emitters. We find similar rates of close-pair interaction between ELGs and non-ELGs. For galaxy surface densities $\Sigma \leq 30$ Mpc$^{-2}$, we find no consistent effect of environment on star formation.

Studying the formation and evolution of galaxies at the earliest cosmic times, and their role in reionization, requires the deepest imaging possible. Ultra-deep surveys like the HUDF and HFF have pushed to mag \mAB$\,\sim\,$30, revealing galaxies at the faint end of the LF to $z$$\,\sim\,$9$\,-\,$11 and constraining their role in reionization. However, a key limitation of these fields is their size, only a few arcminutes (less than a Mpc at these redshifts), too small to probe large-scale environments or clustering properties of these galaxies, crucial for advancing our understanding of reionization. Achieving HUDF-quality depth over areas $\sim$100 times larger becomes possible with a mission like the Wide Field Infrared Survey Telescope (WFIRST), a 2.4-m telescope with similar optical properties to HST, with a field of view of $\sim$1000 arcmin$^2$, $\sim$100$\times$ the area of the HST/ACS HUDF. This whitepaper motivates an Ultra-Deep Field survey with WFIRST, covering $\sim$100$\,-\,$300$\times$ the area of the HUDF, or up to $\sim$1 deg$^2$, to \mAB$\,\sim\,$30, potentially revealing thousands of galaxies and AGN at the faint end of the LF, at or beyond $z$\,$\sim$\,9$\,-\,$10 in the epoch of reionization, and tracing their LSS environments, dramatically increasing the discovery potential at these redshifts. (Note: This paper is a somewhat expanded version of one that was submitted as input to the Astro2020 Decadal Survey, with this version including an Appendix (which exceeded the Astro2020 page limits), describing how the science drivers for a WFIRST Ultra Deep Field might map into a notional observing program, including the filters used and exposure times needed to achieve these depths.)

Recent near-IR power-spectra and panchromatic Extragalactic Background Light measurements provide upper limits on the near-IR surface brightness (SB>31 mag/arcsec^2) that may come from Pop III stars and accretion disks around resulting stellar-mass black holes (BHs) in the epoch of First Light (z=7-17). Physical parameters for zero metallicity Pop III stars at z>7 can be estimated from MESA stellar evolution models through helium-depletion, and for BH accretion disks from quasar microlensing results and multicolor accretion models. Second-generation stars can form at higher multiplicity, so that BH accretion disks may be fed by Roche-lobe overflow from lower-mass companions in their AGB stage. The near-IR SB constraints can be used to calculate the number of caustic transits behind lensing clusters that JWST and the 25~39 m ground-based telescopes may detect for both Pop III stars and stellar mass BH accretion disks. Because Pop III stars and stellar mass BH accretion disks have sizes of a few x 10^-11 arcsec at z>7, typical caustic magnifications can be mu=10^4~10^5, with rise times of hours and decline times of < 1 year for cluster transverse velocities of v<1000 km/s. Microlensing by intracluster medium objects can modify transit magnifications, and lengthen visibility times. Depending on BH masses, accretion-disk radii and feeding efficiencies, stellar-mass BH accretion-disk caustic transits could outnumber those from Pop III stars. To observe Pop III caustic transits directly may require monitoring 3~30 lensing clusters to AB< 29 mag over a decade or more. Such a program must be started with JWST in Cycle 1, and -- depending on the role of microlensing in the Intra Cluster Light -- should be continued for decades with the GMT and TMT, where JWST and the ground-based telescopes each will play a unique and strongly complementary role.

The epoch of reionization (6 < z < 10) marks the period in our universe when the first large galaxies grew to fruition, and began to affect the universe around them. Massive stars, and potentially accreting supermassive black holes, filled the universe with ionizing radiation, burning off the haze of neutral gas that had filled the intergalactic medium (IGM) since recombination (z~1000). The evolution of this process constrains key properties of these earliest luminous sources, thus observationally constraining reionization is a key science goal for the next decade. The measurement of Lyman-alpha emission from photometrically-identified galaxies is a highly constraining probe of reionization, as a neutral IGM will resonantly scatter these photons, reducing detectability. While significant work has been done with 8-10m telescopes, these observations require extremely large telescopes (ELTs); the flux limits available from today's 10m class telescopes are sufficient for only the brightest known galaxies (m < 26). Ultra-deep surveys with the Giant Magellan Telescope (GMT) and Thirty Meter Telescope (TMT) will be capable of detecting Lyman-alpha emission from galaxies 2-3 magnitudes fainter than today's deepest surveys. Wide-field fiber spectroscopy on the GMT combined with narrow-field AO-assisted slit spectroscopy on the TMT will be able to probe the expected size of ionized bubbles throughout the epoch of reionization, following up degree scale deep imaging surveys with the Wide Field Infrared Space Telescope. These data will provide the first resolved Lyman-alpha-based maps of the ionized intergalactic medium throughout the epoch of reionization, constraining models of both the temporal and spatial evolution of this phase change.

Realizing the utility of Ly$\alpha$ emission to trace the evolution of the intergalactic medium (IGM) during the epoch of reionization requires deep spectroscopy across the boundary of optical and near-infrared (NIR) spectrographs at $z\sim7.2$ when Ly$\alpha$ emission is at $\sim$1$\mu$m. Our Texas Spectroscopic Search for Ly$\alpha$ Emission at the End of Reionization includes 18 nights of deep spectroscopic observations using the Keck DEIMOS (optical) and MOSFIRE (NIR) spectrographs. Within this dataset we observe Ly$\alpha$ emission from 183 photometric-redshift selected galaxies at $z =$ 5.5 - 8.3 from the Cosmic Assembly Near infrared Deep Extragalactic Legacy Survey (CANDELS). Our overlapping MOSFIRE observations, over 84 galaxies total, provide the deepest NIR spectroscopic data yet obtained for Ly$\alpha$ from galaxies $z > 7$, with $>16$ hr integration time for four observed galaxies. Here we analyze these four targets, and we report the discovery of a new $z = 7.60$ Ly$\alpha$ detection as well as provide an updated observation of the previously confirmed $z=7.51$ Ly$\alpha$ emission from Finkelstein et al. (2013) with a $\sim$3$\times$ longer exposure time. Our analysis of these Ly$\alpha$ emission line profiles reveal a significant asymmetric shape. The two detected Ly$\alpha$ emission lines from bright sources ($M_{\text{UV}}<-20.25$) could imply that these bright galaxies inhabit ionized bubbles in a partially neutral IGM, although deeper exposures may yet reveal Ly$\alpha$ emission in the fainter sources.

We present the results from the application of a two-dimensional emission line detection method, EMission-line two-Dimensional (EM2D), to the near-infrared G102 grism observations obtained with the Wide-Field Camera 3 (WFC3) as part of the Cycle 22 {\em Hubble Space Telescope} Treasury Program: the Faint Infrared Grism Survey (FIGS). Using the EM2D method, we have assembled a catalog of emission line galaxies (ELGs) with resolved star formation from each of the four FIGS fields. Not only can one better assess the global properties of ELGs, but the EM2D method allows for the analysis and an improved study of the individual emission-line region {\it within} each galaxy. This paper includes a description of the methodology, advantages, and the first results of the EM2D method applied to ELGs in FIGS. The advantage of 2D emission line measurements includes significant improvement of galaxy redshift measurements, approaching the level of accuracy seen in high-spectral-resolution data, but with greater efficiency; and the ability to identify and measure the properties of multiple sites of star-formation and over scales of $\sim$ 1 kpc within individual galaxies out to z $\sim$ 4. The EM2D method also significantly improves the reliability of high-redshift ($z\sim7$) Lyman-$\alpha$ detections. Coupled with the wide field of view and high efficiency of space-based grism observations, EM2D provides a noteworthy improvement on the physical parameters that can be extracted from grism observations.

EUSO-Balloon is a pathfinder mission for the Extreme Universe Space Observatory onboard the Japanese Experiment Module (JEM-EUSO). It was launched on the moonless night of the 25$^{th}$ of August 2014 from Timmins, Canada. The flight ended successfully after maintaining the target altitude of 38 km for five hours. One part of the mission was a 2.5 hour underflight using a helicopter equipped with three UV light sources (LED, xenon flasher and laser) to perform an inflight calibration and examine the detectors capability to measure tracks moving at the speed of light. We describe the helicopter laser system and details of the underflight as well as how the laser tracks were recorded and found in the data. These are the first recorded laser tracks measured from a fluorescence detector looking down on the atmosphere. Finally, we present a first reconstruction of the direction of the laser tracks relative to the detector.

We present full spectral scans from 200-670$\mu$m of 26 Class 0+I protostellar sources, obtained with $Herschel$-SPIRE, as part of the "COPS-SPIRE" Open Time program, complementary to the DIGIT and WISH Key programs. Based on our nearly continuous, line-free spectra from 200-670 $\mu$m, the calculated bolometric luminosities ($L_{\rm bol}$) increase by 50% on average, and the bolometric temperatures ($T_{\rm bol}$) decrease by 10% on average, in comparison with the measurements without Herschel. Fifteen protostars have the same Class using $T_{\rm bol}$ and $L_{\rm bol}$/$L_{\rm submm}$. We identify rotational transitions of CO lines from J=4-3 to J=13-12, along with emission lines of $^{13}$CO, HCO$^+$, H$_{2}$O, and [CI]. The ratios of $^{12}$CO to $^{13}$CO indicate that $^{12}$CO emission remains optically thick for $J_{\rm up}$ < 13. We fit up to four components of temperature from the rotational diagram with flexible break points to separate the components. The distribution of rotational temperatures shows a primary population around 100 K with a secondary population at $\sim$350 K. We quantify the correlations of each line pair found in our dataset, and find the strength of correlation of CO lines decreases as the difference between $J$-level between two CO lines increases. The multiple origins of CO emission previously revealed by velocity-resolved profiles are consistent with this smooth distribution if each physical component contributes to a wide range of CO lines with significant overlap in the CO ladder. We investigate the spatial extent of CO emission and find that the morphology is more centrally peaked and less bipolar at high-$J$ lines. We find the CO emission observed with SPIRE related to outflows, which consists two components, the entrained gas and shocked gas, as revealed by our rotational diagram analysis as well as the studies with velocity-resolved CO emission.

The distribution of Ly$\alpha$ emission is an presently accessible method for studying the state of the intergalactic medium (IGM) into the reionization era. We carried out deep spectroscopic observations in order to search for Ly$\alpha$ emission from galaxies with photometric redshifts $z$ = 5.5 - 8.3 selected from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS). Utilizing data from the Keck/DEIMOS spectrograph, we explore a wavelength coverage of Ly$\alpha$ emission at $z$ ~ 5 - 7 with four nights of spectroscopic observations for 118 galaxies, detecting five emission lines with ~ 5$\sigma$ significance: three in the GOODS-N and two in the GOODS-S field. We constrain the equivalent width (EW) distribution of Ly$\alpha$ emission by comparing the number of detected objects with the expected number constructed from detailed simulations of mock emission lines that account for the observational conditions (e.g., exposure time, wavelength coverage, and sky emission) and galaxy photometric redshift probability distribution functions. The Ly$\alpha$ EW distribution is well described by an exponential form, $\text{dN/dEW}\propto \text{exp(-EW/}W_0)$, characterized by the $e$-folding scale ($W_0$) of ~ 60 - 100$\AA$ at 0.3 < $z$ < 6. By contrast, our measure of the Ly$\alpha$ EW distribution at 6.0 < $z$ < 7.0 rejects a Ly$\alpha$ EW distribution with $W_0$ > 36.4$\AA$ (125.3$\AA$) at 1$\sigma$ (2$\sigma$) significance. This provides additional evidence that the EW distribution of Ly$\alpha$ declines at $z$ > 6, suggesting an increasing fraction of neutral hydrogen in the IGM at that epoch.

We present the results of an unbiased search for Ly{\alpha} emission from continuum-selected 6 < z < 8 galaxies. Our dataset consists of 160 orbits of G102 slitless grism spectroscopy obtained with the Hubble Space Telescope (HST) Wide Field Camera 3 as part of the Faint Infrared Grism Survey (FIGS; PI: Malhotra), which obtains deep slitless spectra of all sources in four fields, and was designed to minimize contamination in observations of previously-identified high-redshift galaxy candidates. The FIGS data can potentially spectroscopically confirm the redshifts of galaxies, and as Ly{\alpha} emission is resonantly scattered by neutral gas, FIGS can also constrain the ionization state of the intergalactic medium (IGM) during the epoch of reionization. These data have sufficient depth to detect Ly{\alpha} emission in this epoch, as Tilvi et al. (2016) have published the FIGS detection of previously known (Finkelstein et al. 2013) Ly{\alpha} emission at z = 7.51. The FIGS data use five separate roll-angles of HST to mitigate the contamination by nearby galaxies. We created a method that accounts for and removes the contamination from surrounding galaxies, and also removes any dispersed continuum light from each individual spectrum (Pirzkal et al. 2017). We searched for significant (> 4{\sigma}) emission lines using two different automated detection methods, free of any visual inspection biases. Applying these methods on photometrically-selected high-redshift candidates between 6 < z < 8 we find two emission lines, one previously published by Tilvi et al. (2016), and a new line at 1.028{\mu}m. We identify this lines as Ly{\alpha} at z = 7.452 +/- 0.003. This newly spectroscopically confirmed galaxy has the highest Ly{\alpha} rest-frame equivalent width (EW Ly{\alpha}) yet published at z > 7 (140.3 +/- 19.0{\AA}).

We study the nature of feedback mechanisms in the 11 CLASH brightest cluster galaxies (BCGs) that exhibit extended ultraviolet and nebular line emission features. We estimate star formation rates (SFRs), dust masses, and starburst durations using a Bayesian photometry-fitting technique that accounts for both stellar and dust emission from the UV through far IR. By comparing these quantities to intracluster medium (ICM) cooling times and free-fall times derived from X-ray observations and lensing estimates of the cluster mass distribution, we discover a tight relationship between the BCG SFR and the ICM cooling time to free-fall time ratio, $t_{cool}/t_{ff}$, with an upper limit on the intrinsic scatter of 0.15 dex. Furthermore, starburst durations may correlate with ICM cooling times at a radius of $0.025R_{500}$, and the two quantities converge upon reaching the Gyr regime. Our results provide a direct observational link between the thermodynamical state of the ICM and the intensity and duration of BCG star formation activity, and appear consistent with a scenario where active galactic nuclei (AGN) induce condensation of thermally unstable ICM overdensities that fuel long-duration ($>$ 1 Gyr) BCG starbursts. This scenario can explain (a) how gas with a low cooling time is depleted without causing a cooling flow and (b) the scaling relationship between SFR and $t_{cool}/t_{ff}$. We also find that the scaling relation between SFR and dust mass in BCGs with SFRs $<100$ M$_{\odot}$ yr$^{-1}$ is similar to star-forming field galaxies; BCGs with large ($>100$ M$_{\odot}$ yr$^{-1}$) SFRs have dust masses comparable to extreme starbursts.

The Faint Infrared Grism Survey (FIGS) is a deep Hubble Space Telescope (HST) WFC3/IR (Wide Field Camera 3 Infrared) slitless spectroscopic survey of four deep fields. Two fields are located in the Great Observatories Origins Deep Survey-North (GOODS-N) area and two fields are located in the Great Observatories Origins Deep Survey-South (GOODS-S) area. One of the southern fields selected is the Hubble Ultra Deep Field. Each of these four fields were observed using the WFC3/G102 grism (0.8$\mu m$-1.15$\mu m$ continuous coverage) with a total exposure time of 40 orbits (~ 100 kilo-seconds) per field. This reaches a 3 sigma continuum depth of ~26 AB magnitudes and probes emission lines to $\approx 10^{-17}\ erg\ s^{-1} \ cm^{-2}$. This paper details the four FIGS fields and the overall observational strategy of the project. A detailed description of the Simulation Based Extraction (SBE) method used to extract and combine over 10000 spectra of over 2000 distinct sources brighter than m_F105W=26.5 mag is provided. High fidelity simulations of the observations is shown to significantly improve the background subtraction process, the spectral contamination estimates, and the final flux calibration. This allows for the combination of multiple spectra to produce a final high quality, deep, 1D-spectra for each object in the survey.

We present the COPS-DIGIT-FOOSH (CDF) Herschel spectroscopy data product archive, and related ancillary data products, along with data fidelity assessments, and a user-created archive in collaboration with the Herschel-PACS and SPIRE ICC groups. Our products include datacubes, contour maps, automated line fitting results, and best 1-D spectra products for all protostellar and disk sources observed with PACS in RangeScan mode for two observing programs: the DIGIT Open Time Key Program (KPOT_nevans_1 and SDP_nevans_1; PI: N. Evans), and the FOOSH Open Time Program (OT1_jgreen02_2; PI: J. Green). In addition, we provide our best SPIRE-FTS spectroscopic products for the COPS Open Time Program (OT2_jgreen02_6; PI: J. Green) and FOOSH sources. We include details of data processing, descriptions of output products, and tests of their reliability for user applications. We identify the parts of the dataset to be used with caution. The resulting absolute flux calibration has improved in almost all cases. Compared to previous reductions, the resulting rotational temperatures and numbers of CO molecules have changed substantially in some sources. On average, however, the rotational temperatures have not changed substantially (< 2%), but the number of warm (Trot ~ 300 K) CO molecules has increased by about 18%.

We utilize observations of sub-millimeter rotational transitions of CO from a Herschel Cycle 2 open time program ("COPS", PI: J. Green) to identify previously predicted turbulent dissipation by magnetohydrodynamic (MHD) shocks in molecular clouds. We find evidence of the shocks expected for dissipation of MHD turbulence in material not associated with any protostar. Two models fit about equally well: model 1 has a density of 10$^{3}$ cm$^{-3}$, a shock velocity of 3 km s$^{-1}$, and a magnetic field strength of 4 ${\mu}$G; model 2 has a density of 10$^{3.5}$ cm$^{-3}$, a shock velocity of $2$ km s$^{-1}$, and a magnetic field strength of 8 $\mu$G. Timescales for decay of turbulence in this region are comparable to crossing times. Transitions of CO up to $J$ of 8, observed close to active sites of star formation, but not within outflows, can trace turbulent dissipation of shocks stirred by formation processes. Although the transitions are difficult to detect at individual positions, our Herschel-SPIRE survey of protostars provides a grid of spatially-distributed spectra within molecular clouds. We averaged all spatial positions away from known outflows near seven protostars. We find significant agreement with predictions of models of turbulent dissipation in slightly denser (10$^{3.5}$ cm$^{-3}$) material with a stronger magnetic field (24 $\mu$G) than in the general molecular cloud.

This brief historical review highlights the early work of Hayashi and his associates on the thermal physics of star-forming clouds, as summarized in the temperature-density diagrams first presented by this group. Some of the more recent developments in this subject, including its application to understanding stellar masses and to understanding the formation of the first stars, are also briefly reviewed.

This contribution contains the introductory historical remarks that I presented at IAU Symposium 270 on "Computational Star Formation" held in Barcelona, Spain, May 31 - June 4, 2010. I give a personal view of some of the early histoy of the subject and I comment on some of the important things that were learned from numerical work on star formation.

The formation of compact objects like stars and black holes is strongly constrained by the requirement that nearly all of the initial angular momentum of the diffuse material from which they form must be removed or redistributed during the formation process. The mechanisms that may be involved and their implications are discussed for (1) low-mass stars, most of which probably form in binary or multiple systems; (2) massive stars, which typically form in clusters; and (3) supermassive black holes that form in galactic nuclei. It is suggested that in all cases, gravitational interactions with other stars or mass concentrations in a forming system play an important role in redistributing angular momentum and thereby enabling the formation of a compact object. If this is true, the formation of stars and black holes must be a more complex, dynamic, and chaotic process than in standard models. The gravitational interactions that redistribute angular momentum tend to couple the mass of a forming object to the mass of the system, and this may have important implications for mass ratios in binaries, the upper stellar IMF in clusters, and the masses of supermassive black holes in galaxies.

The stellar mass in disk galaxies scales approximately with the fourth power of the rotation velocity, and the masses of the central black holes in galactic nuclei scale approximately with the fourth power of the bulge velocity dispersion. It is shown here that these relations can be accounted for if, in a forming galaxy with an isothermal mass distribution, gas with a column density above about 8 Msun/pc^2 goes into stars while gas with a column density above about 2 g/cm^2 (10^4 Msun/pc^2) goes into a central black hole. The lower critical value is close to the column density of about 10 Msun/pc^2 at which atomic gas becomes molecular, and the upper value agrees approximately with the column density of about 1 g/cm^2 at which the gas becomes optically thick to its cooling radiation. These results are plausible because molecule formation is evidently necessary for star formation, and because the onset of a high optical depth in a galactic nucleus may suppress continuing star formation and favour the growth of a central black hole.

We study gravitational instabilities in disks, with special attention to the most massive clumps that form because they are expected to be the progenitors of globular-type clusters. The maximum unstable mass is set by rotation and depends only on the surface density and orbital frequency of the disk. We propose that the formation of massive clusters is related to this largest scale in galaxies not stabilized by rotation. Using data from the literature, we predict that globular-like clusters can form in nuclear starburst disks and protogalactic disks but not in typical spiral galaxies, in agreement with observations.

We study gravitational instabilities in disks, with special attention to the most massive clumps that form because they are expected to be the progenitors of globular-type clusters. The maximum unstable mass is set by rotation and depends only on the surface density and orbital frequency of the disk. We propose that the formation of massive clusters is related to this largest scale in galaxies not stabilized by rotation. Using data from the literature, we predict that globular-like clusters can form in nuclear starburst disks and protogalactic disks but not in typical spiral galaxies, in agreement with observations.

Although the basic physics of star formation is classical, numerical simulations have yielded essential insights into how stars form. They show that star formation is a highly nonuniform runaway process characterized by the emergence of nearly singular peaks in density, followed by the accretional growth of embryo stars that form at these density peaks. Circumstellar disks often form from the gas being accreted by the forming stars, and accretion from these disks may be episodic, driven by gravitational instabilities or by protostellar interactions. Star-forming clouds typically develop filamentary structures, which may, along with the thermal physics, play an important role in the origin of stellar masses because of the sensitivity of filament fragmentation to temperature variations. Simulations of the formation of star clusters show that the most massive stars form by continuing accretion in the dense cluster cores, and this again is a runaway process that couples star formation and cluster formation. Star-forming clouds also tend to develop hierarchical structures, and smaller groups of forming objects tend to merge into progressively larger ones, a generic feature of self-gravitating systems that is common to star formation and galaxy formation. Because of the large range of scales and the complex dynamics involved, analytic models cannot adequately describe many processes of star formation, and detailed numerical simulations are needed to advance our understanding of the subject.

The essential features of the stellar Initial Mass Function are, rather generally, (1) a peak at a mass of a few tenths of a solar mass, and (2) a power-law tail toward higher masses that is similar to the original Salpeter function. Recent work suggests that the IMF peak reflects a preferred scale of fragmentation associated with the transition from a cooling phase of collapse at low densities to a nearly isothermal phase at higher densities, where the gas becomes thermally coupled to the dust. The Salpeter power law is plausibly produced, at least in part, by scale-free accretion processes that build up massive stars in dense environments. The young stars at the Galactic Center appear to have unusually high masses, possibly because of a high minimum mass resulting from the high opacity of the dense star-forming gas.

We review recent advances in our understanding of the origin of the initial mass function (IMF). We emphasize the use of numerical simulations to investigate how each physical process involved in star formation affects the resulting IMF. We stress that it is insufficient to just reproduce the IMF, but that any successful model needs to account for the many observed properties of star forming regions including clustering, mass segregation and binarity. Fragmentation involving the interplay of gravity, turbulence, and thermal effects is probably responsible for setting the characteristic stellar mass. Low-mass stars and brown dwarfs can form through the fragmentation of dense filaments and disks, possibly followed by early ejection from these dense environments which truncates their growth in mass. Higher-mass stars and the Salpeter-like slope of the IMF are most likely formed through continued accretion in a clustered environment. The effects of feedback and magnetic fields on the origin of the IMF are still largely unclear. Lastly, we discuss a number of outstanding problems that need to be addressed in order to develop a complete theory for the origin of the IMF.

The thermal properties of star-forming clouds have an important influence on how they fragment into stars, and it is suggested in this paper that the low-mass stellar IMF, which appears to be almost universal, is determined largely by the thermal physics of these clouds. In particular, it is suggested that the characteristic stellar mass, a little below one solar mass, is determined by the transition from an initial cooling phase of collapse to a later phase of slowly rising temperature that occurs when the gas becomes thermally coupled to the dust. Numerical simulations support the hypothesis that the Jeans mass at this transition point plays an important role in determining the peak mass of the IMF. A filamentary geometry may also play a key role in the fragmentation process because the isothermal case is a critical one for the collapse of a cylinder: the collapse and fragmentation of a cylinder can continue freely as long as the temperature continues to decrease, but not if it begins to increase. The limited available results on the dependence of the thermal properties of clouds on metallicity do not suggest a strong dependence of the IMF on metallicity, but the far-infrared background radiation in starburst regions and in the early universe may significantly shift the peak mass to higher masses in these situations.

Identifying the processes that determine the initial mass function of stars (IMF) is a fundamental problem in star formation theory. One of the major uncertainties is the exact chemical state of the star forming gas and its influence on the dynamical evolution. Most simulations of star forming clusters use an isothermal equation of state (EOS). However, theoretical predictions and observations suggest that the effective polytropic exponent gamma in the EOS varies with density. We address these issues and study the effect of a piecewise polytropic EOS on the formation of stellar clusters in turbulent, self-gravitating molecular clouds using three-dimensional, smoothed particle hydrodynamics simulations. To approximate the results of published predictions of the thermal behavior of collapsing clouds, we increase the polytropic exponent gamma from 0.7 to 1.1 at some chosen density n_c, which we vary. The change of thermodynamic state at n_c selects a characteristic mass scale for fragmentation M_ch, which we relate to the peak of the observed IMF. Our investigation generally supports the idea that the distribution of stellar masses depends mainly on the thermodynamic state of the star-forming gas. The thermodynamic state of interstellar gas is a result of the balance between heating and cooling processes, which in turn are determined by fundamental atomic and molecular physics and by chemical abundances. Given the abundances, the derivation of a characteristic stellar mass can thus be based on universal quantities and constants.

(abridged version) Identifying the processes that determine the initial mass function of stars (IMF) is a fundamental problem in star formation theory. One of the major uncertainties is the exact chemical state of the star forming gas and its influence on the dynamical evolution. Most simulations of star forming clusters use an isothermal equation of state (EOS). We address these issues and study the effect of a piecewise polytropic EOS on the formation of stellar clusters in turbulent, self-gravitating molecular clouds. We increase the polytropic exponent gamma from 0.7 to 1.1 at some chosen density n_c, which we vary from from 4.3x10^4 cm^-3 to 4.3x10^7 cm^-3. The change of thermodynamic state at n_c selects a characteristic mass scale for fragmentation M_ch, which we relate to the peak of the observed IMF. We find a relation M_ch ~ n_c^-0.5, supporting the idea that the distribution of stellar masses largely depends on the thermodynamic state of the star-forming gas.

It is suggested that the thermal physics of star-forming clouds may play a more important role than has usually been recognized in the origin of the stellar IMF and in determining a characteristic mass scale. The importance of the thermal physics has been clearly demonstrated for the formation of the first stars in the universe, where it is well understood and results in cooling to a characteristic minimum temperature at a preferred density, and hence in a characteristic scale for fragmentation. In present-day star-forming clouds, an analogous situation may exist in that at low densities the temperature is expected to decrease with increasing density, reaching a minimum when the gas becomes thermally coupled to the dust and then rising slowly at higher densities. A minimum temperature of about 5 K is predicted to occur at a density of the order of 10^(-18) g cm^(-3), and at this point the Jeans mass is about 0.3 solar masses, similar to the mass at which the IMF peaks. If most of the fragmentation in star-forming clouds occurs in filaments, as is suggested by many simulations as well as by observations, fragmentation seems likely to occur preferentially at the density where the temperature reaches a minimum, and the Jeans mass at this point may then determine a characteristic scale for fragmentation and hence a preferred stellar mass.

Using SPH numerical simulations, we investigate the effects of gas on the inspiral and merger of a massive black hole binary. This study is motivated by the very massive nuclear gas disks observed in the central regions of merging galaxies. Here we present results that expand on the treatment in a previous work (Escala, Larson, Coppi & Mardones 2004; henceforth Paper I), by studying more realistic models. We run a variety of models, ranging from simulations with a relatively smooth gas disk to cases in which the gas has a more clumpy spatial distribution. We also vary the inclination angle between the plane of the binary and the plane of the disk, and the mass ratio between the MBHs and the gaseous disk. We find that as in Paper I, in the early evolution of the system the binary separation diminishes due to gravitational drag, and in the later stages the medium responds by forming an ellipsoidal density enhancement whose axis lags behind the binary axis; this offset produces a torque on the binary that causes continuing loss of angular momentum and is able to reduce the separation to distances where gravitational radiation is efficient. In the simulations that we perform, we find that the coalescence timescale for the MBH binary varies between 5 10^6 yr and 2.5 10^7 yr. For MBHs that satisfy the observed `m-sigma_c' relation, we predict that in a merger of galaxies that have at least 1% of their total mass in gas, the MBHs will coalesce soon after the galaxies merge. We also predict that if the MBHs depart from the `m-sigma_c' relation, the MBH binary can create a circumbinary gap in the disk that stalls the coalescence, but this gap formation can act as a self-regulatory mechanism on MBH growth that can help to explain the existence of the `m-sigma_c' relation.

Using high-resolution SPH numerical simulations, we investigate the effects of gas on the inspiral and merger of a massive black hole binary. This study is motivated by both observational and theoretical work that indicate the presence of large amounts of gas in the central regions of merging galaxies. N-body simulations have shown that the coalescence of a massive black hole binary eventually stalls in a stellar background. However, our simulations suggest that the massive black hole binary will finally merge if it is embedded in a gaseous background. Here we present results in which the gas is assumed to be initially spherical with a relatively smooth distribution. In the early evolution of the binary, the separation dimishes due to the gravitational drag exerted by the background gas. In the later stages, when the binary dominates the gravitational potential in its vicinity, the medium responds by forming an ellipsoidal density enhancement whose axis lags behind the binary axis, and this offset produces a torque on the binary that causes continuing loss of angular momentum and is able to reduce the binary separation to distances where gravitational radiation is efficient. Assuming typical parameters from observations of Ultra Luminous Infrared Galaxies, we predict that a black hole binary will merge within $10^{7}$yrs; therefore these results imply that in a merger of gas-rich galaxies, any massive central black holes will coalescence soon after the galaxies merge. Our work thus supports scenarios of massive black hole evolution and growth where hierarchical merging plays an important role. The final coalescence of the black holes leads to gravitational radiation emission that would be detectable up to high redshift by LISA. We show that similar physical effects are important for the formation of close binary stars.

We review recent theoretical results on the formation of the first stars in the universe, and emphasize related open questions. In particular, we discuss the initial conditions for Population III star formation, as given by variants of the cold dark matter cosmology. Numerical simulations have investigated the collapse and the fragmentation of metal-free gas, showing that the first stars were predominantly very massive. The exact determination of the stellar masses, and the precise form of the primordial initial mass function, is still hampered by our limited understanding of the accretion physics and the protostellar feedback effects. We address the importance of heavy elements in bringing about the transition from an early star formation mode dominated by massive stars, to the familiar mode dominated by low mass stars, at later times. We show how complementary observations, both at high redshifts and in our local cosmic neighborhood, can be utilized to probe the first epoch of star formation.

Our current understanding of the physical processes of star formation is reviewed, with emphasis on processes occurring in molecular clouds like those observed nearby. The dense cores of these clouds are predicted to undergo gravitational collapse characterized by the runaway growth of a central density peak that evolves toward a singularity. As long as collapse can occur, rotation and magnetic fields do not change this qualitative behavior. The result is that a very small embryonic star or protostar forms and grows by accretion at a rate that is initially high but declines with time as the surrounding envelope is depleted. Rotation causes some of the remaining matter to form a disk around the protostar, but accretion from protostellar disks is not well understood and may be variable. Most, and possibly all, stars form in binary or multiple systems in which gravitational interactions can play a role in redistributing angular momentum and driving episodes of disk accretion. Variable accretion may account for some peculiarities of young stars such as flareups and jet production, and protostellar interactions in forming systems of stars will also have important implications for planet formation. The most massive stars form in the densest environments by processes that are not yet well understood but may include violent interactions and mergers. The formation of the most massive stars may have similarities to the formation and growth of massive black holes in very dense environments.

A brief review is given of the basic observed features of the stellar IMF, and also of some of the theoretical ideas and simulations that may help to explain these features. The characteristic stellar mass of order one solar mass may derive from the typical masses of the observed star-forming clumps in molecular clouds, and the typical clump mass may in turn be determined mainly by the Jeans mass, as is true in many numerical simulations of cloud collapse and fragmentation. The Salpeter power law that approximates the IMF at higher masses is less well understood, but it may result from the effects of continuing gas accretion by the more massive forming stars together with the effects of interactions and mergers in dense forming stellar systems. The mass of the most massive star that forms in a cluster is observed to increase systematically with the mass of the cluster, and the most massive clusters may form stars that are massive enough to collapse completely to black holes.

Nearly all of the initial angular momentum of the matter that goes into each forming star must somehow be removed or redistributed during the formation process. The possible transport mechanisms and the possible fates of the excess angular momentum are discussed, and it is argued that transport processes in disks are probably not sufficient by themselves to solve the angular momentum problem, while tidal interactions with other stars in forming binary or multiple systems are likely to be of very general importance in redistributing angular momentum during the star formation process. Most, if not all, stars probably form in binary or multiple systems, and tidal torques in these systems can transfer much of the angular momentum from the gas around each forming star to the orbital motions of the companion stars. Tidally generated waves in circumstellar disks may contribute to the overall redistribution of angular momentum. Stars may gain much of their mass by tidally triggered bursts of rapid accretion, and these bursts could account for some of the most energetic phenomena of the earliest stages of stellar evolution, such as jet-like outflows. If tidal interactions are indeed of general importance, planet-forming disks may often have a more chaotic and violent early evolution than in standard models, and shock heating events may be common. Interactions in a hierarchy of subgroups may play a role in building up massive stars in clusters and in determining the form of the upper IMF. Many of the processes discussed here have analogs on galactic scales, and there may be similarities between the formation of massive stars by interaction-driven accretion processes in clusters and the buildup of massive black holes in galactic nuclei.

Paramount among the processes that ended the cosmic `dark ages' must have been the formation of the first generation of stars (the so-called Population III). We summarize recent progress in constraining its nature, and we discuss the basic physical reasons for why the first stars might have been quite massive, possibly even very massive with M > 100 M_sun.

Recent theoretical investigations have suggested that the formation of the very first stars, forming out of metal-free gas, was fundamentally different from the present-day case. In this paper, we study the effect of metallicity on the evolution of the gas in a collapsing dark matter mini-halo. We model such a system as an isolated 3\sigma peak of mass 2x10^6 M_sun that collapses at z_coll=30, using smoothed particle hydrodynamics. The gas has a supposed level of pre-enrichment of either 10^-4 Z_sun or 10^-3 Z_sun. We find that the evolution proceeds very differently for the two cases. The gas in the lower metallicity simulation fails to undergo continued collapse and fragmentation, whereas the gas in the higher metallicity case dissipatively settles into the center of the dark matter halo. The central gas, characterized by densities n > 10^4 cm^-3, and a temperature, T \sim 90 K, which closely follows that of the CMB, is gravitationally unstable and undergoes vigorous fragmentation. We discuss the physical reason for the existence of a critical metallicity, Z_crit \sim 5x10^-4 Z_sun, and its possible dependence on redshift. Compared to the pure H/He case, the fragmentation of the 10^-3 Z_sun gas leads to a larger relative number of low-mass clumps.

To constrain the nature of the very first stars, we investigate the collapse and fragmentation of primordial, metal-free gas clouds. We explore the physics of primordial star formation by means of three-dimensional simulations of the dark matter and gas components, using smoothed particle hydrodynamics, under a wide range of initial conditions, including the initial spin, the total mass of the halo, the redshift of virialization, the power spectrum of the DM fluctuations, the presence of HD cooling, and the number of particles employed in the simulation. We find characteristic values for the temperature, T ~ a few 100 K, and the density, n ~ 10^3-10^4 cm^-3, characterising the gas at the end of the initial free-fall phase. These values are rather insensitive to the initial conditions. The corresponding Jeans mass is M_J ~ 10^3 M_sun. The existence of these characteristic values has a robust explanation in the microphysics of H2 cooling, connected to the minimum temperature that can be reached with the H2 coolant, and to the critical density at which the transition takes place betweeb levels being populated according to NLTE, and according to LTE. In all cases, the gas dissipatively settles into an irregular, central configuration which has a filamentary and knotty appearance. The fluid regions with the highest densities are the first to undergo runaway collapse due to gravitational instability, and to form clumps with initial masses ~ 10^3 M_sun, close to the characteristic Jeans scale. These results suggest that the first stars might have been quite massive, possibly even very massive with M_star > 100 M_sun.

We briefly review the motivations for studying the formation of the first ``Population III'' stars and present recent results from our numerical simulations in this area. We discuss the new questions raised as a result of the simulations presented by us and others at this meeting.

This contribution summarizes briefly the main topics covered at this wide-ranging conference. Much of the evidence presented indicates that star formation occurs in discrete episodes or bursts that produce stellar groupings of all sizes, the fossil remnants of which make up the present stellar populations of galaxies.

The overall frequency and other statistical properties of binary systems suggest that star formation is intrinsically a complex and chaotic process, and that most binaries and single stars actually originate from the decay of multiple systems. Interactions between stars forming in close proximity to each other may play an important role in the star formation process itself, for example via tidally induced accretion from disks. Some of the energetic activity of newly formed stars could be due to bursts of rapid accretion triggered by interactions with close companions.

The first bound star-forming systems in the universe are predicted to form at redshifts of about 30 and to have masses of the order of 10^6 M_sun. Although their sizes and masses are similar to those of present star-forming regions, their temperatures are expected to be much higher because cooling is provided only by trace amounts of molecular hydrogen. Several recent simulations of the collapse and fragmentation of primordial clouds have converged on a thermal regime where the density is about 10^3-10^4 cm^{-3} and the temperature is about 300 K; under these conditions the Jeans mass is of the order of 10^3 M_sun, and all of the simulations show the formation of clumps with masses of this order. The temperatures in these clumps subsequently rise slowly as they collapse, so little if any further fragmentation is expected. As a result, the formation of predominantly massive or very massive stars is expected, and star formation with a normal present-day IMF seems very unlikely. The most massive early stars are expected to collapse to black holes, and these in turn are predicted to end up concentrated near the centers of present-day large galaxies. Such black holes may play a role in the origin of AGNs, and the heavy elements produced by somewhat less massive stars also formed at early times may play an important role in chemically enriching the inner parts of large galaxies and quasars.

In order to constrain the initial mass function (IMF) of the first generation of stars (Population III), we investigate the fragmentation properties of metal-free gas in the context of a hierarchical model of structure formation. We investigate the evolution of an isolated 3-sigma peak of mass 2x10^6 M_solar which collapses at z_coll=30 using Smoothed Particle Hydrodynamics. We find that the gas dissipatively settles into a rotationally supported disk which has a very filamentary morphology. The gas in these filaments is Jeans unstable with M_J~10^3 M_solar. Fragmentation leads to the formation of high density (n>10^8 cm^-3) clumps which subsequently grow in mass by accreting surrounding gas and by merging with other clumps up to masses of ~10^4 M_solar. This suggests that the very first stars were rather massive. We explore the complex dynamics of the merging and tidal disruption of these clumps by following their evolution over a few dynamical times.

How did star formation begin in the universe? Some of the questions addressed at this first meeting on "The First Stars" are summarized here from a theoretical perspective, and some brief comments are made on what we may have learned so far.

The current status of both the observational evidence and the theory of the stellar initial mass function (IMF) is reviewed, with particular attention to the two basic, apparently universal features shown by all observations of nearby stellar systems: (1) a characteristic stellar mass of the order of one solar mass, and (2) a power-law decline of the IMF at large masses similar to the original Salpeter law. Considerable evidence and theoretical work supports the hypothesis that the characteristic stellar mass derives from a characteristic scale of fragmentation in star-forming clouds which is essentially the Jeans scale as calculated from the typical temperature and pressure in molecular clouds. The power-law decline of the IMF at large masses suggests that the most massive stars are built up by scale-free accretion or accumulation processes, and the observed formation of these stars in dense clusters and close multiple systems suggests that interactions between dense prestellar clumps or protostars in forming clusters will play a role. A simple model postulating successive mergers of subsystems in a forming cluster accompanied by the accretion of a fraction of the residual gas by the most massive protostar during each merger predicts an upper IMF of power-law form and reproduces the Salpeter law with a plausible assumed accretion efficiency.

We summarize recent observational and theoretical progress aimed at understanding the origin of the stellar initial mass function (IMF) with specific focus on galactic star-forming regions. We synthesize data from various efforts to determine the IMF in very young, partially-embedded stellar clusters and find: i) no significant variations in the low-mass IMF have been observed between different star-forming regions; and ii) the mass distributions of young stars just emerging from molecular clouds are consistent with having been drawn from the IMF derived from field stars in the solar neighborhood. These results apply only to gross characterizations of the IMF (e.g. the ratio of high to low mass stars); present observations do not rule out more subtle regional differences. Further studies are required in order to assess whether or not there is evidence for a universal turnover near the hydrogen--burning limit. We also provide a general framework for discussing theories of the IMF, and summarize recent work on several physical mechanisms which could play a role in determining the form of the stellar initial mass function.

In order to constrain the initial mass function of the first generation of stars (Population III), one has to study the fragmentation properties of primordial gas. We present results from 3D simulations, based on Smoothed Particle Hydrodynamics, which explore the idea that small-scale fluctuations in the (cold) dark matter recreate a filamentary and clumpy structure in the gas component on scales smaller than the initial Jeans mass, where all primordial fluctuations would have been wiped out.

It has frequently been suggested that the stellar IMF in galaxies was top-heavy at early times. This would be plausible physically if the IMF depends on a mass scale such as the Jeans mass that was higher at earlier times because of the generally higher temperatures that were present then. In this paper it is suggested, on the basis of current evidence and theory, that the IMF has a universal Salpeter-like form at the upper end but flattens below a characteristic stellar mass that may vary with time. Much of the evidence that has been attributed to a top-heavy early IMF, including the ubiquitous G-dwarf problem, the high abundance of heavy elements in clusters of galaxies, and the high rate of formation of massive stars in high-redshift galaxies, can be accounted for with such an IMF if the characteristic stellar mass was several times higher during the early stages of galaxy evolution. However, significant variations in the mass-to-light ratios of galaxies and large amounts of dark matter in stellar remnants are not as easily explained in this way because they require more extreme and less plausible assumptions about the form and variability of the IMF. Metal-free 'population III' stars are predicted to have an IMF that consists exclusively of massive stars, and they could help to account for some of the evidence that has been attributed to a top-heavy early IMF, as well as contributing importantly to the energetics and chemical enrichment of the early universe.