It is shown that strange quark matter (SQM) objects, stars, and planets, can very efficiently convert the mechanical energy into hadronic energy when they oscillate. This is because the mass density at the edge of SQM objects, $\rho_0{=}4.7{\times}10^{14}\frac{\mathrm{g}}{\mathrm{cm}^3}$, is the critical density below which SQM is unstable with respect to decay into photons, hadrons, and leptons. We consider here radial oscillations of SQM objects that could be induced in stellar or planetary systems where tidal interactions are ubiquitous. Oscillations of $0.1\%$ radius amplitude already result in $1\,$keV per unit baryon number excitation near the surface of SQM stars. The excitation energy is converted into electromagnetic energy in a short time of 1 ms, during a few oscillations. Higher amplitude oscillations result in faster energy release that could lead to fragmentation or dissolution of SQM stars. This would have significant consequences for hypothetical SQM star binaries and planetary systems of SQM planets with regard to gravitational wave emission.

We study cylindrically symmetric steady-state accretion of polytropic test matter spiraling onto the symmetry axis in power-law and logarithmic potentials. The model allows one to qualitatively understand the accretion process in a symmetry different from that of the classical Bondi accretion. We study the integral curves as level lines of some Hamiltonian and also apply this method to Bondi accretion. The isothermal solutions in power-law potentials (as well as in any radius-dependent potential) can be expressed in exact form in terms of the Lambert W function, while in the case of logarithmic potential, exact solutions can be found for any polytropic exponent.

Observations suggest a slower stellar rotation relative to gas rotation in the outer part of the Milky Way Galaxy. This difference could be attributed to an interaction with the interstellar magnetic field. In a simple model, fields of order 10 micro Gauss are then required, consistently with the observed values. This coincidence suggests a tool for estimating magnetic fields in spiral galaxies. A North-South asymmetry in the rotation of gas in the Galaxy could be of magnetic origin too.

Within the disk model framework used to approximately describe flattened galaxies, we develop an iterative method of determining column mass density from rotation curve supplemented with isotropic velocity dispersion profile. This generalizes our previous iterative method to the case when the velocity dispersion becomes important. We show on the example of UGC 6446 galaxy, that taking the velocity dispersion into account results in some observational signatures in the behavior of the local mass-to-light ratio. Along with galactic magnetic fields, this is another factor allowing to substantially reduce the local mass-to-light ratio at galactic outskirts. Taking the velocity dispersion into account may also have some consequences for the division of mass distribution between various mass components in modeling rotation curves.

We study the vertical gradient in azimuthal velocity of spiral galaxy NGC 4244 in a thin disk model. With surface density accounting for the rotation curve, we model the gradient properties in the approximation of quasi-circular orbits and find the predictions to be consistent with the gradient properties inferred from measurements. This consistency may suggest that the mass distribution in this galaxy is flattened.

Ring mass density and the corresponding circular velocity in thin disk model are known to be integral transforms of one another. But it may be less familiar that the transforms can be reduced to one-fold integrals with identical weight functions. It may be of practical value that the integral for the surface density does not involve the velocity derivative, unlike the equivalent and widely known Toomre's formula.

Recently, we determined a lower bound for the Milky Way mass in a point mass approximation. This result was obtained for most general spherically symmetric phase-space distribution functions consistent with a measured radial velocity dispersion. As a stability test of these predictions against a perturbation of the point mass potential, in this paper we make use of a representative of these functions to set the initial conditions for a simulation in a more realistic potential of similar mass and accounting for other observations. The predicted radial velocity dispersion profile evolves to forms still consistent with the measured profile, proving structural stability of the point mass approximation and the reliability of the resulting mass estimate of $2.1\times10^{11}\mathrm{M}_{\odot}$ within $150\,\mathrm{kpc}$. We also find an interesting coincidence with the recent estimates based on the kinematics of the extended Orphan Stream. As a byproduct, we obtain the equations of motion in axial symmetry from a nonstandard Hamiltonian, and derive a formula in the spherical symmetry relating the radial velocity dispersion profile to a directly measured kinematical observable.

