Electrochemical Performance of Aluminum Doped Ni_1−xAl_xCo₂O₄ Hierarchical Nanostructure: Experimental and Theoretical Study

Abstract

For electrochemical supercapacitors, nickel cobaltite (NiCo₂O₄) has emerged as a new energy storage material. The electrocapacitive performance of metal oxides is significantly influenced by their morphology and electrical characteristics. The synthesis route can modulate the morphological structure, while their energy band gaps and defects can vary the electrical properties. In addition to modifying the energy band gap, doping can improve crystal stability and refine grain size, providing much-needed surface area for high specific capacitance. This study evaluates the electrochemical performance of aluminum-doped Ni_1−xAl_xCo₂O₄ (0 ≤ x ≤ 0.8) compounds. The Ni_1−xAl_xCo₂O₄ samples were synthesized through a hydrothermal method by varying the Al to Ni molar ratio. The physical, morphological, and electrochemical properties of Ni_1−xAl_xCo₂O₄ are observed to vary with Al³⁺ content. A morphological change from urchin-like spheres to nanoplate-like structures with a concomitant increase in the surface area, reaching up to 189 m²/g for x = 0.8, was observed with increasing Al³⁺ content in Ni_1−xAl_xCo₂O₄. The electrochemical performance of Ni_1−xAl_xCo₂O₄ as an electrode was assessed in a 3M KOH solution. The high specific capacitance of 512 F/g at a 2 mV/s scan rate, 268 F/g at a current density of 0.5 A/g, and energy density of 12.4 Wh/kg was observed for the x = 0.0 sample, which was reduced upon further Al³⁺ substitution. The as-synthesized Ni_1−xAl_xCo₂O₄ electrode exhibited a maximum energy density of 12.4 W h kg⁻¹ with an outstanding high-power density of approximately 6316.6 W h kg⁻¹ for x = 0.0 and an energy density of 8.7 W h kg⁻¹ with an outstanding high-power density of approximately 6670.9 W h kg⁻¹ for x = 0.6. The capacitance retention of 97% and 108.52% and the Coulombic efficiency of 100% and 99.24% were observed for x = 0.0 and x = 0.8, respectively. First-principles density functional theory (DFT) calculations show that the band-gap energy of Ni_1−xAl_xCo₂O₄ remained largely invariant with the Al³⁺ substitution for low Al³⁺ content. Although the capacitance performance is reduced upon Al³⁺ doping, overall, the Al³⁺ doped Ni_1−xAl_xCo₂O₄ displayed good energy, powder density, and retention performance. Thus, Al³⁺ could be a cost-effective alternative in replacing Ni with the performance trade off.

1. Introduction

Researchers are actively engineering an efficient, clean, and cost-effective electrode for better electrochemical performance [1]. Conventional capacitors and batteries are broadly exploited by energy storage devices [2]. Previously, secondary batteries were suitable electrical energy storage devices, but they were not very satisfactory because of their smaller power density, presence of toxic waste substances, and shorter lifetime [3]. On the other hand, conventional capacitors display high power density compared to a battery, but because of their lower energy density, researchers are interested in advanced energy storing devices that mitigate the demerits of both battery and capacitor. The energy gap between batteries and capacitors can be bridged by a supercapacitor (SC) [4]. SC electrode materials must have a large specific surface area, controlled porosity, high electronic conductivity, suitable electroactive sites, high thermal stability, low raw material costs, and ease of manufacture [5]. Supercapacitors have recently attracted attention due to their quick charging/discharging, high power density, wide range of working temperatures, high efficiency, and extended lifespans [6]. SCs are of three types: (i) electrical double layer capacitance (EDLCs) examples are highly porous carbon materials, (ii) pseudocapacitor (PCs) examples are transition metal oxides and conducting polymers, and (iii) hybrid capacitor examples are battery type hybrid devices such as the lithium-ion capacitor [7]. EDLCs store electrical energy in the electric double layer, which occurs at the interface between the electrode’s active material and the electrolyte by an electrostatic accumulation process. In EDLCs, energy density becomes lower (between 3 and 6 Wh/Kg), and specific capacitance is restricted because of the specific surface area of the active electrode material [8]. The ionic double layer’s fast diffusion provides a very rapid discharge profile and hence high-power densities. PCs, on the other hand, are caused by reversible Faradaic processes at the electrode surface, which cause charge storage [9]. Moreover, in PCs, the electrochemical process occurs near the surface of the active electrode material. Thus, it can give specific capacitance 10–100 times greater than that provided by EDLCs [10]. Moreover, hybrid capacitors are a combination of both EDLCs and PCs. Usually, hybrid capacitor electrodes are made from (1) both EDLCs materials and PCs materials; symmetric geometry (2) with one EDLCs electrode and other PCs or battery type electrode; and asymmetric geometry (3) with an asymmetric structure with one PC electrode and another rechargeable battery type electrode [7].

In a pseudocapacitor, a redox reaction takes place between the active electrode material and the electrolyte. Various kinds of transition metal oxides having higher theoretical specific capacitance, such as RuO₂ [11], MnO₂ [12], NiO [13], Ni(OH)₂ [14], CuO [15], Co₃O₄ [16], and Co(OH)₂ [17], have been investigated as an active electrode material. However, their poor electronic conductivity [18,19], lower cyclic performance [20,21], and substantial volume change during the charge/discharge process [22] demand the development of nanostructured materials to meet the increasing and urgent demands for energy storing devices. Recently, mixed-transition-metal-oxides (MTMOs), A_xB_3−xO₄ materials have emerged with promising applications for electrodes. The MTMOs with A and B metal ions can provide many oxidation states for efficient redox charge reactions [23]. Nickel cobaltite (NiCo₂O₄) is one of the most valuable electrode materials studied because of its high theoretical specific capacitance, lower cost, low toxic nature, good redox activity, higher electrical conductivity, and environmentally benign nature [24]. NiCo₂O₄ offers high specific capacitance due to the presence of Ni and Co ions together, which are responsible for producing higher electrochemical activity and electrical conductivity [25]. NiCo₂O₄ is a spinel-type, AB₂O₄, binary metal oxide where nickel cations occupy octahedral sites and cobalt cations are randomly distributed on both tetrahedral and octahedral sites [26].

Doped NiCo₂O₄ with various transition metals, e.g., Ni_1−xRu_xCo₂O₄ [27], Ni_1−xMn_xCo₂O₄ [28], Ni_1−xCr_xCo₂O₄ [29], Ni_1−xZn_xCo₂O₄ [30], Ni_1−xFe_xCo₂O₄ [31], and Ni_1−xCa_xCo₂O₄ [32], have been studied to understand the effect on morphology and their electrochemical performance. This paper describes a simple hydrothermal method for making Al³⁺ doped NiCo₂O_4, i.e., Ni_1−xAl_xCo₂O₄. The study systematically assesses the effect of Al³⁺ doping on the crystal structure, morphology, chemical composition, and electrochemical performance of Ni_1−xAl_xCo₂O₄. Aluminum ion substitution is favored because of its inexpensive cost, large abundance, non-toxic nature, and superior corrosion resistance. First-principles density functional theory (DFT) calculations were employed to assess the energy band-gap of the doped Ni_1−xAl_xCo₂O₄.

2. Experimental

2.1. Synthesis

The required high purity nitrate salts, viz. Al(NO₃)₂.6H₂O, Co(NO₃)₂.6H₂O, and Ni(NO₃)₂.6H₂O, were purchased from Sigma Aldrich (St. Louis, MO, USA).

Table 1. Stoichiometry of chemicals used in the synthesis of Ni_1−xAl_xCo₂O₄.

Ni1−xAlxCo2O4, x	Co(NO3)2·6H2O (gm)	Ni(NO3)2·6H2O (gm)	Al(NO3)2·6H2O (gm)	Urea (gm)
0.0	2.420	1.209	0.000	1.497
0.2	2.485	0.993	0.320	1.497
0.4	2.550	0.765	0.658	1.497
0.6	2.620	0.525	1.016	1.497
0.8	3.850	0.384	1.985	1.497

lists the required amount of precursor salts needed to prepare Ni_1−xAl_xCo₂O₄. The nitrate salts were dissolved in 35 mL of water with a fixed amount of Urea (CO(NH₂)₂) and stirred for 30 min. The solution was then transferred to a 45 mL Teflon-lined autoclave and maintained at 190 °C for 12 h, followed by a natural cooling to room temperature. The precipitates were washed with distilled water and acetone several times by centrifugation and then dried for 12 h at 80 °C. After this, the as derived metal hydroxide composite particles were calcined at 350 °C in the air for 3.5 h.

2.2. Characterization

