Templating Nanostructured Aromatic Based Materials as Possible Anode Electrodes for Na-ion Batteries: A Computational DFT Approach

Templating Nanostructured Aromatic Based Materials as Possible Anode Electrodes for Na-ion Batteries: A Computational DFT Approach

A. Ahmad A.A.M. Abahussain M.H. Nazir S.Z.J. Zaidi*

Institute of Chemical Engineering and Technology, University of the Punjab, Pakistan

Department of Chemical Engineering, College of Engineering, King Saud University, Riyadh 11421, Saudi Arabia

Faculty of Computing Engineering and Sciences, University of South Wales, United Kingdom

Laboratory for Energy Water and Healthcare Technologies, University of the Punjab, Pakistan

Corresponding Author Email: 
zohaib.icet@pu.edu.pk
Page: 
38-43
|
DOI: 
https://doi.org/10.14447/jnmes.v27i1.a06
Received: 
2 June 2023
|
Revised: 
10 October 2023
|
Accepted: 
17 October 2023
|
Available online: 
31 January 2024
| Citation

© 2024 The authors. This article is published by IIETA and is licensed under the CC BY 4.0 license (http://creativecommons.org/licenses/by/4.0/).

OPEN ACCESS

Abstract: 

Due to their widespread availability and lower costs Na-ion batteries have grabbed attention as a promising substitute for lithium-ion batteries. In this paper, we delve into a computational investigation of the potential use of 2,4-Chloronitrotoluene and pyrazine as electrode materials for SIBs. We utilized Density Functional Theory (DFT) calculations to examine their properties and suitability for energy storage applications. Adsorption energies, highest occupied molecular orbital, lowest unoccupied molecular orbital, HLG energies, total energies, and bond lengths were computed using computational theories for aromatic compounds. This computation is based here on the DFT and computational engineering modelling approach. The results indicate that both 2,4-chloronitrotoluene and pyrazine exhibit favorable adsorption energies for Na-ion adsorption, indicating their potential as electrode materials. Furthermore, the investigation delved into the electronic properties of the materials. Notably, the analysis showed that the EHOMO, ELUMO, and EHLG are facilitative to the efficient transfer of electrons during the operation of the battery. Moreover, the computed bond lengths indicate that stable Na-ion adsorption can appear on both materials, further solidifying their potential for integration into the framework of Na-ion batteries. In essence, this study yields valuable insights into the electrochemical aspects of 2,4- chloronitrotoluene and pyrazine, highlighting their promising candidacy as electrode materials for SIBs.

Keywords: 

2,4- Chloronitrotoluene, pyrazine, Na-ion batteries, DFT, energy storage, electrode materials

1. Introduction

Na-ion batteries have garnered significant attention as a favorable substitute to the prevalent Li-ion batteries, which presently dominate the storing energy landscape. Sodium, being more abundant and cost-effective than lithium, presents a sustainable choice for storing energy on large scale, marking a important step towards a more sustainable environment [1]. Sodium-ion batteries present hope for a greener future. By using sodium ions rather of the more common lithium ions, these batteries are not only cheaper but also better for the environment [2].

A key advantage of sodium-ion batteries lies in their utilization of abundant materials. With sodium ranking as the sixth most abundant element on our planet, its general availability exceeds that of lithium. This intrinsic abundance implies that scaling up the production of sodium-ion batteries could be achieved more seamlessly, without facing the geopolitical challenges associated with lithium-ion battery production [3]. Moreover, the materials used in sodium-ion batteries are generally safer and less harmful than those used in lithium-ion batteries. LIBs often lean on resources like cobalt, which are frequently sourced under environmentally damaging and ethically concerning conditions [4]. Contrariwise, sodium-ion batteries can utilize materials like $\mathrm{Mg}$, which are more abundant and environmentally sustainable. This distinction positions sodium-ion batteries as potentially more cost-effective than their lithium-ion counterparts [5]. With the surge in battery demand, the limitations and increasing prices of lithium are becoming evident. In contrast, the vast availability and affordability of sodium present a more sustainable avenue for large-scale battery production [6]. Furthermore, sodium-ion batteries have the potential to offer higher energy densities and more extended lifetimes when compared to LIBs. However, making high-performance SIBs comes with its own set of challenges, such as dealing with the relatively larger size of $\mathrm{Na}$ ions and the absence of appropriate electrode materials [7].

Sodium-ion batteries (SIBs) function based on the same basic principle as other rechargeable batteries, involving the movement of ions between two electrodes through an electrolyte during charging and discharging cycles. In the case of SIBs, the charge carriers are sodium ions $(\mathrm{Na}+)$, in contrast to lithium ions used in lithium-ion batteries [8]. During the charging process, sodium ions are drawn from the cathode material, generally made of layered transition metal oxides or polyanionic compounds. These ions cross through the electrolyte, usually a salt dissolved in an organic solvent, and get stored in the anode material, generally composed of carbon. Meanwhile, electrons transfer through the external circuit to the anode, creating a voltage difference between the two electrodes [9]. When the battery discharges, the process inverts. The stored $\mathrm{Na}$ ions in the anode travel back to the cathode through the electrolyte, while electrons float in external circuit, generating an electric current to power devices. This cycle can be replicated multiple times, enabling the battery to be recharged and used again [10].

The performance SIBs relies on various factors, including the choice of electrode and electrolyte materials, the structure of the cell, and the operating conditions. One of the primary challenges in SIB development is the quest for materials that can present high capacity, stable cycling, and fast charge/discharge rates while being cost-effective and abundant. Some studies have also concentrated on transition metal oxides for various electrochemical aspects of energy generation and storage [11]. Continued research in this domain seeks to enhance the efficiency and reliability of SIBs for diverse energy storage applications. 2,4-Chloronitrotoluenes and pyrazine have appeared as promising anode materials for $\mathrm{Na}-\mathrm{ion}$ Batteries (SIBs) due to their strong conductivity and capacity. 2,4-Chloronitrotoluenes, a type of carbon with a special cylindrical structure, offer a considerable surface area and superior electrical conductivity [12]. On the other hand, pyrazine is an organic molecule with a flat structure that can host sodium ions between its layers [13]. This study employed density functional theory (DFT) to examine 2,4chloronitrotoluenes and pyrazine, providing valuable insights into their electrochemical properties as anode materials for SIBs. These understandings have played a pivotal role in driving the design and enhancement of these materials to boost battery performance. Key properties such as adsorption energies, HOMO energy, LUMO energy, and the energy difference between HOMO and LUMO orbitals were closely examined.

Figure 1. New sodium-ion battery materials: 2,4-Chloronitrotoluene and pyrazine

2.Computational Details

Our research methodology relied on an ab initio approach, using DFT calculations through computational engineering modeling package along with a molecular editor and visualizer [14, 15]. The main aim of our research was to gain insights into the fundamental electronic properties of electrode materials. We specifically concentrated on analyzing a single atomic layer of these materials. To achieve this, we examined adsorption energies, EHOMO, ELUMO, and EHLG, considering the geometrically optimized structures using density functional theory (DFT). We adjusted the vibrational frequencies to guarantee that the calculated geometries accurately represented the local minimum points [16]. The structures of 2,4-chloronitrotoluenes and Pyrazine with Sodium ion were constructed using the DFT/B3LYP approach. For various tasks such as structure optimization and energy calculations, the method selected was ground state, utilizing DFT with default spin, and the Basis was selected as $6-31 \mathrm{G}$.

3. Results and Discussions

3.1 Geometrically optimized structures of 2,4-Chloronitrotoluenes and pyrazine with $\mathrm{Na}$

This study involved an analysis of the interaction between sodium ions and 2,4-Chloronitrotoluenes and Pyrazine molecules. The investigation entailed an examination of four distinct optimized structures for each compound. Optimized structures were generated by systematically altering the positions of sodium ions within the molecular framework. The optimization was conducted using the 'optimization' job type, and DFT served as the ground-state method. The calculations were performed with default spin settings and a 6-31G basis set, providing insights into the most energetically favorable configurations of sodium ions within the molecules. The results obtained from this study have meaningful implications for the potential use of these molecules as electrode materials in SIBs.

In Figure 2, we showcase a visual representation of our optimized structures of 2,4-chloronitrotoluene with sodium ion configurations, providing insight into their spatial arrangement and geometry.

Figure 2. Optimized structures of 2,4-Chloronitrotoluenes and pyrazine with sodium ions

3.2 Theoretical capacity and energy density

The Theoretical capacity ($Q_{\text {theoretical }}$) of the battery can be calculated using Faraday's law:

$\mathrm{Q}_{\text {theoretical }}=\mathrm{C}=\mathrm{nF} / \mathrm{M}$              (1)

The theoretical capacity, represented as $Q_{\text {theoretical, is }}$ measured in $\mathrm{Ah} / \mathrm{g}$ or $\mathrm{Ah} / \mathrm{cm}^3$. Here, 'n' represents the number of electrons participating in the redox reaction, ' $F$ ' is the Faraday constant $(96,485 \mathrm{C} / \mathrm{mol})$, and ' $\mathbf{M}$ ' represents the molar mass of the active material. Pyrazine, having a molecular weight of $80.10 \mathrm{~g} / \mathrm{mol}$, is capable of engaging in a sodium-ion battery configuration with a distinct charge transfer reaction.

$\mathrm{C}_4 \mathrm{H}_4 \mathrm{~N}_2+\mathrm{Na}^{+}+\mathrm{e}^{-} \rightarrow \mathrm{NaC}_4 \mathrm{H}_4 \mathrm{~N}_2$                  (2)

In this reaction, one electron is transferred per sodium ion, so $\mathrm{n}=1$. Plugging in the values, we get $\mathrm{Q}=1204 \mathrm{Ah} / \mathrm{g}$. The charge transfer mechanism for $2,4-\mathrm{CNT}$ in a sodium ion battery is:

$2,4-\mathrm{CN}+\mathrm{Na}^{+}+e^{-} \rightarrow \text { Sodium 2,4-CN }$                     (3)

In this reaction, one electron is transferred per sodium ion, so $\mathrm{n}=1$. The molar mass of $2,4-\mathrm{CNT}$ is $181.59 \mathrm{~g} / \mathrm{mol}$. Plugging in the values, we get: $\mathrm{Q}=(1\times 96,485 \mathrm{C} / \mathrm{mol}) /$ $181.59 \mathrm{~g} / \mathrm{mol}=531 \mathrm{Ah} / \mathrm{g}$.

In order to convert the units of theoretical Capacity from $\mathrm{Ah} / \mathrm{g}$ to $\mathrm{Wh} / \mathrm{kg}$ we need to multiply it with Vcell and divide it by 1000 . In order to calculate Vcell we use $\mathrm{V}_{\text {cell }}=\Delta \mathrm{G} / \mathrm{nF}$, Where $\mathrm{F}$ is Faradays constant $\mathrm{n}$ is charge on sodium Delta $\mathrm{G}$ is gibbs free energy. We used B3LYP 6-31+G basis for calculations. The values of $\mathrm{V}_{\text {cell }}$ of 2,4-Chloronitrotoluenes and pyrazine were calculated as 5.621 volts and 1.778 volts respectively and the energy densities are calculated as 2.991 $\mathrm{Wh} / \mathrm{kg}$ and $2.140 \mathrm{Wh} / \mathrm{kg}$ respectively.

3.3 Sodium Ion adsorption energies

The adsorption energy of sodium in 2,4-Chloronitrotoluenes and pyrazine plays a crucial role in determining capacity of Na-ion intercalation, stability and durability of the anode material [18]. A high adsorption energy of sodium in 2,4-Chloronitrotoluenes and pyrazine can result in strong binding between the sodium ions and the electrode, which can improve rate capability and cycling stability of the electrode material. However, if the adsorption energy is too high, the sodium ions may become strongly trapped in the electrode material, which can reduce the capacity and voltage of the battery. [18] In our study we calculated the adsorption energies of sodium on different transition structures of 2,4-Chloronitrotoluene and Pyrazine by selecting the job as frequency and the method as ground state, DFT, default spin and the Basis as $6-31 \mathrm{G}$ in computational engineering modeling package. The adsorption energies can be calculated as [19].

$\mathbf{E}_{\text {ad }}=E_{\text {total }}(o p t)-\left(E_{\text {mon } 1}(o p t)+E_{\text {mon } 2}(o p t)\right)$         (4)

Table 1. Adsorption energies of Na-ion on 2,4-CN

Transition State

Ead(eV)

PR1

16.526

PR2

16.526

PR3

16.526

PR4

16.526

Table 2. Adsorption energies of Na-ion on Pyrazine

Tansition State

Bond Length (Å)

CN1

2.170

CN2

2.170

CN3

2.170

CN4

2.170

It can be seen from the above result that the adsorption energies of sodium in 2,4- Chloronitrotoluene varies between 21.621 eV to 30.068 eV the higher Ead value for CN4 is due to the presence of a new lone pair on the 2,4-Chloronitrotoluene molecule. When a sodium ion adsorbs onto the 2,4-Chloronitrotoluene molecule, it can interact with the electron density in the molecule, including any lone pairs. The adsorption energy of sodium in Pyrazine was measured as 16.5266 eV. With respect to adsorption energy pyrazine is favorable electrode material because its adsorption energy is not too high like 2,4-Chloronitrotoluene. Sodium ions may become strongly trapped in the 2,4- Chloronitrotoluene, which can reduce the capacity and voltage of the battery.

3.4 Bond lengths

The bond length between $\mathrm{Na}$ with 2,4-Chloronitrotoluene and pyrazine can have a significant effect on the working of a sodium-ion battery. In general, the shorter the length of bond between $\mathrm{Na}$ and anode material, the stronger the connection between them, which can improve the battery performance [20].

This happens due to sodium $(\mathrm{Na})$ creating a strong connection with the anode material, it stores and carries $\mathrm{Na}$ ions within the material, boosting the battery's capacity and long-term stability. On the flip side, if the bond between them is too tight, it can restrict the movement of $\mathrm{Na}$ ions, leading to slow ion diffusion and ultimately reducing the battery's performance [21]. In our research we calculated the bond length of different transition states of $\mathrm{Na}$ with 2,4Chloronitrotoluene and pyrazine. The results are expressed below.

Table 3. Evaluated bond length of different transition states of sodium with 2,4-Chloronitrotoluene

Transition State

Bond Length Å

PR1

2.170

PR2

2.170

PR3

2.170

PR4

2.170

Table 4. Evaluated bond length of different transition states of sodium with pyrazine

Transition State

EHOMO(eV)

ELUMO(eV)

EHLG(eV)

ETotal (eV)

CN1

-5.596

-4.330

-1.266

-29866.473

CN2

-5.612

-4.368

-1.245

-29867.159

CN3

-5.583

-4.368

-1.215

-29859.199

CN4

-6.297

-4.543

-1.753

-29858.026

There was no difference observed between the bond lengths of $\mathrm{Na}$ with 2,4- chloronitrotoluene and pyrazine because in both cases the bond is between hydrogen and sodium.

3.5 Molecular orbital energies

The performance of SIBs can be significantly affected by the energy levels of the molecular orbitals in the anode materials. In the context of an SIB, the active electrode material undergoes a reversible process of inserting and removing $\mathrm{Na}$ ions during the charge and discharge cycles. This entire mechanism is dictated by the disparity in energy between the molecular orbitals of the electrode material and the Na ions [22]. More specifically, the energy level of the highest occupied molecular orbital (HOMO) of the electrode material affects the ease of $\mathrm{Na}$ ion insertion, while the energy level of the lowest unoccupied molecular orbital (LUMO) impacts the ease of $\mathrm{Na}$ ion removal [23]. If the energy level of the HOMO is too high, it can hinder the insertion of $\mathrm{Na}$ ions into the electrode material, leading to lower capacity and slower rate capability. Conversely, if the energy level of the LUMO is too low, it can make it difficult for the $\mathrm{Na}$ ions to de-insert, resulting in poor cycling stability [24].

In our study, we performed calculations to determine the energy levels of HOMO (Ehomo) and LUMO (Elumo), as well as the energy difference between them, for four different optimized structures of 2,4-Chloronitrotoluene and pyrazine. Figures 3 and 4 visually depict charge transfer processes in these molecules interacting with sodium ions, highlighting the roles of HOMO and LUMO orbitals in understanding electronic structure and energy levels. Additionally, we computed the total energy of the complex formed by the $\mathrm{Na}$ ions interacting with the electrode material. The total energy refers to the overall energy associated with the entire system, which includes both the electrode material and the $\mathrm{Na}$ ions. It represents the thermodynamic stability of the system and is a fundamental quantity in computational chemistry.

Table 5. Energy of HOMO, energy of LUMO, difference energy between HOMO and LUMO and total energy of Na on 2,4-Chloronitrotoluene

Transition State

Ead(eV)

CN1

21.621

CN2

20.935

CN3

20.935

CN4

30.068

Table 6. Energy of HOMO, energy of LUMO, difference energy between HOMO and LUMO and total energy of Na on pyrazine

Transition State

EHOMO (eV)

ELUMO (eV)

EHLG (eV)

ETotal (eV)

PR1

-8.360

-5.311

-3.049

-11589.365

PR2

-8.115

-5.356

-2.759

-11589.365

PR3

-8.361

-5.311

-3.049

-11589.36495

PR4

-8.361

-5.311

-3.049

-11589.365

Figure 3. HOMO and LUMO orbitals of 2,4-Chloronitrotoluene with sodium ions

Figure 4. HOMO and LUMO orbitals of pyrazine with sodium ions

As shown from the results that the $\mathbf{E}_{\text {HLG }}$ of 2,4-Chloronitrotoluene with sodium varies from $-1.215 \mathrm{eV}$ to $1.753 \mathrm{eV}$ and for pyrazine with sodium it varies from -2.759 $\mathrm{eV}$ to $-3.049 \mathrm{eV}$. So, 2,4-Chloronitrotoluene will be favorable because its $\mathbf{E}_{\text {HLG }}$ is not too high unlike pyrazine. In case of pyrazine when $\mathbf{E}_{\text {HLG }}$ is too high because of this the electrodes can lead to a higher risk of unwanted side reactions,  such as electrolyte decomposition, which can degrade the battery's performance over time.

Table 7. Specific capacitance, energy density and adsorption energy of anode material for sodium ion batteries

Anode Material

Specific Capacitance (F/g)

Enery Density (Wh/kg)

Ead(eV)

Reference

CNT loaded MnO2

140

21

-1.83

[25]

LiMn2O4

137

56

N.A

[26]

MnO2

45

43

3.77

[27]

Pyrazine

N.A

2.1

16.52

This Study

Activated Carbon

180

3.6

N.A

[27]

2,4-Chloronitro toluene

N.A

3

20.93

This Study

4. Conclusions

The research highlights the possibility of Na-ion batteries (SIBs) as a viable substitute to Li-ion batteries, citing the abundant and cost-effective nature of $\mathrm{Na}$. Investigating 2,4chloronitrotoluene and pyrazine as anode materials, the study emphasizes the influential role of sodium's adsorption energy, varying from $21.62 \mathrm{eV}$ to $30.07 \mathrm{eV}$ in 2,4-Chloronitrotoluene and $16.53 \mathrm{eV}$ in pyrazine. Additionally, insights into Еномо, ELUMO, and EHLG provide valuable guidance for material optimization, improving energy density and overall SIB performance. This research signifies a notable step forward in the development of efficient and reliable SIBs for diverse applications of energy storage, paving the way for further improvements in the field.

5. Acknowledgment

The authors also would like to extend their sincere appreciation to Researchers Supplorting Project number (RSPD2024R612), King Saud University, Riyadh, Saudi Arabia.

Nomenclature

CN

2,4-chloronitrotoluene

PR

Pyrazine

DFT

Density functional theory

SIBs

Sodium ion batteries

HOMO

Highest occupied molecular orbital

LUMO

Lowest unoccupied molecular orbital

HLG

HOMO LUMO Gap

B3LYP

Becke, 3-parameter, Lee–Yang–Parr

Subscripts

ad

adsorption

homo

highest occupied molecular

lumo

lowest unoccupied molecular orbital

  References

[1] Deng, J., Luo, W. B., Chou, S. L., Liu, H. K., Dou, S.X. (2018). Sodium ‐ ion batteries: from academic research to practical commercialization. Advanced Energy Materials, 8(4): 1701428. https://doi.org/10.3390/batteries9010048

[2] Zaidi, S.Z.J., Hassan, S., Raza, M., Harito, C., Yuliarto, B., Walsh, F.C. (2021). Conceptualized Simulation for Templating Carbon-Based Nano Structures for Li-ion Batteries: A DFT Investigation. Journal of New Materials for Electrochemical Systems, 24(2). https://doi.org/10.14447/jnmes.v24i2.a02

[3] Barpanda, P., Oyama, G., Nishimura, Si. et al. (2014). A 3.8-V earth-abundant sodium battery electrode. Nat Commun, 5: 4358. https://doi.org/10.1038/ncomms5358

[4] Zaidi, S.Z.J., Raza, M., Hassan, S., Harito, C., Walsh, F.C. (2021). A DFT Study of Heteroatom Doped-Pyrazine as an Anode in Sodium ion Batteries. Journal of New Materials for Electrochemical Systems, 24(1): 1-8. https://doi.org/10.14447/jnmes.v24i1.a01

[5] Alharbi, A.F., Abahussain, A.A.M., Nazir, M.H. and Zaidi, S.Z.J. (2022). A High-Energy-Density Magnesium-Air Battery with Nanostructured Polymeric Electrodes. Polymers, 14(15), p.3187. https://doi.org/10.3390/polym14153187

[6] Lukas Mauler, Fabian Duffner, Wolfgang G. Zeie, Jens Leker. (2021). Battery cost forecasting: a review of methods and results with an outlook to 2050. Energy Environ. Sci., 14: 4712-4739. https://doi.org/10.1039/D1EE01530C (Review Article)

[7] Zaidi, S.Z.J., Nazir, M.H., Raza, M. and Hassan, S., 2022. A High Energy Density Li-ion Battery with Lithium Titanium Oxide Anode. Int. J. Electrochem. Sci, 17(2). https://doi.org/10.20964/2022.04.07

[8] Wang, T., Su, D., Shanmukaraj, D. et al. (2018). Electrode Materials for Sodium-Ion Batteries: Considerations on Crystal Structures and Sodium Storage Mechanisms. Electrochem. Energ. Rev., 1: 200-237. https://doi.org/10.1007/s41918-018-0009-9

[9] Nurohmah, A.R., Nisa, S.S., Stulasti, K.N.R. et al. (2022). Sodium- ion battery from sea salt: a review. Mater Renew Sustain Energy, 11: 71-89. https://doi.org/10.1007/s40243-022-00208-1

[10] Jinsoo Park, Jin-Woo Park, Jeong-Hui Han, Sang-Won Lee, Ki-Young Lee, Ho-Suk Ryu, Ki-Won Kim Guoxiu Wang, Jou-Hyeon Ahn, Hyo-Jun Ahn. (2014). Charge–discharge properties of tin dioxide for sodium-ion battery. Materials Research Bulletin, 58: 186-189. https://doi.org/10.1016/j.materresbull.2014.04.051

[11] Abahussain, A.A., Alharbi, A.F., Nazir, M.H., Zaidi, S.Z.J. (2022). Anodic Generation of Hydrogen Peroxide by Using Manganese (III) Oxide Based Catalyst. Int. J. Electrochem. Sci, 17(221130): 2. https://doi.org/10.20964/2022.11.04

[12] Xiong, Z., Yun, Y.S., Jin, H.J. (2013). Applications of Carbon Nanotubes for Lithium Ion Battery Anodes. Materials (Basel), 6(3): 1138-1158. https://doi.org/10.3390/ma6031138

[13] Shi, R., Liu, L., Lu, Y. et al. (2020). Nitrogen-rich covalent organic frameworks with multiple carbonyls for high- performance sodium batteries. Nat Commun, 11: 178. https://doi.org/10.1038/s41467-019-13739-5

[14] Barca, G. M. J., Bertoni, C., Carrington, L., Datta, D., De Silva, N., Deustua, J. E., Fedorov, D. G., Gour, J. R., Gunina, A. O., Guidez, E., Harville, T., Irle, S., Ivanic, J., Kowalski, K., Leang, S. S., Li, H., Li, W., Lutz, J. J., Magoulas, I., Mato, J., Mironov, V., Nakata, H., Pham, B. Q., Piecuch, P., Poole, D., Pruitt, S. R., Rendell, A. P., Roskop, L. B., Ruedenberg, K., Sattasathuchana, T., Schmidt, M. W., Shen, J., Slipchenko, L., Sosonkina, M., Sundriyal, V., Tiwari, A., Galvez Vallejo, J. L., Westheimer, B., Wloch, M., Xu, P., Zahariev, F., & Gordon, M. S. (2020). GAMESS 2023 R1: Recent developments in the general atomic and molecular electronic structure system. The Journal of Chemical Templating nanostructured aromatic based materials as possible anode electrodes .../ J. New Mat. Electrochem. Systems Physics, 152(15): 154102. http://aip.scitation.org/doi/10.1063/5.0005188

[15] Avogadro: an open-source molecular builder and visualization tool. Version 1.2.0n. http://avogadro.cc/.

[16] João Otávio Anhaia-Machado, Artur Caminero Gomes Soares, Claudinéia Aparecida Sales de Oliveira Pinto, Andres Ignacio Ávila Barrera , André Rolim Baby, and Gustavo Henrique Goulart Trossini . Molecular Modeling Based on Time-Dependent Density Functional Theory (TD-DFT) Applied to the UV-Vis Spectra of Natural Compounds. https://doi.org/10.3390/chemistry5010004

[17] Ramirez-de-Arellano, J.M., Canales, M., Magaña, L.F. (2021). Carbon Nanostructures Doped with Transition Metals for Pollutant Gas Adsorption Systems. Molecules, 26: 5346. https://doi.org/10.3390/molecules26175346

[18] Abahussain, A.A., Zaidi, S.Z.J., Nazir, M.H., Raza, M., Hassan, S. (2022). A DFT study of Graphene as a drug carrier for gemcitabine anticancer drug. Journal of New Materials for Electrochemical Systems, 25(4). https://doi.org/10.14447/jnmes.v25i4.a02

[19] Al-Hamdani, Yasmine, Michaelides, Angelos, von Lilienfeld, Anatole. (2017). Exploring water adsorption on isoelectronically doped graphene using alchemical derivatives. The Journal of Chemical Physics, 147. https://doi.org/10.1063/1.4986314

[20] Lu, S., Wang, B., Zhang, P., Jiang, X., Zhao, X., Wang, L., Yin, Z., Wu, J. (2022). First-Principles-Based Optimized Design of Fluoride Electrolytes for Sodium-Ion Batteries. Molecules, 27(20): 6949. https://doi.org/10.3390/molecules27206949

[21] Qiannan Liu, Zhe Hu, Chao Zou, Huile Jin, Shun Wang, Lin Li.  Structural engineering of electrode materials to boost high-performance sodium-ion batteries. Cell Reports Physical Science. https://doi.org/10.1016/j.xcrp.2021.100551

[22] Hwang, J.Y., Myung, S.T., Sun, Y. K. (2017). Sodium-ion batteries: present and future. Chemical Society Reviews, 46(12): 3529-3614. https://doi.org/10.1039/C6CS00776G

[23] Miar, M., Shiroudi, A., Pourshamsian, K., Oliaey, A.R., Hatamjafari, F. (2021). Theoretical investigations on the HOMO– LUMO gap and global reactivity descriptor studies, natural bond orbital, and nucleus-independent chemical shifts analyses of 3- phenylbenzo[d]thiazole-2(3H)-imine and its para-substituted   derivatives: Solvent and substituent effects. Journal of Chemical Research. 2021;45(1-2):147-158. https://doi.org/10.1177/1747519820932091

[24] ACS Omega 2022, 7, 7, 5605–5614 Publication Date:February 11, 2022. https://doi.org/10.1021/acsomega.1c05794

[25] J. Yan, Z. Fan, T. Wei, W. Qian, M. Zhang and F. Wei, Carbon, 48 (2010) 3825-3833. https://doi.org/10.1002/adma.201001029

[26] J. Li, and I. Zhitomirsky, Materials Chemistry and Physics, 112 (2008) 525-530. https://doi.org/10.1016/j.matchemphys.2008.05.086

[27] S.Z.J Zaidi, S.Hassan, M.Raza, F.C Walsh. (2021). Transition metal oxide as possible electrode materials for Li-ion batteries: A DFT Analysis. Int. J. Electrochem. Sci., 16: 210322. https://doi.org/10.20964/2021.03.66

[28] J. Yan, Z. Fan, T. Wei, W. Qian, M. Zhang and F. Wei, Carbon, 48 (2010) 3825-3833.W. Xiao, H. Xia, J.Y.H. Fuh and L. Lu, Journal of The electrochemical Society, 156 (2009) A627.

[29] Kumar, U., Goonetilleke, D., Gaikwad, V., Pramudita, J.C., Joshi, R.K., Sharma, N., Sahajwalla, V. (2019). Activated Carbon from E- Waste Plastics as a Promising Anode for Sodium- Ion Batteries. ACS Sustainable Chem. Eng., 7(12): 10310-10322. https://doi.org/10.1021/acssuschemeng.9b00135