The Tolman--Oppenheimer--Volkoff (TOV) equations are a partially uncoupled system of nonlinear and non-autonomous ordinary differential equations which describe the structure of isotropic spherically symmetric static fluids. Nonlinearity makes finding explicit solutions of TOV systems very difficult and such solutions and very rare. In this paper we introduce the notion of pseudo-asymptotic TOV systems and we show that the space of such systems is at least fifteen-dimensional. We also show that if the system is defined in a suitable domain (meaning the extended real line), then well-behaved pseudo-asymptotic TOV systems are genuine TOV systems in that domain, ensuring the existence of new fourteen analytic solutions for extended TOV equations. The solutions are classified according to the nature of the matter (ordinary or exotic) and to the existence of cavities and singularities. It is shown that at least three of them are realistic, in the sense that they are formed only by ordinary matter and contain no cavities or singularities.

Space and astrophysical plasmas are frequently found in the regime of differential rotation, where the presence of a magnetic field can result in the magnetorotational instability, directly responsible for important phenomena such as turbulent angular momentum transport in accretion disks. In the absence of an imposed magnetic field, a nonlinear dynamo is necessary for this transport mechanism to take place. In protoplanetary disks there are regions with high density and very low temperatures, which are two necessary conditions for the Hall effect to operate, affecting the development of the dynamo and the associated turbulence. In this work we perform local magnetohydrodynamic (MHD) simulations to study transition to weak turbulence in Keplerian shear flows with Hall effect. The Hall effect is shown to lead the system to long-lived turbulent transients whose decay time follows an exponential dependence on the magnetic Reynolds number and the Hall parameter.

In this article we prove three obstruction results on the existence of equations of state in clusters of stellar systems fulfilling mass-radius relationships and some additional bound (on the mass, on the radius or a causal bound). The theorems are proved in large generality. We start with a motivating example of TOV systems and we close by applying our results in stellar systems arising from experimental data.

The classical, relativistic thin-disk theory of Novikov and Thorne (NT) predicts a maximum accretion efficiency of 40% for an optically thick, radiatively efficient accretion disk around a maximally spinning black hole (BH). However, when a strong magnetic field is introduced to numerical simulations of thin disks, large deviations in efficiency are observed, in part due to mass and energy carried by jets and winds launched by the disk or BH spin. The total efficiency of accretion can be significantly enhanced beyond that predicted by NT but it has remained unclear how the radiative component is affected. In order to study the effect of a dynamically relevant large-scale magnetic field on radiatively efficient accretion, we have performed numerical 3D general relativistic - radiative - magnetohydroynamic (GRRMHD) simulations of a disk with scale height to radius ratio of $H/R~0.1$ around a moderately spinning BH (a=0.5) using the code HARMRAD. Our simulations are fully global and allow us to measure the jet, wind, and radiative properties of a magnetically arrested disk (MAD) that is kept thin via self-consistent transport of energy by radiation using the M1 closure scheme. Our fiducial disk is MAD out to a radius of ~16R_g and the majority of the total ~13\% efficiency of the accretion flow is carried by a magnetically driven wind. We find that the radiative efficiency is slightly suppressed compared to NT, contrary to prior MAD GRMHD simulations with an ad hoc cooling function, but it is unclear how much of the radiation and thermal energy trapped in the outflows could ultimately escape.

In this paper we introduce the first results that use data extracted directly from numerical simulations as inputs to the analytic twisted disk model of Zhuravlev & Ivanov. In both numerical and analytic approaches, fully relativistic models of tilted and twisted disks having a moderate effective viscosity around a slowly rotating Kerr black hole are considered. Qualitatively, the analytic model demonstrates the same dynamics as the simulations, although with some quantitative offset. Namely, the GRMHD simulations generally give smaller variations of tilt and twist across the disk. When the black hole and the disk rotate in the same sense, the simulated tilted disk and analytic model show no sign of Bardeen-Petterson alignment, even in the innermost parts of the disk, where the characteristic time for relaxation to a quasi-stationary configuration is of the same order as the computation time. In the opposite case, when the direction of the disk's rotation is opposite to that of black hole, a partial alignment is observed, in agreement with previous theoretical estimates. Thus, both fully numerical and analytic schemes demonstrate that the Bardeen-Petterson effect may not be possible for the case of prograde rotation provided that disk's effective viscosity is sufficiently small. This may have implications in modeling of different astrophysical phenomena, such as disk spectra and jet orientation.

