Cluster Mass Estimate and a Cusp of the Mass Density Distribution in Clusters of Galaxies

We study density cusps in the center of clusters of galaxies to reconcile X-ray mass estimates with gravitational lensing masses. For various mass density models with cusps we compute X-ray surface brightness distribution, and fit them to observations to measure the range of parameters in the density models. The Einstein radii estimated from these density models are compared with Einstein radii derived from the observed arcs for Abell 2163, Abell 2218, and RX J1347.5-1145. The X-ray masses and lensing masses corresponding to these Einstein radii are also compared. While steeper cusps give smaller ratios of lensing mass to X-ray mass, the X-ray surface brightnesses estimated from flatter cusps are better fits to the observations. For Abell 2163 and Abell 2218, although the isothermal sphere with a finite core cannot produce giant arc images, a density model with a central cusp can produce a finite Einstein radius, which is smaller than the observed radii. We find that a total mass density profile which declines as $\sim r^{-1.4}$ produces the largest radius in models which are consistent with the X-ray surface brightness profile. As the result, the extremely large ratio of the lensing mass to the X-ray mass is improved from 2.2 to 1.4 for Abell 2163, and from 3 to 2.4 for Abell 2218. For RX J1347.5-1145, which is a cooling flow cluster, we cannot reduce the mass discrepancy.Comment: 23 pages, 10 figures, Latex, uses aasms4.sty, accepted for publication in ApJ, Part

Star cluster ecology IVa: Dissection of an open star cluster---photometry

The evolution of star clusters is studied using N-body simulations in which the evolution of single stars and binaries are taken self-consistently into account. Initial conditions are chosen to represent relatively young Galactic open clusters, such as the Pleiades, Praesepe and the Hyades. The calculations include a realistic mass function, primordial binaries and the external potential of the parent Galaxy. Our model clusters are generally significantly flattened in the Galactic tidal field, and dissolve before deep core collapse occurs. The binary fraction decreases initially due to the destruction of soft binaries, but increases later because lower mass single stars escape more easily than the more massive binaries. At late times, the cluster core is quite rich in giants and white dwarfs. There is no evidence for preferential evaporation of old white dwarfs, on the contrary the formed white dwarfs are likely to remain in the cluster. Stars tend to escape from the cluster through the first and second Lagrange points, in the direction of and away from the Galactic center. Mass segregation manifests itself in our models well within an initial relaxation time. As expected, giants and white dwarfs are much more strongly affected by mass segregation than main-sequence stars. Open clusters are dynamically rather inactive. However, the combined effect of stellar mass loss and evaporation of stars from the cluster potential drives its dissolution on a much shorter timescale than if these effects are neglected. The often-used argument that a star cluster is barely older than its relaxation time and therefore cannot be dynamically evolved is clearly in error for the majority of star clusters.Comment: reduced abstract, 33 pages (three separate color .jpg figures), submitted to MNRA

Deviations from Berry--Robnik Distribution Caused by Spectral Accumulation

By extending the Berry--Robnik approach for the nearly integrable quantum systems,\cite{[1]} we propose one possible scenario of the energy level spacing distribution that deviates from the Berry--Robnik distribution. The result described in this paper implies that deviations from the Berry--Robnik distribution would arise when energy level components show strong accumulation, and otherwise, the level spacing distribution agrees with the Berry--Robnik distribution.Comment: 4 page

Parallelization, Special Hardware and Post-Newtonian Dynamics in Direct N - Body Simulations

The formation and evolution of supermassive black hole (SMBH) binaries during and after galaxy mergers is an important ingredient for our understanding of galaxy formation and evolution in a cosmological context, e.g. for predictions of cosmic star formation histories or of SMBH demographics (to predict events that emit gravitational waves). If galaxies merge in the course of their evolution, there should be either many binary or even multiple black holes, or we have to find out what happens to black hole multiples in galactic nuclei, e.g. whether they come sufficiently close to merge resulting from emission of gravitational waves, or whether they eject each other in gravitational slingshot interactions

Static cylindrically symmetric spacetimes

We prove existence of static solutions to the cylindrically symmetric Einstein-Vlasov system, and we show that the matter cylinder has finite extension. The same results are also proved for a quite general class of equations of state for perfect fluids coupled to the Einstein equations, extending the class of equations of state considered in \cite{BL}. We also obtain this result for the Vlasov-Poisson system.Comment: Added acknowledgemen

Fringe Fields and Dynamic Aperture in the FNAL Muon Storage Ring

Quadrupole fringe fields can limit the dynamic aperture of muon storage rings.Using the computer code COSY INFINITY for particle tracking and normal-form analysis, we evaluate the importance of fringe fields in the FNAL muon storage ring, and identify the regions of the machine where they are most critical. Dynamic aperture and linear tune shifts with amplitude are calculated considering an ideal machine without any errors or misalignments. We also explore the efficiency of various nonlinear correction schemes, study the momentum acceptance, and evaluate the spin decoherence over the transverse phase space

Spherically symmetric relativistic stellar structures

We investigate relativistic spherically symmetric static perfect fluid models in the framework of the theory of dynamical systems. The field equations are recast into a regular dynamical system on a 3-dimensional compact state space, thereby avoiding the non-regularity problems associated with the Tolman-Oppenheimer-Volkoff equation. The global picture of the solution space thus obtained is used to derive qualitative features and to prove theorems about mass-radius properties. The perfect fluids we discuss are described by barotropic equations of state that are asymptotically polytropic at low pressures and, for certain applications, asymptotically linear at high pressures. We employ dimensionless variables that are asymptotically homology invariant in the low pressure regime, and thus we generalize standard work on Newtonian polytropes to a relativistic setting and to a much larger class of equations of state. Our dynamical systems framework is particularly suited for numerical computations, as illustrated by several numerical examples, e.g., the ideal neutron gas and examples that involve phase transitions.Comment: 23 pages, 25 figures (compressed), LaTe

Harvesting graphics power for MD simulations

We discuss an implementation of molecular dynamics (MD) simulations on a graphic processing unit (GPU) in the NVIDIA CUDA language. We tested our code on a modern GPU, the NVIDIA GeForce 8800 GTX. Results for two MD algorithms suitable for short-ranged and long-ranged interactions, and a congruential shift random number generator are presented. The performance of the GPU's is compared to their main processor counterpart. We achieve speedups of up to 80, 40 and 150 fold, respectively. With newest generation of GPU's one can run standard MD simulations at 10^7 flops/$.Comment: 12 pages, 5 figures. Submitted to Mol. Si