Next Article in Journal
Single-Pixel Photon-Counting Imaging Based on Dual-Comb Interferometry
Next Article in Special Issue
General Molecular Dynamics Approach to Understand the Mechanical Anisotropy of Monocrystalline Silicon under the Nanoscale Effects of Point Defect
Previous Article in Journal
Efficient Visible-Light Photocatalysis of TiO2-δ Nanobelts Utilizing Self-Induced Defects and Carbon Doping
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Interaction between Graphene-Based Materials and Small Ag, Cu, and CuO Clusters: A Molecular Dynamics Study

by
Isabel Lado-Touriño
* and
Alicia Páez-Pavón
School of Architecture, Engineering and Design, Universidad Europea de Madrid, 28670 Villaviciosa de Odón, Spain
*
Author to whom correspondence should be addressed.
Nanomaterials 2021, 11(6), 1378; https://doi.org/10.3390/nano11061378
Submission received: 20 April 2021 / Revised: 16 May 2021 / Accepted: 21 May 2021 / Published: 23 May 2021
(This article belongs to the Special Issue Nanomaterials Investigation by Molecular Dynamics Simulation)

Abstract

:
The excessive use of antibiotics has contributed to the rise in antibiotic-resistant bacteria, and thus, new antibacterial compounds must be developed. Composite materials based on graphene and its derivatives doped with metallic and metallic oxide nanoparticles, particularly Ag, Cu, and Cu oxides, hold great promise. These materials are often modified with polyethylene glycol (PEG) to improve their pharmacokinetic behavior and their solubility in biological media. In this work, we performed molecular dynamics (MD) simulations to study the interaction between small Ag, Cu, and CuO clusters and several graphene-based materials. These materials include pristine graphene (PG) and pristine graphene nanoplatelets (PGN) as well as PEGylated graphene oxide (GO_PEG) and PEGylated graphene oxide nanoplatelets (GO-PEG_N). We calculated the adsorption energies, mean equilibrium distances between the nanoparticles and graphene surfaces, and mean square displacement (MSD) of the nanoclusters. The results show that PEGylation favors the adsorption of the clusters on the graphene surfaces, causing an increase in adsorption energies and a decrease in both distances and MSD values. The strengthening of the interaction could be crucial to obtain effective antibacterial compounds.

Graphical Abstract

1. Introduction

The excessive use of antibiotics has contributed to the rise in antibiotic-resistant bacteria, and new antibacterial compounds must be developed. Graphene is a 2D monolayer material with long edges and large accessible surfaces, which give it an excellent ability to immobilize diverse molecules including nanoparticles or drugs. Hu et al. [1] first reported the antibacterial activity of graphene oxide (GO) and reduced graphene oxide (rGO) nanosheets against E. coli bacteria. Many subsequent studies showed that pristine graphene, GO, and rGO have all antibacterial activity [2,3,4,5,6,7,8,9]. As graphene-based materials tend to aggregate due to strong van der Waals interactions, thus reducing their effective surface and antibacterial activity, they are modified with metals, oxides, polymers, or a combination of these [10,11,12,13,14,15,16]. In particular, the use of Ag, Cu, and their oxides adsorbed on the graphene surface appears promising [17,18,19,20,21,22,23,24,25]. Moreover, graphene-based composites, prepared by incorporating both nanoparticles and polymers, exhibit synergistic antibacterial effects [16,26]. The polymers used include polyethylene glycol (PEG), which is a polymer that helps to obtain stable graphene−PEG materials dispersible in most biocompatible solvents [27,28]. For instance, PEG-GO remains highly dispersed in serum solutions [29] and shows good pharmacokinetic behavior [30] and biocompatibility [31,32].
Molecular dynamics (MD) simulations allow the study of materials from a molecular viewpoint, providing information on an atomic scale usually inaccessible to experimental techniques [33]. The interaction between graphene-based materials and different nanoparticles and polymers has been described in many research works by MD [34,35,36,37,38,39,40,41,42]. Several other authors have performed molecular modeling studies on the interaction between graphene and PEG [43,44,45,46,47,48].
In this work, we carried out MD simulations of the adsorption of small Ag, Cu, and CuO clusters on pristine graphene monolayers (PG), PG nanoplatelets (PGN), PEGylated graphene oxide monolayers (GO_PEG), and PEGylated graphene oxide nanoplatelets (GO-PEG_N). The work focuses on a comparison of the behavior of these three nanoclusters interacting with the above-mentioned graphene-based materials and allows completing our previous results on Ag and Cu nanoclusters interacting with a graphene surface [35]. In our new simulations, much longer simulation lengths are used, and more systems have been included in the calculations: PGN, GO-PEG_N, and CuO clusters. We calculated adsorption energies, mean equilibrium distances between the nanoparticles and graphene surfaces, and mean square displacement (MSD) of the nanoclusters. The results show that PEGylation favors the adsorption of the clusters on the graphene surfaces, causing an increase in adsorption energies and a decrease in both distances and MSD values. Thus, this work can provide important clues on the way that graphene-based materials can be modified to become effective antibacterial compounds.

2. Models and Calculation Method

