Abstract

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.

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 single-layer structures, due to the strong intralayer chemical bonding and weak van der Waals interlayer interaction [58]. 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 [1618]. 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 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é , 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.

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.

Next, we evaluated the cohesive energy for the 2D free-standing 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 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.

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 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 band DOS are somewhat different in the two cases, although both are metallic.

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.

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.

Supplementary Materials

Molecular dynamics movie of freestanding Cu monolayer at 800 K for 0.5 ps.

  1. Supplementary Material