In the thermal spectral range, there are millions of individual absorption lines of water vapor, CO2, and other trace gases. Radiative transfer calculations of wavelength-integrated quantities, such as irradiance and heating rate, are computationally expensive, requiring a high spectral resolution for accurate numerical weather prediction and climate modeling. This paper introduces a method that could highly reduce the cost of integration in the thermal spectrum by employing an optimized wavelength sampling method. Absorption optical thicknesses for various trace gases were calculated from the HITRAN 2012 spectroscopic dataset using the ARTS line-by-line model as input to a fast Schwarzschild radiative transfer model. Using a simulated annealing algorithm, different optimized sets of wavelengths and corresponding weights were identified, which allowed for accurate integrated quantities to be computed as a weighted sum, reducing the computational time by several orders of magnitude. For each set of wavelengths, a lookup table, including the corresponding weights and absorption cross-sections, is created and can be applied to any atmospheric setups for which it was trained. We applied the lookup table to calculate irradiances and heating rates for a large set of atmospheric profiles from the ECMWF 91-level short-range forecast. Ten wavelength nodes are sufficient to obtain irradiances within an average root mean square error (RMSE) of upward and downward radiation at any height below 1 Wm−2 while 100 wavelengths allowed for an RSME of below 0.05 Wm−2. The applicability of this method was confirmed for irradiances and heating rates in clear conditions and for an exemplary cloud at 3.2 km height. Representative spectral gridpoints for integrated quantities in the thermal spectrum (REPINT) is available as absorption parameterization in the libRadtran radiative transfer package, where it can be used as an efficient molecular absorption parameterization for a variety of radiative transfer solvers.