Probing the electric double layer structure at nitrogen-doped graphite electrodes by constant-potential molecular dynamics simulations

Electric double layer(EDL)is a critical topic in electrochemistry and largely determines the working performance of lithium batteries.However,atomic insights into the EDL structures on heteroatom-modified graphite anodes and EDL evolution with electrode potential are very lacking.Herein,a constant-potential molecular dynamics(CPMD)method is proposed to probe the EDL structure under working conditions,taking N-doped graphite electrodes and carbonate electrolytes as an example.An interface model was developed,incorporating the electrode potential and atom electronegativities.As a result,an insightful atomic scenario for the EDL structure under varied electrode potentials has been established,which unveils the important role of doping sites in regulating both the EDL structures and the following electrochemical reactions at the atomic level.Specifically,the negatively charged N atoms repel the anions and adsorb Li~+at high and low potentials,respectively.Such preferential adsorption suggests that Ndoped graphite can promote Li~+desolvation and regulate the location of Li~+deposition.This CPMD method not only unveils the mysterious function of N-doping from the viewpoint of EDL at the atomic level but also applies to probe the interfacial structure on other complicated electrodes.

Three-dimensional thermohaline structure estimation derived from HY-2 satellite data over the Maritime Silk Road and its applications

Estimated ocean subsurface fields derived from satellite observations provide potential data sources for operational marine environmental monitoring and prediction systems.This study employs a statistic regression reconstruction method,in combination with domestic autonomous sea surface height and sea surface temperature observations from the Haiyang-2(HY-2)satellite fusion data,to establish an operational quasi-realtime three-dimensional(3D)temperature and salinity products over the Maritime Silk Road.These products feature a daily temporal resolution and a spatial resolution of 0.25°×0.25°and exhibit stability and continuity.We have demonstrated the accuracy of the reconstructed thermohaline fields in capturing the 3D thermohaline variations through comprehensive statistical evaluations,after comparing them against Argo observations and ocean analysis data from 2022.The results illustrate that the reconstructed fields effectively represent seasonal variations in oceanic subsurface structures,along with structural changes resulting from mesoscale processes,and the upper ocean’s responses to tropical cyclones.Furthermore,the incorporation of HY-2 satellite observations notably enhances the accuracy of temperature and salinity reconstructions in the Northwest Pacific Ocean and marginally improves salinity reconstruction accuracy in the North Indian Ocean when compared to the World Ocean Atlas 2018 monthly climatology thermohaline fields.As a result,the reconstructed product holds promise for providing quasi-real-time 3D temperature and salinity field information to facilitate fast decisionmaking during emergencies,and also offers foundational thermohaline fields for operational ocean reanalysis and forecasting systems.These contributions enhance the safety and stability of ocean subsurface activities and navigation.

Influence of layer thickness on formation quality,microstructure,mechanical properties,and corrosion resistance of WE43 magnesium alloy fabricated by laser powder bed fusion

Laser powder bed fusion(L-PBF)of Mg alloys has provided tremendous opportunities for customized production of aeronautical and medical parts.Layer thickness(LT)is of great significance to the L-PBF process but has not been studied for Mg alloys.In this study,WE43 Mg alloy bulk cubes,porous scaffolds,and thin walls with layer thicknesses of 10,20,30,and 40μm were fabricated.The required laser energy input increased with increasing layer thickness and was different for the bulk cubes and porous scaffolds.Porosity tended to occur at the connection joints in porous scaffolds for LT40 and could be eliminated by reducing the laser energy input.For thin wall parts,a large overhang angle or a small wall thickness resulted in porosity when a large layer thicknesses was used,and the porosity disappeared by reducing the layer thickness or laser energy input.A deeper keyhole penetration was found in all occasions with porosity,explaining the influence of layer thickness,geometrical structure,and laser energy input on the porosity.All the samples achieved a high fusion quality with a relative density of over 99.5%using the optimized laser energy input.The increased layer thickness resulted to more precipitation phases,finer grain sizes and decreased grain texture.With the similar high fusion quality,the tensile strength and elongation of bulk samples were significantly improved from 257 MPa and 1.41%with the 10μm layer to 287 MPa and 15.12%with the 40μm layer,in accordance with the microstructural change.The effect of layer thickness on the compressive properties of porous scaffolds was limited.However,the corrosion rate of bulk samples accelerated with increasing the layer thickness,mainly attributed to the increased number of precipitation phases.

