Topological Defects in Quasi-2D Ferrofluids: Dynamics of Domain Wall Annihilation Mediated by Interfacial Surface Tension Gradients

Introduction

The collective behavior of confined magnetic fluids, such as ferrofluids, provides an accessible physical realization of many complex phenomena observed in condensed matter physics, including phase transitions and defect dynamics. When driven into a metastable state exhibiting magnetic domains, interfaces between these domains—the topological defects—govern the system’s macroscopic evolution toward equilibrium. In our system, a quasi-2D thin film of a specifically engineered ferrofluid exhibits spontaneous magnetization reversal, leading to the formation of extended domain walls separating domains magnetized in opposite directions. The stability and annihilation kinetics of these walls are highly sensitive to boundary conditions and internal energy dissipation mechanisms.

We hypothesize that the annihilation rate of these domain walls is not solely determined by the bulk magnetic energy landscape but is crucially modulated by the interplay between the domain wall velocity and the localized gradients in the interfacial surface tension ($\gamma$). This gradient introduces a non-uniform energy landscape across the wall itself, effectively providing a directional bias for restructuring.

Theoretical Framework and Model Setup

To quantitatively describe this process, we employ a continuum mechanical model coupled with a simplified Landau-Lifshitz description for the localized magnetization vector $\mathbf{M}$. The dynamics are governed by minimizing the total free energy functional, $F$, which includes bulk exchange energy, anisotropy energy, and the interfacial energy term:

$$F = \int dS \left[ \frac{A}{2} (\nabla M)^2 + K (1 - M_z)^2 + \gamma \right]$$

where $A$ relates to exchange stiffness, $K$ accounts for easy-axis anisotropy, and $\gamma$ is the effective surface energy density term. The evolution of the domain wall velocity, $v_w$, is coupled to the driving magnetic field gradient, $\nabla B$, and the local curvature of the interface, $\kappa$, such that the equation of motion balances viscous drag against magnetic driving forces. A simplified expression for the velocity driven by the field gradient is given by:

$$v_w = \mu_0 \frac{M_s}{\eta} (B_{\text{ext}} - B_{\text{depolarizing}}) - \zeta \frac{\partial \gamma}{\partial x} \label{eq:velocity}$$

where $\eta$ is the fluid viscosity and $\zeta$ is the mobility coefficient describing the coupling strength between the magnetic state and the interfacial tension dynamics [1].

Experimental Observations and Data Acquisition

Experiments were conducted using a custom-built microfluidic cell maintaining temperature control within $\pm 0.01 K$. The quasi-2D film thickness ($h$) was held constant at $100 \text{ nm}$. Domain formation was visualized using high-speed Lorentz microscopy capable of resolving domain structures down to $50 \text{ nm}$. We systematically varied the external magnetic field gradient, $\nabla B$, applied orthogonal to the film plane, while monitoring the decay of domain wall density, $\rho_w(t)$.

Table 1. Decay of Domain Wall Density ($\rho_w$) vs. Time ($t$) under Varying Field Gradients ($\nabla B$)

Observations clearly show a non-linear dependence of the annihilation time constant ($\tau_{1/2}$) on the driving gradient, consistent with the viscous drag term becoming increasingly negligible compared to the driving magnetic torque as $\nabla B$ increases [2].

Discussion: The Role of Surface Tension Anisotropy

The critical finding is the quantitative correlation between the imposed $\nabla B$ and the observed annihilation kinetics. At low gradients, the annihilation follows a quasi-diffusive relaxation governed by thermal fluctuations overcoming pinning sites ($\tau_{1/2} \propto 1/(\nabla B)^2$). However, as $\nabla B$ exceeds a threshold value, $B_{\text{crit}}$, the domain wall enters a ballistic regime characterized by rapid, directed motion. This transition correlates precisely with the condition where the kinetic energy imparted by the magnetic field matches the characteristic energy scale associated with the surface tension gradient, $\Delta E_{\gamma} \propto abla (\gamma)$.

This implies that the surface tension gradient acts effectively as a chiral bias along the domain wall trajectory, biasing the annihilation pathway towards low-energy configurations much faster than simple bulk diffusion would predict. This mechanism provides a unique example of how interfacial energy scales can impose global dynamical constraints on mesoscopic structures [1].

Conclusion

We successfully modeled the annihilation dynamics of topological domain walls in quasi-2D ferrofluids, demonstrating that the controlled interplay between external magnetic gradients and anisotropic interfacial tension gradients dictates the relaxation timescale of the system. The ballistic regime observed at high $\nabla B$ suggests a robust mechanism for coherently directing defect motion using tuned boundary conditions. Future work will explore the possibility of reversing this process to induce stable, long-lived domain structures by balancing the driving field with carefully engineered pinning potentials.