Phononic crystals (PCs) are composite materials made of arrays of constituents embedded in host materials 1 2 3 . The interesting characteristics of the PCs are related to the existence of phononic bandgaps (i.e., frequency gaps) due to the Bragg reflections of the phonons with long wavelengths. We can regard the PC as an opaque barrier for the phonons within the phononic bandgaps 4 . This suggests the potential for designing various phonon optic devices, such as phonon filters, mirrors, resonators, etc.

Recent advances in fabrication methods enable realization of one-dimensional hetero-structures, i.e., nanowire superlattices (NWSLs) 5 6 7 8 9 10 . Their electronic and optical properties were studied, and a variety of possible applications utilizing the characteristics were also proposed 11 12 13 14 . In addition, the NWSLs are expected to yield interesting physical effects on phonons, which influence the electronic states and the transport properties via the electron-phonon interaction. These NWSLs can be regarded as wire-type phononic crystals (WPCs), in which the phononic bandgaps are induced by the periodicity along the wire axis.

In a previous paper 15 , we developed a numerical method to derive phonon modes in a free-standing NWSL of anisotropic material with an arbitrary shape of cross-section. As examples, the phonon modes were calculated for the rectangular and square cross-section GaAs/AlAs and InP/InAs NWSLs composed of anisotropic materials 15 16 .

Though above result revealed the important aspects of phonon modes in the NWSLs, it seems to be difficult to design WPCs with complete phononic bandgaps because in the dispersion relations of these NWSLs, many subbands are folded into the mini-Brillouin zone and the frequencies of gaps are different with phonon modes. In addition, the gap widths are narrow in these NWSLs because the difference of the acoustic impedance between the GaAs and AlAs layers (or between the InP and InAs layers) is small. The NWSLs with large acoustic mismatch would be suitable for designing the phonon optic devices.

In the present work, we numerically calculate the dispersion relations and corresponding displacement fields for a circular cross-section NWSL consisting of GaN and AlN, and we determine a set of parameters which gives complete frequency gaps.

The equation giving the eigenfrequencies of phonon modes in a freestanding NWSL composed of anisotropic crystals was formulated in 15 . In this method, the displacement components u_i (i = x, y, z) are expanded in terms of a set of basis functions ϕ _α  (r)

u i ( r ) = ∑ α A αi ϕ α ( r ) .

The expansion coefficients A _αi and the eigenfrequencies ω are determined by solving the generalized eigenvalue equation:

∑ α , ℓ H βi , αℓ A αℓ = ω 2 ∑ α , ℓ S βi , αℓ A αℓ .

For the NWSLs composed of cubic materials, the matrix elements H _βi,αℓ and S _βi,αℓ are written as

H β i , α ℓ = i β C 11 i α + ∑ j ≠ i j β C 44 j α δ i ℓ + i β C 12 ℓ α + ℓ β C 44 i α 1 − δ i ℓ ,

S β i , α ℓ = δ i ℓ β ρ α ,

where

j β C μ ν k α = ∫ V ∂ ϕ β ∗ ( r ) ∂ x j C μ ν ( r ) ∂ ϕ α ( r ) ∂ x k d r ,

β ρ α = ∫ V ϕ β ∗ ( r ) ρ ( r ) ϕ α ( r ) d r .

Here, C _μν and ρ are the stiffness tensor and mass density, respectively, which are dependent on r in the NWSLs.

As basis functions, we adopt the product of powers of the Cartesian coordinates in the xy plane. The z dependence is expressed in the form of the Bloch wave:

ϕ α ( k , r ) = 1 V x R m y R n e i k + G z .

Here, R denotes the radius of the wire; G is the reciprocal lattice vector; and V = Π R ² D is the volume of the unit cell, where D is the length of the unit cell in the z direction. The basis functions are specified with α = (m,n,G). The expressions for the matrix elements H _βi,αℓ and S _βi,αℓ can be analytically obtained for the circular NWSLs.

Based on group theory 17 18 , the phonon modes are classified symmetrically. In the present study, the constituent layers are assumed to be cubic materials, i.e., zinc-blende structure. The group of k is C _4v for 0 < k < Π/D. The irreducible representations are A₁, A₂, B₁, B₂, and E 15 .

