Power dissipation and entropy production rate of high-dimensional optical matter systems

Entropy production is an essential aspect of creating and maintaining nonequilibrium systems. Despite their ubiquity, calculation of entropy production rates is challenging for high-dimensional systems, so it has only been reported for simple (i.e., l-particle) systems. Moreover, there is a dearth of nontrivial experimental systems where precise measurements of entropy production rate and characterization of the nonequilibrium steady state (NESS) are simultaneously possible. We report an approach to calculate the entropy production rate of overdamped, nonconservative, $N$-body systems and demonstrate this on a six-particle triangle optical matter (OM) system as a nontrivial example. OM systems consist of (nano-)particles organized into ordered arrays that are bound by electrodynamic interactions associated with the scattering and interference of light, and the associated induced-polarizations in and among the particles in coherent optical beams. The flux of laser light in OM systems in a solution environment necessitates that they dissipate energy, produce entropy, and relax to a NESS. The NESS may have several ordered particle configurations (i.e., isomers) that can interchange by barrier crossing processes. Understanding the power dissipation and entropy production rate of a NESS in an OM system along different (collective) modes of motion can advance understanding of the relative stability of the NESSs as well as inform design and control of OM structures. Therefore, we compute the components of the entropy production rate and power dissipation along the collective coordinates of the 6 Ag nanoparticle triangle OM system from OM NESS trajectory data and verify the Seifert relation [U. Seifert, Rep. Prog. Phys. 75, 126001 (2012)] for these complex systems with a nuanced interpretation.