Measurement Precision and Thermal and Absorption Properties of Nanostructures in Aqueous Solutions by Transient and Steady-State Thermal-Lens Spectrometry

: A simultaneous steady-state and transient photothermal-lens modality was used for both the thermal and optical parameters of aqueous dispersed systems (carbon and silica nanoparticles, metal iodides, surfactants, heme proteins, albumin, and their complexes). Heat-transfer parameters (thermal diffusivity and thermal effusivity), the temperature gradient of the refractive index, light absorption, and concentration parameters were assessed. To simultaneously measure thermal and optical parameters, the time scale of thermal lensing (characteristic time, t c ) should correspond to an excitation beam size of 60–300 µ m, and the relative time intervals ( 0.5 ÷ 5 ) t c and ( 5 ÷ 20 ) t c should be selected for transient and steady-state measurements, respectively. Dual-beam thermal-lens spectrometers in a mode-mismatched optical schematic at various excitation wavelengths were built. The spectrometers implement back-synchronized detection, providing different measurement conditions for the heating and cooling parts of the thermal-lens cycle. By varying the measurement parameters depending on the dispersed system, the conditions providing the suitable precision (replicability, repeatability, and reproducibility) of thermal-lens measurements were found; setups with a broad excitation beam (waist size, 150 and 300 µ m) provide longer times to attain a thermal equilibrium and, thus, the better precision of measurements of thermal diffusivity.

