# Site-specific carbon isotope measurements of vanillin reference materials by nuclear magnetic resonance spectrometry

Phuong Mai Le, Estelle Martineau, Serge Akoka, Gerald Remaud, Michelle Chartrand, Juris Meija, Zoltán Mester

<figure><img src="image_1.png" /></figure>

To cite this version:

Phuong Mai Le, Estelle Martineau, Serge Akoka, Gerald Remaud, Michelle Chartrand, et al. Site-specific carbon isotope measurements of vanillin reference materials by nuclear magnetic resonance spectrometry. Analytical and Bioanalytical Chemistry, 2022, 414 (24), pp.7153-7165. 10.1007/s00216-022-04292-0. hal-04417081

HAL Id: hal-04417081

https://hal.science/hal-04417081

Submitted on 30 May 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Site-specific carbon isotope measurements of vanillin reference materials by nuclear magnetic resonance spectrometry

Phuong Mai Le, Estelle Martineau, Serge Akoka, Gérald Remaud, Michelle

Chartrand, Juris Meija, Zoltán Mester

<figure><img src="image_1.png" /></figure>

To cite this version:

Phuong Mai Le, Estelle Martineau, Serge Akoka, Gérald Remaud, Michelle Chartrand, et al. Site-specific carbon isotope measurements of vanillin reference materials by nuclear magnetic resonance spectrometry. Analytical and Bioanalytical Chemistry, 2022, 414 (24), pp.7153-7165. 10.1007/s00216-022-04292-0. hal-04417081

HAL Id: hal-04417081

https://hal.science/hal-04417081v1

Submitted on 30 May 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Site-specific carbon isotope measurements of vanillin reference materials by nuclear magnetic resonance spectrometry

Phuong Mai Le<fnref n="1" />☐, Estelle Martineau<fnref n="2" />,<fnref n="3" />, Serge Akoka<fnref n="2" />, Gerald Remaud<fnref n="2" />, Michelle M. G. Chartrand<fnref n="1" />, Juris Meija<fnref n="1" />, Zoltán Mester<fnref n="1" />

Received: 12 July 2022 / Revised: 8 August 2022 / Accepted: 16 August 2022 / Published online: 13 September 2022
© Crown 2022

## Abstract

Vanillin, one of the world's most popular flavor used in food and pharmaceutical industries, is extracted from vanilla beans or obtained (bio)-synthetically. The price of natural vanillin is considerably higher than that of its synthetic alternative which leads increasingly to counterfeit vanillin. Here, we describe the workflow of combining carbon isotope ratio combustion mass spectrometry with quantitative carbon nuclear magnetic resonance spectrometry (¹³C-qNMR) to obtain carbon isotope measurements traceable to the Vienna Pee Dee Belemnite (VPDB) with 0.7‰ combined standard uncertainty (or expanded uncertainty of 1.4‰ at 95% confidence level). We perform these measurements on qualified Bruker 400 MHz instruments to certify site-specific carbon isotope delta values in two vanillin materials, VANA-1 and VANB-1, believed to be the first intramolecular isotopic certified reference material (CRMs).

**Keywords** Site-specific isotope ratios · Carbon · ¹³C-qNMR · Vanillin · Certified reference material

## Introduction

Vanillin is one of the most common flavors used in food, cosmetics, and pharmaceuticals. Due to the high cost and demand for vanillin, synthetic vanillin is common. It is estimated that more than 90% of vanillin is from synthetic origins [1], made largely from petrochemicals and by-products from the paper industry (see Fig. 1) [2]. Consumers are increasingly interested in the authenticity and traceability of food products. Chromatographic techniques such as GC and HPLC rely on composition profiles of impurities to distinguish between natural and synthetic vanillin [3] but are unable to detect sophisticated counterfeit methods which add flavoring substances found in vanilla [4].

Stable isotope analysis is valuable for detecting food counterfeiting [6, 7]. Isotope ratio mass spectrometry (IRMS) is commonly used to determine the bulk isotopic composition for a given molecule whereas isotope ratio NMR is capable to determine the site-specific isotopic composition of a given compound.

Isotopic analysis of hydrogen, carbon, and oxygen has been used for the authentication of vanillin products [4, 8–11], with the power of this tool increasing when hydrogen and carbon isotopic compositions are combined [10–12]. Despite this, IRMS provides access to bulk, not to the intramolecular ¹³C composition of the analyte, resulting in some limitations in the authentication of bio- or chemical-synthetic vanillin obtained from the same precursor [9, 11] and detection of the addition of synthetic vanillin to natural vanillin [13]. In contrast, quantitative NMR reveals all spectrally resolved sites in the compound and provides position-specific isotope ratio measurements [14], which enables a more nuanced look into the formation of a molecule by differentiating the sources of the individual carbon atoms [15]. The absolute carbon isotope (¹³C/¹²C) ratio of VPDB references such as glycine reference materials (USGS64, USGS65, and USGS66) has been recently measured using the signals associated with ¹²C and ¹³C (¹H-¹²C and ¹H-¹³C signals) in the proton nuclear magnetic resonance spectroscopy (¹H

<footnote>☒ Phuong Mai Le
PhuongMai.Le@nrc-cnrc.gc.ca</footnote>

<footnote>☒ Juris Meija
Juris.Meija@nrc-cnrc.gc.ca</footnote>

<footnote n="1">¹ Metrology, National Research Council Canada, 1200
Montreal Road, Ottawa, ON K1A 0R6, Canada</footnote>

<footnote n="2">² Nantes Université, CNRS, CEISAM, UMR6230,
F-44000 Nantes, France</footnote>

<footnote n="3">³ CAPACITÉS SAS, Nantes, France</footnote>
# Site-specific carbon isotope measurements of vanillin reference materials by nuclear magnetic resonance spectrometry

Phuong Mai Le<fnref n="1" />☐, Estelle Martineau<fnref n="2" />,<fnref n="3" />, Serge Akoka<fnref n="2" />, Gerald Remaud<fnref n="2" />, Michelle M. G. Chartrand<fnref n="1" />, Juris Meija<fnref n="1" />, Zoltán Mester<fnref n="1" />

Received: 12 July 2022 / Revised: 8 August 2022 / Accepted: 16 August 2022 / Published online: 13 September 2022
© Crown 2022

## Abstract

Vanillin, one of the world's most popular flavor used in food and pharmaceutical industries, is extracted from vanilla beans or obtained (bio)-synthetically. The price of natural vanillin is considerably higher than that of its synthetic alternative which leads increasingly to counterfeit vanillin. Here, we describe the workflow of combining carbon isotope ratio combustion mass spectrometry with quantitative carbon nuclear magnetic resonance spectrometry (¹³C-qNMR) to obtain carbon isotope measurements traceable to the Vienna Pee Dee Belemnite (VPDB) with 0.7‰ combined standard uncertainty (or expanded uncertainty of 1.4‰ at 95% confidence level). We perform these measurements on qualified Bruker 400 MHz instruments to certify site-specific carbon isotope delta values in two vanillin materials, VANA-1 and VANB-1, believed to be the first intramolecular isotopic certified reference material (CRMs).

**Keywords** Site-specific isotope ratios · Carbon · ¹³C-qNMR · Vanillin · Certified reference material

## Introduction

Vanillin is one of the most common flavors used in food, cosmetics, and pharmaceuticals. Due to the high cost and demand for vanillin, synthetic vanillin is common. It is estimated that more than 90% of vanillin is from synthetic origins [1], made largely from petrochemicals and by-products from the paper industry (see Fig. 1) [2]. Consumers are increasingly interested in the authenticity and traceability of food products. Chromatographic techniques such as GC and HPLC rely on composition profiles of impurities to distinguish between natural and synthetic vanillin [3] but are unable to detect sophisticated counterfeit methods which add flavoring substances found in vanilla [4].

Stable isotope analysis is valuable for detecting food counterfeiting [6, 7]. Isotope ratio mass spectrometry (IRMS) is commonly used to determine the bulk isotopic composition for a given molecule whereas isotope ratio NMR is capable to determine the site-specific isotopic composition of a given compound.

Isotopic analysis of hydrogen, carbon, and oxygen has been used for the authentication of vanillin products [4, 8–11], with the power of this tool increasing when hydrogen and carbon isotopic compositions are combined [10–12]. Despite this, IRMS provides access to bulk, not to the intramolecular ¹³C composition of the analyte, resulting in some limitations in the authentication of bio- or chemical-synthetic vanillin obtained from the same precursor [9, 11] and detection of the addition of synthetic vanillin to natural vanillin [13]. In contrast, quantitative NMR reveals all spectrally resolved sites in the compound and provides position-specific isotope ratio measurements [14], which enables a more nuanced look into the formation of a molecule by differentiating the sources of the individual carbon atoms [15]. The absolute carbon isotope (¹³C/¹²C) ratio of VPDB references such as glycine reference materials (USGS64, USGS65, and USGS66) has been recently measured using the signals associated with ¹²C and ¹³C (¹H-¹²C and ¹H-¹³C signals) in the proton nuclear magnetic resonance spectroscopy (¹H

<footnote>☒ Phuong Mai Le
PhuongMai.Le@nrc-cnrc.gc.ca</footnote>

<footnote>☒ Juris Meija
Juris.Meija@nrc-cnrc.gc.ca</footnote>

<footnote n="1">¹ Metrology, National Research Council Canada, 1200
Montreal Road, Ottawa, ON K1A 0R6, Canada</footnote>

<footnote n="2">² Nantes Université, CNRS, CEISAM, UMR6230,
F-44000 Nantes, France</footnote>

<footnote n="3">³ CAPACITÉS SAS, Nantes, France</footnote>
<figure><img src="image_3.png" /><figcaption>**Fig. 1** Principal sources of vanillin [5]</figcaption></figure>

NMR) [16]. This method provides a new and independent measure of the carbon isotope composition of VPDB, which can calibrate with the IRMS approach to contribute to a better understanding of the isotopic composition of VPDB. However, this method can show some limitations when the overlapping of ¹H-¹²C and ¹H-¹³C signals of a given compound are detected. Furthermore, this method can measure the carbon isotope ratio of a given compound, which has a direct connection between carbon and hydrogen, not for quaternary carbons.

The isotopic analysis of vanillin by ²H-qNMR proved to be efficient to authenticate natural vanillin extracts using the position-specific isotopic ratio of each protonated position of the molecule [17] and enabled the detection and quantification of additions of up to 10% of synthetic or semi-synthetic vanillin in natural vanillin [13]. However, ²H-qNMR requires a long measurement time (5 h acquisition on common spectrometers) and a large sample size (~1 g of vanillin), and is limited to small molecules because of the complexity of the ²H-NMR spectra with a small ²H chemical shift range (~12 ppm at 400 MHz). In addition, hydrogen atoms are easily exchangeable [14]. Furthermore, this method is insufficient to differentiate between the vanillin obtained from ferulic acid by biotechnological processes, and from natural ferulic acid via microorganisms [5].

The analysis of carbon-13 (¹³C-qNMR) presents several advantages over ²H, such as improved sensitivity (natural abundance of ~1.1% and ~0.015%, for ¹³C and ²H, respectively), a larger chemical shift range with fewer overlapping signals, narrower peaks giving a better signal-to-noise

ratio, and robustness against chemical exchange [14, 18]. A disadvantage of measuring ¹³C over ²H, however, is its much smaller range of variations (isotopic variation in natural compounds ranges from ~50 to ~500‰ for ¹³C and ²H, respectively, on the δ scale), which requires an NMR method capable of delivering a precision of a few parts per thousand. In order to attain such a high precision, ¹³C-NMR spectra must be recorded with (i) a high signal-to-noise ratio achieved via extensive signal averaging, (ii) nuclear Overhauser effect (nOe) and partial saturation elimination leading to long experimental times, and (iii) an optimized ¹³C-¹H decoupling sequence in order to perform uniform decoupling over the entire range of ¹H chemical shifts. These requirements were addressed by using (i) inverse-gated decoupled acquisitions with a 90° flip angle [19], (ii) a large inter-pulse delay [20], and (iii) an adiabatic ¹³C-¹H decoupling sequence [21, 22]. The site-specific isotopic compositions of carbon are then calculated from the overall (bulk) ¹³C abundance measured by IRMS and from the peak areas observed by ¹³C-qNMR. The trueness and precision of this single-pulse ¹³C-NMR sequence approach has been demonstrated before [19, 20]. In addition, the reproducibility of this method has been addressed through an inter-laboratory collaboration for which vanillin samples were used as the molecular probe on eight similar spectrometers [23].

In this work, we apply the workflow for measuring carbon isotope ratios in vanillin reference materials by NMR for which we provide certified values for each of the constituent carbon atoms.

<figure><img src="image_4.png" /></figure>
<figure><img src="image_3.png" /><figcaption>**Fig. 1** Principal sources of vanillin [5]</figcaption></figure>

NMR) [16]. This method provides a new and independent measure of the carbon isotope composition of VPDB, which can calibrate with the IRMS approach to contribute to a better understanding of the isotopic composition of VPDB. However, this method can show some limitations when the overlapping of ¹H-¹²C and ¹H-¹³C signals of a given compound are detected. Furthermore, this method can measure the carbon isotope ratio of a given compound, which has a direct connection between carbon and hydrogen, not for quaternary carbons.

The isotopic analysis of vanillin by ²H-qNMR proved to be efficient to authenticate natural vanillin extracts using the position-specific isotopic ratio of each protonated position of the molecule [17] and enabled the detection and quantification of additions of up to 10% of synthetic or semi-synthetic vanillin in natural vanillin [13]. However, ²H-qNMR requires a long measurement time (5 h acquisition on common spectrometers) and a large sample size (~1 g of vanillin), and is limited to small molecules because of the complexity of the ²H-NMR spectra with a small ²H chemical shift range (~12 ppm at 400 MHz). In addition, hydrogen atoms are easily exchangeable [14]. Furthermore, this method is insufficient to differentiate between the vanillin obtained from ferulic acid by biotechnological processes, and from natural ferulic acid via microorganisms [5].

The analysis of carbon-13 (¹³C-qNMR) presents several advantages over ²H, such as improved sensitivity (natural abundance of ~1.1% and ~0.015%, for ¹³C and ²H, respectively), a larger chemical shift range with fewer overlapping signals, narrower peaks giving a better signal-to-noise

ratio, and robustness against chemical exchange [14, 18]. A disadvantage of measuring ¹³C over ²H, however, is its much smaller range of variations (isotopic variation in natural compounds ranges from ~50 to ~500‰ for ¹³C and ²H, respectively, on the δ scale), which requires an NMR method capable of delivering a precision of a few parts per thousand. In order to attain such a high precision, ¹³C-NMR spectra must be recorded with (i) a high signal-to-noise ratio achieved via extensive signal averaging, (ii) nuclear Overhauser effect (nOe) and partial saturation elimination leading to long experimental times, and (iii) an optimized ¹³C-¹H decoupling sequence in order to perform uniform decoupling over the entire range of ¹H chemical shifts. These requirements were addressed by using (i) inverse-gated decoupled acquisitions with a 90° flip angle [19], (ii) a large inter-pulse delay [20], and (iii) an adiabatic ¹³C-¹H decoupling sequence [21, 22]. The site-specific isotopic compositions of carbon are then calculated from the overall (bulk) ¹³C abundance measured by IRMS and from the peak areas observed by ¹³C-qNMR. The trueness and precision of this single-pulse ¹³C-NMR sequence approach has been demonstrated before [19, 20]. In addition, the reproducibility of this method has been addressed through an inter-laboratory collaboration for which vanillin samples were used as the molecular probe on eight similar spectrometers [23].

In this work, we apply the workflow for measuring carbon isotope ratios in vanillin reference materials by NMR for which we provide certified values for each of the constituent carbon atoms.

