Next Article in Journal
Technique for High-Quality Protein Crystal Growth by Control of Subgrain Formation under an External Electric Field
Next Article in Special Issue
Elastodynamic Analysis of a Hollow Cylinder with Decagonal Quasicrystal Properties: Meshless Implementation of Local Integral Equations
Previous Article in Journal
Strategies to Approach Stabilized Plasticity in Metals with Diminutive Volume: A Brief Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Exact Method to Determine the Photonic Resonances of Quasicrystals Based on Discrete Fourier Harmonics of Higher-Dimensional Atomic Surfaces

by
Farhad A. Namin
1 and
Douglas H. Werner
2,*
1
Department of Electrical Engineering, Amirkabir University of Technology (Tehran Polytechnic), Tehran 15875, Iran
2
Department of Electrical Engineering, The Pennsylvania State University, University Park, PA 16802, USA
*
Author to whom correspondence should be addressed.
Crystals 2016, 6(8), 93; https://doi.org/10.3390/cryst6080093
Submission received: 18 May 2016 / Revised: 31 July 2016 / Accepted: 4 August 2016 / Published: 10 August 2016
(This article belongs to the Special Issue Structure and Properties of Quasicrystals 2016)

Abstract

:
A rigorous method for obtaining the diffraction patterns of quasicrystals is presented. Diffraction patterns are an essential analytical tool in the study of quasicrystals, since they can be used to determine their photonic resonances. Previous methods for approximating the diffraction patterns of quasicrystals have relied on evaluating the Fourier transform of finite-sized super-lattices. Our approach, on the other hand, is exact in the sense that it is based on a technique that embeds quasicrystals into higher dimensional periodic hyper-lattices, thereby completely capturing the properties of the infinite structure. The periodicity of the unit cell in the higher dimensional space can be exploited to obtain the Fourier series expansion in closed-form of the corresponding atomic surfaces. The utility of the method is demonstrated by applying it to one-dimensional Fibonacci and two-dimensional Penrose quasicrystals. The results are verified by comparing them to those obtained by using the conventional super-lattice method. It is shown that the conventional super-cell approach can lead to inaccurate results due to the continuous nature of the Fourier transform, since quasicrystals have a discrete spectrum, whereas the approach introduced in this paper generates discrete Fourier harmonics. Furthermore, the conventional approach requires very large super-cells and high-resolution sampling of the reciprocal space in order to produce accurate results leading to a very large computational burden, whereas the proposed method generates accurate results with a relatively small number of terms. Finally, we propose how this approach can be generalized from the vertex model, which assumes identical particles at all vertices, to a more realistic case where the quasicrystal is composed of different atoms.

Graphical Abstract

1. Introduction

A crystal is defined as a solid having an essentially discrete diffraction pattern [1]. In 1984, Dr. Daniel Shechtman reported that the electron diffraction pattern of a rapidly solidified alloy of aluminum (Al) and manganese (Mn) displayed a discrete spectrum and five-fold rotational symmetry [2]. The discrete diffraction peaks satisfied the definition of a crystal; however, mathematically, five-fold rotational symmetry could not be produced by a strictly periodic lattice [1]. His results were highly controversial and met with great resistance in the academic community, since it was a commonly-held belief that all crystals possessed translational symmetry. Over the next decade, several other groups reported similar observations, and this new class of crystals that lacked periodicity came to be known as quasicrystals (QCs). In recent years, there has been considerable interest in utilizing quasicrystalline morphologies in a wide range of applications, such as surface-enhanced Raman scattering (SERS) substrates [3], electronic band gap materials [4,5,6], ultra-wideband antenna arrays [7,8,9] and broadband plasmonic enhancement [10,11,12].
The diffraction patterns of QC lattices consist of discrete, Bragg-like peaks. The presence of discrete Bragg peaks suggests that, similar to periodic lattices, QCs possess photonic resonances. Additionally, due to their unique properties, QC lattices can possess multiple photonic resonances in close proximity. These properties of QC lattices make them excellent candidates for a wide variety of broadband optical applications, such as SERS substrates or broadband absorbers. Recently, plasmon-assisted coatings have been proposed to enhance the efficiency of solar cells [13]. Using QC plasmonic lattices in solar cells can enhance the absorption spectra. Moreover, due to the higher orders of symmetry, the response will be more isotropic and polarization insensitive [14].
A significant obstacle in studying QCs is the lack of analytical tools to accurately and efficiently model them. An indispensable tool for analyzing the resonant properties of QCs is their diffraction pattern. The diffraction patterns of QCs have been successfully used to predict their photonic resonances [12,15,16]. It has been shown that discrete peaks in the diffraction pattern can be associated with vectors in the reciprocal space, referred to as “reciprocal vectors” (RVs). These RVs can be directly associated with photonic resonances of QC lattices [16,17]. Thus, obtaining accurate diffraction patterns for QCs is essential, since they can be used to predict their resonances.
There is a close relationship between QCs and certain aperiodic tilings of the one-dimensional (1D) and two-dimensional (2D) plane [1]. Aperiodic tilings are generated from a set of prototiles and placing the prototiles next to each other according to a given set of matching rules. Alternatively, aperiodic tilings can be generated using a process known as Stone inflation. Stone inflation, also sometimes referred to as inflate-and-substitute, specifies how each prototile is first expanded by a given factor and then subdivided into smaller prototiles [18]. Other examples include the Danzer [1] aperiodic tiling, which possesses seven-fold rotational symmetry and the Ammann–Beenker aperiodic tiling with eight-fold rotational symmetry [18]. Converting an aperiodic tiling to a point lattice simply requires placing points at the vertices of the tiling.
The Fourier diffraction pattern of a QC is traditionally obtained by applying a Fourier transform (FT) to a finite-sized point lattice based on the QC morphology and subsequently sampling the reciprocal space to capture all of the diffraction spots. This traditional approach can often lead to inaccurate results due to the fact that the reciprocal space of QCs is densely filled [19]. To obtain an accurate diffraction pattern, two requirements must be taken into account:
(1)
The approach attempts to capture the properties of an infinite structure with a finite lattice. Thus, the truncated lattice must be large enough to display the long-range order characteristics of the QC.
(2)
The reciprocal space of QCs is densely filled [19]; thus, the sampling of the FT must be performed with a very fine resolution to capture all of the diffraction spots.
Besides the fact that this approach is often computationally expensive, it can produce inaccurate results since it attempts to capture the properties of an infinite structure by only considering a finite subset of it. Another challenge that will be examined later is that using larger and larger lattices as indicated by Requirement #1 will lead to narrower, delta-like peaks in the reciprocal space that can only be detected by increasing the sampling rate in #2, which further compounds the computational burden of this approach. We will demonstrate this problem more clearly in the next section.
In this paper, we present a novel and rigorous method to generate the diffraction patterns of QCs based on the so-called higher dimensional approach, which captures all of the properties and symmetries of the entire infinite structure by a single unit cell in higher dimensional space [20]. The QC lattices are generated by taking 1D or 2D irrational slices of higher dimensional periodic hyper-lattices. Formally, this method is referred to as the cut-and-project method [20]. The dimensionality of the hyper-lattice corresponds to the rotational symmetry of the QC. While mathematically more involved, this approach elegantly captures all of the symmetries of the QC lattice within a unit cell of the hyper-lattice [18].
Here, we present a rigorous method of obtaining the diffraction pattern of a QC based on the Fourier series (FS) expansion of the higher dimensional unit cell, which leads to discrete harmonics. Unlike the traditional methods, which apply a FT to a finite-sized supper-lattice, our method is capable of capturing the exact diffraction pattern corresponding to an infinite QC lattice. To demonstrate the validity of our approach, we present two examples. First, we consider a 1D Fibonacci QC, which can be obtained as an irrational slice of a 2D periodic lattice. Next, we consider a 2D Penrose QC, which can be obtained either as a 2D slice of a 5D hyper-cubic lattice [20] or as a 2D slice of a 4D root lattice A 4 [21]. Here, we employ the 4D root lattice in our derivations, since it reduces the computational burden of the problem.

