六角密堆结构钇的解析型键序作用势

为了构建拟合势需要的数据库,采用密度泛函理论方法计算了六角密堆结构钇(hcp-Y)的晶格参数、弹性常数、内聚能、结构能差以及相关的点、面缺陷性质. 基于解析型键序作用势,构建了hcp-Y的多体作用势模型. hcp-Y势模型是通过拟合Y的晶格参数、弹性常数、体弹模量、内聚能、空位形成能和不同相之间的结构能差而构建.分析发现,所得到的势模型能够很好地描述hcp-Y的自填隙原子形成能、空位形成能、双空位键能以及其它体性质,同时,构建的势模型用来研究Y的热动力学性质的相关结果也比较理想.     

Bayesian ensemble approach to error estimation of interatomic potentials

Using a Bayesian approach a general method is developed to assess error bars on predictions made by models fitted to data. The error bars are estimated from fluctuations in ensembles of models sampling the model-parameter space with a probability density set by the minimum cost. The method is applied to the development of interatomic potentials for molybdenum using various potential forms and databases based on atomic forces. The calculated error bars on elastic constants, gamma-surface energies, structural energies, and dislocation properties are shown to provide realistic estimates of the actual errors for the potentials.     

Construction of embedded-atom-method interatomic potentials for alkaline metals （Li,Na,and K） by lattice inversion

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.     

Thermodynamic properties and structural stability of thorium dioxide

Using density functional theory (DFT) calculations, we have systematically investigated the thermodynamic properties and structural stabilities of thorium dioxide (ThO(2)). Based on the calculated phonon dispersion curves, we have calculated the thermal expansion coefficient, bulk modulus, and heat capacities at different temperatures for ThO(2) under the quasi-harmonic approximation. All the results are in good agreement with corresponding experiments proving the validity of our methods. Our theoretical studies can aid a clearer understanding of the thermodynamic behaviors of ThO(2) at different temperatures. In addition, we have also studied possible defect formations and diffusion behaviors of helium in ThO(2), to discuss its structural stability. It is found that in intrinsic ThO(2) without any Fermi energy shifts, the interstitial Th(i)(4+) defect rather than oxygen or thorium vacancies, interstitial oxygen, or any kinds of Frenkel pairs, is the most probable to form with an energy release of 1.74 eV. However, after upshifting the Fermi energy, the formation of the other defects also becomes possible. Regarding helium diffusion, we find that only through the thorium vacancy can it occur with the small energy barrier of 0.52 eV. Otherwise, helium atoms can hardly incorporate or diffuse in ThO(2). Our results indicate that upward shifting of the Fermi energy of ThO(2) should be prevented to avoid the formation of thorium vacancies so as to avert helium caused damage.     

A potential for Th from inversion of cohesive energy: Elastic constants

An interatomic pair potential for Th was derived by using the Chen–Mobius lattice inversion of cohesive energy for fcc Th as a starting point to develop a free-parameter potential suitable to be used in molecular dynamic calculations for predicting microstructure evolution and thermal properties in multicomponent nuclear fuel. The cohesive energy versus lattice parameter of Th was computed from first principles electronic structure calculations. The elastic constants for fcc Th were calculated by applying different types of strain to the starting crystal. Based on this information, the shear modulus, the Youngs modulus and the Poissons ratio were obtained. The computed elastic constants of fcc Th are found to be in a good agreement with experiments and previous theoretical results.     

Atomistic simulation of ferroelectric–ferroelastic gadolinium molybdate

