Structure sensitive enantioselectivity on surfaces: tartaric acid on all surfaces vicinal to Cu(111) †

Comprehensive mapping of enantiospecific surface reactivity versus the crystallographic orientation of Cu( hkl ) surfaces vicinal to Cu(111) has been conducted using a spherically shaped single crystal on which the surface normal vectors, [ hkl ], span all possible orientations lying with 14 1 of the [111] direction. This has allowed direct measurement on 169 diﬀerent Cu( hkl ) surfaces of the two rate constants, k ( hkl ) i and k ( hkl ) e , that determine the kinetics of the vacancy-mediated, explosive decomposition of tartaric acid (TA). The initiation rate constant, k ( hkl ) i , quantifies the kinetics of an initiation step that creates vacancies in the adsorbed TA monolayer. The explosion rate constant, k ( hkl ) e , quantifies the kinetics of a vacancy-mediated explosion step that results in TA decomposition and product desorption. Enantiospecificity is revealed by the dependence of TA decomposition kinetics on the chirality of the local surface orientation. Diastereomerism is demonstrated by the fact that D -TA is more reactive than L -TA on S surfaces while L -TA is more reactive on R surfaces. The time to reach half coverage, t ( hkl ) 1/2 , during isothermal TA decomposition at 433 K allowed determination of the most enantiospecific surface orientation; Cu(754). The ideal Cu(754) surface structure consists of (111) terraces separated by monoatomic steps formed by the (100) and (110) microfacets.