2. Quasicrystals

As was noted earlier, crystals are identified by their discrete diffraction pattern. Mathematically, there is a close relationship between diffraction patterns and the FT of point lattices. The Dirac delta function introduces a convenient way to represent a set of scatterers as a summation of infinitesimal points. We start by considering Λ = x 1 , x 2 , as a discrete set of points in n-dimensional Euclidean space E n where x k is an n-dimensional vector representing the location of the k-th point in Λ. The FT of this lattice denoted by f ^ ( s ) is obtained from the following integral:
f ^ ( s ) = R n f ( x ) exp 2 π i x · s d x
where i = 1 . It can easily be shown that for an infinite periodic lattice, the FT in Equation (1) becomes a discrete periodic impulse train or in essence, it becomes a FS. A periodic point lattice is said to possess translational symmetry, since a translation mapping by any integer summation of lattice vectors will map the lattice onto itself. Another property associated with crystals is their rotational symmetry. The rotational symmetry of a crystal is defined in terms of the rotational symmetry of its diffraction pattern. Thus, if the diffraction pattern of a point set is unchanged by a 2 π n rotation, the point set is said to possess n-fold rotational symmetry. According to the crystallographic restriction theorem, rotational symmetries of order five and those greater than six are impossible for diffraction patterns of periodic lattices existing in 2D and 3D space [1].
The Fibonacci QC is an example of a 1D QC. A Fibonacci sequence can be obtained iteratively from the alphabet set A , B , using the substitution rule: A A B , B A . Traditionally, F n denotes the Fibonacci sequence obtained after n iterative substitutions starting with F 0 = B [18]:
F 0 = B , F 1 = A , F 2 = A B , F 3 = A B A ,
Generating a Fibonacci point lattice from a Fibonacci sequence is accomplished by assigning a thickness α, β to A and B, respectively, such that α β = τ where τ = 1 + 5 / 2 is the golden ratio [1] and placing points at all junctions. Figure 1 shows the normalized diffraction pattern of a Fibonacci lattice with 4181 point scatterers, which corresponds to F 18 , as defined in Equation (2). The diffraction pattern is obtained as the magnitude of the FT defined in Equation (1) where the reciprocal space region 10 s 10 has been sampled at 100,001 equidistant points.
The Penrose QC, which is obtained from the Penrose aperiodic tiling, is an example of a 2D QC. The Penrose aperiodic tiling was discovered by Sir Roger Penrose in 1974 [22]. The Penrose lattice prototiles consist of a narrow (vertex angles π / 4 and 4 π / 5 ), shown in Figure 2a, and a wide (vertex angles 2 π / 5 and 3 π / 5 ) rhombi tile, shown in Figure 2b, of identical sides. The tiling, shown in Figure 2c, is obtained based on specific matching rules, which requires aligning arrows and double arrows, while pacing prototiles next to each other. As noted earlier, the tiling is converted to a point lattice simply by placing points at the vertices of the tiling. Figure 3 shows the normalized FT (magnitude) of a Penrose lattice of 711 point scatterers obtained from a Penrose tiling with tile-side of a. The diffraction pattern displayed in Figure 3 has been sampled at 501 points along 4 π a < k x < 4 π a and 501 points along 4 π a < k y < 4 π a for a total of 501 2 251 × 10 3 points. As can be seen from Figure 3, the diffraction pattern displays five-fold rotational symmetry, which is forbidden in periodic lattices.

3. Higher Dimensional Approach

As noted earlier, the first observation of QCs was reported in X-ray diffraction patterns of rapidly solidified Al-Mn alloys, which displayed icosahedral symmetry, that is inconsistent with translational symmetry [2]. After this discovery, extensive studies took place to clarify the underlying aperiodic structure. Eventually, the most accurate models were obtained using the higher dimensional approach to derive the structure model of these icosahedral QCs [23], which completely explained their X-ray and neutron diffraction intensity. The atomic configuration was derived from a 3D section of periodic atomic surfaces in 6D space [24]. Thus, the so-called higher dimensional approach is the most rigorous method to study QCs.
The higher dimensional approach is based on the fact that a QC lattice in n-dimensional physical space E n with k-fold rotational symmetry ( k > n ) can be embedded in a k-dimensional hyper-lattice in E k [20]. The higher dimensional space is thus referred to as the embedding space. The cut-and-project method [20], which is based on the higher dimensional approach, thus generates a QC lattice in n-dimensional space as an irrational slice of an integer hyper-lattice in E k . In the cut-and-project method, the hyper-lattice is projected onto two orthogonal subspaces of dimensions n and k n . The n-dimensional space E is referred to as the parallel or external subspace and corresponds to the physical space. The ( k n ) -dimensional subspace E is referred to as the perpendicular or internal subspace and corresponds to the unphysical space. Here, it is important to note that for a QC with k-fold rotational symmetry, k is not always the minimum dimension for the embedding space. The minimum dimension of the embedding space for a QC with k-fold rotational symmetry is φ ( k ) , where φ ( k ) is Euler’s totient function, which is the number of all positive integers less than or equal to k that are relatively prime to k [18].
The higher dimensional approach elegantly restores all hidden symmetries of QCs into the unit cell of the hyper-lattice. This approach has been used to accurately calculate the band diagram of a Fibonacci photonic QC by applying Bloch’s theorem to the higher dimensional unit cell [25]. More recently, this approach has also been used to accurately and efficiently analyze diffraction from QC gratings by applying Floquet’s theorem to a higher dimensional unit cell [26]. In both cases, the results were obtained using the previously-discussed super-cell approximation approach with its stated shortcomings.
Traditionally, the diffraction patterns of QCs are obtained by first taking the FT of finite-sized point lattices and then sampling the reciprocal space to detect the peaks. We discussed the requirements for obtaining accurate diffraction patterns in the previous section and also briefly alluded to the fact that the true increased computational burden of using a larger lattice is much greater than initially expected. This is due to the fact that by definition, a crystal has a discrete diffraction pattern. By a discrete diffraction pattern, we essentially mean that the diffraction pattern is a summation of weighted Dirac delta functions. Thus, as the number of points in a crystalline (or quasicrystalline) lattice gets larger, the diffraction peaks become narrower and sharper, at infinity becoming delta functions, which in turn require sampling the reciprocal space at a higher rate to detect them, further compounding the computational burden of the problem.

