Preparation of Tautomer-Pure Molecular Beams by Electrostatic Deflection

Tautomers are ubiquitous throughout chemistry and typically considered inseparable in solution. Yet (bio)chemical activity is highly tautomer-specific, with common examples being the amino and nucleic acids. While tautomers exist in an equilibrium in solution, in the cold environment of a molecular beam the barrier to tautomerization is typically much too high for interconversion, and tautomers can be considered separate species. Here we demonstrate the first separation of tautomers within a molecular beam and the production of tautomerically pure gas-phase samples. We show this for the 2-pyridone/2-hydroxypyridine system, an important structural motif in both uracil and cytosine. Spatial separation of the tautomers is achieved via electrostatic deflection in strong inhomogeneous fields. We furthermore collect tautomer-resolved photoelectron spectra using femtosecond multiphoton ionization. This paves the way for studying the structure–function–dynamic relationship on the level of individual tautomers, using approaches that typically lack the resolution to do so, such as ultrafast dynamics experiments.

Tautomerization is ubiquitous across organic and biochemistry and commonly associated with any polar molecule containing acidic functional groups.Tautomers and their different tautomer-specific chemical functionalities play a key role in many biochemical processes.From enzyme catalysis [1] and ligand binding [2], to DNA mutation and protein misfolding [3,4] and photo-stability of nucleic acids.While in solution the tautomer equilibrium can be heavily influenced by the solvent and solvent properties, most tautomeric systems are considered non-separable in solution.The situation in typical gas-phase experiments is very different.
Here the extremely cold and collision-free environment of a supersonically expanded molecular beam usually prevents any tautomerization from occurring.However, since the transfer into the gas phase is so rapid, tautomer distributions typically remain nearly unchanged, with structures frozen into their respective local minima, rather than relaxing into the global minimum of the most stable tautomer.Hence gas-phase experiments still face the challenge of molecular samples containing multiple tautomers.While frequently these can subsequently be distinguished using high-resolution spectroscopy techniques, they can not be separated.This implies that properties such as relaxation dynamics can only be studied with tautomer resolution if the respective absorption wavelengths are sufficiently different for the tautomers.In the case of the 2-pyridone /2-hydroxypyridine tautomer system considered here, for example, the dynamics have been studied at specific excitation wavelength [5,6], where it is possible to selectively excite or ionize only a single tautomer.Studies of the relaxation dynamics were both tautomers absorb, e.g. at wavelengths deeper in the UV such as 200 nm [7], are not feasible with this approach.A case in point here are the still much debated electronic relaxation dynamics of cytosine, of crucial relevance for understanding DNA photostability.Different studies have reported vastly different excited-state lifetimes and differing numbers of relaxation channels [8][9][10].These discrepancies can, at least partly, be attributed to different tautomer contributions within the gas-phase samples [11].
Similarly many modern gas-phase methodologies cannot distinguish different tautomer populations.In particular, experiments utilizing ultrafast laser pulses, such as attosecond charge-migration measurements [12] or reaction pathway observations [13], inherently lack the spectral resolution to distinguish tautomers.Moreover, molecular collision studies typically require pure samples of the collision partners [14,15], as do any imaging experiments based on electron or x-ray diffraction, which now have the potential for atomic temporal and spatial resolution [16][17][18].We show here spatial separation of tautomers in the gas-phase with the production of a tautomerically-pure cold molecular beam, by utilizing the established electrostatic deflection technique [19][20][21].In particular we demonstrate this approach for the 2-hydroxypyridine (enol or lactim tautomer)/2-pyridone (keto or lactam tautomer) system (Figure 1), a structural motif for both uracil and cytosine and a common model system for DNA tautomerization.
Both tautomers are stable in the gas phase, with a barrier of 34-38 kcal/mol between them.
The 2-hydroxypyridine tautomer is slightly lower in energy by about 300 cm −1 , which leads to an expected tautomer distribution of about 3:1 (2-hydroxypyridine:2-pyridone) [6] for pure samples at room temperature, whereas in polar solvents and the crystal phase 2-pyridone dominates [22].Additionally, this system is known to form very stable dimers in the gas phase due to strong N-H and O-H hydrogen bonding interactions, making it an ideal model system for DNA hydrogen bonding.
In order to control and separate the two tautomers we use the electrostatic deflection approach [19][20][21].This exploits the difference in the permanent dipole moments of the tautomers, leading to a species-specific Stark interaction with an electrostatic field [21], as shown in Figure 1(b).This shows the calculated Stark energy shift for the rotational levels J = 0 − 2 for both tautomers, calculated using the freely available CMIstark software package [23].It clearly highlights that, due to their distinct permanent electric dipole moments, both tautomers exhibit significantly different Stark shifts.We utilize this by passing the tautomers through an electrostatic deflector with a strong inhomogeneous electric field inside [19], leading to the exertion of a tautomer-specific force and hence acceleration.
This transverse acceleration then leads to a spatial separation of species based on their effective dipole-moment to mass ratio [21].This approach has previously been successfully demonstrated for the separation of, amongst others, individual quantum states of small (triatomic) molecules [24,25], rotational isomers [20,26,27] and clusters [28][29][30].Here we use this approach to separate the two tautomers.We moreover collect photoelectron spectra of the tautomer-separated samples, allowing us to extract and assign tautomer-specific photoelectron spectra using femtosecond multiphoton ionization (fs-MPI).Our separation measurements are supported by trajectory simulations and all spectral assignments are confirmed by quantum-chemical calculations.
Typical mass spectra measured in the molecular beam following fs-MPI are shown in Figure 2(a).These are dominated by the parent ion at m/z=95, although some fragmentation is observed with significant channels at m/z=28, 39 and 67.In panel (b) we show a reference mass spectrum from NIST obtained using electron impact ionization [31], which exhibits identical fragmentation channels, with an overall larger degree of fragmentation.
This confirms our molecular assignment and highlights the relative softness of the fs-MPI process [32].Pyridone is known to form very stable dimers, and these were also present in our mass spectrum.Even under supersonic expansion conditions optimized to minimize dimer formation, we observed a signal in the dimer channel with an intensity of around 5% that of the monomer, as shown in the inset of Figure 2. 2-pyridone can form three different dimers, a homomeric dimer consisting of two 2-hydroxypyridine units, a homomeric dimer of two 2-pyridone units and a mixed dimer.At a temperature of 400 K, the expected relative populations of the dimers are 1.3%, 98.4% and 0.4% for the (2-hyroxypyridine) 2 (enol-enol ), (2-pyridone) 2 (keto-keto) and mixed (keto-enol ) dimers [33], respectively.Hence the 2-pyridone dimers are expected to dominate in the molecular beam, as experimentally shown by Held et al [34].This also follows from computational optimization of the dimer structure, as the mixed dimer optimizes to the structure of the (2-pyridone) 2 dimers through hydrogen shuttling, indicating that this dimer is the lowest energy configuration.However, other studies have previously observed both dimer systems [35,36], without quoting relative abundances.The dimers consisting of a single tautomer do not exhibit a dipole moment due to their inherent symmetry.The mixed dimer is expected to have a dipole moment of about 5 D, however, as outlined above, it is not expected to have any significant population in the molecular beam.
In order to separate the tautomers, the molecular beam was passed through an electrostatic deflector, with a potential difference of 19 kV applied across the electrodes (see S.I. for details).
The resulting deflection of the molecular beam was monitored by translating the tightly focused probe laser beam and recording mass spectra at various positions, yielding a spatial density profile.This is shown in Figure 3 for the m/z=95 channel for both the field-free case (0 kV, black data points) and with the strong field present (19 kV, red data points).In the field-free case a symmetric distribution centered at 0 mm was observed, which in the presence of the strong field deflected to the right and developed a noticeable edge on the more deflected side, around 1.5 mm.Since this data is extracted from fs-MPI and time-of-flight mass spectrometry, it is not tautomer-specific.In order to understand the relative tautomer distributions within the beam, we have performed trajectory simulations for both tautomers under realistic experimental conditions [21].These simulations were then fitted to the experimental data, yielding an initial rotational temperature of 3.4 K and a 2-pyridone fraction of 16% in the molecular beam, as discussed in the supplementary information.Simulated spatial density profiles for the individual tautomers are shown as shaded areas in Figure 3, whilst the sum of the two contributions is shown as a solid line.The simulation results are in excellent agreement with the experimental data points, and show that at spatial positions > 1.2 mm a pure beam of 2-pyridone is produced.We note that the deflected beam is also free from contamination by the dimer system, as expected based on the known dimer distributions and dipole moments.This is confirmed by the mass spectrum around m/z=190, as shown in the inset of Figure 2, with no detectable dimer signal in the deflected beam (see SI for further details).
In order to further confirm the tautomer separation in the deflected beam, we measured photoelectron images in the center of the beam and in the deflected edge, at the positions indicated by the dashed lines in Figure 3, which correspond to 0.1 and 1.4 mm above the center of the molecular beam.These measurements were again conducted using fs-MPI.The ionization energies for 2-pyridone and 2-hydroxypyridine have previously been determined as 8.93 eV and 8.45 eV, respectively [37].Thus necessitating 3 photons at 400 nm (total photon energy 9.3 eV) for ionization.
Collected photoelectron images and spectra are shown in Figure 4. We note that all images show a sharp signal in the center of the image, corresponding to 0 eV electrons, that originate from residual water.In panel (a) we show the data collected at 0.1 mm above the center of the molecular beam, and therefore containing electrons originating from both tautomers, as well as contributions from the dimer systems.The photoelectron spectrum was extracted using Abel inversion [38,39], and shows several distinct peaks in particular in the range 0-2 eV of electron kinetic energy.The corresponding image and photoelectron spectrum for the deflected beam are shown in Figure 4(b), and show only a subset of the features present in (a).We assign this spectrum collected in the deflected beam to the 2-pyridone keto tautomer.
The known ionization energy should yield a maximum electron kinetic energy of 0.9 eV, agreeing well with our data when taking into account the bandwidth of a 3-photon fs-MPI process.We have furthermore calculated the expected vertical ionization energy for the formation of the X 0 0 cation ground state, which yielded an electron kinetic energy of 0.73 eV (indicated by the vertical bar in Figure 4(b)), in good agreement with the experimental data.The experimental spectrum contains 3 closely spaced features approximately 0.11 eV (890 cm −1 ) and 0.20 eV (1610 cm −1 ) apart, most likely arising from vibrational excitation of in-plane bending modes upon ionization [40].
In Figure 4(c) we show a photoelectron spectrum of the mixed beam (a) with the pure-keto spectrum (b) subtracted after normalization on the 2-pyridone signal.The resulting difference spectrum contains two major features.We assign the low energy feature to formation of the X 0 0 cation ground state in the 2-hydroxypyridine tautomer.The calculated vertical ionization energy for this yields electrons with kinetic energies of 0.23 eV (blue vertical bar in Figure 4(c)), whilst the published adiabatic ionization energy corresponds to an energetic cutoff at 0.4 eV, both in agreement with our data.This energetic cut-off also implies that the high energy feature observed cannot correspond to the enol tautomer.Instead we attribute this feature to the dimers present in the molecular beam.We have calculated the vertical ionization energies for the formation of the cation ground state for all the three potential dimer systems (yielding expect electron kinetic energies of enol-enol, 0.83 eV; keto-enol, 1.15 eV and keto-keto, 1.34 eV ), as indicated by the grey vertical bars in Figure 4(c).The major peak observed energetically agrees with the calculated vertical ionization energies of both the keto-enol and keto-keto geometries.However, given that we did not observe any deflection in the dimer channel, we assign this peak to the keto-keto dimer structure, in agreement with previous experimental studies [34] and the relative dimer stabilities, that predict > 98% of dimers to be in the keto-keto geometry [33].From the collected spectrum we obtain values of the vertical and adiabatic ionization energy of the dimer of 8.1 eV and 7.7 eV, respectively.Both the feature attributed to the enol tautomer, as well as that due to the dimer, are absent from the photoelectron spectrum recorded in the deflected molecular beam (Figure 4(b)), confirming the creation of a pure tautomer-selected sample.
In summary, we demonstrated here the separation of tautomers in the gas phase, utilizing electrostatic deflection, for the important 2-pyridone/2-hydroxypyridine system.This enabled the creation of a pure sample of the keto tautomer, which furthermore was free from contamination by the strongly-bound dimers.Experimental results showed excellent agreement with detailed trajectory simulations, and were furthermore corroborated via photoelectron spectroscopy with a universal 400 nm femtosecond multiphoton ionization process.Recorded spectra could be assigned using quantum chemical calculations.The methodology used here is widely applicable and only requires that the molecule of interest can be entrained in a cold molecular beam and that the tautomers have a sufficient difference in dipole moment to mass ratio [21].This makes this approach highly applicable to many important biological building blocks, such as cytosine [41].The spatial dispersion of tautomers can generally add tautomer sensitivity to experimental techniques that otherwise inherently lack the resolution to distinguish them, as we have shown here for femtosecond multiphoton ionization.This paves the way for ultrafast dynamics experiments or time-resolved imaging studies focusing on the biologically relevant tautomers of important biomolecules.