Gadolinium molybdate Gd₂(MoO₄)₃ orthorhombic ferroelectric ferroelastic (β'-phase) is simulated by the method of interatomic potentials. The simulation is performed using the GULP 4.0.1 code (General Utility Lattice Program), which is based on the minimization of the energy of the crystal structure. Parameters of the gadolinium–oxygen interatomic interaction potentials are determined by fitting to the experimental structural data and elastic constants by a procedure available in the GULP code. Atomistic modeling using the effective atomic charges and the system of interatomic potentials made it possible to obtain reasonable estimates of structural parameters, atomic coordinates, and the most important physical, mechanical, and thermodynamic properties of these crystals. Temperature dependences of the heat capacity and vibrational entropy of the crystal are obtained. The calculated parameters of gadolinium–oxygen interaction potentials can be used to simulate more complex gadolinium-containing compounds.     

Molecular Dynamics Simulations of the Elastic Anisotropy of Pd at Extreme Conditions

It is very interesting to discover the elastic properties of engineering material palladium, especially its elastic anisotropy along Hugoniot states. We here investigate the evolution of its high pressure and temperature(PT) elastic ansotropy along Hugoniot using molecular dynamics simulations based on accurate classical interatomic potential. In order to testify the validity of the interatomic potential of Pd in describing the high PT elastic properties, we calculate its isothermal and adiabatic elastic moduli using molecular dynamics method. The obtained data are in good agreement with experimental data. From the isothermal elastic constants, we deduce the Hugoniot acoustic velocities and find that the resulting data are in good agreement with experimental acoustic velocity data. Based on the reliable elastic constants, we further investigate the spacial elastic ansotropy along Hugoniot PT states. It is found that the spacial elastic anisotropy of Pd increases along Hugoniot states.     

Interatomic potential for uranium in a wide range of pressures and temperatures

Using the force-matching method we develop an interatomic potential that allows us to study the structure and properties of α-U, γ-U and liquid uranium. The potential is fitted to the forces, energies and stresses obtained from ab initio calculations. The model gives a good comparison with the experimental and ab initio data for the lattice constants of α-U and γ-U, the elastic constants, the room-temperature isotherm, the normal density isochore, the bond-angle distribution functions and the vacancy formation energies. The calculated melting line of uranium at pressures up to 80 GPa and the temperature of the α-γ transition at 3 GPa agree well with the experimental phase diagram of uranium.     

Lattice static properties of vacancy clusters and interstitials in hcp magnesium: Computer simulation studies

Computer simulation studies have been made to investigate the static properties of mono-, di- and tri-vacancy clusters and of self-interstitials in hcp magnesium in different configurations. Three interatomic potentials have been chosen for which the results have been compared. A crystallite containing about 1500 atoms and a model with the interatomic interaction extending upto the fourth neighbour distance have been used. Relaxation field, defect relaxation and formation energies, strength dipole tensors and relative changes in volume in the above defects have been computed and our final results compared with those of earlier workers. The formation energies of the defects are highly sensitive to the choice of the potential whose detailed structure guides the nature of relaxation and the dipole tensors. Calculations have been done for octahedral, tetrahedral and dumb-bell interstitials of which the last is found to be the most stable.     

金红石结构TiO2晶体点缺陷形成能的经验途径计算

基于经验参数化途径 ,通过对晶体结构、晶格形成能、介电性质和弹性实验数据拟合确定金红石结构氧化物晶体TiO2 的电子壳模型参数和非Coulomb互作用势参数 .计算点缺陷形成能 .论证Schottky缺陷是金红石结构TiO2中的本征缺陷 .     

Study of elastic properties of CeO2 by statistical moment method

The purpose of the present paper is to investigate the temperature and pressure dependences of the elastic properties of cerium dioxide using the statistical moment method (SMM). The equation of states of bulk CeO₂ is derived from the Helmholtz free energy, and the pressure dependences of the elastic moduli like the bulk modulus, B_T, shear modulus, G, Young’s modulus, E, and elastic constants (C11, C12, and C44) are presented taking into account the anharmonicity effects of the thermal lattice vibrations. In the present study, the influence of temperature and pressure on the elastic moduli and elastic constants of CeO₂ has also been studied, using three different interatomic potentials. We compare the results of the present calculations with those of the previous theoretical calculations as well as with the available experiments.     

A continuum thermal stress theory for crystals based on interatomic potentials

This paper presents a new continuum thermal stress theory for crystals based on interatomic potentials.The effect of finite temperature is taken into account via a harmonic model.An EAM potential for copper is adopted in this paper and verified by computing the effect of the temperature on the specific heat,coefficient of thermal expansion and lattice constant.Then we calculate the elastic constants of copper at finite temperature.The calculation results are in good agreement with experimental data.The thermal stress theory is applied to an anisotropic crystal graphite,in which the Brenner potential is employed.Temperature dependence of the thermodynamic properties,lattice constants and thermal strains for graphite is calculated.The calculation results are also in good agreement with experimental data.     

稀土金属间化合物RFe2Zn20-xInx的结构属性研究

RFe2Zn20(R代表稀土元素)型稀土金属间化合物因其低稀土含量和良好的铁磁性,已成为铁磁材料的研究热点之一.添加第四组元对该系列化合物的晶体结构和材料性能会产生一定影响.利用晶格反演方法获得了一系列有效的原子间相互作用势,对三元RFe2Zn20和四元RFe2Zn20-xInx化合物进行原子级模拟计算.研究表明,随着稀土元素原子量的增加,三元体系的晶格参数和体积呈线性下降,第四组元引入与否对该线性关系无直接影响.第四组元In替代Zn时,择优占据16c晶位,占满16c后选择占据96g晶位,始终不占据48f晶位.择优占位的结论符合实验观测,并与晶格反演势分析的结果一致.     

Elastic constants for noncentral interatomic potentials with crystal symmetry

The method of homogeneous deformation has been modified and applied to the calculation of second-order elastic constants for a class of noncentral interatomic potentials with the symmetry of Bravais lattices. The conditions of rotational invariance of the potential energy and the zero initial stress conditions were found and their validity was verified. For cubic symmetry a method was elaborated which applies to the noncentral interatomic potentials proposed byJohnson andWilson. The presented method of homogeneous deformation is compared with the original method ofFuchs and with the method of long waves in calculations of elastic constants.     

金红石结构TiO2晶体点缺陷形成能的经验途径计算

基于经验参数化途径,通过对晶体结构、晶格形成能、介电性质和弹性实验数据拟合确定金红石结构氧化物晶体TiO₂的电子壳模型参数和非Coulomb互作用势参数.计算点 缺陷形成能.论证Schottky缺陷是金红石结构TiO₂中的本征缺陷. 关键词： 2')" href="#">金红石TiO₂ 点缺陷形成能 电子壳模型 经验途径计算     

Cluster centers in forsterite laser crystals Mg2SiO4 : Cr and Mg2SiO4 : Cr : Li

The structures of chromium-containing cluster centers in forsterite laser crystals Mg₂SiO₄ : Cr and Mg₂SiO₄ : Cr : Li have been simulated using the interatomic potential method. A system of position-dependent parameters of interatomic interaction potentials in forsterite has been developed. In the ionic approximation, this system adequately describes the structure, properties, and defects of the crystal with correct representation of the preferred positional arrangement of chromium ions. The preferred mechanisms of chromium dissolution in forsterite crystals have been evaluated from a comparison of the energies of formation of chromium-containing clusters with different configurations. It has been demonstrated that the results of the simulation of interatomic interaction potentials are consistent with the experimental data obtained from electron paramagnetic resonance and optical spectroscopy studies.