Stability and Reactivity of Silicon Magic Numbers Doped with

Dec 4, 2018 - Gabriel F. S. Fernandes , Max Pinheiro Jr. , Francisco Bolivar Correto Machado , and Luiz Fernando de Araujo Ferrão. J. Phys. Chem. A , ...
0 downloads 0 Views 3MB Size
Subscriber access provided by YORK UNIV

A: Spectroscopy, Molecular Structure, and Quantum Chemistry

Stability and Reactivity of Silicon Magic Numbers Doped with Aluminum and Phosphorus Atoms Gabriel F. S. Fernandes, Max Pinheiro Jr., Francisco Bolivar Correto Machado, and Luiz Fernando de Araujo Ferrão J. Phys. Chem. A, Just Accepted Manuscript • DOI: 10.1021/acs.jpca.8b10214 • Publication Date (Web): 04 Dec 2018 Downloaded from http://pubs.acs.org on December 9, 2018

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 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 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.

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 37 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

Stability and Reactivity of Silicon Magic Numbers Doped with Aluminum and Phosphorus Atoms Gabriel F.S. Fernandes1, Max Pinheiro Júnior1, Francisco B. C. Machado1, Luiz F. A. Ferrão1* 1Departamento

de Química, Instituto Tecnológico de Aeronáutica, São José dos Campos, SP

12228-900, Brasil.

__________________________________ *Corresponding author. E-mail address: [email protected]

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

Page 2 of 37

2

Abstract The progressive scaling down of the silicon-based electronics has allowed to develop devices at nanometer scale, requiring new engineering techniques guided by fundamental chemical and physical concepts. Particularly, the nanostructured cluster systems are promising materials since their physical-chemical properties are sensitive to its shape, size and chemical components, such that completely different materials can be produced by the simple addition or removal of a single atom. These size-tunable properties can open a new area in materials science and engineering. In the present work, quantum chemical methods were used to study the chemical substitution effects caused by subvalent (aluminum) and supervalent (phosphorus) atoms in the physical-chemical properties of some small silicon clusters which demonstrate high stability, called magic numbers. The changes in the electronic structure and chemical acceptance to the dopants were evaluated with respect to: ionization potential, electronic excitation energy, stability and reactivity parameters. Taken together, these results enable to identify the most stable silicon-doped clusters. Regarding electrophilic reactions, Si10P is the most favorable system, while for nucleophilic reactions, none of the doped clusters resulted in higher stability.

ACS Paragon Plus Environment

Page 3 of 37 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

3

I. Introduction The tremendous increasing rate of electronic device miniaturization over the past decade, pushing the limits of the Moore’s law1–4, has enabled the industrial production of silicon transistors with tens of nanometers using top-down techniques, such as extreme ultraviolet lithography5,6. In laboratorial scale, 1 nm transistor gate length has become a reality, been synthesized by atomic layer deposition (ALD) technique of ZrO2 in single wall carbon nanotube (SWNCT) followed by pick-and-place dry transfer of MoS2 onto the SWCNT covered with ZrO27. This capacity to synthesize nanomaterials made of few atoms with pre-defined positions allow the development of devices with specific properties and application in many technological areas8, such as chemical sensors9–11, solar cells12,13 and catalysis14–22, for example. This broad range of applications of nanomaterials in industry rise due to the ability to modulate their electronic structure by tailoring the composition, size and topology, thus resulting in a non-linear size dependent behavior of their properties not observed in the bulk material Also, they can be applied as building blocks in the development of new materials23–25. Recently, functionalized nanodevices built by silicon clusters arranged in cage-like structures were proposed because of their unique opto-electronic and magnetic properties26. Silicon nanoclusters, Si2-11, were extensively studied by previous theoretical and experimental works27–37, focusing in the global minimum, electronic structure and magic numbers (n = 4, 6, 7 and 10) characterization. However, none of these works have systematically sought to characterize cage-like silicon clusters, since their experimental formation is difficult to obtain because of their sp3 hybridization38–40. Cage-like structure was experimentally observed on small fullerene cluster (C28) by uranium addition inside the

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

Page 4 of 37

4

carbon cluster41. Motivated by this work41, Jackson and Nellermoe40 studied the chemical substitution on Si20 cage structure with carbon, oxygen, germanium, tin, selenium and zirconium atoms by means of density function theory (DFT) within local density functional approximation (LDA), and found zirconium as the best candidate atom to stabilize the silicon cage. Previous theoretical and experimental works26,42–53 have also demonstrated that the transition metal atoms are encapsulated by silicon clusters, showing a relation between the transition metal valence with the silicon cluster size. Kulshreshtha et al.47 doped Si10 cluster with lithium, sodium, beryllium, magnesium, boron, aluminum and carbon using B3LYP/631G(d) methodology, showing that the beryllium atom is the only one able to form the cagelike structure. The choice of aluminum as dopant atom in the silicon clusters is motivated by at least four observations: firstly, the self-assembled aluminum-silicon nanowires networks on glass and silicon substrates54; secondly, the electrical conductivity enhancement in silicon nanowires (SiNW) as compared to pristine SiNW55. Thirdly, the formation of superlattice on the Si (111) 7x7 surface of identical-size nanoclusters (magic clusters)56 and, lastly, by its electronic structure that has one less electron than the silicon, corresponding to the cation Si. The phosphorus atom was selected based on its application on the semiconductor industry as p-dopant type, its behavior as the anion Si and electron donor in SiNW57. Previous theoretical works51,58 characterized the global minimum of the aluminum doped silicon clusters series (Si1-13Al) through density functional theory (B3LYP/6-31(d)) and molecular dynamics. Energetics properties of these molecular systems such as binding energy, dissociation energy and fragmentation energy were determined. Tam et al.59,60 studied the stability and electronic structure of both aluminum and phosphorus doped silicon

ACS Paragon Plus Environment

Page 5 of 37 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

5

clusters, SinAlm// and SinPm// with n = 1 – 10 and m = 1 – 2, by DFT methodology (B3LYP/6-31+G(d)). The stability analysis was also carried out by dissociation and clustering energies, while the electronic structure was described employing the jellium model61, in which density of states (DOS) was performed for the closed shell clusters. The present work focuses in study the chemical doping effects caused by subvalent (aluminum) and supervalent (phosphorus) atoms in the physical-chemical properties on silicon clusters Sin, n = 6, 7, 10 and 11. These systems were chosen by their remarkable high/low (Si6, Si7, Si10) or low/high (Si11) stability/reactivity, when compared to neighbor clusters (in size). In particular, Si11 was chosen as a counterexample of the most stable clusters (Si6, Si7, Si10) in the silicon series (Si2-11). Therefore, from a quantum chemical perspective, it is an interesting system to be investigated in terms of electronic structure changes induced by the doping species that could eventually lead to improved stabilization properties. To analyze the physical-chemical properties of the doped silicon clusters, simple stability and reactivity descriptors derived from wavefunction methodologies were considered, namely, the Hartree-Fock method and the multiconfigurational approach CASPT2. As a starting point to investigate the electronic properties, the Hartree-Fock orbitals diagram was employed, in special, the highest occupied molecular orbital (HOMO). The multiconfigurational CASPT2 method was applied to compute the vertical excitation energy (Te) and the electronic density variation due to the doping. The chemical stability was probed by means of our previously developed 3 function37, which encompasses both thermodynamical and kinetic aspects of the systems based on the HOMO energy, vertical excitation energy and the atomization free energy. Also, the reactivity was studied using Fukui’s functions62–66 calculated by electronic density of the neutral and ionic species.

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

