Web of Science, ISC, CAS, Google Scholar

Document Type : Original Article


1 Department of Chemistry, College of Basic Sciences, Yadegar-e-Imam Khomeini (RAH) Shahre Rey Branch, Islamic Azad University, Tehran, Iran

2 Department of Chemistry, Faculty of Pharmaceutical Chemistry, Tehran Medical Sciences, Islamic Azad University, Tehran, Iran


This paper examined interaction of Graphene with Amoxicillin antibiotic through density functional theory (DFT) and by using molecular docking method. For this, the structures of Amoxicillin and, Graphene were initially optimized with Gaussian program. Then, by using the molecular docking strategy and its grading system, we computed the arrangement of 10 structures with additional negative binding energy and a fixed state compared with other samples. Finally, for the most fixed arrangement with Graphene, molecular orbitals evaluations were conducted, and binding energy along with thermodynamic evaluated, the results indicated that the adsorption of Amoxicillin antibiotic on Graphene was an exothermic. Finally, the QTAIM calculations were performed to evaluate the type of interaction and bonds created between amoxicillin and graphene.

Graphical Abstract

Interaction of Graphene with Amoxicillin Antibiotic by in Silico Study


Main Subjects


Graphene is a sheeting of genuine carbon, in which each atom is available for chemical reaction from both sides due to its flat structure [1]. It is used to design a new generation of chemical and biological sensors due to its outstanding physical and chemical properties, excellent catalytic activity, and low production cost [2]. The performance of various voltammetry biosensors based on graphene nanocomposites is well defined [3], and the use of graphene in nanocomposites in the application of biosensor due to its unique properties in the field of optimal electrical conductivity [4], flexibility, lightness, and good strength has improved the performance and progress in the properties of biosensor. Graphene also improves the linear periphery and the optimal correlation between concentration and flow. Likewise, the detection limit, the sensitivity of the sensor, and its response time are other things that the use of graphene has led to the improvement of their properties [3]. Among the major environmental problems are drug contaminants, including antibiotics, which are significantly present in the environment due to their high consumption [5]. Because after affecting the body of patients, they enter the wastewater treatment processes along with body excreta. Antibiotics prevent biological wastewater treatment and they cause a lot of problems in the environment [6] and are among the main causes of drinking water [7]. Therefore, it seems necessary to research sensors that detect antibiotics such as amoxicillin. In this paper, by using computational chemistry and molecular docking method [8], we examined the interaction of amoxicillin with graphene to investigate thermodynamic variables and the extent of heat from the absorption of antibiotic amoxicillin through these calculations to see the amoxicillin sensor (Figure 1) [5]. Therefore, we employed Gauss View software to develop the graphene sensor with amoxicillin. To minimize the energy of systems, first, the explored systems were geometrically optimized by the Gauss software via B3LYP model and 6-31G basis set [9]. However, the physical features of the systems were examined by using the Autodock Tools program, and finally, the nature of the resulting bonds was calculated via Quantum Theory of Atoms In the Molecule QTAIM [10]. Physical chemistry indexes and features were estimated, compared, and then the most steady conditions and positions of graphene were identified by the antibiotic amoxicillin.

Figure 1: a) Chemical Amoxicillin structure, b) Graphene structure; Experimental (carbon: gray, nitrogen: blue, oxygen: red, hydrogen: white, and sulfur: yellow)

Computational methods

Initially, amoxicillin and graphene structures were extracted by nanotube modeller [11], and Gauss View 5.0 software [12]. In the next step, geometric optimization was carried out via the density function hypothesis technique and computational level: B3LYP / 6-31G (d). This computational level was chosen since its findings were in good line with empirical data in previously conducted studies. All computations were conducted at a temperature and a pressure of 298 K and 1 atmosphere, respectively by using Gaussian software [13]. Then, it was used by Auto Dock software (Auto Dock, version 4.2). Auto Dock created a sequence of ten conformation designs, indicating the ten greatest anticipating models about the way the antibiotic could possibly cooperate with the Graphene. Auto Dock generated a sequence of energy values (binding energy, ligand proficiency, Van der Waals, hinderance constant, intermolecular energy, electrostatic, and overall inner energy) via a  program [8]. The process under investigation was expressed as the following formula:

Amoxicillin + adsorbent(G)→ Amoxicillin –adsorbent(G)

Results and Discussion

Discussion and results

