Influence of Irradiation Parameters on Structure and Properties of Oak Wood Surface Engraved with a CO2 Laser

The work investigates the effects of CO2 laser parameters (laser power and raster density) on wood mass loss in oak wood and impacts on its morphology, chemical structure, and surface properties (colour and hydrophilicity). The energy amount supplied onto the wood surface with a laser beam under different combinations of the irradiation parameters was expressed through a single variable—total irradiation dose. The mass loss was confirmed as linear-dependent on the irradiation dose. With the mass reduction, the roughness was enhanced. The roughness parameters Ra and Rz increased linearly with the mass loss associated with the increasing irradiation dose. The FTIR (Fourier transform infrared spectroscopy) spectroscopy also detected chemical changes in the main wood components, influencing primarily the wood colour space. Conspicuous discolouration of the engraved wood surface was observed, occurring just at the minimum laser power and raster density. The additional increasing of laser parameters caused a novel colour compared to the original one. The detected dependence of wood discolouration on the total irradiation dose enables us to perform targeted discolouration of the oak wood. The engraved surfaces manifested significantly better wettability with standard liquids, both polar and non-polar, and higher surface energy values. This guarantees appropriate adhesion of film-forming materials to wood. Identification of the changes in wood surface structure and properties, induced by specific CO2 laser-treatments, is important for obtaining targeted discolouration of the wood surface as well as for the gluing or finishing of the surfaces treated in this way.

