Investigation of thermal expansion and compressibility of rare-earth orthovanadates using a dielectric chemical bond method

The chemical bond properties, lattice energies, linear expansion coefficients, and mechanical properties of ReVO 4 (Re = La, Ce, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb, Lu, Sc, Y) are investigated systematically by the dielectric chemical bond theory. The calculated results show that the covalencies of Re-O bonds are increasing slightly from La to Lu and that the covalencies of V-O bonds in crystals are decreasing slightly from La to Lu. The linear expansion coefficients decrease progressively from LaVO 4 to LuVO 4; on the contrary, the bulk moduli increase progressively. Our calculated results are in good agreement with some experimental values for linear expansion coefficients and bulk moduli.     

LaX(X=N,P,As,Sb)晶体的能带结构和化学键性质研究

用LMTO-ASA能带程序计算了LaX(X=N,P,As,Sb)晶体的能带结构,得到的晶体能隙分别为LaN_2.30eV,LaP_2.05eV,LaAs_1.66eV,LaSb_1.34eV,与实验结果基本相符.利用价电子总数在阴阳离子上的分配数之比,给出计算晶体化学键性质的经验关系式,根据该式计算晶体化学键的共价性与文献结果非常吻合,说明了该关系式的合理性.     

Dependence of charge transfer energy on crystal structure and composition in Eu3+-doped compounds

We report a method for estimating the positions of charge transfer (CT) bands in Eu3+-doped complex crystals. The environmental factor (h(e)) influencing the CT energy is presented. h(e) consists of four chemical bond parameters: the covalency, the bond volume polarization, the presented charge of the ligand in the chemical bond, and the coordination number of the central ion. These parameters are calculated with the dielectric theory of complex crystals. The relationship between the experimental CT energies and calculated environmental factors was established by an empirical formula. The calculated values are in good agreement with the experimental results. Such a relationship was confirmed by detailed analysis. In addition, our method is also useful to predict the charge-transfer position of any other rare earth ion.     

Density functional LCAO calculations for solids: Comparison among Hartree–Fock,DFT local density approximation,and DFT generalized gradient approximation structural properties

In this article, the results of a recently implemented DFT a posteriori and Kohn-Sham (KS ) linear combination of atomic orbital computational scheme for solids are presented. The equilibrium lattice parameters, bulk moduli, and lattice energies are calculated for eight crystallized systems. Local density approximation (LDA ) and generalized gradient approximation (GCA ) functionals and potentials are used. The maps of the Hartree-Fock (HF ) and Ks electronic densities and band structures are depicted. The KS results confirm the trend of the a posteriori scheme. Very good agreement between calculated and experimental lattice energies has been found for GGA potentials. © 1995 John Wiley & Sons, Inc.     

Cohesive properties of CeN and LaN from first principles

The effect of electron-correlation on the ground-state properties of CeN and LaN is studied by ab initio quantum-chemical methods. The approach which is used combines two separate steps: (1) the ground-state Hartree-Fock calculations for the crystal; (2) application of the method of increments to the studied system, which allows an expansion of bulk properties using the information from quantum-chemical calculations performed for finite clusters. As can be expected, for CeN correlation plays a significant role: with Hartree-Fock method only 49% of the experimental cohesive energy has been recovered, whereas after correlation corrections (coupled-cluster approach) the ground-state properties were found to be in good agreement with the experimental data found in literature. Thus, we obtained about 90% of the expected cohesive energy; the computed lattice constants and bulk moduli also agree well with the experimental values. For comparison, the equivalent treatment has been performed for LaN, where no f orbital is occupied. There the HF contribution to the ground-state properties is larger and hence the correlation effects weaker.     

Influence of electronic correlations on the ground-state properties of cerium dioxide

The electron-correlation effects on the ground-state properties of CeO(2) are studied by ab initio quantum-chemical methods. For this purpose the method of increments is applied. It combines Hartree-Fock calculations for periodic systems with correlation calculations requiring only information of the corresponding finite-cluster calculations. Using the coupled-cluster approach for the evaluation of the individual increments, we recover 93% of the experimental cohesive energy. The lattice constant and bulk modulus are found to be in good agreement with experimental values. For comparison also the results obtained with density functional methods are presented.     

Bond energy prediction of Curie temperature of lithium niobate crystals

A general expression of the Curie temperature (Tc) and spontaneous polarization (Ps) of lithium niobate (LN) crystals is energetically proposed by employing the viewpoint of the bond energy of constituent chemical bonds within the LN crystallographic frame. The calculated Tc values of various pure and doped LN crystals are in a good agreement with those reported data. Ps values of these LN crystals can also be quantitatively estimated in this work. It is found that the Li site is a sensitive lattice position to dominate the ferroelectricity of LN crystals. This novel method provides us a good understanding of ferroelectric behaviors of LN crystals, which may be applicable to the estimation of ferroelectric behaviors of LN-type solids.     

Ab initio all-electron periodic Hartree-Fock study of hydrostatic compression of pentaerythritol tetranitrate

We present a computational study of hydrostatic compression effects on the pentaerythritol tetranitrate (PETN) energetic material up to 22.7 GPa by means of the ab initio all-electron periodic Hartree-Fock quantum mechanical method with the STO-3G Gaussian basis set. We fitted the calculated volume-energy relation to the energy SJEOS polynomial function from which we obtained the compression dependence of the pressure (P), the bulk modulus (B), and its pressure derivative (B'). We also fitted the experimental volume-pressure relation to the pressure SJEOS polynomial function, which allowed us to calculate the experimental bulk modulus (B(exp)) and its pressure derivative (). Our calculated values, B = 6.73 GPa and B' = 24.63, are in reasonable agreement with the values B(exp) = 8.48 GPa and = 14.42 from our fit to the experimental X-ray data and with the value B(exp) = 9.8 GPa that was derived from the experimental elastic constants. In addition, we present a discussion on how the lattice vectors and the internal coordinates (i.e., bond lengths, bond angles, and torsion angles) of the C(CH(2)ONO(2))(4) molecules in the PETN lattice change during hydrostatic compression of the crystal. Our calculated results suggest that the C(CH(2)ONO(2))(4) molecules cannot be considered as being rigid but are in fact flexible, accommodating lattice compression through torsions, bendings in their bond angles, and contractions in their bond lengths. At pressures higher than about 8 GPa, however, both the C(CH(2)ONO(2))(4) molecules and the c lattice vector seem to stiffen somewhat. The a lattice vector does not exhibit this stiffening. As a consequence, the pressure dependence of the c/a ratio shows a minimum at about 8 GPa.     

High-pressure investigation of Li2MnSiO4 and Li2CoSiO4 electrode materials for lithium-ion batteries

In this work, the high-pressure behavior of Pmn2(1)-Li(2)MnSiO(4) and Pbn2(1)-Li(2)CoSiO(4) is followed by in situ X-ray diffraction at room temperature. Bulk moduli are 81 and 95 GPa for Pmn2(1)-Li(2)MnSiO(4) and Pbn2(1)-Li(2)CoSiO(4), respectively. Regardless of the moderate values of the bulk moduli, there is no evidence of any phase transformation up to a pressure of 15 GPa. Pmn2(1)-Li(2)MnSiO(4) shows an unusual expansion of the a lattice parameter upon compression. A density functional theory investigation yields lattice parameter variations and bulk moduli in good agreement with experiments. The calculated data indicate that expansion of the a lattice parameter is inherent to the crystal structure and independent of the nature of the transition-metal atom (M). The absence of pressure-driven phase transformation is likely associated with the incapability of the Li(2)MSiO(4) composition to adopt denser structures while avoiding large electrostatic repulsions.     

Density functional theory study of pyrophyllite and M-montmorillonites (M = Li, Na, K, Mg, and Ca): role of dispersion interactions

The stacking parameters, lattice constants, bond lengths, and bulk moduli of pyrophyllite and montmorillonites (MMTs), with alkali and alkali earth metal ions, are investigated using density functional theory with and without dispersion corrections. For pyrophyllite, it is found that the inclusion of the dispersion corrections significantly improves the agreement of the calculated values of the lattice parameters and bulk modulus with the experimental values. For the MMTs, the calculations predict that the interlayer spacing varies approximately linearly with the cation radius. The inclusion of dispersion corrections leads to sizable shifts of the interlayer spacings to shorter values. In Li-MMT, compaction of the interlayer distance triggers migration of the Li ion into the tetrahedral sheet and close coordination with basal oxygen atoms. Analysis of electron density distributions shows that the isomorphic octahedral Al(3+)/Mg(2+) substitution in MMT causes an increase of electron density on the basal oxygen atoms of the tetrahedral sheets.     

Accurate lattice energies of organic molecular crystals from periodic turbomole calculations

Accurate lattice energies of organic crystals are important i.e. for the pharmaceutical industry. Periodic DFT calculations with atom‐centered Gaussian basis functions with the Turbomole program are used to calculate lattice energies for several non‐covalently bound organic molecular crystals. The accuracy and convergence of results with basis set size and k‐space sampling from periodic calculations is evaluated for the two reference molecules benzoic acid and naphthalene. For the X23 benchmark set of small molecular crystals accurate lattice energies are obtained using the PBE‐D3 functional. In particular for hydrogen‐bonded systems, a sufficiently large basis set is required. The calculated lattice energy differences between enantiopure and racemic crystal forms for a prototype set of chiral molecules are in good agreement with experimental results and allow the rationalization and computer‐aided design of chiral separation processes. © 2018 Wiley Periodicals, Inc.     

Lattice energy estimation for inorganic ionic crystals

An empirical method based on chemical bond theory for the estimation of the lattice energy for ionic crystals has been proposed. The lattice energy contributions have been partitioned into bond dependent terms. For an individual bond, the lattice energy contribution made by it has been separated into ionic and covalent parts. Our calculated values of lattice energies agree well with available experimental and theoretical values for diverse ionic crystals. This method, which requires detailed crystallographic information and elaborate computation, might be extended and possibly yield further insights with respect to bond properties of materials.     

Low-compressibility and hard material carbon nitride imide C2N2(NH): First principles calculations

First principles calculations are performed to investigate the structural, mechanical, and electronic properties of C₂N₂(NH). Our calculated lattice parameters are in good agreement with the experimental data and previous theoretical values. Orthorhombic C₂N₂(NH) phase is found to be mechanically stable at an ambient pressure. Based on the calculated bulk modulus and shear modulus of polycrystalline aggregate, C₂N₂(NH) can be regarded as a potential candidate of ultra-incompressible and hard material. Furthermore, the elastic anisotropy and Debye temperatures are also discussed by investigating the elastic constants and moduli. Density of states and electronic localization function analysis show that the strong C-N covalent bond in CN₄ tetrahedron is the main driving force for the high bulk and shear moduli as well as small Poisson's ratio of C₂N₂(NH).     

First-principles calculations of hydrostatic pressure effects on the structural,elastic and thermodynamic properties of cubic monocarbides XC (X = Ti,V, Cr,Nb, Mo,Hf)

Six transition metal monocarbides (TiC, VC, CrC, NbC, MoC, HfC) with the rock-salt structure were chosen for a detailed comparative ab initio study of their structural, electronic, elastic, and thermodynamic properties at ambient and elevated up to 50 GPa hydrostatic pressures. Special attention was paid to the relation between the elastic and bonding properties and the number of valence electrons in each compound. Elastic anisotropy of the considered carbides was analyzed; the directions in the crystal lattice corresponding to the greatest and smallest Young's moduli values were identified. The calculated values of the elastic constants were used for further estimations of the Debye temperatures, Grüneisen parameters, specific heat capacities and linear coefficients of thermal expansion. Comparison of the calculated results with available experimental and theoretical data for TiC, VC, NbC, HfC yielded good agreement. The specific heat capacities and thermal expansion coefficients for CrC and MoC were calculated for the first time, to the best of the authors' knowledge.     

Thermodynamical study of the thermoelectric effect for magnesium silicide

The thermoelectric effect of magnesium silicide is studied by using a thermodynamical method in the presence of an electric field. The thermoelectric potential is evaluated from the partial derivative of free energy with respect to charge in which the free energy is calculated at the B3LYP/6-31G(d,p) level of density functional theory. This free energy is also utilized to determine the average dipole moment from which the polarizability, alpha; molar polarization, Psi; and dielectric constant can be computed. The present calculation for the dielectric constant (approximately 24-20) is in very good agreement with the experimental value (20). This accurate dielectric constant can be used to derive the relation of the thermoelectric potential with respect to temperature, from which the thermoelectric power or the Seebeck coefficients are calculated. The present result shows good agreement with experiment measurement for the Seebeck coefficients. In comparison, that calculation from the energy band structure theory is far off from the experimental values.     

NdP5O14晶体的晶场参数

NdP5O14晶体是一种高钕浓度的化学计量比的新型激光材料,它在高的Nd(III)浓度下荧光几乎不发生猝灭.本效计算了它的主要晶场参数,以深入了解它的光学性质.     

High pressure study on the phonon spectra and thermal properties in hafnium nitride and zirconium nitride

We report ab initio calculations of the thermal properties for transition metal nitrides, hafnium and zirconium nitride at ambient and high pressures. The assessment of thermodynamical properties like lattice specific heat, vibrational energy, internal energy and entropy for two nitrides has been carried out. The basic calculations of ingredient phonon density of states for the determination of thermal properties have been done using density functional perturbation theory including external perturbations like strains and electric fields in periodic systems. The ground state properties such as equilibrium lattice constants and bulk modulus obtained for two nitrides are in good agreement with the available experimental value. The calculated pressure variation of the phonon density of states shows trend similar to the experimental pressure dependent Raman spectra. The lattice specific heat, internal energy, entropy and Helmholtz energy increases with pressure.     

Calculation of the transport and relaxation properties of methane. I. Shear viscosity, viscomagnetic effects, and self-diffusion

Transport properties of pure methane gas have been calculated in the rigid-rotor approximation using the recently proposed intermolecular potential energy hypersurface [R. Hellmann et al., J. Chem. Phys. 128, 214303 (2008)] and the classical-trajectory method. Results are reported in the dilute-gas limit for shear viscosity, viscomagnetic coefficients, and self-diffusion in the temperature range of 80-1500 K. Compared with the best measurements, the calculated viscosity values are about 0.5% too high at room temperature, although the temperature dependence of the calculated values is in very good agreement with experiment between 210 and 390 K. For the shear viscosity, the calculations indicate that the corrections in the second-order approximation and those due to the angular-momentum polarization are small, less than 0.7%, in the temperature range considered. The very good agreement of the calculated values with the experimental viscosity data suggests that the rigid-rotor approximation should be very reasonable for the three properties considered. In general, the agreement for the other measured properties is within the experimental error.     

闪锌矿ZnTe高温高压下的弹性及热力学性质

利用第一性原理平面波赝势密度泛函理论并结合准谐德拜模型计算了闪锌矿结构ZnTe在高温高压下的弹性及热力学性质.得到了绝对零度、零压强时的晶格常数为0.6095 nm,仅比实验值(0.6103 nm)小0.1%.计算的体弹模量及弹性常数也与实验值符合较好.根据计算的高压下的弹性常数,得到其相变点约为10 GPa,与已知的实验值一致.通过准谐德拜模型,得到了常温下(T=300 K)的德拜温度为249 K,并得到了不同温度、不同压强下的热容.热容随着压强增加而减小,在高温、高压下,热容接近于Dulong-Petit极限.     

Overlapping shells model applied to diamondlike crystals

The model based on the assumption of the existence of an interatomic distance-dependent, local, effective crystal field applied to the alkaline metals (Int. J. Quantum Chem. 52, 321–328 (1994)) is modified and applied to the diamondlike structure crystals (C, Si, Ge, Sn). In the referred to model, a part of the electron density was missed—not included in the calculation (the density in the spaces between the shells). Such an approach could be used for the alkaline metals, but for the covalent crystals, this is a bad approximation. To avoid that problem, we assumed that the atom shells can overlap in such a way that the entire electron density is taken into the calculation. In this case, the electron density is “moved” from the outside of the shells mostly into the interatomic bond region. We applied the modified model to the calculation of the binding energy and the bulk modulus for the diamondlike crystals. The results show that well-chosen parameters allows one to reproduce the proper values of the binding energy at the equilibrium position. The bulk moduli calculated for these crystals are in quite good agreement with ones calculated as (regular crystal structure) B = 1/3(C₁₁ + 2 C₁₂), where C₁₁ and C₁₂ are elastic constants. © 1997 John Wiley & Sons, Inc.