Recent radio astronomical observations have revealed that HC$_{5}$N, the second shortest cyanopolyyne (HC$_{2n+1}$N), is abundant around some massive young stellar objects (MYSOs), which is not predicted by classical carbon-chain chemistry. For example, the observed HC$_{5}$N abundance toward the G28.28$-$0.36 MYSO is higher than that in L1527, which is one of the warm carbon chain chemistry (WCCC) sources, by more than one order of magnitude (Taniguchi et al., 2017). In this paper, we present chemical simulations of hot-core models with a warm-up period using the astrochemical code Nautilus. We find that the cyanopolyynes are formed initially in the gas phase and accreted onto the bulk and surface of granular ice mantles during the lukewarm phase, which occurs at $25 < T < 100$ K. In slow warm-up period models, the peak abundances occur as the cyanopolyynes desorb from dust grains after the temperature rises above 100 K. The lower limits of the abundances of HC$_{5}$N, CH$_{3}$CCH, and CH$_{3}$OH observed in the G28.28$-$0.36 MYSO can be reproduced in our hot-core models, after their desorption from dust grains. Moreover, previous observations suggested chemical diversity in envelopes around different MYSOs. We discuss possible interpretations of relationships between stages of the star-formation process and such chemical diversity, such as the different warm-up timescales. This timescale depends not only on the mass of central stars but also on the relationship between the size of warm regions and their infall velocity.

The emission line spectra of cyanoacetylene and methanol reveal chemical and physical heterogeneity on very small (< 0.1 pc) scales toward the peak in cyanopolyyne emission in the Taurus Molecular Cloud, TMC-1 (CP). We generate grids of homogeneous chemical models using a three-phase rate equation approach to obtain all time-dependent abundances spanning the physical conditions determined from molecular tracers of compact and extended regions of emission along this line of sight. Each time-dependent abundance is characterized by one of four features: a maximum/minimum, a monotonic increase/decrease, oscillatory behavior, or inertness. We similarly classify the time-dependent agreement between modeled and observed abundances by calculating both the root-mean-square logarithm difference and root-mean-square deviation between the modeled and observed abundances at every point in our grid models for three groups of molecules: (i) a composite group of all species present in both the observations and our chemical network G, (ii) the cyanopolyynes C = {HC$_3$N, HC$_5$N, HC$_7$N, HC$_9$N}, and (iii) the oxygen-containing organic species methanol and acetaldehyde S = {CH$_3$OH, CH$_3$CHO}. We discuss how the Bayesian uncertainties in the observed abundances constrain solutions within the grids of chemical models. The calculated best fit times at each grid point for each group are tabulated to reveal the minimum solution space of the grid models and the effects the Bayesian uncertainties have on the grid model solutions. The results of this approach separate the effect different physical conditions and model-free parameters have on reproducing accurately the abundances of different groups of observed molecular species.