In this paper, we have derived the field equations in an extended theory of gravity in an anisotropic space time background and in the presence of magnetic field. The physical and geometrical parameters of the models are determined with respect to the Hubble parameter using some algebraic approaches. A time varying scale factor has been introduced to analyze the behavior of the model. From some diagnostic approach, we found that the model behaves as $\Lambda_{CDM}$ model at late time of cosmic evolution.

We use global magnetohydrodynamic simulations to study the influence of net vertical magnetic fields on the structure of geometrically thin ($H/r \approx 0.05$) accretion disks in the Newtonian limit. We consider initial mid-plane gas to magnetic pressure ratios $\beta_0 = 1000,\, 300$ and $100$, spanning the transition between weakly and strongly magnetized accretion regimes. We find that magnetic pressure is important for the disks' vertical structure in all three cases, with accretion occurring at $z/R\approx 0.2$ in the two most strongly magnetized models. The disk midplane shows outflow rather than accretion. Accretion through the surface layers is driven mainly by stress due to coherent large scale magnetic field rather than by turbulent stress. Equivalent viscosity parameters measured from our simulations show similar dependencies on initial $\beta_0$ to those seen in shearing box simulations, though the disk midplane is not magnetic pressure dominated even for the strongest magnetic field case. Winds are present but are not the dominant driver of disk evolution. Over the (limited) duration of our simulations, we find evidence that the net flux attains a quasi-steady state at levels that can stably maintain a strongly magnetized disk. We suggest that geometrically thin accretion disks in observed systems may commonly exist in a magnetically "elevated" state, characterized by non-zero but modest vertical magnetic fluxes, with potentially important implications for disk phenomenology in X-ray binaries (XRBs) and active galactic nuclei (AGN).

We perform detailed variability analysis of two-dimensional viscous, radiation hydrodynamic numerical simulations of Shakura-Sunyaev thin disks around a stellar mass black hole. Disk models are initialized on both the gas-, as well as radiation-, pressure-dominated branches of the thermal equilibrium curve, with mass accretion rates spanning the range from $\dot{M} = 0.01 L_\mathrm{Edd}/c^2$ to $10 L_\mathrm{Edd}/c^2$. An analysis of temporal variations of the numerically simulated disk reveals multiple robust, coherent oscillations. Considering the local mass flux variability, we find an oscillation occurring at the maximum radial epicyclic frequency, $3.5\times 10^{-3}\,t_\mathrm{g}^{-1}$, a possible signature of a trapped fundamental ${\it g}$-mode. Although present in each of our simulated models, the trapped ${\it g}$-mode feature is most prominent in the gas-pressure-dominated case. The total pressure fluctuations in the disk suggest strong evidence for standing-wave ${\it p}$-modes, some trapped in the inner disk close to the ISCO, others present in the middle/outer parts of the disk. Knowing that the trapped ${\it g}$-mode frequency and maximum radial epicyclic frequency differ by only $0.01\%$ in the case of a non-rotating black hole, we simulated an additional initially gas-pressure-dominated disk with a dimensionless black hole spin parameter $a_* = 0.5$. The oscillation frequency in the spinning black hole case confirms that this oscillation is indeed a trapped ${\it g}$-mode. All the numerical models we report here also show a set of high frequency oscillations at the vertical epicyclic and breathing mode frequencies. The vertical oscillations show a 3:2 frequency ratio with oscillations occurring approximately at the radial epicyclic frequency, which could be of astrophysical importance in observed twin peak, high-frequency quasi-periodic oscillations.

Recent developments in compact object astrophysics, especially the discovery of merging neutron stars by LIGO, the imaging of the black hole in M87 by the Event Horizon Telescope (EHT) and high precision astrometry of the Galactic Center at close to the event horizon scale by the GRAVITY experiment motivate the development of numerical source models that solve the equations of general relativistic magnetohydrodynamics (GRMHD). Here we compare GRMHD solutions for the evolution of a magnetized accretion flow where turbulence is promoted by the magnetorotational instability from a set of nine GRMHD codes: Athena++, BHAC, Cosmos++, ECHO, H-AMR, iharm3D, HARM-Noble, IllinoisGRMHD and KORAL. Agreement between the codes improves as resolution increases, as measured by a consistently applied, specially developed set of code performance metrics. We conclude that the community of GRMHD codes is mature, capable, and consistent on these test problems.

