Search

Hot QCD physics studies the nuclear strong force under extreme temperature and densities. Experimentally these conditions are achieved via high-energy collisions of heavy ions at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC). In the past decade, a unique and substantial suite of data was collected at RHIC and the LHC, probing hydrodynamics at the nucleon scale, the temperature dependence of the transport properties of quark-gluon plasma, the phase diagram of nuclear matter, the interaction of quarks and gluons at different scales and much more. This document, as part of the 2023 nuclear science long range planning process, was written to review the progress in hot QCD since the 2015 Long Range Plan for Nuclear Science, as well as highlight the realization of previous recommendations, and present opportunities for the next decade, building on the accomplishments and investments made in theoretical developments and the construction of new detectors. Furthermore, this document provides additional context to support the recommendations voted on at the Joint Hot and Cold QCD Town Hall Meeting, which are reported in a separate document., Comment: 190 pages, 69 figures

We present results of hydrodynamic modelling of Au-Au collisions from $\sqrt{s_{NN}}$ = 7.7 to 200 GeV. Our simulations have three novel components. Firstly, we use a Linear EXtrapolation of Ultrarelativistic nucleon-nucleon Scattering to nucleus-nucleus collisions (LEXUS) inspired Monte-Carlo initial state model. Secondly, we use a crossover equation of state at finite baryon densities without a critical point. Finally, we use departure functions derived from the quasiparticle theory of transport coefficients for hadronic matter at non-zero baryon densities., Comment: 19 pages, 12 figures

Whether Quantum Chromodynamics (QCD) exhibits a phase transition at finite temperature and density is an open question. It is important for hydrodynamic modeling of heavy ion collisions and neutron star mergers. Lattice QCD simulations have definitively shown that the transition from hadrons to quarks and gluons is a crossover when the baryon chemical potential is zero or small. We combine the parametric scaling equation of state, usually associated with the 3D Ising model, with a background equation of state based on a smooth crossover from hadrons to quarks and gluons. Comparison to experimental data from the Beam Energy Scan II at the Relativistic Heavy Ion Collider or in heavy ion experiments at other accelerators may allow the critical exponents and amplitudes in the scaling equation of state to be determined for QCD if a critical point exists., Comment: 17 pages, 9 figures

Lattice QCD simulations have shown unequivocally that the transition from hadrons to quarks and gluons is a crossover when the baryon chemical potential is zero or small. Many model calculations predict the existence of a critical point at a value of the chemical potential where current lattice simulations are unreliable. We show how to embed a critical point in a smooth background equation of state so as to yield the critical exponents and critical amplitude ratios expected of a transition in the same universality class as the liquid-gas phase transition and the 3D Ising model. The resulting equation of state has parameters which may be inferred by hydrodynamic modeling of heavy ion collisions in the Beam Energy Scan II at the Relativistic Heavy Ion Collider or in experiments at other accelerators., Comment: 4 pages, 4 figures

The question of whether quark matter exists in neutron stars is a long standing one. Generally one finds that a first order phase transition from baryons to quarks softens the equation of state so much that the star would collapse into a black hole. We consider a crossover equation of state, similar to the crossover that is found in lattice QCD studies at finite temperature and zero or small baryon chemical potentials. We find that with reasonable parameters it may be possible to support neutron stars up to about 2.2 solar masses. In that case 1 to 10% of the pressure would be contributed by quark matter in the central core of the highest mass stars., Comment: 4 pages, 5 figures

We propose a dynamical mechanism which provides an interaction between the spins of hyperons and antihyperons and the vorticity of the baryon current in noncentral high energy nuclear collisions. The interaction is mediated by massive vector and scalar bosons which is well known to describe the nuclear spin-orbit force. It follows from the Foldy-Wouthuysen transformation and leads to a strong-interaction Zeeman effect. The interaction may explain the difference in polarizations of $\Lambda$ and $\bar{\Lambda}$ hyperons as measured by the STAR Collaboration at the BNL Relativistic Heavy Ion Collider. The signs and magnitudes of the meson-baryon couplings are closely connected to the binding energies of hypernuclei and to the abundance of hyperons in neutron stars., Comment: 5 pages, 2 figures, version accepted for publication in PRC Rapids

