Search

Two of the most widely used electronic structure theory methods, namely Hartree-Fock and Kohn-Sham density functional theory, both requires the iterative solution of a set of Schr\"odinger-like equations. The speed of convergence of such self-consistent field process depends on the complexity of the system under investigation and on the initial guess for the density matrix. An initial density matrix close to the ground-state one will effectively allow one to cut out many of the self-consistent steps necessary to achieve convergence. Here, we predict the density matrix of Kohn-Sham density functional theory by constructing a neural network, which uses the atomic positions as only information. Such neural network provides an initial guess for the density matrix far superior to any other recipes available, in particular for molecules with metallic bonds. Furthermore, the quality of such neural-network density matrix is good enough for the evaluation of interatomic forces. This allows us to run accelerated {\it ab-initio} molecular dynamics with little to no self-consistent steps.

\noindent Using results drawn from the pages of PRL\cite{Zhu2014}, the authors comment on \textquotedblleft Semiconducting Layered Blue Phosphorus: A Computational Study\textquotedblright. In recent letter\cite{Zhu2014} unknown phase of phosphorus with high stability and a wide fundamental gap was proposed by Zhu, and Tom\'anek. The first half of the comment questioned elements of phosphorus blue structure, while the second half devoted to van der Waals forces between AB stacked blue phosphorus layers which missed by them., Comment: 1 pages, 1 figures

We present a computational scheme for extracting the energy level alignment of a metal/molecule interface, based on constrained density functional theory and local exchange and correlation functionals. The method, applied here to benzene on Li(100), allows us to evaluate charge transfer energies, as well as the spatial distribution of the image charge induced on the metal surface. We systematically study the energies for charge transfer from the molecule to the substrate as function of the molecule-substrate distance, and investigate the effects arising from image charge confinement and local charge neutrality violation. For benzene on Li(100) we find that the image charge plane is located at about 1.8 \AA above the Li surface, and that our calculated charge transfer energies compare perfectly with those obtained with a classical electrostatic model having the image plane located at the same position. The methodology outlined here can be applied to study any metal/organic interface in the weak coupling limit at the computational cost of a total energy calculation. Most importantly, as the scheme is based on total energies and not on correcting the Kohn-Sham quasiparticle spectrum, accurate results can be obtained with local/semi-local exchange and correlation functionals. This enables a systematic approach to convergence.

In transition metal doped ZnO, the energy position of dopant 3$d$ states relative to host conduction and valence bands is crucial in determining the possibilty of long range ferromagnetism. Density functional theory based estimates of the energy position of Co-3$d$ states in Co doped ZnO differ substantially depending upon the choice of exchange-correlation functional. In this work we investigate many-body $GW$ corrections on top of DFT$+U$ and hybrid-DFT groundstates to provide a theoretical benchmark for the quasiparticle energies in wurtzite ZnO:Co. Both single shot $G_0W_0$ as well as partially self-consistent $GW_0$ wherein the wavefunctions are held fixed at the DFT level but the eigenvalues in G are iterated, are considered. The predicted energy position of the minority spin Co-$t_2$ states is 3.0-3.6 eV above the ZnO conduction band minimum which is closer to hybrid-DFT based estimates., Comment: 6 pages, 2 figures

Octahedral Fe$^{2+}$ molecules are particularly interesting as they often exhibit a spin-crossover transition. In spite of the many efforts aimed at assessing the performances of density functional theory for such systems, an exchange-correlation functional able to account accurately for the energetic of the various possible spin-states has not been identified yet. Here we critically discuss the issues related to the theoretical description of this class of molecules from first principles. In particular we present a comparison between different density functionals for four ions, namely [Fe(H$_2$O)$_6$]$^{2+}$, [Fe(NH$_3$)$_6$]$^{2+}$, [Fe(NCH)$_6$]$^{2+}$ and [Fe(CO)$_6$]$^{2+}$. These are characterized by different ligand-field splittings and ground state spin multiplicities. Since no experimental data are available for the gas phase, the density functional theory results are benchmarked against those obtained with diffusion Monte Carlo, one of the most accurate methods available to compute ground state total energies of quantum systems. On the one hand, we show that most of the functionals considered provide a good description of the geometry and of the shape of the potential energy surfaces. On the other hand, the same functionals fail badly in predicting the energy differences between the various spin states. In the case of [Fe(H$_2$O)$_6$]$^{2+}$, [Fe(NH$_3$)$_6$]$^{2+}$, [Fe(NCH)$_6$]$^{2+}$, this failure is related to the drastic underestimation of the exchange energy. Therefore quite accurate results can be achieved with hybrid functionals including about 50% of Hartree-Fock exchange. In contrast, in the case of [Fe(CO)$_6$]$^{2+}$, the failure is likely to be caused by the multiconfigurational character of the ground state wave-function and no suitable exchange and correlation functional has been identified.

