# Ultrafast probing of isotope-induced explicit symmetry breaking in ethylene

Alessandro Nicola Nardi, Alexie Boyer, Yaowei Hu, Vincent Loriot, Franck

Lépine, Morgane Vacher, Saikat Nandi

► To cite this version:

Alessandro Nicola Nardi, Alexie Boyer, Yaowei Hu, Vincent Loriot, Franck Lépine, et al. Ultrafast probing of isotope-induced explicit symmetry breaking in ethylene. Communications Chemistry, 2025, 8, pp.222. 10.1038/s42004-025-01621-z. http://dx.doi.org/10.1038/s42004-025-01621-z

HAL Id: hal-05263331

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

Submitted on 16 Sep 2025

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_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Ultrafast probing of isotope-induced explicit symmetry breaking in ethylene

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

Check for updates

Alessandro Nicola Nardi Ⓓ1,3, Alexie Boyer Ⓓ2,3, Yaowei Hu², Vincent Loriot Ⓓ2, Franck Lépine Ⓓ2,

Morgane Vacher Ⓓ1✉ & Saikat Nandi Ⓓ2✉

Symmetry governs nature's laws, yet many of the natural phenomena occur due to the breakdown of symmetry. Here, we show how isotope-induced inversion symmetry breaking influences ultrafast photoisomerization processes in ethylene. Using extreme ultraviolet pump-near infrared probe time-of-flight mass spectrometry, we find that replacing one of the carbon atoms in ethylene with a ¹³C isotope leads to twice-faster structural relaxation via ethylene-ethylidene isomerization in the photo-excited molecular cation. Advanced trajectory surface hopping calculations incorporating the nuclear symmetry of the molecular systems, reveal that it arises from the mixing of different normal modes in the isotope-substituted species, interactions otherwise forbidden by symmetry. Although the mixing does not alter the symmetry of the electronic Hamiltonian, it modifies that of the nuclear Hamiltonian, causing explicit symmetry breaking. This facilitates efficient intra-molecular vibrational energy redistribution, lowering the isomerization yield. Our findings offer opportunities to use isotope-induced nuclear symmetry breaking to control the outcome of light-molecule interactions across ultrafast timescales.

Since Pierre Curie's seminal work<fnref n="1" />, symmetry and symmetry breaking have been known to play a central role in contemporary physics. Notable instances include gauge-symmetry breaking in superconductors², or, in the Higgs mechanism³. In these situations, the breakdown of symmetry takes place in systems that have infinite degrees of freedom, and thus an asymmetric ground state⁴. For small molecules, such as diatomic, or triatomic systems, the ground state is highly symmetric due to the reduced number of electronic and nuclear degrees of freedom. In a pioneering experiment with H₂ molecules, it was shown that such a symmetry can be broken by absorption of a linearly polarized photon⁵. With the development of tabletop high-order harmonic generation (HHG) based sources as well as free-electron lasers (FELs) providing attosecond (1 as = 10⁻¹⁸ s) and femtosecond (1 fs = 10⁻¹⁵ s) pulses in the short-wavelength region (few tens of nanometers to few nanometers), there is a growing interest in the ultrafast breakdown of symmetry in polyatomic molecules with high-symmetry point groups⁶,⁷. For instance, femtosecond transient X-ray absorption has been used to study symmetry breaking in tetrahedral molecules, such as methane⁸ and tetrafluoromethane⁹. The signature of symmetry breaking in Jahn-Teller-distorted methane cation was also visualized using ultrafast Coulomb explosion imaging¹⁰. Attosecond charge migration in hydrogen and water molecules has been predicted to break electronic symmetry while preserving the nuclear one¹¹. Carbon dioxide dimers photoionized by

extreme ultraviolet (XUV) pulses from an FEL are shown to undergo femtosecond symmetry-breaking dynamics leading to charge asymmetry¹². These studies collectively demonstrate how ultrashort radiation can induce symmetry-breaking dynamics in various molecular systems.

Here, we investigate a different approach to symmetry breaking: substituting one of the carbon atoms in a neutral ethylene molecule with a ¹³C isotope. This substitution reduces the molecular symmetry from D₂h to C₂v, allowing us to explore how such a change influences light-induced relaxation dynamics across ultrafast timescales. The effect of isotope-induced inversion symmetry breaking is well documented in the frequency domain, especially for hydrogen¹³⁻¹⁵ and nitrogen¹⁶ molecules. Isotopic substitution changes the zero-point energy of the molecule resulting in kinetic isotope effects, which were recently investigated across few-femtosecond¹⁷ and attosecond¹⁸,¹⁹ timescales. The present study examines how the explicit breakdown of the inversion symmetry that preserves the symmetry of the electronic Hamiltonian, but not of the nuclear Hamiltonian, affects photoisomerization in ethylene cation. The electronic Hamiltonian depends only on the positions and charges of the nuclei, making it insensitive to mass differences. However, the nuclear Hamiltonian includes a kinetic energy term that depends on nuclear masses, so substituting one carbon with a heavier isotope explicitly breaks the inversion symmetry in the nuclear motion. Following photoionization by a broadband XUV attosecond pulse

<footnote n="1">¹Nantes Université, CNRS, CEISAM, UMR 6230, F-44000 Nantes, France. ²Universite Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, UMR 5306, F-69100 Villeurbanne, France. ³These authors contributed equally: Alessandro Nicola Nardi, Alexie Boyer. ✉e-mail: morgane.vacher@univ-nantes.fr; saikat.nandi@univ-lyon1.fr</footnote>
train produced via HHG, the resulting C₂H₄⁺ was promoted to an electronic excited state from which it could relax via different pathways¹⁷,²⁰: H-loss, H₂-loss, ethylene-ethylidene (CH₃CH⁺) isomerization. A time-delayed femtosecond near-infrared (NIR) pulse probed the relaxation dynamics in the cation by further ionizing and/or, dissociating it. The CH₃⁺ fragment carried the signature of the relaxation pathway via ethylene-ethylidene conversion²⁰⁻²². We interpreted the corresponding ion yields, measured as a function of the pump-probe delay, using non-adiabatic trajectory surface hopping simulations. These simulations revealed that the difference in relaxation timescales between the isotopologues arises from the mixing of different normal modes, otherwise symmetry-forbidden, enabled by the explicit symmetry breaking. Furthermore, our calculations showed that the isotopic substitution, accompanied by a reduction in molecular symmetry, significantly suppresses the isomerization yield.

