We presented the first particle based, Lagrangian code that can follow the collisional/accretional/dynamical evolution of a large number of km-sized planetesimals through the entire growth process to become planets. We refer to it as the 'Lagrangian Integrator for Planetary Accretion and Dynamics' or LIPAD. LIPAD is built on top of SyMBA, which is a symplectic $N$-body integrator. In order to handle the very large number of planetesimals required by planet formation simulations, we introduce the concept of a `tracer' particle. Each tracer is intended to represent a large number of disk particles on roughly the same orbit and size as one another, and is characterized by three numbers: the physical radius, the bulk density, and the total mass of the disk particles represented by the tracer. We developed statistical algorithms that follow the dynamical and collisional evolution of the tracers due to the presence of one another. The tracers mainly dynamically interact with the larger objects (`planetary embryos') in the normal N-body way. LIPAD's greatest strength is that it can accurately model the wholesale redistribution of planetesimals due to gravitational interaction with the embryos, which has recently been shown to significantly affect the growth rate of planetary embryos . We verify the code via a comprehensive set of tests which compare our results with those of Eulerian and/or direct N-body codes.

The discovery of over 400 extrasolar planets allows us to statistically test our understanding of formation and dynamics of planetary systems via numerical simulations. Traditional N-body simulations of multiple-planet systems without gas disks have successfully reproduced the eccentricity (e) distribution of the observed systems, by assuming that the planetary systems are relatively closely packed when the gas disk dissipates, so that they become dynamically unstable within the stellar lifetime. However, such studies cannot explain the small semi-major axes (a) of extrasolar planetary systems, if planets are formed, as the standard planet formation theory suggests, beyond the ice line. In this paper, we numerically study the evolution of three-planet systems in dissipating gas disks, and constrain the initial conditions that reproduce the observed semi-major axis and eccentricity distributions simultaneously. We adopt the initial conditions that are motivated by the standard planet formation theory, and self-consistently simulate the disk evolution, and planet migration by using a hybrid N-body and 1D gas disk code. We also take account of eccentricity damping, and investigate the effect of saturation of corotation resonances on the evolution of planetary systems. We find that the semi-major axis distribution is largely determined in a gas disk, while the eccentricity distribution is determined after the disk dissipation. We also find that there may be an optimum disk mass which leads to the observed a-e distribution. Our simulations generate a larger fraction of planetary systems trapped in mean-motion resonances (MMRs) than the observations, indicating that the disk's perturbation to the planetary orbits may be important to explain the observed rate of MMRs. We also find much lower occurrence of planets on retrograde orbits than the current observations of close-in planets suggest.

One of the most challenging problems we face in our understanding of planet formation is how Jupiter and Saturn could have formed before the the solar nebula dispersed. The most popular model of giant planet formation is the so-called 'core accretion' model. In this model a large planetary embryo formed first, mainly by two-body accretion. This is then followed by a period of inflow of nebular gas directly onto the growing planet. The core accretion model has an Achilles heel, namely the very first step. We have undertaken the most comprehensive study of this process to date. In this study we numerically integrate the orbits of a number of planetary embryos embedded in a swarm of planetesimals. In these experiments we have included: 1) aerodynamic gas drag, 2) collisional damping between planetesimals, 3) enhanced embryo cross-sections due to their atmospheres, 4) planetesimal fragmentation, and 5) planetesimal driven migration. We find that the gravitational interaction between the embryos and the planetesimals lead to the wholesale redistribution of material - regions are cleared of material and gaps open near the embryos. Indeed, in 90% of our simulations without fragmentation, the region near that embryos is cleared of planetesimals before much growth can occur. The remaining 10%, however, the embryos undergo a burst of outward migration that significantly increases growth. On timescales of ~100,000 years, the outer embryo can migrate ~6 AU and grow to roughly 30 Earth-masses. We also find that the inclusion of planetesimal fragmentation tends to inhibit growth.

The differential migration of two planets due to planet-disk interaction can result in capture into the 2:1 eccentricity-type mean-motion resonances. Both the sequence of 2:1 eccentricity resonances that the system is driven through by continued migration and the possibility of a subsequent capture into the 4:2 inclination resonances are sensitive to the migration rate within the range expected for type II migration. If the migration rate is fast, the resonant pair can evolve into a family of 2:1 eccentricity resonances different from those found by Lee (2004). This new family has outer orbital eccentricity e_2 > 0.4-0.5, asymmetric librations of both eccentricity resonance variables, and orbits that intersect if they are exactly coplanar. Although this family exists for an inner-to-outer planet mass ratio m_1/m_2 > 0.2, it is possible to evolve into this family by fast migration only for m_1/m_2 > 2. Thommes and Lissauer (2003) have found that a capture into the 4:2 inclination resonances is possible only for m_1/m_2 < 2. We show that this capture is also possible for m_1/m_2 > 2 if the migration rate is slightly slower than that adopted by Thommes and Lissauer. There is significant theoretical uncertainty in both the sign and the magnitude of the net effect of planet-disk interaction on the orbital eccentricity of a planet. If the eccentricity is damped on a timescale comparable to or shorter than the migration timescale, e_2 may not be able to reach the values needed to enter either the new 2:1 eccentricity resonances or the 4:2 inclination resonances. Thus, if future observations of extrasolar planetary systems were to reveal certain combinations of mass ratio and resonant configuration, they would place a constraint on the strength of eccentricity damping during migration, as well as on the migration rate itself.

Mean-motion resonances (MMRs) are likely to play an important role both during and after the lifetime of a protostellar gas disk. We study the dynamical evolution and stability of planetary systems containing two giant planets on circular orbits near a 2:1 resonance and closer. We find that by having the outer planet migrate inward, the two planets can capture into either the 2:1, 5:3, or 3:2 MMR. We use direct numerical integrations of ~1000 systems in which the planets are initially locked into one of these resonances and allowed to evolve for up to ~10^7 yr. We find that the final eccentricity distribution in systems which ultimately become unstable gives a good fit to observed exoplanets. Next, we integrate ~500 two-planet systems in which the outer planet is driven to continuously migrate inward, resonantly capturing the inner; the systems are evolved until either instability sets in or the planets reach the star. We find that although the 5:3 resonance rapidly becomes unstable under migration, the 2:1 and 3:2 are very stable. Thus the lack of observed exoplanets in resonances closer than 2:1, if it continues to hold up, may be a primordial signature of the planet formation process.

The growth of Jovian mass planets during migration in their protoplanetary disks is one of the most important problems that needs to be solved in light of observations of the exosolar planets. Studies of the migration of planets in standard gas disk models routinely show that migration is too fast to form Jovian planets, and that such migrating planetary cores generally plunge into the central stars in less than a Myr. In previous work, we have shown that a poorly ionized, less viscous region in a protoplanetary disk called a dead zone slows down the migration of fixed-mass planets. In this paper, we extend our numerical calculations to include dead zone evolution along with the disk, as well as planet formation via accretion of rocky and gaseous materials. Using our symplectic-integrator-gas dynamics code, we find that dead zones, even in evolving disks wherein migrating planets grow by accretion, still play a fundamental role in saving planetary systems. We demonstrate that Jovian planets form within 2.5 Myr for disks that are ten times more massive than a minimum mass solar nebula (MMSN) with an opacity reduction and without slowing down migration artificially. Our simulations indicate that protoplanetary disks with an initial mass comparable to the MMSN only produce Neptunian mass planets. We also find that planet migration does not help core accretion as much in the oligarchic planetesimal accretion scenario as it was expected in the runaway accretion scenario. Therefore we expect that an opacity reduction (or some other mechanisms) is needed to solve the formation timescale problem even for migrating protoplanets, as long as we consider the oligarchic growth. We also point out a possible role of a dead zone in explaining long-lived, strongly accreting gas disks.

The ensemble of now more than 250 discovered planetary systems displays a wide range of masses, orbits and, in multiple systems, dynamical interactions. These represent the end point of a complex sequence of events, wherein an entire protostellar disk converts itself into a small number of planetary bodies. Here, we present self-consistent numerical simulations of this process, which produce results in agreement with some of the key trends observed in the properties of the exoplanets. Analogs to our own solar system do not appear to be common, originating from disks near the boundary between barren and (giant) planet-forming.

