Temperature dependence of tunnel magnetoresistance (TMR) ratio, resistance, and coercivity from 4.2 K to room temperature (RT), applied de bias voltage dependence of the TMR ratio and resistances at 4.2 K and RT, tunn...Temperature dependence of tunnel magnetoresistance (TMR) ratio, resistance, and coercivity from 4.2 K to room temperature (RT), applied de bias voltage dependence of the TMR ratio and resistances at 4.2 K and RT, tunnel current I and dynamic conductance dI/dV as functions of the de bias voltage at 4.2 K, and inelastic electron tunneling (IET) spectroscopy, d(2)I/dV(2) versus V, at 4.2 K for a tunnel junction of Ta(5 nm)/Ni79Fe21(25 nm)/Ir22Mn78(12 nm)/Co75Fe25(4 nm)/Al(0.8 nm)-oxide/Co75Fe25(4 nm)/Ni79Fe21(25 nm)/Ta(5 nm) were systematically investigated. High TMR ratio of 59.2% at 4.2 K and 41.3% at RT were observed for this junction after annealing at 275 degreesC for an hour. The temperature dependence of TMR ratio and resistances from 4.2 to 300 K at 1.0 mV bias and the de bias voltage dependence of TMR ratio at 4.2 K from 0 to 80 mV can be evaluated by a comparison of self-consistent calculations with the experimental data based on the magnon-assisted inelastic excitation model and theory. An anisotropic wavelength cutoff energy of spin-wave spectrum in magnetic tunnel junctions (MTJs) was suggested, which is necessary for self-consistent calculations, based on a series of IET spectra observed in the MTJs.展开更多
In this study,we investigate the ef-ficacy of a hybrid parallel algo-rithm aiming at enhancing the speed of evaluation of two-electron repulsion integrals(ERI)and Fock matrix generation on the Hygon C86/DCU(deep compu...In this study,we investigate the ef-ficacy of a hybrid parallel algo-rithm aiming at enhancing the speed of evaluation of two-electron repulsion integrals(ERI)and Fock matrix generation on the Hygon C86/DCU(deep computing unit)heterogeneous computing platform.Multiple hybrid parallel schemes are assessed using a range of model systems,including those with up to 1200 atoms and 10000 basis func-tions.The findings of our research reveal that,during Hartree-Fock(HF)calculations,a single DCU ex-hibits 33.6 speedups over 32 C86 CPU cores.Compared with the efficiency of Wuhan Electronic Structure Package on Intel X86 and NVIDIA A100 computing platform,the Hygon platform exhibits good cost-effective-ness,showing great potential in quantum chemistry calculation and other high-performance scientific computations.展开更多
β-decay half-lives of some magic and semi-magic nuclei have been studied in a fully self-consistent Skyrme Hartree-Fock(HF) plus charge-exchange random phase approximation(RPA).The self-consistency is addressed,in th...β-decay half-lives of some magic and semi-magic nuclei have been studied in a fully self-consistent Skyrme Hartree-Fock(HF) plus charge-exchange random phase approximation(RPA).The self-consistency is addressed,in that the same Skyrme energy density functional is adopted in the calculation of ground states and Gamow-Teller excited states.First,the impact of J2 terms on the β-decay half-lives is investigated by using the SGII interaction,revealing a large influence.Subsequently,numerical calculations are performed for the selected nuclei with Skyrme energy density functionals SGII,LNS,SKX,and SAMi.Finally,comparisons to available experimental data and predictions of different theoretical models are discussed.展开更多
As cataract surgery progresses from “restoration of sight” to “refractive correction”, precise prediction of intraocular lens (IOL) power is critical for enhancing postoperative visual quality in patients. IOL pow...As cataract surgery progresses from “restoration of sight” to “refractive correction”, precise prediction of intraocular lens (IOL) power is critical for enhancing postoperative visual quality in patients. IOL power calculation methods have evolved and innovated throughout time, from early theoretical and regression formulas to nonlinear formulas for estimating effective lens position (ELP), multivariable formulas, and innovative formulas that use optical principles and AI-based online formulas. This paper thoroughly discusses the development and iteration of traditional IOL calculation formulas, the emergence of new IOL calculation formulas, and the selection of IOL calculation formulas for different patients in the era of refractive cataract surgery, serving as a reference for “personalized” IOL implantation in clinical practice.展开更多
An integrated calculated approach based on weakly coupled finite element(FEM)-viscoplastic selfconsistent(VPSC)model was established to simulate the texture evolution during the variable strain path extrusion process ...An integrated calculated approach based on weakly coupled finite element(FEM)-viscoplastic selfconsistent(VPSC)model was established to simulate the texture evolution during the variable strain path extrusion process of magnesium alloys.The spiral die extrusion(SDE)process with additional circumferential shear deformation was applied to investigate the effect of path control on texture adjustment and verify the accuracy of the model.The results indicated that the additional spiral shear resulting from the overall inclined flow path effectively reduced the intensity of the{0002}//ED fiber texture by suppressing basal slip activation in the core area,while the local shear deformation along the spiral equal channel strain path led to the formation of an inclined{0002}//ND plane texture on the side.Using the modified Hall-Petch relationship,the correlation between texture and yield strength was quantified.Specifically,the weakening of the texture effectively suppressed{10-12}tensile twinning,which compensated for the deficiency of compressive yield strength without significantly sacrificing tensile yield strength,and thus improved the tension-compression asymmetry.Furthermore,the strongly inclined{0002}//ND plane texture inhibited the widespread activation of basal slip during tensile yielding,thereby enhancing the yield strength.展开更多
The effects of pressure on the structural stability,elasticity,electronic properties,and thermodynamic properties of Al,Al_(3)Cu,Al_(2)Cu,Al_(4)Cu_(9),AlCu_(3),and Cu were investigated using first-principles calculati...The effects of pressure on the structural stability,elasticity,electronic properties,and thermodynamic properties of Al,Al_(3)Cu,Al_(2)Cu,Al_(4)Cu_(9),AlCu_(3),and Cu were investigated using first-principles calculations.The experimental results indicate that the calculated equilibrium lattice constant,elastic constant,and elastic modulus agree with both theoretical and experimental data at 0 GPa.The Young's modulus,bulk modulus,and shear modulus increase with increasing pressure.The influence of pressure on mechanical properties is explained from a chemical bond perspective.By employing the quasi-harmonic approximation model of phonon calculation,the temperature and pressure dependence of thermodynamic parameters in the range of 0 to 800 K and 0 to 100 GPa are determined.The findings demonstrate that the thermal capacity and coefficient of thermal expansion increase with increasing temperature and decrease with increasing pressure.This study provides fundamental data and support for experimental investigations and further theoretical research on the properties of aluminum-copper intermetallic compounds.展开更多
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.展开更多
Steel-concrete composite beams,due to their superior mechanical properties,are widely utilized in engineering structures.This study systematically investigates the calculation methods for internal forces and load-bear...Steel-concrete composite beams,due to their superior mechanical properties,are widely utilized in engineering structures.This study systematically investigates the calculation methods for internal forces and load-bearing capacity of composite beams based on elastic theory,with a focus on the transformed section method and its application under varying neutral axis positions.By deriving the geometric characteristics of the transformed section and incorporating a reduction factor accounting for slip effects,a computational model for sectional stress and ultimate load-bearing capacity is established.The results demonstrate that the slip effect significantly influences the flexural load-bearing capacity of composite beams.The proposed reduction factor,which considers the influence of the steel beam’s top flange thickness,offers higher accuracy compared to traditional methods.These findings provide a theoretical foundation for the design and analysis of composite beams,with significant practical engineering value.展开更多
Understanding the adsorption behavior of heavy metals and metalloids on clay minerals is essential for remediating heavy metal-contaminated soils.The adsorption of heavy metals and metalloids on illite(001)and sodium ...Understanding the adsorption behavior of heavy metals and metalloids on clay minerals is essential for remediating heavy metal-contaminated soils.The adsorption of heavy metals and metalloids on illite(001)and sodium montmorillonite(Na-MMT)(001)surfaces was investigated using first-principles calculations in this study,especially As atom and H_(3)AsO_(3) molecule.The adsorption energies of the As atom were−1.94 eV on the illite(001)and−0.56 eV on the Na-MMT(001),whereas,the adsorption energies of the H_(3)AsO_(3) molecule were−1.40 eV on illite(001)and−1.01 eV on Na-MMT(001).The above results indicate that the adsorption was more energetically favorable on illite(001).Additionally,compared to Na-MMT(001),there were more significant interactions between the atoms/molecules on the illite(001).After As atom and H_(3)AsO_(3) molecule adsorption,the electrons were transferred from mineral surface atoms to the adsorbates on both illite(001)and Na-MMT(001)surfaces.Moreover,the adsorption of As atom on illite(001)and Na-MMT(001)surfaces were more energy favorable compared to Hg,Cd,and Cr atoms.Overall,this work provides new insights into the adsorption behavior of As atoms and As molecules on illite and Na-MMT.The results indicate that illite rich soils are more prone to contamination by arsenic compared to soils primarily composed of Na-MMT minerals.展开更多
AIM:To investigate the effect of pharmacological pupil alterations on intraocular lens(IOL)power calculations.METHODS:A systematic review and Meta-analysis of studies published before December 2023 in the PubMed,Embas...AIM:To investigate the effect of pharmacological pupil alterations on intraocular lens(IOL)power calculations.METHODS:A systematic review and Meta-analysis of studies published before December 2023 in the PubMed,Embase,and Cochrane library databases on the accuracy of pharmacological pupil changes on IOL power calculation was performed.The primary outcome was the results of IOL power calculations before and after the use of medications.Subgroup analyses were performed based on participants’basic characteristics,such as age,axial length(AL),and whether miosis or mydriasis were used as classification criteria for further analyses.Each eligible study was evaluated for potential risk of bias by the AHRQ assessment scale.The study was registered on PROSPERO(CRD 42024497535).RESULTS:A total of 3062 eyes from 21 studies were eligible.There was no significant difference in the IOL power calculation before and after pharmacological pupil changes using any of the Hoffer Q(WMD=0.055,95%CI=-0.046–0.156;P=0.29),SRK/T(WMD=0.003,95%CI=-0.073–0.080;P=0.93),Haigis(WMD=-0.030,95%CI=-0.176–0.116;P=0.69),Holladay 2(WMD=-0.042,95%CI=-0.366–0.282;P=0.80),and Barrett Universal Ⅱ(WMD=0.033,95%CI=-0.061–0.127;P=0.49)formulas.On the measurement of parameters related to IOL power calculation,for either miosis or mydriasis AL(P=0.98 and 0.29,respectively),lens thickness(P=0.96 and 0.13,respectively),and mean keratometry(P=0.90 and 0.86,respectively)did not present significant differences,while anterior chamber depth(P=0.07 and<0.01,respectively)and white-to-white distance(P=0.01 and 0.04,respectively)changed significantly between the two measurements prior and posterior.At the same time,despite there being some participants with the difference between the before and after calculations greater than 0.5 diopter,there was no significant difference in the incidence rate between these formulas.CONCLUSION:There is no significant effect of pharmacological pupil changes on the IOL power calculation.It will considerably reduce the visit time burden for patients who require cataract surgery.展开更多
The longitudinal cracks distributed along the dam axis in the corridor of a dam may have potential safety hazards.According to the detection results of crack depth and width and the analysis of monitoring data,a three...The longitudinal cracks distributed along the dam axis in the corridor of a dam may have potential safety hazards.According to the detection results of crack depth and width and the analysis of monitoring data,a three-dimensional finite element model is established for numerical simulation calculation and the influence of cracks on the safety of dam structure is analyzed from different aspects such as deformation,stress value,and distribution range.The calculation results show that the maximum principal tensile stress value and the location of the dam body are basically independent of the change of crack depth(within 1.0 m).Regarding local stress around the corridor,the high upstream water level causes cracks to deepen,resulting in an increase in the maximum tensile stress near the crack tip and an expansion of the tensile stress region.展开更多
The micro-area characterization experiments like scanning Kelvin probe force microscope(SKPFM)and Kernel average misorientation have the defects of complex sample preparation and occasional errors in test results,whic...The micro-area characterization experiments like scanning Kelvin probe force microscope(SKPFM)and Kernel average misorientation have the defects of complex sample preparation and occasional errors in test results,which makes it impossible to accurately and quickly analyze the pitting behavior induced by inclusions in some cases,prompting attempts to turn to simulation calculation research.The method of calculating band structure and work function can be used to replace current-sensing atomic force microscopy and SKPFM to detect the potential and conductivity of the sample.The band structure results show that Al_(2)O_(3) inclusion is an insulator and non-conductive,and it will not form galvanic corrosion with the matrix.Al_(2)O_(3) inclusion does not dissolve because its work function is higher than that of the matrix.Moreover,the stress concentration of the matrix around the inclusion can be characterized by first-principles calculation coupled with finite element simulation.The results show that the stress concentration degree of the matrix around Al_(2)O_(3) inclusion is serious,and the galvanic corrosion is formed between the high and the low stress concentration areas,which can be used to explain the reason of the pitting induced by Al_(2)O_(3) inclusions.展开更多
In this study,6061 aluminum alloy and galvanized steel fusion-brazed lap joints were obtained using a laser-arc hybrid heat source,and the effects of laser power variation on the microstructure,mechanical properties,a...In this study,6061 aluminum alloy and galvanized steel fusion-brazed lap joints were obtained using a laser-arc hybrid heat source,and the effects of laser power variation on the microstructure,mechanical properties,and fracture mechanism of the joints were ana-lyzed.The results showed that the tensile shear load initially increased with rising laser power,followed by a decrease.At a laser power of 240 W,the maximum tensile shear load was 2479.8 N/cm and the weak section of joint was in the Al-Fe reaction layer con-sisting of Fe(Al,Si)_(3),Fe_(2)(Al,Si)_(5),and Fe(Al,Si)intermetallic compounds(IMCs).Computational results showed that the inherently high brittleness and hardness of Fe(Al,Si)_(3) and the high mismatch rates of Fe(Al,Si)_(3)/Al interfaces were the key factor leading to the failure of the joints at lower heat input.展开更多
Neutron-rich boron,carbon,and nitrogen isotopes have garnered extensive experimental and theoretical interest.In the present work,we conducted a comprehensive study of these nuclei by utilizing ab initio valence-space...Neutron-rich boron,carbon,and nitrogen isotopes have garnered extensive experimental and theoretical interest.In the present work,we conducted a comprehensive study of these nuclei by utilizing ab initio valence-space in-medium similarity renormalization group calculations with chiral nucleon-nucleon and three-nucleon interactions.First,we systematically calculated the spectra of nuclei.Our results align well with the available experimental data,which are comparable to phenomenological shell model calculations.Subsequently,the evolution of the N=14 and N=16 shell gaps is discussed based on the calculated spectra and the effective single-particle energies.Our calculations suggest that the N=14 neutron subshell is present in the oxygen isotopes but disappears in the boron,carbon,and nitrogen isotopic chains.Moreover,the N=16 subshell is present in all isotopes but gradually decreases from^(24)O to^(21)B.These results provide valuable information for future studies.展开更多
Two pairs of novel 6/6/6/9 tetracyclic merosesquiterpenoid enantiomers,dauroxonanols A(1)and B(2),possessing an unprecedented 9,15-dioxatetracyclo[8.5.3.0^(4.17).0^(14.18)]octadecane core skeleton,were isolated from R...Two pairs of novel 6/6/6/9 tetracyclic merosesquiterpenoid enantiomers,dauroxonanols A(1)and B(2),possessing an unprecedented 9,15-dioxatetracyclo[8.5.3.0^(4.17).0^(14.18)]octadecane core skeleton,were isolated from Rhododendron dauricum.The nuclear magnetic resonance(NMR)spectra of 1 and 2 showed very broad resonances,and^(13)C NMR spectrum of 1 exhibited only 13 instead of 22 carbon resonances.These broadening or missing NMR resonances led to a great challenge to elucidate their structures using NMR data analysis.Their structures and absolute configurations of 1 and 2 were finally determined by single crystal X-ray diffraction analysis,chiral separation,and electronic circular dichroism(ECD)calculations.Plausible biosynthetic pathways for 1 and 2 are proposed.Conformational analysis,density functional theory(DFT)calculations,and dynamic NMR assigned the coalescent NMR phenomena of 1 and 2 to the conformational changes of the flexible oxonane ring.Dauroxonanols A(1)and B(2)showed potentα-glucosidase inhibitory activities,2-8 times potent than acarbose,an antidiabetic drug targetingα-glucosidase in clinic.展开更多
Oxygen release from Ni-rich cathode is one of the major structural degradations resulting in rapid performance fading in lithium-ion batteries(LIBs).The energy gap between the transition metals(TM)-d band and the O-p ...Oxygen release from Ni-rich cathode is one of the major structural degradations resulting in rapid performance fading in lithium-ion batteries(LIBs).The energy gap between the transition metals(TM)-d band and the O-p band serves as an effective evaluation metric in characterizing the potential for oxygen release.Given that the primary oxidation factors of NCM811 materials vary at different states of charge(SOC),this study employs high-throughput density functional theory(DFT)calculations combined with machine learning(ML)to systematically investigate the regulation mechanism of heteroatoms on the energy gap between the TM-d band(TM=Ni,Co)and O-p band at various SOC levels.Highthroughput DFT calculations were used to study doping thermodynamic stability and complete the database.The results indicate that dopant atoms remain at their original sites even at 50%SOC.Correlation analysis reveals that at 0 SOC,the dopant reduces Ni-O bonding interactions by forming its own bonds with oxygen,thereby preventing lattice oxygen escape and weakening the oxygen binding of the system during Ni redox.At 50%SOC,the dopant and Co atoms synergistically strengthen their bonding interactions with oxygen,thereby maintaining structural stability and inhibiting lattice oxygen escape.Based on R^(2)and root-mean-square error(RMSE),the gradient boosting regression(GBR)algorithm is identified as optimal for predicting the energy gaps between the Ni-d band and O-p band,as well as between the Co-d band and O-p band.Feature importance analysis demonstrates that the magnetic moment(Dma)of the doped atom significantly contributes to the prediction of ΔNi-O and ΔCo-O.In this study,the energy gap regulation mechanisms of Ni-d/O-p and Co-d/O-p are systematically investigated using non-empirical first principle calculations combined with data-driven machine learning,aiming to provide insights into the electrochemical stability of NCM811 and related materials.展开更多
By comprehensively considering the influences of temperature and pressure on fluid density in high temperature and high pressure(HTHP)wells in deepwater fractured formations and the effects of formation fracture defor...By comprehensively considering the influences of temperature and pressure on fluid density in high temperature and high pressure(HTHP)wells in deepwater fractured formations and the effects of formation fracture deformation on well shut-in afterflow,this study couples the shut-in temperature field model,fracture deformation model,and gas flow model to establish a wellbore pressure calculation model incorporating thermo-hydro-mechanical coupling effects.The research analyzes the governing patterns of geothermal gradient,bottomhole pressure difference,drilling fluid pit gain,and kick index on casing head pressure,and establishes a shut-in pressure determination chart for HPHT wells based on coupled model calculation results.The study results show:geothermal gradient,bottomhole pressure difference,and drilling fluid pit gain exhibit positive correlations with casing head pressure;higher kick indices accelerate pressure rising rates while maintaining a constant maximum casing pressure;validation against field case data demonstrates over 95%accuracy in predicting wellbore pressure recovery after shut-in,with the pressure determination chart achieving 97.2%accuracy in target casing head pressure prediction and 98.3%accuracy in target shut-in time.This method enables accurate acquisition of formation pressure after HPHT well shut-in,providing reliable technical support for subsequent well control measures and ensuring safe and efficient development of deepwater and deep hydrocarbon reservoirs.展开更多
This study investigates the application of Gurney and flight of fragment equations,typically used to predict metal fragment velocities,in modeling the water jet behavior.Three shotgun cartridge sizes were used as the ...This study investigates the application of Gurney and flight of fragment equations,typically used to predict metal fragment velocities,in modeling the water jet behavior.Three shotgun cartridge sizes were used as the energy source:2.59 g,5.83 g,and 7.13 g.Two configurations were tested:standard(full-barrel water load)and"negative 8"(partial water load).High-speed footage captured water column velocities,and Gurney models,including infinitely tamped and open-faced configurations,combined with the flight of fragment model were used to assess prediction accuracy.Results showed charge strength significantly affects water column velocity,with higher strengths yielding greater stability and velocity retention over distance.The infinitely tamped Gurney model closely predicted experimental velocities,deviating by as little as 1.4%for standard charges and 2.8% for negative 8 charges.Additionally,interesting dynamics such as a 1-2°rise in jet height and the rear overtaking the front was observed.These findings have significant implications for optimizing PAN disruptors and enhancing performance in high-velocity fluid applications and explosive breaching systems.展开更多
AIM:To evaluate the accuracy of intraocular lens(IOL)power calculation formulas with/without preoperative aphakic anterior chamber depth(aph-ACD)in pediatric aphakia.METHODS:A total of 102 pediatric patients(150 eyes)...AIM:To evaluate the accuracy of intraocular lens(IOL)power calculation formulas with/without preoperative aphakic anterior chamber depth(aph-ACD)in pediatric aphakia.METHODS:A total of 102 pediatric patients(150 eyes)undergoing secondary IOL implantation were divided into two groups(in-the-bag or ciliary sulcus).Prediction error was calculated for 9 IOL power calculation formulas,including:1)not requiring ACD:Hoffer Q,Holladay 1,SRK/T;2)usable without or with entering ACD:Barrett Universal II(BUII),Emmetropia Verifying Optical(EVO)2.0,and Ladas Artificial Intelligence Super(Ladas AI);3)requiring ACD:Haigis,Kane,and Pearl-DGS.Mean prediction error(ME),mean absolute error(MAE),median absolute error(MedAE)and the percentage of eyes within±0.25,±0.50,±0.75,and±1.00 D were calculated.RESULTS:For the BUII,EVO 2.0,and Ladas AI,with aph-ACD demonstrated a higher MedAE compared to without aph-ACD(BUII:1.27 vs 1.13 D,EVO 2.0:1.26 vs 1.06 D,Ladas AI:1.30 vs 1.10 D;all P<0.05).Formulas requiring ACD(Haigis,Kane,and Pearl-DGS)exhibited larger MedAE than those not requiring aph-ACD(Hoffer Q,Holladay 1,and SRK/T;P<0.05).In the capsular group,the percentage of eyes within±1.00 D ranged from 44.83%to 74.14%,and it was 19.57%to 32.61%in the sulcus group.CONCLUSION:The introduction of aph-ACD does not improve the accuracy of IOL calculation for pediatric aphakia,regardless of in-the-bag or sulcus IOL secondary implantation.The relationship between aph-ACD and effective lens position in pediatric aphakia warrants further study.展开更多
Thermodynamic optimization of the AF-BeF_(2)(A=K,Rb,and Cs),KF-CsF,and RbF-CsF systems was performed within the framework of phase diagrams calculation.The model parameters were optimized based on experimental data an...Thermodynamic optimization of the AF-BeF_(2)(A=K,Rb,and Cs),KF-CsF,and RbF-CsF systems was performed within the framework of phase diagrams calculation.The model parameters were optimized based on experimental data and theoretically calculated values.The results show that the thermodynamically calculated values for the AF-BeF_(2)(A=K,Rb,and Cs),KF-CsF,and RbF-CsF systems agree well with the experimental data.Next,a set of reliable and self-consistent thermodynamic databases was built,and the liquidus projections and invariant points of the sub-ternary systems of the KF-RbF-CsF-BeF_(2)system were calculated.Furthermore,the melting temperature with the corresponding composition was predicted using the phase diagrams calculation technique,and the radial distribution functions,coordination numbers,angular distribution functions,and diffusion coefficients of the quaternary KF-RbF-CsF-BeF_(2)system were calculated using ab initio molecular dynamics.The results show that the quaternary KF-RbF-CsF-BeF_(2)system with the proportion 3.50-28.92-21.78-45.80 mol%or 1.80-35.42-52.40-10.38 mol%is one of the most promising candidate coolants for molten salt reactors in terms of thermodynamics and kinetics.This work provides direct guidelines for the screening and optimization of molten salts in the nuclear energy field.展开更多
文摘Temperature dependence of tunnel magnetoresistance (TMR) ratio, resistance, and coercivity from 4.2 K to room temperature (RT), applied de bias voltage dependence of the TMR ratio and resistances at 4.2 K and RT, tunnel current I and dynamic conductance dI/dV as functions of the de bias voltage at 4.2 K, and inelastic electron tunneling (IET) spectroscopy, d(2)I/dV(2) versus V, at 4.2 K for a tunnel junction of Ta(5 nm)/Ni79Fe21(25 nm)/Ir22Mn78(12 nm)/Co75Fe25(4 nm)/Al(0.8 nm)-oxide/Co75Fe25(4 nm)/Ni79Fe21(25 nm)/Ta(5 nm) were systematically investigated. High TMR ratio of 59.2% at 4.2 K and 41.3% at RT were observed for this junction after annealing at 275 degreesC for an hour. The temperature dependence of TMR ratio and resistances from 4.2 to 300 K at 1.0 mV bias and the de bias voltage dependence of TMR ratio at 4.2 K from 0 to 80 mV can be evaluated by a comparison of self-consistent calculations with the experimental data based on the magnon-assisted inelastic excitation model and theory. An anisotropic wavelength cutoff energy of spin-wave spectrum in magnetic tunnel junctions (MTJs) was suggested, which is necessary for self-consistent calculations, based on a series of IET spectra observed in the MTJs.
基金supported by the National Natural Science Foundation of China(No.22373112 to Ji Qi,No.22373111 and 21921004 to Minghui Yang)GH-fund A(No.202107011790)。
文摘In this study,we investigate the ef-ficacy of a hybrid parallel algo-rithm aiming at enhancing the speed of evaluation of two-electron repulsion integrals(ERI)and Fock matrix generation on the Hygon C86/DCU(deep computing unit)heterogeneous computing platform.Multiple hybrid parallel schemes are assessed using a range of model systems,including those with up to 1200 atoms and 10000 basis func-tions.The findings of our research reveal that,during Hartree-Fock(HF)calculations,a single DCU ex-hibits 33.6 speedups over 32 C86 CPU cores.Compared with the efficiency of Wuhan Electronic Structure Package on Intel X86 and NVIDIA A100 computing platform,the Hygon platform exhibits good cost-effective-ness,showing great potential in quantum chemistry calculation and other high-performance scientific computations.
基金Supported by the National Natural Science Foundation of China(11975096,11805280,11775014,11635003,11161130520)National Basic Research Program of China(2010CB832903)+1 种基金European Commission’s 7th Framework Programme(Fp7-PEOPLE-2010-IRSES)under Grant Agreement Project(269131)Continuous Basic Scientific Research Project(WDJC-2019-13)。
文摘β-decay half-lives of some magic and semi-magic nuclei have been studied in a fully self-consistent Skyrme Hartree-Fock(HF) plus charge-exchange random phase approximation(RPA).The self-consistency is addressed,in that the same Skyrme energy density functional is adopted in the calculation of ground states and Gamow-Teller excited states.First,the impact of J2 terms on the β-decay half-lives is investigated by using the SGII interaction,revealing a large influence.Subsequently,numerical calculations are performed for the selected nuclei with Skyrme energy density functionals SGII,LNS,SKX,and SAMi.Finally,comparisons to available experimental data and predictions of different theoretical models are discussed.
文摘As cataract surgery progresses from “restoration of sight” to “refractive correction”, precise prediction of intraocular lens (IOL) power is critical for enhancing postoperative visual quality in patients. IOL power calculation methods have evolved and innovated throughout time, from early theoretical and regression formulas to nonlinear formulas for estimating effective lens position (ELP), multivariable formulas, and innovative formulas that use optical principles and AI-based online formulas. This paper thoroughly discusses the development and iteration of traditional IOL calculation formulas, the emergence of new IOL calculation formulas, and the selection of IOL calculation formulas for different patients in the era of refractive cataract surgery, serving as a reference for “personalized” IOL implantation in clinical practice.
基金supported by the National Natural Science Foundation of China(Grant Nos.51975146,52205344)Shandong Province Natural Science Foundation(Grant No.ZR2020QE171)supported by the National Research Foundation of Korea(NRF)grant funded by the Korea government(MSIP)(Nos.NRF-2021R1A2C3006662,NRF-2022R1A5A1030054).
文摘An integrated calculated approach based on weakly coupled finite element(FEM)-viscoplastic selfconsistent(VPSC)model was established to simulate the texture evolution during the variable strain path extrusion process of magnesium alloys.The spiral die extrusion(SDE)process with additional circumferential shear deformation was applied to investigate the effect of path control on texture adjustment and verify the accuracy of the model.The results indicated that the additional spiral shear resulting from the overall inclined flow path effectively reduced the intensity of the{0002}//ED fiber texture by suppressing basal slip activation in the core area,while the local shear deformation along the spiral equal channel strain path led to the formation of an inclined{0002}//ND plane texture on the side.Using the modified Hall-Petch relationship,the correlation between texture and yield strength was quantified.Specifically,the weakening of the texture effectively suppressed{10-12}tensile twinning,which compensated for the deficiency of compressive yield strength without significantly sacrificing tensile yield strength,and thus improved the tension-compression asymmetry.Furthermore,the strongly inclined{0002}//ND plane texture inhibited the widespread activation of basal slip during tensile yielding,thereby enhancing the yield strength.
基金Funded by the National Key R&D Program of China(No.2021YFB3802300)the Foundation of National Key Laboratory of Shock Wave and Detonation Physics(No.JCKYS2022212004)the National Natural Science Foundation of China(No.52171045),and the Joint Fund(No.8091B022108)。
文摘The effects of pressure on the structural stability,elasticity,electronic properties,and thermodynamic properties of Al,Al_(3)Cu,Al_(2)Cu,Al_(4)Cu_(9),AlCu_(3),and Cu were investigated using first-principles calculations.The experimental results indicate that the calculated equilibrium lattice constant,elastic constant,and elastic modulus agree with both theoretical and experimental data at 0 GPa.The Young's modulus,bulk modulus,and shear modulus increase with increasing pressure.The influence of pressure on mechanical properties is explained from a chemical bond perspective.By employing the quasi-harmonic approximation model of phonon calculation,the temperature and pressure dependence of thermodynamic parameters in the range of 0 to 800 K and 0 to 100 GPa are determined.The findings demonstrate that the thermal capacity and coefficient of thermal expansion increase with increasing temperature and decrease with increasing pressure.This study provides fundamental data and support for experimental investigations and further theoretical research on the properties of aluminum-copper intermetallic compounds.
基金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.
文摘Steel-concrete composite beams,due to their superior mechanical properties,are widely utilized in engineering structures.This study systematically investigates the calculation methods for internal forces and load-bearing capacity of composite beams based on elastic theory,with a focus on the transformed section method and its application under varying neutral axis positions.By deriving the geometric characteristics of the transformed section and incorporating a reduction factor accounting for slip effects,a computational model for sectional stress and ultimate load-bearing capacity is established.The results demonstrate that the slip effect significantly influences the flexural load-bearing capacity of composite beams.The proposed reduction factor,which considers the influence of the steel beam’s top flange thickness,offers higher accuracy compared to traditional methods.These findings provide a theoretical foundation for the design and analysis of composite beams,with significant practical engineering value.
基金Project(22376221)supported by the National Natural Science Foundation of ChinaProject(2024JJ2074)supported by the Natural Science Foundation of Hunan Province,ChinaProject(2023QNRC001)supported by the Young Elite Scientists Sponsorship Program by CAST。
文摘Understanding the adsorption behavior of heavy metals and metalloids on clay minerals is essential for remediating heavy metal-contaminated soils.The adsorption of heavy metals and metalloids on illite(001)and sodium montmorillonite(Na-MMT)(001)surfaces was investigated using first-principles calculations in this study,especially As atom and H_(3)AsO_(3) molecule.The adsorption energies of the As atom were−1.94 eV on the illite(001)and−0.56 eV on the Na-MMT(001),whereas,the adsorption energies of the H_(3)AsO_(3) molecule were−1.40 eV on illite(001)and−1.01 eV on Na-MMT(001).The above results indicate that the adsorption was more energetically favorable on illite(001).Additionally,compared to Na-MMT(001),there were more significant interactions between the atoms/molecules on the illite(001).After As atom and H_(3)AsO_(3) molecule adsorption,the electrons were transferred from mineral surface atoms to the adsorbates on both illite(001)and Na-MMT(001)surfaces.Moreover,the adsorption of As atom on illite(001)and Na-MMT(001)surfaces were more energy favorable compared to Hg,Cd,and Cr atoms.Overall,this work provides new insights into the adsorption behavior of As atoms and As molecules on illite and Na-MMT.The results indicate that illite rich soils are more prone to contamination by arsenic compared to soils primarily composed of Na-MMT minerals.
基金Supported by Beijing Natural Science Foundation from Beijing Municipal Government(No.7202030).
文摘AIM:To investigate the effect of pharmacological pupil alterations on intraocular lens(IOL)power calculations.METHODS:A systematic review and Meta-analysis of studies published before December 2023 in the PubMed,Embase,and Cochrane library databases on the accuracy of pharmacological pupil changes on IOL power calculation was performed.The primary outcome was the results of IOL power calculations before and after the use of medications.Subgroup analyses were performed based on participants’basic characteristics,such as age,axial length(AL),and whether miosis or mydriasis were used as classification criteria for further analyses.Each eligible study was evaluated for potential risk of bias by the AHRQ assessment scale.The study was registered on PROSPERO(CRD 42024497535).RESULTS:A total of 3062 eyes from 21 studies were eligible.There was no significant difference in the IOL power calculation before and after pharmacological pupil changes using any of the Hoffer Q(WMD=0.055,95%CI=-0.046–0.156;P=0.29),SRK/T(WMD=0.003,95%CI=-0.073–0.080;P=0.93),Haigis(WMD=-0.030,95%CI=-0.176–0.116;P=0.69),Holladay 2(WMD=-0.042,95%CI=-0.366–0.282;P=0.80),and Barrett Universal Ⅱ(WMD=0.033,95%CI=-0.061–0.127;P=0.49)formulas.On the measurement of parameters related to IOL power calculation,for either miosis or mydriasis AL(P=0.98 and 0.29,respectively),lens thickness(P=0.96 and 0.13,respectively),and mean keratometry(P=0.90 and 0.86,respectively)did not present significant differences,while anterior chamber depth(P=0.07 and<0.01,respectively)and white-to-white distance(P=0.01 and 0.04,respectively)changed significantly between the two measurements prior and posterior.At the same time,despite there being some participants with the difference between the before and after calculations greater than 0.5 diopter,there was no significant difference in the incidence rate between these formulas.CONCLUSION:There is no significant effect of pharmacological pupil changes on the IOL power calculation.It will considerably reduce the visit time burden for patients who require cataract surgery.
基金Zhejiang Provincial Natural Science Foundation of China for Young Scholars(Project No.:LQ20A020009)National College Students’Innovation and Entrepreneurship Training Program(Project No.:202311842014X)。
文摘The longitudinal cracks distributed along the dam axis in the corridor of a dam may have potential safety hazards.According to the detection results of crack depth and width and the analysis of monitoring data,a three-dimensional finite element model is established for numerical simulation calculation and the influence of cracks on the safety of dam structure is analyzed from different aspects such as deformation,stress value,and distribution range.The calculation results show that the maximum principal tensile stress value and the location of the dam body are basically independent of the change of crack depth(within 1.0 m).Regarding local stress around the corridor,the high upstream water level causes cracks to deepen,resulting in an increase in the maximum tensile stress near the crack tip and an expansion of the tensile stress region.
基金supported by the National Natural Science Foundation of China(Nos.52364044 and 52204364)Central Guidance on Local Science and Technology Development Fund Projects of Inner Mongolia Autonomous Region(No.2022ZY0090)Basic Scientific Research Business Expenses of Colleges and Universities in Inner Mongolia Autonomous Region(Nos.2023QNJS011 and 0406082226).
文摘The micro-area characterization experiments like scanning Kelvin probe force microscope(SKPFM)and Kernel average misorientation have the defects of complex sample preparation and occasional errors in test results,which makes it impossible to accurately and quickly analyze the pitting behavior induced by inclusions in some cases,prompting attempts to turn to simulation calculation research.The method of calculating band structure and work function can be used to replace current-sensing atomic force microscopy and SKPFM to detect the potential and conductivity of the sample.The band structure results show that Al_(2)O_(3) inclusion is an insulator and non-conductive,and it will not form galvanic corrosion with the matrix.Al_(2)O_(3) inclusion does not dissolve because its work function is higher than that of the matrix.Moreover,the stress concentration of the matrix around the inclusion can be characterized by first-principles calculation coupled with finite element simulation.The results show that the stress concentration degree of the matrix around Al_(2)O_(3) inclusion is serious,and the galvanic corrosion is formed between the high and the low stress concentration areas,which can be used to explain the reason of the pitting induced by Al_(2)O_(3) inclusions.
基金supported by the National Key Research and Development Program of China(No.2022YFB4600900).
文摘In this study,6061 aluminum alloy and galvanized steel fusion-brazed lap joints were obtained using a laser-arc hybrid heat source,and the effects of laser power variation on the microstructure,mechanical properties,and fracture mechanism of the joints were ana-lyzed.The results showed that the tensile shear load initially increased with rising laser power,followed by a decrease.At a laser power of 240 W,the maximum tensile shear load was 2479.8 N/cm and the weak section of joint was in the Al-Fe reaction layer con-sisting of Fe(Al,Si)_(3),Fe_(2)(Al,Si)_(5),and Fe(Al,Si)intermetallic compounds(IMCs).Computational results showed that the inherently high brittleness and hardness of Fe(Al,Si)_(3) and the high mismatch rates of Fe(Al,Si)_(3)/Al interfaces were the key factor leading to the failure of the joints at lower heat input.
基金supported by the National Key R&D Program of China(Nos.2024YFE0109800 and 2024YFE0109802)the National Natural Science Foundation of China(Nos.12405141,12347106,12205340,and 12121005)+1 种基金the Gansu Natural Science Foundation(No.25JRRA467)the Key Research Program of the Chinese Academy of Sciences(No.XDPB15)。
文摘Neutron-rich boron,carbon,and nitrogen isotopes have garnered extensive experimental and theoretical interest.In the present work,we conducted a comprehensive study of these nuclei by utilizing ab initio valence-space in-medium similarity renormalization group calculations with chiral nucleon-nucleon and three-nucleon interactions.First,we systematically calculated the spectra of nuclei.Our results align well with the available experimental data,which are comparable to phenomenological shell model calculations.Subsequently,the evolution of the N=14 and N=16 shell gaps is discussed based on the calculated spectra and the effective single-particle energies.Our calculations suggest that the N=14 neutron subshell is present in the oxygen isotopes but disappears in the boron,carbon,and nitrogen isotopic chains.Moreover,the N=16 subshell is present in all isotopes but gradually decreases from^(24)O to^(21)B.These results provide valuable information for future studies.
基金supported by the National Natural Science Foundation of China(Nos.22207036,22277034,22477034,and 22107033)Interdisciplinary Research Program of Huazhong University of Science and Technology(No.2023JCYJ037)International Cooperation Project of Hubei Provincial Key R&D Plan(No.2023EHA040)。
文摘Two pairs of novel 6/6/6/9 tetracyclic merosesquiterpenoid enantiomers,dauroxonanols A(1)and B(2),possessing an unprecedented 9,15-dioxatetracyclo[8.5.3.0^(4.17).0^(14.18)]octadecane core skeleton,were isolated from Rhododendron dauricum.The nuclear magnetic resonance(NMR)spectra of 1 and 2 showed very broad resonances,and^(13)C NMR spectrum of 1 exhibited only 13 instead of 22 carbon resonances.These broadening or missing NMR resonances led to a great challenge to elucidate their structures using NMR data analysis.Their structures and absolute configurations of 1 and 2 were finally determined by single crystal X-ray diffraction analysis,chiral separation,and electronic circular dichroism(ECD)calculations.Plausible biosynthetic pathways for 1 and 2 are proposed.Conformational analysis,density functional theory(DFT)calculations,and dynamic NMR assigned the coalescent NMR phenomena of 1 and 2 to the conformational changes of the flexible oxonane ring.Dauroxonanols A(1)and B(2)showed potentα-glucosidase inhibitory activities,2-8 times potent than acarbose,an antidiabetic drug targetingα-glucosidase in clinic.
基金supported by the National Natural Science Foundation of China(Grant no.52463025,and 52062035)the Major Discipline Academic and Technical Leaders Training Program of Jiangxi Province(Grant no.20213BCJ22056)+2 种基金the Key R&D Program of Jiangxi Province(Grant no.20223BBE51028)the Jiangxi Province Key Laboratory of Lithium-ion Battery Materials and Application(2024SSY05202)the Jiangxi Province Graduate Innovation Special Fund Project(YC2023-B004)。
文摘Oxygen release from Ni-rich cathode is one of the major structural degradations resulting in rapid performance fading in lithium-ion batteries(LIBs).The energy gap between the transition metals(TM)-d band and the O-p band serves as an effective evaluation metric in characterizing the potential for oxygen release.Given that the primary oxidation factors of NCM811 materials vary at different states of charge(SOC),this study employs high-throughput density functional theory(DFT)calculations combined with machine learning(ML)to systematically investigate the regulation mechanism of heteroatoms on the energy gap between the TM-d band(TM=Ni,Co)and O-p band at various SOC levels.Highthroughput DFT calculations were used to study doping thermodynamic stability and complete the database.The results indicate that dopant atoms remain at their original sites even at 50%SOC.Correlation analysis reveals that at 0 SOC,the dopant reduces Ni-O bonding interactions by forming its own bonds with oxygen,thereby preventing lattice oxygen escape and weakening the oxygen binding of the system during Ni redox.At 50%SOC,the dopant and Co atoms synergistically strengthen their bonding interactions with oxygen,thereby maintaining structural stability and inhibiting lattice oxygen escape.Based on R^(2)and root-mean-square error(RMSE),the gradient boosting regression(GBR)algorithm is identified as optimal for predicting the energy gaps between the Ni-d band and O-p band,as well as between the Co-d band and O-p band.Feature importance analysis demonstrates that the magnetic moment(Dma)of the doped atom significantly contributes to the prediction of ΔNi-O and ΔCo-O.In this study,the energy gap regulation mechanisms of Ni-d/O-p and Co-d/O-p are systematically investigated using non-empirical first principle calculations combined with data-driven machine learning,aiming to provide insights into the electrochemical stability of NCM811 and related materials.
基金Supported by the Joint Fund Key Program of the National Natural Science Foundation of China(U21B2069)Key Research and Development Program of Shandong Province(2022CXGC020407)Basic Science Center Program of the National Natural Science Foundation of China(52288101)。
文摘By comprehensively considering the influences of temperature and pressure on fluid density in high temperature and high pressure(HTHP)wells in deepwater fractured formations and the effects of formation fracture deformation on well shut-in afterflow,this study couples the shut-in temperature field model,fracture deformation model,and gas flow model to establish a wellbore pressure calculation model incorporating thermo-hydro-mechanical coupling effects.The research analyzes the governing patterns of geothermal gradient,bottomhole pressure difference,drilling fluid pit gain,and kick index on casing head pressure,and establishes a shut-in pressure determination chart for HPHT wells based on coupled model calculation results.The study results show:geothermal gradient,bottomhole pressure difference,and drilling fluid pit gain exhibit positive correlations with casing head pressure;higher kick indices accelerate pressure rising rates while maintaining a constant maximum casing pressure;validation against field case data demonstrates over 95%accuracy in predicting wellbore pressure recovery after shut-in,with the pressure determination chart achieving 97.2%accuracy in target casing head pressure prediction and 98.3%accuracy in target shut-in time.This method enables accurate acquisition of formation pressure after HPHT well shut-in,providing reliable technical support for subsequent well control measures and ensuring safe and efficient development of deepwater and deep hydrocarbon reservoirs.
基金supported and funded internally through Dr. Catherine Johnson's research funds at Missouri S&T
文摘This study investigates the application of Gurney and flight of fragment equations,typically used to predict metal fragment velocities,in modeling the water jet behavior.Three shotgun cartridge sizes were used as the energy source:2.59 g,5.83 g,and 7.13 g.Two configurations were tested:standard(full-barrel water load)and"negative 8"(partial water load).High-speed footage captured water column velocities,and Gurney models,including infinitely tamped and open-faced configurations,combined with the flight of fragment model were used to assess prediction accuracy.Results showed charge strength significantly affects water column velocity,with higher strengths yielding greater stability and velocity retention over distance.The infinitely tamped Gurney model closely predicted experimental velocities,deviating by as little as 1.4%for standard charges and 2.8% for negative 8 charges.Additionally,interesting dynamics such as a 1-2°rise in jet height and the rear overtaking the front was observed.These findings have significant implications for optimizing PAN disruptors and enhancing performance in high-velocity fluid applications and explosive breaching systems.
基金Supported by the Joint Funding Project of Municipal Schools(Colleges)of Science and Technology Program of Guangzhou,China(No.2023A03J0188)the Construction Project of High-Level Hospitals in Guangdong Province(No.303020102)the Natural Science Fund of Guangdong Province(No.2023A1515011102).
文摘AIM:To evaluate the accuracy of intraocular lens(IOL)power calculation formulas with/without preoperative aphakic anterior chamber depth(aph-ACD)in pediatric aphakia.METHODS:A total of 102 pediatric patients(150 eyes)undergoing secondary IOL implantation were divided into two groups(in-the-bag or ciliary sulcus).Prediction error was calculated for 9 IOL power calculation formulas,including:1)not requiring ACD:Hoffer Q,Holladay 1,SRK/T;2)usable without or with entering ACD:Barrett Universal II(BUII),Emmetropia Verifying Optical(EVO)2.0,and Ladas Artificial Intelligence Super(Ladas AI);3)requiring ACD:Haigis,Kane,and Pearl-DGS.Mean prediction error(ME),mean absolute error(MAE),median absolute error(MedAE)and the percentage of eyes within±0.25,±0.50,±0.75,and±1.00 D were calculated.RESULTS:For the BUII,EVO 2.0,and Ladas AI,with aph-ACD demonstrated a higher MedAE compared to without aph-ACD(BUII:1.27 vs 1.13 D,EVO 2.0:1.26 vs 1.06 D,Ladas AI:1.30 vs 1.10 D;all P<0.05).Formulas requiring ACD(Haigis,Kane,and Pearl-DGS)exhibited larger MedAE than those not requiring aph-ACD(Hoffer Q,Holladay 1,and SRK/T;P<0.05).In the capsular group,the percentage of eyes within±1.00 D ranged from 44.83%to 74.14%,and it was 19.57%to 32.61%in the sulcus group.CONCLUSION:The introduction of aph-ACD does not improve the accuracy of IOL calculation for pediatric aphakia,regardless of in-the-bag or sulcus IOL secondary implantation.The relationship between aph-ACD and effective lens position in pediatric aphakia warrants further study.
基金supported by the National Natural Science Foundation of China(Nos.12205364 and 12375282)Guangdong Provincial Natural Science Foundation(Nos.2024A1515012570 and 2024A1515010885)the Fundamental Research funds for the Central Universities,Sun Yat sen University。
文摘Thermodynamic optimization of the AF-BeF_(2)(A=K,Rb,and Cs),KF-CsF,and RbF-CsF systems was performed within the framework of phase diagrams calculation.The model parameters were optimized based on experimental data and theoretically calculated values.The results show that the thermodynamically calculated values for the AF-BeF_(2)(A=K,Rb,and Cs),KF-CsF,and RbF-CsF systems agree well with the experimental data.Next,a set of reliable and self-consistent thermodynamic databases was built,and the liquidus projections and invariant points of the sub-ternary systems of the KF-RbF-CsF-BeF_(2)system were calculated.Furthermore,the melting temperature with the corresponding composition was predicted using the phase diagrams calculation technique,and the radial distribution functions,coordination numbers,angular distribution functions,and diffusion coefficients of the quaternary KF-RbF-CsF-BeF_(2)system were calculated using ab initio molecular dynamics.The results show that the quaternary KF-RbF-CsF-BeF_(2)system with the proportion 3.50-28.92-21.78-45.80 mol%or 1.80-35.42-52.40-10.38 mol%is one of the most promising candidate coolants for molten salt reactors in terms of thermodynamics and kinetics.This work provides direct guidelines for the screening and optimization of molten salts in the nuclear energy field.