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.     

The coupling coefficients of pulsation for radiative stellar models

The second order theory of coupling is discussed regarding the radial pulsation of stellar models which are constructed ignoring convection. The formula including the nonadiabatic effect is presented. Numerical values given for model classical cepheids are considerably greater than the adiabatic values.     

PULSE: a finite element code for solving adiabatic nonradial pulsation equations

We briefly present the nonradial adiabatic pulsation code PULSE first developped for white dwarf asteroseismology and now used to compute adiabatic oscillation properties for various types of stellar objects. Numerical tests show that the code is able to provide the accuracy (for a given stellar model) required to deal with the precision in frequency expected from the COROT long runs. While the ultimate objective is to compare the output of various pulsation codes (see these proceedings), we already emphasize problems that need to be addressed concerning, in particular, the mesh resolution of the input stellar models and its impact on the accuracy at which frequencies can be computed.     

A New Code for Nonradial Stellar Pulsations and its Application to Low-Mass, Helium White Dwarfs

We present a finite difference code intended for computing linear, adiabatic, non radial pulsations of spherical stars. This code is based on a slight modification of the general Newton-Raphson technique in order to handle the relaxation of the eigenvalue(square of the eigenfrequency) of the modes and their corresponding eigenfunctions. This code has been tested computing the pulsation spectra of polytropic spheres finding a good agreement with previous work. Then, we have coupled this code to our evolutionary code and applied it to the computation of the pulsation spectrum of a low mass, pure-helium white dwarf of 0.3 M _⊙ for a wide range of effective temperatures. In making this calculation we have taken an evolutionary time step short enough such that eigenmodes corresponding to a given model are used as initial approximation to those of the next one. Specifically, we have computed periods, period spacing, eigenfunctions, weight functions, kinetic energies and variational periods for a wide range of modes. To our notice this is the first effort in studying the pulsation properties of helium white dwarfs. The solution we have found working with these realistic white dwarf models are in good accord with the predictions of the asymptotic theory of Tassoul (1980) for high order modes. This indicates that the code presented here is able to work adequately also with realistic stellar models. This revised version was published online in July 2006 with corrections to the Cover Date.     

Bifurcation in hydrodynamic models of stellar pulsation

Phenomena of bifurcation in hydrodynamic stellar models of radial pulsation are reviewed. By changing control parameters of models, we can see qualitatively different pulsation behaviors in hydrodynamic models with transitions due to various types of bifurcation.In weakly dissipative models (classical Cepheids), the bifurcation is induced by modal resonances. Two types of the modal resonances found in models are discussed: The higherharmonic resonances of the second overtone mode in the fundamental mode pulsator and of the fourth overtone mode in the first overtone pulsator are relevant to observations. The subharmonic resonance between the fundamental and first overtone modes is confirmed in classical Cepheid models.In strongly dissipative models (less-massive supergiant stars), the bifurcation of nonlinear pulsation is induced by the hydrodynamics of ionization zones as well as modal resonances. The sequence of the bifurcation sometimes leads to chaotic behaviors in nonlinear pulsation. The transition routes from regular to the chaotic pulsations found in models are discussed with respect to the theory of chaos in simple dynamical systems: The cascade of period-doubling bifurcation is confirmed to cause chaotic pulsation in W Virginis models. For models of higher luminosity, the tangent bifurcation is found to lead intermittent chaos.Finally, hydrodynamic models for chaotic pulsation with small amplitudes observed in the post-AGB stars are briefly discussed.     

An essay on stellar oscillations and evolution

It is cautioned that solar models adjusted in such a way as to achieve a match between theoretical solar oscillation characteristics and observed ones may produce neutrino fluxes inconsistent with the observations and that this is likely to be explicable as a deficiency in modeling that portion of the envelope which is most strongly affected by uncertainties in the treatment of convection. Then follows a summary of how the results of pulsation theory and of stellar evolution theory have been used together to learn about the structure and evolution of RR Lyrae stars, classical Cepheids, and high luminosity AGB stars.     

On properties of the fourier harmonics in the limit-cycle models of classical cepheids

The Fourier coefficients of the hydrodynamic variables are calculated for the limit-cycle models of classical Cepheids having periods from 7.2 to 10.9 days. In adiabatically pulsating layers of the stellar envelope, each Fourier harmonic of orderk 8 is shown to be identified with a corresponding standing wave, so that the pulsation motions of the adiabatic layers may be represented as a superposition of standing waves. Each Fourier harmonic of orderk may also be identified with the eigenfunction of orderl of the linear adiabatic wave equation when the resonance condition _l /₀ =k is fulfilled. The spectra of the oscillatory moment of inertia and the spectra of kinetic energy obey the power law for the Fourier harmonics of orderk 15, the spectrum slope being steeper for shorter pulsation periods. In the helium and hydrogen ionizing regions all of the Fourier harmonics drive the pulsation instability, whereas in the radiative damping region the mechanical work done by each Fourier harmonic is negative. In classical Cepheids having periods shorter than 10 days the period dependence of the secondary bump is due to phase changes of the second order Fourier harmonic in the outer nonadiabatic layers of the stellar envelope. At a pulsation period of II 9.7 days the second order Fourier harmonic is identified with the second overtone. At periods II > 10 days the second order Fourier harmonic tends to be attracted by the fundamental mode in such a way that their phases coincide in the outer layers of the stellar envelope.     

Coupled-oscillators

Modal coupling oscillation models for the stellar radial pulsation and coupled-oscillators are reviewed. Coupled-oscillators with the second-order and third-order terms seemed to behave non-systematically. Using the equation by Schwarzschild and Savedoff (1949) with the dissipation term of van del Pol's type which is third-order, we demonstrate the effect of each term. The effects can be understood by the terms of the nonlinear dynamics, which is recently developing, that is. phase-locking, quasi-periodicity, period doubling, and chaos. As the problem of stellar pulsation, especially of double-mode cepheids on the period-ratio, we examine the dependence on the stellar structure from which the coupling constants in the second-order terms are derived. Eigen functions for adiabatic pulsations had been used for the calculation of the constants. It is noted that only two set of the constants are available, that is, for the polytrope model withn = 3 and a cepheid model without convection. Some examples of nonlinear dynamical effects will be shown.It is shown that if the constants were suitable values, the period-ratio of double-mode cepheids is probably realized. The possibility is briefly suggested.     

Non-adiabatic oscillations in an isothermal atmosphere

The problem of propagation and linear transformation of atmospheric and temperature waves in an isothermal atmosphere with stratified heat exchange has been solved analytically. The cases are discussed where temperature inhomogeneities produced by the waves are either optically thin or optically thick. Formulae have been derived for the absorption, reflection, and transformation of waves when they are transformed from adiabatic into isothermal. The properties of temperature waves are considered. The theory of non-adiabatic atmospheric waves is applied to an analysis of oscillations in a quiet solar atmosphere.     

Mixing-length theory and the excitation of solar acoustic oscillations

The stability of radial solar acoustic oscillations is studied using a time-dependent formulation of mixing-length theory. Though the radiation field is treated somewhat simplistically with the Eddington approximation, and we appreciate that any coupling of the pulsation to the radiation field is important, for the lower frequency radial modes that have been computed this should not produce too serious an error. Instead, we have concentrated upon treating the coupling with convection as accurately as is currently possible with generalized mixing-length theory in order to learn something about its pertinence. Our principal conclusion is that, according to this theory, solar radial acoustic oscillations are expected to be stable and generated by turbulence. Moreover, the theory predicts changes in mode frequency that may, in part, explain the discrepancy between solar observations and the adiabatic pulsation frequencies of theoretical models. We also compute the amplitudes of the modes using a theory of stochastic excitation. These are in good agreement with observed power spectra.     

Linear response of galactic haloes to adiabatic gravitational perturbations

We determine the response of a self-similar isothermal stellar system to small adiabatic gravitational perturbations. For odd spherical harmonics the response is identical to the response of the analogous isothermal fluid system. For even spherical harmonics the response can be regarded as an infinite series of wavetrains in log r , implying alternating compression and rarefaction in equal logarithmic radius intervals. Partly because of the oscillatory nature of the solutions, tidal fields from external sources are not strongly amplified by an intervening isothermal stellar system, except at radii ≲10^−3.5 times the satellite radius; at some radii the stellar system can even screen the external tidal field in a manner analogous to Debye screening. As Weinberg has pointed out, individual resonances in a stellar system can strongly amplify external tidal fields over a limited radial range, but we cannot address this possibility because we examine only adiabatic perturbations. We also discuss the application of our method to the halo response caused by the slow growth of an embedded thin disc.     

Model Atmospheres of AGB Stars: Dynamics, Molecules and Dust

The atmospheres and circumstellar envelopes of AGB stars are characterized by complex physical phenomena like shock waves caused by stellar pulsation or formation of molecules and dust which often lead to a heavy mass loss and have a strong influence on IR properties as observed by ISO. To allow a physical interpretation of various observations we have constructed improved dynamical model atmospheres of long-period variables. In this contribution we mainly investigate the dependence of the atmospheric structure and its variability on stellar pulsation, molecular opacities and time-dependent dust formation. IR spectra resulting from our models are discussed in detail by Loidl et al. (1997b) and compared to ISO-SWS spectra obtained by Hron et al. (1997). This revised version was published online in August 2006 with corrections to the Cover Date.     

Semiclassical approximation for low-degree stellar p modes – I. The classical eigenfrequency equation

A new eigenfrequency equation for low-degree solar-like oscillations in stars is developed, based on the assumption of purely classical propagation in the stellar interior of acoustic waves modified by buoyancy and gravity . Compared with high-frequency asymptotic analysis, the eigenfrequency equation has a new functional form, with expansion in powers of ℓ(ℓ+1) instead of 1/ ω . Basic observable quantities, the 'large' and 'small' frequency separations , are interpreted as the dependence on frequency and refraction angle of a classical action integral for waves propagating close to the stellar diameter. The new eigenfrequency equation gives a significant improvement in accuracy over previous analyses when tested with solar p modes, suggesting this as an alternative and more powerful tool for applications in stellar seismology.     

