Selective detection of naphthalene with nanostructured WO 3 gas sensors prepared by pulsed laser deposition

Pulsed laser deposition (PLD) at room temperature with a nanosecond laser was used to prepare WO3 layers on both MEMS microheater platforms and Si/SiO2 substrates. Structural characterization showed that the layers are formed of nanoparticles and nanoparticle agglomerates. Two types of layers were prepared, one at an oxygen partial pressure of 0.08 mbar and one at 0.2 mbar. The layer structure and the related gas sensing properties were shown to be highly dependent on this deposition parameter. At an oxygen pressure of 0.2 mbar, formation of ε-phase WO3 was found, which is possibly contributing to the observed increase in sensitivity of the sensor material. The gas sensing performance of the two sensor layers prepared via PLD was tested for detection of volatile organic compounds (benzene, formaldehyde and naphthalene) at ppb level concentrations, with various ethanol backgrounds (0.5 and 2 ppm) and gas humidities (30, 50 and 70 % RH). The gas sensors were operated in temperature cycled operation. For signal processing, linear discriminant analysis was performed using features extracted from the conductance signals during temperature variations as input data. Both WO3 sensor layers showed high sensitivity and selectivity to naphthalene compared to the other target gases. Of the two layers, the one prepared at higher oxygen partial pressure showed higher sensitivity and stability resulting in better discrimination of the gases and of different naphthalene concentrations. Naphthalene at concentrations down to 1 ppb could be detected with high reliability, even in an ethanol background of up to 2 ppm. The sensors show only low response to ethanol, which can be compensated reliably during the signal processing. Quantification of ppb level naphthalene concentrations was also possible with a high success rate of more than 99 % as shown by leave-one-out cross validation.


