Hoyle and Folwler showed that there could be Radiation Pressure Supported Stars (RPSS) even in Newtonian gravity. Much later, Mitra found that one could also conceive of their General Relativistic (GR) version, "Relativistic Radiation Pressure Supported Stars" (RRPSSs). While RPSSs have $z\ll 1$, RRPSSs have $z \gg 1$, where $z$ is the surface gravitational redshift. Here we elaborate on the formation of RRPSSs during continued gravitational collapse by recalling that a contracting massive star must start trapping radiation as it would enter its {\em photon sphere}. It is found that, irrespective of the details of the contraction process, the trapped radiation flux should attain the corresponding Eddington value at sufficiently large $z\gg 1$. This means that continued GR collapse may generate an intermediate RRPSS with $z\gg 1$ before a true BH state with $z=\infty$ is formed asymptotically. An exciting consequence of this is that the stellar mass black hole candidates, at present epoch, should be hot balls of quark gluon plasma, as has been discussed by Royzen in a recent article entitled "{\it QCD against black holes?}".

A neutron star in binary orbit with a low-mass non-degenerate companion becomes a source of x-rays with millisecond variability when mass accretion spins it up. Centrifugally driven changes in density profile may initiate a phase transition in a growing region of the core parallel to what may take place in an isolated millisecond pulsar, but in reverse. Such a star will spend a longer time in the spin frequency range over which the transition occurs than elsewhere because the change of phase, paced by the spinup rate, is accompanied by a growth in the moment of inertia. The population of accreters will exhibit a clustering in the critical frequency range. A phase change triggered by changing spin and the accompanying adjustment of moment of inertia has its analogue in rotating nuclei.

Theoretically, the phase transition between the confined and deconfined phases of quarks can have a remarkable effect on the spin properties of millisecond pulsars and on the spin distribution of the population of x-ray neutron stars in low-mass binaries. In the latter class of stars, the effect has already been observed---a strong clustering in the population in a narrow band of spins. The observed clustering cannot presently be uniquely assigned to the phase transition as cause. However, there is another possible signal---not so far observed---in millisecond pulsars that we also discuss. It would have the same origin, and its discovery would tend to confirm the interpretation in terms of a phase transition in the stellar core and that the high density phase exists in the higer mass canonical neutron stars.

We model the evolution of canonical pulsars from the death line to millisecond pulsars through the X-ray neutron star stage of accretion from a low-mass companion. We trace this evolution in magnetic field strength starting at $B=10^{12}$ to $10^8$ G and in period of about 1 s to milliseconds. Important factors are accretion rate and the decay rate of the magnetic field. A broad swathe is traced in the $B-P$ plane according to the value of these factors, which represent different conditions of the binary pair. An important ingredient is the dependence of the stellar moment of inertia on rotation frequency (time).

Uncertainty of the hyperon couplings, in particular, that of the $\Sigma^-$, in dense matter raises the question of the behavior of the electrochemical potential in neutron star matter, which is crucial to the possible presence of the kaon condensed phase. We show that regardless of this uncertainty, the $\Lambda$ hyperon, whose coupling can be constrained by its binding in nuclear matter and other observations, also aided by the $\Xi^-$, introduce a saturation of the electrochemical potential just as the $\Sigma^-$ would otherwise do, which tends to mitigate against kaon condensation. The maximum possible mass of neutron stars appears to be $\sim 1.5 \msun$ independent of the uncertainties.

We study the spin evolution of X-ray neutron stars in binary systems, which are being spun up by mass transfer from accretion disks. Our investigation reveals that a quark phase transition resulting from the changing central density induced by the changing spin, can lead to a pronounced peak in the frequency distribution of X-ray neutron stars. This finding provides one of several possible explanations available in the literature, or at least a contributor to part of the observed anomalous frequency distribution of neutron stars in low-mass X-ray binaries (LMXBs), which lie in a narrow band centered at about 300 Hz, as found by the Rossi Explorer (RXTE).

We give a general thermodynamical argument showing that in neutron stars, the Coulomb structured mixed phase is always favored for any first order phase transition involving systems in equilibrium with baryon number and electric charge as the two independent components. This finding is likely to have important consequences for many neutron star properties, e.g., glitch phenomena, transport and superfluid properties, r-mode instabilities, and the braking index.

We calculate for the first time the surface tension and curvature coefficient of a first order phase transition between two possible phases of cold nuclear matter, a normal nuclear matter phase in equilibrium with a kaon condensed phase, at densities a few times the saturation density. We find the surface tension is proportional to the difference in energy density between the two phases squared. Furthermore, we show the consequences for the geometrical structures of the mixed phase region in a neutron star.

While model-independent limits are always interesting, a limit on neutron star radius as a function of mass attains special interest in light of recent interpretations of the periodic as well as quasi-periodic oscillations (QPOs) in brightness of X-rays emitted from neutron stars that are accreting matter from a low-mass companion. Here we derive such a limit based only on well accepted principles. We discuss our limit in connection with a recent interpretation of X-ray pulsations from SAX J1808.4-3658 as indicating a strange-star candidate, and show that this object can also be a normal neutron star, though one whose central core has very high density. The most plausible high-density phase of hadronic matter, which is also expected to be very compressible, is quark matter. So an alternative to the strange star interpretation of SAX J1808.4-3658 is that it is a hybrid neutron star.

The work of J. A. Wheeler in the mid 1960's showed that for smooth equations of state no stable stellar configurations with central densities above that corresponding to the limiting mass of ``neutron stars'' (in the generic sense) were stable against acoustical vibrational modes. A perturbation would cause any such star to collapse to a black hole or explode. Accordingly, there has been no reason to expect that a stable degenerate family of stars with higher density than the known white dwarfs and neutron stars might exist. We have found a class of exceptions corresponding to certain equations of state that describe a first order phase transition. We discuss how such a higher density family of stars could be formed in nature, and how the promising new exploration of oscillations in the X-ray brightness of accreting neutron stars might provide a means of identifying them. Our proof of the possible existence of a third family of degenerate stars is one of principle and rests on general principles like causality, microstability of matter and General Relativity.

First order Bose condensation in asymmetric nuclear matter and in neutron stars is studied, with particular reference to kaon condensation. We demonstrate explicitly why the Maxwell construction fails to assure equilibrium in multicomponent substances. Gibbs conditions and conservation laws require that for phase equilibrium, the charge density must have opposite sign in the two phases of isospin asymmetric nuclear matter. The mixed phase will therefore form a Coulomb lattice with the rare phase occupying lattice sites in the dominant phase. Moreover, the kaon condensed phase differs from the normal phase, not by the mere presence of kaons in the first, but also by a difference in the nucleon effective masses. The mixed phase region, which occupies a large radial extent amounting to some kilometers in our model neutron stars, is thus highly heterogeneous. It should be particularly interesting in connection with the pulsar glitch phenomenon as well as transport properties.

We discuss the nature of the kaon condensation phase transition. We find several features which, if kaons condense in neutron stars, are not only remarkable, but must surely effect such properties as superfluidity and transport properties, which in turn are relevant to the glitch phenomenon and cooling rates of neutron stars. The mixed phase, because of the extensive pressure range that it spans, will occupy a broad radial extent in a neutron star. This region is permeated with microscopic drops (and other configurations) located at lattice sites of one phase immersed in the background of the other phase. The electric charge on drops is opposite to that of the background phase {\sl and} nucleons have a mass approximately a factor two different depending on whether they are in the drops or the background phase. A large part of the stellar interior has this highly non-homogeneous structure.

We study the influence of a possible H dibaryon condensate on the equation of state and the overall properties of neutron stars whose population otherwise contains nucleons and hyperons. In particular, we are interested in the question of whether neutron stars and their masses can be used to say anything about the existence and properties of the H dibaryon. We find that the equation of state is softened by the appearance of a dibaryon condensate and can result in a mass plateau for neutron stars. If the limiting neutron star mass is about that of the Hulse-Taylor pulsar a condensate of H dibaryons of vacuum mass 2.2 GeV and a moderately attractive potential in the medium could not be ruled out. On the other hand, if the medium potential were even moderately repulsive, the H, would not likely exist in neutron stars. If neutron stars of about 1.6 solar mass were known to exist, attractive medium effects for the H could be ruled out. Certain ranges of dibaryon mass and potential can be excluded by the mass of the Hulse-Taylor pulsar which we illustrate graphically.

A solitary millisecond pulsar, if near the mass limit, and undergoing a phase transition, either first or second order, provided the transition is to a substantially more compressible phase, will emit a blatantly obvious signal---spontaneous spin-up. Normally a pulsar spins down by angular momentum loss to radiation. The signal is trivial to detect and is estimated to be ``on'' for 1/50 of the spin-down era of millisecond pulsars. Presently about 25 solitary millisecond pulsars are known. The phenomenon is analogous to ``backbending'' observed in high spin nuclei in the 1970's.

A phase transition in the nature of matter in the core of a neutron star, such as quark deconfinement or Bose condensation, can cause the spontaneous spin-up of a solitary millisecond pulsar. The spin-up epoch for our model lasts for $2\times 10^7$ years or 1/50 of the spin-down time.

A phase transition in the nature of matter in the core of a neutron star, such as quark deconfinement or Bose condensation, can cause the spontaneous spin-up of a solitary millisecond pulsar. The spin-up epoch for our model lasts for $2\times 10^7$ years or 1/50 of the spin-down time (Glendenning, Pei and Weber in Ref. \cite{glen97:a}). The possibility exists also for future measurements on X-ray neutron stars with low-mass companions for mapping out the tell-tale ``backbending'' behavior of the moment of inertia. Properties of phase transitions in substances such as neutron star matter, which have more than one conserved charge, are reviewed.

The equation of state and the properties of neutron stars are studied for a phase transition to a charged kaon condensate. We study the mixed phase by using Gibbs condition with comparison to the hitherto applied Maxwell construction. Implications for kaon condensation and for the mass-radius relation of condensed neutron stars are examined.

Phase transitions in compact stars are discussed including hyperonization, deconfinement and crystalline phases. Reasons why kaon condensation is unlikely is reviewed. Particular emphasis is placed on the evolution of internal structure with spin-down of pulsars. A signature of a first order phase transition in the timing structure of pulsars which is strong and easy to measure, is identified.

In the first of these two lectures I will discuss the rich constitution of neutron stars as a consequence of the Pauli principle which is engaged by the dominance of gravity over the nuclear force. Three especially interesting phenomena are discussed in this contect--(1) a mechanism for the formation of low-mass black holes distinct in their mass-range from the black holes formed in the prompt collapse of an entire star, (2) a multilayered crystalline structure consisting of confined hadronic matter embedded in a background of deconfined quark matter (or vice versa) which occupies a many kilometer thick inner region, and (3) a clean and pronounced signal of the formation of quark matter in the interior of neutron stars. In the second lecture I will discuss the strange matter hypothesis, its viability as well as its consequences for compact stars and a new family of white dwarfs with dense nuclear matter central regions some orders of magnuitude greater than in ordinary white dwarfs.

We review the differences in first order phase transition of single and multi-component systems, and then discuss the crystalline structure expected to exist in the mixed confined deconfined phase of hadronic matter. The particular context of neutron stars is chosen for illustration. The qualitative results are general and apply for example to the vapor-liquid transition in subsaturated asymmetric nuclear matter.

We demonstrate that the form and location (not the size or spacing) of the energetically preferred geometrical structure of the crystalline quark-hadron mixed phase in a neutron star is very sensitive to finite size terms beyond the surface term. We consider two independent approaches of including further finite size terms, namely the multiple reflection expansion of the bag model and an effective surface tension description. Thus care should be taken in any model requiring detailed knowledge of these crystalline structures.

The conversion of nuclear matter to quark matter in the core of a rotating neutron star alters its moment of inertia. Hence the epoch over which conversion takes place will be signaled in the spin-down "signal_prl.tex" 581 lines, 22203 characters characteristics of pulsars. We find that an observable called the braking index should be easily measurable during the transition epoch and can have a value far removed (by orders of magnitude) from the canonical value of three expected for magnetic dipole radiation, and may have either sign. The duration of the transition epoch is governed by the slow loss of angular momentum to radiation and is further prolonged by the reduction in the moment of inertia caused by the phase change which can even introduce an era of spin-up. We estimate that about one in a hundred pulsars may be passing through this phase. The phenomenon is analogous to ``bachbending'' observed in the moment of inertia of rotating nuclei observed in the 1970's, which also signaled a change in internal structure with changing spin.

For the past 20 years it had been thought that the coexistence phase of the confined hadronic and quark matter phases, assumed to be a first order transition, was strictly excluded from neutron stars. This, however, was due to a seemingly innocuous idealization which has approximated away important physics. The reason is that neutron stars constitute multi-component bodies rather than single-component ones formerly (and incorrectly) used to describe the deconfinement phase transition in neutron stars. So, contrary to earlier claims, `neutron' stars may very well contain quark matter in their cores surrounded by a mixed-phase region of quark and hadronic matter. Such objects are called hybrid stars. The structure of such stars as well as an observable signature that could signal the existence of quark matter in their cores are discussed in this paper.

After giving an overview of the history and idea of neutron stars, I shall introduce, in part one of my lectures, a representative collection of models for the equation of state of superdense neutron star matter, which account for various physical possibilities concerning the unknown behavior of such matter. Part two of my lectures deal with the construction of models of static as well as rapidly rotating neutron stars, and an investigation of the cooling behavior of such objects. Both is performed in the framework of Einstein's theory of general relativity. For this purpose the broad collection of models for the equation of state, derived in part one, will be used as an input. Various computed neutron star properties (e.g., masses, radii, moments of inertia, red and blueshifts, limiting rotational periods, cooling behavior) will be compared with the body of observed pulsar data. Finally, part three deals with an investigation of the properties of stars made up of absolutely stable strange quark matter. Their properties are compared with those of their nonstrange counterparts, i.e., neutron stars.

This paper gives an overview of the properties of all possible equilibrium sequences of compact strange-matter stars with nuclear crusts, which range from strange stars to strange dwarfs. In contrast to their non-strange counterparts, --neutron stars and white dwarfs--, their properties are determined by two (rather than one) parameters, the central star density and the density at the base of the nuclear crust. This leads to stellar strange-matter configurations whose properties are much more complex than those of the conventional sequence. As an example, two generically different categories of stable strange dwarfs are found, which could be the observed white dwarfs. Furthermore we find very-low-mass strange stellar objects, with masses as small as those of Jupiter or even lighter planets. Such objects, if abundant enough in our Galaxy, should be seen by the presently performed gravitational microlensing searches.

This paper gives an overview of the properties of all possible equilibrium sequences of compact strange-matter stars with nuclear crusts, which range from strange stars to strange dwarfs. In contrast to their non-strange counterparts, --neutron stars and white dwarfs--, their properties are determined by two (rather than one) parameters, the central star density and the density at the base of the nuclear crust. This leads to stellar strange-matter configurations whose properties are much more complex than those of the conventional sequence. As an example, two generically different categories of stable strange dwarfs are found, which could be the observed white dwarfs. Furthermore we find very low-mass strange stellar objects, with masses as small as those of Jupiter or even lighter planets. Such objects, if abundant enough in our Galaxy, should be seen by the presently performed gravitational microlensing searches. Further aspects studied in this paper concern the limiting rotational periods and the cooling behavior of neutron stars and their strange counterparts.

A collection of modern, field-theoretical equations of state is applied to the investigation of cooling properties of compact stars. These comprise neutron stars as well as hypothetical strange matter stars, made up of absolutely stable 3-flavor strange quark matter. Various uncertainties in the behavior of matter at supernuclear densities, e.g., hyperonic degrees of freedom, behavior of coupling strengths in matter, pion and meson condensation, superfluidity, transition to quark matter, absolute stability of strange quark matter, and last but not least the many-body technique itself are tested against the body of observed cooling data.