Search

Plasmonic near-field transducers (NFTs) play a key role in administering nanoscale heating for a number of applications ranging from medical devices to next generation data processing technology. We present a novel multi-scale approach, combining quantum many-body perturbation theory with finite-element modelling, to predict the electric and thermal material parameters of various Au-based, noble metal (M) alloys. Specifically, we focus on modelling their performance within an NFT designed to focus high-intensity, sub-diffracted light for technologies such as nanoscale etching, manipulation, sensing, and heat-assisted magnetic recording (HAMR). Elemental Au is the long-standing general-purpose NFT medium due its excellent plasmonic performance at relevant wavelengths. However, elemental Au is a soft, ductile material that tends to extrude and deform in response to extreme temperature gradients. Therefore, alloying Au with other noble metals such as Ag, Cu, Pd or Pt, has attracted considerable interest for improved mechanical and thermal robustness while reaching threshold plasmonic generation at standard optoelectronics operating wavelengths (e.g., $\approx 830$~nm) and approximate high-power NFT temperatures ($\approx 400$~K). We predict that certain Au-Ag alloys may offer improved thermal stability as whole-NFT media compared to elemental Au, alongside plasmonic figures of merit comparable to that of Au. Simulations of certain solid solution Au-Pd/Pt alloys enable us to predict significantly enhanced thermal conductivity. We predict that alloying with Pd at low concentrations $\sim 10\%$ may preserve the NFT performance of Au, while offering the benefits of improved thermal and mechanical stability., Comment: 8 pages, 4 figures, 1 table. Includes Supplemental Material with 10 tables

Burgess et al. have recently introduced the BLOR corrective exchange-correlation functional that is, by construction, the unique simplified rotationally-invariant DFT+U functional that enforces the flat-plane condition separately on each effective orbital of a localized subspace. Detached from the Hubbard model, functionals of this type are both double-counting correction free and, when optimized in situ using appropriate error quantifiers, effectively parameter free. In this work, the extension of the BLOR functional to address many-body errors (mBLOR) is derived. The mBLOR functional is built to enforce the flat-plane condition on the entire subspace, rather than on each orbital individually. In this way inter-orbital errors are corrected on the same footing as the single-particle ones. Focusing on exact test cases with strong inter-orbital interactions, the BLOR and mBLOR functionals were benchmarked against contemporary DFT+U functionals using the total energy extensivity condition on stretched homo-nuclear p-block dimers that represent various self-interaction and static-correlation error regimes. The BLOR functional outperformed all other DFT+$U$ functionals tested, which often act to increase total-energy errors, yet it still yielded large errors in some systems. mBLOR instead yielded low energy errors across all four strongly-correlated dimers, while being constructed using only semi-local approximation ingredients. As mBLOR would not otherwise introduce a band-gap correction in the manner that is a desirable feature of DFT+U, we developed a cost-free technique to reintroduce it automatically by moving the functional's unusual explicit derivative discontinuity into the potential. With this in place, mBLOR is the only known DFT$+U$ functional that opens the bandgap of stretched neutral homo-nuclear dimers without the aid of unphysical spin-symmetry breaking., Comment: 21 pages, 9 figures, and 4 tables

We present the implementation of the Hubbard ($U$) and Hund ($J$) corrected Density Functional Theory (DFT+$U$+$J$) functionality in the Quickstep program, which is part of the CP2K suite. The tensorial and L\"owdin subspace representations are implemented and compared. Full analytical DFT+$U$+$J$ forces are implemented and benchmarked for the tensorial and L\"owdin representations. We also present the implementation of the recently proposed minimum-tracking linear-response method that enables the $U$ and $J$ parameters to be calculated on first principles basis without reference to the Kohn-Sham eigensystem. These implementations are benchmarked against recent results for different materials properties including DFT+$U$ band gap opening in NiO, the relative stability of various polaron distributions in TiO$_2$, the dependence of the calculated TiO$_2$ band gap on +$J$ corrections, and, finally, the role of the +$U$ and +$J$ corrections for the computed properties of a series of the hexahydrated transition metals. Our implementation provides results consistent with those already reported in the literature from comparable methods. We conclude the contribution with tests on the influence of the L\"owdin orthonormalization on the occupancies, calculated parameters, and derived properties.

Members of the DFT+U family of functionals are increasingly prevalent methods of addressing errors intrinsic to (semi-) local exchange-correlation functionals at minimum computational cost, but require their parameters U and J to be calculated in situ for a given system of interest, simulation scheme, and runtime parameters. The SCF linear response approach offers ab initio acquisition of the U and has recently been extended to compute the J analogously, which measures localized errors related to exchange-like effects. We introduce a renovated post-processor, the lrUJ utility, together with this detailed best-practices guide, to enable users of the popular, open-source Abinit first-principles simulation suite to engage easily with in situ Hubbard parameters and streamline their incorporation into material simulations of interest. Features of this utility, which may also interest users and developers of other DFT codes, include $n$-degree polynomial regression, error analysis, Python plotting facilities, didactic documentation, and avenues for further developments. In this technical introduction and guide, we place particular emphasis on the intricacies and potential pitfalls introduced by the projector augmented wave (PAW) method, SCF mixing schemes, and non-linear response, several of which are translatable to DFT+U(+J) implementations in other packages., Comment: 46 pages, 6 figures (+3 figures in appendix), 2 tables (+1 table in appendix)

The Niobium surface is almost always covered by a native oxide layer which greatly influences the performance of superconducting devices. Here we investigate the highly stable Niobium oxide overlayer of Nb(110), which is characterised by its distinctive nanocrystal structure as observed by scanning tunnelling microscopy (STM). Our ab-initio density functional theory (DFT) calculations show that a subtle reconstruction in the surface Niobium atoms gives rise to rows of 4-fold coordinated oxygen separated by regions of 3-fold coordinated oxygen. The 4-fold oxygen rows are determined to be the source of the nanocrystal pattern observed in STM, and the two chemical states of oxygen observed in core-level X-ray photoelectron spectroscopy (XPS) are ascribed to the 3-fold and 4-fold oxygens. Furthermore, we find excellent agreement between the DFT calculated electronic structure with scanning tunnelling spectroscopy and valence XPS measurements., Comment: 8 pages, 4 figures, plus 3 pages of Supporting Information

It has long been postulated that within density-functional theory (DFT) the total energy of a finite electronic system is convex with respect to electron count, so that 2 E_v[N_0] <= E_v[N_0 - 1] + E_v[N_0 + 1]. Using the infinite-separation-limit technique, this article proves the convexity condition for any formulation of DFT that is (1) exact for all v-representable densities, (2) size-consistent, and (3) translationally invariant. An analogous result is also proven for one-body reduced density matrix functional theory. While there are known DFT formulations in which the ground state is not always accessible, indicating that convexity does not hold in such cases, this proof nonetheless confirms a stringent constraint on the exact exchange-correlation functional. We also provide sufficient conditions for convexity in approximate DFT, which could aid in the development of density-functional approximations. This result lifts a standing assumption in the proof of the piecewise linearity condition with respect to electron count, which has proven central to understanding the Kohn-Sham band-gap and the exchange-correlation derivative discontinuity of DFT., Comment: This Communication has been published in the Journal of Chemical Physics 159, 211102 (2023) which can be found at the following link: https://doi.org/10.1063/5.0174159

The piecewise linearity condition on the total energy with respect to the total magnetization of finite quantum systems is derived, using the infinite-separation-limit technique. This generalizes the well-known constancy condition, related to static correlation error, in approximate density functional theory (DFT). The magnetic analog of the DFT Koopmans' theorem is also derived. Moving to fractional electron count, the tilted plane condition is derived, lifting certain assumptions in previous works. This generalization of the flat plane condition characterizes the total energy surface of a finite system for all values of electron count N and magnetization M. This result is used in combination with tabulated spectroscopic data to show the flat plane-structure of the oxygen atom, among others. We find that derivative discontinuities with respect to electron count sometimes occur at non-integer values. A diverse set of tilted plane structures is shown to occur in d-orbital subspaces, depending on chemical coordination. General occupancy-based total-energy expressions are demonstrated thereby to be necessarily dependent on the symmetry-imposed degeneracies., Comment: 13 pages, 7 figures

