In recent years a new interesting research avenue has emerged in non-equilibrium statistical physics, namely studies of collective responses in spatially inhomogeneous systems. Whereas substantial progress has been made in understanding the origins and the often universal nature of cooperative behavior in systems far from equilibrium, it is still unclear whether it is possible to control their global collective stochastic dynamics through local manipulations. Therefore, a comprehensive characterization of spatially inhomogeneous non-equilibrium systems is required.

In the first system, we explore a variant of the Katz–Lebowitz–Spohn (KLS) driven lattice gas in two dimensions, where the lattice is split into two regions that are coupled to heat baths with distinct temperatures T > T_c and T_c respectively, where T_c indicates the critical temperature for phase ordering. The geometry was arranged such that the temperature boundaries are oriented perpendicular or parallel to the external particle drive and resulting net current. For perpendicular orientation of the temperature boundaries, in the hotter region, the system behaves like the (totally) asymmetric exclusion processes (TASEP), and experiences particle blockage in front of the interface to the critical region. This blockage is induced by extended particle clusters, growing logarithmically with system size, in the critical region. We observe the density profiles in both high- and low-temperature subsystems to be similar to the well-characterized coexistence and maximal-current phases in (T)ASEP models with open boundary conditions, which are respectively governed by hyperbolic and trigonometric tangent functions. Yet if the lower temperature is set to T_c, we detect marked fluctuation corrections to the mean-field density profiles, e.g., the corresponding critical KLS power-law density decay near the interfaces into the cooler region.

For parallel orientation of the temperature boundaries, we have explored the changes in the dynamical behavior of the hybrid KLS model that are induced by our choice of the hopping rates across the temperature boundaries. If these hopping rates at the interfaces satisfy particle-hole symmetry, the current difference across them generates a vector flow diagram akin to an infinite flat vortex sheet. We have studied the finite-size scaling of the particle density fluctuations in both temperature regions, and observed that it is controlled by the respective temperature values. If the colder subsystem is maintained at the KLS critical temperature, while the hotter subsystem's temperature is set much higher, the interface current greatly suppresses particle exchange between the two regions. As a result of the ensuing effective subsystem decoupling, strong fluctuations persist in the critical region, whence the particle density fluctuations scale with the KLS critical exponents. However, if both temperatures are set well above the critical temperature, the particle density fluctuations scale according to the totally asymmetric exclusion process. We have also measured the entropy production rate in both subsystems; it displays intriguing algebraic decay in the critical region, while it saturates quickly at a small but non-zero level in the hotter region.

The second system is a lattice gas that simulates the spread of COVID-19 epidemics using the paradigmatic stochastic Susceptible-Infectious-Recovered (SIR) model. In our effort to control the spread of the infection of a lattice, we robustly find that the intensity and spatial spread on the epidemic recurrence wave can be limited to a manageable extent provided release of these restrictions is delayed sufficiently (for a duration of at least thrice the time until the peak of the unmitigated outbreak).