We present a time-dependent study of electron transport through a strongly correlated quantum dot. The time-dependent current is obtained with the multiple-probe battery method, while adiabatic lattice density functional theory in the Bethe ansatz local-density approximation to the Hubbard model describes the dot electronic structure. We show that for a certain range of voltages the quantum dot can be driven into a dynamical state characterized by regular current oscillations. This is a manifestation of a recently proposed dynamical picture of Coulomb blockade. Furthermore, we investigate how the various approximations to the electron-electron interaction affect the line-shapes of the Coulomb peaks and the I-V characteristics. We show that the presence of the derivative discontinuity in the approximate exchange-correlation potential leads to significantly different results compared to those obtained at the simpler Hartree level of description. In particular, a negative differential conductance (NDC) in the I-V characteristics is observed at large bias voltages and large Coulomb interaction strengths. We demonstrate that such NDC originates from the combined effect of electron-electron interaction in the dot and the finite bandwidth of the electrodes., Comment: 10 pages, 7 figures

The electronic origin of a large resistance change in nanoscale junctions incorporating spin crossover molecules is demonstrated theoretically by using a combination of density functional theory and the non-equilibrium Green's functions method for quantum transport. At the spin crossover phase transition there is a drastic change in the electronic gap between the frontier molecular orbitals. As a consequence, when the molecule is incorporated in a two terminal device, the current increases by up to four orders of magnitude in response to the spin change. This is equivalent to a magnetoresistance effect in excess of 3,000 %. Since the typical phase transition critical temperature for spin crossover compounds can be extended to well above room temperature, spin crossover molecules appear as the ideal candidate for implementing spin devices at the molecular level.

Rapid and cost-effective DNA sequencing at the single nucleotide level might be achieved by measuring a transverse electronic current as single-stranded DNA is pulled through a nano-sized pore. In order to enhance the electronic coupling between the nucleotides and the electrodes and hence the current signals, we employ a pair of single-walled close-ended (6,6) carbon nanotubes (CNTs) as electrodes. We then investigate the electron transport properties of nucleotides sandwiched between such electrodes by using first-principles quantum transport theory. In particular we consider the extreme case where the separation between the electrodes is the smallest possible that still allows the DNA translocation. The benzene-like ring at the end cap of the CNT can strongly couple with the nucleobases and therefore both reduce conformational fluctuations and significantly improve the conductance. The optimal molecular configurations, at which the nucleotides strongly couple to the CNTs, and which yield the largest transmission, are first identified. Then the electronic structures and the electron transport of these optimal configurations are analyzed. The typical tunneling currents are of the order of 50 nA for voltages up to 1 V. At higher bias, where resonant transport through the molecular states is possible, the current is of the order of several $\mu$A. Below 1 V the currents associated to the different nucleotides are consistently distinguishable, with adenine having the largest current, guanine the second-largest, cytosine the third and finally thymine the smallest. We further calculate the transmission coefficient profiles as the nucleotides are dragged along the DNA translocation path and investigate the effects of configurational variations. Based on these results we propose a DNA sequencing protocol combining three possible data analysis strategies., Comment: 12 pages, 17 figures, 3 tables

We present a fully variational generalization of the pseudo self-interaction correction (VPSIC) approach previously presented in two implementations based on plane-waves and atomic orbital basis set, known as PSIC and ASIC, respectively. The new method is essentially equivalent to the previous version for what concern the electronic properties, but it can be exploited to calculate total-energy derived properties as well, such as forces and structural optimization. We apply the method to a variety of test cases including both non-magnetic and magnetic correlated oxides and molecules, showing a generally good accuracy in the description of both structural and electronic properties., Comment: 23 pages, 9 tables, 16 figures

Here we show the room temperature integrated detection of the helicity of photons with 1300 nm wavelength, via spin- photodiodes based on fully epitaxial Fe/MgO/Ge(001) heterostructures. These devices convert the photon helicity in a variation of the photocurrent ({\Delta}I) for opposite orientations of the Fe magnetization. We found {\Delta}I on the order of 5.9% for 100% circularly polarized light, both in the forward and the reverse bias, with a contribution of 0.8% from the Fe dichroic light absorption and of 5.1% from the spin filtering of photo-generated carriers crossing the Fe/MgO/Ge tunneling barrier., Comment: This paper has been withdrawn by the author due to some imperfections in the figures

The magnetic properties of the transition metal monoxides MnO and NiO are investigated at equilibrium and under pressure via several advanced first-principles methods coupled with Heisenberg Hamiltonian MonteCarlo. The comparative first-principles analysis involves two promising beyond-local density functionals approaches, namely the hybrid density functional theory and the recently developed variational pseudo-self-interaction correction method, implemented with both plane-wave and atomic-orbital basis sets. The advanced functionals deliver a very satisfying rendition, curing the main drawbacks of the local functionals and improving over many other previous theoretical predictions. Furthermore, and most importantly, they convincingly demonstrate a degree of internal consistency, despite differences emerging due to methodological details (e.g. plane waves vs. atomic orbitals)

The evolution of the magnetic ordering temperature of the 4d3 perovskites RTcO3 (R=Ca,Sr,Ba) and its relation with its electronic and structural properties has been studied by means of hybrid density functional theory and Monte Carlo simulations. When compared to the most widely studied 3d perovskites the large spatial extent of the 4d shells and their relatively strong hybridization with oxygen weaken the tendency to form Jahn-Teller like orbital ordering. This strengthens the superexchange interaction. The resulting insulating G-type antiferromagnetic ground state is characterized by large superexchange coupling constants (26-35 meV) and Neel temperatures (750-1200 K). These monotonically increase as a function of the R ionic radius due to the progressive enhancement of the volume and the associated decrease of the cooperative rotation of the TcO6 octahedra., Comment: 4 pages, 3 figures

We investigate the spin transport properties of molecules belonging to the acenes series by using density functional theory combined with the non-equilibrium Green's function approach to electronic transport. While short acenes are found to be non-magnetic, molecules comprising more than nine acene rings have a spin-polarized ground state. In their gas phase these have a singlet total spin configuration, since the two unpaired electrons occupying the doubly degenerate highest molecular orbital are antiferromagnetically coupled to each other. Such an orbital degeneracy is however lifted once the molecule is attached asymmetrically to Au electrodes via thiol linkers, leading to a fractional magnetic moment. In this situation the system Au/n-acene/Au can act as an efficient spin-filter with interesting applications in the emerging field of organic spintronics., Comment: 13 pages,9 figures

The Bethe-Ansatz local density approximation (LDA) to lattice density functional theory (LDFT) for the one-dimensional repulsive Hubbard model is extended to current-LDFT (CLDFT). The transport properties of mesoscopic Hubbard rings threaded by a magnetic flux are then systematically investigated by this scheme. In particular we present calculations of ground state energies, persistent currents and Drude weights for both a repulsive homogeneous and a single impurity Hubbard model. Our results for the ground state energies in the metallic phase compares favorably well with those obtained with numerically accurate many-body techniques. Also the dependence of the persistent currents on the Coulomb and the impurity interaction strength, and on the ring size are all well captured by LDA-CLDFT. Our study demonstrates that CLDFT is a powerful tool for studying one-dimensional correlated electron systems with high accuracy and low computational costs.

We present a Monte Carlo study of the finite temperature properties of an extended Hubbard-Peierls model describing one dimensional $\pi$-conjugated polymers. The model incorporates electron-phonon and hyperfine interaction and it is solved at the mean field level for half filling. In particular we explore the model as a function of the strength of electron-electron and electron-phonon interactions. At low temperature the system presents a diamagnetic to antiferromagnetic transition as the electron-electron interaction strength increases. At the same time by increasing the electron-phonon coupling there is a transition from a homogeneous to a Peierls dimerized geometry. As expected such a Peierls dimerized phase disappears at finite temperature as a result of thermal vibrations. More intriguing is the interplay between the electron-phonon and the electron-electron interactions at finite temperature. In particular we demonstrate that for a certain region of the parameter space there is a spin-crossover, where the system transits from a low-spin to a high-spin state as the temperature increases. In close analogy to standard spin-crossover in divalent magnetic molecules such a transition is entropy driven. Finally we discuss the role played by the hyperfine interaction over the phase diagram.

We present a theoretical study on the response properties to an external electric field of strongly correlated one-dimensional metals. Our investigation is based on the recently developed Bethe-Ansatz local density approximation (BALDA) to the density functional theory formulation of the Hubbard model. This is capable of describing both Luttinger liquid and Mott-insulator correlations. The BALDA calculated values for the static linear polarizability are compared with those obtained by numerically accurate methods, such as exact (Lanczos) diagonalization and the density matrix renormalization group, over a broad range of parameters. In general BALDA linear polarizabilities are in good agreement with the exact results. The response of the exact exchange and correlation potential is found to point in the same direction of the perturbing potential. This is well reproduced by the BALDA approach, although the fine details depend on the specific parameterization for the local approximation. Finally we provide a numerical proof for the non-locality of the exact exchange and correlation functional., Comment: 8 pages and 8 figures

The bias-dependent transport properties of short poly(G)-poly(C) A-DNA strands attached to Au electrodes are investigated with first principles electronic transport methods. By using the non- equilibrium Green's function approach combined with self-interaction corrected density functional theory, we calculate the fully self-consistent coherent I-V curve of various double-strand polymeric DNA fragments. We show that electronic wave-function localization, induced either by the native electrical dipole and/or by the electrostatic disorder originating from the first few water solvation layers, drastically suppresses the magnitude of the elastic conductance of A-DNA oligonucleotides. We then argue that electron transport through DNA is the result of sequence-specific short-range tunneling across a few bases combined with general diffusive/inelastic processes., Comment: 15 pages, 13 figures, 1 table

The effects of water wetting conditions on the transport properties of molecular nano-junctions are investigated theoretically by using a combination of classical molecular dynamics and first principles electronic transport calculations. These are at the level of the non-equilibrium Green's function method implemented for self-interaction corrected density functional theory. We find that water effectively produces electrostatic gating to the molecular junction, with a gating potential determined by the time-averaged water dipole field. Such a field is large for the polar benzene-dithiol molecule, resulting in a transmission spectrum shifted by about 0.6 eV with respect to that of the dry junction. The situation is drastically different for carbon nanotubes (CNTs). In fact, because of their hydro-phobic nature the gating is almost negligible, so that the average transmission spectrum of wet Au/CNT/Au junctions is essentially the same as that in dry conditions. This suggests that CNTs can be used as molecular interconnects also in water-wet situations, for instance as tips for scanning tunnel microscopy in solution or in biological sensors.

The electronic structure of the neutral and singly charged Mg vacancy in MgO is investigated using density functional theory. For both defects, semilocal exchange correlation functionals such as the local spin density approximation incorrectly predict a delocalized degenerate ground state. In contrast functionals that take strong correlation effects into account predict a localized solution, in agreement with spin resonance experiments. Our results, obtained with the HSE hybrid, atomic self-interaction corrected and LDA+U functionals, provide a number of constraints to the possibility of ferromagnetism in hole doped MgO.

Density functional theory calculations demonstrate that rocksalt MgN is a magnetic material at the verge of half-metallicity, with an electronic structure robust against strong correlations and spin-orbit interaction. Furthermore the calculated heat of formation describes the compound as metastable and suggests that it can be fabricated by tuning the relative Mg and N abundance during growth. Intriguingly the equilibrium lattice constant is close to that of MgO, so that MgN is likely to form as an inclusion during the fabrication of N-doped MgO. We then speculate that the MgO/MgN system may represent a unique materials platform for magnetic tunnel junctions not incorporating any transition metals.

We present results of electronic structure and transport calculations for metallic interfaces, based on density functional theory and the non-equilibrium Green's functions method. Starting from the electronic structure of smooth Al, Cu, Ag, and Au interfaces, we study the effects of different kinds of interface roughness on the transmission coefficient and the I-V characteristic. In particular, we compare prototypical interface distortions, including vacancies and metallic impurities., Comment: 15 pages, 10 figures

We demonstrate that the different magnetic states of a Mn12 molecule can be distinguished in a two-probe transport experiment from a complete knowledge of the current-voltage curve. Our results, obtained with state-of-the-art non-equilibrium transport methods combined with density functional theory, demonstrate that spin configuration-specific negative differential resistances (NDRs) appear in the I-V curves. These originate from the interplay between electron localization and the re-hybridization of the molecular levels in an external electric field and allow the detection of the molecule's spin-state., Comment: 8 pages, 8 figures, 2 page supplementary material

Simulations of scanning tunneling microscopy measurements for molecules on surfaces are traditionally based on a perturbative approach, most typically employing the Tersoff-Hamann method. This assumes that the STM tip is far from the sample so that the two do not interact with each other. However, when the tip gets close to the molecule to perform measurements, the electrostatic interplay between the tip and substrate may generate non-trivial potential distribution, charge transfer and forces, all of which may alter the electronic and physical structure of the molecule. These effects are investigated with the ab initio quantum transport code SMEAGOL, combining non-equilibrium Green's functions formalism with density functional theory. In particular, we investigate alkanethiol molecules terminated with either CH3 or CF3 end-groups on gold surfaces, for which recent experimental data are available. We discuss the effects connected to the interaction between the STM tip and the molecule, as well as the asymmetric charge transfer between the molecule and the electrodes., Comment: 10 pages, 18 figures

The interface structure of Fe/MgO(100) magnetic tunnel junctions predicted by density functional theory (DFT) depends significantly on the choice of exchange and correlation functional. Bader analysis reveals that structures obtained by relaxing the cell with the local spin-density approximation (LSDA) display a different charge transfer than those relaxed with the generalized gradient approximation (GGA). As a consequence, the electronic transport is found to be extremely sensitive to the interface structure. In particular, the conductance for the LSDA-relaxed geometry is about one order of magnitude smaller than that of the GGA-relaxed one. The high sensitivity of the electronic current to the details of the interface might explain the discrepancy between the experimental and calculated values of magnetoresistance., Comment: Submitted to PRL, 5 figures

In this work, we calculate with ab initio methods the current-voltage characteristics for ideal single- and double-barrier Fe/MgO (001) magnetic tunnel junctions. The current is calculated in the phase-coherent limit by using the recently developed SMEAGOL code, combining the nonequilibrium Green function formalism with density-functional theory. In general we find that double-barrier junctions display a larger magnetoresistance, which decays with bias at a slower pace than their single-barrier counterparts. This is explained in terms of enhanced spin filtering from the middle Fe layer sandwiched in between the two MgO barriers. In addition, for double-barrier tunnel junctions, we find a well defined peak in the magnetoresistance at a voltage of V=0.1 V. This is the signature of resonant tunneling across a majority quantum well state. Our findings are discussed in relation to recent experiments.

Predicting magnetism originating from 2$p$ orbitals is a delicate problem, which depends on the subtle interplay between covalency and Hund's coupling. Calculations based on density functional theory and the local spin density approximation fail in two remarkably different ways. On the one hand the excessive delocalization of spin-polarized holes leads to half-metallic ground states and the expectation of room temperature ferromagnetism. On the other hand, in some cases a magnetic ground state may not be predicted at all. We demonstrate that a simple self-interaction correction scheme modifies both these situations via an enhanced localization of the holes responsible for the magnetism and possibly Jahn-Teller distortion. In both cases the ground state becomes insulating and the magnetic coupling between the impurities weak., Comment: 4 pages, 4 figures

Magnetic 3d-ions doped into wide-gap oxides show signatures of room temperature ferromagnetism, although their concentration is two orders of magnitude smaller than that in conventional magnets. The prototype of these exceptional materials is Co-doped ZnO, for which an explanation of the room temperature ferromagnetism is still elusive. Here we demonstrate that magnetism originates from Co2+ oxygen-vacancy pairs with a partially filled level close to the ZnO conduction band minimum. The magnetic interaction between these pairs is sufficiently long-ranged to cause percolation at moderate concentrations. However, magnetically correlated clusters large enough to show hysteresis at room temperature already form below the percolation threshold and explain the current experimental findings. Our work demonstrates that the magnetism in ZnO:Co is entirely governed by intrinsic defects and a phase diagram is presented. This suggests a recipe for tailoring the magnetic properties of spintronics materials by controlling their intrinsic defects.

In transport calculations for molecular junctions based on density functional theory the choice of exchange and correlation functional may dramatically affect the results. In particular local and semi-local functionals tend to over-delocalize the molecular levels thus artificially increasing their broadening. In addition the same molecular levels are usually misplaced with respect to the Fermi level of the electrodes. These shortfalls are reminiscent of the inability of local functionals to describe Mott-Hubbard insulators, but they can be corrected with a simple and computationally undemanding self-interaction correction scheme. We apply such a scheme, as implemented in our transport code Smeagol, to a variety of phenyl-based molecular junctions attached to gold electrodes. In general the corrections reduce the current, since the resonant Kohn-Sham states of the molecule are shifted away from the contact Fermi level. In contrast, when the junction is already described as insulating by local exchange and correlation potentials, the corrections are minimal and the I-V is only weakly modified., Comment: 12 pages, 21 figures

Standard density functional approximations greatly over-estimate the static polarizability of longchain polymers, but Hartree-Fock or exact exchange calculations do not. Simple self-interaction corrected (SIC) approximations can be even better than exact exchange, while their computational cost can scale only linearly with the number of occupied orbitals., Comment: 4 pages, 3 figures

The approximate atomic self-interaction corrections (ASIC) method to density functional theory is put to the test by calculating the exchange interaction for a number of prototypical materials, critical to local exchange and correlation functionals. ASIC total energy calculations are mapped onto an Heisenberg pair-wise interaction and the exchange constants J are compared to those obtained with other methods. In general the ASIC scheme drastically improves the bandstructure, which for almost all the cases investigated resemble closely available photo-emission data. In contrast the results for the exchange parameters are less satisfactory. Although ASIC performs reasonably well for systems where the magnetism originates from half-filled bands, it suffers from similar problems than those of LDA for other situations. In particular the exchange constants are still overestimated. This reflects a subtle interplay between exchange and correlation energy, not captured by the ASIC., Comment: 10 pages

Electronic transport in nickel magnetic point contacts is investigated with a combination of density functional theory and the non-equilibrium Green functions method. In particular we address the possibility of huge ballistic magnetoresistance in impurity-free point contacts and the effects of oxygen impurities. On-site corrections over the local spin density approximation (LSDA) for the exchange and correlation potential, namely the LDA+U method, are applied in order to account for low-coordination and strong correlations. We show that impurity-free point contacts present magnetoresistance never in excess of 50%. This value can raise up to about 450 % in the case of oxygen contamination. These results suggest that magnetoresistance in excess of 1,000 % can not have solely electronic origin., Comment: 9 pages, 12 figures

Density functional theory calculations of electronic transport based on local exchange and correlation functionals contain self-interaction errors. These originate from the interaction of an electron with the potential generated by itself and may be significant in metal-molecule-metal junctions due to the localized nature of the molecular orbitals. As a consequence, insulating molecules in weak contact with metallic electrodes erroneously form highly conducting junctions, a failure similar to the inability of local functionals of describing Mott-Hubbard insulators. Here we present a fully self-consistent and still computationally undemanding self-interaction correction scheme that overcomes these limitations. The method is implemented in the Green's function non-equilibrium transport code Smeagol and applied to the prototypical cases of benzene molecules sandwiched between gold electrodes. The self-interaction corrected Kohn-Sham highest occupied molecular orbital now reproduces closely the negative of the molecular ionization potential and is moved away from the gold Fermi energy. This leads to a drastic reduction of the low bias current in much better agreement with experiments., Comment: 4 pages, 5 figures

We present an atomic orbital based approximate scheme for self-interaction correction (SIC) to the local density approximation of density functional theory. The method, based on the idea of Filippetti and Spaldin [Phys. Rev. B 67, 125109 (2003)], is implemented in a code using localized numerical atomic orbital basis sets and is now suitable for both molecules and extended solids. After deriving the fundamental equations as a non-variational approximation of the self-consistent SIC theory, we present results for a wide range of molecules and insulators. In particular, we investigate the effect of re-scaling the self-interaction correction and we establish a link with the existing atomic-like corrective scheme LDA+U. We find that when no re-scaling is applied, i.e. when we consider the entire atomic correction, the Kohn-Sham HOMO eigenvalue is a rather good approximation to the experimental ionization potential for molecules. Similarly the HOMO eigenvalues of negatively charged molecules reproduce closely the molecular affinities. In contrast a re-scaling of about 50% is necessary to reproduce insulator bandgaps in solids, which otherwise are largely overestimated. The method therefore represents a Kohn-Sham based single-particle theory and offers good prospects for applications where the actual position of the Kohn-Sham eigenvalues is important, such as quantum transport., Comment: 16 pages, 7 figures

We present a study of the effects of inelastic scattering on the transport properties of various nanoscale devices, namely H$_2$ molecules sandwiched between Pt contacts, and a spin-valve made by an organic molecule attached to model half-metal ferromagnetic current/voltage probes. In both cases we use a tight-binding Su-Schrieffer-Heeger Hamiltonian and the inelastic effects are treated with a multi-channel method, including Pauli exclusion principle. In the case of the H$_2$ molecule, we find that inelastic backscattering is responsible for the drop of the differential conductance at biases larger than the excitation energy of the lower of the molecular phonon modes. In the case of the spin-valve, we investigate the different spin-currents and the magnetoresistance as a function of the position of the Fermi level with respect to the spin-polarized band edges. In general inelastic scattering reduces the spin-polarization of the current and consequently the magnetoresistance.

{\it Ab initio} computational methods for electronic transport in nanoscaled systems are an invaluable tool for the design of quantum devices. We have developed a flexible and efficient algorithm for evaluating $I$-$V$ characteristics of atomic junctions, which integrates the non-equilibrium Green's function method with density functional theory. This is currently implemented in the package {\it Smeagol}. The heart of {\it Smeagol} is our novel scheme for constructing the surface Green's functions describing the current/voltage probes. It consists of a direct summation of both open and closed scattering channels together with a regularization procedure of the Hamiltonian, and provides great improvements over standard recursive methods. In particular it allows us to tackle material systems with complicated electronic structures, such as magnetic transition metals. Here we present a detailed description of {\it Smeagol} together with an extensive range of applications relevant for the two burgeoning fields of spin and molecular-electronics., Comment: For additional materials on this topic please follow http://www.smeagol.tcd.ie

All density functional calculations of single-molecule transport to date have used continuous exchange-correlation approximations. The lack of derivative discontinuity in such calculations leads to the erroneous prediction of metallic transport for insulating molecules. A simple and computationally undemanding atomic self-interaction correction greatly improves the agreement with experiment for the prototype Au/dithiolated-benzene/Au junction., Comment: 4 pages. Also available at http://www.smeagol.tcd.ie

Using first principles simulations we perform a detailed study of the structural, electronic and transport properties of monoatomic platinum chains, sandwiched between platinum electrodes. First we demonstrate that the most stable atomic configuration corresponds to a zigzag arrangement that gradually straightens as the chains are stretched. Secondly, we find that the conductance at equilibrium atomic spacing does not oscillate with the number of atoms $n$ in the chain, but instead decreases almost monotonically with $n$. In contrast, the conductances of chains of fixed $n$ oscillate as the end atoms are pulled apart, due to the gradual closing and opening of conductance channels as the chain straightens., Comment: 4 pages. 5 Figures. Version with a new reference

Recent advances in atomic and nano-scale growth and characterization techniques have led to the production of modern magnetic materials and magneto-devices which reveal a range of new fascinating phenomena. The modeling of these is a tough theoretical challenging since one has to describe accurately both electronic structure of the constituent materials, and their transport properties. In this paper I review recent advances in modeling spin-transport at the atomic scale using first-principles methods, focusing both on the methodological aspects and on the applications. The review, which is designed as tutorial for students at postgraduate level, is structured in six main sections: 1) Introduction, 2) General concepts in spin-transport, 3) Transport Theory: Linear Response, 4) Transport Theory: Non-equilibrium Transport, 5) Results, 6) Conclusion. In addition an overview of the computational codes available to date is also included., Comment: To be published in "Handbook of Computational Nanotechnology" (American Scientific Publishers, California 2004). Full text, also available at http://www.smeagol.tcd.ie 117 Pages, 42 figures, 270 references

In view of the recent experimental reports of unexpected ferromagnetism in HfO$_{2}$ thin films, we carried out first principles investigations looking for magnetic order possibly brought about by the presence of small concentrations of intrinsic point defects. {\it Ab initio} electronic structure calculations using density functional theory (DFT) show that isolated cation vacancy sites in HfO$_{2}$ lead to the formation of high spin defect states. Furthermore these appear to be ferromagnetically coupled with a rather short range magnetic interaction, resulting in a ferromagnetic ground state for the whole system. More interestingly, the occurrence of these high spin states and ferromagnetism is in the low symmetry monoclinic phase of HfO$_{2}$. This is radically different from other systems previously known to exhibit point defect ferromagnetism, warranting a closer look at the phenomenon., Comment: 4 pages, 6 figures

We report a detailed theoretical study of the bonding and conduction properties of an hydrogen molecule joining either platinum or palladium electrodes. We show that an atomic arrangement where the molecule is placed perpendicular to the electrodes is unstable for all distances between electrodes. In contrast, the configuration where the molecule bridges the electrodes is stable in a wide range of distances. In this last case the bonding state of the molecule does not hybridize with the leads and remains localized within the junction. As a result, this state does not transmit charge so that electronic transport is carried only through the anti-bonding state. This fact leads to conductances of 1 $G_0$ at most, where $G_0=2e^2/h$. We indeed find that G is equal to 0.9 and 0.6 $G_0$ for Pt and Pd contacts respectively., Comment: 6 pages, 6 figures. Submitted to Phys. Rev. B

We present a theoretical study of the transport properties of magnetic point contacts under bias. Our calculations are based on the Keldish's non-equilibrium Green's function formalism combined with a self-consistent empirical tight-binding Hamiltonian, which describes both strong ferromagnetism and charging effects. We demonstrate that large magnetoresistance solely due to electronic effects can be found when a sharp domain wall forms inside a magnetic atomic-scale point contact. Moreover we show that the symmetry of the $I$-$V$ characteristic depends on the position of the domain wall in the constriction. In particular diode-like curves can arise when the domain wall is placed off-center within the point contact, although the whole structure does not present any structural asymmetry., Comment: 7 figures, submitted to PRB

We present an extensive density functional theory study of the electronic, magnetic and transport properties of GaAs and AlAs digital ferromagnetic heterostructures. These can be obtained by $\delta$-doping with Mn the GaAs layers of a GaAs/AlAs superlattice. Our analysis spans a range of Mn concentrations and considers the presence of compensating defects such as As antisites. In the defect-free case all the heterostructures studied present an half-metallic electronic structure. In contrast when As antisites are present the half-metallic state is destroyed and the heterostructures behave as dirty planar metals. In this case they show a large $p$-type metallic conductance in the Mn plane mainly due to majority spin electrons, and an $n$-type hopping-like conductance in the GaAs planes mainly due to minority spin electrons. This suggests that if the As antisites can be kept far from the Mn planes, spatial separation of the different spin currents can be achieved. Finally we show that in the case of AlAs/(Ga,Mn)As digital ferromagnetic heterostructures the AlAs/GaAs valence band offset produces an additional confining potential for the holes responsible for the ferromagnetism. Therefore the ferromagnetic coupling between the Mn ions becomes larger and more robust to the presence of As antisites., Comment: 14 pages, 14 figures

We investigate the effect of quantum confinement on the ferromagnetism of diluted magnetic semiconductor Ga$_{1-x}$Mn$_x$As using a combination of tight-binding and density functional methods. We observe strong majority-spin Mn $d$-As $p$ hybridization, as well as half metallic behavior, down to sizes as small as 20 \AA in diameter. Below this critical size, the doped holes are self-trapped by the Mn-sites, signalling both valence and electronic transitions. Our results imply that magnetically doped III-V nanoparticles will provide a medium for manipulating the electronic structure of dilute magnetic semiconductors while conserving the ferromagnetic properties and even enhancing it in certain size regime., Comment: 4 pages, 3 figures

In this paper we review the latest achievements of density functional theory in understanding the physics of diluted magnetic semiconductors. We focus on transition metal doped III-V semiconductors, which show spontaneous ferromagnetic order at relatively high temperature and good structural compatibility with existing III-V devices. We show that density functional theory is a very powerful tool for i) studying the effects of local doping defects and disorder on the magnetic properties of these materials, ii) predicting properties of new materials and iii) providing parameters, often not accessible from experiments, for use in model Hamiltonian calculations. Such studies are facilitated by recent advances in numerical implementations of density functional theory, which make the study of systems with a very large number of atoms possible., Comment: Invited article for Journal of Superconductivity and Novel Magnetic Materials

We study spin-dependent electronic transport across ferromagnet/superconductor ballistic junctions modeled using tight-binding Hamiltonians with s, p and d orbitals and material-specific parameters. We show that by accurately modeling the band structure of the bulk materials, one can reproduce the measured differential conductance of Cu/Pb nanocontacts. In contrast the differential conductance of Co/Pb contacts can only be reproduced if an enhanced magnetic moment is present at the interface., Comment: 4 pages, 4 figures

We study the Giant Magnetoresistance (GMR) ratio in magnetic multiayers with a single superconducting contact in the presence of spin-mixing processes. It has been recently shown that the GMR ratio of magnetic multilayers is strongly suppressed by the presence of a superconducting contact when spin-flipping is not allowed. In this Letter we demonstrate that the GMR ratio can be dramatically enhanced by spin-orbit interaction and/or non-collinear magnetic moments. The system is described using a tight-binding model with either s-p-d or s-d atomic orbitals per site., Comment: 5 pages, 5 figures, submitted to Phys. Rev. Lett

Using a scattering technique based on a parametrized linear combination of atomic orbitals Hamiltonian, we calculate the ballistic quantum conductance of multi-wall carbon nanotubes. We find that inter-wall interactions not only block some of the quantum conductance channels, but also redistribute the current non-uniformly over the individual tubes across the structure. Our results provide a natural explanation for the unexpected integer and non-integer conductance values reported for multi-wall nanotubes in Ref.[S. Frank, P. Poncharal, Z.L. Wang, and W. A. de Heer, Science 280, 1744 (1998)], Comment: 16 pages, 5 figures, submitted to PRL