Fluctuational electrodynamics for nonlinear media

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.