<figure><img src="image_4.png" /></figure>
# Material and methods

## Chemicals

Acetone-$d_6$, [$^{13}$C$_2$]ethanol and [$^{13}$C$_2$]acetic acid (both materials containing 0.99 mol/mol enrichment of $^{13}$C), and dimethyl sulfoxide-$d_6$ were obtained from Cambridge Isotope Laboratories (Andover, MA, USA). Chromium acetylacetonate (Cr(Acac)$_3$, 99%) and HPLC/Spectro grade acetone were obtained from Sigma-Aldrich (Oakville, ON, Canada). IRMS reference materials, IAEA-CH-6, USGS65, IAEA-600, NBS22, USGS61, IAEA-603, IAEA-610, IAEA-611, and IAEA-612 were obtained from the International Atomic Energy Agency.

NMR tubes (5 mm inner diameter) were obtained from Wilmad LabGlass (Buena, NJ, USA). Two high-purity synthetic vanillin materials were obtained from NRC [24]. Maleic acid (SRM grade, TraceCERT®, having a purity of $0.9989 \pm 0.0031$ g/g) was obtained from Fluka. VANA-1 and VANB-1 are both synthetic high-purity vanillin materials obtained from Fisher Scientific, Waltham, MA, USA) and from Sigma-Aldrich (St. Louis, MO, USA), respectively.

## IRMS analyses

The bulk isotopic composition of carbon in VANA-1 and VANB-1 was determined by elemental analysis combustion IRMS. Specifically, the carbon isotope delta value relative to the VPDB, $\delta(^{13}\text{C})$, was determined at NRC. Details of the standard procedure for $\delta(^{13}\text{C})$ measurements are described by Chartrand et al. [24]. Briefly, approx. 650 µg of vanillin samples and an appropriate amount of carbon isotope reference materials were weighed into 5 × 3.5 mm tin capsules (Elemental Microanalysis; Okehampton, UK) and were loaded onto an elemental analyzer (Vario EL III; Elementar Americas Inc., Mt. Laurel, NJ, USA) interfaced with a gas flow controller (Conflow III; Thermo Fisher; Bremen, Germany) to an isotope ratio mass spectrometer (Delta$^+$XP; Thermo Fisher; Bremen, Germany). Combustion and reduction reactors were set to 950 °C and 500 °C, respectively. Helium dilution on the Conflow III was set to 0.5 bar (7 psi) pressure. Carbon isotope delta measurements were calibrated on the VPDB scale using nine reference materials: IAEA-CH-6, USGS65, IAEA-600, NBS22, USGS61, IAEA-603, IAEA-610, IAEA-611, and IAEA-612. The carbon isotope delta values for VANA-1 and VANB-1 relative to the VPDB have been determined to be $-31.30 \pm 0.06\%$ and $-25.85 \pm 0.05\%$ (expanded uncertainties are at 95% confidence), respectively. Details regarding the assignment of these values are described by Chartrand et al. [25].

## NMR analyses

### $^{1}$H-qNMR

The chemical purity of vanillin in VANA-1 and VANB-1 was determined by quantitative $^{1}$H-NMR on a Bruker-NMR 400 spectrometer (Bruker Optics Inc; MA, USA) operating at 400.13 MHz for proton ($^{1}$H) resonance frequency equipped with a 5 mm ID NMR probe controlled at 298 K, following a method adapted from Le et al. [26]. In brief, three replicate samples of each vanillin sample (15 mg) and maleic acid (8 mg) as an internal calibrator were weighed on an analytical microbalance (UMT-2; Mettler-Toledo) and co-dissolved in 2 mL DMSO-$d_6$ in a 20-mL vial. The solutions were thoroughly vortex-mixed and 0.8 mL of each sample solution was transferred into an NMR tube. The $^{1}$H-qNMR acquisition was performed with a 90° pulse, a sweep width of 15 ppm (6000 Hz), 32 scans of 32,000 data points preceded by 2 dummy scans, an acquisition time of 2.72 s, and a recovery delay of 45 s ($> 7T_{1,max}$) to ensure accuracy of 0.1% with a 90° pulse. Data processing was performed using 128,000 data points. The FIDs were Fourier transformed and apodized by an exponential function leading to a line broadening of 0.3 Hz (LB). The $^{1}$H-qNMR spectra were manually phased and a third-order polynomial baseline correction was applied manually. The signals used for quantitation of vanillin correspond to the following regions of the spectra: 7.70–7.10 ppm (2H, H$_7$ and H$_5$), 7.07–6.70 ppm (1H, H$_3$), 4.10–3.60 ppm (3H, OCH$_3$), and 6.56–5.98 ppm (2H, maleic acid) as shown in Figure S1 (see Electronic Supplementary Material). The reference frequency was set to the DMSO-$d_6$ signal at 2.50 ppm. All spectra were integrated using TopSpin software (v3.6, Bruker) and included the $^{13}$C satellite signals.

### $^{13}$C-qNMR

Quantitative $^{13}$C-NMR spectra were recorded at 100.62 MHz using two Bruker Avance III 400 spectrometers—one in Ottawa (Canada) and the other in Nantes (France)—both equipped with a 5 mm i.d. NMR-BBFO probe. The temperature of the probe was set to 303 K. Probe tuning and matching were performed at the recording frequency of 100.62 MHz. Both $^{13}$C and $^{1}$H channels were automatically tuned and matched, followed by manual verification to check the obtained results.

**Spectrometer adjustments and qualification** The carbon site-specific isotope delta values of the vanillin molecule, $\delta(^{13}\text{C}_k)$ where $k=1...8$ (see Fig. 2 for carbon numbering), were determined using the single-pulse $^{13}$C-NMR sequence mentioned above. Details of the systematic protocol for
# Material and methods

## Chemicals

Acetone-$d_6$, [$^{13}$C$_2$]ethanol and [$^{13}$C$_2$]acetic acid (both materials containing 0.99 mol/mol enrichment of $^{13}$C), and dimethyl sulfoxide-$d_6$ were obtained from Cambridge Isotope Laboratories (Andover, MA, USA). Chromium acetylacetonate (Cr(Acac)$_3$, 99%) and HPLC/Spectro grade acetone were obtained from Sigma-Aldrich (Oakville, ON, Canada). IRMS reference materials, IAEA-CH-6, USGS65, IAEA-600, NBS22, USGS61, IAEA-603, IAEA-610, IAEA-611, and IAEA-612 were obtained from the International Atomic Energy Agency.

NMR tubes (5 mm inner diameter) were obtained from Wilmad LabGlass (Buena, NJ, USA). Two high-purity synthetic vanillin materials were obtained from NRC [24]. Maleic acid (SRM grade, TraceCERT®, having a purity of $0.9989 \pm 0.0031$ g/g) was obtained from Fluka. VANA-1 and VANB-1 are both synthetic high-purity vanillin materials obtained from Fisher Scientific, Waltham, MA, USA) and from Sigma-Aldrich (St. Louis, MO, USA), respectively.

## IRMS analyses

The bulk isotopic composition of carbon in VANA-1 and VANB-1 was determined by elemental analysis combustion IRMS. Specifically, the carbon isotope delta value relative to the VPDB, $\delta(^{13}\text{C})$, was determined at NRC. Details of the standard procedure for $\delta(^{13}\text{C})$ measurements are described by Chartrand et al. [24]. Briefly, approx. 650 µg of vanillin samples and an appropriate amount of carbon isotope reference materials were weighed into 5 × 3.5 mm tin capsules (Elemental Microanalysis; Okehampton, UK) and were loaded onto an elemental analyzer (Vario EL III; Elementar Americas Inc., Mt. Laurel, NJ, USA) interfaced with a gas flow controller (Conflow III; Thermo Fisher; Bremen, Germany) to an isotope ratio mass spectrometer (Delta$^+$XP; Thermo Fisher; Bremen, Germany). Combustion and reduction reactors were set to 950 °C and 500 °C, respectively. Helium dilution on the Conflow III was set to 0.5 bar (7 psi) pressure. Carbon isotope delta measurements were calibrated on the VPDB scale using nine reference materials: IAEA-CH-6, USGS65, IAEA-600, NBS22, USGS61, IAEA-603, IAEA-610, IAEA-611, and IAEA-612. The carbon isotope delta values for VANA-1 and VANB-1 relative to the VPDB have been determined to be $-31.30 \pm 0.06\%$ and $-25.85 \pm 0.05\%$ (expanded uncertainties are at 95% confidence), respectively. Details regarding the assignment of these values are described by Chartrand et al. [25].

## NMR analyses

### $^{1}$H-qNMR

The chemical purity of vanillin in VANA-1 and VANB-1 was determined by quantitative $^{1}$H-NMR on a Bruker-NMR 400 spectrometer (Bruker Optics Inc; MA, USA) operating at 400.13 MHz for proton ($^{1}$H) resonance frequency equipped with a 5 mm ID NMR probe controlled at 298 K, following a method adapted from Le et al. [26]. In brief, three replicate samples of each vanillin sample (15 mg) and maleic acid (8 mg) as an internal calibrator were weighed on an analytical microbalance (UMT-2; Mettler-Toledo) and co-dissolved in 2 mL DMSO-$d_6$ in a 20-mL vial. The solutions were thoroughly vortex-mixed and 0.8 mL of each sample solution was transferred into an NMR tube. The $^{1}$H-qNMR acquisition was performed with a 90° pulse, a sweep width of 15 ppm (6000 Hz), 32 scans of 32,000 data points preceded by 2 dummy scans, an acquisition time of 2.72 s, and a recovery delay of 45 s ($> 7T_{1,max}$) to ensure accuracy of 0.1% with a 90° pulse. Data processing was performed using 128,000 data points. The FIDs were Fourier transformed and apodized by an exponential function leading to a line broadening of 0.3 Hz (LB). The $^{1}$H-qNMR spectra were manually phased and a third-order polynomial baseline correction was applied manually. The signals used for quantitation of vanillin correspond to the following regions of the spectra: 7.70–7.10 ppm (2H, H$_7$ and H$_5$), 7.07–6.70 ppm (1H, H$_3$), 4.10–3.60 ppm (3H, OCH$_3$), and 6.56–5.98 ppm (2H, maleic acid) as shown in Figure S1 (see Electronic Supplementary Material). The reference frequency was set to the DMSO-$d_6$ signal at 2.50 ppm. All spectra were integrated using TopSpin software (v3.6, Bruker) and included the $^{13}$C satellite signals.

### $^{13}$C-qNMR

Quantitative $^{13}$C-NMR spectra were recorded at 100.62 MHz using two Bruker Avance III 400 spectrometers—one in Ottawa (Canada) and the other in Nantes (France)—both equipped with a 5 mm i.d. NMR-BBFO probe. The temperature of the probe was set to 303 K. Probe tuning and matching were performed at the recording frequency of 100.62 MHz. Both $^{13}$C and $^{1}$H channels were automatically tuned and matched, followed by manual verification to check the obtained results.

**Spectrometer adjustments and qualification** The carbon site-specific isotope delta values of the vanillin molecule, $\delta(^{13}\text{C}_k)$ where $k=1...8$ (see Fig. 2 for carbon numbering), were determined using the single-pulse $^{13}$C-NMR sequence mentioned above. Details of the systematic protocol for
<figure><img src="image_5.png" /><figcaption>**Fig. 2** $^{13}$C-qNMR spectrum of vanillin in acetone-d$_{6}$ recorded at NRC on 400 MHz NMR instrument. Numbering of carbon atoms is in the order of decreasing $^{13}$C chemical shift</figcaption></figure>

NMR spectrometer and experimental setup are described by Chaintreau et al. [23].

Determination of the decoupling power attenuation was performed as follows. A series of $^{13}$C spectra of a 10% CHCl$_3$ sample in acetone-d$_{6}$ were acquired using the decp90 pulse program. The transmitter excitation frequencies O1p (for $^{13}$C) and O2p (for $^{1}$H) were set to $\delta(^{13}C) = \delta(^{13}C,CHCl_3) + 0.025$ ppm and $\delta(^{1}H) = \delta(^{1}H,CHCl_3) + 0.1$ ppm. The first $^{13}$C-qNMR spectrum was recorded with an impulsion of 14.1 µs (p3) with the corresponding power level (parameter PL12), delay of 20 s equal to $J/2$ (with $J=212$ Hz), and $^{1}$H power level (PL2=120 dB). The experiment was repeated by decreasing PL2 until the $^{13}$C doublet signal disappears, which corresponds to the 90° $^{1}$H pulse (PL2dec = 13 W). This decoupling attenuation value was used as the decoupling power in the inverse-gated sequence for the adiabatic decoupling.

NMR spectrometer qualification was performed as follows. The adiabatic decoupling pulse and CPD program for adiabatic decoupling were inserted into the $^{13}$C-NMR (zgig) existing pulse sequence in the NMR spectrometer. The ability of the NMR for the proton decoupling conditions was evaluated by using [ $^{13}$C$_2$ ] ethanol. When adiabatic decoupling was used, the acquisition time was limited to 1 s to control the heating effect of decoupling. The $^{13}$C NMR spectrum of the sealed tube was recorded under quantitative conditions of inverse-gated $^{1}$H decoupling and repetition time $d_1=20$ s ($> 10T_{1,max}$) to avoid partial saturation and build-up of the nOe. Thirteen acquisitions of $^{13}$C with $^{1}$H-decoupling were performed with a systematic variation of the $^{1}$H decoupling offset set to $\nu_0 = \frac{1}{2}\nu(^{1}H,CH_3) + \frac{1}{2}\nu(^{1}H,CH_2)$ for ethanol, ranging between $\nu_0 \pm 3000$ Hz with 500 Hz increments with five spectra recorded for each offset value.

Doubly-labelled [ $^{13}$C$_2$ ] ethanol was used to evaluate the NMR performance as described in detail below. In short,

stoichiometry dictates that the peak areas of the two doublets are identical and this was used to evaluate the accuracy of the isotope ratio measurements by $^{13}$C-qNMR [23].

**Measurements of vanillin samples** Approximately 250 mg of vanillin was dissolved in 400 µL of acetone-d$_{6}$ followed by the addition of 100 µL of 0.1 M Cr(Acac)$_3$ (in acetone) relaxation solution using a Hamilton glass syringe. The resulting solution was vortex-mixed and then equilibrated at room temperature for approx. 1 h. The vanillin solution was carefully filtered to remove any non-dissolved particles of Cr(Acac)$_3$ before transferring it to the NMR tube. The NMR tube was sealed with a plug and parafilm to prevent evaporation of the solvent.

The temperature of the probe and sample was set at 303 K. The offsets for both $^{13}$C and $^{1}$H were set at the middle of the frequency range: O1p at 119 ppm (for $^{13}$C) and O2p at 6.5 ppm (for $^{1}$H). Inverse-gated decoupling was applied to avoid the nOe. The longitudinal relaxation time ($T_1$) of $^{13}$C nuclei was determined by using an inversion recovery sequence, with 8 inversion-time values ranging from $\tau=0.1$ s to 20 s, and by using the $T_1$ processing software of the spectrometer. For most vanillin samples, the 90° $^{13}$C pulse widths were calibrated to 9.9 µs. All $T_1$ values of $^{13}$C nuclei were less than 1.6 s. Thus, the repetition delay between each 90° pulse was fixed at 20 s ($> 10T_{1,max}$) to achieve quantitative relaxation of the magnetization and nOe cancellation. The decoupling sequence used adiabatic full-passage pulses with cosine square amplitude modulation ($\nu_{2,max} = 17.6$ kHz) and offset independent adiabaticity with optimized frequency sweep as described elsewhere [21]. The acquisition conditions were as follows: acquisition time of 0.7 s, repetition delays of 20 s, 4 dummy scans, and 400 sample scans, resulting in the overall acquisition time of 140 min per spectrum. Five spectra were consecutively acquired for each sample. A typical $^{13}$C-qNMR spectrum of vanillin is shown in Fig. 2.

