Thermodynamics of small systems has become an important field of statistical physics. Such systems are driven out of equilibrium by a control, and the question is naturally posed how such a control can be optimized. We show that optimization problems in small system thermodynamics are solved by (deterministic) optimal transport, for which very efficient numerical methods have been developed, and of which there are applications in cosmology, fluid mechanics, logistics, and many other fields. We show, in particular, that minimizing expected heat released or work done during a nonequilibrium transition in finite time is solved by the Burgers equation and mass transport by the Burgers velocity field. Our contribution hence considerably extends the range of solvable optimization problems in small system thermodynamics.

The Rayleigh–Taylor instability of two immiscible fluids in the limit of small Atwood numbers is studied by means of a phase-field description. In this method, the sharp fluid interface is replaced by a thin, yet finite, transition layer where the interfacial forces vary smoothly. This is achieved by introducing an order parameter (the phase-field) continuously varying across the interfacial layers and uniform in the bulk region. The phase-field model obeys a Cahn–Hilliard equation and is two-way coupled to the standard Navier–Stokes equations. Starting from this system of equations we have first performed a linear analysis from which we have analytically rederived the known gravity–capillary dispersion relation in the limit of vanishing mixing energy density and capillary width. We have performed numerical simulations and identified a region of parameters in which the known properties of the linear phase (both stable and unstable) are reproduced in a very accurate way. This has been done both in the case of negligible viscosity and in the case of non-zero viscosity. In the latter situation, only upper and lower bounds for the perturbation growth rate are known. Finally, we have also investigated the weakly nonlinear stage of the perturbation evolution and identified a regime characterized by a constant terminal velocity of bubbles/spikes. The measured value of the terminal velocity is in agreement with available theoretical prediction. The phase-field approach thus appears to be a valuable technique for the dynamical description of the stages where hydrodynamic turbulence and wave-turbulence come into play.

We study the effects of a stable density stratification on the turbulent dynamics of thin fluid layers forced at intermediate scales. By means of a set of high-resolution numerical simulations, performed within the Boussinesq approximation, we investigate how the stratification and confinement affect the mechanisms of kinetic and potential energy transfer. The detailed analysis of the statistics of the energy-dissipation rates and energy-exchange rates and of the spectral fluxes of potential and kinetic energy shows that stratification provides a new channel for the energy transfer towards small scales which reduces the large-scale flux of kinetic energy. We also discuss the role of vortex stretching and enstrophy flux in the transfer of kinetic energy into potential energy at small scales.

Master equations are one of the main avenues to study open quantum systems. When the master equation is of the Lindblad–Gorini–Kossakowski–Sudarshan form, its solution can be “unraveled in quantum trajectories” i.e., represented as an average over the realizations of a Markov process in the Hilbert space of the system. Quantum trajectories of this type are both an element of quantum measurement theory as well as a numerical tool for systems in large Hilbert spaces. We prove that general time-local and trace-preserving master equations also admit an unraveling in terms of a Markov process in the Hilbert space of the system. The crucial ingredient is to weigh averages by a probability pseudo-measure which we call the “influence martingale”. The influence martingale satisfies a 1d stochastic differential equation enslaved to the ones governing the quantum trajectories. We thus extend the existing theory without increasing the computational complexity.