Rh(111) Nanomesh by ... - ACS Publications

Aug 28, 2015 - Another investigated mechanism is the functionalization by low-energy ion implantation. After exposing the hexagonal boron nitride ...
1 downloads 0 Views 4MB Size
Subscriber access provided by UPSTATE Medical University Health Sciences Library

Article

Ar Implantation at the hBN/Rh(111) Nanomesh by ab initio Molecular Dynamics Marcella Iannuzzi J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.5b06774 • Publication Date (Web): 28 Aug 2015 Downloaded from http://pubs.acs.org on September 5, 2015

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry C is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Ar Implantation at the hBN/Rh(111) Nanomesh by ab initio Molecular Dynamics Marcella Iannuzzi∗ Institut für Chemie, Universität Zürich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland E-mail: [email protected]

Abstract The controlled intercalation of selected species underneath graphene or hexagonal boron nitride (hBN) on a substrate opens new ways for the functionalization and the tuning of properties of these systems. In this work, the case of hBN on rhodium exposed to a low energy Ar-ion beam [Cun et al., Nano Letters 2013] is further considered theoretically. With the help of ab initio molecular dynamics, the structural rearrangements induced by the interaction between the impacting Ar ion and the substrate are investigated. It is shown that the ion can be intercalated and trapped between the metal and hBN by breaking some BN bonds and penetrating the over layer. The resulting defective structure relaxes quickly while Ar moves to a stable site where a protrusion appears on top of the characteristic super honeycomb lattice of the nanomesh. The presented results provide a first atomistic description of the complex processes leading to the formation of the, so called, boron nitride nanotents.

Keywords: hexagonal boron nitride, implantation, thermal stability, nanomesh, nanotent, density functional theory, molecular dynamics

∗ To

whom correspondence should be addressed

1 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Introduction The growth of hexagonal boron nitride (hBN) on various transition metals leads to formation of highly periodic superstructures 1 that have been extensively studied 2–8 in the past few years. Of particular relevance is the variation of the hBN interaction with the underlaying metal, which is responsible for the topographic corrugation of the monolayer, the spatial modulation of the BN bond length as well as of the electronic structure. Not only these features can be exploited to template self-assembly of molecules 9–11 and selective adsorption processes, 5,12–16 but also specific chemical reaction and directed functionalization have been proved to be feasible. 17–19 Recently, the intercalation of single atoms (H 20 , Ar 21–23 , O 24 , Co 25 , Si 26 ) and small molecules (CO 27–29 ,O2 30,31 ,H2 O 32 ) under sp2 hybridized hBN and graphene layers grown on transition metals has been achieved. In most of the cases, intercalation is thought to occur under sufficient gas pressure by diffusion through extended defects (e.g., grain boundaries, island edges). The success of such processes depends on the competing binding of the adsorbate with the metal in comparison to the strength of the interaction between the metal and the 2D film. Thereby, the intercalation changes the structural and electronic properties of the system, often decoupling the sp2 layer from the metallic substrate 20,29,33 . On the other hand, the removal of the intercalated species is in general also possible upon annealing. Another investigated mechanism is the functionalization by low energy ion implantation. After exposing to “low energy” Ar ion beam the hexagonal boron nitride nanomesh 1 (NM) at room temperature, scanning tunnelling spectroscopy (STM) reveals the presence of protrusions on top of the characteristic super honeycomb structure. These have been interpreted as Ar atoms intercalated between the sp2 bonded boron nitride layer and the rhodium substrate. 22 Similar implantation phenomena have also been obtained with Ne 22 and Rb 34 atoms, and are expected to be possible for other atoms as well. These new features decorating the NM have been named nanotents. They are found at distinct sites in the 3.2 nm super honeycomb unit cell, are stable at room temperature, and survive exposure to air. Irradiation by ion beam has been applied also to graphene supported on Ir(111) 35,36 and Ru(0001) 21 leading to the controlled creation of defects and the implantation 2 ACS Paragon Plus Environment

Page 2 of 27

Page 3 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

of ions in the interface region. Also in these cases, the modulation of the interaction between the metal and graphene characteristic of the Moiré-like superstructure determines the observed patterns of defects. The implantation process is believed to occur by Ar creating local damage on an otherwise pristine hBN over layer and being trapped underneath. Indeed, the structural analysis carried out after the exposure shows an intact sample, where the NM is still recognisable. The only indication of damage are some relatively small irregularities that sometime appear in the vicinity of the protrusions and have been associated to vacancy defects formed during the dosing. 23 Moreover, processing the sample by several annealing cycles, the NM structure is not further deteriorated, but rather the dynamics of the nanotents, by hopping and clustering, can be observed over an otherwise regular super-honeycomb lattice. 37 Better understanding of this type of implantation process is of great relevance because the controlled functionalization of the sp2 layers on transition metals might improve the templating character of the NM, for example to induce molecular assembly or specific modulation patterns of electronic and magnetic properties. In the previous works on the nanotents, density functional theory (DFT) calculations have been used to study the structural properties of the intercalated Ar atom. The site selectivity observed in experiment has been reproduced by the calculations. 22 Moreover, it has been confirmed that nanotents formed by more Ar atoms aggregated in a cluster are favoured. The stability of single and pair vacancies in the hBN layer has also been investigated by DFT, 23,38 concluding that the interaction with Rh favours the stabilisation of the dangling bonds. This work is a first attempt to describe the dynamical processes that lead to the implantation of Ar and the successive stabilisation of the generated defects and of the nanotent itself. DFT calculations are applied in combination with molecular dynamics (MD) to shed light on the microstructural evolution of the NM in the few femtosecond following the impact of Ar, and in the few picosencods following the effective implantation, leading to the stabilisation of the nanotent.

