Permittivity

In electromagnetism, the absolute permittivity, often simply called permittivity and denoted by the Greek letter ε (epsilon), is a measure of the electric polarizability of a dielectric. A material with high permittivity polarizes more in response to an applied electric field than a material with low permittivity, thereby storing more energy in the material. In electrostatics, the permittivity plays an important role in determining the capacitance of a capacitor.

A dielectric medium showing orientation of charged particles creating polarization effects. Such a medium can have a lower ratio of electric flux to charge (more permittivity) than empty space

In the simplest case, the electric displacement field D resulting from an applied electric field E is

${\displaystyle \mathbf {D} =\varepsilon \mathbf {E} .}$

More generally, the permittivity is a thermodynamic function of state.[1] It can depend on the frequency, magnitude, and direction of the applied field. The SI unit for permittivity is farad per meter (F/m).

The permittivity is often represented by the relative permittivity ε_r which is the ratio of the absolute permittivity ε and the vacuum permittivity ε₀

${\displaystyle \kappa =\varepsilon _{\mathrm {r} }={\frac {\varepsilon }{\varepsilon _{0}}}}$.

This dimensionless quantity is also often and ambiguously referred to as the permittivity. Another common term encountered for both absolute and relative permittivity is the dielectric constant which has been deprecated in physics and engineering[2] as well as in chemistry.[3]

By definition, a perfect vacuum has a relative permittivity of exactly 1 whereas at STP, air has a relative permittivity of κ_air = 1.0006.

Relative permittivity is directly related to electric susceptibility (χ) by

${\displaystyle \chi =\kappa -1}$

otherwise written as

${\displaystyle \varepsilon =\varepsilon _{\mathrm {r} }\varepsilon _{0}=(1+\chi )\varepsilon _{0}}$

Units

The standard SI unit for permittivity is farad per meter (F/m or F·m⁻¹).[4]

${\displaystyle {\frac {\text{F}}{\text{m}}}={\frac {\text{C}}{{\text{V}}{\cdot }{\text{m}}}}={\frac {{\text{C}}^{2}}{{\text{N}}{\cdot }{\text{m}}^{2}}}={\frac {{\text{C}}^{2}{\cdot }{\text{s}}^{2}}{{\text{kg}}{\cdot }{\text{m}}^{3}}}}$

Explanation

In electromagnetism, the electric displacement field D represents the distribution of electric charges in a given medium resulting from the presence of an electric field E. This distribution includes charge migration and electric dipole reorientation. Its relation to permittivity in the very simple case of linear, homogeneous, isotropic materials with "instantaneous" response to changes in electric field is:

${\displaystyle \mathbf {D} =\varepsilon \mathbf {E} }$

where the permittivity ε is a scalar. If the medium is anisotropic, the permittivity is a second rank tensor.

In general, permittivity is not a constant, as it can vary with the position in the medium, the frequency of the field applied, humidity, temperature, and other parameters. In a nonlinear medium, the permittivity can depend on the strength of the electric field. Permittivity as a function of frequency can take on real or complex values.

In SI units, permittivity is measured in farads per meter (F/m or A²·s⁴·kg⁻¹·m⁻³). The displacement field D is measured in units of coulombs per square meter (C/m²), while the electric field E is measured in volts per meter (V/m). D and E describe the interaction between charged objects. D is related to the charge densities associated with this interaction, while E is related to the forces and potential differences.

Vacuum permittivity

The vacuum permittivity ε₀ (also called permittivity of free space or the electric constant) is the ratio D/E in free space. It also appears in the Coulomb force constant,

${\displaystyle k_{\text{e}}={\frac {1}{4\pi \varepsilon _{0}}}}$

Its value is[5]

${\displaystyle \varepsilon _{0}{\stackrel {\mathrm {def} }{=}}\ {\frac {1}{c_{0}^{2}\mu _{0}}}={\frac {1}{35\,950\,207\,149.472\,7056\pi }}{\text{ F/m}}\approx 8.854\,187\,8176\ldots \times 10^{-12}{\text{ F/m }}}$

where

• c₀ is the speed of light in free space,[lower-alpha 1]
• µ₀ is the vacuum permeability.

The constants c₀ and μ₀ were defined in SI units to have exact numerical values until redefinition of SI units in 2019.[7] (The approximation in the second value of ε₀ above stems from π being an irrational number.)

Relative permittivity

The linear permittivity of a homogeneous material is usually given relative to that of free space, as a relative permittivity ε_r (also called dielectric constant, although this term is deprecated and sometimes only refers to the static, zero-frequency relative permittivity). In an anisotropic material, the relative permittivity may be a tensor, causing birefringence. The actual permittivity is then calculated by multiplying the relative permittivity by ε₀:

${\displaystyle \varepsilon =\varepsilon _{\mathrm {r} }\varepsilon _{0}=(1+\chi )\varepsilon _{0},}$

where χ (frequently written χ_e) is the electric susceptibility of the material.

The susceptibility is defined as the constant of proportionality (which may be a tensor) relating an electric field E to the induced dielectric polarization density P such that

${\displaystyle \mathbf {P} =\varepsilon _{0}\chi \mathbf {E} ,}$

where ε₀ is the electric permittivity of free space.

The susceptibility of a medium is related to its relative permittivity ε_r by

${\displaystyle \chi =\varepsilon _{\mathrm {r} }-1.}$

So in the case of a vacuum,

${\displaystyle \chi =0.}$

The susceptibility is also related to the polarizability of individual particles in the medium by the Clausius-Mossotti relation.

The electric displacement D is related to the polarization density P by

${\displaystyle \mathbf {D} =\varepsilon _{0}\mathbf {E} +\mathbf {P} =\varepsilon _{0}(1+\chi )\mathbf {E} =\varepsilon _{\mathrm {r} }\varepsilon _{0}\mathbf {E} .}$

The permittivity ε and permeability µ of a medium together determine the phase velocity v = c/n of electromagnetic radiation through that medium:

${\displaystyle \varepsilon \mu ={\frac {1}{v^{2}}}.}$

Practical applications

Determining capacitance

The capacitance of a capacitor is based on its design and architecture, meaning it will not change with charging and discharging. The formula for capacitance in a parallel plate capacitor is written as

${\displaystyle C=\varepsilon \ {\frac {A}{d}}}$

where ${\displaystyle A}$ is the area of one plate, ${\displaystyle d}$ is the distance between the plates, and ${\displaystyle \varepsilon }$ is the permittivity of the medium between the two plates. For a capacitor with relative permittivity ${\displaystyle \kappa }$, it can be said that

${\displaystyle C=\kappa \ \varepsilon _{0}{\frac {A}{d}}}$

Gauss's law

Permittivity is connected to electric flux (and by extension electric field) through Gauss's law. Gauss's law states that for a closed Gaussian surface, s

${\displaystyle \Phi _{E}={\frac {Q_{\text{enc}}}{\varepsilon _{0}}}=\oint _{s}\mathbf {E} \cdot \mathrm {d} \mathbf {A} }$

where ${\displaystyle \Phi _{E}}$ is the net electric flux passing through the surface, ${\displaystyle Q_{\text{enc}}}$ is the charge enclosed in the Gaussian surface, ${\displaystyle \mathbf {E} }$ is the electric field vector at a given point on the surface, and ${\displaystyle \mathrm {d} \mathbf {A} }$ is a differential area vector on the Gaussian surface.

If the Gaussian surface uniformly encloses an insulated, symmetrical charge arrangement, the formula can be simplified to

${\displaystyle EA\cos(\theta )={\frac {Q_{\text{enc}}}{\varepsilon _{0}}}}$

where ${\displaystyle \theta }$ represents the angle between the electric field vector and the area vector.

If all of the electric field lines cross the surface at 90°, the formula can be further simplified to

${\displaystyle E={\frac {Q_{\text{enc}}}{\varepsilon _{0}A}}}$

Because the surface area of a sphere is ${\displaystyle 4\pi r^{2}}$, the electric field a distance ${\displaystyle r}$ away from a uniform, spherical charge arrangement is

${\displaystyle E={\frac {Q}{\varepsilon _{0}A}}={\frac {Q}{\varepsilon _{0}\left(4\pi r^{2}\right)}}}$

${\displaystyle E={\frac {Q}{4\pi \varepsilon _{0}r^{2}}}={\frac {kQ}{r^{2}}}}$

where ${\displaystyle k}$ is Coulomb's constant (${\displaystyle \sim 9.0\times 10^{9}\ {\text{m}}/{\text{F}}}$). This formula applies to the electric field due to a point charge, outside of a conducting sphere or shell, outside of a uniformly charged insulating sphere, or between the plates of a spherical capacitor.

