Next Article in Journal
Comparative Investigation on the Emission Properties of RAl3(BO3)4 (R = Pr, Eu, Tb, Dy, Tm, Yb) Crystals with the Huntite Structure
Next Article in Special Issue
A Preview of a Construction of a Crystal Lattice Based on Intermolecular Interactions
Previous Article in Journal
Acknowledgement to Reviewers of Crystals in 2018
Previous Article in Special Issue
Cyclodextrin-Driven Formation of Double Six-Ring (D6R) Silicate Cage: NMR Spectroscopic Characterization from Solution to Crystals
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Determination of the Full 207Pb Chemical Shift Tensor of Anglesite, PbSO4, and Correlation of the Isotropic Shift to Lead–Oxygen Distance in Natural Minerals

1
Department of Chemistry, University of Munich (LMU), Butenandtstr. 5-13, 81377 Munich, Germany
2
Mineralogical State Collection Munich (SNSB), Theresienstr. 4, 80333 Munich, Germany
*
Author to whom correspondence should be addressed.
Submission received: 23 December 2018 / Revised: 9 January 2019 / Accepted: 11 January 2019 / Published: 15 January 2019
(This article belongs to the Special Issue NMR Crystallography)

Abstract

:
The full 207 Pb chemical shift (CS) tensor of lead in the mineral anglesite, PbSO 4 , was determined from orientation-dependent nuclear magnetic resonance (NMR) spectra of a large natural single crystal, using a global fit over two rotation patterns. The resulting tensor is characterised by the reduced anisotropy Δ δ = ( 327 ± 4 ) ppm, asymmetry η C S = 0.529 ± 0.002 , and δ i s o = ( 3615 ± 3 ) ppm, with the isotropic chemical shift δ i s o also verified by magic-angle spinning NMR on a polycrystalline sample. The initially unknown orientation of the mounted single crystal was included in the global data fit as well, thus obtaining it from NMR data only. By use of internal crystal symmetries, the amount of data acquisition and processing for determination of the CS tensor and crystal orientation was reduced. Furthermore, a linear correlation between the 207 Pb isotropic chemical shift and the shortest Pb–O distance in the co-ordination sphere of Pb 2 + solely surrounded by oxygen has been established for a large database of lead-bearing natural minerals.

Graphical Abstract

1. Introduction

For the structural characterization of periodic solids, either in single or multicrystalline form, X-ray crystallography is an invaluable tool [1]. As a complementary technique for the elucidation of structure and dynamics in solids, including amorphous ones, nuclear magnetic resonance (NMR) spectroscopy has become an established method [2,3,4], to the extent that, since 2014, a “Commission on NMR Crystallography and Related Methods” exists in the International Union of Crystallography [5]. In the context of NMR crystallography, many questions (e.g., the determination of asymmetric units, assignment of space groups) may already be answered by considering only the isotropic chemical shift δ i s o of the NMR-observed nuclide [2]. The deceptively simple scalar δ i s o is, however, the result of a contraction (of the isotropic part) of the second-rank chemical shift (CS) tensor:
δ i s o = 1 3 δ 11 + δ 22 + δ 33
Here, δ 11 , δ 22 , and δ 33 are the components of the main diagonal of the CS tensor, which are the ’eigenvalues’ in its principal axes system. The full CS tensor δ reflects the spatial distribution of electrons around the observed nucleus, which ’shield’ it from the external magnetic field to a certain extent. Therefore, a more complete picture becomes available when the full chemical shift tensor is known, which provides possible information on co-ordination, the influence of electron lone pairs, etc. The eigenvalues δ 11 , δ 22 , and δ 33 may be determined by measuring the NMR spectrum of a multicrystalline (’powder’) sample under static conditions, as exemplified by Reference [6]. It is also possible to derive these eigenvalues from the rotational side-band pattern of magic-angle spinning (MAS) spectra [7] to even higher precision [8]. In unrivalled accuracy, however, the eigenvalues of the CS tensor may be extracted from orientation-dependent measurements of single crystals [9,10]. Furthermore, the orientation of the CS tensor in the crystal structure, expressed by the corresponding eigenvectors, is unequivocally only accessible via single-crystal NMR spectroscopy. (In certain cases, indirect determination might be possible from the relative orientation of the CS tensor to the dipolar coupling tensor [11].) In spite of these obvious advantages, single-crystal NMR is comparatively seldom employed for tensor determination. This may be traced to (at least) three problems with this method:
(i) 
Data acquisition and evaluation is rather time-consuming when compared to static or MAS NMR of polycrystalline samples;
(ii) 
tensor eigenvector determination in the crystal frame requires knowledge of the exact orientation of the crystal on the goniometer axis;
(iii) 
single crystals of sufficient size must be available because of the comparatively poor signal-to-noise ratio of NMR spectroscopy.
One purpose of this paper is to show that, in many cases, these difficulties may be overcome. As for problem (i), the amount of necessary work could be greatly reduced by exploiting internal crystal symmetries, i.e., considering the NMR resonances of symmetry-related atoms in the unit cell as an additional dataset acquired for an independent but virtual rotation axis [12,13,14]. The same concept can be used to address problem (ii) by calculating the exact orientation of the rotation axis from NMR data alone [15,16,17]. Regarding problem (iii), many systems are amenable to crystal growth by methods available in a standard laboratory. In addition, natural minerals offer a great supply of single crystalline material already prepared by nature, often of sufficient size for NMR investigations [16,17,18,19]. To illustrate how to tackle the above listed problems (i)–(iii), here we present the determination of the full 207 Pb chemical shift tensor of anglesite, PbSO 4 , using a single crystal from a natural mineral deposit, as pictured in Figure 1a. Anglesite crystallises in the orthorhombic space group Pbnm [20], with the location of the 207 Pb in the unit cell shown in Figure 1b.
NMR spectroscopy of 207 Pb (the only stable lead isotope with a nuclear spin, namely, I = 1 / 2 ) is a valuable analytical tool in solid-state research. As a ‘heavy’ nuclide, lead has a very large chemical shift range, extending over approximately 10,000 ppm, which makes 207 Pb a sensitive reporter of the local electronic state [22,23,24]. Recently, there has been intensified interest in 207 Pb-NMR to characterise organolead halide perovskites, which have been identified as promising materials in photovoltaics [25,26,27]. Frequently, 207 Pb NMR studies are augmented with quantum-mechanical calculations of chemical shifts [16,17,24,27]. While, in general, such calculations have matured to the point where their results are very useful for comparison with experimental data [28], calculations for heavy nuclei are still problematic, even when relativistic effects for the electrons are included [29]. Hence, another purpose of the current work is to add to the comparatively small database of 207 Pb CS tensors available in the literature, and to provide more experimental numbers against which calculation methods may be calibrated. Finally, to further advance the use of 207 Pb-NMR as an analytical tool for structure characterisation, we examine an earlier attempt to correlate 207 Pb isotropic shift δ i s o to interatomic distances between lead and the co-ordinating oxygen atoms [22]. In this previous work, a connection was proposed between δ i s o and the average distance of all oxygens in a co-ordination sphere, which, however, is not unequivocally defined. Here, we suggest to instead use the shortest distance between lead and oxygen, which is not only a well-defined parameter, but also improves the observed correlation.
Before discussing our results in detail, we briefly outline the basic principles of our main analytical technique, i.e., NMR spectroscopy of spin I = 1 / 2 in a single crystal.

2. Single-Crystal NMR of 207 Pb with Spin I = 1 / 2

2.1. General Definitions

