Ab initio investigation of the CdTe (001) surface

Based on first-principles density functional pseudopotential calculations, the CdTe (001) surface has been studied. Atomic configurations of different reconstructions are obtained with good accuracy. This reveals the efficiency of the force calculation for the surface relaxations. The surface energies were calculated on relaxed surface slabs as a function of the chemical potentials. The main result is that the energy of the Te-terminated dimerized surface with a (2 × 1) reconstruction is larger than the Cd-terminated c(2×2) reconstructed surface. This is in agreement with what was suspected by the equilibrium model introduced by F. Tinjod et al., which explains the formation of the quantum dots in CdTe/ZnTe superstructures.     

First-principles study of the structural properties and magnetism of R2CoIn8 (R=Y,Pr, Nd,and Dy) intermetallic compounds

The electronic structure of R₂CoIn₈ (R=Y, Pr, Nd, and Dy) intermetallic compounds is calculated from first principles based on the density functional theory (DFT). The Kohn–Sham single-particle equations of the DFT are solved using two independent computational methods, namely APW+lo and FPLO. First the structural properties of Y₂CoIn₈ are studied. Good agreement of calculated equilibrium volume and c/a ratio with the experimental data is found. Also we minimize the forces at equilibrium volume and calculate the symmetry-free structural parameters of the space group P4/mmm for Y₂CoIn₈. In Y₂CoIn₈, the Co 3d states are almost fully occupied and situated below the Fermi level. We applied the fixed-spin-moment method and a stable paramagnetic ground state for Y₂CoIn₈ was found. Finally, the crystal field (CF) parameters were calculated for R=Pr, Nd and Dy from first principles. The microscopic tetragonal CF Hamiltonian was diagonalized and the obtained eigenvalues and eigenfunctions were used to predict the anisotropy of the magnetic properties of R₂CoIn₈ single crystals.     

Elastic and brittle properties of the B2-MgRE (RE = Sc,Y, Ce,Pr, Nd,Gd, Tb,Dy, Ho,Er) intermetallics

The brittle and elastic properties of the B2-MgRE (RE = Sc, Y, Ce, Pr, Nd, Gd, Tb, Dy, Ho, Er) intermetallics have been investigated using first-principles density functional calculations. The calculated equilibrium lattice constants and enthalpies of formation are in overall agreement with the available experiment and theoretical results. The related physical properties of those compounds are compared with that of ductile YCu. The Fermi energy occurs above a peak in the DOS for B2-MgRE intermetallics, whereas for ductile YCu the Fermi energy occurs near a minimum in the DOS. For B2-YCu, the partial density of states of d-states at the Fermi energy is low, while for B2-MgRE the RE d-states are partially occupied, indicating their important roles in the directional bonding for this material. The Cauchy pressure (C₁₂-C₄₄) and the ratio of bulk to shear modulus B/G are used to assess the brittle/ductile behavior of B2-MgRE and YCu compounds. It can be concluded that the B2-MgRE alloys have brittle behavior. MgSc is the most brittle, and MgHo is the least brittle amongst those alloys.     

Ab initio calculations of CO physisorption on ceria(1 1 1)

Applying the method of increments, we have performed MP2 and CCSD(T) calculations for the physisorption of CO on a cerium site on the ceria(1 1 1) surface. Our calculations predict an interaction energy of −0.28 eV. We have compared our calculations to previous CCSD(T) calculations for the physisorption of CO on a cerium site on the ceria(1 1 0) surface and found a difference in the interaction energy that is related to the different structure of the two surfaces. On the ceria(110) surface only 30% of the interaction energy originate from electron correlation effects, but on the ceria(111) surface almost the entire binding energy (80%) is due to electron correlation effects. Analyses of the interaction energy contributions show that most of the electron correlation part originates from the interaction of CO with the O ions in the topmost surface layer.     

Structural relative stabilities and pressure-induced phase transitions for lanthanide trihydrides REH3 (RE=Sm,Gd, Tb,Dy, Ho,Er, Tm,and Lu)

