We integrate numerically axially symmetric stationary Einstein equations describing self-gravitating disks around spinless black holes. The numerical scheme is based on a method developed by Shibata, but contains important new ingredients. We derive a new general-relativistic Keplerian rotation law for self-gravitating disks around spinning black holes. Former results concerning rotation around spin-less black holes emerge in the limit of a vanishing spin parameter. These rotation curves might be used for the description of rotating stars, after appropriate modification around the symmetry axis. They can be applied to the description of compact torus--black hole configurations, including active galactic nuclei or products of coalescences of two neutron stars.

We show that the solutions obtained in the paper `An exact solution for arbitrarily rotating gaseous polytropes with index unity' by Kong, Zhang, and Schubert represent only approximate solutions of the free-boundary Euler-Poisson system of equations describing uniformly rotating, self-gravitating polytropes with index unity. We discuss the quality of such solutions as approximations to the rigidly rotating equilibrium polytropic configurations.

Considering probability distribution as a function of the average density $\bar{\rho}$ computed for 424 extrasolar planets we identify three log-normal Gaussian population components. The two most populous components at $\bar{\rho}\simeq0.7$ g/cc and $\bar{\rho}\simeq7$ g/cc are the ice/gas giants and iron/rock super-Earths, respectively. A third component at $\bar{\rho}\simeq30$ g/cc is consistent with brown dwarfs, i.e., electron degeneracy supported objects. We note presence of several extreme density planetary objects.

Goal of the presented research is to construct simplified model of the core-halo structures in binary systems. Examples are provided by Thorne-Zytkov objects, hot Jupiters, protoplanets with large moons, red supergiants in binaries and globular clusters with central black hole. Instability criteria due to resonance between internal and orbital frequencies in such a systems has been derived. To achieve assumed goals, generalized planar circular restricted three body problem is investigated with one of the point masses, $M$, replaced with spherical body of finite size. Mechanical system under consideration includes two large masses $m$ and $M$ and the test body with small mass $\mu$. Only gravitational interactions are considered. Equations of motion are presented, and linear instability criteria are derived using quantifier elimination. Motion of the test mass $\mu$ is shown to be unstable due to resonance between orbital and internal frequencies if $\frac{M}{d^3} < \frac{4}{3} \pi \rho < \frac{ M + 3 m \left( 1+\mu/M \right)^{-1}}{d^3}$, where $\rho$ is the central density of mass $M$, and $d$ distance between masses $m$ and $M$ (circular orbit diameter). The above result is important for core-collapse supernova theory, with mass $\mu$ identified with helium core of the exploding massive star. The instability cause off-center supernova "ignition" relative to the center-of-mass of the hydrogen envelope. The instability is also inevitable during protoplanet growth, with hypothetical ejection of the rocky core from gas giants and formation of the "puffy planets" due to resonance with orbital frequency. Hypothetical central intermediate black holes of the globular clusters are also in unstable position with respect to perturbations caused by the Galaxy.

We study the evolution of the collapsing core of a 15 Msun blue supergiant supernova progenitor from the core bounce until 1.5 seconds later. We present a sample of hydrodynamic models parameterized to match the explosion energetics of SN 1987A. We find the spatial model dimensionality to be an important contributing factor in the explosion process. Compared to two-dimensional simulations, our three-dimensional models require lower neutrino luminosities to produce equally energetic explosions. We estimate that the convective engine in our models is 4% more efficient in three dimensions than in two dimensions. We propose that the greater efficiency of the convective engine found in three-dimensional simulations might be due to the larger surface-to-volume ratio of convective plumes, which aids in distributing energy deposited by neutrinos. We do not find evidence of the standing accretion shock instability nor turbulence being a key factor in powering the explosion in our models. Instead, the analysis of the energy transport in the post-shock region reveals characteristics of penetrative convection. The explosion energy decreases dramatically once the resolution is inadequate to capture the morphology of convection on large scales. This shows that the role of dimensionality is secondary to correctly accounting for the basic physics of the explosion. We also analyze information provided by particle tracers embedded in the flow, and find that the unbound material has relatively long residency times in two-dimensional models, while in three dimensions a significant fraction of the explosion energy is carried by particles with relatively short residency times.

Aims: We present neutrino light curves and energy spectra for two representative type Ia supernova explosion models: a pure deflagration and a delayed detonation. Methods: We calculate the neutrino flux from $\beta$ processes using nuclear statistical equilibrium abundances convoluted with approximate neutrino spectra of the individual nuclei and the thermal neutrino spectrum (pair+plasma). Results: Although the two considered thermonuclear supernova explosion scenarios are expected to produce almost identical electromagnetic output, their neutrino signatures appear vastly different, which allow an unambiguous identification of the explosion mechanism: a pure deflagration produces a single peak in the neutrino light curve, while the addition of the second maximum characterizes a delayed-detonation. We identified the following main contributors to the neutrino signal: (1) weak electron neutrino emission from electron captures (in particular on the protons Co55 and Ni56) and numerous beta-active nuclei produced by the thermonuclear flame and/or detonation front, (2) electron antineutrinos from positron captures on neutrons, and (3) the thermal emission from pair annihilation. We estimate that a pure deflagration supernova explosion at a distance of 1 kpc would trigger about 14 events in the future 50 kt liquid scintillator detector and some 19 events in a 0.5 Mt water Cherenkov-type detector. Conclusions: While in contrast to core-collapse supernovae neutrinos carry only a very small fraction of the energy produced in the thermonuclear supernova explosion, the SN Ia neutrino signal provides information that allows us to unambiguously distinguish between different possible explosion scenarios. These studies will become feasible with the next generation of proposed neutrino observatories.

Novel method of calculating Nuclear Statistical Equilibrium is presented. Basic equations are carefully solved using arbitrary precision arithmetic. Special interpolation procedure is then used to retrieve all abundances using tabulated results for neutrons and protons, together with basic nuclear data. Proton and neutron abundance tables, basic nuclear data and partition functions for nuclides used in calculations are provided. Simple interpolation algorithm using pre-calculated p and n abundances tabulated as a functions of kT, rho and Ye is outlined. Unique properties of this method are: (1) ability to pick-up out of NSE selected nuclei only (2) computational time scaling linearly with number of re-calculated abundances (3) relatively small amount of stored data: only two large tables (4) slightly faster than solving NSE equations using traditional Newton-Raphson methods for small networks (few tens of species); superior for huge (800-3000) networks (5) do not require initial guess; works well on random input (6) can tailored to specific application (7) ability to use third-party NSE solvers to obtain fully compatible tables (8) encapsulation of the NSE code for bug-free calculations. Range of applications for this approach is possible: coverage test of traditional NSE Newton-Raphson codes, generating starting values, code-to-code verification and possible replacement of the old legacy procedures in supernova simulations.

An algorithm for integration of polynomial functions with variable weight is considered. It provides extension of the Gaussian integration, with appropriate scaling of the abscissas and weights. Method is a good alternative to usually adopted interval splitting.

Spherically symmetric, static model of the cosmological voids is constructed in the framework of the Tolman-Oppenheimer-Volkov equation with the cosmological constant. Extension of the Tooper result (dimensionless form of the TOV equation) is provided for non-zero $\Lambda$. Then, the equation is simplified in $\alpha \to 0$, $\lambda \to 0$, $\lambda/\alpha = const$ regime, suitable for largest structures in $\Lambda$-dominated universe. Voids are treated as an underdensity regions in the static Einstein universe. Both overdensity and underdensity (relative to static universe) solutions exist. They are identified with standard astrophysical spherical objects and voids, respectively. Model is tested against observed properties (the radius - the central density relation) and density profiles of voids. Analytical formulae for radial density contrast profile and radii of the voids are derived. Some consequences for cosmological n-body simulations are suggested. Hints on the dark matter/dark energy EOS filling the voids are provided.

The spectral emission of neutrinos from a plasma in nuclear statistical equilibrium (NSE) is investigated. Particular attention is paid to the possible emission of high energy (>10 MeV) neutrinos or antineutrinos. A newly developed numerical approach for describing the abundances of nuclei in NSE is presented. Neutrino emission spectra, resulting from general Fuller, Fowler, Newman (FFN) conditions, are analyzed. Regions of T-rho-Ye space favoring detectability are selected. The importance of critical Y_e values with zero net rate of neutronization (Ye dot) is discussed. Results are provided for the processing of matter under conditions typical for thermonuclear and core-collapse supernovae, pre-supernova stars, and neutron star mergers.

Spectrum of the neutrinos produced in the massive photon and longitudal plasmon decay process has been computed with four levels of approximation for the dispersion relations. Some analytical formulae in limiting cases are derived. Interesting conclusions related to previous calculations of the energy loss in stars are presented. High energy tail of the neutrino spectrum is shown to be proportional to exp(-E/kT), where E is the neutrino energy and kT is the temperature of the plasma.

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

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.

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.

Approximate analytical formula for density distribution in differentially rotating stars is derived. Any barotropic EOS and conservative rotation law can be handled with use of this method for wide range of differential rotation strength. Results are in good qualitative agreement with comparison to the other methods. Some applications are suggested and possible improvements of the formula are discussed.