Giammarco
Nalin
a,
Kilian
Fehre
a,
Florian
Trinter
bc,
Nikolay M.
Novikovskiy
de,
Nils
Anders
a,
Daniel
Trabert
a,
Sven
Grundmann
a,
Max
Kircher
a,
Arnab
Khan
a,
Raghu
Tomar
a,
Max
Hofmann
a,
Markus
Waitz
a,
Isabel
Vela-Pérez
a,
Gregor
Kastirke
a,
Juliane
Siebert
a,
Dimitrios
Tsitsonis
a,
Hironobu
Fukuzawa
f,
Kiyoshi
Ueda
f,
Joshua B.
Williams
g,
Denis
Kargin
h,
Martin
Maurer
h,
Catmarna
Küstner-Wetekam
d,
Lutz
Marder
d,
Johannes
Viehmann
d,
André
Knie
d,
Till
Jahnke
ai,
Markus
Ilchen
bdi,
Reinhard
Dörner
a,
Rudolf
Pietschnig
h,
Philipp V.
Demekhin
*d and
Markus S.
Schöffler
a
aInstitut für Kernphysik, Goethe-Universität, Max-von-Laue-Straße 1, 60438 Frankfurt am Main, Germany. E-mail: nalin@atom.uni-frankfurt.de; fehre@atom.uni-frankfurt.de; schoeffler@atom.uni-frankfurt.de
bDeutsches Elektronen-Synchrotron (DESY), Notkestraße 85, 22607 Hamburg, Germany
cMolecular Physics, Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195, Berlin, Germany
dInstitut für Physik und CINSaT, Universität Kassel, Heinrich-Plett-Straße 40, 34132 Kassel, Germany. E-mail: demekhin@physik.uni-kassel.de
eInstitute of Physics, Southern Federal University, 344090 Rostov-on-Don, Russia
fInstitute of Multidisciplinary Research for Advanced Materials, Tohoku University, Sendai 980-8577, Japan
gDepartment of Physics, University of Nevada, Reno, Nevada 89557, USA
hInstitut für Chemie und CINSaT, Universität Kassel, Heinrich-Plett-Straße 40, 34132 Kassel, Germany. E-mail: pietschnig@uni-kassel.de
iEuropean XFEL GmbH, Holzkoppel 4, 22869 Schenefeld, Germany. E-mail: markus.ilchen@xfel.eu
First published on 15th July 2021
The photoelectron circular dichroism (PECD) of the O 1s-photoelectrons of trifluoromethyloxirane (TFMOx) is studied experimentally and theoretically for different photoelectron kinetic energies. The experiments were performed employing circularly polarized synchrotron radiation and coincident electron and fragment ion detection using cold target recoil ion momentum spectroscopy. The corresponding calculations were performed by means of the single center method within the relaxed-core Hartree–Fock approximation. We concentrate on the energy dependence of the differential PECD of uniaxially oriented TFMOx molecules, which is accessible through the employed coincident detection. We also compare the results for the differential PECD of TFMOx to those obtained for the equivalent fragmentation channel and similar photoelectron kinetic energy of methyloxirane (MOx), studied in our previous work. Thereby, we investigate the influence of the substitution of the methyl group by the trifluoromethyl group at the chiral center on the molecular chiral response. Finally, the presently obtained angular distribution parameters are compared to those available in the literature.
PECD is a universal chiroptical effect occurring in all regimes,7 from one-photon to multiphoton ionization8,9 and strong-field ionization,7,10,11 with a majority of experiments, so far, reporting on the ionization of valence electrons.12–16 PECD shows a complex dependence on the electronic17 and vibronic18 configuration, as well as the conformation,19,20 dimerization,21 clustering22 of molecules, and the energy of the emitted photoelectrons.23 The typical strength of PECD is on the order of a few percent. Its unique sensitivity has been exploited to study ultrafast molecular dynamics,24,25 especially time-resolved molecular relaxation,26 and it has been employed for real-time determination of the enantiomeric excess in racemic mixtures.27
Several experiments reported on PECD after core ionization of chiral molecules.13,23,28–36 In our recent work,32 we performed a differential study of PECD after O 1s-photoionization of uniaxially oriented methyloxirane molecules (MOx, C3H6O). To this end, the O 1s-photoelectrons were detected as coincident with two charged fragments. The molecular fragmentation axis was assumed to coincide with the emission direction of those ionic fragments relying on the validity of the so-called axial-recoil approximation.37 This assumption is often well justified for cases of rapid Auger decay of an inner-shell vacancy and subsequent Coulomb explosion of a molecular dication. Orientation of achiral molecules in space results in circular dichroism in the angular distribution (CDAD38,39) and in so-called apparent40 PECD. Such extrinsic angular asymmetries caused by the molecular orientation superimpose with the intrinsic asymmetries in the photoelectron emission from chiral molecules. It was found that a selection of certain fragmentation directions of MOx with respect to the light propagation direction yields an increase of the PECD by a factor of 10, as compared to randomly oriented molecules.
In the present work, we extend our previous study32 to a close relative of MOx, the trifluoromethyloxirane (TFMOx, C3H3F3O) molecule, in which all three hydrogen atoms of the methyl group are replaced by much heavier fluorine atoms. For further discussion, it is important to stress that by definition of the Cahn–Ingold–Prelog (CIP41) rules, this substitution exchanges the assigned handedness (S- to R-enantiomer and vice versa) despite the unchanged connectivity of the molecular backbone and the stereo-center in MOx and TFMOx (see Fig. 1). We investigate here the O 1s-photoelectrons of TFMOx for different kinetic energies. One of these energies coincides with that of the study on MOx presented in ref. 32. We selected, furthermore, for our investigation of TFMOx, one of the fragmentation channels which was studied previously for MOx.
Fig. 1 Top: Structure of S- and R-methyloxirane (C3H6O, MOx). Bottom: Structure of R- and S-trifluoromethyloxirane (C3H3F3O, TFMOx). Note that substitution of the methyl group (CH3) with the trifluoromethyl group (CF3) at the chiral center causes renaming of the respective R- and S-enantiomers according to the CIP rule.41 |
Fig. 2 Ratio of enantiomers present in the samples of racemic RS-TFMOx (a) and its S- and R-enantiomers (b and c, respectively) as experimentally determined using HPLC. |
The electron energy calibration was carried out using a series of measurements on argon atoms, ionizing both, the 2p1/2 and 2p3/2 shells. From this series, a calibration function was created which links the measured electron positions-of-impact and times-of-flight with the photoelectron energy. The five photon energies for the TFMOx measurement (541.5, 542.5, 544.5, 546.5, and 550.1 eV) were chosen such that bands of the O 1s-photoelectrons peaked at kinetic energies of 3.1, 4.1, 6.1, 8.1, and 11.7 eV are created. The energy-calibration procedure covered this range of electron kinetic energies. The kinetic energy intervals of ±0.5 eV (within the FWHM = 1.0 eV of the photoelectron bands) were used in the subsequent analysis of the chiral asymmetry. The absolute and relative spatial orientation of the electron and ion detectors was checked comparing the molecular-frame photoelectron angular distributions of N2 K-shell electrons measured during the same beamtime to those reported in ref. 39. The impact positions and times-of-flight of the particles were measured using position- and time-sensitive microchannel plate (MCP) detectors, using delay-line anodes for position read-out (Roentdek GmbH). A HEX90 anode (∅ 90 mm) equipped with an 80 mm diameter MCP49 was used for the electron, and a HEX125L anode (∅ 125 mm) equipped with an efficiency-enhanced 75 mm diameter funnel MCP (Hamamatsu)50 for the ion detection. The momentum vector of each fragment at the instant of ionization is derived from their time-of-flight and impact position on the detector. Data were recorded with typical rates of 6 kHz on the ion and 10 kHz on the electron detector. The combined CF3+(m/z = 69)−C2Hi=1,2,3O+(m/z = 41, 42, 43) breakup channel contains roughly 2.5 × 105 valid events for each enantiomer, light helicity, and photon energy.
In order to cross-check the experimental results presented here, we performed a second short experiment at the permanent COLTRIMS end-station of the PIPE instrument, located at beamline P04 of the synchrotron PETRA III (DESY, Hamburg, Germany).51 In this case, hexagonal delay-line detectors with an active diameter of 80 mm (HEX90 anodes) were used for both, ion and electron detection. The ion arm of that spectrometer was a 17 cm long acceleration region without any meshes, while the electron arm was built in a time-focussing geometry.52 The dichroic parameter β1 was retrieved from the measured data by integrating over all possible fragmentation channels of TFMOx. Therefore, it was sufficient to record datasets for each polarization direction and photon energy for only 20 minutes. The data were recorded at detection rates of 15 kHz for the electron detector. The photoelectron-energy calibrations as well as detector-orientation calibrations were similarly performed as described above for the measurements at synchrotron SOLEIL.
(1) |
(2) |
Fig. 3 Retrieval of the anisotropy and dichroic parameters from the measured PECD of R-TFMOx recorded at the O 1s-photoelectron energy of 11.7 eV. A fit (red line) of the experimental data (circles with error bars) using eqn (2) yielded β1 = 0.76% ± 0.06% and β2 = 1.14 ± 0.15. |
In order to check for possible systematic experimental errors, various subsets of the experimental data were considered separately for this analysis. For example, the high electric fields of the spectrometer degraded the time-of-flight resolution and thus the resolution of the momentum component parallel to the spectrometer axis. As one cross-check, only the central part of the photoelectron time-of-flight distribution was examined and analyzed as a subset. Furthermore, in order to exclude possible influences of inhomogeneous detection probabilities, the electron detector was divided into its left and right half along the light propagation direction, and each subset was analyzed separately. The influence of possible ion feedback from the ion detector was investigated with the help of applying different conditions on the ion time-of-flight and by a selection of different ranges of sums of momenta for various fragmentation channels. All these subsets of the data yielded (within statistical errors) the same angular distribution parameters β1 and β2.
The calculations were performed for different orientations of the molecular frame (MF) with respect to the laboratory frame (LF, ZLF is defined by the propagation direction of the circularly polarized light), which is given by the two Euler angles α and β. Note, that the third angle γ (which defines a rotation around the laboratory ZLF axis) is irrelevant in case circularly polarized light is used. In the experiment, a breakup into two molecular fragments was examined. Thus, only a single fragmentation axis was determined. Therefore, the corresponding differential PECD(θ, β) was obtained by numerical integration of the computational results over the orientation angle α. Here, the orientation interval α ∈ [0, 2π) was covered in steps of Δα = 0.025π. Since the fragmentation dynamics can be rather complicated (and not all of the particles are detected in coincidence), it is not straightforward to relate the fragmentation axis obtained in the experiment from the asymptotic momenta of the two detected ions to a given molecular orientation at the instant of ionization. Therefore, just as in our previous work on uniaxially oriented MOx,32 we used the molecular orientation axis at the instant of the photoionization process as a free parameter in the calculations and searched for the best visual correspondence between the computed and measured differential PECD. Details of this procedure can be found in the supplemental material document of our previous work.32
We first consider the photoelectron kinetic energy of 11.7 eV, which is very similar to that used in ref. 32 for MOx. The measured differential PECDs of uniaxially oriented enantiomers of TFMOx are depicted in the right column of Fig. 4 in panels (b) and (d). For some molecular orientations (i.e., some spatial orientations of the molecular fragmentation), the differential PECD reaches about 20%. Our measurements thus confirm the conclusion of ref. 32 that fixing already one axis of a chiral molecule in space enhances the chiral asymmetry by about an order of magnitude (as compared to about 1% observed for the randomly oriented molecules shown in Fig. 3). The asymmetry switches its sign upon exchanging of the enantiomers, which confirms the chiral origin of the observed effect. This normalized difference also obeys the analytically derived asymmetry equation:32 PECD(π − θ, π − β) = −PECD(θ, β), which implies a flip of its sign upon rotation of this two-dimensional map by 180° around the normal to the picture plane. Please note that this equivalence of changing the sign of the light helicity and exchanging the enantiomer holds only for enantiomeric pure samples and perfectly circularly polarized light (see ref. 10 for more details).
Fig. 4 Measured differential PECD of the uniaxially oriented enantiomers of MOx (left column, data from ref. 32) and TFMOx (right column, present work) as a function of the photoelectron emission angle θ and the molecular orientation angle β. The PECD maps were obtained for O 1s-photoionization yielding similar electron kinetic energies (11.5 eV for MOx and 11.7 eV for TFMOx). Conceptually the same fragmentation channels are chosen, i.e., CH3+(m/z = 15)−C2H2O+ (m/z = 42) for MOx and CF3+ (m/z = 69)−C2Hi=1,2,3O+(m/z = 41, 42, 43) for TFMOx. According to the CIP rule, opposite enantiomers of these two chiral molecules with identical geometrical structure are compared in each row (see Fig. 1 for details). |
For comparison, panels (a) and (c) in the left column of Fig. 4 show the differential PECD of MOx from ref. 32. One can see that for equal fragmentation channels and similar photoelectron kinetic energies the differential PECDs observed for these two molecules are very similar in the strengths and structures. It is surprising that exchanging the three hydrogen atoms in the methyl group by the much heavier scatterers (fluorine atoms) does not significantly change the PECD landscape, observed at this photoelectron kinetic energy. A similar observation was made in a comparative study17 of the valence-shell PECD of MOx and TFMOx, where noticeable differences between the respective chiral asymmetries were found only for kinetic energies below 4 eV. In each row of this figure, we compare the data for MOx and TFMOx molecules of similar geometrical structures (therefore, opposite R- and S-enantiomers). One can see that the differential PECDs of the two molecules of similar geometrical structures have opposite-sign landscapes, whereas for identical R- or S-enantiomers of both molecules the measured sign landscapes are very similar. This fact confirms that the inner-shell PECD effect emerges as a result of a complex multiple scattering of the emitted photoelectron wave by the molecular potential. Indeed, similar stereo-descriptors of enantiomers assume a similar circular arrangement of the stronger-to-weaker scatterers, which couple the rotational motion of a photoelectron, driven by the electric field vector of circularly polarized light, with its translational motion (like a nut on a bolt23,32).
Fig. 6 Dichroic parameter β1 (a) and the anisotropy parameter β2 (b) as given by eqn (1) for the O 1s-photoelectrons of R-TFMOx. Circles: experimental results recorded at the BL13-2 beamline of SSRL using a VMI setup, as reported in ref. 33. Stars: β1 obtained by the present reanalysis of the raw data underlying the values shown by circles from ref. 33 (see Section 3.3). Triangles: results of the present experiment performed at beamline SEXTANTS of SOLEIL using COLTRIMS. Squares: corresponding β1 values measured at beamline P04 of PETRA III. Dashed curve: frozen-core Hartree–Fock calculations performed with the SC method as reported in ref. 33. Solid curve: present relaxed-core Hartree–Fock calculations performed using the SC method. |
In order to calculate the differential PECD of uniaxially oriented TFMOx, we first assumed validity of the axial-recoil approximation,37 for which the momentum difference of two fragments represents an orientation of the molecular axis in space at the instant of ionization (i.e., the molecule does not rotate significantly during the Auger decay and subsequent fragmentation). This approximation worked well in our previous study.32 There, we found that—for the single-bond breaking with the loss of the methyl group—the optimal molecular axis, which provides the best description of the experiment, connects the oxygen atom with the carbon of the methyl group (see the ESI of ref. 32). Starting with this assumption, we searched for the optimal molecular fragmentation axis, which provides the best match between theory and experiment for the differential PECDs of uniaxially oriented TFMOx for all five electron energies.
As a result of our search, we find that the best visual correspondence between the measured and computed differential PECD is provided by a fragmentation axis which connects the oxygen atom with the CH2-group on the oxirane ring. It should be stressed that all five photoelectron energies yielded very similar optimal fragmentation axes which coincide within a solid angle of about 0.2 sr. The fact, that the axis obtained this way does not connect two charged fragments (the trifluoromethyl group and the oxirane ring) but lies almost in the plane of the oxirane ring, indicates a breakdown of the axial-recoil approximation for the TFMOx molecule. Indeed, a considerably larger mass of the trifluoromethyl group (as compared to the methyl group in MOx) and a much richer electronic structure of fluorine atoms (as compared to hydrogens) may lead to the fact that the centers of charge of the fragments do not coincide with their centers of mass. In this case, an intricate rotation of two fragments with respect to each other takes place in the course of a considerably slower dissociation of TFMOx (as compared to the same fragmentation channel of MOx with the much lighter methyl group), and, as a consequence, the differences of their asymptotic momenta do not represent the molecular orientation at the instant of photoionization.
The results of our calculations of the differential PECD of R-TFMOx, obtained for the optimal fragmentation axis, are depicted in the right column of Fig. 5. The calculations provide a semi-quantitative theoretical description of the measured differential PECD. Indeed, the overall landscapes and signs of the experimental differential PECDs, including the arrangement of stripes and islands discussed above, are reproduced by the theory. Detailed distributions of the experimental and theoretical asymmetries in their forms and strengths are, however, somewhat different. This leads to the fact that averaging the theoretical and experimental signals over molecular orientations yields different β1 values. This fact can be seen, e.g., for the photoelectron energies of 8.1 and 11.7 eV in the forward emission direction. In particular, in the experiment, molecular orientations yielding positive chiral asymmetries dominate over those which yield a negative PECD, while this situation is opposite in the theory [cf., in Fig. 5 panels (g) and (h), and separately panels (i) and (j), for the emission angles around θ = 0° and orientation angles around β = 90° with highest contribution weights]. This difference is reflected in the opposite signs of the computed and measured dichroic parameters at 8.1 and 11.7 eV [see panel (a) of Fig. 6 and its discussion in Section 3.3].
In order to understand the considerable deviation between the presently and previously measured dichroic parameters β1 seen in panel (a) of Fig. 6 in the energy range of 3–6 eV, in a joint effort with the authors of ref. 33, we reanalyzed the original Stanford Synchrotron Radiation Lightsource (SSRL) raw data recorded at the BL13-2 beamline with a VMI setup. The reanalysis of raw images was performed independently for every run of a certain photon energy, helicity of the light, and enantiomer, using the following steps. The intensity in the surrounding of the momentum spheres was subtracted as the background. Then, the momentum sphere had to be centered. It turned out that the center varied from run to run on the level of a few percent of the radius of the momentum sphere giving rise to systematic errors if not corrected for. This shift had not been accurately corrected for in the original analysis in ref. 33. An Abel transformation was applied to the centered distributions followed by a transformation into polar coordinates. Finally, only those electrons, which have the proper kinetic energy, are selected for further analysis. This is done by gating on a certain thickness of the momentum sphere of the Abel-inversed data. The resulting subset of data (intensity as function of the emission angle with respect to the light propagation axis) is treated analogous to the COLTRIMS data (see Section 2.3). Results of this reanalysis of the SSRL data from ref. 33 are shown in panel (a) of Fig. 6 by pink stars (note that the reanalysis was performed only in the relevant energy range of the present experiment). The specified uncertainties correspond to the fit error from the extraction of the β1 parameter. A subdivision of each run from the SSRL measurement into two independent data sets confirmed both, the retrieved values of the β1 parameter and the respective uncertainties. As one can see from this figure, the dichroic parameters determined in all three experiments are now consistent.
In order to find a possible reason for the disagreement between the previously computed and all measured dichroic parameters in panel (a) of Fig. 6, we cross-checked the computational results of ref. 33. First of all, we reexamined the equilibrium geometry of the molecule obtained in ref. 33 on the (2,2)-CASSCF/6-31+G(d,p) level of theory and used there in the subsequent calculations. For this purpose, we now performed its optimization with DFT methods using the B3LYP hybrid functional with a 6-311+G(d,p) basis set. The minimal-energy structure was confirmed by the absence of imaginary frequencies. It turned out that both equilibrium structures differ only negligibly and the angular distribution parameters recomputed here in the FCHF approximation basically confirm those computational results of ref. 33 (not shown here for brevity). As a next step, we recomputed the angular distribution parameters of the O 1s-photoelectrons of TFMOx using the RCHF approximation. It is known to be more accurate for core ionization than the FCHF approximation used in the previous work.33 The results of these calculations are shown in Fig. 6 by red solid curves. As one can see, the effect of the core relaxation caused by the created inner-shell vacancy does not significantly change the energy dependence of the computed angular distribution parameters, and the presently computed β1 and β2 parameters agree well with the previously computed and measured data from ref. 33 (blue dashed curve and yellow circles). For the dichroic parameter β1, which is of primary interest here, relaxation of the core has just somewhat reduced its maximal and minimal values. Possible reasons for the remaining disagreement between the theoretical and experimental dichroic parameters are discussed in Section 4.
The molecular fragmentation axis, found theoretically for this breakup channel, does not connect the CF3+ and C2H3O+ fragments, which is a clear indication that for larger molecules, the connection between the orientation of the molecule at the instant of photoionization and the asymptotic momenta of the detected fragments have a rather intricate interconnection. Such theoretically found optimal axes of molecular fragmentation might in principle be used to estimate times for a relative rotation of the fragments during the considered breakup. The computed differential PECDs reproduce the experimental results for all photoelectron kinetic energies on a semi-quantitative level. Having rather similar landscapes, the computed and measured differential PECDs differ in their forms and maximal asymmetry strengths. As a consequence, the averaging of those rather strong asymmetries over all molecular orientations yields different small values of the theoretical and experimental dichroic parameters β1. We performed a reanalysis of the raw data from the experiment published in ref. 33 leading to a revision of some of the data points previously reported. These corrected previous experimental results agree with new results from two different experiments performed at two light sources in the present work.
Being very similar in their geometrical structures, S/R-MOx and R/S-TFMOx molecules are very different in their chemical and, thus, electronic structures. This could be a reason for a very different degrees of agreement between theory and experiment in our earlier study of MOx32 and in the present and earlier studies of TFMOx.33 Because of considerably more electrons in TFMOx, the respective calculations seem to be much more sensitive to minor details of the molecular potential, especially for slow photoelectrons. Experimentally, very different masses of constituent parts in TFMOx make an unambiguous assignment of a molecular orientation axis to a certain fragmentation channel impossible. In order to clarify the remaining disagreement between the theoretical and experimental dichroic parameters β1, the O 1s-photoionization of TFMOx needs to be reinvestigated with alternative theoretical methods (e.g., TDDFT56 or CMS-Xα57). Ref. 28 reports a comprehensive theoretical study of the photoelectron angular distributions from chiral derivatives of oxirane by the former theoretical approach (note that TFMOx was not studied there). As one can learn from Fig. 8 of ref. 28, dichroic parameters computed for O 1s-photoelectrons of all considered oxirane derivatives exhibit a clear enhancement to about 5–8% within the photoelectron kinetic energy interval of 2–8 eV, supporting our theoretical findings on TFMOx. Taking also into account the robustness of the present and previous33 calculations performed by the SC method, a reason for the remaining disagreement between the theoretical and experimental dichroic parameters β1 can lie beyond the (fixed-nuclei one-particle approximation) electronic structure calculations performed here. For instance, a presence of a shape resonance in the O 1s-photoionization of TFMOx just above the ionization threshold may reduce the accuracy of the one-particle Hartree–Fock calculations performed here. In addition, slow photoelectrons may share their chiral asymmetry with the Auger electrons via post-collision interactions (PCI) neglected in the present one-particle approach.
This journal is © the Owner Societies 2021 |