Effective Topological Charge Cancelation Mechanism

Biaxial Structures of Localized Deformations and Line-like Distortions in Effectively 2D Nematic Films

We numerically studied localized elastic distortions in curved, effectively two-dimensional nematic shells. We used a mesoscopic Landau-de Gennes-type approach, in which the orientational order is theoretically considered by introducing the appropriate tensor nematic order parameter, while the three-dimensional shell shape is described by the curvature tensor. We limited our theoretical consideration to axially symmetric shapes of nematic shells. It was shown that in the surface regions of stomatocyte-class nematic shell shapes with large enough magnitudes of extrinsic (deviatoric) curvature, the direction of the in-plane orientational ordering can be mutually perpendicular above and below the narrow neck region. We demonstrate that such line-like nematic distortion configurations may run along the parallels (i.e., along the circular lines of constant latitude) located in the narrow neck regions of stomatocyte-like nematic shells. It was shown that nematic distortions are enabled by the order reconstruction mechanism. We propose that the regions of nematic shells that are strongly elastically deformed, i.e., topological defects and line-like distortions, may attract appropriately surface-decorated nanoparticles (NPs), which could potentially be useful for the controlled assembly of NPs.

Coupling of nematic in-plane orientational ordering and equilibrium shapes of closed flexible nematic shells

The impact of the intrinsic curvature of in-plane orientationally ordered curved flexible nematic molecules attached to closed 3D flexible shells was studied numerically. A Helfrich-Landau-de Gennes-type mesoscopic approach was adopted where the flexible shell's curvature field and in-plane nematic field are coupled and concomitantly determined in the process of free energy minimisation. We demonstrate that this coupling has the potential to generate a rich diversity of qualitatively new shapes of closed 3D nematic shells and the corresponding specific in-plane orientational ordering textures, which strongly depend on the shell's volume-to-surface area ratio, so far not predicted in mesoscopic-type numerical studies of 3D shapes of closed flexible nematic shells.

Stable Assemblies of Topological Defects in Nematic Orientational Order

We considered general mechanisms enabling the stabilization of localized assemblies of topological defects (TDs). There is growing evidence that physical fields represent fundamental natural entities, and therefore these features are of interest to all branches of physics. In general, cores of TDs are energetically costly, and consequently, assemblies of TDs are unfavorable. Owing to the richness of universalities in the physics of TDs, it is of interest to identify systems where they are easily experimentally accessible, enabling detailed and well-controlled analysis of their universal behavior, and cross-fertilizing knowledge in different areas of physics. In this respect, thermotropic nematic liquid crystals (NLCs) represent an ideal experiment testbed for such studies. In addition, TDs in NLCs could be exploited in several applications. We present examples that emphasize the importance of curvature imposed on the phase component of the relevant order parameter field. In NLCs, it is represented by the nematic tensor order parameter. Using a simple Landau-type approach, we show how the coupling between chirality and saddle splay elasticity, which can be expressed as a Gaussian curvature contribution, can stabilize Meron TDs. The latter have numerous analogs in other branches of physics. TDs in 2D curved manifolds reveal that the Gaussian curvature dominantly impacts the assembling and stabilization of TDs. Furthermore, a strong enough curvature that serves as an attractor for TDs is a respective field that could be imposed in a fast enough phase transition. Assemblies of created TDs created in such a disordered environment could be stabilized by appropriate impurities.

History-dependent phase transition character

We consider history-dependent behavior in domain-type configurations in orientational order that are formed in configurations reached via continuous symmetry-breaking phase transitions. In equilibrium, these systems exhibit in absence of impurities a spatially homogeneous order. We focus on cases where domains are formed via (i) Kibble-Zurek mechanism in fast enough quenches or by (ii) Kibble mechanism in strongly supercooled phases. In both cases, domains could be arrested due to pinned topological defects that are formed at domain walls. In systems exhibiting polar or quadrupolar order, point and line defects (disclinations) dominate, respectively. In particular, the disclinations could form complex entangled structures and are more efficient in stabilizing domains. Domain patterns formed by fast quenches could be arrested by impurities imposing a strong enough random-field type disorder, as suggested by the Imry-Ma theorem. On the other hand, domains formed in supercooled systems could be also formed if large enough energy barriers arresting domains are established due to large enough systems' stiffness. The resulting effective interactions in established domain-type patterns could be described by random matrices. The resulting eigenvectors reveal expected structural excitations formed in such structures. The most important role is commonly played by the random matrix largest eigenvector. Qualitatively different behavior is expected if this eigenvector exhibits a localized or extended character. In the former case, one expects a gradual, non-critical-type transition into a glass-type structure. However, in the latter case, a critical-like phase behavior could be observed.

Spatial manipulation of topological defects in nematic shells

It is well known that positions of topological defects (TDs) in liquid crystals can be manipulated experimentally by locally distorting the liquid crystalline (LC) order, as for example by melting induced by optical tweezers. In this work, we study numerically the nematic ordering profiles and the corresponding topological defect configurations in thin nematic liquid crystalline shells controlled by imposed local distortion of LC order. We demonstrate that within curved LC films such manipulations could be strongly affected by local Gaussian curvature if it exhibits strong spatial variations. We use mesoscopic approach in which the shell geometry and LC orientational order are described by curvature of the surface and nematic order parameter tensor. For illustration purposes, we consider LC shells exhibiting spherical topology. We show that on increasing prolateness of shells, which imposes spatially inhomogeneous Gaussian curvature, TDs are relatively strongly "glued" to a local Gaussian curvature.

Effects of orientational order on modulated cylindrical interfaces

Cylindrical interfaces occur in sheared or deformed emulsions and as biological or technological lipid monolayer or bilayer tubules. Like the corresponding spherical droplets and vesicles, these cylinderlike surfaces may host orientational order with n-fold rotational symmetry, for example in the positions of lipid molecules or of spherical nanoparticles. We examine how that order interacts with and induces shape modulations of cylindrical interfaces. While on spherical droplets 2n topological defects necessarily exist and can induce icosahedral droplet shapes, the cylindrical topology is compatible with a defect-free patterning. Nevertheless, once a modulation is introduced by a mechanism such as spontaneous curvature, nontrivial patterns of order, including ones with excess defects, emerge and have nonlinear effects on the shape of the tube. By examining the equilibrium energetics of the system analytically and with a lattice-based Markov chain Monte Carlo simulation, we predict low-temperature morphologies of modulated cylindrical interfaces hosting orientational order. A shape modulation induces a banded pattern of alternatingly isotropic and ordered interfacial material. Furthermore, cylindrical systems can be divided into type I, without defects, and type II, which go through a spectrum of defect states with up to 4n excess defects. The character of the curvature-induced shape transition from unmodulated to modulated cylinders is continuous or discontinuous accordingly.

Cationic vacancies as defects in honeycomb lattices with modular symmetries

Layered materials tend to exhibit intriguing crystalline symmetries and topological characteristics based on their two dimensional (2D) geometries and defects. We consider the diffusion dynamics of positively charged ions (cations) localized in honeycomb lattices within layered materials when an external electric field, non-trivial topologies, curvatures and cationic vacancies are present. The unit (primitive) cell of the honeycomb lattice is characterized by two generators, [Formula: see text] of modular symmetries in the special linear group with integer entries, corresponding to discrete re-scaling and rotations respectively. Moreover, applying a 2D conformal metric in an idealized model, we can consistently treat cationic vacancies as topological defects in an emergent manifold. The framework can be utilized to elucidate the molecular dynamics of the cations in exemplar honeycomb layered frameworks and the role of quantum geometry and topological defects not only in the diffusion process such as prediction of conductance peaks during cationic (de-)intercalation process, but also pseudo-spin and pseudo-magnetic field degrees of freedom on the cationic honeycomb lattice responsible for bilayers.

On the Role of Electrostatic Repulsion in Topological Defect-Driven Membrane Fission

Within a modified Langevin Poisson–Boltzmann model of electric double layers, we derived an analytical expression for osmotic pressure between two charged surfaces. The orientational ordering of the water dipoles as well as the space dependencies of electric potentials, electric fields, and osmotic pressure between two charged spheres were taken into account in the model. Thus, we were able to capture the interaction between the parent cell and connected daughter vesicle or the interactions between neighbouring beads in necklace-like membrane protrusions. The predicted repulsion between them can facilitate the topological antidefect-driven fission of membrane daughter vesicles and the fission of beads of undulated membrane protrusions.

