**Research article**| 22 Nov 2022

# Mechanical ordering of pigment crystallites in oil binder: can electron paramagnetic resonance reveal the gesture of an artist?

Elise Garel, Laurent Binet, and Didier Gourier

**Elise Garel et al.**Elise Garel, Laurent Binet, and Didier Gourier

- Chimie-ParisTech, PSL University, CNRS, Institut de Recherche de Chimie-Paris, 75005 Paris, France

- Chimie-ParisTech, PSL University, CNRS, Institut de Recherche de Chimie-Paris, 75005 Paris, France

**Correspondence**: Laurent Binet (laurent.binet@chimieparistech.psl.eu)

**Correspondence**: Laurent Binet (laurent.binet@chimieparistech.psl.eu)

Received: 05 Jul 2022 – Discussion started: 08 Jul 2022 – Revised: 30 Sep 2022 – Accepted: 24 Oct 2022 – Published: 22 Nov 2022

Is it possible to reconstruct the gesture of an ancient artist applying a paint layer, considering that the orientation distribution of crystallites of an inorganic pigment remains definitively imprinted on the support after drying of the layer? If the pigment contains paramagnetic transition metal ions whose magnetic interactions are themselves anisotropic, then the shape of the electron paramagnetic resonance (EPR) spectrum should reflect the distribution of grain orientations. We have demonstrated this effect in the case of Egyptian blue CaCuSi_{4}O_{10}, a pigment used for at least 3 millennia in antiquity, by reconstructing the probability density of crystallite orientations under various modes of application, such as brush painting, dabbing and droplet deposition.

- Article
(7901 KB) -
Supplement
(1796 KB) - BibTeX
- EndNote

In magnetic resonance, including electron paramagnetic resonance (EPR), the interactions contained in the spin Hamiltonian are essentially anisotropic. For a system in a rigid state, much information on the structure at the molecular scale of the observed species can be obtained from the anisotropic parts of these interactions. The anisotropic interactions can also be helpful for analyzing the orientation distributions of paramagnetic species in disordered materials, thus providing information on the texturation of the sample at a macroscopic scale (Hentschel et al., 1978; Friesner et al., 1979). For instance, EPR has already been used to investigate ordering in liquid crystals (Meirovitch et al., 1982; Imrie et al., 1997; Yankova et al., 2013; Bogdanov and Vorobiev, 2022), in polymer films (Vorobiev and Chumakova, 2005), in muscle fibers (Fajer, 1994), in oriented bacteria (Frank et al., 1979) or in graphene oxide membranes (Chumakova et al., 2022). This approach should also be fruitful in the field of cultural heritage, where complex composite magnetic materials are often to be found in cultural artifacts. As an example, paintings are generally made of layers of
pigment grains dispersed in a polymer binder. Most often pigments contain
paramagnetic transition ions as coloring species. It can be expected that
the orientation of the inorganic grains within the layer is not random but
rather keeps the memory of the painter's gesture. Therefore, analyzing the
orientation distribution of pigment grains within a painting could potentially provide useful information on the making. In this paper, we aim
at testing the potential of EPR in determining the orientation distributions
of pigment grains in a polymer film. The chosen system is cuprorivaite
CaCuSi_{4}O_{10}, also known as Egyptian blue, a pigment widely used in the Mediterranean basin since about 2500 BC (4th Egyptian Dynasty) until the end of the Roman Empire, its manufacturing recipe having been lost around the 7th century AD (Pagès-Camagna and Colinart, 2003). In this work we study by EPR the orientation effects of cuprorivaite crystallites dispersed in dried linseed oil binder. Different deposition modes of the mixture on a substrate were tested to investigate their influence on the orientation distribution.

## 2.1 Sample preparation

The Egyptian blue pigment was synthesized by solid-state reaction from calcium carbonate CaCO_{3}, amorphous silica SiO_{2}, copper oxide CuO
and 3 wt % of sodium carbonate Na_{2}CO_{3}. The powders were mixed, ground together and then pressed into pellets under a uniaxial pressure of 4 t cm^{−2}. The pellets were first sintered in air at 1000 ^{∘}C for 5 h. The resulting samples were then ground, pressed into new pellets, and sintered again for 17 h in air at 1000 ^{∘}C. After this final sintering, the samples were ground again in powder. The purity of the material was tested by X-ray diffraction (XRD) using an X Panalytical X'Pert Pro diffractometer with the *K*_{α1} ray of a copper anticathode (*λ*=0.15406 nm). XRD patterns show that the synthesized powder consists of cuprorivaite CaCuSi_{4}O_{10} with traces of wollastonite CaSiO_{3} and SiO_{2} (Binet et al., 2021). Cuprorivaite is a member of the phyllosilicate group, all minerals of which have the form of platelets. Its structure is made up of double layers of corner-sharing [SiO_{4}] tetrahedra separated by a layer of Ca^{2+} ions. It is tetragonal, with space group P4${/}_{\mathrm{nnc}}$ (Pabst, 1959), which means that the crystallographic *c* axis (*C*_{4} axis) is normal to the silicate layers and thus to the platelets. Cu^{2+} ions occupy plane-square sites with
*D*_{4h} point symmetry (Pabst, 1959), with the *C*_{4} axis (which is also
the ${g}_{\left|\right|}$ axis) parallel to the crystallographic *c* axis.

