We use density functional theory to study a free-standing 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 free-standing 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.
startup fund20060131183004013105independent innovation research fund0118013090Huazhong University of Science and TechnologyNational Natural Science Foundation of China216730871. 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 single-layer 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 three-dimensional (3D) close packed structures. Therefore, monolayer free-standing 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 free-standing 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 oxidation-reduction reactions [12, 13]. Copper complexes at the active site in metalloproteins and important contributions to the biological functions [14, 15]. Cyclic Cun structural motifs (especially Cu3) 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 size-dependent 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 free-standing 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 free-standing Cu monolayers. Motivated by these experimental results on free-standing Fe films, as well as by recent computational predictions of free-standing 2D Au and Ag monolayers [27, 28], we decided to investigate the free-standing Cu monolayer.
In this paper, we present ab initio studies of the free-standing 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 free-standing 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 plane-wave pseudopotential (PWPP) technique as implemented in the Vienna ab initio Simulation Package (VASP) [29]. Ion-electron interactions were implemented with the projector-augmented wave (PAW) [30, 31] method. The generalized gradient approximation (GGA) and the PBE functional [32] were used with a 500 eV limit for the plane-waves. 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 z direction was used. For geometry optimization and electronic structure calculations, the Brillouin zone was sampled with a 24 × 24 × 1 Γ-centered Monkhorst-Pack (MP) [33] K-point grid.
Phonon calculations were performed using the finite displacement method [34] using CASTEP in Materials Studio 7.0 [35]. The local-density 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 Monkhorst-Pack mesh q-points was used.
Ab initio Born-Oppenheimer molecular dynamics simulations were performed to assess the thermal stability of the Cu monolayer. The scalar-relativistic DFT-D and the Tkatchenko-Scheffler (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é Q=2, 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 lowest-enthalpy structures allowed to produce the next generation through heredity (60%), lattice mutation (30%), and atomic permutation (10%). In addition, two lowest-enthalpy 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 free-standing Cu monolayer is shown in Figure 1. Other low-lying 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 in-plane 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 Cu-Cu bond length of 2.44 Å. We note that this free-standing 2D Cu monolayer holds great potential to be realized experimentally over a small hole.
Different motifs of Cu monolayer: (a) hexagonal close packed, (b) trigonal + tetragonal, (c) square, and (d) tetracoordinate.
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) Charge density and (b) ELF of the 2D Cu monolayer and 3D bulk Cu.
Next, we evaluated the cohesive energy for the 2D free-standing monolayer and the bulk 3D material. (1)Ecoh2D=xECu-atom-ECu-monolayerx,Ecoh3D=xECu-atom-ECu-bulkx.ECu-atom, ECu-monolayer, and ECu-bulk 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.
Phonon dispersion and phonon density of states of the 2D Cu monolayer. Γ (0, 0, 0), M (0, 1/2, 0), and K (1/3, 2/3, 0) refer to special points in the first Brillouin zone in reciprocal space.
To verify the thermal stability of this new material at both ambient conditions and also at elevated temperatures, we have performed ab initio Born-Oppenheimer 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 quasi-2D 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.
Snapshots of the final frame of each molecular dynamics simulation from 500 to 1400 K (top and side views). Bonds to atoms outside this section exist but are not shown.
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.
Bond length distribution for the copper monolayer at 500 K, 800 K, and 1200 K. A black vertical line is drawn at 2.41 Å to show the 0 K bond length. We observe that there is a wide distribution of bond lengths due primarily to motion in the Z direction.
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 E=0. 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 Cu-d states. The states at the Fermi level have major contributions from Cu-p and Cu-d states, with p > d > s (right). There is apparent hybridization between the Cu-p and Cu-d states. The Cu-d 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 K-point). 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 d band DOS are somewhat different in the two cases, although both are metallic.
(a) Band structure and DOS of the free-standing 2D Cu monolayer. (b) Band structure and DOS of the bulk 3D Cu. The Fermi level is at 0 eV.
In Table 1, we compare the predicted properties of the free-standing Cu, Ag, and Au monolayers [27, 28]. All of these free-standing 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 free-standing 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 free-standing 2D Au layer is unusually robust, but the Cu layer should also be experimentally accessible.
Comparison of the predicted properties of Cu, Ag, and Au free-standing 2D monolayers.
Bond length, Å
Bond length extension
Bulk
2D cohesive energy,eV/atom
Bulk cohesive energy,eV/atom
In-plane stiffness, N/m
Highest phonon freq, cm−1
500 K
800 K
1200 K
1400 K
Melting Pt, K
Cu
2.41
16%
29%
33%
Melted
1356
2.76
3.48
325
Ag
2.79
29%
36%
Melted
1235
2.01
2.49
31
230
Au
2.73
12%
17%
21%
43%
1338
2.82
3.11
56
220
4. Conclusions
In summary, we have studied the free-standing 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 free-standing 2D gold layer, but more resilient than the free-standing 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 free-standing 2D metal monolayers.
Competing Interests
The authors declare that they have no competing interests.
Acknowledgments
Li-Ming 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). Li-Ming Yang and Eric Ganz thank Matthew Dornfeld for his assistance. The authors thank the Minnesota Supercomputer Institute for support.
MiaoF.WijeratneS.ZhangY.CoskunU. C.BaoW.LauC. N.Phase-coherent transport in graphene quantum billiards200731758441530153310.1126/science.11443592-s2.0-34548672547ZhangY.TanY.-W.StormerH. L.KimP.Experimental observation of the quantum Hall effect and Berry's phase in graphene2005438706520120410.1038/nature042352-s2.0-27744475163BalandinA. A.Thermal properties of graphene and nanostructured carbon materials20111056958110.1038/nmat3064HeH.BourgesP.SidisY.UlrichC.RegnaultL. P.PailhèsS.BerzigiarovaN. S.KolesnikovN. N.KeimerB.Magnetic resonant mode in the single-layer high-temperature superconductor Tl2Ba2CuO6+δ200229555571045104710.1126/science.10678772-s2.0-0037039824ChhowallaM.ShinH. S.EdaG.LiL.-J.LohK. P.ZhangH.The chemistry of two-dimensional layered transition metal dichalcogenide nanosheets20135426327510.1038/nchem.15892-s2.0-84875413255AtacaC.ŞahinH.CiraciS.Stable, single-layer MX2 transition-metal oxides and dichalcogenides in a honeycomb-like structure2012116168983899910.1021/jp212558p2-s2.0-84860361457LebègueS.BjörkmanT.KlintenbergM.NieminenR. M.ErikssonO.Two-dimensional materials from data filtering and ab initio calculations201333703100210.1103/PhysRevX.3.031002PandeyM.VojvodicA.ThygesenK. S.JacobsenK. W.Two-dimensional metal dichalcogenides and oxides for hydrogen evolution: a computational screening approach2015691577158510.1021/acs.jpclett.5b00353CSA—Discovery Guides, A Brief History of Copper, http://csa.comShelefM.Selective catalytic reduction of NOx with N-free reductants199595120922510.1021/cr00033a008HenrichV. E.CoxP. A.1994New York, NY, USACambridge University PressKaimW.SchwederskiB.1991Stuttgart, GermanyTeubnerCottonF. A.WilkinsonG.MurilloC. A.BochmannM.19996thNew York, NY, USAJohn Wiley & SonsSolomonE. I.SundaramU. M.MachonkinT. E.Multicopper oxidases and oxygenases19969672563260510.1021/cr950046o2-s2.0-7744236144HolmR. H.KennepohlP.SolomonE. I.Structural and functional aspects of metal sites in biology19969672239231410.1021/cr95003902-s2.0-0011509370DiasH. V. R.DiyabalanageH. V. K.Rawashdeh-OmaryM. A.FranzmanM. A.OmaryM. A.Bright phosphorescence of a trinuclear copper(I) complex: luminescence thermochromism, solvatochromism, and ‘concentration luminochromism’200312540120721207310.1021/ja036736o2-s2.0-0141958041WangX.-S.ZhaoH.QuZ.-R.YeQ.ZhangJ.XiongR.-G.YouX.-Z.FunH.-K.An unprecedented 8-fold interpenetrating diamondoid-like coordination polymer containing a Cu+4(RCO2)4 cluster as connecting node200342195786578810.1021/ic0344620ErikssonH.HåkanssonM.Mesitylcopper: tetrameric and pentameric199716204243424410.1021/om970462u2-s2.0-4243785018LiX.CaiW.AnJ.KimS.NahJ.YangD.PinerR.VelamakanniA.JungI.TutucE.BanerjeeS. K.ColomboL.RuoffR. S.Large-area synthesis of high-quality and uniform graphene films on copper foils200932459321312131410.1126/science.11712452-s2.0-66749119012JoshiS.EcijaD.KoitzR.IannuzziM.SeitsonenA. P.HutterJ.SachdevH.VijayaraghavanS.BischoffF.SeufertK.BarthJ. V.AuwärterW.Boron nitride on Cu(111): an electronically corrugated monolayer201212115821582810.1021/nl303170m2-s2.0-84869161583KrotoH. W.HeathJ. R.O'BrienS. C.CurlR. F.SmalleyR. E.C60: Buckminsterfullerene1985318604216216310.1038/318162a02-s2.0-1542743726IijimaS.Helical microtubules of graphitic carbon19913546348565810.1038/354056a02-s2.0-0342819025NovoselovK. S.GeimA. K.MorozovS. V.JiangD.ZhangY.DubonosS. V.GrigorievaI. V.FirsovA. A.Electric field in atomically thin carbon films2004306569666666910.1126/science.11028962-s2.0-7444220645TaylorK. J.Pettiette-HallC. L.CheshnovskyO.SmalleyR. E.Ultraviolet photoelectron spectra of coinage metal clusters19929643319332910.1063/1.4619272-s2.0-36449000057McNellisE. R.MeyerJ.ReuterK.Azobenzene at coinage metal surfaces: role of dispersive van der Waals interactions2054802020541410.1103/PhysRevB.80.205414ZhaoJ.DengQ.BachmatiukA.SandeepG.PopovA.EckertJ.RümmeliM. H.Free-standing single-atom-thick iron membranes suspended in graphene pores201434361761228123210.1126/science.12452732-s2.0-84903362414YangL.-M.FrauenheimT.GanzE.The new dimension of silver20151730196951969910.1039/c5cp03465e2-s2.0-84937803462YangL.-M.DornfeldM.FrauenheimT.GanzE.Glitter in a 2D monolayer20151739260362604210.1039/c5cp04222d2-s2.0-84942774562KresseG.HafnerJ.Ab initio molecular dynamics for liquid metals199347155856110.1103/PhysRevB.47.5582-s2.0-12844286241BlöchlP. E.Projector augmented-wave method19945024179531797910.1103/PhysRevB.50.179532-s2.0-25744460922KresseG.JoubertD.From ultrasoft pseudopotentials to the projector augmented-wave method19995931758177510.1103/physrevb.59.17582-s2.0-0011236321PerdewJ. P.BurkeK.ErnzerhofM.Generalized gradient approximation made simple199677183865386810.1103/PhysRevLett.77.38652-s2.0-4243943295MonkhorstH. J.PackJ. D.Special points for Brillouin-zone integrations197613125188519210.1103/PhysRevB.13.51882-s2.0-1842816907KresseG.FurthmüllerJ.HafnerJ.Ab initio force constant approach to phonon dispersion relations of diamond and graphite199532972910.1209/0295-5075/32/9/005ClarkS. J.SegallM. D.PickardC. J.HasnipP. J.ProbertM. I. J.RefsonK.PayneM. C.First principles methods using CASTEP20052205-656757010.1524/zkri.220.5.567.650752-s2.0-20144380685MartynaG. J.KleinM. L.TuckermanM.Nosé-Hoover chains: the canonical ensemble via continuous dynamics19929742635264310.1063/1.4639402-s2.0-36449000062OganovA. R.GlassC. W.Crystal structure prediction using ab initio evolutionary techniques: principles and applications200612424244704MommaK.IzumiF.VESTA 3 for three-dimensional visualization of crystal, volumetric and morphology data20114461272127610.1107/S00218898110389702-s2.0-82055197100BeckeA. D.EdgecombeK. E.A simple measure of electron localization in atomic and molecular systems19909295397540310.1063/1.4585172-s2.0-36549100412SavinA.JepsenO.FladJ.AndersenO. K.PreussH.von SchneringH. G.Electron localization in solid-state structures of the elements: the diamond structure199231218718810.1002/anie.1992018712-s2.0-33745374559SilviB.SavinA.Classification of chemical bonds based on topological analysis of electron localization functions1994371649968368610.1038/371683a02-s2.0-0027946619