Next Article in Journal
Enhancing Performance of Organic Pollutant Degradation via Building Heterojunctions with ZnO Nanowires and Na Doped Conjugated 2,4,6-Triaminopyrimidin-g-C3N4
Previous Article in Journal
High-Throughput Mining of Novel Compounds from Known Microbes: A Boost to Natural Product Screening
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exploring the Structural and Electronic Properties of Niobium Carbide Clusters: A Density Functional Theory Study

by
Hui-Fang Li
1,
Huai-Qian Wang
1,2,* and
Yu-Kun Zhang
1
1
College of Engineering, Huaqiao University, Quanzhou 362021, China
2
College of Information Science and Engineering, Huaqiao University, Xiamen 361021, China
*
Author to whom correspondence should be addressed.
Molecules 2024, 29(13), 3238; https://doi.org/10.3390/molecules29133238
Submission received: 5 June 2024 / Revised: 29 June 2024 / Accepted: 5 July 2024 / Published: 8 July 2024
(This article belongs to the Topic Advances in Computational Materials Sciences)

Abstract

:
This paper systematically investigates the structure, stability, and electronic properties of niobium carbide clusters, NbmCn (m = 5, 6; n = 1–7), using density functional theory. Nb5C2 and Nb5C6 possess higher dissociation energies and second-order difference energies, indicating that they have higher thermodynamic stability. Moreover, ab initio molecular dynamics (AIMD) simulations are used to demonstrate the thermal stability of these structures. The analysis of the density of states indicates that the molecular orbitals of NbmCn (m = 5, 6; n = 1–7) are primarily contributed by niobium atoms, with carbon atoms having a smaller contribution. The composition of the frontier molecular orbitals reveals that niobium atoms contribute approximately 73.1% to 99.8% to NbmCn clusters, while carbon atoms contribute about 0.2% to 26.9%.

Graphical Abstract

1. Introduction

Previous studies have indicated that detailed experimental and theoretical investigations of transition metal clusters often reveal remarkable relationships between their physical properties and chemical reactivity [1,2,3,4]. Niobium clusters are among the most interesting transition metal clusters that have been extensively studied. They possess a long and rich history of experimental and theoretical investigations, attributed to several special properties, including the relative tendency to form clusters, the existence of a single naturally occurring isotope, high melting point, high temperature resistance, and superconductivity [5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24]. The relative propensity to form clusters is experimentally advantageous, and the presence of only one naturally occurring isotope simplifies the interpretation of niobium spectra by avoiding complications from overlapping features caused by different isotopomers. On the experimental side, the ionization potential (IP) [5], dissociation energy (DE) [6,7,8], reactivity [9,10,11], and electron affinity (EA) [12,13,14] have been measured on niobium clusters and their ions. The results have revealed dramatic size-dependent fluctuations of the electronic and chemical properties of niobium clusters. For example, both the reactivity rates of Nbn (n 20) clusters with D2 and the ionization potential show that Nb8, Nb10, and Nb16 are relatively unreactive [15,16]. The strong size dependence of their chemical reactivity further makes them attractive for study, as certain species may possess the required balance of properties for nanotechnology applications. On the theoretical side, niobium has garnered significant attention, with a focus on studying the structural and electronic properties of a series of pure niobium or niobium-related clusters through the application of density functional theory (DFT) [19,20,21]. Many studies on small Nb clusters confirm that Nb8, Nb10, and Nb16 clusters have higher stability than other clusters. Pansini and co-workers studied the electronic structure and electrical properties of Al-doped niobium clusters [21]. By using all-electron density functional theory with Douglas–Kroll–Hess correction, the spin state, geometry, hardness, and mean static dipole polarizability of NbxAly clusters were examined.
Among the studies, niobium carbides, NbmCn, have been the topic of extensive theoretical and experimental research [25,26,27,28,29,30,31,32,33,34,35] due to their potential applications in numerous areas, including catalysis, electronic materials, and energy storage. Mass spectrometry [25,26,27,28], X-ray diffraction spectroscopy [29,30], photoelectron spectroscopy [30,31,32], and electron spin resonance spectroscopy [33], combined with DFT, have been applied to study niobium carbide clusters. The mass spectroscopic studies by Mafuné and co-workers used a double laser ablation technique, and ionization potentials (IPs) of a series of NbnCm (n = 3–10, m = 0–7) were measured and compared with the IPs estimated by density functional theory calculations [28]. Electron spin resonance spectroscopy was performed to study diatomic NbC by Hamrick and Weltner [33]. The elastic properties of niobium carbide NbCy were studied by Gusev and co-workers [34]. Metha and co-workers used photoionization efficiency spectroscopy in combination with DFT to determine the ionization potential of Nb3Cn (n = 1–4) and Nb4Cn (n = 1–6), and Nb3C2 and Nb4C4 clusters exhibited the lowest IPs for the two series, respectively [27].
In this study, we investigate the geometry, stability, and electronic structures of NbmCn (m = 5, 6; n = 1–7) clusters using DFT calculations with the BPW91 functional. Firstly, a global search is performed to obtain the geometric structures of Nb5Cn and Nb6Cn (n = 1–7) clusters. Secondly, the relative stability of the optimized structure is determined by analyzing the average binding energy, dissociation energy, and second-order difference energy, and the stability reasons are explained. Finally, the density of states and compositions of the frontier molecular orbitals are discussed.

2. Results and Discussion

2.1. Geometry

The structures of different isomers for Nb5Cn and Nb6Cn (n = 1–7) are obtained using the BPW91 functional with the 6-311G+(2d) basis set for C and SDD basis set for Nb. Figure 1 displays the lowest-energy geometries and some low-lying isomers of Nb5Cn and Nb6Cn (n = 1–7). According to the total energy from low to high, the isomers of Nb5Cn (n = 1–7) are designated by 5-na, 5-nb, and 5-nc, while the isomers of Nb6Cn (n = 1–7) are designated by 6-na, 6-nb, and 6-nc, where n is the number of C atoms. The Cartesian coordinates of the lowest-energy structures of Nb5Cn and Nb6Cn (n = 1–7) are summarized in Table S1 of the Supplementary Materials. The optimized results show that the most stable structures of Nb5Cn (n = 1–7) exhibit a doublet spin state, whereas those of Nb6Cn (n = 1–7) exhibit a singlet spin state. In the following sections, we discuss these results in detail.

2.1.1. Nb5Cn (n = 1–7)

