String loop vibration around Schwarzschild black hole

String loop vibrations in a central plane of a Schwarzschild black hole are investigated for various string equations of state. We discuss string loop stability and derive frequencies of vibrational modes. Using the vibrating string loop model we ﬁt the quasi-periodic oscillation (QPO) observed in X-ray signal coming from some compact sources. We demonstrate how the string-loop parameters are related to the radial and vertical fundamental vibration modes, and how the vibrational instability can be related to the Q-factor characterizing the observed QPOs

Long-range forces provided by gravitational and electromagnetic (EM) interactions play a vital role in comprehending the proper understanding of high-energy processes around black holes (BHs) [30].Observational evidence strongly supports the existence of magnetic fields (MFs) in the vicinity of BHs, emphasizing their significance in these processes [7,8].Orders of magnitude of MFs around BHs may vary from a few Gs up to 10 8 Gs and more, depending on the source generating the field.For stellar-mass BHs observed in X-ray binaries, the characteristic strength of MFs is of order 10 8 Gs, while for supermassive BHs (SMBHs), the characteristic strength is of order 10 4 Gs.In realistic astrophysical scenarios, the energy densities of MFs at such magnitudes are insufficient to significantly influence the geometry of the background spacetime -the spacetime metric around a BH can be accurately described by the Kerr or Schwarzschild solutions of the Einstein field equations.
In the vicinity of BH where the gravitational and electromagnetic field plays an important role, relativistic magnetohydrodynamics (GRMHD) should be used.It was shown by V. S. Semenov and his colleagues [24] that there is formal equivalence between GRMHD equations and equations for dynamics of relativistic string.Originally, the cosmic strings were introduced as remnants of the phase transitions in the very early universe [32] or strings represented as superconducting vortices were considered in [33].In the case of the cosmic strings general relativistic effect of "conical reduction" of the spacetime occurs -behavior of such strings in relation to black holes was studied e.g. by Larsen [17].For the MHD inspired string loops considered in the present paper, the GR effects are irrelevant regarding their internal structure but influenced significantly the external conditions.However, the current-carrying string loops could represent also plasma exhibiting a string-like behavior due to dynamics of the magnetic field lines [23], or due to the thin flux tubes of magnetized plasma simply described as 1D strings [24].Within this project, we would like to explore connections between complex magnetic flux tube evolution (1D plasma dynamic) and a more simple string loop model.
In this article, we are studying circular string structures (current-carrying with scalar field, "dust" particle loop) threaded onto the symmetry axis of the BH, see Fig. 1.The string loop can oscillate, changing its radius in the plane perpendicular to the symmetry axis, while propagating along the symmetry axis [13].The axisymmetric string loops are governed by their tension and angular momentum.The tension of the string loops prevents their expansion beyond some radius, while their worldsheet current introduces an angular momentum preventing them from collapsing.String loop vibration traveling along the loop can also generate angular moment on the sting loop even for simple Nambu-Goto strings (without a scalar field).It has been shown that frequencies of the string loop quasi-harmonic oscillations (QPOs) can fit twin QPOs observed with resonant frequency ratio 3 : 2 in three Galactic microquasars [28].Resonances and higher frequency string loop vibrational modes and their relation to observed QPOs will be also explored.The oscillating string loops can have frequencies similar to frequencies of epicyclic geodesic motion [29], or epicyclic oscillatory motion of charged particles around magnetized black holes [14].
Based on work [19] we examine string loop vibrations around an equilibrium configuration and its stability against linear perturbations around a Schwarzschild black hole.Due to its relative simplicity, we have chosen the Schwarzschild black hole as just a starting point to a more general Kerr case.A lot of observed black holes are not rotating rapidly and Schwarzschild solution can be used as a good approximation.The string loop can orbit around the black hole without influencing the black hole rotation.String loop is a testing object only, without any direct link to make the black hole rotate.
As established in [19], the string loop around a Schwarzschild black hole is always stable under the polar perturbations, while it is altogether unstable under the equatorial perturbations, due to the instability of the first radial mode.We concentrate on some points that were not discussed in [19], namely on studying equatorial perturbations for different modes separately, and we show that some higher vibrational modes are still stable for a special choice of string parameters.

