Complete Band Gaps in Nano-Piezoelectric Phononic Crystals

Abstract

We study the band structure of elastic waves propagating in a nano-piezoelectric phononic crystal consisting of a polymeric matrix reinforced by BaTiO₃ inclusions in square, rectangular, triangular, honeycomb and Kagomé lattices. We also investigate the influence of inclusion cross section geometry - circular, hollow circular, square and rotated square with a 45º angle of rotation with respect to x and y axes. Plane wave expansion method is used to solve the governing equations of motion of a piezoelectric solid based on classical elasticity theory, ignoring nanoscopic size effects, considering two-dimensional periodicity and wave propagation in the xy plane. Complete band gaps between XY and Z modes are observed for all inclusions and the best performance is for circular inclusion in a triangular lattice. Piezoelectricity influences significantly the band gaps for hollow circular inclusion in lower frequencies. We suggest that nano-piezoelectric phononic crystals are feasible for elastic vibration management in GHz.

Keywords:
nano-piezoelectric phononic crystal; band structure; plane wave expansion method; complete band gaps; vibration control

1. Introduction

Phononic crystals (PCs) are artificial periodic composites designed to exhibit phononic band gaps and they have been quite studied^1-14. There are no mechanical (elastic or acoustic) propagating waves in phononic band gaps, only evanescent waves. These band gaps are created by the periodically mismatch between the constituent materials. This mismatch can be considered to arise either from difference of material properties or geometry (continuum-scale theory), or from interatomic force constants and masses (atomic-scale theory).

The ability of creating phononic band gaps is similar to the electronic and photonic band gaps in semiconductors and photonic crystals^15-16, respectively. The physical origin of phononic and photonic band gaps can be understood at micro-scale using the classical wave theory to describe the Bragg and Mie resonances, respectively, based on the scattering of mechanical and electromagnetic waves propagating within the crystal¹⁷.

PCs have many applications, such as vibration isolation technology^18-22, acoustic barriers/filters^23-25, noise suppression devices^26-27, surface acoustic devices²⁸, architectural design²⁹, sound shields³⁰, acoustic diodes³¹, elastic metamaterials^{21-22,25,27,32} and thermal metamaterials^33-39.

There are also smart PCs that have been studied, such as piezoelectric^40-54, piezomagnetic^55-58 and magnetoelectroelastic^14,59-63 PCs. Among these intelligent PCs, piezoelectric PCs are sensitive to elastic and electric field. Even though the band structure behavior of piezoelectric PCs have been already investigated, to our knowledge only few studies^46-49,51 focused on the influence of inclusion geometry and lattice on band gap formation. Wang and co-workers⁴⁶ considered a piezoelectric PC with square lattice and different inclusion geometries (regular triangle, square, hexagon, circle and oval). They found that the largest complete band gap is obtained by selecting the inclusion with the same symmetry of lattice for the first band gap. Hsu et al.⁴⁷ studied the band structure of a piezoelectric PC with square and triangular lattices using Mindlin-Reissner plate theory and considering only circular inclusion. Qian and co-workers⁴⁸ studied the band structure of piezoelectric PCs with square lattice and circular and square inclusions. They revealed the existence of several very large complete band gaps in PZT rods reinforced polythene composite. Zhou et al.⁴⁹ investigated the band structure of piezoelectric PCs consisting of rectangular inclusion in an epoxy substrate. They considered a PC with square and triangular lattices. Wang and co-workers⁵¹ considered a piezoelectric PC with circular inclusions and square and rectangular lattices. They verified that the largest band gap width is not always obtained for a square lattice. In some cases, rectangular lattice can produce broader gaps. In this context, we extend the studies about piezoelectric PCs^40-54 considering the influence of different inclusion geometries - circular, hollow circular, square and rotated square with 45º angle of rotation with respect to the x and y axes, and different lattices - square, rectangular, triangular, honeycomb and Kagomé on the band structure.

Furthermore, researches have used PCs on the scale µm^10,17,64,65 and mm, resulting in band gaps ranging from GHz and kHz to MHz, respectively. More recently, with the advance of nanomaterials fabrication, nanophononic crystals have been studied and it is possible to control wave propagation in a frequency range from hypersonic^3-6,66-75 to thermal^33-39. However, nano-piezoelectric PCs have not been investigated yet even though studies about nano-piezoelectric materials have been reported^76-79.

The main purpose of this study is to investigate the elastic band structure, also known as dispersion relation, of a nano-piezoelectric phononic crystal with square, rectangular, triangular, honeycomb and Kagomé lattices, composed by a polymeric matrix and BaTiO₃ inclusions, considering wave propagation in the xy plane and two-dimensional periodicity. We consider the following inclusions geometries - circular, hollow circular, square and rotated square with a 45º angle of rotation with respect to the x and y axes. To obtain the elastic band structure we use the semi-analytical plane wave expansion^1,2 (PWE) method.

To the best of our knowledge, it is the first time that the band structure of nano-piezoelectric PCs with different inclusion geometries and different lattices is investigated.

2. 2D Nano-Piezoelectric Phononic Crystal Model by PWE Method

PWE formulation of a 2D PC is presented in this section. PWE method^1,2, also known as ω(k) method, is one of the most used methods to calculate the elastic band structure of PCs and it has been applied in micro¹⁰ and nanophononic crystals^71,75 (NPCs). PWE method is used to solve the constitutive equations of a piezoelectric material based on classical elasticity theory, ignoring nanoscopic size effects, similar to other studies^10,71,75.

We consider two-dimensional periodicity, transversely isotropic elastic solid and wave propagation in the xy plane. Complete band gaps between XY (longitudinal-transverse vibration) and Z (transverse vibration) modes are observed for all inclusion geometry - circular, hollow circular, square and rotated square, considering square, rectangular, triangular, honeycomb and Kagomé lattices. Piezoelectric effect on the Z modes is also investigated.

2.1 2D Phononic Crystal Model

Figure 1 (a-e) sketches the cross section of nano-piezoelectric PC unit cells taking in to account square, rectangular, triangular, honeycomb and Kagomé lattices, respectively, with arbitrary inclusion geometry. Figure 1 (f-h) represents the first Brillouin zone⁸⁰ (FBZ) for square, rectangular, triangular, honeycomb and Kagomé lattices. Note that triangular, honeycomb and Kagomé lattices present the same FBZ shape (Figure 1 (h)). We consider four types of BaTiO₃ inclusions: circular, hollow circular, square and rotated square with a 45º angle of rotation with respect to the x, y axes. We underline that only BaTiO₃ inclusions are rotated for rotated square configuration, not the entire lattice.

Figure 1
Transverse cross-section of the nano-piezoelectric phononic crystal unit cell: BaTiO₃ inclusions distributed in a polymeric matrix for square (a), rectangular (b), triangular (c), honeycomb (d) and Kagomé (e) lattices. First Brillouin zone for square (f), rectangular (g), triangular (h), honeycomb (h) and Kagomé (h) lattices.

In addition, we highlight that we consider 1/4 of the FBZ as the actual first irreducible Brillouin zone⁸⁰ (FIBZ) for square and rotated square inclusions in triangular, honeycomb and Kagomé lattices because the symmetries of phononic crystals are reduced since the phononic crystal basis has different symmetry from the lattice⁸¹. Thus, for these configurations we scan M'-Γ-X-M-X'-M'. It is very important to make this adjustment since it changes the phononic crystal band structures⁸¹. For circular and hollow circular inclusions in triangular, honeycomb and Kagomé lattices, phononic crystal symmetry is not reduced thus we scan only M-Γ-X-M.

FIBZ points in Figure 1 (f-h) are Γ (0,0), X $\left(\frac{\mathit{\pi }}{\mathit{a}},0\right)$ and M $\left(\frac{\mathit{\pi }}{\mathit{a}},\frac{\mathit{\pi }}{\mathit{a}}\right)$ for square lattice, Γ (0,0), X $\left(\frac{\mathit{\pi }}{{\mathit{a}}_1},0\right)$, M $\left(\frac{\mathit{\pi }}{{\mathit{a}}_1},\frac{\mathit{\pi }}{{\mathit{a}}_2}\right)$ and K $\left(0,\frac{\mathit{\pi }}{{\mathit{a}}_2}\right)$ for rectangular lattice, Γ (0,0), X $\left(\frac{4\mathit{\pi }}{3\mathit{a}},0\right)$, M $\left(\frac{\mathit{\pi }}{\mathit{a}},\frac{\mathit{\pi }}{\sqrt{3\mathit{a}}}\right)$, X' $\left(\frac{2\mathit{\pi }}{3\mathit{a}},\frac{2\mathit{\pi }}{\sqrt{3\mathit{a}}}\right)$ and M' $\left(0,\frac{2\mathit{\pi }}{\sqrt{3\mathit{a}}}\right)$ for triangular lattice, Γ (0,0), X $\left(\frac{4\mathit{\pi }}{3\sqrt{3\mathit{a}}},0\right)$, M $\left(\frac{\mathit{\pi }}{\sqrt{3\mathit{a}}},\frac{\mathit{\pi }}{3\mathit{a}}\right)$, X' $\left(\frac{2\mathit{\pi }\sqrt{3}}{9\mathit{a}},\frac{2\mathit{\pi }}{3\mathit{a}}\right)$ and M' $\left(0,\frac{2\mathit{\pi }}{3\mathit{a}}\right)$ for honeycomb lattice, and Γ (0,0), X $\left(\frac{2\mathit{\pi }}{3\mathit{a}},0\right)$, M $\left(\frac{\mathit{\pi }}{2\mathit{a}},\frac{\mathit{\pi }}{2\sqrt{3\mathit{a}}}\right)$, X' $\left(\frac{\mathit{\pi }}{3\mathit{a}},\frac{\mathit{\pi }}{\sqrt{3\mathit{a}}}\right)$ and M' $\left(0,\frac{\mathit{\pi }}{\sqrt{3\mathit{a}}}\right)$ for Kagomé lattice, respectively, where a is the lattice parameter for square, triangular, honeycomb and Kagomé lattices, whereas a₁, a₂ are the lattice parameters for rectangular lattice.

