Abstract
Noniridescent structural colors based on disordered arrangement of monodisperse spherical particles, also called photonic glass, show low color saturation due to gradual transition in the reflectivity spectrum. No significant improvement is usually expected from particles optimization, as Mie resonances are broad for small dielectric particles with moderate refractive index. Moreover, the short range order of a photonic glass alone is also insufficient to cause sharp spectral features. We show here, that the combination of a wellchosen particle geometry with the short range order of a photonic glass has strong synergetic effects. Using a firstorder approximation and an Ewald sphere construction the reflectivity of such structures can be related to the Fourier transform of the permittivity distribution. The Fourier transform required for a highly saturated color can be achieved by tailoring the substructure of the motif. We show that this can be obtained by choosing coreshell particles with a nonmonotonous refractive index distribution from the center of the particle through the shell and into the background material. The firstorder theoretical predictions are confirmed by numerical simulations.
Introduction
Structural color is a color based on selective light scattering and reflection from nanostructures^{1,2,3}. The commercial pigment based color derives from light absorption by electron transitions and is dependent on the presence of a defined chemical structure, which can be altered by UV radiation during later use or high temperature processing during manufacturing^{2,4}. Also, some of the pigments contain toxic materials that can be harmful in production or disposal, initiating the need for alternatives^{2,5,6}. At the same time, structural colors depend on the refractive index distribution, only and thus can be produced from environmentally friendly materials such as silica, alumina, zirconia etc. and therefore bear the potential of high UV and temperature stability. Structural colors can be divided into two classes: iridescent and noniridescent colors. An iridescent color is usually based on periodical structures with the periodical length in the order of visible light, known as photonic crystals (PhCs)^{7}. The noniridescent structural colors are angleindependent which means the color impression is the same for different illumination and observation angles. Historically, many research groups focused on microstructures mimicking biological structures to achieve noniridescent colors. For example, the feathers of many birds can exhibit bright noniridescent structural colors^{8,9,10,11}. Some birds’ feathers have structures similar to random compact arrangement of spherical particles. Such disordered arrangements, also called photonic glass (PhG) in contrast to PhC, can be obtained by selfassembly of monodispersed spherical particles^{12,13}. Recently PhGs have attracted a lot of attention in the field of noniridescent structural colors^{14,15,16,17,18}. Noniridescent structural colors produced by amorphous structures are mostly short wavelength colors such as violet or blue^{19}, since the typical band of scattered wavelengths is situated at the edge of the human eye sensitivity range which leads to the impression of a pure blue or violet as even shorter scattered wavelengths do not contribute to a color mix. Longerwavelength structural colors towards red are difficult to obtain as PhG structures, in addition, always have significant scattering in short wavelength range. These wavelengths will mix with the intended red which spoils the color impression. Still, longer wavelength colors can be produced by introducing a broadband absorber such as carbon black or others^{20,21,22,23,24}. These will take out particularly the shorter multiply scattered wavelengths (blue or green colors), which spent more time (or longer paths) inside the disordered medium and, therefore, experience a larger probability of absorption. The PhG possesses a shortrange order and the Fourier transform (FT) of its permittivity distribution is a spherical shell. The bright noniridescent structural colors directly correlate to the spherical shell shape of the FT such as observed for feathers of the male Plumthroated Cotinga^{11}.
However, the reported transition from the noreflection to backreflection regime is still rather smooth resulting in low color saturation. Some experimental effort was invested in PhG based structural colors with coreshell particles leading to marginal spectral improvement only^{14,25,26,27,28}, which is due to the lack of the theoretical understanding of the influence of coreshell geometry on color purity. Most of the explanations so far are based on a manipulation of Mie resonances in the particles^{1,3,16,22,29}. At the same time, low order Mie resonances in the lowrefractiveindex particles are usually spectrally very wide and thus cannot lead to sharp transitions in the scattering properties. Also low order Mie resonances of adjacent particles will strongly interact in the PhG, which is difficult to take into account. To the best of our knowledge, for the first time, we are providing a comprehensive theoretical and simulation treatment of structural colors employing photonic glasses based on firstorder Born approximation. It helps to explain the main mechanisms of color generation and supplies clear design and synthesis rules to achieve high color saturation.
In this work, we describe the relationships between the spectrum of a noniridescent structural color and the FT of the permittivity distribution via the firstorder approximation. We split the PhG structure into the disordered lattice and the repeating motif. We show that sharp transitions in the FT of the PhG structure can be obtained by tailoring the substructure of the motif which leads to a shift of the first zero position of the motif Fourier transform to smaller wave numbers. Numerical simulations confirm the appearance of sharp transitions in the reflection spectra for the optimized structures.
Results and Discussion
The light scattering properties of disordered structures with small permittivity perturbation \({\rm{\Delta }}\varepsilon (\vec{r})\) with respect to the background level can be estimated from a firstorder Born approximation^{30} which is given in the supplementary materials. The total electric field in the homogeneous material with a perturbation \({\rm{\Delta }}\varepsilon (\vec{r})\) can be expanded in a Taylor series with respect to the amplitude of \({\rm{\Delta }}\hat{\varepsilon }\)^{31}. For a small value of \({\rm{\Delta }}\hat{\varepsilon }\), the total field can be reasonably well approximated by the firstorder term, which represents the firstorder approximation. The firstorder term of the scattered wave contains as a source the excess polarization induced by the incident wave due to the permittivity perturbation. For the validity of the approximation the scattered wave amplitude should be much smaller than the amplitude of the incident wave, which is the nondepleted input criterion. Particularly helpful tool which derives from the firstorder approximation is the Ewald sphere construction which geometrically predicts the wavelength dependence and the directions of the scattered light^{30,32}. Figure 1 shows the schematic representation of the Ewald sphere construction for reflection from PhG. The PhG (Fig. 1a) has the FT of \({\rm{\Delta }}\varepsilon (\vec{r})\) in the shape of a spherical shell (Fig. 1b). The thickness of the shell in reciprocal space is related to the positional order of the spheres and will be discussed later. The wave vector of the incident light is ending at the origin of the reciprocal space. The length of the wave vector \({\vec{k}}_{in}\) is defined by the frequency ω and speed c of light and refractive index of the background material n_{ b }: \({k}_{in}={n}_{b}\omega /c\). The Ewald sphere has the radius of the incident wave number and is centered at the starting point of the incident wave vector. The scattering directions are defined by the scattering wave vectors \({\vec{k}}_{s}\) starting at the center of the Ewald sphere and pointing to the overlap regions between Ewald sphere and FT of \({\rm{\Delta }}\varepsilon (\vec{r})\). When the incident light has a small wave number (long wavelength, right image in Fig. 1a,b), there is no overlap between the Ewald sphere and FT of \({\rm{\Delta }}\varepsilon (\vec{r})\), so that the light cannot be scattered. When the wave number is increasing (intermediate wavelength, the middle image in Fig. 1a,b), the Ewald sphere starts to overlap with the FT of \({\rm{\Delta }}\varepsilon (\vec{r})\) and the incident light will be backscattered only. If we further increase the wave number (short wavelength, left image in Fig. 1a,b), the light will be backscattered into a cone of light with its opening angle spreading as the wavelength is further reduced. According to that, the expected reflection of the structure is schematically shown in Fig. 1c. There will be no light reflection for long wavelengths. Then, when the wavelength decreases such that the Ewald sphere overlaps with the FT of \({\rm{\Delta }}\varepsilon (\vec{r})\), the incident light starts to be reflected. The reflected power (P) from a scattering volume increases proportionally to the square of the absolute value of FT of \({\rm{\Delta }}\varepsilon (\vec{r})\) integrated over the Ewald sphere surface (ESS) (see supplementary materials):
where I_{0} is the intensity of the incident plane wave of light, θ is the angle between scattered \({\vec{k}}_{s}\) and input \({\vec{k}}_{in}\) wavevectors and for unpolarized light \(g(\theta )=(1+{\cos }^{2}\theta )/2\). Besides the dielectric strengths of the individual scatterers it is the overlap of the Ewald sphere with the square of the FT in reciprocal space which governs the reflected power. This way the light scattering can be fully analyzed from the FT of the permittivity. Please note that the FT of \({\rm{\Delta }}\varepsilon (\vec{r})\) alone does not define the scattering directions, and \(\vec{k}\) is not the scattering vector but the difference of the wave vector of the scattered and incident waves. The lightreflection transition between the noreflection and backreflection (shown in Fig. 1c) is determined by the sharpness, i.e. the slope of the FT spectrum of \({\rm{\Delta }}\varepsilon (\vec{r})\) at the inner boundary of the spherical shell. In other words, if we want to achieve a sharp reflection edge, the contributions of the square of the FT of \({\rm{\Delta }}\varepsilon (\vec{r})\) inside the shell should be as little as possible and, most importantly, the transition to large values should be sharp.
In order to produce the sharp reflection edge, the square of the FT of \({\rm{\Delta }}\varepsilon (\vec{r})\) of the PhG structure should be understood and tailored. As can be seen from Fig. 2a, the PhG structure can be seen as the convolution of the disordered lattice function \(l(\vec{r})\) with the motif function \(m(\vec{r})\), where \(l(\vec{r})\) represents the distribution of the spheres’ center points in space and \(m(\vec{r})\) represents the distribution of the permittivity difference in the single spherical particle. Mathematically, the FT of the whole structure is the multiplication of the lattice FT \({ {\mathcal F} }_{l}(\vec{k})\) and the motif FT \({ {\mathcal F} }_{m}(\vec{k})\) as shown in Fig. 2b:
The lattice of the PhG can be characterized by the minimum distance a between lattice points and the average coordination number (which is directly related to the packing density of the PhG constructed from solid spheres^{33}). For the case of hard spheres as shown in Fig. 2a, a is equal to the particle diameter d. We are not interested in the particular realization of the lattice but in the average absolute square of the FTs of such lattice realizations. Such an average quantity normalized for one lattice point is also called the structure factor \({\mathscr{S}}(\vec{k})\)^{33,34}. \({\mathscr{S}}(\vec{k})\) is a dimensionless function which approaches one for complete disorder, meaning that there is no correlation of motif positions and that the scattered intensity per particle in the matrix is equal to the intensity scattered by a single particle. In the case of incomplete disorder the correlation between particles leads to the fact that scattering into some direction is enhanced and in other directions is reduced. Thus in some directions the structure can scatter more than just the sum of intensities from single particles and \({\mathscr{S}}(\vec{k})\) is larger than one. The approximate function of radial distribution of \({\mathscr{S}}(\vec{k})\) can be derived from solving the OrnsteinZernike integral equation by choosing the hard sphere PercusYevick approximation^{35,36,37}. Due to the spherical symmetry of the structure factor \({\mathscr{S}}(\vec{k})\) in \(\vec{k}\) space resulting from the isotropy of the spatial lattice, \({\mathscr{S}}(\vec{k})\) is a function only of the absolute distance \(\vec{k}=k\) from the origin of the \(\vec{k}\) space and we can write \({\mathscr{S}}(\vec{k})={\mathscr{S}}(k)\) and assume \({\mathscr{S}}(k)\) to be real:
where \(\bar{N}\) is the average number density of the spheres, i.e. the number of spheres per unit volume, d is the diameter of the sphere with a = d (hard spheres), C(k) (the full equation of which can be found in supplementary material) is the Fourier transform of the direct correlation function c(r) which represents the direct interactions between particles, and \(\varphi =(\pi \bar{N}{d}^{3})/6\) is the sphere packing density. The coefficients are defined as \({\xi }_{1}={(1+2\varphi )}^{2}/{(1\varphi )}^{4}\), \({\xi }_{2}=\,{(1+\varphi /2)}^{2}/{(1\varphi )}^{4}\). As can be seen from the equations, the structure factor depends on the sphere packing density and the minimum distance between lattice points which is equal to diameter of the sphere for the solid hard spheres. Thus, the average intensity of scattered light per motif in the lattice is proportional to the product of the structure factor \({\mathscr{S}}\) and the form factor \({\mathscr{P}}\) where the latter is the square of the motif FT normalized by the volume of the motif (V) and can be written as \({\mathscr{P}}={{ {\mathcal F} }_{m}}^{2}/{V}^{2}={ {\mathcal F} }_{m}^{2}/{V}^{2}\) because the FT of our spherically symmetric particles is always real. Thus the average square of the FT from N particles is:
In the following examples we will consider PhGs with ultimate packing density of 64%^{38}. The manufacturing packing densities are slightly lower but that does not change the presented approach and conclusions.
The motif can be a homogeneous solid particle, an air hole, a coreshell particle or a hollow particle, etc. The form factor amplitude of a solid sphere with permittivity contrast Δε has the following radial function in reciprocal space^{33,39}:
where d is the particle diameter. Next, we consider a particular example of a direct PhG where the motif is the silica sphere in the background air.
Figure 2b shows the result for the PhG out of homogenous solid silica spheres (\(n=1.46\)) with the diameter d = a and with packing density of 64% embedded in air. The product \({\mathscr{S}}\cdot {\mathscr{P}}\) is a spherically symmetric function and thus only a 1D intensity spectrum along the radial direction is shown in Fig. 2b. For our considerations, the packing density is 64%, so the first peak of the lattice FT is located at about \({k}_{lp}=1.21(2\pi /a)\). For a smaller packing density the main peak of the structure factor will broaden and will slightly shift to smaller wave numbers^{34}. As can be seen, the \({\mathscr{S}}\cdot {\mathscr{P}}\) in the smaller k region has relatively large intensities and first peak (located approximately at k_{ lp }) has a smooth left edge. However, the right side of the peak has a sharper edge. This is because the firstzero point of \({ {\mathcal F} }_{m}\) (at k_{m0}) is located at the right side of the lattice peak at \({k}_{lp}=1.21(2\pi /a)\) (k_{m0} > k_{ lp }) leading to a sharper right peak edge of the product function \({\mathscr{S}}\cdot {\mathscr{P}}\).
Having understood this mechanism, the main idea behind our work is to influence the zero position of the form factor function such that it causes a maximum slope of the product function \({\mathscr{S}}\cdot {\mathscr{P}}\) on the lowk part of the spectrum, thus for the longwavelength edge, as well as a practically empty reciporcal space towards lowk numbers. These two properties eventually will lead to a reflection behavior of the photonic glass which yields a pure and highly saturated structural color. Therefore, if we move k_{m0} to the left side of the peak (k_{m0} < k_{ lp }), we could get a sharper left edge and lower intensity in the smallk region. We will now follow this approach with a coreshell sphere as a motif.
When the packing method and the lattice parameter a are fixed, the lattice function will not change. Accordingly, we can move the position of k_{m0} by modifying substructure of the motif. The coreshell particle structure (Fig. 3a,b) can be described as a sphere with permittivity contrast of the shell material (\({\Delta }{\varepsilon }_{1}={\varepsilon }_{s}{\varepsilon }_{b}\)) plus a smaller sphere with permittivity contrast between core and shell materials (\({\Delta }{\varepsilon }_{2}={\varepsilon }_{c}{\varepsilon }_{s}\)). So that the FT of the whole particle can be written by the following formula:
First, we consider the structure with a monotonous change of the refractive index from core through shell into the background. Thus Δε_{1} and Δε_{2} can be both positive or both negative. Figure 3c shows the FT of the coreshell sphere (Fig. 3a) with the corresponding refractive index of ε_{ c } = 2, ε_{ s } = 1.5 and ε_{ b } = 1 (background is air). The ratio of the core diameter to the diameter of the whole sphere is d_{ c }/d = 0.5. In this case, \({\Delta }{\varepsilon }_{1}=0.5\) and \({\Delta }{\varepsilon }_{2}=0.5\) have the same sign, k_{m0} of the coreshell sphere locates at the k position between the k_{m0} of the solid core and of the solid shell as can be seen from Fig. 3c. Thus, the k_{m0} is always on the right side of the lattice peak.
Now we consider an alternative situation where the refractive index is not changing monotonously. Thus, Δε_{1} and Δε_{2} have different signs. Figure 3d shows the FT of the coreshell sphere (Fig. 3b) with the corresponding refractive of ε_{ c } = 1 (core is air), ε_{ s } = 2 and \({\varepsilon }_{b}=1.5\). The coreshell ratio of the diameter is the same as in Fig. 3a but \({\Delta }{\varepsilon }_{1}=0.5\) and \({\Delta }{\varepsilon }_{2}=\,1\) have different signs. In this case the zero position k_{m0} of the normalized coreshell sphere FT amplitude is located at a k position smaller than k_{m0} of the core or shell as can be seen from Fig. 3d. Thus, by a proper choice of refractive index contrast and shell thickness the motif zero point can be moved to the left side of the lattice peak. We have derived an equation to determine shell thickness at given refractive index contrast in the supplementary materials.
The main finding of the presented investigation of the influence of the varying core diameter in relation to the sphere diameter in coreshell particles with nonmonotonous refractive index distribution is the fact that the zero position of the motif FT can be positioned anywhere from the right to the left side of the peak of the lattice structure factor function. Therefore, we obtain a design degree of freedom which allows us to manipulate the overall FT of the PhG, thus the quality of its structural color properties, by mere changing of the particle geometry.
We now consider a particular coreshell particle system to implement a sharp transition in the reflection spectrum. We first optimize the FT of the structure and later compare finite integration technique (FIT) simulations for PhGs with and without optimization. Figure 4 shows\(\,{\mathscr{S}}\), \({ {\mathcal F} }_{m}/V\) and \({\mathscr{S}}\cdot {\mathscr{P}}\) for the hollow sphere PhG of zirconia particles with aircores and with a background of air. We consider this particular example as such particles can be readily synthesized and such PhGstructures can be obtained by coassembly^{40}. More examples of combinations of different materials for core and shell are presented in the supplementary materials. The refractive indices are n_{ c } = 1, n_{ s } = 2.12, and n_{ b } = 1 for the coreshell sphere. The particle diameter and minimum lattice distance are \(d=a=221.8\) nm to obtain transition for blue color, the considered packing density of the spheres is assumed to be the theoretical limit for PhG of 64%. As can be seen from Fig. 4, the k_{m0}can be moved from the right side to the left side of k_{ lp } when we increase dc/d from 0.4 to 0.9. For d_{ c } = 86.2 nm shown in Fig. 4a, the zero point of \({ {\mathcal F} }_{m}\) is located on the right side of the lattice peak (\({k}_{m0} > {k}_{lp}\)). In this case, the left side edge of the product \({\mathscr{S}}{\mathscr{P}}\) is very smooth and there is a substantial contribution at smaller k values. When we increase the core diameter to 135.5 nm, the zero position shifts to a smaller k value, and just overlaps with the peak (\({k}_{m0}\approx {k}_{lp}\)). As can be seen from Fig. 4b, the first peak of \({\mathscr{S}}\cdot {\mathscr{P}}\) almost vanishes. When we further increase the diameter of the core to 198 nm, k_{m0} is moved to the left side of the lattice peak (\({k}_{m0} < {k}_{lp}\)). As can be seen from Fig. 4c, the left side edge of the first peak of \({\mathscr{S}}\cdot {\mathscr{P}}\,\)is very steep, and now the smallk region has very low intensity. This kind of motif should lead to sharp reflection transition with low scattering for larger wavelengths.
To predict the reflection curves we calculated the effective scattering cross section per single particle in the PhG:
It can be derived as the effective power scattered per single particle divided by the incident light intensity using equations (1) and (5). Figure 5 shows the scattering cross section per particle normalized by the particle geometrical cross section which represents the particle scattering efficiency compared for two PhGs: the homogenous solid silica particle PhG (with the FT shown by Fig. 2b) is shown by Fig. 5a and the coreshell particle PhG (with the FT shown by Fig. 4c) is shown by Fig. 5b, respectively. As can be seen from Fig. 5, the hollow sphere PhG shows a sharper scattering transition from no scattering to strong scattering than the full sphere PhG. Thus, the full sphere PhG shows significant scattering in 500–600 nm range, whereas the hollow sphere PhG demonstrates almost no scattering in this range, which is the basis for superior color saturation. It should be noted that the scattering wavevector \({k}_{s}={n}_{m}\omega /c\) is used as the radius of the Ewald sphere, where n_{ m } is the mean refractive index of the PhG structure. Different sphere diameters are chosen to compensate for the different mean refractive index and by that match the transition wavelengths. It can be seen in Fig. 5 that the effective scattering per particle is much smaller than one in the long wavelength range. Thus, interference of scattering from many particles and thus short range order plays the significant role in PhG in this range. At short wavelengths the firstorder approximation will fail and other approximate methods should be considered. This is beyond the scope of this work.
The applicability of the firstorder approximation is limited to a small refractive index contrast and thin PhG films, e.g. for film thicknesses which do not substantially change the power of the light propagating through the film. For the zirconia based structures considered here, the presented approach thus can be considered only as an approximate solution. Therefore, a very important question we needed to answer was that on the predictive power of our firstorder approach in view of real structures with substantial refractive index contrasts. To address this question we performed bruteforce 3D FIT simulations on randomly packed PhG assemblies of coreshell spheres.
We have numerically simulated the PhGs presented in Fig. 5. The structure realizations of both randomly packed PhGs of solid silica and of hollow zirconia spheres with 64% packing density were obtained by the packing generator MUSEN^{41,42}. The periodic boundary condition was used in the packing algorithm to avoid packing density variations at the edges of the packed volume. Under normal incidence, the light reflectance (R) of the PhG was simulated by using the finite integration technique simulation with CST Microwave Studio^{43}. For the homogeneous solid silica particle PhG, the size of the simulated structure is 2.5 × 2.5 × 12.3 µm^{3} (number of particles is 9122). The corresponding analytical permittivity FT of the structure is shown in Fig. 2b. Here, we assume a sphere diameter of 215.6 nm. For the hollow zirconia sphere PhG, the size of the simulated structure is 1.6 × 1.6 × 7.9 µm^{3} (number of particles is 2961). The simulated volume size here is adjusted to make the maximal reflection comparable to that of solid sphere PhGs. The corresponding analytical permittivity FT of the structure composed of particles with 221.8 nm outer diameter and air cores of 198 nm diameter is shown in Fig. 4c.
The PhGs are then excited by plane wave incident vertically from air. The lateral sides of the simulation volume are mirrors such that light can exit the simulation volume only through the open boundaries at the top and the bottom. The bottom of the PhG is adjusted to the homogeneous substrate material with refractive index equal to the average refractive index of PhG in order to minimize reflections from this boundary. The homogeneous substrate material is then terminated by an open boundary condition. The reflected power is calculated as Poynting vector integration over the upper boundary. We present here simulations of arbitrarily chosen single realizations of the hollow sphere and solid sphere PhGs and do not average over many simulations with different realizations. As light interacts with thousands of particles in a single simulation, we believe the volume averaging is sufficient to represent the reflection spectrum and do not conduct ensemble averaging. The ensemble averaging would also go beyond the currently available simulation capacity.
The light reflection spectra for (a) the solid silica sphere PhG and (b) the hollow zirconia sphere PhG are shown in Fig. 6. These light reflection spectra show the same trend as described in Fig. 1c. In the longer wavelength region, there is nearly no light reflection. When the wavelength of incident light decreases, the Ewald sphere starts to overlap with permittivity FT of the structure, so the incident light starts to be reflected. The remarkable detail is that the light reflection spectrum of the hollow sphere PhG shows a much sharper transition compared to that of the homogenous sphere PhG. There is a good correspondence at small reflection values between the scattering efficiency per particle and the simulated reflectivity functions of the PhGfilms (compare Figs 5 and 6). The enhancing effect of the presented sharp transition in reflectivity on blue color appearance is discussed in the supplementary material. For reflectivity larger than 50% the firstorder approximation is not applicable as the incident wave becomes strongly depleted and the assumption of the same incident intensity on each particle is not justified anymore.
However, we are interested in analyzing the slope of the transition from very low reflection at long wavelengths to high reflection at shorter wavelengths as this slope is responsible for the saturation of the structural color. As it turned out this slope is very well predicted by the firstorder approximation even considering large refractive indices, such as that of zirconia. We can therefore state that the firstorder Ewald sphere approach is a simple and very wellfunctioning technique for predicting the quality of structural colors.
Conclusion
In conclusion, the selective reflectivity of the PhG structure is related to the spherical shellshaped Fourier transform of its permittivity distribution. To explain the connection between the two properties we employ the Ewald sphere construction resulting from firstorder approximation. For sharp spectral selectivity the radial distribution of permittivity Fourier transform should obtain a peak with a sharp edge at lower wave numbers. We have shown that coreshell particles can be used with nonmonotonous refractiveindex distribution to achieve this property. Namely, the permittivity difference Δε_{1} between shell and background should have opposite sign compared to the permittivity difference Δε_{2}. between shell and core. That is, the Fourier transform of the photonic glass can be modified by changing the substructure of the PhG motif. In the optimal situation the zero point of the motif Fourier transform is positioned just at the small wave number ed of the peak corresponding to the Fourier transform of the PhG lattice. Numerical simulations using the finiteintegration timedomain simulation confirm that the structure with optimized motif has sharp reflection transition. A particular example of hollow sphere was presented. But much more combinations are possible when two materials are combined in the coreshell sphere or background porosity of PhG is filled with a third material. The presented examples are based on a PhG packing density of 64%. For a smaller packing density the peak position of the structure factor will shift to smaller wavenumbers. Thus the zero point of the motif FT should be also shifted to smaller wavenumbers, which can still be adjusted by the motif optimization. The proposed innovative approach paves the road for novel structural colors with high color saturation.
References
 1.
Kinoshita, S., Yoshioka, S. & Miyazaki, J. Physics of structural colors. Rep. Prog. Phys. 71, 076401 (2008).
 2.
Faulkner, E. B. & Schwartz, R. J. High Performance Pigments, (John Wiley & Sons, 2009).
 3.
Dumanli, A. G. & Savin, T. Recent advances in the biomimicry of structural colours. Chem. Soc. Rev. 45, 6698–6724 (2016).
 4.
Eppler, R. A. & Eppler, D. R. Which colors can and cannot be produced in ceramic glazes. Ceram. Eng. Sci. Proc. 15, 281 (1994).
 5.
Wendusu, Y. T., Masui, T. & Imanaka, N. Novel environmentally friendly inorganic red pigments based on calcium bismuth oxide. J. Adv. Ceram. 4, 39 (2015).
 6.
Jansen, M. & Letschert, H. P. Inorganic yellowred pigments without toxic metals. Nature 404, 980–982 (2000).
 7.
Joannopoulos, J. D., Johnson, S. G., Winn, J. N. & R. D. Meade, Photonic crystals: Molding the flow of light, 2nd ed. (Princeton University Press: Princeton, 2008).
 8.
Prum, R. O., Torres, R. H., Williamson, S. & Dyck, J. Coherent light scattering by blue feather barbs. Nature 396, 28–29 (1998).
 9.
