Quantum Electrodynamics Vacuum Fluctuation Casimir Force MEMS

# Quantum Electrodynamic Vacuum Fluctuation Forces and Casimir-Polder Phenomena in Nanoscale Systems

## 1. Introduction: Quantum Vacuum and Zero-Point Energy

The quantum vacuum is not empty—it is seething with virtual particle-antiparticle pairs constantly appearing and annihilating. These zero-point fluctuations of the electromagnetic field carry measurable energy density and exert forces on macroscopic objects.

The zero-point energy of a quantum harmonic oscillator (mode $\omega$ ):

E_{\text{vac}}(\omega) = \frac{1}{2}\hbar\omega

Summing over all electromagnetic modes in a cavity:

E_{\text{total}} = \sum_{\vec{k}} \frac{1}{2}\hbar\omega_{\vec{k}} = \frac{\hbar}{2} \int_0^\infty \frac{\rho(\omega)}{\pi} \omega \, d\omega

where $\rho(\omega)$ is the density of electromagnetic states.

In free space, this energy density is formally infinite, but Casimir forces arise from the difference in zero-point energy between configurations—a finite and measurable quantity.

## 2. Casimir Effect: Parallel Plate Geometry

The paradigmatic Casimir effect occurs between two parallel conducting plates separated by distance $$d$$ . The electromagnetic modes between the plates are quantized by boundary conditions (perfect conductors), while outside the plates, modes are unquantized.

### Mode Quantization

Between plates, the wave vector parallel to the surface is arbitrary, but the perpendicular component (z-direction) must satisfy:

k_z = \frac{n\pi}{d}, \quad n = 1, 2, 3, \ldots

The energy difference (per unit area) between the confined region and free space:

\Delta E(d) = E_{\text{inside}} - E_{\text{outside}} = -\frac{\pi^2 \hbar c}{720 d^3}

The Casimir force per unit area (negative = attractive):

F_C = -\frac{\partial E}{\partial d} = -\frac{\pi^2 \hbar c}{240 d^4}

This force is purely relativistic and quantum mechanical—it vanishes in the classical limit ( $\hbar \to 0$ ) and requires relativistic treatment of the electromagnetic field.

### Experimental Verification

The Casimir force was experimentally confirmed by:
- Lamoreaux (1997): Precision measurement with conducting plate and Au sphere, achieved $\Delta F / F \sim 1$ % precision
- Harris et al. (2000): Direct force measurement between polished plates
- MEMS devices (2001-present): Detection of Casimir stiction in micro-scale structures

## 3. Casimir-Polder Interaction: Atom-Surface Potential

For an atom near a dielectric or conducting surface, the interaction arises from polarization of the atom by vacuum electromagnetic fluctuations. The Casimir-Polder potential at large distance $$r$$ :

V_{\text{CP}}(r) = -\frac{\hbar c}{2\pi} \alpha(i\omega) \int_0^\infty d\xi \, G(\xi, r)

where $\alpha(i\omega)$ is the dynamic polarizability and $G(\xi, r)$ is the Green's function for electromagnetic propagation.

For atom above conducting plate (London dispersion limit):

V(r) = -\frac{\hbar c \alpha_0}{16\pi^2 \epsilon_0} \left[\frac{1}{r^4} + \frac{1}{(2r)^4} + \ldots\right] \approx -\frac{C_4}{r^4}

For atom above dielectric interface:

V(r) = -\frac{\hbar c}{32\pi^2 \epsilon_0} \frac{\epsilon(i\omega) - 1}{\epsilon(i\omega) + 1} \frac{\alpha(i\omega)}{r^4}

The $r^{-4}$ power law dominates the potential at nanoscale separations.

## 4. Lifshitz Theory: Dispersion Forces Between Dielectric Media

Lifshitz theory generalizes the Casimir effect to arbitrary dielectric media separated by a finite gap (vacuum or medium). The force per unit area between two semi-infinite dielectrics:

P(d) = \frac{\hbar c}{16\pi^2 d^4} \int_0^\infty d\xi \left[\text{Tr}(\ln[\mathbb{I} - \mathcal{D}])\right]

where $\mathcal{D}$ is the dynamical matrix relating electromagnetic amplitudes in the two media.

For identical media ( $\epsilon_1 = \epsilon_2 = \epsilon$ ) separated by vacuum:

P(d) = -\frac{\hbar c}{8\pi d^4} \int_0^\infty d\xi \, \frac{(\epsilon(i\xi) - 1)^2}{[\epsilon(i\xi) + 1]^2}

The integral is dominated by virtual UV photons ( $\xi \sim \hbar \omega_{\text{UV}} / c$ ), making the force sensitive to dielectric response across the entire spectrum.

## 5. Casimir Force Between Multi-Layer Stacks

For a finite stack of dielectric layers (relevant to photoresist stacks, EUV pellicles), the Lifshitz formula becomes a product over reflection coefficients:

\mathcal{D} = \prod_{i} r_i \exp(2ik_i d_i) r_{i+1} \exp(2ik_{i+1} d_{i+1})

Each layer's contribution is phase-weighted by the thickness $$d_i$$ and wave vector $k_i(\xi)$ in the imaginary frequency domain.

Resonances occur when $2k_i d_i \approx n\pi$ —these enhance or suppress the Casimir force depending on the refractive index profile.

## 6. Retardation Effects and Frequency Dependence

At short separations ( $d \ll \lambda_{\text{plasma}}$ ), the non-retarded (London-van der Waals) limit applies:

V(r) \propto -\frac{1}{r^6}

At large separations ( $d \gg \lambda_{\text{plasma}} = 2\pi c/\omega_p$ ), the retarded limit dominates:

V(r) \propto -\frac{1}{r^4}

The transition occurs around $d \sim 100$ nm for typical materials, making retardation essential for NEMS/MEMS at micrometer scales.

## 7. Temperature Effects on Casimir Force

At finite temperature, thermal photons contribute to the vacuum energy. The temperature correction:

F(T) = F(0) + \frac{k_B T}{d} \sum_{n=1}^\infty' \frac{\zeta(2n)}{(\pi d / \hbar c)^{2n}}

where $\zeta(2n)$ is the Riemann zeta function and the prime denotes exclusion of $$n=0$$ .

At room temperature ( $$T = 300$$ K) and $$d = 100$$ nm:
- Classical term: $\sim 1 \times 10^{-3}$ N/m²
- Quantum vacuum term: $\sim 1 \times 10^{-2}$ N/m² (dominates)

## 8. Stiction and Adhesion in Photoresist Structures

During extreme ultraviolet (EUV) lithography, narrow photoresist lines (width < 20 nm) with large aspect ratios experience van der Waals and Casimir attraction causing:

1. Line collapse: Sidewalls adhere due to capillary forces and Casimir attraction
2. Pattern roughness: Randomly collapsed regions create LWR (line width roughness)
3. Resist scumming: Incompletely developed resist blocks transfer to substrate

For a resist pillar of height $$h$$ , width $$w$$ , and separation $$s$$ from adjacent pillar:

F_{\text{adhesion}} = \frac{C_3 w h}{s^3} + \frac{C_4 w h}{s^4}

where $$C_3$$ dominates at separations >100 nm (van der Waals) and $$C_4$$ at <100 nm (Casimir).

Mitigation strategies:
- Resist chemistry: Lower surface energy polymers reduce van der Waals
- Topcoat layers: Low- $$n$$ polymer overlayers decouple resist-vacuum interaction
- Drying aids: Supercritical $\text{CO}_2$ extraction prevents capillary collapse

## 9. NEMS Stiction and Electrostatic Levitation

In Nano-Electro-Mechanical Systems (NEMS), suspended beams and cantilevers are subject to Casimir-induced stiction when operated at low power:

F_{\text{Casimir}} = -\frac{A \pi^2 \hbar c}{240 d^4}

where $$A$$ is the contact area.

For a 1 μm × 100 nm beam at $$d = 100$$ nm separation:
- Casimir force: $\sim 100$ nN
- Required electrostatic restoring force: $$F_E = C V^2 / d$$ (capacitive)

The pull-in voltage where Casimir attraction overcomes electrostatic repulsion:

V_{\text{pull-in}} = \sqrt{\frac{240 k d^5}{\pi^2 \hbar c C}}

For NEMS actuators operating in vacuum, Casimir forces set a fundamental minimum operating gap—typically 50-200 nm.

## 10. Python Implementation: Casimir Force Calculator

import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import quad, simps
from scipy.optimize import fminbound

def drude_model_permittivity(omega, omega_p=9.0, gamma=0.04):
    """
    Drude model for metal permittivity.
    
    ε(ω) = 1 - ω_p² / (ω² + iγω)
    
    Parameters:
    -----------
    omega : array (eV)
        Angular frequency
    omega_p : float (eV)
        Plasma frequency
    gamma : float (eV)
        Damping coefficient
    
    Returns:
    --------
    epsilon : complex array
        Complex permittivity
    """
    
    omega = np.asarray(omega)
    scalar_input = False
    if omega.ndim == 0:
        omega = omega[None]
        scalar_input = True
    
    epsilon = 1.0 - omega_p**2 / (omega**2 + 1j * gamma * omega)
    
    if scalar_input:
        return epsilon[0]
    return epsilon

def lifshitz_integrand_parallel(xi, d, eps1, eps2, polarization='TM'):
    """
    Compute Lifshitz integrand for parallel plates.
    
    Parameters:
    -----------
    xi : float
        Imaginary frequency variable (ℏ × imaginary freq / eV)
    d : float (m)
        Separation distance
    eps1, eps2 : complex
        Permittivities at imaginary frequency
    polarization : str
        'TM' (transverse magnetic) or 'TE'
    
    Returns:
    --------
    integrand : float
        Contribution to Lifshitz integral
    """
    
    c = 3e8  # m/s
    hbar = 1.055e-34  # J·s
    
    # Wave vectors
    k_perp = np.sqrt(xi**2 / c**2 + 1e-30)  # Avoid singularity
    
    # Reflection coefficients in imaginary frequency domain
    if polarization == 'TM':
        r1 = (eps1 * k_perp - np.sqrt(eps1 * k_perp**2 + xi**2/c**2)) / \
             (eps1 * k_perp + np.sqrt(eps1 * k_perp**2 + xi**2/c**2))
        r2 = (eps2 * k_perp - np.sqrt(eps2 * k_perp**2 + xi**2/c**2)) / \
             (eps2 * k_perp + np.sqrt(eps2 * k_perp**2 + xi**2/c**2))
    else:  # TE
        r1 = (k_perp - np.sqrt(k_perp**2 + xi**2/c**2)) / \
             (k_perp + np.sqrt(k_perp**2 + xi**2/c**2))
        r2 = r1
    
    # Phase factor
    phase = 2.0 * np.sqrt(xi**2/c**2 + eps1*k_perp**2) * d
    
    # Lifshitz integrand
    integrand = np.real(r1 * r2 * np.exp(2*phase)) / (1 - r1*r2*np.exp(2*phase) + 1e-30)
    
    return np.abs(integrand)

def casimir_force_per_area(d, material1='Au', material2='Au', T=300):
    """
    Compute Casimir force per unit area between two materials.
    
    F/A = (ℏc / 8π) ∫₀^∞ dξ Tr[ln(1 - DD†)]
    
    Parameters:
    -----------
    d : float (m)
        Separation distance
    material1, material2 : str
        'Au', 'Ag', 'Cu', 'Vacuum'
    T : float (K)
        Temperature
    
    Returns:
    --------
    F_per_area : float (N/m²)
        Casimir pressure (negative = attractive)
    """
    
    c = 3e8  # m/s
    hbar = 1.055e-34  # J·s
    k_B = 1.381e-23  # J/K
    
    # Plasma frequencies (eV)
    omega_p = {'Au': 9.0, 'Ag': 9.2, 'Cu': 10.8, 'Vacuum': 0.0}
    damping = {'Au': 0.04, 'Ag': 0.03, 'Cu': 0.05, 'Vacuum': 0.0}
    
    xi_array = np.logspace(-3, 3, 500)  # Imaginary frequency range
    dxi = xi_array[1] - xi_array[0]
    
    integrand = np.zeros_like(xi_array)
    
    for i, xi in enumerate(xi_array):
        eps1 = drude_model_permittivity(1j*xi, omega_p.get(material1, 9.0), 
                                       damping.get(material1, 0.04))
        eps2 = drude_model_permittivity(1j*xi, omega_p.get(material2, 9.0), 
                                       damping.get(material2, 0.04))
        
        integrand[i] = lifshitz_integrand_parallel(xi, d, eps1, eps2, 'TM')
    
    # Integrate
    integral = simps(integrand, xi_array)
    
    # Casimir pressure
    F_per_area = -(hbar * c / (8 * np.pi)) * integral / d  # Pressure = F/A
    
    # Temperature correction (Matsubara sum)
    n_max = 20
    T_correction = 0
    for n in range(1, n_max):
        xi_n = 2 * np.pi * k_B * T * n / hbar
        eps1 = drude_model_permittivity(1j*xi_n, omega_p.get(material1, 9.0), 
                                       damping.get(material1, 0.04))
        eps2 = drude_model_permittivity(1j*xi_n, omega_p.get(material2, 9.0), 
                                       damping.get(material2, 0.04))
        T_correction += lifshitz_integrand_parallel(xi_n, d, eps1, eps2, 'TM')
    
    F_per_area += (k_B * T / d) * T_correction
    
    return F_per_area

def stiction_force_photoresist(width, height, gap, rho_resist=1.1e3):
    """
    Compute van der Waals + Casimir adhesion force on photoresist pillar.
    
    Parameters:
    -----------
    width : float (m)
        Pillar width
    height : float (m)
        Pillar height
    gap : float (m)
        Separation from adjacent pillar
    rho_resist : float (kg/m³)
        Resist density
    
    Returns:
    --------
    F_adhesion : float (N)
        Total adhesion force
    """
    
    hbar = 1.055e-34
    c = 3e8
    
    # Van der Waals (Hamaker constant for resist ~10^-19 J)
    A_hamaker = 1e-19  # J
    F_vdw = A_hamaker * width * height / gap**3
    
    # Casimir (C4 coefficient)
    C4 = np.pi**2 * hbar * c / 240
    F_casimir = C4 * width * height / gap**4
    
    return F_vdw + F_casimir

# Calculations
separations = np.logspace(-8, -6, 100)  # 10 nm to 1 μm
F_Au_Au = np.array([casimir_force_per_area(d, 'Au', 'Au') for d in separations])
F_Au_Ag = np.array([casimir_force_per_area(d, 'Au', 'Ag') for d in separations])

# Photoresist stiction
width_resist = 15e-9  # 15 nm
height_resist = 100e-9  # 100 nm
gaps = np.logspace(-8, -7, 50)  # 10-100 nm
F_adhesion = np.array([stiction_force_photoresist(width_resist, height_resist, g) 
                       for g in gaps])

# Plotting
fig, axes = plt.subplots(2, 2, figsize=(13, 10))

# Panel A: Casimir force vs separation
axes[0, 0].loglog(separations*1e9, np.abs(F_Au_Au), 'b-', linewidth=2, label='Au-Au')
axes[0, 0].loglog(separations*1e9, np.abs(F_Au_Ag), 'r--', linewidth=2, label='Au-Ag')
axes[0, 0].axhline(y=1e-5, color='gray', linestyle=':', alpha=0.7, label='Measurement limit')
axes[0, 0].set_xlabel('Separation (nm)', fontsize=11)
axes[0, 0].set_ylabel('Casimir Pressure (Pa)', fontsize=11)
axes[0, 0].set_title('Casimir Force vs. Material Pair', fontsize=12, fontweight='bold')
axes[0, 0].legend(fontsize=10)
axes[0, 0].grid(True, alpha=0.3, which='both')

# Panel B: Force power laws
# Fit to r^-4 and r^-3 in different regimes
short_range_idx = separations < 200e-9
long_range_idx = separations > 200e-9

axes[0, 1].loglog(separations[long_range_idx]*1e9, np.abs(F_Au_Au[long_range_idx]), 
                 'b-', linewidth=2.5, label='Data (retarded)')
fit_r4 = np.abs(F_Au_Au[50]) * (separations[50] / separations)**4
axes[0, 1].loglog(separations*1e9, fit_r4, 'b--', alpha=0.5, linewidth=1.5, label='∝ r⁻⁴')

axes[0, 1].set_xlabel('Separation (nm)', fontsize=11)
axes[0, 1].set_ylabel('Casimir Pressure (Pa)', fontsize=11)
axes[0, 1].set_title('Retarded Regime (r⁻⁴ Scaling)', fontsize=12, fontweight='bold')
axes[0, 1].legend(fontsize=10)
axes[0, 1].grid(True, alpha=0.3, which='both')

