We develop fluctuational electrodynamics for media with nonlinear optical response. In a perturbative manner, we amend the stochastic Helmholtz equation to describe fluctuations in a nonlinear setting, in agreement with the fluctuation dissipation theorem, and identify the local (Rytov) current fluctuations. We show how the linear response (the solution of the scattering problem) of a collection of objects is found from the individual responses, as measured in isolation. As an example, we compute the Casimir force acting between nonlinear objects which approaches the result for linear optics for large separations, and deviates for small distances.

We develop fluctuational electrodynamics for media with nonlinear optical response in and out of thermal equilibrium. Starting from the stochastic nonlinear Helmholtz equation and using the fluctuation dissipation theorem, we obtain perturbatively a deterministic nonlinear Helmholtz equation for the average field, the physical linear response, as well as the fluctuations and Rytov currents. We show that the effects of nonlinear optics, in or out of thermal equilibrium, can be taken into account with an effective, system-aware dielectric function. We discuss the heat radiation of a planar, nonlinear surface, showing that Kirchhoff's must be applied carefully. We find that the spectral emissivity of a nonlinear nanosphere can in principle be negative, implying the possibility of heat flow reversal for specific frequencies.

We study the classical thermal component of Casimir, or van der Waals, forces between point particles with highly anharmonic dipole Hamiltonians when they are subjected to an external electric field. Using a model for which the individual dipole moments saturate in a strong field (a model that mimics the charges in a neutral, perfectly conducting sphere), we find that the resulting Casimir force depends strongly on the strength of the field, as demonstrated by analytical results. For a certain angle between external field and center to center axis, the fluctuation force can be tuned and suppressed to arbitrarily small values. We compare the forces between these anharmonic particles to those between harmonic ones, and also provide a simple formula for asymptotically large external fields, which we expect to be generally valid for the case of saturating dipole moments.