Chemical trends of defect formation in Si quantum dots: The case of group-III and group-V dopants

Using first-principles methods, we have systematically calculated the defect formation energies and transition energy levels of group-III and group-V impurities doped in H passivated Si quantum dots (QDs) as functions of the QD size. The general chemical trends found in the QDs are similar to that found in bulk Si. We show that defect formation energy and transition energy level increase when the size of the QD decreases; thus, doping in small Si QDs becomes more difficult. B-Si has the lowest acceptor transition energy level, and it is more stable near the surface than at the center of the H passivated Si QD. On the other hand, P-Si has the smallest donor ionization energy, and it prefers to stay at the interior of the H passivated Si QD. We explained the general chemical trends and the dependence on the QD size in terms of the atomic chemical potentials and quantum confinement effects.

Quantum master equation scheme of time-dependent density functional theory to time-dependent transport in nanoelectronic devices

In this work a practical scheme is developed for the first-principles study of time-dependent quantum transport. The basic idea is to combine the transport master equation with the well-known time-dependent density functional theory. The key ingredients of this paper include (i) the partitioning-free initial condition and the consideration of the time-dependent bias voltages which base our treatment on the Runge-Gross existence theorem; (ii) the non-Markovian master equation for the reduced (many-body) central system (i.e., the device); and (iii) the construction of Kohn-Sham master equations for the reduced single-particle density matrix, where a number of auxiliary functions are introduced and their equations of motion (EOMs) are established based on the technique of spectral decomposition. As a result, starting with a well-defined initial state, the time-dependent transport current can be calculated simultaneously along with the propagation of the Kohn-Sham master equation and the EOMs of the auxiliary functions.

Electronic, structural, and optical properties of crystalline yttria

The electronic structure of crystalline Y2O3 is investigated by first-principles calculations within the local-density approximation (LDA) of the density-functional theory. Results are presented for the band structure, the total density of states (DOS), the atom-and orbital-resolved partial DOS. effective charges, bond order, and charge-density distributions. Partial covalent character in the Y-O bonding is shown, and the nonequivalency of the two Y sites is demonstrated. The calculated electronic structure is compared with a variety of available experimental data. The total energy of the crystal is calculated as a function of crystal volume. A bulk modulus B of 183 Gpa and a pressure coefficient B' of 4.01 are obtained, which are in good agreement with compression data. An LDA band gap of 4.54 eV at Gamma is obtained which increases with pressure at a rate of dE(g)/dP = 0.012 eV/Gpa at the equilibrium volume. Also investigated are the optical properties of Y2O3 up to a photon energy of 20 eV. The calculated complex dielectric function and electron-energy-loss function are in good agreement with experimental data. A static dielectric constant of epsilon(O)= 3.20 is obtained. It is also found that the bottom of the conduction band consists of a single band, and direct optical transition at Gamma between the top of the valence band and the bottom of the conduction band may be symmetry forbidden.

无机材料的结构与物理性能：第一性原理计算

本论文的研究方向是通过密度泛函理论研究各种材料的物理化学性质。近年来相关理论和数值算法的飞速发展，使得基于密度泛函理论的第一性原理方法成为凝聚态物理、量子化学和材料科学中的常规计算研究手段。本论文对无机材料进行第一性原理研究，研究涉及材料物性包括几何构型、电子结构、磁性能、和力学性质等等。 第一章简要地介绍了密度泛函理论的基本框架和近年来的理论发展。密度泛函理论的发展以及寻找合适的交换相关能量泛函为主线。从最初的局域密度近似（LDA）、广义梯度近似（GGA）到现在的非局域泛函、自相互作用修正，多种泛函形式的出现使得密度泛函理论可以提供越来越精确的计算结果。除了改进交换相关泛函，近年来密度泛函理论向动力学平均场和含时理论等方面扩展也很活跃。这些扩展式的密度泛函理论的应用领域不断扩大。在本章的最后，我们介绍一些密度泛函理论的应用程序。 第二章我们通过第一性原理从头算系统的研究了5d过渡金属二硼化物TMB2 (TM = Hf, Ta, W, Re, Os and Ir, Pt)在假想的Pmmn空间群中的结合能、生成焓、力学常数和电子能带结构。我们的计算结果表明在力学常数和价电子浓度之间存在一种关系：即当价电子浓度在6.8到7.2之间时，体模量和剪切模量达到最大值。再者，这种力学常数与价电子浓度之间的关系可以通过电子能带结构分析，如费米面附近价电子的占据情况，决定了体系的结合能和力学性能。最大的体模量和剪切模量的获得归因于TM d-B p成键态的几乎完全占据和反键态的未占据。依据上述这种关系，我们预测了在正交结构OsB2合金化W和Re将比合金化Ir元素更容易和更硬。事实上，我们的进一步计算证实了我们的期望。 通过第一性原理计算研究了ReB2和WB2的结构、弹性力学，和电子性质。计算结果表明：ReB2的平衡态结构参数和报导的实验结构一致。在常温常压下，WB2在P63/mmc空间群比在P6/mmm更稳定。依据我们计算的多晶聚集体的体模量、剪切模量，ReB2和WB2可以被看作是具有前景的低压缩率和硬材料。再者，化合物的力学各项异性通过计算得到的弹性力学常数来详细的分析讨论。态密度和电子密度分析揭示过渡金属和硼原子之间的共价键是材料具有高的体模量、剪切模量和小的Poisson比率的原因。 第四章，采用WIEN2k程序包中的缀加平面波+局域轨道方法的第一性原理计算，研究了层状钙钛矿化合物Cs2AgF4的结构、电子结构和磁性等性质。我们的计算结果表明Cs2AgF4的基态为正交相，能量比四方相低。我们同时发现Cs2AgF4应该表现出一种强的二维铁磁性，铁磁层之间为弱的反铁磁性耦合，这一结果与实验观察相一致。更主要的是，通过分析态密度图和自旋电子密度等密度面的分析，可以清楚地看出体系中存在 和 轨道的反铁电弹性有序现象。 第五章中，采用WIEN2k程序包中的缀加平面波+局域轨道方法的第一性原理计算，研究了YBaFe2O5的晶体结构、电子和磁性等方面的物理性能，特别是有关电荷和轨道有序的情况。尽管总的3d电荷不均衡程度很小，以Fe2+ 和Fe3+阳离子t2g轨道占据数的差别所定义的轨道有序序参量相当大（0.73），无可非议的显示YBaFe2O5中存在电荷和轨道有序。O 2p和Fe eg轨道之间的强的杂化作用使得Fe2+ 和Fe3+阳离子之间总电荷差别几乎完全消失。此外，我们讨论了轨道有序和电荷有序以及磁有序之间的关系。dxz轨道有序决定了G型反铁磁性自旋有序的稳定性和电荷有序模式。

Room-Temperature Ferromagnetism in Co-Doped In2O3 Nanocrystals

Co-doped In2O3 nanocrystals showing room-temperature ferromagnetism have been successfully prepared by a simple sol-gel synthesis route. The sample displays it clear ferromagnetism behavior above 300 K. Phase and structure analyses reveal that the nanocrystals are crystallized with Co ions substituted for In ions in the In2O3 matrix, and no trace of secondary phases or clusters is detected. The experimental results are explained theoretically by first-principles calculations based on density functional theory, which indicate that the native ferromagnetic behavior of Co-doped In2O3 could be mainly ascribed to the strong d-d coupling of the magnetic ions.

Anomalous optical and electronic properties of dense sodium

Based on the density functional theory, we systematically study the optical and electronic properties of the insulating dense sodium phase (Na-hp4) reported recently (Ma et al., 2009). The structure is found optically anisotropic. Through Bader analysis, we conclude that ionicity exists in the structure and becomes stronger with increasing pressure.

Charge Separation in Wurtzite/Zinc-Blende Heterojunction GaN Nanowires

The electronic properties of wurtzite/zinc-blende (WZ/ZB) heterojunction GaN are investigated using first-principles methods. A small component of ZB stacking formed along the growth direction in the WZ GaN nanowires does not show a significant effect on the electronic property, whereas a charge separation of electrons and holes occurs along the directions perpendicular to the growth direction in the ZB stacking. The later case provides an efficient way to separate the charge through controlling crystal structure. These results have significant implications for most state of the art excitonic solar cells and the tuning region in tunable laser diodes.

ELECTRONIC AND MAGNETIC-STRUCTURE OF THE TERNARY COMPOUND ND2FE17N

