Pub Date : 2026-01-17DOI: 10.1016/j.commatsci.2026.114504
Apurba Sarker, Sourav Saha
The future of space exploration and human settlement beyond Earth hinges on a deeper understanding of in-space manufacturing processes. The unique physical conditions and scarcity of experimental data demand robust computational models to investigate the atomic-scale physics of solidification. This work presents a molecular dynamics (MD) model to examine the solidification behavior of the – binary alloy, focusing on the influence of varying compositions and gravity levels (Earth, Lunar, Martian, and microgravity) on atomistic solidification mechanisms and the resulting mechanical properties — specifically, hardness — of as-solidified nanostructures. Hardness is evaluated via nanoindentation simulations. The study confirms that gravitational forces significantly affect the solidification pathways of – alloys. Notably, by tuning alloy composition, the influence of gravity can be modulated—and in some cases, even reversed. Moreover, hardness exhibits a coupled dependence on both composition and gravity, offering a promising avenue for bottom-up design of components tailored for extraterrestrial environments. The article delves into the nanoscale physical mechanisms underlying these phenomena and outlines future directions for extending this modeling framework to broader applications.
{"title":"Gravity and composition modulated solidification and mechanical properties of Al-Cu nanostructures","authors":"Apurba Sarker, Sourav Saha","doi":"10.1016/j.commatsci.2026.114504","DOIUrl":"10.1016/j.commatsci.2026.114504","url":null,"abstract":"<div><div>The future of space exploration and human settlement beyond Earth hinges on a deeper understanding of in-space manufacturing processes. The unique physical conditions and scarcity of experimental data demand robust computational models to investigate the atomic-scale physics of solidification. This work presents a molecular dynamics (MD) model to examine the solidification behavior of the <span><math><mrow><mi>A</mi><mi>l</mi></mrow></math></span>–<span><math><mrow><mi>C</mi><mi>u</mi></mrow></math></span> binary alloy, focusing on the influence of varying compositions and gravity levels (Earth, Lunar, Martian, and microgravity) on atomistic solidification mechanisms and the resulting mechanical properties — specifically, hardness — of <em>as-solidified</em> nanostructures. Hardness is evaluated via nanoindentation simulations. The study confirms that gravitational forces significantly affect the solidification pathways of <span><math><mrow><mi>A</mi><mi>l</mi></mrow></math></span>–<span><math><mrow><mi>C</mi><mi>u</mi></mrow></math></span> alloys. Notably, by tuning alloy composition, the influence of gravity can be modulated—and in some cases, even reversed. Moreover, hardness exhibits a coupled dependence on both composition and gravity, offering a promising avenue for <em>bottom-up</em> design of components tailored for extraterrestrial environments. The article delves into the nanoscale physical mechanisms underlying these phenomena and outlines future directions for extending this modeling framework to broader applications.</div></div>","PeriodicalId":10650,"journal":{"name":"Computational Materials Science","volume":"265 ","pages":"Article 114504"},"PeriodicalIF":3.3,"publicationDate":"2026-01-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145976387","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"材料科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-17DOI: 10.1016/j.commatsci.2026.114516
Yuanyuan Tian , Jinfeng Du , Zhiyang Li , Xin He , Benxi Zhang
Impingement of deposited droplets with coming ones is regarded as a fundamental process in terms of many practical applications, especially for micro-nano systems. However, the understanding of such systems from dynamic evolution to effect of characteristic parameters remains limited and unsatisfactory. To address this gap, the current study employs molecular dynamics (MD) simulations to unravel the process of impacting deposited droplets with coming on the molecular level. By exploring the effect of various parameter groups, We, Δ, and θ0, systematically, we map all observed outcomes into a series of phase diagrams that encompass four distinct regimes, including deposition, regular bounce, bounce with closed holes, and breakup bounce. The free evolution of targeted systems, critical We between different regimes, and contact time have been investigated and discussed through directly extracting data from numerical simulations and theoretical calculations. Ultimately, we develop a novel system with wrapped nanoparticles to quickly sweep deposited droplets in order to recover hydrophobicity of solid surfaces. The underlying mechanisms are attributable to enhanced energy conversion between merged droplets and wrapped nanoparticles. The current work helps researchers obtain comprehensively insights into the nanoscale collision of droplets with unequal size. Moreover, this approach opens a window to recover the superhydrophobicity of solid surfaces that are contaminated by deposited water droplets.
{"title":"Impingement of nanoscale droplets upon deposited ones at different given conditions: A molecular dynamics study","authors":"Yuanyuan Tian , Jinfeng Du , Zhiyang Li , Xin He , Benxi Zhang","doi":"10.1016/j.commatsci.2026.114516","DOIUrl":"10.1016/j.commatsci.2026.114516","url":null,"abstract":"<div><div>Impingement of deposited droplets with coming ones is regarded as a fundamental process in terms of many practical applications, especially for micro-nano systems. However, the understanding of such systems from dynamic evolution to effect of characteristic parameters remains limited and unsatisfactory. To address this gap, the current study employs molecular dynamics (MD) simulations to unravel the process of impacting deposited droplets with coming on the molecular level. By exploring the effect of various parameter groups, <em>We</em>, <em>Δ</em>, and <em>θ</em><sub>0</sub>, systematically, we map all observed outcomes into a series of phase diagrams that encompass four distinct regimes, including deposition, regular bounce, bounce with closed holes, and breakup bounce. The free evolution of targeted systems, critical <em>We</em> between different regimes, and contact time have been investigated and discussed through directly extracting data from numerical simulations and theoretical calculations. Ultimately, we develop a novel system with wrapped nanoparticles to quickly sweep deposited droplets in order to recover hydrophobicity of solid surfaces. The underlying mechanisms are attributable to enhanced energy conversion between merged droplets and wrapped nanoparticles. The current work helps researchers obtain comprehensively insights into the nanoscale collision of droplets with unequal size. Moreover, this approach opens a window to recover the superhydrophobicity of solid surfaces that are contaminated by deposited water droplets.</div></div>","PeriodicalId":10650,"journal":{"name":"Computational Materials Science","volume":"265 ","pages":"Article 114516"},"PeriodicalIF":3.3,"publicationDate":"2026-01-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145976388","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"材料科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-16DOI: 10.1016/j.commatsci.2026.114514
Yang Sun , Jia Liu , Guangzhao Deng , Shuan Ma , Dengbao Xiao
This study proposes an innovative micromechanics-based deep neural network method to efficiently investigate the effects of fiber shape and interphase on the thermoelastic properties of unidirectional composites. Firstly, this work establishes a micromechanical finite element approach by simulating the internal microstructure of the composite and verifies its rationality by comparing it with experimental results. Subsequently, using DOE sampling method based on global arrangement, data groups for training are obtained through the finite element simulation, and the machine learning model is further constructed utilizing deep neural network algorithm. The effectiveness of the machine learning model is validated by comparing the true values from the finite element simulation with the predicted values from the machine learning. Finally, a comprehensive investigation is conducted to elucidate the effects of fiber concentration and morphology, interphase concentration and characteristics on the thermoelastic behavior of composites. The results show that the established machine learning model provides a fast and accurate prediction for the thermoelastic properties of composites considering microstructural features.
{"title":"A machine learning approach to modeling the effects of fiber shape and interphase on the thermoelastic properties of composites","authors":"Yang Sun , Jia Liu , Guangzhao Deng , Shuan Ma , Dengbao Xiao","doi":"10.1016/j.commatsci.2026.114514","DOIUrl":"10.1016/j.commatsci.2026.114514","url":null,"abstract":"<div><div>This study proposes an innovative micromechanics-based deep neural network method to efficiently investigate the effects of fiber shape and interphase on the thermoelastic properties of unidirectional composites. Firstly, this work establishes a micromechanical finite element approach by simulating the internal microstructure of the composite and verifies its rationality by comparing it with experimental results. Subsequently, using DOE sampling method based on global arrangement, data groups for training are obtained through the finite element simulation, and the machine learning model is further constructed utilizing deep neural network algorithm. The effectiveness of the machine learning model is validated by comparing the true values from the finite element simulation with the predicted values from the machine learning. Finally, a comprehensive investigation is conducted to elucidate the effects of fiber concentration and morphology, interphase concentration and characteristics on the thermoelastic behavior of composites. The results show that the established machine learning model provides a fast and accurate prediction for the thermoelastic properties of composites considering microstructural features.</div></div>","PeriodicalId":10650,"journal":{"name":"Computational Materials Science","volume":"265 ","pages":"Article 114514"},"PeriodicalIF":3.3,"publicationDate":"2026-01-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145976386","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"材料科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-16DOI: 10.1016/j.commatsci.2026.114510
Imane Bezzaoui , Soukaina Elhalimi , Abdelmajid El Badraoui , Abdelmajid Elmansouri , Najim Tahiri , Omar El Bounagui
In this study, we investigate the hydrogen storage potential of KMgH3 using a comprehensive first-principles approach that combines density functional theory (DFT) and ab initio molecular dynamics (AIMD) simulations. We examine the effect of 25% potassium substitution with aluminum and titanium (Al,Ti), the introduction of potassium vacancies, and the application of compressive strain on the doped KMgH3 structures, focusing on their structural stability, thermodynamic, hydrogen storage, and electronic properties. Our results show that both doping and vacancy engineering significantly reduce the formation enthalpy and desorption temperature, improving dehydrogenation performance. Notably, potassium vacancies increase the gravimetric capacity from 4.55 wt% to 5.33 wt%, while compressive strain reduces the formation enthalpy to −40.62 kJ/mol.H₂ and the desorption temperature to 313.46 K. The system exhibits structural and thermal stability under these modifications, as confirmed by elastic constant analysis and AIMD simulations at ambient temperature. To explain this reduction, we examined the electronic structure, which reveals that the reduction in potassium atoms introduces states at the Fermi level, indicating enhanced electrical conductivity and weakened metal‑hydrogen interactions that facilitate hydrogen release. In summary, these findings demonstrate the ability of this combined approach to optimize KMgH3 for advanced hydrogen storage applications.
{"title":"Tuning the hydrogen storage properties of KMgH3 through metal doping, vacancy engineering, and compressive strain: A first principles study","authors":"Imane Bezzaoui , Soukaina Elhalimi , Abdelmajid El Badraoui , Abdelmajid Elmansouri , Najim Tahiri , Omar El Bounagui","doi":"10.1016/j.commatsci.2026.114510","DOIUrl":"10.1016/j.commatsci.2026.114510","url":null,"abstract":"<div><div>In this study, we investigate the hydrogen storage potential of KMgH<sub>3</sub> using a comprehensive first-principles approach that combines density functional theory (DFT) and ab initio molecular dynamics (AIMD) simulations. We examine the effect of 25% potassium substitution with aluminum and titanium (Al,Ti), the introduction of potassium vacancies, and the application of compressive strain on the doped KMgH<sub>3</sub> structures, focusing on their structural stability, thermodynamic, hydrogen storage, and electronic properties. Our results show that both doping and vacancy engineering significantly reduce the formation enthalpy and desorption temperature, improving dehydrogenation performance. Notably, potassium vacancies increase the gravimetric capacity from 4.55 wt% to 5.33 wt%, while compressive strain reduces the formation enthalpy to −40.62 kJ/mol.H₂ and the desorption temperature to 313.46 K. The system exhibits structural and thermal stability under these modifications, as confirmed by elastic constant analysis and AIMD simulations at ambient temperature. To explain this reduction, we examined the electronic structure, which reveals that the reduction in potassium atoms introduces states at the Fermi level, indicating enhanced electrical conductivity and weakened metal‑hydrogen interactions that facilitate hydrogen release. In summary, these findings demonstrate the ability of this combined approach to optimize KMgH<sub>3</sub> for advanced hydrogen storage applications.</div></div>","PeriodicalId":10650,"journal":{"name":"Computational Materials Science","volume":"265 ","pages":"Article 114510"},"PeriodicalIF":3.3,"publicationDate":"2026-01-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145969405","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"材料科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-15DOI: 10.1016/j.commatsci.2026.114492
Michael J.D. Rushton , Michael W.D. Cooper , Ghanshyam Pilania , Blas P. Uberuaga
While density functional theory (DFT) has become the de facto approach for accurate simulation of materials at the atomic scale, there are many aspects of materials that are simply out of reach of DFT methods. In particular, finite temperature properties such as diffusivities, the structure and properties of grain boundaries and interfaces, and the study of defect properties in complex alloys are computationally challenging for DFT methods. Recently, a new class of spinels in which three cations order over two sublattices was discovered. In order to predict the properties of these types of structures, classical potentials are a must. In this work, we derive a new classical potential for Mg-bearing spinels in which the B cations are Al, Ga, and/or In. The potential does well in describing the DFT energetics of various spinel structures as a function of chemistry and inversion. In particular, it reproduces the thermodynamically favorable MgAlGaO4 structure while correctly predicting that neither MgAlInO4 nor MgGaInO4 are stable. Further, it reproduces physical trends in elastic properties as compared against experiment.
{"title":"A new interatomic potential for mixed Mg-Al-Ga-In spinels","authors":"Michael J.D. Rushton , Michael W.D. Cooper , Ghanshyam Pilania , Blas P. Uberuaga","doi":"10.1016/j.commatsci.2026.114492","DOIUrl":"10.1016/j.commatsci.2026.114492","url":null,"abstract":"<div><div>While density functional theory (DFT) has become the de facto approach for accurate simulation of materials at the atomic scale, there are many aspects of materials that are simply out of reach of DFT methods. In particular, finite temperature properties such as diffusivities, the structure and properties of grain boundaries and interfaces, and the study of defect properties in complex alloys are computationally challenging for DFT methods. Recently, a new class of spinels in which three cations order over two sublattices was discovered. In order to predict the properties of these types of structures, classical potentials are a must. In this work, we derive a new classical potential for Mg-bearing spinels in which the B cations are Al, Ga, and/or In. The potential does well in describing the DFT energetics of various spinel structures as a function of chemistry and inversion. In particular, it reproduces the thermodynamically favorable MgAlGaO<sub>4</sub> structure while correctly predicting that neither MgAlInO<sub>4</sub> nor MgGaInO<sub>4</sub> are stable. Further, it reproduces physical trends in elastic properties as compared against experiment.</div></div>","PeriodicalId":10650,"journal":{"name":"Computational Materials Science","volume":"264 ","pages":"Article 114492"},"PeriodicalIF":3.3,"publicationDate":"2026-01-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145974141","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"材料科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-15DOI: 10.1016/j.commatsci.2026.114503
Chaitanya Bhave, Somayajulu L.N. Dhulipala, Mathew Swisher, Jacob A. Hirschhorn, Ryan Terrence Sweet, Stephen R. Novascone
TRistructural ISOtropic (TRISO) particle fuel relies on a silicon carbide (SiC) layer as the primary structural material and barrier to metallic fission products (FPs) release. Accurate prediction of palladium (Pd) transport and penetration into the SiC is therefore critical for qualifying TRISO fuels for advanced reactors. The empirical correlation for Pd penetration in BISON is derived from historical particle-fuel data, but it cannot explain the large scatter in the experimental data that arises from varying experimental conditions. To aid fuel qualification, we previously developed a mechanistic reduced order model (ROM) using BISON that resolves these dependencies (Bhave et al., 2025). In this work we built on that mechanistic ROM, validated it, and quantified its uncertainty using Bayesian uncertainty quantification (UQ). We calibrated against a suite of in-pile and out-of-pile experiments spanning particle compositions, geometries, and operating conditions, and benchmarked the mechanistic ROM against the empirical correlation. We used Bayesian UQ to identify influential parameters and calibrate them to data, which yielded predictive intervals. Results show that while the empirical correlation can be tuned to fit a single experiment type, it transfers poorly; the mechanistic ROM sustains accuracy with credible uncertainty across disparate conditions. This process demonstrates a practical path — via Bayesian UQ applied to mechanistic ROMs — to leverage single-effect experiments for inferring in-reactor behavior and supporting TRISO fuel qualification.
三结构各向同性(TRISO)粒子燃料依靠碳化硅(SiC)层作为主要结构材料和金属裂变产物(FPs)释放的屏障。因此,准确预测钯(Pd)在碳化硅中的传输和渗透对于先进反应堆的TRISO燃料的资格至关重要。BISON中Pd渗透的经验相关性来源于历史颗粒-燃料数据,但它不能解释实验数据中由于不同实验条件而产生的大分散。为了帮助燃料鉴定,我们之前使用BISON开发了一种机制降order模型(ROM)来解决这些依赖关系(Bhave et al., 2025)。在这项工作中,我们建立了机械ROM,验证了它,并使用贝叶斯不确定性量化(UQ)量化了它的不确定性。我们根据一套桩内和桩外实验进行了校准,涵盖了颗粒组成、几何形状和操作条件,并根据经验相关性对机械ROM进行了基准测试。我们使用贝叶斯UQ来识别有影响的参数,并将其校准为数据,从而产生预测区间。结果表明,虽然经验相关性可以调整到适合单一实验类型,但它的转移性很差;机械式只读存储器在不同的条件下保持具有可靠的不确定性的准确性。该过程展示了一种实用的途径——通过将贝叶斯UQ应用于机械rom——利用单效应实验来推断反应堆内行为并支持TRISO燃料鉴定。
{"title":"Bayesian discovery of optimal reduced order models from mechanistic and experimental data: A case study of Pd penetration in TRISO fuels using BISON","authors":"Chaitanya Bhave, Somayajulu L.N. Dhulipala, Mathew Swisher, Jacob A. Hirschhorn, Ryan Terrence Sweet, Stephen R. Novascone","doi":"10.1016/j.commatsci.2026.114503","DOIUrl":"10.1016/j.commatsci.2026.114503","url":null,"abstract":"<div><div>TRistructural ISOtropic (TRISO) particle fuel relies on a silicon carbide (SiC) layer as the primary structural material and barrier to metallic fission products (FPs) release. Accurate prediction of palladium (Pd) transport and penetration into the SiC is therefore critical for qualifying TRISO fuels for advanced reactors. The empirical correlation for Pd penetration in BISON is derived from historical particle-fuel data, but it cannot explain the large scatter in the experimental data that arises from varying experimental conditions. To aid fuel qualification, we previously developed a mechanistic reduced order model (ROM) using BISON that resolves these dependencies (Bhave et al., 2025). In this work we built on that mechanistic ROM, validated it, and quantified its uncertainty using Bayesian uncertainty quantification (UQ). We calibrated against a suite of in-pile and out-of-pile experiments spanning particle compositions, geometries, and operating conditions, and benchmarked the mechanistic ROM against the empirical correlation. We used Bayesian UQ to identify influential parameters and calibrate them to data, which yielded predictive intervals. Results show that while the empirical correlation can be tuned to fit a single experiment type, it transfers poorly; the mechanistic ROM sustains accuracy with credible uncertainty across disparate conditions. This process demonstrates a practical path — via Bayesian UQ applied to mechanistic ROMs — to leverage single-effect experiments for inferring in-reactor behavior and supporting TRISO fuel qualification.</div></div>","PeriodicalId":10650,"journal":{"name":"Computational Materials Science","volume":"264 ","pages":"Article 114503"},"PeriodicalIF":3.3,"publicationDate":"2026-01-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145974140","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"材料科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-14DOI: 10.1016/j.commatsci.2026.114508
Ignacio Palos-Reynoso , Humberto Híjar
We present a neural network-based algorithm for the identification and classification of topological defects in two-dimensional nematic liquid crystals confined within square geometries. The nematic configurations are generated through Nematic-Multiparticle Collision Dynamics, a mesoscopic simulation method that captures both hydrodynamic and orientational fluctuations. Our supervised learning framework is trained on synthetic images labeled with topological defects of positive and negative charge. In the inference stage, we employ local winding number estimations to propose candidate defect locations, which are then evaluated by the neural network to determine their authenticity and type. The algorithm achieves robust classification performance, with a macro-averaged score of 0.92, indicating balanced precision and recall across all defect classes, with most misclassifications arising from director field fluctuations near the corners of the confinement domain. Beyond static identification, our method enables the temporal tracking of defect dynamics, including annihilation events. This work demonstrates the potential of deep learning tools to extract and quantify topological information in fluctuating soft matter systems.
{"title":"Deep learning detection of topological defects in confined two-dimensional nematics","authors":"Ignacio Palos-Reynoso , Humberto Híjar","doi":"10.1016/j.commatsci.2026.114508","DOIUrl":"10.1016/j.commatsci.2026.114508","url":null,"abstract":"<div><div>We present a neural network-based algorithm for the identification and classification of topological defects in two-dimensional nematic liquid crystals confined within square geometries. The nematic configurations are generated through Nematic-Multiparticle Collision Dynamics, a mesoscopic simulation method that captures both hydrodynamic and orientational fluctuations. Our supervised learning framework is trained on synthetic images labeled with topological defects of positive and negative charge. In the inference stage, we employ local winding number estimations to propose candidate defect locations, which are then evaluated by the neural network to determine their authenticity and type. The algorithm achieves robust classification performance, with a macro-averaged <span><math><msub><mrow><mi>F</mi></mrow><mrow><mn>1</mn></mrow></msub></math></span> score of 0.92, indicating balanced precision and recall across all defect classes, with most misclassifications arising from director field fluctuations near the corners of the confinement domain. Beyond static identification, our method enables the temporal tracking of defect dynamics, including annihilation events. This work demonstrates the potential of deep learning tools to extract and quantify topological information in fluctuating soft matter systems.</div></div>","PeriodicalId":10650,"journal":{"name":"Computational Materials Science","volume":"264 ","pages":"Article 114508"},"PeriodicalIF":3.3,"publicationDate":"2026-01-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145974222","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"材料科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-13DOI: 10.1016/j.commatsci.2026.114505
Mary A. Mazannikova , Vladimir I. Anisimov , Dmitry Y. Novoselov
Layered electrides, characterized by anionic electrons confined in interstitial sites, present a unique platform for engineering exotic electronic and magnetic phenomena. This study employs a combination of density functional theory, maximally localized Wannier functions, and dynamical mean-field theory to systematically investigate the emergence and control of magnetism in a family of twelve isostructural electrides (M Ca, Sr, Ba; X N, P, As, Sb). We demonstrate that the magnetic state is governed by the local geometry of the interstitial cavities, specifically by the ratio of intra- to inter-layer metal–metal distances (). A magnetic ground state emerges when this ratio falls below unity, a condition that can be selectively induced by hydrostatic pressure. Electronic structure analysis reveals that this transition is driven by a Stoner-like instability, associated with the flattening of an electride-derived band at the Fermi level. Our DMFT calculations confirm the presence of significant electron correlations and spin fluctuations near the magnetic instability, indicative of a correlated metallic state. The strong coupling between magnetic ordering and the crystal lattice, evidenced by concurrent structural and magnetic phase transitions, underscores a robust magneto-structural coupling. We establish simple empirical criteria based on atomic radii and electronegativities to predict magnetic behavior within this family of compounds. These findings provide a comprehensive microscopic understanding of magnetism in layered electrides and establish design principles for creating and tuning magnetic materials via pressure or chemical substitution from non-magnetic elements.
层状电子,其特征是阴离子电子被限制在间隙位置,为工程奇异的电子和磁现象提供了一个独特的平台。本研究采用密度泛函理论、最大定域万涅尔函数和动力学平均场理论相结合的方法,系统地研究了12种M2X等结构电子(M = Ca, Sr, Ba; X = N, P, As, Sb)中磁性的产生和控制。我们证明了磁性状态是由间隙腔的局部几何形状控制的,特别是由层内与层间金属-金属距离的比率(lintra/linter)控制的。当这个比率低于1时,磁性基态就会出现,这种情况可以由静水压力选择性地诱导。电子结构分析表明,这种转变是由一种类似斯通纳的不稳定性驱动的,这种不稳定性与费米能级上电极衍生带的平坦化有关。我们的DMFT计算证实了磁不稳定性附近存在显著的电子相关性和自旋波动,表明存在相关的金属态。磁有序与晶格之间的强耦合,通过同时发生的结构和磁相变证明,强调了强磁-结构耦合。我们建立了基于原子半径和电负性的简单经验准则来预测这类化合物的磁性行为。这些发现为层状电子中的磁性提供了全面的微观理解,并建立了通过压力或非磁性元素的化学替代来创建和调整磁性材料的设计原则。
{"title":"Electrides: From fundamental concepts to tunable magnetism in layered systems","authors":"Mary A. Mazannikova , Vladimir I. Anisimov , Dmitry Y. Novoselov","doi":"10.1016/j.commatsci.2026.114505","DOIUrl":"10.1016/j.commatsci.2026.114505","url":null,"abstract":"<div><div>Layered electrides, characterized by anionic electrons confined in interstitial sites, present a unique platform for engineering exotic electronic and magnetic phenomena. This study employs a combination of density functional theory, maximally localized Wannier functions, and dynamical mean-field theory to systematically investigate the emergence and control of magnetism in a family of twelve isostructural <span><math><mrow><msub><mrow><mi>M</mi></mrow><mrow><mn>2</mn></mrow></msub><mi>X</mi></mrow></math></span> electrides (M <span><math><mo>=</mo></math></span> Ca, Sr, Ba; X <span><math><mo>=</mo></math></span> N, P, As, Sb). We demonstrate that the magnetic state is governed by the local geometry of the interstitial cavities, specifically by the ratio of intra- to inter-layer metal–metal distances (<span><math><mrow><msub><mrow><mi>l</mi></mrow><mrow><mtext>intra</mtext></mrow></msub><mo>/</mo><msub><mrow><mi>l</mi></mrow><mrow><mtext>inter</mtext></mrow></msub></mrow></math></span>). A magnetic ground state emerges when this ratio falls below unity, a condition that can be selectively induced by hydrostatic pressure. Electronic structure analysis reveals that this transition is driven by a Stoner-like instability, associated with the flattening of an electride-derived band at the Fermi level. Our DMFT calculations confirm the presence of significant electron correlations and spin fluctuations near the magnetic instability, indicative of a correlated metallic state. The strong coupling between magnetic ordering and the crystal lattice, evidenced by concurrent structural and magnetic phase transitions, underscores a robust magneto-structural coupling. We establish simple empirical criteria based on atomic radii and electronegativities to predict magnetic behavior within this family of compounds. These findings provide a comprehensive microscopic understanding of magnetism in layered electrides and establish design principles for creating and tuning magnetic materials via pressure or chemical substitution from non-magnetic elements.</div></div>","PeriodicalId":10650,"journal":{"name":"Computational Materials Science","volume":"264 ","pages":"Article 114505"},"PeriodicalIF":3.3,"publicationDate":"2026-01-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145974145","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"材料科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-13DOI: 10.1016/j.commatsci.2026.114486
Veeramohan Rao M
Aluminosilicate glasses incorporating alkaline earth metals, Zn, and Mn are of significant interest within the realms of material science and geoscience. The structural implications of the MO/Al2O3 ratio, specifically at a value of 3, in aluminosilicate glasses remain incompletely elucidated. In the present study, molecular dynamic simulations are employed to investigate the influence of cation field strength on the structural and elastic characteristics of aluminosilicate glasses. The cation field strength is observed to induce alterations in structural properties attributes such as bond length, coordination number, bond angle distribution, and the presence of oxygen species including FO, NBO, BO, and TBO. The glass transition temperature and elastic constants are determined through these simulations. It is found that an increase in cation field strength correlates with an elevation in elastic constants and a reduction in glass transition temperature. These findings provide atomic-scale insights into the effects of cation field strength on the properties of glasses.
{"title":"Molecular dynamic simulation study on 3MOAl2O3 3SiO2 [M = Ba, Sr, Ca, Mg, Zn and Mn] glasses","authors":"Veeramohan Rao M","doi":"10.1016/j.commatsci.2026.114486","DOIUrl":"10.1016/j.commatsci.2026.114486","url":null,"abstract":"<div><div>Aluminosilicate glasses incorporating alkaline earth metals, Zn, and Mn are of significant interest within the realms of material science and geoscience. The structural implications of the MO/Al<sub>2</sub>O<sub>3</sub> ratio, specifically at a value of 3, in aluminosilicate glasses remain incompletely elucidated. In the present study, molecular dynamic simulations are employed to investigate the influence of cation field strength on the structural and elastic characteristics of aluminosilicate glasses. The cation field strength is observed to induce alterations in structural properties attributes such as bond length, coordination number, bond angle distribution, and the presence of oxygen species including FO, NBO, BO, and TBO. The glass transition temperature and elastic constants are determined through these simulations. It is found that an increase in cation field strength correlates with an elevation in elastic constants and a reduction in glass transition temperature. These findings provide atomic-scale insights into the effects of cation field strength on the properties of glasses.</div></div>","PeriodicalId":10650,"journal":{"name":"Computational Materials Science","volume":"264 ","pages":"Article 114486"},"PeriodicalIF":3.3,"publicationDate":"2026-01-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145974142","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"材料科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-13DOI: 10.1016/j.commatsci.2025.114443
Avik Mahata
We present a second-nearest-neighbor Modified Embedded Atom Method (2NN–MEAM) potential for Scandium (Sc) and Aluminum-Scandium (Al–Sc) alloys that unifies cohesive, thermodynamic, and solidification behavior within a single transferable framework. The Sc component accurately reproduces cohesive energy, lattice constants, defect energetics, and the experimental melting point obtained from two-phase coexistence, demonstrating reliable description of both hcp and liquid phases. The Al–Sc binary interaction parameters were fitted using the L1–AlSc reference and benchmarked against first-principles and calorimetric data. The potential reproduces the strong negative formation enthalpy of AlSc (–0.45 eV atom−1), correct relative stability of competing phases, and realistic elastic properties. Mixing enthalpies of the liquid alloy agree with ideal-associated-solution and CALPHAD models, confirming that the potential captures exothermic Al–Sc association in the melt. Molecular-dynamics simulations of solidification reveal the expected temperature and composition dependence of homogeneous nucleation. Pure Al crystallizes readily, while Al–1 at.% Sc exhibits a longer incubation and slower growth at the same absolute temperature due to reduced undercooling and solute drag. Within the alloy, ordered AlSc-type L1 embryos appear spontaneously, with Sc atoms occupying cube-corner (B) sites surrounded by twelve Al neighbors. Energy–volume trajectories confirm that the potential links thermodynamics to microstructural evolution. Overall, the developed 2NN–MEAM potential provides a quantitatively grounded basis for modeling melting, solidification, and intermetallic ordering in Sc and Al–Sc systems, enabling future multicomponent alloy design and large-scale nucleation studies.
{"title":"Development and validation of interatomic potential for Sc and Al–Sc alloys: Thermodynamics, solidification, and intermetallic ordering","authors":"Avik Mahata","doi":"10.1016/j.commatsci.2025.114443","DOIUrl":"10.1016/j.commatsci.2025.114443","url":null,"abstract":"<div><div>We present a second-nearest-neighbor Modified Embedded Atom Method (2NN–MEAM) potential for Scandium (Sc) and Aluminum-Scandium (Al–Sc) alloys that unifies cohesive, thermodynamic, and solidification behavior within a single transferable framework. The Sc component accurately reproduces cohesive energy, lattice constants, defect energetics, and the experimental melting point obtained from two-phase coexistence, demonstrating reliable description of both hcp and liquid phases. The Al–Sc binary interaction parameters were fitted using the L1<span><math><msub><mrow></mrow><mrow><mn>2</mn></mrow></msub></math></span>–Al<span><math><msub><mrow></mrow><mrow><mn>3</mn></mrow></msub></math></span>Sc reference and benchmarked against first-principles and calorimetric data. The potential reproduces the strong negative formation enthalpy of Al<span><math><msub><mrow></mrow><mrow><mn>3</mn></mrow></msub></math></span>Sc (–0.45 eV atom<sup>−1</sup>), correct relative stability of competing phases, and realistic elastic properties. Mixing enthalpies of the liquid alloy agree with ideal-associated-solution and CALPHAD models, confirming that the potential captures exothermic Al–Sc association in the melt. Molecular-dynamics simulations of solidification reveal the expected temperature and composition dependence of homogeneous nucleation. Pure Al crystallizes readily, while Al–1 at.% Sc exhibits a longer incubation and slower growth at the same absolute temperature due to reduced undercooling and solute drag. Within the alloy, ordered Al<span><math><msub><mrow></mrow><mrow><mn>3</mn></mrow></msub></math></span>Sc-type L1<span><math><msub><mrow></mrow><mrow><mn>2</mn></mrow></msub></math></span> embryos appear spontaneously, with Sc atoms occupying cube-corner (B) sites surrounded by twelve Al neighbors. Energy–volume trajectories confirm that the potential links thermodynamics to microstructural evolution. Overall, the developed 2NN–MEAM potential provides a quantitatively grounded basis for modeling melting, solidification, and intermetallic ordering in Sc and Al–Sc systems, enabling future multicomponent alloy design and large-scale nucleation studies.</div></div>","PeriodicalId":10650,"journal":{"name":"Computational Materials Science","volume":"264 ","pages":"Article 114443"},"PeriodicalIF":3.3,"publicationDate":"2026-01-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145974143","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"材料科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}