## Results

Absorption of a photon from the attosecond pulse train (energies between 17.0 and 23.3 eV) can produce ethylene cation mainly in four different cationic states from D₀ to D₃, with only a modest contribution from the D₄ state (see Supplementary Note 1 for details about the energetics of the photoionization process from our ab initio calculations and previous literature). Typical two-color ion-yields as a function of the XUV pump - NIR probe delay for ¹²CH₃⁺ from ¹²C₂H₄, ¹³CH₃⁺ from ¹³C₂H₄ and ¹³CH₃⁺ from ¹³C¹²CH₄, are shown in Fig. 1a-c. The appearance energy for this fragment was measured to be around 20 eV, with the probability of production reaching a maximum at around 27 eV²³. The interpretation of exponential decays observed in ion-yields in XUV pump - NIR probe experiments has been addressed in multiple studies over the years²⁴,²⁵. Here, to understand the time-dependent signal, we make the distinction between the time it takes to form the ethylidene isomer and the timescale of the dynamics via that relaxation channel¹⁷. The molecule in an excited cationic state relaxes non-adiabatically via certain pathways and the NIR depopulates the state by promoting a part of the population to a higher-lying dissociative state leading to the production of CH₃⁺ ions. If the probe arrives at a delay by which time the molecule has already relaxed electronically and/or, structurally via change in nuclear coordinates, the NIR pulse may no longer promote the system to the dissociative state to generate any CH₃⁺ ion. This causes an exponentially decaying signal as a function of the delay between the two pulses. The same mechanism has been put forward previously for describing pump-probe delay-dependent fragmentation yields in ethylene

and deuterated ethylene¹⁷. Recently, Vismarra et al.²⁶ have developed a similar explanation behind observed time-dependent yields in XUV-ionized donor-π-acceptor systems. In their case, as the nuclear wave-packet evolves, the energy separation between the populated states and the dissociative final state increases, making multi-photon absorption less efficient over time. It progressively reduces the likelihood of promoting the wave-packet to the dissociative state of interest, ultimately leading to the observed decay in ion-yield. This is also conceptually consistent with the way the time-dependent signal is interpreted in the present manuscript. The corresponding experimental time constants (τ) from three different targets (see Methods for details about the data analysis) are shown in Fig. 1(d). The average decay constant for ¹³CH₃⁺ from ¹³C¹²CH₄ is 50 ± 8 fs, almost twice smaller than that of 108 ± 4 fs for ¹²CH₃⁺ from ¹²C₂H₄. For ¹³CH₃⁺ from ¹³C₂H₄, it is 83 ± 4 fs. They do not change significantly with the NIR probe intensity, confirming that the observed time-dependent ion yield is representative of the dynamics initiated by the XUV-pump only. It is interesting to note that the observed trend in the τ-values for different targets cannot be explained by mass effects²⁷. To better understand the trend observed in the measurements, and to identify the relevant electronic and/or nuclear relaxation pathway(s) associated with the two-color signal, we simulated the non-adiabatic relaxation dynamics using trajectory surface hopping calculations within the OpenMolcas code²⁸,²⁹, including the five lowest energy cationic states in all three isotopologues (see Methods for details). The active space was chosen to be large and flexible enough to describe bond breaking and formation. The initial conditions were generated using Newton-X³⁰, with phase space sampling performed in an uncorrelated manner based on the Wigner distribution³¹. While both ¹²C₂H₄ and ¹³C₂H₄ belong to the D₂h symmetry point group, a single ¹³C-isotopic substitution leads to symmetry breaking and the resulting ¹³C¹²CH₄ molecule belongs to the C₂v point group. Due to the stochastic nature of the sampling process, the different sets of initial conditions do not strictly fulfill symmetry. To correct for this, the initial conditions of each isotopologue were symmetrized a posteriori according to its point group. In particular, the geometries explored along the non-adiabatic simulations were replicated by applying the symmetry elements (except for the identity) of the corresponding point group: seven elements (3 symmetry axes, 3 mirror planes, and the inversion) for ¹²C₂H₄⁺ and ¹³C₂H₄⁺, and three elements (the symmetry axis and 2 mirror planes) for ¹³C¹²CH₄⁺, at the equilibrium geometry (see Methods for details).

The Dyson amplitudes (an approximation of photoelectron intensities) were calculated for 100 of the initial conditions for each isotopologue

<figure><img src="image_4.png" /><figcaption>**Fig. 1 |** Experimental results. Typical two-color signal as a function of the XUV pump -- NIR probe delay for a ¹²CH₃⁺ from ¹²C₂H₄, b ¹³CH₃⁺ from ¹³C₂H₄ and c ¹³CH₃⁺ from ¹³C¹²CH₄. The error-bars correspond to statistical fluctuations. In each panel, the solid line denotes fitting with an exponentially decaying function (see Methods for details). d The corresponding time constants, τ, at different NIR-probe intensities for three different isotopologues: C₂H₄ (solid red circles), ¹³C₂H₄ (solid green diamonds) and ¹³C¹²CH₄ (solid blue squares). The uncertainties stem only from the fitting. In each case, the dashed line corresponds to the weighted average and the shaded area represents the 95% confidence interval.</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>**Fig. 2** | Ethylene-ethylidene photoisomerization yield. Theoretical ethylidene population as a function of time for the parent ions *a* C₂H₄⁺, *b* ¹³C₂H₄⁺, and *c* ¹³C¹²CH₄⁺, upon excitation to D₃ electronic state. The dashed-dotted lines correspond to the fitting with a logistic function (see main text for details).</figcaption></figure>

