Electronic Structure of disordered alloys, surfaces and interfaces.”, Kluwer Academic Publishers (1997). A. Chattopadhyay, S. Das Sarma and A.J. Millis, Phys. Rev. Lett. [**87**]{}, 227202 (2001). In the case of $S=1/2$, this corresponds to a singlet state $E_{S}=-\frac{3}{4}J$ with degeneracy 1 and a triplet state $E_{T}=+\frac{1}{4} J$ with degeneracy 3. In GaAs the volume of the unit cell is $v=a_{0}^{3}/4$ where $a_{0}=0.565$ nm. For simplicity, our calculations are performed on a simple cubic lattice, thus the lattice spacing we take is $a_{1}=(a_{0}^{3}/4)^{1/3}$. Assuming a larger $n_{h}$ of order $0.3-0.4$ leads to a value of $J$ weakly affected (only within $10 \%$). J. Okabayashi [*et al.*]{}, Phys. Rev. B [**58**]{}, R4211 (1998).
 ---
abstract: |
    Large scale tensors, including large scale Hankel tensors, have many applications in science and engineering. In this paper, we propose an inexact curvilinear search optimization method to compute Z- and H-eigenvalues of $m$th order $n$ dimensional Hankel tensors, where $n$ is large. Owing to the fast Fourier transform, the computational cost of each iteration of the new method is about $\mathcal{O}(mn\log(mn))$. Using the Cayley transform, we obtain an effective curvilinear search scheme. Then, we show that every limiting point of iterates generated by the new algorithm is an eigen-pair of Hankel tensors. Without the assumption of a second-order sufficient condition, we analyze the linear convergence rate of iterate sequence by the Kurdyka-[Ł]{}ojasiewicz property. Finally, numerical experiments for Hankel tensors, whose dimension may up to one million, are reported to show the efficiency of the proposed curvilinear search method.

    [**Key words:**]{} Cayley transform, curvilinear search, eigenvalue, fast Fourier transform, Hankel tensor, Kurdyka-[Ł]{}ojasiewicz property, large scale tensor.

    [**AMS subject classifications (2010):**]{} 15A18, 15A69, 65F15, 65
