Top Qs
Timeline
Chat
Perspective

Discrete dipole approximation

From Wikipedia, the free encyclopedia

Discrete dipole approximation
Remove ads

The discrete dipole approximation (DDA), also known as the coupled dipole approximation, is a numerical method for computing the scattering and absorption of electromagnetic radiation by particles of arbitrary shape and composition. The method represents a continuum target as a finite array of small, polarizable dipoles, and solves for their interactions with the incident field and with each other. DDA can handle targets with inhomogeneous composition and anisotropic material properties, as well as periodic structures. It is widely applied in fields such as nanophotonics, radar scattering, aerosol physics, biomedical optics, and astrophysics.


Thumb
Magnitude of the electric field strength (colored) and the Poynting vector (black arrows) in the near field of the vertically oriented dipole in the image plane. Blue/red colors indicate an electric field oriented downwards/upwards.
Remove ads

Basic concepts

Summarize
Perspective

The basic idea of the DDA was introduced in 1964 by DeVoe[1] 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[2] 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.[3][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 cubical 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.

There are several reviews of DDA method. [7][6][8][9]

The method was improved by Draine, Flatau, and Goodman, who applied the fast Fourier transform to solve fast convolution problems arising in the discrete dipole approximation (DDA). This allowed for the calculation of scattering by large targets. They distributed an open-source code DDSCAT.[7][10] There are now several DDA implementations,[6] extensions to periodic targets,[11] and particles placed on or near a plane substrate.[12][13] Comparisons with exact techniques have also been published.[14] Other aspects, such as the validity criteria of the discrete dipole approximation, were published.[15] The DDA was also extended to employ rectangular or cuboid dipoles,[16] which are more efficient for highly oblate or prolate particles.


Remove ads

Theory

Summarize
Perspective

In the discrete dipole approximation, a target object is represented as a finite array of N point dipoles located at positions (). The polarization vector of each dipole is related to the local electric field at that dipole by its polarizability tensor :

In anisotropic case (diagonal polarizability)

where is diagonal

This leads to componentwise relations:

For isotropic materials, , so

.

The local electric field acting on the j‑th dipole is given by the sum of the incident field and the fields radiated by all other dipoles:

Here, is the dyadic Green's function describing the field at position due to a unit dipole at the origin.

Dyadic Green’s function

In the discrete dipole approximation, the electric field at dipole due to dipole is given by

where , , and . Here is the identity matrix and is the vacuum wavenumber.

Define

The dyadic Green’s function is:

Notice that it is symmetric: .

Here is the displacement vector from dipole to dipole , is the distance between them, and is the unit vector pointing from to . The components of are defined as:

Polarizability

In the discrete dipole approximation, the electromagnetic response of a target is modeled by replacing the continuous material with a finite array of point dipoles. Each dipole represents a small volume of the material and acts as a polarizable unit that interacts with both the incident field and the fields radiated by all other dipoles. The key parameter that describes how each dipole responds to the local electric field is its polarizability . For a homogeneous material, the polarizability of a dipole is determined by the material’s complex dielectric function , which depends on the wavelength of light in vacuum. The dielectric function is related to the complex refractive index through . The goal in DDA is to assign to each dipole a polarizability such that the array of dipoles reproduces, as accurately as possible, the scattering and absorption behavior of the original continuous medium. For isotropic materials, a common starting point is the Clausius–Mossotti relation, which connects the polarizability to the dielectric function:

In the discrete dipole approximation, the total volume of the target is divided into small cubic cells of volume , where is the lattice spacing. The Clausius–Mossotti polarizability for each dipole is

where is the relative permittivity of the material at the dipole’s position. The dipole volume is constant across all dipoles.

This formula assumes that each dipole occupies a volume embedded in an otherwise uniform dielectric medium. In most implementations of DDA the formulation is expressed in Gaussian units (CGS). In these units, the polarizability has dimensions of volume (cm³). In the discrete dipole approximation, the total volume of the target is divided into small cubic cells of volume , where is the lattice spacing and is the total number of dipoles. The total target volume is thus

To improve the accuracy of the method various corrections to are applied. These include: the lattice dispersion relation (LDR) polarizability (Draine & Goodman, 1993), which adjusts to ensure that the dispersion relation of an infinite lattice of dipoles matches that of the continuous material; the radiative reaction (RR) correction, which compensates for the fact that each dipole radiates energy and is influenced by its own radiation field.


Size parameter

The size parameter is a dimensionless quantity used in scattering theory to characterize the size of a particle relative to the wavelength of the incident light. For a sphere, it is defined as:

where: is the size parameter (dimensionless), is the radius of the sphere, is the wavelength of light in vacuum,

is the wavenumber.

In case of a sphere, the size parameter determines the scattering regime:

  • If , [Rayleigh scattering|Rayleigh scattering]] dominates.
  • If , the scattering is in the regime of Mie scattering.
  • If , the geometric optics approximation becomes valid.



Effective radius and dipole discretization

For nonspherical targets with the same volume as a sphere, the effective radius is often used in place of , with:

where: is the total number of dipoles, is the dipole spacing, is the total volume represented by the dipoles. This gives effective size parameter

One convenient trick in certain DDA accuracy tests is to define wavelength as , in such a case effective radius is the same as effective size parameter.

Dipole-scale size parameter

Each polarizable point (dipole) occupies a cubic volume with side length . Analogous to the global size parameter used for whole particles, one can define a local size parameter for each dipole:

This local parameter quantifies the ratio of the dipole size to the wavelength of light inside the material. For the DDA to be accurate, the field should vary slowly over the size of each dipole. This condition is satisfied when:

This ensures that each dipole is optically small, fields vary slowly over the dipole and the polarizability formula used for each dipole is accurate. Notice that a similar parameter plays a crucial role in the anomalous diffraction theory of van de Hulst, where the total phase shift experienced by light rays traveling through or around the particle is given by:

This describes the optical path difference introduced by the particle (or in the case of DDA by a dipole).


Explicit Matrix Form of the DDA System

The Discrete Dipole Approximation (DDA) linear system is expressed as:

Let N be the number of dipoles. Each dipole has a polarization vector . The total system is a matrix equation of size :

Each block is a complex matrix, defined by:

where: is the inverse polarizability tensor, is the dyadic Green’s tensor for interaction between dipoles and , are the dipole polarization and incident electric field at dipole , respectively.

Typically dipoles are arranged on a regular grid. This implies translational invariance:

Because , the matrix is symmetric:

As a result, the interaction matrix is a block Toeplitz symmetric matrix.

Remove ads

Fast Fourier Transform for fast convolution calculations

The Fast Fourier Transform (FFT) method was introduced in 1991 by Goodman, Draine, and Flatau[17] for the discrete dipole approximation. They utilized a 3D FFT GPFA written by Clive Temperton. The interaction matrix was extended to twice its original size to incorporate negative lags by mirroring and reversing the interaction matrix. Several variants have been developed since then. Barrowes, Teixeira, and Kong[18] in 2001 developed a code that uses block reordering, zero padding, and a reconstruction algorithm, claiming minimal memory usage. McDonald, Golden, and Jennings[19] in 2009 used a 1D FFT code and extended the interaction matrix in the x, y, and z directions of the FFT calculations, suggesting memory savings due to this approach. This variant was also implemented in the MATLAB 2021 code by Shabaninezhad and Ramakrishna[20]. Other techniques to accelerate convolutions have been suggested in a general context[21][22] along with faster evaluations of Fast Fourier Transforms arising in DDA problem solvers.

Conjugate gradient iteration schemes and preconditioning

Some of the early calculations of the polarization vector were based on direct inversion[2] and the implementation of the conjugate gradient method by Petravic and Kuo-Petravic.[23] Subsequently, many other conjugate gradient methods have been tested.[24] Advances in the preconditioning of linear systems of equations arising in the DDA setup have also been reported.[25]

Thermal discrete dipole approximation

Thermal discrete dipole approximation is an extension of the original DDA to simulations of near-field heat transfer between 3D arbitrarily-shaped objects.[26][27]

Discrete dipole approximation codes

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). There are some published comparisons of existing DDA codes.[14]

Remove ads

See also

References

Loading related searches...

Wikiwand - on

Seamless Wikipedia browsing. On steroids.

Remove ads