3 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 27

Computational Details The simulations reported in the present work are based on electronic structure calculations using Kohn–Sham density functional theory within the Gaussian and plane waves (GPW) formalism as implemented in the Quickstep module in the CP2K program package. 39 Dual-space pseudopotentials 40 41 are used to describe the interaction of valence electrons with atomic cores. Exchange and correlation are calculated with the Perdew-Burke-Ernzerhof (PBE) 42,43 GGA exchange-correlation (XC) functional. Long-range dispersion interactions are been computed using DFT-D3 formalisms 44 (more details in SI). We employ a slab model that consists of four 12×12 Rh(111) atomic layers, covered on one side by the 13×13 hBN over layer. The computational set-up used for this work reproduces the NM properties in good agreement with experiment. 15,38 The optimisation leads to the corrugation of hBN and the formation of the NM pore where the registry corresponds to (N-top,B-fcc). The pore is defined as the flat area where hBN lays between 2.2 and 2.6 Å over the Rh surface. In this model the pore occupies more than 40% of the NM unit cell. Here, the interaction with the substrate is dominated by the charge redistribution at the N centers, which is polarised towards the Rh surface. In order to optimise this interaction, the BN bonds in the pore are stretched up to two percent more than the equilibrium bond length. The stretching strain is compensated by the compression of the BN bonds in the surrounding region, the wire. Here, the interaction with the metal is mainly dominated by the dispersion forces. The transition region, where hBN bends to connect pore and wire, is called the rim. The Tersoff-Hamann approximation 45,46 for STM simulations (THSTM) is used to reproduce the iso-current topography above the nanomesh. The THSTM image calculated from the described model reproduces very well size and shape of the pore region. The slight over-estimation of the corrugation along the density profile is a common feature of the this approximation. Figure 1 summarises some results obtained with the described NM model. Panel a is a top view of the replicated NM unit cell, where the color code represents the height of the BN bonds over Rh(111). The wire presents two non equivalent crossing sites differing for registry. WXA is where (B-hcp,N4 ACS Paragon Plus Environment

Page 5 of 27

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

fcc), and it is slightly higher (0.2 Å) over Rh than WXB where (B-top,N-hcp). The deviation from the equilibrium BN bond length (1.453 Å) versus the height over Rh is plotted in panel b. The distribution among pore, rim and wire is clearly distinguished, they correspond to 40%, 20% and 36% of the NM unit cell, respectively. The THSTM image obtained for this NM model is reported in panel c together with the topographic profile as calculated along the red line running over the unit cell through the pore (Fig. 1d). Ab initio molecular dynamics 47 (AIMD) simulations are carried out with the described NM model, always within the spin polarised density functional approximation. The temperature is controlled by coupling the system with an external thermostat. 48 The integration time step is tuned on the fly, by setting a maximum value of 0.8 fs, and reducing it to smaller values when needed, in order to avoid that atomic displacements become larger than 0.1 Å per integration step. This expedient is particularly important when simulating Ar impacting against the NM with a starting kinetic energy of tens of eV. In order to more efficiently observe the possible relaxation processes occurring after the implantation of Ar, the metadynamics approach (MTD) 49,50 is also used. MTD accelerates the exploration of the accessible phase space, thus allowing the simulation of rare events that otherwise would require unfeasibly long simulation times. It is based on the definition of a few collective variables (CVs), which are conveniently selected to describe the most relevant structural changes. In the present work, the CVs need to be associated with the formation or the healing of damage in the hBN mesh. A natural choice of CV is the coordination number (related to the number of BN bonds) of those B and N atoms involved in the defective structure. The nudged elastic band method 51 (NEB) is applied in order to estimate the potential energy barrier for the diffusion of the implanted Ar atom. The initial and final states are the optimised configurations of Ar trapped at wire crossing sites. The minimum energy pathway describes the diffusion of Ar over a distance of about 20 Å and it is described by 48 NEB beads.

6 ACS Paragon Plus Environment

Page 6 of 27

Page 7 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Implantation process by AIMD The implantation occurs by sputtering ions with sufficiently large kinetic energy, such that the hBN can be penetrated, but small enough to minimise the damage on the substrate. The initial kinetic energy, the impact angle and the impact site on the NM are expected to determine the success of the implantation process. Ab initio MD is applied to reproduce the first events after the impact of an Ar ion onto the NM. The choice of such a computational demanding methodologies makes it difficult to provide extended sampling and large statistics, which would be required to adequately describe the physics of the low-energy ion beam interacting with the sample. However, the goal of this study is to shine light on the microstructure evolution occurring after the first impact of the ion, addressing in particular the ability of bond breaking and reforming, the relaxation of the electronic structure at defective sites, the role of the metal below hBN and the charge redistribution. Therefore, simulations based on electronic structure theory are better suited. First the system has been thermalised at 100 K, by a short simulation in the canonical ensemble (NVT). After the equilibration, one Ar ion is added above the hBN surface. A series of 23 simulations has been started always changing the initial position (rAr ) and velocity (vAr ) of the ion. These variables determine approximately where Ar is going to strike the NM and with which impact angle θ , i.e., the angle between vAr and the normal to the Rh(111) surface. We distinguish among impact sites on the wire (W), the pore (P), and the rim (R). Ar can be directed against one boron atom (B), one nitrogen atom (N), one bond center (BN), or the center of one BN hexagonal ring (hex). The sketch in figure 2 illustrates three examples of approaching directions. Hence, the varying parameters are the Ar kinetic energy KAr , measured when the ion is at about 2 Å from the closest atom on hBN, the impact angle θ , and the impact site. KAr is chosen between 40 and 105 eV. These energies are most probably in the energy spectrum of the beam used for the experiments resulting in the formation of nanotents. 22 These runs are not meant to be exhaustive of all the possible events occurring in experiment during the ion beam exposure, but, changing the initial conditions, an approximate picture of the effects on the NM is provided. The implantation attempt is not successful when Ar is reflected back 7 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

Page 8 of 27

Page 9 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

recovered. This is true in particular when the impact occurs in the pore, where the hBN is flat, between 2.2 and 2.5 Å above Rh(111), and interact strongly with Rh, leaving not much space for the inwards deformation. Actually, of the seven attempts carried out sending Ar to the pore, none has resulted in a successful implantation. When Ar is directed against the wire, larger effects are observed. The wire is slightly curved, with height between 3.6 and 4.4 Å above Rh(111), the BN bonds are predominantly shorter than the equilibrium bond length. The low stiffness of the entire wire/rim region, which promotes the topographic adaptation of the layer, has been also demonstrated by measuring the response to strain. 52 It turns out that, on the wire, Ar penetrates more easily the sp2 layer, but, in most of the cases, it is reflected back by Rh. The effective trapping occurs only when the outgoing direction brings Ar in an intercalated region, where the hBN is still intact. After the first impact, the kinetic energy of Ar is not anymore sufficient to break again the sp2 layer, and it remains trapped, rebouncing between Rh and hBN until it relaxes into an equilibrium intercalated position. Figure 3 summarises one run started with KAr = 100 eV and θ = 40◦ , where Ar has been directed towards the center of one BN-hexagon on the wire (run-4 in SI-1). The analysed MD trajectory is only 200 fs long, and it shows the impact of Ar, followed by its penetration and trapping below the hBN membrane. Approaching the substrate, the ion first pushes away a B atom of the wire, which moves inwards into the interlayer region (coloured black starting from the top view taken at 30 fs). While Ar breakes through the hBN layer, also a N atom, is knocked out and moves outwards. Its trajectory is also indicated in black in the figure. The hBN lattice in the vicinity of the generated hole is strongly distorted. A second N moves to the Rh layer (50-80 fs), while Ar pushes the layer outwards giving rise to a protrusion. The potential energy of the whole system is plotted in panel d of figure 3. It increases by several eV, due to the formation of the defect. Meanwhile the Ar kinetic energy quickly decreases (Fig. 3e). In such a short simulation time, the thermostat cannot efficiently dissipate the excess heat, hence locally the system is still hot due to the interaction with Ar (3f). After 200 fs, the structure around the impact site is still defective. Ar is trapped, the nanotent has formed, the knocked-out B atom is in the interlayer

9 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

region, not too far from the original site, while N is free above the layer. The two dangling B atoms form a BB bond of 2 Å, well visible in figure 3c, whereas, the two dangling N atoms bend towards the Rh surface, thus strengthening the interaction with the metal. Indeed, the two closest Rh atoms move by 0.3 Å upwards. Due to the presence of the intercalated B, one corner of the pore is slightly lifted, as indicated by the color change of the ball and stick model reported in figure 3a. The Mulliken charges 53 computed along the trajectory show that the positive Ar ion is sharing electrons with the substrate once it is captured in the interlayer region. The charges of the knocked-out B and N atoms undergo larger fluctuations soon after the impact. Afterwards, B turns out to be positively charged, indicating the interaction with the metal, while the N above the surface behaves as a neutral radical. The highest rate of success for the Ar implantation has been obtained by shooting the ion against the rim, taking advantage of the curvature of the layer. When Ar penetrates the BN membrane, it moves towards the wire, where there is more space between hBN and Rh. This occurs also by transferring a relative low kinetic energy. Moreover, the created hole remains smaller and the relaxation of the hBN soon after the impact prevents that Ar moves back into vacuum after being deflected by the Rh surface. In one case of successful implantation (run-19 in SI-1), the MD run has been extended for more than one ps, in order to observe the further relaxation of the structure. At the impact, Ar removes one N and two B atoms from the rim. While the two B atoms quickly come back in place, after bouncing onto the Rh surface, the N atom moves together with Ar under the wire and it is pushed against the over-laying atomic layer. Here, it is grabbed by another B atom. The presence of Ar causes the deformation of the layer and the weakening of the bonds. Hence, the freed N can be inserted in the layer forming a NN bond in a defective geometry. The four top panels of figure 4 show the first femtoseconds of the path of Ar penetrating hBN. After 150 fs, the two B atoms are bonded again, while from below the N atom and Ar push to form the protrusion causing the lifting of another nitrogen, coloured in red in the figure. Finally, the knocked-out N reappears at the surface, bonded with one N and two B atoms. From this state, the simulation is carried on for more

