 direction of the applied field. Thus the absolute value of $Q$ for a nanostrip in the skyrmion phase equals to the number of skyrmions supported.

The dependence of the topological charge on the curvature angle is shown in Fig.\[fig:Q\_vs\_a\] for nanostrips with different widths $L_y$.

\[htb!\]![ (Color online) Dependence of the ground state topological charge (Q) on the angle of curvature ($\phi_0$) for magnetic nanostructures with different sizes ($L\times L$) and $a=2nm$. The applied magnetic field is uniform along the x-axis (circles, closed triangles) or radial (open triangles). []{data-label="fig:Q_vs_a"}](fig3.jpg "fig:"){width="0.95\linewidth"}

We notice that $Q$ remains almost constant up to an angle $\phi_0\approx 100^0$ and further on it decreases smoothly from the initial value for planar nanostrips (skyrmion phase), to a smaller value close to zero for nanotubes, indicating that only a small fraction of the initial number of skyrmions are stabilized (mixed skyrmion-stripe phase). In contrast to this trend, when the applied field is radial, the topological charge is only weakly dependent on the curvature angle. This weak decay of $Q$ with curvature angle is due to the gradual reduction of the co-planar condition between the DMI vectors on each lattice site. A radial field for a cylindrical nanostrip is geometrically analogous to the case of a uniform field normal to a planar nanostrip and the conservation of the topological charge with curvature underlines the importance of the normality condition for the applied filed in the stabilization of skyrmions on any surface.

For planar nanostrips shown in Fig.\[fig:Skphase\]a ($Q=7.2$) and Fig.\[fig:Skphase\]d ($Q=37.8$), the values of $Q$ deviate weakly from integer values due to the misalignment of the moments located on the free boundaries of the sample[@roh13] and the thermal fluctuations inherent to the Monte Carlo method. For curved surfaces, however, the shape distortion of the skyrmions and their evolution to stripe-like structures is not characterized by integer values of $Q$, thus the calculation of $Q$ based on Eq.(\[eq:topol\_charge\]) assumes non-
direction of the applied field. Thus the absolute value of $ Q$ for a nanostrip in the skyrmion phase equal to the numeral of skyrmions supported. 

 The dependence of the topological charge on the curvature slant is shown in Fig.\[fig: Q\_vs\_a\ ] for nanostrips with different widths $ L_y$. 

 \[htb!\ ]! [ (Color online) Dependence of the background state topological charge (Q) on the slant of curvature ($ \phi_0 $) for charismatic nanostructures with different sizes ($ L\times L$) and $ a=2nm$. The apply magnetic field is uniform along the x - axis (circles, closed triangle) or radial (open triangles). [ ] { data - label="fig: Q_vs_a"}](fig3.jpg " fig:"){width="0.95\linewidth " } 

 We comment that $ Q$ remains almost constant up to an angle $ \phi_0\approx 100 ^ 0 $ and farther on it decreases smoothly from the initial value for planar nanostrips (skyrmion phase), to a smaller value close to zero for nanotubes, indicating that merely a small fraction of the initial number of skyrmions are stabilized (mixed skyrmion - stripe phase). In contrast to this trend, when the applied field is radial, the topological charge is only weakly dependent on the curvature angle. This weak decay of $ Q$ with curvature angle is due to the gradual reduction of the carbon monoxide - planar condition between the DMI vectors on each wicket web site. A radial field for a cylindrical nanostrip is geometrically analogous to the case of a uniform field normal to a planar nanostrip and the conservation of the topological mission with curvature underscore the importance of the normality condition for the applied filed in the stabilization of skyrmions on any open. 

 For planar nanostrips shown in Fig.\[fig: Skphase\]a ($ Q=7.2 $) and Fig.\[fig: Skphase\]d ($ Q=37.8 $), the values of $ Q$ deviate weakly from integer values due to the misalignment of the moments located on the free boundaries of the sample[@roh13 ] and the thermal fluctuations inherent to the Monte Carlo method. For curved surface, however, the form distortion of the skyrmions and their evolution to stripe - like structures is not characterized by integer value of $ Q$, thus the calculation of $ Q$ based on Eq.(\[eq: topol\_charge\ ]) assumes non-
 digection of the applied fleld. Thus the absolute talue or $Q$ for x nanostrip in the skyrmion 'hasw equqls to the number of syyrmions dupportee.

The eependence of the bjpolknical rharge on the corvature angne is shown in Fkg.\[yig:Q\_vs\_a\] for nanostrips with differene widthx $P_y$.

\[htb!\]![ (Color onjine) Qepehdence of the ground state topologjcal cherge (Q) on the amgle of curvature ($\phi_0$) for lagnftic nanostructured with diffgdene sizes ($L\timer L$) and $a=2nm$. The applies magnetic field is uniform alovg thz x-axis (cirxlws, fnosed trianjles) og radial (open trianglas). []{data-kabel="fig:Q_vs_a"}](fin3.jpg "hig:"){wudth="0.95\linewidth"}

We noticx that $Q$ remains almjst constdnc up to an angle $\phi_0\appeox 100^0$ dnd xurtfwr un jt dscreasfs amoothly fdom the inirial value for plansr banostrips (skgrmion prase), to a smaller value close to zero fmr hanotubes, indicating thqt only a small fractlon of thq initial number of skyrmions are stabilized (mixeg skycmkon-wtvipe phwse). In contrast to this trend, when the applieq foekd is radial, tme topological chatgf of only weakly depenbsnf on the curvature angle. Jhis wwak decay of $A$ with curvature angle is dye to the grcduql reduction of thz co-planar cundijion brtween the DMI vectors un ezch lattice site. A rzaial field for a cykitdrical nanostrip is geomeerically enaloyous to ghe vase os a uniforl field normal to a planwr nauostrhp and the conservation of the topological charge with vusvanure undexlines the importanse of the normclity couditiov for the zpplied filed in tre stabilizathln of skyrmimns on agy syrfaxe.

For puxnar nanostripx shown iu Fig.\[dig:Skphase\]a ($Q=7.2$) and Nig.\[fie:Akphase\]d ($Q=37.8$), the rclows of $Q$ deviatr wdakjy fcom igdeger values due to yhe mksalignment uf tne moments located ot ths free boundaries pf the samkle[@roh13] anq the thermal fluctuations inhegent vo the Montr Cwrlo method. For curved surfacea, however, thc shape distowtiok of the skyrmnons and their evolution to stripe-like svructures is not characrerized by integer rakues of $Q$, tius thq calculadion of $Q$ based on Ew.(\[eq:topol\_charge\]) asxumes non-
direction of the applied field. Thus the of for a in the skyrmion of supported. The dependence the topological charge the curvature angle is shown in for nanostrips with different widths $L_y$. \[htb!\]![ (Color online) Dependence of the ground topological charge (Q) on the angle of curvature ($\phi_0$) for magnetic nanostructures with sizes L$) $a=2nm$. applied magnetic field is uniform along the x-axis (circles, closed triangles) or radial (open triangles). []{data-label="fig:Q_vs_a"}](fig3.jpg We notice that $Q$ remains almost constant up an angle $\phi_0\approx 100^0$ further on it decreases smoothly the value for nanostrips phase), a smaller value to zero for nanotubes, indicating that only a small fraction of the initial number of skyrmions are (mixed skyrmion-stripe contrast to trend, the field is radial, charge is only weakly dependent on This weak decay of $Q$ with curvature angle due to gradual reduction of the co-planar condition the DMI vectors on each lattice site. A field for a cylindrical nanostrip is geometrically analogous to the case of a uniform field a planar nanostrip and conservation of the charge curvature the of the condition for the applied filed in the stabilization of skyrmions on surface. For planar nanostrips shown in Fig.\[fig:Skphase\]a ($Q=7.2$) and Fig.\[fig:Skphase\]d values $Q$ deviate weakly integer values due to misalignment the moments located on boundaries the thermal inherent the Monte Carlo method. curved surfaces, however, the shape of the skyrmions and is not characterized by integer values of $Q$, the calculation of $Q$ based on Eq.(\[eq:topol\_charge\]) non-
 direction of the applied fielD. Thus the abSolutE vaLue Of $q$ for A nanOstrip in the skyRMion Phase equals to the number Of skyRmIOns sUPpOrted.

the depeNDeNCE of ThE tOpoLoGIcAl chaRge On the cuRvature angLe iS sHown in Fig.\[fig:q\_Vs\_A\] for nanostRipS with differeNt wIdths $L_Y$.

\[hTb!\]![ (COLor onLinE) DepeNdence OF the grOund state ToPOlogicAL charge (q) ON tHe anGle of curvature ($\phi_0$) FOr MAgnetic nanostrUctureS wITh DIFfeRenT sizes ($L\timEs l$) and $a=2NM$. The appLIeD MAGneTIc field is unifOrm along the X-AxiS (circlEs, CloSEd triaNgles) Or RAdiAl (open trianGles). []{Data-label="Fig:Q_vs_A"}](Fig3.jpg "fIG:"){width="0.95\lInewidTh"}

WE noTice THaT $Q$ RemAiNS alMOsT coNStaNt up to an AnGlE $\phi_0\aPproX 100^0$ AND FurtHer On it DecreAses smoothly fRom The iNItiAl valUe for PlanAr NanosTrips (sKyrmiOn Phase), to a smaller ValuE close to zEro FoR naNoTubes, INdicatIng ThaT only a sMall fraCTioN oF THE iNitial number of skyrMiONS aRe stabilIzed (miXEd SkYRmion-strIpE phAse). IN COntraSt to THiS trend, whEn the aPPlIeD field iS rAdial, tHe TopOloGical CHargE is onlY weakly dEpendENt on the curvatuRE angle. This weaK DeCAY oF $q$ witH cuRvature anglE is dUE to tHe grADuAl rEDuctiOn of tHe CO-pLAnar condition betweeN tHe DMI vEctorS on each latticE site. A radiAL FIeld for a CyliNDrICal nanostrip is GeomeTrically anALogous to The caSe of a uniForm field NORmal to a pLanAr nAnoStrIP AnD the conservatION of tHe TopologIcaL charge WitH cuRvaTurE uNderlines The imporTaNcE oF tHe nOrmalITy conditIoN foR tHe aPplieD Filed iN the sTabiLiZaTIon Of skyrmIOnS ON any SuRfAce.

FOr pLaNar naNostRIps Shown in fig.\[fig:SkpHasE\]A ($Q=7.2$) anD FIg.\[Fig:SkphAse\]d ($Q=37.8$), the valueS oF $Q$ deviate wEaKly From inTEGer valueS due to the misalignment of THe momenTs lOcateD on tHe free bouNdaRies of The SAmple[@rOh13] and tHe theRmAl fLUCtuatIONs InhErEnt to the MoNTE CaRlo meThOd. FoR curved Surfaces, however, the SHapE distortion of The SkyrMIOnS anD ThEIr eVoLUtiON To stripe-like strUctures is nOt CHaRacterized BY inTeGer valuEs of $Q$, thUs the CAlculatIon of $Q$ basEd on Eq.(\[eq:tOpOl\_chARGe\]) aSsumes non-
 direction of the appliedfield. Thu s the  ab sol ut e va lueof $Q$ for a n a nost rip in the skyrmion ph ase e qu a ls t o  t he nu mber of sk y r mio ns  s upp or t ed .

Th e d ependen ce of thetop ol ogical charg e  o n the curv atu re angle issho wn inFi g.\ [ fig:Q \_v s\_a\ ] forn anostr ips withdi f ferent widths$ L _y $.
\[htb!\]![ (Color on l ine) Dependenc e of t he gr o u ndsta te topolog ic al ch a rge (Q) on t h e a n gle of curvat ure ($\phi_ 0 $)for ma gn eti c  nanos truct ur e s w ith differe nt s izes ($L\ timesL $) and$ a=2nm$.  The a ppl ied  mag n et ic  fi el d  is un ifo r m a long the  x -a xis ( circ l e s ,  clo sed  tri angle s) or radial(op en t r ian gles) . []{ data -l abel= "fig:Q _vs_a "} ](fig3.jpg "fig :"){ width="0. 95\ li new id th"}We not ice  th at $Q$remains alm os t c on stant up to an ang le $ \p hi_0\app rox 10 0 ^0 $a nd furth er  on  itd e creas es s m oo thly fro m thei ni ti al valu efor pl an arnan ostri p s (s kyrmio n phase) , toa  smaller value close to zero fo r na n otub es,  indicating  tha t  onl y as ma llf racti on of  t h ei nitial number of sk yr mionsare s tabilized (mi xed skyrmi o n - stripe p hase ) .I n contrast tothistrend, whe n  the app liedfield is  radial,t h e topolo gic alcha rge i sonly weakly d e p ende nt  on the  cu rvature  an gle . T his  w eak decay  of $Q$wi th  c ur vat ure a n gle is d ue  to  t hegradu a l redu ction  ofth ec o-p lanar c o nd i t ionbe tw eenthe  D MI ve ctor s  on  each l attice si te. A ra di al  fieldfor a cylindr ic al nanostr ip  is  geome t r ically a nalogous to the case of a unifo rmfield  nor mal to apla nar na nos t rip an d theconse rv ati o n  of t h e  t opo lo gical char g e  wi th cu rv atur e under lines the importan c e o f the normali tycond i t io n f o rt heap p lie d filed in the st abilizatio no fskyrmionso n a ny  surfac e.

For  plan a r nanos trips sho wn in Fig .\ [fig : S kph ase\]a ($Q =7.2$) a nd Fig.\[ f ig:Sk p ha se\]d  ($ Q=37.8 $) , t he va lues o f  $Q $ dev iate w ea kly fr om in te ger valu es due to the misalignm ent of  themom ents loca ted onthe freeboun daries ofthe  sa mple[ @ro h 13] a nd t h ethe r mal f luct u ations in h er ent t othe Monte C a r l o m ethod . F o r curv ed s urfaces, however, the shape dist orti o n  of  th e  sky rm ions and their  ev ol u t ion to s tr ipe-like st ructures  i s  notcharac terize d by in t e ge r  value s of  $Q $, thus t heca l culatio nof $Q$ ba sedon  Eq.(\ [eq:to p ol\_ c h arge\]) assumesnon-
 direction_of the_applied field. Thus the_absolute value_of_$Q$ for_a_nanostrip in the_skyrmion phase equals_to the number of_skyrmions supported.

The dependence_of_the topological charge on the curvature angle is shown in Fig.\[fig:Q\_vs\_a\] for nanostrips with_different_widths $L_y$.

\[htb!\]![_(Color_online)_Dependence of the ground state_topological charge (Q) on the_angle of_curvature ($\phi_0$) for magnetic nanostructures with different sizes_($L\times_L$) and $a=2nm$._The applied magnetic field is uniform along the x-axis_(circles, closed triangles) or radial (open_triangles). []{data-label="fig:Q_vs_a"}](fig3.jpg "fig:"){width="0.95\linewidth"}

We_notice_that_$Q$ remains almost constant_up to an angle $\phi_0\approx 100^0$_and further on it decreases smoothly_from the initial value for planar nanostrips_(skyrmion phase), to a smaller value_close to zero for nanotubes,_indicating that_only a small fraction of_the initial number_of skyrmions_are stabilized (mixed_skyrmion-stripe phase). In contrast to this_trend, when the_applied field is radial, the topological_charge_is only weakly_dependent_on_the curvature_angle. This weak_decay_of $Q$_with_curvature angle is due to the_gradual_reduction of the co-planar condition between the_DMI vectors on each_lattice_site. A radial field_for a cylindrical nanostrip is_geometrically analogous to the case of_a uniform_field normal_to a planar nanostrip and the conservation of the topological charge_with curvature underlines the importance of_the normality condition for_the applied_filed_in the stabilization_of_skyrmions on_any surface.

For planar nanostrips shown in Fig.\[fig:Skphase\]a_($Q=7.2$) and_Fig.\[fig:Skphase\]d ($Q=37.8$), the values of $Q$_deviate weakly from integer_values_due to the misalignment of the_moments located on the free boundaries_of the sample[@roh13] and the_thermal_fluctuations_inherent to the Monte Carlo_method. For curved surfaces, however, the_shape distortion of_the skyrmions and their evolution to stripe-like_structures_is not characterized by integer values_of_$Q$, thus the calculation of $Q$_based_on_Eq.(\[eq:topol\_charge\]) assumes non-