We discuss low-frequency g modes excited by resonant couplings with weakly unstable oscillatory convective modes in the rotating convective core in early-type main-sequence stars. Our non-adiabatic pulsation analyses including the effect of Coriolis force for $2\,M_\odot$ main-sequence models show that if the convective core rotates slightly faster than the surrounding radiative layers, g modes in the radiative envelope are excited by a resonance coupling. The frequency of the excited g mode in the inertial frame is close to $|m\Omega_{\rm c}|$ with $m$ and $\Omega_{\rm c}$ being the azimuthal order of the g mode and the rotation frequency of the convective core, respectively. These g mode frequencies are consistent with those of photometric rotational modulations and harmonics observed in many early-type main-sequence stars.

We calculate small amplitude gravitational and thermal tides of uniformly rotating hot Jupiters composed of a nearly isentropic convective core and a geometrically thin radiative envelope. We treat the fluid in the convective core as a viscous fluid and solve linearized Navier Stokes equations to obtain tidal responses of the core, assuming that the Ekman number ${\rm Ek}$ is a constant parameter. In the radiative envelope, we take account of the effects of radiative dissipations on the responses. The properties of tidal responses depend on thermal timescales $\tau_*$ in the envelope and Ekman number Ek in the core and on weather the forcing frequency $\omega$ is in the inertial range or not, where the inertial range is defined by $|\omega|\le2\Omega$ for the rotation frequency $\Omega$. If ${\rm Ek}\gtrsim 10^{-7}$, the viscous dissipation in the core is dominating the thermal contributions in the envelope for $\tau_*\gtrsim 1$ day. If ${\rm Ek}\lesssim 10^{-7}$, however, the viscous dissipation is comparable to or smaller than the thermal contributions and the envelope plays an important role to determine the tidal torques. If the forcing is in the inertial range, frequency resonance of the tidal forcing with core inertial modes significantly affects the tidal torques, producing numerous resonance peaks of the torque. Depending on the sign of the torque in the peaks, we suggest that there exist cases in which the resonance with core inertial modes hampers the process of synchronization between the spin and orbital motion of the planets.

We calculate tidal torque due to semi-diurnal thermal tides in rotating hot Jupiters, taking account of the effects of radiative cooling in the envelope and of the planets rotation on the tidal responses. We use a simple Jovian model composed of a nearly isentropic convective core and a thin radiative envelope. To represent the tidal responses of rotating planets, we employ series expansions in terms of spherical harmonic functions $Y_l^m$ with different $l$s for a given $m$. For low forcing frequency, there occurs frequency resonance between the forcing and the $g$- and $r$-modes in the envelope and inertial modes in the core. We find that the resonance enhances the tidal torque, and that the resonance with the $g$- and $r$-modes produces broad peaks and that with the inertial modes very sharp peaks, depending on the magnitude of the non-adiabatic effects associated with the oscillation modes. We also find that the behavior of the tidal torque as a function of the forcing frequency (or period) is different between prograde and retrograde forcing, particularly for long forcing periods because the $r$-modes, which have long periods, exist only on the retrograde side.

We calculate overstable convective modes of uniformly rotating hot Jupiters, which have a convective core and a thin radiative envelope. Convective modes in rotating planets have complex frequency $\omega$ and are stabilized by rapid rotation so that their growth rates $\propto\omega_{\rm I}={\rm Im}(\omega)$ are much smaller than those for the non-rotating planets. The stabilized convective modes excite low frequency gravity waves in the radiative envelope by frequency resonance between them. We find that the convective modes that excite envelope gravity waves remain unstable even in the presence of non-adiabatic dissipations in the envelope. We calculate the heating rates due to non-adiabatic dissipations of the oscillation energy of the unstable convective modes and find that the magnitudes of the heating rates cannot be large enough to inflate hot Jupiters sufficiently so long as the oscillation amplitudes remain in the linear regime.

We calculate axisymmetric magnetic modes of a neutron star possessing a mixed poloidal and toroidal magnetic field, where the toroidal field is assumed to be proportional to a dimensionless parameter $\zeta_0$. Here, we assume an isentropic structure for the neutron star and consider no effects of rotation. Ignoring the equilibrium deformation due to the magnetic field, we employ a polytrope of the index $n=1$ as the background model for our modal analyses. For the mixed poloidal and toroidal magnetic field with $\zeta_0\not=0$, axisymmetric spheroidal and toroidal modes are coupled. We compute axisymmetric spheroidal and toroidal magnetic modes as a function of the parameter $\zeta_0$ from $0$ to $\sim 1$ for the surface field strengths $B_S=10^{14}$G and $10^{15}$G. We find that the frequency $\omega$ of the magnetic modes decreases with increasing $\zeta_0$. We also find that the frequency of the spheroidal magnetic modes is almost exactly proportional to $B_S$ for $\zeta_0\lesssim 1$ but that this proportionality holds only when $\zeta_0\ll 1$ for the toroidal magnetic modes. The wave patterns of the spheroidal magnetic modes and toroidal magnetic modes are not strongly affected by the coupling so long as $\zeta_0\lesssim 1$. We find no unstable modes having $\omega^2<0$.

