Next Article in Journal
Role of the Number of Adsorption Sites and Adsorption Dynamics of Diffusing Particles in a Confined Liquid with Langmuir Kinetics
Previous Article in Journal
Effect of Air Annealing on the Structural, Textural, Magnetic, Thermal and Luminescence Properties of Cerium Fluoride Nanoparticles
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The True Nature of the Energy Calibration for Nuclear Resonant Vibrational Spectroscopy: A Time-Based Conversion

1
SETI Institute, Mountain View, CA 94043, USA
2
Research and Utilization Division, SPring-8/JASRI, 1-1-1 Kouto, Sayo, Hyogo 679-5198, Japan
3
School of Computer Science, Georgia Institute of Technology, Atlanta, GA 30332, USA
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Physchem 2022, 2(4), 369-388; https://doi.org/10.3390/physchem2040027
Submission received: 30 July 2022 / Revised: 12 November 2022 / Accepted: 21 November 2022 / Published: 28 November 2022

Abstract

:
Nuclear resonant vibrational spectroscopy (NRVS) is an excellent synchrotron-based vibrational spectroscopy. Its isotope specificity and other advantages are particularly good to study, for example, iron center(s) inside complicated molecules such as enzymes. In order to investigate some small energy shifts, the energy scale variation from scan to scan must be corrected via an in-situ measurement or with other internal reference peak(s) inside the spectra to be calibrated. On the other hand, the energy re-distribution within each scan also needs attention for a sectional scan which has a different scanning time per point in different sections and is often used to measure weak NRVS signals. In this publication, we: (1) evaluated the point-to-point energy re-distribution within each NRVS scan or within an averaged scan with a time-scaled (not energy-scaled) function; (2) discussed the errorbar contributed from the improper “distribution” of ΔEi or the averaged ΔE within one scan (Eerr1) vs. that due to the different ΔEi from different scans (Eerr2). It is well illustrated that the former (Eerr1) is as important as, or sometimes even more important than, the latter (Eerr2); and (3) provided a procedure to re-calibrate the published NRVS-derived PVDOS spectra in case of need. This article establishes the concept that, at least for sectional NRVS scans, the energy positions should be corrected according to the time scanned rather than be scaled with a universal constant, as in a conventional calibration procedure.

1. Introduction

Calibrating energies is a critical matter for all spectroscopies, including nuclear resonant vibrational spectroscopy (NRVS), which is synchrotron radiation (SR)-based modern vibrational spectroscopy with a lot of distinguished advantages [1,2,3,4,5,6,7,8,9]. In this publication, we provided a brief review on the energy calibration issues for NRVS, discussed the point-by-point energy calibration within one NRVS scan with different data acquisition time per point in different sections, and established the concept that (at least for the sectional scans) the energies should be calibrated according to the time scanned rather than scaled with a universal constant.

1.1. Nuclear Resonant Vibrational Spectroscopy

NRVS is a nuclear resonant inelastic scattering spectroscopy that measures the vibrational modes (or in other words the created or annihilated phonons) associated with the nuclear resonant transition [1,2,3,4,5,6,7,8,9]. As illustrated in Figure 1a, while an extremely monochromic incident X-ray beam (e.g., in case of 57Fe, E1 = 14.41425 keV or 14.4 keV for a simpler description hereafter) scans through the interested nuclear and vibrational transitions, the probed scattering energy is precisely “defined” by the nuclear back radiation at E2 = hν1 which should have the same linewidth as that for the incoming X-ray beam. Meanwhile, both the nuclear fluorescence at hν1 and the converted K-shell electron fluorescence at hν2 are collected as the raw counts. These pure nuclear events (lifetime = nanoseconds) can be distinguished and extracted from the often huge (e.g., 20 million cts/s) electronic scattering background (lifetime = femtoseconds) in the time domain [8,10]. Such a time domain-filtered nuclear signal vs. the vibrational energy (Evib = E1 − E2) forms a raw NRVS spectrum. NRVS becomes available due to the advancement in modern synchrotron radiation (SR) rings, which can provide high intensity/low emittance incoming beams with appropriate time structure(s), advanced high-resolution monochromators (HRM) which reduce the energy bandwidth to 1 meV or narrower (e.g., 0.8 meV at SPring-8 BL35XU, BL19LXU), and fast detectors which distinguish the nuclear events from the huge electronic scattering background. In turn, it presents some distinct advantages in comparison with traditional vibrational spectroscopies such as FTIR [8], resonant Raman spectroscopy [8], and laser-induced fluorescence spectroscopy [11,12,13,14].
Most extraordinarily, NRVS is isotope-specific (site-specific) and is, therefore, an excellent tool to study, for example, Fe-S/P/Cl/O, Fe-CO/CN/NO, and Fe-H/D vibrations [15,16,17,18], etc., inside complicated systems in physics, earth sciences, materials sciences, coordination chemistry, and bioinorganic chemistry [10,16,17,18,19,20,21,22,23,24,25,26]. In particular, when a specific site can be selectively labeled with 57Fe while all the other sites are not labeled [8,27,28], NRVS becomes a pinpoint tool to target the site interested in complicated systems. For example, [2Fe]H site inside several [FeFe] hydrogenases can be selectively enriched with 57Fe and all other iron sites are left unenriched [27,28], making all the NRVS signals originate from the [2Fe]H site. With the same principle, the Fe(II) vs. Fe(III) sites inside a Prussian blue complex or analogs can be alternatively enriched with 57Fe and probed with NRVS [8,29]. Sometimes different elements can be enriched with different isotopes and measured with different NRVS, e.g., (Et4N)3[57Fe4125Te4(SPh)4] can be probed with both 57Fe and 125Te NRVS [30]. On the other hand, NRVS also has coverage for almost all the 57Fe-related vibrations except for a few of them whose 57Fe atoms do not move during the vibrations [16]. In comparison, either IR spectroscopy or Raman spectroscopy has a “restrictive” selection rule.
Due to the nature of nuclear events, which have much longer decaying times (in the order of ns) than electronic scattering (in the order of fs), and due to the fact that nuclear backscattering has a narrow transition bandwidth, NRVS does not use a diffraction spectrometer and thus has a much better photon in and photon out efficiency [6,8,10,20,22]. Due to the same nature, NRVS also has an almost zero background (after filtering out the electronic scattering background), which leads to the observation of some extremely weak signals, such as the 0.1 cts/s Ni−H−Fe [10,22] or X−Fe−H [27,28,31] features inside various hydrogenase samples. NRVS spectra can also start from a real 0 cm−1 while IR or Raman measurement, on the other hand, often starts from 80–100 cm−1.
Due to the simplicity of the NRVS intensity rule (proportional to the displacement of the 57Fe site in a particular vibrational mode), the pure molecule-based/experiment-independent partial vibrational density of state (PVDOS) can be obtained from all raw NRVS spectra. Therefore, NRVS peak positions, as well as its intensities, can be well reproduced via various theoretical calculations, including density function theory (DFT) [27,28]. In contrast, IR or Raman calculation often involves assumptions about molecular properties such as trail dipole moments or polarizabilities in order to obtain an approximate molecular VDOS.
With regard to practical aspects, NRVS has some compelling practical advantages over established methods. For example, it is water transparent in comparison with far IR spectroscopy, and thus well suited for investigating biological samples in their natural aqueous environment; it is free of fluorescence problems in comparison with resonance Raman spectroscopy, and thus suitable for evaluating photosensitive states; it distinguishes well among O, N, and C in comparison with extended X-ray absorption fine structure (EXAFS) or crystallography [32,33]. Due to these advantages, its user community, as well as the research results, grow exponentially [10,16,17,18,19,20,21,22,23,24,25,27,28,31,34]. For example, in the last decade, it has become the third modern X-ray spectroscopy to become popular among biochemical researchers, following crystallography and EXAFS [8,33].
The simplest NRVS scan is an even time scan which measures each data point with the same amount of data acquisition time (also called scanning time), for example five seconds per point (s/p), we call it an even time scan or even scan. For investigating biological samples or other dilute samples, an extremely long scanning time per point must be used in order to obtain spectra with satisfactory statistics. For example, an 89 KDa NiFe hydrogenase molecule has one Fe inside the interested NiFe center, leading to a wt./wt. concentration of ~0.06% Fe in the NiFe center. This does not include the H2O amount used as the solvent. To probe such a low concentration of Fe, the experimentalists often start to use three seconds per point (s/p) for the low energy region (e.g., from −240 cm−1 meV to +400 cm−1) to measure Fe−S-related vibrations and 10 s/p for the higher energy region (e.g., +400 cm−1~+650 cm−1) to measure weak Fe-CO/CN and extremely weak Ni−H−Fe vibrations [22,23]. We call it a sectional scan. When the NRVS features in the Fe−S and Fe−CO/CN regions are well resolved, the experimentalists start to focus on the extremely weak Ni−H−Fe features with 30 s/p scanning time in the corresponding region (e.g., +650 cm−1~+800 cm−1) while leaving 1 s/p or sometimes even 0.5 s/p for the lower energy region (−240 cm−1~+650 cm−1). In some extreme cases, 30 s/p will be used for the X−Fe−H features while 1 s/p is used for the nuclear resonant peak region, and the region between the two will be skipped—we call this a jump scan.

1.2. Energy Calibration in NRVS

Energy calibration is a procedure that uses a correlation function to map the observed raw energies (Eobs) onto the calibrated real energies (Ereal). All spectroscopies need energy calibration, those using SR as light sources are no exceptions [8,35,36,37,38]. For NRVS, while the high-energy resolution beam from an HRM provides the benefit to resolve detailed vibrational spectral features, the extremely high resolving power (>107), which is defined as the exit energy (e.g., 14.4 keV) over its narrow bandwidth (e.g., 0.8 meV), also makes it harder to track and maintain accurate/precise exit energies from an HRM. In other words, it demands more careful energy calibration. For example, in an extreme case, a 0.1 meV energy shift between two spectra can be observed as long as the following two conditions are met: (1) the two spectra in comparison must both have good enough signal-to-noise ratios (S/N) to define and compare the peak centroids with a 0.1 meV difference, and (2) the energy calibration must have a better than 0.1 meV precision/accuracy. While S/N can be improved with more data acquisition time, a good enough energy calibration becomes the key to observing small energy shifts. Although a high energy resolution in the incident beam is also preferable to obtain sharp spectral peaks, this is not a requirement for observing a small energy shift.
According to classical concepts [39], the instrumental scale (ruler) can be calibrated while additional dispersions (instabilities due to environmental differences) are attributed to random/un-calibratable source(s). However, not all the dispersions are imprecisions and at least some energy instabilities are also trackable and “calibratable” [40]. Therefore, there should be two categories of energy calibration for a monochromator, spectrometers, or other instruments. The first one is often performed by the instrument designers or beamline scientists to calibrate a particular device as finely as possible and as repeatably as possible — the aim is to evaluate the design and/or establish the scale for the device. We call this the absolute calibration. To perform an as-fine-as-possible calibration, a high order or more complicated function form is often used to calibrate Eobs to Ereal [41]. For being as repeatable as possible, the calibrations should be performed under a standardized condition, e.g., at 25 K, etc., when the measured correlation between Eobs and Ereal will be stable and repeatable. Although important, this method is not the focus of this publication.
The second is the one that spectroscopists use to calibrate their measured data under whatever the experimental conditions and/or at whatever moment of the measurements—we call this a practical calibration. Such a calibration often omits the higher order corrections and calibrates the energy axis with a linear correlation:
Ereal = (αEobs + E*)
where Ereal, Eobs, E* stand for the calibrated real energy, the uncalibrated observed energy, the correction for the zero energy position, respectively, and α is the energy axis scaling factor which needs to measure a standard calibration spectrum [5,6,20] to calculate. All the calibrations discussed in this publication are practical calibrations.
Figure 1b summarizes examples of a series of the uncalibrated spectra for [FeCl4] measured at four nuclear resonant scattering (NRS) beamlines around the world. All these spectra have different energy scales [6,8]—using α = 0.999 for APS 03ID, 0.986 for ESRF ID18, 0.965 for SPring-8 BL09XU, and 0.920 for Petra-III P01, all the uncalibrated spectra will be aligned to the published data for its IR spectrum [16]—no figures. If each beamline can have a fixed α value as above, the corresponding α value is referred to as the energy calibration scale for that particular beamline or that particular HRM. Unfortunately, as mentioned in an earlier publication [40], and as will be detailed in the following, this is not the case: the measured α values are different from beamtime to beamtime and sometimes even from different calibration measurements within one beamtime. For example, the measured α values for SPring-8 BL09XU were once between 0.938 and 0.976 [5]. Later, the dispersion span becomes narrower: 0.952–0.966. Still, there is >0.01 in its variation span (or >1% vs. the average α value). BL19LXU is a non-dedicated NRS beamline and its energy scales range from 0.918 to 0.937, dispersing over 2%. For most other NRS beamlines, such variations also exist.
Since most of the Fe−X stretching and bending vibrations have 2–4% in their isotopic shifts e.g., X = 13C/12C or 36S/32S [5] and even a larger amount in their redox shifts [5,17,23], the current once per beamtime calibration or quick switch calibration [5] still serves the general purposes for NRVS measurements.
It seems no surprise that the measured α values can be different at different moments because practical calibrations were performed under various experimental conditions (e.g., a minor difference in the temperature surrounding the monochromator’s crystals). A realistic practice to better track the energy variation is to make more frequent calibrations during NRVS beamtimes. In addition to the conventional energy calibrations which can take six hours in total time [5,6,23], a quick switch calibration procedure was also explored and published earlier [5], in which the incident X-ray beam can be altered between to pass through and measure the main NRVS sample inside the cryostat and to pass right over the main sample’s top surface and measure the calibration sample at an room temperature (RT) stage behind the cryostat [5]. This lets calibration measurements be done without moving the main NRVS samples out of the cryostat base, saving the time for sample changing and cryostat temperature cycling. Limited progress has been reported using this calibration: ≤0.4% difference was found between two adjacent calibrations when it is performed in every six (6) hours—using tons of beamtime. Nevertheless, there is no guarantee that the scans between the two calibration measurements will have a fixed α value or a predictable dispersion span in its possible values.

1.3. Using In-Situ Energy Calibrations

The α values measured with practical energy calibrations should include the real instrument-related scale (αins, which could be equal or at least close to the value from an absolute calibration) and the time-dependent portion (αvar): α = αins·αvar. The conventional NRVS calibration procedures are equivalent to treating the energy instabilities from time to time (αvar) as un-calibratable imprecision and using whatever the measured or averaged α value to represent the real energy calibration scale for each beamtime. This issue must be resolved because it limits the accuracy of the energy calibration and the overall benefit of NRVS. For measuring and distinguishing small energy shift(s), the time-dependent portion (αvar) must be tracked, and an in-situ energy calibration procedure becomes necessary. For examples, the 2.4 cm−1 isotopic shift (383.1 → 385.5 cm−1) was observed for [57Fe4S4Cl4]= → [XFe4S4Cl4]= [5] where XFe stands for the mixture of 57Fe or 54Fe. In a different case, a 2–4 cm−1 peak shift was found in the X−Fe−H features between the natural (CH) and 13CD substituted [FeFe] hydrogenase [31], etc. However, in the latter, when the main NRVS sample has a low concentration in the interested iron site [27,28,31], the experimentalists must leave most of the incident X-ray beam for the main NRVS sample and thus have not much to share for the calibration sample, which often leads to quite noisy calibration spectra [40]. This shows that, although useful in some cases, a direct in-situ measurement is still not an excellent solution to obtain a better or more reliable calibration in most cases.
Another option is to use an environment-related parameter, e.g., temperatures, to track the real exit energies in situ. For example, the four-bounce HRM at APS 03ID uses the angle positions as well as the temperatures at its four crystals to calculate the real beam energies. As a result, the exit beam energies at APS 03ID almost always have an excellent scaling factor around α = 1, e.g., 0.999 per the authors’ experience as shown in Figure 1b (red curve). However, most other NRS beamlines do not have a temperature monitoring system. The experimentalists must then look for other parameters to represent the temperature variations. According to a recent publication [40], the zero energy drift (ΔEi, i = scan number index) between adjacent scans becomes the candidate to present and calibrate the variation in the energy scale (αi) although the individual E0 peak positions can just be “used” to correct for the E* value in the conventional calibration procedure [6,8,20]. Instead, the real energies Ereal can be calculated from the observed energies Eabs via the following linear function vs. the accumulated scanning time (Σtk) [40]:
Ereal = [Eobs − (Σtk/Ttot)⋅(ΔEi)]⋅α0
where Eobs, Ereal, and ΔEi are defined as above, while the rest variables/parameters are defined as: (1) k = the data point index inside one particular scan while i = the scan number index; (2) tk = the data acquisition time at each data point, Σtk = the accumulated time scanned from the beginning of the scan to the particular point (1 → k); and (3) Ttot = the total scanning time for the whole scan, Σtk/Ttot is the ratio of the accumulated time scanned (at the point k) to the total time for the whole scan [40]. The Formula (2) can be understood as a stepwise process: (1) Eobs is corrected for the energy variation from scan to scan using E′ = [Eobs − (Σtk/Ttot)⋅(ΔEi)] to create an intermediate energy axis (E′); (2) E′ is further converted to the calibrated real energies Ereal via an additional universal scaling factor (α0) which should be a constant for each beamline or each HRM [40] (or similar to the αins as mentioned above). Such a stepwise calibration leads to a much better energy correction than the ones using (1) [40].
In addition to tracking the αi variation from scan to scan via ΔEi, Formula (2) also includes the nature that Eobs should be mapped onto Ereal according to the accumulated time scanned at each point. In other words, the energy drift amount per scan (ΔEi) needs to be distributed to the individual data point (k) according to the ratio of the accumulated time scanned vs. the total time for one scan (Σtk/Ttot) rather than scaled with a universal factor [using the Equation (1)]. Although this formula was first proposed in an earlier publication [40], the nature and the application of its time-scaled energy distribution were not elaborated. The focus of that publication [40] was to track and correct the energy variations between different scans rather than within one scan. In this publication, instead, we will focus on: (1) evaluating the distributed energies (Eobs → Ereal) with the time-scaled function within one NRVS scan or an average NRVS scan while assuming one representing ΔE value for all NRVS data; (2): analyzing the errorbar contribution due to the improper “distribution” of ΔE within one scan vs. that due to the different ΔEi from different scans. It will be well illustrated later that the former part is as important as or sometimes even more important than the latter part; and (3) establishing a procedure to re-calibrate previous NRVS spectra in case of need. In addition, the case of jump scans where an “unimportant” energy region is skipped will also be evaluated.
Finally, this article aimed to establish the concept that energy positions in one sectional NRVS scan should be tracked and corrected according to the time scanned [as in Equation (2)] rather than be scaled with a universal constant [as in Equation (1)] even though all the NRVS scans have the same zero energy drift value ΔE.

2. Experimental Aspects

The NRVS data cited or tested in this study were previously measured at SPring-8 BL09XU or BL19LXU. The science for these NRVS data has been published elsewhere in various scientific journals [10,17,18,22,23,25,27,28]. The raw NRVS data were converted into PVDOS via a PHOENIX software package [42,43] or a PHOENIX inclusive web tool at http://spectra.tools [19,20], which was also detailed elsewhere. As we focused on the energy calibration process, the details about the samples/NRVS measurements/PVDOS conversion are omitted here.
The calibrated value for the zero energy position should always be zero [Ereal(0) = 0]. However, the observed zero energy position [Eobs(0) or E0 hereafter] often has a continuous drift in the subsequent NRVS scans (see Figure S1 in the supporting information). The amount of the E0 drift per scan ΔEi is defined as ΔEi = E0(i) − E0(i−1), where E0(i) is the E0 value for the current scan and E0(i−1) is the E0 value for the previous scan. For the first order approximation, the energy drift amount for the first scan (ΔE1) is assumed to be the same as that for the second scan (ΔE2). For previous processed PVDOS spectra, original ΔEi values can be further averaged to one ΔE value to represent the energy drift for the converted PVDOS.
Using ΔEi to re-distribute energy drift per scan point by point inside one NRVS scan or using an averaged ΔE to re-calibrate the previous “calibrated” PVDOS is the central topic of this publication and will be discussed in detail in the following sections.

3. Results and Discussions

3.1. Understanding Time-Scaled Calibrations

To understand the mechanism for Formula (2) [40], let us start from the basics—the law of crystal diffraction:
Eobs = hc/(d0·sinθobs)
Ereal = hc/[d(T)·sinθobs]
where T and d(T) are the temperature and the corresponding atomic space of the diffraction crystals inside the HRM at the moment of measurement; d0 is the atomic space at a standard condition (e.g., at 25 K); Eobs is the observed energy corresponding to the real energy Ereal; h and c are constants. The comparison of (3) vs. (4) leads to:
Eobs = [d(T)/d0]·Ereal ∝ d(T)
Therefore, corresponding to one real energy (e.g., at the nuclear resonant peak Ereal = 14.414 keV), the observed energy Eobs is proportional to d(T), or Eobs ∝ d(T). In short, the energy observed at the HRM (Eobs) is a function of the monochromator crystals’ temperature (T), which is mentioned in the introduction and illustrated here as Equation (5). Although Equation (2) is developed from a pure experimental observation [40], a simplified model of single crystal diffraction and their Equations (3)–(5) provide the readers a theoretical outline of the meaning behind it.
From a qualitative perspective, the trend can be understood in this sequence: once the beam shines on the HRM → the heat load at its crystals increases↑ → crystals’ temperature (T) increases ↑ → the crystals’ atomic space [d(T)]↑ → Eobs↑—the longer the time scanned (Σtk), the higher the crystal temperature (T), the larger the energy drift amount (ΔEi). From a more quantitative perspective, Equation (5) tells us Eobs ∝ d while the thermal expansion properties for silicon around RT indicates d ∝ T [44], therefore Eobs, d, and T have a linear relationship with one another (or more particularly Eobs ∝ T) within a small variation of temperature T. Since the temperature changing rate is proportional to the temperature difference [e.g., dT/dt ∝ (TM − T)] [45,46], its integration shows T (thus d or Eobs) has an increasing form of the exponential decay function with time t as below [40]:
T = TM − (TM − T0)·e−kt
d = dM − (dM − d0)·e−kt
Eobs = EM − (EM − E00)·e−kt
where t is the beam-on time from the beginning of the measurement, T is the crystal temperature, d is the atomic spacing of the crystal [=d(T) in (4) and (5)], Eobs is the observed energy at the HRM; T0, d0, E00 are their initial values at RT; and TM, dM, EM are their equilibrium limits when the beam is on for a long time and the k here is the Boltzmann constant. The function (8) shows that the observed energy position is the exponential function of the beam-on time. For a first order approach, these exponential functions [especially (8)] can be approximated with a linear function vs. time (t) when we approximate the exponential function with the first two terms of its Taylor series expansion. This explains the nature that the energy drift amount within one scan (ΔEi) should be “re-distributed” according to a (linear) function of the time scanned rather than a (linear) function of the energy position.
In a practical aspect, Formula (2) can be understood as a stepwise process that includes: (1) Eobs → E′ which re-distributes ΔEi within each scan and (2) E′ → Ereal via an additional universal scale with an instrumental-related α0. Since the measured zero energy point [such as E* in (1)] can always be calibrated in the final data analysis in PHOENIX or so [5,8], we can make any initial alignment between the observed and calibrated energies—it is not necessary to align them to their nuclear resonant peak, as in Formula (1). The NRVS measured at SPring-8 scans in the direction from the higher energy end (let’s call it E2′) to the lower energy end (E1′), it is easier to align the observed E2obs and the intermediate energies E2′, i.e., set E2obs = E2′, as illustrated in Figure 2a vs. Figure 2b. To illustrate the concept, we assumed an exaggerated larges ΔE which is half of the size of (E2′ − E1′)—that is (Eobs2 − Eobs1) = (½)·(E2′ − E1′) in Figure 8. When the scan reaches the ending energy at the left end (the setpoint at E1′), its observed energy position E1obs should have drifted +ΔEi in comparison with its target (E1′) or, in another word, E1′ should have a −ΔEi correction in comparison with the observed E1obs: E1′ = E1obs − ΔEi, as shown in Figure 2b vs. a. A positive energy drift (ΔEi > 0) means the E1′ is lower in energy position than E1obs, and vice versa. Here, we in fact have assumed that the amount of energy drift when HRM scans from E2′ to E1′ is the same as the E0 position drift (ΔEi) between the current and the previous scans, which is almost true.
When the scan is arriving at a middle point between E2′ and E1′, the energy drift between the observed Eobs and the ΔE corrected E′ accumulates gradually and continuously, starting from 0 to ΔE = (½)·(E2′ − E1′). According to (8), this partial amount should be in general following an exponential decay function vs. the beam-on time (t). Within a short period (e.g., one hour or one scan), however, a linear function can almost approximate the relationship between the energy amount drifted (at point k) vs. the beam-on time t from point 1 to point k (Σtk). Therefore, the Eobs energy position at one particular point k should have drifted a part of the—ΔE (or ΔEi if each NRVS scan is processed individually) that is proportional to the beam-on time which is the Equation (2). When an even scan is performed, the result from the Equation (2) is the same as that from the Equation (9) [40]:
E′ = Eobs − {(E2obs − Eobs)/(E2obs − E1obs)}⋅ΔEi
          = Eobs⋅{1 + ΔEi/(E2obs − E1obs)} − {E2obs⋅ΔEi/(E2obs − E1obs)}
Under this condition, it is reasonable to uniformly scale the Eobs according to the Equation (9), rather than transform it point by point with the Equation (2). Therefore, it is the nature of a sectional scan that makes it necessary to correct the energies point to point according to the accumulated beam-on time (Σtk), which is (2).

3.2. Comparing Two Calibration Procedures

Figure 2 assumes an exaggerated 8-data-point NRVS file that scans from the higher energy end E2′ (right) to the lower energy end E1′ (left), mimicking the real situation for the data measured at SPring-8. Figure 2a illustrates the theoretical energy range in E′ in mind when a scan is planned while (b) illustrates the observed raw energies (Eobs). The data points 1, 2, 3, 4 (red) are measured with 3t s/p, and the points 5, 6, 7, 8 (blue) are measured with t s/p, 1/3 of that for the points 1, 2, 3, 4. This conceptual spectrum has a zero vibrational energy position at point 7 and a few assumed “NRVS peaks” at points 5 (grey), 4 (yellow), 3 (green). We also assume an exaggerated but illustrative energy drift amount of ΔE which is equal to half of the targeted energy span (E2′–E1′). When such an NRVS is scanned, we can first align the starting points for the observed and intermediately calibrated NRVS and let E2obs = E2′. When the NRVS reaches its ending points, its observed energy should have an energy drift of +ΔE in comparison with the targeted energy E1′ (or E1′ = Eobs1 − ΔE). This E1′ should have the same energy value no matter whether we use the calibration procedure (9) [equivalent to (1)] or the procedure (2) to calibrate it, such as shown in Figure 2 vs. d: Figure 2c illustrates the E′ which is calibrated according to a universal energy scale (9) and 2d illustrates the E′ which is calibrated according to the accumulated beam-on time [according to (2)]. The re-distributed E′ for the two calibrations in the middle region between E1′ and E2′ are very different: Figure 2c vs. d.
However, the peak position difference between Figure 2c and d is not the real difference between the energy positions obtained from the two calibration methods. To compare, the zero vibrational energy positions for 2c and 2d should first be aligned with each other, e.g., aligned to 0, shifting 2d to 2d’. Comparing Figure 2c vs. d’, the peak position difference is the real difference between the two energy calibration methods. It is obvious that the peak position difference is relatively small near E0 or near the start of the scan ~E2′. On the other hand, the energy difference for the peaks around the changing point for the data acquisition time, e.g., at points 4 and 5, is much larger than for the other peaks. The NRVS experimentalists often care more about the peaks above the time changing point because that is often the region of interest.
Figure 3 illustrates another imaginary NRVS file but with more realistic scan parameters including an assumed total energy drift of 0.8 meV per scan in a total scanning span of 1000 cm−1. Figure 3a is an example of scanning time (s/p) at each data point: the time is not to the scale but the scanning time per point is noted on the grey horizontal line: 30 s/p is used from 800 to 650 cm−1 while 1 s/p is used for the rest (low energy) regions. Although Figure 3 is from pure conceptual data, such a 30 s/p vs. 1 s/p scanning time has often been used to scan the Ni–H–Fe wagging mode in DvMF [NiFe] hydrogenase [10,22] and the X–Fe–H bending features for several [FeFe] hydrogenases [28,31] in a similar region. Figure 3b shows the energy re-distribution of the assumed 0.8 meV energy drift within one NRVS scan: the red curve represents the re-distribution of this 0.8 meV according to the accumulated time scanned [using the Equation (2)]; the blue curve shows the even spread of the 0.8 meV via a procedure similar to the Equation (9). Although the energies are in general different for the two energy correction methods, the two curves in (b) (blue and red) begin at the same starting point (E2′ = E2obs) and reach the same ending point at E1′ = E1obs − 0.8 meV. The difference between the two calibration methods (blue vs. red) reaches its maximum at 650 cm−1 when the scanning time per second changes from 30 s/p to 1 s/p, which is consistent with the conclusion from Figure 2.
In addition to the analysis with formulas, the two calibration processes can also be understood figuratively as the following. The region of 650–800 cm−1 accounts for (150/1000) = 15% of the whole energy range (from −200 cm−1 to 800 cm−1) but takes 30∗150/(30∗150 + 1∗850) = 84% of the total scanning time. Therefore, an energy-scaled traditional calibration procedure [the Equation (9)] leads to 15% of the −0.8 meV (or −6.45 cm−1) being assigned to this region while a time-scaled calculation procedure [the Equation (2)] leads to 84% of the −0.8 meV (or −6.45 cm−1) being assigned to the same region, leaving −4.6 cm−1 difference between the two calibration methods at 650 cm−1, as shown in Figure 3b. However, this −4.6 cm−1 is not the real difference between the two calibration methods. These curves must still align their zero vibrational energies with each other which lead to the curves in Figure 3c. After the zero energies for the two calibrations are re-aligned to 0 cm−1 as shown in Figure 3c, the maximum difference at 650 cm−1 becomes about 3.7 cm−1 instead, which is the real difference between the two “re-distribution” methods.
Figure 3d illustrates one assumed NRVS spectrum in the region of 625–705 cm−1 with two peaks at 650 and 675 cm−1 respectively. The calibrated spectra using the above two calibration methods are shown as in Figure 3e [blue = the spectrum calibrated with a universal energy scale using the Equation (9) vs. red = the spectrum calibrated with time-based re-distribution using the Equation (2)]. Although both methods (e, blue or red) lead to larger peak position shifts in comparison with the original un-calibrated data (d, black) at the higher energy peak (675 cm−1), it is the peak at 650 cm−1 which has the larger difference (error) between the two different calibration methods (blue vs. red). The value of this error varies according to the scanning details of the sectional scans. It is 0.46 meV (or 3.7 cm−1) here.
As mentioned in the introduction [40], the Equation (2) was initially used to account for the αi variation due to the different ΔEi values from different NRVS scans [40]. Then, the errorbar due to the difference between different scans is estimated as (ΔEi − ΔEcal) = 0.8 − 0.3 = 0.5 meV (=4.03 cm−1), where ΔEi corresponds to the energy drift for different scans while ΔEcal stands for the energy drift under which the calibration scan is measured (which is assumed at 0.3 meV here). In comparison, under an assumed 0.8 meV energy drift per scan and the above scan parameters, the Eerr1 (0.46 meV) is almost the same as Eerr2 (0.5 meV). In case that a ΔEi = 0.6 meV is assumed, the Eerr1 will be 0.35 meV while Eerr2 = 0.6 − 0.3 = 0.3 meV: Eerr1 > Eerr2. Therefore, Eerr1 is a critical portion of the total errorbar that needs to be resolved.

3.3. Re-Calibrating Published PVDOS

According to the above discussion, when the energy drift amount per scan (ΔEi) has about the same level, the errorbar due to the energy variation from scan to scan (Eerr2) can be omitted. However, for sectional scans, the “internal” re-distribution of the averaged energy drift (ΔE) still causes an obvious (if not significant) errorbar (Eerr1). How can we “re-calibrate” the old NRVS to include the time-scaled re-distribution using the Equation (2)?
The first task is to “reverse the original calibration process” using the original energy scale constant α: E*obs = E*real/α, where E*real is the energy that was already “calibrated” with the traditional calibration procedure (1). In this step, we need to use the original value of α from the old analysis record. Most published NRVS spectra were already in the form of PVDOS which was averaged and transferred from a series of raw NRVS spectra [6,42,43]. Such spectra (PVDOS in E*real) often have E*real(0) = 0 while the E*obs can in principle have any value including 0, so we do not need to care about it here. In case of need, it will be automatically shifted in the final analysis with PHOENIX [5,43], just like the constant c in an infinity integration.
The second step is to convert the pseudo “raw” data just obtained (in E*obs) to that in a new intermediate energy scale E′ using the first portion of the Equation (2). For processing individual scans, ΔEi is needed to convert it from E*obs to E′. For processing an averaged PVDOS, a uniform ΔE is needed, which can be obtained by averaging the ΔEi per scan from the original raw data corresponding to the PVDOS.
The last step is to convert the data in E′ to that in Ereal using Ereal = E′·α0 which needs a scaling constant α0. This constant is in principle a beamline-dependent constant and must be obtained via a large amount of energy calibration data using a previously proposed stepwise procedure [40]. At the moment, only the α0 value for SPring-8 BL19LXU has been evaluated as 0.918—a preliminary number [40]. Before the α0 values become available for all other beamlines, they can be estimated by using the Equation (1) [or its alternative form (9)] and (2):
α0 = α/{1 + ΔE/(E2*obs − E1*obs)}
which shows that α0 can be estimated from α (for the original energy calibration), ΔE (which is averaged from the ΔEi for the original raw data), and (E2obs − E1obs) (which is the energy range scanned). A combination of all the above three steps leads to:
Ereal = [E*obs − (Σtk/Ttot)⋅(ΔE)]⋅α0
               = [E*real/α − (Σtk/Ttot)⋅(ΔE)]⋅α/{1 + ΔE/(E2*obs − E1*obs)}
where E*real is the energy axis from the previously “calibrated” NRVS (PVDOS) while Ereal is the newly calibrated energy axis processed with the new time-scaled calibration procedure (2), all other variables are defined as before.
When a published PVDOS needs to be re-calibrated, in addition to the published PVDOS (in E*real), the original NRVS data are also needed to find out the tk and Σtk as well as the ΔEi to calculate an averaged energy drift ΔE. The α0 also must be calculated from α value that was used in the previous calibration according to the Equation (10).

3.4. Examples of Re-Calibrated NRVS

Let us start to look at a few real cases. [FeFe] hydrogenases are complex metalloenzymes, which can lead to hydrogen production in numerous organisms [47], and therefore have potential importance in the post-carbon era economy. The mechanism of H+ transfer, Fe−H or Fe−H2 bonding, and H2 bond activation are all associated with the aminodithiolate (ADT) group inside the [2Fe]H sub-cluster which is the enzyme’s active center. Chemically, scientists have successfully developed artificial maturation [26,27,28] which can replace the original NH at the ADT with an O atom to form an oxodithiolate (ODT) variant or to change the functional group at the surrounding position 169 from –SH to –OH (from Cystine to Serine), in addition to enriching a few selected iron sites with 57Fe while leaving other irons unenriched. Spectroscopically, the site-specific NRVS can extract just the 57Fe-related vibrations and thus is better to distinguish the minor differences among different enzymatic variants. For example, two X−Fe−H wagging and bending vibrational peaks are observed from Chlamydomonas reinhardtii [FeFe] hydrogenase (or Cr HydA1 for short) and its variants [27,28], with one at ~675 cm−1 and another in the region of 725–775 cm−1. The first peaks are almost at the same location for all the variants but the positions for the second peaks are different corresponding to the different variants: ODT variant (green) → wild type enzyme (blue) → C169S variant (grey), as illustrated in Figure 4a. The detailed science and meaningfulness are as published [27,28] and will not be repeated here. In short, understanding the variation of these X−Fe−H peaks and their shifts corresponding to the variations in and around the NH site is relevant to the catalytic mechanism of the hydrogenases.
As these peaks are in the energy region where the 30 s/p data acquisition time is used while the front region is scanned with a 1 s/p data acquisition time, the peak positions can have an obvious error in the traditional energy calibration (1) as discussed above. This can be re-corrected with a time-scaled energy re-distribution process (2) or (11). With ΔE and α which are obtained from the original data set and analysis record, Equation (11) produces a new PVDOS with re-distributed energies (Ereal) shown as a red solid curve in Figure 4b. For a direct comparison, the traditionally calibrated PVDOS for the same sample is also duplicated from Figure 4a to Figure 4b as a dashed blue line. A small shift of ~3 cm−1 is observed [Figure 4b blue vs. red]. Although shifts are not great, it is large enough to affect the comparison of some small energy shifts, such as the 2 cm−1 as mentioned in the introduction for the CH-ADT and its C13D substituted counterpart [31]. We also noticed that the peak at 675 cm−1 has a bit more energy difference between the two calibration methods than the peak at 750 cm−1 does, consistent with the analysis in Figure 2 and Figure 3.
In addition to the X−Fe−H bending features, NRVS for the Fe−H stretching mode is also observed but at the current time only for model complexes. Since the mass ratio of H:Fe is 1:57, most of the motion in the Fe−H stretch mode is at H, not at Fe. This nature leads to a very weak NRVS signal for the Fe−H feature although it is still the best method to observe a Fe−H vibration. For example, several Fe−H and Fe−H2 stretching modes in trans-[57Fe(η2-H2)(H)(dppe)2][BPh4] (H2FeH for a simple description) were published [48] and are re-illustrated as in Figure 5a, again with the traditional calibration process (1) (blue, corresponding to E*real). The vibrational modes associated with the NRVS peaks at 1773 and 1915 cm−1 were assigned to the asymmetric Fe–H stretching mode from the Fe(H2) component and the Fe–H stretching from the “local” FeH bond [48]. The peak at 1956 For cm−1 is the Fe–H stretching mode from a byproduct when H2 accidentally dissociates from the complex [48]. Nevertheless, this peak is the NRVS feature with the highest energy observed to date. When we use the original α value, the averaged energy drift per scan ΔE, and Equation (11) to re-calculate the energy transformation E*real → Ereal, the obtained spectrum is illustrated in Figure 5a (red) while the original spectrum is shown n the same figure as blue. To better illustrate the details, Figure 5b,c show the same spectra in two narrow regions.
In this particular case, 20 s/p was used for 1600–2000 cm−1 while 0.1 s/p was used for the lower energy region. It is clear that the closer the peak position is to the scanning time changing point, the higher the energy shift amount is. We observed that there is a −3.4 cm−1 for the peak at 1773 cm−1; a −2.3 cm−1 for the peak at 1915 cm−1; and a −1.8 cm−1 for the peak at 1756 cm−1.
In Figure 2, the peak position lower than the timing changing position also has obvious differences between the two calibration procedures, such as point 5 in Figure 2c vs. Figure 2d’. In a real NRVS, for example, the Fe−CO peaks in the DvMF NiFe hydrogenase can have obvious differences in their peak positions in the spectrum calibrated with the two different calibration procedures. This is well illustrated in Figure 6. On this side, the lower the peak position, the smaller the energy difference between the two calibration procedures: e.g., −3.0 cm−1 for the Fe−CO peak at 606 cm−1; −2.4 cm−1 for the Fe−CO peak at 548 cm−1. People usually do not care much about the lower energy side because the region(s) are usually scanned in much less time than the focused higher energy region. Nevertheless, in cases where the Fe−CO peaks are to be used as internal calibration marks for tracking the X−Fe−H bending peak positions, the issues discussed here become critical to consider.

3.5. Dealing with Jump Scans

A sectional scan does not necessarily include just two sections; multiple sections or other formats of time structure are also possible. In one example, 30 s/p can be used for the interested region, such as the above X−Fe−H features [27,28], while 1–2 s/p can be used for the nuclear resonant peak region. The region between these two scan regions can be skipped without scanning—we call this type of scan a jump scan. Such a jump scan is often used to measure: (1) the NRVS whose general feature is well observed and only a small portion of interested region needs to be measured further, such as the X−Fe−H features [27,28]; (2) a particular feature in unenriched samples, such as or the Fe−CN for the unenriched Fe(CN)6 related complexes which is the only measurable peak [29,49,50,51].
Figure 7a first exhibits the NRVS for 57Fe enriched (NH4)2Mg(II)[57Fe(II)(CN)6] [49] [29] as a reference. It has a prominent feature around 594 cm−1, which corresponds to Fe(II)(CN)6 core vibrations [29]. Other features are much weaker. For evaluating unenriched (similar) complexes, the weak features are even weaker and are probably out of the detection limit [the 2% 57Fe in natural abundance vs. the 100% 57Fe in the fully enriched samples] [29]. For those cases, measuring a full NRVS scan is usually out of question and probing the prominent feature at 594 cm−1 becomes the only hope. Therefore, a jump scan will be useful. For example, measured NRVS for unenriched KEu(III)Fe(II)(CN)6 is shown as in Figure 7b. It is calibrated with the traditional procedure (1) (blue), which scales energies with a universal factor (blue) and the new calibration procedure (2) (red), which re-distributes the energies according to the time-scanned (Σtk). Figure 7c illustrates the seconds used to scan each data point—the text clearly shows 2 s/p for the nuclear resonant peak and 50 s/p for the Fe-CN region. As each scan took about 1.5 h (relatively long), and the beam is not very stable at the moment of the measurement, the energy drifts (the ΔEi values) are ~1.5 meV per scan (a bit higher than the regular scans”). Therefore, the re-calibrated Fe-CN peak has about a −6 cm−1 change towards the lower energy direction.

3.6. Further Discussions

Figure 8 summarizes a few conceptual cases with different scan parameters but the same energy drift per scan (e.g., ΔE = 0.8 meV). The calibrated energy value for each data point (Ereal(k)) is calculated using Equation (2) and distributed according to different scanning parameters. For simpler consideration, we only consider the first calibration step, Eobs → E′ in the following discussions (Figure 8). The ΔEk is defined as the energy correction from the Eobs to E′ at one particular point (k). For different calibration procedures, the ΔEk should have different values at one particular point (k). However, their starting energy (E2′ = E2obs) and ending energy (E1′ = E1obs − ΔE) are always the same no matter how the energies between the two ends are calibrated. The blue curves in Figure 8a,b illustrate the cases with an even time NRVS scan as the reference to discuss other scans. Figure 8a starts with an extreme case: a “jump” scan with 30 s/p scanning time for 800 → 650 cm−1 and 0 s/p for the rest 650 → −200 cm−1 (red). The difference between this scan (red) and an even time scan (blue) is clear with the largest difference occurring at 650 cm−1 where the scanning time changes from 30 s/p to 0 s/p Two additional examples of “jump” scans are also provided in Figure 8a where the scanning time changes from 30 s/p to 0 s/p at 400 and 0 cm−1 respectively (the 2 black curves). As a general rule, the following “conclusions” hold:
(1)
it is obvious that the wider the scanning region (vs. the skipped region) [(800 → 650 cm−1) → (800 → 400 cm−1) → (800 → 0 cm−1)], the closer its energy distribution curve to the case of an even NRVS scan. For example, the left black curve is the closest one to the blue curve in Figure 8a;
(2)
the time ratio of the two adjacent scanning regions dictates the difference between the energies calculated via the Equation (1) and those calculated via the Equation (2): for example, the 5:10 [Figure 8b, light blue] provides a much closer result to the even time scan (blue) than the 1:30 (purple) or 30:0 (red) does. This leads to the concept that a multiple section scan that changes the scanning time at a gradual pace is closer to the even time scan. In practical measurements, it is better to start with an even time scan or a sectional scan but with multiple and stepwise changes in its scanning time parameters. When the major features are all well resolved and calibrated, some extremely weak features need to be probed with heavy counting on one region, such as the 30/1 or the 30/0 s/p scans discussed in Figure 8a. It becomes necessary to use the time-based energy calibration procedure (2) for any sectional scans but it is especially necessary for the ones with large time steps, such as the 30/1 or 30/0 ones;
(3)
no matter which calibration procedure is to be used and no matter what the scanning parameters, no point should have an energy drift amount greater than the total energy drift per scan (ΔEi). Therefore, the final factor for controlling the possible energy calibration errors is to scan NRVS with a lesser ΔEi value. For example, NRVS experimentalists need to avoid the moments right after each hutch opening or so and wait for the beam to be stabilized. For further reduction of ΔEi values, NRVS users also have the option to use less time for each scan and take more scans to average.
Figure 8. The energy shift per data point (ΔEk) distributed from an assumed total 0.8 meV energy drift per scan using the Equation (2) but for NRVS data with different scanning parameters: (a) an even time scan (blue) vs. the 30/0 s/p jump scan (red) which changes the scanning time at 650 (red), 400 (right black) and 0 (left black) cm−1 respectively; and (b) an even time scan (blue) vs. the 10/5/1 s/p sectional scan (light blue), the 30/1 s/p sectional scan (purple) and the 30/0 jump scan (red).
Figure 8. The energy shift per data point (ΔEk) distributed from an assumed total 0.8 meV energy drift per scan using the Equation (2) but for NRVS data with different scanning parameters: (a) an even time scan (blue) vs. the 30/0 s/p jump scan (red) which changes the scanning time at 650 (red), 400 (right black) and 0 (left black) cm−1 respectively; and (b) an even time scan (blue) vs. the 10/5/1 s/p sectional scan (light blue), the 30/1 s/p sectional scan (purple) and the 30/0 jump scan (red).
Physchem 02 00027 g008
Although the examples used in this publication are the data from 57Fe NRVS, NRVS for other isotopes (e.g., 125Te [30], 40K [52], 151Eu [29,53]) should have the same principle discussed here—energy should be calibrated according to the time scanned rather than the energy position travel because the energy drifts according to the beam-on time.

4. Summary

In this publication, we evaluated the energy calibration (Eobs → Ereal) with the time-scaled function within one NRVS scan according to the Equation (2). The errorbar contributed from the improper “distribution” of ΔEi within one scan (or ΔE for the averaged scan) (Eerr1) vs. that due to the different ΔEi (and thus different αi) from scan to scan (Eerr2) were analyzed and compared. It was found that the former (Eerr1) is as important as or sometimes even more important than the latter (Eerr2).
In case of need, a procedure [the Equation (11)] was established to re-calibrate previous NRVS spectra (PVDOS) with E*real as their energy axis to new spectra with Ereal as their energy axis. For practical purposes, one averaged ΔE value calculated from the old raw NRVS data is used to re-calibrate the old NRVS-derived PVDOS. Sectional scans with different scanning conditions as well as jump scans where an “unimportant” energy region is skipped from the scanning are also discussed in detail.
In the end, via this article, the concept can be established that energy positions in one sectional NRVS scan should be corrected or re-distributed according to the time scanned [as in the Equations (2) or (11)] rather than be scaled with a universal constant [as in the Equation (1)].
Although the Formula (2) itself was mentioned in an earlier publication [40], it is this publication that elaborated on the meaningfulness and the applications of the energy re-distribution within one NRVS scan. Several published NRVS spectra were re-calibrated with a few cm−1 energy shifts in comparison with the published spectra (PVDOS) calibrated with a traditional energy calibration (1). These spectra are shown as in Figure 4, Figure 5, Figure 6 and Figure 7.
The transition from the old calibration procedure (1) to the new stepwise calibration practice takes time and some users may prefer to continue using the old calibration procedure. Even without tracking the energy variation from scan to scan due to different ΔEi {Wang, 2022 #55}, using an averaged ΔE for all the NRVS scans and redistributing this averaged ΔE according to the time scanned [using the Equation (2)] rather than according to a universal scale [using the Equation (1)] can fix the Eerr1 and reduce the total errorbar by ~50%.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/physchem2040027/s1, Figure S1: A series of NRVS spectra for 57Fe powder: the E0 position drifts to the higher-energy direction as scanning time increases (red → orange → yellow → green → cyan → dark blue → purple).

Author Contributions

H.W. and Y.Y. measured the NRVS data and proposed the initial concept, J.W. streamlined and tested the details in the mathematical modeling. All the authors contribute equally. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

We thank Stephen P. Cramer at SETI Institute (USA) for his overall support for this project via US NIH grants GM-65440. The NRVS spectra used, cited, mentioned, or tested for this publication were measured at SPring-8 BL09XU or BL19LXU (e.g., proposals 20200013, 20210033 via RIKEN, and 2018A1409, 2018B1379, 2018A1033, 2019A1259 via JASRI).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Seto, M.; Yoda, Y.; Kikuta, S.; Zhang, X.W.; Ando, M. Observation of Nuclear Resonant Scattering Accompanied by Phonon Excitation Using Synchrotron Radiation. Phys. Rev. Lett. 1995, 74, 3828–3831. [Google Scholar] [CrossRef]
  2. Sturhahn, W.; Toellner, T.S.; Alp, E.E.; Zhang, X.; Ando, M.; Yoda, Y.; Kikuta, S.; Seto, M.; Kimball, C.W.; Dabrowski, B. Phonon Density of States Measured by Inelastic Nuclear Resonant Scattering. Phys. Rev. Lett. 1995, 74, 3832–3835. [Google Scholar] [CrossRef] [Green Version]
  3. Yoda, Y.; Yabashi, M.; Izumi, K.; Zhang, X.W.; Kishimoto, S.; Kitao, S.; Seto, M.; Mitsui, T.; Harami, T.; Imai, Y.; et al. Nuclear resonant scattering beamline at SPring-8. Nucl. Inst. Meth. A 2001, 467, 715–718. [Google Scholar] [CrossRef]
  4. Yoda, Y.; Imai, Y.; Kobayashi, H.; Goto, S.; Takeshita, K.; Seto, M. Upgrade of the nuclear resonant scattering beamline, BL09XU in SPring-8. Hyperfine Interact. 2012, 206, 83. [Google Scholar] [CrossRef]
  5. Wang, H.; Yoda, Y.; Dong, W.; Huang, S.D. Energy calibration issues in nuclear resonant vibrational spectroscopy: Observing small spectral shifts and making fast calibrations. J. Synchrotron Radiat. 2013, 20, 683–690. [Google Scholar] [CrossRef] [PubMed]
  6. Wang, H.; Alp, E.; Yoda, Y.; Cramer, S. A Practical Guide for Nuclear Resonance Vibrational Spectroscopy (NRVS) of Biochemical Samples and Model Compounds. In Metalloproteins; Fontecilla-Camps, J.C., Nicolet, Y., Eds.; Humana Press: Totowa, NJ, USA, 2014; Volume 1122, pp. 125–137. [Google Scholar]
  7. Yoda, Y.; Okada, K.; Wang, H.; Cramer, S.P.; Seto, M. High-resolution monochromator for iron nuclear resonance vibrational spectroscopy of biological samples. Jap. J. App. Phys. 2016, 55, 122401. [Google Scholar] [CrossRef] [PubMed]
  8. Wang, H.; Braun, A.; Cramer, S.P.; Gee, L.B.; Yoda, Y. Nuclear Resonance Vibrational Spectroscopy: A Modern Tool to Pinpoint Site-Specific Cooperative Processes. Crystals 2021, 11, 909. [Google Scholar] [CrossRef] [PubMed]
  9. Wang, J.; Li, L.; Wang, H. Machine learning concept in de-spiking process for nuclear resonant vibrational spectra—Automation using no external parameter. Vib. Spectrosc. 2022, 119, 103352. [Google Scholar] [CrossRef]
  10. Wang, H.; Yoda, Y.; Ogata, H.; Tanaka, Y.; Lubitz, W. A strenuous experimental journey searching for spectroscopic evidence of a bridging nickel-iron-hydride in [NiFe] hydrogenase. J. Synchrotron Radiat. 2015, 22, 1334–1344. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Wang, H.; Salzberg, A.P.; Weiner, B.R. Laser ablation of aluminum at 193, 248, and 351 nm. Appl. Phys. Lett. 1991, 59, 935–937. [Google Scholar] [CrossRef]
  12. He, M.; Wang, H.; Weiner, B.R. Production and laser-induced fluorescence spectrum of aluminum sulfide. Chem. Phys. Lett. 1993, 204, 563–566. [Google Scholar] [CrossRef]
  13. Wang, H.; Chen, X.; Weiner, B.R. Laser photodissociation dynamics of thionyl chloride: Concerted and stepwise cleavage of S-Cl bonds. J. Phys. Chem. 1993, 97, 12260–12268. [Google Scholar] [CrossRef]
  14. Chen, X.; Asmar, F.; Wang, H.; Weiner, B.R. Nascent sulfur monoxide (X3.SIGMA.-) vibrational distributions from the photodissociation of sulfur dioxide, sulfonyl chloride, and dimethylsulfoxide at 193 nm. J. Phys. Chem. 1991, 95, 6415–6417. [Google Scholar] [CrossRef]
  15. Guo, Y.; Wang, H.; Xiao, Y.; Vogt, S.; Thauer, R.K.; Shima, S.; Volkers, P.I.; Rauchfuss, T.B.; Pelmenschikov, V.; Case, D.A.; et al. Characterization of the Fe site in iron-sulfur cluster-free hydrogenase (Hmd) and of a model compound via nuclear resonance vibrational spectroscopy (NRVS). Inorg. Chem. 2008, 47, 3969–3977. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Smith, M.C.; Xiao, Y.; Wang, H.; George, S.J.; Coucouvanis, D.; Koutmos, M.; Sturhahn, W.; Alp, E.E.; Zhao, J.; Cramer, S.P. Normal-Mode Analysis of FeCl4- and Fe2S2Cl42- via Vibrational Mössbauer, Resonance Raman, and FT-IR Spectroscopies. Inorg. Chem. 2005, 44, 5562–5570. [Google Scholar] [CrossRef]
  17. Xiao, Y.; Wang, H.; George, S.J.; Smith, M.C.; Adams, M.W.; Jenney, F.E., Jr.; Sturhahn, W.; Alp, E.E.; Zhao, J.; Yoda, Y.; et al. Normal mode analysis of Pyrococcus furiosus rubredoxin via nuclear resonance vibrational spectroscopy (NRVS) and resonance raman spectroscopy. J. Am. Chem. Soc. 2005, 127, 14596–14606. [Google Scholar] [CrossRef]
  18. Xiao, Y.; Fisher, K.; Smith, M.C.; Newton, W.; Case, D.A.; George, S.J.; Wang, H.; Sturhahn, W.; Alp, E.E.; Zhao, J.; et al. How Nitrogenase Shakes—Initial Information about P-Cluster and FeMo-Cofactor Normal Modes from Nuclear Resonance Vibrational Spectroscopy (NRVS). J. Am. Chem. Soc. 2006, 128, 7608–7612. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Gee, L.B.; Pelmenschikov, V.; Wang, H.; Mishra, N.; Liu, Y.-C.; Yoda, Y.; Tamasaku, K.; Chiang, M.-H.; Cramer, S.P. Vibrational characterization of a diiron bridging hydride complex—A model for hydrogen catalysis. Chem. Sci. 2020, 11, 5487–5493. [Google Scholar] [CrossRef] [PubMed]
  20. Gee, L.B.; Wang, H.; Cramer, S.P. Chapter Fourteen—NRVS for Fe in Biology: Experiment and Basic Interpretation. In Methods in Enzymology; David, S.S., Ed.; Academic Press: Cambridge, MA, USA, 2018; Volume 599, pp. 409–425. [Google Scholar]
  21. Guo, Y.; Brecht, E.; Aznavour, K.; Nix, J.C.; Xiao, Y.; Wang, H.; George, S.J.; Bau, R.; Keable, S.; Peters, J.W.; et al. Nuclear resonance vibrational spectroscopy (NRVS) of rubredoxin and MoFe protein crystals. Hyperfine Interact. 2012, 222, 77–90. [Google Scholar] [CrossRef]
  22. Ogata, H.; Kramer, T.; Wang, H.; Schilter, D.; Pelmenschikov, V.; van Gastel, M.; Neese, F.; Rauchfuss, T.B.; Gee, L.B.; Scott, A.D.; et al. Hydride bridge in [NiFe]-hydrogenase observed by nuclear resonance vibrational spectroscopy. Nat. Commun. 2015, 6, 7890. [Google Scholar] [CrossRef] [Green Version]
  23. Kamali, S.; Wang, H.; Mitra, D.; Ogata, H.; Lubitz, W.; Manor, B.C.; Rauchfuss, T.B.; Byrne, D.; Bonnefoy, V.; Jenney, F.E.; et al. Observation of the Fe-CN and Fe-CO Vibrations in the Active Site of [NiFe] Hydrogenase by Nuclear Resonance Vibrational Spectroscopy. Angew. Chem. Int. Ed. 2013, 52, 724–728. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Tinberg, C.E.; Tonzetich, Z.J.; Wang, H.; Do, L.H.; Yoda, Y.; Cramer, S.P.; Lippard, S.J. Characterization of Iron Dinitrosyl Species Formed in the Reaction of Nitric Oxide with a Biological Rieske Center. J. Am. Chem. Soc. 2010, 132, 18168–18176. [Google Scholar] [CrossRef] [PubMed]
  25. Do, L.H.; Wang, H.; Tinberg, C.E.; Dowty, E.; Yoda, Y.; Cramer, S.P.; Lippard, S.J. Characterization of a synthetic peroxodiiron(III) protein model complex by nuclear resonance vibrational spectroscopy. Chem. Commun. 2011, 47, 10945–10947. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Gilbert-Wilson, R.; Siebel, J.F.; Adamska-Venkatesh, A.; Pham, C.C.; Reijerse, E.; Wang, H.; Cramer, S.P.; Lubitz, W.; Rauchfuss, T.B. Spectroscopic Investigations of [FeFe] Hydrogenase Maturated with [57Fe2(adt)(CN)2(CO)4]2−. J. Am. Chem. Soc. 2015, 137, 8998–9005. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Pelmenschikov, V.; Birrell, J.A.; Pham, C.C.; Mishra, N.; Wang, H.X.; Sommer, C.; Reijerse, E.; Richers, C.P.; Tamasaku, K.; Yoda, Y.; et al. Reaction Coordinate Leading to H-2 Production in [FeFe]-Hydrogenase Identified by Nuclear Resonance Vibrational Spectroscopy and Density Functional Theory. J. Am. Chem. Soc. 2017, 139, 16894–16902. [Google Scholar] [CrossRef] [PubMed]
  28. Pham, C.C.; Mulder, D.W.; Pelmenschikov, V.; King, P.W.; Ratzloff, M.W.; Wang, H.; Mishra, N.; Alp, E.E.; Zhao, J.; Hu, M.Y.; et al. Terminal Hydride Species in [FeFe]-Hydrogenases are Vibrationally Coupled to the Active Site Environment. Ang. Chem. 2018, 57, 10605–10609. [Google Scholar] [CrossRef] [PubMed]
  29. Wang, H.; Huang, S.D.; Yan, L.; Hu, M.Y.; Zhao, J.; Alp, E.E.; Yoda, Y.; Petersen, C.M.; Thompson, M.K. Europium-151 and iron-57 nuclear resonant vibrational spectroscopy of naturally abundant KEu(iii)Fe(ii)(CN)6 and Eu(iii)Fe(iii)(CN)6 complexes. Dalton Trans. 2022. [Google Scholar] [CrossRef] [PubMed]
  30. Wittkamp, F.; Mishra, N.; Wang, H.; Wille, H.-C.; Steinbrügge, R.; Kaupp, M.; Cramer, S.P.; Apfel, U.-P.; Pelmenschikov, V. Insights from 125Te and 57Fe nuclear resonance vibrational spectroscopy: A [4Fe–4Te] cluster from two points of view. Chem. Sci. 2019, 10, 7535–7541. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Pelmenschikov, V.; Birrell, J.A.; Gee, L.B.; Richers, C.P.; Reijerse, E.J.; Wang, H.; Arragain, S.; Mishra, N.; Yoda, Y.; Matsuura, H.; et al. Vibrational Perturbation of the [FeFe] Hydrogenase H-Cluster Revealed by 13C2H-ADT Labeling. J. Am. Chem. Soc. 2021, 143, 8237–8243. [Google Scholar] [CrossRef]
  32. Gu, W.W.; Jacquamet, L.; Patil, D.S.; Wang, H.X.; Evans, D.J.; Smith, M.C.; Millar, M.; Koch, S.; Eichhorn, D.M.; Latimer, M.; et al. Refinement of the nickel site structure in Desulfovibrio gigas hydrogenase using range-extended EXAFS spectroscopy. J. Inorg. Biochem. 2003, 93, 41–51. [Google Scholar] [CrossRef]
  33. DeBeer George, S.; Metz, M.; Szilagyi, R.K.; Wang, H.; Cramer, S.P.; Lu, Y.; Tolman, W.B.; Hedman, B.; Hodgson, K.O.; Solomon, E.I. A Quantitative Description of the Ground-State Wave Function of CuA by X-ray Absorption Spectroscopy:  Comparison to Plastocyanin and Relevance to Electron Transfer. J. Am. Chem. Soc. 2001, 123, 5757–5767. [Google Scholar] [CrossRef] [PubMed]
  34. O’Dowd, B.; Williams, S.; Wang, H.X.; No, J.H.; Rao, G.D.; Wang, W.X.; McCammon, J.A.; Cramer, S.P.; Oldfield, E. Spectroscopic and Computational Investigations of Ligand Binding to IspH: Discovery of Non-diphosphate Inhibitors. Chembiochem 2017, 18, 914–920. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Funk, T.; Deb, A.; George, S.J.; Wang, H.; Cramer, S.P. X-ray magnetic circular dichroism—A high energy probe of magnetic properties. Coord. Chem. Rev. 2005, 249, 3–30. [Google Scholar] [CrossRef]
  36. Wang, H.X.; Ralston, C.Y.; Patil, D.S.; Jones, R.M.; Gu, W.; Verhagen, M.; Adams, M.; Ge, P.; Riordan, C.; Marganian, C.A.; et al. Nickel L-edge soft X-ray spectroscopy of nickel-iron hydrogenases and model compounds—Evidence for high-spin nickel(II) in the active enzyme. J. Am. Chem. Soc. 2000, 122, 10544–10552. [Google Scholar] [CrossRef]
  37. Wang, H.X.; Ge, P.H.; Riordan, C.G.; Brooker, S.; Woomer, C.G.; Collins, T.; Melendres, C.A.; Graudejus, O.; Bartlett, N.; Cramer, S.P. Integrated X-ray L absorption spectra. Counting holes in Ni complexes. J. Phys. Chem. B 1998, 102, 8343–8346. [Google Scholar] [CrossRef]
  38. Wang, H.; Friedrich, S.; Li, L.; Mao, Z.; Ge, P.; Balasubramanian, M.; Patil, D.S. L-edge sum rule analysis on 3d transition metal sites: From d10 to d0 and towards application to extremely dilute metallo-enzymes. Phys. Chem. Chem. Phys. 2018, 20, 8166–8176. [Google Scholar] [CrossRef] [PubMed]
  39. Pierce and Rod. Accuracy and Precision. Available online: https://www.mathsisfun.com/accuracy-precision.html (accessed on 1 June 2022).
  40. Wang, J.; Yoda, Y.; Wang, H. Tracking energy scale variations from scan to scan in nuclear resonant vibrational spectroscop y: In situ correction using zero-energy position drifts ΔEi rather than making in situ calibration measurements. Rev. Sci. Instrum. 2022, 93, 095101. [Google Scholar] [CrossRef] [PubMed]
  41. Zhao, J.Y.; Sturhahn, W. High-energy-resolution X-ray monochromator calibration using the detailed-balance principle. J. Synchrotron Radiat. 2012, 19, 602–608. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Sturhahn, W. Nuclear resonant spectroscopy. J. Phys. Condens. Matter 2004, 16, S497–S530. [Google Scholar] [CrossRef]
  43. Sturhahn, W. CONUSS and PHOENIX: Evaluation of nuclear resonant scattering data. Hyperfine Interact. 2000, 125, 149–172. [Google Scholar] [CrossRef]
  44. Watanabe, H.; Yamada, N.; Okaji, M. Linear Thermal Expansion Coefficient of Silicon from 293 to 1000 K. Int. J. Thermophys. 2004, 25, 221–236. [Google Scholar] [CrossRef]
  45. Lumen-Physics. Heat and Heat Transfer Methods. Available online: https://courses.lumenlearning.com/physics/chapter/14-5-conduction/ (accessed on 1 June 2022).
  46. Resnick, R.; Halliday, D.; Walker, J. Fundamentals of Physics Extended; Welly: New York, NY, USA, 2016. [Google Scholar]
  47. Morra, S. Fantastic [FeFe]-Hydrogenases and Where to Find Them. Front. Microbiol. 2022, 13, 853626. [Google Scholar] [CrossRef] [PubMed]
  48. Chiang, M.-H.; Pelmenschikov, V.; Gee, L.B.; Liu, Y.-C.; Hsieh, C.-C.; Wang, H.; Yoda, Y.; Matsuura, H.; Li, L.; Cramer, S.P. High-Frequency Fe–H and Fe–H2 Modes in a trans-Fe(η2-H2)(H) Complex: A Speed Record for Nuclear Resonance Vibrational Spectroscopy. Inorg. Chem. 2021, 60, 555–559. [Google Scholar] [CrossRef] [PubMed]
  49. Chumakov, A.I.; Rüffer, R.; Leupold, O.; Sergueev, I. Insight to Dynamics of Molecules with Nuclear Inelastic Scattering. Struct. Chem. 2003, 14, 109–119. [Google Scholar] [CrossRef]
  50. Shokouhimehr, M.; Soehnlen, E.S.; Hao, J.; Griswold, M.; Flask, C.; Fan, X.; Basilion, J.P.; Basu, S.; Huang, S.D. Dual purpose Prussian blue nanoparticles for cellular imaging and drug delivery: A new generation of T-1-weighted MRI contrast and small molecule delivery agents. J. Mater. Chem. 2010, 20, 5251–5259. [Google Scholar] [CrossRef]
  51. Kandanapitiye, M.S.; Valley, B.; Yang, L.D.; Fry, A.M.; Woodward, P.M.; Huang, S.D. Gallium Analogue of Soluble Prussian Blue KGa Fe(CN)(6) center dot nH(2)O: Synthesis, Characterization, and Potential Biomedical Applications. Inorg. Chem. 2013, 52, 2790–2792. [Google Scholar] [CrossRef] [PubMed]
  52. Seto, M.; Kitao, S.; Kobayashi, Y.; Haruki, R.; Mitsui, T.; Yoda, Y.; Zhang, X.W.; Kishimoto, S.; Maeda, Y. Nuclear Resonant Inelastic and Forward Scattering of Synchrotron Radiation by 40K. Hyperfine Interact. 2002, 141, 99–108. [Google Scholar] [CrossRef]
  53. Kobayashi, H.; Yoda, Y.; Shirakawa, M.; Ochiai, A. 151Eu Nuclear Resonant Inelastic Scattering of Eu4As3 around Charge Ordering Temperature. J. Phys. Soc. Jpn. 2006, 75, 034602. [Google Scholar] [CrossRef]
Figure 1. (a) An illustrative energy diagram for NRVS transitions; (b) the uncalibrated NRVS spectra for (NEt4)[FeCl4] measured at four NRS beamlines around the world: APS 03ID (red), ESRF ID18 (blue); SPring-8 BL09XU (green), and Petra–III P01 (black). The top inserts in (b) are the overviews for the four SR centers.
Figure 1. (a) An illustrative energy diagram for NRVS transitions; (b) the uncalibrated NRVS spectra for (NEt4)[FeCl4] measured at four NRS beamlines around the world: APS 03ID (red), ESRF ID18 (blue); SPring-8 BL09XU (green), and Petra–III P01 (black). The top inserts in (b) are the overviews for the four SR centers.
Physchem 02 00027 g001
Figure 2. A conceptual NRVS measurement which scans from the higher energy end (E2′) to the lower energy end (E1′) and which scans points 1, 2, 3, 4 (red) with 3t s/p and points 5, 6, 7, 8 (blue) with t s/p: (a) a targeted theoretical range for the NRVS spectrum (in E′); (b) the observed raw energies (at points 1−8) when we assume an exaggerated larges ΔE which is half of (E2′ − E1′)—that is (Eobs2 − Eobs1) = (½)·(E2′ − E1′); (c) the corrected intermediate energies (E′) which map the large amount ΔE onto the 8 data points using universal energy scale Formula (1); (d) the corrected intermediate energies (E′) which map the large amount ΔE onto the 8 data points using time-scaled correlation (2); and (d’) the corrected energies in (d) with the E0 position aligned to 0.
Figure 2. A conceptual NRVS measurement which scans from the higher energy end (E2′) to the lower energy end (E1′) and which scans points 1, 2, 3, 4 (red) with 3t s/p and points 5, 6, 7, 8 (blue) with t s/p: (a) a targeted theoretical range for the NRVS spectrum (in E′); (b) the observed raw energies (at points 1−8) when we assume an exaggerated larges ΔE which is half of (E2′ − E1′)—that is (Eobs2 − Eobs1) = (½)·(E2′ − E1′); (c) the corrected intermediate energies (E′) which map the large amount ΔE onto the 8 data points using universal energy scale Formula (1); (d) the corrected intermediate energies (E′) which map the large amount ΔE onto the 8 data points using time-scaled correlation (2); and (d’) the corrected energies in (d) with the E0 position aligned to 0.
Physchem 02 00027 g002
Figure 3. (a) An example scanning time per data point (s/p) used in a conceptual NRVS sectional scan between −200 cm−1 and +800 cm−1: the scanning time is not to the scale but it is noted on the grey horizontal line; (b) the energy shift per data point assuming a total 0.8 meV energy drift per scan: the energy drift at each point calibrated with Equation (9) (blue) vs. that calibrated according to the accumulated time scanned with (2) (red); (c) the results from (b) with the zero vibrational peak position moving to 0 cm−1 in both calibration methods; (d) a conceptual NRVS spectral section between 625 and 705 cm−1; and (e) the NRVS spectra calibrated with the two calibration methods used in (c) (color match).
Figure 3. (a) An example scanning time per data point (s/p) used in a conceptual NRVS sectional scan between −200 cm−1 and +800 cm−1: the scanning time is not to the scale but it is noted on the grey horizontal line; (b) the energy shift per data point assuming a total 0.8 meV energy drift per scan: the energy drift at each point calibrated with Equation (9) (blue) vs. that calibrated according to the accumulated time scanned with (2) (red); (c) the results from (b) with the zero vibrational peak position moving to 0 cm−1 in both calibration methods; (d) a conceptual NRVS spectral section between 625 and 705 cm−1; and (e) the NRVS spectra calibrated with the two calibration methods used in (c) (color match).
Physchem 02 00027 g003
Figure 4. (a) NRVS PVDOS in the X–Fe–H region for Cr [FeFe] hydrogenase in the Hhyd form calibrated with the traditional procedure (1): green = ODT, blue = wild type, grey = C169S; and (b) NRVS PVDOS for the sample spectrum of wild type enzyme re-calibrated with the time-scaled procedure (2) (red) vs. that calibrated with the traditional procedure (1) (dashed blue): the duplication of the wild type spectrum in 4a.
Figure 4. (a) NRVS PVDOS in the X–Fe–H region for Cr [FeFe] hydrogenase in the Hhyd form calibrated with the traditional procedure (1): green = ODT, blue = wild type, grey = C169S; and (b) NRVS PVDOS for the sample spectrum of wild type enzyme re-calibrated with the time-scaled procedure (2) (red) vs. that calibrated with the traditional procedure (1) (dashed blue): the duplication of the wild type spectrum in 4a.
Physchem 02 00027 g004
Figure 5. (a) NRVS PVDOS of complex trans−[57Fe(η2−H2)(H)(dppe)2][BPh4] (H2FeH for short) in the region for Fe–H stretching modes calibrated with the traditional procedure (1) (blue) and with new time-scaled redistribution method (2) (red); (b) PVDOS for (a) in the narrower region of 1720−1820 cm−1; and (c) PVDOS for (a) in the narrower region of 1885–1985 cm−1.
Figure 5. (a) NRVS PVDOS of complex trans−[57Fe(η2−H2)(H)(dppe)2][BPh4] (H2FeH for short) in the region for Fe–H stretching modes calibrated with the traditional procedure (1) (blue) and with new time-scaled redistribution method (2) (red); (b) PVDOS for (a) in the narrower region of 1720−1820 cm−1; and (c) PVDOS for (a) in the narrower region of 1885–1985 cm−1.
Physchem 02 00027 g005
Figure 6. The FeCO feature [in DvMF NiFe hydrogenase NRVS] in the lower energy side of the 30 s/p to 1 s/p timing changing point: (a) the NRVS spectrum calibrated with the traditional calibration procedure; and (b) the NRVS spectrum re-calibrated with the time-scaled calibration procedure (2).
Figure 6. The FeCO feature [in DvMF NiFe hydrogenase NRVS] in the lower energy side of the 30 s/p to 1 s/p timing changing point: (a) the NRVS spectrum calibrated with the traditional calibration procedure; and (b) the NRVS spectrum re-calibrated with the time-scaled calibration procedure (2).
Physchem 02 00027 g006
Figure 7. The NRVS spectra for a series of Fe(CN)6 containing complexes: (a) PVDOS for (NH4)2MgFe(CN)6; (b) raw NRVS spectrum for KEuFe(CN)6 which was measured with a jump scan procedure around the Fe(CN) region and the nuclear resonant peak region. The blue spectrum was calibrated with the tradition procedure (1) and the red one is re-calibrated with a time-based procedure using (2); and (c) the scanning time per point (s/p) used in different energy regions during the NRVS measurement for KEuFe(CN)6 (b).
Figure 7. The NRVS spectra for a series of Fe(CN)6 containing complexes: (a) PVDOS for (NH4)2MgFe(CN)6; (b) raw NRVS spectrum for KEuFe(CN)6 which was measured with a jump scan procedure around the Fe(CN) region and the nuclear resonant peak region. The blue spectrum was calibrated with the tradition procedure (1) and the red one is re-calibrated with a time-based procedure using (2); and (c) the scanning time per point (s/p) used in different energy regions during the NRVS measurement for KEuFe(CN)6 (b).
Physchem 02 00027 g007
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, H.; Yoda, Y.; Wang, J. The True Nature of the Energy Calibration for Nuclear Resonant Vibrational Spectroscopy: A Time-Based Conversion. Physchem 2022, 2, 369-388. https://doi.org/10.3390/physchem2040027

AMA Style

Wang H, Yoda Y, Wang J. The True Nature of the Energy Calibration for Nuclear Resonant Vibrational Spectroscopy: A Time-Based Conversion. Physchem. 2022; 2(4):369-388. https://doi.org/10.3390/physchem2040027

Chicago/Turabian Style

Wang, Hongxin, Yoshitaka Yoda, and Jessie Wang. 2022. "The True Nature of the Energy Calibration for Nuclear Resonant Vibrational Spectroscopy: A Time-Based Conversion" Physchem 2, no. 4: 369-388. https://doi.org/10.3390/physchem2040027

APA Style

Wang, H., Yoda, Y., & Wang, J. (2022). The True Nature of the Energy Calibration for Nuclear Resonant Vibrational Spectroscopy: A Time-Based Conversion. Physchem, 2(4), 369-388. https://doi.org/10.3390/physchem2040027

Article Metrics

Back to TopTop