## Controlling Voigt waves by the Pockels effect 1

Tom G. Mackay 2

School of Mathematics and Maxwell Institute for Mathematical Sciences University of Edinburgh, Edinburgh EH9 3JZ, UK

and

NanoMM - Nanoengineered Metamaterials Group

Department of Engineering Science and Mechanics Pennsylvania State University, University Park, PA 16802-6812, USA

## Abstract

Voigt wave propagation was investigated in a homogenized composite material (HCM) arising from a porous electro-optic host material infiltrated by a fluid of refractive index n a . The constitutive parameters of the HCM were estimated using the extended Bruggeman homogenization formalism. Numerical studies revealed that the directions which support Voigt wave propagation in the HCM could be substantially controlled by means of an applied dc electric field. Furthermore, the extent to which this control could be achieved was found to be sensitive to the porosity of the host material, the shapes, sizes and orientations of the pores, as well as the refractive index n a . These findings may be particularly significant for potential technological applications of Voigt waves, such as in optical sensing.

Keywords: Voigt waves, electro-optic materials, Bruggeman homogenization formalism, Pockels effect

## 1 Introduction

The setting for this numerical study is electromagnetic plane-wave propagation in anisotropic dielectric materials. Usually in such materials two plane waves with different phase speeds and orthogonal polarizations propagate in a given direction [1]. However, circumstances can be such that these two plane waves coalesce to form a single plane wave, known as a Voigt wave [2, 3, 4, 5]. A key characteristic of a Voigt wave is that its amplitude is linearly dependent upon propagation distance. Indeed, this linear dependence on propagation distance could provide the basis for experimental studies and technological applications that exploit Voigt wave propagation. Mathematically, the requirements for Voigt wave propagation are met when an eigenvalue of the corresponding plane-wave propagation matrix has an algebraic multiplicity which exceeds its geometric multiplicity [6].

While certain pleochroic minerals can support Voigt wave propagation [2], greater scope for realizing this singular form of propagation is offered by engineered materials [7, 8]. In particular, homogenized composite materials (HCMs) may be conceptualized which support Voigt wave propagation whereas the component materials from which they arise do not [9, 10]. By judicious design of the structure and selection of the constitutive properties of the component materials, the directions in which Voigt waves propagate in a HCM may be controlled. This ability to control the directions for Voigt wave propagation may be attractive from the point of view of possible technological applications. Recently the potential that Voigt waves offer for optical sensing applications was highlighted [11]. The scenario considered involved an HCM arising from a porous biaxial-dielectric host material which was infiltrated by a fluid of refractive index n a . Numerical studies showed that the directions which supported Voigt wave propagation in the HCM could be acutely sensitive to the refractive index n a ; indeed, sensitivities of up to 300 ◦ per RIU were reported.

1 This manuscript is an extension of the conference paper: T.G. Mackay, 'Voigt waves in electro-optic homogenized composite materials,' SPIE Optics &amp; Photonics 2014: Nanostructured Thin Films VII , San Diego, USA.

2 E-mail: T.Mackay@ed.ac.uk

In the present study, we consider the prospect of controlling Voigt wave propagation in a HCM at the post-fabrication stage, by means of an applied dc electric field. The HCM here arises from two component materials: a porous electro-optic host material which is infiltrated by a fluid of refractive index n a . An extended version [12] of the well-established Bruggeman homogenization formalism [13, 14] is employed to estimate the constitutive parameters of the HCM. The Bruggeman homogenization formalism was chosen chiefly because (i) there is no restriction on volume fraction of the component materials (unlike in the popular Maxwell Garnett formalism which is restricted to small volume fractions [14]); and (ii) the shapes and orientations of the particles which make up the two component materials can be independently varied (unlike in the Maxwell Garnett formalism wherein only the shape and orientation of the designated 'inclusion' particles can be varied, and unlike in the strong-permittivity-fluctuation-theory formalism wherein the particles of both component materials are required to have the same shape and orientation [14]).

As regards the notation adopted in this paper: vectors are denoted by bold typeface, and the addition of the ˆ symbol indicates a unit vector. Accordingly, the unit vectors aligned with the Cartesian axes are written as ˆ, x ˆ, and ˆ. y z Normal typeface combined with double underlining signifies a 3 × 3 dyadic. Thus, I = ˆ ˆ + ˆ ˆ +ˆ ˆ is the identity 3 xx yy zz × 3 dyadic and 0 is the null 3 × 3 dyadic. The dyadic transpose is identified by the superscript T . The determinant of a dyadic is delivered by the dyadic operator 'det'. Blackboard bold typeface combined with double underlining signifies a 6 × 6 dyadic. The permittivity and permeability of free space are represented by the symbols /epsilon1 0 and µ 0 , respectively; the free-space wavenumber is expressed as k 0 = ω √ /epsilon1 0 µ 0 where ω is the angular frequency.

## 2 Voigt waves

Our study concerns the propagation of Voigt waves in a homogenized composite material (HCM). As described later in § 3.2, here the HCM is an anisotropic dielectric material whose electromagnetic properties are characterized by the symmetric dyadic /epsilon1 HCM . In general, there are two directions in which Voigt waves may propagate in the HCM. We determine these two directions in an indirect fashion: Voigt wave propagation is considered along the z coordinate axis for all possible orientations of the HCM. Accordingly, we introduce the HCM permittivity dyadic in a rotated coordinate frame

$$\begin{array} { r c l } \dots \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \rstackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \left ( 1 } \, \\ \end{array} \\ \text{where the orthogonal rotation dyadics}$$

where the orthogonal rotation dyadics

$$\underline { \underline { R } } _ { y } ( \nu ) & = \cos \nu \left ( \hat { x } \hat { x } + \hat { z } \hat { z } \right ) + \sin \nu \left ( \hat { z } \hat { x } - \hat { x } \hat { z } \right ) + \hat { y } \hat { y } \\ \underline { \underline { R } } _ { z } ( \nu ) & = \cos \nu \left ( \hat { x } \hat { x } + \hat { y } \hat { y } \right ) + \sin \nu \left ( \hat { x } \hat { y } - \hat { y } \hat { x } \right ) + \hat { z } \hat { z } \right ) \\ \cdot \cdot \cdot \cdot \cdot$$