High-performance batteries, heterogeneous catalysts and next-generation photovoltaics often centrally involve transition metal oxides (TMOs) that undergo charge or spin-state changes. Demand for accurate DFT modeling of TMOs has increased in recent years, driving improved quantification and correction schemes for approximate DFT's characteristic errors, notably those pertaining to self-interaction and static correlation. Of considerable interest, meanwhile, is the use of DFT-accessible quantities to compute parameters of coarse-grained models such as for magnetism. To understand the interference of error corrections and model mappings, we probe the prototypical Mott-Hubbard insulator NiO, calculating its electronic structure in its antiferromagnetic I/II and ferromagnetic states. We examine the pronounced sensitivity of the first principles calculated Hubbard U and Hund's J parameters to choices concerning Projector Augmented Wave (PAW) based population analysis, we reevaluate spin quantification conventions for the Heisenberg model, and we seek to develop best practices for calculating Hubbard parameters specific to energetically meta-stable magnetic orderings of TMOs. Within this framework, we assess several corrective functionals using in situ calculated U and J parameters, e.g., DFT+U and DFT+U+J. We find that while using a straightforward workflow with minimal empiricism, the NiO Heisenberg parameter RMS error with respect to experiment was reduced to 13%, an advance upon the state-of-the-art. Methodologically, we used a linear-response implementation for calculating the Hubbard U available in the open-source plane-wave DFT code Abinit. We have extended its utility to calculate the Hund's exchange coupling J, however our findings are anticipated to be applicable to any DFT+U implementation., Comment: 19 pages, 5 figures (+2 in SI), 7 tables (+1 in SI). Accepted Nov. 14, 2023 to Physical Review B

We study the modification of the magnetocrystalline anisotropy (MCA) of Co slabs induced by several different conjugated molecular overlayers, i.e., benzene, cyclooctatetraene, naphthalene, pyrene and coronene. We perform first-principles calculations based on Density Functional Theory and the magnetic force theorem. Our results indicate that molecular adsorption tends to favour a perpendicular MCA at surfaces. A detailed analysis of various atom-resolved quantities, accompanied by an elementary model, demonstrates that the underlying physical mechanism is related to the metal-molecule interfacial hybridization and, in particular, to the chemical bonding between the molecular C $p_z$ and the out-of-plane Co $d_{z^2}$ orbitals. This effect can be estimated from the orbital magnetic moment of the surface Co atoms, a microscopic observable accessible to both theory and experiments. As such, we suggest a way to directly assess the MCA modifications at molecule-decorated surfaces, overcoming the limitations of experimental studies that rely on fits of magnetization hysteresis loops. Finally, we also study the interface between Co and both C$_{60}$ and Alq$_3$, two molecules that find widespread use in organic spintronics. We show that the modification of the surface Co MCA is similar upon adsorption of these two molecules, thereby confirming the results of recent experiments., Comment: 10 figures in main text and 3 in the SM, 20 pages

Multi-principal element materials (MPEMs) have been attracting a rapidly growing interest due to their exceptional performance under extreme conditions, from cryogenic conditions to extreme-high temperatures and pressures. Despite the simple conceptual premise behind their formation, computational high-throughput first-principles design of such materials is extremely challenging due to the large number of realizations required for sufficient statistical sampling of their design space. Furthermore, MPEMs are also known to develop short-ranged orderings (SROs) which can play a significant role in their stability and properties. Here, we present an expedient and efficient first-principles computational framework for assessing the compositional and mechanical properties of MPEMs, including SRO effects. This heuristic methodology systematically corrects phase-averaged free-energies of MPEMs to include SRO phases, while imposing constraints for materials design. To illustrate the methodology, we study the stability and mechanical properties of equi-molar refractory-metal high-entropy alloy carbonitrides (RHEA-CNs) such as ZrNbMoHfTaWC3N3. We show that SRO, arising due to preferential neighboring among refractory metals, is necessary for thermodynamic and mechanical stability and to satisfy the imposed design criteria, leading to complex compositions for which their molar fraction and mechanical properties are predicted., Comment: 7 pages, 4 figures, and 1 table. Supplemental Material of 20 pages, 10 figures, and 4 tables. 17MB

A new DFT+U type corrective functional is derived from first principles to enforce the flat plane condition on localized subspaces, thus dispensing with the need for an ad hoc derivation from the Hubbard model. The newly derived functional as given by equation 5 yields relative errors below 0.6% in the total energy of the dissociated s-block dimers as well as the dissociated H5+ ring system. In comparison bare PBE and PBE+U (using Dudarev's 1998 Hubbard functional) yields relative energetic errors as high as 8.0% and 20.5% respectively., Comment: 8 pages, 5 figures, plus 13 pages of Supplementary Information

DFT+U provides a convenient, cost-effective correction for the self-interaction error (SIE) that arises when describing correlated electronic states using conventional approximate density functional theory (DFT). The success of a DFT+U(+J) calculation hinges on the accurate determination of its Hubbard U and Hund's J parameters, and the linear response (LR) methodology has proven to be computationally effective and accurate for calculating these parameters. This study provides a high-throughput computational analysis of the U and J values for transition metal d-electron states in a representative set of over 2000 magnetic transition metal oxides (TMOs), providing a frame of reference for researchers who use DFT+U to study transition metal oxides. In order to perform this high-throughput study, an atomate workflow is developed for calculating U and J values automatically on massively parallel supercomputing architectures. To demonstrate an application of this workflow, the spin-canting magnetic structure and unit cell parameters of the multiferroic olivine LiNiPO4 are calculated using the computed Hubbard U and Hund J values for Ni-d and O-p states, and are compared with experiment. Both the Ni-d U and J corrections have a strong effect on the Ni-moment canting angle. Additionally, including a O-p U value results in a significantly improved agreement between the computed lattice parameters and experiment., Comment: 18 pages, 6 figures

First-principles Hubbard-corrected approximate density-functional theory (DFT+U) is a low-cost, potentially high throughput method of simulating materials, but it has been hampered by empiricism and inconsistent band-gap correction in transition-metal oxides. DFT+U property prediction of non-magnetic systems such as d0 and d10 transition-metal oxides is typically faced with excessively large calculated Hubbard U values, and with difficulty in obtaining acceptable band-gaps and lattice volumes. Meanwhile, Hund's exchange coupling J is an important but often neglected component of DFT+U, and the J parameter has proven challenging to directly calculate by means of linear response. In this work, we provide a revised formula for computing Hund's J using established self-consistent field DFT+U codes. For non-magnetic systems, we introduce a non-approximate technique for calculating U and J simultaneously in such codes, at no additional cost. Using unmodified Quantum ESPRESSO, we assess the resulting values using two different DFT+U functionals incorporating J, namely the widely used DFT+(U-J) and the readily available DFT+U+J. We assess a test set comprising TiO2, ZrO2, HfO2, Cu2O and ZnO, and apply the corrections both to metal and oxygen centered pseudoatomic subspaces. Starting from the PBE functional, we find that DFT+(U-J) is significantly out-performed in band-gap accuracy by DFT+U+J, the RMS band-gap error of which matches that of the hybrid functional HSE06. ZnO, a long-standing challenge case for DFT+U, is addressed by means of Zn 4s instead of Zn 3d correction, in which case the first-principles DFT+U+J band-gap error is half of that reported for HSE06., Comment: 32 pages 7 figures, 10 tables