We find, from our study of binary spiral galaxies in the Sloan Digital Sky Survey Data Release 6, that the relative orientation of disks in binary spiral galaxies is consistent with their being drawn from a random distribution of orientations. For 747 isolated pairs of luminous disk galaxies, the distribution of phi, the angle between the major axes of the galaxy images, is consistent with a uniform distribution on the interval [0 degrees, 90 degrees]. With the assumption that the disk galaxies are oblate spheroids, we can compute cos(beta), where beta is the angle between the rotation axes of the disks. In the case that one galaxy in the binary is face-on or edge-on, the tilt ambiguity is resolved, and cos(beta) can be computed unambiguously. For 94 isolated pairs with at least one face-on member, and for 171 isolated pairs with at least one edge-on member, the distribution of cos(beta) is statistically consistent with the distribution of cos(i) for isolated disk galaxies. This result is consistent with random orientations of the disks within pairs.

We identify, categorize, and quantify alignment effects among host and satellite galaxies using a spectroscopically-confirmed, low-redshift (z<0.23) galaxy sample from the Sloan Digital Sky Survey Data Release 6. Consistent with other recent findings, we find that satellite galaxies (SGs) of red, centrally concentrated (elliptical) host galaxies (HGs) with radial velocity separation |Delta_V|<600 km/s preferentially reside near the projected major axes of their HGs. Among these, this preference is strongest for highly concentrated, red SGs. We find that fractional anisotropy increases with decreasing \Delta_V and Delta_R and is nearly 40% greater among the closest SGs (Delta_R<250 kpc/h) relative to more distant (Delta_R>500 kpc/h) SGs. For highly concentrated SGs at small (<300 kpc/h) projected separation, we observe a strong radial (hostward) alignment signal in isophotal position angles (PAs) due to isophotal twisting and contamination that is not present when using galaxy model PAs. Among objects for which both isophotal and galaxy model PAs agree to <15 degrees, this elongation signal is significantly weaker. We also investigate the "Holmberg Effect," a well-known result wherein nearby (<40 kpc/h) SGs of large, inclined spiral HGs were seen to preferentially reside near the minor axes of their HGs. Survey limitations preclude a strict test of this effect using only SDSS spectroscopic galaxies. By adopting a looser set of cuts than those of Holmberg's study, we recover a comparable preference among faint blue SGs for the HG minor axis at marginal significance (~3 sigma). We conclude that several types of alignment likely exist among different galaxy populations, but that the observed nature and strength of alignment trends depend sensitively on both selection criteria and on the method used to determine galaxy orientation.

We analyze the absolute magnitude (M_r) and color (u-r) of low redshift (z<0.06) galaxies in the Sloan Digital Sky Survey Data Release 6. galaxies with nearly exponential profiles (Sloan parameter fracDeV < 0.1) fall on the blue sequence of the color - magnitude diagram; if, in addition, these exponential galaxies have M_r < -19, they show a dependence of u-r color on apparent axis ratio q expected for a dusty disk galaxy. By fitting luminosity functions for exponential galaxies with different values of q, we find that the dimming is well described by the relation Delta M_r = 1.27 (log q)^2, rather than the Delta M = C log q law that is frequently assumed. When the absolute magnitudes of bright exponential galaxies are corrected to their "face-on" value, M_r^f = M_r - Delta M_r, the average u-r color is linearly dependent on M_r^f for a given value of q. Nearly face-on exponential galaxies (q > 0.9) have a shallow dependence of mean u-r color on M_r^f (0.096 magnitudes redder for every magnitude brighter). By comparison, nearly edge-on exponential galaxies (q < 0.3) are 0.265 magnitudes redder for every magnitude brighter. When the dimming law Delta M = 1.27 (log q)^2 is used to create an inclination-corrected sample of bright exponential galaxies, their apparent shapes are confirmed to be consistent with a distribution of mildly non-circular disks, with median short-to-long axis ratio gamma = 0.22 and median disk ellipticity epsilon = 0.08.

The apparent shapes of spiral galaxies in the 2-Micron All Sky Survey Large Galaxy Atlas are used to constrain the intrinsic shape of their disks. When the distribution of apparent axis ratios is estimated using a nonparametric kernel method, the shape distribution is inconsistent with axisymmetry at the 90% confidence level in the B band and at the 99% confidence level in the K band. If spirals are subdivided by Hubble type, the late-type spirals (Sc and later) are consistent with axisymmetry, while the earlier spirals are strongly inconsistent with axisymmetry. The distribution of disk ellipticity can be fitted adequately with either a Gaussian or a lognormal distribution. The best fits for the late spirals imply a median ellipticity of epsilon = 0.07 in the B band and epsilon = 0.02 in the K band. For the earlier spirals, the best fits imply a median ellipticity of epsilon = 0.18 in the B band and epsilon = 0.30 in the K band. The observed scatter in the Tully-Fisher relation, for both late and early spirals, is consistent with the disk ellipticity measured in the B band. This indicates that excluding spirals of Hubble type earlier than Sc will reduce the intrinsic scatter in the Tully-Fisher relation used as a distance indicator.

Using a sample of galaxies from the Sloan Digital Sky Survey (SDSS) Data Release 4, we study the trends relating surface brightness profile type and apparent axis ratio to the local galaxy environment. We use the SDSS parameter `fracDeV' to quantify the profile type. We find that galaxies with M_r > -18 are mostly described by exponential profiles in all environments. Galaxies with -21 < M_r < -18 mainly have exponential profiles in low density environments and de Vaucouleurs profiles in high density environments. The most luminous galaxies, with M_r < -21, are mostly described by de Vaucouleurs profiles in all environments. For galaxies with M_r < -19, the fraction of de Vaucouleurs galaxies is a monotonically increasing function of local density, while the fraction of exponential galaxies is monotonically decreasing. For a fixed surface brightness profile type, apparent axis ratio is frequently correlated with environment. As the local density of galaxies increases, we find that for -20 < M_r < -18, galaxies of all profile types become slightly rounder, on average; for -22 < M_r < -20, galaxies with exponential profiles tend to become flatter, while galaxies with de Vaucouleurs profiles become rounder; for M_r < -22, galaxies with exponential profiles become flatter, while the de Vaucouleurs galaxies become rounder in their inner regions, yet exhibit no change in their outer regions. We comment on how the observed trends relate to the merger history of galaxies.

For a sample of 96,951 galaxies from the Sloan Digital Sky Survey Data Release 3, we study the distribution of apparent axis ratios as a function of r-band absolute magnitude and surface brightness profile type. We use the parameter fracDeV to quantify the profile type (fracDeV = 1 for a de Vaucouleurs profile; fracDeV = 0 for an exponential profile). When the apparent axis ratio q_{am} is estimated from the moments of the light distribution, the roundest galaxies are very bright (M_r \sim -23) de Vaucouleurs galaxies and the flattest are modestly bright (M_r \sim -18) exponential galaxies. When the apparent axis ratio q_{25} is estimated from the axis ratio of the 25 mag/arcsec^2 isophote, we find that de Vaucouleurs galaxies are flatter than exponential galaxies of the same absolute magnitude. For a given surface brightness profile type, very bright galaxies are rounder, on average, than fainter galaxies. We deconvolve the distributions of apparent axis ratios to find the distribution of the intrinsic short-to-long axis ratio gamma, assuming constant triaxiality T. For all profile types and luminosities, the distribution of apparent axis ratios is inconsistent with a population of oblate spheroids, but is usually consistent with a population of prolate spheroids. Bright galaxies with a de Vaucouleurs profile (M_r < -21.84, fracDeV > 0.9) have a distribution of q_{am} that is consistent with triaxiality in the range 0.4 < T < 0.8, with mean intrinsic axis ratio 0.66 < gamma < 0.69. The fainter de Vaucouleurs galaxies are best fit with prolate spheroids (T = 1) with mean axis ratio gamma = 0.51.

The disks of spiral galaxies are generally elliptical rather than circular. The distribution of ellipticities can be fit with a log-normal distribution. For a sample of 12,764 galaxies from the Sloan Digital Sky Survey Data Release 1 (SDSS DR1), the distribution of apparent axis ratios in the i band is best fit by a log-normal distribution of intrinsic ellipticities with ln epsilon = -1.85 +/- 0.89. For a sample of nearly face-on spiral galaxies, analyzed by Andersen and Bershady using both photometric and spectroscopic data, the best fitting distribution of ellipticities has ln epsilon = -2.29 +/- 1.04. Given the small size of the Andersen-Bershady sample, the two distribution are not necessarily inconsistent. If the ellipticity of the potential were equal to that of the light distribution of the SDSS DR1 galaxies, it would produce 1.0 magnitudes of scatter in the Tully-Fisher relation, greater than is observed. The Andersen-Bershady results, however, are consistent with a scatter as small as 0.25 magnitudes in the Tully-Fisher relation.

We investigate the effect of rotating, triaxial halos on disk galaxies through an extensive set of numerical N-body simulations. Our simulations use a rigid potential field for the halos and bulges and collisionless particles for the disks. The triaxiality and the rotation rate of the halo are varied, as well as the masses of all three galaxy components. We analyze both the bar stability and spiral response of the disk under these conditions. We characterize most of our models by the mass ratio of the disk to the halo at 2.3 disk scale lengths, (M_d/M_h)_Rsun. For models with a mass ratio greater than 0.8, a halo pattern speed Omega_h = 6.7 km/s/kpc, and an intermediate-to-major axis ratio q_b=0.85, a strong bar will develop within 3 Gyr, even for models with a bulge mass M_b=0.3M_d. Models in which the bulge mass is reduced by half develop bars earlier and with lower Omega_h. We create an artificial Hubble sequence of galaxies by varying the bulge-to-disk ratio of our models from 0 to 2.5. The torque induced by a rotating, non-axisymmetric halo creates bisymmetric spiral structure in the disk. We find the the pitch angle of the spiral arms in these models follows the same general trend found in observations of spiral galaxies, namely that later type galaxies have higher pitch angles. Our simulations follow closely the observational relation of spiral pitch angle with maximum rotational velocity of the disk, where galaxies with faster rotation have more tightly wound spiral arms. This relation is followed in our simulations regardless of whether the dominating mass component of the galaxy is the disk, the halo, or the bulge.

We estimate the distribution of apparent axis ratios $q$ for galaxies in the Sloan Digital Sky Survey (SDSS) Early Data Release. We divide the galaxies by profile type (de Vaucouleurs versus exponential) as well as by color ($u^* - r^* \leq 2.22$ versus $u^* - r^* > 2.22$). The axis ratios found by fitting models to the surface photometry are generally smaller than those found by taking the second moments of the surface brightness distribution. Using the axis ratios found from fitting models, we find that galaxies with de Vaucouleurs profiles have axis ratio distributions which are inconsistent, at the 99% confidence level, with their being a population of randomly oriented oblate spheroids. Red de Vaucouleurs galaxies are slightly rounder, on average, than blue de Vaucouleurs galaxies. By contrast, blue galaxies with exponential profiles appear very much flatter, on average, than red galaxies with exponential profiles. The red exponential galaxies are primarily disk galaxies seen nearly edge-on, with reddening due to the presence of dust, rather than to an intrinsically red stellar population.

(Abridged) Motivated by forthcoming data from the Sloan Digital Sky Survey, we present a theoretical framework that can be used to interpret Principal Component Analysis (PCA) of disk galaxy properties. We use the formalism introduced by Mo, Mao, & White to compute the observable properties of galaxies in a number of model populations, varying assumptions about which physical parameters determine structural quantities and star formation histories. We then apply PCA to these model populations. Our baseline model assumes that halo mass, spin parameter, and formation redshift are the governing input parameters and that star formation is determined by surface density through a Schmidt law. In all cases, the first principal component is primarily a measure of the shape of the spectral energy distribution (SED), and it is usually driven by variations in the spin parameter, which influences star formation through the disk surface density. The second and (in some cases) third principal components consist mainly of ``scale'' parameters like luminosity, disk radius, and circular velocity. However, the detailed division of these scale parameters, the disk surface brightness, and the rotation curve slope among the principal components changes significantly from model to model. Our calculations yield predictions of principal component structure for the baseline model of disk galaxy formation, and a physical interpretation of these predictions. They also show that PCA can test the core assumptions of that model and reveal the presence of additional physical parameters that may govern observable galaxy properties.

Elliptical galaxies formed in a major merger tend to become more nearly spherical with time, thanks to the gravitational effect of their central black hole (or black holes). Observational results indicate that elliptical galaxies with older stellar populations (t > 7.5 Gyr) have rounder central isophotes than ellipticals with younger stellar populations. In addition, older ellipticals tend to have core profiles, while younger ellipticals tend to have power-law profiles. Numerical simulations of galaxy mergers indicate that if one or both of the progenitors have a central black hole with a mass around 0.2 percent of the mass in stars, then the effect of the black hole(s) is to make the central regions of the remnant rounder, with a characteristic time of a few gigayears.

We test the hypothesis that the apparent axis ratio of an elliptical galaxy is correlated with the age of its stellar population. We find that old ellipticals (with estimated stellar ages t > 7.5 Gyr) are rounder on average than younger ellipticals. The statistical significance of this shape difference is greatest at small radii; a Kolmogorov-Smirnov test comparing the axis ratios of the two populations at R = Re/16 yields a statistical significance greater than 99.96 percent. The relation between age and apparent axis ratio is linked to the core/powerlaw surface brightness profile dichotomy. Core ellipticals have older stellar populations, on average, than powerlaw ellipticals and are rounder in their inner regions. Our results are consistent with a scenario in which powerlaw ellipticals are formed in gas-rich mergers while core ellipticals form in dissipationless mergers, with cores formed and maintained by the influence of a binary black hole.

The sizes and shapes of voids in a galaxy survey depend not only on the physics of structure formation, but also on the sampling density of the survey and on the algorithm used to define voids. Using an N-body simulation with a CDM power spectrum, we study the properties of voids in samples with different number densities of galaxies, both in redshift space and in real space. When voids are defined as regions totally empty of galaxies, their characteristic volume is strongly dependent on sampling density; when they are defined as regions whose density is 0.2 times the mean galaxy density, the dependence is less strong. We compare two void-finding algorithms, one in which voids are nonoverlapping spheres, and one, based on the algorithm of Aikio and Mahonen, which does not predefine the shape of a void. Regardless of the algorithm chosen, the characteristic void size is larger in redshift space than in real space, and is larger for low sampling densities than for high sampling densities. We define an elongation statistic Q which measures the tendency of voids to be stretched or squashed along the line of sight. Using this statistic, we find that at sufficiently high sampling densities (comparable to the number densities of galaxies brighter than L_*), large voids tend to be slightly elongated along the line of sight in redshift space.

For a sample of 22 core ellipticals and 22 power law ellipticals, I examine the apparent axis ratios of isophotes as a function of the dynamical time at the isophote's mean radius. The two types of elliptical galaxy have a significantly different distribution of shapes in the region where the dynamical time is short (about 1 Myr); the core galaxies are rounder in this region. In the outer regions, the two types have indistinguishable shape distributions. I perform a more detailed analysis of the apparent and intrinsic shapes at two chosen locations: a inner region, where the dynamical time is 0.75 Myr, and an outer region, where the dynamical time is 50 Myr. The inner isophotes of core galaxies are consistent with their being randomly oriented oblate objects, with a peak intrinsic axis ratio of c/a = 0.86. The outer regions of core galaxies are inconsistent with their being randomly oriented oblate or prolate objects. The oblate hypothesis cannot be rejected for the inner and outer regions of power law ellipticals at the 90% confidence level. The observed shapes are consistent with the hypothesis that the inner regions of core galaxies are nearly spherical and nearly axisymmetric, perhaps as the result of the physical process which formed the core.

We present the results of Smoothed Particle Hydrodynamics (SPH) simulations of the formation of a massive counterrotating disk in a spiral galaxy. The current study revisits and extends (with SPH) previous work carried out with sticky particle gas dynamics, in which adiabatic gas infall and a retrograde gas-rich dwarf merger were tested as the two most likely processes for producing such a counterrotating disk. We report on experiments with a cold primary similar to our Galaxy, as well as a hot, compact primary modeled after NGC 4138. We have also conducted numerical experiments with varying amounts of prograde gas in the primary disk, and an alternative infall model (a spherical shell with retrograde angular momentum). The structure of the resulting counterrotating disks is dramatically different with SPH. The disks we produce are considerably thinner than the primary disks and those produced with sticky particles. The time-scales for counterrotating disk formation are shorter with SPH because the gas loses kinetic energy and angular momentum more rapidly. Spiral structure is evident in most of the disks, but an exponential radial profile is not a natural byproduct of these processes. The infalling gas shells that we tested produce counterrotating bulges and rings rather than disks. The presence of a considerable amount of preexisting prograde gas in the primary causes, at least in the absence of star formation, a rapid inflow of gas to the center and a subsequent hole in the counterrotating disk. In general, our SPH experiments yield stronger evidence to suggest that the accretion of massive counterrotating disks drives the evolution of the host galaxies towards earlier (S0/Sa) Hubble types.

The non-Euclidean geometry of spacetime induces an anisotropy in the apparent correlation function of high-redshift quasars. This effect can constrain the cosmological constant \Lambda independent of any assumptions about evolution of luminosities, sizes, or clustering. We examine the prospects for distinguishing between low-density (\Omega_0=0.1-0.4) cosmological models with flat and open space geometry using the quasar samples anticipated from the 2dF and Sloan redshift surveys. We show that even these large quasar surveys are likely to reside in the "sparse sampling" regime, so that measurement errors in the correlation function obey Poisson statistics. As a result: (a) one can devise a simple maximum-likelihood scheme for estimating clustering parameters, (b) one can generate Monte Carlo realizations of correlation function measurements without creating artificial quasar distributions, and (c) for fixed quasar number, a deeper survey over a smaller area has greater statistical power than a shallow, large-area survey. Adopting recent (quite uncertain) estimates of the quasar correlation length, we find that the 2dF and Sloan samples can provide clear discrimination between flat and open geometries for \Omega <= 0.2 but only marginal discrimination for \Omega = 0.4. Clear discrimination is possible for \Omega = 0.4 if the true quasar correlation length is a factor of two larger, and a high-density survey of 30,000 quasars in 200 square degrees would provide clear discrimination even for the lower correlation length.

The Sa(r) galaxy NGC 4138 has been recently found to contain an extensive counterrotating disk which appears to be still forming. Up to a third of the stars in the disk system may be on retrograde orbits. A counterrotating ring of H II regions, along with extended counterrotating H I gas, suggests that the retrograde material has been recently acquired in the gas phase and is still trickling in. Using numerical simulations, we have attempted to model the process by which the counterrotating mass has been accreted by this galaxy. We investigate two possibilities: continuous retrograde infall of gas, and a retrograde merger with a gas-rich dwarf galaxy. Both processes are successful in producing a counterrotating disk of the observed mass and dimensions without heating up the primary significantly. Contrary to our experience with a fiducial cold, thin primary disk, the gas-rich merger works well for the massive, compact primary disk of NGC 4138 even though the mass of the dwarf galaxy is a significant fraction of the mass of the primary disk. Although we have restricted ourselves mainly to coplanar infall and mergers, we report on one inclined infall simulation as well. We also explore the possibility that the H-alpha ring seen in the inner half of the disk is a consequence of counterrotating gas clouds colliding with corotating gas already present in the disk and forming stars in the process.

The shapes of isolated Bok globules and embedded dense cores of molecular clouds are analyzed using a nonparametric method, under the alternate hypotheses that they are randomly oriented prolate objects or that they are randomly oriented oblate objects. In all cases, the prolate hypothesis gives a better fit to the data. If Bok globules are oblate, they must be very flat; the average axis ratio is b/a = 0.3, and few or no globules can have b/a > 0.7. If Bok globules are prolate, then the mean axis ratio is b/a = 0.5. For most data samples of dense cores, the randomly-oriented oblate hypothesis can be rejected at the 99% confidence level. If the dense cores are prolate, their mean axis ratio is approximately 0.4 to 0.5. Dense cores are significantly different in shape from the clouds in which they are embedded; clouds have flatter apparent shapes, and are inconsistent with a population of randomly oriented axisymmetric objects.

Using two-dimensional numerical simulations of gravitational clustering, with initial power spectra of a power-law form with index n, we compare the properties of voids in real space to their properties in redshift space. Both the void probability function (VPF) and the underdense probability function (UPF) are enhanced in redshift space. The enhancement is greatest in the limit n=-2. The VPF and UPF treat voids statistically and assume that two-dimensional voids are circular; we present an algorithm which detects individual voids. Voids found by this algorithm are ellipses whose enclosed density of galaxies falls below a threshold density. When voids are identified using this algorithm, the mean void size and the maximum void size both increase in going from real space to redshift space. The increase is greatest in the limit n=-2. In redshift space, the principal axes of the largest voids in the n=2 and n=0 simulations show a statistically significant tendency (at the 95% confidence level) to be distributed anisotropically, relative to the line of sight from the origin to the void center.

I compute the estimated distribution function f(q) for the apparent axis ratio q of various types of stellar systems, using a nonparametric kernel method. I then invert f(q) to find the distribution of intrinsic axis ratios, using two different hypotheses: first, that the stellar systems are all oblate, and second, that they are all prolate. The shapes of globular clusters in our own galaxy are consistent, at the 99% confidence level, with both the oblate and prolate hypothesis. The shapes of dwarf galaxies in the Virgo cluster are consistent, at the 99% confidence level, with the prolate hypothesis, but inconsistent with the oblate hypothesis. The shapes of star clusters in the Large Magellanic Cloud, of ordinary elliptical galaxies, of brightest cluster ellipticals, and of galaxy clusters are all inconsistent, at the 99% confidence level, with both the oblate and prolate hypotheses. The globular clusters in our galaxy are older than their half-mass relaxation time, and are most likely rotationally flattened oblate spheroids. The other stellar systems considered are generally younger than their half-mass relaxation time, and thus are triaxial bodies flattened by anisotropy of their velocity dispersion.

We present results of numerical simulations of the formation of a massive counterrotating gas disk in a spiral galaxy. Using a hierarchical tree gravity solver combined with a sticky-particle gas dissipation scheme for our simulations, we have investigated three mechanisms: episodic and continuous gas infall, and a merger with a gas-rich dwarf galaxy. We find that both episodic and continuous gas infall work reasonably well and are able to produce a substantial gas counterrotating disk without upsetting the stability of the existing disk drastically, but it is very important for the gas to be well-dispersed in phase-space and not form concentrated clumps prior to its absorption by the disk galaxy. The initial angular momentum of the gas also plays a crucial role in determining the scale length of the counterrotating disk formed and the time it takes to form. The rate of infall, i.e. the mass of gas falling in per unit time, has to be small enough to preclude excessive heating of the preexisting disk. It is much easier in general to produce a smaller counterrotating disk than it is to produce an extensive disk whose scale length is similar to that of the original prograde disk. A gas-rich dwarf merger does not appear to be a viable mechanism to produce a massive counterrotating disk, because only a very small dwarf galaxy can produce a counterrotating disk without increasing the thickness of the existing disk by an order of magnitude, and the time-scale for this process is prohibitively long because it makes it very unlikely that several such mergers can accumulate a massive counterrotating disk over a Hubble time.

Because the transformation from distance to redshift is nonlinear, maps in redshift space become increasingly distorted as the redshift z becomes greater. As noted by Alcock \& Paczy\'nski (1979), observed redshift distortions can be used to estimate the deceleration parameter $q_0$. If $q_0$ is greater than -1, voids in redshift space will be elongated along the line of sight. In addition, distant voids will have a greater volume in redshift space than nearby voids. Accurate measurement of the volume and the axis ratio of voids, as a function of their central redshift, will provide an estimate of $q_0$. To test this method of estimating $q_0$, I create a two-dimensional toy universe, free of peculiar velocities, in which the galaxies are located near the walls of Voronoi cells. The galaxies are then mapped into redshift space, adopting different values of $q_0$. In redshift space, I estimate the area and the axis ratio of the voids by fitting ellipses within the voids, using an algorithm which maximizes the area of the empty ellipses and ensures that ellipses do not overlap. The accuracy of the estimated values of $q_0$ is limited by the intrinsic scatter in the size and shape of the voids. In the toy universe, distinguishing between a $q_0 = -1$ universe and a $q_0 = 1/2$ universe requires a survey which goes to a depth z > 0.1 in redshift space. Peculiar velocities will create an additional source of uncertainty for the values of $q_0$ measured in N-body simulations and in the real universe.