MD simulations were done with the Forcite module of the Materials Studio 9 software [49] in the NVT ensemble (constant number of particles N, constant volume V, and constant temperature T) at 298 K during 10 ns of simulation length. The interactions between atoms were calculated using the COMPASS forcefield, which is a force field parameterized using experimental and ab initio results that enables an accurate prediction of various gas-phase and condensed-phase properties of several organic and inorganic materials [50], including graphene and graphene-based materials [44,51,52]. The temperature was kept constant using a Nosé-Hoover thermostat. [53,54]. Electrostatic and van der Waals interactions were calculated using an atom-based summation method with a cut-off of 12.5 Å.
The starting configurations used for the production trajectories are shown in Figure 1 (graphene monolayers) and Figure 2 (nanoplatelets). A basal graphene sheet with a width of 46.5 Å and a length of 56.2 Å was used for all structures. To build these starting structures, all the systems were subjected to an annealing procedure composed of NVT MD simulations and further steepest descent and conjugate gradient energy minimization cycles at each temperature. This annealing procedure allows exploring the conformational space for low-energy structures, which will be used as starting configurations, by periodically increasing and then decreasing the temperature (from 300 up to 1000 K) of a classical dynamics trajectory, to avoid trapping the structure in local energy minima. For the sake of comparison, PG, PGN, GO_PEG, and GO-PEG_N were constructed. The PEGylated structures were functionalized with eight short chains of PEG (degree of polymerization n = 10). In the interest of saving computational time, only two layers of graphene were used to model the nanoplatelets. In the PEGylated nanoplatelets, both layers had the same grafting density, i.e., the number of PEG chains attached to the edges of the individual GO layers. In these starting structures, the mean final distances between the basal planes of graphene were 3.7 and 6.9 Å for PGN and GO-PEG_N, respectively. The clusters used to represent the metal nanoparticles were composed of thirteen atoms arranged in an icosahedral shape. This number was chosen because of the well-known stability of this regular icosahedron geometry for Cu and Ag nanoclusters [55,56]. For CuO nanoparticles, a Cu6O6 cluster consisting of two hexagonal layers, one above the other, was chosen. This structure was found to be stable by a density functional theory calculation [57]. No constraints were imposed on the systems, so that all molecules could move freely over the whole simulation length. At the beginning of the calculation, the clusters were placed far from the surfaces, at distances of about 8 Å, as can be seen in the lower part of both Figures, and they were left to evolve with time.
The adsorption energy of the clusters on the surface was calculated from the following equation:
E = EGBM+cluster − (EGBM + Ecluster)
where EGBM+cluster is the mean equilibrium potential energy of the graphene-based material (GBM = graphene-based materials = monolayers or nanoplatelets) interacting with the cluster, and EGBM and Ecluster are the mean equilibrium potential energies of the isolated graphene-based materials and clusters, respectively (all the energies were calculated as average values from three independent 10 ns simulations. In turn, the result of each individual simulation is the average of the energy values obtained during the last 5 ns of simulation). As the energy values are negative, the more negative the adsorption energy is, the stronger the interaction. Previous studies have shown that this adsorption energy could be a good estimate of the binding strengths of adsorbate–adsorbent systems [58].
The mean square displacement (MSD) of the clusters, which measures the spatial extent of clusters random motion, was calculated from their positions over simulation length according to:
MSD = Σ < (ri(0) − ri(t))2 > /N
where ri(0) is the reference position of the cluster, ri(t) is its position at time t, and N is the number of atoms of the cluster. Thus, MSD was calculated from the particle positions obtained from the trajectory using a time step of 5 ps.
MSD can be related to time through the following equation:
MSD = Kαtα
where Kα is the generalized diffusion coefficient and α is the diffusion exponent [59,60]. The value of this exponent is related to the mechanism of diffusion.

3. Results and Discussion

3.1. Adsorption Energies

The adsorption energy values calculated from Equation (1) are presented in Table 1.
From the values listed in Table 1, it can be concluded that the strength of the interaction depends on both the chemical nature of the cluster and the functionalization of the graphene-based material. All the clusters interact more weakly with pristine materials than with functionalized materials in both monolayers and nanoplatelets. The Cu6O6 cluster is the most attracted by the surface and the Cu13 cluster is the least attracted. Ag13 lies in between. The introduction of PEG chains favors the adsorption process, increasing the adsorption energy. In addition, the nanoplatelets seem to increase the adsorption energy in all cases. This last trend should be confirmed by increasing the number of layers of the nanoplatelets in future studies.
Here, we must note the qualitative nature of these results. It is well known that cluster geometry and size influence their adsorption energy and reactivity on graphene surfaces, and several authors have shown that quantum mechanical calculations are needed to fully understand and quantify this kind of interaction [61]. Even when using density functional theory calculations, the results vary widely (adsorption energies ranging from a few kcal/mol to a hundred kcal/mol and cluster–graphene surface distances ranging from 2 to 3 Å) depending on the molecular models, the calculation method, or the number of atoms included in the nanocluster [62,63,64]. What is clear is that metal d-orbitals are involved in charge transfer processes between the clusters and the graphene surface, although this charge transfer is weak in some cases, and the bonding can be mainly attributed to van der Waals forces. Cu and Ag are examples of that, as it has been demonstrated by numerous calculations [65,66]. On the other hand, experimental results are diverse and show that the structure of the material after the interaction depends on many factors, including the initial graphene form and the methods used for its synthesis [67]. In addition, due to limited computational resources, our simulations were performed in vacuum, and we are aware that this may influence the results. Thus, the small clusters and the graphene model used in this work represent a first step to compare qualitatively the three types of particles, and we realize that they cannot capture the whole complexity of real systems.
Cu6O6 interacts strongly with both pristine and functionalized surfaces, although the adsorption process is improved by the presence of PEG chains, as already observed for Cu and Ag clusters. We think the atomic charges of the cluster increase the electrostatic attraction between the cluster and surface, although the exact nature of the interaction should be verified by more accurate quantum mechanical calculations. Sun et al. [68] determined that during the process of nucleation of graphene on Cu-based substrates, the adsorption energies of C were larger for copper oxides than for metallic Cu. Ko et al. [69] prepared CuO nanoparticles covered with a monolayer graphene shell, in which strong C-O-Cu links were formed. Using density functional theory calculations, Mohammadi-Manesh et al. [70] calculated the binding energy of different configurations of adsorption of Cu and CuO on graphene surfaces and found that the interaction between Cu and graphene is physical and the interaction between CuO and graphene is chemical.
Even if the calculation method used in this work cannot describe the formation of chemical bonds between the cluster and the graphene surfaces, the results obtained agree with those obtained from more accurate quantum chemical calculations.

3.2. Cluster–Surface Distances

The mean equilibrium distances between the cluster and nearest surface atoms are shown in Table 2. The reduction of the distances used in the starting structures (8 Å) indicates that the clusters approached the graphene surface as the simulation progressed and attained an equilibrium value toward the end of the calculation. The values listed in Table 2 show that the mean distance between the clusters and the pristine surfaces is about 2.9 Å for both monolayers and nanoplatelets. The introduction of PEG chains brings the cluster closer to the graphene surface leading to reduced distances of about 2.4–2.5 Å. Cluster–pristine material distances are larger than cluster–functionalized material distances, which correlates well with improved adsorption energies for GO_PEG-containing materials.
The distances shown in Table 2 for the pristine systems are similar to those found in previous quantum mechanical [65,66] and MD [71,72] simulations of the interaction between graphene and Cu and Ag. DFT results give distances of 2.3–3.9 and 2.8 Å for Cu and Ag, respectively. The distance found by MD between Ag and graphene is 3 Å. MD results for Cu show a range of binding distances varying between 2 and 3.4 Å, depending on the orientation of Cu relative to graphene. They are related to weak bonding (physical adsorption) attributed to the full occupancy of metal d orbitals [73,74].
The radial distribution functions (RDF) of clusters and graphene-based surfaces are shown in Figure 3. These plots give us an idea of the spread of distances between the clusters and surfaces. All plots present maxima between 2.5 and 5–6 Å, and there are no large differences between the different models. It must be noted that the mean equilibrium distances shown in Table 2 represent distances between the cluster and nearest surface atoms, which must not necessarily coincide with maxima in Figure 3. In the case of the Cu6O6 cluster, the plots corresponding to the Cu atom (shown in red) are slight shifted to larger values, which is in accordance with the results shown in Table 2. Figure 3 also shows that in the PEGylated systems, distances are somewhat smaller.
Figure 4 and Figure 5 show the representative positions of the clusters on the monolayer surface (Figure 4) and the nanoplatelet surface (Figure 5). Only the result obtained from one of the simulations is depicted here. In the case of the GO_PEG surfaces, it was observed that the clusters could stick to any lateral PEG chain of the layers, and once stuck, they remained bonded to them at a distance of about 2.5–2.6 Å in all cases (see Table 2). In the case of CuO clusters, this distance corresponds to O-PEG distances. In all cases, the Cu atom in the cluster was located at somewhat larger distances. PEG is frequently used to prepare metal and metal oxide nanoparticles in solution [75,76,77], as it is well known that PEG molecules strongly adsorb on metal nanoparticles surfaces by coordination through the ether bond of PEG and prevent their aggregation [78]. The capping capability of PEG could explain the strong attraction it exerts on the clusters and their subsequent immobilization on the graphene surfaces.

3.3. Mean Square Displacement

Although a detailed analysis of the mechanisms of cluster diffusion is beyond the scope of this work, we think that the MSD of the cluster center of mass can help shed more light on the differences observed in adsorption energies. Total MSD is plotted versus simulation length in Figure 6 and Figure 7 for monolayers and nanoplatelets, respectively. Larger MSD values indicate greater cluster mobility.
Cu13 and Ag13 clusters interacting with the pristine materials present the highest MSD values. Both PG and PGN weakly attract both clusters, which keep constantly moving over the material surface. However, the functionalization of the graphene layers drastically reduces their mobility (see bottom plots in Figure 6 and Figure 7). The PEG chains act by trapping the clusters and blocking their displacement over the surface. The Cu6O6 cluster shows a completely different behavior as its MSD values are low (blue line plots) in both the pristine and functionalized materials.
To show more clearly the cluster movement, the diffusion trajectories of the cluster center of mass over the graphene surface during the last 250 ps of simulation length for both PG and PGN are displayed in Figure 8 and Figure 9, respectively. Some video representative of the different cluster dynamics is included in the Supplementary Materials (Videos S1, S2, and S3).
The MSD values in all directions are shown in Figure 10 and Figure 11 for monolayers and nanoplatelets, respectively. The graphene surfaces are oriented parallel to the xy plane. As can be inferred from these figures, the MSD in the z direction is approximately zero (green line), as once the cluster approaches the basal plane, no vertical displacement from the surface is observed. MSDxx and MSDyy components have non-zero values due to the movement of the cluster over the surface and parallel to it during the whole simulation length. In the functionalized structures, the cluster may carry out slight movements between or along PEG chains (see Figure 12), although it remains bonded to them. Therefore, MSDxx and MSDyy present low non-zero values. An example of this dynamic behavior is included in the Supplementary Material (Video S4).
The MSD values corroborate the trends followed by the adsorption energies: the more dynamic the cluster is, the lower the adsorption energy. In agreement with our results, Gervilla et al. [79,80] studied the diffusion of Cu and Ag adatoms and clusters on graphene, using ab initio and classical molecular dynamics simulations. They found clusters diffused by a super-diffusive mechanism (α > 1 in Equation (3)), which is characterized by a continuous displacement of the clusters over the graphene surface without being trapped by adsorption sites. They attributed this behavior to a flat potential energy landscape on the surface, which facilitated the cluster diffusion. Manade et al. [74] estimated the diffusion energy barrier of Cu and Ag adatoms on graphene and arrived at the same conclusion.
Due to their higher adsorption energy, the MSD values of the Cu6O6 cluster are low on both the pristine and functionalized surfaces. Once it was attracted by the surface, it remained stuck to it (see Video S2).

4. Conclusions

MD simulations were carried out to study the interaction between Cu13, Ag13, and Cu6O6 clusters on pristine and functionalized graphene-based materials with the goal of finding the best structures for adsorption. Adsorption energies, cluster–surface distances, and MSD values show that both the chemical nature of the cluster and the PEGylation of the surface are critical to strengthen the interaction between the clusters and surfaces. The introduction of PEG chains favors the adsorption of the three clusters by acting as trapping sites, thus reducing the distances between the clusters and graphene surfaces and increasing the adsorption energy. MSD values also point in this direction. Pristine materials weakly attract Cu13 and Ag13 clusters, which keep constantly moving over them. PEG chains act by blocking the displacement of the nanoparticles, which are attracted by their functional groups, thus reducing MSD values. Cu6O6 is more strongly attracted than the metal clusters by both pristine and functionalized materials, as can be inferred from adsorption energies and MSD values. The forcefield used in this work does not allow reproducing the formation of chemical bonds between atoms. Previous ab initio results on metallic clusters deposited on pristine graphene and graphene oxide surfaces show that there may be oxidation of the particles and the formation of covalent bonds, which strengthens the interaction and favors the adsorption process. To correctly describe the behavior of these systems, more accurate quantum mechanical calculations are needed. In addition, the use of an experimentally associated medium (e.g., water) would enable to address in more detail the equilibrium configurations of the structures, as well as the energetic affinity between the clusters and the graphene-based surfaces. Although these more demanding simulations allows gain a deeper insight into the mechanisms that control the interaction between graphene-based materials and nanoparticles, we think our MD results can provide useful qualitative information and allow handling a larger number of atoms in shorter computational times (typical CPU running times of 10–13 h using an Intel® Xeon® Quad-core (4 Core) were needed for systems containing 3389 atoms). In future work, DFT calculations will be done to validate these results and be able to consider our molecular dynamics simulations as a quicker reliable method to handle these kinds of systems. Another important aspect that should be considered is the use of a solvent in the calculations. The inclusion of water may have severe effects to several important parameters, such as the equilibrium configurations of the surfaces, the spatial arrangement and the dynamics of the metallic clusters, as well as the energetic affinity between the clusters and the surfaces. Thus, future simulations must include solvent to get a more realistic analysis of real systems.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/nano11061378/s1, Video S1: Cu diffusion on a pristine graphene surface, Video S2: CuO diffusion on a pristine nanoplatelet surface, Video S3: Ag diffusion on a PEGylated nanoplatelet surface, Video S4: Cu cluster movement between two lateral PEG chains in a PEGylated graphene surface.

Author Contributions

Methodology, I.L.-T.; Software, I.L.-T.; Formal Analysis, I.L.-T., A.P.-P.; Investigation, I.L.-T., A.P.-P.; Writing—Original Draft, I.L.-T., A.P.-P.; Writing—Review & Editing, A.P.-P.; Supervision, I.L.-T.; Project Administration, I.L.-T. Both authors have read and agreed to the published version of the manuscript.

Funding

The APC was funded by Universidad Europea de Madrid.

Informed Consent Statement

Not applicable.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Acknowledgments

The authors would like to thank Universidad Europea de Madrid for covering publication costs.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hu, W.; Peng, C.; Luo, W.; Lv, M.; Li, X.; Li, D.; Huang, Q.; Fan, C. Graphene-based antibacterial paper. ACS Nano 2010, 4, 4317–4323. [Google Scholar] [CrossRef] [PubMed]
  2. Kumar, P.; Peipei, H.; Rongzhao, Z.; Bo, L. Antibacterial Properties of Graphene-Based Nanomaterials. Nanomaterials 2019, 9, 737. [Google Scholar] [CrossRef] [Green Version]
  3. Chen, J.; Feng, D.; Yongyou, H.; Jian, S.; Yongyang, Y. Antibacterial activity of graphene-modified anode on Shewanella oneidensis MR-1 biofilm in microbial fuel cell. J. Power Sources 2015, 290, 80–86. [Google Scholar] [CrossRef]
  4. Hegab, H.M.; ElMekawy, A.; Zou, L.; Mulcahy, D.; Saint, C.P.; Ginic-Markovic, M. The controversial antibacterial activity of graphene-based materials. Carbon 2016, 105, 362–376. [Google Scholar] [CrossRef]
  5. Li, X.; Li, F.; Gao, Z.; Fang, L. Toxicology of graphene oxide nanosheets against paecilomyces catenlannulatus. Bull. Environ. Contam. Toxicol. 2015, 95, 25–30. [Google Scholar] [CrossRef] [PubMed]
  6. Pieper, H.; Chercheja, S.; Eigler, S.; Halbig, C.E.; Filipovic, M.R.; Mokhir, A. Toxizität von Graphenoxid: Endoperoxide als Ursache. Angew. Chem. 2016, 128, 413–416. [Google Scholar] [CrossRef]
  7. Jayanthi, S.; Eswar, N.K.; Singh, S.A.; Chatterjee, K.; Madras, G.; Sood, A. Macroporous three-dimensional graphene oxide foams for dye adsorption and antibacterial applications. RSC Adv. 2016, 6, 1231–1242. [Google Scholar] [CrossRef] [Green Version]
  8. Chen, J.; Peng, H.; Wang, X.; Shao, F.; Yuan, Z.; Han, H. Graphene oxide exhibits broad-spectrum antimicrobial activity against bacterial phytopathogens and fungal conidia by intertwining and membrane perturbation. Nanoscale 2014, 6, 1879–1889. [Google Scholar] [CrossRef] [PubMed]
  9. Sengupta, I.; Bhattacharya, P.; Talukdar, M.; Neogi, S.; Pal, S.K.; Chakraborty, S. Bactericidal effect of graphene oxide and reduced graphene oxide: Influence of shape of bacteria. Colloid Interface Sci. Commun. 2019, 28, 60–68. [Google Scholar] [CrossRef]
  10. Trinh, L.T.; Quynh, L.A.B.; Hieu, N.H. Synthesis of zinc oxide/graphene oxide nanocomposite material for antibacterial application. Int. J. Nanotechnol. 2018, 15, 108–117. [Google Scholar] [CrossRef]
  11. Raja, A.; Selvakumar, K.; Rajasekaran, P.; Arunpandian, M.; Ashokkumar, S.; Kaviyarasu, K.; Bahadur, S.A.; Swaminathan, M. Visible active reduced graphene oxide loaded titania for photodecomposition of ciprofloxacin and its antibacterial activity. Colloids Surf. A Physicochem. Eng. Asp. 2019, 564, 23–30. [Google Scholar] [CrossRef]
  12. Yang, S.; Lei, P.; Shan, Y.; Zhang, D. Preparation and characterization of antibacterial electrospun chitosan/poly (vinyl alcohol)/graphene oxide composite nanofibrous membrane. Appl. Surf. Sci. 2018, 435, 832–840. [Google Scholar] [CrossRef]
  13. Tu, Q.; Zhang, Q.; Wang, Y.; Jiao, Y.; Xiao, J.; Peng, T.; Wang, J. Antibacterial properties of poly (dimethylsiloxane) surfaces modified with graphene oxide-catechol composite. Prog. Org. Coat. 2019, 129, 247–253. [Google Scholar] [CrossRef]
  14. Hosseini, M.; Sarafbidabad, M.; Fakhri, A.; NoorMohammadi, Z.; Tahami, S. Preparation and characterization of MnS2/chitosan-sodium alginate and calcium alginate nanocomposites for degradation of analgesic drug: Photocorrosion, mechanical, antimicrobial and antioxidant properties studies. Int. J. Biol. Macromol. 2018, 118, 1494–1500. [Google Scholar] [CrossRef]
  15. Liu, Y.; Wen, J.; Gao, Y.; Li, T.; Wang, H.; Yan, H.; Niu, B.; Guo, R. Antibacterial graphene oxide coatings on polymer substrate. Appl. Surf. Sci. 2018, 436, 624–630. [Google Scholar] [CrossRef]
  16. Teymourinia, H.; Salavati-Niasari, M.; Amiri, O.; Yazdian, F. Application of green synthesized TiO2/Sb2S3/GQDs nanocomposite as high efficient antibacterial agent against E. coli and Staphylococcus aureus. Mater. Sci. Eng. C 2019, 99, 296–303. [Google Scholar] [CrossRef]
  17. Wang, X.; Han, Q.; Yu, N.; Wang, T.; Wang, C.; Yang, R. GO-AgCl/Ag nanocomposites with enhanced visible light-driven catalytic properties for antibacterial and biofilm-disrupting applications. Colloids Surf. B Biointerfaces 2018, 162, 296–305. [Google Scholar] [CrossRef] [PubMed]
  18. Naeem, H.; Ajmal, M.; Qureshi, R.B.; Muntha, S.T.; Farooq, M.; Siddiq, M. Facile synthesis of graphene oxide–silver nanocomposite for decontamination of water from multiple pollutants by adsorption, catalysis and antibacterial activity. J. Environ. Manag. 2019, 230, 199–211. [Google Scholar] [CrossRef]
  19. Moghayedi, M.; Goharshadi, E.K.; Ghazvini, K.; Ahmadzadeh, H.; Ranjbaran, L.; Masoudi, R.; Ludwig, R. Kinetics and mechanism of antibacterial activity and cytotoxicity of Ag-RGO nanocomposite. Colloids Surf. B Biointerfaces 2017, 159, 366–374. [Google Scholar] [CrossRef] [PubMed]
  20. Huo, P.; Liu, C.; Wu, D.; Guan, J.; Li, J.; Wang, H.; Tang, Q.; Li, X.; Yan, Y.; Yuan, S. Fabricated Ag/Ag2S/reduced graphene oxide composite photocatalysts for enhancing visible light photocatalytic and antibacterial activity. J. Ind. Eng. Chem. 2018, 57, 125–133. [Google Scholar] [CrossRef]
  21. Deng, C.H.; Gong, J.L.; Zeng, G.M.; Zhang, P.; Song, B.; Zhang, X.G.; Liu, H.Y.; Huan, S.Y. Graphene sponge decorated with copper nanoparticles as a novel bactericidal filter for inactivation of Escherichia coli. Chemosphere 2017, 184, 347–357. [Google Scholar] [CrossRef]
  22. Ouyang, Y.; Cai, X.; Shi, Q.; Liu, L.; Wan, D.; Tan, S. Poly-l-lysine-modified reduced graphene oxide stabilizes the copper nanoparticles with higher water-solubility and long-term additively antibacterial activity. Colloids Surf. B Biointerfaces 2013, 107, 107–114. [Google Scholar] [CrossRef] [PubMed]
  23. Yang, Z.; Hao, X.; Chen, S.; Ma, Z.; Wang, W.; Wang, C.; Yue, L.; Sun, H.; Shao, Q.; Murugadoss, V. Long-term antibacterial stable reduced graphene oxide nanocomposites loaded with cuprous oxide nanoparticles. J. Colloid Interface Sci. 2019, 533, 13–23. [Google Scholar] [CrossRef] [PubMed]
  24. Feng, Y.; Chang, Y.; Sun, X.; Liu, N.; Cheng, Y.; Feng, Y.; Zhang, H.; Li, X. Understanding the Property–Activity Relationships of Polyhedral Cuprous Oxide Nanocrystals in Terms of Reactive Crystallographic Facets. Toxicol. Sci. 2017, 156, 480–491. [Google Scholar] [CrossRef] [PubMed]
  25. Zavareh, S.; Norouzi, E. Impregnation of GO with Cu2+ for enhancement of aniline adsorption and antibacterial activity. J. Water Process Eng. 2017, 20, 160–167. [Google Scholar] [CrossRef]
  26. Liu, H.; Lv, M.; Deng, B.; Li, J.; Yu, M.; Huang, Q.; Fan, C. Laundering durable antibacterial cotton fabrics grafted with pomegranate-shaped polymer wrapped in silver nanoparticle aggregations. Sci. Rep. 2014, 4, 5920. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Liu, Z.; Joshua, T.R.; Hongjie, D. PEGylated Nanographene Oxide for Delivery of Water-Insoluble Cancer Drugs. J. Am. Chem. Soc. 2008, 130, 10876–10877. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Ghosh, S.; Chatterjee, K. Poly(Ethylene Glycol) Functionalized Graphene Oxide in Tissue Engineering: A Review on Recent Advances. Int. J. Nanomed. 2020, 15, 5991–6006. [Google Scholar] [CrossRef] [PubMed]
  29. Chen, D.; Feng, H.; Li, J. Graphene oxide: Preparation, functionalization, and electrochemical applications. Chem. Rev. 2012, 112, 6027–6053. [Google Scholar] [CrossRef]
  30. Yang, K.; Wan, J.; Zhang, S.; Zhang, Y.; Lee, S.T.; Liu, Z. In vivo pharmacokinetics, long-term biodistribution, and toxicology of PEGylated graphene in mice. ACS Nano. 2011, 5, 516–522. [Google Scholar] [CrossRef]
  31. Demirel, E.; Karaca, E.; Durmaz, Y. Effective PEGylation method to improve biocompatibility of graphene derivatives. Eur. Polym. J. 2020, 124, 10876–109504. [Google Scholar] [CrossRef]
  32. Luo, N.; Weber, J.K.; Wang, S. PEGylated graphene oxide elicits strong immunological responses despite surface passivation. Nat. Commun. 2017, 8, 14537. [Google Scholar] [CrossRef] [Green Version]
  33. Kamberaj, H. Molecular Dynamics Simulations in Statistical Physics: Theory and Applications; Springer Nature Switzerland AG: Cham, Switzerland, 2020. [Google Scholar]
  34. Haghighi, S.; Ansari, R.; Ajori, S. A molecular dynamics study on the interfacial properties of carbene-functionalized graphene/polymer nanocomposites. Int. J. Mech. Mater. Des. 2020, 16, 387–400. [Google Scholar] [CrossRef]
  35. Lado-Touriño, I.; Cerpa-Naranjo, A.; Arce, M.P. A molecular dynamics modelling adsorption study of Cu and Ag nanoparticles on pristine and functionalized graphene surfaces. Mater. Today Proc. 2020, 33, 1830–1834. [Google Scholar] [CrossRef]
  36. Xiunan, W.; Yi, L.; Jingcheng, X.; Shengjuan, L.; Fada, Z.; Qian, Y.; Xiao, Z.; Xinluo, Z. Molecular Dynamics Study of Stability and Diffusion of Graphene-Based Drug Delivery Systems. J. Nanomat. 2015, 16, 872079. [Google Scholar]
  37. Safina, L.L.; Baimov, J.A. Molecular dynamics simulation of fabrication of Ni-graphene composite: Temperature effect. Micro Nano Lett. 2020, 15, 176–180. [Google Scholar] [CrossRef]
  38. Bačová, P.; Rissanou, A.; Harmandaris, V.A. Edge-Functionalized Graphene as a Nanofiller: Molecular Dynamics Simulation Study. Macromolecules 2015, 48, 9024–9038. [Google Scholar] [CrossRef]
  39. Georgios, K.; Rakesh, P.; Soumyadipta, S.; Kostas, K.; Arun, V.; Alexey, V.L. Nanostructure and Dynamics of Humidified Nafion-GrapheneOxide Composites via Molecular Dynamics Simulations. J. Phys. Chem. C 2018, 122, 22864–22875. [Google Scholar]
  40. Georgios, K.; Kostas, K. Temperature dependence of dynamic and mechanical properties in poly(acrylic acid)/graphene oxide nanocomposites. Mater. Today Commun. 2017, 13, 359–366. [Google Scholar]
  41. Chen, J.; Chen, L.; Wang, Y. Molecular dynamics simulations of the adsorption of DNA segments onto graphene oxide. J. Phys. D Appl. Phys. 2014, 47, 505401. [Google Scholar] [CrossRef]
  42. Azimi, M.; Mirjavadi, S.S.; Hamouda, A.M.S. Heterogeneities in polymer structural and dynamic properties in graphene and graphene oxide nanocomposites: Molecular dynamics simulations. Macromol. Theory Simul. 2017, 26, 1600086. [Google Scholar] [CrossRef]
  43. Roza-Eleftheria, R.; Karatasos, K. Graphene/poly(ethylene glycol) nanocomposites as studied by molecular dynamics simulations. Mater. Des. 2016, 97, 163–174. [Google Scholar]
  44. Mirhosseini, M.M.; Seyed, S.Z.; Reza, K. Molecular dynamics simulation of functionalized graphene surface for high efficient loading of doxorubicin. J. Mol. Struct. 2017, 1141, 441–450. [Google Scholar] [CrossRef]
  45. Xu, L.; Xiaoning, Y. Molecular dynamics simulation of adsorption of pyrene–polyethylene glycol onto graphene. J. Colloid Interface Sci. 2014, 418, 66–73. [Google Scholar] [CrossRef]
  46. Mahdavi, M.; Fattahi, A.; Tajkhorshid, E.; Nouranian, S. Molecular insights into the loading and dynamics of doxorubicin on PEGylated graphene oxide nanocarriers. ACS Appl. Bio Mater. 2020, 3, 1354–1363. [Google Scholar] [CrossRef]
  47. Kumar, S.; Parekh, S.H. Linking graphene-based material physicochemical properties with molecular adsorption, structure and cell fate. Commun. Chem. 2020, 3, 8. [Google Scholar] [CrossRef] [Green Version]
  48. Bedrov, D.; Smith, G.D. Molecular Dynamics Simulation Study of the Structure of Poly(ethylene oxide) Brushes on Nonpolar Surfaces in Aqueous Solution. Langmuir 2006, 22, 6189–6194. [Google Scholar] [CrossRef] [PubMed]
  49. Biovia Materials Studio, an Integrated, Multi-Scale Modelling Environment. Available online: https://www.3ds.com/products-services/biovia/products/molecular-modeling-simulation/biovia-materials-studio (accessed on 4 March 2021).
  50. Sun, H.; Pengyu, R.; Fried, J.R. The COMPASS force field: Parameterization and validation for phosphazenes. Comput. Theor. Polym. Sci. 1998, 8, 229–246. [Google Scholar] [CrossRef]
  51. Rahmani, F.; Mahdavi, M.; Nouranian, S.; Al-Ostaz, A. Confinement effects on the thermal stability of poly(ethylene oxide)/graphene nanocomposites: A reactive molecular dynamics simulation study. J. Polym. Sci. B Polym. Phys. 2017, 55, 1026–1035. [Google Scholar] [CrossRef]
  52. Cha, J.; Kyoung, W.; Song, K.; Park, S.; Lim, T.; Lee, J.; Kang, H. Quantitative Evaluation of the Dispersion of Graphene Sheets With and Without Functional Groups Using Molecular Dynamics Simulations. Nanoscale Res. Lett. 2016, 11, 136. [Google Scholar] [CrossRef] [Green Version]
  53. Hoover, W.G. Canonical dynamics: Equilibrium phase-space distributions. Phys. Rev. A 1985, 31, 1695–1697. [Google Scholar] [CrossRef] [Green Version]
  54. Nosé, S. A molecular dynamics method for simulations in the canonical ensemble. Mol. Phys. 1984, 52, 255–268. [Google Scholar] [CrossRef]
  55. Böyükata, M.; Jadson, C.B. Structural and Energetic Analysis of Copper Clusters: MD Study of Cun (n = 2–45). J. Braz. Chem. Soc. 2008, 19, 884–893. [Google Scholar] [CrossRef] [Green Version]
  56. Rao, Y.; Yimin, L.; Xiangyuan, C.; Zongwen, L.; Fuyi, C. Optical and magnetic properties of Cu-doped 13-atom Ag nanoclusters. J. Alloys Compd. 2013, 56, 50–55. [Google Scholar] [CrossRef]
  57. Sriram, S.; Rubalya, V.; Chandiramouli, R. Structural and Electronic Properties of CuO, CuO2 and Cu2O Nanoclusters –a DFT Approach. Mater. Sci. 2015, 21, 173–178. [Google Scholar]
  58. Mirhosseini, M.M.; Haddadi-Asl, V.; Zargarian, S.S. Fabrication and characterization of hydrophilic poly(ε-caprolactone)/pluronic P123 electrospun fibers. Appl. Polym. Sci. 2016, 133, 43345. [Google Scholar] [CrossRef]
  59. Einstein, A. A new determination of molecular dimensions. Ann. Phys. 1906, 19, 289–306. [Google Scholar]
  60. Metzler, R.; Jeon, J.H.; Cherstvy, A.G.; Barkai, E. Anomalous diffusion models and their properties: Non-stationarity, non-ergodicity, and ageing at the centenary of single particle tracking. Phys. Chem. Chem. Phys. 2014, 16, 24128–24164. [Google Scholar] [CrossRef] [Green Version]
  61. Rêgo, C.R.C.; Tereshchuk, P.; Oliveira, L.N.; Da Silv, J.L.F. Graphene-supported small transition metal clusters: A density functional theory investigation within van der Waals corrections. Phys. Rev. B 2017, 95, 235422. [Google Scholar] [CrossRef]
  62. Berahman, M.; Sheikhi, M.; Zarifkar, A.; Nadgaran, H. Structural and electronic properties of zigzag graphene nanoribbon decorated with copper cluster. J. Comput. Electron. 2015, 14, 270–279. [Google Scholar] [CrossRef]
  63. Del Castillo, R.M.; Sansores, L.E. Study of the electronic structure of Ag, Au, Pt and Pd clusters adsorption on graphene and their effect on conductivity. Eur. Phys. J. B 2015, 88, 248. [Google Scholar] [CrossRef]
  64. Martin., A.; Lebègue, S.; Erikssonand, O.; Skorodumova, N.V. Adsorption of Cu, Ag, and Au atoms on graphene including van der Waals interactions. J. Phys. Condens. Matter. 2011, 23, 395001. [Google Scholar]
  65. Chen, Y.J.; Lee, Y.T.; Yeh, P.L.; Wang, B.C. Optimized geometry, electronic structure and Ag adsorption property of nanosheet graphene with different symmetry shapes: A theoretical investigation. Res. Chem. Intermed. 2017, 43, 3613–3620. [Google Scholar] [CrossRef]
  66. Mashhadzadeh, A.H.; Vahedi, A.M.; Ardjmand, A.M.; Ahangari, M.G. Investigation of heavy metal atoms adsorption onto graphene and graphdiyne surface: A density functional theory study. Superlattices Microstruct. 2016, 100, 1094–1102. [Google Scholar] [CrossRef]
  67. Hidalgo-Manrique, P.; Xianzhang Lei, R.X.; Mingyu, Z.; Kinloch, I.A.; Young, R.J. Copper/graphene composites: A review. J. Mater. Sci. 2019, 54, 12236–12289. [Google Scholar] [CrossRef] [Green Version]
  68. Sun, X.; Sun, Z.; Zhang, J.; Liu, X.; Yu, F.; Cheng, X.; Zhao, X. Graphene Nucleation Preference at CuO Defects Rather Than Cu2O on Cu (111): A Combination of DFT Calculation and Experiment. ACS Appl. Mater. Interfaces 2018, 10, 43156–43165. [Google Scholar] [CrossRef]
  69. Ko, Y.; Shim, J.; Lee, C.H.; Lee, K.S.; Cho, H.; Lee, K.T.; Son, D.I. Synthesis and characterization of CuO/graphene (Core/shell) quantum dots for electrochemical applications. Mater. Lett. 2018, 217, 113–116. [Google Scholar] [CrossRef]
  70. Mohammadi-Manesh, E.; Vaezzadeh, M.; Saeidi, M. Cu- and CuO-decorated graphene as a nanosensor for H2S detection at room temperature. Surf. Sci. 2015, 636, 36–41. [Google Scholar] [CrossRef]
  71. Sidorenkov, A.; Kolesnikov, S.; Saletsky, A. Molecular dynamics simulation of graphene on Cu (111) with different Lennard-Jones parameters. Eur. Phys. J. B 2016, 89, 220. [Google Scholar] [CrossRef] [Green Version]
  72. Sunil, K. Spreading and orientation of silver nano-drops over a flat graphene substrate: An atomistic investigation. Carbon 2018, 138, 26–41. [Google Scholar]
  73. Valencia, H.; Gil, A.; Frapper, G. Trends in the Adsorption of 3d Transition Metal Atoms onto Graphene and Nanotube Surfaces: A DFT Study and Molecular Orbital Analysis. J. Phys. Chem. C 2010, 114, 14141–14153. [Google Scholar] [CrossRef]
  74. Manade, M.; Viñes, F.; Illas, F. Transition metal adatoms on graphene: A systematic density functional study. Carbon 2015, 95, 525–534. [Google Scholar] [CrossRef]
  75. Dang, T.M.D.; Le, T.T.T.; Fribourg-Blanc, E.; Dang, M.C. Synthesis and optical properties of copper nanoparticles prepared by a chemical reduction method. Adv. Nat. Sci. Nanosci. Nanotechnol. 2011, 2, 015009. [Google Scholar] [CrossRef]
  76. Tejamaya, M.; Römer, I.; Merrifield, R.C.; Lead, J.R. Stability of Citrate, PVP, and PEG Coated Silver Nanoparticles in Ecotoxicology Media. Environ. Sci. Technol. 2012, 46, 7011–7017. [Google Scholar] [CrossRef]
  77. Hemalatha, T.; Akilandeswari, S. Effect of Poly Ethylene Glycol on CuO Nanoparticles and its Antibacterial Application. Int. Lett. Chem. Phys. Astron. 2016, 63, 111–118. [Google Scholar] [CrossRef] [Green Version]
  78. Nursen, S.; Kahraman, E.; Sari, B.; Özgün, A. Synthesis of Some Polymer-Metal Complexes and Elucidation of their Structures. J. Macromol. Sci. Pure Appl. Chem. 2006, 43, 1227–1235. [Google Scholar]
  79. Gervilla, V.; Zarshenas, M.; Sangiovanni, D.G.; Sarakinos, K. Anomalous versus Normal Room-Temperature Diffusion of Metal Adatoms on Graphene. J. Phys. Chem. Lett. 2020, 11, 8030–8936. [Google Scholar] [CrossRef]
  80. Zarshenas, M.; Gervilla, V.; Sangiovanni, D.G.; Sarakinos, K. Room-Temperature Diffusion of Metal Clusters on Graphene. Available online: https://europepmc.org/article/PPR/PPR254964 (accessed on 4 March 2021).
