Molecular dynamics(MD)is a powerful method widely used in materials science and solid-state physics.The accuracy of MD simulations depends on the quality of the interatomic potentials.In this work,a special class of e...Molecular dynamics(MD)is a powerful method widely used in materials science and solid-state physics.The accuracy of MD simulations depends on the quality of the interatomic potentials.In this work,a special class of exact solutions to the equations of motion of atoms in a body-centered cubic(bcc)lattice is analyzed.These solutions take the form of delocalized nonlinear vibrational modes(DNVMs)and can serve as an excellent test of the accuracy of the interatomic potentials used in MD modeling for bcc crystals.The accuracy of the potentials can be checked by comparing the frequency response of DNVMs calculated using this or that interatomic potential with that calculated using the more accurate ab initio approach.DNVMs can also be used to train new,more accurate machine learning potentials for bcc metals.To address the above issues,it is important to analyze the properties of DNVMs,which is the main goal of this work.Considering only the point symmetry groups of the bcc lattice,34 DNVMs are found.Since interatomic potentials are not used in finding DNVMs,they are exact solutions for any type of potential.Here,the simplest interatomic potentials with cubic anharmonicity are used to simplify the analysis and to obtain some analytical results.For example,the dispersion relations for small-amplitude phonon modes are derived,taking into account interactions between up to the fourth nearest neighbor.The frequency response of the DNVMs is calculated numerically,and for some DNVMs examples of analytical analysis are given.The energy stored by the interatomic bonds of different lengths is calculated,which is important for testing interatomic potentials.The pros and cons of using DNVMs to test and improve interatomic potentials for metals are discussed.Since DNVMs are the natural vibrational modes of bcc crystals,any reliable interatomic potential must reproduce their properties with reasonable accuracy.展开更多
The quest for sustainable energy solutions has intensified the search for alternative feedstocks that can supplement or replace fossil fuels. Obtaining fuels or chemicals through the conversion of renewable biomass is...The quest for sustainable energy solutions has intensified the search for alternative feedstocks that can supplement or replace fossil fuels. Obtaining fuels or chemicals through the conversion of renewable biomass is a promising candidate [1,2]. Some noblemetal-based (e.g., Pt, Pd and Rh) catalysts exhibit significant catalytic activity to the conversion reaction of these biomass.展开更多
To explore atomic-level phenomena in the Cu-Ni-Sn alloy,a second nearest-neighbor modified embedded-atom method(2NN MEAM)potential has been developed for the Cu-Ni-Sn system,building upon the work of other researchers...To explore atomic-level phenomena in the Cu-Ni-Sn alloy,a second nearest-neighbor modified embedded-atom method(2NN MEAM)potential has been developed for the Cu-Ni-Sn system,building upon the work of other researchers.This potential demonstrates remarkable accuracy in predicting the lattice constant,with a relative error of less than 0.5%when compared to density functional theory(DFT)results,and it achieves a 10%relative error in the enthalpy of formation compared to experimental data,marking substantial advancements over prior models.The bulk modulus is predicted with a relative error of 8%compared to DFT.Notably,the potential effectively simulates the processes of melting and solidification of Cu-15Ni-8Sn,with a simulated melting point that closely aligns with the experimental value,within a 7.5%margin.This serves as a foundation for establishing a 2NN MEAM potential for a flawless Cu-Ni-Sn system and its microalloying systems.展开更多
Molecular Dynamics(MD)simulation for computing Interatomic Potential(IAP)is a very important High-Performance Computing(HPC)application.MD simulation on particles of experimental relevance takes huge computation time,...Molecular Dynamics(MD)simulation for computing Interatomic Potential(IAP)is a very important High-Performance Computing(HPC)application.MD simulation on particles of experimental relevance takes huge computation time,despite using an expensive high-end server.Heterogeneous computing,a combination of the Field Programmable Gate Array(FPGA)and a computer,is proposed as a solution to compute MD simulation efficiently.In such heterogeneous computation,communication between FPGA and Computer is necessary.One such MD simulation,explained in the paper,is the(Artificial Neural Network)ANN-based IAP computation of gold(Au_(147)&Au_(309))nanoparticles.MD simulation calculates the forces between atoms and the total energy of the chemical system.This work proposes the novel design and implementation of an ANN IAP-based MD simulation for Au_(147)&Au_(309) using communication protocols,such as Universal Asynchronous Receiver-Transmitter(UART)and Ethernet,for communication between the FPGA and the host computer.To improve the latency of MD simulation through heterogeneous computing,Universal Asynchronous Receiver-Transmitter(UART)and Ethernet communication protocols were explored to conduct MD simulation of 50,000 cycles.In this study,computation times of 17.54 and 18.70 h were achieved with UART and Ethernet,respectively,compared to the conventional server time of 29 h for Au_(147) nanoparticles.The results pave the way for the development of a Lab-on-a-chip application.展开更多
Al,Ca,and Zn are representative commercial alloying elements for Mg alloys.To investigate the effects of these elements on the deformation and recrystallization behaviors of Mg alloys,we develop interatomic potentials...Al,Ca,and Zn are representative commercial alloying elements for Mg alloys.To investigate the effects of these elements on the deformation and recrystallization behaviors of Mg alloys,we develop interatomic potentials for the Al-Ca,Al-Zn,Mg-Al-Ca and Mg-Al-Zn systems based on the second nearest-neighbor modified embedded-atom method formalism.The developed potentials describe structural,elastic,and thermodynamic properties of compounds and solutions of associated alloy systems in reasonable agreement with experimental data and higher-level calculations.The applicability of these potentials to the present investigation is confirmed by calculating the generalized stacking fault energy for various slip systems and the segregation energy on twin boundaries of the Mg-Al-Ca and Mg-Al-Zn alloys,accompanied with the thermal expansion coefficient and crystal structure maintenance of stable compounds in those alloys.展开更多
One of the major tasks in a molecular dynamics （MD） simulation is the selection of adequate potential functions, from which forces are derived. If the potentials do not model the behaviour of the atoms correctly, th...One of the major tasks in a molecular dynamics （MD） simulation is the selection of adequate potential functions, from which forces are derived. If the potentials do not model the behaviour of the atoms correctly, the results produced from the simulation would be useless. Three popular potentials, namely, Lennard-Jones （L J）, Morse, and embedded-atom method （EAM） potentials, were employed to model copper workpiece and diamond tool in nanometric machining. From the simulation results and further analysis, the EAM potential was found to be the most suitable of the three potentials. This is because it best describes the metallic bonding of the copper atoms; it demonstrated the lowest cutting force variation, and the potential energy is most stable for the EAM.展开更多
Abstract The process of γ＇ phase precipitating from Ni75Al14MO11 is studied by a computational simulation technique based on microscopic phase-field kinetics model. We studied the phase transformation with the purpo...Abstract The process of γ＇ phase precipitating from Ni75Al14MO11 is studied by a computational simulation technique based on microscopic phase-field kinetics model. We studied the phase transformation with the purpose of clarifying the influence of the nearest interatomic potential V Ni-Al （the nearest interatomic potential） on the precipitation process of γ＇ phase. The result demonstrates that there are two kinds of ordered phases, respective Llo and L12 in the early stage, and Llo phase transforms into L12 phase subsequently. For L12 phase, Ni atoms mainly occupy α site （face center positions）, while Al and Mo atoms occupy fl sites （the vertex positions）. When VNi-Al is increased by 10 MeV, the occupation probability of Ni atoms on α sites and Al atoms on β sites are enhanced. Enhanced VNi-Al facilitates clustering and ordering of Al atom, which promotes the formation of the γ＇ phase. At last, the simulation result was discussed by employing the thermodynamic stability.展开更多
In this paper, we deduce the analytical form of many-body interatomic potentials based on the Green's function in tight-binding representation. The many-body potentials are expressed as the functions of the hoppin...In this paper, we deduce the analytical form of many-body interatomic potentials based on the Green's function in tight-binding representation. The many-body potentials are expressed as the functions of the hopping integrals which are the physical origin of cohesion of atoms. For thesimple case of s-valent system, the inversion of the many-body potentials are discussed in detail by using the lattice inversion method.展开更多
We investigate the tripartite entanglement dynamics of three two-level atoms in a multi-mode vacuum field. By considering the influences of the interatomic distance and the initial condition on the lower bound of conc...We investigate the tripartite entanglement dynamics of three two-level atoms in a multi-mode vacuum field. By considering the influences of the interatomic distance and the initial condition on the lower bound of concurrence and the tripartite negativity, we show that an optimal interatomic distance can be found to minimize the collective damping. Interestingly, at the same optimal distance, the tripartite entanglement would be maximized in the open dynamics process. In the case of shorter interatomic distance, the tripartite entanglement can display the oscillatory behavior in the initial short-time limit and be trapped in a stationary value in the long-time limit. In addition, the tripartite entanglement for the general situation with different interatomic distances is also discussed.展开更多
The cleavage force F(z) needed to separate parallel atomic planes by a distance z is first discussed for simple s-p metals using density functional theory.For the s-p nearly free-electron metals the linearized Thomas-...The cleavage force F(z) needed to separate parallel atomic planes by a distance z is first discussed for simple s-p metals using density functional theory.For the s-p nearly free-electron metals the linearized Thomas-Fermi equation is solved self-consistently in the cases of (a) semi-infinite planes of jellium (i.e. smeared uniform positive ions) and (b) a semi-infinite cylinder of finite radius, cleaved by a plane perpendicular to its axis. In (a), the elastic region has the form F(z)=Az ∝ Zrs-11/2, where rs is the mean interelectronic distance in the jellium model. Size effects are then considered, with possible relevance to atomic force microscopy.Defect energies are treated, using both electron theory and pair force laws.展开更多
In silico prediction of potential synthetic targets is the prerequisite for function-led discovery of new zeolites. Millions of hypothetical zeolitic structures have been predicted via various computational methods, b...In silico prediction of potential synthetic targets is the prerequisite for function-led discovery of new zeolites. Millions of hypothetical zeolitic structures have been predicted via various computational methods, but most of them are experimentally inaccessible under conventional synthetic conditions.Screening out unfeasible structures is crucial for the selection of synthetic targets with desired functions.The local interatomic distance（LID） criteria are a set of structure rules strictly obeyed by all existing zeolite framework types. Using these criteria, many unfeasible hypothetical structures have been detected. However, to calculate their LIDs, all hypothetical structures need to be fully optimized without symmetry constraints. When evaluating a large number of hypothetical structures, such calculations may become too computationally expensive due to the forbiddingly high degree of freedom. Here, we propose calculating LIDs among structures optimized with symmetry constraints and using them as new structure evaluation criteria, i.e., the LIDsymcriteria, to screen out unfeasible hypothetical structures. We find that the LIDsymcriteria can detect unfeasible structures as many as the original non-symmetric LID criteria do, yet require at least one order of magnitude less computation at the initial geometry optimization stage.展开更多
The lattice-inversion embedded-atom-method interatomic potential developed previously by us is extended to alkaline metals including Li,Na,and K.It is found that considering interatomic interactions between neighborin...The lattice-inversion embedded-atom-method interatomic potential developed previously by us is extended to alkaline metals including Li,Na,and K.It is found that considering interatomic interactions between neighboring atoms of an appropriate distance is a matter of great significance in constructing accurate embedded-atom-method interatomic potentials,especially for the prediction of surface energy.The lattice-inversion embedded-atom-method interatomic potentials for Li,Na,and K are successfully constructed by taking the fourth-neighbor atoms into consideration.These angular-independent potentials markedly promote the accuracy of predicted surface energies,which agree well with experimental results.In addition,the predicted structural stability,elastic constants,formation and migration energies of vacancy,and activation energy of vacancy diffusion are in good agreement with available experimental data and first-principles calculations,and the equilibrium condition is satisfied.展开更多
We applied an approach to the development of many-body interatomic potentials for NiZr alloys,gaining an improved accuracy and reliability.The functional form of the potential is that of the embedded method,but it has...We applied an approach to the development of many-body interatomic potentials for NiZr alloys,gaining an improved accuracy and reliability.The functional form of the potential is that of the embedded method,but it has been improved as follows. (1) The database used for the development of the potential includes both experimental data and a large set of energies of different structures of the alloys generated by Fab initio calculations. (2) The optimum parametrization of the potential for the given database is obtained by fitting. Using this approach we developed reliable interatomic potentials for Ni and Zr. The potential accurately reproduces basic equilibrium properties of the alloys.展开更多
In this paper, a mathematical relation was found between interatomic Hooke’s force constant and both the bulk modulus and interatomic distance in solid crystals, considering that the forces which have effect on an at...In this paper, a mathematical relation was found between interatomic Hooke’s force constant and both the bulk modulus and interatomic distance in solid crystals, considering that the forces which have effect on an atom are only those resulted from the neighboring atoms, and the forces are subject to Hooke’s law as the deflections of atoms from their equilibrium positions are very small. This work has been applied on some solid semiconducting crystals of diatomic primitive cell, including crystals of mono-atomic primitive cell automatically, by using linear statistical fitting with computer programming and, then, using mathematical analysis, proceeding from the vibrational dispersion relation of solid linear lattice, these two methods have been used in the process in order to support each other and for the result to be satisfying and reasonable. This is a contribution to the process of using computer programming in physics to facilitate mathematical analyses and obtain the required relations and functions by designing and developing appropriate computer programs in line with the macro and micro natures of materials. The importance of this is in enhancing our understanding of the interatomic actions in cells and of the crystal structure of materials in general and semiconductors in particular, as it is a step of the initial steps to facilitate the process of calculating energies and extracting mathematical relations between correlation energy and temperature as well as between sub-fusion and fusion energies with temperature.展开更多
Motivated by the special theory of gradient elasticity (GradEla), a proposal is advanced for extending it to construct gradient models for interatomic potentials, commonly used in atomistic simulations. Our focus is o...Motivated by the special theory of gradient elasticity (GradEla), a proposal is advanced for extending it to construct gradient models for interatomic potentials, commonly used in atomistic simulations. Our focus is on London’s quantum mechanical potential which is an analytical expression valid until a certain characteristic distance where “attractive” molecular interactions change character and become “repulsive” and cannot be described by the classical form of London’s potential. It turns out that the suggested internal length gradient (ILG) generalization of London’s potential generates both an “attractive” and a “repulsive” branch, and by adjusting the corresponding gradient parameters, the behavior of the empirical Lennard-Jones potentials is theoretically captured.展开更多
It is shown that electrons forming simple and multiple covalent bonds may have different contribu-tions to the interatomic interactions due to the degeneracy of electron states. A simple relationship between the lengt...It is shown that electrons forming simple and multiple covalent bonds may have different contribu-tions to the interatomic interactions due to the degeneracy of electron states. A simple relationship between the length of covalent bond, its order and atomic numbers of the interacting atoms is de-duced.展开更多
Experimental results of the temperature dependence of critical resolved. shear stresses (CRSS)of Mo, Fe, Al and Mg single crystals are shown. Associating reports in recent years, we point out that the approximate expo...Experimental results of the temperature dependence of critical resolved. shear stresses (CRSS)of Mo, Fe, Al and Mg single crystals are shown. Associating reports in recent years, we point out that the approximate exponential relationship between CRSS and the absolute temperatureat least in the region of the steep temperature dependence range of many materials is more common, even for bcc, fcc, and hcp single crystals. polycrystals and other covalent crystals,provided that the slip plane and slip direction are kept the same. Successful explanation with atomic force law shows that the interatomic forces (electronic structure) play a decisive role in determining the temperature dependence of yield stresses for a large number of materials.展开更多
The known multiferroics(MFs)-ternary oxides with perovskite-type structure PbB′_(0.5)B″_(0.5)O_(3),that undergo successive phase transitions(PTs),ferroelectric(FE)or antiferroelectric(AFE)--at the Curie temperature,...The known multiferroics(MFs)-ternary oxides with perovskite-type structure PbB′_(0.5)B″_(0.5)O_(3),that undergo successive phase transitions(PTs),ferroelectric(FE)or antiferroelectric(AFE)--at the Curie temperature,T_(C),and ferromagnetic(FM),antiferromagnetic(AFM)or ferrimagnetic at the Neel temperature,T_(N) and classical FEs and AFEs are considered.The dependences of the T_(C),T_(N) on the interatomic bond A-O strains in their perovskite structures have been constructed.On constructed dependencies,some ternary MFs are discovered,which have comparatively high temperatures of first FE or AFE and second magnetic PTs but their difference T_(C)-T_(N) values are high comparatively with the binary MFs.展开更多
Machine Learned Interatomic Potentials(MLIPs)combine the predictive power of Density Functional Theory(DFT)with the speed and scaling of interatomic potentials,enabling theoretical spectroscopy to be applied to larger...Machine Learned Interatomic Potentials(MLIPs)combine the predictive power of Density Functional Theory(DFT)with the speed and scaling of interatomic potentials,enabling theoretical spectroscopy to be applied to larger and more complex systems than is possible with DFT.展开更多
So far, it has been a challenge for existing interatomic potentials to accurately describe a wide range of physical properties and maintain reasonable efficiency. In this work, we develop an interatomic potential for ...So far, it has been a challenge for existing interatomic potentials to accurately describe a wide range of physical properties and maintain reasonable efficiency. In this work, we develop an interatomic potential for simulating radiation damage in body-centered cubic tungsten by employing deep potential, a neural network-based deep learning model for representing the potential energy surface. The resulting potential predicts a variety of physical properties consistent with first-principles calculations, including phonon spectrum, thermal expansion, generalized stacking fault energies, energetics of free surfaces, point defects, vacancy clusters, and prismatic dislocation loops. Specifically, we investigated the elasticity-related properties of prismatic dislocation loops, i.e., their dipole tensors, relaxation volumes, and elastic interaction energies. This potential is found to predict the maximal elastic interaction energy between two 1/2 <1 1 1> loops better than previous potentials, with a relative error of only 7.6%. The predicted threshold displacement energies are in reasonable agreement with experimental results, with an average of 128 eV. The efficiency of the present potential is also comparable to the tabulated gaussian approximation potentials and modified embedded atom method potentials, meanwhile, can be further accelerated by graphical processing units. Extensive benchmark tests indicate that this potential has a relatively good balance between accuracy, transferability, and efficiency.展开更多
基金support of the RSF Grant No.24-11-00139(analytics,numerical results,manuscript writing)Daxing Xiong acknowledges the support of the NNSF Grant No.12275116,the NSF Grant No.2021J02051,and the startup fund Grant No.MJY21035For Aleksey A.Kudreyko,this work was supported by the Bashkir StateMedicalUniversity StrategicAcademic Leadership Program(PRIORITY-2030)(analytics).
文摘Molecular dynamics(MD)is a powerful method widely used in materials science and solid-state physics.The accuracy of MD simulations depends on the quality of the interatomic potentials.In this work,a special class of exact solutions to the equations of motion of atoms in a body-centered cubic(bcc)lattice is analyzed.These solutions take the form of delocalized nonlinear vibrational modes(DNVMs)and can serve as an excellent test of the accuracy of the interatomic potentials used in MD modeling for bcc crystals.The accuracy of the potentials can be checked by comparing the frequency response of DNVMs calculated using this or that interatomic potential with that calculated using the more accurate ab initio approach.DNVMs can also be used to train new,more accurate machine learning potentials for bcc metals.To address the above issues,it is important to analyze the properties of DNVMs,which is the main goal of this work.Considering only the point symmetry groups of the bcc lattice,34 DNVMs are found.Since interatomic potentials are not used in finding DNVMs,they are exact solutions for any type of potential.Here,the simplest interatomic potentials with cubic anharmonicity are used to simplify the analysis and to obtain some analytical results.For example,the dispersion relations for small-amplitude phonon modes are derived,taking into account interactions between up to the fourth nearest neighbor.The frequency response of the DNVMs is calculated numerically,and for some DNVMs examples of analytical analysis are given.The energy stored by the interatomic bonds of different lengths is calculated,which is important for testing interatomic potentials.The pros and cons of using DNVMs to test and improve interatomic potentials for metals are discussed.Since DNVMs are the natural vibrational modes of bcc crystals,any reliable interatomic potential must reproduce their properties with reasonable accuracy.
文摘The quest for sustainable energy solutions has intensified the search for alternative feedstocks that can supplement or replace fossil fuels. Obtaining fuels or chemicals through the conversion of renewable biomass is a promising candidate [1,2]. Some noblemetal-based (e.g., Pt, Pd and Rh) catalysts exhibit significant catalytic activity to the conversion reaction of these biomass.
基金sponsored by the Science and Technology Foundation of Guizhou Provincial Education Department(No.QJJ[2024]60)Guizhou Provincial Basic Research Program(Natural Science)(No.QKHJC[2024]Youth 214)+1 种基金Science and Technology Foundation of Guizhou Minzu University(No.GZMUZK[2024]QD21)Research Projects of Anshun University(No.asxybsjj202413).
文摘To explore atomic-level phenomena in the Cu-Ni-Sn alloy,a second nearest-neighbor modified embedded-atom method(2NN MEAM)potential has been developed for the Cu-Ni-Sn system,building upon the work of other researchers.This potential demonstrates remarkable accuracy in predicting the lattice constant,with a relative error of less than 0.5%when compared to density functional theory(DFT)results,and it achieves a 10%relative error in the enthalpy of formation compared to experimental data,marking substantial advancements over prior models.The bulk modulus is predicted with a relative error of 8%compared to DFT.Notably,the potential effectively simulates the processes of melting and solidification of Cu-15Ni-8Sn,with a simulated melting point that closely aligns with the experimental value,within a 7.5%margin.This serves as a foundation for establishing a 2NN MEAM potential for a flawless Cu-Ni-Sn system and its microalloying systems.
文摘Molecular Dynamics(MD)simulation for computing Interatomic Potential(IAP)is a very important High-Performance Computing(HPC)application.MD simulation on particles of experimental relevance takes huge computation time,despite using an expensive high-end server.Heterogeneous computing,a combination of the Field Programmable Gate Array(FPGA)and a computer,is proposed as a solution to compute MD simulation efficiently.In such heterogeneous computation,communication between FPGA and Computer is necessary.One such MD simulation,explained in the paper,is the(Artificial Neural Network)ANN-based IAP computation of gold(Au_(147)&Au_(309))nanoparticles.MD simulation calculates the forces between atoms and the total energy of the chemical system.This work proposes the novel design and implementation of an ANN IAP-based MD simulation for Au_(147)&Au_(309) using communication protocols,such as Universal Asynchronous Receiver-Transmitter(UART)and Ethernet,for communication between the FPGA and the host computer.To improve the latency of MD simulation through heterogeneous computing,Universal Asynchronous Receiver-Transmitter(UART)and Ethernet communication protocols were explored to conduct MD simulation of 50,000 cycles.In this study,computation times of 17.54 and 18.70 h were achieved with UART and Ethernet,respectively,compared to the conventional server time of 29 h for Au_(147) nanoparticles.The results pave the way for the development of a Lab-on-a-chip application.
文摘Al,Ca,and Zn are representative commercial alloying elements for Mg alloys.To investigate the effects of these elements on the deformation and recrystallization behaviors of Mg alloys,we develop interatomic potentials for the Al-Ca,Al-Zn,Mg-Al-Ca and Mg-Al-Zn systems based on the second nearest-neighbor modified embedded-atom method formalism.The developed potentials describe structural,elastic,and thermodynamic properties of compounds and solutions of associated alloy systems in reasonable agreement with experimental data and higher-level calculations.The applicability of these potentials to the present investigation is confirmed by calculating the generalized stacking fault energy for various slip systems and the segregation energy on twin boundaries of the Mg-Al-Ca and Mg-Al-Zn alloys,accompanied with the thermal expansion coefficient and crystal structure maintenance of stable compounds in those alloys.
文摘One of the major tasks in a molecular dynamics （MD） simulation is the selection of adequate potential functions, from which forces are derived. If the potentials do not model the behaviour of the atoms correctly, the results produced from the simulation would be useless. Three popular potentials, namely, Lennard-Jones （L J）, Morse, and embedded-atom method （EAM） potentials, were employed to model copper workpiece and diamond tool in nanometric machining. From the simulation results and further analysis, the EAM potential was found to be the most suitable of the three potentials. This is because it best describes the metallic bonding of the copper atoms; it demonstrated the lowest cutting force variation, and the potential energy is most stable for the EAM.
基金financially supported by the National Natural Science Foundation of China (Nos. 51,204,147 and 51274175)International Cooperation Project Supported by Ministry of Science and Technology of China (No. 2014DFA50320)International Science and Technology Cooperation Project of Shanxi Province (Nos. 2013081017 and 2012081013)
文摘Abstract The process of γ＇ phase precipitating from Ni75Al14MO11 is studied by a computational simulation technique based on microscopic phase-field kinetics model. We studied the phase transformation with the purpose of clarifying the influence of the nearest interatomic potential V Ni-Al （the nearest interatomic potential） on the precipitation process of γ＇ phase. The result demonstrates that there are two kinds of ordered phases, respective Llo and L12 in the early stage, and Llo phase transforms into L12 phase subsequently. For L12 phase, Ni atoms mainly occupy α site （face center positions）, while Al and Mo atoms occupy fl sites （the vertex positions）. When VNi-Al is increased by 10 MeV, the occupation probability of Ni atoms on α sites and Al atoms on β sites are enhanced. Enhanced VNi-Al facilitates clustering and ordering of Al atom, which promotes the formation of the γ＇ phase. At last, the simulation result was discussed by employing the thermodynamic stability.
文摘In this paper, we deduce the analytical form of many-body interatomic potentials based on the Green's function in tight-binding representation. The many-body potentials are expressed as the functions of the hopping integrals which are the physical origin of cohesion of atoms. For thesimple case of s-valent system, the inversion of the many-body potentials are discussed in detail by using the lattice inversion method.
基金supported by the National Natural Science Foundation of China(Grant Nos.61178012,11204156,and 11304179)the Specialized Research Fund for the Doctoral Program of Higher Education,China(Grant Nos.20133705110001 and 20123705120002)the Natural Science Foundation of Shandong Province,China(Grant Nos.BS2013DX034,ZR2012FQ024,and ZR2014AP009)
文摘We investigate the tripartite entanglement dynamics of three two-level atoms in a multi-mode vacuum field. By considering the influences of the interatomic distance and the initial condition on the lower bound of concurrence and the tripartite negativity, we show that an optimal interatomic distance can be found to minimize the collective damping. Interestingly, at the same optimal distance, the tripartite entanglement would be maximized in the open dynamics process. In the case of shorter interatomic distance, the tripartite entanglement can display the oscillatory behavior in the initial short-time limit and be trapped in a stationary value in the long-time limit. In addition, the tripartite entanglement for the general situation with different interatomic distances is also discussed.
文摘The cleavage force F(z) needed to separate parallel atomic planes by a distance z is first discussed for simple s-p metals using density functional theory.For the s-p nearly free-electron metals the linearized Thomas-Fermi equation is solved self-consistently in the cases of (a) semi-infinite planes of jellium (i.e. smeared uniform positive ions) and (b) a semi-infinite cylinder of finite radius, cleaved by a plane perpendicular to its axis. In (a), the elastic region has the form F(z)=Az ∝ Zrs-11/2, where rs is the mean interelectronic distance in the jellium model. Size effects are then considered, with possible relevance to atomic force microscopy.Defect energies are treated, using both electron theory and pair force laws.
基金supported by the National Natural Science Foundation of China(Nos.21622102,21621001 and 21320102001)the National Key Research and Development Program of China(No.2016YFB0701100)
文摘In silico prediction of potential synthetic targets is the prerequisite for function-led discovery of new zeolites. Millions of hypothetical zeolitic structures have been predicted via various computational methods, but most of them are experimentally inaccessible under conventional synthetic conditions.Screening out unfeasible structures is crucial for the selection of synthetic targets with desired functions.The local interatomic distance（LID） criteria are a set of structure rules strictly obeyed by all existing zeolite framework types. Using these criteria, many unfeasible hypothetical structures have been detected. However, to calculate their LIDs, all hypothetical structures need to be fully optimized without symmetry constraints. When evaluating a large number of hypothetical structures, such calculations may become too computationally expensive due to the forbiddingly high degree of freedom. Here, we propose calculating LIDs among structures optimized with symmetry constraints and using them as new structure evaluation criteria, i.e., the LIDsymcriteria, to screen out unfeasible hypothetical structures. We find that the LIDsymcriteria can detect unfeasible structures as many as the original non-symmetric LID criteria do, yet require at least one order of magnitude less computation at the initial geometry optimization stage.
基金Project supported by the National Basic Research Program of China (Grant No. 2011CB606401)
文摘The lattice-inversion embedded-atom-method interatomic potential developed previously by us is extended to alkaline metals including Li,Na,and K.It is found that considering interatomic interactions between neighboring atoms of an appropriate distance is a matter of great significance in constructing accurate embedded-atom-method interatomic potentials,especially for the prediction of surface energy.The lattice-inversion embedded-atom-method interatomic potentials for Li,Na,and K are successfully constructed by taking the fourth-neighbor atoms into consideration.These angular-independent potentials markedly promote the accuracy of predicted surface energies,which agree well with experimental results.In addition,the predicted structural stability,elastic constants,formation and migration energies of vacancy,and activation energy of vacancy diffusion are in good agreement with available experimental data and first-principles calculations,and the equilibrium condition is satisfied.
基金Supported by the National Natural Science Foundation of China(No.2 9892 16 6 ,2 980 30 0 6 ,2 99830 0 1)
文摘We applied an approach to the development of many-body interatomic potentials for NiZr alloys,gaining an improved accuracy and reliability.The functional form of the potential is that of the embedded method,but it has been improved as follows. (1) The database used for the development of the potential includes both experimental data and a large set of energies of different structures of the alloys generated by Fab initio calculations. (2) The optimum parametrization of the potential for the given database is obtained by fitting. Using this approach we developed reliable interatomic potentials for Ni and Zr. The potential accurately reproduces basic equilibrium properties of the alloys.
文摘In this paper, a mathematical relation was found between interatomic Hooke’s force constant and both the bulk modulus and interatomic distance in solid crystals, considering that the forces which have effect on an atom are only those resulted from the neighboring atoms, and the forces are subject to Hooke’s law as the deflections of atoms from their equilibrium positions are very small. This work has been applied on some solid semiconducting crystals of diatomic primitive cell, including crystals of mono-atomic primitive cell automatically, by using linear statistical fitting with computer programming and, then, using mathematical analysis, proceeding from the vibrational dispersion relation of solid linear lattice, these two methods have been used in the process in order to support each other and for the result to be satisfying and reasonable. This is a contribution to the process of using computer programming in physics to facilitate mathematical analyses and obtain the required relations and functions by designing and developing appropriate computer programs in line with the macro and micro natures of materials. The importance of this is in enhancing our understanding of the interatomic actions in cells and of the crystal structure of materials in general and semiconductors in particular, as it is a step of the initial steps to facilitate the process of calculating energies and extracting mathematical relations between correlation energy and temperature as well as between sub-fusion and fusion energies with temperature.
文摘Motivated by the special theory of gradient elasticity (GradEla), a proposal is advanced for extending it to construct gradient models for interatomic potentials, commonly used in atomistic simulations. Our focus is on London’s quantum mechanical potential which is an analytical expression valid until a certain characteristic distance where “attractive” molecular interactions change character and become “repulsive” and cannot be described by the classical form of London’s potential. It turns out that the suggested internal length gradient (ILG) generalization of London’s potential generates both an “attractive” and a “repulsive” branch, and by adjusting the corresponding gradient parameters, the behavior of the empirical Lennard-Jones potentials is theoretically captured.
文摘It is shown that electrons forming simple and multiple covalent bonds may have different contribu-tions to the interatomic interactions due to the degeneracy of electron states. A simple relationship between the length of covalent bond, its order and atomic numbers of the interacting atoms is de-duced.
文摘Experimental results of the temperature dependence of critical resolved. shear stresses (CRSS)of Mo, Fe, Al and Mg single crystals are shown. Associating reports in recent years, we point out that the approximate exponential relationship between CRSS and the absolute temperatureat least in the region of the steep temperature dependence range of many materials is more common, even for bcc, fcc, and hcp single crystals. polycrystals and other covalent crystals,provided that the slip plane and slip direction are kept the same. Successful explanation with atomic force law shows that the interatomic forces (electronic structure) play a decisive role in determining the temperature dependence of yield stresses for a large number of materials.
基金financially supported by the Ministry of Science and Higher Education of the Russian Federation(State assignment in the field of scientific activity,Southern Federal University,2023,Project No.FENW-2023-0015).
文摘The known multiferroics(MFs)-ternary oxides with perovskite-type structure PbB′_(0.5)B″_(0.5)O_(3),that undergo successive phase transitions(PTs),ferroelectric(FE)or antiferroelectric(AFE)--at the Curie temperature,T_(C),and ferromagnetic(FM),antiferromagnetic(AFM)or ferrimagnetic at the Neel temperature,T_(N) and classical FEs and AFEs are considered.The dependences of the T_(C),T_(N) on the interatomic bond A-O strains in their perovskite structures have been constructed.On constructed dependencies,some ternary MFs are discovered,which have comparatively high temperatures of first FE or AFE and second magnetic PTs but their difference T_(C)-T_(N) values are high comparatively with the binary MFs.
基金funding from the EPSRC CDT in Modelling of Heterogeneous Systems funded by EP/S022848/1N.D.M.H.acknowledges support from EPSRC grant number EP/V000136/1Computing facilities were provided by the Scientific Computing Research Technology Platform of the University of Warwick through the use of the High Performance Computing(HPC)cluster Avon,and the Sulis Tier 2 platforms at HPC Midlands+funded by the Engineering and Physical Sciences Research Council(EPSRC),grant number EP/T022108/1.
文摘Machine Learned Interatomic Potentials(MLIPs)combine the predictive power of Density Functional Theory(DFT)with the speed and scaling of interatomic potentials,enabling theoretical spectroscopy to be applied to larger and more complex systems than is possible with DFT.
基金supported by the National Key Research and Development Program of China(Grant No.2022YFE03110000)the National Natural Science Foundation of China(Nos.52171084 and 12192282)the Foundation of President of Hefei Institutes of Physical Science,Chinese Academy of Sciences(Nos.YZJJQY202203 and BJPY2021A05).
文摘So far, it has been a challenge for existing interatomic potentials to accurately describe a wide range of physical properties and maintain reasonable efficiency. In this work, we develop an interatomic potential for simulating radiation damage in body-centered cubic tungsten by employing deep potential, a neural network-based deep learning model for representing the potential energy surface. The resulting potential predicts a variety of physical properties consistent with first-principles calculations, including phonon spectrum, thermal expansion, generalized stacking fault energies, energetics of free surfaces, point defects, vacancy clusters, and prismatic dislocation loops. Specifically, we investigated the elasticity-related properties of prismatic dislocation loops, i.e., their dipole tensors, relaxation volumes, and elastic interaction energies. This potential is found to predict the maximal elastic interaction energy between two 1/2 <1 1 1> loops better than previous potentials, with a relative error of only 7.6%. The predicted threshold displacement energies are in reasonable agreement with experimental results, with an average of 128 eV. The efficiency of the present potential is also comparable to the tabulated gaussian approximation potentials and modified embedded atom method potentials, meanwhile, can be further accelerated by graphical processing units. Extensive benchmark tests indicate that this potential has a relatively good balance between accuracy, transferability, and efficiency.
