The MINiature Exoplanet Radial Velocity Array (MINERVA) is a dedicated observatory of four 0.7m robotic telescopes fiber-fed to a KiwiSpec spectrograph. The MINERVA mission is to discover super-Earths in the habitable zones of nearby stars. This can be accomplished with MINERVA's unique combination of high precision and high cadence over long time periods. In this work, we detail changes to the MINERVA facility that have occurred since our previous paper. We then describe MINERVA's robotic control software, the process by which we perform 1D spectral extraction, and our forward modeling Doppler pipeline. In the process of improving our forward modeling procedure, we found that our spectrograph's intrinsic instrumental profile is stable for at least nine months. Because of that, we characterized our instrumental profile with a time-independent, cubic spline function based on the profile in the cross dispersion direction, with which we achieved a radial velocity precision similar to using a conventional "sum-of-Gaussians" instrumental profile: 1.8 m s$^{-1}$ over 1.5 months on the RV standard star HD 122064. Therefore, we conclude that the instrumental profile need not be perfectly accurate as long as it is stable. In addition, we observed 51 Peg and our results are consistent with the literature, confirming our spectrograph and Doppler pipeline are producing accurate and precise radial velocities.

We report a new analysis protocol for HCN hyperfine data, based on the PYSPECKIT package, and results of using this new protocol to analyse a sample area of seven massive molecular clumps from the Census of High- and Medium-mass Protostars (CHaMP) survey, in order to derive maps of column density for this species. There is a strong correlation between the HCN integrated intensity, IHCN, and previously reported IHCO+ in the clumps, but IN2H+ is not well correlated with either of these other two "dense gas tracers". The four fitted parameters from PYSPECKIT in this region fall in the range of VLSR = 8-10 km/s, {\sigma} V = 1.2-2.2 km/s, Tex = 4-15 K, and {\tau} = 0.2-2.5. These parameters allow us to derive a column density map of these clouds, without limiting assumptions about the excitation or opacity. A more traditional (linear) method of converting IHCN to total mass column gives much lower clump masses than our results based on the hyperfine analysis. This is primarily due to areas in the sample region of low I, low Tex, and high {\tau} . We conclude that there may be more dense gas in these massive clumps not engaged in massive star formation than previously recognized. If this result holds for other clouds in the CHaMP sample, it would have dramatic consequences for the calibration of the Kennicutt-Schmidt star formation laws, including a large increase in the gas depletion time-scale in such regions.

We report the second complete molecular line data release from the {\em Census of High- and Medium-mass Protostars} (CHaMP), a large-scale, unbiased, uniform mapping survey at sub-parsec resolution, of mm-wave line emission from 303 massive, dense molecular clumps in the Milky Way. This release is for all $^{12}$CO $J$=1$\rightarrow$0 emission associated with the dense gas, the first from Phase II of the survey, which includes $^{12}$CO, $^{13}$CO, and C$^{18}$O. The observed clump emission traced by both $^{12}$CO and HCO$^+$ (from Phase I) shows very similar morphology, indicating that, for dense molecular clouds and complexes of all sizes, parsec-scale clumps contain $\Xi$ ~ 75% of the mass, while only 25% of the mass lies in extended (>~ 10 pc) or "low density" components in these same areas. The mass fraction of all gas above a density 10$^9$ m$^{-3}$ is $\xi_9$ >~ 50%. This suggests that parsec-scale clumps may be the basic building blocks of the molecular ISM, rather than the standard GMC concept. Using $^{12}$CO emission, we derive physical properties of these clumps in their entirety, and compare them to properties from HCO$^+$, tracing their denser interiors. We compare the standard X-factor converting $I_{CO}$ to $N_{H_2}$ with alternative conversions, and show that only the latter give whole-clump properties that are physically consistent with those of their interiors. We infer that the clump population is systematically closer to virial equilibrium than when considering only their interiors, with perhaps half being long-lived (10s of Myr), pressure-confined entities which only terminally engage in vigorous massive star formation, supporting other evidence along these lines previously published.