Electric-field induced trends of exchange interactions in transition-metal trilayers

Using density functional theory, this study demonstrates that external electric fields induce a nearly linear, layer-dependent modulation of both pairwise and higher-order exchange interactions in unsupported transition-metal trilayers by altering the spin-dependent local density of states at the Fermi level, while preserving the overall magnetic ground state.

The Gist

Imagine a tiny, three-layered sandwich made of magnetic metals. The bottom slice is Iridium, the middle slice is Iron, and the top slice is a different metal like Palladium, Rhodium, or Ruthenium. This isn't a lunch you can eat; it's a microscopic structure scientists use to study how magnets behave.

The researchers in this paper wanted to see what happens to the "friendship" between the atoms in this sandwich when they zap it with an electric field. In the world of magnetism, atoms have tiny magnetic arrows (spins) that want to point in specific directions relative to their neighbors. Sometimes they want to point the same way (friends), and sometimes they want to point in opposite ways (rivals). The strength of this relationship is called the "exchange interaction."

Here is what the study found, using simple analogies:

1. The Electric Field is Like a Gentle Hand
The scientists applied an electric field (a push or pull on electrons) to this sandwich. They expected the magnetic "friendships" to change drastically, perhaps flipping the whole sandwich from one magnetic state to another.

• The Result: The electric field acted more like a gentle hand adjusting the volume on a radio rather than a hammer smashing the device. The "volume" (the strength of the magnetic connections) went up or down depending on the direction of the field, but the "station" (the fundamental magnetic arrangement) stayed the same. The ground state didn't flip; it just got slightly louder or quieter.

2. The "Volume Knob" Effect
When they turned up the electric field, the magnetic connections changed in a very predictable way, almost like a straight line on a graph.

• The Analogy: Imagine the magnetic bonds are like rubber bands. The electric field stretches or compresses these bands. For the closest neighbors (atoms right next to each other), the stretch was small (a few percent). For neighbors a bit further away, the stretch was much more noticeable (up to 30-40%).
• The Catch: This "stretching" depended heavily on which metal was on the top slice of the sandwich. Changing the top metal from Palladium to Rhodium to Ruthenium changed exactly how the rubber bands reacted to the electric push.

3. The "Team Dynamics" (Higher-Order Interactions)
Usually, we think of magnets as just pairs of atoms talking to each other. But this study looked at more complex conversations where groups of three or four atoms talk at once (called "higher-order interactions").

• The Finding: Even these complex group conversations changed when the electric field was applied. Just like the simple pairs, these group dynamics shifted linearly with the field. This is important because these complex group talks are often what hold special magnetic shapes (like skyrmions, which are tiny, stable magnetic whirlpools) together.

4. Why Did This Happen? (The Electronic Screen)
To understand why the magnetic bonds changed, the researchers looked at the electrons inside the metal.

• The Analogy: Think of the electric field as a strong wind blowing over the surface of the sandwich. The electrons inside the metal act like a crowd of people trying to shield themselves from the wind.
• The Mechanism: The wind pushed the electrons around, specifically changing how many "spin-up" and "spin-down" electrons were hanging out near the surface and in the middle iron layer. It's like the wind rearranged the furniture in the room. Because the magnetic "friendships" depend on how these electrons are arranged, changing the furniture (the electron density) changed the strength of the friendships (the exchange interactions).

5. The Bottom Line
The paper concludes that while the electric field didn't flip the magnetic state of these specific metal sandwiches, it did successfully "tune" the strength of the magnetic connections between atoms.

The authors suggest that because these magnetic connections are the glue holding complex magnetic shapes (like skyrmions) together, being able to tune them with an electric field is a powerful tool. It means we might be able to switch these magnetic shapes on and off or move them around using electricity instead of heat or heavy currents, which is a key goal for future, more efficient data storage devices. However, the paper strictly focuses on the theoretical calculation of these changes in the metal layers and does not claim to have built a working device yet.

Technical Summary

Technical Summary: Electric-field induced trends of exchange interactions in transition-metal trilayers

Problem and Motivation
The manipulation of topological spin structures, such as magnetic skyrmions, is critical for next-generation spintronic and neuromorphic devices. While spin-polarized currents can drive these structures, they suffer from the skyrmion Hall effect (transverse motion) and Joule heating. Electric-field (E-field) assisted switching offers a more energy-efficient alternative by locally tuning magnetic interactions rather than transporting the skyrmion. Although experimental studies have demonstrated E-field control of skyrmion properties in multilayers, a systematic theoretical understanding of how E-fields modify the underlying magnetic exchange interactions—specifically the Heisenberg pairwise exchange and beyond-Heisenberg multi-spin higher-order exchange interactions (HOI)—across different transition-metal systems is lacking. This study addresses this gap by investigating unsupported transition-metal trilayers as model systems for ultrathin films.

Methodology
The authors employed density functional theory (DFT) using the full-potential linearized augmented-plane-wave (FLAPW) method within the fleur code. The study focused on unsupported 4d/Fe/Ir trilayers, specifically Pd/Fe/Ir, Rh/Fe/Ir, and Ru/Fe/Ir, with both face-centered cubic (fcc) and hexagonal close-packed (hcp) stackings of the 4d overlayer. The in-plane lattice constant was fixed to that of Ir (2.70 Å), and interlayer distances were kept consistent with previous studies on Rh/Fe/Ir(111).

To probe the magnetic phase space, the authors calculated the total energy dispersion of spin spirals (1Q states) and multi-Q states (including up-up-down-down/uudd and 3Q states) under external electric fields ranging from 0 to ±1.0 V/Å. The electric field was modeled as a symmetric uniform field generated by charged sheets placed 5.3 Å above and below the trilayer.

Magnetic interactions were extracted by mapping the DFT total energies onto a spin Hamiltonian containing:

1. Heisenberg pairwise exchange constants ($J_{ij}$): Derived from spin spiral energy dispersions.
2. Higher-order exchange (HOI) constants: Including biquadratic ($B_1$), three-site four-spin ($Y_1$), and four-site four-spin ($K_1$) interactions. These were calculated using energy differences between multi-Q states and single-Q (spin spiral) states based on fourth-order perturbation theory.

Key Results

• Stability of Ground States: Across all trilayers and electric field strengths (up to ±1.0 V/Å), the qualitative nature of the energy dispersion remained unchanged. The magnetic ground states (spin spirals for fcc-Pd/Fe/Ir and fcc-Rh/Fe/Ir; ferromagnetic for hcp-Pd/Fe/Ir; row-wise antiferromagnetic for Ru/Fe/Ir variants) were preserved, and no field-induced phase transitions were observed.
• Linear Dependence of Exchange Constants: The pairwise exchange constants ($J_i$) and HOI constants ($B_1, Y_1, K_1$) exhibited a nearly linear dependence on the electric field up to approximately ±0.6 V/Å. Deviations from linearity were observed at higher fields for Rh/Fe/Ir and Ru/Fe/Ir systems.
• Magnitude of Variation:
  • The signs of all interaction constants remained unchanged under the applied fields.
  • Nearest-neighbor exchange constants ($J_1$) showed variations on the order of a few percent.
  • Beyond nearest-neighbor constants and HOI constants exhibited larger variations, reaching up to a few tens of percent (e.g., $J_3$ in Rh/Fe/Ir showed slope magnitudes up to ~100 %/(V/Å)).
  • The specific trends (increase or decrease with field polarity) depended on the specific 4d element (Pd, Rh, Ru) and the stacking order (fcc vs. hcp).
• Electronic Mechanism: The study analyzed spin-dependent screening and local density of states (LDOS). The electric field induces charge redistribution primarily in the $p_z$ orbitals of the surface (4d) and bottom (Ir) layers, extending into the vacuum. Near the atomic layers, the screening is governed by $d$ orbitals. The field-induced shifts in the spin-resolved LDOS, particularly at the Fermi level for the Fe layer, correlate with the modifications in exchange coupling. The Fe magnetic moment remained largely constant (~2.8 $\mu_B$ for Pd/Rh, ~1.75 $\mu_B$ for Ru), indicating that the interaction changes are driven by electronic structure modifications rather than moment quenching.

Significance and Claims
The paper claims to provide a systematic analysis of how external electric fields tune magnetic interactions in transition-metal interfaces. The authors emphasize that:

1. Exchange interactions, including both pairwise and higher-order terms, are significantly modifiable by electric fields in these systems.
2. Since exchange frustration and higher-order exchange are known to influence the stability of skyrmions and antiskyrmions, the demonstrated field-induced tuning of these terms suggests a viable mechanism for controlling topological spin textures.
3. The study serves as a theoretical foundation to motivate experimental investigations into electric-field switching of skyrmions in ultrathin transition-metal films, using the unsupported trilayers as simplified model systems.

The authors explicitly note that while DMI and magnetocrystalline anisotropy (MAE) are crucial for skyrmion stability, they were excluded from this study due to their dependence on substrate layers, focusing instead on the intrinsic exchange interactions of the unsupported trilayers.