  We study integrability properties of a reversible deterministic cellular
automaton (the rule 54 of [Bobenko et al., Commun. Math. Phys. 158, 127
(1993)]) and present a bulk algebraic relation and its inhomogeneous extension
which allow for an explicit construction of Liouvillian decay modes for two
distinct families of stochastic boundary driving. The spectrum of the many-body
stochastic matrix defining the time propagation is found to separate into sets,
which we call orbitals, and the eigenvalues in each orbital are found to obey a
distinct set of Bethe-like equations. We construct the decay modes in the first
orbital (containing the leading decay mode) in terms of an exact inhomogeneous
matrix product ansatz, study the thermodynamic properties of the spectrum and
the scaling of its gap, and provide a conjecture for the Bethe-like equations
for all the orbitals and their degeneracy.
