A theoretical model to describe the thermochemical nonequilibrium and radiation phenomena in the shock layer over a spacecraft entering the planet Neptune is derived by assembling existing information. This model is applied to predict the flow in a constant-area channel behind a normal shock wave in a mixture representing the atmosphere of the planet. Calculation is made for aerocapture flights into Neptune that have been considered in the past. It is shown that the presence of methane in the planet's atmosphere shortens the time for equilibration and typically increases the radiative heat flux by about 50%. Radiative heat flux values are calculated over a range of densities and velocities expected in an aerocapture flight into Neptune. Nomenclature a = reaction rate constant [Eq. (1)], cm 3 =S E = energy feedback factor K F = forward reaction rate coefficient, cm 3 =S N = preexponential power in rate coefficient [Eq. (1)], dimensionless P = pressure, Pa or torr Q R = radiative heat flux, Kw=cm 2 t = heavy particle translational temperature, K t A = effective temperature of reaction, K t D = reaction temperature [Eq. (1)], K uS = shock velocity, M=S or km=S X = distance from shock wave, cm = density, kg=M 3 Subscript 1 = freestream Superscript = excited state