Search

We present analysis of the mass-metallicity relation (MZR) for a sample of 67 [OIII]-selected star-forming galaxies at a redshift range of $z=1.99 - 2.32$ ($z_{\text{med}} = 2.16$) using \emph{Hubble Space Telescope} Wide Field Camera 3 grism spectroscopy from the Quasar Sightline and Galaxy Evolution (QSAGE) survey. Metallicities were determined using empirical gas-phase metallicity calibrations based on the strong emission lines [OII]3727,3729, [OIII]4959,5007 and H$\beta$. Star-forming galaxies were identified, and distinguished from active-galactic nuclei, via Mass-Excitation diagrams. Using $z\sim0$ metallicity calibrations, we observe a negative offset in the $z=2.2$ MZR of $\approx -0.51$ dex in metallicity when compared to locally derived relationships, in agreement with previous literature analysis. A similar offset of $\approx -0.46$ dex in metallicity is found when using empirical metallicity calibrations that are suitable out to $z\sim5$, though our $z=2.2$ MZR, in this case, has a shallower slope. We find agreement between our MZR and those predicted from various galaxy evolution models and simulations. Additionally, we explore the extended fundamental metallicity relation (FMR) which includes an additional dependence on star formation rate (SFR). Our results consistently support the existence of the FMR, as well as revealing an offset of $0.28\pm0.04$ dex in metallicity compared to locally-derived relationships, consistent with previous studies at similar redshifts. We interpret the negative correlation with SFR at fixed mass, inferred from an FMR existing for our sample, as being caused by the efficient accretion of metal-poor gas fuelling SFR at cosmic noon., Comment: 15 pages, 10 figures. This paper has been accepted for publication in MNRAS

We demonstrate that the UV brightest quasars at z=1-2 live in overdense environments. This is based on an analysis of deep Hubble Space Telescope WFC3 G141 grism spectroscopy of the galaxies along the lines-of-sight to UV luminous quasars in the redshift range z=1-2. This constitutes some of the deepest grism spectroscopy performed by WFC3, with 4 roll angles spread over a year of observations to mitigate the effect of overlapping spectra. Of the 12 quasar fields studied, 8 display evidence for a galaxy overdensity at the redshift of the quasar. One of the overdensities, PG0117+213 at z=1.50, has potentially 36 spectroscopically confirmed members, consisting of 19 with secure redshifts and 17 with single-line redshifts, within a cylinder of radius ~700 kpc. Its halo mass is estimated to be log (M/Msol)=14.7. This demonstrates that spectroscopic and narrow-band observations around distant UV bright quasars may be an excellent route for discovering protoclusters. Our findings agree with previous hints from statistical observations of the quasar population and theoretical works, as feedback regulated black hole growth predicts a correlation between quasar luminosity and halo mass. We also present the high signal-to-noise rest-frame optical spectral and photometric properties of the quasars themselves., Comment: Published in MNRAS. This is the author accepted version. A small number of typos were corrected after journal proof stage

The Widefield ASKAP L-band Legacy All-sky Blind surveY (WALLABY) is a next-generation survey of neutral hydrogen (HI) in the Local Universe. It uses the widefield, high-resolution capability of the Australian Square Kilometer Array Pathfinder (ASKAP), a radio interferometer consisting of 36 x 12-m dishes equipped with Phased-Array Feeds (PAFs), located in an extremely radio-quiet zone in Western Australia. WALLABY aims to survey three-quarters of the sky (-90 degr < Dec < +30 degr) to a redshift of z < 0.26, and generate spectral line image cubes at ~30 arcsec resolution and ~1.6 mJy/beam per 4 km/s channel sensitivity. ASKAP's instantaneous field of view at 1.4 GHz, delivered by the PAF's 36 beams, is about 30 sq deg. At an integrated signal-to-noise ratio of five, WALLABY is expected to detect over half a million galaxies with a mean redshift of z ~ 0.05 (~200 Mpc). The scientific goals of WALLABY include: (a) a census of gas-rich galaxies in the vicinity of the Local Group; (b) a study of the HI properties of galaxies, groups and clusters, in particular the influence of the environment on galaxy evolution; and (c) the refinement of cosmological parameters using the spatial and redshift distribution of low-bias gas-rich galaxies. For context we provide an overview of previous large-scale HI surveys. Combined with existing and new multi-wavelength sky surveys, WALLABY will enable an exciting new generation of panchromatic studies of the Local Universe. - First results from the WALLABY pilot survey are revealed, with initial data products publicly available in the CSIRO ASKAP Science Data Archive (CASDA)., Comment: Accepted for publication in ApSS (38 pages, 14 figures), see also https://www.atnf.csiro.au/research/WALLABY/ - Contact email: Baerbel.Koribalski@csiro.au

About 35 years ago a class of galaxies with unusually strong Balmer absorption lines and weak emission lines was discovered in distant galaxy clusters. These objects, alternatively referred to as post-starburst, E+A or k+a galaxies, are now known to occur in all environments and at all redshifts, with many exhibiting compact morphologies and low-surface brightness features indicative of past galaxy mergers. They are commonly thought to represent galaxies that are transitioning from blue to red sequence, making them critical to our understanding of the origins of galaxy bimodality. However, recent observational studies have questioned this simple interpretation. From observations alone, it is challenging to disentangle the different mechanisms that lead to the quenching of star formation in galaxies. Here we present examples of three different evolutionary pathways that lead to galaxies with strong Balmer absorption lines in the EAGLE simulation: classical blue-to-red quenching, blue-to-blue cycle and red-to-red rejuvenation. The first two are found in both post-starburst galaxies and galaxies with truncated star formation. Each pathway is consistent with scenarios hypothesised for observational samples. The fact that 'post-starburst' signatures can be attained via various evolutionary channels explains the diversity of observed properties, and lends support to the idea that slower quenching channels are important at low redshift., Comment: Published in Nature Astronomy, 25th March 2019 (submitted version available on arxiv only). The final accepted version is available for all to read here https://rdcu.be/bsQIC. Corresponding author: V. Wild. 11 pages, 4 figures

We present the first results from a study of OVI absorption around galaxies at $z<1.44$ using data from a near-infrared grism spectroscopic Hubble Space Telescope Large Program, the Quasar Sightline and Galaxy Evolution (QSAGE) survey. QSAGE is the first grism galaxy survey to focus on the circumgalactic medium at $z\sim1$, providing a blind survey of the galaxy population. Using the first of 12 fields, we provide details of the reduction methods, in particular the handling of the deep grism data which uses multiple position angles to minimise the effects of contamination from overlapping traces. The resulting galaxy sample is H$\alpha$ flux limited ($f({\rm H\alpha}) > 2\times10^{-17}$ erg s$^{-1}$ cm$^{-2}$) at 0.680.68, we find 5 OVI absorption systems along the line of sight with identified galaxies lying at impact parameters of $b\approx100-350$ kpc (proper), whilst we find a further 13 galaxies with no significant associated OVI absorption (i.e. $N({\rm OVI})<10^{13.5-14}$ cm$^{-2}$) in the same impact parameter and redshift range. We find a large scatter in the stellar mass and star-formation rates of the closest galaxies with associated OVI. Whilst one of the OVI absorber systems is found to be associated with a low mass galaxy group at $z\approx1.08$, we infer that the detected OVI absorbers typically lie in the proximity of dark matter halos of masses $10^{11.5} {\rm M_\odot}\lesssim M_{\rm halo}\lesssim10^{12} {\rm M_\odot}$., Comment: 22 pages, 18 figures, 4 tables, submitted to MNRAS 14th Sept 2018

This work presents a study of galactic outflows driven by stellar feedback. We extract main sequence disc galaxies with stellar mass $10^9\le$ M$_{\star}/$M$_{\odot} \le 5.7\times10^{10}$ at redshift $z=0$ from the highest resolution cosmological simulation of the Evolution and Assembly of GaLaxies and their Environments (EAGLE) set. Synthetic gas rotation velocity and velocity dispersion ($\sigma$) maps are created and compared to observations of disc galaxies obtained with the Sydney-AAO Multi-object Integral field spectrograph (SAMI), where $\sigma$-values greater than $150$ km s$^{-1}$ are most naturally explained by bipolar outflows powered by starburst activity. We find that the extension of the simulated edge-on (pixelated) velocity dispersion probability distribution depends on stellar mass and star formation rate surface density ($\Sigma_{\rm SFR}$), with low-M$_{\star}/$low-$\Sigma_{\rm SFR}$ galaxies showing a narrow peak at low $\sigma$ ($\sim30$ km s$^{-1}$) and more active, high-M$_{\star}/$high-$\Sigma_{\rm SFR}$ galaxies reaching $\sigma>150$ km s$^{-1}$. Although supernova-driven galactic winds in the EAGLE simulations may not entrain enough gas with T $<10^5$ K compared to observed galaxies, we find that gas temperature is a good proxy for the presence of outflows. There is a direct correlation between the thermal state of the gas and its state of motion as described by the $\sigma$-distribution. The following equivalence relations hold in EAGLE: $i)$ low-$\sigma$ peak $\,\Leftrightarrow\,$ disc of the galaxy $\,\Leftrightarrow\,$ gas with T $<10^5$ K; $ii)$ high-$\sigma$ tail $\,\Leftrightarrow\,$ galactic winds $\,\Leftrightarrow\,$ gas with T $\ge 10^5$ K., Comment: 18 pages, 12 figures, accepted for publication in MNRAS

