1. Introduction
It is usually a simple task to tell by eye whether an object is chiral or not: Achiral objects are superimposable onto their mirror image and, accordingly, they possess a mirror plane [
1]. Recently, chiral scatterers have gained significant interest in nano-optics due to their potential to enhance the weak optical signal of chiral molecules [
2,
3,
4]. Especially, the quantities of optical chirality and optical helicity as well as their relation to duality symmetry are subjects of current research [
5]. The most established experimental technique in this field is the analysis of the circular dichroism (CD) spectrum which equals the differential energy extinction due to the illumination by right- and left-handed circularly polarized light [
6].
In order to observe such chiral electromagnetic response, it seems to be obvious that geometrically chiral scatterers are required. However, it has been shown that extrinsic chirality, that is, a chiral configuration of the illumination and geometric parameters, yields comparable effects as intrinsically chiral objects [
7]. By tuning the far-field polarization of the illumination, large chiral near-fields may even be generated in the viscinity of achiral objects [
8]. In CD measurements, randomly orientied molecules are investigated, which can be classified by their
T-matrix [
9]. The latter has been used for quantifying the electromagnetic (e.m.) chirality, based on a novel definition of it [
10].
However, the quantification of the geometric chirality is an elusive task [
11] and even the unambiguous association of the terms right- and left-handed enantiomer of a general object is impossible [
12]. Different coefficients attempting to rate the chirality of an object are based on the maximal overlap of two mirror images [
13] as well as the Hausdorff distance [
14]. The choice of a specific coefficient determines the most chiral object [
15], that is, there is no natural choice for quantifying geometric chirality. This also holds for the various figures of merit estimating the e.m. chirality. Similar correlations between geometric and optical properties are investigated with respect to the non-sphericity of arbitrary scatterers [
16].
In this study, we start by transferring the simple procedure of finding a mirror plane to optics. Such symmetries are present in different mathematical descriptions as for example, block structures in the Mueller scattering matrix [
17]. Here, we analyse the
T-matrix and its associated geometric mirror symmetries by employing translation and rotation theorems of vector spherial harmonics. We illustrate this concept with numerical simulations of an experimentally realized gold helix. Different quantifications of the e.m. chirality are compared. Furthermore, the symmetry planes found in the optical response by our method are correlated to those of geometric origin. It is shown that the complex optical response, including higher order multipoles, yields mirror planes in the
T-matrix which are not directly related to geometric symmetries.
In the following, we would like to briefly introduce the theory behind the methods described in this study. Further information may be found in
Supplementary Materials.
The most general description of an isolated optical scatterer is the well-known
T-matrix [
18]. It relates an arbitrary incident field with the scattered field caused by the scattering object. The optical response to
any incident field is included in the
T-matrix. Accordingly, the following analysis of
T is independent of
specific illumination parameters such as the direction, polarization and beam shape. The goal of this study is to obtain insights into illumination-
independent symmetries of the scatterer.
Usually, both the incident as well as the scattered field are given in the basis of vector spherical harmonics for computations with the
T-matrix [
19] (see
Supplementary Materials). Physically observable quantities such as the scattered energy, the absorption, as well as the flux of optical chirality are readily computed from
T [
9]. In numerical simulations,
T may be computed with high accuracy [
20]. Knowing the response of the left-handed object
enables the analytic computation of the response of its mirror image
:
where we choose mirroring in the
-plane
without loss of generality (see
Supplementary Materials for further details on notation). Note that the terminology of right-
and left-handed
is ambiguous, as pointed out before, and may be interchanged.
Since we aim to investigate arbitrary mirror planes, we note that an arbitrary plane is given by the three spherical coordinates of its normal: the inclination
and the azimuthal angle
, as well as the distance
d from the origin. We define the according transformation
acting on the object as
where
is the translation of the
T-matrix in the direction given by the angles and the distance and
and
are the rotations around the
z- and
y-axis, respectively [
21] (see
Supplementary Materials).
For a geometrically achiral object (see
Figure 1a) there exists at least one transformation
such that
. On the other hand, the lack of a geometric mirror plane of a chiral object
Figure 1b implies that there exists no such transformation and that
and
do not coincide for any set of transformation parameters
. Note that this does not generally hold in the long wavelength limit, that is, the incident wavelength being much larger than the dimension of the scatterer, due to chiral dispersion [
22].
2. Results
For investigating the role of the geometric shape in nano-optics, it is of interest to identify those planes of highest symmetry of a chiral object: Although there is no mirror plane in a chiral object, a transformation may be identified in which the right- and left-handed
T-matrices are closest to one another. Rating the closeness is done here by calculating the 2-norm of the difference of these two matrices. Accordingly, we introduce the coefficient
which minimizes the difference between the
T-matrices of mirror images as
This means that for the mirror plane corresponding to minimal parameters
of (
3), the optical responses of the two mirror images are as similar as possible. In other words, the mirror images are hardly distinguishable. For an achiral object
vanishes since there exists a transformation for which the mirror images are identical.
Obviously, the choice of the norm is not unique and other quantifications of similarity of the mirror images could be defined (see
Supplementary Materials for the physical relevance of the 2-norm). A recently introduced coefficient
is, for example, based on the singular-value decomposition of the
T-matrix in the helicity basis [
10]. Alternatively, the angular-averaged differential energy extinction
due to illuminating with either right- or left-handed circularly polarized plane waves is experimentally accessible as the CD spectrum.
In order to exemplary introduce our formalism and compare it to previous work, we investigate a nano-optical device numerically. The finite element method is employed to accurately simulate the electromagnetic properties due to incident monochromatic light. Within this study we use the commercial FEM package JCMsuite [
23]. In postprocessing, the
T-matrix is computed by decomposing the scattered field into vector spherical wave functions [
20] from illumination with 150 plane waves with randomly chosen parameters (see
Supplementary Materials).
In
Figure 2, we compare simulations of the aforementioned three coefficients quantifying the e.m. chirality for a gold helix as realized experimentally [
24]. The helix is constructed on the surface of a cylinder with height 230nm and radius 60nm (see
Supplementary Materials). The CD spectrum
shows zero values at incident wavelengths of
and
. If only these wavelengths were analyzed, one could draw the conclusion that an achiral object is investigated. This contradicts the goal of this study to obtain insights into illumination-
independent symmetries of the scatterer—for illuminations with
and
, the scatterer seems to be geometrically achiral which is obviously not the case. Nevertheless, CD makes the chiral geometric nature of the helix visible as a maximum at
and a minimum at
of smaller amplitude. For a helix with an opposite twist—that is, the mirror image—the roles of the extrema are interchanged.
On the other hand, the coefficient
is normalized by the average interaction strength of the
T-matrix at each wavelength. This yields a fairly flat spectrum with two narrow minima below
at the two
for which
. These minima are not present in the minimized
introduced in (
3). However, the maxima of this latter coefficient are in accordance with the experimentally observable CD extrema (
). In the long wavelength regime, all three coefficients tend to zero as expected for point-like particles due to vanishing off-diagonal elements in the
T-matrix.
The minimization in the three-dimensional parameter space in (
3) is carried out using Bayesian optimization [
25] (see
Supplementary Materials). Since the shape of the minimized function highly depends on the actual object, the Bayesian approach is well suited for finding a global minimum. The parameters
of the optimized value are related to geometric mirror planes. In
Figure 3a, the planes following from the respective transformation
of the
-plane are plotted for all incident wavelengths from
to
. The inclination
and azimuthal angle
are given in the shown coordinate system which is centered at the centroid of the helix.
We identify three distinct classes shown in blue, red and green. These correspond to planes which are parallel and perpendicular to the helix axis, as well as tilted by a small angle
from the horizontal position, respectively. The dark grey plane corresponds to the minimal geometric parameters which will be explained in the following paragraphs. Details on the optimization such as challenging flat behaviour for translations from the centroid and, on the obtained minimizing parameters, are given in
Supplementary Materials. Note that the minimization required to obtain the illumination-independent coefficient
involves significantly higher numerical effort than the simple averaging for
for which most information contained in
T is ignored.
Next, we compare the findings on the symmetry based on the optical
T-matrix to those stemming from purely geometric properties. As discussed previously, there is no coefficient which unambiguously rates the geometric chirality of an object. We choose a coefficient
based on the overlap of the left-
and right-handed
object, where the latter results from mirroring
at the
-plane and transformation with
. Namely, the volume
V of the overlap is compared to the volume of the object [
13] (see
Supplementary Materials):
This coefficient vanishes for achiral objects as required for a degree of chirality [
14].
Figure 3b displays the geometric chirality coefficient
for planes rotated around the centroid of the helix as a grey colourmap. Dark regions with large values of
indicate a vanishing overlap between the two mirrored helices. Note that for large distances to the origin
, the mirror images do not overlap and
. However, this is always possible no matter if the object is chiral or not. As in the case of
, the parameter points of interest of
are those corresponding to a minimum: The minimum
in
Figure 3b occurs at
and
which show the intrinsic chiral property of the investigated helix. These two minima are equivalent since a finite helix is
symmetric. The corresponding transformed
-plane is shown in dark grey in
Figure 3a.
Alongside the geometric coefficient
, the planes identified for the minimized
T-matrix difference are shown as colored circles in
Figure 3b. The colors (red, blue and green) of these circles are the same colors used for the planes, that is, a direct comparison of the angle parameters is possible. As seen, the planes are ranked according to their
values: The perpendicular class 1 (blue) has
The flat planes belonging to class 2 (red) show
and
and the tilted class 3 (green) has
and
.
3. Discussion
None of the three optical symmetry planes is directly related to the geometric mirror plane of the helix. However,
Figure 3b enables the comparison of geometric and optical symmetries. In order to further analyze the optical response, we show the wavelength-dependent classification of the symmetry planes on top of
Figure 4. The three classes correspond to sharply separated wavelength ranges: Class 1 is valid for
. For larger wavelengths
, the
T-matrix possesses the symmetry according to planes of class 2. Finally, in the long wavelength regime (
), the symmetry is in class 3.
The analysis in
Figure 3b suggests that class 3 (green) is the closest one to the geometric mirror plane. This is further strengthened by the full angular spectrum of the optical chirality coefficient
(see
Supplementary Materials). Accordingly, we find that the optical response is dominated by the geometric shape for long wavelengths. Obviously, the optics is dominated by the electric dipole moment in this regime which is also shown in
Figure 4. Here, the mean of the diagonal entries of submatrices of the
T-matrix are shown. These are proportional to the electric and magnetic dipole moments as well as to the electric quadrupole moments.
The three symmetry classes of the
T-matrix occur close to three electric dipole peaks (
= 623, 833, and 1473 nm) and are influenced by the anisotropy of the
T-matrix. Truly chiral behaviour, as observed here, however, originates not from anisotropy but from coupling between electric and magnetic multipoles [
26]. In
Supplementary Materials, we elaborate on the complex interplay between these different contributions in the dipolar limit. Here, we limit the discussion to the main aspects of different multipolar contributions.
For large wavelengths with symmetry of class 3, the electric dipoles are much larger then any other induced multipole. In the intermediated regime of symmetry class 2, the magnetic dipole moment significantly increases. For short wavelengths with planes of class 1, the electric quadrupole moment is stronger than the magnetic dipole moment which yields the change in the optical symmetry. Higher order multipoles including mixed electric-magnetic moments are depicted in
Supplementary Materials, in which it is shown that the dominant dipolar moments contribute additionally to the variation of mirror planes. This elaborated study of multipolar resonances underlines again that the chiral response deviates from expectations due to a purely geometrical analysis of the scatterer.