Introduction
The kinetics of surface processes such as catalysis are often dictated by the atomistic structures of the surfaces on which they occur. In principle, this means that atomistic surface structure should be one of the design parameters considered in the development of catalysts. However, rigorous surface design would require the existence and use of quantitative models for prediction of reaction kinetics on surfaces of arbitrary atomic structure. Examples of structure sensitivity in catalytic processes include ammonia synthesis on Fe and Ru, [1][2][3][4][5] hydrocarbon catalysis, [6][7][8][9] and many other reactions. [10][11][12] The basic approach to understanding the origin of the structure sensitivity of these processes requires determination of the reaction mechanism at the level of individual elementary steps and identification of the reaction step (or steps) that are rate controlling. Then one must understand the influence of surface structure on the rate constants for those reaction steps. In this work, a high throughput approach was used to measure the surface structure dependence of the two rate constants, k (hkl) i and k (hkl) e , that determine the overall kinetics of tartaric acid (TA, HO 2 CCH(OH)CH(OH)CO 2 H) decomposition on Cu(hkl) surfaces. The initiation rate constant, k (hkl) i , quantifies the kinetics of an initiation step that creates vacancies in the adsorbed TA monolayer. The explosion rate constant, k (hkl) e , quantifies the kinetics of the vacancy-mediated explosion step that results in auto-catalytic TA decomposition and product desorption. 13 An understanding of the relationship between surface reaction kinetics and surface structure is important to the design and optimization of catalytic processes. This work studies enantiospecificity of the reactions of chiral molecules on chiral surfaces. This is perhaps the subtlest form of surface structure sensitivity or chemical selectivity. Enantioselectivity is solely dictated by the diastereomeric structural relationship between the chiral adsorbate and the chiral surface. 14,15 Chiral surfaces are good candidates for the development and understanding of enantioselective processes because they play an important role in enantioselective chemical processes such as adsorption, catalysis and crystallization. One of the challenges to optimizing enantioselective processes on chiral surfaces is that it requires a fundamental understanding of the enantiospecific interactions between chiral molecules and chiral surfaces. 14,16,17 The optimization of enantioselective processes on chiral surfaces requires an understanding of the influence of surface structure on surface reaction kinetics. This is complicated by the fact that there are an infinite number of surface orientations that can be exposed by any bulk crystal structure. A high throughput approach is needed to study comprehensively the kinetics of enantioselective reactions on chiral surfaces of different crystallographic orientations. Chiral surfaces can be prepared by a variety of means. Many examples of chiral surfaces are made by the adsorption of chiral molecules onto achiral surfaces. The presence of a chiral adsorbate renders the surface chiral by breaking the structural symmetry of the initially achiral surface. 6 Some inorganic materials such as quartz and intermetallic compounds such as GaPd have chiral bulk structures that expose chiral surfaces. 14,18 Intrinsically chiral surfaces can also be obtained from materials with achiral bulk structures, e.g. pure metals. These can be obtained by cleaving the bulk metal lattice along a low-symmetry plane. A face-centered-cubic (FCC) crystal surface is achiral, if any two of its Miller indices (hkl) are equal to one another or, if any of them are zero. This means that, a FCC(hkl) surface is chiral, if h a k a l a h and h Â k Â l a 0. [19][20][21][22] The ideal structures of chiral metal surfaces have low-Miller-index terraces separated by monoatomic kinked step edges as shown for the Cu (3,1,17) surface in the ESI, † Fig. SI 1. 22 These chiral surfaces vary in terms of which low-Miller-index facets form the terraces, the straight portions of the step edges and the kinks. These surfaces also vary in terms of the spacing between the steps separating the terraces and the spacing between the kinks along the step edges. Intrinsically chiral metal surfaces have been shown to exhibit enantiospecific interactions with chiral adsorbates. 13,14,[23][24][25][26] The emergence of the use of single crystal samples in surface science made it possible to measure rates of individual elementary reaction steps and overall catalytic reaction rates on isolated surface orientations. [1][2][3][4][5]27 The ability to perform studies on isolated facets of metal single crystals has revealed surface structure sensitivity. However, that work also reveals the limitations of the flat single crystal approach in being able to provide a comprehensive understanding of structure sensitivity that would allow prediction of surface reaction kinetics across all possible surface orientations. Because of the infinite number of surface orientations, this requires an unrealistic amount of work. For this reason, a high throughput approach is needed to study comprehensively the kinetics of enantioselective reactions on chiral surface orientations spanning crystallographic orientation space. 28,29 Significant development and application of high throughput experimental methods has occurred over the last decade. This allows rapid collection and analysis of datasets spanning broad continuous parameter spaces. [30][31][32][33] Curved single crystal surfaces, which we refer to as surface structure spread single crystals (S 4 Cs), shown in Fig. 1A, expose a continuous distribution of crystallographic surface orientations and, thereby, serve as high throughput libraries for study of structure dependent surface properties. 34 The first use of a curved sample surface was reported by Linder in 1927 in a study of the photoelectric current emitted from the surface of a single crystal cylinder of Zn. 35 The measured photocurrent under UV irradiation varied periodically around the perimeter of the cylinder. In the 1940's Gwathmey used Cu single crystal spheres to reveal surface structure dependent trends in processes such as corrosion and surface oxidation that could be detected visually. [36][37][38][39] In the 1970's Woodruff et al. used Auger spectroscopy to detect structure sensitive oxygen uptake on clean Cu(hkl) surfaces around the perimeter of a single crystal Cu cylinder. 40  a continuous 2D distribution of surface orientations, all treated using the same experimental conditions. This is critical to our study of chiral surfaces. A graphic illustration of a Cu(111)AE141-S 4 C is shown in Fig. 1A. The Cu(111)AE141-S 4 C was used to map the enantiospecific decomposition kinetics of TA decomposition across 169 surface orientations lying within 141 of the [111] direction. The curved dome of the Cu sample is spherical with a 21 mm radius of curvature. The (111) plane is located at the top center point of the dome. As one moves away from the [111] facet, the step density increases with (110) steps being exposed along three symmetry-equivalent directions, and (100) steps along the three opposite directions. The six high symmetry directions divide the S 4 C into six regions containing surfaces alternating between R-and S-chirality. 19,21,48,49 The points along the solid red lines in Fig (111) terraces that are B5 close-packed rows in width separated by monoatomic step edges. The steps are formed by either (100) or (110) microfacets. The use of a S 4 C that exposes both chiral and achiral crystal surfaces, coupled with spatially resolved experimental techniques allows study of enantiospecific surface reaction kinetics across a continuous distribution of surface orientations. This enables the detection of enantiospecificity during surface reactions of chiral adsorbates.
It is important to point out that the real structures of these chiral high-Miller-index surfaces need not be the same as one would predict on the basis of simply cleaving and ideal single crystal bulk lattice. The ideal structures of surfaces are subject to phenomena such as relaxation, reconstruction, step bunching and thermal roughening. The structures of the surfaces of a Cu(111)AE141-S 4 C have been studied using scanning tunneling microscopy. 34 These have shown that the center point is truly a large terrace of Cu(111) orientation. As one moves radially from the (111) center point, monoatomic steps are observed with average spacings that are consistent with expectations based on the lattice spacing of bulk Cu. Furthermore, these steps have roughly the expected orientation as one moves azimuthally around the center point. Although the steps meander with respect to the ideal orientation, there is no sign of gross surface reconstruction or step bunching. High resolution images have also been obtained from the Cu(643) surface. 50 These reveal the expected structure of dense steps with kinks, however, the step edges are thermally roughened such that the kinks are no longer just single atom kinks spaced periodically along the step edge. Although this reduces the areal density of kinks relative to that of the ideal Cu(643) structure, the surface retains its net chirality in the sense that the kinks all have the same orientation of their surrounding low Miller index microfacets. This represents our current understanding of the real structure across the clean Cu(111)AE141-S 4 C surface. That said, it is important to realize also, that adsorbates can complicate this understanding by influencing the underlying structure of the surfaces to which they are adsorbed. 51 TA is an ideal chiral probe for study of enantiospecific surface chemistry because its decomposition mechanism and kinetics are well understood on Cu single crystal surfaces. [52][53][54][55][56][57] After adsorption at 400 K on Cu(110), the saturation coverage has 5/16 TA/Cu atom and with TA singly deprotonated (HO 2 C-CH(OH)CH(OH)CO 2 -Cu). At a lower absolute coverage of 1 4 TA/ Cu atom, the TA is doubly deprotonated copper tartrate (Cu-O 2 CCH(OH)CH(OH)CO 2 -Cu). During heating, the adsorbed TA decomposes into CO 2 , CO, H 2 O, H 2 and leaves some C and O atoms adsorbed to the surface. The decomposition process occurs by a vacancy-mediated explosion reaction in which the reaction rate depends on both the TA coverage, y TA , and the coverage of vacancies or empty sites in the adsorbed monolayer, y vac = (1 À y TA ). 54,55,58 The rate law is given by This surface explosion rate law contains an initiation rate constant, k i , and an explosion rate constant, k e . The objective of this work has been to map the values of k (hkl) i and k (hkl) e across surface orientation space. The exponent in the explosion term describes the order of the reaction with respect to the vacancy concentration, y vac . The second-order suggests that two adjacent vacancies are needed for decomposition of adsorbed TA. The same second-order dependence on vacancy coverage has also been observed for aspartic acid (Asp) decomposition. These non-linear reaction kinetics observed for TA and Asp can result in highly enantiospecific reaction rates on chiral Cu(hkl) surfaces. 13,59 Furthermore, these same explosive reaction kinetics have been observed in studies on Cu(hkl) surfaces with several different orientations: Cu(17,5,1) R&S , Cu(531) R&S , Cu(651) R&S , Cu(643) R&S , and achiral Cu(110) and Cu(111). 13,58,60,61 In order to extract the rate constants k (hkl) i and k (hkl) e we must measure the time-dependent coverage of TA, y (hkl) TA (t;T iso ), at many different surface orientations during isothermal decomposition of adsorbed TA. Herein, we demonstrate the use of a high throughput approach to measure TA coverage, y (hkl) TA (t;T iso ), at T iso = 433 K and at 169 discrete points evenly distributed across the surface of the Cu(111)AE141-S 4 C. From these we have been able to extract values for the surface structure dependence of the rate constants, k (hkl) i and k (hkl) e , and thereby gain insight into the origins of the observed surface structure dependence of TA decomposition and enantiospecificity on chiral Cu(hkl) surfaces that are vicinal to Cu(111).