The as-obtained Egyptian blue pigment was then mixed with boiled linseed oil as a binder in a proportion of 10 wt % of pigment. The liquid pigment–oil mixture was deposited on horizontal supports to allow drying. Either a glass plate or a semi-rigid transparent polymer substrate (transparency for an overhead projector) were used. Four different deposition methods were tested (Fig. 1): (i) a one-way spreading with a film applicator to produce a 250 µm-thick (before drying) film, (ii) a one-way spreading with a paint brush, (iii) a dabbing with a paint brush, and (iv) a droplet deposition with a pipette. The samples were then left to dry for several days until the linseed oil became rigid. Of relevance for further discussion is the fact that, during the drying, the vertical direction, and hence gravity, was perpendicular to the sample plane. For EPR analysis and for samples which had to be rotated about an in-plane direction, small pieces with dimensions 25 mm × 2–3 mm were cut perpendicularly to the spreading direction (when relevant), as shown in Fig. 1. For samples which had to be rotated about the normal to the plane, pieces of approximate size 2 mm ×2 mm were cut to fit into standard EPR tubes.

## 2.2 EPR analysis

EPR spectra were recorded at room temperature with a Bruker Elexsys E500 spectrometer equipped with a Bruker SHQE resonator operating at X-band (microwave frequency ≈ 9.25–9.39 GHz, depending on the sample). The applied magnetic field was modulated at 100 kHz with a 1 mT amplitude for lock-in detection. The microwave power was 2 mW. To analyze theoretically the orientation distribution of pigment grains in a fluid medium, a sample made of pigment powder in fluid linseed oil was submitted to a 20 mT magnetic field in the EPR spectrometer at room temperature. Then the sample was cooled down to 200 K with gaseous nitrogen to freeze the oil and to keep the grain orientations fixed in the sample. EPR spectra at different orientations of the sample in the external field were then recorded at this temperature. Datasets are available in Binet et al. (2022a).

All calculations for EPR spectral analysis were performed with Matlab (script available in Binet et al., 2022b). The Easyspin package (Stoll and Schweiger, 2006) was used for EPR spectrum simulation by taking an axial
**g** matrix with principal values ${g}_{x}={g}_{y}=\mathrm{2.055}={g}_{\perp}$ and ${g}_{z}=\mathrm{2.350}={g}_{\parallel}$ for Cu^{2+} in the Egyptian blue pigment. Here the (*x*, *y*, *z*) molecular frame is such that the *z* axis is the *C*_{4} symmetry axis of the Cu^{2+} coordination site with *D*_{4h} symmetry (Ford and Hitchman, 1979). This axis is also the crystal *c* axis of the structure (Pabst, 1959). A Voigt line shape was considered for the EPR transitions, with an empirical dependence of the peak-to-peak width Δ*B* on the angle *θ* between the applied field and the *z* axis as $\mathrm{\Delta}B=\mathrm{0.86}+\sqrt{\mathrm{0.37}{\mathrm{cos}}^{\mathrm{2}}\mathit{\theta}+\mathrm{0.12}{\mathrm{sin}}^{\mathrm{2}}\mathit{\theta}}$ (mT) for both Lorentzian and Gaussian components.

## 3.1 EPR spectrum of a powder sample with the isotropic orientation distribution of the crystallites

The EPR of the Cu^{2+} ion in cuprorivaite can be described by a spin
Hamiltonian including only the Zeeman interaction with the external applied
magnetic field ** B**:

where *μ*_{B} is the electron Bohr magneton, **g** the electron
*g* factor matrix and ** S** the electron spin. The hyperfine interaction with the central nucleus

^{63}Cu or

^{65}Cu is not to be considered in the Hamiltonian as it is averaged out by exchanged interaction between Cu

^{2+}ions (Binet et al., 2021). The resonance field is then determined by the

*g*factor only and its anisotropy. An EPR spectrum is recorded by scanning the applied field at fixed frequency

*ν*of the electromagnetic wave. For the magnetic field making an angle

*θ*with the molecular

*z*axis, an EPR transition occurs at a resonance magnetic field

with *h* the Planck constant. The resonance field thus varies in a range
defined by the boundary values ${B}_{r}^{\parallel}={B}_{r}\left(\mathit{\theta}=\mathrm{0}\right)=h\mathit{\nu}/{g}_{\parallel}{\mathit{\mu}}_{B}$ and ${B}_{r}^{\perp}={B}_{r}\left(\mathit{\theta}=\frac{\mathit{\pi}}{\mathrm{2}}\right)=h\mathit{\nu}/{g}_{\perp}{\mathit{\mu}}_{B}$, with ${B}_{r}^{\parallel}<{B}_{r}^{\perp}$ in the present case. The EPR transition of a Cu^{2+} ion
thus depends on the orientation of the molecular frame $\left(x,y,z\right)$ within the laboratory frame $\left({X}_{\mathrm{0}},{Y}_{\mathrm{0}},{Z}_{\mathrm{0}}\right)$ defined such as ** B**∥

*Z*

_{0}through the angle

*θ*. In the case of a powder sample where crystallites exhibit a perfectly isotropic distribution of orientations, the EPR spectrum is given by summing over all possible orientations

*θ*:

where *f*(*B*−*B*_{r}(*θ*)) is a normalized lineshape function, and $\mathit{\omega}\left(\mathit{\phi},\mathit{\theta},\mathit{\psi}\right)$ is the transition probability, which depends on the orientation, given by the set of Euler angles $\left(\mathit{\phi},\mathit{\theta},\mathit{\psi}\right)$ of the molecular frame within the laboratory frame. A typical EPR spectrum in the case of an isotropic distribution of crystallite orientations and for an axial **g** matrix is shown in Fig. 2. It exhibits a low-field maximum at ${B}_{r}^{\parallel}$ and a baseline crossing close to (but not exactly at) ${B}_{r}^{\perp}$.

## 3.2 Powder EPR spectrum with the preferred orientation of the crystallites

To describe the EPR spectrum of a sample with a non-uniform distribution of
crystallite orientations, three different frames need to be considered, namely, the molecular frame $(x,y,z)$, the sample frame $\left(X,Y,Z\right)$ with the *X*, *Y* and *Z* axes having a defined position with respect to the sample, and the laboratory frame $\left({X}_{\mathrm{0}},{Y}_{\mathrm{0}},{Z}_{\mathrm{0}}\right)$, as shown in Fig. 3a. The orientation of the molecular frame within the sample
frame is specified by a set of three Euler angles $\mathrm{\Omega}=\left(\mathit{\alpha},\mathit{\beta},\mathit{\gamma}\right)$, its orientation within the laboratory frame
by the set $\mathrm{\Omega}{}^{\prime}=\left(\mathit{\phi},\mathit{\theta},\mathit{\psi}\right)$ and the orientation of the sample frame within the laboratory frame by the set ${\mathrm{\Omega}}_{\mathrm{0}}=\left({\mathit{\alpha}}_{\mathrm{0}},{\mathit{\beta}}_{\mathrm{0}},{\mathit{\gamma}}_{\mathrm{0}}\right)$ (Fig. 3b).

The relevant quantity to characterize the non-random orientations is the
orientation probability density 𝒫(Ω) of the crystallites in the sample frame $\left(X,Y,Z\right)$ such that 𝒫(Ω)d*α*sin *β*d*β*d*γ* is the probability of the molecular frame orientation in the sample frame being found in the range $[\mathit{\alpha},\mathit{\alpha}+\mathrm{d}\mathit{\alpha}[\times [\mathit{\beta},\mathit{\beta}+\mathrm{d}\mathit{\beta}[\times [\mathit{\gamma},\mathit{\gamma}+\mathrm{d}\mathit{\gamma}[$. This density 𝒫(Ω) is directly related to the way the pigment was applied on the sample. It is normalized so that ${\int}_{\mathit{\alpha}=\mathrm{0}}^{\mathrm{2}\mathit{\pi}}{\int}_{\mathit{\beta}=\mathrm{0}}^{\mathit{\pi}}{\int}_{\mathit{\gamma}=\mathrm{0}}^{\mathrm{2}\mathit{\pi}}\mathcal{P}\left(\mathrm{\Omega}\right)\mathrm{d}\mathit{\alpha}\mathrm{sin}\mathit{\beta}\mathrm{d}\mathit{\beta}\mathrm{d}\mathit{\gamma}=\mathrm{1}$. For a specific orientation ${\mathrm{\Omega}}_{\mathrm{0}}=\left({\mathit{\alpha}}_{\mathrm{0}},{\mathit{\beta}}_{\mathrm{0}},{\mathit{\gamma}}_{\mathrm{0}}\right)$ of the sample in the laboratory frame, the EPR spectrum is then given by

In the case of an isotropic orientation distribution, 𝒫(Ω) is constant with $\mathcal{P}\left(\mathrm{\Omega}\right)=\mathrm{1}/\mathrm{8}{\mathit{\pi}}^{\mathrm{2}}$ and can be omitted in Eq. (4), which then reduces to Eq. (3). The transition probability *ω*(Ω,Ω_{0}) scales as the *g* value in the case of a field-swept EPR spectrum (Aasa and Vänngård, 1975). When the anisotropy of the *g* factor is small, as is the case here, the transition probability is almost independent of the orientation (Ω,Ω_{0}). It then appears as an irrelevant scaling factor and will be dropped hereafter for the sake of simplification. The determination of 𝒫(Ω) from experimental EPR spectra is essentially adapted from Hentschel et al. (1978). The unknown probability density 𝒫(Ω) is first expanded as a combination of Wigner matrix elements ${D}_{mn}^{\left(l\right)}\left(\mathrm{\Omega}\right)$:

where

with *Y*_{l,m} and *Y*_{l,n} being spherical harmonics and ${\widehat{L}}_{y}$ the *y* component of the orbital angular momentum operator. An essential simplification can be introduced here because all samples happened to have a revolution symmetry about the *Z* axis of the sample frame (see below) and because the *g* factor is axially symmetric about the molecular *z* axis (*C*_{4} axis), meaning that the spectrum is unchanged by any rotation of the crystallites about the molecular *z* axis. The probability density 𝒫(Ω) then no longer depends on *α* and *γ*, but only on *β*. The dependence on *α* and *γ* arises from terms ${p}_{lmn}{D}_{mn}^{\left(l\right)}\left(\mathrm{\Omega}\right)$ in Eq. 5 with non-zero values of *n* or *m*. The general way to get rid of this
dependence is to set ${p}_{lmn}={p}_{l}\times {\mathit{\delta}}_{n,\mathrm{0}}\times {\mathit{\delta}}_{m,\mathrm{0}}$. This selects terms in Eq. (5) with $m=n=\mathrm{0}$ so that

where ${D}_{\mathrm{00}}^{\left(l\right)}\left(\mathrm{\Omega}\right)={P}_{l}\left(\mathrm{cos}\mathit{\beta}\right)$, with *P*_{l}(*x*) being the *l*th-order Legendre polynomial. The orientation probability density 𝒫(Ω) is then fully determined by the set of coefficients *p*_{l} with *l*∈ℕ, which must be obtained from experimental EPR spectra. Detailed calculations given in the Supplement show that the EPR spectrum for a specific orientation ${\mathrm{\Omega}}_{\mathrm{0}}=\left({\mathit{\alpha}}_{\mathrm{0}},{\mathit{\beta}}_{\mathrm{0}},{\mathit{\gamma}}_{\mathrm{0}}\right)$ of the sample in the laboratory frame is given by

From the above equation, the EPR spectrum only depends on the angle
*β*_{0} between the sample *Z* axis and the applied magnetic field. Therefore, only a set of spectra at different *β*_{0} angles is needed to determine the *p*_{l} coefficients. The numerical implementation of the determination of the *p*_{l} coefficients from a set of experimental spectra at different orientations *β*_{0} of the sample in the laboratory frame is given in the Supplement. It should however be noted that since the EPR phenomenon is unchanged by reverting the orientation of the applied field, only *p*_{l} coefficients with *l* even can be determined so that the actual probability density 𝒫(Ω) cannot be fully obtained from EPR.

The determination of the preferred orientation of pigment crystallites is illustrated in the case of cuprorivaite crystallites dispersed in a fluid medium, such as oil, and oriented by an external magnetic field. Under an external field, the magnetic potential energy per unit volume of a crystallite is

where *N* is the number of Cu^{2+} ions per unit volume and 〈*S*_{z}〉 is the thermal average of the Cu^{2+} electron spin component along the magnetic field direction. In a fluid, a pigment crystallite is free to rotate to minimize its potential energy. It appears from Eq. (8) that the
potential energy is minimal for an orientation *θ* of the molecular
*z* axis with respect to the field corresponding to the maximum value of the *g* factor $g\left(\mathit{\theta}=\mathrm{0}\right)={g}_{\parallel}$. Therefore, we expect pigment crystallites to be mostly oriented so that their crystal *c* axis
(which is parallel to the molecular *z* axis) is along the applied field, which here defines the sample *Z* axis. After orientation in the field, the sample was cooled to freeze the oil and thus the grain orientations, and EPR
spectra were recorded upon varying the angle *β*_{0} between the
measuring field (*Z*_{0} direction) and the sample *Z* axis. These spectra are shown in Fig. 4 for different values of *β*_{0}.

For *β*_{0}=0 (measuring field parallel to the orienting field), the
spectral density is mostly concentrated about the lowest resonance field
${B}_{r}^{\parallel}=\frac{h\mathit{\nu}}{{g}_{\parallel}{\mathit{\mu}}_{B}}$. When the angle
*β*_{0} increases up to $\frac{\mathit{\pi}}{\mathrm{2}}$, the EPR spectral density shifts to higher fields and is concentrated about the highest resonance field ${B}_{r}^{\perp}=h\mathit{\nu}/{g}_{\perp}{\mathit{\mu}}_{B}$ at ${\mathit{\beta}}_{\mathrm{0}}=\frac{\mathit{\pi}}{\mathrm{2}}$
(measuring field perpendicular to the orienting field). This angular dependence of the EPR spectra shows that, as expected, crystallites have
been predominantly oriented with their crystal *c* axis along the orienting
field. The orientation probability density 𝒫(*β*) calculated from the EPR spectra is shown in Fig. 5a, where *β* is also here the angle between the crystal *c* axis and the orienting field. The calculated EPR spectra are given in Fig. S1 in the Supplement. The polar representation of
𝒫(*β*) shows a large maximum of 𝒫(*β*) for *β*=0 and can be compared to the probability density for an isotropic distribution of orientations (red circle in the middle of Fig. 5a). This clearly indicates a strongly anisotropic orientation distribution of the pigment grains. Figure 5b shows the *p*_{l} coefficients in the expansion of 𝒫(*β*) in Eq. (6). Order *l*=0 corresponds to the isotropic contribution to
𝒫(*β*) and values *l*≠0 to the anisotropic components. The highly anisotropic character of 𝒫(*β*) is shown by the fact that the isotropic component *l*=0 is not the dominant one, but rather components with *l* from 2 to 10.

## 5.1 Films deposited with the applicator

Films were deposited on a flat substrate with an applicator by a one-way
spreading along the *X* direction in the substrate plane (Fig. 3a) of a mixture of fluid linseed oil and pigment grains. This application process is
expected to be the most reproducible one among those used in this work.
Since an anisotropy along the *X* direction is induced by the process, it is
meaningful to determine whether a preferred orientation of the grains is induced and how the distribution is oriented with respect to the spreading
*X* direction. EPR spectra corresponding to the rotation of the sample about
the sample *Z* axis and about the sample *Y* axis are shown in Fig. 6a and b,
respectively. It turns out that there is no significant change in the
spectrum shape upon rotation about the *Z* axis (Fig. 6a). This means that the orientation distribution has a revolution symmetry about the *Z* axis.

This is confirmed by the polar plot of the corresponding density 𝒫(*β*) (here *β* is the angle from the *X* axis in the *X*–*Y* plane), which is very close to a circle (Fig. 7a), and by the *p*_{l} coefficients close to zero except for *l*=0 (Fig. 8).

Significant changes in the EPR spectrum shape are observed upon rotation
about the *Y* axis, the magnetic field direction moving from the *Z* direction (*β*_{0}=0) perpendicular to the sample plane to the *X* direction $({\mathit{\beta}}_{\mathrm{0}}=\mathit{\pi}/\mathrm{2})$ in the sample plane (Fig. 6b). The low-field part of the spectrum, corresponding to crystallites with the *c* axis parallel to the field, is enhanced at *β*_{0}≈0, while the high field part (crystallites with the *c* axis perpendicular to the field) is enhanced at ${\mathit{\beta}}_{\mathrm{0}}\approx \mathit{\pi}/\mathrm{2}$. This means that crystallites are mostly oriented with their *c* axis perpendicular to the sample plane. This is confirmed by the probability density 𝒫(*β*) being maximum at *β*=0
(Fig. 7b). Calculated EPR spectra are given in Figs. S2 and S3. The
variability in the preferred orientation of the crystallites has been
analyzed in a set of six samples prepared under the same conditions and the
corresponding 𝒫(*β*) densities shown in Fig. 7b with different colors. The average values of the *p*_{l} coefficients over those six samples are plotted in Fig. 8, with the error bars given by their standard deviation. The maximum value $\mathcal{P}\left(\mathit{\beta}=\mathrm{0}\right)=$ 0.04–0.06 of the probability density for the films is lower than in the case of magnetically oriented grains in fluid oil where the maximum value was $\mathcal{P}\left(\mathit{\beta}=\mathrm{0}\right)\approx \mathrm{0.15}$. For films, the highest value of *p*_{l} is at *l*=0, corresponding to the isotropic contribution to 𝒫(*β*), while for magnetically oriented grains in fluid oil the anisotropic contributions at non-zero orders $l=\mathrm{2}-\mathrm{10}$ dominate the isotropic one (Fig. 5b). All these results indicate that the preferred grain orientation in dried films after mechanical
application is less pronounced than under an external magnetic field in
fluid oil.

## 5.2 Preferred orientations with manual deposition modes

The preferred orientation with more operator-dependent modes of deposition
was also investigated, namely, spreading along the *X* direction in the sample plane with a paint brush, dabbing onto the substrate with a paint brush and droplet deposition with a pipette. In each case two or three samples were
prepared and analyzed. All samples exhibited a revolution symmetry about the
*Z* axis, as in the case of films deposited with the applicator. The EPR spectra under rotation about the *Y* axis are given in Figs. S4 to S6. The orientation probability densities are shown in Fig. 9a and b and the corresponding *p*_{l} coefficients (average values over two to three samples and standard deviations as error bars) in Fig. 10.

In all cases, the probability densities 𝒫(*β*) exhibit a maximum value for *β* in the range 0–17^{∘}. This again shows a preferred orientation with the crystal *c* axis close to the sample *Z* axis. However, the maximum values of 𝒫(*β*), about 0.02–0.03, are significantly lower than for films deposited with the applicator. This is in line with the *p*_{l} coefficients for *l*≠0 in Fig. 10, which
are lower than for films deposited with the applicator showing that the
orientation distribution with paint brushing, dabbing or droplet deposition is less anisotropic than when the film is deposited with an applicator. When
using a paint brush, increasing the strength of the gesture (“strong brush painting” data as compared to “brush painting” data in Figs. 9a and 10) does not seem to induce a significant difference in the orientation distribution, as shown by close 𝒫(*β*) curves in Fig. 9a
and almost equal leading *p*_{l} coefficients for $l=\mathrm{2},\mathrm{4},\mathrm{6}$ in Fig. 10. The
orientation distributions in the cases of dabbing and droplet deposition
appear to be slightly less anisotropic than for brush painting because their maximal values of 𝒫(*β*) are closer to the lower bound 0.02 and because they have lower values of the leading *p*_{l} coefficients for $l=\mathrm{2},\mathrm{4},\mathrm{6}$.

## 5.3 Possible cause of the preferred orientation

Whatever the application mode of the oil–pigment mixture, two common features were revealed by the above analyses: (i) the pigment crystallites exhibit a preferred orientation with their crystal *c* axes oriented almost
perpendicularly to the sample plane, and (ii) the orientation distribution has a revolution symmetry about the normal to the sample plane. More surprising is
that these features are observed even when a mechanical stress is applied
along the sample plane (*X* axis) during the deposition (whether with the
applicator or with a paint brush) and that this mechanical stress seems to induce the most pronounced preferred orientations. It must also be mentioned
that, during the film deposition and drying, the gravitational force was perpendicular to the sample plane. A possible explanation for this preferred
orientation is based on the sedimentation of the pigment crystallites within
the oil while the latter is still viscous. Indeed, microcrystals of cuprorivaite exhibit a tabular morphology (Bloise et al., 2016). Then, when
microcrystals fall through the liquid oil under the gravitational force,
they tend to lay down with their largest face horizontal (Fig. 11a).
Consequently, microcrystals adopt a preferred orientation with the normal to
the largest face perpendicular to the sample plane. The piling of the
crystals is not perfect, some crystals being tilted or aggregated in larger
grains. This is mostly the case with the droplet or dabbing deposition. When a horizontal mechanical stress is applied with an applicator or a paint brush, the polymer chains of the oil become aligned along the stress direction (Fig. 11b). This may have two effects: first, a dissociation of the aggregates into separate microcrystals, second a better alignment of the crystal tablets in
the horizontal plane by polymer–crystal surface interactions. This could
explain why the preferred orientation seems to be enhanced when such a
mechanical stress is applied during film deposition.

We have seen that three competing forces can contribute to the preferential
orientation of pigment grains in an organic binder: a magnetic field,
gravity, and a mechanical stress. The orientation in a magnetic field could
be potentially interesting if it is effective in the Earth field. However,
in the present work on Egyptian blue pigment, we could hardly observe the
orientation effect in a field weaker than 20 mT, which is 1000 times higher
than the Earth field. The effect of gravity and mechanical strain can induce a preferential orientation only if pigment grains have a symmetry lower than
cubic, i.e., in the case of platelets, cylinders, or rods. In the case of Egyptian blue, whose crystallites are in the form of platelets, these two
forces orientate the platelets parallel to each other. Therefore, if the
substrate is horizontal, gravity and mechanical stresses will both orient the platelets parallel to the substrate, as shown in the experiments described
in this work (Fig. 11). Nonetheless, these two situations can be
distinguished by a careful analysis of the EPR line shape, as shown above. However, the platelet shape of the crystallites does not allow us to identify
the direction of the painter's gesture because the *z* axis of the platelets is perpendicular to the spreading gesture, whatever the direction of this
movement. In principle, the frescoes were painted on walls and therefore on vertical supports, whereas the pigment was deposited on horizontal supports
in the present work. Unfortunately, it was not possible to test experimentally the effect of verticality because the binding oil used in this work is a viscous liquid that dries very slowly, so that it flows down when the sample is placed vertically.

The ideal case for determining an artist's gesture by EPR would certainly be the case where pigment crystallites have the shape of elongated parallelepipeds, rods, and cylinders, with an additional condition that the
*z* axis of the molecular frame is along the main axis of the crystallite. In the case of a horizontal substrate, it is anticipated that gravity will
randomly orient the *z* axes of the rods in the *X*–*Y* plane of the sample frame, while application of the pigment will orient the rods in the direction of
the gesture. One would then expect to obtain orientation probability
densities that would precisely reflect the artist's intentions.

An important point to consider is the applicability of the method. With
traditional EPR, which uses a closed resonant cavity containing a small
diameter tube, only small fragments of fresco could be analyzed. This
situation often occurs in archeology, where fragments have spontaneously become detached from walls. In some cases, and only with permission, it is
possible to sample fragments from frescoes or paintings that are sufficiently small to fit into EPR tubes. The ideal situation, towards which
we should move, is the possibility of studying an entire object in situ and non-invasively. A low-frequency (355 MHz) EPR spectrometer for detecting and imaging paramagnetic species on a flat surface has been recently developed
by Hornak's team (Switala et al., 2017). However, at such a low field, the anisotropy of the *g* factor cannot be resolved so that the orientation effects described in the present work cannot be observed. In the project that
finances this work (https://anr.fr/Project-ANR-17-CE29-0002, last access: 20 September 2022), it is planned to build a portable EPR spectrometer, working at 5 GHz, based
on planar micro-resonators of the “microstrip” type. In this case, planar
(or nearly planar) and large objects, which is the case with frescoes and paintings, can be analyzed by moving and rotating the spectrometer over the
surface of the object.

This work shows that the simple act of spreading a paramagnetic pigment dispersed in a binder on a surface introduces a preferential distribution of grain orientations on the paint layer that can be detected by EPR. This effect has been demonstrated with Egyptian blue, a pigment that was used for several millennia in antiquity. By analyzing the variation of the EPR line shape as a function of the angle between the magnetic field and the surface of the layer, the orientation probability density of the grains can be determined and was found to deviate significantly from the isotropic distribution. It is possible to differentiate between orientation distributions induced by mechanical spreading and those in which only gravity is at work (droplet deposition, dabbing). This effect has been demonstrated here with a pigment whose grains are in the form of platelets. It would certainly be greater for rod-shaped grains.

The EPR and graph data are available at https://doi.org/10.17632/bwjy757vfs.1 (Binet et al., 2022a).

The Matlab script for calculations is available at https://doi.org/10.17632/bdxrcmst37.1 (Binet et al., 2022b).

Demonstration of Eq. (7); numerical calculation of the coefficients *p*_{l} from experimental
EPR spectra; Matlab script for the calculation of the orientation
probability density; experimental and calculated EPR spectra. The supplement related to this article is available online at: https://doi.org/10.5194/mr-3-211-2022-supplement.

EG prepared the samples and performed the EPR experiments. EG and LB processed the data and performed the calculations. The three authors discussed and interpreted the results. LB and DG wrote the paper. All the authors have read and agreed to this version of the manuscript.

The contact author has declared that none of the authors has any competing interests.

Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

This research has been supported by the Agence Nationale de la Recherche (grant no. ANR-17-CE29-0002-01).

This paper was edited by Hervé Vezin and reviewed by two anonymous referees.

Aasa, R. and Vänngård, T.: EPR signal intensity and powder shapes: A reexamination, J. Magn. Reson., 19, 308–315, https://doi.org/10.1016/0022-2364(75)90045-1, 1975.

Binet, L., Lizion, J., Bertaina, S., and Gourier, D.: Magnetic and New Optical Properties in the UV–visible Range of the Egyptian Blue Pigment
Cuprorivaite CaCuSi_{4}O_{10}, J. Phys. Chem. C, 15, 25189–25196, https://doi.org/10.1021/acs.jpcc.1c06060, 2021.

Binet, L., Gourier, D., and Garel, E.: Data-MR_paper-Garel_Binet_Gourier, V1, Mendeley Data [data set], https://doi.org/10.17632/bwjy757vfs.1, 2022a.

Binet, L., Garel, E., and Gourier, D.: Matlab_script_MR_paper_Garel_Binet_Gourier, V1, Mendeley Data [data set], https://doi.org/10.17632/bdxrcmst37.1, 2022b.

Bloise, A., Abd El Salam, S., De Luca, R., Crisci, G. M., and Miriello, D.: Flux growth and characterization of cuprorivaite: the influence of temperature, flux, and silica source, Appl. Phys. A 122, 650, https://doi.org/10.1007/s00339-016-0198-1, 2016.

Bogdanov, A. V. and Vorobiev, A. Kh.: Orientation distribution of molecules in a smectic liquid crystal with a distorted director geometry, J. Phys. Chem. B, 126, 2573, https://doi.org/10.1021/acs.jpcb.1c09795, 2022.

Chumakova, N. A., Vorobiev, A. Kh., Astvatsaturov, D. A., Yankova, T. S., Rebrikova, A. T., and Matveev, M. V.: Quantitative determination of the orientational ordering of graphene oxide membranes by the spin probe technique, Phys. Chem. Chem. Phys., 24, 4881, https://doi.org/10.1039/d1cp05417a, 2022.

Fajer, P.: Determination of spin-label orientation within the myosin head, P. Natl. Acad. Sci. USA, 91, 937–941, https://doi.org/10.1073/pnas.91.3.937, 1994.

Ford, R. and Hitchman, M.: Single crystal electronic and EPR spectra of
CaCuSi_{4}O_{10}, a synthetic silicate containing copper(II) in a
four-coordinate, planar ligand environment, Inorg. Chim. Acta, 33,
L167–L170, https://doi.org/10.1016/S0020-1693(00)89436-0, 1979.

Frank, H. A., Friesner, R., Nairn, J. A., Dismukes, G. C., and Sauer, K.: Biochim. Biophys. Acta, 547, 484–501, https://doi.org/10.1016/0005-2728(79)90029-X, 1979.

Friesner, R., Nairn, J. A., and Sauer, K.: Direct calculation of the orientational distribution function of partially ordered ensembles from the EPR line shape, J. Chem. Phys., 71, 358–365, https://doi.org/10.1063/1.438078, 1979.

Hentschel, R., Schlitter, J., and Spiess, H. W.: Orientational distributions in partially ordered solids as determined from NMR and ESR line-shapes, J. Chem. Phys., 68, 56–66, https://doi.org/10.1063/1.435473, 1978.

Imrie, C. T., Ionescu, D., and Luckhurst, G. R.: Molecular organization of the polymer backbone in a side group liquid crystal polymer. An ESR investigation, Macromolecules, 30, 4597–4600, https://doi.org/10.1021/ma970365c, 1997.

Meirovitch, E., Igner, D., Igner, E., Moro, G., and Freed, J. H.: Electron-spin relaxation and ordering in smectic and supercooled nematic liquid crystals, J. Chem. Phys., 77, 3915–3938, https://doi.org/10.1063/1.444346, 1982.

Pabst, A.: Structures of some tetragonal sheet silicates, Acta Crystallogr., 12, 733–739, https://doi.org/10.1107/S0365110X5900216X, 1959.

Pagès-Camagna, S. and Colinart, S.: The Egyptian green pigment: its manufacturing process and links to Egyptian blue, Archaeometry, 45, 637–658, https://doi.org/10.1046/j.1475-4754.2003.00134.x, 2003.

Stoll, S. and Schweiger, A.: Easyspin, a comprehensive software package for spectral simulation and analysis in EPR, J. Magn. Reson., 178, 42–55, https://doi.org/10.1016/j.jmr.2005.08.013, 2006.

Switala, L. E., Black, B. E., Mercovich, C. A., Seshadri, A., and Hornack, J. P.: An electron paramagnetic resonance mobile universal surface explorer, J. Magn. Reson., 285, 18–25, https://doi.org/10.1016/j.jmr.2017.10.004, 2017.

Vorobiev, A. K. and Chumakova, N. A.: Determination of orientation distribution function of anisotropic paramagnetic species by analysis of ESR spectra angular dependence, J. Magn. Reson., 175, 146–157, https://doi.org/10.1016/j.jmr.2005.04.001, 2005.

Yankova, T. S., Chumakova, N. A., Pomogailo, D. A., and Vorobiev, A. K.: Orientational order of guest molecules in aligned liquid crystal as measured by EPR and UV-vis techniques, Liq. Cryst., 40, 1135–1145, https://doi.org/10.1080/02678292.2013.795621, 2013.

- Abstract
- Introduction
- Experimental procedures
- Theoretical background for orientation distribution analysis
- Orientation distribution of pigment crystallites in fluid oil under the magnetic field
- Orientation distribution of pigment grains in dried films: effect of the application process
- Can EPR reveal the gesture of an artist?
- Conclusion
- Data availability
- Author contributions
- Competing interests
- Disclaimer
- Financial support
- Review statement
- References
- Supplement

^{2+}EPR signal reflects the symmetry of the pigment crystallites, we have shown that the orientation distribution of the grains, as revealed by simulation of EPR spectra, is indeed influenced by the nature of the gesture.

- Abstract
- Introduction
- Experimental procedures
- Theoretical background for orientation distribution analysis
- Orientation distribution of pigment crystallites in fluid oil under the magnetic field
- Orientation distribution of pigment grains in dried films: effect of the application process
- Can EPR reveal the gesture of an artist?
- Conclusion
- Data availability
- Author contributions
- Competing interests
- Disclaimer
- Financial support
- Review statement
- References
- Supplement