Issue 
EPJ Appl. Metamat.
Volume 7, 2020
Metamaterial Research Updates from China



Article Number  6  
Number of page(s)  6  
DOI  https://doi.org/10.1051/epjam/2020006  
Published online  23 November 2020 
https://doi.org/10.1051/epjam/2020006
Research Article
Optical bistability in twodimensional nonlinear composites of coated cylinders with nonlinear core and graded shell
School of Science, Jiangnan University, Wuxi 214122, PR China
^{*} email: wxwym2098@163.com
^{**} email: yanghuang@jiangnan.edu.cn
Received:
10
January
2020
Accepted:
21
September
2020
Published online: 23 November 2020
The intrinsic optical bistability of the nonlinear composite consisting of coated cylinders with nonlinear core and graded Drude shell is investigated with selfconsistent meanfield approximation. We derive the nonlinear equation of near field for such graded composite in the quasistatic limit. We demonstrate that the bistability threshold and the bistable width are highly dependent on the core radius, the incident frequency, and the graded coefficient of the coated cylinder inclusion. It is found that the optical bistability appears only when the core radius and the incident frequency satisfy some specific conditions. Therefore, the optical bistability for nonlinear composite materials may be optimized by the suitable adjustment of the physical and geometrical parameters. Our results may be helpful for the design of the nonlinear device with appropriate bistability.
Key words: Optical bistability / cylindrical / graded particle / composite medium
© Y.M. Wu et al., published by EDP Sciences, 2020
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
Optical bistability has attracted widespread interests for decades, due to its novel feature of controlling light with light [1–6]. For twocomponent nonlinear composite system, variational methods [7] and selfconsistent meanfield approximation [8,9] were developed to investigate the nonlinear response of composite system and to explore optical bistable behavior [7–9]. It is wellknown that the local field enhancement of composite system may result in the optical bistability, which can be adjusted by controlling the local microstructure and the physical parameters of components in the composite system. With the development of nanotechnology, preparation of some special microstructures such as coreshell nanoparticle becomes possible [10–12]. For metal/semiconductor coreshell structure composite system, it was found that there exists large enhancement of the optical nonlinear susceptibility [13]. In addition, studies of nonlinear susceptibilities generated by two (three) harmonics in isotropic (or anisotropic) composite systems with coreshell nanoparticles were also widely investigated [14,15]. For instance, researchers found that the twodimensional nonlinear composites exhibit huge enhancement of the thirdorder nonlinear optical susceptibility under the partial resonant condition [16,17]. Moreover, for the threedimensional nonlinear composite consisting of the coreshell nanostructures, the nonlinearity of the whole composite system can be dramatically enhanced by adjusting the material properties of coreshell particle composition ratio [18,19]; Furthermore, Gao et al. [20,21] developed the selfconsistent meanfield approximation to study the nonlinear composite media consisting of coated nanoparticles with nonlinear core and nonlinear shell. Recently, Shewamare et al. [22] studied the optical bistability of the random composites with nonlinear metallic nanoparticles. And Jackson et al. [23] showed that local field can get extremely high magnitude of Raman scattering enhancement in the order of 10^{6} in the nanoparticles with coreshell geometry.
On the other hand, effective linear dielectric constant [24], threeorder nonlinear susceptibility [25] and two(three) harmonics generations for the nonlinear susceptibilities [26] of the random granular composites with the graded coreshell structure were paid a lot of attentions. In this connection, Huang et al. [27,28] found that the gradient structure allows linear and nonlinear responses of the composite to show a wide plasmonic resonance band, which helps to enhance the effective nonlinear response of the composite. Later, Wei et al. [29] studied the local field distribution in nanoparticles with gradient shell and showed that the electric field in the shell is adjustable. More recently, third order nonlinear optical properties of magnetite nanoparticles in colloid were investigated experimentally as a sample of extension from electrical field to magnetic field [30].
To the best of our knowledge, the above works for graded composites only focus on the weak incident light intensity, and only the effective linear and nonlinear responses of the graded system were considered. It is well known that when the light intensity is strong, the composite system may exhibit optical bistability. Therefore, in this paper, we would like to investigate the optical bistable behavior of the graded composite media. Actually, under the surface plasmon resonant conditions, large enhancement of the local electric field takes place and results in the appreciable nonlinear bistability. With further enhancement of light intensity the combined the coreshell microstructure configuration and the gradient characteristic of cylinder components may provide us original optical bistable behavior. We believe our study may be helpful for the design and preparation of optical bistable materials and relative devices.
2 Theoretical model
We assume that the composite is composed of dilute coreshell cylinder inclusions randomly embedded in the linear host medium as shown in Figure 1 in 2D. The distance between inclusions is long enough to neglect the electrostatic dipole interactions among the inclusions in the dilute condition. The coreshell inclusions in this threecomponent composite consist of Kerrtype cores and graded Drudetype shells. Here, we make the assumption that the dielectric function of the graded materials varies along the radial direction r in cylindrical coordinates (r, ϕ, z), and can be written as ϵ _{ α }(r) with the subscript α = c, s, h indicating the core, the shell, and the host media respectively. The inner and outer radii of the cylindrical shell are assumed to be a and b. Within the quasistatic approximation, the electric potentials in whole space satisfy the following equation,(1)where φ _{ α } is electric potential in each area and corresponding local electric field is E _{ α } = −∇ φ _{ α }. For infinitely long cylindrical configuration, Equation (1) reduces to(2)
Therefore, the electric potentials in each region have the following general expressions,(3) (4) (5)where A, B, A _{1}, B _{1} are four unknown coefficients to be determined, P _{ m }(cosϕ) is the mth order Legendre polynomials, and the radial function R _{ m }(r) in shell region satisfies the following equation(6)
and are the two solutions which are regular at the origin and infinity, respectively. Applying the continuous boundary conditions of the normal components of the electric potentials and electric displacements, we can obtain the coefficients [29] as,(7) (8) (9) (10) (11) (12) (13) (14) (15)
where B in equation (8) is equivalent to the dipole factor for the system and F(b) in equation (14) is the equivalent dielectric constant of particles which are determined by the function of the gradation.
To further investigate the distribution of electric potentials in each region, we assume the cylindrical shell is described by the graded Drude model [24,25](16)where ω denotes the frequency of the external electric field, ω _{p}(r) is the graded plasma frequency which varies spatially in radial direction and Γ is the relaxation rate. Without loss of generality, we normalize the powerlawform graded plasma frequency as where h and k are two gradient constants denoting the gradation of the shell [31]. We further normalize the external field frequency ω and the relaxation rate Γ with ω_{p} , thus the Drude model can be rewritten as [29],(17)
Substituting equation (17) into equation (6), we yield(18) (19) where F (α, β, γ, hr ^{ k }/d) is the hypergeometric function which is analytic in the whole complex plane except at the singular points. In addition, , , γ _{±1} = ± 2/k + 1 and d = ω ^{2} + iωΓ − 1.
Therefore, we have analytically derived the potential solutions of graded cylindrical coreshell system with graded Drude model, hence it is easy to achieve the corresponding electric fields throughout the whole space with E _{ α } = − ∇ φ _{ α },(20) (21) (22)
In the dilute limit, we only consider the thirdorder nonlinear response of the core medium, in which the dielectric function of the core has a fielddependent expression as(23)
where χ _{c} is the thirdorder nonlinear susceptibility. Substituting for ϵ _{ c } in equation (7), we yield the nonlinear local field of core,(24)
In the case of strong external field, the second term in equation (23) becomes comparable to the first term, so that the optical bistability in certain conditions may occur. Equation (24) is the key result of this work based on what we will give some numerical calculations and discussions in the following part.
Fig. 1 Schematic diagram of the proposed composite in which the coreshell inclusions have Kerrtype dielectric cores and graded Drudetype shells. 
3 Numerical results and discussions
We are now in a position to study the optical bistability in the graded composite materials. For numerical calculations, we choose the following physical parameters for the nonlinear core, graded shell, and host media: Γ/ω _{p} = 0.02, ϵ _{c} = − 10 + 0.037i, ϵ _{h} = 64, χ _{c} = 10^{−11}(m/V)^{2} [27–29]. Note that the graded materials are rare in nature and most of them in study are artificial with some proper methods. As to the graded shell we prosed here, one possible way to achieve such gradation is to fabricate a graded metallic shell by using different noble metals as different layers covering the core. The large host dielectric constant is essential to achieve a bistable optical response. Such a large value of host medium may be realized by use of the CuCl material near the biexcitonic transition [32]. In reality, the dielectric constants of candidate materials who possess previous physical parameters vary with different frequency. Here we keep them constant for the sake of simplicity. Meanwhile, in the quasistatic approximation, the ratio of the inner and outer radii of the coated particle makes more sense than their exact value in the present model. Therefore, we keep the outer radius b constant and provide it with a dimensionless value 1, and vary the inner radius a from 0 to 1 in the following calculation.
Figure 2a shows the local field in core as a function of the external electric field. It is evident that the Sshaped loop is found which indicates the optical bistable behavior in the near field with strong incident intensity. The abscissa E _{0} of right and left inflection points in the bistable curve are, respectively, the upper and lower thresholds. Figure 2a illustrates that the upper and lower thresholds will simultaneously vary with the core radius a. It is found that both the upper and lower thresholds are increased with increasing the core radius in the small value regime, and then are decreased. Moreover, for some specific values of a, for instance a = 0.31, the bistable threshold has the maximal width with 4.584 × 10^{6} V/m as shown in Figure 2b; When a = 0.93 however, the optical bistable curves will disappear. This indicates a critical value for core radius a _{c} in the determination of optical bistablity within the present physical parameters, i.e., when the particle core radius a < a _{ c }, the bistability appears; The bistability would disappear if a > a _{ c }. Besides that, dash line in Figure 2b shows the dependence of upper and lower thresholds on the core radius for the nongraded shell (h = 0). Variation of bistable thresholds with a in the nongraded case is quite different from that in the graded case. In details, there exists a certain value of a after which the upper and lower thresholds will decrease in the graded case. In nongraded case however, the upper and lower thresholds vary slowly with a, but finally increase rapidly when a is close to 1. Note that all the optical bistability occur under strong field intensity condition and one could not achieve this bistable curve with weak field intensity. These results provide us a good theoretical guidance for the design of lowthreshold optical biastability.
Figure 3a shows the local electric field as a function of the external field with different frequency ω. The maximal width of the threshold is 4.584 × 10^{6} V/m which occurs at ω/ω _{p} = 0.7. Once ω/ω _{p} is down to 0.08, the upper and lower threshold curves begin to degenerate and finally the optical bistability disappear when ω/ω _{p} < 0.08. This indicates that there is a critical incident frequency ω _{c,} above which bistability appears and below which the bistability disappers. This may result from the fact that the average local field of core is weaker when ω < ω _{c}, so the bistablity does not appear. With the increase of ω, the local field also increases, hence leading to bistable phenomenon. Compared to the gradient case (h = 0.9) composite with nongradient media (h = 0) can achieve lower thresholds and larger width of optical bistability thresholds.
To illustrate the dependence of upper and lower thresholds on a and ω/ω _{p} in a boarder parameter space, we define the width of bistable threshold ΔE _{0} and plot it in Figure 4. Again, the width of bistable threshold is found tunable by adjusting the core radius a and external field frequency ω/ω _{p}. When a = 0.01 and ω/ω _{p} = 0.97, the maximal width of threshold is found to be 1.121 × 10^{7}(V/m). Thus, it is possible to achieve the optimal width of threshold by our suitable adjustment of a and ω/ω _{p} simultaneously.
Figure 5 illustrates bistable curves in near field with different host media, it is obviously found that the upper threshold reduces when the value of ϵ _{h} increases. To further demonstrate, we investigate the dependence of thresholds on the host media ϵ _{h} in Figure 6. It shows that the upper and lower thresholds are simultaneously reducing with increasing ϵ _{h}, especially in the low ϵ _{h} regime, i.e., from 1 to 4. After that, they vary slowly with ϵ _{h} and are almost unchanged when ϵ _{h} reaches to high values. The situation is quite different in the nongraded case (dashed line) that the upper threshold is increasing rapidly in the low ϵ _{h} regime from 1 to 6, while the lower thresholds firstly reduce and then slightly increase. In the high ϵ _{h} regime, the dependence of upper and lower thresholds on ϵ _{h} is similar to the graded case.
In the end, we study the dependence of the upper and lower thresholds on the gradient factor h. It can be seen in Figure 7 that there are two peaks for the upper thresholds. The first peak occurs at h = 0.52 that corresponds to the width of threshold of 5.637 × 10^{6} V/m, and the second happens at h = 0.9 with the width of threshold being 4.584 × 10^{6} V/m. When h continues to increase, both the upper and lower thresholds exhibit a downward trend.
Our results suggest that there exists the optical bistability under certain conditions for cylindrical composite consisting of nonlinear coated nanocylinders. This nonlinear feature mainly results from the surface plasmon resonance (SPR) inside the nanostructure, which can dramatically enhance the local electric field inside the nonlinear core, hence boosting the nonlinear response. Therefore the optical bistability occurs. Actually, the SPR enhanced optical nonlinearity resulting in the optical bistability in composite media are the conclusions demonstrated in many previous studies [33–35]. In addition, the introduction of gradient shell give a new degree of freedom for us to obtain desired optical bistability threshold. In the end, we shall remark that although all the parameters in the present manuscript might be realized experimentally, it still has many challenges in precisely preparing the proposed model and testing the bistable behavior in experiment.
Fig. 2 (a) The local electric field as a function of the external field with different core radii. (b) The upper and lower thresholds as a function of core radius a with h = 0.9 (solid line) and h = 0 (dashed line), respectively. 
Fig. 3 (a) The local electric field as a function of the external field with different frequency ω/ω _{p}. (b) The upper and lower thresholds as a function of frequency ω/ω _{p} with h = 0.9 (solid line) and h = 0 (dashed line), respectively. 
Fig. 4 The width of bistable threshold as a function of a and ω/ω _{p}. 
Fig. 5 The local field as a function of the external field with different host media. The gradation constant h = 0.9. 
Fig. 6 The upper and lower thresholds as a function of host medium with h = 0.9 (solid) and h = 0 (dashed). 
Fig. 7 The upper and lower thresholds as a function of h. 
4 Conclusion
Within quasistatic approximation, we derive the nonlinear equations of near field for the graded metaldielectric composite, which is composed of Kerr dielectric coregraded metallic shell nanocylinders randomly embedded in a linear host medium, and study the optical bistability in such composite. We demonstrate the dependences of the bistability threshold and width on the core radius, the frequency, and gradient coefficient. There exist a critical core radius a _{c} and a critical incident frequency ω _{c}, below (above) which the optical bistable behavior takes place. The results indicate that the gradient shell structure is in favor of appearance of optical bistability with lower threshold and larger width of threshold. Moreover, the character of optical bistability can be tunable by changing the geometrical and physical parameters of composite, to achieve maximal width of threshold and lower thresholds. The results can provide a theoretical basis for the design and preparation of optical bistable material.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (grant no. 11704158) and the Natural Science Foundation of Jiangsu Province (grant no. BK20170170).
References
 K.M. Leung, Phys. Rev. A 33, 2461 (1986) [CrossRef] [Google Scholar]
 R.W. Boyd, Nonlinear optics (New York: Academic, 1992) [Google Scholar]
 H.M. Gibbs, S.L. McCall, T.N.C. Venkatesan, Phys. Rev. Lett. 36, 113 (1976) [CrossRef] [Google Scholar]
 H. Chang, H.B. Wu, C.D. Xie, H. Wang, Phys. Rev. Lett. 93, 213901 (2004) [CrossRef] [PubMed] [Google Scholar]
 T.K. Mavrogordatos, et al., Phys. Rev. Lett. 118, 040402 (2017) [CrossRef] [Google Scholar]
 M. Kim, S. Kim, S. Kim, Opt. Express 26, 11620 (2018) [CrossRef] [Google Scholar]
 D.J. Bergman, O. Levy, D. Stroud, Phys. Rev. B 49, 129 (1994) [CrossRef] [Google Scholar]
 L. Gao, L.P. Gu, Z.Y. Li, Phys. Rev. E 68, 066601 (2003) [CrossRef] [Google Scholar]
 Y. Huang, L. Gao, Opt. Lett. 43, 2836 (2018) [CrossRef] [Google Scholar]
 N.A. Nicorovici, R.C. McPhedran, G.W. Miton, Phys. Rev. B 49, 8479 (1994) [CrossRef] [Google Scholar]
 K. Mallik, M. Mandal, N. Pradhan, T. Pal, Nano. Lett. 1, 319 (2001) [CrossRef] [Google Scholar]
 J.J. Schneider, Adv. Mater. 13, 259 (2001) [CrossRef] [Google Scholar]
 A.E. Neeves, M.H. Birnboim, J. Opt. Soc. Am. B 6 (1989) [Google Scholar]
 P.M. Hui, C. Xu, D. Stroud, Phys. Rev. B 69, 014203 (2004) [CrossRef] [Google Scholar]
 X.C. Jiang, Y.W. Zhou, D.L. Gao, Y. Huang, L. Gao, Opt. Express 28, 17384 (2020) [CrossRef] [Google Scholar]
 O. Levy, J. Appl. Phys. 77, 1696 (1995) [CrossRef] [Google Scholar]
 L. Gao, Phys. Status Solidi B 236, 182 (2003) [CrossRef] [Google Scholar]
 X.Y. Liu, Z.Y. Li, Phys. Lett. A 223, 475 (1996) [CrossRef] [Google Scholar]
 C. Xu, Z.Y. Li, Phys. Lett. A 252, 233 (1999) [CrossRef] [Google Scholar]
 L. Gao, Phys. Lett. A 318, 119 (2003) [CrossRef] [Google Scholar]
 L.P. Gu, L. Gao, Physica B 368, 275 (2005) [Google Scholar]
 S. Shewamare, V.N. Mal'nev, Physica B 407, 4837 (2012) [CrossRef] [Google Scholar]
 J.B. Jackson, S.L. Westcott, L.R. Hirsch, J.L. West, N.J. Halas, Appl. Phys. Lett. 82, 257 (2003) [CrossRef] [Google Scholar]
 J.P. Huang, K.W. Yu, G.Q. Gu, M. Karttunen, Phys. Rev. E 67, 054105 (2003) [CrossRef] [Google Scholar]
 L. Gao, J.P. Huang, K.W. Yu, Phys. Rev. B 69, 075105 (2004) [CrossRef] [Google Scholar]
 L. Gao, K.W. Yu, Phys. Rev. B 72, 075111 (2005) [CrossRef] [Google Scholar]
 J.P. Huang, K.W. Yu, J. Opt. Soc. Am. B 22, 1641 (2005) [Google Scholar]
 J.P. Huang, P.M. Hui, K.W. Yu, Phys. Lett. A 342, 484 (2005) [CrossRef] [Google Scholar]
 E.B. Wei, Chin. Phys. B 19, 107802 (2010) [CrossRef] [Google Scholar]
 D.H.G. Espinosa, C.L.P. Oliveira, A.M.F. Neto, J. Opt. Soc. Am. B 35, 346 (2018) [CrossRef] [Google Scholar]
 J.P. Huang, K.W. Yu, Phys. Rep. 431, 87 (2006) [CrossRef] [Google Scholar]
 K.L. Stokes, A. Puri, Opt. Lett. 15, 986 (1990) [CrossRef] [Google Scholar]
 J.W. Haus et al., J. Opt. Soc. Am. B 6, 797 (1989) [CrossRef] [Google Scholar]
 N. Kalyaniwalla, J.W. Haus, R. Inguva, M.H. Bimboim, Phys. Rev. A 42, 5613 (1990) [CrossRef] [PubMed] [Google Scholar]
 D.J. Bergman, O. Levy, D. Stround, Phys. Rev. B 49, 129 (1994) [CrossRef] [Google Scholar]
Cite this article as: Ya M. Wu, Guo Q. Chen, Wen C. Wu, Yi X. Wang, Yang Huang, Optical bistability in twodimensional nonlinear composites of coated cylinders with nonlinear core and graded shell, EPJ Appl. Metamat. 7, 6 (2020)
All Figures
Fig. 1 Schematic diagram of the proposed composite in which the coreshell inclusions have Kerrtype dielectric cores and graded Drudetype shells. 

In the text 
Fig. 2 (a) The local electric field as a function of the external field with different core radii. (b) The upper and lower thresholds as a function of core radius a with h = 0.9 (solid line) and h = 0 (dashed line), respectively. 

In the text 
Fig. 3 (a) The local electric field as a function of the external field with different frequency ω/ω _{p}. (b) The upper and lower thresholds as a function of frequency ω/ω _{p} with h = 0.9 (solid line) and h = 0 (dashed line), respectively. 

In the text 
Fig. 4 The width of bistable threshold as a function of a and ω/ω _{p}. 

In the text 
Fig. 5 The local field as a function of the external field with different host media. The gradation constant h = 0.9. 

In the text 
Fig. 6 The upper and lower thresholds as a function of host medium with h = 0.9 (solid) and h = 0 (dashed). 

In the text 
Fig. 7 The upper and lower thresholds as a function of h. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.