Page 6 of 37

6

II. Computational Methods The initial structures can be built by using sophisticated techniques such as the particle swarm optimization method67–78 or global optimization methods, e.g. basin hopping (BH)79, simulated annealing80, minima hopping (MH)81,82, genetic algorithm (GA)82,83, funneling hopping (FH)84 and random sampling search85. Alternatively, an initial space of isomers can also be obtained by simpler techniques, as it was employed in the present work, by performing substitution23,44,47,50,86 and addition (exohedral) routes 23,42,47,52,87, as shown in Fig.1. Also, point group symmetry was considered to construct structures27–36 which were not taken into account by the substitution or addition routes when incorporating aluminum and phosphorus atoms on the pure silicon clusters37. For each starting structure, the allowed spin multiplicities were also considered. Single doped clusters exhibit an odd number of electrons, then only doublet and quartet spin multiplicity were considered. For the double doped structures, singlet and triplet states were computed. The initial structures were optimized using density functional theory within a hybrid meta-GGA approximation, M0688, combined with 6-311++G(3df,3pd) basis set89. Frequencies analysis was carried out to identify any possible saddle point and the global minimum was determined employing zeropoint energy correction on the total electronic energy. The ground state total energy was further corrected with single-point calculations by perturbation theory to the second order (MS-CASPT2)90–95, which uses as the zero order wave function the complete active self-consistent field (CASSCF)96,97 at the same atomic basis set, 6-311++G(3df,3pd)89. The CASSCF wave function was built in a state-average fashion, including the ground state and the first excited electronic state. For each doped silicon cluster, an analysis of the natural orbitals occupations of a preliminary averaged quadradic coupled cluster (AQCC) calculation was performed to select the most adequate ACS Paragon Plus Environment

Page 7 of 37 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

7

active space orbitals for the CASSCF and CASPT2 wavefunctions. Orbitals with occupations in the range of 0.1 to 1.9 were included in the CASSCF wavefunction and correlated in the perturbation step, the resulting active spaces for all systems are displayed on Table S1. In the CASPT2 calculation, a level shift (0.1 hartree) was included to remove possible intruder states in the excited states calculations98.

Figure 1. Doped silicon clusters addition and substitution scheme, showing some doping possibilities. All the non-equivalent doping positions were considered as starting geometries of a given system. Red and blue colors refers to the dopant. The chemical substitution effects on the electronic structure of the silicon clusters were evaluated through the molecular orbitals energy diagrams, vertical excitation energy and the electronic density variation under the addition (phosphorus) or removal (aluminum) of an electron. Hartree-Fock orbitals energy diagrams were selected to analyze the electronic structure of the silicon doped clusters since the HOMO energy have physical significance, as the ionization energy by Koopman's theorem. From the characterization of the first excited

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

Page 8 of 37

8

electronic states, the vertical excitation energy and electronic density difference between the excited state and ground state were obtained using the single-point CASPT2//M06/6311++G(3df,3pd) methodology. Applying the single-point methodology, the electron addition or removal effects caused by the subvalent and supervalent substitution on the electronic density of the silicon clusters was evaluated performing an electronic density difference between the doped silicon cluster and pure silicon cluster. An estimate of the chemical affinity between the silicon clusters and the doping atoms was performed through the 3 function, recently proposed in our previous work.37 The 3 function37 is a stability ranking function composed by a thermodynamical descriptor, the atomization free energy (∆G°atom, 298K), and two indirect kinetics descriptors, the HOMO energy (homo) and vertical excitation energy (Te). ε3 = |∆𝐺°𝑎𝑡𝑜𝑚, 298𝐾| ∗ | εℎ𝑜𝑚𝑜| ∗ |𝑇𝑒|

(1)

The reactivity of the doped silicon clusters was studied using the Fukui’s functions62– 66, f(r)

and f(r). It has been showed that the Fukui’s functions are useful tools to identify

nucleophilic, electrophilic and radicalar chemical reactions sites62–66, since these functions are based on the variation of the number of electrons in the molecule. Anionic and cationic structures were optimized using the same methodology applied on the neutral clusters (M06/6-311++G(3df,3pd)). From the optimized geometries, the Fukui’s functions were calculated using the electronic density and Mulliken charge difference between the neutral and ionic structures, in both cases. The electronic densities of the neutral and ionics species were calculated separately using the single-point CASPT2//M06/6-311++G(3df,3pd) methodology, as given by the equations 1 - 4 below.

ACS Paragon Plus Environment

Page 9 of 37 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

9

𝑓 + (𝑟) = 𝜌𝑆𝑖𝑛―― 1𝐴𝑙/𝑃 ― 𝜌𝑆𝑖𝑛 ― 1𝐴𝑙/𝑃

(2)

𝑓 + (𝑟) = 𝜌𝑆𝑖𝑛―― 2𝐴𝑙𝑃 ― 𝜌𝑆𝑖𝑛 ― 2𝐴𝑙𝑃

(3)

𝑓 ― (𝑟) = 𝜌𝑆𝑖𝑛 ― 1𝐴𝑙/𝑃 ― 𝜌𝑆𝑖𝑛+― 1𝐴𝑙/𝑃

(4)

𝑓 ― (𝑟) = 𝜌𝑆𝑖𝑛 ― 2𝐴𝑙𝑃 ― 𝜌𝑆𝑖𝑛+― 2𝐴𝑙𝑃

(5)

The electronic density differences were plotted using an isosurface value of 0.003. The difference between the doped and pure silicon clusters was plotted on the doped cluster geometry, and the Fukui’s functions difference was plotted on the neutral cluster geometry. Performing the electronic density calculation in separated geometries did not show any significant changes from the electronic densities calculated in the same geometry. The electronic density colors red and blue shows the electronic redistribution on the molecule. The red color refers to electron poor regions in the molecule and the blue color is related to electron rich regions when compared to the reference density, i.e., the anion compared to the neutral system for f(r) and the neutral compared to the cation system for f(r). The geometry optimization was performed using the Gaussian 09 package99 and the CASSCF and CASPT2 calculations by Molpro 2015100.

III. Results and Discussions The atomization Free Gibbs energy, atomization enthalpy, electronic excitation energy, transition moment, electronic configuration, the cartesian coordinates of the characterized isomers are displayed in the supplementary material on Fig.S1, Table S2, S3, S4, respectively. Reactive parameters such as Fukui’s functions (f f f), global hardness

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

Page 10 of 37

10

(), local softness (sk, sk, sk) (calculated by Mulliken charge) are gathered on Table S5 and the Mulliken charge are collected on Table S6.

A. Structural properties The optimized ground state geometries of Sin-a-bAlaPb characterized by M06/6311++G(3df,3pd) are shown in Fig.2, while the cartesian coordinates are presented in Table S4 in the supplementary material.

