Shape-Dependent Single-Electron Levels for Au Nanoparticles

The shape of metal nanoparticles has a crucial role in their performance in heterogeneous catalysis as well as photocatalysis. We propose a method of determining the shape of nanoparticles based on measurements of single-electron quantum levels. We first consider nanoparticles in two shapes of high symmetry: cube and sphere. We then focus on Au nanoparticles in three characteristic shapes that can be found in metal/inorganic or metal/organic compounds routinely used in catalysis and photocatalysis. We describe the methodology we use to solve the Schrödinger equation for arbitrary nanoparticle shape. The method gives results that agree well with analytical solutions for the high-symmetry shapes. When we apply our method in realistic gold nanoparticle models, which are obtained from Wulff construction based on first principles calculations, the single-electron levels and their density of states exhibit distinct shape-dependent features. Results for clean-surface nanoparticles are closer to those for cubic particles, while CO-covered nanoparticles have energy levels close to those of a sphere. Thiolate-covered nanoparticles with multifaceted polyhedral shape have distinct levels that are in between those for sphere and cube. We discuss how shape-dependent electronic structure features could be identified in experiments and thus guide catalyst design.


Introduction
Although most properties of nanoparticle-containing materials are governed by nanoparticle size, the importance of nanoparticle shape is increasingly appreciated for applications, most notably in optoelectronics and catalysis. Consequently, in addition to controlling size, considerable effort has recently been devoted to determining and tailoring the exact atomic arrangement in nanoparticles. The shape of nanoparticles plays a crucial role in the efficiency of many types of catalysts, such as transition metal heterogeneous catalysts [1][2][3] and titanium oxide photocatalysts [4,5]. For example, for rutile-supported Au, the supporting oxide is found to enhance the catalytic properties of the metal by enforcing specific shapes on the metal clusters, and also by providing a reaction path at the metal oxide interface [6]. The shape of Au nanoparticles can be easily modified during synthesis by means of surfactants and other experimental conditions [7][8][9][10].
The precise characterization of the shape of a metal nanoparticle in a sample is not a trivial task. Despite enormous progress and widespread use during the last few decades, atomic-scale imaging techniques now described in textbooks [11], such as optical, electron, scanning tunneling and atomic force microspcopy, present technical difficulties and inherent resolution limitations. Accurate atomic-level nanostructure determination requires state-of the-art methods which often combine experimental techniques and nanoparticle modeling [12,13]. On the other hand, the quantum energies in metal nanoparticles are clearly affected by their shape [14] and the presence of adsorbates, as shown both theoretically [15] and experimentally [16]. Here we propose an alternative method of shape characterization: probing the single-electron energy levels of the nanoparticles.
In addition to shape-probing, the study of electron confinement in nanoparticles has applications in many other fields of modern materials science. Quantization phenomena are important even in large nanoparticles, where single-electron discrete energy level spacing is on the order of a few meV. There are optical and electrical properties which are determined by quantized electronic states and can be used to probe nanoparticle size and shape. Optical measurements, such as absorption spectra, that provide discrete levels for semiconductor nanoparticles, might be difficult to apply to metal nanoparticles, where plasmon resonance effects dominate in optical frequencies [17,18]. This is especially the case with gold, where bulk plasmon frequencies in the UV shift in the range of visible light for nanoparticles. Interaction of metal nanoparticles with light depends on their composition, size, and shape, as well as on the surrounding medium. Extensive work over the past century, starting with light scattering and absorption from spherical particles [19] and developed in the field of plasmonics, has led to the complete understanding of light nanoparticle interaction and to diverse practical applications [20]. Resonances, local-field enhancement, and absorption are used in optoelectronic devices [21,22]; in photovoltaics, photo-sensitive materials include nanoparticles tuned so that the peak of their absorption coefficient is close to the maximum of the solar spectrum [23]; in the paint industry, nanoparticles give specific colors to colloidal suspensions [24]. Theoretical results connect optical data with the size and the overall nanoparticle shape characteristics [25,26], but not with the details of the atomic arrangement. Electrical measurements are another possibility for investigating metal nanoparticle properties. Discrete electronic levels are probed by single-electron tunneling spectroscopy with single electron transistors (SET) [27][28][29][30].
In this work, we present results for the electronic structure of realistic gold nanoparticles. The nanoparticle models were obtained by an atomistic Wulff construction in inert and reactive environment using density functional theory (DFT) results for the surface energies of different facets. The shape of metal nanoparticles at equilibrium depends on the surface energies of the various (hkl) crystalline planes of the material. Strong interactions between the metallic particle and the surrounding inorganic or organic matter result in significant changes of its shape, which in turn modifies its electronic, optical, and catalytic properties. We show how discrete energy levels and the corresponding density of states are related to nanoparticle shape so that knowledge about nanoparticle electronic structure can be used for shape determination. SET experiments with gold nanoparticles indicate that the independent electron approximation is valid [28]. The eigenenergies of the Schrödinger equation for a free electron moving inside the nanoparticle will be dependent on the nanoparticle shape. Reversing the problem, if one measures a series of electron energies for a given material, they can deduce what nanoparticle shape corresponds to this distribution of energies.
The task in front of us is similar to the problem set by mathematician Mark Kac fifty years ago: Can one hear the shape of a drum [31,32]? In other words, can one deduce the shape of a two-dimensional drum from the vibrational eigenvalues of its membrane? The Schrödinger equation for a single free electron in a nanoparticle is mathematically the same as the Helmholtz equation that gives the eigenfrequencies of vibration, the only difference being the number of dimensions (three vs. two). There is no clear answer to this question in three dimensions, although this is still an active research topic-see reference [33] for another application in nanoparticles.