The electronic and the magnetic structure of the Nd2Fe17N1 phase in the family of Nd-Fe-N ternary compounds have been calculated using the first-principles, spin-polarized orthogonalized linear-combination-of-atomic-orbitals method. Results are presented in the form of site-decomposed and spin-projected partial density of states. The occupation sites of the three N atoms are determined by an average radial distribution of all possible N site configurations. Both cases of N occupying the 3b and the 18g sites are studied. The results indicate that the 6c Fe sites have the maximum and the 18h Fe sites have the minimum local moments. This is in good agreement with experiment. It is concluded that the influence on the local moment due to lattice expansion is less important compared to that due to interatomic interaction between the N atom and its neighbors. The results also show the important role of N atoms in raising the Curie temperature of this compound.

Ab initio calculations for the neutral and charged O vacancy in sapphire

The energetics, lattice relaxation, and the defect-induced states of st single O vacancy in alpha-Al2O3 are studied by means of supercell total-energy calculations using a first-principles method based on density-functional theory. The supercell model with 120 atoms in a hexagonal lattice is sufficiently large to give realistic results for an isolated single vacancy (square). Self-consistent calculations are performed for each assumed configuration of lattice relaxation involving the nearest-neighbor Al atoms and the next-nearest-neighbor O atoms of the vacancy site. Total-energy data thus accumulated are used to construct an energy hypersurface. A theoretical zero-temperature vacancy formation energy of 5.83 eV is obtained. Our results show a large relaxation of Al (O) atoms away from the vacancy site by about 16% (8%) of the original Al-square (O-square) distances. The relaxation of the neighboring Al atoms has a much weaker energy dependence than the O atoms. The O vacancy introduces a deep and doubly occupied defect level, or an F center in the gap, and three unoccupied defect levels near the conduction band edge, the positions of the latter are sensitive to the degree of relaxation. The defect state wave functions are found to be not so localized, but extend up to the boundary of the supercell. Defect-induced levels are also found in the valence-band region below the O 2s and the O 2p bands. Also investigated is the case of a singly occupied defect level (an F+ center). This is done by reducing both the total number of electrons in the supercell and the background positive charge by one electron in the self-consistent electronic structure calculations. The optical transitions between the occupied and excited states of the: F and F+ centers are also investigated and found to be anisotropic in agreement with optical data.

Electronic structure of microporous titanosilicate ETS-10

The electronic structure of a microporous titanosilicate framework, ETS-10 is calculated by means of a first-principles self-consistent method. It is shown that without the inclusion of the alkali atoms whose positions in the framework are unknown, ETS-10 is an electron deficient system with 32 electrons per unit cell missing at the top of an otherwise semiconductor-like band structure. The calculated density of slates are resolved into partial components. It is shown that the states of the missing electrons primarily originate from the Ti-O bond. The local density of states of the Ti-3d orbitals in the ETS-10 framework is quite different from the perovskite BaTiO3. The possibilities of ETS-10 crystal being ferroelectric or having other interesting properties are discussed.

Ab initio study on the physical properties of CoN3 and RhN3 with skutterudite structure

The elastic and electronic properties of hypothetical CoN3 and RhN3 with cubic skutterudite structure were studied by first principles calculations based on density functional theory. By choosing different initial geometries, two local minima or modifications were located on the potential energy surface, termed as modifications I and II. Both compounds are mechanically stable. For each compound, modification I is lower in energy than II. Thermodynamically stable phases can be achieved by applying pressures. Modification II is lower in energy than I at above 50 GPa for both compounds.

Structural stability and electronic properties of Co2N, Rh2N and Ir2N

The structural stability and electronic properties of Co2N, Rh2N and Ir2N were Studied by using the first principles based on the density functional theory. Two Structures were considered for each nitride, orthorhombic Pnnm phase and cubic Pa (3) over bar phase. The results show that they are all mechanically stable. Co2N in both phases are thermodynamically stable due to the negative formation energy, while the remaining two compounds are thermodynamically unstable.

Charge disproportionation in CaCu3Fe4O12

