Popular models for describing the luminosity-density profiles of dynamically hot stellar systems (e.g., Jaffe, Hernquist, Dehnen) were constructed to match the deprojected form of de Vaucouleurs' $R^{1/4}$ light-profile. However, we now know that elliptical galaxies and bulges display a mass-dependent range of structural profiles. To compensate this, the model in Terzic & Graham was designed to closely match the deprojected form of Sersic $R^{1/n}$ light-profiles, including deprojected exponential light-profiles and galaxies with partially depleted cores. It is thus applicable for describing bulges in spiral galaxies, dwarf elliptical galaxies, both ``power-law'' and ``core'' elliptical galaxies, also dark matter halos formed from $\Lambda$CDM cosmological simulations. In this paper, we present a new family of triaxial density-potential-force triplets, which generalizes the spherical model reported in Terzic & Graham to three dimensions. If the (optional) power-law core is present, it is a 5-parameter family, while in the absence of the core it reduces to 3 parameters. The isodensity contours in the new family are stratified on confocal ellipsoids and the potential and forces are expressed in terms of integrals which are easy to evaluate numerically. We provide the community with a suite of numerical routines for orbit integration, which feature: optimized computations of potential and forces for this family; the ability to run simulations on parallel platforms; and modular and easily editable design.

We have recently shown that the 3-parameter density-profile model from Prugniel & Simien provides a better fit to simulated, galaxy- and cluster-sized, dark matter halos than an NFW-like model with arbitrary inner profile slope gamma (Paper I). By construction, the parameters of the Prugniel-Simien model equate to those of the Sersic R^{1/n} function fitted to the projected distribution. Using the Prugniel-Simien model, we are therefore able to show that the location of simulated (10^{12} M_sun) galaxy-sized dark matter halos in the <mu>_e-log(R_e) diagram coincides with that of brightest cluster galaxies, i.e., the dark matter halos appear consistent with the Kormendy relation defined by luminous elliptical galaxies. These objects are also seen to define the new, and equally strong, relation log(rho_e) = 0.5 - 2.5log(R_e), in which rho_e is the internal density at r=R_e. Simulated (10^{14.5} M_sun) cluster-sized dark matter halos and the gas component of real galaxy clusters follow the relation log(rho_e) = 2.5[1 - log(R_e)]. Given the shapes of the various density profiles, we are able to conclude that while dwarf elliptical galaxies and galaxy clusters can have dark matter halos with effective radii of comparable size to the effective radii of their baryonic component, luminous elliptical galaxies can not. For increasingly large elliptical galaxies, with increasingly large profile shapes `n', to be dark matter dominated at large radii requires dark matter halos with increasingly large effective radii compared to the effective radii of their stellar component.

We use techniques from nonparametric function estimation theory to extract the density profiles, and their derivatives, from a set of N-body dark matter halos. We consider halos generated from LCDM simulations of gravitational clustering, as well as isolated, spherical collapses. The logarithmic density slopes gamma = d(log rho)/d(log r) of the LCDM halos are found to vary as power-laws in radius, reaching values of gamma ~ -1 at the innermost resolved radii (~0.01 r_virial). This behavior is significantly different from that of broken power-law models like the NFW profile, but similar to that of models like de Vaucouleurs'. Accordingly, we compare the N-body density profiles with various parametric models to find which provide the best fit. We consider an NFW-like model with arbitrary inner slope; Dehnen & McLaughlin's anisotropic model; Einasto's model (identical in functional form to Sersic's model but fit to the space density); and the density model of Prugniel & Simien that was designed to match the deprojected form of Sersic's R^{1/n} law. Overall, the best-fitting model to the LCDM halos is Einasto's, although the Prugniel-Simien and Dehnen-McLaughlin models also perform well. With regard to the spherical collapse halos, both the Prugniel-Simien and Einasto models describe the density profiles well, with an rms scatter some four times smaller than that obtained with either the NFW-like model or the 3-parameter Dehnen-McLaughlin model. Finally, we confirm recent claims of a systematic variation in profile shape with halo mass.

We have recently shown that both the Prugniel-Simien model and Sersic's function (hereafter referred to as the Einasto model when applied to internal density profiles) describe simulated dark matter halos better than an NFW-like model with equal number of parameters. Here we provide analytical expressions for the logarithmic slopes of these models, and compare them with data from real galaxies. Depending on the Einasto parameters of the dark matter halo, one can expect an extrapolated, inner (0.01--1 kpc), logarithmic profile slope ranging from -0.2 to -1.5, with a typical value at 0.1 kpc around -0.7. Application of this (better fitting) model therefore alleviates some of the past disagreement with observations on this issue. We additionally provide useful expressions for the concentration and assorted scale radii: r_s, r_{-2}, r_e, R_e, r_virial, and r_max -- the radius where the circular velocity profile has its maximum value. We also present the circular velocity profiles and the radial behavior of rho(r)/sigma(r)^3 for both the Einasto and Prugniel-Simien models. We find this representation of the phase-space density profile to be well approximated by a power-law with slope slightly shallower than -2 near r=r_{-2}.

Popular models for describing the luminosity-density profiles of dynamically hot stellar systems (e.g., Jaffe, Hernquist, Dehnen) were constructed with the desire to match the deprojected form of an R^{1/4} light-profile. Real galaxies, however, are now known to have a range of different light-profile shapes that scale with mass. Consequently, although highly useful, the above models have implicit limitations, and this is illustrated here through their application to a number of real galaxy density profiles. On the other hand, the analytical density profile given by Prugniel & Simien (1997) closely matches the deprojected form of Sersic R^{1/n} light-profiles - including deprojected exponential light-profiles. It is thus applicable for describing bulges in spiral galaxies, dwarf elliptical galaxies, and both ordinary and giant elliptical galaxies. Here we provide simple equations, in terms of elementary and special functions, for the gravitational potential and force associated with this density profile. Furthermore, to match galaxies with partially depleted cores, and better explore the supermassive black hole / galaxy connection, we have added a power-law core to this density profile and derived similar expressions for the potential and force of this hybrid profile. Expressions for the mass and velocity dispersion, assuming isotropy, are also given. These models may also prove appropriate for describing the dark matter distribution in halos formed from LCDM cosmological simulations.

Orbits in a three-dimensional potential subjected to periodic driving, V(x^i,t)=[1+m_0 sin(omega t) V_0(x^i), divide naturally into two types, regular and chaotic, between which transitions are seemingly impossible. The chaotic orbits divide in turn into two types, apparently separated by entropy barriers, namely `sticky' orbit segments, which are `locked' to the driving frequency and exhibit little systematic energy diffusion, and `wildly' chaotic segments, which are not so locked and can exhibit significant energy diffusion. Attention focuses on how the relative abundance of these different orbit types and the transition rate between sticky and wildly chaotic orbits depends on amplitude m_0, and on the extent to which these quantities can be altered by weak friction and/or noise and by pseudo-random variations in the driving frequency, idealized as an Ornstein-Uhlenbeck process with (in general) nonzero autocorrelation time t_c. When, in the absence of perturbations, there exist large measures of both regular and chaotic orbits, the primary effect of weak noise is to increase the relative measure of chaotic orbits. Alternatively, when almost all the orbits are already chaotic, noise serves primarily to accelerate transitions from sticky to wildly chaotic behavior. The presence or absence of friction is unimportant and the details of the noise seem largely immaterial. In particular, there is only a weak, roughly logarithmic dependence on amplitude, and additive and multiplicative noise typically have virtually identical effects. By contrast, allowing for high frequency `noisy' variations in omega tends to weaken effects of the driving, decreasing the relative measure of chaotic orbits and suppressing large scale energy diffusion.

This paper focuses on symmetric potentials subjected to periodic driving. Four unperturbed potentials V_0(r) were considered, namely the Plummer potential and Dehnen potentials with \gamma=0.0, 0.5, and 1.0, each subjected to a time-dependence of the form V(r,t)=V_0(r)(1+m_0\sin(\omega t)). In each case, the orbits divide clearly into regular and chaotic, distinctions which appear absolute. In particular, transitions from regularity to chaos are seemingly impossible. Over finite time intervals, chaotic orbits subdivide into what can be termed `sticky' chaotic orbits, which exhibit no large scale secular changes in energy and remain trapped in the phase space region where they started; and `wildly chaotic' orbits, which do exhibit systematic drifts in energy as the orbits diffuse to inhabit different phase space regions. This latter distinction is not absolute, apparently corresponding instead to orbits penetrating a `leaky' phase space barrier. The three different orbit types can be identified simply in terms of the frequencies for which their Fourier spectra have the most power. An examination of the statistical properties of orbit ensembles as a function of driving frequency \omega allows one to identify the specific resonances that determine orbital structure. Attention focuses also on how, for fixed amplitude m_0, such quantities as the mean energy shift, the relative measure of chaotic orbits, and the mean value of the largest Lyapunov exponent vary with driving frequency \omega and how, for fixed \omega, the same quantities depend on m_0.

This paper focuses on the dynamical implications of close supermassive black hole binaries both as an example of resonant phase mixing and as a potential explanation of inversions and other anomalous features observed in the luminosity profiles of some elliptical galaxies. The presence of a binary comprised of black holes executing nearly periodic orbits leads to the possibility of a broad resonant coupling between the black holes and various stars in the galaxy. This can result in efficient chaotic phase mixing and, in many cases, systematic increases in the energies of stars and their consequent transport towards larger radii. Allowing for the presence of a supermassive black hole binary with plausible parameter values near the center of a spherical, or nearly spherical, galaxy characterised initially by a Nuker density profile enables one to reproduce in considerable detail the central surface brightness distributions of such galaxies as NGC 3706.

We use Schwarzschild's orbit superposition method to build self-consistent models of elliptical galaxies with scale-free potentials. Our exhaustive study of all physical shapes and central densities for galaxies with scale-free potentials establishes a relationship between the presence of chaos and self-consistency. The extent and the onset of chaos is studied by varying parameters of the model which are believed to be its inducers, such as the steepness of the central density cusp, flatness and triaxiality. We show that gravitational scattering of the central density cusp plays a dominant role in restricting the shapes of elliptical galaxies.

Recent work has shown that statistical arguments, seemingly well-justified in higher dimensions, can also be used to derive reasonable, albeit less accurate, estimates of the largest Lyapunov exponent ${\chi}$ in lower-dimensional Hamiltonian systems. This letter explores the detailed assumptions incorporated into these arguments. The predicted values of ${\chi}$ are insensitive to most of these details, which can in any event be relaxed straightforwardly, but {\em can} depend sensitively on the nongeneric form of the auto-correlation function characterising the time-dependence of an orbit. This dependence on dynamics implies a fundamental limitation to the application of thermodynamic arguments to such lower-dimensional systems.