By considering the above symmetry, the symmetry-adapted basis function in the present system can be constructed, and the phonon dispersion relations of each mode are independently calculated.

The boundary condition at the free surface of the wire requires that the wave numbers in the lateral direction are discretized. On the other hand, the wave number k in the longitudinal direction has a continuous value. Therefore, even the homogeneous plain nanowire has subband structure. In the NWSL, the subbands are folded into the mini-Brillouin zone. In other words, the size of the mini-Brillouin zone and corresponding phononic bandgaps are determined by the periodicity D of the NWSL.

For the nanowires with smaller R, the maximum wavelengths in the lateral directions become shorter. As a result, the subbands except for the lowest dispersion curve of each mode go up to the higher frequency region, though a lot of subbands are folded into the mini-Brillouin zone. Changing the ratio of R and D, we can control the phonon modes in the lower frequency range.

Figure 1a illustrates the phonon dispersion relations calculated for R = 5.0 nm and D = 20.0 nm (thicknesses of both GaN and AlN layers are 10.0 nm). Other parameters we used are as follows: ρ= 6.15 g/cm³, C₁₁ = 2.96, C₁₂ = 1.54, and C₄₄ = 2.06 (all in units of 10¹² dyn/cm²) for GaN; ρ = 3.26 g/cm³, C₁₁ = 3.04, C₁₂ = 1.52, and C₄₄ = 1.99 (all in units of 10¹² dyn/cm²) for AlN 19 .

In the present frequency range, we can see four different modes, i.e., the A₁, A₂, B₁, and E modes. The dispersion curves corresponding to the B₂ modes exist in higher frequency range. The lowest dispersion curve of the A₂ mode is mostly overlapped with that of the A₁ mode.

The complete frequency gaps are realized for the parameters we selected. For comparison, we show in Figure 1b the phonon dispersion relations calculated for R = 5.0 nm and D = 45.0 nm (thicknesses of GaN and AlN layers are 25.0 and 20.0 nm, respectively). In this example, the broader bandgaps of the A₁, A₂, and E modes are nearly coincident with each other. Here, we note that the complete bandgaps disappear if the isotropic approximation is used for each constituent layer.

The lowest dispersion curves of the A₁ and A₂ modes are linear in k, i.e., ω vanishes at k = 0. On the other hand, the lowest frequencies of the B₁ and B₂ modes have finite values. For thin (thick) NWSLs, the lowest frequencies of B₁ and B₂ modes become higher (lower).

Figure 2a,b shows the displacement patterns corresponding to the lowest A₁ and A₂ modes at k = Π / D, respectively. The lowest A₁ mode has the feature of a dilatational mode, while the lowest A₂ mode shows the feature of a torsional mode.

Figure 2c,d,e shows the displacement patterns corresponding to the B₁, B₂, and E modes at k = Π / D, respectively. The lowest B₁ mode shows the feature of a stretching mode, i.e., the alternating dilatation and contraction in the x and y directions, while the lowest B₂ mode shows the features of a shear mode, i.e., alternating stretching in the two diagonal directions.

All dispersion curves of the E mode are doubly degenerate because the irreducible representation of the E mode is two-dimensional. For the E modes, the lowest dispersion curve near k = 0 is proportional to k². This parabolic behavior is due to the fact that these modes correspond to the bending of the NWSL. Figure 2e clearly shows that the E modes have a feature of flexural mode.

We theoretically studied the acoustic phonon modes in circular cross-section NWSLs consisting of cubic GaN and AlN. We calculated their dispersion relations and phonon displacement fields. These modes are classified into five types, i.e., A₁, A₂, B₁, B₂, and E modes, which have features of dilatational, torsional, stretching, shear, and flexural modes, respectively. We determined a set of parameters which gives complete phononic bandgaps. The results suggest the realization of the optimized phonon devices, such as phonon filters or mirrors in the micro/nano electromechanical systems.

In the present work, we only showed the results for the NWSLs consisting of cubic materials (i.e., zinc-blende structure). The results for the NWSLs with wurtzite structure will be given and the difference will be discussed elsewhere.

The author declares that he has no competing interests.