10 ACS Paragon Plus Environment

Page 10 of 27

Page 11 of 27

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

than one picosecond, setting the thermostat at 300 K. Panel b, c, and d of figure 4 show the potential energy, the Ar temperature, and the temperature of the full system together with the temperature of the B and N subsystems, respectively. Ar looses most of its kinetic energy and stabilises under the protrusion at the wire crossing site. Locally hBN presents an irregular structure with one NN (1.84 Å) and two BB (1.91 and 1.89 Å) bonds, giving rise to 5 and 7 member rings. The defect is stabilised such that all B and N atoms are three fold coordinated (Fig. 4e). At the very top of the protrusion is one N atom, as the apex of a triangular pyramid, below which the Ar atom sits. The simulated STM image obtained from this relaxed structure (Fig. 4f) presents the features observed in experiment, i.e. the protrusion at the wire crossing position and the deformed shape of the pore.

This short simulation time is not sufficient to allow for a complete healing of the defect, and the small hole on the protrusion just above the rim is still present. However, from experiment it is known that by annealing the implanted NM, in a temperature range between 300 and 500 K, the super structure does not deteriorate, on the contrary the number of defects is reduced, and the nanotents are still present. Though, they might diffuse and change in size (becoming larger), suggesting the aggregation of Ar atoms in small clusters.

Lattice reconstruction Whether the knocked-out N and B atom trapped in the interlayer region may contribute to the formation of the protrusions is still an open question. In order to address this point, MD can be used to simulate the processes involving the dangling and freed species generated by the Ar impact. To this purpose, we investigate here the case of an Ar ion impacting to the rim with kinetic energy of about 50 eV, leading to the implantation of the ion. The Ar ion strikes against a BN bond. The B and N atoms are knocked out from hBN and move towards the Rh surface. While the N re-bounces on Rh and soon moves back to bind again with the mesh, the freed B atom is stabilised in the interlayer region at an hollow site over Rh(111), not too far from the vacancy hole left in hBN and from the implanted Ar (at 3.6 Å). After a short equilibration carried out at 12 ACS Paragon Plus Environment

Page 12 of 27

Page 13 of 27

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

low temperature, the local structure around the Ar+vacancy defect does not change anymore. Of the three dangling N atoms at the B vacancy, two bind to the closest Rh atoms underneath, thus inducing the downwards bending of the intact BN bonds. Ar instead pushes upwards, creating the protrusion. These two opposite strains applied to the membrane brake another BN bond near the vacant site. Therefore, in the final configuration there are four N atoms and two B atoms with a reduced coordination due to the breaking of the BN bonds. The first snapshot in figure 5a shows this initial state, where the B and N atoms with reduced coordination are coloured in black. The B atom in the interlayer region is hidden under a N atom of the mesh. Starting from this state, a MD simulation at 300 K is carried out to verify whether the defective sp2 layer easily heals and how the protrusion evolves. Since further processes are expected to occur on a longer time scale, the MD simulation is accelerated by means of metadynamics. Two collective variables are defined: the average coordination number of the four dangling N atoms with respect to all B atoms (CN ), and the average coordination number of two displaced B atoms (one trapped on Rh, the other with one BN-bond less) with respect to all N atoms (CB ). The two coordination numbers have starting values around 2 and 0.9, respectively, and should increase to 3 for the pristine hBN layer. Values smaller than the initial ones are to be expected if the defective area is enlarged and more bonds are broken. Five snapshots of the MTD run taken at 1, 2, 6, 10 and 15 ps are displayed in figure 5a. The plot in panel b, instead, shows the time evolution of the two collective variables, while the MTD potential is activating structural changes. The first rearrangement occurs already within the first 2 ps, leading to a concerted stepwise increase of CN and CB . One BN bond is mended, the Ar atom moves farther from the hole towards the wire crossing site, and the interlayer B atom comes closer to the vacancy site. Three N atoms still form only 2 BN bonds, and only the interlayer B atom is not bonded to any N. During the following 10 ps, the Ar wanders around under the wire, while the B atom does not change substantially its site. After 14 ps, the B atom leaves the interlayer region and binds to the three dangling N atoms, finalising the reconstruction of the nanomesh. The representation of the free energy surface in the reduced space of the two CVs is reported in figure 5c. The contour plot shows that the first minimum is more localised, while the second state, where

14 ACS Paragon Plus Environment

Page 14 of 27

Page 15 of 27

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 27

the system remains for longer time, is indeed broader and deeper. The healing process occurs rather quickly crossing free energy barriers of about 1 eV. Soon after the recombination of B with the mesh the simulation has been interrupted. Hence the most stable state, corresponding to the intact nanomesh, has not been sampled.