Apparatus
TA adsorption and decomposition on the Cu(111)AE141-S 4 C was studied in a ThetaProbet. This apparatus consists of a preparation chamber and an analysis chamber. The preparation chamber is equipped with a homemade sublimation source for TA deposition onto the surface. The analysis chamber is equipped with a focused and monochromated X-ray source, a sputter ion gun to clean the surface, and a hemispherical energy analyzer to perform Spatially Resolved X-ray Photoelectron Spectroscopy (SR-XPS). A Tescan Vega3 SEM equipped with a detector for Electron Back Scatter Diffraction (EBSD) was used to determine the rotational orientation of the crystal lattice vectors. A pole figure assembled from the EBSD pattern is shown in the ESI † Fig. SI 3.

S 4 C preparation
The Cu(111)AE141-S 4 C was prepared from a flat Cu(111) single crystal (Monocrystals Inc.) with dimensions of 10 mm in diameter and 2 mm in thickness. One side was shaped by a commercial vendor of nanomachining services using single point diamond turning. The final shape was a section of a 21 mm radius spherical dome. The S 4 C was then mounted by vacuum brazing with a Cu-Ag braze to a 14 Â 14 Â 3 mm 3 Cu block. It was then mounted onto the manipulator of a UHV chamber and subjected to cycles of Ar + ion sputtering and annealing at 900 K until the surface yielded a low energy electron diffraction (LEED) pattern.
To determine the rotational orientation of the crystal lattice vectors, the sample was inserted into a Tescan Vega3 SEM equipped with an EBSD detector. A pole figure from the sample assembled from the EBSD pattern is shown in Fig. SI 3 of the ESI. † The S 4 C on the Cu block was then mounted to a sample holder that allows heating, cooling, and temperature measurement in the ThetaProbet. This apparatus is used for SR-XPS. Within the ThetaProbet, the S 4 C surface was sputtered and annealed further until the surface was clean of contaminants.