Reversible Mn^(2+)/Mn^(4+)double-electron redox in P3-type layer-structured sodium-ion cathode

The balance between cationic redox and oxygen redox in layer-structured cathode materials is an important issue for sodium batteries to obtain high energy density and considerable cycle stability.Oxygen redox can contribute extra capacity to increase energy density,but results in lattice instability and capacity fading caused by lattice oxygen gliding and oxygen release.In this work,reversible Mn^(2+)/Mn^(4+)redox is realized in a P3-Na_(0.65)Li_(0.2)Co_(0.05)Mn_(0.75)O_(2)cathode material with high specific capacity and structure stability via Co substitution.The contribution of oxygen redox is suppressed significantly by reversible Mn^(2+)/Mn^(4+)redox without sacrificing capacity,thus reducing lattice oxygen release and improving the structure stability.Synchrotron X-ray techniques reveal that P3 phase is well maintained in a wide voltage window of 1.5-4.5 V vs.Na^(+)/Na even at 10 C and after long-term cycling.It is disclosed that charge compensation from Co/Mn-ions contributes to the voltage region below 4.2 V and O-ions contribute to the whole voltage range.The synergistic contributions of Mn^(2+)/Mn^(4+),Co^(2+)/Co^(3+),and O^(2-)/(O_n)^(2-)redox in P3-Na_(0.65)Li_(0.2)Co_(0.05)Mn_(0.75)O_(2)lead to a high reversible capacity of 215.0 m A h g^(-1)at 0.1 C with considerable cycle stability.The strategy opens up new opportunities for the design of high capacity cathode materials for rechargeable batteries.

Case Studies of the Microphysical and Kinematic Structure of Summer Mesoscale Precipitation Clouds over the Eastern Tibetan Plateau

Three cases of microphysical characteristics and kinematic structures in the negative temperature region of summer mesoscale cloud systems over the eastern Tibetan Plateau(TP)were investigated using X-band dual-polarization radar.The time-height series of radar physical variables and mesoscale horizontal divergence δderived by quasi-vertical profiles(QVPs)indicated that the dendritic growth layer(DGL,-20°C to-10°C)was ubiquitous,with large-value zones of K_(DP)(specific differential phase),Z_(DR)(differential reflectivity),or both,and corresponded to various dynamic fields(ascent or descent).Ascents in the DGL of cloud systems with vigorous vertical development were coincident with large-value zones of Z_(DR),signifying ice crystals with a large axis ratio,but with no obvious large values of K_(DP),which differs from previous findings.It is speculated that ascent in the DGL promoted ice crystals to undergo further growth before sinking.If there was descent in the DGL,a high echo top corresponded to large values of K_(DP),denoting a large number concentration of ice crystals;but with the echo top descending,small values of K_(DP) formed.This is similar to previous results and reveals that a high echo top is conducive to the generation of ice crystals.When ice particles fall to low levels(-10℃ to 0℃),they grow through riming,aggregation,or deposition,and may not be related to the kinematic structure.It is important to note that this study was only based on a limited number of cases and that further research is therefore needed.

In-situ additive manufacturing of high strength yet ductility titanium composites with gradient layered structure using N_(2)

It has always been challenging work to reconcile the contradiction between the strength and plasticity of titanium materials.Laser powder bed fusion(LPBF) is a convenient method to fabricate innovative composites including those inspired by gradient layered materials.In this work,we used LPBF to selectively prepare Ti N/Ti gradient layered structure(GLSTi)composites by using different N_(2)–Ar ratios during the LPBF process.We systematically investigated the mechanisms of in-situ synthesis Ti N,high strength and ductility of GLSTi composites using microscopic analysis,TEM characterization,and tensile testing with digital image correlation.Besides,a digital correspondence was established between the N_(2) concentration and the volume fraction of LPBF in-situ synthesized Ti N.Our results show that the GLSTi composites exhibit superior mechanical properties compared to pure titanium fabricated by LPBF under pure Ar.Specifically,the tensile strength of GLSTi was more than 1.5times higher than that of LPBF-formed pure titanium,reaching up to 1100 MPa,while maintaining a high elongation at fracture of 17%.GLSTi breaks the bottleneck of high strength but low ductility exhibited by conventional nanoceramic particle-strengthened titanium matrix composites,and the hetero-deformation induced strengthening effect formed by the Ti N/Ti layered structure explained its strength-plasticity balanced principle.The microhardness exhibits a jagged variation of the relatively low hardness of 245 HV0.2 for the pure titanium layer and a high hardness of 408 HV0.2 for the N_(2) in-situ synthesis layer.Our study provides a new concept for the structure-performance digital customization of 3D-printed Ti-based composites.