For thermal-lens measurements of dispersed solutions, the spatial distribution of thermal properties in the medium is crucial. Thus, the TLS application for dispersed systems is usually based on high-precision microscopic techniques combined with data visualization and sophisticated data treatment [32,[83][84][85][86][87][88][89][90][91][92]. Still, TLS can be implemented with relatively simple dual-beam setups that provide the required sensitivity in rather large-volume bulk samples [2,12,16,93,94]. This approach is considered relatively straightforward for the determination of optical absorption and analyte assessment (similarly to transmission photometry) and is used for the estimation of these properties of solvents. Still, the complex and dynamic nature of photothermal effects provides very informative data, which can be used, even with a non-microscopic schematic, to assess more detailed information: absorption, thermophysical properties, and size parameters [4,50,[95][96][97].
A thermal-lens schematic includes a unit to synchronize the start of the blooming of the thermal-lens effect by the excitation beam and a pulse from the excitation laser or a modulated (chopped) cw beam. The principles of the most common thermal-lens schematic can be described as a directly synchronized (lock-in) TLS technique, as the primary trigger for a new measurement cycle is a new laser pulse or a modulator opening, and the measurements are made at a constant frequency determined by the pulse repetition rate or a modulator [2,89,94]. This schematic is simple and reliable and provides high sensitivity [86,98]. However, for TLS, single-frequency synchronization may not capture information for complex samples, in which the transient and steady-state responses depend on both optical and thermal heterogeneity [2,99]. Furthermore, for dispersed solutions at a constant frequency, the measurement error may accumulate over time due to incomplete heat dissipation from cycle to cycle, providing a lower temperature response [2,28]. This, in turn, leads to thermally induced processes that result in an uncontrolled contribution to probe-beam blooming.
The back-synchronized measurement excitation/data-treatment technique differs at the start of a new measurement cycle; its start is governed by reaching steady-state thermal-lens conditions [100,101]. In this technique, the advantages are (i) the possibility of simultaneously using transient and repeated steady-state measurements; (ii) a linear dynamic range of more than five orders of magnitude of detectable absorbances, including strongly absorbing and scattering samples; and (iii) the possibility of detection under batch and flow conditions with no change in the optical-scheme design of the instrument [102]. Although the detection sensitivity is somewhat lower compared to more common lock-in schemes [2,89,94], the flexibility of back-synchronized TLS provides a considerably larger information volume for studies of complex formation at trace concentrations or transient heat dynamics around absorbing nanoparticles [84].
In this study, we implemented a simultaneous steady-state and transient photothermallens modality for various aqueous dispersions using back-synchronized setups with different time-response scales. The aim was to demonstrate the capability of such a modality for measuring the absorbance and assessing the concentrations of dispersed solutions, along with monitoring the thermal properties (thermal diffusivity and the thermo-optical constant) with the selection of measurement conditions based on precision parameters.
instrument [102]. Although the detection sensitivity is somewhat lower compared to more common lock-in schemes [2,89,94], the flexibility of back-synchronized TLS provides a considerably larger information volume for studies of complex formation at trace concentrations or transient heat dynamics around absorbing nanoparticles [84].
In this study, we implemented a simultaneous steady-state and transient photothermal-lens modality for various aqueous dispersions using back-synchronized setups with different time-response scales. The aim was to demonstrate the capability of such a modality for measuring the absorbance and assessing the concentrations of dispersed solutions, along with monitoring the thermal properties (thermal diffusivity and the thermo-optical constant) with the selection of measurement conditions based on precision parameters.
The excitation beam is directed to a PD L-3DP3C photodiode (Photometry PD) to account for the excitation beam photometric (transmission) signal (Appendix A). Additionally, this mirror prevents the blooming of the thermal lens in the following bandpass filter. The excitation beam is absorbed by the filter (KS-11 color glass, 2 mm depth). The signal is converted and amplified and enters the ADC-DAC unit. The probe beam passes through a pinhole (diameter, 4 mm) centered at the optical axis and hits the Primary PD L-3DP3C photodiode. The signal is converted and amplified by the control unit and enters Channel 1 of the ADC-DAC unit. Through an amplifier that converts a photocurrent into voltage, the signal from this photodiode enters channel 2 of the ADC and DAC converter, connected to a PC.
Lenses L1 and L2 and the sample compartment can be moved along the directions of the beams (step, 0.2 mm), which provides variation in the spectrometer geometry (degree of mode mismatching, see below). The spectrometer implements a secondary channel (an L-3DP3C photodiode, Secondary PD, Figure 1) for gathering scattering or luminescence signals if present.
A homemade ADC-DAC unit based on a C8051F061 board (two 16-bit ADC and two 12-bit DAC channels, ADC time, 2 µs; readout frequency, 1-5 kHz; C8051F060DK development kit, Silicon Labs, Boston, MA, USA) was used in the external-trigger mode from the PC [100]. The homemade software communicates with the unit through an RS-232 interface. The homemade software is written in C++, Borland C++ ver. 5 (Borland Corp., Austin, TX, USA). After obtaining a full curve of repeated heating-cooling cycles, it is divided into separate cycles; each cycle is divided into heating and cooling half-curves, and each half-cycle curve is divided into transient and steady-state parts. A steady-state signal is calculated automatically during measurements, and the transient curves are treated after the experiment.

Figure 2.
Single excitation on/off cycle: a normalized thermal-lens curve, Equation (6), in backsynchronized mode for measurements of transient and steady-state thermal-lens signals exemplified by ferroin in NMP; TLS-60; excitation wavelength, 514.5 nm; excitation power, 25 mW. The selection of times t1, t2, t3, and t4 in this figure is for the illustration of the principle only. For the cooling part of the cycle, times t3 and t4 are measured from the shutter-off trigger.

Transient Thermal-Lens Measurements
A transient (time-resolved) thermal-lens signal ϑ(t) for an excitation on/off cycle of a modulated continuous-wave (cw) beam was calculated as the relative change in the probe-beam intensity at a far-field detector plane ( ) at moment t [2]: where (0) is the intensity of the probe beam at the photodetector plane in the central part of the beam at time t = 0, and ( ) is the intensity of the probe beam at time t from Figure 2. Single excitation on/off cycle: a normalized thermal-lens curve, Equation (6), in backsynchronized mode for measurements of transient and steady-state thermal-lens signals exemplified by ferroin in NMP; TLS-60; excitation wavelength, 514.5 nm; excitation power, 25 mW. The selection of times t 1 , t 2 , t 3 , and t 4 in this figure is for the illustration of the principle only. For the cooling part of the cycle, times t 3 and t 4 are measured from the shutter-off trigger.

Transient Thermal-Lens Measurements
A transient (time-resolved) thermal-lens signal ϑ(t) for an excitation on/off cycle of a modulated continuous-wave (cw) beam was calculated as the relative change in the probe-beam intensity at a far-field detector plane I p (t) at moment t [2]: where I p (0) is the intensity of the probe beam at the photodetector plane in the central part of the beam at time t = 0, and I p (t) is the intensity of the probe beam at time t from the excitation on/off cycle start. According to the diffraction model [104,105], the equation for the thermal-lens signal in relation to light absorption and concentration is REVIEW 6 the excitation on/off cycle start. According to the diffraction model [104,105], the equation for the thermal-lens signal in relation to light absorption and concentration is where is the excitation laser power, ϰ = 2.303 is the linear light-absorption coefficient, l is the sample path length, ε is molar absorptivity, and с is molar concentration of the test absorbing substance in the sample. For absorbances up to ca. 0.005, excitation powers up to ca. 200 mW, and, thus, small temperature changes in the thermal-lens effect, Equation (2) The factor in Equations (2) and (3) is the power-normalized thermo-optical constant (full thermo-optical constant is = ), characterizing the strength of the thermal-lens effect per unit excitation laser power, . Here, k is the thermal conductivity, is the probe-laser wavelength, and the temperature gradient of the refractive index ( ⁄ ) ∝ α (volumetric thermal expansion coefficient).
The parameter ( ) in Equations (2) and (3) is the time-dependent geometrical constant of the optical schematic [104,105], Physchem 2023, 3 161 where P e is the excitation laser power, where is the excitation laser power, ϰ = 2.303 is the linear light-absorption coefficient, l is the sample path length, ε is molar absorptivity, and с is molar concentration of the test absorbing substance in the sample. For absorbances up to ca. 0.005, excitation powers up to ca. 200 mW, and, thus, small temperature changes in the thermal-lens effect, Equation (2) can be simplified due to the small contributions of the quadratic terms and can be defined as a linear equation (see Figures S5 and S6, Supplementary Information): The factor in Equations (2) and (3) is the power-normalized thermo-optical constant (full thermo-optical constant is = ), characterizing the strength of the thermal-lens effect per unit excitation laser power, . Here, k is the thermal conductivity, is the probe-laser wavelength, and the temperature gradient of the refractive index ( ⁄ ) ∝ α (volumetric thermal expansion coefficient).
The parameter ( ) in Equations (2) and (3) is the time-dependent geometrical constant of the optical schematic [104,105], Here, m is the ratio of the cross-sectional areas of the probe and excitation beams at the sample, and V is the relative distance from the excitation waist to the sample, from the traditional notation of TLS geometry parameters proposed by R. Snook, J. Shen, and M. Baesso [104,105], and a, b, and c are related geometry constants of the spectrometer used in the approximation for equation simplicity's sake [106]. The parameter in Equation (4) is the characteristic time of the thermal lens [2] and is calculated from the transient curve of thermal-lens measurements ( Figure 2).
where is the excitation beam waist radius, and is thermal diffusivity. To compare the transient heating and cooling curves under various conditions, we used the relative intensity of the probe beam [107]: where (∞) is the intensity of the probe beam at thermal equilibrium due to the steadystate thermal-lens effect (t→∞). This form of Equation (1) makes it possible to compare the blooming and dissipation of the curves for various media and with various light absorption characteristics, as all of the curves are normalized to a scale from 0 (fully developed steady-state thermal lens) to 1 (no photothermal effects), thus focusing on differences in the vicinity of the characteristic time of the process.

Steady-State Measurements
The steady-state cw thermal-lens signal is and, under the same conditions as linear Equation (3), is defined as: Here, θ ≡ 2.303 (9) = 2.303εc is the linear light-absorption coefficient, l is the sample path length, ε is molar absorptivity, and c is molar concentration of the test absorbing substance in the sample. For absorbances up to ca. 0.005, excitation powers up to ca. 200 mW, and, thus, small temperature changes in the thermal-lens effect, Equation (2) can be simplified due to the small contributions of the quadratic terms and can be defined as a linear equation (see Figures S5 and S6, Supplementary Information): the excitation on/off cycle start. According to the diffraction model [104,105], the equation for the thermal-lens signal in relation to light absorption and concentration is where is the excitation laser power, ϰ = 2.303 is the linear light-absorption coefficient, l is the sample path length, ε is molar absorptivity, and с is molar concentration of the test absorbing substance in the sample. For absorbances up to ca. 0.005, excitation powers up to ca. 200 mW, and, thus, small temperature changes in the thermal-lens effect, Equation (2) can be simplified due to the small contributions of the quadratic terms and can be defined as a linear equation (see Figures S5 and S6, Supplementary Information): The factor in Equations (2) and (3) is the power-normalized thermo-optical constant (full thermo-optical constant is = ), characterizing the strength of the thermal-lens effect per unit excitation laser power, . Here, k is the thermal conductivity, is the probe-laser wavelength, and the temperature gradient of the refractive index ( ⁄ ) ∝ α (volumetric thermal expansion coefficient).
The parameter ( ) in Equations (2) and (3) is the time-dependent geometrical constant of the optical schematic [104,105], Here, m is the ratio of the cross-sectional areas of the probe and excitation beams at the sample, and V is the relative distance from the excitation waist to the sample, from the traditional notation of TLS geometry parameters proposed by R. Snook, J. Shen, and M. Baesso [104,105], and a, b, and c are related geometry constants of the spectrometer used in the approximation for equation simplicity's sake [106].
The parameter in Equation (4) is the characteristic time of the thermal lens [2] and is calculated from the transient curve of thermal-lens measurements ( Figure 2).
where is the excitation beam waist radius, and is thermal diffusivity.
To compare the transient heating and cooling curves under various conditions, we used the relative intensity of the probe beam [107]: where (∞) is the intensity of the probe beam at thermal equilibrium due to the steady-state thermal-lens effect (t→∞). This form of Equation (1) makes it possible to compare the blooming and dissipation of the curves for various media and with various light absorption characteristics, as all of the curves are normalized to a scale from 0 (fully developed steady-state thermal lens) to 1 (no photothermal effects), thus focusing on differences in the vicinity of the characteristic time of the process.

Steady-State Measurements
The steady-state cw thermal-lens signal is and, under the same conditions as linear Equation (3), is defined as: Here, θ ≡ 2.303 The factor E 0 in Equations (2) and (3) is the power-normalized thermo-optical constant (full thermo-optical constant is E = P e E 0 ), characterizing the strength of the thermal-lens effect per unit excitation laser power, E 0 = (−dn/dT)/ λ p k . Here, k is the thermal conductivity, λ p is the probe-laser wavelength, and the temperature gradient of the refractive index (dn/dT) ∝ α T (volumetric thermal expansion coefficient).
The parameter B(t) in Equations (2) and (3) is the time-dependent geometrical constant of the optical schematic [104,105], Here, m is the ratio of the cross-sectional areas of the probe and excitation beams at the sample, and V is the relative distance from the excitation waist to the sample, from the traditional notation of TLS geometry parameters proposed by R. Snook, J. Shen, and M. Baesso [104,105], and a, b, and c are related geometry constants of the spectrometer used in the approximation for equation simplicity's sake [106].
The parameter t c in Equation (4) is the characteristic time of the thermal lens [2] and is calculated from the transient curve of thermal-lens measurements ( Figure 2).
where ω 0e is the excitation beam waist radius, and D T is thermal diffusivity.
To compare the transient heating and cooling curves under various conditions, we used the relative intensity of the probe beam [107]: where I p (∞) is the intensity of the probe beam at thermal equilibrium due to the steadystate thermal-lens effect (t→∞). This form of Equation (1) makes it possible to compare the blooming and dissipation of the curves for various media and with various light absorption characteristics, as all of the curves are normalized to a scale from 0 (fully developed steadystate thermal lens) to 1 (no photothermal effects), thus focusing on differences in the vicinity of the characteristic time t c of the process.

Steady-State Measurements
The steady-state cw thermal-lens signal is and, under the same conditions as linear Equation (3), is defined as: Physchem 2023, 3, FOR PEER REVIEW 6 the excitation on/off cycle start. According to the diffraction model [104,105], the equation for the thermal-lens signal in relation to light absorption and concentration is where is the excitation laser power, ϰ = 2.303 is the linear light-absorption coefficient, l is the sample path length, ε is molar absorptivity, and с is molar concentration The factor in Equations (2) and (3) is the power-normalized thermo-optical constant (full thermo-optical constant is = ), characterizing the strength of the thermal-lens effect per unit excitation laser power, is the probe-laser wavelength, and the temperature gradient of the refractive index ( ⁄ ) ∝ α (volumetric thermal expansion coefficient).
The parameter ( ) in Equations (2) and (3) is the time-dependent geometrical constant of the optical schematic [104,105], ]. (4) Here, θ ≡ 2.303E 0 P e εlc (9) is the steady-state thermal-lens signal corrected for the steady-state geometry constant B(t→∞). The calibration plots for TLS were built as θ ∝ c and as a comparison of thermallens and transmission measurements (photometry), with where κ is a proportionality coefficient, the theoretical value from Equation (9) is 2.303E 0 P e , and A = εlc is absorbance from photometric measurements. The experimental values of the signal θ were corrected, if necessary, for a decrease in the excitation power due to light-scattering losses A s : The recalculation of the absorbance from photothermal measurements (A PT ) was calculated from θ corr , Equation (11), and the equation deduced from Equation (9) and Beer's law: Molar absorptivities ε PT from photothermal measurements were calculated from Beer's law, A PT = ε PT cl. The experimental values of sample absorbance A exp (a spectrophotometry or photometry channel, Appendix A) were corrected for scattering: For known values of absorbance from absorption spectra obtained by spectrophotometry, the theoretical increase in the calibration slope for thermal lensing was calculated as

Parameters of Thermal-Lens Measurements
The data of 100-1000 replicate experimental transient thermal-lens curves I p (t) were measured and averaged. Next, the averaged curves were approximated as [104,105] where t 0 is the delay time, which depends on the transfer open/close times of the mechanical shutter and found empirically from the actual start of the change in the probe-beam intensity. Additionally, we used a previously developed numerical description for signal generation for multipoint-absorbing (disperse) solutions [96,108] (Appendix B). From the viewpoint of non-equilibrium (time-dependent) processes with characteristic times differing from that of the neat solvent, Equation (5) was solved for the effective time-dependent characteristic timet c (t): where the auxiliary function is X(t) ≡ 2 I(t) + 1 /θ. The effective thermal diffusivity as a function of time, D T (t), was calculated from Equations (5) and (16) as t c (t) curves show two types of periodic noises: short-time periodic noise due to fluctuations of the probe-laser-beam intensity and periodic beating with a frequency of 15 Hz caused by the vibrations imposed on shutter opening/closing, which become visible att c (t), especially for concentrated dispersions. Thus, the calculated curvest c (t) and D T (t) were subjected to two-stage smoothing: first, weighted adjacent averaging (5 points) to remove short-scale fluctuations, followed by Savitzky-Golay smoothing to remove the periodic noise. This two-stage algorithm enabled the preservation of the general curve shape but affected the starting point and the maximum value of the curve.
The equation of temperature dependence in the thermal-lens center on the cooling (thermal-lens dissipation) curve for a homogeneous medium is as follows: [2]: REVIEW 6 the excitation on/off cycle start. According to the diffraction model [104,105], the equation for the thermal-lens signal in relation to light absorption and concentration is where is the excitation laser power, ϰ = 2.303 is the linear light-absorption coefficient, l is the sample path length, ε is molar absorptivity, and с is molar concentration of the test absorbing substance in the sample. For absorbances up to ca. 0.005, excitation powers up to ca. 200 mW, and, thus, small temperature changes in the thermal-lens effect, Equation (2) can be simplified due to the small contributions of the quadratic terms and can be defined as a linear equation (see Figures S5 and S6, Supplementary Information): The factor in Equations (2) and (3) is the power-normalized thermo-optical constant (full thermo-optical constant is = ), characterizing the strength of the thermal-lens effect per unit excitation laser power, . Here, k is the thermal conductivity, is the probe-laser wavelength, and the temperature gradient of the refractive index ( ⁄ ) ∝ α (volumetric thermal expansion coefficient).
The parameter ( ) in Equations (2) and (3) is the time-dependent geometrical constant of the optical schematic [104,105], Here, m is the ratio of the cross-sectional areas of the probe and excitation beams at the sample, and V is the relative distance from the excitation waist to the sample, from the traditional notation of TLS geometry parameters proposed by R. Snook, J. Shen, and M. Baesso [104,105], and a, b, and c are related geometry constants of the spectrometer used in the approximation for equation simplicity's sake [106]. The parameter in Equation (4) is the characteristic time of the thermal lens [2] and is calculated from the transient curve of thermal-lens measurements ( Figure 2).
where is the excitation beam waist radius, and is thermal diffusivity. To compare the transient heating and cooling curves under various conditions, we used the relative intensity of the probe beam [107]: where (∞) is the intensity of the probe beam at thermal equilibrium due to the steadystate thermal-lens effect (t→∞). This form of Equation (1) makes it possible to compare the blooming and dissipation of the curves for various media and with various light absorption characteristics, as all of the curves are normalized to a scale from 0 (fully developed steady-state thermal lens) to 1 (no photothermal effects), thus focusing on differences in the vicinity of the characteristic time of the process.

Steady-State Measurements
The steady-state cw thermal-lens signal is and, under the same conditions as linear Equation (3), is defined as: Here, where P e is the power of the excitation laser, k is thermal conductivity, t heat is the time of the heating period ( Figure 2), and hyschem 2023, 3, FOR PEER REVIEW 6 the excitation on/off cycle start. According to the diffraction model [104,105], the equation for the thermal-lens signal in relation to light absorption and concentration is where is the excitation laser power, ϰ = 2.303 is the linear light-absorption coefficient, l is the sample path length, ε is molar absorptivity, and с is molar concentration of the test absorbing substance in the sample. For absorbances up to ca. 0.005, excitation powers up to ca. 200 mW, and, thus, small temperature changes in the thermal-lens effect, Equation (2) can be simplified due to the small contributions of the quadratic terms and can be defined as a linear equation (see Figures S5 and S6, Supplementary Information): The factor in Equations (2) and (3) is the power-normalized thermo-optical constant (full thermo-optical constant is = ), characterizing the strength of the thermal-lens effect per unit excitation laser power, . Here, k is the thermal conductivity, is the probe-laser wavelength, and the temperature gradient of the refractive index ( ⁄ ) ∝ α (volumetric thermal expansion coefficient).
The parameter ( ) in Equations (2) and (3) is the time-dependent geometrical constant of the optical schematic [104,105], Here, m is the ratio of the cross-sectional areas of the probe and excitation beams at the sample, and V is the relative distance from the excitation waist to the sample, from the traditional notation of TLS geometry parameters proposed by R. Snook, J. Shen, and M. Baesso [104,105], and a, b, and c are related geometry constants of the spectrometer used in the approximation for equation simplicity's sake [106]. The parameter in Equation (4) is the characteristic time of the thermal lens [2] and is calculated from the transient curve of thermal-lens measurements ( Figure 2).
where is the excitation beam waist radius, and is thermal diffusivity. To compare the transient heating and cooling curves under various conditions, we used the relative intensity of the probe beam [107]: where (∞) is the intensity of the probe beam at thermal equilibrium due to the steadystate thermal-lens effect (t→∞). This form of Equation (1) makes it possible to compare the blooming and dissipation of the curves for various media and with various light absorption characteristics, as all of the curves are normalized to a scale from 0 (fully developed steady-state thermal lens) to 1 (no photothermal effects), thus focusing on differences in the vicinity of the characteristic time of the process.

Steady-State Measurements
The steady-state cw thermal-lens signal is and, under the same conditions as linear Equation (3), is defined as: Here, θ ≡ 2.303 (9) is the linear absorption coefficient. For calculations of photothermal parameters for dispersed systems, we used the data given in [2], as well as a previously developed theoretical description for signal generation in thermal-lens spectrometry for multipoint-absorbing (disperse) solutions [96].
The volumetric heat capacity ρC p was calculated from the experimentally determined density and isobaric specific heat data. The thermal conductivity k and thermal effusivity e T (thermal inertia, heat penetration coefficient) were calculated from the photothermally measured thermal diffusivity and the specific heat capacity and density:

Precision Parameters
Data analysis was completed using STATISTICA 10 software (StatSoft, Tulsa, OK, USA). The normal distribution of the results was checked with the Shapiro-Wilk test. A distribution was considered to be normal if p ≤ 0.05. Simple variance analysis-the separation of the total standard deviation into repeatability (deviation between replicate measurements) and temporal (day-to-day deviation) parts-was performed according to [109].

Repeatability
A single jth measurement of a sample consisted of a sequence of r = 50-100 excitation on/off cycles, giving a series of results [ϑ j1 , ϑ j2 , . . . , ϑ ji , . . . , ϑ jr ] calculated from Equation (1) or (7). The following parameters were measured: the average signal for the sample, ϑ j , and repeatability RSD rep (measurement stability from excitation on/off cycle to cycle in a short period of measurements, 2-20 min), according to ISO 5725-4 [110]:

Replicability
For a single concentration of the analyte c, n = 10-30 replicate sample analyses were conducted by changing the solution in the cell, giving a series of values θ 1 , θ 2 , . . . θ j . . . θ n calculated from Equation (9). The following parameters were calculated: the average sample signal θ(c) and relative standard deviation (replicability) RSD rpl (c) for a given concentration:  [109,110], which establishes the interconnection of the variance V Y of the function Y = Y(x 1 , x 2 . . . x k ) with variances V x I of its independent variables. The equations for the thermal-lens signal, Equations (1) and (7), for two independent variables I p (t = 0) and I p (t = ∞), can be used to obtain e error of measurements of probe light intensity by a photodetector, ( ) d by Equation (3), θ is determined by Equation (9), and ( ) = 2.303 ( )θ. ady-state signals, the descending reproducibility curve is described by a RSD = θ , where is a constant depending on the spectrometer e, and is a medium-dependent constant in the range 0.8-1.1 [111]. sults of thermal-lens measurements at different concentrations of test ere obtained; 5-7 series of experiments were carried out (by changing the he cell), and replicability, RSD rpl ( ), Equation (22), was determined for 3-5 fferent days (reproducibility, RSD , Equation (23)). The results obtained ed to variance analysis. ults of thermal-lens measurements (cobalt complexes and ferroin) in various at different concentrations were obtained; 5-7 experiments were carried out the solution in the cell, and replicability, RSD rpl ( ) , Equation (22), was for 20-30 consecutive measurements of the signal (each measurement, 10-15 ability, RSD , Equation (21)). The results obtained were subjected to lysis.   [97]; analytical-grade sodium acetate; hydroxide (high-purity grade); potassium iodide; ascorbic acid tical grade); sodium sulfite; copper(I) chloride; copper(II) acetate; lead(II) ium dihydroantimonate hexahydrate; tin(II) chloride; sulfuric acid (highe); nitric acid (analytical grade); sodium perchlorate; hydrochloric acid; ade potassium dichromate; analytical-grade acetic acid, 99.8%; a standard r solution, pH 8.50; chemically pure HCl (1 + 20); and chemically pure ydroxide. Unless otherwise stated, all reagents were of cp grade and from for two independent variables Ip(t = 0) and Ip(t = ¥), can be used to obtain ϒ = This results in an equation suitable for various thermal-lens modalities [95,98]: where δI is the error of measurements of probe light intensity by a photodetector, ( ) is determined by Equation (3), θ is determined by Equation (9), and ( ) = 2.303 ( )θ. For low steady-state signals, the descending reproducibility curve is described by a hyperbola, RSD = θ , where is a constant depending on the spectrometer optical scheme, and is a medium-dependent constant in the range 0.8-1.1 [111]. The results of thermal-lens measurements at different concentrations of test substances were obtained; 5-7 series of experiments were carried out (by changing the solution in the cell), and replicability, RSD rpl ( ), Equation (22), was determined for 3-5 series on different days (reproducibility, RSD , Equation (23)). The results obtained were subjected to variance analysis.
The results of thermal-lens measurements (cobalt complexes and ferroin) in various solvents and at different concentrations were obtained; 5-7 experiments were carried out by changing the solution in the cell, and replicability, RSD rpl ( ) , Equation (22), was determined for 20-30 consecutive measurements of the signal (each measurement, 10-15 min; repeatability, RSD , Equation (21)). The results obtained were subjected to variance analysis.

Performance Parameters
The measurement results are presented as mean values with confidence intervals (p < 0.05 or p = 0.95) under the requirements for the competence of testing and calibration laboratories in ISO/IEC 17025:2005. Calibration parameters were calculated according to IUPAC 1998 recommendations to present the results of the chemical analysis.
The limits of detection (LOD) and quantification (LOQ) were calculated as 3σ and 10σ criteria, respectively, according to IUPAC 1998 recommendations. The variances were compared using the Fisher, Bartlett, and Cochran criteria to determine the significance of contributions to the signal in control experiments.

Reagents and Chemicals
Water from a Milli-Q water purification system (Millipore SAS, Molsheim, France) was used: pH 6.8; specific resistance 18.2 MΩ × cm, Fe, 2 ppt; dissolved SiO2, 3 ppb; total ion amount <0.2 ppb; and TOC <10 ppb. The glassware was washed with acetone, followed by conc. nitric acid. All solutions were prepared daily before the experiments.
The following reagents were used throughout (unless otherwise stated, the products were from Reakhim, Moscow, Russia): ferroin (ferrous tris (1,10- [97]; analytical-grade sodium acetate; bismuth(III) hydroxide (high-purity grade); potassium iodide; ascorbic acid (pharmaceutical grade); sodium sulfite; copper(I) chloride; copper(II) acetate; lead(II) nitrate; sodium dihydroantimonate hexahydrate; tin(II) chloride; sulfuric acid (highpurity grade); nitric acid (analytical grade); sodium perchlorate; hydrochloric acid; analytical-grade potassium dichromate; analytical-grade acetic acid, 99.8%; a standard borate buffer solution, pH 8.50; chemically pure HCl (1 + 20); and chemically pure potassium hydroxide. Unless otherwise stated, all reagents were of cp grade and from Physchem 2023, 3, FOR PEER REVIEW 9 for two independent variables Ip(t = 0) and Ip(t = ¥), can be used to obtain ϒ = This results in an equation suitable for various thermal-lens modalities [95,98]: where δI is the error of measurements of probe light intensity by a photodetector, ( ) is determined by Equation (3), θ is determined by Equation (9), and ( ) = 2.303 ( )θ. For low steady-state signals, the descending reproducibility curve is described by a hyperbola, RSD = θ , where is a constant depending on the spectrometer optical scheme, and is a medium-dependent constant in the range 0.8-1.1 [111]. The results of thermal-lens measurements at different concentrations of test substances were obtained; 5-7 series of experiments were carried out (by changing the solution in the cell), and replicability, RSD rpl ( ), Equation (22), was determined for 3-5 series on different days (reproducibility, RSD , Equation (23)). The results obtained were subjected to variance analysis.
The results of thermal-lens measurements (cobalt complexes and ferroin) in various solvents and at different concentrations were obtained; 5-7 experiments were carried out by changing the solution in the cell, and replicability, RSD rpl ( ) , Equation (22), was determined for 20-30 consecutive measurements of the signal (each measurement, 10-15 min; repeatability, RSD , Equation (21)). The results obtained were subjected to variance analysis.

Performance Parameters
The measurement results are presented as mean values with confidence intervals (p < 0.05 or p = 0.95) under the requirements for the competence of testing and calibration laboratories in ISO/IEC 17025:2005. Calibration parameters were calculated according to IUPAC 1998 recommendations to present the results of the chemical analysis.
The limits of detection (LOD) and quantification (LOQ) were calculated as 3σ and 10σ criteria, respectively, according to IUPAC 1998 recommendations. The variances were compared using the Fisher, Bartlett, and Cochran criteria to determine the significance of contributions to the signal in control experiments.

Reagents and Chemicals
Water from a Milli-Q water purification system (Millipore SAS, Molsheim, France) was used: pH 6.8; specific resistance 18.2 MΩ × cm, Fe, 2 ppt; dissolved SiO2, 3 ppb; total ion amount <0.2 ppb; and TOC <10 ppb. The glassware was washed with acetone, followed by conc. nitric acid. All solutions were prepared daily before the experiments.
The following reagents were used throughout (unless otherwise stated, the products were from Reakhim, Moscow, Russia): ferroin (ferrous tris (1,10- [97]; analytical-grade sodium acetate; bismuth(III) hydroxide (high-purity grade); potassium iodide; ascorbic acid (pharmaceutical grade); sodium sulfite; copper(I) chloride; copper(II) acetate; lead(II) nitrate; sodium dihydroantimonate hexahydrate; tin(II) chloride; sulfuric acid (highpurity grade); nitric acid (analytical grade); sodium perchlorate; hydrochloric acid; analytical-grade potassium dichromate; analytical-grade acetic acid, 99.8%; a standard borate buffer solution, pH 8.50; chemically pure HCl (1 + 20); and chemically pure potassium hydroxide. Unless otherwise stated, all reagents were of cp grade and from I(t=∞) . This results in an equation suitable for various thermal-lens modalities [95,98]: where δ I is the error of measurements of probe light intensity I p by a photodetector, B(t) is determined by Equation (3), θ is determined by Equation (9), and B(t) = 2.303B(t)θ. For low steady-state signals, the descending reproducibility curve is described by a hyperbola, where k 1 is a constant depending on the spectrometer optical scheme, and k 2 is a medium-dependent constant in the range 0.8-1.1 [111].
The results of thermal-lens measurements at different concentrations of test substances were obtained; 5-7 series of experiments were carried out (by changing the solution in the cell), and replicability, RSD rpl (c), Equation (22), was determined for 3-5 series on different days (reproducibility, RSD rpd , Equation (23)). The results obtained were subjected to variance analysis.
The results of thermal-lens measurements (cobalt complexes and ferroin) in various solvents and at different concentrations were obtained; 5-7 experiments were carried out by changing the solution in the cell, and replicability, RSD rpl (c), Equation (22), was determined for 20-30 consecutive measurements of the signal (each measurement, 10-15 min; repeatability, RSD rep , Equation (21)). The results obtained were subjected to variance analysis.

Performance Parameters
The measurement results are presented as mean values with confidence intervals (p < 0.05 or p = 0.95) under the requirements for the competence of testing and calibration laboratories in ISO/IEC 17025:2005. Calibration parameters were calculated according to IUPAC 1998 recommendations to present the results of the chemical analysis. The limits of detection (LOD) and quantification (LOQ) were calculated as 3σ and 10σ criteria, respectively, according to IUPAC 1998 recommendations. The variances were compared using the Fisher, Bartlett, and Cochran criteria to determine the significance of contributions to the signal in control experiments.

Reagents and Chemicals
Water from a Milli-Q water purification system (Millipore SAS, Molsheim, France) was used: pH 6.8; specific resistance 18.2 MΩ × cm, Fe, 2 ppt; dissolved SiO 2 , 3 ppb; total ion amount < 0.2 ppb; and TOC < 10 ppb. The glassware was washed with acetone, followed by conc. nitric acid. All solutions were prepared daily before the experiments.
Aqueous preparations of hemoglobin (SPD RENAM, Moscow, Russia) with a concentration of 159 g/L and p.a. grade were used. Absorption spectra showed that the raw preparations of hemoglobin were mainly in the form of oxyhemoglobin, HbO 2 . Hemoglobin was transformed into MetHbCN via the interaction of oxyhemoglobin with a transforming solution (pH 7.3) containing potassium ferricyanide (200 mg/L), potassium cyanide (5 mg/L), potassium dihydrogen phosphate (150 mg/L), and Triton X-100 (3 mg/L). Absorption spectra within the range 300-700 nm compared to a solution containing the above reagents at the same concentrations were recorded.

Procedures
Procedure 1. Measurement and calculation of thermal parameters of aqueous colloidal solutions. A 3 mL sample of the solution of the test sample was placed in a quartz cell (path length, 10.00 mm), and the transient thermal-lens signal (λ = 532 nm) was measured at a room temperature of 20 • C. Next, 30-60 measurement cycles were averaged, and thermal diffusivity was calculated. Three aqueous solutions of ferroin at different concentrations (10-50 µM) were used to determine the thermal diffusivity of water.
Procedure 2. Measurement and calculation of thermal parameters of aqueous colloidal solutions of silicon oxide. Since silicon oxide practically does not absorb at 532 nm, when diluting the initial solutions of colloidal silicon, we added ferroin (10-50 µM). Next, 3 mL of the solution of the test sample was placed in the quartz cell (path length, 10.00 mm), and the transient thermal-lens signal (λ = 532 nm) was measured at 25 • C. Next, 30-60 measurement cycles were averaged, and thermal diffusivity was calculated.
Procedure 3. Measurement of albumin and its complexes with ferroin and CrEDTA. A small amount of water was placed in a 50 mL glass flask, and 450 mg of NaCl was dissolved. Then, 40 µL of 0.025 M ferroin solution was added and brought to the mark (20 µM). In a 50 mL glass flask, the prepared CrEDTA complex (0.002 M) was placed, 450 mg of NaCl (0.9%) was added, and the mixture was dissolved in a flask. Albumin solutions were prepared by diluting the original 10% solution with 0.9% NaCl.
Auxiliary procedures are given in the Supplementary Materials.

Selection of Samples
As protein dispersions, we used albumin, myoglobin, and cyanmethemoglobin (hemoglobin cyanide, MetHbCN). The determination of hemoglobin as MetHbCN (Drabkin method [112]) has been used for a long time and was approved by the International Committee for Standardization in Hematology (ICSH) in 1963, and this form is the most studied. Currently, this is the main method of hemoglobinometry [113]. The principle of the method is based on the conversion of all forms of hemoglobin into one, MetHbCN. Hemoglobin conversion is carried out when it interacts with a transforming solution containing ferrocyanide, cyanide, components of acid-base buffer systems, and a nonionic detergent.
To compare the results for proteins, we used aqueous dispersions of carbon nanomaterials, fullerenes and nanodiamonds (NDs), and silicon oxide nanoparticles. These carbon nanomaterials were selected as nanoparticles (clustered nanostructures) with their own light absorption properties, fullerene cages, and sp 2 Pandey chains in NDs [114][115][116]. We also used fullerene C 60 solutions in toluene (no aggregate formation, molecular solutions), in N-methyl-2-pyrrolidone (NMP) (in which C 60 aggregates), and in water (aggregates are formed, and the charge-transfer complex stabilizes them) [102].
Silicon oxide aqueous dispersions were selected, as they contain particles with increased thermal conductivity compared to water (heat-conducting nanofluids) and, at the same time, they are low-absorbing samples, being low-absorbing silicon oxide nanoparticles in a non-absorbing medium, water ( Figure S2, Supplementary Materials). The dominant contribution of Rayleigh light scattering to the absorbance recorded by the spectrophotometer is noticeable ( Figure S2, Supplementary Materials).
For comparison with dispersed systems and as reference points of the thermal-lens effect in true solutions, we used solutions of photostable and non-fluorescent dyes that form true molecular solutions that strongly absorb in the wavelength range of excitation lasers (450-540 nm): ferroin, gentian violet, and cobalt nitrate hexahydrate. Ferroin was used as a primary reference sample, as it is a chemically inert reference substance that does not form dispersed systems in water.
To assess the effect of electrolyte solutions on the sensitivity of the determination, photometric and thermal-lens measurements were carried out in saline solutions of chlorides, phosphates, and phosphate buffer saline, as well as in mixtures of concentrated chlorides (2.8 M NaCl + 0.4 M KCl). Surfactant-and polymer-modified aqueous solutions were represented by Triton X-100, Brij-35, Tween-80, and sodium dodecyl sulfate surfactants and polyethylene glycol polymer modifiers.

Thermal-Lens Setups
The back-synchronized measurement modality features different measurement conditions for the blooming and dissipation of the thermal lens [14]. The ADC-DAC unit continuously transfers the converted signal (voltage) from the signal photodiode to the measurement software. The shutter is closed at t = 0, and the probe intensity is maximum (Figure 2). At the measurement start, a clock pulse is sent to the shutter, the latter is opened, and the thermal field starts to bloom. The probe-beam intensity decreases due to the development of the diverging thermal lens, which is displayed as a probe-intensity vs. time chart. In the interval from t = t 1 to t = t 2 (the operator sets the trigger times of t 1 and t 2 ), the software reads the intensity I p (∞). At t 2 , the next clock trigger pulse is sent, the shutter gets closed, and the photothermally induced temperature field dissipates. In the interval from t = t 3 to t = t 4 , which is set up manually as well, the software reads the intensity I p (0). At t 4, the next measurement cycle is started. The software displays and stores the series, calculated by Equations (1) and (21).
To implement dual-beam coaxial measurements, we used dichroic mirrors to combine and to separate the probe beam from the excitation radiation. After the cell penetrated the dichroic mirror, the full intensity of the excitation beam was gathered with a focusing lens by the primary photometric detector (Figure 1), compared with the pre-calibrated power meter used as the absorbance ground signal. All three spectrometers implement a third optical channel ( Figure 1) for gathering scattering (or a luminescence signal, not used in this study), if present, at an angle of 90 • to the incident beam to compensate for the loss of the excitation laser power and the loss in the strength of the thermal-lens effect (Equations (11)-(13)).
The basic setup (TLS-60, Table 1) was built using a narrow excitation beam, with an excitation beam size of ca. 60 µm, to provide short thermal-lens blooming times (characteristic time in water, 6.2 ms) and high-sensitivity thermal-lens measurements. We also made two thermal-lens spectrometers with broad excitation beams (Table 1). TLS-150, with a beam size of 150 µm (characteristic time in water, ca. 40 ms), is more versatile and is designed both for assessing concentrations and for studying thermophysical or optical parameters. TLS-300, with a beam size of 300 µm (characteristic time in water, ca. 155 ms), has a 2-fold broader beam compared to TLS-150 and is designed to study the thermophysical parameters of strongly absorbing media alone because the sensitivity of light absorption is low (approximately 25-fold lower than for TLS-60). The instruments have different working wavelengths, which makes it possible to study different samples with the maximum sensitivity.
In photometry, scattered radiation enters the detector, increasing the signal, which leads to errors in the measurement of light absorption [117]. An approach for minimizing (but not eliminating) the scattered light effect uses high-fluence radiation, making it possible to move the detector to a longer distance from the sample and reduce the solid angle (the field of view) under which the detector is visible. A far-field thermal-lens spectrometer schematic automatically implements this feature ( Figure 1); thus, we added a photometer (transmission) modality to the setup (Appendix A), which accounts for absorbance data in thermal-lens measurements (Equations (10) and (11)).
The spectrometers do not have a single time constant, and irradiation times for an excitation on/off cycle vary from 0.05 to 200 s (depending on the measurement parameters, namely, on the data throughput rate and time and the number of points to be averaged). The spectrometers allow the readjustment of the optical-scheme geometry depending on the task by varying the positions of focusing lenses L1 and L2 ( Figure 1) and the distances between these lenses and the cell and the detector. This feature, along with varying irradiation times, provides an absorbed energy range from 0.1 mJ to 20 J and intensities in the sample cell from 250 to 5 × 10 7 W/cm 2 . The temperature increase due to the photothermal effect is 0.001-2 • C; thus, it is nondestructive and noninvasive.
We achieved rather low detection limits (1-2 orders of magnitude lower than for transmission photometry) for all spectrometers. For TLS-150, the instrumental LOD (IDL) for ferroin is 10 nM (300 mW excitation at 514.5 nm). The corresponding range of linear absorption coefficients for 10 mm optical pathways is from 1 × 10 −6 (300 mW excitation) to 2 × 10 −2 cm −1 (40 mW excitation). For TLS-150, the IDL for ferroin is 300 nM (at 400 mW excitation), which is an order higher than for common TLS setups [14], the linear range is slightly wider than two orders of magnitude (the corresponding range of linear absorption coefficients for 10 mm optical paths is 1 × 10 −5 -2 × 10 −3 cm −1 ).

Accuracy and Precision
The accuracy of steady-state thermal-lens measurements of light-absorption parameters was tested using model systems, ferroin, gentian violet, and cobalt nitrate hexahydrate, which are chemically inert and provide high absorptivities at the working wavelengths of excitation lasers. In all cases, the molar absorptivities of these substances have negligible differences between calculations using the thermal-lens signal (Equation (12)) and those using photometric measurements with the transmission channel of the spectrometer (Appendix A) and measurements of absorbance with a serial spectrophotometer. Transient thermal-lens curves for these systems are described by Equation (3), and the characteristic time of thermal-lens development (Equation (5)) was estimated ( Table 2) and compared to the theoretical calculations using the reference values for water, and a negligible difference was found.
To evaluate the accuracy and precision of thermophysical parameters, we tested true solutions of inert colorants and calculated the characteristic time thermal diffusivities of water and ethanol. In all of these cases, the characteristic time of the thermal lens (Equation (5)) did not differ significantly from the theoretical value (Table 2), and the theoretical curve (Equation (3)) fit the data with good reliability (Figure 3). Table 2. Performance parameters and characteristic times of thermal-lens determination of ferroin as a chelate dye and cyanmethemoglobin (MetHbCN) under various conditions in phosphate buffer solutions (pH 7.0) and phosphate buffer saline (pH 7.4), TLS-60; excitation wavelength, 514.5 nm; excitation power, 40.0 mW (p = 0.95, n = 8). Theoretical characteristic times of the thermal lens were calculated by Equation (5), and theoretical slopes were calculated by Equation (14).  (5), and theoretical slopes were calculated by Equation (14).   Table 3 shows that TLS provides good accuracy for thermal diffusivity valu silicone oxide nanoparticles; however, the accuracy of the thermal conductivity estim is lower, as the accuracy of measuring other parameters contributes to the final v Transient curves at the time of the attainment of the thermal equilibrium for 1-1 fullerene C60 aqueous dispersions (AFDs) fit the theoretical equations well, and we We used colloidal silicon oxide (silica) dispersions and ferroin. In this case, the medium is responsible for absorbing excitation radiation and nanoparticles for the rate of heat transfer. Using Equations (17), (19) and (20), we evaluated thermal diffusivity, thermal conductivity, and thermal effusivity ( Table 3). The thermal conductivity values of colloidal silicon oxide are consistent with the existing data under similar conditions for colloidal silicon oxide with a similar nanoparticle size [118,119]. Table 3. Experimental values of thermophysical parameters by thermal-lens spectrometry (n = 3, p = 0.95) and calculated values for colloidal silicon oxide solutions, AM (ω(SiO 2 ), 18% w/w), SM-30 (ω(SiO 2 ), 17.5% w/w), CL-X (ω(SiO 2 ), 26.5% w/w), TMA (ω(SiO 2 ), 20% w/w), and HS-40 (ω(SiO 2 ), 24% w/w); setup, TLS-60; excitation wavelength, 514.5 nm; excitation power, 40.0 mW.   Table 3 shows that TLS provides good accuracy for thermal diffusivity values for silicone oxide nanoparticles; however, the accuracy of the thermal conductivity estimation is lower, as the accuracy of measuring other parameters contributes to the final value. Transient curves at the time of the attainment of the thermal equilibrium for 1-10 µM fullerene C 60 aqueous dispersions (AFDs) fit the theoretical equations well, and we used the final part of the transient curves to estimate the bulk thermophysical properties of AFDs at the thermal equilibrium. Thermal diffusivities ((1.43 ± 0.03) × 10 −7 m 2 /s), thermal effusivities (1.55 ± 0.04) × 10 3 J/(m 2 ·K·s 1 2 ), and thermal conductivities for AFDs from thermal-lens data (with heat capacity and density obtained from other methods) show negligible changes compared to water, which is in rather good concordance with previous papers on fullerenes solubilized in water [120,121].

Sample
For all of the test systems, there is a significant difference between the components of the root-mean-square deviations of replicability (within a series, Equation (22)) and reproducibility (between series, Equation (23)) (Fisher's test, F exp > F theor ). There is an increase in the root-mean-square deviation of replicability only for the case of the background (water absorption). For all of the studied protein samples, RSD rep behaves similarly to water with inert colorants forming true solutions, with slightly higher values for signals over 0.1.

Hemoglobin
The determination of MetHbCN and ferroin shows low and similar LODs under moderate-power excitation conditions with the most sensitive setup TLS-60, as summarized in Table 2. The theoretical calculations of thermal-lens sensitivity for ferroin in water and various electrolyte solutions are in good agreement with the theoretical values from Equation (14) [122]. Compared to ferroin, the slope of the calibration plot of the thermal-lens signal vs. absorbance for MetHbCN in water is lower (Table 2 and Figure 4).
To prove that this effect results from the dispersive nature of MetHbCN solutions, we used mixtures of ferroin and MetHbCN; as the total thermal-lens signal is additive (as in photometry and Beer's law), the slope should also be an additive function of the molarities of the components of such binary mixtures. For mixtures in 1.0 M PBS (Table 2), the experimental results are in good concordance with this hypothesis. The slope for a mixture can be predicted as S mix = η ferr S ferr + η MetHbCN S MetHbCN , where S values are slopes for individual components, and η values are the molar fractions of ferroin and MetHbCN.
The determination of MetHbCN and ferroin shows low and similar LODs under moderate-power excitation conditions with the most sensitive setup TLS-60, as summarized in Table 2. The theoretical calculations of thermal-lens sensitivity for ferroin in water and various electrolyte solutions are in good agreement with the theoretical values from Equation (14) [122]. Compared to ferroin, the slope of the calibration plot of the thermal-lens signal vs. absorbance for MetHbCN in water is lower (Table 2 and Figure  4).  Measurements in saline solutions were intended, firstly, to show that a change in the environment leads to an increase in the sensitivity of the hemoglobin determination due to a change in thermophysical parameters [122] and, secondly, to find out how salts affect the signal linearity of the dispersed hemoglobin system. A mixture of sodium and potassium chlorides shows a significant increase in sensitivity, whereas phosphate solutions barely change the signal, and the mixture of four salts has the highest impact, which is not reducible to a mere sum of individual salt effects.
A comparison of the slopes of Equations (24) and (25) shows that there is a significant increase in the sensitivity of the determination of MetHbCN in saline media, and for ferroin, the sensitivity coefficient almost does not change compared to true solutions (Appendix C). The LOD of cyanmethemoglobin in saline media is 20 nmol/L (Table 2), which is 1.5 times lower than in the non-saline medium. The ferroin LOD does not change, which is consistent with the previous data for other saline solutions [122]. The ratio of sensitivity factors for cyanmethemoglobin and ferroin in a saline medium is 7.8 (Equations (24) and (25)), which agrees well with the value expected from the molar absorption coefficients. Thus, the influence of electrolytes is contrary to the thermal diffusion nature of such a difference in water [2], since thermal diffusion in electrolyte solutions should not be weaker than in distilled water. On the other hand, the effect of electrolytes on thermal conductivity and thermal effusivity is significant and may hide the weaker effect of dissolved chromophores. In general, the thermal-lens determination of cyanmethemoglobin in electrolyte solutions leads to a gain in the sensitivity coefficient and a decrease in detection limits.

Myoglobin
For myoglobin, the calibration slopes for TLS vs. transmission photometry (Equation (10)) for low powers of excitation are linear and barely change compared to ferroin as a true solution (Figure 4). For concentrations over 0.08% w/w, the concentration dependence of the thermal-lens signal depends on the excitation power in a more sophisticated way, and it changes from a nonlinear ascending curve to a curve with a maximum (Figure 5a

Bovine Serum Albumin
Albumins interact well with small molecules through electrostatic interactions. Thus, to change the properties of albumin, we used ferroin, which interacts well with this protein due to a pronounced positive charge, and the Cr(III) complex with EDTA (CrEDTA), which has a significantly lower charge density, i.e., it weakly interacts with albumin [124]. The modifying effect on the thermophysical parameters of albumin itself was estimated ( Figure S3, Supplementary Materials). It shows a pronounced linear dependence on the concentration of albumin. At the wavelengths of excitation lasers, albumin does not absorb, which means that the signal increases mainly due to changes in thermal parameters.
In the presence of ferroin, there was an increase in the albumin concentration compared to albumin solutions, leading to a decrease in the slope of the calibration curves ( Figure 6a). Moreover, the intersection of the calibration curves is in a narrow concentration range, proving that an increase in the ferroin concentration acts in the opposite direction to the albumin concentration. In the case of CrEDTA, we see linear dependences ( Figure 6b) with the same slope.

Bovine Serum Albumin
Albumins interact well with small molecules through electrostatic interactions. Thus, to change the properties of albumin, we used ferroin, which interacts well with this protein due to a pronounced positive charge, and the Cr(III) complex with EDTA (CrEDTA), which has a significantly lower charge density, i.e., it weakly interacts with albumin [124]. The modifying effect on the thermophysical parameters of albumin itself was estimated ( Figure S3, Supplementary Materials). It shows a pronounced linear dependence on the concentration of albumin. At the wavelengths of excitation lasers, albumin does not absorb, which means that the signal increases mainly due to changes in thermal parameters.
In the presence of ferroin, there was an increase in the albumin concentration compared to albumin solutions, leading to a decrease in the slope of the calibration curves ( Figure 6a). Moreover, the intersection of the calibration curves is in a narrow concentration range, proving that an increase in the ferroin concentration acts in the opposite direction to the albumin concentration. In the case of CrEDTA, we see linear dependences ( Figure 6b) with the same slope.
The ratio of the slope for a given concentration of albumin, E 0,alb , to the slope in water, E 0,w , shows a different character for ferroin, for which it shows a constant decrease (Figure 7a), which depends on the ferroin concentration (Figure 7b). For CrEDTA, this relative slope is somewhat similar to the plots for myoglobin, with a small increase followed by a decrease in the value (compare Figures 5 and 7). To further test albumin, we performed experiments with MetHbCN mixtures with albumin. For Hb concentrations corresponding to the middle of the calibration range (to decrease the overall error), the albumin addition (10-30 g/L) results in a decrease in the signal of more than 20% ( Figure S4, Supplementary Information). thermal parameters.
In the presence of ferroin, there was an increase in the albumin concentration compared to albumin solutions, leading to a decrease in the slope of the calibration curves ( Figure 6a). Moreover, the intersection of the calibration curves is in a narrow concentration range, proving that an increase in the ferroin concentration acts in the opposite direction to the albumin concentration. In the case of CrEDTA, we see linear dependences ( Figure 6b) with the same slope. The ratio of the slope for a given concentration of albumin, E0,alb, to the slope in water, E0,w, shows a different character for ferroin, for which it shows a constant decrease ( Figure  7a), which depends on the ferroin concentration (Figure 7b). For CrEDTA, this relative slope is somewhat similar to the plots for myoglobin, with a small increase followed by a decrease in the value (compare Figures 5 and 7). To further test albumin, we performed experiments with MetHbCN mixtures with albumin. For Hb concentrations corresponding to the middle of the calibration range (to decrease the overall error), the albumin addition (10-30 g/L) results in a decrease in the signal of more than 20% ( Figure  S4, Supplementary Information).

Transient Thermal-Lens Measurements
The transient curves for true solutions of ferroin and dispersed solutions of MetHbCN differ (Figure 8), and the characteristic times of thermal-lens measurements (Equation (5)) decrease significantly ( Table 2). This behavior depends on the concentration (Figure 9a).

Transient Thermal-Lens Measurements
The transient curves for true solutions of ferroin and dispersed solutions of MetHbCN differ (Figure 8), and the characteristic times of thermal-lens measurements (Equation (5)) decrease significantly ( Table 2). This behavior depends on the concentration (Figure 9a).  For MetHbCN and ferroin mixtures in saline solutions, the difference in characteristic time becomes smaller and even negligible (Figure 8b) compared to low ionic strengths. For myoglobin (Figure 8a), transient curves also differ from those of the ferroin solution, but the difference is much less compared to the same concentrations of MetHbCN, and the characteristic times of thermal-lens measurements (Equation (5)) decrease to a lesser degree than for hemoglobin for the same signals. The transient thermal-lens signal ( Figure  8c) for albumin solutions does not show any significant deviations from the shape for most test concentrations; for 50 g/L in the presence of ferroin, a decrease in the relative probe intensity is revealed, which is much less than for heme proteins, as shown above.
For silica nanoparticles and nanodiamonds (Figure 10), transient curves expectedly differ from those of aqueous ferroin solutions, much more strongly compared to MetHbCN at the same absorption level, and the characteristic times of thermal-lens measurements (Equation (5)) decrease, which results from an increase in thermal diffusivity. As shown previously, for a dispersed system, the first part of the curve, up to approximately 10 , corresponds to the appearance of a system of nano-thermal-lenses around the nanoparticles, which results in an apparent increase in the thermal diffusivity For MetHbCN and ferroin mixtures in saline solutions, the difference in characteristic time becomes smaller and even negligible (Figure 8b) compared to low ionic strengths. For myoglobin (Figure 8a), transient curves also differ from those of the ferroin solution, but the difference is much less compared to the same concentrations of MetHbCN, and the characteristic times of thermal-lens measurements (Equation (5)) decrease to a lesser degree than for hemoglobin for the same signals. The transient thermal-lens signal (Figure 8c) for albumin solutions does not show any significant deviations from the shape for most test concentrations; for 50 g/L in the presence of ferroin, a decrease in the relative probe intensity is revealed, which is much less than for heme proteins, as shown above.
For silica nanoparticles and nanodiamonds (Figure 10), transient curves expectedly differ from those of aqueous ferroin solutions, much more strongly compared to MetHbCN at the same absorption level, and the characteristic times of thermal-lens measurements (Equation (5)) decrease, which results from an increase in thermal diffusivity. As shown previously, for a dispersed system, the first part of the curve, up to approximately 10t c , corresponds to the appearance of a system of nano-thermal-lenses around the nanoparticles, which results in an apparent increase in the thermal diffusivity [48,102], and at times over 40t c , the course of the transient curve may be distorted due to thermal diffusion and convection. Thus, both light-absorption and thermophysical parameters were calculated for the range 10 ÷ 40t c ; the actual times differ for different setups. The cooling (dissipation) curves of the thermal-lens measurements were taken separately from the heating curves due to the different behavior of these curves, as the signal development is governed by a balance between the external source fluence and the sample's thermal properties, while dissipation depends on the stored energy and the solvent's thermal effusivity. The curves for proteins showed little difference compared to water (Figure 11), and in all cases, the difference is lower than the error of the experiments and provides the same thermal diffusivity. To the contrary, the dissipation curves at high concentrations of silicon oxide nanoparticles ( Figure 12) start to change with a steeper increase in the intensity due to the higher thermal diffusivities of these samples (Table 3). The cooling (dissipation) curves of the thermal-lens measurements were taken separately from the heating curves due to the different behavior of these curves, as the signal development is governed by a balance between the external source fluence and the sample's thermal properties, while dissipation depends on the stored energy and the solvent's thermal effusivity. The curves for proteins showed little difference compared to water (Figure 11), and in all cases, the difference is lower than the error of the experiments and provides the same thermal diffusivity. To the contrary, the dissipation curves at high concentrations of silicon oxide nanoparticles ( Figure 12) start to change with a steeper increase in the intensity due to the higher thermal diffusivities of these samples (Table 3). signal development is governed by a balance between the external source fluence and the sample's thermal properties, while dissipation depends on the stored energy and the solvent's thermal effusivity. The curves for proteins showed little difference compared to water (Figure 11), and in all cases, the difference is lower than the error of the experiments and provides the same thermal diffusivity. To the contrary, the dissipation curves at high concentrations of silicon oxide nanoparticles ( Figure 12) start to change with a steeper increase in the intensity due to the higher thermal diffusivities of these samples (Table 3).

Back-Synchronized Detection Modality
In TLS, mode-mismatched two-beam (two-laser) schematics are dominant [2,4,31,49,50,125,126]. It requires separate beams: a beam inducing the photothermal effect (the excitation beam) and a beam that is distorted in the photothermal refractive-index field and then probed (the probe beam). The term 'mode mismatching' means that the waists of the focused beams are shifted relative to one another along the propagation axis, which provides the maximum probe sensitivity of the thermal-lens effect [127][128][129][130]. In photothermal microscopy, mode mismatching is necessary, as this technique provides the necessary sensitivity of measurements [3,13,35]. Even for confocal schematics, very slight mode mismatching is needed to provide the probe signal [22,41,131,132].

Back-Synchronized Detection Modality
In TLS, mode-mismatched two-beam (two-laser) schematics are dominant [2,4,31,49,50,125,126]. It requires separate beams: a beam inducing the photothermal effect (the excitation beam) and a beam that is distorted in the photothermal refractive-index field and then probed (the probe beam). The term 'mode mismatching' means that the waists of the focused beams are shifted relative to one another along the propagation axis, which provides the maximum probe sensitivity of the thermal-lens effect [127][128][129][130]. In photothermal microscopy, mode mismatching is necessary, as this technique provides the necessary sensitivity of measurements [3,13,35]. Even for confocal schematics, very slight mode mismatching is needed to provide the probe signal [22,41,131,132].
Unlike lock-in amplifier schematics, the data acquisition and processing system in the back-synchronized mode registers a thermal-lens signal that is directly calculated from the detected intensities of the probe beam Thus, it does not need an empirical conversion (amplification) factor as in lock-in schematics, which makes the thermal-lens signal value more meaningful, since it relies on the assessment and calculation of factors determined by the environment [100], e.g., molar absorptivities (Equation (3)). Thus, the sensitivity of thermal-lens measurements can be compared by using the absolute values of the signals from Equation (3), without involving the ratio of the signal of the test sample and some reference point, which is needed in most lock-in schematics [133][134][135][136].
Thermal-lens spectrometers for high-sensitivity measurements of light absorption usually have a tightly focused excitation beam (waist size of ca. 50-60 µm [14,86] or even less in micro-photothermal techniques [21,22,41,106,[137][138][139][140][141][142]), which provides sufficient heating in the excited zone and, therefore, high optical detection sensitivity. However, such a schematic is not preferable for dispersed systems, especially when assessing the thermophysical parameters of the sample, because the heat flow from the heated zone depends on thermal conductivity, which may vary from the neat base solvents. Hence, the thermal-lens size becomes different for even a narrow concentration range, thus affecting the sensitivity and precision [2]. As a result, the geometry of the optical scheme for each sample should be changed to match the thermal lens and probe-beam sizes. Therefore, we designed spectrometers with a broad excitation beam (waist sizes, 150 and 300 µm), prolonging the time to attain a thermal equilibrium (Figure 1), and compared them to the more tightly focused beam in the 'reference' setup, TLS-60 ( Table 1).
The main parameter of the generic optimization of the optical scheme in thermal lensing is the desired diameter of the excitation beam. This determines the spatial locality (the thermal element diameter) and the sensitivity (power density) [2,130,136,143]. Experimental and theoretical approaches for optimization result in the same conclusion [2,98]. Hence, all other parameters of the spectrometer (frequency) should be selected so that the thermal element diameter is approximately threefold larger than the excitation, and the probe-beam diameter should be equal to the bloomed thermal element diameter. This will give the most extended linear calibration range and the lowest LODs and LOQs [2].
Optimizing a thermal-lens spectrometer involves changing the interposition of the elements of the optical circuit of the thermal-lens spectrometer so that the size of the excitation beam does not change, and the probe beam passes through the entire heated region and includes three steps. This is necessary for the maximum accuracy and sensitivity of measurements. If the probe beam is wider than the heated zone, interference occurs between the part of the beam that does not pass through the solution and the part of the beam that is deflected when passing through the thermal lens (Figure 13a), resulting in noise in the signal, because of which the sensitivity decreases [2,144]. When the probe beam is too narrow (ca. several micrometers), the sensitivity drops because it passes only through a part of the heated region (Figure 13c). Thus, the maximum sensitivity is observed when the probe beam passes through the entire heated area (Figure 13b) [2,144]. sensitivity of measurements. If the probe beam is wider than the heated zone, interference occurs between the part of the beam that does not pass through the solution and the part of the beam that is deflected when passing through the thermal lens (Figure 13a), resulting in noise in the signal, because of which the sensitivity decreases [2,144]. When the probe beam is too narrow (ca. several micrometers), the sensitivity drops because it passes only through a part of the heated region (Figure 13c). Thus, the maximum sensitivity is observed when the probe beam passes through the entire heated area (Figure 13b) [2,144]. The optimization of the optical scheme made it possible to achieve the constancy of the geometric parameter B for thermal-lens spectrometers. We achieved the maximum width of the linearity range of these dependences ( Table 2). Thermal-lens spectrometers have been shown to work well with aqueous solutions and worse with organic media; this is due to the strength of the thermal-lens effect in non-polar solvents, which, for TLS-150, Figure 13. Schematic representation of excitation (boundaries are green lines) and probe beams (boundaries are red lines; dashed lines denote the boundaries without a thermal-lens effect) through the sample in a thermal-lens experiment. The red-to-white gradient denotes the thermally induced field in the sample (the thermal-lens effect). (a) The probe beam is too broad, and only the central part broadens due to the thermal-lens effect and interferes with the undisturbed part at the detector plane (not shown here for simplicity's sake); (b) the correct size of the probe beam, so it is almost equal to the size of the photothermal field (thermal lens); and (c) the probe beam is too narrow, and it does not probe the whole thermal field, degrading the sensitivity of measurements.
The optimization of the optical scheme made it possible to achieve the constancy of the geometric parameter B for thermal-lens spectrometers. We achieved the maximum width of the linearity range of these dependences ( Table 2). Thermal-lens spectrometers have been shown to work well with aqueous solutions and worse with organic media; this is due to the strength of the thermal-lens effect in non-polar solvents, which, for TLS-150, led to the quadratic dependence of the signal on the substance concentration, which is due to a change in the geometric parameters of the thermal lens and the deviation of the thermal-lens model (Equations (1)- (3) and (15)) with concentration [2,104]. There was also a decrease in the determination sensitivity for TLS-300 compared to TLS-150 due to beam broadening; however, this made it possible to work with toluene solutions, while TLS-150 does not allow this.

Transient and Steady-State Measurements
In addition to steady-state signals, which are mainly applied in thermal-lens measurements of light absorption or chemical analyses (concentration measurements), the cycle of the formation and dissipation of a thermal lens is characterized by periods of fall and increase in the intensity of the probe radiation (Figure 2), while transient thermal-lens signals are significantly more dependent on the thermophysical parameters of the medium than on light absorption (somewhat in contrast to steady-state signals) [2]. The transient thermal-lens curve depends on the thermal conductivity, heat capacity, thermal diffusivity, and dn/dT of the sample, absorption, and instrument geometry (Equation (3)). It should be noted that the difference in the transient development and dissipation curves of the thermal lens can be detected using the back-synchronized thermal-lens spectrometer because the thermal equilibrium is attained at every stage.
In true solutions, the absorbed energy is transformed after laser irradiation, which depends only on the energy distribution in the excitation beam. As a result, a gradient of the refractive index, the thermal lens, is formed. Due to the presence of large light-absorbing dispersed particles, part of the energy goes toward heating the particles themselves, and the heat transfer to the solution decreases, which leads to a noticeably weaker total thermal field in the sample irradiated by the laser [145]. Thus, for dispersed systems, the local heating at the beginning is located around particles with a high heat capacity, followed by the transfer of heat to the solvent, resulting in significantly more rapid heating and a reduction in the thermal-lens characteristic time. On the contrary, thermal-lens dissipation exhibits smaller changes compared to homogeneous samples because, for small temperature changes (10 mK, typical for the concentrations of most absorbing compounds in TLS), the dispersity of the medium provides a smaller impact on heat dissipation [2]. According to previous findings and models [41,96,108,146], this behavior (Figures 4-6) is due to a higher heating rate near strongly light-absorbing dispersed particles, which, therefore, decrease the apparent thermal-lens characteristic time. Meanwhile, the heat transfer into the aqueous dispersion is slower [93,108], which leads to more heat dissipation during the development (blooming) of the thermal lens. The average distance between the globules for this concentration range is 30-300 nm, comparable to the Hb globule size, and should affect the thermal properties. If the absorbance of solutions is the same, this effect decreases the steady-state thermal-lens signal from a dispersed system compared to a true solution without a linearity loss (Figure 4).

Measurement Conditions
The microheterogeneity of the medium, e.g., the presence of a surfactant in an aqueous solution, usually increases the thermal-lens signal, and with an increase in its concentration, the signal expectedly increases due to changes in thermal parameters [147][148][149]. However, special attention should be paid to the conditions (excitation power and time selection) of a thermal-lens experiment in dispersed media, particularly for steady-state (stationary) measurements [84,97,150]. Firstly, in two-component media, the role of the temperature gradient in the solution is twofold. It forms both a refractive-index profile (the thermal lens) and a concentration profile of solute molecules in solution due to thermal diffusion (thermodiffusion, the Soret effect) [28,[151][152][153]. This concentration spatial profile may function as another lens-like optical element contributing to probe-beam blooming, together with the thermal-lens action. Thus, the signal (positive or negative) due to the Soret effect may mistakenly be interpreted as an increase or decrease in the thermal-lens signal, which may lead to errors in assessing a light-absorbing analyte [2]. To separate the effects of temperature and concentration gradients, it is necessary to use special techniques, e.g., time selection of the signal [28,154], or the settings of the parameters to efficiently evaluate and exclude the accompanying thermally induced diffusion processes.
Secondly, thermal-lens experiments are accompanied by convective heat transfer, which decreases the thermal-lens strength and disturbs the photothermal effect, with periodic fluctuations resulting in increased measurement errors, i.e., both bias and random errors of photothermal measurements [155][156][157][158]. The periodic nature of thermal-lens measurements facilitates convection and provides the possibility of studying convective phenomena. However, in most tasks of thermal-lens measurements (thermophysical or light-absorption/concentration parameters), convection is a degradation factor for both accuracy and precision, and it should be suppressed if possible. Back-synchronized thermal lensing, by default, uses different times for steady-state heating, and full dissipation usually suppresses convective effects in thermal lensing.
Thus, time selection in thermal lensing is of prime importance because we need to select the time of blooming and the dissipation of the thermal lens to provide (i) the maximum sensitivity of the thermal-lens effect and (ii) the maximum accuracy and precision of the measurements. The second problem depends on the information we need to obtain from the experiments: thermal diffusivity requires short times and the analysis of the transient curves, while absorption measurements require steady-state measurements.
It is worth mentioning that as the thermal-lens effect is scaled depending on the fluence (and the excitation beam size) and thermal diffusivity of the base medium, the characteristic time (Equation (5)) is a scale parameter for time selection. Thus, to generalize the problem, it is expedient to use the relative time scale of measurements, i.e., t/t c values ( Figure 14). Using the example of the most narrowly focused TLS-60 setup (characteristic time in water is 6.1 ms, Table 2), Figure 14 (a single, non-averaged curve to more explicitly show the changing error) shows that in time curves in the interval of steady-state measurements of the heating part (Figure 2), the time (0.5 ÷ 5) × t c corresponds to the transient process (marked by the yellow rectangle), and (5 ÷ 20) × t c corresponds to approaching the steady state (the green rectangle). Longer times correspond to steady-state measurements, and t 1 + t 2 = (20 ÷ 40) × t c provides steady-state measurements with the minimum error, while at longer values (the red rectangle), some other effects start to appear (the Soret effect). As a result, up to 20t c (ca. 125 ms), the curves for water reasonably fit the theoretical equation (Equation (3)). At 100 ms, the time curves approach the plateau, which is in good concordance with the theory. Physchem 2023, 3, FOR PEER REVIEW 24 thermal diffusivity would result in an overestimation. As a result, the straightforward use of the linear part of the cooling transient curve at (0.5 ÷ 10) for dispersed media seems incorrect due to a positive bias, and the use of longer times seems more expedient. However, long times may result in the appearance of conductive and convective heat transfer [158] and also seem inexpedient. Thus, the time selected for data handling depends on the task: for light-absorption parameters, the plateaus should be selected with the maximum number of points to decrease the error, while for thermal parameters, the early stages corresponding to nanolens formation in dispersed media [106] should be excluded, as well as long heating times, which may involve conductive and convective heat transfer as well as thermal diffusion effects.

Reproducibility and Replicability
The first precision parameter is reproducibility, which, according to ISO 5725-1, specifies the closeness of agreement between independent test results obtained with the same method on identical test material but under different conditions (different apparatuses, different operators and laboratories, and a considerable time interval between measurements). Previously, we showed that the reproducibility of thermal-lens measurements in true solutions [111], aqueous surfactant solutions [159], and polymerunmodified aqueous solutions [97] does not differ from that in pure water and mainly depends on the error in sample preparation procedures. This is supported by other findings in thermal lensing [11,88,101,160] and photothermal-lens microscopy [11,91,161,162].
A significant difference between the components of the root-mean-square deviations of replicability (Equation (22)) and reproducibility (Equation (23)) confirms our previous findings [100] that the main contribution to the error is from reproducibility (procedures and protocols for sample preparation), while the repeatability precision of thermal-lens measurements is almost independent of the test substance and its concentration and depends on the thermal-lens setup only and its instrumental optimization. For proteins, the experimental curve seems to fit the theoretical prediction well (Equation (23)).

Repeatability
Another precision factor is repeatability, which, according to ISO 5725-1, specifies the closeness between results obtained under the same experimental conditions (the same It should be noted that t c is the heating characteristic time, and it could be used to characterize the dissipation curve with caution. The physical meaning of t c , the halfblooming of the probe beam, is the relative decrease in intensity, which is 0.5 ( Figure 2) compared to the steady state. In the cooling half-cycle, the relative decrease in intensity of 0.5 has no such meaning ( Figure 14). However, as t c is the generic time scale of a thermal-lens experiment, the normalized time scale seems appropriate for characterizing the dissipation of the thermal lens. Thus, for the cooling period of the excitation on/off cycle, the transient period is the same as for the heating period, (0.5 ÷ 7)t c , and the period of approaching the steady state is longer, (7 ÷ 30)t c . The steady state can be measured in the time interval t 3 + t 4 = (30 ÷ 50)t c or for even longer times.
However, the accuracy of measurements may be affected by the deviation of the transient curve from Equation (3) in dispersed media due to the nanolens formation at early stages [106]. Thus, for nanodiamonds and colloidal silica particles (Figure 10), a steeper decrease in the relative probe-beam intensity means decreasing the apparent characteristic time. However, the direct application of Equation (5) for the calculation of thermal diffusivity would result in an overestimation. As a result, the straightforward use of the linear part of the cooling transient curve at (0.5 ÷ 10)t c for dispersed media seems incorrect due to a positive bias, and the use of longer times seems more expedient. However, long times may result in the appearance of conductive and convective heat transfer [158] and also seem inexpedient. Thus, the time selected for data handling depends on the task: for light-absorption parameters, the plateaus should be selected with the maximum number of points to decrease the error, while for thermal parameters, the early stages corresponding to nanolens formation in dispersed media [106] should be excluded, as well as long heating times, which may involve conductive and convective heat transfer as well as thermal diffusion effects.

Reproducibility and Replicability
The first precision parameter is reproducibility, which, according to ISO 5725-1, specifies the closeness of agreement between independent test results obtained with the same method on identical test material but under different conditions (different apparatuses, different operators and laboratories, and a considerable time interval between measurements). Previously, we showed that the reproducibility of thermal-lens measurements in true solutions [111], aqueous surfactant solutions [159], and polymer-unmodified aqueous solutions [97] does not differ from that in pure water and mainly depends on the error in sample preparation procedures. This is supported by other findings in thermal lensing [11,88,101,160] and photothermal-lens microscopy [11,91,161,162].
A significant difference between the components of the root-mean-square deviations of replicability (Equation (22)) and reproducibility (Equation (23)) confirms our previous findings [100] that the main contribution to the error is from reproducibility (procedures and protocols for sample preparation), while the repeatability precision of thermal-lens measurements is almost independent of the test substance and its concentration and depends on the thermal-lens setup only and its instrumental optimization. For proteins, the experimental curve seems to fit the theoretical prediction well (Equation (23)).

Repeatability
Another precision factor is repeatability, which, according to ISO 5725-1, specifies the closeness between results obtained under the same experimental conditions (the same operator, laboratory, and apparatus, with a short interval between measurements). This characteristic is less significant in a chemical analysis in biomedical studies because it mainly characterizes an instrumental (meter) error.
In thermal-lens spectrometry, this parameter is represented by short-term fluctuations of the signal from excitation on/off cycles. The contribution to the repeatability RSD from the instrument error (laser power fluctuations and jitter of optical elements) is significant but is much lower than the solvent effect (Equations (3) and (9)). The most curious fact is that RSD in some solvents depends much less on the sample absorbance or the excitation power in some signal ranges [93,111]. In other words, a statistical parameter acquires the meaning of the thermo-optical constant of the sample, or at least of a signature parameter. Repeatability depends on thermal diffusivity, non-equilibrium heat transfer (including nanolens formation), convection, and the Soret effect. Still, the separate contributions of thermophysical processes are difficult to find directly. Therefore, a repeatability comparison for two samples (all other conditions being the same) provides the estimation for thermallens detection limits and sensitivity coefficients, and this works for different solvents and additives in solutions [163].
Here, we considered repeatability parameters for the studied proteins and nanoparticle solutions. For pure water with inert colorants, RSD rep (Equation (21)) vs. the steady-state thermal-lens signal (Equation (7)) ( Figure 15a) is similar to the low-absorption branch of the reproducibility curve ( Figure S1, Supplementary Materials [101]), as the instrumental meter error decreases with an increase in the light absorption for absorbances low enough to ensure no significant deviation from the thermal-lens model [2,104,164,165]. For proteins, the actual values increase (Figure 15a), probably due to a change in the size of the thermal lens in these media. Such independence of the relative standard deviation from the concentration (Figure 15a) was also found for surfactants of various origins within the range 1-10 mmol/L. This is in good agreement with the existing data showing that, for these concentrations of nonionic surfactants, the water structure is not significantly changed owing to insignificant changes in the hydrogen-bond structure [147].
of rather finely dispersed solutions of these species. This is confirmed by a gradual increase in RSDrep prior to changes in the thermal-lens signal in the monitoring of spontaneous precipitation [169]. Despite multiple manifestations of this effect in thermallens experiments, its description is still qualitative, and its use is limited. The differences in RSDrep for proteins found in this study are slight compared to true solutions of model systems, which made it impossible to use them for characterization. Still, the use of RSDrep as a signature parameter in thermal-lens studies seems relevant for the further development of the method. Most previous studies estimated the repeatability of steady-state measurements of lock-in measurements with a time constant. In this study, we estimated RSDrep for transient measurements as a function of the transient thermal-lens signal (Equation (3)). By deconvoluting a transient thermal-lens curve, it is possible to describe the thermal equilibrium and to assess the early stages of thermal-lens development that may depend more on the dispersed-phase parameters. As expected, the data corresponding to approaching the thermal equilibrium ( Figure 14) show the same good precision as steadystate signals, while the estimation of RSDrep for earlier stages (0.5 ÷ 10 × ) did not show any valuable data.

Hemoglobin and Myoglobin
Changes in the characteristic time for hemoglobin solutions (Figures 8 and 9) compared to true solutions are more distinct than changes in the slope of a concentration calibration ( Figure 4) and are comparable to those for nanoparticle solutions ( Figure 10). However, as surfactants affect the thermal properties of water much more, the following changes in RSD rep curves for nonionic surfactants and SDS most likely result from a change in the thermal properties of the aqueous solution caused by an increase in the number of micelles rather than in their size. The growth of fluctuations for higher surfactant concentrations (over 100 mmol/L) is probably caused by the increasing non-homogeneity of dispersed solutions and the contribution of light scattering, which can be accounted for by a change in the micelle shape from spherical to cylindrical as well as their aggregation [166]. In Triton X-100, fluctuations grow and reach a plateau at 7 mmol/L. Then, in the range of 0.007-0.15 mol/L Triton X-100, the fluctuations are twofold higher than in pure water (20-25%). For Brij-35, the fluctuations are virtually constant up to 50 mmol/L ( Figure 15a) and start increasing again with a further increase in the surfactant concentration. Thus, the effect of surfactants on the thermal-lens repeatability in aqueous media differs significantly from that of non-surfactant substances and should be taken into account, e.g., when inert modifiers are used along with proteins in light-absorption measurements [99].
It is worth mentioning that dispersed phases also affect the thermal-lens measurements and can be used not only as a measure of thermal properties but also as the signal source. The interference from poorly soluble compounds of metals on the thermal-lens determination of bismuth [167] manifests itself in the relative standard deviation of the reagent mixture (blank) and the analyte signals (Figure 15b). The blank signal does not depend on the concentrations of certain amounts of other metals (without bismuth), and RSD rep also does not change, which is shown as the horizontal parts of the curves in a rather broad range of the molar ratios of the interferent metal and bismuth ions. However, starting from a certain concentration of the interferent metal, RSD rep is increased by a factor of 2-20, as represented by the increasing parts of the curves (Figure 15b). The metal concentrations corresponding to this increase in RSD rep are much higher than the solubility constants of corresponding iodides or oxides that can form under these conditions [168]. Thus, it could be said that the thermal-lens effect reacts to the formation of rather finely dispersed solutions of these species. This is confirmed by a gradual increase in RSD rep prior to changes in the thermal-lens signal in the monitoring of spontaneous precipitation [169]. Despite Physchem 2023, 3 182 multiple manifestations of this effect in thermal-lens experiments, its description is still qualitative, and its use is limited. The differences in RSD rep for proteins found in this study are slight compared to true solutions of model systems, which made it impossible to use them for characterization. Still, the use of RSD rep as a signature parameter in thermal-lens studies seems relevant for the further development of the method.
Most previous studies estimated the repeatability of steady-state measurements of lock-in measurements with a time constant. In this study, we estimated RSD rep for transient measurements as a function of the transient thermal-lens signal (Equation (3)). By deconvoluting a transient thermal-lens curve, it is possible to describe the thermal equilibrium and to assess the early stages of thermal-lens development that may depend more on the dispersed-phase parameters. As expected, the data corresponding to approaching the thermal equilibrium ( Figure 14) show the same good precision as steady-state signals, while the estimation of RSD rep for earlier stages (0.5 ÷ 10 × t c ) did not show any valuable data.

Hemoglobin and Myoglobin
Changes in the characteristic time t c for hemoglobin solutions (Figures 8 and 9) compared to true solutions are more distinct than changes in the slope of a concentration calibration ( Figure 4) and are comparable to those for nanoparticle solutions ( Figure 10). The characteristic time in MetHbCN solutions depends on the protein concentration ( Figure 9); the t c value for ferroin-MetHbCN mixtures is additive (Table 2 and Figure 9b). An equimolar mixture of MetHbCN and ferroin shows a shift in t c , as expected from component molar fractions. This model for thermal lensing in MetHbCN solutions is well simulated by applying the multipoint heat-source theoretical approach (Appendix B) [95,96,108]. This effect of Hb globules with high local light absorption (compared to the much less absorbing solution around them) is responsible for local hot zones previously reported for RBCs [170,171].
The approach based on thermal-lens development as an ensemble of point sources predicts this in a semiquantitative way [95,96]; however, a full solution requires further developments in theory. Moreover, apart from the change in the conductive heat-transfer mechanism in dispersed solutions, such a significant difference in the sensitivity coefficient and characteristic time can result from convection heat transfer [143,158,172,173] and the Soret effect [28,[151][152][153]. However, in such cases, the transient curves would be poorly described by the mono-exponential dependence characteristic of a purely thermal-lens effect and described much better by bi-exponential dependences corresponding to concurrent conductive and convective or thermal and concentration lenses.
If the hypothesis about the influence of dissolved hemoglobin on the thermophysical properties of the medium is correct, then the development and dissipation curves of the thermal lens should remain close to mono-exponential, but the characteristic times of the development of the thermal lens should differ from those of a purely aqueous medium, as discussed above. In addition, the characteristic time of thermal-lens development and dissipation should differ to a greater extent than in a water environment due to the different effects of the solute (in this case, hemoglobin) on the thermal conductivity and thermal effusivity of the medium. To verify this, transient development curves of the thermal-lens signal were recorded at different concentrations of MetHbCN, ferroin, and their mixtures in aqueous and saline media (Figures 8 and 9). Based on the selection of measurement parameters (Section 4.3), we ensure the achievement of a stationary thermal state.
For thermal-lens heating and cooling curves, the functions y = ae −bx and y = a(1 − e −bx ), respectively, gave the worst approximations, with correlation coefficients less than 0.7 in all cases, so they were not used further. For other functions, y = y 0 + ae −bx for a purely thermal-lens effect and y = y 0 + ae −bx + f e −gx , representing two competitive mechanisms, the correlation coefficients were not lower than 0.97 (Table S2, Supplementary Materials). Thus, the thermal-lens development curve is well described by a shifted mono-exponent. It indicates a weak Soret effect or a slight contribution from convective heat transfer to the transient signal under the selected measurement conditions and test samples.
The medium changes when adding electrolytes that affect the sensitivity ( Table 2). As expected, the slope changes for ferroin follow the electrolyte thermal diffusivity and the temperature gradient of the refractive index [122]. However, the changes in the calibration slope for MetHbCN are more drastic and do not directly depend on the ionic strength but rather on the electrolyte nature, with an additional increase in PBS compared to chloride solutions only (Table 2). Additionally, for higher ionic strengths, the difference in the thermal-lens characteristic time between true solutions and protein dispersions becomes smaller and even negligible (Figure 8), and the error in the estimation of changes in the characteristic time using Equation (5) increases significantly. The multipoint response calculation (Appendix B) results in a smooth spatial profile with a higher temperature in the beam center, which correlates entirely with the increase in the slope and characteristic time for MetHbCN (Table 2). Thus, transient curves and characteristic times for various electrolyte-based media containing MetHbCN showed concordance with a decrease in the slope of the calibration θ = f (A) (Equation (10)) responsible for bulk thermal properties and is in good concordance with the result of a multipoint light-absorbing model of thermal-lens measurements [96,108].
A comparison of the slopes of the curves shows a significant increase in the sensitivity of the determination of MetHbCN in saline media, and for ferroin, the sensitivity coefficient barely changes. The detection limit of MetHbCN in the saline medium is 20 nmol/L ( Table 2), which is 1.5 times lower than in the non-saline medium. The ferroin detection limit does not change, which is consistent with previous data for other saline media [122]. The ratio of sensitivity factors for both MetHbCN and ferroin in the case of a saline medium is 7.8, which agrees well with the value expected from the molar absorptivities for these substances. Thus, the impact of electrolytes eliminates the difference between the development of the thermal-lens effect for MetHbCN and ferroin, which is contrary to the thermal diffusion nature of such a difference [2], since thermal diffusion in electrolyte solutions should be no weaker than in distilled water. On the other hand, the effect of electrolytes on the parameters of the medium (thermal conductivity and thermal effusivity) is significant and may hide the weaker effect of dissolved chromophores. In general, the thermal-lens determination of MetHbCN in electrolyte solutions leads to a gain in the sensitivity coefficient and a favorable decrease in detection limits.
In the case of myoglobin, all of the effects discussed in the previous section for hemoglobin are much less pronounced (Figure 8 for transient and Figure 4 for steady-state measurements). This correlates well with the point-source approach's profile calculation. Very characteristic is the slope dependence and the whole behavior of the concentration dependence of myoglobin on laser power ( Figure 5), which shows that local overheating due to nanolens formation depends on the power and fluence. This behavior is remarkably similar to the local overheating and so-called nonlinear mode in photothermal microscopy [21,41,123]. In the latter case, it is due to the formation of nanobubbles in hot zones around the nanostructures, which increases the sensitivity of imaging and quantification. In our case of a rather slow thermal-lens effect and steady-state thermal lensing, the result is quite the opposite: a decrease in the sensitivity down to a decrease in the signal with the chromophore concentration. In our opinion, such behavior can be readily seen with myoglobin, as it allows high excitation powers to be used with no decomposition of this protein.

Bovine Serum Albumin
Contrary to heme proteins, the transient thermal-lens signal for albumin solutions (Figure 8c) does not show any significant deviations in either the shape or the characteristic time. It confirms the absence of thermal diffusion in the solution, which further confirms the thermal character, rather than being concentration-based, of the signal decrease in this case. Figure S3 (Supplementary Information) shows the pronounced linear dependence of the signal on the albumin concentration. Such an increase in the signal of aqueous albumin solutions shows its modifying effect on the thermal parameters. Spectrophotometric measurements showed that in the selected wavelength range, albumin does not absorb, which means that the signal increases due to the photothermal parameters of the solution.
Experiments with ferroin and CrEDTA chelate were used to further prove the effects of point-source absorption on the response of the system under photothermal heating. Ferroin, which, as expected, interacts with albumin to a much higher degree than CrEDTA, is similar to hemoglobin: light-absorbing nanostructures (albumin with ferroin) in a nonabsorbing medium. While albumin itself does not absorb at the excitation wavelength and thus cannot be a heat source, the ferroin-albumin complex leads to the formation of light-absorbing point sources that transfer less heat to the solution compared to aqueous ferroin, and it leads to local overheating. Thus, from the viewpoint of light absorption and photothermal heating, we have a dispersed medium with excess ferroin, with a part of it, upon interaction, being attached to albumin globules. Thus, we have almost the same situation as with MetHbCN, where some of the light-absorption centers are point heat sources in photothermal heating. Thus, the more albumin added (Figure 7a), the greater the fraction of light absorption of point heat sources, and the lower the final signal. If we consider the ferroin concentration, we have a confirmation of this: the higher the ferroin concentration (more point heat sources formed), the steeper the slope (Figure 7b), while the total signal increases due to the higher total absorption of ferroin.
In the case of the CrEDTA complex with a weaker interaction between albumin and CrEDTA, the calibration slopes ( Figure 6b) do not differ significantly from an unmodified aqueous medium. However, in general, it has a similar behavior (Figure 7a) to myoglobin ( Figure 5). This is due to competitive effects: an increase in the signal due to changes in the thermal parameters of the solution upon an increase in the albumin concentration and a decrease in the signal due to point sources on the transient curves. As the pointsource effect is much slower, it is manifested as a lower decrease compared to ferroin and at higher concentrations. As a proof-of-concept, we used samples of albumin with cyanmethemoglobin. These proteins can interact, forming a larger nanostructure [174]. Thus, this effect should be revealed in thermal lensing under selected conditions. As in previous cases, photometric measurements did not show a significant change in absorbance between MetHbCN and MetHbCN/albumin solutions. However, thermal lensing shows a ca. 20% decrease compared to MetHbCN ( Figure S4, Supplementary Information), which is higher than the effect of ferroin/albumin mixtures. Within the framework of this study, this was intended as a demonstration experiment, which proved the interaction effect on thermal lensing, and further studies are needed for a more detailed description.

Comparison with Engineered Nanoparticles
Thus, heating curves for all of the studied proteins show changes compared to aqueous solutions without proteins, while dissipation curves show much smaller changes in their behavior compared to heating curves (compare Figure 11 vs. Figures 8 and 9). This is accounted for by a very slight change in the thermal diffusivity of aqueous protein solutions compared to solutions of nanoparticles and fullerene/nanodiamond aggregate structures ( Figure 11), which, in the case of silica nanoparticles, are heat-conducting nanofluids ( Table 3). This effect ensures more profound nanolens-driven transient curves for nanoparticles. For nanodiamonds, the heating curves seem to show the same situation, lowering t c with the nanoparticle concentration. However, as we have shown previously [106,119], nanodiamonds are rather poor heat-conducting nanofluids, as photothermal heating accumulates the heat during the heating period of the photothermal-excitation cycle. This results in the reversed behavior of the cooling curve ( Figure 11): the dissipation becomes slower due to the action of secondary heating sources. Thus, the cooling curve for a colloidal solution cannot be approximated by the classical temperature equation for a homogeneous medium (Equation (18)), which is experimentally confirmed in the case of water. As the studied protein solutions show overheating and no changes in cooling curves, this shows their somewhat insignificant effect on thermal parameters.
The whole comparison of the results for all three studied proteins and nanoparticle systems allows us to make the following methodological conclusions for measurements in dispersed media. First, tests for accuracy and precision with true and dispersed solutions are demonstrated with known protocols of production (gold, silver, or silica). As protein solution properties may vary due to reaction conditions, prior to thermal lensing, calibration with nanoparticles with different but well-established properties should be performed. In our opinion, the silica nanoparticles used in this study or metal or metal-oxide particles forming heat-conducting nanofluids can be used for this purpose. The advantage of metal or oxide nanoparticles is the possibility of obtaining thermal diffusivity or thermal effusivity values from alternative methods [106,119]. The use of nanodiamonds as calibration samples seems inexpedient due to various properties depending on the nanodiamond brand and the uncommon behavior of the cooling curve [106,107,143,145]. Certainly, if a more precise model for the thermal properties of aqueous dispersions of nanodiamonds (as heat-conducting and simultaneously heat-accumulating solutions) is presented, nanodiamonds can be used as well. Contrary to nanodiamonds, fullerenes may be an expedient candidate for calibration, as they mildly alter the bulk thermal properties compared to nanofluids, and their nanostructural parameters in solution change from true to slightly heat-conducting aggregates depending on the solvent and environment [107,143,145].

Conclusions
This study shows that transient and steady-state thermal lensing should not be performed alone for samples with unknown absorption/concentration and thermal properties, and they should be used together. This study shows that the simultaneous estimation of thermal and optical/concentration parameters of dispersed systems should provide (i) the maximum sensitivity of the thermal-lens effect (absorption measurements) and (ii) the maximum accuracy as well as precision of the measurements (thermal parameters). Thus, time selection for thermal lensing (stages of transient blooming, transient dissipation, and steady states) becomes the key stage, as assessing thermal diffusivity requires the analysis of transient curves and relatively short excitation times, while absorption parameters necessitate steady-state measurements with longer times and prolonged signal accumulation. Therefore, the time scale (thermal-lens characteristic time, t c ) should be selected at the stage of the thermal-lens setup design, and the relative time intervals should be selected: (0.5 ÷ 5)t c for transient and (5 ÷ 20)t c for steady-state measurements. These estimations can be treated as generic for aqueous dispersed systems, as various substances were used (proteins, surfactants, inorganic iodides, and silica oxide and carbon nanoparticles). Still, these parameters should be studied for other dispersed systems and base solvents.
Thus, a thermal-lens spectrometer for dispersed solutions could be implemented as a multifunctional apparatus upon the selection of measurement parameters governed by both accuracy and precision. The balance between the precision of measurements of thermal and optical parameters corresponds to excitation beam sizes of 60-300 µm; a broad excitation beam provides longer times to attain a thermal equilibrium and, thus, the better precision of measurements of thermal diffusivity, though increasing the measurement time and decreasing the accuracy of steady-state measurements. The setups fabricated in this study provide good balance and flexibility for both measuring light-absorption parameters and assessing the concentrations of dispersed solutions with trace levels of constituents and assessing the thermal properties of dispersions. Despite a decrease in the sensitivity of light-absorption measurements of lock-in schematics, the setups retain the features of the sensitivity and instrumental/data-handling simplicity of common dual-beam thermal lensing. In our opinion, an important feature of the back-synchronized schematics used in this study is the possibility of selecting and monitoring both the accuracy and precision of measurements depending on the sample and task (thermal or light-absorption properties). Apart from the multi-signal feature, these joint measurements enable the separation of thermal and concentration parameters and should be accompanied, if possible, by the assessment of either light-absorption or thermophysical properties using other methods.
The obvious disadvantage of conventional thermal lensing is that the thermal-lens effect is ensemble-based, not molecule/structure-based, and non-microscopic measurements always provide the bulk properties of the sample. Still, the very large volume of information from transient thermal lensing simultaneously with the information channel from steady-state measurements provides a way to obtain rather valuable information. For instance, in the case of albumin interaction, thermal lensing can be used for screening various conditions or classifying the interactions by signature effects before more detailed experiments are performed. This feature of thermal lensing can be implemented even in a handheld device for rapid tests of complex samples alongside or before photothermalmicroscopy, light-scattering, or fluorescence measurements. Such an experimental design can be extended to pulsed thermal-lens excitation conditions [52,[175][176][177][178][179][180] not considered here. In that case, short excitation times close to or shorter than the thermal-lens characteristic times may increase the sensitivity due to the more substantial overheating of such entities as Hb globules, mitochondria, etc., and reveal features that are hidden or averaged at the early stages of signal development in the cw excitation conditions used in this study. It may provide, e.g., a means of estimation of the protein globule size in solution, which is beyond the abilities of the setups built in this study due to an insignificant contribution of individual nanolens entities (light-absorbing point sources) in bulk cw thermal lensing.
Certainly, these experiments may (and should) be combined with photothermal microscopic techniques, especially for transient measurements. Transient thermal lensing can overcome the diffraction limit of optical photothermal spectroscopy [41,181,182] to attain the super-resolution of photothermal imaging, which is highly important for various applications, such as studies of mitochondria, erythrocyte pathologies, and similar problems [146,181]. The following studies should integrate all photothermal measurements into a single compact device for micro-measurements [43].
Finally, despite the relatively reliable estimations attained in this study, a theory that combines the size and the nature of dispersed particles and the whole quantitative picture of the thermal-lens effect is not fully developed, and only the estimation of the range of the particle/aggregate/cluster size can be performed; addressing this could be the aim of advanced research in this direction.

Conflicts of Interest:
The authors declare no conflict of interest.

Appendix A. Photometry (Transmission) Modality of the Thermal-Lens Setup
For dye solutions (prepared according to Procedures S2-S4, Supplementary Materials) and nanodiamonds (Procedure S5, Supplementary Materials), the absorbance of solutions in the range of 450-540 nm (scanning step, 0.1 nm; scanning speed, 100 nm/min) was recorded in cells with a 10 mm optical path length using a Shimadzu UVmini 1240 spectrophotometer (Shimadzu Europa GmbH., Duisburg, Germany), and the baseline sample was distilled water. A thermal-lens spectrometer using an Ophir Optronics Nova II power meter (Ophir Optronics Solutions, Jerusalem, Israel) with a highly sensitive thermoelectric sensor recorded the radiation power without a cell (I e0 ), the reference channel, and the radiation power that passed through the cell with the solution under study (I e ). We registered the radiation intensity before and after the cell with gentian violet and cobalt nitrate solutions at all wavelengths of argon-ion and DPSS lasers (Table A1) and calculated absorbance values from transmittance (Equation (A1)). Measurements were carried out at the wavelengths of the coherent Innova 90 argon laser (514.5, 501.7, 496.5, 488.0, 476.5, 472.7, 465.8, 457.9, and 454.5 nm) and the wavelength of the second harmonic of the DPSS laser (532.0 nm). The incident radiation power was maintained as high as possible, determined by the characteristics and the working medium of the laser, and varied depending on the wavelength from 15 mW (454.5 nm) to 1 W (488.0 nm). Absorbance does not play any role, and the lasers worked with a fully open output aperture (>12 mm). The transmission T of the sample was assessed as T = I e /I e0 .
The distance from the cell to the power-meter sensor was about 2 cm. Such a small distance ensured that, with the proper control, all radiation passed onto the sensor crystal. Special attention was paid to concentrated solutions, in which a strong thermal lens is excited, causing a strong divergence of the excitation beam along with the probe beam.
In addition, for a correct comparison with absorbance measured by a spectrophotometer, power measurements were also carried out for a cell with distilled water.
Ferroin solutions were used to compare the linearity ranges of the laser photometer (primary wavelength of 488.0 nm and maximum power of ca. 1 W) and a spectrophotometer. The results are shown in Figure A1, which proves that under the selected conditions, the linear range of the photometer modality of the thermal-lens spectrometer is higher than a commercially manufactured spectrophotometer for routine measurements. Figure A1. Linearity ranges of calibration functions of a spectrophotometer (orange rhombs) and a laser photometer mode of a thermal-lens spectrometer (green squares) for ferroin solutions in water. Setup, TLS-60; wavelength, 488.0 nm; excitation power, 1 W.

PEER REVIEW 33
The correlation coefficients between the values measured on the spectrophotometer and laser spectrometer are above 0.9995. The correlation is high, indicating that the readings of both instruments are not significantly different. A similar procedure was performed for solutions of RUDDM nanodiamonds; the results show coefficients of correlation no less than 0.99; thus, we also observed high correlation coefficients between the spectrophotometer data and the laser photometer and the linearity of the data.

Appendix B. Modeling of Multipoint Secondary Heating
To estimate the thermal field in dispersed media, we used a previously developed numerical description for signal generation in photothermal lensing for multipointabsorbing (dispersed) solutions [96,108]. This approach is based on the calculation of the signal from the probe-laser phase shift. A single measurement of the transient photothermal signal ϑ(t) for an excitation on/off cycle of a cw excitation beam and a pulsed beam was calculated as the phase shift in the probe-beam wavefront Φ at a distance from the laser source z and a distance from the beam center r at time t ( is the probe-laser wavelength, l is the sample path length, dn/dT is the temperature coefficient of the refractive index, and ΔT is a photothermal temperature change) as a relative change in the probe-beam intensity [104,164]. The temperature change ( , , ) is found by summing individual thermal waves from multiple heat sources distributed in N thin layers plus the conventional contribution from the light-absorbing background with due regard to the random movement (shift) of heat sources in the signal generation The correlation coefficients between the values measured on the spectrophotometer and laser spectrometer are above 0.9995. The correlation is high, indicating that the readings of both instruments are not significantly different. A similar procedure was performed for solutions of RUDDM nanodiamonds; the results show coefficients of correlation no less than 0.99; thus, we also observed high correlation coefficients between the spectrophotometer data and the laser photometer and the linearity of the data. thin layers plus the conventional contribution from the light-absorbing background with due regard to the random movement (shift) of heat sources in the signal generation process.
∆T(z, r, t) = ∑ M n=1 ∑ N i=1 ∆T n,i (z, r, t) + ∑ M n=1 ∑ N i=1 ∆T n,i (z + Z, r + R, t) + ..., Here, M is the number of excitation pulses in a pulse train (cw excitation is simulated with a large train of ultrashort pulses). The second double-sum term represents the shifting of the heat element during the heating period by substituting the initial variables (r, z) into (r + R, z + Z), where (R, Z) is the shift vector. Thus, summing the shifted heat functions (Equation (A2)) gives the whole matrix of heat generation without binding it to the radial symmetry of the excitation beam. The model represents heat flow from every individual heat source, thus resulting in the exact spatially resolved thermal picture, which is much closer to the real situation in the description of dilute solutions [96,108]. The thermal parameters for heterogeneity are not considered in the calculation of the heat diffusion, as the contribution of the total heat-source volume to the total sample volume is negligible, and thermal heterogeneity is considered instantaneous; thus, only the photothermal parameters of the medium are taken into account. The application of a quasi-cw function for the description of the solution excitation does not prevent the calculation of transient functions of the temperature profile and signal.
The temperature distribution in layer i, ∆T i , is based on the separate consideration of the axial (below, in Equation (A4), the term in curly brackets) and radial (the term in square brackets) components of the temperature profile: here, I i−1 is the intensity of the excitation radiation incident to layer i. The parameter t 0 is the quasi-pulse duration (expresses the energy transmitted to the absorbing medium only), while the transfer of the total pulse energy is considered instantaneous [95]. The temperature response to cw irradiation is achieved by summing the temperature responses to a series of laser quasi-pulses, with each quasi-pulse starting at t n = (n -1)t 0 , where n is the quasi-pulse ordinal. Heat dissipation functions from different heat sources create the final temperature profile. The heat-wave profiles generated by different heat sources are described with the same function and photothermal response amplitudes that differ according to the particle position within the excitation beam.

Appendix C. Preliminary Optical Measurements of True Solutions
For a solution of MetHbCN in distilled water, the photometric calibration function (concentration in mol/L) is described by the equation: The detection limit of 200 nmol/L is not significantly different from that obtained in the absence of salts.
The corresponding calibration functions (concentration in mol/L) for the photometric determination of ferroin in water and salt (532.0 nm, ε 532 = 7500) are: The ratio of the sensitivity coefficients of Equations (A5)-(A10) with due regard to the measurement error agrees well with the ratio of the molar absorptivities of MetHbCN and ferroin and myoglobin and ferroin at 532 nm, i.e., 5.65 and 0.65, respectively.
Spectrophotometric measurements showed that in the wavelength range of thermallens spectrometers, albumin does not absorb. The absorption of selected silica particles can also be considered low (Supplementary Materials, Figure S2) A comparison of the absorption spectra of MetHbCN, myoglobin, albumin, and mixtures of albumin with dyes and MetHbCN before and after irradiation with laser radiation under the same conditions (excitation and exposure time) as in thermal-lens experiments showed the absence of changes, which verifies the absence of photochemical processes in these samples, as well as the lack of significant light scattering and fluorescence, which fully confirms the available data [183].