Samples of VANA-1 and VANB-1 were measured in Ottawa and additional samples of VANA-1 and VANB-1 were also measured in Nantes with the same $^{13}$C-qNMR parameters except for the $^{13}$C-pw90 for where the value optimized for vanillin on the Nantes instrument was used.

**NMR data processing**

Data were processed using the TopSpin software (v3.6) with the same parameters and the same procedure for free induction decays of 256 K: (i) an exponential multiplication inducing a 2 Hz line broadening, (ii) manual phasing, and (iii) automatic baseline correction using a third-order polynomial function. The carbon signal peak areas were calculated in accordance with a frequency domain mathematical model
<figure><img src="image_5.png" /><figcaption>**Fig. 2** $^{13}$C-qNMR spectrum of vanillin in acetone-d$_{6}$ recorded at NRC on 400 MHz NMR instrument. Numbering of carbon atoms is in the order of decreasing $^{13}$C chemical shift</figcaption></figure>

NMR spectrometer and experimental setup are described by Chaintreau et al. [23].

Determination of the decoupling power attenuation was performed as follows. A series of $^{13}$C spectra of a 10% CHCl$_3$ sample in acetone-d$_6$ were acquired using the decp90 pulse program. The transmitter excitation frequencies O1p (for $^{13}$C) and O2p (for $^{1}$H) were set to $\delta(^{13}C) = \delta(^{13}C,CHCl_3) + 0.025$ ppm and $\delta(^{1}H) = \delta(^{1}H,CHCl_3) + 0.1$ ppm. The first $^{13}$C-qNMR spectrum was recorded with an impulsion of 14.1 µs (p3) with the corresponding power level (parameter PL12), delay of 20 s equal to $J/2$ (with $J=212$ Hz), and $^{1}$H power level (PL2=120 dB). The experiment was repeated by decreasing PL2 until the $^{13}$C doublet signal disappears, which corresponds to the 90° $^{1}$H pulse (PL2dec = 13 W). This decoupling attenuation value was used as the decoupling power in the inverse-gated sequence for the adiabatic decoupling.

NMR spectrometer qualification was performed as follows. The adiabatic decoupling pulse and CPD program for adiabatic decoupling were inserted into the $^{13}$C-NMR (zgig) existing pulse sequence in the NMR spectrometer. The ability of the NMR for the proton decoupling conditions was evaluated by using [ $^{13}$C$_2$ ] ethanol. When adiabatic decoupling was used, the acquisition time was limited to 1 s to control the heating effect of decoupling. The $^{13}$C NMR spectrum of the sealed tube was recorded under quantitative conditions of inverse-gated $^{1}$H decoupling and repetition time $d_1=20$ s ($> 10T_{1,max}$) to avoid partial saturation and build-up of the nOe. Thirteen acquisitions of $^{13}$C with $^{1}$H-decoupling were performed with a systematic variation of the $^{1}$H decoupling offset set to $\nu_0 = \frac{1}{2}\nu(^{1}H,CH_3) + \frac{1}{2}\nu(^{1}H,CH_2)$ for ethanol, ranging between $\nu_0 \pm 3000$ Hz with 500 Hz increments with five spectra recorded for each offset value.

Doubly-labelled [ $^{13}$C$_2$ ] ethanol was used to evaluate the NMR performance as described in detail below. In short,

stoichiometry dictates that the peak areas of the two doublets are identical and this was used to evaluate the accuracy of the isotope ratio measurements by $^{13}$C-qNMR [23].

**Measurements of vanillin samples** Approximately 250 mg of vanillin was dissolved in 400 µL of acetone-d$_6$ followed by the addition of 100 µL of 0.1 M Cr(Acac)$_3$ (in acetone) relaxation solution using a Hamilton glass syringe. The resulting solution was vortex-mixed and then equilibrated at room temperature for approx. 1 h. The vanillin solution was carefully filtered to remove any non-dissolved particles of Cr(Acac)$_3$ before transferring it to the NMR tube. The NMR tube was sealed with a plug and parafilm to prevent evaporation of the solvent.

The temperature of the probe and sample was set at 303 K. The offsets for both $^{13}$C and $^{1}$H were set at the middle of the frequency range: O1p at 119 ppm (for $^{13}$C) and O2p at 6.5 ppm (for $^{1}$H). Inverse-gated decoupling was applied to avoid the nOe. The longitudinal relaxation time ($T_1$) of $^{13}$C nuclei was determined by using an inversion recovery sequence, with 8 inversion-time values ranging from $\tau=0.1$ s to 20 s, and by using the $T_1$ processing software of the spectrometer. For most vanillin samples, the 90° $^{13}$C pulse widths were calibrated to 9.9 µs. All $T_1$ values of $^{13}$C nuclei were less than 1.6 s. Thus, the repetition delay between each 90° pulse was fixed at 20 s ($> 10T_{1,max}$) to achieve quantitative relaxation of the magnetization and nOe cancellation. The decoupling sequence used adiabatic full-passage pulses with cosine square amplitude modulation ($\nu_{2,max} = 17.6$ kHz) and offset independent adiabaticity with optimized frequency sweep as described elsewhere [21]. The acquisition conditions were as follows: acquisition time of 0.7 s, repetition delays of 20 s, 4 dummy scans, and 400 sample scans, resulting in the overall acquisition time of 140 min per spectrum. Five spectra were consecutively acquired for each sample. A typical $^{13}$C-qNMR spectrum of vanillin is shown in Fig. 2.

Samples of VANA-1 and VANB-1 were measured in Ottawa and additional samples of VANA-1 and VANB-1 were also measured in Nantes with the same $^{13}$C-qNMR parameters except for the $^{13}$C-pw90 for where the value optimized for vanillin on the Nantes instrument was used.

**NMR data processing**

Data were processed using the TopSpin software (v3.6) with the same parameters and the same procedure for free induction decays of 256 K: (i) an exponential multiplication inducing a 2 Hz line broadening, (ii) manual phasing, and (iii) automatic baseline correction using a third-order polynomial function. The carbon signal peak areas were calculated in accordance with a frequency domain mathematical model
using Rnmfit package in R [27] using the following integration windows: 191.7–189.7 ppm (C1); 153.7–151.7 ppm (C2); 149.0–147.0 ppm (C3); 130.8–128.8 ppm (C4); 127.4–125.4 ppm (C5); 116.2–114.2 ppm (C6); 111.2–109.2 ppm (C7); and 56.6–54.4 ppm (C8).

## Site-specific carbon isotope delta values

The carbon isotope specific-position for each vanillin was established and traced based on the bulk carbon isotopic composition obtained by IRMS. For subsequent model calculations, we have adopted the ¹³C/¹²C isotope ratio R<sub>VPDB</sub>=0.011108 for the VPDB which is the weighted average of three recent independent estimates [16, 28, 29].

Carbon isotope ratios for each carbon atom were calculated from the processed spectra using the measurement model of Silvestre et al. [30] as outlined below. The input data are the areas corresponding to the main (central) peaks observed in ¹³C-NMR spectra. The established procedure is to adjust these measured peak areas (S) to account for the slight loss of intensity caused by the satellite peaks arising from the ¹³C–¹³C interactions [23]:

$$S_{i,\text{adj}} = S_i(1 + 0.011n_i) \qquad (1)$$

where n<sub>i</sub> is the number of carbon atoms directly bonded to the carbon atom corresponding to each NMR signal. For example, C1 (i=1) is directly bound to C4 therefore n<sub>1</sub>=1 whereas C4 is directly bound to C1, C5, and C7 therefore n<sub>4</sub>=3.

The isotopic distribution of vanillin is first characterized by the relative share of ¹³C for each carbon atom:

$$f_i(^{13}\text{C}) = S_{i,\text{adj}} / (S_{1,\text{adj}} + S_{2,\text{adj}} + \dots + S_{8,\text{adj}}) \qquad (2)$$

The relative share of ¹³C when all eight atoms have identical isotopic composition is:

$$F_i(^{13}\text{C}) = P_i / (P_1 + P_2 + \dots + P_8) \qquad (3)$$

where P<sub>i</sub> is the number of equivalent carbon atoms for the molecular position k. Each peak observed in the ¹³C-NMR spectrum of vanillin presents a single equivalent carbon (Fig. 2); therefore P<sub>i</sub>=1 for i=1...8 and F<sub>i</sub>=1/8.

The average abundance of ¹³C in the vanillin sample, x(¹³C), is obtained from δ(¹³C) measurements by IRMS and the ¹³C/¹²C reference value for VPDB:

$$R_{13} = R_{\text{VPDB}} [1 + \delta(^{13}\text{C})] \qquad (4)$$

$$x(^{13}\text{C}) = R_{13} / (1 + R_{13}) \qquad (5)$$

This average value is used to obtain the isotopic abundance of carbon-13 at each position k:

$$x_i(^{13}\text{C}) = x(^{13}\text{C}) \cdot f_i / F_i \qquad (6)$$

Similarly, the ¹³C/¹²C isotope ratio is calculated for each carbon atom:

$$R_{13,i} = x_i(^{13}\text{C}) / [1 - x_i(^{13}\text{C})] \qquad (7)$$

Lastly, one obtains carbon isotope delta values for each carbon atom in vanillin:

$$\delta_i(^{13}\text{C}) = R_{13,i} / R_{\text{VPDB}} - 1 \qquad (8)$$

In order to further investigate the above measurement model, we note that Eqs. 1–8 (which we term model M0) assume an identical abundance of ¹³C (1.1%) for all carbon atoms in the clumped isotopologues having two or more C<sub>13</sub> atoms, instead of using the actual abundances to perform the correction of signal areas. To overcome this simplification, one can identify the specific isotopologues that make up the central (main) peaks integrated into this study. Figure 3 shows such a schematic for the C1 signal. Taken these considerations together, one can establish a more refined measurement model for vanillin in Eqs. 9–16.

The measurement model for site-specific carbon isotope delta values in vanillin can be written in the form of a generative model which lends itself more suitable for Bayesian analysis (see Supplementary file), with the observed peak areas for vanillin carbon atoms, S<sub>i</sub> (i=1...8), modelled as an outcome of (unknown) NMR sensitivity factor, s, and (unknown) isotopic abundances of carbon-13 at each site, x<sub>13,i</sub>:

$$S_1 = s(a_1 + a_{1,S} - a_{1,4}) \qquad (9)$$

$$S_2 = s(a_2 + a_{2,S} - a_{2,3} - a_{2,6}) \qquad (10)$$

$$S_3 = s(a_3 + a_{3,S} - a_{3,2} - a_{3,7}) \qquad (11)$$

$$S_4 = s(a_4 + a_{4,S} - a_{4,1} - a_{4,5} - a_{4,7}) \qquad (12)$$

$$S_5 = s(a_5 + a_{5,S} - a_{5,4} - a_{5,6}) \qquad (13)$$

$$S_6 = s(a_6 + a_{6,S} - a_{6,2} - a_{6,5}) \qquad (14)$$

$$S_7 = s(a_7 + a_{7,S} - a_{7,3} - a_{7,4}) \qquad (15)$$

$$S_8 = s(a_8 + a_{8,S}) \qquad (16)$$

Parameters a<sub>i</sub> and a<sub>i,j</sub> (i,j=1...8) are as described in Fig. 3 and a<sub>i,S</sub>=sum{a<sub>i,k</sub>} (k=1...8). These equations are solved for s and x<sub>13,i</sub> (i=1...8) using numerical methods along with the constraint imposed on the values of x<sub>13,i</sub> that the individual
using Rnmfit package in R [27] using the following integration windows: 191.7–189.7 ppm (C1); 153.7–151.7 ppm (C2); 149.0–147.0 ppm (C3); 130.8–128.8 ppm (C4); 127.4–125.4 ppm (C5); 116.2–114.2 ppm (C6); 111.2–109.2 ppm (C7); and 56.6–54.4 ppm (C8).

## Site-specific carbon isotope delta values

The carbon isotope specific-position for each vanillin was established and traced based on the bulk carbon isotopic composition obtained by IRMS. For subsequent model calculations, we have adopted the ¹³C/¹²C isotope ratio R<sub>VPDB</sub>=0.011108 for the VPDB which is the weighted average of three recent independent estimates [16, 28, 29].

Carbon isotope ratios for each carbon atom were calculated from the processed spectra using the measurement model of Silvestre et al. [30] as outlined below. The input data are the areas corresponding to the main (central) peaks observed in ¹³C-NMR spectra. The established procedure is to adjust these measured peak areas (S) to account for the slight loss of intensity caused by the satellite peaks arising from the ¹³C–¹³C interactions [23]:

$$S_{i,\text{adj}} = S_i(1 + 0.011n_i) \quad (1)$$

where n<sub>i</sub> is the number of carbon atoms directly bonded to the carbon atom corresponding to each NMR signal. For example, C1 (i=1) is directly bound to C4 therefore n<sub>1</sub>=1 whereas C4 is directly bound to C1, C5, and C7 therefore n<sub>4</sub>=3.

The isotopic distribution of vanillin is first characterized by the relative share of ¹³C for each carbon atom:

$$f_i(^{13}\text{C}) = S_{i,\text{adj}} / (S_{1,\text{adj}} + S_{2,\text{adj}} + \dots + S_{8,\text{adj}}) \quad (2)$$

The relative share of ¹³C when all eight atoms have identical isotopic composition is:

$$F_i(^{13}\text{C}) = P_i / (P_1 + P_2 + \dots + P_8) \quad (3)$$

where P<sub>i</sub> is the number of equivalent carbon atoms for the molecular position k. Each peak observed in the ¹³C-NMR spectrum of vanillin presents a single equivalent carbon (Fig. 2); therefore P<sub>i</sub>=1 for i=1...8 and F<sub>i</sub>=1/8.

The average abundance of ¹³C in the vanillin sample, x(¹³C), is obtained from δ(¹³C) measurements by IRMS and the ¹³C/¹²C reference value for VPDB:

$$R_{13} = R_{\text{VPDB}} [1 + \delta(^{13}\text{C})] \quad (4)$$

$$x(^{13}\text{C}) = R_{13} / (1 + R_{13}) \quad (5)$$

This average value is used to obtain the isotopic abundance of carbon-13 at each position k:

$$x_i(^{13}\text{C}) = x(^{13}\text{C}) \cdot f_i / F_i \quad (6)$$

Similarly, the ¹³C/¹²C isotope ratio is calculated for each carbon atom:

$$R_{13,i} = x_i(^{13}\text{C}) / [1 - x_i(^{13}\text{C})] \quad (7)$$

Lastly, one obtains carbon isotope delta values for each carbon atom in vanillin:

$$\delta_i(^{13}\text{C}) = R_{13,i} / R_{\text{VPDB}} - 1 \quad (8)$$

In order to further investigate the above measurement model, we note that Eqs. 1–8 (which we term model M0) assume an identical abundance of ¹³C (1.1%) for all carbon atoms in the clumped isotopologues having two or more C<sub>13</sub> atoms, instead of using the actual abundances to perform the correction of signal areas. To overcome this simplification, one can identify the specific isotopologues that make up the central (main) peaks integrated into this study. Figure 3 shows such a schematic for the C1 signal. Taken these considerations together, one can establish a more refined measurement model for vanillin in Eqs. 9–16.

The measurement model for site-specific carbon isotope delta values in vanillin can be written in the form of a generative model which lends itself more suitable for Bayesian analysis (see Supplementary file), with the observed peak areas for vanillin carbon atoms, S<sub>i</sub> (i=1...8), modelled as an outcome of (unknown) NMR sensitivity factor, s, and (unknown) isotopic abundances of carbon-13 at each site, x<sub>13,i</sub>:

$$S_1 = s(a_1 + a_{1,S} - a_{1,4}) \quad (9)$$

$$S_2 = s(a_2 + a_{2,S} - a_{2,3} - a_{2,6}) \quad (10)$$

$$S_3 = s(a_3 + a_{3,S} - a_{3,2} - a_{3,7}) \quad (11)$$

$$S_4 = s(a_4 + a_{4,S} - a_{4,1} - a_{4,5} - a_{4,7}) \quad (12)$$

$$S_5 = s(a_5 + a_{5,S} - a_{5,4} - a_{5,6}) \quad (13)$$

$$S_6 = s(a_6 + a_{6,S} - a_{6,2} - a_{6,5}) \quad (14)$$

$$S_7 = s(a_7 + a_{7,S} - a_{7,3} - a_{7,4}) \quad (15)$$

$$S_8 = s(a_8 + a_{8,S}) \quad (16)$$

Parameters a<sub>i</sub> and a<sub>i,j</sub> (i,j=1...8) are as described in Fig. 3 and a<sub>i,S</sub>=sum{a<sub>i,k</sub>} (k=1...8). These equations are solved for s and x<sub>13,i</sub> (i=1...8) using numerical methods along with the constraint imposed on the values of x<sub>13,i</sub> that the individual
<figure><img src="image_6.png" /><figcaption>**Fig. 3** Schematic components of the C1 carbon atom peak in the ¹³C-qNMR spectrum of vanillin. Note that only the ¹³C–¹³C interactions are observed in this acquisition mode. The observed peak area for C1 is proportional to the sum of these components with the exception of ¹³C1¹³C4 isotopologue ($a_{1,4}$) which is outside the integration area of the main peak</figcaption></figure>

isotope values, $\delta_i(^{13}\text{C})$, calculated using Eqs. 7–8, must average to the value provided by IRMS. In the simplest interpretation, Eqs. 9–16 are parametrized in terms of eight unknown ¹³C abundances, $x_{13,i}$ ($i=1...8$). This assumes, for example, that the isotopic abundance of ¹³C in the C1 atom is the same in all isotopologues containing the C1 carbon atom shown in Fig. 3. We term this model M1. This nonlinear measurement model can be solved using standard numerical solver procedures available in either Excel (see Electronic Supplementary Material) or with dedicated mathematical software such as Mathematica or R. Our results show that the sensitivity factor ($s$) varies by approx. 1–3 parts per thousand between the replicate measurements of vanillin samples. It is simply a scaling factor between the (relative) signals predicted by the model and the observed peak areas.

In reality, there might be some variation in the isotopic abundances of ¹³C (in the carbon atoms of the same position in the vanillin molecule) among the various isotopologues. This variability can be addressed using hierarchical random effects models, similar to how between-unit homogeneity is commonly evaluated in the reference materials. This, however, is outside the scope of this work.

# Results and discussion

## ¹H-qNMR for purity determination

Determination of the site-specific carbon isotope profile by ¹³C-qNMR requires high-purity (> 99%) organic materials [31]. In many cases, this can be achieved by the use

of commercial chemicals. We used the ¹H-qNMR technique with maleic acid qNMR internal standard to evaluate the chemical purity of vanillin in both VANA-1 and VANB-1 [32]. The integration regions in the ¹H-qNMR measurements are wide enough to contain the adjacent ¹³C satellites and the integrated signal area: for vanillin material 7.70–7.10 ppm: 2H ($H_5$ and $H_7$); 7.07–6.70 ppm: 1H ($H_3$); and 4.10–3.60: 3H–OCH₃; for maleic acid as IS 3H: 6.56–5.98 ppm (as described in “¹H-qNMR”).

The purity of vanillin in VANA-1 and VANB-1 materials of 0.9982 g/g and 0.9973 g/g, respectively determined using ¹H-qNMR, is thus suitable to use for the determination of the carbon isotope profile of these materials. The assignment of these values is described in Electronic Supplementary Material Table S1.

## ¹³C-NMR pulse sequence for carbon isotope site-specific measurement

The main difficulty with obtaining high-precision isotope ratio values with ¹³C-NMR lies with ensuring uniform ¹H decoupling efficiency over the whole range of ¹H chemical shifts. This problem has been solved with the help of an adiabatic ¹H decoupling sequence using a numerical calculation of the frequency sweep from an amplitude evolution with a cosine shape [21]. The performance qualification of the NMR spectrometer was evaluated using doubly-labelled [¹³C₂]ethanol as a molecular probe to demonstrate uniform and homogeneous ¹H decoupling efficiency over a large range of resonances [23]. The ¹³C-qNMR spectrum of [¹³C₂]ethanol shows four main peaks (Fig. 4)—a doublet for ¹³CH₂
<figure><img src="image_6.png" /><figcaption>**Fig. 3** Schematic components of the C1 carbon atom peak in the ¹³C-qNMR spectrum of vanillin. Note that only the ¹³C–¹³C interactions are observed in this acquisition mode. The observed peak area for C1 is proportional to the sum of these components with the exception of ¹³C1¹³C4 isotopologue ($a_{1,4}$) which is outside the integration area of the main peak</figcaption></figure>

isotope values, $\delta_i(^{13}\text{C})$, calculated using Eqs. 7–8, must average to the value provided by IRMS. In the simplest interpretation, Eqs. 9–16 are parametrized in terms of eight unknown ¹³C abundances, $x_{13,i}$ ($i=1...8$). This assumes, for example, that the isotopic abundance of ¹³C in the C1 atom is the same in all isotopologues containing the C1 carbon atom shown in Fig. 3. We term this model M1. This nonlinear measurement model can be solved using standard numerical solver procedures available in either Excel (see Electronic Supplementary Material) or with dedicated mathematical software such as Mathematica or R. Our results show that the sensitivity factor ($s$) varies by approx. 1–3 parts per thousand between the replicate measurements of vanillin samples. It is simply a scaling factor between the (relative) signals predicted by the model and the observed peak areas.

In reality, there might be some variation in the isotopic abundances of ¹³C (in the carbon atoms of the same position in the vanillin molecule) among the various isotopologues. This variability can be addressed using hierarchical random effects models, similar to how between-unit homogeneity is commonly evaluated in the reference materials. This, however, is outside the scope of this work.

# Results and discussion

## ¹H-qNMR for purity determination

Determination of the site-specific carbon isotope profile by ¹³C-qNMR requires high-purity (> 99%) organic materials [31]. In many cases, this can be achieved by the use

of commercial chemicals. We used the ¹H-qNMR technique with maleic acid qNMR internal standard to evaluate the chemical purity of vanillin in both VANA-1 and VANB-1 [32]. The integration regions in the ¹H-qNMR measurements are wide enough to contain the adjacent ¹³C satellites and the integrated signal area: for vanillin material 7.70–7.10 ppm: 2H ($H_5$ and $H_7$); 7.07–6.70 ppm: 1H ($H_3$); and 4.10–3.60: 3H–OCH₃; for maleic acid as IS 3H: 6.56–5.98 ppm (as described in “¹H-qNMR”).

The purity of vanillin in VANA-1 and VANB-1 materials of 0.9982 g/g and 0.9973 g/g, respectively determined using ¹H-qNMR, is thus suitable to use for the determination of the carbon isotope profile of these materials. The assignment of these values is described in Electronic Supplementary Material Table S1.

## ¹³C-NMR pulse sequence for carbon isotope site-specific measurement

The main difficulty with obtaining high-precision isotope ratio values with ¹³C-NMR lies with ensuring uniform ¹H decoupling efficiency over the whole range of ¹H chemical shifts. This problem has been solved with the help of an adiabatic ¹H decoupling sequence using a numerical calculation of the frequency sweep from an amplitude evolution with a cosine shape [21]. The performance qualification of the NMR spectrometer was evaluated using doubly-labelled [¹³C₂]ethanol as a molecular probe to demonstrate uniform and homogeneous ¹H decoupling efficiency over a large range of resonances [23]. The ¹³C-qNMR spectrum of [¹³C₂]ethanol shows four main peaks (Fig. 4)—a doublet for ¹³CH₂
and a doublet for ¹³CH₃ resonances. The expected area of these two doublets is S<sub>CH2</sub>/S<sub>CH3</sub> = 1, irrespective of the precise isotopic composition of either of the carbon atoms, and this was used to evaluate the accuracy of ¹³C-qNMR. Any deviation from this expected ratio would indicate a lack of accuracy in the NMR configuration used and for our purposes this deviation should be within few parts per thousand (1–2‰) [23].

When applied to the 400 MHz NMR at NRC, this instru-
ment qualification procedure showed a deviation of ≤2‰
over the wide range of ¹H excitation frequencies (Fig. 4),
thus indicating that our NMR spectrometer is suitably quali-
fied to measure ¹³C/¹²C isotope ratios at natural abundances.
We also compared the adiabatic ¹³C-¹H-decoupling scheme
to the WALTZ-16 standard pulse on doubly-labelled ethanol
(Fig. 4). The latter pulse sequence shows significant biases
indicating that the WALTZ-16 ¹³C-¹H-decoupling is sensi-
tive to off-resonance effects and to the presence of ¹H-¹H
couplings and is not suitable for high-precision carbon iso-
tope ratio measurements by NMR [21]. To obtain an esti-
mate of how reliable the peak area measurements are for
vanillin, a similar qualification procedure was performed
using doubly-labelled [¹³C₂]-acetic acid. The average devia-
tion from the 1:1 ratio of the two doublets (–COOH and
–CH₃) is approx. 2‰, identical to the ethanol.

Isotopic composition of vanillin by ¹³C-qNMR

Isotopic.¹³C NMR experiments applied on VANA-1
and VANB-1

Isotope ratio precision of ¹³C-qNMR measurements
depends mainly on the signal-to-noise ratio achievable in
a reasonable period of time which, in turn, is a function
of the acquisition time, the delay between pulses, and the

number of scans [14]. The vanillin molecule presents sev-
eral challenging technical features: a wide range of chemi-
cal shifts for carbon atoms (135 ppm) and hydrogen atoms
(6 ppm), long relaxation times (up to 19 s), and a wide
range of nOe factors. These issues were resolved with the
help of an optimized adiabatic ¹H decoupling sequence [19,
21] and 0.1 M Cr(Acac)₃ paramagnetic relaxation solution.
At this concentration, Cr(Acac)₃ does not affect the sensi-
tivity and accuracy of ¹³C-qNMR [20], while a large reduc-
tion of the T₁ values (from 19 to 1.6 s) for each carbon in
vanillin is observed upon the addition of Cr(Acac)₃. This
provides a significant reduction of the inter-pulse delay
thus allowing a tenfold reduction of the overall acquisi-
tion time [18, 20]. In our study, the longest T₁ of 1.57 s
was observed at C4 carbon in vanillin (Table 1). In such
conditions, the achieved peak area precision of 1–2 parts
per thousand was achieved for the determination of the
position-specific isotopic analysis of VANA-1 and VANB-
1. Setting the acquisition time up to 0.7 s for adiabatic
decoupling ensured probe integrity as well as avoided
the heating effect of decoupling and nOe. The use of 400
scans for 250 mg of vanillin solution in 500 µL acetone-d₆
with 0.1 M Cr(Acac)₃ (in acetone) solution resulted in the
signal-to-noise ratio S:N > 720 for all carbon peaks in the
vanillin samples.

Site-specific ¹³C content in vanillin samples

Using the optimized experimental conditions described
above, we have performed site-specific carbon isotope
measurements in reference materials VANA-1 and VANB-
1. First, we discuss the effect of the measurement model. We
note that the measurements of a single vanillin sample can
provide carbon isotope delta values that differ by up to 1‰
between the two measurement models outlined in this work

<figure><img src="image_7.png" /><figcaption>**Fig. 4** NRC 400 MHz NMR qualification using adiabatic and standard WALTZ-16 ¹³C-¹H decoupling pulses on the doubly-labelled [¹³C₂]-ethanol</figcaption></figure>
and a doublet for ¹³CH₃ resonances. The expected area of these two doublets is S<sub>CH2</sub>/S<sub>CH3</sub> = 1, irrespective of the precise isotopic composition of either of the carbon atoms, and this was used to evaluate the accuracy of ¹³C-qNMR. Any deviation from this expected ratio would indicate a lack of accuracy in the NMR configuration used and for our purposes this deviation should be within few parts per thousand (1–2‰) [23].

When applied to the 400 MHz NMR at NRC, this instru-
ment qualification procedure showed a deviation of ≤2‰
over the wide range of ¹H excitation frequencies (Fig. 4),
thus indicating that our NMR spectrometer is suitably quali-
fied to measure ¹³C/¹²C isotope ratios at natural abundances.
We also compared the adiabatic ¹³C-¹H-decoupling scheme
to the WALTZ-16 standard pulse on doubly-labelled ethanol
(Fig. 4). The latter pulse sequence shows significant biases
indicating that the WALTZ-16 ¹³C-¹H-decoupling is sensi-
tive to off-resonance effects and to the presence of ¹H-¹H
couplings and is not suitable for high-precision carbon iso-
tope ratio measurements by NMR [21]. To obtain an esti-
mate of how reliable the peak area measurements are for
vanillin, a similar qualification procedure was performed
using doubly-labelled [¹³C₂]-acetic acid. The average devia-
tion from the 1:1 ratio of the two doublets (–COOH and
–CH₃) is approx. 2‰, identical to the ethanol.

Isotopic composition of vanillin by ¹³C-qNMR

Isotopic.¹³C NMR experiments applied on VANA-1
and VANB-1

Isotope ratio precision of ¹³C-qNMR measurements
depends mainly on the signal-to-noise ratio achievable in
a reasonable period of time which, in turn, is a function
of the acquisition time, the delay between pulses, and the

number of scans [14]. The vanillin molecule presents sev-
eral challenging technical features: a wide range of chemi-
cal shifts for carbon atoms (135 ppm) and hydrogen atoms
(6 ppm), long relaxation times (up to 19 s), and a wide
range of nOe factors. These issues were resolved with the
help of an optimized adiabatic ¹H decoupling sequence [19,
21] and 0.1 M Cr(Acac)₃ paramagnetic relaxation solution.
At this concentration, Cr(Acac)₃ does not affect the sensi-
tivity and accuracy of ¹³C-qNMR [20], while a large reduc-
tion of the T₁ values (from 19 to 1.6 s) for each carbon in
vanillin is observed upon the addition of Cr(Acac)₃. This
provides a significant reduction of the inter-pulse delay
thus allowing a tenfold reduction of the overall acquisi-
tion time [18, 20]. In our study, the longest T₁ of 1.57 s
was observed at C4 carbon in vanillin (Table 1). In such
conditions, the achieved peak area precision of 1–2 parts
per thousand was achieved for the determination of the
position-specific isotopic analysis of VANA-1 and VANB-
1. Setting the acquisition time up to 0.7 s for adiabatic
decoupling ensured probe integrity as well as avoided
the heating effect of decoupling and nOe. The use of 400
scans for 250 mg of vanillin solution in 500 µL acetone-*d*₆
with 0.1 M Cr(Acac)₃ (in acetone) solution resulted in the
signal-to-noise ratio S:N > 720 for all carbon peaks in the
vanillin samples.

Site-specific ¹³C content in vanillin samples

Using the optimized experimental conditions described
above, we have performed site-specific carbon isotope
measurements in reference materials VANA-1 and VANB-
1. First, we discuss the effect of the measurement model. We
note that the measurements of a single vanillin sample can
provide carbon isotope delta values that differ by up to 1‰
between the two measurement models outlined in this work

<figure><img src="image_7.png" /><figcaption>**Fig. 4** NRC 400 MHz NMR qualification using adiabatic and standard WALTZ-16 ¹³C-¹H decoupling pulses on the doubly-labelled [¹³C₂]-ethanol</figcaption></figure>
(Table 2). In this work, we adopt the arithmetic average of these two models.

The site-specific carbon isotope delta values, δ<sub>t</sub>(<sup>13</sup>C), measured at the NRC by <sup>13</sup>C-qNMR for VANA-1 and VANB-1 samples, are presented in Tables 3–4. The δ<sub>t</sub>(<sup>13</sup>C) values are expressed in permilles relative to the VPDB, five measurements for each bottle. The consistency of δ<sub>t</sub>(<sup>13</sup>C) values from sample to sample range from 0.3 to 1.1‰ for both VANA-1 and VANB-1, as evaluated by the repeatability standard deviations from NRC measurements shown in Tables 3–4.

VANA-1 and VANB-1, have average δ(<sup>13</sup>C) values of –31‰ and –26‰ relative to the VPDB, respectively, consistent with synthetic origin. Vanillin from petrochemical origin typically has δ(<sup>13</sup>C) values from –36 to –25‰ and lignin-based biovanillin has δ(<sup>13</sup>C) values from –29 to –26‰ [2, 5]. Thus, these can be distinguished from natural vanillin from the vanilla pod having δ(<sup>13</sup>C) values from –22 to –14‰ [11]. Although the average isotopic composition of vanillin from VANA-1 and VANB-1 are somewhat similar, the difference in the isotopic composition

of carbon C8 (OCH<sub>3</sub>) between these two materials is extremely large, –54‰ and –31‰, which likely reflects the different synthetic routes for these vanillin samples. In addition, we observe substantial variations in the isotopic composition of carbon atoms in the benzene ring of vanillin for both VANA-1 and VANB-1. Both materials exhibit approximately 10‰ difference between the individual carbon atoms—well in excess of the repeatability standard deviation (<1‰). In both materials, C7 carbon atoms show the highest abundance of <sup>13</sup>C and the C2 atoms the lowest.

## Inter-laboratory comparison

To support the NRC certification work, an inter-laboratory comparison study was conducted with Nantes laboratory (France), which also used an adiabatic pulse sequence for <sup>1</sup>H decoupling on the 400 MHz spectrometer. To avoid variability associated with the NMR data processing, the free induction decays from vanillin samples

<table><caption>**Table 1** Longitudinal relaxation times of each carbon in vanillin in the presence of 0.1 M Cr(Acac)<sub>3</sub><sup>a</sup></caption><thead><tr><th></th><th>C1</th><th>C2</th><th>C3</th><th>C4</th><th>C5</th><th>C6</th><th>C7</th><th>C8</th></tr></thead><tbody><tr><th>T<sub>1</sub>/s</th><td>1.03 (0.07)</td><td>1.26 (0.09)</td><td>1.46 (0.08)</td><td>1.57 (0.05)</td><td>0.85 (0.04)</td><td>0.84 (0.05)</td><td>0.93 (0.03)</td><td>0.86<br>(0.06)</td></tr></tbody></table>

<table><caption>**Table 2** The effect of measurement model on the carbon isotope delta results for a single measurement of vanillin in VANA-1 sample from vial number 170 at NRC</caption><thead><tr><th>Model</th><th>C1</th><th>C2</th><th>C3</th><th>C4</th><th>C5</th><th>C6</th><th>C7</th><th>C8</th></tr></thead><tbody><tr><th>M0</th><td>–20.8</td><td>–30.7</td><td>–33.4</td><td>–28.7</td><td>–28.0</td><td>–30.0</td><td>–24.0</td><td>–54.8</td></tr><tr><th>M1</th><td>–20.4</td><td>–31.0</td><td>–33.6</td><td>–29.2</td><td>–28.2</td><td>–30.2</td><td>–24.3</td><td>–53.6</td></tr></tbody></table>

<table><caption>**Table 3** Site-specific average carbon isotope delta values of vanillin in VANA-1 samples measured at NRC</caption><thead><tr><th></th><th>VANA-1 vial number</th><th>C1</th><th>C2</th><th>C3</th><th>C4</th><th>C5</th><th>C6</th><th>C7</th><th>C8</th></tr></thead><tbody><tr><th>120</th><td></td><td>–22.1</td><td>–31.8</td><td>–32.6</td><td>–28.4</td><td>–29.6</td><td>–29.9</td><td>–23.2</td><td>–52.8</td></tr><tr><th>170</th><td></td><td>–21.1</td><td>–32.1</td><td>–33.1</td><td>–29.2</td><td>–29.4</td><td>–29.3</td><td>–22.8</td><td>–53.5</td></tr><tr><th>290</th><td></td><td>–21.1</td><td>–32.6</td><td>–31.8</td><td>–27.9</td><td>–30.6</td><td>–30.0</td><td>–23.0</td><td>–53.4</td></tr><tr><th>380</th><td></td><td>–21.8</td><td>–32.4</td><td>–31.4</td><td>–27.0</td><td>–30.4</td><td>–29.8</td><td>–24.2</td><td>–53.3</td></tr><tr><th>380</th><td></td><td>–20.9</td><td>–31.9</td><td>–32.7</td><td>–27.7</td><td>–30.4</td><td>–28.6</td><td>–25.1</td><td>–53.1</td></tr><tr><th>510</th><td></td><td>–21.5</td><td>–31.8</td><td>–34.0</td><td>–28.8</td><td>–28.9</td><td>–27.9</td><td>–24.4</td><td>–53.2</td></tr><tr><th>540</th><td></td><td>–21.2</td><td>–31.5</td><td>–32.8</td><td>–28.6</td><td>–29.4</td><td>–30.0</td><td>–23.4</td><td>–53.6</td></tr><tr><th>650</th><td></td><td>–22.0</td><td>–32.3</td><td>–32.5</td><td>–27.7</td><td>–29.8</td><td>–29.2</td><td>–23.6</td><td>–53.4</td></tr><tr><th>870</th><td></td><td>–20.8</td><td>–33.3</td><td>–32.6</td><td>–28.7</td><td>–29.8</td><td>–29.6</td><td>–22.2</td><td>–53.5</td></tr><tr><th>910</th><td></td><td>–21.2</td><td>–31.3</td><td>–31.7</td><td>–28.6</td><td>–30.1</td><td>–29.5</td><td>–25.5</td><td>–52.7</td></tr><tr><th>Mean</th><td></td><td><strong>–21.4</strong></td><td><strong>–32.1</strong></td><td><strong>–32.5</strong></td><td><strong>–28.3</strong></td><td><strong>–29.8</strong></td><td><strong>–29.4</strong></td><td><strong>–23.7</strong></td><td><strong>–53.2</strong></td></tr><tr><th>SD</th><td></td><td>0.5</td><td>0.6</td><td>0.8</td><td>0.7</td><td>0.6</td><td>0.7</td><td>1.0</td><td>0.3</td></tr></tbody></table>

The values in bold are mean values. Standard deviations are given one row below

<footnote><sup>a</sup>Values in the parenthesis are the standard deviations of N=20 replicate measurements</footnote>
(Table 2). In this work, we adopt the arithmetic average of these two models.

The site-specific carbon isotope delta values, δ<sub>t</sub>(<sup>13</sup>C), measured at the NRC by <sup>13</sup>C-qNMR for VANA-1 and VANB-1 samples, are presented in Tables 3–4. The δ<sub>t</sub>(<sup>13</sup>C) values are expressed in permilles relative to the VPDB, five measurements for each bottle. The consistency of δ<sub>t</sub>(<sup>13</sup>C) values from sample to sample range from 0.3 to 1.1‰ for both VANA-1 and VANB-1, as evaluated by the repeatability standard deviations from NRC measurements shown in Tables 3–4.

VANA-1 and VANB-1, have average δ(<sup>13</sup>C) values of –31‰ and –26‰ relative to the VPDB, respectively, consistent with synthetic origin. Vanillin from petrochemical origin typically has δ(<sup>13</sup>C) values from –36 to –25‰ and lignin-based biovanillin has δ(<sup>13</sup>C) values from –29 to –26‰ [2, 5]. Thus, these can be distinguished from natural vanillin from the vanilla pod having δ(<sup>13</sup>C) values from –22 to –14‰ [11]. Although the average isotopic composition of vanillin from VANA-1 and VANB-1 are somewhat similar, the difference in the isotopic composition

of carbon C8 (OCH<sub>3</sub>) between these two materials is extremely large, –54‰ and –31‰, which likely reflects the different synthetic routes for these vanillin samples. In addition, we observe substantial variations in the isotopic composition of carbon atoms in the benzene ring of vanillin for both VANA-1 and VANB-1. Both materials exhibit approximately 10‰ difference between the individual carbon atoms—well in excess of the repeatability standard deviation (<1‰). In both materials, C7 carbon atoms show the highest abundance of <sup>13</sup>C and the C2 atoms the lowest.

## Inter-laboratory comparison

To support the NRC certification work, an inter-laboratory comparison study was conducted with Nantes laboratory (France), which also used an adiabatic pulse sequence for <sup>1</sup>H decoupling on the 400 MHz spectrometer. To avoid variability associated with the NMR data processing, the free induction decays from vanillin samples

<table><caption>**Table 1** Longitudinal relaxation times of each carbon in vanillin in the presence of 0.1 M Cr(Acac)<sub>3</sub><sup>a</sup></caption><thead><tr><th></th><th>C1</th><th>C2</th><th>C3</th><th>C4</th><th>C5</th><th>C6</th><th>C7</th><th>C8</th></tr></thead><tbody><tr><th>T<sub>1</sub>/s</th><td>1.03 (0.07)</td><td>1.26 (0.09)</td><td>1.46 (0.08)</td><td>1.57 (0.05)</td><td>0.85 (0.04)</td><td>0.84 (0.05)</td><td>0.93 (0.03)</td><td>0.86<br>(0.06)</td></tr></tbody></table>

<table><caption>**Table 2** The effect of measurement model on the carbon isotope delta results for a single measurement of vanillin in VANA-1 sample from vial number 170 at NRC</caption><thead><tr><th>Model</th><th>C1</th><th>C2</th><th>C3</th><th>C4</th><th>C5</th><th>C6</th><th>C7</th><th>C8</th></tr></thead><tbody><tr><th>M0</th><td>–20.8</td><td>–30.7</td><td>–33.4</td><td>–28.7</td><td>–28.0</td><td>–30.0</td><td>–24.0</td><td>–54.8</td></tr><tr><th>M1</th><td>–20.4</td><td>–31.0</td><td>–33.6</td><td>–29.2</td><td>–28.2</td><td>–30.2</td><td>–24.3</td><td>–53.6</td></tr></tbody></table>

<table><caption>**Table 3** Site-specific average carbon isotope delta values of vanillin in VANA-1 samples measured at NRC</caption><thead><tr><th></th><th>VANA-1 vial number</th><th>C1</th><th>C2</th><th>C3</th><th>C4</th><th>C5</th><th>C6</th><th>C7</th><th>C8</th></tr></thead><tbody><tr><th>120</th><td></td><td>–22.1</td><td>–31.8</td><td>–32.6</td><td>–28.4</td><td>–29.6</td><td>–29.9</td><td>–23.2</td><td>–52.8</td></tr><tr><th>170</th><td></td><td>–21.1</td><td>–32.1</td><td>–33.1</td><td>–29.2</td><td>–29.4</td><td>–29.3</td><td>–22.8</td><td>–53.5</td></tr><tr><th>290</th><td></td><td>–21.1</td><td>–32.6</td><td>–31.8</td><td>–27.9</td><td>–30.6</td><td>–30.0</td><td>–23.0</td><td>–53.4</td></tr><tr><th>380</th><td></td><td>–21.8</td><td>–32.4</td><td>–31.4</td><td>–27.0</td><td>–30.4</td><td>–29.8</td><td>–24.2</td><td>–53.3</td></tr><tr><th>380</th><td></td><td>–20.9</td><td>–31.9</td><td>–32.7</td><td>–27.7</td><td>–30.4</td><td>–28.6</td><td>–25.1</td><td>–53.1</td></tr><tr><th>510</th><td></td><td>–21.5</td><td>–31.8</td><td>–34.0</td><td>–28.8</td><td>–28.9</td><td>–27.9</td><td>–24.4</td><td>–53.2</td></tr><tr><th>540</th><td></td><td>–21.2</td><td>–31.5</td><td>–32.8</td><td>–28.6</td><td>–29.4</td><td>–30.0</td><td>–23.4</td><td>–53.6</td></tr><tr><th>650</th><td></td><td>–22.0</td><td>–32.3</td><td>–32.5</td><td>–27.7</td><td>–29.8</td><td>–29.2</td><td>–23.6</td><td>–53.4</td></tr><tr><th>870</th><td></td><td>–20.8</td><td>–33.3</td><td>–32.6</td><td>–28.7</td><td>–29.8</td><td>–29.6</td><td>–22.2</td><td>–53.5</td></tr><tr><th>910</th><td></td><td>–21.2</td><td>–31.3</td><td>–31.7</td><td>–28.6</td><td>–30.1</td><td>–29.5</td><td>–25.5</td><td>–52.7</td></tr><tr><th>Mean</th><td></td><td><strong>–21.4</strong></td><td><strong>–32.1</strong></td><td><strong>–32.5</strong></td><td><strong>–28.3</strong></td><td><strong>–29.8</strong></td><td><strong>–29.4</strong></td><td><strong>–23.7</strong></td><td><strong>–53.2</strong></td></tr><tr><th>SD</th><td></td><td>0.5</td><td>0.6</td><td>0.8</td><td>0.7</td><td>0.6</td><td>0.7</td><td>1.0</td><td>0.3</td></tr></tbody></table>

The values in bold are mean values. Standard deviations are given one row below

<footnote><sup>a</sup>Values in the parenthesis are the standard deviations of N=20 replicate measurements</footnote>
were processed by the NRC using the same software and analysts. Given that NRC measurements were obtained in two separate measurement campaigns, we consider our results as if they are from three laboratories. The median standard deviation between the laboratory results is 1‰ (Tables 5–6).

## Isotope measurements on high-field NMR spectrometer

Tables 5–6 indicate that close agreement between the results of site-specific isotope values can be obtained with two qualified Bruker 400 MHz spectrometers in Nantes and Ottawa. In this work, we also compare the results obtained from a Bruker 400 and Bruker Avance HD 700 MHz spectrometers. The uniform decoupling over the entire range of ¹H chemical shifts of the doubly-labelled ethanol on the 400 and 700 MHz spectrometers is 2‰ and 5‰, respectively. This difference can be explained by differences in the coil architecture: a BBFO probe on the 400 MHz and a cryogenically cooled probe on the 700 MHz. Differences in the

shape and the length of each coil (for ¹H and ¹³C) may result in ¹H decoupling deficiency (¹H coil) or edge effects (¹³C coil). In the case of vanillin, we observed that the values of carbon isotopic composition of VANA-1 and VANB-1 on 700 MHz spectrometers are significantly different from those measured on the 400 MHz (Table 7): the value for C8 in VANA-1 observed from 700 MHz instrument is –78‰ whereas the average value from 400 MHz instrument is –53‰. However, the relative difference between the measured values of VANA-1 and VANB-1 remains nearly the same: 23‰ (700 MHz) and 24‰ (400 MHz). This finding further indicates that reliable isotope delta measurements require careful qualification and calibration of the instrument. Additionally, the measurements on 700 MHz instrument can be calibrated against the VANA-1 with the correction factors for each carbon atom (see Electronic Supplementary Material Table S2). The current single-pulse ¹³C-NMR sequence needs to be further investigated for incorporation with the high-field NMR such as 700 MHz. Currently, however, the high-field NMR is not particularly well-suited for these measurements.

<table><caption>**Table 4** Site-specific average carbon isotope delta values of vanillin in VANB-1 samples measured at NRC</caption><thead><tr><th>VANB-1 vial number</th><th>C1</th><th>C2</th><th>C3</th><th>C4</th><th>C5</th><th>C6</th><th>C7</th><th>C8</th></tr></thead><tbody><tr><td>5</td><td>-19.3</td><td>-31.2</td><td>-33.6</td><td>-24.6</td><td>-22.1</td><td>-26.9</td><td>-18.6</td><td>-30.7</td></tr><tr><td>115</td><td>-19.9</td><td>-29.7</td><td>-32.6</td><td>-24.8</td><td>-24.8</td><td>-26.5</td><td>-18.6</td><td>-29.9</td></tr><tr><td>465</td><td>-20.8</td><td>-29.9</td><td>-33.1</td><td>-24.1</td><td>-24.8</td><td>-26.7</td><td>-17.2</td><td>-30.2</td></tr><tr><td>585</td><td>-20.9</td><td>-29.6</td><td>-31.9</td><td>-23.9</td><td>-25.3</td><td>-26.5</td><td>-18.4</td><td>-30.4</td></tr><tr><td>715</td><td>-20.5</td><td>-29.9</td><td>-33.2</td><td>-24.5</td><td>-23.5</td><td>-26.9</td><td>-18.1</td><td>-30.2</td></tr><tr><td>725</td><td>-20.4</td><td>-30.5</td><td>-33.0</td><td>-23.5</td><td>-24.8</td><td>-26.3</td><td>-18.4</td><td>-30.0</td></tr><tr><td>785</td><td>-19.8</td><td>-29.5</td><td>-32.2</td><td>-25.6</td><td>-23.3</td><td>-24.9</td><td>-19.9</td><td>-31.5</td></tr><tr><td>855</td><td>-19.6</td><td>-29.6</td><td>-33.0</td><td>-24.6</td><td>-22.7</td><td>-26.0</td><td>-20.5</td><td>-30.9</td></tr><tr><td>905</td><td>-20.3</td><td>-29.6</td><td>-32.1</td><td>-24.3</td><td>-24.4</td><td>-27.2</td><td>-17.9</td><td>-31.0</td></tr><tr><td>995</td><td>-19.3</td><td>-29.4</td><td>-33.0</td><td>-24.9</td><td>-22.9</td><td>-25.7</td><td>-20.6</td><td>-31.1</td></tr><tr><td>Mean</td><td>-20.1</td><td>-29.9</td><td>-32.8</td><td>-24.5</td><td>-23.9</td><td>-26.4</td><td>-18.8</td><td>-30.6</td></tr><tr><td>SD</td><td>0.6</td><td>0.5</td><td>0.5</td><td>0.6</td><td>1.1</td><td>0.7</td><td>1.1</td><td>0.5</td></tr></tbody></table>

The values in bold are mean values. Standard deviations are given one row below

<table><caption>**Table 5** Inter-laboratory measurements of carbon isotope deltas in VANA-1</caption><thead><tr><th>Laboratory</th><th>C1</th><th>C2</th><th>C3</th><th>C4</th><th>C5</th><th>C6</th><th>C7</th><th>C8</th></tr></thead><tbody><tr><td>NRC (Apr 2021)</td><td>-21.3</td><td>-32.3</td><td>-32.6</td><td>-28.6</td><td>-29.8</td><td>-29.8</td><td>-22.9</td><td>-53.4</td></tr><tr><td>NRC (Jun-Oct 2021)</td><td>-21.5</td><td>-31.9</td><td>-32.5</td><td>-28.0</td><td>-29.9</td><td>-29.0</td><td>-24.5</td><td>-53.1</td></tr><tr><td>Nantes (Feb 2020)</td><td>-22.5</td><td>-31.7</td><td>-32.6</td><td>-31.0</td><td>-29.6</td><td>-26.7</td><td>-25.6</td><td>-50.7</td></tr><tr><td>SD</td><td>0.7</td><td>0.3</td><td>0.1</td><td>1.6</td><td>0.2</td><td>1.6</td><td>1.3</td><td>1.5</td></tr></tbody></table>

<table><caption>**Table 6** Inter-laboratory measurements of carbon isotope deltas in VANB-1</caption><thead><tr><th>Laboratory</th><th>C1</th><th>C2</th><th>C3</th><th>C4</th><th>C5</th><th>C6</th><th>C7</th><th>C8</th></tr></thead><tbody><tr><td>NRC (Mar-Apr 2021)</td><td>-20.0</td><td>-30.0</td><td>-32.7</td><td>-24.7</td><td>-23.9</td><td>-26.5</td><td>-18.4</td><td>-30.7</td></tr><tr><td>NRC (Jun-Oct 2021)</td><td>-20.1</td><td>-29.8</td><td>-32.8</td><td>-24.3</td><td>-23.8</td><td>-26.3</td><td>-19.2</td><td>-30.5</td></tr><tr><td>Nantes (Feb 2020)</td><td>-21.7</td><td>-31.4</td><td>-34.0</td><td>-25.9</td><td>-21.5</td><td>-24.4</td><td>-20.3</td><td>-27.7</td></tr><tr><td>SD</td><td>0.9</td><td>0.9</td><td>0.7</td><td>0.9</td><td>1.4</td><td>1.2</td><td>1.0</td><td>1.7</td></tr></tbody></table>
were processed by the NRC using the same software and analysts. Given that NRC measurements were obtained in two separate measurement campaigns, we consider our results as if they are from three laboratories. The median standard deviation between the laboratory results is 1‰ (Tables 5–6).

## Isotope measurements on high-field NMR spectrometer

Tables 5–6 indicate that close agreement between the results of site-specific isotope values can be obtained with two qualified Bruker 400 MHz spectrometers in Nantes and Ottawa. In this work, we also compare the results obtained from a Bruker 400 and Bruker Avance HD 700 MHz spectrometers. The uniform decoupling over the entire range of ¹H chemical shifts of the doubly-labelled ethanol on the 400 and 700 MHz spectrometers is 2‰ and 5‰, respectively. This difference can be explained by differences in the coil architecture: a BBFO probe on the 400 MHz and a cryogenically cooled probe on the 700 MHz. Differences in the

shape and the length of each coil (for ¹H and ¹³C) may result in ¹H decoupling deficiency (¹H coil) or edge effects (¹³C coil). In the case of vanillin, we observed that the values of carbon isotopic composition of VANA-1 and VANB-1 on 700 MHz spectrometers are significantly different from those measured on the 400 MHz (Table 7): the value for C8 in VANA-1 observed from 700 MHz instrument is –78‰ whereas the average value from 400 MHz instrument is –53‰. However, the relative difference between the measured values of VANA-1 and VANB-1 remains nearly the same: 23‰ (700 MHz) and 24‰ (400 MHz). This finding further indicates that reliable isotope delta measurements require careful qualification and calibration of the instrument. Additionally, the measurements on 700 MHz instrument can be calibrated against the VANA-1 with the correction factors for each carbon atom (see Electronic Supplementary Material Table S2). The current single-pulse ¹³C-NMR sequence needs to be further investigated for incorporation with the high-field NMR such as 700 MHz. Currently, however, the high-field NMR is not particularly well-suited for these measurements.

<table><caption>**Table 4** Site-specific average carbon isotope delta values of vanillin in VANB-1 samples measured at NRC</caption><thead><tr><th>VANB-1 vial number</th><th>C1</th><th>C2</th><th>C3</th><th>C4</th><th>C5</th><th>C6</th><th>C7</th><th>C8</th></tr></thead><tbody><tr><td>5</td><td>-19.3</td><td>-31.2</td><td>-33.6</td><td>-24.6</td><td>-22.1</td><td>-26.9</td><td>-18.6</td><td>-30.7</td></tr><tr><td>115</td><td>-19.9</td><td>-29.7</td><td>-32.6</td><td>-24.8</td><td>-24.8</td><td>-26.5</td><td>-18.6</td><td>-29.9</td></tr><tr><td>465</td><td>-20.8</td><td>-29.9</td><td>-33.1</td><td>-24.1</td><td>-24.8</td><td>-26.7</td><td>-17.2</td><td>-30.2</td></tr><tr><td>585</td><td>-20.9</td><td>-29.6</td><td>-31.9</td><td>-23.9</td><td>-25.3</td><td>-26.5</td><td>-18.4</td><td>-30.4</td></tr><tr><td>715</td><td>-20.5</td><td>-29.9</td><td>-33.2</td><td>-24.5</td><td>-23.5</td><td>-26.9</td><td>-18.1</td><td>-30.2</td></tr><tr><td>725</td><td>-20.4</td><td>-30.5</td><td>-33.0</td><td>-23.5</td><td>-24.8</td><td>-26.3</td><td>-18.4</td><td>-30.0</td></tr><tr><td>785</td><td>-19.8</td><td>-29.5</td><td>-32.2</td><td>-25.6</td><td>-23.3</td><td>-24.9</td><td>-19.9</td><td>-31.5</td></tr><tr><td>855</td><td>-19.6</td><td>-29.6</td><td>-33.0</td><td>-24.6</td><td>-22.7</td><td>-26.0</td><td>-20.5</td><td>-30.9</td></tr><tr><td>905</td><td>-20.3</td><td>-29.6</td><td>-32.1</td><td>-24.3</td><td>-24.4</td><td>-27.2</td><td>-17.9</td><td>-31.0</td></tr><tr><td>995</td><td>-19.3</td><td>-29.4</td><td>-33.0</td><td>-24.9</td><td>-22.9</td><td>-25.7</td><td>-20.6</td><td>-31.1</td></tr><tr><td>Mean</td><td>-20.1</td><td>-29.9</td><td>-32.8</td><td>-24.5</td><td>-23.9</td><td>-26.4</td><td>-18.8</td><td>-30.6</td></tr><tr><td>SD</td><td>0.6</td><td>0.5</td><td>0.5</td><td>0.6</td><td>1.1</td><td>0.7</td><td>1.1</td><td>0.5</td></tr></tbody></table>

The values in bold are mean values. Standard deviations are given one row below

<table><caption>**Table 5** Inter-laboratory measurements of carbon isotope deltas in VANA-1</caption><thead><tr><th>Laboratory</th><th>C1</th><th>C2</th><th>C3</th><th>C4</th><th>C5</th><th>C6</th><th>C7</th><th>C8</th></tr></thead><tbody><tr><td>NRC (Apr 2021)</td><td>-21.3</td><td>-32.3</td><td>-32.6</td><td>-28.6</td><td>-29.8</td><td>-29.8</td><td>-22.9</td><td>-53.4</td></tr><tr><td>NRC (Jun-Oct 2021)</td><td>-21.5</td><td>-31.9</td><td>-32.5</td><td>-28.0</td><td>-29.9</td><td>-29.0</td><td>-24.5</td><td>-53.1</td></tr><tr><td>Nantes (Feb 2020)</td><td>-22.5</td><td>-31.7</td><td>-32.6</td><td>-31.0</td><td>-29.6</td><td>-26.7</td><td>-25.6</td><td>-50.7</td></tr><tr><td>SD</td><td>0.7</td><td>0.3</td><td>0.1</td><td>1.6</td><td>0.2</td><td>1.6</td><td>1.3</td><td>1.5</td></tr></tbody></table>

<table><caption>**Table 6** Inter-laboratory measurements of carbon isotope deltas in VANB-1</caption><thead><tr><th>Laboratory</th><th>C1</th><th>C2</th><th>C3</th><th>C4</th><th>C5</th><th>C6</th><th>C7</th><th>C8</th></tr></thead><tbody><tr><td>NRC (Mar-Apr 2021)</td><td>-20.0</td><td>-30.0</td><td>-32.7</td><td>-24.7</td><td>-23.9</td><td>-26.5</td><td>-18.4</td><td>-30.7</td></tr><tr><td>NRC (Jun-Oct 2021)</td><td>-20.1</td><td>-29.8</td><td>-32.8</td><td>-24.3</td><td>-23.8</td><td>-26.3</td><td>-19.2</td><td>-30.5</td></tr><tr><td>Nantes (Feb 2020)</td><td>-21.7</td><td>-31.4</td><td>-34.0</td><td>-25.9</td><td>-21.5</td><td>-24.4</td><td>-20.3</td><td>-27.7</td></tr><tr><td>SD</td><td>0.9</td><td>0.9</td><td>0.7</td><td>0.9</td><td>1.4</td><td>1.2</td><td>1.0</td><td>1.7</td></tr></tbody></table>
<table><caption>**Table 7** Measurements of VANA-1 and VANB-1 on 700 MHz NMR instrument</caption><thead><tr><th rowspan="2">Atom</th><th>δ<sub>t</sub>(<sup>13</sup>C,<br>VANA-1) vs<br>VPDB</th><th>δ<sub>t</sub>(<sup>13</sup>C,<br>VANB-1) vs<br>VPDB</th><th>δ<sub>t</sub>(<sup>13</sup>C,<br>VANB-1) vs<br>VANA-1 *</th><th>δ<sub>t</sub>(<sup>13</sup>C, VANB-<br>1) vs VANA-1<br>*</th></tr><tr><th>(700 MHz)</th><th>(700 MHz)</th><th>(700 MHz)</th><th>(400 MHz)</th></tr></thead><tbody><tr><td>C1</td><td>-10.5</td><td>-9.5</td><td>0.9</td><td>1.3</td></tr><tr><td>C2</td><td>-17.6</td><td>-12.7</td><td>4.9</td><td>2.3</td></tr><tr><td>C3</td><td>-19.3</td><td>-17.1</td><td>2.2</td><td>-0.3</td></tr><tr><td>C4</td><td>-20.1</td><td>-14.2</td><td>6.1</td><td>3.9</td></tr><tr><td>C5</td><td>-33.9</td><td>-30.0</td><td>4.1</td><td>6.1</td></tr><tr><td>C6</td><td>-30.1</td><td>-28.3</td><td>1.9</td><td>3.1</td></tr><tr><td>C7</td><td>-40.8</td><td>-37.9</td><td>3.0</td><td>5.0</td></tr><tr><td>C8</td><td>-78.1</td><td>-57.1</td><td>22.7</td><td>23.9</td></tr></tbody></table>

## Further validation of the <sup>13</sup>C results: gravimetric sample mixing

At the NRC, we have performed additional analyses of VANA-1 and VANB-1 mixtures. In essence, the gravimetric mixing of the two materials should produce isotopic compositions that can be predicted beforehand from the mixing ratios and the isotopic composition of the two materials. Figure 5 shows the results for two carbon atoms (C5 and C8) whose results show a median residual of 0.85‰ from the linear model fits.

## Measurement uncertainty

We have evaluated the uncertainty associated with the site-specific carbon isotope values using several approaches. First, we have assessed the effect of the measurement model (Table 2), and the laboratory effect (Tables 5–6). Both show uncertainties below 1‰. Thus, we performed a Bayesian evaluation of the overall measurement uncertainty which incorporated uncertainty contributions from the isotopic composition of the

VPDB, the reliability of peak area measurements (set to $u=2‰$ as assessed from the experiments using doubly-labelled acetic acid), and the uncertainty associated with the bulk isotope value from the IRMS measurements (see Electronic Supplementary Material section S3). The analysis is shown in the Electronic Supplementary Material and leads to the measurement combined standard uncertainty of 0.7‰ (or expanded uncertainty of 1.4‰ at a 95% confidence level) for the individual carbon isotope values in either of the two vanillin materials.

## Homogeneity and stability of the materials

The batch homogeneity of carbon isotope values in VANA-1 and VANB-1 was assessed using the NRC measurements (Tables 3–4) and the aforementioned estimate of the measurement uncertainty (0.7‰) by fitting the random effects model to these data using *NIST Consensus Builder* [33]. The median uncertainty due to homogeneity, averaged over all carbon atoms and both materials, is $u_{hom}=0.3‰$.

An isochronous short-term stability of VANA-1 and VANB-1 was conducted to examine the effect of sample storage. For each material, two vials were stored, one in the oven (+40 °C) and one at room temperature. After two and four weeks, the bulk $\delta(^{13}\text{C})$ and carbon site-specific isotope delta values of each vial were measured at the NRC by IRMS and <sup>13</sup>C-qNMR. No significant effect in the $\delta(^{13}\text{C})$ values was observed for the samples stored at +40 °C.

## Summary results

Our evaluation of the uncertainties associated with the site-specific carbon isotope values for all carbon atoms in the two materials can be summarized as follows. The uncertainty due to characterization, $u_{char}=0.3‰$, the uncertainty due to homogeneity, $u_{hom}=0.3‰$, the uncertainty due to short-term stability, $u_{stab}=0.0‰$, and the uncertainty due to laboratory effect, $u_{lab}=0.6‰$ ($1.0‰/\sqrt{3}$). Combining these sources of

<figure><img src="image_8.png" /><figcaption>**Fig. 5** Isotope delta values for vanillin carbon atoms C5 and C8 measured from gravimetric mixtures of VANA-1 and VANB-1 at NRC. Horizontal axis shows the mass fraction of VANA-1 in such mixtures, $w_{VANA-1} = m_{VANA-1}/(m_{VANA-1} + m_{VANB-1})$</figcaption></figure>

<footnote>*Difference between the certified values. Note that $\delta_{VANA}(VANB) = [1 + \delta_{VPDB}(VANB)]/[1 + \delta_{VPDB}(VANA)] - 1$</footnote>
<table><caption>**Table 7** Measurements of VANA-1 and VANB-1 on 700 MHz NMR instrument</caption><thead><tr><th rowspan="2">Atom</th><th>δ<sub>t</sub>(<sup>13</sup>C,<br>VANA-1) vs<br>VPDB</th><th>δ<sub>t</sub>(<sup>13</sup>C,<br>VANB-1) vs<br>VPDB</th><th>δ<sub>t</sub>(<sup>13</sup>C,<br>VANB-1) vs<br>VANA-1 *</th><th>δ<sub>t</sub>(<sup>13</sup>C, VANB-<br>1) vs VANA-1<br>*</th></tr><tr><th>(700 MHz)</th><th>(700 MHz)</th><th>(700 MHz)</th><th>(400 MHz)</th></tr></thead><tbody><tr><td>C1</td><td>-10.5</td><td>-9.5</td><td>0.9</td><td>1.3</td></tr><tr><td>C2</td><td>-17.6</td><td>-12.7</td><td>4.9</td><td>2.3</td></tr><tr><td>C3</td><td>-19.3</td><td>-17.1</td><td>2.2</td><td>-0.3</td></tr><tr><td>C4</td><td>-20.1</td><td>-14.2</td><td>6.1</td><td>3.9</td></tr><tr><td>C5</td><td>-33.9</td><td>-30.0</td><td>4.1</td><td>6.1</td></tr><tr><td>C6</td><td>-30.1</td><td>-28.3</td><td>1.9</td><td>3.1</td></tr><tr><td>C7</td><td>-40.8</td><td>-37.9</td><td>3.0</td><td>5.0</td></tr><tr><td>C8</td><td>-78.1</td><td>-57.1</td><td>22.7</td><td>23.9</td></tr></tbody></table>

## Further validation of the <sup>13</sup>C results: gravimetric sample mixing

At the NRC, we have performed additional analyses of VANA-1 and VANB-1 mixtures. In essence, the gravimetric mixing of the two materials should produce isotopic compositions that can be predicted beforehand from the mixing ratios and the isotopic composition of the two materials. Figure 5 shows the results for two carbon atoms (C5 and C8) whose results show a median residual of 0.85‰ from the linear model fits.

## Measurement uncertainty

We have evaluated the uncertainty associated with the site-specific carbon isotope values using several approaches. First, we have assessed the effect of the measurement model (Table 2), and the laboratory effect (Tables 5–6). Both show uncertainties below 1‰. Thus, we performed a Bayesian evaluation of the overall measurement uncertainty which incorporated uncertainty contributions from the isotopic composition of the

VPDB, the reliability of peak area measurements (set to $u=2‰$ as assessed from the experiments using doubly-labelled acetic acid), and the uncertainty associated with the bulk isotope value from the IRMS measurements (see Electronic Supplementary Material section S3). The analysis is shown in the Electronic Supplementary Material and leads to the measurement combined standard uncertainty of 0.7‰ (or expanded uncertainty of 1.4‰ at a 95% confidence level) for the individual carbon isotope values in either of the two vanillin materials.

## Homogeneity and stability of the materials

The batch homogeneity of carbon isotope values in VANA-1 and VANB-1 was assessed using the NRC measurements (Tables 3–4) and the aforementioned estimate of the measurement uncertainty (0.7‰) by fitting the random effects model to these data using *NIST Consensus Builder* [33]. The median uncertainty due to homogeneity, averaged over all carbon atoms and both materials, is $u_{hom}=0.3‰$.

An isochronous short-term stability of VANA-1 and VANB-1 was conducted to examine the effect of sample storage. For each material, two vials were stored, one in the oven (+40 °C) and one at room temperature. After two and four weeks, the bulk $\delta(^{13}\text{C})$ and carbon site-specific isotope delta values of each vial were measured at the NRC by IRMS and <sup>13</sup>C-qNMR. No significant effect in the $\delta(^{13}\text{C})$ values was observed for the samples stored at +40 °C.

## Summary results

Our evaluation of the uncertainties associated with the site-specific carbon isotope values for all carbon atoms in the two materials can be summarized as follows. The uncertainty due to characterization, $u_{char}=0.3‰$, the uncertainty due to homogeneity, $u_{hom}=0.3‰$, the uncertainty due to short-term stability, $u_{stab}=0.0‰$, and the uncertainty due to laboratory effect, $u_{lab}=0.6‰$ ($1.0‰/\sqrt{3}$). Combining these sources of

<figure><img src="image_8.png" /><figcaption>**Fig. 5** Isotope delta values for vanillin carbon atoms C5 and C8 measured from gravimetric mixtures of VANA-1 and VANB-1 at NRC. Horizontal axis shows the mass fraction of VANA-1 in such mixtures, $w_{VANA-1} = m_{VANA-1}/(m_{VANA-1} + m_{VANB-1})$</figcaption></figure>

<footnote>*Difference between the certified values. Note that $\delta_{VANA}(VANB) = [1 + \delta_{VPDB}(VANB)]/[1 + \delta_{VPDB}(VANA)] - 1$</footnote>
<table><caption>**Table 8** Summary of site-specific carbon isotope delta values of vanillin in VANA-1 and VANB-1 (relative to the VPDB) and their expanded uncertainties at 95% confidence</caption><thead><tr><th>Carbon atom</th><th>VANA-1</th><th>VANB-1</th></tr></thead><tbody><tr><td>C1</td><td>-21.4 ± 1.4‰</td><td>-20.1 ± 1.4‰</td></tr><tr><td>C2</td><td>-32.1 ± 1.4‰</td><td>-29.9 ± 1.4‰</td></tr><tr><td>C3</td><td>-32.5 ± 1.4‰</td><td>-32.8 ± 1.4‰</td></tr><tr><td>C4</td><td>-28.3 ± 1.4‰</td><td>-24.5 ± 1.4‰</td></tr><tr><td>C5</td><td>-29.8 ± 1.4‰</td><td>-23.9 ± 1.4‰</td></tr><tr><td>C6</td><td>-29.4 ± 1.4‰</td><td>-26.4 ± 1.4‰</td></tr><tr><td>C7</td><td>-23.7 ± 1.4‰</td><td>-18.8 ± 1.4‰</td></tr><tr><td>C8</td><td>-53.2 ± 1.4‰</td><td>-30.6 ± 1.4‰</td></tr><tr><td>IRMS average</td><td>-31.30 ± 0.06‰</td><td>-25.85 ± 0.05‰</td></tr></tbody></table>

uncertainty yields combined standard uncertainty $u=0.7‰$ ($k=1$) or expanded uncertainty $U=1.4‰$ ($k=2$) at the 95% confidence level, applicable to all carbon isotope delta values for both materials. Table 8 summarizes the results for both materials with the average NRC results from Tables 3–4.

## Conclusion

The combination of ¹³C-qNMR and IRMS can provide site-specific isotope values for each carbon atom in vanillin. Our measurements for VANA-1 and VANB-1 reference materials provide results that are reliable to within approx. 1‰. While such a precision is an order of magnitude larger than typical IRMS measurements can provide, it is sufficient to show significant differences between the isotopic compositions of the various vanillin carbon atoms in these materials, even showing significant differences within the benzene ring. Our certified carbon isotope values of VANA-1 and VANB-1 were obtained with a qualified Bruker 400 MHz NMR instrument. However, since these values may not be recovered when samples are measured on other NMR spectrometers, one can use VANA-1 as the delta-zero material for intramolecular carbon isotope measurements. To our knowledge, these measurements and the associated reference materials represent the current state of the art in the field and we believe that these two vanillin CRMs, both available from the NRC, will prove useful to the stable isotope community and contribute to detecting counterfeiting of vanillin and other materials of significant economic value.

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1007/s00216-022-04292-0.

**Acknowledgements** The authors from CEISAM acknowledge the French National Infrastructure for Metabolomics and Fluxomics MetaboHUB-ANR-11-INBS-0010 (www.metabohub.fr) and the Corsaire metabolomics core facility (Biogenouest). This work includes NMR experiments carried out on the CEISAM NMR platform.

**Funding** Open Access provided by National Research Council Canada.

## Declarations

**Competing interests** The authors declare no competing interests.

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

## References

1. Bomgarder MM. The problem with vanilla. Chem Eng News. 2016;94:38–42.

2. Gallage NJ, Møller BL. Vanillin-bioconversion and bioengineering of the most popular plant flavor and its de novo biosynthesis in the vanilla orchid. Mol Plant. 2015;8:40–57.

3. Kaunzinger A, Juchelka D, Mosandl A. Progress in the authenticity assessment of vanilla 1 Initiation of authenticity profiles. J Agric Food Chem. 1997;45:1752–7.

4. Lamprecht G, Pichlmayer F, Schmid ER. Determination of the authenticity of vanilla extracts by stable isotope ratio analysis and component analysis by HPLC. J Agric Food Chem. 1994;42:1722–7.

5. Priefert H, Rabenhorst J, Steinbüchel A. Biotechnological production of vanillin. Appl Microbiol Biotechnol. 2001;56:296–314.

6. Muccioza Z, Jackson GP. Isotope ratio mass spectrometry. Analyst. 2009;134:213–22.

7. Jamin E, Thomas F. SNIF-NMR applications in an economic context: fraud detection in food products. In Mod Magn Reson. Edited by Springer C 2017: 1–12

8. Fronza G, Fuganti C, Serra S, Burke A, Guillou C, Reniero F. The positional δ(18O) values of extracted and synthetic vanillin. Helv Chim Acta. 2001;84:351–9.

9. Krueger DA, Krueger HW. Carbon isotopes in vanillin and the detection of falsified “natural” vanillin. J Agric Food Chem. 1983;31:1265–8.

10. Hansen A-M, Fromberg A, Frandsen H. Authenticity and traceability of vanilla flavors by analysis of stable isotopes of carbon and hydrogen. J Agric Food Chem. 2014;62:10326–31.

11. Wilde A, Frandsen H, Fromberg A, Smedsgaard J, Greule M. Isotopic characterization of vanillin ex glucose by GC-IRMS - new challenge for natural vanilla flavour authentication? Food Control. 2019;106:106735.

12. Greule M, Tumino LD, Kronewald T, Hener U, Schleucher Jr, Mosandl A, Keppler F. Improved rapid authentication of vanillin using d13C and d2 H values. Eur Food Res Technol. 2010;231:933–41.

13. Remaud GS, Martin YL, Martin GG, Martin GJ. Detection of sophisticated adulterations of natural vanilla flavors and extracts: application of the SNIF-NMR method to vanillin and p-hydroxybenzaldehyde. J Agric Food Chem. 1997;45:859–66.

14. Akoka S, Remaud GS. NMR-based isotopic and isotopomic analysis. Prog Nucl Magn Reson Spectrosc. 2020;120–121:1–24.

15. Guyader S, Thomas F, Jamin E, Grand M, Akoka S, Silvestre V, Ali E. Combination of 13C and 2H SNIF-NMR isotopic fingerprints
<table><caption>**Table 8** Summary of site-specific carbon isotope delta values of vanillin in VANA-1 and VANB-1 (relative to the VPDB) and their expanded uncertainties at 95% confidence</caption><thead><tr><th>Carbon atom</th><th>VANA-1</th><th>VANB-1</th></tr></thead><tbody><tr><td>C1</td><td>-21.4 ± 1.4‰</td><td>-20.1 ± 1.4‰</td></tr><tr><td>C2</td><td>-32.1 ± 1.4‰</td><td>-29.9 ± 1.4‰</td></tr><tr><td>C3</td><td>-32.5 ± 1.4‰</td><td>-32.8 ± 1.4‰</td></tr><tr><td>C4</td><td>-28.3 ± 1.4‰</td><td>-24.5 ± 1.4‰</td></tr><tr><td>C5</td><td>-29.8 ± 1.4‰</td><td>-23.9 ± 1.4‰</td></tr><tr><td>C6</td><td>-29.4 ± 1.4‰</td><td>-26.4 ± 1.4‰</td></tr><tr><td>C7</td><td>-23.7 ± 1.4‰</td><td>-18.8 ± 1.4‰</td></tr><tr><td>C8</td><td>-53.2 ± 1.4‰</td><td>-30.6 ± 1.4‰</td></tr><tr><td>IRMS average</td><td>-31.30 ± 0.06‰</td><td>-25.85 ± 0.05‰</td></tr></tbody></table>

uncertainty yields combined standard uncertainty $u=0.7‰$ ($k=1$) or expanded uncertainty $U=1.4‰$ ($k=2$) at the 95% confidence level, applicable to all carbon isotope delta values for both materials. Table 8 summarizes the results for both materials with the average NRC results from Tables 3–4.

## Conclusion

The combination of ¹³C-qNMR and IRMS can provide site-specific isotope values for each carbon atom in vanillin. Our measurements for VANA-1 and VANB-1 reference materials provide results that are reliable to within approx. 1‰. While such a precision is an order of magnitude larger than typical IRMS measurements can provide, it is sufficient to show significant differences between the isotopic compositions of the various vanillin carbon atoms in these materials, even showing significant differences within the benzene ring. Our certified carbon isotope values of VANA-1 and VANB-1 were obtained with a qualified Bruker 400 MHz NMR instrument. However, since these values may not be recovered when samples are measured on other NMR spectrometers, one can use VANA-1 as the delta-zero material for intramolecular carbon isotope measurements. To our knowledge, these measurements and the associated reference materials represent the current state of the art in the field and we believe that these two vanillin CRMs, both available from the NRC, will prove useful to the stable isotope community and contribute to detecting counterfeiting of vanillin and other materials of significant economic value.

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1007/s00216-022-04292-0.

**Acknowledgements** The authors from CEISAM acknowledge the French National Infrastructure for Metabolomics and Fluxomics MetaboHUB-ANR-11-INBS-0010 (www.metabohub.fr) and the Corsaire metabolomics core facility (Biogenouest). This work includes NMR experiments carried out on the CEISAM NMR platform.

**Funding** Open Access provided by National Research Council Canada.

## Declarations

**Competing interests** The authors declare no competing interests.

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

## References

1. Bomgarder MM. The problem with vanilla. Chem Eng News. 2016;94:38–42.

2. Gallage NJ, Møller BL. Vanillin-bioconversion and bioengineering of the most popular plant flavor and its de novo biosynthesis in the vanilla orchid. Mol Plant. 2015;8:40–57.

3. Kaunzinger A, Juchelka D, Mosandl A. Progress in the authenticity assessment of vanilla 1 Initiation of authenticity profiles. J Agric Food Chem. 1997;45:1752–7.

4. Lamprecht G, Pichlmayer F, Schmid ER. Determination of the authenticity of vanilla extracts by stable isotope ratio analysis and component analysis by HPLC. J Agric Food Chem. 1994;42:1722–7.

5. Priefert H, Rabenhorst J, Steinbüchel A. Biotechnological production of vanillin. Appl Microbiol Biotechnol. 2001;56:296–314.

6. Muccio Z, Jackson GP. Isotope ratio mass spectrometry. Analyst. 2009;134:213–22.

7. Jamin E, Thomas F. SNIF-NMR applications in an economic context: fraud detection in food products. In *Mod Magn Reson*. Edited by Springer C 2017: 1–12

8. Fronza G, Fuganti C, Serra S, Burke A, Guillou C, Reniero F. The positional δ(18O) values of extracted and synthetic vanillin. Helv Chim Acta. 2001;84:351–9.

9. Krueger DA, Krueger HW. Carbon isotopes in vanillin and the detection of falsified “natural” vanillin. J Agric Food Chem. 1983;31:1265–8.

10. Hansen A-M, Fromberg A, Frandsen H. Authenticity and traceability of vanilla flavors by analysis of stable isotopes of carbon and hydrogen. J Agric Food Chem. 2014;62:10326–31.

11. Wilde A, Frandsen H, Fromberg A, Smedsgaard J, Greule M. Isotopic characterization of vanillin ex glucose by GC-IRMS - new challenge for natural vanilla flavour authentication? Food Control. 2019;106:106735.

12. Greule M, Tumino LD, Kronewald T, Hener U, Schleucher Jr, Mosandl A, Keppler F. Improved rapid authentication of vanillin using d13C and d2 H values. Eur Food Res Technol. 2010;231:933–41.

13. Remaud GS, Martin YL, Martin GG, Martin GJ. Detection of sophisticated adulterations of natural vanilla flavors and extracts: application of the SNIF-NMR method to vanillin and p-hydroxybenzaldehyde. J Agric Food Chem. 1997;45:859–66.

14. Akoka S, Remaud GS. NMR-based isotopic and isotopomic analysis. Prog Nucl Magn Reson Spectrosc. 2020;120–121:1–24.

15. Guyader S, Thomas F, Jamin E, Grand M, Akoka S, Silvestre V, Ali E. Combination of 13C and 2H SNIF-NMR isotopic fingerprints
of vanillin to control its precursors. Flavour and Fragrance Journal.
2019;34:133–44.

16. Hoffman DW, Rasmussen C. Absolute carbon stable isotope ratio in the Vienna PeeDee Belemnite isotope reference determined by 1H NMR spectroscopy. Anal Chem. 2022;94:5240–7.

17. Toulemonde B, Horman I, Egli HMD. Food-related applications of high-resolution NMR. Part II. Differentiation between natural and synthetic vanillin samples using 2H-NMR. Helvetica Chimica Acta. 1983;66:2342–5.

18. Remaud GS, Akoka S. A review of flavors authentication by position-specific isotope analysis by nuclear magnetic resonance spectrometry: the example of vanillin. Flavour Fragr J. 2017;32:77–84.

19. Caytan E, Botosoa EP, Silvestre V, Robins RJ, Akoka S, Remaud GS. Accurate quantitative 13C NMR spectroscopy: repeatability over time of site-specific 13C isotope ratio determination. Anal Chem. 2007;79:8266–9.

20. Caytan E, Remaud GS, Tenailleau E, Akoka S. Precise and accurate quantitative 13C-NMR with reduced experimental time. Talanta. 2007;71:1016–21.

21. Tenailleau E, Akoka S. Adiabatic 1H decoupling scheme for very accurate intensity measurements in 13C NMR. J Magn Reson. 2007;185:50–8.

22. Tenailleau E, Remaud G, Akoka S. Quantification of the 1H-decoupling effects on the accuracy of 13C-NMR measurements. Instrum Sci Technol. 2005;33:391–9.

23. Chaintreau A, Fieber W, Sommer H, Gilbert A, Yamada K, Yoshida N, Pagelot A, Moskau D, Moreno A, Schleucher J, et al. Site-specific 13C content by quantitative isotopic 13C nuclear magnetic resonance spectrometry: a pilot inter-laboratory study. Anal Chim Acta. 2013;788:108–13.

24. Chartrand MMG, Meija J, Kumkrong P, Mester Z. Three certified sugar reference materials for carbon isotope delta measurements. Rapid Commun Mass Spectrom. 2019;33:272–80.

25. Chartrand MMG; Meija J; Hélie J-F; Middlestead P; Ramalingam M; Aziz AA; Mester Z. Development of vanillin carbon isotope delta reference materials Anal Bioanal Chem. 2022, Submitted.

26. Le PM, Ding J, Leek D, Mester Z, Robertson G, Windust A, Meija J. Determination of chemical purity and isotopic composition of natural and carbon-13-labeled arsenobetaine bromide standards by quantitative 1H-NMR. Anal Bional Chem. 2016;408:7413–21.

27. Sokolenko S, Jézéquel T, Hajjar G, Farjon J, Akoka S, Giraudeau P. Robust 1D NMR lineshape fitting using real and imaginary data in the frequency domain. J Meg Res. 2019;298:91–100.

28. Malinovsky D, Dunn PJH, Holcombe G, Cowen S, Goenaga-Infante H. Development and characterisation of new glycine certified reference materials for SI-traceable 13C/12C isotope amount ratio measurements. J Anal At Spectrom. 2019;34:147–59.

29. Fleisher AJ, Yi H, Srivastava A, Polyansky OL, Zobov NF, Hodges JT. Absolute 13C/12C isotope amount ratio for Vienna PeeDee Belemnite from infrared absorption spectroscopy. Nat Phys. 2021;17:889–93.

30. Silvestre V, Mboula VM, Jouitteau C, Akoka S, Robins RJ, Remaud GS. Isotopic 13C NMR spectrometry to assess counterfeiting of active pharmaceutical ingredients: site-specific 13C content of aspirin and paracetamol. J Pharm Biomed Anal. 2009;50:336–41.

31. Robins RJ; Remaud GS; Akoka S. Isotope ratio monitoring 13C nuclear magnetic resonance spectrometry for the analysis of position-specific isotope ratios. Elsevier Inc edn2017.

32. Bharti S, Roy R. The 1H NMR technique is currently used in purity determinations of organic substances. Trends Anal Chem. 2012;35:5–26.

33. Koepke A, Lafarge T, Possolo A, Toman B. Consensus building for interlaboratory studies, key comparisons, and meta-analysis. Metrologia. 2017;54:S34–62.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_9.png" /></figure>

**Phuong Mai Le** is a Research Council Officer at the National Research Council Canada. Her research focuses on the investigation and application of the advanced NMR and MS approaches for (i) analysis of the intramolecular isotope profile of molecules at natural abundance and the development of certified reference materials; (ii) identification and quantification of compounds in complex mixtures, especially in food and medicinal plants.

<figure><img src="image_10.png" /></figure>

**Estelle Martineau** is a research engineer at CAPACITÉS SAS, Nantes (France). She is a specialist in highly precise quantitative analysis and in analyzing complex mixtures by 1D and 2D NMR.

<figure><img src="image_11.png" /></figure>

**Serge Akoka** is a Full Professor of analytical chemistry at Nantes Université (France). His current research interests are (i) isotopic analysis by NMR at natural abundance and (ii) methodological developments in quantitative 1D and 2D NMR.

<figure><img src="image_12.png" /></figure>

**Gerald Remaud** tenured Full Professor at Nantes Université (France), is a specialist in the development of methods for intramolecular isotope measurements and their applications on authenticity, mainly in food. Previously, he held a number of positions in the industry, following research in Sweden into NMR methods for studying nucleotide chemistry.

<figure><img src="image_13.png" /></figure>
of vanillin to control its precursors. Flavour and Fragrance Journal.
2019;34:133–44.

16. Hoffman DW, Rasmussen C. Absolute carbon stable isotope ratio in the Vienna PeeDee Belemnite isotope reference determined by 1H NMR spectroscopy. Anal Chem. 2022;94:5240–7.

17. Toulemonde B, Horman I, Egli HMD. Food-related applications of high-resolution NMR. Part II. Differentiation between natural and synthetic vanillin samples using 2H-NMR. Helvetica Chimica Acta. 1983;66:2342–5.

18. Remaud GS, Akoka S. A review of flavors authentication by position-specific isotope analysis by nuclear magnetic resonance spectrometry: the example of vanillin. Flavour Fragr J. 2017;32:77–84.

19. Caytan E, Botosoa EP, Silvestre V, Robins RJ, Akoka S, Remaud GS. Accurate quantitative 13C NMR spectroscopy: repeatability over time of site-specific 13C isotope ratio determination. Anal Chem. 2007;79:8266–9.

20. Caytan E, Remaud GS, Tenailleau E, Akoka S. Precise and accurate quantitative 13C-NMR with reduced experimental time. Talanta. 2007;71:1016–21.

21. Tenailleau E, Akoka S. Adiabatic 1H decoupling scheme for very accurate intensity measurements in 13C NMR. J Magn Reson. 2007;185:50–8.

22. Tenailleau E, Remaud G, Akoka S. Quantification of the 1H-decoupling effects on the accuracy of 13C-NMR measurements. Instrum Sci Technol. 2005;33:391–9.

23. Chaintreau A, Fieber W, Sommer H, Gilbert A, Yamada K, Yoshida N, Pagelot A, Moskau D, Moreno A, Schleucher J, et al. Site-specific 13C content by quantitative isotopic 13C nuclear magnetic resonance spectrometry: a pilot inter-laboratory study. Anal Chim Acta. 2013;788:108–13.

24. Chartrand MMG, Meija J, Kumkrong P, Mester Z. Three certified sugar reference materials for carbon isotope delta measurements. Rapid Commun Mass Spectrom. 2019;33:272–80.

25. Chartrand MMG; Meija J; Hélie J-F; Middlestead P; Ramalingam M; Aziz AA; Mester Z. Development of vanillin carbon isotope delta reference materials Anal Bioanal Chem. 2022, Submitted.

26. Le PM, Ding J, Leek D, Mester Z, Robertson G, Windust A, Meija J. Determination of chemical purity and isotopic composition of natural and carbon-13-labeled arsenobetaine bromide standards by quantitative 1H-NMR. Anal Bional Chem. 2016;408:7413–21.

27. Sokolenko S, Jézéquel T, Hajjar G, Farjon J, Akoka S, Giraudeau P. Robust 1D NMR lineshape fitting using real and imaginary data in the frequency domain. J Meg Res. 2019;298:91–100.

28. Malinovsky D, Dunn PJH, Holcombe G, Cowen S, Goenaga-Infante H. Development and characterisation of new glycine certified reference materials for SI-traceable 13C/12C isotope amount ratio measurements. J Anal At Spectrom. 2019;34:147–59.

29. Fleisher AJ, Yi H, Srivastava A, Polyansky OL, Zobov NF, Hodges JT. Absolute 13C/12C isotope amount ratio for Vienna PeeDee Belemnite from infrared absorption spectroscopy. Nat Phys. 2021;17:889–93.

30. Silvestre V, Mboula VM, Jouitteau C, Akoka S, Robins RJ, Remaud GS. Isotopic 13C NMR spectrometry to assess counterfeiting of active pharmaceutical ingredients: site-specific 13C content of aspirin and paracetamol. J Pharm Biomed Anal. 2009;50:336–41.

31. Robins RJ; Remaud GS; Akoka S. Isotope ratio monitoring 13C nuclear magnetic resonance spectrometry for the analysis of position-specific isotope ratios. Elsevier Inc edn2017.

32. Bharti S, Roy R. The 1H NMR technique is currently used in purity determinations of organic substances. Trends Anal Chem. 2012;35:5–26.

33. Koepke A, Lafarge T, Possolo A, Toman B. Consensus building for interlaboratory studies, key comparisons, and meta-analysis. Metrologia. 2017;54:S34–62.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_9.png" /></figure>

**Phuong Mai Le** is a Research Council Officer at the National Research Council Canada. Her research focuses on the investigation and application of the advanced NMR and MS approaches for (i) analysis of the intramolecular isotope profile of molecules at natural abundance and the development of certified reference materials; (ii) identification and quantification of compounds in complex mixtures, especially in food and medicinal plants.

<figure><img src="image_10.png" /></figure>

**Estelle Martineau** is a research engineer at CAPACITÉS SAS, Nantes (France). She is a specialist in highly precise quantitative analysis and in analyzing complex mixtures by 1D and 2D NMR.

<figure><img src="image_11.png" /></figure>

**Serge Akoka** is a Full Professor of analytical chemistry at Nantes Université (France). His current research interests are (i) isotopic analysis by NMR at natural abundance and (ii) methodological developments in quantitative 1D and 2D NMR.

<figure><img src="image_12.png" /></figure>

**Gerald Remaud** tenured Full Professor at Nantes Université (France), is a specialist in the development of methods for intramolecular isotope measurements and their applications on authenticity, mainly in food. Previously, he held a number of positions in the industry, following research in Sweden into NMR methods for studying nucleotide chemistry.
<figure><img src="image_14.png" /></figure>

**Michelle M. G. Chartrand** is a Research Council Officer at the National Research Council Canada, specializing in isotope ratio mass spectrometry. Her research interests include the production of certified reference materials and the advancement of carbon isotope delta ratio measurements and calibration techniques.

<figure><img src="image_15.png" /></figure>

**Zoltán Mester** is a Principal Reach Officer at the National Research Council Canada. He conducts research in analytical mass spectrometry and certified reference materials. He is the Chair of the Isotope Ratio Working Group of the Consultative Committee for Amount of Substance: Metrology in Chemistry and Biology (CCQM), home of the SI unit, mole.

<figure><img src="image_16.png" /></figure>

**Juris Meija** is a Senior Research Officer at the National Research Council Canada. His work entails research to help improve the reliability of chemical measurements through the development of certified reference materials, evaluation of measurement uncertainty, and better understanding of isotope ratio measurements.
<figure><img src="image_13.png" /></figure>

**Michelle M. G. Chartrand** is a Research Council Officer at the National Research Council Canada, specializing in isotope ratio mass spectrometry. Her research interests include the production of certified reference materials and the advancement of carbon isotope delta ratio measurements and calibration techniques.

<figure><img src="image_14.png" /></figure>

**Zoltán Mester** is a Principal Reach Officer at the National Research Council Canada. He conducts research in analytical mass spectrometry and certified reference materials. He is the Chair of the Isotope Ratio Working Group of the Consultative Committee for Amount of Substance: Metrology in Chemistry and Biology (CCQM), home of the SI unit, mole.

<figure><img src="image_15.png" /></figure>

**Juris Meija** is a Senior Research Officer at the National Research Council Canada. His work entails research to help improve the reliability of chemical measurements through the development of certified reference materials, evaluation of measurement uncertainty, and better understanding of isotope ratio measurements.