🔬 materials science

Contrasting magnetic anisotropy in CrCl3 and CrBr3: A first-principles study

This first-principles study reveals that the contrasting easy magnetization axes in layered CrCl3 (in-plane) and CrBr3 (out-of-plane) arise from the interplay between shape anisotropy and spin-orbit coupling-induced magnetocrystalline anisotropy, where the distinct spatial distribution, hybridization, and spin-selection rules of the halogen p orbitals determine whether the net anisotropy energy favors an in-plane or out-of-plane orientation.

Original authors: Jiazhuang Si, Shuyuan Liu, Bing Wang, Chongze Wang, Fengzhu Ren, Yu Jia, Jun-Hyung Cho

Published 2026-03-03

📖 4 min read☕ Coffee break read

CC0 1.0

Original authors: Jiazhuang Si, Shuyuan Liu, Bing Wang, Chongze Wang, Fengzhu Ren, Yu Jia, Jun-Hyung Cho

Original paper dedicated to the public domain under CC0 1.0 (http://creativecommons.org/publicdomain/zero/1.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine you have two identical-looking teams of tiny magnets (atoms) arranged in a flat, honeycomb pattern. These are the materials CrCl₃ (Chromium Chloride) and CrBr₃ (Chromium Bromide). They are made of the same "captain" atom (Chromium) but have different "assistant" atoms (Chlorine vs. Bromine).

You might expect them to behave the same way, but they have a secret disagreement about which way they want to point their magnetic "compass needles."

CrCl₃ wants its magnets to lie flat on the table (in-plane).
CrBr₃ wants its magnets to stand straight up like a flagpole (out-of-plane).

This paper is a detective story that uses super-computers to figure out why these two cousins act so differently. Here is the breakdown in simple terms:

The Two Forces at Play

To understand the magnets' direction, we have to look at a tug-of-war between two invisible forces:

The "Shape" Force (Shape-MAE): Imagine a stack of coins. It's much easier to push the stack sideways than to push it straight down through the table. Because these materials are flat layers, their shape naturally pushes the magnets to lie flat. This is the "Shape Force."
The "Atomic Spin" Force (SOC-MAE): This is a quantum mechanical force caused by how electrons spin and orbit around the nucleus. It's like a tiny internal compass inside every atom that wants to point in a specific direction based on the atom's internal structure.

The final direction the magnet points depends on which force wins the tug-of-war.

The Mystery: Why do they choose different winners?

In both materials, the "Shape Force" is trying to pull the magnets flat. So, the question is: Why does the "Atomic Spin Force" win in one but lose in the other?

The answer lies in the "assistants" (the Chlorine vs. Bromine atoms) and how they talk to the "captain" (Chromium).

Case 1: CrCl₃ (The Chlorine Team)

The Assistants: Chlorine atoms are small and hold their electrons very tightly. They are like shy, introverted neighbors who stay close to their own houses.
The Interaction: Because they are so tight and small, they don't mix well with the Chromium captain. When the quantum forces try to align the magnets, the Chlorine atoms create a confusing signal.
The Result: The forces pulling the magnets "up" and the forces pulling them "down" are almost equal but opposite. It's like two people pushing a car from opposite sides with equal strength—the car doesn't move. The "Atomic Spin Force" cancels itself out.
The Winner: Since the internal spin force is weak and canceled out, the Shape Force wins. The magnets lie flat.

Case 2: CrBr₃ (The Bromine Team)

The Assistants: Bromine atoms are larger and their electrons are more spread out and "loose." They are like outgoing, energetic neighbors who hang out in the street and mix with everyone.
The Interaction: Because they are larger and more active, they mix strongly with the Chromium captain. They also have a heavier atomic weight, which makes their internal "spin" much stronger (like a heavier, more powerful gyroscope).
The Result: Instead of canceling each other out, the forces from the Bromine atoms work together. They create a massive, unified push in the "up" direction.
The Winner: The Atomic Spin Force becomes so strong that it overpowers the Shape Force. The magnets stand straight up.

The Big Picture

The scientists discovered that the difference isn't just about the size of the atoms, but about how their electrons behave:

In Chlorine, the electrons are "localized" (stuck in one spot), leading to a confusing mix of signals that cancel out.
In Bromine, the electrons are "delocalized" (spread out), leading to a strong, unified signal that points upward.

Why Does This Matter?

This is like learning the secret recipe for building better computer chips. If we want to build future computers that use magnetism instead of electricity (spintronics), we need to know exactly how to make magnets point the way we want them to.

By understanding that changing the "assistant" atom from Chlorine to Bromine flips the magnetic switch, scientists can now design new materials with custom magnetic properties. It's like realizing that if you swap a wooden gear for a steel one in a clock, the whole clock runs differently. This knowledge helps engineers build faster, smaller, and more efficient devices for the future.

1. Problem Statement

Two-dimensional (2D) van der Waals (vdW) magnetic materials, such as chromium trihalides, are promising candidates for spintronic applications. However, a fundamental challenge in this field is understanding the microscopic origins of Magnetic Anisotropy Energy (MAE), which determines the Easy Magnetization Axis (EMA).

The Paradox: While the Mermin-Wagner theorem suggests long-range magnetic order is forbidden in isotropic 2D systems at finite temperatures, these materials exhibit robust ferromagnetism due to MAE, which opens a spin-wave gap.
The Specific Issue: Layered chromium trihalides CrCl $_3$ and CrBr $_3$ possess nearly identical crystal structures and electronic configurations (high-spin $t_{2g}^3$ $t_{2 g}^{3}$ ). Yet, they exhibit contrasting EMAs:
- CrCl $_3$ : Prefers an in-plane EMA.
- CrBr $_3$ : Prefers an out-of-plane EMA.
Knowledge Gap: Previous studies have largely focused on monolayers or identified the need to include shape anisotropy. However, the specific orbital-resolved mechanisms, particularly the role of halogen $p$ -orbitals, spin-orbit coupling (SOC) strength, and hybridization differences between Cl and Br in the bulk phase, remain incompletely understood.

2. Methodology

The authors employed first-principles Density Functional Theory (DFT) calculations to investigate the bulk phases of CrCl $_3$ and CrBr $_3$ .

Software & Basis: Vienna ab initio Simulation Package (VASP) with a plane-wave basis set and Projector Augmented-Wave (PAW) potentials.
Exchange-Correlation: Generalized Gradient Approximation (GGA-PBE) with the DFT+U method (Dudarev formalism, $U_{eff} = 3$ eV) to correctly describe localized Cr 3 $d$ electrons.
Interactions: Van der Waals (vdW) interactions were accounted for using the DFT-D3 scheme.
Calculations:
- Structural relaxation until forces were $<0.001$ eV/Å.
- Calculation of Spin-Orbit Coupling induced Magnetocrystalline Anisotropy Energy (SOC-MAE) via total energy differences between magnetization directions.
- Calculation of Shape Magnetic Anisotropy Energy (shape-MAE) arising from long-range dipole-dipole interactions (summed over a 1000 Å cutoff).
- Orbital-resolved analysis: Decomposition of MAE contributions by atom (Cr vs. Halogen) and orbital ( $d$ vs. $p$ ), and inter-orbital coupling channels using second-order perturbation theory.

3. Key Contributions

The paper provides a comprehensive microscopic explanation for the contrasting magnetic behaviors, moving beyond simple structural comparisons to an orbital-selective mechanism:

Decoupling MAE Components: The study explicitly separates and quantifies the competition between SOC-MAE (favoring out-of-plane in both, but stronger in CrBr $_3$ ) and shape-MAE (favoring in-plane due to the 2D geometry).
Identification of the Halogen $p$ -Orbital Role: It establishes that while Cr $d$ -orbitals contribute similarly in both compounds, the halogen $p$ -orbitals are the decisive factor.
Spin-Selection Rules: The paper elucidates how spin-flip vs. spin-conserving SOC interactions, mediated by the specific characteristics of Cl 3 $p$ vs. Br 4 $p$ orbitals, dictate the sign and magnitude of the net MAE.

4. Key Results

A. Structural and Electronic Properties

Lattice & Bonding: CrBr $_3$ has larger lattice parameters and interlayer spacing due to the larger Br radius. The Cr–Br bond is more covalent with weaker charge transfer, whereas Cr–Cl is more ionic with stronger charge transfer.
Band Structure: Both are semiconductors. CrCl $_3$ has a larger band gap (2.46 eV) compared to CrBr $_3$ (1.84 eV).
SOC Impact: Including SOC significantly reduces the band gap of CrBr $_3$ (by 0.21 eV) due to the strong atomic SOC of Br, while CrCl $_3$ remains largely unaffected.

B. Total MAE and EMA

CrCl $_3$ :
- SOC-MAE: +15 $\mu$ eV (favors out-of-plane).
- Shape-MAE: -23 $\mu$ eV (favors in-plane).
- Net Result: Shape-MAE dominates $\rightarrow$ In-plane EMA (Total MAE = 8 $\mu$ eV).
CrBr $_3$ :
- SOC-MAE: +97 $\mu$ eV (favors out-of-plane).
- Shape-MAE: -20 $\mu$ eV (favors in-plane).
- Net Result: SOC-MAE dominates $\rightarrow$ Out-of-plane EMA (Total MAE = 77 $\mu$ eV).
Observation: The SOC-MAE in CrBr $_3$ is more than six times larger than in CrCl $_3$ .

C. Microscopic Origin of the Difference

The study reveals that the difference stems from the interorbital coupling of the halogen $p$ -orbitals:

In CrCl $_3$ (Localized 3 $p$ ):
- Weak SOC and localized orbitals promote spin-flip interactions ( $\sigma \neq \sigma'$ ).
- Two dominant channels, $(p_x, p_y)$ and $(p_y, p_z)$ , contribute with opposite signs (negative and positive, respectively) and similar magnitudes.
- Result: Significant cancellation leads to a small net SOC-MAE, allowing shape anisotropy to win.
In CrBr $_3$ (Delocalized 4 $p$ ):
- Stronger SOC and delocalized orbitals promote spin-conserving interactions ( $\sigma = \sigma'$ ).
- The $(p_x, p_y)$ channel provides a large positive contribution, while the $(p_y, p_z)$ channel provides a small negative contribution.
- Result: The cancellation is suppressed; the large positive term dominates, yielding a sizable net SOC-MAE that overcomes shape anisotropy.

5. Significance

Fundamental Insight: The work demonstrates that magnetic anisotropy in layered semiconductors is not solely determined by the transition metal (Cr) but is critically governed by the ligand (halogen) orbital character, specifically the interplay between orbital delocalization, SOC strength, and spin-selection rules.
Design Principles: It provides a roadmap for engineering magnetic properties in 2D materials. By selecting ligands with specific orbital delocalization and SOC strengths, one can tune the EMA from in-plane to out-of-plane.
Application: These findings are crucial for the development of 2D spintronic and magnetoelectronic devices where controlling the direction of magnetization is essential for data storage and logic operations.

In summary, the paper resolves the contrasting magnetic anisotropy in CrCl $_3$ and CrBr $_3$ by showing that CrBr $_3$ 's out-of-plane preference is driven by strong, spin-conserving SOC interactions in delocalized Br 4 $p$ orbitals, whereas CrCl $_3$ 's in-plane preference results from the cancellation of competing spin-flip interactions in localized Cl 3 $p$ orbitals.