On the properties of strange modes

Properties of the so-called strange modes occurring in linear stability calculations of stellar models are discussed. The behaviour of these modes is compared for two different sets of stellar models, for very massive zero-age main-sequence stars and for luminous hydrogen-deficient stars, both with high luminosity-to-mass ratios. We have found that the peculiar behaviour of the frequencies of the strange modes with the change of a control parameter is caused by the pulsation amplitude of a particular eigenmode being strongly confined to the outer part of the envelope, around the density inversion zone. The frequency of a strange mode changes because the depth of the confinement zone changes with the control parameter. Weakly non-adiabatic strange modes tend to be overstable because the amplitude confinement quenches the effect of radiative damping. On the other hand, extremely non-adiabatic strange modes become overstable because the perturbation of radiation force (gradient of radiation pressure) provides a restoring force that can be out of phase with the density perturbation. We discuss this mechanism by using a plane-parallel two-zone model.     

Comparison of the large-scale clustering in the APM and the EDSGC galaxy surveys

We present the results of a numerical study of the fluid f, p and the gravitational w modes for increasingly relativistic non-rotating polytropes. The results for f and w modes are in good agreement with previous data for uniform density stars, which supports an understanding of the nature of the gravitational wave modes based on the uniform density data. We show that the p modes can become extremely long-lived for some relativistic stars. This effect is attributed to the change in the perturbed density distribution as the star becomes more compact.     

On the r-mode spectrum of relativistic stars: the inclusion of the radiation reaction

We consider both mode calculations and time-evolutions of axial r modes for relativistic uniformly rotating non-barotropic neutron stars, using the slow-rotation formalism, in which rotational corrections are considered up to linear order in the angular velocity Ω. We study various stellar models, such as uniform density models, polytropic models with different polytropic indices n , and some models based on realistic equations of state. For weakly relativistic uniform density models and polytropes with small values of n , we can recover the growth times predicted from Newtonian theory when standard multipole formulae for the gravitational radiation are used. However, for more compact models, we find that relativistic linear perturbation theory predicts a weakening of the instability compared to the Newtonian results. When turning to polytropic equations of state, we find that for certain ranges of the polytropic index n , the r mode disappears, and instead of a growth, the time-evolutions show a rapid decay of the amplitude. This is clearly at variance with the Newtonian predictions. It is, however, fully consistent with our previous results obtained in the low-frequency approximation.     

On one classical problem in the radial orbit instability theory

Antonov’s classical problem of stability of a collisionless sphere with a purely radial motion of stars is considered as a limit of the problem in which stars move in nearly radial orbits. We provide the proper limiting equations that take into account the singularity in the density distribution at the sphere center and give their solutions. We show that there is instability for even and odd spherical harmonics, with all unstable modes being not slow. The growth rates of aperiodic even modes increase indefinitely when approaching purely radial models. The physics of the radial orbit instability is discussed.     

Nonlinear Dynamics of Ionization Fronts in HII Regions

Hydrodynamic instability of an accelerating ionization front (IF) is investigated with 2D hydrodynamic simulations, including absorption of incident photoionizing photons, recombination in the HII region, and radiative molecular cooling. When the amplitude of the perturbation is large enough, nonlinear dynamics of the IF triggered by the separation of the IF from the cloud surface is observed. This causes the second harmonic of the imposed perturbation to appear on the cloud surfaces, whereas the perturbation in density of ablated gas in the HII region remains largely single mode. This mismatch of modes between the IF and the density perturbation in the HII region prevents the strong stabilization effect seen in the linear regime. Large growth of the perturbation caused by Rayleigh-Taylor-like instability is observed late in time.     

Early star formation and galaxy formation triggered by the evolution of large-scale density perturbations

The evolution of small-scale density perturbations on the background of increasing large-scale perturbations of supercluster size will be considered. In the case that the characteristic length scales of both perturbation modes differ significantly, the interaction between both modes has to be taken into account already within lowest order of approximation. It will be shown that in this case an effective amplification for the smaller-scale perturbations occurs. For these perturbations the characteristic times of evolution decreases in dependence on the considered mass-scales more or less rapidly. Therefore, the growth of adiabatic density perturbations on mass-scales up to galaxy masses seems to be triggered by the density evolution of superclusters which the smaller-mass perturbations are embedded in. A model for the formation of observed condensed matter distribution will be proposed.     

Coupled galactic density-wave modes in a composite system of thin stellar and gaseous discs

A rotating disc galaxy is modelled as a composite system consisting of thin stellar and gaseous discs, which are described by a two-fluid modal formalism. The composite disc system is assumed to retain axisymmetry in the background equilibrium. General density-wave perturbations in the two discs are coupled through the mutual gravitational interaction. We study the basic properties of open and tight spiral density-wave modes in such a composite disc system. Within the Lindblad resonances, perturbation enhancements of surface mass density in stellar and gaseous discs are in phase; this is also true during the initial growth phase of density-wave perturbations. Outside the Lindblad resonances, there exists a possible spiral density-wave branch for which perturbation enhancements of surface mass density in stellar and gaseous discs are out of phase. We discuss implications of these results on the critical parameters for global star formation in barred and normal spiral galaxies and on magnetohydrodynamic density waves within the Lindblad resonances.     

Alfvén polar oscillations of relativistic stars

We study polar Alfvén oscillations of relativistic stars endowed with a strong global poloidal dipole magnetic field. Here, we focus only on the axisymmetric oscillations which are studied by numerically evolving the two-dimensional perturbation equations. Our study shows that the spectrum of the polar Alfvén oscillations is discrete in contrast to the spectrum of axial Alfvén oscillations which is continuous. We also show that the typical fluid modes, such as the f and p modes, are not significantly affected by the presence of the strong magnetic field.     

Third Order Effect of Rotation on Stellar Oscillations of a B Star

We aim at investigating the effect of rotation up to the third order in the angular velocity of a star on the p and g modes, based on the formalism developed by Soufi et al. Our ultimate goal is the study of oscillations of β Cephei stars which are often rapidly rotating stars. Our results show that the third-order perturbation formalism presented by Soufi et al. should be corrected for some missing terms and some misprints in the equations. As a first step in our study of β Cephei stars, we quantify by numerical calculations the effect of rotation on the oscillation frequencies of a uniformly rotating zero-age main-sequence star with 12 M<,??>. For an equatorial velocity of 100km s-1, it is found that the second-and third-order corrections for (l, m)=(2, 2), for instance, are of the order of 0.01% of the frequency for radial order n=6 and reaches up to 0.5% for n=14.     

Oscillations of rapidly rotating stratified neutron stars

We use time evolutions of the linear perturbation equations to study the oscillations of rapidly rotating neutrons stars. Our models account for the buoyancy due to composition gradients and we study, for the first time, the nature of the resultant g modes in a fast spinning star. We provide detailed comparisons of non-stratified and stratified models. This leads to an improved understanding of the relationship between the inertial modes of a non-stratified star and the g modes of a stratified system. In particular, we demonstrate that each g mode becomes rotation dominated, i.e. approaches a particular inertial mode, as the rotation rate of the star is increased. We also discuss issues relating to the gravitational wave driven instability of the various classes of oscillation modes.     

Non-linear hydrodynamical evolution of rotating relativistic stars: numerical methods and code tests

We present numerical hydrodynamical evolutions of rapidly rotating relativistic stars, using an axisymmetric, non-linear relativistic hydrodynamics code. We use four different high-resolution shock-capturing (HRSC) finite-difference schemes (based on approximate Riemann solvers) and compare their accuracy in preserving uniformly rotating stationary initial configurations in long-term evolutions. Among these four schemes, we find that the third-order piecewise parabolic method scheme is superior in maintaining the initial rotation law in long-term evolutions, especially near the surface of the star. It is further shown that HRSC schemes are suitable for the evolution of perturbed neutron stars and for the accurate identification (via Fourier transforms) of normal modes of oscillation. This is demonstrated for radial and quadrupolar pulsations in the non-rotating limit, where we find good agreement with frequencies obtained with a linear perturbation code. The code can be used for studying small-amplitude or non-linear pulsations of differentially rotating neutron stars, while our present results serve as testbed computations for three-dimensional general-relativistic evolution codes.     

Stability of power-law discs — I. The Fredholm integral equation

