Search

QTRAJ is a computer code that simulates the propagation of quarkonium in the quark-gluon plasma (QGP) based on the quantum-trajectory algorithm. This algorithm solves a master equation in which the quarkonium is treated as an open quantum system (OQS). A major advantage of this approach is that it turns a 3D spatial evolution for a density matrix into a 1D Schr\"odinger equation for a wavefunction with a non-hermitian Hamiltonian, drastically reducing the computational cost. So far, the interaction implemented in the master equation was obtained within the framework of potential non-relativistic QCD (pNRQCD), and restricted to the regime $rT \ll 1$, where $r$ is the size of the color dipole and $T$ is the temperature. In the environment produced in heavy-ion collisions (HIC's) this limit is accurate for $\Upsilon(1S)$, but the applicability to other quarkonium states is dubious. In the present study we generalize the above approach, extending it to the regime $rT\!\sim\! 1$ in the one-gluon exchange approximation, with proper Hard Thermal Loop (HTL) resummation of medium effects. This is done by implementing new jump operators connecting different color states of the $Q\bar Q$ pair and expanding them in plane waves, giving rise to a variation of the algorithm present in QTRAJ 1.0. Here we provide an overview of this approach comparing the $rT \ll 1$ and $rT\sim 1$ cases, and we discuss prospects for phenomenological application to excited states of bottomonium., Comment: Proceedings of the Conference QNP 2024

We study the potential of $X(3872)$ at finite temperature in the Born-Oppenheimer approximation under the assumption that it is a tetraquark. We argue that, at large number of colors, it is a good approximation to assume that the potential consists in a real part plus a constant imaginary term. The real part is then computed adapting an approach by Rothkopf and Lafferty and using as input lattice QCD determinations of the potential for hybrids. This model allows us to qualitatively estimate at which temperature range the formation of a heavy tetraquark is possible, and to propose a qualitative picture for the dissociation of the state in a medium. Our approach can be applied to other suggested internal structures for the $X(3872)$ and to other exotic states., Comment: 12 pages, 7 figures. Numerical mistakes affecting two plots corrected, text changed accordingly; main conclusions unchanged

The evolution of hard probes in a medium is a complex multiscale problem that significantly benefits from the use of Effective Field Theories (EFTs). Within the EFT framework, we aim to define a series of EFTs in a way that addresses each energy scale individually in separate steps. However, studying hard probes in a medium presents challenges. This is because an EFT is typically constructed by formulating the most general Lagrangian compatible with the problem's symmetries. Nevertheless, medium effects may not always be encoded adequately in an effective action. In this paper, we construct an EFT that is valid for studying the evolution of a heavy quark in a QCD plasma containing few other heavy quarks, where degrees of freedom with an energy of the order of the temperature scale are integrated out. Through this example, we explicitly demonstrate how to handle the doubling of degrees that arise in non-equilibrium field theory. As a result, we derive a Fokker-Planck equation using only symmetry and power counting arguments. The methods introduced in this paper will pave the way for future developments in the study of quarkonium suppression., Comment: 34 pages. Matches published version

We demonstrate the importance of quantum jumps in the nonequilibrium evolution of bottomonium states in the quark-gluon plasma. Based on nonrelativistic effective field theory and the open quantum system framework, we evolve the density matrix of color singlet and octet pairs. We show that quantum regeneration of singlet states from octet configurations is necessary to understand experimental results for the suppression of both bottomonium ground and excited states. The values of the heavy-quarkonium transport coefficients used are consistent with recent lattice QCD determinations., Comment: 5 pages, 7 figures, 1 page supplement; v2 - minor changes, updated references, published version

Using the potential non-relativistic quantum chromodynamics (pNRQCD) effective field theory, we derive a Lindblad equation for the evolution of the heavy-quarkonium reduced density matrix that is accurate to next-to-leading order (NLO) in the ratio of the binding energy of the state to the temperature of the medium. The resulting NLO Lindblad equation can be used to more reliably describe heavy-quarkonium evolution in the quark-gluon plasma at low temperatures compared to the leading-order truncation. For phenomenological application, we numerically solve the resulting NLO Lindblad equation using the quantum trajectories algorithm. To achieve this, we map the solution of the three-dimensional Lindblad equation to the solution of an ensemble of one-dimensional Schr\"odinger evolutions with Monte-Carlo sampled quantum jumps. Averaging over the Monte-Carlo sampled quantum jumps, we obtain the solution to the NLO Lindblad equation without truncation in the angular momentum quantum number of the states considered. We also consider the evolution of the system using only the complex effective Hamiltonian without stochastic jumps and find that this provides a reliable approximation for the ground state survival probability at LO and NLO. Finally, we make comparisons with our prior leading-order pNRQCD results and experimental data available from the ATLAS, ALICE, and CMS collaborations., Comment: 40 pages, 8 figures