Single-Electron Levels for Typical Nanoparticle Shapes
The task before us is to find stationary states of an electron that is free to move inside a nanoparticle. No force is acting on the electron, and the only constraint in its motion is that it cannot escape the confinement space, and therefore its wavefunction will be zero along the planes that define the boundary of the nanoparticle. The energy levels corresponding to stationary states for this free electron can be found by solving the time-independent Schrödinger equation: where E i and ψ i (x, y, z) are the i th eigenvalue and eigenstate of the system, respectively; −¯h 2 2m ∇ 2 is the kinetic energy of a particle with mass m; and V(x, y, z) is the potential energy. In our case, potential energy is defined as: The potential described in Equation (2) is the simplest possible confining potential for a metal nanoparticle. In reality, the potential will have a finite depth, which will be of the order of the work functions of the various (hkl) surfaces contained in the nanoparticle. This workfunction will in turn be modified by the presence of adsorbates, as adsorbates will modify the electronic structure, in particular for smaller nanoparticles [16]. The presence of the support might result in modifications of the nanoparticle shape [34]; charge transfer to Au nanoparticles from oxides such as rutile TiO 2 is minimal [35]. Moreover, a realistic nanoparticle will contain structural and chemical defects that will also affect its electronic structure. All these considerations make the problem extremely complicated (and at the same time interesting), making it very difficult to study.
We use a value of 100 a.u. (about 2700 eV) as an approximation of infinity in Equation (2). We tried different values of the confining potential, but found no qualitative difference in results, such as symmetry and degeneracy of the eigenstates or the shape-dependence of the excitation energies. For this reason, we proceed with the simplest possible potential that guarantees that the wave functions are zero outside the nanoparticle. The infinite-potential well of Equation (2) is the simplest possible mathematical model that reproduces the shape-dependence of the excitation energies for nanoparticles of the same volume and the same symmetry point group.
In Section 3, we present the methodology we followed to solve Equation (1). In the following subsections, we describe our findings for five different nanoparticle shapes: 1. Cube; 2. Sphere; 3. Truncated octahedron, formed when Au has very weak interactions with its environment, such as Au nanoparticles bonded to C nanotubes [36]; 4. Complex rounded shape, formed when Au is in equilibrium with small molecules with strong binding to Au, such as CO [37]; 5. Deltoidal icositetrahedron, formed when Au is in equilibrium with large molecules with strong binding to Au, such as alkanethiols [15,38].
The first shape is used as the basis for all other calculations, whereas the second one is used to verify the accuracy of our numerical method. The last three shapes are typical equilibrium shapes for Au nanoparticles.
We only consider shapes that can be formed from face-centered cubic (fcc) Au under thermodynamic equilibrium, as these are the shapes most commonly seen in experiments. Shapes considered here are invariant under all transformations of the fcc lattice space group, O h . The effect of major shape changes, such as from spheres to ellipsoids, to electron and plasmon energies is reviewed in reference [14].

Cubic Nanoparticle
The solution of the three-dimensional "particle-in-a-box" problem with cubic symmetry can be found in any Quantum-Mechanics textbook. The eigefunctions of the Hamiltonian are where a is the length of the cube; and n, m, and l are positive integers, the so-called "quantum numbers" of the problem. Being eigenfunctions of a hermitian operator, the φ n,m,l (x, y, z) functions form a complete and orthonormal set of functions that are vanishing in the borders of the box. The eigenvalues of the Hamiltonian, E n,m,l , are given by: where,h is the reduced Planck constant; and m e is the mass of the electron. The ground state is non-degenerate (all three quantum numbers are equal to one). The next three eigenenergies are triply degenerate. These three triplet states have quantum numbers equal to (n, l, m) = (1, 1, 2), (1, 2, 2), (1, 1, 3), respectively.
We use these eigenfunctions and eigenvalues as the basis for the solution of the Schrödinger equation in other nanoparticle shapes. The details of the calculation are given in Section 3.

Spherical Nanoparticle
The problem of an electron that moves freely in a sphere of radius R can also be solved analytically. The energy levels, E nl , are given by: where χ nl is the n th root of the l th order spherical Bessel functions j l (x). The corresponding eigenfunctions are where j l (x) is the spherical Bessel function of order l and Y lm (θ, φ) is a spherical harmonic. k nl is related to the energy through E nl =¯h 2 k 2 nl 2m e . The symbols n and l represent the principle quantum number and the angular momentum quantum number of each state, respectively. As usual, the third quantum number, m, is smaller in magnitude than l and does not enter Equation (5). This is the reason that each one of the energy levels is (2l + 1)-fold degenerate. In contrast to what is found in the solution of the Schrödinger equation for the hydrogen atom, here l is not necessarily smaller than n. For the spherical nanoparticle, the ground state is again non-degenerate and the first excited state is a triplet. A key difference between spherical and cubic particles occurs in the next six states: we found two triplets for the cube, whereas a quintuplet and a singlet were found for the sphere.
To validate our methodology, we compare the energy levels from Equation (5) to those found by numerical solution of an atomistic model for a spherical nanoparticle. We constructed a spherical particle of volume V = 481.42 nm 3 and sphericity equal to 99.6%. The sphericity of a particle is a measure of how close its shape is to a sphere. It is defined as π 1/3 (6V) 2/3 /A, where V is the volume and A is the area. A sphere has sphericity 1 while a cube has sphericity of about 0.8.
The numerical values are calculated using a set of 32,768 basis functions, with the infinite potential outside the nanoparticle defined as V 0 = 100 a.u. (i.e., about 2700 eV). Both the number of basis functions and the value of V 0 , are chosen so that the ground state and the 100th eigenenergy are well converged. In Table 1, we present the theoretical and calculated energies of the thirty lowest-energy eigenstates for a spherical nanoparticle with volume V = 481.42 nm 3 . The eigenvalues of the analytical solution for different l and n are put together and sorted; we use a new index, i, to count them. The calculated eigenenergies were found to be slightly overestimated compared to the analytically obtained ones. The relative error was 4% or less. The simulated nanoparticle, being a collection of Au atoms in the fcc lattice, is not a perfect sphere; any vertices tend to increase the kinetic energy of the electron. Another factor that affects the agreement between the analytical solution and the simulation could be sought in the vague definition of "volume" for a nanoparticle. Here, for simplicity, we consider that the boundary of the nanoparticle consists of planes that pass through the nuclei of the outermost atoms. Due to the quantum nature of the electron, however, the actual volume that is available will be a bit larger. So, the simulated electron is a bit more confined and hence its energy is a bit higher.
Despite these small and expected discrepancies, overall the quality of the simulation is satisfactory: differences between energy levels are correct to the first decimal, and degeneracies are correctly reproduced. Moreover, the qualitative features of the wave functions are well represented, including the symmetries of the eigenstates.
In Figure 1, we present the probability function, |ψ| 2 , of the eigenfunctions for the first four calculated energy levels. In the analytical solution, the four lowest energy eigenvalues correspond to a singlet (n = 1, l = 0), followed by a triplet (n = 1, l = 1), a quintet (n = 1, l = 2), and a singlet (n = 2, l = 0). The simulations correctly yield all these degeneracies. The density plots shown in Figure 1 correctly resemble the familiar shapes of s, p, d, and s * orbitals, which are the real-valued linear combinations of the spherical harmonics with l = 0, l = 1, and l = 2.
As a last test of the method, we compare the density of states per volume of an electron in a spherical nanoparticle. The density of states (DOS) per volume, ρ V (E), is defined as: where V is the volume of the nanoparticle and the sum runs over all eigenstates of the Hamiltonian. The delta function can be approximated using a Lorentzian distribution with broadening σ: In Figure 2, we plot the ρ V (E) of an electron in a spherical nanoparticle for two sets of eigenenergies: those obtained from our simulation and the analytical values of them using Equation (5). We used two different broadenings for the Lorentzian approximation of the delta function in Equation (8). To make our plot size-independent, we multiply our energies with the volume to the power of 2/3. To change our energy units to ordinary energy values, one needs to divide the values of EV 2/3 by V 2/3 , that is, the volume of the nanoparticle to the power of 2/3. For example, a value of 10 eV nm 2 in the x-axis of Figure 2 corresponds to 38 meV and 154 meV for a spherical nanoparticles with radii of 10 nm and 5 nm, respectively.
The calculated DOS matches the analytically obtained one very well. The first ten peaks are at the same energies, while even for quite large values the difference is small. This difference can be attributed to the underestimation of volume in the atomistic nanoparticle and is small. For instance, the value of 60 eV nm 2 , where a significant discrepancy between the two functions is observed in Figure 2b, corresponds to energy of the order of 1 eV for a nanoparticle of radius 5 nm. For comparison, the difference between the ground-and the first excited state for this nanoparticle is only 15 meV. Therefore, our method provides reliable DOS for all practical purposes.

Single Electron in Equilibrium-Shaped Au Nanoparticles
The shape of Au nanoparticles depends on the concentration and strength of any bonds that are formed between the Au atoms at the surface of the nanoparticle and any ligands or molecules of matter that surround that nanoparticle. Formally, the shape is obtained by means of a Wulff construction where the surface tensions, γ, are substituted by the interfacial tensions, γ int , between Au surfaces and their environment. A simple relationship between the two has been found when there is little interaction between adsorbates [39]: In this equation, θ is the coverage of Au surfaces (θ = 1 corresponds to full coverage) while E ads is the adsorption energy of ligands (exothermic adsorption means negative E ads ). A at is the area per Au atom, that is, the total area divided by the number of surface Au atoms. According to Equation (9), γ int ≈ γ when the coverage, θ, or the adsorption energy, E ads , is small. For this reason, Au nanoparticles are found to have similar shapes in a variety of materials.
We begin by considering an atomistic Wulff construction for Au with no interactions with its environment. The same shape is found for cases where adsorbates exist on Au surfaces, but the ratio of the binding energy over the area of Au is smaller than about 5% of the surface tension of Au [39]. This shape is found to be a truncated octahedron (see Figure 3). The Wulff construction agrees very well with microscopy images of Au nanoparticles in a variety of materials [36,40,41]. In the following, we will use the term "clean" to describe this shape. A set of 32,768 basis functions was used to calculate the eigenenergies and eigenfunctions of the nanoparticle shown in Figure 3 (left). The diameter of this nanoparticle is d = 10.45 nm, it contains 28,806 atoms, its volume is 485.38 nm 3 , its sphericity equals 89.87%. The density of step-edge atoms that typically act as active sites for catalysis is 72 µmol/g.
In Figure 4, we present the probability function, |ψ| 2 , of the eigenstates for the first four eigenenergies in the case of a clean nanoparticle. Apart from the obvious change of the nanoparticle boundary, the energy degeneracy and the symmetry of the wavefunctions for the two lowest eigenenergies is the same as for the cubic nanoparticle. The ground state is non-degenerate and has an eigenfunction with the same sign everywhere. The first excited state, which is triply degenerate, has three eigenfunctions that change sign once and are antisymmetric with respect to the three main mirror planes of the nanoparticle. The next energy level, however, differs from both the cubic nanoparticle (where the third lowest eigenenergy is a triplet) and the spherical one (where the third lowest energy is a quintet). Here, due to the absence of spherical symmetry, the quintet of the l = 2 states breaks into two sets: first, a doubly degenerate state with eigenfunctions similar to the d z 2 orbital; second, a triply degenerate state with eigenfunctions similar to d xy , d yz , and d zx orbitals.