and α β , , and γ are the three Euler angles [15].

For the material characterized by the permittivity dyadic ˜ /epsilon1 HCM ( α, β, γ ), Voigt wave propagation parallel to the z axis is supported provided that [16]:

- (i) Y ( α, β, γ ) = 0 and
- (ii) W α,β,γ ( ) = 0,

/negationslash where the scalar quantities

$$Y ( \alpha, \beta, \gamma ) \ = \ \ & \xi _ { 1 3 } ^ { 4 } + \tilde { \epsilon } _ { 2 3 } ^ { 4 } - 2 \tilde { \epsilon } _ { 2 3 } \tilde { \epsilon } _ { 3 3 } \left [ 2 \tilde { \epsilon } _ { 1 2 } \tilde { \epsilon } _ { 1 3 } - ( \tilde { \epsilon } _ { 1 1 } - \tilde { \epsilon } _ { 2 2 } ) \, \tilde { \epsilon } _ { 2 3 } \right ] + \left [ \left ( \tilde { \epsilon } _ { 1 1 } - \tilde { \epsilon } _ { 2 2 } \right ) ^ { 2 } + 4 \tilde { \epsilon } _ { 1 2 } ^ { 2 } \right ] \tilde { \epsilon } _ { 3 3 } ^ { 2 } \\ + 2 \tilde { \epsilon } _ { 1 3 } \left \{ \tilde { \epsilon } _ { 2 3 } ^ { 2 } \tilde { \epsilon } _ { 1 3 } - [ 2 \tilde { \epsilon } _ { 1 2 } \tilde { \epsilon } _ { 2 3 } + ( \tilde { \epsilon } _ { 1 1 } - \tilde { \epsilon } _ { 2 2 } ) \, \tilde { \epsilon } _ { 1 3 } ] \tilde { \epsilon } _ { 3 3 } \right \} & & ( 4 ) \\ & & 2$$

and

$$W ( \alpha, \beta, \gamma ) & = \tilde { \epsilon } _ { 1 2 } \tilde { \epsilon } _ { 3 3 } - \tilde { \epsilon } _ { 1 3 } \tilde { \epsilon } _ { 2 3 } \,. \quad & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & && & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &  & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &
 & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & \\ & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & \ & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & ; & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &. & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & } & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & - & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & _ & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & = & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & + & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &, & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &	 & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &

 & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & b & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & ^ & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &$$

These two necessary and sufficient conditions emerge from an eigenanalysis of the corresponding propagation matrix for plane-wave propagation [16]. In particular, these conditions can be satisfied by certain biaxial dielectric materials but not by isotropic or uniaxial dielectric materials.

## 3 Extended Bruggeman homogenization formalism

## 3.1 Component materials

Our study is based on the homogenization of two particulate component materials, namely component material a and component material b . Component material a is simply envisaged as an isotropic dielectric fluid characterized by the permittivity dyadic /epsilon1 ( a ) = /epsilon1 0 /epsilon1 ( a ) I in the optical regime; its relativity permittivity /epsilon1 ( a ) = n 2 a . Component material b is an electro-optic material that exhibits the Pockels effect in the optical regime. Its linear electro-optic properties are conventionally characterized via its inverse permittivity dyadic [17, 18]. Retaining only the first-order components of the dc electric field (which is the usual approximation in electro-optics [19]), the permittivity dyadic for component material b may be written as [20]

$$\underline { \epsilon } ^ { ( b ) } & \approx \ \epsilon _ { 0 } \left \{ \sum _ { k = 1 } ^ { 3 } \left [ \epsilon _ { k } ^ { ( b ) } \left ( 1 - \epsilon _ { k } ^ { ( b ) } s _ { k } \right ) \hat { u } _ { k } \hat { u } _ { k } \right ] \\ & - \ \epsilon _ { 2 } ^ { ( b ) } \epsilon _ { 3 } ^ { ( b ) } s _ { 4 } \left ( \hat { u } _ { 2 } \hat { u } _ { 3 } + \hat { u } _ { 3 } \hat { u } _ { 2 } \right ) - \epsilon _ { 1 } ^ { ( b ) } \epsilon _ { 3 } ^ { ( b ) } s _ { 5 } ( \hat { u } _ { 1 } \hat { u } _ { 3 } + \hat { u } _ { 3 } \hat { u } _ { 1 } ) - \epsilon _ { 1 } ^ { ( b ) } \epsilon _ { 2 } ^ { ( b ) } s _ { 6 } ( \hat { u } _ { 1 } \hat { u } _ { 2 } + \hat { u } _ { 2 } \hat { u } _ { 1 } ) \right \}, \\ \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot$$

under the assumption that

