Exploring the Structure–Activity Relationship on Platinum Nanoparticles

The design of active and stable Pt-based nanoscale electrocatalysts for the oxygen reduction reaction (ORR) plays the central role in ameliorating the efficiency of proton exchange membrane fuel-cells towards future energy applications. On that front, theoretical studies have contributed significantly to this research area by gaining deeper insights and understanding of the ongoing processes. In this work, we present an approach capable of characterizing differently-shaped platinum nanoparticles undergoing thermally- and adsorbate-induced restructuring of the surface. Further, by performing ReaxFF-Grand Canonical Molecular Dynamics simulations we explored the water formation on these roughened (“realistic”) nanoparticles in a H2/O2 environment. Taking into consideration the coverage of oxygen-containing intermediates and occurring surface roughening the nanoparticles’ activities were explored. Hereby, we succeeded in locally resolving the water formation on the nanoparticles’ surfaces, allowing an allocation of the active sites for H2O production. We observed that exposed, low-coordinated sites as well as pit-shaped sites originating from roughening of vertices and edges are most active towards H2O formation.


Introduction
Advances in the performance of proton exchange membrane fuel cells (PEMFCs) for economically viable future applications are aimed for in various areas of research, both theoretically and experimentally. As the cathodic oxygen reduction reaction (ORR) on platinum-based electrodes is still considered as one of the bottleneck reactions due to its sluggish kinetics, considerable efforts have been undertaken to its understanding [1]. It has become accepted knowledge that this kinetic hindrance is largely due to oxygenated reaction intermediates expressing too strong binding to the platinum electrode, their removal thereby posing the rate and potential limiting step [2,3]. Following the Sabatier principle, the ORR activity is strongly linked to the adsorption energies of its reaction intermediates. Merging these ideas, socalled volcano plots have been established as guideposts for the catalyst design, correlating surface adsorption energies with predictions of the catalytic activity of materials. Based on these correlations, it has been proposed that an optimal catalyst should express a weaker binding towards hydroxide molecules (OH*) by ~ 0.1 eV compared to Pt(111) [4][5][6][7][8].
Further screening methods, also based on density functional theory (DFT) methods, have been established over the last years, allowing to connect measurable quantities with predictions and trends in the catalytic activity and behavior. One pioneering approach among those methods is surely the d-band model by Hammer and Nørskov, aiming at predicting the reactivity of late transition metals based on the d-band center of the surface atoms [9].
Also, by applying the concept of generalized coordination numbers (gCN) information about the local structure of catalysts can be obtained and correlated to trends in e.g. adsorption energies [8,[10][11][12]. Following intuition, the binding energy of oxygenated species should decrease Electronic supplementary material The online version of this article (https ://doi.org/10.1007/s1124 4-020-01324 -w) contains supplementary material, which is available to authorized users. with increasing coordination of the active site [13][14][15]. The gCNs constitute a first-order extension of the usual coordination numbers by including the second nearest neighbors into consideration. Thus, "structure sensitivity" should be captured more thoroughly; in particular, the comparison between nanoparticles and extended surfaces can be improved by covering "finite-size effects". Now, coupling Sabatier's volcano plots with gCNs results in coordination-activity plots, assigning functional groups of Pt with higher gCN than the model Pt(111) surface higher activities [8,12,16]. Platinum nanoparticles can enhance the performance of the ORR while their high surface area to volume ratios also reduce catalyst loading. However, the challenge remains to find optimum sizes, shapes and morphologies for an enhanced catalytic process [17][18][19][20].
Apart from the previously mentioned criteria for optimized ORR catalysts, long-term stability under the reaction's harsh conditions is of utter importance. Yet, the above-mentioned simple screening models cannot account for thermally-or adsorbate-induced reconstructions of platinum surfaces or nanoparticles and as a consequence thereof lack predictions about possible changes in the catalytic behavior. Alternatively to DFT methods, reactive force field-based approaches such as ReaxFF are predestined to model such structural changes and reactions of extended systems containing several thousand atoms [21][22][23][24]. ReaxFF features the description of electronicstructure-dependent chemical properties of atoms, though without an explicit description of individual electrons, and is highly transferable up to complex reaction chemistry at interfaces. An extension to this method is our recently developed Grand Canonical Molecular Dynamics (GCMD) approach that promises to merge reactive molecular dynamics (MD) simulations with constant composition of the system's components, enabling the simulation of catalytic processes near equilibrium [25].
In this work, selected nanoparticles in a hydrogen (H 2 ) and oxygen (O 2 ) containing environment have been investigated regarding the catalytic formation of water (H 2 O). In addition, the site-specific preference of the nanoparticles' surfaces-both pristine and after catalysis-towards electrophilic and electrophobic species has been calculated. Our GCMD method allows capturing oxygen-containing intermediate's surface coverages, structural rearrangements as well as roughening of the nanoparticles, both thermally activated and adsorption induced. Applying the gCN as geometrical descriptor for platinum surface atoms also produced observable trends in nanoparticle roughening. During the occurring surface rearrangement, the activity of water formation is classified according to the gCN of the platinum surface functional groups, thereby allowing an identification and allocation of active sites.

ReaxFF Reactive Force Field
ReaxFF is a versatile reactive force field approach being successfully applied to a wide range of chemical, physical and material science related research topics [21,23]. It relies on a bond-order-dependent formulation of the potential energy in conjunction with a time-dependent, polarizable charge description. The bond order is updated at every calculated iteration step depending on the local atomic environment, i.e. the respective bond lengths. As such, bond formation as well as bond breaking/dissociation can be captured and therefore entire reactions or chemical processes can be described as well. The ReaxFF potential involves both bonding terms (e.g. bond, angle and torsion contributions) and non-bonding interaction terms (e.g. van der Waals and Coulomb contributions). Atomic charges are derived from the electronegativity equalization method and serve for the description of electrostatic interactions [26]. Our self-developed Pt/O/H force field (details can be found in Refs. [27,28]) has been used for all simulations. As it has not only been trained and intensely tested against bulk platinum phases, oxygen adsorption and initial stages of platinum oxide formation, but also for energies and reaction barriers along various water formation mechanisms, this forcefield is aptly suited for the here investigated systems. Finally, coverages have been evaluated by analysing the bonding table during the simulation runs.

Grand Canonical Molecular Dynamics
Advantageous to simulations of time-dependent properties (such as processes on surfaces) is extending the canonical NVT ensemble (N = number of particles, V = system volume, T = temperature) by a grand canonical treatment. In our GCMD approach the composition as well as the quantities of the different components in the system are maintained by regularly adjusting the amount of H 2 and O 2 molecules (i.e. reactants) as well as H 2 O and H 2 O 2 molecules (i.e. products). In this way, a nearly constant level of partial pressures of H 2 and O 2 is conserved by adding molecules if necessary and removing product molecules (e.g. H 2 O and H 2 O 2 ) from the system. As chemical potentials are dependent on temperature and the respective activity, our GCMD approach aims at keeping the chemical potentials constant by firstly controlling the temperature, and secondly by adjusting the partial pressures or activities. Further information and details concerning this method can be found in Ref. [25].
As even ReaxFF simulations can only cover nanoseconds, applying a higher temperature (e.g. attributing a higher kinetic energy to the particles) is a convenient way to overcome higher activation barriers and to effectively sample reaction events. Thus, different temperature setups have been chosen for the present GCMD simulations, as summarized in Table 1: Firstly, the platinum structures have been appointed 300 K, whereas higher temperatures have been assigned to the H 2 and O 2 molecules of the surrounding atmosphere (T1-T3). These temperature setups have proven to accelerate the establishment of intermediate coverages on the nanoparticle structures, the observation of product formation and the induced structural rearrangement of the nanoparticles. Also, to compensate the high temperature of the gas phase molecules and hence avoiding an excessive preference of desorption processes, the platinum nanoparticles were kept at 300 K, thus providing a corresponding temperature basin. Secondly, temperature setups allocating the same temperatures to the nanoparticle and the H 2 and O 2 environment (T4-T5) have been applied for an overall comparison of the observed processes and reactions.
Prior to the GCMD simulations presented in the following, the nanoparticle systems in contact with a stochiometric 2:1-H 2 /O 2 atmosphere have been equilibrated (NVT ensemble) at 300 K for 125 ps. For the environment a relatively high pressure of 50 bar was used, which is justified by the characteristics and limitations of the theoretical approach: Only by increasing the pressure well above experimental conditions, sufficient particle densities could be achieved (36-120 H 2 and O 2 molecules depending on the system). It should be mentioned that under these simulation conditions, the H 2 /O 2 environment is still gaseous. After the initial equilibration, the subsequent GCMD simulations ran for 3 ns (a total of 12 million iterations). All simulations were performed such that after each 10,000 MD iterations, the system was analysed regarding the occurrence of reactant (i.e. H 2 , O 2 ) and product molecules (i.e. H 2 O, H 2 O 2 ). Based on this, H 2 and O 2 molecules were added to the system to maintain their respective partial pressure and present product molecules have been removed from the system, thus keeping the chemical potentials constant. During the reactive MD simulations presented in the following, a timestep of 0.25 fs was employed, utilizing a velocity-Verlet-algorithm. The simulation temperature has been controlled applying a Berendsen thermostat with a temperature-damping constant of 100 fs [29]. All simulations were performed employing ReaxFF implementation within the Amsterdam Density Functional (ADF) package [30,31]. For visualizations of the nanoparticle structures the program VMD has been utilized [32].

Nanoparticle System Setup
All nanoparticles have been cut from crystalline fcc-Pt using the optimized lattice constant of 3.95 Å. The cube consists of {100} facets, the octahedron possesses {111} and the dodecahedron {110} facets, whereas the cuboctahedron is made up of both {111} and {100} and the sphere consists of mixed facets. A diameter around 3-4 nm has been chosen as a compromise between the observed peak in mass activity at ~ 3 nm and to ensure for sufficient stability, as Pt dissolution is more pronounced on small nanoparticles below 4 nm [33,34]. Therefore, the diameter range between 3 and 4 nm should be catalytically most interesting [35]. The respective diameters of the investigated nanoparticles as well as the respective system sizes and compositions can be found in Table 2. The used ReaxFF force field's transferability from extended surfaces to nanoparticle structures has been confirmed recently [36]. There, the stability of octahedral, cubic, dodecahedral, cuboctahedral and spherical nanoparticles with sizes ranging between 1 and 10 nm has been tested, resulting in the stability sequence: cuboctahedron > octahedron > sphere > dodecahedron > cube. In addition, the formation energy of Pt clusters has been compared to corresponding DFT calculations [37], showing good agreement for the particle's size range studied here.

Results and Discussion
In the following section, we first characterize and discuss the investigated nanoparticles both in their pristine and roughened shapes, enabling the allocation of active sites for water formation in the subsequently performed GCMD simulations. Here, exemplarily only images regarding the cube and

Generalized Coordination Numbers (gCN)
For the characterization of the different nanoparticle shapes and their specific sites the generalized coordination number (gCN) turned out to be a robust and meaningful geometric descriptor. For all systems gCN i , being the generalized coordination number of a surface top site i, is calculated as the sum of all first-nearest neighbours' coordination numbers cn(j) divided by the bulk coordination cn max (e.g. 12 for face centred cubic lattices), see Eq.
(1). The respective coordination number, cn(j) , of each atom j is calculated as the sum of neighbouring platinum atoms expressing a bond order ≥ 0.3 to atom j.
Characterizing the nanoparticles with this gCN yields different values for adsorption sites, depending on their vicinity to vertices, edges and facets. For instance, the gCN of a top site on a {111} surface is calculated to be gCN i = (6·9 + 3·12)/12 = 7.5, summing up six neighbours in the first (each surrounding by 9 atoms) and three neighbours in the second layer (each surrounded by 12 atoms). Hereby, double counting of nearest neighbours is avoided and the gCN are between 0 and 12 for fcc crystals [10, 13]. With the gCN being able to distinguish between adsorption sites on surfaces, it allows estimating the preference of adsorbing species towards certain sites: e.g., hollow sites on the {111} can be separated in hcp (gCN = 7.5) and fcc (gCN = 6.95). Thus, the adsorption of electronegative species such as O* on fcc sites can be identified as adsorption at the three-fold hollow side with lower gCN value. As for the classical coordination numbers, a gCN value smaller than the bulk crystal value (e.g. 12 for fcc crystals) is typical for surfaces and describes atoms with a tendency towards the formation of new bonds, allowing to compensate for a lack in coordination. There exists a proportionality between the undercoordination and the inclination to form new bonds, as bond-order conservation theory suggests [10, [38][39][40][41]. It needs to be noted, however, that though some trends can be predicted the full description of adsorption behaviour still calls for more sensitive calculations or descriptors [10]. (1)

Solvent Accessible Surface Area (SASA)
To link the geometric arrangement of the nanoparticles to their estimated adsorption properties, we utilized a fast and efficient screening method. Therefore, we calculated the nanoparticle's solvent-accessible surface area (SASA) and performed a screening regarding the interaction towards atomic oxygen and hydrogen. Here, the interaction signifies (and is a qualitative measure for) the preference of the respective surface site towards an electrophilic (e.g. chargedonating, tested with the oxygen atom) or electrophobic (e.g. charge-accepting, tested with the hydrogen atom) species. Firstly, the SASA is obtained by distributing a set of analysis points along the van der Waals surface of the nanoparticle. Afterwards, at each point the interaction of a probe atom to the particle was evaluated by ReaxFF single point calculations, whereby the corresponding probe atom was placed at a distance of 0.5 Å from the nanoparticle's van der Waals-surfaces. 1

Connecting Both Descriptors (gCN and SASA)
In the following, these two descriptors for geometry and interaction were applied to the pristine, highly-symmetric nanoparticles.
Starting with the cube, the gCN of its surface as well as its preference towards electrophilic and electrophobic intermediates is depicted in Fig. 1. The lowest gCN values can be found on the cube's vertices (3.17) and on the edges (3.67). According to bond-order preservation theory [40], these sites are expected to express the highest preference towards charge-donating species (e.g. the oxygen probe atom). This is confirmed in Fig. 1b, showing the interaction with the oxygen probe atom.
Subsequently, the inverse relation between gCN and preference towards electrophilic particles can be observed in Fig. 1-that is decreasing gCN should result in increasing bond strengths [8,10]. Further, the preference of the Pt particle towards the hydrogen probe atom (being representative for a charge-accepting species) is more pronounced towards the nanoparticle's facets than to its edges and vertices.
For the octahedron, the increasing gCN values from the edges to the facets can be seen in Fig. S2.1 (see Supporting Information (SI)). The facets are even expressing the character of extended {111} surfaces, as becomes indicative by a gCN value of 7.5. As the vertices are displaying gCN values of 2.33, a strong inclination of bond formation to electrophilic species is expected, which is also confirmed by the vertices' red colouring in Fig. S2.1(b).
On the cuboctahedron, a broader variety of gCN numbers can be observed (see can also be concluded, that the underlying second layer is not accessible for adsorbates or intermediates, as the preference towards the oxygen and hydrogen probe atom is weak. The spherical nanoparticle consists of mixed facets and therefore also displays a wide variety of gCN values, see Fig. S2.4 (a). The atoms that are located at exposed sites are displaying the lowest gCN numbers as well as the highest preference towards electrophilic species, as determined by the oxygen probe atom.

Surface Roughening During Operation
Surface roughening of nanoparticle structures needs to be considered when accounting for the harsh catalytic conditions during the ORR [42]. While the roughening process or degradation may lead to numerous different particle structures, shapes and sizes, whose full analysis is out of the scope of molecular dynamics simulations, it is still possible to study the interplay between surface structure and intermediates. For instance, exposing the pristine, highly-symmetric nanoparticles to the hydrogen and oxygen-containing environment and the different temperature setups during our GCMD simulations, various effects of surface roughening could be observed on the investigated nanoparticle shapes. Again, the gCN as well as the preference of the surface sites towards a charge-donating or -accepting probe atom served as structural descriptors. The GCMD simulations have been performed for 3 ns, during which the distribution of the nanoparticle's gCN numbers have been determined within two timeframes. In the following, t 0 represents the starting point of the GCMD, Δt 1 covers the interval 0.75-1.5 ns and Δt 2 covers the second interval between 2.25 and 3.0 ns. The gCN distributions have then been averaged over the respective intervals and yield information about the time-dependent structural shape of the nanoparticles. As can be seen in Fig. 2a, the pristine cube (prior to the GCMD simulations, e.g. at t 0 ) features distinct gCN values. However, the clear attribution of the gCN to structural features becomes difficult already during the first time interval Δt 1 . Overall, one can see that the thermally-and/or adsorbate-induced local surface rearrangements cause a broadened distribution of gCNs.
Hereby, the peak corresponding to the gCN number of {100} facets (e.g. 7.0) shrinks considerably, while gCN values laying between the identified numbers of facets, vertices and edges appear. Remarkably, roughening of the cube's vertices also causes the appearance of gCN values ≤ 2, suggesting a strong preference towards charge-donating species. During the second time interval ( Δt 2 ) these low gCN values become even more frequent. Also, gCN values attributed to subsurface atoms become blurred. Along with the changes in gCN distribution, changes in the preference towards the oxygen and hydrogen probe atom are of interest when aiming at deeper insights into the Fig. 1 a Generalized coordination numbers of the pristine cubic nanoparticle. The gCN of the surface cover a range of 3.17 to 6.67. b Visual representation of the preference of the surface sites towards electrophilic (probing with an oxygen atom) or electrophobic (probing with a hydrogen atom) intermediates catalytic behaviour. As can be seen in Fig. 3, the exposed and strongly roughened vertices of the cube feature these low gCN values ≤ 2 as well as a stronger preference towards electrophilic and electrophobic species than the few remaining ordered facets. However, pit-shaped local arrangements (indicated by black arrows in Fig. 3b) also display an increased preference towards charge-donating or -accepting species. These shapes represent a concave site, already postulated to be more active towards the ORR by Calle-Vallejo et al. [15] The influence of the different temperature setups on the observed surface roughening is included in the SI exemplary for the cubic nanoparticle. As to be expected, the degree of the observed phenomena correlates with the applied temperature, however, the previously mentioned features subsist generally.
For the octahedron, the peak corresponding to gCN = 7.5 (i.e. the {111} facets) is reduced by ~ 90% until the end of the 3 ns long GCMD simulation (see Fig. 2b). This is mainly caused by adsorbate-induced surface buckling, where a platinum surface atom is pushed out of the surface plane to prevent oxygen-containing intermediates (e.g. OH*, OOH*, O*) laying too close to one another, thus reducing electrostatic repulsion [27]. The surface buckling observed on the {111} facets of the octahedron also creates sites with increased preference towards electrophilic species, see the exposed red atoms in Fig. S3.2(b). Apart from the buckling, the {111} facets are mostly maintained. As observed on the cube, roughening up of the edges causes the appearance of lower gCN numbers.
The {110} facets of the dodecahedron appear to be less stable than the {111} facets during the GCMD simulations: Structural rearrangements affect the anisotropic structure of the first and second layer and are thereby creating a lowsymmetry, rough surface (see Fig. S3.3). Concave sites (as already mentioned on the cubic nanoparticle) can also be observed on the dodecahedron, where the rearrangement of atoms belonging to the {110} facets near vertices is inclined to create pits.
The cuboctahedron, featuring both {111} and {100} facets expresses a similar behavior as observed on the cube and on the octahedron: The roughening of the edges invokes the occurrence of gCN ≤ 2, whereas the facets remain mainly intact, apart from local surface buckling (see Fig. S3.4). Those structural changes have resulted in exposed platinum atoms as well as pit-shaped sites near vertices displaying a strong preference towards charge-donating species.
On the spherical nanoparticle an almost even distribution of gCN values between 0.5 and 11.0 can be observed in Δt 1 (see Fig. S3.1). After 3 ns, only a minor increase of gCN values < 3.5 can be seen. Due to the sphere's mixed facet character, it features many pit-shaped arrangements of platinum surface atoms (see Fig. S3.5).
In general, all investigated nanoparticles express a broadened distribution of gCNs with tails both towards smaller values and higher values caused by surface roughening during the GCMD simulations. On all nanoparticles, the breaking up of the ordered nanoparticle structures in particular starts at vertices and to a lesser degree at edges. The more distinct the vertices (e.g. cube or octahedron) the more intense the roughening.

Coverage of Intermediates
Another feature that needs to be taken into consideration when studying the catalytic behaviour of structures is the coverage of intermediates and adsorbates during the processes, as the dynamically changing local atomic environment is presumably influencing the catalytic properties. However, these cannot be captured by descriptors solely focusing on a single binding energy per species. During our performed GCMD simulations, surface coverages of intermediates (e.g. O*, OH*, OOH*) indeed fluctuated over time and caused structural rearrangements, as previously discussed (e.g. surface roughening and buckling). In the following, the surface composition of the nanoparticles is briefly discussed. Hereby, to allow for a comparison between the different nanoparticles, the coverage is measured in number of adsorbate molecules per surface area (nm 2 ), see Table 3. The respective van der Waals surfaces of the roughened Pt nanoparticles have been obtained by the SASA method. 2 Starting with the cubic nanoparticle, OH* and O* show the highest coverages. Though the coverage of O* is higher than OH* (see Table 3), these O* are mainly incorporated into the roughened edges and vertices and should therefore be seen as surface oxide-species rather than adsorbates. OH* is covering the whole nanoparticle surface, on the exposed sites of strongly roughened vertices one can even observe multiple OH* per platinum atom. Additionally, OH* is encircling the pit-shaped arrangements.
On the octahedron, OH* is covering the nanoparticle in loose chains or ring-like structures. O 2 * is mainly adsorbed on the edges while O* is incorporated in the roughened sites. Also, OOH* is adsorbed on edges and the roughened vertices.
Similar observations can be made on the dodecahedron, with OH* forming chains on the topmost layers of the {110} facets and decorating exposed platinum atoms or island on the surface. No oxygen-containing intermediate is adsorbed on the second layer, as predicted by the SASA descriptor. O 2 * adsorption is more pronounced on the edges than on the facets.
Regarding the cuboctahedron, OH* is covering edges and facets likewise, on some exposed platinum surface atoms multiple OH* are adsorbed. O* is incorporated in the roughened vertices and edges. The same trends can be found for the sphere with O 2 * and OH* decorating the exposed sites while OH* is covering the whole nanoparticle surface. Also, O* is involved in the roughening process.
On all nanoparticle structures, OH* is the predominant intermediate during the GCMD simulations: Its coverage is steadily increasing, however, near the end of the simulations its increase slows. These mentioned observations are in their character observable in all applied temperature setups; however, their degree differs.

Activity of Water Formation
Also, the activity of the nanoparticle differs clearly, as can be seen in Table 4. To allow for a comparison between the different nanoparticles, the activity is measured as number of produced H 2 O molecules per ns and nm 2 . The cubic nanoparticle appears to be the most active one towards water formation, showing the highest value of H 2 O production. The cuboctahedron and sphere perform moderate, expressing These observations can only carefully be compared to experimental findings. Our simplifications made (e.g. the lack of electrolyte) complicates the rating of the nanoparticle's activity: Sánchez-Sánchez et al. [43] discovered the dependence of the {111} and {100} facet activity of the anion present in solution. Other computational approaches, like the density functional theory study performed by Mahata et al. [44], investigating the stability, ORR mechanisms and kinetic aspects, attribute the cuboctahedral nanoparticle a higher activity and selectivity compared to the octahedron. This is mirrored in our observed activities.

Allocation of Water Formation Sites
After applying the geometric descriptors for the characterization of the roughened nanoparticle structures discussed in the previous section, the nanoparticle's actual catalytic behavior was then investigated during the GCMD simulations. Here, the activity of water formation has been classified according to the gCN of the platinum atoms, thereby allowing identifications of active sites. The produced H 2 O molecules have been traced during the MD simulation regarding their metallic binding partners, whose respective gCN values have been calculated. Thereby, the H 2 O production could be allocated to certain gCNs. This evaluation has been performed separately during the time frames Δt 1 (0.75-1.5 ns) and Δt 2 (2.25-3 ns) of the GCMD simulations. For all nanoparticles, three GCMD simulations have been performed at each temperature setup T1-T5. Then, the formation of the water molecules has been traced in all simulations. The combined results of all simulations of one nanoparticle (for the respective time frame) were used to generate the respective locally resolved activity diagram. These diagrams show the number of produced H 2 O molecules per gCN number (scaled by its frequency on the nanoparticle) and allow the allocation and comparison of active sites regarding water formation.
For the cubic nanoparticle, the locally resolved activity diagram is displayed in Fig. 4a. Two distinct trends can be observed: First, water formation is performed rather at the cube's edges and vertices than at facets (as the gCN would suggest), though those sites have already undergone slight structural rearrangements. Second, with increasing roughness of the surface (evaluating Δt 2 ) the H 2 O formation preferentially occurs at surface sites characterized by even lower gCN values ( ≤ 2). Such low gCN values can be found at the exposed, low-symmetric sites of the cube (as depicted in Fig. 3a).
For the octahedron and cuboctahedron, the sites featuring gCN values in the range 1.0-2.5 are also the most active with respect to H 2 O formation. On the octahedron (see Fig. 3b), sites with low gCN values are formed due to surface buckling on the {111} facets, therefore, these buckled platinum atoms seem to be active for water formation as well. Additionally, for gCN values around 7 the resolved water formation activity is slightly enhanced with increasing simulation time: These sites can be attributed to pit-shaped formations, emerging due to the roughening processes, or intact {111} facets. As the roughening of the nanoparticle increases with time, it is more probable for the pit-shaped formations to express the water forming activity as the {111} facets are expected to lose their integrity. On the cuboctahedron (see Fig. S5.1(b)), the water formation activity can be mainly attributed to roughened edges, as those are expressing the lowest gCN values.
However, regarding the dodecahedron, no clear preferences can be observed: Sites characterized by gCN values between 0.5 and 7 contribute to the product formation. The same behavior can be observed for the sphere (see  Fig. S5.1(c)), though here a tendency towards lower gCNs becomes more pronounced when evaluating the second time frame Δt 2 , compared to Δt 1 . Comparing all investigated nanoparticles regarding their resolved local water forming activity, two general observations can be made: Active sites towards H 2 O formation are on the one hand exposed platinum surface atoms, being found on the roughened edges and vertices of the nanoparticles or on buckled surface atoms. These exposed sites feature low gCN values and are decorated with oxygen-containing intermediates (e.g. OH*). Second, pit-shaped sites being characterized by gCN values around 7, mostly located near roughened edges or vertices, also pose active water forming sites. As even the GCMD can only capture a few nanoseconds, the formation of the pit-shaped sites may become more important due to an enhanced surface roughening with time.
These findings are in line with conclusions drawn from experiments: The importance of investigating the effects of heat treatment on platinum nanoparticles has been stressed by Aarons et al. [45], observing surface roughening on Pt nanoparticles going along with a general lowering of the surface gCN. This is further supported by Chung et al. [46], attributing the nanoparticle surficial rearrangement a beneficial contribution to the ORR activity, the degree of which being affected by the heat treatment conditions.

Conclusion
Transferring structural properties to catalytic behaviour plays a key role in designing catalytic systems. In this work, we have illustrated how gas-phase platinum nanoparticles are modified in their surface morphology by the catalytic conditions and based on this, how these rough, low-symmetry nanoparticles behave regarding water formation. Hereby, we have taken into consideration the influence of oxygen containing intermediates on the surfaces. We could observe that edges and vertices are prone to surface roughening, while {111} facets are experiencing surface buckling. These exposed, low-coordinated sites have been detected as the most active regions towards water formation. Another feature, caused by surface roughening, is the occurrence of pit-shaped sites. While being characterized by gCN values around 7, these also express an increased preference towards electrophilic and electrophobic species and pose another region active for water formation. Moreover, the effects associated with adsorbate-induced or thermally-activated processes causing a broadening of the gCN distribution (i.e. transferable to a roughening of the nanoparticle structures), are already occurring at temperatures of 500-600 K. Therefore, it is of utter importance to further elucidate the properties and catalytic behaviour of realistic platinum nanoparticles. Hereby, the theoretical-computational approaches can be envisioned to further enhance the understanding of those catalysts.