Figure 2. Ground state geometries optimized of Sin-a-bAlaPb (n= 6, 7, 10 and 11; a, b =0, 1) with 6-311++G(3df,3pd). Orange refers to the phosphorus atom and beige to the aluminum atom. The chemical substitution behavior of the aluminum and phosphorus atoms on the silicon clusters differs from the one observed in semiconductors of type p and n. In solidstate physics, the chemical doping behavior substitutional (substituting matrix atom for

ACS Paragon Plus Environment

Page 11 of 37 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

11

dopant atom) or interstitial (dopant atom added in structural sites) is predicted by HumeRothery101 rules, but these rules might not be applicable to nanometric systems. If the HumeRothery rules were valid, the aluminum doping would have a substitutional character while the phosphorus doping would present an interstitial behavior in the case of silicon crystal. However, in the nanoclusters systems, the aluminum and phosphorus doping character are correlated with the cluster size, as illustrated in Fig.2. In the case of six atoms clusters, n = 6, aluminum atom substitutes one silicon atom on Si6, exhibiting a substitutional behavior and the phosphorus chemical substitution presents a hollow addition on Si5, characterizing an interstitial doping, which, in this case, agrees with Hume-Rothery rules101. But in the Si4AlP cluster the doping behavior is inverted, the phosphorous doping presents a substitutional character and the aluminum chemical substitution has a bridge addition site. For the remaining doped clusters, n = 7, 10 and 11, the aluminum and phosphorus chemical doping behavior are the same for clusters with the same size. For n = 7, the chemical substitution presents a substitutional character, while on the n = 10 and 11 possess a hollow addition nature. Based on the Mulliken charge and the electronic density distribution on the doped silicon clusters, Table S6 and Fig.4, respectively, the chemical substitution sites on the pure silicon clusters can be associated to the charge transferred between the dopants and the silicon cluster. The aluminum sites tend to be the ones that provide less electrons to the silicon cluster, while the phosphorus sites are the ones that acquire electrons from the silicon atoms. The charge transferred is also related to the dopant coordination number and the cluster size. The coordination number (number of neighboring atoms) can be also associated to the dopant electronegativity. In the subvalent chemical substitution, the aluminum coordination number

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

Page 12 of 37

12

tend to be as minimal as possible, since presents lower values of electronegativity. The opposite occurs in the case of the supervalent chemical substitution, where the phosphorus coordination number tend to be as maximal as possible. It is found that there is an enhancement of the charge transfer effect as the cluster size increases, (probably) due to the large pool of valence electrons provided by the neighboring silicon atoms. Therefore, as bigger is the cluster, more electrons are available and consequently, more negative charge is transferred to the dopant. However, there are some exceptions to the cluster size effect on the charge transfer, such as Si10Al and Si10P. On the Si10Al cluster, the aluminum atom coordinates only three atoms (as in Si9Al) rather than four, as seen in Si5Al and Si6Al. In this way, it acquires more negative charge from the silicon cluster. On the other hand, the phosphorus atom on the Si10P has its lower coordination number and also the lowest negative charge. In contrast to the donor character of the phosphorus atom seen in crystalline silicon or silicon nanoparticles, the phosphorus atoms present an acceptor behavior when doping small silicon clusters such as those studied in this work. Therefore, there is a balance between the cluster size and the doping atom coordination number on the charge transferred to the dopant. B. Electronic structure In general, molecular systems near and in the equilibrium region can be adequately described by a monoconfigurational wavefunction. In this sense, the Hartree-Fock method can be considered as a good starting-point to qualitatively describe the electronic groundstate, providing molecular orbital energies with physical interpretation.

ACS Paragon Plus Environment

Page 13 of 37 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

13

Figure 3. Orbital energy diagram obtained with the RHF//M06/6-311++G(3df,3pd) method for the studied clusters. The red and blue colors refer to the contribution of the aluminum and phosphorous atomic orbitals to each molecular orbital. The physical interpretation of the Hartree-Fock orbital is focused in the frontier orbitals, HOMO and LUMO. As it is well known, the HOMO energy can represent the ionization potential (Koopman’s theorem) being associated to the strength of the weakest electron bonded in the molecular system and indicate, in some extent, the stability of the cluster. Although more inaccurately due to the lack of electron correlation, the LUMO energy can still be correlated to the electron affinity, providing information about the electronacceptor character of the system. Based on this interpretation, the chemical substitution (Al,

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

Page 14 of 37

14

P, or AlP) on all studied silicon clusters (Si6, Si7 and Si10), except for Si11, result in less bonded doped clusters, as suggested by the increase of the HOMO energy of these systems, shown in Fig.3. As expected, the silicon clusters with high stability (Si6, Si7 and Si10) do not accept substitutions, even the isoelectronic substitution (a subvalent and a supervalent substitution). On the other hand, the highly unstable Si11, benefits from any substitution, but mostly from the subvalent and supervalent chemical doping. In addition to providing information on the physical-chemical properties, ionization potential and electron-acceptor behavior, the Hartree-Fock orbitals diagram indicates how occurs the chemical bond among the atoms. Variations in the molecular orbital energy of the pure clusters indicates how the atomic orbitals of the dopants interacts with the silicon clusters molecular orbitals. Comparing the subvalent and supervalent substitution changes on the pure silicon cluster molecular orbitals (Fig.3), the aluminum atomic orbitals participate more of the frontier orbitals formation than the phosphorus orbitals, causing lower variations on the HOMO energy than the phosphorus substitution, suggesting the possibility of the Sin1Al

clusters be more bonded than the Sin-1P. However, the subvalent substitution also lowers

the energy of molecular orbitals that are mostly formed by silicon atoms, indicating an increment on the cohesion energy, and suggests that the phosphorus orbitals contribute more on the internal molecular orbitals.

ACS Paragon Plus Environment

Page 15 of 37 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

15

Figure 4. Electronic density difference between the doped and pure silicon clusters obtained through CASPT2//M06/6-311++G(3df,3pd) with isosurface value of 0.003 𝑒.𝑏𝑜ℎ𝑟 ―3. Yellow and white circles refers to aluminum and phosphorus atoms, respectively. In order to qualitatively evaluate the charge transfer, discussed in the geometry section, the electronic density between the doped and pure silicon clusters was employed to analyze the addition (phosphorus) or removal (aluminum) of one electron in the electronic density, displayed in Fig.4. In the smallest aluminum doped clusters, n = 6 and 7, the dopant atom is under poor regions of electronic density (red color), thus suggesting a charge transfer from the aluminum to the neighboring silicon atoms. However, this charge transfer direction is inverted on the larger clusters (n = 10 and 11), where the aluminum atom is under the rich region of electronic density (blue color), acting as an electron acceptor center on the silicon cluster. On the other hand, the supervalent substitution redistributes a significant part of the electronic density to the dopant surroundings, and therefore, the phosphorus atom behaves as

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

Page 16 of 37

16

electron acceptor site, due its higher electronegativity. When both doping atoms are present on the doped cluster, the electrons are basically transferred from the aluminum atom to the phosphorus atom direction, especially on the Si4AlP and Si5AlP. These results are consistent with the analysis of atomic charge distribution based on Mulliken population. (Table S6 in the supplementary material).