We revisit the "dynamical shakeup" model of Solar System terrestrial planet formation, wherein the whole process is driven by the sweeping of Jupiter's secular resonance as the gas disk is removed. Using a large number of 0.5 Gyr-long N-body simulations, we investigate the different outcomes produced by such a scenario. We confirm that in contrast to existing models, secular resonance sweeping combined with tidal damping by the disk gas can reproduce the low eccentricities and inclinations, and high radial mass concentration, of the Solar System terrestrial planets. At the same time, this also drives the final assemblage of the planets on a timescale of several tens of millions of years, an order of magnitude faster than inferred from previous numerical simulations which neglected these effects, but possibly in better agreement with timescales inferred from cosmochemical data. In addition, we find that significant delivery of water-rich material from the outer Asteroid Belt is a natural byproduct.

Context: Many current and future surveys aim to detect the highest redshift (z >~ 7) sources through their Lyman-alpha (Ly-alpha) emission, using the narrow-band imaging method. However, to date the surveys have only yielded non-detections and upper limits as no survey has reached the necessary combination of depth and area to detect these very young star forming galaxies. Aims: We aim to calculate model luminosity functions and mock surveys of Ly-alpha emitters at z >~ 7 based on a variety of approaches. Methods: We calculate model luminosity functions at different redshifts based on three different approaches: a semi-analytical model based on CDM, a simple phenomenological model, and an extrapolation of observed Schechter functions at lower redshifts. The results of the first two models are compared with observations made at redshifts z ~ 5.7 and z ~ 6.5, and they are then extrapolated to higher redshift. Results: We present model luminosity functions for redshifts between z = 7 - 12.5 and give specific number predictions for future planned or possible narrow-band surveys for Ly-alpha emitters. We also investigate what constraints future observations will be able to place on the Ly-alpha luminosity function at very high redshift. Conclusion: It should be possible to observe z = 7 - 10 Ly-alpha emitters with present or near-future instruments if enough observing time is allocated. In particular, large area surveys such as ELVIS (Emission Line galaxies with VISTA Survey) will be useful in collecting a large sample. However, to get a large enough sample to constrain well the z >= 10 Ly-alpha luminosity function, instruments further in the future, such as an ELT, will be necessary.

We show that interaction with a gas disk may produce young planetary systems with closely-spaced orbits, stabilized by mean-motion resonances between neighbors. On longer timescales, after the gas is gone, interaction with a remnant planetesimal disk tends to pull these configurations apart, eventually inducing dynamical instability. We show that this can lead to a variety of outcomes; some cases resemble the Solar System, while others end up with high-eccentricity orbits reminiscent of the observed exoplanets. A similar mechanism has been previously suggested as the cause of the lunar Late Heavy Bombardment. Thus, it may be that a large-scale dynamical instability, with more or less cataclysmic results, is an evolutionary step common to many planetary systems, including our own.

We compute the impact of dark energy at last scattering on measurements of baryon acoustic oscillations (BAOs). We show that an early dark energy component can contribute a systematic uncertainty to BAO measurements of up to 2.5%. Whilst this effect turns out to only slightly affect current BAO surveys, the results of future BAO surveys might become biased. We find that BAO surveys alone appear unable to resolve this systematic uncertainty, so supplementary measurements are necessary.

The tidal interaction between a disk and a planet leads to the planet's migration. A long-standing question regarding this mechanism is how to stop the migration before planets plunge into their central stars. In this paper, we propose a new, simple mechanism to significantly slow down planet migration, and test the possibility by using a hybrid numerical integrator to simulate the disk-planet interaction. The key component of the scenario is the role of low viscosity regions in protostellar disks known as dead zones, which affect planetary migration in two ways. First of all, it allows a smaller-mass planet to open a gap, and hence switch the faster type I migration to the slower type II migration. Secondly, a low viscosity slows down type II migration itself, because type II migration is directly proportional to the viscosity. We present numerical simulations of planetary migration by using a hybrid symplectic integrator-gas dynamics code. Assuming that the disk viscosity parameter inside the dead zone is (alpha=1e-4-1e-5), we find that, when a low-mass planet (e.g. 1-10 Earth masses) migrates from outside the dead zone, its migration is stopped due to the mass accumulation inside the dead zone. When a low-mass planet migrates from inside the dead zone, it opens a gap and slows down its migration. A massive planet like Jupiter, on the other hand, opens a gap and slows down inside the dead zone, independent of its initial orbital radius. The final orbital radius of a Jupiter mass planet depends on the dead zone's viscosity. For the range of alpha's noted above, this can vary anywhere from 7 AU, to an orbital radius of 0.1 AU that is characteristic of the hot Jupiters.

In the core accretion model, gas giant formation is a race between growth and migration; for a core to become a jovian planet, it must accrete its envelope before it spirals into the host star. We use a multizone numerical model to extend our previous investigation of the "window of opportunity" for gas giant formation within a disk. When the collision cross-section enhancement due to core atmospheres is taken into account, we find that a broad range of protoplanetary disks posses such a window.

In the standard model of gas giant planet formation, a large solid core (~ 10 times the Earth's mass) forms first, then accretes its massive envelope (100 or more Earth masses) of gas. However, inward planet migration due to gravitational interaction with the proto-stellar gas disk poses a difficulty in this model. Core-sized bodies undergo rapid "Type I" migration; for typical parameters their migration timescale is much shorter than their accretion timescale. How, then, do growing cores avoid spiraling into the central star before they ever get the chance to become gas giants? Here, we present a simple model of core formation in a gas disk which is viscously evolving. As the disk dissipates, accretion and migration timescales eventually become comparable. If this happens while there is still enough gas left in the disk to supply a jovian atmosphere, then a window of opportunity for gas giant formation opens. We examine under what circumstances this happens, and thus, what predictions our model makes about the link between proto-stellar disk properties and the likelihood of forming giant planets.

Young planets interact with their parent gas disks through tidal torques. An imbalance between inner and outer torques causes bodies of mass $\ga 0.1$ Earth masses to lose angular momentum and migrate inward rapidly relative to the disk; this is known as ``Type I'' migration. However, protoplanets that grow to gas giant mass, O($10^2) M_\oplus$, open a gap in the disk and are subsequently constrained to migrate more slowly, locked into the disk's viscous evolution in what is called "Type II" migration. In a young planetary system, both Type I and Type II bodies likely coexist; if so, differential migration ought to result in close encounters when the former originate on orbits exterior to the latter. We investigate the resulting dynamics, using two different numerical approaches: an N-body code with dissipative forces added to simulate the effect of the gas disk, and a hybrid code which combines an N-body component with a 1-dimensional viscous disk model, treating planet-disk interactions in a more self-consistent manner. In both cases, we find that sub-gap-opening bodies have a high likelihood of being resonantly captured when they encounter a gap-opening body. A giant planet thus tends to act as a barrier in a protoplanetary disk, collecting smaller protoplanets outside of its orbit. Such behavior has two important implications for giant planet formation: First, for captured protoplanets it mitigates the problem of the migration timescale becoming shorter than the growth timescale. Secondly, it suggests one path to forming systems with multiple giant planets: Once the first has formed, it traps/accretes the future solid core of the second in an exterior mean-motion resonance, and so on. The most critical step in giant planet formation may thus be the formation of the very first one.

Many young stars reside within the central half-parsec from SgrA*, the supermassive black hole in the Galactic Center. The origin of these stars remains a puzzle. Recently, Hansen and Milosavljevic (2003, HM) have argued that an Intermediate-Mass Black Hole (IMBH) could have delivered the young stars to the immediate vicinity of SgrA*. Here we focus on the final stages of the HM scenario. Namely, we integrate numerically the orbits of stars which are initially bound to the IMBH, but are stripped from it by the tidal field of SgrA*. Our numerical algorithm is a symplectic integrator designed specifically for the problem at hand; however, we have checked our results with SYMBA, a version of the widely available SWIFT code. We find that the distribution of the post-inspiral orbital parameters is sensitive to the eccentricity of the inspiraling IMBH. If the IMBH is on a circular orbit, then the inclinations of numerically computed orbits relative to the inspiral plane are almost always smaller than 10 degrees, and therefore (a) the simulations are in good agreement with the observed motions of stars in a clockwise-moving stellar disc, (b) the simulations never reproduce the orbits of stars outside this disc, which include those in the second thick ring of stars and the randomly oriented unrelaxed orbits of some of the S-stars. If the IMBH's orbital eccentricity is e=0.6, then approximately half of the stars end up with orbital inclinations below 10 degrees, and another half have inclinations anywhere between 0 and 180 degrees; this is somewhat closer to what's observed. We also show that if IRS13 cluster is bound by an IMBH, as has been argued by Maillard et. al. 2004, then the same IMBH could not have delivered all of the young stars to their present location.

We present model calculations for the expected surface density of Ly-$\alpha$ emitting primeval galaxies (PGs) at high redshifts. We assume that elliptical galaxies and bulges of spiral galaxies (= spheroids) formed early in the universe and that the Ly-$\alpha$ emitting PGs are these spheroids during their first burst of star formation at high redshift. One of the main assumptions of the models is that the Ly-$\alpha$ bright phase of this first starburst in the spheroids is confined to a short period after its onset due to rapid formation of dust. The models do not only explain the failure of early surveys for Ly-$\alpha $ emitting PGs but are also consistent with the limits of new surveys (e.g. the Calar Alto Deep Imaging Survey - CADIS). At faint detection limits $S_{lim} \le 10^{-20}$W/m$^2$ the surface density of Ly-$\alpha$ emitters is expected to vary only weakly in the redshift range between $z=3$ and $z=6$ with values $> 10^3 / {deg}^2 / \Delta z = 0.1 $ reaching its maximum at $z_0 \approx 4$. At shallower detection limits, $S_{lim} \ge 3\times 10^{-20}$W/m$^2$ the surface density of high-z Ly-$\alpha$ emitters is expected to be a steep function of redshift and detection limit. This explains the low success in finding bright Ly-$\alpha$ galaxies at $z > 5$. We demonstrate how the observed surface densities of Ly-$\alpha $ emitting PGs derived from recent surveys constrain the parameters of our models. Finally, we discuss the possibility that two Ly-$\alpha$ bright phases occur in the formation process of galaxies: An initial -- primeval -- phase in which dust is virtually non-existing, and a secondary later phase in which strong galactic winds as observed in some Lyman break galaxies facilitate the escape of Ly-$\alpha $ photons after dust has already been formed.

The observed orbits of extrasolar planets suggest that many giant planets migrate a considerable distance towards their parent star as a result of interactions with the protoplanetary disk, and that some of these planets become trapped in eccentricity-exciting mean motion resonances with one another during this migration. Using three-dimensional numerical simulations, we find that as long as the timescale for damping of the planets' eccentricities by the disk is close to or longer than the disk-induced migration timescale, and the outer planet is more than half the mass of the inner, resonant inclination excitation will also occur. Neither the addition of a (simple, fixed) disk potential, nor the introduction of a massive inner planetary system, inhibit entry into the inclination resonance. Therefore, such a mechanism may not be uncommon in the early evolution of a planetary system, and a significant fraction of exoplanetary systems may turn out to be non-coplanar.

Runaway growth ends when the largest protoplanets dominate the dynamics of the planetesimal disk; the subsequent self-limiting accretion mode is referred to as ``oligarchic growth.'' Here, we begin by expanding on the existing analytic model of the oligarchic growth regime. From this, we derive global estimates of the planet formation rate throughout a protoplanetary disk. We find that a relatively high-mass protoplanetary disk ($\sim$ 10$\times$ minimum-mass) is required to produce giant planet core-sized bodies ($\sim$ 10 M$_{\oplus}$) within the lifetime of the nebular gas ($\la$ 10 million years). However, an implausibly massive disk is needed to produce even an Earth mass at the orbit of Uranus by 10 Myrs. Subsequent accretion without the dissipational effect of gas is even slower and less efficient. In the limit of non-interacting planetesimals, a reasonable-mass disk is unable to produce bodies the size of the Solar System's two outer giant planets at their current locations on {\it any} timescale; if collisional damping of planetesimal random velocities is sufficiently effective, though, it may be possible for a Uranus/Neptune to form in situ in less than the age of the Solar System. We perform numerical simulations of oligarchic growth with gas, and find that protoplanet growth rates agree reasonably well with the analytic model as long as protoplanet masses are well below their estimated final masses. However, accretion stalls earlier than predicted, so that the largest final protoplanet masses are smaller than those given by the model. Thus the oligarchic growth model, in the form developed here, appears to provide an upper limit for the efficiency of giant planet formation.

Galaxies at high redshift with a strong Ly-Alpha emission line trace massive star formation in the absence of dust, and can therefore be regarded as a prime signature of the first major starburst in galaxies. We report results of the Ly-Alpha search within the Calar Alto Deep Imaging Survey (CADIS). With imaging Fabry-Perot interferometer CADIS can detect emission lines in three waveband windows free of night-sky emission lines at 700nm, 820nm, and 920nm. The typical flux detection limit for Ly-Alpha emission redshifted into these windows, Flim > 3X10^(-20) Wm^(-2), corresponds to (unobscured) star formation rates of >10Msun/yr at z=6. Candidate Ly-Alpha-emitting galaxies are selected from the total emission line sample, which contains more than 97% of objects at z<1.2, by the absence of flux below the Lyman limit (B-band "dropouts"), and the non-detection of secondary emission lines in narrow band filters. We have detected 5 bright Ly-Alpha-emitting galaxy candidates at z ~ 4.8, and 11 candidates at z ~ 5.7. For two of four observed Ly-Alpha candidates, one candidate at z ~ 4.8, and the other at z ~ 5.7, the emission line detected with the Fabry-Perot has been verified spectroscopically at the VLT. When compared to Ly-Alpha surveys at z<3.5 even the upper limits set by our list of candidates show that bright Ly-Alpha galaxies are significantly rarer at z>5 than the assumption of a non-evolving population would predict. Therefore we conclude that the Ly-Alpha bright phase of primeval star formation episodes reached its peak at redshifts between 3 and 6.

A planetary system may undergo significant radial rearrangement during the early part of its lifetime. Planet migration can come about through interaction with the surrounding planetesimal disk and the gas disk--while the latter is still present--as well as through planet-planet interactions. We review the major proposed migration mechanisms in the context of the planet formation process, in our Solar System as well as in others.

In the nucleated instability picture of gas giant formation, the final stage is the rapid accretion of a massive gas envelope by a solid core, bringing about a tenfold or more increase in mass. This tends to trigger the scattering of any nearby bodies, including other would-be giant planet cores; it has been shown in past work that the typical outcome is an outer planetary system very similar to our own. Here, we show that the gravitational scattering accompanying the formation of gas giant planets can also produce, in some cases, outer planets with semimajor axes much larger than those in the Solar System, and eccentricities which remain high for tens of millions of years. Rings, gaps and asymmetries detected in a number of circumstellar dust disks, which hint at the presence of embedded planets at stellocentric distances far beyond where planet formation is expected to occur, may be connected to such a scenario.

Orbits of known extrasolar planets that are located outside the tidal circularization regions of their parent stars are often substantially eccentric. By contrast, planetary orbits in our Solar System are approximately circular, reflecting planet formation within a nearly axisymmetric, circumsolar disk. We propose that orbital eccentricities may be generated by divergent orbital migration of two planets in a viscously accreting circumstellar disk. The migration is divergent in the sense that the ratio of the orbital period of the outer planet to that of the inner planet grows. As the period ratio diverges, the planets traverse, but are not captured into, a series of mean-motion resonances that amplify their orbital eccentricities in rough inverse proportion to their masses. Strong viscosity gradients in protoplanetary disks offer a way to reconcile the circular orbits of Solar System gas giants with the eccentric orbits of currently known extrasolar planets.

The outer giant planets, Uranus and Neptune, pose a challenge to theories of planet formation. They exist in a region of the Solar System where long dynamical timescales and a low primordial density of material would have conspired to make the formation of such large bodies ($\sim$ 15 and 17 times as massive as the Earth, respectively) very difficult. Previously, we proposed a model which addresses this problem: Instead of forming in the trans-Saturnian region, Uranus and Neptune underwent most of their growth among proto-Jupiter and -Saturn, were scattered outward when Jupiter acquired its massive gas envelope, and subsequently evolved toward their present orbits. We present the results of additional numerical simulations, which further demonstrate that the model readily produces analogues to our Solar System for a wide range of initial conditions. We also find that this mechanism may partly account for the high orbital inclinations observed in the Kuiper belt.

We present a weak shear analysis of the Abell 901/902 supercluster, composed of three rich clusters at z=0.16. Using a deep R-band image from the 0.5 x 0.5 degree MPG/ESO Wide Field Imager together with supplementary B-band observations, we build up a comprehensive picture of the light and mass distributions in this region. We find that, on average, the light from the early-type galaxies traces the dark matter fairly well, although one cluster is a notable exception to this rule. The clusters themselves exhibit a range of mass-to-light (M/L) ratios, X-ray properties, and galaxy populations. We attempt to model the relation between the total mass and the light from the early-type galaxies with a simple scale-independent linear biasing model. We find M/L_B=130h for the early type galaxies with zero stochasticity, which, if taken at face value, would imply Omega_m < 0.1. However, this linear relation breaks down on small scales and on scales equivalent to the average cluster separation (approximately 1 Mpc), demonstrating that a single M/L ratio is not adequate to fully describe the mass-light relation in the supercluster. Rather, the scatter in M/L ratios observed for the clusters supports a model incorporating non-linear biasing or stochastic processes. Finally, there is a clear detection of filamentary structure connecting two of the clusters, seen in both the galaxy and dark matter distributions, and we discuss the effects of cluster-cluster and cluster-filament interactions as a means to reconcile the disparate descriptions of the supercluster.

We present K-band number counts for the faint galaxies in the Calar Alto Deep Imaging Survey (CADIS). We covered 4 CADIS fields, a total area of 0.2deg^2, in the broad band filters B, R and K. We detect about 4000 galaxies in the K-band images, with a completeness limit of K=19.75mag, and derive the K-band galaxy number counts in the range of 14.25 < K < 19.75mag. This is the largest medium deep K-band survey to date in this magnitude range. The B- and R-band number counts are also derived, down to completeness limits of B=24.75mag and R=23.25mag. The K-selected galaxies in this magnitude range are of particular interest, since some medium deep near-infrared surveys have identified breaks of both the slope of the K-band number counts and the mean B-K color at K=17\sim18mag. There is, however, a significant disagreement in the K-band number counts among the existing surveys. Our large near-infrared selected galaxy sample allows us to establish the presence of a clear break in the slope at K=17.0mag from dlogN/dm = 0.64 at brighter magnitudes to dlogN/dm = 0.36 at the fainter end. We construct no-evolution and passive evolution models, and find that the passive evolution model can simultaneously fit the B-, R- and K-band number counts well. The B-K colors show a clear trend to bluer colors for K > 18mag. We also find that most of the K=18-20mag galaxies have a B-K color bluer than the prediction of a no-evolution model for an L_* Sbc galaxy, implying either significant evolution, even for massive galaxies, or the existence of an extra population of small galaxies.

Gravitational lensing magnifies the luminosity of galaxies behind the lens. We use this effect to constrain the total mass in the cluster Abell 1689 by comparing the lensed luminosities of background galaxies with the luminosity function of an undistorted field. Since galaxies are assumed to be a random sampling of luminosity space, this method is not limited by clustering noise. We use photometric redshift information to estimate galaxy distance and intrinsic luminosity. Knowing the redshift distribution of the background population allows us to lift the mass/background degeneracy common to lensing analysis. In this paper we use 9 filters observed over 12 hours with the Calar Alto 3.5m telescope to determine the redshifts of 1000 galaxies in the field of Abell 1689. Using a complete sample of 151 background galaxies we measure the cluster mass profile. We find that the total projected mass interior to 0.25h^(-1)Mpc is (0.48 +/- 0.16) * 10^(15)h^(-1) solar masses, where our error budget includes uncertainties from the photometric redshift determination, the uncertainty in the off-set calibration and finite sampling. This result is in good agreement with that found by number count and shear-based methods and provides a new and independent method to determine cluster masses.

We use a multi-color classification method introduced by Wolf, Meisenheimer & Roeser (2000) to reliably identify stars, galaxies and quasars in the up to 16-dimensional color space provided by the filter set of the Calar Alto Deep Imaging Survey (CADIS). The samples of stars, galaxies and quasars obtained this way have been used for dedicated studies published in separate papers. The classification is good enough to detect quasars rather completely and efficiently without confirmative spectroscopy. The multi-color redshifts are accurate enough for most statistical applications, e.g. evolutionary studies of the galaxy luminosity function. We characterize our current dataset on the CADIS 1h-, 9h- and 16h-fields. Using Monte-Carlo simulations we model the classification performance expected for CADIS. We present a summary of the classification results and discuss unclassified objects. More than 99% of the whole catalog sample at R<22 (more than 95% at R<23) are successfully classified matching the expectations derived from the simulations. A small number of peculiar objects challenging the classification are discussed in detail. Spectroscopic observations are used to check the reliability of the multi-color classification (6 mistakes among 151 objects with R<24). We also determine the accuracy of the multi-color redshifts which are rather good for galaxies (sigma_z = 0.03) and useful for quasars. We find the classification performance derived from the simulations to compare well with results from the real survey. Finally, we locate areas for potential improvement of the classification.

We present initial results from a field survey for extremely red objects (EROs, defined here as (R-K') >= 6 mag) covering 154 square arcminutes of sky, from the first of 7 deep, wide-field K' images obtained as part of the Calar Alto Deep Imaging Survey (CADIS). The 5-sigma point source detection limits are K'=20.5 and R=25.0, while extended-source limits are up to 0.5-0.75 mag brighter. We identify a total of 8 bright EROs with K' =< 19.0. Six of these bright EROs are resolved and are likely to be galaxies, while the remaining 2 are unresolved, with colors consistent with their being low-mass galactic stars. We derive a surface density for the 6 bright, extragalactic EROs of 0.039 +/- 0.016 per sq. arcmin, higher by a factor of 4 than previous values. We estimate that the volume density of bright EROs to be as high as that of nearby Seyfert galaxies.

Due to the invention of new powerful instruments in the recent past (e.g. 10m class telescopes) high redshift galaxies are no longer a curiosity. High redshift young star forming galaxies can be effectively discriminated from the much more abundant foreground galaxies by their special spectral properties: a prominent break at the Lyman limit (i.e. a complete absence of flux at wavelength below the Lyman limit), an intrinsically flat spectrum at wavelength long-ward of the Lyman limit and in the very early phase of evolution a strong Ly-alpha emission line with high equivalent width. In the last couple of years several hundred star forming galaxies with 2.5 < z < 3.5 could be identified using deep broad band images which identify the Lyman break (Steidel et al. 1996, 1998). Spectroscopic followup observations confirmed their high redshifts. I summarize the main properties of these Lyman break galaxies (LBGs). Furthermore, the very recent discoveries of very high-z strong Ly-alpha emitting galaxies at z>5 and new large systematic surveys (the Calar Alto Deep Imaging Survey and the survey of Hu et al. 1998) for such objects, suggest that the cosmic star formation rate in strong Ly-alpha emitters does not decrease at redshifts z > 3.5 as suggested for the Lyman break galaxy sample. I discuss the galaxies with the highest redshifts we know today and give an overview over the survey of Hu et al. and the Calar Alto Deep Imaging Survey.

We present the results of the monitoring of the GRB 970111 field that started 19 hours after the event. This observation represents the fastest ground-based follow-up performed for GRB 970111 in all wavelengths. As soon as the detection of the possible GRB 970111 X-ray afterglow was reported by Feroci et al. (1998) we reanalyzed the optical data collected for the GRB 970111 field. Although we detect small magnitude variability in some objects, no convincing optical counterpart is found inside the WFC error box. Any change in brightness 19 hours after the GRB is less than 0.2 mag for objects with B < 21 and R < 20.8. The bluest object found in the field is coincident with 1SAXJ1528.8+1937. Spectroscopic observations revealed that this object is a Seyfert-1 galaxy with redshift z=0.657, which we propose as the optical counterpart of the X-ray source. Further observations allowed to perform multicolour photometry for objects in the GRB 970111 error box. The colour-colour diagrams do not show any object with unusual colours. We applied a photometric classification method to the objects inside the GRB error box, that can distinguish stars from galaxies and estimate redshifts. We were able to estimate photometric redshifts in the range 0.2 < z < 1.4 for several galaxies in this field and we did not find any conspicuous unusual object. We note that GRB 970111 and GRB 980329 could belong to the same class of GRBs, which may be related to nearby sources (z ~1) in which high intrinsic absorption leads to faint optical afterglows.