Photothermal excitation at frequencies below the mechanical resonance of the atomic force microscopy (AFM) microcantilever can be utilized in force modulation microscopy, fast force displacement curve acquisition, and tip-based mass spectroscopy. To understand the microcantilever bending response in these modes, accurate models of the thermoelastic response of the AFM microcantilever are needed. We study the sub-resonance photothermal vibrational response of coated and uncoated AFM microcantilevers as a function of laser modulation frequency and spot location. The sub-resonance microcantilever response shows distinct thermoelastic regimes. Below the microcantilever's thermal roll-off frequency, the vibration amplitude is mostly constant. Past this frequency, the vibration amplitude decreases with increasing frequency. At modulation frequencies below the thermal roll-off frequency, the most efficient photothermal laser spot to excite harmonic motion is near the free end of both coated and uncoated microcantilevers. For the tested coated microcantilevers, the most efficient photothermal laser location migrates from near the free end of the microcantilever to near the fixed end as modulation frequency increases. For the tested uncoated microcantilever, the most efficient photothermal laser location remains unchanged at the tested frequencies. To predict the bending response of the coated microcantilever, a bilayer bending model is implemented. At low frequencies, this model underpredicts the bending response compared to experiments by up to 90%. This may be due to neglecting microcantilever bending contributed by a through-thickness temperature gradient. Our results illustrate different aspects of the frequency-dependent photothermal laser spot optimization that can guide users to maximizing microcantilever response to a given input power. [ABSTRACT FROM AUTHOR]