[Paper Review] Two-Dimensional Modeling of Soft Ferromagnetic Films

Antonio DeSimone, Robert V. Kohn+3 more|arXiv (Cornell University)|May 2, 2000

Advanced Mathematical Modeling in Engineering4 citations

TL;DR

This paper presents a two-dimensional asymptotic model for soft ferromagnetic films that accurately predicts magnetization patterns under in-plane magnetic fields, even when the field penetrates the film. By reducing the micromagnetic energy to a convex variational problem for charge density and using a robust numerical scheme to reconstruct unit-length magnetization, the theory achieves excellent agreement with experimental observations across all field strengths, including the penetration regime where prior models fail.

ABSTRACT

We examine the response of a soft ferromagnetic film to an in-plane applied magnetic field. Our theory, based on asymptotic analysis of the micromagnetic energy in the thin-film limit, proceeds in two steps: first we determine the magnetic charge density by solving a convex variational problem; then we construct an associated magnetization field using a robust numerical method. Experimental results show good agreement with the theory. Our analysis is consistent with prior work by van den Berg and by Bryant and Suhl, but it goes much further; in particular it applies even for large fields which penetrate the sample.

Motivation & Objective

• To develop a theoretical framework that extends prior 2D models of soft ferromagnetic films beyond the weak-field regime where field penetration occurs.
• To provide a mathematically rigorous and numerically robust method for predicting equilibrium magnetization configurations in thin films.
• To clarify the validity regime of 2D models and their connection to full micromagnetics.
• To identify physical quantities—such as charge density, field penetration region, and magnetization in penetrated zones—that should exhibit little or no hysteresis.

Proposed method

• Perform asymptotic analysis of the micromagnetic energy in the thin-film limit, identifying dominant energy contributions at leading and second order.
• Formulate a convex variational problem for the magnetic charge density by minimizing the reduced energy functional under unit-length and divergence-free constraints.
• Use an interior point method with a self-concordant barrier to numerically solve the convex minimization problem for the initial magnetization field.
• Apply the level set method to compute the viscosity solution of the Hamilton-Jacobi equation, ensuring a robust and physically plausible unit-length magnetization field.
• Reconstruct the final magnetization field by solving a Hamilton-Jacobi equation with the potential of the penetrated field, ensuring consistency with field penetration and wall-free regions.
• Validate predictions against experimental Kerr microscopy data on Permalloy films under varying field strengths.

Experimental results

Research questions

• RQ1 How can a 2D micromagnetic model be extended to accurately describe soft ferromagnetic films under strong in-plane magnetic fields that penetrate the film?
• RQ2 What is the mathematical and physical basis for the emergence of domain patterns in thin films, particularly in the regime of field penetration?
• RQ3 Which physical quantities in the magnetization response are expected to be nearly hysteresis-free, and why?
• RQ4 How does the proposed numerical scheme ensure robustness and consistency with experimental observations, especially in the absence of explicit wall energy minimization?

Key findings

• The theory predicts that the magnetic charge density, the region of field penetration, and the magnetization within the penetrated region exhibit little or no hysteresis, consistent with experimental observations.
• The model accurately reproduces experimental domain patterns in Permalloy films across all field strengths, including the transition from field expulsion to penetration, with remarkable quantitative agreement.
• The numerical method successfully reconstructs a physically plausible unit-length magnetization field by solving a Hamilton-Jacobi equation via the level set method, favoring configurations with minimal domain walls.
• The model identifies that the magnetostatic energy and exchange energy dominate at leading and second order, while wall energy and anisotropy contribute only at higher order, justifying their neglect in the principal approximation.
• The theory confirms that field penetration and domain wall expulsion are two manifestations of the same physical phenomenon, both governed by the same underlying energy minimization principle.
• The predicted magnetization in the penetrated region is aligned with the gradient of the potential of the penetrated field, as confirmed by level curves in numerical and experimental images.