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.Comment: AJ in pres

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}.Comment: AJ, in press. (Paper I can be found at astro-ph/0509417

Popular models for describing the luminosity-density profiles of dynamically hot stellar systems (e.g. Jaffe, Hernquist, Dehnen) were constructed with

The genetic algorithm (GA) is a powerful technique that implements the principles nature uses in biological evolution to optimize a multidimensional nonlinear problem. The GA works especially well for problems with a large number of local extrema, where traditional methods (such as conjugate gradient, steepest descent, and others) fail or, at best, underperform. The field of accelerator physics, among others, abounds with problems which lend themselves to optimization via GAs. In this paper, we report on the successful application of GAs in several problems related to the existing Continuous Electron Beam Accelerator Facility nuclear physics machine, the proposed Medium-energy Electron-Ion Collider at Jefferson Lab, and a radio frequency gun-based injector. These encouraging results are a step forward in optimizing accelerator design and provide an impetus for application of GAs to other problems in the field. To that end, we discuss the details of the GAs used, include a newly devised enhancement which leads to improved convergence to the optimum, and make recommendations for future GA developments and accelerator applications.

We have recently shown that the three-parameter density profile model from Prugniel & Simien provides a better fit to simulated galaxy- and cluster-sized dark matter halos than an Navarro-Frenk-White-like model with arbitrary inner profile slope γ (Paper I). By construction, the parameters of the Prugniel-Simien model equate to those of the Sérsic 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¹² M_⊙) galaxy-sized dark matter halos in the μ_e- log R_e diagram coincides with that of the 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 important, relation log(ρ_e) = 0.5 - 2.5 log(R_e), in which ρ_e is the internal density at r = R_e. Simulated (10^14.5 M_⊙) cluster-sized dark matter halos and the gas component of real galaxy clusters follow the relation log(ρ_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 cannot. 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 components.

We present a novel and quite general analysis of the interaction of a high-field chirped laser pulse and a relativistic electron, in which exquisite control of the spectral brilliance of the up-shifted Thomson-scattered photon is shown to be possible. Normally, when Thomson scattering occurs at high field strengths, there is ponderomotive line broadening in the scattered radiation. This effect makes the bandwidth too large for some applications and reduces the spectral brilliance. We show that such broadening can be corrected and eliminated by suitable frequency modulation of the incident laser pulse. Furthermore, we suggest a practical realization of this compensation idea in terms of a chirped-beam-driven free electron laser oscillator configuration and show that significant compensation can occur, even with the imperfect matching to be expected in these conditions.