We present 247-GHz molecular line observations of methanol (CH$_3$OH) toward sixteen massive star-forming regions, using the APEX telescope with an angular resolution of $25''$. The sample covers a range of evolutionary states, including warm molecular cores, hot molecular cores, and ultracompact HII regions. The hot cores, all of which include UC HII regions, show rich molecular line spectra, although the strength of different species and transitions varies from source to source. In contrast, the warm cores do not show significant molecular line emission. Multiple methanol transitions are detected toward nine of the hot cores; eight of these had enough transitions to use the rotation diagram method to estimate rotational temperatures and column densities. The temperatures lie in the range 104$-$168 K and column densities from $3\times10^{16}$ to $7\times10^{18}$ cm$^{-2}$. Using the average methanol line parameters, we estimate virial masses, which fall in the range from 145 to 720 M$_\odot$ and proved to be significantly higher than the measured gas masses. We discuss possible scenarios to explain the chemical differences between hot cores and warm molecular cores. One of the observed methanol lines, $4_{2}-5_{1}A^{+}$ at 247.228 GHz, is predicted to show class II maser emission, similar in intensity to previously reported $J_0-J_{-1}E$ masers at 157 GHz. We did not find any clear evidence for maser emission among the observed sources; however, a weak maser in this line may exist in G345.01+1.79.

We present Submillimeter Array observations, in the 1.3 mm continuum and the CH_3CN(12-11) line of 17 hot molecular cores associated with young high-mass stars. The angular resolution of the observations ranges from 1".0 to 4".0. The continuum observations reveal large (>3500 AU) dusty structures with gas masses from 7 to 375 Msun, that probably surround multiple young stars. The CH_3CN line emission is detected toward all the molecular cores at least up to the K=6-component and is mostly associated with the emission peaks of the dusty objects. We used the multiple K-components of the CH_3CN and both the rotational diagram method and a simultaneous synthetic LTE model with the XCLASS program to estimate the temperatures and column densities of the cores. For all sources, we obtained reasonable fits from XCLASS by using a model that combines two components: an extended and warm envelope, and a compact hot core of molecular gas, suggesting internal heating by recently formed massive stars. The rotational temperatures lie in the range of 40-132 K and 122-485 K for the extended and compact components, respectively. From the continuum and CH_3CN results, we infer fractional abundances from 10^{-9} to 10^{-7} toward the compact inner components, that increase with the rotational temperature. Our results agree with a chemical scenario in which the CH_3CN molecule is efficiently formed in the gas phase above 100-300 K, and its abundance increases with temperature.