3.1. Numerical Simulation in COMSOL
For a better explanation of the deposition phenomena, it is also necessary to understand the thermal phenomena developed on the target surface in the spot of incidence of the laser beam in order to evaluate the magnitude of ablation as well as the fluid phase of the resulting material. Since the temperature developed during the interaction of the laser beam with the material cannot be measured due to the very short action time, it was estimated through the results of numerical simulation in the COMSOL software. The theoretical and mathematical model is according to Cocean et al. 2017 [
18] and Cocean et al. 2024 [
17].
Figure 2 shows the results of the numerical simulation in COMSOL that reflect the thermal effects in space and time produced on the Cu target and on the Ag target in space (
Figure 2a) and in time (
Figure 2b). The numerical simulation takes into account the laser parameters according to those used for experimental deposition and the material constants and parameters of the two materials used as targets (Cu and Ag). The optical parameters of Cu and silver materials were taken from the database [
19,
20], namely the refractive index n(Cu) = 1.1159 and n(Ag) = 0.054007, while the extinction coefficients k(Cu) = 2.5956 and k(Ag) = 3.4290, values to which correspond the absorption coefficients α(Cu) = x 10
5 cm
-1 and α(Ag) = 8.0996 x 10
5 cm
-1. The absorption coefficients are important both for calculating the heat source and the thermal effect resulting from the interaction of the laser beam with the target material, as well as for being able to optimize the meshing for finite element calculation so that no values of interest are lost, but also so that the computer's working memory is not overloaded. The other material parameters were taken from the COMSOL software library as functions of temperature, extrapolating with the closest function for temperature values higher than those provided in the functions supplied by COMSOL.
As a reference for evaluating the thermal effect, the melting points of the two materials were taken into account (M.P. Cu = 1357.77 K; M.P. Ag = 1234.95 K). The melting point is considered in this numerical study as the temperature above which the fluid phase results in the laser ablation process and which is important for the deposition, also influencing the amount of material deposited. For comparing the results obtained for the two materials, the dimensions on which the conditions for the formation of the liquid phase are met on each target were also measured at two time values (t=10 ns and t=20 ns) and which are indicated on the graph by Δx. Also, the results of the thermal effect were compared through the time during which the conditions for the existence of the fluid phase were met (temperature greater than or equal to the melting point), this parameter being indicated on the graph with Δt.
For the Cu target, the simulation results show a maximum of 9394 K at time t=11 ns and values of 9230 K during the pulse (t=10 ns) and 1458 K at the end of the studied interval (t=20 ns). The time interval during which the liquid phase occurs was calculated to be 17.5 ns, and the spot size where the solid phase threshold is reached ranges between Dx=9.9 mm at time t=10 ns and Dx=1.2 mm at time t=20 ns. Regarding the silver target, a maximum temperature resulting from the interaction with the laser radiation of 1811 K is observed at time t=11 ns, and a pulse temperature (t=10 ns) of 1797 K, which ensures the production of the liquid phase. At the end of the studied interval (t=20 ns), the temperature at the center of the laser-irradiated spot drops well below the melting point of silver (T=515 K). From the studied time interval, the time during which the liquid phase occurs is Δt=7 ns. The area of the spot affected by the laser irradiation that reaches temperatures greater than or equal to the melting point is at time t=10 ns equal to Δx=3.2 mm.
These values show an estimated ablation yield relative to the fluid phase much higher for copper compared to silver. Much higher temperatures resulting on the surface of the Cu target during laser irradiation lead to the production of a fluid with lower viscosity, which should be reflected in a deposit with lower amplitude roughness, and the obtained thin film should be more consistent, with a greater thickness than in the case of silver.
3.2. Scanning Electron Microscopy and Energy Dispersive X-Ray (SEM-EDS)
In the SEM images from
Figure 3 a and b, the Cu thin film is smoother and more uniform in its morphology compared to the Ag/Cu double film. For the latter, the aspect of the surface is given by the top layer which is silver and it is similar to the thin film obtained in a previous study reported by Cocean et al., 2024 [
17]. The Cu thin film [
Figure 3a] presents round-shaped particles, similar to pearls of maximum 3 μm diameter which are in small number, and the large majority of the particles with their diameter less than 1 μm, meaning hundreds of nanometers as also Lorusso et al., 2023 reported for the thin films they obtained by PLD even though they used an UV laser beam and different parameters and conditions [
4]. The silver thin film on top of the copper thin film in the Ag/Cu sample (
Figure 3b), shows numerous structures in the range of 10-20 μm. These structures were explained as being the result of the fluid instabilities [
17] and it can be seen that they are ring-shaped (crown-splash), pearl-shaped and aggregated structures. The most of the structures are around 5 μm diameter and less. The difference between the surface morphology of the thin Cu film and that of Ag obtained by laser ablation and deposition lies in the efficiency of producing the fluid phase as well as in the quality of the fluid phase (lower viscosity) of copper compared to silver, as shown by the COMSOL simulation results presented in
Section 3.1. The Cu-RB21 sample presents crystalline structures on all the analyzed areas (
Figure 2c,d) reflecting a good dispersion of the RB21-NaHCO
3 aqueous solution that could be influenced by the copper when compared with the lower dispersion that was obtained on the silver thin film reported by Cocean et al., 2024 [
14]. The good spreading of RB21-NaHCO
3 aqueous solution can be explained by the copper dissolution due to NaHCO
3 resulting copper ions. The copper ions enter into reaction with the bicarbonate, replacing the sodium ions. However, the amount of copper that reacts is small due to the high reducing potential of copper compared to that of sodium [
13]. Also, the morphology of the thin Cu layer shows fewer droplet and crown splash type structures and they are also smaller in size compared to those in the morphology of the thin Ag layer and the double Ag/Cu layer. The RB21-NaHCO
3 on the double thin film Ag/Cu (
Figure 3e,f) appears to be less dispersed, similar to what was observed for the silver thin film [
17].
The EDS spectrum of Cu thin film
Figure 4a) and the elemental composition (
Figure 4c) show a strong adsorption of oxygen on the thin film, the oxygen being in proportion of 95.5% weight on the thin film surface. For the double film Ag/Cu where silver is the top thin film, based on its EDS spectrum (
Figure 4b) and its elemental composition (
Figure 4d), the adsorbed oxygen is reduced to only 52% weight. That makes copper a better oxygen sorbent, but also more sensitive to oxidation. This property may play an important role if the sorption can be reversed by controlled mechanism.
3.3. FT-IR Spectroscopy
In the spectra presented in
Figure 5 with the bands of vibration listed in
Table 1, the bands assigned to O-H free stretching are noticed at 3861 cm
-1, 3810 cm
-1, 3727 cm
-1 (RB21), 3857 cm
-1, 3734 cm
-1 (RB21-NaHCO
3), 3850 cm
-1, 3760 cm
-1, 3647 cm
-1 (R-Cu), 3803 cm
-1, 3740 cm
-1 (Q-Cu), 3857 cm
-1, 3734 cm
-1 (R-Ag/Cu) and 3857 cm
-1, 3743 cm
-1, 3647 cm
-1 (Q-Ag/Cu). The O-H and N-H free and bonded strong stretching vibrations are observed at 3449 cm
-1 (RB21) and 3439 cm
-1 (RB21-NaHCO
3) in the control samples, while bands reduced in intensity are noticed at 3504 cm
-1 -3426 cm
-1 and 3510 cm
-1 for the Q-Cu and Q-Ag/Cu, respectively. They are missing for the R-Cu and R-Ag/Cu samples. These changes in the spectra indicate severe modification in both the chromophore and the chromogen structures. The vibrations assigned to N-H bending and to the pyrrole heteroatomic rings skeleton at 1594 cm
-1 and 1508 cm
-1 in RB21 are missing in RB21-NaHCO
3, R-Cu and Q-Cu, but there are bands at 1552 cm
-1 and 1526 cm
-1 for R-Ag/Cu and Q-Ag/Cu, respectively, which can be assigned to nitro groups (N-O) stretching.
The aromatic rings evidenced by the 3072 cm-1 C-H aromatic stretching for RB21 and 3014 cm-1 for RB21-NaHCO3 are no longer observed to R-Cu, Q-Cu, R-Ag/Cu and Q-Ag/Cu. They may not be visible because they overlap with the broadband in that area, which extends to the 2900 cm-1 -2887 cm-1 aliphatic region in all spectra. This hypothesis is supported by the presence of vibrations at 1652 cm-1, 1641 cm-1, and 1625 cm-1 in the Q-Cu, R-Ag/Cu, and Q-Ag/Cu spectra, respectively, these vibrations being attributed to C-H aromatic bending vibrations. In the R-Cu spectrum, the band corresponding to the C-H bending is missing.
The aliphatic C-H stretching denoted by the medium intensity doublet bands at 2930 cm
-1, 2855 cm
-1 in the RB21 spectrum and 2919 cm
-1, 2849 cm
-1 in the RB21-NaHCO
3 (
Figure 5a), become very weak and single in the tested samples: 2944 cm
-1 (R-Cu), 2882 cm
-1 (Q-Cu) and 2887 cm
-1 (R-Ag/Cu and Q-Ag/Cu). In the fingerprint area, the corresponding bands for CH
3 asymmetric and CH
2 bending are indicated by the bands at 1493 cm
-1, 1463 cm
-1 in the RB21 spectrum, 1455 cm
-1 in the RB21-NaHCO
3; 1458 cm
-1 (R-Cu), 1422 cm
-1 (R-Ag/Cu) and 1494 cm
-1 (Q-Ag/Cu) and for CH
3 symmetric bending, they overlap with the sulfonic and sulfonamides groups for CH
3 symmetric and CH
2 bending are at 1301 cm
-1 (RB21), 1353 cm
-1 (RB21-NaHCO
3), 1390 cm
-1 (R-Cu) and 1396 cm
-1 (Q-Cu).
The Carbonate (CO
32-) groups are evidenced in the RB21-NaHCO
3 spectrum (
Figure 5a) by the 1748 cm
−1 and 1462 cm
−1 bands; the latter also indicating the lattice modes [
17]. The crystalline structure of the NaHCO
3 is denoted by the ascending spectrum baseline, which is also assigned to the carbonates as Mie scattering effect which is also observed for the R-Cu, Q-Cu, R-Ag/Cu and Q-Ag/Cu baselines of the spectra (
Figure 5b,c). In the R-Cu, Q-Cu and R-Ag/Cu spectra, the C=O group in carbonates is denoted by the bands at 1861 cm
-1, 1761 cm
-1, 1730 cm
-1 and 1720 cm
-1, while in the Q-Ag/Cu spectrum, the band is missing.
The sulfonic groups (SO3) in the chromogen and the sulfone groups (SO2) in the sulfonamides of the chromogen are denoted by the 1301 cm-1 in the RB21 spectrum and 1353 cm-1 in the RB21-NaHCO3 spectrum. In the R-Cu and Q-Cu spectra, the shifts of the bands at 1390 cm-1 and 1396 cm-1 (higher wavenumbers) can be assigned to transitions to sulfate groups (SO42-). In the R-Ag/Cu and Q-Ag/Cu spectra, the vibrations in the range of sulfonic and sulfone groups are missing and no vibrations in the range of groups with sulfur and oxygen are observed.
The changes in the RB21 dye IR spectrum (
Figure 5a) in its interaction with the NaHCO
3 denoted by the spectral bands for the RB21-NaHCO
3 are found as being highly enhanced when the aqueous solution of the two compounds enter in contact with the thin films of copper (R-Cu and Q-Cu samples) and with the Ag/Cu double layered thin films (R-Ag/Cu and Q-Ag/Cu). Both thin films induce chemical transformations of high aggressiveness on the RB21 dye, with decomposition results similar to each other and to those observed in the presence of the silver thin film presented earlier [
17]. The catalytic effect of thin copper films and the double-layered Ag/Cu films is proven the same as it was for the thin silver film [
17].
3.4. Profilometry
The topography of the surfaces of thin copper films and double silver/copper thin films produced by the PDL method was studied through the analysis of peak and valley profile measurements and the statistical calculation of characteristic parameters.The measurements were carried out under the following conditions:
| Scan Duration: |
30 s |
| Scan Length: |
3000 µm |
| Scan Resolution: |
0.333259 µm |
| Scan Type: |
Standard Scan |
| Stylus Force: |
10 mg |
| Stylus Scan Range: |
65.5 µm |
| Stylus Type: |
Radius: 12.5 µm |
The thickness of the thin films was evaluated by measuring the profile starting from their edge, being able to observe the step denoting the change in the height of the baseline of the profile. Thus, the copper thin film on area denoted as Cu(1) shows a first step of 75.55 nm in height, and after 1 cm in length in a horizontal direction, the base line rises up to 312.58 nm height along an upward slanting direction (
Figure 6a). This profile indicates that the deposition of the copper on the substrate is more concentrated to the middle of the substrate used for thin film deposition. On a different area, Cu(2), the copper thin film presents the same pattern of the baseline, only less difference between step1 (59.67 nm) and step2 (97.86 nm). Also, when compared the two areas on the copper thin film, Ra=150.78 nm which is the peaks’ arithmetical mean height of Cu(1) area is higher than Ra=60.48 nm of the Cu(2) area, despite the Rp=666.56 nm value of maximum peak height of Cu(1) area being less than Rp=806.30 nm of Cu(2) area (
Figure 6a,c). Regarding Ag/Cu thin film, the profiles of the two analyzed areas show a step of 226.93 nm corresponding to the edge of the thin film, following a horizontal direction of the baseline. The highest peaks of the surface profile of the double thin film Ag/Cu have the values Rp = 6899.45 nm for Ag/Cu(1) and Rp = 5218.49 nm for the area Ag/Cu(2). The highest peaks of the surface profile of the double thin film Ag/Cu have values of Rp = 6899.45 nm for Ag/Cu(1) and Rp = 5218.49 nm for the area Ag/Cu(2), which are approximately 10 times higher than those observed on the surface of the copper thin film. Also, the values of the arithmetic mean on the two areas on the surface of the thin Ag/Cu film are approximately 10 times higher, namely 986.75 nm on Ag/Cu (1) and 817.49 nm Ag/Cu (2).
To better understand the aspect of the roughness of the studied surfaces, peaks from the profile obtained through the measurements performed were analyzed. This way, proportions between the width and height (Width/Height ratio) of the peaks with values of 28.19 for the peak A on the Cu(1), 19.42 for the peak B on the Cu(2), 54.21 for the peak C on the Cu(2), 4.38 for the peak D of the Ag/Cu (1) and 5.37 for the peak D of the Ag/Cu (2) (
Figure 6b,d,e,g,i). The calculated Width/Height ratios are the smallest, because the selected peaks are the maximum peaks on the analyzed surface area. The information obtained consists of the fact that, in general, droplets have Width/Height ratios much smaller than these. The Width/Height ratio shows that the width of the droplets on the surface of the thin copper film is much greater than their height. Comparatively, the Ag/Cu thin films show droplets with Width/Height ratios 4 to 10 times smaller.
The parameters of the surface roughness of the measured areas on the two samples of thin films are presented in
Table 2.
The roughness of all the analyzed areas of the samples – Cu(1), Cu(2), Ag/Cu(1) and Ag/Cu(2) – is characterized by a positive skewness (Rsk>0). In the right-skewed (positively skewed) distribution most data points cluster on the left, with a long tail to the right. The tails represent the probability or frequency of values that are extremely high or low compared to the mean. This means a few unusually high values pull the mean upward and it indicates that a larger number of peaks is below the mean line of the roughness.
Kurtoisis as a measure of the tailedness of a distribution shows how often outliers occur. The kurtosis variation for the Cu thin film surface roughness from Rku =2.76 to 13.48 on the two measured areas indicates non uniform distribution of the peaks and pits (extreme values) from a platykurtic (thin tails) distribution (Rku < 3) to a high leptokurtic (fat tails) distribution (Rku >> 3). That means that on some areas the Cu thin film presents fewer peaks and valleys than a normal distribution and more peaks and valleys than the normal distribution on other areas. The kurtosis of 6.08 and 5.34 of the two areas measured on the Ag/Cu thin film shows a leptokurtic (thin tails) distribution of the extreme values [
25].
The profile analysis of the two thin films complements the analysis of surface morphology using the SEM technique, providing information about the surface relief and their ability to retain the aqueous solutions with which they are treated. The results of the roughness measurements are consistent with the estimates obtained by simulating the laser irradiation process of the two targets, the production of the fluid phase, and its quality in the ablation process.
Subsequently, the analysis of the profile of the two thin films after treatment aims to evaluate how the aqueous RB21-NaHCO3 solution produces corrosion on its surface.
The profiles after treatement with RB21-NaHCO
3 (
Figure 7) show a high degradation of the Cu thin film and of the Ag/Cu thin film. The corrosive effect of the treatment solution on the thin layers of Cu and the double-layer Ag/Cu is highlighted by the baseline that slopes downward to the right for the thin Cu film and by a deep valley in the case of Ag/Cu, where residual peaks are observed, probably from the thin layer prior to treatment with the solution. The changes in the height of the maximum peak after treatment at the value of Rp=2070.93 nm for Cu-RB21-NaHCO3, which shows an increase of approximately three times compared to that determined on the copper thin film, and at the value of Rp=1782.54 nm for Ag/Cu-RB21-NaHCO3, which shows a decrease of almost 4 times compared to the double thin film, are due to the different action of the RB21-NaHCO3 solution on the surface of copper and silver, differently modifying their roughness. Regarding the mean values, a significant change is observed for the Cu-RB21-NaHCO3 sample, for which the arithmetic mean of the height Ra=271.80 nm represents an increase of almost 2 times and more than 4 times compared to the initial thin film, while the mean value Ra=924.59 nm of the Ag/Cu-RB21-NaHCO3 sample remained almost the same as the initial thin film. It can thus be deduced that corrosion occurred in the pits (valleys) also due to changes in the size of the pits. In this regard, increases in the maximum pit height value were recorded, 35 times higher for the Cu-RB21-NaHCO3 sample, for which Rv=730.67 nm, compared to the initial thin copper film, for which Rv=20.61 nm, while for the Ag/Cu-RB21-NaHCO3 sample, the maximum pit height value Rv=2311.69 nm is 70 times higher than Rv=33 nm recorded for the initial thin double Ag/Cu film.
By calculating the Width/Height ratios for the highest peak in the roughness profile for the two thin layers after treatment with the aqueous RB21-NaHCO3 solution, values of 10.50 are obtained for the Cu-RB21-NaHCO3 sample and 12.75 for the Ag/Cu-RB21-NaHCO3 sample. These results show a decrease of 2 to 5 times in the Width/Height ratio for Cu-RB21-NaHCO3 and 2 times increase for Ag/Cu-RB21-NaHCO3.
The left-skewed (negatively skewed) obtained for the roughness of the two films after treatment with RB21 – NaHCO
3 aqueous solution (
Table 3) shows that most data points cluster on the right, with a long tail to the left. A few unusually low values pull the mean downward. The negative skewness (Rsk < 0) resulted for the Cu-RB21-NaHCO
3 sample and for the Ag/Cu-RB21-NaHCO
3 sample shows that both thin films were subjected to notable corrosive degradation, which led to the arrangement of several peaks below the median line, in contrast to the positive skewness noticed for the thin films before their treatment with the RB21-NaHCO
3 solution.
The kurtosis value of Rku =5.65 for the Cu-RB21-NaHCO
3 being greater than 3 indicates that the distribution of the extreme values is leptokurtic (fat tails) [
19], while Rk=1.94 for the Ag/Cu-RB21-NaHCO
3 shows that the distribution of the extreme values is platykurtic (thin tails) [
25].
Comparing the profilometric results of the surface roughness of the studied films for the initial thin films and after their treatment with the RB21-NaHCO3 solution, it can be observed that corrosion has reduced the extreme values, also leading to the lowering of the peaks to a level below the mean line. Also, the fact that the corrosion process led to the increase in the depth of the pits on the surface of the thin films is due to the accumulation and longer retention of the solution in the pits, which allows the chemical process to take place for a longer period of time.
3.6. UV-Vis Spectroscopy
Figure 9 shows the UV-Vis spectra which indicate an advanced decomposition of the dye compared to the RB21-NaHCO3 control sample. In the visible range, a shift of the peaks from 619 nm and 664 nm to 625 nm and 665 nm is observed in the spectra of the Cu-RB21-NaHCO
3 and respectively Ag/Cu-RB21-NaHCO
3 samples, as well as a reduction in intensity. These spectral changes in the visible range show a shift towards higher wavelength values (bathochromic effect), the color of the absorbing material (the solution resulting after the interaction of RB21-NaHCO
3 with the surface of thin films of Cu and Ag/Cu). The decrease in intensity reflects a hypochromic shift. Both effects indicate changes at the level of the chromophore –SO
2NH–Benz.–SO
2CH
2CH
2OH of the RB21 dye, which is consistent with the results of the FTIR analysis.
Changes in conjugation or the conjugation effect are observed in the UV domain by the shifting of the peaks from 261 nm and 337 nm in the spectrum of the RB21-NaHCO3 control sample to 272 nm and 342 nm in the spectra of the solutions resulting after interaction with the two thin films. As with the visible range, a drastic decrease in intensity is observed. The conjugation effect indicates changes within the chromogen of the RB21 dye.
The bathochromic effect observed for the peaks in the UV range indicates changes in the auxochromes of the dye molecule (methyl, hydroxyl, amino groups) regarding their interaction with the π electron conjugation system.