Fast, scalable, and interactive software for Landau-de Gennes numerical modeling of nematic topological defects

Abstract

Numerical modeling of nematic liquid crystals using the tensorial Landau-de Gennes (LdG) theory provides detailed insights into the structure and energetics of the enormous variety of possible topological defect configurations that may arise when the liquid crystal is in contact with colloidal inclusions or structured boundaries. However, these methods can be computationally expensive, making it challenging to predict (meta)stable configurations involving several colloidal particles, and they are often restricted to system sizes well below the experimental scale. Here we present an open-source software package that exploits the embarrassingly parallel structure of the lattice discretization of the LdG approach. Our implementation, combining CUDA/C++ and OpenMPI, allows users to accelerate simulations using both CPU and GPU resources in either single- or multiple-core configurations. We make use of an efficient minimization algorithm, the Fast Inertial Relaxation Engine (FIRE) method, that is well-suited to large-scale parallelization, requiring little additional memory or computational cost while offering performance competitive with other commonly used methods. In multi-core operation we are able to scale simulations up to supra-micron length scales of experimental relevance, and in single-core operation the simulation package includes a user-friendly GUI environment for rapid prototyping of interfacial features and the multifarious defect states they can promote. To demonstrate this software package, we examine in detail the competition between curvilinear disclinations and point-like hedgehog defects as size scale, material properties, and geometric features are varied. We also study the effects of an interface patterned with an array of topological point-defects.

16 pages, 6 figures, 1 youtube link. The full catastrophe

References in corpus (13)

• Reconfigurable knots and links in chiral nematic colloids
• Steady-state hydrodynamic instabilities of active liquid crystals: Hybrid lattice Boltzmann simulations
• Modeling planar degenerate wetting and anchoring in nematic liquid crystals
• Introduction to Q-tensor theory
• Excitable Patterns in Active Nematics
• Triclinic nematic colloidal crystals from competing elastic and electrostatic interactions
• Shape-controlled orientation and assembly of colloids with sharp edges in nematic liquid crystals
• Sculpting stable structures in pure liquids
• Multiparticle Collision Dynamics for Tensorial Nematodynamics
• Stochastic Rotation Dynamics for Nematic Liquid Crystals
• Particle selection through topographic surface patterns in nematic colloids
• Equilibrium state of a cylindrical particle with flat ends in nematic liquid crystals
• The Search direction Correction makes first-order methods faster