4. Diffraction Patterns Based on Fourier Series Expansion of Higher Dimensional Atomic Surfaces

In [23], it was noted that the diffraction pattern of a QC can be regarded as the projection of the Fourier spectrum of the hyper-lattice onto the physical space. Here, we exploit this idea to obtain the diffraction pattern of a QC from the FS expansion of the higher dimensional unit cell rather than the FT of a finite-sized super-cell. We demonstrate the method for a Fibonacci QC and a Penrose QC. The validity of the results is confirmed by comparing them to the FT of very large periodic super-cells.

4.1. Fibonacci QC

The Fibonacci QC lattice can be obtained using the cut-and-project method by taking a 1D irrational slice of a 2D periodic lattice [1], such as that shown in Figure 4a. The primitive lattice vectors are unit vectors along e 1 and e 2 denoted respectively by a 1 and a 1 , as shown in Figure 4a, where without loss of generality, we have assumed that the lattice constant is unity. As a result, every lattice point can be represented by a set of integers. This is referred to as an integer lattice and denoted by I 2 p , where the subscript 2 corresponds to the dimensionality of the lattice. The Voronoï cell at the origin is highlighted and denoted by V ( 0 ) .
We start by defining a new orthonormal basis with unit vectors x ^ and y ^ directed along x and y, as shown in Figure 4a. As can be seen, the unit vectors x ^ and y ^ are obtained by rotating the unit vectors e ^ 1 and e ^ 2 by an angle φ, where tan ϕ = 1 / τ . It can easily be shown that the transformation from the e ^ 1 e ^ 2 basis to the x ^ y ^ basis can be obtained by the unitary mapping M :
M = κ τ κ κ κ τ
where κ = 1 1 + τ 2 . Using this mapping, the subspace x is referred to as the parallel or external subspace and is denoted by E . The parallel subspace corresponds to the physical space. Subspace y is referred to as the perpendicular or internal subspace and is denoted by E . The perpendicular subspace corresponds to the unphysical space. As shown in Figure 4a, for every lattice point, P, its projection onto the parallel subspace E is denoted by Π ( P ) , and its projection onto the perpendicular subspace E is denoted by Π ( P ) . A Fibonacci lattice can now be obtained by the projection of lattice points onto the parallel subspace E . However, not all lattice points are projected onto the parallel subspace, and a selection criterion is used to determine the set of lattice points that should be projected to obtain a Fibonacci QC lattice. The selection window is shown in Figure 4a as the region between two parallel lines above and below the x-axis. This region corresponds to lattice points who’s orthogonal projection onto the perpendicular subspace, denoted by Π ( P ) , falls inside the perpendicular projection of the Voronoï cell at the origin, which is denoted by Π ( V ( 0 ) ) in Figure 4a.
An alternative way to generate QC lattices using the cut-and-project method is by decorating the lattice points with an atomic surface (AS) [18,27]. ASs are projections of the higher dimensional unit cell onto the perpendicular subspace. The points of intersection of the decorated hyper-lattice with the parallel subspace yield the QC lattice. Figure 4b shows the decorated hyper-lattice for the Fibonacci QC. The ASs are line segments of length Δ = sin ϕ + cos ϕ .
Since the decorated hyper-lattice is periodic with lattice vectors a ^ 1 and a ^ 2 , it can be expanded in terms of a FS:
ρ ( e 1 , e 2 ) = G c G exp i G · r
where the sum is over all RVs G defined as:
G = m 1 a 1 * + m 2 a 2 * m 1 , m 2 Z
such that a 1 * and a 2 * are the primitive RVs, which are related to lattice vectors a ^ 1 and a ^ 2 by a m * · a ^ k = 2 π δ m k where δ m k is the Kronecker delta function and r = e 1 a ^ 1 + e 2 a ^ 2 is the position vector. It is customary to use an equal number of terms for each dimension, which we denote by N, leading to a total number of ( 2 N + 1 ) 2 terms in the FS. The value of the Fourier expansion coefficient c G in Equation (4) is given by:
c G = U ρ ( e 1 , e 2 ) exp i G · r d e 1 d e 2
where the integration is performed over the area of a unit cell. It is more convenient to evaluate the integral in Equation (6) in the x y basis using a change of variables:
c G = U ρ ˜ ( x , y ) e 2 i π κ m 1 τ x + m 1 y + m 2 x m 2 τ y d x d y
where ρ ˜ ( x , y ) is the AS in the x y basis and can be expressed as ρ ˜ ( x , y ) = δ ( x ) rect y Δ . Using this expression, we arrive at the following result for c G :
c G = Δ if G = 0 sin Δ π κ m 1 m 2 τ π κ m 1 m 2 τ if G 0
Equation (8) gives the value of the Fourier harmonic c G for the RVs G as defined in Equation (5). The RVs are defined with respect to the lattice vectors a ^ 1 and a ^ 2 . We can also express the RVs in terms of the parallel and perpendicular subspaces as G = G x ^ + G y ^ where G is the projection of G onto the parallel subspace and G is the projection of G onto the perpendicular subspace. Thus, for the RVs G as defined in Equation (5), G and G are:
G = 2 π κ m 1 τ + m 2 G = 2 π κ m 1 m 2 τ
As it was noted earlier, the parallel subspace corresponds to the physical space. Hence, in the same manner that the 1D QC lattice was obtained by taking a slice of the decorated hyper-lattice along the parallel subspace, we can obtain the diffraction pattern of the QC by taking a slice of the reciprocal space of the hyper-lattice. Thus, only the parallel components of the RVs and their corresponding Fourier harmonics are considered. Figure 5 shows the normalized diffraction pattern of the Fibonacci QC using this method. A FS expansion with 441 terms N = 10 has been applied to the hyper-lattice ASs, and the resulting Fourier harmonics have been plotted versus the parallel components of the RVs. The results shown in Figure 5 can be verified by comparing them to Figure 1, which obtained the diffraction pattern by taking the FT of a very large super-cell structure.
Figure 1 was obtained by evaluating the FT of a 2585 lattice at 100,001 points in the reciprocal space. As was noted, when taking the FT of a very large, but finite point lattice, it is critically important to use a very high resolution in sampling the reciprocal space, since otherwise, some strong peaks might be missed. This point is demonstrated in Figure 6. Both Figure 6a and Figure 6b correspond to the FT (magnitude) spectra of the 2585 Fibonacci lattice, which was considered previously. The only difference is that in Figure 6a, the reciprocal space has been sampled at 10,001 points, and in Figure 6b, it has been sampled at 20,001 points. As can be seen, both of these patterns lack many of the peaks shown Figure 1, which was sampled at 100,001 points.

4.2. Penrose QC

Using the cut-and-project method, the Penrose QC lattice can be obtained as a 2D slice of a 5D periodic hyper-lattice. Obtaining the Penrose QC as a 2D slice of the 5D space is known as the canonical projection, which means that the unit cell of the hyper-lattice is a 5D hyper-cube. Using this method, the hyper-cube points are mapped onto a 2D parallel space and a 3D perpendicular space. However, as was noted previously, the minimum dimension of the embedding space for a QC with k-fold rotational symmetry is φ ( k ) , where φ ( k ) is Euler’s totient function and φ ( 5 ) = 4 . Since lower dimensionality reduces the computational burden of the problem, we use the minimal space embedding to obtain our results. An important distinction between the 5D embedding and the 4D embedding is that in the case of the 4D embedding, the hyper-lattice is non-orthogonal. The 4D non-orthogonal hyper-lattice is known as the A 4 root lattice, and the 4D unit cell is a hyper-rhombohedral [20]. A complete discussion of the derivation of the 4D lattice is given in [18,21]. Here, we only state the results, which are important for this work. We consider the 4D Euclidean space with the standard orthonormal Euclidean basis e ^ 1 , e ^ 2 , e ^ 3 , e ^ 4 . The primitive lattice vectors for the A 4 root lattice are given by:
d 1 = 2 5 c 1 1 , s 1 , c 2 1 , s 2 T d 2 = 2 5 c 2 1 , s 2 , c 4 1 , s 4 T d 3 = 2 5 c 3 1 , s 3 , c 6 1 , s 6 T d 4 = 2 5 c 4 1 , s 4 , c 8 1 , s 8 T
where c j = cos 2 j π 5 and s j = sin 2 j π 5 . The parallel subspace E is the plane spanned by e ^ 1 e ^ 2 , and the perpendicular subspace E is the plane spanned by e ^ 3 e ^ 4 . The AS consists of four pentagons in the perpendicular subspace E , which we denote by A , B , C , and D. The A , B , C , and D of the AS are respectively centered at points d s / 5 , 2 d s / 5 , 3 d s / 5 , and 4 d s / 5 , where d s = d 1 + d 2 + d 3 + d 4 is the body diagonal of the unit cell. The mapping matrix M from the lattice basis to the Euclidean basis is:
M = 2 5 c 1 1 c 2 1 c 3 1 c 4 1 s 1 s 2 s 3 s 4 c 2 1 c 4 1 c 6 1 c 8 1 s 2 s 4 s 6 s 8
Figure 7 shows the projection of the hyper-rhombohedral unit cell onto the perpendicular space E . Each vertex is denoted by its coordinate in the d 1 , d 2 , d 3 , d 4 basis. The four ASs have been highlighted and the perpendicular space basis vectors e ^ 3 and e ^ 4 have been shown for reference. The radii of the pentagons are [18]:
R A = R D = 2 τ 2 5 , R B = R C = 2 τ 5
Taking a 2D slice along the parallel space of the 4D hyper-lattice will generate a Penrose QC lattice, where the tile side is 2 τ 2 / 5 [18].
At this point, it is desired to obtain the FS expansion of our higher dimensional unit cell. To accomplish this goal, the periodic AS can be expanded in terms of Fourier harmonics:
ρ ( r 1 , r 2 , r 3 , r 4 ) = G c G exp i G · r
where the sum is over all RVs G defined as:
G = m 1 d 1 * + m 2 d 2 * + m 3 d 3 * + m 4 d 4 * m 1 , , m 4 Z
and r is the 4D position vector given by:
r = r 1 d 1 + r 2 d 2 + r 3 d 3 + r 4 d 4 d 1 , , d 4 R
As before, for practical reasons, the infinite expansion is truncated to the finite set of RVs defined as:
G = m 1 d 1 * + m 2 d 2 * + m 3 d 3 * + m 4 d 4 * : m 1 , m 4 N
which leads to a total of 2 N + 1 4 terms. The value of the Fourier harmonic coefficient c G in Equation (4) can be determined from:
c G = 1 det M ρ ( r 1 , r 2 , r 3 , r 4 ) exp i G · r d r 1 d r 2 d r 3 d r 4
where the integration is performed over the 4D A 4 root lattice. The term det M corresponds to the volume of the unit cell, which due to non-orthogonality, is not unity and must be accounted for. It is more convenient to evaluate the integral in Equation (17) in the e ^ 1 , e ^ 2 , e ^ 3 , e ^ 4 basis. Applying a change of variables to Equation (17), we arrive at:
c G = ρ ˜ ( e 1 , e 2 , e 3 , e 4 ) e 2 π i m 1 d 1 + m 2 d 2 + m 3 d 3 + m 4 d 4 d e 1 d e 2 d e 3 d e 4
where ρ ˜ ( e 1 , e 2 , e 3 , e 4 ) is the AS in the e ^ 1 , e ^ 2 , e ^ 3 , e ^ 4 basis. As shown in Figure 7, the AS surface of the 4D hyper-lattice consists of four pentagons, which lie in the perpendicular space E . The centers for ASs A , B , C and D in the e ^ 1 , e ^ 2 , e ^ 3 , e ^ 4 basis, respectively, are ϑ , 0 , ϑ , 0 , 2 ϑ , 0 , 2 ϑ , 0 , 3 ϑ , 0 , 3 ϑ , 0 and 4 ϑ , 0 , 4 ϑ , 0 , where ϑ = 2 5 . The AS of the 4D hyper-lattice can be represented as the sum of four pentagonal ASs, as shown in Figure 7. The linearity property of the FS expansion can be exploited to evaluate the expansion for each ASs separately. Evaluating individual FS terms requires performing double integration over the area of a pentagon. While this could be performed numerically, we were able to derive analytical expressions, which improved the calculation time by almost three orders of magnitude. The derivations for the analytical results are included in the Appendix. Using these results, we can now evaluate the Fourier harmonic coefficient c G for each RV G in closed-form. However, as was done in the case of the Fibonacci QC, we only need to consider the parallel component of each RV, which corresponds to the physical space. For RV G as defined in Equation (14), the parallel space component, denoted by G , is a 2D vector in the e 1 e 2 plane given by:
G = 2 π 5 m 1 c 1 + m 2 c 2 + m 3 c 3 + m 4 c 4 m 1 s 1 + m 2 s 2 + m 3 s 3 + m 4 s 4
such that the first element of the G corresponds to projection along e ^ 1 and the second element corresponds to projection along e ^ 2 , where at this point, we can associate e ^ 1 and e ^ 2 with the Cartesian dimensions of the physical space.
As noted earlier, a QC possess a discrete diffraction pattern. Figure 8a shows the magnitude of the normalized FT of a Penrose lattice with 731 points obtained from a Penrose tiling with prototile sides of 2 τ 2 / 5 . Moreover, this prototile side corresponds to the lattice obtained by taking a 2D slice of the 4D hyper-lattice. As can be seen from Figure 8a, while there are several bright spots in the pattern, the pattern is not truly discrete, since there are ripples surrounding the bright spots. This is due to the finite size of the super-cell lattice. Figure 8b shows the magnitude of the normalized FT of a Penrose lattice with 1961 points. Compared to Figure 8a, the diffraction pattern in Figure 8b displays a much better discrete nature. However, while increasing the size of the super-lattice increases the accuracy of the diffraction pattern, it will also require much finer sampling in the reciprocal space since the diffraction peaks will become more delta-like with a very narrow width. Thus, a very fine sampling of the reciprocal space is required to capture all of the diffraction spots, which further increases the computational burden of the problem. The diffraction pattern in Figure 8b was obtained by sampling 3001 points along each dimension of the reciprocal space for a total of roughly nine million points.
Figure 9 shows the normalized diffraction pattern of the Penrose QC obtained from the FS expansion of the 4D hyper-lattice with 279,841 terms ( N = 11 ) . Here, we note some important details regarding Figure 9. First, due to the very large number Fourier harmonics, we have performed a thresholding and only plotted Fourier harmonics above a certain threshold. We have used the threshold of 40 dB in Figure 9. Secondly, since it is extremely difficult to discern delta-like peaks in the diffraction pattern due to their infinitesimal width, to better visualize the discrete peaks, we have used Gaussian pulses of the form A e ( α R 2 ) / A , where A is the amplitude of the Fourier harmonic expansion coefficient, R is the Euclidean distance from the location of the Fourier harmonic in the reciprocal space and α is some arbitrary positive decay constant.
Comparing Figure 8b with Figure 9, we can see that many of their features are identical. Both figures have been normalized, and as a result, it is difficult to discern peaks with a smaller magnitude. To overcome this issue, instead of plotting the normalized magnitude of the diffraction peaks, we plot the square root of the normalized magnitude. This non-linear mapping has the effect of enhancing smaller peaks, which allows for a better comparison of finer details. This comparison is illustrated in Figure 10. The square root of the normalized FS expansion harmonics obtained from the expansion of the ASs in the 4D hyper-lattice are displayed in Figure 10a, while Figure 10b plots the square root of the normalized FT of a Penrose super-lattice. As can be seen, these plots display many of the finer features not visible in Figure 8b and Figure 9. One of the most striking features seen in Figure 10a, which is not present in Figure 10b, is the discrete nature of the diffraction pattern. This phenomenon is particularly noticeable in the region surrounding the origin, where in Figure 10a, a series of discrete peaks are present, whereas in Figure 10b, a series of ripples surround the origin. The presence of ripples surrounding the bright spots in FT patterns can lead to inaccurate results in several ways. First, they can mask weaker diffraction peaks, as can be seen in Figure 10b, in the region surrounding the origin. Moreover, they also can lead to the appearance of spurious spots in the pattern (in-phase addition).

