Abstract
We introduce two tractable analytical models to describe dynamic effects at resonant light scattering by subwavelength particles. One of them is based on a generalization of the temporal coupledmode theory, and the other employs the normal mode approach. We show that sharp variations in the envelope of the incident pulse may initiate unusual, counterintuitive dynamics of the scattering associated with interference of modes with fast and slow relaxation. To exhibit the power of the models, we apply them to explain the dynamic light scattering of a squareenvelope pulse by an infinite circular cylinder made of GaP, when the pulse carrier frequency lies in the vicinity of the destructive interference at the Fano resonances. We observe and explain intensive sharp spikes in scattering crosssections just behind the leading and trailing edges of the incident pulse. The latter occurs when the incident pulse is over and is explained by the electromagnetic energy released in the particle at the previous scattering stages. The accuracy of the models is checked against their comparison with results of the direct numerical integration of the complete set of Maxwell’s equations and occurs very high. The models’ advantages and disadvantages are revealed, and the ways to apply them to other types of dynamic resonant scattering are discussed.
1 Introduction
HighQ resonances are of utmost importance in a wide diversity of problems [1]. It is explained by the fact that to obtain strong resonant effects, the corresponding resonance should have a high amplitude, and hence a high Qfactor, at least if a spatially bounded system is a concern. However, the characteristic relaxation time for a resonance is inversely proportional to its Qfactor. Thus, the price one must pay for making use of highQ resonances is longlasting transient effects. On the other hand, the frontier of modern photonics moves toward short and ultrashort pulses. Nowadays, currently, these two factors together make a typical situation where the duration of the laser pulse becomes comparable or even shorter than the relaxation time of the resonance effects initiated by this pulse.
Though in quantum spectroscopy, it is wellknown that nonsteady resonant scattering may qualitatively differ from its steadystate realizations, see, e.g. Ref. [2–15], the corresponding studies in light scattering by subwavelength particles have begun only recently [16–18]. As it could be expected, these studies also reveal qualitatively new effects, which do not exist at the steadystate scattering. However, for the time being, analytical descriptions of the resonant light scattering by particles are still based on the solutions of Maxwell’s equations describing the steadystate scattering, where the processes of gradual “swinging” of resonant modes are not taken into account.
Thus, at the moment, the only possible theoretical description of transient effects at the resonant light scattering by subwavelength particles is made with the help of direct numerical integration of the complete set of Maxwell’s equations. Due to the existence of various numerical methods and many pieces of software, both free and commercial, created to perform this integration, it has become a more or less routine procedure. However, to find in the range of the problem parameters windows, where the desired effect is the most pronounced, the dependence of the scattering on these parameters in a wide domain of their variations is required. To this end, numerical methods are not appropriate, and analytically tractable models are required. To the best of our knowledge, for the time being, none of them exists.
Here, we present two such models, apply them to describe nonsteady Fano resonances, and compare the results with direct numerical integration of the complete set of Maxwell’s equations. The comparison indicates the high accuracy of both models and reveals their mutual advantages and disadvantages. The models unveil the physical grounds for the counterintuitive spikes in the scattered radiation observed behind the leading and trailing edges of the incident pulse.
The first model is based on the temporal coupledmode theory (TCMT) [19] generalized to applications to essentially nonsteady scattering. The second model mimics nonsteady resonant vibrations by a superposition of dynamics of driven harmonic oscillators (HO). The latter approach looks similar to the harmonic inversion, see, e.g. Ref. [20–22]. However, in our case, the mode selection for the approximation and, most importantly, the choice of the values of the model parameters are based on the system in question’s physical properties. Thus, it reveals the role of different excitations in the system dynamic and sheds light on the physical nature of the system as a whole. Besides, this makes it possible to reduce the number of modes to be studied just to a few with nontrivial dynamics. In addition to the purely academic interest, the results obtained may be employed to design a new generation of fast, multifunctional passive nanodevices whose properties vary quantitatively with a variation of the duration of the incident pulse.
The paper has the following structure: in Section 2, the problem formulation is presented. Section 3 is devoted to the problem analysis and discussion of the obtained results. In Section 4, we formulate conclusions. Cumbersome details and expressions and several plots illustrating the developed approach are moved to Appendix.
2 Problem formulation
2.1 Preliminary
The simplest exactly solvable light scattering problems correspond to the steadystate scattering of a linearly polarized plane wave by a homogeneous sphere (the Mie solution) or infinite circular cylinder. In these solutions, the scattered field is presented as an infinite series of partial waves (dipole, quadrupole, etc.), also called multipoles [23]. Here, we generalize this approach to dynamic light scattering, making it possible to study quite intricate transient features.
One of the most typical resonant responses in light scattering by finite obstacles is the Fano resonance [24]. A characterizing it asymmetric lineshape is explained by either constructive (the maximal scattering) or destructive (the minimal scattering) interference occurring close to each other in the frequency domain. The interfering parties are the socalled resonant (narrow line) and background (broad line) partitions of the same multipole. The narrowline and broadline partitions correspond to excitations with slow and fast relaxation time in the timedomain, respectively.
Naturally, the procedure of splitting a single partial wave into the two partitions is not unique. Accordingly, there are two main equivalent approaches to it. In the first approach, a partial wave is presented as a sum of the radiation of the conventional electric and toroidal multipoles [25]. However, in what follows, it will be more convenient for us to employ another approach [26]. In this approach, the resonant partition is associated with the corresponding electromagnetic mode excited in the bulk of the particle (volume polariton). In contrast, the background partition corresponds to the radiation of the surface current induced by the same incident wave, scattered by the particle with the same geometrical shape but made of a hypothetical material called the perfect electric conductor (PEC).
At a point of destructive interference, the two partitions cancel each other. As a result, the contribution of the corresponding multipole to the stationary scattering is totally suppressed. For a subwavelength particle, when just a few first multipoles produce the overwhelming contribution to the overall scattering, the suppression even of one of them may reduce the scattering crosssection dramatically [26]. However, since, as has been mentioned above, the resonant and background partitions are characterized by substantially different relaxation times, it is evident that during transient regimes, the mutual cancelation does not occur. Therefore, the violation of the destructive interference conditions must give rise to a considerable increase in the scattering intensity during the transient and other unusual effects. Some of them have already been discussed in our previous publications [16, 17].
We stress that though this increase in the scattering intensity looks similar to the wellknown overshoot effect when a driven highQ oscillator exhibits oscillatory relaxation to the steadystate, the physical grounds for the former is entirely different: If the overshoot is related to vibrations of a single oscillator, the discussed effect is explained by a superposition of two different oscillations. This distinction in the physical nature of the two cases gives rise to the corresponding difference in their manifestations. In particular, in contrast to the overshoot, the transient at the point of the destructive interference corresponds to relaxation to the steadystate with zero amplitude of the given multipole. Therefore, for a small particle with just a few dominant multipoles, the amplitude of the spike in the overall scattering may be in orders of magnitude larger than that exhibited by the overshoot.
Thus, the dynamic resonant light scattering is a new, rich, and practically untouched subfield. Plenty of exciting effects hidden there are still undiscovered. In the present paper, we continue to explore this appealing topic.
To understand typical main features of the phenomenon, we consider the simplest problem formulation corresponding to the one employed in Ref. [16] for the computer simulation, namely the scattering of a square incident pulse with duration τ (the amplitudes of the electromagnetic fields A(t) = A _{0} = const ≠ 0 inside the pulse and zero outside it), carrier frequency ω, and temporal dependence of the fields ∼exp(−iωt) by an infinite circular highindex cylinder with the base radius R and complex refractive index m = n + iκ (n ≫ 1, κ ≪ 1). The refractive index of the surrounding cylinder medium equals unity. The cylinder is nonmagnetic, so its permeability μ = 1. For further simplification, just the TE polarization and normal incidence are considered, see below Figure 1a.
2.2 Instantaneous scattering crosssection
In the conventional steadystate case, the scattering is quantitatively described by the crosssection C _{sca} calculated per unit of length of the cylinder. C _{sca} is defined as the ratio of the integral power flux through a closed remote surface surrounding the scatterer to the intensity of the incident light. At nonsteady scattering, both these quantities are timedependent, and their ratio is not a constant anymore. Moreover, the ratio depends on the shape and position of the surface used to calculate the flux since the speed of light is finite.
To describe nonsteady scattering, we introduce the instantaneous scattering crosssection C
_{sca}(t) as the ratio of the instantaneous value of the power flux through a cylindrical surface, coaxial with the scattering cylinder and lying in the far wave zone, calculated per a unit of length of the cylinder, to the constant intensity of the square incident pulse I
_{0} [16]. In a more general case, when the pulse envelope has a timedependent shape I(t), the scattered flux may be normalized over, e.g. the maximal value of I(t), i.e.
The corresponding dimensionless scattering efficiency, Q _{sca}(t) is connected with C _{sca}(t) by the usual relation Q _{sca}(t) = C _{sca}(t)/(2R), where R is the radius of the base of the cylinder. We also do not perform the time averaging of the Poynting vector over the period of oscillations.
The exact solution describing the steadystate scattering is built as an infinite series of partial waves (multipoles). For the problem in question, the complex amplitudes of the multipoles (scattering coefficients) associated with the outgoing partial waves a _{ ℓ } and the field within the cylinder d _{ ℓ } are given by the wellknown expressions presented in Appendix Eqs. (A.1) and (A.2), see also, e.g. [23]. Here, ℓ designates the multipole order (dipole, quadrupole, etc.). Note that a _{ ℓ } = a _{−ℓ }, d _{ ℓ } = d _{−ℓ } [23].
For the given problem, the scattering coefficients satisfy the identity [26]:
where x = kR stands for the size parameter; k = ω/c; c is the speed of light in a vacuum; J
_{
ℓ
}(z) and
Routine calculations result in the following expressions [16]:
Here,
2.3 Fano resonances
The Fano resonances [24, 27, 28] are a good example demonstrating unusual, counterintuitive effects in transient processes of resonant light scattering. For the steadystate scattering, their detailed discussion is presented, e.g. in Ref. [26]. Though in that paper, a spherical particle is considered, generalization to the case of a cylinder is a straightforward matter.
As it has been mentioned above, a key point of the Fano resonances is a presentation of the scattered wave as a sum of two partitions: resonant and background. In the proximity of the resonance, the amplitude and phase of the former have a sharp ωdependence, while for the latter its ωdependence is weak. An important conclusion following from the results of Ref. [26, 29] is that the splitting of a
_{
ℓ
} into the two terms, given by Eq. (1), actually, is the singling out the background (
Our goal is to recover the full timedependence Q _{sca}(t). For highQ resonances, which we are interested in, the characteristic time scales of the transients should be large relative to the period of the field oscillations 2π/ω. Then, a quasisteady approximation may be employed. It implies the same structure of the solution like that for steadystate scattering. However, now the scattering coefficients are regarded as slowly varying functions of time. Let us apply this assumption to the TCMT and HO.
2.4 TCMT
In the specified case, the TCMT equations read as follows, see, e.g. Ref. [30]:
Here, s
^{+}(t) and s
^{−}(t) are the amplitudes of the incoming (converging) and outgoing (diverging) cylindrical waves, respectively; B is the background reflection coefficient; κ and ζ are coupling constants; p(t) describes the internal resonant mode excitation; and
The analysis of Eqs. (4) and (5) performed in Ref. [30] for the steadystate scattering indicates that
Meanwhile, there are other ways to fix ϕ, free from this disadvantage. In this paper, we fix ϕ from the condition
The next difficulty is that the employed expression for
Regarding Q
_{sca}(t), in the discussed quasisteady approximation the solution retains the same structure as that for the steadystate. Then, Eqs. (2) and (3) still remain valid but the replacement
at 0 ≤ t ≤ τ and
at t > τ (remember that γ < 0). Note, that
Another point to be stressed is that
2.5 Harmonic oscillators
Another model is based on the wellknown fact that any linear oscillatory dynamic may be approximated by that of a system of driven coupled HO. We just have to apply it to the problem in question. The general solution of the equations for driven coupled HO has the form (see, e.g. [31])
where z
_{
k
}(t) is a complex coordinate of the kth oscillator;
The main idea of the adaptation of Eq. (8) for a drive with a carrier frequency ω is that only the dynamics of the resonant eigenmodes with the frequency mismatch
Next, according to what just has been mentioned, for the problem in question the eigenmode dynamic is described by the equation:
supplemented with the initial conditions
Equation (9) has the following exact solution:
at 0 ≤ t ≤ τ and
at t > τ. Here,
In what follows, we employ Eqs. (10) and (11) to model the dynamics of the scattering coefficients. Naturally, for different coefficients, the values of A
_{0}, ω
_{0}, and γ are also different and will be defined below. The key equation now is Eq. (1), where the steadystate
Let us stress the dramatic difference in the steadystate profiles d
_{
ℓ
}(x) and
To employ Eq. (9) for modeling the dynamics of the eigenmodes, we have to fix the values of the following four parameters: Re[A _{0}], Im[A _{0}], ω _{0}, and γ. The approximation details are as follows: For the resonant dmodes, the values of ω _{0} and γ are given by the corresponding poles of d _{ ℓ }(ω) in the same manner as that in TCMT. Then, we require that at the drive frequency ω, the complex amplitude of the oscillator coincides with that for d _{ ℓ }(ω) of the exact solution.
The profile of a PECmode is far from a narrow resonant line. Therefore, the poles in the complex plane might have nothing to do with the dynamic of the mode. Then, ω _{0} and γ become free parameters, and we need two more conditions to complete the problem. For them, we select (i) the equality of the frequencies maximizing the profile a ^{(PEC)}(ω) and the corresponding profile of the oscillator and (ii) the equality of the maxima themselves. This procedure fixes all four parameters of the HO model unambiguously. For more details, see Appendix.
3 Results and discussions
3.1 Parameters, variables, and plots
To illustrate the accuracy of the two models, their advantages and disadvantages, we perform computer simulations of the scattering based on the direct numerical integration of the complete set of Maxwell’s equations with the help of Lumerical’s FDTD SOLUTIONS and compare their results with the dynamics described by the models. As it has been mentioned above, the mutual orientation of the cylinder and the incident wave corresponds to that shown in Figure 1a. It is convenient to transfer to the dimensionless time: t _{new} = t _{old} c/R. Then, ω _{new} = ω _{old} R/c ≡ x. Since below only the dimensionless quantities are in use, the subscript “new” will be dropped.
The refractive index is selected purely real and equal to m = n = 3.125, which corresponds to that of gallium phosphide (GaP) irradiated in a vacuum at the wavelength λ
^{(GaP)} = 966.4 nm [32]. For the size parameter the value x = 1.702 is selected. At the given wavelength, it implies that the radius of the cylinder equals R = 261.8 nm. The choice is done since this pair of m and x corresponds to a local minimum of
For all other multipoles, x = x _{min} corresponds to the offresonant regions. Thus, only the dynamics of the modes with ℓ = 0, ±2 should be approximated. Moreover, since the scattering coefficients differing only by the sign of ℓ are identical, the modes with ℓ = ±2 may be regarded as a single one.
The pulse duration, τ equals 191.28. This is much larger than any problem’s characteristic relaxation time at the given values of the other parameters. Then, it makes it possible to study both: the transient to the stationary scattering (behind the leading edge of the incident pulse) and its decay (behind the trailing edge).
The poles of the scattering coefficients adjacent to ω = ω _{min} are
Then, the described above procedure gives rise to the following values of the models’ parameters (see Appendix for details):
For TCMT:
For HO:
dmodes:
PECmodes:
Note that, despite the error in the approximation of a ^{(PEC)}, the approximation of a _{ ℓ }(ω) in the vicinity of the Fano resonances is quite accurate, cf. Figures 8 and 9.
The quantitative comparison of the two models with the numerics at the dimensionless ω = ω _{min}, shown in Figure 2, exhibits their high accuracy. The origin of the taxis in Figure 2 (and below, in Figures 3 and 4) is shifted so that t = 0 corresponds to the moment when the scattered radiation for the first time is detected by the measuring monitors, situated in the far wave zone.
To demonstrate that the developed models are applicable to describe with high accuracy the dynamic scattering not only in the vicinity of the minimum of the Fano line but also within the entire Fano resonance line, as a whole, we also show in Figure 2, their comparison with the numerics at the point of the local maximum of the steadystate scattering (ω = ω _{max} = 1.525; λ ^{(GaP)} = 1079 nm) and in the middle of the line (ω = ω _{mid} = 1.589; λ ^{(GaP)} = 1035 nm), see Figure 1b; the meaning of the subscripts min and max is connected with the scattering intensity, therefore, please do not be confused by the fact that numerical value of ω _{max} is smaller than that for ω _{min}.
The number of dominant modes in these cases remains the same as that at ω = ω _{min}. Therefore, the poles and resonant conditions of each mode do not change too. The fitting procedure remains identical to that described above, though, of course, the matching conditions should be implemented for the corresponding new values of the driver frequency.
Now, let us discuss the most interesting, counterintuitive effect in the scattering dynamics, namely the sharp intensive spikes just behind the leading (precursor) and trailing (postcursor) edges of the incident pulse, observed at ω = ω _{min} both in the total and directional scattering, see Figures 2(a) and 4. Note that while the pre and postcursor look similar, their nature is completely different. The precursor is related to the fast excitation of the PECmodes with the low Qfactors. It lasts for the period while the relatively slow resonant modes are not “swung” yet. Since the line of any PECmode is broad, the variation of the carrier frequency of the incident pulse in the range 1.525 ≤ ω ≤ 1.702 affects the value of the PECline profile very little, see Figure 8. Accordingly, it weakly affects the precursor dynamics; see Figure 3, where we show a set of the precursor lines in the discussed range of ω variations with the step Δω = 0.00885. For all curves shown in Figure 3, the maximal value of the precursor is about 1.6 and its shape does not change much.
In contrast, the postcursor is explained by the violation of the balance between the resonant and PECmodes caused by the rapid decay of the latter, when the pumping with the incident pulse is over, and electromagnetic energy of the dark, nonradiating in the steadystate scattering modes, begins to leak out. Thus, the maximal value of the postcursor is bounded by the amplitude of the resonant modes just behind the trailing edge of the incident pulse. Due to the narrowness of the resonant lines, the same variation of ω (1.525 ≤ ω ≤ 1.702) affects the steadystate amplitude of the resonant modes substantially, see Figure 7. Accordingly, at the point of the destructive interference at ω = 1.702, where the steadystate amplitudes of the resonant and background partitions are almost equal, the maximal values of the pre and postcursor are almost equal too. Regarding the duration of the pre and postcursor, both are determined by the relaxation times of the resonant modes (i.e. by 1/γ), which, therefore, are the same for the pre and postcursor. As a result, the shapes of the pre and postcursor in Figure 2a are very similar to each other.
However, a departure from ω = 1.702 toward ω = 1.525 increases the steadystate amplitude of the resonant partition, while the one for the background remains practically unchanged. This gives rise to a decrease in the ratio of the maximal value of the postcursor to the amplitude of the steadystate scattering. Eventually, at the local maximum of the scattering, the “basement” of the stationary scattering eats out the postcursor almost entirely, cf. Figures 2a and c.
It is important to note that both developed models make it possible to describe the dynamic of each multipole independently. Therefore, in addition to the reproduction of the total dynamic scattering, the models can be equally applied to describe the directional dynamic scattering too. The latter may be more practical from the viewpoint of getting experimental evidence of the discussed effects.
Qualitatively, the directional scattering occurs similar to the overall one. In particular, in the vicinity of ω = ω
_{min} = 1.702 (
3.2 RMSE
To get a quantitative criterion of the accuracy of the models in a given window t _{1} ≤ t ≤ t _{2} we do the following: The time in the window is sampled with a step δt. Then, the root mean square error (RMSE) is defined as follows
is calculated (do not mix dummy n here with the real part of the complex refractive index m). Here,
For the two windows with essentially nonsteady dynamics adjacent to the leading and trailing edges of the incident pulse (t _{1} = 0, t _{2} = 38 and t _{1} = τ, t _{2} = τ + 38),^{[1]} the calculated errors are summarized in Table 1 for the total dynamic scattering and in Table 2 for the directional one. The comparison of the results for the precursor RMSE^{(TCMT)} = 0.181, RMSE^{(HO)} = 0.144, and the postcursor RMSE^{(TCMT)} = 0.166, RMSE^{(HO)} = 0.114, in the vicinity of the destructive Fano interference conditions (ω = ω _{min} = 1.702) illustrates that the HO model is more accurate. This is expected since the dynamics of each mode is modeled by the secondorder equation instead of the first one (like in the TCMT case), allowing to capture the transient dynamics in greater detail.
RMSE for total scattering cross section  

HO  TCMT  HO  TCMT  
ω = ω _{min}  0.144  0.181  0.114  0.166 
ω = ω _{mid}  0.282  0.211  0.192  0.273 
ω = ω _{max}  0.295  0.252  0.298  0.212 
Precursor  Postcursor 
RMSE for directional 90° scattering  

HO  TCMT  HO  TCMT  
Forward  0.0295  0.0261  0.0265  0.0285 
Backward  0.0291  0.0236  0.0188  0.0189 
90°  0.0066  0.0130  0.0063  0.0144 
Precursor  Postcursor 
4 Conclusions
Thus, though both models demonstrate high accuracy, the one of the HO is a little better than that of TCMT. In addition, while the TCMTmodel just describes the complicated dynamics of the resonant light scattering, the HO model also elucidates the physical grounds for the sharp intensive spikes in the scattering observed behind the leading and trailing edges of the incident pulse and explained by the violation of the balance between the resonant and background excitations. The violation, in turn, is coursed by the difference in the corresponding relaxation times. Then, the uncompensated part of the interfering modes is exhibited as a spike.
Note that if the “basement” of the scattered pulse is cut off at the level above the intensity of the steadystate scattering, the duration of the remaining parts of the spikes occurs very small. This procedure may be used as a new method of radiation at the nanoscale of sharply directed pulses with duration of a few periods of the field oscillations.
We have also shown that both models can be successfully applied to describe the total and directional dynamic scattering at an arbitrary wavelength. On top of that, importantly, the models also reveal the underlying physical phenomena associated with the fast and slow modes excitation, whose interference leads to new effects, including the formation of the sharp spikes.
Regarding possible extensions of the models, while the governing equations of TCMT are simpler than those of HO, the accuracy of the HO model is higher relative to the one of TCMT. Besides, the implementation of the HO model is more straightforward and does not require a sophisticated procedure to connect the model’s parameters with those of the initial underlying problem.
If just a single resonant excitation is a concern, the HO model may be used based even on fitting an experimentally obtained spectrum of the output signal modules. On the other hand, when the interference of several excitations is essential, its application implies knowledge of the phases, obtained either from an analytical solution or measured experimentally. While the latter is challenging in the optical range, it is a routine procedure at radio frequencies; see, e.g. [33, 34] and references therein. Thus, both models complement each other and may be very useful in descriptions of a wide variety of resonant phenomena.
Funding source: Ministry of Education and Science of the Russian Federation
Award Identifier / Grant number: 02.a03.21.0005
Funding source: Russian Science Foundation
Award Identifier / Grant number: 197230012
Award Identifier / Grant number: 211200151
Funding source: Russian Foundation for Basic Research
Award Identifier / Grant number: 200200086
Funding source: Australian Research Council
Award Identifier / Grant number: DP200101353
Acknowledgment
The authors are very grateful to Boris Y. Rubinstein for his valuable help in symbolic computer calculations.

Author contribution: All the authors have accepted responsibility for the entire content of this submitted manuscript and approved submission.

Research funding: M.I.T. acknowledges the financial support of the Russian Foundation for Basic Research (Projects No. 200200086) for the analytical study, the Moscow Engineering Physics Institute Academic Excellence Project (agreement with the Ministry of Education and Science of the Russian Federation of 27 August 2013, Project No. 02.a03.21.0005) for the modeling of the resonant light scattering, as well as the contribution of the Russian Science Foundation for the computer simulation (Project No. 211200151) and the provision of user facilities (Project No. 197230012).

Conflict of interest statement: The authors declare no conflicts of interest regarding this article.
A.1 Scattering coefficients
The scattering coefficient for the problem under consideration read as follows [23]:
A.2 The approximation procedure
In what follows ω _{sim} designates the carrier frequency of the incident pulse corresponding to a given simulation, while all specific examples are given at ω _{sim} = ω _{min} = 1.702. The application procedure at other values of ω _{sim} discussed in the main text is identical to that at ω _{sim} = ω _{min}.
A.2.1 TCMTmodel
To obtain the value of ϕ, note that according to Ref. [30]
(the same expression may be obtained from Eq. (6), if we formally consider its limit at t → ∞). Then, after some algebra Eq. (A.3) gives rise to the following expression for
where ϵ = −(ω − ω _{0})/γ and q = − cot(ϕ/2). We remind that γ < 0. The last formula in Eq. (A.4) is the conventional Fano profile with the asymmetry parameter q [24, 27, 28].
Following the procedure described in the main text, we have to equalize
where ϵ _{sim} = −(ω _{sim} − ω _{0})/γ. Note that for the problem in question a _{ ℓ } is always smaller than unity, see, e.g. [35]. Therefore, the roots in Eq. (A.5) are always real.
Plots a _{ ℓ }(ω), as well as their approximations at ω _{sim} = ω _{min}, q = q _{(+)}, q = q _{(−)}, and ℓ = 0, 2 for the values of ω _{0}, γ given by Eq. (12), are shown in Figures 5 and 6. Thus, for ℓ = 0, the better overall approximation is at q = q _{(−)}. For the case in question its numerical value is 0.4700, which corresponds to ϕ ≈ −2.263. In contrast, for ℓ = 2, the better overall approximation is at q = q _{(+)} ≈ −4.172, corresponding to ϕ ≈ 0.471.
A.2.2 HO model
A steadystate solution of Eq. (9) with the r.h.s equals A _{0} exp[−iωt] is
Thus, to apply the HO model to the discussed light scattering problem, the following parameters should be fixed: the complex drive amplitude A _{0} ≡ A _{0}e^{iφ }, the frequency ω _{0}, and the damping factor γ. In the case of the dmodes, the values of ω _{0} and γ are given by the position of the nearest to ω _{sim} complex pole of the steadystate coefficient d _{ ℓ }(ω). The remaining undefined parameters are A _{0} and φ. They are readily obtained from the equality F _{ s }(ω _{sim}) = d _{ ℓ }(ω _{sim}). The results of this approximation at ω _{sim} = ω _{min} and the values of the other parameters specified in the main text are shown in Figure 7.
The case of the PECmodes is more tricky since neither ω _{0}, nor γ can be obtained in the same easy manner as that for the dmodes. Once again, it is convenient to present the complex amplitude A _{0} in the form A _{0} = A _{0}e^{iφ }. Then, the maximum of F _{ s }(ω)^{2} equals
and is achieved at
We remind that
The condition
where
Now, the last unfixed parameter φ is readily obtained from the condition that the phase of
The application of this procedure results in the values of the parameters presented in the main text, see Eqs. (15) and (16). Note that while the approximation errors for the PEC modes are large, the final approximations for a _{ ℓ } in the proximity of ω _{sim} are quite accurate, cf. Figures 8 and 9.
Next, there are some points, where
Another point of a special interest is the one, when ω
_{sim} = ω
_{max}. In this case the straightforward application of Eq. (A.9) gives rise to an indeterminate form 0/0. As usual, it means that we have to consider the limit ω
_{sim} → ω
_{max}. Expanding in Eq. (A.9)
where
References
[1] L. Huang, L. Xu, M. Rahmani, D. Neshev, and A. E. Miroshnichenko, “Pushing the limit of highq mode of a single dielectric nanocavity,” Adv. Photonics, vol. 3, p. 016004, 2021. https://doi.org/10.1117/1.ap.3.1.016004.Search in Google Scholar
[2] A. Kaldun, A. Blättermann, V. Stooß, et al.., “Observing the ultrafast buildup of a fano resonance in the time domain,” Science, vol. 354, no. 6313, pp. 738–741, 2016. https://doi.org/10.1126/science.aah6972.Search in Google Scholar
[3] O. Katz, E. Small, Y. Bromberg, and Y. Silberberg, “Focusing and compression of ultrashort pulses through scattering media,” Nat. Photonics, vol. 5, no. 6, pp. 372–377, 2011. https://doi.org/10.1038/nphoton.2011.72.Search in Google Scholar
[4] P. M. Johnson, A. Imhof, B. P. J. Bret, J. G. Rivas, and A. Lagendijk, “Timeresolved pulse propagation in a strongly scattering material,” Phys. Rev. E, vol. 68, p. 016604, 2003. https://doi.org/10.1103/PhysRevE.68.016604.Search in Google Scholar
[5] G. Chaussonnet and P. M. Bardet, “Scattering of an ultrashort laser pulse by a spherical air bubble,” Opt. Express, vol. 28, pp. 27358–27371, 2020. https://doi.org/10.1364/oe.401010.Search in Google Scholar
[6] D. Q. Chowdhury, S. C. Hill, and P. W. Barber, “Time dependence of internal intensity of a dielectric sphere on and near resonance,” J. Opt. Soc. Am. A, vol. 9, pp. 1364–1373, 1992. https://doi.org/10.1364/josaa.9.001364.Search in Google Scholar
[7] E. E. M. Khaled, D. Q. Chowdhury, S. C. Hill, and P. W. Barber, “Internal and scattered timedependent intensity of a dielectric sphere illuminated with a pulsed Gaussian beam,” J. Opt. Soc. Am. A, vol. 11, pp. 2065–2071, 1994. https://doi.org/10.1364/josaa.11.002065.Search in Google Scholar
[8] K. S. Shifrin and l. G. Zolotov, “Quasistationary scattering of electromagnetic pulses by spherical particles,” Appl. Opt., vol. 33, pp. 7798–7804, 1994. https://doi.org/10.1364/ao.33.007798.Search in Google Scholar
[9] L. Mees, G. Gréhan, and G. Gouesbet, “Timeresolved scattering diagrams for a sphere illuminated by plane wave and focused short pulses,” Opt. Commun., vol. 194, no. 1, pp. 59–65, 2001. https://doi.org/10.1016/s00304018(01)012688.Search in Google Scholar
[10] J. A. Lock and P. Laven, “Mie scattering in the time domain. part 1. the role of surface waves,” J. Opt. Soc. Am. A, vol. 28, pp. 1086–1095, 2011. https://doi.org/10.1364/josaa.28.001086.Search in Google Scholar PubMed
[11] P. Tomchuk and D. Butenko, “Single and double ultrashort laser pulse scattering by spheroidal metallic nanoparticles,” J. Nanophotonics, vol. 10, no. 1, pp. 1–14, 2016. https://doi.org/10.1117/1.jnp.10.016018.Search in Google Scholar
[12] V. A. Astapenko, F. B. Rosmej, and E. S. Khramov, “Time dependence of ultrashort laser pulses scattering by atom in high frequency limit,” Atoms, vol. 8, no. 3, p. 41, 2020. https://doi.org/10.3390/atoms8030041.Search in Google Scholar
[13] P. J. Ho, A. E. A. Fouda, K. Li, G. Doumy, and L. Young, “Ultraintense, ultrashort pulse xray scattering in small molecules,” Faraday Discuss., vol. 228, pp. 139–160, 2021. https://doi.org/10.1039/d0fd00106f.Search in Google Scholar PubMed
[14] V. A. Astapenko and N. N. Moroz, “Elastic scattering of ultrashort laser pulses on noble gas atoms,” Russ. Phys. J., vol. 61, no. 7, pp. 1329–1337, 2018. https://doi.org/10.1007/s1118201815388.Search in Google Scholar
[15] D. N. Makarov, “Quantum theory of scattering of ultrashort electromagnetic field pulses by polyatomic structures,” Opt. Express, vol. 27, pp. 31989–32008, 2019. https://doi.org/10.1364/oe.27.031989.Search in Google Scholar
[16] M. I. Tribelsky and A. E. Miroshnichenko, “Dynamics of destructive fano resonances,” Phys. Rev. A, vol. 100, p. 053824, 2019. https://doi.org/10.1103/physreva.100.053824.Search in Google Scholar
[17] S. E. Svyakhovskiy, V. V. Ternovski, and M. I. Tribelsky, “Anapole: its birth, life, and death,” Opt. Express, vol. 27, pp. 23894–23904, 2019. https://doi.org/10.1364/oe.27.023894.Search in Google Scholar PubMed
[18] O. ÁvalosOvando, L. V. Besteiro, Z. Wang, and A. O. Govorov, “Temporal plasmonics: Fano and rabi regimes in the time domain in metal nanostructures,” Nanophotonics, vol. 9, no. 11, pp. 3587–3595, 2020. https://doi.org/10.1515/nanoph20200229.Search in Google Scholar
[19] W. H. Louisell, Coupled Mode and Parametric Electronics, New York and London, Wiley, 1960.Search in Google Scholar
[20] V. A. Mandelshtam and H. S. Taylor, “Harmonic inversion of time signals and its applications,” J. Chem. Phys., vol. 107, no. 17, pp. 6756–6769, 1997. https://doi.org/10.1063/1.475324.Search in Google Scholar
[21] P. Barone, E. Massaro, and A. Polichetti, “The segmented prony method for the analysis of nonstationary time series,” Astron. Astrophys., vol. 209, pp. 435–444, 1989.Search in Google Scholar
[22] A. Roessling and J. Ringwood, “Finite order approximations to radiation forces for wave energy applications,” in Renewable Energies Offshore, vol. 359, London, Taylor and Francis, 2015, pp. 359–366.10.1201/b1897351Search in Google Scholar
[23] C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles, New York, WileyVCH Verlag, 1998.10.1002/9783527618156Search in Google Scholar
[24] A. E. Miroshnichenko, S. Flach, and Y. S. Kivshar, “Fano resonances in nanoscale structures,” Rev. Mod. Phys., vol. 82, no. 3, pp. 2257–2298, 2010. https://doi.org/10.1103/revmodphys.82.2257.Search in Google Scholar
[25] A. E. Miroshnichenko, A. B. Evlyukhin, Y. F. Yu, et al.., “Nonradiating anapole modes in dielectric nanoparticles,” Nat. Commun., vol. 6, no. 1, p. 8069, 2015. https://doi.org/10.1038/ncomms9069.Search in Google Scholar PubMed PubMed Central
[26] M. I. Tribelsky and A. E. Miroshnichenko, “Giant inparticle field concentration and fano resonances at light scattering by highrefractiveindex particles,” Phys. Rev., vol. 93, p. 053837, 2016. https://doi.org/10.1103/physreva.93.053837.Search in Google Scholar
[27] U. Fano, “On the absorption spectrum of noble gases at the arc spectrum limit,” Nuovo Cimento, vol. 12, pp. 154–161, 1935. https://doi.org/10.1007/bf02958288. Available at: http://arXiv.org/abs/condmat/0502210v1.Search in Google Scholar
[28] U. Fano, “Effects of configuration interaction on intensities and phase shifts,” Phys. Rev., vol. 124, pp. 1866–1878, 1961. https://doi.org/10.1103/physrev.124.1866.Search in Google Scholar
[29] M. V. Rybin, K. B. Samusev, I. S. Sinev, et al.., “Mie scattering as a cascade of Fano resonances,” Opt. Express, vol. 21, no. 24, p. 30107, 2013. https://doi.org/10.1364/oe.21.030107.Search in Google Scholar PubMed
[30] Z. Ruan and S. Fan, “Temporal coupledmode theory for Fano resonance in light scattering by a single,” J. Phys. Chem. C, vol. 114, no. 16, pp. 7324–7329, 2010. https://doi.org/10.1021/jp9089722.Search in Google Scholar
[31] L. Landau and E. Lifshitz, Mechanics: Volume 1 (Course of Theoretical Physics Series), vol. 23, Oxford, Oxford Pergamon Press, 2000.Search in Google Scholar
[32] M. Polyanskiy, Refractive Index Database. Available at: http://refractiveindex.info/.Search in Google Scholar
[33] M. I. Tribelsky, J.M. Geffrin, A. Litman, C. Eyraud, and F. Moreno, “Small dielectric spheres with high refractive index as new multifunctional elements for optical devices,” Sci. Rep., vol. 5, no. 1, p. 12288, 2015. https://doi.org/10.1038/srep12288.Search in Google Scholar PubMed PubMed Central
[34] M. I. Tribelsky, J.M. Geffrin, A. Litman, C. Eyraud, and F. Moreno, “Directional fano resonances in light scattering by a high refractive index dielectric sphere,” Phys. Rev. B, vol. 94, p. 121110, 2016. https://doi.org/10.1103/physrevb.94.121110.Search in Google Scholar
[35] M. I. Tribelsky, “Phenomenological approach to light scattering by small particles and directional fano’s resonances,” EPL, vol. 104, p. 34002, 2013. https://doi.org/10.1209/02955075/104/34002.Search in Google Scholar
© 2021 Michael I. Tribelsky and Andrey E. Miroshnichenko, published by De Gruyter, Berlin/Boston
This work is licensed under the Creative Commons Attribution 4.0 International License.