TA/S 4 C surface preparation
To prepare the Cu(111)AE141-S 4 C surface for measurements of the TA decomposition kinetics, the surface was held at 400 K while being exposed to TA vapor in the preparation chamber of the ThetaProbet. TA was adsorbed on the surface from the sublimation source held at B390 K. A saturated monolayer of L-or D-TA was adsorbed by exposure of the surface to an L-or D-TA flux for 25 minutes. Once the S 4 C surface was coated with a uniform, saturated monolayer of TA, it was transferred on its holder into the analysis chamber of the ThetaProbet and mounted onto an automated xyz positioning stage. On the positioning stage, the sample can be heated and the temperature measured.

SR-XPS analysis
After first adsorbing the TA monolayer, the spatio-temporal evolution of TA coverage, y (hkl) TA (t), across the Cu(111)AE141-S 4 C was measured using SR-XPS during isothermal heating at 433 K. The O 1s XPS signal was measured on the circular grid of 169 points shown in the ESI, † Fig. SI 4. Measuring the isothermal decomposition kinetics of TA on the Cu(111)AE141-S 4 C surface consisted of cycles in which the sample was heated to 433 K at 1 K s À1 , held at 433 K for a prescribed period of time, cooled back to room temperature, and then analyzed using spatially resolved XPS, once the temperature had decreased to below 400 K. This process was repeated 29 times for a total time of 4000 s at 433 K. TA decomposition does leave some residual O adsorbed on the surface, but at that point there was no further decrease in the O 1s signal with additional annealing at 433 K. The heating and cooling schedule is shown by the hash marks at the top of Fig. 2 indicating the quench times.
Mapping of the TA coverage across the S 4 C surface was accomplished by measuring the intensity (area) of the O 1s XPS signal at 169 points on the surface. The concentric circular grid of points radiated from the center of the S 4 C to a radial distance of 4.5 mm with a spacing of 0.75 mm between circles. The points around each circle were also spaced by B0.75 mm.
Adsorption of TA on the Cu S 4 C surface at 400 K results in its adsorption in a deprotonated state that then decomposes during heating to yield CO 2 , CO, H 2 and other species in the gas phase. This results in the deposition of B10% of the O atoms onto the Cu surface. [54][55][56]58 This atomically adsorbed O remains on the surface indefinitely during further heating at 433 K. The TA coverage can be estimated from the XPS spectra using eqn (2).
where I (hkl) O1s is the intensity of the O 1s XPS signal at a given surface orientation.

Results
The core dataset collected during the course of this study consists of the O 1s spectra obtained from D-and L-TA adsorbed at 169 points on the Cu(111)AE141-S 4 C surface and measured as  The simplest metric of reactivity to be extracted from these decay curves is the half-time, t (hkl) 1/2 , to reach a coverage of y (hkl) TA = 0.5 ML having started at a saturation coverage of 1 ML. It is important to note that the value of t (hkl) 1/2 can be measured with higher accuracy than either the initiation or explosion rate constants. 62 The optimal fit of the D-TA decomposition rate law (eqn (1)) to the measured data for y (111) D-TA (t) is shown by the black curve in Fig. 2. The optimal fit yields estimates of k . 62 Given the condition that k (hkl) e c k (hkl) i , it can be shown that t 1=2 ffi p=2 ffiffiffiffiffiffiffiffi ffi k i k e p . The relationship between t 1/2 and k i k e serves as a metric for the enantiospecificity of the decomposition kinetics. The high accuracy of estimates of t 1/2 from the fits to the y (hkl) D-TA (t) data can be obtained because of the nature of the solution to eqn (1). 62 In order to assess the impact of intermittent temperature quenching on the values of k (hkl) i , k (hkl) e , and t (hkl) 1/2 obtained from y (hkl) TA (t), we conducted measurements of y (hkl) D-TA (t) at three points on the Cu(111)AE141-S 4 C using continuous isothermal conditions, T iso = 433 K, without quenching of the temperature. Fig. 3 (inset) shows a map of the points (blue, black and red squares) studied on the Cu(111)AE141-S 4 C. This experiment first saturated the surface with D-TA and then heated the sample to T iso = 433 K and held it there continuously while obtaining O 1s XP spectra from each of the three points sequentially until the TA coverage dropped to zero at all three. Fig. 3 shows the D-TA coverage decay curves, y (hkl) D-TA (t), during decomposition at 433 K on the three different Cu(hkl) surfaces for non-quenching (Fig. 3B) and quenching experiments (Fig. 3A). Also shown are the values of k (hkl) i , k (hkl) e , and t (hkl) 1/2 estimated from the XPS data obtained using the two different heating profiles. These data demonstrate that quenching introduces a systematic error leading to over estimation of the rate constants by B2 Â.
The full dataset of rate constants at all 169 points on the Cu(111)AE141-S 4 C was obtained from y (hkl) D-TA (t) decay curves obtained from O 1s XPS spectra obtained during 29 periods of intermittent quenching. The coverage maps, y

