Abstract
The ever increasing availability of supercomputing resources led computer-based materials science into a new era of high-throughput calculations. Recently, Pizzi et al. introduced the AiiDA framework that provides a way to automate calculations while allowing to store the full provenance of complex workflows in a database. We present the development of the AiiDA-KKR plugin that allows to perform a large number of ab initio impurity embedding calculations based on the relativistic full-potential Korringa-Kohn-Rostoker Green function method. The capabilities of the AiiDA-KKR plugin are demonstrated with the calculation of several thousand impurities embedded into the prototypical topological insulator Sb2Te3. The results are collected in the JuDiT database which we use to investigate chemical trends as well as Fermi level and layer dependence of physical properties of impurities. This includes the study of spin moments, the impurity’s tendency to form in-gap states or its effect on the charge doping of the host-crystal. These properties depend on the detailed electronic structure of the impurity embedded into the host crystal which highlights the need for ab initio calculations in order to get accurate predictions.
Similar content being viewed by others
Introduction
In recent years computer-driven materials design has become increasingly important in the field of materials science. The ever increasing availability of supercomputing resources opened up new possibilities towards data-driven condensed matter research. Apart from large collections of crystal structure information1,2,3,4, fully integrated frameworks of tools and databases have arisen that allow for high-throughput investigations using a huge amount of, mainly, density-functional-theory-based calculations5,6,7,8. Here, we present the AiiDA-KKR plugin9 which connects our full-potential relativistic Korringa-Kohn-Rostoker Green function (KKR) method10 to the AiiDA (Automated Interactive Infrastructure and Database for Computational Science) framework8,11.
The AiiDA infrastructure implements the FAIR principle12 of findable, accessible, interoperable and reusable data sharing which provides a flexible plugin-based python environment. Through a common interface, different density-functional theory codes13 can even be used in combination to exploit the strengths of different implementations and realise multi-code workflows within the same framework. The KKR method is an all-electron implementation of density functional theory that allows accurate electronic structure calculations14, the extraction of magnetic response functions15 or gives access to transport properties16,17,18. One of the advantages of the Green function formulation of the KKR method, which we focus on in this work, is the efficient treatment of defective systems (i.e., systems that contain defects and impurities) which can be very expensive to treat with wavefunction-based methods that often require very large supercells for this task. The AiiDA-KKR plugin is used to perform a large number of impurity embedding calculations into the prototypical topological insulator Sb2Te3.
Topological insulators (TIs) have been the center of attention in solid state research since their extraordinary physical properties, that lead to topologically protected surface states, have been discovered19. In the past decade the field around topological materials has evolved steadily and now aims at functionalizing TI materials by interfacing them with other states of matter. For instance, realizing the quantum anomalous Hall (QAH)20,21 insulator state or Majorana zero modes22, that might lead to topological qubits, is pursued. The former requires to combine the topological band structures of TIs with magnetism while the latter needs interfacing topological materials with superconductors. Controlling the interface and understanding the effect defects and imperfections have remains a major challenge in this field to this day.
Apart from the development of the AiiDA-KKR plugin, the outcome of this study is the JuDiT database (Jülich Database of impurities embedded into a Topological insulator) of physical properties of impurities embedded into the surface of Sb2Te3. We study their tendency for charge doping (i.e., to introduce p- or n-doping), their impurity magnetic moments and their density of states (DOS). This collection of impurity properties allows to uncover chemical trends and can help to optimize the next generation of TI-based materials in the future. In particular, we investigate the layer and Fermi level dependence of the spin moment and find that MoSb defects show a high spin moment while introducing only a small charge doping, which is more than 5 times smaller than for magnetic 3d impurities. Furthermore, we find the Mo defect to be a good candidate for future applications since no impurity resonance appears in the bulk band gap region. An impurity resonance would otherwise lead to higher scattering rates of topological surface state electrons off this defect23 and a higher probability to produce impurity bands in the gap.
The paper is structured as follows. The theoretical setting of the impurity embedding problem within the KKR method can be found in the Methods section. In the Results section, First the AiiDA-KKR package is presented where the calculation and workflow plugins, that are implemented in AiiDA-KKR, are discussed. Afterwards the developments are showcased at the example of high-throughput impurity embedding into the topological insulator material Sb2Te3. Finally, the Discussions section concludes the paper with a summary.
Results
The AiiDA-KKR plugin
The AiiDA framework8,11 is a python package that allows to provide a code agnostic interface for different ab initio codes which enables the automation of calculations. Simultaneously, inputs and outputs of complex sequences of calculations are automatically saved in a database24. This ensures the reproducibility of all results due to the stored provenance, which consist of nodes with directed connections, in the database. AiiDA implements the ADES model8 which provides a common layer of data structures that are used by the plugins of different codes13. This enables workflows that use multiple codes, which allows interoperability and exploitation of the individual strengths of different implementations.
To enable high-throughput KKR calculations with AiiDA we developed the open source AiiDA-KKR python plugin9 that provides a set of calculations and workflows (i.e., complex sequences of calculations) and some accompanying tools. A detailed and up-to-date description of the input and output structure of the individual calculations, workflows and tools of AiiDA-KKR is included in its online documentation25 where additionally examples for the usage of the plugin are given. In the following only a brief overview of the features implemented in the AiiDA-KKR plugin is given.
Each calculation plugin comes with the functionality to create code-specific input files from AiiDA objects (e.g., StructureData objects that contain the structural information of a system) and a parser that is able to parse the retrieved output files. The conversion of AiiDA Dict objects to the input file needed by the JuKKR code is facilitated with a python class called kkrparams26 that also contains methods to verify the consistency of input parameter and write the input file. This class also knows about KKR-specific features like dealing with alloys within the coherent potential approximation14. Running a KKR calculation through AiiDA creates an acyclic directed graph in the database which is shown schematically for a KkrCalculation in Fig. 1(a). It can be seen that a calculation requires a few input nodes (input parameter Dict node, a Code node, and a ParentCalculation node), which are all AiiDA objects stored in the AiiDA database. The output files are retrieved once the calculation finished and parsed to extract output parameters, that are then stored as a Dict node in the database. It it worthwhile noting that the AiiDA daemon takes care of automatically storing the resulting nodes with their directed connections to the AiiDA database. This eventually results in a complex graph as it is shown schematically in (Fig. 1b, c) for a complete impurity-embedding sequence including all steps outlined in the Methods section. The full provenance of a complex procedure of calculations is stored in the database which allows to reproduce all inputs and the intermediate steps that have been performed to arrive at the final result. In the following the different calculations provided in the AiiDA-KKR plugin9 are briefly discussed.
-
The VoronoiCalculation plugin allows to use the voronoi code of the JuKKR package10 which constructs the shape functions27,28 needed for the full-potential treatment and generates starting potentials. The graph of a VoronoiCalculation looks similar to the one shown in Fig. 1(a) except that the ParentCalculation input node is replaced by an AiiDA StructureData node that contains all structural information on the crystal (e.g., lattice constant, atom positions and kinds)11.
-
The KkrCalculation plugin provides an interface to the KKRhost code of JuKKR10 which allows to perform self-consistency (scf), density of states (DOS), bandstructure and additional postprocessing calculations (e.g., calculation of Heisenberg exchange interaction parameters15).
-
The KkrimpCalculation plugin connects AiiDA to the KKRimp code of the JuKKR package that solves the Dyson equation for impurity embedding (Eq. (3)). This calculation needs in addition to the usual inputs (input parameter (Dict) or ParentCalculation node) the host Green function in the impurity cluster region (Ghost) which is written out with a special post-processing run-mode of the KKRhost program using a KkrCalculation. The KkrimpCalculation can then be used to perform electronic structure calculations for the impurity problem.
In addition to the calculation plugins, the AiiDA-KKR package provides some workflows that automate complex sequences of voronoi, KKRhost and KKRimp calculations. The workflows contained in AiiDA-KKR have been developed in a modular way and build upon each other. This is illustrated in Fig. 1(c) where the relation of the kkr_imp_wc and kkr_imp_sub_wc workflows is shown. Internally these are AiiDA WorkChains which is indicated by the _wc ending in the names of the workflows. For an in-depth discussion of the input and output structure and their usage we refer to the online documentation25. Here we restrict our discussion to a short overview of the workflows of AiiDA-KKR.
-
The kkr_dos_wc workflow conveniently wraps around a KkrCalculation and provides the necessary inputs to perform a DOS calculation. Additionally the output is parsed and the output DOS data is stored as an array in the database which allows easy access and plotting of the output DOS.
-
The voro_start_wc workflow wraps the VoronoiCalculation and performs some additional verification of input structure and KKR-specific parameters in order to make sure the chosen starting setting is reasonable. One of the checks performed automatically within voro_start_wc makes use of the kkr_dos_wc workflow.
-
The kkr_scf_wc workflow builds upon the voro_start_wc workflow and a sophisticated series of KkrCalculations which is intended to reach convergence of a given host system reliably. This makes sure the starting setup is reasonable before the potential is pre-converged until finally convergence with higher accuracy settings is pursued.
-
The gf_writeout_wc workflow takes care of setting the necessary options for a KkrCalculation in order to write out Ghost in preparation of the impurity embedding step. For a given impurity position and screening cluster size the host Green function can be reused for several impurities that respect this embedding geometry.
-
The kkr_imp_sub_wc workflow performs the self consistency cycle of KKRimp calculations similar to the logic implemented in the kkr_scf_wc workflow. It includes features that deal with possible convergence problems automatically.
-
The kkr_imp_wc workflow combines the voro_start_wc, gf_writeout_wc and kkr_imp_sub_wc steps for the impurity problem. This allows to conveniently start with a converged host calculation and the information on the impurity (e.g., its position in the host crystal, the size of the screening cluster) which defines the problem given by Eq. (3) completely.
-
The kkr_imp_dos_wc workflow gives an easier access to calculate the DOS of an impurity embedded into a host crystal.
Apart from tools that are used internally within the calculation and workflow plugins (e.g., to prepare the real space screening cluster), the AiiDA-KKR plugin contains a plotting tool called plot_kkr. This tool takes a node identifier (the AiiDA node instance, its pk or uuid) or even a list of nodes and plots a standard, yet fully customisable, plot of the respective data. For example, a typical self-consistency workflow for an impurity calculation (kkr_imp_wc node) as input to plot_kkr will produce a plot of the convergence behavior with the scf iteration number or an instance of a kkr_imp_dos_wc workflow will produce the orbital-resolved plot of the total DOS in the impurity cluster. Examples of such plots are given in Fig. 2. The tool conveniently abstracts away the need to extract the relevant data from the nodes in the database and provides a straightforward way to create commonly used plots.
Example usage of AiiDA-KKR: Ag defect embedded into bulk Cu
To illustrate how AiiDA-KKR’s python interface facilitates complex density functional calculations we demonstrate the use of the top-level workflows kkr_scf_wc and kkr_imp_wc at the example of an Ag impurity embedded into bulk fcc Cu. Solving this problem requires setting up the starting potential, converging the bulk electronic structure, writing out the host’s Green function for the impurity embedding and finally performing the impurity embedding step as sketched in Fig. 1(b, c). The AiiDA-KKR plugin conveniently automates this complex series of tasks. The code snippet in Fig. 3 illustrates how the self-consistent calculation for the Cu bulk can be submitted to the AiiDA daemon that takes care of orchestrating the necessary sequence of VoronoiCalculation and KkrCalculations.
Once the calculation for bulk fcc Cu finishes, the impurity embedding step can be done. This is equally simple with AiiDA-KKR’s python interface as the code snippet in Fig. 4 demonstrates.
In this example we define the Ag impurity which replaces a Cu atom and include a screening cluster around the impurity of 4 Å which contains the first two shells of Cu neighbors. From the converged calculation we can investigate the impurity DOS which is shown in Fig. 5. It can be seen that embedding an Ag atom in the Cu crystal locally changes the electronic structure around the defect. This becomes evident in the appearance of a peak in the nearest neighbor Cu DOS at the position in energy where the Ag atoms has a resonance its DOS (highlighted by the black arrow). This resonance disappears already in the second Cu neighbor which shows the localized character of the impurity state. In addition, the inset in Fig. 5 visualizes the resulting database structure of such a series of calculations (bulk scf, impurity scf, impurity DOS) where each node in the graph is an entry in the AiiDA database. The number of nodes and their connections highlight the complexity of the impurity embedding task that is conveniently taken over by the AiiDA-KKR plugin.
JuDiT – A database for impurities embedded into a TI
We apply the AiiDA-KKR plugin to embed a large number of impurities into the topological insulator Sb2Te3. The resulting JuDiT database29 comes with a webinterface30 for convenient access to the included data. The following analysis shows some of the physical insights obtained through this study. Our data analysis does not aim at being comprehensive but is intended to showcase the usefulness of our application. A future data-driven study might give additional insights but is beyond the scope of this work.
For the JuDiT database we considered substitutional defects (denoted by XY for impurity X replacing host atom Y) in the first 3 quintuple layers (QL) of a 6 QL thick film of Sb2Te3. This allows to study the influence of the topological surface state, that is mainly located in the 1st QL, on impurity properties. The band structure and DOS of the clean host system are shown in Fig. 6. In order to take into consideration doping of the host material, we investigated three possible positions of the Fermi level (EF located in the valence band (VB) and in the conduction band (CB), as well as EF in the bulk band gap). The shifted positions of the Fermi level are highlighted in Fig. 6 with red (EF in VB) and green (EF in CB) lines. For the impurity embedding we neglected structural relaxations of the atoms around the impurities but included the first 3-4 shells of neighbors (containing 21–27 atoms and empty cells in the impurity cluster) within a radius of 4.8 Å around the impurities. We used the exact description of the atomic cells27,28 and the local spin density approximation31 (LSDA) for the exchange correlation functional. A cutoff for the angular expansion of \({\ell }_{\max }=3\) was chosen and corrections for the truncation error using Lloyd’s formula have been applied32. Relativistic corrections arising within the scalar-relativistic approximations as well as spin-orbit coupling have been taken into account fully self-consistently for both host and impurity calculations. This approach has been applied successfully in the past to study magnetic and non-magnetic defects in topological insulators where a good agreement between our theoretical predictions and different experimental observations was verified23,33,34,35.
We start our analysis with an overview of the contents of the JuDiT database. In total more than 2100 impurities have been embedded self-consistently into the Sb2Te3 host system. For each impurity we computed physical properties like the spin and orbital moments, the impurity’s DOS as well as the tendency to show impurity resonance in the region of the bulk band gap. Furthermore, we analyzed the charge doping introduced by the defect which we define as
where nimp (nhost) are the electron densities for the impurity (host) atom embedded into the surrounding host crystal integrated in the Voronoi cell of the atom. Here, Zimp (Zhost) is the nuclear charge of the impurity (host) atom. The expression in brackets in the rhs of Eq. (1) therefore contain the information how much charge is transferred to/from the impurity.
The results can conveniently be visualized with the JuDiT webinterface30. The plots from Fig. 6 as well as the ones of Fig. 7 have been created using the tools available there. Figure 7(a) displays the, over the different impurity sites averaged, charge doping given in electrons per impurity. The observed chemical trends partly fit the behavior of the Pauling electronegativity36. This causal relation is highlighted by a Pearson correlation coefficient of 0.5 which is found between the induced charge density and the impurity’s electronegativity (see Supplementary Note 2 and Supplementary Fig. 3 for details). The details of the bonding mechanism are however more subtle and can be quantified using ab initio data. This reflects the intricate physics of the chemical bonding which was recently classified to be metavalent for the Sb2Te3 class of materials37.
Each impurity in JuDiT also has a detail page where the complete output dictionary of the converged calculation is given and from where its provenance can be browsed. It also features a plot of the impurity DOS which is shown exemplary in Fig. 7(b) for a TcSb impurity located in the fourth Sb layer from the surface. The position of the impurity in the Sb2Te3 host crystal is visualized in Fig. 7(c). It can be seen that the d-states of the Tc atom are exchange- and crystal-field-split which results in a magnetic moment of the impurity and a resonance in the impurity DOS around the Fermi level consequently in the bulk band gap region (blue area in Fig. 7(b), see also Fig. 6). This particular defect is therefore expected to lead to strong scattering of the topological surface state electrons23 which can induce significant back-scattering since the Tc defect is a resonant scatterer and is magnetic with a spin moment of 1.86 μB.
In order to quantify the gap-filling nature of all considered defects, we define the number of states introduced by the defect in the bulk band gap region as
where \({E}_{\min }\) and \({E}_{\max }\) are the edges of the bulk band gap region (blue area in Fig. 6), Ωimp is the Voronoi cell around the impurity atom and ρ(r; E) denotes the charge density around the atom in the Voronoi cell. A high gap-filling value consequently signals that scattering off that particular impurity will be increased which could be detrimental to the desired transport properties of TI materials. This is especially the case if the impurity is magnetic and the k → −k backscattering channel reopens due to broken time reversal symmetry.
We now focus our attention to magnetic defects which are found for some transition metal impurities. These systems are interesting in the context of realizing a robust QAH phase. We start by investigating the layer and Fermi level dependence of the spin moment of 3d transition metal impurities shown in Fig. 8(a). We see a gradual increase of the magnetic moment when going from V over Cr up to Mn dopants (blue orange and green symbols, respectively) before a subsequent decrease with Fe and Co impurities (red, violet) is observed. This behavior is expected from Hund’s rule and reflects the, from V to Co, increasing filling of the d-shell (see also Supplementary Note 1 and Supplementary Fig. 1).
The details of the impurity’s electronic structure are largely determined by two factors: (i) the atomic nature of the impurity atom determining its number of electrons in the atomic configuration and (ii) the interaction with the surrounding atoms of the host crystal that affect the hybridization of the atomic states of the impurity with the host’s band structure. This effect is seen in the layer dependence of the size of the spin-moment which is highlighted for VSb and VTe defects with the blue dashed lines in Fig. 8(a). We observe a higher spin-moment for VTe compared to VSb which can be attributed to a larger charge transfer to the impurity in the Sb substitutional site compared to the Te substitutional site. The larger charge transfer to the impurity results in a higher filling of the V d-shell and therefore a higher spin moment which can also be seen for the Cr impurity. The same mechanism leads to a decrease in the spin moment that is found for CoSb defects compared to CoTe.
Investigating the dependence of the impurity spin moment on the Fermi level (different symbols in Fig. 8(a)) reveals that the details of the hybridization with the electronic band structure of the host material can be controlled via the position of the Fermi level in the host material. A shift in the host’s Fermi level can experimentally be achieved by appropriate doping with BiSb impurities38. This change in the impurity moment has been seen previously both theoretical34 as well as experimental35 and can be attributed to the competition between the impurity seeking charge neutrality and the strong change in the hybridization with the host’s electronic structure with varying position of the Fermi level due to the presence of the bulk band gap in TI materials. This leads, for instance, to decreasing (increasing) spin moments for VSb (VTe) with increasing position of the Fermi level.
Overall we observe that for V impurities the spread in the spin moment with the impurity’s surrounding (i.e., its layer dependence) is twice as large as with varying Fermi level. On the contrary, for Mn defects the change in the spin moment in different layers and with varying Fermi level is always rather small. This results from the half-filling of the Mn d orbital that make the spin moment relatively insensitive to small changes in the hybridization with the host’s electronic structure. Nevertheless, the magnetic interactions among multiple magnetic impurity atoms can, even at small changes in the impurity hybridization, be strongly affected34,35.
Figure 8(b) summarizes the charge doping Δnimp for all defects included in JuDiT. The zigzag behavior with the impurity’s core charge reflects the structure of the periodic table with its isoelectronic groups. This is verified by a Pearson correlation coefficient of 0.68 between the impurity’s group index and the induced charge doping (see Supplementary Fig. 4 for details). Each data point in Fig. 8(b) is colored by the impurity’s spin moment which shows that the maximum of the spin-moment is found for Mn and Fe impurities (yellow points). It can however be seen that the charge doping introduced by these 3d defects is fairly large which reflects the significant difference in electronegativity compared to the Sb and Te host atoms36.
Applying magnetic doping to achieve a robust QAH phase needs to fulfill some boundary conditions in order to be feasible in experiments. In order to not tune the Fermi level out of the bulk band gap by magnetic doping, the induced charge doping should be as small as possible. At the same time magnetism is the key ingredient which calls for a sizable spin moment of the impurity. Furthermore, the magnetic impurity should not show a high DOS in the bulk band gap region to reduce the appearance of unwanted impurity bands with increasing magnetic doping. If we apply these conditions of a gap filling of \({n}_{{\rm{gap}}}^{{\rm{imp}}}\, < \,0.02\ {\rm{e}}\), a charge doping of Δnimp < 0.1 e and a spin-moment of ms > 1.5 μB we find that MoSb defects meet all these criteria. Compared to the Te substitutional site, which also shows a small charge doping and high spin moment, the Sb substitution have a gap filling which is an order of magnitude smaller and could therefore be desirable. To be able to use Mo-dopants for the realization of the QAH state the magnetic exchange coupling between Mo atoms needs to be investigated further in the future. This is, however, beyond the scope of this work. The study of the exchange interactions will be especially interesting since Mo-doping of Bi2Se3 showed signatures of antiferromagnetic coupling39 which could be possibly overcome by appropriate band structure and defect engineering as it was seen for Mn and Co doping of Bi2Te334. Here, additional codoping with other defects could open another way to design TI-based materials for future applications40.
In the framework of density functional theory the effect of on-site Coulomb repulsion for localized d or f electrons can be included with the LDA+U method41. We applied this scheme for 153 transition metal defects in the JuDiT database. For the parametrization of the correlations we used the U and J values calculated from the constrained random-phase approximation42. The values for the LDA+U parametrization used in this work for the transition metal impurities are given in Supplementary Table 1. Generally, including correlations within the LDA+U method increases the exchange splitting in magnetic impurities. This is seen, for example, in the comparison of the DOS of 3d transition metal defects with and without U-corrections (see Supplementary Fig. 2). Adding correlation effects also changes other physical properties which manifests in an increase in the median value of the spin moment by 36% (see Supplementary Fig. 5) for the LDA+U calculations compared to the previously discussed LDA results. An even stronger effect is seen on the orbital moment where we find a decrease of 58% in the median value for the subset of impurities where LDA+U calculations have been included. These results indicate that local correlations can strongly affect the outcome and should be considered when comparing calculations and experimental results.
In order to review the accuracy of our calculations and estimate possible shortcomings of our approach we attempt a comparison of physical properties contained in the JuDiT database to already published experimental and theoretical results. We focus on (i) the experimentally observed charge doping in Sb2Te3 and (ii) magnetic properties reported mainly for V and Cr doped (Bi,Sb)2Te3.
In experiments, Sb2Te3 is typically found to be p-doped which is associated to the abundance of intrinsic SbTe anti site defects43,44. This trend is confirmed by results for the impurity induced charge doping in the JuDiT database which is found to be negative ( ≈ − 0.01 e per imp. for p-doped Sb2Te3). In addition, TeSb defects show the same trend which is in line with the experimental observation that n-type doping is not realized in Sb2Te3 even under Te rich growing conditions44. However, the size of the charge doping that results from our calculations is rather small which might indicate that in strongly doped materials collective effects can lead to a more pronounced effect. This can, for instance, result from the long-ranged Coulomb interaction among charged defects that leads to a shift in the Fermi level of the whole crystal45.
In the literature doping of Sb2Te3 is mostly studied in the context of possible realizations of the QAH effect with magnetic dopants. In Table 1 we collected some experimental (e.g., from XMCD data) and theoretical (e.g., DFT supercell calculations) results for the spin moment of transition metal doped Sb2Te3 compounds. We focus on the LDA results since most calculations in the literature do not report values including correlation effects.
We can see that there is a considerable spread in the reported values for the spin moment of different impurities. Taking this methodological variance into account we find a reasonable agreement with our data from the JuDiT database. We suspect that the slight overestimation of the spin moment for V and Cr is a result of the neglected structural relaxations around the defects. The appearance of a resonance in the impurity DOS around the Fermi level (see Supplementary Figs. 1 and 2) could be removed by structural relaxations. In the spirit of the Jahn-Teller effect this could change the d-filling of the impurity and consequently lead to smaller spin moments for V and Cr defects. A future study focusing on this effect might give more insights into the effect of structural relaxations. Some experiments35,46,47 additionally report on the measured impurity DOS for V and Cr defects from resonant photoemission spectroscopy or scanning tunneling spectroscopy. A comparison to our calculated impurity DOS spectra shows good agreement for defects at the substitutional Sb site which is in line with the reasonable agreement of the spin moment reported in Table 1.
Discussion
In conclusion, we have developed the AiiDA-KKR plugin which is an open source python package that connects the JuKKR code family to the AiiDA framework. This allows to perform Korringa-Kohn-Rostoker Green function calculations in an automated high-throughput manner. We concentrated on the ability to perform ab initio impurity embedding into the topological insulator Sb2Te3.
We considered several thousand different impurities embedded into the different layers of the Sb2Te3 host crystal. This procedure allowed us to study the layer and Fermi level dependence of physical properties of defects. Specifically, we studied the chemical trends in terms of the impurity’s charge doping, their tendency to create resonances in the bulk band gap and their magnetic properties. The results have been collected in the JuDiT database which is openly available and comes with online tools for data visualization and export. Throughout our analysis we have seen that the details of the electronic structure of an impurity embedded into a host crystal is very important. The hybridization of the impurity states with its surrounding plays a crucial role for its physical properties. This highlights the relevance and the need for our ab initio calculations which provide predictive power.
In the future the AiiDA-KKR plugin in general and the resulting data of this study in particular can be used in broader high-throughput studies for quantum materials. The capabilities of AiiDA-KKR could be extended to further include the automated calculation of scattering16,17,48 and transport properties18,49 or to investigate magnetic exchange interactions15. Especially the vast space of combinations codoping with other impurities introduces will be of interest in order to find ways to tune physical properties and engineer the behavior of TI-based materials.
METHODS
Ab initio impurity embedding
One of the advantages that arises from the Green function formulation of the KKR method lies in its ability to include impurities efficiently into crystalline solids14,50. This is achieved making use of the Dyson equation
where Ghost is the Green function of the crystalline host system, ΔV = Vimp − Vhost is the difference in the potential introduced due to the presence of the impurity and Gimp is the Green function that describes the impurity embedded into the periodic host crystal. It is important to mention that the change in the potential ΔV occurs only in a small region around the impurity which is why the Dyson equation can be solved in a small real space region around the impurity site. This impurity cluster contains a few neighboring shells of host atoms that are necessary to properly treat the charge screening of the impurity by the neighboring host atoms. It is worthwhile noting that Gimp contains the complete information on physical properties like the DOS which is computed as \(\rho ({\boldsymbol{r}};E)=-\frac{1}{\pi }{\rm{Im}}{\rm{Tr}}{G}^{{\rm{imp}}}({\boldsymbol{r}},{\boldsymbol{r}};E)\) (the trace is implied over spin-, atom- and orbital momentum degrees of freedom of the Green function). This impurity embedding scheme assumes a single impurity embedded into the infinite host crystal and therefore locally breaks the translational invariance. Such a calculation is typically done for the dilute limit of defects where the Fermi level is assumed to be fixed by the host crystal. The embedded defect is then allowed to transfer charge to and from the surrounding atoms of the host crystal that are included in the impurity cluster. The collective effect of changing Fermi level can then be included by using the host Green function computed for shifted Fermi level which can affect the impurity’s charge state and may strongly affect the crystal’s overall charge doping45.
In summary, an impurity embedding calculation in the KKR formalism therefore consists of (i) calculating Ghost (from a converged host calculation), (ii) creating ΔV, (iii) performing a self-consistent field cycle (scf) to converge Gimp.
Data availability
The data generated and analysed during the current study are available in the materialscloud archive29,51. The JuDiT web interface30 is published at https://www.materialscloud.org/discover. The figures of this work can be reproduced with the tools developed for the JuDiT web interface30.
Code availability
The sourcecode of the AiiDA-KKR plugin9 is published as open source software under the MIT license at https://github.com/JuDFTteam/aiida-kkr. The source code of the JuDiT web interface30 with its visualization capabilities is open source under the MIT license as well and can be found at https://github.com/PhilippRue/JuDiT-discover-section.
References
Villars, P., Cenzual, K., Gladyshevskii, R. & Iwata, S. PAULING FILE - towards a holistic view. Chem. Met. Alloys 11, 43–76 (2018).
Hellenbrandt, M. The inorganic crystal structure database (ICSD)-present and future. Crystallography Rev. 10, 17–22 (2014).
Grazulis, S. et al. Crystallography Open Database - an open-access collection ofcrystal structures. J. Appl. Cryst. 42, 726–729 (2009).
Draxl, C. & Scheffler, M. NOMAD: The FAIR concept for big data-driven materials science. MRS Bulletin 43, 676–682 (2018).
Larsen, A. H. et al. The atomic simulation environment-A python library for working with atoms. J. Phys.: Condens. Matter 29, 273002 (2017).
Curtarolo, S. et al. The high-throughput highway to computational materials design. Nature Mater. 12, 191–201 (2013).
Jain, A. et al. Commentary: the Materials Project: a materials genome approach to accelerating materials innovation. APL Materials 1, 011002 (2013).
Pizzi, G., Cepellotti, A., Sabatini, R., Marzari, N. & Kozinsky, B. AiiDA: Automated Interactive Infrastructure and Database for Computational Science. Comp. Mat. Sci. 111, 218 (2016).
The AiiDA-KKR developers. The AiiDA-KKR plugin for the JuKKR codes. https://github.com/JuDFTteam/aiida-kkr; https://doi.org/10.5281/zenodo.3628251 (2020).
The JuKKR developers. The Jülich KKR Codes. https://jukkr.fz-juelich.de (2020).
The AiiDA developers. AiiDA: Automated Interactive Infrastructure and Database for Computational Science. http://www.aiida.net.
Wilkinson, M. D. et al. The FAIR Guiding Principles for scientific data management and stewardship. Sci Data 3, 160018 (2016).
The AiiDA plugin registry. https://aiidateam.github.io/aiida-registry/.
Ebert, H., Ködderitzsch, D. & Minár, J. Calculating condensed matter properties using the KKR-Green’s function method-recent developments and applications. Rep. Prog. Phys. 74, 096501 (2011).
Liechtenstein, A., Katsnelson, M., Antropov, V. & Gubanov, V. Local spin density functional approach to the theory of exchange interactions in ferromagnetic metals and alloys. J. Magn. Magn. Mater. 67, 65 (1987).
Heers, S. Effect of spin-orbit scattering on transport properties of low-dimensional dilute alloys. Ph.D. thesis, RWTH Aachen University (2011).
Long, N. H. et al. Spin relaxation and spin Hall transport in 5d transition-metal ultrathin film. Phys. Rev. B 90, 064406 (2014).
Zimmermann, B. et al. Fermi surfaces, spin-mixing parameter, and colossal anisotropy of spin relaxation in transition metals from ab initio theory. Phys. Rev. B 93, 144403 (2016).
Hasan, M. Z. & Kane, C. L. Colloquium: topological insulators. Rev. Mod. Phys. 82, 3045 (2010).
Yu, R. et al. Quantized anomalous hall effect in magnetic topological insulators. Science 329, 61–64 (2010).
Cui-Zu, C. et al. Experimental observation of the quantum anomalous hall effect in a magnetic topological insulator. Science 340, 6129 (2013).
SankarDas, S., Freedman, M. & Nayak, C. Majorana zero modes and topological quantum computation. npj Quantum Inf 1, 15001 (2015).
Rüßmann, P., Mavropoulos, P. & Blügel, S. Lifetime and surface-to-bulk scattering off vacancies of the topological surface state in the three-dimensional strong topological insulators Bi2Te3 and Bi2Se3. Journal of Physics and Chemistry of Solids https://doi.org/10.1016/j.jpcs.2017.12.009 (2017).
The PostgreSQL Global Development Group. PostgreSQL: The world’s most advanced open source database. http://www.postgresql.org (2020).
The AiiDA-KKR developers. Documentation of the aiida-kkr plugin. https://aiida-kkr.readthedocs.io (2020).
The Masci-tools developers. Tools, utility, parsers useful in daily material science work. https://github.com/judftteam/aiida-kkr (2020).
Stefanou, N., Akai, H. & Zeller, R. An efficient numerical method to calculate shape truncation functions for Wigner-Seitz atomic polyhedra. Comput. Phys. Commun. 60, 231 (1990).
Stefanou, N. & Zeller, R. Calculation of shape-truncation functions for Voronoi polyhedra. J. Phys.: Cond. Matter 3, 7599 (1991).
Rüßmann, P., Bertoldo, F. & Blügel, S. The JuDiT database of impurities embedded into a Topological Insulator. Materials Cloud Archive 2020.0030/v1 https://doi.org/10.24435/materialscloud:nq-ht (2020).
Rüßmann, P., Bertoldo, F. & Blügel, S. JuDiT Database Discover section. https://pgi-jcns.fz-juelich.de/judit (2020).
Vosko, S. H., Wilk, L. & Nusair, M. Accurate spin-dependent electron liquid correlation energies for local spin density calculations: a critical analysis. Can. J. Phys. 58, 1200 (1980).
Zeller, R. An elementary derivation of Lloyd’s formula valid for full-potential multiple-scattering theory. J. Phys. Cond. Matt. 16, 6453 (2004).
Sessi, P. et al. Superparamagnetism-induced mesoscopic electron focusing in topological insulators. Phys. Rev. B 94, 075137 (2016).
Rüßmann, P. et al. Towards microscopic control of the magnetic exchange coupling at the surface of a topological insulator. J. Phys. Mater. 1, 015002 (2018).
Peixoto, T. R. F. et al. Non-local effect of impurity states on the exchange coupling mechanism in magnetic topological insulators. npj Quantum Mater. 5, 87 (2020).
Lide, D. R. Table of electronegativity on the Pauling scale. CRC Handbook of Chemistry and Physics 90 (2010).
Yu, Y., Cagnoni, M., Cojocaru-Mirédin, O. & Wuttig, M. Chalcogenide thermoelectrics empowered by an unconventional bonding mechanism. Adv. Funct. Mater. 1904862 (2019).
Kellner, J. et al. Tuning the Dirac point to the Fermi level in the ternary topological insulator (Bi1−xSbx)2Te3. Appl. Phys. Lett. 107, 251603 (2015).
Deng, B., Liu, F. & Zhu, J. Absence of quantum anomalous Hall state in 4d transition-metal-doped Bi2Se3: An ab initio study. Phys. Rev. B 96, 174404 (2017).
Kim, J., Jhi, S.-H., MacDonald, A. H. & Wu, R. Ordering mechanism and quantum anomalous Hall effect of magnetically doped topological insulators. Phys. Rev. B 96, 140410(R) (2017).
Ebert, H., Perlov, A. & Mankovsky, S. Incorporation of the rotationally invariant LDA + U scheme into the SPR-KKR formalism: application to disordered alloys. Solid State Commun. 127, 443–446 (2003).
Sasioglu, E., Friedrich, C. & Blügel, S. Effective Coulomb interaction in transition metals from constrained random-phase approximation. Phys. Rev. B 83, 121101(R) (2011).
Lostak, P., Stary, Z., Horak, J. & Pancir, J. Substitutional defects in Sb2Te3 crystals. Phys. Stat. Sol. (a) 115, 87 (1989).
Cava, R. J., Ji, H., Fuccillo, M. K., Gibson, Q. D. & Ho, Y. S. Crystal structure and chemistry of topological insulators. J. Mater. Chem. C 1, 3176 (2013).
Van de Walle, C. G. & Neugebauer, J. First-principles calculations for defects and impurities: Applications to III-nitrides. J. Appl. Phys. 95, 3851 (2004).
Peixoto, T. R. F. et al. Impurity states in the magnetic topological insulator V:(Bi,Sb)2Te3. Phys. Rev. B 94, 195140 (2016).
Islam, M. F. et al. Systematics of electronic and magnetic properties in the transition metal doped Sb2Te3 quantum anomalous Hall platform. Phys. Rev. B 97, 155429 (2018).
Rüßmann, P., Mavropoulos, P. & Blügel, S. Ab-initio theory of fourier-transformed quasiparticle interference maps and application to the topological insulator Bi2Te3. phys. stat. sol. (b), in press, https://doi.org/10.1002/pssb.202000031 (2020).
Kosma, A., Rüßmann, P., Blügel, S. & Mavropoulos, P. Strong Spin-Orbit Torque effect on magnetic defects due to topological surface state electrons in Bi2Te3. Phys. Rev. B 102, 144424 (2020).
Bauer, D. S. G. Development of a relativistic full-potential first-principles multiple scattering Green function method applied to complex magnetic textures of nano structures at surfaces. Ph.D. thesis, RWTH Aachen University (2013).
Talirz, L. et al. Materials Cloud, a platform for open computational science. Sci. Data 7, 299 (2020).
Zhang, J.-M. et al. Stability, electronic, and magnetic properties of the magnetically doped topological insulators Bi2Se3, Bi2Te3, and Sb2Te3. Phys. Rev. B 88, 235131 (2013).
Tcakaev, A. et al. Comparing magnetic ground-state properties of the V- and Cr-doped topological insulator (Bi,Sb)2Te3. Phys. Rev. B 101, 045127 (2020).
Acknowledgements
We acknowledge funding from the Priority Programme SPP-1666 Topological Insulators of the Deutsche Forschungsgemeinschaft (DFG) (project MA4637/3-1), from the VITI Programme of the Helmholtz Association, by the MaX Center of Excellence funded by the EU through the H2020-EINFRA-2015-1 project: GA 676598 as well as support by the Joint Lab Virtual Materials Design (JLVMD). PR and SB acknowledge support by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy - Cluster of Excellence Matter and Light for Quantum Computing (ML4Q) EXC 2004/1 - 390534769. This work was supported by computing time granted by the JARA Vergabegremium and provided on the JARA Partition part of the supercomputer CLAIX at RWTH Aachen University. P.R. would like to thank P. Mavropoulos, J. Bröder, and the AiiDA-team for fruitful discussions.
Funding
Open Access funding enabled and organized by Projekt DEAL.
Author information
Authors and Affiliations
Contributions
P.R. and S.B. conceived the project. P.R. and F.B. developed the AiiDA-KKR plugin. P.R. performed the calculations, developed the JuDiT webinterface for the data analysis, and wrote the initial manuscript. All authors contributed to the final manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Rüßmann, P., Bertoldo, F. & Blügel, S. The AiiDA-KKR plugin and its application to high-throughput impurity embedding into a topological insulator. npj Comput Mater 7, 13 (2021). https://doi.org/10.1038/s41524-020-00482-5
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41524-020-00482-5
This article is cited by
-
Shared metadata for data-centric materials science
Scientific Data (2023)
-
Skyrmionic spin structures in layered Fe5GeTe2 up to room temperature
Communications Physics (2022)