This paper presents our latest numerical simulations of accretion disks that are misaligned with respect to the rotation axis of a Kerr black hole. In this work we use a new, fully conservative version of the Cosmos++ general relativistic magnetohydrodynamics (GRMHD) code, coupled with an ad hoc cooling function designed to control the thickness of the disk. Together these allow us to simulate the thinnest tilted accretion disks ever using a GRMHD code. In this way, we are able to probe the regime where the dimensionless stress and scale height of the disk become comparable. We present results for both prograde and retrograde cases. The simulated prograde tilted disk shows no sign of Bardeen-Petterson alignment even in the innermost parts of the disk. The simulated retrograde tilted disk, however, does show modest alignment. The implication of these results is that the parameter space associated with Bardeen-Petterson alignment for prograde disks may be rather small, only including very thin disks. Unlike our previous work, we find no evidence for standing shocks in our simulated tilted disks. We ascribe this to the combination of small black hole spin, small tilt angle, and small disk scale height in these simulations. We also add to the growing body of literature pointing out that the turbulence driven by the magnetorotational instability in global simulations of accretion disks is not isotropic. Finally, we provide a comparison between our moderately thin, untilted reference simulation and other numerical simulations of thin disks in the literature.

We present a new technique for obtaining model fittings to VLBI images of astrophysical jets. The method minimizes a performance function proportional to the sum of the squared difference between the model and observed images. The model image is constructed by summing elliptical Gaussian sources characterized by six parameters: two-dimensional peak position, peak intensity, eccentricity, amplitude and orientation angle of the major axis. We present results for the fitting of two main benchmark jets: the first, constructed from three individual Gaussian sources, the second formed by five Gaussian sources. Both jets were analyzed by our cross-entropy technique in finite and infinite signal-to-noise regimes, the background noise chosen to mimic that found in interferometric radio maps. We show that our technique is capable of recovering the parameters of the sources with a similar accuracy to that obtained from the traditional AIPS task IMFIT when the image is relatively simple (e.g., few components). For more complex maps, our method displays superior performance in recovering the parameters of the jet components. Our methodology is also able to show quantitatively the number of individual components present in an image. An additional application of the cross-entropy technique to a real image of a BL Lac object is shown and discussed. Our results indicate that our cross-entropy technique must be used in situations involving the analysis of complex emission regions having more than three sources, even though it is substantially slower than current model fitting tasks (at least 10,000 times slower for a single processor, depending on the number of sources to be optimized). As in the case of any model fitting performed in the image plane, caution is required in analyzing images constructed from a poorly sampled (u,v) plane.

Several galaxy clusters are known to present multiple and misaligned pairs of cavities seen in X-rays, as well as twisted kiloparsec-scale jets at radio wavelengths. It suggests that the AGN precessing jets play a role in the formation of the misaligned bubbles. Also, X-ray spectra reveal that typically these systems are also able to supress cooling flows, predicted theoretically. The absence of cooling flows in galaxy clusters has been a mistery for many years since numerical simulations and analytical studies suggest that AGN jets are highly energetic, but are unable to redistribute it at all directions. We performed 3D hydrodynamical simulations of the interaction between a precessing AGN jet and the warm intracluster medium plasma, which dynamics is coupled to a NFW dark matter gravitational potential. Radiative cooling has been taken into account and the cooling flow problem was studied. We found that precession is responsible for multiple pairs of bubbles, as observed. The misaligned bubbles rise up to scales of tens of kiloparsecs, where the thermal energy released by the jets are redistributed. After $\sim 150$ Myrs, the temperature of the gas within the cavities is kept of order of $\sim 10^7$ K, while the denser plasma of the intracluster medium at the central regions reaches $T \sim 10^5$ K. The existence of multiple bubbles, at diferent directions, result in an integrated temperature along the line of sight much larger than the simulations of non-precessing jets. This result is in agreement with the observations. The simulations reveal that the cooling flows cessed $\sim 50 - 70$ Myr after the AGN jets are started.

The Perseus galaxy cluster is known to present multiple and misaligned pairs of cavities seen in X-rays, as well as twisted kiloparsec-scale jets at radio wavelengths; both morphologies suggest that the AGN jet is subject to precession. In this work we performed 3D hydrodynamical simulations of the interaction between a precessing AGN jet and the warm intracluster medium plasma, which dynamics is coupled to a NFW dark matter gravitational potential. The AGN jet inflates cavities that become buoyantly unstable and rise up out of the cluster core. We found that under certain circumstances precession can originate multiple pairs of bubbles. For the physical conditions in the Perseus cluster, multiple pairs of bubbles are obtained for a jet precession opening angle > 40 degrees acting for at least three precession periods, reproducing well both radio and X-ray maps. Based on such conditions, assuming that the Bardeen-Peterson effect is dominant, we studied the evolution of the precession opening angle of this system. We were able to constrain the ratio between the accretion disc and black hole angular momenta as 0.7 - 1.4. We were also able to constrain the present precession angle to 30 - 40 degrees, as well as the approximate age of the inflated bubbles to 100 - 150 Myrs.