Arxiv:2102.13497V2 [Cond-Mat.Soft] 7 May 2021

Nematic versus ferromagnetic shells: new insights in curvature-induced eﬀects

Gaetano Napoli∗ Dipartimento di Matematica e Fisica ”E. De Giorgi”, Universit`adel Salento, Lecce (Italy).

Oleksandr V. Pylypovskyi† Helmholtz-Zentrum Dresden-Rossendorf e.V., Institute of Ion Beam Physics and Materials Research, 01328 Dresden, Germany and Kyiv Academic University, 03142 Kyiv, Ukraine

Denis D. Sheka‡ Taras Shevchenko National University of Kyiv, 01601 Kyiv, Ukraine

Luigi Vergori§ Dipartimento di Ingegneria, Universitàdi Perugia, Perugia (Italy) (Dated: May 10, 2021) Within the framework of continuum theory, we draw a parallel between ferromagnetic materials and nematic liquid crystals confined on curved surfaces, which are both characterized by local interaction and anchoring potentials. We show that the extrinsic curvature of the shell combined with the out-of-plane component of the director field gives rise to chirality effects. This interplay produces an effective energy term reminiscent of the chiral term in cholesteric liquid crystals, with the curvature tensor acting as a sort of anisotropic helicity. We discuss also how the different nature of the order parameter, a vector in ferromagnets and a tensor in nematics, yields different textures on surfaces with the same topology as the sphere. In particular, we show that the extrinsic curvature governs the ground state configuration on a nematic spherical shell, favouring two antipodal disclinations of charge +1 on small particles and four +1/2 disclinations of charge located at the vertices of a square inscribed in a great circle on larger particles.

Magnetic materials and liquid crystals are examples effects, for the review see Refs. [1,2,4]. Moreover, re- of materials with orientational order which give rise to cent advances in experimental techniques have also made textures whose complexity is as beautiful as challenging possible the manipulation of these effects for the design of to study. Their confinement to curved layers causes the new functional materials and applications for spintronics, emergence of geometry-induced effects that are not usu- shapeable magnetoelectronics, magnonics, biomedicine, ally observed in flat layers. Recent advances in the theory and soft robotics [1,2,4,5]. of curvilinear magnetism have highlighted a range of fas- Soft matter also provides an area in which the inter- cinating geometry-induced effects in the magnetic prop- play between geometry of the substrate and the order erties of materials [1,2]. When confined to thin curved parameter plays a crucial role. One example is provided domains, effective physical features arise from the inter- by liquid crystal (LC) shells [6]. These are microscopic play between the curved geometry and magnetic texture. colloidal particles coated with a thin layer of nematic According to continuum micromagnetic description LC, and have potential applications as the topological of the ferromagnetic media, the magnetic textures can defects (which may occur on them) can be engineered to be well described by the vector order parameter m = emulate the linear, trigonal and tetrahedral geometries M/ M , which is the normalized magnetization vector. of carbon atoms [7]. This feature opens up the possibil- The| energy| of the ferromagnet typically includes: (i) ity to design meso-atoms with special optical properties a short-range exchange energy that penalizes the non- whose valence and directional-binding can be controlled. uniformity of the magnetization; (ii) an anisotropic term Photonic lattices made of LC shells are the new frontier arXiv:2102.13497v2 [cond-mat.soft] 7 May 2021 that models the existence of directions preferred by the for the manufacture of a new-generation optical crypto- magnetization; (iii) an non-local term describing the graphic devices [8]. long-range magnetostatic interactions. Nematic LCs are aggregates of rodlike molecules. When ferromagnets are confined in thin curvilinear Within the classic theory of nematics[9], the average mi- layers, the magnetic energy can be decomposed [3] to croscopic molecular orientation is described by a sole vec- reveal the emergence of geometry-induced anisotropy tor order parameter n called the director. In the 1960s de and geometry-induced chiral interaction with emergent Gennes introduced the order-tensor theory which bases Dzyaloshinskii–Moriya interaction (DMI) as characteris- on the orientational probability distribution and provides tic example. With this decomposition, a number of new measures of the degree of orientation and biaxiality. In its effects in ferromagnetic spherical shells have been stud- simplest form, this theory uses as state variable a second- ied including topological patterning and magnetochiral order symmetric traceless tensor Q = s(n n 1 I), with ⊗ − 3 2 the scalar parameter s being the degree of orientation with a = 1, is in bijective correspondence with the real | | 2 that vanishes at points where there is no privileged direc- projective RP plane which, as known, is diffeomorphic tion. Contrarily to the director theory, de Gennes theory to the quotient space of S2 under the equivalence relation allows the study of nematic-isotropic phase transitions. induced by the antipodal map. We also note that using Most theories of nematic shells are based on energy N as a state variable for nematics is equivalent to fixing functionals defined on surfaces expressed in terms of vec- s = 1 in the Q tensor theory. − tor [10, 11] or tensor [12–14] order parameters. Theories We shall see that, given an energy density, taking the in which the director field n is purely tangential have the alignment tensor N as the state variable allows a richer flaw that topological defects, i.e. points where the direc- variety of equilibrium textures than on using the director is not uniquely defined, inevitably arise on surfaces tor n as a macroscopic descriptor of the orientation of with the same topology as sphere. Unavoidably, the en- nematic molecules. To justify this assertion, consider a ergy blows up in a neighbourhood of a defect. This flaw spherical shell. Due to the topology of the sphere, at can be overcome by introducing a theory in which the equilibrium, ferromagnets (whose orientation is described tensor order parameter Q is tangential, i.e. the normal by a unit vector) form vortices with integer topological to the surface is an eigenvector of Q corresponding to the charges, while nematics (whose orientation is more ap- zero eigenvalue. In this framework, melting points (points propriately described by a second-order tensor) may form at which the degree of orientation vanishes) occur, and also vortices with half-integer topological charge. The the energy is finite over the entire domain. A different difference in the possible topological charges of a vortex way to avoid singularities within the director theory is between the two theories stems from the different ori- to release the constraint that the director field is tangen- entability of the state variables: arrows for ferromagnets, tial (as suggested by recent experimental observations lines for nematics. [15, 16]), and to assign an energy cost to the out-of-plane We start by considering an ordinary 3D nematics con- components of n. This hypothesis, first proposed in [17], fined in a thin layer. In the bulk the energy due to the allows to deal with smooth fields as topological defects re- distortion of the molecular field is given by the celebrated vert into non-singular disclinations. This alternative will Frank formula[23], while at the boundary a surface en- be explored in this paper and will allow us to compare ergy term promotes a tangential alignment of the director the theories of ferromagnetic and nematic shells. Note field. As a result of a dimensionality reduction conducted that in a more comprehensive framework, fusion of the in section I in the supplementary information, such an en- defect core and out-of-plane escape can be considered si- ergy density approximates to the effective surface energy multaneously [18], as well as biaxial order reconstruction density [19, 20]. k 2 % In theories based on vector order parameters in which ws(N, sN) = sN + N (ν ν), (1) a privileged direction is defined everywhere on the shell, ∇ 4 |∇ | 2 · ⊗ E due to the different orientability of the state variables, ≡ ex there is a fundamental difference between the response of where k > 0 is the reduced| {z elastic} stiffness and % repre- ferromagnets and nematics. For nematics, the head-tail sents the anchoring strength. The sign of % determines symmetry of the molecules is translated at the macro- the direction of the easy axis. The easy axis is normal to scopic level by the assumption that n must be indistin- the shell if % is negative, tangential if % is positive. De- guishable from n. Roughly speaking, it is more appro- noting q (i = 1, .., 6) the components of N with respect priate to regard− the director as a line instead of an arrow i to the local basis B in the Voigt notation (see section [21, 22]. We may then conclude that no vector order pa- II in the supplementary information for details), we can rameter can describe properly the mechanical response decompose the exchange energy density of a nematics as of nematics. A tensor order parameter has to be introduced. For this reason, in our attempt to determining k 2 k k (γ) (γ) Eex = (ðγ qi) + Hijqiqj + D κγ L , (2) the equilibrium textures of a nematics on a curved shell 4 4 4 ij ij S, we choose to use the second-order alignment tensor E0 EA EDM N = n n instead of the director field n. Introduced the ≡ ≡ ≡ ⊗ local orthonormal basis B = e1, e2, ν , where e1 and where γ =| 1, 2,{zi, j} =| 1, ...,{z 6,} and| the{z symbol} ð de- { } e2 are the principal directions on S and ν = e1 e2 is notes the modified tangential derivative [3]. Except the unit normal to S, and expanding the director× field as for the different dimensions of the vector spaces which n = n1e1 + n2e2 + n3ν, the components of N form the q = (q1, ..., q6) and the unit magnetization vector belong symmetric 3 3-matrix (Nij = ninj)i,j=1,2,3. Observe to, the exchange energy density (2) is formally similar to that n and N×are topologically different from each other. the one for magnetic systems[24]. In fact, the set of unit vectors is in one-to-one correspon- This decomposition allows to separate the intrinsic ef- dence with the points of the two-dimensional unit sphere fects from the spurious effects due to the embedding of S2, while the set of second-order tensors in the form a a, a 2D surface into a 3D Euclidean space. In fact, the ⊗ 3 term w0 in the exchange energy density (2) involves only form the 3 3-matrix covariant derivatives and hence represents the intrinsic × part of the exchange interaction. This term is zero if and 0 κ2 0 only if the alignment tensor is uniform in the tangent [τ ] =  κ 0 0  . (4) plane. The other two terms, instead, account also for the − 1   extrinsic curvature of the shell.  0 0 0    Since H is a symmetric 6 6-matrix whose elements   are functions of the principal× curvatures only (see sec- From (4) it is easy to deduce that τ is skew-symmetric tion II in the supplementary information), the energy on spheres, symmetric on minimal surfaces. term EA in (2) couples the extrinsic curvature with the A direct inspection shows that the essential ingredient alignment and induces anisotropic effects. At each point, for the DMI, apart from curvature, is the existence of an EA reaches its absolute minimum when the texture is out-of-plane component of the texture. The DMI term purely tangential and the director field is aligned along favours rotation of the alignment around the normal di- the principal direction with the smallest modulus of cur- rection, where the sense of rotation, in both nematic and vature. To explain this point, observe that EA can be magnetic shells, is determined by the sign of the princi- 2 2 rewritten as EA = k[κ1(1 q2) + κ2(1 q1)]/2, with pal curvatures. This term is instead suppressed in purely 2− − q1, q2 T = (u, v) R : u, v [0, 1], u + v 1 , tangential or purely normal textures. The influence of and check∈ that,{ depending∈ on the magnitude∈ of the≤ mod-} surface geometry and easy axes on local interaction terms uli of the principal curvatures, it reaches its minimum in is sketched in Figure 2 in the paper by Sheka et al. [3]. T for q1 = 1 or q2 = 1. As a particular case, since any In spherical ferromagnetic shells, the interplay between unit vector tangent to a sphere is principal, at each point the DMI term and the normal anisotropy (% < 0) is re- on a spherical shell EA attains its absolute minimum at sponsible of the emergence of skyrmions states as per- any tangential alignment. In any case, the curvature- turbations of hedgehogs-like ground states [27]. Very re- induced anisotropy term is activated whenever the sur- cently, LCs skyrmions have been realized as micron sized face is curved. It combines with anchoring potential, solitons in a chiral nematic material confined between strengthening its effects if % > 0 and weakening them two parallel substrates [28]. if % < 0. Note that, purely tangential field are admissible on Observe that the intrinsic and anisotropic terms are surfaces of genus one, while are not allowed on surface present also when considering purely tangential order of genus zero. Thus, on a spherical shell, except for the parameters [11, 14], while the occurrence of the en- hedgehog configuration, all terms of the local interaction ergy term EDM is completely novel in the theory of ne- are unavoidably activated. matic shells. This term is responsible for a curvature- We now address the equilibrium textures on nematic induced chiral interaction. The Lifshitz invariants of N, spherical shells by the numerical minimization of the en- (γ) ergy functional (1) under appropriate constraints on the Lij = qiðγ qj qjðγ qi, combines with the principal cur- − (γ) variables q’s (see the Supplementary Information for de- vatures κγ , with Dij being skew-symmetric matrices tails). We have used a finite differences scheme with with constant entries, giving rise to an anisotropic DMI constant steps α = π/40 for colatitude and longitude, with geometry-dependent coefficients. A term similar to and ∆r = 2λ, with λ = k/(2%) being a characteristic EDM occurs in the Landau-de Gennes theory of cholesteric length which measures the scale over which the align- p LCs (which are intrinsically chiral materials) and is re- ment tensor varies on the shell at equilibrium, for the sponsible for the emergence of two-dimensional modu- radius r [0.4λ, 10λ]. The texture with the lowest en- lated states and the formation of vortices or skyrmion ergy results∈ from the interplay between the exchange and lattices [25, 26]. In the supplementary information we anchoring energies. For small enough values of the ratio prove that in terms of the components of the alignment r/λ the exchange energy dominates, which leads to ax- tensor EDM reads isymmetric textures that are almost uniform far from two small regions around two antipodal points where the ra- EDM = k εijαεαβLβγ Nihðγ Nhj, (3) dial component (q3) is pronounced (Figures1(a,b,c)). This configuration corresponds to the 3D onion texture where Lβγ are the components of the extrinsic curvature in ferromagnetic spherical shells[29, 30]. tensor. Comparing (3) with the chiral term in the energy As the ratio r/λ increases the role played by the an- density (1) in [25] we realize that the role played by the choring energy becomes more and more important, the (constant) helicity q0 in cholesterics is here substituted by nematics is forced to align tangentially to the shell, the second-order tensor τ = εαβLβγ eα eγ . This tensor and more complex textures occur. We find that for can be regarded as a sort of helicity tensor⊗ which accounts r/λ ξ? 6 the equilibrium texture with four discli- for the anisotropy induced by the extrinsic curvatures. nations≥ has≈ lower energy than the onion configuration The components of τ with respect to the local basis B (Figures1(a,d,e ,h,i)). In particular, we have found two 4

Onion Tetrahedral (a) Trigonal Squared Onion conﬁguration 0 1 Trigonal conﬁguration

(b) q1 (longitudinal) (c) q3 (radial) (f) q1 (longitudinal) (g) q3 (radial) 40

+1/2 +1 30 E/k

20 c3 Tetrahedral conﬁguration Squared conﬁguration

Energy c4 c2 (d) q1 (longitudinal) (e) q2 (colatitudinal) (h) q1 (longitudinal) (i) q3 (radial) 10

c1 +1/2 0 c3 c4 c2 c3 c4 c2 +1/2 0 2 4 6 8 10 +1/2 Sphere radius r/λ

FIG. 1. Equilibrium textures on nematic spherical shells. (a) Energies of stable onion, trigonal, tetrahedral and squared configurations as functions of the ratio r/λ. Insets show the configurations of the nematic director field. Red spots mark the positions of disclinations of the tetrahedral texture. Longitudinal (q1) and radial (q3) components of the onion configuration are represented by Kavrayskiy VII map projection in (b,c), where areas with different colours are separated by isolines. Longitudinal (q1) and colatitudinal (q2) components of the tetrahedral configuration are displayed in (d,e). Longitudinal and radial components of trigonal and squared configurations are sketched in (f,g,h,i). stable configurations with four +1/2 disclinations: one depend on the radius of the shell. in which the disclinations are located at the vertices of Textures with a single +2 disclination or with two an- a square inscribed in a great circle (Figures1(h,i)), and tipodal disclinations of charge +1/2 and +3/2 are config- the other with disclinations situated at the vertices of urations observed in cholesteric LCs shells [37]. Despite a tetrahedron (Figures1(d,e)). The former configura- the presence of the chiral term (3) in the energy, we found tion (which, for brevity, we shall henceforth refer to as them unstable for any r/λ [0.4, 10]. ∈ the squared configuration) has slightly lower energy than Some of these shell textures computed on a sphere the latter, which is instead stable for any r/λ [0.4, 10]. through our model approach agree with those obtained ∈ This result is in disagreement with other theoretical pre- using the Frank 2D theory [10], Landau-de Gennes 2D dictions in the literature on thin nematic shells [10, 31] theory [13, 38], Monte Carlo simulation studies [39, 40]. in which it is claimed that the tetrahedral configuration Numerical results show that, although the energy denis the ground state on spherical shells of any radius. A sity can be regarded as being the same in the two theo- possible experiment with particles of different sizes would ries with n and N, the class of solutions admitted within clarify this controversy. Instead, the metastability of the tensor theory is wider than that of vector theory. these solutions is in agreement with experimental ob- Thus, while ferromagnetics can admit only vortices of servations [6, 32] and Monte Carlo simulations [33]. It integer charge, nematics can also admit those of half- should be noted, that theoretical models with the purely integer charge. This aspect has already been addressed tangential alignment predict that the configuration with in the literature of LCs. Indeed, Ball and Zarnescu [21] two defects is unstable, in favour of the tetrahedral con- discuss the case of a nematic in a planar configuration, figuration [10, 34]. Stable configurations with two discli- confined in a finite plane region with two holes. Report- nations are instead reported for thick spherical shells or ing their words: “...it turns out that there are continu- spherical shells with non-uniform thickness [35, 36]. ous, even smooth, line fields for which it is impossible A fourth metastable solution we have tested numeri- to ‘assign arrows to the lines’; that is, one cannot make cally is the trigonal configuration with three disclinations a choice of a vector out of each line in such a way that located at the vertices of an isosceles triangle. One discli- no discontinuity is created”. It is clear that discontinu- nation has topological charge equal to +1, the other two ous configurations are not allowed as they require infinite equal to +1/2 (Figures1(f,g)). The distance between energy. the two +1/2 disclinations decreases as the ratio r/λ When trying to force a field of vectors or lines on a decreases. The trigonal configuration has been experi- surface S with the topology of the sphere, discontinuities mentally observed and theoretically predicted (within the are unavoidable. Topological defects arise and, according Frank theory) in [6]. Unlike our findings, the theory pro- to the Poincaré-Hopftheorem, the sum of the topolog- posed in [6] explains that the mutual distance between ical charges of all defects on S must be equal to 2, the the disclinations of the trigonal configuration does not Euler-Poincarécharacteristics of a sphere. The topolog- 5 ical charge of each defect is an integer multiple of 1/2. of active materials: From mechanobiology to smart de- However, it must be said that defects might not be the vices”, project n. 2017KL4EF3. only points at which the solution lacks continuity. In fact, while it is always possible to have a continuous field of lines on the whole sphere (of course, excluding the points where the defects are located), the same is not true for fields of vectors. In particular, in the presence of a de- ∗ [email protected] † [email protected] fect with fractional charge there are unavoidably curves ‡ [email protected] through which the vector field is discontinuous. An ex- § [email protected] ample can be found in the work of Vitelli and Nelson [10] [1] R. Streubel, P. Fischer, F. Kronast, V. P. Kravchuk, where such discontinuities are branch cuts in the com- D. D. Sheka, Y. Gaididei, O. G. Schmidt, and plex plane. So, assuming the points where the defects D. Makarov, Magnetism in curved geometries (topical re- are excluded, there would be more continuous configu- view), Journal of Physics D: Applied Physics 49, 363001 rations allowed with lines (nematics) than with arrows (2016). [2] E. Y. Vedmedenko, R. K. Kawakami, D. Sheka, P. Gam- (ferromagnets). When the field of lines or vectors can es- bardella, A. Kirilyuk, A. Hirohata, C. Binek, O. A. cape from the tangent plane, topological defects no longer Chubykalo-Fesenko, S. Sanvito, B. Kirby, J. Grollier, exist. However, since out-of-plane escape has a certain K. Everschor-Sitte, T. Kampfrath, C.-Y. You, and energy cost, the texture tends to lie on the surface when A. Berger, The 2020 magnetism roadmap, Journal of possible. For topological reasons this is not possible at all Physics D: Applied Physics 53, 453001 (2020). points on a sphere and, therefore, uplifts reverts to local- [3] D. Sheka, P. Pylypovskyi, O.V.An Landeros, and et al., ized vortices that can have different structures depending Nonlocal chiral symmetry breaking in curvilinear magnetic shells, Commun Phys 3, 128 (2020). on the theory. In other words, the texture around a lift- [4] P. Fischer, D. Sanz-Hernández, R. Streubel, and ing point mimics that of topological defects. A. Fernández-Pacheco, Launching a new dimension with In conclusion, we have introduced a novel theory for 3D magnetic nanostructures, APL Materials 8, 010701 nematic shells in which the optical axis is assumed not to (2020). be purely tangential. Accordingly, we have observed the [5] D. Makarov, M. Melzer, D. Karnaushenko, and O. G. occurrence of novel curvature-induced effects that have Schmidt, Shapeable magnetoelectronics, Applied Physics been observed within theories for ferromagnetic shells. Reviews 3, 011101 (2016). [6] T. Lopez-Leon, V. Koning, K. B. S. Devaiah, V. Vitelli, In particular, the extrinsic curvature can be seen as a and A. Fernandez-Nieves, Frustrated nematic order in source of chirality responsible for the onset of skyrmion- spherical geometries, Nature Physics 7, 391 (2011). or meron-type nonlinear waves. On the other hand, it [7] D. R. Nelson, Toward a tetravalent chemistry of colloids, must be said that the difference in the orientability of the Nano Letters 2, 1125 (2002). order parameters yields a class of equilibrium solutions [8] Y. Geng, J. Noh, I. Drevensek-Olenik, , R. Rupp, for nematics wider than that for ferromagnets. G. Lenzini, and J. P. F. Lagerwall, High-fidelity spher- We believe that the impact of this work goes beyond ical cholesteric liquid crystal bragg reflectors generating unclonable patterns for secure authentication, Scientific the liquid crystals. In particular, we expect that the Reports 6, 26840 (2016). theory can be extended to curvilinear antiferromagnets. [9] E. G. Virga, Variational Theories For Liquid Crystals These are be described by several vector order parame- Variational Theories For Liquid Crystals (Chapman- ters. In the simplest case, the σ-model allows antiferro- Hall, London, 1994). magnets to be described by the Néelvector[41]. Being [10] V. Vitelli and D. R. Nelson, Nematic textures in spherical a direction by definition, in many cases the Néelvec- shells, Phys. Rev. E 74, 021711 (2006). tor behaves like the nematic director and hence it allows [11] G. Napoli and L. Vergori, Extrinsic curvature effects on nematic shells, Phys. Rev. Lett. 108, 207803 (2012). textures with fractional topological charges [42–45]. The [12] P. Biscari and E. M. Terentjev, Nematic membranes: concept of curvilinear antiferromagnetism has been intro- Shape instabilities of closed achiral vesicles, Phys. Rev. duced in [46]. One can expect that curvilinear antiferro- E 73, 051706 (2006). magnets admit novel textures similar to those predicted [13] S. Kralj, R. Rosso, and E. G. Virga, Curvature control here for nematics. of valence on nematic shells, Soft Matter 7, 670 (2011). [14] G. Napoli and L. Vergori, Surface free energies for nematic shells, Physical Review E 85, 061701 (2012). [15] G. Durey, Y. Ishii, and T. Lopez-Leon, Temperature- ACKNOWLEDGMENTS driven anchoring transitions at liquid crystal/water in- terfaces, Langmuir 36, 9368 (2020). D. S. acknowledges the financial support from the [16] J. Noh, Y. Wang, H.-L. Liang, V. S. R. Jampani, A. Ma- Ministry of Education and Science of Ukraine (project jumdar, and J. P. F. Lagerwall, Dynamic tuning of the 19BF052-01). The work of G. N. and L. V. has been director field in liquid crystal shells using block copoly- funded by the MIUR (Italian Ministry of Education, Uni- mers, Phys. Rev. Research 2, 033160 (2020). versity and Research) project PRIN 2017, ”Mathematics 6

[17] Cladis, P.E. and Kléman,M., Non-singular disclinations [32] T. Lopez-Leon, A. Fernandez-Nieves, M. Nobili, and of strength s = + 1 in nematics, J. Phys. France 33, 591 C. Blanc, Nematic-smectic transition in spherical shells, (1972). Phys. Rev. Lett. 106, 247802 (2011). [18] A. L. Susser, S. Harkai, S. Kralj, and C. Rosenblatt, Tran- [33] M. A. Bates, Nematic ordering and defects on the surface sition from escaped to decomposed nematic defects, and of a sphere: A monte carlo simulation study, J. Chem. vice versa, Soft Matter 16, 4814 (2020). Phys. 128, 104707 (2008). [19] N. Schopohl and T. J. Sluckin, Defect core structure in [34] G. Napoli and L. Vergori, Cooling spherical nematic nematic liquid crystals, Phys. Rev. Lett. 59, 2582 (1987). shells, Submitted (2021). [20] G. Carbone, G. Lombardo, R. Barberi, I. Muˇseviˇc,and [35] A. Fernández-Nieves, V. Vitelli, A. S. Utada, D. R. Link, U. c. v. Tkalec, Mechanically induced biaxial transition in M. Márquez,D. R. Nelson, and D. A. Weitz, Novel de- a nanoconfined nematic liquid crystal with a topological fect structures in nematic liquid crystal shells, Phys. Rev. defect, Phys. Rev. Lett. 103, 167801 (2009). Lett. 99, 157801 (2007). [21] J. M. Ball and A. Zarnescu, Orientable and non- [36] V. Koning, T. Lopez-Leon, A. Fernandez-Nieves, and orientable line field models for uniaxial nematic liquid V. Vitelli, Bivalent defect configurations in inhomoge- crystals, Molecular Crystals and Liquid Crystals 495, neous nematic shells, Soft Matter 9, 4993 (2013). 221/[573] (2008). [37] A. Darmon, M. Benzaquen, S. Copar,ˇ O. Dauchot, and [22] J. M. Ball and A. Zarnescu, Orientability and energy T. Lopez-Leon, Topological defects in cholesteric liquid minimization in liquid crystal models, Archive for Ratio- crystal shells, Soft Matter 12, 9280 (2016). nal Mechanics and Analysis 202, 493 (2011). [38] I. Nitschke, M. Nestler, S. Praetorius, H. Löwen, and [23] F. C. Frank, On the theory of liquid crystals, Discuss. A. Voigt, Nematic liquid crystals on curved surfaces: Faraday Soc. 25, 19 (1958). a thin film limit, Proceedings of the Royal Society A: [24] D. D. Sheka, V. P. Kravchuk, and Y. Gaididei, Curva- Mathematical, Physical and Engineering Sciences 474, ture effects in statics and dynamics of low dimensional 20170686 (2018). magnets, Journal of Physics A: Mathematical and Theo- [39] H. Shin, M. J. Bowick, and X. Xing, Topological de- retical 48, 125202 (2015). fects in spherical nematics, Phys. Rev. Lett. 101, 037802 [25] A. Duzgun, J. V. Selinger, and A. Saxena, Comparing (2008). skyrmions and merons in chiral liquid crystals and mag- [40] M. A. Bates, G. Skaˇcej,and C. Zannoni, Defects and nets, Physical Review E 97, 062706 (2018). ordering in nematic coatings on uniaxial and biaxial col- [26] A. Duzgun, A. Saxena, and J. V. Selinger, Alignment- loids, Soft Matter 6, 655 (2010). induced reconfigurable walls for patterning and assem- [41] B. A. Ivanov, Mesoscopic antiferromagnets: statics, dy- bly of liquid crystal skyrmions, Phys. Rev. Research 3, namics, and quantum tunneling (Review), Low Temper- L012005 (2021). ature Physics 31, 635 (2005). [27] V. P. Kravchuk, U. K. Rößler,O. M. Volkov, D. D. Sheka, [42] I. E. Dzyaloshinski˘i, Domains and dislocations in antfer- J. van den Brink, D. Makarov, H. Fuchs, H. Fangohr, and romagnets, Pis’ma Zh. Exp. Teor. Fiz. 25, 110 (1977). Y. Gaididei, Topologically stable magnetization states on [43] B. A. Ivanov and V. E. Kireev, Spin disclination in a a spherical shell: Curvature-stabilized skyrmions, Phys. layered antiferromagnet with a screw dislocation, Journal Rev. B 94, 144402 (2016). of Experimental and Theoretical Physics Letters 73, 188 [28] A. Duzgun and C. Nisoli, Skyrmion spin ice in liquid (2001). crystals, Phys. Rev. Lett. 126, 047801 (2021). [44] V. E. Kireev and B. A. Ivanov, Localized magnetic non- [29] D. Kong, S. Wang, and C. Chen, Magnetization ground uniformities in an antiferromagnet with a system of dis- states and phase diagrams for a nanosized Co hollow locations, Low Temperature Physics 45, 1256 (2019). sphere: An onion-type magnetization state, Journal of [45] E. G. Galkina, V. E. Kireev, R. V. Ovcharov, R. S. Applied Physics 104, 013923 (2008). Khymyn, and B. A. Ivanov, Spin dynamics in antiferro- [30] M. I. Sloika, D. D. Sheka, V. P. Kravchuk, O. V. Pyly- magnets with domain walls and disclinations, Ukrainian povskyi, and Y. Gaididei, Geometry induced phase tran- Journal of Physics 65, 924 (2020). sitions in magnetic spherical shell, Journal of Magnetism [46] O. V. Pylypovskyi, D. Y. Kononenko, K. V. Yershov, and Magnetic Materials 443, 404 (2017). U. K. Rößler,A. V. Tomilo, J. Fassbender, J. van den [31] T.C. Lubensky and Jacques Prost, Orientational order Brink, D. Makarov, and D. D. Sheka, Curvilinear one- and vesicle shape, J. Phys. II France 2, 371 (1992). dimensional antiferromagnets, Nano Letters 20, 8157 (2020). Supplementary Information: Nematic versus ferromagnetic shells: new insights in curvature-induced effects

1, 2, 3, 4, 5, Gaetano Napoli, ∗ Oleksandr V. Pylypovskyi, † Denis D. Sheka, ‡ and Luigi Vergori § 1Dipartimento di Matematica e Fisica ”E. De Giorgi”, Universitàdel Salento, Lecce (Italy). 2Helmholtz-Zentrum Dresden-Rossendorf e.V., Institute of Ion Beam Physics and Materials Research, 01328 Dresden, Germany 3Kyiv Academic University, 03142 Kyiv, Ukraine 4Taras Shevchenko National University of Kyiv, 01601 Kyiv, Ukraine 5Dipartimento di Ingegneria, Universitàdi Perugia, Perugia (Italy) (Dated: May 10, 2021) In this supplementary material we derive the surface Frank-Rapini-Papoular energy density as a result of a perturbation analysis of the usual 3D energy functional under the assumption that the region occupied by the nematics is very thin. We then decompose the energy density as the sum of the four contributions that are discussed in details in the Letter. Finally, by means of the numerical minimization of the energy density we find equilibrium textures on a spherical nematic shell with two, three and four disclinations.

κ PARAMETRIZATION OF THE REGION e = e Ω + 2s ν e , (S3c) ∇ 2 1 ⊗ 1 ξκ ⊗ 2 OCCUPIED BY THE NEMATIC SHELL. − 2s DIFFERENTIAL OPERATORS where Ω = (Ω e + Ω e ), − 1 1 2 2 Let us assume that the nematic shell occupies a thin Ω1s ξ sκ1s e2s region V of thickness h around a regular compact surface Ω1 = + ∇ · , − 1 ξκ2s (1 ξκ1s)(1 ξκ2s) S. Let νs be the normal unit vector field to S. We − − − (S3d) parametrize points in the bulk through a coordinate set Ω ξ κ e Ω = 2s ∇s 2s · 1s , (u, v, ξ) such that 2 − 1 ξκ − (1 ξκ )(1 ξκ ) − 1s − 1s − 2s p(u, v, ξ) = p (u, v) + ξν (u, v), (S1) S s and Ω1s and Ω2s are the geodesic curvatures of the lines of curvature on S. where p is the normal projection of p onto S, and ξ , S The differential operator in (S3d) is the surface gra- with ξ [ h/2, h/2], is the distance of p from the same| | s dient. To introduce this differential∇ operator in the most surface.∈ Such− a coordinate set is well defined in a finite general setting, let Φ be a smooth scalar, vector or tensor neighbourhood of S. Next, we introduce the principal field defined on S. Then the surface gradient of Φ is the curvatures κ (p ) and κ (p ) of S at p S, and as- 1s S 2s S S tensor field resulting from the composition of the usual sume that ∈ gradient of Φ and the projection onto the tangent plane 1 of S, P = I νs νs, namely h min (max κ1s(pS) , κ2s(pS) )− = `. (S2) − ⊗ pS S {| | | |} ∈ sΦ = ( Φ)P. (S4) For every fixed ξ [ h/2, h/2], equation (S1) defines ∇ ∇ ∈ − a parallel surface Sξ = pS + ξνs(pS): pS S located In analogy with the usual gradient operator, for a given at distance ξ from S with{ the unit normal∈ vector} field | | vector field u we define the surface divergence and the ν : p Sξ νs(pS). In this way, the unit vector field surface curl of u as the trace of su and twice the axial ν is defined∈ 7→ on the entire region V . The second-order ∇ vector corresponding to the skew-symmetric part of su, tensor ν is symmetric. Its eigenvectors are ν (with a respectively. In formulas we have ∇ null eigenvalue)∇ and the unit vector fields arXiv:2102.13497v2 [cond-mat.soft] 7 May 2021 divsu = tr su = su P, curlsu = ε su, ei : p = pS + ξνs V eis(pS)(i = 1, 2), ∇ ∇ · − ∇ ∈ 7→ where ε denotes the Ricci alternator. with e1s and e2s denoting the principal directions fields on S. The spatial gradients of the eigenvectors ν, e1 and e2 are[1]: DERIVATION OF THE SURFACE κ1s κ2s FRANK-RAPINI-PAPOULAR ENERGY ν = e e e e , (S3a) ∇ −1 ξκ 1 ⊗ 1 − 1 ξκ 2 ⊗ 2 − 1s − 2s Let n be a unit vector, called the director, which rep- κ resents the average alignment of the molecules of the ne- e = e Ω + 1s ν e , (S3b) ∇ 1 − 2 ⊗ 1 ξκ ⊗ 1 matics. The celebrated Frank formula[2] for the elastic − 1s 2 energy density (per unit of volume) associated with the where, due to (S9) and following the similar arguments director distortion consists of four terms as in Appendix A in [1],

