Jozef Stefan International Postgraduate School, Jamova 39 , 1000 Ljubljana, Slovenia
Comtron, Trzaska 21, 2000 Maribor, Slovenia
Department of Physics, Case Western Reserve University, Cleveland, Ohio 44106-7079, USA
Department of Physics, Case Western Reserve University, Cleveland, Ohio 44106-7079, USA
Department of Physics, Faculty of Natural Sciences and Mathematics, University of Maribor, Koroska cesta 160, 2000 Maribor, Slovenia
Jozef Stefan Institute, Jamova 39, 1000 Ljubljana, Slovenia
^{1}Jozef Stefan International Postgraduate School, Jamova 39 , 1000 Ljubljana, Slovenia
^{2}Comtron, Trzaska 21, 2000 Maribor, Slovenia
^{3}Department of Physics, Case Western Reserve University, Cleveland, Ohio 44106-7079, USA
^{4}Department of Physics, Faculty of Natural Sciences and Mathematics, University of Maribor, Koroska cesta 160, 2000 Maribor, Slovenia
^{5}Jozef Stefan Institute, Jamova 39, 1000 Ljubljana, Slovenia
Guest Editor: A. Emelyanenko
Beilstein J. Nanotechnol. 2018, 9, 109–118. https://doi.org/10.3762/bjnano.9.13
Received 10 Sep 2017,
Accepted 12 Dec 2017,
Published 10 Jan 2018
Using a Landau–de Gennes approach, we study the impact of confinement topology, geometry and external fields on the spatial positioning of nematic topological defects (TDs). In quasi two-dimensional systems we demonstrate that a confinement-enforced total topological charge of m > 1/2 decays into elementary TDs bearing a charge of m = 1/2. These assemble close to the bounding substrate to enable essentially bulk-like uniform nematic ordering in the central part of a system. This effect is reminiscent of the Faraday cavity phenomenon in electrostatics. We observe that in certain confinement geometries, varying the correlation length size of the order parameter could trigger a global rotation of an assembly of TDs. Finally, we show that an external electric field could be used to drag the boojum fingertip towards the interior of the confinement cell. Assemblies of TDs could be exploited as traps for appropriate nanoparticles, opening several opportunities for the development of functional nanodevices.
Keywords: nanoparticles; nematic liquid crystals; topological charge; topological defects
Topological defects (TDs) [1] represent an interdisciplinary research area [2] that is of high interest for nearly all branches of science. Due to their topological origin they exhibit several universal features that are independent of the microscopic details of the system in which the TDs appear. Their complete understanding might even resolve some of the most intriguing unanswered questions of nature: There are several strong indications that fields represent the basic entities of nature and not fundamental particles [3], which are, in this case, an emergent phenomenon. For example, as far back as 1962, Skyrme [4] developed a theory in which he described hadrons as topological defects in the pion field.
A convenient system in which to study the fundamental behavior of TDs are various liquid crystal (LC) phases [5]. They are relatively easily accessible to various experimental methods [6] due to their unique combination of optical anisotropy and transparency, fluid character, and mechanical softness. In addition, the diversity of LC phases and structures guarantees the existence of many qualitatively different TD structures.
Several recent studies reveal that TDs in LCs could efficiently control the position of assemblies of appropriate nanoparticles (NPs) [7-9]. Note that order parameters, which can host defects, possess two qualitatively different components [10]: an amplitude (also referred to as a hydrodynamic) field, and a symmetry breaking (also referred to as a gauge or nonhydrodynamic) field. If the characteristic size of a nanoparticle is comparable to the amplitude correlation length of an order parameter (which roughly estimates the core size of a defect), and if the nanoparticle does not sufficiently disturb the symmetry breaking field surrounding the core of the defect, then the defect could efficiently trap the NP due to the defect core displacement (DCR) mechanism [11]. In this case, a relatively energetically expensive defect core volume is (at least partially) replaced by the volume of the NP, thereby reducing the overall energy. It has been shown that lattices of orientational (disclinations) [11,12] and translational (dislocations) [9] defects can readily trap such NPs. Furthermore, it was demonstrated that line defects could be exploited to form nanowire-type structures [9,13] consisting of NPs.
In this contribution we study numerically the effects of geometry and an external electric field on the positions of nematic TDs using the Landau–de Gennes mesoscopic approach.
Of interest is the impact of confinement and/or an external electric field on topological defects in a nematic liquid crystal. We use the Landau–de Gennes approach [5] in terms of the tensor order parameter . In its eigenframe it is expressed as , where and λ_{i} are the corresponding eigenvectors and eigenvalues, respectively. We consider uniaxial LCs where the bulk equilibrium ordering is described by the uniaxial tensor
The unit vector points along the local uniaxial direction and is referred to as the nematic director. The uniaxial orientational order parameter, S, quantifies the extent of fluctuations along , and is the unit tensor. If the LC ordering is distorted, the system can exhibit biaxial states. In simulations we study topological defects either in the Cartesian coordinate system (x,y,z) or in the cylindrical coordinate system (r,φ,z). Their coordinate frames are determined by the unit vectors and , respectively. We parameterize the nematic order parameter as
where {q_{1},q_{2},q_{3}} are the variational parameters.
A convenient metric for the degree of biaxiality is the biaxiality parameter [14,15]:
A uniaxial state and configurations with maximum biaxiality are denominated by β^{2} = 0 and β^{2} = 1, respectively.
We write the free energy as the sum of volume and surface integrals:
The condensation (f_{c}), elastic (f_{e}) and external electric field (f_{f}) free energy densities are expressed as [5,15]:
respectively. The quantities A_{0}, B, C are material constants, T* is the supercooling temperature, L is the representative characteristic elastic constant in the single elastic constant approximation, is the external electric field, ε_{0} is the permittivity of free space, and Δε is the dielectric constant anisotropy. We model conditions at the LC confining boundaries either by [15,16]:
or
where w is the surface interaction strength, is the surface normal of the local confinement, and describes preferred nematic ordering of the surface. The surface term given by Equation 8 enforces for w > 0 (w < 0) degenerate tangential (homeotropic) anchoring. On the other hand, the contribution in Equation 9 is minimized for , assuming w > 0. In simulations we consider cases where enforces uniaxial ordering given by Equation 1, where S possesses the bulk equilibrium value, and [5,17]:
This ansatz enforces a topological defect of strength m, where m is an integer multiple of 1/2.
For numerical and presentational convenience [15], we introduce the reduced temperature
where T** is the nematic superheating temperature. In this scaling the bulk phase transition temperature T_{IN} corresponds to
and the bulk degree of uniaxial ordering minimizing Equation 5 can be expressed as
The materials properties of the LC are reflected in various characteristic lengths describing the responses of LC ordering to different perturbations. The relevant lengths for our study are the biaxial correlation length, ξ_{b}, the external field coherence length, ξ_{E,} and the surface extrapolation length d_{e.} We define them as follows:
Note that for scaling purposes we express the latter two distances at T = T**. To estimate their values we choose materials properties of a typical nematic LC, 4-cyano-4'-pentylbiphenyl (5CB), for which [5] A_{0} ≈ 0.3·10^{0} J/(K·m^{3}), B ≈ 4.8·10^{6} J/m^{3}, C ≈ 1.6·10^{7} J/m^{3}, L ≈ 10^{−11} J/m, Δε ≈ 5, T_{IN} ≈ 318 K, T_{IN} − T* ≈ 1.1 K, T** − T_{IN} ≈ 0.18. For temperatures close below T_{IN} it follows ξ_{b} ≈ 20 nm, ξ_{E}(E ≈ 10^{6} V/m) ≈ 1 μm, d_{e}(w ≈ 10^{−4} J/m^{2}) ≈ 1 μm.
We obtained nematic structures for the given boundary conditions by minimizing the total free energy of the system. The resulting Euler–Lagrange equilibrium equations for the variational parameters {q_{1},q_{2},q_{3}} are solved using the standard over-relaxation method, the calculation details of which are given in [15].
We consider thin plane-parallel cells of thickness h. The top and bottom plates are placed at z = 0 and z = h, respectively. We consider the cells either in the Cartesian (x,y,z) or cylindrical geometry (r,φ,z), as illustrated in Figure 1. Accordingly, we use two different boundary conditions, to which we refer as “boundary anchoring condition” (BAC) [18] and “surface anchoring condition” (SAC), respectively.
We use BAC in the “Cartesian” cells. We assume that the cells are relatively thin and the nematic ordering is entirely dominated by conditions at the top “master” plate. At the master plate we impose a strong uniaxial boundary condition of order parameter qiven by Equation 10. The latter is either a circle of radius or a trapezoid that is characterized by a distance R. It is assumed that R is large with respect to the relevant nematic order correlation length. Inside the boundary we allow the nematic tensor frame to freely rotate in the (x,y)-plane. In these simulations the nematic ordering is effectively two-dimensional. Hence, we neglect variations along the z-axis and set .
Regarding SAC we consider thicker cells and permit spatial variations along the z-direction. We perform simulations in the cylindrical coordinate system [19] and we impose the cylindrical symmetry, . At the top plate we enforce uniaxial boundary conditions given by Equation 10, in which we set m = 1. At the bottom plate we enforce homeotropic anchoring conditions using the ansatz in Equation 8. At the lateral boundaries we assume free boundary conditions. These conditions impose a boojum topological defect at the top plate [19,20].
Note that in our simulations we mimic geometric set-ups that could be realized experimentally using, for instance, the atomic force microscope (AFM) scribing method [17]. In a typical experimental set up one confines a nematic LC within a thin plane-parallel cell, where at least one (“master”) surface imposes anchoring conditions inscribed via an AFM stylus [17,18], with a planar degenerate “slave” as the other surface. In Figure 2 we depict an example of a ”master” substrate enforcing 2D topological defects of strength m = ±2, and the corresponding experimentally measured textures using polarized optical microscopy.
In the following we present results of our simulations. We consider structures using the “Cartesian” cells, where we study how different patterns of TDs emerge. Afterwards we focus on an external field defect core structure driven changes in the “cylindrical” cell.
We first study patterns emerging from the BAC boundary condition. We enforce a total topological charge of strength m inside the circular boundary of radius R. At the boundary we strongly impose the nematic ordering defined by Equation 10. The energy-minimized configurations are plotted in Figure 3 and Figure 4.
In Figure 3 the equilibrium biaxiality profiles are plotted, in which cores of TDs are clearly visible. The imposed total charges always decompose into TDs bearing elementary charges m_{0} = 1/2. The fingerprint of the cores of m_{0} = 1/2 TDs is a volcano-like rim where β^{2} = 1 [15,18]. For m = 1, m = 2, m = 4 and m = 6 the patterns exhibit 2, 4, 8 and 12 TDs, respectively. The TDs tend to assemble close to the boundary. The resulting director orientation is plotted in Figure 4: The eigenvectors of with the largest positive eigenvalue (which we set to be ) are plotted. These correspond in the uniaxial limit to the nematic director field. Note that nematic textures are essentially uniaxial, except close to the defects cores. For this reason we henceforth refer to as the nematic director field. The ordering becomes increasingly spatially uniform in the central region with increasing values of m.
This phenomenon is reminiscent of the Faraday cavity effect in conductors. Namely, if one puts electric charges on a conducting body, the charges assemble at its surface and the resulting electric field inside the body vanishes. The Faraday-like behavior in our simulations is clearly visible for cases of R/ξ_{b} >> 1. In our simulations, we set R/ξ_{b} = 30, and for typical LCs ξ_{b} is approximately 20 nm. The absence of an electric field inside the conductor in the electrostatic analogue corresponds in our simulations to a spatially uniform nematic director in the area separated by a distance greater than ξ_{b} from the confining boundary.
We next consider cases in which we change the symmetry of the bounding surface. In Figure 5 and Figure 6 we show nematic configurations for a trapezoid-shaped boundary, through which we enforce m = 3 using Equation 10. In all cases the imposed charge decays into elementary charges m_{0} = 1/2, and the charges assemble at the confining boundary as discussed in the previous subsection. For this specific confinement symmetry we observe changes in the nematic patterns when the ratio η = R/ξ_{b} is decreased, where R corresponds to the bottom length of the trapezoid.
In practice, η can be varied by changing temperature of the sample, which affects ξ_{b}. In this case, a rotation of defect patterns is expected according to our simulations. The changes in patterns reveal that the energy landscape of the system substantially changes with varying η. From the perspective of TDs, the LC configurations reflect the interplay between mutual repulsion among defects and interaction of TDs with the confinement geometry. Note that for sufficiently symmetric confinements, the “rotation” disappears. It is also sensitive to the number of TDs.
We next consider thicker cells and SAC boundary conditions. For sufficiently strong anchoring, a boojum surface defect [19,20] occurs at the top surface. Its structure has been studied in detail in [19]. Its surrounding nematic director field resembles a “classical” half-hedgehog structure. The core structure is relatively complex [19] and is schematically shown in Figure 7. Its core is characterized by a negatively uniaxial finger, surrounded by a shell exhibiting maximal biaxiality, β^{2} = 1. The fingertip is melted due to the topology of the surroundings. To understand this let us consider ideal cylindrically symmetry, which we also adopt in our simulations. The symmetry axis of the defect is uniaxial. Namely, in terms of the parametrization defined by Equation 3 the elastic free energy density includes a term that is linearly proportional with [(3q_{1} + q_{2})^{2} + q_{3}^{3}]/r^{2}. The singularity at r = 0 can be avoided if uniaxial states are introduced, for which 3q_{1} + q_{2} = q_{3} =0. When decreasing the value of z from the top plate, the negative uniaxiality extends to the fingertip (Figure 8). Below the tip the axis is positively uniaxial. The transition from negative to positive uniaxiality requires melting of the nematic ordering in the transition area. A typical spatial variation of the director field, , in the radial direction and order changes of the parameter along the symmetry axis are depicted in Figure 8. In the given case, the fingertip (where the nematic ordering is melted) is located at (r,z) = (0,h − ξ_{f}), where ξ_{f} ≈ 0.8ξ_{b}. For z > h − ξ_{f} the director field is roughly radial everywhere and S(r = 0,z) < 0. For z < h − ξ_{f} the nematic order is positively uniaxial at r = 0 and . With increasing r the director field monotonically increases its departure from the z-axis.
We next apply an external electric field along the z-axis and assume that the LC possesses a negative dielectric anisotropy (Δε < 0). In this case the external field favors the negatively uniaxial part of the boojum. Consequently, it becomes elongated on increasing the field strength, as it is illustrated in Figure 9. To achieve this, the external field must be relatively strong, i.e., ξ_{E} ≈ ξ_{b}, which in a typical LC would correspond to E ≈ 10^{8} V/m.
We studied numerically the impact of geometry, topology, and external fields on patterns and positions of nematic topological defects. In our phenomenological study we used the Landau–de Gennes approach in terms of the nematic tensor order parameter. We found that, in quasi 2D systems, TDs rearrange relatively close to the boundary line, which topologically enforces their presence. Such assemblies of TDs enable the formation of essentially spatially uniform orientational ordering in the central part of the confined nematic. This effect is reminiscent of the Faraday cavity phenomenon in electrostatics. Furthermore, we demonstrated that, for certain confinement geometries (we chose a trapezoid), one could induce collective rotation of a pattern of TDs on changing the defect core size. For example, the latter could be varied by changing the temperature. Furthermore, we demonstrated that one could extend the boojum finger towards the cell interior if an external field is imposed approximately along its symmetry axis for LCs exhibiting negative external field anisotropy.
These mechanisms could be exploited for indirect positional manipulation [7-9,11,12] of certain NPs via positionally controlled TDs. Namely, appropriately surface-decorated NPs could be efficiently trapped within cores of TDs through the defect core displacement mechanism. Therefore, by manipulating the texture of TDs, one could control or reconfigure the positions of trapped NPs, which could be exploited in future nano-scale devices.
SK acknowledges the financial support from the Slovenian Research Agency (research core funding No. P1-0099). PK acknowledges the Slovenian Human Resources and Scholarship Fund. CR was supported by the National Science Foundation’s Condensed Matter Physics program under grant DMR-1505389, and BSM was supported by the National Aeronautics and Space Administration under grant NNX17AC76G.
19. | Kralj, S.; Rosso, R.; Virga, E. G. Phys. Rev. E 2008, 78, 031701. doi:10.1103/PhysRevE.78.031701 |
19. | Kralj, S.; Rosso, R.; Virga, E. G. Phys. Rev. E 2008, 78, 031701. doi:10.1103/PhysRevE.78.031701 |
7. | Kikuchi, H.; Yokota, M.; Hisakado, M.; Yang, H.; Kajiyama, T. Nat. Mater. 2002, 1, 64–68. doi:10.1038/nmat712 |
8. | Karatairi, E.; Rožič, B.; Kutnjak, Z.; Tzitzios, V.; Nounesis, G.; Cordoyiannis, G.; Thoen, J.; Glorieux, C.; Kralj, S. Phys. Rev. E 2010, 81, 041703. doi:10.1103/PhysRevE.81.041703 |
9. | Coursault, D.; Grand, J.; Zappone, B.; Ayeb, H.; Levi, G.; Felidj, N.; Lacaze, E. Adv. Mater. 2012, 24, 1461–1465. doi:10.1002/adma.201103791 |
11. | Rožič, B.; Tzitzios, V.; Karatairi, E.; Tkalec, U.; Nounesis, G.; Kutnjak, Z.; Cordoyiannis, G.; Rosso, R.; Virga, E. G.; Muševič, I.; Kralj, S. Eur. Phys. J. E 2011, 34, 17–27. doi:10.1140/epje/i2011-11017-8 |
12. | Cordoyiannis, G.; Jampani, V. S. R.; Kralj, S.; Dhara, S.; Tzitzios, V.; Basina, G.; Nounesis, G.; Kutnjak, Z.; Tripathi, C. S. P.; Losada-Pérez, P.; Jesenek, D.; Glorieux, C.; Muševič, I.; Zidanšek, A.; Amenitsch, H.; Thoen, J. Soft Matter 2013, 9, 3956–3964. doi:10.1039/c3sm27644a |
1. | Mermin, N. D. Rev. Mod. Phys. 1979, 51, 591–648. doi:10.1103/RevModPhys.51.591 |
5. | Kleman, M.; Lavrentovich, O. D. Soft Matter Physics; Springer: New York, NY, USA, 2003. doi:10.1007/b97416 |
5. | Kleman, M.; Lavrentovich, O. D. Soft Matter Physics; Springer: New York, NY, USA, 2003. doi:10.1007/b97416 |
15. | Kralj, S.; Virga, E. G. J. Phys. A: Math. Gen. 2001, 34, 829–838. doi:10.1088/0305-4470/34/4/309 |
4. | Skyrme, T. H. R. Nucl. Phys. 1962, 31, 556–569. doi:10.1016/0029-5582(62)90775-7 |
15. | Kralj, S.; Virga, E. G. J. Phys. A: Math. Gen. 2001, 34, 829–838. doi:10.1088/0305-4470/34/4/309 |
16. | Nobili, M.; Durand, G. Phys. Rev. A 1992, 46, R6174–R6177. doi:10.1103/PhysRevA.46.R6174 |
3. | Hobson, A. Am. J. Phys. 2013, 81, 211–223. doi:10.1119/1.4789885 |
5. | Kleman, M.; Lavrentovich, O. D. Soft Matter Physics; Springer: New York, NY, USA, 2003. doi:10.1007/b97416 |
2. | Zurek, W. H. Nature 1985, 317, 505–508. doi:10.1038/317505a0 |
14. | Kaiser, P.; Wiese, W.; Hess, S. J. Non-Equilib. Thermodyn. 1992, 17, 153–169. |
15. | Kralj, S.; Virga, E. G. J. Phys. A: Math. Gen. 2001, 34, 829–838. doi:10.1088/0305-4470/34/4/309 |
11. | Rožič, B.; Tzitzios, V.; Karatairi, E.; Tkalec, U.; Nounesis, G.; Kutnjak, Z.; Cordoyiannis, G.; Rosso, R.; Virga, E. G.; Muševič, I.; Kralj, S. Eur. Phys. J. E 2011, 34, 17–27. doi:10.1140/epje/i2011-11017-8 |
9. | Coursault, D.; Grand, J.; Zappone, B.; Ayeb, H.; Levi, G.; Felidj, N.; Lacaze, E. Adv. Mater. 2012, 24, 1461–1465. doi:10.1002/adma.201103791 |
10. | Palffy-Muhoray, P. Phys. Today 2007, 60, 54–60. doi:10.1063/1.2784685 |
9. | Coursault, D.; Grand, J.; Zappone, B.; Ayeb, H.; Levi, G.; Felidj, N.; Lacaze, E. Adv. Mater. 2012, 24, 1461–1465. doi:10.1002/adma.201103791 |
13. | Ok, J. M.; Kim, Y. H.; Lee, T. Y.; Yoo, H.-W.; Kwon, K.; Jung, W.-B.; Kim, S.-H.; Jung, H.-T. Langmuir 2016, 32, 13418–13426. doi:10.1021/acs.langmuir.6b03355 |
7. | Kikuchi, H.; Yokota, M.; Hisakado, M.; Yang, H.; Kajiyama, T. Nat. Mater. 2002, 1, 64–68. doi:10.1038/nmat712 |
8. | Karatairi, E.; Rožič, B.; Kutnjak, Z.; Tzitzios, V.; Nounesis, G.; Cordoyiannis, G.; Thoen, J.; Glorieux, C.; Kralj, S. Phys. Rev. E 2010, 81, 041703. doi:10.1103/PhysRevE.81.041703 |
9. | Coursault, D.; Grand, J.; Zappone, B.; Ayeb, H.; Levi, G.; Felidj, N.; Lacaze, E. Adv. Mater. 2012, 24, 1461–1465. doi:10.1002/adma.201103791 |
6. | Kurik, M. V.; Lavrentovich, O. D. Sov. Phys. - Usp. 1988, 31, 196–224. doi:10.1070/PU1988v031n03ABEH005710 |
11. | Rožič, B.; Tzitzios, V.; Karatairi, E.; Tkalec, U.; Nounesis, G.; Kutnjak, Z.; Cordoyiannis, G.; Rosso, R.; Virga, E. G.; Muševič, I.; Kralj, S. Eur. Phys. J. E 2011, 34, 17–27. doi:10.1140/epje/i2011-11017-8 |
12. | Cordoyiannis, G.; Jampani, V. S. R.; Kralj, S.; Dhara, S.; Tzitzios, V.; Basina, G.; Nounesis, G.; Kutnjak, Z.; Tripathi, C. S. P.; Losada-Pérez, P.; Jesenek, D.; Glorieux, C.; Muševič, I.; Zidanšek, A.; Amenitsch, H.; Thoen, J. Soft Matter 2013, 9, 3956–3964. doi:10.1039/c3sm27644a |
5. | Kleman, M.; Lavrentovich, O. D. Soft Matter Physics; Springer: New York, NY, USA, 2003. doi:10.1007/b97416 |
5. | Kleman, M.; Lavrentovich, O. D. Soft Matter Physics; Springer: New York, NY, USA, 2003. doi:10.1007/b97416 |
17. | Murray, B. S.; Pelcovits, R. A.; Rosenblatt, C. Phys. Rev. E 2014, 90, 052501. doi:10.1103/PhysRevE.90.052501 |
15. | Kralj, S.; Virga, E. G. J. Phys. A: Math. Gen. 2001, 34, 829–838. doi:10.1088/0305-4470/34/4/309 |
15. | Kralj, S.; Virga, E. G. J. Phys. A: Math. Gen. 2001, 34, 829–838. doi:10.1088/0305-4470/34/4/309 |
18. | Kralj, S.; Murray, B. S.; Rosenblatt, C. Phys. Rev. E 2017, 95, 042702. doi:10.1103/PhysRevE.95.042702 |
19. | Kralj, S.; Rosso, R.; Virga, E. G. Phys. Rev. E 2008, 78, 031701. doi:10.1103/PhysRevE.78.031701 |
20. | Volovik, G. E. JETP Lett. 1978, 28, 59–62. |
17. | Murray, B. S.; Pelcovits, R. A.; Rosenblatt, C. Phys. Rev. E 2014, 90, 052501. doi:10.1103/PhysRevE.90.052501 |
17. | Murray, B. S.; Pelcovits, R. A.; Rosenblatt, C. Phys. Rev. E 2014, 90, 052501. doi:10.1103/PhysRevE.90.052501 |
18. | Kralj, S.; Murray, B. S.; Rosenblatt, C. Phys. Rev. E 2017, 95, 042702. doi:10.1103/PhysRevE.95.042702 |
19. | Kralj, S.; Rosso, R.; Virga, E. G. Phys. Rev. E 2008, 78, 031701. doi:10.1103/PhysRevE.78.031701 |
19. | Kralj, S.; Rosso, R.; Virga, E. G. Phys. Rev. E 2008, 78, 031701. doi:10.1103/PhysRevE.78.031701 |
20. | Volovik, G. E. JETP Lett. 1978, 28, 59–62. |
15. | Kralj, S.; Virga, E. G. J. Phys. A: Math. Gen. 2001, 34, 829–838. doi:10.1088/0305-4470/34/4/309 |
18. | Kralj, S.; Murray, B. S.; Rosenblatt, C. Phys. Rev. E 2017, 95, 042702. doi:10.1103/PhysRevE.95.042702 |
© 2018 Kurioz et al.; licensee Beilstein-Institut.
This is an Open Access article under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
The license is subject to the Beilstein Journal of Nanotechnology terms and conditions: (http://www.beilstein-journals.org/bjnano)