Liquid crystals (LCs), which combine fluidity with molecular alignment, exhibit complex phenomena such as defects and reconfigurable structures under electric fields. Modeling these systems requires solving partial differential equations (PDEs) for the Q-tensor—a mathematical object describing LC order via a scalar order strength and a director vector (average molecular alignment). We present a Julia-based software package for simulating LC behavior, leveraging the language’s high-performance capabilities and growing ecosystem. The package streamlines tensor operations critical to the Landau-de Gennes PDE, such as contractions of rank-3 tensors, to simplify code via Einstein notation. This approach enhances readability while maintaining computational efficiency. We demonstrate the framework’s utility by simulating electric field-driven director reorientation near particle surfaces, a precursor to generate 3D solitons—non-dispersive waves with applications in energy and photonic systems. By integrating Julia’s strengths in numerical computing with modular, extensible design, this work advances computational tools for soft matter physics. The package bridges theoretical models (e.g., flexoelectric/dielectric effects) to real-world LC technologies, offering a platform for future innovations.
