Concentration Quantification of TiO2 Nanoparticles Synthesized by Laser Ablation of a Ti Target in Water

In this work, we present a quantitative method for determining the concentration of metal oxide nanoparticles (NP) synthesized by laser ablation in liquid. The case study was performed with titanium dioxide nanoparticles (TiO2 NP), which were synthesized by laser ablation of a Ti target in water. After synthesis, a colloidal solution was analyzed with UV-Vis spectroscopy. At the same time, the craters that remained on the Ti target after ablation were evaluated with an optical microscope to determine the volume of the ablated material. SEM microscopy was used to determine the TiO2 NP size distribution. It was found that synthesized TiO2 NP followed a Log-Normal diameter distribution with a maximum at about 64 nm. From the volume of ablated material and NP size distribution, under the assumption that most of the ablated material is consumed to form nanoparticles, a concentration of nanoparticles can be determined. The proposed method is verified by comparing the calculated concentrations to the values obtained from the Beer–Lambert law using the Mie scattering theory for the NP cross-section calculation.


Introduction
Pulsed laser ablation in liquid (PLAL) is a versatile method for the synthesis of colloidal solutions of nanoparticles (NP) [1]. In this method, the laser pulses hit the target immersed in liquid and cause a laser-induced breakdown, followed by the shock wave and plasma plume formation, its expansion and cooling, cavitation bubble formation, and its expansion and collapse [2]. Nanoparticles are produced during the plasma cooling phase in the processes of nucleation and condensation. After the collapse of the cavitation bubble, nanoparticles are released and diffused into the liquid, leading to the formation of a colloidal solution of NP [2]. Because of the liquid environment, the plasma plume is very strongly confined and therefore has a much larger temperature, pressure, and density compared to the ablation in gas [3]. Unlike ablation in gas, there is no need for a closed chamber for NP collection in PLAL. The liquid layer around the plasma plume is also transformed into the plasma phase, leading to a mixture of two plasmas so chemical reactions can take place between them, thus affecting the NP structure [3]. For instance, the ablation of a Ti target in deionized water with a NS laser at a wavelength of 1064 nm results in the synthesis of TiO 2 NP, as in the present work [4].
A PLAL-synthesized colloidal solution of NP usually has large stability due to the presence of a negative charge on the NP surface [5]. The main advantages of PLAL over chemical methods of NP fabrication are the high purity of synthesized nanoparticles, the absence of unwanted residual byproducts in NP colloidal solution, and the possibility of NP synthesis from a large variety of materials [6]. The properties of NP properties, such as morphology, size distribution, and crystal structure, are affected by PLAL parameters such as laser fluence, wavelength, ablation spot size, pulse duration, type of liquid, and presence of surfactants in it. Therefore, there is a large space for modification and engineering of NP properties [7].
Most of these methods have serious limitations in achieving accurate determinations of NP concentration. For example, the size limitation of NP (in turbidimetry, resistive pulse-sensing, NP tracking analysis), overshadowing the signal of the small nanoparticles (in DLS), or the need for having calibration solution (in DLS and spICMS) [15]. An additional consideration for these methods is the difficulty in sampling, whereby diluting the as-prepared sample over several orders of magnitude is often needed. The main disadvantage of determining concentration from UV-Vis is the need for knowledge of the NP size and refractive indices.
In our previous work [27], a simple model for calculating PLAL-synthesized pure metallic NP concentrations is introduced that considers the volume of an ablated crater and NP size distribution. All disadvantages existing in the previously mentioned methods for concentration calculation do not exist here. The only limitations are those related to the accuracy of crater volume and size distribution determination. This model is successfully tested for Ag NP synthesized by laser ablation of Ag target using Beer-Lambert law.
In the present work, we show that this model can be adapted for the calculation of the concentration of metal oxide NP synthesized by PLAL of the metallic targets that are very prone to oxidation. The proposed method is tested for TiO 2 NP synthesized by the laser ablation of a Ti target immersed in water. TiO 2 NP have stability, corrosion resistance, low reactivity, long agglomeration, and sedimentation time. They have a broad range of industrial and scientific applications, especially photocatalysis [28]. The well-known phases of TiO 2 are anatase, rutile, and brookite, while amorphous TiO 2 , as one synthesized in this work, is much less investigated. Although amorphous TiO 2 is characterized by lower photocatalytic performance than one of the crystal phases, it can act as an active component for visible and near-IR light-harvesting, leading to improved photocatalytic activity in this part of the light energy spectrum, so it can also be applied in photodegradation of organic dyes, hydrogen production, and CO 2 photoreduction [29].
Furthermore, amorphous TiO 2 has lower toxicity and better corrosion resistance when compared to TiO 2 in the crystal phase, which can be an advantage in photoprotective and cosmetic applications [30]. It also exhibits excellent antibacterial properties [31], while a special energy band structure and improved charge separation may be advantageous in SERS [29]. Amorphous TiO 2 NP can be easily transformed to a crystalline phase by heat treatment [29], thus retaining the NP concentration value.
To verify this method, concentrations of TiO 2 NP calculated from crater volume (C V ) were compared to the concentrations obtained from the UV-Vis photoabsorbance data using Beer-Lambert law and Mie scattering theory (C A ). Beer-Lambert law is usually used for the calculation of the concentration of metal nanoparticles with a well-defined resonant frequency, for example, gold [24,32,33] or silver as in [27]. In the case of TiO 2, such calculation is more complicated because it is not plasmonic material. Still, it is a semiconductor that absorbs mainly in the UV range, so it does not possess resonant frequency. Furthermore, TiO 2 refractive index at some wavelength depends on many parameters such as TiO 2 bandgap, crystal structure, defects in it, and NP size [34,35], thus it is not easy to ensure that correct values of refraction index from literature are taken for calculation of C A . Moreover, the dipole approximation that was satisfied in [27] due to the sufficiently small size of synthesized Ag nanoparticles (<λ SPR /10, where λ SPR ≈ 400 nm is silver absorption peak) is not applicable in the present work dealing with NP, which are much larger and do not possess resonant frequency. Thus, in the present paper, model verification by the Beer-Lambert law is adapted when compared to the same method used in [27] to be applicable for the determination of semiconductor NP concentration, and the Mie theory for the calculation of TiO 2 NP extinction cross-section is exactly applied instead of using the dipole approximation.

Materials and Methods
A colloidal solution of TiO 2 NP was synthesized by the pulsed-laser ablation of the Ti target (purity 99.9% and thickness 3 mm, Kurt J. Lesker) immersed in a beaker containing 40 mL of deionized water using Nd:YAG laser (Quantell, Brilliant). Laser specifications are a pulse duration of 4 ns, wavelength of 1064 nm, output energy of 290 mJ and repetition rate of 5 Hz. The laser beam was directed by a system of prisms and focused by a lens (the focal length of 10 cm) onto the target surface. The laser pulse energy in front of the target was 210 mJ while a diameter of a focused pulse on the target surface was 1 mm, which yielded a laser fluence of 27 J/cm 2 . The thickness of a water layer above the target was kept constant at 2 cm during the experiment to keep the ablation efficiency and thus the NP properties, constant [36]. The target was fixed during the experiments to allow the drilling of a crater and thus determine ablated material mass. A scheme of the experimental setup for PLAL is depicted in detail in [37].
After the experiments, the craters created on the targets were studied with an optical microscope (Leica DM2700M, Leica Microsystems, Wetzlar, Germany). Microscopical images were taken at different focal positions with respect to the target surface in order to obtain crater radius dependence on crater depth (i.e., crater profile). Using the obtained crater profile, crater volume is calculated according to the procedure described in [38].
The optical absorbance of laser synthesized colloidal solutions was assessed using a UV-Vis spectrophotometer (Lambda 25, Perkin Elmer, Waltham, MA, USA) in the wavelength range of 200-600 nm.
The crystallinity and crystalline phases were studied by grazing incidence X-ray diffraction (GIXRD). To perform a structural characterization of TiO 2 NP, the produced colloid was dropped onto a silicon substrate and left to air-dry to obtain TiO 2 NP film. The crystalline structure of TiO 2 NP was investigated using a D5000 diffractometer (Siemens, Munich, Germany) in a parallel beam geometry with Cu Kα radiation, a point detector, and a collimator in front of the detector. Grazing incidence X-ray diffraction (GIXRD) scans were acquired with the constant incidence angle α i of 1 • , ensuring that the information contained in the collected signal covers the entire film thickness.
The morphology and size distribution of the synthesized nanoparticles was studied by a field emission scanning electron microscope (FE-SEM, JSM-7600F, Jeol Ltd., Tokyo, Japan). Samples for SEM imaging were prepared by depositing one drop of a suspension on a polished Al sample holder. The holder with the specimens was coated with a 3 nm-thick amorphous carbon layer (PECS 682) prior to the SEM imaging. SEM images were analyzed with ImageJ to determine the NP size distribution.
A Zetasizer Ultra (Malvern Panalytical, Malvern, UK) was used to determine the zetapotential (ζ) by electrophoretic light scattering (ELS). It was calculated from the measured electrophoretic mobility by means of the Henry equation using the Smoluchowski approximation (f(k a ) = 1.5). Results are reported as an average value of three measurements. The data processing was done by the ZS Xplorer 1.20 (Malvern Panalytical).
Numerical calculations of extinction, scattering, and absorption cross-sections of TiO 2 NP using equations from Mie scattering theory are performed by Mätzler's MATLAB code [39].

TiO 2 NP Characterization
The laser ablation of metallic Ti target in water resulted in the synthesis of a colloidal solution of semiconductor TiO 2 NP. The GIXRD measurements shown in Figure 1 revealed that TiO 2 NP are amorphous due to the lack of appearance of any of the major Bragg The ζ-potential of the produced TiO 2 NP colloidal solution was measured to be 30 ± 1 mV, making the solution of moderate stability (aggregation and precipitation of nanoparticles appeared after two days). Therefore, TiO 2 NP concentration can also be considered a homogenous and well-defined value in the as-prepared solution.

TiO2 NP Characterization
The laser ablation of metallic Ti target in water resulted in the synthesis of a colloidal solution of semiconductor TiO2 NP. The GIXRD measurements shown in Figure 1 revealed that TiO2 NP are amorphous due to the lack of appearance of any of the major Bragg peaks of TiO2 phases (theoretical peaks for rutile and anatase are shown for comparison). The ζ-potential of the produced TiO2 NP colloidal solution was measured to be 30 ± 1 mV, making the solution of moderate stability (aggregation and precipitation of nanoparticles appeared after two days). Therefore, TiO2 NP concentration can also be considered a homogenous and well-defined value in the as-prepared solution. The morphology and size distribution of TiO2 NP were obtained using SEM imaging. A typical SEM image of TiO2 NP is shown in Figure 2a, where 5000 pulses were applied in ablation. It can be seen that TiO2 NP are spherical and with a broad range of sizes. For each number of laser pulses applied in ablation, the size distribution was obtained by measuring the diameters of 200 single nanoparticles from SEM images. It was found that the size distribution was similar for all tested numbers of pulses. The size distribution obtained from SEM images is shown in Figure 2b. The morphology and size distribution of TiO 2 NP were obtained using SEM imaging. A typical SEM image of TiO 2 NP is shown in Figure 2a, where 5000 pulses were applied in ablation. It can be seen that TiO 2 NP are spherical and with a broad range of sizes. For each number of laser pulses applied in ablation, the size distribution was obtained by measuring the diameters of 200 single nanoparticles from SEM images. It was found that the size distribution was similar for all tested numbers of pulses. The size distribution obtained from SEM images is shown in Figure 2b. It can be seen that the size distribution is relatively broad, and it was fitted with the Log-Normal function with the maximum at 64 nm. Log-Normal is a very common func tion for fitting the size distribution of nanoparticles synthesized by PLAL [40][41][42][43]. Alt hough non-NP species are, to some degree, also present in synthesized colloidal solutions such as large debris structures, from the obtained size distribution in Figure 2b, it can be assumed that their presence is negligible, so the term "TiO2 NP" is justified. Since the na noparticles are spherical, the average NP volume is calculated as = , where d is nanoparticle diameter and can be obtained from the size distribution. From the given size distribution in Figure 2b, the average volume of nanoparticles is calcu lated as: It can be seen that the size distribution is relatively broad, and it was fitted with the Log-Normal function with the maximum at 64 nm. Log-Normal is a very common function for fitting the size distribution of nanoparticles synthesized by PLAL [40][41][42][43]. Although non-NP species are, to some degree, also present in synthesized colloidal solutions, such as large debris structures, from the obtained size distribution in Figure 2b, it can be assumed that their presence is negligible, so the term "TiO 2 NP" is justified. Since the nanoparticles are spherical, the average NP volume V NP is calculated as V NP = 1 6 d 3 π, where d is nanoparticle diameter and d 3 can be obtained from the size distribution. From the given size distribution in Figure 2b, the average volume V NP of nanoparticles is calculated as: where M is the total number of size ranges in the size distribution, d i is the average diameter that corresponds to size range i (single column bar in size distribution), and n i = N i /N is the ratio between the number N i of nanoparticles corresponding to size range i (relative abundance of NP within corresponding column bar) and the total number of nanoparticles N.

