-
Truncated, Tilted Discs as a Possible Source of Quasi-Periodic Oscillations
Authors:
A. Deepika Bollimpalli,
P. Chris Fragile,
W. Janosz Dewberry,
Włodek Kluźniak
Abstract:
Many accreting black holes and neutron stars exhibit rapid variability in their X-ray light curves, termed quasi-periodic oscillations (QPOs). The most commonly observed type is the low-frequency ($\lesssim 10$ Hz), type-C QPO, while only a handful of sources exhibit high-frequency QPOs ($\gtrsim 60$ Hz). The leading model for the type-C QPO is Lense-Thirring precession of a hot, geometrically thi…
▽ More
Many accreting black holes and neutron stars exhibit rapid variability in their X-ray light curves, termed quasi-periodic oscillations (QPOs). The most commonly observed type is the low-frequency ($\lesssim 10$ Hz), type-C QPO, while only a handful of sources exhibit high-frequency QPOs ($\gtrsim 60$ Hz). The leading model for the type-C QPO is Lense-Thirring precession of a hot, geometrically thick accretion flow that is misaligned with the black hole's spin axis. However, existing versions of this model have not taken into account the effects of a surrounding, geometrically thin disc on the precessing, inner, geometrically thick flow. In Bollimpalli et. al 2023, using a set of GRMHD simulations of tilted, truncated accretion discs, we confirmed that the outer thin disc slows down the precession rate of the precessing torus, which has direct observational implications for type-C QPOs. In this paper, we provide a detailed analysis of those simulations and compare them with an aligned truncated disc simulation. We find that the misalignment of the disc excites additional variability in the inner hot flow, which is absent in the comparable aligned-disc simulations. This suggests that the misalignment may be a crucial requirement for producing QPOs. We attribute this variability to global vertical oscillations of the inner torus at epicyclic frequencies corresponding to the transition radius. This explanation is consistent with current observations of higher frequency QPOs in black hole X-ray binary systems.
△ Less
Submitted 22 December, 2023;
originally announced December 2023.
-
Solar evolution models with a central black hole
Authors:
Earl P. Bellinger,
Matt E. Caplan,
Taeho Ryu,
Deepika Bollimpalli,
Warrick H. Ball,
Florian Kühnel,
R. Farmer,
S. E. de Mink,
Jørgen Christensen-Dalsgaard
Abstract:
Hawking (1971) proposed that the Sun may harbor a primordial black hole whose accretion supplies some of the solar luminosity. Such an object would have formed within the first 1 s after the Big Bang with the mass of a moon or an asteroid. These light black holes are a candidate solution to the dark matter problem, and could grow to become stellar-mass black holes (BHs) if captured by stars. Here…
▽ More
Hawking (1971) proposed that the Sun may harbor a primordial black hole whose accretion supplies some of the solar luminosity. Such an object would have formed within the first 1 s after the Big Bang with the mass of a moon or an asteroid. These light black holes are a candidate solution to the dark matter problem, and could grow to become stellar-mass black holes (BHs) if captured by stars. Here we compute the evolution of stars having such a BH at their center. We find that such objects can be surprisingly long-lived, with the lightest black holes having no influence over stellar evolution, while more massive ones consume the star over time to produce a range of observable consequences. Models of the Sun born about a BH whose mass has since grown to approximately $10^{-6}~\rm{M}_\odot$ are compatible with current observations. In this scenario, the Sun would first dim to half its current luminosity over a span of 100 Myr as the accretion starts to generate enough energy to quench nuclear reactions. The Sun would then expand into a fully-convective star, where it would shine luminously for potentially several Gyr with an enriched surface helium abundance, first as a sub-subgiant star, and later as a red straggler, before becoming a sub-solar-mass BH. We also present results for a range of stellar masses and metallicities. The unique internal structures of stars harboring BHs may make it possible for asteroseismology to discover them, should they exist. We conclude with a list of open problems and predictions.
△ Less
Submitted 11 December, 2023;
originally announced December 2023.
-
Effect of geometrically thin discs on precessing, thick flows: Relevance to type-C QPOs
Authors:
D. A. Bollimpalli,
P. C. Fragile,
W. Kluźniak
Abstract:
Type-C quasi-periodic oscillations (QPOs) are the low-frequency QPOs most commonly observed during the hard spectral state of X-ray binary systems. The leading model for these QPOs is the Lense-Thirring precession of a hot, geometrically thick accretion flow that is misaligned with respect to the black hole spin axis. However, none of the work done to date has accounted for the effects of a surrou…
▽ More
Type-C quasi-periodic oscillations (QPOs) are the low-frequency QPOs most commonly observed during the hard spectral state of X-ray binary systems. The leading model for these QPOs is the Lense-Thirring precession of a hot, geometrically thick accretion flow that is misaligned with respect to the black hole spin axis. However, none of the work done to date has accounted for the effects of a surrounding, geometrically thin disc on this precession, as would be the case in the truncated disc picture of the hard state. To address this, we perform a set of GRMHD simulations of truncated discs misaligned with the spin axes of their central black holes. Our results confirm that the inner-hot flow still undergoes precession, though at a rate that is only 5 percent of what is predicted for an isolated, precessing torus. We find that the exchange of angular momentum between the outer, thin and the inner, thick disc causes this slow-down in the precession rate and discuss its relevance to type-C QPOs.
△ Less
Submitted 6 October, 2022;
originally announced October 2022.
-
Epicyclic frequencies of spheroidal stars with non-uniform density
Authors:
D. A. Bollimpalli
Abstract:
We consider the gravitational potential of a rotating star with non-uniform density to derive the orbital and epicyclic frequencies of the particles orbiting the star. We assume that the star is composed of concentric spheroids of constant density, with a global power-law distribution of density inside the star. At the lowest order approximation, we recover the known result for the Maclaurin spher…
▽ More
We consider the gravitational potential of a rotating star with non-uniform density to derive the orbital and epicyclic frequencies of the particles orbiting the star. We assume that the star is composed of concentric spheroids of constant density, with a global power-law distribution of density inside the star. At the lowest order approximation, we recover the known result for the Maclaurin spheroid that the maximum in the radial epicyclic frequency occurs at $r=\sqrt{2}ae$, for eccentricities $\geq 1/\sqrt{2}$. We find that the nature of these characteristic frequencies differs based on the geometry of the rotating star. For an oblate spheroid, the orbits resemble retrograde-Kerr orbits and the location of the radial epicyclic maximum approaches the stellar surface as the density variation inside the star becomes steeper. On the contrary, orbits around a prolate spheroid resemble prograde-Kerr orbits, but the marginally stable orbit does not exist for prolate-shaped stars. The orbital frequency is larger (smaller) than the Keplerian value for an oblate (prolate) star with the equality attained as $e \rightarrow 0$ or $r \rightarrow \infty$. The radial profiles of the angular velocity and the angular momentum allow for a stable accreting disc around any nature of oblate/prolate spheroid.
△ Less
Submitted 28 July, 2022;
originally announced July 2022.
-
Magneto-rotational instability in magnetically polarized discs
Authors:
Oscar M. Pimentel,
P. Chris Fragile,
F. D. Lora-Clavijo,
Bridget Ierace,
Deepika Bollimpalli
Abstract:
The magneto-rotational instability (MRI) is the most likely mechanism for transportation of angular momentum and dissipation of energy within hot, ionized accretion discs. This instability is produced through the interactions of a differentially rotating plasma with an embedded magnetic field. Like all substances in nature, the plasma in an accretion disc has the potential to become magnetically p…
▽ More
The magneto-rotational instability (MRI) is the most likely mechanism for transportation of angular momentum and dissipation of energy within hot, ionized accretion discs. This instability is produced through the interactions of a differentially rotating plasma with an embedded magnetic field. Like all substances in nature, the plasma in an accretion disc has the potential to become magnetically polarized when it interacts with the magnetic field. In this paper, we study the effect of this magnetic susceptibility, parameterized by $χ_m$, on the MRI, specifically within the context of black hole accretion. We find from a linear analysis within the Newtonian limit that the minimum wavelength of the first unstable mode and the wavelength of the fastest growing mode are shorter in paramagnetic ($χ_m>0$) than in diamagnetic ($χ_m<0$) discs, all other parameters being equal. Furthermore, the magnetization parameter (ratio of gas to magnetic pressure) in the saturated state should be smaller when the magnetic susceptibility is positive than when it is negative. We confirm this latter prediction through a set of numerical simulations of magnetically polarized black hole accretion discs. We additionally find that the vertically integrated stress and mass accretion rate are somewhat larger when the disc is paramagnetic than when it is diamagnetic. If astrophysical discs are able to become magnetically polarized to any significant degree, then our results would be relevant to properly interpreting observations.
△ Less
Submitted 24 May, 2021;
originally announced May 2021.
-
Asymmetric Jet launching
Authors:
Aleksandra Kotek,
Miljenko Čemeljić,
Włodek Kluźniak,
Deepika A. Bollimpalli
Abstract:
In resistive and viscous magnetohydrodynamical (MHD) simulations we obtain axial jets launched from the innermost magnetosphere of a star-disk system. We found that in a part of the parameter space continuous asymmetric jets, which are propagating in opposite directions, are launched. We compare the speed of propagation and rotation of obtained jets with recent observational results.
In resistive and viscous magnetohydrodynamical (MHD) simulations we obtain axial jets launched from the innermost magnetosphere of a star-disk system. We found that in a part of the parameter space continuous asymmetric jets, which are propagating in opposite directions, are launched. We compare the speed of propagation and rotation of obtained jets with recent observational results.
△ Less
Submitted 7 June, 2020;
originally announced June 2020.
-
Looking for the underlying cause of black hole X-ray variability in GRMHD simulations
Authors:
D. A. Bollimpalli,
R. Mahmoud,
C. Done,
P. C. Fragile,
W. Kluźniak,
R. Narayan,
C. J. White
Abstract:
Long-term observations have shown that black hole X-ray binaries exhibit strong, aperiodic variability on time-scales of a few milliseconds to seconds. The observed light curves display various characteristic features like a log-normal distribution of flux and a linear rms-flux relation, which indicate that the underlying variability process is stochastic in nature. It is also thought to be intrin…
▽ More
Long-term observations have shown that black hole X-ray binaries exhibit strong, aperiodic variability on time-scales of a few milliseconds to seconds. The observed light curves display various characteristic features like a log-normal distribution of flux and a linear rms-flux relation, which indicate that the underlying variability process is stochastic in nature. It is also thought to be intrinsic to accretion. This variability has been modelled as inward propagating fluctuations of mass accretion rate, although the physical process driving the fluctuations remains puzzling. In this work, we analyse five exceptionally long duration general relativistic magnetohydrodynamic (GRMHD) simulations of optically thin, geometrically thick, black hole accretion flows to look for hints of propagating fluctuations in the simulation data. We find that the accretion profiles from these simulations do show evidence for inward propagating fluctuations below the viscous frequency by featuring strong radial coherence and positive time lags when comparing smaller to larger radii, although these time lags are generally shorter than the viscous time-scale and frequency independent. Our simulations also support the notion that the fluctuations in $\dot{M}$ build up in a multiplicative manner, as the simulations exhibit linear rms-mass flux relations, as well as log-normal distributions of their mass fluxes. When combining the mass fluxes from the simulations with an assumed emissivity profile, we additionally find broad agreement with observed power spectra and time lags, including a recovery of the frequency dependency of the time lags.
△ Less
Submitted 27 July, 2020; v1 submitted 1 June, 2020;
originally announced June 2020.
-
Atmospheric oscillations provide simultaneous measurement of neutron star mass and radius
Authors:
D. A. Bollimpalli,
M. Wielgus,
D. Abarca,
W. Kluźniak
Abstract:
Neutron stars with near-Eddington observable luminosities were shown to harbor levitating atmospheres, suspended above their surface. We report a new method to simultaneously measure the mass and radius of a neutron star based on oscillations of such atmospheres. In this paper, we present an analytic derivation of a family of relativistic, oscillatory, spherically symmetric eigenmodes of the optic…
▽ More
Neutron stars with near-Eddington observable luminosities were shown to harbor levitating atmospheres, suspended above their surface. We report a new method to simultaneously measure the mass and radius of a neutron star based on oscillations of such atmospheres. In this paper, we present an analytic derivation of a family of relativistic, oscillatory, spherically symmetric eigenmodes of the optically and geometrically thin levitating atmospheres, including the damping effects induced by the radiation drag. We discover characteristic maxima in the frequencies of the damped oscillations and show that using the frequency maxima, one can estimate mass and radius of the neutron star, given the observed frequency and the corresponding luminosity of the star during the X-ray burst. Thus, our model provides a new way to probe the stellar parameters. We also show that the ratio of any two undamped eigenfrequencies depends only on the adiabatic index of the atmosphere, while for the damped eigenfrequencies, this ratio varies with the luminosity. The damping coefficient is independent of the mode number of the oscillations. Signatures of these atmospheres' dynamics will be reflected in the source's X-ray light curves.
△ Less
Submitted 2 July, 2019; v1 submitted 4 December, 2018;
originally announced December 2018.
-
Disc instabilities and nova eruptions in symbiotic systems: RS Ophiuchi and Z Andromedae
Authors:
D. A. Bollimpalli,
J. -M. Hameury,
J. -P. Lasota
Abstract:
Using the disc instability model for dwarf novae and soft X-ray transients, we investigate the stability of accretion discs in long-period binary systems. We simulate outbursts due to this thermal-viscous instability for two symbiotic systems, RS~Ophiuchi and Z~Andromedae. The outburst properties deduced from our simulations suggest that, although the recurrent nova events observed in RS~Oph are d…
▽ More
Using the disc instability model for dwarf novae and soft X-ray transients, we investigate the stability of accretion discs in long-period binary systems. We simulate outbursts due to this thermal-viscous instability for two symbiotic systems, RS~Ophiuchi and Z~Andromedae. The outburst properties deduced from our simulations suggest that, although the recurrent nova events observed in RS~Oph are due to a thermonuclear runaway at the white dwarf surface, these runaways are triggered by accretion disc instabilities. In quiescence, the disc builds up its mass and it is only during the disc-instability outburst that mass is accreted on to the white dwarf at rates comparable to or larger than the mass-transfer rate. For a mass-transfer rate in the range $10^{-8}$ to $10^{-7}~{\rm M}_{\odot}$ yr$^{-1}$, the accretion rate and the mass accreted are sufficient to lead to a thermonuclear runaway during one of a series of a few dwarf nova outbursts, barely visible in the optical, but easily detectable in X-rays. In the case of Z~And, persistent irradiation of the disc by the very hot white-dwarf surface strongly modifies the dwarf-nova outburst properties, making them significant only for very high mass-transfer rates, of the order of $10^{-6}~{\rm M}_{\odot}$ yr$^{-1}$, much higher than the expected secular mean in this system. It is thus likely that the so-called `combination nova' outburst observed in years 2000 to 2002 was triggered not by a dwarf-nova instability but by a mass-transfer enhancement from the giant companion, leading to an increase in nuclear burning at the accreting white-dwarf surface.
△ Less
Submitted 21 October, 2018; v1 submitted 21 April, 2018;
originally announced April 2018.
-
Radial modes of levitating atmospheres around Eddington-luminosity neutron stars
Authors:
Deepika Ananda Bollimpalli,
Włodek Kluźniak
Abstract:
We consider an optically thin radiation-supported levitating atmosphere suspended well above the stellar surface, as discussed recently in the Schwarzschild metric for a star of luminosity close to the Eddington value. Assuming the atmosphere to be geometrically thin and polytropic, we investigate the eigenmodes and calculate the frequencies of the oscillations of the atmosphere in Newtonian forma…
▽ More
We consider an optically thin radiation-supported levitating atmosphere suspended well above the stellar surface, as discussed recently in the Schwarzschild metric for a star of luminosity close to the Eddington value. Assuming the atmosphere to be geometrically thin and polytropic, we investigate the eigenmodes and calculate the frequencies of the oscillations of the atmosphere in Newtonian formalism. The ratio of the two lowest eigenfrequencies is $\sqrt{γ+1}$, i.e., it only depends on the adiabatic index.
△ Less
Submitted 3 October, 2017; v1 submitted 12 March, 2017;
originally announced March 2017.
-
Perturbation of mass accretion rate, associated acoustic geometry and stability analysis
Authors:
Deepika A Bollimpalli,
Sourav Bhattacharya,
Tapas K Das
Abstract:
We investigate the stability of stationary integral solutions of an ideal irrotational fluid in a general static and spherically symmetric background, by studying the profile of the perturbation of the mass accretion rate. We consider low angular momentum axisymmetric accretion flows for three different accretion disk models and consider time dependent and radial linear perturbation of the mass ac…
▽ More
We investigate the stability of stationary integral solutions of an ideal irrotational fluid in a general static and spherically symmetric background, by studying the profile of the perturbation of the mass accretion rate. We consider low angular momentum axisymmetric accretion flows for three different accretion disk models and consider time dependent and radial linear perturbation of the mass accretion rate. First we show that the propagation of such perturbation can be determined by an effective $2\times2$ matrix, which has qualitatively similar acoustic causal properties as one obtains via the perturbation of the velocity potential. Next, using this matrix we analytically address the stability issues, for both standing and travelling wave configurations generated by the perturbation. Finally, based on this general formalism we briefly discuss the explicit example of the Schwarzschild spacetime and compare our results of stability with the existing literature, which instead address this problem via the perturbation of the velocity potential.
△ Less
Submitted 6 September, 2016; v1 submitted 8 July, 2014;
originally announced July 2014.