5. From the Vertex Model to Natural QCs

During our derivations, we assumed identical particles at all vertices. We call this the vertex model. This assumption is valid for many optical and photonic structures based on QC morphologies [12,14,16,17]. However, in natural QCs, this is not the case, and different atoms are placed at different vertices. Hence, an important question arises, which is, can this methodology be extended to natural QCs where the vertex model is no longer valid? Here, we show that our approach can be generalized to all types of QCs where ASs are decorated with different atoms.
To demonstrate the validity of our approach, let us consider the Penrose lattice with icosahedral symmetry. The Penrose lattice prototiles consist of a narrow rhombus tile (vertex angles π / 4 and 4 π / 5 ) shown in Figure 2a and a wide rhombus tile (vertex angles 2 π / 5 and 3 π / 5 ) shown in Figure 2b of identical sides. The tiling, shown in Figure 2c, is obtained by placing prototiles next to each other based on specific matching rules, which requires aligning arrows and double arrows. A close inspection of Figure 2c shows eight distinct types of vertices in the tiling denoted by D 1 3 , J 2 5 , K 4 4 , T 1 6 , V 2 7 , S 5 5 , S 5 and Q 3 3 , as shown in Figure 11. For each vertex, the subscript corresponds to the number of double arrows surrounding it, and the superscript corresponds to the number of tiles surrounding it [28]. In fact, each of these vertices can be associated with a specific portion of the four pentagonal atomic surfaces, as shown in Figure 7. AS segments corresponding to each vertex type are shown in Figure 12, and the frequency of each vertex is directly proportional to its corresponding atomic surface [28].
In our vertex model, we assume identical particles at all vertices. Hence, to evaluate c G in Equation (18), the integration was performed uniformly over all ASs. As noted before, this is valid for many optical and photonic structures based on QC morphologies; however, in natural QCs, this is not the case, and different atomic surfaces are decorated with different atoms. To account for different atoms at different vertices, the integration in Equation (18) must include a weighing factor for different regions of ASs corresponding to different vertices. While this has not been completely investigated, we believe that the weighing factor should be proportional to the mass of the atom occupying the atomic surface.

6. Conclusions

We have presented a rigorous method to obtain and analyze the diffraction properties of QC lattices. Previously, the diffraction patterns of QCs have been obtained by taking the FT of very large super-lattices. The first problem with using super-lattices is that they attempt to capture the properties of a QC by only considering a finite subset of an infinite structure. As a result, very large super-lattices are required to produce accurate results, which significantly increases the computational burden of the problem. Furthermore, as the size of the super-lattice is increased, the diffraction pattern becomes more discrete with delta-like peaks. Thus, the sampling in the reciprocal space must be performed with a very fine resolution to capture all of the diffraction spots.
Our method is based on the cut-and-project technique, which constructs QC lattices as lower dimensional slices of higher dimensional periodic hyper-lattices. The periodicity of the higher dimensional hyper-lattice makes it possible to expand the associated unit cell in terms of Fourier harmonics, which leads to a discrete spectrum. Unlike the super-lattice approach, our method is exact and captures all of the properties of the complete infinite QC structure.
The utility of the proposed method was demonstrated for two important QCs; the 1D Fibonacci and the 2D Penrose QCs. In both cases, we showed that it was possible to derive analytical expressions for the Fourier harmonics of the higher dimensional ASs. These findings were validated by comparing them to results obtained via the conventional approach based on evaluating the FT of very large hyper-lattices. Besides the large computational burden associated with applying a FT to super-lattices, we also demonstrated that such methods can lead to inaccurate results due to the continuous nature of the FT, which can generate spurious modes (due to the coherent addition of ripples) or null some modes (due to the destructive addition of ripples). Finally, we proposed how our approach can be generalized from the vertex model, which assumed identical particles at all vertices, to a more realistic case where different ASs are occupied by different atoms.

Acknowledgments

This work was partially supported by the Penn State MRSEC-Center for Nanoscale Science NSF DMR-1420620.

Author Contributions

Farhad A. Namin conceived of the idea of deriving the diffraction pattern of QC lattices based on the FS of higher dimensional periodic lattices, which leads to a discrete spectrum, rather than traditional FT. Farhad A. Namin developed the theoretical results, which required projecting higher dimensional RVs onto the the physical space. Douglas H. Werner supervised and supported the project. Both authors co-wrote the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix: Fourier Series Expansion of a Pentagonal Atomic Surface in the A4 Root Lattice

We start by considering a pentagonal AS with its center at ( 0 , 0 , 0 , 0 ) and radius of R. The geometry of the AS in the perpendicular subspace is shown in Figure A1. The Fourier series expansion coefficients are determined from the integral in Equation (17). The integral in Equation (17) is written in terms of the lattice basis vectors d 1 , d 2 , d 3 , d 4 . It is more convenient to evaluate the integral in the Euclidean basis e ^ 1 , e ^ 2 , e ^ 3 , e ^ 4 . To evaluate the integral in the Euclidean basis, we employ the mapping:
r 1 , r 2 , r 3 , r 4 T = M 1 e 1 , e 2 , e 3 , e 4 T
where matrix M is shown in Equation (11). The atomic surface in the Euclidean basis is given by:
ρ ˜ ( e 1 , e 2 , e 3 , e 4 ) = δ ( e 1 , e 2 ) 1 e 3 , e 4
where 1 e 3 , e 4 is defined as:
1 e 3 , e 4 = 1 if e 3 , e 4 0 o.w.
and ⬠ represents the pentagon shown in Figure A1. We can now express the exponent in Equation (17) in terms of its parallel and perpendicular subspace components:
G · r = G · r + G · r
The parallel subspace component of Equation (A4) is of trivial importance. However, the perpendicular component is required to evaluate the Fourier expansion harmonics. Using Equation (A1), the expression for G · r is found to be:
G · r = 2 π ζ 1 e 3 + ζ 2 e 4
where:
ζ 1 = m 1 τ 2 5 + m 2 2 5 τ + m 3 2 5 τ m 4 τ 2 5 ζ 2 = m 1 2 5 4 τ m 2 τ 2 5 4 + m 3 τ 2 5 4 m 4 2 5 4 τ
Substituting all of these results into Equation (17), we arrive at:
c G = δ ( e 1 , e 2 ) exp i G · r d e 1 d e 2 1 e 3 , e 4 exp i G · r d e 3 d e 4
Figure A1. AS is a pentagon of radius R in the perpendicular subspace.
Figure A1. AS is a pentagon of radius R in the perpendicular subspace.
Crystals 06 00093 g013
Using the sifting property of the delta function, it follows that the value of the first double integral in Equation (A7) is simply unity. The second integral involves integrating exp i G · r over the pentagon ⬠, as shown in Figure A1. After some manipulation, we arrive at the following closed-from expression, which we denote by c G 0 as:
c G 0 R = 5 R 2 2 sin 2 π 5 if ζ 1 = ζ 2 = 0 F 1 ( ζ 1 , ζ 2 ) + F 2 ( ζ 1 , ζ 2 ) o.w.
where F 1 ( ζ 1 , ζ 2 ) and F 2 ( ζ 1 , ζ 2 ) are given by:
F 1 ( ζ 1 , ζ 2 ) = e i π ζ 1 R / τ γ ζ 2 cos π ζ 2 R β e i 5 π ζ 1 R cos π ζ 2 R α + 5 i ζ 1 sin π ζ 2 R β e i 5 π ζ 1 R sin π ζ 2 R α 2 ζ 2 π 2 5 ζ 1 2 + 2 5 5 ζ 2 2
F 2 ( ζ 1 , ζ 2 ) = e i τ + 1 π ζ 1 R ψ ζ 2 e i π ζ 1 R / τ + 5 i ζ 1 e 2 i π ζ 1 R sin η π ζ 2 R + ψ ζ 2 e 2 i π ζ 1 R cos η π ζ 2 R 2 ζ 2 π 2 5 ζ 1 2 5 + 2 5 ζ 2 2
and the constants α , β , γ , η , and ψ are defined as:
α = 2 + 2 5 τ 1 2 5 β = 2 + 2 5 + 1 τ 1 2 5 (A11) γ = 5 5 2 5 η = 1 τ 5 + 2 5 (A12) ψ = 5 5 + 2 5
Comparing the orientation of the pentagon in Figure A1 and the four ASs shown in Figure 7, we note that ASs A and C have the same orientation as Figure A1; however, they are shifted by vectors v A = ϑ , 0 , ϑ , 0 T and v C = 3 ϑ , 0 , 3 ϑ , 0 T , where ϑ = 2 5 . The FS expansion of shifted ASs may be obtained by employing the well-known shifting property:
ρ r c G ρ r v c G e i G · v
ASs B and D are obtained by rotating the AS shown in Figure A1 by an angle of π in the e 3 e 4 plane. The expansion coefficient for ASs B and D, which we denote as c G π , can be obtained using the rotation property of a FS. The superscript π corresponds to the fact that the ASs B and D have been rotated by an angle of π with respect to the AS of Figure A1, whose expansion coefficient is denoted by c G 0 . The shifting vectors for ASs B and D are v B = 2 ϑ , 0 , 2 ϑ , 0 T and v D = 4 ϑ , 0 , 4 ϑ , 0 T . Combining all of the results together, we can express the FS expansion coefficients of the 4D hyper-lattice as:
c G = c G 0 R A e i G · v A + c G π R B e i G · v B + c G 0 R C e i G · v C + c G π R D e i G · v D
where the values for R A , R B , R C and R D are given in Equation (12).

References

  1. Senechal, M. Quasicrystals and Geometry; Cambridge University Press: Cambridge, UK, 1996. [Google Scholar]
  2. Shechtman, D.; Blech, I.; Gratias, D.; Cahn, J.W. Metallic phase with long-range orientational order and no translational symmetry. Phys. Rev. Lett. 2013, 53, 1951–1953. [Google Scholar] [CrossRef]
  3. Gopinath, A.; Boriskina, S.V.; Reinhard, B.M.; Negro, L.D. Deterministic aperiodic arrays of metal nanoparticles for surface-enhanced Raman scattering (SERS). Opt. Express 2009, 17, 3741–3753. [Google Scholar] [CrossRef] [PubMed]
  4. Villa, A.D.; Galdi, V.; Capolino, F.; Pierro, V.; Enoch, S.; Tayeb, G. A comparative study of representative categories of EBG dielectric quasi-crystals. IEEE Antennas Wirelss Propagat. Lett. 2006, 5, 331–334. [Google Scholar] [CrossRef]
  5. Wang, X.; Werner, D.; Turpin, J. Investigation of scattering properties of large-scale aperiodic tilings using a combination of the characteristic basis function and adaptive integral methods. IEEE Trans. Antenn. Propag. 2013, 61, 3149–3160. [Google Scholar] [CrossRef]
  6. Mittra, R. (Ed.) Computational Electromagnetics: Recent Advances and Engineering Applications; Springer: New York, NY, USA, 2014; Chapter 12.
  7. Spence, T.G.; Werner, D.H. Design of broadband planar arrays based on the optimization of aperiodic tilings. IEEE Trans. Antenn. Propag. 2008, 56, 76–86. [Google Scholar] [CrossRef]
  8. Namin, F.; Petko, J.S.; Werner, D.H. Analysis and design optimization of robust aperiodic micro-UAV swarm-based antenna arrays. IEEE Trans. Antenn. Propag. 2012, 60, 2295–2308. [Google Scholar] [CrossRef]
  9. Gross, F.B. (Ed.) Frontiers in Antennas: Next Generation Design & Engineering; McGraw-Hill: New York, NY, USA, 2011; Chapter 1.
  10. Gopinath, A.; Boriskina, S.V.; Feng, N.; Reinhard, B.M.; Negro, L.D. Photonic-plasmonic scattering resonances in deterministic aperiodic structures. Nano Lett. 2008, 8, 2423–2431. [Google Scholar] [CrossRef] [PubMed]
  11. Namin, F.; Wang, X.; Werner, D.H. Reflection and transmission coefficients for finite-sized aperiodic aggregates of spheres. J. Opt. Soc. Am. B 2013, 30, 1008–1016. [Google Scholar] [CrossRef]
  12. Namin, F.A.; Yuwen, Y.A.; Liu, L.; Panaretos, A.H.; Werner, D.H.; Mayer, T.S. Efficient design , accurate fabrication and effective characterization of plasmonic quasi-crystalline arrays of nano- spherical particles 2016. Sci. Rep. 2016, 6, 22009. [Google Scholar] [CrossRef] [PubMed]
  13. Pala, R.A.; White, J.; Barnard, E.; Liu, J.; Brongersma, M.L. Design of plasmonic thin-film solar cells with broadband absorption enhancements. Adv. Mat. 2009, 21, 3504–3509. [Google Scholar] [CrossRef]
  14. Bauer, C.; Kobiela, G.; Giessen, H. 2D quasiperiodic plasmonic crystals. Sci. Rep. 2012, 2, 681. [Google Scholar] [CrossRef] [PubMed]
  15. Kaliteevski, M.A.; Brand, S.; Abram, R.A.; Krauss, T.F.; Millar, P.; Rue, R.M.D.L. Diffraction and transmission of light in low-refractive index Penrose-tiled photonic quasicrystals. J. Phys. Condens. Matter 2001, 13, 10459–10470. [Google Scholar] [CrossRef]
  16. Matsui, T.; Agrawal, A.; Nahata, A.; Vardeny, Z.V. Transmission resonances through aperiodic arrays of subwavelength apertures. Nature 2007, 446, 517–521. [Google Scholar] [CrossRef] [PubMed]
  17. Agrawal, A.; Matsui, T.; Vardeny, Z.V.; Nahata, A. Terahertz transmission properties of quasiperiodic and aperiodic aperture arrays. J. Opt. Soc. Am. B 2007, 24, 2545–2555. [Google Scholar] [CrossRef]
  18. Steurer, W.; Deloudi, S. Crystallography of Quasicrystals: Concepts, Methods and Structures; Springer: Berlin, Germany, 2009. [Google Scholar]
  19. Achanta, V.G. Plasmonic quasicrystals. Prog. Quant. Electron. 2015, 39, 1–23. [Google Scholar] [CrossRef]
  20. Janssen, T.; Chapuis, G.; de Boissieu, M. Aperiodic Crystals: From Modulated Phases to Quasicrystals; Oxford University Press: Oxford, UK, 2007. [Google Scholar]
  21. Baake, M.; Kramer, P.; Schlottmann, M.; Zeidler, D. Planar patterns with fivefold symmetry as sections of periodic structures in 4-space. Int. J. Mod. Phys. B 1990, 04, 2217–2268. [Google Scholar] [CrossRef]
  22. Penrose, R. The role of aesthetics in pure and applied mathematical research. Inst. Math. Appl. Bull. 1974, 10, 266–271. [Google Scholar]
  23. Yamamoto, A.; Hiraga, K. Structure of an icosahedral Al-Mn quasicrystal. Phys. Rev. B 1988, 37, 6207. [Google Scholar] [CrossRef]
  24. Ishihara, K.N.; Yamamoto, A. Penrose patterns and related structures. I. Superstructure and generalized Penrose patterns. Acta Crystallogr. Sect. A 1988, 44, 508–516. [Google Scholar] [CrossRef]
  25. Rodriguez, A.W.; McCauley, A.P.; Avniel, Y.; Johnson, S.G. Computation and visualization of photonic quasicrystal spectra via Bloch’s theorem. Phys. Rev. B 2008, 77, 104201. [Google Scholar] [CrossRef]
  26. Namin, F.A.; Werner, D.H. Rigorous analysis of diffraction from quasicrystalline gratings via Floquet’s theorem in higher-dimensional space. ACS Photonics 2014, 1, 212–220. [Google Scholar] [CrossRef]
  27. Dubois, J. Useful Quasicrystals; World Scientific: Hackensack, NJ, USA, 2005. [Google Scholar]
  28. Pavlovitch, A.; Kléman, M. Generalised 2D Penrose tilings: structural properties. J. Phys. A 1987, 20, 687. [Google Scholar] [CrossRef]
Figure 1. Normalized FT (magnitude) of a Fibonacci lattice with 4181 points. The reciprocal space region 10 s 10 has been sampled at 100,001 points.
Figure 1. Normalized FT (magnitude) of a Fibonacci lattice with 4181 points. The reciprocal space region 10 s 10 has been sampled at 100,001 points.
Crystals 06 00093 g001
Figure 2. (a) The narrow prototile (vertex angles π / 4 and 4 π / 5 ) and (b) the wide prototile (vertex angles 2 π / 5 and 3 π / 5 ) of the Penrose tiling; (c) the tiling obtained based on specific matching rules, which requires aligning arrows and double arrows, while placing prototiles next to each other.
Figure 2. (a) The narrow prototile (vertex angles π / 4 and 4 π / 5 ) and (b) the wide prototile (vertex angles 2 π / 5 and 3 π / 5 ) of the Penrose tiling; (c) the tiling obtained based on specific matching rules, which requires aligning arrows and double arrows, while placing prototiles next to each other.
Crystals 06 00093 g002
Figure 3. The normalized FT (magnitude) of the Penrose point lattice.
Figure 3. The normalized FT (magnitude) of the Penrose point lattice.
Crystals 06 00093 g003
Figure 4. 2D embedding of the Fibonacci quasicrystal (QC) showing (a) the selection window for the 2D hyper-lattice projection and (b) the decorated hyper-lattice of a Fibonacci QC with atomic surfaces (ASs) of length Δ.
Figure 4. 2D embedding of the Fibonacci quasicrystal (QC) showing (a) the selection window for the 2D hyper-lattice projection and (b) the decorated hyper-lattice of a Fibonacci QC with atomic surfaces (ASs) of length Δ.
Crystals 06 00093 g004
Figure 5. Diffraction pattern of the Fibonacci QC obtained by Fourier series (FS) expansion of the higher dimensional unit cell, as shown in Figure 4b. Fourier harmonics (magnitude) are plotted versus the parallel component of the corresponding reciprocal vectors (RVs).
Figure 5. Diffraction pattern of the Fibonacci QC obtained by Fourier series (FS) expansion of the higher dimensional unit cell, as shown in Figure 4b. Fourier harmonics (magnitude) are plotted versus the parallel component of the corresponding reciprocal vectors (RVs).
Crystals 06 00093 g005
Figure 6. Normalized FT (magnitude) of a Fibonacci lattice with 2585 points. The reciprocal space region 10 s 10 has been sampled at 10,001 points in (a) and at 20,001 points in (b).
Figure 6. Normalized FT (magnitude) of a Fibonacci lattice with 2585 points. The reciprocal space region 10 s 10 has been sampled at 10,001 points in (a) and at 20,001 points in (b).
Crystals 06 00093 g006
Figure 7. Projection of 4D hyper-rhombohedral unit cell onto the perpendicular space E . The four ASs have been highlighted, and the perpendicular space basis vectors e ^ 3 and e ^ 4 have been shown for reference.
Figure 7. Projection of 4D hyper-rhombohedral unit cell onto the perpendicular space E . The four ASs have been highlighted, and the perpendicular space basis vectors e ^ 3 and e ^ 4 have been shown for reference.
Crystals 06 00093 g007
Figure 8. Normalized FT (magnitude) of a Penrose lattice with 371 points in (a) and with 1961 points in (b).
Figure 8. Normalized FT (magnitude) of a Penrose lattice with 371 points in (a) and with 1961 points in (b).
Crystals 06 00093 g008
Figure 9. Diffraction pattern of the Penrose QC obtained from the FS expansion of a 4D hyper-lattice with 2 × 11 + 1 4 = 279,841 terms. The magnitude of the Fourier harmonics is plotted versus the parallel space component of the RVs.
Figure 9. Diffraction pattern of the Penrose QC obtained from the FS expansion of a 4D hyper-lattice with 2 × 11 + 1 4 = 279,841 terms. The magnitude of the Fourier harmonics is plotted versus the parallel space component of the RVs.
Crystals 06 00093 g009
Figure 10. (a) Square root of the normalized FS harmonics (magnitude) of the Penrose QC obtained using the expansion of the 4D hyper-lattice; (b) square root of the normalized FT (magnitude) of a Penrose super-lattice.
Figure 10. (a) Square root of the normalized FS harmonics (magnitude) of the Penrose QC obtained using the expansion of the 4D hyper-lattice; (b) square root of the normalized FT (magnitude) of a Penrose super-lattice.
Crystals 06 00093 g010
Figure 11. Eight distinct types of vertices in the Penrose tiling, as shown in Figure 2c, denoted by D 1 3 , J 2 3 , K 4 4 , T 1 6 , V 2 7 , S 5 5 , S 5 and Q 3 3 .
Figure 11. Eight distinct types of vertices in the Penrose tiling, as shown in Figure 2c, denoted by D 1 3 , J 2 3 , K 4 4 , T 1 6 , V 2 7 , S 5 5 , S 5 and Q 3 3 .
Crystals 06 00093 g011
Figure 12. (a) The portions of atomic surface D corresponding to vertices D 1 3 , S 5 , K 4 4 ; (b) portions of atomic surface C corresponding to vertices Q 3 3 , J 2 5 , T 1 6 , V 2 7 and S 5 . Those in B and A are related by an inversion center.
Figure 12. (a) The portions of atomic surface D corresponding to vertices D 1 3 , S 5 , K 4 4 ; (b) portions of atomic surface C corresponding to vertices Q 3 3 , J 2 5 , T 1 6 , V 2 7 and S 5 . Those in B and A are related by an inversion center.
Crystals 06 00093 g012

Share and Cite

MDPI and ACS Style

Namin, F.A.; Werner, D.H. An Exact Method to Determine the Photonic Resonances of Quasicrystals Based on Discrete Fourier Harmonics of Higher-Dimensional Atomic Surfaces. Crystals 2016, 6, 93. https://doi.org/10.3390/cryst6080093

AMA Style

Namin FA, Werner DH. An Exact Method to Determine the Photonic Resonances of Quasicrystals Based on Discrete Fourier Harmonics of Higher-Dimensional Atomic Surfaces. Crystals. 2016; 6(8):93. https://doi.org/10.3390/cryst6080093

Chicago/Turabian Style

Namin, Farhad A., and Douglas H. Werner. 2016. "An Exact Method to Determine the Photonic Resonances of Quasicrystals Based on Discrete Fourier Harmonics of Higher-Dimensional Atomic Surfaces" Crystals 6, no. 8: 93. https://doi.org/10.3390/cryst6080093

APA Style

Namin, F. A., & Werner, D. H. (2016). An Exact Method to Determine the Photonic Resonances of Quasicrystals Based on Discrete Fourier Harmonics of Higher-Dimensional Atomic Surfaces. Crystals, 6(8), 93. https://doi.org/10.3390/cryst6080093

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