We prove that Einstein's equations coupled to equations of Israel-Stewart-type, describing the dynamics of a relativistic fluid with bulk viscosity and nonzero baryon charge (without shear viscosity or baryon diffusion) dynamically coupled to gravity, are causal in the full nonlinear regime. We also show that these equations can be written as a first-order symmetric hyperbolic system, implying local existence and uniqueness of solutions to the equations of motion. We use an arbitrary equation of state and do not make any simplifying symmetry or near-equilibrium assumption, requiring only physically natural conditions on the fields. These results pave the way for the inclusion of bulk viscosity effects in simulations of gravitational-wave signals coming from neutron star mergers.

We consider a first order formulation of relativistic fluids with bulk viscosity based on a stress-energy tensor introduced by Lichnerowicz. Choosing a barotropic equation of state, we show that this theory satisfies basic physical requirements and, under the further assumption of vanishing vorticity, that the equations of motion are causal, both in the case of a fixed background and when the equations are coupled to Einstein's equations. Furthermore, Lichnerowicz's proposal does not fit into the general framework of first order theories studied by Hiscock and Lindblom, and hence their instability results do not apply. These conclusions apply to the full-fledged non-linear theory, without any equilibrium or near equilibrium assumptions. Similarities and differences between the approach here explored and other theories of relativistic viscosity, including the Mueller-Israel-Stewart formulation, are addressed. Cosmological models based on the Lichnerowicz stress-energy tensor are studied. Finally, it is also shown how the present model can be generalized to a second order formulation.

We examine the cosmological consequences of an alternative to the standard expression for bulk viscosity, one which was proposed to avoid the propagation of superluminal signals without the necessity of extending the space of variables of the theory. The Friedmann equation is derived for this case, along with an expression for the effective pressure. We find solutions for the evolution of the density of a viscous component, which differs markedly from the case of conventional Eckart theory; our model evolves toward late-time phantom-like behavior with a future singularity. Entropy production is addressed, and some similarities and differences to approaches based on the Mueller-Israel-Stewart theory are discussed.