Relativistic jets are intrinsic phenomena of active galactic nuclei (AGN) and quasars. They have been observed to also emanate from systems containing compact objects, such as white dwarfs, neutron stars and black hole candidates. The corresponding Lorentz factors, $\Gamma$, were found to correlate with the compactness of the central objects. In the case of quasars and AGNs, plasmas with $\Gamma$-factors larger than $8$ were detected. However, numerically consistent modelling of propagating shock-fronts with $\Gamma \geq 4$ is a difficult issue, as the non-linearities underlying the transport operators increase dramatically with $\Gamma$, thereby giving rise to a numerical stagnation of the time-advancement procedure or alternatively they may diverge completely. In this paper, we present a unified numerical solver for modelling the propagation of one-dimensional shock fronts with high Lorentz factors. The numerical scheme is based on the finite-volume formulation with adaptive mesh refinement (AMR) and domain decomposition for parallel computation. It unifies both time-explicit and time-implicit numerical schemes within the framework of the pre-conditioned defect-correction iteration solution procedure. We find that time-implicit solution procedures are remarkably superior over their time-explicit counterparts in the very high $\Gamma$-regime and therefore most suitable for consistent modelling of relativistic outflows in AGNs and micro-quasars.

The Crab and Vela are well-studied glitching pulsars and the data obtained sofar should enable us to test the reliability of models of their internal structures. Very recently it was proposed that glitching pulsars are embedded in bimetric spacetime: their incompressible superfluid cores (SuSu-cores) are embedded in flat spacetime, whereas the ambient compressible and dissipative media are enclosed in Schwarzschild spacetime. In this letter we apply this model to the Crab and Vela pulsars and show that a newly born pulsar initially of $1.25\text{M}_{\odot}$ and an embryonic SuSu-core of $0.029\text{M}_{\odot}$ could evolve into a \ARZL{Crab-like}\ARZR pulsar after 1000 years and into a \ARZL{Vela-like}\ARZR pulsar 10000 years later to finally fade away as an invisible dark energy object after roughly 10 Myr. Based thereon we infer that the Crab and the Vela pulsars should have SuSu-cores of $0.15\text{M}_{\odot}$ and $0.55\text{M}_{\odot}$, respectively. Futhermore, the under- and overshootings phenomena observed to accompany the glitch events of the Vela pulsar are rather a common phenomenon of glitching pulsars that can be well-explained within the framework of bimetric spacetime.

We present a modification of our previous model for the mechanisms underlying the glitch phenomena in pulsars and young neutron stars. Accordingly, pulsars are born with embryonic cores comprising of purely incompressible superconducting gluon-quark superfluid (henceforth SuSu-cores). As the ambient medium cools and spins down due to emission of magnetic dipole radiation, the mass and size of SuSu-cores are set to grow discretely with time, in accordance with the Onsager-Feynmann analysis of superfluidity. Presently, we propose that the spacetime embedding glitching pulsars is dynamical and of bimetric nature: inside SuSu-cores the spacetime must be flat, whereas the surrounding region, where the matter is compressible and dissipative, the spacetime is Schwarzschild. It is further proposed that the topological change of spacetime is derived by the strong nuclear force, whose operating length scales is found to increase with time to reach O(1) cm at the end of the luminous lifetimes of pulsars. The model presented here model is in line with the recent radio and GW observations of pulsars and NSs.

The observed recurrence of glitches in pulsars and neutron stars carry rich information about the evolution of their internal structures. In this article I show that the glitch-events observed in pulsars are exact quantum signatures for their metamorphosis into dark super-baryons (SBs), whose interiors are made of purely incompressible superconducting gluon-quark superfluids. Here the quantum nuclear shell model is adopted to describe the permitted energy levels of the SB, which are assumed to be identical to the discrete spinning rates $\Omega_{SB},$ that SBs are allowed to rotate with. Accordingly, a glitch-event corresponds to a prompt spin-down of the superconducting SB from one energy level to the next, thereby expelling a certain number of vortices, which in turn spins-up the ambient medium. The process is provoked mainly by the negative torque of the ambient dissipative nuclear fluid and by a universal scalar field $\phi$ at the background of a supranuclear dense matter. As dictated by the Onsager-Feynman equation, the prompt spin-down must be associated with increase of the dimensions of the embryonic SB to finally convert the entire pulsar into SB-Objects on the scale of Gyrs. Based on our calculations, a Vela-like pulsar should display billions of glitches during its lifetime, before it metamorphoses entirely into a maximally compact SB-object and disappears from our observational windows. The present model predicts the mass of SBs and $\Delta \Omega/\Omega$ in young pulsars to be relatively lower than their older counterparts.

