abstract: We investigate the polarization properties of high harmonics generated with the bichromatic counterrotating circularly polarized (BCCP) laser fields by numerically solving time-dependent Schrödinger equation (TDSE). It is found that, the helicity of the elliptically polarized harmonic emission is reversed at particular harmonic orders. Based on the time-frequency analysis and the classical three-step model, the correspondence between the positions of helicity reversions and the classical trajectories of continuum electrons is established. It is shown that, the electrons ionized at one lobe of laser field can be divided into different groups based on the different lobes they recombine at, and the harmonics generated by adjacent groups have opposite helicities. Our study performs a detailed analysis of high harmonics in terms of electron trajectories and depicts a clear and intuitive physical picture of the HHG process in BCCP laser field.
author: Xiaofan Zhang$^{1}$; Liang Li$^{1}$; Xiaosong Zhu$^{1,}$[^1]; Xi Liu$^{1}$; Qingbin Zhang$^{1}$; Pengfei Lan$^{1}$; Peixiang Lu$^{1,2}$
date: 2024-10-01
title: Helicity reversion in high harmonic generation driven by bichromatic counterrotating circularly polarized laser fields

# Introduction

High-harmonic generation (HHG) occurs in the interaction of intense femtosecond laser field with atoms or molecules . HHG has been intensively studied in experimental and theoretical investigations , because it promises both applications for generating attosecond pulses and for probing the molecular structure and ultrafast dynamics . The HHG process can be nicely explained by the three-step model . First, under the influence of the laser field, an electron escapes from the atom or molecule and reaches the continuum. Then it is accelerated in the laser field. Finally the accelerated electron is driven back to the parent ion and a high energy photon is emitted. According to this picture, only laser fields, which are able to drive the electrons back to the ion, can be used to efficiently generate high harmonics.

Very recently, the bicircualr counter-rotating circularly polarized (BCCP) driving laser field has been a hot spot due to its unique features and applications. A BCCP laser field consists of two coplanar counter-rotating circularly polarized fields with the fundamental field frequency $\omega$ and its second harmonic $2\omega$ . The time-dependent electric field performs a Lissajous figure on its polarization plane with a three-fold symmetry, i.e. it consists of three segments or lobes per cycle. Therefore, the laser field peaks three times one cycle of the fundamental, producing three ionization bursts. The electrons ionized to the continuum near the laser peaks can successfully recombine with the parent ion, radiating an attosecond pulse train with three pulses one cycle of the fundamental .

This attractive laser field has been demonstrated to be an effective tool for fully controlling the high harmonic polarization from linear through elliptical to circular polarization . Importantly, the conversion efficiency for harmonics with arbitrary polarization is comparable to that of the standard HHG process, which yields linearly polarized high harmonics driven by a linearly polarized laser pulse . When the BCCP laser pulse is acting on an isotropic medium, in the frequency domain, a comb of high efficiency harmonic with frequency (3n$\pm$`<!-- -->`{=html}1)$\omega$ is radiated . According to the selection rules , the harmonics with frequency (3n+1)$\omega$ have the same helicity as the fundamental while the (3n-1)$\omega$ harmonics have the same helicity as the second harmonic pulse (3n$\omega$ is parity forbidden). Selecting harmonics with one kind of helicity would generate attosecond pulses with circular polarization. Some efforts have been made for the selection, e.g., optimizing the phase-matching conditions or using crossed beams of counter-rotating circularly polarized driving lasers . This circularly polarized light source is very useful for ultrafast spectroscopy, especially for the ultrafast circular dichroism studies . It was also shown that the bicircular HHG driven by this BCCP laser pulse can be a new method to probe dynamical symmetries of atoms and molecules and their evolution in time . These attractive features of HHG in BCCP laser fields are related to the symmetry of the laser field and the corresponding three-step dynamics of the continuum electrons. To better understand the features of HHG in the BCCP laser field, the correspondence between the harmonic polarization and the electron trajectories needs to be established.

In our study, the polarization properties of high harmonics generated in 2D model Ne atoms by the BCCP laser fields are investigated. We focus on the helicity reversions in the harmonic spectrum, where the rotation direction of the harmonic electric field reverses at specific orders. To explain this phenomenon, we depict the three-step picture of HHG in BCCP fields and reveal the correspondence between the time-frequency distribution of harmonic degree of circular polarization and the classical trajectories. Our study helps people better understand the HHG process and depicts a clear picture for the HHG process in the BCCP laser field.

# Theoretical model

The high harmonic spectrum generated from the interaction between the BCCP laser field and Ne atoms is calculated by numerically solve the two-dimensional time-dependent Schrödinger equation (TDSE) in length-gauge . The 2D model potential of Ne atom is given by (atomic units(a.u.) are used throughout this paper): $$\begin{aligned}
V(\vec{r})= - \frac{Z(\vec{r})}{\sqrt{\vec{r}^{2}+a}},
\end{aligned}$$ where $Z(\vec{r})=1+9\exp(-\vec{r}^{2})$ and a=2.88172 to obtain the ionization potential of Ne atom $\textit{I}_{p}$=0.793 a.u. for 2$\textit{p}$ orbitals. $\vec{r}$$\equiv(x,y)$ denotes the electron position in the two-dimensional x-y plane.

The bicircular electric field is obtained by combining left-circularly polarized and right-circularly polarized laser fields. The laser vector in the $x-y$ polarized plane is defined by $$\begin{aligned}
E(t)=\nonumber E_{0}f(t)(\cos[\omega t]+\cos[\gamma\omega t])\hat{x} \\ +E_{0}f(t)(\sin[\omega t]-\sin[\gamma\omega t])\hat{y},
\end{aligned}$$ where $\gamma$ denotes the integer multiple of the fundamental frequency. The laser field has $\gamma+1$ laser lobes and reflects $\gamma+1$ fold symmetry. The $\omega$ field with helicity +1 rotates counterclockwise while the second field $\gamma\omega$ with helicity -1 rotates clockwise. $f(t)$ is the trapezoidal envelop with 2 cycle rising and 2 cycle falling edges and 5 cycle plateau (in units of fundamental). We set the same laser intensity between the $\omega$ field and the $\gamma\omega$ field. $E_{0}$ is the amplitude of the laser field.