$$\frac { 1 } { \epsilon _ { 0 } } \left \{ | \epsilon _ { 1 } ^ { ( b ) } |, | \epsilon _ { 2 } ^ { ( b ) } |, | \epsilon _ { 3 } ^ { ( b ) } | \right \} _ { \max } \ \{ | s _ { 1 } |, | s _ { 2 } |, | s _ { 3 } |, | s _ { 4 } |, | s _ { 5 } |, | s _ { 6 } | \} _ { \max } \ll 1 \,. \\ ( \delta \, \text{$two$ orientation of the unit} \text{ note}$$

In Eq. (6) the orientations of the unit vectors

$$\hat { u } _ { 1 } = - ( \hat { x } \cos \phi _ { b } + \hat { y } \sin \phi _ { b } ) \cos \theta _ { b } + \hat { z } \sin \theta _ { b } \\ \hat { u } _ { 2 } = \hat { x } \sin \phi _ { b } - \hat { y } \cos \phi _ { b } \quad \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \
\i e \, \text{angles $\theta_{b} \in [0,\pi]$ and $\phi_{b} \in [0,2\pi], are determined by the crystallographic structure of} \\ \i e \, \text{endency of commonent material} \, b ^ { \prime } s \text{ constitutive parameters union a uniform} \, \text{dec} \, \text{electric}$$

as specified by the angles θ b ∈ [0 , π ] and φ b ∈ [0 , 2 π ], are determined by the crystallographic structure of the material. The dependency of component material b 's constitutive parameters upon a uniform dc electric field E dc is captured by the scalar parameters

$$s _ { j } = \sum _ { k = 1 } ^ { 3 } r _ { j k } \hat { \mathfrak { u } } _ { k } \cdot \mathbf E ^ { d c }, \quad \ \ j \in \{ 1, 2, 3, 4, 5, 6 \}, \quad \ \ ( 9 ) \\ \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \colon \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \, \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \dots$$

which are expressed in terms of the electro-optic coefficients r jk ( j ∈ { 1 2 3 4 5 6 , , , , , , } k ∈ (1 , 2 3 , } ). Thus, depending upon the relative values of the principal permittivity scalars /epsilon1 ( b ) 1 2 3 , , , component material b may be an isotropic, uniaxial or biaxial material; furthermore, depending upon the relative values of the electro-optic coefficients, it may belong to one of 20 crystallographic classes of point group symmetry [14].

Both component materials are composed of spheroidal particles. These component spheroids are randomly distributed with the spheroids of component material a occupying the volume fraction f a and those of component material b occupying the volume fraction f b = 1 -f a . It is assumed that all component a spheroids have the same shape and orientation, and all component b spheroids have the same shape and

orientation. The surface of each component spheroid, relative to its centroid, is prescribed by the position vector

$$r _ { \ell } = \eta _ { \ell } \underline { U } _ { \ell } \cdot \hat { \mathbf r }, \quad ( \ell = a, b ). \quad \ \ \ \ \ ( 1 0 )$$

Here the position vector ˆ prescribes the surface of the unit sphere; the shape and orientation of the compor nent spheroid are encapsulated by the real-symmetric dyadic U /lscript ; and the linear dimensions of the component spheroid are characterized by the size parameter η /lscript &gt; 0. The homogenization of component materials a and b is founded on the notion that η /lscript is much smaller than the wavelengths involved. But, in extended homogenization formalisms, η /lscript need not be vanishingly small [12]. Accordingly, in the following numerical studies, the range 0 2 . &lt; k η 0 /lscript &lt; 0 3 is used. . For simplicity, we assume that η a ≡ η b ; and henceforth η is written in lieu of η /lscript ( /lscript = a, b ).

For definiteness, we choose the axis of rotational symmetry for the component b spheroids to be aligned with the x coordinate axis. Thus, the surface dyadic for the component material b spheroids may be expressed as

$$\underline { U } _ { = b } = \frac { 1 } { \frac { 3 } { \sqrt { U _ { x } U ^ { 2 } } } } \left [ U _ { x } \hat { x } \hat { x } + U \left ( \hat { y } \hat { y } + \hat { z } \hat { z } \right ) \right ], \quad ( U _ { x }, U > 0 ).$$

The axis of rotational symmetry for the component material a spheroids is rotated by an angle ϕ in the xy plane relative to the axis of rotational symmetry for the component material b spheroids. Thus, the surface dyadic for the component material a spheroids may be expressed as

$$\underline { U } _ { a } = \underline { R } _ { z } ( \varphi ) \cdot \underline { U } _ { b } \cdot \underline { R } _ { z } ^ { T } ( \varphi ).$$

## 3.2 Homogenized composite material

In the long-wavelength regime, the random mixture of component materials a and b may be regarded as a homogenized composite material (HCM). The electromagnetic properties of this HCM are characterized by its symmetric permittivity dyadic /epsilon1 HCM . We utilize the extended Bruggeman formalism [12, 21] to estimate /epsilon1 HCM . This approach is based on the nonlinear dyadic equation [22]

$$f _ { a } \left \{ \left ( \underline { \epsilon } ^ { ( a ) } - \underline { \epsilon } _ { H C M } \right ) \cdot \left [ \underline { L } + \underline { D } _ { a } \cdot \left ( \underline { \epsilon } ^ { ( a ) } - \underline { \epsilon } _ { H C M } \right ) \right ] ^ { - 1 } \right \} + f _ { b } \left \{ \left ( \underline { \epsilon } ^ { ( b ) } - \underline { \epsilon } _ { H C M } \right ) \cdot \left [ \underline { L } + \underline { D } _ { a } \cdot \left ( \underline { \epsilon } ^ { ( b ) } - \underline { \epsilon } _ { H C M } \right ) \right ] ^ { - 1 } \right \} = \underline { \epsilon } \cdot \\ T h a \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \end{matrix}$$

The depolarization dyadics D a,b herein encapsulate the electromagnetic responses of U a,b -shaped spheroids embedded in the HCM. Further details are provided in the Appendix.

## 4 Numerical investigations

The question that we now turn to is: to what extent can the directions that support Voigt wave propagation be influenced by the application of a dc electric field E dc ? We do so by means of representative numerical calculations.

## 4.1 Component materials

We consider a porous electro-optic host material b , characterized by the permittivity dyadic /epsilon1 ( b ) , which is infiltrated by a fluid of refractive index n a . For component material b , we choose potassium niobate which is specified by [23]: /epsilon1 ( b ) 1 = 4 72, . /epsilon1 ( b ) 2 = 5 20, . /epsilon1 ( b ) 3 = 5 43, . r 13 = 34 × 10 -12 m V -1 , r 23 = 6 × 10 -12 m V -1 , r 33 = 63 4 . × 10 -12 m V -1 , r 42 = 450 × 10 -12 m V -1 , r 51 = 120 × 10 -12 m V -1 , and all other r jk ≡ 0. And for component material a we consider the range n a ∈ [1 05 . , 1 3]. . The eccentricity parameter ρ is introduced as a gauge of the shape of the component spheroids: we let U x = 1 + ρ , U = 1 -( ρ/ 20), and consider the range ρ ∈ [2 , 4].

In order to emphasize the degree of electrical control that may be attained over the Voigt wave directions, we explore cases where the Pockels effect is most easily discernable. To this end, the principal crystallographic

axis of component material b is chosen to be aligned with its component spheroids; in addition, the dc electric field E dc is also chosen to be aligned with the component b spheroids. Accordingly, the crystallographic angles θ b = φ b = 0; and we have E dc = E dc 3 ˆ u 3 with the range for the electric field component E dc 3 being such that the inequality (7) is satisfied.

## 4.2 Extended Bruggeman estimates of /epsilon1 HCM

Before exploring the directions which support Voigt wave propagation in relation to a dc electric field, we first consider how such an electric field effects the constitutive parameters of the HCM. For the constitutive parameter values chosen for our numerical investigations, the HCM is a biaxial dielectric material characterized by a symmetric permittivity dyadic of the general form [24]

$$\underline { \epsilon } _ { H C M } = \epsilon _ { 0 } \left [ \epsilon _ { x } \hat { x } \hat { x } + \epsilon _ { y } \hat { y } + \epsilon _ { t } \left ( \hat { x } \hat { y } + \hat { y } \hat { x } \right ) + \epsilon _ { z } \hat { z } \hat { z } \right ].$$

In Fig. 1 the real and imaginary parts of the extended Bruggeman estimates of the relative permittivity parameters /epsilon1 x,y,z,t are plotted versus the dc electric field component E dc 3 , and the spheroid orientation angle ϕ . For these calculations, the volume faction f a = 0 3, the refractive index . n a = 1 1, the size parameter . η = 0 3 . /k 0 , and the eccentricity parameter ρ = 3. We see that both the real and imaginary parts of /epsilon1 x,y,z decrease in an approximately linear manner as E dc 3 increases from - × 1 10 8 Vm -1 to 1 × 10 8 Vm -1 , albeit this effect is much more obvious for /epsilon1 z than it is for /epsilon1 x,y . The effect of E dc 3 on both the real and imaginary parts of /epsilon1 t is barely discernable in Fig. 1. Both the real and imaginary parts of /epsilon1 x,y vary markedly as ϕ increases from 0 ◦ to 90 ◦ . However, the relative permittivity parameter /epsilon1 z seems to be largely unaffected by the rotation in the xy plane represented by ϕ . The off-diagonal relative permittivity parameter /epsilon1 t = 0 at ϕ = 0 ◦ and 90 ; both the real and imaginary parts of ◦ /epsilon1 t increase and then decrease markedly as ϕ increases from 0 ◦ to 90 ◦ .

In addition to the dc electric field and the orientation of the component spheroids, the electromagnetic properties of HCM are also sensitive to the shapes and sizes of the component spheroids, the volume fraction, and the constitutive parameters of the component materials. These matters have been explored using the extended Bruggeman homogenization formalism in earlier studies, to which the reader is referred for further details [25, 26].

## 4.3 Electrical control over directions that support Voigt wave propagation

/negationslash

/negationslash

Now we investigate the effects of a dc electric field on the directions which support Voigt wave propagation in the HCM. These Voigt wave directions are represented by the Euler angles α β , , and γ for which Y ( α, β, γ ) = 0 and W α,β,γ ( ) = 0. Notice that the angular coordinate γ may be eliminated from our investigations since Voigt wave propagation parallel to the z axis (in the rotated coordinate system) is independent of rotation about that axis. For the HCM specified by a permittivity dyadic of the form given in Eq. (14), there are generally two distinct orientations for which the Voigt wave conditions Y = 0 and W = 0 are satisfied. We label these two orientations with the angular coordinates α = α 1 2 , and β = β 1 2 , ; the corresponding values of the quantity W are written as W 1 2 , .

## 4.3.1 Orientation of component spheroids

We begin by considering the effect of the dc electric field on the Voigt wave directions in relation to the orientation of the component spheroids. In Fig. 2, the angles α 1 2 , and β 1 2 , along with the quantities | W 1 2 , | are plotted versus E dc 3 for three different orientations of the component material a spheroids as given by ϕ = 30 ◦ (green, dashed curves), 60 ◦ (red, solid curves), and 90 ◦ (blue, broken dashed curves). Here the refractive index n a = 1 1, the volume fraction . f a = 0 3, the size parameter . η = 0 3 . /k 0 and the eccentricity parameter ρ = 3. Clearly the directions that support Voigt wave propagation are highly sensitive to the dc electric field, for all orientations of the component material a spheroids considered. Specifically, the α angular coordinate varies by as much as 38 ◦ , and the β angular coordinate varies by as much as 12 ◦ , as the

value of E dc 3 ranges from - × 1 10 8 Vm -1 to 1 × 10 8 Vm -1 . Generally, greater changes in the Voigt wave directions are observed for smaller values of the spheroid orientation angle ϕ . Also, the angular coordinate β varies relatively little for E dc 3 ∈ (0 , 1) × 10 8 Vm -1 , but varies considerably more for E dc 3 ∈ ( -1 0) , × 10 8 Vm -1 . We note that the quantities | W 1 2 , | are nonzero for all calculations, in accordance with requirements for Voigt wave propagation. Other than being nonzero, the quantities | W 1 2 , | are of no obvious physical significance; therefore, plots of these quantities shall be omitted henceforth.

## 4.3.2 Shape of component spheroids

We repeat the calculations of § 4.3.1 but now fixing the orientation of the component material a spheroids at ϕ = 60 ◦ and allowing the eccentricity parameter ρ to vary. The results are presented in Fig. 3 where the angles α 1 2 , and β 1 2 , are plotted versus E dc 3 for ρ = 2 (green, dashed curves), 3 (red, solid curves), and 4 (blue, broken dashed curves). The sensitivity of the Voigt wave directions to the dc electric field is significantly modulated by the degree of eccentricity of the component spheroids. Specifically, the angular coordinates α 1 2 , are moderately sensitive to ρ . While the angular coordinates β 1 2 , are relatively insensitive to ρ for E dc 3 ∈ (0 , 1) × 10 8 Vm -1 , these angles are highly sensitive to ρ for E dc 3 ∈ ( -1 0) , × 10 8 Vm -1 . The influence of the shape parameter upon the sensitivity of the Voigt wave directions in relation to E dc 3 appears to be greatest for the smallest value of ρ considered here.

## 4.3.3 Size of component spheroids

Next we turn to the size of the component spheroids, as gauged by the size parameter η . In Fig. 4 plots similar to those in Fig. 3 are provided but here the eccentricity parameter is fixed at ρ = 3 and three different values of size parameter are employed, namely η = 0 2 . /k 0 (green, dashed curves), 0 25 . /k 0 (blue, broken dashed curves), and 0 3 . /k 0 (red, solid curves). The general trends observable in Fig. 4 are similar to those observable in Fig. 3. That is, the size parameter does significantly modulate the dependency of the Voigt wave directions on the dc electric field; the influence of the size parameter is relatively modest for the α 1 2 , angular coordinates but is considerably stronger for the β 1 2 , angular coordinates, especially at lower values of E dc 3 .

## 4.3.4 Volume fraction

The influence of the volume fraction occupied by component material a (or, equivalently, the porosity of component material b ) is now considered. Plots similar to those in Fig. 4 are provided in Fig. 5 but here the size parameter η = 0 3 . /k 0 ; results are plotted for the volume fractions f a = 0 3 (red, solid curves), 0 35 . . (green, dashed curves), and 0 4 (blue, broken dashed curves). The approximate linearity of the relationship . between the α 1 2 , angular coordinates and E dc 3 is not greatly affected by the volume fraction. Furthermore, in the vicinity of E dc 3 ≈ 0 × 10 8 Vm -1 , the angular coordinates β 1 2 , seem to be insensitive to the volume fraction. But for larger and smaller values of E dc 3 , the volume faction has a strong influence on the angular coordinates β 1 2 , in relation to E dc 3 . This influence is strongest at smaller values of f a when E dc 3 is negative-valued and at larger values of f a when E dc 3 is positive-valued.

## 4.3.5 Refractive index of infiltrating fluid

Lastly, let us consider the influence of the refractive index of component material a , namely n a . We repeat the calculations of Fig. 5 but with the volume fraction fixed at f a = 0 3; results are plotted in Fig. 6 for . the refractive index values n a = 1 05 (green, dashed curves), 1 1 (red, solid curves), and 1 3 (blue, broken . . . dashed curves). It may be seen that the refractive index strongly modulates the sensitivity of the Voigt wave directions to the dc electric field. Reminiscent of the general trends in Fig. 3, the refractive index n a has only a modest influence upon the angular coordinates β 1 2 , for E dc 3 ∈ (0 , 1) × 10 8 Vm -1 , but these angles are highly sensitive to n a for E dc 3 ∈ -( 1 0) , × 10 8 Vm -1 . The effects of the refractive index upon the sensitivity of the Voigt wave directions in relation to E dc 3 are generally greatest for the smallest value of n a considered here.

## 5 Closing comments

Voigt wave propagation has been considered for a HCM arising from a porous electro-optic host material, namely potassium niobate, infiltrated by a fluid of refractive index n a . By means of representative numerical calculations, it has been demonstrated that the directions which support the propagation of Voigt waves may be substantially controlled by the application of a dc electric field. Furthermore, the extent to which this control may be achieved is dependent upon the porosity of the host material, the shapes, sizes and orientations of the pores, as well as the refractive index n a .

In order to be physically realistic, the magnitude of the applied dc electric field should be less than the magnitudes of the breakdown electric fields of the materials involved. However, the breakdown electric fields depend upon the morphologies and thicknesses of the materials, with thin films generally capable of withstanding higher electric fields than bulk materials. Typical breakdown fields for potassium niobate thin films are estimated to be about 8 × 10 7 V m -1 [27], while Kingon et al. reported no breakdown when a field of 1 × 10 8 V m -1 was applied to a potassium niobate thin film in their experiments [28]. Accordingly, it seems reasonable here to present results for applied dc electric fields of magnitudes up to 1 0 . × 10 8 V m -1 . Parenthetically, the inequality (7) is satisfied for such electric fields. In a practical device, the maximum dc electric field that could be applied may well be determined by factors other than the breakdown of the electro-optic material, such as how well the electronic components of the device can be shielded from the dc electric field. Furthermore, potassium niobate was chosen here as a representative example of an electrooptic material which has relatively large electro-optic coefficients. If artificial electro-optic materials were to be developed with larger electro-optic coefficients then this would allow the possibility of using dc electric fields of lower magnitude.

The findings reported here may be particulary significant for potential technological applications of Voigt waves. Notably, Voigt waves appear to present promising opportunities for optical sensing applications, since the directions which support Voigt wave propagation in a porous host material can be highly sensitive to changes in the refractive index of a fluid that infiltrates the host material [11]. Thus, it may be envisaged that electrical control could be harnessed to track Voigt wave propagation in an optical-sensor setting. Other possible applications of electrical control of Voigt waves may involve optical switching [29]. Such Voigt wave applications are matters for future study.

## Appendix

In the extended Bruggeman formalism encapsulated by Eq. (13), the depolarization dyadics comprise the sums [30]

$$\underline { \underline { D } } _ { \ell } = \underline { \underline { D } } _ { \ell } ^ { 0 } + \underline { \underline { D } } _ { \ell } ^ { + }, \quad ( \ell = a, b ).$$

The dyadic term D 0 /lscript represents the depolarization contribution arising in the limit η → 0; it is given by the double integral [31, 32]

$$\underline { \underline { D } } _ { \ell } ^ { 0 } = \frac { 1 } { 4 \pi } \int _ { \phi = 0 } ^ { 2 \pi } \int _ { \theta = 0 } ^ { \pi } \frac { \left ( \underline { U } ^ { - 1 } _ { \ell } \cdot \hat { \mathbf q } \right ) \left ( \underline { U } ^ { - 1 } _ { \ell } \cdot \hat { \mathbf q } \right ) \, \sin \theta } { \left ( \underline { U } ^ { - 1 } _ { \ell } \cdot \hat { \mathbf q } \right ) \cdot \underline { \epsilon } _ { H C M } \cdot \left ( \underline { U } ^ { - 1 } _ { \ell } \cdot \hat { \mathbf q } \right ) } \, d \theta \, d \phi, \quad ( \ell = a, b ), \quad ( 1 6 ) \\ \text{ith the unit vector} \, \hat { q } = \sin \theta \cos \phi \hat { x } + \sin \theta \sin \phi \hat { y } + \cos \theta \hat { z }. \text{ The depolarization contribution associated with}$$

with the unit vector ˆ = sin q θ cos φ ˆ +sin x θ sin φ ˆ +cos y θ ˆ. z The depolarization contribution associated with the nonzero size of the component particles is represented by the dyadic term D + /lscript , which is conveniently given as a 3 × 3 subdyadic of the 6 × 6 dyadic D + /lscript per

$$\left [ \underline { D } ^ { + } _ { \ell } \right ] _ { m n } = \left [ \underline { \mathbb { D } } ^ { + } _ { \ell } \right ] _ { m n }, \quad ( m, n \in \{ 1, 2, 3 \} ).$$

Here the 6 × 6 dyadic [30]

$$\underline { \mathbb { D } } ^ { + } \ = \ \frac { \omega ^ { 4 } } { 4 \pi \mu _ { 0 } } \int _ { \phi = 0 } ^ { 2 \pi } \int _ { \theta = 0 } ^ { \pi } \frac { \sin \theta } { \left [ \left ( \underline { U } ^ { - 1 } \cdot \hat { \mathfrak { q } } \right ) \cdot \underline { \mathfrak { e } } _ { H C M } \cdot \left ( \underline { U } ^ { - 1 } \cdot \hat { \mathfrak { q } } \right ) \right ] \left ( \underline { U } ^ { - 1 } \cdot \hat { \mathfrak { q } } \right ) \cdot \left ( \underline { U } ^ { - 1 } \cdot \hat { \mathfrak { q } } \right ) } \\ & \times \left [ \frac { 1 } { \kappa _ { + } - \kappa _ { - } } \left ( \frac { \exp \left ( i \eta \right ) } { 2 q ^ { 2 } } \left ( 1 - i \eta \right ) \left \{ \det \left [ \underline { \Delta } \left ( \underline { U } ^ { - 1 } \cdot \mathfrak { q } \right ) \right ] \underline { \underline { \mathfrak { e } } ^ { + } \left ( \underline { U } ^ { - 1 } \cdot \mathfrak { q } \right ) } \right ] \underline { \underline { \mathfrak { e } } ^ { + } \left ( \underline { U } ^ { - 1 } \cdot \mathfrak { q } \right ) } \\ & \quad + \det \left [ \underline { \Delta } \left ( - \underline { U } ^ { - 1 } \cdot \mathfrak { q } \right ) \right ] \underline { \underline { \mathfrak { e } } ^ { + } \left ( - \underline { U } ^ { - 1 } \cdot \mathfrak { q } \right ) \right \} \right ) _ { q = \sqrt { \kappa _ { + } } } + \frac { \det \left [ \underline { \Delta } \left ( \theta \right ) \right ] } { \kappa _ { + } + \kappa _ { - } } \underline { \underline { \underline { \mathfrak { e } } ^ { + } \left ( \theta \right ) } } \right ) d \theta \, d \phi, \quad \left ( \ell = a, b \right ), \\ \text{with the vector} \, \mathfrak { q } = q \, \hat { q }, \text{ the scalars} \, \kappa _ { + } \text{ being the} \, \mathfrak { q } ^ { 2 } \text{ roots of} \, \det \left [ \underline { \Delta } \left ( \underline { U } ^ { - 1 } \cdot \mathfrak { q } \right ) \right ] = 0, \text{ and the} \, \mathfrak { \, } \mathfrak { \, } \mathfrak { \, } \mathfrak { \, } \mathfrak { \, } \mathfrak { \, } \mathfrak { \, } \mathfrak { \, } \mathfrak { \, } \mathfrak { \, } \mathfrak { \, } \mathfrak { \, } \mathfrak { \, } \mathfrak { \, } \mathfrak { \, } \mathfrak { \, } \mathfrak { \, } \mathfrak { \, } \mathfrak { \, } \mathfrak { \, } \mathfrak { \, } \mathfrak { \, } \mathfrak { \, } \mathfrak { \, } \mathfrak { \, } \mathfrak { \,                                                                                                                                                                                                         }$$

with the vector q = q ˆ, the scalars q κ ± being the q 2 roots of det [ A ( U -1 · q ) ] = 0, and the 6 × 6 dyadics

$$\underline { \mathbb { A } } ( \mathbf p ) = \left [ \begin{array} { c c } \underline { \epsilon } _ { H C M } & ( \mathbf p / \omega ) \times \underline { I } \\ - ( \mathbf p / \omega ) \times \underline { I } & \mu _ { 0 } \underline { I } \end{array} \right ]$$

and

$$\underline { \underline { G } } ^ { + } ( \mathbf p ) = \underline { \underline { \Delta } } ^ { - 1 } ( \mathbf p ) - \lim _ { | \mathbf p | \to \infty } \underline { \underline { \mathbf A } } ^ { - 1 } ( \mathbf p ).$$

Generally, numerical methods are needed to evaluate the integrals on the right sides of Eqs. (16) and (18), and also to extract the HCM permittivity dyadic /epsilon1 HCM from the nonlinear dyadic equation (13) [33].

## References

- [1] M. Born and E. Wolf, Principles of Optics , 6th Edition, Pergamon Press, Oxford, UK (1980).
- [2] W. Voigt, 'On the behaviour of pleochroitic crystals along directions in the neighbourhood of an optic axis,' Phil. Mag. 4 , 90-97 (1902).
- [3] S. Pancharatnam, 'The optical interference figures of amethystine quartz - Part II,' Proc. Ind. Acad. Sci. A 47 , 210-229 (1958).
- [4] F. I. Fedorov and A. M. Goncharenko, 'Propagation of light along the circular optical axes of absorbing crystals,' Opt. Spectrosc. (USSR) 14 , 51-53 (1963).
- [5] G. S. Ranganath, 'Optics of absorbing anisotropic media,' Curr. Sci. 67 , 231-237 (1994).
- [6] T. G. Mackay and A. Lakhtakia, 'Electromagnetic fields in linear bianisotropic mediums,' Prog. Optics 51 , 121-209 (2008).
- [7] A. Lakhtakia, 'Anomalous axial propagation in helicoidal bianisotropic media,' Opt. Commun. 157 , 193-201 (1998).
- [8] M. V. Berry, 'The optical singularities of bianisotropic crystals,' Proc. R. Soc. A 461 , 2071-2098 (2005).
- [9] T. G. Mackay and A. Lakhtakia, 'Voigt wave propagation in biaxial composite materials,' J. Opt. A: Pure Appl. Opt. 5 , 91-95 (2003).
- [10] T. G. Mackay and A. Lakhtakia, 'Correlation length facilitates Voigt wave propagation,' Waves Random Media 14 , L1-L11 (2004).
- [11] T. G. Mackay, 'On the sensitivity of directions that support Voigt wave propagation in infiltrated biaxial dielectric materials,' J. Nanophotonics 8 , 083993 (2014).

- [12] T. G. Mackay, 'On extended homogenization formalisms for nanocomposites,' J. Nanophotonics 2 , 021850 (2008).
- [13] L. Ward, The Optical Constants of Bulk Materials and Films , 2nd edition, Institute of Physics, Bristol, UK (2000).
- [14] T. G. Mackay and A. Lakhtakia, Electromagnetic Anisotropy and Bianisotropy: A Field Guide , Word Scientific, Singapore (2010).
- [15] G. B. Arfken and H. J. Weber, Mathematical Methods for Physicists , 4th Edition, Academic Press, London, UK (1995).
- [16] J. Gerardin and A. Lakhtakia, 'Conditions for Voigt wave propagation in linear, homogeneous, dielectric mediums,' Optik 112 , 493-495 (2001).
- [17] R. W. Boyd, Nonlinear Optics , 2nd edition, Academic Press, San Diego, CA, USA (2003).
- [18] A. Lakhtakia and T. G. Mackay, 'Electrical control of the linear optical properties of particulate composite materials,' Proc. R. Soc. A 463 , 583-592 (2007).
- [19] A. Yariv and P. Yeh, Photonics: Optical Electronics in Modern Communications , 6th edition, Oxford University Press, New York, NY, USA (2007).
- [20] A. Lakhtakia, 'Electrically tunable, ultranarrowband, circular-polarization rejection filters with electrooptic structurally chiral materials,' J. Eur. Opt. Soc. - Rapid Pubs. 1 , 06006 (2006).
- [21] A. V. Goncharenko, 'Generalizations of the Bruggeman equation and a concept of shape-distributed particle composites,' Phys. Rev. E 68 , 041108 (2003).
- [22] W. S. Weiglhofer, A. Lakhtakia and B. Michel, 'Maxwell Garnett and Bruggeman formalisms for a particulate composite with bianisotropic host medium,' Microwave Opt. Technol. Lett. 15 , 263-266 (1997); Erratum: 22 , 221 (1999).
- [23] M. Zgonik, R. Schlesser, I. Biaggio, E. Volt, J. Tscherry, and P. G¨nter, 'Material constants of KNbO u 3 relevant for electro- and acousto-optics,' J. Appl. Phys. 74 , 1287-1297 (1993).
- [24] T. G. Mackay and W. S. Weiglhofer, 'Homogenization of biaxial composite materials: dissipative anisotropic properties,' J. Opt. A: Pure Appl. Opt. 2 , 426-432 (2000).
- [25] T. G. Mackay and A. Lakhtakia, 'On scattering loss in electro-optic particulate composite materials,' J. Appl. Phys. 101 , 083523 (2007).
- [26] T. G. Mackay, 'Voigt waves in homogenized particulate composites based on isotropic dielectric components,' J. Opt. 13 , 105702 (2011).
- [27] M. J. Madou, Fundamentals of Microfabrication: The Science of Miniaturization (2nd ed) , CRC Press, Boca Raton, FL (2002).
- [28] A. I. Kingon, S. H. Rou, M. S. Ameen, T. M. Graettinger, K. Gifford, and O. Auciello, 'Deposition of electrooptic thin films,' Ceramic Transactions 14 , Electro-Optics and Non-linear Optic Materials, (Am. Cer. Soc., Westerville, OH), 1990, pp. 179-196.
- [29] T. S. El-Bawab, Optical Switching , Springer, New York, NY, USA (2006).
- [30] T. G. Mackay, 'Depolarization volume and correlation length in the homogenization of anisotropic dielectric composites,' Waves Random Media 14 , 485-498 (2004); Erratum: Waves Random Complex Media 16 , 85 (2006).

- [31] B. Michel, 'A Fourier space approach to the pointwise singularity of an anisotropic dielectric medium,' Int. J. Appl. Electromagn. Mech. 8 , 219-227 (1997).
- [32] B. Michel and W. S. Weiglhofer, 'Pointwise singularity of dyadic Green function in a general bianisotropic medium,' Arch. Elekron. Ubertrag. ¨ 51 , 219-223 (1997); Erratum: 52 , 31 (1998).
- [33] B. Michel, A. Lakhtakia and W. S. Weiglhofer, 'Homogenization of linear bianisotropic particulate composite media - Numerical studies,' Int. J. Appl. Electromag. Mech. 9 , 167-178 (1998); Erratum: 10 , 537-538 (1999).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0604v1-figure-ebf28d1b-a8f3-404a-a9a4-01f515723f27.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0604v1-figure-7d5ded96-abbf-4954-bff9-e3b90871707b.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0604v1-figure-6eeb5aa5-468d-4fd4-8281-4998fbb3fd06.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0604v1-figure-03472c0e-5489-4995-89c8-cfbcfd608a72.jpg)

Figure 1: The extended Bruggeman estimates of relative permittivity parameters of the HCM plotted versus E dc 3 ∈ -( 1 1) , × 10 8 Vm -1 and spheroid orientation angle ϕ ∈ (0 , 90) . The refractive index ◦ n a = 1 1, volume . fraction f a = 0 3, size parameter . η = 0 3 . /k 0 , and eccentricity parameter ρ = 3.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0604v1-figure-f3489178-6811-4849-b6ab-07bb2f369e45.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0604v1-figure-3edd64f2-8435-4d0b-a685-4280e6cf17b6.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0604v1-figure-c679b4be-d616-4c89-809e-a55e2bc250b9.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0604v1-figure-b8f9ff77-e91b-4f5e-b1ec-abf66afe831a.jpg)

Figure 2: The angles α 1 2 , and β 1 2 , , along with the quantities | W 1 2 , | , plotted versus E dc 3 ∈ ( -1 1) , × 10 8 Vm -1 for the spheroid orientation angles ϕ = 60 ◦ (red, solid curves), 30 ◦ (green, dashed curves), and 90 ◦ (blue, broken dashed curves). The refractive index n a = 1 1, volume fraction . f a = 0 3, size parameter . η = 0 3 . /k 0 , and eccentricity parameter ρ = 3.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0604v1-figure-5683e37e-190d-4202-bcba-c25b6d61b46d.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0604v1-figure-adaeedb3-d8e0-4126-b276-cd36f0fe398a.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0604v1-figure-be7c36c9-ba25-4569-9a02-9a5eb97a0964.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0604v1-figure-19181eb7-0ccd-49d4-b62f-f6cc950f3342.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0604v1-figure-78005bad-b753-43e2-ac5d-0f8aab7fb3d6.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0604v1-figure-d1067a13-fa62-4e95-ae4b-66ea39bd8e59.jpg)