Electronic Structure of disordered alloys, surfaces and interfaces. ”, Kluwer Academic Publishers (1997). A. Chattopadhyay, S. Das Sarma and A.J. Millis, Phys. Rev. Lett. [ * * 87 * * ] { }, 227202 (2001). In the subject of $ S=1/2 $, this correspond to a singlet state $ E_{S}=-\frac{3}{4}J$ with degeneracy 1 and a triplet country $ E_{T}=+\frac{1}{4 } J$ with degeneracy 3. In GaAs the volume of the whole cellular telephone is $ v = a_{0}^{3}/4 $ where $ a_{0}=0.565 $   nm. For simplicity, our calculations are performed on a simple cubic wicket, thus the lattice spacing we take is $ a_{1}=(a_{0}^{3}/4)^{1/3}$. simulate a larger $ n_{h}$ of order $ 0.3 - 0.4 $ leads to a value of $ J$ weakly affected (only within $ 10 \%$). J. Okabayashi [ * et al. * ] { }, Phys. Rev. B [ * * 58 * * ] { }, R4211 (1998). 
  --- 
 abstract: | 
     bombastic scale tensors, including big scale Hankel tensors, have many applications in science and technology. In this paper, we aim an inexact curvilinear search optimization method to compute Z- and H - eigenvalue of $ m$th order $ n$ dimensional Hankel tensors, where $ n$ is large. Owing to the fast Fourier transform, the computational cost of each iteration of the new method is about $ \mathcal{O}(mn\log(mn))$. Using the Cayley transform, we obtain an effective curvilinear search scheme. Then, we show that every limiting point of iterates generated by the new algorithm is an eigen - couple of Hankel tensor. Without the assumption of a second - decree sufficient circumstance, we analyze the linear convergence rate of iterate sequence by the Kurdyka-[Ł]{}ojasiewicz property. Finally, numerical experiments for Hankel tensor, whose dimension may up to one million, are reported to show the efficiency of the proposed curvilinear search method. 

     [ * * Key actor's line :* * ] { } Cayley transform, curvilinear search, eigenvalue, fast Fourier transform, Hankel tensor, Kurdyka-[Ł]{}ojasiewicz property, large scale tensor. 

     [ * * AMS subject categorization (2010 ): * * ] { } 15A18, 15A69, 65F15, 65
Eleftronic Structure of disurdered alloys, surfacev and jnterfacds.”, Kluwer Academic Publisherd (1997). A. Chqttopadhyay, S. Das Sarmx and A.J. Lillis, Pyys. Cev. Lett. [**87**]{}, 227202 (2001). In tis case of $S=1/2$, tgls coxrxsponds to a sikglet state $E_{S}=-\frac{3}{4}J$ with geeeueracy 1 and a triplet state $E_{T}=+\frac{1}{4} J$ with drgfneracy 3. In GaWs tnq vomlmt of the unit cell is $v=a_{0}^{3}/4$ where $a_{0}=0.565$ nj. For spmplicity, our calvulations are performed on a slmple cubic latticf, thus the oattyxe spacing wd take is $a_{1}=(a_{0}^{3}/4)^{1/3}$. Assuming z larger $n_{h}$ of order $0.3-0.4$ leads to x valbe of $J$ weajlt avxected (only withpn $10 \%$). J. Okabayashi [*et dl.*]{}, Phys. Rev. B [**58**]{}, R4211 (1998).
 ---
absbract: |
    Large scale tensors, including large ssale Hankal tensors, have mant qppliwatimns kb sziehcx ahd englnexring. In thjs paper, we propose an inexact ctgfilinear seadch opeiiization method to compute Z- and H-eigennaluss of $m$th order $n$ dimenwional Hankel tensors, where $n$ ys large. Owing to the fast Fourier transform, the wompuvagiouql corr lf each iteration of the new method is about $\iztncsl{O}(mn\log(mn))$. Usikg the Cayley tranxflrk, we obtain an effeccjvs curvilinear searfh scheie. Thwn, we shor thst every limiting point of uterates genvratwd by the new algoxithm is an zigen-psir og Hankel tensors. Withouc the zssumption lf a secoha-order sufficieng cpngition, wt analyze the linewr converjence rate ow itgrate sqquence by the Kurdyka-[Ł]{}ojasiewicz pgoperjy. Findlly, numerlcal experiments for Hankel tensors, whose dimgnshon may up co one million, are weported to shpw the zfficidncy of thv proposev curvilineaw search methmf.

    [**Key worvs:**]{} Cayley trabsfoem, curvkuinear search, rigenvaluv, yast Fouruer transform, Hankcl tevaor, Kurdyka-[Ł]{}ojasntwixz property, larbe rcaje txnsor.

    [**AMS subjewt cuasroficagions (2010):**]{} 15A18, 15A69, 65N15, 65
Electronic Structure of disordered alloys, surfaces and Academic (1997). A. S. Das Sarma Lett. 227202 (2001). In case of $S=1/2$, corresponds to a singlet state $E_{S}=-\frac{3}{4}J$ degeneracy 1 and a triplet state $E_{T}=+\frac{1}{4} J$ with degeneracy 3. In GaAs volume of the unit cell is $v=a_{0}^{3}/4$ where $a_{0}=0.565$ nm. For simplicity, our are on simple lattice, thus the lattice spacing we take is $a_{1}=(a_{0}^{3}/4)^{1/3}$. Assuming a larger $n_{h}$ of order $0.3-0.4$ to a value of $J$ weakly affected (only $10 \%$). J. Okabayashi al.*]{}, Phys. Rev. B [**58**]{}, (1998). abstract: | scale including scale Hankel tensors, many applications in science and engineering. In this paper, we propose an inexact curvilinear search optimization method compute Z- of $m$th $n$ Hankel where $n$ is to the fast Fourier transform, the each iteration of the new method is about Using the transform, we obtain an effective curvilinear scheme. Then, we show that every limiting point iterates generated by the new algorithm is an eigen-pair of Hankel tensors. Without the assumption second-order sufficient condition, we the linear convergence of sequence the property. Finally, experiments for Hankel tensors, whose dimension may up to one million, reported to show the efficiency of the proposed curvilinear search words:**]{} transform, curvilinear search, fast Fourier transform, Hankel Kurdyka-[Ł]{}ojasiewicz large scale tensor. [**AMS (2010):**]{} 15A69,
Electronic Structure of disoRdered alloYs, surFacEs aNd InteRfacEs.”, Kluwer AcademIC PubLishers (1997). A. Chattopadhyay, S. das SaRmA And A.j. miLlis, PHys. Rev. LETt. [**87**]{}, 227202 (2001). iN The CaSe Of $S=1/2$, ThIS cOrresPonDs to a siNglet state $e_{S}=-\fRaC{3}{4}J$ with degeneRAcY 1 and a triplEt sTate $E_{T}=+\frac{1}{4} J$ wIth DegeneRaCy 3. IN gaAs tHe vOlume Of the uNIt cell Is $v=a_{0}^{3}/4$ where $A_{0}=0.565$ nM. for simPLicity, oUR CaLculAtions are performeD On A Simple cubic latTice, thUs THe LATtiCe sPacing we taKe Is $a_{1}=(a_{0}^{3}/4)^{1/3}$. ASSuming a LArGER $N_{h}$ oF Order $0.3-0.4$ leads to a Value of $J$ weaKLy aFfecteD (oNly WIthin $10 \%$). J. okabaYaSHi [*eT al.*]{}, Phys. Rev. B [**58**]{}, r4211 (1998).
 ---
absTract: |
    LargE scale TEnsors, iNCluding Large sCalE HaNkel TEnSoRs, hAvE ManY ApPliCAtiOns in sciEnCe And enGineERING. In tHis PapeR, we prOpose an inexacT cuRvilINeaR searCh optImizAtIon meThod to CompuTe z- and H-eigenvalueS of $m$Th order $n$ dImeNsIonAl hankeL TensorS, whEre $N$ is largE. Owing tO The FaST fOuRier transform, the coMpUTAtIonal cosT of eacH ItErATion of thE nEw mEthoD IS abouT $\matHCaL{O}(mn\log(mN))$. Using THe caYley traNsForm, we ObTaiN an EffecTIve cUrviliNear searCh schEMe. Then, we show thAT every limitinG PoINT oF IterAteS generated bY the NEw alGoriTHm Is aN Eigen-Pair oF HANkEL tensors. Without the aSsUmptioN of a sEcond-order sufFicient conDITIon, we anaLyze THe LInear convergenCe ratE of iterate SEquence bY the KUrdyka-[Ł]{}oJasiewicz PROperty. FiNalLy, nUmeRicAL ExPeriments for HANKel tEnSors, whoSe dImensioN maY up To oNe mIlLion, are rePorted to ShOw ThE eFfiCiencY Of the proPoSed CuRviLineaR Search MethoD.

    [**Key WoRdS:**]{} cayLey tranSFoRM, CurvIlInEar sEarCh, EigenValuE, FasT FourieR transforM, HaNKel tEnSoR, KurdykA-[Ł]{}ojasiewicz pRoPerty, large ScAle Tensor.

    [**ams subject Classifications (2010):**]{} 15A18, 15A69, 65F15, 65
Electronic Structure of di sordered a lloys , s urf ac es a nd i nterfaces.”, K l uwer  Academic Publishers ( 1997) .A . Ch a tt opadh yay, S. Da s Sar ma  a ndA. J .Milli s,Phys. R ev. Lett.[** 87 **]{}, 22720 2  ( 2001). Inthe  case of $S= 1/2 $, thi scor r espon dsto asingle t  state  $E_{S}=- \f r ac{3}{ 4 }J$ wit h de gene racy 1 and a trip l et state $E_{T}=+ \frac{ 1} { 4} J $ w ith  degenerac y3. In GaAs th e  v o l u meo f the unit ce ll is $v=a_ { 0}^ {3}/4$  w her e  $a_{0 }=0.5 65 $  nm . For simpl icit y, our ca lculat i ons are perform ed ona s imp le c u bi clat ti c e,t hu s t h e l attice s pa ci ng we  tak e i s  $a_ {1} =(a_ {0}^{ 3}/4)^{1/3}$.  As sumi n g a  larg er $n _{h} $of or der $0 .3-0. 4$  leads to a val ue o f $J$ wea kly  a ffe ct ed (o n ly wit hin  $1 0 \%$).  J. Oka b aya sh i [ *e t al.*]{}, Phys. R ev . B[**58**] {}, R4 2 11  ( 1 998).
 - -- 
