%%    TEMPLATE for articles submitted to the full-tex econf proceedings
%%
%%
%%     Please do not remove lines commented out with %+
%%           these are for the editors' use.
%%
%%     Questions?  Send email to :  mpeskin@slac.stanford.edu
%%

\documentclass[12pt]{article}
\usepackage{epsfig}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%  Do not change these:
\textwidth=6.0in  \textheight=8.25in

%%  Adjust these for your printer:
\leftmargin=-0.3in   \topmargin=-0.20in


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%  personal abbreviations and macros
%    the following package contains macros used in this document:
\input econfmacros.tex
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%  To include an item in the INDEX of the conference volume,
%           flag it with    \index{<item name>}
%  The use of this macro is illustrated in the text.
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\def\Title#1{\begin{center} {\Large {\bf #1} } \end{center}}

\begin{document}


\Title{The QCD Phase Diagram and Explosive Astrophysics}

\bigskip\bigskip

%+\addtocontents{toc}{{\it D. Reggiano}}
%+\label{ReggianoStart}

\begin{raggedright}

{\it Stephen D.H. Hsu\index{Hsu}\\
Department of Physics\\
University of Oregon\\
Eugene, OR USA}
\bigskip\bigskip
\end{raggedright}



\section{Introduction}

I was asked by the organizers to give an overview of the QCD phase
diagram, geared towards a mixed audience of astrophysicists and
particle theorists. I chose to emphasize the phase transition(s)
from the normal nuclear phase to the color superconducting and
chirally symmetric phases at high density and low
temperature\footnote{By low temperature here I mean low relative
to the QCD scale, or $\sim 100$ MeV; by astrophysical standards
supernova temperatures of tens of MeV are obviously quite high.}.
I did so because it is possible that these transitions might occur
during violent astrophysical phenomena such as supernovae or
neutron star mergers \cite{HHS}, with potentially dramatic
consequences. If so, these astrophysical events could provide a
window into the physics of QCD.

This contribution is organized as follows. In section 2 I
summarize our current understanding of the QCD phase diagram. In
section 3 I discuss the phase transition characteristics most
relevant to astrophysics, such as critical temperature and latent
heat. In section 4 I discuss supernova core collapse and argue
that conditions at core bounce might lead to the crossing of a QCD
phase boundary. I also briefly discuss neutron star mergers and
hypernovae. I conclude with some discussion from the conference.

\section{QCD phase diagram}

Figure \ref{fig:phase} (taken from Rajagopal and Wilczek's
excellent review in \cite{RW}) shows a possible QCD phase diagram
as a function of temperature and chemical potential. The region at
low temperature and density (lower left) is the normal nuclear
matter phase (the small spur at the bottom is the phase boundary
for the formation of nuclear matter). At sufficiently high density
and low temperature (lower right) it has recently been established
that the ground state of quark matter exhibits Color
Superconductivity (CSC), resulting from the Cooper pairing of
quark quasiparticles near the Fermi surface \cite{RW} --
\cite{Interface}. At high temperature and low density (upper left,
extending to the upper right) we expect to find the quark-gluon
plasma. I will comment on the more specific features of the
diagram below.



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%
%%   use this format to include an .eps figure into your paper
%%
\begin{figure}[htb]
\begin{center}
\epsfig{file=qm99fig3.ps,height=2in} \caption{QCD phase diagram
(simplified) in the density-temperature plane.} \label{fig:phase}
\end{center}
\end{figure}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

Let me try to give you some flavor of the recent results on color
superconductivity. Recall that Cooper pairing (as in ordinary
superconductors) involves excitations near the Fermi surface, but
with equal and opposite momenta (figure \ref{fig:FS}). In QCD the
excitations have quark quantum numbers, and hence there is a
directly attractive channel due to gluon exchange in the
anti-triplet ($\bar{3}$) color representation. In the case of
electrons in a metal photon exchange is repulsive, and the
attractive channel is due to the exchange of phonons (lattice
vibrations).

What is common between quark matter and an ordinary metal is the
existence of a Fermi surface. The Pauli exclusion principle
requires that fermions like quarks or electrons occupy distinct
quantum states, hence at zero temperature the ground state
consists of filled levels up to some energy (this boundary is the
Fermi surface). The lowest energy excitations of this system,
called quasiparticles, are states just above the Fermi surface.
Now, a Fermi surface is unstable with respect to attractive
interactions in the Cooper pairing channel. That is, even an
arbitrarily weak interaction in this channel can lead to pairing
of quasiparticles, which reduces the overall energy of the system.
Heuristically, the spherical symmetry of the Fermi surface (in
momentum space, the energy of a quasiparticle only depends on its
distance from the surface, and not on its angular position)
reduces the pairing dynamics to that of a 1+1 dimensional system,
where even weak interactions lead to pairing.


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%
%%   use this format to include an .eps figure into your paper
%%
\begin{figure}[htb]
\begin{center}
\epsfig{file=FS.ps,height=2in} \caption{Excitations on opposite
sides of the Fermi surface.} \label{fig:FS}
\end{center}
\end{figure}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%



Now, recall that QCD exhibits a property called asymptotic
freedom, which means that short distance (large momentum transfer)
interactions are weak, whereas long distance (small momentum
transfer) interactions are strong. At very high quark densities,
most interactions occur over short distances, and we therefore
have very good control over calculations. However, at quark
densities likely to be found in a neutron star, roughly a few per
cubic Fermi, the typical interactions are quite strong, precluding
reliable quantitative results. Nevertheless, the indications of at
least an {\it instability} to the formation of Cooper pairs is
still evident, even if we can't say more about the details.

In figure \ref{fig:phase} there are two distinct CSC phases
displayed, the 2SC and CFL phases. In the former, the up and down
quarks pair in the isosinglet channel, and the strange quark is
left to pair with itself in an exotic (possibly spin 1) channel.
In the latter, the flavor and color orientations of all three
quarks are correlated in a non-trivial way (see figure
\ref{fig:N3}): pairing in a particular flavor channel corresponds
to a particular orientation in color space. In this talk I will
assume that figure \ref{fig:phase} is correct in that the 2SC
phase has the lowest energy at intermediate density. The
transition from nuclear matter to the 2SC phase is likely to be
first order \cite{Universal}, as we will discuss further below.



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%
%%   use this format to include an .eps figure into your paper
%%
\begin{figure}[htb]
\begin{center}
\epsfig{file=N3.ps,height=1in} \caption{Di-quark orientation in
the CFL phase. The arrow represents a $\bar{3}$ orientation in
color space.} \label{fig:N3}
\end{center}
\end{figure}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


Before leaving the subject of the QCD phase diagram, let me
discuss another transition -- the chiral phase transition -- that
is likely to occur as we increase the baryon number density.
Chiral symmetry has to do with unitary rotations among different
flavors of quarks:
\begin{equation}
q_i^{L,R} \rightarrow U_{ij} q_j^{L,R}~~~,
\end{equation}
where i labels the flavor (i = up, down, strange) and $U$ is a
unitary matrix. L and R label whether the quark is left handed or
right handed. Chiral symmetry reflects the fact that, in the limit
of small quark masses, the strong interactions (gluons) can't tell
the difference between different flavors of left and right handed
quarks.

These symmetries are broken by quark condensation in the QCD
ground state:
$$
\langle \bar{q}_L q_R \rangle \neq 0~~~,
$$
but are restored in the 2SC phase and in the quark-gluon plasma.
(However, in the CFL state chiral symmetries do remain broken at
high density.) Their restoration involves a phase transition,
again at quark densities of roughly a few per cubic Fermi. What
can we say about this transition?

Because it occurs at a density where the effective coupling is
strong, we can't make quantitative predictions. However, we can
learn something about the order of the transition from something
called ``universality''. The logic is as follows: if the
transition is second order (or higher), then at the precise point
of transition there must be very long wavelength excitations in
the system. A heuristic way to understand this is to visualize the
effective potential near the transition. In a second order
transition the potential becomes very flat as the second
derivative changes sign from positive to negative. (In contrast,
in a first order transition the system tunnels or fluctuates out
of the disfavored vacuum while the second derivative is still
positive.)

Technically, a second order transition implies fluctuations of
infinite wavelength (or equivalently, zero mass excitations).
There are not very many models which describe such low-energy
dynamics. To be self-consistent they need to exhibit something
called an infrared fixed point in the renormalization group
evolution of the coupling constant. Roughly speaking, the
low-energy dynamics must become almost scale invariant as the
correlation length diverges. The candidate models describing the
low-energy dynamics must possess the same symmetry properties as
the underlying system, but otherwise can be quite simple. In some
cases, there are no candidate models with the right symmetries and
low-energy particle content that exhibit an infrared fixed point.
In this case the transition is predicted to be first order. Based
on this kind of analysis, both the 2SC and chiral phase transition
at finite density are likely to be first order
\cite{ChPT,Universal}.

So, we have two candidate phase transitions which, as we discuss
below, might have astrophysical implications. Obviously, a number
of issues are still unresolved. For example, are there two
separate phase transitions, or does chiral symmetry restoration
coincide with color superconductivity? The renormalization group
analysis of color superconductivity \cite{ESH,SW} suggests that
once the effective excitations start to look like quarks, there is
a pairing instability and the system is likely to be
superconducting. However, one could also imagine a chirally
restored phase with nucleonic excitations, so chiral symmetry does
not necessarily imply CSC. In any case, both transitions are
likely to occur at temperatures and densities that might be
achieved by astrophysical phenomena.