Implanted species The results of the MD simulations confirm that after implantation and annealing, protrusions are formed at wire crossing positions. Some properties of the implanted structures have been already discussed in previous works. 22,37 Here, some results are summarised, for completeness, and some new features are presented. In particular, details about the structural and electronic features around the implanted species are given, which are considered important to characterise the stabilisation of the defective structure. It should be noticed that the PBE+D3 functional used for this work, which gives better agreement of the NM structure with the experiments, differs from the functional applied for previous simulations reported in Ref. 22,37 This explains some differences in the reported data. Geometry optimisation has been carried out by locating the Ar atom between hBN and Rh(111) at three intercalated sites: WXA, WXB, and center of the pore. In all cases, the protrusion is formed and stabilised, without the formation of any additional defect. The implantation energy is calculated as Eimp = EhBN/Ar/Rh − ENM − EAr ,

(1)

where ENM and EAr are the energies of the pristine NM and of the Ar atom in vacuum, respectively. It turns out that the WX positions are clearly favoured, with a slight preference for WXA (1.38 eV) with respect to WXB (1.45 eV). At the pore center, instead, Eimp increases to 6.62 eV. These results are in agreement with the site selectivity observed in experiment, where protrusions are found only at WX sites, with a prevalence of site A with respect to site B, which depends on coverage and post-annealing treatment. 23 The positive Eimp values indicate that Ar is trapped underneath the mesh. The dominant contributions to this energy are due to the deformation of the hBN layer and 16 ACS Paragon Plus Environment

Page 17 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

its weakened interaction to Rh(111). This latter term is particularly influent when Ar is at the pore, since here the binding interaction to the metal is stronger. Indeed, by running a MD simulation of few ps at 300 K, starting with Ar trapped at the pore, the system relaxes quickly while Ar diffuses to the wire. In a similar MD run, starting with Ar in the WXB position, instead, the atom wanders around, without effectively leaving the initial site. Even though the WXA site would be more stable, the diffusion of Ar underneath the wire requires a rearrangement of the NM atomic structure, which does not occur during such a short sampling. Some additional details on these two MD runs are given in the SI. At WXA, Ar sits on top of a Rh atom, exactly under the center of a BN-hexagon. The energy increase due to the distortion of the NM amounts to 0.67 eV. From the distribution of BN bond length with respect to height (panel a of figure 6) it appears that about 6% of the NM unit cell is involved, and the maximum height from Rh reaches 5.3 Å. The rearrangement of the electronic structure is investigated by means of projected densities of states and the calculation of the electronic density difference, obtained as

∆ρ (r) = ρhBN//Ar/Rh (r) − ρdNM (r) − ρAr (r).

(2)

Here, ρhBN//Ar/Rh is the electronic density of the full system, while ρdNM and ρAr are the densities of the system without Ar and the Ar atom alone, respectively, at exactly the same coordinates as in the full system. The ∆ρ iso-surfaces shown in panel b of figure 6 correspond to the accumulation (red) and the depletion (blue) of electrons at the value of 0.04 el/Å3 . The Ar atom atom appears as compressed between the metal and hBN. Some charge is transferred from Ar to the metal and the remaining positive Ar-center polarises the electronic clouds of the closest B and N forming the protrusion. Hence, the modification remains quite localised at the protrusion site, and it acts as a trapping potential that stabilises the intercalated species. In the plots in figure 6c, the PDOS calculated for the Rh, N, and B atoms directly interacting with Ar are compared to PDOS of the same atoms in the pristine NM. The p and s atomic states of Ar are well below the Fermi energy,

17 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

though some density appears at the Fermi energy due to the re-hybridisation with the d-band of Rh. The Rh-d band has reduced density around the Fermi energy and some features appear at the energy of the Ar p-states (-8 eV). This is another indication of the charge redistribution occurring locally. Also the N-pz and B-pz states distributions are slightly modified. In particular, the broadening of the band towards lower energies indicates the polarisation of the electrons around the positive Ar-center. The experiments show that the functionalised NM survives annealing temperatures up to 600 K. However, the distribution of the protrusions evolves. The number of protrusions decreases but their size increases, suggesting the aggregation of more Ar atoms at the same WX site. 37 The optimisation of the NM structure when two intercalated Ar atoms are present in the unit cell shows that the implantation energy per atom is reduced locating them at the same WX site (structures and STM images are shown in the SI). The main reason is the reduced distortion energy due to the additional strain of the NM, with respect to having two atoms forming two independent protrusions. The implantation energy obtained for different configurations of one and two intercalated Ar are reported in figure 6e. The energy barrier for the diffusion of Ar from WXB to WXA has been also estimated by means on a nudged elastic band simulation. The shortest direct path from B to A runs below the wire and it is about 17 Å long. The computed potential energy barrier is about 1 eV, which is consistent with the fact that experimentally hopping processes are observed. However, the rearrangements of the nanotents is a more complicated process, which most probably involves Ar clusters of larger size and is going to be addressed in future studies.

Conclusions In this work a first description of the dynamic processes leading to the formation of the boron nitride nanotents is presented. By means of AIMD simulations it is possible to follow the local structural evolution induced by the interaction of the hot Ar ion striking against the NM. It is shown that the penetration of Ar through the hBN layer occurs predominantly at the rim and it

18 ACS Paragon Plus Environment

Page 18 of 27

