Motivated by results for the propagation of active-passive interfaces of bacterial Serratia marcescens swarms [Nat. Comm., 9, 5373 (2018)] we use a hydrodynamic multiphase model to investigate the propagation of interfaces of active nematics on substrates. We characterize the active nematic phase of the model and discuss its description of the statistical dynamics of the swarms. We show that the velocity correlation functions and the energy spectrum of the active turbulent phase of the model scale with the active length, for a range of activities. In addition, the energy spectrum exhibits two power-law regimes with exponents close to those reported for other models and bacterial swarms. Although the exponent of the rising branch of the spectrum (small wavevector) appears to be independent of the activity, the exponent of the decay changes with activity systematically, albeit slowly. We characterize also the propagation of circular and flat active-passive interfaces. We find that the closing time of the circular passive domain decays quadratically with the activity and that the structure factor of the flat interface is similar to that reported for swarms, with an activity dependent exponent. Finally, the effect of the substrate friction was investigated. We found an activity dependent threshold, above which the turbulent active nematic forms isolated islands that shrink until the system becomes isotropic and below which the active nematic expands, with a well defined propagating interface. The interface may be stopped by a fricion gradient.

The Boltzmann equation with the Bhatnagar-Gross-Krook collision operator is considered for the Bose-Einstein and Fermi-Dirac equilibrium distribution functions. We show that the expansion of the microscopic velocity in terms of Hermite polynomials must be carried until the fourth order to correctly describe the energy equation. The viscosity and thermal coefficients, previously obtained by J.Y. Yang et al 1,2 through the Uehling-Uhlenbeck approach, are also derived here. Thus the construction of a lattice Boltzmann method for the quantum fluid is possible provided that the Bose-Einstein and Fermi-Dirac equilibrium distribution functions are expanded until fourth order in the Hermite polynomials.

We use numerical simulations to investigate the hydrodynamic behavior of the interface between nematic (N) and isotropic (I) phases of a confined active liquid crystal. At low activities, a stable interface with constant shape and velocity is observed separating the two phases. For nematics in homeotropic channels, the velocity of the interface at the NI transition increases from zero (i) linearly with the activity for contractile systems and (ii) quadratically for extensile ones. Interestingly, the nematic phase expands for contractile systems while it contracts for extensile ones, as a result of the active forces at the interface. Since both activity and temperature affect the stability of the nematic, for active nematics in the stable regime the temperature can be tuned to observe static interfaces, providing an operational definition for the coexistence of active nematic and isotropic phases. At higher activities, beyond the stable regime, an interfacial instability is observed for extensile nematics. In this regime defects are nucleated at the interface and move away from it. The dynamics of these defects is regular and persists asymptotically for a finite range of activities. We used an improved hybrid model of finite differences and lattice Boltzmann method with multirelaxation-time collision operator, the accuracy of which allowed us to characterize the dynamics of the distinct interfacial regimes.