Frequency-Dependent Electric Dipole Polarizability

• Frequency-dependent electric dipole polarizability is a linear-response function that relates an applied electric field at a given frequency to the induced dipole moment in various systems.
• It is characterized by tensorial representations in anisotropic and hyperfine-resolved environments and is analyzed using spectral, sum-over-states, and fluctuation-dissipation approaches.
• Applications span high-precision atomic clocks, optical trapping, plasmonic nanoshells, and periodic systems, with computational methods like coupled-cluster and QED offering precise evaluations.

Frequency-dependent electric dipole-electric dipole polarizability is the linear-response quantity that relates an applied electric field at angular frequency $\omega$ to the induced electric dipole moment. In stationary systems it is the Fourier-domain form of a causal dipole-response kernel; in tensor notation it is written as $\alpha_{ij}(\omega)$ or, in response-theory form, $$\alpha_{AB}(-\omega;\omega)=-\langle\!\langle \mu_A;\mu_B\rangle\!\rangle_\omega$$. In isotropic settings it reduces to a scalar, whereas in anisotropic, hyperfine-resolved, cavity-coupled, periodic, or temporally modulated systems it is intrinsically tensorial and may depend on additional geometric, quantum-number, or time arguments. Recent treatments span high-precision atomic clocks, empirically assembled oscillator-strength distributions, cavity QED response theory, periodic coupled-cluster theory, time-varying particles, Wigner crystals, and multilayer nanoshells (Chakraborty et al., 2023, Babb, 2015, Yuwono et al., 23 Jun 2026, Mirmoosa et al., 2020).

1. Definition as a linear-response function

In the stationary time domain, the induced dipole moment is expressed as a causal convolution,

$$p(t)=\int_0^\infty \alpha(\gamma)\,E(t-\gamma)\,d\gamma,$$

while the frequency-domain relation takes the form

$$p(\omega)=\epsilon_0 \alpha_{ee}(\omega)\,E(\omega)$$

for a stationary, linear, isotropic particle (Mirmoosa et al., 2020). In tensor form for molecular or periodic systems, the induced dipole satisfies

$m_i(\omega)=\sum_j \alpha_{ij}(\omega)E_j(\omega),$

and coupled-cluster linear response uses the convention

$$\alpha_{ij}(\omega)=-\langle\!\langle \mu_i;\mu_j\rangle\!\rangle_\omega$$

(Caricato et al., 2 Aug 2025). In QED-HF, the corresponding working definition is

$$\alpha_{AB}(-\omega;\omega)=-\langle\!\langle \mu_A^e;\mu_B^e\rangle\!\rangle_\omega$$

(Yuwono et al., 23 Jun 2026).

Atomic Stark-shift theory makes the same object directly observable. For a monochromatic field, Chakraborty and Sahoo write the leading light shift of a hyperfine level as

$$\Delta E_{\rm light}=-\frac{1}{2}\,\alpha_{F,M_F}(\omega)\,\mathcal{E}_L^2(\omega),$$

with $\alpha_{F,M_F}(\omega)$ reducing to the static polarizability at $\alpha_{ij}(\omega)$0 (Chakraborty et al., 2023). The same response quantity therefore underlies microscopic spectroscopy, optical trapping, clock-state Stark shifts, and macroscopic constitutive relations.

A persistent source of confusion is not the physics but the notation. The literature represented here uses atomic units, SI units, and the dimensionless relative polarizability $\alpha_{ij}(\omega)$1; it also alternates between scalar, Cartesian-tensor, hyperfine-tensor, and two-time formulations. These are different representations of the same linear electric-dipole response rather than different observables.

2. Spectral, sum-over-states, and correlational representations

A standard stationary representation expresses $\alpha_{ij}(\omega)$2 through discrete and continuum oscillator strengths. For magnesium, the real- and imaginary-frequency forms are

$\alpha_{ij}(\omega)$3

$\alpha_{ij}(\omega)$4

with the Thomas-Reiche-Kuhn sum rule

$\alpha_{ij}(\omega)$5

and the static limit

$\alpha_{ij}(\omega)$6

(Babb, 2015). In that construction, enforcing $\alpha_{ij}(\omega)$7 required inclusion of 3s, 2p, 2s, and 1s strength, and the resulting static polarizability was $\alpha_{ij}(\omega)$8 (Babb, 2015).

The absorptive part can also be derived from fluctuation-dissipation theory. For a Wigner crystal with dipole operator $\alpha_{ij}(\omega)$9, the imaginary part is written as

$$\alpha_{AB}(-\omega;\omega)=-\langle\!\langle \mu_A;\mu_B\rangle\!\rangle_\omega$$0

or equivalently through the dipole spectral density (Singh, 9 Sep 2025). In the high-temperature coarse-grained limit the same work obtains $$\alpha_{AB}(-\omega;\omega)=-\langle\!\langle \mu_A;\mu_B\rangle\!\rangle_\omega$$1 for $$\alpha_{AB}(-\omega;\omega)=-\langle\!\langle \mu_A;\mu_B\rangle\!\rangle_\omega$$2, whereas for $$\alpha_{AB}(-\omega;\omega)=-\langle\!\langle \mu_A;\mu_B\rangle\!\rangle_\omega$$3 the absence of resonant modes causes $$\alpha_{AB}(-\omega;\omega)=-\langle\!\langle \mu_A;\mu_B\rangle\!\rangle_\omega$$4 to drop sharply, implying transparency above the Wigner cutoff (Singh, 9 Sep 2025).

The real part follows from causality through a Kramers-Kronig integral,

$$\alpha_{AB}(-\omega;\omega)=-\langle\!\langle \mu_A;\mu_B\rangle\!\rangle_\omega$$5

a relation stated explicitly for the Wigner-crystal treatment and, in the nonstationary case, satisfied by the temporal complex polarizability for each fixed observation time (Singh, 9 Sep 2025, Mirmoosa et al., 2020). In this sense, frequency-dependent dipole polarizability simultaneously encodes dispersion, absorption, and the oscillator-strength budget of the system.

3. Angular-momentum structure, hyperfine resolution, and Stark shifts

For hyperfine-resolved atomic levels, the polarizability decomposes into scalar, axial-vector, and tensor parts. For the $$\alpha_{AB}(-\omega;\omega)=-\langle\!\langle \mu_A;\mu_B\rangle\!\rangle_\omega$$6Cs clock states,

$$\alpha_{AB}(-\omega;\omega)=-\langle\!\langle \mu_A;\mu_B\rangle\!\rangle_\omega$$7

where $$\alpha_{AB}(-\omega;\omega)=-\langle\!\langle \mu_A;\mu_B\rangle\!\rangle_\omega$$8 is the degree of polarization, $$\alpha_{AB}(-\omega;\omega)=-\langle\!\langle \mu_A;\mu_B\rangle\!\rangle_\omega$$9 is the angle between propagation direction and quantization axis, and $$p(t)=\int_0^\infty \alpha(\gamma)\,E(t-\gamma)\,d\gamma,$$0 is the polarization angle (Chakraborty et al., 2023). Chakraborty and Sahoo further separate each component into hyperfine-independent and hyperfine-induced pieces,

$$p(t)=\int_0^\infty \alpha(\gamma)\,E(t-\gamma)\,d\gamma,$$1

where $$p(t)=\int_0^\infty \alpha(\gamma)\,E(t-\gamma)\,d\gamma,$$2 denotes $$p(t)=\int_0^\infty \alpha(\gamma)\,E(t-\gamma)\,d\gamma,$$3 E1 interactions and $$p(t)=\int_0^\infty \alpha(\gamma)\,E(t-\gamma)\,d\gamma,$$4 M1 hyperfine interactions (Chakraborty et al., 2023).

A central selection rule is that the atomic tensor polarizability $$p(t)=\int_0^\infty \alpha(\gamma)\,E(t-\gamma)\,d\gamma,$$5 vanishes for $$p(t)=\int_0^\infty \alpha(\gamma)\,E(t-\gamma)\,d\gamma,$$6. For the $$p(t)=\int_0^\infty \alpha(\gamma)\,E(t-\gamma)\,d\gamma,$$7 ground state of cesium, this means that any tensor shift at the hyperfine level arises solely from hyperfine-induced third order (Chakraborty et al., 2023). The vector component also acquires hyperfine-level dependence through hyperfine angular factors and the M1 operator.

The same framework directly yields clock-state Stark shifts. The differential clock polarizability is

$$p(t)=\int_0^\infty \alpha(\gamma)\,E(t-\gamma)\,d\gamma,$$8

and the frequency shift is

$$p(t)=\int_0^\infty \alpha(\gamma)\,E(t-\gamma)\,d\gamma,$$9

(Chakraborty et al., 2023). For the static scalar differential shift, the Stark-shift coefficient is

$$p(\omega)=\epsilon_0 \alpha_{ee}(\omega)\,E(\omega)$$0

in excellent agreement with the precise measurement $$p(\omega)=\epsilon_0 \alpha_{ee}(\omega)\,E(\omega)$$1 (Chakraborty et al., 2023). The same study gives $$p(\omega)=\epsilon_0 \alpha_{ee}(\omega)\,E(\omega)$$2 a.u., $$p(\omega)=\epsilon_0 \alpha_{ee}(\omega)\,E(\omega)$$3 a.u., and $$p(\omega)=\epsilon_0 \alpha_{ee}(\omega)\,E(\omega)$$4 a.u., making explicit the strong frequency dependence relevant to optical trapping and clock interrogation (Chakraborty et al., 2023).

4. Computational realizations and many-body methodology

High-accuracy polarizabilities now emerge from several distinct computational pipelines. For cesium clock states, the dominant valence contributions were evaluated by combining relativistic coupled-cluster matrix elements with measurements, while core and mixed terms were treated with random phase approximation. The electronic structure used Dirac-Hartree-Fock with the Dirac-Coulomb Hamiltonian, augmented by Breit and low-order QED corrections, and many-body correlation was included through RCCSD and RCCSDT (Chakraborty et al., 2023). The paper decomposes the response into valence, core, core-core, core-valence, and valence-core sectors and shows that continuum “tail” contributions change the static scalar Stark coefficient by about $$p(\omega)=\epsilon_0 \alpha_{ee}(\omega)\,E(\omega)$$5, making them mandatory for precision work (Chakraborty et al., 2023).

For magnesium, the dynamic function was assembled empirically from discrete resonance oscillator strengths and continuum photoionization cross sections. The adopted principal line strength $$p(\omega)=\epsilon_0 \alpha_{ee}(\omega)\,E(\omega)$$6 brought the oscillator-strength distribution into TRK and $$p(\omega)=\epsilon_0 \alpha_{ee}(\omega)\,E(\omega)$$7 consistency; choosing $$p(\omega)=\epsilon_0 \alpha_{ee}(\omega)\,E(\omega)$$8 or $$p(\omega)=\epsilon_0 \alpha_{ee}(\omega)\,E(\omega)$$9 would have increased $m_i(\omega)=\sum_j \alpha_{ij}(\omega)E_j(\omega),$0 to $m_i(\omega)=\sum_j \alpha_{ij}(\omega)E_j(\omega),$1 or $m_i(\omega)=\sum_j \alpha_{ij}(\omega)E_j(\omega),$2 and $m_i(\omega)=\sum_j \alpha_{ij}(\omega)E_j(\omega),$3 to $m_i(\omega)=\sum_j \alpha_{ij}(\omega)E_j(\omega),$4 or $m_i(\omega)=\sum_j \alpha_{ij}(\omega)E_j(\omega),$5, outside recommended ranges unless other inputs were retuned (Babb, 2015). This construction is not a formal many-body response calculation, but it is an internally constrained reconstruction of the same response function.

In cavity QED Hartree-Fock, the response amplitudes obey a coupled electron-photon linear system,

$m_i(\omega)=\sum_j \alpha_{ij}(\omega)E_j(\omega),$6

with electronic TD-HF blocks, photon blocks at $m_i(\omega)=\sum_j \alpha_{ij}(\omega)E_j(\omega),$7, and a light-matter coupling vector $m_i(\omega)=\sum_j \alpha_{ij}(\omega)E_j(\omega),$8 (Yuwono et al., 23 Jun 2026). The dynamic polarizability is then extracted as

$m_i(\omega)=\sum_j \alpha_{ij}(\omega)E_j(\omega),$9

(Yuwono et al., 23 Jun 2026). The implementation was benchmarked against real-time TD-QED-HF and agreed to $$\alpha_{ij}(\omega)=-\langle\!\langle \mu_i;\mu_j\rangle\!\rangle_\omega$$0 across all tested coupling strengths (Yuwono et al., 23 Jun 2026).

For periodic systems, the first implementation at the CCSD level with periodic boundary conditions was reported for 1D chains in CCResPy, an open-source software based on Python and the NumPy library (Caricato et al., 2 Aug 2025). The response uses a symmetric CC linear-response formulation with perturbed amplitudes at $$\alpha_{ij}(\omega)=-\langle\!\langle \mu_i;\mu_j\rangle\!\rangle_\omega$$1, collective orbital-$$\alpha_{ij}(\omega)=-\langle\!\langle \mu_i;\mu_j\rangle\!\rangle_\omega$$2 indices, and DIIS acceleration. The periodic dipole operator in the length gauge introduces a technical difficulty: “missing integer” terms in the $$\alpha_{ij}(\omega)=-\langle\!\langle \mu_i;\mu_j\rangle\!\rangle_\omega$$3 derivative matrix are usually set to zero, and the paper identifies this as the dominant residual source of discrepancy for longitudinal components in LiH chains (Caricato et al., 2 Aug 2025).

5. Nonstationary, cavity-dressed, and periodic generalizations

When time-translation invariance is broken, the polarizability becomes a two-time kernel,

$$\alpha_{ij}(\omega)=-\langle\!\langle \mu_i;\mu_j\rangle\!\rangle_\omega$$4

and the temporal complex polarizability is defined by

$$\alpha_{ij}(\omega)=-\langle\!\langle \mu_i;\mu_j\rangle\!\rangle_\omega$$5

(Mirmoosa et al., 2020). Under periodic modulation,

$$\alpha_{ij}(\omega)=-\langle\!\langle \mu_i;\mu_j\rangle\!\rangle_\omega$$6

so a single drive frequency produces Floquet sidebands at $$\alpha_{ij}(\omega)=-\langle\!\langle \mu_i;\mu_j\rangle\!\rangle_\omega$$7 (Mirmoosa et al., 2020). A key point made explicitly in the time-varying-particle treatment is that the resulting effective permittivity is fundamentally different from the conventional stationary Drude-Lorentz model with time-dependent parameters simply substituted into it (Mirmoosa et al., 2020).

Cavity dressing modifies the response even for off-resonant perturbing frequencies. In QED-HF the photon enters the electronic response equations through a self-consistent feedback source proportional to

$$\alpha_{ij}(\omega)=-\langle\!\langle \mu_i;\mu_j\rangle\!\rangle_\omega$$8

so the cavity effect depends on the projection $$\alpha_{ij}(\omega)=-\langle\!\langle \mu_i;\mu_j\rangle\!\rangle_\omega$$9 (Yuwono et al., 23 Jun 2026). The same analysis identifies special perturbing frequencies at which the percent change in $$\alpha_{AB}(-\omega;\omega)=-\langle\!\langle \mu_A^e;\mu_B^e\rangle\!\rangle_\omega$$0 crosses zero independent of $$\alpha_{AB}(-\omega;\omega)=-\langle\!\langle \mu_A^e;\mu_B^e\rangle\!\rangle_\omega$$1, consistent with the condition

$$\alpha_{AB}(-\omega;\omega)=-\langle\!\langle \mu_A^e;\mu_B^e\rangle\!\rangle_\omega$$2

(Yuwono et al., 23 Jun 2026). Numerically, p-nitroaniline shows enhancements up to approximately $$\alpha_{AB}(-\omega;\omega)=-\langle\!\langle \mu_A^e;\mu_B^e\rangle\!\rangle_\omega$$3, $$\alpha_{AB}(-\omega;\omega)=-\langle\!\langle \mu_A^e;\mu_B^e\rangle\!\rangle_\omega$$4, and $$\alpha_{AB}(-\omega;\omega)=-\langle\!\langle \mu_A^e;\mu_B^e\rangle\!\rangle_\omega$$5 in $$\alpha_{AB}(-\omega;\omega)=-\langle\!\langle \mu_A^e;\mu_B^e\rangle\!\rangle_\omega$$6, $$\alpha_{AB}(-\omega;\omega)=-\langle\!\langle \mu_A^e;\mu_B^e\rangle\!\rangle_\omega$$7, and $$\alpha_{AB}(-\omega;\omega)=-\langle\!\langle \mu_A^e;\mu_B^e\rangle\!\rangle_\omega$$8 near $$\alpha_{AB}(-\omega;\omega)=-\langle\!\langle \mu_A^e;\mu_B^e\rangle\!\rangle_\omega$$9 a.u. at $$\Delta E_{\rm light}=-\frac{1}{2}\,\alpha_{F,M_F}(\omega)\,\mathcal{E}_L^2(\omega),$$0 a.u. (Yuwono et al., 23 Jun 2026).

Periodic boundary conditions produce a different kind of environment dressing. In the 1D LiH calculations, going from an isolated molecule to a periodic chain inverts the relative magnitude of the polarizability tensor elements: molecular $$\Delta E_{\rm light}=-\frac{1}{2}\,\alpha_{F,M_F}(\omega)\,\mathcal{E}_L^2(\omega),$$1, whereas periodic $$\Delta E_{\rm light}=-\frac{1}{2}\,\alpha_{F,M_F}(\omega)\,\mathcal{E}_L^2(\omega),$$2 (Caricato et al., 2 Aug 2025). The paper interprets this as enhanced longitudinal response and reduced transverse response induced by intermolecular coupling in the chain (Caricato et al., 2 Aug 2025). This suggests that in low-dimensional solids the frequency-dependent dipole polarizability is not merely a localized molecular observable repeated in space, but a collective response per unit cell.

6. Nanostructures, collective modes, applications, and limits of validity

For spherical nanoparticles in the long-wavelength approximation, the frequency-dependent dipole polarizability is already complex because the permittivity is complex and dispersive. For a homogeneous sphere,

$$\Delta E_{\rm light}=-\frac{1}{2}\,\alpha_{F,M_F}(\omega)\,\mathcal{E}_L^2(\omega),$$3

and the Fröhlich condition is approximately $$\Delta E_{\rm light}=-\frac{1}{2}\,\alpha_{F,M_F}(\omega)\,\mathcal{E}_L^2(\omega),$$4 (Ugwuoke et al., 2024). For multilayer nanoshells, the normalized polarizability obeys a recursion,

$$\Delta E_{\rm light}=-\frac{1}{2}\,\alpha_{F,M_F}(\omega)\,\mathcal{E}_L^2(\omega),$$5

with $$\Delta E_{\rm light}=-\frac{1}{2}\,\alpha_{F,M_F}(\omega)\,\mathcal{E}_L^2(\omega),$$6 determined by the inner structure’s scaled polarizability (Ugwuoke et al., 2024). The corresponding Fröhlich function,

$$\Delta E_{\rm light}=-\frac{1}{2}\,\alpha_{F,M_F}(\omega)\,\mathcal{E}_L^2(\omega),$$7

locates dipolar localized surface plasmon resonances (Ugwuoke et al., 2024). The recursive formula reproduces usual long-wavelength-approximation results and was compared with Mie-theory and FEM spectra for structures with up to $$\Delta E_{\rm light}=-\frac{1}{2}\,\alpha_{F,M_F}(\omega)\,\mathcal{E}_L^2(\omega),$$8 layers, but the same paper notes the expected limits of the approximation: under-prediction of absorption magnitudes, smaller redshifts, breakdown at larger size parameter, and possible nonlocal or quantum effects for very thin metallic layers (Ugwuoke et al., 2024).

At the atomic and molecular scale, the same response function is the input for dispersion and surface interactions. The empirically constructed magnesium $$\Delta E_{\rm light}=-\frac{1}{2}\,\alpha_{F,M_F}(\omega)\,\mathcal{E}_L^2(\omega),$$9 yields

$\alpha_{F,M_F}(\omega)$0

with a quadrature value $\alpha_{F,M_F}(\omega)$1 for $\alpha_{F,M_F}(\omega)$2 (Babb, 2015). The magnesium study also shows that incomplete inner-shell modeling is most visible in $\alpha_{F,M_F}(\omega)$3 for $\alpha_{F,M_F}(\omega)$4 a.u., affecting $\alpha_{F,M_F}(\omega)$5 and higher-order dispersion more than $\alpha_{F,M_F}(\omega)$6 (Babb, 2015).

The Wigner-crystal formulation extends the same concept to a lattice of localized electrons. There the imaginary part is computed from the dipole-moment correlator in the harmonic approximation and, with phenomenological damping, becomes a thermally weighted continuum of Lorentz-like oscillators up to the Wigner cutoff $\alpha_{F,M_F}(\omega)$7 (Singh, 9 Sep 2025). A standard ellipsometry set-up is proposed as a test of the predicted low-frequency linear rise of $\alpha_{F,M_F}(\omega)$8 and the transparent regime for applied frequencies greater than $\alpha_{F,M_F}(\omega)$9 (Singh, 9 Sep 2025).

Several limitations recur across these otherwise disparate realizations. The cesium calculations show that continuum tails and mixed core-valence terms, though small, are necessary for percent-level accuracy (Chakraborty et al., 2023). The periodic CCSD implementation avoids explicit damping and therefore omits genuinely resonant regions unless a small imaginary shift is introduced (Caricato et al., 2 Aug 2025). The QED-HF implementation likewise uses Hermitian response without linewidths or explicit Kramers-Kronig analysis (Yuwono et al., 23 Jun 2026). The long-wavelength nanoshell recursion neglects retardation and multipoles by construction (Ugwuoke et al., 2024). These are not contradictions among definitions of polarizability; they are regime-specific truncations applied to the same frequency-dependent electric dipole-electric dipole response function.