Nitrogen defects and ferromagnetism in Cr-doped dilute magnetic semiconductor AlN from first principles

High Curie temperature of 900 K has been reported in Cr-doped AlN diluted magnetic semiconductors prepared by various methods, which is exciting for spintronic applications. It is believed that N defects play important roles in achieving the high-temperature ferromagnetism in good samples. Motivated by these experimental advances, we use a full-potential density-functional-theory method and supercell approach to investigate N defects and their effects on ferromagnetism of (Al,Cr)N with N vacancies (V-N). We investigate the structural and electronic properties of V-N, single Cr atom, Cr-Cr atom pairs, Cr-V-N pairs, and so on. In each case, the most stable structure is obtained by comparing different atomic configurations optimized in terms of the total energy and the force on every atom, and then it is used to calculate the defect formation energy and study the electronic structures. Our total-energy calculations show that the nearest substitutional Cr-Cr pair with the two spins in parallel is the most favorable and the nearest Cr-V-N pair makes a stable complex. Our formation energies indicate that V-N regions can be formed spontaneously under N-poor condition because the minimal V-N formation energy equals -0.23 eV or Cr-doped regions with high enough concentrations can be formed under N-rich condition because the Cr formation energy equals 0.04 eV, and hence real Cr-doped AlN samples are formed by forming some Cr-doped regions and separated V-N regions and through subsequent atomic relaxation during annealing. Both of the single Cr atom and the N vacancy create filled electronic states in the semiconductor gap of AlN. N vacancies enhance the ferromagnetism by adding mu(B) to the Cr moment each but reduce the ferromagnetic exchange constants between the spins in the nearest Cr-Cr pairs. These calculated results are in agreement with experimental observations and facts of real Cr-doped AlN samples and their synthesis. Our first-principles results are useful to elucidate the mechanism for the ferromagnetism and to explore high-performance Cr-doped AlN diluted magnetic semiconductors.

Stone-Wales defects created by low energy recoils in single-walled silicon carbide nanotubes

The defect creation at low energy events was studied using density functional theory molecular dynamics simulations in silicon carbide nanotubes, and the displacement threshold energies determined exhibit a dependence on sizes, which decrease with decreasing diameter of the nanotubes. The Stone-Wales (SW) defect, which is a common defect configurations induced through irradiation in nanotubes, has also been investigated, and the formation energies of the SW defects increase with increasing diameter of the nanotubes. The mean threshold energies were found to be 23 and 18 eV for Si and C in armchair (5,5) nanotubes. (C) 2009 American Institute of Physics. [doi: 10.1063/1.3238307]

Quantum transport from the perspective of quantum open systems

By viewing the non-equilibrium transport setup as a quantum open system, we propose a reduced-density-matrix based quantum transport formalism. At the level of self-consistent Born approximation, it can precisely account for the correlation between tunneling and the system internal many-body interaction, leading to certain novel behavior such as the non-equilibrium Kondo effect. It also opens a new way to construct time-dependent density functional theory for transport through large-scale complex systems. (c) 2006 Elsevier B.V. All rights reserved.

Mechanisms for extra conductance plateaus in quantum wires

Spin-density-functional theory is employed to calculate the conductance G through a quasi-one-dimensional quantum wire. In addition to the usual subband quantization plateaus at G=n(2e(2)/h), we find additional structures at (n+1/2)(2e(2)/h). The extra structures appear whether or not the electrons in the wire spin polarize. However, only the spin-polarized case reproduces the experimental temperature and magnetic field dependences.

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.

Phase equilibria and plate-fluid interfacial tensions for associating hard sphere fluids confined in slit pores

The excess Helmholtz free energy functional for associating hard sphere fluid is formulated by using a modified fundamental measure theory [Y. X. Yu and J. Z. Wu, J. Chem. Phys. 117, 10156 (2002)]. Within the framework of density functional theory, the thermodynamic properties including phase equilibria for both molecules and monomers, equilibrium plate-fluid interfacial tensions and isotherms of excess adsorption, average molecule density, average monomer density, and plate-fluid interfacial tension for four-site associating hard sphere fluids confined in slit pores are investigated. The phase equilibria inside the hard slit pores and attractive slit pores are determined according to the requirement that temperature, chemical potential, and grand potential in coexistence phases should be equal and the plate-fluid interfacial tensions at equilibrium states are predicted consequently. The influences of association energy, fluid-solid interaction, and pore width on phase equilibria and equilibrium plate-fluid interfacial tensions are discussed.

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.