We study the oscillations of an axisymmetric, viscous, radiative, general relativistic hydrodynamical simulation of a geometrically thin disk around a non-rotating, $6.62\,M_\odot$ black hole. The numerical setup is initialized with a Novikov-Thorne, gas-pressure-dominated accretion disk, with an initial mass accretion rate of $\dot{m} = 0.01\,L_\mathrm{Edd}/c^2$ (where $L_\mathrm{Edd}$ is the Eddington luminosity and $c$ is the speed of light). Viscosity is treated with the $\alpha$-prescription. The simulation was evolved for about $1000$ Keplerian orbital periods at three Schwarzschild radii (ISCO radius). Power density spectra of the radial and vertical fluid velocity components, the total (gas $+$ radiation) midplane pressure, and the vertical component of radiative flux from the photosphere, all reveal strong power at the local breathing oscillation frequency. The first, second and third harmonics of the breathing oscillation are also clearly seen in the data. We quantify the properties of these oscillations by extracting eigenfunctions of the radial and vertical velocity components and total pressure. This confirms that these oscillations are associated with breathing motion.

We present results from general relativistic calculations of the tidal disruption of white dwarf stars from near encounters with intermediate mass black holes. We follow the evolution of 0.2 and $0.6 M_\odot$ stars on parabolic trajectories that approach $10^3$ - $10^4 M_\odot$ black holes as close as a few Schwarzschild radii at periapsis, paying particular attention to the effect tidal disruption has on thermonuclear reactions and the synthesis of intermediate to heavy ion elements. These encounters create diverse thermonuclear environments characteristic of Type I supernovae and capable of producing both intermediate and heavy mass elements in arbitrary ratios, depending on the strength (or proximity) of the interaction. Nuclear ignition is triggered in all of our calculations, even at weak tidal strengths $\beta \sim 2.6$ and large periapsis radius $R_P \sim 28$ Schwarzschild radii. A strong inverse correlation exists between the mass ratio of calcium to iron group elements and tidal strength, with $\beta \lesssim 5$ producing predominately calcium-rich debris. At these moderate to weak interactions, nucleosynthesis is not especially efficient, limiting the total mass and outflows of calcium group elements to $< 15$\% of available nuclear fuel. Iron group elements however continue to be produced in greater quantity and ratio with increasing tidal strength, peaking at $\sim 60$\% mass conversion efficiency in our closest encounter cases. These events generate short bursts of gravitational waves with characteristic frequencies 0.1-0.7 Hz and strain amplitudes $0.5\times10^{-22}$ - $3.5\times10^{-22}$ at 10 Mpc source distance.

We present results from two-dimensional, general relativistic, viscous, radiation hydrodynamic numerical simulations of Shakura-Sunyaev thin disks accreting onto stellar mass Schwarzschild black holes. We consider cases on both the gas- and radiation-pressure-dominated branches of the thermal equilibrium curve, with mass accretion rates spanning the range from $\dot{M} = 0.01 L_\mathrm{Edd}/c^2$ to $10 L_\mathrm{Edd}/c^2$. The simulations directly test the stability of this standard disk model on the different branches. We find clear evidence of thermal instability for all radiation-pressure-dominated disks, resulting universally in the vertical collapse of the disks, which in some cases then settle onto the stable, gas-pressure-dominated branch. Although these results are consistent with decades-old theoretical predictions, they appear to be in conflict with available observational data from black hole X-ray binaries. We also find evidence for a radiation-pressure-driven instability that breaks the unstable disks up into alternating rings of high and low surface density on a timescale comparable to the thermal collapse. Since radiation is included self-consistently in the simulations, we are able to calculate lightcurves and power density spectra (PDS). For the most part, we measure radiative efficiencies (ratio of luminosity to mass accretion rate) close to 6%, as expected for a non-rotating black hole. The PDS appear as broken power laws, with a break typically around 100 Hz. There is no evidence of significant excess power at any frequencies, i.e. no quasi-periodic oscillations are observed.

We describe some recent results on the evolution of accretion disks around black holes and neutron stars obtained in magnetohydrodynamic (MHD) simulations in general relativity (GR) with the inclusion of radiation (GRRMHD).

We present results of a set of three-dimensional, general relativistic radiation magnetohydrodynamics simulations of thin accretion discs around a non-rotating black hole to test their thermal stability. We consider two cases, one that is initially radiation pressure dominated and expected to be thermally unstable and another that is initially gas-pressure dominated and expected to remain stable. Indeed, we find that cooling dominates over heating in the radiation pressure dominated model, causing the disc to collapse vertically on roughly the local cooling timescale. We also find that heating and cooling within the disc have a different dependence on the mid-plane pressure, a prerequisite of thermal instability. Comparison of our data with the relevant thin-disc thermal equilibrium curve suggests that our disc may be headed for the thermally stable, gas-pressure-dominated branch. However, because the disc collapses to the point that we are no longer able to resolve it, we had to terminate the simulation. On the other hand, the gas pressure dominated model, which was run for twice as long as the radiation pressure dominated one, remains stable, with heating and cooling roughly in balance. Finally, the radiation pressure dominated simulation shows some evidence of viscous instability. The strongest evidence is in plots of surface density, which show the disc breaking up into rings.