Calculation of TiO 2 NP Concentration from Crater Volume (C V )
The mass of synthesized TiO 2 NP depends on the amount of ablated material by PLAL, represented as a volume of a crater left on the target after ablation. Due to the fact that the target was made of almost pure titanium (thickness of surface oxidation layer can be neglected when compared to crater volume), the amount of titanium in the formed TiO 2 NP is coming solely from the target while the oxygen is coming from the water. The volume of the ablated crater V crat is turned into a fraction of titanium in the synthesized TiO 2 NP. The total number N of synthesized nanoparticles can be determined from the known volume of ablated material V crat by performing the following calculation: where N is the total number of nanoparticles in colloidal solution, V NP is the average volume of nanoparticles, m(TiO 2 ) is the total mass of TiO 2 NP in colloidal solution, ρ(Ti) is Ti target density (4.506 g/cm 3 near R.T.), ρ(TiO 2 ) is the TiO 2 NP density (3.8 ± 0.1 g/cm 3 ) for amorphous TiO 2 [44], while m(Ti) and m(O 2 ) are the total masses of titanium and oxygen atoms in TiO 2 NP in colloidal solution, respectively. A(O) and A(Ti) are atomic mass numbers for oxygen atom and titanium atom, which are 16 and 48, respectively. As shown in [45], TiO 2 NP of 21 nm in size have only 4% larger density than TiO 2 particles 200 nm in size, so it can be concluded that, due to the prevalence of large nanoparticles in the TiO 2 NP size distribution (Figure 2b), it is justified to neglect the NP size dependence of ρ(TiO 2 ) in Equation (2). Number concentration C V [mL −1 ] of TiO 2 nanoparticles is then Equation (3): where total number of synthesized nanoparticles is divided by an amount of liquid V liquid where synthesis was performed (in our case V liquid = 40 mL). In Figure 3a, the profiles of craters left after ablation with 1000, 3000, and 5000 applied laser pulses are shown (note that the x-y axis is not in scale). The dashed line denotes the target surface. Craters are characterized by a relatively large aspect ratio (radius to depth), as the depth of a crater is about five times lower than the crater radius. In Figure 3b, the volumes of crater V crat dependent on the number of pulses are shown, where volumes are calculated from profiles in Figure 3a, as described in [38]. For each number of applied pulses, the total number of nanoparticles in water is calculated by inserting the measured V crat in Equation (2). It provides the following number of nanoparticles dependence on number of pulses: N (1000p) = 7.69 × 10 9 , N (3000p) = 15.94 × 10 9 , N (5000p) = 20.04 × 10 9 . When the number of nanoparticles N is divided by water volume V liquid = 40 mL, as in Equation (3), the concentration C V of TiO 2 NP is obtained for each number of applied pulses. Their values are C V (1000p) = (1.9 ± 0.2) × 10 8 mL −1 , C V (3000p) = (4.0 ± 0.4) × 10 8 mL −1 and C V (5000p) = (5.0 ± 0.5) × 10 8 mL −1 .

Calculation of TiO2 NP Concentration CA from Beer-Lambert Law (CA)
UV-Vis measurements shown in Figure 4 are performed to obtain the concen CA of TiO2 NP in colloidal solution using the Beer-Lambert law. Indirect bandgap en for each TiO2 colloidal solution are calculated from the photoabsorbance measure using the Tauc plot as shown in the inset of Figure 4, and their values are 2.87 eV, 2 and 3.08 eV for 1000, 3000, and 5000 laser pulses, respectively. The calculated in bandgap energies are close to the indirect bandgap energy 3.0 eV obtained in [ amorphous TiO2 thin films but are also close to the bandgap energies expected

Calculation of TiO 2 NP Concentration from Beer-Lambert Law (C A )
UV-Vis measurements shown in Figure 4 are performed to obtain the concentration C A of TiO 2 NP in colloidal solution using the Beer-Lambert law. Indirect bandgap energies for each TiO 2 colloidal solution are calculated from the photoabsorbance measurements using the Tauc plot as shown in the inset of Figure 4, and their values are 2.87 eV, 2.98 eV, and 3.08 eV for 1000, 3000, and 5000 laser pulses, respectively. The calculated indirect bandgap energies are close to the indirect bandgap energy 3.0 eV obtained in [46] for amorphous TiO 2 thin films but are also close to the bandgap energies expected for the most common TiO 2 phases-rutile (3.0 eV) and to a lesser extent, anatase (3.20 eV) [28].  For the calculation of the NP concentration CA from the Beer-Lambert law size-distribution data (Figure 2b) was used for CV calculation. If φin (λ) is the e φout (λ) is the output flux of light with wavelength λ through the TiO2 colloidal (λ) is the absorbance at wavelength λ, defined as in Equation (4): If it is assumed that the colloidal solution contains spherical TiO2 NP wi ferent size ranges, each with a concentration ci (i = 1 to M) and average optica cross-section σi (λ) within the size-range i, the absorbance A can be expressed b Lambert law as in Equation (5): where l = 1 cm is the absorption path length in the cuvette. In our experiment, ance measurements were done immediately after the production of a colloid as such, there was no precipitation and agglomeration of TiO2 NP. From the For the calculation of the NP concentration C A from the Beer-Lambert law, the same size-distribution data (Figure 2b) was used for C V calculation. If φ in (λ) is the entrant flux, φ out (λ) is the output flux of light with wavelength λ through the TiO 2 colloidal solution, A (λ) is the absorbance at wavelength λ, defined as in Equation (4): If it is assumed that the colloidal solution contains spherical TiO 2 NP within M different size ranges, each with a concentration c i (i = 1 to M) and average optical extinction cross-section σ i (λ) within the size-range i, the absorbance A can be expressed by the Beer-Lambert law as in Equation (5): where l = 1 cm is the absorption path length in the cuvette. In our experiment, the absorbance measurements were done immediately after the production of a colloidal solution; as such, there was no precipitation and agglomeration of TiO 2 NP. From the ζ-potential measurements, it was previously concluded that the TiO 2 NP colloidal solution was mod-erately stable. Therefore, the concentration was approximatively the same in the whole solution and can be expressed as: where n i is abundance (∑ M i=1 n i = 1) of TiO 2 NP in the corresponding size range. Abundance n i is calculated from the size distribution (Figure 2b) for each size range. If the expression in Equation (6) for c i is inserted in Equation (5), the following expression for concentration C A is obtained: where σ ext (λ) = ∑ i=M i=1 n i σ i (λ) is the average extinction cross-section of TiO 2 NP at wavelength λ. σ i (λ) is approximated with extinction cross-section σ ext (d i , λ) corresponding to NP diameter d i , which is the arithmetic mean of size-range i. The Mie theory is used for the calculation of σ ext (d i , λ). According to the Mie theory, the extinction cross-section σ ext (d) of a spherical particle with diameter d is the sum of the corresponding absorption and scattering cross-section, as shown in Equation (8): where Q ext is the extinction coefficient, which is the sum of the absorption coefficient Q abs and the scattering coefficient Q scatt . According to [47] Q ext , Q scatt , Q abs can be calculated from Equations (9)-(11): Q abs = Q ext − Q scatt (11) where x is the size parameter defined as: k is the wavenumber of light in a medium (water). n medium is a real number, due to assumed water transparency. Wavelength-dependent values for water refraction indices used in all calculations are comprehensively listed in [48,49] for water temperature 25 • C. The Mie coefficients a n and b n can be calculated from Equations (13) and (14) [47]: where j n and h n are spherical Bessel's functions of order n. m is defined in Equation (15): m = n n medium (15) where n is the TiO 2 refractive index, which is a complex number, defined as in Equation (16): where n and n" are real and imaginary parts of the refractive index, respectively. In order to calculate wavelength dependence of TiO 2 NP average extinction crosssection σ ext (λ) = ∑ i=M i=1 n i σ ext (d i , λ) needed for C A calculation by Equation (7), n i values were taken from NP size distribution (Figure 2b). Equations (8)- (16) were used for the calculation of σ ext (d i , λ) for each d i , with input values of wavelength and wavelengthdependent TiO 2 and H 2 O refractive indices. The refractive index of amorphous TiO 2 is similar to the one of anatase, and rutile to a lesser extent, as obtained in [50]. Therefore, using the values of refractive indices reported in the literature for crystallized TiO 2 is justified for the purpose of estimating the extinction cross-section σ ext (λ) of the TiO 2 NP synthesized in this work.
To make the best estimation of σ ext (λ), calculations were made for three different dependencies of TiO 2 refractive index on wavelength found in the each of the following works: Siefke et al. [51], Sarkar et al. [52], and Bodurov et al. [53,54], which are comprehensively listed on the website [55]. Siefke et al. [51] performed the WGP (wire grid polarizer) technique for the determination of refractive indices in ALD-prepared TiO 2 thin film with a thickness of 350 nm in the wavelength range of 120 nm-125 µm using TiO 2 material with an indirect bandgap at 3.2 eV. Sarkar et al. [52] have used opto-plasmonic sensors for the determination of the refractive indices in TiO 2 rutile thin film with a thickness of 200 nm in a wavelength range of 300 nm-1.69 µm. Bodurov et al. [53,54] used Bruggeman's effective medium approximation for the calculation of refractive indices in TiO 2 anatase nanoparticles (diameter smaller than 35 nm) dispersed in water in the wavelength range 405-635 nm; the measurements were done with a laser micro-refractometer. The calculations of Equations (8)- (16) for the determination of σ ext (d i , λ) were performed numerically, using Matzler's MATLAB code [39]. Figure 5 shows the wavelength-dependence of cross-sections σ ext (λ) calculated for TiO 2 refractive indices from all of the three mentioned works (Siefke et al. [51], Sarkar et al. [52], Bodurov et al. [53,54]) in the wavelength range 390-600 nm. From Figure 5, it can be concluded that extinction cross-sections for all three cases are very similar to each other at wavelengths 390-415 nm, corresponding to a common TiO 2 bandgap energy range 3.0-3.2 eV. For larger wavelengths, the difference between them is much greater. Therefore, Equation (7) will provide the best estimations of TiO 2 NP concentration C A while inserting A(λ) and σ ext (λ) at wavelength λ that corresponds to bandgap energy. Furthermore, both the experimental absorbance A(λ) and Mie extinction cross-section σ ext (λ) have larger values and therefore lower relative error at lower wavelengths, so this is another reason why C A calculation by Equation (7) has the highest accuracy at wavelengths corresponding to bandgap energy. The third advantage of such an approach is the fact that bandgap energy is easily calculated by the Tauc plot (as in Inset of Figure 4), so it is exactly known which experimental absorbances correspond to the bandgap and should be inserted in Equation (7) for the calculation of C A . These absorbances are 0.042 for 1000p, 0.072 for 3000p, and 0.090 for 5000p.  Figure 6 shows the dependence of absorbance at bandgap energy on the umes. It can be seen that this dependence is linear, so it means that the avera cross-section of TiO2 NP does not depend on the number of laser pulses. Th expected due to the similarity in the TiO2 NP size distribution at a different n pulses. Therefore, for each number of pulses, the same average optical cross-s be inserted in Equation (7) for CA calculation. The extinction cross-section i Equation (7) for the calculation of CA is taken from calculations made by usin fractive indices from Siefke et al. [51] and has a value of = 5.42 × 10 −10 cm length 390 nm ( Figure 5), which corresponds to the bandgap of TiO2 used in paper. The selection of the paper by Siefke et al. [51] is made due to the well-de bandgap, which is indirect-like the TiO2 NP in this paper. However, the differ that would occur in the case of using TiO2 refractive indices from the other tw (Sarkar et al. [52], Bodurov et al. [53,54]) is included as the contribution to the u Figure 5. Extinction cross-section dependence on wavelength calculated from the Mie theory for TiO 2 solution of spherical nanoparticles using SEM size-distribution and TiO 2 refractive indices from three different works: Siefke et al. [51], Sarkar et al. [52], and Bodurov et al. [53,54] (listed on the website [55]). Figure 6 shows the dependence of absorbance at bandgap energy on the crater volumes. It can be seen that this dependence is linear, so it means that the average optical cross-section of TiO 2 NP does not depend on the number of laser pulses. This result is expected due to the similarity in the TiO 2 NP size distribution at a different number of pulses. Therefore, for each number of pulses, the same average optical cross-section can be inserted in Equation (7) for C A calculation. The extinction cross-section inserted in Equation (7) for the calculation of C A is taken from calculations made by using TiO 2 refractive indices from Siefke et al. [51] and has a value of σ ext = 5.42 × 10 −10 cm 2 at wavelength 390 nm ( Figure 5), which corresponds to the bandgap of TiO 2 used in the same paper. The selection of the paper by Siefke et al. [51] is made due to the well-defined TiO 2 bandgap, which is indirect-like the TiO 2 NP in this paper. However, the difference in C A that would occur in the case of using TiO 2 refractive indices from the other two papers (Sarkar et al. [52], Bodurov et al. [53,54]) is included as the contribution to the uncertainty of C A . The calculated values of concentrations C A are: C A (1000p) = (1.8 ± 0.3) × 10 8 mL −1 , C A (3000p) = (3.1 ± 0.5) × 10 8 mL −1 and C A (5000p) = (3.8 ± 0.6) × 10 8 mL −1 .
x FOR PEER REVIEW Figure 6. The absorbance at the wavelength corresponding to the bandgap energy vs. th volume for TiO2 colloidal solution synthesized with PLAL. Linear fit.

Discussion
The concentrations of CV and CA as calculated here are both listed in Tabl shown in Figure 7 for each number of applied pulses. CA is close to CV, but slightly 5% for 1000p, 22% for 3000p, and 24% for 5000p. Therefore, our proposed method determination of laser-synthesized TiO2 NP concentration from the size distribut volume of the crater remained on an ablated target is verified, at least within the l obtained uncertainty, which are acceptable for many practical purposes.
There are a few possible explanations why CA is lower than CV. First, the actu of the optical cross-section may be smaller than the calculated value due to the scattering effects that occur on TiO2 NP in colloidal solution. Second, there is a po that the crater volume was slightly larger than the volume of Ti material include formation of TiO2 NP. This was due to the crater expansion that may have occurred the Coulomb explosion in the process of laser ablation or due to the synthesis of str that are not observable by SEM images, like Ti ions. Third, the TiO2 NP analyzed paper are amorphous, so they probably have, according to [50], a slightly smaller tive index than the one used for the calculation of σ ext from the paper by Siefke et The greatest contribution to the calculated concentration uncertainties is related size-distribution uncertainty, which affects the calculation accuracy of the average and average optical cross-section of colloidal TiO2 NP.

Discussion
The concentrations of C V and C A as calculated here are both listed in Table 1 and shown in Figure 7 for each number of applied pulses. C A is close to C V , but slightly lower: 5% for 1000p, 22% for 3000p, and 24% for 5000p. Therefore, our proposed method for the determination of laser-synthesized TiO 2 NP concentration from the size distribution and volume of the crater remained on an ablated target is verified, at least within the limits of obtained uncertainty, which are acceptable for many practical purposes. Table 1. TiO 2 concentrations calculated from the crater volume (C V ) or from the Beer-Lambert law (C A ) per number of pulses. The proposed method for the calculation of TiO2 NP concentration using cra umes has significant advantages over the method using Beer-Lambert law. In the lation of concentration using the Beer-Lambert law, a complex Mie theory is appli therefore, numerical calculations are needed, the NP refractive indices must be k there usually exists the upper limit of solution density for the correct absorbance d nation, the high concentration homogeneity in solution is required, and in order more accurate results, the angular scattering in Mie theory should be considered issues do not exist in the method proposed in this article for the calculation of T concentration from ablated crater volumes.

Conclusions
The method for the calculation of PLAL-synthesized NP concentration, origin troduced in our previous work [27] where it is applied for the calculation of pure m NP concentration (Ag NP), was herein adapted for the calculation of metal oxide N centration (TiO2 NP). It is based solely on the determination of the volume of an material crater and the determination of the size distribution of nanoparticles. In o verify this method, concentrations of TiO2 NP in PLAL-synthesized colloidal so obtained by the presented method are compared to the concentrations calculated the Beer-Lambert law. Concentrations obtained from both methods are similar: C 25% smaller than CA in the TiO2 NP concentration range (1.8-5.0) × 10 8 mL −1 . There can be concluded that the proposed method is verified. This method has many adva compared to traditional methods of concentration determination because it does quire a calibration solution, there is no uncertainty related to the nature of light inte with nanoparticles, and the only requirement for NP size is the possibility of thei There are a few possible explanations why C A is lower than C V . First, the actual value of the optical cross-section may be smaller than the calculated value due to the multiscattering effects that occur on TiO 2 NP in colloidal solution. Second, there is a possibility that the crater volume was slightly larger than the volume of Ti material included in the formation of TiO 2 NP. This was due to the crater expansion that may have occurred during the Coulomb explosion in the process of laser ablation or due to the synthesis of structures that are not observable by SEM images, like Ti ions. Third, the TiO 2 NP analyzed in this paper are amorphous, so they probably have, according to [50], a slightly smaller refractive index than the one used for the calculation of σ ext from the paper by Siefke et al. [51]. The greatest contribution to the calculated concentration uncertainties is related to NP size-distribution uncertainty, which affects the calculation accuracy of the average volume and average optical cross-section of colloidal TiO 2 NP.
The proposed method for the calculation of TiO 2 NP concentration using crater volumes has significant advantages over the method using Beer-Lambert law. In the calculation of concentration using the Beer-Lambert law, a complex Mie theory is applied and, therefore, numerical calculations are needed, the NP refractive indices must be known, there usually exists the upper limit of solution density for the correct absorbance determination, the high concentration homogeneity in solution is required, and in order to have more accurate results, the angular scattering in Mie theory should be considered. These issues do not exist in the method proposed in this article for the calculation of TiO 2 NP concentration from ablated crater volumes.

Conclusions
The method for the calculation of PLAL-synthesized NP concentration, originally introduced in our previous work [27] where it is applied for the calculation of pure metallic NP concentration (Ag NP), was herein adapted for the calculation of metal oxide NP concentration (TiO 2 NP). It is based solely on the determination of the volume of an ablated material crater and the determination of the size distribution of nanoparticles. In order to verify this method, concentrations of TiO 2 NP in PLAL-synthesized colloidal solutions obtained by the presented method are compared to the concentrations calculated using the Beer-Lambert law. Concentrations obtained from both methods are similar: C V is 5-25% smaller than C A in the TiO 2 NP concentration range (1.8-5.0) × 10 8 mL −1 . Therefore, it can be concluded that the proposed method is verified. This method has many advantages compared to traditional methods of concentration determination because it does not require a calibration solution, there is no uncertainty related to the nature of light interaction with nanoparticles, and the only requirement for NP size is the possibility of their detection by microscope. However, to apply this method correctly, the size distribution should be determined with high certainty because it is an important parameter that has a large impact on the calculated value of NP concentration, so this is the main limitation of this method. Furthermore, the method is not applicable when nanoparticles in colloidal solution are very inhomogeneous in terms of their stoichiometry or shape. It can be expected this method will also work for the calculation of the concentration of other metal oxide NP synthesized by PLAL and also while using other types of laser in PLAL, such as a pulsed femtosecond or microsecond laser, but further research may be carried out for confirmation.