1. Introduction
Acousto-optic devices are well known for their broad application in research and technology. To date, more than a dozen acousto-optic devices using different principles and serving various objectives have been introduced and explored [
1].
In designing new acousto-optic devices, particular emphasis is put on the choice of a material in which optical radiation interacts with an acoustic wave. The choice of material depends on the acoustic, optic, and acousto-optic properties of the medium [
2]. Most commonly, priority is given to the acoustic properties of the material because they play a vital role in defining the diffraction efficiency. It is known that the acousto-optic figure of merit of the material is inversely proportional to the cube of the acoustic wave phase velocity [
3]. In the acousto-optic field, it is therefore essential to find and analyze the materials characterized by the low phase velocity of propagating acoustic waves.
A crucial factor in choosing the material is the anisotropy of its acoustic properties. Strong acoustic anisotropy makes it possible to create and implement new configurations of acousto-optic interaction. In particular, a high degree of anisotropy of the acoustic wave phase velocity leads to significantly larger values of the acoustic energy walk-off angle. These high values of the angle between the phase and group velocities allow, for example, efficient implementation of the collinear diffraction mode [
4,
5].
It is known that due to the photoelastic effect, the acoustic perturbation in the medium leads to the induced phase grating in the material. In this case, the polarization vector of the acoustic wave determines the effective photoelastic constant, a value that directly affects the diffraction efficiency [
6]. Information about the polarization composition of the acoustic wave is therefore important for practical applications. It is interesting to consider the conditions under which the acoustic wave becomes purely longitudinal or purely transverse, as well as the conditions under which the wave changes the polarization type.
Monocrystalline materials take a special place in the field of acousto-optics. Among the most common and frequently used media are crystals of tellurium (Te), lithium niobate (LiNbO
3), and TAS (Tl
3AsSe
3), as well as various mercury compounds [
7,
8,
9,
10]. Nowadays though, a paratellurite (TeO
2) crystal is used as a working material in acousto-optic devices most of the time.
A paratellurite crystal is known for its unusual acoustic properties, which distinguish it from other monocrystalline media. As an example,
Figure 1 shows the spatial distributions of the inverse phase velocities of acoustic waves and their acoustic energy walk-off angles. It is evident that paratellurite has extremely low values of the phase velocity of acoustic waves propagating along specific directions. Thus, the minimum phase velocity of a slow shear acoustic mode propagating along one direction is 616 m/s [
11]. To note is the strong acoustic anisotropy of the medium, which leads to record-high angles of acoustic energy walk-off. The maximum value of the acoustic walk-off angle in the paratellurite crystal reaches 74° [
12].
A paratellurite crystal is the standard of modern acousto-optics, and it can be used in creating a wide range of acousto-optic devices. However, with a transparency range of 0.35–6.5 μm [
13], paratellurite cannot be used in creating acousto-optic devices operating in the ultraviolet and far-infrared regions of the electromagnetic spectrum. Moreover, the physical properties of paratellurite are determined and they cannot be changed, so the value of the acousto-optic effect in the medium is fixed. Thus, an imperative task of acousto-optics today is not merely the search for acousto-optically effective materials but also the development of innovative artificial media, whose properties could be controlled, depending on the unique challenges in creating devices. In this study, we turn to acoustic metamaterials, namely phononic crystals [
14,
15,
16].
A phononic crystal is a periodic structure comprising materials with various acoustic properties [
17]. The propagation of acoustic waves in such a crystal is described by the dispersion relation, also known as the band structure. An important feature of the band structure is the presence of band gaps, which are frequency ranges in which the propagation of acoustic waves is forbidden [
18]. In addition, under certain conditions, the dispersion curves have a negative slope, which results in negative refraction [
19,
20]. Moreover, the two-dimensional structures may be applied for tunable broadband polarization converters [
21], cloaking sensors [
22], and plasmonic photodetectors [
23]. The characteristics of phononic crystals are key in studying these materials. However, the periodic nature of phononic crystals predetermines a number of other features, such as strong acoustic anisotropy of the material, low phase velocities of acoustic waves, and large angles of acoustic walk-off [
24,
25,
26]. At the same time, the presence of periodicity and the variation in phononic crystal parameters make it possible to artificially impose various acoustic properties even on isotropic materials.
This paper represents a theoretical study of two-dimensional phononic crystals made of fused silica. The distributions of all the main acoustic characteristics of practical importance in the creation of acousto-optic devices are obtained. The distributions of the inverse phase velocities of acoustic modes of a phononic crystal, the directions of their polarizations, and the acoustic energy walk-off angles are presented. Calculations of the acoustic characteristics are carried out by analyzing the phononic crystal band structure. The dispersion dependences of acoustic waves propagating in a phononic crystal are obtained by solving the elastic wave equation in its variational form, as described in [
17]. The work presents the results of studies conducted to assess the possibility of using the acoustic properties of phononic crystals in creating acousto-optic devices.
2. Theoretical Model and Formulations
This paper analyzes a two-dimensional square-lattice phononic crystal based on fused silica. The model of the phononic crystal considered in the work is shown in
Figure 2. The phononic crystal represents a sample of fused silica with periodically located cylindrical through holes filled with air (
Figure 2a). It is considered an infinite, perfectly periodic medium.
The spatial periodicity of a phononic crystal is characterized by its unit cell, as shown in
Figure 2b. The most important geometric parameters of the unit cell are the diameter of its cylindrical holes, d, and the lattice constant, a, which defines the length of the unit cell; the region Ω denotes the part of the unit cell taken up by fused silica.
Fused silica is chosen as a matrix material forming the phononic crystal, due to its wide availability and ease in manufacturing. In addition, fused silica is known for its low attenuation of acoustic waves [
27]. Thus, when considering the problem of propagation of elastic waves in a phononic crystal, fused silica is considered an ideal medium and losses due to the propagation of acoustic waves are not considered. A possible approach to accounting for the viscoelastic properties of a phononic crystal, leading to the attenuation of acoustic waves, is presented in [
28]. In our work, we also assume that the propagation of acoustic waves in a phononic crystal occurs only in the material of fused silica. The excitation and propagation of acoustic waves in cylindrical holes can be neglected due to the relatively low acoustic impedance of the air.
The theoretical model is based on the elastic wave equation of the form:
where
ui is an elastic wave displacement vector,
Tij is a Cauchy stress tensor characterizing the reaction of the medium to the applied mechanical stress,
fi is the bulk density of external forces, and ρ is the bulk density of the medium in which the elastic wave propagates. The indices take values of
i,j = 1,2,3 and obey the Einstein convention.
The degree of deformation of the medium under the impact of an elastic wave propagating in it is characterized by the strain tensor
Skl, which is defined as:
In the limit of the small strains, the stress tensor
Tij is proportional to the strain tensor
Skl of the medium:
where
cijkl is an elastic constants tensor of the medium. The latest equation represents Hooke’s law in the case of an anisotropic medium, in which stresses and strains are related linearly.
Given that the
cijkl tensor is symmetric with respect to the second pair of indices, in the absence of external forces to the bulk material, the wave equation (Equation (1)) can be written as follows:
In the case of a phononic crystal, the elastic moduli cijkl(r) and the bulk density ρ(r) are step functions of coordinates. Since the propagation of elastic waves occurs only in the material of fused silica, these functions take on the values of the corresponding constants of fused silica in the Ω region, and they are equal to zero in the region of cylindrical holes.
The spatial periodicity of a phononic crystal may be accounted for using Bloch’s theorem. In this case, the general solution of Equation (4) has the form of Bloch waves:
where the amplitude
ũi(
r) is a periodic function of coordinates, defined in the Ω region of the phononic crystal unit cell, and equal to zero in the region of cylindrical holes;
ω = 2πf is the angular frequency of the elastic wave; and
k is the wave vector of the elastic wave.
The vector ũ(r) sets the amplitude of the elastic wave and defines the direction of its polarization. The absolute value of the wave vector |k| = ω/V sets the phase velocity of acoustic waves propagating in a phononic crystal. Based on the axial symmetry of the unit cell, we limit the analysis to the wave vector in the form k = (kx, ky, 0).
The strain tensor of a periodic medium also, according to its definition (Equation (2)), takes the form of a Bloch wave:
where the amplitudes
ij equals:
The explicit form of solution (Equation (5)) allows obtaining the variational formulation of the wave equation (Equation (4)) for infinite periodic media. The resulting variational problem can be solved by the finite element method (FEM), which makes it possible to find the dispersion dependences of acoustic waves. The core idea of the method is to approximate a solution belonging to some infinite-dimensional function space (e.g., Sobolev space) by a linear combination of the fixed-basis vectors of a subspace of the space. The FEM approximation depends on the choice of a subspace and the associated basis. Usually, a space of piecewise polynomial functions is chosen as such a subspace. In our research, we used the Lagrange P
2-elements as a basis. After choosing the basis functions, the FEM comes down to the optimization problem of minimizing the approximation error of some functional (given by a weak formulation of the original differential equation) on the FEM mesh. An example of the FEM mesh used in our calculcations for a phononic crystal with
d/a = 0.6 is shown in
Figure 3. According to [
17], the described variational problem takes the following form:
Here, v(r) is an arbitrary test function that belongs to the same functional space of finite elements as the solution ũ(r). Introduction of this function is necessary for solving the problem using the finite element method. In this case, the Voigt notation is used for tensor quantities, so the indices I ↔ (ij) and J ↔ (kl) take values from 1 to 6 according to the rule 1 ↔ (11), 2 ↔ (22), 3 ↔ (33), 4 ↔ (23), 5 ↔ (13), 6 ↔ (12). Considering the thermodynamic relations, the components of the six-dimensional deformation amplitudes vector are defined as I = (2 − δij) • ij, where δij is the Kronecker symbol.
The integral equation (Equation (8)) is a generalized eigenfunction and eigenvalue problem. The eigenvalue in this case is λ = ω2. The set of eigenvectors ũλ determines the amplitudes and polarization directions of acoustic waves of frequency ω, which can propagate in a phononic crystal. According to Equation (7), the amplitude of the deformation vector is directly affected by the components of the wave vector. Thus, by changing the value of the wave vector within the first zone and solving the integral equation (Equation (8)) for each individual value of the wave vector k, one can obtain the dispersion relation ω(k).
To identify the main features of the solution to the integral equation (Equation (8)), it has to be brought to its dimensionless form. The following dimensionless quantities are introduced:
where ρ′ and
c′
IJ are dimensionless bulk density and elastic constants, respectively, carrying the numerical values of the corresponding physical constants of fused silica, and ρ
0 and c
0 are quantities that characterize the dimensions of those physical constants. In dimensionless quantities, the integral equation (Equation (8)) takes the form:
where the integration is carried out over the region ζ of the unit cell, which is made dimensionless, as shown in
Figure 2c. The explicit form of the integral equation (Equation (10)) allows making several important observations.
First, the eigenvalues of this equation, and therefore also the eigenvalues of Equation (8), are not as dependent on the specific values of a and d as on their ratio d/a. Theoretically, the normalized hole diameter can be infinitely close to unity. However, the maximum d/a ratio considered in this work is 0.8. This choice pertains to the technological difficulties in creating phononic crystals, where the normalized hole diameter exceeds this value. It should be noted that increasing the diameter of the holes reduces the space taken up by the matrix material. This, in turn, leads to a decrease in the acousto-optic interaction region, which results in a lower diffraction efficiency. In addition, for a given geometry of a phononic crystal, that is, for a fixed region of integration ζ, the product of ω•a = const. Thus, a change in the unit cell constant a of a phononic crystal leads to the corresponding scaling of the dispersion dependences. Further in this work, for the sake of presentation clarity, we consider the solution of the integral equation (Equation (8)) in dimensional quantities. The unit cell constant a is considered equal to 10 μm.
3. The Method for Calculations of the Phononic Crystal Acoustic Characteristics
A consistent solution of the generalized eigenfunction and eigenvalue problem (Equation (8)) for different values of the wave vector
k makes it possible to obtain the dispersion dependence ω(k) of acoustic waves propagating in a phononic crystal.
Figure 4 shows the band structure of a phononic crystal with a normalized hole diameter
d/a = 0.8.
Figure 4a shows the dispersion dependences calculated by the finite element method. Only the first three dispersion surfaces are displayed. They correspond to three different acoustic modes with a frequency of up to 200 MHz. Higher-order surfaces are above the indicated frequency. There are no absolute band gaps in the frequency range of up to 200 MHz and this crystal geometry, as evident from
Figure 4a.
The dispersion surfaces of the first three acoustic modes of the phononic crystal originate from the center of the first Brillouin zone. In the lower frequency range, up to 50 MHz, these surfaces are mutually exclusive and monotonic and do not intersect. With higher frequencies exceeding 50 MHz, the surfaces are no longer monotonic. The emerging frequency regions with a negative slope may lead to a negative refraction effect. At frequencies above 100 MHz, the second and third surfaces begin to intersect. Such a complex structure is explained by the relatively high
d/a ratio. It is also evident from
Figure 4a that the dispersion surfaces are symmetric with respect to the first irreducible Brillouin zone Γ-X-M-Γ, which is explained by the symmetry of the unit cell.
The dispersion relation ω(k) connects the frequencies and wave vectors of acoustic waves that can propagate in a medium. By fixing the frequency at ω0 = 2πf0, it is possible to obtain sets of wavenumbers kx(f0) and ky(f0) of all acoustic waves of frequency f0 that can propagate in a phononic crystal. At a given ultrasound frequency f0, the wavenumbers kx and ky are proportional to the components of the inverse phase velocity of acoustic waves, since kx = 2πf0/Vx and ky = 2πf0/Vy. Thus, the line of the frequency contour f0 of the dispersion dependence defines the cross sections of the acoustic slowness surface S(φ) in the XY plane of the phononic crystal, where φ is the polar angle. In an acousto-optic device, the isofrequency f0 is the frequency of ultrasound excited by the piezoelectric transducer.
Figure 4a shows the contour lines of the isofrequency
f0 = 50 MHz, while the projection of those contour lines to the first Brillouin zone can be seen in
Figure 4b. Similar to the dispersion surfaces, the isofrequency contour lines are symmetric with respect to the first irreducible Brillouin zone Γ-X-M-Γ. For isofrequencies
f0 < 50 MHz, the contour line of each individual acoustic mode is a closed curve. For isofrequencies
f0 > 50 MHz, where the dispersion surfaces become nonmonotonic, the contour lines become piecewise interval lines. The calculation of acoustic characteristics at these higher isofrequencies requires a separate analysis. This work presents the case in which the isofrequency is chosen to be
f0 = 50 MHz.
The spatial periodicity of a phononic crystal leads to anisotropy of the physical characteristics of acoustic waves propagating in a phononic crystal. In particular, the vectors of the phase and group velocities of acoustic waves become non-codirectional. The explicit form of the acoustic slowness curves S(φ) allows finding the angle ψ of the acoustic energy walk-off. The acoustic energy walk-off angle between the Pointing vector and the acoustic wave phase velocity, according to [
6], can be found as:
In acousto-optics, one of the most important characteristics of an acoustic wave is the direction of its polarization. The eigenvectors
ũλ(
r) of the integral problem (Equation (8)) define the amplitudes and polarization directions of ω frequency acoustic waves, which can propagate in a phononic crystal. Since eigenvectors
ũλ(
r) are periodic functions of coordinates, the polarization direction of each acoustic mode varies within the unit cell. Therefore, from a practical point of view, it seems necessary to carry out the calculation of the polarization vector
ũ0 averaged over the unit cell region. The components of the averaged polarization vector can be obtained as:
where the index
p =
x,
y,
z.
According to Equation (12), the squared values of the integral equation (Equation (8)) solution components are averaged over the unit cell region of the phononic crystal and then normalized. However, only the squared values of the averaged polarization vector components can be obtained this way. The actual components of the averaged polarization vector ũ0 are defined to within an overall sign, as described in Equation (12). It is clear that in the limit d/a → 0, the averaged solution of Equation (8) should converge to that of solid isotropic fused silica. The requirement of the physical behavior in the mentioned limit makes it possible to choose the sign in Equation (12). We also assume that with a small change in the phononic crystal geometry, that is, with a small change in d/a, the change in the polarization vector direction is continuous. With this assumption, the signs of the components of the averaged polarization vector in Equation (12) can be determined for any values of the ratio d/a.