Rayleigh–Taylor (R–T) instability between plasma species is examined in a kinetic test and near-inertial confinement fusion (ICF) regimes. A transport approximation to the plasma species kinetics is used to represent viscosity and species mass transport within a hydrodynamic fluid code (xRage). R–T simulation results are compared in a kinetic test regime with a fully kinetic particle-in-cell approach [vectorized particle-in-cell (VPIC)] and with an analytic model for the growth rate of R–T instability. Single-mode growth rates from both codes and the analytic model are in reasonable agreement over a range of initial wavelengths including the wavenumber of maximum growth rate. Both codes exhibit similar diffusive mixing fronts. Small code-to-code differences arise from the kinetics, while simulation-analytic model differences arise from several sources dominated by the choice of gradients establishing the hydrostatic equilibrium initial conditions. After demonstrating code agreement in the kinetic test regime, which is practically accessible to the VPIC code, then the xRage code, with the fluid plasma transport approximation, is applied to single mode R–T instability under deceleration conditions closer to an ICF implosion, approximated with a carbon (C) shell imploding on a deuterium (D) fuel. The analytic wavelength of maximum instability is limited by the kinetics, primarily in the viscosity, and is found to be ≈10 μm for an ion temperature near 1 keV at this C–D interface, with the most unstable wavelength increasing as temperature increases. The analytic viscous model agrees with simulation results over a range of initial perturbation wavelengths, provided the simulation results are analyzed over a sufficiently short duration (⪅0.2 ns in this case). Details of the fluid structure evolution during this R–T deceleration are compared between the inviscid Euler equations and cases, which include plasma transport over a range in initial wavelengths and initial perturbation amplitudes. The inviscid Euler solutions show a grid-dependent cascade to smaller scale structures often seen in the R–T instability, while simulations with plasma transport in this deceleration regime develop a single vortex roll-up, as the plasma transport smoothes all hydrodynamic fluid structures smaller than several micrometers. This leads to a grid-converged transient solution for the R–T instability when kinetic effects are included in the simulations, and thus represents a direct numerical simulation of the thermal ions during R–T unstable mixing in ICF relevant conditions.