Punctuated equilibria and 1/f noise in a biological coevolution model with individual-based dynamics

We present a study by linear stability analysis and large-scale Monte Carlo simulations of a simple model of biological coevolution. Selection is provided through a reproduction probability that contains quenched, random interspecies interactions, while genetic variation is provided through a low mutation rate. Both selection and mutation act on individual organisms. Consistent with some current theories of macroevolutionary dynamics, the model displays intermittent, statistically self-similar behavior with punctuated equilibria. The probability density for the lifetimes of ecological communities is well approximated by a power law with exponent near -2, and the corresponding power spectral densities show 1/f noise (flicker noise) over several decades. The long-lived communities (quasi-steady states) consist of a relatively small number of mutualistically interacting species, and they are surrounded by a ``protection zone'' of closely related genotypes that have a very low probability of invading the resident community. The extent of the protection zone affects the stability of the community in a way analogous to the height of the free-energy barrier surrounding a metastable state in a physical system. Measures of biological diversity are on average stationary with no discernible trends, even over our very long simulation runs of approximately 3.4x10^7 generations.Comment: 20 pages RevTex. Minor revisions consistent with published versio

Scanning the phases of QCD with BRAHMS

BRAHMS has the ability to study relativistic heavy ion collisions from the final freeze-out of hadrons all the way back to the initial wave-function of the gold nuclei. This is accomplished by studying hadrons with a very wide range of momenta and angles. In doing so we can scan various phases of QCD, from a hadron gas, to a quark gluon plasma and perhaps to a color glass condensate.Comment: 8 pages, 6 figures, proceedings of plenary talk at Quark Matter 2004 conferenc

High Pt Hadron Spectra at High Rapidity

We report the measurement of charged hadron production at different pseudo-rapidity values in deuteron+gold as well as proton+proton collisions at $sqrt{s_{NN}}$ = 200GeV at RHIC. The nuclear modification factors $R_{dAu}$ and $R_{cp}$ are used to investigate new behaviors in the deuteron+gold system as function of rapidity and the centrality of the collisions respectively.Comment: Nine pages 4 figures to be published in the QM2004 Proceedings, typos corrected and one reference adde

Self-optimization, community stability, and fluctuations in two individual-based models of biological coevolution

We compare and contrast the long-time dynamical properties of two individual-based models of biological coevolution. Selection occurs via multispecies, stochastic population dynamics with reproduction probabilities that depend nonlinearly on the population densities of all species resident in the community. New species are introduced through mutation. Both models are amenable to exact linear stability analysis, and we compare the analytic results with large-scale kinetic Monte Carlo simulations, obtaining the population size as a function of an average interspecies interaction strength. Over time, the models self-optimize through mutation and selection to approximately maximize a community fitness function, subject only to constraints internal to the particular model. If the interspecies interactions are randomly distributed on an interval including positive values, the system evolves toward self-sustaining, mutualistic communities. In contrast, for the predator-prey case the matrix of interactions is antisymmetric, and a nonzero population size must be sustained by an external resource. Time series of the diversity and population size for both models show approximate 1/f noise and power-law distributions for the lifetimes of communities and species. For the mutualistic model, these two lifetime distributions have the same exponent, while their exponents are different for the predator-prey model. The difference is probably due to greater resilience toward mass extinctions in the food-web like communities produced by the predator-prey model.Comment: 26 pages, 12 figures. Discussion of early-time dynamics added. J. Math. Biol., in pres

Charged particle densities from Au+Au collisions at sqrt{s_{NN}}=130 GeV

We present charged particle densities as a function of pseudorapidity and collision centrality for the 197Au+197Au reaction at sqrt{s_{NN}}=130 GeV. An integral charged particle multiplicity of 3860+/-300 is found for the 5% most central events within the pseudorapidity range -4.7 <= eta <= 4.7. At mid-rapidity an enhancement in the particle yields per participant nucleon pair is observed for central events. Near to the beam rapidity, a scaling of the particle yields consistent with the ``limiting fragmentation'' picture is observed. Our results are compared to other recent experimental and theoretical discussions of charged particle densities in ultra-relativistic heavy-ion collisions.Comment: 14 pages, 4 figures; to be published in Phys. Lett.

Nuclear Modification Factor for Charged Pions and Protons at Forward Rapidity in Central Au+Au Collisions at 200 GeV

We present spectra of charged pions and protons in 0-10% central Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV at mid-rapidity ($y=0$) and forward pseudorapidity ($\eta=2.2$) measured with the BRAHMS experiment at RHIC. The spectra are compared to spectra from p+p collisions at the same energy scaled by the number of binary collisions. The resulting nuclear modification factors for central Au+Au collisions at both $y=0$ and $\eta=2.2$ exhibit suppression for charged pions but not for (anti-)protons at intermediate $p_T$. The $\bar{p}/\pi^-$ ratios have been measured up to $p_T\sim 3$ GeV/$c$ at the two rapidities and the results indicate that a significant fraction of the charged hadrons produced at intermediate $p_T$ range are (anti-)protons at both mid-rapidity and $\eta = 2.2$

Forward and midrapidity like-particle ratios from p+p collisions at sqrt(s)=200 GeV

