AV

A. Vasileiadis

info

Please Note

34 records found

Journal article (2026) - Tim Horner, Enis Oğuzhan Eren, Elif Begüm Yılmaz, Jiyong Kim, Ernesto Scoppola, Alexandros Vasileiadis, Nadezda V. Tarakina, Markus Antonietti, Evgeny Senokos, More authors...
Understanding sulfur confinement and chemical transformation in hybrid sulfur-carbon materials is critical for advancing metal-sulfur batteries. Here, we investigate the structural evolution of a sulfur-rich polymer into a hybrid sulfur-carbon via inverse vulcanization and thermal condensation. Multiscale analyses reveal a stepwise transformation, beginning with the emergence of sulfur radicals at ∼175°C, followed by the progressive development of a carbon matrix above 300°C that stabilizes the radical species. Around 450°C, a transitional phase forms, consisting of conjugated carbon clusters covalently bonded to sulfur chains. This hybrid structure confines sulfur within pseudo-graphitic nanodomains, effectively suppressing polysulfide dissolution and enhancing redox stability. DFT simulations show how sulfur confinement modulates Na-S reaction energetics, while electrochemical testing confirms high sulfur utilization, delivering ∼1000 mAh (Formula presented.) and 1200 Wh (Formula presented.), setting a new performance benchmark for room-temperature Na─S batteries. These findings provide critical insights into the correlation between structural evolution and electrochemical performance, offering design principles for next-generation sulfur-based electrodes. ...

A Python toolkit for site-resolved diffusion analysis in solid-state molecular dynamics

Journal article (2026) - Anastasia K. Lavrinenko, Theodosios Famprikis, Victor Landgraf, Jouke R. Heringa, Stef Smeets, Victor Azizi, Simone Ciarella, Marnix Wagemaker, Alexandros Vasileiadis
Molecular dynamics (MD) simulations have become essential for understanding diffusion mechanisms in solid-state materials such as ionic conductors, fuel cells, and gas sensors, yet most existing studies and software tools extract only standard metrics, leaving much of the information contained in the trajectories unused. Here we introduce GEMDAT, a user-friendly Python toolkit for site-resolved diffusion analysis of MD simulations of solid-state materials (https://github.com/GEMDAT-repos/GEMDAT). Beyond mean-squared displacements, radial distribution functions, and Arrhenius-based activation energies, GEMDAT provides jump rates, attempt frequencies, site-specific activation energies, rotational diffusion, and percolation. Our tool provides access to vibrational amplitudes, site geometries, and site occupancies—quantities that are also directly comparable to experimental diffraction data. Migration sites can be defined manually or identified automatically from the trajectory. A built-in caching approach, together with rapid visualization capabilities, makes the workflow fast and interactive. We demonstrate GEMDAT on a series of case studies involving crystalline Li- and Na-ion conductors, plastic crystals, amorphous structures, and surface configurations, showing how the code extracts atomic-level structural features and connects them to macroscopic transport properties, thereby guiding the optimization and development of solid-state materials. ...
Journal article (2026) - Jeffrey Zom, Alexandros Vasileiadis, Shuaishuai Yuan, Erik van der Kolk
Photoluminescence in vanadate compounds has traditionally been attributed to charge-transfer transitions within isolated [VO4]3− centres, with excitation and emission described using molecular orbitals. In previous work, we proposed an alternative mechanism: photoluminescence arises from conventional interband excitation, followed by electron polaron-mediated formation of self-trapped excitons. In this study, we provide further evidence for this model through spectroscopic measurements and ab initio calculations. Using density functional theory calculations on a series of alkali vanadates (MVO3, M = Li, Na, K, Rb, Cs), we show that photoexcited electrons spontaneously localise on V5+ ions, leading to self-trapped exciton formation via hole localisation on neighbouring oxygen atoms. The calculated energies for band gaps and self-trapped excitons closely match experimental values obtained from diffuse reflectance and luminescence spectroscopy. Importantly, temperature- and time-resolved luminescence measurements reveal that quenching predominantly occurs before the formation of the luminescent state, challenging earlier models that assumed quenching to occur from the final emitting state. To explain this behaviour, we note that the width of the conduction band states is expected to govern the rate of electron trapping, and we indeed find a correlation between calculated bandwidths and measured quantum efficiencies. This suggests that non-radiative relaxation of free carriers at defects, prior to self-trapping, is the dominant quenching mechanism. Consequently, the electron self-trapping rate, the self-trapped exciton formation rate, and the defect concentration are expected to critically determine the luminescent efficiency of vanadate phosphors. ...
Silicon (Si) anodes promise high energy density for all-solid-state batteries (ASSBs), while avoiding Li-metal dendrites. However, one of the key limitations is the initial loss of active Li, reflected in a low initial Coulombic efficiency (ICE), arising from solid electrolyte (SE) decomposition, lithium trapping, and native oxide conversion. Here, we report an (electro)chemically guided interfacial-engineering strategy leveraging an irreducible SE to construct a dynamic nanoscale Si|SE interface that enhances ICE. The intrinsic stability of the SE down to the working potentials of Si anodes eliminates decomposition, whereas mechanochemical mixing activates interfacial reactions that realize in situ prelithiation. The engineered interface manifests efficient ion transport and replenishes Li inventory through redox reactions, while stabilizing electrochemical performance. Our approach delivers a record-high ICE approaching 100% in Si half-cells, together with >95% in high-loading LiCoO₂ full cells, demonstrating that precise interface control can unlock the full potential of Si anodes for high-energy ASSBs. ...
Journal article (2026) - A. Gautam, Ruihua Zhou, A.K. Lavrinenko, H.A.A. Al-Kutubi, S. Ganapathy, Xin Zhang, A. Vasileiadis, Shuo Wang, M. Wagemaker
Lithium argyrodite solid electrolytes are promising for all-solid-state batteries due to their high ionic conductivity and low elastic modulus. However, poor chemical stability in humid conditions remains a major hurdle towards large-scale practical implementation. This motivates substituting phosphorus with softer acids to stabilize sulfur, where the challenge is to maintain high Li-ion mobility. In this work, we explore the effects of aliovalent (Ge) substitutions on the structure, ionic transport, and air stability of Li6AsS5Br. The induced structural modification releases the rate-limiting step in long-range Li-ion transport, as demonstrated by Molecular Dynamics and percolation simulations. As a result, the ionic conductivity reaches 13 mS/cm for Li6.5As0.5Ge0.5S5Br, a 6-fold improvement over pristine Li6AsS5Br. The key enabler for combining high Li-ion mobility and improved stability towards moisture is the increased polarity of the Ge─S bonds compared to As─S or P─S bonds. This reduces Li-ion trapping effectively flattening the energy landscape for diffusion and thermodynamically suppresses H2S formation upon exposure to moisture, as experimentally demonstrated. Finally, stable cycling performances in combination with high nickel ternary cathodes are demonstrated. Hereby, this research provides a deeper understanding of the role of composition on the ionic conductivity and moisture stability of lithium argyrodites. ...
Review (2026) - Anastasia K. Lavrinenko, James A. Quirk, James A. Dawson, Marnix Wagemaker, Alexandros Vasileiadis
Solid–solid interfaces play a critical role in the performance of energy storage systems, catalysts, and electronic devices, yet their atomic-scale structure and behaviour remain challenging to probe experimentally. Explicit atomistic modelling provides a powerful approach for investigating interfaces by directly constructing and simulating contacting phases, enabling the calculation of interfacial energetics, electronic structure, charge transfer, stability, and transport processes. This Review summarises the methodology and applications of explicit interface modelling, including interface construction strategies, properties accessible from atomistic simulations, and recent advances enabled by machine-learned interatomic potentials that extend simulations to larger time and length scales. We also discuss how computational results can be interpreted and validated using experimental characterisation techniques, and highlight current challenges and emerging opportunities in predictive interface modelling and materials design. ...
Journal article (2026) - Gangbin Yan, Pierfrancesco Ombrini, Zhichu Tang, Shakul Pathak, Maoyu Wang, Barbara Lavina, Alexandros Vasileiadis, Marnix Wagemaker, Chong Liu, More Authors
Ionic diffusion in solids underpins energy storage, electronics, and catalysis, yet conventional diffusion models often fail to capture complexities arising from confinement, crystallographic disorder, lattice distortions, and coupled transport with phonons or electrons. These challenges are particularly pronounced in battery materials, where ionic and electronic carriers move together, complicating the interpretation of electrochemical measurements. Here we employ tracer exchange as a direct, non-electrochemical probe to reveal rich ion dynamics in the model one-dimensional (1D) conductor olivine LiXFePO4 (0 ≤ X ≤ 1). 6Li-7Li isotope exchange confirms single-file diffusion (SFD), where 1D confinement prevents ion bypassing and preserves spatial order. Kinetic Monte Carlo (KMC) simulations and chronoamperometry further quantify Faradaic and non-Faradaic surface exchange, identifying electron transport as rate-limiting during electrochemical reactions. In contrast, Li-Na exchange exhibits apparent superdiffusion, where the exchange rate increases with Na content. Simulations attribute this behavior to surface-exchange limitation and Na+-enhanced Li+ cross-channel hopping that drives a dimensional crossover from 1D to quasi-2D transport, supported by 4D-STEM and in situ synchrotron XRD. These results establish tracer exchange as a powerful platform for probing coupled multi-ion and electron transport in solids. ...
Phase separation, inducing a miscibility gap and non-monotonic open-circuit potential (OCP), is typical for widespread Li-ion battery electrodes such as LiFePO4, Li4Ti5OPhase separation, inducing a miscibility gap and non-monotonic open-circuit potential (OCP), is typical for widespread Li-ion battery electrodes such as LiFePO4, Li4Ti5O12 and Graphite. Although particle-scale effects of phase separation are well documented, its influence on transport-limited, porous electrodes remains largely overlooked. Here we embed physically consistent non-monotonic OCP profiles in a simplified Doyle–Fuller–Newman framework to compare their behavior against that of solid- solution materials with monotonic OCPs. Our findings provide deeper and general understanding of the different electrode ensemble behavior of solid solution (monotonic OCP) and phase separating (non-monotonic OCP) electrode materials, demonstrating why larger miscibility gaps are associated with decreasing rate capabilities and electrode utilization, amplifying local current heterogeneity and electrolyte depletion. By contrast, simulations employing conventional flat, fitted OCPs mask these effects and overpredict performance—particularly under dynamic cycling protocols such as galvanostatic intermittent titration (GITT). Our results reveal why accounting for realistic OCPs is essential for reliable modelling of high-loading electrodes, providing fundamental understanding and guidance for model-driven design and control of next-generation batteries and Graphite. Although particle-scale effects of phase separation are well documented, its influence on transport-limited, porous electrodes remains largely overlooked. Here we embed physically consistent non-monotonic OCP profiles in a simplified Doyle–Fuller–Newman framework to compare their behavior against that of solid- solution materials with monotonic OCPs. Our findings provide deeper and general understanding of the different electrode ensemble behavior of solid solution (monotonic OCP) and phase separating (non-monotonic OCP) electrode materials, demonstrating why larger miscibility gaps are associated with decreasing rate capabilities and electrode utilization, amplifying local current heterogeneity and electrolyte depletion. By contrast, simulations employing conventional flat, fitted OCPs mask these effects and overpredict performance—particularly under dynamic cycling protocols such as galvanostatic intermittent titration (GITT). Our results reveal why accounting for realistic OCPs is essential for reliable modelling of high-loading electrodes, providing fundamental understanding and guidance for model-driven design and control of next-generation batteries.. ...
All-solid-state batteries receive ample attention due to their promising safety characteristics and energy density. The latter holds true if they are compatible with next-generation high-capacity anodes, but most highly ion-conductive solid electrolytes decompose at low operating potentials, leading to lithium loss and increased cell resistances. Here the dynamic stability of solid electrolytes that can improve all-solid-state battery performance is demonstrated. Halide electrolytes Li3YCl3Br3 and Li2ZrCl6, considered unstable at low potentials, are found to exhibit structurally reversible redox activity beyond their electrochemical stability windows, increasing compatibility with anodes and contributing to capacity without compromising ionic conductivity. The benefit of this dynamic stability window is demonstrated with cost-effective red phosphorus anodes, resulting in high reversible capacities (2,308 mAh g−1), high rate capacity retention (1,024 mAh g−1 at 7.75 mA cm−2) and extended cycle life (61% retention after 1,780 cycles). Furthermore, high areal capacity (7.65 mAh cm−2) and stability (70% retention after 1,000 cycles) are achieved for halide-based full cells with red phosphorous anodes. The beneficial redox activity of halide electrolytes greatly expands their application scenarios and suggests valuable battery design principles to enhance performance. ...
Solid-state batteries currently receive extensive attention due to their potential to outperform lithium-ion batteries in terms of energy density when featuring next-generation anodes such as lithium metal or silicon. However, most highly conducting solid electrolytes decompose at the low operating voltages of next-generation anodes leading to irreversible lithium loss and increased cell resistance. Such performance losses may be prevented by designing electrolytes which are thermodynamically stable at low operating voltages (anolytes). Here, we report on the discovery of a new family of irreducible (i.e., fully reduced) electrolytes by mechanochemically dissolving lithium nitride into the Li2S antifluorite structure, yielding highly conducting crystalline Li2+xS1-xNx phases reaching >0.2 mS cm-1 at ambient temperature. Combining impedance spectroscopy experiments and ab initio density functional theory calculations we clarify the mechanism by which the disordering of the sulfide and nitride ions in the anion sublattice boosts ionic conductivity in Li2+xS1-xNx phases by a factor 105 compared to the Li2S host structure. This advance is achieved through a novel theoretical framework, leveraging percolation analysis with local-environment-specific activation energies and is widely applicable to disordered ion conductors. The same methodology allows us to rationalize how increasing nitrogen content in Li2+xS1-xNx antifluorite-like samples leads to both increased ionic conductivity and lower conductivity-activation energy. These findings pave the way to understanding disordered solid electrolytes and eliminating decomposition-induced performance losses on the anode side in solid-state batteries. ...
By varying the bromine content and cooling method, we are able to induce site disorder in the Li6-xPS5-xBr1+x (x = 0, 0.3, 0.5) system via two routes, allowing us to disentangle the impact of site disorder and chemical composition on conductivity. Through solid-state nuclear magnetic resonance (NMR), we can explore the chemical environment as well as short-range lithium-ion dynamics and compare these to results obtained from neutron diffraction and electrochemical impedance spectroscopy (EIS). We find that the cooling method has a profound effect on the 7Li and 31P environment that cannot be explained through 4d site disorder alone. The configurational entropy (Sconf) is used as a more complete descriptor of structural disorder and linked to distortions in both the phosphorus and lithium environment. These distortions are correlated to increased intercage movement through 7Li T1 spin-lattice relaxation (SLR) NMR. Further analysis of the prefactors obtained from SLR NMR and EIS allows us to obtain the migrational entropy (ΔSm). For short-range SLR movement, the ΔSm correlates well with Sconf, implying that increased intercage movement is related to distortion of the lithium cages as well as a decrease of the intercage distance. Comparison to EIS shows that an increase in short-range movement translates into increased long-range movement in a straightforward manner for slow-cooled samples. However, for quench-cooled samples, this correlation is lost. Lattice softness and phonon-ion interactions are suggested to play an important role in long-range conduction which only becomes apparent when chemical composition and disorder are disentangled. This work shows that by altering one synthesis step, the relationship between site-occupancy-based descriptors (site disorder or Sconf) and lithium dynamics is changed profoundly. Furthermore, it shows that chemical composition and descriptors of site disorder cannot be seen as one and the same, as both play a role that changes with the length scale probed. Finally, it challenges the implicit assumption that increased short-range diffusivity automatically results in increased long-range diffusivity. ...

New Perspectives on Stabilizing High-Capacity Anodes in Solid-State Batteries

Irreducible solid electrolytes (SEs), characterized by non-Li framework ions in their lowest oxidation states, offer intrinsic compatibility with low-reduction-potential, high-capacity negative electrodes, such as lithium metal and silicon. In these SE materials, disorder engineering and vacancy formation reduce lithium-ion diffusion barriers, achieving room-temperature ionic conductivities exceeding 0.1 mS cm–1. Experiments and atomistic simulations confirm that irreducible SEs form decomposition-free interfaces with Li metal. Their limited oxidative stability can be addressed by pairing them with an electrolyte layer stable with practical cathodes yet demanding interface compatibility between the two electrolyte layers. Here we highlight key research directions to accelerate irreducible SE transition from laboratory to practical application, including expanding compositional diversity, optimizing interfaces with cathode-facing electrolytes, developing scalable thin-film processing, and exploring compatibility with other low working potential anodes like silicon. Addressing these challenges is essential to unlock the full potential of irreducible SEs for high-energy, long-life, all-solid-state batteries. ...

Origin of fast charging in hard carbon anodes (Nature Energy, (2024), 9, 2, (134-142), 10.1038/s41560-023-01414-5)

Correction to: Nature Energyhttps://doi.org/10.1038/s41560-023-01414-5, published online 3 January 2024. In the version of this article initially published, lithium (green, “Li”) and sodium (purple, “Na”) color key labels in Fig. 3a,d,e were interchanged and are now amended in the HTML and PDF versions of the article. ...
Journal article (2024) - Evgeny Senokos, Heather Au, Enis Oğuzhan Eren, Tim Horner, Zihan Song, Nadezda V. Tarakina, Elif Begüm Yılmaz, Alexandros Vasileiadis, Paolo Giusto, More authors...
Nanoconfinement is a promising strategy in chemistry enabling increased reaction rates, enhanced selectivity, and stabilized reactive species. Sulfur's abundance and highly reversible two-electron transfer mechanism have fueled research on sulfur-based electrochemical energy storage. However, the formation of soluble polysulfides, poor reaction kinetics, and low sulfur utilization are current bottlenecks for broader practical application. Herein, a novel strategy is proposed to confine sulfur species in a nanostructured hybrid sulfur-carbon material. A microporous sulfur-rich carbon is produced from sustainable natural precursors via inverse vulcanization and condensation. The material exhibits a unique structure with sulfur anchored to the conductive carbon matrix and physically confined in ultra-micropores. The structure promotes Na+ ion transport through micropores and electron transport through the carbon matrix, while effectively immobilizing sulfur species in the nanoconfined environment, fostering a quasi-solid-state redox reaction with sodium. This translates to ≈99% utilization of the 2e− reduction of sulfur and the highest reported capacity for a room temperature Na−S electrochemical system, with high rate capability, coulombic efficiency, and long-term stability. This study offers an innovative approach toward understanding the key physicochemical properties of sulfurcarbon nanohybrid materials, enabling the development of high-performance cathode materials for room-temperature Na-S batteries with efficient sulfur utilization. ...

A unified view on the role of sulfur/halide distribution and local environments

Understanding diffusion mechanisms in solid electrolytes is crucial for advancing solid-state battery technologies. This study investigates the role of structural disorder in Li7−xPS6−xBrx argyrodites using ab initio molecular dynamics, focusing on the correlation between key structural descriptors and Li-ion conductivity. Commonly suggested parameters, such as configurational entropy, bromide site occupancy, and bromine content, correlate with Li-ion diffusivity but do not consistently explain conductivity trends. We find that a uniform distribution of bromine and sulfur ions across the 4a and 4d sublattices is critical for achieving high conductivity by facilitating optimal lithium jump activation energies, anion-lithium distances, and charge distribution. Additionally, we introduce the ionic potential as a simple descriptor that predicts argyrodite conductivity by assessing the interaction strength between cations and anions. By analyzing the correlation between ionic potential and conductivity for a range of argyrodite compositions published over the past decade, we demonstrate its broad applicability. Minimizing and equalizing ionic potentials across both sublattices enhances conductivity by reducing the strength of anion-lithium interactions. Our analysis of local environments coordinating Li jumps reveals that balancing high and low-energy pathways is crucial for enabling macroscopic diffusion, supported by investigating percolating pathways. This study highlights the significance of the anionic framework in lithium mobility and informs the design of solid electrolytes for improved energy storage systems. ...
Due to their high ionic conductivity, lithium-ion conducting argyrodites show promise as solid electrolytes for solid-state batteries. Aliovalent substitution is an effective technique to enhance the transport properties of Li6PS5Br, where aliovalent Si substitution triples ionic conductivity. However, the origin of this experimentally observed increase is not fully understood. Our density functional theory (DFT) study reveals that Si4+ substitution increases Li diffusion by activating Li occupancy in the T4 sites. Redistribution of Li-ions within the lattice results in a more uniform distribution of Li around the T4 and neighboring T5 sites, flattening the energy landscape for diffusion. Since the T4 site is positioned in the intercage jump pathway, an increase in the intercage jump rate is found, which is directly related to the macroscopic diffusion and bulk conductivity. Analysis of neutron diffraction experiments confirms partial T4 site occupancy, in agreement with the computational findings. Understanding the aliovalent substitution effect on interstitials is crucial for improving solid electrolyte ionic conductivity and advancing solid-state battery performance. ...
Solid-state batteries currently receive ample attention due to their potential to outperform lithium-ion batteries in terms of energy density when featuring next-generation anodes such as lithium metal or silicon. One key remaining challenge is identifying solid electrolytes that combine high ionic conductivity with stability in contact with the highly reducing potentials of next-generation anodes. Fully reduced electrolytes, based on irreducible anions, offer a promising solution by avoiding electrolyte decomposition altogether. In this study, we demonstrate the compositional flexibility of the disordered antifluorite framework accessible by mechanochemical synthesis and leverage it to discover irreducible electrolytes with high ionic conductivities. We show that the recently investigated Li 9N 2Cl 3 and Li 5NCl 2 phases are part of the same solid solution of Li-deficient antifluorite phases existing on the LiCl-Li 3N tie line with a general chemical formula of Li 1+2xCl 1−xN x (0.33 < x < 0.5). Using density functional theory calculations, we identify the origin of the 5-order-of-magnitude conductivity increase of the Li 1+2xCl 1−xN x phases compared to the structurally related rock-salt LiCl phase. Finally, we demonstrate that S Cl- and Br Cl-substituted analogues of the Li 1+2xCl 1−xN x phases may be synthesized, enabling significant conductivity improvements by a factor of 10, reaching 0.2 mS cm −1 for Li 2.31S 0.41Br 0.14N 0.45. This investigation demonstrates for the first time that irreducible antifluorite-like phases are compositionally highly modifiable; this finding lays the ground for discovery of new compositions of irreducible antifluorite-like phases with even further increased conductivities, which could help eliminate solid-electrolyte decomposition and decomposition-induced Li losses on the anode side in high-performance next-generation batteries. ...
Journal article (2024) - Alexandros Vasileiadis, Quan Zhou, Yaxiang Lu, Yu Li, Pierfrancesco Ombrini, Zhao Chen, Remco van der Jagt, Swapna Ganapathy, Marnix Wagemaker, More authors...
Transport electrification and grid storage hinge largely on fast-charging capabilities of Li- and Na-ion batteries, but anodes such as graphite with plating issues drive the scientific focus towards anodes with slopped storage potentials. Here we report fast charging of ampere-hour-level full Na-ion batteries within about 9 minutes for continuous 3,000 cycles based on hard carbon anodes. These anodes, in addition to displaying a sloped storage voltage, provide capacity at a nearly constant voltage just above the plating potential, without observing Na-metal plating under high areal capacity. Comparing the electrochemical behaviour of Li and Na in hard carbon through experimental and computational techniques, a unified storage mechanism relying on the dimensions of wedge nanopores and drawing parallels with underpotential deposition for metals is brought forward, providing a rational guide for achieving fast storage in hard carbon anodes. ...
Journal article (2023) - Alexandros Vasileiadis, Yuqi Li, Yaxiang Lu, Yong Sheng Hu, Marnix Wagemaker
There are several questions and controversies regarding the Na storage mechanism in hard carbon. This springs from the difficulty of probing the vast diversity of possible configurational environments for Na storage, including surface and defect sites, edges, pores, and intercalation morphologies. In the effort to explain the observed voltage profile, typically existing of a voltage slope section and a low-voltage plateau, several experimental and computational studies have provided a variety of contradicting results. This work employs density functional theory to thoroughly examine Na storage in hard carbon in combination with electrochemical experiments. Our calculation scheme disentangles the possible interactions by evaluating the enthalpies of formation, shedding light on the storage mechanisms. Parallel evaluation of the Li and K storage, and comparison with experiments, put forward a unified reaction mechanism for the three alkali metals. The results underline the importance of exposed metal surfaces and metal-carbon interfaces for the stability of the pore-filling mechanism responsible for the low-voltage plateau, in excellent agreement with the experimental voltage profiles. This generalized understanding provides insights into hard carbons as negative electrodes and their optimized properties. ...

From an extended regular solution theory to a phase-field model of LiMnyFe1-yPO4

Journal article (2023) - Pierfrancesco Ombrini, Martin Z. Bazant, Marnix Wagemaker, Alexandros Vasileiadis
Phase separation during the lithiation of redox-active materials is a critical factor affecting battery performance, including energy density, charging rates, and cycle life. Accurate physical descriptions of these materials are necessary for understanding underlying lithiation mechanisms, performance limitations, and optimizing energy storage devices. This work presents an extended regular solution model that captures mutual interactions between sublattices of multi-sublattice battery materials, typically synthesized by metal substitution. We apply the model to phospho-olivine materials and demonstrate its quantitative accuracy in predicting the composition-dependent redox shift of the plateaus of LiMnyFe1-yPO4 (LFMP), LiCoyFe1-yPO4 (LFCP), LiCoxMnyFe1-x-yPO4 (LFMCP), as well as their phase separation behavior. Furthermore, we develop a phase-field model of LFMP that consistently matches experimental data and identifies LiMn0.4Fe0.6PO4 as a superior composition that favors a solid solution phase transition, making it ideal for high-power applications. ...