Discussion
The objective of this work is to map enantiospecific surface reaction rate constants, k (hkl) i and k (hkl) e , versus surface crystallographic orientation. The resulting data can serve as the basis for future development of predictive models for the origins of structure sensitive surface reaction kinetics. The focus of this work is enantioselectivity which is a purely structural form of selectivity. Enantioselectivity is perhaps one of the most subtle forms of reaction selectivity, in the sense that the enantiospecificities of surface reaction energetics are typically only on the order of a few kJ per mole; i.e. DDE RÀS B k B T. 23,24,64,65 This work provides insight into the origins of enantioselectivity in the decomposition of L-and D-TA. A comprehensive study of enantiospecific surface structure sensitivity requires the use of curved samples such as our S 4 C. These span both dimensions of surface orientation. The spherical S 4 C sample spans surface orientations that lack mirror symmetry and, therefore, expose chiral surface structures. The origins of structure sensitivity and enantioselectivity in surface chemistry lie in the surface structure dependence of the rate constants. Estimates of the surface structure dependence of the rate constants k (hkl) i and k (hkl) e were obtained from the coverage decay curves of D-and L-TA on Cu(hkl) surfaces vicinal to Cu(111). Bear in mind that, as indicated by the data in Fig. 3, the rate constants are overestimated systematically by B2Â. This does not, however, influence the relative values of the rate constants obtained from different points on the S 4 C.

Values of the halftime for TA decomposition
The analytical solution for y (hkl) TA (t) allows one to relate features of the decay curves directly to the reaction rate constants. The time to reach a coverage of y TA = 0.5 is given by . 62 The values of t (hkl) 1/2 measured across the Cu(111)AE141-S 4 C span the range t (hkl) 1/2 = 350 to 1300 s. The dependence of t 1/2 on surface orientation, (hkl), is a clear indicator of surface structure sensitivity and enantiospecificity. Rather than using t 1/2 , the values of t 1/2 À1 are used as a metric of enantiospecificity on chiral surfaces because it has the units of a pseudo first-order rate constant and it can be estimated more accurately from the data than either k i or k e . The rotationally averaged values of t 1/2 À1 across the Cu(111)AE141-S 4 C are shown in Fig. 5C for D-and L-TA. The maps of t 1/2 À1 show both structure sensitivity and enaniospecificity. The values of t 1/2 À1 range from 1.0 Â 10 À3 to 3.0 Â 10 À3 s À1 . The regions vicinal to the Cu(111) center of the S 4 C reveal minima in t 1/2 À1 for both enantiomers; i.e. surfaces with no or low step density are relatively unreactive for TA decomposition. Reactivity increases with step density, but enantiospecifically. The values of t 1/2 À1 become significantly enantiospecific for surfaces with polar angles \91 relative to the (111) plane. Enantiospecificity is revealed by the fact that the values of t 1/2 À1 are greater for L-TA in the R-regions than in the S-regions. Consistent with the expected diastereomerism, the values of t 1/2 À1 are greater for D-TA in the S-regions than in the R-regions.