Page 19 of 27

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

induces a relatively limited damage of the NM, also thanks to the efficient self healing dynamics of the supported sp2 layer. A few B and N atoms happen to be knocked-out from their lattice sites, but in general remain in the vicinity of the defect, because of the interaction with Rh. It is demonstrated that the displaced species can be easily re-located in the hBN mesh by self-structural healing. At higher temperatures it turns out that the hBN layer undergoes fast rearrangements, which facilitate the reconstruction after the damage but also the flexible adaptation of the mesh at the presence of Ar. A clear example is the wave-like movement of hBN that allows the displacement of Ar from less favourable sites (e.g. at the pore) to more stable positions. Finally, the stability and properties of the single-Ar nanotent are quickly reviewed, providing additional details on the electronic structure effects. All these results are consistent with the impressive stability and reversibility of the intercalated-NM as demonstrated in experiments. On the other hand, there is still a lot to understand about the dynamics of the Ar clusters and of the defects, which might be responsible for the formation of the nano voids observed upon annealing to temperatures higher than 600K. Moreover, the intercalation of other species, like He and Rb, also leads to peculiar defect structures that need to be better understood. Hence, we are going to extend the study on the funtionalization of the NM in these directions. Acknowledgements: I would like to thank Prof. Thomas Greber and Dr. Huanyao Cun for fruitful discussions and valuable input from the experimental perspective, and Prof. Jürg Hutter for the advise and support . The Swiss National Supercomputer Centre (CSCS) is acknowledged for the generous allocation of computer time.

References 1. Corso, M.; Auwarter, W.; Muntwiler, M.; Tamai, A.; Greber, T.; Osterwalder, J. Boron Nitride Nanomesh. Science 2004, 303, 217–220. 2. Berner, S.; Corso, M.; Widmer, R.; Groening, O.; Laskowski, R.; Blaha, P.; Schwarz, K.; Goriachko, A.; Over, H.; Gsell, S. et al. Boron Nitride Nanomesh: Functionality from a Corrugated Monolayer. Angew Chem Int Edit 2007, 46, 5115–5119. 20 ACS Paragon Plus Environment

Page 20 of 27

Page 21 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

3. Preobrajenski, AB,; Vinogradov, A. S.; Martensson, N. Monolayer of h-BN Chemisorbed on Cu(111) and Ni(111): The Role of the Transition Metal 3d States. Surf Sci 2005, 582, 21–30. 4. Preobrajenski, A. B.; Vinogradov, A. S.; Ng, M. L.; Cavar, E.; Westerstrom, R.; Mikkelsen, A.; Lundgren, E.; Martensson, N. Influence of Chemical Interaction at the Lattice-mismatched hBN/Rh(111) and h-BN/Pt(111) Interfaces on the Overlayer Morphology. Phys Rev B 2007, 75, 245412. 5. Joshi, S.; Ecija, D.; Koitz, R.; Iannuzzi, M.; Seitsonen, A. P.; Hutter, J.; Sachdev, H.; Vijayaraghavan, S.; Bischoff, F.; Seufert, K. et al. Boron Nitride on Cu(111): An Electronically Corrugated Monolayer. Nano Lett 2012, 12, 5821–5828. 6. Brugger, T.; Guenther, S.; Wang, B.; Dil, J. H.; Bocquet, M.-L.; Osterwalder, J.; Wintterlin, J.; Greber, T. Comparison of Electronic Structure and Template Function of Single-layer Graphene and a Hexagonal Boron Nitride Nanomesh on Ru(0001). Phys Rev B 2009, 79. 7. Morscher, M.; Corso, M.; Greber, T.; Osterwalder, J. Formation of Single Layer h-BN on Pd(111). Surf Sci 2006, 600, 3280–3284. 8. Mueller, F.; Huefner, S.; Sachdev, H.; Gsell, S.; Schreck, M. Epitaxial Growth of Hexagonal Boron Nitride Monolayers by a Three-step Boration-oxidation-nitration Process. Phys Rev B 2010, 82, 075405. 9. Schulz, F.; Drost, R.; Hämäläinen, S. K.; Liljeroth, P. Templated Self-assembly and Local Doping of Molecules on Epitaxial Hexagonal Boron Nitride. ACS Nano 2013, 7, 11121–11128. 10. Urgel, J. I.; Schwarz, M.; Garnica, M.; Stassen, D.; Bonifazi, D.; Ecija, D.; Barth, J. V.; Auwarter, W. Controlling Coordination Reactions and Assembly on a Cu(111) Supported Boron Nitride Monolayer. J Am Chem Soc 2015, 137, 2420–2423. 11. Morchutt, C.; Bjoerk, J.; Krotzky, S.; Gutzler, R.; Kern, K. Covalent Coupling via Dehalogenation on Ni(111) Supported Boron Nitride and Graphene. Chem Commun 2015, 51, 2440–2443. 21 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