The constitutive equations of an elastic piezoelectric material are⁴¹:

(1) ${\mathit{\sigma }}_{\mathit{ij}}={\mathit{c}}_{\mathit{ijkl}}{\mathit{u}}_{\mathit{k}.\mathit{l}}-{\mathit{e}}_{\mathit{lij}}{\mathit{E}}_{\mathit{l}},$

(2) ${\mathit{D}}_{\mathit{i}}={\mathit{e}}_{\mathit{ikl}}{\mathit{u}}_{\mathit{k}.\mathit{l}}-{\mathit{\varepsilon }}_{\mathit{il}}{\mathit{E}}_{\mathit{l}},$

where i,j,k,l = 1,2,3, σ_ij is the elastic stress tensor, D_i is the electric displacement vector, u_i is the elastic displacement vector, E_l is the electric field vector, c_ijkl is the elastic stiffness tensor, e_lij is the piezoelectric tensor and ϵ_il is the dielectric tensor. The standard tensor notation is used with Latin indices running from 1 to 3. They obey Einstein's summation convention when repeated.

We restricted the treatment to linear media, thus the elastic strain tensor ε_kl is simplified:

(3) ${\epsilon }_{\mathit{kl}}=\frac{1}{2}\left({\mathit{u}}_{\mathit{k},\mathit{l}}+{\mathit{u}}_{\mathit{l},\mathit{k}}\right).$

Furthermore, based on the quasi-static approximation, there is no electro source (i.e. since there is no voltage source, this configuration is also known as an open circuit) and the curl is zero, thus the electric field is taken as gradient of scalar potential and one can write:

(4) ${\mathit{E}}_{\mathit{l}}\Phi -{\mathit{\Phi }}_{,\mathit{l}}.$

The differential equations of motion in absence of body forces are given by:

(5) ${\mathit{\sigma }}_{\mathit{ij},\mathit{i}}=\mathit{\rho }{\ddot{\mathit{u}}}_{\mathit{j}},$

(6) ${\mathit{D}}_{\mathit{i},\mathit{i}}=0,$

where ρ is the mass density and dot denotes differentiation with respect to time. Substituting Eqs. (1-2) in Eqs. (5-6), applying the simplifications of Eqs. (3-4), considering a transversely isotropic elastic solid and for a two-dimensional problem, ∂/∂x₃=0, results:

(7) $\mathit{\rho }{\ddot{\mathit{u}}}_1={\left({\mathit{c}}_{11}{\mathit{u}}_{1,1}+{\mathit{c}}_{12}{\mathit{u}}_{2,2}\right)}_{,1}+{\left[{\mathit{c}}_{66}\left({\mathit{u}}_{1,2}+{\mathit{u}}_{2,1}\right)\right]}_{,2},$

(8) $\mathit{\rho }{\ddot{\mathit{u}}}_2={\left({\mathit{c}}_{12}{\mathit{u}}_{1,1}+{\mathit{c}}_{11}{\mathit{u}}_{2,2}\right)}_{,2}+{\left[{\mathit{c}}_{66}\left({\mathit{u}}_{1,2}+{\mathit{u}}_{2,1}\right)\right]}_{,1},$

(9) $\mathit{\rho }{\ddot{\mathit{u}}}_3={\left({\mathit{c}}_{44}{\mathit{u}}_{3,1}+{\mathit{e}}_{15}{\mathit{\Phi }}_{,1}\right)}_{,1}+{\left({\mathit{c}}_{44}{\mathit{u}}_{3,2}+{\mathit{e}}_{15}{\mathit{\Phi }}_{,2}\right)}_{,2},$

(10) $0={\left({\mathit{e}}_{15}{\mathit{u}}_{3,1}-{\mathit{\varepsilon }}_{11}{\mathit{\Phi }}_{,1}\right)}_{,1}+{\left({\mathit{e}}_{15}{\mathit{u}}_{3,2}-{\mathit{\varepsilon }}_{11}{\mathit{\Phi }}_{,2}\right)}_{,2},$

where c₆₆= $\frac{1}{2}$ (c₁₁ - c₁₂). Note that Eqs. (7-10) are written using Voigt notation and from now on this notation is adopted.

In addition, considering a piezoelectric PC, one can note that c₁₁ = c₁₁ (r), c₁₂ = c₁₂ (r), c₆₆ = c₆₆ (r), c₄₄ = c₄₄ (r), e₁₅ = e₁₅ (r), ϵ₁₁ = ϵ₁₁ (r), ρ = ρ(r), because we consider two different materials - BaTiO₃ inclusions and a polymeric matrix, and u_i = u_i (r,t). For a two-dimensional periodicity (the system has translational symmetry in z direction and material parameters depend only on the x and y coordinates), then r = xe₁ + ye₂ (x,y ∈ ℝ) is the two-dimensional spatial vector and e_i (i = 1,2) are the basis vectors in real space.

In order to eliminate the factor time in Eqs. (7-10), we apply the temporal Fourier transform. Applying Floquet-Bloch's theorem^82-83, expanding Bloch wave amplitude as Fourier series in reciprocal space and considering wave propagation in the xy plane (k₃ = 0), we can write:

(11) $\begin{array}{c}{\mathit{u}}_{\mathit{i}}\left(\mathbf{r}\right)={\mathit{e}}^{\mathit{j}\mathbf{k}.\mathbf{r}}{\mathit{u}}_{\mathit{i}\mathbf{k}}\left(\mathbf{r}\right)={\mathit{e}}^{\mathit{j}\mathbf{k}.\mathbf{r}}\hfill \\ \sum _{\mathbf{g}=-\infty }^{+\infty }{\mathit{u}}_{\mathit{i}\mathbf{k}}\left(\mathbf{g}\right){\mathit{e}}^{\mathit{j}\mathbf{g}.\mathbf{r}}=\hfill \\ \sum _{\mathbf{g}=-\infty }^{+\infty }{\mathit{u}}_{\mathit{i}\mathbf{k}}\left(\mathbf{g}\right){\mathit{e}}^{\mathit{j}\left(\mathbf{k}\mathit{+}\mathbf{g}\right).\mathbf{r}},\hfill \end{array}$

where u_ik (r) is the Bloch wave amplitude, note that u_ik (r) = u_ik (r + r) and u_i(r + r) = u_i(r)e^jk.r, e^jk.r is called Bloch periodic boundary condition, k = ub₁ + 𝑣b₂ is the Bloch wave vector, u, 𝑣 ∈ ℚ are the symmetry points within the FIBZ in reciprocal space, or we may write k = k₁ e₁ + k₂ e₂, k₁, k₂ ∈ ℝ are the point coordinates within the FIBZ in Figure 1 (f-h) for the reciprocal space. The basis vectors in reciprocal space b_i (i = 1,2) are defined as a_i ∙ b_j = 2πδ_ij, δ_ij = 0 if i ≠ j or δ_ij = 1 if i = j is the Kronecker delta, b₁ = 2π $\frac{{\mathbf{a}}_2\times {\mathbf{a}}_3}{{\mathbf{a}}_1·\left({\mathbf{a}}_2\times {\mathbf{a}}_3\right)}$b₂ = 2π $\frac{{\mathbf{a}}_3\times {\mathbf{a}}_1}{{\mathbf{a}}_2·\left({\mathbf{a}}_3\times {\mathbf{a}}_1\right)}$a_i (i = 1,2) are the components of the lattice vector r = (pa₁ + qa₂), p, q ∈ ℤ.

The lattice vector components are a_i = ae_i (i = 1,2) for square lattice, a₁ = a₁ e₁, a₂ = a₂ e₂ for rectangular lattice, a₁ = ae₁, a₂ = $\frac{\mathit{a}}{2}$e₁ + $\frac{\mathit{a}\sqrt{3}}{2}$e₂ for triangular lattice, a₁ = $\frac{\mathit{a}\sqrt{3}}{2}$e₁ + $\frac{3\mathit{a}}{2}$e₂, a₂ = -$\frac{\mathit{a}\sqrt{3}}{2}$e₁ + $\frac{3\mathit{a}}{2}$e₂ for honeycomb lattice and a₁ = ae₁ + $\mathit{a}\sqrt{3}$e₂, a₂ = -ae₁ + $\mathit{a}\sqrt{3}$e₂ for Kagomé lattice.

