Yoshiteru Itagaki*a, Ramakant M. Kadama, Anders Lunda, Einar Sagstuenb and Janina Goslarc
aDepartment of Physics and Measurement Technology, Linköping University, S581 83 Linköping, Sweden. E-mail: yoita@ifmoliu.se
bDepartment of Physics, University of Oslo, N-0316 Oslo, Norway
cInstitute of Molecular Physics, Polish Academy of Sciences, PL 60-179, Poznan, Poland
First published on UnassignedUnassigned22nd December 1999
X-irradiated Li(N2H5)SO4 single crystals were investigated using ESR and ENDOR spectroscopy at several temperatures. The hydrazine radical cation N2H4•+ was selectively produced by irradiation at room temperature. From the analysis of the orientation dependent ENDOR spectra, the 1H-hfc tensors of the cation radical were precisely obtained and the radical structure was supported by theoretical calculations. It is suggested that the cation radical has a planar π* structure D2h (2B2g) in the crystal down to 230 K. By using the evaluated 1H-hf tensor the powder ESR line shape was successfully simulated. Concomitant with the radical formation, the N–N bond of N2H4•+ is suggested to reorient so as to optimize hydrogen bond interactions. 1H-ENDOR line splitting for the N2H4•+ radical was observed at temperatures below 230 K. Apparently this splitting is due to a reversible structural change where one of the NH2 moieties in N2H4•+ becomes slightly bent out of the molecular plane, whereas the other one remains planar. This deformation evidently arises from interactions between the cation radical, adjacent H2N–NH3+ molecules and the SO4–LiO4 framework. Interacting N2H5+···N2H4•+···N2H5+ molecules along the c-axis are proposed to explain the deformation mechanism.
![]() | ||
Fig. 1 Illustration of the arrangement of LiHzS molecules in the crystal highlighting the N–H···N chain of hydrogen bonds believed to be important for the proton conductivity of these crystals. |
Significant reduction of the conductivity was reported by Hilczer et al.14 in a crystal with damage induced by ionizing irradiation. Morawski et al.15 carried out ESR and ESE studies on irradiated LiHzS single crystals to deduce the nature of the radiation defects and their dynamic behavior. These authors suggested that the decrease in conductivity was a result of γ-irradiation producing N˙H–NH3+ radicals and consequently that the hydrogen-bond between adjacent N2H5+ molecules was broken.15 These studies further confirmed that tunneling of protons along the N–H···N chain activated by the –NH2 group rotation is the preferred mechanism for the conductivity in these solids.15–17
Generally, ESR (electron spin resonance) studies of single crystals containing paramagnetic centers give very useful information about anisotropic parameters such as g and hfc (hyperfine coupling) tensors, and the free radical structure and orientation in the crystalline lattice is amenable to be determined from those parameters. A paramagnetic center interacts with surrounding molecules or nuclei, and due to these interactions, the radical species in a crystal may have a structure slightly different from that of the energetically most stable one. This difference, however, is usually indistinguishable by ESR because of line broadening. The ENDOR (electron nuclear double resonance) technique, however, is a quite powerful technique to resolve hf couplings hidden within the ESR line width, and a more accurate evaluation of hfc tensors is possible. Sagstuen et al.18 observed 1H-ENDOR spectra of irradiated N2H5HC2O4 (abbreviated as HzHO hereafter) and successfully obtained 1H-hfc tensors of the N2H4•+ radical species. These authors also detected a small degree of asymmetry of the species.
So far no ENDOR investigation on the structure of the radical produced in LiHzS has been carried out. It is desirable to obtain accurate hfc tensors and information about orientation and small structural changes of the radical in LiHzS for the purpose of understanding the mechanisms governing proton transfer processes in this system.
The crystals were irradiated at room temperature for 1 h using X-ray tube with a W anode operated at 70 kV and 20 mA (approximate dose rate 250 Gy min−1). The irradiated crystals were mounted to the goniometer head of a Weissenberg X-ray diffraction camera, and the crystallographic axes were aligned to rotation axis of the camera within 1.0° using oscillation diagrams. The crystals were transferred to quartz rod for ESR and ENDOR measurements. As described below, some misalignment took place during the transfer process. This was corrected for during data analysis.
X-band ESR and ENDOR spectra were recorded using a Bruker ELEXSYS E580 spectrometer in CW-mode equipped with an EN 801 ENDOR cavity. The ENDOR spectra were obtained with a Bruker E250 ENDOR accessory and a ENI RF 500 amplifier. The radiofrequency (rf) modulation depth for ENDOR measurements was 190 kHz. The sample temperature was regulated by a Bruker VT4111 temperature controller using a cold nitrogen gas flow.
The ENDOR data were analyzed to obtain experimental hfc tensors by using the program ‘‘Magres ’’ written by W. H. Nelson.22
The theoretical calculations were performed using the Gaussian 94 program.23
![]() | ||
Fig. 2 X-band ESR spectra of an X-irradiated Li(N2H5)SO4 single crystal, recorded at 240 K. The magnetic field is parallel to a, b and c axes of the crystal, respectively. Dashed curves are calculated spectra based on the proton hfc parameter presented in Table 1 together with the reported 14N hfc.19 |
![]() | ||
Fig. 3 ENDOR spectra corresponding to the ESR in Fig. 2 recorded at 240 K. The external magnetic field was locked to the central field of the ESR (ca. 340 mT). The lines marked by * is the 2nd harmonic of a real line due to the non-linearity of the rf power amplifier. |
All sets of angular dependent ENDOR spectra were recorded at 240 K by rotating the irradiated LiHzS crystal about the three crystallographic axes. Fig. 4 shows a plot of the ENDOR frequencies of the A and B lines observed at the dif ferent orientations. The angular variations for both A and B lines were accompanied with site splitting. There are four symmetry-related molecular sites in an orthorhombic structure and by rotation about a crystallographic axis, sets of two are always magnetically equivalent and cannot be distinguished by EPR. For rotation about the a and b axes, there were, however, slight misalignments of the crystal on the sample holder, and for this reason further line splitting (all four molecular sites resolved) occurred, and thus coinciding data points at 〈a〉 and 〈b〉 in the figure were not obtained. The ENDOR data were successfully analyzed by taking these misalignments into account and experimental hfc tensors for the A and B resonance lines were evaluated and are reported in Table 1.
![]() | ||
Fig. 4 The angular variation of the high frequency ENDOR resonance lines from the strongly bonded protons in the N2H4•+ cation radical. Solid curves are theoretical fit calculated using the 1H-hfc tensor parameters in Table 1. |
Eigenvectors | |||||
---|---|---|---|---|---|
Principal value/MHz | a | b | c | ||
A | −59.64 | 0.406 | 0.647 | −0.645 | This work |
−37.34 | 0.234 | −0.756 | −0.611 | 240 K | |
+2.28 | 0.883 | −0.097 | 0.459 | ||
B | −0.963 | −0.255 | −0.087 | ||
0.258 | −0.779 | −0.572 | |||
−0.078 | 0.573 | −0.816 | |||
I | −64.20 | 0.389 | 0.642 | −0.661 | This work |
−42.13 | 0.241 | −0.763 | −0.600 | 200 K | |
−2.38 | 0.889 | −0.074 | 0.452 | ||
II | −54.51 | 0.391 | 0.666 | −0.636 | |
−33.63 | 0.191 | −0.733 | −0.652 | ||
4.28 | 0.901 | −0.133 | 0.413 | ||
I | −61.11 | 0.315 | −0.497 | 0.809 | Ref. 18 |
−39.26 | 0.935 | 0.013 | −0.355 | 295 K | |
2.01 | 0.166 | 0.868 | 0.468 | ||
II | −59.33 | 0.306 | −0.515 | 0.800 | |
−37.52 | 0.938 | 0.029 | −0.345 | ||
3.09 | 0.160 | 0.857 | 0.490 |
The evaluated principal values are identical for the A and B. Moreover, they are in good agreement with the proton hfc tensors of N2H4•+ reported for the HzHO single crystal,18,19 and the isotropic 1H-hf (31.57 MHz) is very close to that of N2H4•+ detected in aqueous sulfonic acid solution (30.8 MHz).24 Thus it is proposed that the radical formed in LiHzS is the N2H4•+ species.
The principal values are −59.64, −37.34 and +2.28 MHz, and are expected to be along the perpendicular to the N–H bond, vertical to the NHH plane and along the N–H bond direction, respectively. Sagstuen et al.18 reported two different hfc tensors for N2H4•+ in the HzHO single crystal due to two physically inequivalent sets of protons at 295 K. Since only one set of 1H-hfc tensor was evaluated in the present work, this means that all four protons are physically equivalent even at 240 K; it is striking that the principal values obtained for LiHzS crystal are intermediate between two different values in the HzHO which suggests that N2H4•+ is more flexible to average the hfc in the former crystalline lattice than in the latter. The eigenvectors relative to the crystallographic axes of each principal value are also given in Table 1. Both hfc tensors have one common axis (z[cmb.circ]1, z[cmb.circ]2) and two other axes (x[cmb.circ]1, x[cmb.circ]2 and ŷ1, ŷ2) make an angle of 120.0° corresponding to the HNH sp2 bonding angle. This means that N2H4•+ species exhibits a planar structure in the present crystal, whereas a slightly bent structure was suggested in the HzHO crystal. Thus, ENDOR signals A and B are each attributable to two of the four protons in N2H4•+, respectively, which are chemically equivalent but magnetically inequivalent; if the N2H4•+ molecule exhibits a planar D2h structure, the hfc tensors of the two protons in trans sites are identical.
Theoretical calculations for N2H4•+ have been performed by Almlöf et al.25 by using unrestricted Hartree–Fock level of theory (UHF) and showed good agreement with our experimental values. In the present work, the 1H-hf tensor of N2H4•+ cation radical was theoretically evaluated using DFT (density functional theory) methods and is shown in Fig. 5. The single point calculation was performed employing the B3LYP density function with the 6-31++G(d,p) basis set based on an optimized structure evaluated using the MP2 method with the 6-31++G(d,p) basis set, i.e. B3LYP/6-31++G(d,p)//MP2/6-31++G(d,p). The N2H4•+ molecule was optimized to be a D2h planar structure in which the unpaired electron is delocalized in nitrogen 2pπ* orbital (2B2g). The theoretical principal values of hfc tensor based on the optimized structure are also in good agreement with the experimental ones. The principal axes of the tensor were evaluated to be perpendicular to the N–H bond, vertical to the NHH plane and along the bond direction; the theoretical hfc tensor results are included in Fig. 5 together with experimental values. The nitrogen hfc tensor was also calculated and showed good agreement with that previously reported.19
![]() | ||
Fig. 5 Optimized geometry and calculated hfc tensor for the N2H4•+ cation radical (B3LYP/6-31++G(d,p)//MP2/6-31++G(d,p)). Experimental values of 1H and 14N21 couplings are given in parentheses.19 |
Powder ESR spectra of irradiated sulfate compounds are shown in Fig. 6. The ESR line shapes of LiHzS and HzHO are very similar, and both are in good agreement with the simulated spectrum of N2H4•+ using the proton hfc parameter presented in Table 1 together with the reported 14N hfc.19 The ESR pattern of irradiated N2H6SO4 single crystal has been also assigned to N2H4•+,21 although the powder pattern shown in Fig. 6 has a slightly different line shape due to differences in anisotropic parameters. The ESR spectra of the single crystal in Fig. 2 were also simulated by using the same hfc tensors as the powder case: small degree of inconsistencies are due to the nitrogen couplings.
![]() | ||
Fig. 6 ESR spectra of irradiated powder of (a) Li(N2H5)SO4, (b) N2H5HC2O4 and (c) N2H6SO4. Dashed curve is the ESR simulation for the N2H4•+ radical obtained using the 1H-hfc tensors obtained in this work (Table 1) and the 14N-hfc tensor reported in ref. 19. |
As commented in the Introduction, the radical produced in LiHzS by γ-ray irradiation was previously suggested to exhibit the structure H N˙–NH3+,15 although the powder ESR line shape is similar to that in Fig. 6(a). In the present study, the ENDOR investigation using single crystal and powder ESR line simulations clearly shows that it is the N2H4•+ species which is predominantly produced.
![]() | ||
Fig. 7 Changes in the ENDOR line shape with decreasing temperature, recorded with the magnetic field parrallel with the c-axis. |
![]() | ||
Fig. 8 A plot of the temperature dependent ENDOR frequencies of the strongly bonded protons in the N2H4•+ cation radical. (Upper plot is for the A resonance line, and the lower plot for the B resonance line.) |
The ENDOR lines of the weakly bonded protons did not change to any large extent. The complexity of those signals at 110 K is mainly due to the appearance of intense nitrogen resonance lines accompanied by quadrupole splitting. To further investigate the splitting of the A and B lines, the angular dependent ENDOR spectra were measured at 200 K, where the split lines are well resolved. Thus, two different hfc tensors were obtained for each line (Table 1). The isotropic hfc constants at 200 K are −36.24 and −27.95 MHz and Δaiso, referred to the values at 240 K, are thus −4.67 for I and+3.62 MHz for II (Table 1). Clearly at 200 K two kinds of pairs of physically equivalent protons exist in the radical which most probably is due to asymmetric deformation by lowering the temperature. Comparing the evaluated eigenvectors of the 1H-hfc tensors at 200 K with the corresponding eigenvectors at 240 K, the angles of deviation are 1.51, 1.08 and 1.55° for I (the numerically larger hf), and 0.93, 4.28 and 3.58° for II (the numerically smaller hf). It is obvious that the deformation is larger for the smaller hf II and the observations suggest that the –NH2 moiety has become bent out of the molecular plane. The angles of deviation 1.51° for I and 0.93° for II are in a direction perpendicular to the N–H bond. Apparently, this direction is less affected by the deformation process than the other two principal directions. This deformation will be further discussed in conjunction with theoretical calculations in the next section.
The N–N bond of the N2H5+ fragment in the undamaged LiHzS crystal is directed along (0.581, 0.813, 0.042), whereas the bond of the radical is along (0.804, 0.476, −0.357). Those vectors make angle of 33.0°: the N–N bond direction of the N2H4•+ has changed also in the LiHzS crystal. This reorientation results in the 2pz orbital of N2H4•+ (perpendicular to the molecular plane) becoming oriented to lone-pair orbital (LPO) direction: the angles between the direction of the 2pz orbital of the N2H4•+ radical and that of the N–Hα and N–Hβ bonds of the neighboring N2H5+ molecule (NHαHβ–NH3+) are 108.3 and 111.7, respectively. These angles are very close to the sp3 bonding angle ∠HNH. Thus, it is possible that nitrogen 2pz orbital of the N2H4•+ molecule is oriented so as to become parallel with the LPO of the N2H5+ molecule. Furthermore, the N–N direction of the N2H4•+ molecule makes almost same angle with the N–Hα and N–Hβ bond directions of the N2H5+ molecule (54.3 and 59.6°): this means the N2H4•+ molecule is symmetrically located with respect to the –NH2 group of N2H5+ molecule through the LPO. Assuming that the nitrogen atom of N2H4•+ is located 3.0 Å away from N(1) of N2H5+ molecule based on X-ray diffraction data,4,5 the orientation in the N2H5+···N2H4•+···N2H5+ part is proposed as shown in Fig. 9. Considering that the LPO of the N2H5+ molecule originally is partially oriented towards one of the two protons to form an effective N–H···N chain (Fig. 1), the generation of N2H4•+ consequently might contribute to the proton conductivity loss in irradiated LiHzS.14
![]() | ||
Fig. 9 Schematic representations of the interacting N2H5+···N2H4•+···N2H5+ models. Detailed explanations are given in the text. |
Theoretical calculations were carried out to attain further understanding of the deformation mechanism. As noted in the previous section, the 1H-hfc of the bending part (II) decreased in absolute value compared to the planar structure, on the other hand the absolute value of the 1H-hfc value of the other part increased. This hfc change can be interpreted as follows. The nitrogen spin orbital of the bending –NH2 part in the N2H4•+ contributes a weak hydrogen bond to the neighboring N2H5+ molecule, and consequently the nitrogen acquires some sp3-like orbital nature. The absolute value of the 1H-hfc consequently decreases by a positive spin contribution from the sp3-like spin orbital. By the interaction between 2pz of the N(2) and the LPO of N(3), spin density on the N(2) is partially shifted to the N(1) and then 1H-hfc of the planar –NH2 part numerically increases (Fig. 9). We propose a simple model shown in Fig. 9 in which the N2H4•+ cation radical is put in between two N2H5+ molecules in the c-axis direction (N2H5+···N2H4•+···N2H5+) and the geometry of the two N2H5+ molecules is fixed based on the molecular orientation in the unirradiated crystal (Fig. 1). The site geometry of the N2H4•+ is obtained as follows. The interacting nitrogen N(2) of the N2H4•+ molecule was located 3.0 Å from N(3) as discussed above. The rotation of the N(1)–N(2) bond about any axes was frozen (assuming those motions are restricted by the hydrogen bonds (Fig. 9). Furthermore, the non-interacting –N(1)H2 group was kept planar. Next, the bending angle θ was evaluated by optimizing the remaining variable geometrical parameters and the corresponding 1H-hfc tensors were calculated using the B3LYP/6-31++G(d,p)//B3LYP/6-31G(d,p) density functional method. The optimized geometry exhibited a bent angle θ of +22.1° (positive direction of the c-axis). The lower N2H5+ molecule (H2N(3)–NH3) has a large contribution to the structural deformation compared to the upper molecule through the weak bonding between 2pz orbital of the N(2) and the LPO. Consequently, the N2H4•+ and the lower N2H5+ form a dimer like structure combined through the LPO and 2pz orbital of the N(2) and this is reflected in the positive bent angle. Thus, the hydrogen bonding between the N2H4•+ cation radical and the two N2H5+ molecules leads the cation radical into a bent structure. Probably the theoretically evaluated bending angle is too large, because the calculation of isotropic 1H-hfc for the bending –N(2)H2 group resulted in 1.25 and 4.17 MHz and these are very different from the experimental value of −27.95 MHz. The bending is likely to be somewhat restricted by hydrogen bonding with oxygen atoms in the molecular framework: the hydrogen bonding could also be a reason why the –N(1)H2 part remains planar. To obtain reasonable 1H-hfc's we fixed the bending angle to the experimental value of +2.5° and optimized the remaining variable parameters. As a result, the geometry obtained exhibits reasonable theoretical isotropic hfcs −40.70 and −40.85 MHz (−40.78 MHz in average) for the bending –N(2)H2 group and −45.60 and −45.58 MHz (−45.59 MHz in average) for the planar –N(1)H2 group. The difference from experimental 1H-isotropic hfc (Δaiso) between planar (240 K) and bending N2H4•+ (200 K) was −4.67 (planar part) and +3.62 (bending part). Whereas the Δaiso between the optimized structure of the isolated N2H4•+ (−43.86 MHz) (Fig. 5) and the bending N2H4•+ in the N2H5+–N2H4•+–N2H5+ model (Fig. 9) are −1.73 (planar part) and +3.09 MHz (bending part), and these are close to the experimental values. We conclude that the bent structure of the N2H4•+ cation radical is brought about by the weak bonding with the LPO of one of neighboring N2H5+ molecules, and that the degree of bending is restricted by the surrounding oxygen atoms. At temperatures above 230 K, molecular vibration of the N2H4•+ radical and/or the N2H5+ molecules to some extent can break these interactions, and for this reason the N2H4•+ radical exhibits a planar structure is above 230 K.
This journal is © the Owner Societies 2000 |