X-ray diffraction (XRD) was used to examine the crystalline structure of calcined powder using a Bruker D8 Advance X-ray diffractometer (Bruker, Billerica, MA, USA). The surface morphology of prepared samples was studied by a scanning electron microscope (SEM) (at 10 keV using Phenom, Oak Park, CA, USA). Autosorb (Quantachrome, Boynton Beach, FL 33426, USA, model No. AS1MP) was used to quantify specific surface area and pore volume at 77K using nitrogen as the adsorbing gas. The Brunauer–Emmett–Teller (BET) method was used to measure the specific surface area of as-prepared samples. Also, surface area, average pore volume, and pore radius were obtained from the Barret–Joyner–Halenda (BJH) model. Between 400 and 1000 cm⁻¹, FTIR spectra were acquired using an FTIR spectrometer (Nicolet iS10, Thermo Fisher Scientific, Waltham, MA, USA). The x-ray photoelectron spectrometer (Thermo Scientific K-alpha, XPS) equipped with a monochromatic X-ray source at 1486.6 eV, corresponding to the Al K_α line, was employed to collect XPS data (for the brevity XPS spectra not shown). The atomic weight percent of elements was determined using XPS.

The working electrode was prepared using nickel foam. Nickel foam was first cleaned with acetone and 3 M HCl solution for 10 min, which helped remove the oxide layer from the nickel foam’s surface. After washing and drying the nickel foam, the working electrode was made by combining 80 wt.% of the synthesized powder, 10 wt. % of acetylene black, and 10 wt. % of polyvinylidene difluoride (PVDF) with N-methyl pyrrolidinone (NMP). The slurry paste was obtained after thoroughly mixing the components and pasted onto the nickel foam. The prepared electrode on nickel foam was placed in a vacuum oven at a temperature of 60 °C for 10 h. The actual sample loading on nickel foam was determined by weighing the nickel foam before and after loading the nickel foam sample using an analytical balance (MS105DU, Mettler Toledo, max. 120 g, 0.01 mg of resolution). A platinum wire counter electrode, a saturated calomel electrode (SCE) as a reference electrode, and pasted samples on nickel foam as a working electrode were utilized in the electrochemical cell. The electrochemical performance of the electrode was assessed in 3 M KOH as an electrolyte. The electrochemical performance of the as-prepared electrode was evaluated using cyclic voltammetry (CV), galvanostatic charge-discharge techniques, and electrochemical impedance spectroscopy (EIS) by using a Versastat4-500 electrochemical workstation (Princeton Applied Research, Oak Ridge, TN, USA).

2.3. DFT Study

First-principles density functional theory (DFT) calculations were performed using the Perdew–Burke–Ernzerhof (PBE) [33] exchange-correlation functional for self-consistent calculations and geometry optimization along with the advanced hybrid functional HSE06 method to find the electronic density of states of the Al_xNi_1−xCo₂O₄ structure. The variable x ranges from 0 to 1 with a step of 0.25. Calculations were performed using VASP [34] (Vienna ab initio simulation package) with the projected augmented wave (PAW) [35] type pseudo-potential and the plane-wave basis set. The exchange and correlation part for the HSE calculations was described by hybrid functional [36] containing 40% Hartree–Fock exchange. We used a 4 × 4 × 2 k-points mesh centered at the Gamma point with a plane wave cutoff of 400 eV, which is enough for the self-consistent calculation. The global break condition for the electronic and ionic self-consistent (SC) loop was set to be 10⁻⁴ eV and 10⁻³ eV, respectively.

3. Results and Discussion

The room temperature powder XRD patterns of Ni_1−xAl_xCo₂O₄ samples are shown in Figure 1. The XRD patterns of the as-synthesized doped Ni_1−xAl_xCo₂O₄ are almost identical and can be assigned to NiCo₂O₄ (ICDD card no.02-0770). The diffraction peaks are located at 18.9°, 31.3°, 36.8°, 38.4°, 44.5°, 55.4°, 59.4°, and 65.2° can be assigned to the (111), (220), (311), (222), (400), (422), (511), and (440) agreeing with the standard crystallographic pattern of spinel NiCo₂O₄ phase respectively [37]. The absence of impurities in the XRD patterns confirms the complete dissolution of Al³⁺ into the Ni_1−xAl_xCo₂O₄ compound. The lattice parameter of Ni_1−xAl_xCo₂O₄ is listed in

Table 2. Surface area parameters derived from the adsorption-desorption curve, crystallite size obtained using Scherrer’s formula, lattice parameter obtained from XRD refinement, and Al³⁺ content was obtained from the XPS analysis of Ni_1−xAl_xCo₂O₄ compound.

Ni1−xAlxCo2O4, x	Al Content Determined via XPS	BET Surface Area (m2/g)	BJH Surface Area (m2/g)	BJH Avg. Pore Radius (nm)	BJH Avg. Pore Volume (cc/g)	Lattice Parameter, a (Å)	Crystallite Size (nm)
0.0	0.00	81.74	89.09	1.419	0.239	8.0983	13.212
0.2	0.37	83.16	126.13	1.423	0.183	8.0804	16.525
0.4	0.57	107.19	120.26	1.522	0.192	8.0666	13.804
0.6	0.79	106.29	128.81	2.144	0.194	8.0585	14.776
0.8	0.89	189.00	354.83	1.593	0.276	8.0449	12.286

. As expected, the lattice parameter of Ni_1−xAl_xCo₂O₄ decreases with the content x due to Al³⁺ (r_ion ~0. 535 nm) replacing a bigger Ni ion (r_ionic ~0.6 nm). The average crystallite size of the NiCo₂O₄ calculated using Scherrer’s formula [38] is listed in Table 2. The average crystallite size of Ni_1−xAl_xCo₂O₄ falls within the range of 12.28 nm for x = 0.8 to 16.52 nm for x = 0.2. Moreover, sharp diffraction peaks indicated the crystalline nature of the calcined Ni_1−xAl_xCo₂O₄ nanostructures. The aluminum content in Ni_1−xAl_xCo₂O₄ as determined via XPS (XPS study is not reported here) is listed in Table 2 and is within the accepted range of stoichiometry of the compound.

Figure 2 shows the SEM images of Ni_1−xAl_xCo₂O₄ as a function of Al³⁺ content. Figure 2 shows the gradual evolution of Ni_1−xAl_xCo₂O₄ nanostructures from an urchin-like architecture for x = 0.0 to nanoplates for x = 0.8. It is evident that the initial urchin whiskers slowly merge and morph into a plate-like structure with increasing Al content. The morphology of Ni_1−xAl_xCo₂O₄ is dependent on the Al³⁺ content in the sample. Usually, the morphology of hierarchal nanostructures depends on the synthesis techniques and urea concentration [4]. As a dispersing agent, urea helps avoid agglomeration by finely separating nanoparticles and plays a vital role in controlling the particles’ morphology [39,40,41].

Furthermore, NH₃ reacts with water to form NH₄⁺ and OH⁻, and CO₂ hydrolyzes to form Co₃²⁻ ions [42]. Tiny crystalline nuclei are created during the oxide crystal development phase, and nanoparticles of this oxide are precipitated by a rise in pH caused by NH₄⁺ ions generated from NH₃ owing to urea breakdown at higher temperatures [43,44]. Urea’s hydrolysis leads to a rise in the pH due to the increased release of NH₄⁺ in the solution [45]. The urea hydrolysis progresses slowly at milder circumstances and at a certain urea level, and the basic solution undergoes supersaturation of the metal-hydroxide species [46].

Figure 3a shows N₂ adsorption/desorption measurements were performed at 77 K between relative pressures of P/P_o ~ 0.029 to 0.99. From these curves, BET-specific surface areas and corresponding pore sizes were calculated by the BJH method for all Ni_1−xAl_xCo₂O₄ (0 ≤ x ≤ 0.8). Figure 2c inset shows the pore distribution. These curves show the largest number of pores distribution at around 3 nm for all Ni_1−xAl_xCo₂O₄. This pore distribution lies between ~2.0 and 6.0 nm, a range best suited to electrode materials [46]. The larger pore volume helps boost the diffusion kinetics inside the electrode material [40,46,47]. Table 2 lists the measured BJH pore radius, volume, and BET surface area of Ni_1−xAl_xCo₂O₄ samples. The high surface area of samples, in the range of 81.74–189 m²/g, indicates that the hydrated electrolyte ions could have an increased contact area at the electrolyte/electrode surface for the Faradaic redox reaction [48,49].

The FTIR spectrum Figure 3b shows two different bands arising from the stretching vibrations of the metal-oxygen bonds at 534.2 (ν₁) and 628.7 (ν₂) cm⁻¹ [1]. For the crystal structure of NiCo₂O₄, the ν₁ band is attributed to M-O (M = Ni) bond vibrations in octahedral coordination, whereas the ν₂ band is attributed to M-O (M-Co) bond vibrations in tetrahedral coordination. The stretching vibration of Co³⁺-O²⁻ in the octahedral sites is shown by the FTIR spectrum at 534.2 cm⁻¹, and the vibration of Ni²⁺-O⁻ at the tetrahedral sites is indicated by the FTIR spectrum at 628.7 cm⁻¹ for NiCo₂O₄, respectively [50]. The presence of these vibration bands indicates that pure phase spinel NiCo₂O₄ nanostructures have formed.

Table 3. Two distinct band positions from FTIR spectrum for Ni_1−xAl_xCo₂O₄ samples.

Ni1−xAlxCo2O4, x	ν1 (cm−1)	ν2 (cm−1)
0.0	534.2	628.7
0.2	544.0	645.2
0.4	551.0	652.7
0.6	553.6	653.9
0.8	556.1	665.1

lists the band positions for the tetrahedral and octahedral sites of the compound. Aluminum doping brings the shift of the stretching peak towards the right with an increase in Al³⁺ content. As the vibrational frequency is inversely proportional to the mass, the increase in vibrational frequency with Al³⁺ is expected due to the lower atomic weight of Al³⁺ replacing the Ni ion.

Cyclic voltammetry (CV) and galvanostatic charge-discharge (GCD) measurements give the electrochemical performance of Ni_1−xAl_xCo₂O₄ nanostructures. Figure 4a–e shows the cyclic voltammograms obtained at different scan rates (2–300 mV/s) in a voltage window of 0–0.6 V (vs. SCE). The Faradaic reaction for NiCo₂O_4, x = 0 is given by the equation below [51,52].

MO + OH⁻ ↔ MOOH + e⁻

NiCo₂O₄ + OH⁻ + H₂O ↔ NiOOH + 2CoOOH + e⁻

In both the cathodic and anodic scans, there are noticeable oxidation and reduction peaks [53]. Obtained non-rectangular and unsymmetric CV curves indicate the pseudocapacitive nature of electrodes [54]. The CV curves of Ni_1−xAl_xCo₂O₄ show no extra redox peaks, implying that Ni_1−xAl_xCo₂O₄ redox processes are very similar to NiCo₂O₄ redox processes. The peaks are due to the redox reaction related to M-O/M-O-OH, where M stands for Ni or Co ions [55]. A positive shift in oxidation peak potential and a negative shift in reduction peak potential are found as scan rate increases [56]. This PC characteristic is derived from the Faradaic redox reaction related to the reversible reaction of Ni²⁺/Ni³⁺ and Co³⁺/Co⁴⁺ transitions. Figure 5a–e shows the cyclic stability curves of Ni_1−xAl_xCo₂O₄ measured up to 2000 cycles without any change in the position of redox peaks.

The Randles–Sevcik plots of the Ni_1−xAl_xCo₂O₄ samples are shown in Figure 5f. When the scan rate is increased from 2 to 300 mV s⁻¹, the anodic peak current (I_pa) and the cathodic peak current (I_pc) both increase. More precisely, the Randles–Sevcik equation shows that both I_pa and I_pc change linearly with the square root of the scan rate (ν^1/2) [57,58]. These results suggest that the electrode–electrolyte interface’s redox events are rapid, quasi-reversible, and only limited by electrolyte diffusion [59]. The positive shift of oxidation peak and negative shift of reduction peak potential indicate relatively low electrode materials resistance and good electrochemical reversibility [1].

The calculation of specific capacitance was done from CV curves using the equation below [60]:

$$C_{\mathit{sp}}=\frac{{{\displaystyle \int }}_{V_1}^{V_2}I\ast V\ast dV}{m\ast v\ast \left(V_2-V_1\right)}$$

(1)

where I (A) is current, m (g) is the mass of the active material, v (mV/s) is the scan rate, V₁ and V₂ (V) are the start and end voltage of the CV scan.

Figure 4f shows the C_sp of Ni_1−xAl_xCo₂O₄ samples as a function of scan rates. The obtained values of C_sp for Ni_1−xAl_xCo₂O₄ nanostructures are 512, 368, 371, 380, and 356 F/g at a scan rate of 2 mV/s for Ni_1−xAl_xCo₂O₄ x = 0.0, 0.2, 0.4, 0.6, and 0.8, respectively. Compared to all Ni_1−xAl_xCo₂O₄ nanostructures, NiCo₂O₄ x = 0.0 exhibited the highest specific capacitance at all scan rates with a maximum value of 512 F/g at 2 mV/s. The higher electrochemical utilization and high electroactive surface area of the synthesized nanostructure occur for Ni_1−xAl_xCo₂O₄, x = 0. The observed decrease in C_sp with Al³⁺ substitution suggests that in the presence of Al³⁺, the Ni and Co in Ni_1−xAl_xCo₂O₄ may not behave synergistically as an effective Faradaic charging site as compared with the Ni and Co sites in NiCo₂O₄. Moreover, it can be deduced that, with the increase in Al³⁺ content, overall effective redox sites have decreased, which leads to a decrease in the specific capacitance.

The total charge storage process is affected by factors such as the Faradaic contribution from the insertion process of electrolyte ions, the Faradaic contribution from the charge-transfer process, and due to the high surface area, the contribution from both pseudocapacitance and the non-Faradaic contribution from the double layer effects [61]. CV data are analyzed at different scan rates and can characterize capacitive effects by equation [62,63].

I = av^b

(2)

where the peak current is I (A), the scan rate is v (mV/s), and a and b are adjustable parameters. The value of b specifies the charge storage mechanism. If b = 1, it indicates a surface-controlled reaction that the capacitive surface mechanism is dominant, and if b = 0.5, it indicates a diffusion-controlled Faradaic reaction during the charge storage mechanism [64]. Figure 5f shows the peak current vs. square root of the scan rate curve. These curves are fitted with the Equation (2) which gives values of b ~ 0.6239, 0.6398, 0.6218, 0.6134, and 0.5894 for x = 0.0, 0.2, 0.4, 0.6, and 0.8, respectively. That b > 0.5 for all electrodes indicates both diffusion-controlled and surface-controlled reactions [54]. Moreover, in Figure 5f, the lower value of the scan rate region has a linear dependence relation with the square root of the scan rate; this indicates a quasi-reversible electrochemical reaction [65]. Figure 5f also indicates that, at a higher scan rate > 100 mV/s, the ion diffusion is limited to the surface of the active material of the electrode, i.e., EDLC dominates the pseudocapacitor, and diffusion of OH⁻ ions can adhere only to the outer layer of the nanostructure, which contributes less to the electrochemical capacitive behavior [66]. On the other hand, the Faradaic redox reaction dominates at scan rates < 100 mV/s due to more effective usage of the working electrode’s active material [67]. Moreover, the diffusion of OH^- ions can easily penetrate deep into the nanostructure’s interlayer, which leads to adsorption of more ions and hence ends up with higher specific capacitance [68].

The capacitive and diffusion processes are the two mechanisms for current response at a fixed potential [69,70]. The equation below gives their contributions [71].

C_sp = k₁ + k₂^−1/2

(3)

k₁ and k₂ are determined from the C_sp vs.V^−1/2 linear plot, where k₂ is the slope and k₁ is the intercept. k₁ indicates diffusion and k₂ capacitance contribution to the total specific capacitance for a given voltage. For the calculation, the specific capacitance is plotted against the slow scan rate up to a value of 20 mV/s and performed a regression fit using Equation (3). The obtained value for k₁ and k₂ was used to determine the fractional contribution in terms of diffusion and capacitance from total specific capacitance [72], as given in Figure 5g. Figure 5g shows that the contribution to the current response at a fixed potential is more diffusive than capacitive and decreases with the Al content.

To further estimate the potential application of as-prepared electrodes for supercapacitors, the galvanostatic charge-discharge (GCD) measurement was performed. GCD was measured within the voltage window of 0.0 to 0.6 at different current densities between 1 A/g to 30 A/g in 3 M KOH solution are given in Figure 6a–e. These figures indicate that prepared electrodes of Ni_1−xAl_xCo₂O₄ can charge and discharge rapidly with good electrochemical reversibility at different constant current densities. As in CV curves, the non-linear relationship between the potential and time in both the charge-discharge cycle indicates the capacitance of studied nanomaterials is not constant in between the potential range and reflects typical PCs behavior [73]. The non-linearity of the GCD curve is a consequence of the Co³⁺/Co⁴⁺, and Ni²⁺/Ni³⁺ ions redox reactions with OH⁻. The discharge process occurs in three distinct steps. First, there is a rapid potential drop that occurs owing to the internal resistance. The second is a slow potential decay at an intermediate time, which is due to the Faradaic redox reaction. The third is a fast potential decay which is created by electric double layer capacitance [74]. The specific capacitance values for electrodes of Ni_1−xAl_xCo₂O₄ were calculated by using Equation (4) [75].

$$C_{\mathit{sp}}=\frac{I\ast t}{m\ast \Delta V}$$

(4)

where C_sp is the specific capacitance (F/g), I (A) is the charge-discharge current, ΔV (V) is the potential range, m (g) is the mass of the electroactive materials and t (s) is the discharging time.

The pseudocapacitance behavior of electrodes with respect to their discharging time is confirmed by the presence of voltage plateau on GCD curves [76]. The specific capacitance values of Ni_1−xAl_xCo₂O₄ were calculated using Equation (4) at 0.5 A/g are 268, 194, 192, 176, and 167 F/g for x = 0.0, 0.2, 0.4, 0.6, and 0.8, respectively, and are listed in

Table 4. Data of specific capacitance obtained from cyclic voltammetry, specific capacitance obtained from charge-discharge curves, energy density, and power density for Ni_1−xAl_xCo₂O₄.

x, Content	Specific Capacitance at 2 mv/s (F/g)	Specific Capacitance at 0.5 A/g (F/g)	Energy Density (Wh/Kg)	Power Density (W/Kg)
0.0	512	268	12.4	6316.6
0.2	368	194	8.9	6289.1
0.4	371	192	9.1	5281.2
0.6	380	176	8.7	6670.9
0.8	356	167	8.2	4660.4

. Figure 6f shows the dependence of current density with specific capacitance. A decrease in the specific capacitance with the increase of the discharge current is caused by the insufficient time available for the diffusion of the electrolyte ions into the inner electrode surface and the increase of the potential drop towards higher discharge currents [77].

The cyclic stability of Ni_1−xAl_xCo₂O₄, x = 0.0, and x = 0.8, electrodes were evaluated by the repeated charge-discharge measurements up to 5000 cycles at a constant current density of 10 A/g and 5 A/g in the potential range between 0.0 and 0.6 V in 3 M KOH, shown in Figure 7a,b respectively. The percentage retention in specific capacitance was calculated using Equation (5) [78].

% retention in specific capacitance = (C_#/C₁) × 100

(5)

where C_# and C₁ are specific capacitance at various cycles and the 1st cycle, respectively. Hence, the % retention in specific capacitance after 5000 cycles for Ni_1−xAl_xCo₂O₄ x = 0.0 and x = 0.8 is 97% and 108.5% respectively. This shows that even after 5000 cycles, the electrode shows outstanding cyclic stability, and hence there is not much high specific capacitance fade from its initial values in x = 0.0. Similarly, in the case of x = 0.8, the specific capacitance improves to 108.5% of its initial value, and this suggests an improvement in the surface moistened by electrolyte ions, access of electrolyte ions up to the pore of electrode materials, and activation of electrodes by electrolyte ions upon cycle repetition [79]. Furthermore, the Coulombic efficiency for x = 0.0 and x = 0.8 is 100% and 99.2%, as shown in Figure 7a,b, respectively.

Ragone plots of synthesized Ni_1−xAl_xCo₂O₄ are given in Figure 8a. The energy (E) and power (P) densities are determined using the equations below [32],

E = (1/2)CV²

(6)

P = E/t

(7)

where C is the electrodes’ specific capacitance, V is the operating voltage, and the letter t denotes the discharge time in seconds. An excellent-performance supercapacitor must exhibit a high value of both energy density and power density together. Figure 8a shows the superior performance of electrodes with an energy density of up to 12.4 Wh/kg for x = 0.0, with up to a peak power density of 6670.9 W/kg for x = 0.6. Furthermore,

Table 5. Comparison of specific capacitance of electrodes derived from doping of transition elements on NiCo₂O₄ using the hydrothermal method.

Electrode Material	Scan Rate (mV/s)	Specific Cap. (Csp)	Current Density (A/g)	Specific Cap. (Csp)	Reference
Ni1−xRuxCo2O4 (x = 0.00, 0.02, 0.05, 0.10, 0.20)	----	-----	1	831 F/g, 913 F/g, 1004 F/g, 1527 F/g, 1483 F/g	[27]
Ni1−xMnxCo2O4 (x =1.0, 0.7, 0.5, 0.3, 0.0)	-----	----	1	243 F/g (x = 0.3), 366 F/g (x = 0.5, max), 162 F/g (x = 0.7)	[28]
Ni1−xZnxCo2O4 (x = 0.5)	---	----	2	1228 F/g	[30]
Ni1−xFexCo2O4 (x = 1)	----	----	1	2237 F/g at 1 h,2208 F/g at 3 h, 1976 F/g at 5 h.	[31]
Ni1−xCaxCo2O4 (x = 0.0, 0.2, 0.4, 0.6, 0.8)	2	506 F/g,508 F/g, 355 F/g, 934 F/g, 571 F/g	0.5	284 F/g, 365 F/g, 317 F/g, 558 F/g, 253 F/g	[32]
Ni1−xAlxCo2O4 (x = 0.0, 0.2, 0.4, 0.6, 0.8)	2	512 F/g,368 F/g, 371 F/g, 380 F/g, 356 F/g	0.5	268 F/g, 194 F/g, 192 F/g, 176 F/g, 167 F/g	[This study]

compares the electrodes’ specific capacitances derived from the doping of transition elements on NiCo₂O₄ using the hydrothermal method.

The superior electrochemical performance of as-synthesized Ni_1−xAl_xCo₂O₄ nanostructure electrode materials was further confirmed by the electrochemical impedance spectroscopy (EIS) measurement held at 10 mV AC perturbation in the frequency range from 10 kHz to 0.05 Hz. Figure 8b shows the tendency of the real part of the impedance (Z′) to decrease with increasing frequency for all samples. When x = 0.0, it has the lowest real part of the impedance. The magnitude of the real part of the impedance (Z’) increases as Al³⁺ content increases, which affects the rise of grains, grain boundaries, and electrode interface resistance. Figure 8c shows the Nyquist plots for all of the Ni_1−xAl_xCo₂O₄ nanostructures. In this plot, the NiCo₂O₄ nanostructure electrode exhibits a small semicircle that suggests low internal resistance and charge transfer resistance, and hence this resistance is smaller than Al³⁺ doped NiCo₂O₄. Towards higher frequencies, the real part (Z_real) of the impedance represents a combined resistance of contact resistance at active material/current collector interface, ionic resistance of the electrolyte, and intrinsic resistance of the nanomaterials used in the electrode [80]. The semicircle at the high-frequency range is due to the Faradaic charge transfer resistance (R_ct) of the redox reaction between the electrode and electrolyte. The decrease in the diameter of the semicircle indicates lower charge transfer resistance [81]. These variations of the EIS curve from semicircle towards higher frequencies and linearity towards lower frequencies relate to reversible Faradaic redox reaction and access of OH⁻ ions into the pore of the electrode material. This behavior of the EIS curve describes the PCs properties of the prepared Ni_1−xAl_xCo₂O₄ nanostructure electrode [82]. The slope of the 45° portion of the curve is called Warburg resistance (Z_w) and results from the frequency dependence of ion diffusion/transport in the electrolyte to the electrode surface [83,84]. The Nyquist plot is a vertical line for an ideal electrode material, and more is the vertical line represents better electrolyte diffusion and capacitive behavior [85]. Figure 8c shows that R_ct is smaller for NiCo₂O₄ and has a more vertical line, which offers better rate capability.

First-principles density functional theory (DFT) calculations were performed using VASP² [34] (Vienna ab initio simulation package) with the projected augmented wave (PAW) [35] type pseudo-potential and the plane-wave basis set. The exchange and correlation part for the HSE calculations was described by hybrid functional [36] containing a 40% Hartree-Fock exchange. Here, we used a 4 × 4 × 2 k-points mesh centered at the Gamma point with a plane wave cutoff of 400 eV, which is enough for the self-consistent calculation. The global break condition for the electronic and ionic self-consistent (SC) loop was set to be 10⁻⁴ eV and 10⁻³ eV, respectively. In our calculation, we dealt with a cell containing 28 atoms as a total. The value of gaussian smearing was adjusted to be 0.20. Even though the ratio of x is 20% in the experiment, we chose 25% in the calculations because changing it to 20% would require a much larger supercell, which eventually increases the computing time and memory excessively. Figure 9a shows the side view along the x-axis of the crystal structure of Ni_1−xAl_xCo₂O₄, while Figure 9b shows the band-gap variation from 2.42 eV to 4.82 eV for the different values of x. For the system without Al atom, which is Ni₁Co₂O_4, the band gap was found at 2.42 eV, which agrees with the previously reported experimental values of 2.06 eV [86] and 2.64 eV [87]. This value is the minimum among all the cases of adding Al to the system. Adding an Al atom above 0.8, the band-gap reached its maximum value of 4.82 eV for the case where Al completely replaces Ni, as shown in Figure 9a. These results suggest that the energy band-gap of Ni_1−xAl_xCo₂O₄ remains largely unaltered (x < 0.8) with the Al³⁺ substitution, which translates to the fact that the electrical conductivity of Ni_1−xAl_xCo₂O₄ could remain unaltered with the Al³⁺ substitution below x = 0.8. In the absence of the contribution of electrical conductivity to the electrochemical performance of Ni_1−xAl_xCo₂O_4, at least for x < 0.8, it can be concluded that the electrochemical performance of Ni_1−xAl_xCo₂O₄ could be dictated by other factors, such as morphology, the oxidation potential of participating ions, active metal sites, and the internal resistance of the compound.

4. Conclusions

Finally, a simple and cost-effective technique was developed to report the effect of Al³⁺ doping on the structural and electrochemical performance of the NiCo₂O₄ nanostructure. The experimental results revealed the dependence of surface area, specific capacitance, and electrochemical performance of NiCo₂O₄ on Al³⁺ content. The doped compound’s overall electrochemical performance, viz. energy and power density, was observed to be significantly superior to other doped NiCo₂O₄ compounds. The DFT calculation suggests that the energy band-gap of Ni_1−xAl_xCo₂O₄ (<0.8) remains largely invariant, and hence the electrical conductivity. Thus, the observed reduction in the specific capacitance of doped Ni_1−xAl_xCo₂O₄ could be attributed to the decrease in active redox sites with Al³⁺ substitution. The study indicates that Al³⁺ can be an excellent cost-effective substitution for expensive cobalt in NiCo₂O₄ to provide the desired electrochemical performance. Hence, this work highlights the development of cheaper and more promising electrode doped materials without sacrificing the performance of NiCo₂O₄.