Dispersion and causality

In general, a material cannot polarize instantaneously in response to an applied field, and so the more general formulation as a function of time is

${\displaystyle \mathbf {P} (t)=\varepsilon _{0}\int _{-\infty }^{t}\chi \left(t-t'\right)\mathbf {E} \left(t'\right)\,dt'.}$

That is, the polarization is a convolution of the electric field at previous times with time-dependent susceptibility given by χ(Δt). The upper limit of this integral can be extended to infinity as well if one defines χ(Δt) = 0 for Δt < 0. An instantaneous response would correspond to a Dirac delta function susceptibility χ(Δt) = χδ(Δt).

It is convenient to take the Fourier transform with respect to time and write this relationship as a function of frequency. Because of the convolution theorem, the integral becomes a simple product,

${\displaystyle \mathbf {P} (\omega )=\varepsilon _{0}\chi (\omega )\mathbf {E} (\omega ).}$

This frequency dependence of the susceptibility leads to frequency dependence of the permittivity. The shape of the susceptibility with respect to frequency characterizes the dispersion properties of the material.

Moreover, the fact that the polarization can only depend on the electric field at previous times (i.e. effectively χ(Δt) = 0 for Δt < 0), a consequence of causality, imposes Kramers–Kronig constraints on the susceptibility χ(0).

Complex permittivity

A dielectric permittivity spectrum over a wide range of frequencies. ε′ and ε″ denote the real and the imaginary part of the permittivity, respectively. Various processes are labeled on the image: ionic and dipolar relaxation, and atomic and electronic resonances at higher energies.[8]

As opposed to the response of a vacuum, the response of normal materials to external fields generally depends on the frequency of the field. This frequency dependence reflects the fact that a material's polarization does not change instantaneously when an electric field is applied. The response must always be causal (arising after the applied field), which can be represented by a phase difference. For this reason, permittivity is often treated as a complex function of the (angular) frequency ω of the applied field:

${\displaystyle \varepsilon \rightarrow {\hat {\varepsilon }}(\omega )}$

(since complex numbers allow specification of magnitude and phase). The definition of permittivity therefore becomes

${\displaystyle D_{0}e^{-i\omega t}={\hat {\varepsilon }}(\omega )E_{0}e^{-i\omega t},}$

where

• D₀ and E₀ are the amplitudes of the displacement and electric fields, respectively,
• i is the imaginary unit, i² = −1.

The response of a medium to static electric fields is described by the low-frequency limit of permittivity, also called the static permittivity ε_s (also ε_DC):

${\displaystyle \varepsilon _{\mathrm {s} }=\lim _{\omega \rightarrow 0}{\hat {\varepsilon }}(\omega ).}$

At the high-frequency limit (meaning optical frequencies), the complex permittivity is commonly referred to as ε_∞ (or sometimes ε_opt[9]). At the plasma frequency and below, dielectrics behave as ideal metals, with electron gas behavior. The static permittivity is a good approximation for alternating fields of low frequencies, and as the frequency increases a measurable phase difference δ emerges between D and E. The frequency at which the phase shift becomes noticeable depends on temperature and the details of the medium. For moderate field strength (E₀), D and E remain proportional, and

${\displaystyle {\hat {\varepsilon }}={\frac {D_{0}}{E_{0}}}=|\varepsilon |e^{-i\delta }.}$

Since the response of materials to alternating fields is characterized by a complex permittivity, it is natural to separate its real and imaginary parts, which is done by convention in the following way:

${\displaystyle {\hat {\varepsilon }}(\omega )=\varepsilon '(\omega )-i\varepsilon ''(\omega )=\left|{\frac {D_{0}}{E_{0}}}\right|\left(\cos \delta -i\sin \delta \right).}$

where

• ε′ is the real part of the permittivity;
• ε″ is the imaginary part of the permittivity;
• δ is the loss angle.

The choice of sign for time-dependence, e^−iωt, dictates the sign convention for the imaginary part of permittivity. The signs used here correspond to those commonly used in physics, whereas for the engineering convention one should reverse all imaginary quantities.

The complex permittivity is usually a complicated function of frequency ω, since it is a superimposed description of dispersion phenomena occurring at multiple frequencies. The dielectric function ε(ω) must have poles only for frequencies with positive imaginary parts, and therefore satisfies the Kramers–Kronig relations. However, in the narrow frequency ranges that are often studied in practice, the permittivity can be approximated as frequency-independent or by model functions.

At a given frequency, the imaginary part, ε″, leads to absorption loss if it is positive (in the above sign convention) and gain if it is negative. More generally, the imaginary parts of the eigenvalues of the anisotropic dielectric tensor should be considered.

In the case of solids, the complex dielectric function is intimately connected to band structure. The primary quantity that characterizes the electronic structure of any crystalline material is the probability of photon absorption, which is directly related to the imaginary part of the optical dielectric function ε(ω). The optical dielectric function is given by the fundamental expression:[10]

${\displaystyle \varepsilon (\omega )=1+{\frac {8\pi ^{2}e^{2}}{m^{2}}}\sum _{c,v}\int W_{c,v}(E){\bigl (}\varphi (\hbar \omega -E)-\varphi (\hbar \omega +E){\bigr )}\,dx.}$

In this expression, W_c,v(E) represents the product of the Brillouin zone-averaged transition probability at the energy E with the joint density of states,[11][12] J_c,v(E); φ is a broadening function, representing the role of scattering in smearing out the energy levels.[13] In general, the broadening is intermediate between Lorentzian and Gaussian;[14][15] for an alloy it is somewhat closer to Gaussian because of strong scattering from statistical fluctuations in the local composition on a nanometer scale.

Tensorial permittivity

According to the Drude model of magnetized plasma, a more general expression which takes into account the interaction of the carriers with an alternating electric field at millimeter and microwave frequencies in an axially magnetized semiconductor requires the expression of the permittivity as a non-diagonal tensor.[16] (see also Electro-gyration).

${\displaystyle \mathbf {D} (\omega )={\begin{vmatrix}\varepsilon _{1}&-i\varepsilon _{2}&0\\i\varepsilon _{2}&\varepsilon _{1}&0\\0&0&\varepsilon _{z}\\\end{vmatrix}}\operatorname {\mathbf {E} } (\omega )}$

If ε₂ vanishes, then the tensor is diagonal but not proportional to the identity and the medium is said to be a uniaxial medium, which has similar properties to a uniaxial crystal.

Classification of materials

Classification of materials based on permittivity

εr″/εr′	Current conduction	Field propagation
0		perfect dielectric lossless medium
≪ 1	low-conductivity material poor conductor	low-loss medium good dielectric
≈ 1	lossy conducting material	lossy propagation medium
≫ 1	high-conductivity material good conductor	high-loss medium poor dielectric
∞	perfect conductor

Materials can be classified according to their complex-valued permittivity ε, upon comparison of its real ε′ and imaginary ε″ components (or, equivalently, conductivity, σ, when accounted for in the latter). A perfect conductor has infinite conductivity, σ = ∞, while a perfect dielectric is a material that has no conductivity at all, σ = 0; this latter case, of real-valued permittivity (or complex-valued permittivity with zero imaginary component) is also associated with the name lossless media.[17] Generally, when σ/ωε′ ≪ 1 we consider the material to be a low-loss dielectric (although not exactly lossless), whereas σ/ωε′ ≫ 1 is associated with a good conductor; such materials with non-negligible conductivity yield a large amount of loss that inhibit the propagation of electromagnetic waves, thus are also said to be lossy media. Those materials that do not fall under either limit are considered to be general media.

Lossy medium

In the case of a lossy medium, i.e. when the conduction current is not negligible, the total current density flowing is:

${\displaystyle J_{\text{tot}}=J_{\mathrm {c} }+J_{\mathrm {d} }=\sigma E+i\omega \varepsilon 'E=i\omega {\hat {\varepsilon }}E}$

where

• σ is the conductivity of the medium;
• ${\displaystyle \varepsilon '=\varepsilon _{0}\varepsilon _{r}}$ is the real part of the permittivity.
• ${\displaystyle {\hat {\varepsilon }}=\varepsilon '-i\varepsilon ''}$is the complex permittivity

The size of the displacement current is dependent on the frequency ω of the applied field E; there is no displacement current in a constant field.

In this formalism, the complex permittivity is defined as:[18][19]

${\displaystyle {\hat {\varepsilon }}=\varepsilon '\left(1-i{\frac {\sigma }{\omega \varepsilon '}}\right)=\varepsilon '-i{\frac {\sigma }{\omega }}}$

In general, the absorption of electromagnetic energy by dielectrics is covered by a few different mechanisms that influence the shape of the permittivity as a function of frequency:

• First are the relaxation effects associated with permanent and induced molecular dipoles. At low frequencies the field changes slowly enough to allow dipoles to reach equilibrium before the field has measurably changed. For frequencies at which dipole orientations cannot follow the applied field because of the viscosity of the medium, absorption of the field's energy leads to energy dissipation. The mechanism of dipoles relaxing is called dielectric relaxation and for ideal dipoles is described by classic Debye relaxation.
• Second are the resonance effects, which arise from the rotations or vibrations of atoms, ions, or electrons. These processes are observed in the neighborhood of their characteristic absorption frequencies.

The above effects often combine to cause non-linear effects within capacitors. For example, dielectric absorption refers to the inability of a capacitor that has been charged for a long time to completely discharge when briefly discharged. Although an ideal capacitor would remain at zero volts after being discharged, real capacitors will develop a small voltage, a phenomenon that is also called soakage or battery action. For some dielectrics, such as many polymer films, the resulting voltage may be less than 1–2% of the original voltage. However, it can be as much as 15–25% in the case of electrolytic capacitors or supercapacitors.

Quantum-mechanical interpretation

In terms of quantum mechanics, permittivity is explained by atomic and molecular interactions.

At low frequencies, molecules in polar dielectrics are polarized by an applied electric field, which induces periodic rotations. For example, at the microwave frequency, the microwave field causes the periodic rotation of water molecules, sufficient to break hydrogen bonds. The field does work against the bonds and the energy is absorbed by the material as heat. This is why microwave ovens work very well for materials containing water. There are two maxima of the imaginary component (the absorptive index) of water, one at the microwave frequency, and the other at far ultraviolet (UV) frequency. Both of these resonances are at higher frequencies than the operating frequency of microwave ovens.

At moderate frequencies, the energy is too high to cause rotation, yet too low to affect electrons directly, and is absorbed in the form of resonant molecular vibrations. In water, this is where the absorptive index starts to drop sharply, and the minimum of the imaginary permittivity is at the frequency of blue light (optical regime).

At high frequencies (such as UV and above), molecules cannot relax, and the energy is purely absorbed by atoms, exciting electron energy levels. Thus, these frequencies are classified as ionizing radiation.

While carrying out a complete ab initio (that is, first-principles) modelling is now computationally possible, it has not been widely applied yet. Thus, a phenomenological model is accepted as being an adequate method of capturing experimental behaviors. The Debye model and the Lorentz model use a first-order and second-order (respectively) lumped system parameter linear representation (such as an RC and an LRC resonant circuit).

Measurement

The relative permittivity of a material can be found by a variety of static electrical measurements. The complex permittivity is evaluated over a wide range of frequencies by using different variants of dielectric spectroscopy, covering nearly 21 orders of magnitude from 10⁻⁶ to 10¹⁵ hertz. Also, by using cryostats and ovens, the dielectric properties of a medium can be characterized over an array of temperatures. In order to study systems for such diverse excitation fields, a number of measurement setups are used, each adequate for a special frequency range.

Various microwave measurement techniques are outlined in Chen et al..[20] Typical errors for the Hakki-Coleman method employing a puck of material between conducting planes are about 0.3%.[21]

• Low-frequency time domain measurements (10⁻⁶ to 10³ Hz)
• Low-frequency frequency domain measurements (10⁻⁵ to 10⁶ Hz)
• Reflective coaxial methods (10⁶ to 10¹⁰ Hz)
• Transmission coaxial method (10⁸ to 10¹¹ Hz)
• Quasi-optical methods (10⁹ to 10¹⁰ Hz)
• Terahertz time-domain spectroscopy (10¹¹ to 10¹³ Hz)
• Fourier-transform methods (10¹¹ to 10¹⁵ Hz)

At infrared and optical frequencies, a common technique is ellipsometry. Dual polarisation interferometry is also used to measure the complex refractive index for very thin films at optical frequencies.

See also

• Acoustic attenuation
• Density functional theory
• Electric-field screening
• Green–Kubo relations
• Green's function (many-body theory)
• Linear response function
• Rotational Brownian motion
• Electromagnetic permeability

Notes

1. Current practice of standards organizations such as NIST and BIPM is to use c₀, rather than c, to denote the speed of light in a vacuum according to ISO 31. In the original Recommendation of 1983, the symbol c was used for this purpose.[6]

References

1. Landau, L. D.; Lifshitz, E. M.; Pitaevskii, L. P. (2009). Electrodynamics of continuous media. Elsevier Butterworth-Heinemann. ISBN 978-0-7506-2634-7. OCLC 756385298.
2. IEEE Standards Board (1997). "IEEE Standard Definitions of Terms for Radio Wave Propagation". p. 6.
3. Braslavsky, S.E. (2007). "Glossary of terms used in photochemistry (IUPAC recommendations 2006)" (PDF). Pure and Applied Chemistry. 79 (3): 293–465. doi:10.1351/pac200779030293. S2CID 96601716.
4. International Bureau of Weights and Measures (2006), The International System of Units (SI) (PDF) (8th ed.), ISBN 92-822-2213-6, archived (PDF) from the original on 2017-08-14, p. 119
5. "2018 CODATA Value: vacuum electric permittivity". The NIST Reference on Constants, Units, and Uncertainty. NIST. 20 May 2019. Retrieved 2019-05-20.
6. NIST Special Publication 330, Appendix 2, p. 45 .
7. "Latest (2006) values of the constants (NIST)". Physics.nist.gov. 2017-07-01. Retrieved 2018-11-20.
8. "Dielectric Spectroscopy". Archived from the original on 2006-01-18. Retrieved 2018-11-20.
9. Hofmann, Philip (2015-05-26). Solid State Physics (2 ed.). Wiley-VCH. p. 194. ISBN 978-3527412822.
10. Peter Y. Yu; Manuel Cardona (2001). Fundamentals of Semiconductors: Physics and Materials Properties. Berlin: Springer. p. 261. ISBN 978-3-540-25470-6.
11. José García Solé, Jose Solé, Luisa Bausa (2001). An introduction to the optical spectroscopy of inorganic solids. Wiley. Appendix A1, pp, 263. ISBN 978-0-470-86885-0.CS1 maint: multiple names: authors list (link)
12. John H. Moore; Nicholas D. Spencer (2001). Encyclopedia of chemical physics and physical chemistry. Taylor and Francis. p. 105. ISBN 978-0-7503-0798-7.
13. Solé, José García; Bausá, Louisa E; Jaque, Daniel (2005-03-22). Solé and Bausa. p. 10. ISBN 978-3-540-25470-6.
14. Hartmut Haug; Stephan W. Koch (1994). Quantum Theory of the Optical and Electronic Properties of Semiconductors. World Scientific. p. 196. ISBN 978-981-02-1864-5.
15. Manijeh Razeghi (2006). Fundamentals of Solid State Engineering. Birkhauser. p. 383. ISBN 978-0-387-28152-0.
16. Prati E. (2003) "Propagation in gyroelectromagnetic guiding systems", J. of Electr. Wav. and Appl. 17, 8, 1177
17. http://www.ece.rutgers.edu/~orfanidi/ewa/ch01.pdf
18. John S. Seybold (2005), Introduction to RF Propagation, John Wiley & Sons, p. 22, eq. (2.6).
19. Kenneth L. Kaiser, Electromagnetic Shielding, CRC Press, 2005, p. 1·28, eqs. (1.80) and (1.81).
20. Linfeng Chen; V. V. Varadan; C. K. Ong; Chye Poh Neo (2004). "Microwave theory and techniques for materials characterization". Microwave electronics. Wiley. p. 37. ISBN 978-0-470-84492-2.
21. Mailadil T. Sebastian (2008). Dielectric Materials for Wireless Communication. Elsevier. p. 19. ISBN 978-0-08-045330-9.

Further reading

• C. J. F. Bottcher, O. C. von Belle & Paul Bordewijk (1973) Theory of Electric Polarization: Dielectric Polarization, volume 1, (1978) volume 2, Elsevier ISBN 0-444-41579-3.
• Arthur R. von Hippel (1954) Dielectrics and Waves ISBN 0-89006-803-8
• Arthur von Hippel editor (1966) Dielectric Materials and Applications: papers by 22 contributors ISBN 0-89006-805-4.

External links

• Electromagnetism, a chapter from an online textbook
• What's all this trapped charge stuff . . ., A different approach to some capacitor problems
• Complex Permittivity and Refractive Index for Metals
• DrudeLorentz.com Online plotting and parameterisation database of the Drude-Lorentz permittivity models of common metals