First-principles investigation of the electronic structure and magnetism of eskolaite

The electronic structure and magnetism of eskolaite are studied by using first-principles calculations where the on-site Coulomb interaction and the exchange interaction are taken into account and the LSDA+U method is used.The calculated energies of magnetic configurations are very well fitted by the Heisenberg Hamiltonian with interactions in five neighbour shells; interaction with two nearest neighbours is found to be dominant. The Neel temperature is calculated in the spin-3/2 pair-cluster approximation. It is found that the measurements are in good agreement with for the values of U and J that are close to those obtained within the constrained occupation method.The band gap is of the Mott-Hubbard type.

Structural, electronic and mechanical properties of ReN2 from first principles

First principles calculations were performed to study the structural, electronic and mechanical properties of hypothetical rhenium dinitride ReN2 for various space groups. It was found that cubic Fm-3m and Pa-3, tetragonal P4(2)/mnm, and orthorhombic Pmmn structures are mechanically stable and metallic. P4(2)/mnm structure is thermodynamically stable at ambient conditions and up to 76 GPa. It has the shortest Re-N bond (1.964 angstrom).

Crystal structures and elastic properties of superhard IrN2 and IrN3 from first principles

First principles calculations were performed to investigate the structural, elastic, and electronic properties of IrN2 for various space groups: cubic Fm-3m and Pa-3, hexagonal P3(2)21, tetragonal P4(2)/mnm, orthorhombic Pmmn, Pnnm, and Pnn2, and monoclinic P2(1)/c. Our calculation indicates that the P2(1)/c phase with arsenopyrite-type structure is energetically more stable than the other phases. It is semiconducting (the remaining phases are metallic) and contains diatomic N-N with the bond distance of 1.414 A. These characters are consistent with the experimental facts that IrN2 is in lower symmetry and nonmetallic. Our conclusion is also in agreement with the recent theoretical studies that the most stable phase of IrN2 is monoclinic P2(1)/c. The calculated bulk modulus of 373 GPa is also the highest among the considered space groups. It matches the recent theoretical values of 357 GPa within 4.3% and of 402 GPa within 7.8%, but smaller than the experimental value of 428 GPa by 14.7%. Chemical bonding and potential displacive phase transitions are discussed for IrN2. For IrN3, cubic skutterudite structure (Im-3) was assumed.

Elastic and electronic properties of CoFe3N, RhFe3N, and IrFe3N from first principles

First principles calculations are performed to investigate the elastic and electronic properties of MFe3N (M=Co,Rh,Ir) at Pm-3m space group. The authors' calculation indicates that the three MFe3N phases are metallic and mechanically stable. For RhFe3N, the calculated lattice parameter of 3.826 A is in excellent agreement with the experimental value of 3.8292 A. The three phases are ferromagnetic with the calculated magnetic moments per f.u. being 8.92 mu(B) for CoFe3N, 9.04 mu(B) for RhFe3N, and 8.50 mu(B) for IrFe3N. The unusually large B/G ratio from 2.47 for CoFe3N and 2.45 for RhFe3N to 1.81 for IrFe3N indicates that they are ductile.

First-principles investigation on the elastic, magnetic and electronic properties of MFe3N (M = Fe, Ru, Os)

The elastic, magnetic and electronic properties of MFe3N (M = Fe, Ru, Os) are investigated via first-principles calculations. The calculated results are in agreement with the experimental and other theoretical data. The high ratios of bulk modulus to shear modulus 2.7, 2.0, and 1.8 for gamma'-Fe4N, RuFe3N, and OsFe3N, respectively, indicate that they have good ductility. gamma'-Fe4N possesses the largest B/C-44 (3.41) ratio, which suggests that it is much prone to shearing. The net magnetic moment per formula unit decreases from 9.90 for gamma'-Fe4N, 7.66 for RuFe3N, to 6.80 mu(B) for OsFe3N.

Orbital ordering in Cs2AgF4 from first principles

First-principles calculations using the augmented plane wave plus local orbital method, as implemented in the WIEN2K code, have been used to investigate the structural, electronic, and magnetic properties of the layered perovskite Cs2AgF4. Our calculations indicate that an orthorhombic ground state for Cs2AgF4 is energetically favored over tetragonal. We also find that Cs2AgF4 should be a strong two-dimensional ferromagnet, with very weak antiferromagnetic coupling between the layers, in agreement with the experiment. More importantly, an antiferrodistortive ordering of z(2)-x(2) and z(2)-y(2) orbitals is inferred from the density of states and from a spin density isosurface analysis.

地幔压力下CaSiO3 Perovskite状态方程和弹性常数的第一性原理研究

In the past decade density functional theory (DFT) has made its way from a peripheral position in quantum chemistry to center. Of course the often excellent accuracy of the DFT based methods has provided the primary driving force of this development. This dissertation is devoted to the study of physical and chemical properties of planetary materials by first-principle calculation. The concerned properties include the geometry, elastic constants and anisotropy. In the first chapter, we give a systematic introduction to theoretical background and review its progress. Development of quantum chemistry promotes the establishment of DFT. Theorem of Hohenberg-Kohn is the fundament of DFT and is developed to Kohn-Sham equation, which can be used to perform real calculations. Now, new corrections and extensions, together with developed exchange-correlation, have made DFT more accurate and suitable for larger systems. In the second chapter, we focus on the calculational methods and technical aspects of DFT. Although it is important to develop methods and program, external package are still often used. At the end of this chapter, we briefly some widely used simulation package and the application of DFT. In the third chapter, we begin to focus on properties of real materials by first principles calculation. We study a kind of minerals named Ca perovskite, investigate its possible structure and anisotropy at Earth’s mental condition. By understanding and predicting geo-physically important materials properties at extreme conditions, we can get the most accurate information to interpret seismic data in the context of likely geophysical processes.

压力和铁含量对铁方镁石(Mg1-xFex)O(x<=0.25)弹性常数影响的第一性原理研究

By first principle methods based on density functional theory (DFT),the equation of state(EOS) and elastic constants of both periclase and ferropericlase are calculated. The pressure and iron doping effects on the elastic constants of ferropericlase are investigated systematically. Firstly, we calculate the elastic constants of periclase and compare the obtained results with experimental data and other theoretical calculations, which shows a encouraging consistence and demonstrates the practicability of first-principle methods. Secondly, by adding iron into periclase crystal model, we build up ferropericlase with iron contents ranging from 0% to 25% mole percent. The corresponding elastic constants are calculated in a large pressure range(0~120GPa). Emphatically, the strong correlation of 3d electrons in transitional elements, such as iron, is difficult to treat in first-principle methods for a long time. The current solution is to make additional correction. During the initial stage of this study, the strong correlation of 3d electrons in iron is not considered, and we observed that addition of iron decreases the volume of ferropericlase, which is totally contradictory to the experimental data. By applying LDA+U approximation in order to solve the strongly correlated 3d electron of iron, we observed the expansion of volume by iron as expected. On the basis of the LDA+U approximation, the elastic constants of ferropericlase are calculated. After a detailed analysis of data obtained from theoretical calculations, we have reached the following conclusions:(1)pressure imposes positive effects on all elastic constants, and the degree of effects is C11>C12>C44. (2) Iron has no distinctive effects on C11 and C12, although some fluctuations are observed around 60GPa. However, iron has obvious softening effects on C44 The softening effects on C44 are intensified as pressure increases. Above the 100GPa, the effects increase greatly, even surpasses the pressure's positive effects in ferropericlase crystal models with iron mole percent of having 12.5%, 18.75% and 25% iron content. (3)As to the modulus deprived from elastic constants, iron has no effect on the adiabatic bulk module BS, only a little fluctuation around 60GPa. We find iron's softening effects on shear modulus G. (4)We find out that, compared with low iron content, elastic constants with iron content approaching 25mole% is consistently fluctuated，which may be caused by the limitations of the LDA+U approximation method itself. (5)We investigate the pressure and Fe doping effects on elastic anisotropy factor（A=(2C44+C12-C11)/C11） of ferropericlase and find out that iron contents will lower the critical isotropic pressure. At the same pressure, when the pressure is below the isotropic pressure, iron softens the anisotropy factor ; when pressure surpasses the isotropic pressure, iron increases the anisotropy factor.