1ST-PRINCIPLES CALCULATIONS FOR QUASI-PARTICLE ENERGIES OF GAP AND GAAS

We have applied the Green-function method in the GW approximation to calculate quasiparticle energies for the semiconductors GaP and GaAs. Good agreement between the calculated excitation energies and the experimental results was achieved. We obtained calculated direct band gaps of GaP and GaAs of 2.93 and 1.42 eV, respectively, in comparison with the experimental values of 2.90 and 1.52 eV, respectively. An ab initio pseudopotential method has been used to generate basis wave functions and charge densities for calculating the dielectric matrix elements and self-enegies. To evaluate the dynamical effects of the screened interaction, the generalized-plasma-pole model has been utilized to extend the dielectric matrix elements from static results to finite frequencies. We presen the calculated quasiparticle energies at various high-symmetry points of the Brillouin zone and compare them with the experimental results and other calculations.

1ST PRINCIPLE CALCULATIONS OF QUASI-PARTICLE ENERGIES FOR BAND STRUCTURES OF SEMICONDUCTORS

We successfully applied the Green function theory in GW approximation to calculate the quasiparticle energies for semiconductors Si and GaAs. Ab initio pseudopotential method was adopted to generate basis wavefunctions and charge densities for calculating dielectric matrix elements and electron self-energies. To evaluate dynamical effects of screened interaction, GPP model was utilized to extend dieletric matrix elements from static results to finite frequencies. We give a full account of the theoretical background and the technical details for the first principle pseudopotential calculations of quasiparticle energies in semiconductors and insulators. Careful analyses are given for the effective and accurate evaluations of dielectric matrix elements and quasiparticle self-energies by using the symmetry properties of basis wavefunctions and eigenenergies. Good agreements between the calculated excitation energies and fundamental energy gaps and the experimental band structures were achieved.

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.

First principles based predictions of the toughness of a metal/oxide interface

We describe a first-principles-based strategy to predict the macroscopic toughness of a gamma-Ni(Al)/alpha-Al2O3 interface. Density functional theory calculations are used to ascertain energy changes upon displacing the two materials adjacent to the interface, with relaxation conducted over all atoms located within adjoining rows. Traction/displacernent curves are obtained from derivatives of the energy. Calculations are performed in mode I (opening), mode II (shear) and at a phase angle of 45 degrees. The shear calculations are conducted for displacements along < 110 > and < 112 > of the Ni lattice. A generalized interface potential function is used to characterize the results. Initial fitting to both the shear and normal stress results is required to calibrate the unknowns. Thereafter, consistency is established by using the potential to predict other traction quantities. The potential is incorporated as a traction/displacement function within a cohesive zone model and used to predict the steady-state toughness of the interface. For this purpose, the plasticity of the Ni alloy must be known, including the plasticity length scale. Measurements obtained for a gamma-Ni superalloy are used and the toughness predicted over the full range of mode mixity. Additional results for a range of alloys are used to demonstrate the influences of yield strength and length scale.

First principles based predictions of the toughness of a metal/oxide interface

We describe a first-principles-based strategy to predict the macroscopic toughness of a gamma-Ni(Al)/alpha-Al2O3 interface. Density functional theory calculations are used to ascertain energy changes upon displacing the two materials adjacent to the interface, with relaxation conducted over all atoms located within adjoining rows. Traction/displacernent curves are obtained from derivatives of the energy. Calculations are performed in mode I (opening), mode II (shear) and at a phase angle of 45 degrees. The shear calculations are conducted for displacements along < 110 > and < 112 > of the Ni lattice. A generalized interface potential function is used to characterize the results. Initial fitting to both the shear and normal stress results is required to calibrate the unknowns. Thereafter, consistency is established by using the potential to predict other traction quantities. The potential is incorporated as a traction/displacement function within a cohesive zone model and used to predict the steady-state toughness of the interface. For this purpose, the plasticity of the Ni alloy must be known, including the plasticity length scale. Measurements obtained for a gamma-Ni superalloy are used and the toughness predicted over the full range of mode mixity. Additional results for a range of alloys are used to demonstrate the influences of yield strength and length scale.