Mass and Magnetic distributions in Self Gravitating Super Alfvenic Turbulence with AMR

In this work, we present the mass and magnetic distributions found in a recent Adaptive Mesh Refinement (AMR) MHD simulation of supersonic, \sa, self gravitating turbulence. Powerlaw tails are found in both volume density and magnetic field probability density functions, with $P(\rho) \propto \rho^{-1.67}$ and $P(B)\propto B^{-2.74}$. A power law is also found between magnetic field strength and density, with $B\propto \rho^{0.48}$, throughout the collapsing gas. The mass distribution of gravitationally bound cores is shown to be in excellent agreement with recent observation of prestellar cores. The mass to flux distribution of cores is also found to be in excellent agreement with recent Zeeman splitting measurements.Comment: 9 pages, 10 figures (3 color). Submitted to the Astrophysical Journa

A modified parallel tree code for N-body simulation of the Large Scale Structure of the Universe

N-body codes to perform simulations of the origin and evolution of the Large Scale Structure of the Universe have improved significantly over the past decade both in terms of the resolution achieved and of reduction of the CPU time. However, state-of-the-art N-body codes hardly allow one to deal with particle numbers larger than a few 10^7, even on the largest parallel systems. In order to allow simulations with larger resolution, we have first re-considered the grouping strategy as described in Barnes (1990) (hereafter B90) and applied it with some modifications to our WDSH-PT (Work and Data SHaring - Parallel Tree) code. In the first part of this paper we will give a short description of the code adopting the Barnes and Hut algorithm \cite{barh86} (hereafter BH), and in particular of the memory and work distribution strategy applied to describe the {\it data distribution} on a CC-NUMA machine like the CRAY-T3E system. In the second part of the paper we describe the modification to the Barnes grouping strategy we have devised to improve the performance of the WDSH-PT code. We will use the property that nearby particles have similar interaction list. This idea has been checked in B90, where an interaction list is builded which applies everywhere within a cell C_{group} containing a little number of particles N_{crit}. B90 reuses this interaction list for each particle $p \in C_{group}$ in the cell in turn. We will assume each particle p to have the same interaction list. Thus it has been possible to reduce the CPU time increasing the performances. This leads us to run simulations with a large number of particles (N ~ 10^7/10^9) in non-prohibitive times.Comment: 13 pages and 7 Figure

Initial data transients in binary black hole evolutions

We describe a method for initializing characteristic evolutions of the Einstein equations using a linearized solution corresponding to purely outgoing radiation. This allows for a more consistent application of the characteristic (null cone) techniques for invariantly determining the gravitational radiation content of numerical simulations. In addition, we are able to identify the {\em ingoing} radiation contained in the characteristic initial data, as well as in the initial data of the 3+1 simulation. We find that each component leads to a small but long lasting (several hundred mass scales) transient in the measured outgoing gravitational waves.Comment: 18 pages, 4 figure

Unambiguous determination of gravitational waveforms from binary black hole mergers

Gravitational radiation is properly defined only at future null infinity (\scri), but in practice it is estimated from data calculated at a finite radius. We have used characteristic extraction to calculate gravitational radiation at \scri for the inspiral and merger of two equal mass non-spinning black holes. Thus we have determined the first unambiguous merger waveforms for this problem. The implementation is general purpose, and can be applied to calculate the gravitational radiation, at \scri, given data at a finite radius calculated in another computation.Comment: 4 pages, 3 figures, published versio

Head-on collisions of unequal mass black holes in D=5 dimensions

We study head-on collisions of unequal mass black hole binaries in D=5 space-time dimensions, with mass ratios between 1:1 and 1:4. Information about gravitational radiation is extracted by using the Kodama-Ishibashi gauge-invariant formalism and details of the apparent horizon of the final black hole. For the first time, we present waveforms, total integrated energy and momentum for this process. Our results show surprisingly good agreement, within 5% or less, with those extrapolated from linearized, point-particle calculations. Our results also show that consistency with the area theorem bound requires that the same process in a large number of spacetime dimensions must display new features.Comment: 10 pages, 5 figures, RevTex4. v2: Published versio

A Comparison of Simple Mass Estimators for Galaxy Clusters

High-resolution N-body simulations are used to investigate systematic trends in the mass profiles and total masses of clusters as derived from 3 simple estimators: (1) the weak gravitational lensing shear field under the assumption of an isothermal cluster potential, (2) the dynamical mass obtained from the measured velocity dispersion under the assumption of an isothermal cluster potential, and (3) the classical virial estimator. The clusters consist of order 2.5e+05 particles of mass m_p \simeq 10^{10} \Msun, have triaxial mass distributions, and significant substructure exists within their virial radii. Not surprisingly, the level of agreement between the mass profiles obtained from the various estimators and the actual mass profiles is found to be scale-dependent. The virial estimator yields a good measurement of the total cluster mass, though it is systematically underestimated by of order 10%. This result suggests that, at least in the limit of ideal data, the virial estimator is quite robust to deviations from pure spherical symmetry and the presence of substructure. The dynamical mass estimate based upon a measurement of the cluster velocity dispersion and an assumption of an isothermal potential yields a poor measurement of the total mass. The weak lensing estimate yields a very good measurement of the total mass, provided the mean shear used to determine the equivalent cluster velocity dispersion is computed from an average of the lensing signal over the entire cluster (i.e. the mean shear is computed interior to the virial radius). [abridged]Comment: Accepted for publication in The Astrophysical Journal. Complete paper, including 3 large colour figures can also be obtained from http://bu-ast.bu.edu/~brainerd/preprints