We report detections of three z ~ 2.5 submillimeter-selected galaxies (SMGs; SMM J14011+0252, SMM J14009+0252, SMM J04431+0210) in the lowest rotational transition of the carbon monoxide molecule (CO J = 1-0) and one nondetection (SMM J04433+0210). For the three galaxies we detected, we find a line-integrated brightness temperature ratio of the J = 3-2 and 1-0 lines of 0.68 +/- 0.08; the 1-0 line is stronger than predicted by the frequent assumption of equal brightnesses in the two lines and by most single-component models. The observed ratio suggests that mass estimates for SMGs based on J = 3-2 observations and J = 1-0 column density or mass conversion factors are low by a factor of 1.5. Comparison of the 1-0 line intensities with intensities of higher-J transitions indicates that single-component models for the interstellar media in SMGs are incomplete. The small dispersion in the ratio, along with published detections of CO lines with J_upper > 3 in most of the sources, indicates that the emission is from multi-component interstellar media with physical structures common to many classes of galaxies. This result tends to rule out the lowest scaling factors between CO luminosity and molecular gas mass, and further increases molecular mass estimates calibrated against observations of galaxies in the local universe. We also describe and demonstrate a statistically sound method for finding weak lines in broadband spectra that will find application in searches for molecular lines from sources at unknown redshifts.

We consider the effect of rapid rotation on the light curves of neutron stars with hot polar caps. For $P \approx 3$ms spin periods, the pulse fractions can be as much as an order of magnitude larger than with simple slowly-rotating (Schwarzschild) estimates. Doppler boosting, in particular, leads to characteristic distortion and ``soft lags'' in the pulse profiles, which are easily measurable in light curves with moderate energy resolution. With $\sim 10^5$ photons it should also be possible to isolate the more subtle distortions of light travel time variations and frame dragging. Detailed analysis of high quality millisecond pulsar data from upcoming X-ray missions must include these effects.

We examine the non-linear stability of the Wisdom-Holman (WH) mapping applied to the integration of perturbed, highly eccentric two-body orbits. We find that the method is unstable and introduces artificial chaos into the computed trajectories, unless the step size is chosen small enough to resolve pericenter. The origin of the instability is analyzed analytically using the Stark problem as a fiducial test case. We similarly examine the robustness of several alternative methods: a regularized WH map due to Mikkola (1997); the potential-splitting (PS) approach of Lee et al. (1997); and two types of Stark-based methods. Comparative simulations show the regularized WH map to be stable at high eccentricities, and an enhanced PS algorithm to perform well when close encounters are additionally present. We find Stark-based schemes to be of marginal use in N-body type integrations.

It has recently been shown that the rate of angular momentum relaxation in nearly-Keplerian star clusters is greatly increased by a process termed resonant relaxation (Rauch & Tremaine 1996), who also argued that tidal disruption of stars in galactic nuclei containing massive black holes could be noticeably enhanced by this process. We describe here the results of numerical simulations of resonant tidal disruption which quantitatively test the predictions made by Rauch & Tremaine. The simulation method is based on an N-body routine incorporating cloning of stars near the loss cone and a semi-relativistic symplectic integration scheme. We also briefly describe the discovery of chaos in the Wisdom-Holman symplectic integrator applied to highly eccentric orbits and propose a modified integration scheme that remains robust under these conditions. We find that resonant disruption rates exceed their non-resonant counterparts by an amount consistent with the predictions; in particular, we estimate the net tidal disruption rate for a fully resonant cluster to be about twice that of its non-resonant counterpart. No significant enhancement in rates is observed outside the critical radius. Relativistic quenching of the effect is found to occur for hole masses M>8*10^7 solar masses. The numerical results combined with the observed properties of galactic nuclei indicate that for most galaxies the resonant enhancement to tidal disruption rates will be very small.

We demonstrate the existence of an enhanced rate of angular momentum relaxation in nearly Keplerian star clusters, such as those found in the centers of galactic nuclei containing massive black holes. The enhanced relaxation arises because the radial and azimuthal orbital frequencies in a Keplerian potential are equal, and hence may be termed {\em resonant relaxation}. We explore the dynamics of resonant relaxation using both numerical simulations and order-of-magnitude analytic calculations. We find that the resonant angular momentum relaxation time is shorter than the non-resonant relaxation time by of order $M_\star/M$, where $M_\star$ is the mass in stars and $M$ is the mass of the central object. Resonance does not enhance the energy relaxation rate. We examine the effect of resonant relaxation on the rate of tidal disruption of stars by the central mass; we find that the flux of stars into the loss cone is enhanced when the loss cone is empty, but that the disruption rate averaged over the entire cluster is not strongly affected. We show that relativistic precession can disable resonant relaxation near the main-sequence loss cone for black hole masses comparable to those in galactic nuclei. Resonant dynamical friction leads to growth or decay of the eccentricity of the orbit of a massive body, depending on whether the distribution function of the stars is predominantly radial or tangential. The accelerated relaxation implies that there are regions in nuclear star clusters that are relaxed in angular momentum but not in energy; unfortunately, these regions are not well-resolved in nearby galaxies by the Hubble Space Telescope.