from the neutral ground state to the first five electronic states of the cation. The results (see Supplementary Fig. 1) show the independence of the estimated photoelectron intensities on the isotopic substitution and the low photoionization intensity for the D₄ state. Only about 15% of the initial conditions for each isotopologue showed a Dyson norm greater than 0.25 for the latter cationic state, suggesting that this state contributes only moderately to the experimental signal. We calculated the time-dependent yields for all possible relaxation channels, namely, H-loss, H₂ and/or, 2H-loss, and ethylene-ethyldiene isomerization. In general, for relaxation from the excited cationic states, D₁, D₂, D₃, and D₄, the corresponding trends observed for three different isotopologues approximately follow the change in mass (see Supplementary Fig. 2). However, for ethylene-ethylidene conversion, the calculations show that the ethylidene population saturates faster and to a lower yield in the symmetry-broken species compared to the other two for relaxation from the D₃ state (see Fig. 2a-c, solid lines). Note that the ethylidene population curves increase in a non-monotonic manner; this is because the isomerization reaction is reversible. To quantitatively compare the calculations with the experiment, the theoretical ethylidene populations were fitted with a logistic function, $f(t) = \frac{P}{1+e^{-k(t-t_{1/2})}}$, where P

denotes the asymptotic population, *k* is the population growth rate, and *t*₁/₂ signifies the time to form half of the asymptotic ethylidene population. The fits reported in Fig. 2 (dashed-dotted lines) were performed assuming *P*, *k* and *t*₁/₂ as free parameters. To verify the robustness of the fitting procedure, we checked the fitting with normalized *P* as well as fixed *k* for three isotopologues. While the *R*² value reduces by less than 5% for constrained fitting, *t*₁/₂ do not change significantly (see Supplementary Fig. 3a, b and Supplementary Table 1). Although ¹²CH₃⁺ fragments from ¹³C¹²CH₄⁺ could not be detected due to overlapping signal from ¹³CH₂⁺, we have calculated the ethylidene population as a function of time for the ¹³C¹²CH₄⁺ → ¹²CH₃⁻·¹³CH⁺ isomerization channel. The yield for this channel saturates at a similar rate compared to that for the ¹³C¹²CH₄⁺ → ¹³CH₃⁻·¹²CH⁺ channel (see Supplementary Fig. 4), indicating both relaxation channels to be sensitive towards the effect of isotope-induced symmetry breaking.

The fits reflect a significantly lower asymptotic ethylidene isomer population in ¹³C¹²CH₄⁺: *P* = 9%, compared to 14% in ¹²C₂H₄, and 13% in ¹³C₂H₄. The population saturates faster: *k* = 0.082 fs⁻¹ in ¹³C¹²CH₄, in contrast with 0.032 fs⁻¹ in ¹²C₂H₄. For ¹³C₂H₄, an intermediate value of 0.063 fs⁻¹ is obtained. More importantly, the inflection point is reached earlier in time in the symmetry-broken molecule. Indeed, *t*₁/₂ is 62 fs in

¹³C¹²CH₄, 103 fs in ¹²C₂H₄, and 80 fs in ¹³C₂H₄, which matches well experimental τ-values (see, Table 1). This also implies that if the NIR-probe arrives at a delay by which the relevant excited state population in D₃ has relaxed via the formation of around 50% of the asymptotic ethylidene population, it can no longer efficiently produce the CH₃⁺ ions from the ethylene cation, leading to the exponential decay observed in the experimental two-color signal. Our finding that CH₃⁺ yield primarily reflects relaxation from a higher-lying excited cationic state, such as D₃, is supported by previous observation³². As in the measurements, while the differences in the fitted parameters between iso-symmetric ¹²C₂H₄ and ¹³C₂H₄ may be explained in terms of changes in mass, the predicted behavior of the symmetry-broken ¹³C¹²CH₄ cannot be. Nevertheless, the trend observed in the experimental decay constants for three isotopologues agree with that observed in the theoretical timescales for the ethylene-ethylidene isomerization channel. This allows us to assign the time-dependent ¹²CH₃⁺ (or, ¹³CH₃⁺) yield with the structural relaxation pathway via ethylene-ethylidene isomerization.

## Discussion

For ethylene and ethylene-¹³C₂, the representation of the twelve normal modes in D₂h symmetry point group is 3Ag + Au + 2B₁u + B₂g + 2B₂u + 2B₃g + B₃u. As the inversion symmetry is broken, the representation in the resulting C₂v point group becomes 5A₁ + A₂ + 2B₁ + 4B₂. The correlation relations are as follows: Ag, B₁u → A₁; Au → A₂; B₂g, B₃u → B₁; B₂u, B₃g → B₂. The normal modes within symmetries A₁, B₁ and B₂ in ¹³C¹²CH₄ are allowed, from a symmetry point of view, to couple and exchange vibrational energies with more modes compared to those having D₂h point group, thus affecting the subsequent relaxation dynamics in the excited cation.

To investigate the effect of the breakdown of the inversion symmetry on the coupling of normal modes, we conducted a principal component analysis (PCA) on the geometries visited along the trajectories initiated upon excitation to the D₃ electronic state for each isotopologue (see Supplementary Note 2 for details about the PCA). The frames where H-loss and H₂ (or, 2H)-loss occurred were excluded from the analysis to focus only on the ethylidene conversion channel. The resulting principal components indicate the essential motions for this reaction pathway. For each isotopologue, the components were projected onto the corresponding normal modes (see Supplementary Fig. 5 for pictorial depiction of the normal modes in ethylene) and analyzed to identify which modes contribute to the ethylene-ethylidene isomerization.
The outcomes of the analysis for D₃ electronic state are presented in Fig. 3a, b for ¹²C₂H₄ and ¹³C₂H₄ (results for other initial electronic states are shown in Supplementary Fig. 6). These maps show the absolute values of the coefficients of each principal component (columns) in terms of normal modes (rows). The principal component (PC) analysis shows that the essential motions for the ethylidene isomerization in ¹²C₂H₄ are along CH₂ twisting normal mode (PC0), CH₂ scissoring, C = C stretching and CH₂ stretching (PC1 and PC3), CH₂ rocking (PC2) and CH₂ scissoring and CH₂ stretching (PC4). No significant differences were observed between the two iso-symmetric isotopologues, except for permutation of some principal components with similar associated variance. However, the analysis in symmetry-broken ¹³C¹²CH₄ (see Fig. 3c) displays significant mixing between normal modes (now symmetry-allowed) with the appearance of three different clusters. For example, two CH₂-scissoring modes with A₉ and B₁₀ symmetries in D₂h that transforms into A₁ symmetry in C₂v now mix with each other. Looking at the principal components 1, 2 and 4 (circled region in Fig. 3c), which individually explains a variance ratio of about 0.4, 0.1, and 0.1, respectively, we found that in ¹³C¹²CH₄ they contain significant contributions from five different normal modes, all of them with the A₁

<table><caption>**Table 1 | Comparison between experimental and theoretical timescales**</caption><thead><tr><th>Parent</th><th>Fragment</th><th>τ (fs)</th><th>t<sub>1/2</sub> (fs)</th><th>P (%)</th><th>k (fs<sup>-1</sup>)</th></tr></thead><tbody><tr><td>¹²C₂H₄</td><td>¹²CH<sub>3</sub><sup>+</sup></td><td>108 ± 4</td><td>103 ± 1</td><td>13.5 ± 0.1</td><td>0.032 ± 0.001</td></tr><tr><td>¹³C₂H₄</td><td>¹³CH<sub>3</sub><sup>+</sup></td><td>83 ± 4</td><td>80 ± 1</td><td>12.9 ± 0.1</td><td>0.063 ± 0.002</td></tr><tr><td>¹³C¹²CH₄</td><td>¹³CH<sub>3</sub><sup>+</sup></td><td>50 ± 8</td><td>62 ± 1</td><td>8.8 ± 0.1</td><td>0.082 ± 0.003</td></tr></tbody></table>

The experimental and theoretical timescales, along with the asymptotic population and the population growth rate. The uncertainties originate from the fitting only.

symmetry. The PCA was performed also on the dynamics upon excitation to the D₁, D₂, and D₄ state, for the three isotopologues, and the results are reported in Supplementary Fig. 6 and the explained variance ratio for all principal components are shown in Supplementary Fig. 7.

For the isotope-substituted molecule, the symmetry-breaking kinetic energy term in the Hamiltonian of the molecule primarily depends on the mass difference as $\propto \Delta m(\nabla_R \cdot \sum_j \nabla_j)$, where R is the position of the substituted nuclei and j denotes the relevant electronic coordinates¹³. For small value of the mass difference $\Delta m$, such as in the present case, it approximately preserves the overall Hamiltonian symmetry. In this situation, the pure electronic relaxation still respects the approximate symmetry; this is seen in the time-dependent electronic populations for different excited states, which do not change for the three isotopologues (Supplementary Fig. 8). However, the structural relaxation mediated by several nuclear coordinates does not respect the approximate symmetry. Hence, the isotope-induced inversion symmetry breaking is an example of explicit symmetry breaking⁴, as opposed to spontaneous symmetry breaking typically associated with Jahn-Teller effects⁸,¹⁰. The XUV-induced ultrafast ethylene-ethylidene conversion process serves as an ideal proving ground for this mechanism. Using the PCA, we clearly identified normal-mode mixing as a key consequence of this explicit symmetry breaking (Fig. 3). The mixing causes an efficient intra-molecular vibrational energy redistribution³³, decreasing the overall isomerization yield and thus, the asymptotic population of ethylidene, P. It is worth noting that the differences between the ethylidene yields of the different isotopologues upon excitation to the D₄ state (see Supplementary Fig. 2) can also be rationalized considering the intramolecular vibrational energy redistribution that follows from the D₂h to C₂v point-group symmetry-lowering induced by a single ¹³C-atom substitution, suggested by the normal mode mixing in the PCs reported in Supplementary Fig. 6g-i.

To get a qualitative picture behind how the mixing between different normal modes can influence the ethylene to ethylidene conversion as

<figure><img src="image_6.png" /><figcaption>**Fig. 3 | Principal component analysis.** The absolute values of the principal component coefficients with D₃ initial electronic state for **a** C₂H₄, **b** ¹³C₂H₄, and **c**, ¹³C¹²CH₄. Note the mixing between different normal modes (denoted by</figcaption></figure>

dashed circle) in ethylene-¹³C following the loss of the g/u-symmetry. The explained variance ratio for different principal components are shown in Supplementary Fig. 7.
**Fig. 4** | **Temporal variation of nuclear coordinates.**
*a* The temporal evolution of the dihedral angle, |*Φ*|, defined by H³C¹C²H⁶ (solid lines) and H³C¹C²H⁵ (dashed lines) in different isotopologues (¹²C₂H₄ in red, ¹³C₂H₄ in green, and ¹³C¹²CH₄ in blue) upon excitation to the D₃ electronic state. *b* The same as *a* but for the distance, *d*, between carbon atom C¹ and hydrogen atom H⁵. *c* Once again the same as *a* but for the angle, α, identified by atoms H⁵C²C¹. The inset in each case shows the corresponding definition for the coordinate. *d* The time dependent variation of the average normal mode (CH₂ scissoring) displacements. In all the analyses the frames in which H-loss and H₂ or, 2H-loss takes place are excluded.

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

