We study the growth of stellar discs of Milky Way-sized galaxies using a suite of cosmological simulations. We calculate the half-mass axis lengths and axis ratios of stellar populations split by age in isolated galaxies with stellar mass $M_* = 10^7 - 10^{10} M_{\odot}$ at redshifts $z$ > 1.5. We find that in our simulations stars always form in relatively thin discs, and at ages below 100 Myr are contained within half-mass height $z_{1/2}$ ~ 0.1 kpc and short-to-long axis ratio $z_{1/2}/x_{1/2}$ ~ 0.15. The disc thickness increases with the age of stellar population, reaching median $z_{1/2}$ ~ 0.8 kpc and $z_{1/2}/x_{1/2}$ ~ 0.6 for stars older than 500 Myr. We trace the same group of stars over the simulation snapshots and show explicitly that their intrinsic shape grows more spheroidal over time. The observed increase in disc thickness is due to the gradual vertical expansion of the originally thin stellar population, combined with the rapidly changing orientation of the galactic plane which mixes the configuration of young stars. The frequently mentioned "upside-down" formation scenario of galactic discs, which posits that young stars form in already thick discs at high redshift, is describing only the final product of this quick disc inflation while missing the actual formation of stars within a fairly thin plane. The location of very young stars in thin discs is due to the correspondingly flat configuration of dense molecular gas that participates in star formation.

The nature of the Sombrero galaxy (M 104 = NGC 4594) has remained elusive despite many observational studies at a variety of wavelengths. Here we present Hubble Space Telescope imaging of two fields at $\sim$16 and 33 kpc along the minor axis to examine stellar metallicity gradients in the extended spheroid. We use this imaging, extending more than 2 mag below the tip of the red giant branch (TRGB), in combination with artificial star tests to forward model observed color-magnitude diagrams (CMDs), measuring metallicity distribution functions (MDFs) at different radii along the minor axis. An important and unexpected result is that the halo of the Sombrero is strikingly metal-rich: even the outer field, located at $\sim$17 effective radii of the bulge, has a median metallicity [Z/H]$\sim$-0.15 and the fraction of stars with [Z/H]<-1.0 is negligible. This is unprecedented among massive galaxy halos studied to date, even among giant ellipticals. We find significant radial metallicity gradients, characterized by an increase in the fraction of metal-poor stars with radius and a gradient in median metallicity of $\sim$-0.01 dex/kpc. The density profile is well fit by power laws with slopes that exhibit a dependence on metallicity, with flatter slopes for more metal-poor stars. We discuss our results in the context of recent stellar MDF studies of other nearby galaxies and potential formation scenarios for the Sombrero galaxy.

Sky survey telescopes and powerful targeted telescopes play complementary roles in astronomy. In order to investigate the nature and characteristics of the motions of very faint objects, a flexibly-pointed instrument capable of high astrometric accuracy is an ideal complement to current astrometric surveys and a unique tool for precision astrophysics. Such a space-based mission will push the frontier of precision astrometry from evidence of earth-massed habitable worlds around the nearest starts, and also into distant Milky way objects up to the Local Group of galaxies. As we enter the era of the James Webb Space Telescope and the new ground-based, adaptive-optics-enabled giant telescopes, by obtaining these high precision measurements on key objects that Gaia could not reach, a mission that focuses on high precision astrometry science can consolidate our theoretical understanding of the local universe, enable extrapolation of physical processes to remote redshifts, and derive a much more consistent picture of cosmological evolution and the likely fate of our cosmos. Already several missions have been proposed to address the science case of faint objects in motion using high precision astrometry ESA missions: NEAT for M3, micro-NEAT for S1 mission, and Theia for M4 and M5. Additional new mission configurations adapted with technological innovations could be envisioned to pursue accurate measurements of these extremely small motions. The goal of this white paper is to address the fundamental science questions that are at stake when we focus on the motions of faint sky objects and to briefly review quickly instrumentation and mission profiles.

We examine the nature of kpc-scale clumps seen in high-redshift galaxies using a suite of cosmological simulations of galaxy formation. We identify rest-frame UV clumps in mock HST images smoothed to 500 pc resolution, and compare them with the intrinsic 3D clumps of young stars identified in the simulations with 100 pc resolution. According to this comparison, we expect that the stellar masses of the observed clumps are overestimated by as much as an order of magnitude, and that the sizes of these clumps are also overestimated by factor of several, due to a combination of spatial resolution and projection. The masses of young stars contributing most of the UV emission can also be overestimated by factor of a few. We find that most clumps of young stars present in a simulation at one time dissolve on a timescale shorter than $\sim$150 Myr. Some clumps with dense cores can last longer but eventually disperse. Most of the clumps are not bound structures, with virial parameter $\alpha_{\rm vir}$ > 1. We find similar results for clumps identified in mock maps of H$\alpha$ emission measure. We examine the predictions for effective clump sizes from the linear theory of gravitational perturbations and conclude that they are inconsistent with being formed by global disc instabilities. Instead, the observed clumps represent random projections of multiple smaller star-forming regions.

In a series of three papers, we introduced a novel cluster formation model that describes the formation, growth, and disruption of star clusters in high-resolution cosmological simulations. We tested this model on a Milky Way-sized galaxy and found that various properties of young massive clusters, such as the mass function and formation efficiency, are consistent with observations in the local universe. Interestingly, most massive clusters -- globular cluster candidates -- are preferentially formed during major merger events. We follow the dynamical evolution of clusters in the galactic tidal field. Due to tidal disruption, the cluster mass function evolves from initial power law to a peaked shape. The surviving clusters at $z=0$ show a broad range of metallicity [Fe/H] from -3 to -0.5. A robust prediction of the model is the age--metallicity relation, in which metal-rich clusters are systematically younger than metal-poor clusters by up to 3 Gyr.

We observe two metal-poor main sequence stars that are members of the recently-discovered Sylgr stellar stream. We present radial velocities, stellar parameters, and abundances for 13 elements derived from high-resolution optical spectra collected using the Magellan Inamori Kyocera Echelle spectrograph. The two stars have identical compositions (within 0.13 dex or 1.2 sigma) among all elements detected. Both stars are very metal poor ([Fe/H] = -2.92 +/- 0.06). Neither star is highly enhanced in C ([C/Fe] < +1.0). Both stars are enhanced in the alpha elements Mg, Si, and Ca ([alpha/Fe] = +0.32 +/- 0.06), and ratios among Na, Al, and all Fe-group elements are typical for other stars in the halo and ultra-faint and dwarf spheroidal galaxies at this metallicity. Sr is mildly enhanced ([Sr/Fe] = +0.22 +/- 0.11), but Ba is not enhanced ([Ba/Fe] < -0.4), indicating that these stars do not contain high levels of neutron-capture elements. The Li abundances match those found in metal-poor unevolved field stars and globular clusters (log epsilon (Li) = 2.05 +/- 0.07), which implies that environment is not a dominant factor in determining the Li content of metal-poor stars. The chemical compositions of these two stars cannot distinguish whether the progenitor of the Sylgr stream was a dwarf galaxy or a globular cluster. If the progenitor was a dwarf galaxy, the stream may originate from a dense region such as a nuclear star cluster. If the progenitor was a globular cluster, it would be the most metal-poor globular cluster known.

Energetic feedback from star clusters plays a pivotal role in shaping the dynamical evolution of giant molecular clouds (GMCs). To study the effects of stellar feedback on the star formation efficiency of the clouds and the dynamical response of embedded star clusters, we perform a suite of isolated GMC simulations with star formation and momentum feedback subgrid models using the moving-mesh hydrodynamics code \textsc{Arepo}. The properties of our simulated GMCs span a wide range of initial mass, radius, and velocity configurations. We find that the ratio of the final stellar mass to the total cloud mass, $\epsilon_{\rm int}$, scales strongly with the initial cloud surface density and momentum feedback strength. This correlation is explained by an analytic model that considers force balancing between gravity and momentum feedback. For all simulated GMCs, the stellar density profiles are systematically steeper than that of the gas at the epochs of the peaks of star formation, suggesting a centrally concentrated stellar distribution. We also find that star clusters are always in a sub-virial state with a virial parameter $\sim0.6$ prior to gas expulsion. Both the sub-virial dynamical state and steeper stellar density profiles prevent clusters from dispersal during the gas removal phase of their evolution. The final cluster bound fraction is a continuously increasing function of $\epsilon_{\rm int}$. GMCs with star formation efficiency smaller than 0.5 are still able to form clusters with large bound fractions.

In previous papers of this series, we developed a new algorithm for modeling the formation of star clusters in galaxy formation simulations. Here we investigate how dissolution of bound star clusters affects the shape of the cluster mass function and the metallicity distribution of surviving clusters. Cluster evolution includes the loss of stars that become unbound due to tidal disruption as well as mass-loss due to stellar evolution. We calculate the tidal tensor along cluster trajectories and use it to estimate the instantaneous mass-loss rate. The typical tidal tensor exhibits large variations on a time-scale of $\sim100$~Myr, with maximum eigenvalue of $10^7$~Gyr$^{-2}$, and median value of $10^4$~Gyr$^{-2}$ for the first Gyr after cluster formation. As a result of dynamical disruption, at the final available output of our simulations at redshift $z\approx1.5$, the cluster mass function has an approximately log-normal shape peaked at $\sim10^{4.3}M_\odot$. Extrapolation of the disruption to $z=0$ results in too many low-mass clusters compared to the observed Galactic globular clusters (GCs). Over 70\% of GC candidates are completely disrupted before the present; only 10\% of the total GC candidate mass remains in surviving clusters. The total mass of surviving clusters at $z=0$ varies from run to run in the range $(2-6)\times10^7M_\odot$, consistent with the observed mass of GC systems in Milky Way-sized galaxies. The metallicity distributions of all massive star clusters and of the surviving GCs have similar shapes but different normalization because of cluster disruption. The model produces a larger fraction of very metal-poor clusters than observed. A robust prediction of the model is the age-metallicity relation, in which metal-rich clusters are systematically younger than metal-poor clusters by up to 3~Gyr.

Globular cluster (GC) systems demonstrate tight scaling relations with the properties of their host galaxies. In previous work, we developed an analytic model for GC formation in a cosmological context and showed that it matches nearly all of the observed scaling relations across 4 orders of magnitude in host galaxy mass. Motivated by the success of this model, we investigate in detail the physical origins and evolution of these scaling relations. The ratio of the combined mass in GCs $M_{\rm GC}$ to the host dark matter halo mass $M_h$ is nearly constant at all redshifts, but its normalization evolves by a factor of $\sim$10 from birth to $z=0$. The relation is steeper than linear at halo masses $M_h \lesssim 10^{11.5} M_{\odot}$, primarily due to non-linearity in the stellar mass-halo mass relation. The near constancy of the ratio $M_{\rm GC}/M_h$, combined with the shape of the stellar mass-halo mass relation, sets the characteristic $U-$shape of the GC specific frequency as a function of host galaxy mass. The contribution of accreted satellite galaxies to the buildup of GC systems is a strong function of the host galaxy mass, ranging from $\approx$0% at $M_h \approx 10^{11} M_{\odot}$ to 80% at $M_h \approx 10^{15} M_{\odot}$. The metal-poor clusters are significantly more likely to form ex-situ relative to the metal-rich clusters, but a substantial fraction of metal-poor clusters still form in-situ in lower mass galaxies. Similarly, the fraction of red clusters increases from $\approx 10$% at $M_h = 10^{11} M_{\odot}$ to $\approx 60$% at $M_h \approx 10^{13} M_{\odot}$, and flattens at higher $M_h$. Clusters formation occurs essentially continuously at high redshift, while at low redshift galactic mergers become increasingly important for cluster formation.

Gaia is currently revolutionizing modern astronomy. However, much of the Galactic plane, center and the spiral arm regions are obscured by interstellar extinction, rendering them inaccessible because Gaia is an optical instrument. An all-sky near infrared (NIR) space observatory operating in the optical NIR, separated in time from the original Gaia would provide microarcsecond NIR astrometry and millimag photometry to penetrate obscured regions unraveling the internal dynamics of the Galaxy.

Observations of young star clusters reveal that the high-mass end of the cluster initial mass function (CIMF) deviates from a pure power-law and instead truncates exponentially. We investigate the effects of this truncation on the formation of globular cluster (GC) systems by updating our analytic model for cluster formation and evolution, which is based on dark matter halo merger trees coupled to empirical galactic scaling relations, and has been shown in previous work to match a wide array of observational data. The cutoff masses of $M_c=10^{6.5} M_{\odot}$ or $10^{7}M_{\odot}$ match many scaling relations: between the GC system mass and host halo mass, between the average metallicity of the GC system and host halo mass, and the distribution of cluster masses. This range of $M_c$ agrees with indirect measurements from extragalactic GC systems. Models with $M_c<10^{6.5}M_{\odot}$ cannot reproduce the observed GC metallicity and mass distributions in massive galaxies. The slope of the mass-metallicity relation for metal-poor clusters (blue tilt) for all $M_c$ models is consistent with observations within their errors, when measured using the same method. We introduce an alternative, more robust fitting method, which reveals a trend of increasing tilt slope for lower $M_c$. In our model the blue tilt arises because the metal-poor clusters form in relatively low-mass galaxies which lack sufficient cold gas to sample the CIMF at highest masses. Massive blue clusters form in progressively more massive galaxies and inherit their higher metallicity. The metal-rich clusters do not exhibit such a tilt because they form in significantly more massive galaxies, which have enough cold gas to fully sample the CIMF.

