The elevated Curie temperature and half-metallicity in the ferromagnetic semiconductor La$_{x}$Eu$_{1-x}$O

Here we study the effect of La doping in EuO thin films using SQUID magnetometry, muon spin rotation ($mu$SR), polarized neutron reflectivity (PNR), and density functional theory (DFT). The $mu$SR data shows that the La$_{0.15}$Eu$_{0.85}$O is homogeneously magnetically ordered up to its elevated $T_{rm C}$. It is concluded that bound magnetic polaron behavior does not explain the increase in $T_{rm C}$ and an RKKY-like interaction is consistent with the $mu$SR data. The estimation of the magnetic moment by DFT simulations concurs with the results obtained by PNR, showing a reduction of the magnetic moment per La$_{x}$Eu$_{1-x}$O for increasing lanthanum doping. This reduction of the magnetic moment is explained by the reduction of the number of Eu-4$f$ electrons present in all the magnetic interactions in EuO films. Finally, we show that an upwards shift of the Fermi energy with La or Gd doping gives rise to half-metallicity for doping levels as high as 3.2 %.

Simulations of Nanocrystals Under Pressure: Combining Electronic Enthalpy and Linear-Scaling Density-Functional Theory

We present an implementation in a linear-scaling density-functional theory code of an electronic enthalpy method, which has been found to be natural and efficient for the ab initio calculation of finite systems under hydrostatic pressure. Based on a definition of the system volume as that enclosed within an electronic density isosurface [Phys. Rev. Lett., 94, 145501 (2005)], it supports both geometry optimizations and molecular dynamics simulations. We introduce an approach for calibrating the parameters defining the volume in the context of geometry optimizations and discuss their significance. Results in good agreement with simulations using explicit solvents are obtained, validating our approach. Size-dependent pressure-induced structural transformations and variations in the energy gap of hydrogenated silicon nanocrystals are investigated, including one comparable in size to recent experiments. A detailed analysis of the polyamorphic transformations reveals three types of amorphous structures and their persistence on depressurization is assessed.

Linear-scaling time-dependent density-functional theory in the linear response formalism

We present an implementation of time-dependent density-functional theory (TDDFT) in the linear response formalism enabling the calculation of low energy optical absorption spectra for large molecules and nanostructures. The method avoids any explicit reference to canonical representations of either occupied or virtual Kohn-Sham states and thus achieves linear-scaling computational effort with system size. In contrast to conventional localised orbital formulations, where a single set of localised functions is used to span the occupied and unoccupied state manifold, we make use of two sets of in situ optimised localised orbitals, one for the occupied and one for the unoccupied space. This double representation approach avoids known problems of spanning the space of unoccupied Kohn-Sham states with a minimal set of localised orbitals optimised for the occupied space, while the in situ optimisation procedure allows for efficient calculations with a minimal number of functions. The method is applied to a number of medium sized organic molecules and a good agreement with traditional TDDFT methods is observed. Furthermore, linear scaling of computational cost with system size is demonstrated on a system of carbon nanotubes.

Bromophenyl functionalization of carbon nanotubes : an ab initio study

We study the thermodynamics of bromophenyl functionalization of carbon nanotubes with respect to diameter and metallic/insulating character using density-functional theory (DFT). On one hand, we show that the activation energy for the grafting of a bromophenyl molecule onto a semiconducting zigzag nanotube ranges from 0.73 eV to 0.76 eV without any clear trend with respect to diameter within numerical accuracy. On the other hand, the binding energy of a single bromophenyl molecule shows a clear diameter dependence and ranges from 1.51 eV for a (8,0) zigzag nanotube to 0.83 eV for a (20,0) zigzag nanotube. This is in part explained by the transition from sp2 to sp3 bonding occurring to a carbon atom of a nanotube when a phenyl is grafted to it and the fact that smaller nanotubes are closer to a sp3 hybridization than larger ones due to increased curvature. Since a second bromophenyl unit can attach without energy barrier next to an isolated grafted unit, they are assumed to exist in pairs. The para configuration is found to be favored for the pairs and their binding energy decreases with increasing diameter, ranging from 4.34 eV for a (7,0) nanotube to 2.27 eV for a (29,0) nanotube. An analytic form for this radius dependence is derived using a tight binding hamiltonian and first order perturbation theory. The 1/R^2 dependance obtained (where R is the nanotube radius) is verified by our DFT results within numerical accuracy. Finally, metallic nanotubes are found to be more reactive than semiconducting nanotubes, a feature that can be explained by a non-zero density of states at the Fermi level for metallic nanotubes.

Calculating dispersion interactions using maximally-localized Wannier functions

We investigate a recently developed approach [P. L. Silvestrelli, Phys. Rev. Lett. 100, 053002 (2008); J. Phys. Chem. A 113, 5224 (2009)] that uses maximally localized Wannier functions to evaluate the van der Waals contribution to the total energy of a system calculated with density-functional theory. We test it on a set of atomic and molecular dimers of increasing complexity (argon, methane, ethene, benzene, phthalocyanine, and copper phthalocyanine) and demonstrate that the method, as originally proposed, has a number of shortcomings that hamper its predictive power. In order to overcome these problems, we have developed and implemented a number of improvements to the method and show that these modifications give rise to calculated binding energies and equilibrium geometries that are in closer agreement to results of quantum-chemical coupled-cluster calculations.

Importance of many body effects in the kernel of hemoglobin for ligand binding

We propose a mechanism for binding of diatomic ligands to heme based on a dynamical orbital selection process. This scenario may be described as bonding determined by local valence fluctuations. We support this model using linear-scaling first-principles calculations, in combination with dynamical mean-field theory, applied to heme, the kernel of the hemoglobin metalloprotein central to human respiration. We find that variations in Hunds exchange coupling induce a reduction of the iron 3d density, with a concomitant increase of valence fluctuations. We discuss the comparison between our computed optical absorption spectra and experimental data, our picture accounting for the observation of optical transitions in the infrared regime, and how the Hunds coupling reduces, by a factor of five, the strong imbalance in the binding energies of heme with CO and O_2 ligands.

Vanadium dioxide : A Peierls-Mott insulator stable against disorder

Vanadium dioxide undergoes a first order metal-insulator transition at 340 K. In this work, we develop and carry out state of the art linear scaling DFT calculations refined with non-local dynamical mean-field theory. We identify a complex mechanism, a Peierls-assisted orbital selection Mott instability, which is responsible for the insulating M$_1$ phase, and furthermore survives a moderate degree of disorder.

Linear-scaling DFT+U with full local orbital optimization

We present an approach to the DFT+U method (Density Functional Theory + Hubbard model) within which the computational effort for calculation of ground state energies and forces scales linearly with system size. We employ a formulation of the Hubbard model using nonorthogonal projector functions to define the localized subspaces, and apply it to a local-orbital DFT method including in situ orbital optimization. The resulting approach thus combines linear-scaling and systematic variational convergence. We demonstrate the scaling of the method by applying it to nickel oxide nano-clusters with sizes exceeding 7,000 atoms.