The energy levels of nuclear spin I in an external magnetic field B 0 are quantum mechanically described by Hamiltonians H i ^ , with the following general form [30,31,32]:
H ^ i = C i I ^ T i B 0
Here, C i is a characteristic constant for the i interaction, and I ^ is the spin operator in vector form. B 0 is the magnetic field vector of an external magnetic field B 0 . The second-rank tensor T i , which can be written as a 3 × 3 matrix, connects the two vector quantities, and therefore encapsulates the orientation dependence of the interaction. For spin–spin interactions, such as direct (dipolar) or indirect (J) couplings, Hamiltonians take a form different from Equation (2). However, for the 207 Pb NMR spectroscopy of anglesite, such couplings result only in some unspecified line broadening, and are not evaluated in detail here. Therefore, the NMR behaviour of 207 Pb (with spin I = 1 / 2 , for which no quadrupolar coupling exists) in our anglesite crystal is primarily determined by the Zeeman interaction H ^ Z and its modification by the electron cloud surrounding the nucleus, the chemical shift H ^ C S [30,31,32]:
H ^ N M R = H ^ Z + H ^ C S = γ n I ^ z B 0 γ n I ^ δ B 0
The nuclear Zeeman interaction determines the resonance (or Larmor) frequency ν 0 for a nuclide with gyromagnetic ratio γ n . The relevant interaction tensor T Z according to Equation 2 is simply the identity matrix, which means that ν 0 solely scales with magnetic field strength. In solids, the chemical shift is orientation-dependent, which is reflected by the existence of the CS tensor T C S = δ . Since the antisymmetric part of δ is practically unobservable, only the symmetric part is considered [33]. The CS tensor, thus, takes the following form in the laboratory co-ordinate system ( L A B ), where the B 0 field vector defines the z axis of a Cartesian co-ordinate system ( x y z ):
δ LAB = δ x x δ x y δ x z δ x y δ y y δ y z δ x z δ y z δ z z
For the NMR spectroscopy of single crystals, it is useful to define two other co-ordinate systems, i.e., the frame of the crystal lattice ( C R Y ), which here is the a b c system of the orthorhombic unit cell of anglesite, and the principal axis system ( P A S ) of the tensor. The CS tensor may be transformed between these three co-ordinate systems using the Euler angles Ω i j relating one frame to the other [34,35]:
δ LAB Ω L C Ω C L δ CRY Ω C P Ω P C δ PAS
Symmetric tensors take diagonal form in their own P A S , where the tensor eigenvalues ( δ 11 , δ 22 , δ 33 ) are the diagonal tensor elements:
δ PAS = δ 11 0 0 0 δ 22 0 0 0 δ 33 = δ i s o Δ δ 2 ( 1 + η C S ) 0 0 0 δ i s o Δ δ 2 ( 1 η C S ) 0 0 0 Δ δ + δ i s o
In the notation on the right-hand side, we used the definition of the isotropic chemical shift δ i s o (Equation (1)), and two alternative tensor parameters according to the Haeberlen convention [31], namely, asymmetry parameter η C S and reduced anisotropy Δ δ :
η C S = δ 22 δ 11 Δ δ ; Δ δ = δ 33 δ i s o
For the above definitions to work properly, the eigenvalues of δ need to be sorted according to:
| δ 33 δ i s o | | δ 11 δ i s o | | δ 22 δ i s o |

2.2. Single-Crystal NMR of Anglesite, PbSO 4

When evaluating the single-crystal NMR of anglesite, it is advantageous to determine the CS tensor of 207 Pb in the orthorhombic crystal frame a b c since the crystallographic symmetry requirements markedly simplify the tensor. The lead atoms occupy Wyckoff position 4c in the crystal structure, which means that they are located on mirror planes parallel to a b , and generated by glides parallel to a c and b c (see also Figure 1b). Since translational elements and inversion do not affect NMR spectra, the four symmetry-related lead atoms in the unit cell form two groups of magnetically inequivalent 207 Pb, which can be considered to be related by a 180 ° rotation about the crystallographic b axis:
Pb ( 1 ) ( x , 1 4 , z ) 180 ° Pb ( 2 ) ( x + 1 2 , 1 4 , z ¯ + 1 2 )
Their respective CS tensors in the C R Y frame are described by only four independent tensor components, P, Q, R, and S:
δ Pb ( 1 ) CRY = P Q 0 Q R 0 0 0 S δ Pb ( 2 ) CRY = P Q 0 Q R 0 0 0 S
To determine the tensor components P, Q, R, and S from our NMR experiments, we now consider what happens when a chemical-shift tensor is rotated in the L A B frame around an axis g , perpendicular to the external magnetic field B 0 (described by normal vector b 0 here), by an angle φ . To this end, a single crystal of anglesite (approximate size 3 × 3 × 2 mm) was fixed on a rod and installed in a goniometer mechanics, as pictured in Figure 2a. The mechanics allows a defined change of rotation angle φ around the goniometer axis g . For such a rotation, the orientation dependence of the resonance frequency (in units of ppm, i.e., scaled by Larmor frequency ν 0 ) is given by Volkoff harmonics [9]:
ν n ( φ ) / ν 0 = b 0 T · δ LAB ( φ ) · b 0 = A n + B n cos 2 φ + C n sin 2 φ
Plots of the resonance positions over the rotation interval of 180 ° , also known as ‘rotation patterns’, are shown in Figure 3. The contribution of each magnetically inequivalent nuclide n in these rotation patterns follows the type of harmonic function shown in Equation (11). The factors A n , B n , and C n are linear combinations of CS tensor components, their precise form depending on the relative orientation of the P A S frame to the L A B frame. In the general case, three rotation experiments about three nonparallel axes with known orientation are required to derive the full CS tensor [9,10]. However, anglesite is one of the cases where internal crystal symmetry may be exploited [12,13,14] to obtain two rotation patterns from just one rotation experiment, corresponding to the two magnetically inequivalent 207 Pb atoms shown in Equation (9). The use of internal symmetry makes the C R Y frame the preferred choice for carrying out all calculations needed to determine the CS tensors of Equation (10). In this frame, the magnetic field vector b 0 has a general orientation depending on the rotation angle φ , and the line position observed for Pb(1) can be expressed as:
ν P b ( 1 ) ( φ ) / ν 0 = b 0 T ( φ ) · δ Pb ( 1 ) CRY · b 0 ( φ ) = P b x 2 + R b y 2 + S b z 2 + Q 2 b x b y
The components of b 0 in the C R Y frame can be described in terms of the (initially unknown) orientation of goniometer axis g :
g = sin θ g cos ϕ g sin θ g sin ϕ g cos θ g
To describe the stepwise rotation of b 0 around g , it is useful to define two auxiliary unit vectors, v and u , in the plane perpendicular to g . To unequivocally set v and u , a reference vector is needed (that obviously needs to be nonparallel to g ), the choice of which determines the definition of offset angle φ Δ in the plane, as shown in Figure 2b. For our data, we chose ( 0 , 0 , 1 ) , i.e., the unit vector along the crystallographic c-axis, so that v and u are given by [15]:
v = 1 sin θ g g × c ; u = v × g = 1 sin θ g g × c × g = 1 sin θ g c g cos θ g
The orientation of b 0 in the C R Y frame may then be written as:
b 0 ( φ ) = v sin φ φ Δ + u cos φ φ Δ
Using this formalism, we can now proceed to determine the chemical shift tensor δ CRY for 207 Pb in anglesite, PbSO 4 .

3. Results and Discussion

3.1. 207 Pb NMR of Anglesite: Single Crystal

From fitting the experimental line positions in Figure 3, we obtained two sets of prefactors ( A 1 , B 1 , C 1 and A 2 , B 2 , C 2 ) from the Volkoff harmonics, see Equation (11). Since the rotation was carried out around one physical goniometer axis g , the two harmonics of the magnetically inequivalent lead atoms are linked by one constraint [14], and only five linear independent parameters may be extracted from one rotation pattern. The chemical shift tensor δ CRY , which we want to calculate from the rotation pattern, possesses four independent tensor elements (Equation (10)) due to constraints imposed by crystal symmetry. However, there are also two parameters ( θ , ϕ ) describing the unknown orientation of the goniometer axis in the C R Y frame, plus the unknown initial orientation of the crystal expressed by offset angle φ Δ , which need to be included into the fit routine. Overall, therefore, our system is underdetermined, with seven fit parameters versus five experimental ones. The obvious solution to this problem is to acquire two rotation patterns about two different rotation axes, g 1 (Figure 3b) and g 2 (Figure 3d), yielding enough free parameters to also include the orientations of g 1 and g 2 into the fit. Thus, the two rotation patterns were simultaneously fitted according to Equation (12) (and its symmetry-related counterpart), with the motion of field vector b 0 in the C R Y frame described by Equation (15). The fit converged on a global solution with the following values for the ten parameters (the errors reflect the fit residuals):
g 1 : θ g = ( 76.8 ± 0.4 ) ° ϕ g = ( 149.8 ± 0.2 ) ° φ Δ = ( 17.5 ± 0.2 ) ° g 2 : θ g = ( 126.6 ± 0.2 ) ° ϕ g = ( 228.3 ± 0.7 ) ° φ Δ = ( 178.4 ± 0.6 ) ° P = ( 3873 ± 4 ) ppm Q = ( 152.0 ± 0.5 ) ppm R = ( 3607 ± 3 ) ppm S = ( 3365 ± 1 ) ppm
The diagonalization of the CS tensor defined by the P , Q , R , S values listed above transforms it to its own P A S frame, giving the eigenvalues and eigenvectors shown in Table 1. Inspecting the properties of δ CRY , we find the anisotropy of the CS tensor ( Δ δ = 327 ppm) to be comparatively small, which indicates that the electron lone pair at the lead atoms has predominantly s-character [36]. The asymmetry parameter η C S = 0.529 , on the other hand, reflects the comparatively low symmetry of the polyhedron defined by the surrounding oxygens, with the co-ordination number being 10 (see Section 3.3 below). Both derived tensor parameters agree reasonably well with those previously reported from static spectra of a polycrystalline sample: Δ δ = 320 ppm and η C S = 0.545 [6].
The orientation of the scaled eigenvectors of the CS tensor in the unit cell of anglesite, as listed in Table 1, may be visualised from Figure 4. These orientations are in agreement with the crystal symmetry, with the eigenvector of one of the lesser tensor components ( δ 11 P A S ) aligning exactly along the c axis in the crystal frame. This is a consequence of the symmetry enforced by the mirror plane in the orthorhombic unit cell at Wyckoff position 4c, on which the 207 Pb atoms are situated. This generates CS tensors of the shape shown in Equation (10), where the principal component S does not change when transforming from the C R Y to the P A S frame. The eigenvectors belonging to the two other principal components are free to orient according to the electronic environment, generated chiefly by the surrounding oxygen atoms. The eigenvector with the second largest corresponding eigenvalue, i.e., d 22 , points exactly along the bisector of the next two nearest oxygen atoms with the same Pb–O distance, as may be best seen in the lower-left corner of Figure 4b.

3.2. 207 Pb NMR of Anglesite: Magic-Angle Spinning

For the single-crystal NMR experiments, the isotropic chemical shift is computed from the diagonal elements of δ according to Equation (1), to δ i s o = 3615 ppm. As a useful comparison, δ i s o may also be obtained from MAS NMR spectroscopy. To this end, we acquired the NMR spectra of a polycrystalline sample (prepared by crushing larger anglesite crystals in a mortar) at various spinning speeds, ranging from 22.5 to 12.5 kHz. The 207 Pb MAS NMR spectrum at 22.5 kHz spinning speed, shown in Figure 5a, shows the isotropic band at 3604.5 ppm, deviating significantly from the single-crystal value. It is, however, well known that for many lead-containing compounds, the 207 Pb chemical shift markedly varies with temperature, to the extent that lead nitrate is routinely used for temperature calibration of NMR probes [37,38,39]. Our anglesite sample also already showed a noticeable change in δ i s o within the comparatively small temperature deviations produced by friction effects of the MAS technique. In order to be able to compare the isotropic shift derived from our single-crystal experiments at static (and therefore room-temperature) conditions, we extrapolated the temperature dependence of the isotropic shift to zero spinning speed, taking into account that the resonance position quadratically varies with MAS frequency [6]. The fit (see Figure 5b) resulted in an ‘static’ isotropic chemical shift of 3614 ± 3 ppm, in near perfect agreement with the single-crystal NMR value. Our δ i s o for anglesite is also in good agreement with the value of 3613 ppm, published in 1997 [22]. Curiously, the isotropic shift determined from a static polycrystalline sample [6] is significantly different at 3505 ppm. All principal tensor components given in that work are about + 110 ppm off ours, so that δ i s o is strongly affected, but derived parameters Δ δ and η C S are not (see above). Since the obvious pitfall of referencing can be ruled out, we currently have no convincing explanation for this large difference in tensor values.
As already mentioned in the introduction, it is also possible to extract all three eigenvalues of the CS tensor from the rotational side-band pattern of MAS spectra. We have performed a Herzfeld–Berger analysis [7] on our MAS spectrum measured at the slowest speed ( 12.5 kHz). The results are listed in Table 1, and agree remarkably well with the single-crystal values. The comparatively large errors of the MAS-derived values could be improved by investing extensive amounts of measurement time to acquire spectra at lower spinning speeds and with better signal-to-noise ratio, but that is outside the scope of this work.

3.3. Correlation of 207 Pb Isotropic Shifts and Pb–O Distances