We explore the properties of dense star clusters that are likely to be the nurseries of stellar black holes pairing in close binaries. We combine a cosmological model of globular cluster formation with analytic prescriptions for the dynamical assembly and evolution of black hole binaries to constrain which types of clusters are most likely to form binaries tight enough to coalesce within a Hubble time. We find that black hole binaries which are ejected and later merge ex-situ form in clusters of a characteristic mass $M_{\rm cl} \sim 10^{5.3}M_{\odot}$, whereas binaries which merge in-situ form in more massive clusters, $M_{\rm cl} \sim 10^{5.7}M_{\odot}$. The clusters which dominate the production of black hole binaries are similar in age and metallicity to the entire population. Finally, we estimate an approximate cosmic black hole merger rate of dynamically assembled binaries using the mean black hole mass for each cluster given its metallicity. We find an intrinsic rate of $\sim 6\,\mathrm{Gpc\,yr^{-1}}$ at $z=0$, a weakly increasing merger rate out to $z=1.5$, and then a decrease out to $z=4$. Our results can be used to provide a cosmological context and choose initial conditions in numerical studies of black hole binaries assembled in star clusters.

The Fermi Gamma-Ray Space Telescope has provided evidence for diffuse gamma-ray emission in the central parts of the Milky Way and the Andromeda galaxy. This excess has been interpreted either as dark matter annihilation emission or as emission from thousands of millisecond pulsars (MSPs). We have recently shown that old massive globular clusters may move towards the center of the Galaxy by dynamical friction and carry within them enough MSPs to account for the observed gamma-ray excess. In this paper we revisit the MSP scenario for the Andromeda galaxy, by modeling the formation and disruption of its globular cluster system. We find that our model predicts gamma-ray emission $\sim 2-3$ times larger than for the Milky Way, but still nearly an order of magnitude smaller than the observed Fermi excess in the Andromeda. Our MSP model can reproduce the observed excess only by assuming $\sim 8$ times larger number of old clusters than inferred from galaxy scaling relations. To explain the observations we require either that Andromeda deviates significantly from the scaling relations, or that a large part of its high-energy emission comes from additional sources.

We investigate the structure of galaxies formed in a suite of high-resolution cosmological simulations. Consistent with observations of high-redshift galaxies, our simulated galaxies show irregular, prolate shapes with thick stellar disks, which are dominated by turbulent motions instead of rotation. Yet molecular gas and young stars are restricted to relatively thin disks. We examine the accuracy of applying the Toomre linear stability analysis to predict the location and amount of gas available for star formation. We find that the Toomre criterion still works for these irregular galaxies, after correcting for multiple gas and stellar components: the $Q$ parameter in $\rm{H_2}$ rich regions is in the range $0.5-1$, remarkably close to unity. Due to the violent stellar feedback from supernovae and strong turbulent motions, young stars and molecular gas are not always spatially associated. Neither the $Q$ map nor the $\rm{H_2}$ surface density map coincide with recent star formation exactly. We argue that the Toomre criterion is a better indicator of future star formation than a single $\rm{H_2}$ surface density threshold because of the smaller dynamic range of $Q$. The depletion time of molecular gas is below 1~Gyr on kpc scale, but with large scatter. Centering the aperture on density peaks of gas/young stars systematically biases the depletion time to larger/smaller values and increases the scatter.

Globular cluster (GC) systems around galaxies of a vast mass range show remarkably simple scaling relations. The combined mass of all GCs is a constant fraction of the total galaxy mass and the mean metallicity and metallicity dispersion of the GC system scale up weakly with galaxy mass. The metallicity of massive, metal-poor ("blue") clusters increases with cluster mass, while that of metal-rich ("red") clusters does not. A significant age-metallicity relation emerges from analysis of resolved stellar populations in Galactic GCs and unresolved populations in nearby galaxies. Remarkably, all these trends can be explained by a simple merger-based model developed in previous work and updated here using recent observations of galaxy scaling relations at high redshift. We show that the increasing dispersion of GC metallicity distributions with galaxy mass is a robust prediction of the model. It arises from more massive galaxies having more mergers that combine satellite GC systems. The average metallicity also increases by 0.6~dex over 3~dex in halo mass. The models show a non-linear trend between the GC system mass and host galaxy mass which is consistent with the data. The model does not consider GC self-enrichment, yet predicts a correlation between cluster mass and metallicity for massive blue clusters. The age-metallicity relation is another robust prediction of the model. Half of all clusters are predicted to form within the redshift range $5<z<2.3$, corresponding to ages of $10.8-12.5$~Gyr, in halos of masses $10^{11}-10^{12.5} M_{\odot}$.

We investigate the possible connection between the most massive globular clusters, such as $\omega$ Cen and M54, and nuclear star clusters of dwarf galaxies that exhibit similar spreads in age and metallicity. We examine galactic nuclei in cosmological galaxy formation simulations at $z \approx 1.5$ to explore whether their age and metallicity spreads could explain these massive globular clusters. We derive structural properties of these nuclear regions, including mass, size, rotation, and shape. By using theoretical supernova yields to model the supernova enrichment in the simulations, we obtain individual elemental abundances for Fe, O, Na, Mg, and Al. Our nuclei are systematically more metal-rich than their host galaxies, which lie on the expected mass-metallicity relation. Some nuclei have a spread in Fe and age comparable to the massive globular clusters of the Milky Way, lending support to the hypothesis that nuclear star clusters of dwarf galaxies could be the progenitors of these objects. None of our nuclear regions contain the light element abundance spreads that characterize globular clusters, even when a large age spread is present. Our results demonstrate that extended star formation history within clusters, with metal pollution provided solely by supernova ejecta, is capable of replicating the metallicity spreads of massive globular clusters, but still requires another polluter to produce the light element variations.

The implementation of star formation and stellar feedback in cosmological simulations plays a critical role in shaping galaxy properties. In the first paper of the series, we presented a new method to model star formation as a collection of star clusters. In this paper, we improve the algorithm by eliminating accretion gaps, boosting momentum feedback, and introducing a subgrid initial bound fraction, $f_i$, that distinguishes cluster mass from stellar particle mass. We perform a suite of simulations with different star formation efficiency per freefall time $\epsilon_{\rm ff}$ and supernova momentum feedback intensity $f_{\rm boost}$. We find that the star formation history of a Milky Way-sized galaxy is sensitive to $f_{\rm boost}$, which allows us to constrain its value, $f_{\rm boost}\approx5$, in the current simulation setup. Changing $\epsilon_{\rm ff}$ from a few percent to 200\% has little effect on global galaxy properties. However, on smaller scales, the properties of star clusters are very sensitive to $\epsilon_{\rm ff}$. We find that $f_i$ increases with $\epsilon_{\rm ff}$ and cluster mass. Through the dependence on $f_i$, the shape of the cluster initial mass function varies strongly with $\epsilon_{\rm ff}$. The fraction of clustered star formation and maximum cluster mass increase with the star formation rate surface density, with the normalization of both relations dependent on $\epsilon_{\rm ff}$. The cluster formation timescale systematically decreases with increasing $\epsilon_{\rm ff}$. Local variations in the gas accretion history lead to a 0.25~dex scatter for the integral cluster formation efficiency. Joint constraints from all the observables prefer the runs that produce a median integral efficiency of 16\%.

We discuss some of the key open questions regarding the formation and evolution of globular clusters (GCs) during galaxy formation and assembly within a cosmological framework. The current state-of-the-art for both observations and simulations is described, and we briefly mention directions for future research. The oldest GCs have ages $\ge$ 12.5 Gyr and formed around the time of reionisation. Resolved colour-magnitude diagrams of Milky Way GCs and direct imaging of lensed proto-GCs at z $\sim$ 6 with JWST promise further insight. Globular clusters are known to host multiple populations of stars with variations in their chemical abundances. Recently, such multiple populations have been detected in $\sim$2 Gyr old compact, massive star clusters. This suggests a common, single pathway for the formation of GCs at high and low redshift. The shape of the initial mass function for GCs remains unknown, however for massive galaxies a power-law mass function is favoured. Significant progress has been made recently modelling GC formation in the context of galaxy formation, with success in reproducing many of the observed GC-galaxy scaling relations.

The Fermi Large Area Telescope has provided the most detailed view toward the Galactic Centre (GC) in high-energy gamma rays. Besides the interstellar emission and point-source contributions, the data suggest a residual diffuse gamma-ray excess. The similarity of its spatial distribution with the expected profile of dark matter has led to claims that this may be evidence for dark matter particle annihilation. Here, we investigate an alternative explanation that the signal originates from millisecond pulsars (MSPs) formed in dense globular clusters and deposited at the GC as a consequence of cluster inspiral and tidal disruption. We use a semi-analytical model to calculate the formation, migration, and disruption of globular clusters in the Galaxy. Our model reproduces the mass of the nuclear star cluster and the present-day radial and mass distribution of globular clusters. For the first time, we calculate the evolution of MSPs from disrupted globular clusters throughout the age of the Galaxy and consistently include the effect of the MSP spin-down due to magnetic-dipole breaking. The final gamma-ray amplitude and spatial distribution are in good agreement with the Fermi observations and provide a natural astrophysical explanation for the GC excess.

In the context of the ESA M5 (medium mission) call we proposed a new satellite mission, Theia, based on relative astrometry and extreme precision to study the motion of very faint objects in the Universe. Theia is primarily designed to study the local dark matter properties, the existence of Earth-like exoplanets in our nearest star systems and the physics of compact objects. Furthermore, about 15 $\%$ of the mission time was dedicated to an open observatory for the wider community to propose complementary science cases. With its unique metrology system and "point and stare" strategy, Theia's precision would have reached the sub micro-arcsecond level. This is about 1000 times better than ESA/Gaia's accuracy for the brightest objects and represents a factor 10-30 improvement for the faintest stars (depending on the exact observational program). In the version submitted to ESA, we proposed an optical (350-1000nm) on-axis TMA telescope. Due to ESA Technology readiness level, the camera's focal plane would have been made of CCD detectors but we anticipated an upgrade with CMOS detectors. Photometric measurements would have been performed during slew time and stabilisation phases needed for reaching the required astrometric precision.

Motivated by theoretical expectations that Nuclear Star Clusters (NSCs) in galactic centers may provide a favorable environment for super-massive black holes to form and/or efficiently grow, we set out to measure the fraction of nearby nucleated galaxies that also host an Active Galactic Nucleus (AGN). We targeted a distance-limited sample of 98 objects with the Chandra X-ray Telescope, down to a uniform X-ray luminosity threshold of $\sim$10$^{38}$ erg s$^{-1}$. The sample is composed of 47 late-types and 51 early-types, enabling us to further investigate the active fraction as a function of galactic morphology. After correcting for contamination to the nuclear X-ray signal from bright X-ray binaries, we measure an active fraction $f$=11.2$\%^{+7.4}_{-4.9}$ (1$\sigma$ C.L.) across the whole sample, in agreement with previous estimates based on an heterogeneous combination of optical, X-ray and radio diagnostics, by Seth et al. (2008). After accounting for the different stellar mass distributions in our samples, we find no statistically significant difference in the active fraction of early- vs. late-type nucleated galaxies, with $f$=10.6$\%^{+11.9}_{-4.9}$ and 10.8$\%^{+11.3}_{-6.3}$, respectively. For the early-type nucleated galaxies, we are able to carry out a controlled comparison with a parent sample of non-nucleated galaxies covering the same stellar mass range, finding again no statistically significant difference in the active fraction. Taken at face value, our findings suggest that the presence of a NSC does not facilitate nor enhance accretion-powered emission from a nuclear super-massive black hole. This is true even for late-type nucleated galaxies, home to bluer NSCs and arguably larger gas reservoirs.

We present new deep photometry of the rich globular cluster (GC) systems around the Brightest Cluster Galaxies UGC 9799 (Abell 2052) and UGC 10143 (Abell 2147), obtained with the HST ACS and WFC3 cameras. For comparison, we also present new reductions of similar HST/ACS data for the Coma supergiants NGC 4874 and 4889. All four of these galaxies have huge cluster populations (to the radial limits of our data, comprising from 12000 to 23000 clusters per galaxy). The metallicity distribution functions (MDFs) of the GCs can still be matched by a bimodal-Gaussian form where the metal-rich and metal-poor modes are separated by ~0.8 dex, but the internal dispersions of each mode are so large that the total MDF becomes very broad and nearly continuous from [Fe/H] = -2.4 to Solar. There are, however, significant differences between galaxies in the relative numbers of \emph{metal-rich} clusters, suggesting that they underwent significantly different histories of mergers with massive, gas-rich halos. Lastly, the proportion of metal-poor GCs rises especially rapidly outside projected radii R > 4 R_eff, suggesting the importance of accreted dwarf satellites in the outer halo. Comprehensive models for the formation of GCs as part of the hierarchical formation of their parent galaxies will be needed to trace the systematic change in structure of the MDF with galaxy mass, from the distinctly bimodal form in smaller galaxies up to the broad continuum that we see in the very largest systems.

We present a new implementation of star formation in cosmological simulations, by considering star clusters as a unit of star formation. Cluster particles grow in mass over several million years at the rate determined by local gas properties, with high time resolution. The particle growth is terminated by its own energy and momentum feedback on the interstellar medium. We test this implementation for Milky Way-sized galaxies at high redshift, by comparing the properties of model clusters with observations of young star clusters. We find that the cluster initial mass function is best described by a Schechter function rather than a single power law. In agreement with observations, at low masses the logarithmic slope is $\alpha\approx 1.8-2$, while the cutoff at high mass scales with the star formation rate. A related trend is a positive correlation between the surface density of star formation rate and fraction of stars contained in massive clusters. Both trends indicate that the formation of massive star clusters is preferred during bursts of star formation. These bursts are often associated with major merger events. We also find that the median timescale for cluster formation ranges from 0.5 to 4 Myr and decreases systematically with increasing star formation efficiency. Local variations in the gas density and cluster accretion rate naturally lead to the scatter of the overall formation efficiency by an order of magnitude, even when the instantaneous efficiency is kept constant. Comparison of the formation timescale with the observed age spread of young star clusters provides an additional important constraint on the modeling of star formation and feedback schemes.