Introduction
In order to evaluate and assess indoor air quality (IAQ), different types of gaseous chemical compounds have to be considered.In addition to carbon monoxide (CO), carbon dioxide (CO 2 ) and nitrogen dioxide (NO 2 ), low concentrations of volatile organic compounds (VOCs) play a significant role in deteriorating the quality of breathing air in buildings (World Health Organization, 2010;Bernstein et al., 2008).Exposure to these substances, even at low concentrations, can lead to severe negative effects on human health.For VOCs, health problems mainly include damage to the respiratory system and skin irritations (Jones, 1999).Additionally, some VOCs are proven to be carcinogenic (e.g., benzene, World Health Organization, 2010) or are suspected to be carcinogenic (e.g., formaldehyde, Guo et al., 2004).Based on toxicity and prevalence, according to the World Health Organization (WHO) and the INDEX project (Koistinen et al., 2008), the highest priority VOCs for IAQ are formaldehyde, benzene and naphthalene.For naphthalene, the WHO guidelines suggest values below 0.01 mg m −3 corresponding to 1.9 ppb (World Health Organization, 2010).The main health concerns for naphthalene are lesions in the respiratory tract, including tumors in the upper airways (World Health Organization, 2010).
In order to detect such small concentrations of VOCs without the need for expensive and time-consuming analytical measurements (e.g., GC-MS analysis, Wu et al., 2004), metal oxide semiconductor gas sensors can be applied.Detection of VOCs in the ppb range with such sensors has been successfully demonstrated using temperature cycled operation (TCO) and pattern recognition signal processing with ceramic-based thick film sensors (Leidinger et al., 2014); however, significant cross-sensitivity to ethanol was found mainly for SnO 2 -based sensors.In order to reduce this crosssensitivity, WO 3 layers were investigated.To obtain highly sensitive sensors with small thermal time constants, the MOS thin film layers were produced on microheater substrates by pulsed laser deposition (PLD).We found that these sensors show a high response to naphthalene in the relevant concentration range with high selectivity compared to other gases, especially relevant interferent gases for indoor air quality assessment, e.g., ethanol.
Pulsed laser deposition is a method for depositing a variety of materials ranging from epitaxial thin films (Hussain et al., 2005) to highly porous nanostructured layers (Balandeh et al., 2015).Porous nanostructured layers have been studied especially in the context of gas sensing materials (Caricato et al., 2009;Nam et al., 2006).PLD offers many advantages compared to other deposition methods, for example easily controllable film composition by deposition parameters, and a good repetition of stoichiometry of the target material in the films deposited on the substrate.When using nanosecond laser PLD, as in this study, with a high oxygen partial pressure in the deposition chamber, nanoparticle formation starts during the deposition process leading to a highly porous nanostructured layer (Harilal et al., 2003;Infortuna et al., 2008;Huotari et al., 2015).These types of layers are very suitable for gas sensing purposes because of their high specific surface area.
WO 3 as a material has been widely studied as it offers a large range of possibilities in practical applications, e.g., in gas sensing (Kohl et al., 2000;Wang et al., 2008;Balazsi et al., 2008), and photocatalytic water splitting (Pihosh et al., 2015).There are several methods to produce WO 3 layers ranging from thick film and thin films technologies to chemical methods (Zheng et al., 2011).In this study, PLD was utilized for depositing WO 3 layers on MEMS microhotplates to produce low-cost and high-performance gas sensor devices.
The performance of the PLD sensor layers has been evaluated in test gas measurements.The three high-priority VOCs, benzene, formaldehyde and naphthalene, have been applied in concentrations below, at and above the respective guideline values, and ethanol has been added as an interferent gas in much higher concentrations in order to simulate typical IAQ applications with background gases from, e.g., cleaning agents or alcoholic beverages.The sensors were operated in dynamic operation using temperature cycled operation (TCO), which is a well-known method for increasing sensitivity and selectivity of gas sensor systems based on MOS sensors (Heilig et al., 1997;Lee and Ready, 1999;Paczkowski et al., 2013;Baur et al., 2014).The resulting signals, after pre-processing and feature extraction, were analyzed by linear discriminant analysis (LDA), a multivariate pattern recognition method which separates different classes of input data while trying to group data sets of the same assigned group (Klecka, 1980;Gutierrez-Osuna, 2002).The combination of TCO and LDA has shown to improve selectivity and sensitivity of gas sensors, both MOS sensors (Gramm et al., 2003;Meier et al., 2007;Reimann and Schütze, 2012;Leidinger et al., 2014) and other types, e.g., GasFET devices (Bur et al., 2012).

Sensor layer deposition
A XeCl laser with a wavelength of λ = 308 nm was used to produce WO 3 layers on both Si/SiO 2 substrates and commercial microheater MEMS platforms from a ceramic WO 3 pellet.The laser pulse length was 25 ns and pulse fluence was I = 1.25 J cm −2 .In all depositions the substrate temperature was kept at room temperature (RT).Two types of samples have been prepared, one at a low O 2 partial pressures of p(O 2 ) = 0.08 mbar, designation PLD0.08O2, and a second type with a higher partial pressure of p(O 2 ) = 0.2 mbar, designation PLD0.2O2.All samples were annealed in a furnace at 400 • C for 1 h after deposition.The samples deposited on Si/SiO 2 substrates were used as reference samples in structural characterization of the layers, and samples with the MEMS heaters were used in gas sensing measurements.
A Bruker D8 Discover device was used in X-ray diffraction studies, and Raman spectroscopy studies were performed with a HORIBA Jobin Yvon LabRAM HR800 in order to study the crystal structure and symmetry of the layers.The surface morphology and the film composition of the samples were studied with a Veeco Dimension 3100 atomic force microscope (AFM) and with Zeiss Sigma FESEM device.