We use the split-operator method to solve Eq.(1). To avoid spurious reflections from the spatial boundaries, the electron wave function $\Psi(\vec{r},t)$ is multiplied by a "mask function" $g$ with the following form $$\begin{aligned}
g=g(x)g(y)
\end{aligned}$$ at each time step, in which $$\begin{aligned}
g(x)=\{^{1\qquad\qquad\qquad\ \ \ |x|<R_{x}-L_{x}} _{\sin^{1/8}(\frac{R_{x}-|x|}{2\pi L_{x}})\qquad|x|\geq R_{x}-L_{x}}
\end{aligned}$$ and $$\begin{aligned}
g(y)=\{^{1\qquad\qquad\qquad\ \ \ |y|<R_{y}-L_{y}} _{\sin^{1/8}(\frac{R_{y}-|y|}{2\pi L_{y}})\qquad|y|\geq R_{y}-L_{y}}.
\end{aligned}$$ For all results reported here, we set the width of the "absorbing" area as $L_{x}$=60 a.u. and $L_{y}$=60 a.u.. In our calculation, we set the range of the "mask function" from -240 a.u. to 240 a.u. in the Cartesian coordinates, i.e., $R_{x}$=$R_{y}$=240 a.u.. The space step is 0.1 a.u.. The time step is $\Delta t=0.05$ a.u..

The time-dependent dipole acceleration can be obtained by means of the Ehrenfest theorem $$\begin{aligned}
A_{q}(t)=-<\Psi(t)|[H(t),[H(t),q]]|\Psi(t)>, q=x,y .
\end{aligned}$$

Then, the harmonic spectrum is obtained by Fourier transforming the dipole acceleration $$\begin{aligned}
a_{q}(\omega)=\int A_{q}(t)\exp(-iq\omega t)dt.
\end{aligned}$$ The intensity of left- and right- polarized harmonic components can be obtained by $$\begin{aligned}
D_{\pm}=|a_{\pm}|^{2},
\end{aligned}$$ where $a_{\pm}=\frac{1}{\sqrt{2}}(a_{x}\pm ia_{y})$. The polarization state of high harmonics can be characterized by the degree of circular polarization (DCP) $$\begin{aligned}
 \zeta=\frac{D_{+}-D_{-}}{D_{+}+D_{-}},
\end{aligned}$$ which can be obtained from the intensity measurements. The ellipticity $\epsilon$ has a one-to-one correspondence with the DCP as $\zeta=\frac{2\epsilon}{1+\epsilon^2}$ . The two parameters both vary in the interval from -1 to 1 and have the same sign.

The rotation direction of the electric field can be quantitatively described by the helicity, which is defined as: $$\begin{aligned}
h=sgn(\epsilon)=sgn(\zeta).
\end{aligned}$$ The helicity $h$ takes the values -1 and +1, indicating the two opposite rotation directions.

# Results and discussions

Firstly, we present the harmonic intensity and the DCP as a function of the harmonic order with different laser parameters in Fig.1. We choose the example of Ne atoms with the 2$\textit{p}$ ground state. In this case, the two components of the bicircular field with $\gamma$=2 have equal intensities. According to the selection rules, the harmonics 3n$\omega$ are absent, and the harmonics (3n+1)$\omega$ and (3n-1)$\omega$ have opposite helicities. The harmonic spectra for left-circularly polarized harmonics (3n+1)$\omega$ and right-circularly polarized harmonics (3n-1)$\omega$ are presented with light red lines and blue lines respectively. The envelopes of these harmonic spectra are presented with the dark red solid lines and blue dashed lines correspondingly. For equal intensities of the two laser field components, the position of the harmonic cutoff is obtained by $$\begin{aligned}
n\omega= \frac{3.17U_{p}}{\sqrt{2}}+1.2|I_{p}|.
\end{aligned}$$ The ponderomotive energy of the entire field is defined by $\textit{U}_{p}=\frac{2E_{0}^2}{4\omega^2}+\frac{2E_{0}^2}{4(2\omega)^2}=\textit{U}_{p1}+\textit{U}_{p2}$. The cutoff obtained by Eq.(11) is almost coincident with that from the simulated results.

The harmonic spectrum driven by the bicircular laser field with the fundamental laser wavelength $\lambda_{1}=1300$ nm and laser intensity $\textit{I}_{1}=1.0\times 10^{14}$ W/cm$^{2}$ is shown in Fig. 1(a). From this figure, one can see that, the intensities of the harmonics (3n-1)$\omega$, having the helicity $h=-1$, are different from that of the harmonics (3n+1)$\omega$, for which $h=+1$. There are two helicity reversions where the red-solid and blue-dashed lines intersect. One is at about the ionization threshold (22nd-order harmonic), the other is at about 107th-order harmonic (labelled by the dash-dot line). Below the ionization threshold, the intensities of (3n+1)th-order harmonics are higher than those of (3n-1)th-order. Between the ionization threshold 22nd-order and the 107th-order harmonic, the intensities of (3n-1)th-order harmonics turn to be higher than those of the (3n+1)th-order. While between 107th-order and the cutoff 120th-order obtained by the Eq.(11), an opposite result to that between the ionization threshold 22nd-order and the 107th-order harmonic order appears. That is, the intensities of the harmonics (3n+1)th-order co-rotating with the $\omega$ field are higher than those of the (3n-1)th-order counter-rotating with the $\omega$ field. In order to clearly show the helicity reversion phenomenon in the harmonic spectra, we display the DCP in Fig.1(e). According to the helicity definition, the helicity reversion appears at the position where the DCP changes its sign. From Fig.1(e), The DCP crosses zero at ionization threshold 22nd-order and 107th-order harmonic, which is coincident with the discussion above.

To check whether the phenomenon is general for the bicircular field interaction with Ne atoms and uncover the reason for the reversion, we apply different laser parameters for fundamental laser wavelength $\lambda_{1}=1300$ nm , 1600 nm and laser intensity $\textit{I}_{1}=1.0\times 10^{14}$ W/cm$^2$, $1.2\times 10^{14}$ W/cm$^2$. The harmonic spectra are presented in Figs.1(b-d) and the DCP are plotted in Figs.1(f-h). It is shown that the intensities of (3n+1)th-order harmonics are higher than those of the (3n-1)th-order harmonics below the ionization threshold, and correspondingly, the values of the DCP are negative. There is a helicity reversion at the ionization threshold for different laser parameters. We have also calculated the harmonic spectra for other model atoms with different two-dimensional soft-core potentials and different ionization potentials. It is shown that the reversion at the threshold is not a general phenomenon and is dependent on the target chosen. Except the helicity reversion at the ionization threshold, the helicity reversion above the threshold is found in all these spectra. The reversion position is labelled by the dash-dot line. In our study, we mainly investigate the helicity reversions in the above-threshold region.

In order to explain the phenomenon of the helicity reversion and elucidate the physics behind it, we present the time-frequency distribution of the harmonic DCP and the classical ionization and recombination paths in Fig.2. The laser parameters are the same as in Fig.1(a). Here we analyze the time-frequency properties of the harmonic DCP with the Gabor transform, which has been proved to be a very powerful tool to analyze the emission times of HHG and to clearly discriminate the electron trajectories. Considering the dipole acceleration $A_{q}(t)$ of Eq.(6), the Gabor transform is defined as

$$\begin{aligned}
GT_{q}[\Omega,t_0]=\frac{1}{2\pi}\int dtA_{q}(t)e^{-i\Omega t}e^{-(t-t_{0})^2/2\sigma^2},
\end{aligned}$$ where $\Omega$ is the frequency of the high harmonics. In our study, we use $\sigma=1/3\omega$. The DCP distribution is obtained by $\zeta=\frac{|GT_{+}|^2-|GT_{-}|^2}{|GT_{+}|^2+|GT_{-}|^2}$, in which $GT_{\pm}=\sqrt{\frac{1}{2}}(GT_{x}\pm iGT_{y})$. In Fig.2(a), we present the DCP distribution of harmonics between 1.5 optical cycle and 8 optical cycle. Due to the symmetry of the BCCP laser pulse, three radiation bursts per fundamental laser pulse cycle can be found in this figure, as predicted in previous works . There are two helicity reversions (labelled by the sold line and the dashed line respectively) and three distinct regions can be identified: () Below threshold 22nd-order (labelled by the solid line), the DCP is mostly negative. The harmonic emissions counter-rotating with the $\omega$ laser field are radiated. () The middle region between the threshold (22nd-order) and the second reversion (about 107th-order, labelled by the dashed line), the DCP is high and positive. The harmonic emissions are emitted with the left-circular polarization. () The upper region between the second reversion (about 107th-order) and the cutoff (120th-order), the DCP is almost negative opposed to the second region. In the second region, the highest DCP is at about the 30th-order harmonic, corresponding to the position of the maximum intensity difference between the left-circularly and right-circularly polarized harmonic components in Fig.1(a). Then the pulse DCP decreases gradually along with the increase of the harmonic order. At about 107th-order denoted by the dashed line, the DCP transits from positive value to negative value. These are directly reflected in Fig.1(e).

To give a clear interpretation of the quantum results above, we explore the classical behavior of the electron path by solving Newton's equation under the BCCP laser field driven . In Fig.2(b), we present the electrons classical ionization times labelled by the blue diamonds, as well as, the recombination times labelled by the red pentagrams. The solid markers are used for the short trajectories and the hollow markers are used for the long trajectories. Note that HHG in bicircular fields is dominated by short trajectories . We focus on the short trajectories of the electrons ionized at one laser field lobe. As shown in this figure, three ionization bursts and three recombination bursts in one fundamental laser cycle reflect the threefold symmetry of the $\omega-2\omega$ counter-rotating circularly polarized laser field . The electrons ionized in one laser lobe can recombine with the parent ion in the following two lobes and two peaks marked as $Q_{1}$ and $Q_{2}$ are generated respectively at laser lobe $P_{1}$ and $P_{2}$ (see Fig.3(b)). The maxima of the peak $Q_{1}$ and peak $Q_{2}$ are 111th-order and 122nd-order harmonics respectively, which approximately correspond to the reversion position and the cutoff position. The slight deviation between the quantum and classical model originates from the quantum effect and the absence of harmonic 3n$\omega$ in quantum model. But it does not influence greatly on our conclusion. By comparing the time-frequency analysis and the classical calculations, one can see that, the DCP between the ionization threshold and about 40th-order harmonic is positive, which is only contributed by $Q_{1}$. However, between the maximum of $Q_{1}$ and the harmonic cutoff, the DCP contributed by only $Q_{2}$ is negative. Besides, at the overlapping part of the two peaks $Q_{1}$ and $Q_{2}$ (40th-order harmonic to 111th-order harmonic), the DCP decreases along with the increase of harmonic order. The correspondence relationship between the time-frequency analysis and the classical calculations indicates that the harmonics generated by $Q_{1}$ are left-circularly polarized while the harmonics generated by $Q_{2}$ are right-circularly polarized. Therefore, the harmonics, which are contributed by the electrons ionized at one field lobe and recombining at the following adjacent lobes, have opposite helicities.

To further illustrate our idea, we present the electron classical trajectories and the plot of the laser field vector of the $\omega-2\omega$ field in Fig.3. The dash-dot green line in Fig.3(a) is the amplitude of the laser field with three peaks per optical cycle. Every peak corresponds to one laser lobe. In order to increase the readability of the figure, the second lobe marked as $P_{1}$ and the third lobe marked as $P_{2}$ are superimposed by the light red background and light blue background respectively. We calculate the short trajectories of the electrons ionized at first laser lobe. It is found that, these trajectories can be classified by the laser lobe they recombine at. That is, the trajectories of the electrons recombining at $P_{1}$ are classified as group $R_{1}$ (the red lines), and the trajectories of the electrons recombining at $P_{2}$ are classified as group $R_{2}$ (the blue dot lines). To give a clear description of the trajectory characteristics, the polar plot of the electric field vector of the $\omega$-$2\omega$ field is shown in Fig.3(b). A clover-leaf structure that has a threefold symmetry is presented. The laser field symmetry is reflected by an attosecond pulse train with three pulses generated in one optical cycle in Fig.2. The electronic field vector at the ionization time at the first laser lobe is marked as bold purple line $d_{i}$. The bold red line at $P_{1}$ and bold blue line at $P_{2}$ denote the laser field vector at the recombination times for electrons of groups $R_{1}$ and $R_{2}$ respectively. The electrons ionized at $d_{i}$ can recombine with the nucleus at $P_{1}$ and $P_{2}$ through the path $R_{1}$ and $R_{2}$. Note that the trajectories $R_{1}$ and $R_{2}$ correspondingly contribute to short trajectories of $Q_{1}$ and $Q_{2}$ in Fig.2(b). Then, based on the discussion about Fig.2, the harmonics radiated from $P_{1}$ and $P_{2}$ have opposite helicities.

Next, we also calculate the relative phase of the $x$,$y$ harmonic components using the strong-field approximation (SFA) model . Fig.4(a) and (b) respectively exhibit the calculation results of the harmonic emissions generated by electrons recombining at $P_{1}$ and $P_{2}$. Based on the classical calculations in Fig.2(b), we choose the recombination time period and integral over a ionized time period for each recombination time chosen. From this figure, one can see that the relative phase of the two components radiated at the laser lobe $P_{1}$ varies between 0.4$\pi$ to $\pi$, which indicates that the radiated harmonics have the positive DCP. However, the relative phase of the two components radiated at the lobe $P_{2}$ changes from 0 to -0.5$\pi$ which indicates that the harmonics have the negative DCP. The result of Fig.4 coincides with that of Fig.2. Therefore, our conclusion that the harmonics contributed by electrons ionized at the first lobe and recombining at the following two lobes have opposite helicities is confirmed by using the classical model and the SFA model.

Finally, we have also calculated the HHG driven by BCCP laser field with $\gamma=3$. The fundamental laser wavelength and laser intensity are $\lambda_{1}=1300$ nm and $\textit{I}_{1}=1.0\times 10^{14}$ W/cm$^2$, respectively. In Fig.5(a), the harmonic spectra for left-circularly polarized harmonics (4n+1)$\omega$ and right-circularly polarized harmonics (4n-1)$\omega$ are presented with light red and blue lines respectively. The envelopes of these harmonic spectra are presented with the dark red solid lines and blue dashed lines correspondingly. In this spectrum, four helicity reversions are obviously observed. The first one is also at the ionization threshold, which is determined by the target. The remaining three reversions are at about 69th-order (the dashed line), 80th-order and 95th-order (the dash-dot line). In order to clearly display the phenomenon of helicity reversion of the harmonics, we plot the time-frequency distribution of the harmonic DCP. Four attosecond pulses are generated per optical cycle, which reflects the fourfold symmetry of the $\omega-3\omega$ laser field. The helicity reversions at about 69th-order and 95th-order are labelled by the dashed line and the dash-dot line. Similarly, the classical electron paths in one optical cycle are exhibited in Fig.5(c). We only study the trajectories of the electrons ionized at one laser lobe and recombining at the following three lobes. The recombination times of these electrons are labelled by the dark red pentagrams. Others are labelled by the light red pentagrams. Three peaks $Q_{1}, Q_{2}$ and $Q_{3}$ are presented in Fig.5(c). Comparing Fig.5(b) with Fig.5(c), it is shown that, the DCP in the region between the ionization threshold 22nd-order and the reversion position 69th-order is positive and decreases with harmonic order. This region is dominated by $Q_{1}$. The region between the reversion position 95th-order (the maximum of peak $Q_{3}$) and the harmonic cutoff 109th-order (the maximum of peak $Q_{2}$) is only contributed by $Q_{2}$. The helicity of $Q_{2}$ is negative, which is opposite to that of $Q_{1}$. The helicity reversions at the 69th-order and the 95th-order in Figs.5(a) and 5(b) appear at the maxima of $Q_{1}$ and $Q_{3}$ in Fig.5(c). In the overlapping region of $Q_{2}$ and $Q_{3}$ between 69th-order and 95th-order, the harmonic emissions have both positive and negative DCP and an additional helicity reversion appears. This will be discussed later. For a clear analysis of the classical electron trajectories in the $\omega-3\omega$ laser field, the laser field vector is plotted in Fig.5(d). The electronic field vector at the ionization time at the first laser lobe is marked as bold purple line $d_{i}$. The recombination times of the short trajectories of $Q_{1}$, $Q_{2}$ and $Q_{3}$ are marked as the red bold line in laser lobe $P_{1}$, cyan bold line in lobe $P_{2}$ and blue bold line in lobe $P_{3}$ respectively. It is found that the peaks $Q_{1}$, $Q_{2}$ and $Q_{3}$ are contributed by the electrons recombining at different laser lobes. Based on the above comparison between Fig.5(b) and Fig.5(c), the harmonic emissions generated in $P_{1}$ have opposite helicities with those in $P_{2}$, and the helicities of harmonics generated in $P_{2}$ are opposite to those in $P_{3}$, which is consistent with the conclusion discussed above for $\omega-2\omega$ laser field.

To understand the appearance of the additional helicity reversion in the overlapping region of $Q_{2}$ and $Q_{3}$, the weights of the contributing electron trajectories should be obtained. Fig. 6 shows the weights of the short trajectories in both the $\omega-3\omega$ and $\omega-2\omega$ laser fields. For each trajectory, the weight is calculated by $W(t_0, v_\perp)=W_0(t_0)W_1(v_\perp)$. $W_0(t_0)$ is the ionization rate at time $t_0$ obtained from the TDSE simulation by calculating the derivative of the time-dependent ionization probability. $W_1(v_\perp)$ is the weight with a Gaussian distribution for the tunnel electron with transverse initial velocity $v_\perp$ . For the $\omega-3\omega$ laser field as shown in Fig. 6(a), the spectrum range can be divided into four regions. In the first region below the 69th-order, the weight of peak $Q_1$ is much greater than those of peaks $Q_2$ and $Q_3$. In the second region from the 69th-order to the 79th-order, peak $Q_1$ disappears and the weight of $Q_2$ is greater than $Q_3$. Then from about 79th-order to 95th-order, the weight of $Q_3$ becomes greater than $Q_2$ in the third region. In the last region above the 95th-order, only peak $Q_2$ remains contributing to the harmonic emissions. In each region, the helicity of the generated harmonics is determined by helicity of the strongest contributing peak, which can well explain the results in Figs. 5(a-b). Owing to the intersection between the weights of peaks $Q_2$ and $Q_3$, the additional helicity reversion in the overlapping region of $Q_{2}$ and $Q_{3}$ appears. For the $\omega-2\omega$ laser field as shown in Fig. 6(b), the spectrum range can be divided by the 111th-order into two regions. The weight of peak $Q_1$ is much greater than that of peak $Q_2$ in the first region. In the second region, only peak $Q_2$ contributes to the harmonic emissions. Since the weights of the two peaks have no intersection, there is no extra helicity reversion in the $\omega-2\omega$ BCCP laser field.

Based on all the discussions above, a picture of the electron classical trajectories in HHG driven by BCCP laser field is revealed. Besides, the correspondence between the harmonic polarization and the classical trajectories has also been established. Specifically, the electron ionized at one laser lobe can recombine with the nucleus at the following laser lobes. These recombining trajectories can be classified into different groups based on the laser lobe they recombining at. Moreover, the harmonics emitted from these groups at adjacent laser lobes have opposite helicities. For instance, in $\omega-2\omega$ laser field, groups $R_{1}$ and $R_{2}$ are both ionized at one laser lobe $d_{i}$. And they recombine with the nucleus at the following laser lobes $P_{1}$ and $P_{2}$ and form the peaks $Q_1$ and $Q_2$, respectively. The harmonics from $Q_{1}$ and $Q_{2}$ are radiated with opposite helicities and the helicity reversion occurs at the maximum of $Q_1$. In $\omega-3\omega$ laser field, electrons ionized at one laser lobe $d_{i}$ recombine at the following lobes $P_{1}$, $P_{2}$, $P_{3}$ and form the peaks $Q_1$, $Q_2$, $Q_3$, respectively. The harmonics from peaks $Q_1$ and $Q_3$ have positive helicity while those from $Q_2$ has negative helicity. The helicity reversion occurs at the position where the dominantly contributing peak is changed to its neighbouring peak, i.e. at the maxima of $Q_1$ and $Q_2$ and at the intersection between the weights of $Q_2$ and $Q_3$.

# Conclusion

In conclusion, under the interaction of bicircular laser field with the 2D model Ne atoms, a harmonic spectrum with different helicities in different spectrum ranges is generated. It is shown that, the helicity of the elliptically polarized harmonic emission is reversed at particular harmonic orders. Based on the time-frequency analysis and the classical three-step model, the relationship between the reversion position and the electron trajectories is obtained. The electrons ionized at one lobe of laser field can recombine with the nucleus at the following laser lobes radiating harmonics. These harmonic emissions generated from adjacent laser lobes have opposite helicities. Besides, we confirm our conclusion through the calculation of relative phase by the SFA model. Our study performs a detailed analysis for the HHG driven by the BCCP laser field in terms of the classical trajectories and established the correspondence between the harmonic polarization and the classical trajectories.

# Acknowledgement

This work was supported by the National Natural Science Foundation of China under Grants No. 11404123, 11234004, 61275126, 11422435 and 11574101. Numerical simulations presented in this paper were carried out using the High Performance Computing Center experimental testbed in SCTS/CGCL (see http://grid.hust.edu.cn/hpcc).

[^1]: Corresponding author: firstname.lastname@example.com


author: L. P. Chitta, S. K. Solanki, H. Peter, R. Aznar Cuadrado, L. Teriaca, U. Schühle, F. Auchère, D. Berghmans, E. Kraaikamp, S. Gissot; C. Verbeeck
date: Received ; accepted
title: Capturing transient plasma flows and jets in the solar corona

Intensity bursts in ultraviolet (UV) to X-ray wavelengths and plasma jets are typical signatures of magnetic reconnection and the associated impulsive heating of the solar atmospheric plasma. To gain new insights into the process, high-cadence observations are required to capture the rapid response of plasma to magnetic reconnection as well as the highly dynamic evolution of jets. Here, we report the first 2 s cadence extreme-UV observations recorded by the 174 Å High Resolution Imager of the Extreme Ultraviolet Imager on board the Solar Orbiter mission. These observations, covering a quiet-Sun coronal region, reveal the onset signatures of magnetic reconnection as localized heating events. These localized sources then exhibit repeated plasma eruptions or jet activity. Our observations show that this spatial morphological change from localized sources to jet activity could occur rapidly on timescales of about 20 s. The jets themselves are intermittent and are produced from the source region on timescales of about 20 s. In the initial phases of these events, plasma jets are observed to exhibit speeds, as inferred from propagating intensity disturbances, in the range of 100 km s$^{-1}$ to 150 km s$^{-1}$. These jets then propagate to lengths of about 5 Mm. We discuss examples of bidirectional and unidirectional jet activity observed to have been initiated from the initially localized bursts in the corona. The transient nature of coronal bursts and the associated plasma flows or jets along with their dynamics could provide a benchmark for magnetic reconnection models of coronal bursts and jets.

# Introduction

Magnetic reconnection is a universal process that is invoked to explain phenomena ranging from astrophysical jets to granular-scale explosive events on the Sun. In the solar atmosphere, signatures of magnetic reconnection are prevalent. In particular, observations of spectral line broadening and double-peaked spectral profiles support the idea of bidirectional jets during magnetic reconnection at the site of megameter-scale explosive events in the transition region at the base of the quiet-Sun corona , as well as ultraviolet (UV) bursts in active regions . High-cadence (1.4 s to 5 s) spectroscopic observations of such bursts show that the observed intensity and broadening of spectral lines fluctuate on timescales as short as 10 s, suggesting that the plasma responds very dynamically to magnetic reconnection .

In chromospheric images that sample temperatures around $10^4$ K, inverse Y-shaped anemone jets are observed in active regions. Such jets outline plasma eruptions on spatial scales of 150 km to 300 km, with speeds of 10 km s$^{-1}$ to 20 km s$^{-1}$ . Signatures of similar small-scale chromospheric jets are also observed at the footpoints of active region coronal loops . Small-scale plasma jets at higher speeds on the order of 100 km s$^{-1}$ are observed at higher temperatures of $\sim$`<!-- -->`{=html}0.1 MK to 1 MK in the transition region and corona . Evidence of magnetic reconnection is further supported by observations that these small-scale jets erupt when patches of small-scale opposite-polarity magnetic elements in the solar photosphere interact and undergo flux cancellation . There is also some evidence for a direct association of explosive events identified in spectroscopic observations and plasma jets observed in imaging observations .

The dot-like, loop-like, and surge- or jet-like EUV brightenings in active regions, observed during the second flight of Hi-C, have provided new insights into the spatial evolution of plasma response to magnetic reconnection by revealing bi- and uni-directional plasma flows at small spatial scales of 2–5. However, the rapid response of high temperature plasma (above 0.1 MK) to magnetic reconnection in terms of the spatial evolution and internal dynamics of jets are not well known due to the lack of extreme ultraviolet (EUV) coronal observations at cadences as high as 1 s to 2 s, which are comparable to high-cadence spectroscopic observations of the transition region . We expect such high-cadence coronal images to reveal the transience of plasma jets as a direct response to the ephemeral nature of magnetic reconnection. We also expect such imaging observations to provide insights into the morphological evolution of jets. In this study, we report the first 2 s high cadence observations of plasma flows or jets using data obtained by the Extreme Ultraviolet Imager on board the Solar Orbiter mission during its cruise phase on 2021 February 23. The reported plasma flow events are not classical or traditional coronal jets as described in the literature .

# Observations

We employ a high-cadence image sequence recorded by the 17.4 nm EUV High Resolution Imager (HRI$_{\textrm{EUV}}$) of the Extreme Ultraviolet Imager on board the Solar Orbiter mission , which is currently in its cruise phase. EUI/HRI$_{\textrm{EUV}}$ records the EUV emission on the central $2048\times2048$ pixel area of its CMOS sensor with an image scale of 0.492 pixel$^{-1}$. The observations used in this study were obtained on 2021 February 23 between UT 17:13:25 and UT 17:20:59 (i.e., total duration of 454 s).[^1] At the time of observations, the sun-bound Solar Orbiter was at a distance of 0.526 astronomical units and hence, the linear image scale of EUI/HRI$_{\textrm{EUV}}$ was about 188 km pixel$^{-1}$. These observations were recorded at 2 s cadence and we used level-2 calibrated data. Using a cross-correlation technique, we removed the remaining jitter in the image sequence.

EUI/HRI$_{\textrm{EUV}}$ samples emission from roughly 1 MK plasma due to Fe ix (at 17.11 nm) and Fe x (at 17.45 nm and 17.72 nm). In Figure , we display the mean emission map (obtained by averaging all the frames in the image sequence) covering the full field of view of HRI$_{\textrm{EUV}}$. The region covers the quiet-Sun corona that is close to the disk center as observed from the vantage point of Solar Orbiter. The Carrington longitude of Solar Orbiter was about 152, while that of Earth was about 292. This means that solar disk center seen from Solar Orbiter could not be seen from Earth and thus no coordinated solar observations from Earth orbit or by ground based telescopes are available.

# Signatures of dynamic magnetic reconnection in the solar corona

The HRI$_{\textrm{EUV}}$ images revealed typical quiet-Sun coronal features such as coronal bright points on scales of 10 Mm interspersed with patches of diffuse EUV emission (Fig. ). Coronal bright points are widely studied features in the literature . The data also revealed a number of transient roundish and loop-like elongated brightenings that are distributed throughout the field of view (not marked in Fig. ). Numerical simulations show that such transient brightenings in the corona can be triggered by magnetic (component) reconnection . There is also observational evidence to support magnetic reconnection as the driver of these transient coronal events based on the detection of double-peaked spectral lines (indicative of bidirectional jets from the site of reconnection) emitted by those structures at transition region temperatures of around 0.1 MK .

These transient events exhibit rapid (internal) dynamics that can be seen as intensity disturbances along the brightenings in the data analyzed here. We focus on four examples of elongated brightenings detected in the 2 s high-cadence EUI data in this study. These appear to be collimated plasma flows, which we refer to as jets in the rest of the paper. The location of these events is marked in Fig.  (boxes A–C). The examples were selected based on their clear visibility, which allows us to track their evolution with minimal confusion coming from the rather diffuse and weak local fore- and background coronal emission.

## Bidirectional jet

Here, we analyze a bidirectional jet event that is a typical signature of magnetic reconnection. In Fig. a-d, we show the spatial morphological evolution of a jet event. The general structure of this jet is discernible also at the center of box-A overlaid on the mean emission map in Fig. . The event started as a weak, spatially localized burst, seen as a compact bright source (indicated by an arrow in Fig. a). The intensity of this compact source increased over a period of 20 s (from UT 17:17:00 to UT 17:17:20; see white-colored light curve of the compact source in panel-e). After this phase of localized brightening, the event rapidly developed into an elongated structure (panels c-d). On one side of this elongation, the flows appears to exhibit complexity and then bifurcate further into multiple segments, diverging away from a Y-type junction (green arrow in panel d). This could reflect the internal magnetic complexity associated with this event. While the jet was visible, the compact source region exhibited further intensity fluctuations on timescales of 10 s to 20 s, pointing to intermittency in the magnetic reconnection. Away from the compact source, the jets are observed to exhibit intensity fluctuations on similar timescales (see cyan coloured light curve of the compact source in panel-e). The propagating nature of disturbances is apparent from the time-lag between the initial intensity rise exhibited by the compact source and that displayed by the jet.

The episodic nature of jet activity and its relation to the initial compact source are demonstrated using the space-time stack plot shown in Fig. e, derived by plotting the intensities along a strip placed along the structure. In particular, intermittent jets are observed to originate on either side of the initial compact source (located between distances of 4 Mm and 6 Mm from the start of the selected strip). These disturbances apparently propagate over distances (i.e., lengths of the jets) of 2 Mm to 4 Mm. The intermittency of jets is further established using smooth-subtracted image of the space-time stack plot (panel-f).[^2]

Initially, these jets exhibit speeds on the order of 100 km s$^{-1}$. Each episode lasts for about 20 s, the timescale that is also observed in the intensity fluctuations of the light curve of the central bright source (panel-e). Assuming that the intensity disturbances are plasma motions, the jet speeds are comparable to the local sound speed of about 150 km s$^{-1}$ in a million Kelvin hot plasma in the corona. The event in its entirety lasted for about 200 s.

## Unidirectional jets

In the following, we analyze intermittent unidirectional jets observed by HRI$_{\textrm{EUV}}$. In Fig. a-d, we present the morphological evolution of another event that exhibits characteristics that are similar to the example in Fig.  (see also box-B in Fig. ). The event in Fig.  also starts as a compact bright source (see arrow in Fig. a). But compared to the example in Fig. , which is more isolated, the event in Fig.  originates in a region of enhanced activity (see bright regions adjacent to the source marked by an arrow in Fig. a). The source exhibits a rapid intensity increase for a period of roughly 20 s (between UT 17:17:20 to UT 17:17:40; see light curve in panel-e). After this phase, repeated plasma eruptions are observed to originate from the source (see panels c, e-f). The propagation speeds of intensity disturbances of about 150 km s$^{-1}$ are comparable to the local sound speed. This event also lasts for about 200 s.

While the intensity disturbances propagate to relatively similar distances on either side of the central source in Fig. f, the detected jets in Fig. f are predominantly in one direction. On one side from the source (located about a distance of 2 Mm from the zero position of the strip), jets have a projected length of nearly 6 Mm. Opposite in direction to these longer jets (shown in Fig. f), there are indications of repeated, but shorter propagating perturbations with projected lengths of only 1 Mm. If these shorter-scale perturbations are also associated with small-scale jet activity, then the difference in the projected lengths of plasma flow in opposite directions can be understood in terms of jets propagating in closed magnetic loops, where the source region is situated closer to one of their footpoints. However, we cannot rule out other plausible explanations, such as magnetohydrodynamic waves at the footpoints of the longer jet driving these intensity perturbations on shorter length-scales of 1 Mm.

Adjacent to the jet activity discussed in Fig. , we observed a pre-existing loop-like plasma feature (green arrow in Fig. c) with one of its apparent footpoints situated closer to the source region marked in Fig. b. This loop-like structure gradually fades away by UT 17:19:20. Its footpoint region exhibited another localized brightening (see Fig. a-d and box-B in Fig. ) that led to plasma eruptions. This burst episode exhibits brightening between UT 17:19:40 and UT 17:20:00, (see Fig. e), after which point a jet is ejected at speeds of about 150 km s$^{-1}$. Unlike jets seen in previous examples, that did not show signatures of speed changes as they propagate, this unidirectional jet exhibits some deceleration after UT 17:20:20. This jet also propagates to distances of about 6 Mm from the initial bright source. There are signatures of subsequent jet activity from this location seen as intensity ridges in the space-time stack plots (between UT 17:20:00 and UT 17:21:00; Fig. e-f). The whole event evolves over a period of about 80 s. The footpoint regions of the events shown in Figs.  and appear to be close to each other (see Fig. b). It is possible that the jet in Fig.  is causally linked to the eruption in Fig. , which is still ongoing when the new jet starts.

An example of a jet that is fainter than the other three is presented in Fig.  (box-C in Fig. ). The bright source becomes clearly distinguishable starting around UT 17:18:20 and a jet is ejected from this region around UT 17:18:40, which then propagates at speeds of about 150 km s$^{-1}$ (Fig. e-f, with at least two additional episodes between UT 17:19:20 and UT 17:20:00). Although the morphology of the jet is quite simple, its source region footpoint exhibits a dynamic evolution. The duration of this event is about 100 s.

# Discussion and conclusions

In the initial phases, the observed transient jets propagate with speeds in the range of 100 km s$^{-1}$ to 150 km s$^{-1}$. Using 5 s cadence EUI data from the 2020-May campaign, investigated the internal dynamics of EUV brightenings and found signatures of propagating disturbances with speeds in the range of 25 km s$^{-1}$ and 60 km s$^{-1}$. This is a factor of two to six smaller than the speeds of jets investigated here. In this regard, the study of complements our work in highlighting a variety of plasma dynamics in coronal structures that could be revealed, thanks to higher cadence EUI data.

The jet speeds reported in our study are on the higher end of the distribution of speeds exhibited by type-II spicules and jetlet-like features from magnetic network lanes observed adjacent to an active region . The speeds, however, are generally on the lower end of transition region network jets (seen as fast propagating intensity disturbances) in coronal hole regions , but comparable to transition region jets observed in the quiet-Sun .

The main difference, however, is that the jets observed in these earlier studies appear to be collimated outflows from the Sun along (locally) open magnetic fields. In contrast, both the bi- and unidirectional jet features we presented here are apparent flows confined to closed magnetic loops in the quiet-Sun . Our quiet-Sun events could be similar to the previously reported active region bi- and unidirectional surges and jet-like features .

Another related difference is the general appearance or spatial morphology of jets itself. Typical coronal jets consist of a bright base along-with a spire that extends outward from the bright jet-base into the corona . However, the collimated flows we investigated here do not show these typical characteristics of classical coronal jets. Why the studied jets differ in terms of spatial morphology in comparison with the classical coronal jets remains to be answered.

A clear imaging of bi-directional plasma jets in the corona was previously reported by using EUV data from the Atmospheric Imaging Assembly on board the Solar Dynamics Observatory . However, the transient and repetitive nature of quiet-Sun jets on shorter timescales of 20 s that we present in all examples here had not been observed before. This is because (as discussed in Appendix ) such short-lived transients could not be observed with routine SDO/AIA EUV observations.

Fast-propagating disturbances in the upper chromosphere and transition region were also detected along bright elongated structures using Lyman-$\alpha$ slit-jaw images recorded by the CLASP sounding rocket experiment . These CLASP Lyman-$\alpha$ propagating disturbances might be transition region counterparts to events similar to those studied here. Whether or not they are visible also in the EUI/HRI Lyman-$\alpha$ Imager requires further investigation with new observations.

Furthermore, a common feature we observed among all four jets is that they erupt from an initial compact bright source. The intensity of the source increases for about 20 s. The jet material is then observed as repetitive eruptions from this local bright source. The short timescales associated with these jets resemble fluctuations observed in the nonthermal broadening of transition region spectral lines (formed at around 0.1 MK) from reconnection-driven microflares that are likely produced by internal turbulent motions .

In three of the four jet events we investigated (i.e., Figs. , , and ), there is no clear signature of acceleration or deceleration of ejected material. Nevertheless, the jets fade away after traversing distances of a few Mm from the source region with almost uniform speeds. One possible explanation is that the jets are further heated or cooled down such that the resulting emission is not captured by the HRI$_{\textrm{EUV}}$ 174 Å filter. In one case (Fig. ), however, the initial jet is observed to slow down. Whether this deceleration reflects the true kinematics of the jet plasma or whether it is an apparent effect (i.e., produced by the loop geometry that curves perpendicular to the plane of sky, as in a downward curved loop) requires further investigation with more such events, coupled with magnetic field extrapolations.

Because these jets are characterized by their intensity enhancement with respect to the local background, it is either the segments of these jets or the whole structures themselves that could be classified as quiet-Sun EUV bursts or campfires , which are likely akin to small-scale bursts detected in active regions . Our spatio-temporal analysis provides further insights into the dynamics and transient nature of these small-scale events.

Magnetic reconnection models of solar bursts predict that the current sheets break into copious magnetic islands or plasmoids that are then ejected outward. These plasmoids then hit the ambient medium and thermalize to produce repeated bursts of observable plasma emission . It is possible that the jets we observed are related to the dynamics of such plasmoids along reconnection current sheets in coronal bursts. The observed speeds of jets, their transient nature (short lifetimes) and their origination from initially localized bursts could be used to test magnetic reconnection models of coronal bursts and jets.

The data we used here are among the highest cadence EUV imaging observations of the solar corona ever recorded. These novel 2 s cadence data revealed the repetitive and transient nature of coronal jets on timescales of about 20 s. There are hints of similar plasma flows (with repetition) also in coronal bright points distributed throughout the field of view. However, most of these regions appear over-exposed in this experimental observational set, limiting the extraction of local flow properties. Similar high-cadence EUV observations (with better exposure settings to capture bright regions) during the perihelion of Solar Orbiter (around 0.3 AU) with roughly two times better spatial resolution than the data used here will further reveal and establish the details of the internal dynamics and the transient nature of these jets and their relation to plasmoids. Combining such observations with photospheric magnetic field information at a similar spatial resolution that is expected from the Polarimetric and Helioseismic Imager on Solar Orbiter will shed light on the magnetic drivers of these transient jets.

[^1]: <https://doi.org/10.24414/z2hf-b008>

[^2]: To enhance the fluctuations on timescales of 20 s observed in the light curves in panel (e), we first smoothened the space-time stack plot in the time domain using a smoothing window of 10 time steps (corresponding to a duration of 20 s). The smoothened image is then subtracted from the original space-time map. The features remain clearly discernible over a range of smoothing windows from 5 to 25 time steps. No smooth-subtraction is done in the spatial direction. This is the case also for the smooth-subtracted images in Figs. , , and . The intermittent nature of jets is evident in the associated online animations.