Negative Energy And Angular Momentum Modes Of Thin Accretion Disks

This work derives the linearized equations of motion, the Lagrangian density, the Hamiltonian density, and the canonical angular momentum density for general perturbations [∝ exp (imφ) with m = 0, ± 1, ...] of a geometrically thin self-gravitating, homentropic fluid disk including the pressure. The theory is applied to “eccentric,” m = ± 1 perturbations of a geometrically thin Keplerian disk. We find m = 1 modes at low frequencies relative to the Keplerian frequency. Further, it is shown that these modes can have negative energy and negative angular momentum. The radial propagation of these low-frequency m = 1 modes can transport angular momentum away from the inner region of a disk and thus increase the rate of mass accretion. Depending on the radial boundary conditions there can be discrete low-frequency, negative-energy, m = 1 modes.     

Dynamical instability of laminar axisymmetric flows of ideal incompressible fluid

The instability of axisymmetric flows of ideal incompressible fluid with respect to infinitesimal perturbations with the nonconservation of angular momentum is investigated by numerically integrating the differential equations of hydrodynamics. The problem has been solved for two types of rotation profiles of an unperturbed flow: with zero and nonzero pressure gradients at the flow boundaries. Both rigid and free boundary conditions have been considered. The stability of axisymmetric flows with free boundaries is of great importance in disk accretion problems. Our calculations have revealed a crucial role of the flow pattern near the boundaries in the instability of the entire main flow. When the pressure gradient at the boundaries is zero, there is such a limiting scale of perturbations in azimuthal coordinate that longer-wavelength perturbations grow, while growing shorter-wavelength perturbations do not exit. In addition, for a fixed radial flow extent, there exists a nonzero minimum amplitude of the deviation of the angular velocity from the Keplerian one at which the instability vanishes. For a nonzero pressure gradient at the boundaries, the flow is unstable with respect to perturbations of any scale and at any small deviation of the angular velocity from the Keplerian one.     

Optimal growth of small perturbations in thin gaseous disks

A thin gaseous disk with a nearly Keplerian rotation profile and free boundaries in the external gravitational field of a point gravitating object does not generate any growing perturbation eigenmodes. In spite of this, a significant transient growth of linear perturbations measured by the evolution of their total acoustic energy is possible in such a disk. This is shown within the framework of the simplest model of an inviscid polytropic thin disk with a finite radial extent in which small adiabatic perturbations that are a linear combination of neutral eigenmodes with a corotation radius beyond the outer flow boundary are considered.     

Thermomagnetic and magnetorotational instabilities of thin Keplerian discs in poloidal magnetic fields

Stability of thin hot Keplerian discs is investigated asymptotically in small disc's aspect ratio, ε. The study is carried out in the local approximation for short vertical waves in the disc‐thickness scale. Besides the radial rotation shear and the vertical magnetic field, the background configuration is characterized by a vertically near‐constant temperature profile with a small vertical gradient. The temperature‐gradient term in Ohm's law, which characterizes the thermomagnetic transport is found to be of the order of ε. The effect of the thermomagnetic transport slightly modifies the conventional magnetorotational instability (MRI), while a new thermomagnetic instability (TMI) emerges in regions of the wavenumber space where MRI is absent. Explicit solutions are obtained for a wide range of values of plasma beta, β, and thermomagnetic transport coefficient, λ. In particular, it is shown for λ ≪ 1 that the MRI dominates in weak magnetic fields, β ≫ 1, while the TMI is exhibited in strong magnetic fields, β ∼ 1, also with the growth rate of the order of inverse rotation period (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Dynamical instability of laminar axisymmetric flows of ideal compressible fluid

The instability of axisymmetric flows of inviscid compressible fluid with respect to two-dimensional infinitesimal perturbations with the nonconservation of angular momentum is investigated by numerically integrating the differential equations of hydrodynamics. The compressibility is taken into account for a homentropic flow with an adiabatic index varying over a wide range. The problem has been solved for two angular velocity profiles of an initial axisymmetric flow. In the first case, a power-law rotation profile with a finite enthalpy gradient at the flow edges has been specified. For this angular velocity profile, we show that the instability of sonic and surface gravity modes in a nearly Keplerian flow, when a radially variable vorticity exists in the main flow, can be explained by the combined action of the Landau mechanism and mode coupling. We also show that including a radially variable vorticity makes the limiting exponent in the rotation law at which the unstable surface gravity modes vanish dependent on the fluid compressibility. In the second case, a Keplerian rotation law with a quasi-sinusoidal deviation has been specified in such a way that the enthalpy gradient vanished at the flow edges. We have found than the sonic modes are then stabilized and the flow is unstable only with respect to the perturbations that also exist in an incompressible fluid.     

Turbulence‐driven angular momentum transport in modulated Kepler flows

The velocity fluctuations in a spherical shell arising from sinusoidal perturbations of a Keplerian shear flow with a free amplitude parameter ε are studied numerically by means of fully 3D nonlinear simulations. The investigations are performed at high Reynolds numbers, i.e. 3000 < Re < 5000. We find Taylor‐Proudman columns of large eddies parallel to the rotation axis for sufficiently strong perturbations. An instability sets in at critical amplitudes with ε_crit ∝ Re^—1. The whole flow turns out to be almost axisymmetric and nonturbulent exhibiting, however, a very rich radial and latitudinal structure. The Reynolds stress 〈u′_ru′_ϕ〉 is positive in the entire computational domain, from its Gaussian radial profile a positive viscosity‐alpha of about 10^—4 is derived. The kinetic energy of the turbulent state is dominated by the azimuthal component 〈u^′2_ϕ whereas the other components are smaller by two orders of magnitude. Our simulations reveal, however, that these structures disappear as soon as the perturbations are switched off. We did not find an “effective” perturbation whose amplitude is such that the disturbance is sustained for large times (cf. Dauchot & Daviaud 1995) which is due to the effective violation of the Rayleigh stability criterion. The fluctuations rapidly smooth the original profile towards to pure Kepler flow which, therefore, proves to be stable in that sense.     

Axisymmetry vs. nonaxisymmetry of a Taylor‐Couette flow with azimuthal magnetic fields

The instability of a supercritical Taylor‐Couette flow of a conducting fluid with resting outer cylinder under the influence of a uniform axial electric current is investigated for magnetic Prandtl number Pm = 1. In the linear theory the critical Reynolds number for axisymmetric perturbations is not influenced by the current‐induced axisymmetric magnetic field but all axisymmetric magnetic perturbations decay. The nonaxisymmetric perturbations with m = 1 are excited even without rotation for large enough Hartmann numbers (“Tayler instability”). For slow rotation their growth rates scale with the Alfvén frequency of the magnetic field but for fast rotation they scale with the rotation rate of the inner cylinder. In the nonlinear regime the ratio of the energy of the magnetic m = 1 modes and the toroidal background field is very low for the non‐rotating Tayler instability but it strongly grows if differential rotation is present. For super‐Alfv´enic rotation the energies in the m = 1 modes of flow and field do not depend on the molecular viscosity, they are almost in equipartition and contain only 1.5 % of the centrifugal energy of the inner cylinder. The geometry of the excited magnetic field pattern is strictly nonaxisymmetric for slow rotation but it is of the mixed‐mode type for fast rotation – contrary to the situation which has been observed at the surface of Ap stars. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Boundary layer emission in luminous LMXBs

We show that aperiodic and quasiperiodic variability of bright LMXBs – atoll and Z‐sources – on ∼ sec‐msec time scales is caused primarily by variations of the luminosity of the boundary layer. The emission of the accretion disk is less variable on these time scales and its power density spectrum follows P _disk(f ) ∝ f ^–1 law, contributing to observed flux variation at low frequencies and low energies only. The kHz QPOs have the same origin as variability at lower frequencies, i.e. independent of the nature of the “clock”, the actual luminosity modulation takes place on the neutron star surface. The boundary layer spectrum remains nearly constant in the course of the luminosity variations and is represented to certain accuracy by the Fourier frequency resolved spectrum. In the investigated range of ∼ (0.1 – 1) _Edd it depends weakly on the global mass accretion rate and in the limit ∼ _Edd is close toWien spectrum with kT ∼ 2.4 keV. Its independence on the global value of lends support to the theoretical suggestion by Inogamov & Sunyaev (1999) that the boundary layer is radiation pressure supported. Based on the knowledge of the boundary layer spectrum we attempt to relate the motion along the Z‐track to changes of physically meaningful parameters. Our results suggest that the contribution of the boundary layer to the observed emission decreases along the Z‐track from conventional ∼50% on the horizontal branch to a rather small number on the normal branch. This decrease can be caused, for example, by obscuration of the boundary layer by the geometrically thickened accretion disk at ∼ _Edd. Alternatively, this can indicate significant change of the structure of the accretion flow at ∼ _Edd and disappearance of the boundary layer as a distinct region of the significant energy release associated with the neutron star surface. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Non-axisymmetric instabilities in shocked adiabatic accretion flows

We investigate the linear stability of a shocked accretion flow on to a black hole in the adiabatic limit. Our linear analyses and numerical calculations show that, despite the post-shock deceleration, the shock is generally unstable to non-axisymmetric perturbations. The simulation results of Molteni, Tóth & Kuznetsov can be well explained by our linear eigenmodes. The mechanism of this instability is confirmed to be based on the cycle of acoustic waves between the corotation radius and the shock. We obtain an analytical formula to calculate the oscillation period from the physical parameters of the flow. We argue that the quasi-periodic oscillation should be a common phenomenon in accretion flows with angular momentum.     

Angular momentum transport by thermal emission in black hole accretion disks

We calculate the amount of angular momentum that thermal photons carry out of a viscous black hole accretion disk, due to the strong Doppler shift imparted to them by the high orbital velocity of the radiating disk material. While thermal emission can not drive accretion on its own, we show that along with disk heating it does nonetheless result in a loss of specific angular momentum, thereby contributing to an otherwise viscosity‐driven accretion flow. In particular, we show that the fraction of the angular momentum that is lost to thermal emission at a radius r in a standard, multi‐color disk is ∼0.4r_s/r, where r_s is the Schwarzschild radius of the black hole. We briefly highlight the key similarties between this effect and the closely related Poynting‐Robertson effect (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Current‐driven instabilities in weakly ionized disks

Cool weakly ionized gaseous rotating disk form the basis for many models in astrophysics objects. Instabilities against perturbations in such disks play an important role in the theory of the formation of stars and planets. Traditionally, axisymmetric magnetohydrodynamic (MHD) and recently Hall‐MHD instabilities have been thoroughly studied as providers of an efficient mechanism for radial transfer of angular momentum, and of density radial stratification. In the current work, the Hall instability against axisymmetric perturbations in incompressible rotating fluid in external poloidal and toroidal magnetic field is considered. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Relativistic magnetofluid disc with nonzero radial velocity in Schwarzschild geometry

In this paper we consider the equilibrium of a magnetofluid disc in Schwarzschild background with an external magnetic field, having the azimuthal and the radial components of the flow velocity nonzero. The electrical conductivityσ of the fluid is taken to be finite and thus the solution for the electromagnetic field is required to satisfy the Ohm’s law too with the four-current having onlyJ ^ϕ andJ ^t nonzero. The various physical parameters that have to correlate for possible equilibrium configurations are identified and their respective magnitudes estimated. It is found that for a given angular momentum distribution the inner edge of the disc can reach well within the usual6m limit only when the surface magnetic field of the central object is not too high when the matter density at the outer edge of the disc and the accretion rate are taken with reasonable limits     

Low‐temperature primordial gas in merging halos

The thermal regime of the baryons behind shock waves arising in the process of virialization of dark matter halos is governed at certain conditions by radiation of HD lines. A small fraction of the shocked gas can cool down to the temperature of the cosmic microwave background (CMB). We estimate an upper limit for this fraction: at z = 10 it increases sharply from about q_T ∼ 10^–3 for dark halos of M = 5 × 10⁷ M_⊙ to ∼ 0.1 for halos with M = 10⁸ M_⊙. Further increase of the halo mass does not lead however to a significant growth of q_T – the asymptotic value for M ≫ 10⁸ M_⊙ is 0.3. We estimate the star formation rate associated with such shock waves, and show that they can provide a small but not negligible fraction of the star formation. We argue that extremely metal‐poor low‐mass stars in the Milky Way may have been formed from primordial gas behind such shocks. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Oscillation modes of relativistic slender tori

Accretion flows with pressure gradients permit the existence of standing waves which may be responsible for observed quasi-periodic oscillations (QPO's) in X-ray binaries. We present a comprehensive treatment of the linear modes of a hydrodynamic, non-self-gravitating, polytropic slender torus, with arbitrary specific angular momentum distribution, orbiting in an arbitrary axisymmetric space–time with reflection symmetry. We discuss the physical nature of the modes, present general analytic expressions and illustrations for those which are low order, and show that they can be excited in numerical simulations of relativistic tori. The mode oscillation spectrum simplifies dramatically for near Keplerian angular momentum distributions, which appear to be generic in global simulations of the magnetorotational instability. We discuss our results in light of observations of high frequency QPO's, and point out the existence of a new pair of modes which can be in an approximate 3:2 ratio for arbitrary black hole spins and angular momentum distributions, provided the torus is radiation pressure dominated. This mode pair consists of the axisymmetric vertical epicyclic mode and the lowest order axisymmetric breathing mode.     

The Bochum survey of the southern Galactic disk: I. Survey design and first results on 50 square degrees monitored in 2011

We are monitoring a 6° wide stripe along the southern Galactic disk simultaneously in the r and i bands, using a robotic 15‐cm twin telescope of the Universitätsternwarte Bochum near Cerro Armazones in Chile. Utilising the telescope's 2.7° field of view, the survey aims at observing a mosaic of 268 fields once per month and to monitor dedicated fields once per night. The survey reaches a sensitivity from 10^m down to 18^m (AB system), with a completeness limit of r ∼ 15.5^m and i ∼ 14.5^m which – due to the instrumental pixel size of 2.″4 – refers to stars separated by >3″. This brightness range is ideally suited to examine the intermediately bright stellar population supposed to be saturated in deep variability surveys with large telescopes. To connect to deep surveys or to explore faint long term variables, coadded images of several nights reach a depth of ∼ 20^m. The astrometric accuracy is better than 1″, as determined from the overlap of neighbouring fields. We describe the survey design, the data properties and our procedures to derive the light curves and to extract variable stars. We present a list of ∼2200 variable stars identified in 50 square degrees with 50‐80 observations between May and October 2011. For bright stars the variability amplitude A reaches down to A ∼ 0.05^m, while at the faint end variations of A > 1^m are detected. About 200 stars were known tobe variable, and their amplitudes and periods – as far as determinable from our six month monitoring – agree with literature values, demonstrating the performance of the Bochum Galactic Disk Survey (© 2012 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Nonlocal density wave theory for gravitational instability of protoplanetary disks without sharp boundaries

The stability of a self‐gravitating infinitesimally thin gaseous disk rotating around a central mass is studied. Our global linear analysis concerns marginal stability, i.e. it yields the critical temperature for the onset of instability for any given ratio of the disk mass to the central mass. Both axisymmetric and low‐m nonaxisymmetric excitations are analysed. When the fractional disk mass increases, the symmetry character of the instability changes from rings (m = 0) to one‐armed trailing spirals (m = 1). The distribution of the surface density along the spiral arms is not uniform, but describes a sequence of maxima that might be identified with forming planets. The number of the mass concentrations decreases with increasing fractional disk mass. We also obtain solutions in the form of global nonaxisymmetric vortices, which are, however, never excited.     

The applications of the MHD Alfvén wave oscillation model for kHz quasi‐periodic oscillations

In this paper, we improve the previous work on the MHD Alfvén wave oscillation model for the neutron star (NS) kHz quasi‐periodic oscillations (QPOs), and compare the model with the updated twin kHz QPO data. For the 17 NS X‐ray sources with the simultaneously detected twin kHz QPO frequencies, the stellar mass M and radius R constraints are given by means of the derived parameter A in the model, which is associated with the averaged mass density of the star as 〈ρ 〉 = 3M /(4πR³) ≃ 2.4 × 10¹⁴ (A /0.7)² g/cm³, and we also compare the M ‐R constraints with the stellar equations of state. Moreover, we also discuss the theoretical maximum kHz QPO frequency and maximum twin peak separation, and some expectations on SAX J1808.4–3658 are mentioned, such as its highest kHz QPO frequency ∼ 870 Hz, which is about 1.4–1.5 times less than those of the other known kHz QPO sources. The estimated magnetic fields for both Z sources (about Eddington accretion rate ) and Atoll sources (∼ 1% ) are approximately ∼10⁹ G and ∼10⁸ G, respectively. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Oscillatory motions in an active prominence

Different types of oscillatory motions were detected in the late phases of eruption of a prominence. We found oscillations of the prominence axis and diameter with periods of 4.3 and 9.1 min, corresponding to the eigenmodes m = 4 and m = 8 with a damping factor 4.6 × 10^–3 s^–1. A period about 4.5 min was found for oscillations of the pitch angle of the helically twisted filaments. The m = 2 and m = 3 eigenmodes could be also identified and they led to the final relaxation of the prominence axis. The observations are compared with a model in which we consider forces acting in a curved, cylindrical magnetic tube anchored at both ends in the photosphere and carrying an electric current. The stability of the prominence is discussed.