In this paper we model the vertical structure of mass distribution of the Milky Way galaxy in the framework of a finite-width global disk model. Assuming the Galactic rotation curve only, we test inside the solar orbit the predictions of the model for two measurable and unrelated to each other processes: the gravitational microlensing that allows to fix the disk width-scale by the best fit to measurements, and the vertical gradient of rotation modelled in the quasi-circular orbits approximation. The former is sensitive to the gravitating mass in compact objects and the latter is sensitive to all kinds of gravitating matter. The analysis points to a small width-scale of the considered disks and, at most, insignificant contribution of non-baryonic dark mater in the solar circle. The predicted high vertical gradient values in the rotation are consistent with the gradient measurements.

The magnetic fields that are present in the galaxy NGC 253 are exceptionally strong. This means that they can influence the rotation of matter and hence the mass-to-light ratio. In this context, we address the issue of the presence of a non-baryonic dark matter halo in this galaxy.

We model the phase-space of the kinematic tracers using general, smooth distribution functions to derive a conservative lower bound on the total mass within 150-200 kpc. By approximating the potential as Keplerian, the phase-space distribution can be simplified to that of a smooth distribution of energies and eccentricities. Our approach naturally allows for calculating moments of the distribution function, such as the radial profile of the orbital anisotropy. We construct a family of phase-spaces with the resulting radial velocity dispersion overlapping with that of distant kinematic tracers, while making no assumptions about the density of the tracers and the radial profile of the velocity anisotropy (beta). While there is no apparent upper bound for the Milky Way mass, at least as long as only the radial motions are concerned, we find a sharp lower bound for the mass that is small. In particular, a mass value of $2.4 \times 10^{11}$ of solar masses, is still consistent with the dispersion profile at larger radii. Compared with much greater mass values in the literature, this result shows that determining the Milky Way mass is strongly model dependent. We expect a similar reduction of mass estimates in models assuming more realistic mass profiles.

The aim of this work is to show, in the framework of a simple finite-width disk model, that the amount of mass seen through gravitational microlensing measurements in the region $0<R<R0$ is consistent with the dynamical mass ascertained from Galaxy rotation after subtracting gas contribution. Since microlensing only detects compact objects, this result suggests that a non-baryonic mass component may be negligible in this region.

We address the problem of the influence of large-scale magnetic fields on galactic rotation for the example of the spiral galaxy NGC 891. Based on its rotation curve and the surface density of HI we determine, in the framework of the global disc model, the surface density of matter. Then, based on the surface brightness, we determine the corresponding profile of the local mass-to-light ratio. We also model the vertical gradient of azimuthal velocity in the quasi-circular-orbit approximation, and compare it with measurements. We discuss what factors may influence the rotation of matter in NGC 891 and how this can translate to changes in the profile of the local mass-to-light ratio. In particular, we discuss the possible effect of magnetic fields on the motion of ionized gas, and, consequently, on the determination of the profile of the local mass-to-light ratio. Finally, we put forward the hypothesis that the asymmetry in magnetic fields observed in NGC 891 might be responsible for the observed anomalous behaviour of the vertical gradient.

We estimate the vertical gradient of rotational velocity for several spiral galaxies in the framework of a global thin-disc model, using the approximation of quasi-circular orbits. We obtain gradients having a broad range of values, in agreement with measurements, for galaxies with both low and high gradients. To model the gradient, it suffices to know the rotation curve only. We illustrate, using the example of galaxy NGC 4302 with particularly high gradients, that mass models of galactic rotation curves that assume a significant spheroidal mass component reduce the predicted gradient value, which may suggest that the mass distribution is dominated by a flattened disc-like component. We conclude that the value and behaviour of the vertical gradient in rotational velocity can be used to study the mass distribution in spiral galaxies.

