Next Article in Journal
Evaluation of Functionalized Porous Titanium Implants for Enhancing Angiogenesis in Vitro
Next Article in Special Issue
Using Central Composite Experimental Design to Optimize the Degradation of Tylosin from Aqueous Solution by Photo-Fenton Reaction
Previous Article in Journal
Temperature-Driven Structural and Morphological Evolution of Zinc Oxide Nano-Coalesced Microstructures and Its Defect-Related Photoluminescence Properties
Previous Article in Special Issue
Photocatalytic Properties of g-C3N4–TiO2 Heterojunctions under UV and Visible Light Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Shape-Dependent Single-Electron Levels for Au Nanoparticles

by
Georgios D. Barmparis
1,2,
Georgios Kopidakis
1,3 and
Ioannis N. Remediakis
1,3,*
1
Department of Materials Science and Technology, University of Crete, 71003 Heraklion, Greece
2
Crete Center for Quantum Complexity and Nanotechnology, Department of Physics, University of Crete, 71003 Heraklion, Greece
3
Institute of Electronic Structure and Laser, Foundation for Research and Technology-Hellas, 70013 Heraklion, Greece
*
Author to whom correspondence should be addressed.
Materials 2016, 9(4), 301; https://doi.org/10.3390/ma9040301
Submission received: 28 February 2016 / Revised: 8 April 2016 / Accepted: 12 April 2016 / Published: 21 April 2016
(This article belongs to the Special Issue Advancement of Photocatalytic Materials 2016)

Abstract

:
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.

1. 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 ( h k l ) 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.

2. Results

2.1. 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:
- 2 2 m 2 + V ( x , y , z ) ψ i ( x , y , z ) = E i ψ i ( x , y , z )
where E i and ψ i ( x , y , z ) are the i t h eigenvalue and eigenstate of the system, respectively; - 2 2 m 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:
V ( x , y , z ) = 0 inside the nanoparticle otherwise
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 ( h k l ) 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:
  • Cube;
  • Sphere;
  • Truncated octahedron, formed when Au has very weak interactions with its environment, such as Au nanoparticles bonded to C nanotubes [36];
  • Complex rounded shape, formed when Au is in equilibrium with small molecules with strong binding to Au, such as CO [37];
  • 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].

2.1.1. 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
ϕ n , m , l ( x , y , z ) = 8 a 3 sin ( n π x a ) sin ( m π y a ) sin ( l π z a )
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:
E n m l = 2 π 2 2 m e a 2 ( n 2 + m 2 + l 2 )
where, 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.

2.1.2. 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 n l , are given by:
E n l = 2 χ n l 2 2 m e R 2
where χ n l is the n t h root of the l t h order spherical Bessel functions j l ( x ) . The corresponding eigenfunctions are
ψ n l m ( r , θ , ϕ ) j l ( k n l r ) Y l m ( θ , ϕ )
where j l ( x ) is the spherical Bessel function of order l and Y l m ( θ , ϕ ) is a spherical harmonic. k n l is related to the energy through E n l = 2 k n l 2 2 m 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 ( 2 l + 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 ( 6 V ) 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:
ρ V ( E ) = 1 V i δ ( E - E i )
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 σ:
δ ( x ) = 1 π σ ( 1 + x 2 σ 2 )
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 E V 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.

2.1.3. 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, γ i n t , between Au surfaces and their environment. A simple relationship between the two has been found when there is little interaction between adsorbates [39]:
γ i n t = γ + θ E a d s A a t
In this equation, θ is the coverage of Au surfaces ( θ = 1 corresponds to full coverage) while E a d s is the adsorption energy of ligands (exothermic adsorption means negative E a d s ). A a t is the area per Au atom, that is, the total area divided by the number of surface Au atoms. According to Equation (9), γ i n t γ when the coverage, θ, or the adsorption energy, E a d s , 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 x y , d y z , and d z x orbitals.

2.1.4. 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.

2.2. 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.
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.

3. Materials and Methods

3.1. Atomistic Wulff Construction

The Au nanoparticles used in this study are generated by means of our atomistic-Wulff- construction method. The Wulff theorem states that the equilibrium shape of a crystal is a polyhedron. The distance of each ( h k l ) face of this polyhedron from the center is proportional to the surface tension of the ( h k l ) 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].

3.2. 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:
ψ i ( x , y , z ) = j N c i , j ϕ j ( x , y , z )
where c i , j is a scalar coefficient describing the weight of the contribution of the j t h 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:
ϕ n , m , l ( x , y , z ) = 8 a b c sin ( n π x a ) sin ( m π y b ) sin ( l π z c ) E n , m , l = 2 π 2 2 m e ( n 2 a 2 + m 2 b 2 + l 2 c 2 )
These are solutions of the Schrödinger equation for the potential energy V ( x , y , z ) defined as:
V ( x , y , z ) = 0 0 x a , 0 y b , and 0 z c otherwise
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 )
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 d V ,
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 r e g i o n , in each of them using the form:
V r e g i o n = 8 V 0 a b c x i x f y i y f 0 z f sin ( n π x a ) sin ( m π y b ) sin ( l π z c ) sin ( n π x a ) sin ( m π y b ) sin ( l π z c ) d z d y d x
Using the Euler’s formula, we rewrite each sine with exponential functions:
sin ( x ) = e i x - e - i x 2 i
and thus Equation (15) is written as:
V r e g i o n = 8 V 0 a b c x i x f y i y f 0 z f e i n π a x - e - i n π a x 2 i e i m π b y - e - i m π b y 2 i e i l π c z - e - i l π c z 2 i × e i n π a x - e - i n π a x 2 i e i m π b y - e - i m π b y 2 i e i l π c z - e - i l π c z 2 i d z d y d x
Expanding this integral we obtain that V r e g i o n is given by:
V r e g i o n = - 8 V 0 2 6 a b c x i x f y i y f 0 z f ( e + i π ( n + n ) a x e + i π ( m + m ) b y e + i π ( l + l ) c z - e + i π ( n + n ) a x e + i π ( m + m ) b y e + i π ( l - l ) c z - e + i π ( n + n ) a x e + i π ( m + m ) b y e - i π ( l - l ) c z + e + i π ( n + n ) a x e + i π ( m + m ) b y e - i π ( l + l ) c z - e + i π ( n + n ) a x e + i π ( m - m ) b y e + i π ( l + l ) c z + e + i π ( n + n ) a x e + i π ( m - m ) b y e + i π ( l - l ) c z + e + i π ( n + n ) a x e + i π ( m - m ) b y e - i π ( l - l ) c z - e + i π ( n + n ) a x e + i π ( m - m ) b y e - i π ( l + l ) c z - e + i π ( n + n ) a x e - i π ( m - m ) b y e + i π ( l + l ) c z + e + i π ( n + n ) a x e - i π ( m - m ) b y e + i π ( l - l ) c z + e + i π ( n + n ) a x e - i π ( m - m ) b y e - i π ( l - l ) c z - e + i π ( n + n ) a x e - i π ( m - m ) b y e - i π ( l + l ) c z + e + i π ( n + n ) a x e - i π ( m + m ) b y e + i π ( l + l ) c z - e + i π ( n + n ) a x e - i π ( m + m ) b y e + i π ( l - l ) c z - e + i π ( n + n ) a x e - i π ( m + m ) b y e - i π ( l - l ) c z + e + i π ( n + n ) a x e - i π ( m + m ) b y e - i π ( l + l ) c z - e + i π ( n - n ) a x e + i π ( m + m ) b y e + i π ( l + l ) c z + e + i π ( n - n ) a x e + i π ( m + m ) b y e + i π ( l - l ) c z + e + i π ( n - n ) a x e + i π ( m + m ) b y e - i π ( l - l ) c z - e + i π ( n - n ) a x e + i π ( m + m ) b y e - i π ( l + l ) c z + e + i π ( n - n ) a x e + i π ( m - m ) b y e + i π ( l + l ) c z - e + i π ( n - n ) a x e + i π ( m - m ) b y e + i π ( l - l ) c z - e + i π ( n - n ) a x e + i π ( m - m ) b y e - i π ( l - l ) c z + e + i π ( n - n ) a x e + i π ( m - m ) b y e - i π ( l + l ) c z + e + i π ( n - n ) a x e - i π ( m - m ) b y e + i π ( l + l ) c z - e + i π ( n - n ) a x e - i π ( m - m ) b y e + i π ( l - l ) c z - e + i π ( n - n ) a x e - i π ( m - m ) b y e - i π ( l - l ) c z + e + i π ( n - n ) a x e - i π ( m - m ) b y e - i π ( l + l ) c z - e + i π ( n - n ) a x e - i π ( m + m ) b y e + i π ( l + l ) c z + e + i π ( n - n ) a x e - i π ( m + m ) b y e + i π ( l - l ) c z + e + i π ( n - n ) a x e - i π ( m + m ) b y e - i π ( l - l ) c z - e + i π ( n - n ) a x e - i π ( m + m ) b y e - i π ( l + l ) c z - e - i π ( n - n ) a x e + i π ( m + m ) b y e + i π ( l + l ) c z + e - i π ( n - n ) a x e + i π ( m + m ) b y e + i π ( l - l ) c z + e - i π ( n - n ) a x e + i π ( m + m ) b y e - i π ( l - l ) c z - e - i π ( n - n ) a x e + i π ( m + m ) b y e - i π ( l + l ) c z + e - i π ( n - n ) a x e + i π ( m - m ) b y e + i π ( l + l ) c z - e - i π ( n - n ) a x e + i π ( m - m ) b y e + i π ( l - l ) c z - e - i π ( n - n ) a x e + i π ( m - m ) b y e - i π ( l - l ) c z + e - i π ( n - n ) a x e + i π ( m - m ) b y e - i π ( l + l ) c z + e - i π ( n - n ) a x e - i π ( m - m ) b y e + i π ( l + l ) c z - e - i π ( n - n ) a x e - i π ( m - m ) b y e + i π ( l - l ) c z - e - i π ( n - n ) a x e - i π ( m - m ) b y e - i π ( l - l ) c z + e - i π ( n - n ) a x e - i π ( m - m ) b y e - i π ( l + l ) c z - e - i π ( n - n ) a x e - i π ( m + m ) b y e + i π ( l + l ) c z + e - i π ( n - n ) a x e - i π ( m + m ) b y e + i π ( l - l ) c z + e - i π ( n - n ) a x e - i π ( m + m ) b y e - i π ( l - l ) c z - e - i π ( n - n ) a x e - i π ( m + m ) b y e - i π ( l + l ) c z + e - i π ( n + n ) a x e + i π ( m + m ) b y e + i π ( l + l ) c z - e - i π ( n + n ) a x e + i π ( m + m ) b y e + i π ( l - l ) c z - e - i π ( n + n ) a x e + i π ( m + m ) b y e - i π ( l - l ) c z + e - i π ( n + n ) a x e + i π ( m + m ) b y e - i π ( l + l ) c z - e - i π ( n + n ) a x e + i π ( m - m ) b y e + i π ( l + l ) c z + e - i π ( n + n ) a x e + i π ( m - m ) b y e + i π ( l - l ) c z + e - i π ( n + n ) a x e + i π ( m - m ) b y e - i π ( l - l ) c z - e - i π ( n + n ) a x e + i π ( m - m ) b y e - i π ( l + l ) c z - e - i π ( n + n ) a x e - i π ( m - m ) b y e + i π ( l + l ) c z + e - i π ( n + n ) a x e - i π ( m - m ) b y e + i π ( l - l ) c z + e - i π ( n + n ) a x e - i π ( m - m ) b y e - i π ( l - l ) c z - e - i π ( n + n ) a x e - i π ( m - m ) b y e - i π ( l + l ) c z + e - i π ( n + n ) a x e - i π ( m + m ) b y e + i π ( l + l ) c z - e - i π ( n + n ) a x e - i π ( m + m ) b y e + i π ( l - l ) c z - e - i π ( n + n ) a x e - i π ( m + m ) b y e - i π ( l - l ) c z + e - i π ( n + n ) a x e - i π ( m + m ) b y e - i π ( l + l ) c z ) d z d y d x
The above integral can be reduced using the fact that its 64 terms are all of the form:
s ( t ) e i π N ( t ) a x e i π M ( t ) b y e i π L ( t ) c z
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 r e g i o n 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:
V r e g i o n = - 8 V 0 2 6 a 3 t = 1 64 s ( t ) x i x f y i y f 0 z f e i π N ( t ) a x e + i π M ( t ) a y e + i π L ( t ) a z d z d y d x
With limits defined by the borders of each region given as: x i < x < x f , the lines, y i = a i x + b i < y < y f = a f x + b f , and the surfaces z i = 0 < z < z f = d - h x - k y l , where d = h 2 + k 2 + l 2 d h k l . This triple integral of each one of these terms can be calculated analytically and is given by:
x i x f y i y f 0 z f e i π N ( t ) a x e + i π M ( t ) a y e + i π L ( t ) a z d z d y d x = - 2 a 3 l 2 π 3 L ( t ) [ M ( t ) l - L ( t ) k ] [ N ( t ) l - L ( t ) h + [ M ( t ) l - L ( t ) k ] a f ] { sin [ π d a l L ( t ) + π b f a l [ M ( t ) l - L ( t ) k ] + π x f a l [ N ( t ) l - L ( t ) h + [ M ( t ) l - L ( t ) k ] a f ] ] - sin [ π d a l L ( t ) + π b f a l [ M ( t ) l - L ( t ) k ] + π x i a l [ N ( t ) l - L ( t ) h + [ M ( t ) l - L ( t ) k ] a f ] ] } + 2 a 3 l 2 π 3 L ( t ) [ M ( t ) l - L ( t ) k ] [ N ( t ) l - L ( t ) h + [ M ( t ) l - L ( t ) k ] a i ] { sin [ π d a l L ( t ) + π b i a l [ M ( t ) l - L ( t ) k ] + π x f a l [ N ( t ) l - L ( t ) h + [ M ( t ) l - L ( t ) k ] a i ] ] - sin [ π d a l L ( t ) + π b i a l [ M ( t ) l - L ( t ) k ] + π x i a l [ N ( t ) l - L ( t ) h + [ M ( t ) l - L ( t ) k ] a i ] ] } + 2 a 3 π 3 L ( t ) M ( t ) [ N ( t ) + a f M ( t ) ] { sin [ π b f a M ( t ) + π x f a [ N ( t ) + a f M ( t ) ] ] - s i n [ π b f a M ( t ) + π x i a [ N ( t ) + a f M ( t ) ] ] } + 2 a 3 π 3 L ( t ) M ( t ) [ N ( t ) + a i M ( t ) ] { sin [ π b i a M ( t ) + π x f a [ N ( t ) + a i M ( t ) ] ] - s i n [ π b i a M ( t ) + π x i a [ N ( t ) + a i M ( t ) ] ] }
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 r e g i o n of each one of the regions.
Finally, the expectation value for the potential energy is given by the sum of V r e g i o n over all the regions as:
( ϕ n , m , l , V ( x , y , z ) ϕ n , m , l ) = r e g i o n s V r e g i o n
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.