12. Dil, J. H.; Lobo-Checa, J.; Laskowski, R.; Blaha, P.; Berner, S.; Osterwalder, J.; Greber, T. Surface Trapping of Atoms and Molecules with Dipole Rings. Science 2008, 319, 1824–1826. 13. Koitz, R.; Seitsonen, A. P.; Iannuzzi, M.; Hutter, J. Structural and Electronic Properties of a Large-scale Moiré pattern of Hexagonal Boron Nitride on Cu(111) Studied with Density Functional Theory. Nanoscale 2013, 5, 5589–5595. 14. Koch, H. P.; Laskowski, R.; Blaha, P.; Schwarz, K. Adsorption of Small Gold Clusters on the h-BN/Rh(111) Nanomesh. Phys Rev B 2012, 86, 155404. 15. Iannuzzi, M.; Tran, F.; Widmer, R.; Dienel, T.; Radican, K.; Ding, Y.; Hutter, J.; Gröning, O. Site-selective Adsorption of Phthalocyanine on h-BN/Rh(111) Nanomesh. Phys Chem Chem Phys 2014, 16, 12374–12384. 16. Ma, H.; Ding, Y.; Iannuzzi, M.; Brugger, T.; Berner, S.; Hutter, J.; Osterwalder, J.; Greber, T. Chiral distortion of confined ice oligomers (N = 5,6). Langmuir 2012, 28, 15246–15250. 17. Lyalin, A.; Nakayama, A.; Uosaki, K.; Taketsugu, T. Functionalization of Monolayer h-BN by a Metal Support for the Oxygen Reduction Reaction. J Phys Chem C 2013, 117, 21359–21370. 18. Dienel, T.; Gómez-Díaz, J.; Seitsonen, A. P.; Widmer, R.; Iannuzzi, M.; Radican, K.; Sachdev, H.; Müllen, K.; Hutter, J.; Gröning, O. Dehalogenation and Coupling of a Polycyclic Hydrocarbon on an Atomically Thin Insulator. ACS Nano 2014, 8, 6571–6579. 19. Koitz, R.; Norskov, J. K.; Studt, F. A Systematic Study of Metal-supported Boron Nitride Materials for the Oxygen Reduction Reaction. Phys Chem Chem Phys 2015, 17, 12722–12727. 20. Brugger, T.; Ma, H.; Iannuzzi, M.; Berner, S.; Winkler, A.; Hutter, J.; Osterwalder, J.; Greber, T. Nanotexture Switching of Single-layer Hexagonal Boron Nitride on Rhodium by Intercalation of Hydrogen Atoms. Angew Chem Int Edit 2010, 49, 6120–6124. 21. Cun, H.; Iannuzzi, M.; Hemmi, A.; Osterwalder, J.; Greber, T. Ar Implantation beneath Graphene on Ru(0001): Nanotents and “can-opener” Effect. Surf Sci 2015, 634, 95–102. 22 ACS Paragon Plus Environment

Page 22 of 27

Page 23 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

22. Cun, H.; Iannuzzi, M.; Hemmi, A.; Roth, S.; Osterwalder, J.; Greber, T. Immobilizing Individual Atoms beneath a Corrugated Single Layer of Boron Nitride. Nano Lett 2013, 13, 2098–2103. 23. Cun, H.; Iannuzzi, M.; Hemmi, A.; Osterwalder, J.; Greber, T. Two-Nanometer Voids in Single-Layer Hexagonal Boron Nitride: Formation via the "Can-Opener" Effect and Annihilation by Self-Healing. ACS Nano 2014, 8, 7423–7431. 24. Grånäs, E.; Knudsen, J.; Schroeder, U. A.; Gerber, T.; Busse, C.; Arman, M. A.; Schulte, K.; Andersen, J. N.; Michely, T. Oxygen Intercalation under Graphene on Ir(111): Energetics, Kinetics, and the Role of Graphene Edges. ACS Nano 2012, 6, 9951–9963. 25. Preobrajenski, A. B.; Ng, M. L.; Vinogradov, N. A.; Vinogradov, A. S.; Lundgren, E.; Mikkelsen, A.; Martensson, N. Impact of Oxygen Coadsorption on Intercalation of Cobalt under the h-BN Nanomesh. Nano Lett 2009, 9, 2780–2787. 26. Cui, Y.; Gao, J.; Jin, L.; Zhao, J.; Tan, D.; Fu, Q.; Bao, X. An Exchange Intercalation Mechanism for the Formation of a Two-dimensional Si Structure underneath Graphene. Nano Res 2012, 5, 352–360. 27. Jin, L.; Fu, Q.; Dong, A.; Ning, Y.; Wang, Z.; Bluhm, H.; Bao, X. Surface Chemistry of CO on Ru(0001) under the Confinement of Graphene Cover. J Phys Chem C 2014, 118, 12391– 12398. 28. Grånäs, E.; Andersen, M.; Arman, M. A.; Gerber, T.; Hammer, B.; Schnadt, J.; Andersen, J. N.; Michely, T.; Knudsen, J. CO Intercalation of Graphene on Ir(111) in the Millibar Regime. J Phys Chem C 2013, 117, 16438–16447. 29. Ng, M. L.; Shavorskiy, A.; Rameshan, C.; Mikkelsen, A.; Lundgren, E.; Preobrajenski, A.; Bluhm, H. Reversible modification of the structural and electronic properties of a boron nitride monolayer by CO intercalation. ChemPhysChem 2015, 16, 923–927.

23 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