We present new deep photometry of the globular cluster system (GCS) around NGC 6166, the central supergiant galaxy in Abell 2199. HST data from the ACS and WFC3 cameras in F475W, F814W are used to determine the spatial distribution of the GCS, its metallicity distribution function (MDF), and the dependence of the MDF on galactocentric radius and on GC luminosity. The MDF is extremely broad, with the classic red and blue subpopulations heavily overlapped, but a double-Gaussian model can still formally match the MDF closely. The spatial distribution follows a Sersic-like profile detectably to a projected radius of at least $R_{gc} = 250$ kpc. To that radius, the total number of clusters in the system is N_{GC} = 39000 +- 2000, the global specific frequency is S_N = 11.2 +- 0.6, and 57\% of the total are blue, metal-poor clusters. The GCS may fade smoothly into the Intra-Cluster Medium of A2199; we see no clear transition from the core of the galaxy to the cD halo or the ICM. The radial distribution, projected ellipticity, and mean metallicity of the red (metal-richer) clusters match the halo light extremely well for R > 15 kpc, both of them varying as \sigma_{MRGC} ~ \sigma_{light} ~ R^-1.8. By comparison, the blue (metal-poor) GC component has a much shallower falloff \sigma_{MPGC} ~ R^-1.0 and a more nearly spherical distribution. This strong difference in their density distributions produces a net metallicity gradient in the GCS as a whole that is primarily generated by the population gradient. With NGC 6166 we appear to be penetrating into a regime of high enough galaxy mass and rich enough environment that the bimodal two-phase description of GC formation is no longer as clear or effective as it has been in smaller galaxies.

We measure proper motions with the Hubble Space Telescope for 16 extreme radial velocity stars, mostly unbound B stars in the Milky Way halo. Twelve of these stars have proper motions statistically consistent with zero, and thus have radial trajectories statistically consistent with a Galactic center "hypervelocity star" origin. The trajectory of HE 0437-5439 is consistent with both Milky Way and Large Magellanic Cloud origins. A Galactic center origin is excluded at 3-sigma confidence for two of the lowest radial velocity stars in our sample, however. These two stars are probable disk runaways and provide evidence for ~500 km/s ejections from the disk. We also measure a significant proper motion for the unbound sdO star US 708. Its 1,000 km/s motion is in some tension with proposed supernova ejection models, but can be explained if US 708 was ejected from the stellar halo. In the future, we expect Gaia will better constrain the origin of these remarkable unbound stars.

The mass distribution and chemical composition of globular cluster (GC) systems preserve fossil record of the early stages of galaxy formation. The observed distribution of GC colors within massive early-type galaxies in the ACS Virgo Cluster Survey (ACSVCS) reveals a multi-modal shape, which likely corresponds to a multi-modal metallicity distribution. We present a simple model for the formation and disruption of GCs that aims to match the ACSVCS data. This model tests the hypothesis that GCs are formed during major mergers of gas-rich galaxies and inherit the metallicity of their hosts. To trace merger events, we use halo merger trees extracted from a large cosmological N-body simulation. We select 20 halos in the mass range of $2\times 10^{12}$ to $7\times 10^{13}M_\odot$ and match them to 19 Virgo galaxies with K-band luminosity between $3\times 10^{10}$ and $3\times 10^{11}L_\odot$. To set the [Fe/H] abundances, we use an empirical galaxy mass-metallicity relation. We find that a minimal merger ratio of 1:3 best matches the observed cluster metallicity distribution. A characteristic bimodal shape appears because metal-rich GCs are produced by late mergers between massive halos, while metal-poor GCs are produced by collective merger activities of less massive hosts at early times. The model outcome is robust to alternative prescriptions for cluster formation rate throughout cosmic time, but a gradual evolution of the mass-metallicity relation with redshift appears to be necessary to match the observed cluster metallicities. We also affirm the age-metallicity relation, predicted by an earlier model, in which metal-rich clusters are systematically several billion years younger than their metal-poor counterparts.

We present the first results from our HST Brightest Cluster Galaxy (BCG) survey of seven central supergiant cluster galaxies and their globular cluster (GC) systems. We measure a total of 48000 GCs in all seven galaxies, representing the largest single GC database. We find that a log-normal shape accurately matches the observed luminosity function (LF) of the GCs down to the GCLF turnover point, which is near our photometric limit. In addition, the LF has a virtually identical shape in all seven galaxies. Our data underscore the similarity in the formation mechanism of massive star clusters in diverse galactic environments. At the highest luminosities (log L > 10^7 L_Sun) we find small numbers of "superluminous" objects in five of the galaxies; their luminosity and color ranges are at least partly consistent with those of UCDs (Ultra-Compact Dwarfs). Lastly, we find preliminary evidence that in the outer halo (R > 20 kpc), the LF turnover point shows a weak dependence on projected distance, scaling as L_0 ~ R^-0.2, while the LF dispersion remains nearly constant.

We revisit the hypothesis that dense galactic nuclei are formed from inspiraling globular clusters. Recent advances in understanding of the continuous formation of globular clusters over cosmic time and the concurrent evolution of the galaxy stellar distribution allow us to construct a simple model that matches the observed spatial and mass distributions of clusters in the Galaxy and the giant elliptical galaxy M87. In order to compare with observations, we model the effects of dynamical friction and dynamical evolution, including stellar mass loss, tidal stripping of stars, and tidal disruption of clusters by the growing galactic nucleus. We find that inspiraling globular clusters form a dense central structure, with mass and radius comparable to the typical values in observed nuclear star clusters (NSCs) in late-type and low-mass early-type galaxies. The density contrast associated with the NSC is less pronounced in giant elliptical galaxies. Our results indicate that the NSC mass as a fraction of mass of the galaxy stellar spheroid scales as M_NSC/M_* = 0.0025 M_{*,11}^{-0.5}. Thus disrupted globular clusters could contribute most of the mass of NSCs in galaxies with stellar mass below 10^{11} Msun. The inner part of the accumulated cluster may seed the growth of a central black hole via stellar dynamical core collapse, thereby relieving the problem of how to form luminous quasars at high redshift. The seed black hole may reach about 10^5 Msun within 1 Gyr of the beginning of globular cluster formation.

We investigate the transition from primordial Population III (Pop III) star formation to normal Pop II star formation in the first galaxies using new cosmological hydrodynamic simulations. We find that while the first stars seed their host galaxies with metals, they cannot sustain significant outflows to enrich the intergalactic medium, even assuming a top-heavy initial mass function. This means that Pop III star formation could potentially continue until z~6 in different unenriched regions of the universe, before being ultimately shut off by cosmic reionization. Within an individual galaxy, the metal production and stellar feedback from Pop II stars overtake Pop III stars in 20-200 Myr, depending on galaxy mass.

We revisit the formation and evolution of the first galaxies using new hydrodynamic cosmological simulations with the ART code. Our simulations feature a recently developed model for H2 formation and dissociation, and a star formation recipe that is based on molecular rather than atomic gas. Here, we develop and implement a recipe for the formation of metal-free Population III stars in galaxy-scale simulations that resolve primordial clouds with sufficiently high density. We base our recipe on the results of prior zoom-in simulations that resolved the protostellar collapse in pre-galactic objects. We find the epoch during which Pop III stars dominated the energy and metal budget of the first galaxies to be short-lived. Galaxies which host Pop III stars do not retain dynamical signatures of their thermal and radiative feedback for more than 10^8 yr after the lives of the stars end in pair-instability supernovae, even when we consider the maximum reasonable efficiency of the feedback. Though metals ejected by the supernovae can travel well beyond the virial radius of the host galaxy, they typically begin to fall back quickly, and do not enrich a large fraction of the intergalactic medium. Galaxies with total mass in excess of 3 x 10^6 Msun re-accrete most of their baryons and transition to metal-enriched Pop II star formation.

Star formation plays a central role in the evolution of galaxies and of the Universe as a whole. Studies of star-forming regions in the local universe have shown that star formation typically occurs in a clustered fashion. Building a coherent picture of how star clusters form and evolve is therefore critical to our overall understanding of the star formation process. Most clusters disrupt after they form, thus contributing to the field star population. However, the most massive and dense clusters remain bound and survive for a Hubble time. These globular clusters provide unique observational probes of the formation history of their host galaxies. In particular, the age and metallicity can be determined for each globular cluster individually, allowing the distribution of ages and metallicities within host galaxies to be constrained. We show how space-based UV-to-near-IR imaging covering a wide field of view (>= 20 arcmin per axis) and deep UV/Optical multi-object spectroscopy of globular cluster systems in nearby galaxies would allow one to place important new constraints on the formation history of early-type galaxies and their structural subcomponents (e.g., bulge, halo).

The condensation of gas and stars in the inner regions of dark matter halos leads to a more concentrated dark matter distribution. While this effect is based on simple gravitational physics, the question of its validity in hierarchical galaxy formation has led to an active debate in the literature. We use a collection of several state-of-the-art cosmological hydrodynamic simulations to study the halo contraction effect in systems ranging from dwarf galaxies to clusters of galaxies, at high and low redshift. The simulations are run by different groups with different codes and include hierarchical merging, gas cooling, star formation, and stellar feedback. We show that in all our cases the inner dark matter density increases relative to the matching simulation without baryon dissipation, at least by a factor of several. The strength of the contraction effect varies from system to system and cannot be reduced to a simple prescription. We present a revised analytical model that describes the contracted mass profile to an rms accuracy of about 10%. The model can be used to effectively bracket the response of the dark matter halo to baryon dissipation. The halo contraction effect is real and must be included in modeling of the mass distribution of galaxies and galaxy clusters.

We study the detailed structure of galaxies at redshifts z > 2 using cosmological simulations with improved modeling of the interstellar medium and star formation. The simulations follow the formation and dissociation of molecular hydrogen, and include star formation only in cold molecular gas. The molecular gas is more concentrated towards the center of galaxies than the atomic gas, and as a consequence, the resulting stellar distribution is very compact. For halos with total mass above 10^{11} Mo, the median half-mass radius of the stellar disks is 0.8 kpc at z = 3. The vertical structure of the molecular disk is much thinner than that of the atomic neutral gas. Relative to the non-radiative run, the inner regions of the dark matter halo change shape from prolate to mildly oblate and align with the stellar disk. However, we do not find evidence for a significant dark disk of dark matter around the stellar disk. The outer halo regions retain the orientation acquired during accretion and mergers, and are significantly misaligned with the inner regions. The radial profile of the dark matter halo contracts in response to baryon dissipation, establishing an approximately isothermal profile throughout most of the halo. This effect can be accurately described by a modified model of halo contraction. The angular momentum of a fixed amount of inner dark matter is approximately conserved over time, while in the dissipationless case most of it is transferred outward during mergers. The conservation of the dark matter angular momentum provides supporting evidence for the validity of the halo contraction model in a hierarchical galaxy formation process.

A basic property of objects, like galaxies and halos that form in cosmological structure formation simulations, is their shape. Here, we critically investigate shape determination methods that are commonly used in the literature. It is found that using an enclosed integration volume and weight factors r^{-2} and r_{ell}^{-2} (elliptical radius) for the contribution of each particle or volume element in the shape tensor leads to biased axis ratios and smoothing of details when calculating the local shape as a function of distance from the center. To determine the local shape of matter distributions as a function of distance for well resolved objects (typically more than O(10^4) particles), we advocate a method that (1) uses an ellipsoidal shell (homoeoid) as an integration volume without any weight factors in the shape tensor and (2) removes subhalos.

Modern hydrodynamic simulations of galaxy formation are able to predict accurately the rates and locations of the assembly of giant molecular clouds in early galaxies. These clouds could host star clusters with the masses and sizes of real globular clusters. I describe current state-of-the-art simulations aimed at understanding the origin of the cluster mass function and metallicity distribution. Metallicity bimodality of globular cluster systems appears to be a natural outcome of hierarchical formation and gradually declining fraction of cold gas in galaxies. Globular cluster formation was most prominent at redshifts z>3, when massive star clusters may have contributed as much as 20% of all galactic star formation.

The Hypervelocity Star survey presents the currently largest sample of radial velocity measurements of halo stars out to 80 kpc. We apply spherical Jeans modeling to these data in order to derive the mass profile of the Galaxy. We restrict the analysis to distances larger than 25 kpc from the Galactic center, where the density profile of halo stars is well approximated by a single power law with logarithmic slope between -3.5 and -4.5. With this restriction, we also avoid the complication of modeling a flattened Galactic disk. In the range 25 < r < 80 kpc, the radial velocity dispersion declines remarkably little; a robust measure of its logarithmic slope is between -0.05 and -0.1. The circular velocity profile also declines remarkably little with radius. The allowed range of V_c(80 kpc) lies between 175 and 231 km/s, with the most likely value 193 km/s. Compared with the value at the solar location, the Galactic circular velocity declines by less than 20% over an order of magnitude in radius. Such a flat profile requires a massive and extended dark matter halo. The mass enclosed within 80 kpc is 6.9(+3.0-1.2) 10^11 Msun. Our sample of radial velocities is large enough that the biggest uncertainty in the mass is not statistical but systematic, dominated by the density slope and anisotropy of the tracer population. Further progress requires modeling observed datasets within realistic simulations of galaxy formation.

We use Hubble Space Telescope imaging to measure the absolute proper motion of the hypervelocity star (HVS) HE 0437-5439, a short-lived B star located in the direction of the Large Magellanic Cloud (LMC). We observe (\mu_\alpha, \mu_\delta)=(+0.53+-0.25(stat)+-0.33(sys), +0.09+-0.21(stat)+-0.48(sys)) mas/yr. The velocity vector points directly away from the center of the Milky Way; an origin from the center of the LMC is ruled out at the 3-sigma level. The flight time of the HVS from the Milky Way exceeds its main-sequence lifetime, thus its stellar nature requires it to be a blue straggler. The large space velocity rules out a Galactic-disk ejection. Combining the HVS's observed trajectory, stellar nature, and required initial velocity, we conclude that HE 0437-5439 was most likely a compact binary ejected by the Milky Way's central black hole.

