Spontaneous natural optical activity in disordered media

Spontaneous natural optical activity in disordered media

F. A. Pinheiro, V. A. Fedotov, N. Papasimakis, N. I. Zheludev Instituto de Física, Universidade Federal do Rio de Janeiro, Rio de Janeiro-RJ, 21941-972, Brazil Optoelectronics Research Centre and Centre for Photonic Metamaterials, University of Southampton, United Kingdom Centre for Disruptive Photonic Technologies, School of Physical and Mathematical Sciences and The Photonics Institute, Nanyang Technological University, Singapore

We demonstrate natural optical activity in disordered ensembles of non-chiral plasmonic resonators. We show that the statistical distributions of rotatory power and spatial dichroism are strongly dependent on the scattering mean free path in diffusive random media. This result is explained in terms of the intrinsic geometric chirality of disordered media, as they lack mirror symmetry. We argue that chirality and natural optical activity of disordered systems can be quantified by the standard deviation of both rotatory power and spatial dichroism. Our results are based on microscopic electromagnetic wave transport theory coupled to vectorial Green’s matrix method for pointlike scatterers, and are independently confirmed by full-wave simulations.

The concept of chirality, introduced by Lord Kelvin to designate any geometrical object or ensemble of points that lack mirror symmetry, pervades the natural world. The DNA double-helix structure and the related long-standing enigma of the origin of homochirality of life are just some examples of this fundamental concept in science wagniere (). Since the pioneer work by Pasteur on the resolution of tartaric acid, the development of efficient synthesis processes of naturally occurring chiral molecules is of utmost importance in chemistry, biology, and medicine since chiral systems have widespread applications in pharmaceuticals, asymmetric synthesis, and catalysis wagniere (). In addition to natural chiral systems, designing and characterizing artificial chiral structures constitute an intense research activity in recent years guerrero (). Chiral metamaterials are striking examples of artificial media without mirror symmetry rogacheva (); they have been applied to generate novel optical properties such as negative index media pendry (); plum2009 (), broadband circular polarizers gansel (), and enantioselective photochemistry tang2011 (). Another possibility of artificially generating chiral systems is to assembly achiral plasmonic nanoparticles in chiral geometries, such as helices, pyramids, and twisted nanorod pairs zheludev2003 (); zheludev2006 (); valev2013 (); kuzuk (); shen2013 (); hentschel (); zan (); ferry (). Remarkably, natural optical activity in these systems can be considerably larger than in molecular ones. On a more fundamental side, the conditions for observing natural optical activity, involving the preservation of the helicity of light, have been recently revisited corbaton1 (); corbaton2 (); corbaton3 ().

Despite all efforts to understand the optical properties of naturally occurring chiral media and to design artificial ones, there is an important yet overlooked class of chiral systems: disordered structures. Such systems have neither center nor plane of symmetry, so that they should exhibit natural optical activity as any other chiral medium. Interestingly, there is some previous experimental evidence of natural optical activity in random media, but this effect has never been attributed to the presence of disorder hadley (); silverman (); drachev (); instead, alternative explanations, such as surface contamination by unwanted chiral substances, have been put forward hadley (); silverman ().

Figure 1: Examples of chiral ensembles: dipole scatterers arranged in a helix (a) and in a random configuration (b); (c) Schematic representation of polarisation rotation in a disordered medium.

Here, we demonstrate natural optical activity due to intrinsic geometric chirality in disordered, diffusive scattering systems. By means of microscopic electromagnetic wave transport theory, we derive an expression for the rotatory power and spatial dichroism of a medium consisting of randomly distributed pointlike scatterers (see Fig. 1). We perform a statistical analysis of natural optical activity in random media to conclude that the standard deviation of both rotatory power and spatial dichroism are strongly dependent on the density of scatterers and the scattering mean free path. Finally, we show that natural optical activity in disordered ensembles can be more that an order of magnitude higher than in the archetypical helical configurations zan ().

In the following we describe the derivation of the expression for natural optical activity in disordered media using the well established microscopic theory of vector wave diffusion pingboek (), which finds its equivalent in electronic transport theory mahan (). The second-rank tensor represents the specific intensity of the diffuse radiation in the direction . It can be expressed as mahan ():


with the Dyson Green’s tensor. For a low density of the scatterers the (second rank) mass-operator is related to the T-matrix of one independent scatterer pingboek (),


where we have implicitly applied rotational averaging of Eq. (2). In Eq. (1) the second-rank tensor describes the angular dependence of the diffuse energy flow. It obeys the integral equation mahan ()


where the current tensor is ( being the unit matrix), and the (four-rank) irreducible vertex is related, after rotational averaging, to the T-matrix of one independent scatterer according to pingboek ()


with . For isotropic, achiral media subject to weak disorder, the well known Boltzmann results applies, with .

In isotropic chiral media, the tensor (neglecting longitudinal terms) takes the form josa ():


where is the projector upon the space of transverse polarization (normal to ), is the antisymmetric Hermitian tensor ( being the Levi-Civita tensor). In chiral media obeys the symmetry relation .

Our objective is to solve Eq. (3) and to explicitly obtain the coefficients , in Eq. (5), and a expression for natural optical activity. To this end, let us separate the mass operator in a symmetric (S) and an antisymmetric (A) part:


where is given by Eq. (2). Equation (6) defines the pseudoscalar , whose real and imaginary parts determines the rotatory power and the spatial dichroism of the effective medium, respectively. Inserting Eq. (5) into Eq. (3) and performing the operations and results in the two following coupled equations for and :


where we have put and used the fact that . It can be shown that , , and can be explicitly written in terms of the (off-shell) scattering total T-matrix of an ensemble of point dipoles josa (), which describes multiple scattering from the direction to


with a symmetric matrix describing scattering of light from particle to that is obtained by diagonalization dipolemodel (). As a result, multiple scattering is treated exactly, i.e. all scattering orders are included; the only approximation within this model is at the level of single scattering, as scatterers are treated as point dipoles. Hence the knowledge of the scattering T-matrix for the ensemble of pointlike dipoles enables one to calculate , , , and, consequently, and  josa (). In terms of and after averaging over all solid rotations, the final expression for the pseudoscalar , which describes natural optical activity, is josa ()


where is the spherical Bessel function of the first kind, and is the relative position between scatterers and . We have verified that Eq. (10) is indeed a pseudoscalar, i.e. it changes sign by a mirror operation. We have also verified that for , identically vanishes, reflecting the fact that a chiral system must be composed of at least 4 particles.

In Fig. 2 the distributions of , which gives the optical rotatory power, is calculated using Eq. (10) for 1000 different disorder realizations and two different particle densities. The scatterers are randomly distributed inside a cylinder of fixed volume (height and radius , where is the incident wavelength). In Fig. 2(a) is calculated for scatterers, which corresponds to (with the wavenumber and the mean free path). In Fig. 2(b) is calculated for scatterers, in which case multiple light scattering is stronger (). For both densities light transport in the random medium is diffusive. The values of have been normalized by the ones corresponding to scatterers distributed along a helix inside a cylinder of equal volume (see Fig. 1a). Hence Fig. 2b reveals that for a given configuration of random scatterers the value of rotatory power can be up to 60 times larger than for a system with the same density, where the particles are distributed along a helix, which is the hallmark of a chiral system. Figure 2 also shows that , which can be explained by the fact that for large number of disorder realizations the mirror image of any configuration is equally probable. However, Fig. 2 demonstrates that the standard deviation of the distribution of is strongly dependent on , which is a measure of how strong light scattering is inside the medium. We have verified that the standard deviation of exhibits a similar dependence on .

Figure 2: Histogram of that gives the rotatory power in disordered medium for (a) 10 and (b) 700 pointlike scatterers randomly distributed inside a cylinder of fixed volume (height and radius ).

