Research Article  Open Access
Properties of the FreeStanding TwoDimensional Copper Monolayer
Abstract
We use density functional theory to study a freestanding 2D copper monolayer. We find that the Cu monolayer is stable in 15 ps ab initio molecular dynamics simulations up to 1200 K. Due to the smaller number of bonds per atom in the 2D layer compared to the 3D bulk, we observe a significantly enhanced energy per bond (0.92 versus 0.58 eV/bond). This is similar to the increase in bond strength going from 3D diamond to 2D graphene. We predict various properties of this material, including band structure and density of states. The freestanding 2D Cu monolayer is hexagonal close packed and is the global minimum structure. One valence electron from each atom is delocalized and is donated into a 2D nearly free electron gas.
1. Introduction
Ultrathin materials including graphene, translation metal dichalcogenides, and boron nitride have generated great interest and exhibit fascinating properties including high carrier mobility [1], quantum Hall effect [2], extraordinary thermal conduction [3], magnetic resonance, and superconductivity [4]. Generally, these ultrathin materials have been based on single layers of lamellar van der Waals materials. These materials are stable in few or even singlelayer structures, due to the strong intralayer chemical bonding and weak van der Waals interlayer interaction [5–8]. In contrast to van der Waals materials, metal atoms have a strong preference for threedimensional (3D) close packed structures. Therefore, monolayer freestanding metallic structures with numerous unsaturated atoms are difficult to stabilize and their synthesis has been challenging. In this paper, we will discuss the stability and properties of freestanding monolayer Cu films.
Bulk copper is estimated to have been discovered in 9000 BC in the Middle East [9]. Copper is a ductile and malleable metal with very high thermal and electrical conductivity. Copper catalysts have numerous applications, including oxidation of alcohols, NO, and CO [10]. Copper oxides include high temperature superconductors with scientific and technological importance [11]. Copper is a key cofactor in many biological oxidationreduction reactions [12, 13]. Copper complexes at the active site in metalloproteins and important contributions to the biological functions [14, 15]. Cyclic Cu_{n} structural motifs (especially Cu_{3}) have been reported to be building blocks of various organic and inorganic systems [16–18]. Copper is essential to all living organisms as a trace dietary mineral because it is a key constituent of the respiratory enzyme complex cytochrome c oxidase. Cu(111) is widely used as substrate to fabricate numerous 2D sheets including graphene [19] and boron nitride [20].
Dimensionality is crucial for material properties. The same chemical compounds can exhibit dramatically different properties depending on whether they are arranged in 0D, 1D, 2D, or 3D crystal structures. This is demonstrated by different carbon materials which can vary from 0D clusters and fullerenes [21] to 1D nanotubes [22], 2D graphene [23], and 3D graphite and diamond. We see that structures of different dimensionality have unique properties.
Copper clusters have also demonstrated sizedependent properties that are completely different from the bulk phase [24]. This demonstrates the importance of quantum confinement and the size or shape of the object on its properties, functionalities, and applications. There have also been interesting experiments on the adsorption of small molecules (such as azobenzene and molecular switches) on the Cu(111) surface [25].
In 2014, Zhao et al. reported the experimental fabrication of small freestanding atomically thick iron patches of up to ten atoms in graphene pores [26]. One can imagine making such patches out of other metals as well. For example, this method could inspire the experimental fabrication of small freestanding Cu monolayers. Motivated by these experimental results on freestanding Fe films, as well as by recent computational predictions of freestanding 2D Au and Ag monolayers [27, 28], we decided to investigate the freestanding Cu monolayer.
In this paper, we present ab initio studies of the freestanding Cu monolayer. We demonstrate that this material is quite stable, as it can be annealed for 15 ps up to 1200 K. Molecular dynamics movies are included in the Supplemental Information in Supplementary Material available online at http://dx.doi.org/10.1155/2016/8429510, and we also evaluate the bond length distribution. We present calculations of the band structure and density of states and compare the 2D material to bulk 3D Cu. We will also compare our results to calculations of freestanding 2D gold and silver monolayers [27, 28]. We hope that these results will inform experimental efforts in this direction. We expect that these materials will exhibit new and unusual properties.
2. Computational Methods
We calculated the electronic structure and total energy using density functional theory (DFT) with the planewave pseudopotential (PWPP) technique as implemented in the Vienna ab initio Simulation Package (VASP) [29]. Ionelectron interactions were implemented with the projectoraugmented wave (PAW) [30, 31] method. The generalized gradient approximation (GGA) and the PBE functional [32] were used with a 500 eV limit for the planewaves. An energy threshold of 10^{−6} eV and a force threshold of 10^{−2} eV/Å were used. The 2D atom sheet was located in the xy plane. To avoid the influence of neighboring layers in the periodic simulation, a 20 Å spacing in the direction was used. For geometry optimization and electronic structure calculations, the Brillouin zone was sampled with a 24 × 24 × 1 Γcentered MonkhorstPack (MP) [33] Kpoint grid.
Phonon calculations were performed using the finite displacement method [34] using CASTEP in Materials Studio 7.0 [35]. The localdensity approximation with ultrasoft pseudopotentials was used. The SCF tolerance was 5 × 10^{−7} eV/atom, and the force tolerance was 10^{−3} eV/Å, and the energy cutoff was 330 eV, at the ultrafine setting. The Brillouin zone was sampled with a 24 × 24 × 1 MP grid for both phonon dispersion and phonon density of states. A 7 × 7 supercell was used with an 8.0 Å cutoff radius. For the Brillouin zone, an average distance of 0.003 Å^{−1} between MonkhorstPack mesh qpoints was used.
Ab initio BornOppenheimer molecular dynamics simulations were performed to assess the thermal stability of the Cu monolayer. The scalarrelativistic DFTD and the TkatchenkoScheffler (TS) method were used in CASTEP [35] in Materials Studio 7.0. The NVT ensemble was used. Typical simulations ran for up to 10 ps with a time step of 1.0 fs (parameters: accuracy fine, SCF = 3 × 10^{−6}, smearing = 0.04, DIIS = 20, NoséHoover method [36], Nosé , Nosé chain length = 2, and fixed center of mass). Materials Studio was also used to create the initial structures and visualize the results. We use periodic boundary conditions for all calculations.
Crystal structure predictions were performed using the evolutionary algorithm as implemented in the USPEX code [37]. In these calculations, initial structures are randomly produced using planar group symmetry. All the newly produced structures are relaxed, and the relaxed energies are used to select structures as parents for the next generation. Each subsequent generation is produced by carefully designed variation operators such as heredity and soft mutation. We considered systems with up to 18 atoms in the unit cell and used 30 structures in each generation, with 60% of the lowestenthalpy structures allowed to produce the next generation through heredity (60%), lattice mutation (30%), and atomic permutation (10%). In addition, two lowestenthalpy structures were allowed to survive into the next generation. The structure relaxations during the evolutionary algorithm were performed using the PBE functional as implemented in VASP. The Visualization for Electronic and Structural Analysis software (VESTA, series 3) [38] was used for visualization and plotting.
3. Results and Discussion
We use the USPEX search program to find the global minimum and low energy isomers of a 2D Cu monolayer [37]. The hexagonal close packed ground state structure of the freestanding Cu monolayer is shown in Figure 1. Other lowlying isomers shown in Figure 1 include trigonal + tetragonal, square, and tetracoordinate configurations. These isomers have relative energies of 0.15, 0.30, and 0.71 eV/atom, respectively. We also observed that the open honeycomb lattice of the Cu monolayer transformed to the square lattice during geometry optimization. This close packed structure maximizes the number of bonds in the plane and is the expected structure. These inplane bonds act to stabilize and maintain the planar configuration of the sheet. The model has P6/mmm (#191) space group symmetry. The unit cell consists of a single Cu atom with lattice constant a = b = 2.44 Å and CuCu bond length of 2.44 Å. We note that this freestanding 2D Cu monolayer holds great potential to be realized experimentally over a small hole.
(a)
(b)
(c)
(d)
The electron charge density is shown in Figure 2(a) and can be used to evaluate the chemical bonding between copper atoms of the monolayer. Further insight into the bonding interaction comes from the analysis of the electron localization function (ELF) shown in Figure 2(b). This provides a good description of electron delocalization in molecules [39] and solids [40] and is a useful tool for chemical bond classification [41]. The ELF values are very low. The distribution of electrons in the 2D Cu monolayer is similar to that for the 3D Cu(111) surface as shown in Figure 2.
(a)
(b)
Next, we evaluated the cohesive energy for the 2D freestanding monolayer and the bulk 3D material. , , and are the total energies of a single Cu atom, one unit cell of the 2D Cu monolayer, and one unit cell of 3D bulk Cu, respectively. Based on our results, the 2D Cu monolayer and the 3D bulk have cohesive energies of 2.76 and 3.48 eV/atom, respectively. For comparison, using the same computational method, the cohesive energies of a 2D and bulk 3D Au layer are 2.82 and 3.11 eV per atom, respectively [28]. The bond strength is increased in comparison to the bulk due to the smaller number of bonds per atom (6 versus 12). Therefore each bond is stronger in the sheet than in the bulk (0.92 eV/bond versus 0.58 eV/bond, resp.).
We tested the stability of this system by calculating the phonon dispersion, which is shown in Figure 3 [34]. All the frequencies are real, indicating kinetic stability. The largest observed frequency is 325 cm^{−1}.
To verify the thermal stability of this new material at both ambient conditions and also at elevated temperatures, we have performed ab initio BornOppenheimer molecular dynamics simulations. Periodic 5 × 6 and 8 × 8 supercells were used in the MD simulations. A series of individual 4–15 ps simulations were carried out at temperatures of 500, 800, 1200, and 1400 K. Snapshots taken at the end of each of these simulations are shown in Figure 4. A survey of bond length extensions shows an increase of up to 16%, 29%, and 33% during last 1 ps of the 500, 800, and 1200 K runs. These results demonstrate that the Cu monolayer has good thermal stability and can maintain its structural integrity during brief 15 ps annealing up to 1200 K. At 1400 K, the system is starting to melt (holes appearing and diffusion) already on these short time scales. These numbers can be compared to the melting point of 1358 K for bulk Cu. At elevated temperatures, we may consider that this is a quasi2D monolayer. The liquid regime for this material will be explored in a separate paper. A movie of one of the molecular dynamics simulations at 800 K for 0.5 ps is available in the Supplemental Information.
In Figure 5, we show the results of the bond length distribution for the copper monolayer at temperatures of 500 K, 800 K, and 1200 K. We observe a significant range of bond lengths, primarily due to motion in the Z direction. Bond length extensions up to 3.3 Å are observed.
To get insight into the electronic properties of this system, we have computed the band structure and the density of states (DOS). As shown in Figure 6(a), the 2D Cu monolayer shows a band structure typical for metals. The metallic character of the Cu monolayer is demonstrated by the Fermi level (E = 0) being located inside the bands and no observation of a bandgap at energy . Therefore, the Cu monolayer is a nonmagnetic metal. The partial density of states (PDOS) is also plotted in Figure 6(a) to visualize the contributions of individual orbitals. There is a sharp peak under the Fermi level mainly contributed by the Cud states. The states at the Fermi level have major contributions from Cup and Cud states, with p > d > s (right). There is apparent hybridization between the Cup and Cud states. The Cud states are close to and below the Fermi level. From the band structure in Figure 6(a), we see that one band crosses the Fermi level. This band is highly dispersive and crosses the Fermi level halfway up (between the Γ point and the Kpoint). This band corresponds to a similar free electron band in bulk 3D Cu and indicates a donation of one electron to the free electron gas. We can see that one valence electron from each atom is delocalized and participates in a nearly free homogeneous electron gas. In Figure 6(b), we show the band structure of bulk 3D Cu for comparison. We see that the band DOS are somewhat different in the two cases, although both are metallic.
(a)
(b)
In Table 1, we compare the predicted properties of the freestanding Cu, Ag, and Au monolayers [27, 28]. All of these freestanding 2D monolayers have been studied using the same parameters and computational methods. At the lowest temperature of 500 K, we observe that the copper and gold systems show very low bond length extensions. Cu has the highest bulk cohesive energy space (by 10%) and the highest bulk melting point (by 1%). Therefore, it is quite surprising that, as a 2D layer, copper ends up somewhat weaker than gold. We see that Au 2D freestanding monolayer has the highest cohesive energy (by 2%). Furthermore, we see that the 2D Au monolayer can be annealed for 10 ps up to 1400 K without destruction of the framework in that time, while the Cu layer only went up to 1200 K. We see that silver has the lowest bulk cohesive energy, which corresponds to the lowest 2D cohesive energy and the lowest observed 2D annealing temperature ranges. We conclude that the freestanding 2D Au layer is unusually robust, but the Cu layer should also be experimentally accessible.

4. Conclusions
In summary, we have studied the freestanding 2D Cu monolayer. We have used molecular dynamics simulations to demonstrate that the bonding network survives 15 ps annealing runs up to 1200 K. This material turns out to be slightly less strong than the freestanding 2D gold layer, but more resilient than the freestanding Ag monolayer. The 2D Cu monolayer is HCP with planar hexacoordinate bonding. With only six atoms per bond, the individual bond energy is substantially larger than the corresponding bulk value. The material is metallic with one electron donated into the free electron gas. We used an evolutionary algorithm to search the 2D space to confirm that this is a global minimum structure. The phonon spectrum was tested to ensure stability against negative modes. A wide range of bond lengths are observed. We hope that researchers will be inspired to study this system or other similar freestanding 2D metal monolayers.
Competing Interests
The authors declare that they have no competing interests.
Acknowledgments
LiMing Yang gratefully acknowledges the support from startup fund (2006013118 and 3004013105) and independent innovation research fund (0118013090) from Huazhong University of Science and Technology and National Natural Science Foundation of China (Grant no. 21673087). LiMing Yang and Eric Ganz thank Matthew Dornfeld for his assistance. The authors thank the Minnesota Supercomputer Institute for support.
Supplementary Materials
Molecular dynamics movie of freestanding Cu monolayer at 800 K for 0.5 ps.
References
 F. Miao, S. Wijeratne, Y. Zhang, U. C. Coskun, W. Bao, and C. N. Lau, “Phasecoherent transport in graphene quantum billiards,” Science, vol. 317, no. 5844, pp. 1530–1533, 2007. View at: Publisher Site  Google Scholar
 Y. Zhang, Y.W. Tan, H. L. Stormer, and P. Kim, “Experimental observation of the quantum Hall effect and Berry's phase in graphene,” Nature, vol. 438, no. 7065, pp. 201–204, 2005. View at: Publisher Site  Google Scholar
 A. A. Balandin, “Thermal properties of graphene and nanostructured carbon materials,” Nature Materials, vol. 10, pp. 569–581, 2011. View at: Publisher Site  Google Scholar
 H. He, P. Bourges, Y. Sidis et al., “Magnetic resonant mode in the singlelayer hightemperature superconductor Tl_{2}Ba_{2}CuO_{6+δ},” Science, vol. 295, no. 5557, pp. 1045–1047, 2002. View at: Publisher Site  Google Scholar
 M. Chhowalla, H. S. Shin, G. Eda, L.J. Li, K. P. Loh, and H. Zhang, “The chemistry of twodimensional layered transition metal dichalcogenide nanosheets,” Nature Chemistry, vol. 5, no. 4, pp. 263–275, 2013. View at: Publisher Site  Google Scholar
 C. Ataca, H. Şahin, and S. Ciraci, “Stable, singlelayer MX_{2} transitionmetal oxides and dichalcogenides in a honeycomblike structure,” Journal of Physical Chemistry C, vol. 116, no. 16, pp. 8983–8999, 2012. View at: Publisher Site  Google Scholar
 S. Lebègue, T. Björkman, M. Klintenberg, R. M. Nieminen, and O. Eriksson, “Twodimensional materials from data filtering and ab initio calculations,” Physical Review X, vol. 3, no. 3, Article ID 031002, 7 pages, 2013. View at: Publisher Site  Google Scholar
 M. Pandey, A. Vojvodic, K. S. Thygesen, and K. W. Jacobsen, “Twodimensional metal dichalcogenides and oxides for hydrogen evolution: a computational screening approach,” The Journal of Physical Chemistry Letters, vol. 6, no. 9, pp. 1577–1585, 2015. View at: Publisher Site  Google Scholar
 CSA—Discovery Guides, A Brief History of Copper, http://csa.com.
 M. Shelef, “Selective catalytic reduction of NO_{x} with Nfree reductants,” Chemical Reviews, vol. 95, no. 1, pp. 209–225, 1995. View at: Publisher Site  Google Scholar
 V. E. Henrich and P. A. Cox, The Surface Science of Metal Oxides, Cambridge University Press, New York, NY, USA, 1994.
 W. Kaim and B. Schwederski, Bioanorgansiche Chemie, Teubner, Stuttgart, Germany, 1991.
 F. A. Cotton, G. Wilkinson, C. A. Murillo, and M. Bochmann, Advanced Inorganic Chemistry, John Wiley & Sons, New York, NY, USA, 6th edition, 1999.
 E. I. Solomon, U. M. Sundaram, and T. E. Machonkin, “Multicopper oxidases and oxygenases,” Chemical Reviews, vol. 96, no. 7, pp. 2563–2605, 1996. View at: Publisher Site  Google Scholar
 R. H. Holm, P. Kennepohl, and E. I. Solomon, “Structural and functional aspects of metal sites in biology,” Chemical Reviews, vol. 96, no. 7, pp. 2239–2314, 1996. View at: Publisher Site  Google Scholar
 H. V. R. Dias, H. V. K. Diyabalanage, M. A. RawashdehOmary, M. A. Franzman, and M. A. Omary, “Bright phosphorescence of a trinuclear copper(I) complex: luminescence thermochromism, solvatochromism, and ‘concentration luminochromism’,” Journal of the American Chemical Society, vol. 125, no. 40, pp. 12072–12073, 2003. View at: Publisher Site  Google Scholar
 X.S. Wang, H. Zhao, Z.R. Qu et al., “An unprecedented 8fold interpenetrating diamondoidlike coordination polymer containing a ${{\text{Cu}}^{+}}_{4}$(RCO_{2})_{4} cluster as connecting node,” Inorganic Chemistry, vol. 42, no. 19, pp. 5786–5788, 2003. View at: Publisher Site  Google Scholar
 H. Eriksson and M. Håkansson, “Mesitylcopper: tetrameric and pentameric,” Organometallics, vol. 16, no. 20, pp. 4243–4244, 1997. View at: Publisher Site  Google Scholar
 X. Li, W. Cai, J. An et al., “Largearea synthesis of highquality and uniform graphene films on copper foils,” Science, vol. 324, no. 5932, pp. 1312–1314, 2009. View at: Publisher Site  Google Scholar
 S. Joshi, D. Ecija, R. Koitz et al., “Boron nitride on Cu(111): an electronically corrugated monolayer,” Nano Letters, vol. 12, no. 11, pp. 5821–5828, 2012. View at: Publisher Site  Google Scholar
 H. W. Kroto, J. R. Heath, S. C. O'Brien, R. F. Curl, and R. E. Smalley, “C_{60}: Buckminsterfullerene,” Nature, vol. 318, no. 6042, pp. 162–163, 1985. View at: Publisher Site  Google Scholar
 S. Iijima, “Helical microtubules of graphitic carbon,” Nature, vol. 354, no. 6348, pp. 56–58, 1991. View at: Publisher Site  Google Scholar
 K. S. Novoselov, A. K. Geim, S. V. Morozov et al., “Electric field in atomically thin carbon films,” Science, vol. 306, no. 5696, pp. 666–669, 2004. View at: Publisher Site  Google Scholar
 K. J. Taylor, C. L. PettietteHall, O. Cheshnovsky, and R. E. Smalley, “Ultraviolet photoelectron spectra of coinage metal clusters,” The Journal of Chemical Physics, vol. 96, no. 4, pp. 3319–3329, 1992. View at: Publisher Site  Google Scholar
 E. R. McNellis, J. Meyer, and K. Reuter, “Azobenzene at coinage metal surfaces: role of dispersive van der Waals interactions,” Physical Review B, vol. 80, no. 20, Article ID 205414, 2054. View at: Publisher Site  Google Scholar
 J. Zhao, Q. Deng, A. Bachmatiuk et al., “Freestanding singleatomthick iron membranes suspended in graphene pores,” Science, vol. 343, no. 6176, pp. 1228–1232, 2014. View at: Publisher Site  Google Scholar
 L.M. Yang, T. Frauenheim, and E. Ganz, “The new dimension of silver,” Physical Chemistry Chemical Physics, vol. 17, no. 30, pp. 19695–19699, 2015. View at: Publisher Site  Google Scholar
 L.M. Yang, M. Dornfeld, T. Frauenheim, and E. Ganz, “Glitter in a 2D monolayer,” Physical Chemistry Chemical Physics, vol. 17, no. 39, pp. 26036–26042, 2015. View at: Publisher Site  Google Scholar
 G. Kresse and J. Hafner, “Ab initio molecular dynamics for liquid metals,” Physical Review B, vol. 47, no. 1, pp. 558–561, 1993. View at: Publisher Site  Google Scholar
 P. E. Blöchl, “Projector augmentedwave method,” Physical Review B, vol. 50, no. 24, pp. 17953–17979, 1994. View at: Publisher Site  Google Scholar
 G. Kresse and D. Joubert, “From ultrasoft pseudopotentials to the projector augmentedwave method,” Physical Review B, vol. 59, no. 3, pp. 1758–1775, 1999. View at: Publisher Site  Google Scholar
 J. P. Perdew, K. Burke, and M. Ernzerhof, “Generalized gradient approximation made simple,” Physical Review Letters, vol. 77, no. 18, pp. 3865–3868, 1996. View at: Publisher Site  Google Scholar
 H. J. Monkhorst and J. D. Pack, “Special points for Brillouinzone integrations,” Physical Review B, vol. 13, no. 12, pp. 5188–5192, 1976. View at: Publisher Site  Google Scholar
 G. Kresse, J. Furthmüller, and J. Hafner, “Ab initio force constant approach to phonon dispersion relations of diamond and graphite,” Europhysics Letters, vol. 32, no. 9, p. 729, 1995. View at: Publisher Site  Google Scholar
 S. J. Clark, M. D. Segall, C. J. Pickard et al., “First principles methods using CASTEP,” Zeitschrift fur Kristallographie, vol. 220, no. 56, pp. 567–570, 2005. View at: Publisher Site  Google Scholar
 G. J. Martyna, M. L. Klein, and M. Tuckerman, “NoséHoover chains: the canonical ensemble via continuous dynamics,” The Journal of Chemical Physics, vol. 97, no. 4, pp. 2635–2643, 1992. View at: Publisher Site  Google Scholar
 A. R. Oganov and C. W. Glass, “Crystal structure prediction using ab initio evolutionary techniques: principles and applications,” The Journal of Chemical Physics, vol. 124, no. 24, Article ID 244704, 2006. View at: Google Scholar
 K. Momma and F. Izumi, “VESTA 3 for threedimensional visualization of crystal, volumetric and morphology data,” Journal of Applied Crystallography, vol. 44, no. 6, pp. 1272–1276, 2011. View at: Publisher Site  Google Scholar
 A. D. Becke and K. E. Edgecombe, “A simple measure of electron localization in atomic and molecular systems,” The Journal of Chemical Physics, vol. 92, no. 9, pp. 5397–5403, 1990. View at: Publisher Site  Google Scholar
 A. Savin, O. Jepsen, J. Flad, O. K. Andersen, H. Preuss, and H. G. von Schnering, “Electron localization in solidstate structures of the elements: the diamond structure,” Angewandte Chemie (International Edition in English), vol. 31, no. 2, pp. 187–188, 1992. View at: Publisher Site  Google Scholar
 B. Silvi and A. Savin, “Classification of chemical bonds based on topological analysis of electron localization functions,” Nature, vol. 371, no. 6499, pp. 683–686, 1994. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2016 LiMing Yang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.