In this paper, we aim to estimate the vertical gradients in the rotational velocity of the Galaxy. This is carried out in the framework of a global thin disc model approximation. The predicted gradient values coincide with the observed vertical fall-off in the rotation curve of the Galaxy. The gradient is estimated based on a statistical analysis of trajectories of test bodies in the gravitational field of the disc and in an analytical way using a quasi-circular orbit approximation. The agreement of the results with the gradient measurements is remarkable in view of other more complicated, non-gravitational mechanisms used for explaining the observed gradient values. Finally, we find that models with a significant spheroidal component give worse vertical gradient estimates than the simple disc model. In view of these results, we can surmise that, apart from the central spherical bulge and Galactic halo, the gross mass distribution in the Galaxy forms a flattened rather than spheroidal figure.

Spiral galaxies are studied using a simple global disc model as a means for approximate determination of mass profiles. Based on rotation curves and the amount of gas (HI+He), we find global surface mass densities consistent with measurements and compare them with B-band surface brightness profiles. As a result we obtain mass-to-light ratio profiles. We give some arguments for why our approach is reliable and sometimes better than those assuming ad hoc the presence of a massive non-baryonic dark matter halo. Using this model, we study galaxies NGC 7793, 1365, 6946 and UGC 6446. Based on a rotation curve from The HI Nearby Galaxy Survey (THINGS) we also study galaxy NGC 4536 and compare the results with those we published elsewhere for the same galaxy.

Following Wolszczan's landmark discovery of planets in orbit around pulsar PSR B1257+12 in 1991, over 300 planets in more than 200 planetary systems have been found. Therefore, the meaning of Wolszczan's discovery cannot be overestimated. In this paper we aim to convince the reader that the objects accompanying pulsar PSR B1257+12 are more exotic than thought so far. They might not be ordinary planets but dwarf strange quark stars, whereas the pulsar might be a quark star with standard mass, not a neutron star. If this was the case, it would indicate that strange quark matter is the ground state of matter.

Propagation of relativistic particles in the Schwarzschild gravitational field is studied. Particles emitted radially outward with speed at infinity exceeding c/2^0.5 are observed to be accelerated in the gravitational field by a distant observer. This is the Shapiro effect for relativistic particles. Slower particles are decelerated, as in Newtonian gravity. A speed-dependent potential for relativistic particles corresponding to the speed measured in terms of the coordinate time is derived to be $V=GM gamma(gamma^2 -2)/r which is repulsive for v>c/2^0.5. The gravitational repulsion could be revealed in satellite experiments with beams of relativistic particles subject to very precise time measurements. Principles of laboratory measurements able to test kinetic energy changes of relativistic particles in the Earth gravitational field are discussed.

Non-monotonic features of rotation curves, and also the related gravitational effects typical of thin disks -- like backward-reaction or amplification of rotation by negative surface density gradients -- which are characteristic imprints of disk-like mass distributions, are discussed in the axisymmetric thin disk model. The influence of the data cutoff in rotational velocity measurements on the determination of the mass distribution in flattened galaxies is studied. It has also been found that the baryonic matter distribution in the spiral galaxy NGC 5475, obtained in the axisymmetric thin disk approximation, accounts for the rotation curve of the galaxy. To obtain these results, the iteration method developed recently by the authors has been applied.

An iterative method for reconstructing mass distribution in spiral galaxies using a thin disk approximation is developed. As an example, the method is applied to galaxy NGC 4736; its rotation curve does not allow one to employ a model with a massive spherical halo. We find a global mass distribution in this galaxy (without non-baryonic dark matter) that agrees perfectly with the high resolution rotation curve of the galaxy. This mass distribution is consistent with the $I$-band luminosity profile with the mean mass-to-light ratio $M/L_I=1.2$, and also agrees with the amount of hydrogen observed in the outermost regions of the galaxy. We predict the total mass of the galaxy to be only $3.43\times10^{10}M_{\sun}$. It is very close to the value predicted by the modified gravity models and much less than the currently accepted value of $5.0\times10^{10}M_{\sun}$ (with $\approx70%$ of the mass in a dark matter halo).

Effects of the neutrino trapping and symmetry energy behavior are investigated in the framework of the chiral Kaplan-Nelson model with kaon condensation. Decrease in the condensation threshold during deleptonization if found to be generic regardless uncertainties in the nucleon-kaon interactions and symmetry energy. Quantitatively however, differences are shown to be important

