In an effort to better understand collisionless relaxation processes in gravitational systems, we investigate one-dimensional models. Taking advantage of a Hermite-Legendre expansion of relevant distribution functions, we present analytical and numerical behaviors of Maxwell-Boltzmann entropy. In particular, we modestly perturb systems about a separable-solution equilibrium and observe their collisionless evolution to a steady state. We verify the time-independence of fine-grained entropy in these systems before turning our attention to the behavior of coarse-grained entropy. We also verify that there is no analogue to the collisional H-theorem for these systems. Competing terms in the second-order coarse-grained entropy make it impossible to guarantee continuously increasing entropy. However, over dynamical time-scales the coarse-grained entropy generally increases, with small oscillations occurring. The lack of substantive differences between the entropies in test-particle and self-gravitating cases suggests that phase mixing, rather than violent relaxation associated with potential changes, more significantly drives the coarse-grained entropy evolution. The effects of violent relaxation can be better quantified through analysis of energy distributions rather than phase-space distributions.

Building on the development of a Hermite-Legendre analysis of one-dimensional gravitating collisionless systems, we present a technique for determining the steady states of such systems. This provides an important component for understanding the physics involved in the relaxation of these kinds of systems. As the dark matter structures in the universe should have traits in common with these systems, insight into this relaxation can provide clues to larger astrophysical questions. For large perturbation strengths, we determine physically motivated parameter ranges for the simplest families of steady states as well as their stability. We also demonstrate that any set of initial conditions in the linear regime can be resolved into unique time-independent and time-dependent modes. Combinations of time-independent modes then describe the steady state of any system linearly perturbed from equilibrium. These results highlight the importance of initial conditions over relaxation mechanisms in the evolution of these systems.

Numerous previous studies have investigated the phenomenon wherein initially spherical N-body systems are distorted to triaxial shapes. We report on an investigation of a previously described orbital instability that should oppose triaxiality. After verifying the instability with numerical orbit integrations that extend the original analysis, we search for evidence of the instability in N-body systems that become triaxial. Our results highlight the difficulty in separating dynamical process from finite-N effects. While we argue that our analysis points to the presence of the instability in simulated triaxial systems, discreteness appears to play a role in mimicking the instability. This suggests that predicting the shapes of real-world systems, such as dark matter halos around galaxies, based on such simulations involves more uncertainty than previously thought.

Using an $N$-body evolution code that does not rely on softened potentials, I have created a suite of unbound interacting cluster pair simulations. The motions of the centers of mass of the clusters have been tracked and compared to the trajectories of point masses interacting via one of four different softened potential prescriptions. I find that the relationship between the impact parameter of the cluster interaction and the point-mass softening length that best approximates each cluster's center-of-mass motion depends on the adopted prescription. In general, the range of allowed softening lengths grows roughly linearly with the impact parameter, but zero softening is acceptable in the majority of situations. In an $N$-body simulation that adopts a fixed softening length, such relationships lead to the possibility of two-body effects, like dynamical friction, being either larger or smaller than the corresponding cluster situation. Further consideration of more specific $N$-body situations leads estimating that a very small fraction of point-mass encounters experience two-body effects significantly different than those of equivalent clusters.

The equilibria formed by the self-gravitating, collisionless collapse of simple initial conditions have been investigated for decades. We present the results of our attempts to describe the equilibria formed in $N$-body simulations using thermodynamically-motivated models. Previous work has suggested that it is possible to define distribution functions for such systems that describe maximum entropy states. These distribution functions are used to create radial density and velocity distributions for comparison to those from simulations. A wide variety of $N$-body code conditions are used to reduce the chance that results are biased by numerical issues. We find that a subset of initial conditions studied lead to equilibria that can be accurately described by these models, and that direct calculation of the entropy shows maximum values being achieved.

The current paradigm for understanding galaxy formation in the universe depends on the existence of self-gravitating collisionless dark matter. Modeling such dark matter systems has been a major focus of astrophysicists, with much of that effort directed at computational techniques. Not surprisingly, a comprehensive understanding of the evolution of these self-gravitating systems still eludes us, since it involves the collective nonlinear dynamics of many-particle systems interacting via long-range forces described by the Vlasov equation. As a step towards developing a clearer picture of collisionless self-gravitating relaxation, we analyze the linearized dynamics of isolated one-dimensional systems near thermal equilibrium by expanding their phase space distribution functions f(x,v) in terms of Hermite functions in the velocity variable, and Legendre functions involving the position variable. This approach produces a picture of phase-space evolution in terms of expansion coefficients, rather than spatial and velocity variables. We obtain equations of motion for the expansion coefficients for both test-particle distributions and self-gravitating linear perturbations of thermal equilibrium. N-body simulations of perturbed equilibria are performed and found to be in excellent agreement with the expansion coefficient approach over a time duration that depends on the size of the expansion series used.

We studied the application of statistical reconstruction algorithms, namely maximum likelihood and least squares methods, to the problem of event reconstruction in a dual phase liquid xenon detector. An iterative method was developed for in-situ reconstruction of the PMT light response functions from calibration data taken with an uncollimated gamma-ray source. Using the techniques described, the performance of the ZEPLIN-III dark matter detector was studied for 122 keV gamma-rays. For the inner part of the detector (R<100 mm), spatial resolutions of 13 mm and 1.6 mm FWHM were measured in the horizontal plane for primary and secondary scintillation, respectively. An energy resolution of 8.1% FWHM was achieved at that energy. The possibility of using this technique for improving performance and reducing cost of scintillation cameras for medical applications is currently under study.

Since globular clusters (GCs) are old, low-N systems their dynamics is widely believed to be fully dominated by collisional two-body processes, and their surface brightness profiles are fit by King models. However, for many GCs, especially those with HST-resolved central regions, and `extra-tidal' features, King models provide poor fits. We suggest that this is partly because collisionless dynamics is also important and contribute to shaping the cluster properties. We show using time-scale and length-scale arguments that except for the very centers of clusters, collisionless dynamics should be more important than collisional. We then fit 38 GCs analyzed by Noyola and Gebhardt (2006) with (collisional) King and (collisionless) DARKexp models over the full available radial range, and find that the latter provide a better fit to 29 GCs; for six of these the fit is at least ~5x better in term of rms. DARKexp models are theoretically derived maximum entropy equilibrium states of self-gravitating collisionless systems and have already been shown to fit the results of dark matter N-body simulations. (We do not attempt fits with ad hoc fitting functions.)

We present an analysis of two thermodynamic techniques for determining equilibria of self-gravitating systems. One is the Lynden-Bell entropy maximization analysis that introduced violent relaxation. Since we do not use the Stirling approximation which is invalid at small occupation numbers, our systems have finite mass, unlike Lynden-Bell's isothermal spheres. (Instead of Stirling, we utilize a very accurate smooth approximation for $\ln{x!}$.) The second analysis extends entropy production extremization to self-gravitating systems, also without the use of the Stirling approximation. In addition to the Lynden-Bell (LB) statistical family characterized by the exclusion principle in phase-space, and designed to treat collisionless systems, we also apply the two approaches to the Maxwell-Boltzmann (MB) families, which have no exclusion principle and hence represent collisional systems. We implicitly assume that all of the phase-space is equally accessible. We derive entropy production expressions for both families, and give the extremum conditions for entropy production. Surprisingly, our analysis indicates that extremizing entropy production rate results in systems that have maximum entropy, in both LB and MB statistics. In other words, both thermodynamic approaches lead to the same equilibrium structures.

We report experimental upper limits on WIMP-nucleon elastic scattering cross sections from the second science run of ZEPLIN-III at the Boulby Underground Laboratory. A raw fiducial exposure of 1,344 kg.days was accrued over 319 days of continuous operation between June 2010 and May 2011. A total of eight events was observed in the signal acceptance region in the nuclear recoil energy range 7-29 keV, which is compatible with background expectations. This allows the exclusion of the scalar cross-section above 4.8E-8 pb near 50 GeV/c^2 WIMP mass with 90% confidence. Combined with data from the first run, this result improves to 3.9E-8 pb. The corresponding WIMP-neutron spin-dependent cross-section limit is 8.0E-3 pb. The ZEPLIN programme reaches thus its conclusion at Boulby, having deployed and exploited successfully three liquid xenon experiments of increasing reach.

Plastic scintillators are widely used in industry, medicine and scientific research, including nuclear and particle physics. Although one of their most common applications is in neutron detection, experimental data on their response to low-energy nuclear recoils are scarce. Here, the relative scintillation efficiency for neutron-induced nuclear recoils in a polystyrene-based plastic scintillator (UPS-923A) is presented, exploring recoil energies between 125 keV and 850 keV. Monte Carlo simulations, incorporating light collection efficiency and energy resolution effects, are used to generate neutron scattering spectra which are matched to observed distributions of scintillation signals to parameterise the energy-dependent quenching factor. At energies above 300 keV the dependence is reasonably described using the semi-empirical formulation of Birks and a kB factor of (0.014+/-0.002) g/MeVcm^2 has been determined. Below that energy the measured quenching factor falls more steeply than predicted by the Birks formalism.

ZE3RA is the software package responsible for processing the raw data from the ZEPLIN-III dark matter experiment and its reduction into a set of parameters used in all subsequent analyses. The detector is a liquid xenon time projection chamber with scintillation and electroluminescence signals read out by an array of 31 photomultipliers. The dual range 62-channel data stream is optimised for the detection of scintillation pulses down to a single photoelectron and of ionisation signals as small as those produced by single electrons. We discuss in particular several strategies related to data filtering, pulse finding and pulse clustering which are tuned to recover the best electron/nuclear recoil discrimination near the detection threshold, where most dark matter elastic scattering signatures are expected. The software was designed assuming only minimal knowledge of the physics underlying the detection principle, allowing an unbiased analysis of the experimental results and easy extension to other detectors with similar requirements.

The ZEPLIN-III experiment is operating in its second phase at the Boulby Underground Laboratory in search of dark matter WIMPs. The major upgrades to the instrument over its first science run include lower background photomultiplier tubes and installation of a plastic scintillator veto system. Performance results from the veto detector using calibration and science data in its first six months of operation in coincidence with ZEPLIN-III are presented. With fully automated operation and calibration, the veto system has maintained high stability and achieves near unity live time relative to ZEPLIN-III. Calibrations with a neutron source demonstrate a rejection of 60% of neutron-induced nuclear recoils in ZEPLIN-III that might otherwise be misidentified as WIMPs. This tagging efficiency reduces the expected untagged nuclear recoil background from neutrons during science data taking to a very low rate of ~0.2 events per year in the WIMP acceptance region. Additionally, the veto detector provides rejection of 28% of gamma-ray induced background events, allowing the sampling of the dominant source of background in ZEPLIN-III - multiple scatter gamma-rays with rare topologies. Since WIMPs will not be tagged by the veto detector, and tags due to gamma-rays and neutrons are separable, this population of multiple scatter events may be characterised without biasing the analysis of candidate WIMP signals in the data.

Certain thermal non-equilibrium situations, outside of the astrophysical realm, suggest that entropy production extrema, instead of entropy extrema, are related to stationary states. In an effort to better understand the evolution of collisionless self-gravitating systems, we investigate the role of entropy production and develop expressions for the entropy production rate in two particular statistical families that describe self-gravitating systems. From these entropy production descriptions, we derive the requirements for extremizing the entropy production rate in terms of specific forms for the relaxation function in the Boltzmann equation. We discuss some implications of these relaxation functions and point to future work that will further test this novel thermodynamic viewpoint of collisionless relaxation.

The design, optimisation and construction of an anti-coincidence veto detector to complement the ZEPLIN-III direct dark matter search instrument is described. One tonne of plastic scintillator is arranged into 52 bars individually read out by photomultipliers and coupled to a gadolinium-loaded passive polypropylene shield. Particular attention has been paid to radiological content. The overall aim has been to achieve a veto detector of low threshold and high efficiency without the creation of additional background in ZEPLIN-III, all at a reasonable cost. Extensive experimental measurements of the components have been made, including radioactivity levels and performance characteristics. These have been used to inform a complete end-to-end Monte Carlo simulation that has then been used to calculate the expected performance of the new instrument, both operating alone and as an anti-coincidence detector for ZEPLIN-III. The veto device will be capable of rejecting over 65% of coincident nuclear recoil events from neutron background in the energy range of interest in ZEPLIN-III. This will reduce the background in ZEPLIN-III from ~0.4 to ~0.14 events per year in the WIMP acceptance region, a significant factor in the event of a non-zero observation. Furthermore, in addition to providing valuable diagnostic capabilities, the veto is capable of tagging over 15% for gamma-ray rejection, all whilst contributing no significant additional background. In conjunction with the replacement of the internal ZEPLIN-III photomultiplier array, the new veto is expected to improve significantly the sensitivity of the ZEPLIN-III instrument to dark matter, allowing spin independent WIMP-nucleon cross sections below 1E-8 pb to be probed.

The ZEPLIN-III experiment in the Palmer Underground Laboratory at Boulby uses a 12kg two-phase xenon time projection chamber to search for the weakly interacting massive particles (WIMPs) that may account for the dark matter of our Galaxy. The detector measures both scintillation and ionisation produced by radiation interacting in the liquid to differentiate between the nuclear recoils expected from WIMPs and the electron recoil background signals down to ~10keV nuclear recoil energy. An analysis of 847kg.days of data acquired between February 27th 2008 and May 20th 2008 has excluded a WIMP-nucleon elastic scattering spin-independent cross-section above 8.1x10(-8)pb at 55GeV/c2 with a 90% confidence limit. It has also demonstrated that the two-phase xenon technique is capable of better discrimination between electron and nuclear recoils at low-energy than previously achieved by other xenon-based experiments.

Using a suite of self-gravitating, collisionless N-body models, we systematically explore a parameter space relevant to the onset and behavior of the radial orbit instability (ROI), whose strength is measured by the systemic axis ratios of the models. We show that a combination of two initial conditions, namely the velocity anisotropy and the virial ratio, determines whether a system will undergo ROI and exactly how triaxial the system will become. A third initial condition, the radial shape of the density profile, plays a smaller, but noticeable role. Regarding the dynamical development of the ROI, the instability a) begins after systems collapse to their most compact configuration and b) evolves fastest when a majority of the particles have radially anisotropic orbits while there is a lack of centrally-concentrated isotropic orbits. We argue that this is further evidence that self-reinforcing torques are the key to the onset of the ROI. Our findings support the idea that a separate orbit instability plays a role in halting the ROI.

For a decade, N-body simulations have revealed a nearly universal dark matter density profile, which appears to be robust to changes in the overall density of the universe and the underlying power spectrum. Despite its universality, the physical origin of this profile has not yet been well understood. Semi--analytic models by Barnes et al. (2005) have suggested that the density structure of dark matter halos is determined by the onset of the radial orbit instability (ROI). We have tested this hypothesis using N-body simulations of collapsing dark matter halos with a variety of initial conditions. For dynamically cold initial conditions, the resulting halo structures are triaxial in shape, due to the mild aspect of the instability. We examine how variations in initial velocity dispersion affect the onset of the instability, and find that an isotropic velocity dispersion can suppress the ROI entirely, while a purely radial dispersion does not. The quantity sigma^2/vc^2 is a criterion for instability, where regions with sigma^2/vc^2 <~1 become triaxial due to the ROI or other perturbations. We also find that the radial orbit instability sets a scale length at which the velocity dispersion changes rapidly from isotropic to radially anisotropic. This scale length is proportional to the radius at which the density profile changes shape, as is the case in the semi--analytic models; however, the coefficient of proportionality is different by a factor of ~2.5. We conclude that the radial orbit instability is likely to be a key physical mechanism responsible for the nearly universal profiles of simulated dark matter halos.

We present a non-parametric method for decomposition of the light of disk galaxies into disk, bulge and bar components. We have developed and tested the method on a sample of 68 disk galaxies for which we have acquired I-band photometry. The separation of disk and bar light relies on the single assumption that the bar is a straight feature with a different ellipticity and position angle from that of the projected disk. We here present the basic method, but recognise that it can be significantly refined. We identify bars in only 47% of the more nearly face-on galaxies in our sample. The fraction of light in the bar has a broad range from 1.3% to 40% of the total galaxy light. If low-luminosity galaxies have more dominant halos, and if halos contribute to bar stability, the luminosity functions of barred and unbarred galaxies should differ markedly; while our sample is small, we find only a slight difference of low significance.

Using two-dimensional velocity maps and I-band photometry, we have created mass models of 40 spiral galaxies using the Milgrom relation (the basis of modified Newtonian dynamics, or MOND) to complement previous work. A Bayesian technique is employed to compare several different dark matter halo models to Milgrom and Newtonian models. Pseudo-isothermal dark matter halos provide the best statistical fits to the data in a majority of cases, while the Milgrom relation generally provides good fits as well. We also find that Milgrom models give mass-to-light ratios that roughly correlate with galaxy color, as predicted by stellar population models. A subsample of galaxies in the Hydra cluster follow a tight relation between mass-to-light and color, but one that is significantly different from relations found in previous studies. Ruling out the Milgrom relation with rotational kinematics is difficult due to systematic uncertainties in the observations as well as underlying model assumptions. We discuss in detail two galaxies for which the Milgrom relation appears to fail and find that relaxing the assumption of constant stellar mass-to-light ratio can maintain Milgrom models' viability.

There is no accepted mechanism that explains the equilibrium structures that form in collisionless cosmological N-body simulations. Recent work has identified nonextensive thermodynamics as an innovative approach to the problem. The distribution function that results from adopting this framework has the same form as for polytropes, but the polytropic index is now related to the degree of nonextensiveness. In particular, the nonextensive approach can mimic the equilibrium structure of dark matter density profiles found in simulations. We extend the investigation of this approach to the velocity structures expected from nonextensive thermodynamics. We find that the nonextensive and simulated N-body rms-velocity distributions do not match one another. The nonextensive rms-velocity profile is either monotonically decreasing or displays little radial variation, each of which disagrees with the rms-velocity distributions seen in simulations. We conclude that the currently discussed nonextensive models require further modifications in order to corroborate dark matter halo simulations. (adapted from TeX)

It has long been realized that dark matter halos formed in cosmological N-body simulations are characterized by density profiles rho(r) that, when suitably scaled, have similar shapes. Additionally, combining the density and velocity dispersion profiles sigma(r), each of which have decidedly nonpower-law shapes, leads to quantity rho/(sigma^3) that is a power-law in radius over 3 orders of magnitude in radius. Halos' velocity anisotropy profiles beta(r) vary from isotropic near the centers of halos to quite radially anisotropic near the virial radius. Finally, there appears to be a nearly linear correlation between beta and the logarithmic density slope gamma for a wide variety of halos. This work is part of a continuing investigation of the above interrelationships and their origins using analytical and semi-analytical techniques. Our finding suggest that the nearly linear beta--gamma relationship is not just another expression of scale-free rho/(sigma^3) behavior. We also note that simultaneously reproducing density and anisotropy profiles like those found in simulations requires beta(r) and gamma(r) to have similar shapes, leading to nearly linear beta--gamma correlations. This work suggests that the beta--gamma and power-law rho/(sigma^3) relations have distinct physical origins.

We investigate the connection between collisionless equilibria and the phase-space relation between density $\rho$ and velocity dispersion $\sigma$ found in simulations of dark matter halo formation, $F=\psd \propto r^{-\alpha}$. Understanding this relation will shed light on the physics relevant to collisionless collapse and on the subsequent structures formed. We show that empirical density profiles that provide good fits to N-body halos also happen to have nearly scale-free \psd distributions when in equilibrium. We have also done a preliminary investigation of variables other than $r$ that may match or supercede the correlation with $F$. In the same vein, we show that $\rho/\sigma^m$, where $m=3$ is the most appropriate combination to use in discussions of the power-law relationship. Since the mechanical equilibrium condition that characterizes the final systems does not by itself lead to power-law $F$ distributions, our findings prompt us to posit that dynamical collapse processes (such as violent relaxation) are responsible for the radial power-law nature of the \psd distributions of virialized systems.

Although N-body studies of dark matter halos show that the density profiles, rho(r), are not simple power-laws, the quantity rho/sigma^3, where sigma(r) is the velocity dispersion, is in fact a featureless power-law over ~3 decades in radius. In the first part of the paper we demonstrate, using the semi-analytic Extended Secondary Infall Model (ESIM), that the nearly scale-free nature of rho/sigma^3 is a robust feature of virialized halos in equilibrium. By examining the processes in common between numerical N-body and semi-analytic approaches, we argue that the scale-free nature of rho/sigma^3 cannot be the result of hierarchical merging, rather it must be an outcome of violent relaxation. The empirical results of the first part of the paper motivate the analytical work of the second part of the paper, where we use rho/sigma^3 proportional to r^{-alpha} as an additional constraint in the isotropic Jeans equation of hydrostatic equilibrium. Our analysis shows that the constrained Jeans equation has different types of solutions, and in particular, it admits a unique ``periodic'' solution with alpha=1.9444. We derive the analytic expression for this density profile, which asymptotes to inner and outer profiles of rho ~ r^{-0.78}, and rho ~ r^{-3.44}, respectively.

We investigate a hypothesis regarding the origin of the scalelength in halos formed in cosmological N-body simulations. This hypothesis can be viewed as an extension of an earlier idea put forth by Merritt and Aguilar. Our findings suggest that a phenomenon related to the radial orbit instability is present in such halos and is responsible for density profile shapes. This instability sets a scalelength at which the velocity dispersion distribution changes rapidly from isotropic to radially anisotropic. This scalelength is reflected in the density distribution as the radius at which the density profile changes slope. We have tested the idea that radially dependent velocity dispersion anisotropy leads to a break in density profile shape by manipulating the input of a semi-analytic model to imitate the velocity structure imposed by the radial orbit instability. Without such manipulation, halos formed are approximated by single power-law density profiles and isotropic velocity distributions. Halos formed with altered inputs display density distributions featuring scalelengths and anisotropy profiles similar to those seen in cosmological N-body simulations.

High resolution N-body simulations have all but converged on a common empirical form for the shape of the density profiles of halos, but the full understanding of the underlying physics of halo formation has eluded them so far. We investigate the formation and structure of dark matter halos using analytical and semi-analytical techniques. Our halos are formed via an extended secondary infall model (ESIM); they contain secondary perturbations and hence random tangential and radial motions which affect the halo's evolution at it undergoes shell-crossing and virialization. Even though the density profiles of NFW and ESIM halos are different their phase-space density distributions are the same: \rho/\sigma^3 ~ r^{-\alpha}, with \alpha=1.875 over ~3 decades in radius. We use two approaches to try to explain this ``universal'' slope: (1) The Jeans equation analysis yields many insights, however, does not answer why \alpha=1.875. (2) The secondary infall model of the 1960's and 1970's, augmented by ``thermal motions'' of particles does predict that halos should have \alpha=1.875. However, this relies on assumptions of spherical symmetry and slow accretion. While for ESIM halos these assumptions are justified, they most certainly break down for simulated halos which forms hierarchically. We speculate that our argument may apply to an ``on-average'' formation scenario of halos within merger-driven numerical simulations, and thereby explain why \alpha=1.875 for NFW halos. Thus, \rho/\sigma^3 ~ r^{-1.875} may be a generic feature of violent relaxation.

We model the mass distributions of 40 high surface brightness spiral galaxies inside their optical radii, deriving parameters of mass models by matching the predicted velocities to observed velocity maps. We use constant mass-to-light disk and bulge models, and we have tried fits with no halo and with three different halo density profiles. The data require a halo in most, but not all, cases, while in others the best fit occurs with negligible mass in the luminous component, which we regard as unphysical. All three adopted halo profiles lead to fits of about the same quality, and our data therefore do not constrain the functional form of the halo profile. The halo parameters display large degeneracies for two of the three adopted halo functions, but the separate luminous and dark masses are better constrained. However, the fitted disk and halo masses vary substantially between the adopted halo models, indicating that even high quality 2-D optical velocity maps do not provide significant constraints on the dark matter content of a galaxy. We demonstrate that data from longslit observations are likely to provide still weaker constraints. We conclude that additional information is needed in order to constrain the separate disk and halo masses in a galaxy.

We investigate the impact of nonaxisymmetric structure on estimates of galaxy inclinations and position angles. A new minimization technique is used to obtain estimates of inclination and position angle from a global fit to either photometric or kinematic data. We discuss possible systematic uncertainties which are much larger than statistical uncertainties. Our investigation reveals that systematic uncertainties associated with fitting photometric data dominate the formal statistical uncertainties. For our sample of inclined galaxies, we estimate that nonaxisymmetric features introduce inclination and position angle uncertainties of approximately 5 degrees, on average. The magnitudes of these uncertainties weaken the arguments for intrinsically elliptical galaxy disks.

In an effort to more fully understand the variety of stellar distribution functions that can be used to construct models of realistic galaxies, the correlation integral method for orbit characterization that previously has been introduced by Grassberger & Procaccia (1983) and Carnevali & Santangelo (1984) is examined in considerable detail. The broad utility of the method is validated and demonstrated by applying it to orbits in a number of different, previously studied test cases (1D, 2D, and 3D; nonrotating and rotating). At the same time, the correlation integral method is compared and contrasted with other more traditional characterization tools such as Lyapunov exponents and surfaces of section. The method is then extended to orbits in a previously unexamined rotating, 3D bar potential. The correlation integral method is found to be a simple and reliable way to quantitatively categorize orbits in virtually any potential. It is recommended that it be broadly adopted as a tool for characterizing the properties of orbits and, by extension, stellar distribution functions, in all Hamiltonian dynamical systems.

In an effort to better understand the formation and evolution of barred galaxies, we have examined the properties of equatorial orbits in the effective potential of one model of a rapidly rotating, steady-state gas-dynamical bar that has been constructed via a self-consistent hydrodynamical simulation. Using a ``Restriction Hypothesis'' to determine initial conditions, we find that a significant fraction of orbits in this potential are quasi-ergodic and that regular orbits have a ``bowtie'' shape in contrast to the more typical x1 orbits. This bowtie orbit should give a boxy-peanut shape to such systems.