Self-consistent numerical modeling of radiatively damped Lorentz oscillators
AbstractRecent progress towards realizing quantum emitters (QEs) suitable for integration in quantum information networks stimulates the demand for a self-consistent numerical approach to describe scattering of radiatively limited QEs in complex dielectric environments. As the QE response has to be characterized without the use of phenomenological damping parameters, the divergent nature of the pointlike emitter's in-phase self-field has to be carefully dealt with to avoid unphysical frequency shifts. Here we provide a solution to this problem and show two ways to obtain accurate results in the weak excitation limit using finite-difference time-domain algorithms. One of these approaches lays important groundwork needed for future simulations of nonlinear QE networks. The solution for dealing with the frequency shift reveals that dynamical contributions to the resonant depolarization field of arbitrarily small dielectric objects make crucial contributions to the net dipole moment induced by an external field when radiative scattering is the only source of damping.