ab stra c t : |
   L a rg e scaletensor s ,in cluding  l arge s ca leHan kel t e nsor s, hav e many a pplic a tions in scien c e and enginee r in g .  I n  thi s p aper, we pr opos e  aninex a ct  cu r vilin ear s ea r ch optimization method  t o comp ute Z - and H-eigen values of$ m $ th order  $n$ di m ensional Hanke l ten sors, wher e  $n$ islarge . Owingto the fa s t  Fourier  tr ans for m,t h ecomputational c ostof  each i ter ation o f t henew  me th od is abo ut $\mat hc al {O }( mn\ log(m n ))$. Usi ng  th eCay ley t r ansfor m, we  obt ai na n e ffectiv e  c u r vili ne ar  sea rch  s cheme . Th e n,we show  that eve ryl imit in gpoint o f iterates ge ne rated by t he  ne w algo r i thm is a n eigen-pair of Hankelt ensors.  Wi thout  the  assumpti onof a s eco n d-orde r suff icien tcon d i tion, w eana ly ze the lin e a r c onver ge ncerate of  iterate sequenceb y t he Kurdyka-[Ł ]{} ojas i e wi czp ro p ert y. Fin a l ly, numerical e xperiments  f o rHankel ten s ors ,whose d imensio n may up to o ne millio n, are re po rted t o s how the ef ficiency  of the p r opose d  c urvil ine ar sea rc h m ethod .

  [** Key w ords:* *] {} Cay ley t ra nsform,curvilinear search, eig envalu e, fa stFourier t ran s for m, Hankel  ten sor, Kurdy ka- [Ł] {}oja sie w icz p rope r ty , l a rge s cale tensor.
  [* * A MS  subject cl a s s ifi catio ns( 2010): **]{ } 15A18, 15A69, 6 5 F15, 65
Electronic Structure_of disordered_alloys, surfaces and interfaces.”,_Kluwer Academic_Publishers_(1997). A._Chattopadhyay,_S. Das Sarma_and A.J. Millis,_Phys. Rev. Lett. [**87**]{},_227202 (2001). In_the_case of $S=1/2$, this corresponds to a singlet state $E_{S}=-\frac{3}{4}J$ with degeneracy 1 and_a_triplet state_$E_{T}=+\frac{1}{4}_J$_with degeneracy 3. In GaAs_the volume of the unit_cell is_$v=a_{0}^{3}/4$ where $a_{0}=0.565$ nm. For simplicity, our calculations are_performed_on a simple_cubic lattice, thus the lattice spacing we take is_$a_{1}=(a_{0}^{3}/4)^{1/3}$. Assuming a larger $n_{h}$ of_order $0.3-0.4$ leads_to_a_value of $J$ weakly_affected (only within $10 \%$). J._Okabayashi [*et al.*]{}, Phys. Rev. B_[**58**]{}, R4211 (1998).
 ---
abstract: |
  _ Large scale tensors, including large_scale Hankel tensors, have many_applications in_science and engineering. In this_paper, we propose_an inexact_curvilinear search optimization_method to compute Z- and H-eigenvalues_of $m$th order_$n$ dimensional Hankel tensors, where $n$_is_large. Owing to_the_fast_Fourier transform,_the computational cost_of_each iteration_of_the new method is about $\mathcal{O}(mn\log(mn))$._Using_the Cayley transform, we obtain an effective_curvilinear search scheme. Then,_we_show that every limiting_point of iterates generated by_the new algorithm is an eigen-pair_of Hankel_tensors. Without_the assumption of a second-order sufficient condition, we analyze the linear_convergence rate of iterate sequence by_the Kurdyka-[Ł]{}ojasiewicz property. Finally,_numerical experiments_for_Hankel tensors, whose_dimension_may up_to one million, are reported to show_the efficiency_of the proposed curvilinear search method.

_   [**Key_words:**]{}_Cayley transform, curvilinear search, eigenvalue, fast_Fourier transform, Hankel tensor, Kurdyka-[Ł]{}ojasiewicz property,_large scale tensor.

  __[**AMS_subject classifications (2010):**]{} 15A18, 15A69,_65F15, 65