Figure 5. Electronic difference between the ground state and the first excited state of the doped silicon clusters obtained using CASPT2//M06/6-311++G(3df,3pd) with isosurface value of 0.003 𝑒.𝑏𝑜ℎ𝑟 ―3. A quantitative analysis of the electronic structure (Tables S2 and S3) was carried out by vertical excitations on the electronic ground state through the single-point CASPT2 methodology. For the phosphorus doped nanoclusters, Sin-1P and Sin-2AlP, the first excited state corresponds to one electron promotion characterized as the HOMO–LUMO transition, observed previously for the Si2-11 series27–37, except for Si5P cluster, in which the transition

ACS Paragon Plus Environment

Page 17 of 37 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

17

corresponds to one electron excitation from the HOMO-1 to singly occupied molecular orbital (SOMO), as seen for the aluminum doped clusters (Sin-1Al). This HOMO-1–SOMO transition has lower excitation energy than the HOMO–LUMO transition. Concerning to the electronic ground state, the wave function of almost all doped (Al or P and AlP) structures preservers the monoconfigurational character observed on the pure silicon clusters37, except for Si5P and Si9Al clusters. As one goes to the excited states, there is a significant wave function multiconfigurational character for the Si5P, Si5AlP, Si9Al, Si8AlP, Si10Al and Si9AlP clusters, as seen in Table S3. To evaluate the consistency of the Hartree-Fock method to describe the electronic structure, we can compare the electron displacement on the cluster provided by the electronic density difference between the ground state and excited state calculated by our single-point CASPT2 methodology, with the atomic orbitals contribution on the molecular orbital construction obtained by the RHF methodology. The electronic transitions occurs among the HOMO-1, HOMO (isoelectronic) or SOMO (subvalent and supervalent) and LUMO, these molecular orbitals are considered in the analysis of Fig.3. Taken together, the results of Fig.3 and Fig.5, one can see that the atoms largely contributing for the electronic density variations upon excitation are the same which have strong participation in the construction of the Hartree-Fock molecular orbitals. The Hartree-Fock electronic structure differs from the CASPT2 for the Si10Al and Si10P clusters. Therefore, the Hartree-Fock wave function is in qualitatively good agreement with the CASPT2 methodology.

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

Page 18 of 37

18

C. Stability and reactivity The stability and reactivity of a molecular system can be used to estimate which chemical substitution and silicon cluster combination produce favorable doped clusters. Variations in the molecular system stability can be used as a parameter to evaluate the chemical affinity of a given molecular system with respect to the introduction of a dopant atom. Which is investigated here by means of the 3 function (Fig.6). The ability of a molecule to react is related to the electronic density of the frontier orbitals, since these regions represent the reactive regions to an electrophilic and nucleophilic reaction. To determine these reactive sites, Fukui’s functions (f(r) and f(r))62–66 (Fig.7) use the variation of the spatial electronic distribution (density) in respect to the number of electrons, as given in equations 1 - 4.

Figure 6. Stability ranking function (3), applied on doped silicon clusters. Color black, red, blue and green refers to n =6, 7, 10 and 11, respectively.

ACS Paragon Plus Environment

Page 19 of 37 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

19

As discussed in detail in our previous study37, small values of 3 indicates low stability, while large values represents high stability. In this way, in Fig.6, the 3 function37 demonstrates that the Si11 has a higher acceptance degree to the chemical substitutions (Al, P or AlP) than the Si6, Si7 and Si10 clusters, since the doped silicon clusters (Si10Al, Si10P and Si9AlP) possess higher values of 3 than the Si11, indicating a higher stability. The Si6, Si7 and Si10 clusters have lower acceptance degree because they are more stable than their respective doped structure. These results could be expected, based on the fact that, the Si11 is one of the less stable on the Si2-11 series and the Si6, Si7 and Si10 are the most stable ones37. In general, the double substitution (AlP) produces more stable clusters than the monodoping (Al or P), once they combine both aluminum and phosphorus stability properties; the increase on the strength of the weakest chemical bond and the increase on the cohesion energy, respectively.

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

Page 20 of 37

20

Figure 7. Fukui's functions calculated using electronic density difference between the neutral and ionics species obtained using CASPT2//M06/6-311++G(3df,3pd) with isosurface value of 0.003 𝑒.𝑏𝑜ℎ𝑟 ―3. Yellow and white circles refers to aluminum and phosphorus atoms, respectively. From the Fukui's functions, Fig.7, in general, one can infer that the chemical substitution tends to increase the Si6, Si7 and Si10 reactivity due to the significant increment of the electronic density volume. However, no significant variation on the electronic density volume were observed for the Si6Al and Si9P (electrophilic reactions, f(r)), and Si4AlP, Si9Al and Si9P (nucleophilic reactions, f(r)). The chemical substitutions effects on the Si11 reactivity, suggest to produce high reactive species (Si10Al, Si10P), since their electronic density volume are equal or higher than the Si11. This extremely reactive character occurs by the presence of the electronic density covering the whole doped structure, except for Si10P,

ACS Paragon Plus Environment

Page 21 of 37 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

21

whose reactivity drastically decreases, indicating a kinetic stability for electrophilic reactions. Combining the stability and reactivity analysis, one can conclude that the clusters with higher stability and lower reactivity are prompter to be formed. Magic numbers (n= 6, 7 and 10) still are the stable ones among the studied clusters, considering clusters with same size. And, the Si11 possess, exclusively, a chemical substitution acceptor character, if stability analysis is considered individually. However, if the reactivity is also considered, the results indicate that the synthesis of the Si10P can be achieved. The others studied doped silicon clusters are more unstable and reactive than their pure silicon counterparts.

IV. Conclusions The doping effects in the physical-chemistry properties caused by subvalent and supervalent atoms, aluminum and phosphorus, on the Si6, Si7, Si10 and Si11 silicon clusters were studied using molecular quantum chemistry methods. The applicability of the HumeRothery rules, were tested in nanometric systems, comparing the aluminum and phosphorus substitution character observed on the doped clusters with the Hume-Rothery predictions. It was observed that the chemical doping character changes according to the cluster size, i.e., a size-dependent property, differing from the Hume-Rothery predictions. The most probable substitution sites were determined combining the Mulliken charge and the electronic density difference between the doped and pure silicon cluster. We noticed that aluminum and phosphorus preferable sites are related to the negative charge transferred between the dopant and the silicon cluster, where the aluminum position are the ones that transfer less electrons

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

Page 22 of 37

22