To detect the most steady arrangement of amoxicillin with graphene, as can be seen in Figure 2, we employed the docking technique via Auto dock software [14]. Docking is the process of determining the orientation and binding energy of two compounds. In this work, the interaction was examined between amoxicillin and graphene as a receptor in its active site. Docking is a molecular mechanic computation. That automated molecular algorithm connects a smaller compound and specifies the (ligand) to the active site of the larger molecule (target). This method includes determining the composition orientation, the geometric structure of the conformation (formulation), and the ranking. Ranking can be a measure of connection energy, free energy, or a numerical qualitative measure. Each automatic docking algorithm attempts to put the combination in active position in different orientations and formulations and calculate a score for each [15].

Figure 2: Ten steady complexes (AG) created from molecular docking computations, Amoxicillin (A) with Graphene(G), because of their higher negative binding energy compared with other instances (Amoxicillin: pink, graphene: green)

Table 1 and Figure 3 indicate the results of amoxicillin derivatives with graphene, suggesting that the energy changes were not positive, so it reflects the exothermic of the adsorption procedure in these arrangements. Hence, we expect the physical adsorption of amoxicillin to occur in the interplay of graphene with amoxicillin. Furthermore, this nanostructure can be used for producing novel thermal sensors to calculate amoxicillin. These sensors generally measure the temperature variations resulting from the development of a procedure by using an extremely precise and sensitive thermistor that is used as an indicator to calculate the volume of analytes.

Table 1: In normal circumstance (T= 298.15 K, p=1 atm); The structures computed by docking method for adsobtion amoxicillin on graphene (AG)

Figure 3: The results of docking for A G

Analysing the findings of molecular orbitals computations

HOMO in chemistry refers to the High Occupied Molecular Orbital, while LUMO indicates Low Unoccupied Molecular Orbital or vacant molecular orbital. In addition, energy gap refers to the energy differentiation between the two orbitals, which is generally characterized by the HLG sign, and formula 2 is applied to measure it. In this formula, EH and EL refer to the HOMO and LUMO orbital energies, respectively. The energy gap is straight associated with the electric conductivity of molecules. Indeed, mixtures with little energy gaps can transfer electrons more effortlessly via the prohibited band into the conduction band. Thus, compounds with smaller energy gaps exhibit greater electric conductivity compared with materials with profound energy gaps. The findings illustrated in Table 2 demonstrate that the level of energy gap after absorbing amoxicillin reduced dramatically on the graphene surface. Indeed, the amoxicillin conductivity enhanced sharply after interplay with graphene. Another noteworthy issue is that we can utilize the rise in conductivity induced by the absorption on the graphene surface to detect and calculate them. Otherwise, graphene is employed to produce novel electrochemical sensors for calculating amoxicillin. The next explored variable is chemical stiffness (η) as the amount drawn from equation 3. Chemical stiffness can be a reliable account of the reactivity of a novel material since molecules with softer structures and lower chemical stiffness can more effortlessly alter their electron thickness. Consequently, the electron transfer, which is necessary for chemical interactions, takes place = easier and more suitable in soft materials. Table 2 demonstrates that the amoxicillin reactivity is enhanced after being absorbed on graphene since all derivatives extracted via interaction with graphene possess less chemical stiffness compared with intact amoxicillin. The value of chemical potential (µ) utilized for obtaining the other variables was also measured based on equation 4. Electrophilicity (ω) and the highest charge transmitted to the system (∆Nmax) are both acceptable features that show the inclination of a material to attract electrons. These two features were measured via equations 5 and 6, respectively. When two molecules start reacting with each other, one molecule plays the role of an electrophile whereas the other functions as a nucleophile. Moreover, a molecule with greater electrophile capacity and charge potential is tended to serve as a receptor electron. However, a compound with little electrophile capacity and charge potential will be more inclined to deliver electrons to the system. According to the results illustrated in the table, amoxicillin is inclined to act as an electron donor in reaction with the nanostructure as its electrophile potential amounts to 0 electron volts. Conversely, intact graphene acts as an electron receiver since its electrophile capacity amounts to 0.01 electron volts. Hence, we can conclude that graphene is able to be involved in electrochemical reactions with amoxicillin. Furthermore, Table 2 reveals that the amoxicillin electrophilicity raised after absorption on the graphene surface. Consequently, we conclude that the amoxicillin tendency towards adsorbing electron was enhanced after reacting with nanostructure. Also, bipolar time of the investigated structures was examined. This feature is a suitable benchmark for measuring the solvency level of molecules in polar solvents. Compounds with greater dipole time will have superior solvency in water, while molecules with less dipole time would have less solvency in polar solvents. As it can be observed, the dipole time of amoxicillin rises after absorption on the graphene surface. Accordingly, graphene derivatives with amoxicillin will manifest greater solvency in water compared with amoxicillin.

