Orbital-polarization terms: From a phenomenological to a first-principles description of orbital magnetism in density-functional theory

Phenomenological orbital-polarizition (OP) terms have been repeatedly introduced in the single-particle equations of spin-density-functional theory, in order to improve the description of orbital magnetic moments in systems containing transition metal ions. Here we show that these ad hoc corrections can be interpreted as approximations to the exchange-correlation vector potential A(xc) of current-density functional theory (CDFT). This connection provides additional information on both approaches: phenomenological OP terms are connected to first-principles theory, leading to a rationale for their empirical success and a reassessment of their limitations and the approximations made in their derivation. Conversely, the connection of OP terms with CDFT leads to a set of simple approximations to the CDFT potential A(xc), with a number of desirable features that are absent from electron-gas-based functionals. (C) 2008 Wiley Periodicals, Inc.

Orbital Magnetization of Quantum Spin Hall Insulator Nanoparticles

Both spin and orbital degrees of freedom contribute to the magnetic moment of isolated atoms. However, when inserted in crystals, atomic orbital moments are quenched because of the lack of rotational symmetry that protects them when isolated. Thus, the dominant contribution to the magnetization of magnetic materials comes from electronic spin. Here we show that nanoislands of quantum spin Hall insulators can host robust orbital edge magnetism whenever their highest occupied Kramers doublet is singly occupied, upgrading the spin edge current into a charge current. The resulting orbital magnetization scales linearly with size, outweighing the spin contribution for islands of a few nm in size. This linear scaling is specific of the Dirac edge states and very different from Schrodinger electrons in quantum rings. By modeling Bi(111) flakes, whose edge states have been recently observed, we show that orbital magnetization is robust with respect to disorder, thermal agitation, shape of the island, and crystallographic direction of the edges, reflecting its topological protection.

Giant diamagnetism of gold nanorods

The presence of giant diamagnetism in Au nanorods, NRs, is shown to be a possible consequence of field induced currents in the surface electrons. The distance, Delta , between quantum surface energy levels has been calculated as a function of the NRs radius. Note that those electrons occupying states for which Delta > k(B)T are steadily orbiting with constant orbital moment. The diamagnetic response induced when a field is turned on remains constant during the time the field is acting. As the NRs radius increases, Delta decreases and accordingly the electron fraction available to generate constant currents decreases, consequently the surface diamagnetic susceptibility decreases towards its bulk value. The surface electronic motion induced by the axial applied field on electrons confined into a cylindrical surface accounts with extremely good quantitative agreement for the giant diamagnetism recently measured and reported.

EXACT AND APPROXIMATE RELATIONS FOR THE SPIN-DEPENDENCE OF THE EXCHANGE ENERGY IN HIGH MAGNETIC FIELDS

The exchange energy of an arbitrary collinear-spin many-body system in an external magnetic field is a functional of the spin-resolved charge and current densities, E(x)[n(up arrow), n(down arrow), j(up arrow), j(down arrow)]. Within the framework of density-functional theory (DFT), we show that the dependence of this functional on the four densities can be fully reconstructed from either of two extreme limits: a fully polarized system or a completely unpolarized system. Reconstruction from the limit of an unpolarized system yields a generalization of the Oliver-Perdew spin scaling relations from spin-DFT to current-DFT. Reconstruction from the limit of a fully polarized system is used to derive the high-field form of the local-spin-density approximation to current-DFT and to magnetic-field DFT.

KO2: Realization of Orbital Ordering in a p-Orbital System

KO2 is a molecular solid consisting of oxygen dimers. K present in the lattice donates an electron which goes on to occupy the O p levels.As the basic electronic structure is similar to that of an oxygen molecule, except for broadening due to solid state effects, KO2 represents the realization of the doping of oxygen molecules arranged in a lattice. These considerations alone result in magnetism with high ordering temperatures as our calculations reveal. However, we find that the high temperature structure is unstable to an orbital ordering (OO) transition. The microscopic considerations driving the OO transition, however, are electrostatic interactions instead of the often encountered superexchange driven ordering within the Kugel-Khomskii model often used to describe the OO. This OO transition is also found to preclude any possibility of high magnetic ordering temperatures, which otherwise seemed possible.

Origins of magnetism in transition metal doped Cul

Cupric iodide is a p-type semiconductor and has a large band gap. Doping of Mn, Co, and Ni are found to make gamma-CuI ferromagnetic ground state, while Cr-doped and Fe-doped CuI systems are stabilized in antiferromagnetic configurations. The origins of the magnetic ordering are demonstrated successfully by the phenomenological band coupling model based on d-d level repulsions between the dopant ions. Furthermore, using a molecular-orbital bonding model, the electronic structures of the doped CuI are well understood. According to Heisenberg model, high-T-C may be expected for CuI:Mn and CuI:Ni if there are no native defects or other impurities.

Competitions of magnetism and superconductivity in FeAs-based materials

Using the numerical unrestricted Hartree-Fock approach, we study the ground state of a two-orbital model describing newly discovered FeAs-based superconductors. We observe the competition of a (0, π) mode spin-density wave and the superconductivity as the doping concentration changes. There might be a small region in the electron-doping side where the magnetism and superconductivity coexist. The superconducting pairing is found to be spin singlet,orbital even, and coexisting sxy + dx~2-y~2 wave (even parity).

Magnetism of Co overlayers and nanostructures on W(001): A first-principles study

The magnetic properties of Co nanostructures and a Co monolayer on W(001) have been studied in the framework of density functional theory. Different geometries such as planar and three-dimensional clusters have been considered, with cluster sizes varying between 2 and 13 atoms. The calculations were performed using the real-space linear muffin-tin orbital method (RS-LMTO-ASA). With respect to the stability of the magnetic state, we predict an antiferromagnetic (AFM) structure for the ground state of the planar Co clusters and a ferromagnetic (FM) state for the three-dimensional clusters. For the three-dimensional clusters, one of the AFM arrangements leads to frustration due to the competing FM and AFM exchange interactions between different atoms in the cluster, and gives rise to a non-collinear state with energy close to that of the FM ground state. The relative role of the Co-Co and Co-W exchange interactions is also investigated. (C) 2007 Elsevier B.V. All rights reserved.

Magnetism of C adatoms on BN nanostructures : implications for functional nanodevices

Spin-polarized density functional calculations reveal that magnetism can be induced by carbon adatoms on boron nitride nanotubes (BNNTs) and BN hexagonal sheets. As a result of the localization of impurity states, these hybrid sp-electron systems are spin-polarized, with a local magnetic moment of 2.0 μB per C adatom regardless of the tube diameter and the bonding between the C atom and the BNNTs/BN sheets. An analysis of orbital hybridization indicates that two valence electrons participate in the bonding and the remaining two electrons of the C adatom are confined at the adsorption site and contribute to the magnetism accordingly. The effective interaction distance between the C-induced magnetic moments is evaluated. In terms of the diffusion barrier and the adsorption energy of C adatoms on the BN nanotubes/ sheets, a fabrication method for BN-C-based functional nanodevices is proposed, and a series of virtual building blocks for functional devices are illustrated.

First-principles studies of complex magnetism in Mn nanostructures on the Fe(001) surface

The magnetic properties of Mn nanostructures on the Fe(001) surface have been studied using the noncollinear first-principles real space-linear muffin-tin orbital-atomic sphere approximation method within density-functional theory. We have considered a variety of nanostructures such as adsorbed wires, pyramids, and flat and intermixed clusters of sizes varying from two to nine atoms. Our calculations of interatomic exchange interactions reveal the long-range nature of exchange interactions between Mn-Mn and Mn-Fe atoms. We have found that the strong dependence of these interactions on the local environment, the magnetic frustration, and the effect of spin-orbit coupling lead to the possibility of realizing complex noncollinear magnetic structures such as helical spin spiral and half-skyrmion.

Materials properties from electron density distributions: Molecular magnetism and linear optical properties

The general goal of this thesis is correlating observable properties of organic and metal-organic materials with their ground-state electron density distribution. In a long-term view, we expect to develop empirical or semi-empirical approaches to predict materials properties from the electron density of their building blocks, thus allowing to rationally engineering molecular materials from their constituent subunits, such as their functional groups. In particular, we have focused on linear optical properties of naturally occurring amino acids and their organic and metal-organic derivatives, and on magnetic properties of metal-organic frameworks. For analysing the optical properties and the magnetic behaviour of the molecular or sub-molecular building blocks in materials, we mostly used the more traditional QTAIM partitioning scheme of the molecular or crystalline electron densities, however, we have also investigated a new approach, namely, X-ray Constrained Extremely Localized Molecular Orbitals (XC-ELMO), that can be used in future to extracted the electron densities of crystal subunits. With the purpose of rationally engineering linear optical materials, we have calculated atomic and functional group polarizabilities of amino acid molecules, their hydrogen-bonded aggregates and their metal-organic frameworks. This has enabled the identification of the most efficient functional groups, able to build-up larger electric susceptibilities in crystals, as well as the quantification of the role played by intermolecular interactions and coordinative bonds on modifying the polarizability of the isolated building blocks. Furthermore, we analysed the dependence of the polarizabilities on the one-electron basis set and the many-electron Hamiltonian. This is useful for selecting the most efficient level of theory to estimate susceptibilities of molecular-based materials. With the purpose of rationally design molecular magnetic materials, we have investigated the electron density distributions and the magnetism of two copper(II) pyrazine nitrate metal-organic polymers. High-resolution X-ray diffraction and DFT calculations were used to characterize the magnetic exchange pathways and to establish relationships between the electron densities and the exchange-coupling constants. Moreover, molecular orbital and spin-density analyses were employed to understand the role of different magnetic exchange mechanisms in determining the bulk magnetic behaviour of these materials. As anticipated, we have finally investigated a modified version of the X-ray constrained wavefunction technique, XC-ELMOs, that is not only a useful tool for determination and analysis of experimental electron densities, but also enables one to derive transferable molecular orbitals strictly localized on atoms, bonds or functional groups. In future, we expect to use XC-ELMOs to predict materials properties of large systems, currently challenging to calculate from first-principles, such as macromolecules or polymers. Here, we point out advantages, needs and pitfalls of the technique. This work fulfils, at least partially, the prerequisites to understand materials properties of organic and metal-organic materials from the perspective of the electron density distribution of their building blocks. Empirical or semi-empirical evaluation of optical or magnetic properties from a preconceived assembling of building blocks could be extremely important for rationally design new materials, a field where accurate but expensive first-principles calculations are generally not used. This research could impact the community in the fields of crystal engineering, supramolecular chemistry and, of course, electron density analysis.

Sedimentology, mineral magnetism, andoxygen isotope stratigraphy of ODP Hole 133-819A

To investigate late Quaternary paleoclimatic and paleoceanographic change in the sedimentary record, preserved on the Australian Continental Margin during the late Quaternary, core material was collected from Ocean Drilling Program, Leg 133, Site 819. An expanded sequence of late Quaternary, rhythmically bedded, predominantly hemipelagic sediments were recovered from Hole 819A. The foraminiferal d18O record preserved at Hole 819A suggests that the late Quaternary section is incomplete. Both benthic and planktonic d18O stratigraphies can be traced tentatively downcore to stage 6 at about 32.5 mbsf, where a major hiatus occurs. At this level, a slump detachment surface has been identified (Shipboard Scientific Party, 1991). This slump has removed marine oxygen isotope stages 7 to 13. Below 32.5 mbsf, continuous correlation can be achieved in the planktonic d18O curve, with existing deep-sea foraminiferal oxygen isotope stratigraphies from stage 14 through stage 28. The major hiatus at 32.5 mbsf marks the position of a significant change in the character of the sedimentation at Site 819. Sediments below 32.5 mbsf, relative to those above 32.5 mbsf, are characterized by less variation in mean particle size; lower percentages of carbonate content in the coarse fraction (>63 µm); a stronger relationship between the percentage of fine fraction and magnetic mineral concentration, and lower foraminiferal abundances. Above the hiatus, large fluctuations in mean particle size occurred, which have been interpreted to be the result of high foraminiferal abundances. Early highstands show high terrigenous influx in the fine fraction above the hiatus. This is the opposite of the general idea of high terrigenous influx during lowstands of sea level on siliciclastic dominated continental margins. We are far from understanding the origin of this material and further investigation will be required (see also Glenn et al., this volume). All our records, except the planktonic foraminiferal oxygen isotope record, indicate that the major hiatus marks the position of a significant change in the environment at Site 819. The planktonic foraminiferal d18O record suggests that environmental change occurred prior to the formation of the hiatus (i.e., near the Brunhes/Matuyama [B/M] boundary). The interval between the B/M boundary and the hiatus represents a transitional period between two different patterns of ocean circulation. Throughout most of the lower part of the sequence, Site 819 was at a shallow-water depth and local oceanographic conditions were dominated by sluggish Subtropical Central Water (SCW) flow. However, near the B/M boundary, ocean circulation patterns intensified, reflecting a worldwide change in paleoenvironment. Enhanced ocean circulation patterns were possibly aided by tectonic subsidence. During this period Site 819 became progressively more under the influence of Antarctic Intermediate Water (AAIW), than SCW. In the upper part of the sequence at Hole 819 A, we see a continuation of the pattern of oceanographic reorganization suggested during stages 21 through 14. Intensification of the subsurface oceanographic circulation was also accompanied by the progressive wedging southward of surface waters associated with the East Australian Current (EAC). The change in the nature of the records in the lower and upper parts of the sequence at Site 819 are thought to reflect perturbations by the orbital eccentricity cycle.

Magnetism in graphene nanoislands

We study the magnetic properties of nanometer-sized graphene structures with triangular and hexagonal shapes terminated by zigzag edges. We discuss how the shape of the island, the imbalance in the number of atoms belonging to the two graphene sublattices, the existence of zero-energy states, and the total and local magnetic moment are intimately related. We consider electronic interactions both in a mean-field approximation of the one-orbital Hubbard model and with density functional calculations. Both descriptions yield values for the ground state total spin S consistent with Lieb’s theorem for bipartite lattices. Triangles have a finite S for all sizes whereas hexagons have S=0 and develop local moments above a critical size of ≈1.5  nm.

Time-Dependent Density Functional Molecular Orbital and Excited State Calculations on Bis(porphyrinyl)butadiynes in the Monocationic, Neutral, Monoanionic, and Dianionic Oxidation States

We report a theoretical study of the multiple oxidation states (1+, 0, 1−, and 2−) of a meso,meso-linked diporphyrin, namely bis[10,15,20-triphenylporphyrinatozinc(II)-5-yl]butadiyne (4), using Time-Dependent Density Functional Theory (TDDFT). The origin of electronic transitions of singlet excited states is discussed in comparison to experimental spectra for the corresponding oxidation states of the close analogue bis{10,15,20-tris[3‘,5‘-di-tert-butylphenyl]porphyrinatozinc(II)-5-yl}butadiyne (3). The latter were measured in previous work under in situ spectroelectrochemical conditions. Excitation energies and orbital compositions of the excited states were obtained for these large delocalized aromatic radicals, which are unique examples of organic mixed-valence systems. The radical cations and anions of butadiyne-bridged diporphyrins such as 3 display characteristic electronic absorption bands in the near-IR region, which have been successfully predicted with use of these computational methods. The radicals are clearly of the “fully delocalized” or Class III type. The key spectral features of the neutral and dianionic states were also reproduced, although due to the large size of these molecules, quantitative agreement of energies with observations is not as good in the blue end of the visible region. The TDDFT calculations are largely in accord with a previous empirical model for the spectra, which was based simplistically on one-electron transitions among the eight key frontier orbitals of the C4 (1,4-butadiyne) linked diporphyrins.

Super-orbital re-entry in Australia : laboratory measurement, simulation and flight observation

There are large uncertainties in the aerothermodynamic modelling of super-orbital re-entry which impact the design of spacecraft thermal protection systems (TPS). Aspects of the thermal environment of super-orbital re-entry flows can be simulated in the laboratory using arc- and plasma jet facilities and these devices are regularly used for TPS certification work [5]. Another laboratory device which is capable of simulating certain critical features of both the aero and thermal environment of super-orbital re-entry is the expansion tube, and three such facilities have been operating at the University of Queensland in recent years[10]. Despite some success, wind tunnel tests do not achieve full simulation, however, a virtually complete physical simulation of particular re-entry conditions can be obtained from dedicated flight testing, and the Apollo era FIRE II flight experiment [2] is the premier example which still forms an important benchmark for modern simulations. Dedicated super-orbital flight testing is generally considered too expensive today, and there is a reluctance to incorporate substantial instrumentation for aerothermal diagnostics into existing missions since it may compromise primary mission objectives. An alternative approach to on-board flight measurements, with demonstrated success particularly in the ‘Stardust’ sample return mission, is remote observation of spectral emissions from the capsule and shock layer [8]. JAXA’s ‘Hayabusa’ sample return capsule provides a recent super-orbital reentry example through which we illustrate contributions in three areas: (1) physical simulation of super-orbital re-entry conditions in the laboratory; (2) computational simulation of such flows; and (3) remote acquisition of optical emissions from a super-orbital re entry event.