Most enantiospecific surface orientation
One of the most interesting and potentially valuable characteristics of enantiospecific surface chemistry that can be garnered from data obtained using an S 4 C is the orientation of the surface that is most enantiospecific for the reaction of interest. This requires identification of some reaction metric that can be shown to exhibit a maximum in its difference on the R-and S-enantiomorphs of a given surface orientation.
In this work, the most enantiospecific surface orientation vicinal to Cu(111) was determined using the half time, t (hkl) 1/2 , needed to reach a coverage of y (hkl) TA (t 1/2 ) = 0.5. All surface orientations showed that k (hkl) is useful because it serves as a metric for the enantiospecificity of the TA decomposition kinetics. 62 The nature of the solution to eqn (1) allows high accuracy estimates of t (hkl) 1/2 from the fits to the y  Fig. 6A and B, respectively. The values of Dt (hkl) 1/2 in Fig. 6A span the range À400 to 400 s with positive (red) values in the S-chirality zones and negative (blue) values in the R-chirality zones. Note that the white curves in Fig. 6A denote points along which Dt (hkl) 1/2 = 0 and the surfaces exhibit no enantioselectivity. In an ideal world, these white curves would be aligned with the solid and dashed white lines along which the surfaces are achiral.
By determining the locations exhibiting highest magnitude of 1/2 , the most enantiospecific surface orientation was determined. The largest difference is Dt (hkl) 1/2 = 430 s and most enantiospecific surface orientation is Cu(754) R&S . The various equivalent points denoting the positions of the family of Cu{754} R&S surfaces are marked on Fig. 6A with black crosses. At these points, the magnitude of Dt (754) 1/2 reaches B430 seconds which is a reasonably large fraction of a process that takes B2500 s to reach completion. One obvious point to make is that the Cu(754) orientation lies at the edge of the area sampled on the Cu(111)AE141-S 4 C and the data clearly suggest that at polar angles 4121 one will find surfaces with Dt (hkl) 1/2 4 430 s. The ideal chiral surface structures of Cu(754) R&S are depicted in Fig. 6B. These are formed of (111) terraces separated by (310) monoatomic step edges formed by two (100) unit cells separated by kinks formed by one (110) unit cell. The terraces are B4 Cu spacings across or, in other words, B1 nm in width.
In microfacet notation, the ideal surface structure can be described as 4(111) Â (310). 66 The chirality is dictated by the rotational ordering of the (111), (100) and (110) microfacets about the surface normal. It is interesting to note that this structure is similar to that of the Cu(643) R&S surfaces that can be described as having 3(111) Â (310) structure. In other words, (643) is formed by (111) terraces that are B3 atoms wide and separated by (310) step edges. It lies at the same azimuthal angle as (754) but at a longer radial distance from the (111) point; i.e. just off the edge of our Cu(111)AE141-S 4 C. Cu(643) shows very high enantiospecificity for D-and L-TA decomposition. 13 More importantly, that study on a flat Cu(643) R&S single crystal yielded the same enantiospecificity as the current study on the curved Cu(111)AE141-S 4 C; i.e. D-TA exhibits higher reactivity, t 1/2 À1 , on the S-surfaces than on the R-surfaces. The highest enantiospecificity for TA decomposition on our Cu(111)AE141-S 4 C is observed on chiral regions having high step density at large polar angles. We can only speculate on the origin of the observed enantiospecificity. One microscopic explanation could be that the packing density of D-or L-TA on the narrow terraces at the edge of the sample is influenced enantiospecifically by the chirality of the nearby kinks in the (310) step edges. The difference in the packing densities of the two TA enantiomers may influence the rates of formation of the vacancies needed for the explosion reaction.

Conclusions
Both rate constants, k (hkl) i and k (hkl) e , for D-and L-TA decomposition on surfaces vicinal to Cu(111) reveal clear signs of structure sensitivity and enantiospecificity. The TA decomposition kinetics are well-described by the two rate constants across the entire range of surface orientations probed on the Cu(111)AE141-S 4 C. Enantiospecificity is most pronounced in the chiral zones of the S 4 C having high step density; i.e. at polar angles \ 91 from the Cu(111) surface. The values of k (hkl) e show enantiospecificity and structure sensitivity for both TA enantiomers. The clearest manifestation of enantiospecificity and diastereoemerism are observed in the values of t (hkl) 1/2 , with D-TA exhibiting greatest reactivity (lowest t (hkl) 1/2 ) on the surface orientations with S-chirality and the least reactivity on the surface orientation with R-chirality. As dictated by diastereomerism, L-TA exhibits the opposite trend in reactivity. Using Dt (hkl) 1/2 as a metric, we have identified surfaces having (754) orientation as the 'most' enantioselective for TA decomposition on the Cu(111)AE141-S 4 C.

Conflicts of interest
There are no conflicts to declare.