Asteroseismic inference from pressure modes (p modes) and buoyancy, or gravity, modes (g modes) is ubiquitous for stars across the Hertzsprung--Russell diagram. Until now, however, discussion of r modes (global Rossby waves) has been rare. Here we derive the expected frequency ranges of r modes in the observational frame by considering the visibility of these modes. We find that the frequencies of r modes of azimuthal order $m$ appear as groups at slightly lower frequency than $m$ times the rotation frequency. Comparing the visibility curves for r modes with Fourier amplitude spectra of Kepler light curves of upper main-sequence B, A and F stars, we find that r modes are present in many $\gamma$ Dor stars (as first discovered by Van Reeth et al. 2016), spotted stars, and so-called Heartbeat stars, which are highly eccentric binary stars. We also find a signature of r modes in a frequently bursting Be star observed by Kepler. In the amplitude spectra of moderately to rapidly rotating $\gamma$ Dor stars, r-mode frequency groups appear at lower frequency than prograde g-mode frequency groups, while in the amplitude spectra of spotted early A to B stars, groups of symmetric (with respect to the equator) r-mode frequencies appear just below the frequency of a structured peak that we suggest represents an approximate stellar rotation rate. In many Heartbeat stars, a group of frequencies can be fitted with symmetric $m=1$ r modes, which can be used to obtain rotation frequencies of these stars.

We calculate axisymmetric spheroidal modes of neutron stars magnetized with an axisymmetric poloidal magnetic field. We use polytropes of the indices $n\sim1$ as background equilibrium models of neutron stars where we ignore the deformation due to the magnetic fields. For a poloidal magnetic field, axisymmetric normal modes of non-rotating stars are decoupled into spheroidal modes and toroidal modes, and we can treat spheroidal modes separately from toroidal modes. For the surface field strength $B_S$ ranging from $10^{14}$G to $10^{16}$G, we calculate axisymmetric spheroidal magnetic modes whose oscillation frequency is proportional to $B_S$. The typical oscillation frequency of the magnetic modes is $\sim 0.01\times\sqrt{GM/R^3}$ for $B_S\sim 10^{15}$G, where $M$ and $R$ are respectively the mass and radius of the star and $G$ is the gravitational constant. For $M=1.4M_\odot$ and $R=10^6$cm, this frequency is $\sim 20$Hz, which may explain low frequency QPOs found for SGR 1806-204 and SGR 1900+14. We also find modes of frequency $> \sqrt{GM/R^3}$ corresponding to the radial fundamental and first harmonic modes. No unstable magnetic modes are found for axisymmetric spheroidal oscillations of magnetized stars.

We calculate tidally driven mean flows in a slowly and uniformly rotating massive main sequence star in a binary system. We treat the tidal potential due to the companion as a small perturbation to the primary star. We compute tidal responses of the primary as forced linear oscillations, as a function of the tidal forcing frequency $\omega_{\rm tide}=2(\Omega_{\rm orb}-\Omega)$, where $\Omega_{\rm orb}$ is the mean orbital angular velocity and $\Omega$ is the angular velocity of rotation of the primary star. The amplitude of the tidal responses is proportional to the parameter $f_0\propto (M_2/M)(a_{\rm orb}/R)^{-3}$, where $M$ and $M_2$ are the masses of the primary and companion stars, $R$ is the radius of the primary and $a_{\rm orb}$ is the mean orbital separation between the stars. For a given $f_0$, the amplitudes depend on $\omega_{\rm tide}$ and become large when $\omega_{\rm tide}$ is in resonance with natural frequencies of the star. Using the tidal responses, we calculate axisymmetric mean flows, assuming that the mean flows are non-oscillatory flows driven via non-linear effects of linear tidal responses. We find that the $\phi$-component of the mean flow velocity dominates. We also find that the amplitudes of the mean flows are large only in the surface layers where non-adiabatic effects are significant and that the amplitudes are confined to the equatorial regions of the star. Depending on $M_2/M$ and $a_{\rm orb}/R$, the amplitudes of mean flows at the surface become significant.

Zonal and meridional axisymmetric flows can deeply impact the rotational and chemical evolution of stars. Therefore, momentum exchanges between waves propagating in stars, differential rotation, and meridional circulation must be carefully evaluated. In this work, we study axisymmetric mean flows in rapidly and initially uniformly rotating massive stars driven by small amplitude non-axisymmetric $\kappa$-driven oscillations. We treat them as perturbations of second-order of the oscillation amplitudes and derive their governing equations as a set of coupled linear ordinary differential equations. This allows us to compute 2-D zonal and meridional mean flows driven by low frequency $g$- and $r$-modes in slowly pulsating B (SPB) stars and $p$-modes in $\beta$ Cephei stars. Oscillation-driven mean flows usually have large amplitudes only in the surface layers. In addition, the kinetic energy of the induced 2-D zonal rotational motions is much larger than that of the meridional motions. In some cases, meridional flows have a complex radial and latitudinal structure. For SPB stars, we find that there are low frequency retrograde $g$-modes and $r$-modes that drive mean flows with positive velocity components in the radial and azimuthal directions at the stellar surface. This suggests that low frequency retrograde modes can transport angular momentum to the surface layers, which may help to form a circumstellar gaseous disc. Moreover, pulsation-driven and rotation-driven meridional flows can have similar amplitudes. These results show the importance of taking wave-- mean flow interactions into account when studying the evolution of massive stars.

We calculate non-axisymmetric oscillations of neutron stars magnetized by purely poloidal magnetic fields. We use polytropes of index $n=1$ and 1.5 as a background model, where we ignore the equilibrium deformation due to the magnetic field. Since separation of variables is not possible for the oscillation of magnetized stars, we employ finite series expansions for the perturbations using spherical harmonic functions. Solving the oscillation equations as the boundary and eigenvalue problem, we find two kinds of discrete magnetic modes, that is, stable (oscillatory) magnetic modes and unstable (monotonically growing) magnetic modes. For isentropic models, the frequency or the growth rate of the magnetic modes is exactly proportional to $B_{\rm S}$, the strength of the field at the surface. The oscillation frequency and the growth rate are affected by the buoyant force in the interior, and the stable stratification tends to stabilize the unstable magnetic modes.

We calculate non-axisymmetric oscillations of uniformly rotating polytropes magnetized with a purely toroidal magnetic field, taking account of the effects of the deformation due to the magnetic field. As for rotation, we consider only the effects of Coriolis force on the oscillation modes, ignoring those of the centrifugal force, that is, of the rotational deformation of the star. Since separation of variables is not possible for the oscillation of rotating magnetized stars, we employ finite series expansions for the perturbations using spherical harmonic functions. We calculate magnetically modified normal modes such as $g$-, $f$-, $p$-, $r$-, and inertial modes. In the lowest order, the frequency shifts produced by the magnetic field scale with the square of the characteristic Alfv\'en frequency. As a measure of the effects of the magnetic field, we calculate the proportionality constant for the frequency shifts for various oscillation modes. We find that the effects of the deformation are significant for high frequency modes such as $f$- and $p$-modes but unimportant for low frequency modes such as $g$-, $r$-, and inertial modes.

We calculate axisymmetric toroidal modes of magnetized neutron stars with a solid crust in the general relativistic Cowling approximation. We assume that the interior of the star is threaded by a poloidal magnetic field, which is continuous at the surface with an outside dipole field. We examine the cases of the field strength $B_{{\rm{S}}}\sim10^{16}$ G at the surface. Since separation of variables is not possible for the oscillations of magnetized stars, we employ finite series expansions for the perturbations using spherical harmonic functions. We find discrete normal toroidal modes of odd parity, but no toroidal modes of even parity are found. The frequencies of the toroidal modes form distinct mode sequences and the frequency in a given mode sequence gradually decreases as the number of radial nodes of the eigenfunction increases. From the frequency spectra computed for neutron stars of different masses, we find that the frequency is almost exactly proportional to $B_{{\rm{S}}}$ and is well represented by a linear function of $R/M$ for a given $B_{{\rm{S}}}$, where $M$ and $R$ are the mass and radius of the star. The toroidal mode frequencies for $B_{{\rm{S}}}\sim 10^{15}$ G are in the frequency range of the quasi-periodic oscillations (QPOs) detected in the soft-gamma-ray repeaters, but we find that the toroidal normal modes cannot explain all the detected QPO frequencies.

We estimate the amount of angular momentum transferred by the low-frequency oscillations detected in the rapidly rotating hot Be star HD 51452. Here, we assume that the oscillations detected are stochastically excited by convective motions in the convective core of the star, that is, we treat the oscillations as forced oscillations excited by the periodic convective motions of the core fluids having the frequencies observationally determined. With the observational amplitudes of the photometric variations, we determine the oscillation amplitudes, which makes it possible to estimate the net amount of angular momentum transferred by the oscillations using the wave-meanflow interaction theory. Since we do not have any information concerning the azimuthal wavenumber $m$ and spherical harmonic degree $l$ for each of the oscillations, we assume that all the frequencies detected are prograde or retrograde in the observer's frame and they are all associated with a single value of $m$ both for even modes ($l=|m|$) and for odd modes ($l=|m|+1$). We estimate the amount of angular momentum transferred by the oscillations for $|m|=1$ and 2, which are typical $|m|$ values for Be stars, and find that the amount is large enough for a decretion disc to form around the star. Therefore, transport of angular momentum by waves stochastically excited in the core of Be stars might be responsible for the Be phenomenon.

We discuss candidates for non-radial modes excited in a mass accreting and rapidly rotating neutron star to explain the coherent frequency identified in the light curves of a millisecond X-ray pulsar XTE J1751-305. The spin frequency of the pulsar is $\nu_{\rm spin}\cong435$Hz and the identified coherent frequency is $\nu_{\rm osc}=0.5727595\times\nu_{\rm spin}$. Assuming the frequency identified is that observed in the corotating frame of the neutron star, we find that the surface $r$-modes of $l^\prime=m=1$ and 2 excited by $\epsilon$-mechanism due to helium burning in the thin shell can give the frequency ratio $\kappa=\nu_{\rm osc}/\nu_{\rm spin}\simeq0.57$ at $\nu_{\rm spin}=435$Hz. As another candidate for the observed ratio $\kappa$, we also suggest a toroidal crustal mode that has penetrating amplitudes in the fluid core and is destabilized by gravitational wave emission. Since the surface fluid layer is separated from the fluid core by a solid crust, the amplitudes of an $r$-mode in the core, which is destabilized by emitting gravitational waves, can be by a large factor different from those in the fluid ocean. We find that the amplification factor defined as $f_{\rm amp}=\alpha_{\rm surface}/\alpha_{\rm core}$ is as large as $f_{\rm amp}\sim 10^2$ for the $l^\prime=m=2$ $r$-mode at $\nu_{\rm spin}=435$Hz for a $M=1.4M_\odot$ neutron star model. Because of this significant amplification of the $r$-mode amplitudes in the surface fluid layer, we suggest that, when proper corrections to the $r$-mode frequency such as due to the general relativistic effects are taken into consideration, the core $r$-mode of $l^\prime=m=2$ can be a candidate for the detected frequency, without leading to serious contradictions to, for example, the spin evolution of the underlying neutron star.

The Be phenomenon, that is the ejection of matter from Be stars into a circumstellar disk, has been a long lasting mystery. In the last few years, the CoRoT (Convection, Rotation and planetary Transits) satellite brought clear evidence that Be outbursts are directly correlated with pulsations. We found that it may be the transport of angular momentum by waves or pulsation modes that brings the already rapid stellar rotation to its critical value at the surface, and allows the star to eject material. The recent discovery of stochastically excited gravito-inertial modes by CoRoT in a hot Be star strengthens this scenario. We present the CoRoT observations and modeling of several Be stars and describe the new picture of the Be phenomenon which arose from these results.

We discuss steady viscous Keplerian decretion discs around rapidly rotating stars. We assume that low frequency modes, which may be excited by the opacity bump mechanism, convective motion in the core, or tidal force if the star is in a binary system, can transport an enough amount of angular momentum to the region close to the stellar surface. Under this assumption, we construct a star-disc system, in which there forms a viscous decretion disc around a rapidly rotating star because of the angular momentum supply. We find a series of solutions of steady viscous decretion discs around a rapidly rotating star that extend to $R_{\rm disc}\gtsim~10R_*$ with $R_*$ being the equatorial radius of the star, depending on the amount of angular momentum supply.

Using weakly non-linear theory of oscillation, we estimate the amplitudes of low frequency modes in a slowly pulsating B (SPB) star, taking account of the effects of rotation on the modes. Applying the formulation by Schenk et al (2002), we compute non-linear coupling coefficient between the low frequency modes and estimate the equilibrium amplitudes of the modes excited in the star, assuming the amplitudes of the unstable modes are saturated as a result of non-linear coupling with stable modes, that is, as a result of parametric instability expected between one unstable mode and two stable modes. We use the traditional approximation to calculate adiabatic and non-adiabatic oscillations in a rotating star. We find $r$-modes in a rapidly rotating star play a significant role in the amplitude determination through non-linear coupling. We also find that for low $m$ modes, the fractional amplitudes of the radiative luminosity caused by the low frequency modes are of order $10^{-4}$ to $10^{-3}$ at the surface.

We have studied the stability of low degree $g$-modes in uniformly rotating B-type stars, taking into account the effects of the Coriolis force and the rotational deformation. From an analysis treating rotation frequency as a small parameter it is found that slow rotation tends to $destabilize$ high radial-order $retrograde$ $g$-modes, although the effect is very small or absent for relatively low order modes. Calculating eigenfrequencies at selected rotation rates, we find, on the other hand, that rapid rotation tends to $stabilize$ $retrograde$ $g$-modes. The stabilizing effect appears stronger for less massive B-type stars having low effective temperatures. If we change rotation rate continuously, the frequency of a $g$-mode belonging to ($l', m$) crosses frequencies of other $g$-modes belonging to ($l', m$). If the parity of the two encountering modes are the same, they interact each other and the stability (i.e., imaginary part of eigenfrequency) of each mode is modified. Using an asymptotic method we discuss the property of such mode crossings and couplings. For rapidly rotating stars mode couplings are important for the stability of low degree $g$-modes. In particular, we find that the stabilization of retrograde $g$-modes in rapidly rotating stars is due to many strong mode couplings, while %prograde sectoral% modes are exceptionally immune to the damping effects from the mode couplings.

We calculate light curves produced by a hot spot of a rapidly rotating neutron star, assuming that the spot is perturbed by a core $r$-mode, which is destabilized by emitting gravitational waves. To calculate light curves, we take account of relativistic effects such as the Doppler boost due to the rapid rotation and light bending assuming the Schwarzschild metric around the neutron star. We assume that the core $r$-modes penetrate to the surface fluid ocean to have sufficiently large amplitudes to disturb the spot. For a $l'=m$ core $r$-mode, the oscillation frequency $\omega\approx2m\Omega/[l'(l'+1)]$ defined in the co-rotating frame of the star will be detected by a distant observer, where $l'$ and $m$ are respectively the spherical harmonic degree and the azimuthal wave number of the mode, and $\Omega$ is the spin frequency of the star. In a linear theory of oscillation, using a parameter $A$ we parametrize the mode amplitudes such that ${\rm max}\left(|\xi_\theta|,|\xi_\phi|\right)/R=A$ at the surface, where $\xi_\theta$ and $\xi_\phi$ are the $\theta$ and $\phi$ components of the displacement vector of the mode and $R$ is the radius of the star. For the $l'=m=2$ $r$-mode with $\omega=2\Omega/3$, we find that the fractional Fourier amplitudes at $\omega=2\Omega/3$ in light curves depend on the angular distance $\theta_s$ of the spot centre measured from the rotation axis and become comparable to or even larger than $A\sim0.001$ for small values of $\theta_s$.

We investigate non-axisymmetric low frequency modes of a rotating and magnetized neutron star, assuming that the star is threaded by a dipole magnetic field whose strength at the stellar surface, $B_0$, is less than $\sim 10^{12}$G, and whose magnetic axis is aligned with the rotation axis. For modal analysis, we use a neutron star model composed of a fluid ocean, a solid crust, and a fluid core, where we treat the core as being non-magnetic assuming that the magnetic pressure is much smaller than the gas pressure in the core. Here, we are interested in low frequency modes of a rotating and magnetized neutron star whose oscillation frequencies are similar to those of toroidal crust modes of low spherical harmonic degree and low radial order. For a magnetic field of $B_0\sim 10^7$G, we find Alfv\'en waves in the ocean have similar frequencies to the toroidal crust modes, and we find no $r$-modes confined in the ocean for this strength of the field. We calculate the toroidal crustal modes, the interfacial modes peaking at the crust/core interface, and the core inertial modes and $r$-modes, and all these modes are found to be insensitive to the magnetic field of strength $B_0\ltsim10^{12}$G. We find the displacement vector of the core $l^\prime=|m|$ $r$-modes have large amplitudes around the rotation axis at the stellar surface even in the presence of a surface magnetic field $B_0\sim10^{10}$G, where $l^\prime$ and $m$ are the spherical harmonic degree and the azimuthal wave number of the $r$-modes, respectively. We suggest that millisecond X-ray variations of accretion powered X-ray millisecond pulsars can be used as a probe into the core $r$-modes destabilized by gravitational wave radiation.

We calculate axisymmetric toroidal modes of magnetized neutron stars with a solid crust. We assume the interior of the star is threaded by a poloidal magnetic field that is continuous at the surface with the outside dipole field whose strength $B_p$ at the magnetic pole is $B_p\sim 10^{16}$G. Since separation of variables is not possible for oscillations of magnetized stars, we employ finite series expansions of the perturbations using spherical harmonic functions to represent the angular dependence of the oscillation modes. For $B_p\sim 10^{16}$G, we find distinct mode sequences, in each of which the oscillation frequency of the toroidal mode slowly increases as the number of radial nodes of the eigenfunction increases. The frequency spectrum of the toroidal modes for $B_p\sim 10^{16}$G is largely different from that of the crustal toroidal modes of the non-magnetized model, although the frequency ranges are overlapped each other. This suggests that an interpretation of the observed QPOs based on the magnetic toroidal modes may be possible if the field strength of the star is as strong as $B_p\sim 10^{16}$G.

We calculate axisymmetric oscillations of rotating neutron stars composed of the surface fluid ocean, solid crust, and fluid core, taking account of a dipole magnetic field as strong as $B_S\sim 10^{15}$G at the surface. The adiabatic oscillation equations for the solid crust threaded by a dipole magnetic field are derived in Newtonian dynamics, on the assumption that the axis of rotation is aligned with the magnetic axis so that perturbations on the equilibrium can be represented by series expansions in terms of spherical harmonic functions $Y_l^m(\theta,\phi)$ with different degrees $l$ for a given azimuthal wave number $m$ around the the magnetic axis. Although the three component models can support a rich variety of oscillation modes, axisymmetric ($m=0$) toroidal $_{l}t_n$ and spheroidal $_ls_n$ shear waves propagating in the solid crust are our main concerns, where $l$ and $n$ denote the harmonic degree and the radial order of the modes, respectively. In the absence of rotation, axisymmetric spheroidal and toroidal modes are completely decoupled, and we consider the effects of rotation on the oscillation modes only in the limit of slow rotation. We find that the oscillation frequencies of the fundamental toroidal torsional modes $_{l}t_n$ in the crust are hardly affected by the magnetic field as strong as $B_S\sim 10^{15}$G at the surface. As the radial order $n$ of the shear modes in the crust becomes higher, however, both spheroidal and toroidal modes become susceptible to the magnetic field and their frequencies in general get higher with increasing $B_S$. We also find that the surface $g$ modes and the crust/ocean interfacial modes are suppressed by a strong magnetic field, and that there appear magnetic modes in the presence of a strong magnetic field.

The Microvariability and Oscillations of stars (MOST) satellite has detected low-amplitude light variations ($\Delta m\sim$1 mmag) in the Be star $\beta$ CMi (B8Ve). The observations lasted 41 days and the variations have typical periods $\sim 0.3$ days. We demonstrate that the dominant frequencies are consistent with prograde high-order g-modes of $m=-1$ excited by the Fe-bump of opacity in an intermediate-mass ($\approx 3.5 M_\odot$) star with a nearly critical rotation period of 0.38 days. This is the first detection of nonradial g-mode pulsations in a Be star later than B6 leading to the possibility that pulsations are excited in all classical Be stars.

We examine pulsational stability of low $m$ $r$ modes in SPB stars by calculating fully nonadiabatic oscillations of uniformly rotating stars, where $m$ is an integer representing the azimuthal wave number around the rotation axis. $R$ modes are rotationally induced, non-axisymmetric, oscillation modes, whose oscillation frequency strongly depends on the rotation frequency $\Omega$ of the star. They are conveniently classified by using two integer indices $m$ and $l^\prime\ge |m|$ that define the asymptotic oscillation frequency $2m\Omega/[l^\prime(l^\prime+1)]$ in the limit of $\Omega\to 0$. We find low $m$, high radial order, odd $r$ modes with $l^\prime=m$ in SPB stars are excited by the same iron opacity bump mechanism that excites low frequency $g$ modes of the variables, when the rotation frequency $\Omega$ is sufficiently high. No even $r$ modes with low $m$ are found to be pulsationally unstable. Since the surface pattern of the temperature perturbation of odd modes is antisymmetric about the equator of the star, observed photometric amplitudes caused by the unstable odd $r$ modes with $l^\prime=m$ are strongly dependent on the inclination angle between the axis of rotation and the line of sight. Applying the wave-meanflow interaction formalism to nonadiabatic $r$ modes in rapidly rotating SPB models, we find that because of the $r\phi$ component of the Reynolds stress and the radial transport of the eddy fluctuation of density in the rotating star, the surface rotation is accelerated by the forcing due to the low $l^\prime=m$ unstable $r$ modes.

We calculate light curves produced by $r$-modes with small azimuthal wavenumbers, $m$, propagating in the surface fluid ocean of rotating neutron stars. We include relativistic effects due to rapid rotation, and propagate photons from the stellar surface to a distant observer using the Schwarzschild metric. The wave motions of the surface $r$-modes are confined to the equatorial region of the star, and the surface pattern of the temperature variation can be either symmetric (for even modes) or anti-symmetric (for odd modes) with respect to the equator. Since for the surface $r$-modes the oscillation frequency in the corotating frame of the star is much smaller than the rotation frequency, $\Omega$, we employ the approximation in which the oscillation frequency in the inertial frame, $\sigma$, is given by $\sigma= -m\Omega$. We find that the $even$, $m = 1$ $r$-mode produces the largest light variations. The dominant Fourier component in the light curves of these modes is the fundamental having $\sigma=-\Omega$, and the first harmonic component having $\sigma=-2\Omega$ is always negligible in comparison. The dominant Fourier component of the even, $m=2$ $r$-modes is the first harmonic. Although the $odd$ $r$-modes produce smaller amplitude light variations compared with the $even$ modes, the light curves of the former have a stronger first harmonic component. If both $m=1$ and 2 $r$-modes are excited simultaneously, a rich variety of light curves is possible, including those having an appreciable first harmonic component. We show that the phase difference, $\delta-\delta_E$, between the bolometric light curve and that at a particular photon energy can possibly be used as a probe of the stellar compactness, $R/M$, where $R$ and $M$ are the radius and mass of the star.

We study $r$-modes of a rotating magnetized neutron star, assuming a magnetic dipole field whose axis is aligned with the axis of rotation. We approach the problem by applying a singular perturbation theory to the oscillations of rotating stars. In this treatment, we divide the star into a thin surface magnetic layer and a non-magnetic core. We integrate linearized ideal MHD equations in the surface magnetic layer and non-magnetic oscillation equations in the core, and match the two integrations at the interface to obtain a complete solution. For a polytropic neutron star model of mass $M=1.4M_\odot$ and radius $R=10^6$cm, the magnetic dipole field becomes effective on the modal properties of the $r$-modes only when the field strength $B_S$ is much greater than $10^{14}$G. We also find that the damping effects caused by very short magnetic perturbations in the surface layer are not important for the $r$-mode instability of rapidly rotating neutron stars if the field strength $B_S$ is smaller than $10^{12}$G.

We study the $r$-modes propagating in steadily mass accreting, nuclear burning, and geometrically thin envelopes on the surface of rotating neutron stars. For the modal analysis, we construct the envelope models which are fully radiaitive or have a convective region. As the angular rotation frequency $\Omega$ is increased, the oscillation frequency $\omega$ of the $r$-modes in the thin envelopes deviates appreciably from the asymptotic frequency $\omega=2m\Omega/l^\prime(l^\prime+1)$ defined in the limit of $\Omega\to 0$, where $\omega$ is the frequency observed in the corotating frame of the star, and $m$ and $l^\prime$ are the indices of the spherical harmonic function $Y_{l^\prime}^m$ representing the angular dependence of the modes. We find that the fundamental $r$-modes in the convective models are destabilized by strong nuclear burning in the convective region. Because of excessive heating by nuclear buring, the corotating-frame oscillation frequency $\omega$ of the $r$-modes in the convective models becomes larger, and hence the inertial-frame oscillation frequency $|\sigma|$ becomes smaller, than those of the corresopnding $r$-modes in the radiative models, where $\sigma=\omega-m\Omega$ is negative for the $r$-modes of positive $m$. We find that the relative frequency change $f=-(\sigma_{conv}-\sigma_{rad})/\sigma_{rad}$ is always positive and becomes less than $\sim$0.01 for the fundamental $r$-modes of $l^\prime>|m|+1$ at $|\sigma_{rad}|/2\pi\sim$300Hz for $m=1$ or at $|\sigma_{rad}|/2\pi\sim$600Hz for $m=2$, where $\sigma_{conv}$ and $\sigma_{rad}$ denote the oscillation frequencies for the convective and the radiative envelope models, respectively.

We investigate the modal properties of inertial modes of rotating neutron stars with the core filled with neutron and proton superfluids, taking account of entrainment effects between the superfluids. In this paper, the entrainment effects are modeled by introducing a parameter $\eta$ so that no entrainment state is realized at $\eta=0$. We find that inertial modes of rotating neutron stars with the superfluid core are split into two families, which we call ordinary fluid inertial modes ($i^o$-mode) and superfluid inertial modes ($i^s$-mode). The two superfluids in the core counter-move for the $i^s$-modes. For the $i^o$-modes, $\kappa_0=\lim_{\Omega\to 0}\omega/\Omega$ is only weakly dependent on the entrainment parameter $\eta$, where $\Omega$ and $\omega$ are the angular frequency of rotation and the oscillation frequency observed in the corotating frame of the star, respectively. For the $i^s$-modes, on the other hand, $|\kappa_0|$ almost linearly increases as $\eta$ increases. Avoided crossings as functions of $\eta$ are therefore quite common between $i^o$- and $i^s$-modes. We find that some of the $i^s$-modes that are unstable against the gravitational radiation reaction at $\eta=0$ become stable when $\eta$ is larger than $\eta_{crit}$, the value of which depends on the mode. Since the radiation driven instability associated with the current multipole radiation is quite weak for the inertial modes and the mutual friction damping in the superfluid core is strong, the instability caused by the inertial modes will be easily suppressed unless the entrainment parameter $\eta$ is extremely small and the mutual friction damping is sufficiently weak.

We discuss the modal properties of the $r$-modes of relativistic superfluid neutron stars, taking account of the entrainment effects between superfluids. In this paper, the neutron stars are assumed to be filled with neutron and proton superfluids and the strength of the entrainment effects between the superfluids are represented by a single parameter $\eta$. We find that the basic properties of the $r$-modes in a relativistic superfluid star are very similar to those found for a Newtonian superfluid star. The $r$-modes of a relativistic superfluid star are split into two families, ordinary fluid-like $r$-modes ($r^o$-mode) and superfluid-like $r$-modes ($r^s$-mode). The two superfluids counter-move for the $r^s$-modes, while they co-move for the $r^o$-modes. For the $r^o$-modes, the quantity $\kappa\equiv\sigma/\Omega+m$ is almost independent of the entrainment parameter $\eta$, where $m$ and $\sigma$ are the azimuthal wave number and the oscillation frequency observed by an inertial observer at spatial infinity, respectively. For the $r^s$-modes, on the other hand, $\kappa$ almost linearly increases with increasing $\eta$. It is also found that the radiation driven instability due to the $r^s$-modes is much weaker than that of the $r^o$-modes because the matter current associated with the axial parity perturbations almost completely vanishes.

We investigate the modal properties of the $r$-modes of rotating neutron stars with the core filled with neutron and proton superfluids, taking account of entrainment effects between the superfluids. The stability of the $r$-modes against gravitational radiation reaction is also examined considering viscous dissipation due to shear and a damping mechanism called mutual friction between the superfluids in the core. We find the $r$-modes in the superfluid core are split into ordinary $r$-modes and superfluid $r$-modes, which we call, respectively, $r^o$- and $r^s$-modes. The two superfluids in the core flow together for the $r^o$-modes, while they counter-move for the $r^s$-modes. For the $r^o$-modes, the coefficient $\kappa_0\equiv\lim_{\Omega\to 0}\omega/\Omega$ is equal to $2m/[l^\prime(l^\prime+1)]$, almost independent of the parameter $\eta$ that parameterizes the entrainment effects between the superfluids, where $\Omega$ is the angular frequency of rotation, $\omega$ the oscillation frequency observed in the corotating frame of the star, and $l^\prime$ and $m$ are the indices of the spherical harmonic function representing the angular dependence of the $r$-modes. For the $r^s$-modes, on the other hand, $\kappa_0$ is equal to $2m/[l^\prime(l^\prime+1)]$ at $\eta=0$ (no entrainment), and it almost linearly increases as $\eta$ is increased from $\eta=0$. The mutual friction in the superfluid core is found ineffective to stabilize the $r$-mode instability caused by the $r^o$-mode except in a few narrow regions of $\eta$. The $r$-mode instability caused by the $r^s$-modes, on the other hand, is extremely weak and easily damped by dissipative processes in the star.

Nonradial oscillations of relativistic neutron stars with a solid crust are computed in the relativistic Cowling approximation, in which all metric perturbations are ignored. For the modal analysis, we employ three-component relativistic neutron star models with a solid crust, a fluid core, and a fluid ocean. As a measure for the relativistic effects on the oscillation modes, we calculate the relative frequency difference defined as $\delta\sigma/\sigma\equiv(\sigma_{GR}-\sigma_N)/\sigma_{GR}$, where $\sigma_{GR}$ and $\sigma_R$ are, respectively, the relativistic and the Newtonian oscillation frequencies. The relative difference $\delta\sigma/\sigma$ takes various values for different oscillation modes of the neutron star model, and the value of $\delta\sigma/\sigma$ for a given mode depends on the physical properties of the models. We find that $|\delta\sigma/\sigma|$ is less than $\sim0.1$ for most of the oscillation modes we calculate, although there are a few exceptions such as the fundamental (nodeless) toroidal torsional modes in the crust, the surface gravity modes confined in the surface ocean, and the core gravity modes trapped in the fluid core. We also find that the modal properties, represented by the eigenfunctions, are not strongly affected by introducing general relativity. It is however shown that the mode characters of the two interfacial modes, associated with the core/crust and crust/ocean interfaces, have been interchanged between the two through an avoided crossing when we move from Newtonian dynamics to general relativistic dynamics.

This paper has been withdrawn by the authors, due to inappropriate discussions in Discussions.

We investigate the properties of relativistic $r$-modes of slowly rotating neutron stars by using a relativistic version of the Cowling approximation. In our formalism, we take into account the influence of the Coriolis like force on the stellar oscillations, but ignore the effects of the centrifugal like force. For three neutron star models, we calculated the fundamental $r$-modes with $l'=m=2$ and 3. We found that the oscillation frequency $\bar\sigma$ of the fundamental $r$-mode is in a good approximation given by $\bar\sigma\approx \kappa_0 \Omega$, where $\bar\sigma$ is defined in the corotating frame at the spatial infinity, and $\Omega$ is the angular frequency of rotation of the star. The proportional coefficient $\kappa_0$ is only weakly dependent on $\Omega$, but it strongly depends on the relativistic parameter $GM/c^2R$, where $M$ and $R$ are the mass and the radius of the star. All the fundamental $r$-modes with $l'=m$ computed in this study are discrete modes with distinct regular eigenfunctions, and they all fall in the continuous part of the frequency spectrum associated with Kojima's equation (Kojima 1998). These relativistic $r$-modes are obtained by including the effects of rotation higher than the first order of $\Omega$ so that the buoyant force plays a role, the situation of which is quite similar to that for the Newtonian $r$-modes.

We investigate the properties of $r$-mode oscillations of a slowly rotating neutron star with a solid crust, by taking account of the effects of the Coriolis force. For the modal analysis we employ three-component neutron star models that are composed of a fluid core, a solid crust and a surface fluid ocean. For the three-component models, we find that there exist two kinds of $r$-modes, that is, those confined in the surface fluid ocean and those confined in the fluid core, which are most important for the $r$-mode instability. The $r$-modes do not have any appreciable amplitudes in the solid crust if rotation rate of the star is sufficiently small. We find that the core $r$-modes are strongly affected by mode coupling with the crustal torsional (toroidal) modes and lose their simple properties of the eigenfunction and eigenfrequency as functions of the angular rotation velocity $\Omega$. This indicates that the extrapolation formula, which is obtained in the limit of $\Omega\to 0$, cannot be used to examine the $r$-mode instability of rapidly rotating neutron stars with a solid crust unless the effects of mode coupling with the crustal torsional modes are correctly taken into account.

We investigate the properties of r-mode and inertial mode of slowly rotating, non-isentropic, Newtonian stars, by taking account of the effects of the Coriolis force and the centrifugal force. For the non-isentropic models we consider only two cases, that is, the models with the stable fluid stratification in the whole interior and the models that are fully convective. For simplicity we call these two kinds of models "radiative" and "convective" models in this paper. For both cases, we assume the deviation of the models from isentropic structure is small. Examining the dissipation timescales due to the gravitational radiation and several viscous processes for the polytropic neutron star models, we find that the gravitational radiation driven instability of the r-modes remains strong even in the non-isentropic models. Calculating the rotational modes of the radiative models as functions of the angular rotation frequency $\Omega$, we find that the inertial modes are strongly modified by the buoyant force at small $\Omega$, where the buoyant force as a dominant restoring force becomes comparable with or stronger than the Coriolis force. Because of this property we obtain the mode sequences in which the inertial modes at large $\Omega$ are identified as g-modes or the r-modes with l=|m| at small $\Omega$. We also note that as $\Omega$ increases from $\Omega=0$ the retrograde g-modes become retrograde inertial modes, which are unstable against the gravitational radiation reaction.

We investigate inertial mode oscillations of slowly and uniformly rotating, isentropic, Newtonian stars. Inertial mode oscillations are induced by the Coriolis force due to the star's rotation, and their characteristic frequencies are comparable with the rotation frequency $\Omega$ of the star. So called r-mode oscillations form a sub-class of the inertial modes. In this paper, we use the term ``r-modes'' to denote the inertial modes for which the toroidal motion dominates the spheroidal motion, and the term ``inertial modes'' to denote the inertial modes for which the toroidal and spheroidal motions have comparable amplitude to each other. Using the slow rotation approximation consistent up to the order of $\Omega^3$, we study the properties of the inertial modes and r-modes, by taking account of the effect of the rotational deformation of the equilibrium on the eigenfrequencies and eigenfunctions. The eigenfrequencies of the r-modes and inertial modes calculated in this paper are in excellent agreement with those obtained by Lindblom et al (1999) and Lockitch & Friedman (1998). We also estimate the dissipation timescales due to the gravitational radiation and several viscous processes for polytropic neutron star models. We find that for the inertial modes, the mass multipole gravitational radiation dominates the current multipole radiation, which is dominating in the case of the r-modes. It is also found that the inertial mode instability is more unstable than previously reported by Lockitch & Friedman (1998), and survives the viscous damping processes relevant in neutron stars.