A brief introduction of the studied system:
Optimized geometry of small Al clusters
A series of Aluminum chain containing 2 to 9 atoms along with a triangular and a square shaped small clusters have been depicted in Fig. 1. All of these clusters are optimized structure with zero or more imaginary frequencies. For Al2 cluster interatomic distance is 2.531 Å, it is 2.472 Å for Al3 cluster, 2.728 Å for Al4 cluster, 2.734 Å for Al5 cluster, 2.738 Å for Al6 cluster, 2.693 Å for Al7 cluster, 2.723 Å for Al8 cluster, 2.721 Å for Al9 cluster. For triangular Al3 cluster, it is 2.509 Å and for square shaped Al4 cluster, it is 2.600 Å. Length of the Al chain with with number of atom is proportional which has been plotted in Fig. 2.
Now the turn to discuss the set of planar and quasi-planar Al13+ cluster whose optimized structure and electronic character have already been reported by us [19, 20]. Among these three clusters one is perfectly planar whereas two other are quasi-planar structures. Among these two quasi- At these structures, one atom has been replaced from any symmetry position by alkali atoms (Li, Na and K) [28]. Among, all these set of structures, doping at position 4 are the minimum energy structures which are depicted in Fig. S1.These structures have been chosen to see the plasmonic effect.
Time Dependent Density Functional Theory (Tddft):
TDDFT study for small Al clusters
TDDFT calculations provide a nice idea about TDM and OS of the structures. For, pure linear Al2 cluster, bright plasmons have been found at state 11 with TDM value and Oscillator strength are − 2.79 debye and 0.87 respectively (see Table 1).For, Al3 cluster this TDM and oscillator strength have been increased significantly to 4.75 Debye and 2.23 respectively (State 28). For, Al4 cluster, bright plasmons have been found at state 17 with TDM value and Oscillator strength are − 6.72 Debye and 1.95 respectively. In case of Al5 cluster, TDM and oscillator strength have also been enhanced significantly to -6.19 Debye and 3.31 respectively (State 55). For, Al6 cluster, bright plasmons have been found at state 33 with TDM value and Oscillator strength are − 6.89 Debye and 3.97 respectively. For, Al7 cluster this TDM and oscillator strength have been increased significantly to -7.48 Debye and 4.40 respectively (State 28). For, Al4 cluster, bright plasmons have been found at state 17 with TDM value and Oscillator strength are − 6.72 Debye and 1.95 respectively. Whereas, for Al5 this TDM and oscillator strength have also been enhanced to -6.19 debye and 3.31 respectively (State 55). So, there is a systematic variation of OS and TDM with Al chain has been plotted in Fig. 2.
Table 1
Table for highest plasmonics transition of small Al clusters obtained from TDDFT calculation.
Cluster | Length of the chain | Excited State | Optimized energy in Hartree | Absorbed wavelength in nm | Oscillator strength | TDM |
Linear Chain |
Al2 | 2.531 | 11 | -484.5209875 | 270.54 | 0.87 | -2.79 |
Al3 | 4.943 | 28 | -726.8232394 | 307.94 | 2.23 | 4.75 |
Al4 | 7.980 | 17 | -969.1440582 | 348.10 | 1.95 | -4.72 |
Al5 | 10.371 | 55 | -1211.4548099 | 351.94 | 3.31 | -6.19 |
Al6 | 12.784 | 33 | -1453.7481465 | 363.42 | 3.97 | -6.89 |
Al7 | 15.245 | 80 | -1696.0604769 | 386.67 | 4.40 | -7.48 |
Al8 | 17.636 | 49 | -1938.344392 | 395.60 | 4.37 | 7.54 |
Al9 | 20.057 | 114 | -2180.6589412 | 415.83 | 4.83 | -8.13 |
Triangle |
Al3 | 2.509 | 73 | -726.880578 | 182.43 | 1.14 | 6.8180 |
Square |
Al4 | 2.600 | 74 | -969.186197 | 173.90 | 0.98 | 5.5902 |
Tddft Study For Undoped And Doped Al Clusters
The TDDFT study reveals that for the undoped planar Al13+ cluster the bright plasmonic state is the state 189 with TDM value 4.37 Debye and Oscillator strength 2.59. Details of this transition have been jotted down in Table 2 and Table 3. For quasi-planar cluster, this TDM and oscillator strength has been reduced significantly to 2.27 Debye and 0.68 respectively which is observed for the state 185. Alkali atom doping has also significant effect on the TDM and OS and generally the TDM and OS decreases upon doping. For the Li doped planar Al13+ cluster state 174 correspond to the bright plasmon state which has TDM value and Oscillator strength 3.07 Debye and 1.17 respectively. For Li doped quasi-planar cluster this TDM and oscillator strength have been reduced significantly to 1.17 Debye and 0.72 respectively (State 188). For, planar Na doped Al13+ cluster, bright plasmons have been found at state 178 with TDM value and Oscillator strength are 2.65 Debye and 0.88 respectively. Similarly, for Na doped quasi-planar cluster this TDM and oscillator strength have been reduced significantly to 3.07 debye along Y-axis and 1.30 respectively (State 196). For, planar K doped Al13+ cluster, bright plasmons have been found at state 179 with TDM value and Oscillator strength are 3.25 Debye and 1.29 respectively. But in case of K doped quasi-planar cluster this TDM and oscillator strength have been enhanced to 3.29 Debye along Y-axis (-0.46 Debye along z-axis due to quasi-planarity) and 1.70 respectively (State 202). Detail study of all plasmonic parameters for different clusters has been given in Table 4. Undoped Transition state Al13+ cluster exhibits more plasmonic character than planar cluster but less plasmonic behaviour than quasiplanar cluster (Table S1). In each cases, major contribution of transition dipole moment has been found along Y-axis, in which bilateral symmetry axis of each cluster passes. But, in case of quasi-planar clusters these major contribution of TDM value along Y- axis is shared towards Z-axis due to non-planarity of the cluster. Naturally, dipole moment along a specific direction (Y-axis/ symmetry axis) has been reduced significantly, in case of corresponding quasi planar clusters. Oscillator strength which is proportional to TDM is also reduces, it indicate plasmonicity also reduces.
Table 2
Table for highest plasmonics transition of planar Al13+ clusters obtained from TDDFT calculation.
Clusters | State | Transition Energy ∆E in eV | Transition (Contributions) | ∆n | Transition Weight | Electric Transition Dipole Moments (TDM) | Oscillator Strength(f) |
µx | µy | µz |
Planar Clusters |
Pure Al13+ Cluster | 189 | 5.5253 | H -> L + 25 (21%) | 26 | 0.32205 | 0.0004 | 4.3749 | 0.0000 | 2.5909 |
H-5-> L + 14 (11%) | 20 | 0.23239 |
H-12 -> L (7%) | 13 | 0.18053 |
H-3 ->L + 20 (7%) | 24 | 0.18526 |
H-3 ->L + 23 (6%) | 27 | 0.17912 |
H-5 ->L + 19 (4%) | 25 | -0.13323 |
H-7 -> L + 10 (4%) | 18 | -0.13308 |
H-2 ->L + 22 (3%) | 25 | -0.12511 |
H-5 ->L + 17 (3%) | 23 | 0.11841 |
H-2 ->L + 19 (2%) | 22 | -0.10041 |
Li doping | 174 | 5.0860 | H-8->L + 6 (11%) | 15 | 0.23204 | 3.0729 | 0.0020 | 0.0000 | 1.17 |
H-7->L + 12 (25%) | 20 | 0.11406 |
H-6->L + 14 (16%) | 21 | 0.35035 |
H-8->L + 8 (3%) | 17 | -0.15449 |
H-6->L + 10 (5%) | 17 | 0.27844 |
H-6->L + 16 (5%) | 23 | -0.15522 |
H-3->L + 21 (3%) | 25 | -0.12555 |
H->L + 26 (8%) | 27 | 0.20037 |
Na doping | 178 | 5.0972 | H-7->L + 12 (23%) | 20 | 0.33977 | 0.0001 | 2.6508 | 0.0000 | 0.88 |
H-6->L + 14 (11%) | 21 | 0.23146 |
H-6->L + 16 (34%) | 16 | 0.41436 |
H-3->L + 21 (4%) | 16 | -0.14426 |
H->L + 25 (4%) | 26 | -0.14494 |
K doping | 179 | 4.9901 | H-8->L + 6(14%) | 15 | 0.26760 | 0.0002 | 3.2478 | 0.0000 | 1.29 |
H-6->L + 16 (25%) | 23 | 0.35385 |
H-7->L + 14 (6%) | 22 | -0.17591 |
H-6->L + 11 (4%) | 18 | -0.14795 |
H-6->L + 18 (3%) | 25 | 0.13107 |
H-4->L + 17 (5%) | 22 | 0.15866 |
H-3->L + 10 (2%) | 14 | -0.10910 |
H-3->L + 21 (6%) | 25 | 0.16746 |
H-2->L + 18 (2%) | 21 | -0.10478 |
H->L + 27 (4%) | 28 | -0.13594 |
H->L + 35 (4%) | 36 | 0.14602 |
Table 3
Table for highest plasmonics transition of quasi-planar Al13+ clusters obtained from TDDFT calculation.
Clusters | State | Transition Energy ∆E in eV | Transition (Contributions) | ∆n | Transition Weight | Electric Transition Dipole Moments (TDM) | Oscillator Strength(f) |
µx | µy | µz |
Quasi-planar Clusters |
Pure Al13+ Cluster | 185 | 5.3768 | H-1->L + 22 (18%) | 24 | 0.30168 | 2.2707 | -0.0011 | -0.0668 | 0.68 |
H-9 -> L + 7 (17%) | 17 | -0.29278 |
H-8 -> L + 11 (9%) | 20 | -0.20618 |
H-2 -> L + 21 (8%) | 24 | -0.19914 |
H-2 ->L + 20 (5%) | 23 | 0.15992 |
H-5 ->L + 19 (4%) | 25 | 0.14463 |
H-3 ->L + 19 (3%) | 23 | -0.12052 |
H -> L + 23 (3%) | 24 | 0.12147 |
Li doping | 188 | 5.6113 | H-4->L + 20 (31%) | 25 | -0.10528 | 0.0169 | -2.3003 | -0.0038 | 0.72 |
H-11->L + 5 (7%) | 17 | 0.19042 |
H-7->L + 13 (6%) | 21 | 0.17952 |
H-6-> L + 15 (2%) | 22 | -0.10981 |
H-6-> L + 17 (9%) | 24 | 0.21201 |
H-5-> L + 18 (4%) | 24 | -0.14332 |
H-5-> L + 20 (2%) | 26 | 0.39438 |
H-3-> L + 19 (3%) | 23 | 0.12151 |
Na doping | 196 | 5.6381 | H-10->L + 7 (13%) | 18 | -0.25106 | -3.0653 | 0.0004 | 0.0699 | 1.30 |
H -2-> L + 21 (16%) | 24 | -0.21331 |
H-7->L + 14 (9%) | 22 | -0.14356 |
H-6->L + 20 (4%) | 27 | -0.11638 |
H-5->L + 14 (3%) | 20 | 0.11596 |
H-5->L + 18 (3%) | 24 | -0.12056 |
H-5->L + 21 (3%) | 27 | -0.12384 |
H-4->L + 18 (3%) | 23 | 0.27847 |
H-1->L + 25 (3%) | 27 | -0.12663 |
H-1->L + 26 (3%) | 18 | -0.11838 |
K doping | 202 | 5.6072 | H-11->L + 2 (3%) | 14 | -0.11799 | 3.4921 | 0.0000 | -0.4644 | 1.70 |
H-10->L + 7 (6%) | 18 | 0.16918 |
H-8 -> L + 8 (9%) | 17 | 0.21429 |
H-5-> L + 14 (4%) | 20 | 0.14802 |
H-5-> L + 18 (4%) | 24 | -0.13577 |
H-5->L + 21 (5%) | 27 | 0.16383 |
H-4->L + 18 (2%) | 23 | 0.10386 |
H-1->L + 27 (8%) | 29 | -0.19779 |
H-> L + 26 (4%) | 27 | -0.13618 |
Table 4
Variation of absorbed wavelength at resonance with dimension of the clusters.
Planar/ Quasi-planar | Pure/Doping | Dimensions | States at maximum oscillator strength | Absorbed wavelength at resonance in nm | Absorbed energy at resonance in eV | Oscillator strength |
Max horizontal length in Å | Max vertical length in Å | Average Length in Å |
Planar | Pure Al13+ Cluster | 7.95 | 7.53 | 7.74 | 189 | 224.39 | 5.53 | 2.59 |
Li doping | 7.83 | 8.96 | 8.39 | 174 | 243.78 | 5.09 | 1.17 |
Na doping | 7.83 | 9.38 | 8.60 | 178 | 243.24 | 5.10 | 0.88 |
K doping | 7.81 | 10.03 | 8.92 | 179 | 248.46 | 4.99 | 1.29 |
Quasi-planar | Pure Al13+ Cluster | 7.51 | 7.36 | 7.43 | 185 | 230.59 | 5.38 | 0.68 |
Li doping | 6.92 | 6.68 | 6.80 | 188 | 220.95 | 5.61 | 0.73 |
Na doping | 7.08 | 6.69 | 6.88 | 196 | 219.90 | 5.63 | 1.30 |
K doping | 7.11 | 6.72 | 6.91 | 202 | 221.11 | 5.60 | 1.70 |
Uv-vis Spectra Study:
UV-Vis spectra for small Al clusters
UV-Vis spectra [29] of these three clusters have been depicted in Fig. 3. For linear Al2 cluster three consecutive distinct sharp resonance peaks have been observed, among them peak located at 270.54 nm with oscillator strength 0.87 is highest. Whereas for Al3 cluster, it is located at 307.94 nm which exhibits significant oscillator strength of 2.23. For, Al4 chain, a sharp resonance peak has been observed at 348.10 nm with oscillator strength 1.95. Whereas, for Al5 cluster, it is located at 351.94 nm with oscillator strength 3.31. For Al6 chain, this sharp resonance peak has been observed at 363.42 nm with oscillator strength 3.97. Whereas, for Al7 cluster, it is located at 386.67 nm with oscillator strength 4.40. For Al8 cluster, this sharp resonance peak has been observed at 395.60 nm with oscillator strength 4.97. Whereas, for Al9 cluster, it is found at 415.83 nm having oscillator strength 4.83. Now turns for 2d clusters. For triangular Al3 cluster, highest plasmonic emission is observed at 182.43 nm with oscillator strength 1.14. For square shaped Al4 cluster, it is found at 173.90 nm with oscillator strength 0.98. So, resonance absorbance coefficient which is directly proportional to number of plasmons associated with this resonance or oscillator strength is higher for planar clusters than its corresponding quasi-planar cluster. Also, resonance frequency is higher for planar clusters than its corresponding quasi-planar clusters. So, for these pure and doped Aluminum clusters, the plasmon band is located in the ultraviolet region. Generally, for usually available Ag, Au metal nanoparticles, it is seen within IR region. So, resonance absorbance frequency of Al nanoparticle is high. Also, this resonance absorption frequency can be finely tuned by distorting and doping the clusters.
UV-Vis spectra for Al 13 + clusters.
UV-Vis spectra of planar clusters and quasi-planar clusters have been depicted in Fig. 4. For pure planar Al13+ cluster, sharp resonance peak has been observed at 224.39 nm with oscillator strength 2.59.Whereas, for undoped true minima quasi-planar Al13+ cluster, it is located at 230.59 nm with oscillator strength 0.6798. For undoped transition state Al13+ cluster, sharp resonance peak has been observed at 228.21 nm with oscillator strength 1.10 (Fig. S2). For planar Li doped Al13+ cluster, this resonance peak has been observed at 243.78 nm with oscillator strength 1.17. Whereas, for quasi-planar Li doped Al13+ cluster, it is located at 220.95 nm with oscillator strength 0.72. For planar Na doped Al13+cluster, this sharp resonance peak has been observed at 243.24 nm with oscillator strength 0.88. Whereas, for quasi-planar Na doped Al13+ cluster, it is located at 219.90 nm with oscillator strength 1.30. For planar K doped Al13+cluster, a sharp resonance peak has been observed at 229.70 nm with oscillator strength 1.8362. Whereas, for quasi-planar K doped Al13+ clusters, it is located at 221.11 nm with oscillator strength 1.70. On the other hand, resonance absorbance coefficient which is directly proportional to number of plasmons associated with this resonance is higher for planar clusters than its corresponding quasi-planar clusters. Also, resonance frequency is higher for planar clusters than its corresponding quasi-planar clusters. So, for these pure and doped Aluminum clusters, the plasmon band is located in the ultraviolet region. Generally, for usually available Ag, Au metal nanoparticles, it is seen within IR region. So, resonance absorbance frequency of Al nanoparticle is high. Also, this resonance absorption frequency can be finely tuned by restructuring and doping the clusters.
Canonical Molecular Orbitals:
Canonical Molecular Orbital for small Al clusters
Canonical molecular orbitals [30] for Al chain with Triangular and square shaped Al clusters have been shown in Fig. 5. In case of Al2 and Al3, HOMO is delocalized due to \(\pi\) electrons but LUMO is localized due to\(\sigma\)electrons. In case of Al4 system, both HOMO and LUMO are delocalized due to \(\pi\) electrons. In Al5 cluster, \(\pi\) electron contributes to HOMO and it is delocalized, but is localized by \(\sigma\) electron. In case of Al6 and Al7 clusters, HOMO is delocalized due to \(\pi\) electron, whereas LUMO are localized due to \(\sigma\) electron. In case of Al8 and Al9 clusters, both localized HOMO and LUMO have contribution to \(\pi\) electron. In case of triangular Al3 and Al4 clusters, HOMO are delocalized and this comes due to \(\pi\) electron. LUMO is not delocalized and it has contribution to \(\sigma\) electron. So, most of the HOMO orbital are delocalized due to \(\pi\) electron may also the origin of exhibiting plasmonic character.
Canonical Molecular Orbital For Small Al Clusters
Canonical Molecular Orbital of planar and quasi-planar clusters has been shown in Fig. 6. Here, canonical molecular orbital of highest oscillator strength with greater transition weight have been depicted. For pure planar Al13+ cluster, delocalized orbitals of transition from HOMO to LUMO + 25 have probability 21%. In this case, delocalized \(\pi\) orbitals have been confined in a shorter region after transition. Next higher transition from HOMO-5 to LUMO + 14 has probability 11%. In this case, σ delocalized orbitals have been confined in a shorter region after transition. For planar Li doped Al13+ cluster, delocalized orbitals of transition from HOMO-7 to LUMO + 12 has probability 25%. In this case, contributions of both orbitals are due to \(\sigma\) electrons. Transition from HOMO-6 to LUMO + 16 have probability 11%. In this case, transition has been taken place from σ delocalized orbitals to another σ delocalized orbital. For planar Na doped Al13+ cluster, delocalized orbitals of transition from HOMO-6 to LUMO + 16 has probability 34%. In this case, a small contribution of \(\pi\) orbitals is found after transition from delocalized σ orbital. Next higher transition from HOMO-7 to LUMO + 12 has probability 23%. In this both cases, contribution of delocalized σ orbital is found. For planar K doped Al13+ cluster, delocalized orbital of transition from HOMO-6 to LUMO + 16 has probability 25%. In this case, transition occurs from delocalized σ orbital to \(\pi\) orbitals. Next higher transition from HOMO-8 to LUMO + 6 has probability 14%.Now the turns come for quasi-planar clusters. For pure quasi-planar Al13+ cluster, delocalized orbitals of transition from HOMO-1 to LUMO + 22 has probability 18%. In this case, transition has been taken place from delocalized \(\pi\) orbitals to σ orbitals. Next higher transition from HOMO-9 to LUMO + 7 has probability 17%. In this case, area of σ delocalized orbitals have been confined in a shorter region after transition. For quasi-planar Li doped Al13+ cluster, delocalized orbitals of transition from HOMO-4 to LUMO + 20 has probability 31%. In this case, transition has been taken from delocalized \(\pi\)orbitals \(\sigma\) orbitals.Transition from HOMO-6 to LUMO + 17 have probability 9%. In this case, delocalized \(\pi\) orbitals have been confined in a shorter region after transition from delocalized σ orbitals. For quasi-planar Na doped Al13+ cluster, delocalized orbitals of transition from HOMO-2 to LUMO + 21 has probability 16%. In this case, transition occurs from delocalized \(\sigma\) orbitals to similar σ orbitals. Next higher transition from HOMO-10 to LUMO + 7 has probability 13%. In this case, delocalized \(\sigma\) orbitals have been confined in a shorter region after transition. For quasi-planar K doped Al13+ cluster, delocalized orbitals of transition from HOMO-8 to LUMO + 8 have probability 9%. In this case, delocalized \(\sigma\) orbitals have been confined in a shorter region after transition. Next higher transition from HOMO-1 to LUMO + 27 has probability 8%. In this case, transition occurs from delocalized \(\sigma\) orbitals to similar σ orbitals. For undoped transition state Al13+ cluster, delocalized orbitals of transition from HOMO-5 to LUMO + 19 have probability 23%. In this case, delocalized \(\sigma\) orbitals have been confined in a shorter region after transition. Next higher transition from HOMO-4 to LUMO + 18 has probability 17%. In this case, transition occurs from delocalized \(\sigma\) orbitals to another σ orbitals. So, most of the contribution of electronic transition which is generally plasmonic transition is found in delocalized orbitals.
Polarizability Study:
Polarizability is a measure of how easily an electron cloud is distorted by an electric field [31]. Generally, the electron cloud will being to an atom or molecule or ion. But, in this case electron cloud is due to plasmonic nanopaticle of Al13+ nanoclusters. It can be represented in tensor or matix form for any arbitrary body. Sometimes due to symmetry of the system, its off-diagonal vanish. This is called anisotropic polarizability, it is represented by-
αiso=\(\frac{1}{3}\) (αxx + αyy + αzz)
This expression is applicable in our planar system.
However, during the dynamics, when the planar clusters transform into quasi-planar clusters, the nuclear displacements due to the activated bending mode cause a change in the polarizability that cannot be described in terms of its diagonal elements. The plane deviates from the planar symmetry and for this reason the off-diagonal polarizability contributions have to be included in the molecular polarizability. The anisotropic term, including the off-diagonal terms, can be represented by
αaniso= \(\frac{1}{\sqrt{2}}\sqrt{{({{\alpha }}_{\text{x}\text{x}}-{{\alpha }}_{yy})}^{2}+{({{\alpha }}_{\text{x}\text{x}}-{{\alpha }}_{zz})}^{2}+{({{\alpha }}_{yy}- {{\alpha }}_{zz})}^{2}+6({{\alpha }}_{xy}^{2}+{{\alpha }}_{\text{y}\text{z}}^{2}+{{\alpha }}_{zx}^{2})}\)
Polarizability of these reported undoped and doped clusters has been provided in Table 5. For pure planar Al13+ cluster, y-component of polarizability (αyy )is largest (837.28 unit) and z-component of polarizabilty (αzz ) is smallest (279.04 unit).There is no cross/off-diagonal components of polarizability. Isotropic and anisotropic polarizabilities are respectively 642.99 and 546.34 units. For pure quasi-planar Al13+ cluster, y-component of polarizability (αyy )is largest(788.16 unit) and z-component of polarizabilty (αzz ) is smallest (332.14 unit).There also exist little amount of cross/off-diagonal components of polarizability due to quasi-planarity and principal component of polarizability reduces simultaneously. Isotropic and anisotropic polarizabilities are respectively 641.75 and 492.86 units. For undoped transition state Al13+ cluster, y-component of polarizability (αyy )is largest (802.84 unit) and z-component of polarizabilty (αzz ) is smallest (310.02 unit). Like quasi-planar true minima cluster, there also exist little amount of cross/off-diagonal components of polarizability due to quasi-planarity and principal component of polarizability reduces simultaneously. Isotropic and anisotropic polarizabilities are respectively 645.19 and 507.98 units (Table S2). For planar Li doped Al13+ cluster, y-component of polarizability (αyy ) is largest (864.39 unit) and z-component of polarizabilty (αzz ) is smallest (271.61 unit).There is no cross/off-diagonal components of polarizability. Isotropic and anisotropic polarizabilities are respectively 644.70 and 378.14 units. For Li doped quasi-planar Al13+ cluster, y-component of polarizability (αyy ) is also largest (677.51 unit) and z-component of polarizabilty (αzz ) is smallest (364.02 unit).There also exist little amount of cross/off-diagonal components of polarizability due to quasi-planarity and eventually principal component of polarizability reduces simultaneously. Isotropic and anisotropic polarizabilities are respectively 568.58 and 400.74 units. For planar Na doped Al13+ cluster, y-component of polarizability (αyy) is largest (907.09 unit) and z-component of polarizabilty (αzz ) is smallest (277.41 unit).There is no cross/off-diagonal components of polarizability. Isotropic and anisotropic polarizabilities are respectively 663.85 and 828.87 units. For Na doped quasi-planar Al13+ cluster, y-component of polarizability (αyy` ) is largest (693.96 unit) and z-component of polarizabilty (αzz ) is smallest (370.51 unit).There also exist little amount of cross/off-diagonal components of polarizability due to quasi-planarity and eventually principal component of polarizability reduces simultaneously. Isotropic and anisotropic polarizabilities are respectively 580.00 and 256.29 units. For K doped planar Al13+ cluster, y-component of polarizability (αyy ) is largest (706.30 unit) and z-component of polarizabilty (αzz ) is smallest (379.80 unit).There also exist little amount of cross/off-diagonal components of polarizability due to quasi-planarity and eventually principal component of polarizability reduces simultaneously. Isotropic and anisotropic polarizabilities are respectively 591.47 and 345.68 units. For K doped quasi-planar Al13+ cluster, y-component of polarizability (αyy ) is largest (706.30 unit) and z-component of polarizabilty (αzz ) is smallest (379.80 unit).There also exist little amount of cross/off-diagonal components of polarizability due to quasi-planarity and eventually principal component of polarizability reduces simultaneously. Isotropic and anisotropic polarizabilities are respectively 591.47 and 345.68 units. These interesting result also nicely explains the plasmonic responds of the clusters. As these clusters has larger vertical dimension and bilateral symmetry axis passes through y-axis, polarizabilities of all the clusters have been found along y-axis. Due to transition from planar structure to quasi-planar structure its off-diagonal or cross component of polarizability (αxy or αyz or αzx )is enhanced and principal component (αxx or αyy or αzz ) of polarizability is discreased. It is found due to charge separation of the cluster from x-y plane.
Table 5
Polarizabilities for pure planar and quasi-planar Al13+ clusters.
Clusters | Pure/ Doping | αxx | αxy | αyy | αxz | αyz | αzz | αiso | αaniso |
Planar | Pure Al13+ cluster | 812.64 | -0.01 | 837.28 | 0.00 | 0.00 | 279.04 | 642.99 | 546.34 |
Li doping | 798.11 | -0.02 | 864.39 | 0.00 | 0.00 | 271.61 | 644.70 | 378.14 |
Na doping | 807.05 | -0.01 | 907.09 | 0.00 | 0.00 | 277.41 | 663.85 | 828.87 |
K doping | 814.58 | -0.02 | 929.39 | 0.00 | 0.00 | 283.14 | 675.70 | 597.18 |
Quasi-planar | Pure Al13+ cluster | 804.96 | -13.10 | 788.16 | -82.47 | 45.08 | 332.14 | 641.75 | 492.86 |
Li doping | 664.21 | 9.32 | 677.51 | -131.04 | 69.60 | 364.02 | 568.58 | 400.74 |
Na doping | 675.53 | 13.33 | 693.96 | -129.70 | 68.80 | 370.51 | 580.00 | 256.29 |
K doping | 688.32 | 12.98 | 706.30 | -130.47 | 69.21 | 379.80 | 591.47 | 345.68 |
Density Of States:
Plot of density of state (Fig. 7) is also the evidence of permanent electric dipole moment (which has also been found from Tables 2 & 3) having larger Transition dipole moment (TDM) [32]. These plots for all these structures show no mirror symmetry between valance band and conduction band - the states within valance band and conduction band are not symmetrical. Also, none of the dotted line passes through 0 eV energy. Interestingly, these dotted line where it intersects on x-axis i.e. energy axis provides the energy corresponds to the excitation energy with highest oscillator strength (Tables 2 & 3). For planar Al13+ cluster it is 5.53 eV, for Li doped cluster it is 5.09 eV, for Na doped cluster it is 5.10 eV, for K doped cluster it is 4.99 eV. Similarly, for quasi-planar Al13+ cluster it is 5.37 eV, for Li doped cluster it is 5.61 eV, for Na doped cluster it is 5.63 eV, for K doped cluster it is 5.60 eV. For transition state Al13+ cluster it is 5.42 eV (Fig. S4).