We study the dynamics of radiation pressure supported tori around Schwarzschild black holes, focusing on their oscillatory response to an external perturbation. Using KORAL, a general relativistic radiation hydrodynamics code capable of modeling all radiative regimes from the optically thick to the optically thin, we monitor a sample of models at different initial temperatures and opacities, evolving them in two spatial dimensions for $\sim 165$ orbital periods. The dynamics of models with high opacity is very similar to that of purely hydrodynamics models, and it is characterized by regular oscillations which are visible also in the light curves. As the opacity is decreased, the tori quickly and violently migrate towards the gas-pressure dominated regime, collapsing towards the equatorial plane. When the spectra of the $L_2$ norm of the mass density are considered, high frequency inertial-acoustic modes of oscillations are detected (with the fundamental mode at a frequency $68 M_{\rm BH}^{-1}\,\rm Hz$), in close analogy to the phenomenology of purely hydrodynamic models. An additional mode of oscillation, at a frequency $129 M_{\rm BH}^{-1}\,\rm Hz$, is also found, which can be unambiguously attributed to the radiation. The spectra extracted from the light curves are typically more noisy, indicating that in a real observation such modes would not be easily detected.

We simulate an oscillating purely hydrodynamical torus with constant specific angular mo- mentum around a Schwarzschild black hole. The goal is to search for quasi-periodic oscil- lations (QPOs) in the light curve of the torus. The initial torus setup is subjected to radial, vertical and diagonal (combination of radial and vertical) velocity perturbations. The hydro- dynamical simulations are performed using the general relativistic magnetohydrodynamics code Cosmos++ and ray-traced using the GYOTO code. We found that a horizontal velocity perturbation triggers the radial and plus modes, while a vertical velocity perturbation trig- gers the vertical and X modes. The diagonal perturbation gives a combination of the modes triggered in the radial and vertical perturbations.

We present the spectrum of eigenfrequencies of axisymmetric acoustic-inertial oscillations of thin accretion disks for a Schwarzschild black hole modeled with a pseudo-potential. There are nine discrete frequencies, corresponding to trapped modes. Eigenmodes with nine or more radial nodes in the inner disk belong to the continuum, whose frequency range starts somewhat below the maximum value of the radial epicyclic frequency. The results are derived under the assumption that the oscillatory motion is parallel to the midplane of the disk.

We investigated a semi-analytic and numerical model to study the geometrically thin and optically thick accretion disk around Maclaurin spheroid (MS). The main interest is in the inner region of the so called {\alpha}-disk, {\alpha} being the viscosity parameter. Analytical calculations are done assuming radiation pressure and gas pressure dominated for close to Eddington mass accretion rate and $\dot{M}\lesssim 0.1\dot{M_{Edd}}$ respectively. We found that the change in eccentricity of MS gives a change at high frequency region in the emitted spectra. We found that disk parameters are dependent on eccentricity of MS. Our semi-analytic results show that qualitatively an increase in eccentricity of MS has same behavior as decrease in mass accretion rate. Numerical work has been carried out to see the viscous time evolution of the accretion disk around MS. In numerical model we showed that if the eccentricity of the object is high the matter will diffuse slowly during its viscous evolution. This gives a clue that how spin-up or spin-down can change the time evolution of the accretion disk using a simple Newtonian approach. The change in spectra can be used to determine the eccentricity of MS and thus period of the MS.

We present eigenfrequencies and eigenfunctions of trapped acoustic-inertial oscillations of thin accretion disks for a Schwarzschild black hole and a rapidly rotating Newtonian star (a Maclaurin spheroid). The results are derived in the formalism of Nowak and Wagoner (1991) with the assumption that the oscillatory motion is parallel to the midplane of the disk. The first four radial modes for each of five azimuthal modes $m = 0$ through $m = 4$ are presented. The frequencies and wavefunctions of the lowest modes may be accurately approximated by the Airy function.

Optically thin coronae around neutron stars suffering an X-ray burst can be ejected as a result of rapid increase in stellar luminosity. In general relativity (GR), radiation pressure from the central luminous star counteracts gravitational attraction more strongly than in Newtonian physics. However, motion near the neutron star is very effectively impeded by the radiation field. We discuss coronal ejection in a general relativistic calculation of the motion of a test particle in a spherically symmetric radiation field. At every radial distance from the star larger than that of the ISCO, and any initial luminosity of the star, there exists a luminosity change which leads to coronal ejection. The luminosity required to eject from the system the inner parts of the optically thin neutron-star corona is very high in the presence of radiation drag and always close to the Eddington luminosity. Outer parts of the corona, at a distance of ~20 $R_G$ or more, will be ejected by a sub-Eddington outburst. Mildly fluctuating luminosity will lead to dissipation in the plasma and may explain the observed X-ray temperatures of coronae in low mass X-ray binaries (LMXBs). At large radial distances from the star ($3\cdot 10^3 R_G$ or more) the results do not depend on whether or not Poynting-Robertson drag is included in the calculation.