As light intensity increases in photoacoustic imaging, the saturation of optical absorption and the temperature dependence of the thermal expansion coefficient result in a measurable nonlinear dependence of the photoacoustic (PA) signal on the excitation pulse fluence. Here, under controlled conditions, we investigate the intensity-dependent photoacoustic signals from oxygenated and deoxygenated hemoglobin at varied optical wavelengths and molecular concentrations. The wavelength and concentration dependencies of the nonlinear PA spectrum are found to be significantly greater in oxygenated hemoglobin than in deoxygenated hemoglobin. These effects are further influenced by the hemoglobin concentration. These nonlinear phenomena provide insights into applications of photoacoustics, such as measurements of average inter-molecular distances on a nm scale or with a tuned selection of wavelengths, a more accurate quantitative PA tomography. Photoacoustic microscopy (PAM) is an effective in vivo functional and molecular imaging tool. In the PA phenomenon, light is absorbed by molecules and converted to heat. Subsequent thermoelastic expansion generates an acoustic wave, termed the PA wave. 1,2 Detection of PA waves sequentially excited at multiple optical wavelengths provides quantitative information about the concentrations of multiple chromophores such as oxygenated and deoxygenated hemoglobin molecules in red blood cells. Thus, the relative concentration and the oxygen saturation (sO 2 ) of hemoglobin can be extracted. [3][4][5] Generally, the amplitude of the PA signal is assumed to be linearly proportional to the excitation pulse fluence. However, as the excitation laser intensity increases, both the saturation of the optical absorption 6,7 and the temperature dependence of thermal expansion 8-10 result in a measurable nonlinear dependence of the PA signal on the excitation pulse fluence. PA nonlinearity has recently been used in several applications such as quantifying picosecond absorption relaxation times with a nanosecond laser, 6 differentiating optical absorbers, 10 measuring oxygen saturation in vivo, 7 and performing label-free PA nanoscopy of biological structures having undetectable fluorescence. 11 In the presence of nonlinearity, quantitative PA measurements require a detailed analysis of the intensity-dependence of the PA signal, particularly for hemoglobin, the major intrinsic absorber in tissue for PA imaging. Based on thermal nonlinearity, analysis of the nonlinear PA signal from a single point source has been previously discussed. 8 However, the wavelengthand concentration-dependent effects of both optical saturation and thermal nonlinearity on the PA signals for multiple absorbers have not been reported. Here, we investigate nonlinear PA effects in oxygenated and deoxygenated hemoglobin using a PA spectrometer with a flat-top beam illumination, which effectively reduces uncertainty in our measurements that would otherwise arise from inhomogeneous spatial distribution of the optical fluence. In ox...