Quarkonium suppression is one of the more useful observables to obtain information about the hot medium created in ultrarelativistic heavy-ion collisions. In this manuscript, we discuss a simple way to implement both the initial-state effects and the hot-medium evolution, and to compute the quarkonium nuclear modification factor if the survival probability for a bound state at a given energy density is known. Based on the Glauber model, the temperature of the evolving medium and the centrality dependence of the nuclear modification factor will be analysed., Comment: 22 pages, 4 figure. Updated to match published version

We report predictions for the suppression and elliptic flow of the $\Upsilon(1S)$, $\Upsilon(2S)$, and $\Upsilon(3S)$ as a function of centrality and transverse momentum in ultra-relativistic heavy-ion collisions. We obtain our predictions by numerically solving a Lindblad equation for the evolution of the heavy-quarkonium reduced density matrix derived using potential nonrelativistic QCD and the formalism of open quantum systems. To numerically solve the Lindblad equation, we make use of a stochastic unraveling called the quantum trajectories algorithm. This unraveling allows us to solve the Lindblad evolution equation efficiently on large lattices with no angular momentum cutoff. The resulting evolution describes the full 3D quantum and non-abelian evolution of the reduced density matrix for bottomonium states. We expand upon our previous work by treating differential observables and elliptic flow; this is made possible by a newly implemented Monte-Carlo sampling of physical trajectories. Our final results are compared to experimental data collected in $\sqrt{s_{NN}} = 5.02$ TeV Pb-Pb collisions by the ALICE, ATLAS, and CMS collaborations., Comment: 14 pages, 11 figures

We introduce an open-source package called QTraj that solves the Lindblad equation for heavy-quarkonium dynamics using the quantum trajectories algorithm. The package allows users to simulate the suppression of heavy-quarkonium states using externally-supplied input from 3+1D hydrodynamics simulations. The code uses a split-step pseudo-spectral method for updating the wave-function between jumps, which is implemented using the open-source multi-threaded FFTW3 package. This allows one to have manifestly unitary evolution when using real-valued potentials. In this paper, we provide detailed documentation of QTraj 1.0, installation instructions, and present various tests and benchmarks of the code., Comment: 35 pages, 10 figures; v2 - minor typos fixed; published version

The study of heavy quarkonium suppression in heavy-ion collisions represents an important source of information about the properties of the quark-gluon plasma produced in such collisions. In a previous paper, we have considered how to model the evolution of a quarkonium in such a way that the solution of the resulting equations evolves toward the correct thermal equilibrium distribution for an homogeneous and static medium. We found that it is crucial to take into account the energy gap between singlet and octet configurations when the temperature is not much greater than this energy gap. In this manuscript, we explore in more detail the phenomenological consequences of this observation in the more realistic situation of an expanding system. We consider two different scenarios, based on the same approximation scheme, but on different choices of parameters. In the first case, we rely on a Hard Thermal Loop approximation, while the second case is based on a recent determination of the static potential in lattice QCD. In both cases, we compute the decay width and the nuclear modification factor, both taking the energy gap into account and ignoring it. We find that the impact on the predictions is as large in the expanding medium as it is in the static case. Our conclusion is that this energy gap should be taken into account in phenomenological studies., Comment: Published version. 25 pages, 11 figures

We solve the Lindblad equation describing the Brownian motion of a Coulombic heavy quark-antiquark pair in a strongly coupled quark-gluon plasma using the highly efficient Monte Carlo wave-function method. The Lindblad equation has been derived in the framework of pNRQCD and fully accounts for the quantum and non-Abelian nature of the system. The hydrodynamics of the plasma is realistically implemented through a 3+1D dissipative hydrodynamics code. We compute the bottomonium nuclear modification factor and compare with the most recent LHC data. The computation does not rely on any free parameter, as it depends on two transport coefficients that have been evaluated independently in lattice QCD. Our final results, which include late-time feed down of excited states, agree well with the available data from LHC 5.02 TeV PbPb collisions., Comment: 47 pages, 20 figures, corrections to feed down matrix leading to slight change in phenomenological predictions, journal version

We study the transitions between the different color states of a static quark-antiquark pair, singlet and octet, in a thermal medium. This is done non-perturbatively exploiting the infinite mass limit of QCD. This study is interesting because it can be used for future developments within the framework of Effective Field Theories (EFTs) and because it can be combined with other techniques, like lattice QCD or AdS/CFT, to gain non-perturbative information about the evolution of quarkonium in a medium. We also study the obtained expressions in the large $N_{c}$ limit. This allows us to learn lessons that are useful to simplify phenomenological models of quarkonium in a plasma., Comment: 32 pages, 14 figures. Matches published version

We study exclusive quarkonium production in the dipole picture at next-to-leading order (NLO) accuracy, using the non-relativistic expansion for the quarkonium wavefunction. This process offers one of the best ways to obtain information about gluon distributions at small $x$, in ultraperipheral heavy ion collisions and in deep inelastic scattering. The quarkonium light cone wave functions needed in the dipole picture have typically been available only at tree level, either in phenomenological models or in the nonrelativistic limit. In this paper, we discuss the compatibility of the dipole approach and the non-relativistic expansion and compute NLO relativistic corrections to the quarkonium light-cone wave function in light-cone gauge. Using these corrections we recover results for the NLO decay width of quarkonium to $e^{+}e^{-}$ and we check that the non-relativistic expansion is consistent with ERBL evolution and with B-JIMWLK evolution of the target. The results presented here will allow computing the exclusive quarkonium production rate at NLO once the one loop photon wave function with massive quarks, currently under investigation, is known., Comment: 33 pages. 9 figures. Updated to incorporate the changes made to the published version

Hard Probes are an essential tool to discover the properties of the quark-gluon plasma created in heavy-ion collisions. The study of hard probes always involves taking into account very different energy scales, and this is precisely the situation in which Effective Fields Theories (EFTs) are useful. EFTs can be used to separate the short-distance and perturbative physics from the long-distance and non-perturbative. This method combined with Lattice QCD evaluations of the long-distance effects can provide accurate and first principles results. In this proceeding, I will report recent advances in this direction. Results from an EFT computation of quarkonium $R_{AA}$ at $\sqrt{s_{NN}}=5.02\,\textit{TeV}$ are shown for the first time here., Comment: 8 pages, 2 figures, for the International Conference on Hard and Electromagnetic Probes of High-Energy Nuclear Collisions, 30 September - 5 October 2018 in Aix-Les-Bains, Savoie, France

We derive equations of motion for the reduced density matrix of a heavy quarkonium in contact with a quark-gluon plasma in thermal equilibrium. These equations allow in particular a proper treatment of the regime when the temperature of the plasma is comparable to the binding energy of the quarkonium. These equations are used to study how the quarkonium approaches equilibrium with the plasma, and we discuss the corresponding entropy increase, or free energy decrease, depending on the temperature regime. The effect of collisions can be accounted for by the generalization of the imaginary potential introduced in previous studies, and from which collision rates are derived. An important outcome of the present analysis is that this imaginary potential has a sizeable dependence on the energy of the relevant transitions., Comment: 47 pages, 8 figures

We derive equations for the time evolution of the reduced density matrix of a collection of heavy quarks and antiquarks immersed in a quark gluon plasma. These equations, in their original form, rely on two approximations: the weak coupling between the heavy quarks and the plasma, the fast response of the plasma to the perturbation caused by the heavy quarks. An additional semi-classical approximation is performed. This allows us to recover results previously obtained for the abelian plasma using the influence functional formalism. In the case of QCD, specific features of the color dynamics make the implementation of the semi-classical approximation more involved. We explore two approximate strategies to solve numerically the resulting equations in the case of a quark-antiquark pair. One involves Langevin equations with additional random color forces, the other treats the transition between the singlet and octet color configurations as collisions in a Boltzmann equation which can be solved with Monte Carlo techniques., Comment: 62 pages, 14 figures

In the early stages of heavy-ion collisions, the hot QCD matter expands more longitudinally than transversely. This imbalance causes the system to become rapidly colder in the longitudinal direction and a local momentum anisotropy appears. In this paper, we study the heavy-quarkonium spectrum in the presence of a small plasma anisotropy. We work in the framework of pNRQCD at finite temperature. We inspect arrangements of non-relativistic and thermal scales complementary to those considered in the literature. In particular, we consider temperatures larger and Debye masses smaller than the binding energy, which is a temperature range relevant for presently running LHC experiments. In this setting we compute the leading thermal corrections to the binding energy and the thermal width induced by quarkonium gluo-dissociation., Comment: 6 pages, 5 figures

In a model where Majorana neutrinos heavier than the electroweak scale couple to Standard Model Higgs bosons and leptons, we compute systematically thermal corrections to the direct and indirect CP asymmetries in the Majorana neutrino decays. These are key ingredients entering the equations that describe the thermodynamic evolution of the induced lepton-number asymmetry eventually leading to the baryon asymmetry in the universe. We compute the thermal corrections in an effective field theory framework that assumes the temperature smaller than the masses of the Majorana neutrinos and larger than the electroweak scale, and we provide the leading corrections in an expansion of the temperature over the mass. In this work, we consider the case of two Majorana neutrinos with nearly degenerate masses., Comment: Equations (B.25), (B.26) and (B.27) corrected in the parts depending on the SU(2) gauge coupling; final results modified accordingly

Heavy quarkonium suppression was proposed long ago as a signal of the formation of a deconfined phase in heavy ion collisions. Originally the mechanism responsible for this suppression was thought to be color screening. However perturbative computations and recent lattice studies suggest the existence of an imaginary part of the potential which could have a more important role for suppression than screening. In this work we review some general aspects of effective field theories for heavy quarkonium in the medium and discuss the physical phenomena behind the imaginary part of the potential and the decay width of heavy quarkonium and their corresponding cross sections., Comment: 5 pages, no figures; Contribution for FAIRNESS 2013 proceedings

Heavy Majorana neutrinos enter in many scenarios of physics beyond the Standard Model: in the original seesaw mechanism they provide a natural explanation for the small masses of the Standard Model neutrinos and in the simplest leptogenesis framework they are at the origin of the baryonic matter of the universe. In this paper, we develop an effective field theory for non-relativistic Majorana particles, which is analogous to the heavy-quark effective theory. Then, we apply it to the case of a heavy Majorana neutrino decaying in a hot and dense plasma of Standard Model particles, whose temperature is much smaller than the mass of the Majorana neutrino but still much larger than the electroweak scale. The neutrino width gets zero-temperature contributions that can be computed from in-vacuum matrix elements, and thermal corrections. Only the latter will be addressed. Symmetry and power counting arguments made manifest by the effective field theory restrict the form of the thermal corrections and simplify their calculation. The final result agrees with recent determinations obtained with different methods. The effective field theory presented here is suitable to be used for a variety of different models involving non-relativistic Majorana fermions., Comment: Published version. 25 pages, 7 figures, references added

By means of effective field theory techniques, we study the modifications of some properties of weakly coupled heavy quarkonium states propagating through a quark-gluon plasma at temperatures much smaller than the heavy quark mass, m_Q. Two different cases are considered, corresponding to two different hierarchies between the typical size of the bound state, r, the binding energy, E, the temperature, T, and the screening mass, m_D. The first case corresponds to the hierarchy m_Q >> 1/r >> T >> E >> m_D, relevant for moderate temperatures, and the second one to the hierarchy m_Q >> T >> 1/r, m_D >> E, relevant for studying the dissociation mechanism. In the first case we determine the perturbative correction to the binding energy and to the decay width of states with arbitrary angular momentum, finding that the width is a decreasing function of the velocity. A different behavior characterizes the second kinematical case, being the width of s-wave states a non-monotonic function of the velocity, increasing at moderate velocities and decreasing in the ultra-relativistic limit. We obtain a simple analytical expression of the decay width for T >> 1/r >> m_D >> E at moderate velocities, and we derive the s-wave spectral function for the more general case T >> 1/r, m_D >> E. A brief discussion of the possible experimental signatures as well as a comparison with the relevant lattice data are also presented., Comment: Published version, title changed and improved bibliography, 21 pages. 8 figures

In a weak-coupling effective field theory framework we study quarkonium dissociation induced by inelastic scattering with partons in the medium. This is the dominant dissociation process for temperatures such that the Debye mass is larger than the binding energy. We evaluate the dissociation cross section and the corresponding thermal decay width. At leading order we derive a convolution formula relating the two, which is consistent with the optical theorem and QCD at finite temperature. Bound state effects are systematically included. They add contributions to the cross section and width that are beyond a quasi-free approximation, whose validity is critically reviewed. For temperatures such that the Debye mass is smaller than the binding energy, the dominant dissociation mechanism is gluo-dissociation consisting in quarkonium dissociation induced by the absorbtion of a gluon from the medium. We calculate the gluo-dissociation cross section and width at next-to-leading-order accuracy., Comment: 33 pages. 11 figures

In the framework of Soft-Collinear Effective Theory, the jet quenching parameter, $\hat{q}$, has been evaluated by adding the effect of Glauber gluon interactions to the propagation of a highly-energetic collinear parton in a medium. The result, which holds in covariant gauges, has been expressed in terms of the expectation value of two Wilson lines stretching along the direction of the four-momentum of the parton. In this paper, we show how that expression can be generalized to an arbitrary gauge by the addition of transverse Wilson lines. The transverse Wilson lines are explicitly computed by resumming interactions of the parton with Glauber gluons that appear only in non-covariant gauges. As an application of our result, we discuss the contribution to $\hat{q}$ coming from transverse momenta of order $g^2T$ in a medium that is a weakly-coupled quark-gluon plasma., Comment: 31 pages, 7 figures; journal version

The thermal width of heavy-quarkonium bound states in a quark-gluon plasma has been recently derived in an effective field theory approach. Two phenomena contribute to the width: the Landau damping phenomenon and the break-up of a colour-singlet bound state into a colour-octet heavy quark-antiquark pair by absorption of a thermal gluon. In the paper, we investigate the relation between the singlet-to-octet thermal break-up and the so-called gluo-dissociation, a mechanism for quarkonium dissociation widely used in phenomenological approaches. The gluo-dissociation thermal width is obtained by convoluting the gluon thermal distribution with the cross section of a gluon and a 1S quarkonium state to a colour octet quark-antiquark state in vacuum, a cross section that at leading order, but neglecting colour-octet effects, was computed long ago by Bhanot and Peskin. We will, first, show that the effective field theory framework provides a natural derivation of the gluo-dissociation factorization formula at leading order, which is, indeed, the singlet-to-octet thermal break-up expression. Second, the singlet-to-octet thermal break-up expression will allow us to improve the Bhanot--Peskin cross section by including the contribution of the octet potential, which amounts to include final-state interactions between the heavy quark and antiquark. Finally, we will quantify the effects due to final-state interactions on the gluo-dissociation cross section and on the quarkonium thermal width., Comment: 17 pages, 6 figures

Heavy quarkonium at finite temperature has been the subject of intense theoretical studies, for it provides a potentially clean probe of the quark-gluon plasma. Recent studies have made use of effective field theories to exploit in a systematic manner the hierarchy of energy scales that characterize the system. In the case of a quarkonium in a medium whose temperature is smaller than the typical momentum transfer in the bound state but larger than its energy, the suitable effective field theory is pNRQCD_HTL, where degrees of freedom with energy or momentum larger than the binding energy have been integrated out. Thermal effects are expected to break Poincar\'e invariance, which, at zero temperature, manifests itself in a set of exact relations between the matching coefficients of the effective field theory. In the paper, we evaluate the leading-order thermal corrections to the spin-orbit potentials of pNRQCD_HTL and show that Poincar\'e invariance is indeed violated., Comment: 17 page, 4 figures. Version published on JHEP

We study the propagation of non-relativistic bound states moving at constant velocity across a homogeneous thermal bath and we develop the effective field theory which is relevant in various dynamical regimes. We consider values of the velocity of the bound state ranging from moderate to highly relativistic and temperatures at all relevant scales smaller than the mass of the particles that form the bound state. In particular, we consider two distinct temperature regimes, corresponding to temperatures smaller or higher than the typical momentum transfer in the bound state. For temperatures smaller or of the order of the typical momentum transfer, we restrict our analysis to the simplest system, a hydrogen-like atom. We build the effective theory for this system first considering moderate values of the velocity and then the relativistic case. For large values of the velocity of the bound state, the separation of scales is such that the corresponding effective theory resembles the soft collinear effective theory (SCET). For temperatures larger than the typical momentum transfer we also consider muonic hydrogen propagating in a plasma which contains photons and massless electrons and positrons, so that the system resembles very much heavy quarkonium in a thermal medium of deconfined quarks and gluons. We study the behavior of the real and imaginary part of the static two-body potential, for various velocities of the bound state, in the hard thermal loop approximation. We find that Landau damping ceases to be the relevant mechanism for dissociation from a certain "critical" velocity on in favor of screening. Our results are relevant for understanding how the properties of heavy quarkonia states produced in the initial fusion of partons in the relativistic collision of heavy ions are affected by the presence of an equilibrated quark-gluon plasma., Comment: 25 pages, 4 figures. Minor corrections, reference added. Closely matches version to be published in Physical Review D

We illustrate how to apply modern effective field theory techniques and dimensional regularization to factorise the various scales which appear in QED bound states at finite temperature. We focus here on the muonic hydrogen atom. Vacuum polarization effects make the physics of this atom at finite temperature very close to that of heavy quarkonium states. We comment on the implications of our results for these states in the quark gluon plasma. In particular, we estimate the effects of a finite charm quark mass in the dissociation temperature of bottomonium., Comment: 22 pages, 8 figures. Journal version, reference added

We calculate the heavy quarkonium energy levels and decay widths in a quark-gluon plasma, whose temperature T and screening mass m_D satisfy the hierarchy m alpha_s >> T >> m alpha_s^2 >> m_D (m being the heavy-quark mass), at order m alpha_s^5. We first sequentially integrate out the scales m, m alpha_s and T, and, next, we carry out the calculations in the resulting effective theory using techniques of integration by regions. A collinear region is identified, which contributes at this order. We also discuss the implications of our results concerning heavy quarkonium suppression in heavy ion collisions., Comment: 25 pages, 2 figures

We study the recently proposed effective field theory for the phonon of an arbitrary non-relativistic superfluid. After computing the one-loop phonon self-energy, we obtain the low temperature T contributions to the phonon dispersion law at low momentum, and see that the real part of those can be parametrized as a thermal correction to the phonon velocity. Because the phonons are the quanta of the sound waves, at low momentum their velocity should agree with the speed of sound. We find that our results match at order T^4ln(T) with those predicted by Andreev and Khalatnikov for the speed of sound, derived from the superfluid hydrodynamical equations and the phonon kinetic theory. We get also higher order corrections of order T^4, which are not reproduced pushing naively the kinetic theory computation. Finally, as an application, we consider the cold Fermi gas in the unitarity limit, and find a universal expression for the low T relative correction to the speed of sound for these systems., Comment: 14 pages, 2 figures. References added

We compute the coefficients of bulk viscosity for a non-relativistic superfluid corresponding to a fermionic system close to the unitarity limit. We consider the low temperature regime assuming that the transport properties of the system are dominated by phonons. To compute the coefficients of bulk viscosity we use kinetic theory in the relaxation time approximation and the low energy effective field theory of the corresponding system. We show that the three independent bulk viscosity coefficients, $\zeta_1, \zeta_2, \zeta_3$, associated with irreversible flows vanish for phonons with a linear dispersion law. Considering a phonon dispersion law with a cubic term in momentum we find that in the conformal limit $\zeta_1 = \zeta_2=0$, while $\zeta_3$ is non-zero. Including a conformal breaking term which arises for a large but finite s-wave scattering length, $a$, at the leading order in $1/a$ we obtain that $\zeta_1 \propto 1/a$ and $\zeta_2 \propto 1/a^2$., Comment: 14 pages, minor corrections, references added, closely matches published version

We illustrate how to apply modern effective field theory techniques and dimensional regularization to factorize the various scales which appear in non-relativistic bound states at finite temperature. We focus here on the simplest case: the hydrogen atom. We discuss in detail the interplay of the hard, soft and ultrasoft scales of the non-relativistic system at zero temperature with the additional scales induced at finite temperature. We also comment on the implications of our results for heavy quarkonium bound states in the quark gluon plasma., Comment: 28 pages, 1 table. Numerical factors in formulas (23), (25), (26) and (28), as well as a few misprints elswhere corrected. Supersedes journal version