The lowest-energy structure of Nb5C is found to be a distorted prism with a doublet spin state, which can be viewed as the one obtained from replacing one Nb atom by a C atom in the Nb6 structure. The corresponding quartet states are higher in energy than the doublet state by 0.26 eV. Another two similar structures, 5-1b and 5-1c, are obtained, and their relative energies to the ground state are 0.53 and 0.82 eV, respectively.
Nb5C2 is an irregular pentagonal bipyramid with similar features as the structure of Nb7. The structure of Nb5C2 can be described as the one obtained from replacing two Nb atoms by C atoms sitting on the pentagonal ring. The point group symmetry is C2v and the calculated equatorial bond lengths of the pentagonal ring are 2.04, 2.05, and 2.84 Å, whereas the polar-equatorial bond lengths are 2.19, 2.63, and 2.55 Å. The pole-to-pole bond length is 3.13 Å. Structure 5-2b is also a pentagonal bipyramid with C2v point group symmetry. Comparing structures 5-2a and 5-2b, the two C atoms occupy opposite positions in the pentagonal ring in structure 5-2a and adjacent positions in structure 5-2b. Structure 5-2b lies 0.69 eV above structure 5-2a. Structure 5-2c possesses low Cs symmetry and is less stable than 5-2a by 0.93 eV.
A previous study predicted that Nb5C3 has a similar cubic structure as Nb4C4. The BPW91 calculation determined that Nb5C3 is a distorted cubic structure and has low C1 symmetry. Structure 5-3a consists of two quadrangles. The upper quadrangle contains two C atoms in opposite positions and the lower quadrangle has one C atom. Isomer 5-3b presents a similar geometry as structure 5-3a, with a slight distortion, but with two C atoms occupying adjacent positions. Isomers 5-3b and 5-3c are less stable than 5-3a by 0.30 and 0.53 eV, respectively.
The three low-lying isomers of Nb5C4 listed in this paper are all capped structures, and the main difference between the three isomers is the position of the four C atoms in the cluster. It should be pointed out that the energies of structures 4a and 4b are very close, and the energy of structure 5-4b is only 0.02 eV higher than that of 5-4a. In view of such a small relative energy and in order to ensure the accuracy of calculations, we conducted DFT computations with the Def2-DZVPP basis set instead of the SDD basis set. It was found that the energy difference between 5-4a and 5-4b was very small, but the energy level order was reversed. The energy of structure 5-4b was 0.01 eV lower than that of the 5-4a structure (see Table S2). Additionally, DFT calculations were also performed using the B3LYP functional. The B3LYP functional calculations showed that isomer 5-4b is the lowest energy structure, in which each of the two C atoms occupies the opposite position of the upper and lower quadrangle. In this case, we calculated the vertical ionization potentials (VIPs) of both isomers 5-4a and 5-4b. The results found that the calculated VIP of isomer 5-4b is 5.15 eV (B3LYP) or 5.13 eV (BPW91), and the calculated VIP of isomer 5-4a is 4.85 eV (B3LYP) or 4.90 eV (BPW91). No matter the B3LYP or BPW91 functional, the calculated VIP value for structure 5-4b is in better agreement with the experimental result (5.12 ± 0.12 eV) than for structure 5-4a. Therefore, isomer 5-4b should be the ground state structure. Isomer 5-4c is less stable than 5-4a by 0.15 eV.
From the point of view of growth, the lowest-energy structure of Nb5C5 can be viewed as the one obtained from capping of the C atom to the lowest-energy structure of Nb5C4 with a little local distortion. Another isomer, 5-5b, which possesses similar shape but a different C atom position with isomer 5-5a, is less stable than 5-5a by 0.45 eV. Structure 5-5c, which has relatively high C2v point group symmetry and can be obtained by adding one C atom to the bottom of structure 5-4c, is 1.03 eV higher in energy than structure 5-5a.
The lowest-energy structure for Nb5C6 is found to be a ring, consisting of four C atoms and three Nb atoms capped with one C atom and one Nb atom on each side of the ring. The point group symmetry is C2v. Another two structures, 5-6b and 5-6c, with lower C1 symmetry, are less stable than structure 5-6a by 1.27 and 2.22 eV, respectively. For Nb5C7, when a C atom is added to the lowest-energy structure of 5-6a at different positions, isomers 5-7a, 5-7b, and 5-7c are formed. Isomer 5-7a is more stable than 5-7b and 5-7c by 0.07 and 0.25 eV, respectively.

2.1.2. Nb6Cn (n = 1–7)

For Nb6C, the lowest-energy structure is a pentagonal bipyramid with four Nb atoms and one C atom in a plane, and the other two Nb atoms are located on each side of the pentagonal ring. The energies of the other two structures, 6-1b and 6-1c, are higher than those of structure 6-1a by 0.66 eV and 0.99 eV, respectively. We found the lowest-energy structure of Nb6C2 is approximately C2 point group symmetry. Three Nb atoms and two C atoms make up a ring that is not in the plane, and one Nb atom and a dimer Nb2 occupy each side of the ring, respectively. Another structure, 6-2b, with Cs point group symmetry, is less stable than 6-2a by 0.67 eV. The lowest-energy structure of Nb6C3 can be obtained by capping of an atom to Nb5C3 with different positions of C atoms. Structures 6-3b and 6-3c are less stable than 6-3a by 0.50 and 0.64 eV, respectively. The lowest-energy structure 6-4a has relatively high D2d point group symmetry. The structure can be viewed as the one obtained from capping of a Nb atom to the bottom of structure 5-4c of Nb5C4, with small structure distortion. Structures 6-4b and 6-4c are 0.75 and 0.95 eV higher in energy than structure 6-4a. As for Nb6C5, Nb6C6, and Nb6C7 clusters, three low-lying isomers of each system are listed in this paper; isomer 6-5a is the structure with the lowest energy, and the energies of the other two isomers are 0.24 and 0.43 eV higher than that of structure 6-5a, respectively. Isomer 6-6a is the lowest-energy structure for the Nb6C6 cluster, and isomers 6-6b and 6-6c are less stable than isomer 6-6a by 0.73 and 0.84 eV, respectively. Isomers 6-7b and 6-7c are less stable than the lowest-energy isomer 6-7a by 1.75 and 2.31 eV, respectively.

2.2. Comparison between Calculated and Experimental Vertical Ionization Potential

The VIPs of the stable structures of Nb5Cn and Nb6Cn (n = 1–7), determined with BPW91 calculations, are shown in Table 1 and Figure 2, along with the experimental results [28]. The calculated VIP values are in excellent agreement with the experiments for Nb5Cn (n = 1–6); however, for Nb5C7, the calculated result likely overestimates this parameter, and the relative error is 8.6%. We tried to obtain more of the various structures and calculate the VIPs of different isomers, but the results were not satisfactory. For Nb6Cn (n = 1–7) clusters, the employed DFT method predicted VIPs with acceptable agreement with the experiments. The calculation results showed that the relative error of Nb6Cn (n = 1–7) clusters is within 3.5%. The good agreement between the calculated VIPs and the experimental values gives confidence in the assigned the ground state for the complexes considered in the present paper.

2.3. Relative Stability

To determine the structural stabilities of these Nb5Cn and Nb6Cn (n = 1–7) clusters, the average binding energy per atom (Eb) can be calculated using the following formulas:
Eb (Nb5Cn) = [5E(Nb) + nE(C) − E(Nb5Cn)]/(5 + n)
Eb (Nb6Cn) = [6E(Nb) + nE(C) − E(Nb6Cn)]/(6 + n)
Here, E(Nb), E(C), E(Nb5Cn), and E(Nb6Cn) represent the total energy of Nb, C, Nb5Cn, and Nb6Cn, respectively, where “n” denotes the number of C atoms. The average binding energy is a measurement of the thermodynamic stability of the clusters, and the results are listed in Table 2 and plotted in Figure 3 as a function of the number of C atoms. As shown in Figure 3, the Eb values of these Nb5Cn and Nb6Cn (n = 1–7) clusters are quite large, ranging from approximately 4 to 6 eV, and are close to each other. For Nb5Cn (n = 1–7), the Eb becomes a growing function of the number of C atoms and has two inflection points, corresponding to n = 2 and 6, indicating that Nb5C2 and Nb5C6 clusters are more stable than others. For the Nb6Cn cluster, the Eb shows a monotonic increasing trend with increasing cluster size. Starting from n = 3, the rate of increase in average binding energy slightly decelerates, indicating that Nb6C3 exhibits slightly higher but not significantly pronounced stability.
The stability of Nb5Cn and Nb6Cn (n = 1–7) clusters was further investigated by calculating the dissociation energy (DE), defined as follows:
DE (Nb5Cn) = E (Nb5Cn−1) + E(C) − E(Nb5Cn)
DE (Nb6Cn) = E (Nb6Cn−1) + E(C) − E(Nb6Cn)
The dissociation energies to remove a C atom from the clusters are illustrated in Table 2 and Figure 3. For the Nb5Cn (n = 1–7) cluster, the dissociation energies of Nb5C2 and Nb5C6 corresponding to dissociation paths N b 5 C 2 N b 5 C + C and N b 5 C 6 N b 5 C 5 + C yield the local maxima of all calculated clusters, indicating that Nb5C2 and Nb5C6 clusters are more stable than their neighbors. The dissociation energy of Nb5C3 corresponding to the dissociation path N b 5 C 3 N b 5 C 2 + C yields the local minima of all calculated clusters, indicating that Nb5C3 has relatively weak stability. For Nb6Cn, the dissociation energy of Nb6C3 corresponding to the dissociation path N b 6 C 3 N b 6 C 2 + C yields a small local maximum, indicating slightly higher stability compared to neighboring clusters. This result is consistent with the results of the average binding energy.
The second-order difference (Δ2E) refers to taking the difference again based on the first-order difference of energy. Δ2E, a more sensitive quantity, reflects the relative stability of the clusters. Positive peaks of Δ2E indicate greater stability of the cluster, while a smaller Δ2E suggests weaker stability of the cluster. The Δ2E values of Nb5Cn and Nb6Cn (n = 1–7) clusters were determined using the formulas:
Δ2E (Nb5Cn) = E(Nb5Cn−1) + E(Nb5Cn+1) − 2E(Nb5Cn)
Δ2E (Nb6Cn) = E(Nb6Cn−1) + E(Nb6Cn+1) − 2E(Nb6Cn)
where E(Nb5/6Cn−1), E(Nb5/6Cn), and E(Nb5/6Cn+1) represent the total energy of Nb5/6Cn−1, Nb5/6Cn, and Nb5/6Cn+1, respectively. The Δ2E is presented in Table 2 and plotted in Figure 3 as a function of the number of C atoms. The minima of Nb5Cn (n = 1–7) are found at Nb5C3 and Nb5C5, indicating that these clusters show obviously weak stability. The maxima of Nb5Cn are found at n = 2 and 6, indicating that Nb5C2 and Nb5C6 possess higher stability than their neighbors. The curve for DE shows a similar behavior in Figure 3, confirming the thermodynamic stability of these clusters. For Nb6Cn (n = 1–7), Δ2E initially decreases, then increases, and gradually decreases again from n = 3 to 6. The curve of Δ2E reaches a small local maximum at n = 3, showing a trend similar to the dissociation energy, implying that the stability of Nb6C3 clusters is slightly higher than that of neighboring clusters.
Comparing the analysis results of average binding energy, dissociation energy, and second-order difference energy, it can be observed that Nb5C2 and Nb5C6 exhibit significantly higher stability. The stability of Nb6C3 is slightly higher among Nb6Cn clusters but not notably so. Next, we use ab initio molecular dynamics (AIMD) analysis to explain the reasons for the higher stability of Nb5C2 and Nb5C6. Furthermore, when n = 1 and n = 7, due to the limitations imposed by the cluster sizes computed in this study, Nb5C and Nb5C7 as well as Nb6C and Nb6C7 are positioned at the endpoints of the line. Their stability will also be further discussed in conjunction with AIMD analysis.

2.4. Stability Analysis via AIMD Simulations

In the above discussion, through a comprehensive comparison of the results of average binding energy, dissociation energy, and second-order difference energy, we identified that Nb5C2 and Nb5C6 clusters exhibit higher stability. To analyze the reasons for the stability of these clusters, AIMD simulations, implemented in Born–Oppenheimer MD (BOMD) format, were conducted to analyze their stability. Additionally, AIMD calculations were performed on clusters positioned at the endpoints, namely Nb5C, Nb5C7, Nb6C, and Nb6C7. Each system was set at three different temperatures (100 K, 200 K, and 300 K) under vacuum conditions and maintained a certain temperature stability every 15 fs. The simulations allowed us to observe the thermal motion of the atoms and the geometric fluctuations in these systems. Figure 4 presents the root mean square deviation (RMSD) curves for the three AIMD trajectories at these temperatures. Here, RMSD is defined as the square root of the mean square difference between the position where each atom moves and its initial position. The trajectory from the beginning to the end is represented by three colors: red, white, and blue. The more overlapping of the atoms represented by these three colors, the smaller the range of atomic fluctuations, corresponding to smaller RMSD values and better thermal stability.
As can be seen in Figure 4, by comparing the RMSD of these systems, it is evident that the order of fluctuation amplitude is Nb5C2(Nb5C6) < Nb6C < Nb5C7 < Nb6C7 < Nb5C. Neither isomerization nor dissociation of these systems was observed during 2000 fs simulations at the three temperatures, which implies that the thermal stability should never be underestimated, at least in vacuum and not too high temperature. As the temperature increases, the differences in the ability of different systems to maintain equilibrium gradually manifest. For Nb5C, Nb5C7, and Nb6C7, the fluctuation amplitude significantly increases with the increase in temperature, and no regular vibration period forms. Especially, the RMSD values of Nb5C do not show a convergence trend, indicating that it is likely to undergo dissociation. By contrast, for Nb5C2, Nb5C6, and Nb6C, the variation amplitude with temperature change is not significant (usually less than 0.1 Å), and it has a relatively obvious vibration period (about 250 fs) with a clear convergence trend, indicating that these structures have a certain degree of flexibility and can also maintain good stability.
The reasons for the differences in stability among the different systems can also be analyzed based on the atom colors (red–white–blue) in Figure 5. For Nb5C2 and Nb5C6, they both exhibit a certain degree of symmetry and balanced bond lengths. The horizontal movements have relatively small amplitudes, and the movements perpendicular to the paper surface are almost stationary. This results in these structures having a better ability to resist interference. On the other hand, the other candidates show a reduction or loss of symmetry, leading to a decrease in their ability to resist interference. One example is the loss of partial symmetry, as seen in Nb6C. Although the positioning of the doped atoms on the pentagonal ring distorts the structure of the pentagonal bipyramid, the overall structure remains largely unchanged, allowing it to maintain a certain degree of stability. The second example is the loss of symmetry through the addition of atoms, as observed in Nb5C7 and Nb6C7. The additional atoms have a greater impact on their neighboring atoms, while their impact on other atoms is relatively small, resulting in fluctuations occurring only between a portion of the atoms and minimal overall decreases in stability. The worst case is Nb5C, where the uneven bond lengths and fewer atoms result in weak binding forces between the atoms. Consequently, this leads to the highest volatility at high temperatures.
When comparing Figure 3, Figure 4 and Figure 5, we can observe that Nb5C2, Nb5C6, and Nb6C clusters, which exhibit high thermal stability, also perform well in terms of dissociation energy (Ed) and second-order difference energy (Δ2E). Additionally, Nb5C2 and Nb5C6 correspond to points where the slope of the average binding energy increases rapidly, indicating that these structures are possible magic number clusters.

2.5. Density of States

To better understand the molecular orbitals of Nb5Cn and Nb6Cn (n = 1–7) clusters, the total density of states (TDOS) and the partial densities of states (PDOS) of Nb5/Nb6 and Cn atoms are plotted in Figure 6, Figure 7, Figure 8 and Figure 9 using the BPW91 functional. TDOS refers to the density of electronic states or the number of electronic energy levels per unit energy within a given energy range. It describes the distribution of electronic energy levels within the cluster and can be used to analyze the electronic structure of the material. The calculation of PDOS involves decomposing DOS to obtain contributions from different atomic orbitals or interatomic interactions. The calculation of Nb5Cn (n = 1–7) is based on an unrestricted open-shell treatment, and both alpha-TDOS and beta-TDOS are obtained. For Nb6Cn (n = 1–7), the calculation predicts it to be a closed shell and the spin state is singlet. From these graphs, it can be observed that for both Nb5Cn and Nb6Cn (n = 1–7) series, the PDOS profiles of the C atoms are commonly lower, compared to the PDOS curves of Nb. With the increase in the number of C atoms, the contribution of C atoms to the molecular orbitals increases slightly. That is to say, Nb atoms play an important role during the formation of molecular orbitals compared with the usually smaller contribution of C atoms. For low energy regions (the first peak), C atoms are active and contribute significantly to the molecular orbitals. However, as the energy increases, the contribution of C atoms decreases. In particular, the contribution of C atoms to the frontier molecular orbitals is very small. We provide the compositions of the frontier molecular orbitals (alpha-HOMO, alpha-LUMO, beta-HOMO, beta-LUMO) for the open-shell structures Nb5Cn (n = 1–7) and closed-shell structures Nb6Cn (n = 1–7) in Table 3 and Figure 10 and Figure 11. The results indicated that the HOMOs and LUMOs are primarily derived from the contribution of Nb atoms, and the contribution from C atoms is commonly small. For example, for α-HOMO, 79.7–99.7% is contributed by Nb atoms, while 0.3–20.3% is contributed by C atoms. For α-LUMO, 87.4–99.6% is contributed by Nb atoms, while 0.4–12.6% is contributed by C atoms. For β-HOMO, 73.1–98.2% is contributed by Nb atoms, while 1.8–26.9% is contributed by C atoms. For β-LUMO, 82.1–99.5% is contributed by Nb atoms, while 0.5–17.9% is contributed by C atoms. For the Nb6Cn (n = 1–7) cluster, 83.4–99.8% is contributed to HOMO by Nb atoms, while 0.2–16.6% is contributed to HOMO by C atoms. Additionally, 88.1–99.2% is contributed to LUMO by Nb atoms, while 0.8–11.9% is contributed to LUMO by C atoms.

3. Computational Methods

The ground states of Nb5Cn and Nb6Cn (n = 1–7) clusters were determined using DFT with the generalized gradient approximation by employing GAUSSIAN programs [36]. In recent years, various algorithms and strategies have played an important role in exploring the potential energy surface (PES) of clusters [37,38,39]. We first performed the search for the global minima on the PES for Nb5Cn and Nb6Cn (n = 1–7) systems using the stochastic kicking (SK) method [40,41,42]. In recent years, our research group has effectively predicted the ground-state structures and electronic properties of binary mixed clusters utilizing the SK-DFT method [4,19,35,43]. All of the atoms are placed at the same point initially and then are “kicked” randomly with a sphere of some radius. The kick method was run approximately 500 times at the BPW91/3-21G level until no new minima appeared. Subsequently, several relevant lower-lying isomers were selected for further optimization utilizing the larger basis set. The triply-split basis set with polarization and diffuse functions, 6-311+G(2d) and SDD pseudopotential basis set, which accounts for relativistic effects, were used for C atoms and Nb atoms, respectively. The BPW91 functional was employed in these calculations [44,45], which has been successfully used in the calculations of Nb clusters [46,47,48]. Different possible spin multiplicities were also considered for each of these structural isomers to determine the preferred spin states of these complexes. For Nb5Cn (n = 1–7) clusters, spin multiplicities of 2, 4, 6, and 8 were considered during the calculation process, while for Nb5Cn (n = 1–7) clusters, spin multiplicities of 1, 3, 5, and 7 were considered. Spin-restricted DFT calculations were employed for the singlet state, while spin-unrestricted DFT calculations were employed for all other electronic states. The calculated results showed that the most stable structures of Nb5Cn and Nb6Cn (n = 1–7) prefer the lowest spin state. Vibrational frequency calculations were performed to ascertain the stability of lowest energy isomers. Hence, the structures depicted in Figure 1 represent true local minima as they exhibit positive frequencies. In addition, the low-lying energy isomers identified in Figure 1 were recalculated with higher precision using the Def2-DZVPP basis set instead of the SDD basis set, as shown in Table S2 of the Supplementary Materials. VIP is defined as the energy difference between the cation clusters at optimized neutral geometry clusters and optimized neutral clusters: VIP = E(cation clusters at optimized neutral geometry) − E(optimized neutral). AIMD simulations were performed with ORCA 5.0 code [49], and the dynamic trajectories were visualized by Visual Molecular Dynamics (VMD 1.9.3) software [50]. In the simulations, the BPW91 functional was selected, the step size was set to 1 fs, and a CSVR thermostat with a time constant of 15 fs was employed to maintain the temperature. Density of states and orbital composition analysis were performed by the Multiwfn program [51] and visualized by VMD 1.9.3 software.

4. Conclusions

The geometries, stabilities, and electronic properties of niobium carbide clusters, NbmCn (m = 5, 6; n = 1–7), have been investigated within the framework of density functional theory. The average binding energy, dissociation energy, and second-order difference energy, combined with AIMD simulations, are used to discuss the stability of NbmCn (m = 5, 6; n = 1–7) clusters. The results reveal that Nb5C2, Nb5C6, and Nb6C exhibit relatively higher thermodynamic stability among all clusters investigated in this study. In addition, the molecular orbitals of NbmCn (m = 5, 6; n = 1–7) are primarily contributed by niobium atoms, and niobium atoms contribute approximately 73.1–99.8% to NbmCn (m = 5, 6; n = 1–7) clusters. By contrast, carbon atoms make a modest contribution, ranging from 0.2% to 26.9%.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules29133238/s1, Table S1: Cartesian coordinates for the lowest energy structure of NbmCn (m = 5,6; n =1–7) at the BPW91/Nb/SDD//C/6-311+G(2d) level; Table S2: The relative energies of NbmCn (m = 5,6; n =1–7) at the BPW91/Nb/SDD//C/6-311+G(2d) and BPW91/Nb/Def2-TZVPP//C/6-311+G(2d) level.

Author Contributions

Conceptualization, H.-F.L. and H.-Q.W.; methodology, H.-F.L. and H.-Q.W.; software, H.-Q.W.; validation, H.-F.L. and H.-Q.W.; formal analysis, H.-Q.W. and Y.-K.Z.; investigation, H.-F.L.; resources, H.-Q.W.; data curation, H.-F.L.; writing—original draft preparation, H.-F.L.; writing—review and editing, H.-F.L., H.-Q.W. and Y.-K.Z.; visualization, H.-F.L.; supervision, H.-Q.W.; project administration, H.-Q.W.; funding acquisition, H.-Q.W. All authors have read and agreed to the published version of the manuscript.

Funding

The project was supported by the Natural Science Foundation of Xiamen (Grant No. 3502Z202373051) and the Natural Science Foundation of Fujian Province of China (Grant No. 2023J01141).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available upon request from the corresponding authors.

Acknowledgments

The authors would like to express their gratitude to the School of Engineering, Huaqiao University, which is affiliated with Huaqiao University. This research was conducted at the School of Engineering, Huaqiao University.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Alonso, J.A. Electronic and atomic structure, and magnetism of transition-metal clusters. Chem. Rev. 2000, 100, 637–678. [Google Scholar] [CrossRef] [PubMed]
  2. Knickelvein, M.B. Experimental observation of superparamagnetism in manganese clusters. Phys. Rev. Lett. 2001, 86, 5255–5257. [Google Scholar] [CrossRef] [PubMed]
  3. Xu, X.S.; Yin, S.Y.; Moro, R.; De Heer, W.A. Magnetic moments and adiabatic magnetization of free cobalt clusters. Phys. Rev. Lett. 2005, 95, 237209. [Google Scholar] [CrossRef] [PubMed]
  4. Fan, Y.W.; Kong, X.Y.; Zhao, L.J.; Wang, H.Q.; Li, H.F.; Zhan, Q.; Xie, B.; Xu, H.G.; Zheng, W.J. A joint experimental and theoretical study on structural, electronic, and magnetic properties of MnGen (n = 3–14) clusters. J. Chem. Phys. 2021, 154, 204302. [Google Scholar] [CrossRef] [PubMed]
  5. Brock, L.; Duncan, M.A. Near-threshold photoionization to probe neutral “met-car” clusters. J. Phys. Chem. 1996, 100, 5654–5659. [Google Scholar] [CrossRef]
  6. Loh, S.K.; Lian, L.; Armentrout, P.B. Oxidation reactions at variably sized transition metal centers: Fen+ and Nbn+ + O2 (n = 1–3). J. Chem. Phys. 1989, 91, 6148–6163. [Google Scholar] [CrossRef]
  7. Hales, D.A.; Lian, L.; Armentrout, P.B. Collision-induced dissociation of Nbn+ (n = 2–11): Bond energies and dissociation pathways. Int. J. Mass Spectrom. Ion Process. 1990, 102, 269–301. [Google Scholar] [CrossRef]
  8. Cartier, S.F.; May, B.D.; Castleman, A.W., Jr. The delayed ionization and atomic ion emission of binary metal metallocarbohedrenes TixMyC12 (M = Zr, Nb; 0 ≤ y ≤ 4; x + y = 8). J. Chem. Phys. 1996, 104, 3423–3432. [Google Scholar] [CrossRef]
  9. Deng, H.T.; Guo, B.C.; Kerns, K.P.; Castleman, A.W., Jr. Gas phase reactions of the met-cars Ti8C12+, Nb8C12+ and Ti7Nb12+ with acetone and methyl lodide. J. Phys. Chem. 1994, 98, 13373–13378. [Google Scholar] [CrossRef]
  10. Yeh, C.S.; Byun, Y.G.; Afzaal, S.; Kan, S.Z.; Lee, S.; Freiser, B.S.; Hay, P.J. Experimental and theoretical studies on Nb4C40/+: Reactivity and structure of the smallest cubic niobium-carbon cluster. J. Am. Chem. Soc. 1995, 117, 4042–4048. [Google Scholar] [CrossRef]
  11. Byun, Y.G.; Lee, S.A.; Kan, S.Z.; Freiser, B.S. Reactivities of metallocarbohedrenes: Nb8C12+. J. Phys. Chem. 1996, 100, 14281–14288. [Google Scholar] [CrossRef]
  12. Yang, D.S.; Zgierski, M.Z.; Berces, A.; Hackett, P.A.; Roy, P.N.; Martinez, A.; Carrington, T.; Salahub, D.R.; Fournier, R.; Pang, T.; et al. Vibrational and geometric structures of Nb3C2 and Nb3C2+ from pulsed field ionization-zero electron kinetic energy photoelectron spectra and density functional calculations. J. Chem. Phys. 1996, 105, 10663–10671. [Google Scholar] [CrossRef]
  13. Li, S.; Wu, H.; Wang, L. Probing the electronic structure of metallocarbohedrenes: M8C12 (M = Ti, V, Cr, Zr, Nb). J. Am. Chem. Soc. 1997, 119, 7417–7422. [Google Scholar] [CrossRef]
  14. Knappenberger, K.L., Jr.; Clayborne, P.A.; Reveles, J.U.; Sobhy, M.A.; Jones, C.E.; Gupta, U.U.; Khanna, S.N.; Iordanov, I.; Sofo, J.; Castleman, A.W., Jr. Anion photoelectron spectroscopy and density functional investigation of diniobium-carbon clusters. ACS Nano 2007, 1, 319–326. [Google Scholar] [CrossRef] [PubMed]
  15. Geusic, M.E.; Morse, M.D.; Smalley, R.E. Hydrogen chemisorption on transition metal clusters. J. Chem. Phys. 1985, 82, 590–591. [Google Scholar] [CrossRef]
  16. Whetten, R.L.; Zakin, M.R.; Cox, D.M.; Trevor, D.J.; Kaldor, A. Electron binding and chemical inertness of specific Nbx clusters. J. Chem. Phys. 1986, 85, 1697–1698. [Google Scholar] [CrossRef]
  17. Eduok, U. Niobia Nanofiber-Reinforced Protective Niobium Oxide/Acrylate Nanocomposite Coatings. ACS Omega 2020, 5, 30716–30728. [Google Scholar] [CrossRef] [PubMed]
  18. Korzyński, M.D.; Xie, L.S.; Dincă, M. Structural Characterization of a High-Nuclearity Niobium(V) Carboxylate Cluster Based on Pivalic Acid. Helv. Chim. Acta 2020, 103, e2000186. [Google Scholar] [CrossRef]
  19. Wang, H.Q.; Li, H.F. Probing the structural and electronic properties of al-doped small niobium clusters. Chem. Phys. Lett. 2012, 554, 231–235. [Google Scholar] [CrossRef]
  20. Grönbeck, H.; Rosén, A.; Andreoni, W. Structural electronic, and vibrational properties of neutral and charged Nbn (n = 8, 9, 10) clusters. Phys. Rev. A 1998, 58, 4630–4636. [Google Scholar] [CrossRef]
  21. Pansinia, F.N.N.; Camposb, M.D.; Netoa, A.C.; Sergiob, C.S. Theoretical study of the electronic structure and electrical properties of Al-doped niobium clusters. Chem. Phys. 2020, 535, 110778. [Google Scholar] [CrossRef]
  22. Prudnikava, A.; Tamashevich, Y.; Babenkov, S.; Makarova, A.; Smirnov, D.; Aristov, V.; Molodtsova, O.; Kugeler, O.; Viefhaus, J.; Foster, B. Systematic study of niobium thermal treatments for superconducting radio frequency cavities employing x-ray photoelectron spectroscopy. Supercond. Sci. Technol. 2022, 35, 065019. [Google Scholar] [CrossRef]
  23. Oliveira, L.; Pereira, M.; Heitman, A.P.; Filho, J.; Olivira, C.; Ziolek, M. Niobium: The focus on catalytic application in the conversion of biomass and biomass derivatives. Molecules 2023, 28, 1527. [Google Scholar] [CrossRef]
  24. Luo, P.; Zhao, Y.H. Niobium nitride preparation for superconducting single-photon detectors. Molecules 2023, 28, 6200. [Google Scholar] [CrossRef] [PubMed]
  25. Wei, S.; Guo, B.C.; Deng, H.T.; Kerns, K.; Purnell, J.; Buzza, S.A.; Castleman, A.W., Jr. Formation of met-cars and face-centered cubic structures: Thermodynamically or kinetically controlled? J. Am. Chem. Soc. 1994, 116, 4475–4476. [Google Scholar] [CrossRef]
  26. Pilgrim, J.S.; Brock, L.R.; Duncan, M.A. Photodissociation of niobium-carbon clusters and nanocrystals. J. Phys. Chem. 1995, 99, 544–550. [Google Scholar] [CrossRef]
  27. Dryza, V.; Addicoat, M.A.; Gascooke, J.R.; Buntine, M.A.; Metha, G.F. Threshold photoionization and density functional theory studies of the niobium carbide clusters Nb3Cn (n = 1–4) and Nb4Cn (n = 1–6). J. Phys. Chem. A 2008, 112, 5582–5592. [Google Scholar] [CrossRef] [PubMed]
  28. Fukushima, N.; Miyajima, K.; Mafune, F. Ionization energies of niobium carbide clusters NbnCm (n = 3–10, m = 0–7). J. Phys. Chem. A 2009, 113, 2309–2315. [Google Scholar] [CrossRef]
  29. Gusev, A.I. Anisotropy of microstructure and elastic properties of niobium carbide nanopowders. Solid State Sci. 2020, 100, 106092. [Google Scholar] [CrossRef]
  30. Chebanenko, M.I.; Danilovich, D.P.; Lobinsky, A.A.; Popkov, V.I.; Rempel, A.A.; Valeeva, A.A. Novel high stable electrocatalyst based on non-stoichiometric nanocrystalline niobium carbide toward effective hydrogen evolution. Int. J. Hydrogen Energy 2021, 46, 16907–16916. [Google Scholar] [CrossRef]
  31. Zhai, H.J.; Liu, S.R.; Wang, L.S. Photoelectron spectroscopy of mono-niobium carbide clusters NbCn (n= 2–7): Evidence for a cyclic to linear structural transition. J. Chem. Phys. 2001, 115, 5170–5178. [Google Scholar] [CrossRef]
  32. Il’in, E.G.; Parshakov, A.S.; Teterin, Y.A.; Maslakov, K.I.; Teterin, A.Y. Surface morphology and composition of a NbC/C composite studied by scanning electron microscopy and X-ray photoelectron spectroscopy. Inorg. Mater. 2020, 56, 443–450. [Google Scholar] [CrossRef]
  33. Hamrick, Y.M.; Weltner, W., Jr. Quenching of angular momentum in the ground states of VC, NbC, VSi, and NbSi molecules. J. Chem. Phys. 1991, 94, 3371–3380. [Google Scholar] [CrossRef]
  34. Valeeva, A.A.; Gusev, A.I. Effect of nonstoichiometry on elastic properties of niobium carbide NbCy. Int. J. Refract. Met. H. 2021, 95, 105435. [Google Scholar] [CrossRef]
  35. Li, H.F.; Wang, H.Q.; Zhang, J.M.; Qin, L.X.; Zheng, H.; Zhang, Y.H. Investigation of structures, stabilities, and electronic and magnetic properties of niobium carbon clusters Nb7Cn (n = 1–7). Molecules 2024, 29, 1692. [Google Scholar] [CrossRef] [PubMed]
  36. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G.A.; et al. Gaussian 09, Revision C.01; Gaussian Inc.: Wallingford, CT, USA, 2010. [Google Scholar]
  37. Diego, I.; Alejandro, V.S.; Luis, L.P.; Williams, G.A.; Maria, L.C.; Osvaldo, Y.; William, T. Revisiting the potential-energy surface of CnBe3n+2H2n+22+ (n = 2–4) clusters: Are planar pentacoordinate carbon structures the global minima? Phys. Chem. Chem. Phys. 2023, 25, 20235–20240. [Google Scholar]
  38. Venkatesan, S.T.; Aland, S.; Diego, I.; Pothiappan, V.; Krishnan, T.; Saikat, R.; Anakuthil, A.; William, T. Why an integrated approach between search algorithms and chemical intuition is necessary? Phys. Chem. Chem. Phys. 2022, 24, 11680–11686. [Google Scholar]
  39. Osvaldo, Y.; Diego, I.; Brandon, U.A.; Alejandro, V.E.; Ricardo, P.R.; Mauricio, T.S.; Jorge, G.; Jorge, B.; Gabriel, M.; William, T. Evaluation of restricted probabilistic cellular automata on the exploration of the potential energy surface of Be6B11. Theor. Chem. Acc. 2020, 139, 41. [Google Scholar]
  40. Roy, D.; Corminboeuf, C.; Wannere, C.S.; King, R.B.; Schleyer, P.V.R. Planar tetracoordinate carbon atoms centered in bare four-membered rings of late transition metals. Inorg. Chem. 2006, 45, 8902–8906. [Google Scholar] [CrossRef]
  41. Bera, P.P.; Sattelmeyer, K.W.; Saunders, M.; Schaefer, H.F.; Schleyer, P.R. Mindless chemistry. J. Phys. Chem. A 2006, 110, 4287–4290. [Google Scholar] [CrossRef]
  42. Saunders, M. Stochastic search for isomers on a quantum mechanical surface. J. Comput. Chem. 2004, 25, 621–626. [Google Scholar] [CrossRef] [PubMed]
  43. Li, H.F.; Wang, H.Q. Stabilization of golden cages by encapsulation of a single transition metal atom. R. Soc. Open Sci. 2018, 5, 171019. [Google Scholar] [CrossRef] [PubMed]
  44. Becke, A.D. Density-functional exchange-energy approximation with correct asymptotic behavior. Phys. Rev. A 1988, 38, 3098–3100. [Google Scholar] [CrossRef]
  45. Perdew, J.P.; Wang, Y. Accurate and simple analytic representation of the electron-gas correlation energy. Phys. Rev. B 1992, 45, 13244–13249. [Google Scholar] [CrossRef]
  46. Huang, B.B.; Zhang, H.Y.; Gan, W.; Yang, M.Z.; Luo, Z.X.; Yao, J.N. Nb12+—Niobespherene: A full-metal hollow-cage cluster with superatomic stability and resistance to CO attack. Natl. Sci. Rev. 2023, 10, nwac197. [Google Scholar] [CrossRef] [PubMed]
  47. Pham, V.N.; Vu, T.N.; Truong, B.T.; Minh, T.N. Electronic Structures, Vibrational and Thermochemical Properties of Neutral and Charged Niobium Clusters Nbn, n = 7–12. J. Phys. Chem. A 2011, 115, 3523–3535. [Google Scholar]
  48. Pham, V.N.; Vu, T.N.; Minh, T.N. A New Look at the Structure and Vibrational Spectra of Small Niobium Clusters and Their Ions. J. Phys. Chem. C 2010, 114, 13210–13218. [Google Scholar]
  49. Neese, F. Software update: The ORCA program system-Version 5.0. Wiley Isnterdiscip. Rev. Comput. Mol. Sci. 2022, 12, e1606. [Google Scholar] [CrossRef]
  50. Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual molecular dynamics. J. Mol. Graph. 1996, 14, 33–38. [Google Scholar] [CrossRef]
  51. Lu, T.; Chen, F.W. Multiwfn: A multifunctional wavefunction analyzer. J. Comput. Chem. 2012, 33, 580–592. [Google Scholar] [CrossRef]
Figure 1. Geometric structures, point group symmetries, as well as relative energies of Nb5Cn and Nb6Cn (n = 1–7) clusters. The relative energies are in eV.
Figure 1. Geometric structures, point group symmetries, as well as relative energies of Nb5Cn and Nb6Cn (n = 1–7) clusters. The relative energies are in eV.
Molecules 29 03238 g001
Figure 2. The experimental and calculated VIPs of Nb5Cn and Nb6Cn (n = 1–7) clusters.
Figure 2. The experimental and calculated VIPs of Nb5Cn and Nb6Cn (n = 1–7) clusters.
Molecules 29 03238 g002
Figure 3. The average binding energy per atom (Eb), dissociation energy (DE), and second-order difference energy ( Δ 2 E ) for Nb5Cn and Nb6Cn (n = 1–7) clusters.
Figure 3. The average binding energy per atom (Eb), dissociation energy (DE), and second-order difference energy ( Δ 2 E ) for Nb5Cn and Nb6Cn (n = 1–7) clusters.
Molecules 29 03238 g003
Figure 4. RMSD of (a) Nb5C; (b) Nb5C2; (c) Nb5C6; (d) Nb5C7; (e) Nb6C; and (f) Nb6C7 at three different temperatures. The trajectories have been aligned to the first frame prior to the RMSD calculation.
Figure 4. RMSD of (a) Nb5C; (b) Nb5C2; (c) Nb5C6; (d) Nb5C7; (e) Nb6C; and (f) Nb6C7 at three different temperatures. The trajectories have been aligned to the first frame prior to the RMSD calculation.
Molecules 29 03238 g004
Figure 5. AIMD simulation trajectories of (a) Nb5C; (b) Nb5C2; (c) Nb5C6; (d) Nb5C7; (e) Nb6C; and (f) Nb6C7. The left and right sides of each small image represent the trajectories at 100 K and 300 K, respectively. The structures are extracted every 100 fs, the color corresponds to the time step and varies as red–white–blue.
Figure 5. AIMD simulation trajectories of (a) Nb5C; (b) Nb5C2; (c) Nb5C6; (d) Nb5C7; (e) Nb6C; and (f) Nb6C7. The left and right sides of each small image represent the trajectories at 100 K and 300 K, respectively. The structures are extracted every 100 fs, the color corresponds to the time step and varies as red–white–blue.
Molecules 29 03238 g005
Figure 6. The total density of state (TDOS) and partial density of state (PDOS) of Nb5Cn (n = 1–4) with a full width at half maximum (FWHM) of 0.5 eV.
Figure 6. The total density of state (TDOS) and partial density of state (PDOS) of Nb5Cn (n = 1–4) with a full width at half maximum (FWHM) of 0.5 eV.
Molecules 29 03238 g006
Figure 7. The total density of state (TDOS) and partial density of state (PDOS) of Nb5Cn (n = 5–7) with a full width at half maximum (FWHM) of 0.5 eV.
Figure 7. The total density of state (TDOS) and partial density of state (PDOS) of Nb5Cn (n = 5–7) with a full width at half maximum (FWHM) of 0.5 eV.
Molecules 29 03238 g007
Figure 8. The total density of state (TDOS) and partial density of state (PDOS) of Nb6Cn (n = 1–4) with a full width at half maximum (FWHM) of 0.5 eV.
Figure 8. The total density of state (TDOS) and partial density of state (PDOS) of Nb6Cn (n = 1–4) with a full width at half maximum (FWHM) of 0.5 eV.
Molecules 29 03238 g008
Figure 9. The total density of state (TDOS) and partial density of state (PDOS) of Nb6Cn (n = 5–7) with a full width at half maximum (FWHM) of 0.5 eV.
Figure 9. The total density of state (TDOS) and partial density of state (PDOS) of Nb6Cn (n = 5–7) with a full width at half maximum (FWHM) of 0.5 eV.
Molecules 29 03238 g009
Figure 10. Molecular orbitals (isosurfaces = 0.05) and energy levels of Nb5Cn (n = 1–7).
Figure 10. Molecular orbitals (isosurfaces = 0.05) and energy levels of Nb5Cn (n = 1–7).
Molecules 29 03238 g010
Figure 11. Molecular orbitals (isosurfaces = 0.05) and energy levels of Nb6Cn (n = 1–7).
Figure 11. Molecular orbitals (isosurfaces = 0.05) and energy levels of Nb6Cn (n = 1–7).
Molecules 29 03238 g011
Table 1. The vertical ionization potential (VIP) for Nb5Cn and Nb6Cn (n = 1–7) clusters, all energies are in eV, the values in parentheses are relative error (%).
Table 1. The vertical ionization potential (VIP) for Nb5Cn and Nb6Cn (n = 1–7) clusters, all energies are in eV, the values in parentheses are relative error (%).
IsomersVIPExpt. [28]IsomersVIPExpt. [28]
5-1a5.30 (1.0%)5.25 ± 0.086-1a5.11 (2.7%)5.25 ± 0.07
5-2a4.45 (3.3%)4.6 ± 0.156-2a5.20 (0.2%)5.21 ± 0.07
5-3a4.60 (1.5%)4.53 ± 0.056-3a4.99 (3.5%)4.82 ± 0.06
5-4b5.13 (0.2%)5.12 ± 0.126-4a4.94 (0.8%)4.9 ± 0.1
5-5a5.02 (0.6%)5.05 ± 0.106-5a5.18 (1.6%)5.1 ± 0.1
5-6a5.00 (0.4%)5.02 ± 0.086-6a5.28 (2.5%)5.15 ± 0.12
5-7a5.54 (8.6%) 5.10 ± 0.126-7a5.28 (3.5%)5.1 ± 0.12
Table 2. The average binding energy per atom (Eb), dissociation energy (DE), and second-order difference energy (Δ2E) for Nb5Cn and Nb6Cn (n = 1–7) clusters, all energies are in eV.
Table 2. The average binding energy per atom (Eb), dissociation energy (DE), and second-order difference energy (Δ2E) for Nb5Cn and Nb6Cn (n = 1–7) clusters, all energies are in eV.
IsomersEbDEΔ2EIsomersEbDEΔ2E
5-1a4.4087.956−0.3846-1a4.6398.9521.652
5-2a4.9708.3401.1516-2a4.9717.300−0.160
5-3a5.2477.190−0.3336-3a5.2487.4590.243
5-4b5.5007.523−0.0946-4a5.4447.2160.111
5-5a5.7127.617−0.6286-5a5.5957.105−0.279
5-6a5.9428.2452.1036-6a5.7457.384−1.232
5-7a5.9596.142 6-7a5.9658.617
Table 3. The compositions of the frontier molecular orbitals for open-shell structures of Nb5Cn and closed-shell structures of Nb6Cn (n = 1–7).
Table 3. The compositions of the frontier molecular orbitals for open-shell structures of Nb5Cn and closed-shell structures of Nb6Cn (n = 1–7).
IsomerAtomαβIsomerAtomHOMOLUMO
HOMOLUMOHOMOLUMO
5-1aNb97.9%99.5%96.3%98.2%6-1aNb99.8%99.2%
C2.1%0.5%3.7%1.8%C0.2%0.8%
5-2aNb99.7%99.6%98.2%99.5%6-2aNb97.4%97.8%
C0.3%0.4%1.8%0.5%C2.6%2.2%
5-3aNb96.8%98.4%93.5%92.0%6-3aNb94.4%94.5%
C3.2%1.6%6.5%8.0%C5.6%5.5%
5-4bNb93.8%90.7%95.4%91.9%6-4aNb84.3%98.3%
C6.2%9.3%4.6%8.1%C15.7%1.7%
5-5aNb92.6%94.9%73.1%86.2%6-5aNb83.4%90.8%
C7.4%5.1%26.9%13.8%C16.6%9.2%
5-6aNb91.0%92.3%94.7%90.2%6-6aNb84.2%93.4%
C9.0%7.7%5.3%9.8%C15.8%6.6%
5-7aNb79.7%87.4%84.7%82.1%6-7aNb88.8%88.1%
C20.3%12.6%15.3%17.9%C11.2%11.9%
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

Li, H.-F.; Wang, H.-Q.; Zhang, Y.-K. Exploring the Structural and Electronic Properties of Niobium Carbide Clusters: A Density Functional Theory Study. Molecules 2024, 29, 3238. https://doi.org/10.3390/molecules29133238

AMA Style

Li H-F, Wang H-Q, Zhang Y-K. Exploring the Structural and Electronic Properties of Niobium Carbide Clusters: A Density Functional Theory Study. Molecules. 2024; 29(13):3238. https://doi.org/10.3390/molecules29133238

Chicago/Turabian Style

Li, Hui-Fang, Huai-Qian Wang, and Yu-Kun Zhang. 2024. "Exploring the Structural and Electronic Properties of Niobium Carbide Clusters: A Density Functional Theory Study" Molecules 29, no. 13: 3238. https://doi.org/10.3390/molecules29133238

APA Style

Li, H. -F., Wang, H. -Q., & Zhang, Y. -K. (2024). Exploring the Structural and Electronic Properties of Niobium Carbide Clusters: A Density Functional Theory Study. Molecules, 29(13), 3238. https://doi.org/10.3390/molecules29133238

Article Metrics

Back to TopTop