Nonlinear Methods for Shape Optimization Problems in Liquid Crystal Tactoids

James H. Adler,Anca S. Andrei,Timothy J. Atherton

2024-06-06

Abstract:Anisotropic fluids, such as nematic liquid crystals, can form non-spherical equilibrium shapes known as tactoids. Predicting the shape of these structures as a function of material parameters is challenging and paradigmatic of a broader class of problems that combine shape and order. Here, we consider a discrete shape optimization approach with finite elements to find the configuration of two-dimensional and three-dimensional tactoids using the Landau de Gennes framework and a Q-tensor representation. Efficient solution of the resulting constrained energy minimization problem is achieved using a quasi-Newton and nested iteration algorithm. Numerical validation is performed with benchmark solutions and compared against experimental data and earlier work. We explore physically motivated subproblems, whereby the shape and order are separately held fixed, respectively, to explore the role of both and examine material parameter dependence of the convergence. Nested iteration significantly improves both the computational cost and convergence of numerical solutions of these highly deformable materials.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is how to predict the shapes of liquid - crystal tactoids under different material parameters. Specifically, the article focuses on nematic liquid crystals, which can form non - spherical equilibrium shapes, namely tactoids. Predicting the shapes of these structures as a function of material parameters is a challenging problem and an example of the broader problem of combining shape and order.

### Detailed Explanation:

1. **Research Background**:
- Nematic liquid crystals are anisotropic fluids that can form non - spherical equilibrium shapes, called tactoids, when suspended in the surrounding isotropic fluid.
- These tactoids can exhibit various shapes and orientation - field configurations, depending on their size, elastic properties, and anisotropic surface - tension strength.

2. **Research Objectives**:
- The goal of the article is to develop effective numerical methods based on nonlinear optimization techniques to predict the shape and order problems of tactoids.
- Use the discrete - shape optimization method combined with the finite - element method to find the optimal configurations of tactoids in two - dimensional and three - dimensional cases.
- Describe the order of liquid crystals by using the Landau - de Gennes framework and Q - tensor representation.

3. **Research Methods**:
- **Energy Minimization Problem**: Determine the shape and order of tactoids by solving the constrained energy minimization problem.
- **Optimization Algorithm**: Adopt the quasi - Newton method and the nested - iteration algorithm to efficiently solve the numerical solutions of highly deformable materials.
- **Verification and Comparison**: Conduct numerical verification through benchmark solutions and compare with experimental data and earlier works.

4. **Main Contributions**:
- Propose an integrated optimization method that simultaneously determines the shape and orientation configuration, ensuring physical validity.
- Use the nested - iteration technique to significantly improve computational efficiency and convergence speed.
- Explore the influence of shape and order on the dependence of material parameters.

5. **Application Prospects**:
- Tactoids have broad application potential in emerging technologies due to their ability to change shape and adapt to complex geometric structures, such as enhancing the performance of liquid - crystal displays, drug - carrier design, and developing adaptive - stiffness materials.

### Formula Summary:

- **Q - tensor Representation**:
\[
Q = S(n\otimes n-\frac{I}{d})
\]
where \(S\) is a scalar field, \(n\) is a unit vector, and \(I\) is an identity matrix.

- **Free - Energy Functional**:
\[
F(x, Q)=\int_M[\text{tr}(Q^2)+\frac{2b}{3}\text{tr}(Q^3)+\frac{c}{2}(\text{tr}(Q^2))^2+\frac{L_1}{2}|\nabla Q|^2]dx+\int_{\partial M}[\sigma-\frac{W}{2}\text{tr}((Q - Q_s)^2)]dS
\]

- **Volume Constraint**:
\[
C(x)=\int_Mdx - V_0 = 0
\]

Through these methods and formulas, this paper aims to provide an effective method to predict and understand the shape and order of nematic - liquid - crystal tactoids, providing theoretical support for research and applications in related fields.