Table 2: Energy levels of HOMO and LUMO orbitals, chemical potential, electrophilicity, chemical stiffness, energy gap, highest load transmitted to the system, dipole time for amoxicillin, and its most fixed complex with graphene











Dipole moment


































E HOMO(eV), E LUMO (eV), Gap Energy (eV), ƞ:Hardness(eV), µ: Chemical Potential(eV), S: Softness (eV), ω: Electrophilicity(eV), ΔNmax (eV), Dipole Moment (Debye)

Calculations of computed HOMO and LUMO frontier molecular orbitals show that electron transfer does not occur, so the interaction is of weak van der Waals type (see Figure 4).

Figure 4: a) DFT-computed HOMO and LUMO frontier molecular orbitals and DOS maps for amoxicillin, b) complex of graphene with amoxicillin: AG




η = (EL − EH)/2


µ = (EL + EH)/2


ω = µ2/2η


∆Nmax=− µ/ η


Electrostatic capacity diagrams are exceptionally valuable three-dimensional maps of molecules. They make it possible to visualize the charge dispersions of molecules and charge-related features of molecules. In addition, they lead to observe both dimension and form of molecules. In chemistry, electrostatic potential diagrams are valuable in anticipating activity of complex molecules. To facilitate the data interpretation related to electrostatic potential energy, a colour spectrum, with red and blue as the least and the highest electrostatic potential energy values, respectively, was used to represent the different levels of electrostatic potential energy values. Therefore, red indicates negative charge, while blue represents positive charges. In other words, the red colour with negative charge suggests the lowest electrostatic capacity (i.e. it is loose or additional electrons) and plays the role of an electrophilic attacker. The blue colour demonstrates the highest electrostatic capacity and functions conversely [16].

Due to the shapes of electron surfaces, contour diagrams of the molecular ESP [17] and the positions of the HOMO and LUMO orbitals on the surface of the drug, it appears that at the position of benzene and the oxygen ring on the ring, a higher electron cloud density [5] is located. Therefore, in electrophilic reactions are occurred through this position. on the other hand, the illustrations of more stable structures produced from molecular docking calculations, and the theoretical findings derived from this research are consistent with each other (see Figure 5).

Quantum theory of atoms in molecules (QTAIM)

We applied the AIM analysis to identify the existence of bond critical points (BCPs) of the intramolecular bonds and to calculate their energies, as represented in Table 3. These features are listed in Table 2 for the intramolecular bonds in the explored molecules. A significant relationship was found between the values of (rc) and L(rc). The positive Laplacian  (rc) values in Table 2 demonstrate electronic charge draining along the bond way, which could characterize closed shell interactions bonds. Table 3 illustrates the measured bonds energies as follows: Vc indicates the density of nearby potential electron energy, and Gc denotes the density of nearby kinetic electron energy. Moreover, the Gc/Vc ratio, where Vc denotes the density of nearby potential electron energy, and Gc refers to the density of nearby kinetic electron energy. Finally, we used the Gc/Vc ratio as an indicator of the nature of bonds: for Gc/Vc>1, the bond is noncovalent, while for 0.5<Gc/Vc<1, it is somewhat covalent. In addition, for Gc/Vc<0.5, the bond is covalent (EA., B, kcal.mol-1), and bond critical position data (in a.u.) from quantum theory of atoms in molecules analysis [18-31].

Figure 5: Diagram of a molecule surface or electrostatic potential (ESP) presents incomplete dispersion of alteration along the molecule surface. These diagrams are exceptionally worthwhile and can describe molecular polarity and make the positioning of the dipole arrow possible. The ESP surface is made with red areas and blue areas indicating negative and positive, respectively. Notice: Graphene hydrogens are avoided to better understand the shape

Table 3: Bond critical points (BCPs) of the intramolecular bonds and to calculate their energies for A CG




Ñ2r (rc)



Bond type


























































Quantum theory of atoms in molecules

The most often used criteria of the existence of bonding interactions are the electron density r(rc) and the Laplacian of electron density Ñ2r (rc) at the BCPs. These parameters for the intramolecular bonds in the studied molecules are presented in Table 3. There is a good correlation between the r(rc) and Ñ2r (rc) values. The positive values of Laplacian Ñ2r (rc) in Table 3 indicate depletion of electronic charge along the bond path, which is a characteristic of closed shell interactions bonds. Table 2 lists the bonds energies calculated similar to the following equations:

EA...B  =1/2Vc


Vc  =1/4Ñ2r (rc)- 2Gc    


Where, Vc is the local potential electron energy density and Gc is the local kinetic electron energy density. Finally, the Gc/Vc ratio where, Vc is the local potential electron energy density and Gc is the local kinetic electron energy density. Finally, the Gc/Vc ratio, was used as a criterion of the nature of bonds: for _Gc/Vc>1, the bond is noncovalent, whereas for 0.5<Gc/Vc<1, it is partially covalent and Gc/Vc<0.5 it is covalent. bonds (EA.B, kcal.mol-1), and bond critical point data (in a.u.) from quantum theory of atoms in molecules analysis Equation 7 and 8 (see Figure 6 and 7).

Figure 6: Contour diagrams of the molecular electrostatic potential (ESP) of amoxicillin and AG in its ground state. Contours of places shadowed in dark grey indicate places of negative ESP, while the bright red places indicate the ESP positive state. Calculations were done at the B3LYP/6-31G* level of hypothesis. Electrostatic potential from Total SCF Density (npts = 118, 111, 76; res (A) = 0.176392, 0.176392, and 0.176392), P (0,0,1,0)). Electrostatic potential from Total SCF Density, Electrostatic potential from Total SCF Density (isoval = 0.0004)

Figure 7: Critical points of the intramolecular bonds


Adsorbing Amoxicillin on the Graphene surface was examined by using the molecular docking and density function theories. The computed binding energy indicated the exothermic and spontaneous absorption procedure of this drug and this procedure was measured at room temperature. In addition, the analysis of molecular orbitals suggested that complex Graphene with amoxicillin had more electrophile potential and were more conductive and reactive compared with intact Amoxicillin. It was also found that Graphene could be used to create novel electrochemical sensors to identify and measure Amoxicillin. Graphene is a chemical sensor to detect compounds that structurally bear a close resemblance to Amoxicillin. Based on the measured results, it is recommended that the function of Graphene in removing and measuring Amoxicillin and impact of these nanostructures on their energy qualities be explored empirically. Likewise, since the molecular bond between amoxicillin and graphene is weak, we can look at graphene as a suitable sensor for amoxicillin.


This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.

Authors' contributions

All authors contributed to data analysis, drafting, and revising of the paper and agreed to be responsible for all the aspects of this work.

Conflict of Interest

There are no conflicts of interest in this study.


Roya Ahmadi



Somayeh Pour Karim, Roya Ahmadi, Mohammad Yousefi, khadijeh kalateh, Goldasteh Zaree. Interaction of Graphene with Amoxicillin Antibiotic by in silico study. Chem. Methodol., 2022, 6(11) 861-871


URL: http://www.chemmethod.com/article_154886.html

[1]. Mustafa A., Zahran A.A., Etryl, Pentyl, Hexyl, and Nonyl. Preparation and Explosive Properties, Journal of Chemical and Engineering Data, 1963, 8:135 [Crossref], [Google Scholar], [Publisher]
[2]. Dai Lam T., Van Chat N., Bach V.Q., Loi V.D., Van Anh N., Simultaneous degradation of 2, 4, 6-trinitrophenyl-N-methylnitramine (Tetryl) and hexahydro-1, 3, 5-trinitro-1, 3, 5 triazine (RDX) in polluted wastewater using some advanced oxidation processes, Journal of Industrial and Engineering Chemistry, 2014, 20:1468 [Crossref], [Google Scholar], [Publisher]
[3]. Reddy T.V., Olson G.R., Wiechman B., Reddy G., Torsella J., Daniel F.B., Leach G.J., oxicity of Tetryl (N-Methyl-N, 2, 4, 6–Tetranitroaniline) in F344 Rats, International Journal of Toxicology, 1999, 18:97 [Crossref], [Google Scholar], [Publisher]
[4]. Bai S., Shen X., Graphene–inorganic nanocomposites, RSC Advances, 2012, 2:64 [Crossref], [Google Scholar], [Publisher]
[5]. Stringer R.C., Gangopadhyay S., Grant S.A., Detection of nitroaromatic explosives using a fluorescent-labeled imprinted polymer, Analytical Chemistry, 2010, 82:4015 [Crossref], [Google Scholar], [Publisher]
[6]. Heidary R., simultaneous determination of carbazoles in water samples by cloud point extraction coupled to HPLC, Journal of Applied Chemical Research, 2013, 7:21 [Crossref], [Google Scholar], [Publisher]
[7]. Han G., Gou R.J., Zhang S.H., Wu C.L., Zhu S.F., Theoretical investigation into the influence of molar ratio on binding energy, mechanical property and detonation performance of 1, 3, 5, 7-tetranitro-1, 3, 5, 7-tetrazacyclo octane (HMX)/1-methyl-4, 5-dinitroimidazole (MDNI) cocrystal explosive, Computational and Theoretical Chemistry, 2017, 1109:27 [Crossref], [Google Scholar], [Publisher]
[8]. Wang K., Fu X.L., Tang Q.F., Li H., Shu Y.J., Li J.Q., Pang W.P.,  Theoretical investigations on novel energetic salts composed of 4-nitro-7-(4-nitro-1, 2, 3-triazol-1-olate)-furazano [3, 4-d] pyridazine-based anions and ammonium-based cations, Computational Materials Science, 2018, 146:230 [Crossref], [Google Scholar], [Publisher]
[9]. Allis D.G., Zeitler J.A., Taday P.F., Korter T.M., Theoretical analysis of the solid-state terahertz spectrum of the high explosive RDX, Chemical Physics Letters, 2008, 463:84 [Crossref], [Google Scholar], [Publisher]
[10]. Fayet G., Joubert L., Rotureau P., Adamo C. On the use of descriptors arising from the conceptual density functional theory for the prediction of chemicals explosibility, Chemical Physics Letters, 2009, 467:407 [Crossref], [Google Scholar], [Publisher]
[11]. Ahmadi R., Ebrahimikia M., Calculation of thermodynamic parameters of [2.4. 6] three nitro toluene (TNT) with nanostructures of fullerene and boron nitride nano-cages over different temperatures, using density functional theory, Physical Chemistry Research, 2017, 5:617 [Crossref], [Google Scholar], [Publisher]
[12]. Jalali Sarvestani M.R., Boroushaki T., Ezzati M., The Effect of B12N12 Substitution on the Properties of TEX Energetic Materials in Different Temperature Conditions: A DFT Study, International Journal of New Chemistry, 2018, 5:18 [Crossref], [Google Scholar], [Publisher]
[13]. Jalali Sarvestani M.R., Ahmadi R., Investigating the Complexation of a Recently Synthesized Phenothiazine with Different Metals by Density Functional Theory, International Journal of New Chemistry, 2017, 4:101 [Crossref], [Google Scholar], [Publisher]
[14]. Frisch M.J., Trucks G.W., Schlegel H.B., Scuseria G.E., Robb M.A., Cheeseman J.R., Scalman G., Barone V., Mennucci B., Petersson G.A., Nakatsuji H., Caricato M., Li X., Hratchian H.P., Izmaylov A.F., Bloino J., Zheng G., Sonnenberg J.L., Hada M., Ehara M., Toyota K., Fukuda R., Hasegawa J., Ishida M., T. Nakajima, Honda Y., Kitao O., Nakai H., Vreven T., Montgomery J.A., Jr., Peralta J.E., Ogliaro F., Bearpark M., Heyd J.J., Brothers E., Kudin K.N., Staroverov V.N., Kobayashi R., Normand J., Raghavachari K., Rendell A., Burant J.C., Iyengar S.S., Tomasi J., Cossi M., Rega N., Millam J.M., Klene M., Knox J.E., Cross J.B., Bakken V., Adamo C., Jaramillo J., Gomperts R., Stratmann R.E., Yazyev O., Austin A.J., Cammi R., Pomelli C., Ochterski J.W., Martin R.L., Morokuma K., Zakrzewski V.G., Voth G.A., Salvador P., Dannenberg J.J., Dapprich S., Daniels A.D., Farkas O., Foresman J.B., Ortiz J.V., Cioslowski J., Fox D. J.: Gaussian 09. Revision A.02 ed.; Gaussian, Inc.: Wallingford CT, 2009. [Publisher]
[15]. Becke A.D., Density-functional exchange-energy approximation with correct asymptotic behavior, Physical Review A, 1998, 38:3098 [Crossref], [Google Scholar], [Publisher]
[16]. Lee C., Yang W., Parr R.G., Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density, Physical Review B, 1998, 37:785 [Crossref], [Google Scholar], [Publisher]
[17]. Hazrati M.K., Hadipor N.L., Adsorption behavior of 5-fluorouracil on pristine, B-, Si-, and Al-doped C60 fullerenes: A first-principles study, Physics Letters A, 2016, 380:937 [Crossref], [Google Scholar], [Publisher]
[18]. Vessaly E., Behmagham F., Massoumi B., Hosseinian A., Edjlali L., Carbon nanocone as an electronic sensor for HCl gas: quantum chemical analysis, Vaccum, 2016, 134:40 [Crossref], [Google Scholar], [Publisher]
[19]. Jalali Sarvestani M.R., Gholizadeh Arashti M., Mohasseb M. Quetiapine Adsorption on the Surface of Boron Nitride Nanocage (B12N12): A Computational Study, International Journal of New Chemistry, 2020, 7:87 [Crossref], [Google Scholar], [Publisher]
[20]. Ahmadi R., Jalali Sarvestani M.R., Sadeghi B., Computational study of the fullerene effects on the properties of 16 different drugs: A review, International Journal of Nano Dimension, 2018, 9:325 [Crossref], [Google Scholar], [Publisher]
[21]. Ahmadi R., Jalali Sarvestani M.R., Quarterly Journal of Iranian Chemical Communication, 2019, 7:344 [Google Scholar], [Publisher]
[22]. Jalali Sarvestani M.R., Ahmadi R., Trinitroanisole adsorption on the surface of boron nitride nanocluster (B12N12): A computational study, Journal of Water and Environmental Nanotechnology, 2020, 5:34 [Crossref], [Google Scholar], [Publisher]
[23]. Itodo A., Itodo O., Iornumbe E., Fayomi M.O., Sorptive chelation of metals by inorganic functionalized organic WOx–EDA nanowires: adsorbent characterization and isotherm studies, Progress in Chemical and Biochemical Research, 2019, 1:50 [Crossref], [Google Scholar], [Publisher]
[24]. Adole V.A., Computational Chemistry Approach for the Investigation of Structural, Electronic, Chemical and Quantum Chemical Facets of Twelve Biginelli Adducts, Journal of Applied Organometallic Chemistry, 2021, 1:29 [Crossref], [Google Scholar], [Publisher]
[25]. Duru C.E, Duru, I.A., Computational Investigation of Infrared Vibrational Frequency
Shift Modes in Schiff Base-Transition Metal Complexes, Journal of Applied Organometallic Chemistry, 2022, 2:54 [Crossref], [Google Scholar], [Publisher]
[26]. Farhami, N., A Computational Study of Thiophene Adsorption on Boron Nitride Nanotube, Journal of Applied Organometallic Chemistry, 2022, 2:163 [Crossref], [Google Scholar], [Publisher]
[27]. Oji Moghanlou A., Salimi F., The investigation of antibacterial activity and cell viability of rGO/Cu2O nanocomposite, Asian Journal of Nanoscience and Materials, 2022, 5:63 [Crossref], [Publisher]
[28] Jalali Sarvestani M.R., The effect of doping graphene with silicon on the adsorption of cadmium(II): theoretical investigations, Asian Journal of Nanoscience and Materials, 2020, 3:280 [Crossref], [Google Scholar], [Publisher]
[29] Abedin Dargoush Sh., Hatamie S., Irani Sh., Soliemani M., Hanaee-Ahvaz H., Naderi Sohi A.,  Strontium doped nanohydroxy apatite/reduced graphene oxide nanohybrid is speed up osteogenic differentiation of human mesenchymal stem cells, Asian Journal of Nanoscience and Materials, 2020, 3:226 [Crossref], [Google Scholar], [Publisher]
[30] Jalali Sarvestani M.R., Ahmadi R,. A comprehensive DFT study on the adsorption of tetryl on the surface of graphene, Asian Journal of Green Chemistry, 2020, 4:269 [Crossref], [Google Scholar], [Publisher]
[31] Samadi N., Ansari R., Khodavirdilo B., Synthesis nanoparticles derivations of graphene oxide and poly (Styrene–alternative-maleic anhydride) for removing zinc(II) ions from aqueous solutions, Asian Journal of Green Chemistry, 2019, 3:288 [Crossref], [Google Scholar], [Publisher]