reported in Fig. 2, we studied the temporal evolution of several internal (nuclear) coordinates of the photoionized molecule, upon excitation to the D₃ electronic state. Like in the PCA, the frames where H-loss and H₂ (or, 2H)-loss occurred were excluded. The three relevant internal coordinates are the absolute dihedral angle, |Φ|, the average distance *d*, between carbon atom C¹ and hydrogen atom H⁵, and the angle, α identified by atoms H⁵C¹C²H⁵ (see Fig. 4 insets). For all the panels in Fig. 4, the temporal variations were smoothed using a Gaussian kernel with a full width at half maximum of 36 fs (same as the XUV-NIR cross-correlation) to better reveal the underlying trend. Ethylene-ethylidene isomerization requires the combination of a rotation around |Φ|, a closure of the angle α, and a shortening of the bond length *d*. Fig. 4a shows the temporal variation of |Φ|, identified by the atoms H⁵C¹C²H⁵ and H⁵C¹C²H⁶: the differences between the three different isotopologues are not significant. This is expected since motion around |Φ| corresponds to a normal mode of A_u symmetry in D₂h, A₂ in C₂v, and is not allowed to couple to more modes upon breakdown of the inversion symmetry. However, for *d* and α internal coordinates, clear differences between the different isotopologues can be observed (see, Fig. 4b, c). At longer times (>100 fs), the values of the coordinates *d* and α decrease further in the isosymmetric systems, in particular for ¹³C₂H₄. To understand these evolutions, it is important to realize that the time evolutions of average nuclear coordinates are representative of both the yields (*P*) and the rates (*k*) of the processes. The larger distortion of *d* and α in the ¹³C species results in a relatively large asymptotic population, *P*, and an intermediate rate, *k*. Interestingly, ¹²C₂H₄ and ¹³C¹²CH₄ show a comparable change in the *d* and α coordinates within the simulation time. This may appear strange at first glance but can actually be explained by two effects compensating each other. On one hand, the ¹²C₂H₄ is the isotopologue with the lowest rate, *k*, but the highest yield. On the other hand, the symmetry-broken isotopologue shows the highest rate but lowest yield (see Table 1). These two aspects of the dynamics compensate each other and lead to equivalent variations of these coordinates within the simulated time window. Now, comparing the two iso-symmetric species, they have almost the same asymptotic yield but ¹²C₂H₄ isotope has a slower rate, hence the smaller distortion as a function of time. A similar trend was also found for the average of the relevant normal mode (CH₂ scissoring) displacements (see Fig. 4d). We performed the same

analysis for all the other normal modes by filtering everything but ethylene to ethylidene conversion as well as trajectories ending in ethylidene; for both filters the effect of symmetry breaking was only observed for the relevant modes with A_g symmetry, further reinforcing our interpretation (see Supplementary Figs. 9 and 10). After 100 fs, the coordinates change more slowly in the symmetry-broken system, suggesting that much of the excited-state population has already relaxed structurally by this time. For the two isosymmetric molecules, especially ¹³C₂H₄ species, even at longer times, there is still a significant excited-state population left to decay, causing noticeable changes in the coordinates. As pointed out earlier, this is reflected in the experimentally measured faster decay timescale of the CH₃⁺ ion yield (τ) with respect to the symmetry-intact molecules.

We have shown that isotope-induced explicit symmetry breaking can significantly alter structural relaxation timescales, despite minor changes in the overall mass, due to mixing of normal modes in the photo-excited molecular cation. While the decay of electronic populations remains unaffected by the mixing, it can significantly suppress photo-isomerization yields via intra-molecular vibrational energy redistribution. Isotope substitution has previously been used to change the outcome of light-induced electron transfer processes in molecules³⁴. Our results highlight the potential of symmetry breaking to control similar phenomena in molecules. With the advent of brighter light sources having shorter wavelengths and ultrashort pulse duration, the present study can act as a benchmark not only for future experiments on isotope-induced symmetry-breaking in more complex systems but also to design new photoresponsive molecules based on specific symmetries.

## Methods

### Experimental setup

The experiment was performed at the Circé attosecond beamline at the Institut Lumière Matière, Lyon³⁵. The output of a commercially available Ti:Saph amplifier laser system (Coherent Legend), providing NIR pulses (wavelength: centered around 800 nm) at a repetition rate of 5 kHz with 2 mJ energy per pulse, was divided into two parts along the arms of a Mach-Zehnder interferometer. In the pump arm, the XUV attosecond pulse train was generated by focusing the NIR beam onto a jet of noble gas atoms (either
xenon or, krypton). The attosecond pulse trains produced via HHG represent a comb of odd-order high harmonics in the frequency domain, out of which only three, namely 11, 13 and 15 were selected using a tin filter (see Supplementary Fig. 11 for the HHG spectra). The Sn-filter also removes the co-propagating NIR-pulse. Using the chirp-scan technique³⁶, we determined the temporal duration of the NIR-pulse to be 27.5 ± 0.2 fs (see Supplementary Fig. 12 for the NIR temporal profile). In the probe arm, the remaining NIR beam is focused with a 1-m lens onto the interaction region of a velocity map imaging (VMI) spectrometer. The ionic fragments were collected across 4π solid angle by the VMI spectrometer operating in the time-of-flight mode. We procured the ¹³C-isotope substituted samples from Sigma-Aldrich (CAS 6145-18-2 and 51915-19-6) and used them without further purification. The molecules were introduced in the interaction region as an effusive jet via a 200 µm hole in the repeller plate of the VMI spectrometer. The values of the probe intensity reported in Fig. 1d were estimated to have an uncertainty up to 20%. The cross-correlation between the XUV pump and NIR probe was 35 ± 1 fs (see Supplementary Fig. 13 for cross-correlation measurements in argon).

## Data analysis

The data reported in Fig. 1a–c consists of an average over 10 consecutive sets of forward (–160 fs to + 400 fs) and backward scans (+400 fs to –160 fs), thus, in total 20 scans across the entire range of the pump-probe delay, at each probe intensity. The error-bars corresponds to the statistical fluctuations for these 10 sets of scan, obtained as the standard error. Pure two-color signal was obtained by subtracting the XUV-only and NIR-only yield from the time dependent signal at each pump-probe delay. The XUV-only time-of-flight mass spectra normalized on the parent-ion-yield for all three isotopologues are shown in Supplementary Fig. 14. The two-color signal as a function of the pump-probe delay was fitted with an exponential decay function convoluted (*) with a normalized Gaussian denoting the instrument response function, $g_{\text{IRF}}(t) = \frac{1}{\sigma_p\sqrt{2\pi}}\exp[-(\frac{t}{\sigma_p\sqrt{2}})^2]$. The fitting function can, therefore, be written as, $g(t) = [(A \exp(-\frac{t-t_0}{\tau}) + B)\Theta(t-t_0)] * g_{\text{IRF}}(t) + C$. Here, $\tau$ is the time-constant, 2.355$\sigma_p$ provides a measure of the cross-correlation, the temporal overlap between the pump and the probe is given by $t_0$, and $\Theta(t)$ denotes a Heaviside step function. See Supplementary Note 3 for details about the precise determination of $t_0$ and $\sigma_p$. The linear background, C, was subtracted from the two-color signal for each measurement. The error-bar in $\tau$ for an individual measurement shown in Fig. 1d represents uncertainties obtained from the fitting procedure using the least squares method. The dashed line is the average, weighted against these uncertainties. The error on the weighted average was determined by the corresponding standard deviations of the four independent measurements.

## Theoretical methods

We performed non-adiabatic dynamics simulations of all the ¹³C-containing isotopologues of ethylene cation examined in this study (¹²C₂H⁺₄, ¹³C₂H⁺₄, and ¹³C¹²CH⁺₄). The simulations were conducted using the surface hopping method, with the Tully fewest switches algorithm³⁷, within the OpenMolcas²⁸,²⁹ software (version 22.10-354-g7f2c128³⁸). The trajectories were integrated with a nuclear time step of 20 a.u. for ca. 203 fs. Each nuclear time step was divided into 96 substeps for the electronic wave function integration, carried out using the Hammes-Schiffer Tully³⁹ scheme with biorthonormalization, as implemented in OpenMolcas⁴⁰. To account for decoherence effects, the Granucci and Persico energy-based correction⁴¹ was applied with a decay factor of 0.1 Hartree. After a hop, the entire classical velocity vector is rescaled to ensure energy conservation along the trajectory, and frustrated hops are ignored. The electronic structure of all the studied isotopologues was treated at the Complete Active Space SCF⁴² (CASSCF) level, with state-averaging performed over the five lowest-energy cationic states. See Supplementary Table 2 for the theoretical (vertical) ionization energies of ethylene from its neutral ground state to the first five cationic

states. The selection of the active space, consisting of 11 electrons in 12 orbitals (the $\sigma$ and $\sigma^*$ orbitals of the four C–H bonds, and the $\sigma$, $\sigma^*$, $\pi$, and $\pi^*$ orbitals of the C = C double bond), was based on previous works by some of the present authors¹⁷,²⁰. Indeed, the employed active space was demonstrated to be flexible enough to describe bond breaking and formation. While dynamic correlation may affect the XUV-induced dynamics⁴³, previous works investigating the influence of the dynamical correlation introduced at XMS-CASPT2 level²⁰ showed only a modest reduction of the D₀–D₁ energy gap at Franck-Condon geometry. Moreover, a semi-quantitative agreement between CASSCF and CASPT2 potential energy scans were found along the different relaxation pathways including the ethylene-ethylidene isomerization. The calculations employed the ANO-RCC-VDZP basis set⁴⁴ and Cholesky decomposition¹⁵ to reduce computational cost.

The frequencies used for this sampling were calculated at the ground-state equilibrium geometries of the respective neutral isotopologues. For the ¹³C₂H⁺₄ and ¹³C¹²CH⁺₄ isotopologues, 300 pairs of coordinates and velocities were generated and propagated starting from the D₁ and D₃ electronic states, 200 trajectories were initiated on D₀ and D₂. For the ¹²C₂H⁺₄ isotopologue, the ensemble of trajectories obtained in a previous work²⁰, which followed the same computational protocol, was used and analyzed. The set includes 300 trajectories started on the D₁, D₂, and D₃ electronic states, and 200 initiated on the D₀ state. Finally, 120 additional trajectories were initiated on the D₄ electronic state, even if this state showed very low photo-ionization probability (estimated as Dyson intensities on 100 of the initial conditions for each isotopologue) compared to the other three states (see Supplementary Fig. 1). A small number of trajectories, specifically 0.0%, 0.3%, 3.3%, 7.3%, and 10.9 % of those initiated on D₀, D₁, D₂, D₃, and D₄, respectively, for the ¹²C₂H⁺₄ isotopologue, exhibited a change in total energy exceeding 0.5 eV during the dynamics. Similar percentages were observed for the other two isotopologues, and these trajectories were discarded in the analysis process.

To gain insights into the reactive events in the collected ensemble of the trajectories, accordingly to a previous work²⁰, the ones with exactly three hydrogen atoms within 2.9 Å (value corresponding to the sum of the van der Waals radii of H and C atoms) from their nearest carbon atom were considered as having undergone H-loss. Similarly, trajectories with exactly two hydrogen atoms within 2.9 Å of their nearest carbon atom were considered to undergo through H₂-loss or 2H-loss. Finally, trajectories where all four hydrogen atoms are less than 2.9 Å from their nearest carbon atom, with three hydrogen atoms being closer to one carbon than the other, were classified as the ethylidene form. To account for the symmetries in the isotopologues, the geometries explored along the non-adiabatic simulations were replicated by applying the symmetry elements of the corresponding point group at the equilibrium geometry. For each explored geometry of ¹²C₂H⁺₄ and ¹³C₂H⁺₄, seven additional geometries were generated by performing the symmetry operations corresponding to the three symmetry axes, three mirror planes, and the inversion. In case of ¹³C¹²CH⁺₄, only three additional geometries were produced by performing the symmetry operations corresponding to the symmetry axis and two mirror planes. This enabled us to account for the nuclear symmetry of the investigated isotopologues. The statistical analyses were performed keeping the ensembles of trajectories that started upon excitation on D₀, D₁, D₂, D₃, and D₄ separated for each isotopologue. In the case of the principal component analysis, conducted on each ensemble separately, the frames of the corresponding trajectories were joined, discarding the ones in which H-loss or 2H and/or H₂-loss took place.

