🔬 mesoscale physics

An expandable kinetic Monte Carlo platform for modelling electron transport through chiral molecules

This paper presents an expandable kinetic Monte Carlo platform designed to model and compare competing theories of electron transport through chiral molecules, specifically investigating the relationship between Electronic Magnetochiral Anisotropy (eMChA) and Chirality Induced Spin Selectivity (CISS) by simulating voltage-dependent spin filtering effects.

Original authors: Silvia Giménez-Santamarina, Andrés Mora Martínez, Gérliz M. Gutiérrez-Finol, Alejandro Gaita-Ariño

Published 2026-07-21

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Silvia Giménez-Santamarina, Andrés Mora Martínez, Gérliz M. Gutiérrez-Finol, Alejandro Gaita-Ariño

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.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 a world where the tiny particles that power our devices have a secret personality trait: a spin. In the field of spintronics, scientists are trying to use this spin, in addition to the particle's electric charge, to store and process information. It's like upgrading a simple on/off light switch to a dimmer that can also change color, potentially leading to faster, smaller, and more energy-efficient electronics. But there's a catch: most current devices need heavy, magnetic materials to control this spin, which makes them hard to shrink down for use in flexible or biological systems.

Enter the concept of "chirality." In everyday language, chirality is just a fancy word for "handedness." Think of your hands: your left hand is a mirror image of your right, but you can't stack them perfectly on top of each other. Many molecules in nature, like the DNA in your cells, are chiral. Scientists have discovered something magical: when electrons travel through these twisted, chiral molecules, their spin gets sorted out, almost like a bouncer at a club letting only one type of person in. This phenomenon, known as Chirality Induced Spin Selectivity (CISS), suggests we might be able to control electron spins using just the shape of a molecule, without needing any magnets at all. However, there is another related effect called electric Magnetochiral Anisotropy (eMChA) that behaves slightly differently, and scientists are still arguing over whether these two effects are cousins or twins.

To solve this mystery, the authors of this paper have built a new digital playground: a kinetic Monte Carlo simulation. Think of this as a highly sophisticated video game where they program thousands of virtual electrons to run through a virtual chiral molecule. Instead of building a physical lab setup, they created a code that tracks how these electrons hop from one spot to another, assigning them a spin (either "alpha" or "beta") and seeing how the molecule's twist affects their journey.

The team's main finding is that their simulation successfully reproduces the behavior of eMChA. In their virtual world, when they apply a voltage, the electrons don't just flow; they get filtered based on their spin and the direction they are moving. Crucially, the simulation shows that this effect is voltage-dependent. At very low voltages, the effect vanishes, and the electrons flow almost equally in both directions. But as the voltage increases, the "bouncer" gets stricter, creating a clear difference in how the electrons behave depending on whether the voltage is positive or negative. This matches what experimentalists see in real labs: the effect is a nonlinear response that grows with the push of the voltage.

The paper also explores the relationship between CISS and eMChA. While CISS is often observed as a steady spin filter that works even without an external magnetic field, the authors' simulation suggests that eMChA is a specific type of response that emerges when there is a strong, directional flow of charge. They found that in their model, the "magnetic" part of the effect (which usually requires an external magnet) can be mimicked by the motion of the electrons themselves. When electrons rush through the chiral twist, they generate a tiny, internal magnetic influence that interacts with their spin. This explains why the effect is weak or non-existent at low voltages (where electrons just jitter randomly) but becomes strong when the voltage pushes them into a steady stream.

Importantly, the authors are careful to note that this is a simulation, not a final proof of how nature works. They have created a "frugal" and expandable tool that captures the essential physics without getting bogged down in overly complex calculations. They demonstrate that their code can recover standard electrical laws (like Ohm's law) and then layer on the spin effects to see what happens. They explicitly rule out the idea that this effect is a simple, linear response that happens at any voltage; their results show it requires a certain level of driving force to kick in.

Ultimately, this work doesn't claim to have solved the entire mystery of chiral spin transport. Instead, it offers a flexible, efficient way to test different theories. By adjusting the parameters in their code, the researchers can see which combinations of rules produce the experimental results we see in the real world. They suggest that the interplay between the electron's motion, its spin, and the molecule's twist is the key, and their tool allows them to quantify exactly how strong that interaction needs to be to create the observed effects. It's a step forward in understanding how the shape of a molecule can act as a traffic cop for spinning electrons, paving the way for future devices that are lighter, smarter, and magnet-free.

Technical Summary: An Expandable Kinetic Monte Carlo Platform for Modelling Electron Transport Through Chiral Molecules

Problem Statement
The paper addresses the lack of a clear microscopic consensus regarding two prominent phenomena in chiral molecular systems: Chirality Induced Spin Selectivity (CISS) and Electric Magnetochiral Anisotropy (eMChA). While both effects involve the coupling of electron spin, linear momentum, and molecular chirality, it remains an open question whether they stem from fundamentally different mechanisms or represent different experimental manifestations of the same underlying physics. Existing theoretical models (ranging from tight-binding and non-equilibrium Green's function approaches to correlated electron models) often struggle to simultaneously reproduce the magnitude, voltage dependence, and material generality observed in experiments. Furthermore, many models are computationally demanding or fail to capture the specific nonlinear transport signatures associated with eMChA, particularly its behavior at low bias.

Methodology
The authors developed a core kinetic Monte Carlo (kMC) code designed to model electron transport through chiral molecules under an applied voltage. The methodology is characterized by the following features:

Microscopic Phenomenological Approach: The model simulates an ensemble of $N$ identical molecular strands, each hosting a single conducting electron. The system is discretized into spatial regions (partitions) along the molecular backbone.
Stochastic Transport: Electron motion is treated as a stochastic hopping process. Transition probabilities are determined by Boltzmann populations evaluated at the working temperature and the energy differences ( $\Delta E$ ) arising from the applied electric field.
Spin and Chirality Coupling: The model explicitly treats $\alpha$ and $\beta$ spin channels independently. A phenomenological factor is introduced to the hopping probabilities to account for the interplay between spin ( $s$ ), molecular helicity ( $h$ ), and transport direction ( $v$ ). This factor follows a symmetry requirement where the transport response changes sign upon reversal of any of these three quantities.
Voltage-Dependent Mechanism: To model eMChA, the authors assume that the spin-orbit coupling contribution arises from the magnetic field generated by charge motion along a chiral trajectory. Since a net magnetic field requires a net directional current, the model modulates the spin-dependent hopping probability with a hyperbolic tangent function of the voltage ( $\tanh(V)$ ). This ensures the effect vanishes at zero bias (diffusion-dominated regime) and grows with applied voltage.
Validation and Scaling: The code was validated by demonstrating that it naturally recovers Ohm's law and standard resistivity scaling (linear dependence on voltage, inverse dependence on length) without explicitly imposing these macroscopic relations. The model operates in a linear response regime by choosing appropriate time steps and partition densities.

Key Contributions and Results

Development of an Expandable kMC Framework: The authors present a modular, computationally efficient code that treats spin and charge transport on discrete networks. Unlike previous approaches that may require complex Hamiltonians or explicit leads, this model uses a "frugal" phenomenological approach that can be systematically expanded to include additional mechanisms (e.g., spin-flip, vibronic coupling).
Reproduction of eMChA Signatures: The simulation successfully reproduces the qualitative features of electric magnetochiral anisotropy. Specifically, the model generates a voltage-dependent spin filtering effect that:
- Vanishes at low bias (linear response limit).
- Exhibits asymmetry between positive and negative voltages.
- Shows distinct resistance behaviors for $\alpha$ and $\beta$ spins depending on the molecular handedness and bias polarity.
Distinction from CISS: The results highlight a key difference between the simulated eMChA behavior and typical CISS observations. While the model shows a vanishing effect at low voltage (consistent with eMChA theory), the authors note that CISS experiments often report robust spin filtering even at low bias. The paper posits that this discrepancy may require future model iterations to resolve.
Quantitative Parameter Extraction: The framework allows for the calculation of the eMChA interaction parameters ( $\gamma$ and $\Delta R$ ) as defined by Pop et al. The simulations confirm that the anisotropy scales linearly with the magnetic field ( $B$ ) at low fields and inversely with temperature ( $T$ ) at high temperatures, consistent with Boltzmann population statistics of spin states.

Significance and Claims
The paper claims that this work provides a foundational, "microscopic, frugal, and expandable" tool for investigating chiral spin transport. Its primary significance lies in its ability to:

Bridge the gap between microscopic assumptions and experimentally accessible observables (such as I-V curves and resistance anisotropy).
Demonstrate that the specific nonlinear voltage dependence of eMChA can emerge naturally from stochastic transport dynamics coupled with a phenomenological spin-chirality interaction, without requiring complex external magnetic field terms in the simulation logic.
Offer a platform for systematically testing theoretical mechanisms (such as exchange interactions or vibronic effects) against experimental data.

The authors remain modest regarding the current scope, acknowledging that the code represents a "core" implementation. They emphasize that while the model successfully recovers the expected behavior for eMChA and standard resistive transport, further development is needed to fully reconcile the low-voltage behavior of the simulation with the robust low-bias filtering observed in CISS experiments. The tool is presented not as a final solution, but as a flexible framework to facilitate the systematic incorporation of emerging theoretical mechanisms.

Technical Summary: An Expandable Kinetic Monte Carlo Platform for Modelling Electron Transport Through Chiral Molecules

More like this