Adsorption of chemical warfare agents over C24 fullerene: Effects of decoration of cobalt

Adsorption of chemical warfare agents over C24 fullerene: Effects of decoration of cobalt

Accepted Manuscript Adsorption of chemical warfare agents over C24 fullerene: Effects of decoration of cobalt Alireza Soltani, Masoud Bezi Javan, Moha...

3MB Sizes 0 Downloads 19 Views

Accepted Manuscript Adsorption of chemical warfare agents over C24 fullerene: Effects of decoration of cobalt Alireza Soltani, Masoud Bezi Javan, Mohammad T. Baei, Zivar Azmoodeh PII:

S0925-8388(17)34132-4

DOI:

10.1016/j.jallcom.2017.11.350

Reference:

JALCOM 44043

To appear in:

Journal of Alloys and Compounds

Received Date: 16 May 2017 Revised Date:

16 October 2017

Accepted Date: 28 November 2017

Please cite this article as: A. Soltani, M.B. Javan, M.T. Baei, Z. Azmoodeh, Adsorption of chemical warfare agents over C24 fullerene: Effects of decoration of cobalt, Journal of Alloys and Compounds (2018), doi: 10.1016/j.jallcom.2017.11.350. This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.

ACCEPTED MANUSCRIPT Adsorption of chemical warfare agents over C24 fullerene: effects of decoration of Cobalt Alireza Soltania,*, Masoud Bezi Javanb, Mohammad T. Baeic, Zivar Azmoodehd Golestan Rheumatology Research Center, Golestan University of Medical Science, Gorgan, Iran b Physics Department, Faculty of Sciences, Golestan University, Gorgan, Iran c Department of Chemistry, Azadshahr Branch, Islamic Azad University, Azadshahr, Golestan, Iran d Young Researchers and Elite Club, Behshar Branch, Islamic Azad University, Sari, Iran

SC

RI PT

a

Abstract

M AN U

The present theoretical study suggests the adsorption of the chemical warfare agents (Soman, Chlorosoman, Sarin, and Chlorosarin) over the pure and cobalt-decorated C24 fullerenes by using PBE functional. For all four species, adsorption on the C24 fullerene takes place via a Co_O=P dative covalent bond between a cobalt-decorated site. It has been indicated that the

TE D

interaction between the chemical warfare agents and the cobalt-decorated C24 fullerene is more stable than that of the pure C24 fullerene. We found that the cobalt-decorated C24 fullerene is the most favorable site for Sarin and Soman in comparison with Chlorosarin and

EP

Chlorosoman. In addition, the interactions of chemical warfare agents with cobalt-decorated

AC C

C24 fullerene can be responsible to alter of the structural and electronic properties.

Keywords: C24 fullerene, Co adsorption, Chemical toxic agents, Density functional theory, Thermodynamic stability

E-mail addresses: [email protected], [email protected]; Tel: +98938-4544921

1

ACCEPTED MANUSCRIPT 1. Introduction The detection and removing from chemical warfare agents (CWAs) are a substantial issue for our safety and health. The C4H10FO2P (isopropyl methylphosphonofluoridate) and C7H16FO2P (3, 3-dimethyl-2-buthyl methylphosphonoflouridate) or Sarin and Soman,

RI PT

respectively, are extremely toxic substances belong to the class of organophosphorus compounds which are highly threatening for the living organisms. The first utilization of Sarin was reported when Iraqi army bombarded Iranian citizens in August 1988. Japanese

SC

cult used liquid Sarin to attack a building in Matsumoto and the Tokyo subway in 1994 and 1995, respectively [1-4]. To this time, several investigations into the adsorption of chemical

M AN U

warfare agents have been reported computationally [5-7].

With the discovery of C60 fullerene by Korto et al. in 1985 [8], search for highly stable cages with broad range of applications have been significantly accelerated and reported. Several experimental and theoretical reports have been directed towards the implementation of

TE D

carbon nanostructures in applications like chemical sensor [9-11], hydrogen storage [12, 13], and biomedical applications [14, 15]. Carbon nano-cages can be utilized to improve and advance toxic gas detection technologies; it is substantial to comprehend how each type of

EP

device may be appropriate [16, 17]. Rather and Wael investigated fullerene-C60 sensor for

AC C

ultra-high sensitive detection of bisphenol molecule by green technology [18]. For example, Shariatinia and Shahidi performed a DFT study upon the physical adsorption of cyclophosphamide derivatives on fullerene C60 surface [19]. Very recently, Sun and coworkers [20] have shown the interactions of low molecular weight natural organic acids with C60 in vacuum and water using Molecular dynamics simulations (MD). Mahdy et al. reported that the transition metal doped fullerene C60 systems introduced as a promising candidate for the detection of CO and NO gases compared to the pure C60 [21]. Xu and co-workers theoretically investigated the structural and electronic properties of C24-glycine and Gd@C24-

2

ACCEPTED MANUSCRIPT glycine using the hybrid DFT-B3LYP functional [22]. It was found that the adsorption of glycine from the amino nitrogen is energetically favorable to interact with C24 nano-cage. Baei and co-workers reported a theoretical study of nicotine adsorption in different configurations over C24 nano-cage using B3LYP functional [23]. They found that the C24 nano-cage can be

RI PT

used as electronic sensor for nicotine detection. Recently, the interaction between CWAs and materials was studied mainly upon the adsorptions of Soman and Sarin over Amorphous silica [24], γAl2O3 [25], Clay mineral fragments [26], and hydroxide brucite [27]. To the best of our

SC

knowledge, few reports have been performed on the adsorption of Soman (SOF), Chlorosoman (SOC), Sarin (SAF), and Chlorosarin (SAC) over carbon nanostructures and also Co decorated

M AN U

carbon nanostructure so in this research we are studying the possibility of application of C24 and its Co decorated as an adsorbents or sensor materials for the removal or detection of these toxic agents.