We present the evolution of galaxy sizes, from redshift 2 to 0, for actively star forming and passive galaxies in the cosmological hydrodynamical 1003 cMpc3 simulation of the EAGLE project. We find that the sizes increase with stellar mass , but that the relation weakens with increasing redshift. Separating galaxies by their star formation activity, we find that passive galaxies are typically smaller than active galaxies at fixed stellar mass. These trends are consistent with those found in observations and the level of agreement between the predicted and observed size - mass relation is of order 0.1 dex for z < 1 and 0.2-0.3 dex from redshift 1 to 2. We use the simulation to compare the evolution of individual galaxies to that of the population as a whole. While the evolution of the size-stellar mass relation for active galaxies provides a good proxy for the evolution of individual galaxies, the evolution of individual passive galaxies is not well represented by the observed size - mass relation due to the evolving number density of passive galaxies. Observations of z \approx 2 galaxies have revealed an abundance of massive red compact galaxies, that depletes below z \approx 1. We find that a similar population forms naturally in the simulation. Comparing these galaxies to their z = 0 descendants, we find that all compact galaxies grow in size due to the high-redshift stars migrating outwards. Approximately 60% of the compact galaxies increase in size further due to renewed star formation and/or mergers., Comment: 14 pages, 10 figures

We study high-resolution hydrodynamic simulations of Milky Way type galaxies obtained within the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) project, and identify the those that best satisfy observational constraints on the Milky Way total stellar mass, rotation curve, and galaxy shape. Contrary to mock galaxies selected on the basis of their total virial mass, the Milky Way analogues so identified consistently exhibit very similar dark matter profiles inside the solar circle, therefore enabling more accurate predictions for indirect dark matter searches. We find in particular that high resolution simulated haloes satisfying observational constraints exhibit, within the inner few kiloparsecs, dark matter profiles shallower than those required to explain the so-called Fermi GeV excess via dark matter annihilation., Comment: 33 pages, 15 figures, 4 tables. Matches published version

We investigate the evolution of galaxy masses and star formation rates in the Evolution and Assembly of Galaxies and their Environment (EAGLE) simulations. These comprise a suite of hydrodynamical simulations in a $\Lambda$CDM cosmogony with subgrid models for radiative cooling, star formation, stellar mass loss, and feedback from stars and accreting black holes. The subgrid feedback was calibrated to reproduce the observed present-day galaxy stellar mass function and galaxy sizes. Here we demonstrate that the simulations reproduce the observed growth of the stellar mass density to within 20 per cent. The simulation also tracks the observed evolution of the galaxy stellar mass function out to redshift z = 7, with differences comparable to the plausible uncertainties in the interpretation of the data. Just as with observed galaxies, the specific star formation rates of simulated galaxies are bimodal, with distinct star forming and passive sequences. The specific star formation rates of star forming galaxies are typically 0.2 to 0.4 dex lower than observed, but the evolution of the rates track the observations closely. The unprecedented level of agreement between simulation and data makes EAGLE a powerful resource to understand the physical processes that govern galaxy formation., Comment: Accepted for publication by MNRAS; 18 pages, 7 figures; comments welcome

Feedback from energy liberated by gas accretion onto black holes (BHs) is an attractive mechanism to explain the exponential cut-off at the massive end of the galaxy stellar mass function (SMF). Semi-analytic models of galaxy formation in which this form of feedback is assumed to suppress cooling in haloes where the gas cooling time is large compared to the dynamical time do indeed achieve a good match to the observed SMF. Furthermore, hydrodynamic simulations of individual halos in which gas is assumed to accrete onto the central BH at the Bondi rate have shown that a self-regulating regime is established in which the BH grows just enough to liberate an amount of energy comparable to the thermal energy of the halo. However, this process is efficient at suppressing the growth not only of massive galaxies but also of galaxies like the Milky Way, leading to disagreement with the observed SMF. The Bondi accretion rate, however, is inappropriate when the accreting material has angular momentum. We present an improved accretion model that takes into account the circularisation and subsequent viscous transport of infalling material and include it as a "subgrid" model in hydrodynamic simulations of the evolution of halos with a wide range of masses. The resulting accretion rates are generally low in low mass ($\lsim 10^{11.5} \msun$) halos, but show outbursts of Eddington-limited accretion during galaxy mergers. During outbursts these objects strongly resemble quasars. In higher mass haloes, gas accretion occurs continuously, typically at $~10$ % of the Eddington rate, which is conducive to the formation of radio jets. The resulting dependence of the accretion behaviour on halo mass induces a break in the relation between galaxy stellar mass and halo mass in these simulations that matches observations., Comment: 24 pages, 14 figures. Accepted to MNRAS for publication. Updated version with typos corrected

We present the results of a search for molecular gas emission from a star-forming galaxy at z = 4.9. The galaxy benefits from magnification of 22 +/- 5x due to strong gravitational lensing by the foreground cluster MS1358+62. We target the CO(5-4) emission at a known position and redshift from existing Hubble Space Telescope/ACS imaging and Gemini/NIFS [OII]3727 imaging spectroscopy, and obtain a tentative detection at the 4.3sigma level with a flux of 0.104 +/- 0.024Jkm/s. From the CO line luminosity and assuming a CO-to-H2 conversion factor alpha=2, we derive a gas mass M_gas ~ 1^{+1}_{-0.6} x 10^9 M_sun. Combined with the existing data, we derive a gas fraction Mgas/(Mgas + M*) = 0.59^{+0.11}_{-0.06}. The faint line flux of this galaxy highlights the difficulty of observing molecular gas in representative galaxies at this epoch, and suggests that routine detections of similar galaxies in the absence of gravitational lensing will remain challenging even with ALMA in full science operations., Comment: 14 pages, 4 figures, published in ApJL

We present HST/WFC3 narrowband imaging of the H-alpha emission in a sample of eight gravitationally-lensed galaxies at z = 1 - 1.5. The magnification caused by the foreground clusters enables us to obtain a median source plane spatial resolution of 360pc, as well as providing magnifications in flux ranging from ~10x to ~50x. This enables us to identify resolved star-forming HII regions at this epoch and therefore study their H-alpha luminosity distributions for comparisons with equivalent samples at z ~ 2 and in the local Universe. We find evolution in the both luminosity and surface brightness of HII regions with redshift. The distribution of clump properties can be quantified with an HII region luminosity function, which can be fit by a power law with an exponential break at some cut-off, and we find that the cut-off evolves with redshift. We therefore conclude that `clumpy' galaxies are seen at high redshift because of the evolution of the cut-off mass; the galaxies themselves follow similar scaling relations to those at z = 0, but their HII regions are larger and brighter and thus appear as clumps which dominate the morphology of the galaxy. A simple theoretical argument based on gas collapsing on scales of the Jeans mass in a marginally unstable disk shows that the clumpy morphologies of high-z galaxies are driven by the competing effects of higher gas fractions causing perturbations on larger scales, partially compensated by higher epicyclic frequencies which stabilise the disk., Comment: 16 pages, 9 figures, 2 tables, accepted for publication in MNRAS

The observed stellar mass function (SMF) is very different to the halo mass function predicted by Lambda-CDM, and it is widely accepted that this is due to energy feedback from supernovae and black holes. However, the strength and form of this feedback is not understood. In this paper, we use the phenomenological model GALFORM to explore how galaxy formation depends on the strength and halo mass dependence of feedback. We focus on 'expulsion' models in which the wind mass loading, beta, is proportional to 1/\vdisk^n, with n=0,1,2 and contrast these models with the successful Bower et al.\ 2008 model (B8W7). A crucial development is that our code explicitly accounts for the recapture of expelled gas as the system's halo mass (and thus gravitational potential) increases. We find that a model with modest wind speed but high mass loading matches the flat portion of the SMF. When combined with AGN feedback, the model provides a good description of the observed SMF above 10^9 h^-1 Msol. However, in the expulsion models, the brightest galaxies are assembled more recently than in B8W7, and the specific star formation rates of galaxies decrease strongly with decreasing stellar mass. The expulsion models also tend to have a cosmic star formation density that is dominated by lower mass galaxies at z=1-3, and dominated high mass galaxies at low redshift. These trends are in conflict with observational data, but the comparison highlights some deficiencies of the B8W7 model also. The experiments in this paper give us important physical insight to the impact of the feedback process on the formation histories of galaxies, but the strong mass dependence of feedback adopted in B8W7 still appears to provide the most promising description of the observed universe., Comment: 27 pages (31 color figures). Accepted for publication in MNRAS. E-mail to obtain Model catalogues

We compare the results of thirteen cosmological gasdynamical codes used to simulate the formation of a galaxy in the LCDM structure formation paradigm. The various runs differ in their hydrodynamical treatment (SPH, moving-mesh and AMR) but share the same initial conditions and adopt their latest published model of cooling, star formation and feedback. Despite the common halo assembly history, we find large code-to-code variations in the stellar mass, size, morphology and gas content of the galaxy at z=0, due mainly to the different implementations of feedback. Compared with observation, most codes tend to produce an overly massive galaxy, smaller and less gas-rich than typical spirals, with a massive bulge and a declining rotation curve. A stellar disk is discernible in most simulations, though its prominence varies widely from code to code. There is a well-defined trend between the effects of feedback and the severity of the disagreement with observation. Models that are more effective at limiting the baryonic mass of the galaxy come closer to matching observed galaxy scaling laws, but often to the detriment of the disk component. Our conclusions hold at two different numerical resolutions. Some differences can also be traced to the numerical techniques: more gas seems able to cool and become available for star formation in grid-based codes than in SPH. However, this effect is small compared to the variations induced by different feedback prescriptions. We conclude that state-of-the-art simulations cannot yet uniquely predict the properties of the baryonic component of a galaxy, even when the assembly history of its host halo is fully specified. Developing feedback algorithms that can effectively regulate the mass of a galaxy without hindering the formation of high-angular momentum stellar disks remains a challenge., Comment: 28 pages, 16 figures, MNRAS accepted. High resolution figures and additional images can be found at http://www.aip.de/People/cscannapieco/aquila/

(Abridged) We study the phase-space distribution of satellite galaxies associated with late-type galaxies in the GIMIC suite of simulations. GIMIC consists of re-simulations of 5 cosmologically representative regions from the Millennium simulation, which have higher resolution and incorporate baryonic physics. Whilst the disc of the galaxy is well aligned with the inner regions (r ~ 0.1r_200) of the dark matter halo, both in shape and angular momentum, there can be substantial misalignments at larger radii (r ~ r_200). Misalignments of > 45 deg are seen in ~ 30% of our sample. We find that the satellite population aligns with the shape (and angular momentum) of the outer dark matter halo. However, the alignment with the galaxy is weak owing to the mismatch between the disc and dark matter halo. Roughly 20% of the satellite systems with ten bright galaxies within r_200 exhibit a polar spatial alignment with respect to the galaxy --- an orientation reminiscent of the classical satellites of the Milky Way. We find that a small fraction (~10%) of satellite systems show evidence for rotational support which we attribute to group infall. There is a bias towards satellites on prograde orbits relative to the spin of the dark matter halo (and to a lesser extent with the angular momentum of the disc). This preference towards co-rotation is stronger in the inner regions of the halo where the most massive satellites accreted at relatively early times are located. We attribute the anisotropic spatial distribution and angular momentum bias of the satellites at z=0 to their directional accretion along the major axes of the dark matter halo. The satellite galaxies have been accreted relatively recently compared to the dark matter mass and have experienced less phase-mixing and relaxation --- the memory of their accretion history can remain intact to z=0., Comment: 19 pages, 17 Figures. Accepted by MNRAS

(Abridged) The thermal state of the intracluster medium results from a competition between gas cooling and heating. The heating comes from two distinct sources: gravitational heating from the collapse of the dark matter halo and thermal input from galaxy/black hole formation. However, a long standing problem has been that cosmological simulations based on smoothed particle hydrodynamics (SPH) and Eulerian mesh codes predict different results even when cooling and galaxy/black hole heating are switched off. Clusters formed in SPH simulations show near powerlaw entropy profiles, while those formed in mesh simulations develop a core and do not allow gas to reach such low entropies. Since the cooling rate is closely connected to the minimum entropy of the gas, the differences are of potentially key importance. In this paper, we investigate the origin of this discrepancy. By comparing simulations run using the GADGET-2 SPH code and the FLASH adaptive Eulerian mesh code, we show that the discrepancy arises during the idealised merger of two clusters. The difference is not sensitive to the resolution of our simulations, nor is it is due differences in the gravity solvers, Galilean non-invariance of the mesh code, or an effect of unsuitable artificial viscosity in the SPH code. Instead, we find that the difference is inherent to the treatment of eddies and fluid instabilities. These are suppressed in the SPH simulations, while the cluster mergers generate strong vortices in the mesh simulations that efficiently mix the fluid and erase the low entropy gas. Consequently, particles in the SPH simulations retain a close connection to their initial entropy, while this connection is much weaker in the mesh simulations. We discuss the potentially profound implications of these results., Comment: 19 pages, 10 figures, Accepted for publication in MNRAS

We present Hubble Space Telescope/Wide Field Camera 3 narrow-band imaging of the Hα emission in a sample of eight gravitationally lensed galaxies at z = 1-1.5. The magnification caused by the foreground clusters enables us to obtain a median source plane spatial resolution of 360 pc, as well as providing magnifications in flux ranging from approximately 10× to approximately 50×. This enables us to identify resolved star-forming HII regions at this epoch and therefore study their Hα luminosity distributions for comparisons with equivalent samples at z approximately 2 and in the local Universe. We find evolution in the both luminosity and surface brightness of HII regions with redshift. The distribution of clump properties can be quantified with an HII region luminosity function, which can be fit by a power law with an exponential break at some cut-off, and we find that the cut-off evolves with redshift. We therefore conclude that 'clumpy' galaxies are seen at high redshift because of the evolution of the cut-off mass; the galaxies themselves follow similar scaling relations to those at z = 0, but their HII regions are larger and brighter and thus appear as clumps which dominate the morphology of the galaxy. A simple theoretical argument based on gas collapsing on scales of the Jeans mass in a marginally unstable disc shows that the clumpy morphologies of high-z galaxies are driven by the competing effects of higher gas fractions causing perturbations on larger scales, partially compensated by higher epicyclic frequencies which stabilize the disc.

Appendix A, "Data Analysis When Students Change from Smaller to Larger Schools between the Pre-test and Post-test of the First Year," discusses a methodological problem involved in planning the approximate number of questionnaires to be administered in program. Appendix B, "Estimates of Sample Sizes and Data Collection Costs," includes five tables. Tables 1 and 2 show estimated numbers of schools, students, and teachers for each year of the program. Table 3 shows an approximate cost estimate for data collection, and Table 4 indicates how sample sizes would change if the experimental version of the study were used. Table 5 gives the costs of implementation of the experimental innovations. Appendix C contains basic questionnaires and test batteries to be used in the various surveys. These include student questionnaires, noncognitive outcomes, teacher questionnaire, principal questionnaire, superintendent questionnaire, parent questionnaire, community leadership questionnaires, bibliography, and classroom observation instruments. Appendix D describes Jessen's Method 2 for unequal probability sampling as well as for probability lattice sampling. Also discussed are options that enhance flexibility, involving cell size, zero cells, multidimensions, unequal margins, and noninteger margins. The "Replicability of Results" problem is also examined. (Author/JM)

The U.S. Commission on Civil Rights asked the Rand Corporation to prepare a design for an exhaustive research program on school desegregation. This is Volume 2 of the final report. The 11 sections of this volume describe a general design for a research program in school desegregation. Section 6 describes a community reaction survey. In Section 7 we outline a one-year longitudinal study of student attitudes and achievement in a sample of 1600 elementary, junior high and high schools. In Sections 8 and 9 we describe four substudies: each takes a subsample of schools and pursues a particular topic more intensively. Section 10 discusses the sample. Section 11 is a design for a study of the impact of desegregation policies on the exodus of whites from central cities. Section 12 discusses how the studies of innovative schools can be converted to the more powerful methodology of the randomized experiment. Section 13 discusses how three existing panels of young adults can be used to test hypotheses about the impact of different types of desegregation on social-psychological laboratory experiments to test hypotheses about the dynamics of race relations in schools. Finally, Sections 15 and 16 discuss general problems of the program. (Author/JM)