We present a detailed analysis of the behavior of the triaxial Schwarzschild orbit superposition method near the axisymmetric limit. Orbit superposition modeling is the primary method used to determine dynamical masses of supermassive black holes ($M_\mathrm{BH}$) in nearby galaxies; however, prior studies have reported conflicting results when comparing the outcome from axisymmetric orbit codes with that from a triaxial orbit code in the axisymmetric limit. We show that in order to achieve (oblate) axisymmetry in a triaxial code, care needs to be taken to axisymmetrize the short-axis tube orbits and to exclude both the long-axis tube and box orbits from the orbit library. Using up to 12 Gauss-Hermite moments of the line-of-sight velocity distributions as constraints, we demonstrate the effects of orbit types on the best-fit $M_\mathrm{BH}$ in orbit modeling of the massive elliptical galaxy NGC 1453 reported in Liepold et al. (2020). We identify a subset of slowly precessing quasi-planar orbits for which the typical integration times can be insufficient to fully capture the equilibrium orbital behavior in both axisymmetric and triaxial systems with central black holes. Further investigation is needed for a more reliable treatment of these orbits.

The inclination distribution of circumbinary planets (CBPs) is an important scientific issue. It is of great significance in estimating the occurrence rate of CBPs and studying their formation and evolution. Although the CBPs currently discovered by the transit method are nearly coplanar, the true inclination distribution of CBPs is still unknown. Previous researches on CBPs mostly regarded them as an isolated binary-planet system, without considering the birth environment of their host binaries. It is generally believed that almost all stars are born in clusters. Therefore, it is necessary to consider the impact of the close encounters of stars on the CBP systems. This article discusses how the close encounters of fly-by stars affect the inclination of CBP. Based on extensive numerical simulations, we found that CBPs in close binary with a spacing of about 0.2 au are almost unaffected by star fly-bys. Their orbits remain coplanar. However, when the spacing of the binary stars is greater than 1 au, the 2-3 fly-bys of the intruding star can excite a considerable inclination even for the CBP near the unstable boundary of the binary. For the planets in the outer region, a single star fly-by can excite an inclination to more than 5 degrees. Especially, CBPs in near polar or retrograde orbits can naturally form through bianry-star encounters. If close binaries are born in open clusters, our simulations suggest that there may be high-inclination CBPs in binaries with spacing > 1 au.

We present a new stellar dynamical mass measurement ($M_\mathrm{BH}$) of the supermassive black hole (SMBH) in NGC 1453, a fast-rotating massive elliptical galaxy in the MASSIVE survey. We measure stellar kinematics in 135 spatial bins in the central 1.5 kpc by 2 kpc region of the galaxy using high signal-to-noise ($S/N \sim 130$) spectra from the Gemini-North GMOS integral-field spectrograph (IFS). Combining with wide-field IFS kinematics out to $\sim 3$ effective radii and stellar light distributions from Hubble Space Telescope (HST) WFC3 images, we perform Schwarzschild orbit-based mass modeling in the axisymmetric limit to constrain the mass components in NGC 1453. The best-fit black hole mass is $M_\mathrm{BH} =(2.9 \pm 0.4) 10^9 M_\odot$; the mass models without a central black hole are excluded at the $8.7\sigma$ level. The NGC 1453 black hole lies within the intrinsic scatter of the SMBH and galaxy scaling relations, unlike three other galaxies hosting $\gtrsim 10^{10} M_\odot$ SMBHs in the MASSIVE sample. The high-$S/N$ GMOS spectra enable us to determine 8 moments of the Gauss-Hermite expansion of the line-of-sight velocity distributions (LOSVDs), which are used as constraints in the orbit modeling. The stellar orbits in the mass models are further constrained to produce negligible $h_9$ through $h_{12}$ to minimize spurious behavior in the LOSVDs. We show that truncating the series at $h_4$, as was often done in prior work, leads to a much weaker constraint on the inferred $M_\mathrm{BH}$ for NGC 1453. Furthermore, we discuss precautions and modifications that are needed to achieve axisymmetry in triaxial orbit codes that use the Schwarzschild method to sample the start space of stellar orbits in triaxial gravitational potentials.

We use high spatial resolution stellar velocity maps from the Gemini GMOS integral-field spectrograph (IFS) and wide-field velocity maps from the McDonald Mitchell IFS to study the stellar velocity profiles and kinematic misalignments from $\sim 200$ pc to $\sim 20$ kpc in 20 early-type galaxies with stellar mass $M_* > 10^{11.7} M_\odot$ in the MASSIVE survey. While 80% of the galaxies have low spins ($\lambda < 0.1$) and low rotational velocities ($< 50$ km/s) in both the central region and the main body, we find a diverse range of velocity features and misalignment angles. For the 18 galaxies with measurable central kinematic axes, 10 have well aligned kinematic axis and photometric major axis, and the other 8 galaxies have misalignment angles that are distributed quite evenly from $15^\circ$ to the maximal value of $90^\circ$. There is a strong correlation between central kinematic misalignment and galaxy spin, where all 4 galaxies with significant spins have well aligned kinematic and photometric axes, but only 43% of the low-spin galaxies are well aligned. The central and main-body kinematic axes within a galaxy are not always aligned. When the two kinematic axes are aligned ($\sim 60$% of the cases), they are either also aligned with the photometric major axis or orthogonal to it. We find 13 galaxies to also exhibit noticeable local kinematic twists, and one galaxy to have a counter-rotating core. A diverse assembly history consisting of multiple gas-poor mergers of a variety of progenitor mass ratios and orbits is likely to be needed to account for the predominance of low spins and the wide range of central and main-body velocity features reported here for local massive ETGs.

We identify multi-wavelength counterparts to 1,147 submillimeter sources from the S2COSMOS SCUBA-2 survey of the COSMOS field by employing a recently developed radio$+$machine-learning method trained on a large sample of ALMA-identified submillimeter galaxies (SMGs), including 260 SMGs identified in the AS2COSMOS pilot survey. In total, we identify 1,222 optical/near-infrared(NIR)/radio counterparts to the 897 S2COSMOS submillimeter sources with S$_{850}$>1.6mJy, yielding an overall identification rate of ($78\pm9$)%. We find that ($22\pm5$)% of S2COSMOS sources have multiple identified counterparts. We estimate that roughly 27% of these multiple counterparts within the same SCUBA-2 error circles very likely arise from physically associated galaxies rather than line-of-sight projections by chance. The photometric redshift of our radio$+$machine-learning identified SMGs ranges from z=0.2 to 5.7 and peaks at $z=2.3\pm0.1$. The AGN fraction of our sample is ($19\pm4$)%, which is consistent with that of ALMA SMGs in the literature. Comparing with radio/NIR-detected field galaxy population in the COSMOS field, our radio+machine-learning identified counterparts of SMGs have the highest star-formation rates and stellar masses. These characteristics suggest that our identified counterparts of S2COSMOS sources are a representative sample of SMGs at z<3. We employ our machine-learning technique to the whole COSMOS field and identified 6,877 potential SMGs, most of which are expected to have submillimeter emission fainter than the confusion limit of our S2COSMOS surveys (S$_{850}$<1.5mJy). We study the clustering properties of SMGs based on this statistically large sample, finding that they reside in high-mass dark matter halos ($(1.2\pm0.3)\times10^{13}\,h^{-1}\,\rm M_{\odot}$), which suggests that SMGs may be the progenitors of massive ellipticals we see in the local Universe.

We study black hole - host galaxy correlations, and the relation between the over-massiveness (the distance from the average $M_{BH}-\sigma$ relation) of super-massive black holes (SMBHs) and star formation histories of their host galaxies in the Illustris and TNG100 simulations. We find that both simulations are able to produce black hole scaling relations in general agreement with observations at $z=0$, but with noticeable discrepancies. Both simulations show an offset from the observations for the $M_{BH}-\sigma$ relation, and the relation between $M_{BH}$ and the Sersic index. The relation between $M_{BH}$ and stellar mass $M_*$ is tighter than the observations, especially for TNG100. For massive galaxies in both simulations, the hosts of over-massive SMBHs (those above the mean $M_{BH}-\sigma$ relation) tend to have larger Sersic indices and lower baryon conversion efficiency, suggesting a multidimensional link between SMBHs and properties of their hosts. In Illustris, the hosts of over-massive SMBHs have formed earlier and have lower present-day star formation rates, in qualitative agreement with the observations for massive galaxies with $\sigma>100 \rm km/s$. For low-mass galaxies, such a correlation still holds in Illustris but does not exist in the observed data. For TNG100, the correlation between SMBH over-massiveness and star formation history is much weaker. The hosts of over-massive SMBHs generally have consistently larger star formation rates throughout history. These galaxies have higher stellar mass as well, due to the strong $M_{BH}-M_*$ correlation. Our findings show that simulated SMBH scaling relations and correlations are sensitive to features in the modeling of SMBHs.

Most galaxies are hosted by massive, invisible dark matter haloes, yet little is known about the scatter in the stellar mass--halo mass relation for galaxies with host halo masses $M_{h}\le 10^{11}M_{\odot}$. Using mock catalogues based on dark matter simulations, we find that two observable signatures are sensitive to scatter in the stellar mass--halo mass relation even at these mass scales; i.e., conditional stellar mass functions and velocity distribution functions for neighbouring galaxies. We compute these observables for 179,373 galaxies in the Sloan Digital Sky Survey (SDSS) with stellar masses $M_{\ast} > 10^9 M_{\odot}$ and redshifts 0.01 $< z <$ 0.307. We then compare to mock observations generated from the $\textit{Bolshoi-Planck}$ dark matter simulation for stellar mass--halo mass scatters ranging from 0 to 0.6 dex. The observed results are consistent with simulated results for most values of scatter ($<$0.6 dex), and SDSS statistics are insufficient to provide firm constraints. However, this method could provide much tighter constraints on stellar mass--halo mass scatter in the future if applied to larger data sets, especially the anticipated Dark Energy Spectroscopic Instrument Bright Galaxy Survey. Constraining the value of scatter could have important implications for galaxy formation and evolution.

We use observations from the GEMINI-N/GMOS integral-field spectrograph (IFS) to obtain spatially resolved stellar kinematics of the central $\sim 1$ kpc of 20 early-type galaxies (ETGs) with stellar masses greater than $10^{11.7} M_\odot$ in the MASSIVE survey. Together with observations from the wide-field Mitchell IFS at McDonald Observatory in our earlier work, we obtain unprecedentedly detailed kinematic maps of local massive ETGs, covering a scale of $\sim 0.1-30$ kpc. The high ($\sim 120$) signal-to-noise of the GMOS spectra enable us to obtain two-dimensional maps of the line-of-sight velocity, velocity dispersion $\sigma$, as well as the skewness $h_3$ and kurtosis $h_4$ of the stellar velocity distributions. All but one galaxy in the sample have $\sigma(R)$ profiles that increase towards the center, whereas the slope of $\sigma(R)$ at one effective radius ($R_e$) can be of either sign. The $h_4$ is generally positive, with 14 of the 20 galaxies having positive $h_4$ within the GMOS aperture and 18 having positive $h_4$ within $1 R_e$. The positive $h_4$ and rising $\sigma(R)$ towards small radii are indicative of a central black hole and velocity anisotropy. We demonstrate the constraining power of the data on the mass distributions in ETGs by applying Jeans anisotropic modeling (JAM) to NGC~1453, the most regular fast rotator in the sample. Despite the limitations of JAM, we obtain a clear $\chi^2$ minimum in black hole mass, stellar mass-to-light ratio, velocity anisotropy parameters, and the circular velocity of the dark matter halo.

Supermassive black holes are located at the center of most, if not all, massive galaxies. They follow close correlations with global properties of their host galaxies (scaling relations), and are thought to play a crucial role in galaxy evolution. Yet, we lack a complete understanding of fundamental aspects of their growth across cosmic time. In particular, we still do not understand: (1) whether black holes or their host galaxies grow faster and (2) what is the maximum mass that black holes can reach. The high angular resolution capability and sensitivity of 30-m class telescopes will revolutionize our understanding of the extreme end of the black hole and galaxy mass scale. With such facilities, we will be able to dynamically measure masses of the largest black holes and characterize galaxy properties out to redshift $z \sim 1.5$. Together with the evolution of black hole-galaxy scaling relations since $z \sim 1.5$, the maximum mass black hole will shed light on the main channels of black hole growth.

In this paper we study the molecular gas content of a representative sample of 67 of the most massive early-type galaxies in the local universe, drawn uniformly from the MASSIVE survey. We present new IRAM-30m telescope observations of 30 of these galaxies, allowing us to probe the molecular gas content of the entire sample to a fixed molecular-to-stellar mass fraction of 0.1%. The total detection rate in this representative sample is 25$^{+5.9}_{-4.4}$%, and by combining the MASSIVE and ATLAS$^{\rm 3D}$ molecular gas surveys we find a joint detection rate of 22.4$^{+2.4}_{-2.1}$%. This detection rate seems to be independent of galaxy mass, size, position on the fundamental plane, and local environment. We show here for the first time that true slow rotators can host molecular gas reservoirs, but the rate at which they do so is significantly lower than for fast-rotators. Objects with a higher velocity dispersion at fixed mass (a higher kinematic bulge fraction) are less likely to have detectable molecular gas, and where gas does exist, have lower molecular gas fractions. In addition, satellite galaxies in dense environments have $\approx$0.6 dex lower molecular gas-to-stellar mass ratios than isolated objects. In order to interpret these results we created a toy model, which we use to constrain the origin of the gas in these systems. We are able to derive an independent estimate of the gas-rich merger rate in the low-redshift universe. These gas rich mergers appear to dominate the supply of gas to ETGs, but stellar mass loss, hot halo cooling and transformation of spiral galaxies also play a secondary role.

Pulsar timing arrays (PTAs) are on the verge of detecting low-frequency gravitational waves (GWs) from supermassive black hole binaries (SMBHBs). With continued observations of a large sample of millisecond pulsars, PTAs will reach this major milestone within the next decade. Already, SMBHB candidates are being identified by electromagnetic surveys in ever-increasing numbers; upcoming surveys will enhance our ability to detect and verify candidates, and will be instrumental in identifying the host galaxies of GW sources. Multi-messenger (GW and electromagnetic) observations of SMBHBs will revolutionize our understanding of the co-evolution of SMBHs with their host galaxies, the dynamical interactions between binaries and their galactic environments, and the fundamental physics of accretion. Multi-messenger observations can also make SMBHBs 'standard sirens' for cosmological distance measurements out to $z\simeq0.5$. LIGO has already ushered in breakthrough insights in our knowledge of black holes. The multi-messenger detection of SMBHBs with PTAs will be a breakthrough in the years $2020-2030$ and beyond, and prepare us for LISA to help complete our views of black hole demographics and evolution at higher redshifts.

We measure the stellar populations as a function of radius for 90 early-type galaxies (ETGs) in the MASSIVE survey, a volume-limited integral-field spectroscopic (IFS) galaxy survey targeting all northern-sky ETGs with absolute K-band magnitude M_K < -25.3 mag, or stellar mass M* 4x10^11 M_sun, within 108 Mpc. We are able to measure reliable stellar population parameters for individual galaxies out to 10-20 kpc (1-3 R_e) depending on the galaxy. Focusing on ~R_e (~10 kpc), we find significant correlations between the abundance ratios, sigma, and M* at large radius, but we also find that the abundance ratios saturate in the highest-mass bin. We see a strong correlation between the kurtosis of the line of sight velocity distribution (h4) and the stellar population parameters beyond R_e. Galaxies with higher radial anisotropy appear to be older, with metal-poorer stars and enhanced [alpha/Fe]. We suggest that the higher radial anisotropy may derive from more accretion of small satellites. Finally, we see some evidence for correlations between environmental metrics (measured locally and on >5 Mpc scales) and the stellar populations, as expected if satellites are quenched earlier in denser environments.

We simulate recoiling black hole trajectories from $z=20$ to $z=0$ in dark matter halos, quantifying how parameter choices affect escape velocities. These choices include the strength of dynamical friction, the presence of stars and gas, the accelerating expansion of the universe (Hubble acceleration), host halo accretion and motion, and seed black hole mass. $\Lambda$CDM halo accretion increases escape velocities by up to 0.6 dex and significantly shortens return timescales compared to non-accreting cases. Other parameters change orbit damping rates but have subdominant effects on escape velocities; dynamical friction is weak at halo escape velocities, even for extreme parameter values. We present formulae for black hole escape velocities as a function of host halo mass and redshift. Finally, we discuss how these findings affect black hole mass assembly as well as minimum stellar and halo masses necessary to retain supermassive black holes.

We test the distance--duality relation $\eta \equiv d_L / [ (1 + z)^2 d_A ] = 1$ between cosmological luminosity distance ($d_L$) from the JLA SNe Ia compilation (arXiv:1401.4064) and angular-diameter distance ($d_A$) based on Baryon Oscillation Spectroscopic Survey (BOSS; arXiv:1607.03155) and WiggleZ baryon acoustic oscillation measurements (arXiv:1105.2862, arXiv:1204.3674). The $d_L$ measurements are matched to $d_A$ redshift by a statistically consistent compression procedure. With Monte Carlo methods, nontrivial and correlated distributions of $\eta$ can be explored in a straightforward manner without resorting to a particular evolution template $\eta(z)$. Assuming independent constraints on cosmological parameters that are necessary to obtain $d_L$ and $d_A$ values, we find 9% constraints consistent with $\eta = 1$ from the analysis of SNIa + BOSS and an 18% bound results from SNIa + WiggleZ. These results are contrary to previous claims that $\eta < 1$ has been found close to or above the $1 \sigma$ level. We discuss the effect of different cosmological parameter inputs and the use of the apparent deviation from distance--duality as a proxy of systematic effects on cosmic distance measurements. The results suggest possible systematic overestimation of SNIa luminosity distances compared with $d_A$ data when a Planck {\Lambda}CDM cosmological parameter inference (arXiv:1502.01589) is used to enhance the precision. If interpreted as an extinction correction due to a gray dust component, the effect is broadly consistent with independent observational constraints.

We use spatially resolved two-dimensional stellar velocity maps over a $107"\times 107"$ field of view to investigate the kinematic features of 90 early-type galaxies above stellar mass $10^{11.5}M_\odot$ in the MASSIVE survey. We measure the misalignment angle $\Psi$ between the kinematic and photometric axes and identify local features such as velocity twists and kinematically distinct components. We find 46% of the sample to be well aligned ($\Psi < 15^{\circ}$), 33% misaligned, and 21% without detectable rotation (non-rotators). Only 24% of the sample are fast rotators, the majority of which (91%) are aligned, whereas 57% of the slow rotators are misaligned with a nearly flat distribution of $\Psi$ from $15^{\circ}$ to $90^{\circ}$. 11 galaxies have $\Psi \gtrsim 60^{\circ}$ and thus exhibit minor-axis ("prolate") rotation in which the rotation is preferentially around the photometric major axis. Kinematic misalignments occur more frequently for lower galaxy spin or denser galaxy environments. Using the observed misalignment and ellipticity distributions, we infer the intrinsic shape distribution of our sample and find that MASSIVE slow rotators are consistent with being mildly triaxial, with mean axis ratios of $b/a=0.88$ and $c/a=0.65$. In terms of local kinematic features, 51% of the sample exhibit kinematic twists of larger than $20^{\circ}$, and 2 galaxies have kinematically distinct components. The frequency of misalignment and the broad distribution of $\Psi$ reported here suggest that the most massive early-type galaxies are mildly triaxial, and that formation processes resulting in kinematically misaligned slow rotators such as gas-poor mergers occur frequently in this mass range.

Nuclear rings are excellent laboratories for probing diverse phenomena such as the formation and evolution of young massive star clusters (YMCs), nuclear starbursts, as well as the secular evolution and dynamics of their host galaxies. We have compiled a sample of 17 galaxies with nuclear rings, which are well resolved by high-resolution {\sl Hubble} and {\sl Spitzer Space Telescope} imaging. For each nuclear ring, we identified the ring star cluster population, along with their physical properties (ages, masses, extinction values). We also determined the integrated ring properties, including the average age, total stellar mass, and current star-formation rate (SFR). We find that Sb-type galaxies tend to have the highest ring stellar mass fraction with respect to the host galaxy, and this parameter is correlated with the ring's SFR surface density. The ring SFRs are correlated with their stellar masses, which is reminiscent of the main sequence of star-forming galaxies. There are striking correlations between star-forming properties (i.e., SFR and SFR surface density) and non-axisymmetric bar parameters, appearing to confirm previous inferences that strongly barred galaxies tend to have lower ring SFRs, although the ring star-formation histories turn out to be significantly more complicated. Nuclear rings with higher stellar masses tend to be associated with lower cluster mass fractions, but there is no such relation with the ages of the rings. The two youngest nuclear rings in our sample, NGC 1512 and NGC 4314, which have the most extreme physical properties, represent the young extremity of the nuclear ring age distribution.

We present near-infrared observations of 35 of the most massive early-type galaxies in the local universe. The observations were made using the infrared channel of the Hubble Space Telescope Wide Field Camera 3 in the F110W (1.1 $\mu$m) filter. We measured surface brightness profiles and elliptical isophotal fit parameters from the nuclear regions out to a radius of ~10 kpc in most cases. We find that 37% (13) of the galaxies in our sample have isophotal position angle rotations greater than 20 degrees over the radial range imaged by WFC3/IR, which is often due to the presence of neighbors or multiple nuclei. Most galaxies in our sample are significantly rounder near the center than in the outer regions. This sample contains six fast rotators and 28 slow rotators. We find that all fast rotators are either disky or show no measurable deviation from purely elliptical isophotes. Among slow rotators, significantly disky and boxy galaxies occur with nearly equal frequency. The galaxies in our sample often exhibit changing isophotal shapes, sometimes showing both significantly disky and boxy isophotes at different radii. The fact that parameters vary widely between galaxies and within individual galaxies is evidence that these massive galaxies have complicated formation histories, and some of them have experienced recent mergers and have not fully relaxed. These data demonstrate the value of high spatial resolution IR imaging of galaxies and provide measurements necessary for determining stellar masses, dynamics, and black hole masses in high mass galaxies.

The scaling relations between the mass of supermassive black holes ($M_{\bullet}$) and host galaxy properties (stellar mass, $M_{\star}$, and velocity dispersion, $\sigma$), provide a link between the growth of black holes (BHs) and that of their hosts. Here we investigate if and how the BH-galaxy relations are established in the high-$z$ universe using \textsc{BlueTides}, a high-resolution large volume cosmological hydrodynamic simulation. We find the $M_{\bullet}-M_{\star}$ and $M_{\bullet}-\sigma$ relations at $z=8$: $\log_{10}(M_{\bullet}) = 8.25 + 1.10 \ \log_{10}(M_{\star}/10^{11}M_{\odot})$ and $\log_{10}(M_{\bullet}) = 8.35 + 5.31 \ \log_{10}(\sigma/200kms^{-1})$ at $z=8$, both fully consistent with the local measurements. The slope of the $M_{\bullet}-\sigma$ relation is slightly steeper for high star formation rate and $M_{\star}$ galaxies while it remains unchanged as a function of Eddington accretion rate onto the BH. The intrinsic scatter in $M_{\bullet}-\sigma$ relation in all cases ($\epsilon \sim 0.4$) is larger at these redshifts than inferred from observations and larger than in $M_{\bullet}-M_{\star}$ relation ($\epsilon \sim 0.14$). We find the gas-to-stellar ratio $f=M_{\rm gas}/M_{\star}$ in the host (which can be very high at these redshifts) to have the most significant impact setting the intrinsic scatter of $M_{\bullet}-\sigma$. The scatter is significantly reduced when galaxies with high gas fractions ($\epsilon = 0.28$ as $f<10$) are excluded (making the sample more comparable to low-$z$ galaxies); these systems have the largest star formation rates and black hole accretion rates, indicating that these fast-growing systems are still moving toward the relation at these high redshifts. Examining the evolution (from $z=10$ to 8) of high mass black holes in $M_{\bullet}-\sigma$ plane confirms this trend.

