Scaling limit of stationary coupled Sasamoto-Spohn models

We introduce a family of stationary coupled Sasamoto-Spohn models and show that, in the weakly asymmetric regime, they converge to the energy solution of coupled Burgers equations. Moreover, we show that any system of coupled Burgers equations satisfying the so-called trilinear condition ensuring stationarity can be obtained as the scaling limit of a suitable system of coupled Sasamoto-Spohn models. The core of our proof, which avoids the use of spectral gap estimates, consists in a second order Boltzmann-Gibbs principle for the discrete model.