Figure 1. Models of PG, GO_PEG, and clusters. The starting structures used in the simulations before cluster adsorption on both surfaces are shown in the lower part of the figure.
Figure 1. Models of PG, GO_PEG, and clusters. The starting structures used in the simulations before cluster adsorption on both surfaces are shown in the lower part of the figure.
Nanomaterials 11 01378 g001
Figure 2. Models of PGN and GO_PEG_N. The starting structures used in the simulations before Cu cluster adsorption on both surfaces are shown in the lower part of the figure. Each basal graphene plane in the nanoplatelets is displayed in a different color for a clearer view.
Figure 2. Models of PGN and GO_PEG_N. The starting structures used in the simulations before Cu cluster adsorption on both surfaces are shown in the lower part of the figure. Each basal graphene plane in the nanoplatelets is displayed in a different color for a clearer view.
Nanomaterials 11 01378 g002
Figure 3. RDF of metallic clusters and graphene-based surfaces. For the PEGylated structures (bottom), only the RDF of the clusters and the nearest chains of PEG are shown.
Figure 3. RDF of metallic clusters and graphene-based surfaces. For the PEGylated structures (bottom), only the RDF of the clusters and the nearest chains of PEG are shown.
Nanomaterials 11 01378 g003
Figure 4. Representative position of the nanoclusters for pristine and functionalized graphene monolayers after 10 ns of simulation length. The PEG chains closest to the nanoparticles are shown in blue and green.
Figure 4. Representative position of the nanoclusters for pristine and functionalized graphene monolayers after 10 ns of simulation length. The PEG chains closest to the nanoparticles are shown in blue and green.
Nanomaterials 11 01378 g004
Figure 5. Representative position of the nanoclusters for pristine and functionalized graphene nanoplatelets after 10 ns of simulation length. The PEG chains closest to the nanoparticles are shown in blue and green.
Figure 5. Representative position of the nanoclusters for pristine and functionalized graphene nanoplatelets after 10 ns of simulation length. The PEG chains closest to the nanoparticles are shown in blue and green.
Nanomaterials 11 01378 g005
Figure 6. MSD of the clusters versus simulation length: PG (top) and GO_PEG (bottom).
Figure 6. MSD of the clusters versus simulation length: PG (top) and GO_PEG (bottom).
Nanomaterials 11 01378 g006
Figure 7. MSD of the clusters versus simulation length: PGN (top) and GO_PEG_N (bottom).
Figure 7. MSD of the clusters versus simulation length: PGN (top) and GO_PEG_N (bottom).
Nanomaterials 11 01378 g007
Figure 8. Diffusion trajectories of Ag13, Cu13, and Cu6O6 over PG during the last 250 ps of simulation length. The diffusion trajectory of the Cu6O6 cluster is magnified for a clearer view.
Figure 8. Diffusion trajectories of Ag13, Cu13, and Cu6O6 over PG during the last 250 ps of simulation length. The diffusion trajectory of the Cu6O6 cluster is magnified for a clearer view.
Nanomaterials 11 01378 g008
Figure 9. Diffusion trajectories of Ag13, Cu13, and Cu6O6 over PGN during the last 250 ps of simulation length. The diffusion trajectory of the Cu6O6 cluster is magnified for a clearer view.
Figure 9. Diffusion trajectories of Ag13, Cu13, and Cu6O6 over PGN during the last 250 ps of simulation length. The diffusion trajectory of the Cu6O6 cluster is magnified for a clearer view.
Nanomaterials 11 01378 g009
Figure 10. Total MSD (black) and MSD in x, y, and z directions of the clusters interacting with pristine and functionalized graphene monolayers. Axis units: time in ps (x), MSD in Å2 (y).
Figure 10. Total MSD (black) and MSD in x, y, and z directions of the clusters interacting with pristine and functionalized graphene monolayers. Axis units: time in ps (x), MSD in Å2 (y).
Nanomaterials 11 01378 g010
Figure 11. Total MSD (black) and MSD in x, y, and z directions of the clusters interacting with pristine and functionalized graphene nanoplatelets. Axis units: Time in ps (x), MSD in Å2 (y).
Figure 11. Total MSD (black) and MSD in x, y, and z directions of the clusters interacting with pristine and functionalized graphene nanoplatelets. Axis units: Time in ps (x), MSD in Å2 (y).
Nanomaterials 11 01378 g011
Figure 12. Motion in the xy plane of a Cu13 cluster interacting with two PEG chains of a GO_PEG surface.
Figure 12. Motion in the xy plane of a Cu13 cluster interacting with two PEG chains of a GO_PEG surface.
Nanomaterials 11 01378 g012
Table 1. Average ΔE values for the adsorption of Cu, Ag, and CuO nanoclusters on pristine and functionalized graphene-based materials. The values obtained from each individual simulation are shown in parentheses.
Table 1. Average ΔE values for the adsorption of Cu, Ag, and CuO nanoclusters on pristine and functionalized graphene-based materials. The values obtained from each individual simulation are shown in parentheses.
ClusterSystemE (kcal/mol)
Cu13PG −13.72 (−9.11, −14.27, −17.78)
GNP −18.55 (−16.65, −17.91, −21.09)
GO_PEG −27.69 (−22.17, −29.32, −31.58)
GO_PEG_ N −89.54 (−79.11, −85.76, −103.75)
Ag13PG −29.11 (−25.37, −28.83, −33.13)
GNP −37.53 (−30.54, −37.44, −44.61)
GO_PEG −64.05 (−57.89, −60.45, −73.81)
GO_PEG_ N −116.96 (−99.54, −109.67, −141.67)
Cu6O6PG −77.83 (−59.54, −82.37, −91.78)
GNP −114.18 (−102.54, −118.02, −121.98)
GO_PEG −136.01 (−127.27, −130.32, −150.44)
GO_PEG_ N −220.64 (−191.09, −227.43, −243.40)
Table 2. Mean equilibrium distances between the cluster and nearest graphene surface atoms. The values obtained from each individual simulation are shown in parentheses.
Table 2. Mean equilibrium distances between the cluster and nearest graphene surface atoms. The values obtained from each individual simulation are shown in parentheses.
ClusterSystemd (Å)
Cu13PG 2.92 (2.85, 2.90, 3.01)
GNP 2.92 (2.84, 2.91, 3.01)
GO_PEG 2.45 (2.37, 2.46, 2.52)
GO_PEG_ N 2.44 (2.20, 2.46, 2.66)
Ag13PG 3.01 (2.97, 3.02, 3.04)
GNP 2.99 (2.95, 2.98, 3.04)
GO_PEG 2.51 (2.38, 2.56. 2.59)
GO_PEG_ N 2.44 (2.34, 2.40, 2.58)
Cu6O6PG 2.98 (2.94, 2.99, 3.01) (O)
3.45 (3.21, 3.33. 3.45) (Cu)
GNP 2.93 (2.78, 2.97, 3.04) (O)
3.01 (2.92, 3.07, 3.04) (Cu)
GO_PEG 2.55 (2.34, 2.52, 2.79) (O)
3.11 (2.94, 3.07, 3.32) (Cu)
GO_PEG_ N 2.63 (2.54, 2.61, 2.74) (O)
2.94 (2.72, 2.99, 3.11) (Cu)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lado-Touriño, I.; Páez-Pavón, A. Interaction between Graphene-Based Materials and Small Ag, Cu, and CuO Clusters: A Molecular Dynamics Study. Nanomaterials 2021, 11, 1378. https://doi.org/10.3390/nano11061378

AMA Style

Lado-Touriño I, Páez-Pavón A. Interaction between Graphene-Based Materials and Small Ag, Cu, and CuO Clusters: A Molecular Dynamics Study. Nanomaterials. 2021; 11(6):1378. https://doi.org/10.3390/nano11061378

Chicago/Turabian Style

Lado-Touriño, Isabel, and Alicia Páez-Pavón. 2021. "Interaction between Graphene-Based Materials and Small Ag, Cu, and CuO Clusters: A Molecular Dynamics Study" Nanomaterials 11, no. 6: 1378. https://doi.org/10.3390/nano11061378

APA Style

Lado-Touriño, I., & Páez-Pavón, A. (2021). Interaction between Graphene-Based Materials and Small Ag, Cu, and CuO Clusters: A Molecular Dynamics Study. Nanomaterials, 11(6), 1378. https://doi.org/10.3390/nano11061378

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop