The coincidence rate in two detectors that register signal and idler photons in parametric down-conversion is calculated using a quantum description of the radiation field that extends the semiclassical theory to the resonance regime. The signal is given by a sum over paths in the joint field plus matter space, each involving a pair of sites (molecules) and six radiation field modes. It may not be factorized into amplitudes representing one-site, three-mode paths. All effects of multiple vacuum modes are accounted for by theses paths without adding a Langevin noise source. The molecular information is given by a time-ordered superoperator Green's function rather than the susceptibility ${\ensuremath{\chi}}^{(2)}$ used in the semiclassical theory.