Comprehensive search for new phases and compounds in binary alloy systems based on platinum-group metals, using a computational first-principles approach

We report a comprehensive study of the binary systems of the platinum-group metals with the transition metals, using high-throughput first-principles calculations. These computations predict stability of new compounds in 28 binary systems where no compounds have been reported in the literature experimentally and a few dozen of as-yet unreported compounds in additional systems. Our calculations also identify stable structures at compound compositions that have been previously reported without detailed structural data and indicate that some experimentally reported compounds may actually be unstable at low temperatures. With these results, we construct enhanced structure maps for the binary alloys of platinum-group metals. These maps are much more complete, systematic, and predictive than those based on empirical results alone.

First-principles study of ferroelectricity and isotope effects in H-bonded KH2PO4 crystals

By means of extensive first-principles calculations we studied the ferroelectric phase transition and the associated isotope effect in KH2PO4 (KDP). Our calculations revealed that the spontaneous polarization of the ferroelectric phase is due to electronic charge redistributions and ionic displacements which are a consequence of proton ordering, and not vice versa. The experimentally observed double-peaked proton distribution in the paraelectric phase cannot be explained by a dynamics of only protons. This requires, instead, collective displacements within clusters that include also the heavier ions. These tunneling clusters can explain the recent evidence of tunneling obtained from Compton scattering measurements. The sole effect of mass change upon deuteration is not sufficient to explain the huge isotope effect. Instead, we find that structural modifications deeply connected with the chemistry of the H bonds produce a feedback effect on tunneling that strongly enhances the phenomenon. The resulting influence of the geometric changes on the isotope effect agrees with experimental data from neutron scattering. Calculations under pressure allowed us to analyze the issue of universality in the disappearance of ferroelectricity upon compression. Compressing DKDP so that the distance between the two peaks in the deuteron distribution is the same as for protons in KDP, corresponds to a modification of the underlying double-well potential, which becomes 23 meV shallower. This energy difference is what is required to modify the O-O distance in such a way as to have the same distribution for protons and deuterons. At the high pressures required experimentally, the above feedback mechanism is crucial to explain the magnitude of the geometrical effect.

Rational screening low-cost counter electrodes for dye-sensitized solar cells

Dye-sensitized solar cells have attracted intense research attention owing to their ease of fabrication, cost-effectiveness and high efficiency in converting solar energy. Noble platinum is generally used as catalytic counter electrode for redox mediators in electrolyte solution. Unfortunately, platinum is expensive and non-sustainable for long-term applications. Therefore, researchers are facing with the challenge of developing low-cost and earth-abundant alternatives. So far, rational screening of non-platinum counter electrodes has been hamstrung by the lack of understanding about the electrocatalytic process of redox mediators on various counter electrodes. Here, using first-principle quantum chemical calculations, we studied the electrocatalytic process of redox mediators and predicted electrocatalytic activity of potential semiconductor counter electrodes. On the basis of theoretical predictions, we successfully used rust (alpha-Fe2O3) as a new counter electrode catalyst, which demonstrates promising electrocatalytic activity towards triiodide reduction at a rate comparable to platinum.

High-resolution electron energy-loss spectroscopy of Ba Ti O 3∕ Sr Ti O 3 multilayers

The dielectric properties of BaTiO3 thin films and multilayers are different from bulk materials because of nanoscale dimensions, interfaces, and stress-strain conditions. In this study, BaTiO3/SrTiO3 multilayers deposited on SrTiO3 substrates by pulsed laser deposition have been investigated by high-energy-resolution electron energy-loss spectroscopy. The fine structures in the spectra are discussed in terms of crystal-field splitting and the internal strain. The crystal-field splitting of the BaTiO3 thin layer is found to be a little larger than that of bulk BaTiO3, which has been interpreted by the presence of the internal strain induced by the misfit at the interface. This finding is consistent with the lattice parameters of the BaTiO3 thin layer determined by the selected area diffraction pattern. The near-edge structure of the oxygen K edge in BaTiO3 thin layers and in bulk BaTiO3 are simulated by first-principle self-consistent full multiple-scattering calculations. The results of the simulations are in a good agreement with the experimental results. Moreover, the aggregation of oxygen vacancies at the rough BaTiO3/SrTiO3 interface is indicated by the increased [Ti]/[O] element ratio, which dominates the difference of dielectric properties between BaTiO3 layer and bulk materials.

On the anharmonic, multi-phonon, Debye-Waller contributions to the phonon-limited resistivity of metals : applications to Na and K

The anharmonic, multi-phonon (MP), and Oebye-Waller factor (OW) contributions to the phonon limited resistivity (;0) of metals derived by Shukla and Muller (1979) by the doubletime temperature dependent Green function method have been numerically evaluated for Na and K in the high temperature limit. The anharmonic contributions arise from the cubic and quartic shift of phonons (CS, QS), and phonon width (W) and the interference term (1). The QS, MP and OW contributions to I' are also derived by the matrix element method and the results are in agreement with those of Shukla and Muller (1979). In the high temperature limit, the contributions to;O from each of the above mentioned terms are of the type BT2 For numerical calculations suitable expressions are derived for the anharmonic contributions to ~ in terms of the third and fourth rank tensors obtained by the Ewald procedure. The numerical calculation of the contributions to;O from the OW, MP term and the QS have been done exactly and from the CS, Wand I terms only approximately in the partial and total Einstein approximations (PEA, TEA), using a first principle approach (Shukla and Taylor (1976)). The results obtained indicate that there is a strong pairwise cancellation between the: OW and MP terms, the QS and CS and the Wand I terms. The sum total of these contributions to;O for Na and K amounts to 4 to 11% and 2 to 7%, respectively, in the PEA while in the TEA they amount to 3 to 7% and 1 to 4%, respectively, in the temperature range.

The first real-time worldwide ionospheric predictions network: an advance in support of spaceborne experimentation, on-line model validation, and space weather

We report on the first realtime ionospheric predictions network and its capabilities to ingest a global database and forecast F-layer characteristics and "in situ" electron densities along the track of an orbiting spacecraft. A global network of ionosonde stations reported around-the-clock observations of F-region heights and densities, and an on-line library of models provided forecasting capabilities. Each model was tested against the incoming data; relative accuracies were intercompared to determine the best overall fit to the prevailing conditions; and the best-fit model was used to predict ionospheric conditions on an orbit-to-orbit basis for the 12-hour period following a twice-daily model test and validation procedure. It was found that the best-fit model often provided averaged (i.e., climatologically-based) accuracies better than 5% in predicting the heights and critical frequencies of the F-region peaks in the latitudinal domain of the TSS-1R flight path. There was a sharp contrast however, in model-measurement comparisons involving predictions of actual, unaveraged, along-track densities at the 295 km orbital altitude of TSS-1R In this case, extrema in the first-principle models varied by as much as an order of magnitude in density predictions, and the best-fit models were found to disagree with the "in situ" observations of Ne by as much as 140%. The discrepancies are interpreted as a manifestation of difficulties in accurately and self-consistently modeling the external controls of solar and magnetospheric inputs and the spatial and temporal variabilities in electric fields, thermospheric winds, plasmaspheric fluxes, and chemistry.

Electronic, vibrational and related properties of group IV metal oxides by ab initio calculations

We present our theoretical results for the structural, electronic, vibrational and optical properties of MO(2) (M = Sn, Zr, Hf and Ti) obtained by first-principles calculations. Relativistic effects are demonstrated to be important for a realistic description of the detailed structure of the electronic frequency-dependent dielectric function, as well as of the carrier effective masses. Based on our results, we found that the main contribution of the high values calculated for the oxides dielectric constants arises from the vibrational properties of these oxides, and the vibrational static dielectric constant values diminish with increasing pressure. (c) 2008 Elsevier B.V. All rights reserved.