The reciprocal lattice vector is defined as $\mathbf{g}=\frac{2\mathit{\pi }}{\mathit{a}}\left(\mathit{m}{\mathbf{e}}_1+\mathit{n}{\mathbf{e}}_2\right)$ (m, n ∈ ℤ) for square lattice, $\mathbf{g}=2\mathit{\pi }\left(\frac{\mathit{m}}{{\mathit{a}}_1}{\mathbf{e}}_1+\frac{\mathit{n}}{{\mathit{a}}_2}{\mathbf{e}}_2\right)$ for rectangular lattice, $\mathbf{g}=\frac{2\mathit{\pi }}{\mathit{a}}\left[\mathit{m}{\mathbf{e}}_1+\frac{\left(-\mathit{m}+2\mathit{n}\right)}{\sqrt{3}}{\mathbf{e}}_2\right]$ for triangular lattice, $\mathbf{g}=\frac{2\mathit{\pi }}{\mathit{a}\sqrt{3}}\left[\left(\mathit{m}-\mathit{n}\right){\mathbf{e}}_1+\frac{\mathit{m}+\mathit{n}}{\sqrt{3}}{\mathbf{e}}_2\right]$ for honeycomb lattice and $\mathbf{g}=\frac{\mathit{\pi }}{\mathit{a}}\left[\left(\mathit{m}-\mathit{n}\right){\mathbf{e}}_1+\frac{\left(\mathit{m}+\mathit{n}\right)}{\sqrt{3}}{\mathbf{e}}_2\right]$ for Kagomé lattice. Note that g is a two-dimensional vector because we consider two-dimensional periodicity.

Furthermore, we may expand c₁₁,c₁₂,c₆₆,c₄₄,e₁₅,ϵ₁₁, ρ in Fourier series on the reciprocal space as:

(12) $\mathit{P}\left(\mathbf{r}\right)=\sum _{\stackrel{─}{\mathbf{g}}=-\infty }^{+\infty }\mathit{P}\left(\stackrel{─}{\mathbf{g}}\right){\mathit{e}}^{\mathit{j}\stackrel{─}{\mathbf{g}}.\mathbf{r}},$

where P is one of the c₁₁,c₁₂,c₆₆,c₄₄,e₁₅,ϵ₁₁, ρ and g has the same expressions of g with m, n ∈ ℤ. We use g instead of g to highlight the difference between the Fourier series expansion of material properties and displacements.

Substituting Eqs. (11) and (12) in Eqs. (7-10), with $\tilde{\mathbf{g}}$ = g + g, multiplying by e^-jg.r and integrating over the unit cell, we may write:

(13) $\left(\mathbf{K}-{\mathit{\omega }}^2\mathbf{M}\right)\mathbf{q}=\mathbf{0}$

where

(14) $\mathbf{K}=\left[\begin{array}{c}{\mathbf{K}}_{11}{\mathbf{K}}_{12}\mathbf{0}\mathbf{0}\hfill \\ {\mathbf{K}}_{21}{\mathbf{K}}_{22}\mathbf{0}\mathbf{0}\hfill \\ \mathbf{0}\mathbf{0}{\mathbf{K}}_{33}{\mathbf{K}}_{34}\hfill \\ \mathbf{0}\mathbf{0}{\mathbf{K}}_{43}{\mathbf{K}}_{44}\hfill \end{array}\right]$

The sub-matrices in Eq. (14) are given by:

(15) $\begin{array}{c}{\mathbf{K}}_{11}={\mathit{c}}_{11}\left(\stackrel{─}{\mathbf{g}}-\mathbf{g}\right){\left(k+g\right)}_1{\left(k+\stackrel{─}{g}\right)}_1+\hfill \\ {\mathit{c}}_{66}\left(\stackrel{─}{\mathbf{g}}-\mathbf{g}\right){\left(k+g\right)}_2{\left(k+\stackrel{─}{g}\right)}_2,\hfill \end{array}$

(16) $\begin{array}{c}{\mathbf{K}}_{12}={\mathit{c}}_{12}\left(\stackrel{─}{\mathbf{g}}-\mathbf{g}\right){\left(k+g\right)}_2{\left(k+\stackrel{─}{g}\right)}_1+\hfill \\ {\mathit{c}}_{66}\left(\stackrel{─}{\mathbf{g}}-\mathbf{g}\right){\left(k+g\right)}_1{\left(k+\stackrel{─}{g}\right)}_2,\hfill \end{array}$

(17) $\begin{array}{c}{\mathbf{K}}_{21}={\mathit{c}}_{12}\left(\stackrel{─}{\mathbf{g}}-\mathbf{g}\right){\left(k+g\right)}_1{\left(k+\stackrel{─}{g}\right)}_2+\hfill \\ {\mathit{c}}_{66}\left(\stackrel{─}{\mathbf{g}}-\mathbf{g}\right){\left(k+g\right)}_2{\left(k+\stackrel{─}{g}\right)}_1,\hfill \end{array}$

(18) $\begin{array}{c}{\mathbf{K}}_{22}={\mathit{c}}_{11}\left(\stackrel{─}{\mathbf{g}}-\mathbf{g}\right){\left(k+g\right)}_2{\left(k+\stackrel{─}{g}\right)}_2+\hfill \\ {\mathit{c}}_{66}\left(\stackrel{─}{\mathbf{g}}-\mathbf{g}\right){\left(k+g\right)}_1{\left(k+\stackrel{─}{g}\right)}_1,\hfill \end{array}$

(19) $\begin{array}{c}{\mathbf{K}}_{33}={\mathit{c}}_{44}\left(\stackrel{─}{\mathbf{g}}-\mathbf{g}\right){\left(k+g\right)}_1{\left(k+\stackrel{─}{g}\right)}_1+\hfill \\ {\mathit{c}}_{44}\left(\stackrel{─}{\mathbf{g}}-\mathbf{g}\right){\left(k+g\right)}_2{\left(k+\stackrel{─}{g}\right)}_2,\hfill \end{array}$

(20) $\begin{array}{c}{\mathbf{K}}_{34}={\mathit{e}}_{15}\left(\stackrel{─}{\mathbf{g}}-\mathbf{g}\right){\left(k+g\right)}_1{\left(k+\stackrel{─}{g}\right)}_1+\hfill \\ {\mathit{e}}_{15}\left(\stackrel{─}{\mathbf{g}}-\mathbf{g}\right){\left(k+g\right)}_2{\left(k+\stackrel{─}{g}\right)}_2,\hfill \end{array}$

(21) $\begin{array}{c}{\mathbf{K}}_{44}=-{\mathit{\varepsilon }}_{11}\left(\stackrel{─}{\mathbf{g}}-\mathbf{g}\right){\left(k+g\right)}_1{\left(k+\stackrel{─}{g}\right)}_1-\hfill \\ {\mathit{\varepsilon }}_{11}\left(\stackrel{─}{\mathbf{g}}-\mathbf{g}\right){\left(k+g\right)}_2{\left(k+\stackrel{─}{g}\right)}_2,\hfill \end{array}$

(22) ${\mathbf{K}}_{43}={\mathbf{K}}_{34}.$

The matrix M in Eq. (13) is expressed by:

(23) $\mathbf{M}=\left[\begin{array}{c}\mathit{\rho }\left(\stackrel{─}{\mathbf{g}}-\mathbf{g}\right)\mathbf{0}\mathbf{0}\mathbf{0}\hfill \\ \mathbf{0}\mathit{\rho }\left(\stackrel{─}{\mathbf{g}}-\mathbf{g}\right)\mathbf{0}\mathbf{0}\hfill \\ \mathbf{0}\mathbf{0}\mathit{\rho }\left(\stackrel{─}{\mathbf{g}}-\mathbf{g}\right)\mathbf{0}\hfill \\ \mathbf{0}\mathbf{0}\mathbf{0}\mathbf{0}\hfill \end{array}\right]$

The vector q in Eq. (13) is given by:

(24) $\mathbf{q}=\left\{\begin{array}{c}{\mathit{u}}_{1\mathbf{k}}\left(\mathbf{g}\right)\hfill \\ {\mathit{u}}_{2\mathbf{k}}\left(\mathbf{g}\right)\hfill \\ {\mathit{u}}_{3\mathbf{k}}\left(\mathbf{g}\right)\hfill \\ {\mathit{\Phi }}_{\mathbf{k}}\left(\mathbf{g}\right)\hfill \end{array}\right\}.$

Equation (13) is an infinite system of equations, thus the Fourier series needs to be truncated. We choose m,m,n,n = [-M, ..., M] and the total number of plane waves is (2M + 1)². Equation (13) represents a generalized eigenvalue problem of ω² (k) and should be solved for each k into the FIBZ for square, rectangular, triangular, honeycomb and Kagomé lattices (Figure 1(f-h)).

The Fourier coefficients are defined as:

(25) $\mathit{P}\left(\stackrel{─}{\mathbf{g}}\right)=\left\{\begin{array}{c}{\mathit{fP}}_{\mathit{A}}+\left(1-\mathit{f}\right){\mathit{P}}_{\mathit{B}}\mathrm{for}\stackrel{─}{\mathbf{g}}=\mathbf{0}\hfill \\ \left({\mathit{P}}_{\mathit{A}}-{\mathit{P}}_{\mathit{B}}\right)\mathit{F}\left(\stackrel{─}{\mathbf{g}}\right)\mathrm{for}\stackrel{─}{\mathbf{g}}\ne \mathbf{0}\hfill \end{array}\right.,$

where the indexes A and B in Eq. (25) are related to inclusion (BaTiO₃) and polymeric matrix, respectively, F(g) is the structure function and ƒ is the filling fraction of each type of inclusion, considering circular section of radius r̃, square section of width 2l, rotated square section of width 2l with a 45º angle of rotation with respect to x and y axes and hollow circular section with external radius R̃ and internal radius r̃, R̃ > r̃. The hollow cylinder inclusion has an internal radius r̃ of the same polymer at matrix and a thickness R̃ - r̃ of BaTiO₃, similar to Anjos & Arantes¹⁰.

The filling fraction is defined as ƒ = S_A/S_C for square, rectangular and triangular lattices, whereas for honeycomb and Kagomé lattices it is defined as ƒ = 2S_A/S_C and ƒ = 3S_A/S_C, respectively, where S_A is the cross section area of the BaTiO₃ inclusions and S_C = ‖a₁ × a₂‖ is the cross section area of the unit cell. Note that differently from other lattices honeycomb and Kagomé lattices have two and three inclusions per unit cell, respectively.

The filling fractions for square, rectangular, triangular, honeycomb and Kagomé lattices are:

(26) $\mathit{f}=\left\{\begin{array}{c}\mathit{\pi }{\tilde{\mathit{r}}}^2/{\mathit{a}}^2\mathrm{for}\mathrm{circular}\sec \mathrm{tion}\hfill \\ 4{\mathit{l}}^2/{\mathit{a}}^2\mathrm{for}\mathrm{square}\sec \mathrm{tion}\hfill \\ 4{\mathit{l}}^2/{\mathit{a}}^2\mathrm{for}\mathrm{rotated}\mathrm{square}\sec \mathrm{tion}\hfill \\ \mathit{\pi }\left({\tilde{\mathit{R}}}^2-{\tilde{\mathit{r}}}^2\right)/{\mathit{a}}^2\mathrm{for}\mathrm{hollow}\sec \mathrm{tion}\hfill \end{array}\right.,$

(27) $\mathit{f}=\left\{\begin{array}{c}\mathit{\pi }{\tilde{\mathit{r}}}^2/\left({\mathit{a}}_1{\mathit{a}}_2\right)\mathrm{for}\mathrm{circular}\sec \mathrm{tion}\hfill \\ 4{\mathit{l}}^2/\left({\mathit{a}}_1{\mathit{a}}_2\right)\mathrm{for}\mathrm{square}\sec \mathrm{tion}\hfill \\ 4{\mathit{l}}^2/\left({\mathit{a}}_1{\mathit{a}}_2\right)\mathrm{for}\mathrm{rotated}\mathrm{square}\sec \mathrm{tion}\hfill \\ \mathit{\pi }\left({\tilde{\mathit{R}}}^2-{\tilde{\mathit{r}}}^2\right)/\left({\mathit{a}}_1{\mathit{a}}_2\right)\mathrm{for}\mathrm{hollow}\sec \mathrm{tion}\hfill \end{array}\right.,$

(28) $\mathit{f}=\left\{\begin{array}{c}2\mathit{\pi }{\tilde{\mathit{r}}}^2/\sqrt{3}{\mathit{a}}^2\mathrm{for}\mathrm{circular}\sec \mathrm{tion}\hfill \\ 8{\mathit{l}}^2/\sqrt{3}{\mathit{a}}^2\mathrm{for}\mathrm{square}\sec \mathrm{tion}\hfill \\ 8{\mathit{l}}^2/\sqrt{3}{\mathit{a}}^2\mathrm{for}\mathrm{rotated}\mathrm{square}\sec \mathrm{tion}\hfill \\ 2\mathit{\pi }\left({\tilde{\mathit{R}}}^2-{\tilde{\mathit{r}}}^2\right)/\sqrt{3}{\mathit{a}}^2\mathrm{for}\mathrm{hollow}\sec \mathrm{tion}\hfill \end{array}\right.,$

(29) $\mathit{f}=\left\{\begin{array}{c}4\mathit{\pi }{\tilde{\mathit{r}}}^2/3\sqrt{3}{\mathit{a}}^2\mathrm{for}\mathrm{circular}\sec \mathrm{tion}\hfill \\ 16{\mathit{l}}^2/3\sqrt{3}{\mathit{a}}^2\mathrm{for}\mathrm{square}\sec \mathrm{tion}\hfill \\ 16{\mathit{l}}^2/3\sqrt{3}{\mathit{a}}^2\mathrm{for}\mathrm{rotated}\mathrm{square}\sec \mathrm{tion}\hfill \\ 4\mathit{\pi }\left({\tilde{\mathit{R}}}^2-{\tilde{\mathit{r}}}^2\right)/3\sqrt{3}{\mathit{a}}^2\mathrm{for}\mathrm{hollow}\sec \mathrm{tion}\hfill \end{array}\right.,$

(30) $\mathit{f}=\left\{\begin{array}{c}\mathit{\pi }\sqrt{3}{\tilde{\mathit{r}}}^2/2{\mathit{a}}^2\mathrm{for}\mathrm{circular}\sec \mathrm{tion}\hfill \\ 2\sqrt{3}{\mathit{l}}^2/{\mathit{a}}^2\mathrm{for}\mathrm{square}\sec \mathrm{tion}\hfill \\ 2\sqrt{3}{\mathit{l}}^2/{\mathit{a}}^2\mathrm{for}\mathrm{rotated}\mathrm{square}\sec \mathrm{tion}\hfill \\ \mathit{\pi }\sqrt{3}\left({\tilde{\mathit{R}}}^2-{\tilde{\mathit{r}}}^2\right)/2{\mathit{a}}^2\mathrm{for}\mathrm{hollow}\sec \mathrm{tion}\hfill \end{array}\right.,$

respectively.

The structure function for square, rectangular and triangular lattices F(g) is given by:

(31) $\mathit{F}\left(\stackrel{─}{\mathbf{g}}\right)=\frac{1}{{\mathit{S}}_{\mathit{c}}}\iint {\mathit{e}}^{-\mathit{j}\stackrel{─}{\mathbf{g}}.\mathbf{r}}{\mathit{d}}^2\mathit{r}.$

The integral in Eq. (31) is performed over the unit cell cross section. The structure functions for square, rectangular and triangular lattices are calculated by:

(32) $\mathit{F}\left(\stackrel{─}{\mathbf{g}}\right)=\left\{\begin{array}{c}2{\mathit{fJ}}_1\left(\stackrel{─}{\mathit{g}}\tilde{\mathit{r}}\right)/\stackrel{─}{\mathit{g}}\tilde{\mathit{r}}\mathrm{for}\mathrm{circular}\sec \mathrm{tion}\hfill \\ \mathit{f}\left[\sin \left({\stackrel{─}{\mathit{g}}}_1\mathit{l}\right)/{\stackrel{─}{\mathit{g}}}_1\mathit{l}\right]\left[\sin \left({\stackrel{─}{\mathit{g}}}_2\mathit{l}\right)/{\stackrel{─}{\mathit{g}}}_2\mathit{l}\right]\mathrm{for}\mathrm{square}\sec \mathrm{tion}\hfill \\ \mathit{f}\left\{\frac{\sin \left[\left(\mathit{l}/\sqrt{2}\right)\left({\stackrel{─}{\mathit{g}}}_1+{\stackrel{─}{\mathit{g}}}_2\right)\right]}{\left[\left(\mathit{l}/\sqrt{2}\right)\left({\stackrel{─}{\mathit{g}}}_1+{\stackrel{─}{\mathit{g}}}_2\right)\right]}\right\}\left\{\frac{\sin \left[\left(\mathit{l}/\sqrt{2}\right)\left(-{\stackrel{─}{\mathit{g}}}_1+{\stackrel{─}{\mathit{g}}}_2\right)\right]}{\left[\left(\mathit{l}/\sqrt{2}\right)\left(-{\stackrel{─}{\mathit{g}}}_1+{\stackrel{─}{\mathit{g}}}_2\right)\right]}\right\}\mathrm{for}\mathrm{square}\mathrm{rotated}\sec \mathrm{tion}\hfill \\ 2\mathit{f}\left[{\mathit{J}}_1\left(\stackrel{─}{\mathit{g}}\tilde{\mathit{R}}\right)-\left(\tilde{\mathit{r}}/\tilde{\mathit{R}}\right){\mathit{J}}_1\left(\stackrel{─}{\mathit{g}}\tilde{\mathit{R}}\right)\right]/\left(\stackrel{─}{\mathit{g}}\tilde{\mathit{R}}\right)\mathrm{for}\mathrm{hollow}\sec \mathrm{tion}\hfill \end{array}\right.$

where g = ‖g‖ and g_1,2 = ‖g_1,2‖.

The structure function for honeycomb lattice F_H(g) is defined as⁸⁴:

(33) $\begin{array}{c}{\mathit{F}}_{\mathit{H}}\left(\stackrel{─}{\mathbf{g}}\right)=\cos \left(\stackrel{─}{\mathbf{g}}·{\stackrel{─}{\mathbf{u}}}_1\right)\frac{1}{{\mathit{S}}_{\mathit{c}}}\iint {\mathit{e}}^{-\mathit{j}\stackrel{─}{\mathbf{g}}·\mathbf{r}}{\mathit{d}}^2\mathit{r}=\hfill \\ \cos \left(\stackrel{─}{\mathbf{g}}·{\stackrel{─}{\mathbf{u}}}_1\right)\mathit{F}\left(\stackrel{─}{\mathbf{g}}\right),\hfill \end{array}$

where ū₁ = -ū₂ = a(0,1/2) are vectors that define central position of the two BaTiO₃ inclusions within the honeycomb unit cell. We choose these vectors similar to Gao et al.⁸⁵. Thus, the structure functions of BaTiO₃ inclusions for honeycomb lattice are the same of Eq. (31) multiplied by cos(g ∙ u₁), considering ƒ from Eq. (29).

The structure function for Kagomé lattice F_Kg (g) is defined as:

(34) $\begin{array}{c}{\mathit{F}}_{\mathit{Kg}}\left(\stackrel{─}{\mathbf{g}}\right)=\frac{1}{3}\sum _{\mathit{i}=1}^3{\mathit{e}}^{\mathit{j}\left(\stackrel{─}{\mathbf{g}}·{\stackrel{─}{\mathbf{u}}}_i\right)}\frac{1}{{\mathit{S}}_{\mathit{c}}}\iint {\mathit{e}}^{-\mathit{j}\stackrel{─}{\mathbf{g}}·\mathbf{r}}{\mathit{d}}^2\mathit{r}=\hfill \\ \frac{1}{3}\sum _{\mathit{i}=1}^3{\mathit{e}}^{-\mathit{j}\left(\stackrel{─}{\mathbf{g}}·{\stackrel{─}{\mathbf{u}}}_i\right)}\mathit{F}\left(\stackrel{─}{\mathbf{g}}\right)\hfill \end{array}$

where ${\overline{\mathbf{u}}}_1=-\frac{\mathit{a}}{2}{\mathbf{e}}_1-\frac{\sqrt{3\mathit{a}}}{6}{\mathbf{e}}_2,{\overline{\mathbf{u}}}_2=\frac{\mathit{a}}{2}{\mathbf{e}}_1-\frac{\sqrt{3\mathit{a}}}{6}{\mathbf{e}}_2,{\overline{\mathbf{u}}}_3=\frac{\sqrt{3\mathit{a}}}{3}{\mathbf{e}}_1$ are vectors that define the three BaTiO₃ inclusion positions within the Kagomé unit cell. It is possible to split Eq. (34) in its real and imaginary parts:

(35) ${\mathit{F}}_{\mathit{kg}}\left(\stackrel{─}{\mathbf{g}}\right)=\mathit{Re}\left[{\mathit{F}}_{\mathit{kg}}\left(\stackrel{─}{\mathbf{g}}\right)\right]+\mathit{j}\mathrm{Im}\left[{\mathit{F}}_{\mathit{kg}}\left(\stackrel{─}{g}\right)\right]$

(36) $\mathit{Re}\left[{\mathit{F}}_{\mathit{kg}}\left(\stackrel{─}{\mathbf{g}}\right)\right]=\left[2\cos \left({\stackrel{─}{g}}_1{\stackrel{─}{\mathit{u}}}_{1_1}\right)\cos \left({\stackrel{─}{g}}_2{\stackrel{─}{\mathit{u}}}_{1_2}\right)+\cos \left(2{\stackrel{─}{g}}_2{\stackrel{─}{\mathit{u}}}_{1_2}\right)\right]\mathit{F}\left(\stackrel{─}{\mathbf{g}}\right),$

(37) $\mathrm{Im}\left[{\mathit{F}}_{\mathit{kg}}\left(\stackrel{─}{\mathbf{g}}\right)\right]=\left[-2\cos \left({\stackrel{─}{g}}_1{\stackrel{─}{\mathit{u}}}_{1_1}\right)\sin \left({\stackrel{─}{g}}_2{\stackrel{─}{\mathit{u}}}_{1_2}\right)+\sin \left(2{\stackrel{─}{g}}_1{\stackrel{─}{\mathit{u}}}_{1_2}\right)\right]\mathit{F}\left(\stackrel{─}{\mathbf{g}}\right),$

where ū_j1,2 = ‖ū_j1,2‖.

3. Results and Discussion

The physical parameters of BaTiO₃ inclusion (A) and polymeric matrix (B) are listed in Table 1. We calculate the elastic band structure considering initially a fixed filling fraction, 0.335, for the four inclusion cross section geometries considered - circular, hollow circular, square and rotated square with a 45º angle of rotation with respect to the x, y axes, in a square, rectangular, triangular, honeycomb and Kagomé lattices. In the course of numerical calculations, we consider 441 plane waves for Fourier series expansion (i.e. M = 10), which results in a good convergence. Band structure plots are limited until a maximum frequency of 0.12 THz.

Table 1
Physical parameters of BaTiO₃ inclusion (A) and polymeric matrix (B).

Geometry/Property	Value
Lattice parameter (ɑ)	30 x 10-9 m
Lattice parameter (ɑ1)	30 x 10-9 m
Lattice parameter (ɑ2)	25 x 10-9 m
Filling fraction (ƒ)	0.335
Mass density (ρA,ρB)	5730 kg/m3, 1150 kg/m3
Elastic constant (c11A,c11B)	166 x 109 N/m2, 7.8 x 109 N/m2
Elastic constant (c12A,c12B)	77 x 109 N/m2, 4.7 x 109 N/m2
Elastic constant (c44A,c44B)	43 x 109 N/m2, 1.6 x 109 N/m2
Elastic constant (c66A,c66B)	44.5 x 109 N/m2, 1.55 x 109 N/m2
Piezoelectric coefficient(e15A,e15B)	11.6 C/m2, 0 C/m2
Dielectric coefficient (𝜖11A,𝜖11B)	11.2 x 10-9 C2/Nm2, 0.0398 x 10-9 C2/Nm2

• Note that ais the lattice parameter for square, triangular, honeycomb and Kagomé lattices, and a₁, a₂are the lattice parameters for rectangular lattice.

Figure 2 (a-d) shows the band structures of the nano-piezoelectric PC for square lattice, considering the four types of inclusions and XY (red) and Z (blue) modes. We plot the band structures in the principal symmetry directions of FIBZ (see Figure 1 (f-h)). Plots are given in terms of frequency in Hz versus the reduced Bloch wave vector, k = kɑ/2π. The reduced Bloch wave vector for rectangular lattice is calculated taking into account $\mathit{a}=\sqrt{{\mathit{a}}_1^2+{\mathit{a}}_2^2}$. In Figure 2 (a), two complete band gaps are obtained for circular inclusion.

Figure 2
Elastic band structures of XY (red) and Z (blue) modes of BaTiO₃ inclusions in a polymeric matrix for square lattice. The following types of inclusions are considered - circular (a), hollow circular (b), square (c) and rotated square with a 45º angle of rotation with respect to the x, y axes.

From now on the relation between parameters R̃ and r̃ for hollow circular inclusion is r̃ = 0.2R̃ and we do not investigate the influence of BaTiO₃ thickness, i.e. R̃ - r̃, in the band structure. The thickness influence of a carbon microstructure inclusion for square lattice was investigated by Anjos & Arantes¹⁰. Figure 2 (b) presents one complete band gap for hollow circular inclusion and we can observe that first branches occur in higher frequencies compared to the other inclusions.

Figure 2 (c) shows two complete band gaps for square inclusion and its first complete band gap is broader than the complete band gaps for circular and hollow circular inclusions. When these square inclusions are rotated 45º with respect to x and y axes (Figure 2 (d)) three complete gaps are opened up. The broadest complete band gap for square lattice is for rotated square inclusion with a bandwidth of 11.97 GHz (43.5 GHz is the center frequency).

Figure 3 (a-d) shows the band structures for rectangular lattice. For circular and hollow circular inclusions, the band gap behavior of rectangular lattice is similar to square lattice, however they are narrower than the complete band gaps in a square lattice. For square and rotated square inclusions, only one and two complete band gaps are opened up, respectively. Similar to square lattice, the broadest complete band gap for rectangular lattice is opened up for rotated square inclusion with a bandwidth of 8.51 GHz (44.59 GHz is the center frequency).

Figure 3
Elastic band structures of XY (red) and Z (blue) modes of BaTiO₃ inclusions in a polymeric matrix for rectangular lattice. The following types of inclusions are considered - circular (a), hollow circular (b), square (c) and rotated square with a 45º angle of rotation with respect to the x, y axes.

Figure 4 (a-d) illustrates the band structures for triangular lattice. For circular, square and rotated square inclusions, a wide complete band gap is created. The broadest complete band gap for triangular lattice is for circular inclusion with a bandwidth of 12.98 GHz (36.31 GHz is the center frequency). A narrow complete band gap is also opened up for circular inclusion in high bands. For hollow inclusion, one complete band gap is created between 58.69 GHz-49.94 GHz. Comparing the band structures of triangular lattice to square and rectangular lattices, we observe that triangular lattice with circular inclusion presents the broadest complete band gap.

Figure 4
Elastic band structures of XY (red) and Z (blue) modes of BaTiO₃ inclusions in a polymeric matrix for triangular lattice. The following types of inclusions are considered - circular (a), hollow circular (b), square (c) and rotated square with a 45º angle of rotation with respect to the x, y axes.

In Figures 5 (a-d) and 6 (a-d), honeycomb and Kagomé lattices, respectively, the band structures of square and rotated square inclusions present similar band gap behavior. The band structure of circular inclusion in honeycomb and Kagomé lattices presents seven and five complete band gaps, respectively. For Kagomé lattice, all the complete band gaps are narrower than the broadest complete band gap from square, rectangular, triangular and honeycomb lattices. The broadest complete band gap for honeycomb lattice is for circular inclusion with a bandwidth of 9.77 GHz (31.19 GHz is the center frequency). For hollow inclusion in honeycomb and Kagomé lattices, we can observe two and one complete band gaps, respectively, and the first bands appear in higher frequencies compared to the other inclusions, similar as observed to other lattices - square, rectangular and triangular.

Figure 5
Elastic band structures of XY (red) and Z (blue) modes of BaTiO₃ inclusions in a polymeric matrix for honeycomb lattice. The following types of inclusions are considered - circular (a), hollow circular (b), square (c) and rotated square with a 45º angle of rotation with respect to the x, y axes.

Figure 6
Elastic band structures of XY (red) and Z (blue) modes of BaTiO₃ inclusions in a polymeric matrix for Kagomé lattice. The following types of inclusions are considered - circular (a), hollow circular (b), square (c) and rotated square with a 45º angle of rotation with respect to the x, y axes.

Comparing the band structures for honeycomb and Kagomé lattices, we observe that honeycomb lattice presents the best performance (more complete band gaps and broader complete band gaps) for all inclusions.

In Figures 7-11, we show the comparison between Z modes of the nano-piezoelectric PC with (blue asterisks) and without (e_15A = 0) (black circles) piezoelectricity for square, rectangular, triangular, honeycomb and Kagomé lattices, respectively. The influence of piezoelectricity is significant on the band gaps for all lattices and inclusions. This influence is more evident in high frequencies⁶⁰ for circular, square and rotated square inclusions. However, for hollow circular inclusion in all lattices (Figures 7-11 (b)), the piezoelectricity is also significant in lower frequencies.

Figure 7
Elastic band structures of BaTiO₃ inclusions in a polymeric matrix considering Z mode with (blue asterisks) and without (black circles) piezoelectricity for square lattice. The following types of inclusions are considered - circular (a), hollow circular (b), square (c) and rotated square with a 45º angle of rotation with respect to the x, y axes.

Figure 8
Elastic band structures of BaTiO₃ inclusions in a polymeric matrix considering Z mode with (blue asterisks) and without (black circles) piezoelectricity for rectangular lattice. The following types of inclusions are considered - circular (a), hollow circular (b), square (c) and rotated square with a 45º angle of rotation with respect to the x, y axes.

Figure 9
Elastic band structures of BaTiO₃ inclusions in a polymeric matrix considering Z mode with (blue asterisks) and without (black circles) piezoelectricity for triangular lattice. The following types of inclusions are considered - circular (a), hollow circular (b), square (c) and rotated square with a 45º angle of rotation with respect to the x, y axes.

Figure 10
Elastic band structures of BaTiO₃ inclusions in a polymeric matrix considering Z mode with (blue asterisks) and without (black circles) piezoelectricity for honeycomb lattice. The following types of inclusions are considered - circular (a), hollow circular (b), square (c) and rotated square with a 45º angle of rotation with respect to the x, y axes.

Figure 11
Elastic band structures of BaTiO₃ inclusions in a polymeric matrix considering Z mode with (blue asterisks) and without (black circles) piezoelectricity for Kagomé lattice. The following types of inclusions are considered - circular (a), hollow circular (b), square (c) and rotated square with a 45º angle of rotation with respect to the x, y axes.

We illustrate in Figures 12-16 the complete band gap widths between XY and Z modes as a function of filling fraction for all inclusions and lattices. Piezoelectricity is included in this analysis.

Figure 12
Elastic complete band gap widths between XY and Z modes of BaTiO₃ inclusions in a polymeric matrix as a function of filling fraction for square lattice. The following types of inclusions are considered - circular (a), hollow circular (b), square (c) and rotated square with a 45º angle of rotation with respect to the x, y axes.

In a square lattice (Figure 12), circular and square inclusions provide the best behavior with high complete band gap widths in a broad range of filling fraction. Considering the first band gap, square inclusion shows the highest one in a broad range of filling fraction. Furthermore, square inclusion also presents four complete band gaps. Hollow circular inclusion presents only one complete band gap, whereas rotated square inclusion does not present better performance than square inclusion.

For rectangular lattice, Figure 13, hollow circular inclusion opens up just one narrow complete band gap, whereas rotated square inclusion creates three complete band gaps, however, only the second one appear in a broad range of filling fraction. Similar as square lattice (Figure 12), circular and square inclusions present the broadest complete band gaps for rectangular lattice. Complete band gap widths in rectangular lattice are lower than in square lattice for all inclusions.

Figure 13
Elastic complete band gap widths between XY and Z modes of BaTiO₃ inclusions in a polymeric matrix as a function of filling fraction for rectangular lattice. The following types of inclusions are considered - circular (a), hollow circular (b), square (c) and rotated square with a 45º angle of rotation with respect to the x, y axes.

Figure 14 shows complete band gap widths for triangular lattice. Circular inclusion presents the best performance. Comparing triangular lattice (Figure 14) to square lattice (Figure 12), triangular lattice presents the best performance for circular (first complete band gap), hollow circular and rotated square inclusions. However, for square inclusion, square lattice shows better performance than triangular lattice. There are two important variables to conclude which is the best band gap performance, that is to say the band gap width (the higher, the better) and the filling fraction range in which the band gaps are opened (the broader, the better).

Figure 14
Elastic complete band gap widths between XY and Z modes of BaTiO₃ inclusions in a polymeric matrix as a function of filling fraction for triangular lattice. The following types of inclusions are considered - circular (a), hollow circular (b), square (c) and rotated square with a 45º angle of rotation with respect to the x, y axes.

Comparing honeycomb (Figure 15) and Kagomé (Figure 16) lattices, the best performance is found for honeycomb lattice considering circular, hollow circular and rotated square inclusions. However, the difference among these lattices is subtle. The highest band gap width for honeycomb and Kagomé lattices is for circular inclusion. Comparing honeycomb and Kagomé lattices to the other lattices, we can observe that for all inclusions, the best behavior is found for square (Figure 12) and triangular (Figure 14) lattices. The worst performance is found for rectangular lattice with hollow circular inclusion (Figure 13 (b)).

Figure 15
Elastic complete band gap widths between XY and Z modes of BaTiO₃ inclusions in a polymeric matrix as a function of filling fraction for honeycomb lattice. The following types of inclusions are considered - circular (a), hollow circular (b), square (c) and rotated square with a 45º angle of rotation with respect to the x, y axes.

Figure 16
Elastic complete band gap widths between XY and Z modes of BaTiO₃ inclusions in a polymeric matrix as a function of filling fraction for Kagomé lattice. The following types of inclusions are considered - circular (a), hollow circular (b), square (c) and rotated square with a 45º angle of rotation with respect to the x, y axes.

4. Conclusion

Broad complete band gaps are obtained for a nano-piezoelectric PC, consisting of BaTiO₃ inclusions embedded in a polymeric matrix, for different inclusion geometries and different lattices. We firstly study the elastic band structure for a fixed filling fraction of 0.335 considering XY and Z modes. For square lattice, the broadest complete band gap is found for rotated square inclusion with 11.97 GHz of bandwidth. This inclusion geometry is also the best one for rectangular lattice. Triangular lattice presents the broadest complete band gap than the other lattices for circular inclusion. The elastic band structure for honeycomb and Kagomé lattices presents several complete band gaps. However, for Kagomé lattice, all of them are narrower than the complete band gaps from square, rectangular, triangular and honeycomb lattices.

The Z mode with and without piezoelectric effect is also investigated for a fixed filling fraction of 0.335. This effect is significant on the band gaps for all lattices and inclusions, and it is more evident in higher frequencies for circular, square and rotated square inclusions. However, piezoelectricity is significant in lower frequencies for hollow circular inclusion.

We also analyze the complete band gap bandwidths as a function of filling fraction for all inclusions and lattices with piezoelectric effect. Among square, rectangular and triangular lattices, triangular lattice presents the best performance for circular inclusion (first complete band gap), hollow circular and rotated square inclusions. For square inclusion, square lattice shows the best performance. Among honeycomb and Kagomé lattices, the best performance is found for honeycomb lattice considering circular, hollow circular and rotated square inclusions.

The best band gap performance (i.e. highest band gap width and band gap opened up in a broad range of filling fraction) of the nano-piezoelectric PC investigated is for circular inclusion in a triangular lattice. However, it also depends on the application, for instance the frequency range of interest. Finally, we suggest that elastic complete band gaps in nano-piezoelectric PCs enlarge the potential applications for elastic vibration control in GHz.

5. Acknowledgments

The authors gratefully acknowledge the financial support of this investigation by the Brazilian research funding agency FAPEMA (State Funding Agency of Maranhão) and by IFMA (Federal Institute of Education, Science and Technology of Maranhão).

6. References

• ¹
  Sigalas M, Economou EN. Band structure of elastic waves in two dimensional systems. Solid State Communications 1993;86(3):141-143.
• ²
  Kushwaha MS, Halevi P, Martínez G, Dobrzynski L, Djafari-Rouhani B. Theory of acoustic band structure of periodic elastic composites. Physical Review B 1994;49(4):2313-2322.
• ³
  Huang GL, Sun CT. Modeling heterostructures of nanophononic crystals by continuum model with microstructures. Applied Physics Letters 2006;88(26):261908.
• ⁴
  Hepplestone SP, Srivastava GP. Hypersonic modes in nanophononic semiconductors. Physical Review Letters 2008;101(10):105502.
• ⁵
  Schneider D, Liaqat F, El Boudouti EH, El Hassouani Y, Djafari-Rouhani B, Tremel W, et al. Engineering the hypersonic phononic band gap of hybrid Bragg stacks. Nano Letters 2012;12(6):3101-3108.
• ⁶
  Graczykowski B, Sledzinska M, Kehagias N, Alzina F, Reparaz JS, Sotomayor Torres CM. Hypersonic phonon propagation in one-dimensional surface phononic crystal. Applied Physics Letters 2014;104(12):123108.
• ⁷
  Pennec Y, Vasseur JO, Djafari-Rouhani B, Dobrzynski L, Deymier PA. Two-dimensional phononic crystals: Examples and applications. Surface Science Reports 2010;65(8):229-291.
• ⁸
  Huang J, Shi Z. Attenuation zones of periodic pile barriers and its application in vibration reduction for plane waves. Journal of Sound and Vibration 2013;332(19):4423-4439.
• ⁹
  Yu K, Chen T, Wang X. Band gaps in the low-frequency range based on the two-dimensional phononic crystal plates composed of rubber matrix with periodic steel stubs. Physica B: Condensed Matter 2013;416:12-16.
• ¹⁰
  Anjos V, Arantes A. Phononic band structure in carbon microtube composites. RSC Advances 2015;5(15):11248-11253.
• ¹¹
  Miranda EJP Jr, Dos Santos JMC. Flexural wave band gaps in metamaterial elastic beam. In: Proceedings of the 23^rd ABCM International Congress of Mechanical Engineering (COBEM 2015); 2015 Dec 6-11; Rio de Janeiro, RJ, Brazil. 8 p.
• ¹²
  Miranda EJP Jr, Dos Santos JMC. Phononic band gaps in Al₂O₃/epoxy composite. Materials Science Forum 2018;912:112-117.
• ¹³
  Miranda EJP Jr, Dos Santos JMC. Flexural wave band gaps in Al₂O₃/epoxy composite rectangular plate using Mindlin theory. In: Proceedings of the 3^rd Brazilian Conference on Composite Materials (BCCM-3); 2016 Aug 28-31; Gramado, RS, Brazil. 8 p.
• ¹⁴
  Miranda EJP Jr, Dos Santos JMC. Elastic wave band gaps in a two-dimensional magnetoelectroelastic phononic crystal. Revista Interdisciplinar de Pesquisa em Engenharia 2016;2(13):13-26.
• ¹⁵
  Yablonovitch E. Inhibited spontaneous emission in solid-state physics and electronics. Physical Review Letters 1987;58(20):2059-2062.
• ¹⁶
  John S. Strong localization of photons in certain disordered dielectric superlattices. Physical Review Letters 1987;58(23):2486-2489.
• ¹⁷
  Olsson III RH, El-Kady I. Microfabricated phononic crystal devices and applications. Measurement Science and Technology 2009;20:012002.
• ¹⁸
  Jensen JS. Phononic band gaps and vibrations in one- and two-dimensional mass-spring structures. Journal of Sound and Vibration 2003;266(5):1053-1078.
• ¹⁹
  Wang G, Wen JH, Wen XS. Quasi-one-dimensional phononic crystals studied using the improved lumped-mass method: Application to locally resonant beams with flexural wave band gap. Physical Review B 2005;71(10):104302.
• ²⁰
  Casadei F, Beck BS, Cunefare KA, Ruzzene M. Vibration control of plates through hybrid configurations of periodic piezoeletric shunts. Journal of Intelligent Material Systems and Structures 2012;23(10):1169-1177.
• ²¹
  Miranda EJP Jr, Dos Santos JMC. Flexural wave band gaps in elastic metamaterial thin plate. In: Proceedings of the IX Mechanical Engineering Brazilian Congress (CONEM 2016); 2016 Aug 21-25; Fortaleza, CE, Brazil; 10 p.
• ²²
  Miranda EJP Jr, Dos Santos JMC. Flexural wave band gaps in elastic metamaterial beam. In: Proceedings of the International Conference on Noise and Vibration Engineering (ISMA2016). 2016 Sep 19-23; Leuven, Belgium. p. 2099-2113.
• ²³
  Ho KM, Cheng CK, Yang Z, Zhang XX, Sheng P. Broadband locally resonant sonic shields. Applied Physics Letters 2003;83(26):5566.
• ²⁴
  Qiu C, Liu Z, Mei J, Shi J. Mode-selecting acoustic filter by using resonant tunneling of two-dimensional double phononic crystals. Applied Physics Letters 2005;87(10):104101.
• ²⁵
  Yang Z, Dai HM, Chan NH, Ma GC, Sheng P. Acoustic metamaterial panels for sound attenuation in the 50-1000 Hz regime. Applied Physics Letters 2010;96(4):041906.
• ²⁶
  Casadei F, Dozio L, Ruzzene M, Cunefare KA. Periodic shunts arrays for the control of noise radiation in an enclosure. Journal of Sound and Vibration 2010;329(18):3632-3646.
• ²⁷
  Xiao Y, Wen J, Wen X. Sound transmission loss of metamaterial-based thin plates with multiple subwavelength arrays of attached resonators. Journal of Sound and Vibration 2012;331(25):5408-5423.
• ²⁸
  Benchabane S, Khelif A, Robert L, Rauch JY, Pastureaud T, Laude V. Elastic band gaps for surface modes in an ultrasonic lithium niobate phononic crystal. In: Proceedings of SPIE Photonics Europe; 2006; Strasbourg, France. 2006. 618216.
• ²⁹
  Comerio MC. Can buildings be made earthquake-safe? Science 2006;312(5771):204-206.
• ³⁰
  Gorishnyy T, Ullal CK, Maldovan M, Fytas G, Thomas EL. Hypersonic phononic crystals. Physical Review Letters 2005;94(11):115501.
• ³¹
  Cheng W, Sainidou R, Burgardt P, Stefanou N, Kiyanova A, Efremov M, et al. Elastic properties and glass transition of supported polymer thin films. Macromolecules 2007;40(20):7283-7290.
• ³²
  Liu Z, Zhang X, Mao Y, Zhu YY, Yang Z, Chan CT, et al. Locally resonant sonic materials. Science 2000;289(5485):1734-1736.
• ³³
  Davis BL, Hussein M. Nanophononic metamaterial: Thermal conductivity reduction by local resonance. Physical Review Letters 2014;112(5-7):055505.
• ³⁴
  Wagner MR, Graczykowski B, Reparaz JS, Sachat AE, Sledzinska M, Alzina F, et al. Two-dimensional phononic crystals: Disorder matters. Nano Letters 2016;16(9):5661-5668.
• ³⁵
  Yang L, Yang N, Li B. Extreme low thermal conductivity in nanoscale 3D Si phononic crystal with spherical pores. Nano Letters 2014;14(4):1734-1738.
• ³⁶
  Zen N, Puurtinen TA, Isotalo TJ, Chaudhuri S, Maasilta IJ. Engineering thermal conductance using a two-dimensional phononic crystal. Nature Communications 2014;5:3435.
• ³⁷
  Anufriev R, Nomura M. Thermal conductance boost in phononic crystal nanostructures. Physical Review B 2015;91:245417.
• ³⁸
  Lacatena V, Haras M, Robillard JF, Monfray S, Skotnicki T, Dubois E. Toward quantitative modeling of silicon phononic thermocrystals. Applied Physics Letters 2015;106(11):114104.
• ³⁹
  Honarvar H, Yang L, Hussein MI. Thermal transport size effects in silicon membranes featuring nanopillars as local resonators. Applied Physics Letters 2016; 108(263101):1-5.
• ⁴⁰
  Wilm M, Ballandras S, Laude V, Pastureaud T, Pierre G, Steichen W. A plane-wave-expansion approach for modelling acoustic propagation in 2D and 3D piezoelectric periodic structures. In: IEEE Ultrasonics Symposium; 2001 Oct 7-10; Atlanta, GA, USA. p. 977-980.
• ⁴¹
  Wilm M, Ballandras S, Laude V, Pastureaud T. A full 3D plane-wave-expansion model for 1-3 piezoelectric composite structures. Journal of the Acoustical Society of America 2002;112(3 Pt 1):943-952.
• ⁴²
  Wilm M, Khelif A, Ballandras S, Laude V, Djafari-Rouhani B. Out-of-plane propagation of elastic waves in two-dimensional phononic band-gap materials. Physical Review E 2003;67(6):065602.
• ⁴³
  Hou Z, Wu F, Liu Y. Phononic crystals containing piezoelectric material. Solid State Communications 2004;130(11):745-749.
• ⁴⁴
  Wu TT, Hsu ZC, Huang ZG. Band gaps and the electromechanical coupling coefficient of a surface acoustic wave in a two-dimensional piezoelectric phononic crystal. Physical Review B 2005;71(6):064303.
• ⁴⁵
  Hsu JC, Wu TT. Bleustein-Gulyaev-Shimizu acoustic waves in two-dimensional piezoelectric phononic crystals. IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control 2006;53(6):1169-1176.
• ⁴⁶
  Wang YZ, Li FM, Huang WH, Wang YS. Effects of inclusion shapes on the band gaps in two-dimensional piezoelectric phononic crystals. Journal of Physics: Condensed Matter 2007;19(49):496204.
• ⁴⁷
  Hsu JC, Wu TT. Calculations of Lamb wave band gaps and dispersions for piezoelectric phononic plates using Mindlin's theory-based plane wave expansion method. IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control 2008;55(2):431-441.
• ⁴⁸
  Qian ZH, Jin F, Li FM, Kishimoto K. Complete band gaps in two-dimensional piezoelectric phononic crystals with {1-3} connectivity family. International Journal of Solids and Structures 2008;45(17):4748-4755.
• ⁴⁹
  Zou XY, Chen Q, Liang B, Cheng JC. Control of the elastic wave bandgaps in two-dimensional piezoelectric periodic structures. Smart Materials and Structures 2008;17(1):015008.
• ⁵⁰
  Wang YZ, Li FM, Kishimoto K, Wang YW, Huang WH. Wave band gaps in three-dimensional periodic piezoelectric structures. Mechanics Research Communications 2009;36(4):461-468.
• ⁵¹
  Wang Y, Li F, Wang Y, Kishimoto K, Huang W. Tuning of band gaps for a two-dimensional piezoelectric phononic crystal with a rectangular lattice. Acta Mechanica Sinica 2009;25(1):65-71.
• ⁵²
  Huang Y, Zhang CL, Chen WQ. Elastic wave band structures and defect states in a periodically corrugated piezoelectric plate. Journal of Applied Mechanics 2014;81(8):081005.
• ⁵³
  Lian Z, Jiang S, Hu H, Dai L, Chen X, Jiang W. An enhanced plane wave expansion method to solve piezoelectric phononic crystal with resonant shunting circuits. Shock and Vibration 2016;2016:4015363.
• ⁵⁴
  Li L, Guo Y. Analysis of longitudinal waves in rod-type piezoelectric phononic crystals. Crystals 2016;6(4):45.
• ⁵⁵
  Robillard JF, Bou Matar O, Vasseur JO, Deymier PA, Stippinger M, Hladky-Hennion AC, et al. Tunable magnetoelastic phononic crystals. Applied Physics Letters 2016;95(12):124104.
• ⁵⁶
  Vasseur JO, Bou Matar O, Robillard JF, Hladky-Hennion AC, Deymier PA. Band structures tunability of bulk 2D phononic crystals made of magneto-elastic materials. AIP Advances 2011;1(4):041904.
• ⁵⁷
  Bou Matar O, Robillard JF, Vasseur JO, Hladky-Hennion AC, Deymier PA, Preobrazhensky V. Band gap tunability of magneto-elastic phononic crystal. Journal of Applied Physics 2012;111(5):054901.
• ⁵⁸
  Zhou C, Sai Y, Chen J. Tunable Lamb wave band gaps in two-dimensional magnetoelastic phononic crystal slabs by an applied external magnetostatic field. Ultrasonics 2016;71:69-74.
• ⁵⁹
  Wang YZ, Li FM, Huang WH, Jiang X, Wang YS, Kishimoto K. Wave band gaps in two-dimensional piezoelectric/piezomagnetic phononic crystals. International Journal of Solids and Structures 2008;48(14-15):4203-4210.
• ⁶⁰
  Wang YZ, Li FM, Kishimoto K, Wang YS, Huang WH. Elastic wave band gaps in magnetoelectroelastic phononic crystals. Wave Motion 2009;46(1):47-56.
• ⁶¹
  Wang YZ, Li FM, Kishimoto K, Wang YS, Huang WH. Band gap behaviours of periodic magnetoelectroelastic composite structures with kagome lattices. Waves in Random and Complex Media 2009;19(3):509-520.
• ⁶²
  Lan M, Wei P. Band gap of piezoelectric/piezomagnetic phononic crystal with graded interlayer. Acta Mechanica 2014;225(6):1779-1794.
• ⁶³
  Guo X, Wei P, Lan M, Li L. Dispersion relations of elastic waves in one-dimensional piezoelectric/piezomagnetic phononic crystal with functionally graded interlayers. Ultrasonics 2016;70:158-171.
• ⁶⁴
  El-Kady I, Olsson III RH, Fleming JG. Phononic band-gap crystals for radio frequency communications. Applied Physics Letters 2008;92(23):233504.
• ⁶⁵
  Soliman YM, Su MF, Leseman ZC, Reinke CM, El-Kady I, Olsson III RH. Phononic crystals operating in the gigahertz range with extremely wide band gaps. Applied Physics Letters 2010;97(19):193502.
• ⁶⁶
  Robillard JF, Devos A, Roch-Jeune I. Time-resolved vibrations of two-dimensional hypersonic phononic crystals. Physical Review B 2007;76(9):092301.
• ⁶⁷
  Ezzahri Y, Grauby S, Rampnoux JM, Michel H, Pernot G, Claeys W, et al. Coherent phonons in Si/Si Ge superlattices. Physical Review B 2007;75(19):195309.
• ⁶⁸
  Olsson III RH, Ziaei-Moayyed M, Kim B, Reinke C, Su MF, Hopkins P, et al. Micro and nano fabricated phononic crystals: Technology and applications. In: Proceedings of the IEEE International Ultrasonics Symposium; 2011 Oct 18-21; Orlando, FL, USA. p. 983-988.
• ⁶⁹
  Graczykowski B, Sledzinska M, Alzina F, Gomis-Bresco J, Reparaz JS, Wagner MR, et al. Phonon dispersion in hypersonic two-dimensional phononic crystal membranes. Physical Review B 2015;91(7):075414.
• ⁷⁰
  Travagliati M, Nardi D, Giannetti C, Gusev V, Pingue P, Piazza V, et al. Interface nano-confined acoustic waves in polymeric surface phononic crystals. Applied Physics Letters 2015;106(2):021906.
• ⁷¹
  Sadat SM, Wang RY. Colloidal nanocrystal superlattices as phononic crystals: Plane wave expansion modeling of phonon band structure. RSC Advances 2016;6(50):44578-44587.
• ⁷²
  Sgouros AP, Neupane MR, Sigalas MM, Aravantinos-Zafiris N, Lake RK. Nanoscale phononic interconnects in THz frequencies. Physical Chemistry Chemical Physics 2014;16(42):23355-23364.
• ⁷³
  Zhen N, Wang YS, Zhang C. Surface/interface effect on band structures of nanosized phononic crystals. Mechanics Research Communications 2012;46:81-89.
• ⁷⁴
  Zhen N, Wang YS, Zhang C. Bandgap calculation of in-plane waves in nanoscale phononic crystals taking in account surface/interface effects. Physica E: Low-dimensional Systems and Nanostructures 2013;54:125-132.
• ⁷⁵
  Abdollahi A, Moravvej-Farshi MK. Designing a hypersonic filter on a nanophononic crystal platform. In: Proceedings of the 22^nd Iranian Conference on Electrical Engineering (ICEE 2014); 2014 May 20-22; Tehran, Iran. p. 26-28.
• ⁷⁶
  Wang ZL, Song J. Piezoelectric nanogenerators based on zinc oxide nanowire arrays. Science 2006;312(5771):242-246.
• ⁷⁷
  Wang ZL. The new field of nanopiezotronics. Materialstoday 2007;10(5):20-28.
• ⁷⁸
  Yang F, Li JCM. Micro and Nano Mechanical Testing of Materials and Devices New York: Springer; 2008.
• ⁷⁹
  Wang X. Piezoelectric nanogenerators - Harvesting ambient mechanical energy at the nanometer scale. Nano Energy 2012;1(1):13-24.
• ⁸⁰
  Brillouin L. Wave Propagation in Periodic Structures Mineola: Dover Publications; 1946.
• ⁸¹
  Kuang W, Hou Z, Liu Y. The effects of shapes and symmetries of scatterers on the phononic band gap in 2D phononic crystals. Physics Letters A 2004;332(5-6):481-490.
• ⁸²
  Floquet G. Sur les équations différentielles linéaires à coefficients périodiques. Annales Scientifiques de L'École Normale Supérieure 1883;12:47-88.
• ⁸³
  Bloch F. Über die quantenmechanik der electron in kristallgittern. Zeitschrift für Physik 1929;52(7-8):555-600.
• ⁸⁴
  Cassagne D, Jouanin C, Bertho D. Hexagonal photonic-band gap structures. Physical Review B 1996;53(11):7134-7142.
• ⁸⁵
  Gao Z, Fang J, Zhang Y, Jiang L. Band structure research of a 2D honeycomb lattice phononic crystal. International Journal of Electrochemical Science 2013;8:7918-7925.

Publication Dates

• Publication in this collection
  30 Oct 2017
• Date of issue
  2017

History

• Received
  17 Mar 2017
• Reviewed
  16 Aug 2017
• Accepted
  17 Sept 2017