A Simplified Method for the Stress Analysis of Underground Transfer Structures Crossing Multiple Subway Tunnels

According to the design specifications,the construction of extended piles involves traversing the tunnel’s upper region and extending to the underlying rock layer.To address this challenge,a subterranean transfer structure spanning multiple subway tunnels was proposed.Deliberating on the function of piles in the transfer structure as springs with axial and bending stiffness,and taking into account the force balance and deformation coordination conditions of beams and plates within the transfer structure,we established a simplified mechanical model that incorporates soil stratification by combining it with the Winkler elastic foundation beam model.The resolved established simplifiedmechanicalmodel employed finite difference technology and the Newton-Simpsonmethod,elucidating the mechanical mechanism of the transfer structure.The research findings suggest that the load carried by the upper structural columns can be transferred to the pile foundation beneath the beams through the transfer structure,subsequently reaching the deep soil layer and ensuring minimal impact on adjacent tunnels.The established simplified analysis method can be used for stress analysis of the transfer structure,concurrently considering soil stratification,pile foundation behavior,and plate action.The pile length,pile section size,and beam section size within the transfer structure should account for the characteristics of the upper load,ensuring an even distribution of the beam bending moment.

Influence of topography on the fine structures of stratospheric gravity waves:An analysis using COSMIC-2 temperature data

We derive the potential energy of gravity waves(GWs)in the upper troposphere and stratosphere at 45°S-45°N from December 2019 to November 2022 by using temperature profiles retrieved from the Constellation Observing System for Meteorology,Ionosphere,and Climate-2(COSMIC-2)satellite.Owing to the dense sampling of COSMIC-2,in addition to the strong peaks of gravity wave potential energy(GWPE)above the Andes and Tibetan Plateau,we found weak peaks above the Rocky,Atlas,Caucasus,and Tianshan Mountains.The land-sea contrast is responsible for the longitudinal variations of the GWPE in the lower and upper stratosphere.At 40°N/S,the peaks were mainly above the topographic regions during the winter.At 20°N/S,the peaks were a slight distance away from the topographic regions and might be the combined effect of nontopographic GWs and mountain waves.Near the Equator,the peaks were mainly above the regions with the lowest sea level altitude and may have resulted from convection.Our results indicate that even above the local regions with lower sea level altitudes compared with the Andes and Tibetan Plateau,the GWPE also exhibits fine structures in geographic distributions.We found that dissipation layers above the tropopause jet provide the body force to generate secondary waves in the upper stratosphere,especially during the winter months of each hemisphere and at latitudes of greater than 20°N/S.

Impact of Dietary Lactobacillus plantarum Postbiotics on the Performance of Layer Hens under Heat Stress Conditions

This experiment was conducted to determine the performance of heat-stressed layers fed a diet containing the probiotic Lactobacillus plantarum RS5 or its products of fermentation (postbiotics). Twenty-week-old Isa White layers, were subdivided into six treatments of 32 individually caged birds. Half of the birds were reared under regular temperature conditions, while the other half was subjected to cyclic daily heat stress. Layers were offered one of three diets: 1) Control;2) Control + Lactobacillus plantarum RS5 probiotic;3) Control + Lactobacillus plantarum RS5 postbiotics. Birds were tested for performance and visceral organ development for 5 months. Heat stress negatively affected the birds’ feed intake, egg weight, shell weight percentage, Haugh unit, shell thickness, yolk color, body weight and spleen weight percentage. Postbiotics significantly increased egg production (p < 0.05) in comparison to the control and the probiotic fed group (94.8% vs 92.6% vs 93.1%, respectively). Birds under probiotic or postbiotic diet showed a significantly higher (p < 0.05) feed intake and egg weight, although the probiotic had a more pronounced and gradual effect. Specific gravity, yolk weight percentage and shell thickness didn’t show differences among dietary groups. The Haugh Unit was significantly higher (p < 0.05) in probiotic group which also showed a significantly lower yolk color index (p < 0.05). The different feed treatments did not impact the bird’s viscera weight percentage, except for the ileum that was significantly lower (p < 0.05) under postbiotic supplementation. Both probiotics and postbiotics could be used as a potential growth promoters and might alleviate heat stress impact in poultry industry.

Layered metastructure containing freely-designed local resonators for wave attenuation

Combining periodic layered structure with three-dimensional cylindrical local resonators,a hybrid metastructure with improved wave isolation ability was designed and investigated through theoretical and numerical approaches.The metastructure is composed of periodic rubber layers and concrete layers embedded with three-dimensional resonators,which can be freely designed with multi local resonant frequencies to attenuate vibrations at required frequencies and widen the attenuation bandgap.The metastructure can also effectively attenuate seismic responses.Compared with layered rubber-based structures,the metastructure has more excellent wave attenuation effects with greater attenuation and wider bandgap.

Sedimentary Environment and Mineralization of the Black Shale Polymetallic Layer in the Early Cambrian,SW China:Constraints from in situ LA-ICP-MS Analysis of Pyrite

A polymetallic layer is usually developed at the bottom of the early Cambrian black shale in Guizhou Province.The mineral that makes up the polymetallic layer is related to the sedimentary facies.To analyze the differentiation mechanism between polymetallic deposits(Ni-Mo and V),the Zhijin Gezhongwu profile located in the outer shelf and the Sansui Haishan V deposit located in the lower slope are selected to study the in situ sulfur isotopes and trace elements of pyrite.The results show that δ^(34)S values of pyrite vary widely from−7.8‰to 28‰in the Gezhongwu profile,while the δ^(34)S values are relatively uniform(from 27.8‰to 38.4‰)in the Haishan profile.The isotopic S composition is consistent with the transition that occurs in the sedimentary phase from the shelf to the deep sea on the transgressive Yangtze platform;this indicates that the δ^(34)SO_(4)^(2−)values in seawater must be differently distributed in depositional environments.The sulfur in the Ni-Mo layer is produced after the mixing of seawater and hydrothermal fluid,while the V layer mainly originates from seawater.Overall,the Ni-Mo and V deposits have been differentiated primarily on the basis of the combined effect of continental weathering and hydrothermal fluid.

OptoGPT: A foundation model for inverse design in optical multilayer thin film structures

Optical multilayer thin film structures have been widely used in numerous photonic applications.However,existing inverse design methods have many drawbacks because they either fail to quickly adapt to different design targets,or are difficult to suit for different types of structures,e.g.,designing for different materials at each layer.These methods also cannot accommodate versatile design situations under different angles and polarizations.In addition,how to benefit practical fabrications and manufacturing has not been extensively considered yet.In this work,we introduce OptoGPT(Opto Generative Pretrained Transformer),a decoder-only transformer,to solve all these drawbacks and issues simultaneously.

The Clumpy Structure of Five Star-bursting Dwarf Galaxies in the MaNGA Survey

The star-forming clumps in star-bursting dwarf galaxies provide valuable insights into understanding the evolution of dwarf galaxies.In this paper,we focus on five star-bursting dwarf galaxies featuring off-centered clumps in the Mapping Nearby Galaxies at Apache Point Observatory survey.Using the stellar population synthesis software Fitting Analysis using Differential evolution Optimization,we obtain the spatially resolved distribution of the star formation history,which allows us to construct the g-band images of the five galaxies at different ages.These images can help us to probe the evolution of the morphological structures of these galaxies.While images of a stellar population older than 1 Gyr are typically smooth,images of a stellar population younger than 1 Gyr reveal significant clumps,including multiple clumps which appear at different locations and even different ages.To study the evolutionary connections of these five galaxies to other dwarf galaxies before their star-forming clumps appear,we construct the images of the stellar populations older than three age nodes,and define them to be the images of the"host"galaxies.We find that the properties such as the central surface brightness and the effective radii of the hosts of the five galaxies are in between those of dwarf ellipticals(dEs)and dwarf irregulars(dIrrs),with two clearly more similar to dEs and one more similar to dIrrs.Among the five galaxies,8257-3704 is particularly interesting,as it shows a previous starburst event that is not quite visible from its gri image,but only visible from images of the stellar population at a few hundred million years.The star-forming clump associated with this event may have appeared at around 600 Myr ago and disappeared at around 40 Myr ago.

Relationship of Retinal Nerve Fiber Layer Thickness and Retinal Vessel Calibers with Cognitive Impairment in the Asymptomatic Polyvascular Abnormalities Population

