Ferronematics: integrality of the limiting interface in the strong coupling regime

Technical Summary: Integrality of the Limiting Interface in Ferronematics

Problem Statement
This paper investigates the asymptotic behavior of a vectorial energy functional modeling thin films of ferronematics—composite materials consisting of magnetic nanoparticles dispersed in a nematic liquid crystal host. The model, proposed in [8], involves two order parameters: a reduced $Q$-tensor field $Q$ (describing the nematic component) and a vector field $M$ (representing the average polarization of the magnetic particles). The energy functional $F_\varepsilon$ contains a coupling term that promotes alignment between $Q$ and $M$.

The study focuses on the limit as a small parameter $\varepsilon \to 0$, representing a large domain limit relative to the correlation length. Previous works [12, 13] established that for critical pairs $(Q_\varepsilon, M_\varepsilon)$, the energy concentrates on distinct singular sets: a finite set of points for the $Q$-component and a union of segments (an $H^1$-rectifiable varifold) for the $M$-component. However, the multiplicity of this limiting varifold remained undetermined. While integrality (integer multiplicity) is a characteristic feature of scalar Allen-Cahn and Ginzburg-Landau theories, it is generally false for vectorial problems unless specific structural conditions are met. The central question addressed is whether the limiting varifold associated with the $M$-component possesses integer multiplicity in this coupled system.

Methodology
The authors refine the analysis of non-minimizing critical points of the functional $F_\varepsilon$ under specific boundary conditions. The methodology proceeds through several key analytical steps:

1. Decomposition and Auxiliary Functionals: The vector field $M_\varepsilon$ is decomposed relative to the eigenframe of the $Q$-tensor field. Specifically, $M_\varepsilon$ is projected onto the eigenvector $n_\varepsilon$ (associated with the positive eigenvalue of $Q_\varepsilon$) and its orthogonal complement. This yields components $u_{1,\varepsilon} = M_\varepsilon \cdot n_\varepsilon$ and $u_{2,\varepsilon} = M_\varepsilon \cdot n_\varepsilon^\perp$.
2. Quantitative Convergence of $Q$: A crucial preliminary result (Lemma 1.1) establishes a quantitative rate of convergence for the modulus of $Q_\varepsilon$ and its gradient in $W^{1,2}_{loc}$ away from the singular set of $Q$. This improves upon previous estimates in [12] and is essential for controlling the perturbation terms in the $M$-equation.
3. Structural Analysis of the Coupling: The authors analyze the equation satisfied by the orthogonal component $u_{2,\varepsilon}$. They demonstrate that under the condition $\sqrt{2\beta} > 1$ (where $\beta$ is the coupling constant), the equation for $u_{2,\varepsilon}$ undergoes a structural change that forces a rapid decay of this component. Specifically, Proposition 2.1 proves that the $W^{1,2}$-norm of $u_{2,\varepsilon}$ decays at a power rate with respect to $\varepsilon$ away from the singular set of $Q$.
4. Reduction to Scalar Theory: Due to the rapid decay of $u_{2,\varepsilon}$, the system effectively reduces to a scalar problem for the component $u_{1,\varepsilon}$. The authors rewrite the Euler-Lagrange equation for $u_{1,\varepsilon}$ as a perturbed Allen-Cahn equation. The perturbation terms are shown to be sufficiently small (in $L^2$ and $L^1$ senses) to satisfy the assumptions of the integrality results established by De Philippis and Pigati [26] for scalar Allen-Cahn systems with $L^2$-type perturbations.
5. Integrality Inference: By applying the machinery from [26] to the rescaled scalar component, the authors deduce that the limiting density of the energy concentration is an integer multiple of a material-dependent constant $\sigma_\beta$.

Key Results
The main result, Theorem 1, states that if the coupling constant satisfies $\sqrt{2\beta} > 1$, the limiting varifold $V_\star$ associated with the $M$-component has integer multiplicity after rescaling by the constant $\sigma_\beta$. Specifically, the density function $v_\star$ of the varifold satisfies:
$$\sigma_\beta^{-1} v_\star(x) \in \mathbb{N} \quad \text{for } H^1\text{-almost every } x.$$
Here, $\sigma_\beta = \frac{\sqrt{2}}{3}(1 + \sqrt{2\beta})^{3/2}$.

The proof relies on showing that the orthogonal component $u_{2,\varepsilon}$ does not concentrate energy in the limit, effectively decoupling the vectorial nature of the problem in the strong coupling regime. The limiting potential energy measure is shown to be absolutely continuous with respect to the $H^1$-measure on the singular set, with a density that is an integer multiple of $\sigma_\beta$.

Significance and Claims
The paper claims to provide the first integrality result for a limiting varifold obtained in a vectorial Allen-Cahn problem involving a non-trivial, physically relevant potential energy that is not decoupled at the $\varepsilon$-level. Furthermore, this result holds for non-minimizing critical points, extending beyond the scope of standard minimization arguments.

The authors emphasize that the condition $\sqrt{2\beta} > 1$ is critical; it ensures a fast, quantitative decay of the component of $M_\varepsilon$ orthogonal to the preferred direction defined by $Q_\varepsilon$. This decay allows the vectorial problem to be treated quantitatively as a perturbed scalar problem. The paper notes that the behavior in the opposite regime ($\sqrt{2\beta} \leq 1$) remains an open question, where integrality might fail, potentially offering insights into the sufficient conditions for integrality in general vectorial Allen-Cahn systems.

The work builds upon the foundational analysis of ferronematics in [12, 13] and leverages recent breakthroughs in the analysis of vectorial Allen-Cahn systems [4] and scalar perturbed Allen-Cahn systems [26]. It does not propose new experimental setups or applications but rather provides a rigorous mathematical characterization of the topological singularities in this specific physical model.