Figure 3 shows the real () and imaginary () parts of the standard deviation of for 1000 disorder realizations as a function of for a diffusive scattering system () composed of resonant pointlike scatterers randomly distributed inside a cylinder of fixed volume (height and radius ). The results are normalized by the value of for a system composed of 100 scatterers distributed along a helix, a system that is clearly chiral, contained in a cylinder of equal volume. This allows us to obtain the order-of-magnitude of natural optical activity of a disordered medium, as natural optical activity for nanoparticles oriented along a helix has been calculated in Ref. zan (). Figure 3 reveals that monotonically increases as decreases, showing that the magnitude of natural optical activity increases as the density of scatterers increases. More importantly, for can be approximately 15 times larger than the case of a helix. It is also important to emphasize that both and are small for ballistic and weakly scattering systems, for which . This demonstrates that multiple light scattering is necessary in order to produce significant values of natural optical activity. In particular, we have verified that at least four scattering events are required to generate a non-vanishing value of , so that light can “probe” the chiral configuration associated to random scatterers. Furthermore, in the diffusive regime, in which the system size is much larger than , the statistical distribution of is not expected to depend on the system size for a large set of disorder realizations.

To further investigate the interplay between the chirality associated with the random configuration of scatterers and natural optical activity, we calculate a geometrical chiral index defined in Ref. harris () according to


with the Levi-Civita tensor. In Eq. (11) the tensors , , and are given by


where is the position vector of the scatterers relative to the system center of mass, and and are the eigenvalues corresponding to the matrices and , respectively. The chiral index , which only depends on the scatterers positions, is a pseudoscalar invariant under rotations that vanishes for achiral configurations harris (). The choice of chiral index is not unique and other chiral measures do exist chiralindexes (). The chiral index (11) has been chosen here by its simplicity and, more importantly, because it has been shown to be related to physical observables, such as the magneto-chiral scattering cross-section pinheiro2002 (), the electronic current in chiral carbon nanotubes pinheiro2010 (), and the pitch of a cholesteric liquid crystal harris (). For a random ensemble of scatterers the average value of over many disorder realizations is zero since the mirror image of any configuration is equally probable. However, the standard deviation is strongly dependent on the particle density. This can be seen in Fig. 3, where is calculated using Eq. (11) for 1000 disorder realizations. Remarkably, Fig. 3 shows that is proportional to , confirming that the calculated optical activity is indeed related to the intrinsic chiral nature of disordered systems. Together with previous experimental evidence of optical activity in disordered systems in the absence of any chiral substance hadley (); silverman (); drachev (), these results strongly suggest that large optical activity should be observed in a disordered medium due to its intrinsic chirality. As a possible application, Fig. 3 suggests that one could determine the particle density in solution by measuring natural optical activity.

Figure 3: The standard deviation of the distributions of and for 1000 configurations of randomly distributed identical pointlike scatterers inside a inside a cylinder of fixed volume (height and radius ) as function of , with the wavenumber and the scattering mean free path. The values of and have been normalized by and , which correspond to rotatory power and spatial dichroism of 100 pointlike dipoles distributed along a helix inside a cylinder of the same volume. Also shown is the standard deviation of the distributions of chiral geometrical parameter , defined in Ref. harris (), of the corresponding 1000 random configurations. The insets show a scheme of characteristic random configurations for high particle density (small ) and for low particle density (large ).

In order to investigate whether the calculated natural optical activity is related to the anisotropy of the disordered medium, we have calculated the real and imaginary parts of for a single, fixed configuration of random scatterers inside a cylinder with the same dimensions of Fig. 2, which is continuously rotated by an angle around the main axis of the cylinder. The result, shown in Fig. 4, reveals that does not depend on . Hence we rule out the possibility that the reported results are due to the anisotropy of the random medium as Fig. 4 proves. Indeed, polarization change due to anisotropy does not require multiple light scattering. On the contrary, here optical activity is related to the intrinsic chirality of random media, which light can only probe when multiple scattering occurs (or, more precisely, when at least four scattering events occurs), as Fig. 3 demonstrates.

Figure 4: and for a fixed random configuration of 100 dipoles inside a cylinder of same volume as a function of the rotation angle around the cylinder main axis. The values of and have been normalized by and that correspond to the corresponding values of for 100 scatterers distributed along a helix inside a cylinder of the same volume.
Figure 5: Circular dichroism (CD) spectra of random and helical Au nanoparticle configurations. Blue line represents the CD of a helical arrangement with radius nm & pitch nm. The average, , over 310 realizations of random nanoparticle configurations is represented by the red line. Magenta line represents the standard deviation of CD. The radius of the nanoparticles is considered to be 25 nm. All CD spectra have been obtained by finite element modeling.

To further confirm the validity of the results from the microscopic theory of light transport, we performed finite element full-wave simulations of random and helical configurations of achiral Au nanoparticles. In Fig. 5, we show circular dichroism (CD) as a function of wavelength for random (red and magenta lines) and helical (blue lines) ensembles of 25 nanoparticles. Both random and helical configurations occupy the same volume at a density of nanoparticles/. The helical configuration exhibits strong circular dichroism in the vicinity of the plasmonic nanoparticle resonance ( nm). On the other hand, in the case of random nanoparticle ensembles, the average over 310 realizations is close to zero, as predicted by the previous microscopic analysis. However, the corresponding standard deviation () exceeds the CD of the helical configuration by a factor of 2 at resonance, while for specific realisations this factor can be as high as 10 (not shown here). These results are in qualitative agreement with the predictions of the point-dipole model for diffusive scattering systems, and demonstrate unambiguously that polarization effects in disordered scattering systems can be much stronger than in a typical chiral medium (helix) due to its intrinsic chirality.

In conclusion, we have investigated natural optical activity of diffusive disordered systems composed of non-chiral plasmonic resonators. Using microscopic electromagnetic wave transport theory coupled to vectorial Green’s matrix model for pointlike scatterers, we have demonstrated that the distributions of the rotatory power and spatial dichroism in diffusive disordered media is strongly dependent on the scattering mean free path. In particular, despite the fact the mean value of rotatory power and spatial dichroism in random media vanishes after averaging over many disorder realizations, their standard deviation exhibits a strong monotonic increase as the scattering mean free path decreases. As a result we argue that the standard deviation of both rotatory power and spatial dichroism is the appropriate quantity to probe chirality and natural optical activity in disordered systems of non-chiral plasmonic oscillators. This finding is independently corroborated by full-wave electromagnetic calculations. We have found that the standard deviation of a purely geometrical chiral parameter, defined in Ref. harris () and that only depends on the scatterers’ positions, is proportional to the standard deviation of natural optical activity in random media. This result corroborates that the latter is the appropriate quantity to probe natural optical in disordered media, which is intrinsically related to the fact that random systems lack mirror symmetry.

