SE) for neutron star magnetospheres [@Glampedakis2014; @Fujisawa2014; @Pili2015; @Akgun2016; @Kojima2017; @Kojima2018; @Kojima2018b; @Akgun2018a]. @Akgun2017 performed magneto-thermal evolutions coupling the crustal magnetic field at the stellar surface with an exterior equilibrium solution. The results showed that large twists grow in the magnetosphere up to a critical point beyond which no stable equilibrium solutions where found. A more detailed analysis by @Akgun2018a showed that, for sufficiently large twists, the solutions of the GSE are degenerate with several possible configurations of different energies but matching boundary conditions at the surface. This suggests the possibility of an unstable branch of the solutions and, thus, a possible explanation for the occurrence of bursts and GFs. In this work we explore this possibility by performing three-dimensional (3D) numerical simulations of the equilibrium models in @Akgun2018a. We asses their stability properties and their potential as candidates for transient magnetar phenomenology.

This work is organised as follows. In section \[sec:physics\] we review and discuss the physics involved in magnetars relevant to the processes that we want to study. In section \[sec:ff\_equations\] we briefly review the equations of force-free electrodynamics (FFE) implemented for simulations conducted on the infrastructure of the `Einstein Toolkit` (supplemented by appendix \[sec:augmented\_system\]). A detailed description of the derivation of initial models according to @Akgun2018a is given in section \[sec:twisted\_models\]. In section \[sec:simulations\] we present the numerical setup of our simulations as well as the outcome of the conducted 3D force-free simulations of twisted magnetospheres (reviewing details on maintaining the force-free regime in appendix \[sec:ff\_breakdown\]). The observed rapid dissipation of electromagnetic energy through the magnetar crust is interpreted and related to observable quantities, such as luminosity estimates, shear stresses on the stellar crust, and opacity models, in section \[sec:discussion\]. Along this paper we use Gaussian units in CGS, except for section \[sec:ff\_equations\] in which we use Heaviside-Lorentz with geometrised units ($G=c=M_\odot=1$). For convenience we express current densities in A m$^{-2}$
SE) for neutron star magnetospheres [ @Glampedakis2014; @Fujisawa2014; @Pili2015; @Akgun2016; @Kojima2017; @Kojima2018; @Kojima2018b; @Akgun2018a ]. @Akgun2017 performed magneto - thermal evolutions coupling the crustal charismatic sphere at the stellar surface with an exterior chemical equilibrium solution. The solution showed that large wind develop in the magnetosphere up to a critical point beyond which no static chemical equilibrium solutions where find. A more detailed analysis by @Akgun2018a testify that, for sufficiently large twists, the solutions of the GSE are degenerate with several possible configuration of different energies but matching boundary condition at the surface. This suggests the possibility of an unstable outgrowth of the solutions and, thus, a possible explanation for the occurrence of bursts and GFs. In this workplace we explore this possibility by performing three - dimensional (3D) numerical simulations of the equilibrium models in @Akgun2018a. We asses their stability properties and their potential as candidates for transient magnetar phenomenology. 

 This work is organised as follows. In section   \[sec: physics\ ] we review and discuss the purgative involve in magnetars relevant to the processes that we desire to analyze. In incision   \[sec: ff\_equations\ ] we briefly review the equations of force - free electrodynamics (FFE) implemented for simulations conducted on the infrastructure of the 'Einstein Toolkit '(append by appendix \[sec: augmented\_system\ ]). A detailed description of the derivation of initial models according to @Akgun2018a is given in section   \[sec: twisted\_models\ ]. In section   \[sec: simulations\ ] we award the numerical setup of our simulations as well as the outcome of the conducted 3D force - complimentary simulations of twisted magnetospheres (reviewing details on keep the force - free regime in appendix \[sec: ff\_breakdown\ ]). The ascertained rapid dissipation of electromagnetic department of energy through the magnetar crust is interpreted and related to observable quantities, such as luminosity estimates, shear stress on the stellar crust, and opacity model, in section   \[sec: discussion\ ]. Along this composition we use Gaussian unit in CGS, except for section   \[sec: ff\_equations\ ] in which we practice Heaviside - Lorentz with geometrised units ($ G = c = M_\odot=1 $). For convenience we press out current densities in A m$^{-2}$
SE) vor neutron star magnetorpheres [@Glampedcjis2014; @Fuoisawa2014; @Lili2015; @Akgjn2016; @Kojima2017; @Kojima2018; @Kojima2018b; @Akgnn2018a]. @Qkgun2017 performed magneto-theroal evolunions couplinj the crustal majhetic flzld af the wtellar surfacg with an exderior equilibsijm solution. The results showed that lwrge twoshs grow in the magmqtoslhere up to a critical point beyons which no stable eqiilibrium solutions where vounf. A more detailed wnalysis by @Akgtb2018a showed thxt, for sufficiently latge twists, the solutions of the ESE axe degenerajz wihv several pissibje configurabpons of differrnt energies bmt mavchibg boundary conditionv at the surface. Tris suggevta the possibility od an onstatle ceanzh kf tge solktikns and, thhs, a possiboe explanation for uhe iccurrence of burstf wnd GFs. In this work we explore this povsigility by performing theee-dimensional (3D) numetical simujations of the equilibrium models in @Akgun2018a. We asves tiekr wtwcullty properties and their potential as candidaess fpr transient mcgnetar phenomenpllgu.

Jhis work is otganiseb aa follows. In sectiln \[sec:phisics\] qe review and discuss the physics involvwd in magnetcrs relevant to the pxocesses thac we wsnt tp study. In section \[sec:ff\_zquatikns\] we brievly revies the equations ow fprwe-free electrodynamics (FFE) implemenved fpr simuuatipns cogducted on the ltfrastructure of tje `Eiusteit Toolkit` (dupplemented by appendix \[sec:augmented\_system\]). A datapled descxiptiok of the derivaeion of initiak modelf accurding to @Zkgun2018a ms given in fection \[sec:twivjed\_models\]. In vection \[sqc:sinularions\] wd present the nimerical setup of oyr simulations as xell zs the outcome if uhw conducted 3D gorze-fwev smmulaehons of twisded oagvrtospferes (rtriecivg drtails on maintaininc ths force-free regime ik appendiz \[sec:ff\_bweakdown\]). The pbserved rapid disdipatmon of elecyroiagnetic energy through the mafnetar crkst is interpreeed qnd related co observable quantities, such as luminosmty estimates, shear sttesses on the stellar grust, and opecity iodels, in section \[sec:discussiob\]. Along this papev we use Gaussian unita in CCS, exfept for section \[sec:ff\_equations\] in which we use Heaviside-Lorentz with geomwtrisev tnits ($G=c=M_\odkt=1$). Fpr cotvzniznce we expcers current densinies in A m$^{-2}$
SE) for neutron star magnetospheres [@Glampedakis2014; @Fujisawa2014; @Kojima2017; @Kojima2018b; @Akgun2018a]. performed magneto-thermal evolutions at stellar surface with exterior equilibrium solution. results showed that large twists grow the magnetosphere up to a critical point beyond which no stable equilibrium solutions found. A more detailed analysis by @Akgun2018a showed that, for sufficiently large twists, solutions the are with several possible configurations of different energies but matching boundary conditions at the surface. This suggests possibility of an unstable branch of the solutions thus, a possible explanation the occurrence of bursts and In work we this by three-dimensional (3D) numerical of the equilibrium models in @Akgun2018a. We asses their stability properties and their potential as candidates for magnetar phenomenology. is organised follows. section we review and physics involved in magnetars relevant to we want to study. In section \[sec:ff\_equations\] we review the of force-free electrodynamics (FFE) implemented for conducted on the infrastructure of the `Einstein Toolkit` by appendix \[sec:augmented\_system\]). A detailed description of the derivation of initial models according to @Akgun2018a in section \[sec:twisted\_models\]. In \[sec:simulations\] we present numerical of simulations well as outcome of the conducted 3D force-free simulations of twisted magnetospheres (reviewing on maintaining the force-free regime in appendix \[sec:ff\_breakdown\]). The observed of energy through the crust is interpreted and to quantities, such as luminosity stresses the opacity in \[sec:discussion\]. Along this paper use Gaussian units in CGS, for section \[sec:ff\_equations\] in geometrised units ($G=c=M_\odot=1$). For convenience we express current in A m$^{-2}$
SE) for neutron star magnetospHeres [@GlampEdakiS2014; @FuJisAwA2014; @PilI2015; @AkgUn2016; @Kojima2017; @Kojima2018; @kOjimA2018b; @Akgun2018a]. @Akgun2017 performed MagneTo-THermAL eVolutIons couPLiNG The CrUsTal MaGNeTic fiEld At the stEllar surfaCe wItH an exterior eQUiLibrium solUtiOn. The results ShoWed thaT lArgE TwistS grOw in tHe magnETospheRe up to a crItICal poiNT beyond WHIcH no sTable equilibrium sOLuTIons where found. a more dEtAIlED AnaLysIs by @Akgun2018a ShOwed tHAt, for suFFiCIENtlY Large twists, thE solutions oF The gSE are DeGenERate wiTh sevErAL poSsible confiGuraTions of diFferenT EnergieS But matcHing boUndAry CondITiOnS at ThE SurFAcE. ThIS suGgests thE pOsSibilIty oF AN UNstaBle BranCh of tHe solutions anD, thUs, a pOSsiBle exPlanaTion FoR the oCcurreNce of BuRsts and GFs. In thiS worK we explorE thIs PosSiBilitY By perfOrmIng Three-diMensionAL (3D) nUmERICaL simulations of the eQuILIbRium modeLs in @AkGUn2018A. WE Asses theIr StaBiliTY PropeRtieS AnD their poTentiaL As CaNdidateS fOr tranSiEnt MagNetar PHenoMenoloGy.

This woRk is oRGanised as folloWS. In section \[sec:PHySICs\] WE revIew And discuss tHe phYSics InvoLVeD in MAgnetArs reLeVAnT To the processes that wE wAnt to sTudy. IN section \[sec:ff\_Equations\] wE BRIefly revIew tHE eQUations of force-Free eLectrodynaMIcs (FFE) imPlemeNted for sImulationS COnducted On tHe iNfrAstRUCtUre of the `EinstEIN TooLkIt` (supplEmeNted by aPpeNdiX \[seC:auGmEnted\_systEm\]). A detaiLeD dEsCrIptIon of THe derivaTiOn oF iNitIal moDEls accOrdinG to @AKgUn2018A Is gIven in sECtION \[sec:TwIsTed\_mOdeLs\]. in secTion \[SEc:sImulatiOns\] we presEnt THe nuMeRiCal setuP of our simulatIoNs as well as ThE ouTcome oF THe conducTed 3D force-free simulationS Of twistEd mAgnetOsphEres (revieWinG detaiLs oN MaintaIning tHe forCe-FreE REgime IN ApPenDiX \[sec:ff\_breaKDOwn\]). the obSeRved Rapid diSsipation of electroMAgnEtic energy thrOugH the MAGnEtaR CrUSt iS iNTerPREted and related tO observablE qUAnTities, such AS luMiNosity eStimateS, sheaR StresseS on the steLlar crust, AnD opaCITy mOdels, in secTion \[sec:dIscussion\]. aLong tHIs Paper We uSe GausSiAn uNits iN CGS, exCEpt For seCtion \[sEc:Ff\_equaTions\] In Which we uSe Heaviside-Lorentz with gEometrIsed uNitS ($G=c=M_\odot=1$). FOr cONveNience we eXpreSs current dEnsItiEs in A M$^{-2}$
SE) for neutron star magne tospheres[@Gla mpe dak is 2014 ; @F ujisawa2014; @ P ili2 015; @Akgun2016; @Koji ma201 7; @Koj i ma 2018;  @Kojim a 20 1 8 b;@A kg un2 01 8 a] . @Ak gun 2017 pe rformed ma gne to -thermal evo l ut ions coupl ing  the crustal  ma gnetic  f iel d  at t hestell ar sur f ace wi th an ext er i or equ i librium s ol utio n. The results sh o we d  that large tw ists g ro w  i n the  ma gnetospher eup to a criti c al p o int beyond whichno stable e q uil ibrium  s olu t ions w herefo u nd.  A more det aile d analysi s by @ A kgun201 8 a showe d that , f orsuff i ci en tly  l a rge tw ist s , t he solut io ns  of t he G S E a re d ege nera te wi th several po ssi blec onf igura tions  ofdi ffere nt ene rgies  b ut matching bou ndar y conditi ons  a t t he  surf a ce. Th issug gests t he poss i bil it y o fan unstable branch  o f th e soluti ons an d ,th u s, a pos si ble  exp l a natio n fo r  t he occur renceo fbu rsts an dGFs. I nthi s w ork w e  exp lore t his poss ibili t y by performin g  three-dimens i on a l  ( 3 D) n ume rical simul atio n s of  the eq uil i brium  mode ls in @Akgun2018a. We ass es  their  stab ility propert ies and th e i r  potenti al a s  c a ndidates for t ransi ent magnet a r phenom enolo gy.

Thi s work is o rganised  as  fo llo ws. I nsection \[sec : p hysi cs \] we r evi ew anddis cus s t heph ysics inv olved in  m ag ne ta rsrelev a nt to th epro ce sse s tha t  we wa nt to  stu dy .I n s ection\ [s e c :ff\ _e qu atio ns\ ]we br iefl y  re view th e equatio nso f fo rc e- free el ectrodynamics  ( FFE) imple me nte d fors i mulation s conducted on the infr a structu reof th e `E instein T ool kit` ( sup p lement ed byappen di x \ [ s ec:au g m en ted \_ system\]). A  de taile ddesc ription  of the derivation ofinitial model s a ccor d i ng  to @A k gun 20 1 8ai s  given in secti on \[sec:t wi s te d\_models\ ] . I nsection  \[sec: simul a tions\]  we prese nt the nu me rica l set up of oursimulati ons as we l l ast he  outc ome  of th econ ducte d 3D f o rce -free  simul at ions o f twi st ed magne tospheres (reviewing de tailson ma int aining th e f o rce -free reg imein appendi x \ [se c:ff\ _br e akdow n\]) .  T heo bserv ed r a pid dissi p at ion o felectromagn e t i c e nergy  th r ough t he m agnetar crust isi nterpreted and  rel a t edtoo bser va ble quantities , s uc h as lumin os ity estimat es, shea rs tress es onthe st ellar c r u st ,  and o paci tymodels, i n s ec t ion \[s ec :d i scussi on\] .Alongthis p a perw e  use Gaussian un its i n CGS,e xce pt fo rsection \[se c:ff\_equa tions\] inwhichwe u se He aviside -L orentz  wi th  geometris e d units ( $G=c= M_\odot =1 $).For  conve nien c e  we e xpre ss  cu rrent den s i ti e sin A m $^{- 2}$
SE) for_neutron star_magnetospheres [@Glampedakis2014; @Fujisawa2014; @Pili2015;_@Akgun2016; @Kojima2017;_@Kojima2018;_@Kojima2018b; @Akgun2018a]._@Akgun2017_performed magneto-thermal evolutions_coupling the crustal_magnetic field at the_stellar surface with_an_exterior equilibrium solution. The results showed that large twists grow in the magnetosphere up_to_a critical_point_beyond_which no stable equilibrium solutions_where found. A more detailed_analysis by_@Akgun2018a showed that, for sufficiently large twists, the_solutions_of the GSE_are degenerate with several possible configurations of different energies_but matching boundary conditions at the_surface. This suggests_the_possibility_of an unstable branch_of the solutions and, thus, a_possible explanation for the occurrence of_bursts and GFs. In this work we_explore this possibility by performing three-dimensional_(3D) numerical simulations of the_equilibrium models_in @Akgun2018a. We asses their_stability properties and_their potential_as candidates for_transient magnetar phenomenology.

This work is organised_as follows. In_section \[sec:physics\] we review and discuss the_physics_involved in magnetars_relevant_to_the processes_that we want_to_study. In_section \[sec:ff\_equations\]_we briefly review the equations of_force-free_electrodynamics (FFE) implemented for simulations conducted on_the infrastructure of the_`Einstein_Toolkit` (supplemented by appendix_\[sec:augmented\_system\]). A detailed description of_the derivation of initial models according_to @Akgun2018a_is given_in section \[sec:twisted\_models\]. In section \[sec:simulations\] we present the numerical setup of our_simulations as well as the outcome_of the conducted 3D_force-free simulations_of_twisted magnetospheres (reviewing_details_on maintaining_the force-free regime in appendix \[sec:ff\_breakdown\]). The_observed rapid_dissipation of electromagnetic energy through the_magnetar crust is interpreted_and_related to observable quantities, such as_luminosity estimates, shear stresses on the_stellar crust, and opacity models,_in_section \[sec:discussion\]._Along this paper we use_Gaussian units in CGS, except for_section \[sec:ff\_equations\] in which_we use Heaviside-Lorentz with geometrised units ($G=c=M_\odot=1$)._For_convenience we express current densities in_A_m$^{-2}$