Lauren
Dalecky
a,
Francesco
Sottile
b,
Linda
Hung
c,
Laure
Cazals
a,
Agnès
Desolneux
d,
Aurélia
Chevalier
e,
Jean-Pascal
Rueff
fg and
Loïc
Bertrand
*a
aUniversité Paris-Saclay, ENS Paris-Saclay, CNRS, Photophysique et Photochimie Supramoléculaires et Macromoléculaires, 91190 Gif-sur-Yvette, France. E-mail: loic.bertrand@ens-paris-saclay.fr
bETSF and LSI, CNRS, CEA/DRF/IRAMIS, École Polytechnique, Institut Polytechnique de Paris, F-91120 Palaiseau, France
cEnergy and Materials Division, Toyota Research Institute, Los Altos, CA 94022, USA
dUniversité Paris-Saclay, ENS Paris-Saclay, CNRS, Centre Borelli, 91190 Gif-sur-Yvette, France
eConservation of Cultural Heritage – Aurélia Chevalier Sàrl, Route des Jeunes 4bis, 1227 Les Acacias, Genève, Switzerland
fSynchrotron SOLEIL, L'Orme des Merisiers, BP 48 Saint-Aubin, 91192 Gif-sur-Yvette, France
gLaboratoire de Chimie Physique – Matière et Rayonnement (LCPMR), Sorbonne Université, CNRS, 75005 Paris, France
First published on 26th December 2023
Inelastic X-ray scattering (IXS) spectroscopy has been used in many fields of solid-state physics and theoretical chemistry as an accurate and quantitative probe of elementary excitations. We show that non-resonant IXS spectra in the energy loss range below 100 eV exhibit a strong contrast across a wide range of commercially available pigments, opening new routes for their discrimination. These signatures combine plasmonic transitions, collective excitations and low energy absorption edges. We have performed IXS to discriminate different artists' pigments within complex mixtures and to quantitatively determine rutile and anatase polymorphs of TiO2. The integration of experimental data on pigment powders with suitable ab initio simulations shows a precise fit of the spectroscopic data both in the position of the resonances and in their relative intensity.
X-ray absorption spectroscopy (XAS) methods, often combined with other X-ray techniques, have become key in gaining electronic and geometric structure information (oxidation state, coordination number, and site symmetry) for inorganic paint materials in both crystalline and amorphous forms.10,11 For example, X-ray absorption near-edge structure (XANES) spectroscopy combined with X-ray fluorescence (XRF) has been used to identify oxidised arsenic species resulting from the migration of arsenic sulphide-based yellows (As2S3), reds and oranges (As4S4) in photodegraded 19th century artworks.12 The combination of μ-XANES and μ-XRF with micro-X-ray diffraction (XRD) has been used to evaluate the darkening of chromium-based yellows (PbCrO4, PbCr1−xSxO4) through photoreduction in one of the paintings in Van Gogh's Sunflowers series (1888–1889) by determining the oxidation state and spatial distribution of the paint compounds on the microscale.7 These methods were also coupled to understand the reactivity of crystal defects in cadmium-based (CdS) nanoparticles in yellows from Pablo Picasso's Femme (1907) by correlating elemental and structural distributions with optical emissions to identify degradation mechanisms.13 In these studies, analysis was performed on prepared samples from the painting, usually in the form of cross-sections. Alternatively, some studies only explored the sample surface due to the limited depth of information constrained by the appropriate edge energy.
XAS has several limitations for the study of artistic materials: (i) The depth probed is imposed by the energy of the studied edge, which requires the study of cross-sections or sample surfaces that may not be representative of the bulk material or whose preparation may alter heterogeneous or reactive samples.10 (ii) XAS on artistic materials is mainly performed in the hard or tender X-ray regime, which is not suitable for studying elements with a low atomic number (Z).10,14 This is particularly problematic for recent paint formulations, which may include a range of synthetic organic dyes or extenders made of light elements.15 (iii) XAS yields only a partial view of the electronic structure as it primarily reflects the projected density of empty states of the excited atom, accessible mainly through the dipole selection rule. It therefore lacks information about the symmetry-forbidden density of empty states, the density of occupied states, and, importantly, about elementary excitations that are beyond the X-ray range but might be valuable for material characterisation.
Here, we apply non-resonant inelastic X-ray scattering (IXS) as a powerful method for the identification and analysis of mineral pigments. IXS is a technique capable of revealing electronic excitations over a wide range of energies and momenta. It not only addresses different degrees of freedom compared to infrared or Raman spectroscopy (successfully applied to titanium dioxide16,17), but it is also complementary to optical techniques, like absorption, optical reflectivity, or photoluminescence.18 Specifically, IXS, as a scattering technique, allows analysis of a new dimension, namely, the momentum transfer (q). This capability unveils excitations that would remain undetected in the dipole limit intrinsic to optical absorption experiments. Another noteworthy emerging technique is the resonant counterpart of IXS, known as RIXS. Similar to IXS, RIXS is a scattering technique that examines both energy and momentum, and it has been proven highly effective in the description of titanium compounds.19–21 RIXS offers a different excitation pathway22 and complements IXS. It is more focused on specific aspects of the excitation, due to its resonant condition of the incoming photon, like electron–phonon coupling,23 d–d excitations, and pre-peak features of X-ray absorption edges. The resonant condition, even though it increases the signal/noise ratio of the experiment, also limits the momentum range that can be investigated. Probably the closest technique to IXS is electron energy loss spectroscopy (EELS), which enables the study of excitations over a wide range of energies with unrivalled spatial resolution. EELS has also been proven effective in highlighting differences in bulk polymorphs of titanium,24 but the spectral differences between rutile, anatase, and brookite are too small to be useful in powder compounds.
For these reasons, we conducted a combined experimental–theoretical study focusing on non-resonant IXS, an approach that allows us to address a broader momentum range and explore various aspects of electronic excitations. In particular, we demonstrate the ability of low energy loss IXS to discriminate between common commercial pigments, by exploiting resonances resulting from the electron structure of the materials, including electron–hole or charge transfer excitations, collective excitations related to the plasmonic response, gap transitions, in addition to absorption edges. We show that significant contrasts within spectroscopic signatures hold promise for material differentiation. The dependence of these features on the momentum transfer provides additional signatures not observed by any other analysis. Finally, we show that titanium dioxide, TiO2, which constitutes the most commonly used white pigment since the 20th century, can be readily identified by IXS with a number of characteristic spectral features. In particular, the two main crystalline phases encountered (rutile and anatase) show spectroscopic differences in very good agreement with theoretical calculations of the electronic dynamic structure factor (DSF) as a function of the momentum transfer. We therefore present a state-of-the-art integration of experimental and theoretical methodologies for studying mineral pigments. At the same time, we wish to provide a practical illustration that highlights the applicability of IXS to a wider spectrum of complex materials than those that have been previously studied using this approach.
The pigment powder samples were mounted by pressing them into the 200 μm diameter holes of a 2.5 cm × 2 cm steel sample holder. The setup for spectroscopy was configured in the forward scattering mode with the spectrometer aligned at a scattering angle 2θ = 55° with respect to the centre of the spectrometer, corresponding to q ranging from 1.9 to 4.0 Å−1 from end to end. Spectra for all samples were recorded in the low-resolution mode (Fig. S2.1a, ESI†) for an energy loss range of ΔE = 0–100 eV with an energy step size of 200 meV for ΔE = 0–30 eV and an energy step size of 500 meV for ΔE = 30–98 eV. The acquisition time was 4 s per energy point (∼20 min per spectrum). High-resolution spectra (Fig. S2.1b, ESI†) were obtained for titanium white pigment samples at the plasmonic resonance ΔE = −2 to 18 eV, with a step size of 200 meV and an acquisition time of 20 s per energy point (∼67 min per spectrum). For spectral imaging, the configuration was changed to the backscattering mode with the spectrometer aligned at a scattering angle of 2θ = 130° for higher spatial resolution, and spectra were acquired in the low-resolution mode (1.2 eV) (Fig. S2.1c, ESI†). Spectra were acquired by scanning horizontally across the sample (y) at a fixed energy for 100 ms per spatial point, for all energies from ΔE = −5 to 25 eV with a step size of 1 eV. The vertical position (z) was modified and the scans were repeated.
No. | Pigment designation | Diagnostic phases |
---|---|---|
1 | Titanium white | TiO2 (rutile) |
2 | Turquoise Clair | CoAl2O4o |
3 | Blanc de Zinc | ZnO |
4 | Noir de Vigne | Fe3O4, Fe2O3s,o |
5 | Noir de Mars | Fe3O4s,o |
6 | Terre Ombre Brûlée | Fe2O3, Fe3O4s,o |
7 | Bleu de Prusse | Fe4[Fe(CN)6]3 |
8 | Noir d'Ivoire | Ca(PO4)3(OH) |
9 | Bleu Outremer | Sodalite-type phasesc,o |
10 | Jaune Indien Imitation | NiC8H6N6O6 |
11 | Rouge de Quinacridone | C22H16N2O2 |
12 | Titanium white | TiO2 (anataseb) |
13 | Titanium white | TiO2 (anataseb) |
14 | Titanium white | TiO2 (rutile) |
15 | Titanium white | TiO2 (rutileb) |
16 | “Rutile” | TiO2 (rutileb) |
We further experimentally examined the plasmon resonance dependence on the momentum transfer q for the two crystal structures of TiO2 found in titanium white pigments, anatase and rutile, focusing on pigment 12 and 1, respectively. For this purpose, we used the specific geometry of the MATRIXS spectrometer at the GALAXIES beamline at the SOLEIL synchrotron facility. The 40 Si(110) analyser crystals are arranged in a compact 8 × 5 array (H × V), divided into 4 panels of 2 × 5 analysers (H × V). Each V column corresponds to a specific q value, whose contribution to the IXS spectrum can be isolated by selecting the corresponding signal on the PSD. Four spectra were plotted, each corresponding to one q value (Fig. 2). The most notable effect, observed for both samples, is the decreasing intensity of the plasmon peak as q increases. Two notable differences between the crystal polymorphs were observed: (1) at lower q (1.9–2.2 Å−1), two broad peaks, centred around 25 eV and 35 eV are present for 1, and observed as one broad peak centred around 30 eV for 12, and (2) a broad peak around 8 eV is highlighted for q at 2.6–2.9 Å−1 for 12 and q at 3.2–3.5 Å−1 for 1.
We then calculated the DSF of rutile and anatase crystals at different q values, both in the XY plane and in the Z direction (both structures are tetragonal). For each q, we then qualitatively compared the experimental results on the powder with the arithmetic mean of the amplitude of the simulated spectra as a function of energy (0–50 eV) in the three directions X, Y and Z. This theory–experiment comparison can only be qualitative since the exact distribution in orientation of TiO2 crystallites in the powder samples is unknown. However, the anisotropy in both structures, although present, is not huge (see Fig. S4.1, ESI†), and the averaging operation permits one to carefully analyse the spectra and recognise the characteristic features. We observed a very good agreement between theoretical and experimental spectra, inspite of the important difference between the experimental sample and the theoretically ideal infinite bulk anatase and rutile. In addition, there is a quantitative juxtaposition between theory and experiments in different spectral features, like the sharp plasmon peak at 14 eV or the broad structure around 50 eV. At all investigated q values, we noticed a correctly matching relative intensity of the different peaks. This permitted us to carry out a full analysis of the different spectral features (see Discussion). According to the results of this comparison, we can suggest a criterion to use low-loss DSF to disentangle the contribution of rutile from anatase, thus permitting to distinguish between the two.
High energy resolution spectra (Eres ∼ 300 meV) were collected for energy losses below 16 eV for rutile and anatase-based pigments (1, 12–16). The differences observed experimentally were most notable at q ∼ 2.6–2.9 Å−1, where a broad peak at ∼8 eV was observed for the two anatase-based pigments, 12 and 13, only (see Fig. S3.3 and S3.4, ESI†). The DSF calculations show good agreement with the experimental spectra for both pigments in terms of the relative intensities of the plasmon peak at ∼14 eV and interband transitions O-2p → Ti-3d at ∼8 eV (Fig. 3). Both theory and experiments agree in highlighting that the peak at 8 eV is the most discriminant between anatase and rutile phases.
Furthermore, this peak is broader, observed as a shoulder in pigments 4–6 Fe d3. Compton scattering dominates in pigments containing light elements, and no spectral features are observed in synthetic organic samples, as expected. Although outside the scope of this contribution, IXS spectroscopy would enable the study of their speciation by extending the energy loss range to the C and O edges.26 Additional spectral features were observed that could not be identified at this stage. Indeed, for many samples, references in low energy-loss ranges were not available in the EELS databases, such as the one from Hitchcock and Mancini,27 nor mentioned in the literature; the identification of these resonances would require the acquisition of new reference data and/or future simulations.
We observed a decrease in the intensity of the plasmonic peak of 1 (titanium white, rutile-based TiO2) with increasing q from ∼2 Å−1, until it was no longer observed at ∼4 Å−1. This phenomenon results from q exceeding the critical value qc for TiO2, where the plasmon merges with the particle–hole continuum that dominates the DSF for a wide range of frequencies. The same behaviour was reported for many materials, including silicon28 or SrVO3.29 IXS relies on the energy and momentum transfers measured from a high-energy incident photon in the X-ray range that is defined by its wave vector, energy, and polarisation (k1, ħω1, ε1) by the electron system and the scattered photon (k2, ħω2, ε2) defined by the same parameters.30,31 The measured energy transfer ħω is given by the law of conservation of energy,
ħω = ħω1 − ħω2 |
The IXS spectrum is accurately represented by the S(q,ω) of a material, i.e., encompassing all the observed excitations in the system, for a wide range of energies and momenta. The DSF contains multiple pieces of information, ranging from core excitations and edges to collective plasmons, from excitons to the Compton regime. It does not suffer from the limitations of other X-ray spectroscopy techniques, like the low signal for low-Z elements, or the necessity to have a perfect crystal. In addition, being a scattering technique, IXS can (i) go beyond the dipole regime, giving access, e.g., to d–d or orbital excitations and (ii) be extremely bulk-sensitive, with the use of hard X-rays.32 Some known limitations include the difficulties at low energy loss (typically in the phonon excitation regime). Here, we particularly want to emphasise the dispersion analysis, made possible by acquiring spectra with varying q, which opens the way to an entirely additional dimension for the discrimination of different structures.
Theoretical ab initio calculations of IXS have been used for the study of materials found in pigments, such as titanium dioxide (TiO2); these materials are chosen due to the simplicity of their structure and their wide range of applications in industrial methods.33 Past results for titanium dioxide using ab initio quantum chemistry approaches showed the quality and accuracy of the methods for rutile at vanishing momentum transfer, in comparison with EELS experiments.34 Here we extend this analysis to more polymorphs and to high momentum, in order to compare with IXS.
The dynamical structure factor can be related to the inverse dielectric function
![]() | (1) |
Theoretical and computational details are shown in Text S6 (ESI†). Here, we can discuss the results and analyse the experiment–theory comparison in more detail. First, by focusing on a specific momentum transfer, e.g., for rutile at q = 2.03 Å−1, we can identify the contributions to the spectrum, via the imaginary part of the dielectric function, Im(ε), that is directly related to the interband transitions, and so to band-structure calculations35 (Fig. 4). We can distinguish four zones: (i) the 0–13 eV range originates from transitions close to the Fermi energy, involving O 2p states, slightly hybridised with Ti 3d's; (ii) excitation of oxygen 2s states is responsible for the region 18–30 eV and (iii) Ti 3p excitations are involved above 35 eV; (iv) the sharp peak at 14 eV serves as a reminder of the plasmon peak at zero-momentum transfer (see Fig. S8.2, ESI†) for the full dispersion DSF, from 0 to high momentum transfer. This analysis permits us to identify all the features in the experimental spectra. Furthermore, our calculations are not affected by the elastic peak, which means that we can unambiguously identify even barely visible shoulders in the experimental peak, like at 5–10 eV in Fig. 1. Here we do not report band structure calculations, band gaps, and bandwidths, that are in line with the already published results.35–39 The DSF is also a pertinent tool to distinguish the anatase from rutile phase, but only if we know, in advance, the q and energy loss range to consider. Unlike optical properties, in fact, which clearly distinguish rutile and anatase (an absorption spectrum strongly depends on the optical bandgap,35 for instance), the DSF spectra of the different pigments might appear, at first glance, very similar. Theory can help us to explain the reason. The two allotropic structures are overall very similar, with the same number of atoms in the unit cell, and the same coordination for the titanium centres. If the experiment is conducted on powders, in addition, the anisotropy of the sample cannot be revealed. A comprehensive study of the dispersion of the DSF with the momentum transfer over a wide range of energy can suggest where the contrast between rutile and anatase signals is higher (see also Fig. S5.1, S6.1 and S6.2, ESI†). The small differences between theory and experiments, notably the missing double-plasmon peak in anatase are due to the particular sample used (powder) here and are explained by highlighting the role of the local field effects (Text S9, ESI†).
Despite a low cross-section, recent instrumental advancements now enable the collection of 2D or 3D volumes, providing valuable insights into speciation and chemical identification in real-life systems.40,41 Noteworthy progress has been achieved in X-ray Raman spectroscopy for the analysis of cultural heritage materials.26,42 Specifically, spectral variations observed in transitions at the carbon K-edge (280–320 eV) could be used to segment voxels containing diverse organics within the bulk of palaeontological samples.43
In the case of 3d t-metals with five or fewer unoccupied states, a broad peak representing the collective excitations – a combination of band structure effects – is observed between 20–30 eV.25 Theoretical calculations of the dynamical structure factor show the broadening of his peak with an increase in the number of electrons in the 3d band, due to the increase the number of electrons available for transitions from 3d states to the unoccupied 3p states of the t-metal structure.25 This seems to be in agreement with the spectral results where the band around 25 eV of TiO2-rutile 1 is seen more clearly as a broad peak and for CoAl2O42 and the Fe-based pigments 4–7, this is observed more as a shoulder covered by a broad band from 30–40 eV. For ZnO 3, the absence of this band is observed, given its nature as a d10 metal where no electrons are available for transition from the occupied t-metal 3d bands to the unoccupied 3p bands.
High-dose X-ray irradiation in the hard and soft X-ray domains can induce damage to paint samples.44,45 In our experiments, we did not observe any photodegradation of the unsupported pigments, either through optical effects or signal modification during data collection. This is particularly remarkable for the (sensitive) organic samples, which did not evolve towards a graphitic signature. Although most significant changes were identified for pigments associated with binding media,46,47 alteration was also reported for unsupported powder pigments, particularly when using micro-focused low-energy X-ray beams.48,49 The absence of alteration observed in our experiments is attributed to the use of a large beam, with correspondingly low fluences (Table S2.1, ESI†) and high energy X-rays, which interact significantly less with mater than low energy photons because of the reduced photoelectric cross-section. However, a colour change was observed post-analysis on the raster-scanned paint. Future studies investigating paints should take into account the increased risks due to the use of smaller beams and the presence of a binder that can play a critical role in redox exchanges. For example, studies on the photoreactivity of pigments, such as TiO2, showed the effect of organic media.50–54 As a highly sensitive probe of chemical speciation, IXS spectroscopy could also be a useful approach for studying the mechanisms induced by irradiation in situ and quantitatively.
![]() | ||
Fig. 5 2D spectral image of a canvas painted with 1, 7, 9, 10, and 11 and image segmented by unsupervised K-means clustering based on their IXS spectra (energy losses: 2–25 eV, q: 5.3–5.5 Å−1). |
Inelastic X-ray scattering is a powerful technique that can be applied to a wide range of sample types. It is applicable not only to powders (like the pigments in our study), but also to monocrystals (as we report in the Text S7, ESI†) and even liquids. The method can provide insights under various temperature and pressure conditions. Its broad applicability makes it a valuable tool for a wide spectrum of scientific investigations. In this work, we supplemented the experimental measurements with theoretical calculations using an ab initio approach based on TDDFT. This predictive tool does not rely on experimental results or on models, and permits a full description of the optical and dielectric properties of a vast range of materials with high accuracy. In conclusion, we showcase a powerful and cutting-edge integration of experimental and theoretical methodologies to study pigments. Simultaneously, we provide a practical illustration that underscores its applicability to a broader spectrum of complex materials.
Footnote |
† Electronic supplementary information (ESI) available: manufacturer and chemical information of commercial pigments, experimental configurations, additional experimental and simulated DSF spectra for anatase and rutile-based pigments, anisotropy of rutile and anatase, theoretical approach, computational details of the numerical calculations of TiO2, experimental and calculated IXS spectra of TiO2 crystals, dispersion of bulk rutile and anatase along the Z direction, comparison of the experiment with theory, and the role of local field effects. See DOI: https://doi.org/10.1039/d3cp04753a |
This journal is © the Owner Societies 2024 |