You could discover where the boundary conditions come in (and they don't have to be at t=0) by solving the same DEs the usual way :)
What you are doing is exploiting that you know the solutions to the DEs with the appropriate boundary conditions already because you know the relationship between PD and dQ/dt.
