Exploring Doping Characteristics of Various Adatoms on Single-Layer

Mar 23, 2017 - based on density functional theory to investigate the doping .... Ernzerhof formulation of generalized gradient approximation.43 ..... ...
0 downloads 0 Views 6MB Size
Subscriber access provided by University of Newcastle, Australia

Article

Exploring Doping Characteristics of Various Adatoms on Single-layer Stanene Syeda Rabab Naqvi, Tanveer Hussain, Wei Luo, and Rajeev Ahuja J. Phys. Chem. C, Just Accepted Manuscript • Publication Date (Web): 23 Mar 2017 Downloaded from http://pubs.acs.org on March 23, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry C is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 43

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Exploring Doping Characteristics of Various Adatoms on Single-layer Stanene S. R. Naqvi,1 T. Hussain*,3 W. Luo1 and R. Ahuja1,2 1

Condensed Matter Theory Group, Department of Physics and Astronomy,

Box 516, Uppsala University, S-75120 Uppsala, Sweden 2

Applied Materials Physics, Department of Materials and Engineering,

Royal Institute of Technology (KTH), S-100 44 Stockholm, Sweden 3

Centre for Theoretical and Computational Molecular Science,

Australian Institute for Bioengineering and Nanotechnology, The University of Queensland, Brisbane, Qld 4072, Australia [email protected] +61 7 33463976

 ABSRACT We have performed first-principles calculations based on density functional theory to investigate the doping characteristics of 31 different adatoms on stanene monolayer, which includes the elements of alkali metals (AM), alkaline earth metals (AEM), transition metals (TMs) and group III-VII. The most stable configurations of all the dopants have been explored by calculating and comparing binding energies of all the possible binding sites. In order to comment on the uniform distribution of adatoms on stanene, the adsorption energies (Eads) of adatoms have been compared with their experimental cohesive energies (Ec) in bulk phase. A further comparison reveals that the binding energies of the most of studied adatoms on stanene are much stronger than other group IV monolayers. Apart from structural and binding characteristics, bond lengths, adatom-adatom distance, charge transfer mechanism, electronic properties and work function has also been explored in pristine and doped monolayers. The strong adsorption of adatoms on stanene, tunable electronic properties and formation of dumbbell structures in the case of AEM and TM shows that doped stanene sheets are worth further exploration.

ACS Paragon Plus Environment

1

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 43

 INTRODUCTION The discovery of free-standing monolayer of carbon atoms, graphene,1 was followed by a sparkling upsurge of interest in condensed-matter nature of its pristine and functionalized forms for future applications.2 This stable sp2 hybridized -bonded two-dimensional (2D) network exhibited a fascinating electronic behavior. For example, linearly dispersed energies close to the Fermi level ( ) endowed Dirac-type nature of quasiparticles,3 and excellent transportation at room temperature.4, 5 However, zero band gap was a major hindrance for practical implementation of graphene in nanoscale electronics/optoelectronic tools such as transistors and photovoltaic cells.6-8 To tackle this issue, researchers explored the potential of other 2D materials possessing fascinating properties. Presently, 2D honeycomb systems of group IV elements have attracted an extensive interest due to their structural stability and distinctive nature of bonding, examples include, silicene,9-12, 33, 38, 58, 59 germanene13-17, 58, 59 and stanene.18-21, 37, 56-59 Stanene, a 2D honeycomb layer of Sn atoms is structurally analogous to silicene and germanene with a comparatively larger degree of buckling.22 Higher buckling parameter weakens the  −  orbitals interaction between Sn atoms and enhances an overlay of  and  orbitals that consequently lead to a higher dynamic stability of Sn monolayer.23 The inclusion of spin-orbit coupling confirmed the existence of Quantum Anomalous Hall (QAH) effect in monolayer stanene at room temperature.24 Additionally, the coexistence of Dirac-type behavior and a small band gap (72 meV)22 suggested a possibility of stanene as a potential candidate for nano-sized optoelectronic devices. Studies revealed that functionalization with foreign adatoms provided several possible routes to tune the electronic and magnetic behavior of pristine 2D monolayers.25-28,30-32,51 Functionalization of graphene with AM and AEM lead to the transfer of electrons from metal species to graphene.25 Pang et al. reported the successful band gap engineering of monolayer germanene via AM adsorption.27 TMs functionalized graphene sheets were capable of serving as a catalyst during particular chemical activities.34,35 Lin et al. investigated the functionalization of silicene with 15 different metal adatoms and found an improved binding energies of elements compared to graphene.39 Recently, Kadioglo et al. reported that electronic and magnetic properties of stanene monolayer were affected by AM) and AEM functionalization.40 Subsequently, Xiong et al. reported the tunable magnetic behavior of TMs embedded stanene monolayers.51 The unique electronic structure of stanene monolayer with

ACS Paragon Plus Environment

2

Page 3 of 43

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

enhanced sp3-bonding nature and higher dynamic stability distinguished it from the other group IV monolayers. Moreover, the intriguing properties upon AM and AEM metal doping provided an overwhelming urge to explore the bonding and charge transfer characteristics of functionalized stanene monolayers with other metal and non-metal adatoms. In the present work, we have performed a detailed theoretical investigation on the structural, energetics and electronic properties upon the adsorption of 31 different adatoms including AM, AEM and TM atoms, metals and non-metals of group III-VII on stanene monolayer. Most of the studied dopants not only form uniformed distribution over the monolayer but also modified the charge concentration and consequently the electronic properties of stanene monolayer.

 METHODOLOGY Electronic structure investigations were accomplished in the framework of spin-polarized density functional theory (DFT) using VASP code.41,42 Exchange and correlation energies of electrons were considered according to Perdew-Burke-Ernzerhof formulation of generalized gradient approximation.43 GGA is expected to underestimate binding energies of adatoms over stanene, thus a semi−empirical van der Waals corrections of Grimme (DFT-D2) were included throughout the calculations. Ion−electron interactions were examined according to projector-augmented wave (PAW)44 approach with an energy cut-off set to 500 eV. Moreover, the calculations performed were not spin-polarized. The reason for not including spin-polarization was that, we did not find considerable difference in total energies while preforming test calculations for pure and doped system. However for obtaining the density of states for specific systems, which resulted into highest binding, we performed spin-polarized calculations. The sampling of Brillouin zone (BZ) was acquired with Monkhorst-pack using a 3× 3× 1 mesh for structural optimization and a denser 7× 7× 1 mesh for the density of states calculations according to tetrahedron method.45 In order to avoid interactions between recurring periodic images of the surface, a vacuum thickness of 15 Å was inserted perpendicular to the plane of the sheet. For structural optimization, total energies were allowed to converge until the difference between energy of two consecutive self-consistency steps was less than 10-5 eV. The component of Hellmann-Feynman force on each individual atom was less than 0.01 eV/Å for ground state configuration. Bader charge analysis46 was

ACS Paragon Plus Environment

3

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 43

performed to understand the charge-transfer phenomena and nature of adatom−host binding. The optimized 4× 4× 1 supercell of buckled stanene monolayer is shown in Figure 1.

Figure 1. (a) Top and side views of optimized stanene monolayer (4 × 4). The dotted rhombus indicates the 1 × 1 unit cell. Buckling parameter (δ) is the difference of heights between sub-lattices A and B. Green balls represent Sn atoms.

We considered four symmetry sites for adsorption of adatoms, those include H (hollow of hexagon), T (Top, Sn atom in the upper sub-lattice), V (Valley, Sn atom in the lower sublattice) and B (the bridge between two Sn atoms), as indicated in Figure 2. Adsorption energies at each site were calculated using the following relation, E   E   −  E  E   

(1)

Here, the first term denotes the total energy of the functionalized sheet and the second term in the parenthesis represents a sum of individual energies of pristine sheet and the adatom. The adsorption of adatoms on stanene can be used as a tool to effectively tune stanene work function. The work function of adatom on stanene monolayer is calculated as ϕ  Vϕ − E

(2)

Here ϕ, Vϕ, and E denote work function, electrostatic vacuum potential and Fermi energy of the functionalized stanene sheet, respectively. The surface charge density differences δρ describe the changes in electron density during adatom adsorption to the sheet. Charge densities were acquired using the relation "ρ = ρ (sheet@adatom) – ρ (sheet)- ρ (adatom)

(3)

 RESULTS AND DISCUSSIONS We report a detailed investigation on the adsorption geometries, binding energies, chargetransfer mechanism, work-function and electronic structure of stanene monolayer functionalized with AM (Li−K), AEM (Be−Ca), TM (Sc−Zn) and group III−VII atoms such as metals (Al, Ga), metalloids (B, Si, Ge, As) and non-metals (C, N, O, P, S, Se, F, Cl, Br).

ACS Paragon Plus Environment

4

Page 5 of 43

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

For relaxed 1 × 1 unit-cell, the calculated lattice constants 4.62 Å are in agreement with reported values.40, 47 For adatom adsorption, we optimized a 4× 4 supercell (shown in Figure 1) that was larger enough to avoid the unwanted effects of long-range electrostatic interactions between the recurring images of adatoms.40 The obtained values of lattice parameters, Sn−Sn bond length and buckling parameter for supercell were 18.48 Å, 2.81 Å, and 0.883 Å, respectively. All geometrical parameters showed a reasonable agreement with the previously reported values.40,47-49, 59 Bader charge investigation depicted that each Sn atom in pristine monolayer carries a charge equivalent to +0.01é.

Figure 2. (Left) Possible adsorption sites for adatoms on stanene monolayer are indicated as B (bridge), H (hollow), T (Top, above Sn atom in sub-lattice A) and V (Valley, above Sn atom in sub-lattice B). (Right) Spin-polarized total DOS and the partial DOS (multiplied by a constant) for stanene monolayer.

ACS Paragon Plus Environment

5

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 43

Table 1. Energy and structural parameters for the stanene sheet with different adatoms. Most stable adsorption sites are indicated as H (Hollow), B (Bridge), T (Top), and V (Valley) site. The adsorption energy (eV) for the most stable site, adatom-Sn distance (Å), distortion in the monolayer stanene δSn

(Å), Bader charge on adatom (e), and the work function (eV) are listed. For comparison with previously reported data,41 we listed adsorption energy (eV) and dadatom-Sn Å) below.

Atom Site

Eads

Eads/Ecoh dadatom-Sn

(eV)

Li Na K Be Mg Ca Sc Ti V Cr Mn Fe Co Ni Cu Zn B Al Ga N P As O S Se F Cl Br

H H H V H H H H H H H H V V H V V B H V V T B B B T T T

-3.220 -2.643 -2.694 -3.152 -1.756 -3.351 -4.938 -5.277 -4.214 -3.019 -3.075 -3.759 -4.281 -4.839 -3.282 -1.162 -5.41 -3.39 -3.24 -4.12 -3.67 -3.68 -5.70 -4.27 -3.85 -4.88 -3.49 -3.06

(Å

1.98 2.37 2.88 0.94 1.16 1.82 1.27 1.09 0.79 0.73 1.05 0.88 0.97 1.09 0.94 0.86 0.93 1.0 1.15 0.84 1.08 1.24 2.19 1.49 1.56 0.84 2.49 2.51

2.8 3.13 3.5 2.4 2.96 3.12 2.85 2.77 2.76 2.79 2.75 2.68 2.37 2.37 2.71 2.79 2.26 2.76 3.15 2.13 2.59 2.74 2.01 2.44 2.58 2.00 2.44 2.83

Eads41

dadatom-

(eV)

Sn

-2.687 -2.121 -2.303 -3.058 -1.330 -2.769 ---------------------------------------------------------------------------------------------------------------

41

(Å

2.79 3.12 3.45 2.39 2.98 3.14 ---------------------------------------------------------------------------------------------------------------

δSn

Bader

10-1(Å Charge (e)

0.0 0.2 0.4 0.2 0.2 0.2 0.0 0.1 0.1 0.0 0.0 0.2 0.2 0.2 0.0 0.1 0.3 0.1 0.3 0.3 0.0 0.3 0.0 0.2 0.2 0.0 0.1 0.3

+0.851 +0.801 +0.808 +1.065 +0.902 +1.314 +1.298 +0.976 +0.666 +0.835 +0.435 +0.122 -0.508 -0.555 -0.011 +0.070 -0.781 +0.684 +0.224 -1.380 -0.664 -0.508 -1.164 -0.803 -0.659 -0.725 -0.549 -0.524

Work Function (eV)

4.21 4.14 3.89 4.42 4.35 4.00 4.08 4.22 4.21 5.14 4.10 4.14 4.39 4.40 4.37 4.44 4.36 4.36 4.28 4.19 4.39 4.34 4.47 4.46 4.45 4.60 4.59 4.48

AM, AEM and TM metal adatoms. For convenience we have divided the dopants into two sections, one that consists of AM, AEM, TMs and other section with the remaining dopants. In this section, we discuss the

ACS Paragon Plus Environment

6

Page 7 of 43

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

adsorption of AM, AEM and TM metals on pure stanene sheet. Adsorption energies were calculated for each adatom for all the four adsorptions sites and the most preferential adsorption sites were the ones corresponding to the highest binding energy of the dopants in functionalized systems. The differences of adsorption energies among the adsorption sites, in turn, yield the diffusion energies of adatoms. For example, the energy barrier for migration of Li adatom between H and T sites is obtained as Eads(Li@H)-Eads(Li@T). Similar to the work,40 AM adatoms (Li−K) show a stronger tendency to bind to H sites, which is in accordance with adsorption of such atoms on graphene, silicene, and germanene.25,28,39 Similarly, B site is a transition configuration and AM adatom placed at B site migrates to a nearby V site. The similar behavior of AM adatoms was previously reported for silicene sheet.26 Adsorption energies of AM adatoms on V sites are lower as compared to H site but higher than T site. All the three AM adatoms follow a similar kind of adsorption trend, however, this is not true for the adsorption behavior of AEM adatoms (Be−Ca), where V site is the most preferential substitutional binding site for Be. Due to the electrostatic repulsion, Be exerts an outward thrust on the Sn atom, which is initially located at V site and finally Be dopant substitutes it. Similarly, Be atom placed at B site migrates towards V site and substitutes the Sn atom while propelling it outwards, indicating that B site is a transition site for Be atom. The similar adsorption behavior of Be on silicene was previously reported26. In the most stable configuration, Mg atom prefers to adsorb on H site. V site is the next stable site, where Mg slightly propels Sn out of the sheet and forms a dumbbell structure. These kinds of structures were formed when group IV elements were doped on silicene and germanene.40,50 It is observed that B site does not prove to be a stable adsorption configuration for Mg and Ca atoms. Adsorption energy of Ca is highest at H site compared to other adsorption sites, however Ca placed on B and V sites forms a dumbbell structure at V site similar to other AEM atoms. We have tabulated the distances between adatoms and the closest Sn atom in the sheet. For comparison, adsorption energies and adsorption distances for AM, AEM, TM and group III-VII adatoms are presented in Table 1 along with the results of Kadioglu et al.40 Distortion δSn in doped single-layer stanene is quantified as a deviation of zcoordinate of Sn atoms in doped stanene compared to their average position in pure stanene. Distortion in doped stanene for all the elements (see Table. 1) remains less than 0.05 Å.

ACS Paragon Plus Environment

7

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 43

Figure 3. Adsorption energies of (a) AM, (b) AEM, (c) and (d) TM adatoms for the corresponding adsorption sites on 4× 4× 1 stanene monolayers.

Adsorption energies of TM adatoms on all the four adsorption sites on stanene as plotted in Figure 3c are higher than those of graphene, silicene and germanene.25, 28, 39, 52-53 Like AM adatoms, TMs prefer H site for adsorption in all cases except for Co, Ni, and V dopants who preferred to sit on V site. For all TM adatoms except Zn, B is a transition site and the adatom placed at B site migrates towards V site, which is the minimum energy configuration for Co, Ni, and Zn, while the second preferred site for the other TM adatoms. Each TM adatom on the V site comes closer to the sheet and propels Sn atom out of the plane of sheet to a small distance, whilst forming a dumbbell structure. The adsorption energies (Eads) of Sc, Ti, Mn, and Ni are higher than their respective cohesive energies (Ec) in bulk phase, which depicts that these atoms are capable of forming two-dimensional layers on the surface of stanene. The binding of Ti to stanene is the strongest amongst all other investigated TM adatoms, which is consistent to its binding on graphene.25 The binding enthalpy is highly site-dependent and there is a large diffusion barrier between H and other adsorption sites (> 0.57 eV). The T site preserves the weakest binding affinity for all TM atoms. V atom, which was initially introduced at T site showed a different behavior from other TMs and substituted the Sn atom propelling it out of the plane of stanene sheet to a small distance. The binding energy of Co

ACS Paragon Plus Environment

8

Page 9 of 43

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

adatom on stanene is higher than graphene but lower than its value on silicene where the ratio Eads/Ec was greater than unity.39,

52

Adsorption energies of AM, AEM and TM adatoms

corresponding to different adsorption sites and adsorption pathways are plotted in Figure 3. Bonding nature of the stanene monolayer upon the introduction of various dopants has been investigated by means of Bader charge analysis. Due to the difference of electronegativities of Sn and the dopants considered here, there is a transfer of charge from the Sn-sheet to the dopants and vice versa. In case of AMs and AEMs, comparatively lower electronegativities cause these dopants to lose a portion of their electronic charge to the stanene sheet and acquire a partial positive state. TM with higher electron affinities compared to Sn, such as Ni, Co, and Cu gain electronic charge from the Sn-sheet. Isosurface charge densities in Figure 4(a-d) clearly show the depletion and accumulation of charge on selected dopants Li, Ca, Sc and Ti. Similarly, the elements of group III-VII would lose (electronegativities < Sn) or gain (electronegativities > Sn) electronic charge when doped to stanene sheet. Complete results of the amount of charge transfer between Sn sheet and the dopants are given in Table 1. Similarly, isosurface charge densities in Figure 6(a-d) display an accumulation of charge on selected dopants of B, N, O and F, which correspond to the highest binding energies among the dopants of group III, V, VI and VII respectively. The variation of charge concentration in the functionalized stanene monolayers is coupled with the change in their electronic properties, which has been studied by plotting total (TDOS) and the partial density of states (PDOS). For convenience, we have reported DOS plots for only those dopants, which corresponds to the highest binding energies among their respective groups. For comparison to the doped systems, the TDOS and PDOS of pristine stanene monolayer has been given in Figure 4(a-d), which shows its semi-metallic nature with the valence and conduction band dominated by Sn (p) orbitals. The Sn (p) orbitals, which in the case of pristine sheet (Figure 2) touch at the Fermi level combining the valence and conduction bands, moved towards the left of Fermi in the case of both Li and Ca doping and touch each other at -0.3 eV as shown in Figure 4(a,b). The main contribution is still from Sn (p) on either side of the Fermi with a small contribution of Li (s) at -1.6 eV on the left and 1.5 eV on right. However in the case of Ca doping, more pronounced Ca (s) contributions occurred at 0.3 eV and 1.3 eV on right of the Fermi. For TMs doping, we have reported the cases of Sc and Ti as shown in Figure 4(c,d). It is evident from DOS that in contrast to AM and AEM, TMs dopants have significant contribution with Sc (d) and Ti (d) not only overlapping with Sn (p) at the Fermi but also at -1.60 eV on the its left. Major peaks can also be seen in Figure 4(c, d) at 1.1 eV and 1.5 eV on right of the Fermi level

ACS Paragon Plus Environment

9

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 43

for Sc (d) and Ti (d) respectively. From above discussion, it can be concluded that upon functionalization all the metal dopants (AM, AEM, TMs) transform semi-metallic stanene to a pure metallic one.

ACS Paragon Plus Environment

10

Page 11 of 43

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 4. (Left) Total (TDOS) and the partial density of states (PDOS) are plotted for (a) Li, (b) Ca, (c) Sc, and (d) Ti. (Right) Isosurface charge density plots of the similar systems. The elements Li, Ca, Sc, and Ti are chosen as representative elements of the corresponding group due to the strongest adsorption energies on stanene.

Adsorption of Group III−VII atoms on Stanene. In this section, we report the adsorption characteristics of group III−VII atoms on stanene. The structural properties, adsorption energies, charge transfer, and work function are listed in Table 1. For B and Al atoms, B site is a transition site as adatoms introduced at B site migrate towards V site and substitutes the Sn atom by pushing it out of the stanene’s plane to a small distance. The V site is the most preferred for both of these dopants followed by the H site being second most preferred one. However, the diffusion barrier for B atom (-1.98 eV) is quite large as compared to Al (0.08 eV). Unlike B and Al, Ga atom favors the H site, where V and B sites come after that in adsorption preference. The Eads for Al and Ga on stanene is much higher than their values reported for graphene and silicene.29,39 Moreover, the ratio Eads/Ec ensures a potential for cluster-free adsorption of these dopants on stanene, which was previously not achievable. Previously, Ozcelik et al. reported that the adsorption of group IV adatoms C, Si, Ge, and Sn on silicene and germanene forms the dumbbell structures, however, stanene sheets become unstable upon adsorption of these adatoms.50 We performed MD calculations using Nose-Thermostat algorithm to test the stability of stanene monolayers with C, Si and Ge adatoms. Our results indicate that C, Si and Ge-doped stanene monolayers become highly distorted at temperature higher than 150 K, which is in agreement with the reference.50 Adsorption energies of group VA elements in the most stable configurations are strong enough to avoid formation of three-dimensional clusters over the stanene. Particularly for P and As, the ratio of Eads to Ec is higher than 1. For the lowest energy configurations, N and P

ACS Paragon Plus Environment

11

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 43

atoms prefer to bind at V site whereas N atom initially placed at the V site propels the Sn atom out of the stanene’s plane and substitutes it. The P atom fails to completely substitute Sn and forms a dumbbell structure at V site. However, both P and As initially placed at B site migrate towards V site. The lowest Eads of N and P adatoms is observed at T site. Unlike N and P, As prefers to adsorb at T site and forms a dumbbell. V is the second and H is the third favored site for As. Group VIA elements O, S, and Se prefer B site in their most stable configurations with V, H and T are the second, third and least favored adsorption sites respectively. The O dopant propels the Sn atom at V site and substitutes it. S atom slightly pushes the Sn atom out of the plane of stanene but fails to substitute it and ends up forming a dumbbell structure. The Eads for each of these atoms is higher than Ec, which is desirable to avoid the clustering and facilitating uniform dispersion of dopants over the monolayer. Halogens F, Cl, and Br prefer to adsorb at T site, while V, B and H are next favored sites for adsorption of each of them. The adsorption energies for F are weaker compared to Ec. However, Cl and Br adsorb with a sufficiently large Eads.

ACS Paragon Plus Environment

12

Page 13 of 43

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 5. Adsorption energies of (a) Group III (b) Group V (c) Group VI, and (d) Group VII (Halogens) adatoms for the corresponding adsorption sites on 4× 4× 1 stanene monolayers.

Charge analysis reveals that each atom in the group IIIA−VIIA gets a bulk of electronic charge from the stanene monolayer except Al and Ga. The difference of electron affinities explains the amount of charge transfer and donor/acceptor elements. Sn being more electronegative compared to Al and Ga gets electronic charge from the adatoms. Other adatoms are more electronegative than Sn, thus the adsorbed elements are in the partial negative state. To study the electronic properties of group III-VII, we have selected the cases of B, N, O and F as these dopants yielded the highest binding energies among the elements of their respective groups. Similar to the case of metals (AM, AEM, TMs) doping, even here we have plotted TDOS and PDOS for pure and doped monolayers to have a better understanding of the effect of foreign adatoms in the electronic properties of stanene in Figure 6(a-d). In the case of Bdoped stanene, there is a strong hybridization between Sn (p) and B (p) at -1.85 eV on the left of Fermi as shown in Figure 6a. A small overlap between the p orbitals of both Sn and B can also be seen right at Fermi. This overlap also occurred at 1.40 eV and 2.0 eV on the right of the Fermi level. For N-doped stanene, there is mainly Sn (p) contribution at and on the left of the Fermi with N (p) states appearing relatively deep in the valence band between -3.5 eV to -3.0 eV. However, Sn (p) and N (p) do overlap at 1.0 eV on the right of the Fermi. When it comes to O doping, several states from O (p) appear on the left of the Fermi and few of them overlap with Sn (p) between -3.0 eV and -0.90 eV as shown in Figure 6c. Nevertheless, the conduction band is mainly dominated by Sn (p) with a small contribution of O (p) at 1.1 eV. Finally, we describe the case of F-doped stanene, which also constitutes a metallic system. Similar to N and O dopants mentioned above, even here bulk of the F (p) states appear quite deep in the valence band between -3.55 eV to -2.90 eV and then two overlapping F (p) and Sn (p) states at -1.20 eV and -0.75 eV close to the Fermi.

ACS Paragon Plus Environment

13

The Journal of Physical Chemistry

Page 14 of 43

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

14

Page 15 of 43

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 6. (Left) Total (TDOS) and the partial density of states (PDOS) are plotted for (a) B, (b) N, (c) O, and (d) F. (Right) Isosurface charge density plots of the similar systems. The elements B, N, O, and F are chosen as representative metals of the corresponding group due to the strongest adsorption energies on stanene.

Adsorption of adatoms on stanene leads to the variation of charge concentration as well as surface potential, which could be depicted in terms of the variation of work function (Φ) of pure and functionalized systems.54,

55

Work function of pure stanene sheet is calculated as

4.466 eV whereas Cr and K are the dopants that resulted into the highest and the lowest values of 5.14 eV and 3.89 eV respectively. The complete results of Φ for all the functionalized systems are presented in Table 1. Adsorption energies of adatoms on the host materials typically become weaker at higher doping concentrations. We have plotted the adsorption energies of adatoms on stanene at two different coverages Θ =1:32 and 1:18 in Figure 7. Adatoms having the strongest binding on stanene sheet were chosen as representative elements of the corresponding group for the investigation of adsorption behavior at different coverages. For Θ =1:18, the adsorption energies of all the adatoms except B and N still remain higher than experimental cohesive energies.

ACS Paragon Plus Environment

15

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 43

Figure 7. Adsorption energies (eV) of adatoms on the stanene sheet for two different coverages (Θ  1: 32 and 1: 18 are plotted. Adsorption energies at higher coverage are obtained only for the representative elements of corresponding group having the strongest binding on stanene.

CONCLUSION: In summary, we have studied the electronic structure, adsorption energies, charge transfer mechanism and work function variations of 31 different adatoms on stanene sheet by means density functional theory. We have also established the possibilities to form a uniform distribution of the dopants over the monolayer and the probability of cluster formation by drawing a comparison between binding and cohesive energies of the studied dopants. Our findings have yielded that stanene monolayer exhibits strong adsorption energies for most of the adatoms, which ensured that the formation of two-dimensional layers of these adatoms on the surface of stanene is possible. It is observed that elements of AM, AEM, and TMs being less electronegativities transfer a portion of their electronic charge to the stanene sheet whereas majority of the remaining dopants being more electronegative behave opposite to the formers. The accumulation or depletion of charges modified the charge concentrations in the doped systems and hence electronic properties, which were studied by plotting the TDOS and PDOS. Interestingly the group IV elements led to a large amount of distortion in the stanene sheet, indicating the instability of these doped systems, which was consistent with the

ACS Paragon Plus Environment

16

Page 17 of 43

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

literature. Upon structural optimization AEM, TMs and some other group III-VII elements form dumbbell structures on stanene.

Thus, the captivating properties explored by

introducing various foreign adatoms on stanene sheets opened up the doors for several practical explorations.

 ACKNOWLEDGEMENT: The Swedish Research Council (VR), StandUp, Carl Tryggers Stiftelse För Vetenskaplig Forskning are acknowledged for financial support. Authors are grateful to SNIC and UPPMAX for provided computing time. T. Hussain is indebted to UQ for financial support under UQ postdoctoral fellowship scheme and the resources at NCI National Facility systems at the Australian National University through National Computational Merit Allocation Scheme supported by the Australian Government and the University of Queensland Research Computing Centre.

 REFERENCES [1] Novoselov, K. S.; Geim, A. K.; Morozov, S. V.; Jiang, D.; Zhang, Y.; Dubonos, S. V.; Grigorieva, I. V.; Firsov, A. A. Electric field effect in atomically thin carbon films. Science

2004, 306, 666-669. [2] Geim, A. K.; Novoselov, K. S. The rise of graphene. Nat. Mater. 2007, 6, 183-191. 
 [3] Castro Neto, A. H.; Guinea, F.; Peres, N. M. R.; Novoselov, K. S.; Geim, A. K. The electronic properties of graphene. Rev. Mod. Phys. 2009, 81, 109-162. 
 [4] Zhang, Y.; Tan, Y. W.; Stormer, H. L.; Kim, P. Experimental observation of the quantum hall effect and berry’s phase in graphene. Nature 2005, 438, 201-204. [5] Novoselov, K. S.; Jiang, Z.; Zhang, Y.; Morozov, S. V.; Stormer, H. L.; Zeitler, U.; Maan, J. C.; Boebinger, G. S.; Kim, P.; Geim, A. K. Room-temperature quantum hall effect in graphene. Science 2007, 315,1397.
 [6] Huang, B.; Deng, H. X.; Lee, H.; Yoon, M.; Sumpter, B. G.; Liu, F.; Smith, S. C.; Wei, S. H. Exceptional optoelectronic properties of hydrogenated bilayer silicene. Phys. Rev. X 2014, 4, 021029. 
 [7] Bui, V. Q.; Pham, T. T.; Nguyen, H. V. S.; Le, H. M. Transition metal (Fe and Cr) adsorptions on buckled and planar silicene monolayers: a density functional theory investigation. J. Phys. Chem. C 2013, 117, 23364-23371.

ACS Paragon Plus Environment

17

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 43

[8] Shkrebtii, A. I.; Heritage, E.; McNelles, P.; Cabellos, J. L.; Mendoza, B. S. Graphene and graphane functionalization with hydrogen: electronic and optical signatures. Phys. Status

Solidi C 2012, 9, 1378-1383. 
 [9] Cahangirov, S.; Topsakal, M.; Akturk, E.; Sahin, H.; Ciraci, S. Two- and one- dimensional honeycomb structures of silicon and germanium. Phys. Rev. Lett. 2009, 102, 236804. 
 [10] Lalmi, B.; Oughaddou, H.; Enriquez, H.; Kara, A.; Vizzini, S.; Ealet, B.; Aufray, B. Epitaxial growth of a silicene sheet. Appl. Phys. Lett. 2010, 97, 223109. 
 [11] Feng, B.; Ding, Z.; Meng, S.; Yao, Y.; He, X.; Cheng, P.; Chen, L.; Wu, K. Evidence of silicene in honeycomb structures of silicon on Ag (111). Nano Lett. 2012, 12, 3507-3511. 
 [12] Fleurence, A.; Friedlein, R.; Ozaki, T.; Kawai, H.; Wang, Y.; Yamada-Takamura, Y. Experimental evidence for epitaxial silicene on diboride thin films. Phys. Rev. Lett. 2012,108, 245501. 
 [13] Ozcelik, V. O.; Durgun, E.; Ciraci, S. New phases of germanene. J. Phys. Chem. Lett.

2014, 5, 2694-2699. 
 [14] Rachel, S.; Ezawa, M. Giant magnetoresistance and perfect spin filter in silicene, germanene and stanene. Phys. Rev. B 2014, 89,195303. 
 [15] Ma, Y.; Dai, Y.; Lu, Y. B.; Huang, B. Effective band gap engineering in wrinkled germanane via tiny electric field. J. Mater. Chem. C 2014, 2, 1125-1130. 
 [16] Li, L.; Lu, S.; Pan, J.; Qin, Z.; Wang, Y.; Wang, Y.; Cao, G.; Du, S.; Gao, H. Buckled germanene formation on Pt (111). Adv. Mater. 2014, 26, 4820-4824. 
 [17] Nijamudheen, A.; Bhattacharjee, R.; Choudhury, S.; Datta, A. Electronic and chemical properties of germanene: the crucial role of buckling. J. Phys. Chem. C 2015, 119, 38023809. 
 [18] Xu, Y.; Yan, B.; Zhang, H. J.; Wang, J.; Xu, G.; Tang, P.; Duan, W.; Zhang, S. C. Largegap quantum spin hall insulators in tin films. Phys. Rev. Lett. 2013, 111, 136804. 
 [19] Tang, P.; Chen, P.; Cao, W.; Huang, H.; Cahangirov, S.; Xian, L.; Xu, Y.; Zhang, S. C.; Duan, W.; Rubio, A. Stable two-dimensional dumbbell stanene: a quantum spin hall insulator. Phys. Rev. B 2014, 90, 121408. 
 [20] Ma, Y.; Dai, Y.; Guo, M.; Niu, C.; Huang, B. Intriguing behavior of halogenated twodimensional tin. J. Phys. Chem. C 2012, 116, 12977-12981. 
 [21] Rivero, P.; Yan, J. A.; Garcia-Suarez, V. M.; Ferrer, J.; Barraze-Lopez, S. Stability and properties of high-buckled two-dimensional tin and lead. Phys. Rev. B 2014, 90, 241408. 


ACS Paragon Plus Environment

18

Page 19 of 43

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

[22] Cai, B.; Zhang, S.; Hu, Z.; Hu, Y.; Zou, Y.; Zeng, H. Tinene: a two-dimensional Dirac material with a 72 meV band gap. Phys. Chem. Chem. Phys. 2015, 17, 12634-12638. [23] Chen, X.; Tan, C.; Yang, Q.; Meng, R.; Liang, Q.; Cai, M.; Zhang, S.; Jiang, J. Ab Initio study of the adsorption of small molecules on stanene. J. Phys. Chem. C 2016,120, 1398713984. [24] Wu, S. C.; Shan, G.; Yan, B. Prediction of near-room-temperature quantum anomalous hall effect on honeycomb materials. Phys. Rev. Lett. 2014, 113, 256401. [25] Chan, K. T.; Neaton, J. B.; Cohen, M. L. First-principles study of metal adatom adsorption on graphene. Phys. Rev. B 2008, 77, 235430. [26] Sahin, H.; Peeters, F. M. Adsorption of alkali, alkaline-earth and 3d transition metal atoms on silicene. Phys. Rev. B 2013, 87, 085423. [27] Pang, Q.; Zhang, C. L.; Li, L.; Fu, Z. Q.; Wei, X. M.; Song, Y. L. Adsorption of alkali metal atoms on germanene: a first-principles study. Appl. Surf. Sci. 2014, 314, 15-20. [28] Pang. Q.; Li, L.; Zhang, L. L.; Zhang, C. L.; Song, Y. L., Functionalization of germanene by metal atoms adsorption: a first-principles study. Can. J. Phys. 2015, 93, 1310-1318.
 [29] Legesse, M.; Mellouhi, F. E.; Bentria, E. T.; Madjet, M. E.; Fisher, T. S.; Kais, S.; Alharbi, F. H. Reduced work function of graphene by metal adatoms. Appl. Surf. Sci. 2017, 394, 98-107. [30] Profeta, G.; Calandra, M.; Mauri, F. Phonon-mediated superconductivity in graphene by lithium deposition. Nat. Phys. 2012, 8, 131-134. [31] Cabria, I.; Lopez, M. J.; Alonso, J. A. Enhancement of hydrogen physisorption on graphene and carbon nanotubes by Li doping. J. Chem. Phys. 2005, 123, 204721.
 [32] Ataca, C.; Akturk, E.; Ciraci, S.; Ustunel, H. High-capacity hydrogen storage by metallized graphene. Appl. Phys. Lett. 2008, 93, 043123.
 [33] Hussain, T.; Sudip, C.; Sarkar, A. D.; Johansson, B.; Ahuja, R. Enhancement of energy storage capacity of Mg functionalized silicene and silicane under external strain. Appl. Phys.

Lett. 2014, 105, 123903. [34] Lu, Y. H.; Zhou, M.; Zhang, C.; Feng, Y. P. Metal-embedded graphene: a possible catalyst with high activity. J. Phys. Chem. C 2009, 113, 20156-20160. [35] Li, Y.; Zhou, Z.; Yu, G.; Chen, W.; Chen, Z. CO catalytic oxidation on iron-embedded graphene: computational quest for low-cost nanocatalysts. J. Phys. Chem. C 2010, 114, 62506254. [36] Chapman, J.; Su, Y.; Howard, C. A.; Kundys, D.; Grigorenko, A.; Guinea, F.; Geim, A.

ACS Paragon Plus Environment

19

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 43

K.; Grigorieva, I. V.; Nair, R. R. Superconductivity in Ca-doped graphene laminates. Sci.

Rep. 2016, 6, 23254. [37] Balendhran, S.; Walia, S.; Nili, H.; Sriram, S.; Bhaskaran, M. Elemental analogues of graphene: silicene, germanene, stanene and phosphorene. Small 2015, 11, 640-652. [38] Hussain, T.; Chakraborty, S.; Ahuja, R. Metal-functionalized silicene for efficient hydrogen storage. ChemPhysChem 2013, 14, 3463-3466.
 [39] Lin, X.; Ni, J. Much stronger binding of metal adatoms to silicene than to graphene: a first-principles study. Phys. Rev. B 2012, 86, 075440. [40] Kadioglu, Y.; Ersan, F.; Gokoglu, G.; Akturk, O. U.; Akturk. E. Adsorption of alkali and alkaline-earth metal atoms on stanene: A first-principles study. Mater. Chem. Phys. 2016, 180, 326-331. [41] Kresse, G.; Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. Rev. B 1993, 47, 558-561.
 [42] Kresse, G.; Hafner, J. Ab initio molecular-dynamics simulation of the liquid-metal– amorphous-semiconductor transition in germanium. Phys. Rev. B 1994, 49, 14251-14269. [43] Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 1996, 77, 3865-3868. [44] Blochl, P. E. Projector augmented-wave method. Phys. Rev. B 1994, 50, 17953-17979.
 [45] Monkhorst, H. J.; Pack, J. D. Special points for brillouin-zone integrations. Phys. Rev. B

1976, 13, 5188-5192. [46] Bader, R. F. W., Atoms in molecules - a quantum Theory; Oxford University Press: Oxford, 1990. [47] Peng, B.; Zhang, H.; Shao, H.; Xu, Y.; Zhang, X.; Zhu, H. Low lattice thermal conductivity of stanene. Sci. Rep. 2016, 6, 20225. [48] Fang, Y.; Huang, Z. Q.; Hsu, C. H.; Li, X.; Xu, Y.; Zhou, Y.; Wu, S.; Chuang, F. C.; Zhu. Z. Z. Quantum spin hall states in stanene/Ge (111). Sci. Rep. 2015, 5, 14196. [49] Xu, Y.; Yan, B.; Zhang, H. -J.; Wang, J.; Xu, G.; Tang, P.; Duan, W.; Zhang, S. C. Large-gap quantum spin hall insulators in tin films. Phys. Rev. Lett. 2013, 111, 136804. [50] Ozcelik, V. O.; Kecik, D.; Durgun, E.; Ciraci, S. Adsorption of group IV elements on graphene, silicene, germanene, and stanene: Dumbbell formation. J. Phys. Chem. C 2014, 119, 845-853.

ACS Paragon Plus Environment

20

Page 21 of 43

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

[51] Xiong, W.; Xia, C.; Wang, T.; Du, J.; Peng, Y.; Zhao, X.; Jia, Y. Tuning electronic structures of the stanene monolayer via defects and transition-metal-embedding: spin–orbit coupling. Phys. Chem. Chem. Phys. 2016, 18, 28759-28766. [52] Cao, C.; Wu, M.; Jiang, J.; Cheng, H. P. Transition metal adatom and dimer adsorbed on graphene: Induced magnetization and electronic structures. Phys. Rev. B 2010, 81, 205424. [53] Sevincli, H.; Topsakal, M.; Durgun, E.; Ciraci, S. Electronic and magnetic properties of 3d transition-metal atom adsorbed graphene and graphene nanoribbons. Phys. Rev. B 2008, 77, 195434. [54] Hussain, T.; Panigrahi, P.; Ahuja, R. Enriching physisorption of H2S and NH3 gases on a graphane sheet by doping with Li adatoms. Phys. Chem. Chem. Phys. 2014, 16, 8100-8105. [55] Hussain, T.; Panigrahi, P.; Ahuja, R. Sensing propensity of a defected graphane sheet towards CO, H2O and NO2. Nanotechnology 2014, 25, 325501. [56] Zhu, F.; Chen, W. J; Xu, Y.; Gao, C. L.; Guan, D. D.; Liu, C.; Qian, D.; Zhang, S. C.; Jia, J. F. Epitaxial growth of two-dimensional stanene. Nat. Mater. 2015, 14, 1020-1025. [57] Takahashi, L.; Takahashi, K. Low temperature pollutant trapping and dissociation over two-dimensional tin. Phys. Chem. Chem. Phys. 2015, 17, 21394-21396. [58] Mortazavi, B.; Rahaman, O.; Makaremi, M.; Dianat, A.; Cuniberti, G.; Rabczuk, T. Firstprinciples investigation of mechanical properties of silicene, germanene and stanene. Phys. E

(Amsterdam, Neth.) 2017, 87, 228-232. [59] Mortazavi, B.; Dianat, A.; Cuniberti, G.; Rabczuk, T. Application of silicene, germanene and stanene for Na or Li ion storage: A theoretical investigation, Electrochim.

Acta 2016, 213, 865-870.

ACS Paragon Plus Environment

21

The Journal of Physical Chemistry

Page 22 of 43

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

22

Page 23 of 43

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42

The Journal of Physical Chemistry

Li Mg Sc Ga Mn S As Ni Cl

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42

ACS Paragon Plus Environment

Page 24 of 43

Page 25 of 43

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Paragon Plus Environment

Page 26 of 43

Page 27 of 43

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Paragon Plus Environment

Page 28 of 43

Page 29 of 43

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Paragon Plus Environment

Page 30 of 43

Page 31 of 43

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Paragon Plus Environment

Page 32 of 43

Page 33 of 43

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Paragon Plus Environment

Page 34 of 43

Page 35 of 43

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Paragon Plus Environment

Page 36 of 43

Page 37 of 43

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Paragon Plus Environment

Page 38 of 43

Page 39 of 43

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Paragon Plus Environment

Page 40 of 43

Page 41 of 43

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Paragon Plus Environment

Page 42 of 43

Page 43 of 43

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Paragon Plus Environment