To estimate the photoelectron intensities, we calculated the Dyson amplitudes for 100 of the initial conditions for each isotopologue from the neutral ground state to the first five electronic states of the cation at the same level of theory used for the non-adiabatic dynamics simulations, i.e., SA5-CASSCF(11e,12o). The geometries obtained from the symmetrization of the initial conditions exhibit the same Dyson norm of the corresponding Wigner-sampled structure. Finally, initial conditions were not selected on the basis of its Dyson norm.
**Data availability**

Any additional data related to the findings are available from the corresponding authors upon request.

**Code availability**

Codes developed for this study are available from the corresponding authors upon request.

Received: 21 January 2025; Accepted: 21 July 2025;
Published online: 31 July 2025

**References**

1. Currie, P. Sur la symétrie dans les phénomènes physiques, symétrie d'un champ électrique et d'un champ magnétique. *J. Phys. Théo. Appliquée* **3**, 393 (1894).

2. Bardeen, J., Cooper, L. N. & Schrieffer, J. R. Microscopic theory of superconductivity. *Phys. Rev.* **106**, 162 (1957).

3. Higgs, P. W. Broken symmetries and the masses of gauge bosons. *Phys. Rev. Lett.* **13**, 508 (1964).

4. Gross, D. J. The role of symmetry in fundamental physics. *Proc. Natl. Acad. Sci. USA* **93**, 14256 (1996).

5. Martin, F. et al. Single photon-induced symmetry breaking of H₂ dissociation. *Science* **315**, 629 (2007).

6. González-Castrillo, A., Palacios, A., Bachau, H. & Martín, F. Clocking ultrafast wave packet dynamics in molecules through UV-induced symmetry breaking. *Phys. Rev. Lett.* **108**, 063009 (2012).

7. Palacios, A., Rivière, P., González-Castrillo, A. & Martín, F. XUV lasers for ultrafast electronic control in H₂. in *Ultrafast Phenomena in Molecular Sciences* Vol. 107 (eds de Nalda, R. & Bañares, L.) (Springer, 2014).

8. Ridente, E. et al. Femtosecond symmetry breaking and coherent relaxation of methane cations via X-ray spectroscopy. *Science* **380**, 713 (2023).

9. Pertot, Y. et al. Time-resolved X-ray absorption spectroscopy with a water window high-harmonic source. *Science* **355**, 264 (2017).

10. Li, M. et al. Ultrafast imaging of spontaneous symmetry breaking in a photoionized molecular system. *Nat. Comm.* **12**, 4233 (2021).

11. Haase, D., Manz, J. & Tremblay, J. C. Attosecond charge migration can break electron symmetry while conserving nuclear symmetry. *J. Phys. Chem. A* **124**, 3329 (2020).

12. Livshits, E. et al. Symmetry-breaking dynamics of a photoionized carbon dioxide dimer. *Nat. Comm.* **15**, 6322 (2024).

13. Bunker, P. R. Forbidden transitions in homopolar isotopically unsymmetric diatomic molecules and the dipole moment of HD. *J. Mol. Spectrosc.* **46**, 119 (1973).

14. Reinhold, E., Hogervorst, W. & Ubachs, W. Complete g–u symmetry breaking in highly excited valence states of the HD molecule. *Chem. Phys. Lett.* **296**, 411 (1998).

15. Hölsch, N. et al. Ionization and dissociation energies of HD and dipole-induced g/u-symmetry breaking. *Phys. Rev. A* **108**, 022811 (2023).

16. Rolles, D. et al. Isotope-induced partial localization of core electrons in the homonuclear molecule N₂. *Nature* **437**, 711–715 (2005).

17. Vacher, M., Boyer, A., Loriot, V., Lépine, F. & Nandi, S. Few-femtosecond isotope effect in polyatomic molecules ionized by extreme ultraviolet attosecond pulse trains. *J. Phys. Chem. A* **126**, 5692 (2022).

18. Gong, X. et al. Attosecond delays between dissociative and non-dissociative ionization of polyatomic molecules. *Nat. Comm.* **14**, 4402 (2023).

19. Ertel, D. et al. Influence of nuclear dynamics on molecular attosecond photoelectron interferometry. *Sci. Adv.* **9**, adh7747 (2023).

20. Fransén, L., Tran, T., Nandi, S. & Vacher, M. Dissociation and isomerization following ionization of Ethylene: insights from nonadiabatic dynamics simulations. *J. Phys. Chem. A* **128**, 1457 (2024).

21. Ludwig, A. et al. Ultrafast relaxation dynamics of the ethylene cation C₂H₄. *J. Phys. Chem. Lett.* **7**, 1901 (2016).

22. van Tilborg, J. et al. Femtosecond isomerization dynamics in the ethylene cation measured in an EUV-pump NIR-probe configuration. *J. Phys. B At. Mol. Opt. Phys.* **42**, 081002 (2009).

23. Ibuki, T., Cooper, G. & Brion, C. E. Absolute dipole oscillator strengths for photoabsorption and the molecular and dissociative photoionization of ethylene. *Chem. Phys.* **129**, 295 (1989).