Stationary, axisymmetric and asymptotically flat spacetimes of dust of which trajectories are integral curves of the time translation Killing vector are investigated. The flow has no Newtonian limit. Asymptotic flatness implies the existence of singularities of the curvature scalar that are distributions and that are not isolated from regularity regions of the flow. The singularities are closely related to the presence of additional stresses that contribute negative active mass to the total (Komar) mass, which is zero for asymptotically flat spacetimes. Several families of solutions are constructed and a multipole expansion is presented.

Our goal is to study the effects of the UV radiation from the first stars, quasars and decays of the hypothetical Super Heavy Dark Matter (SHDM) particles on the formation of primordial bound objects in the Universe. We trace the evolution of a spherically symmetric density perturbation in the Lambda Cold Dark Matter ($\Lambda$CDM) and MOND models, solving the frequency-dependent radiative transfer equation, non-equilibrium chemistry, and one-dimensional gas hydrodynamics. We concentrate on the destruction and formation processes of the H$_2$ molecule, which is the main coolant in the primordial objects.

Our goal is to study the effects of the UV radiation from the first stars, quasars and hypothetical Super Heavy Dark Matter (SHDM) particle decays on the formation of primordial bound objects in the Universe. We trace the evolution of a spherically symmetric density perturbation in the Lambda Cold Dark Matter and MOND model, solving the frequency-dependent radiative transfer equation, non-equilibrium chemistry, and one-dimensional gas hydrodynamics. We concentrate on the destruction and formation processes of the $H_{2}$ molecule, which is the main coolant in the primordial objects.

The age of the Universe in the $\Lambda$CDM cosmology with $\Omega_{matter}=0.26$ and $\Omega_{\Lambda}=0.74$ is the same as in the Milne cosmology which correspods to an almost empty universe. In both cases it is a reciprocal Hubble constant, $1/H_0$, that for now preferred value $H_0=71 km/s/Mpc$ is 13.7 billion years. The most curious coincidence is that at the present time, in the $\Lambda$CDM model the decelerated expansion is exactly compensated by the accelerated expansion, as if the Universe coast for 13.7 billion years.

We would like to discuss prospects for neutrino observations of the core-collapse supernova progenitor during neutrino-cooled stage. We will present new theoretical results on thermal neutrino and antineutrino spectra produced deep inside the pre-supernova core. Three competing processes: pair-, photo and plasma-neutrino production, are taken into account. The results will be used to estimate signal in existing and future neutrino detectors. Chance for supernova prediction is estimated, with possible aid to core-collapse neutrino and gravitational wave detectors in the form of early warning.

An new method of calculating the energy spectrum of neutrinos and antineutrinos produced in the electron-positron annihilation processes in hot stellar plasma is presented. Detection of these neutrinos, produced copiously in the presupernova which is evolutionary advanced neutrino-cooled star, may serve in future as a trigger of pre-collapse early warning system. Also, observation of neutrinos will probe final stages of thermonuclear burning in the presupernova. The spectra obtained with the new method are compared to Monte Carlo simulations. To achieve high accuracy in the energy range of interest, determined by neutrino detector thresholds, differential cross-section for production of the antineutrino, previously unknown in an explicit form, is calculated as a function of energy in the plasma rest frame. Neutrino spectrum is obtained as a 3-dimensional integral, computed with the use of the Cuhre algorithm of at least 5% accuracy. Formulae for the mean neutrino energy and its dispersion are given as a combination of Fermi-Dirac integrals. Also, useful analytical approximations of the whole spectrum are shown.

We find a class of exact solutions of differentially rotating dust in the framework of General Relativity. There exist asymptotically flat space-times of the flow with positive mass function that for radii sufficiently large is monotonic and tends to zero at infinity. Some of the space-times may have non-vanishing total angular momentum. The flow is essentially different from another exactly solvable flow described by van Stockum line element.