We re-examine the properties of the star cluster population in the circumnuclear starburst ring in the face-on spiral galaxy NGC 7742, whose young cluster mass function has been reported to exhibit significant deviations from the canonical power law. We base our reassessment on the clusters' luminosities (an observational quantity) rather than their masses (a derived quantity), and confirm conclusively that the galaxy's starburst-ring clusters---and particularly the youngest subsample, $\log(t \mbox{ yr}^{-1}) \le 7.2$---show evidence of a turnover in the cluster luminosity function well above the 90\% completeness limit adopted to ensure the reliability of our results. This confirmation emphasises the unique conundrum posed by this unusual cluster population.

We measure the radial profiles of the stellar velocity dispersions, $\sigma(R)$, for 90 early-type galaxies (ETGs) in the MASSIVE survey, a volume-limited integral-field spectroscopic (IFS) galaxy survey targeting all northern-sky ETGs with absolute $K$-band magnitude $M_K < -25.3$ mag, or stellar mass $M_* > 4 \times 10^{11} M_\odot$, within 108 Mpc. Our wide-field 107" $\times$ 107" IFS data cover radii as large as 40 kpc, for which we quantify separately the inner (2 kpc) and outer (20 kpc) logarithmic slopes $\gamma_{\rm inner}$ and $\gamma_{\rm outer}$ of $\sigma(R)$. While $\gamma_{\rm inner}$ is mostly negative, of the 56 galaxies with sufficient radial coverage to determine $\gamma_{\rm outer}$ we find 36% to have rising outer dispersion profiles, 30% to be flat within the uncertainties, and 34% to be falling. The fraction of galaxies with rising outer profiles increases with $M_*$ and in denser galaxy environment, with 10 of the 11 most massive galaxies in our sample having flat or rising dispersion profiles. The strongest environmental correlations are with local density and halo mass, but a weaker correlation with large-scale density also exists. The average $\gamma_{\rm outer}$ is similar for brightest group galaxies, satellites, and isolated galaxies in our sample. We find a clear positive correlation between the gradients of the outer dispersion profile and the gradients of the velocity kurtosis $h_4$. Altogether, our kinematic results suggest that the increasing fraction of rising dispersion profiles in the most massive ETGs are caused (at least in part) by variations in the total mass profiles rather than in the velocity anisotropy alone.

Supermassive black hole binaries (SMBHBs) in the 10 million to 10 billion $M_\odot$ range form in galaxy mergers, and live in galactic nuclei with large and poorly constrained concentrations of gas and stars. There are currently no observations of merging SMBHBs--- it is in fact possible that they stall at their final parsec of separation and never merge. While LIGO has detected high frequency GWs, SMBHBs emit GWs in the nanohertz to millihertz band. This is inaccessible to ground-based interferometers, but possible with Pulsar Timing Arrays (PTAs). Using data from local galaxies in the 2 Micron All-Sky Survey, together with galaxy merger rates from Illustris, we find that there are on average $91\pm7$ sources emitting GWs in the PTA band, and $7\pm2$ binaries which will never merge. Local unresolved SMBHBs can contribute to GW background anisotropy at a level of $\sim20\%$, and if the GW background can be successfully isolated, GWs from at least one local SMBHB can be detected in 10 years.

The relation between X-ray luminosity (L_X) and ambient gas temperature (T) among massive galactic systems is an important cornerstone of both observational cosmology and galaxy-evolution modeling. In the most massive galaxy clusters, the relation is determined primarily by cosmological structure formation. In less massive systems, it primarily reflects the feedback response to radiative cooling of circumgalactic gas. Here we present a simple but powerful model for the L_X-T relation as a function of physical aperture R within which those measurements are made. The model is based on the precipitation framework for AGN feedback and assumes that the circumgalactic medium is precipitation-regulated at small radii and limited by cosmological structure formation at large radii. We compare this model with many different data sets and show that it successfully reproduces the slope and upper envelope of the L_X-T-R relation over the temperature range from ~0.2 keV through >10 keV. Our findings strongly suggest that the feedback mechanisms responsible for regulating star formation in individual massive galaxies have much in common with the precipitation-triggered feedback that appears to regulate galaxy-cluster cores.

We analyse the environmental properties of 370 local early-type galaxies (ETGs) in the MASSIVE and ATLAS3D surveys, two complementary volume-limited integral-field spectroscopic (IFS) galaxy surveys spanning absolute $K$-band magnitude $-21.5 > M_K > -26.6$, or stellar mass $8 \times 10^{9} < M_* < 2 \times 10^{12} M_\odot$. We find these galaxies to reside in a diverse range of environments measured by four methods: group membership (whether a galaxy is a brightest group/cluster galaxy, satellite, or isolated), halo mass, large-scale mass density (measured over a few Mpc), and local mass density (measured within the $N$th neighbour). The spatially resolved IFS stellar kinematics provide robust measurements of the spin parameter $\lambda_e$ and enable us to examine the relationship among $\lambda_e$, $M_*$, and galaxy environment. We find a strong correlation between $\lambda_e$ and $M_*$, where the average $\lambda_e$ decreases from $\sim 0.4$ to below 0.1 with increasing mass, and the fraction of slow rotators $f_{\rm slow}$ increases from $\sim 10$% to 90%. We show for the first time that at fixed $M_*$, there are almost no trends between galaxy spin and environment; the apparent kinematic morphology-density relation for ETGs is therefore primarily driven by $M_*$ and is accounted for by the joint correlations between $M_*$ and spin, and between $M_*$ and environment. A possible exception is that the increased $f_{\rm slow}$ at high local density is slightly more than expected based only on these joint correlations. Our results suggest that the physical processes responsible for building up the present-day stellar masses of massive galaxies are also very efficient at reducing their spin, in any environment.

The integrated properties of nuclear rings are correlated with their host galaxy's secular evolution and its dynamics, as well as with the formation and evolution of the ring's star cluster population(s). Here we present a new method to accurately measure the spectral energy distribution and current star-formation rate (SFR) of the nuclear ring in the barred spiral galaxy NGC 1512 based on high-resolution {\sl Hubble} and {\sl Spitzer Space Telescope} images. Image degradation does not have a significant negative effect on the robustness of the results. To obtain the ring's SFR for the period spanning $\sim$3--10 Myr, we apply our method to the continuum-subtracted H$\alpha$ and 8 $\mu$m images. The resulting SFR surface density, $\Sigma_{\rm SFR}$=$0.09\, {M}_{\odot} \,{\rm yr}^{-1}$ ${\rm kpc}^{-2}$, which is much higher than the disk-averaged SFR densities in normal galaxies. We also estimate the ring's total stellar mass, log (${M}/{M}_{\odot}$) = 7.1 $\pm$ 0.11 for an average age of $\sim$40 Myr.

We present a multi-wavelength analysis of 52 sub-millimeter galaxies (SMGs), identified using ALMA 870$\mu$m continuum imaging in a pilot program to precisely locate bright SCUBA2-selected sub-mm sources in the UKIDSS Ultra Deep Survey (UDS) field. Using the available deep (especially near-infrared), panoramic imaging of the UDS field at optical-to-radio wavelengths we characterize key properties of the SMG population. The median photometric redshift of the bright ALMA/SCUBA-2 UDS (AS2UDS) SMGs that are detected in a sufficient number of wavebands to derive a robust photometric redshift is $z$=2.65$\pm$0.13. However, similar to previous studies, 27% of the SMGs are too faint at optical-to-near-infrared wavelengths to derive a reliable photometric redshift. Assuming that these SMGs lie at z$\gtrsim$3 raises the median redshift of the full sample to $z$=2.9$\pm$0.2. A subset of 23, unlensed, bright AS2UDS SMGs have sizes measured from resolved imaging of their rest-frame far-infrared emission. We show that the extent and luminosity of the far-infrared emission are consistent with the dust emission arising from regions that are optically thick, on average, at a wavelength of $\lambda_0$$\ge$75$\mu$m (1-$\sigma$ dispersion of 55-90$\mu$m). Using the dust masses derived from our optically-thick spectral energy distribution models we determine that these galaxies have a median hydrogen column density of $N_{H}$=9.8$_{-0.7}^{+1.4}$$\times$10$^{23}$cm$^{-2}$, or a corresponding median $V$-band obscuration of $A_\mathrm{v}$=540$^{+80}_{-40}$mag, averaged along the line of sight to the source of their restframe $\sim$200$\mu$m emission. We discuss the implications of this extreme attenuation by dust for the multiwavelength study of dusty starbursts and reddening-sensitive tracers of star formation.

Mergers and the spin of the dark matter halo are factors traditionally believed to determine the morphology of galaxies within a $\Lambda$CDM cosmology. We study this hypothesis by considering approximately 18,000 central galaxies at $z=0$ with stellar masses $M_{\ast} = 10^{9}-10^{12} \, {\rm M}_{\odot}$ selected from the Illustris cosmological hydrodynamic simulation. The fraction of accreted stars -- which measures the importance of massive, recent and dry mergers -- increases steeply with galaxy stellar mass, from less than 5 per cent in dwarfs to 80 per cent in the most massive objects, and the impact of mergers on galaxy morphology increases accordingly. For galaxies with $M_{\ast} \gtrsim 10^{11} \, {\rm M}_{\odot}$, mergers have the expected effect: if gas-poor they promote the formation of spheroidal galaxies, whereas gas-rich mergers favour the formation and survivability of massive discs. This trend, however, breaks at lower masses. For objects with $M_{\ast} \lesssim 10^{11} \, {\rm M}_{\odot}$, mergers do not seem to play any significant role in determining the morphology, with accreted stellar fractions and mean merger gas fractions that are indistinguishable between spheroidal and disc-dominated galaxies. On the other hand, halo spin correlates with morphology primarily in the least massive objects in the sample ($M_{\ast} \lesssim 10^{10} \, {\rm M}_{\odot}$), but only weakly for galaxies above that mass. Our results support a scenario where (1) mergers play a dominant role in shaping the morphology of massive galaxies, (2) halo spin is important for the morphology of dwarfs, and (3) the morphology of medium-sized galaxies -- including the Milky Way -- shows little dependence on galaxy assembly history or halo spin, at least when these two factors are considered individually.

We present the first systematic investigation of the existence, spatial distribution, and kinematics of warm ionized gas as traced by the [O II] 3727AA emission line in 74 of the most massive galaxies in the local Universe. All of our galaxies have deep integral field spectroscopy from the volume- and magnitude-limited MASSIVE survey of early-type galaxies with stellar mass log(M_*/M_sun) > 11.5 (M_K < -25.3 mag) and distance D < 108 Mpc. Of the 74 galaxies in our sample, we detect warm ionized gas in 28, which yields a global detection fraction of 38\pm6% down to a typical [O II] equivalent width limit of 2AA. MASSIVE fast rotators are more likely to have gas than MASSIVE slow rotators with detection fractions of 80\pm10% and 28\pm6%, respectively. The spatial extents span a wide range of radii (0.6 - 18.2 kpc; 0.1 - 4R_e), and the gas morphologies are diverse, with 17/28 = 61\pm9% being centrally concentrated, 8/28 = 29\pm9% exhibiting clear rotation out to several kpc, and 3/28 = 11\pm6% being extended but patchy. Three out of four fast rotators show kinematic alignment between the stars and gas, whereas the two slow rotators with robust kinematic measurements available exhibit kinematic misalignment. Our inferred warm ionized gas masses are roughly ~10^5M_sun. The emission line ratios and radial equivalent width profiles are generally consistent with excitation of the gas by the old underlying stellar population. We explore different gas origin scenarios for MASSIVE galaxies and find that a variety of physical processes are likely at play, including internal gas recycling, cooling out of the hot gaseous halo, and gas acquired via mergers.