We present a measurement of pi-\pi+, K-/K+ and pbar/p from p+p collisions at sqrt(s) = 20 0GeV over the rapidity range 0<y<3.4. For pT < 2.0 GeV/c we see no significant transverse momentum dependence of the ratios. All three ratios are independent of rapidity for y ~< 1.5 and then steadily decline from y ~ 1.5 to y ~ 3. The pi-\pi+ ratio is below unity for y > 2.0. The pbar/p ratio is very similar for p+p and 20% central Au+Au collisions at all rapidities. In the fragmentation region the three ratios seem to be independent of beam energy when viewed from the rest frame of one of the protons. Theoretical models based on quark-diquark breaking mechanisms overestimate the pbar/p ratio up to y ~< 3. Including additional mechanisms for baryon number transport such as baryon junctions leads to a better description of the data.Comment: 15 pages, 4 figures, uses elsart.sty. Changes to references and discussion based on referee comments, resubmitted to Phys. Lett.

The New Physics at RHIC. From Transparency to High pt_t Suppression

Heavy ion collisions at RHIC energies (Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV) exhibit significant new features as compared to earlier experiments at lower energies. The reaction is characterized by a high degree of transparency of the collisions partners leading to the formation of a baryon-poor central region. In this zone, particle production occurs mainly from the stretching of the color field. The initial energy density is well above the one considered necessary for the formation of the Quark Gluon Plasma, QGP. The production of charged particles of various masses is consistent with chemical and thermal equilibrium. Recently, a suppression of the high transverse momentum component of hadron spectra has been observed in central Au+Au collisions. This can be explained by the energy loss experienced by leading partons in a medium with a high density of unscreened color charges. In contrast, such high $p_t$ jets are not suppressed in d+Au collisions suggesting that the high $p_t$ suppression is not due to initial state effects in the ultrarelativistic colliding nuclei.Comment: 15 pages, 11 figures. to appear in Nucl. Physics A. Invited talk at 'Nucleus-Nucleus Collisions 2003' conference, Mosco

Anisotropic flow of charged hadrons, pions and (anti-)protons measured at high transverse momentum in Pb-Pb collisions at sNN=2.76\sqrt{s_{\rm NN}}=2.76 TeV

The elliptic, $v_2$, triangular, $v_3$, and quadrangular, $v_4$, azimuthal anisotropic flow coefficients are measured for unidentified charged particles, pions and (anti-)protons in Pb-Pb collisions at $\sqrt{s_{\rm NN}} = 2.76$ TeV with the ALICE detector at the Large Hadron Collider. Results obtained with the event plane and four-particle cumulant methods are reported for the pseudo-rapidity range $|\eta|<0.8$ at different collision centralities and as a function of transverse momentum, $p_{\rm T}$, out to $p_{\rm T}=20$ GeV/$c$. The observed non-zero elliptic and triangular flow depends only weakly on transverse momentum for $p_{\rm T}>8$ GeV/$c$. The small $p_{\rm T}$ dependence of the difference between elliptic flow results obtained from the event plane and four-particle cumulant methods suggests a common origin of flow fluctuations up to $p_{\rm T}=8$ GeV/$c$. The magnitude of the (anti-)proton elliptic and triangular flow is larger than that of pions out to at least $p_{\rm T}=8$ GeV/$c$ indicating that the particle type dependence persists out to high $p_{\rm T}$.Comment: 16 pages, 5 captioned figures, authors from page 11, published version, figures at http://aliceinfo.cern.ch/ArtSubmission/node/186

Centrality dependence of charged particle production at large transverse momentum in Pb-Pb collisions at sNN=2.76\sqrt{s_{\rm{NN}}} = 2.76 TeV

The inclusive transverse momentum ($p_{\rm T}$) distributions of primary charged particles are measured in the pseudo-rapidity range $|\eta|<0.8$ as a function of event centrality in Pb-Pb collisions at $\sqrt{s_{\rm{NN}}}=2.76$ TeV with ALICE at the LHC. The data are presented in the $p_{\rm T}$ range $0.15<p_{\rm T}<50$ GeV/$c$ for nine centrality intervals from 70-80% to 0-5%. The Pb-Pb spectra are presented in terms of the nuclear modification factor $R_{\rm{AA}}$ using a pp reference spectrum measured at the same collision energy. We observe that the suppression of high-$p_{\rm T}$ particles strongly depends on event centrality. In central collisions (0-5%) the yield is most suppressed with $R_{\rm{AA}}\approx0.13$ at $p_{\rm T}=6$-7 GeV/$c$. Above $p_{\rm T}=7$ GeV/$c$, there is a significant rise in the nuclear modification factor, which reaches $R_{\rm{AA}} \approx0.4$ for $p_{\rm T}>30$ GeV/$c$. In peripheral collisions (70-80%), the suppression is weaker with $R_{\rm{AA}} \approx 0.7$ almost independently of $p_{\rm T}$. The measured nuclear modification factors are compared to other measurements and model calculations.Comment: 17 pages, 4 captioned figures, 2 tables, authors from page 12, published version, figures at http://aliceinfo.cern.ch/ArtSubmission/node/284