Experimental Methods
All experiments were performed in a custom-build molecular beam setup, equipped with an Even-Lavie valve [42], 30 cm long electrostatic deflector [19][20][21] and velocity-map imaging (VMI) spectrometer [43], further technical details are given in the supporting information.Chemical samples were commercially available (97% purity, Sigma-Aldrich) and used without further purification.Molecules were ionized via femtosecond multiphoton ionization (fs-MPI) using 400 nm laser pulses of 100 fs duration.
Experimental deflection measurements were complemented by trajectory simulations.For these the Stark effect of both isomers was calculated using the freely available CMIstark 9 software package [23].Trajectories of individual quantum states (J = 0 − 15, 2000 trajectories per state) were then propagated through the experimental setup [21].Individual trajectories were combined to a simulated deflection profile through Boltzmann weighting, with the rotational temperature and the tautomer distribution fitted to the experimental data.

FIG. 1 :
FIG. 1: Structures and dipole moments of 2-hydroxypyridine and 2-pyridone (a), and the associated Stark shift as a function of the electrostatic field (b).These are shown for rotational states J = 0 − 2 and clearly show the different extent of the Stark interaction for the two tautomers.

FIG. 2 :
FIG.2: Mass spectrum of 2-hydroxypyridine/2-pyridone collected with 400 nm fs-MPI (a), and a reference spectrum from NIST collected with electron impact ionization (b)[31].The inset shows an enlarged view of the mass region containing dimer signal.In the center of the molecular beam (black trace) the dimer signal amounts to around 5% of the monomer signal, whereas in the deflected beam (red trace) no dimers are observed, see text for further details.

FIG. 3 :
FIG.3: Spatial density profiles of the molecular beam without electrostatic field (black) and with electrostatic field present (red).Measurements are shown as data points with error bars, while solid lines correspond to trajectory simulations.The shaded areas correspond to the tautomer-resolved densities extracted from the simulations.These show that for the deflected beam, at positions > 1.2 mm a pure sample of 2-pyridone is created.

FIG. 4 :
FIG. 4: Photoelectron images and extracted photoelectron spectra in (a) the center of the beam and (b) the deflected edge of the molecular beam.The latter corresponds to a spectrum of pure 2-pyridone.(c) Mixed photoelectron spectrum with the contribution from the 2-pyridone subtracted.This contains features from both 2-hydroxypyridine, as well as from dimers.Calculated vertical ionization energies are shown as solid bars throughout.For the dimer system calculated energies for all three stoichiometries (enol-enol (ee), keto-enol (ke) and keto-keto(kk)) are shown.