Crystal structure characterization of the sensing layers
The grazing incidence diffraction (GID) method of the Xray diffraction was used to characterize the WO 3 layers annealed at 400 • C for 1 h.The results are shown in Fig. 1.A clear difference in the crystal structure can be seen.The phase composition of layers deposited at p(O 2 ) = 0.08 mbar is mostly of the monoclinic γ phase of WO 3 , but in the samples deposited at p(O 2 ) = 0.2 mbar, also the ferroelectric monoclinic ε phase of WO 3 is present.This is emphasized especially by the ( 110) and (−112) reflections located at 2θ ≈ 24.0 and 33.3 • , respectively (Johansson, 2012).However, one must remember that the crystal structures of the γ phase and the ε phase are quite similar, and thus both of them have XRD reflections either in the same 2θ angels or very close to each other.The average grain size of both types of samples after annealing was determined to be around 30 nm using the Warren-Averbach method for XRD data (Marinkovic et al., 2001).In Fig. 2, the Raman spectroscopy studies performed to the WO 3 layers are presented.The Raman spectra of the as-deposited samples immediately after deposition without any heat treatment are shown in Fig. 2a, and in Fig. 2b the spectra of the layers after annealing at 400 • C for 1 h are presented.An interesting property of the deposition process can be identified in the non-annealed samples.When the O 2 partial pressure is 0.08 mbar, the samples seem to be in an amorphous state after deposition, but when the O 2 pressure is 0.2 mbar, some crystallization is already evident during the deposition process at RT, even before any heat treatment to the layers.However, from Fig. 2b, showing the Raman spectra of the samples after the annealing process, it is clearly seen that after heat treatment in a furnace both films have a more crystallized structure.It is also again evident that the layers PLD0.08O2 are composed mostly of γ phase, but the samples PLD0.2O2 have also the ε phase in their crystal structure, verified from the Raman modes at wavenumbers 67,97,144,183,203,272,303,370,425,644, and 680 cm −1 (Wang et al., 2008;Johansson et al., 2012;Souza Filho et al., 2000).Similarly as the reflections in XRD measurements, both γ phase and ε phase have Raman modes either at same wavenumbers or very close to each other.
At this point it should be noted that usually the ferroelectric ε phase only exists in temperatures below −40 • C.However, different studies (Wang et al., 2008;Righettoni et al., 2010;Johansson et al., 2012) show that the ε phase can exist in a solid-state form also at temperatures above RT.The reason for this is believed to be the small particle size of the samples, similarly as in the samples presented in this study.Also, the existence of the ε phase in the WO 3 composition has been proven to enhance WO 3 structures sensitivity to acetone (Wang et al., 2008;Righettoni et al., 2010;Sood and Gouma, 2013).The reason was suggested to be the ferroelectricity of the ε phase, namely the spontaneous electric dipole moments it possesses, which are then highly contributing to the chemical reaction between the WO 3 surface and the target gas.

Film composition characterization of the sensing layers
The surface morphology of the as-deposited and annealed WO 3 samples was studied by atomic force microscopy, and the results are shown in Fig. 3.The surface micrographs of the as-deposited layer and the annealed layer of sample PLD0.08O2 are shown in Fig. 3a and b, respectively.The sample surfaces consist of small agglomerates of nanopar- ticles.In Fig. 3c and d, the surface micrographs of the asdeposited layer and the annealed layer of sample PLD0.2O2 are shown, respectively.The layers consist also of small nanoparticles, but agglomerated to bigger clusters.Also, the layer structure is much rougher and more porous than on the samples deposited at 0.08 mbar O 2 .It can also be clearly seen that the annealing process at 400 • C for 1 h does not have a great effect on the surfaces of the samples.In both cases, the average surface roughness value R q was the same before and after the annealing process, being R q = 5.5 nm for sample PLD0.08O2, and R q = 42.2 nm for sample PLD0.2O2.The crystallization, which is observed for the PLD0.08O2sample seems to be a local process that does not involve larger-scale material transport.
Scanning electron microscopy was used to further study the film composition of the samples.Both surface micrographs and cross-section micrographs were taken from the samples.In Fig. 4a and b the cross-section micrograph and surface micrograph of the as-deposited sample fabricated at p(O 2 ) = 0.08 mbar are shown, respectively.The crosssection graph shows that the film is composed of small nanoparticle agglomerates formed as pillar-like structures, with some porosity in between the columns.The surface graph shows that the film surface is formed of small nanoparticle agglomerates and thus verifies the measurements made with AFM.The cross-section micrograph and surface micrograph of the as-deposited sample fabricated at p(O 2 ) = 0.2 mbar are shown in Fig. 4b and c, respectively.Now the film composition is much more porous and rough compared to the PLD0.08O2film, and also the nanoparticle agglomerate size is larger.The agglomerates form clearer pillar-like morphology on top of the substrate.The surface of the film is highly porous, concurrent with the AFM measurements performed.
3 Gas sensor performance

Gas test measurement setup
The gas sensing performance of the two sensor layers was evaluated in an extensive test measurement.The three target VOC gases were applied in three concentrations each, the middle concentrations representing the respective WHO guideline values of 0.1 mg m −3 (81 ppb) for formaldehyde and 0.01 mg m −3 (1.9 ppb) for naphthalene (World Health Organization, 2010), as well as the European Union guideline value of 5 µg m −3 (1.6 ppb) for benzene (European Parliament, Council of the European Union, 2008).Additionally, ethanol was introduced as a background gas in two concentrations, both much higher than the target gas concentrations.As the third varied parameter, the gas humidity was set in three steps.Table 1 shows all gases and concentrations.
The test gases were generated in a gas mixing system specifically designed for trace gases by Helwig et al. (2014).The gases were mixed into zero air produced by two cascaded zero air generators.Ethanol and formaldehyde test gases were taken from gas cylinders and diluted into the zero air carrier gas stream, either with a one-step dilution (ethanol) or a two-step dilution (formaldehyde).Benzene and naphthalene test gases were generated from permeation tubes in permeation ovens.Each target VOC concentration was set twice during each combination of humidity and ethanol background, first from highest to lowest concentration, then back to the highest concentration.The length of each VOC run was 30 min; between two trace VOC applications the sensors were purged with zero air with the respective ethanol and humidity configuration.In total, 90 different gas mixtures were generated and tested with the sensors; the total length of the measurement was 123 h.
For temperature cycled operation of the sensors, a rampup-ramp-down approach was chosen; the temperature of the microheaters was increased from 200 to 400 • C in 20 s and Table 1.Test gas setup.Each gas concentration was applied at each EtOH background and humidity level for 30 min; between gas exposures sensors were exposed to background for 30 min.then reduced back to 200 • C in the same time, creating a 40 s cycle (see Fig. 5).The two sensor types clearly show differing behavior during the temperature cycle, especially during increase of the heater temperature.

Gas measurement results
For a first signal evaluation, quasi-static sensor signals were extracted from the raw sensor signal data sets.These signals were generated by plotting the signal value of certain points in the TCO cycle for each cycle, i.e., over the course of the complete measurement.An example of the PLD0.2O2sensor is given in Fig. 6.Four points of the temperature cycle were selected, indicated in Fig. 5.A section of the measurement was chosen in which all three test VOCs are applied at 30 % RH gas humidity and without ethanol background.
The sensor response to all concentrations of naphthalene is clearly visible, as well as the much lower responses to the other VOCs.By normalizing the signals, i.e., calculating the relative change of conductance G/G 0 , the sensor response to naphthalene at the different points in the cycle can be determined (cf.Fig. 7).In this plot, it can be seen that the highest sensitivity, of the chosen points in the cycle, is during cooling of the sensor, 30 s after start of the cycle.The  1) in air at 30 % RH without EtOH background.For later normalization, the value G 0 was extracted during background before the first exposure to naphthalene as indicated.change of conductance is approx.15 % for 1 ppb of naphthalene.The reason for the high selectivity to naphthalene is not known in detail and should be studied more closely.However, it was reported earlier that WO 3 has a specific response to aromatic compounds (Sauerwald, 2008) and that in general higher molecular weight leads to an increased sensitivity for this material (Sauerwald, 2008;Kohl et al., 2000).
The dynamic sensor signal patterns were analyzed using LDA.The whole procedure of generating data sets for LDA from dynamic gas sensor signals and the options for LDA application were described by Bur et al. (2014).As input data for each sensor, a limited number of features were extracted from the respective sensor signal of each temperature cycle.In this case, the cycle was divided into 20 segments of equal length (2 s).For each segment, the mean value and the slope of the sensor signal were calculated and used as features.This generates a feature vector with 40 values for each temperature cycle.These data sets were grouped into different groups, depending on the desired data analysis.
For the first analysis, the complete data set, containing 2799 feature vectors, was used as input data for the LDA.All cycles which contain a certain VOC were grouped together, regardless of VOC concentration, humidity and ethanol background.This results in one group for each target VOC and a "background" group, which contains the data sets of the TCO cycles which ran when no trace VOC was applied.This analysis checks the performance of the sensors to discriminate the target gases in varying humidity and background conditions.The LDA result plots for the two sensors with the PLD layers are shown in Fig. 8, left for the PLD0.08O2sensor and right for the PLD0.2O2type.For both sensors, the background, formaldehyde and benzene groups are overlapping strongly, while the naphthalene group is more separated.Especially for the PLD0.2O2sensor the naphthalene group is nearly completely split from the other gases.
As a quantitative measure of the discrimination result, leave-one-out cross validation was performed on the data (LOOCV; Gutierrez-Osuna, 2002), with k nearest neighbors (kNN, k = 5) as classifier.This method calculates how many  Of the two sensors, the one with the PLD0.2O2layer shows the better discrimination result overall, with 71.9 % correctly classified data points.The PLD0.08O2 achieves 66.5 %.The same order is given for classification of the naphthalene data points.The layer deposited at 0.2 mbar of O 2 shows a much better classification performance (99.2 % correct classifications) compared to the second sensor.The LOOCV results are summarized in Table 2, LDA 1a and 1b.
In the next step, it was checked if a quantification of the naphthalene concentration was possible.The full data set was used again.Each naphthalene concentration was assigned a group (with all humidities and ethanol backgrounds), and the other VOCs were assigned to the background group.The results for both tested sensors are plotted in Fig. 9.In this analysis, the PLD0.2O2sensor shows the best separation of the groups again.There are some data points from the 1 ppb group located in the 2 ppb group.Otherwise, the three naphthalene concentrations are well lined up along the first discriminant axis.The PLD0.08O2 sensor shows less discrimination of the groups, which is also shown in LOOCV results (see Table 2, LDA 2).
The same analysis was attempted with a reduced data set, which included only the segments of the measurement without the ethanol background (950 sensor cycles in to-tal).This significantly increases the quality of discrimination and thus naphthalene quantification for the sensors (see Fig. 10).Especially the PLD0.2O2sensor has excellent separation of the naphthalene concentrations along the first discriminant function; the PLD0.08O2sensor layer has much wider groups.The LOOCV results, listed in Table 2, LDA 3, show nearly 100 % correct classifications for the better sensor (PLD0.2O2); the second sensor layer also has over 99 % success rate.
Another evaluation of the sensor data was performed regarding quantification of gas humidity.The full data set was split into three groups, for the three gas humidities set during the measurement.This LDA run checks the sensors' crosssensitivity to humidity, which also shows if the sensor would be able to measure the gas humidity.See Fig. 11 for the LDA result plots.These plots and the corresponding LOOCV results (Table 2, LDA 4) show very good discrimination of the humidities, which means that the sensor layers deposited by PLD have considerable sensitivity to water.However, this also means that the sensors could be used to measure the gas humidity or that the sensor performance could be monitored by comparing the predicted gas humidity with the value of a reference humidity sensor.
If the gas humidity is known, either from the gas sensor itself or from an additional humidity sensor, the quality of the signal processing can be improved by calculating different  LDA projections for several humidity ranges.A simple version of this approach has been tested by reducing the data set to the sensor signals acquired in one humidity, in this case all cycles measured in 50 % RH.The LDA result for gas identification for this is shown in Fig. 12.Compared to the full data set with three humidity levels (Fig. 8), the separation of the gases is clearly improved, especially for the PLD0.08O2sensor.LOOCV also shows significant improvement of the classification results (Table 1, LDA 5).For all gases, the ratio of correct classifications was raised from 66.5 to 85.4 % for sensor PLD0.08O2 and from 71.9 to 86.7 % for sensor PLD0.2O2.For the naphthalene group, perfect classification was achieved for the sensor with 0.2 mbar O 2 partial pressure (Table 1, LDA 5b).A hierarchical data processing approach, which in the first step determines the humidity and in the second step classifies the gas, seems promising.

Conclusions
Nanoporous WO 3 gas sensing layers have successfully been prepared via nanosecond pulsed laser deposition, characterized, and tested for their gas sensing performance, especially for use in IAQ applications.
The characterization of the PLD sensing layers showed films formed of nanoparticle agglomerates with pillar-like morphology.The films were highly porous in their structure, especially when higher oxygen partial pressure was used during the deposition process.The crystal structure of the films was also dependent on the oxygen partial pressure with higher O 2 pressure resulting in the formation of WO 3 ε phase, which in bulk WO 3 samples is only stable at temperatures below −40 • C. The ε phase in the PLD films was observed to withstand annealing at 400 • C for 1 h probably due to its monocrystalline structure.
The two compared PLD sensor layer samples showed significant differences in their gas sensing performance.The layer deposited at higher oxygen partial pressure displayed improved response and excellent selectivity to naphthalene.
Both discrimination of different target VOCs and quantification of naphthalene were more successful with this sensor.
Naphthalene concentrations down to 1 ppb could be quantified with the sensor layer deposited at 0.2 mbar of O 2 with nearly 100 % success rate as determined by leave-one-out cross validation when no ethanol background was present.Even in varying ethanol background of up to 2 ppm, quantification was still successful for 94 % of all temperature cycles for this sensor.Detection of the presence of naphthalene in concentrations of 1 ppb or more had a success rate of more than 99 %, again determined by LOOCV.
The highly porous structure and possibly the formation of the WO 3 ε phase resulting from the higher oxygen pressure during PLD thus boost the gas sensing performance considerably.With the selected parameters, suitable gas sensing layers for detection and quantification of naphthalene have been obtained.Future investigations will address further improvement of the gas sensitive layers, e.g., by further variations of the deposition parameters as well as introduction of additional nanoparticles for doping and catalytic activation.We are also planning to study PLD deposition based on picosecond laser pulses which allows further optimization of the deposition parameters.

Figure 1 .
Figure 1.X-ray diffraction spectra of the deposited WO 3 layers after annealing in air (400 • C for 1 h).

Figure 3 .
Figure 3. Atomic force microscopy surface micrographs of the deposited WO 3 layers showing the influence of oxygen partial pressure during deposition (a, b) as-deposited samples, (c, d) after annealing at 400 • C for 1 h.

Figure 5 .
Figure 5. Gas sensor temperature cycle (green) and corresponding sensor signals of the two PLD WO 3 sensors in air at 30 % RH.The dashed lines indicate the selected points for quasi-static sensor signal analysis (cf.Fig. 6).

Figure 6 .
Figure 6.Quasi-static sensor signals for the four points indicated in Fig.5during exposure of the PLD0.2O2sensor to formaldehyde, naphthalene and benzene at relevant ppb levels (Table1) in air at 30 % RH without EtOH background.For later normalization, the value G 0 was extracted during background before the first exposure to naphthalene as indicated.

Figure 7 .
Figure 7. Normalized quasi-static signals for sensor PLD0.2O2 for the four selected points in the temperature cycle (cf.Fig. 5).