In the context of NMR crystallography, the aspiration is always to find connections between NMR parameters and structural features of the compounds under investigation. As already mentioned in the introduction, it has previously been attempted to correlate the 207 Pb isotropic shift δ i s o to interatomic distances between lead and oxygen for compounds where lead occurs solely in oxygen co-ordination. In this work by Fayon et al. in 1997 [22], a connection was proposed between δ i s o and the mean Pb–O distance in a given co-ordination sphere. However, this average distance is obviously strongly dependent on the chosen co-ordination sphere and the co-ordination number (CN) following from it. In fact, the question of how to define a co-ordination sphere around a given atom in a periodic structure is complex [40], and still a matter of debate among inorganic chemists. When checking the suggested correlation [22] for our anglesite data, classical counting within the co-ordination polyhedra gives CN = 10 for 207 Pb, whereas the more convoluted concept of effective co-ordination number (eCoN) according to Hoppe [40] produces the almost identical value of 9.8. Thus, for lead in PbSO 4 , the co-ordination sphere is comparatively well-defined, and the calculation of a mean Pb–O distance from the atomic co-ordinates is a straightforward task. However, when we tried to reproduce the plot of Reference [22], we encountered many compounds where the definition of the co-ordination sphere was more ambiguous. For a number of those, our best effort at determining the co-ordination number, and the average distance following from it, gave a result that deviated from the earlier work. The CN’s and corresponding Pb–O distances of our assessment are listed in Table 2, where we restricted ourselves to lead compounds that also occur as natural minerals, with the respective mineral names specified. The plot of δ i s o versus mean Pb–O distance is shown in Figure 6(left), evidently correlating these two parameters only poorly.
However, we found that using the shortest Pb–O distance, which, in contrast to the mean distance, can be established unequivocally, results in much improved correlation with the 207 Pb isotropic chemical shift, as shown in Figure 6(right). Thus, the evolution of δ i s o with the shortest Pb–O distance in the co-ordination sphere, d short , may be described by:
δ i s o [ ppm ] = 23925 102780 · d short [ nm ]
Here, it should be emphasised that the above correlation only applies to Pb 2 + . The absence of the lone pair in Pb 4 + obviously strongly changes the electronic surrounding of the 207 Pb nucleus. As one example, we included the Pb 4 + at the Pb(1) site of the mixed-valence oxide Pb 3 O 4 (minium), which clearly lies far away from the suggested correlation lines in both plots of Figure 6. Regarding the well-documented temperature dependence of the 207 Pb isotropic shift (see above), all δ i s o values in Table 2 were reported at ’room temperature’. While this specification is obviously somewhat vague, we note that the possible errors introduced by temperature deviations of a few degrees Celsius are so small that the error bars would be inside the size of the symbols used in the plots of Figure 6.
Furthermore, in the plot using the shortest Pb–O distances, the data points seem to appear in two clusters. Broadly, this clustering follows chemical intuition in the sense that shorter Pb–O distances are possible for the compounds with lower co-ordination number (for higher co-ordination, the ’crowding’ of oxygen around lead makes close contact difficult). There are, however, exceptions to that rule of thumb, in particular for changbaiite, PbNb 2 O 6 . Compounds with CN = 6 also show up in either cluster. It can therefore not be decided from the currently available database whether the observed clustering is a real effect or purely coincidental.

4. Conclusions

In this work, the full 207 Pb chemical shift tensor in the orthorhombic crystal structure of anglesite, PbSO 4 , has been determined using orientation-dependent NMR measurements on a natural single crystal. With this study, we also attempted to demonstrate that problems usually associated with single-crystal NMR, i.e., problems (i)–(iii) listed in the introduction, may be overcome:
(i) 
The presumed extensive effort for data acquisition and evaluation is clearly manageable for the anglesite system. Despite the relatively low natural abundance of 207 Pb of 22.1%, the spectra for the rotation pattern shown in Figure 3b were acquired with only eight scans and 60 s recycle delay, meaning the absolute measurement time was only about 2.5 hours. It can be seen from Figure 3a that the signal-to-noise ratio of the single-crystal spectra is satisfactory under these conditions. In contrast to static spectra of polycrystalline material (see, for example, Reference [6]), spectral intensity is focused on only few resonance lines, making detection less demanding. Our total experimental time was, of course, longer than 2.5 hours, with reorientation of the crystal and especially the subsequent necessary tuning and matching of the probe causing additional delays (the procedure of reorient/tune/match, however, may in principle be automated, considerably lessening the work load of the operator [56]). To further improve the efficiency of both data acquisition and processing, we utilised crystallographic symmetry relations [12,13,14] between two groups of magnetically inequivalent 207 Pb atoms, reducing the necessary data for determining the full CS tensor to two rotation patterns, about two nonparallel rotation axes. A generous estimate of the optimised total time for data acquisition would therefore be 7–8 h. It took three hours to record the MAS spectrum of a polycrystalline sample shown in Figure 5a, yet, for 207 Pb spectra of systems with large Δ δ , two such MAS spectra are usually required to unequivocally assign the isotropic peak. In terms of pure experimental time, therefore, single-crystal NMR compares not unfavourably to MAS NMR, although the necessary effort for data processing and evaluation is obviously more extensive. At the end of this process, however, the full CS tensor is gained, and not only the value of δ i s o available from fast MAS experiments.
(ii) 
The necessity to know the orientation of the rotation axis in the crystal frame is generally solvable in a number of ways, including optical orientation using well-developed crystal faces. The most precise method is certainly the use of X-ray diffraction on the mounted crystal on which NMR measurements are to be performed. However, the large crystals needed for NMR tend to absorb most of the X-ray beam, in particular when compounds with high lead content, such as anglesite, are being investigated. For these systems, resolving the orientation of the rotation axis by relying only on NMR data is a very useful approach, as shown here and in earlier work [15,16,17].
(iii) 
For the scarce availability of single crystals of sufficient size for NMR experiments, no general solution appears to exist. For the area of 207 Pb-NMR, however, a considerable number of compounds also occur as natural minerals (see Table 2). For many lead compounds, comparatively large single crystals are therefore available from mineral collections.
Furthermore, in the context of establishing connections between NMR parameters and crystal structure, we here have reported on a linear correlation between 207 Pb isotropic chemical shift and the shortest Pb–O distance in the co-ordination sphere, for compounds with Pb 2 + solely surrounded by oxygen. This correlation appears to work better than the previously suggested connection between δ i s o and an average distance of all oxygens in a co-ordination sphere [22], in particular because the precise definition of this sphere is debatable.
As an outlook, the strategy of deriving full chemical shift tensors from single crystals with an initially unknown rotation axis discussed in this paper may also be extended to crystal systems with lower symmetry. In such systems, two rotation axes are insufficient due to the lack of magnetically independent but crystallographically related atoms. Further work along these lines is in progress in our laboratory with the aim to determine chemical-shift eigenvalues and eigenvectors of nuclei with spin I = 1 2 with high precision and without the need for additional crystal-alignment techniques, thus bringing NMR spectroscopy and X-ray diffraction crystallography in close relation.

5. Materials and Methods

207 Pb NMR spectra were acquired on a Bruker Avance-III 500 ( B 0 = 11.7 T) spectrometer at LMU Munich, with a Larmor frequency of ν 0 ( 207 Pb) = 104.63 MHz. A spin–echo sequence [57] was employed to minimise baseline roll. A recycle delay of 60 s was used for both single-crystal and MAS spectra. The rotation pattern of Figure 3b was acquired with 8 scans for each spectrum, the pattern in Figure 3d mostly with 16 scans (one exception is the spectrum at 40 degrees shown in Figure 3c, for which 60 scans were added). The orientation change of the single crystal was realised by a goniometer mechanics build by NMR Service GmbH (Erfurt, Germany), which was ’clipped on’ (see Figure 2a) a wide-bore Bruker static NMR probe with a home-built 6 mm solenoid coil. For the MAS spectra, a polycrystalline sample was prepared by crushing single crystals of anglesite with an agate mortar and measured using a 2.5 mm rotor. All spectra were indirectly referenced to 1 H in 100% TMS at 0.1240 ppm. According to our measurements, this is equivalent to the commonly used reference of a static Pb(NO 3 ) 2 powder spectrum at room temperature, with δ i s o = 3492 ppm, or 3474 ppm at the highest point of the line shape, both in good agreement with previously reported values [6]. The global fit of the two rotation patterns was performed with the program Igor Pro 7 from WaveMatrics Inc., which delivers excellent nonlinear fitting performance.

Author Contributions

Conceptualization, O.E.O.Z., R.H., and T.B.; methodology, O.E.O.Z. and T.B.; formal analysis, O.E.O.Z. and J.S.; investigation, O.E.O.Z. and J.S.; resources, R.H. and T.B.; writing—original draft preparation, O.E.O.Z. and T.B.; writing—review and editing, O.E.O.Z., J.S., R.H., and T.B.; visualization, O.E.O.Z. and J.S.; supervision, T.B.

Funding

This research received no external funding.

Acknowledgments

The authors would like to thank Igor Moudrakovski (MPI-FKF, Stuttgart) for support in both the theory and practice of NMR. The photograph of the goniometer in Figure 2a was provided by Otakar Zeman (Ampfing).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ewald, P.P.; Fifty Years of X-ray Diffraction. International Union of Crystallography. 1999. Available online: http://www.iucr.org/publ/50yearsofxraydiffraction (accessed on 14 January 2019).
  2. Harris, R.K. NMR crystallography: The use of chemical shifts. Solid State Sci. 2004, 6, 1025–1037. [Google Scholar] [CrossRef]
  3. Martineau, C.; Senker, J.; Taulelle, F. NMR Crystallography. Ann. Rep. NMR Spectrosc. 2014, 82, 1–57. [Google Scholar] [CrossRef]
  4. Bryce, D.L. NMR crystallography: Structure and properties of materials from solid-state nuclear magnetic resonance observables. IUCrJ 2017, 4, 350–359. [Google Scholar] [CrossRef] [PubMed]
  5. Available online: https://www.iucr.org/iucr/commissions/nmr-crystallography (accessed on 14 January 2019).
  6. Neue, G.; Dybowski, C.; Smith, M.L.; Hepp, M.A.; Perry, D.L. Determination of 207Pb2+ chemical shift tensors from precise powder lineshape analysis. Solid State Nucl. Magn. Reson. 1996, 6, 241–250. [Google Scholar] [CrossRef]
  7. Herzfeld, J.; Berger, A.E. Sideband intensities in NMR spectra of samples spinning at the magic angle. J. Chem. Phys. 1980, 73, 6021–6030. [Google Scholar] [CrossRef]
  8. Hodgkinson, P.; Emsley, L. The reliability of the determination of tensor parameters by solid-state nuclear magnetic resonance. J. Chem. Phys. 1997, 107, 4808–4816. [Google Scholar] [CrossRef]
  9. Volkoff, G.M.; Petch, H.E.; Smellie, D.W.L. Nuclear electric quadrupole interactions in single crystals. Can. J. Phys. 1952, 30, 270–289. [Google Scholar] [CrossRef]
  10. Vosegaard, T.; Hald, E.; Langer, V.; Skov, H.J.; Daugaard, P.; Bildsøe, H.; Jakobsen, H.J. Improved hardware and software for single-crystal NMR spectroscopy. J. Magn. Reson. 1998, 135, 126–132. [Google Scholar] [CrossRef]
  11. Hou, G.; Paramasivam, S.; Byeon, I.J.L.; Gronenborn, A.M.; Polenova, T. Determination of relative tensor orientations by γ-encoded chemical shift anisotropy/heteronuclear dipolar coupling 3D NMR spectroscopy in biological solids. Phys. Chem. Chem. Phys. 2010, 12, 14873–14883. [Google Scholar] [CrossRef]
  12. Weil, J.A. Use of symmetry-related crystal sites for measuring tensor properties in magnetic resonance. J. Magn. Reson. 1973, 10, 391–393. [Google Scholar] [CrossRef]
  13. Berglund, B.; Tegenfeldt, J. The determination of quadrupole coupling tensors from single-crystal NMR data. J. Magn. Reson. 1978, 30, 451–455. [Google Scholar] [CrossRef]
  14. Tesche, B.; Zimmermann, H.; Poupko, R.; Haeberlen, U. The Single-Rotation Method of Determining Quadrupole Coupling Tensors in Monoclinic Crystals. Error Analysis and Application to Bullvalene. J. Magn. Reson. Ser. A 1993, 140, 68–77. [Google Scholar] [CrossRef]
  15. Bräuniger, T.; Poupko, R.; Luz, Z.; Zimmermann, H.; Haeberlen, U. Quantification of the orientational disorder in ortho-dichlorotetramethylbenzene: A single crystal deuterium nuclear magnetic resonance and x-ray study of the site populations. J. Chem. Phys. 2001, 115, 8049–8059. [Google Scholar] [CrossRef]
  16. Zeman, O.E.O.; Moudrakovski, I.L.; Hoch, C.; Hochleitner, R.; Schmahl, W.W.; Karaghiosoff, K.; Bräuniger, T. Determination of the 31P and 207Pb Chemical Shift Tensors in Pyromorphite, Pb5(PO4)3Cl, by Single-Crystal NMR Measurements and DFT Calculations. Z. Anorg. Allg. Chem. 2017, 643, 1635–1641. [Google Scholar] [CrossRef]
  17. Zeman, O.E.O.; Hoch, C.; Hochleitner, R.; Bräuniger, T. NMR interaction tensors of 51V and 207Pb in vanadinite, Pb5(VO4)3Cl, determined from DFT calculations and single-crystal NMR measurements, using only one general rotation axis. Solide State Nucl. Magn. Reson. 2018, 89, 11–20. [Google Scholar] [CrossRef] [PubMed]
  18. Hafner, S.; Raymond, M. The nuclear quadrupole coupling tensors of Al27 in kyanite. Am. Mineral. 1967, 52, 1632–1642. [Google Scholar]
  19. Bryant, P.L.; Harwell, C.R.; Wu, K.; Fronczek, F.R.; Hall, R.W.; Butler, L.G. Single-crystal 27Al NMR of andalusite and calculated electric field gradients: The first complete NMR assignment for a 5-coordinated aluminum site. J. Phys. Chem. A 1999, 103, 5246–5252. [Google Scholar] [CrossRef]
  20. Antao, S.M. Structural trends for celestite (SrSO4), anglesite (PbSO4), and barite (BaSO4): Confirmation of expected variations within the SO4 groups. Am. Mineral. 2012, 97, 661–665. [Google Scholar] [CrossRef]
  21. Momma, K.; Izumi, F. VESTA 3 for three-dimensional visualization of crystal, volumetric and morphology data. J. Appl. Cryst. 2011, 44, 1272–1276. [Google Scholar] [CrossRef]
  22. Fayon, F.; Farnan, I.; Bessada, C.; Coutures, J.; Massiot, D.; Coutures, J.P. Empirical Correlations between 207Pb NMR Chemical Shifts and Structure in Solids. J. Am. Chem. Soc. 1997, 119, 6837–6843. [Google Scholar] [CrossRef]
  23. Mason, H.E.; Hirner, J.J.; Xu, W.; Parise, J.B.; Phillips, B.L. Solid-state NMR spectroscopy of Pb-rich apatite. Magn. Reson. Chem. 2009, 47, 1062–1070. [Google Scholar] [CrossRef] [PubMed]
  24. Dmitrenko, O.; Bai, S.; Beckmann, P.A.; van Bramer, S.; Vega, A.J.; Dybowski, C. 207Pb The Relationship between 207Pb NMR Chemical Shift and Solid-State Structure in Pb(II) Compounds. J. Phys. Chem. A 2008, 112, 3046–3052. [Google Scholar] [CrossRef] [PubMed]
  25. Senocrate, A.; Moudrakovski, I.; Kim, G.Y.; Yang, T.Y.; Gregori, G.; Grätzler, M.; Maier, J. The Nature of Ion Conduction in Methylammonium Lead Iodide: A Multimethod Approach. Angew. Chem. Int. Ed. 2017, 56, 7755–7759. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Hanrahan, M.P.; Men, L.; Rosales, B.A.; Vela, J.; Rossini, A.J. Sensitivity-Enhanced 207Pb Solid-State NMR Spectroscopy for the Rapid, Non-Destructive Characterization of Organolead Halide Perovskites. Chem. Mater. 2018, 30, 7005–7015. [Google Scholar] [CrossRef]
  27. Senocrate, A.; Moudrakovski, I.; Maier, J. Short-range ion dynamics in methylammonium lead iodide by multinuclear solid state NMR and 127I NQR. Phys. Chem. Chem. Phys. 2018, 20, 20043–20055. [Google Scholar] [CrossRef] [PubMed]
  28. Charpentier, T. The PAW/GIPAW approach for computing NMR parameters: A new dimension added to NMR study of solids. Solide State Nucl. Magn. Reson. 2011, 40, 1–20. [Google Scholar] [CrossRef]
  29. Alkan, F.; Dybowski, C. Chemical-shift tensors of heavy nuclei in network solids: A DFT/ZORA investigation of 207Pb chemical-shift tensors using the bond-valence method. Phys. Chem. Chem. Phys. 2015, 17, 25014–25026. [Google Scholar] [CrossRef]
  30. Mehring, M. Principles of High Resolution NMR in Solids; Springer: New York, NY, USA, 1983; ISBN 3-540-11852-7. [Google Scholar]
  31. Haeberlen, U. High Resolution NMR in Solids: Selective Averaging. In Advances in Magnetic Resonance; Waugh, J., Ed.; Academic Press: New York, NY, USA, 1976; ISBN 0-12-025561-8. [Google Scholar]
  32. Laws, D.D.; Bitter, H.-M.L.; Jerschow, A. Solid-State NMR Spectroscopic Methods in Chemistry. Angew. Chem. Int. Ed. 2002, 41, 3096–3129. [Google Scholar] [CrossRef]
  33. Anet, F.A.L.; O’Leary, D.J. The shielding tensor. Part I: Understanding its symmetry properties. Concepts Magn. Reson. 1991, 3, 193–214. [Google Scholar] [CrossRef]
  34. Rose, M.E. Elementary Theory of Angular Momentum; John Wiley & Sons: New York, NY, USA, 1957; ISBN 0-486-68480-6. [Google Scholar]
  35. Millot, Y.; Man, P.P. Active and passive rotations with Euler angles in NMR. Concepts Magn. Reson. 2012, 40A, 215–252. [Google Scholar] [CrossRef]
  36. Greer, G.J.; Michaelis, V.K.; Katz, M.J.; Leznoff, D.B.; Schreckenbach, G.; Kroeker, S. Characterising lone-pair activity of lead(II) by 207Pb solid-state NMR spectroscopy: coordination polymers of [N(CN)2] and [Au(CN)2] with terpyridine ancillary ligands. Chem. Eur. J. 2011, 17, 3609–3618. [Google Scholar] [CrossRef] [PubMed]
  37. Bielecki, A.; Burum, D.P. Temperature dependence of 207Pb MAS spectra of solid lead nitrate: An accurate, sensitive thermometer for variable-temperature MAS. J. Magn. Reson. A 1995, 116, 215–220. [Google Scholar] [CrossRef]
  38. Van Gorkom, L.C.M.; Hook, J.M.; Logan, M.B.; Hanna, J.V.; Wasylishen, R.E. Solid-state lead-207 NMR of lead(II) nitrate: Localized heating effects at high magic angle spinning speeds. Magn. Reson. Chem. 1995, 33, 791–795. [Google Scholar] [CrossRef]
  39. Milder, T.; Ernst, H.; Freude, D. 207Pb NMR detection of spinning-induced temperature gradients in MAS rotors. Solid State Nucl. Magn. Reson. 1995, 5, 269–271. [Google Scholar] [CrossRef]
  40. Hoppe, R. Effective coordination numbers (ECoN) and mean Active ionic radii (MEFIR). Z. Kristallog. Cryst. Mater. 1979, 150, 23–52. [Google Scholar] [CrossRef]
  41. Leciejewicz, J. On the crystal structure of tetragonal (red) PbO. Acta Crystallogr. 1961, 14, 1304. [Google Scholar] [CrossRef]
  42. Hill, R.J. Refinement of the structure of orthorhombic PbO (massicot) by Rietveld analysis of neutron powder diffraction data. Acta Crystallogr. 1985, C41, 1281–1284. [Google Scholar] [CrossRef]
  43. Boucher, M.L.; Peacor, D.R.Z. The crystal structure of alamosite, PbSiO3. Krist 1968, 126, 98–111. [Google Scholar] [CrossRef]
  44. Gavarri, J.R.; Weigel, D.J. Oxydes de plomb. I. Structure cristalline du minium Pb3O4, á température ambiante (293 K). Solid State Chem. 1975, 13, 252–257. [Google Scholar] [CrossRef]
  45. Nelmes, R.J.; Kuhs, W.F. The crystal structure of tetragonal PbTiO3 at room temperature and at 700 K. Solid State Commun. 1985, 54, 721–723. [Google Scholar] [CrossRef]
  46. Okudera, H. Relationships among channel topology and atomic displacements in the structures of Pb5(BO4)3Cl with B = P (pyromorphite), V (vanadinite), and As (mimetite). Am. Mineral. 2013, 98, 1573–1579. [Google Scholar] [CrossRef]
  47. Moreau, J.M.; Galez, P.; Peigneux, J.P.; Korzhik, M.V. Structural characterization of PbWO4 and related new phase Pb7W8O(32−x). J. Alloy. Compd. 1996, 238, 46–48. [Google Scholar] [CrossRef]
  48. Nolle, A. 207Pb Magnetic Shielding Anisotropy in Pb(NO3)2, PbCO3, PbCrO4, PbMoO4 and PbWO4 by Fourier Transform NMR. Z. Naturforsch. 1977, 32a, 964–967. [Google Scholar] [CrossRef]
  49. Sleight, A.W. Accurate cell dimensions for ABO4 molybdates and tungstates. Acta Crystallogr. 1972, B28, 2899–2902. [Google Scholar] [CrossRef]
  50. Effenberger, H.; Pertlik, F. Four monazite type structures: Comparison of SrCrO4, SrSeO4, PbCrO4(crocoite), and PbSeO4. Z. Kristallogr. 1986, 176, 75–83. [Google Scholar] [CrossRef]
  51. Basso, R.; Lucchetti, G.; Zefiro, L.; Palenzona, A. Rosiaite, PbSb2O6, a new mineral from the Cetine mine, Siena, Italy. Eur. J. Mineral. 1996, 8, 487–492. [Google Scholar] [CrossRef]
  52. Yoko, T.; Tadanaga, K.; Miyaji, F.; Sakka, S. A 207Pb MAS-NMR study of Pb-containing glasses. J. Non-Cryst. Solids 1992, 150, 192–196. [Google Scholar] [CrossRef]
  53. Chevrier, G.; Giester, G.; Heger, G.; Jarosch, D.; Wildner, M.; Zemann, J. Neutron single-crystal refinement of cerussite, PbCO3, and comparison with other aragonite-type carbonates. Z. Kristallog. Cryst. Mater. 1992, 199, 67–74. [Google Scholar] [CrossRef]
  54. Brusset, H.; Gillier-Pandraud, H.; Mahe, R. Structure cristalline du metaniobate de plomb rhomboedrique. C. R. Acad. Sci. C 1966, 263, 217–219. [Google Scholar]
  55. Zhao, P.; Prasad, S.; Huang, J.; Fitzgerald, J.J.; Shore, J.S. Lead-207 NMR Spectroscopic Study of Lead-Based Electronic Materials and Related Lead Oxides. J. Phys. Chem. B 1999, 103, 10617–10626. [Google Scholar] [CrossRef]
  56. Pecher, O.; Halat, D.M.; Lee, J.; Liu, Z.; Griffith, K.J.; Braun, M.; Grey, C.P. Enhanced efficiency of solid-state NMR investigations of energy materials using an external automatic tuning/matching (eATM) robot. J. Magn. Reson. 2017, 275, 127–136. [Google Scholar] [CrossRef] [PubMed]
  57. Kunwar, A.C.; Turner, G.L.; Oldfield, E. Solid-state spin-echo Fourier transform NMR of 39K and 67Zn salts at high field. J. Magn. Reson. 1986, 69, 124–127. [Google Scholar] [CrossRef]
Figure 1. (a) Single crystal of anglesite, PbSO 4 , from Monteponi Mine, Iglesias, Sardinia/Italy (Mineralogical State Collection, Munich inventory MSM 6140). (b) Unit cell of anglesite according to Reference [20]. The lead atoms (gray, located at Wyckoff position 4c) are related by glides ( a 1 4 c ) , ( a 3 4 c ) and ( 1 4 b c ) , ( 3 4 b c ) , and are co-ordinated by the ten nearest oxygen atoms (red). The sulfur atoms (yellow) are located in the ( a b 1 4 ) and ( a b 3 4 ) planes, with their covalent bonds to oxygen in tetrahedral co-ordination also shown. (Drawing generated with the Vesta program [21]).
Figure 1. (a) Single crystal of anglesite, PbSO 4 , from Monteponi Mine, Iglesias, Sardinia/Italy (Mineralogical State Collection, Munich inventory MSM 6140). (b) Unit cell of anglesite according to Reference [20]. The lead atoms (gray, located at Wyckoff position 4c) are related by glides ( a 1 4 c ) , ( a 3 4 c ) and ( 1 4 b c ) , ( 3 4 b c ) , and are co-ordinated by the ten nearest oxygen atoms (red). The sulfur atoms (yellow) are located in the ( a b 1 4 ) and ( a b 3 4 ) planes, with their covalent bonds to oxygen in tetrahedral co-ordination also shown. (Drawing generated with the Vesta program [21]).
Crystals 09 00043 g001
Figure 2. (a) ‘Clip-on’ goniometer built by NMR Service GmbH (Erfurt, Germany) for a wide-bore Bruker static NMR probe, with an arbitrary single crystal mounted on a wooden rotation axis with 5 mm diameter. The goniometer axis (red) is perpendicular to the applied magnetic field B 0 (blue). (b) Co-ordinates used for the interpretation of the rotation patterns (Figure 3b,d); a, b, and c are the unit vectors of the orthorhombic crystal system ( C R Y ) of anglesite, g is the unit vector along the goniometer axis, θ g and ϕ g its polar angles in the C R Y , b 0 is the unit vector of the magnetic field, φ is the rotation angle, and u and v are auxiliary unit vectors necessary to define b 0 perpendicular to g [15].
Figure 2. (a) ‘Clip-on’ goniometer built by NMR Service GmbH (Erfurt, Germany) for a wide-bore Bruker static NMR probe, with an arbitrary single crystal mounted on a wooden rotation axis with 5 mm diameter. The goniometer axis (red) is perpendicular to the applied magnetic field B 0 (blue). (b) Co-ordinates used for the interpretation of the rotation patterns (Figure 3b,d); a, b, and c are the unit vectors of the orthorhombic crystal system ( C R Y ) of anglesite, g is the unit vector along the goniometer axis, θ g and ϕ g its polar angles in the C R Y , b 0 is the unit vector of the magnetic field, φ is the rotation angle, and u and v are auxiliary unit vectors necessary to define b 0 perpendicular to g [15].
Crystals 09 00043 g002
Figure 3. 207 Pb NMR spectra of a single crystal of anglesite, PbSO 4 , rotated counterclockwise by the indicated angle φ around rotation axes g 1 (a) and g 2 (c), both perpendicular to the external magnetic field B 0 (varying signal intensity is due to the limited excitation width of the used spin–echo sequence). Full rotation pattern over 180 ° for the two magnetically inequivalent 207 Pb at Wyckoff position 4c, acquired by rotating the anglesite crystal stepwise by 10 ° around g 1 (b) and g 2 (d). Lines represent the fit of the experimental points to the CS tensor elements (see main text for details).
Figure 3. 207 Pb NMR spectra of a single crystal of anglesite, PbSO 4 , rotated counterclockwise by the indicated angle φ around rotation axes g 1 (a) and g 2 (c), both perpendicular to the external magnetic field B 0 (varying signal intensity is due to the limited excitation width of the used spin–echo sequence). Full rotation pattern over 180 ° for the two magnetically inequivalent 207 Pb at Wyckoff position 4c, acquired by rotating the anglesite crystal stepwise by 10 ° around g 1 (b) and g 2 (d). Lines represent the fit of the experimental points to the CS tensor elements (see main text for details).
Crystals 09 00043 g003
Figure 4. Orientation of the experimentally determined 207 Pb chemical-shift eigenvectors (blue) in the unit cell of anglesite, scaled according to the magnitude of the associated eigenvalues, such that an absolute value of 750 ppm corresponds to a length of 1 Å, with the lead atoms (gray) shown with their co-ordination to ten oxygen atoms (red): (a) View down the c axis; (b) view down the b axis. (Drawings generated with the Vesta program [21]).
Figure 4. Orientation of the experimentally determined 207 Pb chemical-shift eigenvectors (blue) in the unit cell of anglesite, scaled according to the magnitude of the associated eigenvalues, such that an absolute value of 750 ppm corresponds to a length of 1 Å, with the lead atoms (gray) shown with their co-ordination to ten oxygen atoms (red): (a) View down the c axis; (b) view down the b axis. (Drawings generated with the Vesta program [21]).
Crystals 09 00043 g004
Figure 5. (a) 207 Pb magic-angle spinning NMR spectra of polycrystalline anglesite, acquired at the indicated MAS frequencies, with isotropic bands marked by a triangle. (b) Plot of the 207 Pb isotropic chemical shift of anglesite versus the squared MAS frequency, with the linear fit given by the dashed line. An extrapolation to zero spinning speed results in an isotropic chemical shift of 3614 ± 3 ppm.
Figure 5. (a) 207 Pb magic-angle spinning NMR spectra of polycrystalline anglesite, acquired at the indicated MAS frequencies, with isotropic bands marked by a triangle. (b) Plot of the 207 Pb isotropic chemical shift of anglesite versus the squared MAS frequency, with the linear fit given by the dashed line. An extrapolation to zero spinning speed results in an isotropic chemical shift of 3614 ± 3 ppm.
Crystals 09 00043 g005
Figure 6. (Left) Evolution of 207 Pb isotropic chemical shift ( δ i s o ) versus mean Pb–O distance for compounds where lead occurs solely in oxygen co-ordination, with anglesite shown in red. (Right) δ i s o versus the shortest Pb–O distance. All data given in Table 2. (Dashed lines show the least-squares linear fit, with fit parameters and their errors given in the upper-right corners. The data point for Pb 4 + (1) in Pb 3 O 4 (Minium) is shown as a gray triangle and was not included in the least-square fits, as the correlation was restricted to only Pb 2 + ).
Figure 6. (Left) Evolution of 207 Pb isotropic chemical shift ( δ i s o ) versus mean Pb–O distance for compounds where lead occurs solely in oxygen co-ordination, with anglesite shown in red. (Right) δ i s o versus the shortest Pb–O distance. All data given in Table 2. (Dashed lines show the least-squares linear fit, with fit parameters and their errors given in the upper-right corners. The data point for Pb 4 + (1) in Pb 3 O 4 (Minium) is shown as a gray triangle and was not included in the least-square fits, as the correlation was restricted to only Pb 2 + ).
Crystals 09 00043 g006
Table 1. Chemical shift tensor δ CRY of 207 Pb in the anglesite structure. Left: Determined from single-crystal NMR experiments about two rotation axes at room temperature. The orientations of the corresponding eigenvectors are listed in spherical co-ordinates ( θ , φ ) in the orthorhombic a b c crystal system C R Y and refer to the atom closest to the origin, i.e., Pb(1). Error values are derived from the fit residuals. Right: Determined from a Herzfeld–Berger analysis [7] of the rotational side-band pattern at ν r = + 12.5  kHz magic-angle spinning (MAS) (i.e., at slightly elevated temperature), with  I 1 / I 0 = 0.52 , I 2 / I 0 = 0.24 , and I 3 / I 0 = 0.054 , leading to ρ = 0.3 ± 0.04 and μ = 4.5 ± 0.3 . The error values of the tensor components are derived from those of ρ and μ .
Table 1. Chemical shift tensor δ CRY of 207 Pb in the anglesite structure. Left: Determined from single-crystal NMR experiments about two rotation axes at room temperature. The orientations of the corresponding eigenvectors are listed in spherical co-ordinates ( θ , φ ) in the orthorhombic a b c crystal system C R Y and refer to the atom closest to the origin, i.e., Pb(1). Error values are derived from the fit residuals. Right: Determined from a Herzfeld–Berger analysis [7] of the rotational side-band pattern at ν r = + 12.5  kHz magic-angle spinning (MAS) (i.e., at slightly elevated temperature), with  I 1 / I 0 = 0.52 , I 2 / I 0 = 0.24 , and I 3 / I 0 = 0.054 , leading to ρ = 0.3 ± 0.04 and μ = 4.5 ± 0.3 . The error values of the tensor components are derived from those of ρ and μ .
Single-Crystal NMRMAS NMR
δ 11 P A S ( 3365 ± 3 ) ppm ( 3369 ± 15 ) ppm
δ 22 P A S ( 3538 ± 3 ) ppm ( 3557 ± 10 ) ppm
δ 33 P A S ( 3942 ± 4 ) ppm ( 3907 ± 30 ) ppm
d 11 0.00 , 0.00
d 22 90.00 , 65.59
d 33 90.00 , 155.59
Δ δ ( 327 ± 4 ) ppm ( 296 ± 30 ) ppm
η C S 0.529 ± 0.002 0.64 ± 0.06
δ i s o ( 3611 ± 4 ) ppm
( 3615 ± 3 ) ppm ( 3614 ± 3 ) ppm a
a From extrapolation of MAS spectra to zero spinning, see Figure 5b.
Table 2. Crystallographic data and 207 Pb isotropic chemical shifts δ i s o of lead-bearing minerals. The first co-ordination sphere of at least one lead side for each compound is solely determined by the next neighbouring oxygen atoms.
Table 2. Crystallographic data and 207 Pb isotropic chemical shifts δ i s o of lead-bearing minerals. The first co-ordination sphere of at least one lead side for each compound is solely determined by the next neighbouring oxygen atoms.
MineralSpace Group (Reference)Pb Coord. No.Pb–O Mean Bond Length (nm)Pb–O Shortest Bond Length (nm) δ iso (ppm)NMR Reference
red PbO [Litharge] P 4 / n m m [41]40.23090.23091939 ± 5[22]
yellow PbO [Massicot] P b c m [42]40.23580.22211515 ± 5[22]
PbSiO 3 : Pb(1) [Alamosite] Pb(2) [Alamosite] Pb(3) P 2 / n [43]5 5 60.2518 0.2492 0.26330.2291 0.2256 0.2299 166 ± 5 93 ± 5 −366 ± 5[22]
[22]
[22]
Pb 3 O 4 : Pb 4 + (1) [Minium] Pb(2) P 4 2 / m b c [44]6 40.2176 0.23740.2131 0.2215−1105 ± 1 795 ± 2[22]
[22]
PbTiO 3 [Macedonite] P 4 m m [45]12 (8+4)0.28460.2517−1395.1 ± 8.4[6]
Pb 5 (VO 4 ) 3 Cl: Pb(1) [Vanadinite] Pb(2) a P 6 3 / m [46]9 60.2719 0.26640.2483 0.2324−1619 ± 2 −1729 ± 9[17]
[17]
PbWO 4 [Stolzite] I 4 1 / a [47]80.25960.2579−2003[48]
PbMoO 4 [Wulfenite] I 4 1 / a [49]80.26570.2608−2009 ± 2[22]
PbCrO 4 [Crocoite] P 2 1 / n [50]80.26340.2532−2236.2 ± 2.9[6]
PbSb 2 O 6 [Rosiaite] P 3 ¯ m [51]60.25540.2554−2581 ± 3[52]
PbCO 3 [Cerussite] P m c n [53]90.26920.2594−2630 ± 2[22]
Pb 5 (PO 4 ) 3 Cl: Pb(1) [Pyromorphite] Pb(2) a P 6 3 / m [46]9 60.2710 0.27610.2556 0.2353−2810 ± 14 −2172 ± 11[16]
[16]
PbNb 2 O 6 : Pb(1) [Changbaiite] Pb(2) [Changbaiite] Pb(3) R 3 m [54]6 6 (3+3) 12 (6+6)0.2577 0.2962 0.27380.2560 0.2571 0.2343−2829 −2829 134[55]
[55]
[55]
PbSO 4 [Anglesite] P b n m [20]100.27850.2612−3615 ± 3this work
a Pb(2) is located on Wyckoff position 6h in vanadinite and pyromorphite, and is not included into the Pb–O distance correlation due to chlorine also being in the first co-ordination sphere.

Share and Cite

MDPI and ACS Style

Zeman, O.E.O.; Steinadler, J.; Hochleitner, R.; Bräuniger, T. Determination of the Full 207Pb Chemical Shift Tensor of Anglesite, PbSO4, and Correlation of the Isotropic Shift to Lead–Oxygen Distance in Natural Minerals. Crystals 2019, 9, 43. https://0-doi-org.brum.beds.ac.uk/10.3390/cryst9010043

AMA Style

Zeman OEO, Steinadler J, Hochleitner R, Bräuniger T. Determination of the Full 207Pb Chemical Shift Tensor of Anglesite, PbSO4, and Correlation of the Isotropic Shift to Lead–Oxygen Distance in Natural Minerals. Crystals. 2019; 9(1):43. https://0-doi-org.brum.beds.ac.uk/10.3390/cryst9010043

Chicago/Turabian Style

Zeman, Otto E. O., Jennifer Steinadler, Rupert Hochleitner, and Thomas Bräuniger. 2019. "Determination of the Full 207Pb Chemical Shift Tensor of Anglesite, PbSO4, and Correlation of the Isotropic Shift to Lead–Oxygen Distance in Natural Minerals" Crystals 9, no. 1: 43. https://0-doi-org.brum.beds.ac.uk/10.3390/cryst9010043

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop