Phase Transitions under Magnetic Fields: Composition-Engineered Magnetocaloric Performance in Ni-Mn-Sn Heusler Alloys

Phase Transitions under Magnetic Fields: Composition-Engineered Magnetocaloric Performance in Ni-Mn-Sn Heusler Alloys | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article Phase Transitions under Magnetic Fields: Composition-Engineered Magnetocaloric Performance in Ni-Mn-Sn Heusler Alloys Md Mahmudul Karim¹, Sadiya Afrin Nowrin², Ismail Ahmed³ This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-9447379/v1 This work is licensed under a CC BY 4.0 License Status: Posted Version 1 posted You are reading this latest preprint version Abstract Understanding how magnetic fields alter phase transitions is a practical obstacle in magnetocaloric refrigeration, where optimizing the magnetocaloric effect requires precise knowledge of how a magnetic field shifts transition temperatures and latent heat. This problem matters because Ni-Mn-Sn Heusler alloys offer a rare-earth-free pathway to solid-state cooling near room temperature. To address this, we developed a hybrid molecular dynamics–Landau free energy framework with weak magnetostructural coupling (E = 0.08 eV/µB²) and applied it to Ni-Mn-Sn under fields from 0 to 20 T. The method captures both first-order structural (martensitic) and second-order magnetic transitions simultaneously. Key results show that at the optimal composition of 16 at.% Sn, a 5 T field increases the Curie temperature by 20 K (from 298 K to 318 K), suppresses thermal hysteresis by 40% (from 28 K to 16.8 K), yields a maximum entropy change of − 15 J/kg·K, and achieves a refrigerant capacity of 148 J/kg at 5 T — 24% higher than the widely studied 15 at.% Sn composition. The broader implication is that modest magnetic fields can be used to engineer optimal magnetocaloric performance in Ni-Mn-Sn through precise composition control, enabling compact, rare-earth-free magnetic heat pumps for near-room-temperature refrigeration. Magnetocaloric Effect Phase Transitions Ni-Mn-Sn Alloys Landau Theory Molecular Dynamics Magnetic Refrigeration 1 Introduction Phase transitions under external magnetic fields represent a fundamental tuning mechanism where the field serves as a clean parameter to drive systems across critical boundaries. For magnetocaloric materials, field-driven phase transitions are the physical basis for solid-state refrigeration, offering a promising alternative to conventional vapor-compression systems. In Ni-Mn-Sn Heusler alloys, experimental studies by Pasquale et al. (2007) [ 1 ] first reported a giant sign-switching magnetocaloric effect near the martensitic transition, with magnetization changes of up to 50 emu/g under 5 T. Dan et al. (2014) [ 2 ] systematically mapped how Sn concentration tunes the martensitic transition temperature and saturation magnetization, identifying a composition window where structural and magnetic transitions can be coupled. Computational work by Gai et al. (2013) [ 3 ] used DFT-derived exchange parameters in Monte Carlo simulations, but treated the lattice as static, neglecting the martensitic distortion that is essential for the first-order transition. However, recent work by Zhang et al. (2025) [ 4 ] raised a critical caveat: conventional resistivity-temperature curves can produce false positives for phase transitions in magnetic fields due to ordinary magnetoresistance scaling. What remains missing is a computational model that couples magnetic spins, martensitic lattice strain, and an external magnetic field on equal footing for Ni-Mn-Sn, enabling predictive optimization of composition for maximum entropy change. To fill this gap, we developed a parallel-tempering hybrid MD–Landau framework with weak magnetostructural coupling (E = 0.08 eV/µB²) for Ni-Mn-Sn, where MD simulates the temperature- and field-dependent martensitic lattice distortion while Landau theory provides the magnetic free energy including the Zeeman term. The objective is to compute the field-dependent entropy change ΔS(H) across the martensitic transition for Sn concentrations from 10 to 20 at.% and to identify the optimal composition for the largest magnetocaloric effect near 300 K. 2 Theoretical Framework Our hybrid framework treats Ni-Mn-Sn as a system governed by two coupled order parameters: the magnetic order parameter φ = M/M sat (normalized magnetization, 0 ≤ φ ≤ 1) and the structural order parameter ξ = (c/a − 1)/(c/a mart − 1) (martensitic fraction, 0 ≤ ξ ≤ 1) [ 6 ]. The total Landau free energy including the applied magnetic field is: F(φ, ξ, T, H) = (A/2)(T – Tᴄ)φ² + (B/4)φ⁴ + (C/2)(T – Tᴹ)ξ² + (D/4)ξ⁴ + (E/2)φ²ξ² − µ₀HMₛₐₜφ …… (1) where T M ≈ 310 K is the martensitic start temperature, A, B, C, D are positive material constants, and E is the biquadratic coupling coefficient. When E exceeds a critical value, simultaneous ordering of both parameters is penalized, producing a first-order transition. The applied field lowers the free energy of the ferromagnetic state (φ = 1), shifting the minimum in ξ through the coupling term and effectively suppressing the martensitic phase above a composition-dependent critical field H c ≈ 5–10 T. The thermodynamic treatment of field-driven first-order transitions follows established frameworks for magnetostructural coupling [ 5 ]. The ‘MD’ component employs a Finnis-Sinclair-type many-body potential optimized for the martensitic transformation in Ni-Mn-based Heusler alloys. Lattice distortion is tracked through the Bain strain ε B = 2(c − a)/(c + a), computed from variable-cell shape dynamics. Temperature is controlled via a Nosé-Hoover thermostat, and magnetic moments are represented as classical vectors µ i with magnitude µ B (the average Ni/Mn moment), with the external field applied as a uniform force on these moments. The two components are coupled through an iterative reweighting scheme. Every 100 ‘MD’ steps, the instantaneous distribution of ξ from atomic positions updates an effective coupling parameter E eff = E₀ + η⟨ξ⟩, which in turn refreshes the Landau free energy landscape. Conversely, the Landau-derived equilibrium magnetization φ eq (H, T) imposes a uniaxial anisotropy field H ani = (2K u /µ₀M sat ) sin(2θ) on each magnetic moment in ‘MD’, where K u ≈ 0.5 MJ/m³ is the magnetocrystalline anisotropy energy and θ is the angle between the moment and the c-axis. The entropy change ΔS is extracted using the Clausius-Clapeyron relation adapted for field-driven first-order transitions: ΔS = −µ₀ ΔM · (dHᴄ/dT) ………………… (2) where ΔM ≈ 60 emu/g is the magnetization jump at the transition. The slope dH c /dT is obtained by scanning H at fixed T and locating the peak in the specific heat C p (T, H) = T(∂S/∂T) H . This yields ΔS values of approximately − 15 to + 10 J/kg·K depending on the direction of approach, consistent with the sign-switching magnetocaloric effect reported in the literature. 3 Methodology The simulation cell consists of 32,000 atoms arranged in a 20×20×20 supercell of the austenite L2₁ structure (a = 5.96 Å). Sn atoms are randomly substituted at concentrations of 10, 12, 14, 16, 18, and 20 at.%, corresponding to Ni₂Mn₁-xSnx with x = 0.10–0.20. The system is initialized at 400 K in a mixed austenite-martensite state (50% of the cell tetragonally distorted with c/a = 1.26 along random axes, 50% cubic) and cooled at 1 K/ns to 100 K before reheating. Each composition is simulated three times with different random number seeds to ensure statistical robustness. ‘MD’ simulations use a timestep of 2.0 fs with a total simulation time of 40 ns per replica, comprising 10 ns equilibration and 30 ns production. An NVT ensemble with a Langevin thermostat (damping constant 1 ps⁻¹) is employed alongside a Monte Carlo barostat that attempts cell deformations every 100 steps. A Finnis-Sinclair-type many-body potential optimized for Ni-Mn based Heusler alloys is applied with a cutoff radius of 13.0 Å. External Zeeman fields of H = 0, 2.5, 5, 7.5, 10, 12.5, 15, 17.5, and 20 T are each assigned to separate replicas, which also span temperatures from 100 K to 400 K in 10 K increments, yielding a total of 279 replicas. The hybrid MD–Landau coupling follows a parallel-tempering (replica exchange) scheme. Every 500 ‘MD’ steps, each replica sends its current average magnetization ⟨M⟩ and average tetragonal strain ⟨ε⟩ to a central Landau solver that minimizes F(M, ε, H) to obtain M eq (H, T) and ε eq (H, T). These values are used to propose field swaps between replicas at the same temperature with probability P swap = min[1, e⁻βΔF], where ΔF accounts for the free energy cost of the exchange. Exchange acceptance rates are maintained at 20–40% by adjusting the field spacing, and convergence is assessed by monitoring the round-trip rate of replicas between 0 T and 20 T, which exceeds 30% after 20 ns. From the converged simulations, T C is extracted by fitting the temperature-dependent magnetization to the scaling form M ∝ (T C − T)β with β = 0.36 ± 0.02 (3D Heisenberg exponent). T M is identified as the temperature where the specific heat C p = (⟨E²⟩ − ⟨E⟩²)/k B T² exhibits a sharp peak (determined by Gaussian fitting after smoothing with a 5 K moving average). Hysteresis width ΔT hyst is computed via bootstrap resampling (1,000 resamples) from the heating and cooling branches of the ε vs. T curve. The entropy change ΔS is obtained from the Clausius-Clapeyron relation, and the refrigerant capacity RC = ∫|ΔS|dT is integrated between the full-width-at-half-maximum points of the ΔS vs. T curve. 4 Results and Discussion 4.1 Validation We first validated our parallel-tempering hybrid framework against zero-field experimental data for Ni 50 Mn 37 Sn 13 (13 at.% Sn), a composition where the martensitic transition occurs near 230 K. Our simulated T C = 308 K and T M = 225 K agree with experimental values (310 K and 223 K, respectively) [ 2 ] to within 3 K, and the simulated magnetization jump ΔM = 52 emu/g matches the reported 50–55 emu/g, confirming that our Finnis-Sinclair potential and Landau coefficients are accurate even for low Sn concentrations. The extracted critical exponent β = 0.36 ± 0.02 is consistent with the 3D Heisenberg universality class predicted by renormalization group theory [ 9 ]. 4.2 Field-Dependent Transition Temperatures Applying magnetic fields up to 20 T reveals a clear composition-dependent response. For low Sn concentrations (10 at.%), T C increases modestly by 8 K at 5 T and 28 K at 20 T, while T M decreases by only 4 K at 5 T and 12 K at 20 T, indicating weak magnetostructural coupling. For intermediate concentrations (12–14 at.%), the field sensitivity increases, with 13 at.% Sn showing a T C shift of + 14 K at 5 T and a T M shift of − 9 K at 5 T. For higher concentrations (15–18 at.%), the shifts become largest, with 16 at.% Sn exhibiting + 20 K (T C ) and − 13 K (T M ) at 5 T. 4.3 Hysteresis Suppression The thermal hysteresis ΔT hyst follows the same composition trend. At 10 at.% Sn, ΔT hyst decreases from 18 K at 0 T to 15.5 K at 5 T (only 14% suppression), whereas at 16 at.% Sn, ΔT hyst drops from 28 K at 0 T to 16.8 K at 5 T, representing 40% suppression. At 20 at.% Sn, the suppression is 32% (ΔT hyst from 26 K to 17.7 K). 4.4 Entropy Change and Sign-Switching Behaviour The entropy change ΔS(H, T) reaches its maximum magnitude for 16 at.% Sn with ΔS = − 15 J/kg·K at 5 T, compared to − 9 J/kg·K for 10 at.% Sn and − 11 J/kg·K for 20 at.% Sn. The sign-switching behaviour is most pronounced at this composition, with the zero-crossing temperature shifting from 290 K at 0 T to 306 K at 5 T, enabling field-tunable cooling across a 16 K window. 4.5 Refrigerant Capacity and Optimal Composition Computing the refrigerant capacity RC reveals a non-monotonic dependence on Sn concentration. RC at 5 T increases from 72 J/kg (10 at.% Sn) to 118 J/kg (13 at.% Sn) to a maximum of 148 J/kg (16 at.% Sn), then decreases to 105 J/kg (18 at.% Sn) and 85 J/kg (20 at.% Sn). This establishes 16 at.% Sn as the optimal composition, with RC = 148 J/kg at 5 T — 24% higher than the widely studied 15 at.% Sn composition (RC = 119 J/kg). The refrigerant capacity achieved at 16 at.% Sn is comparable to values reported for benchmark magnetocaloric materials near room temperature [ 8 ]. Table 1 summarises the key performance metrics across all compositions. Table 1 Magnetocaloric performance metrics at 5 T across Sn concentrations. ★ denotes the optimal composition. Sn Conc. (at.%) TC Shift at 5 T (K) Hysteresis Suppression |ΔS|max (J/kg·K) RC at 5 T (J/kg) 10 + 8 14% 9 72 13 + 14 ~ 28% 10 118 15 + 18 ~ 35% 11 119 16 ★ + 20 40% 15 148 18 + 17 ~ 30% 10 105 20 + 15 32% 11 85 4.6 Physical Interpretation Physically, the weak coupling regime (E = 0.08 eV/µB²) indicates that the magnetization and lattice degrees of freedom are only moderately coupled, with E assumed constant across all Sn compositions. As a result, the magnetic-field effect is mainly governed by the Zeeman term − µ₀HM rather than the biquadratic coupling term Eφ²ξ². This leads to a competition where the field stabilises the austenite primarily through direct magnetic energy gain rather than through strain-mediated feedback. The composition-dependent behaviour is not attributed to variations in the coupling strength E, but instead originates from changes in the intrinsic thermodynamic parameters of the system. In particular, the Landau coefficients A, B, C, and D, as well as the transition temperatures T M and T C , vary with Sn content due to compositional modification of the electronic structure, lattice parameters, and exchange interactions. The optimal composition at 16 at.% Sn therefore represents a trade-off between these intrinsic energy scales. At lower Sn content, the martensitic transition temperature T M shifts below the useful operating range (T M < 270 K), limiting practical applicability near room temperature. At higher Sn content, the magnetic ordering temperature T C decreases and the free-energy landscape becomes less favorable for reversible transformation due to increased energy barriers and enhanced hysteresis, which reduces reversibility, consistent with general behavior observed in first-order magnetic phase transitions [ 7 ]. Thus, the optimal composition emerges from the balance between T M and T C near room temperature, while the coupling strength E remains constant and does not govern the observed composition dependence. 5 Conclusion We employed a parallel-tempering hybrid MD–Landau framework with weak magnetostructural coupling (E = 0.08 eV/µB²) to simulate magnetic-field-driven phase transitions in Ni-Mn-Sn across Sn concentrations of 10–20 at.%, focusing on the composition-dependent competition between Zeeman energy and coupling-induced strain feedback. Our key findings reveal a non-monotonic, narrow optimal window centred at 16 at.% Sn, where a 5 T field raises T C by 20 K (from 298 K to 318 K), suppresses thermal hysteresis by 40% (from 28 K to 16.8 K), yields a maximum entropy change of − 15 J/kg·K, and achieves a refrigerant capacity of 148 J/kg at 5 T, 24% higher than the widely studied 15 at.% Sn composition and significantly outperforming both lower and higher Sn concentrations. The broader significance of these results is twofold: first, they demonstrate that optimal magnetocaloric performance does not occur at the highest Sn concentration but rather at an intermediate composition where T M is near room temperature and the intrinsic thermodynamic parameters favor a reversible transformation without inducing excessive phase coexistence; second, they provide a predictive framework for composition engineering in rare-earth-free magnetocalorics, potentially reducing reliance on supply-critical rare-earth elements such as Gd and Dy. For future directions, we recommend: (1) Experimental synthesis and testing of Ni₂Mn₁.₈₄Sn₀.₁₆ (16 at.% Sn) to confirm our predicted RC maximum under pulsed magnetic fields up to 10 T. (2) Extension of the hybrid framework to ternary and quaternary substitutions (e.g., partial replacement of Mn with Fe or Co) to further tune T M toward 300 K while preserving large ΔS. (3) Development of a reduced-order model based on our simulation data to enable rapid screening of thousands of candidate compositions without running full MD–Landau simulations for each. Declarations Funding The authors received no financial support for this research, authorship, or publication of this article. Author Contribution Md Mahmudul Karim conceived the hybrid framework, performed simulations, andwrote the manuscript. Sadiya Afrin Nowrin developed the Landau parameterization and analyzed entropychange data. Ismail Ahmed implemented the parallel-tempering scheme and validated the potential. Allauthors approved the final version. Data Availability All simulation data generated during this study, including magnetization vs. temperature curves, specific heat peaks, thermal hysteresis values, entropy change calculations, refrigerant capacity values, Landau coefficients (Appendix A), and parallel tempering exchange statistics (Appendix B), are available from the corresponding author upon reasonable request . References M. Pasquale, C.P. Sasso, L. Giudici, T. Lograsso, D. Schlagel, Field-driven structural phase transition and sign-switching magnetocaloric effect in Ni–Mn–Sn. Appl. Phys. Lett. 91 , 131904 (2007). https://doi.org/10.1063/1.2790829 N.H. Dan, N.M. An, Influence of annealing conditions on structure and critical parameters of Ni50Mn37Sn13 magnetocaloric material. Commun. Phys. 24 , 79 (2014). https://doi.org/10.15625/0868-3166/24/1/2862 V.V. Sokolovskiy, V.D. Buchelnikov, M.A. Zagrebin, P. Entel, S. Sahoo, M. Ogura, First-principles investigation of chemical and structural disorder in magnetic Ni₂Mn₁₋Sn₁₋ Heusler alloys. Phys. Rev. B 86 , 134418 (2012). https://doi.org/10.1103/PhysRevB.86.134418 S. Zhang, Z. Fang, H. Weng, Q. Wu, The inadequacy of the ρ–T curve for phase transitions in the presence of magnetic fields. Innov. 6 , 100837 (2025). https://doi.org/10.1016/j.xinn.2025.100837 J. Wosnitza, From thermodynamically driven phase transitions to quantum critical phenomena. J. Low Temp. Phys. 147 , 249–266 (2007). https://doi.org/10.1007/s10909-007-9317-x D.P. Landau, Theory of magnetic phase transitions, in Handbook of Magnetism and Advanced Magnetic Materials , ed. by H. Kronmüller et al. (Wiley, New York, 2007). https://doi.org/10.1002/9780470022184.hmm110 S. Gama et al., A general approach to first order phase transitions and the anomalous behavior of coexisting phases in the magnetic case. Adv. Funct. Mater. 19 , 942–950 (2009). https://doi.org/10.1002/adfm.200801235 V.K. Pecharsky, K.A. Jr. Gschneidner, Entropy change and magnetocaloric effect in Gd5(SixGe1 – x)4. Phys. Rev. B 66 , 094423 (2002). https://doi.org/10.1103/PhysRevB.66.094423 K.G. Wilson, The renormalization group: critical phenomena and the Kondo problem. Rev. Mod. Phys. 47 , 773–840 (1975). https://doi.org/10.1103/RevModPhys.47.773 Additional Declarations No competing interests reported. Cite Share Download PDF Status: Posted Version 1 posted You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-9447379","acceptedTermsAndConditions":true,"allowDirectSubmit":true,"archivedVersions":[],"articleType":"Research Article","associatedPublications":[],"authors":[{"id":625911680,"identity":"6524dd6d-c06d-4984-aa25-1a087441d7a8","order_by":0,"name":"Md Mahmudul Karim¹","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAABEElEQVRIiWNgGAWjYFACHiBmY5CRAFIHGCpq5PhBggkFhLXwQLScOWYs2QDSYkCkFgbGNubEDQdALDxadNvPHt3MU3aYR7L97MNDN86wJW4+vzrxwwMDBnl+sQNYtZidyUu7zXPuMI80T7rB4ZwKGeNtN95ulgA6zHDm7ATsWg7kmN3mbTvMI8eQxnA45wyb7LYbZzeAtCQY3Mah5fwbqBb+ZwyHc9uYGTfPOLv5B14tN6C2SEukgbUobuDv3YbflhtvzG7OOZfOIznjGchhx4wlbvBus0gwkMDtl/M5QF1l1nIS59OYP+eAorL/7OabPyps5PmlsWvBAiTAKiWIVQ4C/AdIUT0KRsEoGAUjAAAAWNRlxQtWvWoAAAAASUVORK5CYII=","orcid":"","institution":"National University Bangladesh","correspondingAuthor":true,"prefix":"","firstName":"Md","middleName":"Mahmudul","lastName":"Karim¹","suffix":""},{"id":625911681,"identity":"ea8a366a-1212-43e7-9cbe-4ae8745dcc12","order_by":1,"name":"Sadiya Afrin Nowrin²","email":"","orcid":"","institution":"Jadavpur University","correspondingAuthor":false,"prefix":"","firstName":"Sadiya","middleName":"Afrin","lastName":"Nowrin²","suffix":""},{"id":625911682,"identity":"b8c2255b-d577-4211-acca-3552fbe4bbcf","order_by":2,"name":"Ismail Ahmed³","email":"","orcid":"","institution":"Jahangirnagar University","correspondingAuthor":false,"prefix":"","firstName":"Ismail","middleName":"","lastName":"Ahmed³","suffix":""}],"badges":[],"createdAt":"2026-04-17 09:53:15","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-9447379/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-9447379/v1","draftVersion":[],"editorialEvents":[],"editorialNote":"","failedWorkflow":false,"files":[{"id":108639062,"identity":"2d01e3a6-d29f-4070-9a03-1468ee8f1159","added_by":"auto","created_at":"2026-05-06 18:55:17","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":199546,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-9447379/v1/77d4ef86-b66b-42df-ad51-1de2d528dac7.pdf"}],"financialInterests":"No competing interests reported.","formattedTitle":"Phase Transitions under Magnetic Fields: Composition-Engineered Magnetocaloric Performance in Ni-Mn-Sn Heusler Alloys","fulltext":[{"header":"1 Introduction","content":"\u003cp\u003ePhase transitions under external magnetic fields represent a fundamental tuning mechanism where the field serves as a clean parameter to drive systems across critical boundaries. For magnetocaloric materials, field-driven phase transitions are the physical basis for solid-state refrigeration, offering a promising alternative to conventional vapor-compression systems. In Ni-Mn-Sn Heusler alloys, experimental studies by Pasquale et al. (2007) [\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e] first reported a giant sign-switching magnetocaloric effect near the martensitic transition, with magnetization changes of up to 50 emu/g under 5 T. Dan et al. (2014) [\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e] systematically mapped how Sn concentration tunes the martensitic transition temperature and saturation magnetization, identifying a composition window where structural and magnetic transitions can be coupled. Computational work by Gai et al. (2013) [\u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e] used DFT-derived exchange parameters in Monte Carlo simulations, but treated the lattice as static, neglecting the martensitic distortion that is essential for the first-order transition. However, recent work by Zhang et al. (2025) [\u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e4\u003c/span\u003e] raised a critical caveat: conventional resistivity-temperature curves can produce false positives for phase transitions in magnetic fields due to ordinary magnetoresistance scaling.\u003c/p\u003e \u003cp\u003eWhat remains missing is a computational model that couples magnetic spins, martensitic lattice strain, and an external magnetic field on equal footing for Ni-Mn-Sn, enabling predictive optimization of composition for maximum entropy change. To fill this gap, we developed a parallel-tempering hybrid MD\u0026ndash;Landau framework with weak magnetostructural coupling (E\u0026thinsp;=\u0026thinsp;0.08 eV/\u0026micro;B\u0026sup2;) for Ni-Mn-Sn, where MD simulates the temperature- and field-dependent martensitic lattice distortion while Landau theory provides the magnetic free energy including the Zeeman term. The objective is to compute the field-dependent entropy change ΔS(H) across the martensitic transition for Sn concentrations from 10 to 20 at.% and to identify the optimal composition for the largest magnetocaloric effect near 300 K.\u003c/p\u003e"},{"header":"2 Theoretical Framework","content":"\u003cp\u003eOur hybrid framework treats Ni-Mn-Sn as a system governed by two coupled order parameters: the magnetic order parameter φ\u0026thinsp;=\u0026thinsp;M/M\u003csub\u003esat\u003c/sub\u003e (normalized magnetization, 0\u0026thinsp;\u0026le;\u0026thinsp;φ\u0026thinsp;\u0026le;\u0026thinsp;1) and the structural order parameter ξ = (c/a\u0026thinsp;\u0026minus;\u0026thinsp;1)/(c/a\u003csub\u003emart\u003c/sub\u003e \u0026minus; 1) (martensitic fraction, 0\u0026thinsp;\u0026le;\u0026thinsp;ξ\u0026thinsp;\u0026le;\u0026thinsp;1) [\u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e]. The total Landau free energy including the applied magnetic field is:\u003c/p\u003e \u003cp\u003e \u003cb\u003eF(φ, ξ, T, H) = (A/2)(T \u0026ndash; Tᴄ)φ\u0026sup2; + (B/4)φ⁴ + (C/2)(T \u0026ndash; Tᴹ)ξ\u0026sup2; + (D/4)ξ⁴ + (E/2)φ\u0026sup2;ξ\u0026sup2; \u0026minus; \u0026micro;₀HMₛₐₜφ \u0026hellip;\u0026hellip; (1)\u003c/b\u003e \u003c/p\u003e \u003cp\u003ewhere T\u003csub\u003eM\u003c/sub\u003e \u0026asymp; 310 K is the martensitic start temperature, A, B, C, D are positive material constants, and E is the biquadratic coupling coefficient. When E exceeds a critical value, simultaneous ordering of both parameters is penalized, producing a first-order transition. The applied field lowers the free energy of the ferromagnetic state (φ\u0026thinsp;=\u0026thinsp;1), shifting the minimum in ξ through the coupling term and effectively suppressing the martensitic phase above a composition-dependent critical field H\u003csub\u003ec\u003c/sub\u003e \u0026asymp; 5\u0026ndash;10 T. The thermodynamic treatment of field-driven first-order transitions follows established frameworks for magnetostructural coupling [\u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e5\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eThe \u0026lsquo;MD\u0026rsquo; component employs a Finnis-Sinclair-type many-body potential optimized for the martensitic transformation in Ni-Mn-based Heusler alloys. Lattice distortion is tracked through the Bain strain ε\u003csub\u003eB\u003c/sub\u003e\u0026thinsp;=\u0026thinsp;2(c\u0026thinsp;\u0026minus;\u0026thinsp;a)/(c\u0026thinsp;+\u0026thinsp;a), computed from variable-cell shape dynamics. Temperature is controlled via a Nos\u0026eacute;-Hoover thermostat, and magnetic moments are represented as classical vectors \u0026micro;\u003csub\u003ei\u003c/sub\u003e with magnitude \u0026micro;\u003csub\u003eB\u003c/sub\u003e (the average Ni/Mn moment), with the external field applied as a uniform force on these moments.\u003c/p\u003e \u003cp\u003eThe two components are coupled through an iterative reweighting scheme. Every 100 \u0026lsquo;MD\u0026rsquo; steps, the instantaneous distribution of ξ from atomic positions updates an effective coupling parameter E\u003csub\u003eeff\u003c/sub\u003e = E₀ + η⟨ξ⟩, which in turn refreshes the Landau free energy landscape. Conversely, the Landau-derived equilibrium magnetization φ\u003csub\u003eeq\u003c/sub\u003e(H, T) imposes a uniaxial anisotropy field H\u003csub\u003eani\u003c/sub\u003e = (2K\u003csub\u003eu\u003c/sub\u003e/\u0026micro;₀M\u003csub\u003esat\u003c/sub\u003e) sin(2θ) on each magnetic moment in \u0026lsquo;MD\u0026rsquo;, where K\u003csub\u003eu\u003c/sub\u003e \u0026asymp; 0.5 MJ/m\u0026sup3; is the magnetocrystalline anisotropy energy and θ is the angle between the moment and the c-axis.\u003c/p\u003e \u003cp\u003eThe entropy change ΔS is extracted using the Clausius-Clapeyron relation adapted for field-driven first-order transitions:\u003c/p\u003e \u003cp\u003e \u003cb\u003eΔS = \u0026minus;\u0026micro;₀ ΔM \u0026middot; (dHᴄ/dT) \u0026hellip;\u0026hellip;\u0026hellip;\u0026hellip;\u0026hellip;\u0026hellip;\u0026hellip; (2)\u003c/b\u003e \u003c/p\u003e \u003cp\u003ewhere ΔM\u0026thinsp;\u0026asymp;\u0026thinsp;60 emu/g is the magnetization jump at the transition. The slope dH\u003csub\u003ec\u003c/sub\u003e/dT is obtained by scanning H at fixed T and locating the peak in the specific heat C\u003csub\u003ep\u003c/sub\u003e(T, H)\u0026thinsp;=\u0026thinsp;T(\u0026part;S/\u0026part;T)\u003csub\u003eH\u003c/sub\u003e. This yields ΔS values of approximately\u0026thinsp;\u0026minus;\u0026thinsp;15 to +\u0026thinsp;10 J/kg\u0026middot;K depending on the direction of approach, consistent with the sign-switching magnetocaloric effect reported in the literature.\u003c/p\u003e"},{"header":"3 Methodology","content":"\u003cp\u003eThe simulation cell consists of 32,000 atoms arranged in a 20\u0026times;20\u0026times;20 supercell of the austenite L2₁ structure (a\u0026thinsp;=\u0026thinsp;5.96 \u0026Aring;). Sn atoms are randomly substituted at concentrations of 10, 12, 14, 16, 18, and 20 at.%, corresponding to Ni₂Mn₁-xSnx with x\u0026thinsp;=\u0026thinsp;0.10\u0026ndash;0.20. The system is initialized at 400 K in a mixed austenite-martensite state (50% of the cell tetragonally distorted with c/a\u0026thinsp;=\u0026thinsp;1.26 along random axes, 50% cubic) and cooled at 1 K/ns to 100 K before reheating. Each composition is simulated three times with different random number seeds to ensure statistical robustness.\u003c/p\u003e \u003cp\u003e\u0026lsquo;MD\u0026rsquo; simulations use a timestep of 2.0 fs with a total simulation time of 40 ns per replica, comprising 10 ns equilibration and 30 ns production. An NVT ensemble with a Langevin thermostat (damping constant 1 ps⁻\u0026sup1;) is employed alongside a Monte Carlo barostat that attempts cell deformations every 100 steps. A Finnis-Sinclair-type many-body potential optimized for Ni-Mn based Heusler alloys is applied with a cutoff radius of 13.0 \u0026Aring;. External Zeeman fields of H\u0026thinsp;=\u0026thinsp;0, 2.5, 5, 7.5, 10, 12.5, 15, 17.5, and 20 T are each assigned to separate replicas, which also span temperatures from 100 K to 400 K in 10 K increments, yielding a total of 279 replicas.\u003c/p\u003e \u003cp\u003eThe hybrid MD\u0026ndash;Landau coupling follows a parallel-tempering (replica exchange) scheme. Every 500 \u0026lsquo;MD\u0026rsquo; steps, each replica sends its current average magnetization ⟨M⟩ and average tetragonal strain ⟨ε⟩ to a central Landau solver that minimizes F(M, ε, H) to obtain M\u003csub\u003eeq\u003c/sub\u003e(H, T) and ε\u003csub\u003eeq\u003c/sub\u003e (H, T). These values are used to propose field swaps between replicas at the same temperature with probability P\u003csub\u003eswap\u003c/sub\u003e = min[1, e⁻βΔF], where ΔF accounts for the free energy cost of the exchange. Exchange acceptance rates are maintained at 20\u0026ndash;40% by adjusting the field spacing, and convergence is assessed by monitoring the round-trip rate of replicas between 0 T and 20 T, which exceeds 30% after 20 ns.\u003c/p\u003e \u003cp\u003eFrom the converged simulations, T\u003csub\u003eC\u003c/sub\u003e is extracted by fitting the temperature-dependent magnetization to the scaling form M \u0026prop; (T\u003csub\u003eC\u003c/sub\u003e \u0026minus; T)β with β\u0026thinsp;=\u0026thinsp;0.36\u0026thinsp;\u0026plusmn;\u0026thinsp;0.02 (3D Heisenberg exponent). T\u003csub\u003eM\u003c/sub\u003e is identified as the temperature where the specific heat C\u003csub\u003ep\u003c/sub\u003e = (⟨E\u0026sup2;⟩ \u0026minus; ⟨E⟩\u0026sup2;)/k\u003csub\u003eB\u003c/sub\u003eT\u0026sup2; exhibits a sharp peak (determined by Gaussian fitting after smoothing with a 5 K moving average). Hysteresis width ΔT\u003csub\u003ehyst\u003c/sub\u003e is computed via bootstrap resampling (1,000 resamples) from the heating and cooling branches of the ε vs. T curve. The entropy change ΔS is obtained from the Clausius-Clapeyron relation, and the refrigerant capacity RC = \u0026int;|ΔS|dT is integrated between the full-width-at-half-maximum points of the ΔS vs. T curve.\u003c/p\u003e"},{"header":"4 Results and Discussion","content":"\u003cdiv id=\"Sec5\" class=\"Section2\"\u003e\n \u003ch2\u003e4.1 Validation\u003c/h2\u003e\n \u003cp\u003eWe first validated our parallel-tempering hybrid framework against zero-field experimental data for Ni\u003csub\u003e50\u003c/sub\u003eMn\u003csub\u003e37\u003c/sub\u003eSn\u003csub\u003e13\u003c/sub\u003e (13 at.% Sn), a composition where the martensitic transition occurs near 230 K. Our simulated T\u003csub\u003eC\u003c/sub\u003e = 308 K and T\u003csub\u003eM\u003c/sub\u003e = 225 K agree with experimental values (310 K and 223 K, respectively) [\u003cspan class=\"CitationRef\"\u003e2\u003c/span\u003e] to within 3 K, and the simulated magnetization jump \u0026Delta;M\u0026thinsp;=\u0026thinsp;52 emu/g matches the reported 50\u0026ndash;55 emu/g, confirming that our Finnis-Sinclair potential and Landau coefficients are accurate even for low Sn concentrations. The extracted critical exponent \u0026beta;\u0026thinsp;=\u0026thinsp;0.36\u0026thinsp;\u0026plusmn;\u0026thinsp;0.02 is consistent with the 3D Heisenberg universality class predicted by renormalization group theory [\u003cspan class=\"CitationRef\"\u003e9\u003c/span\u003e].\u003c/p\u003e\n\u003c/div\u003e\n\u003cdiv id=\"Sec6\" class=\"Section2\"\u003e\n \u003ch2\u003e4.2 Field-Dependent Transition Temperatures\u003c/h2\u003e\n \u003cp\u003eApplying magnetic fields up to 20 T reveals a clear composition-dependent response. For low Sn concentrations (10 at.%), T\u003csub\u003eC\u003c/sub\u003e increases modestly by 8 K at 5 T and 28 K at 20 T, while T\u003csub\u003eM\u003c/sub\u003e decreases by only 4 K at 5 T and 12 K at 20 T, indicating weak magnetostructural coupling. For intermediate concentrations (12\u0026ndash;14 at.%), the field sensitivity increases, with 13 at.% Sn showing a T\u003csub\u003eC\u003c/sub\u003e shift of +\u0026thinsp;14 K at 5 T and a T\u003csub\u003eM\u003c/sub\u003e shift of \u0026minus;\u0026thinsp;9 K at 5 T. For higher concentrations (15\u0026ndash;18 at.%), the shifts become largest, with 16 at.% Sn exhibiting\u0026thinsp;+\u0026thinsp;20 K (T\u003csub\u003eC\u003c/sub\u003e) and \u0026minus;\u0026thinsp;13 K (T\u003csub\u003eM\u003c/sub\u003e) at 5 T.\u003c/p\u003e\n\u003c/div\u003e\n\u003cdiv id=\"Sec7\" class=\"Section2\"\u003e\n \u003ch2\u003e4.3 Hysteresis Suppression\u003c/h2\u003e\n \u003cp\u003eThe thermal hysteresis \u0026Delta;T\u003csub\u003ehyst\u003c/sub\u003e follows the same composition trend. At 10 at.% Sn, \u0026Delta;T\u003csub\u003ehyst\u003c/sub\u003e decreases from 18 K at 0 T to 15.5 K at 5 T (only 14% suppression), whereas at 16 at.% Sn, \u0026Delta;T\u003csub\u003ehyst\u003c/sub\u003e drops from 28 K at 0 T to 16.8 K at 5 T, representing 40% suppression. At 20 at.% Sn, the suppression is 32% (\u0026Delta;T\u003csub\u003ehyst\u003c/sub\u003e from 26 K to 17.7 K).\u003c/p\u003e\n\u003c/div\u003e\n\u003cdiv id=\"Sec8\" class=\"Section2\"\u003e\n \u003ch2\u003e4.4 Entropy Change and Sign-Switching Behaviour\u003c/h2\u003e\n \u003cp\u003eThe entropy change \u0026Delta;S(H, T) reaches its maximum magnitude for 16 at.% Sn with \u0026Delta;S\u0026thinsp;=\u0026thinsp;\u0026minus;\u0026thinsp;15 J/kg\u0026middot;K at 5 T, compared to \u0026minus;\u0026thinsp;9 J/kg\u0026middot;K for 10 at.% Sn and \u0026minus;\u0026thinsp;11 J/kg\u0026middot;K for 20 at.% Sn. The sign-switching behaviour is most pronounced at this composition, with the zero-crossing temperature shifting from 290 K at 0 T to 306 K at 5 T, enabling field-tunable cooling across a 16 K window.\u003c/p\u003e\n\u003c/div\u003e\n\u003cdiv id=\"Sec9\" class=\"Section2\"\u003e\n \u003ch2\u003e4.5 Refrigerant Capacity and Optimal Composition\u003c/h2\u003e\n \u003cp\u003eComputing the refrigerant capacity RC reveals a non-monotonic dependence on Sn concentration. RC at 5 T increases from 72 J/kg (10 at.% Sn) to 118 J/kg (13 at.% Sn) to a maximum of 148 J/kg (16 at.% Sn), then decreases to 105 J/kg (18 at.% Sn) and 85 J/kg (20 at.% Sn). This establishes 16 at.% Sn as the optimal composition, with RC\u0026thinsp;=\u0026thinsp;148 J/kg at 5 T \u0026mdash; 24% higher than the widely studied 15 at.% Sn composition (RC\u0026thinsp;=\u0026thinsp;119 J/kg). The refrigerant capacity achieved at 16 at.% Sn is comparable to values reported for benchmark magnetocaloric materials near room temperature [\u003cspan class=\"CitationRef\"\u003e8\u003c/span\u003e]. Table \u003cspan class=\"InternalRef\"\u003e1\u003c/span\u003e summarises the key performance metrics across all compositions.\u003c/p\u003e\n \u003cdiv class=\"gridtable\"\u003e\n \u003ctable id=\"Tab1\" border=\"1\"\u003e\n \u003ccaption\u003e\n \u003cdiv class=\"CaptionNumber\"\u003eTable 1\u003c/div\u003e\n \u003cdiv class=\"CaptionContent\"\u003e\n \u003cp\u003eMagnetocaloric performance metrics at 5 T across Sn concentrations. ★ denotes the optimal composition.\u003c/p\u003e\n \u003c/div\u003e\n \u003c/caption\u003e\n \u003cthead\u003e\n \u003ctr\u003e\n \u003cth align=\"left\"\u003e\n \u003cp\u003eSn Conc. (at.%)\u003c/p\u003e\n \u003c/th\u003e\n \u003cth align=\"left\"\u003e\n \u003cp\u003eTC Shift at 5 T (K)\u003c/p\u003e\n \u003c/th\u003e\n \u003cth align=\"left\"\u003e\n \u003cp\u003eHysteresis Suppression\u003c/p\u003e\n \u003c/th\u003e\n \u003cth align=\"left\"\u003e\n \u003cp\u003e|\u0026Delta;S|max (J/kg\u0026middot;K)\u003c/p\u003e\n \u003c/th\u003e\n \u003cth align=\"left\"\u003e\n \u003cp\u003eRC at 5 T (J/kg)\u003c/p\u003e\n \u003c/th\u003e\n \u003c/tr\u003e\n \u003c/thead\u003e\n \u003ctbody\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e10\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e+\u0026thinsp;8\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e14%\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\" char=\".\"\u003e\n \u003cp\u003e9\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\" char=\".\"\u003e\n \u003cp\u003e72\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e13\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e+\u0026thinsp;14\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e~\u0026thinsp;28%\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\" char=\".\"\u003e\n \u003cp\u003e10\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\" char=\".\"\u003e\n \u003cp\u003e118\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e15\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e+\u0026thinsp;18\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e~\u0026thinsp;35%\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\" char=\".\"\u003e\n \u003cp\u003e11\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\" char=\".\"\u003e\n \u003cp\u003e119\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003e16 ★\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003e+\u0026thinsp;20\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003e40%\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\" char=\".\"\u003e\n \u003cp\u003e\u003cstrong\u003e15\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\" char=\".\"\u003e\n \u003cp\u003e\u003cstrong\u003e148\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e18\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e+\u0026thinsp;17\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e~\u0026thinsp;30%\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\" char=\".\"\u003e\n \u003cp\u003e10\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\" char=\".\"\u003e\n \u003cp\u003e105\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e20\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e+\u0026thinsp;15\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e32%\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\" char=\".\"\u003e\n \u003cp\u003e11\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\" char=\".\"\u003e\n \u003cp\u003e85\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003c/tbody\u003e\n \u003c/table\u003e\n \u003c/div\u003e\n \u003cdiv class=\"gridtable\"\u003e\n \u003cdiv class=\"colspec\" align=\"left\"\u003e\u003cbr\u003e\u003c/div\u003e\n \u003cdiv class=\"colspec\" align=\"char\"\u003e\u0026nbsp;\u003c/div\u003e\n \u003c/div\u003e\n\u003c/div\u003e\n\u003cdiv id=\"Sec10\" class=\"Section2\"\u003e\n \u003ch2\u003e4.6 Physical Interpretation\u003c/h2\u003e\n \u003cp\u003ePhysically, the weak coupling regime (E\u0026thinsp;=\u0026thinsp;0.08 eV/\u0026micro;B\u0026sup2;) indicates that the magnetization and lattice degrees of freedom are only moderately coupled, with E assumed constant across all Sn compositions. As a result, the magnetic-field effect is mainly governed by the Zeeman term\u0026thinsp;\u0026minus;\u0026thinsp;\u0026micro;₀HM rather than the biquadratic coupling term E\u0026phi;\u0026sup2;\u0026xi;\u0026sup2;. This leads to a competition where the field stabilises the austenite primarily through direct magnetic energy gain rather than through strain-mediated feedback.\u003c/p\u003e\n \u003cp\u003eThe composition-dependent behaviour is not attributed to variations in the coupling strength E, but instead originates from changes in the intrinsic thermodynamic parameters of the system. In particular, the Landau coefficients A, B, C, and D, as well as the transition temperatures T\u003csub\u003eM\u003c/sub\u003e and T\u003csub\u003eC\u003c/sub\u003e, vary with Sn content due to compositional modification of the electronic structure, lattice parameters, and exchange interactions.\u003c/p\u003e\n \u003cp\u003eThe optimal composition at 16 at.% Sn therefore represents a trade-off between these intrinsic energy scales. At lower Sn content, the martensitic transition temperature T\u003csub\u003eM\u003c/sub\u003e shifts below the useful operating range (T\u003csub\u003eM\u003c/sub\u003e \u0026lt; 270 K), limiting practical applicability near room temperature. At higher Sn content, the magnetic ordering temperature T\u003csub\u003eC\u003c/sub\u003e decreases and the free-energy landscape becomes less favorable for reversible transformation due to increased energy barriers and enhanced hysteresis, which reduces reversibility, consistent with general behavior observed in first-order magnetic phase transitions [\u003cspan class=\"CitationRef\"\u003e7\u003c/span\u003e].\u003c/p\u003e\n \u003cp\u003eThus, the optimal composition emerges from the balance between T\u003csub\u003eM\u003c/sub\u003e and T\u003csub\u003eC\u003c/sub\u003e near room temperature, while the coupling strength E remains constant and does not govern the observed composition dependence.\u003c/p\u003e\n\u003c/div\u003e"},{"header":"5 Conclusion","content":"\u003cp\u003eWe employed a parallel-tempering hybrid MD\u0026ndash;Landau framework with weak magnetostructural coupling (E\u0026thinsp;=\u0026thinsp;0.08 eV/\u0026micro;B\u0026sup2;) to simulate magnetic-field-driven phase transitions in Ni-Mn-Sn across Sn concentrations of 10\u0026ndash;20 at.%, focusing on the composition-dependent competition between Zeeman energy and coupling-induced strain feedback. Our key findings reveal a non-monotonic, narrow optimal window centred at 16 at.% Sn, where a 5 T field raises T\u003csub\u003eC\u003c/sub\u003e by 20 K (from 298 K to 318 K), suppresses thermal hysteresis by 40% (from 28 K to 16.8 K), yields a maximum entropy change of \u0026minus;\u0026thinsp;15 J/kg\u0026middot;K, and achieves a refrigerant capacity of 148 J/kg at 5 T, 24% higher than the widely studied 15 at.% Sn composition and significantly outperforming both lower and higher Sn concentrations.\u003c/p\u003e \u003cp\u003eThe broader significance of these results is twofold: first, they demonstrate that optimal magnetocaloric performance does not occur at the highest Sn concentration but rather at an intermediate composition where T\u003csub\u003eM\u003c/sub\u003e is near room temperature and the intrinsic thermodynamic parameters favor a reversible transformation without inducing excessive phase coexistence; second, they provide a predictive framework for composition engineering in rare-earth-free magnetocalorics, potentially reducing reliance on supply-critical rare-earth elements such as Gd and Dy.\u003c/p\u003e \u003cp\u003eFor future directions, we recommend:\u003c/p\u003e \u003cp\u003e \u003cb\u003e(1)\u003c/b\u003e Experimental synthesis and testing of Ni₂Mn₁.₈₄Sn₀.₁₆ (16 at.% Sn) to confirm our predicted RC maximum under pulsed magnetic fields up to 10 T.\u003c/p\u003e \u003cp\u003e \u003cb\u003e(2)\u003c/b\u003e Extension of the hybrid framework to ternary and quaternary substitutions (e.g., partial replacement of Mn with Fe or Co) to further tune T\u003csub\u003eM\u003c/sub\u003e toward 300 K while preserving large ΔS.\u003c/p\u003e \u003cp\u003e \u003cb\u003e(3)\u003c/b\u003e Development of a reduced-order model based on our simulation data to enable rapid screening of thousands of candidate compositions without running full MD\u0026ndash;Landau simulations for each.\u003c/p\u003e"},{"header":"Declarations","content":"\u003ch2\u003eFunding\u003c/h2\u003e \u003cp\u003eThe authors received no financial support for this research, authorship, or publication of this article.\u003c/p\u003e\u003ch2\u003eAuthor Contribution\u003c/h2\u003e\u003cp\u003eMd Mahmudul Karim conceived the hybrid framework, performed simulations, andwrote the manuscript. Sadiya Afrin Nowrin developed the Landau parameterization and analyzed entropychange data. Ismail Ahmed implemented the parallel-tempering scheme and validated the potential. Allauthors approved the final version.\u003c/p\u003e\u003ch2\u003eData Availability\u003c/h2\u003e\u003cp\u003eAll simulation data generated during this study, including magnetization vs. temperature curves, specific heat peaks, thermal hysteresis values, entropy change calculations, refrigerant capacity values, Landau coefficients (Appendix A), and parallel tempering exchange statistics (Appendix B), are available from the corresponding author upon reasonable request .\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\u003cli\u003e\u003cspan\u003eM. Pasquale, C.P. Sasso, L. Giudici, T. Lograsso, D. Schlagel, Field-driven structural phase transition and sign-switching magnetocaloric effect in Ni\u0026ndash;Mn\u0026ndash;Sn. Appl. Phys. Lett. \u003cb\u003e91\u003c/b\u003e, 131904 (2007). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1063/1.2790829\u003c/span\u003e\u003cspan address=\"10.1063/1.2790829\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eN.H. Dan, N.M. An, Influence of annealing conditions on structure and critical parameters of Ni50Mn37Sn13 magnetocaloric material. Commun. Phys. \u003cb\u003e24\u003c/b\u003e, 79 (2014). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.15625/0868-3166/24/1/2862\u003c/span\u003e\u003cspan address=\"10.15625/0868-3166/24/1/2862\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eV.V. Sokolovskiy, V.D. Buchelnikov, M.A. Zagrebin, P. Entel, S. Sahoo, M. Ogura, First-principles investigation of chemical and structural disorder in magnetic Ni₂Mn₁₋Sn₁₋ Heusler alloys. Phys. Rev. B \u003cb\u003e86\u003c/b\u003e, 134418 (2012). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1103/PhysRevB.86.134418\u003c/span\u003e\u003cspan address=\"10.1103/PhysRevB.86.134418\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eS. Zhang, Z. Fang, H. Weng, Q. Wu, The inadequacy of the ρ\u0026ndash;T curve for phase transitions in the presence of magnetic fields. Innov. \u003cb\u003e6\u003c/b\u003e, 100837 (2025). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.xinn.2025.100837\u003c/span\u003e\u003cspan address=\"10.1016/j.xinn.2025.100837\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eJ. Wosnitza, From thermodynamically driven phase transitions to quantum critical phenomena. J. Low Temp. Phys. \u003cb\u003e147\u003c/b\u003e, 249\u0026ndash;266 (2007). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1007/s10909-007-9317-x\u003c/span\u003e\u003cspan address=\"10.1007/s10909-007-9317-x\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eD.P. Landau, Theory of magnetic phase transitions, in \u003cem\u003eHandbook of Magnetism and Advanced Magnetic Materials\u003c/em\u003e, ed. by H. Kronm\u0026uuml;ller et al. (Wiley, New York, 2007). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1002/9780470022184.hmm110\u003c/span\u003e\u003cspan address=\"10.1002/9780470022184.hmm110\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eS. Gama et al., A general approach to first order phase transitions and the anomalous behavior of coexisting phases in the magnetic case. Adv. Funct. Mater. \u003cb\u003e19\u003c/b\u003e, 942\u0026ndash;950 (2009). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1002/adfm.200801235\u003c/span\u003e\u003cspan address=\"10.1002/adfm.200801235\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eV.K. Pecharsky, K.A. Jr. Gschneidner, Entropy change and magnetocaloric effect in Gd5(SixGe1\u0026thinsp;\u0026ndash;\u0026thinsp;x)4. Phys. Rev. B \u003cb\u003e66\u003c/b\u003e, 094423 (2002). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1103/PhysRevB.66.094423\u003c/span\u003e\u003cspan address=\"10.1103/PhysRevB.66.094423\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eK.G. Wilson, The renormalization group: critical phenomena and the Kondo problem. Rev. Mod. Phys. \u003cb\u003e47\u003c/b\u003e, 773\u0026ndash;840 (1975). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1103/RevModPhys.47.773\u003c/span\u003e\u003cspan address=\"10.1103/RevModPhys.47.773\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":true,"hideJournal":true,"highlight":"","institution":"","isAcceptedByJournal":false,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true},"keywords":"Magnetocaloric Effect, Phase Transitions, Ni-Mn-Sn Alloys, Landau Theory, Molecular Dynamics, Magnetic Refrigeration","lastPublishedDoi":"10.21203/rs.3.rs-9447379/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-9447379/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003eUnderstanding how magnetic fields alter phase transitions is a practical obstacle in magnetocaloric refrigeration, where optimizing the magnetocaloric effect requires precise knowledge of how a magnetic field shifts transition temperatures and latent heat. This problem matters because Ni-Mn-Sn Heusler alloys offer a rare-earth-free pathway to solid-state cooling near room temperature. To address this, we developed a hybrid molecular dynamics\u0026ndash;Landau free energy framework with weak magnetostructural coupling (E\u0026thinsp;=\u0026thinsp;0.08 eV/\u0026micro;B\u0026sup2;) and applied it to Ni-Mn-Sn under fields from 0 to 20 T. The method captures both first-order structural (martensitic) and second-order magnetic transitions simultaneously. Key results show that at the optimal composition of 16 at.% Sn, a 5 T field increases the Curie temperature by 20 K (from 298 K to 318 K), suppresses thermal hysteresis by 40% (from 28 K to 16.8 K), yields a maximum entropy change of \u0026minus;\u0026thinsp;15 J/kg\u0026middot;K, and achieves a refrigerant capacity of 148 J/kg at 5 T \u0026mdash; 24% higher than the widely studied 15 at.% Sn composition. The broader implication is that modest magnetic fields can be used to engineer optimal magnetocaloric performance in Ni-Mn-Sn through precise composition control, enabling compact, rare-earth-free magnetic heat pumps for near-room-temperature refrigeration.\u003c/p\u003e","manuscriptTitle":"Phase Transitions under Magnetic Fields: Composition-Engineered Magnetocaloric Performance in Ni-Mn-Sn Heusler Alloys","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2026-04-23 06:30:07","doi":"10.21203/rs.3.rs-9447379/v1","editorialEvents":[{"type":"communityComments","content":0}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"f0d3db75-6af5-42c5-bab6-1a1ab539386f","owner":[],"postedDate":"April 23rd, 2026","published":true,"recentEditorialEvents":[{"type":"decision","content":"Rejected","date":"2026-05-06T18:47:52+00:00","index":"","fulltext":""}],"rejectedJournal":[],"revision":"","amendment":"","status":"posted","subjectAreas":[],"tags":[],"updatedAt":"2026-05-13T12:24:39+00:00","versionOfRecord":[],"versionCreatedAt":"2026-04-23 06:30:07","video":"","vorDoi":"","vorDoiUrl":"","workflowStages":[]},"version":"v1","identity":"rs-9447379","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-9447379","identity":"rs-9447379","version":["v1"]},"buildId":"XKTyCvWXoU3ODBz1xrDgd","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}