On the Role of Curved Membrane Nanodomains, and Passive and Active Skeleton Forces in the Determination of Cell Shape and Membrane Budding

Biological membranes are composed of isotropic and anisotropic curved nanodomains. Anisotropic membrane components, such as Bin/Amphiphysin/Rvs (BAR) superfamily protein domains, could trigger/facilitate the growth of membrane tubular protrusions, while isotropic curved nanodomains may induce undulated (necklace-like) membrane protrusions. We review the role of isotropic and anisotropic membrane nanodomains in stability of tubular and undulated membrane structures generated or stabilized by cyto- or membrane-skeleton. We also describe the theory of spontaneous self-assembly of isotropic curved membrane nanodomains and derive the critical concentration above which the spontaneous necklace-like membrane protrusion growth is favorable. We show that the actin cytoskeleton growth inside the vesicle or cell can change its equilibrium shape, induce higher degree of segregation of membrane nanodomains or even alter the average orientation angle of anisotropic nanodomains such as BAR domains. These effects may indicate whether the actin cytoskeleton role is only to stabilize membrane protrusions or to generate them by stretching the vesicle membrane. Furthermore, we demonstrate that by taking into account the in-plane orientational ordering of anisotropic membrane nanodomains, direct interactions between them and the extrinsic (deviatoric) curvature elasticity, it is possible to explain the experimentally observed stability of oblate (discocyte) shapes of red blood cells in a broad interval of cell reduced volume. Finally, we present results of numerical calculations and Monte-Carlo simulations which indicate that the active forces of membrane skeleton and cytoskeleton applied to plasma membrane may considerably influence cell shape and membrane budding.

Minimizing isotropic and deviatoric membrane energy - An unifying formation mechanism of different cellular membrane nanovesicle types

Tiny membrane-enclosed cellular fragments that can mediate interactions between cells and organisms have recently become a subject of increasing attention. In this work the mechanism of formation of cell membrane nanovesicles (CNVs) was studied experimentally and theoretically. CNVs were isolated by centrifugation and washing of blood cells and observed by optical microscopy and scanning electron microscopy. The shape of the biological membrane in the budding process, as observed in phospholipid vesicles, in erythrocytes and in CNVs, was described by an unifying model. Taking the mean curvature h and the curvature deviator d of the membrane surface as the relevant parameters, the shape and the distribution of membrane constituents were determined theoretically by minimization of membrane free energy. Considering these results and previous results on vesiculation of red blood cells it was interpreted that the budding processes may lead to formation of different types of CNVs as regards the compartment (exo/endovesicles), shape (spherical/tubular/torocytic) and composition (enriched/depleted in particular kinds of molecules). It was concluded that the specificity of pinched off nanovesicles derives from the shape of the membrane constituents and not primarily from their chemical identity, which explains evidences on great heterogeneity of isolated extracellular vesicles with respect to composition.

One of the amazing properties of a biological membrane is the ability to undergo dramatic changes of its shape. It may exhibit very high curvature and thereby enclose nano-sized compartments that pinch off from the mother membrane and become freely moving cellular nanovesicles (CNVs). CNVs externalize the pieces of the cell and make them available to other cells within the same organism or other organisms. Therefore they have been acknowledged as mediators of communication between microorganisms, plants, animals and human. Furthernore, they dwell on the border between living and non-living things. Recent findings report on heterogeneity of the size and composition of CNVs found in isolates from different biological samples. As communication between cells is involved in many physiological and patophysiological processes, it is of importance to understand the mechanisms of CNVs formation and recognize the natural laws that mainly govern them. We point to an unifying mechanism that explains stability of differently shaped and composed CNVs by taking into account that the biological membrane tends to attain the minimum of its relevant energy. Conveniently, the procedure can be described by a mathematical model which allows for transparent comparison between experimentally induced shapes of membrane-enclosed vesicular structures and numerical calculations.

Topological defects of dipole patchy particles on a spherical surface

We investigate the assembly of dipole-like patchy particles confined to a spherical surface by Brownian dynamics simulations. The surface property of the spherical particle is described by the spherical harmonic Y10, and the orientation of the particle is defined as the uniaxial axis. On a flat space, we observe a defect-free square lattice with nematic order. On a spherical surface, defects appear due to the topological constraint. As for the director field, four defects of winding number +1/2 are observed, satisfying the Euler characteristic. We have found many configurations of the four defects lying near a great circle. Regarding the positional order for the square lattice, eight grain boundary scars proliferate linearly with the sphere size. The positions and orientations of the eight grain boundary scars are strongly related to the four +1/2 defect cores.

Normal red blood cells' shape stabilized by membrane's in-plane ordering

Red blood cells (RBCs) are present in almost all vertebrates and their main function is to transport oxygen to the body tissues. RBCs' shape plays a significant role in their functionality. In almost all mammals in normal conditions, RBCs adopt a disk-like (discocyte) shape, which optimizes their flow properties in vessels and capillaries. Experimentally measured values of the reduced volume (v) of stable discocyte shapes range in a relatively broad window between v ~ 0.58 and 0.8. However, these observations are not supported by existing theoretical membrane-shape models, which predict that discocytic RBC shape is stable only in a very narrow interval of v values, ranging between v ~ 0.59 and 0.65. In this study, we demonstrate that this interval is broadened if a membrane's in-plane ordering is taken into account. We model RBC structures by using a hybrid Helfrich-Landau mesoscopic approach. We show that an extrinsic (deviatoric) curvature free energy term stabilizes the RBC discocyte shapes. In particular, we show on symmetry grounds that the role of extrinsic curvature is anomalously increased just below the nematic in-plane order-disorder phase transition temperature.

Tubulation pattern of membrane vesicles coated with biofilaments

Narrow membrane tubes are commonly pulled out from the surface of phospholipid vesicles using forces applied either through laser or magnetic tweezers or through the action of processive motor proteins. Recent examples have emerged in which an array of such tubes grows spontaneously from vesicles coated with bioactive cytoskeletal filaments (e.g., FtsZ, microtubule) in the presence GTP or ATP. We show how a soft vesicle deforms as a result of the interplay between its topology, local curvature, and the forces due to filament bundles. We present results from dynamically triangulated Monte Carlo simulations of a closed membrane vesicle coated with a nematic field (the filaments), and we show how the intrinsic curvature of the filaments and their bundling interactions drive membrane tubulation. We predict interesting patterns consisting of a large number of nematic defects that accompany tubulation. A common theme emerges: defect locations on vesicle surfaces are hot spots of membrane deformation activity, which could be useful for vesicle origami. Although our equilibrium model is not applicable to the nonequilibrium shape dynamics exhibited by active microtubule-coated vesicles, we show that some of the features, such as the size-dependent vesicle shape and the number of tubes, can still be understood from our equilibrium model.

Bistable curvature potential at hyperbolic points of nematic shells

Nematic shells are colloidal particles coated with nematic liquid crystal molecules which may freely glide and rotate on the colloid's surface while keeping their long axis on the local tangent plane. We describe the nematic order on a shell by a unit director field on an orientable surface. Equilibrium fields can then be found by minimising the elastic energy, which in general is a function of the surface gradient of the director field. We learn how to extract systematically out of this energy a fossil component, related only to the surface and its curvatures, which expresses a curvature potential for the molecular torque. At hyperbolic points on the colloid's surface, and only there, the alignment preferred by the curvature potential may fail to be a direction of principal curvature. There the fossil energy becomes bistable.

Curvature-driven stability of defects in nematic textures over spherical disks

Stabilizing defects in liquid-crystal systems is crucial for many physical processes and applications ranging from functionalizing liquid-crystal textures to recently reported command of chaotic behaviors of active matters. In this work, we perform analytical calculations to study the curvature-driven stability mechanism of defects based on the isotropic nematic disk model that is free of any topological constraint. We show that in a growing spherical disk covering a sphere the accumulation of curvature effect can prevent typical +1 and +1/2 defects from forming boojum textures where the defects are repelled to the boundary of the disk. Our calculations reveal that the movement of the equilibrium position of the +1 defect from the boundary to the center of the spherical disk occurs in a very narrow window of the disk area, exhibiting the first-order phase-transition-like behavior. For the pair of +1/2 defects by splitting a +1 defect, we find the curvature-driven alternating repulsive and attractive interactions between the two defects. With the growth of the spherical disk these two defects tend to approach and finally recombine towards a +1 defect texture. The sensitive response of defects to curvature and the curvature-driven stability mechanism demonstrated in this work in nematic disk systems may have implications towards versatile control and engineering of liquid-crystal textures in various applications.