Get Discrete Dipole Approximation essential facts below. View Videos or join the Discrete Dipole Approximation discussion. Add Discrete Dipole Approximation to your PopFlock.com topic list for future reference or share this resource on social media.
Discrete Dipole Approximation
In the discrete dipole approximation, a larger object is approximated in terms of discrete radiating electric dipoles.
Discrete dipole approximation (DDA), also known as coupled dipole approximation,^{[1]} is a method for computing scattering of radiation by particles of arbitrary shape and by periodic structures. Given a target of arbitrary geometry, one seeks to calculate its scattering and absorption properties by an approximation of the continuum target by a finite array of small polarizabledipoles. This technique is used in a variety of applications including nanophotonics, radar scattering, aerosol physics and astrophysics.
Basic concepts
The basic idea of the DDA was introduced in 1964 by DeVoe^{[2]} who applied it to study the optical properties of molecular aggregates; retardation effects were not included, so DeVoe's treatment was limited to aggregates that were small compared with the wavelength. The DDA, including retardation effects, was proposed in 1973 by Purcell and Pennypacker^{[3]}
who used it to study interstellar dust grains. Simply stated, the DDA is an approximation of the continuum target by a finite array of polarizable points. The points acquire dipole moments in response to the local electric field. The dipoles interact with one another via their electric fields, so the DDA is also sometimes referred to as the coupled dipole approximation.^{[1]}^{[4]}
Nature provides the physical inspiration for the DDA - in 1909 Lorentz^{[5]}
showed that the dielectric properties of a substance could be directly related to the polarizabilities of the individual atoms of which it was composed, with a particularly simple and exact relationship, the Clausius-Mossotti relation (or Lorentz-Lorenz), when the atoms are located on a cubic lattice. We may expect that, just as a continuum representation of a solid is appropriate on length scales that are large compared with the interatomic spacing, an array of polarizable points can accurately approximate the response of a continuum target on length scales that are large compared with the interdipole separation.
For a finite array of point dipoles the scattering problem may be solved exactly, so the only approximation that is present in the DDA is the replacement of the continuum target by an array of N-point dipoles. The replacement requires specification of both the geometry (location of the dipoles) and the dipole polarizabilities. For monochromatic incident waves the self-consistent solution for the oscillating dipole moments may be found; from these the absorption and scattering cross sections are computed. If DDA solutions are obtained for two independent polarizations of the incident wave, then the complete amplitude scattering matrix can be determined.
Alternatively, the DDA can be derived from volume integral equation for the electric field.^{[6]} This highlights that the approximation of point dipoles is equivalent to that of discretizing the integral equation, and thus decreases with decreasing dipole size.
With the recognition that the polarizabilities may be tensors, the DDA can readily be applied to anisotropic materials. The extension of the DDA to treat materials with nonzero magnetic susceptibility is also straightforward, although for most applications magnetic effects are negligible.
Extensions
The method was improved by Draine, Flatau, and Goodman who applied fast Fourier transform to calculate convolution problem arising in the DDA which allowed to calculate scattering by large targets. They distributed discrete dipole approximation open source code DDSCAT.^{[7]}^{[8]}
There are now several DDA implementations,^{[6]} extensions to periodic targets^{[9]} and particles placed on or near a plane substrate.^{[10]}^{[11]} and comparisons with exact technique were published.^{[12]}
Other aspects such as the validity criteria of the discrete dipole approximation^{[13]} was published. The DDA was also extended to employ rectangular or cuboid dipoles ^{[14]} which is more efficient for highly oblate or prolate particles.
Discrete dipole approximation codes
There are reviews^{[7]}^{[6]} as well as published comparison of existing codes.^{[12]}
Most of the codes apply to arbitrary-shaped inhomogeneous nonmagnetic particles and particle systems in free space or homogeneous dielectric host medium. The calculated quantities typically include the Mueller matrices, integral cross-sections (extinction, absorption, and scattering), internal fields and angle-resolved scattered fields (phase function).
General-purpose open-source DDA codes
These codes typically use regular grids (cubic or rectangular cuboid), conjugate gradient method to solve large system of linear equations, and FFT-acceleration of the matrix-vector products which uses convolution theorem. Complexity of this approach is almost linear in number of dipoles for both time and memory.^{[6]}
Implements fast and rigorous consideration of a plane substrate, and allows rectangular-cuboid voxels for highly oblate or prolate particles. Can also calculate emission (decay-rate) enhancement of point emitters.Near-fields calculation is not very efficient. Uses Message Passing Interface (MPI) parallelization and can run on GPU (OpenCL).
Also calculates near fields and material absorption. Named differently, but the algorithms are very similar to the ones used in the mainstream DDA.
Specialized DDA codes
These list include codes that do not qualify for the previous section. The reasons may include the following: source code is not available, FFT acceleration is absent or reduced, the code focuses on specific applications not allowing easy calculation of standard scattering quantities.
Name
Authors
References
Language
Updated
Features
DDSURF, DDSUB, DDFILM
Schmehl, Nebeker, and Zhang
^{[10]}^{[23]}^{[24]}
Fortran
2008
Rigorous handling of semi-infinite substrate and finite films (with arbitrary particle placement), but only 2D FFT acceleration is used.
DDMM
Mackowski
^{[25]}
Fortran
2002
Calculates T-matrix, which can then be used to efficiently calculate orientation-averaged scattering properties.
Simulates electron-energy loss spectroscopy and cathodoluminescence. Handles substrate through image approximation, but no FFT acceleration is used.
T-DDA
Edalatpour
^{[30]}
Fortran
2015
Simulates near-field radiative heat transfer. The computational bottleneck is direct matrix inversion (no FFT acceleration is used). Uses OpenMP and MPI parallelization.
Gallery of shapes
Scattering by periodic structures such as slabs, gratings, of periodic cubes placed on a surface, can be solved in the discrete dipole approximation.
Scattering by infinite object (such as cylinder) can be solved in the discrete dipole approximation.
^ ^{a}^{b}Singham, Shermila B.; Salzman, Gary C. (1986). "Evaluation of the scattering matrix of an arbitrary particle using the coupled dipole approximation". The Journal of Chemical Physics. AIP Publishing. 84 (5): 2658-2667. doi:10.1063/1.450338. ISSN0021-9606.
^DeVoe, Howard (1964-07-15). "Optical Properties of Molecular Aggregates. I. Classical Model of Electronic Absorption and Refraction". The Journal of Chemical Physics. AIP Publishing. 41 (2): 393-400. doi:10.1063/1.1725879. ISSN0021-9606.
^
E. M. Purcell; C. R. Pennypacker (1973). "Scattering and absorption of light by nonspherical dielectric grains". Astrophysical Journal. 186: 705. Bibcode:1973ApJ...186..705P. doi:10.1086/152538.
^Singham, Shermila Brito; Bohren, Craig F. (1987-01-01). "Light scattering by an arbitrary particle: a physical reformulation of the coupled dipole method". Optics Letters. The Optical Society. 12 (1): 10-12. doi:10.1364/ol.12.000010. ISSN0146-9592.
^H. A. Lorentz, Theory of Electrons (Teubner, Leipzig, 1909)
^Chaumet, Patrick C.; Rahmani, Adel; Bryant, Garnett W. (2003-04-02). "Generalization of the coupled dipole method to periodic structures". Physical Review B. American Physical Society (APS). 67 (16): 165404. arXiv:physics/0305051. doi:10.1103/physrevb.67.165404. ISSN0163-1829.
^ ^{a}^{b}Schmehl, Roland; Nebeker, Brent M.; Hirleman, E. Dan (1997-11-01). "Discrete-dipole approximation for scattering by features on surfaces by means of a two-dimensional fast Fourier transform technique". Journal of the Optical Society of America A. The Optical Society. 14 (11): 3026-3036. doi:10.1364/josaa.14.003026. ISSN1084-7529.
^Groth, Samuel P and Polimeridis, Athanasios G and White, Jacob K (2020). "Accelerating the discrete dipole approximation via circulant preconditioning". Journal of Quantitative Spectroscopy and Radiative Transfer. 240: 106689.CS1 maint: multiple names: authors list (link)
^M. Zimmermann; A. Tausendfreund; S. Patzelt; G. Goch; S. Kieß; M. Z. Shaikh; M. Gregoire; S. Simon (2012). "In-process measuring procedure for sub-100 nm structures". J. Laser Appl. 24 (4): 042010. Bibcode:2012JLasA..24d2010Z. doi:10.2351/1.4719936.
^B. M. Nebeker (1998). Modeling of light scattering from features above and below surfaces using the discrete-dipole approximation (PhD). Tempe, AZ, USA: Arizona State University.
^V. L. Y. Loke; P. M. Mengüç; Timo A. Nieminen (2011). "Discrete dipole approximation with surface interaction: Computational toolbox for MATLAB". J. Quant. Spectrosc. Radiat. Transfer. 112 (11): 1711-1725. Bibcode:2011JQSRT.112.1711L. doi:10.1016/j.jqsrt.2011.03.012.
^N. W. Bigelow; A. Vaschillo; V. Iberi; J. P. Camden; D. J. Masiello (2012). "Characterization of the electron- and photon-driven plasmonic excitations of metal nanorods". ACS Nano. 6 (8): 7497-7504. doi:10.1021/nn302980u. PMID22849410.
^N. Geuquet; L. Henrard (2010). "EELS and optical response of a noble metal nanoparticle in the frame of a discrete dipole approximation". Ultramicroscopy. 110 (8): 1075-1080. doi:10.1016/j.ultramic.2010.01.013.