Observed metallicities of globular clusters reflect physical conditions in the interstellar medium of their high-redshift host galaxies. Globular cluster systems in most large galaxies display bimodal color and metallicity distributions, which are often interpreted as indicating two distinct modes of cluster formation. The metal-rich and metal-poor clusters have systematically different locations and kinematics in their host galaxies. However, the red and blue clusters have similar internal properties, such as the masses, sizes, and ages. It is therefore interesting to explore whether both metal-rich and metal-poor clusters could form by a common mechanism and still be consistent with the bimodal distribution. We present such a model, which prescribes the formation of globular clusters semi-analytically using galaxy assembly history from cosmological simulations coupled with observed scaling relations for the amount and metallicity of cold gas available for star formation. We assume that massive star clusters form only during mergers of massive gas-rich galaxies and tune the model parameters to reproduce the observed distribution in the Galaxy. A wide, but not entire, range of model realizations produces metallicity distributions consistent with the data. We find that early mergers of smaller hosts create exclusively blue clusters, whereas subsequent mergers of more massive galaxies create both red and blue clusters. Thus bimodality arises naturally as the result of a small number of late massive merger events. This conclusion is not significantly affected by the large uncertainties in our knowledge of the stellar mass and cold gas mass in high-redshift galaxies. The fraction of galactic stellar mass locked in globular clusters declines from over 10% at z>3 to 0.1% at present.

The disruption of the M33 galaxy is evident from its extended gaseous structure. We present new data from the Galactic Arecibo L-Band Feed Array HI (GALFA-HI) Survey that show the full extent and detailed spatial and kinematic structure of M33's neutral hydrogen. Over 18% of the HI mass of M33 (M_HI(tot) =1.4 x 10^9 Msun) is found beyond the star forming disk as traced in the far-ultraviolet (FUV). The most distinct features are extended warps, an arc from the northern warp to the disk, diffuse gas surrounding the galaxy, and a southern cloud with a filament back to the galaxy. The features extend out to 22 kpc from the galaxy center (18 kpc from the edge of the FUV disk) and the gas is directly connected to M33's disk. Only five discrete clouds (i.e., gas not directly connected to M33 in position-velocity space) are catalogued in the vicinity of M33, and these clouds show similar properties to Galactic and M31 halo clouds. M33's gaseous features most likely originate from the tidal disruption of M33 by M31 1-3 Gyr ago as shown from an orbit analysis which results in a tidal radius < 15 kpc in the majority of M33's possible orbits. M33 is now beyond the disruptive gravitational influence of M31 and the gas appears to be returning to M33's disk and redistributing its star formation fuel. M33's high mean velocity dispersion in the disk (~18.5 km/s) may also be consistent with the previous interaction and high rate of star formation. M33 will either exhaust its star formation fuel in the next few Gyrs or eventually become star formation fuel for M31. The latter represents the accretion of a large gaseous satellite by a spiral galaxy, similar to the Magellanic Clouds' relationship to the Galaxy.

Globular cluster systems in most large galaxies display bimodal color and metallicity distributions, which are frequently interpreted as indicating two distinct modes of cluster formation. The metal-rich (red) and metal-poor (blue) clusters have systematically different locations and kinematics in their host galaxies. However, the red and blue clusters have similar internal properties, such as the masses, sizes, and ages. It is therefore interesting to explore whether both metal-rich and metal-poor clusters could form by a common mechanism and still be consistent with the bimodal distribution. We show that if all globular clusters form only during mergers of massive gas-rich protogalactic disks, their metallicity distribution could be statistically consistent with that of the Galactic globulars. We take galaxy assembly history from cosmological dark matter simulations and couple it with the observed scaling relations for the amount of cold gas available for star formation. In the best-fit model, early mergers of smaller hosts create exclusively blue clusters, whereas subsequent mergers of progenitor galaxies with a range of masses create both red and blue clusters. Thus bimodality arises naturally as the result of a small number of late massive merger events. We calculate the cluster mass loss, including the effects of two-body scattering and stellar evolution, and find that more blue clusters than red clusters are disrupted by the present time, because of their smaller initial masses and larger ages. The present-day mass function in the best-fit model is consistent with the Galactic distribution. However, the spatial distribution of model clusters is much more extended than observed and is independent of the parameters of our model.

Large extra dimensions have been proposed as a possible solution to the hierarchy problem in physics. One of the suggested models, the RS2 braneworld model, makes a prediction that black holes evaporate by Hawking radiation on a short timescale that depends on the black hole mass and on the asymptotic radius of curvature of the extra dimensions. Thus the size of the extra dimensions can be constrained by astrophysical observations. Here we point out that the black hole, recently discovered in a globular cluster in galaxy NGC 4472, places the strongest constraint on the maximum size of the extra dimensions, L < 0.003 mm. This black hole has the virtues of old age and relatively small mass. The derived upper limit is within an order of magnitude of the absolute limit afforded by astrophysical observations of black holes.

Great strides have been made in the last two decades in determining how galaxies evolve from their initial dark matter seeds to the complex structures we observe at z=0. The role of mergers has been documented through both observations and simulations, numerous satellites that may represent these initial dark matter seeds have been discovered in the Local Group, high redshift galaxies have been revealed with monstrous star formation rates, and the gaseous cosmic web has been mapped through absorption line experiments. Despite these efforts, the dark matter simulations that include baryons are still unable to accurately reproduce galaxies. One of the major problems is our incomplete understanding of how a galaxy accretes its baryons and subsequently forms stars. Galaxy formation simulations have been unable to accurately represent the required gas physics on cosmological timescales, and observations have only just begun to detect the star formation fuel over a range of redshifts and environments. How galaxies obtain gas and subsequently form stars is a major unsolved, yet tractable problem in contemporary extragalactic astrophysics. In this paper we outline how progress can be made in this area in the next decade.

The fact that dark matter (DM), thus far, has revealed itself only on scales of galaxies and larger, again thrusts onto astrophysics the opportunity and the responsibility to confront the age old mystery "What is the nature of matter?" By deriving basic data on the nature of DM - e.g., mass of its particle(s), present mean temperature, distribution in galaxies and other structures in the universe, and capacity for dissipational collapse - we will be uncovering the properties of the dominant species of matter in the universe and significantly extending the standard models of particle physics. Determining the mass of the DM particle to an order of magnitude would help to sort out the particle family to which it (or they) belongs. Beyond mass, there are issues of stability. The DM particle may be unstable with a measurable half-life, or it may become unstable after absorbing a certain amount of energy from collisions. In both cases it would contribute to the present hot dark matter component. Some key parameters of DM can most accurately be measured in the very nearby universe because DM dominates the mass in the outer Milky Way (MW), in other galaxies in the Local Group, and in the Local Group in its entirety. The presence and amount of DM can be quantified by study of dynamical processes observable in fine detail within these entities. Precise measurements of 3-D velocities for stars, coherent star streams, and stars in satellite stellar systems out to the edge of the Galaxy can reveal "what is the shape, orientation, density law, and lumpiness of the dark matter halo" as well as "what is the total mass of the Galaxy?"

The concordance Lambda Cold Dark Matter (Lambda-CDM) model for the formation of structure in the Universe, while remarkably successful at describing observations of structure on large scales, continues to be challenged by observations on galactic scales. Fortunately, CDM models and their various proposed alternatives make a rich variety of testable predictions that make the Local Group and our own Milky Way Galaxy key laboratories for exploring dark matter (DM) in this regime. However, some of the most definitive tests of local DM require microarcsecond astrometry of faint sources, an astrometric regime that is a unique niche of SIM Lite. This chapter explores the important and distinct contributions that can be made by SIM Lite in the exploration of galaxy dynamics and DM on galaxy scales and that have cosmological consequences. Key areas of potential SIM Lite exploration include (1) measuring the shape, orientation, density law, and lumpiness of the dark halo of the Milky Way and other nearby galaxies, (2) determining the orbits of Galactic satellites, which may be representatives of late infall from the hierarchical formation of the Milky Way, (3) ascertaining the distribution of angular momentum and orbital anisotropy of stars and globular clusters to the outer reaches of the Galactic halo, dynamical properties that hold clues to the early hierarchical formation of the Galaxy, (4) measuring the physical nature of DM by placing strong constraints on the phase space density in the cores of nearby dSph galaxies, and (5) reconstructing the dynamical history of the Local Group through the determination of orbits and masses of its constituent galaxies.

We examine the accretion and merger histories of central and satellite galaxies in a smoothed particle hydrodynamics (SPH) cosmological simulation that resolves galaxies down to 7e9 M[Sun]. Most friends-of-friends halos in the simulation have a distinct central galaxy, typically two to five times more massive than the most massive satellite. As expected, satellites have systematically higher assembly redshifts than central galaxies of the same baryonic mass, and satellites in more massive halos form earlier. However, contrary to the simplest expectations, satellite galaxies continue to accrete gas and convert it to stars; the gas accretion declines steadily over a period of 0.5-1 Gyr after the satellite halo merges with a larger parent halo. Satellites in a cluster mass halo eventually begin to lose baryonic mass. Since z=1, 27% of central galaxies (above 3e10 M[Sun]) and 22% of present-day satellite galaxies have merged with a smaller system above a 1:4 mass ratio; about half of the satellite mergers occurred after the galaxy became a satellite and half before. In effect, satellite galaxies can remain "central" objects of halo substructures, with continuing accretion and mergers, making the transition in assembly histories and physical properties a gradual one. Implementing such a gradual transformation in semi-analytic models would improve their agreement with the observed colour distributions of satellite galaxies in groups and with the observed colour dependence of galaxy clustering.

The paucity of observed dwarf galaxies in the Local Group relative to the abundance of predicted dark matter halos remains one of the greatest puzzles of the LCDM paradigm. Solving this puzzle now requires not only matching the numbers of objects but also understanding the details of their star formation histories. We present a summary of such histories derived from the HST data using the color-magnitude diagram fitting method. To reduce observational uncertainties, we condense the data into five cumulative parameters -- the fractions of stellar mass formed in the last 1, 2, 5, and 10 Gyr, and the mean stellar age. We interpret the new data with a phenomenological model based on the mass assembly histories of dark matter halos and the Schmidt law of star formation. The model correctly predicts the radial distribution of the dwarfs and the fractions of stars formed in the last 5 and 10 Gyr. However, in order to be consistent with the observations, the model requires a significant amount of recent star formation in the last 2 Gyr. Within the framework of our model, this prolonged star formation can be achieved by adding a stochastic variation in the density threshold of the star formation law. The model results are not sensitive to late gas accretion, the slope of the Schmidt law, or the details of cosmic reionization. A few discrepancies still remain: our model typically predicts too large stellar masses, only a modest population of ultra-faint dwarfs, and a small number of dwarfs with anomalously young stellar populations. Nevertheless, the observed star formation histories of Local Group dwarfs are generally consistent the expected star formation in cold dark matter halos.

We review current state of neutron star cooling theory and discuss the prospects to constrain the equation of state, neutrino emission and superfluid properties of neutron star cores by comparing the cooling theory with observations of thermal radiation from isolated neutron stars.

We probe the evolution of globular clusters that could form in giant molecular clouds within high-redshift galaxies. Numerical simulations demonstrate that the large and dense enough gas clouds assemble naturally in current hierarchical models of galaxy formation. These clouds are enriched with heavy elements from earlier stars and could produce star clusters in a similar way to nearby molecular clouds. The masses and sizes of the model clusters are in excellent agreement with the observations of young massive clusters. Do these model clusters evolve into globular clusters that we see in our and external galaxies? In order to study their dynamical evolution, we calculate the orbits of model clusters using the outputs of the cosmological simulation of a Milky Way-sized galaxy. We find that at present the orbits are isotropic in the inner 50 kpc of the Galaxy and preferentially radial at larger distances. All clusters located outside 10 kpc from the center formed in the now-disrupted satellite galaxies. The spatial distribution of model clusters is spheroidal, with a power-law density profile consistent with observations. The combination of two-body scattering, tidal shocks, and stellar evolution results in the evolution of the cluster mass function from an initial power law to the observed log-normal distribution.

We present stellar velocity dispersion profiles for seven Milky Way dwarf spheroidal (dSph) satellite galaxies. We have measured 8394 line-of-sight velocities (+/- 2.5 km/s) for 6804 stars from high-resolution spectra obtained at the Magellan and MMT telescopes. We combine these new data with previously published velocities to obtain the largest available kinematic samples, which include more than 5500 dSph members. All the measured dSphs have stellar velocity dispersion of order 10 km/s that remains approximately constant with distance from the dSph center, out to and in some cases beyond the radius at which the mean surface brightness falls to the background level. Assuming dSphs reside within dark matter halos characterized by the NFW density profile, we obtain reasonable fits to the empirical velocity dispersion profiles. These fits imply that, among the seven dSphs, M_vir ~ 10^[8-9] M_sun. The mass enclosed at a radius of 600 pc, the region common to all data sets, ranges from (2-7) x 10^7 M_sun .

The host galaxies of the five local, z<0.25, long-duration gamma-ray bursts, each of which had a well-documented associated supernova, are all faint and metal-poor compared to the population of local star-forming galaxies. We quantify this statement by using a previous analysis of star-forming galaxies (z<0.2) from the Sloan Digital Sky Survey to estimate the fraction of local star formation as a function of host galaxy oxygen abundance. We find that only a small fraction (<25%) of current star formation occurs in galaxies with oxygen abundance 12+log(O/H)<8.6, i.e., about half that of the Milky Way. However, all five low-z GRB hosts have oxygen abundance below this limit, in three cases very significantly so. If GRBs traced local star formation independent of metallicity, the probability of obtaining such low abundances for all five hosts would be P~0.1%. We conclude that GRBs trace only low-metallicity star formation, and that the Milky Way has been too metal rich to host long GRBs for at least the last several billion years. This result has implications for the potential role of GRBs in mass extinctions, for searches for recent burst remnants in the Milky Way and other large galaxies, for non-detections of late radio emission from local core-collapse supernovae, and for the production of cosmic rays in the local Universe. We also find that the isotropic energy release of these five GRBs, E_iso, steeply decreases with increasing host oxygen abundance, suggesting an upper metallicity limit for ``cosmological'' GRBs at ~0.15 Z_solar. (Abridged)

The concordance cosmological model based on cold dark matter makes definitive predictions for the growth of galaxies in the Universe, which are being actively studied using numerical simulations. These predictions appear to contradict the observations of dwarf galaxies. Dwarf dark matter halos are more numerous and have steeper central density profiles than the observed galaxies. The first of these small-scale problems, the "missing satellites problem", can be resolved by accounting for the low efficiency of gas cooling and star formation in dwarf halos. A newly-discovered class of HyperVelocity Stars will soon allow us to test another generic prediction of CDM models, the triaxial shapes of dark matter halos. Measuring the proper motions of HVS will probe the gravitational potential out to 100 kpc and will constrain the axis ratios and the orientation of the Galactic halo.

If the dark matter particle is a neutralino then the first structures to form are cuspy cold dark matter (CDM) haloes collapsing after redshifts z ~ 100 in the mass range 10^{-6} - 10^{-3} Msun. We carry out a detailed study of the survival of these micro-haloes in the Galaxy as they experience tidal encounters with stars, molecular clouds, and other dark matter substructures. We test the validity of analytic impulsive heating calculations using high resolution N-body simulations. A major limitation of analytic estimates is that mean energy inputs are compared to mean binding energies, instead of the actual mass lost from the system. This energy criterion leads to an overestimate of the stripped mass and underestimate of the disruption timescale since CDM haloes are strongly bound in their inner parts. We show that a significant fraction of material from CDM micro-haloes can be unbound by encounters with Galactic substructure and stars, however the cuspy central regions remain relatively intact. Furthermore, the micro-haloes near the solar radius are those which collapse significantly earlier than average and will suffer very little mass loss. Thus we expect a fraction of surviving bound micro-haloes, a smooth component with narrow features in phase space, which may be uncovered by direct detection experiments, as well as numerous surviving cuspy cores with proper motions of arc-minutes per year, which can be detected indirectly via their annihilation into gamma-rays.

We present five new satellites of the Milky Way discovered in Sloan Digital Sky Survey (SDSS) imaging data, four of which were followed-up with either the Subaru or the Isaac Newton Telescopes. They include four probable new dwarf galaxies -- one each in the constellations of Coma Berenices, Canes Venatici, Leo and Hercules -- together with one unusually extended globular cluster, Segue 1. We provide distances, absolute magnitudes, half-light radii and color-magnitude diagrams for all five satellites. The morphological features of the color-magnitude diagrams are generally well described by the ridge line of the old, metal-poor globular cluster M92. In the last two years, a total of ten new Milky Way satellites with effective surface brightness mu_v >~ 28 mag/sq. arcsec have been discovered in SDSS data. They are less luminous, more irregular and appear to be more metal-poor than the previously-known nine Milky Way dwarf spheroidals. The relationship between these objects and other populations is discussed. We note that there is a paucity of objects with half-light radii between ~40 pc and ~ 100 pc. We conjecture that this may represent the division between star clusters and dwarf galaxies.

We test the hypothesis that metal-poor globular clusters form within disk galaxies at redshifts z>3. Numerical simulations demonstrate that giant gas clouds, which are cold and dense enough to produce massive star clusters, assemble naturally in hierarchical models of galaxy formation at high redshift. Do model clusters evolve into observed globular clusters or are they disrupted before present as a result of the dynamical evolution? To address this question, we calculate the orbits of model clusters in the time-variable gravitational potential of a Milky Way-sized galaxy, using the outputs of a cosmological N-body simulation. We find that at present the orbits are isotropic in the inner 50 kpc of the Galaxy and preferentially radial at larger distances. All clusters located outside 10 kpc from the center formed in satellite galaxies, some of which are now tidally disrupted and some of which survive as dwarf galaxies. The spatial distribution of model clusters is spheroidal and the fit to the density profile has a power-law slope of 2.7. The combination of two-body relaxation, tidal shocks, and stellar evolution drives the evolution of the cluster mass function from an initial power law to a peaked distribution, in agreement with observations. However, not all initial conditions and not all evolution scenarios are consistent with the observed mass function of the Galactic globular clusters. The successful models require the average cluster density to be constant initially for clusters of all mass and to remain constant with time. Synchronous formation of all clusters at a single epoch (z=4) and continuous formation over a span of 1.6 Gyr (between z=9 and z=3) are both consistent with the data. For both formation scenarios, we provide online catalogs of the main physical properties of model clusters.

We investigate structural properties of dark matter halos of disk galaxies in LCDM cosmology, using a well-defined sample of 81 disk-dominated galaxies from the SDSS redshift survey. We model the mass-velocity and fundamental plane relations of these galaxies, which are constructed from the galaxy stellar mass, disk scale length, and optical Halpha rotation velocity at 2.2 scale lengths. We calculate a sequence of model galaxy populations, defined by the distribution of the stellar disk-to-total mass fraction, m_d. We include the effect of adiabatic contraction of dark matter halos in response to condensation of baryons. We find that models with constant m_d underpredict the intrinsic scatter of the TF and FP relations and predict an (unobserved) strong correlation between TF residuals. Introducing a scatter of disk mass fractions and allowing the mean value m_d to scale with the stellar surface density significantly improves observational match of both the slope and intercept of the TF relation and reduces the predicted residual correlation enough to be consistent with the data. Our best-fit models with a Kroupa stellar IMF over-produce the galaxy stellar mass function and predict the virial r-band mass-to-light ratios, M_vir/L_r, systematically lower than those inferred from galaxy-galaxy weak lensing and satellite dynamics. We investigate three possible solutions to these problems: (1) ignoring the effects of adiabatic contraction, (2) adopting a ``light'' stellar IMF, or (3) considering the lower halo concentrations predicted for a low cosmological power spectrum normalization. Any of these solutions yields acceptable residual correlations and relieves most of the observational tension between the TF relation and the galaxy stellar mass function (abridged).

In this Letter, we study a localized stellar overdensity in the constellation of Ursa Major, first identified in Sloan Digital Sky Survey (SDSS) data and subsequently followed up with Subaru imaging. Its color-magnitude diagram (CMD) shows a well-defined sub-giant branch, main sequence and turn-off, from which we estimate a distance of ~30 kpc and a projected size of ~250 x 125 pc. The CMD suggests a composite population with some range in metallicity and/or age. Based on its extent and stellar population, we argue that this is a previously unknown satellite galaxy of the Milky Way, hereby named Ursa Major II (UMa II) after its constellation. Using SDSS data, we find an absolute magnitude of M_V \~ -3.8, which would make it the faintest known satellite galaxy. UMa II's isophotes are irregular and distorted with evidence for multiple concentrations; this suggests that the satellite is in the process of disruption.

We test the hypothesis that globular clusters form in supergiant molecular clouds within high-redshift galaxies. Numerical simulations demonstrate that such large, dense, and cold gas clouds assemble naturally in current hierarchical models of galaxy formation. These clouds are enriched with heavy elements from earlier stars and could produce star clusters in a similar way to nearby molecular clouds. The masses and sizes of the model clusters are in excellent agreement with the observations of young massive clusters. Do these model clusters evolve into globular clusters that we see in our and external galaxies? In order to study their dynamical evolution, we calculate the orbits of model clusters using the outputs of the cosmological simulation of a Milky Way-sized galaxy. We find that at present the orbits are isotropic in the inner 50 kpc of the Galaxy and preferentially radial at larger distances. All clusters located outside 10 kpc from the center formed in the now-disrupted satellite galaxies. The spatial distribution of model clusters is spheroidal, with a power-law density profile consistent with observations. The combination of two-body scattering, tidal shocks, and stellar evolution results in the evolution of the cluster mass function from an initial power law to the observed log-normal distribution. However, not all initial conditions and not all evolution scenarios are consistent with the observed mass function.

We study thermal structure and evolution of magnetars as cooling neutron stars with a phenomenological heat source in a spherical internal layer. We explore the location of this layer as well as the heating rate that could explain high observable thermal luminosities of magnetars and would be consistent with the energy budget of neutron stars. We conclude that the heat source should be located in an outer magnetar's crust, at densities rho < 5e11 g/cm^3, and should have the heat intensity of the order of 1e20 erg/s/cm^3. Otherwise the heat energy is mainly emitted by neutrinos and cannot warm up the surface.

We use measurements of the stellar velocity dispersion profile of the Fornax dwarf spheroidal galaxy to derive constraints on its dark matter distribution. Though the data are unable to distinguish between models with small cores and those with cusps, we show that a large > 1 kpc dark matter core in Fornax is highly implausible. Irrespective of the origin of the core, reasonable dynamical limits on the mass of the Fornax halo constrain its core radius to be no larger than 700 pc. We derive an upper limit core radius of 300 pc by demanding that the central phase space density of Fornax not exceed that directly inferred from the rotation curves of low-mass spiral galaxies. Further, if the halo is composed of warm dark matter then phase-space constraints force the core to be quite small in order to avoid conservative limits from the Ly alpha forest power spectrum, implying a core radius < 85 pc. We discuss our results in the context of the idea that the extended globular cluster distribution in Fornax can be explained by the presence of a large 1.5 kpc core. A self-consistent core of this size would be drastically inconsistent with the expectations of standard warm or cold dark matter models, and would also require an unreasonably massive dark matter halo, with a maximum circular velocity of 200 km/s.

Precise proper motion measurements (sigma_mu ~ 10 mkas/yr) of the recently discovered hyper-velocity star (HVS) SDSS J090745.0+024507 would yield significant constraints on the axis ratios and orientation of a triaxial model for the Galactic halo. Triaxiality of dark matter halos is predicted by Cold Dark Matter models of galaxy formation and may be used to probe the nature of dark matter. However, unless the distance to this star is determined to better than 10%, these constraints suffer from one-dimensional degeneracies, which we quantify. We show how proper motion measurements of several HVSs could simultaneously resolve the distance degeneracies of all such stars and produce a detailed picture of the triaxial halo. Additional HVSs may be found from radial velocity surveys or from parallax/proper-motion data derived from GAIA. High-precision proper-motion measurements of these stars using the Space Interferometry Mission (SIM PlanetQuest) would substantially tighten the constraints they yield on the Galactic potential.

We report on the angular correlation function of Lyman-break galaxies (LBGs) at z~4 and 5 from deep samples obtained from the Great Observatories Deep Origins Survey (GOODS). Similar to LBGs at z~3, the shape of w(theta) of the GOODS LBGs is well approximated by a power-law with slope beta~0.6 at angular separation theta > 10 arcsec. The clustering strength of z~4, 5 LBGs also depends on the rest-frame UV luminosity, with brighter galaxies more strongly clustered than fainter ones, implying a general correlation between halos' mass and LBGs' star-formation rate. At smaller separations, w(theta) of deep samples significantly exceeds the extrapolation of the large-scale power-law fit, implying enhanced spatial clustering at scales r < 1 Mpc. We also find that bright LBGs statistically have more faint companions on scales theta < 20 arcsec than fainter ones, showing that the enhanced small-scale clustering is very likely due to sub-structure, namely the fact that massive halos can host multiple galaxies. A simple model for the halo occupation distribution and the CDM halo mass function reproduce well the observed w(theta). The scaling relationship of the clustering strength with volume density and with redshift is quantitatively consistent with that of CDM halos. A comparison of the clustering strength of three samples of equal luminosity limit at z ~ 3, 4 and 5 shows that the LBGs at z~5 are hosted in halos about one order of magnitude less massive than those in the lower redshift bins, suggesting that star-formation was more efficient at higher-redshift.

We study the "minimal" cooling scenario of superfluid neutron stars with nucleon cores, where the direct Urca process is forbidden and the enhanced cooling is produced by the neutrino emission due to Cooper pairing of neutrons. Extending our previous consideration (Gusakov et al. 2004a), we include the effects of accreted envelopes of light elements. We employ phenomenological density-dependent critical temperatures T_{cp}(\rho) and T_{cnt}(\rho) of singlet-state proton and triplet-state neutron pairing in a stellar core, as well as the critical temperature T_{cns}(\rho) of singlet-state neutron pairing in a stellar crust. We show that the presence of accreted envelopes simplifies the interpretation of observations of thermal radiation from isolated neutron stars in the scenario of Gusakov et al. (2004a) and widens the class of models for nucleon superfluidity in neutron star interiors consistent with the observations.

The cooling of baryons in the centers of dark matter halos leads to a more concentrated dark matter distribution. This effect has traditionally been calculated using the model of adiabatic contraction, which assumes spherical symmetry, while in hierarchical formation scenarios halos grow via multiple violent mergers. We test the adiabatic contraction model in high-resolution cosmological simulations and find that the dissipation of gas indeed increases the density of dark matter and steepens its radial profile compared to the case without cooling. Although the standard model systematically overpredicts the increase of dark matter density, a simple modification of the assumed invariant from M(r)r to M(<r>)r, where <r> is the orbit-averaged particle position, reproduces the simulated profiles within 10%.

We have derived a set of equations to describe the thermal evolution of a neutron star which undergoes small-amplitude radial pulsations. We have taken into account, in the frame of the General Theory of Relativity, the pulsation damping due to the bulk and shear viscosity and the accompanying heating of the star. The neutrino emission of a pulsating non-superfluid star and its heating due to the bulk viscosity are calculated assuming that both processes are determined by the non-equilibrium modified Urca process. Analytical and numerical solutions to the set of equations of the stellar evolution are obtained for linear and strongly non-linear deviations from beta-equilibrium. It is shown that a pulsating star may be heated to very high temperatures, while the pulsations damp very slowly with time (a power law damping for 100-1000 years), as long as the damping is determined by the bulk viscosity. The contribution of the shear viscosity to the damping becomes important in a rather cool star with a low pulsation energy.

We study the cooling of superfluid neutron stars whose cores consist of nucleon matter with the Akmal-Pandharipande-Ravenhall equation of state. This equation of state opens the powerful direct Urca process of neutrino emission in the interior of most massive neutron stars. Extending our previous studies (Gusakov et al. 2004a, Kaminker et al. 2005), we employ phenomenological density-dependent critical temperatures T_{cp}(\rho) of strong singlet-state proton pairing (with the maximum T_{cp}^{max} \sim 7e9 K in the outer stellar core) and T_{cnt}(\rho) of moderate triplet-state neutron pairing (with the maximum T_{cnt}^{max} \sim 6e8 K in the inner core). Choosing properly the position of T_{cnt}^{max} we can obtain a representative class of massive neutron stars whose cooling is intermediate between the cooling enhanced by the neutrino emission due to Cooper pairing of neutrons in the absence of the direct Urca process and the very fast cooling provided by the direct Urca process non-suppressed by superfluidity.

We present a study of the spatial distribution of subhalos in galactic dark matter halos using dissipationless cosmological simulations of the concordance LCDM model. We find that subhalos are distributed anisotropically and are preferentially located along the major axes of the triaxial mass distributions of their hosts. The Kolmogorov-Smirnov probability for drawing our simulated subhalo sample from an isotropic distribution is P_KS \simeq 1.5 \times 10^{-4}. An isotropic distribution of subhalos is thus not the correct null hypothesis for testing the CDM paradigm. The nearly planar distribution of observed Milky Way (MW) satellites is marginally consistent (probability \simeq 0.02) with being drawn randomly from the subhalo distribution in our simulations. Furthermore, if we select the subhalos likely to be luminous, we find a distribution that is consistent with the observed MW satellites. In fact, we show that subsamples of the subhalo population with a centrally-concentrated radial distribution, similar to that of the MW dwarfs, typically exhibit a comparable degree of planarity. We explore the origin of the observed subhalo anisotropy and conclude that it is likely due to (1) preferential accretion of subhalos along filaments, often closely aligned with the major axis of the host halo, and (2) evolution of satellite orbits within the prolate, triaxial potentials typical of CDM halos. Agreement between predictions and observations requires the major axis of the outer dark matter halo of the Milky Way to be nearly perpendicular to the disk. We discuss possible observational tests of such disk-halo alignment with current large galaxy surveys.

We discuss pycnonuclear burning of highly exotic atomic nuclei in deep crusts of neutron stars, at densities up to 1e13 g/cc. As an application, we consider pycnonuclear burning of matter accreted on a neutron star in a soft X-ray transient (SXT, a compact binary containing a neutron star and a low-mass companion). The energy released in this burning, while the matter sinks into the stellar crust under the weight of newly accreted material, is sufficient to warm up the star and initiate neutrino emission in its core. The surface thermal radiation of the star in quiescent states becomes dependent of poorly known equation of state (EOS) of supranuclear matter in the stellar core, which gives a method to explore this EOS. Four qualitatively different model EOSs are tested against observations of SXTs. They imply different levels of the enhancement of neutrino emission in massive neutron stars by (1) the direct Urca process in nucleon/hyperon matter; (2) pion condensates; (3) kaon condensates; (4) Cooper pairing of neutrons in nucleon matter with the forbidden direct Urca process. A low level of the thermal quiescent emission of two SXTs, SAX J1808.4-3658 and Cen X-4, contradicts model (4). Observations of SXTs test the same physics of dense matter as observations of thermal radiation from cooling isolated neutron stars, but the data on SXTs are currently more conclusive.

We study the formation of globular clusters (GCs) in a Milky Way-size galaxy using a high-resolution cosmological simulation. The clusters in our model form in the dense cores of supergiant molecular clouds in the gaseous disks of high-redshift galaxies. The properties of clusters are estimated using a physically-motivated subgrid model of the isothermal cloud collapse. The first clusters in the simulation form at z ~ 12, while we conjecture that the best conditions for GC formation appear to be at z ~ 3-5. Most clusters form in the progenitor galaxies of the virial mass >10^9 Msun and the total mass of the cluster population is strongly correlated with the mass of its host galaxy with a fraction ~2x10^-4 of the galactic baryons being in the form of GCs. In addition, the mass of the GC population and the maximum cluster mass in a given region strongly correlate with the local average star formation rate. We find that the mass, size, and metallicity distributions of the cluster population identified in the simulation are remarkably similar to the corresponding distributions of the Milky Way globulars. We find no clear mass-metallicity or age-metallicity correlations for the old clusters. The zero-age cluster mass function can be approximated by a power-law, dN/dM ~ M^-alpha, with alpha ~ 2, in agreement with the mass function of young stellar clusters in starbursting galaxies. However, the shape of the zero-age mass function may be better described by the high-mass tail of a lognormal distribution which peaks at ~10^3 Msun. We discuss in detail the origin, universality, and dynamical evolution of the globular cluster mass function. Our results indicate that globular clusters with properties similar to those of observed clusters can form naturally within young dense gaseous disks at z >~ 3 in the LCDM cosmology.

The cooling of gas in the centers of dark matter halos is expected to lead to a more concentrated dark matter distribution. The response of dark matter to the condensation of baryons is usually calculated using the model of adiabatic contraction, which assumes spherical symmetry and circular orbits. In contrast, halos in the hierarchical structure formation scenarios grow via multiple violent mergers and accretion along filaments, and particle orbits in the halos are highly eccentric. We study the effects of the cooling of gas in the inner regions of halos using high-resolution cosmological simulations which include gas dynamics, radiative cooling, and star formation. We find that the dissipation of gas indeed increases the density of dark matter and steepens its radial profile in the inner regions of halos compared to the case without cooling. For the first time, we test the adiabatic contraction model in cosmological simulations and find that the standard model systematically overpredicts the increase of dark matter density in the inner 5% of the virial radius. We show that the model can be improved by a simple modification of the assumed invariant from M(r)r to M(r_av)r, where r and r_av are the current and orbit-averaged particle positions. This modification approximately accounts for orbital eccentricities of particles and reproduces simulation profiles to within 10-20%. We present analytical fitting functions that accurately describe the transformation of the dark matter profile in the modified model and can be used for interpretation of observations.

We investigate the combined effect of neutron and proton superfluidities on the cooling of neutron stars whose cores consist of nucleons and electrons. We consider singlet-state pairing of protons and triplet-state pairing of neutrons in the cores of neutron stars. The critical superfluid temperatures T_c are assumed to depend on the density of matter. We study two types of neutron pairing with different components of the total angular momentum of Cooper pairs along the quantization axis (|m_J| =0 or 2). Our calculations are compared with observations of thermal emission from isolated neutron stars. We show that the observations can be interpreted by using two classes of superfluidity models: (1) strong proton superfluidity with a maximum critical temperature in the stellar core T_c^{max} > 4 \times 10^9 K and weak neutron superfluidity of any type (T_c^{max} < 2 \times 10^8 K); (2) strong neutron superfluidity (pairing with |m_J|=0) and weak proton superfluidity. The two types of models reflect an approximate symmetry with respect to an interchange of the critical temperatures of neutron and proton pairing.

The impact of nuclear physics theories on cooling of isolated neutron stars is analyzed. Physical properties of neutron star matter important for cooling are reviewed such as composition, the equation of state, superfluidity of various baryon species, neutrino emission mechanisms. Theoretical results are compared with observations of thermal radiation from neutron stars. Current constraints on theoretical models of dense matter, derived from such a comparison, are formulated.

During the APPLES parallel campaign, the HST Advanced Camera for Surveys has resolved a distant stellar system, which appears to be an isolated dwarf galaxy. It is characterized by a circularly symmetric distribution of stars with an integrated magnitude m(F775W) = 20.13 +- 0.02, a central surface brightness of ~ 21.33 +- 0.18 mag/arcsec^2 and a half-light radius of ~ 1.8 arcsec. The ACS and VLT spectra show no evidence of ionized gas and appear dominated by a 3 Gyr old stellar population. The OB spectral type derived for two resolved stars in the grism data and the systemic radial velocity of ~ 670 km/s measured from the VLT data give a fiducial distance of ~ 9 +- 2 Mpc. These findings, with the support of the spatial morphology, would classify the system among the dwarf spheroidal (dSph) galaxies. Following the IAU rules, we have named this newly discovered galaxy APPLES 1. An intriguing peculiarity of APPLES 1 is that the properties (age and metallicity) of the stellar content so far detected are similar to those of dSph galaxies in the Local Group, where star formation is thought to be driven by galaxy interactions and mergers. Yet, APPLES 1 seems not to be associated with a major group or cluster of galaxies. Therefore, APPLES 1 could be the first example of a field dSph galaxy with self-sustained and regulated star formation and, therefore, would make an interesting test case for studies of the formation and evolution of unperturbed dSph galaxies.

We describe a quasi-equilibrium profile of dark matter particles in the inner parsec of the Galaxy, rho ~ r^-3/2. This "mini-cusp" profile is caused by scattering with the dense stellar cluster around the supermassive black hole in Sgr A* and is independent of the initial conditions. The implications for detection of gamma rays from WIMP dark matter annihilation in the Galactic center are a mild enhancement of the flux, and a characteristic central feature in the angular distribution which could be detectable by high-resolution Atmospheric Cherenkov Telescopes.

Hierarchical Cold Dark Matter (CDM) models predict that Milky Way (MW) sized halos contain hundreds of dense low-mass dark satellites, an order of magnitude more than the number of observed satellites in the Local Group (LG). If the CDM paradigm is correct, we need to understand why most of these halos failed to form stars and become galaxies. We analyze the dynamical evolution of the satellite halos in a high-resolution cosmological simulation of MW sized halos in the LCDM cosmology. We find that about 10% of the substructure halos with the present masses <10^8-10^9 Msun (circular velocities Vmax<30 km/s) had considerably larger masses and circular velocities when they formed at redshifts z>2. After the initial period of mass accretion in isolation, these objects experience dramatic mass loss due to tidal stripping, in some cases even before they are accreted by their host halo. This can explain how the smallest dwarf spheroidal galaxies of the LG were able to build up a sizable stellar mass in their seemingly shallow potential wells. We propose a new model in which all of the luminous dwarf spheroidals in the Local Group are descendants of the relatively massive (>10^9 Msun) high-redshift systems, in which the gas could cool efficiently by atomic line emission and which were not significantly affected by the extragalactic ultraviolet radiation. We present a simple galaxy formation model based on the trajectories extracted from the simulation, which accounts for the bursts of star formation after strong tidal shocks and the inefficiency of gas cooling in halos with virial temperatures Tvir<~10^4 K. Our model reproduces the abundance, spatial distribution, and morphological segregation of the observed Galactic satellites. The results are insensitive to the redshift of reionization.

We show that the direct Urca process of neutrino emission is allowed in two possible phases of nonspherical nuclei (inverse cylinders and inverse spheres) in the mantle of a neutron star near the crust-core interface. The process is open because neutrons and protons move in a periodic potential created by inhomogeneous nuclear structures. In this way the nucleons acquire large quasimomenta needed to satisfy momentum-conservation in the neutrino reaction. The appropriate neutrino emissivity in a nonsuperfluid matter is about 2--3 orders of magnitude higher than the emissivity of the modified Urca process in the stellar core. The process may noticeably accelerate the cooling of low-mass neutron stars.

We simulate cooling of superfluid neutron stars with nucleon cores where direct Urca process is forbidden. We adopt density dependent critical temperatures $T_{cp}(\rho)$ and $T_{cn}(\rho)$ of singlet-state proton and triplet-state neutron pairing in a stellar core and consider a strong proton pairing (with maximum $T_{cp}^{max} \ga 5 \times 10^9$ K) and a moderate neutron pairing ($T_{cn}^{max} \sim 6 \times 10^8$ K). When the internal stellar temperature $T$ falls below $T_{cn}^{max}$, the neutrino luminosity $L_{CP}$ due to Cooper pairing of neutrons behaves $\propto T^8$, just as that produced by modified Urca process (in a non-superfluid star) but is higher by about two orders of magnitude. In this case the Cooper-pairing neutrino emission acts like an enhanced cooling agent. By tuning the density dependence $T_{cn}(\rho)$ we can explain observations of cooling isolated neutron stars in the scenario in which direct Urca process or similar process in kaon/pion condensed or quark matter are absent.

We calculate the thermal structure and quiescent thermal luminosity of accreting neutron stars (warmed by deep crustal heating in accreted matter) in soft X-ray transients (SXTs). We consider neutron stars with nucleon and hyperon cores and with accreted envelopes. It is assumed that an envelope has an outer helium layer (of variable depth) and deeper layers of heavier elements, either with iron or with much heavier nuclei (of atomic weight A > 100) on the top (Haensel & Zdunik 1990, 2003, astro-ph/0305220). The relation between the internal and surface stellar temperatures is obtained and fitted. The quiescent luminosity of the hottest (low-mass) and coldest (high-mass) neutron stars is calculated, together with the ranges of its possible variations due to variable thickness of the helium layer. The results are compared with observations of SXTs, particularly, containing the coldest (SAX J1808.4-3658) and the hottest (Aql X-1) neutron stars. The observations of SAX J1808.4-3658 in a quiescent state on March 24, 2001 (Campana et al. 2002, astro-ph/0206376) can be explained only if this SXT contains a massive neutron star with a nucleon/hyperon core; a hyperon core with a not too low fraction of electrons is preferable. Future observations may discriminate between the various models of hyperon/nucleon dense matter. The thermal emission of SAX J1808.4-3658 is also sensitive to the models of plasma ionization in the outermost surface layers and can serve for testing such models.

