No Arabic abstract
We present calculations of the magnetic ground states of Cr trimers in different geometries on top of a Au(111) surface. By using a least square fit method based on a fully relativistic embedded-cluster Greens function method first we determined the parameters of a classical vector-spin model consisting of second and fourth order interactions. The newly developed method requires no symmetry constraints, therefore, it is throughout applicable for small nanoparticles of arbitrary geometry. The magnetic ground states were then found by solving the Landau-Lifshitz-Gilbert equations. In all considered cases the configurational energy of the Cr trimers is dominated by large antiferromagnetic nearest neighbor interactions, whilst biquadratic spin-interactions have the second largest contributions to the energy. We find that an equilateral Cr trimer exhibits a frustrated 120$^circ$ Neel type of ground state with a small out-of-plane component of the magnetization and we show that the Dzyaloshinsky-Moriya interactions determine the chirality of the magnetic ground state. In cases of a linear chain and an isosceles trimer collinear antiferromagnetic ground states are obtained with a magnetization lying parallel to the surface.
The ground state band structure, magnetic moments, charges and population numbers of electronic shells of Cu and Fe atoms have been calculated for chalcopyrite CuFeS2 using density functional theory. The comparison between our calculation results and experimental data (X-ray photoemission, X-ray absorption and neutron diffraction spectroscopy) has been made. Our calculations predict a formal oxidation state for chalcopyrite as Cu$^{1+}$Fe$^{3+}$S$_{2}$$^{2-}$. However, the assignment of formal valence state to transition metal atoms appears to be oversimplified. It is anticipated that the valence state can be confirmed experimentally by nuclear magnetic and nuclear quadrupole resonance and Mossbauer spectroscopy methods.
We investigate equilibrium and transport properties of a copper phthalocyanine (CuPc) molecule adsorbed on Au(111) and Ag(111) surfaces. The CuPc molecule has essentially three localized orbitals close to the Fermi energy resulting in strong local Coulomb repulsion not accounted for properly in density functional calculations. Hence, they require a proper many-body treatment within, e.g., the Anderson impurity model (AIM). The occupancy of these orbitals varies with the substrate on which CuPc is adsorbed. Starting from density functional theory calculations, we determine the parameters for the AIM embedded in a noninteracting environment that describes the residual orbitals of the entire system. While correlation effects in CuPc on Au(111) are already properly described by a single orbital AIM, for CuPc on Ag(111) the three orbital AIM problem can be simplified into a two orbital problem coupled to the localized spin of the third orbital. This results in a Kondo effect with a mixed character, displaying a symmetry between SU(2) and SU(4). The computed Kondo temperature is in good agreement with experimental values. To solve the impurity problem we use the recently developed fork tensor product state solver. To obtain transport properties, a scanning tunneling microscope (STM) tip is added to the CuPc molecule absorbed on the surface. We find that the transmission depends on the detailed position of the STM tip above the CuPc molecule in good agreement with differential conductance measurements.
In a recent preprint Kong et al, arXiv:0902.0642v1 (2009) claimed to calculate the lattice thermal conductivity of single and bi-layer graphene from first principles. The main findings were that the Umklapp-limited thermal conductivity is only slightly higher than that of high-quality bulk graphite along the basal plane, and that it does not strongly depend on the number of atomic layers. Here we explain that the calculation of Kong et al used a truncation procedure with a hidden parameter, a cut-off frequency for the long-wavelength acoustic phonons, which essentially determined the final result. Unlike in bulk graphite, there is no physical justification for introducing the cut-off frequency for the long wavelength phonons in graphene. It leads to substantial underestimation of graphenes lattice thermal conductivity and a wrong conclusion about the dependence on the number of atomic layers. We outline the proper way for calculating the lattice thermal conductivity of graphene, which requires an introduction of other scattering mechanisms to avoid a logarithmic divergence of the thermal conductivity integral.
We present an ab initio study of the interface energies of cubic yttria-stabilized zirconia (YSZ) epitaxial layers on a (0001)_{alpha-Al_2O_3} substrate. The interfaces are modelled using a supercell geometry and the calculations are carried out in the framework of density-functional theory (DFT) and the local-density approximation (LDA) using the projector-augmented-wave (PAW) pseudopotential approach. Our calculations clearly demonstrate that the (111)_{YSZ} || (0001)_{alpha-Al_2O_3} interface energy is lower than that of (100)_{YSZ} || (0001)_{alpha-Al_2O_3}. This result is central to understanding the behaviour of YSZ thin solid film islanding on (0001)_{alpha-Al_2O_3} substrates, either flat or in presence of defects.
We present a systematic study of the atomic and electronic structure of the Si(111)-(5x2)-Au reconstruction using first-principles electronic structure calculations based on the density functional theory. We analyze the structural models proposed by Marks and Plass [Phys. Rev. Lett.75, 2172 (1995)], those proposed recently by Erwin [Phys. Rev. Lett.91, 206101 (2003)], and a completely new structure that was found during our structural optimizations. We study in detail the energetics and the structural and electronic properties of the different models. For the two most stable models, we also calculate the change in the surface energy as a function of the content of silicon adatoms for a realistic range of concentrations. Our new model is the energetically most favorable in the range of low adatom concentrations, while Erwins 5x2 model becomes favorable for larger adatom concentrations. The crossing between the surface energies of both structures is found close to 1/2 adatoms per 5x2 unit cell, i.e. near the maximum adatom coverage observed in the experiments. Both models, the new structure and Erwins 5x2 model, seem to provide a good description of many of the available experimental data, particularly of the angle-resolved photoemission measurements.