2.1. Reagents and chemicals
The study used dispersions of SiO
2 nanoparticles of three different series of Ludox (SM-30, HS-40, TM-50) purchased from Merck (certificate descriptions are presented in
Table 1). Aqueous dispersions were prepared using deionized water of the Milli-Q Reference system (Millipore SAS, Molsheim, France). We used the following solvents: Ethanol (C
2H
5OH, CAS: 64-17-5, chemically pure, Khimreaktiv, Klin, Russia), Chloroform (CHCl
3, CAS: 67-66-3, analytical grade, Khimreaktiv, Klin, Russia), Toluene (C
6H
5CH
3, CAS: 108-88-3, analytical grade, Lenreaktiv, St. Petersburg, Russia). Photostable indicator Ferroin (1,10-phenanthroline complex of iron(II), Fe(C
12H
8N
2)
3SO
4,
Mr = 596.27 g/mol, CAS: 14634-91-4, analytical grade, Lenreaktiv, St. Petersburg, Russia) were used as chromophores) and dye Sudan 1 (C
16H
12N
2O,
Mr = 248.28 g/mol, CAS: 842-07-9, chemically pure, Reakhim, Moscow, Russia). Mixing of solutions was carried out using a laboratory shaker model PE-6410 (EKROSKHIM, St. Petersburg, Russia).
For TLS, several samples of different brands were prepared with the same mass concentration of the solid phase of 1.60, 3.97, 8.79, 14.39, and 22.4 mg/mL (for TM series, 0.78, 2.80, 5.60, 11.2, and 16.8 mg/mL were additionally prepared). Ferroin was added to each sample. To reduce the systematic error caused by high absorbance, as described [
23], the concentration of ferroin in the dispersions did not exceed 1 µmol/L. After that, the dispersions were mixed for 2 h on a laboratory shaker and used in the analysis.
Based on the data in [
23], we used highly diluted solutions as reference samples, which were a solution of ferroin in water with a concentration of 1 µmol/L, and solutions of the dye Sudan 1 in chloroform, toluene, and ethanol with a concentration of 50 nmol/L.
2.2. Thermal-lens measurements
We apply the Shen–Snook model for a dual-beam TLS with mode-mismatch in the stationary state [
22]. We do not consider here the main conditions for the applicability of the model and briefly consider only the equations necessary to the work. The main equation for the intensity of the probe beam on the detector (transient curve equation) is:
in which the intensity of the probe beam
at each moment of the development of the thermal field (time
t) depends on the characteristic time
, the geometric parameters
m and
V, and the thermooptical signal
θ (
I(
0) is also present in the equation, the intensity of the probe beam at
t = 0). The characteristic time
is related to thermal diffusivity
D as:
where
is the radius of the excitation beam in the waist. Thus, knowing the characteristic time and the excitation beam waist, we can find the thermal diffusivity of the analyte. However, to apply Eq. (1), one must know three parameters:
m is the mode mismatch factor, which is described by the equation:
is the parameter representing the ratio of the radii of the probe (
) and excitation beams at the center of the sample;
V is the dimensional parameter of the spectrometer, which is described by the equation:
, which takes into account the distance of the probe beam waist to the sample (
), from the sample to the detector (
and confocal distance (Rayleigh distance) for the probe beam
); and the last parameter,
θ, which includes the excitation beam power (
P), linear absorption coefficient (
α), optical path length (
l), thermal conductivity (
k), excitation laser wavelength (
λe), refractive index temperature coefficient (
dn/
dT):
In Eq. (1), most parameters are constants and for simplicity, it can be represented as:
where
a,
b and
c are constants:
а = 2
mV,
b = (1 + 2
m)
2 +
V2,
c = 1 + 2
m +
V2.
The characteristic time from Eq. (6) can be represented as a function of time:
where
(
t) is the effective characteristic time at each point on the transient curve,
= 2[1–
I(∞)/
I(0)]/tan
–1(
a/
c), where
is the probe laser intensity in the stationary state. In the case of a homogeneous solution, the steady state means the average value of the last points of the transient curve (in our case, the average of the last 300 ms). Thermal diffusivity
, so-called effective thermal diffusivity, is calculated from Eq. (7) for each value
at time
t. The transition from effective
and
D to true values occurs by averaging the values of the first 100 ms of the functions
and
using Eq. (2).
Before each measurement, we measured the characteristic time and thermal diffusivity of reference samples, which are pure solvents (with a small chromophore content) with precisely known thermophysical and optical parameters. All necessary constants in the thermooptical parameter θ are taken from reference data and absorbance measurements.
We use the full development of the thermal field (a steady state). In this case, the thermal-lens signal is:
Here, in the case of the homogeneous reference sample, we also use the average value of the last 300 ms of the transient curve as the steady state, but in the case of a heterogeneous system, we use the minimum value of I(t), which is the inflection point. For comparisons of transient curves, we normalize them to the range 0÷1.
The development of the thermal-lens signal
was carried out using the calculation of the thermal-lens signal at each point. Eq. (8) was used for this, but
was used instead of
, where
is the probe beam intensity at each time point. Also, the development of the signal in time was normalized to the range of 0 ÷ 1:
where
is the signal in the equilibrium state, where is the inflection point of a transient curve,
is the thermal-lens signal at the initial moment of time.
To find the characteristic time and thermal diffusivity of heterogeneous systems, we used to fit the stationary state of the transient curve. This means that instead of
, the corrected intensity of the probe beam
is used, which is obtained by fitting the first 100 ms of the transient curve so that the last, 100
th, point of the experiment is fully consistent with the theoretical one. In this case, the characteristic time can be found according to Eq. (7) using the corrected intensity
:
The transition from effective characteristic time and thermal diffusivity (’ and D’, respectively) to true ones was carried out in the same way as in the case of homogeneous solutions described above, by averaging the data for first 100 ms.
Transient curves are presented in three different normalized forms. In the first form, normalization is carried out on the largest value,
. In the second and third forms, they are normalized according to the following equation:
where
in the second form is the intensity of the probe beam in the stationary state, which is found by averaging the last 300 ms of the transient curve. In the third form,
is the intensity of the probe beam in the steady state (intensity value of the probe beam at ca. 100 ms). In this case, for the second form, we denote the normalized intensity of the probe beam by
as in Eq. (11), without changes, and for the third form
(with an apostrophe).
The transient curve for thermal lens dissipation
was also normalized to two forms using the following equation:
where
is the intensity of the probe beam at time
t,
is the intensity of the probe beam at
t = 0, at the moment the shutter closes, is the smallest value;
in the first form is the intensity of the probe beam at stationary state, which is obtained by averaging the last 300 ms of the transient curve in the state of complete dissipation of the thermal lens and return of the beam to the level of pure water. In the second form, the intensity of the probe beam in the equilibrium state (intensity value at ca. 80–120 ms). In this case, for the first form, we denote the normalized intensity of the probe beam by
as in Eq. (12), without changes, and for the second form
(with an apostrophe).
The error is calculated according to the following equation:
where
is the measured value,
is the true value calculated theoretically or from a reference.
2.3. Thermal-lens spectrometer
The thermal lens spectrometer is described in detail in [
23].
Figure 1 shows the scheme of the dual-beam thermal-lens spectrometer. Radiation from the MGL-FN-532 solid-state laser (wavelength of 532 nm, TEM
00; Changchun New Industries Optoelectronics Tech. Co., Ltd, Changchun, PRC) passes through a shutter (model SH05, ThorLabs, Newton, NJ, USA), which is controlled by an analog-to-digital and digital-to-analog converter (ADC–DAC) model c8051Fx-DK (Silicon Labs, Boston, MA, USA) connected to a personal computer (PC), and enters the sample in a quartz cell (
l = 10.00 mm), in which a thermal lens is generated. A helium–neon laser HNL050L (wavelength of 632.8 nm, TEM
00; ThorLabs, Newton, NJ, USA,) was used as the probe laser. A photodiode was used as a detector.
The signal from the detector is recorded every cycle. One cycle starts when the shutter opens, continues when the shutter closes, and ends when the shutter is reopened, after which the cycle repeats. The PC receives data from the ADC and the detector, which are processed in the original program (C++ programming, Borland Corp., Austin, TX, USA), where the measurement cycles are formed, displayed, and stored in the form of transient curves (signal intensity vs. time). The laser power was measured using an Optronics Nova II power meter (Ophir Optronics Solutions, Jerusalem, Israel). The operating parameters of the measurements are summed up in
Table 2.
The measurement parameters for the thermal-lens spectrometer are optimized. The shutter frequency corresponds to the time of full development of the thermal lens and the transient curve reach the steady state. The selection of the number of transient curves for averaging is based on the recommendations made in [
23]. The mode-mismatch factor was selected based on the condition of the minimum systematic error identified in [
23] and corresponding to the conditions of the model [
22]. The sample-to-detector distance also corresponds to the conditions of the model [
22]. All systematic errors caused by the periodic divergence of the excitation beam, the relative displacement of two beams in the cell, the bias of the sample cell along the beam propagation, the position of the maximum intensity of the probe beam in the center of the detector are taken into account and reduced to zero, according to the recommendations made in [
23].