30. Yang, Y.; Fu, Q.; Wei, M.; Bluhm, H.; Bao, X. Stability of BN/metal Interfaces in Gaseous Atmosphere. Nano Res 2014, 8, 227–237. 31. Blume, R.; Kidambi, P. R.; Bayer, B. C.; Weatherup, R. S.; Wang, Z.-J.; Weinberg, G.; Willinger, M.-G.; Greiner, M.; Hofmann, S.; Knop-Gericke, A. et al. The Influence of Intercalated Oxygen on the Properties of Graphene on Polycrystalline Cu Under Various Environmental Conditions . Phys Chem Chem Phys 2014, 16, 25989–26003. 32. Feng, X.; Maier, S.; Salmeron, M. Water splits epitaxial graphene and intercalates. J Am Chem Soc 2012, 134, 5662–5668. 33. Dong, A.; Fu, Q.; Wei, M.; Liu, Y.; Ning, Y.; Yang, F.; Bluhm, H.; Bao, X. Facile Oxygen Intercalation between Full Layer Graphene and Ru(0001) under Ambient Conditions. Surf Sci 2015, 634, 37–43. 34. de Lima, L. H.; Cun, H. Y.; Hemmi, A.; Kaelin, T.; Greber, T. Note: An Ion Source for Alkali Metal Implantation Beneath Graphene and Hexagonal Boron Nitride Monolayers on Transition Metals. Review of Scientific Instruments 2013, 84. 35. Standop, S.; Lehtinen, O.; Herbig, C.; Lewes-Malandrakis, G.; Craes, F.; Kotakoski, J.; Michely, T.; Krasheninnikov, A. V.; Busse, C. Ion Impacts on Graphene/Ir(111): Interface Channeling, Vacancy Funnels, and a Nanomesh. Nano Lett 2013, 13, 1948–1955. 36. Ahlgren, E. H.; Hämäläinen, S. K.; Lehtinen, O.; Liljeroth, P.; Kotakoski, J. Structural Manipulation of the Graphene/metal Interface with Ar+ Irradiation. Phys Rev B 2013, 88. 37. Cun, H.; Iannuzzi, M.; Hemmi, A.; Osterwalder, J.; Greber, T. Implantation Length and Thermal Stability of Interstitial Ar Atoms in Boron Nitride Nanotents. ACS Nano 2013, 1014–1021. 38. Bacle, P.; Seitsonen, A. P.; Iannuzzi, M.; Hutter, J. Chemical Reactions on Metal-supported Hexagonal Boron Nitride Investigated with Density Functional Theory. Chimia 2014, 68, 596– 601. 24 ACS Paragon Plus Environment

Page 24 of 27

Page 25 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

39. Hutter, J.; Iannuzzi, M.; Schiffmann, F.; VandeVondele, J. CP2K: Atomistic Simulations of Condensed Matter Systems. WIREs Comput Mol Sci 2014, 4, 15–25. 40. Goedecker, S.; Teter, M. P.; Hutter, J. Separable Dual-space Gaussian Pseudopotentials. Phys Rev B 1996, 54, 1703–1710. 41. Krack, M. Pseudopotentials for H to Kr Optimized for Gradient-corrected Exchangecorrelation Functionals. Theor Chem Acc 2005, 114, 145–152. 42. Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys Rev Lett 1996, 77, 3865–3868. 43. Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple (vol 77, pg 3865, 1996). Phys Rev Lett 1997, 78, 1396–1396. 44. Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A Consistent and Accurate ab initio Parametrization of Density Functional Dispersion Correction (DFT-D) for the 94 Elements H-Pu. J Chem Phys 2010, 132, 154104. 45. Tersoff, J.; Hamann, D. Theory and Application for the Scanning Tunneling Microscope. Phys Rev Lett 1983, 50, 1998–2001. 46. Tersoff, J.; Hamann, D. Theory of the Scanning Tunneling Microscope. Phys Rev, B Condens Matter 1985, 31, 805–813. 47. Marx, D.; Hutter, J. Ab Initio Molecular Dynamics; Basic Theory and Advanced Methods; Cambridge University Press, 2009. 48. Bussi, G.; Donadio, D.; Parrinello, M. Canonical Sampling through Velocity Rescaling. J Chem Phys 2007, 126, 014101. 49. Laio, A.; Parrinello, M. Escaping Free-energy Minima. PNAS 2002, 99, 12562–12566.

25 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

50. Iannuzzi, M.; Laio, A.; Parrinello, M. Efficient Exploration of Reactive Potential Energy Surfaces Using Car-Parrinello Molecular Dynamics. Phys Rev Lett 2003, 90, 238302. 51. Henkelman, G.; Uberuaga, B. P. A Climbing Image Nudged Elastic Band Method for Finding Saddle Points and Minimum Energy Paths. J Chem Phys 2000, 113, 9901. 52. Herden, T.; Ternes, M.; Kern, K. Lateral and Vertical Stiffness of the Epitaxial h-BN Monolayer on Rh(111). Nano Lett 2014, 14, 3623–3627. 53. MULLIKEN, R. S. Electronic Population Analysis on Lcao-Mo Molecular Wave Functions .1. J Chem Phys 1955, 23, 1833–1840.

26 ACS Paragon Plus Environment

Page 26 of 27

Page 27 of 27

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment