the quark structure of the nucleons , suggested by quantum chromodynamics , opens the possibility of a hadron - quark phase transition at high densities and/or temperatures , as suggested by @xcite . in the theories of strong interaction
quark bag models suppose that breaking of physical vacuum takes place inside hadrons .
if the hypothesis of the quark matter is true , then some compact objects identified with neutron stars could actually be strange stars , built entirely of strange matter @xcite . for a review of strange star properties , see @xcite .
most of the investigations of quark star properties have been done within the framework of the so - called mit bag model .
assuming that interactions of quarks and gluons are sufficiently small , neglecting quark masses and supposing that quarks are confined to the bag volume ( in the case of a bare strange star , the boundary of the bag coincides with stellar surface ) , the energy density @xmath0 and pressure @xmath1 of a quark - gluon plasma are related , in the mit bag model , by the equation of state ( eos ) @xcite @xmath2 where @xmath3 is the difference between the energy density of the perturbative and non - perturbative qcd vacuum ( the bag constant ) .
equation ( [ 1 ] ) is essentially the equation of state of a gas of massless particles with corrections due to the qcd trace anomaly and perturbative interactions .
more sophisticated investigations of quark - gluon interactions have shown that eq .
( [ 1 ] ) represents a limiting case of more general equations of state . for example mit bag models with massive strange quarks and lowest order qcd interactions lead to some corrections terms in the equation of state of quark matter .
models incorporating restoration of chiral quark masses at high densities and giving absolutely stable strange matter can no longer be accurately described by using eq .
( [ 1 ] ) . on the other hand in models in which quark interaction is described by an interquark potential originating from gluon exchange and by a density dependent scalar potential which restores the chiral symmetry at high densities @xcite
, the equation of state @xmath4 can be well approximated by a linear function in the energy density @xmath5 @xcite .
@xcite has studied the linear approximation of the equation of state , obtaining all parameters of the eos as polynomial functions of strange quark mass , qcd coupling constant and bag constant .
a complete description of static strange stars has been obtained based on numerical integration of mass continuity and tov ( hydrostatic equilibrium ) equations for different values of the bag constant @xcite . using numerical methods the maximum gravitational mass @xmath6 and the maximum radius @xmath7 of the strange star ,
have been obtained , as a function of the bag constant , in the form @xcite @xmath8 where @xmath9 mev @xmath10 .
there are several proposed mechanisms for the formation of quark stars .
quark stars are expected to form during the collapse of the core of a massive star after the supernova explosion as a result of a first or second order phase transition , resulting in deconfined quark matter @xcite .
the proto - neutron star core or the neutron star core is a favorable environment for the conversion of ordinary matter to strange quark matter @xcite .
another possibility is that some neutron stars in low - mass x - ray binaries can accrete sufficient mass to undergo a phase transition to become strange stars @xcite .
this mechanism has also been proposed as a source of radiation emission for cosmological @xmath11-ray bursts @xcite .
strange quark matter is filled with electromagnetic waves in thermodynamic equilibrium with quarks .
photon emissivity is the basic parameter for determining macroscopic properties of objects made of strange matter . @xcite
have shown that , because of very high plasma frequency near the strange matter edge , photon emissivity of strange matter is very low . for temperatures @xmath12 , where @xmath13 mev is the characteristic transverse plasmon cuttof energy ,
the equilibrium photon emissivity of strange matter is negligible small , compared to the black body one .
the dispersion relation of the electromagnetic waves in quark matter can be written as @xmath14 , where @xmath15 is the frequency of the waves , @xmath16 is their wave number , @xmath17 is the plasma frequency , with @xmath18 the baryon number density of strange quark matter and @xmath19 is the chemical potential of the quarks ( usov 1998 ) .
propagating modes exist only for @xmath20 .
hence the lower limit of the energy @xmath21 of the electromagnetic quanta , which are in equilibrium with quarks is given by @xmath22 mev , where @xmath23 @xmath24 is normal nuclear matter density .
at the quark star surface @xmath25 and the spectrum of equilibrium photons is very hard , with @xmath26 mev @xcite .
the bremsstrahlung emissivity of quark matter has been estimated by @xcite .
they have found that the surface radiation is about four orders of magnitude weaker than the equilibrium black body radiation , even they both have the same temperature dependence .
the structure of a realistic strange star is very complicated , but its basic properties can be described as follows @xcite .
beta - equilibrated strange quark - star matter consists of an approximately equal mixture of up @xmath27 , down @xmath28 and strange @xmath29 quarks , with a slight deficit of the latter .
the fermi gas of @xmath30 quarks constitutes a single color - singlet baryon , with baryon number @xmath31 .
this structure of the quarks leads to a net positive charge inside the star .
since stars in their lowest energy state are supposed to be charge neutral , electrons must balance the net positive quark charge in strange matter stars .
the electrons , being bounded to the quark matter by the electromagnetic interaction and not by the strong force , are able to move freely across the quark surface , but clearly can not move to infinity because of the electrostatic attraction of quarks .
the electron distribution extends up to @xmath32 fm above the quark surface .
the coulomb barrier at the quark surface of a hot strange star may also be a powerful source of @xmath33 pairs , which are created in the extremely strong electric field of the barrier . at surface temperatures of around @xmath34
k , the luminosity of the outflowing plasma may be of the order @xmath35 ergs@xmath36 @xcite .
moreover , as shown by @xcite , the thermal luminosity from the star surface , due to both photon emission and @xmath33 pair production may be , for about one day for normal quark matter and for up to a hundred years for superconducting quark matter , orders of magnitude higher than the eddington limit .
it is the purpose of the present paper to reconsider the problem of the photon emissivity , via bremsstrahlung radiation , of the quark stars .
equilibrium radiation , transmitted through the surface , is dominant at temperatures @xmath37 k , while below this temperature the bremsstrahlung radiation from the surface layer prevails .
hence bremsstrahlung is the main source of radiation for cold quark stars .
this electromagnetic radiation is generated in the dense medium at the stellar surface in which extremely relativistic quarks move . in the 1950 s @xcite
predicted that the cross section for bremsstrahlung from highly relativistic particles in dense media is suppressed , due to the interference between amplitudes of nearby interactions .
the suppression has its roots in the uncertainty principle .
the kinematics of the bremsstrahlung requires that the longitudinal momentum transfer between the fixed and the scattered particles must be small . on the other hand
the uncertainty principle requires that the interaction must occur over a large longitudinal distance scale ( formation zone ) . if the charged particle coulomb scatters while traversing this zone , the bremsstrahlung amplitude from before and after the scattering can interfere , reducing the amplitude for photon emission @xcite .
this effect has been recently confirmed experimentally for the bremsstrahlung emission of electrons in different materials at a 5% accuracy level @xcite .
since the surface of the quark star consists of a very dense medium , with density of the order @xmath38 , multiple collisions of the electromagnetically radiating particle lead to a significant suppression of the quark - quark bremsstrahlung . by adopting a simple model for the elastic scattering of quarks , we derive the spectrum of the classical bremsstrahlung radiation and the emissivity of the quark matter , by considering the multiple scattering effects in the electromagnetic radiation of quark matter . for the range of densities existing at the surfaces of strange stars , which are higher than the nuclear density
, the landau - pomeranchuk effect could reduce the bremsstrahlung emissivity of quark matter by an order of magnitude .
photon emissivity of quark - gluon plasma ( qgp ) , which is conjectured to be formed in ultrarelativistic heavy ion collisions , has been extensively investigated recently ( for a recent review of direct photon emission from qgp , including comparisons of theoretical predictions with experiments see peitzmann & thoma ( 2002 ) ) .
photons and dilepton pairs only interact electromagnetically and their mean free paths are much larger than the size of the qgp .
hence these electromagnetic probes leave the hot and dense qgp without further scattering , thus providing important informations about the early stages of the collision and the structure of the qgp . as sources of direct photons
one can consider quark annihilation , compton scattering and bremsstrahlung following the initial hard scattering of partons of the nuclei , as well as thermal photons from the qgp and from hadronic interactions in the hot hadronic gas after the hadronization of the plasma ( peitzmann & thoma 2002 ) .
the first observation of direct photon production in ultrarelativistic heavy ion collisions has been reported by the wa98 collaboration in @xmath39 collisions at @xmath40 gev at the super proton synchrotron at cern ( aggarwal et al . 2000 ) .
the results display a clear excess of direct photons above the expected background from hadronic decays in the range of transverse momentum @xmath41 gev@xmath42 in the most central collisions .
these experimental findings provide a confirmation of the feasibility of direct photons as reliable probes in heavy ion collisions and may pave the way for the understanding of the formation and evolution of the qgp .
the evaluation of the photo - emission rate from the qgp was initiated by kapusta et al .
( 1991 ) and by baier et al .
. however , as shown by aurenche et al .
( 2000 ) , these initial estimations were incomplete , since the bremsstrahlung and inelastic pair annihilation processes contains collinear enhancements which cause them to contribute at the same parametric order in the coupling as the two to two processes even for large photon energies .
but these calculations are also incomplete , as they did not incorporate the suppression of photon emission due to multiple scatterings during the photon emission process , which limits the coherence length of the emitted radiation ( the landau - pomeranchuk - migdal effect ) . the rate of photo - production by bremsstrahlung and inelastic pair annihilation in a hot , equilibrated plasma at zero chemical potential , fully including the lpm effect , was calculated , to leading orders in both the electromagnetic and strong coupling constants , by ( arnold , moore & jaffe 2001 , arnold , moore & jaffe 2002 ) .
the emission of hard photons from the qgp plasma , using a model based on the thermodynamics of qcd , and the determination of the initial temperature of the expanding fireball has been considered recently by renk ( 2003 ) .
the radiation emitted by the quark layer at the surface of the quark star must also cross the thin electron layer extending beyond the star s surface .
we also consider the effect of this layer on the outgoing radiation and show that it can significantly reduce the intensity of the radiation coming from the quark star .
the present paper is organized as follows . in section
ii we discuss the landau - pomeranchuk - migdal effect for quark matter .
the bremsstrahlung emissivity of quark matter is considered in section iii .
the effect of the electron layer on the radiation outgoing from the quark star surface is analyzed in section iv . in section
v we discuss and conclude our results . throughout this paper we use units so that @xmath43 , with @xmath44 the boltzmann constant .
the dispersion relation shows that only photons with frequency @xmath45 can propagate inside strange matter .
propagation of electromagnetic waves of frequencies lower than @xmath46 is exponentially damped , with the damping coefficient @xmath47 depending on @xmath15 and @xmath46 .
for @xmath48 the value of @xmath47 is given by @xmath49 @xcite .
this means that photons with @xmath50 emitted in various processes inside strange matter are absorbed .
the mean free path of such a photon is of the order of @xmath51 fm @xcite .
therefore , only photons produced just below the surface with momenta pointing outwards can leave strange matter .    in the following
we shall consider the photon emissivity of quark matter , by also taking into account the effect of multiple scatterings on the emitted radiation .
this effect is called the landau - pomeranchuk - migdal ( lpm ) effect , and was considered initially by @xcite , who used classical electrodynamics to demonstrate bremsstrahlung suppression due to multiple scattering : the suppression comes from the interference between photons emitted by different elements of the charged particle trajectory . to study the astrophysical implications of this effect on quark matter we shall take into account the bremsstrahlung photons produced in the quark - quark interactions , that is in the process @xmath52 .
the lpm effect can be discussed in a simple qualitative way as follows .
let s consider a bremsstrahlung photon of energy @xmath15 produced by an ultrarelativistic quark of energy @xmath53 and mass @xmath54 , where @xmath55 .
in this kinematic regime the average angle @xmath56 between the incident quark and the produced photon is small , @xmath57 .
we assume that the angle between the incident and scattered quark is also small .
moreover , we shall assume that all quarks have the same density dependent effective mass @xmath58 , where @xmath59 is the qcd coupling constant and @xmath60 , with @xmath61 the quark particle number densities .
neglecting the scattering and quark - photon angles , it follows that the longitudinal momentum transfer @xmath62 to the static quark is @xmath63 , where @xmath64 @xcite .
the uncertainty principle then requires that the spatial position of the bremsstrahlung process has a longitudinal uncertainty of @xmath65 . in alternate language ,
the quark and photon slowly split apart over the distance @xmath66 . in a sufficiently dense medium
the quark mean free path is much less than @xmath66 ; in the lpm effect the relevant interaction is multiple scattering .
bremsstrahlung is suppressed when the mean square multiple scattering angle over the distance @xmath66 , @xmath67 @xmath68 , where @xmath69 and @xmath70 ^{-1}$ ] , with @xmath71 , is the radiation length , is greater than or equal to @xmath72 , @xmath73 @xcite .
therefore the bremsstrahlung differential cross section for the production of a photon is suppressed when @xmath74 , where @xmath75 @xcite . by assuming that the temperature of the star @xmath76 , with @xmath77 the fermi momentum , the fermi distribution factors force all the quarks to have @xmath78 . for an extreme relativistic quark gas @xmath79 , where @xmath80 is the statistical weight for quark matter .
therefore the critical lpm frequency for bremsstrahlung photon emission in quark matter is given by @xmath81 where @xmath82 .
for strange quark matter @xmath83 .
emission of the bremsstrahlung radiation with frequency smaller than @xmath84 is suppressed , due to the effect of multiple scattering .
as expected , the lpm frequency increases with the density of the strange matter .
the variation of @xmath85 as a function of the particle number density is represented , for different values of the strong coupling constant @xmath59 , in fig .  1 .
the most important parameter characterizing the lpm effect is the average time between two collisions , denoted @xmath86 .
it can be calculated from @xmath87 , where @xmath88 is the density of the medium which depends on its composition and on temperature @xcite .
@xmath89 denotes the cross section of a test particle having velocity @xmath90 and interacting with the particles of the medium .
thus generally @xmath86 is a function of the temperature and density of the medium and of the momentum of the test particle .
the value of @xmath86 is essential for the importance of the lpm effect : if @xmath91 is zero ( @xmath92 ) , the effect is absent @xcite .    the starting point for the rigorous implementation of the lpm effect
is the following textbook equation for the energy radiated per unit momentum by the charged particle @xcite , @xmath93 } \vec{n}\times \vec{v}_{l}dt\right| ^{2 } ,   \label{1a}\ ] ] where @xmath94 describes the trajectory of the charged particle , @xmath95 is its velocity , @xmath96 , @xmath97 is the fine structure constant and @xmath98 is the charge of @xmath99 quarks .
if only one collision occurs and if the velocity is a constant ( @xmath100 ) before and after the collision ( @xmath101 ) , then the radiated energy is given by @xcite @xmath102 which is the standard textbook form for the calculation of the intensity .
if , however , a series of collisions occurs , one has to change eqs .
( [ 1a ] ) and ( [ 2a ] ) to a sum over all pieces of the trajectory @xcite , @xmath103 } \vec{n}\times \vec{v}_{1l}dt+\int_{t_{1}}^{t_{2}}e^{i\left [ \omega t-\vec{k}\cdot \vec{r}_{l}\left ( t\right ) \right ] } \vec{n}\times \vec{v}_{2l}dt+ ... +\int_{t_{n-1}}^{t_{n}}e^{i\left [ \omega t-\vec{k}\cdot \vec{r}_{l}\left ( t\right ) \right ] } \vec{n}\times \vec{v}_{nl}dt+ .... \right| ^{2}\end{aligned}\ ] ]    assuming that the velocity is a constant between two collisions leads to @xmath104 where @xmath105 is the total number of collisions .
the calculation of the square gives @xmath106 +   \nonumber \\ & & 2\frac{\alpha } { \left ( 2\pi \right ) ^{2}}\sum_{m = u ,
d , s}\left ( \frac{q_{m}}{e}\right ) ^{2}\textrm{re}\sum_{j > l}^{n}\frac{\vec{v}_{jm}\cdot \vec{v}_{lm}-\left ( \vec{n}\cdot \vec{v}_{jm}\right ) \left ( \vec{n}\cdot \vec{v}_{lm}\right ) } { \left ( \omega -\vec{k}\cdot \vec{v}_{jm}\right ) \left ( \omega -\vec{k}\cdot \vec{v}_{lm}\right ) } \exp \left [ i\sum_{i = l}^{j-1}\xi _ { i}\left ( \omega -\vec{k}\cdot \vec{v}_{im}\right ) \right]\nonumber \\ & & \left ( 1-e^{-i\left ( \omega -\vec{k}\cdot \vec{v}_{jm}\right ) \xi _ { j}}\right ) \times   \left ( 1-e^{-i\left ( \omega -\vec{k}\cdot \vec{v}_{lm}\right ) \xi _ { l}}\right ) , \end{aligned}\ ] ] where we denoted @xmath107 .
as the quark gas is taken to be in thermal equilibrium , many possible velocities can result in between two collisions .
for this reason the non - diagonal terms in the above equation will give zero contribution after one averages over all velocities , providing of course that no correlation exists between the velocity before and after the collision .
only the diagonal terms remain in this case .
the time between two successive collisions is given by @xmath108 .
taking an average over the time between two collisions can be obtained using the following distribution @xcite @xmath109    then the electromagnetic energy radiated away by the quarks is given by @xmath110 d\xi , \ ] ] where we have introduced the average velocities by means of the definitions @xmath111 , @xmath112 and we also denoted @xmath113 .
the averaging over the time between two collisions can be done with the use of the mathematical identity @xmath114 .
therefore we obtain for the intensity of the radiation emitted by a quark moving in a dense medium the expression @xmath115    since @xmath116 , the intensity can also be written as @xmath117    the angular integration is standard , and finally we obtain for the intensity of the radiation emitted by the quark moving in a dense medium the expression @xcite @xmath118 where the function @xmath119 is given by @xmath120 \nonumber \\ & & + \omega \ln \frac{1+\omega ^{2}\tau ^{2}\left ( 1+v_{l}\right ) ^{2}}{1+\omega ^{2}\tau ^{2}\left ( 1-v_{l}\right ) ^{2}}-2\omega v_{l}.\end{aligned}\ ] ]    in order to compute the numerical values of the intensity we need to know the mean velocities of the quarks randomly moving in a dense medium .
the multiple scattering of the quarks change the quark path and affect the emission of the radiation .
if @xmath121 is the initial velocity of the quark at the moment @xmath122 , then at the moment @xmath123 its velocity will be @xmath124 , where @xmath125 is a unit vector perpendicular to the initial direction of motion and @xmath126 is the electron multiple scattering in the time interval @xmath127 to @xmath123 @xcite . @xcite took @xmath128 .
therefore , in this approximation we obtain @xmath129 . over a distance @xmath130 the multiple scattering
is given by @xmath131 .
for @xmath28 we shall assume that it is the suppression distance , that is , the distance over the interaction actually spread : @xmath132 . in the most general case , with quantum corrections taken into account ,
@xmath133 is given as a root of a quadratic equation : @xmath134 , where @xmath135 is the classical suppression length @xcite .
in the case of the standard approach to the bremsstrahlung radiation of charged particles , one assumes that the energy is emitted in a single collision .
then , from eq .
( [ single ] ) , by taking @xmath136 and @xmath137 , it immediately follows that the single collision spectral distribution of the bremsstrahlung radiation @xmath138 is given by @xcite @xmath139    the total intensity is proportional to the frequency of the emitted radiation and is independent on the medium in which the motion of the particle takes place .
the same result can be derived by taking the limit @xmath140 in eq .
( [ int ] ) : @xmath141    in the opposite limit of very high density , corresponding to @xmath142 , from eq .
( [ int ] ) we obtain @xmath143 and @xmath144 .
therefore the electromagnetic radiation from a quark moving in a very dense medium is completely suppressed .
to calculate the bremsstrahlung emissivity a sensible model for quark interactions has to be used .
we adopt the basic assumption that the details of the interaction should be relatively unimportant , as long as the overall strength is correct and the relevant symmetries are respected . hence we use the well - known nambu - jona - lasinio ( njl ) model in its @xmath145 version in our calculations @xcite .
it has the same symmetries as qcd and describes an effective pointlike interaction with a constant coupling strength .
usually the njl model and its extensions are used in mean field calculations . as a consequence of adopting this model all the aspects related to color superconductivity
are neglected .
let @xmath146 be the cross - section for a given process of scattering of charged particles , which may be accompanied by the emission of a certain number of photons .
for example , @xmath146 could refer to the scattering of a quark by an other quark , with the possible emission of hard photons .
together with this process one could consider another process which differs from it only in that one extra photon is emitted . in this case
the total cross - section @xmath147 can be represented as a product of two independent factors , the cross - section @xmath146 and the probability @xmath148 of emission of a single photon in the collision @xcite .
the emission of a soft photon is a quasi - classical process .
the probability of emission is the same as the classically calculated number of quanta emitted in the collision , that is the same as the classical intensity ( total energy ) of emission @xmath149 , divided by the frequency of the radiation @xmath15 @xcite .
hence as a first step in obtaining the bremsstrahlung emissivity of quark matter we have to calculate the classical radiation intensity emitted by a quark moving in a dense medium in which many inter - particle collisions occur .
the probability @xmath148 of emitting a photon of energy @xmath15 by any of the scattered quarks is @xmath150 .
therefore the total cross section for the emission of soft bremsstrahlung photons is given by @xcite @xmath151    to calculate the total cross section for a photon emission we need to know @xmath146 , describing the scattering of quarks .
within the two - flavor njl model the elementary quark - quark cross sections have been calculated by @xcite.the calculations have been done nonperturbatively in the coupling constant , using the so - called @xmath152 expansion , where @xmath153 is the number of colors . in the first order of the @xmath152 expansion and for temperature @xmath154 , where @xmath155 gev , the quark - quark scattering cross - sections have a remarkably simple dependence on the collision energy @xmath29 , given by @xmath156 @xcite .
the cross section decreases with increasing energy . for a more general parametrization of the quark - quark scattering cross section , with the inclusion of the double scattering terms
see also @xcite .
we assume that the quarks form a degenerate fermi gas , with the particle number density @xmath88 given by @xmath157 p^{2}dp=\frac{g}{2\pi ^{2}}\int_{0}^{\infty } \frac{p^{2}dp}{\exp \left [ \frac{e(p)-\mu } { t}\right ] + 1},\ ] ] where @xmath158 is the quark kinetic energy and @xmath159 = 1/\left ( \exp \left [ \frac{e(p)-\mu } { t}\right ] + 1\right ) $ ] is the fermi - dirac distribution function @xcite .    in order to estimate the importance of the lpm effect on the bremsstrahlung emissivity of quark matter we have to estimate @xmath86 , the average time between two collisions . in the case of
the interaction of two quarks , denoted 1 and 2 , the average time can be obtained from @xcite @xmath160 , \ ] ] where @xmath161    the kinematic invariant @xmath29 is given by @xmath162 @xcite . by assuming that quarks are ultrarelativistic , with @xmath163 , the average collision time
can be approximated by the following analytic representation : @xmath164 , \ ] ] where @xmath165 is the polylogarithm function defined as @xmath166 .
the variation of the collision time as a function of the quark number density is presented , for different values of the temperature , in fig .  2 .
in the framework of our model and for the range of temperatures considered , the collision time is almost independent of @xmath167
. as expected , @xmath86 decreases with increasing density . to calculate @xmath86
we have assumed again that the test quarks , moving in the dense medium , have energies of the same order as the fermi energy , @xmath168 .
a rough estimate of the mean collision time can be also obtained from the relation @xmath169 . by assuming that the cross section of quark - quark interactions is around @xmath170 mb @xcite and the quarks are travelling with the speed of light
, we obtain @xmath171 fm @xmath172 for @xmath173 @xmath24 and @xmath174 fm @xmath175 for @xmath176 @xmath24 .
the knowledge of @xmath86 allows the calculation of the spectrum of the photons emitted via the bremsstrahlung mechanism from the quark star surface .
the applicability of the formulas for the spectral distribution and intensity of the bremsstrahlung radiation we have derived in the previous section is limited by the quantum condition that @xmath15 be small as compared with the total kinetic energy of the quark , which is of the order of @xmath177 : @xmath178 .
since at the star s surface the kinetic energy of the quarks is of the order of @xmath179 mev , the radiation formulae are valid only for frequencies @xmath180 mev . in order to find the maximum frequency @xmath181 of the radiation in the ultrarelativistic case
we assume that it is of the order @xmath182 , that is it is inversely proportional to the duration of the collision and to the square of the energy of the radiating particle @xcite .
if we assume that quarks have energies near the fermi energy , then the corresponding velocities are the fermi velocities , @xmath183 which are of the order of @xmath184 for @xmath185 and @xmath186 for @xmath187 .
hence the maximum allowable frequency up to which the formalism can be applied is given by @xmath188    for a quark - quark interaction cross sections of the order of @xmath189 mb and for @xmath185 we obtain @xmath190 mev . in the following
we shall restrict our study only to quark bremsstrahlung frequencies satisfying the condition @xmath191 mev .
the frequency distribution of the radiation emitted by a single quark moving in the dense matter at the surface of the quark star is presented in fig .
the total intensity @xmath192 of the bremsstrahlung radiation emitted by a single quark , @xmath193 obtained by integrating eq .
( [ int ] ) over all possible frequencies is represented , as a function of the collision time @xmath86 , in fig .
4 .    for the same range of frequencies the calculation of the bremsstrahlung intensity by using eq .
( [ bremst ] ) , which does not take into account the effects of the multiple collisions , gives @xmath194 mev . by taking into account that in the high density quark - gluon plasma the collision time is finite and of the order of @xmath195 fm leads to an intensity of the order @xmath196 mev .
hence the lpm effect could reduce the intensity of the photon bremsstrahlung of quarks at the surface of the star by a factor of @xmath197 to @xmath198 for @xmath199 @xmath24 .
for collision times of the order of @xmath200 fm , which , due to the incertitudes in the quark - quark cross - sections , are also possible , the decrease in the intensity of the bremsstrahlung is of an order of magnitude .
the variation of the intensity of the electromagnetic radiation is represented , as a function of the quark density , and for different values of the temperature , in fig .  5    the intensity of
the radiation emitted by a single particle is decreasing with increasing density , due to the effects of the multiple scatterings of the quarks .
the number of photons is proportional to the number of collisions , but if these happen too often , due to the high density of the medium , the photons emitted at different points of the trajectory start to interfere with each other and the intensity of the radiation is accordingly reduced . on the other hand @xmath192 increases with the temperature of the dense quark matter .
if the quarks are extremely relativistic , with @xmath201 , the intensity of the total radiation emitted by a single quark can be expressed in the following closed form @xmath202 +   \nonumber \\ & & \frac{1}{\tau } \left [ \arctan 2\left ( 1 + 3\pi /4\alpha _
{ s}\right ) -\arctan 2\tau an_{b}^{1/3}\right ] - \nonumber\\ & & 4\left ( \frac{1 + 3\pi /4\alpha _ { s}}{\tau }
-an_{b}^{1/3}\right ) + \frac{1}{2\tau } \left [ d\left ( 1 + 3\pi /4\alpha
_ { s}\right ) -d\left ( \tau an_{b}^{1/3}\right ) \right ] , \end{aligned}\ ] ] where the function @xmath203 is given by @xmath204 $ ] .
@xmath203 is a real function for all @xmath205 , and its power series expansion is given by @xmath206 .
the approximation @xmath207 is appropriate for quarks having energies near the fermi energy and moving with velocities @xmath208 , which are of the order of @xmath189 .
the emission of a photon in the quark matter is the result of the scattering of a quark from a state @xmath209 to a state @xmath210 by collision with another particle of momentum @xmath211 .
the rate of collisions in the quark - gluon plasma in which a soft photon of energy @xmath55 is emitted can be defined as @xcite @xmath212 \left [ 1-f_{2}\left ( \vec{p}_{2}^{\prime } \right ) \right ] \prod_{i=1,2,1^{\prime } , 2^{\prime } } d^{3}\vec{p}_{i}d\omega \left ( \theta , \phi \right ) , \ ] ] where @xmath213 is the solid angle element in the direction @xmath214 , @xmath215 are the statistical weights and @xmath216 , @xmath217 are the velocities of the interacting particles .
the bremsstrahlung emissivity per unit volume of the quark matter is given by @xmath218    withe the use of the explicit form of the cross - section for the photon emission , given by eq .
( [ cross ] ) , the photon emissivity can be written in the following form : @xmath219 where @xmath220 is a function of the quark density and of the temperature .    the bremsstrahlung energy flux from the star
, @xmath221 is given by the emissivity of a thin surface layer , with thickness of the order of the photon mean free path @xmath222 @xcite .
taking into account only photons emitted outwards and produced not deeper than @xmath222 , the energy flux from the surface of the quark star is given by @xmath223 where @xmath224 .
since the temperature dependence of @xmath225 is not significant , it follows that the quark bremsstrahlung spectrum has the same temperature dependence as the black body radiation
.    for arbitrary temperatures the coefficient of the bremsstrahlung emissivity can be represented , in an approximate analytical form , by @xmath226 i\left ( n_{b},t\right ) \lambda .\ ] ]    the calculation of @xmath227 requires , in the general case , the estimation of integrals of the form @xmath228 $ ] , which can be expressed usually in terms of the polylogarithm function . in the limit of small temperatures , with @xmath229 ,
the integral can be evaluated by means of the substitution @xmath230 . for small
@xmath167 one can extend the lower limit of integration to @xmath231 , so that the integral becomes @xmath232 = \left ( \pi ^{2}/6\right ) t^{2}$ ] .
therefore in the limit of small temperatures @xmath233 , but there is no explicit temperature dependence of the coefficient @xmath225 : @xmath234    in this case the only @xmath167-dependence of @xmath225 is via the intensity of the radiation .
the variation of the ratio of the bremsstrahlung emissivity coefficient and the stefan - boltzmann constant @xmath235 is represented , as a function of the mean collision time @xmath86 , and in the limit of small temperatures , in fig .
we have assumed that the thickness of the emitting layer is @xmath236 fm .    in the limit of large collisions times @xmath237 , and thus we recover the single - collision results of @xcite . for collision times specific to the quark - gluon plasma at the surface of the quark stars , of the order of @xmath238 fm , the ratio @xmath239 .
therefore , taking into account multiple collisions in the high density plasma , reduces with an order of magnitude the value of the bremsstrahlung coefficient @xmath225 , making the soft photon radiation from the surface of the quark star about @xmath198 orders of magnitude weaker than the equilibrium black body radiation .
the radiation emitted by the thin quark surface layer of the strange star must propagate through the electron layer formed due to the electrostatic attraction of the quarks . in the followings we shall consider the effect of the electron layer on the electromagnetic radiation emitted , via the bremsstrahlung mechanism , by the quarks in the layer of thickness @xmath222 . as a first step it is necessary to obtain a quantitative analytical description of the properties of the electrons outside the quark star .    to find the electron distribution near the quark star surface we follow the approach of @xcite , but also taking into account the finite temperature effects as discussed by @xcite .
the chemical equilibrium implies that the electron chemical potential @xmath240 is constant , where @xmath241 is the electrostatic potential per unit charge and @xmath242 is the electron s chemical potential .
since far outside the star both @xmath241 and @xmath242 tend to zero , it follows that @xmath243 and @xmath244 .    for the number density @xmath245 of the electrons at the quark star surface we use the expression @xcite @xmath246 where @xmath167 is the temperature of the electron layer , which can be taken as a constant , since we assume the electrons are in thermodynamic equilibrium with the constant temperature quark matter .
the poisson equation for the electrostatic potential @xmath247 generated by the finite temperature electron distribution reads @xcite @xmath248 , z\leq 0 , \label{2}\ ] ] @xmath249    in eqs .
( [ 2])-([3 ] ) , @xmath250 is the space coordinate measuring height above the quark surface , @xmath97 is the fine structure constant and @xmath251 is the quark charge density inside the quark matter
. the boundary conditions for eqs .
( [ 2])-([3 ] ) are @xmath252 as @xmath253 and @xmath254 for @xmath255 . in the case of the zero temperature electron distribution at the boundary @xmath256
we have the condition @xmath257 @xcite .    in the following
we will be interested mainly in the properties of the electron gas extending outside the quark surface of the star , described by eq .
( [ 3 ] ) . in order to find the exact solution of eq .
( [ 3 ] ) , we introduce first a new variable @xmath258 .
then eq .
( [ 3 ] ) can be written in the form @xmath259 immediately leading to the first integral @xmath260 where an arbitrary integration constant has been set to zero .
the general solution of eq .
( [ eqv ] ) is given by @xmath261 } { \exp \left [ 4\sqrt{\frac{\alpha \pi } { 3}}t\left ( z+z_{0}\right ) \right ] -1},\ ] ] where @xmath262 is a constant of integration .
its value can be obtained from the condition of the continuity of the potential across the star s surface , requiring @xmath263 , where @xmath264 is the value of the electrostatic potential in the region @xmath265 , described by eq .
( [ 2 ] ) .
therefore @xmath266 .\ ] ]    hence the electron number density @xmath267 is given , at finite temperatures , by @xmath268 } { \exp \left [ 4\sqrt{\frac{\alpha \pi } { 3}}t\left ( z+z_{0}\right ) \right ] -1}\left\ { 8\left (   \frac{\exp \left [ 2\sqrt{\frac{\alpha \pi } { 3}}t\left ( z+z_{0}\right ) \right ] } { \exp \left [ 4\sqrt{\frac{\alpha \pi } { 3}}t\left ( z+z_{0}\right )
\right ] -1}\right ) ^{2}+1\right\ } t^{3},t\neq 0.\ ] ]    in the limit of zero temperature , @xmath269 , from eq .
( [ eqv ] ) we obtain @xmath270 where @xmath271 and @xmath272 is an integration constant .
@xmath272 can be determined from the boundary condition @xmath257 , which gives @xmath273 . therefore , in this case for the electron particle number distribution , extending several thousands fermis above the quark matter surface , we find the expression : @xmath274    the variation of the electron number density @xmath275 is represented , for different values of the temperature , in fig .
7 .    in the absence of a crust of the quark star
, the electron layer can extend to several thousands fermis outside the star s surface .    the intensity of the photon beam coming from the thin outer layer of the quark star decreases due to the energy absorption in the electron layer .
the variation of the intensity of the beam can be obtained by solving the equation of the radiative transfer @xcite . for a steady state ,
the general solution of the transport equation is by given @xmath276 where the optical depth @xmath277 is defined as @xmath278 , with @xmath279 the rosseland mean opacity ( giving the energy absorbed per unit time per unit volume from the beam with intensity @xmath280 ) , @xmath281 is the intensity of the radiation incident on the electron layer and @xmath282 is the emissivity of the electron . since in the present paper we are interested to study only the effect of the outer electrons on the quark matter emission , we shall take @xmath283 , that is we shall neglect the emissivity of the electron layer .
hence , the variation of the intensity , due to absorption in the electron layer , of the photon beam emitted by the quarks at the strange star surface can be described by the following equation : @xmath284    in the case of the scattering of photons by free non - relativistic electrons , the opacity is given by @xmath285 , where @xmath286 is the thomson cross - section @xcite .
this formula is valid only in the case of a non - degenerate gas at low - temperatures . for systems at high densities and/or temperatures , the corrections to the thomson scattering cross - section
become important .    for high - temperature non - degenerate electrons
the opacity has been calculated first in @xcite .
these calculations have been extended to the semi - degenerate ( @xmath287 ) case by @xcite and @xcite ( see also @xcite and @xcite ) . for a degenerate electron gas at @xmath287 , the rosseland mean opacity
is given by @xcite @xmath288    @xmath289 ( the inverse of the rosseland mean ) is a function of the temperature and of the degeneracy parameter @xmath290 , where @xmath291 is the fermi energy .
@xmath292 can be represented in the form @xcite @xmath293 where @xmath294 .
the variation of the function @xmath295 as a function of the parameters @xmath296 and @xmath167 is represented in fig .
8 .    for a non - degenerate electron gas ,
obeying the maxwell - boltzmann distribution , the rosseland mean opacity is given by @xcite @xmath297 where @xmath298 and the temperature is expressed in units of @xmath299 .    in order to obtain a simple estimate of the absorption of the surface radiation of quarks due to the presence of the electron layer , we consider first the case of the
completely degenerate electron gas at a small temperature @xmath167 .
complete degeneracy corresponds to @xmath300 and @xmath301 .
we also neglect the temperature dependent terms in eq .
( [ g ] ) , which are not important for small temperatures .
hence from eq .
( [ g ] ) we obtain , as a first approximation , @xmath302 .
therefore the rosseland opacity of the electron layer at the quark star surface can be approximated , for low temperatures of the degenerate electron gas , by the expression @xmath303    consequently , for the optical depth of the electron distribution above the quark star surface we find the expression @xmath304 .
\label{tau0}\ ] ]    for large @xmath250 , corresponding to the limit @xmath255 , from eq .
( [ tau0 ] ) we obtain @xmath305    at zero temperature the optical depth is proportional to the square of the value of the quark charge density at the surface of the quark star @xmath306 .
the absorption of the bremsstrahlung radiation by the electron layer essentially depends on the quark charge density @xmath307 inside the quark star , where @xmath27 , @xmath28 and @xmath29 denotes the different quark flavors @xcite .
the condition of the charge neutrality requires @xmath308 , where @xmath309 is the charge density of the electrons inside the star .
@xmath309 is related to the baryon number density @xmath310 by means of the relation @xmath311 , where for the electron abundance @xmath312 we adopt the range @xmath313 @xcite . near the surface of the quark star
the density can be approximated by @xmath314 g/@xmath315 . for a strange quark mass of the order of @xmath316 mev , for the baryon number density near the surface of the star we find @xmath317 mev@xmath318 .
hence the corresponding quark charge density is @xmath319 mev , leading to @xmath320 mev .
the corresponding bremsstrahlung luminosity @xmath321 of the quark star is given by @xmath322 , where @xmath323 is the radius of the star .
for a @xmath324 star and for a degenerate electron layer , @xmath325 .\ ] ]    the variation of the coefficient @xmath326 is represented , as a function of @xmath306 , in fig .
9 .    for an increasing surface quark charge density @xmath326 is a rapidly decreasing function .
therefore in the case of a large @xmath306 the electron layer significantly decreases the intensity of the bremsstrahlung radiation emitted by the thin quark layer at the surface of the star .
the spectrum of the @xmath287 quark star essentially depends on the charge density @xmath327 at the surface of the star , on the degeneracy parameter @xmath296 and on the temperature .
the numerical values of @xmath328 must be obtained by solving eq .
( [ 2 ] ) in the region @xmath329 for different values of @xmath167 , together with the boundary condition @xmath252 for @xmath253 .
therefore even for the non - zero temperature quark star , the spectrum of the bremsstrahlung radiation is determined by the quark charge density inside the star . however , in order to obtain a simple estimate of the temperature effects on the optical length of the radiation in the electron layer , we shall assume that for a suitable range of temperatures the charge density at the surface can be approximated by @xmath330 mev @xcite .
for this value of the charge density the variation of @xmath326 as a function of @xmath296 and @xmath167 is represented in fig .  10 .
since the increase in the temperature also increases the electron number density of the electron layer , the absorption of the bremsstrahlung radiation is much stronger at high temperatures of the semi - degenerate electron gas .
since the proposal of the existence of strange stars , much effort has been devoted to find observational properties that could distinguish strange stars from neutron stars .
one of the most common and simplest method could be to find some specific signatures in the photon emissivity of strange stars . in the first estimate of the surface photon emissivity of quark matter @xcite
it has been concluded that the most important photon radiation mechanism is the bremsstrahlung of quarks from a thin layer at the star s surface , and that the intensity of the radiation is about four orders of magnitude weaker than the equilibrium black body radiation , although both have the same temperature dependence .    in the present paper
we have reconsidered the electromagnetic emission from the surface of bare quark stars , by also taking into account the important effect of the multiple collisions of quarks in the very high density quark - gluon plasma existing at the surface of the quark stars . to obtain the bremsstrahlung spectrum of the photon emission of the quarks we have used several assumptions .
the most important is the adoption of a specific model for the description of quark - quark elastic scattering , which is essential in obtaining the bremsstrahlung emissivity .
the details of the interactions between quarks are generally poorly known , and therefore there are a lot of uncertainties in the description of the strong interactions in a quark - gluon plasma . in our calculations
we have used a classical soft photon approximation , neglecting quantum effects .
we have generally assumed that the quarks are semi - degenerate , with temperature @xmath331 and @xmath332 , but with the kinetic energies of the interacting particles of the order of the fermi energy .    in order to give a complete description of the quark - quark bremsstrahlung we have considered , in a systematic way , the effects of the dense medium on the radiation spectrum .
a dense medium in which many collisions occur will reduce the bremsstrahlung radiation , and as high the density is , as lower is the bremsstrahlung emission . on the other hand
the emission of photons with frequency smaller than the critical lpm frequency is completely suppressed .
@xmath85 depends on the qcd coupling constant and on the density at the quark star s surface . for strange stars
@xmath85 is of the order of @xmath333 mev , so photons with @xmath334 can not be emitted from the surface of the star .
hence only low frequencies are suppressed due to the lpm effect . on the other hand
the multiple collision effect reduces with one order of magnitude the bremsstrahlung coefficient @xmath225 , making the soft photon radiation from the quark star @xmath198 orders of magnitude weaker than the black body radiation .    a second effect , which essentially influences the radiation of photons by quark stars , is the absorption of the radiation in the electron layer formed outside the quark star . the electron layer could extend to several hundreds or thousands fermis , the electron number density being an increasing function of the temperature .
the high - density semi - degenerate electron layer at the surface of a strange star plays an important role in the propagation of the electromagnetic radiation from the star s surface .
the optical depth of the degenerate or semi - degenerate electron layer depends on the electron number density , temperature , the degeneracy parameter and the charge density at the surface of the quark star . for high surface charge densities the intensity of the quark bremsstrahlung radiation is significantly reduced by absorption in the electron layer . for the low - temperature star , with @xmath335 , for a surface potential @xmath336 mev ,
the electron layer reduces the intensity of the radiation by an order of magnitude . in this case
the combined effects of the multiple collisions and absorption in the electron layer give @xmath337 .
hence the energy radiated away by the quark star could be @xmath338 orders of magnitude smaller than the equilibrium black body radiation .    for the sake of comparison we also present the expression for the photon emission rate of an equilibrated
, hot qcd plasma , at zero chemical potential , with the lpm effect fully included .
the spontaneous emission rate @xmath339 of photons with a given momentum @xmath340 can be represented , for two flavor qcd , as ( arnold et al .
2001 , arnold et al . 2002 ) @xmath341 , \label{ar1}\ ] ] where @xmath342 is the fermi - dirac distribution function , @xmath343 ^{-1}$ ] , @xmath344 is the thermal quark mass , given by @xmath345 , and @xmath346 with @xmath347 and @xmath348 @xmath349 $ ] . for @xmath350
, bremsstrahlung becomes the most important process , while for @xmath351 , pair annihilation dominates .
the lpm effect suppresses both bremsstrahlung and pair annihilation processes , but the suppression is not severe ( @xmath352 or less ) , except for photons with @xmath353 or very hard pair annihilation .    in the limit of small temperatures , @xmath269 and @xmath354 .
hence @xmath355 and @xmath356 .
by neglecting the logarithmically divergent terms , the integration over all possible ranges of frequencies in eq .
( [ ar1 ] ) gives @xmath357 where @xmath358 /9 $ ] , with @xmath359 is the riemann zeta function @xmath360 .
therefore the perturbative qcd approach gives , in the low temperature limit , the same temperature dependence of the bremsstrahlung spectrum as in eq .
( [ arbr ] ) .
however , we must point out that eqs.([ar1 ] ) and ( [ ar2 ] ) have been derived for a temperature of the quark - gluon plasma higher than @xmath361 mev , since the temperature of the fireball formed in heavy ion collisions is of the order of @xmath362 mev ( renk 2003 ) .
recently , it has been argued that the strange quark matter in the color - flavor locked(cfl ) phase of qcd , which occurs for large gaps @xmath363 , is rigorously electrically neutral , despite the unequal quark masses , and even in the presence of electron chemical potential @xcite .
however , @xcite pointed out that for sufficiently large @xmath364 the low density regime is rather expected to be in the  2-color - flavor superconductor  phase in which only @xmath27 and @xmath28 quarks of two color are paired in single condensates , while the ones of the third color and @xmath29 quarks of all three colors are unpaired . in this phase , electrons are present . in other words ,
electrons may be absent in the core of strange stars but present , at least , near the surface where the density is lowest .
nevertheless , the presence of cfl effect can reduce the electron density at the surface and hence increases the bremsstrahlung emissivity .    in conclusion , in the present paper we have pointed out two effects that could significantly reduce or even fully suppress the bremsstrahlung radiation from the surface of quark stars : the significant decrease of the intensity due to the multiple collisions in a dense medium and the absorption by the outer electron layer of the stars . the possible observational significance for the detection of quark stars of these effects
will be considered in a future publication .
the authors would like to thank to the anonymous referee , whose comments helped to improve the manuscript .
we are grateful for the useful comments of prof . v. v. usov .
this work is supported by a rgc grant of hong kong government .
aggarwal , m. m. et al .
2000 , phys .
lett . , 85 , 3595 alcock , c. , farhi , e. & olinto , a. 1986 , apj , 310 , 261 alford , m. , rajagopal , k. & wilczek , f. 1998 , phys .
b , 422 , 247 alford , m. , rajagopal , k. & wilczek , f. 1999 , nucl . phys .
b , 537 , 443 anthony , p. l. , becker - szendy , r. , cavalli - sforza , m. , keller , l. p. , kelley , l. a. , klein , s. r. , niemi , g. , perl , m. l. , rochester , l. s. & white , j. l. 1995 , phys .
lett . , 75 , 1949 arnold , p. , moore , g. d. & yaffe , l. g. 2001 , jhep , 0112 , 009 arnold , p. , moore , g. d. & yaffe , l. g. 2002 , jhep , 0206 , 030 aurenche , p. , gelis , f. & zaraket , h. 2000 , phys .
d , 62 , 096012 baier , r. , nakkagawa , h. , niegawa , a. & redlich , k. 1992 z. phys .
c , 53 , 433 berestetskii , v. b. , lifshitz , e. m. & pitaevskii , l. p. 1982
, quantum electrodynamics ( oxford : pergamon press ) boercker , d. b. 1987 , apj , 316 , l95 buchler , j. r. & yueh , w. r. 1976 , apj , 210 , 440 cheng , k. s. & dai , z. g. 1996 , phys . rev .
, 77 , 1210 cheng , k. s. & dai , z. g. 1998 , phys . rev . lett . , 80 , 18 cheng , k. s. , dai , z. g. , wei , d. m. & lu , t. 1998a , science , 280 , 407 cheng , k. s. , dai , z. g. & lu , t. 1998 , int . j. mod .
d , 7 , 139 cheng , k. s. & harko , t. 2000 , phys .
d , 62 , 083001 chin , c. w. 1965 , apj , 142 , 1481 chiu , h .- y .
1968 , stellar physics ( waltham , massachussets : blaisdell publishing company ) chmaj , t. , haensel , p. & slominski , w. 1991 , nucl .
b ( proc . suppl . ) , 24 , 40 cleymans , j. , goloviznin , v. v. & redlich , k. 1993 , phys . rev .
d , 47 , 173 dai , z. g. , peng , q. h. & lu , t. 1995 , apj , 440 , 815 dey , m. , bombacci , i. , dey , j. , ray , s. & samanta , b. c. 1998 , phys .
b , 438 , 123 gondek - rosinska , d. , bulik , t. , zdunik , l. , gourgoulhon , e. , ray , s. , dey , j. & dey , m. 2000 , a&a , 363 , 1005 gourgoulhon , e. , haensel , p. , livine , r. , paluch , e. , bonazzola , s. & marck , j .- a .
1999 , a&a , 349 , 851 haensel , p. , zdunik , j. l. & schaeffer , r. 1986 , a&a , 160 , 121 haglin , k. & pratt , s. 1994 , phys .
b , 328 , 255 harko , t. & cheng , k. s. 2002 , a&a , 385 , 947 hatsuda , t. & kunihiro , t. 1994 , phys .
rept . , 247 , 221 jackson , j. d. 1999 , classical electrodynamics ( new york : john wiley and sons ) kapusta , j. , lichard , p. & seibert , d. 1991 phys .
d , 44 , 2774 kettner , ch . , weber , f. , weigel , m. k. & glendenning , n. k. 1995 , phys .
d , 51 , 1440 klein , s. 1999 , rev .
phys . , 71 , 1501 landau , l. d. & lifshitz , e. m. 1975 , the classical theory of fields ( oxford : pergamon press ) landau , l. d. & pomeranchuk , i. j. 1953 , dokl .
sssr , 92 , 735 manley , d. m. 1997 , phys .
d , 56 , 7325 page , d. & usov , v. v. 2002 , phys .
lett . , 89 , 131101 peitzmann , t. & thoma , m. h. 2002 , phys . rept . , 364 , 175 rajagopal , k. & wilczek , f. 2001 , phys . rev .
, 86 , 3492 renk , t. 2003 , hep - ph/0301133 , to appear in phys .
c rose , s. j. 1995 , apj , 453 , l45 rybicki , g. b. & lightman , a. p. 1979
, radiative processes in astrophysics ( new york : john wiley and sons ) sampson , d. h. 1959 , apj , 129 , 734 shapiro , s. l. & teukolsky , s. a. 1983 , black holes , white dwarfs , and neutron stars ( new york : john wiley and sons ) usov , v. v. 1998 , phys .
, 80 , 230 witten , e. 1984 , phys .
d , 30 , 272 zdunik , j. l. 2000 , a&a , 359 , 311 zhuang , p. , hufner , j. , klevansky , s. p. & neise , l. 1995 , phys .
d , 51 , 3728