Open Access Article
Tahereh Azargoshasb
,
Paria Kamkar and
Hadi Heidari
*
Microelectronics Lab, James Watt School of Engineering, University of Glasgow, Glasgow, UK
First published on 28th October 2025
Double-barrier magnetic tunnel junctions (DB-MTJs) present a promising platform for biosensing applications due to their improved sensitivity and capability for label-free, non-invasive detection. This work models a high-performance iron-free DB-MTJ configuration: Co2MnSi–MgO–MoTe2–MgO–Co2MnSi. By Utilizing density functional theory (DFT) and the nonequilibrium Green's function (NEGF) technique, we computed the tunnelling magnetoresistance (TMR) ratio, transmission spectra, density of states (DOS), and bandgaps for the Co2MnSi–MgO–MoTe2–MgO–Co2MnSi system. The DB-MTJ features iron-free Co2MnSi Heusler alloy electrodes and a MoTe2 transition metal dichalcogenide (TMDC) spacer, achieving a tunneling magnetoresistance (TMR) ratio of 1226%, alongside other configurations with TMR ratios ranging from 22% to 667%. The high TMR, driven by the tunable electronic properties of MoTe2 and the half-metallic nature of Co2MnSi, enables exceptional sensitivity to interfacial modifications, making it ideal for detecting biomarkers. The application of iron-free materials increases the electrodes' stability and spin polarization and improves the sensitivity of the DB-MTJs upon interface variations. This research points out the high-potential applications of Heusler-based and iron-free DB-MTJ's as sensitive and label-free biosensors.
MTJ devices have garnered considerable research interest in recent years due to their potential in spintronic applications. An MTJ typically consists of two ferromagnetic (FM) electrodes separated by a thin non-magnetic insulating or semiconducting (SC) spacer layer, which facilitates quantum mechanical tunneling of spin-polarized electrons.9,10 The tunneling magnetoresistance (TMR) effect in these structures is strongly dependent on the relative spin alignment of the FM layers parallel (P) or antiparallel (AP) enabling the electrical readout of magnetically encoded information.11 The TMR ratio is defined as the relative difference in electrical conductance between the P and AP configurations, normalized to the lower conductance state. Remarkably, TMR ratios exceeding 3700% have been experimentally demonstrated in Fe/MgO/Fe and FeCo/MgO/FeCo systems, establishing these junctions as leading candidates for practical device applications.12–14 Theoretically, TMR ratios can approach infinity when half-metallic ferromagnetic electrodes are employed, due to their complete spin polarization at the Fermi level, making them highly desirable for enhancing device sensitivity.
Heusler alloys, with the general chemical formula X2YZ, where X and Y are typically transition metals and Z is a main-group element exhibit a wide range of tunable electrical and magnetic properties, including half-metallic, semiconducting, and metallic behaviors, depending on their specific composition.15,16 Among them, cobalt-based full Heusler alloys, such as Co2Fe(Si,Al), Co2(Cr,Fe)Al, and Co2MnSi, have attracted significant attention as promising electrode materials due to their high spin polarization and mechanical and thermodynamic stability in spintronic applications.17 These alloys are often integrated with thin MgO insulating barriers to form MTJs.18,19 The strong appeal of Co-based full Heusler alloys arises from their high spin polarization and elevated Curie temperatures (Tc), which are essential for maintaining magnetic stability and spin coherence at room temperature.20 These features make them promising candidates for usage in spintronic devices at room temperature (RT) such as next-generation spintronic biosensors, enabling enhanced sensitivity, reduced signal noise, and cost-effective detection of cancer biomarkers in clinical diagnostics.21 Recent studies have explored the diagnostic potential of volatile organic compounds (VOCs) as non-invasive biomarkers for cancer detection. Emerging evidence suggests that specific VOC signatures can reliably distinguish breast cancer patients from healthy individuals, reflecting alterations in metabolic pathways associated with tumorigenesis.22,23 To detect these disease-specific VOCs, a range of sensing technologies has been developed, including electronic noses, gas chromatography mass spectrometry (GC-MS), and nanomaterial-based sensor platforms.24–26 These approaches have demonstrated improved sensitivity and specificity in identifying VOC profiles associated with cancer, thereby offering promising avenues for early, non-invasive diagnostics. Despite the promise of VOC analysis in cancer diagnostics, significant challenges remain in the standardization of breath sampling protocols and the precise identification of disease-specific VOC biomarkers.27 In parallel, advancements in spintronic biosensing devices have further enhanced detection capabilities. Conventional single-barrier magnetic tunnel junctions (SB-MTJs), typically fabricated using Fe/MgO/Fe configurations, have demonstrated notable TMR effects and compatibility with established nanofabrication techniques.28 However, double-barrier magnetic tunnel junctions (DB-MTJs) have recently emerged as promising architectures to surpass the limitations of SB-MTJs. These multilayered structures have shown substantial improvements in TMR performance, with reported increases from 122% to 802% in Fe/MgO/Fe-based systems.29 The enhancement is primarily attributed to quantum well states and resonant tunneling phenomena that arise within the DB-MTJ structure.30,31 Furthermore, DB-MTJs exhibit increased sensitivity to interfacial modifications, making them particularly well-suited for biosensing applications where detection of biomolecular interactions at the surface is critical. Notably, the development of iron-free DB-MTJs may also offer additional benefits in biocompatibility and device stability. All these achievements open up avenues for future spintronic device developments, including prospective applications in biosensors.32
We explore in this article the possibility of DB-MTJs featuring transition metal dichalcogenides (TMDCs) like MoS2, and MoTe2 as composite tunnel barriers coupled with Co2MnSi Heusler alloy electrodes for biosensing applications of high sensitivity. By applying density functional theory (DFT) and nonequilibrium Green's function (NEGF) formalisms, we compute important properties like TMR ratios, transmission spectra, density of states (DOS), and bandgaps for Co2MnSi–MgO–MoTe2–MgO–Co2MnSi devices. The obtained TMR ratios is as large as 1226%. These TMR values coupled with the electronically tunable features of TMDCs make DB-MTJs a promising device to identify cancer biomarkers and VOCs in different forms of cancer. The iron-free nature helps to increase longevity through reduced degradation compared to classical Fe-based MTJs, and the increased sensitivity to interfacial variations allows for label-free, non-invasive sensing. The research opens ways to create sophisticated biosensors for the diagnosis of cancer at earlier stages, a versatile, high-performance tool to enhance healthcare and environmental monitoring.
![]() | (1) |
DB-MTJs offer an advanced architecture to enhance the TMR effect compared to single-barrier MTJs (SB-MTJs). Schematic representations of SB-MTJ and DB-MTJ structures are shown in Fig. 2(a and b) respectively. In the DB-MTJ configuration, two insulating barriers are separated by an intermediate ferromagnetic layer, which may be composed of either the same or different material as the outer ferromagnetic leads. In SB-MTJs, the P and AP magnetic configurations are established by maintaining the magnetization of the left electrode while reversing the magnetization of the right electrode. In contrast, for DB-MTJs, the P and AP states are realized by fixing the magnetization directions of the two outer electrodes and inverting the magnetization of the central ferromagnetic layer. As a result, in the AP configuration of a DB-MTJ, electrons experience spin misalignment at two interfaces during transport, whereas in the P configuration, the spin alignment is preserved throughout the structure. The transmission probability, which quantifies the likelihood of an electron tunneling from the left to the right electrode, is typically higher in the P configuration than in the AP configuration. This leads to a higher TMR ratio in DB-MTJs compared to SB-MTJs. The tunneling process in DB-MTJs is effectively composed of two sequential and independent tunneling events through each barrier. This mechanism is often referred to as sequential independent tunneling, and it significantly enhances the spin-filtering effect, contributing to the superior TMR performance of DB-MTJs.29
![]() | ||
| Fig. 2 Schematic of (a and b) single-barrier and (c and d) double-barrier magnetic tunnel junctions (MTJs) with corresponding spin-dependent energy band diagrams. | ||
The suggested DB-MTJ has a symmetric layered configuration, FM/barrier/spacer/barrier/FM (Fig. 1(b)). Wherein the FM electrodes consist of iron-free Heusler alloys (Co2MnSi), and the MgO as a tunnel barriers comprise insulators merged with semiconducting 2D materials MoS2, or MoTe2. These DB-MTJ devices in biosensing applications can be coupled with a microfluidic system for the delivery of the analyte toward the sensor surface. As illustrated in Fig. 3, biomarkers binding events linked with specific ligands or antibodies cause local variations in the potential, orientation of dipoles, or the magnetic field, resulting in spin-dependent transport variations in the junction, and this may be measured quantitatively through conductance and TMR variations.
All calculations in this work were performed using DFT employing the linear combination of atomic orbitals (LCAO) method as in Quantum ATK software (version 2023.12). The exchange–correlation potentials were calculated through the generalized gradient approximation (GGA) and the Perdew–Burke–Ernzerhof (PBE) functional.29 Structural relaxation was attained using geometry optimizations utilizing the limited memory Broyden–Fletcher–Goldfarb–Shanno (LBFGS) method, with a force tolerance of 10−4 eV Å−1 and a stress tolerance of 0.1 GPa. Transport characteristics were computed using a mixture of DFT and the non-equilibrium Green's function (NEGF) approach, according to the Landauer–Büttiker formalism. Consequently, the MTJ structures were thoroughly optimized at the interfaces between the barrier and electrode, achieving a total energy convergence of less than 10−5 eV, with the maximum force on each atom remaining below 0.05 eV Å−1. The optimization was conducted via the LBFGS optimizer approach and SGGA exchange–correlation. The k-point grid was established at 7 × 7 × 100, and a linear combination of atomic orbitals (LCAO) utilizing single zeta polarization was chosen as the basis set (see Fig. S1 for more details). To calculate band structure, we incorporated GGA-PBE and DFT-1/2 method. Here, 140 hartree was adjusted for cut-off energy of the density mesh with the Monkhorst–pack k-point of 13 × 13 × 13.
The physical parameters employed for simulations of parallel and anti-parallel configuration simulations are presented in Table 1. The energy convergence criterion for transmission spectrum computations was established at 10−5 eV. The obtained transmission spectrum plots are used to calculate the transmission coefficients at the level of the Fermi energy.
| Physical parameters | |
|---|---|
| Broadening temperature | 1200 K |
| Spin | Polarized |
| Exchange correlation | SGGA (PBE functional) |
| Density mesh cut-off | 140 h |
| Pseudopotential | SG15 |
| k-Point sampling (A, B, C) | (7, 7, 100) |
Transmission spectra were computed for both parallel and antiparallel magnetic configurations and showed spin polarized transport behaviour and determined the most active tunnelling states. Moreover, the TMR ratio was calculated at the Fermi level. These findings provide the basis for the examination of molecular interaction on the sensor interfaces in modulating spin current transmission.
In addition, DOS and bandgap calculations were also performed for MoS2, MoTe2, and Co2MnGe materials. These calculations give insight into the level of spin polarization, the transport channel nature, and the band alignment, which can promote or suppress electron tunnelling. These electronic properties are equally important because they determine how biomarker attachment make changes in dipole orientation, strain, or electrostatic potential, and their effect on the tunnelling probabilities and TMR signal.
This simulation platform develops a firm base for designing and optimizing TMR biosensors through leveraging the increased sensitivity of DB-MTJs toward interface variations. These biosensors will identify biomarkers by quantifiable variations in electronic conductance, spin transport, and pave the way for the development of low-power diagnostic devices for medical and environmental applications.
Specifically, the transmission coefficients at the Fermi level for both P and AP states were extracted from the k-resolved transmission spectra (Fig. 4 and S2) and are summarized in Table 2. Notably, spin-dependent scattering at the interfaces between the FM layers and the insulating barriers significantly influences the tunnelling behaviour. These interfacial effects modulate the spin-polarized transmission pathways, resulting in distinct differences in electron transmission between P and AP configurations, which in turn underlie the observed variations in TMR performance.
| Device | Parallel configuration | Anti-parallel configuration | ||||
|---|---|---|---|---|---|---|
| Spin-up | Spin-down | Total | Spin-up | Spin-down | Total | |
| Co2MnSi–MgO–Co2MnGe–MgO–Co2MnSi | 3.72 × 10−5 | 1.49 × 10−20 | 3.72 × 10−5 | 4.85 × 10−6 | 1.49 × 10−20 | 4.86 × 10−6 |
| Co2MnSi–MgO–MoS2–MgO–Co2MnSi | 6.28 × 10−6 | 1.49 × 10−18 | 9.57 × 10−8 | 9.52 × 10−8 | 1.42 × 10−20 | 9.52 × 10−8 |
| Co2MnSi–MgO–MoTe2–MgO–Co2MnSi | 4.06 × 10−6 | 2.79 × 10−9 | 4.062 × 10−6 | 7.31 × 10−7 | 1.63 × 10−18 | 3.06 × 10−7 |
| Co2MnSi–MgO–Cu–MgO–Co2MnSi | 7.12 × 10−6 | 6.75 × 10−13 | 7.12 × 10−6 | 3.78 × 10−6 | 2.12 × 10−18 | 3.78 × 10−6 |
The calculated transmission values derived using eqn (2) are consistently higher for the P configurations compared to the AP ones, which aligns with the expected lower electrical resistance in the P alignment due to more efficient spin-aligned electron tunneling.
![]() | (2) |
Furthermore, Fig. 5 illustrates the electron transmission spectra at the Fermi level for Co2MnSi–MgO–Co2MnGe–MgO–Co2MnSi and Co2MnSi–MgO–MoTe2–MgO–Co2MnSi DB-MTJ structures under total, P, and AP magnetic configurations. The transmission profiles clearly indicate that the P configuration exhibits significantly higher electron transmission, particularly at the Fermi level compared to the AP configuration in both structures. These findings are consistent with the expected spin-dependent tunneling behavior, where the alignment of spin-polarized states in the P configuration facilitates more efficient electron transport. Quantitative values of electron transmission at the Fermi level for each configuration are summarized in Table 2.
The TMR values were calculated using the transmission coefficients obtained for both P and AP magnetic configurations via QuantumATK simulations. Eqn (1) was employed to determine the TMR ratio based on these spin-dependent transmission values. The resulting TMR ratios for the DB-MTJs investigated in this study are presented in Table 3. Notably, the Co2MnSi–MgO–MoTe2–MgO–Co2MnSi structure exhibited the highest TMR ratio of 1226.20%, significantly surpassing that of Co2MnSi–MgO–MoS2–MgO–Co2MnSi (TMR = 22%), and Co2MnSi–MgO–Co2MnGe–MgO–Co2MnSi (TMR = 667.43%). Additionally, as shown in Fig. 5, the transmission probability at the Fermi level in the P configuration exhibits a substantial increase, reaching 4.06 × 10−6 for spin-up electrons, in contrast to the AP configuration, which yields a significantly lower value of 7.31 × 10−7. This pronounced difference aligns with the high TMR ratio obtained in this study, further confirming the enhanced spin-polarized transport efficiency in the P state. These findings highlight the critical influence of the spacer material on the spin-dependent transport properties and overall device performance. The substantial variation in resistance between magnetic configurations, as captured by the TMR ratio, serves as a key metric for evaluating and comparing the efficiency of different DB-MTJ architectures (see Fig. S3 for more details of script calculates the TMR).
| Materials | TMR% |
|---|---|
| Co2MnSi–MgO–Co2MnGe–MgO–Co2MnSi | 667.43 |
| Co2MnSi–MgO–MoS2–MgO–Co2MnSi | 22 |
| Co2MnSi–MgO–MoTe2–MgO–Co2MnSi | 1226.20 |
| Co2MnSi–MgO–Cu–MgO–Co2MnSi | 421 |
The integration of MoTe2 as a two-dimensional spacer material within the double tunnel barrier structure (MgO–MoTe2–MgO) appears to enhance spin-dependent tunneling, attributed to its optimal bandgap (∼1.1 eV) and intrinsic spin-filtering capabilities, as previously demonstrated for other transition metal dichalcogenides (TMDCs).33 In comparison, the Co2MnSi–MgO–Cu–MgO–Co2MnSi structure yielded a lower TMR ratio of 412%, reinforcing the superior performance of TMDC-based barriers in achieving high TMR values. Notably, a record TMR ratio of 4220.59% was reported by Sinha et al. for a Fe–MgO–MoTe2–MgO–Fe-based DB-MTJ, underscoring the beneficial role of MoTe2 in enhancing spin transport due to its wide bandgap.33 However, the use of iron-free Co2MnSi electrodes in the present study offers distinct advantages, including mitigation of oxidation-related degradation commonly associated with Fe-based systems, thereby improving the long-term stability and reliability of MTJ devices for biosensing applications.
This trend is consistent with the findings reported by Zheng et al., in which Fe/MgO/Fe DB-MTJs exhibited a comparable disparity in electron transmission between the P and AP states, resulting in a TMR ratio of 802%.29 The enhanced transmission observed in the P configuration can be attributed to the sequential independent tunneling mechanism across the two insulating barriers, which reinforces spin filtering and significantly amplifies the spin-dependent transport characteristics. This mechanism underscores the superiority of DB-MTJs over conventional single-barrier MTJs in achieving higher TMR values and improved sensitivity for spintronic applications. The comparison in TMR ratio of the suggested material structures and reference structures are shown in Table 4.
| Structure (this work) | Reference structure | TMR (%) | Ref. |
|---|---|---|---|
| Co2MnSi–MgO–MoTe2–MgO–Co2MnSi (TMR% = 1226.20) | Fe–MgO–MoTe2–MgO–Fe | 4220.59 | 33 |
| Fe–MgO–MoSe2–MgO–CoFeB | 2000 | 33 | |
| CoFe–MgO–MoS2–MgO–CoFe | 1300 | 33 | |
| Co2MnSi–MgO–Cu–MgO–Co2MnSi (TMR% = 421) | — | — | — |
| Co2MnSi–MgO–Co2MnGe–MgO–Co2MnSi (TMR% = 667.43) | Co2TiGe/MgO/Co2TiGe | 645 | 34 |
| Co2FeAl–2L–MgO–Co2FeAl | 373.27 | 19 | |
| Co2TiSi/MgO/Co2TiSi | 1894 | 34 |
Strong hybridization is observed between the d states of the transition metals and the s states of the chalcogen atoms below the Fermi energy level in all structures. Furthermore, when comparing Mo-based chalcogenides from sulfur (S) to tellurium (Te) (Fig. 6(a and b)), a reduction in the bandgap is noted. This trend is attributed to the increasing atomic radii and decreasing electronegativity from S to Te. The electronic states near the Fermi level are primarily derived from Mo d orbitals, and as the ligand field strength of the chalcogen increases (moving from Te to S), broader d-bands are formed, leading to an increased bandgap.
The electrical structure significantly influences the magnetic characteristics of crystalline materials. Consequently, we determined the electronic band structure of the stable phase of the Heusler compound Co2MnGe utilising polarised spin computations inside the GGA-PBE (Table S1) and DFT-1/2 method. Fig. 6(c) displays the computed dispersion curves along the high symmetry directions of the first Brillouin zone. This figure illustrates that the majority-spin band exhibits overlapping valence and conduction bands, resulting in metallic behaviour, whereas the minority-spin band demonstrates semiconductor characteristics, with the valence band maximum and conduction band minimum located at the Γ and X points, respectively.35 Fig. 6(d) presents the calculated electronic band structure and calculated DOS for the structure of MgO. The Fermi level located in a broad bandgap, with significant any electronic states immediately near the Fermi energy was observed. This confirms the highly insulating nature of MgO, consistent with its experimentally known bandgap energy of about 6.7 eV. The DOS plot (Fig. 6(d)) illustrates that the occupied states (valence band) and unoccupied states (conduction band) possess a large bandgap, with negligible hybridization close to the Fermi level. The valence band is constituted mainly of oxygen 2p orbitals, whereas the conduction band comprises mainly magnesium 3s orbitals. This bandgap results (Table 5) in the large dielectric breakdown strength and stability of MgO and makes it a favorable material for application as an insulating barrier in MTJs. On the other hand, the Heusler alloy spacer, Co2MnGe (Fig. 6(c)) is a ferromagnet exhibiting a spin-down bandgap of around ∼0 eV.18 This characteristic guarantees substantial spin polarisation, enhancing the favourable TMR ratio of 667%. The findings underscore the pivotal influence of the spacer material's bandgap on tunnelling efficiency and TMR, with reduced bandgaps often promoting elevated TMR ratios in DB-MTJ systems. Moreover, TMR is not controlled by the barrier bandgap itself, yet greatly tuned by interfacial resonance states, e.g., in (111)-oriented Co/MgO/Co junctions, high TMR ratios are induced by the resonance enhancement of tunneling through one spin channel by interfacial antibonding states which are created close to the Fermi level.36
| Structure | Energy bandgap |
|---|---|
| MoS2 | 1.78 eV |
| MoTe2 | 1.11 eV |
| MgO | 6.7 eV |
| Co2MnGe | ∼0 eV |
Supplementary information: simulation parameters and output data. See DOI: https://doi.org/10.1039/d5ra04973c.
| This journal is © The Royal Society of Chemistry 2025 |