We study the thermal structure of neutron stars with magnetized envelopes composed of accreted material, using updated thermal conductivities of plasmas in quantizing magnetic fields, as well as equation of state and radiative opacities for partially ionized hydrogen in strong magnetic fields. The relation between the internal and local surface temperatures is calculated and fitted by an analytic function of the internal temperature, magnetic field strength, angle between the field lines and the normal to the surface, surface gravity, and the mass of the accreted material. The luminosity of a neutron star with a dipole magnetic field is calculated for various values of the accreted mass, internal temperature, and magnetic field strength. Using these results, we simulate cooling of superfluid neutron stars with magnetized accreted envelopes. We consider slow and fast cooling regimes, paying special attention to very slow cooling of low-mass superfluid neutron stars. In the latter case, the cooling is strongly affected by the combined effect of magnetized accreted envelopes and neutron superfluidity in the stellar crust. Our results are important for interpretation of observations of isolated neutron stars hottest for their age, such as RX J0822-43 and PSR B1055-52.

Numerical simulations of a Milky Way-size galaxy demonstrate that globular clusters with the properties similar to observed can form naturally at z > 3 in the concordance Lambda-CDM cosmology. The clusters in our model form in the strongly baryon-dominated cores of supergiant molecular clouds. The first clusters form at z = 12, while the peak formation appears to be at z = 3-5. The zero-age mass function of globular clusters can be approximated by a power-law dN/dM ~ M^-2, in agreement with observations of young massive star clusters.

The cooling theory of isolated neutron stars is reviewed. The main cooling regulators are discussed, first of all, operation of direct Urca process (or similar processes in exotic phases of dense matter) and superfluidity in stellar interiors. The prospects to constrain gross parameters of supranuclear matter in neutron-star interiors by confronting cooling theory with observations of isolated neutron stars are outlined. A related problem of thermal states of transiently accreting neutron stars with deep crustal heating of accreted matter is discussed in application to soft X-ray transients.

High-redshift clusters of galaxies show an over-abundance of spirals by a factor of 2-3, and the corresponding under-abundance of S0 galaxies, relative to the nearby clusters. This morphological evolution can be explained by tidal interactions with neighboring galaxies and with the hierarchically growing cluster halo. The efficiency of tidal interactions depends on the size and structure of the cluster, as well as on the epoch of its formation. I simulate the formation and evolution of Virgo-type clusters in three cosmologies: a critical density model Omega_0=1, an open model Omega_0=0.4, and a flat model Omega_0=0.4 with a cosmological constant. The orbits of identified halos are traced with a high temporal resolution (~10^7 yr). Halos with low relative velocities merge only shortly after entering the cluster; after virialization mergers are suppressed. The dynamical evolution of galaxies is determined by the tidal field along their trajectories. The maxima of the tidal force do not always correspond to closest approach to the cluster center. They are produced to a large extent by the local density structures, such as the massive galaxies and the unvirialized remnants of infalling groups of galaxies. Collisions of galaxies are intensified by the substructure, with about 10 encounters within 10 kpc per galaxy in the Hubble time. These very close encounters add an important amount (10-50%) of the total heating rate. The integrated effect of tidal interactions is insufficient to transform a spiral galaxy into an elliptical, but can produce an S0 galaxy. Overall, tidal heating is stronger in the low Omega_0 clusters.

Tidal forces acting on galaxies in clusters lead to a strong dynamical evolution. In order to quantify the amount of evolution, I run self-consistent N-body simulations of disk galaxies for a variety of models in the hierarchically forming clusters. The tidal field along the galactic orbits is extracted from the simulations of cluster formation in the Omega_0=1; Omega_0=0.4; and Omega_0=0.4, Omega_Lambda=0.6 cosmological scenarios. For large spiral galaxies with the rotation speed of 250 km/s, tidal interactions truncate massive dark matter halos at 30 +- 6 kpc, and thicken stellar disks by a factor 2 to 3, increasing Toomre's parameter to Q > 2 and halting star formation. Low density galaxies, such as the dwarf spheroidals with the circular velocity of 20 km/s and the extended low surface brightness galaxies with the scale length of 10-15 kpc, are completely disrupted by tidal shocks. Their debris contribute to the diffuse intracluster light. The tidal effects are significant not only in the core but throughout the cluster and can be parametrized by the critical tidal density. The tidally-induced evolution results in the transformation of the infalling spirals into S0 galaxies and in the depletion of the LSB population. In the low Omega_0 cosmological models, clusters form earlier and produce stronger evolution of galaxies.

Despite the rapid observational progress in the study of young massive star clusters, the formation of globular clusters still remains poorly understood. Yet, it is emerging that globular cluster formation is intimately linked to the formation of the Galaxy. I discuss a generic scenario of the cluster formation within progenitor galaxies, based on the available observational constraints. The oldest clusters formed around redshift z = 7, but the process continued at least until z = 3. Because of their high density, globular clusters survived when their progenitor galaxies were disrupted by the Galactic tidal field.

We extend our modeling of cooling of superfluid neutron stars by including the production of muons in the core, in addition to neutrons, protons, and electrons. The results are confronted with observations of middle-aged isolated NSs. Muons have little effect on the hydrostatic structure of NSs, on the slow cooling of low-mass NSs (RX J0822--43 and PSR 1055--52 in our model) and on the rapid cooling of massive NSs. They affect, however, the moderately fast cooling of medium-mass NSs (1E 1207--52, RX J0002+62, PSR 0656+14, Vela, and Geminga) and shift appreciably the mass range of these NSs to lower masses, which is important for correct interpretation of the observations. Moreover, the effects of muons can accurately be reproduced by a simple renormalization of NS models with no muons in the NS cores.

We consider cooling of neutron stars (NSs) with superfluid cores composed of neutrons, protons, and electrons (assuming singlet-state pairing of protons, and triplet-state pairing of neutrons). We mainly focus on (nonstandard) triplet-state pairing of neutrons with the $|m_J| = 2$ projection of the total angular momentum of Cooper pairs onto quantization axis. The specific feature of this pairing is that it leads to a power-law (nonexponential) reduction of the emissivity of the main neutrino processes by neutron superfluidity. For a wide range of neutron critical temperatures $T_{cn}$, the cooling of NSs with the $|m_J| = 2$ superfluidity is either the same as the cooling with the $m_J = 0$ superfluidity, considered in the majority of papers, or much faster. The cooling of NSs with density dependent critical temperatures $T_{cn}(\rho)$ and $T_{cp}(\rho)$ can be imitated by the cooling of the NSs with some effective critical temperatures $T_{cn}$ and $T_{cp}$ constant over NS cores. The hypothesis of strong neutron superfluidity with $|m_J| = 2$ is inconsistent with current observations of thermal emission from NSs, but the hypothesis of weak neutron superfluidity of any type does not contradict to observations.

The upper limit of the effective surface temperature of the neutron star (NS) PSR J0205+6449 in the supernova remnant 3C 58 obtained recently by Slane et al. (2002) is analyzed using a modern theory of NS cooling (Kaminker et al. 2002). The observational limit can be explained by cooling of a superfluid NS with the core composed of neutrons, protons and electrons, where direct Urca process is forbidden. However, combined with the data on the surface temperatures of other isolated NSs, it gives evidence (emphasized by Slane et al.) that direct Urca process is open in the inner cores of massive NSs. This evidence turns out to be less stringent than the evidence provided by the well known observations of Vela and Geminga.

Cooling of neutron stars (NSs) with the cores composed of neutrons, protons, and electrons is analyzed. The main cooling regulators are discussed: opening of direct Urca process in a NS central kernel; superfluidity of nucleons in NS interiors; surface layers of light (accreted) elements; strong surface magnetic fields. An emphasis is paid on the cooling scenario with strong $^1$S$_0$ pairing of protons and weak $^3$P$_2$ pairing of neutrons in the NS core, as well as strong $^1$S$_0$ pairing of neutrons in the NS crust. The theory predicts three types of isolated cooling middle-aged NSs with distinctly different properties: low-mass, slowly cooling NSs; medium-mass, moderately cooling NSs; massive, rapidly cooling NSs. The theory is compared with observations of isolated NSs -- pulsars and radio quiet NSs in supernova remnants. The constraints on physical properties of NSs which can be inferred from such a comparison are outlined.

Using current observational data and simple dynamical modeling, we demonstrate that Omega Cen is not special among the Galactic globular clusters in its ability to produce and retain the heavy elements dispersed in the AGB phase of stellar evolution. The multiple stellar populations observed in Omega Cen cannot be explained if it had formed as an isolated star cluster. The formation within a progenitor galaxy of the Milky Way is more likely, although the unique properties of Omega Cen still remain a mystery.

The observed rotation curves of dark matter-dominated dwarf galaxies indicate low density cores, contrary to the predictions of CDM models. A possible solution of this problem involves stellar feedback. A strong baryonic wind driven by vigorous star formation can remove a large fraction of the gas, causing the dark matter to expand. We explore the maximum effect of the feedback with an instantaneous removal of the gaseous disk. The energy input depends on the compactness of the disk, hence the specific angular momentum of the disk. For the plausible cosmological parameters and a large range of the disk angular momentum, the feedback is insufficient to destroy the central halo cusp, while the inner density is lowered only by a modest factor of 2 to 6. Any realistic modeling of the feedback would have even lesser impact on dark matter. We find that no star formation effect can resolve the problems of CDM cusps.

Cooling of neutron stars (NSs) with the cores composed of neutrons, protons, and electrons is simulated assuming $^1$S$_0$ pairing of neutrons in the NS crust, and also $^1$S$_0$ pairing of protons and weak $^3$P$_2$ pairing of neutrons in the NS core, and using realistic density profiles of the superfluid critical temperatures $T_{\rm c}(\rho)$. The theoretical cooling models of isolated middle-aged NSs can be divided into three main types. (I) {\it Low-mass}, {\it slowly cooling} NSs where the direct Urca process of neutrino emission is either forbidden or almost fully suppressed by the proton superfluidity. (II) {\it Medium-mass} NSs which show {\it moderate} cooling via the direct Urca process suppressed by the proton superfluidity. (III) {\it Massive} NSs which show {\it fast} cooling via the direct Urca process weakly suppressed by superfluidity. Confronting the theory with observations we treat RX J0822--43, PSR 1055--52 and RX J1856--3754 as slowly cooling NSs. To explain these sufficiently warm sources we need a density profile $T_{\rm c}(\rho)$ in the crust with a rather high and flat maximum and sharp wings. We treat 1E 1207--52, RX J0002+62, PSR 0656+14, Vela, and Geminga as moderately cooling NSs. We can determine their masses for a given model of proton superfluidity, $T_{\rm cp}(\rho)$, and the equation of state in the NS core. No rapidly cooling NS has been observed so far.

As shown recently by Kaminker et al. (2001), current observations of thermal emission of isolated middle-aged neutron stars (NSs) can be explained by cooling of NSs of different masses with the cores composed of neutrons, protons and electrons, assuming rather strong superfluidity (SF) of protons, weak triplet-state SF of neutrons and neglecting singlet-state SF of neutrons. We show that this explanation remains correct in the presence of singlet-state SF of neutrons in the NS crust and outermost core but under stringent constraints on the density profile of the SF critical temperature $T_{\rm cns}(\rho)$. In order to explain observations of (young and hot) RX J0822--43 and (old and warm) PSR 1055--52 and RX J1856--3754 as cooling not too massive NSs, the maximum $T_{\rm cns}^{\rm max}$ should be rather high ($\ga 5 \times 10^9$ K) and the decrease of $T_{\rm cns}(\rho)$ outside the maximum should be sharp. These results place important constraints on the models of nucleon SF in dense matter.

The existence of quasars at redshift z > 5 indicates that supermassive black holes were present since the very early times. If they grew by accretion, the seeds of mass ~ 10^5 Msun must have formed at z ~ 9. These seed black holes may result from the collapse and dissipation of primordial gas during the early stages of galaxy formation. I discuss the effects of magnetic fields on the fragmentation of cold gas clouds embedded into a hot diffuse phase and a virialized dark matter halo. The field of 10^-4 G ejected by supernova remnants can halt cloud break-up at 10^4 Msun. High star formation rates keep the clouds partially ionized, making ambipolar diffusion inefficient. The magnetically-supported clouds collapse into black holes, which later spiral via dynamical friction into a central cluster with the total mass Mbh ~ 6 10^6 Msun. As the cluster collapses, the black holes merge emitting gravitational radiation that should be detectable by LISA.

The dynamical evolution of galaxies in clusters is modified if dark matter is self-interacting. Heat conduction from the hot cluster halo leads to evaporation of the relatively cooler galactic halos. The stellar distribution would adiabatically expand as it readjusts to the loss of dark matter, reducing the velocity dispersion and increasing the half-light radius. If the dark matter content within that radius was f_dm = 25-50% of the total, as indicated by current observations, the ellipticals in clusters would be offset from the fundamental plane relation beyond the observational scatter. The requirement that their halos survive for a Hubble time appears to exclude just that range of the dark matter cross-section, 0.3 < sigma/m < 10^4 cm^2/g, thought to be optimal for reducing central halo cusps, unless f_dm < 15%. If the cross-section is allowed to vary with the relative velocity of dark matter particles, sigma \propto v^{-2 delta}, a new problem of evaporation of dark matter arises in the dwarf galaxies with low velocity dispersion. The halos of large galaxies in clusters and dwarf galaxies in the Local Group can both survive only if delta < 1.1 or delta > 1.8. In either case the problem of central density cusps remains.

Estimating the age of the Universe is an old problem. Rapid progress in observational cosmology in recent years has led to more accurate values of the fundamental parameters. The current most popular model is a flat Universe, with about 30% of the critical density in the form of matter (baryonic and non-baryonic) and 70% in the form of `dark energy'. These two densities, together with the Hubble constant (estimated to be about 70 km/s/Mpc), constrain the age of the Universe to be approximately 13 Gyr. This expansion age is uncomfortably close to the age of the oldest globular clusters (approximately 12.5 Gyr), in particular if they formed relatively recently. We review here proposed models of globular cluster formation and point out possible conflicts with cosmology.

Matter of subnuclear density in the inner crusts of neutron stars consists of neutron-rich atomic nuclei immersed in strongly degenerate relativistic gas of electrons and strongly nonideal liquid of neutrons. We give a brief account of thermodynamic and kinetic properties of this matter which are greatly affected by Coulomb and nuclear interactions (a companion paper on the outer envelopes is astro-ph/0012316) and show how they can be studied from observations of thermal radiation of young (age < 100 yr) neutron stars (for a discussion of the latter possibility, see astro-ph/0012306).

The internal properties of the neutron star crust can be probed by observing the epoch of thermal relaxation. After the supernova explosion, powerful neutrino emission quickly cools the stellar core, while the crust stays hot. The cooling wave then propagates through the crust, due to its finite thermal conductivity. When the cooling wave reaches the surface (age 10-100 yr), the effective temperature drops sharply from 250 eV to 30 or 100 eV, depending on the cooling model. The crust relaxation time is sensitive to the (poorly known) microscopic properties of matter of subnuclear density, such as the heat capacity, thermal conductivity, and superfluidity of free neutrons. We calculate the cooling models with the new values of the electron thermal conductivity in the inner crust, based on a realistic treatment of the shapes of atomic nuclei. Superfluid effects may shorten the relaxation time by a factor of 4. The comparison of theoretical cooling curves with observations provides a potentially powerful method of studying the properties of the neutron superfluid and highly unusual atomic nuclei in the inner crust.

We review the main neutrino emission mechanisms in neutron star crusts and cores. Among them are the well-known reactions such as the electron-positron annihilation, plasmon decay, neutrino bremsstrahlung of electrons colliding with atomic nuclei in the crust, as well as the Urca processes and neutrino bremsstrahlung in nucleon-nucleon collisions in the core. We emphasize recent theoretical achievements, for instance, band structure effects in neutrino emission due to scattering of electrons in Coulomb crystals of atomic nuclei. We consider the standard composition of matter (neutrons, protons, electrons, muons, hyperons) in the core, and also the case of exotic constituents such as the pion or kaon condensates and quark matter. We discuss the reduction of the neutrino emissivities by nucleon superfluidity, as well as the specific neutrino emission produced by Cooper pairing of the superfluid particles. We also analyze the effects of strong magnetic fields on some reactions, such as the direct Urca process and the neutrino synchrotron emission of electrons. The results are presented in the form convenient for practical use. We illustrate the effects of various neutrino reactions on the cooling of neutron stars. In particular, the neutrino emission in the crust is critical in setting the initial thermal relaxation between the core and the crust. Finally, we discuss the prospects of exploring the properties of supernuclear matter by confronting cooling simulations with observations of the thermal radiation from isolated neutron stars.

We solve the Boltzmann equation for cosmological neutrinos around the epoch of the electron-positron annihilation in order to verify the freeze-out approximation and to compute accurately the cosmological neutrino distribution function. We find the radiation energy density to be about 0.3% higher than the one predicted by the freeze-out approximation. As a result, the spectrum of the Cosmic Microwave Background anisotropies changes by 0.3-05%, depending on the angular scale, and the amplitude of the mass fluctuations on scales below about 100 h^{-1} Mpc decreases by about 0.2-0.3%.

We present new Fokker-Planck models of the evolution of globular clusters, including gravitational tidal shocks. We extend our calculations beyond the core collapse by adopting three-body binary heating. Effects of the shocks are included by adding the tidal shock diffusion coefficients to the ordinary Fokker-Planck equation: the first order heating term, <dE>, and the second order energy dispersion term, <dE^2>. As an example, we investigate the evolution of models for the globular cluster NGC 6254. Using the Hipparcos proper motions, we are now able to construct orbits of this cluster in the Galaxy. Tidal shocks accelerate significantly both core collapse and the evaporation of the cluster and shorten the destruction time from 24 Gyr to 18 Gyr. We examine various types of adiabatic corrections and find that they are critical for accurate calculation of the evolution. Without adiabatic corrections, the destruction time of the cluster is twice as short. We examine cluster evolution for a wide range of the concentration and tidal shock parameters, and determine the region of the parameter space where tidal shocks dominate the evolution. We present fitting formulae for the core collapse time and the destruction time, covering all reasonable initial conditions. In the limit of strong shocks, the typical value of the core collapse time decreases from 10 t_{rh} to 3 t_{rh} or less, while the destruction time is just twice that number. The effects of tidal shocks are rapidly self-limiting: as clusters lose mass and become more compact, the importance of the shocks diminishes. This implies that tidal shocks were more important in the past.

We derive expressions for the tidal field exerted by a spherically symmetric galaxy having an extended mass distribution, and use our analysis to calculate tidal perturbations and heating of stars in a globular cluster or a satellite galaxy orbiting in the external potential. We consider a Hernquist model and several power-law models for the primary. We integrate the tidal force over accurate orbits in the given potential by using the impulse approximation and by introducing adiabatic corrections that allow for the time variation of the force. We compare these results with the conventional straight path approximation. We also check the results with direct N-body simulations and find very good agreement. Heating on highly eccentric orbits dominates as the adiabatic corrections strongly damp the energy change on low eccentricity orbits. The results are illustrated for the example of globular cluster NGC 6712. For the orbital eccentricities higher than 0.9 the future lifetime of NGC 6712 is less than 10^{10} yr. Our analysis can be used to study tidal interactions of star clusters and galaxies in a semi-analytical manner.

We study the reaction of a globular star cluster to a time-varying tidal perturbation (gravitational shock) using self-consistent N-body simulations and address two questions. First, to what extent is the cluster interior protected by adiabatic invariants. Second, how much further energy change does the postshock evolution of the cluster potential produce and how much does it affect the dispersion of stellar energies. We introduce the adiabatic correction} as ratio of the energy change, <dE>, to its value in the impulse approximation. When the potential is kept fixed, the numerical results for the adiabatic correction for stars with orbital frequency \omega can be approximated as (1 + \omega^2 \tau^2)^{-\gamma}. For shocks with the characteristic duration of the order the half-mass dynamical time of the cluster, \tau < t_{dyn,h}, the exponent \gamma = 5/2. For more prolonged shocks, \tau > 4 t_{dyn,h}, the adiabatic correction is shallower, \gamma = 3/2. When we allow for self-gravity and potential oscillations which follow the shock, the energy of stars in the core changes significantly, while the total energy of the system is conserved. Paradoxically, the postshock potential fluctuations reduce the total amount of energy dispersion, <dE^2>. The effect is small but real and is due to the postshock energy change being statistically anti-correlated with the shock induced heating. These results are to be applied to Fokker-Planck models of the evolution of globular clusters.

Recent observations by the Hubble Space Telescope suggest that a strong dynamical and morphological evolution of disk galaxies has occurred in clusters of galaxies. I present the result of the self-consistent high-resolution simulations of galactc dynamics in the three cosmologically formed clusters. The disks of large galaxies survive but thicken by a factor of two and cease new star formation. The dark matter halos are truncated in agreement with analytical models. The overall effect is the morphological transformation of spirals into S0 galaxies.

The semi-analytic theory of tidal shocks proves to be a powerful tool to study tidal interactions of star clusters and satellite galaxies with their massive hosts. New models of the globular cluster evolution employ a combination of analytic estimates, solutions of the Fokker-Planck equation and direct N-body simulations. The models predict large destruction rates for the Galactic globular clusters. Those on the highly eccentric orbits around the Galactic center are much more likely to be disrupted than the ones on nearly circular orbits. The destruction rates are largely increased near the bulge. Disruption of the low-mass clusters changes the Luminosity Function of the Globular Cluster System, shifting the peak of the Luminosity Function to the brighter end.

Anticipating new spectroscopic studies of quasar companions, we calculate expected emission line spectra for the companion galaxies. Under the assumption of isotropic photoionizing radiation, we predict equivalent widths for Halpha, Hbeta, [O III], [O II], [Ne III], and Mg II lines. We consider 6 quasars, two of which (PKS 2135-14 and PKS 2349-014) have very close companions with projected separation < 4 h^{-1} kpc. For these companions the calculated emission lines should be easily observable, with W_lambda > 20 AA for the [O III], [O II], and Halpha lines. Presence of strong emission lines from the companions occupying only a tiny solid angle would suggest that 1) quasar radiation is isotropic; 2) companion galaxies have measurable amount of gas, which can be assessed using the photoionization models; 3) host galaxies, if present, have been depleted of gas or do not interact with the companions.

We have previously found in the Hubble Deep Field a significant angular correlation of faint, high color-redshift objects on scales below one arcsecond, or several kiloparsecs in metric size. We examine the correlation and nearest neighbor statistics to conclude that 38% of these objects in the HDF have a companion within one arcsecond, three times the number expected in a random distribution with the same number of objects. We examine three dynamical scenarios for these object multiplets: 1) the objects are star-forming regions within normal galaxies, whose disks have been relatively dimmed by K-correction and surface brightness dimming; 2) they are fragments merging into large galaxies; 3) they are satellites accreting onto normal L_* galaxies. We find that hypothesis 1 is most tenable. First, large galaxies in the process of a merger formation would have accumulated too much mass in their centers (5e12 M_sun inside 2 kpc) to correspond to present day objects. Second, accretion by dynamical friction occurs with a predictable density vs. radius slope, not seen among the faint HDF objects. Since the dynamical friction time is roughly (1 Gyr), a steady-state should have been reached by redshift z < 5. Star-forming regions within galaxies clearly present no dynamical problems. Since large spirals would still appear as such in the HDF, we favor a scenario in which the faint compact sources in the HDF are giant starforming regions within small normal galaxies, such as Magellanic irregulars. Finally we checked that reduction in mass-to-light from induced star-formation cannot alone explain the luminosity overdensity.

The dynamical evolution of the globular cluster systems in galaxies is predicted, based on the standard dynamical theory normalized to the example of the Milky Way. The major processes varying with the galactocentric distance are the tidal shocks and dynamical friction. Our simple model explains, on a quantitative basis, the observed differences of the inner and outer populations of globular clusters. We can thus calculate corrections for dynamical evolution for the luminosity function of globular clusters with the assumption that the initial luminosity function is identical in all galaxies (and we can test this assumption as well, in certain cases). Then we can compute the expected distribution of absolute magnitudes and compare it with the observed distribution of apparent magnitudes to estimate the distance moduli for M31 and M87. Using this new method we find dm(M31)=24.05 +- 0.23, dm(M87)=30.83 +- 0.17, as compared to current best estimates using other methods of dm(M31)=24.30 +- 0.20, dm(M87)=31.0 +- 0.1. As a check on the method we compute, and compare with observations, the differences between the inner and outer globular clusters in all three galaxies. This new method, coupled with HST observations, promises to provide an independent method of estimating distances to galaxies with recession velocities < 10,000 km/s, or D < 100 Mpc.

We consider the observational signature of the dynamical effects on the luminosity function of globular clusters. For the three best studied systems, in Milky Way, M31, and M87, there is a statistically significant difference between the inner and outer population of globular clusters. In all cases the inner clusters are on average brighter than the outer clusters ($0.26 < \Delta m_0 < 0.84$) and have a smaller dispersion in magnitudes ($0.04 < \Delta \sigma < 0.53$), with the larger differences for the local, better observed samples. The differences are of the type that would be expected if the inner population had been depleted by tidal shocks. The results suggest that the inner population suffers substantial evolution from its initial distribution and cannot not be used as a standard candle without correction for dynamical evolution.

We investigate the dynamical evolution of the Galactic Globular Cluster System in considerably greater detail than has been done hitherto, finding that destruction rates are significantly larger than given by previous estimates. More than half of the present clusters (52%-58% for the OC galactic model, and 75%-86% for the BSS model) will be destroyed in the next Hubble time. For the evolution of individual clusters we use a Fokker-Planck code including two-body relaxation, tidal truncation of clusters, compressive gravitational shocks while clusters pass through the Galactic disk, and tidal shocks due to passage close to the bulge. Gravitational shocks are treated comprehensively, using a recent result by Kundic & Ostriker (1995) that the shock-induced relaxation term, driving an additional dispersion of energies, is generally more important than the usual energy shift term. We discuss some implications for a past history of the Globular Cluster System, and the initial distribution of the destruction times raising the possibility that the current population is but a very small fraction of the initial population with the remnants of the destroyed clusters constituting presently a large fraction of the spheroid (bulge+halo) stellar population.

Spiral arms, if they are massive, exert gravitational torques that transport angular momentum radially within galactic disks. These torques depend not on the pattern speed or permanence of the arms but only on the nonaxisymmetric mass distribution. Hence the torques can be measured from photometry. We demonstrate this using $gri$ CCD data for M100 (NGC 4321). Since we find consistency among the three bands, we believe that dust and young stars in the arms do not seriously bias our results. If the present epoch is representative, the timescale for redistribution of angular momentum in M100 is $5-10$ Gyr, the main uncertainty being the mass-to-light ratio of the disk.