We study the localization of protons in the core of neutron stars for ten realistic nuclear models that share a common behaviour of nuclear symmetry energy which saturates and eventually decreases at high densities. This results in the low proton fraction of beta-stable neutron star matter. Protons form a small admixture in the neutron star core, which is localized at sufficiently high densities. For every model we calculate the density $n_{loc}$ above which the localization effect is present. Our results indicate that localization occurs at densities above $0.5-1.0 fm^{-3}$. The phase with localized protons occupies a spherical shell or a core region inside neutron stars which contains significant fraction of all nucleons. Proton localization is of great importance for astrophysical properties of neutron stars as it strongly affects transport coefficients of neutron star matter and can produce spontaneous magnetization in neutron stars.

Neutron stars are discussed as laboratories of physics of strong gravitational fields. The mass of a neutron star is split into matter energy and gravitational field energy contributions. The energy of the gravitational field of neutron stars is calculated with three different approaches which give the same result. It is found that up to one half of the gravitational mass of maximum mass neutron stars is comprised by the gravitational field energy. Results are shown for a number of realistic equations of state of neutron star matter.

We study the evolution of a spherically symmetric density perturbation in the Modified Newtonian Dynamics (MOND) model applied to the net acceleration over Hubble flow. The background cosmological model is a $\Lambda$-dominated, low-$\Omega_b$ Friedmann model with no Cold Dark Matter. We include thermal processes and non-equilibrium chemical evolution of the collapsing gas. We find that under these assumptions the first low-mass objects ($M \le 3\times 10^4 M_{\odot}$) may collapse already for $z\sim 30$, which is in quite good agreement with the recent WMAP results. A lower value of $a_0$ would lead to much slower collapse of such objects.

Physical mechanisms that can influence rotation curves of spiral galaxies are discussed. For dark matter studies, possible contributions due to magnetic fields and non-Newtonian gravitational accelerations should be carefully accounted for. We point out that magnetic fields are particularly important in outermost parts of the disk. In the framework of general relativity the physical reason of an enhanced gravity in spiral galaxies depends on the assumed metric. The additional gravity is provided for Schwarzschild metric by nonluminous mass, whereas for Vaidya metric [1] by emission of radiative energy. In the latter case the non-Newtonian acceleration displays 1/r behaviour. Also matter flows contribute to non-Newtonian gravity.

Analysis of the massive star properties during C, Ne, O and Si burning i.e. the neutrino-cooled stage, leads to the simplified neutrino emission model. In the framework of this model we have simulated spectrum of the antineutrinos. Flux normalized according to the massive star model with explicitly given neutrino luminosity allow us to predict signal produced in water Cherenkov detectors. The results are discussed from the point of view of the possibility of the core-collapse supernova event prediction in advance of a few days.

The signal produced in neutrino observatories by the pair-annihilation neutrinos emitted from a 20 $M_{\odot}$ pre-supernova star at the silicon burning phase is estimated. The spectrum of the neutrinos with an average energy $\sim$2 MeV is calculated with the use of the Monte Carlo method. A few relevant reactions for neutrinos and anti-neutrinos in modern detectors are considered. The most promising results are from $\bar{\nu}_e + p \longrightarrow n + e^{+}$ reaction. During the Si-burning phase we expect 1.27 neutrons/day/kiloton of water to be produced by neutrinos from a star located at a distance of 1 kpc. Small admixture of effective neutron-absorbers as e.g. NaCl or GdCl$_{3}$ makes these neutrons easily visible because of Cherenkov light produced by electrons which were hit by $\sim$8 MeV photon cascade emitted by Cl or Gd nuclei. The estimated rate of neutron production for SNO and Super-Kamiokande is, respectively, 2.2 and 41 events per day for a star at 1 kpc. For future detectors UNO and Hyper-Kamiokande we expect 5.6 and 6.9 events per day even for a star 10 kpc away. This would make it possible to foresee a massive star death a few days before its core collapse. Importance of such a detection for theoretical astrophysics is discussed.

An attempt is made to assess the significance of rotation in the core-collapse supernova phenomenon, from both observational and theoretical point of view. The data on supernovae particularly indicative of the role of rotation in the collapse-triggered explosion is emphasized. The problem of including the rotation of presupernova core into the supernova theory is considered. A two-dimensional classification scheme of core-collapse supernovae is proposed which unifies 'classical' supernovae of type Ib/c and type II, 'hypernovae' and some GRB events.

