Elliptic flow of resonances at RHIC: probing final state interactions and the structure of resonances

We propose the measurement of the elliptic flow of hadron resonances at the Relativistic Heavy Ion Collider as a tool to probe the amount of hadronic final state interactions for resonances at intermediate and large transverse momenta. This can be achieved by looking at systematic deviations of the measured flow coefficient $v_2$ from the scaling law given by the quark recombination formalism. Our method can be generalized to explore the structure of exotic particles, such as the recently found pentaquark $\Theta^+ (1540)$.Comment: 5 pages, 2 figures; v2: accepted version for publication in Physical Review C rapid communication

Jet modification in three dimensional fluid dynamics at next-to-leading twist

The modification of the single inclusive spectrum of high transverse momentum ($p_T$) pions emanating from an ultra-relativistic heavy-ion collision is investigated. The deconfined sector is modelled using a full three dimensional (3-D) ideal fluid dynamics simulation. Energy loss of high $p_T$ partons and the ensuing modification of their fragmentation is calculated within perturbative QCD at next-to-leading twist, where the magnitude of the higher twist contribution is modulated by the entropy density extracted from the 3-D fluid dynamics simulation. The nuclear modification factor ($R_{AA}$) for pions with a $p_T \geq 8$ GeV as a function of centrality as well as with respect to the reaction plane is calculated. The magnitude of contributions to the differential $R_{AA}$ within small angular ranges, from various depths in the dense matter is extracted from the calculation and demonstrate the correlation of the length integrated density and the $R_{AA}$ from a given depth. The significance of the mixed and hadronic phase to the overall magnitude of energy loss are explored.Comment: 5 pages, 4 figures, Revte

MAESTRO, CASTRO, and SEDONA -- Petascale Codes for Astrophysical Applications

Performing high-resolution, high-fidelity, three-dimensional simulations of Type Ia supernovae (SNe Ia) requires not only algorithms that accurately represent the correct physics, but also codes that effectively harness the resources of the most powerful supercomputers. We are developing a suite of codes that provide the capability to perform end-to-end simulations of SNe Ia, from the early convective phase leading up to ignition to the explosion phase in which deflagration/detonation waves explode the star to the computation of the light curves resulting from the explosion. In this paper we discuss these codes with an emphasis on the techniques needed to scale them to petascale architectures. We also demonstrate our ability to map data from a low Mach number formulation to a compressible solver.Comment: submitted to the Proceedings of the SciDAC 2010 meetin

Efficient Variable-Coefficient Finite-Volume Stokes Solvers

We investigate several robust preconditioners for solving the saddle-point linear systems that arise from spatial discretization of unsteady and steady variable-coefficient Stokes equations on a uniform staggered grid. Building on the success of using the classical projection method as a preconditioner for the coupled velocity-pressure system [B. E. Griffith, J. Comp. Phys., 228 (2009), pp. 75657595], as well as established techniques for steady and unsteady Stokes flow in the finite-element literature, we construct preconditioners that employ independent generalized Helmholtz and Poisson solvers for the velocity and pressure subproblems. We demonstrate that only a single cycle of a standard geometric multigrid algorithm serves as an effective inexact solver for each of these subproblems. Contrary to traditional wisdom, we find that the Stokes problem can be solved nearly as efficiently as the independent pressure and velocity subproblems, making the overall cost of solving the Stokes system comparable to the cost of classical projection or fractional step methods for incompressible flow, even for steady flow and in the presence of large density and viscosity contrasts. Two of the five preconditioners considered here are found to be robust to GMRES restarts and to increasing problem size, making them suitable for large-scale problems. Our work opens many possibilities for constructing novel unsplit temporal integrators for finite-volume spatial discretizations of the equations of low Mach and incompressible flow dynamics.Comment: Submitted to CiC

Hadronization in heavy ion collisions: Recombination and fragmentation of partons

We argue that the emission of hadrons with transverse momentum up to about 5 GeV/c in central relativistic heavy ion collisions is dominated by recombination, rather than fragmentation of partons. This mechanism provides a natural explanation for the observed constant baryon-to-meson ratio of about one and the apparent lack of a nuclear suppression of the baryon yield in this momentum range. Fragmentation becomes dominant at higher transverse momentum, but the transition point is delayed by the energy loss of fast partons in dense matter.Comment: 4 pages, 2 figures; v2: reference [8] added; v3: Eq.(2) corrected, two references added, version to appear in PR

Possible Resolutions of the D-Paradox

We propose possible ways of explaining the net charge event-by-event fluctuations in Au+Au collisions at the Relativistic Heavy Ion Collider within a quark recombination model. We discuss various methods of estimating the number of quarks at recombination and their implications for the predicted net charge fluctuations. We also discuss the possibility of diquark and quark-antiquark clustering above the deconfinement temperature.Comment: 5 pages, 2 figure