and the phosphorus sites the ones that acquires electrons. Using the energy diagram of the Hartree-Fock orbitals, we could, in a qualitative approach, evaluate changes on the weakest bond and electron acceptor behavior of the systems, showing that the chemical substitution, except for Si11, produced less bonded clusters. With respect to the wave function character, the chemical substitutions maintain the ground state monoconfigurational behavior observed on the silicon clusters. Regarding the vertical excitation energies, the subvalent substitution lowered the excitation energy by changing the molecular orbitals involved on the electronic transition. Combining the stability and reactivity analysis, the most favorable chemical substitution on the studied silicon clusters can be pointed out. Considering an electrophilic environment, the Si10P is the most probable cluster to be found, while for nucleophilic reactions, none of the doped clusters resulted in higher stability. Supporting Information The atomization Free Gibbs energy, atomization enthalpy, electronic excitation energy, transition moment, electronic configuration, the cartesian coordinates of the characterized isomers are displayed in the supplementary material on Fig.S1, Table S2, S3, S4, respectively. Reactive parameters such as Fukui’s functions (f f f), global hardness (), local softness (sk, sk, sk) (calculated by Mulliken charge) are gathered on Table S5 and the Mulliken charge are collected on Table S6.

Acknowledgements This work has been supported by Brazilian agencies Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) under grants 2017/07707-3 and 2017/01359-3, Conselho

ACS Paragon Plus Environment

Page 23 of 37 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

23

Nacional de Desenvolvimento Científico e Tecnológico (CNPq) under grants 307052/20168,

404337/2016-3,

309051/2016-9

and

406107/2016-5

and

Coordenação

de

Aperfeiçoamento de Pessoal de Nível Superior (CAPES) under Project CAPES/ITA Project No. 88882.161993/2017-01.

References (1)

Waldrop, M. M. The Chips Are down for Moore’s Law. Nature 2016, 530, 144–147.

(2)

Mack, C. A. Fifty years of Moore's law. IEEE Transactions on Semiconductor Manufacturing 2011, 24, 202207.

(3)

Schaller, R. R. Moore’s Law: Past, Present and Future. IEEE Spectrum 1997, 34, 52– 59.

(4)

Keyes, R. . The Impact of Moore ’ S Law. IEEE Solid-State Circuits Society. Newsletter 2006, 11, 25–27.

(5)

Bohr, M. 14 nm Transistor Explained—Following the Path of Moore’s Law https://www.intel.com.br/content/www/br/pt/silicon-innovations/standards-14nmexplained-video.html (accessed Feb 7, 2017).

(6)

Anthony, S. IBM unveils world’s first 5nm chip https://arstechnica.com/gadgets/2017/06/ibm-5nm-chip/ (accessed Apr 17, 2018).

(7)

Desai S.B; Madhvapathy, S. R.; Sachid, A. B.; Llinas, J. P.; Wang, Q.; Ahn, G. H.; Pitner, G.; Kim, M. J.; Bokor, J.; Hu, C.; et al. MoS2 Transistors with 1-Nanometer Gate Lengths. Science 2016, 354, 99–102.

(8)

Fahlman, B. D. Materials Chemistry; Springer, 2007.

(9)

Hussain, A. M. P.; Sarangi, S. N.; Kesarwani, J. A.; Sahu, S. N. Au-Nanocluster Emission Based Glucose Sensing. Biosens. Bioelectron. 2011, 29, 60–65.

(10)

Zhang, L.; Wang, E. Metal Nanoclusters: New Fluorescent Probes for Sensors and Bioimaging. Nano Today. 2014, 9, 132–157.

(11)

Franke, M. E.; Koplin, T. J.; Simon, U. Metal and Metal Oxide Nanoparticles in Chemiresistors: Does the Nanoscale Matter? Small 2006, 2, 36–50.

(12)

Chen, Y. S.; Choi, H.; Kamat, P. V. Metal-Cluster-Sensitized Solar Cells. A New Class of Thiolated Gold Sensitizers Delivering Efficiency Greater than 2%. J. Am. Chem. Soc. 2013, 135, 8822–8825.

(13)

Choi, H.; Chen, Y. S.; Stamplecoskie, K. G.; Kamat, P. V. Boosting the Photovoltage of Dye-Sensitized Solar Cells with Thiolated Gold Nanoclusters. J. Phys. Chem. Lett. 2015, 6, 217–223.

(14)

Ferreira, H. S.; Rangel, M. do C. Nanotecnologia: Aspectos Gerais e Potencial de

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

Page 24 of 37

24

Aplicação em Catálise. Quim. Nova 2009, 32, 1860–1870. (15)

Li, G.; Abroshan, H.; Liu, C.; Zhuo, S.; Li, Z.; Xie, Y.; Kim, H. J.; Rosi, N. L.; Jin, R. Tailoring the Electronic and Catalytic Properties of Au25 Nanoclusters via Ligand Engineering. ACS Nano 2016, 10, 7998–8005.

(16)

Nilsson, A.; Pettersson, L. G. M.; Hammer, B.; Bligaard, T.; Christensen, C. H.; Nørskov, J. K. The Electronic Structure Effect in Heterogeneous Catalysis. Catal. Letters 2005, 100, 111–114.

(17)

Vajda, S.; White, M. G. Catalysis Applications of Size-Selected Cluster Deposition. ACS Catal. 2015, 5, 7152–7176.

(18)

Zhang, J.; Wang, H.; Dalai, A. K. Effects of Metal Content on Activity and Stability of Ni-Co Bimetallic Catalysts for CO2 reforming of CH4. Appl. Catal. A 2008, 339 , 121–129.

(19)

Ding, F. The Importance of Strong-Carbon-Metal Adhesion for Catalytic Nucleation of Single-Walled Carbon Nanotubes. Nano Lett. 2008, 8, 463–468.

(20)

Li, G.; Jin, R. Atomically Precise Gold Nanoclusters as New Model Catalysts. Acc. Chem. Res. 2013, 46, 1749–1758.

(21)

Guerreiro, E. D.; Gorriz, O. F.; Larsen, G.; Arrúa, L. A. Cu/SiO2 Catalysts for Methanol to Methyl Formate Dehydrogenation A Comparative Study Using Different Preparation Techniques. Appl. Catal. A 2000, 204, 33–48.

(22)

Ouyang, R.; Jiang, D. Understanding Selective Hydrogenation of Α,β-Unsaturated Ketones to Unsaturated Alcohols on the Au25(SR)18 Cluster. ACS Catal. 2015, 5, 6624–6629.

(23)

Khanna, S. N.; Jena, P. Atomic Clusters: Building Blocks for a Class of Solids. Phys. Rev. B 1995, 51, 13705–13716.

(24)

Castleman, A. W.; Khanna, S. N. Clusters, Superatoms, and Building Blocks of New Materials. J. Phys. Chem. C 2009, 113, 2664–2675.

(25)

Castleman Jr., A. W.; Khanna, S. N. Superatoms: Building Blocks of New Materials. J. Phys. Chem. C 2009, 113, 2664–2675.

(26)

Ohara, M.; Koyasu, K.; Nakajima, A.; Kaya, K. Geometric and Electronic Structures of Metal (M)-Doped Silicon Clusters (M = Ti, Hf, Mo and W). Chem. Phys. Lett. 2003, 371, 490–497.

(27)

Bloomfield, L. A.; Freeman, R. R.; Brown, W. L. Photofragmentation of MassResolved Si2-12+ Clusters. Phys. Rev. Lett. 1985, 54, 2246–2249.

(28)

Heath, J. R.; Liu, Y.; O’Brien, S. C.; Zhang, Q.-L.; Curl, R. F.; Tittel, F. K.; Smalley, R. E. Semiconductor Cluster Beams: One and Two Color Ionization Studies of SiX and GeX. J. Chem. Phys. 1985, 83, 5520–5526.

(29)

Cheshnovsky, O.; Yang, S. H.; Pettiette, C. L.; Craycraft, M. J.; Liu, Y.; Smalley, R. E. Ultraviolet Photoelectron Spectroscopy of Semiconductor Clusters: Silicon and Germanium. Chem. Phys. Lett. 1987, 138, 119–124.

(30)

Raghavachari, K.; Rohlfing, C. M. Bonding and Stabilities of Small Silicon Clusters: A Theoretical Study of Si7Si10. J. Chem. Phys. 1988, 89, 2219–2234.

ACS Paragon Plus Environment

Page 25 of 37 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

25

(31)

Raghavachari, K.; Rohlfing, C. M. Electronic Structures of the Negative Ions Si2Si10: Electron Affinities of Small Silicon Clusters. J. Chem. Phys. 1991, 94, 3670–3678.

(32)

Curtiss, L. A; Deutsch, P. W.; Raghavachari, K. Binding Energies and Electron Affinities of Small Silicon Clusters (n=2-5). J. Chem. Phys. 1992, 96, 6868–6872.

(33)

Tománek, D.; Schluter, M. A. Calculation of Magic Numbers and the Stability of Small Si Clusters. Phys. Rev. Lett. 1986, 56, 10551058.

(34)

Zhu, X.; Zeng, X. C. Structures and Stabilities of Small Silicon Clusters: Ab Initio Molecular-Orbital Calculations of Si7Si11. J. Chem. Phys. 2003, 118, 3558–3570.

(35)

Pouchan, C.; Bégué, D.; Zhang, D. Y. Between Geometry, Stability, and Polarizability: Density Functional Theory Studies of Silicon Clusters Sin (n=3-10). J. Chem. Phys. 2004, 121, 4628–4634.

(36)

Yang, J.; Xu, W.; Xiao, W. The Small Silicon Clusters Sin (n=2-10) and Their Anions: Structures, Themochemistry, and Electron Affinities. J. Mol. Struct. THEOCHEM. 2005, 719, 89–102.

(37)

Fernandes, G. F. S.; Machado, F. B. C.; Ferrão, L. F. A. A Quantitative Tool to Establish Magic Number Clusters, ε3, Applied in Small Silicon Clusters, Si2-11. J. Mol. Model. 2018, 24, 203212.

(38)

Huda, M. N.; Ray, A. K. Carbon Dimer in Silicon Cage: A Class of Highly Stable Silicon Carbide Clusters. Phys. Rev. A 2004, 69, 1–4.

(39)

Li, B. X.; Cao, P. L.; Que, D. L. Distorted Icosahedral Cage Structure of Si60 Clusters. Phys Rev B 2000, 61, 1685–1687.

(40)

Jackson, K.; Nellermoe, B. Zr@Si20: A Strongly Bound Si Endohedral System. Chem. Phys. Lett. 1996, 254, 249–256.

(41)

Guo, T.; Diener, M. D.; Chai, Y.; Alford, M. J.; Haufler, R. E.; Mcclure, S. M.; Weaver, J. H.; Scuseria, G. E.; Smalley, R. E. Uranium Stabilization of C28 : A Tetravalent Fullerene. Science 1992, 257, 1661–1664.

(42)

Chuang, F. C.; Hsieh, Y. Y.; Hsu, C. C.; Albao, M. A. Geometries and Stabilities of Ag-Doped Sin (n=1-13) Clusters: A First-Principles Study. J. Chem. Phys. 2007, 127, 1–6.

(43)

Kawamura, H.; Kumar, V.; Kawazoe, Y. Growth Behavior of Metal-Doped Silicon Clusters SinM(M=Ti,Zr,Hf;n=8–16). Phys. Rev. B 2005, 71, 1–12.

(44)

Jaeger, J. B.; Jaeger, T. D.; Duncan, M. A. Photodissociation of Metal-Silicon Clusters: Encapsulated versus Surface-Bound Metal. J. Phys. Chem. A 2006, 110, 9310–9314.

(45)

Xu, H. G.; Zhang, Z. G.; Feng, Y.; Yuan, J.; Zhao, Y.; Zheng, W. Vanadium-Doped Small Silicon Clusters: Photoelectron Spectroscopy and Density-Functional Calculations. Chem. Phys. Lett. 2010, 487, 204–208.

(46)

Hagelberg, F.; Yanov, I.; Leszczynski, J. Theoretical Investigations on Closed-Shell Silicon Clusters Doped with Cu Atoms. J. Mol. Struct. THEOCHEM. 1999, 487, 183– 192.

(47)

Majumder, C.; Kulshreshtha, S. K. Impurity-Doped Si10 Cluster: Understanding the

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

Page 26 of 37

26

Structural and Electronic Properties from First-Principles Calculations. Phys. Rev. B 2004, 70, 245426-1–245426-7. (48)

Karamanis, P.; Marchal, R.; Carbonniére, P.; Pouchan, C. Doping-Enhanced Hyperpolarizabilities of Silicon Clusters: A Global Ab Initio and Density Functional Theory Study of Si10(Li, Na, K)n (n = 1, 2) Clusters. J. Chem. Phys. 2011, 135, 0445111–044511-15.

(49)

Neukermans, S.; Wang, X.; Veldeman, N.; Janssens, E.; Silverans, R. E.; Lievens, P. Mass Spectrometric Stability Study of Binary MSn Clusters (S=Si, Ge, Sn, Pb, and M=Cr, Mn, Cu, Zn). Int. J. Mass Spectrom. 2006, 252, 145–150.

(50)

Guo, L.; Zhao, G.; Gu, Y.; Liu, X.; Zeng, Z. Density-Functional Investigation of Metal-Silicon Cage Clusters MSin (M= Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn; N= 816). Phys. Rev. B 2008, 77, 195417-1195417-1.

(51)

Majumder, C.; Kulshreshtha, S. K. Influence of Al Substitution on the Atomic and Electronic Structure of Si Clusters by Density Functional Theory and Molecular Dynamics Simulations. Phys. Rev. B 2004, 69, 115432-1115432-8.

(52)

Li, Y.; Tam, N. M.; Claes, P.; Woodham, A. P.; Lyon, J. T.; Ngan, V. T.; Nguyen, M. T.; Lievens, P.; Fielicke, A.; Janssens, E. Structure Assignment, Electronic Properties, and Magnetism Quenching of Endohedrally Doped Neutral Silicon Clusters, SinCo (N = 10–12). J. Phys. Chem. A 2014, 118, 8198–8203.

(53)

Alonso, J. A. Electronic and Atomic Structure, and Magnetism of Transition-Metal Clusters. Chem. Rev. 2000, 100, 637–678.

(54)

Paulose, M.; Grimes, C. A.; Varghese, O. K.; Dickey, E. C. Self-Assembled Fabrication of Aluminum-Silicon Nanowire Networks. Appl. Phys. Lett. 2002, 81, 153–155.

(55)

Landman, U.; Barnett, R. N.; Scherbakov, A. G.; Avouris, P. Metal-Semiconductor Nanocontacts: Silicon Nanowires. Phys. Rev. Lett. 2000, 85, 1958–1961.

(56)

Kotlyar, V. G.; Zotov, A. V.; Saranin, A. A.; Kasyanova, T. V.; Cherevik, M. A.; Pisarenko, I. V.; Lifshits, V. G. Formation of the Ordered Array of Al Magic Clusters on Si(111)7×7. Phys. Rev. B. 2002, 66, 165401-1–165401-4.

(57)

Yan, B.; Rurali, R.; Gali, Á. Ab Initio Study of Phosphorus Donors Acting as Quantum Bits in Silicon Nanowires. Nano Lett. 2012, 12, 3460–3465.

(58)

Nigam, S.; Majumder, C.; Kulshreshtha, S. K. Structural and Electronic Properties of Sin, Sin+, and AlSin-1 (n=2–13) Clusters: Theoretical Investigation Based on Ab Initio Molecular Orbital Theory. J. Chem. Phys. 2004, 121, 77567763.

(59)

Tam, N. M.; Tai, T. B.; Ngan, V. T.; Nguyen, M. T. Structure, Thermochemical Properties, and Growth Sequence of Aluminum-Doped Silicon Clusters SinAlm (n = 1-11, m = 1-2) and Their Anions. J. Phys. Chem. A 2013, 117, 6867–6882.

(60)

Pham, H. T.; Tam, N. M.; Jeilani, Y. A.; Nguyen, M. T. Structural Evolution and Bonding of Phosphorus-Doped Silicon Clusters SinPm−/0/+with n=1–10, m=1, 2. Comput. Theor. Chem. 2017, 1107, 115–126.

(61)

Brack, M. The Physics of Simple Metal Clusters: Self-Consistent Jellium Model and

ACS Paragon Plus Environment

Page 27 of 37 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

27

Semiclassical Approaches. Rev. Mod. Phys. 1993, 65, 677–732. (62)

Balawender, R.; Komorowski, L. Atomic Fukui Function Indices and Local Softness Ab Initio. J. Chem. Phys. 1998, 109, 5203–5211.

(63)

Berkowitz, M.; Parr, R. G. Molecular Hardness and Softness, Local Hardness and Softness, Hardness and Softness Kernels, and Relations among These Quantities. J. Chem. Phys. 1988, 88, 2554–2557.

(64)

Roy, R. K.; Krishnamurti, S.; Geerlings, P.; Pal, S. Local Softness and Hardness Based Reactivity Descriptors for Predicting Intra- and Intermolecular Reactivity Sequences: Carbonyl Compounds. J. Phys. Chem. A 1998, 102, 3746–3755.

(65)

Yang, W.; Parr, R. G. Hardness, Softness, and the Fukui Function in the Electronic Theory of Metals and Catalysis. Proc. Natl. Acad. Sci. 1985, 82, 6723–6726.

(66)

Fukui, K. The Role of Frontier Orbitals in Chemical Reactions (Nobel Lecture). Angew. Chemie Int. Ed. Engl. 1982, 21, 801–809.

(67)

Wang, Y.; Lv, J.; Zhu, L.; Ma, Y. CALYPSO: A Method for Crystal Structure Prediction. Comput. Phys. Commun. 2012, 183, 2063–2070.

(68)

Lv, J.; Wang, Y.; Zhu, L.; Ma, Y. Particle-Swarm Structure Prediction on Clusters. J. Chem. Phys. 2012, 137, 084104-1084104-8 .

(69)

Wang, Y.; Miao, M.; Lv, J.; Zhu, L.; Yin, K.; Liu, H.; Ma, Y. An Effective Structure Prediction Method for Layered Materials Based on 2D Particle Swarm Optimization Algorithm. J. Chem. Phys. 2012, 137, 224108-1224108-6.

(70)

Wang, Y.; Lv, J.; Zhu, L.; Ma, Y. Crystal Structure Prediction via Particle-Swarm Optimization. Phys. Rev. B 2010, 82, 1–8.

(71)

Lu, S.; Wang, Y.; Liu, H.; Miao, M. S.; Ma, Y. Self-Assembled Ultrathin Nanotubes on Diamond (100) Surface. Nat. Commun. 2014, 5, 1–6.

(72)

Le Chen, B.; Sun, W. G.; Kuang, X. Y.; Lu, C.; Xia, X. X.; Shi, H. X.; Maroulis, G. Structural Stability and Evolution of Medium-Sized Tantalum-Doped Boron Clusters: A Half-Sandwich-Structured TaB12-Cluster. Inorg. Chem. 2018, 57, 343–350.

(73)

Jin, Y.; Tian, Y.; Kuang, X.; Zhang, C.; Lu, C.; Wang, J.; Lv, J.; Ding, L.; Ju, M. Ab Initio Search for Global Minimum Structures of Pure and Boron Doped Silver Clusters. J. Phys. Chem. A 2015, 119, 6738–6745.

(74)

Xia, X. X.; Hermann, A.; Kuang, X. Y.; Jin, Y. Y.; Lu, C.; Xing, X. D. Study of the Structural and Electronic Properties of Neutral and Charged Niobium-Doped Silicon Clusters: Niobium Encapsulated in Silicon Cages. J. Phys. Chem. C 2016, 120, 677– 684.

(75)

Ju, M.; Lv, J.; Kuang, X. Y.; Ding, L. P.; Lu, C.; Wang, J. J.; Jin, Y. Y.; Maroulis, G. Systematic Theoretical Investigation of Geometries, Stabilities and Magnetic Properties of Iron Oxide Clusters (FeO)nμ (n = 1-8, μ = 0, ±1): Insights and Perspectives. RSC Adv. 2015, 5, 6560–6570.

(76)

Jin, Y.; Maroulis, G.; Kuang, X.; Ding, L.; Lu, C.; Wang, J.; Lv, J.; Zhang, C.; Ju, M. Geometries, Stabilities and Fragmental Channels of Neutral and Charged Sulfur Clusters: SnQ (n = 3-20, Q = 0, ±1). Phys. Chem. Chem. Phys. 2015, 17, 13590–13597.

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

Page 28 of 37

28

(77)

Xing, X.; Hermann, A.; Kuang, X.; Ju, M.; Lu, C.; Jin, Y.; Xia, X.; Maroulis, G. Insights into the Geometries, Electronic and Magnetic Properties of Neutral and Charged Palladium Clusters. Sci. Rep. 2016, 6, 1–11.

(78)

Sun, W. G.; Wang, J. J.; Lu, C.; Xia, X. X.; Kuang, X. Y.; Hermann, A. Evolution of the Structural and Electronic Properties of Medium-Sized Sodium Clusters: A Honeycomb-Like Na20 Cluster. Inorg. Chem. 2017, 56, 1241–1248.

(79)

Wales, D. J.; Doye, J. P. K. Global Optimization by Basin-Hopping and the Lowest Energy Structures of Lennard-Jones Clusters Containing up to 110 Atoms. J. Phys. Chem. A 1997, 101, 5111–5116.

(80)

Kirkpatrick, S.; Gelatt, C. D.; Vecchi, M. P. Optimization by Simulated Annealing. Science 1983, 220, 671–680.

(81)

Goedecker, S. Minima Hopping: An Efficient Search Method for the Global Minimum of the Potential Energy Surface of Complex Molecular Systems. J. Chem. Phys. 2004, 120, 9911–9917.

(82)

Schönborn, S. E.; Goedecker, S.; Roy, S.; Oganov, A. R. The Performance of Minima Hopping and Evolutionary Algorithms for Cluster Structure Prediction. J. Chem. Phys. 2009, 130, 144108-1144108-9.

(83)

Hartke, B. Global Geometry Optimization of Clusters Using a Growth Strategy Optimized by a Genetic Algorithm. Chem. Phys. Lett. 1995, 240, 560–565.

(84)

Cheng, L.; Feng, Y.; Yang, J.; Yang, J. Funnel Hopping: Searching the Cluster Potential Energy Surface over the Funnels. J. Chem. Phys. 2009, 130, 14112-1 214112-7.

(85)

Needs, R. J. Ab Initio Random Structure Searching. Journal of Physics: Condensed Matter. 2011, 23, 05320105322045.

(86)

Kumar, V.; Bhattacharjee, S.; Kawazoe, Y. Silicon-Doped Icosahedral, Cuboctahedral, and Decahedral Clusters of Aluminum. Phys. Rev. B 2000, 61, 8541– 8547.

(87)

Lu, J.; Nagase, S. Metal-Doped Germanium Clusters MGens at the Sizes of n = 12 and 10: Divergence of Growth Patterns from the MSin Clusters. Chem. Phys. Lett. 2003, 372, 394–398.

(88)

Zhao, Y.; Truhlar, D. G. The M06 Suite of Density Functionals for Main Group Thermochemistry, Thermochemical Kinetics, Noncovalent Interactions, Excited States, and Transition Elements: Two New Functionals and Systematic Testing of Four M06-Class Functionals and 12 Other Function. Theor. Chem. Acc. 2008, 120, 215–241.

(89)

Dunning T.H., J.; Peterson, K. A.; Wilson, A. K. Gaussian Basis Sets for Use in Correlated Molecular Calculations. X. The Atoms Aluminum through Argon Revisited. J. Chem. Phys. 2001, 11, 9244–9253.

(90)

Finley, J.; Malmqvist, P.-Å.; Roos, B. O.; Serrano-Andrés, L. The Multi-State CASPT2 Method. Chem. Phys. Lett. 1998, 288, 299–306.

(91)

Celani, P.; Werner, H.-J. Multireference Perturbation Theory for Large Restricted and

ACS Paragon Plus Environment

Page 29 of 37 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

29

Selected Active Space Reference Wave Functions. J. Chem. Phys. 2000, 112, 5546– 5557. (92)

Shiozaki, T.; Gyroffy, W.; Celani, P.; Werner, H. J. Communication: Extended MultiState Complete Active Space Second-Order Perturbation Theory: Energy and Nuclear Gradients. J. Chem. Phys. 2011, 135, 081106-1081106-4.

(93)

Shiozaki, T.; Werner, H. J. Communication: Second-Order Multireference Perturbation Theory with Explicit Correlation: CASPT2-F12. J. Chem. Phys. 2010, 133, 41103-141103-5.

(94)

Andersson, K.; Malmqvist, P. A.; Roos, B. O.; Sadlej, A. J.; Wolinski, K. SecondOrder Perturbation Theory with a CASSCF Reference Function. J. Phys. Chem. 1990, 94, 5483–5488.

(95)

Roos, B. O.; Linse, P.; Siegbahn, P. E. M.; Blomberg, M. R. A. A Simple Method for the Evaluation of the Second-Order-Perturbation Energy from External DoubleExcitations with a CASSCF Reference Wavefunction. Chem. Phys. 1982, 66, 197– 207.

(96)

Werner, H.-J.; Knowles, P. J. A Second Order Multiconfiguration SCF Procedure with Optimum Convergence. J. Chem. Phys. 1985, 82, 5053.

(97)

Knowles, P. J.; Werner, H.-J. An efficient second-order MC SCF method for long configuration expansions. Chem. Phys. Lett. 1985, 115, 259–267.

(98)

Roos, B. O.; Andersson, K. Multiconfigurational Perturbation Theory with Level Shift — the Cr2 Potential Revisited. Chem. Phys. Lett. 1995, 245, 215–223.

(99)

Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. A.; et al. Gaussian 09 Revision D.01.

(100) Werner, H.-J.; Knowles, P. J.; Knizia, G.; Manby, F. R.; Schütz, M.; Celani, P.; Györffy, W.; Kats, D.; Korona, T.; Lindh, R.; et al. MOLPRO, Version 2015.1, a Package of Ab Initio Programs. Cardiff, UK 2015. (101) Hume-Rothery, W.; Smallman, R. E.; Haworth, C. W. The Structure of Metals and Alloys. Inst. Met. 1 Carlt. House Terrace, London SW 1 Y 5 DB, UK, 1988. 1988.

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

Page 30 of 37

30

TOC Graphic

ACS Paragon Plus Environment

Page 31 of 37 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

Doped silicon clusters addition and substitution scheme, showing some doping possibilities. All the nonequivalent doping positions were considered as starting geometries of a given system. Red and blue colors refers to the dopant. 338x190mm (96 x 96 DPI)

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

Ground state geometries optimized of Sin-a-bAlaPb (n= 6, 7, 10 and 11; a, b =0, 1) with 6311++G(3df,3pd). Orange refers to the phosphorus atom and beige to the aluminum atom. 338x190mm (96 x 96 DPI)

ACS Paragon Plus Environment

Page 32 of 37

Page 33 of 37 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

Orbital energy diagram obtained with the RHF//M06/6-311++G(3df,3pd) method for the studied clusters. The red and blue colors refer to the contribution of the aluminum and phosphorous atomic orbitals to each molecular orbital. 169x131mm (300 x 300 DPI)

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

Electronic density difference between the doped and pure silicon clusters obtained through CASPT2//M06/6311++G(3df,3pd) with isosurface value of 0.003 e.bohr-3. Yellow and white circles refers to aluminum and phosphorus atoms, respectively. 338x190mm (96 x 96 DPI)

ACS Paragon Plus Environment

Page 34 of 37

Page 35 of 37 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

Electronic difference between the ground state and the first excited state of the doped silicon clusters obtained using CASPT2//M06/6-311++G(3df,3pd) with isosurface value of 0.003 e.bohr-3. 338x190mm (96 x 96 DPI)

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

Stability ranking function (ε3), applied on doped silicon clusters. Color black, red, blue and green refers to n =6, 7, 10 and 11, respectively. 169x130mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 36 of 37

Page 37 of 37 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

Fukui's functions calculated using electronic density difference between the neutral and ionics species obtained using CASPT2//M06/6-311++G(3df,3pd) with isosurface value of 0.003 e.bohr-3. Yellow and white circles refers to aluminum and phosphorus atoms, respectively. 338x190mm (96 x 96 DPI)

ACS Paragon Plus Environment