Einstein's general relativity predicts that pressure, in general stresses, play a similar role to energy density in generating gravity. The source of gravitational field, the active gravitational mass density, sometimes referred to as Whittaker's mass density, is not conserved, hence its changes can propagate as monopole gravitational waves. Such waves can be generated only by astrophysical sources with varying gravitational mass. Here we show that relativistic fireballs, considered in modelling gamma-ray burst phenomena, are likely to radiate monopole gravitational waves from high-pressure plasma with varying Whittaker's mass. Also, ejection of a significant amount of initial mass-energy of the progenitor contributes to the monopole gravitational radiation. We identify monopole waves with h^11+h^22 waves of Eddington's classification which propagate (in the z-direction) together with the energy carried by massless fields. We show that the monopole waves satisfy Einstein's equations, with a common stress-energy tensor for massless fields. The polarization mode of monopole waves is Phi_22, i.e. these are perpendicular waves which induce changes of the radius of a circle of test particles only (breathing mode). The astrophysical importance of monopole gravitational waves is discussed.

The cooling of neutron stars by URCA processes in the kaon-condensed neutron star matter for various forms of nuclear symmetry energy is investigated. The kaon-nucleon interactions are described by a chiral lagrangian. Nuclear matter energy is parametrized in terms of the isoscalar contribution and the nuclear symmetry energy in the isovector sector. High density behaviour of nuclear symmetry energy plays an essential role in determining the composition of the kaon-condensed neutron star matter which in turn affects the cooling properties. We find that the symmetry energy which decreases at higher densities makes the kaon-condensed neutron star matter fully protonized. This effect inhibits strongly direct URCA processes resulting in slower cooling of neutron stars as only kaon-induced URCA cycles are present. In contrast, for increasing symmetry energy direct URCA processes are allowed in the almost whole density range where the kaon condensation exists.

We calculate the evolution of a low-mass ($M \le 10^5 M_{\odot}$) spherically symmetric density perturbation in the $\Omega_b h^2=0.02$, $\Omega_M =0.35$, $\Omega_{\Lambda}=0.65$, $h=0.72$ Universe. The results are compared with the ones that assume no cosmological constant and the flat, dark matter dominated Universe. We include thermal processes and non-equilibrium chemical evolution of the collapsing gas. We find that direct formation of bound objects with such masses by $z=8$ is unlikely so in fact they may form only through fragmentation of greater objects. This is in stark contrast to the $\Omega=1$ pure CDM cosmology, where low-mass objects form abundantly at redshits $z>10$.

We study properties of the proton component of neutron star matter for a number of realistic nuclear models. Protons which form a few percent admixture tend to be localized in potential wells corresponding to neutron matter inhomogeneities created by the protons in the neutron medium. We calculate the energy of the Wigner-Seitz cell enclosing a single localized proton. The neutron background is treated in the Thomas-Fermi approximation and the localized proton is described by the Gaussian wave function. The neutron density profile is obtained by solving the appropriate variational equation. This approach gives lower energies of localized protons than obtained previously with less sophisticated methods.

We trace the evolution of a spherically symmetric density perturbation in the MOdified Newtonian Dynamics (MOND) model. The background cosmological model is a $\Lambda$-dominated, low-$\Omega_b$ Friedmann model with no Cold Dark Matter. We include thermal processes and non-equilibrium chemical evolution of the collapsing gas. We find that the first density perturbations which collapse to form luminous objects have mass $\sim 10^5 M_{\odot}$. The time of the final collapse of these objects depends mainly on the value of the MOND acceleration $a_0$ and also on the baryon density $\Omega_b$. For the "standard" value $a_0=1.2\times 10^{-8}$ cm/s$^2$ the collapse starts at redshift $z \sim 160$ for $\Omega_b=0.05$ and $z \sim 110$ for $\Omega_b=0.02$.

The influence of the nuclear symmetry energy on the formation of a mixed quark-nucleon phase in neutron star cores is studied. We use simple parametrizations of the nuclear matter equation of state, and the bag model for the quark phase. The behavior of nucleon matter isobars, which is responsible for the existence of the mixed phase, is investigated. The role of the nuclear symmetry energy changes with the value of the bag constant B. For lower values of B the properties of the mixed phase do not depend strongly on the symmetry energy. For larger B we find that a critical pressure for the first quark droplets to form is strongly dependent on the nuclear symmetry energy, but the pressure at which last nucleons disappear is independent of it.

A model of the ferromagnetic origin of magnetic fields of neutron stars is considered. In this model, the magnetic phase transition occurs inside the core of neutron stars soon after formation. However, owing to the high electrical conductivity the core magnetic field is initially fully screened. We study how this magnetic field emerges for an outside observer. After some time, the induced field that screens the ferromagnetic field decays enough to uncover a detectable fraction of the ferromagnetic field. We conjecture that weak fields of millisecond pulsars of 10^8-10^9 G could be identified with ferromagnetic fields of unshielded fraction f=10^-4 resulting from the decay of screening fields by a factor 1-f in 10^8 yr since their birth.

We study the influence of a high density behaviour of the nuclear symmetry energy on a kaon condensation in neutron stars. We find that the symmetry energy typical for several realistic nuclear potentials, which decreases at high densities, inhibits kaon condensation for weaker kaon-nucleon couplings. There exists a threshold coupling above which the kaon condensate forms at densities exceeding some critical value. This is in contrast to the case of rising symmetry energy, as e.g. for relativistic mean field models, when the kaon condensate can form for any coupling at a sufficiently high density. Properties of the condensate are also different in both cases.

We study the equation of state of beta-stable dense matter and models of neutron stars in the relativistic mean field theory with the isovector scalar mean field corresponding to the delta-meson [a_0(980)]. A range of values of the delta-meson coupling compatible with the Bonn potentials is explored. Parameters of the model in the isovector sector are constrained to fit the nuclear symmetry energy, E_s=30 MeV. We find that the quantity most sensitive to the delta-meson coupling is the proton fraction of neutron star matter. It increases significantly in the presence of the delta-field. The energy per baryon also increases but the effect is smaller. The equation of state becomes slightly stiffer and the maximum neutron star mass increases for stronger delta-meson coupling.

A short introduction is given to astrophysics of neutron stars and to physics of dense matter in neutron stars. Observed properties of astrophysical objects containing neutron stars are discussed. Current scenarios regarding formation and evolution of neutron stars in those objects are presented. Physical principles governing the internal structure of neutron stars are considered with special emphasis on the possible spin ordering in the neutron star matter.

Relativistic mean field (RMF) theory of nuclear matter with the isovector scalar mean field corresponding to the delta-meson [a_0(980)] is studied. While the delta-meson mean field vanishes in symmetric nuclear matter, it can influence properties of asymmetric nuclear matter in neutron stars. The RMF contribution due to delta-field to the nuclear symmetry energy is negative. To fit the empirical value, E_s=30 MeV, a stronger rho-meson coupling is required than in the absence of the delta-field. The energy per particle of neutron matter is then larger at high densities than the one with no delta-field included. Also, the proton fraction of beta-stable matter increases. Splitting of proton and neutron effective masses due to the delta-field can affect transport properties of neutron star matter.

The equation of state of hadronic matter in neutron stars is briefly reviewed. Uncertainties regarding the stiffness and composition of hadronic matter are discussed. Importance of poorly known short range interactions of nucleons and hyperons is emphasized. Condensation of meson fields and the role of subhadronic degrees of freedom is considered. Empirical constraints on the equation of state emerging from observations of neutron stars are discussed. Implications of possible black hole remnant of SN1987A for the equation of state are also discussed.

Recently, Cottingham, Kalafatis and Vinh Mau proposed, in the framework of the Lee-Wick model, a scenario of formation of quark stars in the cosmological quark-hadron phase transition with high degree of supercooling. We show here that the cosmological quark-hadron phase transition in the Lee-Wick model with such a high degree of supercooling cannot be completed. No quark stars could be formed in this scenario.