Objective Cognitive impairment(CI)in older individuals has a high morbidity rate worldwide,with poor diagnostic methods and susceptible population identification.This study aimed to investigate the relationship between different retinal metrics and CI in a particular population,emphasizing polyvascular status.Methods We collected information from the Asymptomatic Polyvascular Abnormalities Community Study on retinal vessel calibers,retinal nerve fiber layer(RNFL)thickness,and cognitive function of 3,785participants,aged 40 years or older.Logistic regression was used to analyze the relationship between retinal metrics and cognitive function.Subgroups stratified by different vascular statuses were also analyzed.Results RNFL thickness was significantly thinner in the CI group(odds ratio:0.973,95%confidence interval:0.953–0.994).In the subgroup analysis,the difference still existed in the non-intracranial arterial stenosis,non-extracranial carotid arterial stenosis,and peripheral arterial disease subgroups(P<0.05).Conclusion A thin RNFL is associated with CI,especially in people with non-large vessel stenosis.The underlying small vessel change in RNFL and CI should be investigated in the future.

Mg/Fe site-specific dual-doping to boost the performance of cobalt-free nickle-rich layered oxide cathode for high-energy lithium-ion batteries

Layer-type LiNi0.9Mn0.1O2is promising to be the primary cathode material for lithium-ion batteries(LIBs)due to its excellent electrochemical performance.Unfortunately,the cathode with high nickel content suffers from severely detrimental structural transformation that causes rapid capacity attenuation.Herein,site-specific dual-doping with Fe and Mg ions is proposed to enhance the structural stability of LiNi0.9Mn0.1O2.The Fe3+dopants are inserted into transition metal sites(3b)and can favorably provide additional redox potential to compensate for charge and enhance the reversibility of anionic redox.The Mg ions are doped into the Li sites(3a)and serve as O_(2)^(-)-Mg^(2+)-O_(2)^(-)pillar to reinforce the electrostatic cohesion between the two adjacent transition-metal layers,which further suppress the cracking and the generation of harmful phase transitions,ultimately improving the cyclability.The theoretical calculations,including Bader charge and crystal orbital Hamilton populations(COHP)analyses,confirm that the doped Fe and Mg can form stable bonds with oxygen and the electrostatic repulsion of O_(2)^(-)-O_(2)^(-)can be effectively suppressed,which effectively mitigates oxygen anion loss at the high delithiation state.This dual-site doping strategy offers new avenues for understanding and regulating the crystalline oxygen redox and demonstrates significant potential for designing high-performance cobalt-free nickel-rich cathodes.

Mechanical responses of underground carriageway structures due to construction of metro tunnels beneath the existing structure:A case study

To understand the mechanical response pattern of the existing structure and ground due to the construction of metro tunnels underneath,the finite difference method is adopted to study the torsional deformation and stress variation of the existing structure and the effect of underground carriageway structures on the surface subsidence.The curves of the maximum differential subsidence,torsion angle,and distortion of the cross-section of the existing structure show two peaks in succession during traversing of two metro tunnels beneath it.The torsion angle of the existing structure changes when the two tunnels traverse beneath it in opposite directions.The first traversing of the shield tunnel mainly induces the magnitude variation in torsional deformation of the existing structure,but the second traversing of the subsurface tunnel may cause a dynamic change in the magnitude and form of torsional deformation in the existing structure.The shielding effect can reduce the surface subsidence caused by metro tunnel excavation to a certain extent,and the development trend of subsidence becomes slower as the excavation continues.

Laccase/caffeic acid-catalyzed crosslinking coupled with galactomannan alters the conformational structure of ovalbumin and alleviates Th2-mediated allergic asthma

Ovalbumin(OVA)is the major allergenic protein that can induce T helper 2(Th2)-allergic reactions,for which current treatment options are inadequate.In this study,we developed a polymerized hypoallergenic OVA product via laccase/caffeic acid(Lac/CA)-catalyzed crosslinking in conjunction with galactomannan(Man).The formation of high molecular weight crosslinked polymers and the Ig G-binding were analyzed by sodium dodecyl sulfate-polyacrylamide gel electrophoresis(SDS-PAGE)and Western blotting.The study indicated that Lac/CA-catalyzed crosslinking plus Man conjugation substantially altered secondary and tertiary structures of OVA along with the variation in surface hydrophobicity.Gastrointestinal digestion stability assay indicated that crosslinked OVA exhibited less resistance in simulated gastric fluid(SGF)and simulated intestinal fluid(SIF).Mouse model study indicated that Lac-Man/OVA ameliorated eosinophilic airway inflammatory response and efficiently downregulated the expression of Th2-related cytokines(interleukin(IL)-4,IL-5,and IL-13),and upregulated IFN-γand IL-10 expression.Stimulation of bone marrow-derived dendritic cells with Lac-Man/OVA suppressed the expression of phenotypic maturation markers(CD80 and CD86)and MHC class II molecules,and suppressed the expression levels of proinflammatory cytokines.The knowledge obtained in the present study offers an effective way to acquire a hypoallergenic OVA product that can have a therapeutic effect in alleviating OVA-induced allergic asthma.

Frequency-dependent viscoelasticity effects on the wave attenuation performance of multi-layered periodic foundations

In this paper,layered periodic foundations(LPFs)are numerically examined for their responses to longitudinal and transverse modes in the time and frequency domains.Three different unit-cells,i.e.,2-layer,4-layer,and 6-layer unit-cells,comprising concrete/rubber,concrete/rubber/steel/rubber,and concrete/rubber/steel/rubber/lead/rubber materials,respectively,are taken into account.Also,the viscoelasticity behavior of the rubber is modeled with two factors,i.e.,a frequency-independent(FI)loss factor and a linear frequency-dependent(FD)loss factor.Following the extraction of the complex dispersion curves and the identification of the band gaps(BGs),the simulations of wave transmission in the time and frequency domains are performed using the COMSOL software.Subsequent parametric studies evaluate the effects of the rubber viscoelasticity models on the dispersion curves and the wave transmission for the longitudinal and transverse modes.The results show that considering the rubber viscoelasticity enhances the wave attenuation performance.Moreover,the transverse-mode damping is more sensitive to the viscoelasticity model than its longitudinal counterpart.The 6-layer unit-cell LPF exhibits the lowest BG,ranging from 4.8 Hz to 6.5 Hz.

Coated boron layers by boronization and a real-time boron coating using an impurity powder dropper in the LHD

In the Large Helical Device(LHD),diborane(B2H6)is used as a standard boron source for boronization,which is assisted by helium glow discharges.In 2019,a new Impurity Powder Dropper(IPD)system was installed and is under evaluation as a real-time wall conditioning technique.In the LHD,which is a large-sized heliotron device,an additional helium(He)glow discharge cleaning(GDC)after boronization was operated for a reduction in hydrogen recycling from the coated boron layers.This operational time of 3 h was determined by spectroscopic data during glow discharges.A flat hydrogen profile is obtained on the top surface of the coated boron on the specimen exposed to boronization.The results suggest a reduction in hydrogen at the top surface by He-GDC.Trapped oxygen in coated boron was obtained by boronization,and the coated boron,which has boron-oxide,on the first wall by B-IPD was also shown.Considering the difference in coating areas between B2H6 boronization and B-IPD operation,it would be most effective to use the IPD and B2H6 boronization coating together for optimized wall conditioning.

A modified smoothed particle hydrodynamics method considering residual stress for simulating failure and its application in layered rock mass

Residual strength is an indispensable factor in evaluating rock fracture,yet the current Smoothed Particle Hydrodynamics(SPH)framework rarely considers its influence when simulating fracture.An improved cracking strategy considering residual stress in the base bond SPH method was proposed to simulate failures in layered rocks and slopes and verified by experimental results and other simulation methods(i.e.,the discrete element method).Modified Mohr–Coulomb failure criterion was applied to distinguish the mixed failure of tensile and shear.Bond fracture markψwas introduced to improve the kernel function after tensile damage,and the calculation of residual stress after the damage was derived after shear damage.Numerical simulations were carried out to evaluate its performance under different stress and scale conditions and to verify its effectiveness in realistically reproducing crack initiation and propagation and coalescence,even fracture and separation.The results indicate that the improved cracking strategy precisely captures the fracture and failure pattern in layered rocks and rock slopes.The residual stress of brittle tock is correctly captured by the improved SPH method.The improved SPH method that considers residual strength shows an approximately 13%improvement in accuracy for the safety factor of anti-dip layered slopes compared to the method that does not consider residual strength,as validated against analytical solutions.We infer that the improved SPH method is effective and shows promise for applications to continuous and discontinuous rock masses.