24. Hervé, M. et al. Ultrafast dynamics of correlation bands following XUV molecular photoionization. *Nat. Phys.* **17**, 327 (2021).

25. Murillo-Sánchez, M. L. et al. Femtosecond XUV-IR induced photodynamics in the methyl iodide cation. *New J. Phys.* **23**, 073023 (2021).

26. Vismarra, F. et al. Few-femtosecond electron transfer dynamics in photoionized donor-π-acceptor molecules. *Nat. Chem.* **16**, 2017 (2024).

27. Felker, P. M. & Zewail, A. H. Rates of photoisomerization of trans-stilbene in isolated and solvated molecules: experiments on the deuterium isotope effect and RRKM behavior. *J. Phys. Chem.* **89**, 5402 (1985).

28. Fdez Galván, I. et al. OpenMolcas: from source code to insight. *J. Chem. Theo. Comput.* **15**, 5925 (2019).

29. Li Manni, G. et al. The OpenMolcas web: a community-driven approach to advancing computational chemistry. *J. Chem. Theory Comput.* **19**, 6933 (2023).

30. Barbatti, M. et al. Newton-X: a surface-hopping program for nonadiabatic molecular dynamics. *Wiley Interdiscip. Rev. Comput. Mol. Sci.* **4**, 26 (2014).

31. Wigner, E. On the quantum correction for thermodynamic equilibrium. *Phys. Rev.* **40**, 749 (1932).

32. Lucchini, M. et al. Few-femtosecond C₂H₄ internal relaxation dynamics accessed by selective excitation. *J. Phys. Chem. Lett.* **13**, 48 (2022).

33. Boyer, A. et al. Ultrafast vibrational relaxation dynamics in XUV-excited polycyclic aromatic hydrocarbon molecules. *Phys. Rev. X* **11**, 041012 (2021).

34. Delor, M. et al. Directing the path of light-induced electron transfer at a molecular fork using vibrational excitation. *Nat. Chem.* **9**, 1099 (2017).

35. Loriot, V. et al. Attosecond metrology of the two-dimensional charge distribution in molecules. *Nat. Phys.* **20**, 765 (2024).

36. Loriot, V., Gitzinger, G. & Forget, N. Self-referenced characterization of femtosecond laser pulses by chirp scan. *Opt. Exp.* **21**, 24879 (2013).

37. Tully, J. C. Molecular dynamics with electronic transitions. *J. Chem. Phys.* **93**, 1061 (1990).

38. OpenMolcas v22.10. https://gitlab.com/Molcas/OpenMolcas/-/releases/v22.1web0 (October, 2022).

39. Hammes-Schiffer, S. & Tully, J. C. Proton transfer in solution: molecular dynamics with quantum transitions. *J. Chem. Phys.* **101**, 4657 (1994).

40. Merritt, I. C. D., Jacquemin, D. & Vacher, M. Nonadiabatic coupling in trajectory surface hopping: how approximations impact excited-state reaction dynamics. *J. Chem. Theory Comput.* **19**, 1827 (2023).

41. Granucci, G. & Persico, M. Critical appraisal of the fewest switches algorithm for surface hopping. *J. Chem. Phys.* **126**, 134114 (2007).

42. Roos, B. O., Taylor, P. R. & Sigbahn, P. E. A Complete Active Space SCF Method (CASSCF) using a density matrix formulated super-Cl approach. *Chem. Phys.* **48**, 157 (1980).

43. Ibele, L. M., Memhood, A., Levine, B. G. & Avagliano, D. Ab initio multiple spawning nonadiabatic dynamics with different CASPT2 flavors: a fully open-source Py Spawn/OpenMolcas interface. *J. Chem. Theo. Comput.* **20**, 8140 (2024).

44. Roos, B. O., Lindh, R., Malmqvist, P.-Å., Veryazov, V. & Widmark, P.-O. Main group atoms and dimers studied with a new relativistic ANO basis set. *J. Phys. Chem. A* **108**, 2851 (2004).

45. Aquilante, F., Gagliardi, L., Pedersen, T. B. & Lindh, R. Atomic cholesky decompositions: a route to unbiased auxiliary basis sets for density fitting approximation with tunable accuracy and efficiency. *J. Chem. Phys.* **130**, 154107 (2009).
## Acknowledgements

We thank Emilien Prost for assistance during experiments. S.N. thanks Center National de la Recherche Scientifique (CNRS) and Fédération de Recherche André Marie Ampère, Lyon for financial support. We acknowledge support from Agence National de la Recherche (ANR-20-CE29-0021, ANR-20-CE29-0021 and ANR-21-CE30-0052). The simulations in this work were performed using HPC resources from CCIPL (Le center de calcul intensif des Pays de la Loire) and from GENCI-IDRIS (grant 2021-101353). The project is partly funded by the European Union (ERC, 101040356, ATTOP, M.V. and A.N.N.). Views and opinions expressed are those of the authors only and do not necessarily reflect those of the European Union or the European Research Council Executive Agency. Neither the European Union nor the granting authority can be held responsible for them. A.N.N. and M.V. wish to thank Anthony Ferté and Lina Fransén for stimulating discussions.

## Author contributions

A.B., Y.H., V.L., F.L., and S.N. carried out the experiment. S.N. analyzed the experimental data. The theoretical calculations were performed by A.N.N. under the supervision of M.V. A.N.N., A.B., F.L., M.V., and S.N. interpreted the results. A.N.N., M.V., and S.N. wrote the manuscript, which all authors discussed. S.N. conceived and led the project.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s42004-025-01621-z.

**Correspondence** and requests for materials should be addressed to Morgane Vacher or Saikat Nandi.

**Peer review information** *Communications Chemistry* thanks Patricia Vindel-Zandbergen, Amit Samanta and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. A peer review file is available.

**Reprints and permissions information** is available at http://www.nature.com/reprints

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

**Open Access** This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, 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 you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. 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-nc-nd/4.0/.

© The Author(s) 2025