in this paper we introduce a general lagrangian framework to study the dissipative properties of two component systems composed of a high - loss and lossless components which can have gyroscopic properties .
this framework covers any linear lagrangian system provided it has a finite number of degrees of freedom , a nonnegative hamiltonian , and losses accounted by the rayleigh dissipative function , ( * ? ? ?
10.11 , 10.12 ) , ( * ? ? ?
8 , 9 , 46 ) .
such physical systems include , in particular , many different types of damped mechanical systems or electric networks .
[ [ motivation . ] ] motivation .
+ + + + + + + + + + +    we are  looking to design and study two - component composite dielectric media consisting of a high - loss and low - loss component in which the lossy component has a useful property ( or functionality ) such as magnetism .
we want to understand what is the trade - off between the losses and useful properties inherited by the composite from its components .
our motivation comes from a major problem in the design of such structures where a component which carries a useful property , e.g. , magnetism , has prohibitively strong losses in the frequency range of interest .
often this precludes the use of such a lossy component with otherwise excellent physically desirable properties .
then the question stands : is it possible to design a composite system having a useful property at a level comparable to that of its lossy component but with significantly reduced losses over a broad frequency range ?
an important and guiding example of a two - component dielectric medium composed of a high - loss and lossless components was constructed in @xcite .
the example was a simple layered structure which had magnetic properties comparable with a natural bulk material but with 100 times lesser losses in a wide frequency range .
this example demonstrated that it is possible to design a composite material / system which can have a desired property comparable with a naturally occurring bulk substance but with significantly reduced losses .    in order to understand the general mechanism for this phenomenon ,
the authors in @xcite considered a general dynamical system as the model for such a medium where the high - loss component of the medium was represented as a significant fraction of the entire system .
we have found that for such a system the losses of the entire structure become small provided that the lossy component is sufficiently lossy .
the general mechanism of this phenomenon , which we proved in @xcite , is the _ modal dichotomy _ when the entire set of modes of the system splits into two distinct classes , _ high - loss _ and _ low - loss modes _ , based on their dissipative properties .
the higher loss modes for a wide range of frequencies contribute very little to losses and this is how the entire structure can be low loss , whereas the useful property is still present . in that work the way we accounted for the presence of a useful property was by simply demanding that the lossy component was a significant fraction of the entire composite structure without explicitly correlating the useful property and the losses .    in this paper
, we continue our studies on the modal dichotomy and overdamping , which began in @xcite , but now we focus are attention on gyroscopic - dissipative systems .
although our results on modal dichotomy apply to the full generality of the dynamical systems considered here , our results on overdamping in this paper will be restricted to systems without gyroscopy .
consideration of overdamping in systems with gyroscopy will be considered in a future work .
[ [ overview - of - results . ] ] overview of results .
+ + + + + + + + + + + + + + + + + + + +    now in order to account in a general form for the physical properties of the two - component composite system with a high - loss and lossless components we introduce a lagrangian framework with dissipation . this framework can be a basis for studying the interplay of physical properties of interest and dissipation .
these general physical properties include in particular gyroscopy ( or gyrotropy ) , which is intimately related to magnetic properties , different symmetries and other phenomena such as overdamping that are not present in typical dynamical systems .    in the lagrangian
setting , the dissipation is taken into account by means of the rayleigh dissipation function ( * ? ? ?
10.11 , 10.12 ) , ( * ? ?
* sec . 8 , 9 , 46 ) . in the case of two - component composite with lossy and lossless components ,
the rayleigh function is assumed to affect only a fraction @xmath0 of the total degrees of freedom of the system . in a rough sense ,
the loss fraction @xmath1 signifies the fraction of the degrees of freedom that are effected by losses .
we now give a concise qualitative description of the main results of this paper which are discussed next and organized under the following three topics : overdamping and selective overdamping in section [ sodan ] ; virial theorem for dissipative systems and equipartition of energy in section [ svir ] ; standard spectral theory vs. krein theory in sections [ smodel ] and [ sspan ] .
* overdamping and selective overdamping*. overdamping is a regime in dissipative lagrangian systems in which some of the systems eigenmodes are overdamped , i.e. , have exactly zero frequency or , in other words , are non - oscillatory when losses are sufficiently large .
the phenomenon of overdamping has been studied thoroughly but only in the case when the entire system is overdamped @xcite , ( * ? ? ?
9 ) , @xcite meaning all its modes are non - oscillatory . the lagrangian framework and the methods we develop in this paper are applicable to the more general case allowing for only a fraction of the modes to be susceptible to overdamping , a phenomenon we call _ selective overdamping_. our analysis of overdamping and this selective overdamping phenomenon is carried out in section [ sodan ] .
our main results are theorems [ tcovd ] , [ tpovd ] , [ cpodr ] , [ thmod ] , and [ tseod ] .
our interest in selective overdamping is threefold .
first , as we prove in this paper the selective overdamping of the system is a rather _ universal phenomenon _ that can occur for a two - component system when the losses of the lossy component are sufficiently large .
moreover , we show that _ the large losses in the lossy component _ cause not only the modal dichotomy , but also _ results in overdamping of all of the high - loss modes while a positive fraction of the low - loss modes remain oscillatory_. second , since mode excitation is efficient only when its frequency matches the frequency of the excitation force , if these high - loss modes go into the overdamping regime having exactly zero frequency these modes can not be excited efficiently and therefore their associated losses are essentially eliminated .
this explains the absorption suppression for systems composed of lossy and lossless components .
third , as we will show , as the losses in the lossy component increase the overdamped high - loss modes are more suppressed while all the low - loss oscillatory modes are more enhanced with increasingly high quality factor .
_ this provides a mechanism for selective enhancement of these high quality factor , low - loss oscillatory modes and selective suppression of the high - loss non - oscillatory modes_.    in fact , we show that _ the fraction of all overdamped modes is exactly the loss fraction _ @xmath1 _ _  which satisfies _ _ @xmath0 , and half of these overdamped modes are the high - loss modes .
thus it is _ exactly the positive fraction _ @xmath2 _ _  of modes which are low - loss oscillatory modes__. in addition , _ these latter modes have extremely high quality factor which increases as the losses in the lossy component increase_. an example of this behavior using the electric circuit in section [ scircuits ] was shown numerically in @xcite .    * virial theorem for dissipative systems and equipartition of energy*. recall that the virial theorem from classical mechanics ( * ? ? ?
* ;  3.4 ) ( which we review in appendix [ apdxvirthm ] ) is about the oscillatory transfer of energy from one form to another and its equipartition . for conservative linear systems
whose lagrangian is the difference between the kinetic and potential energy , e.g. , a spring - mass system or a parallel lc circuit , the virial theorem says that for any time - periodic state of the system , the time - averaged kinetic energy equals the time - averaged potential energy .
this result is known as the equipartition of energy because of the fact that the total energy of the system , i.e. , the hamiltonian , is equal to the sum of the kinetic and potential energy .
it is this result that we generalize in this paper for dissipative lagrangian systems and their damped oscillatory modes .
in particular , we show that _ for any damped oscillatory mode of the system the kinetic energy equals the potential energy _ and so _
there is an equipartition of the system energy for these damped oscillatory modes into kinetic and potential energy_. our precise statements of these results can be found in section [ svir ] in theorem [ tvir ] and corollary [ tvir_eqp ] .
we have also found that the transition to overdamping is characterized by a breakdown of the virial theorem .
more specifically , as the losses of the lossy component increase some of the modes become overdamped with complete ceasing of any oscillations and with breaking of the equality between kinetic and potential energy .
such a _ breaking of the equipartition of the energy can viewed as a dynamical characterization of the overdamped modes ! _
so we now have two different characterizations of overdamping : spectral and dynamical .
the effectiveness of this dynamical characterization is demonstrated in our study of the selective overdamping phenomenon and in proving that a positive fraction of modes always remain oscillatory no matter how large the losses become in the lossy component of the composite system .    * standard spectral theory vs. krein spectral theory*. the study of the eigenmodes of a lagrangian system often relies on the fact the system evolution can be transformed into the hamiltonian form as the first - order linear differential equations .
however , it is not emphasized enough that the hamiltonian setting does not lead directly to the standard spectral theory of self - adjoint or dissipative operators but that only it might be reduced to one in some important and well known cases by a proper transformation as , for instance , in the case of a simple oscillator .
the general spectral theory of hamiltonian systems is known as the _ krein spectral theory _
@xcite , @xcite , @xcite .
this theory is far more complex than the standard spectral theory and it is much harder to apply .
we have found , as discussed in section [ smodel ] , a general transformation that reduces the evolution equation to the standard spectral theory under the condition of positivity of the hamiltonian .
the standard spectral theory is complete , well understood , and allows for an elaborate and effective perturbation theory .
we used that all in our analysis of the modal dichotomy and symmetries of the spectrum in section [ sspan ] as well as the selective overdamping phenomenon in section [ sodan ] .
our main results on spectral symmetries are proposition [ pssym ] and corollary [ cssym ] .
our principal results on the modal dichotomy can be found in theorems [ tmdic ] , [ cpodr ] and corollaries [ cpfsp ] , [ cmdic ] .
in fact , it is important to point out that the key result of this paper that leads to the modal dichotomy which is essential in the analysis of overdamping is proposition [ pevbd ] which relates to certain eigenvalue bounds .
moreover , these bounds and the dichotomy come from some general results that we derive in appendix [ apxebd ] on the perturbation theory for matrices with a large imaginary part .    [ [ organization - of - the - paper . ] ] organization of the paper .
+ + + + + + + + + + + + + + + + + + + + + + + + + +    the rest of the paper is organized as follows .
section [ smodel ] sets up and discusses the lagrangian framework approach and model used in this paper to study the dissipative properties of two - component composite systems with a high - loss and a lossless component . in section [ scircuits ]
we apply the developed approach to an electric circuit showing all key features of the method .
sections [ svir][sodan ] are devoted to a precise formulation of all significant results in the form of theorems , propositions , and so on ( with section [ sproofs ] providing the proofs of these results ) .
more specifically , section [ svir ] discusses our extension of the virial theorem and equipartition of energy from classical mechanics for conservative lagrangian systems to dissipative systems .
section [ sspan ] is on the spectral analysis of the system eigenmodes including spectral symmetries , modal dichotomy , and perturbation theory in the high - loss regime .
section [ sodan ] contains our analysis of the overdamping phenomena including selective overdamping .
appendices [ apxsc ] and [ apxebd ] contain results we need that we believe will be of use more generally in studying dissipative dynamical systems , but especially for composite systems with high - loss and lossless components .
in particular , appendix [ apxebd ] is on the perturbation theory for matrices with a large imaginary part .
finally , appendices [ apenerg ] and [ apdxvirthm ] review the virial theorem from classical mechanics from the energetic point of view .
the lagrangian and hamiltonian frameworks have numerous well known advantages in describing the evolution of physical systems .
these advantages include , in particular , the universality of the mathematical structure , the flexibility in the choice of variables and the incorporation of the symmetries with the corresponding conservations laws .
consequently , it seems quite natural and attractive to integrate dissipation into the lagrangian and hamiltonian frameworks .
we are particular interested in integrating dissipation into systems with gyroscopic features such as dielectric media with magnetic components and electrical networks .
we start with setting up the lagrangian and hamiltonian structures similarly to those in @xcite and @xcite . in particular
, we use matrices to efficiently handle possibly large number of degrees of freedom .
since we are interested in systems evolving linearly the lagrangian @xmath3 is assumed to be a quadratic function ( bilinear form ) of the coordinates @xmath4   _ { r=1}^{n}$ ] ( column vector ) and their time derivatives @xmath5 , that is @xmath6{l}\dot{q}\\ q \end{array } \right ]   ^{\mathrm{t}}m_{\mathrm{l}}\left [ \begin{array } [ c]{l}\dot{q}\\ q \end{array } \right ]   , \qquad m_{\mathrm{l}}=\left [ \begin{array } [ c]{ll}\alpha & \theta\\ \theta^{\mathrm{t } } & -\eta \end{array } \right ]   , \label{dlag1}\ ] ] where @xmath7 denotes the matrix transposition operation , and @xmath8 and @xmath9 are @xmath10-matrices with real - valued entries .
in addition to that , we assume @xmath8 to be symmetric matrices which are positive definite and positive semidefinite , respectively , and we assume @xmath9 is skew - symmetric , that is @xmath11 notice that only the skew - symmetric part of @xmath9 in ( [ dlag1 ] ) matters for the system dynamics because of the form of the euler - lagrange equations ( [ dlag4 ] ) and so without loss of generality we can assume as we have that the matrix @xmath9 is skew - symmetric .
in fact , the symmetric part of the matrix @xmath9 if any alters the lagrangian by a complete time derivative and consequently can be left out as we have done .
the lagrangian ( [ dlag1 ] ) can be written in the energetic form@xmath12 where @xmath13 are interpreted as the kinetic and potential energy , respectively . by hamilton
s principle the dynamics of the system are governed by the euler - lagrange equations @xmath14 which are the following second - order ordinary differential equation ( odes ) @xmath15    the above second - order odes can be turned in the first - order odes with help of the hamiltonian function defined by the lagrangian through the legendre transformation @xmath16 hence,@xmath17 and , consequently @xmath18   = \frac{1}{2}\dot{q}^{t}\alpha\dot{q}+\frac{1}{2}q^{t}\eta q. \label{dlag7}\ ] ] in particular , the hamiltonian is the sum of the kinetic and potential energy , that is@xmath19 we remind that the hamiltonian function @xmath20 is interpreted as the system energy which is a conserved quantity , that is@xmath21 the function @xmath20 defined by ( [ dlag7 ] ) is a quadratic form associated with a matrix @xmath22 through the relation @xmath23{l}p\\ q \end{array } \right ]   ^{\mathrm{t}}m_{\mathrm{h}}\left [ \begin{array } [ c]{l}p\\ q \end{array } \right ]   , \label{dlag8}\ ] ] where @xmath22 is the @xmath24 matrix having the block form @xmath25{ll}\alpha^{-1 } & -\alpha^{-1}\theta\\ -\theta^{\mathrm{t}}\alpha^{-1 } & \theta^{\mathrm{t}}\alpha^{-1}\theta+\eta \end{array } \right ]   = \left [ \begin{array } [ c]{ll}\mathbf{1 } & 0\\ -\theta^{\mathrm{t } } & \mathbf{1}\end{array } \right ]   \left [ \begin{array } [ c]{ll}\alpha^{-1 } & 0\\ 0 & \eta \end{array } \right ]   \left [ \begin{array } [ c]{ll}\mathbf{1 } & -\theta\\ 0 & \mathbf{1}\end{array } \right ]   , \label{dlag7a}\ ] ] where @xmath26 is the identity matrix .
notice that the representations ( [ dlag8 ] ) , ( [ dlag7a ] ) combined with the inequalities ( [ dlag2 ] ) imply@xmath27 the hamiltonian evolution equations then take the form@xmath28{l}p\\ q \end{array } \right ]   = jm_{\mathrm{h}}\left [ \begin{array } [ c]{l}p\\ q \end{array } \right ]   , \quad j=\left [ \begin{array } [ c]{ll}0 & -\mathbf{1}\\ \mathbf{1 } & 0 \end{array } \right ]   \text { , } \label{dlag9}\ ] ] observe that the symplectic matrix @xmath29 satisfies the following relations@xmath30    * the rayleigh dissipation function .
* up to now we have dealt with a conservative system  satisfying the energy conservation ( [ dlag7b ] ) .
let us introduce now dissipative forces using rayleigh s method described in ( * ? ? ?
10.11 , 10.12 ) , ( * ? ? ?
* sec . 8 , 9 , 46 ) .
the _ rayleigh dissipation function _
@xmath31 is defined as a quadratic function of the generalized velocities , namely @xmath32 where the scalar @xmath33 is a dimensionless _ loss parameter _ which we introduce to scale the intensity of dissipation . in particular
, @xmath34 is an @xmath10 symmetric matrix with real - valued entries and , most importantly , is positive semidefinite .
the dissipation is then introduced through the following general euler - lagrange equations of motion with forces@xmath35 where @xmath36 are generalized dissipative forces and @xmath37 is an external force , yielding the following second - order odes @xmath38      now by a linear ( dissipative ) _ lagrangian system _ we mean a system whose state is described by a time - dependent @xmath39 taking values in the hilbert space @xmath40 with the standard inner product @xmath41 ( i.e. , @xmath42 , where @xmath43 denotes the conjugate transpose , i.e. , @xmath44 ) whose dynamics are governed by the odes ( [ radis2a ] ) .
the energy balance equation for any such state , which follows from ( [ radis2a ] ) , is @xmath45 where @xmath46  as in ( [ dlag7c ] ) but now instead of ( [ dlag3_1 ] )  we have@xmath47 we continue to interpret @xmath48 , @xmath49 as the kinetic and potential energies , @xmath50 as the system energy , and @xmath51 as the dissipated power . the term @xmath52 is interpreted as the rate of work done by the force @xmath53 .
the physical significance of the various energetic terms @xmath54 , where @xmath55 , is that for a complex - valued state @xmath56 with real - valued @xmath57 both @xmath58 and @xmath59 are also states representing physical solutions and @xmath60 .
these latter two terms in the sum reduce to the previous definitions of the energetic term @xmath54 above in ( [ dlag3_1 ] ) , ( [ dlag7 ] ) , or ( [ radis1 ] ) for real - valued vector quantities .
we now introduce an important quantity  the _ loss fraction _ @xmath61 .
it is defined as the ratio of the rank of the matrix @xmath34 to the total degrees of freedom @xmath62 of the system@xmath63 from our hypothesis @xmath64 it follows that @xmath65 .
we consider the dissipative lagrangian system to be a model of a two - component composite with a lossy and a lossless component whenever the loss fraction condition @xmath66 is satisfied .
we then associate the range of the operator @xmath34 , i.e. , @xmath67 , with the lossy component of the system and consider the lossy component to be highly lossy when @xmath68 , i.e. , the high - loss regime .
our paper is focused on this case .
* on the eigenmodes and the quality factor .
* to study the dissipative properties of the dissipative lagrangian system ( [ radis2a ] ) in the high - loss regime @xmath68 , one can consider the eigenmodes of the system and their quality factor in this regime .
an eigenmode of this system ( [ radis2a ] ) is defined as a nonzero solutions of the odes ( [ radis2a ] ) with no forcing , i.e. , @xmath69 , having the form @xmath70 , where @xmath71 and @xmath72 are called its frequency and damping factor , respectively . since
the system is dissipative , i.e. , dissipates energy as interpreted from the energy balance equation ( [ radis4 ] ) , the damping factor must satisfy @xmath73 which implies that @xmath74 .
thus , an eigenmode  is a state of the dissipative lagrangian system which is in damped harmonic motion and the motion is oscillatory provided @xmath75 .
an important quantity which characterizes the quality of such a damped oscillation is the quality factor @xmath76 ( i.e. , @xmath77-factor ) of the eigenmode which is defined as the reciprocal of the relative rate of energy dissipation per temporal cycle , i.e. , @xmath78 it follows immediately from this definition and ( [ radis4_1 ] ) that @xmath76 depends only on the complex frequency @xmath79 of the eigenmode , in particular,@xmath80 with the convention @xmath81 if @xmath82 .
now the matrix hamilton form of the euler - lagrange equation involving rayleigh dissipative forces and external force ( [ radis2 ] ) reads@xmath83{ll}\beta r & 0\\ 0 & 0 \end{array } \right ]   \right )   m_{\mathrm{h}}u+\left [ \begin{array } [ c]{l}f\\ 0 \end{array } \right ]   , \quad u=\left [ \begin{array } [ c]{l}p\\ q \end{array } \right ]   .
\label{radis3}\ ] ]    it is important to recognize that the evolution equation ( [ radis3 ] ) are not quite yet of the desired form which is the most suitable for the spectral analysis . advancing general ideas of the hamiltonian treatment of dissipative systems developed in @xcite we can factor the matrix @xmath22 as @xmath84 where the matrix @xmath85 is the block matrix@xmath86{ll}k_{\mathrm{p } } & 0\\ 0 & k_{\mathrm{q}}\end{array } \right ]   \left [ \begin{array } [ c]{ll}\mathbf{1 } & -\theta\\ 0 & \mathbf{1}\end{array } \right ]   = \left [ \begin{array } [ c]{ll}k_{\mathrm{p } } & -k_{\mathrm{p}}\theta\\ 0 & k_{\mathrm{q}}\end{array } \right ]   , \label{radis6}\\ k_{\mathrm{p}}=\sqrt{\alpha}^{-1},\text { } k_{\mathrm{q}}=\sqrt{\eta } \label{radis6a1}\ ] ] which manifestly takes into account the gyroscopic term @xmath9 . here
@xmath87 and @xmath88 denote the unique positive semidefinite square roots of the matrices @xmath89 and @xmath90 , respectively .
in particular , it follows from the properties ( [ dlag2 ] ) and the proof of ( * ? ? ? * , theorem vi.9 )
that @xmath91 , @xmath92 are @xmath10 matrices with real - valued entries with the properties @xmath93    the introduction of the matrix @xmath85 according to @xcite is intimately related to the introduction of force variables@xmath94 consequently , recasting the evolution equation ( [ radis3 ] ) in terms of the force variables @xmath95 yields the desired canonical form , namely @xmath96 and@xmath97{cc}\omega_{\mathrm{p } } & -\mathrm{i}\phi^{\mathrm{t}}\\ \mathrm{i}\phi & 0 \end{array } \right ]   , \text { \ } b = k\left [ \begin{array } [ c]{ll}r & 0\\ 0 & 0 \end{array } \right ]   k^{\mathrm{t}}=\left [ \begin{array } [ c]{ll}\mathsf{\tilde{r } } & 0\\ 0 & 0 \end{array } \right ]   , \text { \ } f=\left [ \begin{array } [ c]{l}k_{\mathrm{p}}f\\ 0 \end{array } \right ]   , \label{radis8a}\\
\omega_{\mathrm{p}}=-\mathrm{i}2k_{\mathrm{p}}\theta k_{\mathrm{p}}^{\mathrm{t}},\text { \ \ } \phi = k_{\mathrm{q}}k_{\mathrm{p}}^{\mathrm{t}},\text { \ \ } \mathsf{\tilde{r}}=k_{\mathrm{p}}rk_{\mathrm{p}}^{\mathrm{t}}. \label{radis8b}\ ] ] observe that the @xmath24 matrices @xmath98 , @xmath99 in block form are hermitian and positive semidefinite , respectively , that is @xmath100 moreover , the matrix @xmath99 does not have full rank and its rank is that of @xmath34 , i.e. @xmath101 the _ system operator _ @xmath102 has the following important properties @xmath103 where the latter property comes from the fact that @xmath104 and @xmath105 . also ,
recall that for any square matrix @xmath106 , one can write @xmath107 , where @xmath108 and @xmath109 denote the real and imaginary parts of the matrix @xmath106 , respectively .
now by the _ canonical system _ we mean the system whose state is described by a time - dependent @xmath110 taking values in the hilbert space @xmath111 with the standard inner product @xmath112 whose dynamics are governed by the odes ( [ radis8 ] ) .
these odes will be referred to as the _ canonical evolution equation_.    in this paper we will focus on the case when there is no forcing , i.e. , @xmath113 ( or , equivalently , @xmath69 ) . in this case , the canonical system represents a dissipative dynamical system ( since @xmath114 ) with the evolution governed by the semigroup @xmath115 . as we shall see in this paper and as mentioned in the introduction
, there are some serious advantages in studying this dissipative dynamical system , i.e. , the canonical system , over the lagrangian or hamiltonian systems .      in a previous work @xcite of ours , we studied canonical systems whose states are solutions of a linear evolution equation in the canonical form ( [ radis8 ] ) with system operator @xmath102@xmath116 @xmath117 , @xmath118 having exactly the properties ( [ radis10 ] ) in which the matrix @xmath99 did not have full rank .
in that paper , the canonical system was a simplified version of an abstract model of an oscillator damped retarded friction that modeled a two - component composite system with a lossy and lossless component .
we showed that the state @xmath110 of such a system satisfied the energy balance equation@xmath119   = -w_{\mathrm{dis}}\left [ v(t)\right ]   + w\left [   v\left (   t\right )   \right ]   , \label{radis11}\ ] ] in which the system energy @xmath120   $ ] , dissipated power @xmath121   $ ] , and @xmath122   $ ] the rate of work done by the force @xmath123 were given by @xmath124   = \frac{1}{2}\left (   v\left (   t\right ) , v\left (   t\right )   \right )   , \text { } w_{\mathrm{dis}}\left [   v(t)\right ] = \beta\left (   v\left (   t\right )   , bv\left (   t\right )   \right )   , \text { } w\left [   v\left (   t\right )   \right ]   = \operatorname{re}\left (   v\left ( t\right )   , f\left (   t\right )   \right )   .
\label{radis11a}\ ] ]    an important result of this paper which follows immediately from the block form ( [ radis8a ] ) , ( [ radis8b ] ) and the relation between the variables @xmath95 , @xmath125 , @xmath126 from ( [ radis7 ] ) , ( [ radis3 ] ) , ( [ dlag6 ] ) ] is that a state @xmath77 of the lagrangian system , whose dynamics are governed by second - order odes ( [ radis2a ] ) , and the corresponding state @xmath95 of the canonical system , whose dynamics are governed by the canonical evolution equation ( [ radis8 ] ) , the energetics are equivalent in the sense @xmath127   = \mathcal{t}\left (   \dot{q},q\right )   + \mathcal{v}\left ( \dot{q},q\right )   , \text { \ } w_{\mathrm{dis}}\left [   v\right ]   = 2\mathcal{r}\left (   \dot{q}\right )   , \text { \ } w\left [   v\right ]   = \operatorname{re}\left (   \dot{q},f\right )   \text { , } \label{radis11b}\ ] ] where , as you will recall from the energy balance equation ( [ radis4 ] ) , @xmath128 was the system energy for the lagrangian system .
* on the eigenmodes and the quality factor .
* we will be interested in the eigenmodes of the canonical system ( [ radis8 ] ) and their quality factor .
an eigenmode of the canonical system is defined as a nonzero solutions of the odes ( [ radis8 ] ) with no forcing , i.e. , @xmath113 , having the form @xmath129 .
the quality factor @xmath130 $ ] of such an eigenmode is defined as in @xcite by @xmath131=2\pi\frac{\text{energy stored in system}}{\text{energy lost per cycle}}=\left\vert \operatorname{re}\zeta\right\vert \frac{u\left [   v\left ( t\right )   \right ]   } { w_{\mathrm{dis}}\left [   v(t)\right ]   } .
\label{radis12}\ ] ] by the energy balance equation ( [ radis11 ] ) it follows that @xmath131=-\frac{1}{2}\frac{\left\vert \operatorname{re}\zeta\right\vert } { \operatorname{im}\zeta } , \label{radis12_1}\ ] ] with the convention @xmath130=+\infty$ ] if @xmath82 .
now given an eigenmode @xmath132 of the lagrangian system [ radis2a ] ) , it follows that the corresponding state @xmath95 of the canonical system ( [ radis8 ] ) , related by ( [ radis7 ] ) , is an eigenmode of the form @xmath133 ( excluding the case @xmath134 , which can only occur if @xmath135 ) .
importantly , the energetic equivalence ( [ radis11b ] ) holds and their quality factors ( [ radis4_3 ] ) , ( [ radis12_1 ] ) are equal , i.e.,@xmath136   = q_{\zeta}\text{. } \label{radis13}\ ] ]    a more detailed discussion on the eigenmodes of the two systems and the relationship between them can be found in section [ sspan ] .
one of the important applications of our methods described above is electric circuits and networks involving resistors representing losses .
a general study of electric networks with losses can be carried out with the help of the lagrangian approach , and that systematic study has already been carried out in this paper . for lagrangian treatment of electric networks and circuits we refer to ( * ? ? ?
9 ) , ( * ? ? ?
2.5 ) , @xcite .
we illustrate the idea and give a flavor of the efficiency of our methods by considering below a rather simple example of an electric circuit as in fig .
[ figc1 ] with the assumptions@xmath137 this example  has the essential features of two component systems incorporating high - loss and lossless components .    , @xmath138 , @xmath139 , two inductances @xmath140 , @xmath141 , a resistor @xmath142 , and two sources @xmath143 , @xmath144 .
this electric circuit example fits within the framework of our model .  indeed , since resistors represent losses , this two component system consists of a lossy component and a lossless component  the right and left circuits , respectively.,width=459,height=251 ]    * the lagrangian system . * to derive evolution equations for the electric circuit in fig .
[ figc1 ] we use a general method for constructing lagrangians for circuits , ( * ? ? ?
9 ) , that yields@xmath145 where @xmath48 and @xmath49 are respectively the kinetic and the potential energies , @xmath146 is the lagrangian , and @xmath31 is the rayleigh dissipative function .
notice that @xmath147 and @xmath148 are the currents .
the general euler - lagrange equations of motion with forces are , ( * ? ? ?
8),@xmath149 where @xmath77 are the charges and @xmath53 the sources@xmath150{c}q_{1}\\ q_{2}\end{array } \right ]   , \quad f=\left [ \begin{array } [ c]{c}e_{1}\\ e_{2}\end{array } \right ]   , \label{circ3}\ ] ] yielding from ( [ circ0])([circ3 ] ) the following second - order odes@xmath151 with the dimensionless loss parameter@xmath152 that scales the intensity of losses in the system , and@xmath153{cc}l_{1 } & 0\\ 0 & l_{2}\end{array } \right ]   , \text { \ \ } r=\left [ \begin{array } [ c]{cc}0 & 0\\ 0 & \ell \end{array } \right ]   , \text { \ \ } \eta=\left [ \begin{array } [ c]{cc}\frac{1}{c_{1}}+\frac{1}{c_{12 } } & -\frac{1}{c_{12}}\\ -\frac{1}{c_{12 } } & \frac{1}{c_{2}}+\frac{1}{c_{12}}\end{array } \right ]   .
\label{circ6}\ ] ] recall , the loss fraction @xmath1 defined in ( [ radis4a ] ) is the ratio of the rank of the matrix @xmath34 to the total degrees of freedom @xmath62 of the system which in this case is @xmath154 thus the lagrangian system ( [ circ4 ] ) has all the properties described in sections [ slagrangian ] , [ slagrmodeleigenmqfac ] and since the loss fraction condition ( [ radis4b ] ) , i.e. , @xmath0 , is satisfied then it is model of a two - component composite with a lossy and a lossless component .    * the canonical system . *
following the method in section [ scansys ] , we introduce the variables @xmath95 defined by ( [ dlag6 ] ) , ( [ radis3 ] ) , ( [ radis7 ] )  in terms of the matrix @xmath85 defined in ( [ radis6 ] ) , ( [ radis6a1 ] ) as@xmath155{ll}k_{\mathrm{p } } & 0\\ 0 & k_{\mathrm{q}}\end{array } \right ]   , \text { \ \ } k_{\mathrm{p}}=\sqrt{\alpha}^{-1}=\left [ \begin{array } [ c]{ll}\frac{1}{\sqrt{l_{1 } } } & 0\\ 0 & \frac{1}{\sqrt{l_{2}}}\end{array } \right ]   , \text { \ \ } k_{\mathrm{q}}=\sqrt{\eta}.\nonumber\end{gathered}\ ] ] as @xmath156 is a @xmath157 matrix in this case , then we can compute its square root explicitly using the formula@xmath158 where @xmath159 denotes the positive square root , which holds for any @xmath157 matrix @xmath160 with @xmath161 .
consequently , in the lagrangian system ( [ circ4 ] ) in terms of the variable @xmath95 is the desired canonical evolution equation from ( [ radis8 ] ) , namely , @xmath162 and@xmath163{cc}0 & -\mathrm{i}\phi^{\mathrm{t}}\\ \mathrm{i}\phi & 0 \end{array } \right ]   , \text { \ } b=\left [ \begin{array } [ c]{ll}\mathsf{\tilde{r } } & 0\\ 0 & 0 \end{array } \right ]   , \text { \ } f=\left [ \begin{array } [ c]{l}k_{\mathrm{p}}f\\ 0 \end{array } \right ]   , \label{circ11}\\ \phi = k_{\mathrm{q}}k_{\mathrm{p}}^{\mathrm{t}},\text { \ \ } \mathsf{\tilde{r}}=k_{\mathrm{p}}rk_{\mathrm{p}}^{\mathrm{t}}=\left [ \begin{array } [ c]{ll}0 & 0\\ 0 & \frac{\ell}{l_{2}}\end{array } \right ]   .\nonumber\end{gathered}\ ] ]
in this section we will introduce a new virial theorem for dissipative lagrangian systems in terms of the eigenmodes of the system .
this result generalizes the virial theorem from classical mechanics ( * ? ? ?
* ;  3.4 ) , as discussed in appendix [ apdxvirthm ] , for conservative lagrangian systems for these types of modes .
recall that the virial theorem is about the oscillatory transfer of energy from one form to another and its equipartition .
we can now state precisely and prove our generalization of the virial theorem to dissipative ( @xmath118 ) lagrangian systems  ( [ radis2a ] ) and the equipartition of energy for the oscillatory eigenmodes .
in fact , the proof in section [ sproofs ] shows that our theorem is true for more general system of odes in the form ( [ radis2a ] ) since essentially all that really matters is that the energy balance equation ( [ radis4 ] ) holds for the eigenmodes .
[ virial theorem][tvir]if @xmath70 is an eigenmode of the lagrangian system ( [ radis2a ] ) then the following identity holds if @xmath75 : @xmath164 where @xmath165 and @xmath166 are the kinetic and potential energy , respectively , defined in ( [ radis4_0 ] ) . on the other hand ,
if @xmath167  then the identity ( [ tvirs1 ] ) no longer holds and either @xmath135 or the identity @xmath168 must hold .
[ energy equipartition][tvir_eqp]for systems with @xmath169 , if @xmath70 is an eigenmode of the lagrangian system ( [ radis2a ] ) with @xmath75 , then the following identity holds@xmath170 in other words , for the oscillatory eigenmodes there is an equipartition of the system energy , i.e. , @xmath128 , between their kinetic energy @xmath48 and their potential energy @xmath49 .
in this section we study time - harmonic solutions to the euler - lagrange equation ( [ radis2a ] ) which constitutes a subject of the spectral theory .
an important subject of this section is the study of the relations between standard spectral theory and the quadratic pencil formulation as it arises naturally as the time fourier transformation of the euler - lagrange evolution equation .
as introduced in section [ slagrmodeleigenmqfac ] , the eigenmodes of the lagrangian system are nonzero solutions of the odes ( [ radis2a ] )  with @xmath69 having the form @xmath70 . for a fixed @xmath33 ,
these modes correspond to the solutions of the quadratic eigenvalue problem ( qep)@xmath171 for the quadratic matrix pencil in @xmath79@xmath172 the set of eigenvalues ( spectrum ) of the pencil @xmath173 is the set@xmath174 which are exactly those values @xmath79 for which a solution to the qep ( [ pfsp1 ] ) exists .
the spectral theory of polynomial operator pencils @xcite can be applied to study the eigenmodes but it has its disadvantages such as being more complicated than standard spectral theory .
thus an alternative approach to the spectral theory is desirable .
often the alternative approach is to use the hamiltonian system and consider its eigenmodes , that is , the nonzero solutions of the hamiltonian equations ( [ radis3 ] ) with @xmath69 having the form @xmath175 .
these modes correspond to the solutions of the eigenvalue problem@xmath176 for the matrix@xmath177{ll}\beta r & 0\\ 0 & 0 \end{array } \right ]   \right )   m_{\mathrm{h}}. \label{pfsp5}\ ] ] an advantage to this approach is the simple correspondence via ( [ dlag6 ] ) , ( [ radis3 ] ) between the set of modes of the two systems , namely , the eigenmodes of the hamiltonian system are solutions of ( [ radis3 ] ) having the block form @xmath178{l}p\\ q \end{array } \right ]   e^{-\mathrm{i}\zeta t}$ ] in which @xmath70 an eigenmode of the lagrangian system and @xmath179 . in particular , this means the matrix @xmath180 and the pencil @xmath181 have the same eigenvalues and hence the same spectrum , i.e.,@xmath182 a major disadvantage of this approach is that @xmath183 is a non - self - adjoint matrix such that the standard theory of self - adjoint or dissipative operators does not apply without further transformation of the system , even in the absence of losses , i.e. , @xmath184 . instead , in this case it is the krein spectral theory @xcite , @xcite , @xcite that is often used .
but this theory is far more complex than the standard spectral theory and much harder to apply .
our approach to these spectral problems which overcomes the disadvantages in the pencil or krein spectral theory is to use the canonical system and consider its eigenmodes , that is , the nonzero solutions of the canonical evolution equation ( [ radis8 ] ) with @xmath113 having the form @xmath185 .
these modes correspond to the solutions of the eigenvalue problem@xmath186 for the system operator @xmath187 with hermitian matrices @xmath98 , @xmath99 with @xmath188 .
the advantage of this approach is that the standard spectral theory can be used since @xmath189 is a dissipative operator and when losses are absent @xmath190 is self - adjoint .
this is a serious advantage since the spectral theory is significantly simpler and allows the usage of the deep and effective results from perturbation theory such as those developed in @xcite to study the modes of such canonical evolution equations having the general form ( [ radis8 ] ) .
this is the approach we take in this paper to study spectral symmetries and the dissipative properties of the eigenmodes of the lagrangian system such as the modal dichotomy and overdamping phenomenon .    *
correspondence between spectral problems .
* we now conclude this section by summarizing the correspondence between the two main spectral problems of this paper , namely , between the standard eigenvalue problem ( [ pfsp7 ] ) and the quadratic eigenvalue problem ( [ pfsp1 ] ) .
we do this in the next corollary which uses the following proposition that tells us the characteristic matrix of the system operator @xmath191 can be factored in terms of the quadratic matrix pencil @xmath192 .
first we introduce some notation that will be useful .
the hilbert space @xmath111 with standard inner product @xmath112 can be decomposed as @xmath193 into the orthogonal subspaces @xmath194 , @xmath195 with orthogonal matrix projections@xmath196{cc}\mathbf{1 } & 0\\ 0 & 0 \end{array } \right ]   , \text { \ \ } p_{\mathrm{q}}=\left [ \begin{array } [ c]{cc}0 & 0\\ 0 & \mathbf{1}\end{array } \right ]   . \label{pfsp8a}\ ] ] in particular , the matrices @xmath98 , @xmath99 , and @xmath102 defined in ( [ radis8])([radis8b ] ) are block matrices already partitioned with respect to the decomposition @xmath197 and any vector @xmath198 can be represented uniquely in the block form@xmath199{c}\varphi\\
\psi \end{array } \right ]   , \text { \ \ where } \varphi = p_{\mathrm{p}}w,\text { \ \ } \psi = p_{\mathrm{q}}w . \label{pfsp8b}\ ] ] then with respect to this decomposition we have the following results .    [ ppfsp]if @xmath200 then @xmath201{cc}k_{\mathrm{p } } & \zeta^{-1}\mathrm{i}\phi^{\mathrm{t}}\\ 0 & \mathbf{1}\end{array } \right ]   \left [ \begin{array } [ c]{cc}\zeta^{-1}\mathbf{1 } & 0\\ 0 & \zeta\mathbf{1}\end{array } \right ]   \left [ \begin{array } [ c]{cc}c(\zeta,\beta ) & 0\\ 0 & \mathbf{1}\end{array } \right ]   \left [ \begin{array } [ c]{cc}k_{\mathrm{p}}^{\mathrm{t } } & 0\\ -\zeta^{-1}\mathrm{i}\phi & \mathbf{1}\end{array } \right ]   .
\label{pfsp8}\ ] ]    [ spectral equivalence][cpfsp]for any @xmath202 , @xmath203 in particular , the system operator @xmath102 and quadratic matrix pencil @xmath204 have the same spectrum , i.e. , @xmath205 moreover , if @xmath200 then the following statements are true :    1 .   if @xmath206 and @xmath207 then @xmath199{c}-\mathrm{i}\zeta\sqrt{\alpha}q\\ \sqrt{\eta}q \end{array } \right ]   , \text { \ \ where } c(\zeta,\beta)q=0,\text { \ \ } q\not = 0 . \label{cpfsp3}\ ] ] 2 .   if @xmath208 and @xmath209 then @xmath210{c}-\mathrm{i}\zeta\sqrt{\alpha}q\\ \sqrt{\eta}q \end{array } \right ]   \not = 0 . \label{cpfsp4}\ ] ]      in section [ stvspen ] it was shown that the study of the eigenmodes of the lagrangian system ( [ radis2a ] ) reduces to the quadratic eigenvalue problem ( [ pfsp1 ] ) and this motivates a study of the spectrum @xmath211 of the quadratic matrix pencil @xmath204 in ( [ pfsp2 ] ) .
corollary [ cpfsp ] makes it clear that we can instead study the eigenmodes of the canonical system ( [ radis8 ] ) and , in particular , the system operator @xmath102 spectrum satisfies @xmath212 .
the purpose of this section is to give a detailed analysis of the set @xmath213 .
recall , the system operator @xmath214 , @xmath118 from ( [ radis8 ] ) with @xmath24 matrices @xmath98 , @xmath99 has the fundamental properties ( [ radis10a ] ) , ( [ radis10b])@xmath215 where @xmath216 .
these properties are particularly important in describing the spectrum @xmath213 of the system operator @xmath102 .
for instance , the next proposition on the spectrum for nondissipative ( @xmath184 ) lagrangian systems ( [ radis2a ] ) follows immediately from these properties which would otherwise not be exactly obvious for gyroscopic systems ( i.e. , @xmath217 ) .
[ real eigenfrequencies]for nondissipative lagrangian systems ( [ radis2a ] ) , that is when @xmath218 , all the eigenfrequencies @xmath79 are real , and consequently any eigenmode evolution is of the form @xmath219 with @xmath82 .    if @xmath184 ( i.e. , no dissipation ) then @xmath220 is a hermitian matrix
thus , the spectrum @xmath221 is a subset of @xmath222 .
hence , if @xmath70 is an eigenmode of the lagrangian systems ( [ radis2a ] ) with @xmath184 then by corollary [ cpfsp ] we have @xmath223 and so @xmath82 .
this completes the proof .    in the next few sections we will give a deeper analysis of the spectral properties of @xmath102 including spectral symmetries in section [ sspsy ] and the modal dichotomy in sections [ sevmd ] , [ smdhlr ] . in order to do
so we must first introduce some notation . in the hilbert space @xmath111 ,
denote by @xmath224 , @xmath225 the nonzero eigenvalues of @xmath99 ( counting multiplicities ) with the smallest denoted by @xmath226 in particular , the spectrum of @xmath99 is @xmath227 where @xmath228 .
denote the largest eigenvalue of @xmath98 by @xmath229 .
it follows from the fact that @xmath98 is a hermitian matrix which is skew - symmetric that @xmath230 where @xmath231 denotes the operator norm on square matrices .
the next proposition describes the spectral symmetries of the system operator @xmath232 which follow from the property @xmath233 .
[ spectral symmetry][pssym]the following statements are true :    1 .
the characteristic polynomial of @xmath102 satisfies @xmath234 for every @xmath202 .
in particular , the spectrum @xmath235 of the system operator @xmath232 has the symmetry@xmath236 2 .   if @xmath237 is an eigenvector of the system operator @xmath232 with corresponding eigenvalue @xmath79 then @xmath238 is an eigenvector of @xmath232 with corresponding eigenvalue @xmath239 .
3 .   if @xmath184 ( i.e. , no dissipation ) then @xmath240 @xmath241 for every @xmath202 .
we will denote the discs centered at the eigenvalues of @xmath242 with radius @xmath229 by@xmath243 two  subsets of the spectrum @xmath244 of the system operator @xmath245 @xmath117 which play a central role in our analysis are@xmath246    [ eigenvalue bounds][pevbd]the following statements are true :    1 .
the eigenvalues of the system operator @xmath102 lie in the union of the closed discs whose centers are the eigenvalues of @xmath242 with radius @xmath229 , that is,@xmath247 2 .   if @xmath207 and @xmath206 then@xmath248 3 .
if @xmath79 is an eigenvalue of @xmath102 and @xmath249 then @xmath250    [ spectral clustering][cssym]the eigenvalues of the system operator @xmath187 , @xmath118 lie in the closed lower half of the complex plane , are symmetric with respect to the imaginary axis , and lie in the union of the closed discs whose centers are the eigenvalues of @xmath242 with radius @xmath229 .  moreover , if @xmath184 ( i.e. , no dissipation ) then the eigenvalues of @xmath220 are real and symmetric with respect to the origin .
[ modal dichotomy i][tmdic]if @xmath251 then @xmath252 furthermore , there exists unique  invariant subspaces @xmath253 , @xmath254 of the system operator @xmath255 @xmath117 with the properties@xmath256 where @xmath111 .
moreover , the dimensions of these subspaces satisfy@xmath257    [ high - loss susceptible subspace][dmdic]for the system operator @xmath258 @xmath117 with @xmath259 we will call its @xmath260-dimensional invariant subspace @xmath254 the high - loss susceptible subspace .  we will call its @xmath261-dimensional invariant subspace @xmath262 the low - loss susceptible subspace .    our reasoning for the definitions of these subspaces is clarified with the following corollary .
[ high - loss subspace : dissipative properties][cmdic]if @xmath263 then@xmath264 furthermore , the quality factor ( [ radis12_1 ] ) of any eigenmode of the canonical system ( [ radis8 ] ) in the high - loss susceptible subspace @xmath254 satisfies @xmath265\leq\frac{1}{2}\frac{\omega_{\text{max}}}{\beta b_{\text{min}}-\omega_{\text{max}}}<\frac{1}{2}. \label{cmdic2}\ ] ] in particular , as the losses go to @xmath266 the damping factor and quality factor of any such eigenmode goes to @xmath267 and @xmath268 , respectively , that is,@xmath269=0 . \label{cmdic3}\ ] ]    we conclude this section with the following remarks .
the spectrum of @xmath270 restricted to the low - loss susceptible subspace @xmath262 , that is , the set @xmath271 in ( [ cmdic1 ] ) , is close to the real axis and actually coalesces to a finite set of real numbers as losses @xmath272 .
this statement is made precise in the following section with theorem [ cpodr ] and the asymptotic expansions in ( [ pod14 ] ) .
on the other hand , results on quality factor for the eigenmodes in @xmath273 is far more subtle than the results in corollary [ cmdic ] for the eigenmodes in the high - loss susceptible subspace @xmath254 . for gyroscopic - dissipative systems considered in this paper , corollary [ cmdic ] above and
proposition [ pqaollm ] below give a partial description of the nature of the quality factor for large losses , i.e. , @xmath68 , in the lossy component of the composite system .
when gyroscopy is absent , i.e. , @xmath169 , then a more complete analysis for quality factor can be carried out , which we have done in section [ sodan ] of this paper in connection to our studies on the overdamping phenomenon .
consideration of quality factor and overdamping in dissipative systems with gyroscopy , however , requires a more subtle and detailed analysis that will be considered in a future work .
we are interested in describing the spectrum @xmath274 of the system operator @xmath275 , @xmath118 in the high - loss regime , i.e. , @xmath68 .
we do this in this section by giving an asymptotic characterization , as @xmath272 , of the modal dichotomy as described in theorem [ tmdic ] and corollary [ cmdic ] .  in order to do
so we need to give a spectral perturbation analysis of the matrix @xmath270 as @xmath272 .
fortunately , this analysis has already been carried out in @xcite .
we now introduce the necessary notion and describe the results .
the hilbert space @xmath111 with standard inner product @xmath112 is decomposed into the direct sum of orthogonal invariant subspaces of the operator @xmath99 , namely @xmath276 where @xmath277 ( the range of @xmath99 ) is the loss subspace of dimension @xmath278 with orthogonal projection @xmath279 and its orthogonal complement , @xmath280 ( the nullspace of @xmath99 ) , is the no - loss subspace of dimension @xmath281 with orthogonal projection @xmath282 .
the operators @xmath98 and @xmath99 with respect to the direct sum ( [ pod5 ] ) are the @xmath157 block operator matrices@xmath283{cc}\omega_{2 } & \theta\\ \theta^{\ast } & \omega_{1}\end{array } \right ]   , \quad b=\left [ \begin{array } [ c]{cc}b_{2 } & 0\\ 0 & 0 \end{array } \right ]   , \label{pod7}\ ] ] where @xmath284 and @xmath285 are restrictions of the operators @xmath98 and @xmath99 respectively to loss subspace @xmath286 whereas @xmath287 is the restriction of @xmath98 to complementary subspace @xmath288 .
also , @xmath289 is the operator @xmath290 whose adjoint is given by @xmath291 .    the perturbation analysis in the high - loss regime @xmath68 for the system operator @xmath292 described in ( * ?
* a , theorem 5 & proposition 11 ) introduces an orthonormal basis @xmath293 diagonalizing the self - adjoint operators @xmath294 and @xmath295 from ( [ pod7 ] ) with@xmath296 where@xmath297 then for @xmath68 the system operator @xmath298 is diagonalizable with basis of eigenvectors @xmath299 satisfying@xmath300 which split into two distinct classes @xmath301 with the following properties .    * the high - loss class : * the eigenvalues have poles at @xmath302 whereas their eigenvectors are analytic at @xmath302 , having the asymptotic expansions@xmath303 the vectors @xmath304 , @xmath305 form an orthonormal basis of the loss subspace @xmath286 and@xmath306 in particular , @xmath224 , @xmath225 are all the nonzero eigenvalues of @xmath99 ( counting multiplicities ) .    * the low - loss class : * the eigenvalues and eigenvectors are analytic at @xmath302 ,
having the asymptotic expansions@xmath307 the vectors @xmath304 , @xmath308 form an orthonormal basis of the no - loss subspace @xmath309 and@xmath310    by ( * ? ? ?
* a , proposition 7 ) we know that the asymptotic formulas for the real and imaginary parts of the complex eigenvalues @xmath311 as @xmath272 are given by @xmath312    observe that the expansions ( [ pod16 ] ) imply @xmath313 justifying the names high - loss and low - loss .
the following theorem is the goal of this section .
it characterizes the spectrum @xmath244 of the system operator @xmath187 and the modal dichotomy from theorem [ tmdic ] and corollary [ cmdic ] in the high - loss regime @xmath68 in terms of the high - loss and low - loss eigenvectors .
[ modal dichotomy ii][cpodr]for @xmath33 sufficiently large , the modal dichotomy occurs as in theorem [ tmdic ] and corollary [ cmdic ] with the following equalities holding : @xmath314 and@xmath315 in particular , @xmath254 and @xmath253 , the high - loss and low - loss susceptible subspaces of @xmath270 , respectively , have as a basis the high - loss eigenvectors @xmath316 and the low - loss eigenvectors @xmath317 , respectively .    in the next section
we will study overdamping phenomena for lagrangian systems ( [ radis2a ] ) .
an important role in our analysis will be played by the eigenmodes @xmath318 , @xmath319 of the canonical system ( [ radis8 ] ) which by the modal dichotomy split into the two distinct classes based on their dissipative properties:@xmath320    it should be emphasized here that the classification of these modes into high - loss and low - loss is based solely on the behavior of their damping factor in ( [ pod17 ] ) as losses become large , i.e. , as @xmath272 , and not necessarily on their quality factor which we will discuss at the end of this section .
moreover , the damping factor is related to energy loss by the energy balance equation ( [ radis11 ] ) , ( [ radis11a ] ) satisfied by these modes which implies the dissipated power is@xmath321   &   = -2\operatorname{im}\zeta_{j}\left (   \beta\right )   u[v_{j}\left (   t,\beta\right )   ] , \text { \ \ } 1\leq j\leq2n,\label{pod20a}\\ u[v_{j}\left (   t,\beta\right )   ]   &   = \frac{1}{2}\left (   w_{j}\left ( \beta\right )   , w_{j}\left (   \beta\right )   \right )   e^{2\operatorname{im}\zeta_{j}\left (   \beta\right )   t}\nonumber\end{aligned}\ ] ] where @xmath322 $ ] defined in ( [ radis11a ] ) was interpreted as the system energy . in particular , by ( [ pod17 ] ) , ( [ pod20a ] ) , and the fact @xmath323 it follows that their system energy for any _ fixed _
satisfies@xmath325=0,\text { \ \ } 1\leq j\leq n_{r};\label{pod20b}\\ \text{low - loss }   &   \text { : \ \ } \lim_{\beta\rightarrow\infty}u[v_{j}\left ( t,\beta\right )   ] = \frac{1}{2},\text { \ \ } n_{r}+1\leq j\leq2n.\nonumber\end{aligned}\ ] ] this combined with ( [ pod17 ] ) , ( [ pod20a ] ) is the justification for the names high - loss and low - loss modes .    * asymptotic formulas for the quality factor . *  the perturbation analysis in the high - loss regime @xmath68 for the quality factor @xmath326   $ ] , @xmath319 of the eigenmodes from ( [ pod20 ] ) as given by ( [ radis12 ] ) , ( [ radis12_1 ] ) has already been carried out in ( * ? ? ?
* a , prop .
we will now describe the results .
the quality factor @xmath327   $ ] , @xmath328 for each high - loss eigenmode has a series expansion containing only odd powers of @xmath329 implying for @xmath68 it is a nonnegative decreasing function in the loss parameter @xmath33 which has the asymptotic formula as @xmath272@xmath330   = \frac { 1}{2}\frac{\left\vert \rho_{j}\right\vert } { b_{j}}\beta^{-1}+o\left ( \beta^{-3}\right )   , \ \
1\leq j\leq n_{r}. \label{pod21}\ ] ] it follows from this discussion that as @xmath272 @xmath330 \searrow0,\ \
1\leq j\leq n_{r } \label{pod21_1}\ ] ] ( here we will use the notation @xmath331 @xmath268 or @xmath332 to denote that a function of @xmath33 is decreasing to @xmath268 or increasing to @xmath267 , respectively , as @xmath272 ) .    the quality factor @xmath327   $ ] , @xmath308 for each low - loss eigenmode has two possibilities : ( i ) @xmath327   = + \infty$ ] for @xmath68 ; ( ii ) @xmath327   $ ] has a series expansion containing only odd powers of @xmath329 implying either @xmath326   \searrow0 $ ] or @xmath333   \nearrow+\infty$ ] as @xmath272 . finding necessary and sufficient conditions for when the cases occur is a subtle problem which is still open in general , but is solved completely in the next section on overdamping when @xmath169 under the nondegeneracy condition @xmath334
. we will now consider the cases @xmath335 [ cf .
( [ pod22 ] ) , ( [ pod22_1 ] ) ] or @xmath336 ( cf . prop .
[ pqaollm ] ) .  in the former case ( which is the typical case see ( * ? ? ? * a , remark 9 ) ) , we have as @xmath272 the asymptotic formula@xmath337   = \frac{1}{2}\frac{\left\vert \rho_{j}\right\vert } { d_{j}}\beta+o\left (   \beta ^{-1}\right )   , \ \
\text{if } d_{j}\not = 0 \label{pod22}\ ] ] and so it follows in this case that as @xmath272 @xmath337 \searrow0\text { , if } \rho_{j}=0\text { and } q\left [   w_{j}\left (   \beta\right ) \right ]   \nearrow+\infty,\text { if } \rho_{j}\not = 0\text{. } \label{pod22_1}\ ] ]    * asymptotic oscillatory low - loss modes . *
now the low - loss modes from ( [ pod20 ] ) with @xmath336 will play a key role in the study of the selective overdamping phenomenon in section [ spad ] .
we call such modes the asymptotic oscillatory low - loss modes since if @xmath336 then the limiting function @xmath338 ( for fix @xmath339 ) is an oscillatory function in @xmath339 with period @xmath340 .
thus as @xmath272 we denote this by @xmath341 from the quality factor discussion above for the low - loss modes  the next proposition follows immediately , regardless of whether @xmath335 or not .
[ q - factor : asymptotic oscillatory low - loss modes][pqaollm ] any low - loss eigenmodes from ( [ pod20 ] ) with @xmath336 has a quality factor @xmath327   $ ] with the property that either @xmath327   = + \infty$ ] for all @xmath68 or @xmath327   \nearrow+\infty$ ] as @xmath272 ( i.e. , is an unbounded increasing function of the loss parameter @xmath33 ) .
the phenomenon of _ overdamping _ ( also called _ heavy damping _ ) is best known for a simple damped oscillator .
namely , when the damping exceeds certain critical value all oscillations cease entirely , see , for instance , ( * ? ?
2 ) . in other words , if the damped oscillations are described by the exponential function @xmath342 with a complex constant @xmath79 then in the case of overdamping ( heavy damping ) @xmath343 .
our interest to overdamping is motivated by the fact that if an eigenmode becomes overdamped it can not resonate at any finite frequency .
consequently , the contribution of such a mode to losses at finite frequencies becomes minimal , and that provides a mechanism for the absorption suppression for systems composed of lossy and lossless components .
the treatment of overdamping for systems with many degrees of freedom involves a number of subtleties particularly in our case when the both lossy and lossless degrees of freedom are present , i.e. ,  when the loss fraction condition ( [ radis4b ] ) is satisfied .
we will show that any lagrangian system ( [ radis2a ] ) with @xmath169 can be overdamped with selective overdamping occurring whenever @xmath34 does not have full rank , i.e. , @xmath344 , which corresponds to a model of a two - component composite with a high - loss and a low - loss component .
let us be clear that we study overdamping in this paper only for the case @xmath169 as the case @xmath217 will differ significantly in the analysis and it is not expected that there exists such a critical value for damping when all oscillations cease entirely for any damped oscillation .
a study of a proper generalization of overdamping in the case @xmath217 will be carried out in a future publication .
we make here statements and provide arguments for overdamping for lagrangian systems .
we will study the phenomenon known as overdamping in this section under the following condition :    [ no gyrotropy][cndod]for our study of overdamping we assume henceforth that@xmath345    the following proposition is a key result that will be referenced often in our study of overdamping .
it describes the spectrum of the @xmath24 matrices @xmath99 , @xmath98 in terms of the spectrum of the @xmath10 matrices @xmath346 , @xmath347 .
[ spectra relations][pspre]for the @xmath24 matrices @xmath99 , @xmath98 and the @xmath348 matrices @xmath89 , @xmath90 , @xmath34 we have @xmath349 ( @xmath350 ) and @xmath351 for every @xmath202 .  in particular ,
if @xmath352 and @xmath229 denote the smallest nonzero eigenvalue and the largest eigenvalue of @xmath99 and @xmath98 , respectively , and @xmath353 denotes the smallest positive eigenvalue of @xmath98 then @xmath354 and@xmath355      to study whether overdamping is possible we need to determine conditions under which@xmath356 has a solution for the system operator @xmath357 .
but we know from corollary [ cpfsp ] that @xmath358 if and only if @xmath359 hence the solutions to ( [ lod1 ] ) are the solution of@xmath360 thus to study overdamping it suffices to study conditions under which the quadratic eigenvalue problem@xmath361 has a solution .  to be explicit we have@xmath362 as it turns out the key to studying overdamping will be to study zero sets of the family of quadratic forms@xmath363 on certain subspaces of the hilbert space @xmath194 .
in particular , we find that@xmath364   ^{2}-\frac{\left (   q,\eta q\right ) } { \left (   q,\alpha q\right )   } } .\nonumber\end{aligned}\ ] ] thus we find that a necessary and sufficient condition for an eigenvalue @xmath79 of the system operator @xmath102 to satisfy @xmath343
is@xmath365    now there are some fundamental inequalities that play a key role in our overdamping analysis which are described in the following proposition
. we will prove part of this proposition , the rest is proved in section [ sproofs ] , since the proof is enlightening showing how the energy conservation law ( [ radis4 ] ) and the virial theorem [ tvir ] can be used to derive these inequalities .
[ fundamental inequalities][pfuineq]the following statements are true .    1 .   for any @xmath366 with @xmath209,@xmath367 moreover , if @xmath368 exists then@xmath369 2 .
if the matrix @xmath34 has full rank , i.e. , @xmath370 , then for any @xmath366 with @xmath209 , @xmath371 3 .
for any @xmath372 @xmath373with @xmath75,@xmath374 4 .
if @xmath375 [ or equivalently , @xmath376 then @xmath343 whenever any one of the following two inequalities is satisfied : ( i ) @xmath377 or ( ii ) @xmath368 exists and @xmath378 . 5 .
if @xmath70 is an eigenmode of the lagrangian system ( [ radis2a ] ) then energy equipartition , i.e. , @xmath379 , can only hold if @xmath380 or , in the case @xmath368 exists , if @xmath381 .
the first two statements are proved in section [ sproofs ] .
the third statement will now be proved using the energy conservation law ( [ radis4 ] ) and the virial theorem [ tvir ] , or more specifically , corollary [ tvir_eqp ] on the equipartition of energy .
@xmath373with @xmath382 .
then @xmath70 is an eigenmode of the lagrangian system ( [ radis2a ] ) .
it follows from ( [ radis4 ] ) and corollary [ tvir_eqp ] that for the system energy @xmath383 we have@xmath384 where @xmath385 and @xmath386 are the kinetic and potential energy , respectively , of this state @xmath77 and @xmath387 is the dissipated power .
notice that for @xmath55 we have @xmath388 .
it now follows from this and ( [ lod7_0c ] ) , ( [ lod7_0d ] ) that@xmath389 the proof of the third statement now follows from these two equalities and the inequality ( [ lod7_0 ] ) .
the fourth statement follows immediately from the inequality ( [ lod7_0b ] ) and the inequality ( [ lod7_00 ] ) .
the fifth statement follows from the fact that if @xmath390 then the equality ( [ lod7_0f ] ) holds , regardless of whether @xmath75 or not as long as @xmath200 , and so the result follows immediately from the inequalities ( [ lod7_0 ] ) and ( [ lod7_00 ] ) if @xmath200 . and
since the existence of @xmath368 prohibits @xmath135 , the proof now follows .
this completes the proof .    from the inequalities ( [ lod7_0a ] ) , ( [ lod7_0b ] ) one can see clearly
that if the matrix @xmath34 has full rank then all the eigenvalues of @xmath270 are purely imaginary once @xmath391 in other words , ( [ lod7_3 ] ) is a sufficient condition for all of the eigenmodes of the canonical system ( [ radis8 ] ) [ or , equivalently , the lagrangian system ( [ radis2a ] ) ] to be overdamped , provided @xmath34 has full rank .
but it is also clear that this argument must be refined if @xmath34 does not have full rank , i.e. , @xmath344 , since then @xmath392 and hence @xmath393 .
the refinement we use is perturbation theory . in particular
, we use our results in sections [ sspsy ] , [ smdhlr ] , [ svir ] on the modal dichotomy and the virial theorem to derive our main results on overdamping .
our argument is essentially based on determining which of the eigenmodes can not maintain the energy equipartition in corollary [ tvir_eqp ] , i.e. , the equality between the kinetic and potential energy , when @xmath33 is sufficiently large by either using the inequalities in corollary [ cmdic ] or using the asymptotic expansions in section [ smdhlr ] for the eigenmodes in ( [ pod20 ] ) .
this is where the fifth statement in proposition [ pfuineq ] is relevant .
we derive the following results on overdamping using the inequalities in proposition [ pfuineq ] .
[ complete overdamping][tcovd]if @xmath394 and @xmath395 has full rank , i.e. , @xmath370 , then all the eigenmodes of the canonical system ( [ radis8 ] ) are overdamped , i.e. , for the system operator @xmath187 we have@xmath396    [ overdamped regime is optimal][eodop]the following example shows for the class of system operators @xmath397 under consideration in this paper , the regime @xmath251 is optimal for guaranteeing overdamping in the sense that we can always find an example of a system operator satisfying @xmath398 with @xmath33 as close to @xmath399 as we like such that the set @xmath400 is empty .
in particular , it is proven using the following family of system operators@xmath163{cc}0 & -\mathrm{i}\\ \mathrm{i } & 0 \end{array } \right ]   , \text { \ \ } b=\left [ \begin{array } [ c]{cc}1 & 0\\ 0 & 0 \end{array } \right ]   , \label{eodop1}\\ a\left (   \beta\right )   = \omega-\mathrm{i}\beta b,\text { \ \ } \beta \geq0,\nonumber\end{gathered}\ ] ] derived from the system with one degree of freedom having lagrangian and rayleigh dissipation function @xmath401{l}\dot{q}\\ q \end{array } \right ]   ^{\mathrm{t}}m_{\mathrm{l}}\left [ \begin{array } [ c]{l}\dot{q}\\ q \end{array } \right ]   , \quad m_{\mathrm{l}}=\left [ \begin{array } [ c]{ll}1 & 0\\ 0 & -1 \end{array } \right ]   , \label{eodop2}\\ r=\frac{1}{2}\left (   \dot{q},\beta r\dot{q}\right )   , \text { \ \ } r=1\text{,}\nonumber\end{gathered}\ ] ] which is just is the lagrangian for a one degree - of - freedom spring - mass - damper system with mass @xmath402 , spring constant @xmath402 , and viscous damping coefficient @xmath33 .
then we have@xmath403 and therefore@xmath404 in this case , at the boundary of the overdamping regime , where @xmath405 , we have critical damping of the system .
now we want to answer the question of which modes of the canonical system ( [ radis8 ] ) are overdamped in the case when @xmath34 does not have full rank , i.e. , @xmath344 .
the next theorem gives us a partial answer to the question .
[ partial overdamping][tpovd]if @xmath251 and @xmath344 then the modal dichotomy occurs as in theorem [ tmdic ] and corollary [ cmdic ] .
in addition , any eigenmode of the canonical system ( [ radis8 ] ) in the invariant subspace @xmath406 ( the high - loss susceptible subspace ) is overdamped , i.e. , @xmath407    now we ask the important and natural question : can any of the modes of the system operator @xmath102 in the low - loss susceptible subspace @xmath273 be overdamped when @xmath34 from the rayleigh dissipative function @xmath31 in ( [ radis1 ] ) does not have full rank , i.e. , in the case when the lagrangian system is a model of a two - component composite with a lossy and a lossless component ?
this question will be investigated in the next section .
we are interested in solving the problem of describing all the modes of the canonical system ( [ radis8 ] ) [ or , equivalently , the lagrangian system ( [ radis2a ] ) ] , with system operator @xmath245 @xmath117 , that are overdamped when losses are sufficiently large , i.e. , @xmath68 .  in the previous section we solved the problem for lagrangian systems with complete rayleigh dissipative , i.e.
, @xmath34 has full rank , for then theorem [ tcovd ] tells us that all the modes will become overdamped . on the other hand , for lagrangian systems with incomplete rayleigh dissipative , i.e. , @xmath34 does not have full rank ,
then theorem [ tpovd ] is a partial solution to the problem since it tells us that , due to sufficiently large losses , modal dichotomy occurs splitting the modes into the two @xmath102-invariant subspaces @xmath273 and @xmath406 ( the high - loss and low - loss susceptible subspaces , respectively ) with the modes belonging to @xmath254 overdamped .
thus to solve the problem completely we need to describe which of the modes in the low - loss susceptible subspace @xmath253 can be overdamped when losses are sufficiently large . in this section
we solve the problem asymptotically as losses become extremely large , i.e. , as @xmath272 , using the virial theorem [ tvir ] and , in particular , energy equipartition [ tvir_eqp ] for the eigenmodes and the perturbation theory developed in @xcite which we described in section [ smdhlr ] .
our description of the selective overdamping phenomenon is based on the modal dichotomy in the high - loss regime @xmath68 as described in section [ smdhlr ] .
in particular , we use theorem [ cpodr ] to characterize overdamping in terms of the high - loss and low - loss eigenmodes @xmath408 , @xmath319 of the canonical system ( [ radis8 ] ) which split into the two distinct classes based on their dissipative properties ( [ pod20 ] ) .
we say that such an eigenmode @xmath409 is _ overdamped _ for @xmath68 if @xmath410 for all
@xmath33 sufficiently large .
we say that such an eigenmode _ remains oscillatory _ for @xmath68 if @xmath411 for all
@xmath33 sufficiently large .
our first result tells us that all of the high - loss eigenmodes @xmath408 , @xmath305 from ( [ pod20 ] ) are overdamped regardless of whether @xmath34 has full rank or not .
it is an immediate corollary of theorem [ tpovd ] and theorem [ cpodr ] .
[ overdamped high - loss modes][thmod]all of the high - loss eigenmodes from ( [ pod20 ] ) are overdamped for @xmath68 .
now we will discuss and give a solution to the problem of determining which of the low - loss eigenmodes @xmath412 , @xmath308 in ( [ pod20 ] ) will be overdamped .  in order to do this
we work under the following nondegeneracy condition :    [ nondegeneracy][cndnd]for our study of the overdamping of the low - loss modes in this section , we assume henceforth that @xmath413    we have two main goals in this section . our first is to find necessary and sufficient conditions to have @xmath414 for @xmath68 .
let @xmath415 be the number of low - loss eigenmodes that are overdamped for @xmath68 .
our second goal is to calculate this number @xmath415 .
we begin by recalling some of the important properties of the low - loss eigenpairs @xmath416 , @xmath417 , @xmath308 in ( [ pod11 ] ) for the system operator @xmath418 that we will need in this section .  by ( [ pod14 ] ) and ( [ pod16 ] ) the limits@xmath419 exist and from ( [ pod8 ] ) , ( [ pod15 ] ) the limiting vectors @xmath420 are an orthonormal basis for @xmath421 which also diagonalize the self - adjoint operator @xmath294 from ( [ pod7 ] ) satisfying@xmath422    * overview of main results . *
before we proceed , let us now summarize the results in this section .
first , its obvious from this discussion that a necessary condition for the low - loss eigenmode @xmath409 to be overdamped for @xmath68 is @xmath423 ( even though it does so with a slow damping of order @xmath329 , as seen in ( [ llod1 ] ) and ( [ pod14 ] ) ) , but what is not obvious is that this is in fact a sufficient condition .
although in the case @xmath368 exists ( i.e. , @xmath424 ) , it just follows from the fourth statement in proposition [ pfuineq ] since if @xmath423 then by ( [ llod1 ] ) we have @xmath425 which implies @xmath426 for @xmath427 and hence @xmath428 for @xmath68 . but this argument fails when  @xmath429 .
thus we must find an alternative method to show that @xmath423 is equivalent to @xmath410 for @xmath68 .
the key idea that leads to an alternative method is to realize that by the virial theorem [ tvir ] and the fifth statement in proposition [ pfuineq ] , when @xmath368 exists the modes with @xmath423 are exactly the modes that break the energy equipartition .
hence our general method is to show , even in the case @xmath424 , that for a low - loss modes with @xmath423 , there is a breakdown of the virial theorem [ tvir ] since as losses increase , i.e. , as @xmath272 , the equality between the kinetic and potential energy ( see corollary [ tvir_eqp ] ) of the corresponding eigenmode of the lagrangian system ( [ radis2a ] ) can no longer be maintained and must eventually break which forces the modes to be non - oscillatory for @xmath68 , i.e. , overdamped .
this proves that @xmath423 is equivalent to @xmath428 for @xmath68 .
thus it follows that the low - loss eigenmodes in ( [ pod20 ] ) of the canonical system ( [ radis8 ] ) which are overdamped for @xmath68 , their limiting vectors ( satisfy @xmath430 @xmath431 for @xmath339 fixed ) form a orthonormal basis for @xmath431 .
this implies of course that @xmath432 is the number of low - loss eigenmodes that are overdamped for @xmath68 . a major result of this section
is that in fact @xmath433 , where @xmath434 .
but @xmath435 is also the number of high - loss eigenpairs in ( [ pod11 ] ) .
thus we conclude using these results and theorem [ thmod ] that the number of eigenmodes in ( [ pod20 ] ) which are overdamped in the high - loss regime @xmath68 is exactly @xmath436 with @xmath260 of these being low - loss eigenmodes and the rest being all the high - loss eigenmodes .
moreover , in the case @xmath344 , the remaining @xmath437 modes are low - loss oscillatory modes with an extremely high quality factor that actually increases as the losses increase , i.e. , as @xmath272 .
* main results . *
we now give the statements of our main results .
the proof of the first theorem we give in this section as it is enlightening .
the proofs of the rest of the statements below are found in section [ sproofs ] .
[ overdamped low - loss modes][tlmod]a necessary and sufficient condition for a low - loss eigenmode @xmath438 from ( [ pod20 ] ) to be overdamped for @xmath68 is @xmath423 ( or equivalently , @xmath439 ) .
let @xmath440 be a low - loss eigenmode from ( [ pod20 ] ) .
if it is overdamped for @xmath68 then @xmath410 for @xmath68 ( i.e. , for all @xmath33 sufficiently large ) and so by ( [ llod1 ] ) , ( [ llod2 ] ) we have @xmath423 and @xmath441 ( the nullspace of @xmath294 ) .
lets now prove the converse .
suppose @xmath423 .
then we must show that @xmath428 for @xmath68 .
suppose that this was not true .
then since @xmath416 was analytic at @xmath302 we must have @xmath442 for @xmath68 .
thus by corollary [ cpfsp ] there is a corresponding eigenmode @xmath443 of the lagrangian system ( [ radis2a ] ) such that with respect to the block representation ( [ pfsp8b ] ) the vector @xmath417 is @xmath444{c}-\mathrm{i}\zeta_{j}\left (   \beta\right )   \sqrt{\alpha}q_{j}\left ( \beta\right ) \\
\sqrt{\eta}q_{j}\left (   \beta\right ) \end{array } \right ]   .
\label{llod3}\ ] ] in particular , @xmath445 and @xmath446 are related by ( [ radis7 ] ) , namely,@xmath447 by section [ seneq ] on the energetic equivalence between the canonical system ( [ radis8 ] ) and the lagrangian system ( [ radis2a ] ) it follows that@xmath448   = \mathcal{t}\left (   \dot{q},q\right )   + \mathcal{v}\left (   \dot{q},q\right )   .
\label{llod5}\ ] ] the term on the left is the system energy of the eigenmode @xmath445 of the canonical system defined in ( [ radis11a ] ) whereas the sum on the right is the system energy of the eigenmode @xmath446 of the lagrangian system ( [ radis2a ] ) , i.e. , the sum of its kinetic and potential energy .
as proved in section [ pspad ] , for any fixed time @xmath339 we have @xmath449   = \lim_{\beta\rightarrow\infty}\mathcal{v}\left (   \dot{q}_{j}\left ( t\,,\beta\right )   , q_{j}\left (   t\,,\beta\right )   \right )   , \label{llod6}\\ \lim_{\beta\rightarrow\infty}\mathcal{t}\left (   \dot{q}_{j}\left ( t\,,\beta\right )   , \dot{q}_{j}\left (   t\,,\beta\right )   \right )   = 0.\nonumber\end{gathered}\ ] ] this proves that @xmath450 and , consequently , by the energy equipartition ( see corollary [ tvir_eqp ] ) that@xmath451 this completes the proof .
[ clmod]the collection of limiting values as @xmath272 ( with @xmath339 fixed ) of the low - loss eigenmodes in ( [ pod20 ] ) of the canonical system ( [ radis8 ] ) which are overdamped for @xmath68 is @xmath452 @xmath453 and is an orthonormal basis for @xmath431 .
[ plmod]let @xmath415 denote the number low - loss eigenmodes in ( [ pod20 ] ) of the canonical system ( [ radis8 ] ) which are overdamped for @xmath68 .
then @xmath454 where @xmath455 .
[ number of overdamped modes][tnlm]the total number of eigenmodes of the canonical system ( [ radis8 ] ) from ( [ pod20 ] ) which are overdamped for @xmath68 is @xmath436 , where @xmath455 .
furthermore , @xmath260 of these are low - loss eigenmodes and the rest are all of the high - loss eigenmodes . moreover ,
if @xmath344 then the @xmath456 modes which are not overdamped as @xmath68 are low - loss eigenmodes which remain oscillatory for @xmath68 .
[ selective overdamping][tseod]if @xmath344 then there is exactly @xmath456 low - loss eigenmodes of the canonical system ( [ radis8 ] ) from ( [ pod20 ] ) which remain oscillatory for @xmath68 ( i.e. , @xmath411 for @xmath68 ) and their quality factors as given by ( [ radis12 ] ) , ( [ radis12_1 ] ) have the following property @xmath457 = + \infty\text { for } \beta\gg1\text { or } q\left [   w_{j}\left (   \beta\right ) \right ]   \nearrow+\infty\text { as } \beta\rightarrow\infty .
\label{tseod1}\ ] ]
this section contains the proofs of the results from sections [ svir ] , [ sspan ] and [ sodan ] .
proofs are organized in the order in which statements appeared and grouped according to the sections and subsections in which they are located in the paper .
all assumptions , notation , and conventions used here will adhere to that previously introduced in this paper .
the proofs of statements in section [ svir ] are contained here .
[ proof of theorem [ tvir]]for @xmath70 an eigenmode of the lagrangian system ( [ radis2a ] ) , we consider the virial@xmath458 in our proof we will use the real quantities @xmath459 and @xmath460 . if @xmath135 we are done .
thus , suppose that @xmath200 .
the virial ( [ vir1 ] ) satisfies@xmath461 taking the time - derivative of ( [ vir9 ] ) and using ( [ radis2a ] ) yields the identity@xmath462 where @xmath463 .
the energy balance equation ( [ radis4 ] ) implies@xmath464 combining the identities ( [ vir10 ] ) and ( [ vir11 ] ) yields immediately the identity@xmath465 from ( [ vir12 ] ) and since @xmath466 we derive the identity@xmath467 if @xmath75 or @xmath468 if @xmath343 .
suppose now that @xmath75 .
then , by taking the imaginary and real parts of the identity ( [ vir14 ] ) , we find that@xmath469 the rest of the proof of theorem [ tvir ] now immediately follows from ( [ vir16 ] ) and ( [ radis4_0 ] ) .
[ proof of corollary [ tvir_eqp]]corollary [ tvir_eqp ] follows immediately from the virial theorem [ tvir ] .
the proofs of statements in section [ sspan ] are contained here .
[ proof of proposition [ ppfsp]]the factorization ( [ pfsp8 ] ) follows from formula ( [ mab2 ] ) in appendix [ apxsc ] and ( [ radis6a1 ] ) , ( [ radis6b ] ) , ( [ radis8])([radis8b ] ) , ( [ pfsp2 ] )  since@xmath470{cc}\zeta\mathbf{1}-\omega_{\mathrm{p}}+\mathrm{i}\mathsf{\tilde{r } } & \mathrm{i}\phi^{\mathrm{t}}\\ -\mathrm{i}\phi & \zeta\mathbf{1}\end{array } \right ] \\ = \left [ \begin{array } [ c]{cc}\mathbf{1 } & \zeta^{-1}\mathrm{i}\phi^{\mathrm{t}}\\ 0 & \mathbf{1}\end{array } \right ]   \left [ \begin{array } [ c]{cc}\zeta\mathbf{1}-\omega_{\mathrm{p}}+\mathrm{i}\mathsf{\tilde{r}}-\zeta ^{-1}\mathrm{i}\phi^{\mathrm{t}}\left (   -\mathrm{i}\phi\right )   & 0\\ 0 & \zeta\mathbf{1}\end{array } \right ]   \left [ \begin{array } [ c]{cc}\mathbf{1 } & 0\\ -\zeta^{-1}\mathrm{i}\phi & \mathbf{1}\end{array } \right ] \\ = \left [ \begin{array } [ c]{cc}k_{\mathrm{p } } & \zeta^{-1}\mathrm{i}\phi^{\mathrm{t}}\\ 0 & \mathbf{1}\end{array } \right ]   \left [ \begin{array } [ c]{cc}\zeta^{-1}\mathbf{1 } & 0\\ 0 & \zeta\mathbf{1}\end{array } \right ]   \left [ \begin{array } [ c]{cc}c(\zeta,\beta ) & 0\\ 0 & \mathbf{1}\end{array } \right ]   \left [ \begin{array } [ c]{cc}k_{\mathrm{p}}^{\mathrm{t } } & 0\\ -\zeta^{-1}\mathrm{i}\phi & \mathbf{1}\end{array } \right ]   .\end{gathered}\ ] ] this completes the proof .    [ proof of corollary [ cpfsp]]the proof of this corollary follows immediately from the factorization ( [ pfsp8 ] ) using the facts @xmath471 , @xmath472 and@xmath473{cc}k_{\mathrm{p}}^{\mathrm{t } } & 0\\ -\zeta^{-1}\mathrm{i}\phi & \mathbf{1}\end{array } \right ]   ^{-1}=\left [ \begin{array } [ c]{cc}\sqrt{\alpha } & 0\\ \zeta^{-1}\mathrm{i}\phi\sqrt{\alpha } & \mathbf{1}\end{array } \right ]   .\ ] ]      [ proof of proposition [ pssym]]the proof of this proposition follows immediately from elementary properties of the determinant and using the fundamental property of the system operator @xmath275 in ( [ sasp1 ] ) , namely , @xmath474 which implies @xmath475 .      [ proof of proposition [ pevbd]]the second statement follows immediately from ( [ sasp1 ] ) since if @xmath206 and @xmath207 then@xmath476 from this and using ( [ sasp1])([sasp3 ] ) , the first and third statements follow immediately from proposition [ apxpebd ] in appendix [ apxebd ] .
this completes the proof .
[ proof of corollary [ cssym]]the proof of this corollary follows immediately from propositions [ pssym ] and [ pevbd ] .
[ proof of theorem [ tmdic]]suppose that @xmath251 .
it follows from theorem [ apxtmd ] in appendix [ apxebd ] that the @xmath24 matrix @xmath102 satisfies the hypotheses of the theorem and together with ( [ sasp1 ] ) , ( [ ebmd1 ] ) implies @xmath477 and the sets @xmath478 , @xmath479 in ( [ mdss ] ) satisfy@xmath480 moreover , theorem [ apxtmd ] gives the existence of unique projection matrices @xmath481 , @xmath482 such that in the space @xmath111 the subspaces @xmath483 and @xmath484 have exactly the desired properties in theorem [ tmdic ] with their dimensions satisfying@xmath485 this completes the proof .    [ proof of corollary [ cmdic]]first , if @xmath237 is an eigenvector of @xmath270 with corresponding eigenvalue @xmath79 then proposition [ pevbd ] and ( [ sasp3 ] ) imply@xmath486 assume now that @xmath251 .
then by ( [ sasp1])([sasp3 ] ) , ( [ pcmdic1 ] ) , and theorem [ tmdic ] the proof of the first part of corollary [ cmdic ] follows immediately from corollary [ apxcmd ] in appendix [ apxebd ] .
the rest of the proof now follows immediately from ( [ pcmdic1 ] ) , formula ( [ radis12_1 ] ) , and the first part of corollary [ cmdic ] .
this completes the proof .
[ proof of theorem [ cpodr]]consider the high - loss and low - loss eigenpairs @xmath416 , @xmath417 , @xmath487 for @xmath68 introduced in section [ smdhlr ] .
the have the properties that @xmath488 are a basis of eigenvectors for the system operator @xmath102 satisfying @xmath489 the proof of theorem [ cpodr ] follows from these facts , theorem [ tmdic ] , and corollary [ cmdic ] .
the proofs of statements in section [ sodan ] are contained here .
recall , that we are assuming from now on that condition [ cndod ] is true i.e. , @xmath490    [ proof of proposition [ pspre]]using the block representation of the matrices @xmath98 , @xmath99 in ( [ radis8a ] ) and the result ( [ mab5 ] ) in appendix [ apxsc ] on the determinant of block matrices with commuting block entries we find that for any @xmath202 we have@xmath491{cc}\zeta\mathbf{1 } & \mathrm{i}\phi^{\mathrm{t}}\\ -\mathrm{i}\phi & \zeta\mathbf{1}\end{array } \right ]   \right )   = \det\left (   \zeta^{2}\mathbf{1-}\phi^{\mathrm{t}}\phi\right )   = \det\left (   \zeta^{2}\mathbf{1}-\alpha^{-1}\eta\right )   , \\
\det\left (   \zeta\mathbf{1}-b\right )    &   = \det\left (   \left [ \begin{array } [ c]{cc}\zeta\mathbf{1-}\mathsf{\tilde{r } } & 0\\ 0 & \zeta\mathbf{1}\end{array } \right ]   \right )   = \zeta^{n}\det\left (   \zeta\mathbf{1-}\mathsf{\tilde{r}}\right )   = \zeta^{n}\det\left (   \zeta\mathbf{1-}\alpha^{-1}r_{\mathrm{h}}\right )   , \end{aligned}\ ] ] since by ( [ radis6a1 ] ) , ( [ radis6b ] ) , ( [ radis8b ] ) we have @xmath492 and @xmath493 .
the proof now follows from these facts .
[ proof of proposition [ pfuineq]]the last three statements have already been proved .
the first two statements are proved using proposition [ pspre ] and the following series of inequalities which just use elementary properties of the operator norm @xmath231 for hermitian matrices @xmath494 where @xmath495 or @xmath34 and the inequality ( [ lod7_1a ] )  requires the hypothesis that @xmath90 or @xmath34 , respectively , be invertible so that @xmath496 .
the completes the proof .
[ proof of theorem [ tcovd]]suppose @xmath394and @xmath497 , where @xmath216 . if @xmath79 was an eigenvalue of @xmath102 with @xmath75 then by proposition [ pfuineq ] we have @xmath498 a contradiction .
therefore , all the eigenvalues of the system operator @xmath102 are real .
this completes the proof .
[ proof of theorem [ tpovd]]suppose @xmath394and @xmath499 .
then by theorem [ tmdic ] and corollary [ cmdic ] , if @xmath500 then @xmath501 and so by proposition [ pfuineq ] we must have @xmath343 .
therefore , all the eigenvalues of @xmath270 in @xmath502 are real .
[ proof of theorem [ thmod]]if @xmath438 is a high - loss eigenmode from ( [ pod20 ] ) then by the modal dichotomy ( [ cpodr ] ) we must have latexmath:[$\zeta_{j}\left (   \beta\right )   \in\sigma\left (   a\left (   \beta\right )    @xmath68 ( i.e. , for @xmath33 sufficiently large ) . by theorem [ tpovd ] it follows that @xmath428 for @xmath68 .
therefore , @xmath409 is overdamped for @xmath68 .
this completes the proof .
we will now assume in the rest of the proofs that condition [ cndnd ] is true i.e. , @xmath504 .
[ proof of theorem [ tlmod]]the proof of this theorem was started in section [ spad2 ] .
it remains to show that ( [ llod6 ] ) holds .
first , we have that @xmath505 with @xmath506 and @xmath507 by the properties of the low - loss eigenmodes in ( [ pod20 ] ) and since ( [ llod1 ] ) , ( [ llod2 ] ) is true for these modes together with the hypothesis @xmath423 .
now we use the representation ( [ llod3 ] ) and ( [ plmodp6 ] ) for @xmath417 and @xmath431 , respectively , to imply that as @xmath272 @xmath444{c}-\mathrm{i}\zeta_{j}\left (   \beta\right )   \sqrt{\alpha}q_{j}\left ( \beta\right ) \\
\sqrt{\eta}q_{j}\left (   \beta\right ) \end{array } \right ]   \rightarrow\mathring{w}_{j}=\begin{bmatrix } 0\\ \mathring{\psi}\end{bmatrix } , \text { } -\mathrm{i}\zeta_{j}\left (   \beta\right )   \sqrt{\alpha}q_{j}\left ( \beta\right )   \rightarrow0,\text { } \sqrt{\eta}q_{j}\left (   \beta\right ) \rightarrow\mathring{\psi } \label{tlmod1}\ ] ] for some @xmath508 .
it follows from ( [ tlmod1 ] ) and the fact @xmath509 that for any fixed time @xmath339 we have as @xmath272@xmath510    &   = \frac{1}{2}\left ( w_{j}\left (   \beta\right )   , w_{j}\left (   \beta\right )   \right ) e^{2\operatorname{im}\zeta_{j}\left (   \beta\right )   t}\rightarrow\frac{1}{2}\left (   \mathring{w}_{j},\mathring{w}_{j}\right )   = \frac{1}{2}.\end{aligned}\ ] ] this proves ( [ llod6 ] ) since by ( [ tlmod1 ] ) we have @xmath511 .
this completes the proof .
[ proof of corollary [ clmod]]this corollary follows immediately from theorem [ tlmod ] and the fact that the @xmath512 are an orthonormal basis for @xmath421 which diagonalize the self - adjoint operator @xmath294 such that ( [ llod2 ] ) is satisfied .
[ proof of theorem [ tnlm]]this theorem follows immediately from theorems [ thmod ] and [ tlmod ] , corollary [ clmod ] , and proposition [ pqaollm ] ( which is proved below ) .
[ proof of theorem [ tseod]]this theorem follows immediately from theorems [ tnlm ] , [ tlmod ] and proposition [ pqaollm ] .
[ proof of proposition [ plmod]]let @xmath415 denote the number of low - loss eigenmodes of the canonical system ( [ radis8 ] ) in ( [ pod20 ] ) which are overdamped for @xmath68 .  then by corollary [ clmod ] we know that@xmath513 we will now prove that @xmath514 where by ( [ radis10a ] ) , @xmath455 .
we will use the elementary facts from linear algebra that since @xmath90 , @xmath515 are @xmath10 matrices which are positive semidefinite then @xmath516 another elementary fact from linear algebra which we need is that if @xmath517 and @xmath518 are subspaces in the same hilbert space then@xmath519 where @xmath520 .
thus with @xmath521 , @xmath522 we have@xmath523 where the sum is direct by condition [ cndnd ] , i.e. , the hypothesis @xmath334 . in particular
, this implies @xmath524    the proof of ( [ plmodp0 ] ) will be carried out in a series of steps which we outline now before we begin .
first , we will prove that @xmath525{c}0\\ \psi \end{array } \right ]   : \sqrt{\eta}\psi\in\operatorname{ran}\eta\cap\operatorname{ran}r\right\ }   , \label{plmodp6}\ ] ] where the block form is with respect to the decomposition @xmath193 , @xmath111 , @xmath526 described in section [ stvspen ] , i.e. , the form ( [ pfsp8b ] ) .
next , it follows immediately from ( [ plmodp1 ] ) that @xmath527 finally , we prove that @xmath528 the proof of ( [ plmodp0 ] ) follows immediately from ( [ plmodp6])([plmodp8 ] ) .
we begin by computing the representation ( [ plmodp6 ] ) for @xmath431 .
first , by its definition from ( [ pod7 ] ) , @xmath294 is the restriction of the operator @xmath98 to the no - loss subspace @xmath280 , that is,@xmath529 where @xmath282 is the orthogonal projection onto @xmath288 in the hilbert space @xmath111 with standard inner product @xmath112 .
denote by @xmath530 the orthogonal projection onto @xmath531 in the hilbert space @xmath194 with the standard inner product .
then it follows from the block representation of @xmath98 and @xmath99 in ( [ pod7 ] ) with respect to the decomposition @xmath193 that @xmath282 and @xmath532 with respect to this decomposition are the block operators @xmath533 from ( [ plmodp9 ] ) , ( [ plmodp10 ] ) it follows that@xmath525{c}\varphi\\ \psi \end{array } \right ]   \in\operatorname{ker}\mathsf{\tilde{r}}\oplus h_{\mathrm{q}}:-\mathrm{i}p_{\mathsf{\tilde{r}}}^{\bot}\phi^{\mathrm{t}}\psi=0\text { and } \mathrm{i}\phi p_{\mathsf{\tilde{r}}}^{\bot}\varphi=0\right\ }   .\ ] ] but if @xmath534 with @xmath535 then @xmath536 by definition of @xmath530 and hence @xmath537 and @xmath538 which implies @xmath539 and  @xmath540 since by ( [ radis8b ] ) and ( [ radis6a1 ] ) we have @xmath541 , @xmath542 , @xmath543 , and @xmath544 . by the hypothesis of condition [ cndnd ] this implies that @xmath545 . also , since @xmath546 is a hermitian matrix it follows that @xmath547 is the orthogonal projection onto the @xmath548 so that @xmath549 is equivalent to @xmath550 and this is equivalent to @xmath551 @xmath67 by properties ( [ radis6b ] ) . but @xmath551 @xmath67 is equivalent to @xmath551 @xmath552 which by ( [ plmodp1 ] ) is equivalent to @xmath551 @xmath553 .
this proves ( [ plmodp6 ] ) .
finally , we will now prove ( [ plmodp8 ] ) .
it follows from the fact that @xmath90 is a hermitian matrix and ( [ plmodp1 ] ) that the restriction of the operator@xmath88 to @xmath554 , i.e. , @xmath555 , is an invertible operator and that @xmath556 from this it follows that@xmath557 the proof of ( [ plmodp8 ] ) immediately follows from this fact and ( [ plmodp5 ] ) .
this completes the proof of the theorem .
let @xmath558{cc}p & q\\ r & s \end{array } \right ]   \label{mab1}\ ] ] be a square matrix represented in block form where @xmath559 and @xmath560 are square matrices with the latter invertible , that is , @xmath561 . then the following _ aitken block - diagonalization formula
_ holds @xcite , @xcite,@xmath558{cc}p & q\\ r & s \end{array } \right ]   = \left [ \begin{array } [ c]{cc}\mathbf{1 } & qs^{-1}\\ 0 & \mathbf{1}\end{array } \right ]   \left [ \begin{array } [ c]{cc}p - qs^{-1}r & 0\\ 0 & s \end{array } \right ]   \left [ \begin{array } [ c]{cc}\mathbf{1 } & 0\\ s^{-1}r & \mathbf{1}\end{array } \right ]   , \label{mab2}\ ] ] where the matrix@xmath562 is known as the _ schur complement _ of @xmath560 in @xmath106 . the aitken formula ( [ mab2 ] ) readily implies@xmath563 from this we can conclude that@xmath564 but the latter statement is true even if @xmath565 which is proved using a limiting argument .  indeed , by substituting for @xmath560 in ( [ mab1 ] ) the small perturbation @xmath566 which is invertible for @xmath567 and commutes with any matrix that @xmath560 commutes with then statement [ mab5 ] is true for @xmath568 and the limit as @xmath569 yields the desired result .
in this appendix we discuss results pertaining to bounds on the eigenvalues of non - hermitian matrices in terms of their real and imaginary parts .
we also discuss a phenomenon , which we call _ modal dichotomy _ , that occurs for any non - hermitian matrix with a non - invertible imaginary part which is large in comparison to the real part of the matrix .
the bounds and dichotomy described here are extremely useful tools in studying the spectral properties of dissipative systems especially for composite systems with high - loss and lossless components .
recall , any square matrix @xmath106 can be written as the sum of a hermitian matrix and a skew - hermitian matrix,@xmath570 the real and imaginary parts of @xmath106 , respectively , and @xmath571 is the conjugate transpose of the matrix @xmath106 .    for the set of @xmath572 matrices , denote by @xmath231 the operator norm@xmath573 where @xmath574 and @xmath112 denotes the standard euclidean inner product on @xmath575 column vectors with entries in @xmath576 .
denote the spectrum of a square matrix @xmath106 by @xmath577 .
recall in the case @xmath106 is a hermitian matrix , i.e. , @xmath578 , or a normal matrix , i.e. , @xmath579 , that @xmath580 where @xmath581 and by convention @xmath582 if @xmath583 @xmath577 .    the first result which we will prove here
is the following proposition which gives bounds on the eigenvalues of a matrix @xmath106 in terms of @xmath584 and @xmath585 .
[ eigenvalue bounds][apxpebd]let @xmath106 be a square matrix .
then@xmath586 in other words , the eigenvalues of @xmath106 lie in the union of the closed discs whose centers are the eigenvalues of @xmath587 with radius @xmath588 .
moreover , the imaginary part of any eigenvalue @xmath79 of @xmath106 satisfies the inequality@xmath589 for some @xmath590 depending on @xmath79 .    before we proceed with the proof we will need the following perturbation result from ( * ? ? ?
v.4 , p. 291
4.8 ) for non - hermitian matrices which we prove here for completeness .    if @xmath591 is a square matrix with @xmath592 normal then@xmath593    let @xmath591 be a square matrix with @xmath592 normal .
then@xmath594 let @xmath595 , i.e. , the resolvent set of @xmath592 .
then for any @xmath596 we have@xmath597 thus if @xmath598 then the matrix @xmath599 is invertible since the neumann series @xmath600 converges to @xmath601 for any square matrix @xmath602 satisfying @xmath603 .
this implies@xmath604 but the inequality@xmath605 and the fact @xmath592 is a normal matrix implies@xmath606 therefore,@xmath607 this completes the proof .
[ proof of proposition [ apxpebd]]let @xmath106 be any square matrix .
then @xmath608 and @xmath609 is a normal matrix . by the previous lemma
this implies@xmath610 now let @xmath79 be any eigenvalue of @xmath106 .
then for any corresponding eigenvector @xmath611 of unit norm@xmath612 moreover , we know there exists an eigenvalue @xmath613 of @xmath585 ( depending on  @xmath79 ) such that @xmath614 hence @xmath615 the proof now follows from ( [ papxpebd1 ] ) and ( [ papxpebd2 ] ) .
the next result we prove is a proposition on dichotomy of the spectrum of a non - hermitian matrix @xmath106 which we mentioned in the introduction of this appendix .
in particular , the following proposition tells us that the spectrum of a matrix @xmath106 will split into two disjoint parts when the imaginary part @xmath585 is non - invertible and large `` in comparison to the real part @xmath584 .
the term large '' means that the bottom of the nonzero spectrum of @xmath616 must be greater than twice the top of the spectrum of @xmath617 ( where for a square matrix @xmath232 , @xmath618 as defined in ( * ? ? ?
* , p. 196 ) ) , that is , @xmath619 this of course implies @xmath620 is large " in comparison to @xmath588 since it follows from this inequality that@xmath621    [ modal dichotomy][apxtmd]let @xmath106 be any @xmath572 matrix which is non - hermitian such that its imaginary part @xmath585 is non - invertible with rank @xmath622 .
denote by @xmath623 , @xmath624
, @xmath625 the nonzero eigenvalues of @xmath585 .
if @xmath626 then @xmath627 where@xmath628 furthermore , there exists unique matrices @xmath629 , @xmath630 with the properties@xmath631 where @xmath632 if @xmath633 and @xmath634 if @xmath635 ( @xmath636 ) . moreover , these projection matrices have rank satisfying@xmath637    all the eigenvalues of @xmath585 by hypothesis are @xmath623 , @xmath624  , @xmath638 and @xmath639 .
it follows from proposition [ apxpebd ] that@xmath640 and the imaginary part of any eigenvalue @xmath79 of @xmath106 satisfies the inequality@xmath641 for some @xmath642 depending on @xmath79 .
in particular , it follows that@xmath643 where @xmath644 now suppose that @xmath645 .
we will now show that @xmath646 .
suppose this were not true then we could find @xmath647 which would imply for some @xmath648 with @xmath649 we have @xmath650 a contradiction .
thus @xmath651 .
now it follows from the spectral theory of matrices that the finite - dimensional vector space @xmath652 of @xmath575 column vectors with entries in @xmath576 can be written as the direct sum of two invariant subspaces for the matrix @xmath106,@xmath653 that is , @xmath654 is the union of the generalized eigenspaces of @xmath106 corresponding to eigenvalues in @xmath655 and @xmath656 is the union of generalized eigenspaces of @xmath106 corresponding to the eigenvalues in @xmath657 .
denote the projection matrix onto @xmath654 along @xmath656 by @xmath629 , that is , the matrix satisfying @xmath658 , @xmath659 , and @xmath660 .
then it follows from the spectral theory of matrices that @xmath629 and @xmath661 are the unique matrices with the properties@xmath662 now it follows that @xmath663 so that to complete the proof of this proposition we need only show that@xmath664    in order to prove this we begin by giving an explicit representation of the matrix @xmath629 in terms of the resolvent of @xmath106 which is defined as @xmath665 to do this we define @xmath666 we will now prove that@xmath667 first , we notice that @xmath668 which proves @xmath669 .
next , we have that if @xmath670 then for any @xmath623 with @xmath671 , @xmath672 implying that @xmath673  @xmath674 .
this proves @xmath675 .
it now follows from the cauchy - riesz functional calculus @xcite that@xmath676 where the contour integral is over the simply closed positively oriented path @xmath677 , @xmath678 .
now we consider the family of matrices@xmath679 notice that for @xmath680 we have@xmath681 it follows from this that the results proven so far in this theorem for the matrix @xmath106 apply to the matrix @xmath682 for each @xmath683   $ ] .  in particular
, for each @xmath684   $ ] we have @xmath685 where@xmath686 furthermore , there exists unique matrices @xmath687 , @xmath688 with the properties@xmath689 moreover , @xmath690 it follows from these facts that @xmath687 is a continuous matrix projection - valued function of @xmath339 in @xmath691   $ ] and@xmath692 where @xmath693 denotes the trace of a square matrix .
but @xmath694 is a continuous function of @xmath339 since @xmath687 is a continuous matrix - valued function implying @xmath695 is a continuous function of @xmath339 taking values in the nonnegative integers only .  from this
we conclude @xmath695 is constant for @xmath696 @xmath691   $ ] .  in particular ,
@xmath697 moreover,@xmath698 but this implies that@xmath699 this completes the proof .
we conclude this appendix by giving an alternative characterization of the subsets @xmath655 and @xmath674 from theorem [ apxtmd ] in terms of the magnitude of the imaginary parts of the eigenvalues of @xmath106 .
in particular , the next corollary tells us that if the @xmath585 is associated with losses and if the imaginary part of the eigenvalues of @xmath106 are associated with damping then only the eigenmodes with eigenvalues in @xmath674 are susceptible to large damping , relative to the norm of @xmath584 , when losses are large .
[ apxcmd]if @xmath645 then@xmath700    by theorem [ apxtmd ] we know that@xmath701 where@xmath702 and @xmath623 , @xmath624
, @xmath638 are the nonzero eigenvalues of @xmath585 .  in particular
, these facts imply immediately that@xmath703 now if @xmath704 then there exists a nonzero eigenvalue @xmath613 of @xmath585 such that@xmath705 which by the hypothesis @xmath645 implies@xmath706 the corollary follows immediately now from these facts .
the term of _ energetics _ refers to a fundamental set up for a system evolution when two energies are defined for any system configuration , namely the _ kinetic energy _ and the _
potential energy_. here is a concise description of the energetics given by h. poincar , @xcite :    _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ the difficulties inherent in the classic mechanics have led certain minds to prefer a new system they call _
energetics_.    energetics took its rise as an outcome of the discovery of the principle of the conservation of energy .
helmholtz gave it its final form .
it begins by defining two quantities which play the fundamental role in this theory .
they are _ kinetic energy , or vis viva _ , and _
potential energy_.    all the changes which bodies in nature can undergo are regulated by two experimental laws :    1@xmath707 .
the sum of kinetic energy and potential energy is constant .
this is the principle of the conservation of energy .
if a system of bodies is at @xmath232 at the time @xmath708 and at @xmath99 at the time @xmath709 , it always goes from the first situation to the second in such a way that the mean value of the difference between the two sorts of energy , in the interval of time which separates the two epochs @xmath708 and @xmath709 , may be as small as possible .
this is hamilton s principle , which is one of the forms of the principle of least action . " _
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _    in other words according to energetics the kinetic and potential energies @xmath48 and @xmath49 that are related to the lagrangian @xmath3 and the total energy @xmath50 by the following relations@xmath710 implying@xmath711 we assume the total energy to be equal to the hamiltonian , that is to be defined by the lagrangian through the legendre transformation@xmath712 then we arrive at the following uniquely defined expressions for the kinetic and potential energies@xmath713 we are particularly interested in the case of the quadratic lagrangian@xmath714 where @xmath89 and @xmath90 are real symmetric matrices which are positive definite and positive semidefinite , respectively , and @xmath9 is a real skew - symmetric matrix , that is@xmath715 then according to ( [ lagham3 ] ) the total energy @xmath50 the takes the form@xmath716 consequently , as it follows from ( [ lagham2 ] ) the kinetic and the potential energies are defined by@xmath717 the euler - lagrange equations corresponding to the quadratic lagrangian ( [ lagham5 ] ) are @xmath718 multiplying the above equation by @xmath719 from the left yields@xmath720 and this equation can be recast recast as@xmath721 combining the above equation with ( [ lagham2 ] ) and ( [ lagham5 ] ) we arrive at the following important identity that holds for any solution @xmath77 to the euler - lagrange equation ( [ lagham9])@xmath722 where the term @xmath723 is often referred to as the _ virial_.
let us introduce now the time average @xmath724 for a time dependent quantity @xmath725 defined by@xmath726 and observe that if @xmath123 is bounded for @xmath727 then @xmath728 notice now that in the case of interest when @xmath729 and consequently @xmath730 , all eigenfrequencies of the system are real and hence all solutions to the euler - lagrange equation ( [ lagham10 ] ) are bounded . consequently , applying the relation ( [ ftaver2 ] ) , using the identity ( [ lagham12 ] ) we obtain the virial theorem:@xmath731    * acknowledgment and disclaimer : * the research of a. figotin was supported through dr .
a. nachman of the u.s .
air force office of scientific research ( afosr ) , under grant number fa9550 - 11 - 1 - 0163 .
both authors are indebted to the referee for his valuable comments on our original paper .
v. i. arnold , _ mathematical methods of classical mechanics _ , translated from the russian by k. vogtmann and a. weinstein , second edition , graduate texts in mathematics , 60 , springer - verlag , new york , 1989 .