\section{Phase transition parameters}

At asymptotic density, the nature of the CSC order parameter, the
binding energy density and the critical temperature are all
precisely calculable. Less is known about intermediate densities
of several times nuclear density. However, there are strong
indications of a Cooper pairing instability, and estimates of the
resulting gap are of order $\Delta \sim$ 40-120 MeV at a quark
chemical potential $\mu$ of $\approx$ 400 MeV. The nature of the
chiral phase transition is similarly poorly determined, although
as we discussed in the previous section it is likely to be first
order. Rather than discuss the two transitions separately, we will
focus on the CSC transition, keeping in mind that if a first order
chiral boundary is crossed, the consequences should be similar.

Consider the possible implications of CSC for the collapse and
explosion of massive stars ($M > 8 M_{\odot}$). We will argue that
it is quite likely that at the moment of maximum compression of
the collapsing Fe core, the densest part of the star crosses the
critical density into the phase where CSC is energetically
favored. The release of latent energy has the potential to
generate an explosive shockwave which powers the resulting
supernova (SN). Current simulations of supernovae are generally
unable to reproduce the explosive behavior observed in nature: the
shockwave generated by the mechanical bounce of the nuclear core
stalls before reaching the surface, unless an appeal is made to
neutrino reheating combined with non-spherical phenomena such as
rotation or convection \cite{SNcollapse,SNsim}. We also note that
the energy liberated in a CSC phase transition is potentially
sufficient to power hypernovae (HN) \cite{Hyper}, which have been
linked to gamma ray burst events \cite{GRB}.



First, let us summarize some results on CSC from the recent
literature. Precise results are only valid at asymptotic densities
where the effective QCD coupling is small, however they should
still be useful guide when dealing with intermediate densities. In
any case, our qualitative results will be insensitive to factors
of 2 in these formulas:

\bigskip

\noindent $\bullet$ Gap size: $\Delta \sim$ 40 - 120 MeV

\smallskip

\noindent $\bullet$ Critical temperature: $T_c \simeq .57 \Delta$

\smallskip

\noindent $\bullet$ At asymptotic density the binding energy
density is $E_{CSC} = {\frac{5.8}{4 \pi^2}} \Delta^2 \mu^2$.
Simple dimensional analysis (given the absence of any small
parameter) also suggests a value in this range. Note that we are
interested here in the latent {\it energy} associated with the
first order transition to the CSC phase (or any other transition
that occurs as the baryon density is increased beyond several
times nuclear density). The baryon density changes on
astrophysical timescales, or very slowly on the timescale of QCD
dynamics, and the vacuum state at zero temperature (or at $T <<
T_c$) is found by minimizing the energy in the sector of the
Hilbert space with fixed baryon number. Studies such as
\cite{Interface}, which involve the free energy $\Omega$ at finite
chemical potential (but not fixed density) are appropriate for
determining pressure equilibrium between nuclear and CSC matter in
circumstances in which baryon number can flow across a boundary
(e.g. in a neutron star), but do not address a possible SN
transition. (See discussion for further remarks.)

\smallskip

\noindent $\bullet$ Phase diagram: the normal nuclear phase is
separated from the CSC phase (most likely the 2SC two flavor
condensate phase, although possibly the CFL phase
\cite{ChPT,Universal,Interface}) by a first order boundary.

\bigskip

\section{Astrophysics of Core Collapse}

Now let us review the standard scenario of Fe core collapse which
is believed to lead to type II supernovae \cite{SNcollapse}.
Nuclear burning during the $10^7$ year lifetime of the star leads
to a shell structure, with the inner core eventually consisting of
Fe ash. Because iron cannot participate in further exothermic
nuclear reactions, there is an eventual cooling and collapse of
the Fe core, whose mass is likely to be $(1-2)~ M_{\odot}$ (or,
roughly the Chandresekhar mass). The collapse of this core is only
halted by neutron degeneracy, which leads to a stiffening of the
equation of state. The resulting bounce produces a shockwave,
whose energy of $\sim 10^{51}$ ergs is a small fraction of the
total available gravitational binding energy released by the
collapse:
\begin{equation}
E_b \sim 3 \cdot 10^{53} {\rm ergs} \left( {\frac{M_{core}
}{M_{\odot}}} \right)^2 \left( {\frac{R }{10 {\rm km}}}
\right)^{-1} \ . \end{equation} Most of this energy escapes in the
form of neutrinos during the supernova, as was observed in the
case of SN1987a.

The pressure in the collapsed core at the instant of the bounce is
most likely {\it greater} than the corresponding pressure in any
remnant neutron star. In order to cause a bounce, the kinetic
energy $E_b$ of the infalling material (which is a sizeable
fraction of a solar mass!) must be momentarily stored as
compressional potential energy in the (sub)nuclear matter. This
additional mechanical squeezing at the bounce suggests that if the
critical density for CSC is ever reached in a neutron star, it
will be reached at this instant.

Simulations of the core bounce result in densities of at least
several times nuclear density ($5-10 \cdot 10^{14} {\rm g/cm^3}$),
and temperatures of roughly 10-20 MeV \cite{SNcollapse}. This
temperature is likely less than $T_c$ for CSC, possibly much
less\footnote{It is conceptually easier to think about the case
where $T$ is much less than $T_c$, since in this case the Free
energy ($F = U - TS$) liberated by the transition is predominantly
energy, with only a small component related to entropy. The
relevant dynamics is governed by energetics rather than Free
energetics.}, and hence the core of the star traverses the phase
diagram horizontally in the density-temperature plane, crossing
the critical density boundary into the CSC phase. It is important
to note that the core region at bounce is probably {\it cooler}
than post-bounce, since degenerate neutrinos tend to heat the
proto-neutron star as they diffuse out \cite{Lattimer}. Studies
quoting larger SN temperatures such as $T \sim 30$ MeV generally
refer to this later stage \cite{Carter}.

Once the core crosses into the CSC part of the phase diagram, the
transition proceeds rapidly, on hadronic timescales. Because the
transition is first order, it proceeds by nucleation of bubbles of
the CSC phase in the normal nuclear background. The rate of bubble
nucleation will be of order $({\rm fm})^4$ (${\rm
fm}=10^{-13}$cm.), due to strong coupling. (In a system governed
by a dimensional scale $\Lambda$, the nucleation rate is given by
$\Gamma \sim \Lambda^4 e^{-S}$, where $S$ is the action of the
Euclidean bounce solution interpolating between the false (normal
nuclear) vacuum and a bubble of critical size. At strong coupling,
$S$ is of order one, so there is no large exponential suppression
of the nucleation rate. The scale $\Lambda$ is of the order of
$\Delta$ or $\mu = 400$ MeV.)


Causality requires that the mechanical bounce of the core happen
over timescales larger than the light crossing time of the core,
or at least $10^{-4} s$. Hence, the phase transition occurs
instantaneously on astrophysical timescales. A nucleated bubble of
CSC phase expands relativistically -- liberated latent heat is
converted into its kinetic energy -- until it collides with other
bubbles. Because the system is strongly coupled, these collisions
lead to the rapid production of all of the low energy excitations
in the CSC phase, including (pseudo)Goldstone bosons and other
hadrons. The resulting release of energy resembles an explosion of
hadronic matter.

To estimate the total CSC energy released in the bounce, we use
the result that the ratio of CSC binding energy density to quark
energy density is of order $\left(\frac{ \Delta }{\mu }\right)^2$.
For $\mu \sim 400$ MeV, and $\Delta \sim$ 40-120 MeV, this ratio
is between .01 and .08, or probably a few percent.
\begin{equation}
E \sim \left(\frac{ \Delta }{\mu }\right)^2 M_{core} \ .
\end{equation}
In other words, the total energy release could be a few percent of
a solar mass, or $10^{52}$ ergs! This is significantly larger than
the energy usually attributed to the core shockwave, and possibly
of the order of the gravitational collapse energy $E_b$. The
implications for SN simulations are obviously quite intriguing.

In \cite{BH} it was suggested that strange matter formation might
overcome the energetic difficulties in producing type-II supernova
explosions. While there are strong arguments that the CSC
transition should be first order, and reasonable order of
magnitude estimates of the latent heat \cite{RW,Interface}, it is
not clear to us why there would be supercooling in a strange
matter transition. The conversion of up quarks to strange quarks
must proceed by the weak interactions, but the rate is still much
faster than any astrophysical timescale. Thus, the population of
strange quarks is likely to track its chemical equilibrium value
as the pressure of the core increases. There may be an important
effect on the nuclear equation of state from strangeness (e.g. a
softening of the pressure-density relationship), but we do not see
why there should be explosive behavior.

Our results are also relevant to hypernovae \cite{Hyper}, which
are observed to have ejecta kinetic energies 10-100 times larger
(of order $10^{52-53}$ ergs) than those of ordinary type II
supernovae. Accounting for this extra kinetic is extremely
challenging in standard scenarios. However, for exceptionally
massive stars with $M < 35  M_{\odot}$ (for $M > 35 M_{\odot}$ the
hydrogen envelope is lost during H-shell burning, and the core
size actually decreases~\cite{chlee}) there is a large core mass
which leads to a larger release of CSC binding energy. In fact,
the released energy might depend nonlinearly and sensitively on
the star's mass at the upper range. For example, the fraction of
$M_{core}$ which achieves critical density might be a sensitive
function of the mass of the star.