TE D

2. Methodologies

In the present work, all the calculations for the toxic agents (Soman, Chlorosoman, Sarin, and Chlorosarin) adsorption over pure and Co-doped C24 fullerenes have been carried

EP

out using GAUSSIAN 03 program package [28]. Subsequently, all the structures have been

AC C

optimized using density functional theory (DFT) with the Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional [29] and 6-311++G** standard basis set. For the Co atom, the standard LANL2DZ basis set was used throughout the calculations [30]. The PBE density functional has been previously reported to study carbon nanostructures [31-33]. Natural bond orbital (NBO) analysis was calculated by using PBE functional. Molecular electrostatic potential (MEP), HOMO and LUMO orbitals, and spin density plots of molecular configurations was generated with the GaussView [34] program (Gaussian, Inc.

3

ACCEPTED MANUSCRIPT Pittsburgh, PA). The corresponding binding energies (Eb) of the mentioned systems were determined through the following equation: Eb= Eadsorbent-C24– (Eadsorbent+ E C24) + BSSE

(1)

Eb= Eadsorbent-decorated-C24 – (Eadsorbent+ E Decorated-C24) + BSSE

(2)

RI PT

where Eadsorbent-C24 and Eadsorbent-doped-C24 are the total energies of soman, chlorosoman, sarin and chlorosarin as adsorbates over the pure and Co-decorated C24 fullerenes as adsorbents. Eadsorbate is the total energies of the soman, chlorosoman, sarin and chlorosarin molecules and

SC

E C24 and EDecorated-C24 is the total energies of the pure and Co-decorated C24 fullerenes. The basis set superposition error (BSSE) for the binding energy was computed by implementing

M AN U

the counterpoise method. HOMO-LUMO gap energies, Natural bond orbital (NBO), molecular electrostatic potential (MEP), partial density of states (PDOS) analyses, and quantum molecular descriptors were calculated accordingly at the same level of theory

3. Results and discussion

TE D

according to previous literature reports [35-40].

3.1. Adsorption of the chemical toxic agents on C24 fullerene

EP

We have considered here several adsorption configurations for per one of the toxic

AC C

agents including Sarin (SAF), Chlorosarin (SAC), Soman (SOF), Chlorosoman (SOC) which participate in the interaction with the pure C24 nano-cage. After full structural optimization without any constraints, the most stable structures are obtained upon the relaxation processes which are shown in Figs. 1 (See Figure 1). We performed full structural relaxation in the C24 fullerene with each of the toxic agents to examine the equilibrium geometries, energetic and electronic properties in Tables 1 and 2. After full optimization, we found the binding energies of Sarin and Chlorosarin on the surface of C24 fullerene is equal to -0.02 and -0.03 eV, respectively. In contrast, Soman and Chlorosoman adsorptions on the surface of C24

4

ACCEPTED MANUSCRIPT fullerene predict values of -0.03 and -0.04 eV, respectively. For all complexes, there is a little different from binding energy and the equilibrium height. As can be seen in Table 2, the equilibrium height of toxic agents with the C24 fullerene is 3.122 Å (SAC), 3.307 Å (SAF), 3.217 Å (SOC), and 3.318 Å (SOF) which predicts weak physisorption owing to van der

RI PT

Waals forces between participating specious in the interaction process. Our calculations present a small change in isometric structure of C24 fullerene. Vaise and co-workers [27] calculated the physisorption of Sarin on brucite with an amount of -0.26 eV. Javan et al.

SC

studied the physisorption of Soman and Chlorosoman over the outer surface BN nanotube with the energy values of -0.052 and -0.098 eV at B3LYP functional, while these values at

M AN U

M06-2X functional are -0.21 and -0.14 eV, respectively [41]. The NBO analysis shows a charge about 0.016 and 0.010 electrons are transferred from C24 fullerene to SOC and SOF and the charges about 0.015 and 0.010 electrons are transferred from C24 fullerene to SAC and SAF, respectively, indicating that the chemical toxic agents do not interact strongly

TE D

(physical adsorption) with the exterior surface of C24 fullerene.


EP



AC C

3.2. Adsorption of Sarin and Chlorosarin on Co-doped C24 The adsorptions of Sarin (SAF) and Chlorosarin (SAC) on the outer surface of Codecorated C24 fullerene were studied at the GGA level (PBE) as described in Figure 2. Geometrical parameters and adsorption energies are calculated and reported in Tables 1 and 3. Addition of cobalt impurity to nano-cluster leads to symmetry change from C4h (in the case of C24) to C4v (in the case of Co-decorated C24) splits the 3d shell into a doublet of eg (dxz, dyz) symmetry and three singlets of a1g (dz2), b2g (dxy) and b1g (dx2-y2) symmetry. The calculated value of binding energy of cobalt adsorbed on the surface of C24 fullerene is negative (-0.53

5

ACCEPTED MANUSCRIPT eV), showing the interaction between two specious is exothermic in nature. Baei and coworkers [42] reported the doping of Co atom at the center of C20 fullerene is endothermic with an energy value of +2.81 eV. Manna and co-workers [43] have shown the encapsulation of a plutonium atom within C24 leads to highly stable Pu@C24 fullerene. After structural

RI PT

optimization, we found the chemisorption energy of one SAC by using O atom upon the Codecorated C24 about -1.29 eV and 0.11 electrons charges were transferred from the adsorbate to the fullerene as Co atom act as Lewis acid, accepting electron pair from Lewis base (O

SC

atom of SAC). Therefore, the interactions between Co 3d electrons and O 2p lone pair electrons both raise the energy of the 3d sublevel and lift the degeneracy of the 3d orbitals.

M AN U

The natural electron configuration of cobalt atom in Co-decorated C24 fullerene is altered from [Ar] 4s0.07 3d3.41 to [Ar] 4s0.30 3d7.93 and [Ar] 4s0.303d8.00 after the interaction with SAF and SAC chemical agents, respectively, while the natural electron configurations of oxygen and phosphorous atoms are found to be [He] 2s0.89 2p2.64 and [He] 3s0.42 3p0.87 on the

TE D

adsorption of SAF with Co-decorated C24 fullerene and He] 2s0.89 2p2.64 and [He] 3s0.47 3p0.98 on the adsorption of SAC with Co-decorated C24 fullerene. The ∆G and ∆H values for SAC adsorbed from its phosphonyl oxygen atom (P=O) are calculated to be -1.60 and -2.10 eV,

EP

respectively. The values of ∆G and ∆H are negative. As results exhibited the value of ∆G is

AC C

slightly reduced in comparison with ∆H value that is owing to entropic effect. The interaction distance of SAC with the C24 fullerene decorated with Co atom is 1.949 Ả. In contrast, when a single SAF adsorbed on Co-decorated C24, the binding energy and interaction distance were about -1.60 eV and 1.991 Ả, respectively. The ∆G and ∆H values for SAF adsorbed on Codecorated C24 fullerene are calculated to be -2.06 and -2.61 eV, respectively, showing that the adsorption of SAF is much stronger than that of the SAC. NBO analysis demonstrated the charges about 0.10 electrons transferred from the SAF to the Co-decorated C24 fullerene. With respect to the Co-decorated C24 fullerene, the point charges upon the Co and C atoms

6

ACCEPTED MANUSCRIPT were 0.702 and -0.154 e, respectively. When SAC adsorbed on the Co-decorated C24 fullerene, the point charges of C and Co atoms of the nano-cage altered to 0.681 and -0.151 electrons compared to the basic values of 0.759 and -0.150 electrons in the SAF adsorbed on the fullerene, respectively. Sharma and Kakkar found the chemisorption energy of

RI PT

approximately -2.27 eV for SAF on Ti-doped MgO nanotube [44]. The infrared (IR) spectrum (Figure 3) of the SAC indicates the peaks at 502 and 1137 cm-1 which are in related to P-Cl and P=O bond stretching, which is in agreement with the results mentioned by Vaiss

SC

and co-workers [27]. The P=O vibrational frequency of a single SAF varies from 1249 cm-1 to 1185 cm-1 after adsorption process of Co-decorated C24 nano-cage, and result was in good

M AN U

agreement with the experimental data (P=O: 1277 cm-1 at liquid and 1280 or 1240 cm-1 at dilute aqueous solutions) as previously mentioned by Kuiper, Elkine, and Braue et al. [4547]. We observed the medium peak at around 906 cm-1 that is related to antisymmetric C-C stretching vibration between the 6-numbered ring and 12-numbered ring of the Co-decorated

TE D

C24 fullerene, which is close to the mentioned previous result by Xu et al. [22]. The C-H band of SAF adsorbed on Co-decorated C24 fullerene is observed in the range of 2950-3102 cm-1. Quantum molecular descriptors analysis also revealed that increase in the values of

EP

global hardness (η) from 0.60 and 0.36 eV (before adsorption) to 2.92 and 0.46 eV (after

AC C

adsorption) for the pure and Co-decorated C24 fullerenes, respectively, lead to increase in stability and decrease in reactivity of the system [48, 49]. Moreover, the same adsorption forms of SAF and SAC resulted to the reduction of values of electron affinity (A) and electrophilicity index (ω) which is in well agreement with the results obtained by the ∆Nmax calculation whereas electron transfers occur from the adsorbates to adsorbents by SAF and SAC interactions with the pure and Co-decorated C24 fullerenes (see Tables 2 and 3).


7

ACCEPTED MANUSCRIPT
3.3. Adsorption of Soman and Chlorosoman on Co-doped C24 After structure relaxations, we calculated two kinds of energetically favorable configurations, in which the Soman (SOF) and Chlorosoman (SOC) adsorb over the Co-decorated C24

RI PT

fullerene as displayed in Figure 2. In Table 3, For SOF and SOC upon the Co-decorated C24 fullerene, the calculated binding energy for both states were about -1.31 and -1.65 eV, and distance of SOF and SOC to the Co-decorated C24 were discovered to be 1.950 and 1.996 Ǻ,

SC

respectively. The BSSE-corrected binding energies for interaction between boron-doped C34 fullerene and Soman molecule was found to be -1.83 eV (from its O head) and -1.65 eV

M AN U

(from it’s OF head) [50]. The obtained results demonstrated a strong chemical interaction (covalent) with a small distance in SOF or SOC on Co-decorated C24 fullerene. The natural electron configuration of cobalt atom in Co-decorated C24 fullerene interacting with SOF and SOC are found to be [Ar] 4s0.19 3d4.540 and [Ar] 4s0.143d4.68, respectively, while the natural

TE D

electron configurations of oxygen and phosphorous atoms are found to be [He] 2s0.89 2p2.63 and [He] 3s0.42 3p0.87 on the adsorption of Soman with Co-decorated C24 fullerene and [He] 2s0.85 2p2.56 and [He] 3s0.47 3p0.98 on the adsorption of SOC with Co-decorated C24 fullerene.

EP

The values of ∆G and ∆H for SOC adsorbed on Co-decorated C24 fullerene are calculated to

AC C

be -1.51 and -2.06 eV, while these values for Soman are about -2.15 and -2.68 eV, respectively. When SOF chemisorbed upon Co-decorated C24 surface, the length of P=O and P-F bonds are equal to 1.513 and 1.611 Ả, respectively, which are different of the isolated SOF molecule with the values of 1.489 (P=O) and 1.631 Ả (P-F). Compared to Soman molecule, the length of P=O and P-Cl bonds of SOC chemisorbed on Co-decorated C24 surface were approximately 1.526 and 2.064 Ả, respectively, while it was 1.493 (P=O) and 2.102 Ả (P-Cl) in isolated SOC molecule. The NBO analysis demonstrates that the charges about 0.301 and 0.116 electrons were transferred from SOF and SOC to Co-decorated C24

8

ACCEPTED MANUSCRIPT fullerene. The infrared spectra for SOF and SOC are very resembling, showing the similarities in their molecular structure and surface bonding energy. After adsorption of SOF and SOC on the surface of Co-decorated C24 fullerene, the bands observed at around 1168 and 1121 cm-1 can be attributed to the P=O vibrations and compared they with the values of

RI PT

the SOF and SOC appears at around 1234 and 1227 cm-1 , respectively [4]. The bands at around 495 and 789 cm-1 are assigned to the P-Cl and P-F stretching vibrations (Figure 3). Davis et al. calculated the band of P-F stretching vibration in related to Sarin adsorbed on

SC

silica in the region of 803 cm-1 [4]. The global hardness values decreased in SOC/C24 and SOF/C24 systems and increased in SOC/Co-C24 and SOF/Co-C24 systems, indicating that

M AN U

stability of the recently mentioned systems increased after adsorptions of SOF and SOC over C24 fullerene while it dwindled after adsorptions of SOF and SOC over Co-C24 fullerene which as the consequence resulted to the diminishment and growth of reactivity of C24 and Co-decorated C24 species, respectively. Concerning electron transfer studies, both electron

TE D

affinity (A) and electrophilicity index (ω) values reduced upon the adsorption processes except in case of SOF/C24 system where the ω increased from 21.59 to 43.13 eV. This along with the values of ∆Nmax as total electron transfer index of the studied systems represent that

EP

electrons were generally transferred from SOF and SOC to the adsorbents which coincides

AC C

well with our previous report [41].

3.4. Electronic properties of the Co-decorated C24 complexes 3.4.1. HOMO-LUMO energy of the complexes In this section, we are reporting on the frontier molecular orbitals (FMO) of the applied adsorbents before (see Figure 4) and after adsorption (see Figures 5 and 6), as it is represented in Figure 5, the occupied (HOMO) and virtual (LUMO) orbitals in both α (alpha) and β (Beta) forms are more situated over the Co-decorated C24 fullerene after interactions

9

ACCEPTED MANUSCRIPT with SAF and SAC molecules except in β forms of LUMO where the more orbitals concentrated on Co atom and its surrounding atoms. Moreover, the calculations of electronic properties revealed that the HOMO energy (EHOMO) and LUMO energy (ELUMO) of C24 has some negligible changes after SAF and SAC adsorptions, while EHOMO and ELUMO of Co-

RI PT

decorated C24 changed appreciably after adsorption of the same molecules (see Table 2). The value of Eg for C24 fullerene is found to be 1.24 eV at PBE functional, which is close to the obtained values by other literatures [51, 52]. When the nano-cage linked to cobalt atom, the

SC

Eg value is reduced to 0.71 eV at alpha form and 0.31 eV at Beta form which this push down expect is related in the change of HOMO energy in comparison with the LUMO energy.

M AN U

Moreover, the differences in EHOMO and ELUMO values (∆Eg) in α forms were 28.17 (SAF) and 29.58% (SAC) and in β forms were 8.57 (SAF) and 42.86% (SAC) whereas the EHOMO and ELUMO values in the same parameter is about zero concerning SAF and SAC adsorptions over C24 fullerene. Figure 5 demonstrates that β form of LUMO is more localized over the

TE D

Co atom and C atoms close to it (after SAF interaction with Co-decorated C24), while the occupied and virtual orbitals of both α and β forms were uniformly distributed throughout the adsorbents.

EP

Concerning HOMO-LUMO of Co-decorated C24 fullerene before and after interactions

AC C

with SOF and SOC, Figure 6 represents that, except β forms of LUMO in both SOF and SOC/Co-decorated C24 systems which electrons are more exist on Co atom and areas close to it, occupied and virtual orbitals of both α and β forms were uniformly distributed throughout the adsorbents. Again, the frontier molecular orbitals (FMOs) values of SOF and SOC adsorbed over C24 were not changes upon the adsorption process while the HOMO and LUMO in both α and β forms reduced in both SOC/CoC24 and SOF/CoC24 systems. As displayed in Table 3, the Eg values of SOC/CoC24 and SOF/CoC24 systems are larger than that of the C24 fullerene (0.71 eV), being 0.91 and 0.92 eV, respectively. It is noticed that, a

10

ACCEPTED MANUSCRIPT large Eg can be related to kinetic and structural stability [53-58]. The values of ∆Eg increased from both α and β forms of SOC/CoC24 and SOF/CoC24 systems while it was again close to zero in case of SOC/C24 and SOF/C24 systems. As shown in Table 4, in the cobalt-decorated C24 nano-cage, five 3d subshells are highly

RI PT

polarized with the dominance of spin-up states. Among five 3d subshells, 3dyz and 3dx2-y2 are the most polarized orbital, while we also notice significant spin polarizations in 3dxz and 3dxy. For all the other species, high positive spin polarization terms were found in all 3d subshells

SC

which resulted in a more severe distortion for the species (see Table 4). Therefore, the planewave calculations predicted that the structures gave antiferromagnetic contributions. The

M AN U

spin-up states in the 3d subshells and polarizations show the effect of strong chemical toxic agent’s dispersion interactions. According to the data presented we can see that the total magnetic moment of the SAF/Co-C24 and SAC/Co-C24 complexes have slightly increased in comparison with Co-C24 host cage. However in the case of the SAF/Co-C24 and SOC/Co-C24

TE D

the total magnetic moment have also a slightly decline (Figure 7). As presented in Figure 5 the Co atom spin polarization is significant and the Co-3d orbitals charge polarization has

EP

remarkable role in induced magnetic moment to the systems.


AC C

The molecular electrostatic potential (MEP) maps for all chemical toxic agents (in the most stable states) interacting with Co-decorated C24 fullerene are displayed in Figures 4 and 7. MEP maps represent nucleophilic sites with blue color for the most positive electrostatic potential (chemical toxic agents) and the positions of electrophilic site with red color is defined as the most electronegative electrostatic potential (Co-decorated C24 fullerene). Hence, MEP maps indicate that the all chemical toxic agents acts as an electron donor and Co-decorated C24 fullerene acts as an electron acceptor.


11

ACCEPTED MANUSCRIPT
3.4.2. Total densities of states (TDOSs) of the complexes

RI PT

The total density of state (TDOS), partial density of state (PDOS), and overlap population density of state (OPDOS) were computed for stable configurations of considered complexes adsorbed upon the Co-decorated C24 fullerene. All the TDOS, PDOS, and OPDOS

SC

plots for each complex are displayed in Figure 8. PDOS plots of first (fullerene) and second (chemical toxic agents) fragments indicate with the red and blue colors of specified

M AN U

complexes. According to the partial density of state demonstrated in the figure we can also see that the HOMO and LUMO states of the Co-decorated C24 states interacting with chemical toxic agents are significantly altered. This suggests the chemical toxic agents are chemically adsorbed (chemisorption) upon the outer surfaces of Co-decorated C24 fullerene

TE D

and consequently, charge transfer does happen between adsorbent and adsorbate. As a result decoration of C24 surface with cobalt can leads to a significant reduces in Eg of the fullerene after the interacting with SAF and SOF in both alpha and beta forms, which can resulting

EP

from overlapping of HOMO or LUMO orbital of chemical toxic agents with Co-C24 orbital

AC C

hybridization [59, 60]. Also in all considered systems, the overlap population density of states (OPDOS) indicates bonding interaction (green color) between adsorbate and adsorbent.


3.4.3. Electron charge density differences We demonstrate in Figures 9 and 10 contour plots of the electron charge density distribution belong to SAC (I), SAF (II), SOC (III), and SOF (IV) interacting with Codecorated C24 fullerene. As it is explicit from Figures 9 and 10, the electron charge density distribution is obviously centralized around the cobalt atom and the oxygen atom of the

12

ACCEPTED MANUSCRIPT complex and indicated a significant overlapped charge density in agreement with covalent interaction scheme leading to the chemical adsorption and high binding energy of the chemical warfare agents to the Co-decorated C24 fullerene.

3.4.4. Electron localization function differences

RI PT



To study further the performance of cobalt atom, we have distinguished the isosurface

SC

plots of the electron localization function (ELF) for the SAC (I), SAF (II), SOC (III), and SOF (IV) adsorbed on cobalt-decorated C24 fullerene. Figures 11 and 12 demonstrate in all

M AN U

cases, the electron localization is significant between oxygen atom of chemical warfare agents and Co atom of adsorbent surface as the electrons more localized on the oxygen atom of adsorbate interact considerably with the vacant orbitals of the Co atom, suggesting that Co facilities the adsorption of chemical warfare agents over the surface of C24 fullerene. In

TE D

Figures 11 and 12 for all cases demonstrate that the interaction between the cobalt atom of the nano-cage and the oxygen atom of chemical warfare agents (Co-O bond) is a strong negative covalent potential in nature, exhibited of high polarization, and in agreement with

AC C

EP

the high dipole moment of all the systems.



4. Summary

The interaction of Soman, Chlorosoman, Sarin, and Chlorosarin over the pure and cobaltdecorated C24 fullerenes have been studied theoretically using DFT calculations to investigate the effect of cobalt atom in the binding energies and the electronic properties in the mentioned systems. The results are as follows:

13

ACCEPTED MANUSCRIPT (1) In the interaction between the chemical warfare agents and the cobalt-decorated C24 fullerene, the adsorption resulted in chemical strong Co-O bonding and charge transfer from adsorbate to adsorbent. (2) The results demonstrated that the energy gap between the HOMO and LUMO orbitals of

RI PT

the mentioned systems decreases quite rapidly as the chemical warfare agents approaches the cobalt-decorated C24 fullerene.

(3) The reported harmonic frequencies of the chemical warfare agents interacting with cobalt-

SC

decorated C24 fullerene indicated a positive vibration mode, suggesting that the complexes are stable. Also, the obtained results by thermodynamic parameters demonstrated the interaction

M AN U

between adsorbate and adsorbent is exothermic and energetically favorable.

Acknowledgments

We would like to thank the clinical Research Development Unit (CRDU), Sayad Shirazi

TE D

Hospital, Golestan University of Medical Sciences, Gorgan, Iran.

References

EP

[1] C. P. Holstege, M. Kirk, F. R Sidell, Crit. Care Clin. 13 (1997) 923-942.

AC C

[2] Y. Tokuda, M. Kikuchi, O. Takahashi, G. H. Stein, Resuscitation 68 (2006) 193-202. [3] T.M. Alam, C.J. Pearce, J.E. Jenkins, Computational and Theoretical Chemistry 995 (2012) 24-35.

[4] E.D. Davis, W.O. Gordon, A.R. Wilmsmeyer, D. Troya, J.R. Morris, J. Phys. Chem. Lett. 5 (2014) 1393-1399. [5] V. Stengl, J. Bludska, F. Oplustil b, T. Nemec, Materials Research Bulletin 46 (2011) 2050-2056. [6] V. M. Bermudez, J. Phys. Chem. C 111 (2007) 3719-3728.

14

ACCEPTED MANUSCRIPT [7] Q. Zhao, G. Liu, H. Zhang, F. Zhou, Y. Li, W. Cai, Journal of Hazardous Materials 324 (2017) 194-202. [8] H. Kroto, J. Heath, S. Obrien, R. Curl, R. Smalley, Nature 318 (1985) 162-163. [9] S.W. de Silva, A. Du, W. Senadeera, Y. Gu, Beilstein J. Nanotechnol. 5 (2014) 413-418.

RI PT

[10] Z. Zhang, C.H. Turner, J. Phys. Chem. C 117(17) (2013) 8758-8766. [11] W. An, C.H. Turner, Chemical Physics Letters 482 (2009) 274-280.

Wu, Chemical Physics Letters 511 (2011) 393-398. [13] S.S. Han, H.M. Lee, Carbon 42 (2004) 2169-2177.

SC

[12] C. Tian, Z. Wang, M. Jin, W. Zhao, Y. Meng, F. Wang, W. Feng, H. Liu, D. Ding, D.

M AN U

[14] S.H. Friedman, D.L. Decamp, R.P. Sijbesma, G. Srdanov, F. Wudl, G.L. Kenyon, Inhibition of the HIV-1 protease by fullerene derivatives: model building studies and experimental verification, J. Am. Chem. Soc. 115 (1993) 6506-6509. [15] H. Kato, Y. Kanazawa, M. Okumura, A. Taninaka, T. Yokawa, H. Shinohara, J. Am.

TE D

Chem. Soc. 125 (2003) 4391-4397.

[16] M.T. Baei, Heteroatom Chemistry 24 (2013) 516-522. [17] A. Soltani, M.T. Baei, M. Ramezani Taghartapeh, E. Tazikeh Lemeski, S. Shojaee,

EP

Struct Chem 26 (2015) 685-693.

AC C

[18] J.A. Rather, K. De Wael, Sensors and Actuators B 176 (2013) 110-117. [19] Z. Shariatini, S. Shahidi, Journal of Molecular Graphics and Modelling 52 (2014) 71-81. [20] Q. Sun, H-B. Xie, J. Chen, X. Li, Z. Wang, L. Sheng, Chemosphere 92 (2013) 429-434. [21] A.M. El Mahdy, Applied Surface Science 383 (2016) 353-366. [22] L. Xu, C. Li, F. Li, X. Li, S. Tao, Spectrochimica Acta Part A 98 (2012) 183-189. [23] M.T. Baei, P. Torabi, H. Mohammadian, H. Saeedeh, Fullerenes, Nanotubes and Carbon Nanostructures 23 (2015) 874-877.

15

ACCEPTED MANUSCRIPT [24] E.D. Davis, W.O. Gordon, A.R. Wilmsmeyer, D. Troya, J.R. Morris, J. Phys. Chem. Lett. 2014, 5, 1393-1399. [25] V. M. Bermudez, J. Phys. Chem. C 113 (2009) 1917-1930.

Leszczynski, J. Phys. Chem. B 110 (2006) 21175-21183.

RI PT

[26] A. Michalkova, J. Martinez, O. A. Zhikol, L. Gorb, O. V. Shishkin, D. Leszczynska, J.

[27] V.S. Vaiss, I. Borges Jr., A.A. Leitao, J. Phys. Chem. C 115 (2011) 24937-24944.

[28] M. Frisch, G. Trucks, H. Schlegel, G. Scuseria, M. Robb, J. Cheeseman, et al., Gaussian

SC

03, revision D. 01, Gaussian Inc, Wallingford, CT, (2004).

[29] J.P. Perdew, K. Burke, M. Ernzerhof, Phys. Rev. Lett. 77 (1996) 3865-3868.

M AN U

[30] H-S. Wu, F-Q. Zhang, X-H. Xu, C-J. Zhang, H. Jiao, J. Phys. Chem. A 107 (2003) 204209.

[31] A. Soltani, M. Bezi Javan, RSC Adv. 5 (2015) 90621-90631.

[32] A. Soltani, A. Sousaraei, M. Bezi Javan, M. Eskandaric, H. Balakheyli, New J. Chem.

TE D

40 (2016) 7018-7026.

[33] D.S. Sabirov, R.G. Bulgakov, Computational and Theoretical Chemistry 963 (2011) 185-190.

EP

[34] GaussView 4.1.2, Gaussian Inc., Wallingford, CT, 2004.

AC C

[35] R.G. Parr, R.A. Donnelly, M. Levy, W.E. Palke, J. Chem. Phys. 68 (1978) 3801–3807. [36] R.G. Parr, L. Szentpaly, S. Liu, J. Am. Chem. Soc. 121 (1999) 1922. [37] T. Koopmans, Über die Zuordnung von Wellenfunktionen und Eigenwerten zuden Einzelnen Elektronen Eines Atoms, Physica 1 (1933) 104. [38] Mohammad T. Baei, A. S. Ghasemi, E. Tazikeh Lemeski, Alireza Soltani, Niloofar Gholami, J Clust Sci 27 (2016) 1081-1096. [39] M. Bezi Javan, A. Soltani, E. Tazikeh Lemeski, A. Ahmadi, S. Moazen Rad, Superlattices and Microstructures 100 (2016) 1-14.

16

ACCEPTED MANUSCRIPT [40] A. Soltani, S. Ghafouri Raz, V. Joveini Rezaei, A. Dehno Khalaji, M. Savar, Applied Surface Science 263 (2012) 619-625. [41] M. Bezi Javan, A. Soltani, A.S. Ghasemi, E. Tazikeh Lemeski, N. Gholami, H. Balakheyli, Applied Surface Science 411 (2017) 1-10.

RI PT

[42] M.T. Baei, A. Soltani, P. Torabi, F. Hosseini, Monatsh Chem 145 (2014)1401-1405. [43] D. Manna, A. Sirohiwal, T.K. Ghanty, J. Phys. Chem. C 118 (2014) 7211-7221. [44] N. Sharma, R. Kakkar, Journal of Computational Science 10 (2015) 225-236.

SC

[45] A.E.T. Kuiper, J.J.G.M. van Brokhoven, J. Medema, J. Catal. 43 (1976) 154-167.

[46] M.D. Elkin, S.I. Tatarinov, E.Y. Stepanovich, J. Appl. Spectrosc. 77 (4) (2010) 479-483.

M AN U

[47] E.H. Braue Jr., M.G. Panella, Appl. Spectrosc. 44 (9) (1990) 1513-1520. [48] N. Saikia, R.C. Deka, Comput. Theor. Chem. 964 (2011) 257-261. [49] N. Saikia, R.C. Deka, Chem. Phys. Lett. 500 (2010) 65-70.

[50] S. Arshadi, A.R. Bekhradnia, F. Abbasi Alaei, Phosphorus, Sulfur, and Silicon 190

TE D

(2015) 2051-2061.

[51] Y. Yue, X.J. Li, Int. J. Quantum. Chem. 111 (2009) 96-102. [52] S. Peng, Y. Zhang, X.J. Li, Y. Ren, D.X. Zhang, Spectrochim. Acta A 74 (2009) 553.

EP

[53] D. Bucher, Chem. Phys. Lett. 534 (2012) 38-42.

AC C

[54] S. Peng, X.J. Li, Spectrochim. Acta A 73 (2009) 67. [55] X.J. Li, J. Mol. Struct. (THEOCHEM) 896 (2009) 25. [56] S. Peng, Y. Zhang, X.J. Li, Y. Ren, D.X. Zhang, Spectrochim. Acta A 74 (2009) 553. [57] Z. Sun, X. Li, M. Tian, G. Zhao, J. Li, B. Ma, J. Mol. Struct. (THEOCHEM) 913 (2009) 265. [58] D.E. Manolopoulos, J.C. May, S.E. Down, Chem. Phys. Lett. 181 (1991) 105. [59] O. Leenaerts, B. Partoens, F. Peeters, Phys. Rev. B 77 (2008) 125416-125421.

17

ACCEPTED MANUSCRIPT [60] H. Pinto, R. Jones, J.P. Goss, P.R. Briddon, J. Phys.: Condens. Matter 21 (2009)

AC C

EP

TE D

M AN U

SC

RI PT

402001-402003.

18

ACCEPTED MANUSCRIPT Figure 1. Adsorption configurations of Sarin (a), Chlorosarin (b), Soman (c), Chlorosoman

AC C

EP

TE D

M AN U

SC

RI PT

(d) on the pure C24 fullerene.

19

ACCEPTED MANUSCRIPT Figure 2. Adsorption configurations of Sarin (e), Chlorosarin (f), Soman (g), Chlorosoman

AC C

EP

TE D

M AN U

SC

RI PT

(h) on the Co-doped C24 nano-cage.

20

ACCEPTED MANUSCRIPT

AC C

EP

TE D

M AN U

SC

RI PT

Figure 3. IR spectra of all the complexes.

21

ACCEPTED MANUSCRIPT Figure 4. Illustrations of FMO, MEP plots and spin density from the cobalt-decorated

AC C

EP

TE D

M AN U

SC

RI PT

C24 fullerene.

22

ACCEPTED MANUSCRIPT Figure 5. The HOMO and LUMO orbitals of the Sarin and Chlorosarin adsorbed on the

AC C

EP

TE D

M AN U

SC

RI PT

surface of Co-decorated C24 fullerene.

23

ACCEPTED MANUSCRIPT Figure 6. The HOMO and LUMO orbitals of the Soman and Chlorosoman adsorbed on the

AC C

EP

TE D

M AN U

SC

RI PT

surface of cobalt-decorated C24 fullerene.

24

ACCEPTED MANUSCRIPT Figure 7. Illustrations of MEP plots and spin density transfer from chemical toxic agents to

AC C

EP

TE D

M AN U

SC

RI PT

the cobalt-decorated C24 fullerene.

25

ACCEPTED MANUSCRIPT Figure 8. Parital density of states spectrums for (a) chlorosarin and (c) chlorosoman (left)

AC C

EP

TE D

M AN U

SC

RI PT

and (b) Sarin and (d) Soman (right) interacting with the Co-decorated C24 fullerene.

26

ACCEPTED MANUSCRIPT Figure 9. Contour plots of charge density for the Sarin and Chlorosarin adsorbed on the

AC C

EP

TE D

M AN U

SC

RI PT

surface of cobalt-decorated C24 fullerene.

27

ACCEPTED MANUSCRIPT Figure 10. Contour plots of charge density for the Soman and Chlorosoman adsorbed on the

AC C

EP

TE D

M AN U

SC

RI PT

surface of cobalt-decorated C24 fullerene.

28

ACCEPTED MANUSCRIPT Figure 11. ELF plots of Sarin (I) and Chlorosarin (II) adsorbed on cobalt-decorated C24

AC C

EP

TE D

M AN U

SC

RI PT

fullerene.

29

ACCEPTED MANUSCRIPT Figure 12. ELF plots of Soman, and Chlorosoman adsorbed on cobalt-decorated C24

AC C

EP

TE D

M AN U

SC

RI PT

fullerene.

30

ACCEPTED MANUSCRIPT Table 1 Optimized structural geometries of Sarin, Chlorosarin, Soman and Chlorosoman chemical agents interacting with the pure and Co-decorated C24 fullerenes using the PBE functional. Lengths are in angstrom (Å).

O=P-O/°

C-Co-O/°

-

117.4

-

-

117.7

-

-

117.2

-

-

-

117.9

-

-

-

-

-

-

-

117.5

-

1.490

-

-

117.6

-

-

1.491

-

-

117.5

-

1.495

2.099

1.492

-

-

117.8

-

-

-

-

1.542

1.952

-

-

-

1.614

1.510

-

2.156

1.827

1.991

117.4

153.1

1.520

2.071

1.540

1.940

1.950

116.7

150.3

1.525

2.063

1.542

1.938

1.950

116.8

148.2

1.513

-

2.158

1.828

1.995

117.4

150.4

P=O/Ǻ

P-Cl/Ǻ

C-C/Ǻ

Co-C/Ǻ

Sarin

1.631

1.489

-

-

-

-

1.492

2.103

-

-

1.616

1.492

-

-

-

Chlorosoman

-

1.493

2.102

-

C24

-

-

-

1.492

1.630

1.491

-

1.491

-

1.495

2.101

1.629

1.491

Chlorosoman/C24

-

Co-C24

Sarin/C24 Chlorosarin/C24 Soman/C24

Sarin/Co-C24 Chlorosarin/Co-C24

-

Soman/Co-C24

1.611

AC C

EP

Chlorosoman/Co-C24

M AN U

Soman

TE D

Chlorosarin

SC

P-F/Ǻ

RI PT

Co-O/Ǻ

Property

31

ACCEPTED MANUSCRIPT Table 2 The structural and electronic parameters and quantum molecular descriptor of chemical warfare agents on the pure C24 fullerene.

C24

Sa(Cl)

Sa(F)

So(Cl)

So(F)

Sa(Cl)/C24

Sa(F)/ C24

So(Cl)/ C24

So(F)/ C24

EHOMO/eV

-5.69

-6.61

-7.00

-6.89

-6.60

-5.46

-5.46

-5.45

-5.52

ELUMO /eV

-4.49

-0.77

-0.79

-1.26

0.73

-4.26

-4.26

-4.24

-4.32

Eg/eV

1.20

5.84

6.21

5.63

7.33

1.20

1.20

1.21

1.20

D/Debye

0.0

3.13

3.00

3.30

2.94

3.74

3.68

3.80

3.13

EF/eV

-5.09

-3.69

-3.89

-4.07

-2.94

-4.86

-3.69

-4.85

-4.92

Ead/eV

-

-

-

-

-

-0.03

-0.02

-0.04

-0.03

D/Ǻ

-

-

-

-

-

-0.03

-0.02

-0.04

-0.03

I /eV

5.69

6.61

7.00

6.89

6.60

5.46

5.46

5.45

5.52

A /eV

4.49

0.77

0.79

1.26

0.73

4.26

4.26

4.24

4.32

η /eV

0.60

2.92

3.11

2.82

µ /eV

-5.09

-3.69

-3.90

-4.08

S /eV-1

0.83

0.17

0.16

0.18

ω /eV

21.59

2.33

2.44

2.95

∆Nmax

0.84

1.26

1.25

1.45

M AN U

SC

RI PT

Property

2.92

2.92

0.47

0.19

-3.67

-3.69

-3.69

-4.18

-4.05

0.17

0.117

0.117

1.063

2.63

2.29

1.593

1.593

18.573

43.13

1.25

1.26

1.26

8.89

21.32

AC C

EP

TE D

2.94

32

ACCEPTED MANUSCRIPT Table 3 The thermodynamic and electronic parameters and quantum molecular descriptor of chemical warfare agents on the Co-decorated C24 fullerene.

CoC24

Sa(Cl)/CoC24

Sa(F)/CoC24

So(Cl)/CoC24

So(F)/CoC24

αEHOMO/eV

-5.12

-4.61

-4.50

-4.60

-4.48

αELUMO /eV

-4.41

-3.70

-3.58

-3.69

-3.56

αEg/eV

0.71

0.91

0.92

0.91

0.92

∆Eg/%

-

28.17

29.58

28.17

29.58

EF/eV

-4.77

-4.16

-4.04

-4.15

-4.02

βEHOMO/eV

-4.79

-4.26

-4.40

-4.25

-4.39

βELUMO/eV

-4.44

-3.94

-3.90

-3.94

-3.92

βEg/eV

0.35

0.32

0.50

0.31

0.47

∆Eg/%

-

8.57

42.86

11.43

34.28

EF/eV

-4.62

-4.10

D/Debye

5.76

12.57

-

-1.29

1.950

1.939

∆G/eV

-

-1.60

∆H/eV

-

-2.10

I /eV

5.12

4.61

A /eV

4.41

η /eV

0.36

µ /eV

-4.77

S /eV-1

1.41

ω /eV

SC

M AN U

-4.16

12.63

12.54

12.84

-1.60

-1.31

-1.65

1.830

1.938

1.828

-2.06

-1.51

-2.15

-2.61

-2.06

-2.68

4.50

4.60

4.48

3.70

3.58

3.69

3.56

0.46

0.46

0.46

0.46

-4.16

-4.04

-4.15

-4.02

1.10

1.09

1.10

1.09

31.98

18.97

17.74

18.88

17.56

13.25

9.04

8.78

9.02

8.74

AC C

∆Nmax

-4.10

TE D

D/Ǻ

-4.15

EP

Ead/eV

RI PT

Property

33

ACCEPTED MANUSCRIPT Table 4 Spin polarization terms (µB) of the Co 3d orbitals in four investigated C24-Cochemical agent complexes obtained from Plane-Wave calculations using σ = 0.002 Ry.

3dz2

3dxz

3dyz

3dx2-y2

3dxy

Spin density

Co-C24

1.798

1.217

1.847

1.840

1.445

1.413

Sa(Cl)/Co-C24

1.771

1.438

1.693

1.600

1.498

1.463

Sa(F)/Co-C24

1.342

1.837

1.739

1.451

1.559

1.106

So(Cl)/Co-C24

1.841

1.552

1.343

1.491

1.707

0.994

So(F)/Co-C24

1.751

1.518

1.632

1.487

1.614

1.448

AC C

EP

TE D

M AN U

SC

RI PT

Property

34

ACCEPTED MANUSCRIPT

Highlights ● We have investigated the C24 fullerene for interacting with chemical toxic agents using DFT. ● Adsorption behavior of chemical toxic agents on Co-decorated C24 fullerene was studied.

RI PT

● Structural and electronic properties of the on Co-decorated C24 fullerene after adsorption was

AC C

EP

TE D

M AN U

SC

investigated.