We use the Illustris simulation to study the relative contributions of in situ star formation and stellar accretion to the build-up of galaxies over an unprecedentedly wide range of masses ($M_{\ast} = 10^9-10^{12} \, {\rm M_{\odot}}$), galaxy types, environments, and assembly histories. We find that the `two-phase' picture of galaxy formation predicted by some models is a good approximation only for the most massive galaxies in our simulation -- namely, the stellar mass growth of galaxies below a few times $10^{11} \, {\rm M_{\odot}}$ is dominated by in situ star formation at all redshifts. The fraction of the total stellar mass of galaxies at $z=0$ contributed by accreted stars shows a strong dependence on galaxy stellar mass, ranging from about 10 per cent for Milky Way-sized galaxies to over 80 per cent for $M_{\ast} \approx 10^{12} \, {\rm M_{\odot}}$ objects, yet with a large galaxy-to-galaxy variation. At a fixed stellar mass, elliptical galaxies and those formed at the centres of younger haloes exhibit larger fractions of ex situ stars than disc-like galaxies and those formed in older haloes. On average, $\sim$50 per cent of the ex situ stellar mass comes from major mergers (stellar mass ratio $\mu > 1/4$), $\sim$20 per cent from minor mergers ($1/10 < \mu < 1/4$), $\sim$20 per cent from very minor mergers ($\mu < 1/10$), and $\sim$10 per cent from stars that were stripped from surviving galaxies (e.g. flybys or ongoing mergers). These components are spatially segregated, with in situ stars dominating the innermost regions of galaxies, and ex situ stars being deposited at larger galactocentric distances in order of decreasing merger mass ratio.

We analyse the cluster luminosity functions (CLFs) of the youngest star clusters in three galaxies exhibiting prominent circumnuclear starburst rings. We focus specifically on NGC 1512 and NGC 6951, for which we have access to H$\alpha$ data that allow us to unambiguously identify the youngest sample clusters. To place our results on a firm statistical footing, we first explore in detail a number of important technical issues affecting the process from converting the observational data into the spectral-energy distributions of the objects in our final catalogues. The CLFs of the young clusters in both galaxies exhibit approximate power-law behaviour down to the 90 per cent observational completeness limits, thus showing that star cluster formation in the violent environments of starburst rings appears to proceed similarly as that elsewhere in the local Universe. We discuss this result in the context of the density of the interstellar medium in our starburst-ring galaxies.

We present spatially-resolved two-dimensional stellar kinematics for the 41 most massive early-type galaxies (MK <~ -25.7 mag, stellar mass M* >~ 10^11.8 Msun) of the volume-limited (D < 108 Mpc) MASSIVE survey. For each galaxy, we obtain high-quality spectra in the wavelength range of 3650 to 5850 A from the 246-fiber Mitchell integral-field spectrograph (IFS) at McDonald Observatory, covering a 107" x 107" field of view (often reaching 2 to 3 effective radii). We measure the 2-D spatial distribution of each galaxy's angular momentum (lambda and fast or slow rotator status), velocity dispersion (sigma), and higher-order non-Gaussian velocity features (Gauss-Hermite moments h3 to h6). Our sample contains a high fraction (~80% ) of slow and non-rotators with lambda <~ 0.2. When combined with the lower-mass ETGs in the ATLAS3D survey, we find the fraction of slow-rotators to increase dramatically with galaxy mass, reaching ~50% at MK ~ -25.5 mag and ~90% at MK <~ -26 mag. All of our fast rotators show a clear anti-correlation between h3 and V/sigma, and the slope of the anti-correlation is steeper in more round galaxies. The radial profiles of sigma show a clear luminosity and environmental dependence: the 12 most luminous galaxies in our sample (MK <~ -26 mag) are all brightest cluster/group galaxies (except NGC 4874) and all have rising or nearly flat sigma profiles, whereas five of the seven "isolated" galaxies are all fainter than MK = -25.8 mag and have falling sigma. All of our galaxies have positive average h4; the most luminous galaxies have average h4 ~ 0.05 while less luminous galaxies have a range of values between 0 and 0.05. Most of our galaxies show positive radial gradients in h4, and those galaxies also tend to have rising sigma profiles. We discuss the implications for the relationship among dynamical mass, sigma, h4, and velocity anisotropy for these massive galaxies.

Bayesian graphical models are an efficient tool for modelling complex data and derive self-consistent expressions of the posterior distribution of model parameters. We apply Bayesian graphs to perform statistical analyses of Type Ia supernova (SN Ia) luminosity distance measurements from the joint light-curve analysis (JLA) data set. In contrast to the $\chi^2$ approach used in previous studies, the Bayesian inference allows us to fully account for the standard-candle parameter dependence of the data covariance matrix. Comparing with $\chi^2$ analysis results, we find a systematic offset of the marginal model parameter bounds. We demonstrate that the bias is statistically significant in the case of the SN Ia standardization parameters with a maximal 6 $\sigma$ shift of the SN light-curve colour correction. In addition, we find that the evidence for a host galaxy correction is now only 2.4 $\sigma$. Systematic offsets on the cosmological parameters remain small, but may increase by combining constraints from complementary cosmological probes. The bias of the $\chi^2$ analysis is due to neglecting the parameter-dependent log-determinant of the data covariance, which gives more statistical weight to larger values of the standardization parameters. We find a similar effect on compressed distance modulus data. To this end, we implement a fully consistent compression method of the JLA data set that uses a Gaussian approximation of the posterior distribution for fast generation of compressed data. Overall, the results of our analysis emphasize the need for a fully consistent Bayesian statistical approach in the analysis of future large SN Ia data sets.

Studies of the physical properties of local elliptical galaxies (e.g., gas temperatures, halo masses, stellar kinematics) are shedding new light on galaxy formation. Here we present the hot X-ray gas properties of 33 early-type systems within the MASSIVE galaxy survey sample that have archival Chandra X-ray observations. Through careful X-ray spectral modeling, we derive X-ray luminosities (L_X) and plasma temperatures (T_gas) for the diffuse gas components in these galaxies. We combine the MASSIVE sample with 41 galaxies from the ATLAS^3D survey to investigate the X-ray and optical properties of a statistically significant sample of nearby early-type galaxies across a wide-range of environments. We deduce that all early-type galaxies (independent of galaxy mass and rotational support) follow a universal scaling law such that L_X~T_gas^4.5. When X-ray measurements are performed consistently in apertures set by the galaxy stellar content, the wide-scale environment does not contribute to the intrinsic scatter (~0.5dex) within the scaling relation. We further demonstrate that the scatter in L_X around both K-band luminosity (L_K) and the galaxy stellar velocity dispersion is primarily driven by T_gas, with no clear trends with halo mass, radio power, or angular momentum of the stars. It is not trivial to tie the origin of the gas directly to either the stellar mass or the galaxy potential. Indeed, our data require a steeper relation between L_X, L_K, and sigma_e than predicted by standard mass-loss models. Finally, we find a statistically significant correlation between sigma_e and T_gas, suggesting that T_gas is set by the galaxy potential inside the optical effective radius. We conclude that within the inner-most 10-30kpc region, early-types maintain pressure-supported hot gas, with a minimum T_gas set by the virial temperature, but the majority show evidence for some additional heating.

Galaxy comoving number density is commonly used to forge progenitor/descendant links between observed galaxy populations at different epochs. However, this method breaks down in the presence of galaxy mergers, or when galaxies experience stochastic growth rates. We present a simple analytic framework to treat the physical processes that drive the evolution and diffusion of galaxies within comoving number density space. The evolution in mass rank order of a galaxy population with time is influenced by the galaxy coagulation rate and galaxy "mass rank scatter" rate. We quantify the relative contribution of these two effects to the mass rank order evolution. We show that galaxy coagulation is dominant at lower redshifts and stellar masses, while scattered growth rates dominate the mass rank evolution at higher redshifts and stellar masses. For a galaxy population at $10^{10} M_\odot$, coagulation has been the dominant effect since $z=2.2$, but a galaxy population at $10^{11} M_\odot$ was dominated by mass rank scatter until $z=0.6$. We show that although the forward and backward median number density evolution tracks are asymmetric, the backward median number density evolution can be obtained by convolving the descendant distribution function with progenitor relative abundances. We tabulate fits for the median number density evolution and scatter which can be applied to improve the way galaxy populations are linked in multi-epoch observational datasets.

Pulsar timing arrays (PTAs) are placing increasingly stringent constraints on the strain amplitude of continuous gravitational waves emitted by supermassive black hole binaries on subparsec scales. In this paper, we incorporate independent information about the dynamical masses $M_{bh}$ of supermassive black holes in specific galaxies at known distances and use this additional information to further constrain whether or not those galaxies could host a detectable supermassive black hole binary. We estimate the strain amplitudes from individual binaries as a function of binary mass ratio for two samples of nearby galaxies: (1) those with direct dynamical measurements of $M_{bh}$ in the literature, and (2) the 116 most massive early-type galaxies (and thus likely hosts of the most massive black holes) within 108 Mpc from the MASSIVE Survey. Our exploratory analysis shows that the current PTA upper limits on continuous waves (as a function of angular position in the sky) can already constrain the mass ratios of hypothetical black hole binaries in many galaxies in our samples. The constraints are stronger for galaxies with larger $M_{bh}$ and at smaller distances. For the black holes with $M_{bh} \gtrsim 5\times 10^9 M_\odot$ at the centers of NGC 1600, NGC 4889, NGC 4486 (M87) and NGC 4649 (M60), any binary companion in orbit within the PTA frequency bands would have to have a mass ratio of a few percent or less.

Quasars are associated with and powered by the accretion of material onto massive black holes; the detection of highly luminous quasars with redshifts greater than z = 6 suggests that black holes of up to ten billion solar masses already existed 13 billion years ago. Two possible present-day dormant descendants of this population of active black holes have been found in the galaxies NGC 3842 and NGC 4889 at the centres of the Leo and Coma galaxy clusters, which together form the central region of the Great Wall - the largest local structure of galaxies. The most luminous quasars, however, are not confined to such high-density regions of the early Universe; yet dormant black holes of this high mass have not yet been found outside of modern-day rich clusters. Here we report observations of the stellar velocity distribution in the galaxy NGC 1600 - a relatively isolated elliptical galaxy near the centre of a galaxy group at a distance of 64 Mpc from Earth. We use orbit superposition models to determine that the black hole at the centre of NGC 1600 has a mass of 17 billion solar masses. The spatial distribution of stars near the centre of NGC 1600 is rather diffuse. We find that the region of depleted stellar density in the cores of massive elliptical galaxies extends over the same radius as the gravitational sphere of influence of the central black holes, and interpret this as the dynamical imprint of the black holes.

We present multiwavelength identifications for the counterparts of 1088 submillimeter sources detected at 850$\mu$m in the SCUBA-2 Cosmology Legacy Survey study of the UKIDSS-UDS field. By utilising an ALMA pilot study on a subset of our bright SCUBA-2 sample as a training set, along with the deep optical-near-infrared data available in this field, we develop a novel technique, Optical-IR Triple Color (OIRTC), using $z-K$, $K-[3.6]$, $[3.6]-[4.5]$ colors to select the candidate submillimeter galaxy (SMG) counterparts. By combining radio identification and the OIRTC technique, we find counterpart candidates for 80% of the Class = 1 $\geq4\,\sigma$ SCUBA-2 sample, defined as those that are covered by both radio and OIR imaging and the base sample for our scientific analyses. Based on the ALMA training set, we expect the accuracy of these identifications to be $82\pm20$%, with a completeness of $69\pm16$%, essentially as accurate as the traditional $p$-value technique but with higher completeness. We find that the fraction of SCUBA-2 sources having candidate counterparts is lower for fainter 850$\mu$m sources, and we argue that for follow-up observations sensitive to SMGs with $S_{850}\gtrsim 1$ mJy across the whole ALMA beam, the fraction with multiple counterparts is likely to be $>40$% for SCUBA-2 sources at $S_{850} \gtrsim 4$ mJy. We find that the photometric redshift distribution for the SMGs is well fit by a lognormal distribution, with a median redshift of $z=2.3\pm0.1$. After accounting for the sources without any radio and/or OIRTC counterpart, we estimate the median redshift to be $z=2.6\pm0.1$ for SMGs with $S_{850} >1$ mJy. We also use this new large sample to study the clustering of SMGs and the the far-infrared properties of the unidentified submillimeter sources by stacking their Herschel SPIRE far-infrared emission.

Massive quiescent galaxies have much smaller physical sizes at high redshift than today. The strong evolution of galaxy size may be caused by progenitor bias, major and minor mergers, adiabatic expansion, and/or renewed star formation, but it is difficult to test these theories observationally. Herein, we select a sample of 35 massive, compact galaxies ($M_* = 1-3 \times 10^{11}$ M$_\odot$, $M_*/R^{1.5} > 10^{10.5}$ M$_\odot$/kpc$^{1.5}$) at $z=2$ in the cosmological hydrodynamical simulation Illustris and trace them forward to $z=0$ to uncover their evolution and identify their descendants. By $z=0$, the original factor of 3 difference in stellar mass spreads to a factor of 20. The dark matter halo masses similarly spread from a factor of 5 to 40. The galaxies' evolutionary paths are diverse: about half acquire an ex-situ envelope and are the core of a more massive descendant, a third survive undisturbed and gain very little mass, 15% are consumed in a merger with a more massive galaxy, and a small remainder are thoroughly mixed by major mergers. The galaxies grow in size as well as mass, and only $\sim$10% remain compact by $z=0$. The majority of the size growth is driven by the acquisition of ex-situ mass. The most massive galaxies at $z=0$ are the most likely to have compact progenitors, but this trend possesses significant dispersion which precludes a direct linkage to compact galaxies at $z=2$. The compact galaxies' merger rates are influenced by their $z=2$ environments, so that isolated or satellite compact galaxies (which are protected from mergers) are the most likely to survive to the present day.

The goal of this presentation is to report the latest progress in creation of the next generation of VLBI-based International Celestial Reference Frame, ICRF3. Two main directions of ICRF3 development are improvement of the S/X-band frame and extension of the ICRF to higher frequencies. Another important task of this work is the preparation for comparison of ICRF3 with the new generation optical frame GCRF expected by the end of the decade as a result of the Gaia mission.

In this work we present CO(1-0) and CO(2-1) observations of a pilot sample of 15 early-type galaxies (ETGs) drawn from the MASSIVE galaxy survey, a volume-limited integral-field spectroscopic study of the most massive ETGs ($M_* >10^{11.5}M_\odot$) within 108 Mpc. These objects were selected because they showed signs of an interstellar medium and/or star formation. A large amount of gas ($>$2$\times$10$^8$ M$_{\odot}$) is present in 10 out of 15 objects, and these galaxies have gas fractions higher than expected based on extrapolation from lower mass samples. We tentatively interpret this as evidence that stellar mass loss and hot halo cooling may be starting to play a role in fuelling the most massive galaxies. These MASSIVE ETGs seem to have lower star-formation efficiencies (SFE=SFR/M$_{\rm H2}$) than spiral galaxies, but the SFEs derived are consistent with being drawn from the same distribution found in other lower mass ETG samples. This suggests that the SFE is not simply a function of stellar mass, but that local, internal processes are more important for regulating star formation. Finally we used the CO line profiles to investigate the high-mass end of the Tully-Fisher relation (TFR). We find that there is a break in the slope of the TFR for ETGs at high masses (consistent with previous studies). The strength of this break correlates with the stellar velocity dispersion of the host galaxies, suggesting it is caused by additional baryonic mass being present in the centre of massive ETGs. We speculate on the root cause of this change and its implications for galaxy formation theories.

The cumulative comoving number-density of galaxies as a function of stellar mass or central velocity dispersion is commonly used to link galaxy populations across different epochs. By assuming that galaxies preserve their number-density in time, one can infer the evolution of their properties, such as masses, sizes, and morphologies. However, this assumption does not hold in the presence of galaxy mergers or when rank ordering is broken owing to variable stellar growth rates. We present an analysis of the evolving comoving number density of galaxy populations found in the Illustris cosmological hydrodynamical simulation focused on the redshift range $0\leq z \leq 3$. Our primary results are as follows: 1) The inferred average stellar mass evolution obtained via a constant comoving number density assumption is systematically biased compared to the merger tree results at the factor of $\sim$2(4) level when tracking galaxies from redshift $z=0$ out to redshift $z=2(3)$; 2) The median number density evolution for galaxy populations tracked forward in time is shallower than for galaxy populations tracked backward in time; 3) A similar evolution in the median number density of tracked galaxy populations is found regardless of whether number density is assigned via stellar mass, stellar velocity dispersion, or dark matter halo mass; 4) Explicit tracking reveals a large diversity in galaxies' assembly histories that cannot be captured by constant number-density analyses; 5) The significant scatter in galaxy linking methods is only marginally reduced by considering a number of additional physical and observable galaxy properties as realized in our simulation. We provide fits for the forward and backward median evolution in stellar mass and number density and discuss implications of our analysis for interpreting multi-epoch galaxy property observations.

We examine stellar population gradients in ~100 massive early type galaxies spanning 180 < sigma* < 370 km/s and M_K of -22.5 to -26.5 mag, observed as part of the MASSIVE survey (Ma et al. 2014). Using integral-field spectroscopy from the Mitchell Spectrograph on the 2.7m telescope at McDonald Observatory, we create stacked spectra as a function of radius for galaxies binned by their stellar velocity dispersion, stellar mass, and group richness. With excellent sampling at the highest stellar mass, we examine radial trends in stellar population properties extending to beyond twice the effective radius (~2.5 R_e). Specifically, we examine trends in age, metallicity, and abundance ratios of Mg, C, N, and Ca, and discuss the implications for star formation histories and elemental yields. At a fixed physical radius of 3-6 kpc (the likely size of the galaxy cores formed at high redshift) stellar age and [alpha/Fe] increase with increasing sigma* and depend only weakly on stellar mass, as we might expect if denser galaxies form their central cores earlier and faster. If we instead focus on 1-1.5 R_e, the trends in abundance and abundance ratio are washed out, as might be expected if the stars at large radius were accreted by smaller galaxies. Finally, we show that when controlling for \sigmastar, there are only very subtle differences in stellar population properties or gradients as a function of group richness; even at large radius internal properties matter more than environment in determining star formation history.

We report the results of a deep SCUBA-2 850- and 450-$\mu$m survey for dust-obscured ultra-luminous infrared galaxies (U/LIRGs) in the field of the z=1.46 cluster XCS J2215.9-1738. We detect a striking overdensity of sub-millimeter sources coincident with the core of this cluster: $\sim 3-4 \times$ higher than expected in a blank field. We use the likely radio and mid-infrared counterparts to show that the bulk of these sub-millimeter sources have spectroscopic or photometric redshifts which place them in the cluster and that their multi-wavelength properties are consistent with this association. The average far-infrared luminosities of these galaxies are $(1.0\pm0.1) \times 10^{12} L_{\odot}$, placing them on the U/LIRG boundary. Using the total star formation occurring in the obscured U/LIRG population within the cluster we show that the resulting mass-normalized star-formation rate for this system supports previous claims of a rapid increase in star-formation activity in cluster cores out to $z\sim1.5$, which must be associated with the on-going formation of the early-type galaxies which reside in massive clusters today.

We have constructed merger trees for galaxies in the Illustris Simulation by directly tracking the baryonic content of subhalos. These merger trees are used to calculate the galaxy-galaxy merger rate as a function of descendant stellar mass, progenitor stellar mass ratio, and redshift. We demonstrate that the most appropriate definition for the mass ratio of a galaxy-galaxy merger consists in taking both progenitor masses at the time when the secondary progenitor reaches its maximum stellar mass. Additionally, we avoid effects from `orphaned' galaxies by allowing some objects to `skip' a snapshot when finding a descendant, and by only considering mergers which show a well-defined `infall' moment. Adopting these definitions, we obtain well-converged predictions for the galaxy-galaxy merger rate with the following main features, which are qualitatively similar to the halo-halo merger rate except for the last one: a strong correlation with redshift that evolves as $\sim (1+z)^{2.4-2.8}$, a power law with respect to mass ratio, and an increasing dependence on descendant stellar mass, which steepens significantly for descendant stellar masses greater than $\sim 2 \times 10^{11} \, {\rm M_{\odot}}$. These trends are consistent with observational constraints for medium-sized galaxies ($M_{\ast} \gtrsim 10^{10} \, {\rm M_{\odot}}$), but in tension with some recent observations of the close pair fraction for massive galaxies ($M_{\ast} \gtrsim 10^{11} \, {\rm M_{\odot}}$), which report a nearly constant or decreasing evolution with redshift. Finally, we provide a fitting function for the galaxy-galaxy merger rate which is accurate over a wide range of stellar masses, progenitor mass ratios, and redshifts.

Massive, quiescent galaxies at high redshift have been found to be considerably more compact than galaxies of similar mass in the local universe. How these compact galaxies formed has yet to be determined, though several progenitor populations have been proposed. Here we investigate the formation processes and quantify the assembly histories of such galaxies in Illustris, a suite of hydrodynamical cosmological simulations encompassing a sufficiently large volume to include rare objects, while simultaneously resolving the internal structure of galaxies. We select massive (~10^11 solar masses) and compact (stellar half-mass radius < 2 kpc) galaxies from the simulation at z=2. Within the Illustris suite, we find that these quantities are not perfectly converged, but are reasonably reliable for our purposes. The resulting population is composed primarily of quiescent galaxies, but we also find several star-forming compact galaxies. The simulated compact galaxies are similar to observed galaxies in star formation activity and appearance. We follow their evolution at high redshift in the simulation and find that there are multiple pathways to form these compact galaxies, dominated by two mechanisms: (i) intense, centrally concentrated starbursts generally triggered by gas-rich major mergers between z~2-4, reducing the galaxies' half-mass radii by a factor of a few to below 2 kpc, and (ii) assembly at very early times when the universe was much denser; the galaxies formed compact and remained so until z~2.

Using the longitudinal expression of Hubble expansion rate for the general Lema\^itre-Tolman-Bondi (LTB) metric as a function of cosmic time, we examine the scale on which the Copernican Principle holds in the context of a void model. By way of performing parameter estimation on the CGBH void model, we show that the Hubble parameter data favors a void with characteristic radius of 2 ~ 3 Gpc. This brings the void model closer, but not yet enough, to harmony with observational indications given by the background kinetic Sunyaev-Zel'dovich effect and the normalization of near-infrared galaxy luminosity function. However, the test of such void models may ultimately lie in the future detection of the discrepancy between longitudinal and transverse expansion rates, a touchstone of inhomogeneous models. With the proliferation of observational Hubble parameter data and future large-scale structure observation, a definitive test could be performed on the question of cosmic homogeneity. Particularly, the spherical LTB void models have been ruled out, but more general non-spherical inhomogeneities still need to be tested by observation. In this paper, we utilise a spherical void model to provide guidelines into how observational tests may be done with more general models in the future.

During galaxy mergers the gas falls to the center, triggers star formation, and feeds the rapid growth of supermassive black holes (SMBHs). SMBHs respond to this fueling by supplying energy back to the ambient gas. Numerical studies suggest that this feedback is necessary to explain why the properties of SMBHs and the formation of bulges are closely related. This intimate link between the SMBH's mass and the large scale dynamics and luminosity of the host has proven to be a difficult issue to tackle with simulations due to the inability to resolve all the relevant length scales simultaneously. In this paper we simulate SMBH growth at high-resolution with {\it FLASH}, accounting for the gravitational focusing effects of nuclear star clusters (NSCs), which appear to be ubiquitous in galactic nuclei. In the simulations, the NSC core is resolved by a minimum cell size of about 0.001 pc or approximately $10^{-3}$ of the cluster's radius. We discuss the conditions required for effective gas funneling to occur, which are mainly dominated by a relationship between NSC velocity dispersion and the local sound speed, and provide a sub-grid prescription for the augmentation of central SMBH accretion rates in the presence of NSCs. For the conditions expected to persist in the centers of merging galaxies, the resultant large central gas densities in NSCs should produce drastically enhanced embedded SMBH accretion rates - up to an order of magnitude increase can be achieved for gas properties resembling those in large-scale galaxy merger simulations. This will naturally result in faster black hole growth rates and higher luminosities than predicted by the commonly used Bondi-Hoyle-Lyttleton accretion formalism.

We analyse HST WFC3/$H_{160}$-band observations of a sample of 48 ALMA-detected submillimeter galaxies (SMGs) in the Extended Chandra Deep Field South field, to study their stellar morphologies and sizes. We detect 79$\pm$17% of the SMGs in the $H_{160}$-band imaging with a median sensitivity of 27.8 mag, and most (80%) of the non-detections are SMGs with 870$\mu$m fluxes of $S_{870} < $3 mJy. With a surface brightness limit of $\mu_H \sim $26 mag arcsec$^{-2}$, we find that 82$\pm$9% of the $H_{160}$-band detected SMGs at $z =$ 1-3 appear to have disturbed morphologies, meaning they are visually classified as either irregulars or interacting systems, or both. By determining a S\'ersic fit to the $H_{160}$ surface-brightness profiles we derive a median S\'ersic index of $n = $1.2$\pm$0.3 and a median half-light radius of $r_e = $4.4$^{+1.1}_{-0.5}$ kpc for our SMGs at $z = $1-3. We also find significant displacements between the positions of the $H_{160}$-component and 870$\mu$m emission in these systems, suggesting that the dusty star-burst regions and less-obscured stellar distribution are not co-located. We find significant differences in the sizes and the S\'ersic index between our $z = $2-3 SMGs and $z \sim $2 quiescent galaxies, suggesting a major transformation of the stellar light profile is needed in the quenching processes if SMGs are progenitors of the red-and-dead $z\sim$2 galaxies. Given the short-lived nature of SMGs, we postulate that the majority of the $z = $2-3 SMGs with $S_{870} \gtrsim $2 mJy are early/mid-stage major mergers.

We present high-resolution (0.3'') ALMA 870um imaging of 52 sub-millimeter galaxies (SMGs) in the Ultra Deep Survey (UDS) field and investigate the size and morphology of the sub-millimeter (sub-mm) emission on 2-10kpc scales. We derive a median intrinsic angular size of FWHM=0.30$\pm$0.04'' for the 23 SMGs in the sample detected at a signal-to-noise ratio (SNR) >10. Using the photometric redshifts of the SMGs we show that this corresponds to a median physical half-light diameter of 2.4$\pm$0.2kpc. A stacking analysis of the SMGs detected at an SNR <10 shows they have sizes consistent with the 870um-bright SMGs in the sample. We compare our results to the sizes of SMGs derived from other multi-wavelength studies, and show that the rest-frame ~250um sizes of SMGs are consistent with studies of resolved 12CO (J=3-2 to 7-6) emission lines, but that sizes derived from 1.4GHz imaging appear to be approximately two times larger on average, which we attribute to cosmic ray diffusion. The rest-frame optical sizes of SMGs are around four times larger than the sub-millimeter sizes, indicating that the star formation in these galaxies is compact relative to the pre-existing stellar distribution. The size of the starburst region in SMGs is consistent with the majority of the star formation occurring in a central region, a few kpc in extent, with a median star formation rate surface density of 90$\pm$30Msol/yr/kpc$^2$, which may suggest that we are witnessing an intense period of bulge growth in these galaxies.

Massive early-type galaxies represent the modern-day remnants of the earliest major star formation episodes in the history of the universe. These galaxies are central to our understanding of the evolution of cosmic structure, stellar populations, and supermassive black holes, but the details of their complex formation histories remain uncertain. To address this situation, we have initiated the MASSIVE Survey, a volume-limited, multi-wavelength, integral-field spectroscopic (IFS) and photometric survey of the structure and dynamics of the ~100 most massive early-type galaxies within a distance of 108 Mpc. This survey probes a stellar mass range M* > 10^{11.5} Msun and diverse galaxy environments that have not been systematically studied to date. Our wide-field IFS data cover about two effective radii of individual galaxies, and for a subset of them, we are acquiring additional IFS observations on sub-arcsecond scales with adaptive optics. We are also acquiring deep K-band imaging to trace the extended halos of the galaxies and measure accurate total magnitudes. Dynamical orbit modeling of the combined data will allow us to simultaneously determine the stellar, black hole, and dark matter halo masses. The primary goals of the project are to constrain the black hole scaling relations at high masses, investigate systematically the stellar initial mass function and dark matter distribution in massive galaxies, and probe the late-time assembly of ellipticals through stellar population and kinematical gradients. In this paper, we describe the MASSIVE sample selection, discuss the distinct demographics and structural and environmental properties of the selected galaxies, and provide an overview of our basic observational program, science goals and early survey results.

All gravitationally bound clusters expand, due to both gas loss from their most massive members and binary heating. All are eventually disrupted tidally, either by passing molecular clouds or the gravitational potential of their host galaxies. However, their interior evolution can follow two very different paths. Only clusters of sufficiently large initial population and size undergo the combined interior contraction and exterior expansion that leads eventually to core collapse. In all other systems, core collapse is frustrated by binary heating. These clusters globally expand for their entire lives, up to the point of tidal disruption. Using a suite of direct N-body calculations, we trace the "collapse line" in r_v-N space that separates these two paths. Here, r_v and N are the cluster's initial virial radius and population, respectively. For realistic starting radii, the dividing N-value is from 10^4 to over 10^5. We also show that there exists a minimum population, N_min, for core collapse. Clusters with N < N_min tidally disrupt before core collapse occurs. At the Sun's Galactocentric radius, R_G = 8.5 kpc, we find N_min >~ 300. The minimum population scales with Galactocentric radius as R_G^{-9/8}. The position of an observed cluster relative to the collapse line can be used to predict its future evolution. Using a small sample of open clusters, we find that most lie below the collapse line, and thus will never undergo core collapse. Most globular clusters, on the other hand, lie well above the line. In such a case, the cluster may or may not go through core collapse, depending on its initial size. We show how an accurate age determination can help settle this issue.

In this paper we outline several problems related to the realization of the international celestial and terrestrial reference frames ICRF and ITRF at the millimeter level of accuracy, with emphasis on ICRF issues. The main topics considered are: analysis of the current status of the ICRF, mutual impact of ICRF and ITRF, and some considerations for future ICRF realizations.

We present a catalog of 1921 spectroscopic redshifts measured in the fields of the massive galaxy clusters MACSJ0416.1--2403 ($z=0.397$), MACSJ0717.5+3745 ($z=0.546$), and MACSJ1149.5+2223 ($z=0.544$), i.e., three of the four clusters selected by STScI as the targets of the Frontier Fields (FF) initiative for studies of the distant Universe via gravitational lensing. Compiled in the course of the MACS project (Massive Cluster Survey) that detected the FF clusters, this catalog is provided to the community for three purposes: (1) to allow the identification of cluster members for studies of the galaxy population of these extreme systems, (2) to facilitate the removal of unlensed galaxies and thus reduce shear dilution in weak-lensing analyses, and (3) to improve the calibration of photometric redshifts based on both ground- and spacebased observations of the FF clusters.

In the standard excursion set model for the growth of structure, the statistical properties of halos are governed by the halo mass and are independent of the larger scale environment in which the halos reside. Numerical simulations, however, have found the spatial distributions of halos to depend not only on their mass but also on the details of their assembly history and environment. Here we present a theoretical framework for incorporating this "assembly bias" into the excursion set model. Our derivations are based on modifications of the path integral approach of Maggiore & Riotto (2010) that models halo formation as a non-Markovian random walk process. The perturbed density field is assumed to evolve stochastically with the smoothing scale and exhibits correlated walks in the presence of a density barrier. We write down conditional probabilities for multiple barrier crossings, and derive from them analytic expressions for descendant and progenitor halo mass functions and halo merger rates as a function of both halo mass and the linear overdensity of the larger-scale environment of the halo. Our results predict a higher halo merger rate and higher progenitor halo mass function in regions of higher overdensity, consistent with the behavior seen in N-body simulations.

We present a $250\,$ks Chandra observation of the cluster merger A2034 with the aim of understanding the nature of a sharp edge previously characterized as a cold front. The new data reveal that the edge is coherent over a larger opening angle and is significantly more bow-shock-shaped than previously thought. Within $\sim 27\,$degrees about the axis of symmetry of the edge the density, temperature and pressure drop abruptly by factors of $1.83^{+0.09}_{-0.08}$, $1.85^{+0.41}_{-0.41}$ and $3.4^{+0.8}_{-0.7}$, respectively. This is inconsistent with the pressure equilibrium expected of a cold front and we conclude that the edge is a shock front. We measure a Mach number $M = 1.59^{+0.06}_{-0.07}$ and corresponding shock velocity $v_{\rm shock}\simeq 2057\,$km/s. Using spectra collected at the MMT with the Hectospec multi-object spectrograph we identify 328 spectroscopically confirmed cluster members. Significantly, we find a local peak in the projected galaxy density associated with a bright cluster galaxy which is located just ahead of the nose of the shock. The data are consistent with a merger viewed within $\sim 23\,$degrees of the plane of the sky. The merging subclusters are now moving apart along a north-south axis approximately $0.3\,$Gyr after a small impact parameter core passage. The gas core of the secondary subcluster, which was driving the shock, appears to have been disrupted by the merger. Without a driving 'piston' we speculate that the shock is dying. Finally, we propose that the diffuse radio emission near the shock is due to the revival of pre-existing radio plasma which has been overrun by the shock.

We study the holographic and agegraphic dark energy models without interaction using the latest observational Hubble parameter data (OHD), the Union2.1 compilation of type Ia supernovae (SNIa), and the energy conditions. Scenarios of dark energy are distinguished by the cut-off of cosmic age, conformal time, and event horizon. The best-fit value of matter density for the three scenarios almost steadily located at $\Omega_{m0}=0.26$ by the joint constraint. For the agegraphic models, they can be recovered to the standard cosmological model when the constant $c$ which presents the fraction of dark energy approaches to infinity. Absence of upper limit of $c$ by the joint constraint demonstrates the recovery possibility. Using the fitted result, we also reconstruct the current equation of state of dark energy at different scenarios, respectively. Employing the model criteria $\chi^2_{\textrm{min}}/dof$, we find that conformal time model is the worst, but they can not be distinguished clearly. Comparing with the observational constraints, we find that SEC is fulfilled at redshift $0.2 \lesssim z \lesssim 0.3$ with $1\sigma$ confidence level. We also find that NEC gives a meaningful constraint for the event horizon cut-off model, especially compared with OHD only. We note that the energy condition maybe could play an important role in the interacting models because of different degeneracy between $\Omega_m$ and constant $c$.

While the star formation rates and morphologies of galaxies have long been known to correlate with their local environment, the process by which these correlations are generated is not well understood. Galaxy groups are thought to play an important role in shaping the physical properties of galaxies before entering massive clusters at low redshift, and transformations of satellite galaxies likely dominate the buildup of local environmental correlations. To illuminate the physical processes that shape galaxy evolution in dense environments, we study a sample of 116 X-ray selected galaxy groups at z=0.2-1 with halo masses of 10^13-10^14 M_sun and centroids determined with weak lensing. We analyze morphologies based on HST imaging and colors determined from 31 photometric bands for a stellar mass-limited population of 923 satellite galaxies and a comparison sample of 16644 field galaxies. Controlling for variations in stellar mass across environments, we find significant trends in the colors and morphologies of satellite galaxies with group-centric distance and across cosmic time. Specifically at low stellar mass (log(M_stellar/M_sun) = 9.8-10.3), the fraction of disk-dominated star-forming galaxies declines from >50% among field galaxies to <20% among satellites near the centers of groups. This decline is accompanied by a rise in quenched galaxies with intermediate bulge+disk morphologies, and only a weak increase in red bulge-dominated systems. These results show that both color and morphology are influenced by a galaxy's location within a group halo. We suggest that strangulation and disk fading alone are insufficient to explain the observed morphological dependence on environment, and that galaxy mergers or close tidal encounters must play a role in building up the population of quenched galaxies with bulges seen in dense environments at low redshift.

We have tested the effect of spatial gradients in stellar mass-to-light ratio (Y) on measurements of black hole masses (MBH) derived from stellar orbit superposition models. Such models construct a static gravitational potential for a galaxy and its central black hole, but typically assume spatially uniform Y. We have modeled three giant elliptical galaxies with gradients alpha = d(log Y)/d(log r) from -0.2 to +0.1. Color and line strength gradients suggest mildly negative alpha in these galaxies. Introducing a negative (positive) gradient in Y increases (decreases) the enclosed stellar mass near the center of the galaxy and leads to systematically smaller (larger) MBH measurements. For models with alpha = -0.2, the best-fit values of MBH are 28%, 27%, and 17% lower than the constant-Y case, in NGC 3842, NGC 6086, and NGC 7768, respectively. For alpha = +0.1, MBH are 14%, 22%, and 17% higher than the constant-Y case for the three respective galaxies. For NGC 3842 and NGC 6086, this bias is comparable to the statistical errors from individual modeling trials. At larger radii, negative (positive) gradients in Y cause the total stellar mass to decrease (increase) and the dark matter fraction within one effective radius to increase (decrease).

New kinematic data and modeling efforts in the past few years have substantially expanded and revised dynamical measurements of black hole masses (Mbh) at the centers of nearby galaxies. Here we compile an updated sample of 72 black holes and their host galaxies, and present revised scaling relations between Mbh and stellar velocity dispersion (sigma), V-band luminosity (L), and bulge stellar mass (Mbulge), for different galaxy subsamples. Our best-fitting power law relations for the full galaxy sample are log(Mbh) = 8.32 + 5.64*log(sigma/200 kms), log(Mbh) = 9.23 + 1.11*log(L/10^{11} Lsun), and log(Mbh) = 8.46 + 1.05*log(Mbulge/10^{11} Msun). A log-quadratic fit to the Mbh-sigma relation with an additional term of beta_2*[log(sigma/200 kms)]^2 gives beta_2 = 1.68 +/- 1.82 and does not decrease the intrinsic scatter in Mbh. When the early- and late-type galaxies are fit separately, we obtain similar slopes of 5.20 and 5.06 for the Mbh-sigma relation but significantly different intercepts -- Mbh in early-type galaxies are about 2 times higher than in late types at a given sigma. Within early-type galaxies, our fits to Mbh(sigma) give Mbh that is about 2 times higher in galaxies with central core profiles than those with central power-law profiles. Our Mbh-L and Mbh-Mbulge relations for early-type galaxies are similar to those from earlier compilations. When the conventional quadrature method is used to determine the intrinsic scatter in Mbh, our dataset shows weak evidence for increased scatter at Mbulge < 10^{11} Msun or L_V < 10^{10.3} Lsun, while the scatter stays constant for 10^{11} < Mbulge < 10^{12.3} Msun and 10^{10.3} < L_V < 10^{11.5} Lsun. A Bayesian analysis indicates a larger sample of Mbh measurements would be needed to detect any statistically significant trend in the scatter with galaxy properties.

We estimate the average radio-AGN (mechanical) power deposited into the hot atmospheres of galaxy clusters over more than three quarters of the age of the Universe. Our sample was drawn from eight major X-ray cluster surveys, and includes 685 clusters in the redshift range 0.1 < z < 0.6 that overlap the area covered by the NVSS. The radio AGN mechanical power was estimated from the radio luminosity of central NVSS sources, using the relation of Cavagnolo et al. 2010 that is based on mechanical powers determined from the enthalpies of X-ray cavities. We find only a weak correlation between radio luminosity and cluster X-ray luminosity, although the most powerful radio sources resides in luminous clusters. The average AGN mechanical power of 3x10^{44} erg/s exceeds the X-ray luminosity of 44% of the clusters, indicating that the accumulation of radio-AGN energy is significant in these clusters. Integrating the AGN mechanical power to redshift z=2, using simple models for its evolution and disregarding the hierarchical growth of clusters, we find that the AGN energy accumulated per particle in low luminosity X-ray clusters exceeds 1 keV per particle. This result represents a conservative lower limit to the accumulated thermal energy. The estimate is comparable to the level of energy needed to "preheat" clusters, indicating that continual outbursts from radio-AGN are a significant source of gas energy in hot atmospheres. Assuming an average mass conversion efficiency of $\eta=0.1$, our result implies that the supermassive black holes that released this energy did so by accreting an average of ~10^9 M_{sun} over time, which is comparable to the level of growth expected during the quasar era.

Locating the centers of dark matter halos is critical for understanding the mass profiles of halos as well as the formation and evolution of the massive galaxies that they host. The task is observationally challenging because we cannot observe halos directly, and tracers such as bright galaxies or X-ray emission from hot plasma are imperfect. In this paper we quantify the consequences of miscentering on the weak lensing signal from a sample of 129 X-ray selected galaxy groups in the COSMOS field with redshifts 0<z<1 and halo masses in the range 10^13 - 10^14 M_sun. By measuring the stacked lensing signal around eight different candidate centers (such as the brightest member galaxy, the mean position of all member galaxies, or the X-ray centroid), we determine which candidates best trace the center of mass in halos. In this sample of groups, we find that massive galaxies near the X-ray centroids trace the center of mass to <~75 kpc, while the X-ray position and centroids based on the mean position of member galaxies have larger offsets primarily due to the statistical uncertainties in their positions (typically ~50-150 kpc). Approximately 30% of groups in our sample have ambiguous centers with multiple bright or massive galaxies, and these groups show disturbed mass profiles that are not well fit by standard models, suggesting that they are merging systems. We find that halo mass estimates from stacked weak lensing can be biased low by 5-30% if inaccurate centers are used and the issue of miscentering is not addressed.

We report the first weak-lensing detection of a large-scale filament funneling matter onto the core of the massive galaxy cluster MACSJ0717.5+3745. Our analysis is based on a mosaic of 18 multi-passband images obtained with ACS aboard the HST, covering an area of \sim 10x20 arcmin^2. We use a weak-lensing pipeline developed for the COSMOS survey, modified for the analysis of galaxy clusters, to produce a weak-lensing catalogue. A mass map is then computed by applying a weak-gravitational-lensing multi-scale reconstruction technique designed to describe irregular mass distributions such as the one investigated here. We test the resulting mass map by comparing the mass distribution inferred for the cluster core with the one derived from strong-lensing constraints and find excellent agreement. The filament is detected within the 3 sigma detection contour of the lensing mass reconstruction, and underlines the importance of filaments for theoretical and numerical models of the mass distribution in the Cosmic Web. We measure the filament's projected length as \sim 4.5 h_{74}^{-1} Mpc, and its mean density as (2.92 \pm 0.66)10^8 h_{74} M_{\odot} kpc^{-2}. Combined with the redshift distribution of galaxies obtained after an extensive spectroscopic follow-up in the area, we can rule out any projection effect resulting from the chance alignment on the sky of unrelated galaxy group-scale structures. Assuming plausible constraints concerning the structure's geometry based on its galaxy velocity field, we construct a 3D model of the large-scale filament. Within this framework, we derive the three-dimensional length of the filament to be 18 h_{74}^{-1} Mpc, and a deprojected density in terms of the critical density of the Universe of (206 \pm 46) \rho_{crit}, a value that lies at the very high end of the range predicted by numerical simulations.

We investigate the gravitational interactions between live stellar disks and their dark matter halos, using LCDM haloes similar in mass to that of the Milky Way taken from the Aquarius Project. We introduce the stellar disks by first allowing the haloes to respond to the influence of a growing rigid disk potential from z = 1.3 to z = 1.0. The rigid potential is then replaced with star particles which evolve self-consistently with the dark matter particles until z = 0.0. Regardless of the initial orientation of the disk, the inner parts of the haloes contract and change from prolate to oblate as the disk grows to its full size. When the disk normal is initially aligned with the major axis of the halo at z=1.3, the length of the major axis contracts and becomes the minor axis by z=1.0. Six out of the eight disks in our main set of simulations form bars, and five of the six bars experience a buckling instability that results in a sudden jump in the vertical stellar velocity dispersion and an accompanying drop in the m=2 Fourier amplitude of the disk surface density. The bars are not destroyed by the buckling but continue to grow until the present day. Bars are largely absent when the disk mass is reduced by a factor of two or more; the relative disk-to-halo mass is therefore a primary factor in bar formation and evolution. A subset of the disks is warped at the outskirts and contains prominent non-coplanar material with a ring-like structure. Many disks reorient by large angles between z=1 and z=0, following a coherent reorientation of their inner haloes. Larger reorientations produce more strongly warped disks, suggesting a tight link between the two phenomena. The origins of bars and warps appear independent: some disks with strong bars show no disturbances at the outskirts, while the disks with the weakest bars show severe warps.

We present stellar kinematics and orbit superposition models for the central regions of four Brightest Cluster Galaxies (BCGs), based upon integral-field spectroscopy at Gemini, Keck, and McDonald Observatories. Our integral-field data span radii from < 100 pc to tens of kpc. We report black hole masses, M_BH, of 2.1 +/- 1.6 x 10^10 M_Sun for NGC 4889, 9.7 + 3.0 - 2.6 x 10^9 M_Sun for NGC 3842, and 1.3 + 0.5 - 0.4 x 10^9 M_Sun for NGC 7768. For NGC 2832 we report an upper limit of M_BH < 9 x 10^9 M_Sun. Stellar orbits near the center of each galaxy are tangentially biased, on comparable spatial scales to the galaxies' photometric cores. We find possible photometric and kinematic evidence for an eccentric torus of stars in NGC 4889, with a radius of nearly 1 kpc. We compare our measurements of M_BH to the predicted black hole masses from various fits to the relations between M_BH and stellar velocity dispersion, luminosity, or stellar mass. The black holes in NGC 4889 and NGC 3842 are significantly more massive than all dispersion-based predictions and most luminosity-based predictions. The black hole in NGC 7768 is consistent with a broader range of predictions.

We present results of a strong-lensing analysis of MACS J0717.5+3745 (hereafter MACS J0717), an extremely X-ray luminous galaxy cluster at z=0.55. Observations at different wavelengths reveal a complex and dynamically very active cluster, whose core is connected to a large scale filament extended over several Mpc. Using multi-passband imaging data obtained with the Hubble Space Telescope Advanced Camera for Surveys (ACS), we identify 15 multiply imaged systems across the full field of view of ACS, five of which we confirmed spectroscopically in groundbased follow-up observations with the Keck telescope. We use these multiply imaged systems to constrain a parametric model of the mass distribution in the cluster core, employing a new parallelized version of the Lenstool software. The main result is that the most probable description of the mass distribution comprises four cluster-scale dark matter haloes. The total mass distribution follows the light distribution but strongly deviates from the distribution of the intra-cluster gas as traced by the X-ray surface brightness. This confirms the complex morphology proposed by previous studies. We interpret this segregation of collisional and collisionless matter as strong evidence of multiple mergers and ongoing dynamical activity. MACS J0717 thus constitutes one of the most disturbed clusters presently known and, featuring a projected mass within the ACS field of view (R=150"=960 kpc) of 2.11 +- 0.23 x 10^{15} M_{sun}, the system is also one of the most massive known.

We present an analysis of the structures and dynamics of the merging cluster Abell~1201, which has two sloshing cold fronts around a cooling core, and an offset gas core approximately 500kpc northwest of the center. New Chandra and XMM-Newton data reveal a region of enhanced brightness east of the offset core, with breaks in surface brightness along its boundary to the north and east. This is interpreted as a tail of gas stripped from the offset core. Gas in the offset core and the tail is distinguished from other gas at the same distance from the cluster center chiefly by having higher density, hence lower entropy. In addition, the offset core shows marginally lower temperature and metallicity than the surrounding area. The metallicity in the cool core is high and there is an abrupt drop in metallicity across the southern cold front. We interpret the observed properties of the system, including the placement of the cold fronts, the offset core and its tail in terms of a simple merger scenario. The offset core is the remnant of a merging subcluster, which first passed pericenter southeast of the center of the primary cluster and is now close to its second pericenter passage, moving at ~1000 km/s. Sloshing excited by the merger gave rise to the two cold fronts and the disposition of the cold fronts reveals that we view the merger from close to the plane of the orbit of the offset core.

We study the likelihood of energy condition violations in the history of the Universe. Our method is based on a set of functions that characterize energy condition violation. FLRW cosmological models are built around these "indication functions". By computing the Fisher matrix of model parameters using type Ia supernova and Hubble parameter data, we extract the principal modes of these functions' redshift evolution. These modes allow us to obtain general reconstructions of energy condition violation history independent of the dark energy model. We find that the data suggest a history of strong energy condition violation, but the null and dominant energy conditions are likely to be fulfilled. Implications for dark energy models are discussed.

Observational work conducted over the last few decades indicates that all massive galaxies have supermassive black holes at their centres. Although the luminosities and brightness fluctuations of quasars in the early Universe suggest that some are powered by black holes with masses greater than 10 billion solar masses, the remnants of these objects have not been found in the nearby Universe. The giant elliptical galaxy Messier 87 hosts the hitherto most massive known black hole, which has a mass of 6.3 billion solar masses. Here we report that NGC 3842, the brightest galaxy in a cluster at a distance from Earth of 98 megaparsecs, has a central black hole with a mass of 9.7 billion solar masses, and that a black hole of comparable or greater mass is present in NGC 4889, the brightest galaxy in the Coma cluster (at a distance of 103 megaparsecs). These two black holes are significantly more massive than predicted by linearly extrapolating the widely-used correlations between black hole mass and the stellar velocity dispersion or bulge luminosity of the host galaxy. Although these correlations remain useful for predicting black hole masses in less massive elliptical galaxies, our measurements suggest that different evolutionary processes influence the growth of the largest galaxies and their black holes.

We present hydrodynamical simulations of major mergers of galaxies and study the effects of winds produced by active galactic nuclei (AGN) on interstellar gas in the AGN's host galaxy. We consider winds with initial velocities ~ 10,000 km/s and an initial momentum (energy) flux of ~ tau_w L/c (~ 0.01 tau_w L), with tau_w ~ 1-10. The AGN wind sweeps up and shock heats the surrounding interstellar gas, leading to a galaxy-scale outflow with velocities ~ 1000 km/s, peak mass outflow rates comparable to the star formation rate, and a total ejected gas mass ~ 3 x 10^9 M_sun. Large momentum fluxes, tau_w > 3, are required for the AGN-driven galactic outflow to suppress star formation and accretion in the black hole's host galaxy. Less powerful AGN winds (tau_w < 3) still produce a modest galaxy-scale outflow, but the outflow has little global effect on the ambient interstellar gas. We argue that this mechanism of AGN feedback can plausibly produce the high velocity outflows observed in post-starburst galaxies and the massive molecular and atomic outflows observed in local ultra-luminous infrared galaxies. Moreover, the outflows from local ultra-luminous infrared galaxies are inferred to have tau_w ~ 10, comparable to what we find is required for AGN winds to regulate the growth of black holes and set the M_BH-sigma relation. We conclude by discussing theoretical mechanisms that can lead to AGN wind mass-loading and momentum/energy fluxes large enough to have a significant impact on galaxy formation.

We examine atmospheric heating by radio active galactic nuclei (AGN) in distant X-ray clusters by cross correlating clusters selected from the 400 Square Degree (400SD) X-ray Cluster survey with radio sources in the NRAO VLA Sky Survey. Roughly 30% of the clusters show radio emission above a flux threshold of 3 mJy within a projected radius of 250 kpc. The radio emission is presumably associated with the brightest cluster galaxy. The mechanical jet power for each radio source was determined using scaling relations between radio power and cavity (mechanical) power determined for nearby clusters, groups, and galaxies with hot atmospheres containing X-ray cavities. The average jet power of the central radio AGN is approximately $2\times 10^{44}$\ergs. We find no significant correlation between radio power, hence mechanical jet power, and the X-ray luminosities of clusters in the redshift range 0.1 -- 0.6. This implies that the mechanical heating rate per particle is higher in lower mass, lower X-ray luminosity clusters. The jet power averaged over the sample corresponds to an atmospheric heating of approximately 0.2 keV per particle within R$_{500}$. Assuming the current AGN heating rate does not evolve but remains constant to redshifts of 2, the heating rate per particle would rise by a factor of two. We find that the energy injected from radio AGN contribute substantially to the excess entropy in hot atmospheres needed to break self-similarity in cluster scaling relations. The detection frequency of radio AGN is inconsistent with the presence of strong cooling flows in 400SD clusters, but does not exclude weak cooling flows. It is unclear whether central AGN in 400SD clusters are maintained by feedback at the base of a cooling flow. Atmospheric heating by radio AGN may retard the development of strong cooling flows at early epochs.

We use a suite of cosmological hydrodynamic simulations to quantify the accretion rates of baryons into dark matter halos and the resulting baryon mass fractions, as a function of halo mass, redshift, and baryon type (including cold and hot gas). We find that the net baryonic accretion rates through the virial radius are sensitive to galactic outflows and explore a range of outflow parameters to illustrate the effects. We show that the cold gas accretion rate is in general not a simple universal factor of the dark matter accretion rate, and that galactic winds can cause star formation rates to deviate significantly from the external gas accretion rates, both via gas ejection and re-accretion. Furthermore, galactic winds can inject enough energy and momentum in the surrounding medium to slow down accretion altogether, especially in low-mass halos and at low redshift. By resolving the accretion rates versus radius from the halo centers, we show how cold streams penetrate the hot atmospheres of massive halos at z>2, but gradually disappear at lower redshift. The total baryon mass fraction is also strongly suppressed by outflows in low-mass halos, but is nearly universal in the absence of feedback in halos above the UV background suppression scale. The transition halo mass, at which the gas mass in halos is equal for the cold and hot components, is roughly constant at ~10^11.5 Msun and does not depend sensitively on the wind prescription. We provide simple fitting formulae for the cold gas accretion rate into halos in the no-wind case. Finally, we show that cold accretion is broadly consistent with driving the bulk of the highly star-forming galaxies observed at z~2, but that the more intense star formers likely sample the high end of the accretion rate distribution, and may be additionally fueled by a combination of gas recycling, gas re-accretion, hot mode cooling, and mergers.

We use simulated Hubble parameter data in the redshift range 0 \leq z \leq 2 to explore the role and power of observational H(z) data in constraining cosmological parameters of the {\Lambda}CDM model. The error model of the simulated data is empirically constructed from available measurements and scales linearly as z increases. By comparing the median figures of merit calculated from simulated datasets with that of current type Ia supernova data, we find that as many as 64 further independent measurements of H(z) are needed to match the parameter constraining power of SNIa. If the error of H(z) could be lowered to 3%, the same number of future measurements would be needed, but then the redshift coverage would only be required to reach z = 1. We also show that accurate measurements of the Hubble constant H_0 can be used as priors to increase the H(z) data's figure of merit.

This paper is a review on the observational Hubble parameter data that have gained increasing attention in recent years for their illuminating power on the dark side of the universe --- the dark matter, dark energy, and the dark age. Currently, there are two major methods of independent observational H(z) measurement, which we summarize as the "differential age method" and the "radial BAO size method". Starting with fundamental cosmological notions such as the spacetime coordinates in an expanding universe, we present the basic principles behind the two methods. We further review the two methods in greater detail, including the source of errors. We show how the observational H(z) data presents itself as a useful tool in the study of cosmological models and parameter constraint, and we also discuss several issues associated with their applications. Finally, we point the reader to a future prospect of upcoming observation programs that will lead to some major improvements in the quality of observational H(z) data.

We present the first direct measurement of the central black hole mass, M_BH, in NGC 6086, the Brightest Cluster Galaxy (BCG) in Abell 2162. Our investigation demonstrates for the first time that stellar dynamical measurements of M_BH in BCGs are possible beyond the nearest few galaxy clusters. We observed NGC 6086 with laser guide star adaptive optics and the integral-field spectrograph (IFS) OSIRIS at the W.M. Keck Observatory, and with the seeing-limited IFS GMOS-N at Gemini Observatory North. We combined the two IFS data sets with existing major-axis kinematics, and used axisymmetric stellar orbit models to determine M_BH and the R-band stellar mass-to-light ratio, M*/L_R. We find M_BH = 3.6(+1.7)(-1.1) x 10^9 M_Sun and M*/L_R = 4.6(+0.3)(-0.7) M_Sun/L_Sun (68% confidence), from models using the most massive dark matter halo allowed within the gravitational potential of the host cluster. Models fitting only IFS data confirm M_BH ~ 3 x 10^9 M_Sun and M*/L_R ~ 4 M_Sun/L_Sun, with weak dependence on the dark matter halo structure. When data out to 19 kpc are included, the unrealistic omission of dark matter causes the best-fit black hole mass to decrease dramatically, to 0.6 x 10^9 M_Sun, and the best-fit stellar mass-to-light ratio to increase to 6.7 M_Sun/L_Sun. The latter value is at further odds with stellar population studies favoring M*/L ~ 2 M_Sun/L_Sun,R. Biases from dark matter omission could extend to dynamical models of other galaxies with central stellar cores, and new measurements of M_BH from models with dark matter could steepen the empirical scaling relationships between black holes and their host galaxies.

We propose two efficient numerical methods of evaluating the luminosity distance in the spatially flat {\Lambda}CDM universe. The first method is based on the Carlson symmetric form of elliptic integrals, which is highly accurate and can replace numerical quadratures. The second method, using a modified version of Hermite interpolation, is less accurate but involves only basic numerical operations and can be easily implemented. We compare our methods with other numerical approximation schemes and explore their respective features and limitations. Possible extensions of these methods to other cosmological models are also discussed.

The excursion set theory based on spherical or ellipsoidal gravitational collapse provides an elegant analytic framework for calculating the mass function and the large-scale bias of dark matter haloes. This theory assumes that the perturbed density field evolves stochastically with the smoothing scale and exhibits Markovian random walks in the presence of a density barrier. Here we derive an analytic expression for the halo bias in a new theoretical model that incorporates non-Markovian extension of the excursion set theory with a stochastic barrier. This model allows us to handle non-Markovian random walks and to calculate perturbativly these corrections to the standard Markovian predictions for the halo mass function and halo bias. Our model contains only two parameters: kappa, which parameterizes the degree of non-Markovianity and whose exact value depends on the shape of the filter function used to smooth the density field, and a, which parameterizes the degree of stochasticity of the barrier. Appropriate choices of kappa and a in our new model can lead to a closer match to both the halo mass function and halo bias in the latest N-body simulations than the standard excursion set theory.

Different methodologies lead to order-of-magnitude variations in predicted galaxy merger rates. We examine and quantify the dominant uncertainties. Different halo merger rates and subhalo 'destruction' rates agree to within a factor ~2 given proper care in definitions. If however (sub)halo masses are not appropriately defined or are under-resolved, the major merger rate can be dramatically suppressed. The dominant differences in galaxy merger rates owe to baryonic physics. Hydrodynamic simulations without feedback and older models that do not agree with the observed galaxy mass function propagate factor ~5 bias in the resulting merger rates. However, if the model matches the galaxy mass function, properties of central galaxies are sufficiently converged to give small differences in merger rates. But variations in baryonic physics of satellites also have dramatic effects. The known problem of satellite 'over-quenching' in most semi-analytic models (SAMs), whereby SAM satellites are too efficiently stripped of gas, could lead to order-of-magnitude under-estimates of merger rates for low-mass, gas-rich galaxies. Fixing the satellite properties to observations tends to predict higher merger rates, but with factor ~2 empirical uncertainties. Choice of mass ratio definition matters: at low masses, most true major mergers (in baryonic/dynamical galaxy mass) will appear to be minor mergers in their stellar or luminosity mass ratio. Observations and models using these criteria may underestimate major merger rates by factors ~5. Orbital parameters and gas fractions also introduce factor ~3 differences in amount of bulge formed by mergers, even for fixed mass ratio encounters.

In the standard Friedmann-Lemaitre-Robertson-Walker (FLRW) cosmological model, the energy conditions provides model-independent bounds on the behavior of the distance modulus. However, this method can not provide us the detailed information about the violation between the energy conditions and the observation. In this paper, we present an extended analysis of the energy conditions based upon the entropy density of the universe. On the one hand, we find that these conditions imply that entropy density $s$ depends on Hubble parameter H(z). On the other hand, we compare the theoretical entropy density from the conservation law of energy-momentum tensor with that from the energy conditions using the observational Hubble parameter. We find that the two kinds of entropy density are in agreement, only when the present-day entropy density satisfies 0.0222 <= s_0 <= 0.7888. We also obtain that the strong energy condition (SEC) accords with the first law of thermodynamics in the redshift range z < 2.7, the null energy condition (NEC) at z<3.2, and the dominant energy condition (DEC) at z > 2.6. In addition, the energy conditions gives the deceleration parameter 0 <= q(z) <= 2, which is in a predicament of the accelerated expansion of the universe. In particular, the NEC suggests q(z) >= 5/3.

We use an 18' x 9' mosaic of ACS images covering the entire large-scale structure around the X-ray luminous cluster MACSJ0717.5 (z=0.545) to study the morphology of galaxies at the cluster redshift. We find the global fraction of morphological types of galaxies to be consistent with results in the literature. Interestingly, we find that the fraction of S0s also correlates with local galaxy density, in contrast to the findings of a study of the cores of 10 clusters at similar redshift by Dressler et al. We suggest that this apparent inconsistency is due to differences in the spatial coverage around clusters, which is supported by the fact that the correlation disappears for S0s within a radius of 0.6R_200 of MACSJ0717. We interpret this result as evidence of the morphology-density relation being caused by a combination of morphological transformation triggered by galaxy-galaxy interactions, and effects related to the formation and evolution of large-scale structure. In environments of low to intermediate density, where galaxy-galaxy interactions are frequent and efficient, the observed pronounced morphology-density relation for S0s reflects the density dependence of the interaction cross section. In clusters, however, the correlation disappears for S0s, as the much higher galaxy velocities in clusters not only lower the interaction cross section, but also cause a spatial redistribution of galaxies that all but destroys such a correlation. This argument does not hold for elliptical galaxies in clusters which, having formed much earlier, have settled into the large-scale cluster potential; hence the morphology-density relation for cluster ellipticals may reflect primarily the state of advanced dynamical relaxation of this population within the cluster rather than a causal link to the environment responsible for the morphological transformation of galaxies.

Neutron-rich isotopes with masses near that of iron are produced in type Ia and II supernovae. Traces of such nucleosynthesis are found in primitive meteorites in the form of variations in the isotopic abundance of 54Cr, the most neutron-rich stable isotope of chromium. The hosts of these isotopic anomalies must be presolar grains that condensed in the outflows of supernovae, offering the opportunity to study the nucleosynthesis of iron-peak nuclei in ways that complement spectroscopic observations and can inform models of stellar evolution. However, despite almost two decades of extensive search, the carrier of 54Cr anomalies is still unknown, presumably because it is fine-grained and is chemically labile. Here we identify in the primitive meteorite Orgueil the carrier of 54Cr-anomalies as nanoparticles, most likely spinels that show large enrichments in 54Cr relative to solar composition (54Cr/52Cr ratio >3.6xsolar). Such large enrichments in 54Cr can only be produced in supernovae. The mineralogy of the grains supports condensation in the O/Ne-O/C zones of a type II supernova, although a type Ia origin cannot be excluded. We suggest that planetary materials incorporated different amounts of these nanoparticles, possibly due to late injection by a nearby supernova that also delivered 26Al and 60Fe to the solar system. This idea explains why the relative abundance of 54Cr and other neutron-rich isotopes vary between planets and meteorites. We anticipate that future isotopic studies of the grains identified here will shed new light on the birth of the solar system and the conditions insupernovae.

The increasing abundance of passive "red-sequence" galaxies since z=1-2 is mirrored by a coincident rise in the number of galaxies with spheroidal morphologies. In this paper, however, we show that in detail the correspondence between galaxy morphology and color is not perfect, providing insight into the physical origin of this evolution. Using the COSMOS survey, we study a significant population of red sequence galaxies with disk-like morphologies. These passive disks typically have Sa-Sb morphological types with large bulges, but they are not confined to dense environments. They represent nearly one-half of all red-sequence galaxies and dominate at lower masses (log Mstar < 10) where they are increasingly disk-dominated. As a function of time, the abundance of passive disks with log Mstar < 11 increases, but not as fast as red-sequence spheroidals in the same mass range. At higher mass, the passive disk population has declined since z~1, likely because they transform into spheroidals. We estimate that as much as 60% of galaxies transitioning onto the red sequence evolve through a passive disk phase. The origin of passive disks therefore has broad implications for understanding how star formation shuts down. Because passive disks tend to be more bulge-dominated than their star-forming counterparts, a simple fading of blue disks does not fully explain their origin. We explore several more sophisticated explanations, including environmental effects, internal stabilization, and disk regrowth during gas-rich mergers. While previous work has sought to explain color and morphological transformations with a single process, these observations open the way to new insight by highlighting the fact that galaxy evolution may actually proceed through several separate stages.

Supernova (SN) rates are potentially powerful diagnostics of metal enrichment and SN physics, particularly in galaxy clusters with their deep, metal-retaining potentials and relatively simple star-formation histories. We have carried out a survey for supernovae (SNe) in galaxy clusters, at a redshift range 0.5<z<0.9, using the Advanced Camera for Surveys (ACS) on the Hubble Space Telescope. We reimaged a sample of 15 clusters that were previously imaged by ACS, thus obtaining two to three epochs per cluster, in which we discovered five likely cluster SNe, six possible cluster SNe Ia, two hostless SN candidates, and several background and foreground events. Keck spectra of the host galaxies were obtained to establish cluster membership. We conducted detailed efficiency simulations, and measured the stellar luminosities of the clusters using Subaru images. We derive a cluster SN rate of 0.35 SNuB +0.17/-0.12 (statistical) \pm0.13 (classification) \pm0.01 (systematic) [where SNuB = SNe (100 yr 10^10 L_B_sun)^-1] and 0.112 SNuM +0.055/-0.039 (statistical) \pm0.042 (classification) \pm0.005 (systematic) [where SNuM = SNe (100 yr 10^10 M_sun)^-1]. As in previous measurements of cluster SN rates, the uncertainties are dominated by small-number statistics. The SN rate in this redshift bin is consistent with the SN rate in clusters at lower redshifts (to within the uncertainties), and shows that there is, at most, only a slight increase of cluster SN rate with increasing redshift. The low and fairly constant SN Ia rate out to z~1 implies that the bulk of the iron mass in clusters was already in place by z~1. The recently observed doubling of iron abundances in the intracluster medium between z=1 and 0, if real, is likely the result of redistribution of existing iron, rather than new production of iron.

We study the growth of massive black holes (BH) in galaxies using smoothed particle hydrodynamic simulations of major galaxy mergers with new implementations of BH accretion and feedback. The effect of BH accretion on gas in its host galaxy is modeled by depositing momentum at a rate ~ tau L/c into the ambient gas, where L is the luminosity produced by accretion onto the BH and tau is the wavelength-averaged optical depth of the galactic nucleus to the AGN's radiation (a free parameter of our model). The accretion rate onto the BH is relatively independent of our subgrid accretion model and is instead determined by the BH's dynamical impact on its host galaxy: BH accretion is thus self-regulated rather than `supply limited.' We show that the final BH mass and total stellar mass formed during a merger are more robust predictions of the simulations than the time dependence of the star formation rate or BH accretion rate. In particular, the latter depend on the assumed interstellar medium physics, which determines when and where the gas fragments to form star clusters; this in turn affects the fuel available for further star formation and BH growth. Simulations over a factor of ~ 30 in galaxy mass are consistent with the observed M_BH-sigma relation for a mean optical depth of tau ~ 25. This requires that most BH growth occur when the galactic nucleus is optically thick to far-infrared radiation, consistent with the hypothesized connection between ultra-luminous infrared galaxies and quasars. We find tentative evidence for a shallower M_BH-sigma relation in the lowest mass galaxies, sigma < 100 km/s. Our results demonstrate that feedback-regulated BH growth and consistency with the observed M_BH-sigma relation do not require that BH feedback terminate star formation in massive galaxies or unbind large quantities of cold gas.