A. Elastic law and speeds of sound
The string worldsheet is described by the spacetime coordinates X α (τ, λ) with α = 0, 1, 2, 3, given as functions of two worldsheet coordinates: τ parameterizes string evolution, λ is running along the string.We assume that the parameter λ is the arclength in the string unstretched configuration.The embedding X induces a metric on the worldsheet, where g αβ is the background metric, |a = ∂ /∂a and a, b ∈ {τ, λ}.The dynamic of the string is described by the action Following [19], we characterize the string by two parameters: energy density ρ and pressure p (since the worldsheet is two-dimensional, this pressure is actually a force -tension or compression of the string).
For an elastic string whose internal energy density ρ depends only on its stretching: where |σ| is the factor by which the string is stretched according to an observer comoving with it, the Lagrangian density is (see for instance [2,11]) where h ≡ det(h ab ).Note that the stretching is related to the induced metric by The choice of the so-called elastic law ρ = F (σ 2 ), corresponding to a particular kind of elastic string, defines also string pressure p as The relation between the energy density ρ and the string pressure p, which is called an equation of state, is therefore determined by the internal structure of the elastic string.It can also be described by the following two parameters: the physical speed of sound for longitudinal waves having the same expression as the speed of sound for a perfect fluid (see for instance [4,18]), and the speed of sound for transverse waves which generalizes the well-known classical result.The speeds of sound c and s should be real.Therefore, are necessary conditions for the stability of the stretched string, as noted in [3] (otherwise there would exist exponentially growing modes in the limit of small wavelengths).
It is worth mentioning cases for string composition given by different equations of state (elastic law).For example, there is a "dust" string, i.e. ring created of test particles on a circular orbit.Such "dust string" is completely without tension, particles are not interacting, and hence c = s = 0. On the other hand, there is an extremely stiff current-carrying string [10,13] for which we have c = 1, s = 0. We have also a very famous Nambu-Goto string model [32] for which c 2 = −1, s = 1.

B. Equilibrium conditions
The line element describing the spacetime of nonrotating Schwarzschild BH with mass M is given as (10) where the function f (r) takes the form FIG. 2. Graphical representation of radial and vertical modes for perturbed circular string loop.Famous ringworld instability can be seen as 1st radial mode (j = 1).
As in [19], we consider an axially symmetric elastic string loop rotating in equilibrium at a constant radius r in the equatorial plane with angular velocity ω, described by an embedding where r, ω and k are constants and λ ∈ 0, 2π k .Here where r 0 is the radius of the unstretched string loop in flat spacetime, according to the definition of λ.Note that the rotating string loop radius relates to the unstretched radius by the combined effects of stretching and length contraction as where the velocity v of the string with respect to the static observer is given by The equations of motion for the rotating string loop in the considered case lead to equation [19] relating the four unknown quantities v, ω, r and σ, with r 0 being a known and fixed parameter of the loop.

III. STRING LOOP VIBRATIONS AND THEIR STABILITY
We now assume that the string loop has well-defined longitudinal and transverse speeds of sound c and s (which is a necessary condition for stability, as noted in [3]), herewith c, s ≤ 1, and also satisfies the weak energy condition In particular, the weak energy condition ( 17) is satisfied by such string models as "warm" cosmic string model, non-prestressed strings with constant longitudinal speed of sound (ranging from the so-called rigid string to the incoherent, or dust, string) and strings with constant transverse speed of sound (including the Nambu-Goto string).Condition ( 17) can be violated when some exotic matter is included but this could lead to nonphysical sound velocities.We first introduce, following [19], the decomposition of general vibrations into vibrational modes, and later we treat the vibrational frequencies and stability of vibrational modes.

A. Vibrational modes
Analysis of the linear stability of the rotating string loops around an equilibrium configuration (satisfying Eq. ( 16)) leads to the linearized equations of motion [19] A δr F δr + G δϕ ′ + H δr ′′ + I δ φ + J δ ṙ′ + L δr = 0 , (20) with the coefficients depending on r, M , c, s and given explicitly in Appendix A. Since the string loop is closed, the perturbations can be Fourier expanded.Looking for solutions proportional to e iΩτ , one obtains that the frequencies Ω of the j-th mode are the roots of the corresponding characteristic polynomial: for the equatorial perturbations, they can be found from the equation [19] where while for the polar perturbations, they can be found from the equation [19] where The embedding ( 12) is linearly stable under equatorial (polar) perturbations only if, for each mode j ∈ Z, all the roots of pj (q j ) have a non-negative imaginary part.Since complex roots come in conjugate pairs, this necessary condition amounts to requiring all the roots to be real.
For the fundamental (j = 0) mode characteristic equa- tion for the equatorial perturbations takes the form and for the polar perturbations take the form Let us first look at the polar Eq. ( 27), which is a quadratic reduced one.It has two distinct real roots since the coefficients L and N are positive for r > 3M , that is outside the photon sphere, which is a natural requirement for the string loop to fulfill.
The equatorial equation ( 26) is a bi-quadratic reduced one.It has two coinciding roots Ω = 0 and two roots of the remaining quadratic equation Due to definition of the coefficients, whether these latter two roots are real or complex depends on the relation between longitudinal c and transverse s speeds of sound of the string loop under consideration.

B. Frequencies of fundamental modes
Using the vibration theory, we can find that the frequencies of the equatorial and polar perturbations for the fundamental mode are given correspondingly by In the case of "dust" string loops (with c = 0, s = 0) they read while for the current carrying string loops (with c = 1, These results coincide with the earlier results for the test particle motion (which form the "dust" string) and current-carrying string loops [13].
There exists an innermost stable circular loop (ISCL) position, which can be calculated from the Ω 2 r = 0 condition.The smallest circular sting loop radius can be found t see Fig. 4. Interestingly similar value can be obtained for charged particle ISCO around magnetized BH [14].

C. Higher vibration modes
For the higher modes j ≥ 1 the characteristic equations have a more general form, containing all the powers of the variable Ω. Equation ( 22) for the equatorial perturbations represents a characteristic polynomial of degree 4. The number of the real roots of this equation can be inferred due to its discriminant and coefficients.In particular, if the discriminant is positive, there are four distinct roots (all real or all complex), if the discriminant is zero, there is at least one multiple root, and if the discriminant is negative, the equation has two distinct real roots and two complex roots.This means that to have four distinct real roots we need a positive discriminant.Even if we know only one real root and the discriminant is positive, we conclude that the equation has all real distinct roots.Real roots depending on the radius r within the chosen range are plotted in Fig. 6 for different cases of c and s values and the number of the mode j.In the cases where we evidently have three or even all four real roots, as for c = 0.2, s = 0.5, j = 20, r = 20 or c = 0.001, s = 0.999, j = 3, r = 10, the discriminant is positive, as it should be, and the string loop is stable.On the other hand, for c = 0.01, s = 0.01, j = 2, r = 10 or c = 0.999, s = 0.001, j = 1, r = 20, the discriminant is negative, which means that the real roots we see on the plot are the only real roots in this case and the other two roots are complex, which means instability.
The equation for the polar perturbations with characteristic polynomial ( 25) is a quadratic one and has two distinct real roots when its discriminant is positive.Since this discriminant is positive for all the cases depicted in Fig. 6, the absence of the corresponding real roots on the plots is due to the range chosen for the plots.
It is worth mentioning that in [19] the altogether stability of the circular string loop is studied, that is for all the modes at once.In particular, it has been established that the discriminant of the polar characteristic equation is always positive, and therefore the string loop is stable under the polar perturbations for all the cases.Therefore we concentrate our studies on the equatorial perturbations and investigate them for different modes j separately.
In Fig. 6 we plot the real roots of Eq. ( 22) for radial perturbations depending on the radius r (thick curves) with different mode numbers j for some values of the string speeds of sound in between the limiting cases of the dust string (c = s = 0) and current-carrying string (c = 1, s = 0).As can be seen from the second and the third rows of the plots, for the strings with the speeds of sound c = 0.2, s = 0.5 and c = 0.5, s = 0.2 for the mode number j = 2 or j = 3 we already can distinguish three distinct real roots for some values of the radius r.This means stability for these values of the radius, which is also confirmed by the positivity of the discriminant in these cases.For the mode number j = 20 we can see all four distinct real roots, and the range of the stability in r is increasing rapidly as the mode number is growing.The case of the string with the speeds of sound c = 0.001, s = 0.999 presented in the fourth row of the plots is similar: for the mode number j = 2 we already can see three distinct real roots for some interval of the radius r values, and the range of the stability in r is increasing with the mode number growing.
In the near-limiting cases c = 0.01, s = 0.01 and c = 0.999, s = 0.001 (the first and the last rows of the plots) we can see at most two distinct real roots and since the discriminant turns out to be negative for these cases, this means instability up to the mode number j = 20 and even further, with stability emerging after j = 60 for the almost dust string (c = 0.01, s = 0.01) and after j = 300 for the almost current-carrying string (c = 0.999, s = 0.001).
Note that for the first mode, we observe instability in all considered cases of string speeds of sound, since we can see at most two distinct real roots, which corresponds to the negativity of the discriminant in these cases.

D. Vibrations quality factor
The frequencies Ω as the roots of eq. ( 22) can be generally complex.Real parts Re(Ω) are responsible for the mode harmonic oscillations, while imaginary parts Im(Ω) are making the mode unstable -the mode will exponentially grow for Im(Ω) < 0, or decay for Im(Ω) > 0. We can define the parameter giving the quality of oscillations related to its instability using the standard Q-factor [9] which will describe the relationship between the initial energy stored within the resonator and the energy dissipated during a single radian of the oscillation cycle.Such a definition is similar to the ratio of a resonator's central frequency ν 0 to its bandwidth δ, which describes how narrowly peaked the signal is, see Fig. 5.
The string loop vibration problem is reduced to the eigenvalue problem and as such is similar to the problem of finding the quasi-normal modes of the compact objects.In both cases, the solutions to the problem form a discrete infinity of complex eigenfrequencies, labeled by the mode number j.The real part of the solution represents the frequency of the oscillations and the imaginary part is responsible for the damping rate of the signal.Obviously in the case when the imaginary part causes an exponential growth of the corresponding eigenfunction, the object under consideration turns out to be unstable under perturbations with a particular mode number.The quasi-normal ringing for different black hole and wormhole spacetimes has been studied by some of us in [5,6].Here we concentrate on the possible role of string loop oscillatory motion for another astrophysically very important case of oscillations observed in accretion disk orbiting BHs or neutron stars [29].

IV. FITTING HF QPOS
The presented string loop model is motivated by a simplified model for magnetic flux tube dynamic around BH [24].Such linear plasma structures are prone to various instabilities, as we have seen in our string loop model as well.Moreover, it is interesting that there is a connection between current carrying string loop vibration and frequency observed for particle motion in strong magnetic field [14].It is then attractive to explore string loop vibrations and apply the calculated frequencies to QPOs observed in microquasars.
Dimensional harmonic oscillation frequencies Ω β , related to the distant static observer, should be expressed in the physical units using conversion factor c 3 /GM .Then the frequencies of the string loop radial and latitudinal harmonic oscillations ν r , ν θ measured by the distant static observers are given by  I. Observed twin HF QPO data for the three microquasars, and the restrictions on mass and spin of the black holes located in them, based on measurements independent of the HF QPO measurements given by the spectral continuum fitting [22,25].
The quasi-harmonic character of the string loop dynamics trapped in a toroidal region around a given radius in the black hole equatorial plane suggests an interesting astrophysical application related to the HF QPOs observed in the LMXB systems containing a black hole in binary systems, or in active galactic nuclei.Some of the HF QPOs come in pairs of the upper and lower frequencies (ν U , ν L ) of twin peaks in the Fourier power spectra.Since the peaks of high frequencies are close to the orbital frequency of the marginally stable circular orbit representing the inner edge of Keplerian discs orbiting black holes (or neutron stars), the strong gravity effects must be relevant in explaining HF QPOs [31].Here we focus on the case of well-known microquasars, summarized in Table I.
The twin peak HF QPOs and their the 3 : 2 ratio has been pointed first out in [12] suggested that these QPOs should have rational ratios, because of the resonances in oscillations of nearly Keplerian accretion disks; see also [1].It seems that the resonance hypothesis is now well supported by observations and that the 3 : 2 ratio (2ν U = 3ν L ) is seen most often in the twin HF QPOs in the LMXB containing black holes (microquasars).
We will assume applicability of the parametric resonance, discussed in [16], focusing attention to the case of the frequency ratios ν θ : ν r = 3 : 2 or ν θ : ν r = 2 : 3, as the observed values of the twin HF QPO frequencies for GRO 1655-40, XTE 1550-564 and GRS 1915+105 sources show clear ratio for the upper ν U and lower ν L frequencies, see Tab.I. We use the epicyclic resonance QPOs model that is modified to string loop vibrations and identify directly the frequencies ν U , ν L with ν θ , ν r or ν r , ν θ frequencies of the oscillating string.In contrast to the standard resonance epicyclic model based on uncharged particle dynamic, the oscillating string loop around Schwarzschild black hole allows both frequency ratios Since r 3:2 < r 2:3 , see Figs. 7, we call the first resonance radius, where ν θ : ν r = 3 : 2, the inner one, and the second resonance radius, where ν θ : ν r = 2 : 3, the FIG. 6. Frequencies for higher modes of circular string loop vibrations.Thick curves are for radial perturbations Ωr, dashed for vertical one Ω θ ; here we put M = 1.For dust string the higher modes are identical to the 0th mode (explored already in Fig. 3), for with low stiffness (almost dust) we see deviations from 0th mode only for really high modes.For stiff strings, where c → 1 or s → 1 the higher mode behavior has really complicated shape.outer one.Note that for the oscillating string loop the 1 : 1, 1 : 2, 2 : 1, or other, resonant frequency ratios can enter the play.There could be also resonance in just one radial or vertical direction only, but between different individual modes of the vibrating string.
The procedure of fitting the string loop oscillation frequencies to the observed frequencies is presented in Fig. 8, for all the three microquasars GRS 1915+105, XTE 1550-564, and GRO 1655-40.From the restrictions on the spacetime mass parameter M for each of the sources, see Tab.I., we obtain simultaneous restrictions on string loop parameters c, s.As we can see from Fig. 8, the string loop should be sufficiently stiff c, s > 0.5 to fit observed HF QPOs using the epicyclic resonance model.So far we have not considered the source rotation that can modify the radial profiles of the vertical and horizontal frequencies -the effect of the black hole rotation is going to be studied in future work, but our preliminary results show that the fitting can be done even for the extremely fast rotating microquasar GRS 1915+105.
There are also observational restrictions on QPOs Qfactor, i.e. the harmonic oscillation should last long enough to have a strong peak in frequency spectra.As pointed out in [9], the Q-factor for BH QPOs should be typically between 5 and 10 or higher.The instability of the string loop 1st radial mode (ring-world instability) restricts harmonic oscillations and Q-factor dramatically.As we already discussed in the previous section, string loop higher modes have generally complex frequency and the related Q-factor is low (Q ∼ 2.5).Although there are some stable higher modes, the 1st radial mode is always unstable, and only for c, s < 0.01 we do have the Q-factor strong enough to fit observed QPOs.Parameter c, s ∼ 0.01 can give us restrictions on maximal plasma (magnetic flux tube) stiffness within our string loop model.
Unfortunately for c, s ∼ 0.01 the string loop vibrational frequencies are unable to fit observed HF QPOs using the presented simple epicyclic resonance model.More complex sting loop QPOs model with different frequency combination should be used or one should consider resonances between basic and higher modes in separate radial and vertical directions.

V. CONCLUSIONS
Based on work [19] we examined string loop vibrations around an equilibrium configuration and its stability for different radii around a Schwarzschild black hole in linear approach.While in [19] the altogether stability of the circular string loop was already studied, that is for all the modes at once, here we provided frequency analysis for different modes of vibrations, and also for different string loop composition, in dependence on speeds of sound c, s within the string.As established in [19], the string loop around a Schwarzschild black hole is always stable under the polar perturbations, while it is altogether unstable under the equatorial perturbations, since the 1st radial mode is always unstable (ringworld instability).We study equatorial perturbations for different modes separately and show that some higher vibrational modes are still stable but only for a special choice of c, s.
As we already pointed out, there is a connection between one-dimensional plasma structures (magnetic flux tubes) and relativistic string which is based on the similarity of their equation of motion [24,26].Contrary to the dust particle orbiting the BH freely along a circular geodesic, the string loop circular structure (plasma flux tube) has some internal tension.On the other hand, the thick accretion torus formed around the BH needs some pressure force to maintain its size and if the torus pressure is removed, the torus will deflate into a circular orbit as well.Hence one can see the presented string loop model as another accretion disk model, approached from different side and introducing stiffness (tension) into play.Both thick torus and string loop are prone to radial instabilities which have geometrical reasons in both cases: for the thick torus it is Papaloizou-Pringle instability [20] and for string loop it is ring-world instability [19].
Our article is basically a theoretical study of a structure with tension (string loop) orbiting in the accretion disk around a nonrotating BH.However, the string loop vibration frequencies and stability of different vibrational modes (Q-factor) can be easily confronted with observed QPO signals.A single test particle does not have the strength to produce observed X-ray flux change and some extended object (accretion torus, string loop) or some plasma collective phenomena are needed.Since the observed QPOs have a limited lifespan, we need a mechanism to create the vibration and to let them die -and for this effect the sting loop radial instability is convenient.In this work, we demonstrated, that zero mode frequency of string loop oscillations, and resonances between vertical and radial mode can be used for fitting observed HF QPOs in X-ray microquasars quite well using the epicyclic resonance model.However, the string loops are unstable throughout ringworld instability in the 1st radial mode, with quite a short lifetime for high sound speeds c and s.Larger and more reasonable Q factors (more stable oscillations) can be obtained for strings with low stiffness c, s < 0.1 (dust-like).Observed black hole HF QPOs have Q factor typically between 5 and 10 or higher [9] and such Q factor condition restricts string models to weak dust string only.Unfortunately for such almost dust string the HF QPOs frequency fit using the epicyclic resonance model is not working in the pure Schwarzschild spacetime.For almost dust string, the frequencies are very close to those obtained for test particles, and a more complex oscillation model should be used.We can summarize the results of our string loop model in the Schwarzschild geometry for different string equations of state: when we are able to fit the frequency data, the Q-factor is too low; and when a good Q-factor is obtained, frequencies will not match the QPOs data.Therefore, modifications of the string equations of state are needed, or modification of the background spacetime (Kerr geometry) could be useful.
It is worth noting that both pressure and tension effects can be combined into a single model of thick accretion torus with a toroidal magnetic field as introduced by Komissarov [15] and applied e.g. in [21].One of the results of this study is that when tension is included in accretion disk structure the vibrational frequencies will be increased, see how string loop frequencies are increased with string tension and compare them with free test particle frequencies.This is the opposite effect of thick torus pressure forces which will make frequencies go down as the structure becomes thicker [27].One can conclude that the frequencies of structure vibrations are going up as the object becomes stiffer, and going down when the object becomes softer, which is in accordance with our daily experience.It will be interesting to see if the stiffer accretion torus with a toroidal magnetic field will be more resilient to various instabilities, which are very common in GRMHD simulations of oscillating structures.In some of our future work, we would like to give realistic estimates on string loop speeds of sound c, s using plasma flux tube model [26].Also, the closed string loop evolution for any general shape beyond the linear perturbation model using numerical integration of the equation of motion could be fruitful with application to magnetic flux tube dynamic [24].

FIG. 1 .
FIG. 1. Schematic representation of vibrating circular string loop in equatorial plane of Schwarzschild BH.

FIG. 3 .
FIG.3.Radial profiles of the sting loop fundamental frequencies measured by distant measured Ωr and Ω θ for different string composition.We compare frequencies with two limiting cases: the case of particle motion c = s = 0 ("dust" loop) with a current carrying string ("ccs" loop) c = 1, s = 0 (same as c = 0, s → 1 case).

FIG. 4 .
FIG. 4. Position of last stable circular loop radially perturbed as a function of longitudinal and traverse speed.We see the innermost stable circular orbit (ISCO) for dust string located at rISCO = 6, while current-carrying sting loop we have the innermost stable circular loop (ISCL) at rISCL .= 4.3.

FIG. 5 .
FIG. 5. Signal generated by a harmonic oscillator with complex frequency for Q-factor Q = 250 (dotted) and Q = 2.5 (solid) from the point of time (left) and frequency (right, Log-Log plot) domain.

FIG. 7 .FIG. 8 .
FIG. 7. Frequencies ν θ , νr of small harmonic oscillations measured by static distant observers are given for the string loop around the Schwarzschild black hole.The positions of the resonant radii are also presented.

FIG. 9 .
FIG. 9. Stability of 3:2 and 2:3 resonance for string loop zeroth mode.The Q-factor is calculated from ratio of upper resonant frequency given by stable (real) zeroth mode to the unstable (imaginary) value of first radial mode, see definition (35).The most stable oscillations (highest Q-factor) can be found for small c, s.