Another alternative is that hypernovae are the result of neutron
star mergers rather than the explosion of an individual star. This
possibility has been examined in relation to hypernovae as the
engines of gamma ray bursts (GRBs) \cite{GRB}. It seems quite
plausible that in the merger of two cold neutron stars a
significant fraction of the stars' mass undergoes the CSC
transition (i.e. crosses the critical pressure boundary for the
first time; in this case the temperature is probably negligible
relative to the CSC scale $\Delta$). This provides a substantial
new source of energy beyond gravitational binding, and may solve
the ``energy crisis'' problem for this model of GRBs \cite{GRB}.


Finally, we note that the trajectory of the SN core in the
temperature-density phase diagram might be rather complicated. The
parameters suggest a density transition (at $T < T_c$), but
subsequent reheating of the core due to the explosion, or to
neutrino diffusion \cite{Lattimer} might raise the temperature
above $T_c$, and lead to additional transitions across the
temperature boundary \cite{Carter}. When $T \sim T_c$, the Free
energies ($F = U - TS$) of the normal and CSC phases are
comparable, due to the larger entropy of the normal phase. The
evolution of bubbles in this regime is governed by relative Free
energies rather than energetics alone, and the transition is
presumably less explosive than the pressure transition at $T <<
T_c$.









\section{Summary}

Our understanding of QCD at high density has evolved dramatically
over the past two years, leading to remarkable progress in
understanding the QCD phase diagram and the color superconducting
state of quark matter. We have argued here that phase transitions
at high density and low temperature (relative to the QCD scale)
might play a role in violent astrophysical phenomena. In
particular, the densities and temperatures associated with core
collapse in massive stellar evolution suggest that CSC could play
an important role in type II supernovae, or possibly hypernovae
(GRBs).

Our assumptions concerning key parameters are conservative, and
taken from distinct (and heretofore independent) regimes of
inquiry: stellar astrophysics and dense quark matter. Yet, they
point to the interesting possibility that supernova explosions are
powered by CSC binding energy. It is well established that the
shockwave energy from core collapse is insufficient to produce an
explosion, and recent results incorporating Boltzman transport of
neutrinos show that neutrino reheating is also insufficient unless
non-spherical phenomena such as rotation or convection are taken
into account \cite{SNsim}. We are optimistic that future progress
in simulations will tell us much about whether and how latent
energy from QCD plays a role in stellar explosions.

\bigskip \noindent {\bf Acknowledgements}

\noindent I would like to thank F. Sannino and R. Ouyed for
organizing this very stimulating workshop. The work of S.H. was
supported in part under DOE contract DE-FG06-85ER40224 and by the
NSF through the Korean-USA Cooperative Science Program, 9982164.


%\newpage

\vskip 1 in \baselineskip=1.6pt

\begin{thebibliography}{99}

\bibitem{HHS} D.K. Hong, S.D.H. Hsu and F. Sannino, Phys. Lett. B516,
362,2001.


\bibitem{RW}  For a review see K.~Rajagopal, F.~Wilczek ;
S.D.H. Hsu, Proceedings of TMU - Yale Symposium on Dynamics of
Gauge Fields, ; D.~K.~Hong, ;
M.~Alford, .

\bibitem{SCF}  S.~C.~Frautschi, Asymptotic freedom and color
superconductivity in dense quark matter, in: Proceedings of the
Workshop on Hadronic Matter at Extreme Energy Density, N.~Cabibbo,
Editor, Erice, Italy (1978).

\bibitem{B}  B.~Barrois, Nucl.~Phys.~B, 129, 390 (1977).

\bibitem{BL}  D.~Bailin, and A.~Love, Phys.~Rept., 107, 325 (1984).

\bibitem{ARW}  M.~Alford, K.~Rajagopal, and F.~Wilczek, Phys.~Lett.~B, 422,
247 (1998).

\bibitem{RSSV}  R.~Rapp, T.~Sch\"afer, E.V.~Shuryak, M.~Velkovsky,
Phys.~Rev.~Lett., 81, 53 (1998).

\bibitem{ESH}  N.~Evans, S.~Hsu and M.~Schwetz, Nucl.~Phys.~B{\bf 551}, 275
(1999); Phys.~Lett.~B{\bf 449}, 281 (1999). S.~Hsu and M.~Schwetz,
Nucl.Phys.B{\bf 572}, 211 (2000).

\bibitem{EHSH} N. Evans, J. Hormuzdiar, S. Hsu and M. Scwhetz, Nucl.Phys.B{\bf
581}, 391 (2000).

\bibitem{Son}  D.T.~Son, Phys.~Rev.~D{\bf 59}, 094019 (1999).