Au Nanoparticles with Strong Interactions with Their Surroundings
The shape of Au nanoparticles changes dramatically when Au has strong interactions with materials around it. However, only very strong binders and/or molecules that can bind Au at high areal densities are capable of modifying the equilibrium shape of Au particles. Here, we consider two characteristic examples: first, the shape of Au at equlibrium with CO gas at CO oxidation conditions, taken from reference [39], and second, the shape of Au at equilibrium with methanethiol molecules at a density typical for Au encapsulated in polymer matrices, taken from reference [42]. The simulated nanoparticles are shown in Figure 3.
The CO-covered one contains 27,981 atoms, its diameter is d = 10.03 nm, its volume is V = 500.36 nm 3 , its sphericity is 98.6%, and its density of active sites is 181 µmol/g. The thiolate-covered nanoparticle has a polyhedral shape close to that of a deltoidal icositetrahedron. A particle with 28,502 atoms, d = 10.87 nm, V = 505.68 nm 3 , sphericity of 95.3%, and density of active sites 273 µmol/g was simulated.
In Figure 5, we present the probability function, |ψ| 2 , of the eigenstates for the lowest single-electron eigenenergies for the thiolate-covered gold nanoparticle. As was the case with a clean nanoparticle, the symmetry of the eigenstates lies in between the symmetry of states in a sphere and the symmetry of states in a cube. Again, the fivefold d state of the sphere has split into a triplet and a doublet. The main qualitative difference between protected and clean nanoparticles is that in the former, the doublet and triplet states are very close in energy, making these five states almost degenerate. The geometric characteristics of the eigenfunctions are the same as those for the clean Au nanoparticle.

Shape-Dependent Electronic Structure
We now compare our findings for the different nanoparticle shapes, trying to address whether is could be possible to probe the nanoparticle shape, and therefore its catalytic activity, without microscopy. In the following, we summarize the key differences in the electronic states of the shapes considered in this work and discuss possible connections to state-of-the-art experimental techniques that could probe these electronic states.
In Figure 6, we present the density of states for all shapes considered in this study. In order to eliminate size-dependent properties, we use dimensional analysis to find that the density of states of a free electron is proportional to the volume while its energy is proportional to the volume to the power of − 2 3 . The x-axis can be transformed to usual energy units by multiplying with V −2/3 , where V is the volume of the nanoparticle; see Section 2.1.2 for examples.
Interestingly, the shape with the highest sphericity, which is the CO-covered nanoparticle, has a DOS that is very close to that of a sphere throughout the range of energies considered. The clean nanoparticle has features that resemble more those of the cube than those of the sphere. Finally, the SCH 3 -covered nanoparticle has distinct features in its DOS that make it distinguishable from both the sphere and the cube. Let us now focus on the first few excited states and the corresponding excitation energies of the electron, which might be a more reliable probe of the shape. In Table 2 we present the ten lowest eigenvalues of the Schrödinger equation, taking the ground state energy to be zero in every case. The nanoparticles are ranked according to the first excitation energy. Table 2.
Energy eigenvalues for five different nanoparticle shapes. The scaled energy, A clean nanoparticle is found to have its first excitation energy 3% higher than the nanoparticles covered with CO or thiolates. Moreover, the faceted shape can be clearly differentiated from the spherical or the cubic shapes. Even the difference between CO-and thiolate-covered particles, which is only 1.9%, is undoubtedly within the resolution of modern state-of-the-art instruments.
In the second excitation energy, the differences are more prominent. A sphere has a five-fold degenerate state; cube and CO-covered nanoparticles have three-fold degenerate states; thiol-covered and clean nanoparticles have doublet states. Different degeneracies yield different peak heights in absorption experiments. For the second excitation energy, peak height might be a more important differentiation tool than peak position, as CO-and SCH 3 -covered nanoparticles have very similar energies: Energy times volume to the 2/3 power equals 2.402 for CO-covered and 2.388 eV nm 2 for SCH 3 -covered nanoparticles (see Table 2).
These differences in excitation energies should in principle be probed by an appropriate experiment. It might be a great challenge to design an optical experiment to distinguish these single-electron levels through the great background of light scattering and absorption. However, recent developments in optical trapping allow for detailed studies of light-matter interaction [43]. In such state-of-the-art optics, one could in principle differentiate between plasmonic and electronic peaks, as the latter will be much narrower. Another possibility for experimental observation of these states would be through tunnelling experiments and other electron-transfer measurements, such as Scanning Tunneling Spectroscopy (STS) or Electron Energy Loss Spectroscopy (EELS). Recent examples of probing single-electron states can be seen in Refs. [16,29,30]. Particularly for the lowest-lying electron states, X-ray techniques might be able to probe them-for instance, X-ray absorption spectroscopy can probe states lying hundreds of eVs below the Fermi level with resolution of 100 meV; see for example [44,45].
If such a method exists, it will allow for quick and reliable characterization of catalysts, as it can provide a direct link between single-electron states and the density of active sites for catalysis. For example, when a Au nanoparticle with a volume of 500 nm 3 is covered by thiols, the first excitation energy drops from 18 to 17 meV while the tenth one rises from 60 meV to 64 meV. These small changes in single-elctron energies correspond to a huge change in the density of active sites, which rises from 72 µmol/g for a clean nanoparicle to 273 µmol/g for a thiol-covered one. Probing single-electron states might be a very accurate tool for catalyst characterization.

Atomistic Wulff Construction
The Au nanoparticles used in this study are generated by means of our atomistic-Wulffconstruction method. The Wulff theorem states that the equilibrium shape of a crystal is a polyhedron. The distance of each (hkl) face of this polyhedron from the center is proportional to the surface tension of the (hkl) surface of the material [46][47][48].
Wulff construction has been a valuable tool for nanoparticle shapes, and has been used for nanoparticles formed by a variety of materials, including supported Au [34,49], diamond [50], TiO 2 [51], Si in amorphous SiO 2 [52], diamond in amorphous C [53], Rh and Pd in oxidizing conditions [54], Cu in N gas [55], Au in oxidizing conditions [56], noble metals with an environment [57], complex metal hydrides [58], iron carbides [59], and dawsonites [60,61], just to name a few. The atomistic Wulff construction takes it one step further and fills the Wulff polyhedron with atoms, taking into account the finite size effects at the nanoscale. It has been used successfully in the context of theoretical catalysis for nanoparticles made of Ru [62][63][64], Ag [65], and Au [39,42,61,66,67].
The atomistic models of the nanoparticles used in this work are constructed using the methodology of references [39,42]; for a review of this method, see reference [68]. The stability of these nanoparticles has been checked in a previous work of the authors by means of Molecular Dynamics (MD) simulations using a many-body Effective-Medium Theory (EMT) potential [69].

