Graphene has a two dimensional crystal without volume and with one surface, that increase the influence of surface dopants and hold new magnetic and electronic properties. In fact, the single atomic plane of graphite is regarded as graphene. The synthesis of garaphene initiated a novel revolution in nanotechnology . The special two-dimensional structure, noticeable mechanical, electrical and thermal properties of graphene exited many investigations [2-5]. Accessible production and economical, a large surface for adsorption of hydrogen are the advantages of graphene. Also, theoretical researches excited attempts for illustration of graphene properties [6-12]. Also, the organometallic chemistry of graphitic compounds has excessive probable for the understanding of novel materials and h6-complexation reactions of chromium with numerous forms of graphene, graphite and carbon nanotubes have been reported .
Cisplatin is used extensively as anticancer drug  which has been predominantly fruitful in treating small cell lung, testicular, ovarian, neck tumors, and head 15. Then the finding of this drug [16, 17] much growth has been made in comprehending its manner of achievement and numerous facts of the mechanism of its antitumor activity are currently perfect [18, 19]. The employment of platinum-based drugs in the clinic has their serious side effects (for example: vomiting. nausea, myelosuppression, neurotoxicity, and ototoxicity) . Consequently, growth of novel delivery systems that can decrease the severity of the drugs’ side-effects is valuable [21, 22]. For instance, a cisplatin slow-release hydrogel drug delivery system based on a formulation of the macrocycle cucurbituril, gelatin and polyvinyl alcohol has been reported . Furthermore, theoretical studied have been reported about the structure, bonding, properties and action mechanism of the drugs [24-28].
In their attempts to progress a drug delivery system, numerous researchers have assessed the capability of carbon nanotubes (CNT) and nanocages . Drug delivery can be achieved either through filling the inner space of the tubes with porphyrins , metals , and biomolecules  or by attaching proteins and compounds (via a specific adsorption or covalent bond) on their external surface. In other research, loading of a phenanthroline-based Platinum(II) complex onto the surface of a carbon nanotube via π–π Stacking has been studied .
Theoretical calculations using the DFT formalism are very general in order to realize thorough electronic structures of complicated coordination compounds from a microscopic viewpoint [26, 34-50]. The current calculations did not contain micro-solvent effect of solvent molecules with solute cations as solvent molecules were not used explicitly.
To realize the solvation effect on the interaction of cis-PtCl2(NH3)2 complex and graphene in solution more quantitatively from the microscopic viewpoint, more sophisticated calculations are essential such as QM/MM (quantum mechanics/molecular mechanics) approach.
The basis goal of this investigation was to reveal the interaction of cis-PtCl2(NH3)2 complex and garaphene by the quantum mechanics method in gas and solvent phases. The stability, dipole moment, structural parameters, of this complex was illustrated Moreover, EDA and QTAIM methods were used for illustration of the interactions in this complex.
All calculations were carried out with the Gaussian 09 suite of program . The calculations of systems contain main group elements described by the standard 6-311G(d,p) basis set [52-55]. Calculation related to Pt element was performed using the element standard Def2-TZVPPD basis set .The pseudo-potential effective core potential (ECP) using the Def2-TZVPPD basis set, was applied to describe Pt . The Modified Perdew-Wang Exchange and Correlation (MPW1PW91) method was employed for the purpose of geometry optimization. Harmonic vibrational frequencies were calculated to verify that the optimization structures have no imaginary frequency.
For the solvation effect study, the structure of complex was reoptimized in selected solvents by a self-consistent reaction field (SCRF) approach, in particular using the polarizable continuum model (PCM) .
The bonding interactions between the graphyne and cis-PtCl2(NH3)2 complex fragment were evaluated considering energy decomposition analysis (EDA) in the Multiwfn 3.3.5 software package . Between the two fragments, the instantaneous interaction energy (Eint) was calculated as:
DEint= DEpolar +DE els+ DE Ex
Where Epolar is the electron density polarization term (the induction term) calculated by subtracting E (SCF last) from E (SCF 1st). Eels and EEx are the electrostatic interaction and the exchange repulsion terms, respectively.
The Multiwfn 3.3.9 software package was also used for the topological analysis of electron density .
RESULTS AND DISCUSSION
Fig. 1 present the optimized structure of cis-Pt(NH3)2Cl2…Graphene complex in gas and solution phases. To evaluate the interaction of cis-Pt(NH3)2Cl2 and graphene, we have selected coronene cluster (C24H12) as a large planar polycyclic aromatic hydrocarbons (PAH) that is accepted to be a proper model for theoretical examination on graphene .It can be observed, square planar of cis-Pt(NH3)2Cl2 complex is oriented approximately parallel with graphene layer.
The absolute energy and solvation energy values of the studied complex in gas and various solvent are gathered in Table 1. The negative value of solvation energy values reveals the more stability of cis-Pt(NH3)2Cl2…Graphene complex in solution in compared to gas phase. The larger negative values of solvation energy values in more polar solvent illustrates the more stability of this complex in the more polar solvent. It can be observed a good linear correlation between solvation energy values and dielectric constant values of solvents (Fig. 2):
Esolv = -0.8799 e - 14.908; R² = 0.9648
Dipole moment values of the of cis-Pt(NH3)2Cl2…Graphene complex are listed in Table 1. The larger dipole moment value of various solvents attribute to the additional dipole moment induced by the solvents. There is a good linear relationship between solvation energy values and dipole moment in different solvents (Fig. 3):
Esolv = -4.3868 m + 51.641; R² = 0.9163
Energy decomposition analysis (EDA)
The chemical bond nature between the graphene fragment with cis-Pt(NH3)2Cl2 complex is explored using an energy decomposition analysis (EDA). Results of these calculations reveal the interaction energy value of graphene fragment with cis-Pt(NH3)2Cl2 complex is -7.57 kcal/mol, in gas. The interaction energy values in different solvent are listed in Table 1.The negative total energy of interaction values between graphene fragment with cis-Pt(NH3)2Cl2 complex show interactions between two fragments. It can be seen, the interaction strength between two fragments is stronger in gas phase in compared to solution phase.
Fig. 3 presents the shortest distances of Pt … C(Graphene) and H(NH3) …. C(Garaphene) in the cis-Pt(NH3)2Cl2…Graphene complex in gas. These values for this complex in various solvents are gathered in Table 1. It can be seen, these distances are longer in solution phase than gas phase. As a result, these distances predict stronger interaction between cis-Pt(NH3)2Cl2 and graphene in gas phase. These distances are compatible with interaction energy values of the complex in two phases.
Molecular orbital analysis
Frontier orbital energy values and HOMO-LUMO gap of cis-Pt(NH3)2Cl2…Graphene complex in the and selected solvents are shown in Table 2. The HOMO and LUMO energy values of the cis-Pt(NH3)2Cl2…Graphene complex are -6.22 and -2.05 eV in gas. The comparison of these values with corresponding values in solvent system show that HOMO are stabilized in solution phase. In contrast, the stability of LUMO decrease in solution phase. The HOMO-LUMO gap values in the studied complex are 4.17 eV in gas. It can be observed HOMO-LUMO gap increase in solution phase. This means that electron transfer between frontier orbitals (HOMO→LUMO) is easy in gas phase in compared to solution phase. The increase in the dielectric constant of the solvent and possible interactions between the solute and solvent molecules causes to increase the energy gap.
Plots of frontier orbital in this complex are presented in Fig. 4. It can be observed, only cis-Pt(NH3)2Cl2 fragment contributes in HOMO of complex. In contrast, only graphene fragment contributes in LUMO complex.
The absolute free energy and enthalpy values of the cis-Pt(NH3)2Cl2…Graphene are summarized in Table 4. The solvation free energy and enthalpy values are evaluated as:
ΔXsolvation = Xsolv -Xgas; X = G, H
It can be found, increasing the polarity of solvent resulted in decreased Gsolv and Hsolv values of
the investigated complex. Additionally, Gsolv and Hsolv values show good relationships with the dielectric constant values:
Gsolv = -0.6243 ε-15.945; R2 = 0.9966
Hsolv = -0.8496 ε - 14.545; R² = 0.9682
Quantum theory of atoms in molecules (QTAIM) analysis
Understanding of the chemical and physical properties of molecules is possible by using Quantum theory of atoms in molecules (QTAIM) analysis [63-67]. Therefore, we study topological parameters of the studied complexes. QTAIM analysis reveals bond critical point (BCP) between Pt and H(NH3) with Cgraphyne atoms. The results of these analysis for Pt-Cgaraphyne and H-C bonds are gathered in Table 4. It can be found, the total electron density is larger on the BCP of the H-C bonds in compared to Pt-C bond.
It is obvious that Pt-C bonds are polar, as is the situation regularly in the coordinate bonds. This is detected in the positive Laplacian values, which illustrates charge depletion in the bond critical points. Laplacian of electron density of in bond critical point (BCP) of Pt-Cgraphyne and H-C bonds bonds in the basis of the QTAIM calculations are listed in Table 4. As seen, values are positive at their corresponding BCP. These values are expected for closed-shell interactions.
Energy density of in bond critical point (BCP) of Pt-C(graphene) and H-C interactions in the basis of the QTAIM calculations are reported in Table 4 the total energy density (H) is defined as:
H= G + V
In this equation G and V are Lagrangian kinetic energy and Virial energy density, respectively. QTAIM results reveal V values are negative. This negative value of H is compatible with shared interactions in Pt-Cgraphyne and H-C bonds bonds. Consequently, these results are compatible with preceding findings for the M–C bonds in organometallic complexes  and transition metal carbonyl clusters , where the characteristics of the metal–ligand bonding is intermediate between closed-shell and shared interactions. On the other hand, the positive (r) and H(r) topological descriptors and -G(r)/V(r) descriptor is also greater than one at BCP of C-H, implying that this interaction should be classified as weak and non-covalent.
Theoretical study on interaction of cis-Pt(NH3)2Cl2 with graphene in gas phase and selected solvents was revealed :
1. Negative value solvation energy of the between cis-Pt(NH3)2Cl2 …. graphene complex signifies the more stability of the complex in solution in comparison to gas phase.
2. There are good linear correlations between Esol with dipole moment and dielectric constant values.
3. EDA and structural parameters predict stronger interaction between cis-Pt(NH3)2Cl2 and graphene in gas phase.
4. cis-Pt(NH3)2Cl2 and graphene fragments contribute in HOMO and LUMO of complex, respectively.
5. In the basis of QTAIM analysis, the characteristics of the Pt…C(Graphene) and H(NH3)… C(Graphene) interactions was intermediate between closed-shell and shared interactions.
CONFLICT OF INTEREST
The authors declare that there is no conflict of interests regarding the publication of this manuscript.