2 2 2 2wOZF = K1(divn) + K2(n curln) + K3 n curln sθs e1s sθs e2s · | × | θ = ∇ · e1 + ∇ · e2, (S11) + (K + K )[( n) ( n)T (divn)2] (S5) ∇ 1 ξκ1s 1 ξκ2s 2 24 ∇ · ∇ − − − and where the constants K1, K2, K3, and K24 are called the splay, twist, bend, and saddle-splay moduli, respectively. sϕs e1s sϕs e2s To ensure a stable undistorted configuration of a nematic ϕ = ∇ · e1 + ∇ · e2. (S12) ∇ 1 ξκ1s 1 ξκ2s liquid crystal in the absence of external fields or confine- − − ments, the three moduli Ki(i = 1, 2, 3) must be non- The scalar fields θs and ϕs on the right-hand sides of negative, whereas the elastic saddle-splay constant must equations (S11) and (S12) are the restrictions of θ and ϕ obey Ericksen’s inequalities[3]: to the surface S, namely θ : p S θ(p ) and ϕ : s S ∈ 7→ S s pS S ϕ(pS). From (S3a), (S3d) and (S10)–(S12) K K ,K + K 2K . ∈ 7→ | 24| ≤ 2 2 24 ≤ 1 it is evident that although the director field is assumed to be constant along directions normal to S, its spatial In the attempt to capture also the effects of the bound- gradient does depend on the spatial variable ξ in the aries ∂S h/2 on the orientation of the nematic director, ± normal direction ν. we consider the energy functional resulting from the sum At ξ = 0 the parametrization of n (S8) and (S9) give of the Frank energy and the Rapini-Papoular surface an- the restriction n of the director field to the surface S, choring energy[4] s

ns = sin θs cos ϕse1s + sin θs sin ϕse2s + cos θsνs, (S13) W = w (n, n)dV OZF ∇ ZV while evaluating (S3a)-(S3d), (S11) and (S12) at ξ = 0 % + (n υ)2da, (S6) gives the surface gradient of ns, 4 Sh/2 S h/2 · Z ∪ − n = (ν n ) ( ϕ Ω ) + ν Ln ∇s s s × s ⊗ ∇s s − s s ⊗ s where the constant % represents the anchoring strength νs (νs ns)ns and υ is the unit outward normal to S . The anchor- (νs ns)L − · sθs, (S14) h/2 − · − νs ns ⊗ ∇ ing energy in (S6) penalizes the deviation± of the molecules | × | where at the boundaries S h/2 from tangential directions if % is positive, from the normal± direction if % is negative. Ω = Ω e Ω e , (S15a) Under the assumption (S2) on the smallness of the s − 1s 1s − 2s 2s thickness of the region V occupied by the nematics, the energy functional (S6) approximates to a surface inte- L = κ1se1s e1s + κ2se2s e2s, (S15b) gral over S, with the integrand representing the surface ⊗ ⊗ Frank-Rapini-Papular energy density. To prove this, we are the vector parametrizing the spin connection on S, first assume that the director field n does not vary along and the extrinsic curvature tensor of S, respectively. In- directions normal to S, that is it satisfies the property cidentally, we recall that the mean and Gaussian curvatures of S are defined as n(p) = n(p ) (S7) S 1 1 H = trL and K = (trL)2 tr(L2) , (S16) at each point p = p + ξν V . We then parametrize 2 2 − S s ∈ the director field as respectively. Recalling the definition of the parameter ε (c.f. (S2)), n = sin θ cos ϕe1 + sin θ sin ϕe2 + cos θν, (S8) from (S10)–(S14) we deduce that (see Napoli and Vergori where, in view of ansatz (S7), the angles θ and ϕ satisfy [1] for details) the properties n = n + o() as ε 0, (S17) ∇ ∇s s → θ(p) = θ(pS) and ϕ(p) = ϕ(pS) (S9) whence for all p = p + ξν V . S s ∈ From (S3a)-(S3c) and (S8) the gradient of n is found divn = divsns + o(), curln = curlsns + o() (S18) to be as ε 0. → n = (ν n) ( ϕ Ω) ν ( ν)n With the parametrization (S1) the volume element dV ∇ × ⊗ ∇ − − ⊗ ∇ ν (ν n)n reads +(ν n) ν − · θ, (S10) · ∇ − ν n ⊗ ∇ dV = (1 2Hξ + Kξ2)dAdξ (S19) | × | − 3 where dA is the area element on S. On the other hand, respectively. In view of (S2) the volume and area ele- on the boundaries S h/2 the unit outward normal and ments(S19) and (S20) approximate to the area element are ±υ = ν and ± s dV = dAdξ + o(ε) and dA h/2 = dA + o(ε) (S21) ± as ε 0. Finally,→ from (S7), (S13), (S17), (S18), (S20) and (S21) h2 dA h/2 = 1 hH + K dA, (S20) we deduce that ± ∓ 4

1 W = K (divn)2 + K (n curln)2 + K n curln 2 + (K + K ) ( n) ( n)T (divn)2 dV 2 1 2 · 3| × | 2 24 ∇ · ∇ − ZV % 2 % 2 + (n νs) dAh/2 + [n ( νs)] dA h/2 − Sh/2 4 · S h/2 4 · − Z Z − h/2 1 2 2 2 2 = K1(divn) + K2(n curln) + K3 n curln (1 2Hξ + Kξ )dA dξ 2 h/2 S · | × | − Z− ( Z ) h/2 1 T 2 2 + (K2 + K24) ( n) ( n) (divn) (1 2Hξ + Kξ )dA dξ 2 h/2 S ∇ · ∇ − − Z− ( Z ) % h2 % h2 + (n ν )2 1 hH + K dA + (n ν )2 1 + hH + K dA 4 · s − 4 4 · s 4 ZS ZS h/2 1 2 2 2 = K1(divsns) + K2(ns curlsns) + K3 ns curlsns dA dξ 2 h/2 S · | × | Z− ( Z ) h/2 1 T 2 % 2 + (K2 + K24) ( sns) ( sns) (divsns) dA dξ + (ns νs) dA + o(ε) 2 h/2 S ∇ · ∇ − S 2 · Z− ( Z ) Z h = K (div n )2 + K (n curl n )2 + K n curl n 2 dA 2 1 s s 2 s · s s 3| s × s s| ZS 1 + h(K + K )[( n ) ( n )T (div n )2] + %(n ν )2 dA + o(ε) as ε 0. (S22) 2 2 24 ∇s s · ∇s s − s s s · s → ZS

SURFACE ENERGY DENSITY IN TERMS OF We have then proven that W Ws = wsdA, with ≈ THE ALIGNMENT TENSOR ZS

2w = k (div n )2 + k (n curl n )2 + k n curl n 2 In what follows we shall consider the surface energy s 1 s s 2 s · s s 3| s × s s| T 2 density (S24), and, for simplicity of notation, we shall + (k2 + k24)[( sns) ( sns) (divsns) ] ∇ · ∇ − omit the subscript s to n, e1, e2 and ν. + %(n ν )2, (S23) s · s We now show that the moduli of the surface gradients of the director n and alignment tensor N = n n are ⊗ where ki = hKi (i = 1, 2, 3, 24) are, respectively, the re- related each other. scaled splay, twist, bend and saddle-splay constants with From the deﬁnition of surface gradient (S4), the com- the same physical dimensions as energy per unit area. ponents of the surface gradient of the alignment tensor Within the one-constant approximation (K1 = K2 = are found to be K3 = K and K24 = 0 in (S5)) the surface energy density (S23) reduces to ∂Nij ∂ni ∂nj ( sN)ijh = Plh = Plhnj + ni Plh ∇ ∂xl ∂xl ∂xl k % n n + n n . (S25) w = n 2 + (n ν )2, (S24) ≡ i;h j i j;h s 2 |∇s s| 2 s · s (In (S25), and in the subsequent formulas, the Einstein 2 with k = hK. summation convention is used.) Since ni = 1 (the direc- 4 tor is a unit vector), we deduce that ðγ N3β = γ N3β + εβζ N3ζ Ωγ , (S32b) ∇ 2 nini;h = (ni );h nini;h = nini;h = nini;h = 0, − − ⇒ ðγ N33 = γ N33. (S32c) ∇ which combined with (S25) yields As known, the space of symmetric second-order ten- 2 2 2 2 2 2 sors is isomorphic to the six-dimensional Euclidean vec- sN = N = n n + 2nini;hnjnj;h + n n |∇ | ijh i;h j i j;h tor space. Speciﬁcally, on setting 2 2 ∂ni 2 = 2ni;h = 2 Plh = 2 sn . (S26) ∂xl |∇ | N11 = q1,N22 = q2,N33 = q3, (S33)

In view of this identity, the surface energy density (S24) √2N23 = q4, √2N13 = q5, √2N12 = q6, can be rewritten in terms of the alignment tensor as we can identify N with the vector q = (q1, q2, q3, q4, q5, q6). Relations (S33) are known as k 2 % ws = sN + N (νs νs), (S27) the Voigt notation. On using the Voigt notation it 4 |∇ | 2 · ⊗ can be proven that the space of second-order tensors and the six-dimensional Euclidean vector space are also On the other hand, with respect to the local basis isometric, whence N 2 = q 2, N 2 = q 2 and e , e , e = ν on S, the surface gradient of the align- s s 1 2 3 ( N )2 = ( q )2.| | | | |∇ | |∇ | {ment tensor reads} ðγ ij ðγ k From (S31), (S32) and (S33) we obtain (after lengthy algebra) N = [ N (e e e )] (e e e ), (S28) ∇s ∇s · i ⊗ j ⊗ γ i ⊗ j ⊗ γ N 2 = [ N (e e )]2 (S34) where, Latin indices range from 1 to 3, Greek ones from |∇s | ∇γ · i ⊗ j 2 (γ) (γ) 1 to 2 ( sN (ei ej ν) = 0 for all i, j = 1, 2, 3), and = (ðγ qi) + Hijqiqj + D κγ L , ∇ · ⊗ ⊗ ij ij N (e e e ) = N (e e ) (S29) where ∇s · i ⊗ j ⊗ γ ∇γ · i ⊗ j = γ Nij + Nhj( γ eh) ei + Nik( γ ek) ej. ðγ q1 = γ q1 + √2Ωγ q6, (S35a) ∇ ∇ · ∇ · ∇ In view of the symmetry of N we have that ðγ q2 = γ q2 √2Ωγ q6, (S35b) N (e e ) = N (e e ), (S30) ∇ − ∇γ · i ⊗ j ∇γ · j ⊗ i whence, with the aid of (S3) evaluated at ξ = 0, we ðγ q3 = γ q3, (S35c) ∇ obtain

ðγ q4 = γ q4 Ωγ q5, (S35d) γ N (eα eβ) = ðγ Nαβ (N3βδαγ +Nα3δβγ )κγ , (S31a) ∇ · ⊗ − ∇ −

ðγ q5 = γ q5 + Ωγ q4, (S35e) γ N (ν eβ) = ðγ N3β + κγ (Nγβ N33δβγ ), (S31b) ∇ ∇ · ⊗ −

ðγ q6 = γ q6 √2Ωγ (q1 q2), (S35f) γ N (ν ν) = ðγ N33 + 2Nγ3κγ , (S31c) ∇ − − ∇ · ⊗ (γ) where L = (qiðγ qj qjðγ qi), (S36) ij −

ðγ Nαβ = γ Nαβ + (εαζ Nζβ + εβζ Nαζ )Ωγ , (S32a) H = (Hij)i,j=1,...,6 is the symmetric sixth order matrix ∇ 5

2κ2 0 2κ2 0 0 0 1 − 1  0 2κ2 2κ2 0 0 0  2 − 2  2 2 2 2   2κ1 2κ2 2(κ1 + κ2) 0 0 0   − −  H =   ,, (S37)  0 0 0 κ2 + 4κ2 0 0   1 2     2 2   0 0 0 0 4κ1 + κ2 0       0 0 0 0 0 κ2 + κ2   1 2   

(γ) (γ) and D = (Dij )i,j=1,...,6 (γ = 1, 2) are skew-symmetric SURFACE ENERGY DENSITY ON A matrices the independent non-zero entries of which are SPHERICAL SHELL

D (1) = D (1) = √2, D (1) = 1, (S38a) 15 − 35 46 − On a spherical shell of radius r the local basis formed and by the principal directions and the unit normal is eϕ, eθ, er , with er, eθ and eϕ being the radial, co- D (2) = D (2) = √2, D (2) = 1. (S38b) { − } 24 − 34 56 − latitudinal and longitudinal directions, respectively, the Finally, in view of (S26), on using the Voigt notation principal curvatures are κϕ = κθ = 1/r and the vector the surface energy density (S24) becomes parametrizing the spin connection is Ω = (cot θ/r)eϕ. Therefore, for a spherical shell the energy density terms k 2 k k (γ) (γ) % E0, EA and EDM, and the whole energy density (S39) ws = (ðγ qi) + Hijqiqj + D κγ L + q3. 4 4 4 ij ij 2 read, respectively,

=E0 =EA =EDM (S39) | {z } | {z } | {z } k 2 To find stable equilibrium configurations one has to √ 2 E0 = 2 q1,ϕ csc θ + 2q6 cot θ + q1,θ minimize the energy functional Ws = wsdA subject 4r ( S to some constraints on the scalar fields q’s. In fact, since 2 R + q csc θ √2q cot θ + q2 q1, ..., q6 are the components of the alignment tensor N = 2,ϕ − 6 2,θ n n they must satisfy the following inequalities and + q2 csc2 θ + q2 + (q csc θ q cot θ)2 equations⊗ 3,ϕ 3,θ 4,ϕ − 5 2 2 2 + q4,θ + (q5,ϕ csc θ + q4 cot θ) + q5,θ qi 0 (i = 1, 2, 3), (S40a) ≥ 2 2 + q6,ϕ csc θ √2(q1 q2) cot θ + q6,θ , q + q + q = 1, (S40b) − − ) 1 2 3 h i

q2 q2 q2 q q 4 = q q 5 = q q 6 = 0, (S40c) 2 3 − 2 1 3 − 2 1 2 − 2 k 2 2 2 2 2 2 EA = 2q1 + 2q2 + 4q3 + 5q4 + 5q5 + 2q6 Inequalities (S40a) impose the non-negativeness of q1, q2 4r2 and q3 as these components of N are the squares of the 4q1q3 4q2q3 , components of the nematic director n. Since n is a unit − − vector, the trace of N must be equal to unity, which justiﬁes (S40b). Finally, equations (S40c) together with (S40a) and (S40b) imply that the components of N form k a matrix with rank equal to one as it has to be because √ EDM = 2 2(q1q5,ϕ q5q1,ϕ q3q5,ϕ + q5q3,ϕ) the alignment tensor is of the form n n. 2r ( − − We conclude this section by observing⊗ that by using h q4q6,ϕ + q6q4,ϕ csc θ (S31) the energy density term EDM in (S39) can be writ- − ten also as i + √2q (2q q q ) 3q q cot θ 4 1 − 2 − 3 − 5 6 E = kε ε L N N . (S41) DM ijα αβ βγ ikðγ kj + √h 2(q q + q q q q iq q ) 2 4,θ 4 3,θ − 4 2,θ − 3 4,θ This expression allows a comparison between the Dzyaloshinskii–Moriya interaction term E and the chi- q5q6,θ + q6q5,θ , DM − ral term in the energy density (1) in [5] (see the Letter). ) 6 and SIMULATIONS

To determine the equilibrium configuration of nematic molecules, we numerically minimize the energy functional Ws with energy density as in (S42). To do this, the con- tinual fields ql, l = 1, 6 are replaced by their discrete counterparts on a finite difference mesh (FDM) with constant steps along the longitudinal and colatitudinal di- 2 k 2 rections ∆φ = 2π/M and ∆θ = π/N, respectively, where ws = (q1,θ) + q1,ϕ csc θ √2q5 +√2q6 cot θ) 4r2 − − N > 1 and M > 1 are integer numbers determining h 2 discretization of the mesh. If f ij represents the FDM + q6,ϕ csc θ q4 √2q1 cot θ + √2q2 cot θ − − counterpart of the function f(θ, φ) on a sphere, then 2 + q csc θ + √2q √2q + q cot θ 5,ϕ 1 − 3 4 f(i∆θ, j∆φ) = f ij, 2 2 + q √2q + q csc θ √2q cot θ f(θ = 0) = f 0, f(θ = π) = f N , (S43) 2,θ − 4 2,ϕ − 6 2 2 i = 1,N 1, j = 0,M 1. + q + √2q √2q + q + √2q − − 4,θ 2 − 3 3,θ 4 2 Integration of Ws is performed using the double Riemann 2 + q3,ϕ csc θ + √2q5 + (q4,ϕ csc θ q5 cot θ + q6) sum method. We minimize the function W = W +W − tot s p % using the conjugate gradients method with the penalty + (q + q )2 + (q q )2 + q . (S42) 5,θ 6 6,θ − 5 2 3 function i

2 2 2 2 2 2 W = Λ (1 qij qij qij)2 + 2qijqij qij + 2qijqij qij + 2qijqij qij p − 1 − 2 − 3 1 2 − 6 2 3 − 4 3 1 − 5 ij ( ) X 2 2 2 2 2 2 + Λ (1 q0 q0 q0)2 + 2q0q0 q0 + 2q0q0 q0 + 2q0q0 q0 (S44) − 1 − 2 − 3 1 2 − 6 2 3 − 4 3 1 − 5 h i h i h i 22 2 2 2 2 +(1 q0 q0 q0)2 + 2qN qN qN + 2qN qN qN + 2qN qN qN , − 1 − 2 − 3 1 2 − 6 2 3 − 4 3 1 − 5 h i h i h i

3 where Λ = 10 . To avoid saddle points of Wtot, the The relaxed state is shown in Fig. S2. multiplier q3 is replaced by 1 q1 q2 in the anchoring The initial state with four +1/2 disclinations at the term. − − vertices of a square inscribed in a great circle is given by To test stability of different configurations, the strictly in-surface initial state is given 1 7 σ1/2 4 = arg sin θ (cos φ + i cos θ sin φ) . (S48) × 2 2 q = cos2 σ, q = sin2 σ, q = √2 sin σ cos σ, "r # 1 2 6 (S45) q3 = q4 = q5 = 0 The relaxed state is shown in Fig. S3. The initial state with with three disclinations, one with for the sphere of each radius with the concrete configu- topological charge equal to +1 and the other two equal ration determined by the angle σ. For the onion config- to +1/2 is in the form uration it reads 1 2iφ 4 θ 1 2 θ σonion = π/2. (S46) σ1+1/2 2 = arg 2ie− sin sin . (S49) × 2 2 − 2 2 The relaxed state is shown in Fig. S1. Note, that in gen- The exchange and anchoring energies, and penalty W eral case, the axis of symmetry not necessary to coincide p of the stable equilibrium configurations are shown in with the direction along the sphere poles. Fig. S5. The stability of the nematic configurations has The initial state with four +1/2 disclinations at the been tested with step ∆(r/λ) = 0.4 for r/λ [0.4, 4], vertices of a tetrahedron is given by with step ∆(r/λ) = 2 for r/λ [4, 10]. ∈ ∈ 1 iφ 3 θ θ 2iφ 4 θ We conclude by observing that for r/λ . 2 the ex- σ1/2 4 = arg 4ie cos sin √2e− sin . × 2 2 2 − 2 change energy of the onion configuration is small and (S47) tends to zero as r/λ approaches zero. This means that 7 for small values of r/λ the director field tends to be uni- [2] F. C. Frank, On the theory of liquid crystals, Discuss. form in space. Faraday Soc. 25, 19 (1958). [3] J. L. Ericksen, Inequalities in liquid crystal theory, Phys. Fluids 9, 1205 (1966). [4] A. Rapini and M. Paapoular, Distortion d’une lamelle nématiquesous champ magnétique.conditions d’angrage aux paroix, Le Journal de Physique Colloques 30, C4 ∗ [email protected] (1969). † [email protected] [5] A. Duzgun, J. V. Selinger, and A. Saxena, Comparing ‡ [email protected] skyrmions and merons in chiral liquid crystals and mag- § [email protected] [1] G. Napoli and L. Vergori, Surface free energies for nematic nets, Physical Review E 97, 062706 (2018). shells, Physical Review E 85, 061701 (2012). 8

FIG. S1. Onion conﬁguration of the nematic director on the sphere of radius r = 4λ (λ = k/(2%). (a) Flux lines of the director ﬁeld on the sphere. (b–d) Components of alignmenttensor tensor N. p

FIG. S2. Tetrahedral conﬁguration on a sphere of radius r = 10λ. The four +1/2 disclinations are located at vertices of tetrahedron (a–d). The components of the alignment tensor N are sketched in (e–g). 9

FIG. S3. Squared conﬁguration on a sphere of radius r = 10λ. The four +1/2 disclinations located at the vertices of a square inscribed in a great circle (a–d). The components of the alignment tensor N are sketched in (e–g).

FIG. S4. Trigonal conﬁguration on a sphere of radius r = 8λ. One +1 disclination and two +1/2 disclinations are located at the vertices of an isosceles triangle inscribed in a great circle (a,b,c). (d) Topologically trivial conﬁguration at the south pole. (e–g) The components of the alignment tensor N are sketched in (e,f,g). 10

Onion Tetrahedral Trigonal Squared

(a) (b) 10 30

5 20

Anchoring energy 0 10 Exchange energy 0.15 (c) 0.1 0.05 0 Penalty 0 0 2 4 6 8 10 0 2 4 6 8 10 Sphere radius r/λ Sphere radius r/λ

FIG. S5. Energy contributions for the diﬀerent equilibrium conﬁgurations. (a) Exchange energy. (b) Anchoring energy. (c) Penalty Wp to satisfy constraints (S40).