Aiming for a new and more comprehensive DIB catalog between 4,000 and 9,000{\AA}, we revisited the Atlas Catalog based on the observations of HD 183143 and HD 204827 (Hobbs et al. 2008 and 2009). Twenty-five medium-to-highly reddened sight lines were selected, sampling a variety of spectral types of the background star and the interstellar environments. The median SNR of these spectra is ~ 1,300 around 6,400{\AA}. Compared to the Atlas Catalog, 22 new DIBs were found, and the boundaries of 27 (sets) of DIBs were adjusted, resulting in an updated catalog containing 559 DIBs that we refer to as the Apache Point Observatory Catalog of Optical Diffuse Interstellar Bands. Measurements were then made based on this catalog. We found our survey most sensitive between 5,500 and 7,000{\AA} due largely to the local SNR of the spectra, the relative absence of interfering stellar lines, and the weakness of telluric residuals. For our data sample, the number of DIBs detected in a given sight line is mostly dependent on EB-V and less on the spectral type of the background star. Some dependence on fH2 is observed, but less well-determined due to the limited size of the data sample. The variations of the wavelengths of each DIB in different sight lines are generally larger than those of the known interstellar lines CH+, CH, and K I. Those variations could be due to the inherent error in the measurement, or to differences in the velocity components among sight lines.

We study the behavior of eight diffuse interstellar bands (DIBs) in different interstellar environments, as characterized by the fraction of hydrogen in molecular form [$f$(H$_2$)], with comparisons to the corresponding behavior of various known atomic and molecular species. The equivalent widths of the five "normal" DIBs ($\lambda\lambda$5780.5, 5797.1, 6196.0, 6283.8, and 6613.6), normalized to $E(B-V)$, show a "Lambda-shaped" behavior: they increase at low $f$(H$_2$), peak at $f$(H$_2$) ~ 0.3, and then decrease. The similarly normalized column densities of Ca, Ca$^+$, Ti$^+$, and CH$^+$ also decline for $f$(H$_2$) > 0.3. In contrast, the normalized column densities of Na, K, CH, CN, and CO increase monotonically with $f$(H$_2$), and the trends exhibited by the three C$_2$ DIBs ($\lambda\lambda$4726.8, 4963.9, and 4984.8) lie between those two general behaviors. These trends with $f$(H$_2$) are accompanied by cosmic scatter, the dispersion at any given $f$(H$_2$) being significantly larger than the individual errors of measurement. The Lambda-shaped trends suggest the balance between creation and destruction of the DIB carriers differs dramatically between diffuse atomic and diffuse molecular clouds; additional processes besides ionization and shielding are needed to explain those observed trends. Except for several special cases, the highest $W$(5780)/$W$(5797) ratios, characterizing the so-called "sigma-zeta effect", occur only at $f$(H$_2$) < 0.2. We propose a sequence of DIBs based on trends in their pair-wise strength ratios with increasing $f$(H$_2$). In order of increasing environmental density, we find the $\lambda$6283.8 and $\lambda$5780.5 DIBs, the $\lambda$6196.0 DIB, the $\lambda$6613.6 DIB, the $\lambda$5797.1 DIB, and the C$_2$ DIBs.