We perform hydrodynamical simulations of major galaxy mergers using new methods for calculating the growth of massive black holes (BH) in galactic nuclei and their impact on the surrounding galaxy. We model BH growth by including a subgrid model for accretion produced by angular momentum transport on unresolved scales. The impact of the BHs radiation on surrounding gas is approximated by depositing momentum into the ambient gas, which produces an outward force away from the BH. We argue that these phenomenological models for BH growth and feedback better approximate the interaction between the BH and dense gas in galaxies than previous models. We show that this physics leads to self-regulated black hole growth: during the peak of activity, the accretion rate onto the BH is largely determined by the physics of BH feedback, not the subgrid accretion model. The BH significantly modifies the gas dynamics in the galactic nucleus (< 300 pc), but does not generate large-scale galactic outflows. Integrated over an entire galaxy merger, BH feedback has little effect on the total number of stars formed, but is crucial for setting the BHs mass.

Using the observations of the star formation rate and HI densities to z ~ 4, with measurements of the Molecular Gas Depletion Rate (MGDR) and local density of H_2 at z = 0, we derive the history of the gas consumption by star formation to z ~ 4. We find that closed-box models in which H_2 is not replenished by HI require improbably large increases in rho(H_2) and a decrease in the MGDR with lookback time that is inconsistent with observations. Allowing the H_2 used in star formation to be replenished by HI does not alleviate the problem because observations show that there is very little evolution of rho(HI) from z = 0 to z = 4. We show that to be consistent with observational constraints, star formation on cosmic timescales must be fueled by intergalactic ionized gas, which may come from either accretion of gas through cold (but ionized) flows or from ionized gas associated with accretion of dark matter halos. We constrain the rate at which the extraglactic ionized gas must be converted into HI and ultimately into H_2. The ionized gas inflow rate roughly traces the SFRD: about 1 - 2 x 10^8 M_sun Gyr^-1 Mpc^-3 from z ~ 1 - 4, decreasing by about an order of magnitude from z=1 to z=0 with details depending largely on MGDR(t). All models considered require the volume averaged density of rho(H_2) to increase by a factor of 1.5 - 10 to z ~ 1.5 over the currently measured value. Because the molecular gas must reside in galaxies, it implies that galaxies at high z must, on average, be more molecule rich than they are at the present epoch, which is consistent with observations. These quantitative results, derived solely from observations, agree well with cosmological simulations.