The power-law discs are a family of infinitesimally thin, axisymmetric stellar discs of infinite extent. The rotation curve can be rising, falling or flat. The self-consistent power-law discs are scale-free, so that all physical quantities vary as a power of radius. They possess simple equilibrium distribution functions depending on the two classical integrals, energy and angular momentum. While maintaining the scale-free equilibrium force law, the power-law discs can be transformed into cut-out discs by preventing stars close to the origin (and sometimes also at large radii) from participating in any disturbance. This paper derives the homogeneous Fredholm integral equation for the in-plane normal modes in the self-consistent and the cut-out power-law discs. This is done by linearizing the collisionless Boltzmann equation to find the response density corresponding to any imposed density and potential. The normal modes — that is, the self-consistent modes of oscillation — are found by requiring the imposed density to equal the response density. In practice, this scheme is implemented in Fourier space, by decomposing both imposed and response densities in logarithmic spirals. The Fredholm integral equation then relates the transform of the imposed density to the transform of the response density. Numerical strategies to solve the integral equation and to isolate the growth rates and the pattern speeds of the normal modes are discussed.     

Seismic signatures of strange stars with crust

We study acoustic oscillations (eigenfrequencies, velocity distributions, damping times) of normal crusts of strange stars. These oscillations are very specific because of huge density jump at the interface between the normal crust and the strange matter core. The oscillation problem is shown to be self-similar. For a low (but non-zero) multipolarity l , the fundamental mode (without radial nodes) has a frequency of ∼300 Hz and mostly horizontal oscillation velocity; other pressure modes have frequencies ≳20 kHz and almost radial oscillation velocities. The latter modes are similar to radial oscillations (having approximately the same frequencies and radial velocity profiles). The oscillation spectrum of strange stars with crust differs from the spectrum of neutron stars. If detected, acoustic oscillations would allow one to discriminate between strange stars with crust and neutron stars and constrain the mass and radius of the star.     

Surface trapping and leakage of low-frequency g modes in rotating early-type stars – II. Global analysis

A global analysis of the surface trapping of low-frequency non-radial g modes in rotating early-type stars is undertaken within the Cowling, adiabatic and traditional approximations. The dimensionless pulsation equations governing these modes are reviewed, and the boundary conditions necessary for solution of the equations are considered; in particular, an outer mechanical boundary condition, which does not enforce complete wave trapping at the stellar surface, is derived and discussed in detail. The pulsation equations are solved for a 7-M_⊙ model star over a range of rotation rates, using a numerical approach.
The results of the calculations confirm the findings of the preceding paper in the series: modes with eigenfrequencies below a cut-off cannot be fully trapped within the star, and exhibit leakage in the form of outwardly propagating waves at the surface. The damping rates resulting from leakage are calculated for such 'virtual' modes, and found to be appreciably larger than typical growth rates associated with opacity-driven pulsation. Furthermore, it is demonstrated that the surface perturbations generated by virtual modes are significantly changed from those caused by fully trapped modes; the latter result suggests differences in the line-profile variations exhibited by these two types of mode.
The findings are discussed in the context of the 53 Per, SPB and pulsating Be classes of variable star. Whilst wave leakage will probably not occur for overstable g modes in the 53 Per and slowly rotating SPB stars, the adoption of the new outer mechanical boundary condition may still affect the pulsational stability of these systems. Wave leakage for overstable modes remains a possibility in Be stars and the more rapidly rotating SPB stars.     

Accretion-induced overstability of density waves in a self-gravitating disk

A two-component differentially rotating disk of self-gravitating particles is considered in the hydrodynamical framework. This system is shown to sustain two pairs of density waves, corresponding to the familiar Jeans modes and an acoustic type of modes (similar to ion-acoustic modes in plasmas). As a result of mass and momentum transfer from the gaseous to the stellar component (an accretion process), the acoustic modes suffer a strong damping, whereas the Jeans modes which were oscillatory, now become overstable provided the thermal velocity of stars is larger than that of gaseous component. The waves with frequencies near the corotation have a rather large growth rate. This amplification can explain the maintenance of spiral structure and a selective amplification could even determine the wave-frequency (or pattern velocity).     

Chemical spots and oscillatory diffusion modes in magnetic stars

The stars of the middle main sequence often have spot‐like chemical structures at their surfaces. We consider diffusion caused by electric currents and argue that such current‐driven diffusion can form chemical inhomogeneities in a plasma. The considered mechanism can contribute to a formation of element spots in Hg‐Mn and Ap‐stars. Due to the Hall effect, diffusion in the presence of electric currents can be accompanied by the propagation of a particular type of magnetohydrodynamic modes in which only the impurity number density oscillates. Such modes exist if the magnetic pressure is much greater than the gas pressure and can be the reason for variations of the abundance peculiarities in stars. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Finite larmor-radius effects on Rayleigh-Taylor instability of a dusty magnetized plasma

This paper discusses the Rayleigh-Taylor (RT) instability of an infinitely conducting medium having an exponential density distribution which includes the effects of finite ion Larmor-radius (FLR) corrections and suspended particles in the presence of a uniform horizontal magnetic field. The relevant equations of the problem are linearized and from the linearized perturbation equations a dispersion relation is obtained, using appropriate boundary conditions. It has been found that the criterion for the stable density stratification remains uninfluenced by the simultaneous inclusion of the FLR corrections and suspended particles. The stability of the medium has been proved for the case of stable stratification when the FLR corrections are not considered in the analysis. The growth rate of unstable RT modes with increasing relaxation frequency of the suspended particles is evaluated analytically. It has been shown that the presence of suspended particles in the medium suppresses the growth rate of the unstable RT modes, thereby implying a stabilizing influence of the particles on the considered configuration.     

Detecting ionized bubbles in redshifted 21-cm maps

Linear transient phenomena induced by flow non-normality in thin self-gravitating astrophysical discs are studied using the shearing sheet approximation. The considered system includes two modes of perturbations: vortex and (spiral density) wave. It is shown that self-gravity considerably alters the vortex mode dynamics; its transient (swing) growth may be several orders of magnitude stronger than in the non-self-gravitating case and two to three times larger than the transient growth of the wave mode. Based on this finding, we comment on the role of vortex mode perturbations in a gravitoturbulent state. We also describe the linear coupling of the perturbation modes, caused by the differential character of disc rotation. The coupling is asymmetric: vortex mode perturbations are able to excite wave mode perturbations, but not vice versa. This asymmetric coupling lends additional significance to the vortex mode as a participant in spiral density waves and shock manifestations in astrophysical discs.     

A numerical study of the r-mode instability of rapidly rotating nascent neutron stars

The first results of numerical analysis of classical r-modes of rapidly rotating compressible stellar models are reported. The full set of linear perturbation equations of rotating stars in Newtonian gravity is solved numerically without the slow rotation approximation. A critical curve of gravitational wave emission induced instability, which restricts the rotational frequencies of hot young neutron stars, is obtained. Taking the standard cooling mechanisms of neutron stars into account, we also show the 'evolutionary curves' along which neutron stars are supposed to evolve as cooling and spinning down proceed. Rotational frequencies of 1.4-M_⊙ stars suffering from this instability decrease to around 100 Hz when the standard cooling mechanism of neutron stars is employed. This result confirms the results of other authors, who adopted the slow rotation approximation.     

Measuring the Cosmic Equation of State with Counts of Galaxies

We show that the gas in growing density perturbations is vulnerable to the influence of winds outflowing from nearby collapsed galaxies that have already formed stars. This suggests that the formation of nearby galaxies with masses less, similar10(9) M( middle dot in circle) is likely to be suppressed, irrespective of the details of galaxy formation. An impinging wind may shock-heat the gas of a nearby perturbation to above the virial temperature, thereby mechanically evaporating the gas, or the baryons may be stripped from the perturbation entirely if they are accelerated to above the escape velocity. We show that baryonic stripping is the most effective of these two processes, because shock-heated clouds that are too large to be stripped are able to radiatively cool within a sound crossing time, limiting evaporation. The intergalactic medium temperatures and star formation rates required for outflows to have a significant influence on the formation of low-mass galaxies are consistent with current observations, but may soon be examined directly via associated distortions in the cosmic microwave background and with near-infrared observations from the Next Generation Space Telescope, which may detect the supernovae from early-forming stars.