In simulations of high energy heavy ion collisions that employ viscous hydrodynamics, single particle distributions are distorted from their thermal equilibrium form due to gradients in the flow velocity. These are closely related to the formulas for the shear and bulk viscosities in the quasi-particle approximation. Distorted single particle distributions are now commonly used to calculate the emission of photons and dilepton pairs, and in the late stage to calculate the conversion of a continuous fluid to individual particles. We show how distortions of the single particle distribution functions due to both shear and bulk viscous effects can be done rigorously in the quasi-particle approximation and illustrate it with the linear $\sigma$ model at finite temperature., Comment: 5 pages, 3 figures. This is the published version which takes into account comments from the referee

We have constructed an equation of state which smoothly interpolates between an excluded volume hadron resonance gas at low energy density to a plasma of quarks and gluons at high energy density. This crossover equation of state agrees very well with lattice calculations at both zero and nonzero baryon chemical potential. We use it to compute the variance, skewness, and kurtosis of fluctuations of baryon number, and compare to measurements of proton number fluctuations in central Au-Au collisions as measured by the STAR collaboration in a beam energy scan at the Relativistic Heavy Ion Collider. The crossover equation of state can reproduce the data if the fluctuations are frozen out at temperatures well below than the average chemical freeze-out., Comment: 5 pages, 7 figures. arXiv admin note: substantial text overlap with arXiv:1506.03408

We develop a flexible quasiparticle theory of transport coefficients of hot hadronic matter at finite baryon density. We begin with a hadronic quasiparticle model which includes a scalar and a vector mean field. Quasiparticle energies and the mean fields depend on temperature and baryon chemical potential. Starting with the quasiparticle dispersion relation, we derive the Boltzmann equation and use the Chapman-Enskog expansion to derive formulas for the shear and bulk viscosities and thermal conductivity. We obtain both relaxation time approximation formulas and more general integral equations. Throughout the work, we explicitly enforce the Landau-Lifshitz conditions of fit and ensure the theory is thermodynamically self-consistent. The derived formulas should be useful for predicting the transport coefficients of the hadronic phase of matter produced in heavy-ion collisions at the Relativistic Heavy Ion Collider (RHIC) and at other accelerators., Comment: 25 pages

Nuclei colliding at very high energy create a strong, quasi-classical gluon field during the initial phase of their interaction. We present an analytic calculation of the initial space-time evolution of this field in the limit of very high energies using a formal recursive solution of the Yang-Mills equations. We provide analytic expressions for the initial chromo-electric and chromo-magnetic fields and for their energy-momentum tensor. In particular, we discuss event-averaged results for energy density and energy flow as well as for longitudinal and transverse pressure of this system. For example, we find that the ratio of longitudinal to transverse pressure very early in the system behaves as $p_L/p_T = -[1-\frac{3}{2a}(Q\tau)^2]/[1-\frac{1}{a}(Q\tau)^2]+\mathcal{O}(Q\tau)^4$ where $\tau$ is the longitudinal proper time, $Q$ is related to the saturation scales $Q_s$ of the two nuclei, and $a = \ln (Q^2/\hat{m}^2)$ with $\hat m$ a scale to be defined later. Our results are generally applicable if $\tau \lesssim 1/Q$. As already discussed in a previous paper, the transverse energy flow $S^i$ of the gluon field exhibits hydrodynamic-like contributions that follow transverse gradients of the energy density $\nabla^i \varepsilon$. In addition, a rapidity-odd energy flow also emerges from the non-abelian analog of Gauss' Law and generates non-vanishing angular momentum of the field. We will discuss the space-time picture that emerges from our analysis and its implications for observables in heavy ion collisions., Comment: 26 pages, 9 figures

Fluctuations of the proton number distribution in central Au-Au collisions have been measured by the STAR collaboration in a beam energy scan at the Relativistic Heavy Ion Collider (RHIC). The motivation is a search for evidence of a critical point in the equation of state. It was found that the skewness and kurtosis display an interesting energy dependence. We compare these measurements to an equation of state which smoothly interpolates between an excluded volume hadron resonance gas at low energy density to a perturbative plasma of quarks and gluons at high energy density. This crossover equation of state agrees very well with the lattice QCD equation of state. The crossover equation of state can reproduce the data if the fluctuations are frozen at a temperature significantly lower than the average chemical freeze-out., Comment: 17 pages, 10 figures

We construct a model for baryon diffusion which has the desired properties of analyticity and causality. The model also has the desired property that the noise correlation function is not a Dirac delta function in space and time. The model depends on three time constants in addition to the diffusion constant. This description can be incorporated into 2nd order viscous hydrodynamical models of heavy ion collisions., Comment: 5 pages, 2 figures, Proceedings of Quark Matter 2014

The fluctuation-dissipation theorem requires the presence of thermal noise in viscous fluids. The time and length scales of heavy ion collisions are small enough so that the thermal noise can have a measurable effect on observables. Thermal noise is included in numerical simulations of high energy lead-lead collisions, increasing average values of the momentum eccentricity and contributing to its event by event fluctuations., Comment: 13 pages, 4 figures

We match three hadronic equations of state at low energy densities to a perturbatively computed equation of state of quarks and gluons at high energy densities. One of them includes all known hadrons treated as point particles, which approximates attractive interactions among hadrons. The other two include, in addition, repulsive interactions in the form of excluded volumes occupied by the hadrons. A switching function is employed to make the crossover transition from one phase to another without introducing a thermodynamic phase transition. A chi-square fit to accurate lattice calculations with temperature $100 < T < 1000$ MeV determines the parameters. These parameters quantify the behavior of the QCD running gauge coupling and the hard core radius of protons and neutrons, which turns out to be $0.62 \pm 0.04$ fm. The most physically reasonable models include the excluded volume effect. Not only do they include the effects of attractive and repulsive interactions among hadrons, but they also achieve better agreement with lattice QCD calculations of the equation of state. The equations of state constructed in this paper do not result in a phase transition, at least not for the temperatures and baryon chemical potentials investigated. It remains to be seen how well these equations of state will represent experimental data on high energy heavy ion collisions when implemented in hydrodynamic simulations., Comment: 19 pages, 9 figures

We construct a model of baryon diffusion which has the desired properties of causality and analyticity. The model also has the desired property of colored noise, meaning that the noise correlation function is not a Dirac delta function in space and time; rather, it depends on multiple time and length constants. The model can readily be incorporated in 3+1 dimensional second order viscous hydro-dynamical models of heavy ion collisions, which is particularly important at beam energies where the baryon density is large., Comment: 20 pages, 5 figures

Using the Bjorken expansion model we study the effect of intrinsic hydrodynamic noise on the correlations observed in heavy-ion collisions., Comment: Proceedings of Quark Matter 2012, (August 13-18, 2012, Washington DC), 4 pages, 3 figures

We develop the relativistic theory of hydrodynamic fluctuations for application to high energy heavy ion collisions. In particular, we investigate their effect on the expanding boost-invariant (Bjorken) solution of the hydrodynamic equations. We discover that correlations over a long rapidity range are induced by the propagation of the sound modes. Due to the expansion, the dispersion law for these modes is non-linear and attenuated even in the limit of zero viscosity. As a result, there is a non-dissipative wake behind the sound front which is generated by any instantaneous point-like fluctuation. We evaluate the two-particle correlators using the initial conditions and hydrodynamic parameters relevant for heavy-ion collisions at RHIC and LHC. In principle these correlators can be used to obtain information about the viscosities because the magnitudes of the fluctuations are directly proportional to them., Comment: 39 pages, 6 figures; references added

A theoretical framework for the calculation of shear and bulk viscosities of hadronic matter at finite temperature is presented. The framework is based on the quasi-particle picture. It allows for an arbitrary number of hadron species with point-like interactions, and allows for both elastic and inelastic collisions. Detailed balance is ensured. The particles have temperature dependent masses arising from mean field or potential effects, which maintains self-consistency between the equation of state and the transport coefficients. As an example, we calculate the shear and bulk viscosity in the linear $\sigma$ model. The ratio of shear viscosity to entropy density shows a minimum in the vicinity of a rapid crossover transition, while the ratio of bulk viscosity to entropy density shows a maximum., Comment: 45 pages

Soft-wall models in AdS/QCD generally have dilaton and scalar fields that vary with the fifth-dimension coordinate. These fields can be parameterized to yield hadron mass spectra with linear radial trajectories and to incorporate spontaneous breaking of chiral symmetry. We show how to construct scalar potentials which lead to such solutions., Comment: 21 pages. Version 3 matches that accepted for publication. Typos have been corrected and references added. Substantial discussion has been added in order to clarify our approach and the issue of the dilaton mass in Sec. 4. Central results and conclusions remain unchanged.