We construct merger trees of dark matter haloes and quantify their merger rates and mass growth rates using the joint dataset from the Millennium and Millennium-II simulations. The finer resolution of the Millennium-II Simulation has allowed us to extend our earlier analysis of halo merger statistics to an unprecedentedly wide range of descendant halo mass (10^10 < M0 < 10^15 Msun), progenitor mass ratio (10^-5 < xi < 1), and redshift (0 < z < 15). We update our earlier fitting form for the mean merger rate per halo as a function of M_0, xi, and z. The overall behavior of this quantity is unchanged: the rate per unit redshift is nearly independent of z out to z~15; the dependence on halo mass is weak (M0^0.13); and it is nearly a power law in the progenitor mass ratio (xi^-2). We also present a simple and accurate fitting formula for the mean mass growth rate of haloes as a function of mass and redshift. This mean rate is 46 Msun/yr for 10^12 Msun haloes at z=0, and it increases with mass as M^{1.1} and with redshift as (1+z)^2.5 (for z > 1). When the fit for the mean mass growth rate is integrated over a halo's history, we find excellent match to the mean mass assembly histories of the simulated haloes. By combining merger rates and mass assembly histories, we present results for the number of mergers over a halo's history and the statistics of the redshift of the last major merger.

We present a statistically complete sample of very X-ray luminous galaxy clusters detected in the MAssive Cluster Survey (MACS). This second MACS release comprises all 34 MACS clusters with nominal X-ray fluxes in excess of 2x10^(-12) erg/s/cm^2 (0.1-2.4 keV) in the ROSAT Bright Source Catalogue; two thirds of them are new discoveries. Extending over the redshift range from 0.3 to 0.5, this subset complements the complete sample of the 12 most distant MACS clusters (z>0.5) published in 2007 and further exemplifies the efficacy of X-ray selection for the compilation of samples of intrinsically massive galaxy clusters. Extensive follow-up observations with Chandra/ACIS led to three additional MACS cluster candidates being eliminated as (predominantly) X-ray point sources. For another four clusters --- which, however, remain in our sample of 34 --- the point-source contamination was found to be about 50%. The median X-ray luminosity of 1.3x10^45 erg/s (0.1-2.4 keV, Chandra, within r_500) of the clusters in this subsample demonstrates the power of the MACS survey strategy to find the most extreme and rarest clusters out to significant redshift. A comparison of the optical and X-ray data for all clusters in this release finds a wide range of morphologies with no obvious bias in favour of either relaxed or merging systems.

We present results of an extensive morphological, spectroscopic, and photometric study of the galaxy population of MACS J0025.4$-$1225 (z=0.586), a major cluster merger with clear segregation of dark and luminous matter, to examine the impact of mergers on galaxy evolution. Based on 436 galaxy spectra obtained with Keck DEIMOS, we identified 212 cluster members within 4 Mpc of the cluster centre, and classified them using three spectroscopic types; we find 111 absorption-line, 90 emission-line (including 23 e(a) and 11 e(b)), and 6 E+A galaxies. The fraction of absorption(emission)-line galaxies is a monotonically increasing(decreasing) function of both projected galaxy density and radial distance to the cluster center. More importantly, the 6 observed E+A cluster members are all located between the dark-matter peaks of the cluster and within ~0.3Mpc radius of the X-ray flux peak, unlike the E+A galaxies in other intermediate-redshift clusters which are usually found to avoid the core region. In addition, we use Hubble Space Telescope imaging to classify cluster members according to morphological type. We find the global fraction of spiral and lenticular galaxies in MACS J0025 to be among the highest observed to date in clusters at z>0.5. The observed E+A galaxies are found to be of lenticular type with Sersic indices of ~2, boosting the local fraction of S0 to 70 per cent between the dark-matter peaks. Combing the results of our analysis of the spatial distribution, morphology, and spectroscopic features of the galaxy population, we propose that the starburst phase of these E+A galaxies was both initiated and terminated during the first core-passage about 0.5--1Gyr ago, and that their morphology has already been transformed into S0 due to ram pressure and/or tidal forces near the cluster core. By contrast, ongoing starbursts are observed predominantly in infalling galaxies, and thus appears to be unrelated to the cluster merger.

We have measured the sub-milli-arcsecond structure of 274 extragalactic sources at 24 and 43 GHz in order to assess their astrometric suitability for use in a high frequency celestial reference frame (CRF). Ten sessions of observations with the Very Long Baseline Array have been conducted over the course of $\sim$5 years, with a total of 1339 images produced for the 274 sources. There are several quantities that can be used to characterize the impact of intrinsic source structure on astrometric observations including the source flux density, the flux density variability, the source structure index, the source compactness, and the compactness variability. A detailed analysis of these imaging quantities shows that (1) our selection of compact sources from 8.4 GHz catalogs yielded sources with flux densities, averaged over the sessions in which each source was observed, of about 1 Jy at both 24 and 43 GHz, (2) on average the source flux densities at 24 GHz varied by 20%-25% relative to their mean values, with variations in the session-to-session flux density scale being less than 10%, (3) sources were found to be more compact with less intrinsic structure at higher frequencies, and (4) variations of the core radio emission relative to the total flux density of the source are less than 8% on average at 24 GHz. We conclude that the reduction in the effects due to source structure gained by observing at higher frequencies will result in an improved CRF and a pool of high-quality fiducial reference points for use in spacecraft navigation over the next decade.

We present [OIII 5007A] observations of the star forming galaxy HDF-BMZ1299 (z=1.598) using Keck Observatory's Adaptive Optics system with the near-infrared integral field spectrograph OSIRIS. Using previous Halpha and [NII] measurements of the same source, we are able for the first time to use spatially resolved observations to place a high-redshift galaxy's substructure on a traditional HII diagnostic diagram. We find that HDF-BMZ1299's spatially concentrated nebular ratios in the central ~1.5 kiloparsec (0."2) are best explained by the presence of an AGN: log([NII]/Halpha)=-0.22+/-0.05 and 2sigma limit of log([OIII]/Hbeta)>0.26. The dominant energy source of this galaxy is star formation, and integrating a single aperture across the galaxy yields nebular ratios that are composite spectra from both AGN and HII regions. The presence of an embedded AGN in HDF-BMZ1299 may suggest a potential contamination in a fraction of other high-redshift star forming galaxies, and we suggest that this may be a source of the "elevated" nebular ratios previously seen in seeing-limited metallicity studies. HDF-BMZ1299's estimated AGN luminosity is L_Halpha = 3.7e41 erg/s and L_[OIII] = 5.8e41 erg/s, making it one of the lowest luminosity AGN discovered at this early epoch.

We present results of a gravitational-lensing and optical study of MACS ,J1423.8+2404 (z=0.545, MACS, J1423). Our analysis uses high-resolution images taken with the Hubble Space Telescope in the F555W and F814W passbands, ground based imaging in eight optical and near-infrared filters obtained with Subaru and CFHT, as well as extensive spectroscopic data gathered with the Keck telescopes. At optical wavelengths the cluster exhibits no sign of substructure and is dominated by a cD galaxy that is 2.1 magnitudes (K-band) brighter than the second brightest cluster member, suggesting that MACS, J1423 is close to be fully virialized. Analysis of the redshift distribution of 140 cluster members reveals a Gaussian distribution, mildly disturbed by the presence of a loose galaxy group that may be falling into the cluster along the line of sight. Combining strong-lensing constraints from two spectroscopically confirmed multiple-image systems near the cluster core with a weak-lensing measurement of the gravitational shear on larger scales, we derive a parametric mass model for the mass distribution. All constraints can be satisfied by a uni-modal mass distribution centred on the cD galaxy and exhibiting very little substructure. The derived projected mass of M(<65\arcsec [415 kpc])=(4.3\pm0.6)\times 10^{14} M_sun is about 30% higher than the one derived from X-ray analyses assuming spherical symmetry, suggesting a slightly prolate mass distribution consistent with the optical indication of residual line-of-sight structure. The similarity in shape and excellent alignment of the centroids of the total mass, K-band light, and intra-cluster gas distributions add to the picture of a highly evolved system [ABRIDGED]

We present Advanced Camera for Surveys observations of MACSJ1149.5+2223, an X-ray luminous galaxy cluster at z=0.544 discovered by the Massive Cluster Survey. The data reveal at least seven multiply-imaged galaxies, three of which we have confirmed spectroscopically. One of these is a spectacular face-on spiral galaxy at z=1.491, the four images of which are gravitationally magnified by ~8<mu<~23. We identify this as an L* (M_B=-20.7), disk-dominated (B/T<~0.5) galaxy, forming stars at ~6Msol/yr. We use a robust sample of multiply-imaged galaxies to constrain a parameterized model of the cluster mass distribution. In addition to the main cluster dark matter halo and the bright cluster galaxies, our best model includes three galaxy-group-sized halos. The relative probability of this model is P(N_halo=4)/P(N_halo<4)>=10^12 where N_halo is the number of cluster/group-scale halos. In terms of sheer number of merging cluster/group-scale components, this is the most complex strong-lensing cluster core studied to date. The total cluster mass and fraction of that mass associated with substructures within R<=500kpc, are measured to be M_tot=(6.7+/-0.4)x10^14Msol and f_sub=0.25+/-0.12 respectively. Our model also rules out recent claims of a flat density profile at >~7sigma confidence, thus highlighting the critical importance of spectroscopic redshifts of multiply-imaged galaxies when modeling strong lensing clusters. Overall our results attest to the efficiency of X-ray selection in finding the most powerful cluster lenses, including complicated merging systems.

Dark matter haloes in Lambda CDM simulations grow by mergers with other haloes as well as accretion of "diffuse" non-halo material. We quantify the mass growth rates via these two processes, dM_mer/dt and dM_dif/dt, and their dependence on halo environment using the ~500,000 haloes in the Millennium simulation. Adopting a local mass density parameter as a measure of halo environment, we find the two rates show strong but opposite environmental dependence, with mergers playing an increasingly important role for halo growths in overdense regions and diffuse accretion dominating growth in voids. This behaviour is independent of the mass cuts used to define haloes vs non-haloes. For galaxy-scale haloes, these two opposite correlations largely cancel out, but a weak environmental dependence remains that results in a slightly lower mean total growth rate, and hence an earlier mean formation redshift, for haloes in denser regions. The mean formation redshift of the ~5000 cluster-mass haloes, on the other hand, appears to have no correlation with halo environment. The origin of the positive correlation of dM_mer/dt with local density can be traced to the surrounding mass reservoir outside the haloes, where more progenitor haloes are available in denser regions. The negative correlation of dM_dif/dt with density, however, is not explained by the available diffuse mass in the reservoir, which is in fact larger in denser regions. The non-halo component may therefore be partially comprised of truly diffuse dark matter particles that are dynamically hotter and are accreted at a suppressed rate in denser regions. We also discuss the implications of these results for how to modify the Extended Press-Schechter model of halo growth, which in its original form does not predict environmental dependence.

This abstract has been withdrawn. The proper replacement is arXiv:0909.3306

Traditionally, secular evolution is defined as evolution of systems where the internal growth of structure and instabilities dominates the growth via external drivers (e.g. accretion/mergers). Most study has focused on 'isolated' galaxies, where seed asymmetries may represent realistic cosmological substructure, but subsequent evolution ignores galaxy growth. Large-scale modes in the disk then grow on a timescale of order a disk rotation period (0.1-1 Gyr). If, however, galaxies evolve cosmologically on a shorter timescale, then it may not be appropriate to consider them 'isolated.' We outline simple scalings to ask whether the timescale for secular evolution is shorter than the timescale for cosmological accretion and mergers. This is the case in a narrow, but important range of perturbation amplitudes corresponding to substructure or mode/bar fractional amplitudes 0.01-0.1, a range of interest for observed strong bars and pseudobulges. At smaller amplitudes <<0.1, systems are not isolated: typical disks will grow by accretion at a comparable level over even a single dynamical time. At larger amplitudes >>0.1, the evolution is no longer secular; direct gravitational evolution of the seed swamps the internal disk response. We derive criteria for when disks can be well-approximated as 'isolated' as a function of mass, redshift, and disk stability. The relevant parameter space shrinks at higher mass, higher disk stability, and higher-z as accretion rates increase. Cosmological rates of galaxy evolution also define a maximum bar/mode lifetime of practical interest, of ~0.1/H(z). Longer-lived modes will de-couple from their drivers (if driven) and encounter cosmological effects.

Observations suggest that effective radii of high-z massive spheroids are as much as a factor ~6 smaller than low-z galaxies of comparable mass. Given the apparent absence of low-z counterparts, this has often been interpreted as indicating that the high density, compact red galaxies must be 'puffed up' by some mechanism. We compare the ensemble of high-z observations with large samples of well-observed low-z ellipticals. At the same physical radii, the stellar surface mass densities of low and high-z systems are comparable. Moreover, the abundance of high surface density material at low redshift is comparable to or larger than that observed at z>1-2, consistent with the continuous buildup of spheroids over this time. The entire population of compact, high-z red galaxies may be the progenitors of the high-density cores of present-day ellipticals, with no need for a decrease in stellar density from z=2 to z=0. The primary difference between low and high-z systems is thus the observed low-density material at large radii in low-z spheroids (rather than the high-density material in high-z spheroids). Such low-density material may either (1) assemble at z<2 or (2) be present, but not yet detected, at z>2. Mock observations of low-z massive systems show that the high-z observations do not yet probe sufficiently low surface brightness material to detect the low surface density 'wings' (if present). Thus, if the high-z galaxies resemble the most massive systems today, their inferred effective radii could be under-estimated by factors ~2-4. This difference arises because massive systems at low redshift are not well-fit by single Sersic profiles. We discuss implications of our results for physical models of galaxy evolution.

Transformation of disks into spheroids via mergers is a well-accepted element of galaxy formation models. However, recent simulations have shown that bulge formation is suppressed in increasingly gas-rich mergers. We investigate the global implications of these results in a cosmological framework, using independent approaches: empirical halo-occupation models (where galaxies are populated in halos according to observations) and semi-analytic models. In both, ignoring the effects of gas in mergers leads to the over-production of spheroids: low and intermediate-mass galaxies are predicted to be bulge-dominated (B/T~0.5 at <10^10 M_sun), with almost no bulgeless systems), even if they have avoided major mergers. Including the different physical behavior of gas in mergers immediately leads to a dramatic change: bulge formation is suppressed in low-mass galaxies, observed to be gas-rich (giving B/T~0.1 at <10^10 M_sun, with a number of bulgeless galaxies in good agreement with observations). Simulations and analytic models which neglect the similarity-breaking behavior of gas have difficulty reproducing the strong observed morphology-mass relation. However, the observed dependence of gas fractions on mass, combined with suppression of bulge formation in gas-rich mergers, naturally leads to the observed trends. Discrepancies between observations and models that ignore the role of gas increase with redshift; in models that treat gas properly, galaxies are predicted to be less bulge-dominated at high redshifts, in agreement with the observations. We discuss implications for the global bulge mass density and future observational tests.

We use the extensive catalog of dark matter haloes from the Millennium simulation to investigate the statistics of the mass accretion histories (MAHs) and accretion rates of ~500,000 haloes from redshift z=0 to 6. We find only about 25% of the haloes to have MAHs that are well described by a 1-parameter exponential form. For the rest of the haloes, between 20% (Milky-Way mass) to 50% (cluster mass) experience late-time growth that is steeper than an exponential, whereas the remaining haloes show plateau-ed late-time growth that is shallower than an exponential. The haloes with slower late-time growth tend to reside in denser environments, suggesting that either tidal stripping or the "hotter" dynamics are suppressing the accretion rate of dark matter onto these haloes. These deviations from exponential growth are well fit by introducing a second parameter: M(z) \propto (1+z)^beta exp(-gamma z). The full distribution of beta and gamma as a function of halo mass is provided. From the analytic form of M(z), we obtain a simple formula for the mean accretion rate of dark matter, dM/dt, as a function of redshift and mass. At z=0, this rate is 42 Msun/yr for 1e12 Msun haloes, which corresponds to a mean baryon accretion rate of dMbaryon/dt=7 Msun/yr. This mean rate increases approximately as (1+z)^1.5 at low z and (1+z)^2.5 at high z, reaching dMbaryon/dt = 27, 69, and 140 Msun/yr at z=1, 2, and 3. The specific rate depends on halo mass weakly: dlogM/dt \propto M^0.127. Results for the broad distributions about the mean rates are also discussed.

The latest cosmological N-body simulations find two intriguing properties for dark matter haloes: (1) their radial density profile, rho, is better fit by a form that flattens to a constant at the halo center (the Einasto profile) than the widely-used NFW form; (2) the radial profile of the pseudo-phase-space density, rho/sigma3, on the other hand, continues to be well fit by a power law, as seen in earlier lower-resolution simulations. In this paper we use the Jeans equation to argue that (1) and (2) cannot both be true at all radii. We examine the implied radial dependence of rho/sigma3 over 12 orders of magnitude in radius by solving the Jeans equation for a broad range of input rho and velocity anisotropy beta. Independent of beta, we find that rho/sigma3 is approximately a power law only over the limited range of halo radius resolvable by current simulations (down to ~0.1% of the virial radius), and rho/sigma3 deviates significantly from a power-law below this scale for both the Einasto and NFW rho. The same conclusion also applies to a more general density-velocity relation rho/sigma_D^epsilon. Conversely, when we enforce rho/sigma^3 r^{-eta} as an input, none of the physically allowed rho (occurring for the narrow range 1.8<eta<1.9444) follows the Einasto form. We expect the next generation simulations with better spatial resolution to settle the debate: either the Einasto profile will continue to hold and rho/sigma3 will deviate from a power law, or rho/sigma3 will continue as a power law and rho will deviate from its current parameterizations.

Using CHANDRA, we investigate the spatial temperature distribution of the intracluster medium (ICM) within 700 kpc of the center of the massive merging cluster MACSJ0717.5+3745 at z=0.55. Combining the X-ray evidence with information about the distribution and velocities of the cluster galaxies near the core provides us with a snapshot of the three-dimensional geometry and dynamics of one of the most complex cluster studied to date. We find MACSJ0717.5+3745 to be an active triple merger with ICM temperatures exceeding 20 keV. Although radial velocity information and X-ray/optical offsets indicate that all three mergers proceed along distinctly different directions, the partial alignment of the merger axes points to a common origin in the large-scale filament south-east of the cluster core. Clear decrements in the ICM temperature observed near two of these subclusters identify the respective X-ray surface brightness peaks as remnants of cool cores; the compactness and low temperature of 5.7 keV of one of these features suggest that the respective merger, a high-velocity collision at 3,000 km/s, is still in its very early stages. Looking beyond the triple merger, we find the large-scale filament to not only provide a spatial as well as temporal arrow for the interpretation of the dynamics of the merger events near the cluster core; we also find tantalizing, if circumstantial, evidence for direct, large-scale heating of the ICM by contiguous infall of low-density gas from the filament.

In an earlier paper we quantified the mean merger rate of dark matter haloes as a function of redshift z, descendant halo mass M0, and progenitor halo mass ratio xi using the Millennium simulation of the LCDM cosmology. Here we broaden that study and investigate the dependence of the merger rate of haloes on their surrounding environment. A number of local mass overdensity variables, both including and excluding the halo mass itself, are tested as measures of a halo's environment. The simple functional dependence on z, M0, and xi of the merger rate found in our earlier work is largely preserved in different environments, but we find that the overall amplitude of the merger rate has a strong positive correlation with the environmental densities. For galaxy-mass haloes, we find mergers to occur ~2.5 times more frequently in the densest regions than in voids at both z=0 and higher redshifts. Higher-mass haloes show similar trends. We present a fitting form for this environmental dependence that is a function of both mass and local density and is valid out to z=2. The amplitude of the progenitor (or conditional) mass function shows a similarly strong correlation with local overdensity, suggesting that the extended Press-Schechter model for halo growth needs to be modified to incorporate environmental effects.

We discuss the X-ray and optical properties of the massive galaxy cluster MACSJ1206.2-0847 (z=0.4385), discovered in the Massive Cluster Survey (MACS). Our Chandra observation of the system yields a total X-ray luminosity of 2.4 x 10^45 erg/s (0.1-2.4 keV) and a global gas temperature of (11.6 +/- 0.7) keV, very high values typical of MACS clusters. In both optical and X-ray images MACSJ1206 appears close to relaxed in projection, with a pronounced X-ray peak at the location of the brightest cluster galaxy (BCG); we interpret this feature as the remnant of a cold core. A spectacular giant gravitational arc, 15" in length, bright (V~21) and unusually red (R-K=4.3), is seen 20" west of the BCG; we measure a redshift of z=1.036 for the lensed galaxy. From our HST image of the cluster we identify the giant arc and its counter image as a seven-fold imaged system. An excess of X-ray emission in the direction of the arc coincides with a mild galaxy overdensity and could be the remnant of a minor merger with a group of galaxies. We derive estimates of the total cluster mass as well as of the mass of the cluster core using X-ray, dynamical, and gravitational-lensing techniques. For the mass enclosed by the giant arc (r<119 kpc) our strong-lensing analysis based on HST imaging yields a very high value of 1.1 x 10^14 M_sun, inconsistent with the much lower X-ray estimate of 0.5 x 10^14 M_sun. Similarly, the virial estimate of 4 x 10^15 M_sun for the total cluster mass, derived from multi-object spectroscopy of 38 cluster members, is significantly higher than the corresponding X-ray estimate of 1.7 x 10^15 M_sun. We take the discrepant mass estimates to be indicative of substructure along the line of sight during an ongoing merger event, an interpretation that is supported by the system's very high velocity dispersion of 1580 km/s.

The Yuan-Tseh Lee Array for Microwave Background Anisotropy (AMiBA) is the first interferometer dedicated to studying the cosmic microwave background (CMB) radiation at 3mm wavelength. The choice of 3mm was made to minimize the contributions from foreground synchrotron radiation and Galactic dust emission. The initial configuration of seven 0.6m telescopes mounted on a 6-m hexapod platform was dedicated in October 2006 on Mauna Loa, Hawaii. Scientific operations began with the detection of a number of clusters of galaxies via the thermal Sunyaev-Zel'dovich effect. We compare our data with Subaru weak lensing data in order to study the structure of dark matter. We also compare our data with X-ray data in order to derive the Hubble constant.

We investigate seven Monte Carlo algorithms -- four old and three new -- for constructing merger histories of dark matter halos using the extended Press-Schechter (EPS) formalism based on both the spherical and ellipsoidal collapse models. We compare, side-by-side, the algorithms' abilities at reproducing the analytic EPS conditional (or progenitor) mass function over a broad range of mass and redshift (z=0 to 15). Among the four old algorithms (Lacey & Cole 1993, Kauffmann & White 1993, Somerville & Kolatt 1999, Cole et al 2000), we find that only KW93 produces a progenitor mass function that is consistent with the EPS prediction for all look-back redshifts. The origins of the discrepancies in the other three algorithms are discussed. Our three new algorithms are designed to generate the correct progenitor mass function at each timestep. We show that this is a necessary and sufficient condition for consistency with EPS at any look-back time. We illustrate the differences among the three new algorithms and KW93 by investigating two other conditional statistics: the mass function of the i_{th} most massive progenitors and the mass function for descendants with N_p progenitors.

We present maps of the cosmic large-scale structure around the twelve most distant galaxy clusters from the Massive Cluster Survey (MACS) as traced by the projected surface density of galaxies on the cluster red sequence. Taken with the Suprime-Cam wide-field camera on the Subaru telescope, the images used in this study cover a 27x27 arcmin^2 area around each cluster, corresponding to 10 x 10 Mpc^2 at the median redshift of z = 0.55 of our sample. We directly detect satellite clusters and filaments extending over the full size of our imaging data in the majority of the clusters studied, supporting the picture of mass accretion via infall along filaments suggested by numerical simulations of the growth of clusters and the evolution of large-scale structure. A comparison of the galaxy distribution near the cluster cores with the X-ray surface brightness as observed with Chandra reveals, in several cases, significant offsets between the gas and galaxy distribution, indicative of ongoing merger events. The respective systems are ideally suited for studies of the dynamical properties of gas, galaxies, and dark matter. In addition, the large-scale filaments viewed at high contrast in these MACS clusters are prime targets for the direct detection and study of the warm-hot intergalactic medium (WHIM).

We present the results of a wide-field spectroscopic analysis of the galaxy population of the massive cluster MACSJ0717.5+3745 and the surrounding filamentary structure (z=0.55), as part of our systematic study of the 12 most distant clusters in the MACS sample. Of 1368 galaxies spectroscopically observed in this field, 563 are identified as cluster members; of those, 203 are classified as emission-line galaxies, 260 as absorption-line galaxies, and 17 as E+A galaxies (defined by $\frac{H_{\delta}+H_{\gamma}}{2}>6$\AA and no detection of [OII] and $H_{\beta}$ in emission). The variation of the fraction of emission- and absorption-line galaxies as a function of local projected galaxy density confirms the well-known morphology-density relation, and becomes flat at projected galaxy densities less than $\sim 20Mpc^{-2}. Interestingly, 16 out of 17 E+A galaxies lie (in projection) within the ram-pressure stripping radius around the cluster core, which we take to be direct evidence of ram-pressure stripping being the primary mechanism that terminates star-formation in the E+A population of galaxy clusters. This conclusion is supported by the rarity of E+A galaxies in the filament which rules out galaxy mergers as the dominant driver of evolution for E+A galaxies in clusters. In addition, we find the 42 e(a) and 27 e(b) member galaxies, i.e., the dusty-starburst and starburst galaxies respectively, to be spread out across almost the entire study area. Their spatial distribution, which shows a strong preference for the filament region, suggests that starbursts are triggered in relatively low-density environments as galaxies are accreted from the field population.

We present a study of the mass-metallicity (M-Z) relation and H II region physical conditions in a sample of 20 star-forming galaxies at 1.0<z<1.5 drawn from the DEEP2 Galaxy Redshift Survey. Our analysis indicates that the zero point of the M-Z relationship evolves with redshift, in the sense that galaxies at fixed stellar mass become more metal-rich at lower redshift. Measurements of [O III]/H-beta and [N II]/H-alpha emission-line ratios show that, on average, objects in the DEEP2 1.0<z<1.5 sample are significantly offset from the excitation sequence observed in nearby H II regions and SDSS emission-line galaxies. In order to fully understand the causes of this offset, which is also observed in z~2 star-forming galaxies, we examine in detail the small fraction of SDSS galaxies that have similar diagnostic ratios to those of the DEEP2 sample. Some of these galaxies indicate evidence for AGN and/or shock activity, which may give rise to their unusual line ratios, and contribute to Balmer emission lines at the level of ~20%. Others indicate no evidence for AGN or shock excitation, yet are characterized by higher electron densities and temperatures, and therefore, interstellar gas pressures, than typical SDSS star-forming galaxies of similar stellar mass. These anomalous objects also have higher concentrations and star-formation rate surface densities, which are directly connected to higher interstellar pressure. Higher star-formation rate surface densities, interstellar pressures, and HII region ionization parameters may also be common at high-redshift. These effects must be taken into account when using strong-line indicators to understand the evolution of heavy elements in galaxies. When such effects are included, the inferred evolution of the M-Z relation out to z~2 is more significant than previous estimates. (abridged)

Analytic models based on spherical and ellipsoidal gravitational collapse have been used to derive the mass functions of dark matter halos and their progenitors (the conditional mass function). The ellipsoidal model generally provides a better match to simulation results, but there has been no simple analytic expression in this model for the conditional mass function that is accurate for small time steps, a limit that is important for generating halo merger trees and computing halo merger rates. We remedy the situation by deriving accurate analytic formulae for the first-crossing distribution, the conditional mass function, and the halo merger rate in the ellipsoidal collapse model in the limit of small look-back times. We show that our formulae provide a closer match to the Millennium simulation results than those in the spherical collapse model and the ellipsoidal model of Sheth & Tormen (2002).

We construct merger trees from the largest database of dark matter haloes to date provided by the Millennium simulation to quantify the merger rates of haloes over a broad range of descendant halo mass (1e12 < M0 < 1e15 Msun), progenitor mass ratio (1e-3 < xi < 1), and redshift (0 < z < 6). We find the mean merger rate per halo, B/n, to have very simple dependence on M0, xi, and z, and propose a universal fitting form for B/n that is accurate to 10-20%. Overall, B/n depends very weakly on the halo mass (proportional to M0^0.08) and scales as a power law in the progenitor mass ratio (proportional to xi^-2) for minor mergers (xi < 0.1) with a mild upturn for major mergers. As a function of time, we find the merger rate per Gyr to evolve as (1+z)^n with n=2-2.3, while the rate per unit redshift is nearly independent of z. Several tests are performed to assess how our merger rates are affected by, e.g. the time interval between Millennium outputs, binary vs multiple progenitor mergers, and mass conservation and diffuse accretion during mergers. In particular, we find halo fragmentations to be a general issue in merger tree construction from N-body simulations and compare two methods for handling these events. We compare our results with predictions of two analytical models for halo mergers based on the Extended Press-Schechter (EPS) model and the coagulation theory. We find the EPS model to overpredict the major merger rates and underpredict the minor merger rates by up to a factor of a few.

The timescale for galaxies within merging dark matter halos to merge with each other is an important ingredient in galaxy formation models. Accurate estimates of merging timescales are required for predictions of astrophysical quantities such as black hole binary merger rates, the build-up of stellar mass in central galaxies, and the statistical properties of satellite galaxies within dark matter halos. In this paper, we study the merging timescales of extended dark matter halos using N-body simulations. We compare these results to standard estimates based on the Chandrasekhar theory of dynamical friction. We find that these standard predictions for merging timescales, which are often used in semi-analytic galaxy formation models, are systematically shorter than those found in simulations. The discrepancy is approximately a factor of 1.7 for $M_sat/M_host \approx 0.1$ and becomes larger for more disparate satellite-to-host mass ratios, reaching a factor of $\sim 3.3$ for $M_sat/M_host\approx 0.01$. Based on our simulations, we propose a new, easily implementable fitting formula that accurately predicts the timescale for an extended satellite to sink from the virial radius of a host halo down to the halo's center for a wide range of $M_sat/M_host$ and orbits. Including a central bulge in each galaxy changes the merging timescale by $\la 10%$. To highlight one concrete application of our results, we show that merging timescales often used in the literature overestimate the growth of stellar mass by satellite accretion by $\approx 40 %$, with the extra mass gained in low mass ratio mergers.

Whether or not deconfined quark phase exists in neutron star cores and represents the ground state of neutral matter at moderate densities are open questions. We use two realistic effective quark models, the three-flavor Nambu-Jona-Lasinio model and the modified quark-meson coupling model, to describe the neutron star matter. After constructing possible hybrid equations of state (EOSes) with unpaired or color superconducting quark phase, we systematically discuss the observational constraints of neutron stars on the EOSes. It is found that the neutron star with pure quark matter core is unstable and the hadronic phase with hyperons is denied, while hybrid EOSes with two-flavor color superconducting phase or unpaired quark matter phase are both allowed by the tight and most reliable constraints from two stars Ter 5 I and EXO 0748-676. And the hybrid EOS with unpaired quark matter phase is allowed even compared with the tightest constraint from the most massive pulsar star PSR J0751+1807. Therefore, we conclude that the ground state of neutral matter at moderate densities is in deconfined quark phase likely.

We present the statistically complete and cosmologically most relevant subset of the twelve most distant galaxy clusters detected at z>0.5 by the MAssive Cluster Survey (MACS). Ten of these systems are new discoveries; only two (MACSJ0018.5+1626 aka CL0016+1609, and MACSJ0454.1-0300 aka MS0451.6-0305) were previously known. We provide fundamental cluster properties derived from our optical and X-ray follow-up observations as well as the selection function in tabulated form to facilitate cosmological studies using this sample.

We analyze the fundamental plane projections of elliptical galaxies as a function of luminosity, using a sample of approximately 80,000 galaxies drawn from Data Release 4 (DR4) of the Sloan Digital Sky Survey (SDSS). We separate brightest cluster galaxies (BCGs) from our main sample and reanalyze their photometry due to a problem with the default pipeline sky subtraction for BCGs. The observables we consider are effective radius (R_e), velocity dispersion (sigma), dynamical mass (M_dyn ~ R_e sigma2), effective density (sigma2/R_e2), and effective surface brightness (mu_e). With the exception of the L-M_dyn correlation, we find evidence of variations in the slope (i.e. the power-law index) of the fundamental plane projections with luminosity for our normal elliptical galaxy population. In particular, the radius-luminosity and Faber-Jackson relations are steeper at high luminosity relative to low luminosity, and the more luminous ellipticals become progressively less dense and have lower surface brightnesses than lower luminosity ellipticals. These variations can be understood as arising from differing formation histories, with more luminous galaxies having less dissipation. Data from the literature and our reanalysis of BCGs show that BCGs have radius-luminosity and Faber-Jackson relations steeper than the brightest non-BCG ellipticals in our sample, consistent with significant growth of BCGs via dissipationless mergers. The variations in slope we find in the Faber-Jackson relation of non-BCGs are qualitatively similar to that reported in the black hole mass-velocity dispersion (M_BH-sigma) correlation. This similarity is consistent with a roughly constant value of M_BH/M_star over a wide range of early type galaxies, where M_star is the stellar mass.

Minor mergers of galaxies are expected to be common in a hierarchical cosmology such as $\Lambda$CDM and have the potential to significantly affect galactic structure. In this paper we dissect the case-by-case outcome from a set of numerical simulations of a single satellite elliptical galaxy accreting onto a massive elliptical galaxy. We take care to explore cosmologically relevant orbital parameters and to set up realistic initial galaxy models that include all three relevant dynamical components: dark matter halos, stellar bulges, and central massive black holes. The effects of several different parameters are considered, including orbital energy and angular momentum, satellite density and inner density profile, satellite-to-host mass ratio, and presence of a black hole at the center of the host. Black holes play a crucial role in protecting the shallow stellar cores of the hosts, as satellites merging onto a host with a central black hole are more strongly disrupted than those merging onto hosts without black holes. Orbital parameters play an important role in determining the degree of disruption: satellites on less bound or more eccentric orbits are more easily destroyed than those on more bound or more circular orbits as a result of an increased number of pericentric passages and greater cumulative effects of gravitational shocking and tidal stripping. In addition, satellites with densities typical of faint elliptical galaxies are disrupted relatively easily, while denser satellites can survive much better in the tidal field of the host. Over the range of parameters explored, we find that the accretion of a single satellite elliptical galaxy can result in a broad variety of changes, in both signs, in the surface brightness profile and color of the central part of an elliptical galaxy.

Several recent observations suggest that gas-poor (dissipationless) mergers of elliptical galaxies contribute significantly to the build-up of the massive end of the red sequence. We perform a series of major merger simulations to investigate the spatial and velocity structure of the remnants of such mergers. Regardless of orbital energy or angular momentum, we find that the stellar remnants lie on the fundamental plane defined by their progenitors, a result of virial equilibrium with a small tilt due to an increasing central dark matter fraction. However, the locations of merger remnants in the projections of the fundamental plane -- the Faber-Jackson and R_e-M_* relations -- depend strongly on the merger orbit, and the relations steepen significantly from the canonical scalings (L sigma^4 and R_e M_*^0.6) for mergers on radial orbits. Our results imply that the projections of the fundamental plane -- but not necessarily the plane itself -- provide a powerful way of investigating the assembly history of massive elliptical galaxies, including the brightest cluster galaxies at or near the centers of galaxy clusters. We argue that most massive ellipticals are formed by anisotropic merging and that their fundamental plane projections should thus differ noticeably from those of lower mass ellipticals even though they should lie on the same fundamental plane. Current observations are consistent with this conclusion. The steepening in the L-sigma relation for luminous ellipticals may also be reflected in a corresponding steepening in the M_BH-sigma relation for massive black holes.

We study the clustering properties of groups and of galaxies in groups in the DEEP2 Galaxy Redshift Survey dataset at z~1. Four clustering measures are presented: 1) the group correlation function for 460 groups with estimated velocity dispersions of sigma>200 km/s, 2) the galaxy correlation for the full galaxy sample, using a flux-limited sample of 9800 objects between 0.7<z<1.0, 3) the galaxy correlation for galaxies in groups, and 4) the group-galaxy cross-correlation function. Using the observed number density and clustering amplitude of the groups, the estimated minimum group dark matter halo mass is M_min~6 10^12 h^-1 M_Sun for a flat LCDM cosmology. Groups are more clustered than galaxies, with a relative bias of b=1.7 +/-0.04 on scales r_p=0.5-15 Mpc/h. Galaxies in groups are also more clustered than the full galaxy sample, with a scale-dependent relative bias which falls from b~2.5 +/-0.3 at r_p=0.1 Mpc/h to b~1 +/-0.5 at r_p=10 Mpc/h. The correlation functions for all galaxies and galaxies in groups can be fit by a power-law on scales r_p=0.05-20 Mpc/h. We empirically measure the contribution to the projected correlation function for galaxies in groups from a `one-halo' term and a `two-halo' term by counting pairs of galaxies in the same or in different groups. The projected cross-correlation between shows that red galaxies are more centrally concentrated in groups than blue galaxies at z~1. DEEP2 galaxies in groups appear to have a shallower radial distribution than that of mock galaxy catalogs made from N-body simulations, which assume a central galaxy surrounded by satellite galaxies with an NFW profile. We show that the clustering of galaxies in groups can be used to place tighter constraints on the halo model than can be gained from using the usual galaxy correlation function alone.

We present the results of near-infrared spectroscopic observations for a sample of 12 star-forming galaxies at 1.0<z<1.5, drawn from the DEEP2 Galaxy Redshift Survey. H-beta, [OIII], H-alpha, and [NII] emission-line fluxes are measured for these galaxies. Application of the O3N2 and N2 strong-line abundance indicators implies average gas-phase oxygen abundances of 50-80% solar. We find preliminary evidence of luminosity-metallicity (L-Z) and mass-metallicity (M-Z) relationships within our sample, which spans from M_B=-20.3 to -23.1 in rest-frame optical luminosity, and from 4x10^9 to 2x10^11 Msun in stellar mass. At fixed oxygen abundance, these relationships are displaced from the local ones by several magnitudes towards brighter absolute B-band luminosity and more than an order of magnitude towards larger stellar mass. If individual DEEP2 galaxies in our sample follow the observed global evolution in the B-band luminosity function of blue galaxies between z~1 and z~0 they will fade on average by ~1.3 magnitudes in M_B. To fall on local (L-Z) and (M-Z) relationships, these galaxies must increase by a factor of 6-7 in M/L_B between z~1 and z~0, and by factor of two in both stellar mass and metallicity. Emission line diagnostic ratios indicate that the z>1 DEEP2 galaxies in our sample are significantly offset from the excitation sequence observed in nearby HII regions and SDSS emission-line galaxies. This offset implies that physical conditions are different in the HII regions of distant galaxies hosting intense star formation, and may affect the chemical abundances derived from strong-line ratios for such objects (Abridged).

We explore the effect of substructure in dark matter halos on the power spectrum and bispectrum of matter fluctuations and weak lensing shear. By experimenting with substructure in a cosmological N = 512^3 simulation, we find that when a larger fraction of the host halo mass is in subhalos, the resulting power spectrum has less power at 1 < k < 100 h Mpc^{-1} and more power at k > 100 h Mpc^{-1}. We explain this effect using an analytic halo model including subhalos, which shows that the 1 < k < 100 h Mpc^{-1} regime depends sensitively on the radial distribution of subhalo centers while the interior structure of subhalos is important at k > 100 h Mpc^{-1}. The corresponding effect due to substructures on the weak lensing power spectrum is up to about 11% at angular scale l < 10^4. Predicting the nonlinear power spectrum to a few percent accuracy for future surveys would therefore require large cosmological simulations that also have exquisite numerical resolution to model accurately the survivals of dark matter subhalos in the tidal fields of their hosts.

We carry out numerical simulations of dissipationless major mergers of elliptical galaxies using initial galaxy models that consist of a dark matter halo and a stellar bulge with properties consistent with the observed fundamental plane. By varying the density profile of the dark matter halo (standard NFW versus adiabatically contracted NFW), the global stellar to dark matter mass ratio, and the orbit of the merging galaxies, we are able to assess the impact of each of these factors on the structure of the merger remnant. Our results indicate that the properties of the remnant bulge depend primarily on the angular momentum and energy of the orbit; for a cosmologically motivated orbit, the effective radius and velocity dispersion of the remnant bulge remain approximately on the fundamental plane. This indicates that the observed properties of elliptical galaxies are consistent with significant growth via late dissipationless mergers. We also find that the dark matter fraction within the effective radius of our remnants increases after the merger, consistent with the hypothesis that the tilt of the fundamental plane from the virial theorem is due to a varying dark matter fraction as a function of galaxy mass.

We study the effects of ellipticity in lens galaxies and external tidal shear from neighboring objects on the statistics of strong gravitational lenses. For isothermal lens galaxies normalized so that the Einstein radius is independent of ellipticity and shear, ellipticity {\it reduces} the lensing cross section slightly, and shear leaves it unchanged. Ellipticity and shear can significantly enhance the magnification bias, but only if the luminosity function of background sources is steep. Realistic distributions of ellipticity and shear {\it lower} the total optical depth by a few percent for most source luminosity functions, and increase the optical depth only for steep luminosity functions. The boost in the optical depth is noticeable (>5%) only for surveys limited to the brightest quasars (L/L_* > 10). Ellipticity and shear broaden the distribution of lens image separations but do not affect the mean. Ellipticity and shear naturally increase the abundance of quadruple lenses relative to double lenses, especially for steep source luminosity functions, but the effect is not enough (by itself) to explain the observed quadruple-to-double ratio. With such small changes to the optical depth and image separation distribution, ellipticity and shear have a small effect on cosmological constraints from lens statistics: neglecting the two leads to biases of just Delta Omega_M = 0.00 \pm 0.01 and Delta Omega_Lambda = -0.02 \pm 0.01 (where the errorbars represent statistical uncertainties in our calculations).

The current realization of the International Celestial Reference Frame (ICRF) comprises a total of 717 extragalactic radio sources distributed over the entire sky. An observing program has been developed to densify the ICRF in the northern sky using the European VLBI network (EVN) and other radio telescopes in Spitsbergen, Canada and USA. Altogether, 150 new sources selected from the Jodrell Bank-VLA Astrometric Survey were observed during three such EVN+ experiments conducted in 2000, 2002 and 2003. The sources were selected on the basis of their sky location in order to fill the "empty" regions of the frame. A secondary criterion was based on source compactness to limit structural effects in the astrometric measurements. All 150 new sources have been successfully detected and the precision of the estimated coordinates in right ascension and declination is better than 1 milliarcsecond (mas) for most of them. A comparison with the astrometric positions from the Very Long baseline Array Calibrator Survey for 129 common sources indicates agreement within 2 mas for 80% of the sources.

When applied to deep photometric catalogs, the two-point angular correlation function, w(theta), is a sensitive probe of the evolution of galaxy clustering properties. Here we present measurements of w(theta) as a function of I_AB magnitude and (R-I) color to a depth of I_AB=24, using a sample of ~350,000 galaxies covering 5 degrees^2 in total over 5 separate fields. Using redshifts of 3319 galaxies in early DEEP2 Galaxy Redshift Survey data, we construct robust galaxy redshift distributions as a function of I_AB and R_AB magnitude and (R-I) color for galaxies between 0<z<2. We constrain models of the redshift evolution of galaxy clustering and find that significant growth of clustering has occurred from z>1 to z=0. A model in which the comoving scale-length, x_0, evolves linearly with redshift, x_0(z)=x_0(0)(1-Bz), fits the data better than the 'epsilon' model proposed by Groth and Peebles (1977). The clustering properties depend strongly on observed (R-I) color, with both the reddest and bluest galaxies exhibiting large clustering amplitudes and steeper slopes. Different observed (R-I) color ranges are sensitive to very disparate redshift regimes. Red galaxies with (R-I)~1.5 lie in a narrow redshift range centered at z~0.85 and have a comoving scale length of clustering at z=0.85 of x_0=5.0 +/-0.3 Mpc/h. These galaxies have early-type spectra and are likely progenitors of massive local ellipticals. The bluest galaxies with (R-I)~0 appear to be a mix of star-forming galaxies, both relatively local (z~0.3-0.6) dwarfs and bright z>1.4 galaxies, and broad-line AGN. We find that local blue dwarfs are relatively unclustered, with x_0=1.6 +/-0.2 Mpc/h. The z>1.4 blue galaxies have a larger clustering scale-length, x_0>5 Mpc/h.

Anisotropic gravitational radiation from a coalescing black hole binary can impart a recoil velocity of up to several hundred km/s to the remnant black hole. We examine the effects of recoiling massive black holes on their host stellar bulges, both for holes that escape their host and those that return to the galactic center via dynamical friction. We show that removal of a black hole via radiation recoil generally results in a rapidly-formed central core in the stellar system, with the effect being largest when the hole stays bound to the bulge and the recoil velocity is comparable to the bulge velocity dispersion. Black hole recoil therefore provides a mechanism for producing cores in some early type galaxies, but it is expected to be most efficient in faint ellipticals that are known to have steep density profiles. We argue that these results may hint at a significant role for gas in facilitating the coalescence of binary black holes in faint (power-law) early-type galaxies.

Much recent discussion about dark matter has been centered on two seemingly independent problems: the abundance of substructure in dark matter halos, and the cuspiness of the halos' inner density profile. We explore possible connections between the two problems by studying the gravitational scattering effects due to subhalos on the phase-space distribution of dark matter particles in the main halos. Our series of controlled numerical experiments indicates that the number and mass density of subhalos can be high enough to cause the collisionless dark matter particles in the inner part of a main halo to diffuse, flattening the main halo's inner cusp within a few dynamical times. Depending on the masses and concentration of the subhalos, the inner density profile of the whole system (main plus sub halos) can either steepen or flatten. Subhalo accretion can therefore introduce significant scatter in the inner density profiles of dark matter halos, offering a possible explanation for the range of profiles seen in both observations and cosmological simulations.

We present a kinetic theory for the evolution of the phase-space distribution of dark matter particles in galaxy halos in the presence of a cosmological spectrum of fluctuations. This theory introduces a new way to model the formation and evolution of halos, which traditionally have been investigated by analytic gravitational infall models or numerical N-body methods. Unlike the collisionless Boltzmann equation, our kinetic equation contains nonzero terms on the right-hand side arising from stochastic fluctuations in the gravitational potential due to substructures in the dark matter mass distribution. Using statistics for constrained Gaussian random fields in standard cosmological models, we show that our kinetic equation to second-order in perturbation theory is of the Fokker-Planck form, with one scattering term representing drift and the other representing diffusion in velocity-space. The drift is radial, and the drift and diffusion coefficients depend only on positions and not velocities; our relaxation process in the quasilinear regime is therefore different from the standard two-body relaxation. We provide explicit expressions relating these coefficients to the linear power spectrum of mass fluctuation and present results for the currently favored cold dark matter model with a nonzero cosmological constant. Solutions to this kinetic equation will provide a complete description of the cold dark matter spatial and velocity distributions for the average halo during the early phases of galaxy halo formation.

We use the two-dimensional Kolmogorov-Smirnov (KS) test to study the correlation between the 60 cosmic ray events above 4x10^19 eV from the AGASA experiment and the positions of infrared luminous galaxies from the IRAS PSCz catalog. These galaxies are expected to be hosts to gamma ray bursts (GRB) and magnetars, both of which are associated with core collapse supernovae and have been proposed as possible acceleration sites for ultra high energy cosmic rays. We find consistency between the models and the AGASA events to have been drawn from the same underlying distribution of positions on the sky with KS probabilities ~50%. Application of the same test to the 11 highest AGASA events above 10^20 eV, however, yields a KS probability of < 0.5%, rejecting the models at >99.5% significance level. Taken at face value, these highest energy results suggest that the existing cosmic ray events above 10^20 eV do not owe their origin to long burst GRBs, rapidly rotating magnetars, or any other events associated with core collapse supernovae. The larger data set expected from the AUGER experiment will test whether this conclusion is real or is a statistical fluke that we estimate to be at the 2 sigma level.

We present the results from a series of collisionless N-body simulations of major mergers of galaxy dark matter haloes with density profiles having either inner cusps or cores. Our simulations range from 2x10^5 to 10^7 particles, allowing us to probe the phase-space distribution of dark matter particles in the innermost regions (less than 0.005 virial radii) of cold dark matter haloes, a subject of much recent debate. We find that a major merger of two cored haloes yields a cored halo and does not result in a cuspy profile seen in many cosmological simulations. This result is unchanged if we consider mergers with parent mass ratios of 3:1 instead of 1:1. Mergers of a cuspy halo with either a cored halo or a second cuspy halo of equal mass, on the other hand, produce cuspy haloes with a slightly reduced inner logarithmic slope. Cuspy haloes, once formed, therefore appear resilient to major mergers. We find the velocity structure of the remnants to be mildly anisotropic, with a Maxwellian velocity distribution near the centre but not in the outer portions of the final haloes. Violent relaxation is effective only during the early phase of mergers, with phase mixing likely to be the dominant relaxation process at late times.

Strong gravitational lensing is a useful probe of both the intrinsic properties of the lenses and the cosmological parameters of the universe. The large number of model parameters and small sample of observed lens systems, however, have made it difficult to obtain useful constraints on more than a few parameters from lensing statistics. Here we examine how the recent WMAP measurements help improve the constraining power of statistics from the radio lensing survey JVAS/CLASS. We find that the absence of theta>3'' lenses in CLASS places an upper bound of beta<1.25 (1.60) at 68% (95%) CL on the inner density profile, rho \propto r^{-beta}, of cluster-sized halos. Furthermore, the favored power spectrum normalization is sigma_8 >= 0.7 (95% CL). We discuss two possibilities for stronger future constraints: a positive detection of at least one large-separation system, and next-generation radio surveys such as LOFAR.

We study contributions from inhomogeneous (patchy) reionization to arcminute scale ($1000 < \ell < 10,000$) cosmic microwave background (CMB) anisotropies. We show that inhomogeneities in the ionization fraction, rather than in the mean density, dominate both the temperature and the polarization power spectra. Depending on the ionization history and the clustering bias of the ionizing sources, we find that rms temperature fluctuations range from 2 $\mu$K to 8 $\mu$K and the corresponding values for polarization are over two orders of magnitude smaller. Reionization can significantly bias cosmological parameter estimates and degrade gravitational lensing potential reconstruction from temperature maps but not from polarization maps. We demonstrate that a simple modeling of the reionization temperature power spectrum may be sufficient to remove the parameter bias. The high-$\ell$ temperature power spectrum will contain some limited information about the sources of reionization.

This paper discusses the IAU Resolutions B1.3, B1.4, B1.5 and B1.9 (2000) that were adopted during the 24th General Assembly in Manchester, 2000 and provides details and explanations for these Resolutions. It is explained why they present significant progress over the corresponding IAU 1991 Resolutions and why they are necessary in the light of present accuracies in astrometry, celestial mechanics and metrology. In fact most of these Resolutions are consistent with astronomical models and software already in use. The metric tensors and gravitational potentials of both the Barycentric Celestial Reference System and Geocentric Celestial Reference System are defined and discussed. The necessity and relevance of the two celestial reference systems are explained. The transformations of coordinates and gravitational potentials are discussed. Potential coefficients parameterizing the post-Newtonian gravitational potentials are expounded. Simplified versions of the time transformations suitable for modern clock accuracies are elucidated. Various approximations used in the Resolutions are explicated and justified. Some models (e.g. for higher spin moments) that serve the purpose for estimating orders of magnitude have actually never been published before.

Sub-arcsecond lensing statistics depend sensitively on the inner mass profiles of low-mass objects and the faint-end slopes of the Schechter luminosity function and the Press-Schechter mass function. By requiring the luminosity and mass functions to give consistent predictions for the distribution of image separation below 1'', we show that dark matter halos with masses below 10^12 M_sun cannot have a single type of profile, be it the singular isothermal sphere (SIS) or the shallower ``universal'' dark matter profile. Instead, consistent results are achieved if we allow a fraction of the halos at a given mass to be luminous with the SIS profile, and the rest be dark with an inner logarithmic slope shallower than -1.5 to compensate for the steeper faint-end slope of the mass function compared with the luminosity function. We quantify how rapidly the SIS fraction must decrease with decreasing halo mass, thereby providing a statistical measure for the effectiveness of feedback processes on the baryon content in low-mass halos.

We develop a method based on the collisionless Boltzmann equation to calculate the gravitational clustering of relic neutrinos in realistic cosmological models dominated by cold dark matter (CDM) and the cosmological constant. This method can be used to estimate the phase-space distribution of any light particles in CDM halos. We find that neutrinos with masses > 0.3 eV cluster appreciably in dark matter halos above the galactic size. The resulting neutrino overdensity above the cosmic mean neutrino density increases with both the neutrino mass and the halo mass, ranging from about 10 for 0.3 eV neutrinos in 10^13 M_sun halos, to about 1500 for 1.8 eV neutrinos in 10^15 M_sun halos. We examine the implications of neutrino clustering for the Z-burst model of ultra high energy cosmic rays (UHECR), which interprets the observed events at E > 4 x 10^19 eV as decay products of Z-bosons from the resonant scattering between relic and high energy neutrinos and anti-neutrinos. We estimate the UHECR energy spectrum for various neutrino masses towards five of the most massive clusters in the local universe (within 100 Mpc): Virgo, Perseus-Pisces, Hydra, Centaurus, and Coma. The UHECR flux in the Z-burst model is expected to be significantly higher towards these clusters if m_nu > 0.3 eV and nearly isotropic otherwise.

Although it has previously been recognized that the CMB polarization measurement on a sky ring is able to separate the B and E-mode patterns, we show that the rotational symmetry of the CMB polarization measurement is indeed unique, in that it enables the B and E modes to be separately measured in a simple and clean manner. The separation of B and E modes can in principle be achieved even with a single-polarization detector of arbitrary beam pattern. Based on this premise, a specific observing strategy is suggested, where the telescope scans along a sky circle with the polarimeter axis always oriented to a fixed direction with respect to the radial direction of the sky circle. The observational strategy for B-mode measurements can further be refined by choosing appropriate beam sizes to suppress the accidental E-mode leakage. We present the expected variances of the measured E and B-mode surface-brightness fluctuations as functions of the sky annulus radius and beam size. The measured peak power is expected to be up to 10% of the total power contained in the CMB polarization fluctuations for both B and E modes. The temperature-polarization, or T-E, correlation can be measured by the present method as well. By adopting this strategy, the AMiBA experiment is expected to detect, in a single sky annulus, a 1\sigma$ B-mode signal in 340 hours.

We derive fully nonlinear expressions for temperature fluctuations from the kinetic Sunyaev-Zeldovich (SZ) effect, the scattering of cosmic microwave background photons off hot electrons in bulk motion. Our result reproduces the Ostriker-Vishniac effect to second order in perturbation theory but contains nonlinear corrections to the electron velocities and densities that were neglected previously. We use the recently developed halo model for nonlinear gravitational clustering to compute the nonlinear kinetic SZ power spectrum, which dominates the primary anisotropy on small angular scales.

A catalog containing milliarcsecond--accurate positions of 1332 extragalactic radio sources distributed over the northern sky is presented - the Very Long Baseline Array Calibrator Survey (VCS1). The positions have been derived from astrometric analysis of dual-frequency 2.3 and 8.4 GHz VLBA snapshot observations; in a majority of cases, images of the sources are also available. These radio sources are suitable for use in geodetic and astrometric experiments, and as phase-reference calibrators in high-sensitivity astronomical imaging. The VCS1 is the largest high-resolution radio survey ever undertaken, and triples the number of sources available to the radio astronomy community for VLBI applications. In addition to the astrometric role, this survey can be used in active galactic nuclei, Galactic, gravitational lens and cosmological studies. The VCS1 catalog will soon be available at www.nrao.edu/vlba/VCS1 .

Studies of the growth of cosmic perturbations are typically focused on galactic scales and above. In this paper we investigate the evolution of perturbations in baryons, photons, and dark matter for masses below 10^6 M_\odot (or wavenumbers above 100 Mpc^-1). Fluctuations on these scales are of interest and importance because they grow to become the earliest collapsed objects and provide the first light sources in the so called dark ages. We investigate both the linear evolution and the second-order nonlinear effects arising from the coupling of large-scale velocity fields to small-scale perturbations in the baryon density and the electron ionization fraction. We find that this second order nonlinear coupling dominates the growth of perturbations with masses \la 10^3 M_\odot immediately after recombination, enhancing the baryon fluctuation amplitudes by a factor of ~ 5, but the nonlinear effect does not persist at late times.

We investigate the statistics of gravitational lenses in flat, low-density cosmological models with different cosmic equations of state w. We compute the lensing probabilities as a function of image separation \theta using a lens population described by the mass function of Jenkins et al. and modeled as singular isothermal spheres on galactic scales and as Navarro, Frenk & White halos on cluster scales. It is found that COBE-normalized models with w > - 0.4 produce too few arcsecond-scale lenses in comparison with the JVAS/CLASS radio survey, a result that is consistent with other observational constraints on w. The wide-separation (\theta > 4'') lensing rate is a particularly sensitive probe of both w and the halo mass concentration. The absence of these systems in the current JVAS/CLASS data excludes highly concentrated halos in w < -0.7 models. The constraints can be improved by ongoing and future lensing surveys of > 10^5 sources.

Most gravitational lens systems consist of two or four observable images. The absence of detectable odd images allows us to place a lower limit on the power-law slope of the inner mass profile of lensing galaxies. Using a sample of six two-image radio-loud lens systems and assuming a singular power-law surface density (Sigma proportional to r^{-beta}) for the inner several kpc of the mass distribution, we find that there is less than a 10% probability that the data are consistent with profile slopes beta < 0.80. Furthermore, individual mass modeling yields beta > 0.85 for B0739+366 and beta > 0.91 for B1030+074. Modeling central black holes as additional point masses changes the constraints in these systems to beta > 0.84 and beta > 0.83, respectively. The inner mass profiles of lensing galaxies are therefore not much shallower than isothermal.

As part of a 4-year Cosmology and Particle Astrophysics (CosPA) Research Excellence Initiative in Taiwan, AMiBA $-$ a 19-element dual-channel 85-105 GHz interferometer array is being specifically built to search for high redshift clusters of galaxies via the Sunyaev-Zeldovich Effect (SZE). In addition, AMiBA will have full polarization capabilities, in order to probe the polarization properties of the Cosmic Microwave Background. AMiBA, to be sited on Mauna Kea in Hawaii or in Chile, will reach a sensitivity of $\sim 1$ mJy or 7$\mu$K in 1 hour. The project involves extensive international scientific and technical collaborations. The construction of AMiBA is scheduled to starting operating in early 2004.

An analytical understanding of the strongly nonlinear regime of gravitational collapse has been difficult to achieve. The only insight has been the stable clustering hypothesis, which assumes that the number of neighbors for objects averaged over small length scales is constant in time. Our recently proposed analytic halo model for N-point correlation functions now provides a tool for calculating gravitational clustering properties in the strongly nonlinear regime. This model also provides a new physical framework for an independent evaluation of the validity of the stable clustering hypothesis. We derive the asymptotic nonlinear behavior of the N-point correlation functions and pairwise peculiar velocities in terms of dark matter halo properties. We show that these statistics exhibit stable clustering only when the halo mass function and halo density profile obey specific relations. The long-cherished stable clustering hypothesis therefore is not necessarily realized in practice.

We present an analytic model for the fully nonlinear power spectrum P and bispectrum Q of the cosmological mass density field. The model is based on physical properties of dark matter halos, with the three main model inputs being analytic halo density profiles, halo mass functions, and halo-halo spatial correlations, each of which has been well studied in the literature. We demonstrate that this new model can reproduce the power spectrum and bispectrum computed from cosmological simulations of both an n=-2 scale-free model and a low-density cold dark matter model. To enhance the dynamic range of these large simulations, we use the synthetic halo replacement technique of Ma & Fry (2000a), where the original halos with numerically softened cores are replaced by synthetic halos of realistic density profiles. At high wavenumbers, our model predicts a slope for the nonlinear power spectrum different from the often-used fitting formulas in the literature based on the stable clustering assumption. Our model also predicts a three-point amplitude Q that is scale dependent, in contrast to the popular hierarchical clustering assumption. This model provides a rapid way to compute the mass power spectrum and bispectrum over all length scales where the input halo properties are valid. It also provides a physical interpretation of the clustering properties of matter in the universe.

High resolution N-body simulations of cold dark matter (CDM) models predict that galaxies and clusters have cuspy halos with excessive substructure. Observations reveal smooth halos with central density cores. One possible resolution of this conflict is that the dark matter is warm (WDM); this will suppress the power spectrum on small scales. The Lyman-alpha forest is a powerful probe of the linear power spectrum on these scales. We use collisionless N-body simulations to follow the evolution of structure in WDM models, and analyze artificial Lyman-alpha forest spectra extracted from them. By requiring that there is enough small-scale power in the linear power spectrum to reproduce the observed properties of the Lyman-alpha forest in quasar spectra, we derive a lower limit to the mass of the WDM particle of 750 eV. This limit is robust to reasonable uncertainties in our assumption about the temperature of the mean density gas (T0) at z=3. We argue that any model that suppresses the CDM linear theory power spectrum more severely than a 750 eV WDM particle cannot produce the Lyman-alpha forest.

We investigate the nonlinear two- and three-point correlation functions of the cosmological density field in Fourier space and test the popular hierarchical clustering model, that the three-point amplitude Q is independent of scale. In high-resolution N-body simulations of both n=-2 scale-free and cold dark matter models, we find that Q at late times is not constant but increases with wavenumber far into the nonlinear regime. Self-similar scaling also does not hold as rigorously for the three-point function as for the two-point function in the n=-2 simulation, perhaps a manifestation of the finite simulation volume. We suggest that a better understanding of the two- and three-point correlation functions in the nonlinear regime lies in the link to the density profiles of dark matter halos. We demonstrate and quantify how the slopes of the correlation functions are affected by the slope of the halo profiles using simple halo shapes and analytic clustering models.

We present new optical observations of the z=0.282 cluster Abell 697 from the Keck II telescope. Images show an unusual disturbed structure in the cD halo and a previously unknown faint gravitational lens arc. A spectrum of the arc did not yield a redshift, but its spectrum and colors suggest it lies at z>1.3. We construct models to reproduce the arc that show the potential is likely to be highly elliptical. We suggest that this cluster may have undergone a recent merger and is in the process of forming its cD galaxy. Analysis of X-ray data from ROSAT and ASCA suggests that the merging process is sufficiently advanced that the gas in the cluster has relaxed, and A697 lies near the L_x-T_x relation for normal clusters.

We present simple analytic approximations for the linear and fully evolved nonlinear mass power spectrum for spatially flat cold dark matter (CDM) cosmological models with quintessence (Q). Quintessence is a time evolving, spatially inhomogeneous energy component with negative pressure and an equation of state w_Q < 0. It clusters gravitationally on large length scales but remains smooth like the cosmological constant on small length scales. We show that the clustering scale is determined by the Compton wavelength of the Q-field and derive a shape parameter, \Gamma_Q, to characterize the linear mass power spectrum. The growth of linear perturbations as functions of redshift, w_Q, and matter density \Omega_m is also quantified. Calibrating to N-body simulations, we construct a simple extension of the formula by Ma (1998) that closely approximates the nonlinear power spectrum for a range of plausible QCDM models.

(abridged) We present a study of the nature and evolution of large-scale structure based on two independent redshift surveys of faint field galaxies conducted with the 176-fiber Norris Spectrograph on the Palomar 200-inch telescope. The two surveys together sparsely cover ~20 sq. degrees and contain 835 r < 21 mag galaxies with redshifts 0.2 < z < 0.5. Both surveys have a median redshift of z = 0.30. In order to obtain a rough estimate of the cosmic variance, we analyze the two surveys independently. We measure the comoving correlation length to be 3.70 +/- 0.13 h^-1 Mpc at z = 0.30 with a power-law slope gamma = 1.77 +/- 0.05. Dividing the sample into low (0.2 < z < 0.3) and high (0.32 < z < 0.5) redshift intervals, we see no evidence for a change in the comoving correlation length over the redshift range 0.2 < z < 0.5. Similar to the well-established results in the local universe, we find that intrinsically bright galaxies are more strongly clustered than intrinsically faint galaxies and that galaxies with little ongoing star formation, as judged from the rest-frame equivalent width of the [OII]3727, are more strongly clustered than galaxies with significant ongoing star formation. The rest-frame pairwise velocity dispersion of the sample is 326^+67_-52 km s^-1, ~25% lower than typical values measured locally. The appearance of the galaxy distribution, particularly in the more densely sampled Abell 104 field, is quite striking. The pattern of sheets and voids which has been observed locally continues at least to z ~ 0.5. A friends-of-friends analysis of the galaxy distribution supports the visual impression that > 90% of all galaxies at z < 0.5 are part of larger structures with overdensities of > 5.

In these lectures I highlight some key features of massive neutrinos in the context of cosmology. I first review the thermal history and the free-streaming kinematics of the uniform cosmic background neutrinos. I then describe how fluctuations in the phase space distributions of neutrinos and other particles arise and evolve after neutrino decoupling according to the linear perturbation theory of gravitational instability. The different clustering properties of massive neutrinos (aka hot dark matter) and cold dark matter are contrasted. The last part discusses the nonlinear stage of gravitational clustering and highlights the effects of massive neutrinos on the formation of cosmological structure.

A non-relativistic Fokker-Planck equation for the electron distribution function is formulated incorporating the effects of stochastic acceleration by whistler-mode waves and Coulomb collisions. The stationary solution $f$ to the equation, subject to a zero-flux boundary condition, is found to be a generalized Lorentzian (or kappa) distribution, which satisfies $f\propto v^{-2(\kappa+1)}$ for large velocity $v$, where $\kappa$ is the spectral index. The parameter $\kappa$ depends strongly on the relative wave intensity $R$. Taking into account the critical energy required for resonance of electrons with whistlers, we calculate a range of values of $R$ for each of a number of different space plasmas for which kappa distributions can be expected to be formed. This study is one of the first in the literature to provide a theoretical justification for the formation of generalized Lorentzian (or kappa) particle distribution functions in space plasmas.

I present a general discussion of the evolution and model-dependence of both the linear and nonlinear power spectrum of density fluctuations. The features of the linear power spectrum in cosmological models with cold dark matter (CDM) and cold+hot dark matter (C+HDM) are reviewed, and useful analytical approximations are summarized. Cosmological numerical simulation results are then used to illustrate the process of nonlinear gravitational clustering and to compute the nonlinear power spectrum. An analytical approximation relating the linear and nonlinear power spectrum in C+HDM models is presented.

This paper presents an analytical formula that closely approximates the fully nonlinear power spectrum of matter fluctuations for redshift $z\approx 5$ to 0 over a wide range of cosmologically interesting flat models with varying matter density $\Omega_m$ and neutrino fraction $\Omega_\nu$. The functional form is motivated by analytical solutions in asymptotic regimes, but in order to obtain accurate approximations, the coefficients are calculated from fits to the nonlinear power spectrum computed from numerical simulations of four cosmological models. The transformation from the linear to nonlinear power spectrum depends on $\Omega_m, \Omega_\nu$, and time. A simple scaling rule is introduced, which greatly simplifies the construction of the functional form and allows the formula to depend directly on the rms linear mass fluctuation $\sigma_8$ instead of on an effective spectral index as in earlier work.

This paper presents a detailed theoretical study of the two-point correlation function $\xi$ for both dark matter halos and the matter density field in five cosmological models with varying matter density $\Omega_m$ and neutrino fraction $\Omega_\nu$. The objectives of this systematic study are to evaluate the nonlinear gravitational effects on $\xi$, to contrast the behavior of $\xi$ for halos vs. matter, and to quantify the redshift evolution of $\xi$ and its dependence on cosmological parameters. Overall, $\xi$ for halos exhibits markedly slower evolution than $\xi$ for matter, and its redshift dependence is much more intricate than the single power-law parameterization used in the literature. Of particular interest is that the redshift evolution of the halo-halo correlation length $r_0$ depends strongly on $\Omega_m$ and $\Omega_\nu$, being slower in models with lower $\Omega_m$ or higher $\Omega_\nu$. Measurements of $\xi$ to higher redshifts can therefore be a potential discriminator of cosmological parameters. The evolution rate of $r_0$ for halos within a given model increases with time, passing the phase of fixed comoving clustering at $z\sim 1$ to 3 toward the regime of stable clustering at $z\sim 0$. The shape of the halo-halo $\xi$, on the other hand, is well approximated by a power law with slope -1.8 in all models and is not a sensitive model discriminator.

We present a summary report on a detailed study of the mass and dynamical state of the Corona Borealis supercluster using the largest redshift survey of this object to date. Extensive $N$-body experiments are performed on simulated superclusters in both critically-flat and low-density cosmological models to test the robustness and errors of the mass estimators used on the data. From the redshifts and dynamics of 528 galaxies in the supercluster, the mass of Corona Borealis is estimated to be at least $3\times 10^{16} h^{-1} M_\odot$, and the mass-to-light ratio is $560 h (M/L)_\odot$, yielding a matter density parameter of $\Omega_m\approx 0.4$. We also discuss the prospects for mapping the mass distribution on supercluster scales with gravitational lensing.

We add the $\alpha-$ effect in the dynamo driven accretion disk model proposed by Tout & Pringle (1992), i.e., a dynamo model depends on the physical processes such as Parker instability, Balbus-Hawley instability, magnetic field reconnection and $\alpha-\omega$ mean field dynamo as well. The $\alpha-$ effect in the dynamo mechanism is determined by the strength of turbulence of the accretion flow. When the turbulent Mach number $M_t$ is less than 0.25, the solutions of the magnetic fields oscillate around their equilibrium values. The increase of the value of $M_t$ makes the amplitude of the oscillation smaller and the period longer, but does not affect the equilibrium values. The Shakura-Sunyaev viscosity parameter $\alpha_{SS}$ oscillates around the equilibrium value of 0.33. When the turbulent Mach number $M_t$ is larger than 0.25, the magnetic field components reach a stable state. In the non-linear dynamo region, the critical turbulent Mach number $M_t$ is 0.44 rather than 0.25. The oscillating magnetic fields and viscosity parameter can explain the basic properties of the dwarf nova eruptions and some properties of quiescent disks (Armitage et al. 1996).

We use the cyclotron-maser theory to explain the core emission from the magnetosphere of pulsars. As a kind of direct and efficient maser type of emission, it can give rise to escaping radiation with extremely high brightness temperature and narrow angle with respect to the magnetic axis. We find that the growth rates and real frequencies of the O-mode electromagnetic wave propagating parallel to the magnetic fields depend on the ratio of the plasma frequency $\omega_p$ and the gyrofrequency $\omega_b$ rather than the plasma frequency alone, as described by other models. The emission takes place in the region where the magnitude of $\omega_p/\omega_b$ is $10^{-2}$. The corresponding altitude is about a few decades of neutron star radius, where the magnetic field strength is about $10^6-10^8 G$. The qualitative spectrum and the lower frequency cut-off of the radio emission is obtained by this model.

(Abridged) We present a study of the structure and dynamics of the Corona Borealis Supercluster (z ~ 0.07) based on the redshifts for 528 galaxies in the supercluster. The galaxy distribution within Corona Borealis is clumpy and appears overall to be far from relaxed. A background supercluster at z ~ 0.11 makes a substantial contribution to the projected surface density of galaxies in the Corona Borealis field. In order to estimate the mass of the supercluster, we have assumed that the mass of the supercluster is proportional to v^2r, where v and r are suitable scale velocity and radius, respectively, and we have used N-body simulations of both critical- and low-density universes to determine the applicability of standard mass estimators based on this assumption. Although superclusters are obviously not in equilibrium, our simulations demonstrate that the virial mass estimator yields mass estimates with an insignificant bias and a dispersion of only ~25% for objects with overdensities >~ 5. Non-uniform spatial sampling can, however, cause systematic underestimates by as much as 30%. All of our simulated superclusters turn out to be bound, and based on the overdensity of the Corona Borealis supercluster, we believe it is also very likely to be bound and may well have started to collapse. The mass of Corona Borealis is at least 3 x 10^16 solar masses, which yields a B(AB)-band mass-to-light ratio of 564h (M/L)_solar on scales of ~20 h^-1 Mpc. The background supercluster has a mass-to-light ratio of 726h (M/L)_solar. By comparing the supercluster mass-to-light ratios with the critical mass-to-light ratio required to close the universe, we determine that Omega_0 >~ 0.4 on supercluster scales.

Any viable cosmological model must produce enough structure at early epochs to explain the amount of gas associated with high-redshift damped Ly$\alpha$ systems. We study the evolution of damped Ly$\alpha$ systems at redshifts $z\ge 2$ in cold dark matter (CDM) and cold+hot dark matter (CDM+HDM) models using both N-body and hydrodynamic simulations. Our approach incorporates the effects of gas dynamics, and we find that all earlier estimates which assumed that all the baryons in dark matter halos would contribute to damped Ly$\alpha$ absorption have overestimated the column density distribution $f(N)$ and the fraction of neutral dense gas $\Omega_g$ in damped Ly$\alpha$ systems. The differences are driven by ionization of hydrogen in the outskirts of galactic halos and by gaseous dissipation near the halo centers, and they tend to exacerbate the problem of late galaxy formation in CDM+HDM models. We only include systems up to the highest observed column density $N\sim 10^{21.8}$ cm$^{-2}$ in the estimation of $\Omega_g$ for a fair comparison with data. If the observed $f(N)$ and $\Omega_g$ inferred from a small number of confirmed and candidate absorbers are robust, the amount of gas in damped Ly$\alpha$ systems at high redshifts in the $\Omega_\nu=0.2$ CDM+HDM model falls well below the observations.

This paper presents simple analytic approximations to the linear power spectra, linear growth rates, and rms mass fluctuations for both components in a family of cold+hot dark matter (CDM+HDM) models that are of current cosmological interest. The formulas are valid for a wide range of wavenumber, neutrino fraction, redshift, and Hubble constant: $k\lo 10\,h$ Mpc$^{-1}$, $0.05\lo \onu\lo 0.3$, $0\le z\lo 15$, and $0.5\lo h \lo 0.8$. A new, redshift-dependent shape parameter $\Gamma_\nu=a^{1/2}\onu h^2$ is introduced to simplify the multi-dimensional parameter space and to characterize the effect of massive neutrinos on the power spectrum. The physical origin of $\Gamma_\nu$ lies in the neutrino free-streaming process, and the analytic approximations can be simplified to depend only on this variable and $\onu$. Linear calculations with these power spectra as input are performed to compare the predictions of $\onu\lo 0.3$ models with observational constraints from the reconstructed linear power spectrum and cluster abundance. The usual assumption of an exact scale-invariant primordial power spectrum is relaxed to allow a spectral index of $0.8\lo n\le 1$. It is found that a slight tilt of $n=0.9$ (no tensor mode) or $n=0.95$ (with tensor mode) in $\onu\sim 0.1-0.2$ CDM+HDM models gives a power spectrum similar to that of an open CDM model with a shape parameter $\Gamma=0.25$, providing good agreement with the power spectrum reconstructed by Peacock and Dodds (1994) and the observed cluster abundance.

This paper presents a systematic treatment of the linear theory of scalar gravitational perturbations in the synchronous gauge and the conformal Newtonian (or longitudinal) gauge. It differs from others in the literature in that we give, in both gauges, a complete discussion of all particle species that are relevant to any flat cold dark matter (CDM), hot dark matter (HDM), or CDM+HDM models (including a possible cosmological constant). The particles considered include CDM, baryons, photons, massless neutrinos, and massive neutrinos (an HDM candidate), where the CDM and baryons are treated as fluids while a detailed phase-space description is given to the photons and neutrinos. Particular care is applied to the massive neutrino component, which has been either ignored or approximated crudely in previous works. Isentropic initial conditions on super-horizon scales are derived. The coupled, linearized Boltzmann, Einstein and fluid equations that govern the evolution of the metric and density perturbations are then solved numerically in both gauges for the standard CDM model and two CDM+HDM models with neutrino mass densities $\onu=0.2$ and 0.3, assuming a scale-invariant, adiabatic spectrum of primordial fluctuations. We also give the full details of the cosmic microwave background anisotropy, and present the first accurate calculations of the angular power spectra in the two CDM+HDM models including photon polarization, higher neutrino multipole moments, and helium recombination. The numerical programs for both gauges are available at http://arcturus.mit.edu/cosmics/ .

We report the most recent results from high-resolution numerical simulations of structure formation in two flat cold+hot dark matter models with neutrino mass densities $\onu=0.2$ and 0.3. We find that structure forms too late in all CDM+HDM models with $\onu>0.2$ to account for the amount of dense neutral gas in high-redshift damped Lyman-$\alpha$ systems. The $\onu=0.2$ model at $z\approx0$ provides a better match to observations than the pure CDM model.

The abundance of galactic systems at high redshifts can impose a strong constraint on the cold+hot dark matter (CDM+HDM) models. The hot component reduces the excessive small-scale power in the COBE-normalized CDM model but also delays the epoch of galaxy formation. We present results from the first numerical simulations that have enough dynamic range to address accurately the issue of high-redshift halo abundances in CDM+HDM models. Equivalent high-resolution particle-particle/particle-mesh $N$-body simulations are performed for spatially flat models with $\Omega_\nu =0.3$ and 0.2 (with $H_0=50$ km s$^{-1}$ Mpc$^{-1}$ and $\Omega_b=0.05$). We study the constraints placed on the models by the high-redshift quasar space density and by the mass fraction in neutral dense gas associated with damped Ly$\alpha$ systems. We find that even with optimistic assumptions, the much-studied $\Omega_\nu=0.3$ model does not produce enough massive halos to account for the observed abundance of quasars at $z>4$. The model passes this test if $\Omega_\nu$ is decreased to 0.2. Both models do not produce enough high column-density halos to account for the amount of gas in damped Ly$\alpha$ systems at $z\go 3$: the $\Omega_\nu=0.3$ model falls short by a factor $\sim$80; the $\Omega_\nu=0.2$ model by a factor $\sim$3. We conclude that only CDM+HDM models with $\Omega_\nu\lo 0.2$ can match observations at high redshift, implying an upper bound of 4.7 eV on the most massive light neutrino (presumably the $\tau$).

We present a systematic treatment of the linear theory of scalar gravitational perturbations in the synchronous gauge and the conformal Newtonian (or longitudinal) gauge. We first derive the transformation law relating the two gauges. We then write down in parallel in both gauges the coupled, linearized Boltzmann, Einstein and fluid equations that govern the evolution of the metric perturbations and the density fluctuations of the particle species. The particle species considered include cold dark matter (CDM), baryons, photons, massless neutrinos, and massive neutrinos (a hot dark matter or HDM candidate), where the CDM and baryon components are treated as fluids while a detailed phase-space description is given to the photons and neutrinos. The linear evolution equations presented are applicable to any $\Omega=1$ model with CDM or a mixture of CDM and HDM. Isentropic initial conditions on super-horizon scales are derived. The equations are solved numerically in both gauges for a CDM+HDM model with $\Omega_{\rm cold}=0.65,$ $\Omega_{\rm hot}=0.3$, and $\Omega_{\rm baryon}=0.05$. We discuss the evolution of the metric and the density perturbations and compare their different behaviors outside the horizon in the two gauges. In a companion paper we integrate the geodesic equations for the neutrino particles in the perturbed conformal Newtonian background metric computed here. The purpose is to obtain an accurate sampling of the neutrino phase space for the HDM initial conditions in $N$-body simulations of the CDM+HDM models.

This paper presents a general-relativistic N-body technique for evolving the phase space distribution of massive neutrinos in linear perturbation theory. The method provides a much more accurate sampling of the neutrino phase space for the HDM initial conditions of N-body simulations in a cold+hot dark matter universe than previous work. Instead of directly sampling the phase space at the end of the linear era, we first compute the evolution of the metric perturbations by numerically integrating the coupled, linearized Einstein, Boltzmann, and fluid equations for all particle species. We then sample the phase space shortly after neutrino decoupling at redshift z=10^9 when the distribution is Fermi-Dirac. To follow the trajectory of each neutrino, we subsequently integrate the geodesic equations for each neutrino in the perturbed background spacetime from z=10^9 to z=13.55, using the linearized metric found in the previous calculation to eliminate discreteness noise. The positions and momenta resulting from this integration represent a fair sample of the full neutrino phase space and can be used as HDM initial conditions for N-body simulations of nonlinear structure evolution in this model. A total of 21 million neutrino particles are used in a 100 Mpc box, with Omega_cdm=0.65, Omega_hdm=0.30, Omega_baryon=0.05, and Hubble constant H_0=50. We find that correlations develop in the neutrino densities and momenta which are absent when only the zeroth-order Fermi-Dirac distribution is considered.

SO(10) cosmic strings formed during the phase transition Spin(10) $\rightarrow$ SU(5) $\times{\cal Z}_2$ are studied. Two types of strings --- one effectively Abelian and one non-Abelian --- are constructed and the string solutions are calculated numerically. The non-Abelian string can catalyze baryon number violation via the ``twisting'' of the scalar field which causes mixing of leptons and quarks in the fermion multiplet. The non-Abelian string is also found to have the lower energy possibly for the entire range of the parameters in the theory. Scattering of fermions in the fields of the strings is analyzed, and the baryon number violation cross section is calculated. The role of the self-adjoint parameters is discussed and the values are computed.