The structures, structural relative stabilities, pressure-induced phase transitions, and equations of state for lanthanide trihydrides REH₃ (RE=Sm, Gd, Tb, Dy, Ho, Er, Tm, and Lu) are systematically studied using ab initio calculations under a core state model (CSM). The obtained ground-state parameters, such as lattice constants and bulk modulus, agree well with the available data. Among the P6₃/mm, P3?c1, and P6₃cm structures, the P6₃cm structure is found to be the most stable structure for lanthanide trihydride via the comparison of the calculated total energies. With the help of Birch–Murnaghan equation of state, the structural transitions from hexagonal to cubic for REH₃ (RE=Sm, Gd, Ho, Er, and Lu) under pressure are affirmed; especially, the similar behavior of REH₃ (RE= Tb, Dy, and Tm) is reasonably predicted for the first time by this means. For the transitions, the repulsive interactions of H–H atoms may play an important role in terms of the analysis of the structures in the vicinity of the theoretical phase transition.     

Ab initio study of some fundamental properties of the M3X (M=Cr,V; X=Si,Ge) compounds

M₃X (M=Cr, V; X=Si, Ge) compounds are studied using first-principles calculations based on the Density Functional Theory (DFT). It is found that the bulk of Cr₃X (X=Si, Ge) compounds are comparable to those of Al₂O₃, the nearest-neighbor distance D_M−M and D_M−X in these compounds increase and the bulk modulus decrease, there is a strong interaction between M and M (M=Cr the interaction is stronger). Also the interaction between M (M=Cr, V) and X (X=Ge) is negative, an anti-bonding-type interaction is dominant between these atoms.     

Half-metallic magnetism of Co2CrX (X=As, Sb) Heusler compounds: An ab initio study

In this study, we present the electronic, magnetic, and structural properties of two novel half-metallic full-Heusler compounds, Co₂CrAs and Co₂CrSb, in cubic L2₁ geometry. The calculations are based on the density functional theory within plane-wave pseudopotential method and spin-polarized generalized gradient approximation of the exchange-correlation functional. The electronic band structures and density of states of the systems indicate half-metallic behavior with vanishing electronic density of states of minority spins at Fermi level, which yields perfect spin polarization. The calculated magnetic moments of both systems in L2₁ structure are 5.00 μ_B, which are largely localized on the chromium site. The energy gaps in minority spin states are restricted by the 3d-states of cobalt atoms on two different sublattices. The formation enthalpies for both structures are negative indicating stability of these systems against decomposition into stable solid compounds.     

Ab initio molecular dynamics study on Ag (n = 4, 5, 6)

Ab initio Molecular Dynamics (MD) method, based on density functional theory (DFT) with planewaves and pseudopotentials, was used to study the stability and internal motion in silver cluster Ag_n, with n =4-6. Calculations on the neutral, cationic and anionic silver dimer Ag₂ show that the bond distance and vibrational frequency calculated by DFT are of good quality. Simulations of Ag₄, Ag₅, and Ag₆ in canonical ensemble reveal distinct characteristics and isomerization paths for each cluster. At a temperature of 800 K, an Ag₄ has no definite structure due to internal motion, while for Ag₅ and Ag₆the clusters maintain the planar structure, with atomic rearrangement observed for Ag₅ but not for Ag₆. At a temperature of 200 K, Ag₄ can exist in two planar structures whilst Ag₅ is found to be stable only in the planar form. In contrast Ag₆ is stable in both planar trigonal and 3D pentagonal structures. Micro-canonical MD simulation was performed for all three clusters to obtain the vibrational density of states (DOS). Received 5 May 1999 and Received in final form 20 August 1999     

Ab initio chemisorption studies of H on Fe(110)

Ab initio configuration interaction calculations are performed to study the chemisorption of atomic H on a Fe(110) surface. The lattice is modeled as an embedded three-layer, 40-atom cluster with the Fe atoms fixed at the bulk position. Fe 3d orbitals are explicitly included on five Fe atoms on the surface. Hydrogen strongly binds to the Fe(110) surface at the long-bridge, short-bridge, and quasi three-fold sites. The calculated adsorption energies are 2.76, 2.73, and 2.71 eV, respectively. H-surface bonding at the on-top Fe site is more than 0.4 eV weaker. The calculated H-surface distances are 0.89, 1.03, and 0.87 Å for H at the long-bridge, short-bridge, and quasi three-fold sites, respectively, which agrees well with the LEED value of 0.9 ± 0.1 Å. The H-surface stretching vibrational frequencies are calculated to be 1070, 1066, and 1073 cm⁻¹, at the long-bridge, short-bridge, and quasi three-fold sites, respectively. The work function of Fe(110) decreases on H adsorption. The present calculations indicate that H diffusion into the bulk through the short-bridge site will have a much higher activation barrier than via the long-bridge and quasi three-fold sites.     

Ab initio calculations of B<Subscript>2</Subscript> type RHg (R = Ce,Pr, Eu and Gd) intermetallic compounds

Spin polarized ab initio calculations have been carried out to study the structural, electronic, elastic and thermal properties of RHg (R = Ce, Pr, Eu and Gd) intermetallic compounds in B₂ structure. The calculations have been performed by using both generalized gradient approximation (GGA) and local spin density approximation (LSDA). The calculated value of lattice constant (a ₀) for these compounds with GGA is in better agreement with the experimental data than those with LSDA. Bulk modulus (B), first-order pressure derivative of bulk modulus and magnetic moment (μ _B ) are also presented. The energy band structure and electron density of states show the occupancy of 4f states for light as well as heavy rare earth atom. The elastic constants are predicted from which all the related mechanical properties like Poisson’s ratio (σ), Young’s modulus (E), shear modulus (G _H ) and anisotropy factor (A) are calculated. The ductility or brittleness of these compounds is predicted from Pugh’s rule (B/G _H ) and Cauchy pressure (C ₁₂ ? C ₄₄). The Debye temperature (θ _D ) is estimated from the average sound velocity, which have not been calculated and measured yet.     

Magnetic properties of RE2Au compounds (RE = Pr,Nd, Gd,Tb, Dy,Ho, Er,Tm and Y)

Magnetic properties of nine RE₂Au compounds have been studied in fields of up to 19 kOe in the temperature range 4.2K–300K. It has been found that all compounds are paramagnetic at room temperature except Gd₂Au. The compounds with Pr, Nd, Ho, Er and Tm exhibit Curie-Weiss behaviour with paramagnetic moments in close agreement with those expected for the free RE³⁺ ion. The moment of gold was found to be zero. The compounds with Pr, Nd, Tb, Dy, Er and Tm are antiferromagnetic at low temperatures. It appears that Ho₂Au is ferromagnetically ordered below 4.5 K. No evidence for magnetic ordering was found for Y₂Au. The compound with Tb exhibits metamagnetic behaviour.     

The third-order elastic moduli and pressure derivatives for AlRE (RE=Y,Pr, Nd,Tb, Dy,Ce) intermetallics with B2-structure: A first-principles study

The third-order elastic moduli and pressure derivatives of the second-order elastic constants of novel B2-type AlRE (RE=Y, Pr, Nd, Tb, Dy, Ce) intermetallics are presented from first-principles calculations. The elastic moduli are obtained from the coefficients of the polynomials from the nonlinear least-squares fitting of the energy–strain functions. The calculated second-order elastic constants of AlRE intermetallics are consistent with the previous calculations. To judge that our computational accuracy is reasonable, the calculated third-order constants of Al are compared with the available experimental data and other theoretical results and are found to have good agreements. In comparison with the theory of the linear elasticity, the third-order effects are very important with the finite strains which are lager than approximately 3.5%. Finally, the pressure derivatives have been discussed.     

Ab initio comparative study of C54 and C49 TiSi2 surfaces

A theoretical comparison of C54 and C49 TiSi₂ surfaces is presented, using ab initio plane-wave ultrasoft pseudopotential method based on generalized gradient approximation (GGA). The different surface energies of TiSi₂ have not only been calculated out, but the preferential formation of C49 phase in solid-state reaction could be explained by smaller surface energies and Poisson's ratio of C49 TiSi₂ as well. As for polar C54 TiSi₂(1 0 0) and C49 TiSi₂(0 1 0) surfaces, the Si termination surfaces are more stable.     

He+2团簇结构和解析势能函数的从头计算研究

采用从头计算的耦合簇方法CCSD(T)和He原子Dunning's相关调和基函数组对He+2团簇的结构参数、势能曲线进行计算.利用Murrell-Sorbie函数和最小二乘法拟合出了解析势能函数,并以此为基础计算出光谱常数.通过比较发现:分子结构和光谱常数计算结果均与实验值符合良好,优于文献报道的结果.说明本文所得势能函数解析表达式准确反映了分子中原子间相互作用,可用于更广泛的研究.     

Comment on: Ab initio calculations of B<Subscript>2</Subscript> type RHg (R = Ce,Pr,Eu and Gd) intermetallic compounds

In a recent article by Devi et al. [Eur. Phys. J. B 87, 268 (2014)], thestructural, electronic, elastic and some thermal properties of B2 type RHg(R = Ce, Pr, Eu and Gd)intermetallic compounds have been studied by ab initio calculations. After the study oftheir article I found that there are some mistakes in predicted crystal density,longitudinal, transverse and average elastic wave velocities, and Debye temperature data.The crystal density has been found multiplied per 4. Also the longitudinal, transverse andaverage elastic wave velocities and Debye temperature are different from my reexaminedvalues (all results represented by Devi et al. have been found divided per 2). Althoughthese small mistakes do not influence their conclusion, it is better to correct them. Inthe present work, I reexamined all data again by using the right formulas, based on thelattice parameters and the elastic constants obtained in the work of Devi et al.     

He_2~+团簇结构和解析势能函数的从头计算研究

采用从头计算的耦合簇方法CCSD(T)和He原子Dunning’s相关调和基函数组对He+2 团簇的结构参数、势能曲线进行计算。利用Murrell Sorbie函数和最小二乘法拟合出了解析势能函数 ,并以此为基础计算出光谱常数。通过比较发现 :分子结构和光谱常数计算结果均与实验值符合良好 ,优于文献报道的结果。说明本文所得势能函数解析表达式准确反映了分子中原子间相互作用 ,可用于更广泛的研究     

Magnetic properties and structure of hydrogen-amorphized intermetallic compounds RFe2Hx (R=Y, Gd, Tb, Dy, Ho, Er)

The magnetic properties and structure of hydrogen-amorphized RFe₂H_x compounds (R=Y, Gd, Tb, Dy, Ho, Er) have been studied. It is shown that amorphization of the RFe₂H_x hydrides results in an increase of Fe-Fe, and a decrease of R-Fe exchange interaction energy compared to their crystalline hydride counterparts. The magnetic structure of amorphous RFe₂H_x compounds, with the exception of those with R=Y and Gd is apparently noncollinear ferrimagnetic, as in the crystalline hydrides. A model of the heterophase state of amorphous RFe₂H_x hydrides based on x-ray diffraction and magnetic data is proposed. Fiz. Tverd. Tela (St. Petersburg) 39, 908–912 (May 1997)     

第一性原理计算Fe从bcc到hcp结构的相变路径及其磁性边界

本文采用基于密度泛函理论的第一性原理方法,计算了压力作用下Fe从bcc到hcp结构相变的势能面、相变路径以及相变过程中的磁性相边界.结果表明:与Burgers路径不同,相变过程中bcc结构(110)bcc面的剪切和相对滑移相互耦合,并伴随有(110)bcc面间距的减小;这一相变机制可以解释Fe的静高压实验中在相变初期观察到的hcp结构异常.因此,并不需要像Wang和Ingalls提出的那样,在相变过程中引入一个亚稳定的fcc相来解释这些实验结果.对相变势能面的计算表明剪切对相变的发生有激活作用.此外,分析表明相变过程中涉及复杂的磁性转变,相变过渡态位置正好位于磁性相边界上,并对原子磁性对结构转变影响的物理机制进行了讨论.