Direct conversion of syngas to light olefins(STO)on bifunctional catalysts has garnered significant attention,yet a comprehensive understanding of the reaction pathway and reaction kinetics remains elusive.Herein,we t...Direct conversion of syngas to light olefins(STO)on bifunctional catalysts has garnered significant attention,yet a comprehensive understanding of the reaction pathway and reaction kinetics remains elusive.Herein,we theoretically addressed the kinetics of the direct STO reaction on typical ZnAl_(2)O_(4)/zeolite catalysts by establishing a complete reaction network,consisting of methanol synthesis and conversion,water gas shift(WGS)reaction,olefin hydrogenation,and other relevant steps.The WGS reaction occurs very readily on ZnAl_(2)O_(4) surface whereas which is less active towards alkane formation via olefin hydrogenation,and the latter can be attributed to the characteristics of the H_(2) heterolytic activation and the weak polarity of olefins.The driving effect of zeolite component towards CO conversion was demonstrated by microkinetic simulations,which is sensitive to reaction conditions like space velocity and reaction temperature.Under a fixed ratio of active sites between oxide and zeolite components,the concept of the“impossible trinity”of high CO conversion,high olefin selectivity,and high space velocity can thus be manifested.This work thus provides a comprehensive kinetic picture on the direct STO conversion,offering valuable insights for the design of each component of bifunctional catalysts and the optimization of reaction conditions.展开更多
The development of high-performance solid electrolytes is pivotal for advancing solid-state battery technologies.In this work,we design an oxysulfide-based solid electrolyte Na MgPO_(3)S by combining bond valence theo...The development of high-performance solid electrolytes is pivotal for advancing solid-state battery technologies.In this work,we design an oxysulfide-based solid electrolyte Na MgPO_(3)S by combining bond valence theory and density functional theory calculations.The material features a wide band gap of 4.0 eV and a considerable reduced Na^(+)migration barrier of 0.44 eV,a 1.26-eV decrease compared to pristine Na MgPO_(4)(~1.70 eV).Ab initio molecular dynamics simulations further reveal significantly enhanced ionic conductivity in the oxysulfide-based system compared to the pristine oxide structure.In addition,the calculated decomposition energy indicates that the modified material exhibits good moisture stability.Our findings suggest that sulfur-doping strategy can simultaneously achieve improved ionic conductivity and high moisture stability in oxide solid electrolytes,which could pave the way for designing high-performance solid electrolytes.展开更多
High-performance intelligent protective materials are vital for nuclear energy systems exposed to extreme irradiation.Among them,tungsten-based alloys demonstrate exceptional potential owing to their superior irradiat...High-performance intelligent protective materials are vital for nuclear energy systems exposed to extreme irradiation.Among them,tungsten-based alloys demonstrate exceptional potential owing to their superior irradiation resistance.Recent experimental studies have demonstrated that W-Ta-Cr-V alloys exhibit excellent irradiation resistance under helium(He)ion irradiation.However,the underlying mechanisms,especially the migration behavior of He atoms,remain unclear.In this work,the influences of uniaxial tensile and compressive strain on He migration in W-Ta-Cr-V complex alloys have been systematically studied through first-principles calculations.Our results demonstrate that He atoms preferentially occupy the tetrahedral interstitial sites,with interstitial formation energies significantly reduced compared to pure W.The introduction of Ta,Cr,and V alloying elements markedly increases the He migration barriers,effectively suppressing He diffusion.Compressive strain increases the migration barriers,inhibiting He bubbles nucleation and growth,while tensile strain decreases the barriers,facilitating bubble formation.Compared to pure W,the W-Ta-Cr-V alloys exhibit both lower He interstitial formation energies and higher migration barriers,with further enhancement under compressive strain.Specifically,compressive strain of 6%increases the He migration barrier of the W-Ta-Cr-V alloy by 0.166 e V,which further widens the difference relative to pure W.These findings provide a theoretical explanation for the superior irradiation resistance of tungsten-based alloys observed experimentally and promote the understanding of irradiation damage in these alloys under strain.展开更多
The rapid advancement of machine learning based tight-binding Hamiltonian(MLTB)methods has opened new avenues for efficient and accurate electronic structure simulations,particularly in large-scale systems and long-ti...The rapid advancement of machine learning based tight-binding Hamiltonian(MLTB)methods has opened new avenues for efficient and accurate electronic structure simulations,particularly in large-scale systems and long-time scenarios.This review begins with a concise overview of traditional tight-binding(TB)models,including both(semi-)empirical and first-principles approaches,establishing the foundation for understanding MLTB developments.We then present a systematic classification of existing MLTB methodologies,grouped into two major categories:direct prediction of TB Hamiltonian elements and inference of empirical parameters.A comparative analysis with other ML-based electronic structure models is also provided,highlighting the advancement of MLTB approaches.Finally,we explore the emerging MLTB application ecosystem,highlighting how the integration of MLTB models with a diverse suite of post-processing tools from linear-scaling solvers to quantum transport frameworks and molecular dynamics interfaces is essential for tackling complex scientific problems across different domains.The continued advancement of this integrated paradigm promises to accelerate materials discovery and open new frontiers in the predictive simulation of complex quantum phenomena.展开更多
Vitrimers belong to a class of polymeric materials capable of bond exchange reactions,showing great promise for environmental protection and sustainable development.However,studies on the coupling mechanism between th...Vitrimers belong to a class of polymeric materials capable of bond exchange reactions,showing great promise for environmental protection and sustainable development.However,studies on the coupling mechanism between the bond exchange kinetics and segmental dynamics near the glass transition temperature(T_(g))remain scarce.Herein,we employed molecular dynamics simulations to investigate the dynamic heterogeneity of the segment motion and bond exchange in vitrimers.The simulation results revealed that the bond exchange energy barrier exerts a much stronger influence on the bond exchange kinetics than on the segmental dynamics.At lower temperatures,slower segmental relaxation further constraind the bond exchange rate.Additionally,increasing the bond exchange energy barrier markedly enhanced the dynamic heterogeneity of segment motion.A close correlation was observed between heterogeneity and bond exchange.This study elucidated the coupling mechanism between bond exchange and segmental dynamics at the molecular scale,thereby providing a theoretical basis for designing vitrimer materials with tunable dynamic properties.展开更多
Compared to the well-studied two-dimensional(2D)ferroelectricity,the appearance of 2D antiferroelectricity is much rarer,where local dipoles from the nonequivalent sublattices within 2D monolayers are oppositely orien...Compared to the well-studied two-dimensional(2D)ferroelectricity,the appearance of 2D antiferroelectricity is much rarer,where local dipoles from the nonequivalent sublattices within 2D monolayers are oppositely oriented.Using NbOCl_(2) monolayer with competing ferroelectric(FE)and antiferroelectric(AFE)phases as a 2D material platform,we demonstrate the emergence of intrinsic antiferroelectricity in NbOCl_(2) monolayer under experimentally accessible shear strain,along with new functionality associated with electric field-induced AFE-to-FE phase transition.Specifically,the complex configuration space accommodating FE and AFE phases,polarization switching kinetics,and finite temperature thermodynamic properties of 2D NbOCl_(2) are all accurately predicted by large-scale molecular dynamics simulations based on deep learning interatomic potential model.Moreover,room temperature stable antiferroelectricity with low polarization switching barrier and one-dimensional collinear polarization arrangement is predicted in shear-deformed NbOCl_(2) monolayer.The transition from AFE to FE phase in 2D NbOCl_(2) can be triggered by a low critical electric field,leading to a double polarization–electric(P–E)loop with small hysteresis.A new type of optoelectronic device composed of AFE-NbOCl_(2) is proposed,enabling electric“writing”and nonlinear optical“reading”logical operation with fast operation speed and low power consumption.展开更多
The first-principles calculations based on density functional theory combined with cluster expansion techniques and Monte Carlo （MC） simulations were used to study the phase diagrams of both wurtzite （WZ） and zinc...The first-principles calculations based on density functional theory combined with cluster expansion techniques and Monte Carlo （MC） simulations were used to study the phase diagrams of both wurtzite （WZ） and zinc-blende （ZB） Cdl_xZnxS alloys. All formation energies are positive for WZ and ZB Cdl-xZnxS alloys, which means that the Cdl-xZnxS alloys are unstable and have a tendency to phase separation. For WZ and ZB Cdl_xZnxS alloys, the consolute temperatures are 655 K and 604 K, respectively, and they both have an asymmetric miscibility gap. We obtained the spatial distributions of Cd and Zn atoms in WZ and ZB Cd0.sZn0.sS alloys at different temperatures by MC simulations. We found that both WZ and ZB phases of Cdo.sZn0.sS alloy exhibit phase segregation of Cd and Zn atoms at low temperature, which is consistent with the phase diagrams.展开更多
The mechanical and thermodynamic properties of W-Ti alloys(including W_(15)Ti_(1),W_(14)Ti_(2),W_(12)Ti_(4) and W_(8)Ti_(8) alloys)were investigated by the first-principles approach based on density functional theory....The mechanical and thermodynamic properties of W-Ti alloys(including W_(15)Ti_(1),W_(14)Ti_(2),W_(12)Ti_(4) and W_(8)Ti_(8) alloys)were investigated by the first-principles approach based on density functional theory.The results indicate that W-Ti alloys except W_(8)Ti_(8) are thermodynamically stable.The modulus and hardness of W-Ti alloys are smaller than those of pure tungsten and gradually decrease with increasing Ti concentration.However,their B/G ratios and Poisson's ratios exceed those of pure tungsten,suggesting that the introduction of Ti decreases the mechanical strength while enhancing the ductility of W-Ti alloys.The thermal expansion coefficients for W-Ti alloys all surpass those of pure tungsten,indicating that the introduction of titanium exacerbates the thermal expansion behavior of W-Ti alloys.Nevertheless,elevated pressure has the capacity to suppress the thermal expansion tendencies in titanium-doped tungsten alloys.This study offers theoretical insights for the design of nuclear materials by exploring the mechanical and thermodynamic properties of W-Ti alloys.展开更多
Deeply buried mountain tunnels are often exposed to the risk of rock bursts,which always cause serious damage to the supporting structures and threaten the safety of the engineers.Due to the limited data available,a s...Deeply buried mountain tunnels are often exposed to the risk of rock bursts,which always cause serious damage to the supporting structures and threaten the safety of the engineers.Due to the limited data available,a suitable approach to predict the rockburst tendency at the preliminary stage becomes very important.In this study,an integrated methodology combining 3D initial stress inversion and rockburst tendency prediction was developed and subsequently applied to a case study of the Sangzhuling Tunnel on the Sichuan–Tibet Railway.The numerical modelling involved inverting the initial stress field using a multiple linear regression method.The tunnel excavation was simulated separately by FDM and DEM,based on a stress boundary condition from the inverted stress field.The comparative analysis demonstrates that the rockburst ratio calculated using DEM(76.70%)exhibits a slight increase compared to FDM(75.38%),and the rockburst location is consistent with the actual situation.This suggests that DEM is more suitable for simulating the stress redistribution during excavation in a jointed rock mass.The numerical simulation combined with the deviatoric stress approach effectively predicts rockburst tendency,meeting the engineering requirements.Despite its limitations,numerical simulation remains a reliable method for predicting rock bursts.展开更多
The frontal edge of the Makran accretionary wedge is characterized by the development of multiple imbricate thrust faults trending E-W and relatively parallel.However,the mechanisms underlying their formation and the ...The frontal edge of the Makran accretionary wedge is characterized by the development of multiple imbricate thrust faults trending E-W and relatively parallel.However,the mechanisms underlying their formation and the factors controlling their development remain subjects of debate.This paper,based on seismic profile analysis,employs physical simulation experiments to establish a'wedge'type subduction model.The study explores the influence of the initial wedge angle,horizontal sand layer thickness,and the presence or absence of a decollement layer on the structural styles of the thrust wedge.Experimental results indicate that as the initial wedge angle decreases from 11°to 8°,the lateral growth of the thrust wedge increases,whereas vertical growth diminishes.When the horizontal sand layer thickness is reduced from 4.5 cm to 3.0 cm,the spacing between the frontal thrusts decreases and the number of thrust faults increases.Both lateral and vertical growth are relatively reduced,resulting in a smaller thrust wedge.When a decollement layer is present,the structural style exhibits layered deformation.The decollement layer constrains the development of back thrusts and promotes the localized formation of frontal thrusts.In conclusion,the imbricate thrust faults at the frontal edge of the Makran accretionary wedge are primarily controlled by the characteristics of the wedge itself and the presence of the decollement layer.展开更多
Eutectic high entropy alloys are noted for their excellent castability and comprehensive mechanical properties.The excellent mechanical properties are closely related to the activation and evolution of deformation mec...Eutectic high entropy alloys are noted for their excellent castability and comprehensive mechanical properties.The excellent mechanical properties are closely related to the activation and evolution of deformation mechanisms at the atomic scale.In this work,AlCoCrFeNi2.1 alloy is taken as the research object.The mechanical behaviors and deformation mechanisms of the FCC and B2 single crystals with different orientations and the FCC/B2 composites with K-S orientation relationship during nanoindentation processes are systematically studied by molecular dynamics simulations.The results show that the mechanical behaviors of FCC single crystals are significantly orientation-dependent,meanwhile,the indentation force of[110]single crystal is the lowest at the elastic-plastic transition point,and that for[100]single crystal is the lowest in plastic deformation stage.Compared with FCC,the stress for B2 single crystals at the elastic-plastic transition point is higher.However,more deformation systems such as stacking faults,twins and dislocation loops are activated in FCC single crystal during the plastic deformation process,resulting in higher indentation force.For composites,the flow stress increases with the increase of B2 phase thickness during the initial stage of deformation.When indenter penetrates heterogeneous interface,the significantly increased deformation system in FCC phase leads to a significant increase in indentation force.The mechanical behaviors and deformation mechanisms depend on the component single crystal.When the thickness of the component layer is less than 15 nm,the heterogeneous interfaces fail to prevent the dislocation slip and improve the indentation force.The results will enrich the plastic deformation mechanisms of multi-principal eutectic alloys and provide guidance for the design of nanocrystalline metallic materials.展开更多
The strategic dispersion of carbon nanotubes(CNTs)within triblock copolymer matrix is key to fabricating nanocomposites with the desired electrical properties.This study investigated the self-assembly and electrical b...The strategic dispersion of carbon nanotubes(CNTs)within triblock copolymer matrix is key to fabricating nanocomposites with the desired electrical properties.This study investigated the self-assembly and electrical behavior of a polystyrene-polybutadiene-polystyrene(SBS)matrix with CNTs of different aspect ratios using hybrid particle-field molecular dynamics simulations.Structural factor analysis of the nanocomposites indicated that CNTs with higher aspect ratios promoted the transition of the SBS matrix from a bicontinuous to a lamellar phase.The resistor network algorithm method showed that the electrical conductivity of SBS and CNTs nanocomposites was influenced by the interplay between the CNTs aspect ratios,concentrations,and domain sizes of the triblock copolymer SBS.Our research sheds light on the relationship between CNTs dispersion and the electrical behavior of SBS/CNTs nanocomposites,guiding the engineering of materials to achieve desired electrical properties through the modulation of CNTs aspect ratios and tailored sizing of triblock copolymer domains.展开更多
To address the limitations of existing coupling methods in aero-engine system simulation,which fail to adaptively adjust iterative parameters and coupling relationships,which can result in low efficiency and in⁃stabil...To address the limitations of existing coupling methods in aero-engine system simulation,which fail to adaptively adjust iterative parameters and coupling relationships,which can result in low efficiency and in⁃stability,this study introduces a‘Dynamic Event-Driven Co-Simulation’algorithm integrated with decision tree algorithms.This algorithm separates the overall coupling relationships and the main solver from the primary mod⁃el,utilizing a dynamic event monitoring module to adaptively adjust simulation strategies,including iteration pa⁃rameters,coupling relationships,and convergence criteria.This facilitates efficient adaptive simulations of dy⁃namic events while balancing solution accuracy and computational efficiency.The research focuses on a twinshaft turbofan engine,establishing six system-level models that encompass overall performance and various sub⁃systems based on three coupling methods,along with a multidisciplinary multi-fidelity simulation framework in⁃corporating a 3D CFD nozzle model.The study tests both model exchange and coupled simulation methods under a 14 s transient acceleration and deceleration scenario.In a 100%throttle condition,a high-fidelity nozzle model is used to analyze the sensitivity of different convergence criteria on computational efficiency and accuracy.Re⁃sults indicate that the accuracy and efficiency achieved with this method are comparable to those of PROOSIS soft⁃ware(18 s and 35 s,respectively),while being 71%more efficient than Simulink software(62 s and 120 s,re⁃spectively).Furthermore,appropriately relaxing the convergence criteria for the 0D model(from 10-6 to 10-4)while enhancing those for the 3D model(from 3000 steps to 6000 steps)can effectively balance computational accuracy and efficiency.展开更多
Accurately determining the Flade potential(E_(Flade))is of significant importance in the design of novel corrosion-resisting alloys.However,due to the complex nature of the E_(Flade)influenced by several factors inclu...Accurately determining the Flade potential(E_(Flade))is of significant importance in the design of novel corrosion-resisting alloys.However,due to the complex nature of the E_(Flade)influenced by several factors including compositions of the alloys and corrosive solutions,there is currently a lack of truly predictive ab initio model.Here,we established the critical potential condition required for passivation in acidic solutions containing chloride ions(Cl^(-))by developing an ab initio model that incorporates the potential drop from the metal electrode to the solution,considering tunneling of electrons at metal/film interface,breakdown of the film,and electrochemical adsorption reactions at film/solution interface.These parameters were derived from the work function of the alloy substrate and passivation film,the band gap of the passivation film,and the Gibbs free energy of adsorption on the passivation film,all of which can be obtainable from first-principles calculations.This theoretical model has been successfully validated for alloyed stainless steel,exhibiting a remarkable agreement with experimental results.Importantly,enabled by the model,we have identified several alloying elements(i.e.,Ta,W,Os,and Ir)that can effec-tively lower the EFlade of the stainless steel.This work constitutes an important step forward in modeling complex passivation behaviors from first-principles,providing a useful tool for the design of corrosion-resisting alloys.展开更多
Based on first-principles calculation framework,the surface model,anodic dissolution,cathodic oxygen absorption reaction,and other related electrochemical corrosion models of Fe-Ce system were constructed,and the infl...Based on first-principles calculation framework,the surface model,anodic dissolution,cathodic oxygen absorption reaction,and other related electrochemical corrosion models of Fe-Ce system were constructed,and the influencing mechanism Ce doping on the corrosion resistance of Fe-Ce system in the Cl medium environment was analyzed.The results show that Ce doping on the first surface and subsurface inhibits the ionization of Fe atoms and greatly promotes the repassivation process of Fe matrix.Ce doping on the first layer is conducive to preventing the detachment of surface Fe atoms from Fe matrix and delaying the occurrence of corrosion.Ce atoms in the subsurface effectively increase the difficulty of Fe atoms detaching from the matrix at high Cl concentrations.When O diffusion is the controlling link of oxygen absorption reaction,Ce doping has no effects on the reaction rate of cathodic oxygen absorption.Ce doping enhances the electrochemical stability of Fe(100)1and reduces the anodic dissolution rate of Fe matrix,thereby improving its corrosion resistance.展开更多
The vacuum reactive wetting and brazing of Er_(2)Si_(2)O_(7)/MoSi_(2) coatings were investigated using a (CoFeNiCrMn)_(88)Nb_(12) high-entropy alloy (HEA) brazing filler. The microstructural evolution and wettability ...The vacuum reactive wetting and brazing of Er_(2)Si_(2)O_(7)/MoSi_(2) coatings were investigated using a (CoFeNiCrMn)_(88)Nb_(12) high-entropy alloy (HEA) brazing filler. The microstructural evolution and wettability of the HEA filler were analyzed, with particular attention to the surface energy, interfacial stability, and electronic properties of the HEA filler/rare earth silicate coating system, as determined by density functional theory (DFT). As Nb diffused into the interface and the ErNbO_(4) phase formed, the wetting angle gradually decreased to 23.12° The effective wetting and spreading of the HEA brazing filler on the rare earth silicate coating surface are strongly correlated with the formation of the ErNbO_(4) phase at the interface. Furthermore, DFT calculations reveal that the interfacial bonding energy between the BCC' and FCC' phases and the ErNbO_(4) phase, after the wetting reaction, is significantly higher than the bonding energy between the initial filler and Er_(2)Si_(2)O_(7). This finding suggests that the formation of the ErNbO_(4) phase improves the wetting and spreading behavior of the filler.展开更多
The adsorption properties of a magnesium porphyrin(MgP)molecule on Au(111)surface covered with up to three lay-ers of sodium chloride(NaCl)were investigated by means of first-principles calculations.The most stable ad...The adsorption properties of a magnesium porphyrin(MgP)molecule on Au(111)surface covered with up to three lay-ers of sodium chloride(NaCl)were investigated by means of first-principles calculations.The most stable adsorption configuration of MgP on the NaCl/Au(111)heterosurfaces was found to be at the Cl-top site with a 20°angle between the[110]lattice direction of NaCl and the Mg–N bond of the molecule.Compared with MgP molecule adsorbed on bare Au(111),the inclusion of NaCl lay-ers can lead to a significant decrease in the adsorption energy of the MgP molecule.The exis-tence of NaCl layers also reduced the charge transfer between the molecule and the surface.For heterosurfaces with two or three monolayers of NaCl,the charge transfer was almost com-pletely suppressed.The obtained partial density of states(PDOS)showed that hybridization between the electronic structures of the adsorbed MgP molecule and the metal surface can be significantly suppressed when NaCl layers were added.For the heterosurface with three lay-ers of NaCl,the PDOS around the Fermi level was almost identical with that of the free molecule,suggesting the electronic structure of the MgP molecule was nicely preserved.Influ-ence of the NaCl layers on the electronic structure of the MgP molecule was mainly found for molecular orbitals(MOs)away from the Fermi level as a result of the large band gap of the NaCl layers.展开更多
Control of hyperfine interaction strength of shallow donors in Si is one of the central issues in realizing Kane quantum computers.First-principles calculations on the hyperfine Stark shift of shallow donors are chall...Control of hyperfine interaction strength of shallow donors in Si is one of the central issues in realizing Kane quantum computers.First-principles calculations on the hyperfine Stark shift of shallow donors are challenging since large supercells are needed to accommodate the delocalized donor wave functions.In this work,we investigated the hyperfine Stark shift and its strain tunability for shallow donors P and As in Si using the potential patching method based on first-principles density functional theory calculations.The good agreement between our calculations and experimental results confirms that the potential patching method is a feasible and accurate first-principles approach for studying wave-function-related properties of shallow impurities,such as the Stark shift parameter.It is further shown that the application of strain expands the range of hyperfine Stark shift and helps improve the response of shallow donor based qubit gates.The results could be useful for developing quantum computing architectures based on shallow donors in Si.展开更多
文摘Direct conversion of syngas to light olefins(STO)on bifunctional catalysts has garnered significant attention,yet a comprehensive understanding of the reaction pathway and reaction kinetics remains elusive.Herein,we theoretically addressed the kinetics of the direct STO reaction on typical ZnAl_(2)O_(4)/zeolite catalysts by establishing a complete reaction network,consisting of methanol synthesis and conversion,water gas shift(WGS)reaction,olefin hydrogenation,and other relevant steps.The WGS reaction occurs very readily on ZnAl_(2)O_(4) surface whereas which is less active towards alkane formation via olefin hydrogenation,and the latter can be attributed to the characteristics of the H_(2) heterolytic activation and the weak polarity of olefins.The driving effect of zeolite component towards CO conversion was demonstrated by microkinetic simulations,which is sensitive to reaction conditions like space velocity and reaction temperature.Under a fixed ratio of active sites between oxide and zeolite components,the concept of the“impossible trinity”of high CO conversion,high olefin selectivity,and high space velocity can thus be manifested.This work thus provides a comprehensive kinetic picture on the direct STO conversion,offering valuable insights for the design of each component of bifunctional catalysts and the optimization of reaction conditions.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.22473010,22303114,and 12474372)the Fundamental Research Funds for the Central Universities,Jilin University,the National Key Research and Development Program of China(Grant No.SQ2023YFB2805600)+4 种基金the Natural Science Foundation of Beijing Municipality(Grant No.Z210004)the Fund from the State Key Laboratory of Information Photonics and Optical Communications(Grant No.IPOC2021ZT01)Beijing Nova Program from Beijing Municipal Science and Technology Commission(Grant No.20230484433)Beijing University of Posts and Telecommunications Excellent Ph.D.Students Foundation(Grant No.CX20241078)Beijing Natural Science Foundation(Undergraduate Program)(Grant No.QY24218)。
文摘The development of high-performance solid electrolytes is pivotal for advancing solid-state battery technologies.In this work,we design an oxysulfide-based solid electrolyte Na MgPO_(3)S by combining bond valence theory and density functional theory calculations.The material features a wide band gap of 4.0 eV and a considerable reduced Na^(+)migration barrier of 0.44 eV,a 1.26-eV decrease compared to pristine Na MgPO_(4)(~1.70 eV).Ab initio molecular dynamics simulations further reveal significantly enhanced ionic conductivity in the oxysulfide-based system compared to the pristine oxide structure.In addition,the calculated decomposition energy indicates that the modified material exhibits good moisture stability.Our findings suggest that sulfur-doping strategy can simultaneously achieve improved ionic conductivity and high moisture stability in oxide solid electrolytes,which could pave the way for designing high-performance solid electrolytes.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11505003 and 52325103)the Anhui Provincial Natural Science Foundation(Grant No.2108085MA25)。
文摘High-performance intelligent protective materials are vital for nuclear energy systems exposed to extreme irradiation.Among them,tungsten-based alloys demonstrate exceptional potential owing to their superior irradiation resistance.Recent experimental studies have demonstrated that W-Ta-Cr-V alloys exhibit excellent irradiation resistance under helium(He)ion irradiation.However,the underlying mechanisms,especially the migration behavior of He atoms,remain unclear.In this work,the influences of uniaxial tensile and compressive strain on He migration in W-Ta-Cr-V complex alloys have been systematically studied through first-principles calculations.Our results demonstrate that He atoms preferentially occupy the tetrahedral interstitial sites,with interstitial formation energies significantly reduced compared to pure W.The introduction of Ta,Cr,and V alloying elements markedly increases the He migration barriers,effectively suppressing He diffusion.Compressive strain increases the migration barriers,inhibiting He bubbles nucleation and growth,while tensile strain decreases the barriers,facilitating bubble formation.Compared to pure W,the W-Ta-Cr-V alloys exhibit both lower He interstitial formation energies and higher migration barriers,with further enhancement under compressive strain.Specifically,compressive strain of 6%increases the He migration barrier of the W-Ta-Cr-V alloy by 0.166 e V,which further widens the difference relative to pure W.These findings provide a theoretical explanation for the superior irradiation resistance of tungsten-based alloys observed experimentally and promote the understanding of irradiation damage in these alloys under strain.
基金supported by the Advanced Materials-National Science and Technology Major Project(Grant No.2025ZD0618401)the National Natural Science Foundation of China(Grant No.12504285)+1 种基金the Natural Science Foundation of Jiangsu Province(Grant No.BK20250472)NFSG grant from BITS-Pilani,Dubai campus。
文摘The rapid advancement of machine learning based tight-binding Hamiltonian(MLTB)methods has opened new avenues for efficient and accurate electronic structure simulations,particularly in large-scale systems and long-time scenarios.This review begins with a concise overview of traditional tight-binding(TB)models,including both(semi-)empirical and first-principles approaches,establishing the foundation for understanding MLTB developments.We then present a systematic classification of existing MLTB methodologies,grouped into two major categories:direct prediction of TB Hamiltonian elements and inference of empirical parameters.A comparative analysis with other ML-based electronic structure models is also provided,highlighting the advancement of MLTB approaches.Finally,we explore the emerging MLTB application ecosystem,highlighting how the integration of MLTB models with a diverse suite of post-processing tools from linear-scaling solvers to quantum transport frameworks and molecular dynamics interfaces is essential for tackling complex scientific problems across different domains.The continued advancement of this integrated paradigm promises to accelerate materials discovery and open new frontiers in the predictive simulation of complex quantum phenomena.
基金financially supported by the National Natural Science Foundation of China(Nos.52173020 and 52573023)。
文摘Vitrimers belong to a class of polymeric materials capable of bond exchange reactions,showing great promise for environmental protection and sustainable development.However,studies on the coupling mechanism between the bond exchange kinetics and segmental dynamics near the glass transition temperature(T_(g))remain scarce.Herein,we employed molecular dynamics simulations to investigate the dynamic heterogeneity of the segment motion and bond exchange in vitrimers.The simulation results revealed that the bond exchange energy barrier exerts a much stronger influence on the bond exchange kinetics than on the segmental dynamics.At lower temperatures,slower segmental relaxation further constraind the bond exchange rate.Additionally,increasing the bond exchange energy barrier markedly enhanced the dynamic heterogeneity of segment motion.A close correlation was observed between heterogeneity and bond exchange.This study elucidated the coupling mechanism between bond exchange and segmental dynamics at the molecular scale,thereby providing a theoretical basis for designing vitrimer materials with tunable dynamic properties.
基金supported by the National Natural Science Foundation of China (Grant No.11574244 for G.Y.G.)the XJTU Research Fund for AI Science (Grant No.2025YXYC011 for G.Y.G.)the Hong Kong Global STEM Professorship Scheme (for X.C.Z.)。
文摘Compared to the well-studied two-dimensional(2D)ferroelectricity,the appearance of 2D antiferroelectricity is much rarer,where local dipoles from the nonequivalent sublattices within 2D monolayers are oppositely oriented.Using NbOCl_(2) monolayer with competing ferroelectric(FE)and antiferroelectric(AFE)phases as a 2D material platform,we demonstrate the emergence of intrinsic antiferroelectricity in NbOCl_(2) monolayer under experimentally accessible shear strain,along with new functionality associated with electric field-induced AFE-to-FE phase transition.Specifically,the complex configuration space accommodating FE and AFE phases,polarization switching kinetics,and finite temperature thermodynamic properties of 2D NbOCl_(2) are all accurately predicted by large-scale molecular dynamics simulations based on deep learning interatomic potential model.Moreover,room temperature stable antiferroelectricity with low polarization switching barrier and one-dimensional collinear polarization arrangement is predicted in shear-deformed NbOCl_(2) monolayer.The transition from AFE to FE phase in 2D NbOCl_(2) can be triggered by a low critical electric field,leading to a double polarization–electric(P–E)loop with small hysteresis.A new type of optoelectronic device composed of AFE-NbOCl_(2) is proposed,enabling electric“writing”and nonlinear optical“reading”logical operation with fast operation speed and low power consumption.
基金supported by the National Natural Science Foundation of China(Grant Nos.11164014 and 11364025)Gansu Science and Technology Pillar Program,China(Grant No.1204GKCA057)
文摘The first-principles calculations based on density functional theory combined with cluster expansion techniques and Monte Carlo （MC） simulations were used to study the phase diagrams of both wurtzite （WZ） and zinc-blende （ZB） Cdl_xZnxS alloys. All formation energies are positive for WZ and ZB Cdl-xZnxS alloys, which means that the Cdl-xZnxS alloys are unstable and have a tendency to phase separation. For WZ and ZB Cdl_xZnxS alloys, the consolute temperatures are 655 K and 604 K, respectively, and they both have an asymmetric miscibility gap. We obtained the spatial distributions of Cd and Zn atoms in WZ and ZB Cd0.sZn0.sS alloys at different temperatures by MC simulations. We found that both WZ and ZB phases of Cdo.sZn0.sS alloy exhibit phase segregation of Cd and Zn atoms at low temperature, which is consistent with the phase diagrams.
基金Funded by National Key R&D Program of China(No.2021YFB3802300)the National Natural Science Foundation of China(No.52171045)the Joint Fund(No.8091B022108)。
文摘The mechanical and thermodynamic properties of W-Ti alloys(including W_(15)Ti_(1),W_(14)Ti_(2),W_(12)Ti_(4) and W_(8)Ti_(8) alloys)were investigated by the first-principles approach based on density functional theory.The results indicate that W-Ti alloys except W_(8)Ti_(8) are thermodynamically stable.The modulus and hardness of W-Ti alloys are smaller than those of pure tungsten and gradually decrease with increasing Ti concentration.However,their B/G ratios and Poisson's ratios exceed those of pure tungsten,suggesting that the introduction of Ti decreases the mechanical strength while enhancing the ductility of W-Ti alloys.The thermal expansion coefficients for W-Ti alloys all surpass those of pure tungsten,indicating that the introduction of titanium exacerbates the thermal expansion behavior of W-Ti alloys.Nevertheless,elevated pressure has the capacity to suppress the thermal expansion tendencies in titanium-doped tungsten alloys.This study offers theoretical insights for the design of nuclear materials by exploring the mechanical and thermodynamic properties of W-Ti alloys.
基金financially supported by the State Key Laboratory of Geohazard Prevention and Geoenvironment Protection(Chengdu University of Technology)(Grant No.SKLGP2020Z007)。
文摘Deeply buried mountain tunnels are often exposed to the risk of rock bursts,which always cause serious damage to the supporting structures and threaten the safety of the engineers.Due to the limited data available,a suitable approach to predict the rockburst tendency at the preliminary stage becomes very important.In this study,an integrated methodology combining 3D initial stress inversion and rockburst tendency prediction was developed and subsequently applied to a case study of the Sangzhuling Tunnel on the Sichuan–Tibet Railway.The numerical modelling involved inverting the initial stress field using a multiple linear regression method.The tunnel excavation was simulated separately by FDM and DEM,based on a stress boundary condition from the inverted stress field.The comparative analysis demonstrates that the rockburst ratio calculated using DEM(76.70%)exhibits a slight increase compared to FDM(75.38%),and the rockburst location is consistent with the actual situation.This suggests that DEM is more suitable for simulating the stress redistribution during excavation in a jointed rock mass.The numerical simulation combined with the deviatoric stress approach effectively predicts rockburst tendency,meeting the engineering requirements.Despite its limitations,numerical simulation remains a reliable method for predicting rock bursts.
基金the National Natural Science Foundation of China(No.42076069)。
文摘The frontal edge of the Makran accretionary wedge is characterized by the development of multiple imbricate thrust faults trending E-W and relatively parallel.However,the mechanisms underlying their formation and the factors controlling their development remain subjects of debate.This paper,based on seismic profile analysis,employs physical simulation experiments to establish a'wedge'type subduction model.The study explores the influence of the initial wedge angle,horizontal sand layer thickness,and the presence or absence of a decollement layer on the structural styles of the thrust wedge.Experimental results indicate that as the initial wedge angle decreases from 11°to 8°,the lateral growth of the thrust wedge increases,whereas vertical growth diminishes.When the horizontal sand layer thickness is reduced from 4.5 cm to 3.0 cm,the spacing between the frontal thrusts decreases and the number of thrust faults increases.Both lateral and vertical growth are relatively reduced,resulting in a smaller thrust wedge.When a decollement layer is present,the structural style exhibits layered deformation.The decollement layer constrains the development of back thrusts and promotes the localized formation of frontal thrusts.In conclusion,the imbricate thrust faults at the frontal edge of the Makran accretionary wedge are primarily controlled by the characteristics of the wedge itself and the presence of the decollement layer.
基金supported by the Natural Science Foundation of Hebei Province(E2024209052)the Youth Scholars Promotion Plan of North China University of Science and Technology(QNTJ202307).
文摘Eutectic high entropy alloys are noted for their excellent castability and comprehensive mechanical properties.The excellent mechanical properties are closely related to the activation and evolution of deformation mechanisms at the atomic scale.In this work,AlCoCrFeNi2.1 alloy is taken as the research object.The mechanical behaviors and deformation mechanisms of the FCC and B2 single crystals with different orientations and the FCC/B2 composites with K-S orientation relationship during nanoindentation processes are systematically studied by molecular dynamics simulations.The results show that the mechanical behaviors of FCC single crystals are significantly orientation-dependent,meanwhile,the indentation force of[110]single crystal is the lowest at the elastic-plastic transition point,and that for[100]single crystal is the lowest in plastic deformation stage.Compared with FCC,the stress for B2 single crystals at the elastic-plastic transition point is higher.However,more deformation systems such as stacking faults,twins and dislocation loops are activated in FCC single crystal during the plastic deformation process,resulting in higher indentation force.For composites,the flow stress increases with the increase of B2 phase thickness during the initial stage of deformation.When indenter penetrates heterogeneous interface,the significantly increased deformation system in FCC phase leads to a significant increase in indentation force.The mechanical behaviors and deformation mechanisms depend on the component single crystal.When the thickness of the component layer is less than 15 nm,the heterogeneous interfaces fail to prevent the dislocation slip and improve the indentation force.The results will enrich the plastic deformation mechanisms of multi-principal eutectic alloys and provide guidance for the design of nanocrystalline metallic materials.
基金financially supported by the National Natural Science Foundation of China(Nos.52273019,62173065,22133002,22273031,and 12274056)Fundamental Research Funds for the Central Universities(No.04442024074)+2 种基金NationalKey R&D Program of China(No.2022YFB3707300)Beijing Natural Science Foundation(No.4242040)Scientific Research Funds Project of Liaoning Provincial Department of Education(No.LJKZ0034)。
文摘The strategic dispersion of carbon nanotubes(CNTs)within triblock copolymer matrix is key to fabricating nanocomposites with the desired electrical properties.This study investigated the self-assembly and electrical behavior of a polystyrene-polybutadiene-polystyrene(SBS)matrix with CNTs of different aspect ratios using hybrid particle-field molecular dynamics simulations.Structural factor analysis of the nanocomposites indicated that CNTs with higher aspect ratios promoted the transition of the SBS matrix from a bicontinuous to a lamellar phase.The resistor network algorithm method showed that the electrical conductivity of SBS and CNTs nanocomposites was influenced by the interplay between the CNTs aspect ratios,concentrations,and domain sizes of the triblock copolymer SBS.Our research sheds light on the relationship between CNTs dispersion and the electrical behavior of SBS/CNTs nanocomposites,guiding the engineering of materials to achieve desired electrical properties through the modulation of CNTs aspect ratios and tailored sizing of triblock copolymer domains.
文摘To address the limitations of existing coupling methods in aero-engine system simulation,which fail to adaptively adjust iterative parameters and coupling relationships,which can result in low efficiency and in⁃stability,this study introduces a‘Dynamic Event-Driven Co-Simulation’algorithm integrated with decision tree algorithms.This algorithm separates the overall coupling relationships and the main solver from the primary mod⁃el,utilizing a dynamic event monitoring module to adaptively adjust simulation strategies,including iteration pa⁃rameters,coupling relationships,and convergence criteria.This facilitates efficient adaptive simulations of dy⁃namic events while balancing solution accuracy and computational efficiency.The research focuses on a twinshaft turbofan engine,establishing six system-level models that encompass overall performance and various sub⁃systems based on three coupling methods,along with a multidisciplinary multi-fidelity simulation framework in⁃corporating a 3D CFD nozzle model.The study tests both model exchange and coupled simulation methods under a 14 s transient acceleration and deceleration scenario.In a 100%throttle condition,a high-fidelity nozzle model is used to analyze the sensitivity of different convergence criteria on computational efficiency and accuracy.Re⁃sults indicate that the accuracy and efficiency achieved with this method are comparable to those of PROOSIS soft⁃ware(18 s and 35 s,respectively),while being 71%more efficient than Simulink software(62 s and 120 s,re⁃spectively).Furthermore,appropriately relaxing the convergence criteria for the 0D model(from 10-6 to 10-4)while enhancing those for the 3D model(from 3000 steps to 6000 steps)can effectively balance computational accuracy and efficiency.
基金supported by the Natural Science Foundation Project of Liaoning province(Nos.2023-MS-017 and 2022-MS-005)the National Science and Technology Major Project(No.J2019-VI-0019-0134)+2 种基金the National Natural Science Foundation of China(Nos.52201028 and 52188101)the Project funded by China Postdoctoral Science Foundation(No.2021M700153)the Special Projects of the Central Government in Guidance of Local Science and Technology Development(No.2024010859-JH6/1006).
文摘Accurately determining the Flade potential(E_(Flade))is of significant importance in the design of novel corrosion-resisting alloys.However,due to the complex nature of the E_(Flade)influenced by several factors including compositions of the alloys and corrosive solutions,there is currently a lack of truly predictive ab initio model.Here,we established the critical potential condition required for passivation in acidic solutions containing chloride ions(Cl^(-))by developing an ab initio model that incorporates the potential drop from the metal electrode to the solution,considering tunneling of electrons at metal/film interface,breakdown of the film,and electrochemical adsorption reactions at film/solution interface.These parameters were derived from the work function of the alloy substrate and passivation film,the band gap of the passivation film,and the Gibbs free energy of adsorption on the passivation film,all of which can be obtainable from first-principles calculations.This theoretical model has been successfully validated for alloyed stainless steel,exhibiting a remarkable agreement with experimental results.Importantly,enabled by the model,we have identified several alloying elements(i.e.,Ta,W,Os,and Ir)that can effec-tively lower the EFlade of the stainless steel.This work constitutes an important step forward in modeling complex passivation behaviors from first-principles,providing a useful tool for the design of corrosion-resisting alloys.
基金Project supported by the National Natural Science Foundation of China(52364044,52204364)Scientific Research Special Project for First-Class Disciplines of Education Department of Inner Mongolia Autonomous Region(YLXKZX-NKD-001,YLXKZX-NKD-011)Basic Scientific Research Business Expenses of Colleges and Universities of Inner Mongolia Autonomous Region(2023QNJS011)。
文摘Based on first-principles calculation framework,the surface model,anodic dissolution,cathodic oxygen absorption reaction,and other related electrochemical corrosion models of Fe-Ce system were constructed,and the influencing mechanism Ce doping on the corrosion resistance of Fe-Ce system in the Cl medium environment was analyzed.The results show that Ce doping on the first surface and subsurface inhibits the ionization of Fe atoms and greatly promotes the repassivation process of Fe matrix.Ce doping on the first layer is conducive to preventing the detachment of surface Fe atoms from Fe matrix and delaying the occurrence of corrosion.Ce atoms in the subsurface effectively increase the difficulty of Fe atoms detaching from the matrix at high Cl concentrations.When O diffusion is the controlling link of oxygen absorption reaction,Ce doping has no effects on the reaction rate of cathodic oxygen absorption.Ce doping enhances the electrochemical stability of Fe(100)1and reduces the anodic dissolution rate of Fe matrix,thereby improving its corrosion resistance.
基金support from the National Natural Science Foundation of China(No.52374402)the National Key Research and Development Program(No.2022YFB3402200)+2 种基金the National Science and Technology Major Project(No.J2022-VII-0003-0045)the Project of Key areas of innovation team in Shaanxi Province(No.2024RS-CXTD-20)the Innovation Foundation for Doctor Dissertation of Northwestern Polytechnical University(No.CX2024055).
文摘The vacuum reactive wetting and brazing of Er_(2)Si_(2)O_(7)/MoSi_(2) coatings were investigated using a (CoFeNiCrMn)_(88)Nb_(12) high-entropy alloy (HEA) brazing filler. The microstructural evolution and wettability of the HEA filler were analyzed, with particular attention to the surface energy, interfacial stability, and electronic properties of the HEA filler/rare earth silicate coating system, as determined by density functional theory (DFT). As Nb diffused into the interface and the ErNbO_(4) phase formed, the wetting angle gradually decreased to 23.12° The effective wetting and spreading of the HEA brazing filler on the rare earth silicate coating surface are strongly correlated with the formation of the ErNbO_(4) phase at the interface. Furthermore, DFT calculations reveal that the interfacial bonding energy between the BCC' and FCC' phases and the ErNbO_(4) phase, after the wetting reaction, is significantly higher than the bonding energy between the initial filler and Er_(2)Si_(2)O_(7). This finding suggests that the formation of the ErNbO_(4) phase improves the wetting and spreading behavior of the filler.
基金supported by the National Natural Science Foundation of China(No.22373084,No.62201494)Hebei Natural Science Foundation(B2022203007)the Cultivation Project for Basic Research and Innovation of Yanshan University(2024LGZD002).
文摘The adsorption properties of a magnesium porphyrin(MgP)molecule on Au(111)surface covered with up to three lay-ers of sodium chloride(NaCl)were investigated by means of first-principles calculations.The most stable adsorption configuration of MgP on the NaCl/Au(111)heterosurfaces was found to be at the Cl-top site with a 20°angle between the[110]lattice direction of NaCl and the Mg–N bond of the molecule.Compared with MgP molecule adsorbed on bare Au(111),the inclusion of NaCl lay-ers can lead to a significant decrease in the adsorption energy of the MgP molecule.The exis-tence of NaCl layers also reduced the charge transfer between the molecule and the surface.For heterosurfaces with two or three monolayers of NaCl,the charge transfer was almost com-pletely suppressed.The obtained partial density of states(PDOS)showed that hybridization between the electronic structures of the adsorbed MgP molecule and the metal surface can be significantly suppressed when NaCl layers were added.For the heterosurface with three lay-ers of NaCl,the PDOS around the Fermi level was almost identical with that of the free molecule,suggesting the electronic structure of the MgP molecule was nicely preserved.Influ-ence of the NaCl layers on the electronic structure of the MgP molecule was mainly found for molecular orbitals(MOs)away from the Fermi level as a result of the large band gap of the NaCl layers.
基金supported by the National Natural Science Foun-dation of China(Grant Nos.12393831 and 12088101).
文摘Control of hyperfine interaction strength of shallow donors in Si is one of the central issues in realizing Kane quantum computers.First-principles calculations on the hyperfine Stark shift of shallow donors are challenging since large supercells are needed to accommodate the delocalized donor wave functions.In this work,we investigated the hyperfine Stark shift and its strain tunability for shallow donors P and As in Si using the potential patching method based on first-principles density functional theory calculations.The good agreement between our calculations and experimental results confirms that the potential patching method is a feasible and accurate first-principles approach for studying wave-function-related properties of shallow impurities,such as the Stark shift parameter.It is further shown that the application of strain expands the range of hyperfine Stark shift and helps improve the response of shallow donor based qubit gates.The results could be useful for developing quantum computing architectures based on shallow donors in Si.