We thank B. A. van Tiggelen, L. Dal Negro, and O. Buchnev for fruitful discussions. We thank The Royal Society-Newton Advanced Fellowship (Grant no. NA150208). F.A.P. thanks the Optoelectronics Research Centre and Centre for Photonic Metamaterials, University of Southampton, for the hospitality, and CAPES for funding his visit (Grant No. BEX 1497/14-6). F.A.P. also acknowledges CNPq (Grant No. 303286/2013-0) for financial support. N. P. acknowledges the support of the UK’s Engineering and Physical Sciences Research Council (Grant no. EP/M008797/1).


  • (1) G. H. Wagniére, On Chirality and the Universal Asymmetry: Reflections on Image and Mirror Image, Wiley-VCH, Zurich, 2007.
  • (2) A. Guerrero-Martinez et al., Nanotoday 6, 381 (2011).
  • (3) A.V. Rogacheva, V.A. Fedotov, A.S. Schwanecke, and N.I. Zheludev, Phys. Rev. Lett. 97, 177401 (2006).
  • (4) B. Pendry, Science 306, 1353 (2004).
  • (5) E. Plum et al., Phys. Rev. B 79, 035407 (2009).
  • (6) J.K. Gansel et al., Science 325,1513 (2009).
  • (7) Y. Tang and A. Cohen, Science 332, 333 (2011).
  • (8) A. Papakostas et al., Phys. Rev. Lett. 90, 107404 (2003).
  • (9) S. Boruhovich, S. L. Prosvirnin, A. S. Schwanecke, and N. I. Zheludev., Proc. Eur. Microw. Assoc. 2, 89 (2006).
  • (10) V.K. Valev, J.J. Baumberg, C. Sibilia, and T. Verbiest, Adv. Mater. 25, 2517 (2013).
  • (11) A. Kuzyk et al., Nature 483, 311 (2012).
  • (12) X. Shen et al., Nano Lett. 13, 2128 (2013).
  • (13) M. Hentschel et al. Nano Lett. 12, 2542 (2012).
  • (14) Z. Fan and A.O. Govorov, Nano Lett. 10, 2580 (2010).
  • (15) V.E. Ferry, J.M. Smith, and A. P. Alivisatos, ACS Photonics 1, 1189 (2014).
  • (16) I. Fernandez-Corbaton, X. Vidal, N. Tischler, G. Molina-Terriza, J. Chem. Phys. 38, 21431 (2013).
  • (17) N. Tischler et al., Light: Science Applications 3, e183 (2014).
  • (18) I. Fernandez-Corbaton, M. Fruhnert, and C. Rockstuhl, ACS Photonics 2, 376 (2015).
  • (19) X. Guo, M. F. G. Wood, and I. Alex Vitkin, Applied Optics 46, 4491 (2007).
  • (20) K. C. Hadley, I.A. Vitkin, J. Biomed. Opt. 7, 291 (2002).
  • (21) M. E Silverman, W. Strange, J. Badoz, I.A. Vitkin, Opt. Comm. 132, 410 (1996).
  • (22) V. P. Drachev et al., J. Opt. Soc. Am. B 18, 1896 (2001).
  • (23) G. D. Mahan, Many-Particle Physics (Plenum, New York, 1981).
  • (24) P. Sheng, Wave Scattering, Localization and Mesoscopic Phenomena (Academic, San Diego, 1995).
  • (25) F. A. Pinheiro and B. A. van Tiggelen, J. Opt. Soc. Am. A 20, 99 (2003).
  • (26) F. A. Pinheiro, M. Rusek, A. Orlowski, and B. A. van Tiggelen, Phys. Rev. E 69, 026605 (2004); F. A. Pinheiro, Phys. Rev. A 78, 023812 (2008); M. Rusek, J. Mostowski, and A. Orlowski, Phys. Rev. A 61, 022704 (2000); A. Goetschy and S.E. Skipetrov, Phys. Rev E 84, 011150 (2011).
  • (27) A.B. Harris, R. D. Kamien, and T. C. Lubensky, Rev. Mod. Phys. 71, 1745 (1999).
  • (28) A Potts, D. M. Bagnall, and N. I. Zheludev, J. Opt. A 6, 193 (2004); G. Gilat, J. Phys. A 22, L545 (1989); M.A. Osipov, B.T. Pickup, and D.A. Dunmur, Mol. Phys. 84,1193 (1995); M. P. Neal, M. Solymosi, M. R. Wilson, and D. J. Earl, The Journal of Chem. Phys. 119, 3567 (2003).
  • (29) F. A. Pinheiro and B. A. van Tiggelen, Phys. Rev. E 66, 016607 (2002).
  • (30) F. A. Pinheiro, S. J. S. da Silva, E. R. Granhen, and J. Del Nero, Phys. Rev. E 81, 115456 (2010).
  • (31) X. Duan, S. Kamin, F. Sterl, H. Giessen, and N. Liu, Nano Lett. 16 1462 (2016).
  • (32) M. Decker, M. W. Klein, M. Wegener, and S. Linden, Opt. Express 32, 856 (2007).
  • (33) B. P. Ablitt et al., Waves in Random Media 9, 561 (1999).
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description