Diffusion of multiple conserved charges from entropy production

We derive dissipative relativistic hydrodynamic equations in the presence of multiple conserved charges, i.e., baryon number ($B$), electric charge ($Q$), and strangeness ($S$), using the Chapman-Enskog (CE) method within the kinetic theory approach. The relativistic Boltzmann equation is solved within the relaxation-time approximation with a momentum-independent relaxation time in the collision term. We derive both first-order (Navier-Stokes limit) and second-order dissipative hydrodynamic equations. Within the kinetic theory framework, using the Boltzmann's H-theorem, and by demanding that for a dissipative system, the entropy must be produced, we find different transport coefficients at the first-order and second-order gradient expansion of the out-of-equilibrium distribution function around the local equilibrium. Apart from the well-known transport coefficients, the shear ($\eta$) and the bulk ($\zeta$) viscosities , we also find the diffusion matrix elements ($\kappa_{qq^{\prime}}$) for the conserved charges $B$, $Q$ and $S$. The diffusion matrix elements ($\kappa_{qq^{\prime}}$) are important to model the multi-component diffusion dynamics sourced by inhomogeneous baryon stopping in the initial state of heavy-ion collisions. We estimate the temperature ($T$) and chemical potential dependence of diagonal and off-diagonal elements of the diffusion matrix elements for the (2+1) flavor quark-gluon plasma. We further estimate the ratio $\kappa_{qq^{\prime}}T/\eta$ for a wide range of temperature and chemical potentials to show the relative importance of the diffusion matrix elements compared to other transport coefficients.