Prum, R. O., Torres, R., Williamson, S. & Dyck, J. Twodimensional Fourier analysis of the spongy medullary keratin of structurally coloured feather barbs. P. Roy. Soc. B: Biol. Sci. 266, 1322 (1999).
 10.
Noh, H. et al. How noniridescent colors are generated by quasiordered structures of bird feathers. Adv. Mater. 22, 2871–2880 (2010).
 11.
Saranathan, V. et al. Structure and optical function of amorphous photonic nanostructures from avian feather barbs: a comparative small angle Xray scattering (SAXS) analysis of 230 bird species. J. R. Soc. Interface 9, 2563–2580 (2012).
 12.
García, P. D., Sapienza, R., Blanco, Á. & López, C. Photonic glass: A novel random material for light. Adv. Mater. 19, 2597–2602 (2007).
 13.
García, P. D., Sapienza, R. & López, C. Photonic glasses: A step beyond white paint. Adv. Mater. 22, 12–19 (2010).
 14.
Park, J. G. et al. Fullspectrum photonic pigments with noniridescent structural colors through colloidal assembly. Angew. Chem. Int. Ed. 53, 2899–2903 (2014).
 15.
Forster, J. D. et al. Biomimetic isotropic nanostructures for structural coloration. Adv. Mater. 22, 2939–2944 (2010).
 16.
Ge, D. et al. A robust smart window: Reversibly switching from high transparency to angleindependent structural color display. Adv. Mater. 27, 2489–2495 (2015).
 17.
Dong, B. Q. et al. Structural coloration and photonic pseudogap in natural random closepacking photonic structures. Opt. Express 18, 14430–14438 (2010).
 18.
Magkiriadou, S., Park, J. G., Kim, Y. S. & Manoharan, V. N. Disordered packings of core shell particles with angleindependent structural colors. Opt. Mater. Express 2, 1343–1352 (2012).
 19.
Magkiriadou, S., Park, J. G., Kim, Y. S. & Manoharan, V. N. Absence of red structural color in photonic glasses, bird feathers, and certain beetles. Phys. Rev. E 90, 062302 (2014).
 20.
Zhang, Y. et al. Using cuttlefish ink as an additive to produce noniridescent structural colors of high color visibility. Adv. Mater. 27, 47194724 (2015).
 21.
Wang, F. et al. Rapid fabrication of angleindependent structurally colored films with a superhydrophobic property. Dyes Pigments 130, 202–208 (2016).
 22.
Takeoka, Y. et al. Production of colored pigments with amorphous arrays of black and white colloidal particles. Angew. Chem. Int. Ed. 52, 7261–7265 (2013).
 23.
Shawkey, M. D., Morehouse, N. I. & Vukusic, P. A protean palette: Colour materials and mixing in birds and butterflies. J. R. Soc. Interface 6, S221–S231 (2009).
 24.
Shawkey, M. D. & Hill, G. E. Carotenoids need structural colours to shine. Biol. Letters 1, 121–124 (2005).
 25.
Kawamura, A. et al. Fullcolor biomimetic photonic materials with iridescent and noniridescent structural colors. Sci. Rep. 6, 33984 (2016).
 26.
Iwata, M., Teshima, M., Seki, T., Yoshioka, S. & Takeoka, Y. Bioinspired bright structurally colored colloidal amorphous array enhanced by controlling thickness and black background. Adv. Mater. 29, 1605050 (2017).
 27.
Kim, S. H. et al. Inverse photonic glasses by packing bi disperse hollow microspheres with uniform cores. ACS Appl. Mater. Interfaces 9, 24155–24160 (2017).
 28.
Kohri, M., Yamazaki, S., Kawamura, A., Taniguchi, T. & Kishikawa, K. Bright structural color films independent of background prepared by the dipcoating of biomimetic melaninlike particles having polydopamine shell layers. Colloids Surf. A: Physicochem. Eng. Aspects 532, 564–569 (2017).
 29.
Montesdeoca, D. et al. Monodisperse silica spheres ensembles with tailored optical resonances in the visible. Part. Part. Syst. Charact. 33, 871–877 (2016).
 30.
Born, M. & Wolf, E. Principles of optics: Electromagnetic theory of propagation, interference and diffraction of light, 7th ed. (Cambridge Univ. Press, 1999).
 31.
Lang, S. et al. A quantitative firstorder approach for the scattering of light by structured thin films, arXiv:1802.07058.
 32.
Maiwald, L. et al. Ewald sphere construction for structural colors. arXiv:1803.01755.
 33.
Pedersen, J. S. Modelling of smallangle scattering data from colloids and polymer systems in neutrons, Xrays and light, (Elsevier Science B.V, 2002).
 34.
lópez de Haro, M. & Roble, M. The structure factor and equation of state of hardsphere fluids. J. Phys.: Condens. Matter. 16, S2089–S2096 (2004).
 35.
Wertheim, M. S. Exact solution of the PercusYevick integral equation for hard spheres. Phys. Rev. Lett. 10, 321–323 (1963).
 36.
Ashcroft, N. W. & Lekner, J. Structure and resistivity of liquid metals. Phys. Rev. 145, 83–90 (1966).
 37.
Doicu, A., Wriedt, T. & Eremin, Y. A. Light scattering by systems of particles: Nullfield method with discrete sources: Theory and programs, (Springer, 2006).
 38.
Torquato, S., Truskett, T. M. & Debenedetti, P. G. Is random close packing of spheres well defined? Phys. Rev. Lett. 84, 2064–2067 (2000).
 39.
Vembu, S. Fourier transformation of the ndimensional radial delta function. Quart. J. Math. 12, 165–168 (1961).
 40.
do Rosario, J. J. et al. Facile deposition of YSZinverse photonic glass films. ACS Appl. Mater. Interfaces 6, 12335–12345 (2014).
 41.
Dosta, M., Antonyuk, S. & Heinrich, S. Multiscale simulation of agglomerate breakage in fluidized beds. Ind. Eng. Chem. Res. 52, 11275–11281 (2013).
 42.
Dosta, M. et al. Numerical and experimental analysis of influence of granule microstructure on its compression breakage. Powder Technol. 299, 87–97 (2016).
 43.
CST, Darmstadt Germany, Microwave Studio software. www.cst.com.
Acknowledgements
The authors gratefully acknowledge financial support from the German Research Foundation (DFG) via SFB 986 “TailorMade MultiScale Materials Systems: M³”, projects C2, A3 and priority program SPP 1839 “Tailored Disorder”, project 278744289, and the Hamburg University of Technology (TUHH) in the funding program Open Access Publishing, DFG project 392323616. The authors also acknowledge the support from CST, Darmstadt, Germany, with their Microwave Studio software.
Author information
Affiliations
Contributions
G.S. performed the simulations, A.P. and M.E. supervised the project, G.S., L.M., H.R., D.J., M.D., S.H., A.P. and M.E. analysed the results and wrote the paper.
Corresponding author
Ethics declarations
Competing Interests
The authors declare no competing interests.
Additional information
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Shang, G., Maiwald, L., Renner, H. et al. Photonic glass for high contrast structural color. Sci Rep 8, 7804 (2018). https://doi.org/10.1038/s41598018261198
Received:
Accepted:
Published:
Further reading

Transparency induced in opals via nanometer thick conformal coating
Scientific Reports (2019)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.