Pulsars and neutron stars are generally more massive than the Sun, whereas black holes have unlimited mass-spectrum, though the mass-gap between 2 - 5$M_\odot$, which applies for both classes, is evident and remains puzzling. Based on the solution of the TOV equation modified to include a universal scalar field $\phi$ at the background of supranuclear densities, we claim that pulsars must be born with embryonic super-baryons (SBs), that form through the merger of individual neutrons at their centers. The cores of SBs are made of purely incompressible superconducting gluon-quark superfluids (henceforth SuSu-fluids). Such quantum fluids have a uniform supranuclear density and governed by the critical EOSs for baryonic matter $P_b = \mathcal{E}_b$ and for $\phi-$induced dark energy $P_\phi= -\mathcal{E}_\phi.$ The incompressibility here ensures that particles communicate on the shortest possible timescale, superfluidity and superconductivity enforce SBs to spin-down promptly as dictated by the Onsager-Feynman equation, whereas their lowest energy state grants SBs lifetimes that are comparable to those of protons. The extraordinary long lifetimes suggests that conglomeration of SuSu-objects would evolve over several big bang events to possibly form dark matter halos that embed the galaxies in the observable universe. Having pulsars been converted into SuSu-objects, which is predicted to last for one Gyr or even shorter, then they become extraordinary compact and turn invisible. It turns out that recent observations on the quantum, stellar and cosmological scales remarkably support the present scenario.

Astronomical observations reveal a gap in the mass spectrum of relativistic objects: neither black holes nor neutron stars with 2 - 5 solar masses have ever been observed. In this article I proceed in presenting the scenario which discloses a possible hidden connection between massive neutron stars (MANSs), dark matter and dark energy in cosmology. Accordingly, when the curved spacetime embedding MANSs compresses the nuclear matter to beyond a critical supranuclear density $n_{cr},$ mesons, generally transmitting the residual nuclear forces between neutrons, could gain energy by frequently interacting with a scalar field $\phi$ at the background. When the effective energy of mesons becomes comparable to the bag energy enclosing the quarks, the neutrons merge together and form a super-baryon (SB), whose interior is made of incompressible gluon-quark superfluid. It turns out that the process has a runaway-character: it enables the super-baryon to grow in mass and volume from inside-to-outside to finally metamorphose the entire object into a completely invisible dark gluon-quark object, practically indistinguishable from isolated stellar black holes. The inability of these objects to merge with other objects whilst agglomerating in clusters makes them excellent candidates both for black holes and for dark matter halos in cosmology.

Astronomical observations reveal a gap in the mass spectrum of relativistic objects: neither black holes nor neutron stars having masses in the range of 2 - 5$\,\MSun$ have ever been observed. Based on the solution of the TOV equation modified to include a universal scalar field $\cal{H},$ we argue that all moderate and massive neutron stars should end invisible dark energy objects (DEOs). Triggered by the $\cal{H}-$baryonic matter interaction, a phase transition from normal compressible nuclear matter into an incompressible quark-superfluid is shown to occur at roughly $3$ times the nuclear density. At the transition front, the scalar field is set to inject energy at the maximum possible rate via a non-local interaction potential $V_\phi = a_0 r^2 + b_0.$ This energy creates a global confining bag, inside which a sea of freely moving quarks is formed in line with the asymptotic freedom of quantum chromodynamics. The transition front, $r_f,$ creeps from inside-to-outside to reach the surface of the object on the scale of Gyrs or even shorter, depending on its initial compactness. Having $r_f$ reached $R_\star,$ then the total injected dark energy via $V_\phi $ turns NSs into invisible DEOs. While this may provide an explanation for the absence of stellar BHs with $M_{BH}\leq 5 \MSun$ and NSs with $M_{NS}\geq 2 \MSun$, it also suggests that DEOs might have hidden connection to dark matter and dark energy in cosmology.

Jet-plasmas emanating from the vicinity of relativistic objects and in gamma-ray bursts have been observed to propagate with Lorentz factors laying in the range between one and several hundreds. On the other hand, the numerical studies of such flows have been focussed so far mainly on the lowest possible range of Lorentz factors $\Gamma,$ specifically, on the regime $1\leq \Gamma \leq 5.$ Therefore, relativistic flows with high $\Gamma-$factors have poorly studied, as most numerical methods are found to encounter severe numerical difficulties or even become numerically unstable for $\Gamma \gg 1. In this paper we present an implicit numerical advection scheme for modeling the propagation of relativistic plasmas with shocks, discuss its consistency with respect to both the internal and total energy formulation in general relativity. Using the total energy formulation, the scheme is found to be viable for modeling moving shocks with moderate Lorentz factors, though with relatively small Courant numbers. In the limit of high Lorentz factors, the internal energy formulation in combination with a fine-tuned artificial viscosity is much more robust and efficient. We confirm our conclusions by performing test calculations and compare the results with analytical solutions of the relativistic shock tube problem. The aim of the present modification is to enhance the robustness of the general relativistic implicit radiative MHD solver: GR-I-RMHD (http://www1.iwr.uni-heidelberg.de/groups/compastro/home/gr-i-mhd-solver) and extend its range of applications into the high $\Gamma-$regime.}

Black holes are inevitable mathematical outcome of spacetime-energy coupling in general relativity. Currently these objects are of vital importance for understanding numerous phenomena in astrophysics and cosmology. However, neither theory nor observations have been capable of unequivocally prove the existence of black holes or granting us an insight of what their internal structures could look like, therefore leaving researchers to speculate about their nature. In this paper the reliability of supermassive Bose-Einstein condensates (henceforth SMBECs) as alternative to supermassive black holes is examined. Such condensates are found to suffer of a causality problem that terminates their cosmological growth toward acquiring masses typical for quasars and triggers their self-collapse into supermassive black holes (SMBHs). It is argued that a SMBEC-core most likely would be subject to an extensive deceleration of its rotational frequency as well as to vortex-dissipation induced by the magnetic fields that thread the crust, hence diminishing the superfluidity of the core. Given that rotating superfluids between two concentric spheres have been verified to be dynamically unstable to non axi-symmetric perturbations, we conclude that the remnant energy stored in the core would be sufficiently large to turn the flow turbulent and dissipative and subsequently lead to core collapse bosonova. The feasibility of a conversion mechanism of normal matter into bosonic condensates under normal astrophysical conditions is discussed as well. We finally conclude that in lack of a profound theory for quantum gravity, BHs should continue to be the favorite proposal for BH candidates.

A general relativistic model for the formation and acceleration of low mass-loaded jets from systems containing accreting black holes is presented. The model is based on previous numerical results and theoretical studies in the Newtonian regime, but modified to include the effects of space-time curvature in the vicinity of the event horizon of a spinning black hole. It is argued that the boundary layer between the Keplerian accretion disk and the event horizon is best suited for the formation and acceleration of the accretion-powered jets in active galactic nuclei and micro-quasars. The model is based on matching the solutions of three different regions: i- a weakly magnetized Keplerian accretion disk in the outer part, where the transport of angular momentum is mediated through the magentorotational instability, ii- a strongly magnetized, advection-dominated and turbulent-free boundary layer (BL) between the outer cold accretion disk and the event horizon, where the plasma rotates sub-Keplerian and iii- a transition zone (TZ) between the BL and the overlying corona, where the electrons and protons are thermally uncoupled, highly dissipative and rotate super-Keplerian. Our model predicts the known correlation between the Lorentz-factor and the spin parameter of the BH. It also shows that the effective surface of the BL, through which the baryons flow into the TZ, shrinks with increasing the spin parameter, implying therefore that low mass-loaded jets most likely originate from around Kerr black holes. When applying our model to the jet in the elliptical galaxy M87, we find a spin parameter in the range 0.99 - 0.998, a transition radius of about 30 gravitational radii and a fraction of 0.05 - 0.1 of the mass accretion rate goes into the TZ, where the plasma speeds up its outward-oriented motion to reach a Lorentz factor of 2.5 - 5.0 at the transition radius.

The distribution of angular momentum of the matter during X-ray bursts on neutron stars is studied by means of 3D axi-symmetric general relativistic hydrodynamics. The set of fully general relativistic Navier-Stokes equations is solved implicitly using the implicit solver GR-I-RMHD in combination with a third order spatial and second order temporal advection scheme. The viscous operators are formulated using a Kerr-like metric in the fixed background of a slowly rotating neutron star whose radius coincides with the corresponding last stable orbit. The importance of these operators and their possible simplifications are discussed as well. In the rotating case and depending on the viscosity parameter, $\alpha_\mm{tur}$, it is found that the viscously-initiated fronts at the center of bursts propagate at much faster speed than the fluid motion. These fast fronts act to decouple angular momentum from matter: angular momentum is transported outwards while matter sinks inwards into the deep gravitational well of the neutron star, thereby enhancing the compression of matter necessary for initiating ignition, that subsequently spreads over the whole surface of the neutron star on the viscous time scale. Based on the numerical simulations, we find that a viscosity parameter $\alpha_\mm{tur} = \mathcal{O}(0.1)$ is most suitable for fitting observations of neutron stars during X-ray bursts. It is argued that the spin up observed in the cooling tails of X-ray bursts is a transient phase, which eventually should be followed by a spin down phase. This delay can be attributed to a significant lengthening of the viscous time scale due to rapid cooling of matter in the outer layers.

An implicit algorithm for solving the equations of general relativistic hydrodynamics in conservative form in three-dimensional axi-symmetry is presented. This algorithm is a direct extension of the pseudo-Newtonian implicit radiative magnetohydrodynamical solver -IRMHD- into the general relativistic regime. We adopt the Boyer-Lindquist coordinates and formulate the hydrodynamical equations in the fixed background of a Kerr black hole. The set of equations are solved implicitly using the hierarchical solution scenario (HSS). The HSS is efficient, robust and enables the use of a variety of solution procedures that range from a purely explicit up to fully implicit schemes. The discretization of the HD-equations is based on the finite volume formulation and the defect-correction iteration strategy for recovering higher order spatial and temporal accuracies. Depending on the astrophysical problem, a variety of relaxation methods can be applied. In particular the vectorized black-white Line-Gauss-Seidel relaxation method is most suitable for modeling accretion flows with shocks, whereas the Approximate Factorization Method is for weakly compressible flows. The results of several test calculations that verify the accuracy and robustness of the algorithm are shown. Extending the algorithm to enable solving the non-ideal MHD equations in the general relativistic regime is the subject of our ongoing research.

Computational gas dynamics has become a prominent research field both in astrophysics and cosmology. In the first part of this review we intend to briefly describe several of the numerical methods used in this field, discuss their range of application and present strategies for converting conditionally-stable numerical methods into unconditionally-stable solution procedures. The underlying aim of the conversion is to enhance the robustness and unification of numerical methods and subsequently enlarge their range of applications considerably. In the second part Fabian Heitsch presents and discusses the implementation of a time-explicit MHD Boltzmann solver.

Internal flows inside gravitationally stable astrophysical objects, such as the Sun, stars and compact stars are compressed and extremely subsonic. Such low Mach number flows are usually encountered when studying for example dynamo action in stars, planets, the hydro-thermodynamics of X-ray bursts on neutron stars and dwarf novae. Treating such flows is numerically complicated and challenging task. We aim to present a robust numerical tool that enables modeling the time-evolution or quasi-stationary of stratified low Mach number flows under astrophysical conditions. It is argued that astrophysical low Mach number flows cannot be considered as an asymptotic limit of incompressible flows, but rather as highly compressed flows with extremely stiff pressure terms. Unlike the pseudo-pressure in incompressible fluids, a Possion-like treatment for the pressure would smooth unnecessarily the physically induced acoustic perturbations, thereby violating the conservation character of the compressible equations. Moreover, classical dimensional splitting techniques, such as ADI or Line-Gauss-Seidel methods are found to be unsuited for modeling compressible flows with low Mach numbers. In this paper we present a nonlinear Newton-type solver that is based on the defect-correction iteration procedure and in which the Approximate Factorization Method (AFM) is used as a preconditioner. This solver is found to be sufficiently robust and is capable of capturing stationary solutions for viscous rotating flows with Mach number as small as $\mcal{M} \approx 10^{-3},$ i.e., near the incompressibility limit.

A method for relaxing the CFL-condition, which limits the time step size in explicit methods in computational fluid dynamics, is presented. The method is based on re-formulating explicit methods in matrix form, and considering them as a special-Jacobi iteration scheme that converge efficiently if the CFL- number is less than unity. By adopting this formulation, one can design various solution methods in arbitrary dimensions that range from explicit to unconditionally stable implicit methods in which CFL-number could reach arbitrary large values. In addition, we find that adopting a specially varying time stepping scheme accelerates convergence toward steady state solutions and improves the efficiently of the solution procedure.

A self-similar solution for the 3D axi-symmetric radiative MHD equations, which revisits the formation and acceleration of accretion-powered jets in AGNs and microquasars, is presented. The model relies primarily on electromagnetic extraction of rotational energy from the disk plasma and forming a geometrically thin super-Keplerian layer between the disk and the overlying corona. The outflowing plasma in this layer is dissipative, two-temperature, virial-hot, advective and electron-proton dominated. The innermost part of the disk in this model is turbulent-free, sub-Keplerian rotating and advective-dominated. This part ceases to radiate as a standard disk, and most of the accretion energy is converted into magnetic and kinetic energies that go into powering the jet. The corresponding luminosities of these turbulent-truncated disks are discussed.

Based on theoretical arguments and quasi-stationary radiative MHD calculations, a model for an accretion-powered jet is presented. It is argued that accretion disks around BHs consist of 1) a cold, Keplerian-rotating and weakly magnetized mediumin the outer part, 2) a highly advective and turbulent-free plasma inside $r_\mm{tr}=10-20$ Schwarzschild radii, where magnetic fields are predominantly of large scale topology and in excess of thermal equipartition, and 3) an ion-dominated torus in the vicinity of the hole, where magnetic fields undergo a topological change into a monopole like-configuration. The action of magnetic fields interior to $r_\mm{tr}$ is to initiate torsional {\Alfven} waves that extract angular momentum from the disk-plasma and deposit it into the transition layer (-TL) above the disk, where the plasma is dissipative and tenuous. A significant fraction of the shear-generated toroidal magnetic field reconnects in the TL, thereby heating the plasma up to the virial-temperature and forming a super-Keplerian rotating, and hence centrifugally accelerated outflow. The strong magnetic field in the TL forces the electrons to cool rapidly yielding thereby an ion-dominated outflow. The toroidal magnetic field in the TL is in thermal equipartition with the ions, whereas the poloidal component is in equipartition with the electrons. Such a strong toroidal magnetic field is essential for increasing the jet-disk luminosity in the radio regime.

Most powerful jets are observed to emanate from accreting black hole systems. Recent quasi-stationary radiative MHD calculations reveal that jet-launching start in the innermost region of a transition layer (:TL), which is located between the disk and the corona. The plasma in the TL is dominated by virial-hot protons, advective, highly magnetic-diffusive and rotates super-Keplerian. The amplified toroidal magnetic field (TMF) in the TL reach values beyond equipartition with respect to the electron thermal energy, thereby considerably enhancing the Poynting energy flux and yielding gravitationally unbound electron-proton outflows with a positive Bernoulli number. We speculate that beaming jets in GRBs may possess similar properties, provided that the GRB-progenitors have rotational energy sufficient for forming a nearly Keplerian and massive disk.

Low cooling plasmas associated with large kinetic energies are likely to be the origin of the kpc-extended and well collimated extra-galactic jets. It is proposed that jets are launched from a layer, governed by a highly diffusive, super-Keplerian rotating and thermally dominated by virial-hot and magnetized ion-plasma. The launching layer is located between the accretion disk and the corona surrounding the nucleus. The matter in the layer is causally connected to both the disk and to the central engine. Moreover we find that coronae, in the absence of heating from below, are dynamically unstable to thermal ion-conduction, and that accretion disks become intrinsically advection-dominated. We confirm the capability of this multi-layer model to form jets by carrying out 3D axisymmetric quasi-stationary MHD calculations with high spatial resolution, and taking into account turbulent and magnetic diffusion. The new multi-layer topology accommodates several previously proposed elements for jet-initiation, in particular the ion-torus, the magneto-centrifugal and the truncated disk - advective tori models.