Numerical Solution of Schrödinger's Equation
The difficulty in solving the eigenvalue problem defined in Equations (1) and (2) lies in the nature of the potential. In the case of faceted nanoparticles, like those considered here, the potential V(x, y, z) in Equation (2) is a nonseparable one. Furthermore, the eigenstates of the electron must be vanishing in a very difficult boundary condition, i.e., the borders of the nanoparticle.
We overcome these difficulties by using a suitable set of basis functions and expand wave functions on this basis. The choice of the set of the basis functions is crucial, since it has to be a complete and orthonormal set. An extra advantage of a good basis set is that if it is appropriately related to the problem to be solved, good accuracy is achieved with a limited number of basis functions. Here, the basis functions we use are the solutions for the cubic nanoparticle, φ j (x, y, z). Each eigenfunction of the system, ψ i (x, y, z), is expanded as a linear combination of basis functions as follows: where c i,j is a scalar coefficient describing the weight of the contribution of the j th function of the base to the eigenfunction i. Using a set of basis functions in order to numerically solve the eigenvalue problem, one has firstly to construct the matrix of the Hamiltonian of the system. Then, to diagonalize this matrix, commonly using Linear Algebra PACKage (LAPACK) subroutines in order to get the eigenenergies and the c i,j coefficients.
We start by finding the bounding box for each nanoparticle. This is a rectangle with dimensions (a, b, c), as shown in Figure 7 (left). We then use the eigenstates of the well-known "particle-in-a-box" problem as a set of basis functions:  These are solutions of the Schrödinger equation for the potential energy V(x, y, z) defined as: The advantage of these eigenfunctions is that they are also eigenstates of the kinetic energy operator of our system. We then write each matrix element of the Hamiltonian of the system, (φ n ,m ,l , Hφ n,m,l ), as: (φ n ,m ,l , Hφ n,m,l ) = E n,m,l δ n ,n δ m ,m δ l ,l + (φ n ,m ,l , V(x, y, z)φ n,m,l ) (13) where (φ n ,m ,l , V(x, y, z)φ n,m,l ) is the expectation value of the potential energy for the potential defined in Equation (2). The problem of calculating the (φ n ,m ,l , Hφ n,m,l ) is now reduced to calculating the (φ n ,m ,l , V(x, y, z)φ n,m,l ) integral. We substitute the infinity value of the potential with a large enough value, V 0 . The matrix element then takes the following form: (φ n ,m ,l , V(x, y, z)φ n,m,l ) = V 0 φ * n ,m ,l φ n,m,l dV, where the integration is over all volume expept fot the interior of the nanoparticle. The nonzero integral of V(x, y, z) outside the nanoparticle can be calculated either numerically, using for example a three-dimensional trapezoidal rule, or analytically. The disadvantage of a numerical solution is that it needs a high-density grid in order to be as adaptive as possible to the boundary condition of the potential and also to be able to catch the fluctuations of the functions of the basis. This high-density grid makes the calculation too slowly and not accurately enough.
In this study, both numerical and analytical techniques were performed. In the following, we describe in detail the method used to calculate the expectation value of the potential energy using an analytical integration method.
Due to symmetry, we can use only the region where: 0 ≤ x ≤ a 2 , 0 ≤ y ≤ b 2 , and 0 ≤ z ≤ c 2 , as shown in Figure 7 (right). The integration volume is spilt into a sum of regions. Each region is defined by corner atoms of the nanoparticle and the space between the nanoparticle and the box. After defining these regions, we can solve the potential integral, V region , in each of them using the form: Using the Euler's formula, we rewrite each sine with exponential functions: and thus Equation (15) is written as: Expanding this integral we obtain that V region is given by: The above integral can be reduced using the fact that its 64 terms are all of the form: where s(t) is the sign of each term; and N(t), M(t), and L(t) represent each one of the 64 sets of combinations of (±n , ±n), (±m , ±m), and (±l , ±l). Substituting N(t), M(t) and L(t) into V region and taking into account that the box used is a cubic one (a = b = c), Equation (18) can be written as a sum over 64 terms in form: x e + iπM(t) a y e + iπL(t) a z dz dy dx (20) With limits defined by the borders of each region given as: This triple integral of each one of these terms can be calculated analytically and is given by: There are a few cases where one or more of the denominators of Equation (21) are equal to zero. For example, when L(t) = 0, all of the denominators are equal to zero. In such a case, it is easy to show that the nominator of the same fraction also equals zero and thus the fraction takes the indeterminate form of 0/0. This indeterminate form can be addressed either by applying the l'Hopital rule, or (more easily) by defining the integral again without the terms that are equal to zero and solving it from scratch. In the code this is implemented using if statements, where each one uses a formula to return the value of the integral for each specific case. Equation (21) was used in order to calculate the V region of each one of the regions.
Finally, the expectation value for the potential energy is given by the sum of V region over all the regions as: (φ n ,m ,l , V(x, y, z)φ n,m,l ) = ∑ regions V region (22) With the expectation values of the kinetic and the potential energy known, we construct the Hamiltonian matrix. This matrix is a real symmetric one and can be diagonalized using the TRED2 and TQL2 subroutines from the Fortran LAPACK library. The former reduces the real symmetric matrix of the Hamiltonian to a symmetric tridiagonal one, and the latter finds the eigenvalues and eigenfunctions of the reduced matrix. These eigenvalues and eigenfunctions are the eigenenergies and eigenstates of a free electron in the nanoparticle.

Conclusions
We presented eigenenergies and eigenstates as well as the density of states for several distinct shapes of gold nanoparticles. We find that both the absolute values of eigenenergies and their degeneracies depend crucially on the shape of the nanoparticle. The density of states clearly shows a distinguishable shape-dependent electronic structure. The first excitation energy alone can be used as a safe criterion for the determination of the shape of Au nanoparticles, which in turn can yield the density of active sites for catalysis. We hope that this work might be the first step towards the identification of nanoparticle shapes using single-electron levels. As such, it might be a step towards catalyst design as an extra characterization tool for catalysts by means of the energy levels of their electrons.