We report the detection of far-IR CO rotational emission from the prototypical Seyfert 2 galaxy NGC 1068. Using Herschel-PACS, we have detected 11 transitions in the J_upper=14-30 (E_upper/k_B = 580-2565 K) range, all of which are consistent with arising from within the central 10" (700 pc). The detected transitions are modeled as arising from 2 different components: a moderate excitation (ME) component close to the galaxy systemic velocity, and a high excitation (HE) component that is blueshifted by ~80 km s^{-1}. We employ a large velocity gradient (LVG) model and derive n_H2~10^{5.6} cm^{-3}, T_kin~170 K, and M_H2~10^{6.7} M_sun for the ME component, and n_H2~10^{6.4} cm^{-3}, T_kin~570 K, and M_H2~10^{5.6} M_sun for the HE component, although for both components the uncertainties in the density and mass are plus/minus (0.6-0.9) dex. We compare the CO line profiles with those of other molecular tracers observed at higher spatial and spectral resolution, and find that the ME transitions are consistent with these lines arising in the ~200 pc diameter ring of material traced by H_2 1-0 S(1) observations. The blueshift of the HE lines may also be consistent with the bluest regions of this H_2 ring, but a better kinematic match is found with a clump of infalling gas ~40 pc north of the AGN. We discuss the prospects of placing the HE component near the AGN, and conclude that while the moderate thermal pressure precludes an association with the ~1 pc radius H_2O maser disk, the HE component could potentially be located only a few parsecs more distant from the AGN, and might then provide the N_H~10^{25} cm^{-2} column obscuring the nuclear hard X-rays. Finally, we also report sensitive upper limits extending up to J_upper=50, which place constraints on a previous model prediction for the CO emission from the X-ray obscuring torus. [Abridged]

There is now convincing evidence that the intensity of HCN molecular line emission is enhanced around active galactic nuclei. In this paper we examine the specific case of the Seyfert galaxy NGC3227, for which there are subarcsecond resolution data for the HCN (1-0) 88 GHz and CO (2-1) 230 GHz rotational lines, enabling us to spatially separate a circumnuclear ring at a radius of 140pc and an inner nuclear region within 40pc of the AGN. The HCN(1-0)/CO(2-1) flux ratio differs by more than an order of magnitude between these two regions. We carry out large velocity gradient (LVG) computations to determine the range of parameters (gas temperature and density, HCN/CO abundance ratio, column densities and velocity gradients) that yield physically plausible solutions for the observed flux ratio in the central 100pc. The observed HCN/CO intensity ratio in the nucleus is consistent with very optically thick thermalized emission in very dense (>=10^5cm^{-3}) gas, in which case the HCN/CO abundance ratio there is unconstrained. Alternatively, the HCN/CO intensity ratio could be due to optically thinner emissions but with very high (~10^{-2}) HCN/CO abundance ratios. This possibility is more consistent with the CO and HCN emissions seen in the nuclei of the Seyfert galaxies NGC1068 and NGC6951. It would imply the velocity gradients are large and the clouds may be gravitationally unbound. We estimate that the X-ray ionisation rate at radii less than 20pc in the centre of NGC3227 exceeds 10^{-13}s^{-1}. X-ray ionisation and heating may lead to high HCN/CO ratios in warm gas in a high-ionisation molecular phase near the AGN.