4. 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.

Acknowledgments

Georgios D. Barmparis acknowledges support from EU/FP7-REGPOT-2012-2013-1 under Grant Agreement No. 316165. This work was supported by the Research Committee of the University of Crete K.A. 3474 and by the by Greek GSRT project ERC02-EXEL Grant No. 6260.

Author Contributions

Georgios D. Barmparis made the atomistic Wulff constructions, wrote and run the code for the energy eigenvalues in the nanoparticles. Georgios Kopidakis helped in analyzing the data and produced the densities of states. Ioannis N. Remediakis conceived and designed the simulations and helped in analyzing the data.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
DFT
Density-functional Theory
DOS
Density of states

References

  1. Sun, Y.; Xia, Y. Shape-Controlled Synthesis of Gold and Silver Nanoparticles. Science 2002, 298, 2176–2179. [Google Scholar] [CrossRef] [PubMed]
  2. Grzelczak, M.; Pérez-Juste, J.; Mulvaney, P.; Liz-Marzán, L. Shape control in gold nanoparticle synthesis. Chem. Soc. Rev. 2008, 37, 1783–1791. [Google Scholar] [CrossRef] [PubMed]
  3. Mostafa, S.; Behafarid, F.; Croy, J.R.; Ono, L.K.; Li, L.; Yang, J.C.; Frenkel, A.I.; Cuenya, B.R. Shape-Dependent Catalytic Properties of Pt Nanoparticles. J. Am. Chem. Soc. 2010, 132, 15714–15719. [Google Scholar] [CrossRef] [PubMed]
  4. Kovács, G.; Fodor, S.; Vulpoi, A.; Schrantz, K.; Dombi, A.; Hernádi, K.; Danciu, V.; Pap, Z.; Baia, L. Polyhedral Pt vs. spherical Pt nanoparticles on commercial titanias: Is shape tailoring a guarantee of achieving high activity? J. Catal. 2015, 325, 156–167. [Google Scholar] [CrossRef]
  5. Vajda, K.; Saszet, K.; Kedves, E.; Kása, Z.; Danciu, V.; Baia, L.; Magyari, K.; Hernádi, K.; Kovács, G.; Pap, Z. Shape-controlled agglomeration of TiO2 nanoparticles. New insights on polycrystallinity vs. single crystals in photocatalysis. Ceram. Int. 2016, 42, 3077–3087. [Google Scholar] [CrossRef]
  6. Remediakis, I.N.; Lopez, N.; Nørskov, J.K. CO Oxidation on Rutile-Supported Au Nanoparticles. Angew. Chem. Int. Ed. 2005, 44, 1824–1826. [Google Scholar] [CrossRef] [PubMed]
  7. Barnard, A.S.; Lin, X.M.; Curtiss, L.A. Equilibrium Morphology of Face-Centered Cubic Gold Nanoparticles >3 nm and the Shape Changes Induced by Temperature. J. Phys. Chem. B 2005, 109, 24465–24472. [Google Scholar] [CrossRef] [PubMed]
  8. Fasi, A.; Palinko, I.; Hernadi, K.; Kiricsi, I. Formation of Au nanorods and nanoforks over MgO support. React. Kinet. Catal. Lett. 2006, 87, 263–268. [Google Scholar] [CrossRef]
  9. Ahmadi, T.S.; Wang, Z.L.; Green, T.C.; Henglein, A.; Sayed, M.A.E. Shape-Controlled Synthesis of Colloidal Platinum Nanoparticles. Science 2007, 272, 1924–1926. [Google Scholar] [CrossRef]
  10. Ruffino, F.; Bongiorno, C.; Giannazzo, F.; Roccaforte, F.; Raineri, V.; Grimaldi, M. Effect of surrounding environment on atomic structure and equilibrium shape of growing nanocrystals: Gold in/on SiO2. Nanoscale Res. Lett. 2007, 2, 240–247. [Google Scholar] [CrossRef] [PubMed]
  11. Kittel, C. Introduction to Solid State Physics, 8th ed.; John Wiley & Sons, Inc.: New York, NY, USA, 2005; Chapter 18. [Google Scholar]
  12. Van Aert, S.; Batenburg, K.J.; Rossell, M.D.; Erni, R.; Van Tendeloo, G. Three-dimensional atomic imaging of crystalline nanoparticles. Nature 2011, 470, 374–377. [Google Scholar] [CrossRef] [PubMed]
  13. Sivaramakrishnan, S.; Wen, J.; Scarpelli, M.E.; Pierce, B.J.; Zuo, J.M. Equilibrium shapes and triple line energy of epitaxial gold nanocrystals supported on TiO2(110). Phys. Rev. B 2010, 82. [Google Scholar] [CrossRef]
  14. De Heer, W.A. The physics of simple metal clusters: Experimental aspects and simple models. Rev. Mod. Phys. 1993, 65, 611–676. [Google Scholar] [CrossRef]
  15. Walter, M.; Akola, J.; Lopez-Acevedo, O.; Jadzinsky, P.D.; Calero, G.; Ackerson, C.J.; Whetten, R.L.; Grönbeck, H.; Häkkinen, H. A unified view of ligand-protected gold clusters as superatom complexes. Proc. Natl. Acad. Sci. USA 2008, 105, 9157–9162. [Google Scholar] [CrossRef] [PubMed]
  16. Stiehler, C.; Calaza, F.; Schneider, W.D.; Nilius, N.; Freund, H.J. Molecular Adsorption Changes the Quantum Structure of Oxide-Supported Gold Nanoparticles: Chemisorption versus Physisorption. Phys. Rev. Lett. 2015, 115. [Google Scholar] [CrossRef] [PubMed]
  17. Pelton, M.; Aizpurua, J.; Bryant, G. Metal-nanoparticle plasmonics. Laser Photonics Rev. 2008, 2, 136–159. [Google Scholar] [CrossRef]
  18. West, P.; Ishii, S.; Naik, G.; Emani, N.; Shalaev, V.; Boltasseva, A. Searching for better plasmonic materials. Laser Photonics Rev. 2010, 4, 795–808. [Google Scholar] [CrossRef]
  19. Mie, G. Beiträge zur Optik trüber Medien, speziell kolloidaler Metallösungen. Ann. Phys. 1908, 330, 377–445. [Google Scholar] [CrossRef]
  20. Maier, S.A. Plasmonics: Fundamentals and Applications, 1st ed.; Springer-Verlag US: New York, NY, USA, 2007. [Google Scholar]
  21. Klimov, V.I.; Mikhailovski, A.A.; Xu, S.; Malko, A.; Hollingsworth, J.A.; Leatherdale, C.A.; Eisler, H.J.; Bawendi, M.G. Optical Gain and Stimulated Emission in Nanocrystal Quantum Dots. Science 2000, 290, 314–317. [Google Scholar] [CrossRef] [PubMed]
  22. Lal, S.; Link, S.; Halas, N.J. Nano-optics from sensing to waveguiding. Nat. Photonics 2007, 1, 641–648. [Google Scholar] [CrossRef]
  23. Kim, S.K.; Day, R.W.; Cahoon, J.F.; Kempa, T.J.; Song, K.D.; Park, H.G.; Lieber, C.M. Tuning Light Absorption in Core/Shell Silicon Nanowire Photovoltaic Devices through Morphological Design. Nano Lett. 2012, 12, 4971–4976. [Google Scholar] [CrossRef] [PubMed]
  24. Beqa, L.; Singh, A.K.; Khan, S.A.; Senapati, D.; Arumugam, S.R.; Ray, P.C. Gold Nanoparticle-Based Simple Colorimetric and Ultrasensitive Dynamic Light Scattering Assay for the Selective Detection of Pb(II) from Paints, Plastics, and Water Samples. ACS Appl. Mater. Interfaces 2011, 3, 668–673. [Google Scholar] [CrossRef] [PubMed]
  25. Hourahine, B.; Papoff, F. The geometrical nature of optical resonances: From a sphere to fused dimer nanoparticles. Meas. Sci. Technol. 2012, 23, 084002. [Google Scholar] [CrossRef]
  26. Hourahine, B.; Papoff, F. Optical control of scattering, absorption and lineshape in nanoparticles. Opt. Express 2013, 21, 20322–20333. [Google Scholar] [CrossRef] [PubMed]
  27. Von Delft, J.; Ralph, D. Spectroscopy of discrete energy levels in ultrasmall metallic grains. Phys. Rep. 2001, 345, 61–173. [Google Scholar] [CrossRef]
  28. Bolotin, K.I.; Kuemmeth, F.; Pasupathy, A.N.; Ralph, D.C. Metal-nanoparticle single-electron transistors fabricated using electromigration. Appl. Phys. Lett. 2004, 84, 3154–3156. [Google Scholar] [CrossRef]
  29. Kuemmeth, F.; Bolotin, K.I.; Shi, S.F.; Ralph, D.C. Measurement of Discrete Energy-Level Spectra in Individual Chemically Synthesized Gold Nanoparticles. Nano Lett. 2008, 8, 4506–4512. [Google Scholar] [CrossRef] [PubMed]
  30. Lin, X.; Nilius, N.; Freund, H.J.; Walter, M.; Frondelius, P.; Honkala, K.; Häkkinen, H. Quantum Well States in Two-Dimensional Gold Clusters on MgO Thin Films. Phys. Rev. Lett. 2009, 102. [Google Scholar] [CrossRef] [PubMed]
  31. Kac, M. Can One Hear the Shape of a Drum? Am. Math. Mon. 1966, 73, 1–23. [Google Scholar] [CrossRef]
  32. Hearing the Shape of a Drum. Available online: https://en.wikipedia.org/wiki/Hearing_the_shape_of_a_drum (accessed on 30 January 2016).
  33. Tolea, F.; Tolea, M. Hearing shapes of few electrons quantum drums: A configuration-interaction study. Physica B 2015, 458, 85–91. [Google Scholar] [CrossRef]
  34. Lopez, N.; Nørskov, J.; Janssens, T.; Carlsson, A.; Puig-Molina, A.; Clausen, B.; Grunwaldt, J.D. The adhesion and shape of nanosized Au particles in a Au/TiO2 catalyst. J. Catal. 2004, 225, 86–94. [Google Scholar] [CrossRef]
  35. Remediakis, I.N.; Lopez, N.; Nørskov, J.K. CO oxidation on gold nanoparticles: Theoretical studies. Appl. Catal. A Gen. 2005, 291, 13–21. [Google Scholar] [CrossRef]
  36. Quintana, M.; Ke, X.; Tendeloo, G.V.; Meneghetti, M.; Bittencourt, C.; Prato, M. Light-Induced Selective Deposition of Au Nanoparticles on Single-Wall Carbon Nanotubes. ACS Nano 2010, 4, 6105–6113. [Google Scholar] [CrossRef] [PubMed]
  37. Ueda, K.; Kawasaki, T.; Hasegawa, H.; Tanji, T.; Ichihashi, M. First observation of dynamic shape changes of a gold nanoparticle catalyst under reaction gas environment by transmission electron microscopy. Surf. Interface Anal. 2008, 40, 1725–1727. [Google Scholar] [CrossRef]
  38. Häkkinen, H.; Walter, M.; Grönbeck, H. Divide and Protect: Capping Gold Nanoclusters with Molecular Gold-Thiolate Rings. J. Phys. Chem. B 2006, 110, 9927–9931. [Google Scholar] [CrossRef] [PubMed]
  39. Barmparis, G.D.; Remediakis, I.N. Dependence on CO adsorption of the shapes of multifaceted gold nanoparticles: A density functional theory. Phys. Rev. B 2012, 86. [Google Scholar] [CrossRef]
  40. Bittencourt, C.; Felten, A.; Douhard, B.; Colomer, J.F.; Tendeloo, G.V.; Drube, W.; Ghijsen, J.; Pireaux, J.J. Metallic nanoparticles on plasma treated carbon nanotubes: Nano2hybrids. Surf. Sci. 2007, 601, 2800–2804. [Google Scholar] [CrossRef]
  41. Espinosa, E.; Ionescu, R.; Bittencourt, C.; Felten, A.; Erni, R.; Tendeloo, G.V.; Pireaux, J.J.; Llobet, E. Metal-decorated multi-wall carbon nanotubes for low temperature gas sensing. Thin Solid Films 2007, 515, 8322–8327. [Google Scholar] [CrossRef]
  42. Barmparis, G.D.; Honkala, K.; Remediakis, I.N. Thiolate adsorption on Au(hkl) and equilibrium shape of large thiolate-covered gold nanoparticles. J. Chem. Phys. 2013, 138. [Google Scholar] [CrossRef] [PubMed]
  43. Marago, O.M.; Jones, P.H.; Gucciardi, P.G.; Volpe, G.; Ferrari, A.C. Optical trapping and manipulation of nanostructures. Nat. Nano 2013, 8, 807–819. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Plummer, E. Electronic states and phases of KxC60 from photoemission and X-ray absorption spectroscopy. Nature 1991, 352, 603–605. [Google Scholar]
  45. Liu, W.; Hu, E.; Jiang, H.; Xiang, Y.; Weng, Z.; Li, M.; Fan, Q.; Yu, X.; Altman, E.I.; Wang, H. A highly active and stable hydrogen evolution catalyst based on pyrite-structured cobalt phosphosulfide. Nat. Commun. 2016, 7. [Google Scholar] [CrossRef] [PubMed]
  46. Wulff, G. Zur Frage der Geschwindigkeit des Wachstums und der Auflösung der Krystallflächen. Z. Kristallogr. 1901, 34, 449–530. [Google Scholar]
  47. Herring, C. Some Theorems on the Free Energies of Crystal Surfaces. Phys. Rev. 1951, 82. [Google Scholar] [CrossRef]
  48. Rosakis, P. Continuum Surface Energy from a Lattice Model. Netw. Heterog. Media 2014, 9, 453–476. [Google Scholar] [CrossRef]
  49. Molina, L.M.; Hammer, B. Active Role of Oxide Support during CO Oxidation at Au/MgO. Phys. Rev. Lett. 2003, 90. [Google Scholar] [CrossRef] [PubMed]
  50. Barnard, A.; Zapol, P. A model for the phase stability of arbitrary nanoparticles as a function of size and shape. J. Chem. Phys. 2004, 121, 4276–4283. [Google Scholar] [CrossRef] [PubMed]
  51. Barnard, A.S.; Curtiss, L.A. Prediction of TiO2 Nanoparticle Phase and Shape Transitions Controlled by Surface Chemistry. Nano Lett. 2005, 5, 1261–1266. [Google Scholar] [CrossRef] [PubMed]
  52. Hadjisavvas, G.; Remediakis, I.N.; Kelires, P.C. Insights into the Shape and Faceting of Embedded Si/α-SiO2 Nanocrystals. Phys. Rev. B 2006, 74. [Google Scholar] [CrossRef]
  53. Kopidakis, G.; Remediakis, I.N.; Fyta, M.G.; Kelires, P.C. Atomic and electronic structure of crystalline-amorphous carbon interfaces. Diam. Relat. Mater. 2007, 16, 1875–1881. [Google Scholar] [CrossRef]
  54. Mittendorfer, F.; Seriani, N.; Dubay, O.; Kresse, G. Morphology of mesoscopic Rh and Pd nanoparticles under oxidizing conditions. Phys. Rev. B 2007, 76. [Google Scholar] [CrossRef]
  55. Soon, A.; Wong, L.; Delley, B.; Stampfl, C. Morphology of copper nanoparticles in a nitrogen atmosphere: A first-principles investigation. Phys. Rev. B 2008, 77. [Google Scholar] [CrossRef]
  56. Shi, H.; Stampfl, C. Shape and surface structure of gold nanoparticles under oxidizing conditions. Phys. Rev. B 2008, 77. [Google Scholar] [CrossRef]
  57. Cortes-Huerto, R.; Goniakowski, J.; Noguera, C. An efficient many-body potential for the interaction of transition and noble metal nano-objects with an environment. J. Chem. Phys. 2013, 138. [Google Scholar] [CrossRef] [PubMed]
  58. Kim, K.C.; Dai, B.; Johnson, J.K.; Sholl, D.S. Assessing nanoparticle size effects on metal hydride thermodynamics using the Wulff construction. Nanotechnology 2009, 20. [Google Scholar] [CrossRef] [PubMed]
  59. Pham, T.H.; Duan, X.; Qian, G.; Zhou, X.; Chen, D. CO Activation Pathways of Fischer Tropsch Synthesis on χ-Fe5C2 (510): Direct versus Hydrogen-Assisted CO Dissociation. J. Phys. Chem. C 2014, 118, 10170–10176. [Google Scholar] [CrossRef]
  60. Lodziana, Z.; Stoica, G.; Perez-Ramirez, J. Reevaluation of the Structure and Fundamental Physical Properties of Dawsonites by DFT Studies. Inorg. Chem. 2011, 50, 2590–2598. [Google Scholar] [CrossRef] [PubMed]
  61. Honkala, K.; Lodziana, Z.; Remediakis, I.N.; Lopez, N. Expanding and Reducing Complexity in Materials Science Models with Relevance in Catalysis and Energy. Top. Catal. 2014, 57, 14–24. [Google Scholar] [CrossRef]
  62. Honkala, K.; Hellman, A.; Remediakis, I.N.; Logadottir, A.; Carlsson, A.; Dahl, S.; Christensen, C.; Nørskov, J.K. Ammonia synthesis from first-principles calculations. Science 2005, 307, 555–558. [Google Scholar] [CrossRef] [PubMed]
  63. Hellman, A.; Honkala, K.; Remediakis, I.N.; Logadottir, A.; Carlsson, A.; Dahl, S.; Christensen, C.H.; Nørskov, J.K. Insights into ammonia synthesis from first-principles. Surf. Sci. 2006, 600, 4264–4268. [Google Scholar] [CrossRef]
  64. Hellman, A.; Honkala, K.; Remediakis, I.N.; Logadottir, A.; Carlsson, A.; Dahl, S.; Christensen, C.H.; Nørskov, J.K. Ammonia synthesis and decomposition on a Ru-based catalyst modeled by first-principles. Surf. Sci. 2009, 603, 1731–1739. [Google Scholar] [CrossRef]
  65. Vile, G.; Baudouin, D.; Remediakis, I.N.; Coperet, C.; Lopez, N.; Perez-Ramirez, J. Silver Nanoparticles for Olefin Production: New Insights into the Mechanistic Description of Propyne Hydrogenation. ChemCatChem 2013, 5, 3750–3759. [Google Scholar] [CrossRef]
  66. Gómez-Graña, S.; Goris, B.; Altantzis, T.; Fernández-López, C.; Carbó-Argibay, E.; Guerrero-Martínez, A.; Almora-Barrios, N.; López, N.; Pastoriza-Santos, I.; Pérez-Juste, J.; et al. Au@Ag nanoparticles: Halides stabilize 100 facets. J. Phys. Chem. Lett. 2013, 4, 2209–2216. [Google Scholar] [CrossRef]
  67. Almora-Barrios, N.; Novell-Leruth, G.; Whiting, P.; Liz-Marzán, L.; López, N. Theoretical description of the role of halides, silver, and surfactants on the structure of gold nanorods. Nano Lett. 2014, 14, 871–875. [Google Scholar] [CrossRef] [PubMed]
  68. Barmparis, G.D.; Lodziana, Z.; Lopez, N.; Remediakis, I.N. Nanoparticle shapes by using Wulff constructions and first-principles calculations. Beilstein J. Nanotechnol. 2015, 6, 361–368. [Google Scholar] [CrossRef] [PubMed]
  69. Barmparis, G.D.; Maniadaki, A.E.; Kopidakis, G.; Remediakis, I.N. Wulff construction and molecular dynamics simulations for Au nanoparticles. J. Chem. Eng. Chem. Res. 2016. accepted. [Google Scholar]
Figure 1. The probability function, | ψ | 2 , of the first ten eigenstates of an electron inside a sphere. Yellow represents positive and cyan represents negative values of the wavefunction.
Figure 1. The probability function, | ψ | 2 , of the first ten eigenstates of an electron inside a sphere. Yellow represents positive and cyan represents negative values of the wavefunction.
Materials 09 00301 g001
Figure 2. Comparison of the ρ V ( E ) of an electron in a spherical nanoparticle of volume V = 481.42 nm 3 using the calculated (red line) and the analytical eigenenergies (black dotted line). Density of States (DOS) has been calculated using Lorentzian broadening of 1 meV in (a); and 5 meV in (b).
Figure 2. Comparison of the ρ V ( E ) of an electron in a spherical nanoparticle of volume V = 481.42 nm 3 using the calculated (red line) and the analytical eigenenergies (black dotted line). Density of States (DOS) has been calculated using Lorentzian broadening of 1 meV in (a); and 5 meV in (b).
Materials 09 00301 g002
Figure 3. Ball-and-stick models of typical equilibrium shapes for Au nanoparticles. Step and kink atoms are shown in darker color. (Left): clean nanoparticle; (center): particle in CO-rich gas at typical CO oxidation conditions; (right): particle in equilibrium with thiol-containing molecules at typical polymer concentrations.
Figure 3. Ball-and-stick models of typical equilibrium shapes for Au nanoparticles. Step and kink atoms are shown in darker color. (Left): clean nanoparticle; (center): particle in CO-rich gas at typical CO oxidation conditions; (right): particle in equilibrium with thiol-containing molecules at typical polymer concentrations.
Materials 09 00301 g003
Figure 4. The probability function, | ψ | 2 , of the lowest-energy eigenstates of an electron inside a clean Au nanoparticle at equilibrium. Yellow represents positive and cyan negative values of the wavefunction.
Figure 4. The probability function, | ψ | 2 , of the lowest-energy eigenstates of an electron inside a clean Au nanoparticle at equilibrium. Yellow represents positive and cyan negative values of the wavefunction.
Materials 09 00301 g004
Figure 5. The probability function, | ψ | 2 , of the first ten eigenstates of an electron inside a thiolate-protected Au nanoparticle. Yellow represents positive and cyan negative values of the wavefunction.
Figure 5. The probability function, | ψ | 2 , of the first ten eigenstates of an electron inside a thiolate-protected Au nanoparticle. Yellow represents positive and cyan negative values of the wavefunction.
Materials 09 00301 g005
Figure 6. The DOS per volume of all shapes considered in this study: spherical (dotted black line), cubic (red line), clean (dashed green line), CO-covered (long-dashed blue line), and thiolate-covered (dot-dashed purple line). DOS has been calculated using Lorentzian broadening of either 1 meV (a) or 5 meV (b). The energies have been multiplied by the volume of the nanoparticle to the exponent 2/3.
Figure 6. The DOS per volume of all shapes considered in this study: spherical (dotted black line), cubic (red line), clean (dashed green line), CO-covered (long-dashed blue line), and thiolate-covered (dot-dashed purple line). DOS has been calculated using Lorentzian broadening of either 1 meV (a) or 5 meV (b). The energies have been multiplied by the volume of the nanoparticle to the exponent 2/3.
Materials 09 00301 g006
Figure 7. (Left): The polyhedron that forms the boundary of the nanoparticle enclosed in the simulated box with dimensions ( a , b , c ); (Right): The octant used to calculate the potential integral. Grey color represents the inner part of the nanoparticle.
Figure 7. (Left): The polyhedron that forms the boundary of the nanoparticle enclosed in the simulated box with dimensions ( a , b , c ); (Right): The octant used to calculate the potential integral. Grey color represents the inner part of the nanoparticle.
Materials 09 00301 g007
Table 1. Comparison of calculated, E i C a l c , and analytically obtained values, E i , of the first thirty eigenenergies of a spherical nanoparticle with volume V = 481.42 nm 3 .
Table 1. Comparison of calculated, E i C a l c , and analytically obtained values, E i , of the first thirty eigenenergies of a spherical nanoparticle with volume V = 481.42 nm 3 .
i E i c a l c (meV) E i (meV)i E i c a l c (meV) E i (meV)i E i c a l c (meV) E i (meV)
116.57715.9071181.60578.70421111.900107.914
233.69032.5431281.69378.70422112.031107.914
333.95332.5431381.79178.70423112.122107.914
433.95332.5431481.79178.70424112.122107.914
555.53453.5381582.01778.70425112.259107.914
655.53453.5381682.34378.70426112.404107.914
755.65253.5381782.34378.70427112.404107.914
855.84253.5381899.60596.18928112.851107.914
956.07053.53819100.37096.18929113.054107.914
1066.30763.63020100.37096.18930138.328133.323
Table 2. Energy eigenvalues for five different nanoparticle shapes. The scaled energy, Δ E i = ( E i - E 0 ) V 2 3 is shown, where E i is the i-th energy eigenvalue and V is te volume of each nanoparticle.
Table 2. Energy eigenvalues for five different nanoparticle shapes. The scaled energy, Δ E i = ( E i - E 0 ) V 2 3 is shown, where E i is the i-th energy eigenvalue and V is te volume of each nanoparticle.
i Δ E i (eV nm 2 )
SphereCOSCH 3 CleanCube
00.0000.0000.0000.0000.000
11.0221.0681.0891.1051.128
21.0221.0681.0891.1051.128
31.0221.0691.0891.1051.128
42.3122.4022.3882.1832.256
52.3122.4022.3892.1832.256
62.3122.4022.5172.7212.256
72.3122.4462.5172.7213.008
82.3122.4462.5172.7213.008
92.9313.0713.1203.0303.008
103.8573.9814.0413.7743.384

Share and Cite

MDPI and ACS Style

Barmparis, G.D.; Kopidakis, G.; Remediakis, I.N. Shape-Dependent Single-Electron Levels for Au Nanoparticles. Materials 2016, 9, 301. https://doi.org/10.3390/ma9040301

AMA Style

Barmparis GD, Kopidakis G, Remediakis IN. Shape-Dependent Single-Electron Levels for Au Nanoparticles. Materials. 2016; 9(4):301. https://doi.org/10.3390/ma9040301

Chicago/Turabian Style

Barmparis, Georgios D., Georgios Kopidakis, and Ioannis N. Remediakis. 2016. "Shape-Dependent Single-Electron Levels for Au Nanoparticles" Materials 9, no. 4: 301. https://doi.org/10.3390/ma9040301

APA Style

Barmparis, G. D., Kopidakis, G., & Remediakis, I. N. (2016). Shape-Dependent Single-Electron Levels for Au Nanoparticles. Materials, 9(4), 301. https://doi.org/10.3390/ma9040301

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop