Forced oscillations in a hydrodynamical accretion disk and QPOs

This is the second of a series of papers aimed to look for an explanation on the generation of high frequency quasi-periodic oscillations (QPOs) in accretion disks around neutron star, black hole, and white dwarf binaries. The model is inspired by the general idea of a resonance mechanism in the accretion disk oscillations as was already pointed out by Abramowicz and Klu’zniak (2001). In a first paper (P'etri, 2005a, paper I), we showed that a rotating misaligned magnetic field of a neutron star gives rise to some resonances close to the inner edge of the accretion disk. In this second paper, we suggest that this process does also exist for an asymmetry in the gravitational potential of the compact object. We prove that the same physics applies, at least in the linear stage of the response to the disturbance in the system. This kind of asymmetry is well suited for neutron stars or white dwarfs possessing an inhomogeneous interior allowing for a deviation from a perfectly spherically symmetric gravitational field. After a discussion on the magnitude of this deformation applied to neutron stars, we show by a linear analysis that the disk initially in a cylindrically symmetric stationary state is subject to {three kinds of resonances: a corotation resonance, a Lindblad resonance due to a driven force and a parametric resonance}. In a second part, we focus on the linear response of a thin accretion disk in the 2D limit. {Waves are launched at the aforementioned resonance positions and propagate in some permitted regions inside the disk, according to the dispersion relation obtained by a WKB analysis}. In a last part, these results are confirmed and extended via non linear hydrodynamical numerical simulations performed with a pseudo-spectral code solving Euler's equations in a 2D cylindrical coordinate frame. {We found that for a weak potential perturbation, the Lindblad resonance is the only effective mechanism producing a significant density fluctuation}. In a last step, we replaced the Newtonian potential by the so called logarithmically modified pseudo-Newtonian potential in order to take into account some general-relativistic effects like the innermost stable circular orbit (ISCO). The latter potential is better suited to describe the close vicinity of a neutron star or a black hole. However, from a qualitative point of view, the resonance conditions remain the same. The highest kHz QPOs are then interpreted as the orbital frequency of the disk at locations where the response to the resonances are maximal. It is also found that strong gravity is not required to excite the resonances.     

Jet Deflection by a Quasi-Steady-State Side Wind in the Laboratory

We present experimental data on the steady state deflection of a highly supersonic jet by a side-wind in the laboratory. The use of a long interaction region enables internal shocks to fully cross the jet, leading to the development of significantly more structure in the jet than in previous work with a similar setup (Lebedev et al., 2004). The ability to control the length of the interaction region in the laboratory allows the switch between a regime representing a clumpy jet or wind and a regime similar to a slowly varying mass loss rate. The results indicate that multiple internal oblique shocks develop in the jet and the possible formation of a second working surface as the jet attempts to tunnel through the ambient medium.     

估算有中子星作伴星的黑洞质量

在本文中，我们讨论了由中子星和黑洞构成的密近双星系统的引力效应。首先，我们清楚地解释了引力辐射存在引起的轨道圆化，并估算了圆化的时间尺度。求解了围绕史瓦西黑洞作圆轨道运动的矢量进动方程，得到了一个精确解。由进动周期和中心黑洞的质量的关系，我们得到了估算中子星伴随的黑洞质量的方法。随后我们讨论了以中子星脉冲信号随自旋进动变化来寻找黑洞的方法。     

Modelling the advance–retreat cycle of a tidewater glacier with simple sediment dynamics

We present a simple coupled glacier-sediment model to simulate the evolution of a tidewater glacier. The model is based on a consideration of the total mass budget of a glacier, whereas ice mechanics are fully parameterized. The calving rate at the glacier terminus is assumed to be proportional to the water depth. We show that the formation of a morainic shoal has a profound influence on the response of a tidewater glacier to climatic forcing. For a slow periodic forcing of the equilibrium-line altitude, the model glacier advances steadily into the estuary, builds up a submarine terminal moraine, and then retreats rapidly when a critical point is hit. For a slowly increasing equilibrium-line altitude, at a rate as low as 1 m per 5 years, we find terminus retreat rates of up to 10 km in 50 years. Our model suggests that, although the response of tidewater glaciers to external forcing is strongly nonlinear, the onset of retreat is controlled by climate change.     

Beam Mode Expansion of Corrugated Conical Horns with Phase Correcting Lens: Application to Radioastronomy Receivers

A classical radioastronomy receiver is fed with a corrugated horn and an independent lens, both placed in a cryostat to lower the noise temperature. The beam is focused and directed using a combination of elliptical and plane mirrors. This paper proposes modifying the initial feeding system by placing the lens onto the horn aperture, thereby allowing a size reduction of the horn and lens, and a simplification of their mechanical design. The profiled lens is shaped to correct the phase error on the horn aperture. A quasi-optical model of the horn-plus-lens system has been developed using a Beam Mode Expansion (BME). Results using both a hyperbolic-planar lens and a spherical-elliptical lens, as well as results obtained by using Geometrical Optics (GO) with a Kirchoff–Huygens integration to get the far-field pattern, have been compared with measurements. As a direct application, a full focusing system for the new 40-m radiotelescope at the “Centro Astronómico de Yebes” is presented for the 22, 30 and 45 GHz bands. This paper has developed a QO model for a corrugated conical horn with a phase-correcting lens. This revised version was published online in July 2006 with corrections to the Cover Date.     

The impact of a supernova explosion in a very massive binary

We consider the effect of a supernova (SN) explosion in a very massive binary that is expected to form in a portion of Population III stars with the mass higher than  100 M_⊙  . In a Population III binary system, a more massive star can result in the formation of a black hole (BH) and a surrounding accretion disc. Such BH accretion could be a significant source of the cosmic reionization in the early Universe. However, a less massive companion star evolves belatedly and eventually undergoes a SN explosion, so that the accretion disc around a BH might be blown off in a lifetime of companion star. In this paper, we explore the dynamical impact of a SN explosion on an accretion disc around a massive BH, and elucidate whether the BH accretion disc is totally demolished or not. For the purpose, we perform three-dimensional hydrodynamic simulations of a very massive binary system, where we assume a BH of  10³ M_⊙  that results from a direct collapse of a very massive star and a companion star of  100 M_⊙  that undergoes a SN explosion. We calculate the remaining mass of a BH accretion disc as a function of time. As a result, it is found that a significant portion of gas disc can survive through three-dimensional geometrical effects even after the SN explosion of a companion star. Even if the SN explosion energy is higher by two orders of magnitude than the binding energy of gas disc, about a half of disc can be left over. The results imply that the Population III BH accretion disc can be a long-lived luminous source, and therefore could be an important ionizing source in the early Universe.     

Discovery of the 'missing' mode in HR 1217 by the Whole Earth Telescope

HR 1217 is a prototypical rapidly oscillating Ap star that has presented a test to the theory of non-radial stellar pulsation. Prior observations showed a clear pattern of five modes with alternating frequency spacings of 33.3 and 34.6 μHz, with a sixth mode at a problematic spacing of 50.0 μHz (which equals  1.5×33.3 μHz)  to the high-frequency side. Asymptotic pulsation theory allowed for a frequency spacing of 34 μHz, but Hipparcos observations rule out such a spacing. Theoretical calculations of magnetoacoustic modes in Ap stars by Cunha predicted that there should be a previously undetected mode 34 μHz higher than the main group, with a smaller spacing between it and the highest one. In this Letter, we present preliminary results from a multisite photometric campaign on the rapidly oscillating Ap star HR 1217 using the 'Whole Earth Telescope'. While a complete analysis of the data will appear in a later paper, one outstanding result from this run is the discovery of a newly detected frequency in the pulsation spectrum of this star, at the frequency predicted by Cunha.     

Evolution of Polytropic Reimann Disks in Halos

The equilibrium of elliptical Riemann disks with a polytropic equation of state and their evolution under the influence of viscosity and gravitational radiation inside spheroidal halos with a relative surface mass density k. The evolutionary trajectory of a disk inside a halo with k<0.5, which is analogous to the evolution of an isolated disk, differs from that of a disk inside a denser halo.     

Evolution of astrophysical jets generated by a vortex mechanism

The nonlinear dynamics of a rotating jet is examined following its ejection from a compact gravitating object by a vortex mechanism. A scenario is described in which a dense stream expands and is subsequently transformed into a nonstationary vortex consisting of a cylindrical core and a “sheath.” At this stage of development, a converging radial flow of matter in the differentially rotating nonuniform sheath collimates the jet and speeds up the rotation of the core, as well as the flow of matter along the jet, in accordance with a power law or “explosive” instability, until the velocity discontinuity at the surface of the core approaches the sound speed. Flows of this type have low energy dissipation and can serve as unique channels for the acceleration and collimation of jet eruptions from young stars, active galactic nuclei, and quasars. Translated from Astrofizika, Vol. 52, No. 1, pp. 135–145 (February 2009).     

Quasi-viscous accretion flow – I. Equilibrium conditions and asymptotic behaviour

In a novel approach to studying viscous accretion flows, viscosity has been introduced as a perturbative effect, involving a first-order correction in the α-viscosity parameter. This method reduces the problem of solving a second-order non-linear differential equation (Navier–Stokes equation) to that of an effective first-order equation. Viscosity breaks down the invariance of the equilibrium conditions for stationary inflow and outflow solutions, and distinguishes accretion from wind. Under a dynamical systems classification, the only feasible critical points of this 'quasi-viscous' flow are saddle points and spirals. On large spatial scales of the disc, where a linearized and radially propagating time-dependent perturbation is known to cause a secular instability, the velocity evolution equation of the quasi-viscous flow has been transformed to bear a formal closeness with Schrödinger's equation with a repulsive potential. Compatible with the transport of angular momentum to the outer regions of the disc, a viscosity-limited length-scale has been defined for the full spatial extent over which the accretion process would be viable.     

Solar wind-magnetosphere interaction as simulated by a 3-D em particle code

We have studied the solar wind-magnetosphere interaction using a 3-D electromagnetic particle code. The results for an unmagnetized solar wind plasma streaming past a dipole magnetic field show the formation of a magnetopause and a magnetotail, the penetration of energetic particles into cusps and radiation belt and dawn-dusk asymmetries. The effects of interplanetary magnetic field (IMF) have been investigated in a similar way as done by MHD simulations. The simulation results with a southward IMF show the shrunk magnetosphere with great particle entry into the cusps and nightside magnetosphere. This is a signature of a magnetic reconnection at the dayside magnetopause. After a quasi-stable state is established with an unmagnetized solar wind we switched on a solar wind with an northward IMF. In this case the significant changes take place in the magnetotail. The waving motion was seen in the magnetotail and its length was shortened. This phenomena are consistent with the reconnections which occur at the high latitude magnetopause. In our simulations kinetic effects will determine the self-consistent anomalous resistivity in the magnetopause that causes reconnections.Deceased January 24, 1993; R. Bunemanet al. 1993.     

Scattered emission from a relativistic outflow and its application to gamma-ray bursts

We investigate a scenario of photon scattering by electrons within a relativistic outflow. The outflow is composed of discrete shells with different speeds. One shell emits radiation for a short duration. Some of this radiation is scattered by the shell(s) behind. We calculate in a simple two-shell model the observed scattered flux density as a function of the observed primary flux density, the normalized arrival time delay between the two emission components, the Lorentz factor ratio of the two shells and the scattering shell's optical depth. Thomson scattering in a cold shell and inverse Compton scattering in a hot shell are both considered. The results of our calculations are applied to the gamma-ray bursts and the afterglows. We find that the scattered flux from a cold slower shell is small and likely to be detected only for those bursts with very weak afterglows. A hot scattering shell could give rise to a scattered emission as bright as the X-ray shallow decay component detected in many bursts, on a condition that the isotropically equivalent total energy carried by the hot electrons is large, ∼10⁵²–10⁵⁶ erg. The scattered emission from a faster shell could appear as a late short γ-ray/MeV flash or become part of the prompt emission depending on the delay of the ejection of the shell.     

Effect of Cold Plasma Injection on Whistler Mode Instability Triggered by Perpendicular Ac Electric Field at Uranus

The effect of cold plasma injection on whistler mode instability has been studied separately for a bi-Maxwellian and a loss-cone hackground plasma with perpendicular AC electric field. The cold plasma is described by a simple Maxwellian distribution, whereas a generalized distribution function with index j that reduces to a bi-Maxwellian for j = 0 and to a loss-cone for j = 1 has been derived for a plasma in the presence of a perpendicular AC electric field, to form a hot/warm background. The dispersion relation is obtained using the method of characteristic solutions and kinetic approach. An expression for the growth rate of a system with added cold plasma injection has been calculated. Results of sample theoretical calculations for representative values of parameters suited to the magnetosphere of Uranus has been obtained. The salient features of the analysis and the results obtained in both cases have been compared and discussed. It is inferred that it is not the magnitude but the frequency of the AC field which influences the growth rate and a loss-cone background plasma has a triggering effect on the growth rate, increasing the value of the real frequency and maximum growth rate by an order of magnitude. These results may go a long way to enable one to get a better understanding of whistlers and diagnostics of plasma parameters in the Uranian magnetosphere.     

Measuring galaxy potentials using shell kinematics

We show that the kinematics of the shells seen around some elliptical galaxies provide a new, independent means for measuring the gravitational potentials of elliptical galaxies out to large radii. A numerical simulation of a set of shells formed in the merger between an elliptical and a smaller galaxy reveals that the shells have a characteristic observable kinematic structure, with the maximum line-of-sight velocity increasing linearly as one moves inward from a shell edge. A simple analytic calculation shows that this structure provides a direct measure of the gradient of the gravitational potential at the shell radius. In order to extract this information from attainable data, we have also derived a complete distribution of line-of-sight velocities for material within a shell; comparing the observed spectra of a shell to a stellar template convolved with this distribution will enable us to measure the gradient of the potential at this radius. Repeating the analysis for a whole series of nested shells in a galaxy allows the complete form of the gravitational potential as a function of radius to be mapped out. The requisite observations lie within reach of the up-coming generation of large telescopes.     

Dark matter in dwarf spheroidals – II. Observations and modelling of Draco

We present stellar radial velocity data for the Draco dwarf spheroidal (dSph) galaxy obtained using the AF2/WYFFOS instrument combination on the William Herschel Telescope. Our data set consists of 186 member stars, 159 of which have good quality velocities, extending to a magnitude   V ≈19.5  with a mean velocity precision of ≈2 km s⁻¹. As this survey is based on a high-precision photometric target list, it contains many more Draco members at large radii. For the first time, this allows a robust determination of the radial behaviour of the velocity dispersion in a dSph.
We find statistically strong evidence of a rising velocity dispersion consistent with a dark matter halo that has a gently rising rotation curve. There is a <2 σ signature of rotation about the long axis, inconsistent with tidal disruption as the source of the rising dispersion. By comparing our data set with earlier velocities, we find that Draco probably has a binary distribution and fraction comparable to those in the solar neighbourhood.
We apply a novel maximum likelihood algorithm and fit the velocity data to a two parameter spherical model with an adjustable dark matter content and velocity anisotropy. Draco is best fit by a weakly tangentially anisotropic distribution of stellar orbits in a dark matter halo with a very slowly rising rotation law  ( v _circ∝ r ^0.17)  . We are able to rule out both a mass-follows-light distribution and an extended halo with a harmonic core at the 2.5 to 3 σ significance level, depending on the details of our assumptions about Draco's stellar binary population. Our modelling lends support to the idea that the dark matter in dwarf spheroidals is distributed in the form of massive, nearly isothermal haloes.     

Implementation of an optimized Cassegrain system for radio telescopes

We present the antenna design for a radio interferometer, the Arcminute Microkelvin Imager, together with its beam pattern measurement. Our aim was to develop a low-cost system with high aperture efficiency and low ground spill across the frequency range 12–18 GHz. We use a modified Cassegrain system consisting of a commercially available paraboloidal primary mirror with a diameter of 3.7 m, and a shaped secondary mirror. The secondary mirror is oversized with respect to a ray optics design and has a surface that is bent towards the primary near its outer edge using a square term for the shaping. The antennas are simple to manufacture and therefore their cost is low. The design increased the antenna gain by approximately 10 per cent compared to a normal Cassegrain system while still maintaining low contamination from ground spill and using a simple design for the horn.     

Light propagation in a universe with spatial inhomogeneities

A light propagation in a universe which is homogeneous only in average systematically differs from a light propagation in a strictly homogeneous Universe. We demonstrate a link between this effect and the general theory of transport phenomenon in random media. The effective spatial curvature of a universe which is homogeneous only in average is introduced. This curvature is governing a light propagation in such universe. We show that the effective spatial curvature is lower than the average curvature. It implies that a universe with critical mean density looks like a space with negative curvature.