Shear-Wave Splitting and Mantle Flow Beneath the Colorado Plateau and its Boundary with the Great Basin

Shear-wave splitting measurements from SKS and SKKS phases show fast polarization azimuths that are subparallel to North American absolute plate motion within the central Rio Grande Rift (RGR) and Colorado Plateau (CP) through to the western rim of the CP, with anisotropy beneath the CP and central RGR showing a remarkably consistent pattern with a mean fast azimuth of 4 degrees +/- degrees 6 E of N. Approaching the rim from the southeast, fast anisotropic directions become north-northeast-south-southwest (NNE-SSW), rotate counter clockwise to north-south in the CP-GB transition, and then to NNW-SSE in the western Great Basin ( GB). This change is coincident with uppermost mantle S-wave velocity perturbations that vary from +4% beneath the western CP and the eastern edge of the Marysvale volcanic field to about -8% beneath the GB. Corresponding delay times average 1.5 sec beneath the central CP, decrease to approximately 0.8 sec near the CP-GB transition, and increase to about 1.2 sec beneath the GB. For the central CP, we suggest anisotropy predominantly controlled by North American plate motion above the asthenosphere. The observed pattern of westward-rotating anisotropy from the western CP through the CP-GB transition may be influenced to asthenospheric flow around a CP lithospheric keel and/or by vertical flow arising from edge-driven small-scale convection. The anisotropic transition from the CP to the GB thus marks a first-order change from absolute plate motion dominated lithosphere-asthenosphere shear to a new regime controlled by regional flow processes. The NNW-SSE anisotropic fast directions of split SKS waves in the eastern GB area are part of a broad circular pattern of seismic anisotropic fast direction in the central GB that has recently been hypothesized to be due to toroidal flow around the sinking Juan de Fuca-Gorda slab.National Science Foundation EAR 9706094, 9707188, 9707190, 0207812Los Alamos National Laboratory Institute of Geophysics and Planetary PhysicsNational Science Foundation Cooperative EAR-000430Department of Energy National Nuclear Security AdministrationGeological Science

Antibaryon density in the central rapidity region of a heavy ion collision

We consider (anti-)baryons production in heavy ion collisions as production of topological defects during the chiral phase transition. Non-zero quark masses which explicitly break chiral symmetry supress the (anti-)baryon density. Hardly any (anti-)baryons will be produced in the central rapidity region of a heavy ion collision.Comment: 3 pages in RevTex, 3 .ps file

Routine Microsecond Molecular Dynamics Simulations with AMBER on GPUs. 1. Generalized Born

We present an implementation of generalized Born implicit solvent all-atom classical molecular dynamics (MD) within the AMBER program package that runs entirely on CUDA enabled NVIDIA graphics processing units (GPUs). We discuss the algorithms that are used to exploit the processing power of the GPUs and show the performance that can be achieved in comparison to simulations on conventional CPU clusters. The implementation supports three different precision models in which the contributions to the forces are calculated in single precision floating point arithmetic but accumulated in double precision (SPDP), or everything is computed in single precision (SPSP) or double precision (DPDP). In addition to performance, we have focused on understanding the implications of the different precision models on the outcome of implicit solvent MD simulations. We show results for a range of tests including the accuracy of single point force evaluations and energy conservation as well as structural properties pertainining to protein dynamics. The numerical noise due to rounding errors within the SPSP precision model is sufficiently large to lead to an accumulation of errors which can result in unphysical trajectories for long time scale simulations. We recommend the use of the mixed-precision SPDP model since the numerical results obtained are comparable with those of the full double precision DPDP model and the reference double precision CPU implementation but at significantly reduced computational cost. Our implementation provides performance for GB simulations on a single desktop that is on par with, and in some cases exceeds, that of traditional supercomputers