Next Article in Journal
Effect of Mineral Admixtures on the Mechanical and Shrinkage Performance of MgO Concrete
Previous Article in Journal
Review of the Interactions between Conventional Cementitious Materials and Heavy Metal Ions in Stabilization/Solidification Processing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Rational Design, Stabilities and Nonlinear Optical Properties of Non-Conventional Transition Metalides; New Entry into Nonlinear Optical Materials

1
Department of Chemistry, College of Science, King Faisal University, Al-Ahsa 31982, Saudi Arabia
2
Chemical Sciences, Faculty of Science, Universiti Brunei Darussalam, Jalan Tungku Link, Gadong BE1410, Brunei
3
University of Rostock, Institute of Chemistry, Physical and Theoretical Chemistry, Albert-Einstein-Straße 27, 18059 Rostock, Germany
4
Department of Science and Technology of Life, Light and Matter, Faculty of Interdisciplinary Research, University of Rostock, 18059 Rostock, Germany
5
Leibniz-Institut für Katalyse an der Universität Rostock, Albert-Einstein-Strasse 29a, 18059 Rostock, Germany
6
Department of Chemistry, COMSATS University Islamabad, Abbottabad Campus, Abbottabad 22060, KPK, Pakistan
*
Authors to whom correspondence should be addressed.
Materials 2023, 16(9), 3447; https://doi.org/10.3390/ma16093447
Submission received: 16 February 2023 / Revised: 14 April 2023 / Accepted: 20 April 2023 / Published: 28 April 2023

Abstract

:
Electronic and nonlinear optical properties of endohedral metallofullerenes are presented. The endohedral metallofullerenes contain transition metal encapsulated in inorganic fullerenes X12Y12 (X = B, Al & Y = N, P). The endohedral metallofullerenes (endo-TM@X12Y12) possess quite interesting geometric and electronic properties, which are the function of the nature of the atom and the size of fullerene. NBO charge and frontier molecular orbital analyses reveal that the transition metal encapsulated Al12N12 fullerenes (endo-TM@Al12N12) are true metalides when the transition metals are Ni, Cu and Zn. Endo-Cr@Al12N12 and endo-Co@Al12N12 are at the borderline between metalides and electrides with predominantly electride characteristics. The other members of the series are excess electron systems, which offer interesting electronic and nonlinear optical properties. The diversity of nature possessed by endo-TM@Al12N12 is not prevalent for other fullerenes. Endo-TM@Al12P12 are true metalides when the transition metals are (Cr-Zn). HOMO-LUMO gaps (EH-L) are reduced significantly for these endohedral metallofullerenes, with a maximum percent decrease in EH-L of up to 70%. Many complexes show odd–even oscillating behavior for EH-L and dipole moments. Odd electron species contain large dipole moments and small EH-L, whereas even electron systems have the opposite behavior. Despite the decrease in EH-L, these systems show high kinetic and thermodynamic stabilities. The encapsulation of transition metals is a highly exergonic process. These endo-TM@X12Y12 possess remarkable nonlinear optical response in which the first hyperpolarizability reaches up to 2.79 × 105 au for endo-V@Al12N12. This study helps in the comparative analysis of the potential nonlinear optical responses of electrides, metalides and other excess electron systems. In general, the potential nonlinear optical response of electrides is higher than metalides but lower than those of simple excess electron compounds. The higher non-linear optical response and interesting electronic characteristics of endo-TM@Al12N12 complexes may be promising contenders for potential NLO applications.

1. Introduction

Alkalides are the compounds where anionic positions are occupied by alkali metal anions [1]. Alkalides are the focus of research in the current decade due to their fascinating properties including the theoretically predicted large nonlinear optical response [2,3,4,5,6,7,8,9,10]. The concept of alkalides has been well known for more than five decades now both in solution [11,12,13,14] and crystalline phase [15,16]; however, the room-temperature stable alkalides were reported by Dye and co-workers in 1999 [17]. A number of reports are available in the literature which illustrate the large nonlinear optical response of alkalides. For example, Wei Chen et al. [18] reported that the hyperpolarizability of AdzLi+Na is remarkably high (βo = 6.20681 × 105 au). Similarly, an increase in the nonlinear optical response of Li+(calix[4]pyrrole)M (M = Li, Na and K) is also reported by the introduction of alkalides [7]. The hyperpolarizability values of these alkalides are larger than their corresponding electrides [7]. Similarly, the alkalides from organic-cage-like compounds are also reported where alkali metal cation resides inside the cage, whereas the alkalide resides outside the cage [7]. The nonlinear optical response of these alkalides [9] is even higher than the nonlinear optical response of inverse metal hydride and Li+(calix[4]pyrrole)M (M = Li, Na and K). More recently, Janus alkalides are reported from facially polarized Janus all-cis-1,2,3,4,5,6-hexafluorocyclohexane. The hyperpolarizability values of these alkalides [10] (calculated at MP2/6-31+G(d,p) reach up to 1.45 × 106 au.
All of the alkalides reported in the literature share a common design principle. A complexant interacts with two atoms (either both alkali metal atoms or an alkali and electropositive atom) where one atom receives a positive charge (alkali or other electropositive atom), while the second alkali metal atom receives a negative charge. For example, in the reverse metal hydride [18,19], a hydrogen atom receives a positive charge upon interaction with nitrogen atoms, whereas the alkali metal atom receives negative charge. For a complexant having coordination sites available, the alkali metal entrapped in the coordination sphere is positively charged, whereas the free alkali metal atom gains a negative charge. The only example of an alkalide where the alkali metal cation is not bound in a cavity is of Janus alkalide [6,10], where the alkali metal cation and anion are present on the open faces of a molecule. These structural variations have a significant impact on the nonlinear optical properties. For example, the βo value of Li+(NH3)4M with a cage-like complexant (NH3)4 is about six times larger than that of the corresponding Li+(calix[4]pyrrole)M- with cup-like complexant calix[4]-pyrrole [4]. Moreover, the nonlinear optical response of the Janus alkalide is much larger than other cage-like alkalides. Similarly, borates consisting of boron oxyanions are a fast-growing field of solid-state chemistry with well-established and important functionality [20]. In recent decades, freshly discovered borates gained the spotlight in research due to the great variety of their molecular structures (0D, 1D, 2D and 3D) [21,22].
Despite these advancements in the rational design of alkalides, one does not find any parallel rational design strategy for transition metal anions. The literature reveals a number of molecules/complexes where transition metals exist as anions. However, for transition metal anions, the difference in electropositivity of combining atoms is generally taken as the guiding principle. For example, platinum and gold atoms in Cs2Pt [23] and CsAu [24,25,26,27] bear a −2 and −1 charge, respectively. In these compounds, the transition metal atom acts as a reservoir to hold 12 electrons in 6s2 5d10 configuration for Pt2− and Au1−. A similar situation is observed for semiconducting oxides containing 18 electron octahedral cations [PtIn6]10+ and [IrIn6]9+ [28,29]. A quick glance at the structure of these cations reveals that the platinum and iridium in these complexes are in −2 and −3 oxidation states.
From the discussion above, one can clearly see that the transition metal anions are all covalent compounds where transition metal is covalently bound to some other electropositive atom (vide supra). Moreover, a design strategy of the simple coordination compound containing the transition metal anion is lacking. The only example reported in the literature for a coordination-based metalide is by Sun and co-workers, who have recently shown that the judicious placement of coinage metal and alkali metals on opposite faces of all-cis-1,2,3,4,5,6 hexfluorcyclohexane can generate metalides with a enhanced nonlinear optical response [30]. Moreover, density functional theory (DFT) calculations have also been extensively reported in the literature on the geometric, thermodynamic stability and electronic properties of heterofullerenes for nonlinear optical response [31,32].
With these targets, a new design strategy is required. If the structures of the alkalides is analyzed, one can clearly see that the alkali metal atom inside the cage gains a positive charge because the positive charge is balanced by an electron donation from electron-rich nitrogen or oxygen atoms. It is proposed here that if the entrapped alkali atom is surrounded by electropositive atom, then it may receive a negative charge without the assistance of any other compensating alkali metal cation. In this regard, the molecules of choice are inorganic fullerenes such as X12Y12 where X = Al or B and Y = N or P. These inorganic fullerenes have electropositive and electronegative atoms available. Although these contain electronegative atoms present in the sphere, they also have electropositive atoms. If the effect of electropositive atoms dominates over the electronegative atoms, then the encapsulation of transition metals in these fullerenes may provide an easy approach to generate metalides. The diversity of atoms present in these fullerenes may help in finding a suitable candidate in this regard.
These inorganic fullerenes can exist in a variety of sizes such as X12Y12, X16Y16, X28Y28, etc. Theoretical studies illustrate that X12Y12 nanoclusters are more stable than higher analogues [33,34]. X12Y12 nanocages are the subject of extensive exploration in the recent past for their application in sensors [35,36,37], storage [38,39], batteries and nonlinear optical materials [40,41,42,43,44]. Several studies reveal that the doping of alkali metal, whether exohedral or substitutional, increases the hyperpolarizability of the nanocage [40,41,43,44,45,46,47,48,49]. Xuri Huang and coworkers have shown that the exohedral doping of alkali metal significantly enhances the hyperpolarizability of Al12N12 nanocages. The authors argued that diffuse excess electrons are responsible for boosting the hyperpolarizability. Similarly, Maria et al. have shown that substitutional doping of alkali metal for X or Y increases the hyperpolarizability by two to three orders of magnitude [40,41,42]. However, this increase is less when compared with those for exohedral doping. Endohedral doping of alkali metals in these fullerenes is also reported for their nonlinear optical response. However, the alkalide characteristic of these EMFs is not discussed. The focus of the current study is to examine the metalide character of transition-metal-encapsulated fullerenes. Encapsulation of the metal atom imparts some unique characteristics to the nanocages, which are non-existent in pure nanocages. Endohedral doping of nanocages have found possible applications in many areas of scientific research including renewable energy, photonics and optoelectronics [50,51,52]. Moreover, the hyperpolarizability of EMF is discussed for the first time. A complete systematic spin-polarized DFT study is carried out.

2. Methodology

All calculations are performed with Gaussian 09 [53], whereas the visualization of the results was achieved through GaussView [54]. The geometry of all structures were optimized without any symmetry constraints at B3LYP/6-31G(d,p). It has been shown by us and others that the geometries of these nanocages are predicted reliably with the hybrid B3LYP method [35,55,56,57,58]. Moreover, the geometries are not affected to any appreciable extent when dispersion-corrected methods (ωB97XD and M05-2X) are applied. For transition metals, the spin-polarized DFT study was carried out in which the four lowest possible spin states are evaluated. The geometries of these endohedral metallofullerenes were confirmed at true minima through frequency analysis (lack of any imaginary frequency). The charge analysis was performed through natural bonding orbital (NBO) approach to realize the metalide nature of these complexes. S2 values for all complexes are evaluated to realize any spin contamination. The binding energies are calculated for these complexes in order to evaluate their stabilities. The binding energies are calculated through the following equation.
Eb = EM@X12Y12 − (EX12Y12 + EM)
where EX12Y12, EM and EM@X12Y12 are the total energies of the undoped X12Y12 molecule, the metal atom M and the corresponding doped system M@X12Y12, respectively. The bond order is also calculated through NBO approach for all these complexes.
Since polarizability is a measure of the distortion of the electronic cloud under the influence of an applied electric field, the relationship between the energy of a system in a weak and homogenous electric field is as follows.
E = EoµαFα − 0.5ααβFαFβ − 1/6βαβγFαFβFγ − …
where Fα is the electric field component along the α direction and Eo is the molecular total energy without electric field. µα, ααβ, and βαβγ are the dipole, the polarizability and the first hyperpolarizability, respectively.
Hyperpolarizability is calculated at the long-range-corrected CAM-B3LYP and LC-BLYP methods because the conventional DFT methods overestimate the hyperpolarizability due to incorrect electric field dependence provided by exchange correlation functional. The CAM-B3LYP method provides long range correction through Coulomb attenuation method, and it is superior for calculating NLO properties.
The mean polarizability and hyperpolarizability are described as
α = 1/3(αxx + αyy + αzz)
β0 = (βx2 + βy2 + βz2)1/2
βvec, which is more relevant to a measurable quantity, is the projection of β on the dipole moment vector and is given by the following equation [59]:
β v e c = μ i β i     μ
where µi, βi and |µ| are the dipole moments along the i direction, polarizability along the i direction and the total dipole moment, respectively.
UV-Vis calculations are performed through the time-dependent DFT approach (TD-DFT) at LC-BLYP method (TD-LC-BYLP) along with 6-31+G(d) basis set. A total of twenty states were computed with the 50:50 singlet triplet. These computations were performed with unrestricted formalism. The well-parametrized Bethe–Salpeter Green’s function approach has an advantage that it is parameter free and has been shown to provide reliable results in a large variety of organic or inorganic systems, with an insulating or metallic character, in extended (bulk) system or finite (reduced) size systems [60,61,62,63,64]. However, the Green’s function formalism cannot be considered to provide higher accuracy than DFT/TDDFT simulations [65].
A two-level model proposed by Oudar and Chemla [66] is also applied to understand the contributing factors for hyperpolarizability.

3. Results and Discussion

3.1. Spin State Analysis

For transition-metal-doped systems, it is better to discuss the spin states prior to discussing the optimized geometries. As presented above, the geometries of transition-metal-doped nanocages are optimized for the four lowest possible spin states. These include doublet, quartet, sextet and octet for Sc, V, Mn, Co and Cu. The spin states studied for Ti, Cr, Fe, Ni and Zn are singlet, triplet, quintet and septet. The results (Table 1) do not deliver any obvious trends at first glance because the most stable spin states are changing with the nature of the nanocage, except for late transition metals Ni, Cu and Zn, where stable spin states are singlet, doublet and singlet, respectively, regardless of the nature of the nanocage. The endohedral cobalt complexes endo-Co@X12Y12 are not very different where the most stable spin state is doublet for all nanocages, except for Al12N12 (endo-Co@Al12N12), where the most stable spin state is quartet. There are some other generalizations available for the data presented in Table 1. For example, transition metals prefer a higher spin in Al nanocages (Al12N12 and Al12P12). Additionally, the most stable spin states of Ti, V, Cr, Mn and Fe in Al12N12 nanocage are quintet, quartet, septet, sextet and quintet, respectively. The stable spin states of Ti, V, Cr, Mn and Fe in Al12P12 are quintet, sextet, septet, sextet and quintet, respectively.
The results in Table 1 also give information about the energy differences of higher spin states from the lowest energy spin state. The energy difference between the lowest energy spin state and the next spin state is significant in all cases except for endo-Co@Al12N12 and endo-Co@Al12P12, where the differences are 0.92 and 0.43 kcal mol−1, respectively. For Al12N12, small differences between the two lowest spin states are also observed for endo-Fe@Al12N12 and endo-Ni@Al12N12: 1.41 and 0.44 kcal mol−1, respectively. The highest difference in energy between the two states is observed for endo-Cu@X12Y12 (for all systems). These values are 44.17 and 45.35 kcal mol−1 for endo-Cu@Al12N12 and endo-Cu@Al12P12, respectively.
The symmetry group of both optimized nanocage Al12N12 and Al12P12 is Th. In the case of endo-TM@Al12N12, the Th is retained for Cr and Zn metal, whereas for Mn the symmetry observed is C3, while the rest of the complexes show C1 symmetry. Similarly, in case of endo-TM@Al12P12, the Th is retained for Mn and Zn metal, whereas C3 symmetry is observed for Ti and Ni. Moreover, C2V symmetry in obtained for Fe, while all the rest of the complexes exhibit C1 symmetry. Even when the symmetry is not C1 for the most stable spin state (as mentioned above), the higher spin states show C1 state in these cases, probably due to the distortion faced in the higher spin state.

3.2. Optimized Geometries

The optimized geometries of transition-metal-doped Al12N12 nanocages are shown in Figure 1. Sideview images are shown for the proper position of transition metal in the nanocage. The optimized geometries of all structures (Figure 1) indicate that most of the transition metals are not present at the center of the nanocages except for certain middle transition metal atoms (Cr, Mn and Fe) and Zn. These results are somewhat perplexing because early transition metal atoms are more electron deficient, and they tend to bind with high coordination number in order to overcome their electron deficiency. The central orientation of Zn in Al12N12 nanocage may be attributed to the lack of a binding tendency in Zn due to the completely filled d orbital, as it is also evident from binding energies analysis (vide infra). The change in the bonding orientation of different transition metals in the nanocage has a strong impact on the charges of transition metals (see NBO section).
An identical situation is observed for endo-TM@Al12P12 complexes, where early transition metals Sc, Ti and V tend to sit on a side of the nanocage, but Cr, Mn and Fe reside almost at the center of the nanocage (Figure 2). Co, Ni and Cu are also oriented towards a side of the nanocage, but Zinc is again present at the center of the nanocage. It is worth mentioning that the complexes where transition metal is not present at the center are not similar to each other. In some cases, the transition metal is facing towards a hexagon of the nanocage, whereas, in others, it is facing a tetragon. Among all complexes, endo-Zn@X12Y12 complexes only contain the metal atom at the center of the nanocages for all inorganic fullerenes, which mainly attributed to lack of any tendency for coordination due to completely filled d orbital.

3.3. Binding Energies

The binding energies of transition metals with the nanocages are analyzed (Table 2). Obviously, no trend could be traced in the interaction energies. The interaction of early transition metals with the Al12N12 is an exothermic process. The Gibbs free energies of the interaction of early transition metals show an odd–even effect. The interaction energy of Sc with Al12N12 is −12.56 kcal mol−1 increases to −41.72 for endo-Ti@Al12N12, but again shows decrease and increase to −8.78 and −13.95 kcal mol−1 for endo-Vi@Al12N12 and endo-Cr@Al12N12, respectively. The Gibbs free energies are positive for endo-Mn@Al12N12 and endo-Fe@Al12N12. The Gibbs free energies are 38.83 and 9.79 kcal mol−1 for endo-Mn@Al12N12 and endo-Fe@Al12N12, respectively. These are the complexes where transition metal resides at the center of the nanocage (vide supra). From the results, it is obvious that central orientation of transition metal in Al12N12 is because of the weak interaction of the metal with the nanocage, except for endo-Cr@Al12N12, where the interaction is exothermic despite the fact that the metal resides in the center of the nanocage. The geometric parameters revealed that Zn in endo-Zn@Al12N12 is also oriented in the center of the nanocage. The interaction energy is also positive for this complex which supports the above-mentioned notion. For these metal-encapsulated nanocages, two factors are operational: (a) attractive force between metal and nanocage and (b) repulsive interaction between nanocage and metal due to steric interaction. For weakly coordinating metals such as Mn and Zn, the repulsive interactions dominate over weak attractive forces, which results in overall endergonic process. For Al12N12, the strongest interaction is observed when nickel atom is encapsulated in the nanocage (−77.84 kcal mol−1). The binding of Cu with Al12N12 is also strongly exergonic by −38.17 kcal mol−1.
Al12P12 nanocage is larger in size; therefore, the repulsive interaction between nanocage and metal atoms is expected to be negligible. Analysis of the interaction energies reveals that the interaction of all transition metals with Al12P12 nanocage is favorable. There are some common features between the endo-TM@Al12N12 and endo-TM@Al12P12 series. For example, the interaction energies of endo-Mn@Al12P12, endo-Fe@Al12P12 and endo-Zn@Al12P12 are lower than other complexes, but fortunately, these interaction energies are exothermic for Al12P12 when compared with those of Al12N12 complexes. The interaction energies in endo-Mn@Al12P12, endo-Fe@Al12P12 and endo-Zn@Al12P12 are −7.83, −19.14 and −17.01 kcal mol−1, respectively. The exergonic energies of the interaction are mainly because of the lack of repulsive interactions from the nanocage with the encapsulated transition metals. In the absence of repulsive forces, attractive forces are only playing their role and the interaction energies are negative. Moreover, the interaction energy of Ni with Al12P12 (endo-Ni@Al12P12) is the highest in the series, a characteristic common with Al12N12 nanocage. Despite these similarities, a marked difference is also observed. The odd–even oscillation of interaction energies of early transition metals with Al12N12 nanocages is not observed for Al12P12 nanocage. The interaction energies of endo-Sc@Al12P12, endo-Ti@Al12P12, endo-V@Al12P12 and endo-Cr@Al12P12 are −58.36, −56.00, −55.83 −53.10 kcal mol−1, respectively. For Al12P12, the highest interaction energy is observed for endo-Ni@Al12P12, followed by endo-Co@Al12P12 and endo-Cu@Al12P12.

3.4. Expansion of the Nanocage

It was of quite interest to see whether encapsulation causes expansion of the nanocage. The expansion of the nanocage is calculated by diameter of the nanocages between two face-to-face hexagons (Figure 1) before and after encapsulation. Analysis of the data reveals that encapsulation of transition metal causes expansion of the nanocages; however, an irregular trend is observed. For endo-TM@Al12N12 complexes, the maximum increase in the diameter is observed for Sc encapsulation (0.157 Å), and it decreases slightly for endo-Ti@Al12N12 (0.154 Å). The diameters of endo-Sc@Al12N12 and endo-Ti@Al12N12 are 4.597 and 4.594 Å, respectively, compared to 4.44 for the bare nanocage, and the maximum increase in the diameter of these nanocage (endo-Sc@Al12N12 and endo-Ti@Al12N12) is attributed to their larger size. The expansion of the nanocage decreases with the increase in the atomic number of the transition metal up to endo-Ni@Al12N12, except for endo-Mn@Al12N12, where slightly pronounced increase in the expansion is observed. This irregular pronounced expansion of endo-Mn@Al12N12 is attributed to half-filled nature of Mn, which resists in coordination with other atoms. The data reveal that the Al12N12 nanocage is not affected by nickel encapsulation which indicates that the size of the nickel atom is ideal to fit in the cavity of the nanocage. After endo-Ni@Al12N12, the ring expansion increases again. The diameters of endo-Cu@Al12N12 and endo-Zn@Al12N12 are 4.480 and 4.503, respectively, compared to 4.44 for the bare nanocage.
Quite similar to Al12N12 nanocage, an irregular trend is also seen for endo-TM@Al12P12. However, there are certain common features. For example, the maximum expansion (0.082) is seen for endo-Sc@Al12P12. The diameter of endo-Sc@Al12P12 is 5.532 Å compared to 5.45 Å for the bare Al12P12 nanocage. Another common characteristic between the two nanocages is the irregular increase in the diameter for endo-Zn@Al12N12 and endo-Mn@Al12N12, which is again attributed to lack of any tendency for coordination due to completely filled and half-filled d orbitals. Surprisingly, the Al12P12 nanocage shrinks slightly to accommodate copper atoms. The diameter of endo-Cu@Al12P12 is 5.44 compared to 5.45 for the bare nanocage. Similar shrinkage of the nanocage is also seen for endo-Co@Al12P12, endo-V@Al12P12 and endo-Ni@Al12P12. The maximum shrinkage in the diameter of the nanocage is seen for endo-Ni@Al12P12. The data reveal that the Al12P12 nanocage is ideal for copper atom to fit in.

3.5. Charge Analysis and Dipole Moment

Next, attention was paid to a very important question: what is the charge of encapsulated metal in these nanocages? Since transition metal offers more diversity in their bonding, their behavior towards charges is worth studying. It was initially thought that placing a transition metal inside inorganic fullerenes will lead to the formation of metalides. Charge analysis reveals that the charge of encapsulated transition metal in the nanocage is strongly dependent on the nature of transition metal and the nanocage. For endo-TM@Al12N12, the general trend is a decrease in positive charge (or increase in the negative charge) while moving from Sc to Zn. (A few exceptions do exist.) The early transition metals are electron deficient; they tend to receive positive charge when encapsulated in the nanocage. The late transition metals are relatively less electron deficient, and they gain negative charge upon encapsulation. The positive charge of early transition metals decreases as the atomic number increases. The NBO charges on Sc, Ti, V and Cr in their respective complexes of Al12N12 endo-Sc@Al12N12, endo-Ti@Al12N12, endo-V@Al12N12 and endo-Cr@Al12N12 are 0.53, 0.41, 0.37 and −0.04, respectively. It is interesting to note that the positive charges on these metal atoms are not as pronounced as expected for early transition metals. Cr in endo-Cr@Al12N12 is negatively charged, although at a very small number. These findings are quite interesting and justify the proposed notion. The negative (or lower positive) charges on the transition metals in these complexes are mainly due to the presence of a large numbers of aluminum atoms in the vicinity, which take most of the positive charge. Moreover, the aluminum atoms in these inorganic fullerenes are bound to interact with transition metals at distances much shorter than their van der Waals radii, as compared to electronegative nitrogen atoms. Therefore, the effect of the aluminum atom on the charges of transition metals are more pronounced than the nitrogen atoms. The anionic character of the transition metal is also noticed for the complexes of the late transition metals. The charges on Co, Ni, Cu and Zn in endo-Co@Al12N12, endo-Ni@Al12N12, endo-Cu@Al12N12 and endo-Zn@Al12N12 are −0.11, −0.27, −0.27 and −0.17, respectively. The charges increase from Co to Cu, and then a slight decrease is observed for Zn. The decrease in anionic character for Zn may be attributed to the weak bonding of Zn with the nanocage. The strongest anionic character is seen for Ni and Cu (0.27 charge each) in their respective complexes. Although the charges on the transition metals are not very strong negatives, at least they provide a new way to generate negative charges on transition metals by encapsulating in inorganic fullerenes.
Interestingly, the negative character of the transition metals is pronounced in Al12P12 nanocages. All encapsulated transition metals except Sc, Ti and V show a negative charge. Sc, Ti and V are positively charged in their complexes; however, their charges are lesser than their corresponding charges in Al12N12 nanocages. For example, the charge of Sc in endo-Sc@Al12P12 is 0.36 compared to 0.53 for endo-Sc@Al12N12. For endo-Cr@Al12P12, endo-Fe@Al12P12 and endo-Co@Al12P12, the negative charges are quite low (in the range of −0.02 to −0.06. The highest negative charge of −0.41 is calculated for endo-Zn@Al12P12. For Al nanocages, the phosphides are better for metalide character, whereas among B nanocages, nitrides are better. It can be rationalized on soft and hard acid base (SHAB) concept.
Another characteristic closely associated with the charges is the dipole moment. The dipole moments of endo-TM@X12Y12 are also analyzed. The dipole moments do not show any regular trend, quite similar to the charges. However, one thing is quite obvious, which is that the complexes where transition metal is present at the center of the nanocage (such as endo-Cr@Al12N12, endo-Mn@Al12N12 and endo-Zn@Al12N12) do not show any appreciable dipole moments. Dipole moments of these complexes are in the range of 0.05–0.07 D. This is quite logical because the dipole moment depends on the charges as well as separation between them. The transition metal atoms in these complexes bear low charges. Moreover, individual dipole moment vectors are cancelled out due to the centrosymmetric nature of the complex. The highest dipole moment (6.11 D) is observed for endo-Sc@Al12N12 nanocages followed by 4.64 D for endo-Co@Al12N12. The dipole moments of endo-Ti@Al12N12 and endo-V@Al12N12 (3.14–3.15 D) are also significantly high. From the results, one more obvious fact is that the magnitude of the dipole moment is not related to the charge only. For example, the Ni and Cu both have similar charges (−0.27) in endo-Ni@Al12N12 and endo-Cu@Al12N12, respectively, but the dipole moments are remarkably different. The dipole moments of endo-Ni@Al12N12 and endo-Cu@Al12N12 are 0.24 and 4.64 D, respectively. This is mainly due to the separation of charges in the nanocage between Al and N atoms. For endo-Cu@Al12N12, the charge on the nitrogen and aluminum atoms of the nanocages is higher compared to those in endo-Ni@Al12N12 (see Table 3). The symmetry-group-related discussion for both optimized nanocage Al12N12 and Al12P12 is delineated in the spin-state analysis (Section 3.1).
As shown in Table 3 and Figure 3, the dipole moments of endo-TM@Al12P12 follow a trend similar to those of endo-TM@Al12N12, where early transition-metals-based complexes show higher dipole moments, whereas lower values are seen for complexes of middle transition metals endo-TM@Al12P12 (M = Cr, Mn and Fe). The dipole moments of these complexes (endo-TM@Al12P12 (M = Cr, Mn and Fe)) and endo-Zn@Al12P12 are low due to the same reason as seen above for endo-TM@Al12N12; the transition metals are present exactly at the center of the fullerenes, and individual dipole moments are cancelled out. In this series, the highest dipole moment is seen for endo-V@Al12P12.

3.6. Frontier Molecular Orbital Analysis

Transition metals offer much diversity for their charges and provide small incremental steps for changes in properties with changes in atomic number and radius. Since some of the metals show metalide character while others do not, their frontier orbital analysis will deliver useful information. The highest occupied molecular orbitals of all transition metal doped nanocages are analyzed. For endo-TM@Al12N12, the early transition metals (Sc-V) show positive charge in Al12N12 nanocage, whereas late transition metals (Ni, Cu, Zn) show metalide character.
The HOMOs of these complexes show quite interesting chemistry lying underneath. For endo-Sc@Al12N12, endo-Ti@Al12N12 and endo-V@Al12N12, the HOMOs are localized mainly on the transition metal (Figure 4). The orbitals plots show that the highest occupied molecular orbitals in these complexes are d orbitals of transition metals. The presence of HOMO on d orbitals indicates that these complexes are excess electron systems. These d orbital electrons are pushed out under the influence of lone pairs of nitrogen in Al12N12 nanocage. These excess electrons are well known to impart large nonlinear optical responses [67,68,69,70,71,72]. We and others have shown previously that such excess electrons on alkali and transition metals can lead to significantly large nonlinear optical response in these systems, particularly for Al12N12 cluster.
The situation is quite interesting for endo-Cr@Al12N12, where HOMO has density in empty space in close proximity of the transition metal, which indicates that this complex is an electride. The NBO charge analysis reveals that transition metal in this case (Cr) has a slight negative charge (metalide), and one should expect HOMO density on transition metal for metalides. Since the HOMO density of this electride is present in empty space near the transition metal, it is easy to understand that the small negative charge in the transition metal is due to the electride characteristic where some part of the density is on the metal. Therefore, the endo-Cr@Al12N12 complex is at the borderline between the metalide and the electride; however, the latter effect is dominating. A similar borderline (between metalide and electride) situation is also seen for endo-Co@Al12N12, where the HOMO indicates electride character; however, the NBO charges show some metalide characteristics. A difference between endo-Cr@Al12N12 and endo-Co@Al12N12 is the positioning of the density in HOMO. In the former, the density is slightly away from the metal center, but in the latter case, it is much closer to the metal center.
Endo-Ni@Al12N12 and endo-Cu@Al12N12 complexes are true metalides because in these cases, the HOMOs are centered on the transition metals, which when combined with NBO charge analysis (negative charge on metal) confirms the metalide character. The HOMO is a dxy type orbital in endo-Ni@Al12N12, whereas it is a dz2 orbital in endo-Cu@Al12N12. Endo-Zn@Al12N12 complex is also a metalide, but the HOMO in this complex is not present on a d orbital; rather, it is an s orbital on Zn. In summary, endo-TM@Al12N12 complexes offer quite a diversity where the first three members have excessive electron phenomenon (without any metalide or electride characteristic), while the next three members are at the borderline between metalide and electride. The last four members of the series are true metalides.
Next, the frontier orbitals of endo-TM@Al12P12 complexes are analyzed. First, three members of the series (endo-Sc@Al12P12, endo-Ti@Al12P12 and endo-V@Al12P12) show behavior very similar to those of their corresponding endo-TM@Al12N12 complexes (Figure 5); these are simply excess electron systems which are characterized by the presence of HOMO on transition metal (positive charge on metal from NBO analysis). Although the HOMOs are present on transition metals in these complexes, the d orbital involved in HOMO (dz2) is different than those of the corresponding endo-TM@Al12N12 complexes where dxy type orbital is involved in HOMOs. Although the first three members in both series show excess electron nature, the complexes from middle transition metals show quite opposite behavior. Endo-TM@Al12P12 (TM = Cr-Zn) are all metalides except endo-Fe@Al12P12, which shows electride characteristics. An interesting feature of all these endo-TM@Al12P12 (TM = Cr-Zn) complexes is the involvement of dz2 orbital in HOMOs, except endo-Zn@Al12P12, where s orbital is involved.
Next, the energies of the frontier molecular orbitals are analyzed (Table 4). For endo-TM@Al12N12, the HOMO-LUMO gaps are reduced significantly by encapsulation of transition metals when compared with the bare Al12N12 nanocage (3.93 eV). The HOMO-LUMO gap (EH-L) for endo-Sc@Al12N12 is 1.67 eV, which shows a reduction of 57.6% compared to that of the undoped system. An odd–even effect on the EH-L is seen in this series. However, the lowest EH-L is seen for endo-Fe@Al12N12, where the EH-L is 1.55 eV and the percent change in EH-L is about 60.5. The next lowest EH-L is seen for endo-V@Al12N12. The highest EH-L in this series is seen for endo-Zn@Al12N12. Doping of transition metals in Al12P12 nanocages also reduces the EH-L; however, the effect is not as much pronounced as for the endo-TM@Al12N12. For endo-TM@Al12P12, the lowest EH-L is seen for endo-Sc@Al12P12, where the EH-L is reduced to 1.53 eV compared to 3.39 eV for the bare nanocage. The next pronounced change in EH-L is seen for endo-Co@Al12P12, where the EH-L is 1.72 eV (49.4% decrease). For endo-TM@Al12P12, an odd–even effect is seen for EH-L, where the odd electron-containing species havelower EH-L, whereas the even electron-containing species have higher EH-L, except endo-Mn@Al12P12. The odd behavior of endo-Mn@Al12P12 may be attributed to the high stability of the Mn because of half-filled d orbital. A similar but more pronounced effect is seen for endo-Zn@Al12P12, where EH-L is, in fact, increased to 3.62 eV (quite contrary to the rest of the series). The EH-L analysis reveals that the gaps are reduced significantly by doping. This reduction in the energy gap is quite crucial for their potential application in optoelectronic devices because the large gaps of the bare nanocages hinder their applications. It is also worth mentioning that the EH-L are reduced, though they are not reduced too much to the extent that these systems are kinetically unstable. The EH-L are in reasonable range for enough electronic stability. For most of the optical applications’ higher transmittance, the visible region is of great interest. The NLO materials having an energy gap in the visible region finds application in material devices such as sensors, light emitters, photonics, optical communication and many other newly emerging applications [73,74,75]. Overall, the higher transmission in the visible range enables these NLO materials to be promising candidates for optoelectronic applications.

3.7. Nonlinear Optical Response

The literature reveals that alkalides are the materials with a high nonlinear optical response. It is worth investigating whether these “metalides” or other transition-metal-doped endohedral fullerenes have potential for nonlinear optical materials or not. The polarizability and first hyperpolarizability (nonlinear optical response) of the endohedral “metalides” endo-TM@X12Y12 is discussed in Table 5 and compared with those of the endohedral alkalides and others from the literature.
The polarizability of bare nanocages is 285 and 598 a.u. for Al12N12 and Al12P12, respectively. For endo-TM@Al12N12, an interesting odd–even oscillation in polarizability values is observed. For example, the polarizability of endo-Sc@Al12N12 is 486 au, and it decreases to 341 au for endo-Ti@Al12N12, but again increases to 522 for endo-V@Al12N12. The highest polarizability in this series is seen for endo-Mn@Al12N12, whereas the lowest is for endo-Ni@Al12N12. The odd–even oscillation seen for endo-TM@Al12N12 is not observed for other complexes.
The bare nanocages (Al12N12, Al12P12) are centrosymmetric and therefore have zero hyperpolarizability. Introduction of a transition metal changes the hyperpolarizability significantly. From the results, it appears that the hyperpolarizability of endo-TM@X12Y12 complexes is independent of the nature of the species whether the complex is a transition metalide or doped transition metal with a positive charge. For example, for endo-TM@Al12N12, the hyperpolarizability follows an odd–even effect, quite similar to polarizability. In this odd–even effect, the highest hyperpolarizability is seen for endo-V@Al12N12 (2.80 × 105 au). The next highest hyperpolarizability is calculated for endo-Sc@Al12N12. The highest hyperpolarizability for any transition metalide is calculated for endo-Co@Al12N12 (1.17 × 104). Many other complexes with metalide character such as endo-Cr@Al12N12, endo-Ni@Al12N12, endo-Cu@Al12N12 and endo-Zn@Al12N12 have relatively low hyperpolarizability values. Interestingly, all these species except endo-Cu@Al12N12 have metals with an even number of electrons. The hyperpolarizability of endo-Cu@Al12N12 is 2342 au. From the results for endo-Cr@Al12N12, it is clear that the hyperpolarizability of odd electron systems is higher, an effect similar to alkali metal doping. Alkali metals have an odd number of electrons, and they impart large nonlinear optical response to inorganic fullerenes upon doping. It is also worth mentioning that the systems where transition metal resides in the center of the nanocage have relatively lower hyperpolarizability (as expected from the trend). For example, the hyperpolarizability of endo-Mn@Al12N12 is merely 1.1 × 103 au, whereas it is 2.80 × 105 au for endo-V@Al12N12.
The odd–even oscillating trend of hyperpolarizability is also seen for endo-TM@Al12P12 to some extent. Moreover, the hyperpolarizability values of complexes in this series are lower than those of endo-TM@Al12N12. For example, the highest values of hyperpolarizability in this series are calculated for endo-Sc@Al12P12 (3954 au), which is about an order of magnitude smaller than that of endo-Sc@Al12N12 (1.76 × 104 au). The hyperpolarizability drops to 2345 au for endo-Ti@Al12P12, but again increases to 2773 au for endo-V@Al12P12 followed by another decrease to 403 au for endo-Cr@Al12P12. In this series, the lowest value is calculated for endo-Zn@Al12P12, followed by endo-Mn@Al12P12. Quite similar to AlN nanocage, the only metalide with reasonably high hyperpolarizability is endo-Co@Al12P12 (probably mainly due to odd–even effect). Quite similar to AlN nanocages, the early transition-metal-doped AlP nanocages have relatively high hyperpolarizability. Both of these complexes have transition metal either half-filled or completely filled. Moreover, the transition metals are present exactly at the center of the nanocage. Both of these factors are probably responsible for low hyperpolarizabilities of these complexes.
For comparison, the polarizability and hyperpolarizability of endo-TM@X12Y12 complexes are compared with some other reported metalides. Table 6 reveals that the first static hyperpolarizabilities of the studied complexes designed in this investigation are very much comparable or greater relative to the complexes reported in the literature, especially in case of endo-V@Al12N12.

4. Conclusions

Electronic and nonlinear optical properties of endohedral metallofullerenes are studied through density functional theory calculations. The endohedral metallofullerenes contain transition metal encapsulated in inorganic fullerenes X12Y12 (X = Al & Y = N, P). The endohedral metallofullerenes (endo-TM@X12Y12) possess quite interesting geometric and electronic properties, which are the function of the nature of atom and the size of fullerene. Early and late transition metals except endo-Zn@X12Y12 tend to sit at the side of the cavity, whereas the middle transition metal atoms reside at the center of the fullerenes. NBO charge and frontier molecular orbital analyses reveal that the endo-TM@Al12N12 are true metalides when transition metals are Ni, Cu and Zn. Additionally, endo-Cr@Al12N12 and endo-Co@Al12N12 are at the borderline between metalide and electride with predominant electride characteristics. The other members of the series are excess electron systems, which offer interesting electronic and nonlinear optical properties. The diversity of nature possessed by endo-TM@Al12N12 is not prevalent for other fullerenes. Additionally, endo-TM@Al12P12 are true metalides when transition metals are (Cr-Zn). Among all complexes studied, the highest metalide character is seen for endo-Zn@Al12P12, where the charge on the metal reaches up to −0.41. HOMO-LUMO gaps (EH-L) are reduced significantly for these endohedral metallofullerenes. Many complexes show odd–even oscillating behavior for EH-L and dipole moments. Odd electron species contain large dipole moments and small EH-L, whereas even electron systems have the opposite behavior. Despite a decrease in EH-L, these systems show high kinetic and thermodynamic stabilities. The encapsulation of transition metals is a highly exergonic process. For all series, nickel-doped complexes are the most stable thermodynamically. These endo-TM@X12Y12 possess remarkable nonlinear optical response where the first hyperpolarizability reaches up to 2.79 × 105 au for endo-V@Al12N12. With few exceptions, the polarizabilities of the complexes are inversely related to the HOMO-LUMO gaps (EH-L). This study helps in the comparative analysis of nonlinear optical response of electrides, metalides and other excess electron systems. In general, the nonlinear optical response of electrides is higher than metalides but lower than those of simple excess electron compounds. In summary, the interesting electronic and large nonlinear optical response of the studied metalides shows that these materials are promising candidates for optoelectronic applications.

Author Contributions

Data curation, N.S.S. and K.A.; Formal analysis, N.S.S., M.A.A., Y.S.S.A.-F. and I.B.; Investigation, M.A.A., Y.S.S.A.-F., I.B. and N.S.S.; Methodology, N.S.S., M.A.A., Y.S.S.A.-F. and I.B.; Project administration, N.S.S. and K.A.; Resources, M.A.A., N.S.S., R.L. and K.A.; Software, N.S.S. and K.A.; Validation, N.S.S., M.A.A. and K.A.; Visualization, N.S.S., I.B. and K.A.; Writing—original draft, N.S.S., M.A.A. and K.A.; Writing—review and editing, N.S.S. and K.A. All authors have read and agreed to the published version of the manuscript.

Funding

Deputyship for Research & Innovation, Ministry of Education, Saudi Arabia (INST208).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data are provided in the manuscript.

Acknowledgments

The authors extend their appreciation to the Deputyship for Research and Innovation, Ministry of Education in Saudi Arabia for funding this research work (Project number INST208).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dye, J.L.; DeBaker, M.G. Physical and Chemical Properties of Alkalides And Electrides. Ann. Rev. Phys. Chem. 1987, 38, 271–301. [Google Scholar] [CrossRef]
  2. Sun, W.-M.; Li, Y.; Li, X.-H.; Wu, D.; He, H.-M.; Li, C.-Y.; Chen, J.-H.; Li, Z.-R. Stability and Nonlinear Optical Response of Alkalides that Contain a Completely Encapsulated Superalkali Cluster. ChemPhysChem 2016, 17, 2672–2678. [Google Scholar] [CrossRef] [PubMed]
  3. Redko, M.Y.; Huang, R.H.; Jackson, J.E.; Harrison, J.F.; Dye, J.L. Barium Azacrypt and Sodide, the First Alkalide with an Alkaline Earth Cation, Also Contains a Novel Dimer, (Na2)2−. J. Am. Chem. Soc. 2003, 125, 2259–2263. [Google Scholar] [CrossRef] [PubMed]
  4. Sun, W.-M.; Wu, D.; Li, Y.; Li, Z.-R. Substituent Effects on the Structural Features and Nonlinear Optical Properties of the Organic Alkalide Li+(calix[4]pyrrole)Li. ChemPhysChem 2013, 14, 408–416. [Google Scholar] [CrossRef] [PubMed]
  5. Sun, W.-M.; Fan, L.-T.; Li, Y.; Liu, J.-Y.; Wu, D.; Li, Z.-R. On the Potential Application of Superalkali Clusters in Designing Novel Alkalides with Large Nonlinear Optical Properties. Inorg. Chem. 2014, 53, 6170–6178. [Google Scholar] [CrossRef] [PubMed]
  6. Jing, Y.-Q.; Li, Z.-R.; Wu, D.; Li, Y.; Wang, B.-Q.; Gu, F.L.; Aoki, Y. Effect of the Complexant Shape on the Large First Hyperpolarizability of Alkalides Li+(NH3)4M. ChemPhysChem 2006, 7, 1759–1763. [Google Scholar] [CrossRef]
  7. Chen, W.; Li, Z.-R.; Wu, D.; Li, Y.; Sun, C.-C.; Gu, F.L.; Aoki, Y. Nonlinear Optical Properties of Alkalides Li+(calix[4]pyrrole)M (M = Li, Na, and K): Alkali Anion Atomic Number Dependence. J. Am. Chem. Soc. 2006, 128, 1072–1073. [Google Scholar] [CrossRef]
  8. Sun, W.-M.; Li, X.-H.; Wu, J.; Lan, J.-M.; Li, C.-Y.; Wu, D.; Li, Y.; Li, Z.-R. Can Coinage Metal Atoms Be Capable of Serving as an Excess Electron Source of Alkalides with Considerable Nonlinear Optical Responses? Inorg. Chem. 2017, 56, 4594–4600. [Google Scholar] [CrossRef]
  9. Wang, F.-F.; Li, Z.-R.; Wu, D.; Wang, B.-Q.; Li, Y.; Li, Z.-J.; Chen, W.; Yu, G.-T.; Gu, F.L.; Aoki, Y. Structures and Considerable Static First Hyperpolarizabilities: New Organic Alkalides (M+@n6adz)M′(M,M′ = Li, Na, K; n = 2, 3) with Cation Inside and Anion Outside of the Cage Complexants. J. Phys. Chem. B 2008, 112, 1090–1094. [Google Scholar] [CrossRef]
  10. Sun, W.-M.; Ni, B.-L.; Wu, D.; Lan, J.-M.; Li, C.-Y.; Li, Y.; Li, Z.-R. Designing Alkalides with Considerable Nonlinear Optical Responses and High Stability Based on the Facially Polarized Janus all-cis-1,2,3,4,5,6-Hexafluorocyclohexane. Organometallics 2017, 36, 3352–3359. [Google Scholar] [CrossRef]
  11. Matalon, S.; Golden, S.; Ottolenghi, M. Nature of the visible absorption bands in metal-amine solutions. J. Phys. Chem. 1969, 73, 3098–3101. [Google Scholar] [CrossRef]
  12. DeBacker, M.G.; Dye, J.L. Metal-ethylenediamine solutions. Extinction coefficients and equilibriums. J. Phys. Chem. 1971, 75, 3092–3096. [Google Scholar] [CrossRef]
  13. Lok, M.T.; Tehan, F.J.; Dye, J.L. Spectra of Na-, K-, and e-solv in amines and ethers. J. Phys. Chem. 1972, 76, 2975–2981. [Google Scholar] [CrossRef]
  14. Ceraso, J.M.; Dye, J.L. 23 Na NMR spectrum of the sodium anion. J. Chem. Phys. 1974, 61, 1585–1587. [Google Scholar] [CrossRef]
  15. Dye, J.L.; Ceraso, J.M.; Lok, M.; Barnett, B.L.; Tehan, F.J. Crystalline salt of the sodium anion (Na-). J. Am. Chem. Soc. 1974, 96, 608–609. [Google Scholar] [CrossRef]
  16. Tehan, F.J.; Barnett, B.L.; Dye, J.L. Alkali anions. Preparation and crystal structure of a compound which contains the cryptated sodium cation and the sodium anion. J. Am. Chem. Soc. 1974, 96, 7203–7208. [Google Scholar] [CrossRef]
  17. Kim, J.; Ichimura, A.S.; Huang, R.H.; Redko, M.; Phillips, R.C.; Jackson, J.E.; Dye, J.L. Crystalline Salts of Na-and K-(Alkalides) that Are Stable at Room Temperature. J. Am. Chem. Soc. 1999, 121, 10666–10667. [Google Scholar] [CrossRef]
  18. Chen, W.; Li, Z.-R.; Wu, D.; Li, Y.; Li, R.-Y.; Sun, C.-C. Inverse Sodium Hydride: Density Functional Theory Study of the Large Nonlinear Optical Properties. J. Phys. Chem. A 2005, 109, 2920–2924. [Google Scholar] [CrossRef]
  19. Redko, M.; Vlassa, M.; Jackson, J.E.; Misiolek, A.W.; Huang, R.H.; Dye, J.L. Inverse Sodium Hydride”: A Crystalline Salt that Contains H+ and Na. J. Am. Chem. Soc. 2002, 124, 5928–5929. [Google Scholar] [CrossRef]
  20. Mutailipu, M.; Poeppelmeier, K.R.; Pan, S. Borates: A Rich Source for Optical Materials. Chem. Rev. 2021, 121, 1130–1202. [Google Scholar] [CrossRef]
  21. Mutailipu, M.; Yang, Z.; Pan, S. Toward the Enhancement of Critical Performance for Deep-Ultraviolet Frequency-Doubling Crystals Utilizing Covalent Tetrahedra. Acc. Mater. Res. 2021, 2, 282–291. [Google Scholar] [CrossRef]
  22. Huang, J.; Shu, S.; Cai, G.-M. Review of Heteroleptic Tetrahedra as Birefringent or Nonlinear Optical Motifs. Cryst. Growth Des. 2022, 22, 1500–1514. [Google Scholar] [CrossRef]
  23. Karpov, A.; Nuss, J.; Wedig, U.; Jansen, M. Cs2Pt: A Platinide(-II) Exhibiting Complete Charge Separation. Angew. Chemie Int. Ed. 2003, 42, 4818–4821. [Google Scholar] [CrossRef]
  24. Sommer, A. Alloys of Gold with Alkali Metals. Nature 1943, 152, 215. [Google Scholar] [CrossRef]
  25. Wood, V.E.; Reitz, J.R. Electronic band structure of cesium gold. J. Phys. Chem. Solids. 1962, 23, 229–235. [Google Scholar] [CrossRef]
  26. Liu, T.L. Optical properties of RbAu and CsAu. Phys. Rev. B 1975, 12, 3008–3012. [Google Scholar] [CrossRef]
  27. Knecht, J.; Fischer, R.; Overhof, H.; Hensel, F. ESCA study of compounds of gold in the oxidation state −1. J. Chem. Soc. Chem. Commun. 1978, 21, 905–906. [Google Scholar] [CrossRef]
  28. Köhler, J.; Chang, J.-H. [PtIn6]10+ Octahedra in PtIn7F13: The First Compound of a New Class of Metal-Cluster Fluorides. Angew. Chemie Int. Ed. 2000, 39, 1998–2000. [Google Scholar] [CrossRef]
  29. Köhler, J.; Friedrich, C.; Lee, H.A.; Whangbo, M.-H. IrIn7GeO8 = [IrIn6](GeO4)(InO4) und Verbindungen der Mischkristallreihe [IrIn6](Ge1+xIn1−4x/3O8) (0 ≤ x ≤ 0.75)–erste Oxide mit [IrIn6]-Oktaedern. Z. Anorg. Allg. Chem. 2007, 633, 1464–1471. [Google Scholar] [CrossRef]
  30. Li, X.-H.; Zhang, X.-L.; Chen, Q.-H.; Zhang, L.; Chen, J.-H.; Wu, D.; Sun, W.-M.; Li, Z.-R. Coinage metalides: A new class of excess electron compounds with high stability and large nonlinear optical responses. Phys. Chem. Chem. Phys. 2020, 22, 8476–8484. [Google Scholar] [CrossRef]
  31. Malloci, G.; Cappellini, G.; Mulas, G.; Satta, G. Quasiparticle effects and optical absorption in small fullerene like GaP clusters. Phys. Rev. B 2004, 70, 205429. [Google Scholar] [CrossRef]
  32. Tozzini, V.; Buda, F.; Fasolino, A. Fullerene-like III−V clusters: A density functional theory prediction. J. Phys. Chem. B 2001, 105, 12477–12480. [Google Scholar] [CrossRef]
  33. Wu, H.-S.; Zhang, F.-Q.; Xu, X.-H.; Zhang, C.-J.; Jiao, H. Geometric and Energetic Aspects of Aluminum Nitride Cages. J. Phys. Chem. A 2003, 107, 204–209. [Google Scholar] [CrossRef]
  34. Seifert, G.; Fowler, P.W.; Mitchell, D.; Porezag, D.; Frauenheim, T. Boron-nitrogen analogues of the fullerenes: Electronic and structural properties. Chem. Phys. Lett. 1997, 268, 352–358. [Google Scholar] [CrossRef]
  35. Shokuhi Rad, A.; Ayub, K. Adsorption of pyrrole on Al12N12, Al12P12, B12N12, and B12P12 fullerene-like nano-cages; a first principles study. Vacuum 2016, 131, 135–141. [Google Scholar] [CrossRef]
  36. Shokuhi Rad, A.; Ayub, K. A comparative density functional theory study of guanine chemisorption on Al12N12, Al12P12, B12N12, and B12P12 nano-cages. J. Alloys Compd. 2016, 672, 161–169. [Google Scholar] [CrossRef]
  37. Rad, A.S.; Ayub, K. How can nickel decoration affect H2 adsorption on B12P12 nano-heterostructures? J. Mol. Liq. 2018, 255, 168–175. [Google Scholar] [CrossRef]
  38. Ayub, K. Transportation of hydrogen atom and molecule through X12Y12nano-cages. Int. J. Hydrog. Energy 2017, 42, 11439–11451. [Google Scholar] [CrossRef]
  39. Ayub, K. Binding affinity and permeation of X12Y12 nanoclusters for helium and neon. J. Mol. Liq. 2017, 244, 124–134. [Google Scholar] [CrossRef]
  40. Iqbal, J.; Ludwig, R.; Ayub, K. Phosphides or nitrides for better NLO properties? A detailed comparative study of alkali metal doped nano-cages. Mater. Res. Bull. 2017, 92, 113–122. [Google Scholar] [CrossRef]
  41. Iqbal, J.; Ayub, K. Enhanced electronic and non-linear optical properties of alkali metal (Li, Na, K) doped boron nitride nano-cages. J. Alloys Compd. 2016, 687, 976–983. [Google Scholar] [CrossRef]
  42. Iqbal, J.; Ayub, K. Theoretical study of the non linear optical properties of alkali metal (Li, Na, K) doped aluminum nitride nanocages. RSC Adv. 2016, 6, 94228–94235. [Google Scholar] [CrossRef]
  43. Munsif, S.; Khan, S.; Ali, A.; Gilani, M.A.; Iqbal, J.; Ludwig, R.; Ayub, K. Remarkable nonlinear optical response of alkali metal doped aluminum phosphide and boron phosphide nanoclusters. J. Mol. Liq. 2018, 271, 51–64. [Google Scholar] [CrossRef]
  44. Ayub, K. Are phosphide nano-cages better than nitride nano-cages? A kinetic, thermodynamic and non-linear optical properties study of alkali metal encapsulated X12Y12 nano-cages. J. Mater. Chem. C 2016, 4, 10919–10934. [Google Scholar] [CrossRef]
  45. Ahsan, A.; Khan, S.; Gilani, M.A.; Ayub, K. Endohedral metallofullerene electrides of Ca12O12with remarkable nonlinear optical response. RSC Adv. 2021, 11, 1569–1580. [Google Scholar] [CrossRef] [PubMed]
  46. Sajid, H.; Ullah, F.; Khan, S.; Ayub, K.; Arshad, M.; Mahmood, T. Remarkable static and dynamic NLO response of alkali and superalkali doped macrocyclic [hexa-]thiophene complexes; a DFT approach. RSC Adv. 2021, 11, 4118–4128. [Google Scholar] [CrossRef] [PubMed]
  47. Kosar, N.; Tahir, H.; Ayub, K.; Mahmood, T. DFT studies of single and multiple alkali metals doped C24 fullerene for electronics and nonlinear optical applications. J. Mol. Graph. Model. 2021, 105, 107867. [Google Scholar] [CrossRef] [PubMed]
  48. Sohail, M.; Khaliq, F.; Mahmood, T.; Ayub, K.; Tabassum, S.; Gilani, M.A. Influence of bi-alkali metals doping over Al12N12 nanocage on stability and optoelectronic properties: A DFT investigation. Radiat. Phys. Chem. 2021, 184, 109457. [Google Scholar] [CrossRef]
  49. Baloach, R.; Ayub, K.; Mahmood, T.; Asif, A.; Tabassum, S.; Gilani, M.A. A New Strategy of bi-Alkali Metal Doping to Design Boron Phosphide Nanocages of High Nonlinear Optical Response with Better Thermodynamic Stability. J. Inorg. Organomet. Polym. Mater. 2021, 31, 3062–3076. [Google Scholar] [CrossRef]
  50. Perez-Jimenez, A.J.; Sancho-García, J.C. Using circumacenes to improve organic electronics and molecular electronics: Design clues. Nanotechnology 2009, 20, 475201. [Google Scholar] [CrossRef]
  51. Mocci, P.; Cardia, R.; Cappellini, G. Si-atoms substitutions effects on the electronic and optical properties of coronene and ovalene. New J. Phys. 2018, 20, 113008. [Google Scholar] [CrossRef]
  52. Mocci, P.; Cardia, R.; Cappellini, G. Inclusions of Si-atoms in Graphene nanostructures: A computational study on the ground-state electronic properties of Coronene and Ovalene. J. Phys. Conf. Ser. 2018, 956, 012020. [Google Scholar] [CrossRef]
  53. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. Gaussian, Inc.: Wallingford, CT, USA, 2016.
  54. Dennington, R.; Keith, T.A.; Millam, J.M. GaussView, version 5.0; Gaussian, Inc.: Wallingford, CT, USA, 2009. [Google Scholar]
  55. Becke, A.D. Density-functional thermochemistry. III. The role of exact exchange. J. Chem. Phys. 1993, 98, 5648–5652. [Google Scholar] [CrossRef]
  56. Lee, C.; Yang, W.; Parr, R.G. Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. Phys. Rev. B 1988, 37, 785–789. [Google Scholar] [CrossRef] [PubMed]
  57. Hussain, R.; Hussain, A.I.; Chatha, S.A.S.; Mansha, A.; Ayub, K. Density functional theory study of geometric and electronic properties of full range of bimetallic AgnYm (n + m = 10) clusters. J. Alloys Compd. 2017, 705, 232–246. [Google Scholar] [CrossRef]
  58. Rad, A.S.; Ayub, K. Coordination of nickel atoms with Al12X12 (X = N, P) nanocages enhances H2 adsorption: A surface study by DFT. Vacuum 2016, 133, 70–80. [Google Scholar] [CrossRef]
  59. Khaliq, F.; Afzaal, A.; Tabassum, S.; Mahmood, T.; Ayub, K.; Khan, A.L.; Yasin, M.; Gilani, M.A. Surface functionalization of Si6Li6 cluster with superalkalis to achieve high nonlinear optical response: A DFT study. Colloids Surf. A Physicochem. Eng. Asp. 2022, 653, 129985. [Google Scholar] [CrossRef]
  60. Godby, R.W.; Schlüter, M.; Sham, L. Self-energy operators and exchange-correlation potentials in semiconductors. Phys. Rev. B 1988, 37, 10159. [Google Scholar] [CrossRef]
  61. Cappellini, G.; Furthmüller, J.; Cadelano, E.; Bechstedt, F. Electronic and optical properties of cadmium fluoride: The role of many-body effects. Phys. Rev. B 2013, 87, 075203. [Google Scholar] [CrossRef]
  62. Cardia, R.; Malloci, G.; Rignanese, G.M.; Blase, X.; Molteni, E.; Cappellini, G. Electronic and optical properties of hexathiapentacene in the gas and crystal phases. Phys. Rev. B 2016, 93, 235132. [Google Scholar] [CrossRef]
  63. Stagi, L.; Chiriu, D.; Scholz, M.; Carbonaro, C.M.; Corpino, R.; Porcheddu, A.; Rajamaki, S.; Cappellini, C.; Cardia, R.; Ricci, P.C. Vibrational and optical characterization of s-triazine derivatives. Spectrochim. Acta Part A Mol. Biomol. Spectrosc. 2017, 183, 348–355. [Google Scholar] [CrossRef] [PubMed]
  64. Pinna, E.; Melis, C.; Antidormi, A.; Cardia, R.; Sechi, E.; Cappellini, G.; D’Ischia, M.; Colombo, L.; Mula, G. Deciphering molecular mechanisms of interface buildup and stability in porous Si/eumelanin hybrids. Int. J. Mol. Sci. 2017, 18, 1567. [Google Scholar] [CrossRef] [PubMed]
  65. Faber, C.; Boulanger, P.; Attaccalite, C.; Duchemin, I.; Blase, X. Excited states properties of organic molecules: From density functional theory to the GW and Bethe–Salpeter Green’s function formalisms. Philos. Trans. R. Soc. A Math. Phys. Eng. Sci. 2014, 372, 20130271. [Google Scholar] [CrossRef]
  66. Oudar, J.L.; Chemla, D.S. Hyperpolarizabilities of the nitroanilines and their relations to the excited state dipole moment. J. Chem. Phys. 1977, 66, 2664–2668. [Google Scholar] [CrossRef]
  67. Ullah, F.; Kosar, N.; Ayub, K.; Mahmood, T. Superalkalis as a source of diffuse excess electrons in newly designed inorganic electrides with remarkable nonlinear response and deep ultraviolet transparency: A DFT study. Appl. Surf. Sci. 2019, 483, 1118–1128. [Google Scholar] [CrossRef]
  68. Ahsan, A.; Ayub, K. Extremely large nonlinear optical response and excellent electronic stability of true alkaline earthides based on hexaammine complexant. J. Mol. Liq. 2020, 297, 111899. [Google Scholar] [CrossRef]
  69. Ahsan, A.; Ayub, K. Adamanzane based alkaline earthides with excellent nonlinear optical response and ultraviolet transparency. Opt. Laser Technol. 2020, 129, 106298. [Google Scholar] [CrossRef]
  70. Ahsan, A.; Sarfaraz, S.; Fayyaz, F.; Asghar, M.; Ayub, K. Enhanced non-linear optical response of calix[4]pyrrole complexant based earthides in the presence of oriented external electric field. J. Mol. Liq. 2022, 350, 118504. [Google Scholar] [CrossRef]
  71. Ahsan, A.; Sarfaraz, S.; Gilani, M.A.; Mahmood, T.; Ahmad, Z.; Ayub, K. Alkaline earthides based on 15-crown-5 ether with remarkable NLO response. Eur. Phys. J. Plus 2022, 137, 1149. [Google Scholar] [CrossRef]
  72. Ahsan, A.; Sarfaraz, S.; Fayyaz, F.; Asghar, M.; Ayub, K. Theoretical study of 36Adz based alkaline earthides M+(36Adz)M (M+ = Li & Na; M = Be, Mg & Ca) with remarkable nonlinear optical response. Mater. Sci. Semicond. Process. 2023, 153, 107119. [Google Scholar] [CrossRef]
  73. Zawadzka, A.; Płóciennik, P.; Strzelecki, J.; Sahraoui, B. Transparent amorphous zinc oxide thin films for NLO applications. Opt. Mater. 2014, 37, 327–337. [Google Scholar] [CrossRef]
  74. Komine, T.; Nakagawa, M. Fundamental analysis for visible-light communication system using LED lights. IEEE Trans. Consum. Electron. 2004, 50, 100–107. [Google Scholar] [CrossRef]
  75. Du, Y.; Ai, X.; Li, Z.; Sun, T.; Huang, Y.; Zeng, X.; Chen, X.; Rao, F. Visible-to-Ultraviolet Light Conversion: Materials and Applications. Adv. Photonics Res. 2021, 2, 2000213. [Google Scholar] [CrossRef]
  76. Rad, A.S.; Ayub, K. Nonlinear optical and electronic properties of Cr-, Ni-, and Ti-substituted C20 fullerenes: A quantum-chemical study. Mater. Res. Bull. 2018, 97, 399–404. [Google Scholar] [CrossRef]
  77. Arshad, Y.; Khan, S.; Hashmi, M.A.; Ayub, K. Transition metal doping: A new and effective approach for remarkably high nonlinear optical response in aluminum nitride nanocages. New J. Chem. 2018, 42, 6976–6989. [Google Scholar] [CrossRef]
  78. Irshad, S.; Ullah, F.; Khan, S.; Ludwig, R.; Mahmood, T.; Ayub, K. First row transition metals decorated boron phosphide nanoclusters as nonlinear optical materials with high thermodynamic stability and enhanced electronic properties; A detailed quantum chemical study. Opt. Laser Technol. 2021, 134, 106570. [Google Scholar] [CrossRef]
Figure 1. Optimized geometries and the most stable spin states of endo-TM@Al12N12 complexes.
Figure 1. Optimized geometries and the most stable spin states of endo-TM@Al12N12 complexes.
Materials 16 03447 g001
Figure 2. Optimized geometries and the most stable spin states of endo-TM@Al12P12 complexes.
Figure 2. Optimized geometries and the most stable spin states of endo-TM@Al12P12 complexes.
Materials 16 03447 g002
Figure 3. Illustration of dipole moments of endo-TM@X12Y12 fullerenes.
Figure 3. Illustration of dipole moments of endo-TM@X12Y12 fullerenes.
Materials 16 03447 g003
Figure 4. Frontier molecular orbitals of endo-TM@Al12N12 complexes plotted at an isovalue of 0.045.
Figure 4. Frontier molecular orbitals of endo-TM@Al12N12 complexes plotted at an isovalue of 0.045.
Materials 16 03447 g004
Figure 5. Frontier molecular orbitals of endo-TM@Al12P12 complexes plotted at isovalue of 0.045.
Figure 5. Frontier molecular orbitals of endo-TM@Al12P12 complexes plotted at isovalue of 0.045.
Materials 16 03447 g005
Table 1. Illustration of ordering and relative energies of different spin states for endo-TM@X12Y12 complexes. All relative energies for a metal complex are with respect to its most stable spin state at 0 kcal mol−1.
Table 1. Illustration of ordering and relative energies of different spin states for endo-TM@X12Y12 complexes. All relative energies for a metal complex are with respect to its most stable spin state at 0 kcal mol−1.
Spin StatesScTiVCrMnFeCoNiCuZn
AlNStable spinDoublet (0)
C1
Quintet (0)
C1
Quartet
(0)
C1
Septet
(0)
Th
Sextet
(0)
C1
Quintet
(0)
C1
Quartet
(0)
C1
Singlet
(0)
C1
Doublet
(0)
C1
Singlet
(0)
Th
2nd stableQuartet (2.72)
C1
Triplet (2.13)
C1
Sextet (5.22)
C1
Quintet (3.18)
Th
Quartet (3.41)
C1
Triplet (1.41)
C1
Doublet (0.92)
C1
Triplet (0.44)
C1
Quartet
(44.17)
C1
Triplet (21.09)
Cs
3rd stableSextet (57.0)
C1
Singlet (9.65)
C1
Doublet (14.46)
C1
Triplet (11.46)
C1
Octet (4.40)
C1
Septet (17.18)
C1
Sextet (23.13)
C1
Quintet (12.32)
C1
Sextet (83.66)
C1
Quintet (76.46)
C1
4th stableOctet (122.9)
C1
Septet (66.64)
C1
Octet (70.54)
C1
Singlet (44.37)
C1
Doublet (19.68)
C1
Singlet (28.09)
C1
Octet (78.48)
C1
Septet
C1
Octet (154.63)
C1
Septet
C1
AlPStable spinQuartet (0)
C1
Quintet (0)
C3
Sextet
(0)
C1
Septet
(0)
C1
Sextet
(0)
Th
Quintet
(0)
C2v
Doublet
(0)
C1
Singlet
(0)
C3
Doublet (0)
C1
Singlet (0)
Th
2nd stableDoublet (4.38)
C1
Triplet (4.82)
C1
Quartet (8.93)
C1
Quintet (11.89)
C1
Quartet (19.64)
C1
Triplet (6.74)
C1
Quartet (0.43)
C1
Triplet (12.76)
C1
Quartet (45.35)
C1
Triplet (41.78)
C1
3rd stableSextet (43.88)
C1
Singlet (9.58)
C1
Doublet (12.55)
C1
Triplet (36.26)
C1
Octet (26.89)
C1
Septet (13.70)
C1
Sextet (29.09) C1Quintet (55.88)
C1
Sextet (90.70)
C1
Quintet (89.62)
C1
4th stableOctet (92.19)
C1
Septet (47.18)
Cs
Octet (38.12)
C1
Singlet (61.57)
C1
Doublet (43.58)
C1
Singlet (41.41)
C1
Octet (78.85)
C1
Septet (112.37)
C1
Octet (144.85)
C1
Septet (132.56)
C3
Table 2. Zero point corrected and Gibbs free binding energies (kcal mol−1) of transition metals with X12Y12 nanocages.
Table 2. Zero point corrected and Gibbs free binding energies (kcal mol−1) of transition metals with X12Y12 nanocages.
Transition MetalNanocage
Al12N12Al12P12
ZPEGibbsZPEGibbs
Sc−20.72−12.56−64.47−58.36
Ti−49.68−41.72−62.63−56.00
V−16.74−8.78−62.74−55.83
Cr−20.93−13.95−58.07−53.10
Mn30.6838.83−14.42−7.83
Fe2.009.79−23.84−19.14
Co−12.10−4.45−81.03−73.76
Ni−85.57−77.84−135.88−127.54
Cu−44.66−38.17−69.07−61.24
Zn21.3528.11−24.34−17.01
Table 3. NBO Charges, diameters (in Å), change in diameters (in Å), dipole moments (in Debye), bond order, orientation of metal and symmetry of endo-TM@X12Y12.
Table 3. NBO Charges, diameters (in Å), change in diameters (in Å), dipole moments (in Debye), bond order, orientation of metal and symmetry of endo-TM@X12Y12.
NanocagePropertyTransition Metals
ScTiVCrMnFeCoNiCuZn
Al12N12Charge0.530.410.37−0.040.170.07−0.11−0.27−0.27−0.17-
Diameter4.5974.5944.5114.5174.5274.504.4924.4304.4804.5034.440
Change0.1570.1540.0710.0770.0870.060.052−0.010.040.063
Dipole Moment6.113.143.150.050.073.434.640.242.840.0370.0
Bond order4.274.393.312.132.572.602.623.352.562.48-
SymmetryC1C1C1ThC3C1C1C1C1ThTh
OrientationSideSideSideCenterCenterSideSideSideSideCenter
Al12P12Charge0.360.140.16−0.02−0.06−0.12−0.03−0.17−0.2−0.41-
Diameter5.5325.4605.4235.4905.5045.4905.4285.3975.4425.5005.45
Change0.0820.010−0.0270.0400.0540.04−0.022−0.053−0.0080.05
Dipole Moment3.254.715.310.130.020.081.791.202.500.000.0
Bond order3.222.802.762.552.913.092.913.502.873.00-
SymmetryC1C3C1C1ThC2VC1C3C1ThTh
OrientationSideSideSideCenterCenterCenterSideSideSideCenter
Table 4. Energies of HOMO (EHOMO in eV), LUMO (ELUMO in eV), HOMO-LUMO gaps (EH-L in eV), percent decrease in HOMO-LUMO gap (% ΔEH-L) and first hyperpolarizabilities (βo in au) endo-TM@X12Y12.
Table 4. Energies of HOMO (EHOMO in eV), LUMO (ELUMO in eV), HOMO-LUMO gaps (EH-L in eV), percent decrease in HOMO-LUMO gap (% ΔEH-L) and first hyperpolarizabilities (βo in au) endo-TM@X12Y12.
NanocagePropertyTransition Metals
ScTiVCrMnFeCoNiCuZnBare
Al12N12ELUMO−2.67−2.05−2.68−1.67−1.78−2.62−2.13−2.99−2.08−2.22−2.54
EHOMO−4.34−4.56−4.23−3.78−4.06−4.18−4.04−4.81−4.24−4.96−6.47
EH-L1.672.511.562.122.281.551.911.822.162.743.93
% ΔEH-L57.636.060.446.242.060.551.453.745.030.40.0
βo17,590480279,4956811131963111,6911002342130
Al12P12ELUMO−3.42−3.09−3.26−2.94−3.06−3.12−3.40−3.49−3.06−3.09−3.36
EHOMO−4.95−5.06−5.20−5.62−6.06−6.04−5.12−5.96−5.64−6.71−6.75
EH-L1.531.971.942.683.002.921.722.472.593.623.39
% ΔEH-L54.841.942.820.811.413.749.427.323.7−6.70.0
βo395423452773403372 3027020.820
Table 5. Polarizability (α in au) and first hyperpolarizability (βo in au) of endo-TM@X12Y12 complexes, calculated at LC-BLYP/6-31+G(d).
Table 5. Polarizability (α in au) and first hyperpolarizability (βo in au) of endo-TM@X12Y12 complexes, calculated at LC-BLYP/6-31+G(d).
Transition MetalEndo-TM@Al12N12Endo-TM@Al12P12
αβoαβo
Sc48617,5906643954
Ti3414806602345
V522279,4956242773
Cr473681662403
Mn55911316633.0
Fe327963164372
Co34811,6916482442
Ni294100586302
Cu3222342618702
Zn326136260.82
Table 6. Comparison of nonlinear optical response of current work with already reported values over different metalide systems.
Table 6. Comparison of nonlinear optical response of current work with already reported values over different metalide systems.
Complexαo (au)βo (au)Ref.
endo-V@Al12N125222.8 × 105This work
endo-Sc@Al12P126643.9 × 103This work
Ti-C19-2.5 × 103[76]
Cu@r6-Al12N12-1.8 × 104[77]
Ni@Mg12O122611.2 × 103[76]
Sc@r4-B12P126124.4 × 104[78]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Alkhalifah, M.A.; Sheikh, N.S.; Al-Faiyz, Y.S.S.; Bayach, I.; Ludwig, R.; Ayub, K. Rational Design, Stabilities and Nonlinear Optical Properties of Non-Conventional Transition Metalides; New Entry into Nonlinear Optical Materials. Materials 2023, 16, 3447. https://doi.org/10.3390/ma16093447

AMA Style

Alkhalifah MA, Sheikh NS, Al-Faiyz YSS, Bayach I, Ludwig R, Ayub K. Rational Design, Stabilities and Nonlinear Optical Properties of Non-Conventional Transition Metalides; New Entry into Nonlinear Optical Materials. Materials. 2023; 16(9):3447. https://doi.org/10.3390/ma16093447

Chicago/Turabian Style

Alkhalifah, Mohammed A., Nadeem S. Sheikh, Yasair S. S. Al-Faiyz, Imene Bayach, Ralf Ludwig, and Khurshid Ayub. 2023. "Rational Design, Stabilities and Nonlinear Optical Properties of Non-Conventional Transition Metalides; New Entry into Nonlinear Optical Materials" Materials 16, no. 9: 3447. https://doi.org/10.3390/ma16093447

APA Style

Alkhalifah, M. A., Sheikh, N. S., Al-Faiyz, Y. S. S., Bayach, I., Ludwig, R., & Ayub, K. (2023). Rational Design, Stabilities and Nonlinear Optical Properties of Non-Conventional Transition Metalides; New Entry into Nonlinear Optical Materials. Materials, 16(9), 3447. https://doi.org/10.3390/ma16093447

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop