Solvent influence on the interaction of cis-PtCl2(NH3)2 complex and graphene: A theoretical study
Reza Fazaeli
1
(
Department of Chemistry, South Tehran Branch, Islamic Azad University, Tehran, IRAN
)
E Ebrahimi Mokarram
2
(
Department of Chemistry, Science and Research Branch, Islamic Azad University, Tehran, IRAN
)
H Aghaei
3
(
Department of Chemistry, Science and Research Branch, Islamic Azad University, Tehran, IRAN
)
K Zare
4
(
Department of Chemistry, Science and Research Branch, Islamic Azad University, Tehran, IRAN
)
Mohammad Yousefi
5
(
Department of Chemistry, Science and Research Branch, Islamic Azad University, Tehran, Iran
)
کلید واژه: solvent effect, Graphene, Cis-PtCl2(NH3)2 Complex, energy decomposition analysis (EDA), polarizable continuum model (PCM),
چکیده مقاله :
In this investigation the interaction of cis-PtCl2(NH3)2 complex and graphene were investigated with MPW1PW91 method in gas and solvent phases. The solvent effect were examined by the self-consistent reaction field theory (SCRF) based on Polarizable Continuum Model (PCM). The selected solvents were chloroform, chlorobenzene, bromoethane, dimethyldisulfide, and dichloroethane. The solvent effect on the frontier orbital energy and HOMO-LUMO gap were studied. Negative value solvation energy of the between cis-Pt(NH3)2Cl2 …. graphene complex signifies the more stability of the complex in solution in compared to gas phase.The characterization of the interaction between two fragments was clarified with energy decomposition analysis (EDA). Pt-C(Graphene) and H(NH3)...C(Garaphen) interactions in the graphene … cis-PtCl2(NH3)2 complex were analyzed using quantum theory of atoms in molecules analysis (QTAIM). cis-Pt(NH3)2Cl2 and graphene fragments contribute in HOMO and LUMO of complex, respectively.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.EDA and structural parameters predict stronger interaction between cis-Pt(NH3)2Cl2 and graphene in gas phase.
چکیده انگلیسی :
Solvent influence on the interaction of cis-PtCl2(NH3)2 complex and graphene: A theoretical study
Elham Ebrahimi Mokarram1, Reza Fazaeli2,*, Hossein Aghaei1, Mohammad Yousefi1, Karim Zare1
1 Department of Chemistry, Science and Research Branch, Islamic Azad University, Tehran, IRAN
2 Department of Chemistry, South Tehran Branch, Islamic Azad University, Tehran, IRAN
Email: el1985e@gmail.com
Abstract:
In this investigation the interaction of cis-PtCl2(NH3)2 complex and graphene were investigated with MPW1PW91 method in gas and solvent phases. The solvent effect were examined by the self-consistent reaction field theory (SCRF) based on Polarizable Continuum Model (PCM). The selected solvents were chloroform, chlorobenzene, bromoethane, dimethyldisulfide, and dichloroethane. The solvent effect on the frontier orbital energy and HOMO-LUMO gap were studied. The characterization of the interaction between two fragments was clarified with energy decomposition analysis (EDA). Pt-C(Graphene) and H(NH3)...C(Garaphen) interactions in the graphene … cis-PtCl2(NH3)2 complex were analyzed using quantum theory of atoms in molecules analysis (QTAIM).
Keywords: Graphene, cis-PtCl2(NH3)2 complex, , solvent effect, energy decomposition analysis (EDA), polarizable continuum model (PCM), Quantum theory of atoms in molecules analysis (QTAIM).
Introduction:
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, single atomic plane of graphite is regarded as graphene. The synthesis of garaphene initiated a novel revolution in nanotechnology 1. 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 has been reported 13.
Cisplatin is used extensively as anticancer drug 14 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) 20. 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 cucurbit[7]uril, gelatin and polyvinyl alcohol has been reported 23. 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 29. Drug delivery can be achieved either through filling the inner space of the tubes with porphyrins 30, metals 31, and biomolecules 32 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 33.
Theoretical calculations using the DFT formalism are very general in order to realize thorough electronic structures of complicated coordination compounds from a microscopic viewpoint 27, 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 were illustrated Moreover, EDA and QTAIM methods were used for illustration of the interactions in this complex.
Computational method
All calculations were carried out with the Gaussian 09 suite of program 51. 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 were performed using the element standard Def2-TZVPPD basis set 56. The pseudo-potential effective core potential (ECP) using the Def2-TZVPPD basis set, was applied to described Pt 57. The Modified Perdew-Wang Exchange and Correlation (MPW1PW91) method was employed for the purpose of geometry optimization58. 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) 59.
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 60. 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 61 .
Results and discussion
1. Energetic aspect
Figure 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 62.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 (Figure 2):
Esolv = -0.8799 e - 14.908; R² = 0.9648
2. Dipole moment
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 (Figure 3):
Esolv = -4.3868 m + 51.641; R² = 0.9163
3. 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.
4. Structural parameters.
Figure 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 distance 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.
5. 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 correspond values in solvent system show that HOMO are stabilized in solution phase. In contrast, the stability of LUMO is decreases 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 these complex are presented in Figure 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.
6. Thermodynamic parameters
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
7. 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 reveal 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, Ñ2r 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 values 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 68 and transition metal carbonyl clusters 69, where the characteristics of the metal–ligand bonding is intermediate between closed-shell and shared interactions. On the other hand, the positive Ñ2r(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.
Conclusion:
Theoretical study on interaction of cis-Pt(NH3)2Cl2 with graphene in gas phase and five 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 compared 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.
References:
2. Majidi, R.; Karami, A. R., Physica E 2014, 59, 169.
4. Geim, A., Science 2009, 324, 1530.
5. Jalili, S.; Majidi, R., J Comp Theor Nanosci 2006, 3, 664.
6. Majidi, R., Theor Chem Acc 2017, 136, 109.
8. Yaraghi, A.; Ozkendir, O. M.; Mirzaei, M., Superlattices and Microstructures 2015, 85, 784–788.
9. Saha, B.; Bhattacharyya, P. K., Computational and Theoretical Chemistry 2016, 106, 45-51.
13. Sarkar, S.; Niyogi, S.; Bekyarova, E.; Haddon, R. C., Chem. Sci. 2011, 2, 1326.
15. Go, R. S.; Adjei, A. A., J. Clin. Oncol 1999, 17, 409-422.
16. Rosenberg, B.; Camp, L. V.; Krigas, T., Nature 1965, 205, 698.
17. Rosenberg, B.; Camp, L. V.; Trosko, J. E.; H.Mansour, V., Nature 1969, 222, 385.
18. Jamieson, E. R.; Lippard, S. J., Chem. Rev. 1999, 99, 2467.
19. Fuertes, M. A.; Alonso, C.; Perez, J. M., Chem. Rev. 2003, 103, 645.
20. Tsang, R. Y.; Al-Fayea, T.; Au, H.-J., Drug Saf. 2009, 323, 1109–1122.
21. Barry, N. P. E.; Sadler, P. J., ACS Nano 2013, 7, 5654–5659.
22. Gabano, E.; Ravera, M.; Osella, D., Curr. Med. Chem. 2009, 34, 4544–4580.
23. Oun, R.; Plumb, J. A.; Wheate, N. J., Journal of Inorganic Biochemistry 2014, 134, 100–105.
24. Valizadeh, A.; Ghiasi, R., J. Struct. Chem 2017, 58, 1307-1317.
25. Ghiasi, R.; Sadeghi, N., J. Chin. Chem. Soc. 2017, 64, 934–939.
26. Ghiasi, R.; Sadeghi, N., J. Struc. Chem. 2017, in press.
27. Rezazadeh, M.; Ghiasi, R.; Jamehbozorgi, S., J. Struc. Chem 2018, 59, 245-251.
28. Sarraf, S.; R.Ghiasi, Structural Chemistry 2018, 29, 435–440.
29. Ghanbari, H.; Cousins, B. G.; Seifalian, A. M., Macromol Rapid Commun. 2011, 32, 1032–1046.
30. Khlobystov, A. N.; Britz, D. A.; Briggs, G. A. D., Acc. Chem. Res. 2005, 38, 901-909.
32. Yanagi, K.; Miyata, Y.; Kataura, H., Adv. Mater. 2006, 18, 437-441
34. Ghiasi, R.; F.Zafarniya; Ketabi, S., Russian Journal of Inorganic Chemistry 2017, 62, 1371-1378.
35. Alavi, H.; Ghiasi, R., J. Struc. Chem 2017, 58, 30-37.
36. Zafarniya, F.; R.Ghiasi; Jameh-Bozorghi, S., Physics and Chemistry of liquids 2017, 55, 444-456.
37. Zafarnia, F.; Ghiasi, R.; Jamehbozorgi, S., J. Struc. Chem 2017, 58, 1324-1331.
39. Ghiasi, R.; Peikari, A., Physical and Chemistry of Liquids 2017, 55, 421-431.
40. Ghiasi, R.; Peikari, A., Russian Journal of Physical Chemistry A 2016, 90, 2211-2216.
41. Ghiasi, R.; Peikari, A., Journal of Applied Spectroscopy 2017, 84, 148-155.
42. Ghiasi, R.; Pasdar, H.; Fereidoni, S., Russian Journal of Inorganic Chemistry 2016, 61, 327–333.
43. Ghiasi, R.; Nemati, M.; H.Hakimioun, A., J. Chil. Chem. Soc 2016, 61, 2921-2928.
44. Peikari, A.; Ghiasi, R.; Pasdar, H., Russian Journal of Physical Chemistry A 2015, 89, 250–255.
45. Ghiasi, R.; Amini, E., Journal of Structural Chemistry 2015, 56, 1483-1494.
46. Fashami, M. Z.; Ghiasi, R., Journal of Structural Chemistry 2015, 56, 1474-1482.
47. Rezazadeh, M.; Ghiasi, R.; Jamehbozorgi, S., Journal of Applied Spectroscopy 2018, 85, 926-934.
49. Rahimi, M.; Ghiasi, R., Journal Molecular Liquid 2018, 265, 164-171.
50. Ghiasi, R., J. Mol. Liqu 2018, 264, 616-620.
52. Krishnan, R.; Binkley, J. S.; Seeger, R.; Pople, J. A., J. Chem. Phys. 1980, 72, 650-654.
53. Wachters, A. J. H., J. Chem. Phys. 1970, 52, 1033.
54. Hay, P. J., J. Chem. Phys. 1977, 66, 4377-4384.
55. McLean, A. D.; Chandler, G. S., J. Chem. Phys. 1980, 72, 5639-5648.
56. Rappoport, D.; Furche, F., J. Chem. Phys. 2010, 133, 134105.
57. Andrae, D.; Haeussermann, U.; M.Dolg; H.Stoll; H.Preuss, Theor.Chim.Acta 1990, 77, 123-141.
58. Perdew, J. P.; Burke, K.; Ernzerhof, M., Phys. Rev. Lett. 1996, 77, 3865-3868.
59. Tomasi, J.; Mennucci, B.; Cammi, R., Chem. Rev. 2005, 105, 2999.
60. Lu, T.; Chen, F., J. Mol.Graphics. Model 2012, 38, 314.
61. Lu, T.; Chen, F., J. Comp. Chem. 2012, 33, 580-592.
62. Duley, W. W.; Seahra, S. S., Astrophys. J 1999, 522, L129.
63. Sobczyk, L.; Grabowski, S. J.; Krygowski, T. M., Chem. Rev. 2005, 105, 3513–3560.
64. Bader, R. F. W.; Matta, C. F.; Corte´s-Guzman, F., Organometallics 2004, 23, 6253–6263.
65. Fradera, X.; Austen, M. A.; Bader, R. F. W., J. Phys. Chem. A 1999, 103, 304–314.
66. Bader, R. F. W.; Fang, D.-F., J. Chem. Theor. Comput. 2005, 1.
67. Mitrasinovic, P. M., Can. J. Chem. 2003, 81, 542.
68. Palusiak, M., J. Organometallic. Chem 2005, 692, 3866-3873
69. Macchi, P.; Sironi, A., Coordination Chemistry Reviews 2003, 239, 383.