The energy concentrated in a laser beam and supplied onto a specific spot on the engraved wood surface is converted to heat. This heat causes thermal degradation of wood surface structures in the area of the heat-affected zone. The high energy values cause the wood surface layer to sublimate, which is accompanied by the formation of a carbonized or melted layer. The wood surface treated in this way exhibits changes in its chemistry and structure, reflected subsequently in the changes in morphology, colour, and wetting performance with liquids [18,[21][22][23][24]. This energy can be controlled by the laser power, the rate of laser head movement, focal distance, and raster density [11,12,16,22,[25][26][27][28]. Apart from the energy supplied amount and its concentration, the thickness of the heat-affected zone depends considerably on the wood species concerned [11,12,19,22,23]. the wood surface. The differences between the authors can be explained by using lasers with different technical parameters based on the different methodical approaches. Despite this fact, the results clearly show that the appropriate adjustment of CO 2 laser radiation parameters can be a means for targeted wood surface modification, in accordance with the current requirements. For this purpose, it is necessary to study the whole range of changes in wood surface properties, using the same methods.
The aim of this work was to evaluate comprehensively the influence of infrared laser radiation generated with a CO 2 laser on the surface structure and properties of oak wood. This aim comprehended the following partial tasks: to determine the total amount of energy supplied with a laser beam onto the wood surface under varying CO 2 laser parameters (laser power, raster density); -to assess the impact of the total radiation dose on the wood mass loss; -to assess the impact of the specific method of surface treatment on the wood surface chemical structure and morphology quantified through roughness parameters; -to inspect colour changes depending on different laser parameters and to identify how the total energy supplied influenced the changes in the colour coordinates; -to find out the influence of engraving on the wetting process and on the surface free energy.

Experimental Material and Irradiation Parameters
The experiments were carried out on specimens prepared from ring-porous oak wood. The specimen dimensions were 50 mm × 100 mm × 15 mm (width × length × thickness). Prior to the grinding, the surfaces were sanded with a sandpaper with a grain size of P180. The irradiated (engraved) surfaces were radial, sized 50 mm × 100 mm (Figure 1a).
During the engraving, the specimens were placed under the focus of a CO 2 focusing lens of laser equipment CM-1309 (Shenzhen Reliable Laser Tech, Shenzhen, China). The specimen distance from the lens was 17 mm. The laser head moved over the specimen surface parallel to the grain, at a constant speed of 350 mm·s −1 (Figure 1b,c). The radiation energy varied with varying laser power and raster density (number of paths per one millimetre of width). Three specimen sets were prepared. The first was engraved at a laser power of 8% (derived based on the maximal laser power, representing the resonator output value of 137.5 W), the second at 12% laser power, and the third at 16% power. The laser power was measured with measuring equipment FieldMaxII-TOP (Coherent, Wilsonville, OR, USA). The equipment sensor PM150-50C placed at the level of the specimen surface measured the power of laser beams perpendicular to the surface. The analysis of the whole summary of the published results shows that the energy absorbed into the surface during CO2 laser treatment induces chemical changes in the main components of wood, in its morphology, colour, and hydrophilicity or hydrophobicity of the wood surface. The differences between the authors can be explained by using lasers with different technical parameters based on the different methodical approaches. Despite this fact, the results clearly show that the appropriate adjustment of CO2 laser radiation parameters can be a means for targeted wood surface modification, in accordance with the current requirements. For this purpose, it is necessary to study the whole range of changes in wood surface properties, using the same methods.
The aim of this work was to evaluate comprehensively the influence of infrared laser radiation generated with a CO2 laser on the surface structure and properties of oak wood. This aim comprehended the following partial tasks: to determine the total amount of energy supplied with a laser beam onto the wood surface under varying CO2 laser parameters (laser power, raster density); -to assess the impact of the total radiation dose on the wood mass loss; -to assess the impact of the specific method of surface treatment on the wood surface chemical structure and morphology quantified through roughness parameters; -to inspect colour changes depending on different laser parameters and to identify how the total energy supplied influenced the changes in the colour coordinates; -to find out the influence of engraving on the wetting process and on the surface free energy.

Experimental Material and Irradiation Parameters
The experiments were carried out on specimens prepared from ring-porous oak wood. The specimen dimensions were 50 mm × 100 mm × 15 mm (width × length × thickness). Prior to the grinding, the surfaces were sanded with a sandpaper with a grain size of P180. The irradiated (engraved) surfaces were radial, sized 50 mm × 100 mm ( Figure  1a).
During the engraving, the specimens were placed under the focus of a CO2 focusing lens of laser equipment CM-1309 (Shenzhen Reliable Laser Tech, Shenzhen, China). The specimen distance from the lens was 17 mm. The laser head moved over the specimen surface parallel to the grain, at a constant speed of 350 mm·s −1 (Figure 1b,c). The radiation energy varied with varying laser power and raster density (number of paths per one millimetre of width). Three specimen sets were prepared. The first was engraved at a laser power of 8% (derived based on the maximal laser power, representing the resonator output value of 137.5 W), the second at 12% laser power, and the third at 16% power. The laser power was measured with measuring equipment FieldMaxII-TOP (Coherent, Wilsonville, OR, USA). The equipment sensor PM150-50C placed at the level of the specimen surface measured the power of laser beams perpendicular to the surface. four specimens, plus four referential ones (0 mm −1 ). Under the experimental conditions, all the specimens were irradiated uniformly over their length and width ( Figure 1). The laser parameters (laser power, head movement rate, focal distance and raster density) were decided based on our experience with previous measurements (Kúdela et al. 2019(Kúdela et al. , 2021) and on the knowledge collected from the literature referred to in the Introduction.
The amount of the energy supplied was expressed as the amount of irradiation dose H'. The total irradiation dose per unit area of irradiated specimens was determined in accordance with Kubovský et al. [32]. The authors calculated the irradiation dose for one route of the laser beam according to their equation: where P e is the effective laser power on the specimen surface, τ is irradiation time during one path (ratio between the specimen dimension x and speed v), and A is area irradiated at one laser beam path. The irradiation dose H' related to a unit area of the irradiated specimen was subsequently calculated according to the equation: where n is the raster density, and y is the width of the track left by the laser beam after one path on the wood surface. The track width was measured with the aid of a camera Leica EC3 (Leica Microsystems, Heerbrugg, Switzerland).

Determining the Weight Ratio of the Wood Mass Loss
For determining the ratio of wood mass loss after engraving the wood surface with a CO 2 laser, new specimens were prepared from the same experimental material. The specimens were dried out to a zero-moisture content at a temperature of 103 ± 2 • C. The dry specimens were removed from the drying chamber and placed into an exicator with silica gel. After the cooling, the specimens were weighed on a laboratory scale with a precision of 0.001 g (mass m 01 ). The dried and weighed specimens were engraved at specified laser powers and raster densities. After this process, the specimens were weighed again, the result of which was the mass m 02 . Then, the equation was used for calculating the weight ration for the wooden mass loss.

Evaluation of Surface Morphology
Wood surface morphology after the laser treatment was evaluated from the viewpoints of anatomy and physics. The roughness profiles were measured with a profilometer Surfcom 130A (Carl Zeiss, Oberkochen, Germany) consisting of a measuring unit and an evaluation unit. The changes to the wood surface morphology induced by CO 2 laser engraving were evaluated through the following roughness parameters: Ra (mean arithmetic deviation), Rz (the maximum peak height plus the maximum depression depth within the cut-off, or sampling length), Rt (the maximum peak height plus the maximum depression depth within the entire evaluation length) and RSm (mean distance between the trenches-arithmetic mean calculated from distances between the profile unevenness within the sampling length).
The roughness was measured on the irradiated radial surfaces, parallel with and perpendicular to the grain. The roughness evaluation started by filtering off the waviness from the basic profile, and then the roughness curve was transferred onto the basic line. The entire transverse length consisted of run-up segment l r , five sampling length (cutoff) segments λ c , and the over travel segment; l p . The basic lines were chosen from the interval 0.025-8 mm, based on the preliminary measured roughness parameters Ra and Rz. For studying the structure of the engraved oak wood surfaces, a light microscope Leica MZ 9.5, camera Leica EC3 was used, and a digital microscope Keyence VHX 7000 (Keyence International, Mechelen, Belgium). Microscopic slides of the transverse cuts were prepared following the methods described in [12].

Detecting Chemical Changes after Engraving
Chemical changes were inspected in specimens irradiated at all the tree laser power and raster density values 10, 20, and 30 mm −1 . For this purpose, FTIR Analysis (Fourier Transform-Infrared Spectroscopy) was used. FTIR spectra of the engraved wood surfaces were recorded on a Nicolet iS10 FTIR spectrometer (Thermo Fisher Scientific, Waltham, MA, USA), equipped with Smart iTR using an attenuated total reflectance (ATR) sampling accessory attached to a diamond crystal. The spectra were acquired by accumulating 64 scans at a spectral resolution of 4 cm −1 in an absorbance mode from 4000 to 650 cm −1 standardised using the baseline method. The obtained results were analysed using OMNIC 9.0 software (Thermo Fisher Scientific, Waltham, MA, USA). Measurements were performed on four replicates per a sample. Each specimen was measured on the radial surface.

Colour Measurement
The colorimetric values L*, a* and b* on the referential and engraved specimens were measured with a spectrophotometer Spectro-guide 45/0 gloss (BYK-GARDNER GmbH, Geretsried, Germany). The measurements were taken at ten spots per one specimen. The colour differences ∆L*, ∆a*, ∆b* under different irradiation conditions and the total colour difference ∆E were determined according to the following equations: where the index "1" represents the colour coordinate value in the referential specimen (not engraved) and the index "2" indicates the value of the colour coordinate of the wood specimen irradiated with the CO 2 laser.

Wood Surface Wetting with Liquids and Determining the Surface Free Energy
Wood wetting with standard liquids was performed, together with measuring the contact angles during the entire wetting process until the complete drop soaking into the substrate were realised with the aid of a goniometer Krüss DSA30 Standard (Krüss, Hamburg, Germany). The wetting process as such was evaluated using a software package DSA3 (Krüss, Hamburg, Germany). Two testing liquids with different polarities were used-redistilled water and diiodomethane. The reasons for using these two liquids follow from [40].
The wood wetting was realised by applying a drop of the relevant liquid with a volume of 0.0018 mL onto the substrate. After the drop contact with the wood surface, a camera scanned the time history of the drop shape in the fibre direction, from the first contact up to the complete soaking. The scanning frequency was adjusted according to the wetting process duration. The drop shape analysis and determining the contact angle were performed by a circle method.
Immediately after applying the drop onto the surface, the contact angle value θ 0 was measured. Based on the parameter d (drop width) variation, the moment on the advancing contact angle conversion into the receding one was identified. The contact angle measured at this moment was defined as the "equilibrium" contact angle-θ e . The contact angle values were measured at twelve different spots on each specimen. The values of contact angles θ 0 and θ e , were used, following the methods described in [41], for the calculation of the contact angle corresponding to an ideally smooth surface θ w . This angle was subsequently used for determining the surface free energy γ S and its disperse and polar shares γ S d and γ S p . The disperse share was obtained from the wood wetting with diiodomethane, the polar share with water. The total surface free energy was defined as the sum of the polar and non-polar components, according to [40].

Results and Discussion
The data concerning the wood mass loss, chemical changes, and the changes in the morphology and performance of wood surface engraved at variable laser power and raster density values are presented and analysed in the following sections.

Determining the Total Irradiation Dose
The values H of the irradiation dose delivered over a unit area were calculated according to Equation (2), with parameters expressing conditions for the relevant CO 2 laser irradiation. These results are listed in Table 1. The total irradiation dose ranged from 6 to 185 J·cm −2 , with the lowest value obtained at the minimum (8%) laser power and the lowest raster density (2 mm −1 ), and the highest at the respective values of 16% and 30 mm −1 . The increasing irradiation dose was attained by a linear increase in the weight rate of the wood mass loss ( Figure 2). Quantitatively similar wood mass loss could also be attained by adjusting the focal distance of the lens and the rate of the laser head movement [16,42]. It was confirmed that all the irradiation parameters impacting s the radiation energy significantly (laser power, head shifting rate, focal distance, and raster density) could be substituted with a single variable-irradiation dose H . Non-uniform mass loss in oak wood was accordingly also reflected in the wood morphology. Immediately after applying the drop onto the surface, the contact angle value θ0 was measured. Based on the parameter d (drop width) variation, the moment on the advancing contact angle conversion into the receding one was identified. The contact angle measured at this moment was defined as the "equilibrium" contact angle-θe. The contact angle values were measured at twelve different spots on each specimen.
The values of contact angles θ0 and θe, were used, following the methods described in [41], for the calculation of the contact angle corresponding to an ideally smooth surface θw. This angle was subsequently used for determining the surface free energy S and its disperse and polar shares S d and S p . The disperse share was obtained from the wood wetting with diiodomethane, the polar share with water. The total surface free energy was defined as the sum of the polar and non-polar components, according to [40].

Results and Discussion
The data concerning the wood mass loss, chemical changes, and the changes in the morphology and performance of wood surface engraved at variable laser power and raster density values are presented and analysed in the following sections.

Determining the Total Irradiation Dose
The values H´ of the irradiation dose delivered over a unit area were calculated according to Equation (2), with parameters expressing conditions for the relevant CO2 laser irradiation. These results are listed in Table 1. The total irradiation dose ranged from 6 to 185 J•cm −2 , with the lowest value obtained at the minimum (8%) laser power and the lowest raster density (2 mm −1 ), and the highest at the respective values of 16% and 30 mm −1 . The increasing irradiation dose was attained by a linear increase in the weight rate of the wood mass loss ( Figure 2). Quantitatively similar wood mass loss could also be attained by adjusting the focal distance of the lens and the rate of the laser head movement [16,42]. It was confirmed that all the irradiation parameters impacting s the radiation energy significantly (laser power, head shifting rate, focal distance, and raster density) could be substituted with a single variable-irradiation dose H´. Non-uniform mass loss in oak wood was accordingly also reflected in the wood morphology.

Morphology of Wood Surfaces Engraved with CO 2 Laser
The surface morphology of laser-engraved oak wood was evaluated based on experimentally obtained roughness profiles, parallel with and perpendicular to the grain, in the radial direction. The basic statistical characteristics for the roughness parameters Ra, Rz, Rt, and RSm for the relevant raster density range, power values, and the two anatomical directions are in Table 2. The results of a three-way variance analysis confirmed a significant influence of all three evaluated factors (raster density, laser power, and anatomical direction) and their interactions on the evaluated roughness variables. In the case of parameters Ra, Rz, and RSm, the number of measurements n for each variant was 60. In the case of Rt, n = 12.
The values of all the roughness parameters in the referential specimens were lower than in the engraved specimens, in both anatomical directions, parallel with and perpendicular to the grain. In the referential specimens, significantly higher roughness parameter values were measured perpendicular to the grain, consistently with the orientation of cell elements.
In the engraved specimens, the roughness parameter values Ra, Rz, and Rt parallel to the grain significantly increased with increasing laser power, and mainly with increasing raster density. At the maximum laser power and the maximum raster density, the roughness parameter increase was several-fold. The roughness parameter values RSm, parallel to the grain, increased with increasing raster density as far as 5 mm −1 ; there followed a moderate decrease ( Figure 3). In the case of parameters Ra, Rz, and RSm, the number of measurements n for each variant was 60. In the case of Rt, n = 12.
The values of all the roughness parameters in the referential specimens were lower than in the engraved specimens, in both anatomical directions, parallel with and perpendicular to the grain. In the referential specimens, significantly higher roughness parameter values were measured perpendicular to the grain, consistently with the orientation of cell elements.
In the engraved specimens, the roughness parameter values Ra, Rz, and Rt parallel to the grain significantly increased with increasing laser power, and mainly with increasing raster density. At the maximum laser power and the maximum raster density, the roughness parameter increase was several-fold. The roughness parameter values RSm, parallel to the grain, increased with increasing raster density as far as 5 mm −1 ; there followed a moderate decrease (Figure 3).  Perpendicular to the grain, the roughness was significantly higher over the whole engraving range. With increasing raster density, the parameter values Ra, Rz, and Rt increased much more than in the longitudinal direction (Figure 3). At the maximum raster density, the Ra values were by order higher compared to the referential ones. The parameter RSm exhibited the most variability compared to the others. The raster-density-dependent change in parameter RSm perpendicular to the grain differed from the corresponding change in the parallel direction, in quantity equally as in quality (Figure 3).
Energy concentrated in the laser beam and supplied on a specific spot was converted to heat. The checking of temperature with a thermo-camera, revealed that the wood surface temperature at the moment of contact between wood and beam fluctuated close to the upper measuring threshold of the camera (1000 • C) just at an 8% laser power. Such extreme heat, concentrated on the wood surface within the laser beam with a tiny diameter, caused immediate demolition and sublimation of the thin surface wood layer. Under the same conditions, the thickness of the sublimated layer was not uniform over the irradiated area, due to the wood heterogeneous structure. This was primarily evident in measuring the roughness perpendicular to the grain course. As the direction was radial, the differences observed were a consequence mainly of the differences in the qualitative and quantitative presence of cell elements between the early and late wood.
Experimental measurements of roughness and the microscopical observations showed that, apart from the energy amount and concentration, the thickness of the sublimated layer was to a large extent affected by the differences in the structure and properties between the earlywood and the latewood, which was also confirmed for other types of wood [12]. Nonuniform destruction of the wood mass within growth rings, due to the density differences between the early and the late wood, was mainly reflected in more pronounced trenches in the early wood and partly also in distances between the profile unevenness. The altered irradiated wood surface morphology is shown in Figures 4 and 5a. Perpendicular to the grain, the roughness was significantly higher over the whole engraving range. With increasing raster density, the parameter values Ra, Rz, and Rt increased much more than in the longitudinal direction (Figure 3). At the maximum raster density, the Ra values were by order higher compared to the referential ones. The parameter RSm exhibited the most variability compared to the others. The raster-density-dependent change in parameter RSm perpendicular to the grain differed from the corresponding change in the parallel direction, in quantity equally as in quality (Figure 3).
Energy concentrated in the laser beam and supplied on a specific spot was converted to heat. The checking of temperature with a thermo-camera, revealed that the wood surface temperature at the moment of contact between wood and beam fluctuated close to the upper measuring threshold of the camera (1000 °C) just at an 8% laser power. Such extreme heat, concentrated on the wood surface within the laser beam with a tiny diameter, caused immediate demolition and sublimation of the thin surface wood layer. Under the same conditions, the thickness of the sublimated layer was not uniform over the irradiated area, due to the wood heterogeneous structure. This was primarily evident in measuring the roughness perpendicular to the grain course. As the direction was radial, the differences observed were a consequence mainly of the differences in the qualitative and quantitative presence of cell elements between the early and late wood.
Experimental measurements of roughness and the microscopical observations showed that, apart from the energy amount and concentration, the thickness of the sublimated layer was to a large extent affected by the differences in the structure and properties between the earlywood and the latewood, which was also confirmed for other types of wood [12]. Non-uniform destruction of the wood mass within growth rings, due to the density differences between the early and the late wood, was mainly reflected in more pronounced trenches in the early wood and partly also in distances between the profile unevenness. The altered irradiated wood surface morphology is shown in Figures 4 and  5a. Oak wood is ring-porous, with early wood vessels having big lumens and thin cell walls, and with libriform fibres exhibiting small lumen diameters and thick cell walls. The roughness profiles measured perpendicular to the fibre course in the radial direction were Oak wood is ring-porous, with early wood vessels having big lumens and thin cell walls, and with libriform fibres exhibiting small lumen diameters and thick cell walls. The roughness profiles measured perpendicular to the fibre course in the radial direction were affected by cell wall degradation, especially in early wood vessels. This resulted in forming deeper depressions in such spots. The paper [12] suggests that the same is true in the case of irradiated tangential surfaces, which is evident on the microscopic slides in Figure 5b-d. If the irradiated oak wood surfaces exhibit dominant tangential shares, the heat supplied onto the surface during engraving penetrates deeper into the wood through its pith rays. The pith rays are composed of parenchyma thin-walled cells that under heat impact degrade easily, up to carbonisation. The carbonised cells in micro-cuts are evidenced as black strips (Figure 5b-d).
Non-uniform engraving of the oak wood surface was primarily responded with higher values of parameters Ra and Rz, increasing linearly with the wood mass loss rate from the oak wood surface ( Figure 6). The results also indicate that the roughness variation was more affected by the raster density than by the laser power, which has also been confirmed by [11,27]. As for the impacts on the wood surface geometry, in this case, all the engraving parameters are possible to substitute with a single variable-the total irradiation dose. With an increasing irradiation dose, the roughness parameters increased linearly. After the engraving, not only the wood surface geometry was changed but also the wood surface chemical structure. affected by cell wall degradation, especially in early wood vessels. This resulted in forming deeper depressions in such spots. The paper [12] suggests that the same is true in the case of irradiated tangential surfaces, which is evident on the microscopic slides in Figure  5b-d. If the irradiated oak wood surfaces exhibit dominant tangential shares, the heat supplied onto the surface during engraving penetrates deeper into the wood through its pith rays. The pith rays are composed of parenchyma thin-walled cells that under heat impact degrade easily, up to carbonisation. The carbonised cells in micro-cuts are evidenced as black strips (Figure 5b-d).
Non-uniform engraving of the oak wood surface was primarily responded with higher values of parameters Ra and Rz, increasing linearly with the wood mass loss rate from the oak wood surface (Figure 6). The results also indicate that the roughness variation was more affected by the raster density than by the laser power, which has also been confirmed by [11,27]. As for the impacts on the wood surface geometry, in this case, all the engraving parameters are possible to substitute with a single variable-the total irradiation dose. With an increasing irradiation dose, the roughness parameters increased linearly. After the engraving, not only the wood surface geometry was changed but also the wood surface chemical structure.   affected by cell wall degradation, especially in early wood vessels. This resulted in forming deeper depressions in such spots. The paper [12] suggests that the same is true in the case of irradiated tangential surfaces, which is evident on the microscopic slides in Figure  5b-d. If the irradiated oak wood surfaces exhibit dominant tangential shares, the heat supplied onto the surface during engraving penetrates deeper into the wood through its pith rays. The pith rays are composed of parenchyma thin-walled cells that under heat impact degrade easily, up to carbonisation. The carbonised cells in micro-cuts are evidenced as black strips (Figure 5b-d).
Non-uniform engraving of the oak wood surface was primarily responded with higher values of parameters Ra and Rz, increasing linearly with the wood mass loss rate from the oak wood surface ( Figure 6). The results also indicate that the roughness variation was more affected by the raster density than by the laser power, which has also been confirmed by [11,27]. As for the impacts on the wood surface geometry, in this case, all the engraving parameters are possible to substitute with a single variable-the total irradiation dose. With an increasing irradiation dose, the roughness parameters increased linearly. After the engraving, not only the wood surface geometry was changed but also the wood surface chemical structure.   The engraved wood surfaces, especially those subjected to higher irradiation doses, exhibited a carbonised layer (Figures 4 and 5) characterised by poor stability, and also a melted layer, especially on the late wood surface (Figure 7). These changes were induced mainly by the chemical changes in the surface layers of the studied oak wood. The engraved wood surfaces, especially those subjected to higher irradiation doses, exhibited a carbonised layer (Figures 4 and 5) characterised by poor stability, and also a melted layer, especially on the late wood surface (Figure 7). These changes were induced mainly by the chemical changes in the surface layers of the studied oak wood.

Changes in the FTIR Spectra
FTIR spectra were measured on the native wood surface and the surface engraved with CO2 laser (for used raster densities of 10, 20, and 30 mm −1 ). All the spectra were evaluated altogether, as the differences between them were small. We can see (Figure 8) that the variance in the band from 3100 to 3600 cm −1 was minimum. This very broad peak has been assigned to the O−H vibrations in the lignin structure, in cellulose fibres, and in hemicelluloses [43]. The changes in the band absorbance with a peak at 3350-3360 cm −1 could be observed at all three laser power values used. Increasing absorbance was especially evident for lower raster density values; conversely, at higher raster densities, the absorbance decreased. This phenomenon may be backed-up with a higher rate of splitting the O−H bonds in the molecules of the water bound to the wood, due to higher amounts of the energy supplied [44]. The interval from 2900 to 2850 cm −1 (symmetric and asymmetric C−H stretching vibrations in aliphatic compounds) [45] exhibits similar effects, however, with lower variance in the absorbance density values.

Changes in the FTIR Spectra
FTIR spectra were measured on the native wood surface and the surface engraved with CO 2 laser (for used raster densities of 10, 20, and 30 mm −1 ). All the spectra were evaluated altogether, as the differences between them were small. We can see (Figure 8) that the variance in the band from 3100 to 3600 cm −1 was minimum. This very broad peak has been assigned to the O−H vibrations in the lignin structure, in cellulose fibres, and in hemicelluloses [43]. The changes in the band absorbance with a peak at 3350-3360 cm −1 could be observed at all three laser power values used. Increasing absorbance was especially evident for lower raster density values; conversely, at higher raster densities, the absorbance decreased. This phenomenon may be backed-up with a higher rate of splitting the O−H bonds in the molecules of the water bound to the wood, due to higher amounts of the energy supplied [44]. The interval from 2900 to 2850 cm −1 (symmetric and asymmetric C−H stretching vibrations in aliphatic compounds) [45] exhibits similar effects, however, with lower variance in the absorbance density values. In the range from 1800 to 800 cm −1 (the so-called fingerprint region), bands assigned to stretching and deformation vibrations of all wood components, more noteworthy changes in absorbance were recorded. The specimens treated at the lowest laser power (8%) displayed enhanced absorbance within the band of 1730 cm −1 (C=O stretching in unconjugated carbonyl groups), depending on the raster density used in engraving ( Figure  8). At the two higher laser power values (12% and 16%), the absorbance of this band increased with increasing raster density (except for the highest density value, at which a moderate decrease was observed). These changes point to the changes in several functional groups in lignin and in hemicelluloses (carbonyls, aldehydes, ketones, and carboxylic acids) [46,47]. Increasing absorbance in the band assigned to the non-conjugated C=O groups indicates increasing amounts of acetyl and carboxyl groups in the lignin and in polysaccharides [48]. Reduction in this band could be, in our case, the result of hemicellulose deacetylation during the exposure of the wood surface structures [24].
Heat-induced processes are associated with the degradation of carbonyl groups in lignin and in hemicelluloses where the cleavage of bonds C=O affects the changes in the In the range from 1800 to 800 cm −1 (the so-called fingerprint region), bands assigned to stretching and deformation vibrations of all wood components, more noteworthy changes in absorbance were recorded. The specimens treated at the lowest laser power (8%) displayed enhanced absorbance within the band of 1730 cm −1 (C=O stretching in unconjugated carbonyl groups), depending on the raster density used in engraving (Figure 8). At the two higher laser power values (12% and 16%), the absorbance of this band increased with increasing raster density (except for the highest density value, at which a moderate decrease was observed). These changes point to the changes in several functional groups in lignin and in hemicelluloses (carbonyls, aldehydes, ketones, and carboxylic acids) [46,47]. Increasing absorbance in the band assigned to the non-conjugated C=O groups indicates increasing amounts of acetyl and carboxyl groups in the lignin and in polysaccharides [48]. Reduction in this band could be, in our case, the result of hemicellulose deacetylation during the exposure of the wood surface structures [24].
Heat-induced processes are associated with the degradation of carbonyl groups in lignin and in hemicelluloses where the cleavage of bonds C=O affects the changes in the chromophores content. The chromophores are structures determining the wood surface colouring (this is also evident from the values related to the discolouration).
The band at 1600 cm −1 (C=C stretching vibration conjugated with an aromatic ring in lignin) practically copied the performance of the band assigned to the carbonyl groups. The absorbance of the band around 1500 cm −1 (C=C stretching vibration and aromatic skeleton vibration) showed only small changes. This indicates decrease in amount of methoxyl groups confirming the reduction in the lignin content [49,50]. At higher temperatures, condensation reactions in lignin are possible [51,52]. Moderate changes, primarily dependent on the raster density were also recorded for the bands 1460 cm −1 (asymmetric CH 3 bonding in methoxyl groups in lignin), 1370 cm −1 (symmetric and asymmetric CH 3 bonding), 1320 cm −1 (C−O vibration in syringyl derivatives), and 1230 cm −1 (C−O stretching vibration in xylan and syringyl ring) assigned to the lignin and hemicelluloses [53][54][55][56]. The decrease in these bands has confirmed the supposition about bond-cleavage-caused degradation processes in lignin and the subsequent structural decomposition [57]. Contrarily to the hitherto observed trends, we observed a permanent decrease for the band 1030 cm −1 (C−O deformation vibrations in cellulose) (Figure 8). The absorbance on this band exhibited a permanent decrease indicating the degradation processes in the cellulose [58].

Discolouration of Engraved Surfaces
The interactions at the phase boundary wood-laser beam are complex. These interactions not only influence the changes in the wood surface structure and chemistry, but they also impact wood discolouration. The results of the two-way variance analysis confirmed the important effects of the two relevant laser parameters (laser power and raster density) on the colour variation in the engraved surfaces. The basic statistical characteristics of the colour coordinates L*, a*, and b* are listed in Table 3. The highest average lightness L* was measured in the referential specimens (L* = 68). With increasing raster density, at the given power values, the lightness decreased significantly down to the final values from 26.6 to 21.8, representing the lightness reduction by 61-68%. The colour coordinates a* and b* varied too. The differences in the individual colour coordinates ∆L*, ∆a*, and ∆b*, together with the total colour difference ∆E*, are in Figure 9.
The coordinate a* values were increasing with increasing raster density until 5 mm −1 ; after this raster value, a decrease by 25 to 33 % was recorded compared to the referential specimens (corresponding to the laser power). With increasing laser power and raster density, the coordinate a* exhibited a progressive loss in red saturation, because of a moderate shift towards an achromatic colour. The coordinate b* values decreased proportionally, over the entire raster range, at all the laser power values. The final decrease in this coordinate was from 41 to 63%.
As for the discussed engraving parameters, the changes in colour coordinates were primarily a consequence of the varying the raster density. The effect of laser power has also been found statistically significant; however, in this case (Figure 9), the discolouration was less conspicuous. In addition, the colour-change effect of the laser power was partially over masked by the oak wood structure heterogeneity underlying the wood colour variability, as illustrated in Figure 10. Table 3. Basic statistical characteristics of colour coordinates for oak wood surfaces engraved with CO 2 laser at different values of raster density and laser power (x represents the mean and s is the standard deviation). The coordinate a* values were increasing with increasing raster density until 5 mm −1 ; after this raster value, a decrease by 25 to 33 % was recorded compared to the referential specimens (corresponding to the laser power). With increasing laser power and raster density, the coordinate a* exhibited a progressive loss in red saturation, because of a moderate shift towards an achromatic colour. The coordinate b* values decreased proportionally, over the entire raster range, at all the laser power values. The final decrease in this coordinate was from 41 to 63%. Figure 9. Raster density-dependent values of colour coordinates L*, a*, b* and total colour difference ΔE* at varying CO2 laser power.

Colour
As for the discussed engraving parameters, the changes in colour coordinates were primarily a consequence of the varying the raster density. The effect of laser power has also been found statistically significant; however, in this case (Figure 9), the discolouration was less conspicuous. In addition, the colour-change effect of the laser power was partially over masked by the oak wood structure heterogeneity underlying the wood colour variability, as illustrated in Figure 10. The impacts of the changes in the individual colour coordinates on the total colour change ΔE* are illustrated in Figure 9. Just at the minimum laser power (8%) and the minimum raster density (2 mm −1 ), the ΔE values calculated using Equation (7) ranged from 9 to 12, corresponding, according to the six-degree scale proposed by [33], to the degree five, expressing an easily visible change in relation to the original colour. With grading the laser power and raster density, the total colour difference values were far beyond 12, corresponding, in accordance with the reported scale, to the total new colour compared to the original. Qualitatively, these changes were like the ones reported by [10,11,16,20,24]. Quantitatively, our results were different, as we studied discolouration on different wood species and under varying laser parameters (head movement rate, focal distance and similar). This fact was an incentive for exploring the dependence of colour change on the total irradiation dose. The result was finding that all the radiation parameters (laser power, head movement rate, focal distance, and raster density) significantly influencing irradiation energy and wood colour could be substituted with a single variable H´, representing the total irradiation dose. The detected dependences of colour coordinates and total colour difference on irradiation dose are displayed in Figure 11. The results confirmed a strong dependence on the discussed colour coordinates on the total irradiation dose. The impacts of the changes in the individual colour coordinates on the total colour change ∆E* are illustrated in Figure 9. Just at the minimum laser power (8%) and the minimum raster density (2 mm −1 ), the ∆E values calculated using Equation (7) ranged from 9 to 12, corresponding, according to the six-degree scale proposed by [33], to the degree five, expressing an easily visible change in relation to the original colour. With grading the laser power and raster density, the total colour difference values were far beyond 12, corresponding, in accordance with the reported scale, to the total new colour compared to the original. Qualitatively, these changes were like the ones reported by [10,11,16,20,24]. Quantitatively, our results were different, as we studied discolouration on different wood species and under varying laser parameters (head movement rate, focal distance and similar). This fact was an incentive for exploring the dependence of colour change on the total irradiation dose. The result was finding that all the radiation parameters (laser power, head movement rate, focal distance, and raster density) significantly influencing irradiation energy and wood colour could be substituted with a single variable H , representing the total irradiation dose. The detected dependences of colour coordinates and total colour difference on irradiation dose are displayed in Figure 11. The results confirmed a strong dependence on the discussed colour coordinates on the total irradiation dose.
The quantification of the energy amount supplied onto the oak wood surface with a laser beam and the identification of irradiation-dose influence on discolouration provide conditions for targeted change in colour space in the relevant wood species treated with a CO 2 laser. The detected dependence of colour changes on the total irradiation dose could enable engraving plastic patterns with a wider range of colour hues on the laser-treated surfaces for particular wood species, as has been discussed in [14]. Figure 11. Colour coordinates and total colour difference depending on the irradiation dose generated by CO2 laser.
The quantification of the energy amount supplied onto the oak wood surface with a laser beam and the identification of irradiation-dose influence on discolouration provide conditions for targeted change in colour space in the relevant wood species treated with a CO2 laser. The detected dependence of colour changes on the total irradiation dose could enable engraving plastic patterns with a wider range of colour hues on the laser-treated surfaces for particular wood species, as has been discussed in [14].

Wetting and Surface Free Energy in Engraved Wood Surfaces
Compared to the original oak wood surfaces, the engraved ones demonstrated much better wettability with both testing liquids. The results for oak wood wetted with water and diiodomethane are summarised in Table 3. Figure 12 illustrates the variability of the contact angles during the wetting process at a raster density of 20 mm −1 , at all three laser power values. Each curve in the graph represents one particular characteristic measurement corresponding to the maximum contact angle at the moment of drop application onto the surface for the relevant treatment way. In many cases, the observed wetting was perfect, equally for water and for diiodomethane. This fact was also reflected in significantly higher values of surface free energy (Figure 13), ensuing from increased disperse and polar components (Table 4). These results are opposite to the results reported in [22], the authors of which observed enhanced hydrophobicity on pine wood surfaces treated with CO2 laser. The authors report pine wood surface temperature reaching up to 200 °C, without the surface carbonisation. The work [59] reports a significant increase in hydrophobicity of a heat-treated beech wood surface at 200 °C, which was reflected in contact angle values bigger than 90°. Worsened wetting of heat-treated modified beech wood was also observed by other authors [60,61]. Figure 11. Colour coordinates and total colour difference depending on the irradiation dose generated by CO 2 laser.

Wetting and Surface Free Energy in Engraved Wood Surfaces
Compared to the original oak wood surfaces, the engraved ones demonstrated much better wettability with both testing liquids. The results for oak wood wetted with water and diiodomethane are summarised in Table 3. Figure 12 illustrates the variability of the contact angles during the wetting process at a raster density of 20 mm −1 , at all three laser power values. Each curve in the graph represents one particular characteristic measurement corresponding to the maximum contact angle at the moment of drop application onto the surface for the relevant treatment way. In many cases, the observed wetting was perfect, equally for water and for diiodomethane. This fact was also reflected in significantly higher values of surface free energy (Figure 13), ensuing from increased disperse and polar components (Table 4). These results are opposite to the results reported in [22], the authors of which observed enhanced hydrophobicity on pine wood surfaces treated with CO 2 laser. The authors report pine wood surface temperature reaching up to 200 • C, without the surface carbonisation. The work [59] reports a significant increase in hydrophobicity of a heat-treated beech wood surface at 200 • C, which was reflected in contact angle values bigger than 90 • . Worsened wetting of heat-treated modified beech wood was also observed by other authors [60,61]. The wetting values also differed from the results obtained for engraved beech wood surfaces by [11]. These different oak wood wetting values were caused by higher irradiation doses per unit area inducing more advanced destruction of lignin and cellulose, more frequent occurrence of carbonised layer, more pronounced roughness, and, consequently, enhanced porosity on the oak wood surface. On such surface, the liquid applied immedi-  The wetting values also differed from the results obtained for engraved beech wood surfaces by [11]. These different oak wood wetting values were caused by higher irradiation doses per unit area inducing more advanced destruction of lignin and cellulose, more frequent occurrence of carbonised layer, more pronounced roughness, and, consequently, enhanced porosity on the oak wood surface. On such surface, the liquid applied immediately soaked into the substrate.
The wetting values also differed from the results obtained for engraved beech wood surfaces by [11]. These different oak wood wetting values were caused by higher irradiation doses per unit area inducing more advanced destruction of lignin and cellulose, more frequent occurrence of carbonised layer, more pronounced roughness, and, consequently, enhanced porosity on the oak wood surface. On such surface, the liquid applied immediately soaked into the substrate. Figure 13. Surface free energy of oak wood surface, dependent on laser power and raster density.
The results of wetting and of surface free energy in oak wood surfaces treated with laser beam allow us to expect an appropriate adhesion of film-forming materials to wood. The pronounced roughness, however, may play a negative role in gluing such surfaces. Moreover, in the case of surfaces engraved at higher raster densities, it is necessary to consider an unstable carbonised layer ( Figure 5), with weak adhesion, possible to peel off from the substrate easily.

Conclusions
The experimental results demonstrate that the studied laser engraving parameters (laser power and raster density) affected the wood mass loss, wood surface chemical structure, and morphology significantly, as well as the surface properties (colour and hydrophilicity) of the analysed oak wood.
The energy amount supplied onto the wood surface under varying irradiation conditions was expressed through the total irradiation dose H'. This variable has been confirmed as an appropriate one for substituting all the irradiation parameters altogether (laser power, laser head movement rate, focal distance, and raster density).
During the engraving, thin surface layers manifested immediate destruction and sublimation. The wood mass loss increased linearly with the irradiation dose. The wood mass destruction was not uniform over the whole surface. In the early wood with dominant early vessels, the trenches created by the laser beam were deeper than in the late wood, and this was reflected in the more pronounced roughness of the early wood.
The major impact on the roughness variability was detected for the raster density. With increasing raster density over the whole studied range the roughness parameters Ra, Figure 13. Surface free energy of oak wood surface, dependent on laser power and raster density.
The results of wetting and of surface free energy in oak wood surfaces treated with laser beam allow us to expect an appropriate adhesion of film-forming materials to wood. The pronounced roughness, however, may play a negative role in gluing such surfaces. Moreover, in the case of surfaces engraved at higher raster densities, it is necessary to consider an unstable carbonised layer ( Figure 5), with weak adhesion, possible to peel off from the substrate easily.

Conclusions
The experimental results demonstrate that the studied laser engraving parameters (laser power and raster density) affected the wood mass loss, wood surface chemical structure, and morphology significantly, as well as the surface properties (colour and hydrophilicity) of the analysed oak wood.
The energy amount supplied onto the wood surface under varying irradiation conditions was expressed through the total irradiation dose H'. This variable has been confirmed as an appropriate one for substituting all the irradiation parameters altogether (laser power, laser head movement rate, focal distance, and raster density).
During the engraving, thin surface layers manifested immediate destruction and sublimation. The wood mass loss increased linearly with the irradiation dose. The wood mass destruction was not uniform over the whole surface. In the early wood with dominant early vessels, the trenches created by the laser beam were deeper than in the late wood, and this was reflected in the more pronounced roughness of the early wood.
The major impact on the roughness variability was detected for the raster density. With increasing raster density over the whole studied range the roughness parameters Ra, Rz, and Rt increased parallel with and similarly perpendicular to the grain. The more pronounced roughness as well as more distinct variability were recorded perpendicular to the grain course. The roughness parameter values Ra and Rz increased linearly with increasing irradiation dose.
The wood surface treatment with the laser induced wood discolouration. Already at the minimum laser power and the minimum raster density, evident discolouration could be declared. The next laser parameter rises caused a totally novel wood surface colour compared with the original one.
The discolouration was primarily due to heat-induced processes associated with the degradation of carbonyl groups in lignin and hemicelluloses in which the broken bonds C=O caused changes in the contents of chromophores determining the colour.
A close dependence of discolouration on the irradiation dose was identified, regardless of the values of the individual irradiation parameters. This dependence can be declared as very important, because it serves a means for targeted discolouration of wood engraved with a CO 2 laser.
The engraved surfaces exhibited significantly enhanced wettability with the standard polar and non-polar liquids, which was reflected in the higher surface free energy. The results concerning the wetting and the surface free energy values obtained for the laserengraved oak wood surfaces allow us to suppose an appropriate spreading of film-forming material on wood surface, and, correspondingly, appropriate adhesion of film-forming materials to wood.
Author Contributions: J.K. conceiving and designing the experiments; I.K. measuring and analysis of the FTIR data; M.A. measuring and analysis of the laser power data; J.K. and M.A. measuring and analysis of the colour, roughness, contact angle, and calculating the surface free energy; J.K., I.K. and M.A. wrote the paper. All authors have read and agreed to the published version of the manuscript.