 persistent for the modes with $m=0,1,2$ and $l=1,2$. This observation raises the important question: What are the electromagnetic QNMs for which one has to look in astrophysical data. Also interesting is that while the $\epsilon$-method leads to significant changes of the frequencies $\omega_{m,n}$, it affects much less the separation parameter $E_{m,n}$ which here for the first time was obtained directly as a solution of the two-dimensional system without any prior approximations for it.

Another general result is that the confluent Heun function proved to be an effective tool for physical problems. Even though its <span style="font-variant:small-caps;">maple</span> realization still has many flaws, its precision proved to be good enough to repeat the already published results, and also studying the solutions, we were able to reveal new properties of the numerical stability of the EM QNMs with respect to changes in the phase-condition.

An interesting question is the results obtained using this method for $a>M$ or the so called naked singularity regime. Preliminary results show that the method is applicable in this case as well and the results will be published elsewhere.

Acknowledgements
================

The authors would like to thank Emanuele Berti for discussion of the numerical values of the EM QNM frequencies obtained through Leaver’s method which was important for the comparison of our method with this already well-established one.

P.F. is deeply indebted to Irina Aref’ieva, Alexei Starobinsky, Bruno Coppi and Michail Sazhin for interesting comments and suggestions about the relations of EM QNM to BH physics and the possibilities to discover such modes in astrophysical observations.

The authors are also thankful to Luciano Rezzolla for drawing our attention to the references [@headon0] and [@GW_] and to Shahar Hod for drawing our attention to the references [@extr1] and [@extr2].

This article was supported by the Foundation “Theoretical and Computational Physics and Astrophysics”, by the Bulgarian National Scientific Fund under contracts DO-1-872, DO-1-895, DO-02-136, and Sofia University Scientific Fund, contract 185/26.04.2010.

Author Contributions
====================

P.F. posed the problem of the evaluation of the EM QNMs of
persistent for the modes with $ m=0,1,2 $ and $ l=1,2$. This observation raises the authoritative doubt: What are the electromagnetic QNMs for which one has to look in astrophysical data. besides interesting is that while the $ \epsilon$-method leads to significant changes of the frequency $ \omega_{m, n}$, it affects much less the interval parameter $ E_{m, n}$ which here for the first prison term was obtained directly as a solution of the two - dimensional organization without any anterior approximations for it. 

 Another general resultant role is that the confluent Heun function proved to be an effective instrument for physical problems. Even though its < span style="font - variant: little - caps;">maple</span > realization still has many flaws, its precision prove to be good enough to repeat the already published results, and also study the solutions, we were able to reveal new properties of the numerical constancy of the EM QNMs with esteem to changes in the phase - condition. 

 An interesting question is the results obtained using this method for $ a > M$ or the so called bare singularity government. Preliminary results usher that the method is applicable in this case as well and the results will be publish elsewhere. 

 Acknowledgements 
 = = = = = = = = = = = = = = = = 

 The authors would like to thank Emanuele Berti for discussion of the numerical value of the EM QNM frequencies obtained through Leaver ’s method which was important for the comparison of our method with this already well - install one. 

 P.F. is deeply indebted to Irina Aref’ieva, Alexei Starobinsky, Bruno Coppi and Michail Sazhin for interesting comments and suggestions about the relations of EM QNM to bohrium physics and the possibility to detect such modes in astrophysical observations. 

 The authors are also thankful to Luciano Rezzolla for drawing our attention to the reference point [ @headon0 ] and [ @GW _ ] and to Shahar Hod for reap our attention to the references [ @extr1 ] and [ @extr2 ]. 

 This article was hold by the Foundation “ Theoretical and Computational Physics and Astrophysics ”, by the Bulgarian National Scientific Fund under contracts DO-1 - 872, DO-1 - 895, DO-02 - 136, and Sofia University Scientific Fund, contract 185/26.04.2010. 

 Author Contributions 
 = = = = = = = = = = = = = = = = = = = = 

 P.F. posed the problem of the evaluation of the EM QNMs of
 pegsistent for the modes wlth $m=0,1,2$ and $l=1,2$. This observetion rzises thd important question: What arx thw elextromagnetic QNMs for dhich one has to oook un astrophbaical data. Alal incecesting is that while the $\epsilon$-method ldabs to significant changes of the frezuenciex $\lmega_{m,n}$, it affgcts ktch mvsw the separation parameter $E_{m,h}$ which here for the first time was obtained dlrechly as a solution lf the two-dumenfuonal system without any prior apptoximations for it.

Another generau resblt is that tye fmnfluent Henn funbtion proved bp be at effecyive tool for ihysiral problems. Even though mts <span style="font-vatiant:small-wa's;">maple</span> realizatiin stiln hav mavt fuawa, mts precidioi proved to be good eniugh to repeat the slwvsdy publishes resujtf, and also studying the solutions, we wege agle to reveal new propeeties of the numericap stabiliey of the EM QNMs with respect to changes in the [hase-rovdiuiin.

An kbtfresting question is the results obtained usigf uhix method for $a>I$ or the so cwlkgd naked singuuarity refime. Preliminary rfsults fhow rhat the iethpd is applicable in this cawe as well aud rhe results will bz published zlsewhgre.

Ackmowledgements
================

The authorr wohld like to thank Emzvuele Berti for aisbusshon of the numerical valuef of the XM QNK frequdncigs obtayned throuhh Leaver’s method which aas ilpmrtant for the comparison of our method wmvh this alreacy wenl-establnshed pne.

P.F. is deepjy indebted to Irina Cref’iexa, Alexei Atarobiisky, Bruno Cjppi and Michdll Sazhin foc intereseing comnents ava suggestions sbout the relations of EM QNM to BH pmysicr and the possibnuuties to discovrr rucr loves ig astrophysicdl ocseffatiovs.

The authovs xre slso thankful to Luchano Rezzolla for drawonn our attgntion to the referencrs [@headon0] and [@GW_] ajd to Shehar Hpd sor drawing our attention to tge referejcef [@extr1] and [@evtr2].

Tmis wrticle wax supported by the Foundation “Theoreticel and Computational Phtsics and Astrophysncx”, by the Bukgariwn Nationdl Scientific Fund ubder contracts DO-1-872, DO-1-895, DO-02-136, and Sofia Univsrsity Sciejtific Fund, contract 185/26.04.2010.

Author Contributions
====================

P.F. posed the problem of the evalyation of the EM QHMs pf
persistent for the modes with $m=0,1,2$ and observation the important What are the has look in astrophysical Also interesting is while the $\epsilon$-method leads to significant of the frequencies $\omega_{m,n}$, it affects much less the separation parameter $E_{m,n}$ which for the first time was obtained directly as a solution of the two-dimensional without prior for Another general result is that the confluent Heun function proved to be an effective tool for problems. Even though its <span style="font-variant:small-caps;">maple</span> realization still many flaws, its precision to be good enough to the published results, also the we were able reveal new properties of the numerical stability of the EM QNMs with respect to changes in the An interesting the results using method $a>M$ or the naked singularity regime. Preliminary results show is applicable in this case as well and results will published elsewhere. Acknowledgements ================ The authors like to thank Emanuele Berti for discussion of numerical values of the EM QNM frequencies obtained through Leaver’s method which was important for of our method with already well-established one. is indebted Irina Alexei Starobinsky, Coppi and Michail Sazhin for interesting comments and suggestions about the of EM QNM to BH physics and the possibilities to modes astrophysical observations. The are also thankful to Rezzolla drawing our attention to [@headon0] [@GW_] Hod drawing attention to the references and [@extr2]. This article was by the Foundation “Theoretical by the Bulgarian National Scientific Fund under contracts DO-1-895, DO-02-136, and Sofia University Scientific Fund, 185/26.04.2010. Author Contributions ==================== P.F. posed the problem of the evaluation of EM QNMs
 persistent for the modes with $M=0,1,2$ and $l=1,2$. This oBservAtiOn rAiSes tHe imPortant questioN: what Are the electromagnetic QnMs foR wHIch oNE hAs to lOok in asTRoPHYsiCaL dAta. alSO iNtereStiNg is thaT while the $\ePsiLoN$-method leads TO sIgnificant ChaNges of the freQueNcies $\oMeGa_{m,N}$, It affEctS much Less thE SeparaTion paramEtER $E_{m,n}$ whICh here fOR ThE firSt time was obtained DIrECtly as a solutioN of the TwO-DiMENsiOnaL system witHoUt any PRior appROxIMATioNS for it.

Another General resuLT is That thE cOnfLUent HeUn funCtIOn pRoved to be an EffeCtive tool For phySIcal proBLems. EveN thougH itS <spAn stYLe="FoNt-vArIAnt:SMaLl-cAPs;">mAple</span> ReAlIzatiOn stILL HAs maNy fLaws, Its prEcision proved To bE gooD EnoUgh to RepeaT the AlReady PublisHed reSuLts, and also studyIng tHe solutioNs, wE wEre AbLe to rEVeal neW prOpeRties of The numeRIcaL sTABIlIty of the EM QNMs with ReSPEcT to changEs in thE PhAsE-ConditioN.

AN inTereSTIng quEstiON iS the resuLts obtAInEd Using thIs Method FoR $a>M$ Or tHe so cALled Naked sIngulariTy regIMe. Preliminary rESults show that THe METhOD is aPplIcable in thiS casE As weLl anD ThE reSUlts wIll be PuBLiSHed elsewhere.

AcknowlEdGementS
================

The aUthors would liKe to thank EMANUele BertI for DIsCUssion of the numEricaL values of tHE EM QNM frEquenCies obtaIned throuGH leaver’s mEthOd wHicH waS IMpOrtant for the cOMPariSoN of our mEthOd with tHis AlrEadY weLl-EstablishEd one.

P.F. iS dEePlY iNdeBted tO irina AreF’iEva, alExeI StarOBinsky, bruno coppI aNd mIchAil SazhIN fOR InteReStIng cOmmEnTs and SuggEStiOns abouT the relatIonS Of EM qNm tO BH physIcs and the possIbIlities to dIsCovEr such MODes in astRophysical observations.

THE authorS arE also ThanKful to LucIanO RezzoLla FOr drawIng our AttenTiOn tO THe refEREnCes [@HeAdon0] and [@GW_] aND To SHahar hoD for Drawing Our attention to the rEFerEnces [@extr1] and [@eXtr2].

this ARTiCle WAs SUppOrTEd bY THe Foundation “TheOretical anD COMpUtational PHYsiCs And AstrOphysicS”, by thE bulgariAn NationaL ScientifIc fund UNDer Contracts Do-1-872, DO-1-895, DO-02-136, and sofia UnivERsity sCiEntifIc FUnd, conTrAct 185/26.04.2010.

authoR ContrIButIons
====================

P.f. posed ThE problEm of tHe EvaluatiOn of the EM QNMs of
 persistent for the modeswith $m=0, 1,2$and  $l =1 ,2$.  Thi s observationr aise s the important questi on: W ha t  are th e ele ctromag n et i c  QN Ms  f orwh i ch  onehas  to loo k in astro phy si cal data. Al s ointerestin g i s that while  th e $\ep si lon $ -meth odleads  to si g nifica nt change so f thef requenc i e s$\om ega_{m,n}$, it af f ec t s much less th e sepa ra t io n par ame ter $E_{m, n} $ whi c h heref or t h e f i rst time wasobtained di r ect ly asasol u tion o f the  t w o-d imensionalsyst em withou t anyp rior ap p roximat ions f orit. 

An o th er  ge ne r alr es ult isthat the  c on fluen t He u n f unct ion  pro ved t o be an effec tiv e to o l f or ph ysica l pr ob lems.  Eventhoug hits <span style ="fo nt-varian t:s ma ll- ca ps;"> m aple</ spa n>realiza tion st i llha s m an y flaws, its preci si o n  p roved to  be go o den o ugh to r ep eat  the a lread y pu b li shed res ults,a nd  a lso stu dy ing th esol uti ons,w e we re abl e to rev eal n e w properties o f  the numerica l  s t a bi l ityofthe EM QNMs  wit h  res pect to  ch a ngesin th ep ha s e-condition.

An in te restin g que stion is theresults ob t a i ned usin g th i sm ethod for $a>M $ orthe so cal l ed naked  sing ularityregime. P r e liminary  re sul tssho w th at the method i s ap pl icableinthis ca seaswel l a nd  the resu lts will  b epu bl ish ed el s ewhere.
A ckn ow led gemen t s
==== ===== ==== == =The  author s  w o u ld l ik eto t han kEmanu eleB ert i for d iscussion  of thenu me rical v alues of theEM  QNM frequ en cie s obta i n ed throu gh Leaver’s method whic h  was im por tantforthe compa ris on ofour method  withthisal rea d y  well - e st abl is hed one.
P . F.is de ep ly i ndebted  to Irina Aref’iev a , A lexei Starobi nsk y, B r u no  Co p pi and  M i cha i l  Sazhin for int eresting c om m en ts and sug g est io ns abou t the r elati o ns of E M QNM toBH physic sandt h e p ossibiliti es to di scover su c h mod e sin as tro physic al  ob serva tions. 
Th e aut hors a re  alsothank fu l to Luc iano Rezzolla for drawi ng our  atte nti on to the  re f ere nces [@he adon 0] and [@G W_]  an d toSha h ar Ho d fo r  d raw i ng ou r at t ention to th e r e f er ences [@ext r 1 ]  an d [@e xtr 2 ].

Th is a rticle was suppor t ed by the Foun dati o n  “T heo r etic al  and Computati ona lP h ysics an dAstrophysic s”, by t he Bulga rian N ationa l Scien t i fi c  Fundunde r c ontractsDO- 1- 8 72, DO- 1- 89 5 , DO-0 2-13 6,  and S ofia U n iver s i ty Scientific Fu nd, c o n tract 185 /26.0 4. 2010.
A utho r Contribu tions
===== ====== ==== ===== 

P.F.po sed th e p ro blem of th e  evaluati on of  the EM  Q NMsof
 persistent_for the_modes with $m=0,1,2$ and_$l=1,2$. This_observation_raises the_important_question: What are_the electromagnetic QNMs_for which one has_to look in_astrophysical_data. Also interesting is that while the $\epsilon$-method leads to significant changes of the_frequencies_$\omega_{m,n}$, it_affects_much_less the separation parameter $E_{m,n}$_which here for the first_time was_obtained directly as a solution of the two-dimensional_system_without any prior_approximations for it.

Another general result is that the confluent_Heun function proved to be an_effective tool for_physical_problems._Even though its <span_style="font-variant:small-caps;">maple</span> realization still has many flaws,_its precision proved to be good_enough to repeat the already published results,_and also studying the solutions, we_were able to reveal new_properties of_the numerical stability of the_EM QNMs with_respect to_changes in the_phase-condition.

An interesting question is the results_obtained using this_method for $a>M$ or the so_called_naked singularity regime._Preliminary_results_show that_the method is_applicable_in this_case_as well and the results will_be_published elsewhere.

Acknowledgements
================

The authors would like to thank_Emanuele Berti for discussion_of_the numerical values of_the EM QNM frequencies obtained_through Leaver’s method which was important_for the_comparison of_our method with this already well-established one.

P.F. is deeply indebted to_Irina Aref’ieva, Alexei Starobinsky, Bruno Coppi_and Michail Sazhin for_interesting comments_and_suggestions about the_relations_of EM_QNM to BH physics and the possibilities_to discover_such modes in astrophysical observations.

The authors_are also thankful to_Luciano_Rezzolla for drawing our attention to_the references [@headon0] and [@GW_] and_to Shahar Hod for drawing_our_attention_to the references [@extr1] and_[@extr2].

This article was supported by the_Foundation “Theoretical and_Computational Physics and Astrophysics”, by the Bulgarian_National_Scientific Fund under contracts DO-1-872, DO-1-895,_DO-02-136,_and Sofia University Scientific Fund, contract_185/26.04.2010.

Author_Contributions
====================

P.F._posed the problem of the_evaluation of the EM QNMs of