First-principles calculations using the APW+lo method, as implemented in the WIEN2K code, have been used to investigate the structural, electronic, and magnetic properties of the perovskite CaCu3Fe4O12, including the high-temperature Im-3 and low-temperature Pn-3 phase. The high-temperature phase presents a homogeneous valence and an orbital degenerate half-metallic behavior, which is consistent with the previous theoretical result. Instead orbital ordering, charge ordering, or disproportionation on Fe sites occur in the low-temperature phase, leading to the insulating character. More importantly, the charge disproportionation is of 2d(5)L -> d(5)L(2) + d(5) type (where L denotes an oxygen hole or a ligand hole), and the origin for the phenomenon is discussed in detail.

通过第一性原理计算研究内地核的结构与成分

The determination of the composition and structure of the Earth’s inner core has long been the major subject in the study of the Earth’s deep interior. It’s widely believed that the Earth’s core is formed by iron with a fraction of nickel. However, light elements must exist in the inner core because the earth core is less dense than pure iron-nickel alloy (~2-3% in the solid inner core and ~6-7% in the liquid outer core). The questions are what and how much light element is there in the iron-nickel alloy. Besides the composition, the crystal structure of the iron with or without light element is also not well known. According to the seismological observations, the sound waves propagate 3-4% faster along the spin axis than in the equatorial plane. That means the inner core is anisotropic. The densest structure of iron-nickel alloy should be h.c.p structure under the very high pressures. However, the h,c,p structure does not propagate waves anisotropic ally. Then what is the structure of the iron-nickel alloy or the iron-nickle-light element alloy. In this study, we tried to predict the composition and the structure of the inner core through ab initio calculation of the Gibbs free energy, which is a function of internal energy, density and entropy. We conclude that the h.c.p structure is more stable than the b.c.c structure under high pressure and 0 K, but with the increase of temperature, the free energy of the b.c.c structure is decreasing much faster than the h.c.p structure caused by the vibration of the atomics, so the b.c.c structure is more stable at high temperatures. With the addition of light elements (S or Si or both), the free energy of b.c.c. decreases even faster, about 3at% of Si not only explains why the inner core is about 2-3 % lighter than the iron-nickle alloy, but also reasons why the inner core is anisotropic, since the b.c.c. structure becomes more stable than the h.c.p structure at 5500-6000K and b.c.c. is anisotropic in propagating seismic waves. Therefore, we infer that the inner core of the earth is formed by b.c.c iron and a fraction of nickel plus ~3at.% Si, with a temperature higher than 5500K, which is consistent with the studies from other approaches.

分子水平上研究地质流体的物理化学性质

Geological fluids exist in every geosphere of the Earth and play important roles in many processes of material transformations, energetic interchanges and geochemical interactions. To study the physicochemical properties and geochemical behaviors of geological fluids turn Girt to be one of the challenging issues in geosciences. Compared with conventional approaches of experiments and semi-theoretical modeling, computer simulation on molecular level shows its advantages on quantitative predictions of the physicochemical properties of geological fluids under extreme conditions and emerges as a promising approach to find the characteristics of geological fluids and their interactions in different geospheres of the Earth interior.This dissertation systematically discusses the physicochemical properties of typical geological fluids with state-of-the-art computer simulation techniques. The main results can be summarized as follows: (1) The experimental phase behaviors of the systems CH4-C2H6 and. CO2 have been successfully reproduced with Monte Carlo simulations. (2) Through comprehensive isothermal-isobaric molecular dynamics simulations, the PVT data of water hia^e been extended beyond experimental range to about 2000 K and 20 GPa and an improved equation of state for water has been established. (3) Based on extensive computer simulations, am optimized molecular potential for carbon dioxide have been proposed, this model is expected to predict different properties of carbon dioxide (volumetric properties, phase equilibria, heat of vaporization, structural and dynamical properties) with improved accuracies. (4) On the basis of the above researches of the end-members, a set of parameters for unlike interactions has been proposed by non-linear fitting to the ab initio potential surface of CO2-H2O and is superior to the common used mixing rule and the results of prior workers vs/Ith remarkable accuracies, then a number of simulations of the mixture have been carried out to generate data under high temperatures and pressures as an important complement to the limited experiments. (5) With molecular dynamics simulations, various structural, dynamical and thermodynamical properties of ionic solvations and associations have been oomprehensively analyzed, these results not only agree well with experimental data and first principle calculation results, but also reveal some new insights into the microscopic ionic solvation and association processes.