The thermodynamics and critical exponents and amplitudes of high temperature and dense matter near the chiral critical point are studied. The parameterized equation of state matches on to that calculated with lattice QCD at zero chemical potential and to the known properties of nuclear matter at zero temperature. The extent to which finite size effects wash out the phase separation near the critical point is determined., Comment: 32 pages, 8 figures

We use the linear sigma model with two flavors of quarks to study the phase diagram at finite temperature and baryon chemical potential as a function of the vacuum pion mass. Our calculations include thermal fluctuations of both the bosonic and fermionic fields. An interesting phase structure occurs that results in two, one, or no critical points depending on the value of the pion mass., Comment: Proceedings of "Critical Point and Onset of Deconfinement 2009", Brookhaven National Laboratory. 7 pages, 5 figures

We employ the linear sigma model to study the chiral dynamics of two flavors of quarks at finite temperature and density. Our calculations include thermal fluctuations of both the bosonic and fermionic fields. In particular, we determine the phase diagram in the plane of temperature and baryon chemical potential as a function of the pion mass. An interesting phase structure occurs that results in zero, one, or two critical points depending on the value of the vacuum pion mass., Comment: 4 pages, Quark Matter 2009 proceedings

We study the shear mode in the gauge/gravity correspondence at finite temperature. First, we confirm the general formula for the shear viscosity in an arbitrary background metric which includes a black hole in the fifth dimension. We then derive a general formula for the shear mode relaxation time which appears in the theory of relativistic dissipative fluid dynamics; it agrees with known expressions in the limit of conformal fields. These results may be useful in relativistic viscous fluid descriptions of high energy nuclear collisions at RHIC and LHC., Comment: Modifications based on helpful feedback from colleagues and referee

Electric and color neutral solutions, and the critical conditions for the formation of gapless color superconductors, are investigated in $K^0$ condensed color-flavor locked quark matter for nonzero strange quark mass. We show that as the strange quark mass increases, gapless modes for up-strange quark pairing occur first, followed by down-strange quark pairing. The behavior of the gaps, the dispersion relations, and the thermodynamic potential are all found as functions of the strange quark mass on the basis of a Nambu Jona-Lasinio type model. To a high degree of accuracy, they are presented as relatively simple elementary functions. This allows for easy computation for any reasonable range of baryon chemical potential and strange quark mass., Comment: 29 pages, 3 figures, 4 tables

Photons and dileptons are being used to probe the properties of nuclear and quark-gluon matter at high energy densities. This is an area where theory and experiment are driving each other to obtain solid results. However, it is important to clearly separate the assumptions and conclusions concerning the correlation and response functions of a system in thermal equilibrium from the space-time dynamics used to model the evolution of the matter created in the nuclear collision., Comment: 9 pages, 10 figures, Proceedings of Hard Probes 2006

The electric/color neutral solution and the critical conditions for gapless formation are investigated in the $K^0$ condensed color-flavor locked matter. We point out that there exist no longer gapless modes for down-strange quark pairing while the gapless phenomenon for up-strange one is dominated in the $K^0$ condensed environment. In a model-independent way, the phase transition to the resulting gapless phase is found to be of first-order. The novel phase structure implies that the chromomagnetic instability happens in the previous-predicted gapless phase might be removed at least partly., Comment: 2 figures

We calculate the classical gluon field created at early times in collisions of large nuclei at high energies. We find that the field is dominated by the longitudinal chromoelectric and chromomagnetic components. We estimate the initial energy density of this gluon field to be approximately 260 GeV/fm$^3$ at RHIC., Comment: 4 pages, 1 fig; v2: some new material included, discussion improved, major conclusions unchanged

The spectrum of positrons and antiprotons produced by the outflow of high temperature matter surrounding microscopic black holes is calculated for energies between five GeV and the Planck energy. The results may be applicable for the last few hours and minutes of a microscopic black hole's lifetime., Comment: 7 pages, 3 figures

We derive a general formula for the correlation function of two identical particles with the inclusion of multiple elastic scatterings in the medium in which the two particles are produced. This formula involves the propagator of the particle in the medium. As illustration of the effect we apply the formula to the special case where the scatterers are static, localized 2-body potentials. In this illustration both $R^2_{\rm side}$ and $R^2_{\rm out}$ are increased by an amount proportional to the square of the spatial density of scatterers and to the differential cross section. Specific numbers are used to show the expected magnitude of the rescattering effect on kaon interferometry., Comment: 15 pages, 4 figures

It is shown that Dirac-type neutrinos display BCS superfluidity for any nonzero mass. The Cooper pairs are formed by attractive scalar Higgs boson exchange between left- and right-handed neutrinos; in the standard SU(2)xU(1) theory, right-handed neutrinos do not couple to any other boson. The value of the gap, the critical temperature, and the Pippard coherence length are calculated for arbitrary values of the neutrino mass and chemical potential. Although such a superfluid could conceivably exist, detecting it would be a major challenge., Comment: This is the version published in PRL

We develop and solve a network of rate equations for the production of baryons and anti-baryons in high energy nuclear collisions. We include all members of the baryon octet and decuplet and allow for transformations among them. This network is solved during a relativistic 2+1 hydrodynamical expansion of the of the hot matter created in the collision. As an application we compare to the number of protons, lambdas, negative cascades, and omega baryons measured at mid-rapidity in central collisions of gold nuclei at 65 GeV per nucleon at the Relativistic Heavy Ion Collider (RHIC)., Comment: 24 pages, 9 figures

The spectrum of neutrinos produced by the outflow of high temperature matter surrounding microscopic black holes is calculated for neutrino energies between one GeV and the Planck energy. The results may be applicable for the last few hours and minutes of a microscopic black hole's lifetime., Comment: 13 pages, 4 figures

The results of the X-ray, microscopic, magnetic, and dielectric measurements performed for PFN single crystals and ceramic samples are presented. The influence of the technological conditions on selected physical properties of the materials examined is shown., Comment: 10 pages with 1 table (html), 7 figures (png), submitted to Ferroelectrics

Radiation of dileptons and photons from high energy nuclear collisions provides information on the space-time evolution of the hot dense matter produced therein. We compute this radiation using relativistic hydrodynamics and a coarse-grained version of the microscopic event generator UrQMD, both of which provide a good description of the hadron spectra. The currently most accurate dilepton and photon emission rates from perturbative QCD and from experimentally-based hadronic calculations are used. Comparisons are made to data on central Pb-Pb and Pb-Au collisions taken at the CERN SPS at a beam energy of 158 A GeV. Both hydrodynamics and UrQMD provide very good descriptions of the photon transverse momentum spectrum measured between 1 and 4 GeV, but slightly underestimate the low mass spectrum of e+e- pairs, even with greatly broadened rho and omega vector mesons. Predictions are given for the transverse momentum distribution of dileptons., Comment: 35 pages, 17 figures

No conventional picture of nucleus-nucleus collisions has yet been able to explain the abundance of Omega and AntiOmega in central collisions between Pb nuclei at 158 A GeV at the CERN SPS. We argue that such a deviation from predictions of statistical thermal models and numerical simulations is evidence that they are produced as topological defects in the form of skyrmions arising from the formation of disoriented chiral condensates. The estimated domain size falls in the right range to be consistent with the so far non-observation of DCC from the distribution of neutral pions., Comment: paper presented at the ICPAQGP-2001, Jaipur, India

This paper reviews the field of gamma-ray astronomy and describes future experiments and prospects for advances in fundamental physics and high-energy astrophysics through gamma-ray measurements. We concentrate on recent progress in the understanding of active galaxies, and the use of these sources as probes of intergalactic space. We also describe prospects for future experiments in a number of areas of fundamental physics, including: searches for an annihilation line from neutralino dark matter, understanding the energetics of supermassive black holes, using AGNs as cosmological probes of the primordial radiation fields, constraints on quantum gravity, detection of a new spectral component from GRBs, and the prospects for detecting primordial black holes., Comment: 28 pages, 18 figures, From Snowmass 2001 "The Future of Particle Physics"

No conventional picture of nucleus-nucleus collisions has yet been able to explain the abundance of Omega and AntiOmega hadrons in central collisions between Pb nuclei at 158 A GeV at the CERN SPS. We argue that this is evidence that they are produced from topological defects in the form of skyrmions arising from the formation of disoriented chiral condensates. It is shown that the excess Omega or AntiOmega produced could not be easily washed out in the hadronic medium and therefore would survive in the final state., Comment: 3 embedded PS figures, to appear in J. Phys. G

The statics and dynamics of a surface separating two phases of a relativistic quantum field theory at or near the critical temperature typically make use of a free energy as a functional of an order parameter. This free energy functional also affords an economical description of states away from equilibrium. The similarities and differences between using a scalar field as the order parameter versus the energy density are examined, and a peculiarity is noted. We also point out several conceptual errors in the literature dealing with the dynamical prefactor in the nucleation rate., Comment: 12 pages plus 5 figures

When a potential for a scalar field has two local minima there arise spherical shell-type solutions of the classical field equations due to gravitational attraction. We establish such solutions numerically in a space which is asymptotically de Sitter. It generically arises when the energy scale characterizing the scalar field potential is much less than the Planck scale. It is shown that the mirror image of the shell appears in the other half of the Penrose diagram. The configuration is smooth everywhere with no physical singularity., Comment: 24 pages, 12 figures

We propose a unified space-time picture of baryon stopping and quark-gluon plasma creation in ultrarelativistic heavy-ion collisions. It is assumed that the highly Lorentz contracted nuclei are decelerated by the coherent color field which is formed between them after they pass through each other. This process continues until the field is neutralized by the Schwinger mechanism. Conservation of energy and momentum allow us to calculate the energy losses of the nuclear slabs and the initial energy density of the quark-gluon plasma., Comment: 11 pages in revtex, 2 eps figures

The relativistic viscous fluid equations describing the outflow of high temperature matter created via Hawking radiation from microscopic black holes are solved numerically for a realistic equation of state. We focus on black holes with initial temperatures greater than 100 GeV and lifetimes less than 6 days. The spectra of direct photons and photons from $\pi^0$ decay are calculated for energies greater than 1 GeV. We calculate the diffuse gamma ray spectrum from black holes distributed in our galactic halo. However, the most promising route for their observation is to search for point sources emitting gamma rays of ever-increasing energy., Comment: 33 pages, 13 figures, to be submitted to PRD

The mass shift, width broadening, and spectral density for the $\rho$ and $\omega$ mesons in a heat bath of nucleons and pions are calculated using a general formula which relates the self-energy to the real and imaginary parts of the forward scattering amplitude. We use experimental data to saturate the scattering amplitude at low energies with resonances and include a background Pomeron term, while at high energies a Regge parameterization is used. The real part obtained directly is compared with the result of a dispersion integral over the imaginary part. The peaks of the spectral densities are little shifted from their vacuum positions, but the widths are considerably increased due to collisional broadening. Where possible we compare with the UrQMD model and find quite good agreement. At normal nuclear matter density and a temperature of 150 MeV the spectral density of the $\rho$ meson has a width of 345 MeV, while that for the $\omega$ is in the range 90--150 MeV., Comment: 21 pages revtex + 9 postscript figures

One and two loop self-energies are worked out explicitly for a heavy scalar field interacting weakly with a light self-interacting scalar field at finite temperature. The ring/daisy diagrams and a set of necklace diagrams can be summed simultaneously. This simple model serves to illustrate the connection between multi-loop self-energy diagrams and multiple scattering in a medium., Comment: 15 pages including 7 figures; v2. with appendix added, accepted by PRD

No conventional picture of nucleus-nucleus collisions has yet been able to explain the abundance of Omega and anti-Omega hyperons in central collisions between Pb nuclei at 158 A GeV at the CERN SPS. We argue that this is evidence that they are produced as topological defects arising from the formation of disoriented chiral condensates (DCC) with an average domain size of about 2 fm., Comment: version 2 containing formulas, accepted by PRL

Microscopic black holes explode with their temperature varying inversely as their mass. Such explosions would lead to the highest temperatures in the present universe, all the way to the Planck energy. Whether or not a quasi-stationary shell of matter undergoing radial hydrodynamic expansion surrounds such black holes is been controversial. In this paper relativistic viscous fluid equations are applied to the problem. It is shown that a self-consistent picture emerges of a fluid just marginally kept in local thermal equilibrium; viscosity is a crucial element of the dynamics., Comment: 11 pages, revtex