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

Abstract 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.


Introduction
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.

Sample preparation
The Egyptian blue pigment was synthesized by solid-state reaction from calcium carbonate CaCO 3 , amorphous silica Published by Copernicus Publications on behalf of the Groupement AMPERE. 212 E. Garel et al.: EPR of a mechanically oriented pigment 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 cornersharing [SiO 4 ] tetrahedra separated by a layer of Ca 2+ ions. It is tetragonal, with space group P4/ 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 || 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.

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 = 2.055 = g ⊥ and g z = 2.350 = g 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 B = 0.86 + 0.37cos 2 θ + 0.12sin 2 θ (mT) for both Lorentzian and Gaussian components.

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 = B r (θ = 0) = hν/g µ B and B ⊥ r = B r θ = π 2 = hν/g ⊥ µ B , with B r < B ⊥ r in the present case. The EPR transition of a Cu 2+ ion thus depends on the orientation of the molecular frame (x, y, z) within the laboratory frame (X 0 , Y 0 , Z 0 ) 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 θ : ) is a normalized lineshape function, and ω (ϕ, θ, ψ) is the transition probability, which depends on the orientation, given by the set of Euler angles (ϕ, θ, ψ) 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 and a baseline crossing close to (but not exactly at) B ⊥ r .

Powder EPR spectrum with the preferred orientation of the crystallites
To describe the EPR spectrum of a sample with a nonuniform distribution of crystallite orientations, three different frames need to be considered, namely, the molecular frame (x, y, z), the sample frame (X, Y, Z) with the X, Y and Z axes having a defined position with respect to the sample, and the laboratory frame (X 0 , Y 0 , Z 0 ), as shown in Fig. 3a. The orientation of the molecular frame within the sample frame is specified by a set of three Euler angles = (α, β, γ ), its orientation within the laboratory frame by the set = (ϕ, θ, ψ) and the orientation of the sample frame within the laboratory frame by the set 0 = (α 0 , β 0 , γ 0 ) (Fig. 3b).
The relevant quantity to characterize the non-random orientations is the orientation probability density P ( ) of the crystallites in the sample frame (X, Y, Z) such that P( )dα sin βdβdγ is the probability of the molecular frame orientation in the sample frame being found in the range This density P( ) is directly related to the way the pigment was applied on the sample. It is normalized so that 2π α=0 π β=0 2π γ =0 P( )dα sin βdβdγ = 1. For a specific orientation 0 = (α 0 , β 0 , γ 0 ) of the sample in the laboratory frame, the EPR spectrum is then given by In the case of an isotropic orientation distribution, P( ) is constant with P( ) = 1/8π 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 fieldswept 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 P( ) from experimental EPR spectra is essentially adapted from Hentschel et al. (1978). The unknown probability density P( ) is first expanded as a combination of Wigner matrix elements D (l) mn ( ): where with Y l,m and Y l,n being spherical harmonics andL 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 P( ) then no longer depends on α and γ , but only on β. The dependence on α and γ arises from terms p lmn D (l) mn ( ) 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 × δ n,0 × δ m,0 . This selects terms in Eq. (5) with m = n = 0 so that where D (l) 00 ( ) = P l (cos β), with P l (x) being the lth-order Legendre polynomial. The orientation probability density P( ) is then fully determined by the set of coefficients p l with l ∈ N, which must be obtained from experimental EPR spectra. Detailed calculations given in the Supplement show that the EPR spectrum for a specific orientation 0 = (α 0 , β 0 , γ 0 ) 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 P( ) cannot be fully obtained from EPR.

Orientation distribution of pigment crystallites in fluid oil under the magnetic field
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 (θ = 0) = g . 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 = hν g µ B . When the angle β 0 increases up to π 2 , the EPR spectral density shifts to higher fields and is concentrated about the highest resonance field B ⊥ r = hν/g ⊥ µ B at β 0 = π 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 P(β) 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 P(β) shows a large maximum of P(β) 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 P(β) in Eq. (6). Order l = 0 corresponds to the isotropic contribution to P(β) and values l = 0 to the anisotropic components. The highly anisotropic character of P(β) 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.

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 cor- responding 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 P(β) (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 (β 0 = π/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 Figure 7. Orientation probability densities P(β) for six films deposited under the same conditions with an applicator calculated from (a) a rotation of the sample about the Z axis (β being the angle from the X axis in the X-Y plane) and (b) a rotation about the Y axis (β being the angle from the Z axis in the X-Y plane). The black circles correspond to an isotropic orientation.
(crystallites with the c axis perpendicular to the field) is enhanced at β 0 ≈ π/2. This means that crystallites are mostly oriented with their c axis perpendicular to the sample plane. This is confirmed by the probability density P(β) 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 P(β) 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 P (β = 0) = 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 P (β = 0) ≈ 0.15. For films, the highest value of p l is at l = 0, corresponding to the isotropic contribution to P(β), while for magnetically oriented grains in fluid oil the anisotropic contributions at non-zero orders l = 2−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.

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 P(β) 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 P(β), 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 sig- Figure 9. Orientation probability densities P(β) for films made by (a) brush painting (full blue curves) and strong brush painting (dashed red curves) and (b) dabbing (full blue curves) and droplet deposition (dashed red curves). The black circles correspond to an isotropic orientation distribution. The sample setting in the laboratory frame was identical to the case of Fig. 6b. nificant difference in the orientation distribution, as shown by close P(β) curves in Fig. 9a and almost equal leading p l coefficients for l = 2, 4, 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 P(β) are closer to the lower bound 0.02 and because they have lower values of the leading p l coefficients for l = 2, 4, 6.

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 Figure 11. Models for the sedimentation of pigment microcrystals in oil leading to (a) poorly oriented crystals upon dabbing or droplet deposition and (b) an enhanced orientation when the oil polymers are stretched by an in-plane mechanical stress in the case of brush painting or deposition with an applicator. the preferred orientation seems to be enhanced when such a mechanical stress is applied during film deposition.

Can EPR reveal the gesture of an artist?
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.