# Panel C: Photoresist stiction
axes[1, 0].loglog(gaps*1e9, F_adhesion*1e9, 'g-o', linewidth=2, markersize=5)
axes[1, 0].set_xlabel('Gap (nm)', fontsize=11)
axes[1, 0].set_ylabel('Adhesion Force (nN)', fontsize=11)
axes[1, 0].set_title('Photoresist Pillar Stiction (15×100 nm)', fontsize=12, fontweight='bold')
axes[1, 0].grid(True, alpha=0.3, which='both')

# Panel D: Temperature dependence
temps = np.array([100, 300, 500])
d_temp = 100e-9
F_temps = [casimir_force_per_area(d_temp, 'Au', 'Au', T=T) for T in temps]

axes[1, 1].bar(temps, np.abs(F_temps), width=100, color='orange', alpha=0.7, edgecolor='black')
axes[1, 1].set_xlabel('Temperature (K)', fontsize=11)
axes[1, 1].set_ylabel('Casimir Pressure @ 100 nm (Pa)', fontsize=11)
axes[1, 1].set_title('Temperature Dependence', fontsize=12, fontweight='bold')
axes[1, 1].grid(True, alpha=0.3, axis='y')

plt.tight_layout()
plt.savefig('casimir_force_analysis.png', dpi=150, bbox_inches='tight')
plt.show()

print("=== Casimir Force Analysis Complete ===")
print(f"Au-Au Casimir pressure @ 100 nm: {casimir_force_per_area(100e-9, 'Au', 'Au'):.2e} Pa")
print(f"Photoresist stiction force @ 20 nm: {stiction_force_photoresist(15e-9, 100e-9, 20e-9)*1e9:.2f} nN")

## 11. Multilayer Thin-Film Stacks and Resonance Enhancement

In photoresist-on-substrate stacks, constructive interference can enhance the Casimir force by factors of 2-10. This occurs when layer thicknesses satisfy:

d_{\text{layer}} \approx \frac{n\lambda_{\text{plasma}}}{4}

where $\lambda_{\text{plasma}} = 2\pi c / \omega_p$ .

For polymeric resists ( $\omega_p \sim 10$ eV), this resonance occurs at layer thicknesses ~50-100 nm—exactly the regime of high-resolution patterning.

## 12. EUV Pellicle Dynamics

Pellicles (protective membranes atop EUV masks) experience:
- Electrostatic attraction to mask substrate (tens of nN)
- Casimir attraction at small gaps (<100 nm)
- Thermal expansion stresses during 250 W exposure

The pellicle sag (deflection) under combined loads:

w_{\text{max}} = \frac{P d^4}{384 EI}

where $$P$$ includes Casimir pressure and thermal gradients.

Finite-element simulations show that Casimir forces can cause 20-50 nm sag for 0.5 μm pellicles, reducing optical transmission by ~5%.

## 13. Comparison with Other Adhesion Mechanisms

Mechanism	Range	Strength	Temperature Dependence
Chemical bonding	0.1-1 Å	~eV	Strong (activation)
van der Waals	1-10 nm	pJ scale	Weak (∝ 1/T)
Casimir (retarded)	10-1000 nm	fJ scale	Weak (∝ 1/T)
Capillary (wet)	<100 nm	Variable	Very strong
Electrostatic	>100 nm	Tunable	Very weak

Casimir forces are most relevant in the 10-100 nm regime for dry MEMS/NEMS.

## 14. Experimental Techniques and Measurement Challenges

Modern Casimir force measurements use:
- Atomic Force Microscopy (AFM): Measure deflection of cantilever spring
- Laser interferometry: Track separation distance with pm precision
- Electrostatic levitation: Balance Casimir with opposing electrostatic force

Key challenges:
- Surface roughness increases effective distance ( $\sim 10$ nm typical)
- Patch potentials create spurious electrostatic artifacts
- Thermal drift limits long-term stability

## 15. Implications for Nanofabrication and Device Design

Casimir forces set fundamental limits on:
- Photoresist line collapse in EUV lithography (<22 nm designs)
- NEMS actuator stability (minimum gap ~50 nm)
- Thermal conductance (gap optimization)
- Quantum capacitors (nonlinear gap dependence)

Future mitigation includes surface modification (coatings), dynamic actuation, and material engineering to suppress attractive forces in critical applications.

Quantum Electrodynamics Vacuum Fluctuation Casimir Force MEMS

Take quantum electrodynamics vacuum fluctuation casimir force MEMS further