Figure 3: The angles α 1 2 , and β 1 2 , plotted versus E dc 3 ∈ ( -1 1) , × 10 8 Vm -1 for the eccentricity parameter values ρ = 3 (red, solid curves), 2 (green, dashed curves), and 4 (blue, broken dashed curves). The refractive index n a = 1 1, volume fraction . f a = 0 3, size parameter . η = 0 3 . /k 0 , and spheroid orientation angle ϕ = 60 . ◦

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0604v1-figure-85d75f62-4d01-4f8b-88a6-db6b7c66f61e.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0604v1-figure-5a91fe9a-2167-4c7a-acb6-2b1709bf5288.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0604v1-figure-c057e92d-aeb6-413e-b573-ba63ed43a23b.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0604v1-figure-13fe7f19-ff53-4446-a56c-33bd79472b74.jpg)

Figure 4: The angles α 1 2 , and β 1 2 , plotted versus E dc 3 ∈ ( -1 1) , × 10 8 Vm -1 for the size parameter values η = 0 3 . /k 0 (red, solid curves), 0 2 . /k 0 (green, dashed curves), and 0 25 . /k 0 (blue, broken dashed curves). The refractive index n a = 1 1, volume fraction . f a = 0 3, spheroid orientation angle . ϕ = 60 , and eccentricity ◦ parameter ρ = 3.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0604v1-figure-a9983115-b800-41ee-94bb-b5052d3bb67f.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0604v1-figure-8975de80-375a-4303-85f2-4053592c0cf5.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0604v1-figure-095f4915-ba52-4886-a11c-49d140bd82e0.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0604v1-figure-8ce831fe-fada-4603-ada9-097d67460694.jpg)