O-bearing Molecules in Carbon-rich Proto-Planetary Objects

We use high-precision multiband photometric data of the first-overtone RR Lyrae star U Comae to investigate the predictive capability of full-amplitude, nonlinear, convective hydrodynamical models. The main outcome of this investigation is that theoretical predictions properly account for the luminosity variations along a full pulsation cycle. Moreover, we find that this approach, because of the strong dependence of this observable and of the pulsation period on stellar parameters, supplies tight constraints on stellar mass, effective temperature, and distance modulus. Pulsational estimates of these parameters appear in good agreement with empirical ones. Finally, a well-defined bump just before the luminosity maximum gave the unique opportunity to calibrate the turbulent convection model adopted for handling the coupling between pulsation and convection.     

Non-adiabatic discrete Alfvén waves in coronal loops and prominences

Discrete Alfvén waves in coronal loops and prominences are investigated in non-ideal magnetohydrodynamics. The non-ideal effects included are anisotropic, thermal conduction, and optically thin radiation. The classic ideal Alfvén continuum is not altered by these non-ideal effects, but the discrete Alfvén modes, which exist under certain conditions above or below the Alfvén continuum in ideal MHD, are shown to be influenced by non-adiabatic effects.The existence of discrete, non-adiabatic Alfvén waves, and their damping and overstability are examined for 1D cylindrical equilibrium states with twisted magnetic fields. First, analytic results are obtained for modes of high radial order by means of a WKB-analysis. The subspectrum of discrete Alfvén modes is computed with a numerical code, with particular emphasis on the modes of low radial order. The results show that discrete Alfvén waves are of potential importance for solar applications and also that the information obtained with the WKB-analysis is of limited use in this context.Research Assistant of the Belgian National for Scientific Research.     

Invariant mass approach to the spiral structure of galaxies

Following an earlier suggestion that spiral density waves in galaxies primarily are a response phenomenon to the increase with time of a central quadrupole moment. The object of this note is to describe a simple procedure adequate for estimating the shape of the density response in a stellar disk for given structures of the background and perturbative gravitational fields, respectively, under the premises of first order epicyclic and adiabatic approximation.     

Pulsar radiation belts and transient radio emission

It is proposed that radiation belts similar to the ones in the planetary magnetosphere can exist for a pulsar with a relatively long period and a strong magnetic field. In the belts located in the closed field line region near the light cylinder relativistic pairs are trapped and maintained at a density substantially higher than the local Goldreich–Julian corotation density. The trapped plasma can be supplied and replenished by either direct injection of relativistic pairs from acceleration of externally supplied particles in a dormant outer gap or in situ ionization of the accreted neutral material in the trapping region. The radiation belts can be disrupted by waves that are excited in the region as the result of plasma instabilities or emitted from the surface due to starquakes or stellar oscillations. The disruption can cause an intermittent particle precipitation towards the star producing radio bursts. It is suggested that such bursts may be seen as rotating radio transients.     

Solar models with low opacity

Solar Physics - Evolutionary models of the present Sun with standard and artifically low opacity of stellar matter are obtained and adiabatic nonradial oscillations of the models are computed. The...     

Reflection of Alfvén waves by non-uniform ionospheres

Magnetospheric Alfvén waves are reflected by the ionosphere. We investigate the effect of horizontally varying ionospheric conductivity on the process of Alfvén wave reflection. Four idealised ionospheric models are considered in detail. We find that the reflection process is strongly dependent on the orientation of the wave electric field vector with respect to the boundary between high and low conductivities, and under certain conditions subsidiary Alfvén waves are generated. The field-aligned currents in the subsidiary Alfvén waves serve to close divergent horizontal currents resulting from the non-uniform ionospheric conductivity. The implications for ground-based pulsation studies are discussed.     

Darkening and visibility functions for the global five-minute oscillations

A theory for the brightness fluctuations of the Sun as a star under the effect of its global oscillations has been developed. Formulas for the darkening and visibility of p-modes are derived and their calculations are performed in the local approximation for adiabatic oscillations. Observational data from the DIFOS multichannel photometer onboard the CORONAS-F satellite are used to solve the inverse problem of determining the amplitude of the five-minute temperature fluctuations in the solar photosphere as a function of the height. Analysis of the solution and comparison with the results of other authors suggest that the predicted temperature waves resulting from a linear transformation of p-modes in the photosphere exist in the photosphere. The wavelength and phase velocity of the temperature waves are considerably smaller than those of acoustic waves. It turns out that the solar brightness fluctuations should be produced mainly by the temperature waves in the photosphere, not by the p-modes themselves. The darkening function for the brightness fluctuations is oscillatory in behavior, while the visibility function can differ markedly from that for the Doppler shifts of spectral lines produced by p-modes. These properties are important for interpreting the observations of stellar oscillations based on stellar brightness fluctuations.