\bibitem{SW}  T.~Sch\"{a}fer and F.~Wilczek, Phys.~Rev.~Lett.~{\bf 82} ,
3956 (1999).

\bibitem{DH}  D.~K.~Hong, Phys. Lett. B {\bf 473}, 118 (2000);
Nucl. Phys. B {\bf 582}, 451 (2000). D.~K.~Hong, T.~Lee and
D.~Min, Phys.~Lett.~B{\bf 477}, 137 (2000). D.~K.~Hong,
V.~A.~Miransky, I.~A.~Shovkovy and L.~C.~R.~Wiejewardhana,
Phys.~Rev.~D{\bf 61}, 056001 (2000).

\bibitem{PR}  R.D.~Pisarski and D.H.~Rischke, Phys.~Rev.~{\bf D61}, 051501
(2000);Phys.~Rev.~{\bf D61}, 074017 (2000).


\bibitem{CG}  R.~Casalbuoni and R.~Gatto, Phys.~Lett.~B{\bf 464}, 11 (1999);
Phys.~Lett.~B{\bf 469}, 213 (1999).

\bibitem{S}  F.~Sannino, Phys.~Lett.~B{\bf 480}, 280, (2000). S.~Hsu,
F.~Sannino and M.~Schwetz, .

\bibitem{CDS}  R.~Casalbuoni, Z.~Duan and F.~Sannino, Phys.~Rev.~D{\bf 62},
094004, (2000). R.~Casalbuoni, Z.~Duan and F.~Sannino,
. R.~Casalbuoni, Z.~Duan and F.~Sannino,
Phys.~Rev.~D63, 114026, 2001. R.~Ouyed, F.~Sannino 
to appear in Phys. Lett. B.

\bibitem{ChPT} S.D.H. Hsu and M. Schwetz,
Phys.Lett.B432:203,1998.

\bibitem{Universal} R. Pisarski and D. Rischke,
Phys.Rev.Lett.83:37,1999.

\bibitem{Interface} M. Alford, K. Rajagopal, S. Reddy and F. Wilczek,
.

\bibitem{SNcollapse}  S. L. Shapiro, and S.A. Teukolsky, Black Holes,
White Dwarfs, and Neutron Stars: the physics of compact objects
(Wiley-Interscience, 1983); B. Cooperstein and E. Baron, in
Supernovae, edited by A. Petschek (Springer-Verlag, 1990).

\bibitem{SNsim}  For recent reviews on the status of supernova simulations,
see H.T. Janka, K. Kifonidis and M. Rampp, ; A.
Mezzacappa, .

\bibitem{Hyper}  T. Galama et al, Nature 395, 670 (1998); Q. Daniel Wang,
ApJ 517, L27 (1999).

\bibitem{GRB}  Piran, T. 1999b, in Gamma-Ray Bursts: The first three
minutes, ASP Conference Series, Vol. 190, eds. J. Poutanen and R.
Svensson, p3.

\bibitem{Lattimer} J.M. Lattimer, in Supernova Shells and Their Birth Events,
edited by W. Kundt (Springer-Verlag, 1988). A. Burrows, T.J.
Mazurek and J.M. Lattimer, Ap. J. {\bf 251}, 325 (1981); A.
Burrows and J.M. Lattimer, Ap. J. {\bf 307}, 178 (1986).

\bibitem{Carter} G.W. Carter and S. Reddy, .


\bibitem{BH} O.G.~Benvenuto and J.~Horvath, Phys.~Rev.~Lett.{\bf
63}, 716 (1989).

\bibitem{chlee} G.~E.~Brown, C.~H.~Lee and H.~A.~Bethe,
New Astron.\  {\bf 4}, 313 (1999).








\end{thebibliography}

\newpage

\def\Discussion{
\setlength{\parskip}{0.3cm}\setlength{\parindent}{0.0cm}
     \bigskip\bigskip      {\Large {\bf Discussion}} \bigskip}
\def\speaker#1{{\bf #1:}\ }
\def\endDiscussion{}

\Discussion

\speaker{Student} How does one identify the order of a phase
transition? What is the difference between first and second order?

\speaker{Hsu} Formally speaking, the order refers to the level of
discontinuity exhibited by thermodynamic quantities (in the
infinite volume limit) at the transition. In an Nth order
transition, the Nth derivative of the free energy (with respect
to, e.g., temperature or density) is discontinuous, while lower
order derivatives are continuous. So, in a first order transition
the derivative of the free energy itself is discontinuous, whereas
in a second order transition the first derivative of the free
energy is continuous and the second derivative is not. One can
imagine even smoother transitions where only some high derivative
of the free energy is discontinuous. Alternatively, one can
distinguish by noting that the order parameter exhibits a jump in
the first order case, but evolves smoothly in the 2nd and higher
order cases. Of course, it is only first order transitions (with
nonzero latent heat) that can lead to the explosive phenomena
considered here.

\bigskip


\speaker{Mark Alford (University of Glasgow)}  (paraphrasing) In
our recent investigations of a model of the nuclear matter--quark
matter interface in neutron stars, we find that while the free
energy of the quark phase is lower, the energy density of the
quark core is actually higher than in the nuclear phase. If this
is the case, how does the phase transition proceed?


\speaker{Hsu} In this discussion I assumed that there is a high
density state (either the CSC or chirally restored phase) with
lower energy than the normal nuclear phase. Note that I really
mean energy here, not free energy: $F = H - \mu N = - P$. If we
consider the phase diagram as a function of baryon density (not
baryon chemical potential), the existence of a new phase implies a
state with lower energy than the normal nuclear phase at the same
density. If the transition is first order, there is the
possibility of latent energy release. In a core bounce, the baryon
density is increased on a timescale too short for bulk baryon
number flow. So, the relevant question is what happens when you
cross a boundary in the baryon density direction. The system will
try to tunnel into the lower energy phase, and the release of
latent energy powers the explosion.

On the other hand, if the compression is slow enough for
significant rearrangement of baryon number, then it is governed by
free energetics, or $F$. In the model you considered, the quark
matter phase has lower free energy, but higher energy and baryon
density than the nuclear phase. The system can only form the quark
phase by aggregating baryon number, and there is probably not
enough time for this to happen during a supernova core bounce. The
conversion from nuclear to quark phases also requires energy
input, which presumably has to come from gravitational collapse.







\endDiscussion

\end{document}



\section{Observations}

My observations began when I was called on an emergency basis to
aid two Albanian noblemen, Guglielmo B. and Ferrando di N., who
had ingested the poison in an apparent attempt to commit suicide.
(The facts of the case were rather mysterious, since these young
men had been seen in the harbor earlier in the day, affecting a
carefree attitude.  However, this need not concern us here.)  The
standard remedy would have been to induce vomiting; however, this
would have offended the sensibilities \index{sensibilities}  of
two noble ladies present, the daughters of the duke in whose
household this unfortunate incident took place.  Thus, I thought
it wise to begin by administering magneto-therapy using the
devices that I had at hand.

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%
%%   use this format to include an .eps figure into your paper
%%
\begin{figure}[htb]
\begin{center}
\epsfig{file=magnet.eps,height=1.5in} \caption{Plan of the magnet
used in the Mesmeric studies.} \label{fig:magnet}
\end{center}
\end{figure}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

The source of magnetic influence used in this therapy was the
simple iron pole magnet shown in Figure~\ref{fig:magnet}.  The
device, roughly the size of human palm, was held three fingers'
length away from each of the patients, while the blood cyanide
level was monitored.  As this operator proceeded in silence, only
a very small effect on the poison was observed.  In a moment,
however, the sisters began to weep and ululate in the most painful
manner.  Under this blanket of sound, the effects of the cyanide
reaction were observed to be completely reversed.  Indeed, even
before the end of a single session of treatment, the Albanians
were seen to wake from their comas and offer thanks for their
rescue.  Some sample blood cyanide levels are reported in
Table~\ref{tab:blood}.



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%
%%   use this format to include a LaTeX table  into your paper
%%
\begin{table}[b]
\begin{center}
\begin{tabular}{l|ccc}
Patient &  Initial level($\mu$g/cc) &  w. Magnet & w. Magnet and
Sound \\ \hline
 Guglielmo B.  &   0.12     &     0.10      &     0.001  \\
 Ferrando di N. &  0.15     &     0.11      &  $< 0.0005$ \\ \hline
\end{tabular}
\caption{Blood cyanide levels for the two patients.}
\label{tab:blood}
\end{center}
\end{table}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\section{Interpretation}

A description of these observations, obviously based on
second-hand reports, has already been reported by L. da
Ponte~\cite{daPonte}. \index{da Ponte}
 In his paper,
da Ponte attributes the mode of action of the magnetic force to
direct magnetic attraction.  This would be implausible with the
relatively small magnet used in this study; thus da Ponte was
forced to propose an absurdly large size for the  magnetic device.
His theory does not account for the acoustic influence, which is
known not to be present in magnetic action on inanimate objects.
A second hypothesis would state that the acoustic fields perturb
the chemistry of the cyanide reactions, allowing the reagents to
couple magnetically.   As with the first hypothesis, we are
skeptical, since again this phenomenon has no inorganic
counterpart.  Thus, we are forced to conclude, with Mesmer, that
magnetism couples directly to the life force, which, as is well
known, can be modulated by musical influences.

In conclusion, we note that the unusual mode of action of
magnetism is only one of many which have been reported in the
literature.  It is likely that there are further unusual magnetic
phenomena to be discovered.  All studies to date have involved
male patients.  But it is likely that the prevailing belief that
`women are all the same' would be overturned in a similar
investigation.




\bigskip
I am grateful to Don Alfonso d'Alba for certain services essential
to this investigation.



\begin{thebibliography}{99}

%%
%%  bibliographic items can be constructed using the LaTeX format in SPIRES:
%%    see    http://www.slac.stanford.edu/spires/hep/latex.html
%%  SPIRES will also supply the CITATION line information; please include it.
%%

\bibitem{Mesmer}
F. A. Mesmer, Proc. Wien. Acad. Sci. {\bf 13}, 1564, 1593 (1762).
%%CITATION = PWASA,13,1564;%%

\bibitem{diCenzo}
A. L. di Cenzo, Trans. Acad. Ducal.  Milan., {\bf 23}, 2647
(1771).
%%CITATION = TAADD,23,2647;%%

\bibitem{Muller}
For an exhaustive review of the Prussian literature, see A. D.
M\"uller, in Workshop on Action at a Distance, F. Eisenschmidt,
ed. (Springer, Berlin, 1774).

\bibitem{daPonte}

L. da Ponte, Trans. N. Y. Acad. Sci., {\bf 3}, 27 (1795).
%%CITATION = TNYAA,3,27;%%


%%%  useful macros for econf proceedings articles:

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%  basic formatting macros:
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%  single-line equations:

\def\beq{\begin{equation}}
\def\eeq#1{\label{#1}\end{equation}}
\def\eeqn{\end{equation}}

%%  multiple-line equations  (use \CR as the carriage return):

\def\beqa{\begin{eqnarray}}
\def\eeqa#1{\label{#1}\end{eqnarray}}
\def\eeqan{\end{eqnarray}}
\def\CR{\nonumber \\ }

%%  reference to an equation number:

\def\leqn#1{(\ref{#1})}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%  bibliographic entries   (use this or the SPIRES LaTeX output)

%%   Journal or bibliographic formatting macros are obsolete!
%%   The SLAC/SPIRES database supplies properly formatted citations;
%%      click on  LaTeX(US) or LaTeX(EU)   

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%  sizing and bars

\def\st{\scriptstyle}
\def\sst{\scriptscriptstyle}

\def\overbar#1{\overline{#1}}
\let\littlebar=\bar
\let\bar=\overbar

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%  text-mode macros:

\def\etal{{\it et al.}}
\def\ie{{\it i.e.}}
\def\eg{{\it e.g.}}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%  expectation values:

\def\VEV#1{\left\langle{ #1} \right\rangle}
\def\bra#1{\left\langle{ #1} \right|}
\def\ket#1{\left| {#1} \right\rangle}
\def\vev#1{\langle #1 \rangle}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%  caligraphic letters (for matrix elements, luminosity, etc.)

\def\D{{\cal D}}
\def\L{{\cal L}}
\def\M{{\cal M}}
\def\O{{\cal O}}
\def\W{{\cal W}}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%  matrix operations and fractions:

\def\One{{\bf 1}}
\def\hc{{\mbox{\rm h.c.}}}
\def\tr{{\mbox{\rm tr}}}
\def\half{\frac{1}{2}}
\def\thalf{\frac{3}{2}}
\def\third{\frac{1}{3}}
\def\tthird{\frac{2}{3}}

\def\Dslash{\not{\hbox{\kern-4pt $D$}}}
\def\dslash{\not{\hbox{\kern-2pt $\del$}}}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%  high-energy physics terminology:

\def\Pl{{\mbox{\scriptsize Pl}}}
\def\eff{{\mbox{\scriptsize eff}}}
\def\CM{{\mbox{\scriptsize CM}}}
\def\BR{\mbox{\rm BR}}
\def\ee{e^+e^-}
\def\sstw{\sin^2\theta_w}
\def\cstw{\cos^2\theta_w}
\def\mz{m_Z}
\def\gz{\Gamma_Z}
\def\mw{m_W}
\def\mt{m_t}
\def\gt{\Gamma_t}
\def\mh{m_h}
\def\gmu{G_\mu}
\def\GF{G_F}
\def\alphas{\alpha_s}
\def\msb{{\bar{\ssstyle M \kern -1pt S}}}
\def\lmsb{\Lambda_{\msb}}
\def\ELER{e^-_Le^+_R}
\def\EREL{e^-_Re^+_L}
\def\ELEL{e^-_Le^+_L}
\def\ERER{e^-_Re^+_R}
\def\eps{\epsilon}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%  supersymmetry:

\def\ch#1{\widetilde\chi^+_{#1}}
\def\chm#1{\widetilde\chi^-_{#1}}
\def\neu#1{\widetilde\chi^0_{#1}}
\def\s#1{\widetilde{#1}}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5