Figure 5: The angles α 1 2 , and β 1 2 , plotted versus E dc 3 ∈ ( -1 1) , × 10 8 Vm -1 for the volume fractions f a = 0 3 . (red, solid curves), 0 35 (green, dashed curves), and 0 4 (blue, broken dashed curves). . . The refractive index n a = 1 1, size parameter . η = 0 3 . /k 0 , spheroid orientation angle ϕ = 60 , and eccentricity parameter ◦ ρ = 3.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0604v1-figure-7c2b2560-0e5b-4c3c-8457-df2e4940c24c.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0604v1-figure-262882ba-911c-4b95-8a21-e87eee25e53f.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0604v1-figure-75ed32d1-0cd5-446a-a1ed-6c00d48f27ff.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0604v1-figure-9879e9fc-793b-456c-8d05-1facc4e1b8c0.jpg)

Figure 6: The angles α 1 2 , and β 1 2 , plotted versus E dc 3 ∈ ( -1 1) , × 10 8 Vm -1 for the refractive index values n a = 1 1 (red, solid curves), 1 05 (green, dashed curves), and 1 3 (blue, broken dashed curves). The volume . . . fraction f a = 0 3, size parameter . η = 0 3 . /k 0 , eccentricity parameter ρ = 3, and spheroid orientation angle ϕ = 60 . ◦

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0604v1-figure-15a32265-33e8-4991-8d2f-ad2043502310.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0604v1-figure-e0efb6ca-6548-45fe-b0b3-a77cb9854513.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0604v1-figure-34392c30-8314-4a2a-92f6-82a3af933250.jpg)
Tom G. Mackay
Controlling Voigt waves by the Pockels effect
Voigt wave propagation was investigated in a homogenized composite material
(HCM) arising from a porous electro--optic host material infiltrated by a fluid
of refractive index $n_a$. The constitutive parameters of the HCM were
estimated using the extended Bruggeman homogenization formalism. Numerical
studies revealed that the directions which support Voigt wave propagation in
the HCM could be substantially controlled by means of an applied dc electric
field. Furthermore, the extent to which this control could be achieved was
found to be sensitive to the porosity of the host material, the shapes, sizes
and orientations of the pores, as well as the refractive index $n_a$. These
findings may be particularly significant for potential technological
applications of Voigt waves, such as in optical sensing.