D.
Catone
a,
N.
Zema
a,
T.
Prosperi
a,
M.
Stener
bc,
P.
Decleva
bc,
P.
Nitti
b and
S.
Turchini
*a
aIstituto Struttura della Materia-CNR (ISM-CNR), Via del Fosso del Cavaliere 100, 00133, Roma, Italy. E-mail: stefano.turchini@ism.cnr.it
bDipartimento di Scienze Chimiche, Universith di Trieste, Via L. Giorgieri 1, I-34127 Trieste, Italy
cIstituto Officina dei Materiali-CNR (CNR-IOM), Area Science Park Strada Statale 14 km 163, 5 - 34149 Basovizza, Trieste, Italy
First published on 20th January 2022
In this work we address a specific experimental and theoretical question regarding the influence of a conformational population in the modeling of photoelectron circular dichroism (PECD) spectroscopy. In the past two decades, PECD has revealed a rich and complex phenomenology in molecular processes with unprecedented insight, especially in molecular geometry sensitivity. Since the early development of this spectroscopy, theory has pointed out the importance of conformer influence on PECD; in particular, the rotation of methyl groups was surprisingly found to be responsible for strong modulation of the PECD signal. Here, to advance understanding of the effect of rotations, we have chosen to study norcamphor, a single-conformer molecule, as a benchmark for a PECD comparison between experiment and theory at the density functional theory (DFT) and time-dependent density functional theory (TDDFT) levels. The excellent agreement between experimental data and theory for norcamphor sheds light on the influence of rotations and gives a solid explanation for the reasonable qualitative agreement in the PECD of camphor, where three methyl groups are added to the same molecular structure.
In the past two decades, a new approach to the conformer study of chiral molecules has been explored using photoelectron circular dichroism (PECD), which is an extension of classical UV-IR circular dichroism applied to the VUV-soft X-ray range where photoionization is present. The final state is a photoelectron with a defined kinetic energy and momentum, which allows us to measure the circular dichroism as a two-dimensional dispersion, while circular dichroism in absorption is only measured as a function of the photon energy.
PECD has been measured in almost every sample target compatible with photoelectron spectroscopy, effusive and supersonic molecular beams in the gas phase, and thin films. The asymmetry has the same order of magnitude for valence and core PECD, and it is one order of magnitude larger than the absorption one, because the matrix element is at the electric dipole level.10 Moreover, PECD has been measured using photoelectron photoion spectroscopy,11–14 photoelectron–Auger electron coincindence,15 resonant photoemission,16 multiphotonics photoelectron spectroscopy,17 and femtosecond time-resolved pump–probe experiments.18
Since the beginning it has been clear that PECD is extremely sensitive to conformer geometry19,20 although this result has turned out to be ambivalent: on the one hand, PECD has revealed itself to be a powerful form of spectroscopy that is capable of detecting elusive electronic and structural properties; on the other hand, this sensitivity deeply challenges the theoretical modeling. The dependence of PECD on tiny distortions of geometry is attributed to the partial wave interference of the form sin(ηl − ηl′) in the matrix element,21 where ηl are the generalized phase shifts of the asymptotic wave of photoelectrons for a non-central potential and l and l′ are the angular momenta. For a phase shift of the same order of magnitude the sine can rapidly change sign. It is worth noticing that classical photoemission depends on the cosine of the difference of the phase shifts and near zero is almost constant.
We underline that the rich and complex oscillatory behavior stemming from each state of the molecule represents a serious benchmark for theory. To better focus the aim and results of our work, we summarize the results so far achieved in the study of conformers using PECD. A comparison between valence band spectra and theory in the methyl-oxirane molecule showed excellent agreement between LCAO–B spline theory22 and experiment.23,24 In a subsequent theoretical study the focus was the effect of the conformational geometry on PECD. The most dramatic variation is related to the rotation of the methyl group of methyl-oxirane with variation of the intensity and sign.25 The same behavior has been reported even for core levels, where the picture of local character of the initial state suggests a less intense effect. However, work on the conformation of methyl-oxirane confirms the previous result: performing a Boltzmann average over conformations, the obtained dispersion is reasonably close to that of the equilibrium geometry used for the previous work.
At this point the most urgent question is how reliable is the equilibrium geometry in the presence of rotations? Subsequently, if the accuracy of the calculations in the case of the equilibrium geometry is not satisfactory, how can the effect of the conformers be efficiently included? The cases of some organic molecules studied gave us some hints at interpretation. In the case of alaninol26 we tried to interpret the data with the population at room temperature formed by the two conformers belonging to the most representative population at equilibrium. For some states we successfully improved the agreement considering the presence of the second conformer. A PECD study of the valence band of 3-methyl-cyclopentanone27 revealed a change in the PECD spectra upon heating the molecular gas at a temperature of 370 K. On the basis of a two-conformer population analysis we found a good agreement in the comparison between the theoretical predictions and experimental data for the difference in the PECD spectra obtained at 370 K and room temperature (300 K). The agreement in the difference of the spectra taken at the above-mentioned temperatures is better than that of the terms of the difference. Our interpretation was that in this case the subtraction of the PECD taken at different temperatures lifted the contribution of the rotating groups and the result was closer to the difference of the calculations at equilibrium geometry. Similar conformational effects were observed for alanine11,28 and proline.29
We planned an experimentum crucis choosing (+)-(1S,4R)-norcamphor (Fig. 1), a molecule with a single conformational geometry. Norcamphor (2-norbornanone, bicyclo[2.2.1]heptan-2-one) might be thought of as being the structural parent compound of the large bornane family, whose more versatile chiral compound is the naturally occurring bicyclic ketone (+)-(1R,4R)-camphor, which differs from norcamphor by the presence of two methyl groups bonded at the C7 atom and one methyl group substituting the hydrogen bonded to C1, adjacent to the carbonyl group. The multiple-site reactivity of the bornane skeleton and its rigidity has made these molecules highly useful in different fields for the introduction of chiral groups into pharmacologically active molecules30 and in enantioselective hydrogenation.31 Finally, norcamphor has been confirmed as a suitable chemical for rare functionalisations at carbons of the bornane moiety that provide an efficient approach for the preparation of important norcamphor scaffolds.32
The experimental PECD of camphor has been carefully studied33,34 through comparison using two different state-of-the-art theoretical methods multiple scattering-Xa (MS-Xa)35,36 and LCAO B-spline (B-spline). The outcome of both calculation methods is in accord and describes qualitatively the experimental dispersion for the highest occupied molecular orbital (HOMO). Our aim is to find evidence of whether or not the average of the conformational rotamer configurations distorts PECD calculations at the equilibrium geometry. The idea at the heart of the experiment is to study norcamphor to suppress the effects of the rotamers and to critically compare the outcome of the theory with the experimental results.
Finally, it is worth mentioning that vibrations can also affect PECD via ion displacement.
Actually, a dramatic breakdown of the Franck–Condon (FC) approximation has been found for the valence band of the methyl-oxirane molecule.37 Apart the obvious consideration that vibrations are related to different geometrical configurations, the complex nature of the matrix element, together with the consideration of pseudo-asymmetries, creates an outcome of the sum of the contributions that is extremely sensitive to tiny phase change between them. Indeed, the current landscape, supported by theoretical work,38 points out that weak FC transitions associated with small ion displacements are the transition that is more amenable to break the FC approximation. Consequently, since in unresolved PECD the transitions with large FC factors are more weighted in the overall signal, it seems to be safe to interpret in this way the success of vibrationally unrestricted PECD calculations, although it is not possible to exclude violations as a result of more complex behavior.39 Clearly, a final goal for theoretical simulations will include the accurate treatment of both the conformational flexibility and the remaining vibrational modes.
In PECD data analysis each i-th electronic state of the molecule gives rise to a dispersion as a function of the photoelectron kinetic energy. Phenomenologically, PECD is described via three dynamical parameters in the angle-resolved differential cross-section of the photoelectrons. In the following the equation for elliptically polarized light is reported:
(+)-(1S,2R,4R)-endo-2-Norborneol was purchased from Sigma-Aldrich. Optical rotations were measured using a PerkinElmer Model 241 polarimeter. TLC was performed using Polygram® Sil G/UV254 silica gel pre-coated plastic sheets (eluent: light petroleum–ethyl acetate, 4:1). Flash chromatography was run on silica gel 230–400 mesh ASTM (Kieselgel 60, Merck).
(+)-(1S,4R)-Norcamphor was synthetized through the oxidation of (+)-(1S,2R,4R)-endo-2-norborneol following the procedure indicated in the literature.41 (+)-endo-2-Norborneol (1.0 g, 8.9 mmol) was dissolved in 50 mL of anhydrous dichloromethane at 0 °C under Ar, and 3 Å molecular sieves (9 g) and pyridinium chlorochromate (6.0 g, 27.8 mmol) were added and the suspension was vigorously stirred for 3 h. Diethyl ether (100 mL) was added and, after decantation, the solution was filtered through a glass filter filled with silica gel with 13% calcium sulphate (Fluka). The black residue was washed three times with diethyl ether and the combined filtrates were carefully concentrated via rotary evaporation. The resulting (+)-(1S,4R)-2-norbornanone was purified by flash chromatography (eluent: light petroleum:diethyl ether, 9:1), 40% total yield with [α]25D = +26.3 (c 1.0, CHCl3) ([α]25D = + 29.1 (c 1.5, CHCl3),42 [α]25D = + 29.8 (c 1.1, CHCl3),43 [α]25D = +27.2 (c 1.8, CHCl3)44); all spectroscopic data are identical to a commercially available sample of racemic 2-norbornanone.
To acquire data in linearly polarized radiation a commercially supplied (purity 98%, Sigma–Aldrich) racemic mixture of norcamphor was used. The norcamphor was introduced into an effusive source exploiting the room temperature vapor tension.
To ease the comparison between norcamphor and camphor we reversed the sign of the dichroism of (1S,4R)-norcamphor to obtain the PECD of the (1R,4S)-norcamphor enantiomer.
The pressure during the measurements was about 5 × 10−6 mbar. PECD measurements were recorded alternating the helicity of the radiation at a 0.05 Hz frequency. A complete description of the experimental D extraction and analysis and the correction for the ellipticity of the light present at the Circular Polarization beamline are reported elsewhere.26
The B-spline LCAO method was employed to calculate the photoelectron wave function in the continuum.47 This is based on the use of basis functions of the type
χilm = 1/rBi(r)Ylm(θ,φ) |
Ionization energies were separately computed to confirm assignment of the spectrum using the OVGF approach and the Gaussian program.48
The isolated HOMO peak allows us to neglect any mixing of other molecular orbitals in the experimental data analysis and, consequently, to avoid convoluting theoretical data with a poor reliable assumption about the energy position, intensity and width.
Fig. 3 displays the experimental valence band D dispersion of the HOMO of (−)-(1R,4S)-norcamphor as a function of kinetic energy along with the DFT and TDDFT calculated dispersions. The theoretical D parameters are shifted 1.6 eV towards higher kinetic energy values to align the minima, where this phenomenological shift is justified by the attractive behavior of the LB94 functional. The D parameter is always negative, and the shape of the dispersion is a valley with a marked minimum at 9 eV and the full shape presents excellent quantitative agreement for both DFT and TDDFT. In particular, the TDDFT dispersion is closer to the experimental one and displays a more structured behaviour; the calculated dispersion is also in agreement with the fine structure of the experimental dispersion in the range of 12–14 eV. This result is of paramount importance because it straightforwardly addresses the issue we posed at the beginning, i.e., the absence of rotating groups has a remarkable influence on the comparison between calculation at the equilibrium geometry and the experiment.
To better understand this crucial aspect, we compare our results with those obtained in the case of the HOMO of camphor and those reported in the literature.34Fig. 4 displays a comparison between the HOMO D dispersion of (+)-(1R,4R)-camphor (data digitalized and plotted from ref. 34) and (−)-(1R,4S)-norcamphor in the energy range of 5–20 eV; to make the comparison easier, the camphor data were multiplied by a factor 2 and shifted 2 eV towards lower kinetic energies to align the minima. The HOMO dispersions of the two molecules present the same sign and profile; the same asymmetric shape around the minimum is present with a longer tail towards higher kinetic energies. This result clearly proves from a phenomenological point of view that structural similarities are reflected on the PECD dispersion, and we believe that the same sign and the resemblance of the two minima are an indication of the similar molecular structures.
Fig. 4 Experimental D dispersions of the HOMO of (+)-(1R,4R)-camphor (red dots) and (−)-(1R,4S)-norcamphor (blue dots). The camphor HOMO dispersion has been multiplied by a factor 2 and shifted 2 eV toward lower kinetic energies to facilitate the visual comparison. Camphor data are digitalized and plotted from ref. 34. Error bars are omitted. The (1S,4R)-norcamphor data have been negated to obtain the (1R,4S)-norcamphor D. |
In the camphor HOMO study, both MS-Xa and B-spline theories reproduce in a satisfactory way the shape of the experimental dispersion in the intermediate energy range, in particular the oscillation around the minimum at 12.2 eV. The MS-Xa theory also follows the oscillations closer to the thresholds. Preliminary calculations using the B-spline code show that this behavior is strongly dependent on the choice of the local potential. At the quantitative level, the comparison of the sign and intensity of the theoretical and experimental dispersion fails in some region of the investigated range, in fact the sign of the experimental data is always negative while the theoretical values oscillates between positive and negative values. The different agreement of the theory of norcamphor and camphor, together with the close similarity of a region of the dispersion, directly indicate that the influence of the rotating methyl groups cannot be considered negligible in the theoretical modeling of the experimental data.
To strengthen this consideration we repeat here that a similar conclusion was drawn for the PECD study of the structural isomers alaninol and isopropanolamine,49 where a strong resemblance between alaninol and isopropanolamine was found for the experimental HOMO and HOMO−1 dispersion, in spite of a qualitative agreement with the DFT B-spline calculations. The calculations at the equilibrium geometry do not fully reproduce the experimental data; nevertheless, the experimental dispersion exhibits similar shapes.
Although at the present stage it is not possible to convert the line of reasoning exposed above into a quantitative result, like a data analysis procedure or a prescription in the calculation, we believe it grasps the important behavior of the PECD conformational sensitivity.
The findings of the C 1s PECD of norcamphor were also studied. Fig. 5 reports the 1s PES core levels acquired at the photon energy of 320 eV. The spectrum is composed of two separate peaks that are assigned to the C 1s signal of the carbonyl group (IE = 293.0 eV) and to the other alkylic C 1s contributions (IE = 290.6 eV). The PECD measurements were performed on the carbonyl-related peak to investigate a single atomic C 1s signal. Fig. 6 shows the comparison between the experimental data and the theoretical DFT and TDDFT calculations. The calculations were shifted 1.2 eV towards higher kinetic energies. The calculations are in good agreement with the descending part of the dispersion and also the shoulder. However, for kinetic energies greater than 20 eV the agreement is hampered by the poor signal-to-noise ratio, and the calculated dispersion is in satisfying agreement close to the experimental error bar.
Fig. 5 PES C 1s of norcamphor taken at a photon energy 305 eV. The ionization energy of the C atom of the carbonyl group is 293 eV. |
In the literature, both MS-Xa and B-spline LCAO also reproduce well the core levels in molecules where the valence band agreement is qualitative or not satisfactory. The prevalent atomic nature of the C 1s core level provides a sharply defined initial state, since for the valence band the theory assignment is more dependent on the chosen theoretical method. This simplification is certainly helpful in the accuracy of the theoretical calculation. Indeed, the conformational geometry effects on the PECD have a different origin in valence band and in the core level states. Despite the local character of the core level excitation, the sensitivity is neither fully local nor fully extended, as in the case of the valence band, and the balance between the local excitation and the whole molecular potential probed by the photoelectrons gives to the core level PECD a sufficient sensitivity to the molecule geometry. As reported above, in methyl oxirane C 1s calculations, the rotations of the methyl group exert their influence also for C atoms that are not directly connected with methyl groups.25
Moreover, it is important to note that the pure p-wave continuum generated by core 1s ionization limits the number of accessible angular momentum channels of the photoelectron continuum wavefunction, and PECD directly probes the rescattering generated by the neighboring atoms.
It is worth noticing that for camphor the experimental C 1s PECD dispersion of the carbonyl group is well reproduced using MS-Xa theory.50 To qualitatively understand the different agreement of the calculations for the valence band and C 1s core of camphor, it could be likely that the methyl groups are far enough to have a tiny influence on the photoelectron wavefunction in the proximity of the C absorbing atom and the effect of the rotations could be considered negligible, although it is not possible to rule out that the Boltzmann average of the configurations approximates the equilibrium calculation. As reported in Fig. 4 for the HOMO dispersion, Fig. 7 shows the experimental D dispersion for the C 1s of the carbonyl group for (+)-(1R,4R)-camphor and (−)-(1R,4S)-norcamphor. The camphor data are extracted and plotted from ref. 50, where the reported asymmetries were divided by cos(54.7°) to obtain the D parameter. As previously reported for the HOMO, there are also similar features in the C 1s dispersion, where the descending part and the negative region of the minimum have the same shape, although the minimum is not clearly outlined for norcamphor. The zero crossing of norcamphor stays at higher electron kinetic energies with respect to camphor.
Fig. 7 Experimental D dispersions of C 1s of (+)-(1R,4R)-camphor (red dots) and (−)-(1R,4S)-norcamphor (blue dots). The camphor C 1s dispersion has been multiplied by a factor 0.5 to facilitate the visual comparison. Camphor data are extracted and plotted from ref. 50, where the reported asymmetries were divided by cos(54.7°) to obtain the D parameter. Error bars are omitted. The (1S,4R)-norcamphor data have been negated to obtain the (1R,4S)-norcamphor D. |
From the comparison of experiment and theory in norcamphor and camphor, we outline the general conclusions that help us to understand some open issues on PECD. State-of-the-art DFT and TDDFT quantum theory at the equilibrium geometry can correctly describe the PECD of molecules in the absence of rotating groups. Rotamers can exert a strong influence on PECD and the average on the configurations can resolve a theoretical dispersion that is different from the dispersion calculated at the equilibrium geometry. Another important result is that the comparison of the experimental HOMO dispersions for camphor and norcamphor shows that there are similarities, due to the identical skeleton structure. When equilibrium geometry calculations are employed in molecules where conformers and/or rotamers are present the outcome can differ from the experimental dispersion although in a not predictable way. However, it is also evident from previous studies that equilibrium geometry calculations always describe the essential features, especially in determining the maxima and minima of the D dispersion as a function of the kinetic energy, and can be posed as a solid foundation of PECD spectroscopy. Due to the possible slight energy misalignment between the experimental and calculated PECD profiles, it is important that the experimental data cover a substantial energy region, i.e., at least several eV, beyond the threshold region.
In this context, future efforts will be devoted to incorporating a treatment to take into account the distribution of the conformational geometries, providing an efficient choice of the ensemble of configurations to avoid the unsustainable burden of brute-force calculations.
A fair trade off between the calculation accuracy and the number of configurations is desirable to avoid to wasting calculation time that can be caused by attempting in a not-needed accuracy or oversampling the configurations. A subsequent question that cannot be ducked is the influence of the density functional on the PECD calculations; in this and previous studies performed by the authors, the density functional LB94 was the preferred choice due to the correct tail for continuum calculations, but it could be of paramount importance to explore the correlation with other functionals. PECD intertwines the structure and electronic properties, which is a gift and a curse because it is not easy to isolate and recognize which approximation is not corresponding to the actual description of the molecular system. From this perspective, the effort of the past two decades has given a solid foundation, although loose ends still remain if the power of PECD is to be fully harnessed.
Given the importance of theoretical simulations for extracting maximum information from experimental data, the development of predictive computational tools is mandatory. The typical size of chiral molecules, and possible conformational averaging, makes higher-level treatments hardly affordable, and spurs the optimization of local and possibly non-local single-particle potentials, especially for the lower kinetic energy region where important PECD signatures are found. To this end, a set of accurate experimental data on rigid systems, to avoid ambiguities in conformational distribution, is highly desirable to serve as a benchmark for theoretical development.
Although circularly polarized core-level spectroscopy is traditionally rooted in X-ray magnetic circular dichroism, the increasing interest in PECD fosters the demand for chemical physics investigations of chiral molecules.
Exploiting core-level PECD, profiting from the reduced number of angular momentum channels and the local nature of the excitation could be of practical importance to access both the electronic and structural properties of chiral molecules in the gas phase.
This journal is © the Owner Societies 2022 |