Chemical short-range order(SRO),a phenomenon at the atomic scale resulting from inhomogeneities in the local chemical environment,is usually studied using machine learning force field-based molecular dynamics simulati...Chemical short-range order(SRO),a phenomenon at the atomic scale resulting from inhomogeneities in the local chemical environment,is usually studied using machine learning force field-based molecular dynamics simulations due to the limitations of experimental methods.To promote the reliable application of machine potentials in high-entropy alloy simulations,first,this work uses NEP models trained on two different datasets to predict the SRO coefficients of NbMoTaW.The results show that within the same machine learning framework,there are significant differences in the prediction of SRO coefficients for the Nb-Nb atomic pair.Subsequently,this work predicts the SRO coefficients of NbMoTaW using the NEP model and the SNAP model,both of which are trained on the same dataset.The results reveal significant discrepancies in SRO predictions for like-element pairs(e.g.,Nb-Nb and W-W)between the two potentials,despite the identical training data.The findings of this study indicate that discrepancies in the prediction results of SRO coefficients can arise from either the same machine learning framework trained on different datasets or different learning frameworks trained on the same dataset.This reflects possible incompleteness in the current training set's coverage of local chemical environments at the atomic scale.Future research should establish unified evaluation standards to assess the capability of training sets to accurately describe complex atomic-scale behaviors such as SRO.展开更多
We present multi-threading and SIMD optimizations on short-range potential calculation kernel in Molecular Dynamics.For the multi-threading optimization,we design a partition-and-two-steps(PTS)method to avoid write co...We present multi-threading and SIMD optimizations on short-range potential calculation kernel in Molecular Dynamics.For the multi-threading optimization,we design a partition-and-two-steps(PTS)method to avoid write conflicts caused by using Newton’s third law.Our method eliminates serialization bottle-neck without extra memory.We implement our PTS method using OpenMP.Afterwards,we discuss the influence of the cutoff if statement on the performance of vectorization in MD simulations.We propose a pre-searching neighbors method,which makes about 70%atoms meet the cutoff check,reducing a large amount of redundant calculation.The experiment results prove our PTS method is scalable and efficient.In double precision,our 256-bit SIMD implementation is about 3×faster than the scalar version.展开更多
Using the asymptotic iteration method, we obtain the S-wave solution for a short-range three-parameter central potential with 1/r singularity and with a non-orbital barrier. To the best of our knowledge, this is the f...Using the asymptotic iteration method, we obtain the S-wave solution for a short-range three-parameter central potential with 1/r singularity and with a non-orbital barrier. To the best of our knowledge, this is the first attempt at calculating the energy spectrum for this potential, which was introduced by H. Bahlouli and A. D. Alhaidari and for which they obtained the “potential parameter spectrum”. Our results are also independently verified using a direct method of diagonalizing the Hamiltonian matrix in the J-matrix basis.展开更多
This paper presents recurrence spectra of highly excited lithium atoms with M = 1 state in parallel electric and magnetic fields at a fixed scaled energy ε = -0.03. Short-ranged potentials including ionic core potent...This paper presents recurrence spectra of highly excited lithium atoms with M = 1 state in parallel electric and magnetic fields at a fixed scaled energy ε = -0.03. Short-ranged potentials including ionic core potential and centrifugal barrier are taken into account. Their effects on the states and photo-absorption spectrum are analysed in detail. This demonstrates that the geometric features of classical orbits are of special importance for modulations of the spectral pattern. Thus the weak polarization as well as the reduction of correlation of electrons induced by short-ranged potentials give rise to the recurrence spectra of lithium M = 1 atoms more compact than that of the M = 0 one, which is in good agreement with the experimental prediction.展开更多
Global optimization of Morse clusters with shortrange potential is a great challenge.Here,we apply our recently developed unbiased fuzzy global optimization method to systematically study Morse clusters with the poten...Global optimization of Morse clusters with shortrange potential is a great challenge.Here,we apply our recently developed unbiased fuzzy global optimization method to systematically study Morse clusters with the potential rangeρ=14 and the number of atoms N up to 400.All the putative global minima reported in the literature have been successfully reproduced with relatively high success ratios.Compared to the available results for N≤240 and several larger Morse clusters,new global minima(and local minima)with lower energies have been found out for N=164,175,188,193,194,197,239,246,260,318,and 389.Clusters with magic numbers are figured out through fitting the size-dependent global minimum energies.The cluster structures tend to be close-packed for short-range potential with large N.展开更多
Quantum energies which are used in applications are usually composed of repulsive and attractive terms. The objective of this study is to use an accurate and efficient fitting of the repulsive energy instead of using ...Quantum energies which are used in applications are usually composed of repulsive and attractive terms. The objective of this study is to use an accurate and efficient fitting of the repulsive energy instead of using standard parametrizations. The investigation is based on Density Functional Theory and Tight Binding simulations. Our objective is not only to capture the values of the repulsive terms but also to efficiently reproduce the elastic properties and the forces. The elasticity values determine the rigidity of a material when some traction or load is applied on it. The pair-potential is based on an exponential term corrected by B-spline terms. In order to accelerate the computations, one uses a hierarchical optimization for the B-splines on different levels. Carbon graphenes constitute the configurations used in the simulations. We report on some results to show the efficiency of the B-splines on different levels.展开更多
Quantum energies which are used in applications are usually composed of repulsive and attractive terms. The objective of this study is to use an accurate and efficient fitting of the repulsive energy instead of using ...Quantum energies which are used in applications are usually composed of repulsive and attractive terms. The objective of this study is to use an accurate and efficient fitting of the repulsive energy instead of using standard parametrizations. The investigation is based on Density Functional Theory and Tight Binding simulations. Our objective is not only to capture the values of the repulsive terms but also to efficiently reproduce the elastic properties and the forces. The elasticity values determine the rigidity of a material when some traction or load is applied on it. The pair-potential is based on an exponential term corrected by B-spline terms. In order to accelerate the computations, one uses a hierarchical optimization for the B-splines on different levels. Carbon graphenes constitute the configurations used in the simulations. We report on some results to show the efficiency of the B-splines on different levels.展开更多
Classical scattering trajectories are known to form a Lagrangian manifold in euclidean phase space,which allows the classification of local focal points for sufficiently small dimensions. For the case of a short-range...Classical scattering trajectories are known to form a Lagrangian manifold in euclidean phase space,which allows the classification of local focal points for sufficiently small dimensions. For the case of a short-range potential, we show that the natural description of focal points at infinity is a Lagrangian manifold in the cotangent bundle of the sphere and establish the relationship between focal points at infinity and the projection singularities of that manifold.展开更多
We study the confinement of a spinless charged particle to a spherical quantum dot under the influence of a linear electric field.The spherical quantum dot is described by a short-range potential given by the power-ex...We study the confinement of a spinless charged particle to a spherical quantum dot under the influence of a linear electric field.The spherical quantum dot is described by a short-range potential given by the power-exponential potential.Then,by analysing the region near the spherical quantum dot centre,we discuss two cases where the energy levels can be obtained for s-waves and how the linear electric field modifies the spectrum of energy of the spherical quantum dot.展开更多
High/medium entropy alloys(H/MEAs)have shown unique strengthening behavior and mechanical prop-erties because of the presence of massive local chemical orderings.Nevertheless,dynamic interactions between chemical shor...High/medium entropy alloys(H/MEAs)have shown unique strengthening behavior and mechanical prop-erties because of the presence of massive local chemical orderings.Nevertheless,dynamic interactions between chemical short-range orders(CSROs)and dislocations,and the underlying atomic strengthening mechanism remain elusive.In this work,we first developed a novel machine learning-embedded atom method(ML-EAM)potential of the CoNiV system,trained on a comprehensive first-principles dataset,which enables accurate and efficient modeling of CSRO formation and dislocation dynamics.Then,we in-vestigated the strengthening mechanisms of CSROs in CoNiV MEA through machine learning-augmented molecular dynamics(MD)simulations.Hybrid MD/Monte Carlo simulations reveal that CSRO domains possess an L1_(2)(NiCo)_(3) V structure,whose size increases with lowering annealing temperatures.These domains significantly enhance strength by impeding dislocation motion through complex energy path-ways,increasing depinning forces,and reducing mobility.Moreover,the MD simulations combined with theoretical analysis elucidate the competition between CSRO-assisted strengthening(via antiphase bound-ary formation)and solid solution weakening(via reduced atomic misfit volume).Phonon-drag effects are also amplified by CSROs,further resisting dislocation glide.Our results demonstrate that L1_(2)-CSROs strengthen CoNiV MEA primarily through antiphase boundary and phonon-drag contributions,providing new insights for designing high-performance multi-principal-element alloys via tailoring CSROs.展开更多
We consider a pair of Hamiltonians （H, H0） on L2（R^n）, where H0=p^2 -x^2 is a SchrSdinger operator with a repulsive potential, and H = H0＋V（x）. We show that, under suitable assumptions on the decay of the elect...We consider a pair of Hamiltonians （H, H0） on L2（R^n）, where H0=p^2 -x^2 is a SchrSdinger operator with a repulsive potential, and H = H0＋V（x）. We show that, under suitable assumptions on the decay of the electric potential, V is uniquely determined by the high energy limit of the scattering operator.展开更多
An accurate total energy calculation is essential in materials computation.To date,many tight-binding(TB)approaches based on parameterized hopping can produce electronic structures comparable to those obtained using f...An accurate total energy calculation is essential in materials computation.To date,many tight-binding(TB)approaches based on parameterized hopping can produce electronic structures comparable to those obtained using first-principles calculations.However,TB approaches still have limited applicability for determining material properties derived from the total energy.That is,the predictive power of the TB total energy is impaired by an inaccurate evaluation of the repulsive energy.The complexity associated with the parametrization of TB repulsive potentials is the weak link in this evaluation.In this study,we propose a new method for obtaining the pairwise TB repulsive potential for crystalline materials by employing the Chen-Möbius inversion theorem.We show that the TB-based phonon dispersions,calculated using the resulting repulsive potential,compare well with those obtained by first-principles calculations for various systems,including covalent and ionic bulk materials and twodimensional materials.The present approach only requires the first-principles total energy and TB electronic band energy as input and does not involve any parameters.This striking feature enables us to generate repulsive potentials programmatically.展开更多
基金Project supported by the Hunan Provincial Natural Science Foundation(Grant Nos.2024JJ6190 and 2024JK2007-1)。
文摘Chemical short-range order(SRO),a phenomenon at the atomic scale resulting from inhomogeneities in the local chemical environment,is usually studied using machine learning force field-based molecular dynamics simulations due to the limitations of experimental methods.To promote the reliable application of machine potentials in high-entropy alloy simulations,first,this work uses NEP models trained on two different datasets to predict the SRO coefficients of NbMoTaW.The results show that within the same machine learning framework,there are significant differences in the prediction of SRO coefficients for the Nb-Nb atomic pair.Subsequently,this work predicts the SRO coefficients of NbMoTaW using the NEP model and the SNAP model,both of which are trained on the same dataset.The results reveal significant discrepancies in SRO predictions for like-element pairs(e.g.,Nb-Nb and W-W)between the two potentials,despite the identical training data.The findings of this study indicate that discrepancies in the prediction results of SRO coefficients can arise from either the same machine learning framework trained on different datasets or different learning frameworks trained on the same dataset.This reflects possible incompleteness in the current training set's coverage of local chemical environments at the atomic scale.Future research should establish unified evaluation standards to assess the capability of training sets to accurately describe complex atomic-scale behaviors such as SRO.
文摘We present multi-threading and SIMD optimizations on short-range potential calculation kernel in Molecular Dynamics.For the multi-threading optimization,we design a partition-and-two-steps(PTS)method to avoid write conflicts caused by using Newton’s third law.Our method eliminates serialization bottle-neck without extra memory.We implement our PTS method using OpenMP.Afterwards,we discuss the influence of the cutoff if statement on the performance of vectorization in MD simulations.We propose a pre-searching neighbors method,which makes about 70%atoms meet the cutoff check,reducing a large amount of redundant calculation.The experiment results prove our PTS method is scalable and efficient.In double precision,our 256-bit SIMD implementation is about 3×faster than the scalar version.
文摘Using the asymptotic iteration method, we obtain the S-wave solution for a short-range three-parameter central potential with 1/r singularity and with a non-orbital barrier. To the best of our knowledge, this is the first attempt at calculating the energy spectrum for this potential, which was introduced by H. Bahlouli and A. D. Alhaidari and for which they obtained the “potential parameter spectrum”. Our results are also independently verified using a direct method of diagonalizing the Hamiltonian matrix in the J-matrix basis.
基金Project supported by the National Natural Science Foundation of China(Grant Nos10774093 and 10374061)
文摘This paper presents recurrence spectra of highly excited lithium atoms with M = 1 state in parallel electric and magnetic fields at a fixed scaled energy ε = -0.03. Short-ranged potentials including ionic core potential and centrifugal barrier are taken into account. Their effects on the states and photo-absorption spectrum are analysed in detail. This demonstrates that the geometric features of classical orbits are of special importance for modulations of the spectral pattern. Thus the weak polarization as well as the reduction of correlation of electrons induced by short-ranged potentials give rise to the recurrence spectra of lithium M = 1 atoms more compact than that of the M = 0 one, which is in good agreement with the experimental prediction.
基金supported by the National Natural Science Foundation of China(No.21803053)the Natural Science Foundation of Zhejiang Province,China(No.LY20B030005)the Open Project Fund of Key Laboratory of Excited-State Materials of Zhejiang Province。
文摘Global optimization of Morse clusters with shortrange potential is a great challenge.Here,we apply our recently developed unbiased fuzzy global optimization method to systematically study Morse clusters with the potential rangeρ=14 and the number of atoms N up to 400.All the putative global minima reported in the literature have been successfully reproduced with relatively high success ratios.Compared to the available results for N≤240 and several larger Morse clusters,new global minima(and local minima)with lower energies have been found out for N=164,175,188,193,194,197,239,246,260,318,and 389.Clusters with magic numbers are figured out through fitting the size-dependent global minimum energies.The cluster structures tend to be close-packed for short-range potential with large N.
文摘Quantum energies which are used in applications are usually composed of repulsive and attractive terms. The objective of this study is to use an accurate and efficient fitting of the repulsive energy instead of using standard parametrizations. The investigation is based on Density Functional Theory and Tight Binding simulations. Our objective is not only to capture the values of the repulsive terms but also to efficiently reproduce the elastic properties and the forces. The elasticity values determine the rigidity of a material when some traction or load is applied on it. The pair-potential is based on an exponential term corrected by B-spline terms. In order to accelerate the computations, one uses a hierarchical optimization for the B-splines on different levels. Carbon graphenes constitute the configurations used in the simulations. We report on some results to show the efficiency of the B-splines on different levels.
文摘Quantum energies which are used in applications are usually composed of repulsive and attractive terms. The objective of this study is to use an accurate and efficient fitting of the repulsive energy instead of using standard parametrizations. The investigation is based on Density Functional Theory and Tight Binding simulations. Our objective is not only to capture the values of the repulsive terms but also to efficiently reproduce the elastic properties and the forces. The elasticity values determine the rigidity of a material when some traction or load is applied on it. The pair-potential is based on an exponential term corrected by B-spline terms. In order to accelerate the computations, one uses a hierarchical optimization for the B-splines on different levels. Carbon graphenes constitute the configurations used in the simulations. We report on some results to show the efficiency of the B-splines on different levels.
文摘Classical scattering trajectories are known to form a Lagrangian manifold in euclidean phase space,which allows the classification of local focal points for sufficiently small dimensions. For the case of a short-range potential, we show that the natural description of focal points at infinity is a Lagrangian manifold in the cotangent bundle of the sphere and establish the relationship between focal points at infinity and the projection singularities of that manifold.
文摘We study the confinement of a spinless charged particle to a spherical quantum dot under the influence of a linear electric field.The spherical quantum dot is described by a short-range potential given by the power-exponential potential.Then,by analysing the region near the spherical quantum dot centre,we discuss two cases where the energy levels can be obtained for s-waves and how the linear electric field modifies the spectrum of energy of the spherical quantum dot.
基金financially supported by the National Natural Science Foundation of China(Nos.52071024,52271003 and 52101188)the National Science Fund for Distinguished Young Scholars of China(No.52225103)+2 种基金the Funds for Creative Research Groups of China(No.51921001)the Projects of International Cooperation and Exchanges NSFC(Nos.51961160729 and 52061135207)the Fundamental Research Fund for the Central Universities of China,and the State Key Laboratory for Advanced Metals and Materials.
文摘High/medium entropy alloys(H/MEAs)have shown unique strengthening behavior and mechanical prop-erties because of the presence of massive local chemical orderings.Nevertheless,dynamic interactions between chemical short-range orders(CSROs)and dislocations,and the underlying atomic strengthening mechanism remain elusive.In this work,we first developed a novel machine learning-embedded atom method(ML-EAM)potential of the CoNiV system,trained on a comprehensive first-principles dataset,which enables accurate and efficient modeling of CSRO formation and dislocation dynamics.Then,we in-vestigated the strengthening mechanisms of CSROs in CoNiV MEA through machine learning-augmented molecular dynamics(MD)simulations.Hybrid MD/Monte Carlo simulations reveal that CSRO domains possess an L1_(2)(NiCo)_(3) V structure,whose size increases with lowering annealing temperatures.These domains significantly enhance strength by impeding dislocation motion through complex energy path-ways,increasing depinning forces,and reducing mobility.Moreover,the MD simulations combined with theoretical analysis elucidate the competition between CSRO-assisted strengthening(via antiphase bound-ary formation)and solid solution weakening(via reduced atomic misfit volume).Phonon-drag effects are also amplified by CSROs,further resisting dislocation glide.Our results demonstrate that L1_(2)-CSROs strengthen CoNiV MEA primarily through antiphase boundary and phonon-drag contributions,providing new insights for designing high-performance multi-principal-element alloys via tailoring CSROs.
文摘We consider a pair of Hamiltonians （H, H0） on L2（R^n）, where H0=p^2 -x^2 is a SchrSdinger operator with a repulsive potential, and H = H0＋V（x）. We show that, under suitable assumptions on the decay of the electric potential, V is uniquely determined by the high energy limit of the scattering operator.
基金supported by the National Natural Science Foundation of China(Grant Nos.12274035,and 11874088)supported by the Fundamental Research Funds for the Central Universities。
文摘An accurate total energy calculation is essential in materials computation.To date,many tight-binding(TB)approaches based on parameterized hopping can produce electronic structures comparable to those obtained using first-principles calculations.However,TB approaches still have limited applicability for determining material properties derived from the total energy.That is,the predictive power of the TB total energy is impaired by an inaccurate evaluation of the repulsive energy.The complexity associated with the parametrization of TB repulsive potentials is the weak link in this evaluation.In this study,we propose a new method for obtaining the pairwise TB repulsive potential for crystalline materials by employing the Chen-Möbius inversion theorem.We show that the TB-based phonon dispersions,calculated using the resulting repulsive potential,compare well with those obtained by first-principles calculations for various systems,including covalent and ionic bulk materials and twodimensional materials.The present approach only requires the first-principles total energy and TB electronic band energy as input and does not involve any parameters.This striking feature enables us to generate repulsive potentials programmatically.
