diff --git a/Electrodynamics.bib b/Electrodynamics.bib new file mode 100644 index 0000000..4961e18 --- /dev/null +++ b/Electrodynamics.bib @@ -0,0 +1,451 @@ + +@article{sauvan_theory_2013, + title = {Theory of the {Spontaneous} {Optical} {Emission} of {Nanosize} {Photonic} and {Plasmon} {Resonators}}, + volume = {110}, + url = {http://link.aps.org/doi/10.1103/PhysRevLett.110.237401}, + doi = {10.1103/PhysRevLett.110.237401}, + abstract = {We provide a self-consistent electromagnetic theory of the coupling between dipole emitters and dissipative nanoresonators. The theory that relies on the concept of quasinormal modes with complex frequencies provides an accurate closed-form expression for the electromagnetic local density of states of any photonic or plasmonic resonator with strong radiation leakage, absorption, and material dispersion. It represents a powerful tool to calculate and conceptualize the electromagnetic response of systems that are governed by a small number of resonance modes. We use the formalism to revisit Purcell’s factor. The new formula substantially differs from the usual one; in particular, it predicts that a spectral detuning between the emitter and the resonance does not necessarily result in a Lorentzian response in the presence of dissipation. Comparisons with fully vectorial numerical calculations for plasmonic nanoresonators made of gold nanorods evidence the high accuracy of the predictions achieved by our semianalytical treatment.}, + number = {23}, + urldate = {2015-04-12}, + journal = {Physical Review Letters}, + author = {Sauvan, C. and Hugonin, J. P. and Maksymov, I. S. and Lalanne, P.}, + month = jun, + year = {2013}, + pages = {237401}, + annote = {Definition of mode volume in lossy systems, cf. eqs (9), (10) and the text around.}, + file = {APS Snapshot:/u/46/necadam1/unix/.zotero/zotero/9uf64zmd.default/zotero/storage/Q4P33DFT/PhysRevLett.110.html:text/html;PhysRevLett.110.237401.pdf:/u/46/necadam1/unix/.zotero/zotero/9uf64zmd.default/zotero/storage/P6ZH92PB/PhysRevLett.110.237401.pdf:application/pdf;Sauvan_supplementary_material.pdf:/u/46/necadam1/unix/.zotero/zotero/9uf64zmd.default/zotero/storage/SSCXRFQ8/Sauvan_supplementary_material.pdf:application/pdf} +} + +@article{pustovit_plasmon-mediated_2010, + title = {Plasmon-mediated superradiance near metal nanostructures}, + volume = {82}, + url = {http://link.aps.org/doi/10.1103/PhysRevB.82.075429}, + doi = {10.1103/PhysRevB.82.075429}, + abstract = {We develop a theory of cooperative emission of light by an ensemble of emitters, such as fluorescing molecules or semiconductor quantum dots, located near a metal nanostructure supporting surface plasmon. The primary mechanism of cooperative emission in such systems is resonant energy transfer between emitters and plasmons rather than the Dicke radiative coupling between emitters. We identify two types of plasmonic coupling between the emitters, (i) plasmon-enhanced radiative coupling and (ii) plasmon-assisted nonradiative energy transfer, the competition between them governing the structure of system eigenstates. Specifically, when emitters are removed by more than several nanometers from the metal surface, the emission is dominated by three superradiant states with the same quantum yield as a single emitter, resulting in a drastic reduction of ensemble radiated energy, while at smaller distances cooperative behavior is destroyed by nonradiative transitions. The crossover between two regimes can be observed in distance dependence of ensemble quantum efficiency. Our numerical calculations incorporating direct and plasmon-assisted interactions between the emitters indicate that they do not destroy the plasmonic Dicke effect.}, + number = {7}, + urldate = {2015-04-27}, + journal = {Physical Review B}, + author = {Pustovit, Vitaliy N. and Shahbazyan, Tigran V.}, + month = aug, + year = {2010}, + pages = {075429}, + annote = {Velmi zajímavý článek}, + file = {APS Snapshot:/u/46/necadam1/unix/.zotero/zotero/9uf64zmd.default/zotero/storage/GP8ZCDES/PhysRevB.82.html:text/html;PhysRevB.82.075429.pdf:/u/46/necadam1/unix/.zotero/zotero/9uf64zmd.default/zotero/storage/6M953H6A/PhysRevB.82.075429.pdf:application/pdf} +} + +@book{bohren_absorption_1983, + title = {Absorption and scattering of light by small particles}, + url = {http://adsabs.harvard.edu/abs/1983asls.book.....B}, + abstract = {Not Available}, + urldate = {2014-05-09}, + author = {Bohren, Craig F. and Huffman, Donald R.}, + year = {1983}, + keywords = {ABSORPTION, DUST, LIGHT SCATTERING, Particles, THEORY}, + file = {(Wiley science paperback series) Craig F. Bohren, Donald R. Huffman-Absorption and scattering of light by small particles-Wiley-VCH (1998).djvu:/u/46/necadam1/unix/.zotero/zotero/9uf64zmd.default/zotero/storage/HES6WJTP/(Wiley science paperback series) Craig F. Bohren, Donald R. Huffman-Absorption and scattering of light by small particles-Wiley-VCH (1998).djvu:image/vnd.djvu} +} + +@book{taylor_optical_2011, + title = {Optical {Binding} {Phenomena}: {Observations} and {Mechanisms}}, + isbn = {978-3-642-21195-9}, + shorttitle = {Optical {Binding} {Phenomena}}, + abstract = {This thesis addresses optical binding - a new area of interest within the field of optical micromanipulation. It presents, for the first time, a rigorous numerical simulation of some of the key results, along with new experimental findings and also physical interpretations of the results. In an optical trap particles are attracted close to areas of high optical intensities and intensity gradients. So, for example, if two lasers are pointed towards each other (a counter propagating trap) then a single particle is trapped in the centre of the two beams – the system is analogous to a particle being held by two springs in a potential well. If one increases the number of particles in the trap then naively one would expect all the particles to collect in the centre of the well. However, the effect of optical binding means that the presence of one particle affects the distribution of light experienced by another particle, resulting in extremely complex interactions that can lead to unusual 1D and 2D structures to form within the trap. Optical binding is not only of theoretical interest but also has applications in micromanipulation and assembly.}, + language = {en}, + publisher = {Springer Science \& Business Media}, + author = {Taylor, Jonathan M.}, + month = jul, + year = {2011}, + keywords = {Science / Physics / Atomic \& Molecular, Science / Physics / Electricity, Science / Physics / General, Science / Physics / Mathematical \& Computational, Science / Physics / Optics \& Light, Technology \& Engineering / Electrical, Technology \& Engineering / Optics}, + file = {(Springer Theses) Jonathan M. Taylor (auth.)-Optical Binding Phenomena_ Observations and Mechanisms -Springer-Verlag Berlin Heidelberg (2011).pdf:/u/46/necadam1/unix/.zotero/zotero/9uf64zmd.default/zotero/storage/7XKKCD9X/(Springer Theses) Jonathan M. Taylor (auth.)-Optical Binding Phenomena_ Observations and Mechanisms -Springer-Verlag Berlin Heidelberg (2011).pdf:application/pdf} +} + +@article{epton_multipole_1995, + title = {Multipole {Translation} {Theory} for the {Three}-{Dimensional} {Laplace} and {Helmholtz} {Equations}}, + volume = {16}, + issn = {1064-8275}, + url = {http://epubs.siam.org/doi/abs/10.1137/0916051}, + doi = {10.1137/0916051}, + abstract = {The mathematical theory of multipole translation operators for the three-dimensional Laplace and Helmholtz equations is summarized and extended. New results for the Laplace equation include an elementary proof of the inner-to-inner translation theorem, from which follows the definition of a far-field signature function analogous to that of the Helmholtz equation. The theory for the Helmholtz equation is developed in terms of a new convolutional form of the translation operator, which is connected to Rokhlin’s diagonal form by means of Wigner 3-j symbols., The mathematical theory of multipole translation operators for the three-dimensional Laplace and Helmholtz equations is summarized and extended. New results for the Laplace equation include an elementary proof of the inner-to-inner translation theorem, from which follows the definition of a far-field signature function analogous to that of the Helmholtz equation. The theory for the Helmholtz equation is developed in terms of a new convolutional form of the translation operator, which is connected to Rokhlin’s diagonal form by means of Wigner 3-j symbols.}, + number = {4}, + urldate = {2015-08-20}, + journal = {SIAM Journal on Scientific Computing}, + author = {Epton, M. and Dembart, B.}, + month = jul, + year = {1995}, + pages = {865--897}, + file = {epton1995.pdf:/u/46/necadam1/unix/.zotero/zotero/9uf64zmd.default/zotero/storage/7READ52S/epton1995.pdf:application/pdf;Snapshot:/u/46/necadam1/unix/.zotero/zotero/9uf64zmd.default/zotero/storage/XI5S8S66/0916051.html:text/html} +} + +@article{coifman_fast_1993, + title = {The fast multipole method for the wave equation: a pedestrian prescription}, + volume = {35}, + issn = {1045-9243}, + shorttitle = {The fast multipole method for the wave equation}, + doi = {10.1109/74.250128}, + abstract = {A practical and complete, but not rigorous, exposition of the fact multiple method (FMM) is provided. The FMM provides an efficient mechanism for the numerical convolution of the Green's function for the Helmholtz equation with a source distribution and can be used to radically accelerate the iterative solution of boundary-integral equations. In the simple single-stage form presented here, it reduces the computational complexity of the convolution from O(N/sup 2/) to O(N/sup 3/2/), where N is the dimensionality of the problem's discretization.{\textless}{\textgreater}}, + number = {3}, + journal = {IEEE Antennas and Propagation Magazine}, + author = {Coifman, R. and Rokhlin, V. and Wandzura, S.}, + month = jun, + year = {1993}, + keywords = {Acceleration, boundary-integral equations, computational complexity, Convolution, Electromagnetic scattering, electromagnetic wave scattering, fast multiple method, Green's function, Green's function methods, Hardware, Helmholtz equation, iterative solution, Message-oriented middleware, Moment methods, numerical convolution, Partial differential equations, Physics computing, source distribution, Surface waves, wave equation, wave equations}, + pages = {7--12}, + file = {IEEE Xplore Abstract Record:/u/46/necadam1/unix/.zotero/zotero/9uf64zmd.default/zotero/storage/KN5GUPCG/abs_all.html:text/html;IEEE Xplore Full Text PDF:/u/46/necadam1/unix/.zotero/zotero/9uf64zmd.default/zotero/storage/9Z8788XV/Coifman et al. - 1993 - The fast multipole method for the wave equation a.pdf:application/pdf} +} + +@incollection{bostrom_transformation_1991, + title = {Transformation properties of plane, spherical and cylindrical scalar and vector wave functions}, + volume = {1}, + url = {http://lup.lub.lu.se/record/1174356}, + language = {eng}, + urldate = {2014-05-19}, + booktitle = {Acoustic, {Electromagnetic} and {Elastic} {Wave} {Scattering}, {Field} {Representations} and {Introduction} to {Scattering}}, + publisher = {Elsevier Science Publishers}, + author = {Boström, Anders and Kristensson, Gerhard and Ström, Staffan}, + year = {1991}, + keywords = {Technology and Engineering}, + pages = {165--210}, + file = {Chapter[04].pdf:/u/46/necadam1/unix/.zotero/zotero/9uf64zmd.default/zotero/storage/PNMMGQA2/Chapter[04].pdf:application/pdf;Snapshot:/u/46/necadam1/unix/.zotero/zotero/9uf64zmd.default/zotero/storage/9HDC2IWZ/1174356.html:text/html} +} + +@article{xu_electromagnetic_1997, + title = {Electromagnetic scattering by an aggregate of spheres: far field}, + volume = {36}, + issn = {0003-6935, 1539-4522}, + shorttitle = {Electromagnetic scattering by an aggregate of spheres}, + url = {https://www.osapublishing.org/ao/abstract.cfm?uri=ao-36-36-9496}, + doi = {10.1364/AO.36.009496}, + language = {en}, + number = {36}, + urldate = {2015-08-21}, + journal = {Applied Optics}, + author = {Xu, Yu-lin}, + month = dec, + year = {1997}, + pages = {9496}, + file = {ao-36-36-9496.pdf:/u/46/necadam1/unix/.zotero/zotero/9uf64zmd.default/zotero/storage/EIKKM5ZP/ao-36-36-9496.pdf:application/pdf} +} + +@article{xu_electromagnetic_1995, + title = {Electromagnetic scattering by an aggregate of spheres}, + volume = {34}, + issn = {0003-6935, 1539-4522}, + url = {https://www.osapublishing.org/ao/abstract.cfm?uri=ao-34-21-4573}, + doi = {10.1364/AO.34.004573}, + language = {en}, + number = {21}, + urldate = {2015-08-21}, + journal = {Applied Optics}, + author = {Xu, Yu-lin}, + month = jul, + year = {1995}, + pages = {4573}, + file = {ao-34-21-4573.pdf:/u/46/necadam1/unix/.zotero/zotero/9uf64zmd.default/zotero/storage/Z77F8CGC/ao-34-21-4573.pdf:application/pdf} +} + +@article{mackowski_analysis_1991, + title = {Analysis of {Radiative} {Scattering} for {Multiple} {Sphere} {Configurations}}, + volume = {433}, + issn = {1364-5021, 1471-2946}, + url = {http://rspa.royalsocietypublishing.org/content/433/1889/599}, + doi = {10.1098/rspa.1991.0066}, + abstract = {An analysis of radiative scattering for an arbitrary configuration of neighbouring spheres is presented. The analysis builds upon the previously developed superposition solution, in which the scattered field is expressed as a superposition of vector spherical harmonic expansions written about each sphere in the ensemble. The addition theorems for vector spherical harmonics, which transform harmonics from one coordinate system into another, are rederived, and simple recurrence relations for the addition coefficients are developed. The relations allow for a very efficient implementation of the `order of scattering' solution technique for determining the scattered field coefficients for each sphere.}, + language = {en}, + number = {1889}, + urldate = {2015-08-21}, + journal = {Proceedings of the Royal Society of London A: Mathematical, Physical and Engineering Sciences}, + author = {Mackowski, Daniel W.}, + month = jun, + year = {1991}, + pages = {599--614}, + file = {mackowski1991.pdf:/u/46/necadam1/unix/.zotero/zotero/9uf64zmd.default/zotero/storage/254TXAN3/mackowski1991.pdf:application/pdf;Snapshot:/u/46/necadam1/unix/.zotero/zotero/9uf64zmd.default/zotero/storage/QV6MH2N9/599.html:text/html} +} + +@article{cantrell_numerical_2013, + title = {Numerical methods for the accurate calculation of spherical {Bessel} functions and the location of {Mie} resonances}, + url = {http://libtreasures.utdallas.edu/jspui/handle/10735.1/2583}, + urldate = {2015-09-01}, + author = {Cantrell, Cyrus D. and {others}}, + year = {2013}, + file = {ECS-TR-EE-Cantrell-310125.24.pdf:/u/46/necadam1/unix/.zotero/zotero/9uf64zmd.default/zotero/storage/TSJ8T9GS/ECS-TR-EE-Cantrell-310125.24.pdf:application/pdf;Snapshot:/u/46/necadam1/unix/.zotero/zotero/9uf64zmd.default/zotero/storage/7AP37Z92/2583.html:text/html} +} + +@article{gumerov_fast_2001, + title = {Fast, {Exact}, and {Stable} {Computation} of {Multipole} {Translation} and {Rotation} {Coefficients} for the 3-{D} {Helmholtz} {Equation}}, + url = {http://drum.lib.umd.edu/handle/1903/1141}, + abstract = {We develop exact expressions for translations and rotations of local and +multipole fundamental solutions of the Helmholtz equation in spherical +coordinates. These expressions are based on recurrence relations that we +develop, and to our knowledge are presented here for the first time. The +symmetry and other properties of the coefficients are also examined, and +based on these efficient procedures for calculating them are presented. Our +expressions are direct, and do not use the Clebsch-Gordan coefficients or +the Wigner 3-j symbols, though we compare our results with methods that use +these, to prove their accuracy. We test our expressions on a number of +simple calculations, and show their accuracy. For evaluating a \$N\_t\$ term +truncation of the translation (involving \$O(N\_t{\textasciicircum}2)\$ multipoles), compared to +previous exact expressions that rely on the Clebsch-Gordan coefficients or +the Wigner \$3-j\$ symbol that require \$O(N\_t{\textasciicircum}5)\$ operations, our expressions require \$O(N\_t{\textasciicircum}4)\$) evaluations, with a small constant multiplying the order +term. + +The recent trend in evaluating such translations has been to use approximate +"diagonalizations," that require \$O(N\_t{\textasciicircum}3)\$ evaluations with a large +coefficient for the order term. For the Helmholtz equation, these +translations in addition have stabilty problems unless the accuracy of the +truncation and approximate translation are balanced. We derive explicit +exact expressions for achieving "diagonal" translations in \$O(N\_t{\textasciicircum}3)\$ +operations. Our expressions are based on recursive evaluations of multipole +coefficients for rotations, and are accurate and stable, and have a much +smaller coeffiicient for the order term, resulting practically in much fewer +operations. Future use of the developed methods in computational acoustic +scattering, electromagnetic scattering (radar and microwave), optics and +computational biology are expected. + +Cross-referenced as UMIACS-TR-2001-44}, + language = {en\_US}, + urldate = {2015-09-10}, + author = {Gumerov, Nail A. and Duraiswami, Ramani}, + month = sep, + year = {2001}, + file = {Full Text PDF:/u/46/necadam1/unix/.zotero/zotero/9uf64zmd.default/zotero/storage/TDEVFZBV/Gumerov and Duraiswami - 2001 - Fast, Exact, and Stable Computation of Multipole T.pdf:application/pdf;Snapshot:/u/46/necadam1/unix/.zotero/zotero/9uf64zmd.default/zotero/storage/V4SZJT43/1141.html:text/html} +} + +@article{moneda_dyadic_2007, + title = {Dyadic {Green}'s function of a cluster of spheres}, + volume = {24}, + issn = {1084-7529, 1520-8532}, + url = {https://www.osapublishing.org/josaa/abstract.cfm?uri=josaa-24-11-3437}, + doi = {10.1364/JOSAA.24.003437}, + language = {en}, + number = {11}, + urldate = {2015-09-10}, + journal = {Journal of the Optical Society of America A}, + author = {Moneda, Angela P. and Chrissoulidis, Dimitrios P.}, + year = {2007}, + pages = {3437}, + file = {josaa-24-11-3437 (1).pdf:/u/46/necadam1/unix/.zotero/zotero/9uf64zmd.default/zotero/storage/NRM37FIF/josaa-24-11-3437 (1).pdf:application/pdf} +} + +@article{moneda_dyadic_2007-1, + title = {Dyadic {Green}'s function of a sphere with an eccentric spherical inclusion}, + volume = {24}, + issn = {1084-7529, 1520-8532}, + url = {https://www.osapublishing.org/josaa/abstract.cfm?uri=josaa-24-6-1695}, + doi = {10.1364/JOSAA.24.001695}, + language = {en}, + number = {6}, + urldate = {2015-09-10}, + journal = {Journal of the Optical Society of America A}, + author = {Moneda, Angela P. and Chrissoulidis, Dimitrios P.}, + year = {2007}, + pages = {1695}, + file = {josaa-24-6-1695.pdf:/u/46/necadam1/unix/.zotero/zotero/9uf64zmd.default/zotero/storage/33I3IGX6/josaa-24-6-1695.pdf:application/pdf} +} + +@article{pellegrini_interacting_2007, + series = {{EMRS} 2006 {Symposium} {A}: {Current} {Trends} in {Nanoscience} - from {Materials} to {Applications}}, + title = {Interacting metal nanoparticles: {Optical} properties from nanoparticle dimers to core-satellite systems}, + volume = {27}, + issn = {0928-4931}, + shorttitle = {Interacting metal nanoparticles}, + url = {http://www.sciencedirect.com/science/article/pii/S0928493106002657}, + doi = {10.1016/j.msec.2006.07.025}, + abstract = {Prompted by the growing interest in the optical properties of coupled metal nanoclusters, we implemented a code in the framework of Generalized Multiparticle Mie theory (GMM) to simulate far-field properties of strongly interacting spherical particles. In order to validate the code different case studies, including systems modeled for the first time, have been treated. The extinction properties of noble metal nanocluster dimers, chains and core-satellite structures have been computed. Influence of parameters like interparticle distance, incident field polarization, number of multipolar expansions and chain length has been studied. The code provided reliable results in agreement with previous works and proved to be affordable and robust in any of the treated case.}, + number = {5–8}, + urldate = {2015-11-18}, + journal = {Materials Science and Engineering: C}, + author = {Pellegrini, G. and Mattei, G. and Bello, V. and Mazzoldi, P.}, + month = sep, + year = {2007}, + keywords = {Coupled plasmons, Interacting nanoparticles, Optical properties}, + pages = {1347--1350}, + file = {ScienceDirect Full Text PDF:/u/46/necadam1/unix/.zotero/zotero/9uf64zmd.default/zotero/storage/77R8E8NQ/Pellegrini et al. - 2007 - Interacting metal nanoparticles Optical propertie.pdf:application/pdf;ScienceDirect Snapshot:/u/46/necadam1/unix/.zotero/zotero/9uf64zmd.default/zotero/storage/IXCSKSFT/S0928493106002657.html:text/html} +} + +@article{xu_efficient_1998, + title = {Efficient {Evaluation} of {Vector} {Translation} {Coefficients} in {Multiparticle} {Light}-{Scattering} {Theories}}, + volume = {139}, + issn = {0021-9991}, + url = {http://www.sciencedirect.com/science/article/pii/S0021999197958678}, + doi = {10.1006/jcph.1997.5867}, + abstract = {Vector addition theorems are a necessary ingredient in the analytical solution of electromagnetic multiparticle-scattering problems. These theorems include a large number of vector addition coefficients. There exist three basic types of analytical expressions for vector translation coefficients: Stein's (Quart. Appl. Math.19, 15 (1961)), Cruzan's (Quart. Appl. Math.20, 33 (1962)), and Xu's (J. Comput. Phys.127, 285 (1996)). Stein's formulation relates vector translation coefficients with scalar translation coefficients. Cruzan's formulas use the Wigner 3jm symbol. Xu's expressions are based on the Gaunt coefficient. Since the scalar translation coefficient can also be expressed in terms of the Gaunt coefficient, the key to the expeditious and reliable calculation of vector translation coefficients is the fast and accurate evaluation of the Wigner 3jm symbol or the Gaunt coefficient. We present highly efficient recursive approaches to accurately evaluating Wigner 3jm symbols and Gaunt coefficients. Armed with these recursive approaches, we discuss several schemes for the calculation of the vector translation coefficients, using the three general types of formulation, respectively. Our systematic test calculations show that the three types of formulas produce generally the same numerical results except that the algorithm of Stein's type is less accurate in some particular cases. These extensive test calculations also show that the scheme using the formulation based on the Gaunt coefficient is the most efficient in practical computations.}, + number = {1}, + urldate = {2015-11-18}, + journal = {Journal of Computational Physics}, + author = {Xu, Yu-lin}, + month = jan, + year = {1998}, + pages = {137--165}, + annote = {N.B. erratum regarding eqs (50,52,53) +http://www.sciencedirect.com/science/article/pii/S0021999197956874}, + file = {ScienceDirect Full Text PDF:/u/46/necadam1/unix/.zotero/zotero/9uf64zmd.default/zotero/storage/STV5263F/Xu - 1998 - Efficient Evaluation of Vector Translation Coeffic.pdf:application/pdf;ScienceDirect Snapshot:/u/46/necadam1/unix/.zotero/zotero/9uf64zmd.default/zotero/storage/VMGZRSAA/S0021999197958678.html:text/html} +} + +@article{xu_calculation_1996, + title = {Calculation of the {Addition} {Coefficients} in {Electromagnetic} {Multisphere}-{Scattering} {Theory}}, + volume = {127}, + issn = {0021-9991}, + url = {http://www.sciencedirect.com/science/article/pii/S0021999196901758}, + doi = {10.1006/jcph.1996.0175}, + abstract = {One of the most intractable problems in electromagnetic multisphere-scattering theory is the formulation and evaluation of vector addition coefficients introduced by the addition theorems for vector spherical harmonics. This paper presents an efficient approach for the calculation of both scalar and vector translational addition coefficients, which is based on fast evaluation of the Gaunt coefficients. The paper also rederives the analytical expressions for the vector translational addition coefficients and discusses the strengths and limitations of other formulations and numerical techniques found in the literature. Numerical results from the formulation derived in this paper agree with those of a previously published recursion scheme that completely avoids the use of the Gaunt coefficients, but the method of direct calculation proposed here reduces the computing time by a factor of 4–6.}, + number = {2}, + urldate = {2015-11-24}, + journal = {Journal of Computational Physics}, + author = {Xu, Yu-lin}, + month = sep, + year = {1996}, + pages = {285--298}, + file = {ScienceDirect Snapshot:/u/46/necadam1/unix/.zotero/zotero/9uf64zmd.default/zotero/storage/H98A3TTE/S0021999196901758.html:text/html;xu1996.pdf:/u/46/necadam1/unix/.zotero/zotero/9uf64zmd.default/zotero/storage/CBABI5M4/xu1996.pdf:application/pdf} +} + +@misc{kristensson_spherical_2014, + title = {Spherical {Vector} {Waves}}, + url = {http://www.eit.lth.se/fileadmin/eit/courses/eit080f/Literature/book.pdf}, + urldate = {2014-05-20}, + author = {Kristensson, Gerhard}, + month = jan, + year = {2014} +} + +@book{jackson_classical_1998, + address = {New York}, + edition = {3 edition}, + title = {Classical {Electrodynamics} {Third} {Edition}}, + isbn = {978-0-471-30932-1}, + abstract = {A revision of the defining book covering the physics and classical mathematics necessary to understand electromagnetic fields in materials and at surfaces and interfaces. The third edition has been revised to address the changes in emphasis and applications that have occurred in the past twenty years.}, + language = {English}, + publisher = {Wiley}, + author = {Jackson, John David}, + month = aug, + year = {1998}, + file = {John David Jackson-Classical Electrodynamics-Wiley (1999).djvu:/u/46/necadam1/unix/.zotero/zotero/9uf64zmd.default/zotero/storage/3BWPD4BK/John David Jackson-Classical Electrodynamics-Wiley (1999).djvu:image/vnd.djvu} +} + +@article{mandelstam_quantum_1962, + title = {Quantum electrodynamics without potentials}, + volume = {19}, + issn = {0003-4916}, + url = {http://www.sciencedirect.com/science/article/pii/0003491662902324}, + doi = {10.1016/0003-4916(62)90232-4}, + abstract = {A scheme is proposed for quantizing electrodynamics in terms of the electromagnetic fields without the introduction of potentials. The equations are relativistically covariant and do not require the introduction of unphysical states and an indefinite metric. Calculations carried out according to current quantization methods in the Coulomb or Lorentz gauges are justified in the new formalism. The theory exhibits an analogy between phases of operators and electromagnetic fields on the one hand, and coordinate systems and space curvature on the other. It is suggested that this analogy may be useful in quantizing the gravitational field.}, + number = {1}, + urldate = {2014-11-23}, + journal = {Annals of Physics}, + author = {Mandelstam, Stanley}, + month = jul, + year = {1962}, + pages = {1--24} +} + +@article{philbin_canonical_2010, + title = {Canonical quantization of macroscopic electromagnetism}, + volume = {12}, + issn = {1367-2630}, + url = {http://iopscience.iop.org/1367-2630/12/12/123008}, + doi = {10.1088/1367-2630/12/12/123008}, + abstract = {Application of the standard canonical quantization rules of quantum field theory to macroscopic electromagnetism has encountered obstacles due to material dispersion and absorption. This has led to a phenomenological approach to macroscopic quantum electrodynamics where no canonical formulation is attempted. In this paper macroscopic electromagnetism is canonically quantized. The results apply to any linear, inhomogeneous, magnetodielectric medium with dielectric functions that obey the Kramers–Kronig relations. The prescriptions of the phenomenological approach are derived from the canonical theory.}, + language = {en}, + number = {12}, + urldate = {2014-03-28}, + journal = {New Journal of Physics}, + author = {Philbin, T. G.}, + month = dec, + year = {2010}, + keywords = {\_tablet}, + pages = {123008}, + file = {Philbin_2010_Canonical quantization of macroscopic electromagnetism.pdf:/u/46/necadam1/unix/.zotero/zotero/9uf64zmd.default/zotero/storage/MHJA3DSX/Philbin_2010_Canonical quantization of macroscopic electromagnetism.pdf:application/pdf;Snapshot:/u/46/necadam1/unix/.zotero/zotero/9uf64zmd.default/zotero/storage/R5AHRZRR/123008.html:text/html} +} + +@article{huttner_quantization_1992, + title = {Quantization of the electromagnetic field in dielectrics}, + volume = {46}, + url = {http://link.aps.org/doi/10.1103/PhysRevA.46.4306}, + doi = {10.1103/PhysRevA.46.4306}, + abstract = {We present a fully canonical quantization scheme for the electromagnetic field in dispersive and lossy linear dielectrics. This scheme is based on a microscopic model, in which the medium is represented by a collection of interacting matter fields. We calculate the exact eigenoperators for the coupled system and express the electromagnetic field operators in terms of them. The dielectric constant of the medium is explicitly derived and is shown to satisfy the Kramers-Kronig relations. We apply these results to treat the propagation of light in dielectrics and obtain simple expressions for the electromagnetic field in the medium in terms of space-dependent creation and annihilation operators. These operators satisfy a set of equal-space commutation relations and obey spatial Langevin equations of evolution. This justifies the use of such operators in phenomenological models in quantum optics. We also obtain two interesting relationships between the group and the phase velocity in dielectrics.}, + number = {7}, + urldate = {2014-03-28}, + journal = {Physical Review A}, + author = {Huttner, Bruno and Barnett, Stephen M.}, + month = oct, + year = {1992}, + keywords = {\_tablet}, + pages = {4306--4322}, + file = {APS Snapshot:/u/46/necadam1/unix/.zotero/zotero/9uf64zmd.default/zotero/storage/VP7HX7MC/PhysRevA.46.html:text/html;Huttner_Barnett_1992_Quantization of the electromagnetic field in dielectrics.pdf:/u/46/necadam1/unix/.zotero/zotero/9uf64zmd.default/zotero/storage/EQ6HCUDJ/Huttner_Barnett_1992_Quantization of the electromagnetic field in dielectrics.pdf:application/pdf} +} + +@article{huttner_canonical_1991, + title = {Canonical {Quantization} of {Light} in a {Linear} {Dielectric}}, + volume = {16}, + issn = {0295-5075}, + url = {http://iopscience.iop.org/0295-5075/16/2/010}, + doi = {10.1209/0295-5075/16/2/010}, + abstract = {Quantization of the macroscopic electromagnetic field via effective susceptibilities leads to inconsistencies if the medium is dispersive. A canonical quantization scheme has to take explicit account of the matter field. By introducing a simple model for the matter, we are able to resolve some of the difficulties highlighted in the recent literature. We demonstrate the fundamental significance of the electromagnetic energy flux (rather than the density) and justify the use of temporal modes of the field. Our analysis leads to an apparently unknown relationship between the group and phase velocity in a linear dielectric medium.}, + language = {en}, + number = {2}, + urldate = {2015-05-08}, + journal = {EPL (Europhysics Letters)}, + author = {Huttner, B. and Baumberg, J. J. and Barnett, S. M.}, + month = sep, + year = {1991}, + pages = {177}, + file = {Full Text PDF:/u/46/necadam1/unix/.zotero/zotero/9uf64zmd.default/zotero/storage/5FJP7FZI/Huttner et al. - 1991 - Canonical Quantization of Light in a Linear Dielec.pdf:application/pdf;Snapshot:/u/46/necadam1/unix/.zotero/zotero/9uf64zmd.default/zotero/storage/CCZPG7NS/010.html:text/html} +} + +@article{hopfield_theory_1958, + title = {Theory of the {Contribution} of {Excitons} to the {Complex} {Dielectric} {Constant} of {Crystals}}, + volume = {112}, + url = {http://link.aps.org/doi/10.1103/PhysRev.112.1555}, + doi = {10.1103/PhysRev.112.1555}, + abstract = {It is shown that the ordinary semiclassical theory of the absorption of light by exciton states is not completely satisfactory (in contrast to the case of absorption due to interband transitions). A more complete theory is developed. It is shown that excitons are approximate bosons, and, in interaction with the electromagnetic field, the exciton field plays the role of the classical polarization field. The eigenstates of the system of crystal and radiation field are mixtures of photons and excitons. The ordinary one-quantum optical lifetime of an excitation is infinite. Absorption occurs only when "three-body" processes are introduced. The theory includes "local field" effects, leading to the Lorentz local field correction when it is applicable. A Smakula equation for the oscillator strength in terms of the integrated absorption constant is derived.}, + number = {5}, + urldate = {2015-05-08}, + journal = {Physical Review}, + author = {Hopfield, J. J.}, + month = dec, + year = {1958}, + pages = {1555--1567}, + file = {APS Snapshot:/u/46/necadam1/unix/.zotero/zotero/9uf64zmd.default/zotero/storage/PUNM4SES/PhysRev.112.html:text/html;PhysRev.112.1555 (1).pdf:/u/46/necadam1/unix/.zotero/zotero/9uf64zmd.default/zotero/storage/7G8J8PTA/PhysRev.112.1555 (1).pdf:application/pdf} +} + +@book{cohen-tannoudji_atom-photon_1998, + address = {New York}, + title = {Atom-{Photon} {Interactions}: {Basic} {Processes} and {Applications}}, + isbn = {978-0-471-29336-1}, + shorttitle = {Atom-{Photon} {Interactions}}, + abstract = {Atom-Photon Interactions: Basic Processes and Applications allows the reader to master various aspects of the physics of the interaction between light and matter. It is devoted to the study of the interactions between photons and atoms in atomic and molecular physics, quantum optics, and laser physics. The elementary processes in which photons are emitted, absorbed, scattered, or exchanged between atoms are treated in detail and described using diagrammatic representation. The book presents different theoretical approaches, including: * Perturbative methods * The resolvent method * Use of the master equation * The Langevin equation * The optical Bloch equations * The dressed-atom approach Each method is presented in a self-contained manner so that it may be studied independently. Many applications of these approaches to simple and important physical phenomena are given to illustrate the potential and limitations of each method.}, + language = {English}, + publisher = {Wiley-VCH}, + author = {Cohen-Tannoudji, Claude and Dupont-Roc, Jacques and Grynberg, Gilbert}, + month = mar, + year = {1998} +} + +@book{cohen-tannoudji_photons_1997, + address = {Weinheim}, + edition = {1st THUS edition}, + title = {Photons and {Atoms}: {Introduction} to {Quantum} {Electrodynamics}}, + isbn = {978-0-471-18433-1}, + shorttitle = {Photons and {Atoms}}, + abstract = {Photons and Atoms Photons and Atoms: Introduction to Quantum Electrodynamics provides the necessary background to understand the various physical processes associated with photon-atom interactions. It starts with elementary quantum theory and classical electrodynamics and progresses to more advanced approaches. A critical comparison is made between these different, although equivalent, formulations of quantum electrodynamics. Using this format, the reader is offered a gradual, yet flexible introduction to quantum electrodynamics, avoiding formal discussions and excessive shortcuts. Complementing each chapter are numerous examples and exercises that can be used independently from the rest of the book to extend each chapter in many disciplines depending on the interests and needs of the reader.}, + language = {English}, + publisher = {Wiley-VCH}, + author = {Cohen-Tannoudji, Claude and Dupont-Roc, Jacques and Grynberg, Gilbert}, + month = mar, + year = {1997}, + file = {[Claude_Cohen-Tannoudji\;_Jacques_Dupont-Roc\;_Gilbe(BookZZ.org).djvu:/u/46/necadam1/unix/.zotero/zotero/9uf64zmd.default/zotero/storage/C8UTI3EA/[Claude_Cohen-Tannoudji\;_Jacques_Dupont-Roc\;_Gilbe(BookZZ.org).djvu:image/vnd.djvu} +} + +@article{gruner_green-function_1996, + title = {Green-function approach to the radiation-field quantization for homogeneous and inhomogeneous {Kramers}-{Kronig} dielectrics}, + volume = {53}, + url = {http://link.aps.org/doi/10.1103/PhysRevA.53.1818}, + doi = {10.1103/PhysRevA.53.1818}, + abstract = {A quantization scheme for the radiation field in dispersive and absorptive linear dielectrics is developed, which applies to both bulk material and multilayer dielectric structures. Starting from the phenomenological Maxwell equations, where the properties of the dielectric are described by a permittivity consistent with the Kramers-Kronig relations, an expansion of the field operators is performed that is based on the Green function of the classical Maxwell equations and preserves the equal-time canonical field commutation relations. In particular, in frequency intervals with approximately vanishing absorption the concept of quantization through mode expansion for dispersive dielectrics is recognized. The theory further reveals that weak absorption gives rise to space-dependent mode operators that spatially evolve according to quantum Langevin equations in the space domain. To illustrate the applicability of the theory to inhomogeneous structures, the quantization of the radiation field in a dispersive and absorptive one-interface dielectric is performed. © 1996 The American Physical Society.}, + number = {3}, + urldate = {2016-05-13}, + journal = {Physical Review A}, + author = {Gruner, T. and Welsch, D.-G.}, + month = mar, + year = {1996}, + pages = {1818--1829}, + file = {APS Snapshot:/u/46/necadam1/unix/.zotero/zotero/9uf64zmd.default/zotero/storage/KNE8THEZ/PhysRevA.53.html:text/html;PhysRevA.53.1818.pdf:/u/46/necadam1/unix/.zotero/zotero/9uf64zmd.default/zotero/storage/7PDC6W3U/PhysRevA.53.1818.pdf:application/pdf} +} \ No newline at end of file diff --git a/README.rst b/README.rst index a167c02..cee538a 100644 --- a/README.rst +++ b/README.rst @@ -2,3 +2,15 @@ Quantum photonic multiple scattering ==================================== TODO description + +Installation +============ +The package depends on numpy, scipy, cython and customized version of py_gmm. +The first three can be obtained by pip, the last one can be obtained from github: + +git clone --branch standalone_mie https://github.com/texnokrates/py_gmm.git + +After all dependencies are installed, install qpms to your local python library using + +python3 setup.py install --user + diff --git a/Scattering and Shit.lyx b/Scattering and Shit.lyx index 2e39e81..1cdefc7 100644 --- a/Scattering and Shit.lyx +++ b/Scattering and Shit.lyx @@ -79,8 +79,21 @@ \begin_body +\begin_layout Standard +\begin_inset FormulaMacro +\newcommand{\vect}[1]{\mathbf{#1}} +\end_inset + + +\begin_inset FormulaMacro +\newcommand{\ud}{\mathrm{d}} +\end_inset + + +\end_layout + \begin_layout Title -Electromagnetic multiple scattering, spherical waves and shit +Electromagnetic multiple scattering, spherical waves and **** \end_layout \begin_layout Author @@ -95,33 +108,781 @@ Zillion conventions for spherical vector waves Legendre polynomials and spherical harmonics: messy from the very beginning \end_layout +\begin_layout Standard +\begin_inset Marginal +status open + +\begin_layout Plain Layout +FIXME check the Condon-Shortley phases. +\end_layout + +\end_inset + + +\end_layout + +\begin_layout Standard +Associated Legendre polynomial of degree +\begin_inset Formula $l\ge0$ +\end_inset + + and order +\begin_inset Formula $m,$ +\end_inset + + +\begin_inset Formula $l\ge m\ge-l$ +\end_inset + +, is given by the recursive relation +\begin_inset Formula +\[ +P_{l}^{-m}=\underbrace{\left(-1\right)^{m}}_{\mbox{Condon-Shortley phase}}\frac{1}{2^{l}l!}\left(1-x^{2}\right)^{m/2}\frac{\ud^{l+m}}{\ud x^{l+m}}\left(x^{2}-1\right)^{l}. +\] + +\end_inset + +There is a relation between the positive and negative orders, +\end_layout + +\begin_layout Standard +\begin_inset Formula +\[ +P_{l}^{-m}=\underbrace{\left(-1\right)^{m}}_{\mbox{C.-S. p.}}\frac{\left(l-m\right)!}{\left(l+m\right)!}P_{l}^{m}\left(\cos\theta\right),\quad m\ge0. +\] + +\end_inset + +The index +\begin_inset Formula $l$ +\end_inset + + (in certain notations, it is often +\begin_inset Formula $n$ +\end_inset + +) is called +\emph on +degree +\emph default +, index +\begin_inset Formula $m$ +\end_inset + + is the +\emph on +order +\emph default +. + These two terms are then transitively used for all the object which build + on the associated Legendre polynomials, i.e. + spherical harmonics, vector spherical harmonics, spherical waves etc. +\end_layout + +\begin_layout Subsection +Kristensson +\end_layout + +\begin_layout Standard +Kristensson uses the Condon-Shortley phase, so (sect. + [K]D.2) +\end_layout + +\begin_layout Standard +\begin_inset Formula +\[ +Y_{lm}\left(\hat{\vect r}\right)=\left(-1\right)^{m}\sqrt{\frac{2l+1}{4\pi}\frac{\left(l-m\right)!}{\left(l+m\right)!}}P_{l}^{m}\left(\cos\theta\right)e^{im\phi} +\] + +\end_inset + + +\begin_inset Formula +\[ +Y_{lm}^{\dagger}\left(\hat{\vect r}\right)=Y_{lm}^{*}\left(\hat{\vect r}\right) +\] + +\end_inset + + +\begin_inset Formula +\[ +Y_{l,-m}\left(\hat{\vect r}\right)=\left(-1\right)^{m}Y_{lm}^{\dagger}\left(\hat{\vect r}\right) +\] + +\end_inset + + +\end_layout + +\begin_layout Standard +Orthonormality: +\begin_inset Formula +\[ +\int Y_{lm}\left(\hat{\vect r}\right)Y_{l'm'}^{\dagger}\left(\hat{\vect r}\right)\,\ud\Omega=\delta_{ll'}\delta_{mm'} +\] + +\end_inset + + +\end_layout + \begin_layout Section -Pi and tau (?), spherical Bessel functions +Pi and tau +\end_layout + +\begin_layout Subsection +Taylor +\end_layout + +\begin_layout Standard +\begin_inset Formula +\begin{eqnarray*} +\tilde{\pi}_{mn}\left(\cos\theta\right) & = & \sqrt{\frac{2n+1}{4\pi}\frac{\left(n-m\right)!}{\left(n+m\right)!}}\frac{m}{\sin\theta}P_{n}^{m}\left(\cos\theta\right)\\ +\tilde{\tau}_{mn}\left(\cos\theta\right) & = & \sqrt{\frac{2n+1}{4\pi}\frac{\left(n-m\right)!}{\left(n+m\right)!}}\frac{\ud}{\ud\theta}P_{n}^{m}\left(\cos\theta\right) +\end{eqnarray*} + +\end_inset + + \end_layout \begin_layout Section Vector spherical harmonics (?) \end_layout -\begin_layout Section -All the conventions +\begin_layout Subsection +Kristensson \end_layout \begin_layout Standard +Original formulation, sect. + [K]D.3.3 +\end_layout + +\begin_layout Standard +\begin_inset Formula +\begin{eqnarray*} +\vect A_{1lm}\left(\hat{\vect r}\right) & = & \frac{1}{\sqrt{l\left(l+1\right)}}\left(\hat{\vect{\theta}}\frac{1}{\sin\theta}\frac{\partial}{\partial\phi}Y_{lm}\left(\hat{\vect r}\right)-\hat{\vect{\phi}}\frac{\partial}{\partial\theta}Y_{lm}\left(\hat{\vect r}\right)\right)\\ +\vect A_{2lm}\left(\hat{\vect r}\right) & = & \frac{1}{\sqrt{l\left(l+1\right)}}\left(\hat{\vect{\theta}}\frac{\partial}{\partial\phi}Y_{lm}\left(\hat{\vect r}\right)-\hat{\vect{\phi}}\frac{1}{\sin\theta}\frac{\partial}{\partial\theta}Y_{lm}\left(\hat{\vect r}\right)\right)\\ +\vect A_{3lm}\left(\hat{\vect r}\right) & = & \hat{\vect r}Y_{lm}\left(\hat{\vect r}\right) +\end{eqnarray*} + +\end_inset + +Normalisation: +\begin_inset Formula +\[ +\int\vect A_{n}\left(\hat{\vect r}\right)\cdot\vect A_{n'}^{\dagger}\left(\hat{\vect r}\right)\,\ud\Omega=\delta_{nn'} +\] + +\end_inset + +Here +\begin_inset Formula $\mbox{ }^{\dagger}$ +\end_inset + + means just complex conjugate, apparently (see footnote on p. + 89). +\end_layout + +\begin_layout Subsection +Jackson +\end_layout + +\begin_layout Standard +\begin_inset CommandInset citation +LatexCommand cite +after "(9.101)" +key "jackson_classical_1998" + +\end_inset + +: +\begin_inset Formula +\[ +\vect X_{lm}(\theta,\phi)=\frac{1}{\sqrt{l(l+1)}}\vect LY_{lm}(\theta,\phi) +\] + +\end_inset + +where +\begin_inset CommandInset citation +LatexCommand cite +after "(9.119)" +key "jackson_classical_1998" + +\end_inset + + +\begin_inset Formula +\[ +\vect L=\frac{1}{i}\left(\vect r\times\vect{\nabla}\right) +\] + +\end_inset + +for its expression in spherical coordinates and other properties check Jackson's + book around the definitions. +\end_layout + +\begin_layout Standard +Normalisation +\begin_inset CommandInset citation +LatexCommand cite +after "(9.120)" +key "jackson_classical_1998" + +\end_inset + +: +\begin_inset Formula +\[ +\int\vect X_{l'm'}^{*}\cdot\vect X_{lm}\,\ud\Omega=\delta_{ll'}\delta_{mm'} +\] + +\end_inset + + +\end_layout + +\begin_layout Standard +Local sum rule +\begin_inset CommandInset citation +LatexCommand cite +after "(9.153)" +key "jackson_classical_1998" + +\end_inset + +: +\begin_inset Formula +\[ +\sum_{m=-l}^{l}\left|\vect X_{lm}(\theta,\phi)^{2}\right|=\frac{2l+1}{4\pi} +\] + +\end_inset + + +\end_layout + +\begin_layout Section +Spherical Bessel functions +\begin_inset CommandInset label +LatexCommand label +name "sec:Spherical-Bessel-functions" + +\end_inset + + +\end_layout + +\begin_layout Standard +Cf. + [DLMF] §10.47–60. +\end_layout + +\begin_layout Standard +The radial dependence of spherical vector waves is given by the spherical + Bessel functions and their first derivatives. + Commonly, the following notation is adopted +\begin_inset Formula +\begin{eqnarray*} +z_{n}^{(1)}(x) & = & j_{n}(x),\\ +z_{n}^{(2)}(x) & = & y_{n}(x),\\ +z_{n}^{(3)}(x) & = & h_{n}^{(1)}(x)=j_{n}(x)+iy_{n}(x),\\ +z_{n}^{(4)}(x) & = & h_{n}^{(2)}(x)=j_{n}(x)-iy_{n}(x). +\end{eqnarray*} + +\end_inset + +Here, +\begin_inset Formula $j_{n}$ +\end_inset + + is the spherical Bessel function of first kind (regular), +\begin_inset Formula $y_{j}$ +\end_inset + + is the spherical Bessel function of second kind (singular), and +\begin_inset Formula $h_{n}^{(1)},h_{n}^{(2)}$ +\end_inset + + are the Hankel functions a.k.a. + spherical Bessel functions of third kind. + In spherical vector waves, +\begin_inset Formula $j_{n}$ +\end_inset + + corresponds to regular waves, +\begin_inset Formula $h^{(1)}$ +\end_inset + + corresponds (by the usual convention) to outgoing waves, and +\begin_inset Formula $h^{(2)}$ +\end_inset + + corresponds to incoming waves. + To describe scattering, we need two sets of waves with two different types + of spherical Bessel functions +\begin_inset Formula $z_{n}^{(J)}$ +\end_inset + +. + Most common choice is +\begin_inset Formula $J=1,3$ +\end_inset + +, because if we decompose the field into spherical waves centered at +\begin_inset Formula $\vect r_{0}$ +\end_inset + +, the field produced by other sources (e.g. + spherical waves from other scatterers or a plane wave) is always regular + at +\begin_inset Formula $\vect r_{0}$ +\end_inset + +. + Second choice which makes a bit of sense is +\begin_inset Formula $J=3,4$ +\end_inset + + as it leads to a nice expression for the energy transport. +\end_layout + +\begin_layout Subsection +Limiting Forms +\end_layout + +\begin_layout Standard +[DLMF] §10.52: +\end_layout + +\begin_layout Subsection +\begin_inset Formula $z\to0$ +\end_inset + + +\end_layout + +\begin_layout Standard +\begin_inset Formula +\begin{eqnarray*} +j_{n}(z) & \sim & z^{n}/(2n+1)!!\\ +h_{n}^{(1)}(z)\sim iy(z) & \sim & -i\left(2n+1\right)!!/z^{n+1} +\end{eqnarray*} + +\end_inset + + +\end_layout + +\begin_layout Section +Spherical vector waves +\end_layout + +\begin_layout Standard +TODO \begin_inset Formula $M,N,\psi,\chi,\widetilde{M},\widetilde{N},u,v,w,\dots$ \end_inset , sine/cosine convention (B&H), ... \end_layout +\begin_layout Standard +There are two mutually orthogonal types of divergence-free (everywhere except + in the origin for singular waves) spherical vector waves, which I call + electric and magnetic, given by the type of multipole source to which they + correspond. + This is another distinction than the regular/singular/ingoing/outgoing + waves given by the type of the radial dependence (cf. + section +\begin_inset CommandInset ref +LatexCommand ref +reference "sec:Spherical-Bessel-functions" + +\end_inset + +). + Oscillating electric current in a tiny rod parallel to its axis will generate + electric dipole waves (net dipole moment of magnetic current is zero) moment + , whereas oscillating electric current in a tiny circular loop will generate + magnetic dipole waves (net dipole moment of electric current is zero). +\end_layout + +\begin_layout Standard +In the usual cases we encounter, the part described by the magnetic waves + is pretty small. +\end_layout + +\begin_layout Subsection +Taylor +\end_layout + +\begin_layout Standard +Definition [T](2.40); +\begin_inset Formula $\widetilde{\vect N}_{mn}^{(j)},\widetilde{\vect M}_{mn}^{(j)}$ +\end_inset + + are the electric and magnetic waves, respectively: +\end_layout + +\begin_layout Standard +\begin_inset Formula +\begin{eqnarray*} +\widetilde{\vect N}_{mn}^{(j)} & = & \frac{n(n+1)}{kr}\sqrt{\frac{2n+1}{4\pi}\frac{\left(n-m\right)!}{\left(n+m\right)!}}P_{n}^{m}\left(\cos\theta\right)e^{im\phi}z_{n}^{j}\left(kr\right)\hat{\vect r}\\ + & & +\left[\tilde{\tau}_{mn}\left(\cos\theta\right)\hat{\vect{\theta}}+i\tilde{\pi}_{mn}\left(\cos\theta\right)\hat{\vect{\phi}}\right]e^{im\phi}\frac{1}{kr}\frac{\ud\left(kr\, z_{n}^{j}\left(kr\right)\right)}{\ud(kr)}\\ +\widetilde{\vect M}_{mn}^{(j)} & = & \left[i\tilde{\pi}_{mn}\left(\cos\theta\right)\hat{\vect{\theta}}-\tilde{\tau}_{mn}\left(\cos\theta\right)\hat{\vect{\phi}}\right]e^{im\phi}z_{n}^{j}\left(kr\right) +\end{eqnarray*} + +\end_inset + + +\end_layout + +\begin_layout Subsection +Kristensson +\end_layout + +\begin_layout Standard +Definition [K](2.4.6); +\begin_inset Formula $\vect u_{\tau lm},\vect v_{\tau lm},\vect w_{\tau lm}$ +\end_inset + + are the waves with +\begin_inset Formula $j=3,1,4$ +\end_inset + + respectively, i.e. + outgoing, regular and incoming waves. + The first index distinguishes between the electric ( +\begin_inset Formula $\tau=2$ +\end_inset + +) and magnetic ( +\begin_inset Formula $\tau=1$ +\end_inset + +). + Kristensson uses a multiindex +\begin_inset Formula $n\equiv(\tau,l,m)$ +\end_inset + + to simlify the notation. +\begin_inset Formula +\begin{eqnarray*} +\left(\vect{u/v/w}\right)_{2lm} & = & \frac{1}{kr}\frac{\ud\left(kr\, z_{l}^{(j)}\left(kr\right)\right)}{\ud\, kr}\vect A_{2lm}\left(\hat{\vect r}\right)+\sqrt{l\left(l+1\right)}\frac{z_{l}^{(j)}(kr)}{kr}\vect A_{3lm}\left(\hat{\vect r}\right)\\ +\left(\vect{u/v/w}\right)_{1lm} & = & z_{l}^{(j)}\left(kr\right)\vect A_{1lm}\left(\hat{\vect r}\right) +\end{eqnarray*} + +\end_inset + + +\end_layout + +\begin_layout Subsection +Relation between Kristensson and Taylor +\begin_inset CommandInset label +LatexCommand label +name "sub:Kristensson-v-Taylor" + +\end_inset + + +\end_layout + +\begin_layout Standard +Kristensson's and Taylor's VSWFs seem to differ only by an +\begin_inset Formula $l$ +\end_inset + +-dependent normalization factor, and notation of course (n.b. + the inverse index order) +\begin_inset Formula +\begin{eqnarray*} +\left(\vect{u/v/w}\right)_{2lm} & = & \frac{\widetilde{\vect N}_{ml}^{(3/1/4)}}{\sqrt{l\left(l+1\right)}}\\ +\left(\vect{u/v/w}\right)_{1lm} & = & \frac{\widetilde{\vect M}_{ml}^{(3/1/4)}}{\sqrt{l\left(l+1\right)}} +\end{eqnarray*} + +\end_inset + + +\end_layout + \begin_layout Section Plane wave expansion \end_layout +\begin_layout Subsection +Taylor +\end_layout + +\begin_layout Standard +\begin_inset Formula $x$ +\end_inset + +-polarised, +\begin_inset Formula $z$ +\end_inset + +-propagating plane wave, +\begin_inset Formula $\vect E=E_{0}\hat{\vect x}e^{i\vect k\cdot\hat{\vect z}}$ +\end_inset + + (CHECK): +\begin_inset Formula +\begin{eqnarray*} +\vect E & = & -i\left(p_{mn}\widetilde{\vect N}_{mn}^{(1)}+q_{mn}\widetilde{\vect M}_{mn}^{(1)}\right)\\ +p_{mn} & = & E_{0}\frac{4\pi i^{n}}{n(n+1)}\tilde{\tau}_{mn}(1)\\ +q_{mn} & = & E_{0}\frac{4\pi i^{n}}{n(n+1)}\tilde{\pi}_{mn}(1) +\end{eqnarray*} + +\end_inset + +while it can be shown that +\begin_inset Formula +\begin{eqnarray*} +\tilde{\pi}_{mn}(1) & = & -\frac{1}{2}\sqrt{\frac{\left(2n+1\right)\left(n\left(n+1\right)\right)}{4\pi}}\left(\delta_{m,1}+\delta_{m,-1}\right)\\ +\tilde{\tau}_{mn}(1) & = & -\frac{1}{2}\sqrt{\frac{\left(2n+1\right)\left(n\left(n+1\right)\right)}{4\pi}}\left(\delta_{m,1}-\delta_{m,-1}\right) +\end{eqnarray*} + +\end_inset + + +\end_layout + +\begin_layout Subsection +Kristensson +\end_layout + +\begin_layout Standard +\begin_inset Formula $x$ +\end_inset + +-polarised, +\begin_inset Formula $z$ +\end_inset + +-propagating plane wave, +\begin_inset Formula $\vect E=E_{0}\hat{\vect x}e^{i\vect k\cdot\hat{\vect z}}$ +\end_inset + + (CHECK, ): +\begin_inset Formula +\[ +\vect E=\sum_{n}a_{n}\vect v_{n} +\] + +\end_inset + + +\begin_inset Formula +\begin{eqnarray*} +a_{1lm} & = & E_{0}i^{l+1}\sqrt{\left(2l+1\right)\pi}\left(\delta_{m,1}+\delta_{m,-1}\right)\\ +a_{2lm} & = & E_{0}i^{l+1}\sqrt{\left(2l+1\right)\pi}\left(\delta_{m,1}+\delta_{m,-1}\right) +\end{eqnarray*} + +\end_inset + + +\end_layout + \begin_layout Section Radiated energy \end_layout +\begin_layout Standard +In this section I summarize the formulae for power +\begin_inset Formula $P$ +\end_inset + + radiated from the system. + For an absorbing scatterer, this should be negative (n.b. + sign conventions can be sometimes confusing). + If the system is excited by a plane wave with intensity +\begin_inset Formula $E_{0}$ +\end_inset + +, this can be used to calculate the absorption cross section (TODO check + if it should be multiplied by the 2), +\begin_inset Formula +\[ +\sigma_{\mathrm{abs}}=-\frac{2P}{\varepsilon\varepsilon_{0}\left|E_{0}\right|^{2}}. +\] + +\end_inset + + +\end_layout + +\begin_layout Subsection +Kristensson +\begin_inset CommandInset label +LatexCommand label +name "sub:Radiated enenergy-Kristensson" + +\end_inset + + +\end_layout + +\begin_layout Standard +Sect. + [K]2.6.2; here this form of expansion is assumed: +\begin_inset Formula +\begin{equation} +\vect E\left(\vect r,\omega\right)=k\sqrt{\eta_{0}\eta}\sum_{n}\left(a_{n}\vect v_{n}\left(k\vect r\right)+f_{n}\vect u_{n}\left(k\vect r\right)\right).\label{eq:power-Kristensson-E} +\end{equation} + +\end_inset + +Here +\begin_inset Formula $\eta_{0}=\sqrt{\mu_{0}/\varepsilon_{0}}$ +\end_inset + + is the wave impedance of free space and +\begin_inset Formula $\eta=\sqrt{\mu/\varepsilon}$ +\end_inset + + is the relative wave impedance of the medium. + +\end_layout + +\begin_layout Standard +The radiated power is then (2.28): +\begin_inset Formula +\[ +P=\frac{1}{2}\sum_{n}\left(\left|f_{n}\right|^{2}+\Re\left(f_{n}a_{n}^{*}\right)\right). +\] + +\end_inset + +The first term is obviously the power radiated away by the outgoing waves. + The second term must then be minus the power sucked by the scatterer from + the exciting wave. + If the exciting wave is plane, it gives us the extinction cross section +\begin_inset Formula +\[ +\sigma_{\mathrm{tot}}=-\frac{\sum_{n}\Re\left(f_{n}a_{n}^{*}\right)}{\varepsilon\varepsilon_{0}\left|E_{0}\right|^{2}} +\] + +\end_inset + + +\end_layout + +\begin_layout Subsection +Taylor +\end_layout + +\begin_layout Standard +Here I derive the radiated power in Taylor's convention by applying the + relations from subsection +\begin_inset CommandInset ref +LatexCommand ref +reference "sub:Kristensson-v-Taylor" + +\end_inset + + to the Kristensson's formulae (sect. + +\begin_inset CommandInset ref +LatexCommand ref +reference "sub:Radiated enenergy-Kristensson" + +\end_inset + +). +\end_layout + +\begin_layout Standard +Assume the external field decomposed as (here I use tildes even for the + expansion coefficients in order to avoid confusion with the +\begin_inset Formula $a_{n}$ +\end_inset + + in +\begin_inset CommandInset ref +LatexCommand eqref +reference "eq:power-Kristensson-E" + +\end_inset + +) +\begin_inset Formula +\[ +\vect E\left(\vect r,\omega\right)=\sum_{mn}\left[-i\left(\tilde{p}_{mn}\vect{\widetilde{N}}_{mn}^{(1)}+\tilde{q}_{mn}\widetilde{\vect M}_{mn}^{(1)}\right)+i\left(\tilde{a}_{mn}\widetilde{\vect N}_{mn}^{(3)}+\tilde{b}_{mn}\widetilde{\vect M}_{mn}^{(3)}\right)\right] +\] + +\end_inset + +(there is minus between the regular and outgoing part!). + The coefficients are related to those from +\begin_inset CommandInset ref +LatexCommand eqref +reference "eq:power-Kristensson-E" + +\end_inset + + as +\begin_inset Formula +\[ +\tilde{p}_{mn}=\frac{k\sqrt{\eta_{0}\eta}}{-i\sqrt{n(n+1)}}a_{2nm},\quad\tilde{q}_{mn}=\frac{k\sqrt{\eta_{0}\eta}}{-i\sqrt{n(n+1)}}a_{1nm}, +\] + +\end_inset + + +\begin_inset Formula +\[ +\tilde{a}_{mn}=\frac{k\sqrt{\eta_{0}\eta}}{i\sqrt{n(n+1)}}f_{2nm},\quad\tilde{b}_{mn}=\frac{k\sqrt{\eta_{0}\eta}}{i\sqrt{n(n+1)}}f_{1nm}. +\] + +\end_inset + +The radiated power is then +\begin_inset Formula +\[ +P=\frac{1}{2}\sum_{m,n}\frac{n\left(n+1\right)}{k^{2}\eta_{0}\eta}\left(\left|a_{mn}\right|^{2}+\left|b_{mn}\right|^{2}-\Re\left(a_{mn}p_{mn}^{*}\right)-\Re\left(b_{mn}q_{mn}^{*}\right)\right). +\] + +\end_inset + +If the exciting wave is a plane wave, the extinction cross section is +\begin_inset Formula +\[ +\sigma_{\mathrm{tot}}=\frac{1}{\varepsilon\varepsilon_{0}\left|E_{0}\right|^{2}k^{2}\eta_{0}\eta}\sum_{m,n}n(n+1)\left(\Re\left(a_{mn}p_{mn}^{*}\right)+\Re\left(b_{mn}q_{mn}^{*}\right)\right) +\] + +\end_inset + + +\end_layout + +\begin_layout Subsection +Jackson +\end_layout + +\begin_layout Standard +\begin_inset CommandInset citation +LatexCommand cite +after "(9.155)" +key "jackson_classical_1998" + +\end_inset + +: +\begin_inset Formula +\[ +P=\frac{Z_{0}}{2k^{2}}\sum_{l,m}\left[\left|a_{E}(l,m)\right|^{2}+\left|a_{M}(l,m)\right|^{2}\right] +\] + +\end_inset + + +\end_layout + \begin_layout Section Limit solutions \end_layout @@ -139,7 +900,7 @@ key "pustovit_plasmon-mediated_2010" \end_inset - + and Jackson (9.169) and around. \end_layout \begin_layout Subsection @@ -147,15 +908,166 @@ Near field limit \end_layout \begin_layout Chapter -Mie Theory +Single particle scattering and Mie theory +\end_layout + +\begin_layout Standard +The basic idea is simple. + For an exciting spherical wave (usually the regular wave in whatever convention +) of a given frequency +\begin_inset Formula $\omega$ +\end_inset + +, type +\begin_inset Formula $\zeta'$ +\end_inset + + (electric or magnetic), degree +\begin_inset Formula $l'$ +\end_inset + + and order +\begin_inset Formula $m'$ +\end_inset + +, the particle responds with waves from the complementary set (e.g. + outgoing waves), with the same frequency +\begin_inset Formula $\omega$ +\end_inset + +, but any type +\begin_inset Formula $\zeta$ +\end_inset + +, degree +\begin_inset Formula $l$ +\end_inset + + and order +\begin_inset Formula $m$ +\end_inset + +, in a way that the Maxwell's equations are satisfied, with the coefficients + +\begin_inset Formula $T_{l,m;l',m'}^{\zeta,\zeta'}(\omega)$ +\end_inset + +. + This yields one row in the scattering matrix (often called the +\begin_inset Formula $T$ +\end_inset + +-matrix) +\begin_inset Formula $T(\omega)$ +\end_inset + +, which fully characterizes the scattering properties of the particle (in + the linear regime, of course). + Analytical expression for the matrix is known for spherical scatterer, + otherwise it is computed numerically (using DDA, BEM or whatever). + So if we have the two sets of spherical wave functions +\begin_inset Formula $\vect f_{lm}^{J_{1},\zeta}$ +\end_inset + +, +\begin_inset Formula $\vect f_{lm}^{J_{2},\zeta}$ +\end_inset + + and the full +\begin_inset Quotes sld +\end_inset + +exciting +\begin_inset Quotes srd +\end_inset + + wave has electric field given as +\begin_inset Formula +\[ +\vect E_{\mathrm{inc}}=\sum_{\zeta'=\mathrm{E,M}}\sum_{l',m'}c_{l'm'}^{\zeta'}\vect f_{l'm'}^{\zeta'}, +\] + +\end_inset + +the +\begin_inset Quotes sld +\end_inset + +scattered +\begin_inset Quotes srd +\end_inset + + field will be +\begin_inset Formula +\[ +\vect E_{\mathrm{scat}}=\sum_{\zeta',l',m'}\sum_{\zeta,l,m}T_{l,m;l',m'}^{\zeta,\zeta'}c_{l'm'}^{\zeta'}\vect f_{lm}^{\zeta}, +\] + +\end_inset + +and the total field around the scaterer is +\begin_inset Formula $\vect E=\vect E_{\mathrm{ext}}+\vect E_{\mathrm{scat}}$ +\end_inset + +. \end_layout \begin_layout Section -Full version +Mie theory – full version +\end_layout + +\begin_layout Standard +\begin_inset Formula $T$ +\end_inset + +-matrix for a spherical particle is type-, degree- and order- diagonal, + that is, +\begin_inset Formula $T_{l',m';l,m}^{\zeta',\zeta}(\omega)=0$ +\end_inset + + if +\begin_inset Formula $l\ne l'$ +\end_inset + +, +\begin_inset Formula $m\ne m'$ +\end_inset + + or +\begin_inset Formula $\zeta\ne\zeta'$ +\end_inset + +. + Moreover, it does not depend on +\begin_inset Formula $m$ +\end_inset + +, so +\begin_inset Formula +\[ +T_{l,m;l',m'}^{\zeta,\zeta'}(\omega)=T_{l}^{\zeta}\left(\omega\right)\delta_{\zeta'\zeta}\delta_{l'l}\delta_{m'm} +\] + +\end_inset + +where for the usual choice +\begin_inset Formula $J_{1}=1,J_{2}=3$ +\end_inset + + +\begin_inset Formula +\begin{eqnarray*} +T_{l}^{E}\left(\omega\right) & = & TODO,\\ +T_{l}^{M}(\omega) & = & TODO. +\end{eqnarray*} + +\end_inset + + \end_layout \begin_layout Section -Long wave approximation +Long wave approximation for spherical nanoparticle \end_layout \begin_layout Standard @@ -170,6 +1082,140 @@ key "pustovit_plasmon-mediated_2010" and around. \end_layout +\begin_layout Section +Note on transforming T-matrix conventions +\end_layout + +\begin_layout Standard +T-matrix depends on the used conventions as well. + This is not apparent for the Mie case as the T-matrix for a sphere is +\begin_inset Quotes sld +\end_inset + +diagonal +\begin_inset Quotes srd +\end_inset + +. + But for other shapes, dipole incoming field can induce also higher-order + multipoles in the nanoparticle, etc. + The easiest way to determine the transformation properties is to write + down the total scattered electric field for both conventions in the form +\begin_inset Formula +\[ +\vect E_{\mathrm{scat}}=\sum_{n'}\sum_{n}T_{n'}^{n}c^{n'}\vect f_{n}=\sum_{n'}\sum_{n}\widetilde{T}_{n'}^{n}\widetilde{c}^{n'}\widetilde{\vect f}_{n} +\] + +\end_inset + +where we merged all the indices into single multiindex +\begin_inset Formula $n$ +\end_inset + + or +\begin_inset Formula $n'$ +\end_inset + +. + This way of writing immediately suggest how to transform the T-matrix into + the new convention if we know the transformation properties of the base + waves and expansion coefficients, as it reminds the notation used in geometry + – +\begin_inset Formula $c^{\alpha}$ +\end_inset + + are +\begin_inset Quotes sld +\end_inset + +vector coordinates +\begin_inset Quotes srd +\end_inset + +, +\begin_inset Formula $\vect f_{\alpha}$ +\end_inset + + are +\begin_inset Quotes sld +\end_inset + +base vectors +\begin_inset Quotes srd +\end_inset + +. + Obviously, T-matrix is then +\begin_inset Quotes sld +\end_inset + +tensor of type (1,1) +\begin_inset Quotes srd +\end_inset + +, and it transforms as vector coordinates (i.e. + wave expansion coefficients) in the +\begin_inset Formula $n$ +\end_inset + + (outgoing wave) indices and as form coordinates in the +\begin_inset Formula $n'$ +\end_inset + + (regular/illuminating wave) indices. + Form coordinates change in the same waves as base vectors +\end_layout + +\begin_layout Subsection +Kristensson to Taylor +\end_layout + +\begin_layout Standard +For instance, let us transform between from the Kristensson's to Taylor's + convention. + We know that the Taylor's base vectors are +\begin_inset Quotes sld +\end_inset + +larger +\begin_inset Quotes srd +\end_inset + +: +\begin_inset Formula $\widetilde{\vect N}_{ml}^{(3/1/4)}=\sqrt{l(l+1)}\left(\vect{u/v/w}\right)_{2lm}$ +\end_inset + + etc, so the coefficients must be smaller by the reciprocal factor, e.g. + +\begin_inset Formula $\tilde{a}_{ml}=f_{2lm}/\sqrt{l(l+1)}$ +\end_inset + + (now we assume that there are no other prefactors in the expansion of the + field, they are already included in the coefficients). + Then the T-matrix in the Taylor's convention (tilded) can be calculated + from the T-matrix in the Kristensson's convention as +\begin_inset Formula +\[ +\underbrace{\widetilde{T}_{\zeta'l'm'}^{\zeta lm}}_{\mbox{Taylor}}=\frac{\sqrt{l'(l'+1)}}{\sqrt{l(l+1)}}\underbrace{T_{\zeta'l'm'}^{\zeta lm}}_{\mbox{Krist.}}\,_{\leftarrow\mbox{illuminating}}^{\leftarrow\mbox{outgoing}}. +\] + +\end_inset + + +\end_layout + +\begin_layout Subsubsection +scuff-tmatrix output +\end_layout + +\begin_layout Standard +Indices of the outgoing wave (without primes) come first, illuminating regular + wave (with primes) second in the output files of scuff-tmatrix. + It seems that it at least in the electric part, the output of scuff-tmatrix + is equivalent to the Kristensson's convention. + Not sure whether it is also true for the E-M cross terms. +\end_layout + \begin_layout Chapter Green's functions \end_layout @@ -183,17 +1229,17 @@ Mie decomposition of Green's function for single nanoparticle \end_layout \begin_layout Chapter -Translation of spherical waves: shit gets insane +Translation of spherical waves: getting insane \end_layout \begin_layout Chapter -Multiple scattering: nice linear algebra born from all the shit +Multiple scattering: nice linear algebra born from all the mess \end_layout \begin_layout Standard \begin_inset CommandInset bibtex LatexCommand bibtex -bibfiles "dipdip" +bibfiles "Electrodynamics" options "plain" \end_inset diff --git a/dipdip.bib b/dipdip.bib index 28c46d2..46f379a 100644 --- a/dipdip.bib +++ b/dipdip.bib @@ -200,6 +200,38 @@ They observe exciton conductance (defined as loss of energy from the last molecu file = {Michael I. Mishchenko, Joachim W. Hovenier, Larry D. Travis-Light Scattering by Nonspherical Particles-Academic Press (1999).pdf:/home/necadam1/.zotero/zotero/9uf64zmd.default/zotero/storage/9HIG2UUN/Michael I. Mishchenko, Joachim W. Hovenier, Larry D. Travis-Light Scattering by Nonspherical Particles-Academic Press (1999).pdf:application/pdf} } +@article{goldstein_dipole-dipole_1997, + title = {Dipole-dipole interaction in optical cavities}, + volume = {56}, + url = {http://link.aps.org/doi/10.1103/PhysRevA.56.5135}, + doi = {10.1103/PhysRevA.56.5135}, + abstract = {At the most fundamental level of quantum electrodynamics, there is no such thing as two-body interactions between atoms. The potentials that describe these interactions are effective potentials resulting from a series of approximations whose validity depends on the precise situation at hand. Considering a one-dimensional geometry for simplicity, we discuss under which conditions the familiar form of the near-resonant dipole-dipole interaction is valid, paying particular attention to the effects of interatomic propagation of light, and to what extent it can be modified in the tailored electromagnetic environments provided by optical resonators. We find that once the atoms have established that they are inside a resonator, the dipole-dipole potential may or may not remain a useful concept, depending upon the strength of the atom-field interaction. In the weak-coupling regime, one finds a two-body dipole-dipole interaction that can be enhanced or inhibited by varying the atom-field detuning. In the strong-coupling regime, by contrast, the two-body dipole-dipole potential ceases to be meaningful.}, + number = {6}, + urldate = {2016-04-04}, + journal = {Physical Review A}, + author = {Goldstein, E. V. and Meystre, P.}, + month = dec, + year = {1997}, + pages = {5135--5146}, + file = {APS Snapshot:/home/necadam1/.zotero/zotero/9uf64zmd.default/zotero/storage/DSPRNX3I/PhysRevA.56.html:text/html;Full Text PDF:/home/necadam1/.zotero/zotero/9uf64zmd.default/zotero/storage/45J3SJBS/Goldstein and Meystre - 1997 - Dipole-dipole interaction in optical cavities.pdf:application/pdf} +} + +@article{dung_resonant_2002, + title = {Resonant dipole-dipole interaction in the presence of dispersing and absorbing surroundings}, + volume = {66}, + url = {http://link.aps.org/doi/10.1103/PhysRevA.66.063810}, + doi = {10.1103/PhysRevA.66.063810}, + abstract = {Within the framework of quantization of the macroscopic electromagnetic field, equations of motion and an effective Hamiltonian for treating both the resonant dipole-dipole interaction between two-level atoms and the resonant atom-field interaction are derived, which can suitably be used for studying the influence of arbitrary dispersing and absorbing material surroundings on these interactions. The theory is applied to the study of the transient behavior of two atoms that initially share a single excitation, with special emphasis on the role of the two competing processes of virtual- and real-photon exchange in the energy transfer between the atoms. In particular, it is shown that for weak atom-field interaction there is a time window, where the energy transfer follows a rate regime of the type obtained by ordinary second-order perturbation theory. Finally, it is shown that the resonant dipole-dipole interaction can change the singlet line of the emitted light to a doublet spectrum for weak atom-field interaction and the doublet spectrum to a triplet spectrum for strong atom-field interaction.}, + number = {6}, + urldate = {2016-04-04}, + journal = {Physical Review A}, + author = {Dung, Ho Trung and Knöll, Ludwig and Welsch, Dirk-Gunnar}, + month = dec, + year = {2002}, + pages = {063810}, + file = {APS Snapshot:/home/necadam1/.zotero/zotero/9uf64zmd.default/zotero/storage/P2VKBU4Q/PhysRevA.66.html:text/html;PhysRevA.66.063810.pdf:/home/necadam1/.zotero/zotero/9uf64zmd.default/zotero/storage/GVUTKSQQ/PhysRevA.66.063810.pdf:application/pdf} +} + @article{lehmberg_radiation_1970, title = {Radiation from an \${N}\$-{Atom} {System}. {II}. {Spontaneous} {Emission} from a {Pair} of {Atoms}}, volume = {2}, @@ -336,6 +368,16 @@ They observe exciton conductance (defined as loss of energy from the last molecu file = {APS Snapshot:/home/necadam1/.zotero/zotero/9uf64zmd.default/zotero/storage/HX4M79VA/PhysRevA.55.html:text/html;get (1).pdf:/home/necadam1/.zotero/zotero/9uf64zmd.default/zotero/storage/I5D2XQ7A/get (1).pdf:application/pdf} } +@book{smith_cython_????, + title = {Cython}, + isbn = {978-1-4919-0155-7}, + url = {http://shop.oreilly.com/product/0636920033431.do}, + abstract = {Build software that combines Python’s expressivity with the performance and control of C (and C++). It’s possible with Cython, the compiler and hybrid programming language used by foundational packages such as NumPy. In this practical...}, + urldate = {2016-03-24}, + author = {Smith, Kurt W.}, + file = {Kurt W. Smith-Cython - A guide for Python programmers-O'Reilly (2015).pdf:/home/necadam1/.zotero/zotero/9uf64zmd.default/zotero/storage/5D4IVZKZ/Kurt W. Smith-Cython - A guide for Python programmers-O'Reilly (2015).pdf:application/pdf;Snapshot:/home/necadam1/.zotero/zotero/9uf64zmd.default/zotero/storage/NJ4C4SSJ/0636920033431.html:text/html} +} + @article{nicolosi_dissipation-induced_2004, title = {Dissipation-induced stationary entanglement in dipole-dipole interacting atomic samples}, volume = {70}, @@ -665,6 +707,22 @@ and 3.4 (p, 389+) too.}, file = {APS Snapshot:/home/necadam1/.zotero/zotero/9uf64zmd.default/zotero/storage/WQJ5AIGU/PhysRevLett.111.html:text/html;PhysRevLett.111.166802.pdf:/home/necadam1/.zotero/zotero/9uf64zmd.default/zotero/storage/ZJZXN4VA/PhysRevLett.111.166802.pdf:application/pdf} } +@article{pirruccio_coherent_2016, + title = {Coherent {Control} of the {Optical} {Absorption} in a {Plasmonic} {Lattice} {Coupled} to a {Luminescent} {Layer}}, + volume = {116}, + url = {http://link.aps.org/doi/10.1103/PhysRevLett.116.103002}, + doi = {10.1103/PhysRevLett.116.103002}, + abstract = {We experimentally demonstrate the coherent control, i.e., phase-dependent enhancement and suppression, of the optical absorption in an array of metallic nanoantennas covered by a thin luminescent layer. The coherent control is achieved by using two collinear, counterpropagating, and phase-controlled incident waves with wavelength matching the absorption spectrum of dye molecules coupled to the array. Symmetry arguments shed light on the relation between the relative phase of the incident waves and the excitation efficiency of the optical resonances of the system. This coherent control is associated with a phase-dependent distribution of the electromagnetic near fields in the structure which enables a significant reduction of the unwanted dissipation in the metallic structures.}, + number = {10}, + urldate = {2016-03-29}, + journal = {Physical Review Letters}, + author = {Pirruccio, Giuseppe and Ramezani, Mohammad and Rodriguez, Said Rahimzadeh-Kalaleh and Rivas, Jaime Gómez}, + month = mar, + year = {2016}, + pages = {103002}, + file = {APS Snapshot:/home/necadam1/.zotero/zotero/9uf64zmd.default/zotero/storage/IESGABCT/PhysRevLett.116.html:text/html;PhysRevLett.116.103002.pdf:/home/necadam1/.zotero/zotero/9uf64zmd.default/zotero/storage/ABU3ENFF/PhysRevLett.116.103002.pdf:application/pdf} +} + @article{kurizki_suppression_1988, title = {Suppression of {Molecular} {Interactions} in {Periodic} {Dielectric} {Structures}}, volume = {61}, @@ -880,6 +938,38 @@ and 3.4 (p, 389+) too.}, file = {ACS Full Text PDF w/ Links:/home/necadam1/.zotero/zotero/9uf64zmd.default/zotero/storage/UVH7KEPI/Zeman and Schatz - 1987 - An accurate electromagnetic theory study of surfac.pdf:application/pdf;ACS Full Text Snapshot:/home/necadam1/.zotero/zotero/9uf64zmd.default/zotero/storage/GJFAX4IF/j100287a028.html:text/html} } +@article{fredkin_resonant_2003, + title = {Resonant {Behavior} of {Dielectric} {Objects} ({Electrostatic} {Resonances})}, + volume = {91}, + issn = {0031-9007, 1079-7114}, + url = {http://adsabs.harvard.edu/abs/2003PhRvL..91y3902F}, + doi = {10.1103/PhysRevLett.91.253902}, + language = {en}, + number = {25}, + urldate = {2014-03-21}, + journal = {Physical Review Letters}, + author = {Fredkin, D. and Mayergoyz, I.}, + month = dec, + year = {2003}, + file = {NASA ADS\: Resonant Behavior of Dielectric Objects (Electrostatic Resonances):/home/necadam1/.zotero/zotero/9uf64zmd.default/zotero/storage/DIW3Q4DW/2003PhRvL..html:text/html} +} + +@article{gruner_green-function_1996, + title = {Green-function approach to the radiation-field quantization for homogeneous and inhomogeneous {Kramers}-{Kronig} dielectrics}, + volume = {53}, + url = {http://link.aps.org/doi/10.1103/PhysRevA.53.1818}, + doi = {10.1103/PhysRevA.53.1818}, + abstract = {A quantization scheme for the radiation field in dispersive and absorptive linear dielectrics is developed, which applies to both bulk material and multilayer dielectric structures. Starting from the phenomenological Maxwell equations, where the properties of the dielectric are described by a permittivity consistent with the Kramers-Kronig relations, an expansion of the field operators is performed that is based on the Green function of the classical Maxwell equations and preserves the equal-time canonical field commutation relations. In particular, in frequency intervals with approximately vanishing absorption the concept of quantization through mode expansion for dispersive dielectrics is recognized. The theory further reveals that weak absorption gives rise to space-dependent mode operators that spatially evolve according to quantum Langevin equations in the space domain. To illustrate the applicability of the theory to inhomogeneous structures, the quantization of the radiation field in a dispersive and absorptive one-interface dielectric is performed. © 1996 The American Physical Society.}, + number = {3}, + urldate = {2016-05-13}, + journal = {Physical Review A}, + author = {Gruner, T. and Welsch, D.-G.}, + month = mar, + year = {1996}, + pages = {1818--1829}, + file = {APS Snapshot:/home/necadam1/.zotero/zotero/9uf64zmd.default/zotero/storage/KNE8THEZ/PhysRevA.53.html:text/html;PhysRevA.53.1818.pdf:/home/necadam1/.zotero/zotero/9uf64zmd.default/zotero/storage/7PDC6W3U/PhysRevA.53.1818.pdf:application/pdf} +} + @book{rother_electromagnetic_2014, address = {Berlin, Heidelberg}, series = {Springer {Series} in {Optical} {Sciences}}, @@ -1088,6 +1178,24 @@ and 3.4 (p, 389+) too.}, file = {(Springer Theses) Jonathan M. Taylor (auth.)-Optical Binding Phenomena_ Observations and Mechanisms -Springer-Verlag Berlin Heidelberg (2011).pdf:/home/necadam1/.zotero/zotero/9uf64zmd.default/zotero/storage/7XKKCD9X/(Springer Theses) Jonathan M. Taylor (auth.)-Optical Binding Phenomena_ Observations and Mechanisms -Springer-Verlag Berlin Heidelberg (2011).pdf:application/pdf} } +@article{de_vries_point_1998, + title = {Point scatterers for classical waves}, + volume = {70}, + url = {http://link.aps.org/doi/10.1103/RevModPhys.70.447}, + doi = {10.1103/RevModPhys.70.447}, + abstract = {The authors present a closed formulation of resonant point scatterers for classical-wave propagation problems. A Green’s-function approach is employed in which all the small-distance singularities are regularized. Application of point scatterers considerably simplifies multiple-scattering calculations needed, for instance, for understanding the optical properties of dense cold gases and optical lattices. In the case of the vector description of light, it is shown that two different regularization parameters are required in order to obtain physically meaningful results. One parameter is related to the physical size of the pointlike scattering particle, while the other is connected to its dynamic properties. All parameters involved are defined in terms of physical observables leading to a complete and self-consistent treatment. The applicability of the point-scatterer model to several physical models is demonstrated. We calculate the local density of states of waves in the presence of one resonant point scatterer. For the vector case, the bare polarizability is shown to enter the local density of states. For a collection of resonant point dipoles, the Lorentz-Lorenz relation for the dielectric constant is derived for cubic lattices and for disordered arrangements.}, + number = {2}, + urldate = {2016-02-15}, + journal = {Reviews of Modern Physics}, + author = {de Vries, Pedro and van Coevorden, David V. and Lagendijk, Ad}, + month = apr, + year = {1998}, + pages = {447--466}, + annote = {Very nice article about light scattering on point dipole. +Renormalisation, relation between the regularisation parameters and measurable physical parameters of the scatterer.}, + file = {APS Snapshot:/home/necadam1/.zotero/zotero/9uf64zmd.default/zotero/storage/7KCG8EN3/RevModPhys.70.html:text/html;RevModPhys.70.447.pdf:/home/necadam1/.zotero/zotero/9uf64zmd.default/zotero/storage/CGE7Z2CZ/RevModPhys.70.447.pdf:application/pdf} +} + @article{epton_multipole_1995, title = {Multipole {Translation} {Theory} for the {Three}-{Dimensional} {Laplace} and {Helmholtz} {Equations}}, volume = {16}, diff --git a/qpms/qpms_p.py b/qpms/qpms_p.py index 473aba6..4d907e9 100644 --- a/qpms/qpms_p.py +++ b/qpms/qpms_p.py @@ -215,6 +215,9 @@ def zJn(n, z, J=1): # The following 4 funs have to be refactored, possibly merged + +# FIXME: this can be expressed simply as: +# $$ -\frac{1}{2}\sqrt{\frac{2n+1}{4\pi}n\left(n+1\right)}(\delta_{m,1}+\delta_{m,-1}) $$ def π̃_zerolim(nmax): # seems OK """ lim_{θ→ 0-} π̃(cos θ) @@ -249,6 +252,8 @@ def π̃_pilim(nmax): # Taky OK, jen to možná není kompatibilní se vzorečky π̃_y = prenorm * π̃_y return π̃_y +# FIXME: this can be expressed simply as +# $$ -\frac{1}{2}\sqrt{\frac{2n+1}{4\pi}n\left(n+1\right)}(\delta_{m,1}-\delta_{m,-1}) $$ def τ̃_zerolim(nmax): """ lim_{θ→ 0-} τ̃(cos θ) @@ -419,7 +424,7 @@ def Ã(m,n,μ,ν,kdlj,θlj,φlj,r_ge_d,J): +math.lgamma(n+ν+1) - math.lgamma(2*(n+ν)+1)) presum = math.exp(exponent) presum = presum * np.exp(1j*(μ-m)*φlj) * (-1)**m * 1j**(ν+n) / (4*n) - qmax = floor(q_max(-m,n,μ,ν)) #nemá tu být +m? + qmax = math.floor(q_max(-m,n,μ,ν)) #nemá tu být +m? q = np.arange(qmax+1, dtype=int) # N.B. -m !!!!!! a1q = a_q(-m,n,μ,ν) # there is redundant calc. of qmax @@ -451,7 +456,7 @@ def B̃(m,n,μ,ν,kdlj,θlj,φlj,r_ge_d,J): presum = math.exp(exponent) presum = presum * np.exp(1j*(μ-m)*φlj) * (-1)**m * 1j**(ν+n+1) / ( (4*n)*(n+1)*(n+m+1)) - Qmax = floor(q_max(-m,n+1,μ,ν)) + Qmax = math.floor(q_max(-m,n+1,μ,ν)) q = np.arange(Qmax+1, dtype=int) if (μ == ν): # it would disappear in the sum because of the factor (ν-μ) anyway ã2q = 0 @@ -500,6 +505,7 @@ def mie_coefficients(a, nmax, #ω, ε_i, ε_e=1, J_ext=1, J_scat=3 FIXME test the magnetic case TODO description + RH concerns the N ("electric") part, RV the M ("magnetic") part # Parameters @@ -670,6 +676,14 @@ def G_Mie_scat_cart(source_cart, dest_cart, a, nmax, k_i, k_e, μ_i=1, μ_e=1, J return G_Mie_scat_precalc_cart_new(source_cart, dest_cart, RH, RV, a, nmax, k_i, k_e, μ_i, μ_e, J_ext, J_scat) +#TODO +def cross_section_Mie_precalc(): + pass + +def cross_section_Mie(a, nmax, k_i, k_e, μ_i, μ_e,): + pass + + # In[9]: # From PRL 112, 253601 (1) diff --git a/setup.py b/setup.py index b34257b..7507c5b 100644 --- a/setup.py +++ b/setup.py @@ -15,7 +15,8 @@ setup(name='qpms', version = "0.1", packages=['qpms'], # setup_requires=['setuptools_cython'], - install_requires=['cython>=0.21','quaternion','spherical_functions'], + install_requires=['cython>=0.21','quaternion','spherical_functions','py_gmm'], + dependency_links=['https://github.com/texnokrates/py_gmm','https://github.com/moble/quaternion','https://github.com/moble/spherical_functions'], ext_modules=[qpms_c], cmdclass = {'build_ext': build_ext}, ) diff --git a/worknotes.lyx b/worknotes.lyx index f657533..c0f7203 100644 --- a/worknotes.lyx +++ b/worknotes.lyx @@ -1717,12 +1717,34 @@ Numerics: \end_layout \begin_layout Section -TODO +Misc \end_layout \begin_layout Itemize -Päivi's suggestion: suppress the dipole and let it interact only with the - higher multipoles. +The +\begin_inset Quotes eld +\end_inset + +zero limits +\begin_inset Quotes erd +\end_inset + + of +\begin_inset Formula $\tilde{\pi},\tilde{\tau}$ +\end_inset + + functions in Taylor's normalisation can be expressed as +\lang finnish + +\begin_inset Formula +\begin{eqnarray*} +\lim_{\theta\to0}\tilde{\pi}_{mn}\left(\cos\theta\right) & = & -\frac{1}{2}\sqrt{\frac{2n+1}{4\pi}n\left(n+1\right)}(\delta_{m,1}+\delta_{m,-1})\\ +\lim_{\theta\to0}\tilde{\tau}_{mn}\left(\cos\theta\right) & = & -\frac{1}{2}\sqrt{\frac{2n+1}{4\pi}n\left(n+1\right)}(\delta_{m,1}-\delta_{m,-1}) +\end{eqnarray*} + +\end_inset + + \end_layout \begin_layout Standard