Neumann heat content asymptotics with singular initial temperature and singular specific heat

We study the asymptotic behavior of the heat content on a compact Riemannian manifold with boundary and with singular specific heat and singular initial temperature distributions imposing Robin boundary conditions. Assuming the existence of a complete asymptotic series we determine the first three terms in that series. In addition to the general setting, the interval is studied in detail as are recursion relations among the coefficients and the relationship between the Dirichlet and Robin settings.