Abstract
1 Introduction
Skyrmions1 are topological defects in vector fields, which can be created, moved, and annihilated. They are ubiquitous and important in condensed matter physics.2,3 Recently, the idea of skyrmions was introduced into the communities of photonics and optics.4 The optical analogies of skyrmions have found advanced applications in various fields such as topological Hall devices5 and deep-subwavelength microscopy.6 An optical skyrmion can be realized when the values of electromagnetic vector fields4,6
High harmonic generation (HHG) is a typical highly nonlinear process in ultrafast science,14 through which the fundamental driving photons can be up-converted to a wide spectral region from vacuum ultraviolet to the X-ray domain.15 In HHG, the local states of driving laser fields are inherently resolved with sub-nanoscale spatial resolution of light–atom interaction. This nature promotes the emerging HHG studies to produce novel spatially structured light beams,16
In this work, we show how to create and control EUV free-space optical skyrmions. By employing HHG as the highly nonlinear interaction process and full Poincaré beams with different wavelengths as the driving fields, we demonstrate that the generated EUV fields reveal explicit skyrmionic structures in Stokes vectors. Their skyrmion numbers are determined by monitoring the polarity and vorticity of Stokes vector distributions, which manifest a close correction with the harmonic orders. We reveal that the underlying mechanism is related to spatial-resolved spin constrain in HHG. More importantly, we further propose, by employing the concept of hybrid-order Poincaré sphere pairs (HyOPSP), a novel control scenario for the topological texture of the EUV skyrmions. In this way, we show that with controlling the geometrical parameters of the driving full Poincaré beams, one can form arbitrary crucial EUV skyrmionic textures (e.g., Bloch-type, Néel-type, anti-type, and higher-order skyrmions) in HHG. This controlling strategy can be extended to other nonlinear optical processes.
Sign up for Advanced Photonics Nexus TOC. Get the latest issue of Advanced Photonics Nexus delivered right to you!Sign up now
2 Results
2.1 Signatures of the Driving Structured Fields
The full Poincaré beam is a typical vector field, whose polarization state covers the full surface of the Poincaré sphere. It is constructed by an arbitrary-order Laguerre–Gaussian beam and a fundamental Gaussian mode with an orthogonal polarization state. Its electric field can be given by: , where is the time, is the angular frequency, denotes the spatial coordinate for the incident driving fields [the notation “” indicates the variables in the incident plane, as shown in the inset of Fig. 1(a)], () is the spatial complex amplitude, and () is the topological charge ( or equals zero). Without loss of generality, () can be approximated as a complex constant for the incident laser. For conciseness, we introduce a geometrical parameter defined as , . As shown later, directly determines the topology of EUV fields.
Figure 1.(a) Schematic of generating EUV skyrmions with HHG. A synthesized two-color (
To demonstrate our scheme, we start by considering one typical kind of full Poincaré beams, namely, cylindrical vortex vector beams (CVVBs). We illustrate the physical scenario in Fig. 1(a), where the incident field is spatiotemporally superposed by two CVVBs with different wavelengths (). The geometrical parameters of the 800- and 400-nm CVVBs are given by and . The notations and represent the angular frequencies of 800 and 400 nm, respectively. As shown in Fig. 1(a), for each CVVB, it carries both radial polarization state (i.e., optical polarization singularity) and spiral phase (optical phase singularity).
In HHG, atoms directly interact with the focal field of the incident laser. It is important to know the signatures of focal field distribution. We, therefore, simulate the focal field of CVVB using the Richards–Wolf vectorial diffraction method.25 We use the 20th order Gauss–Legendre integral formula to numerically calculate the integral, and the grid size is chosen as in the focal plane with a spacing of . In our simulation, the numerical aperture (NA) is selected to be . It is a typical focusing condition of HHG experiments. In this way, we can obtain the focal intensity distributions of both 800- and 400-nm fields [upper panel in Fig. 1(b)]. One can notice that the focal intensity distributions manifest as Gaussian-like structures, which differ from the typical donut-shaped structure of individual phase or polarization singularity. Such exotic distribution largely enriches the highly nonlinear interaction in the spatial domain. Furthermore, we calculate the normalized spin angular momentum (SAM) density of focal fields.26,27 As shown in the bottom panel of Fig. 1(b), spatial-varying polarization is generated. It is interesting to find the presence of circular polarizations in the focal plane despite the incident light being linearly polarized throughout. Here, we emphasize that both the peculiar intensity distribution and SAM distribution of the focal field originate from spin-orbit interaction of light.28 The underlying origin of such spin-orbit interaction comes from the time-varying polarization distribution in the wavefront of incident CVVB.28,29 When the 800- and 400-nm CVVBs are spatiotemporally synthesized at focus, three-leaf polarizations show up [Fig. 1(c)]. This polarization has played an important role in the studies of HHG, which allows us to emit circularly polarized EUV light in HHG.30 Most of the relevant works only involve laser fields with a spatially homogeneous three-leaf polarization. In our case, this two-color synthesized field, known as a Lissajous beam,31,32 carries a versatile polarization state, and the polarization varies rapidly in space.
2.2 Generation of EUV Skyrmions with CVVBs
To explore the HHG driven by the synthesized fields, we perform numerical simulation by solving the Schrödinger equation for hydrogen within strong-field approximation.33 The high harmonic field with angular frequency is calculated from the Fourier components of the dipole moment, . Here, the dipole moment depends on both space () and time (), given by , in which is an arbitrary small positive constant, is the saddle point momentum calculated by saddle point equation , is the quasi-classical action, and contribute to the photoelectron recombination and ionization process, is the focal electric field obtained by the Richards–Wolf vectorial diffraction method, and is the vector potential of light. We insert the focal electric field into the expression of the HHG dipole moment and obtain the spatially resolved HHG field of different orders by changing the harmonic frequency .
In order to elucidate the advantages of the two-color scenario in generating EUV skyrmions, we compare the HHG driven by three different fields, i.e., the 800-nm CVVB, 400-nm CVVB, and their synthesized two-color light field. In simulation, the peak intensities of each driving field are taken to be , and their temporal envelopes are selected to be a function of duration time. In Figs. 2(a)–2(c), we show the HHG spectra for different drivers. Due to the joint restriction taken by parity, energy, and momentum conservations in HHG, harmonic spectra reveal the orders of ()th, ()th, and th for the three different drivers, respectively.
Figure 2.Spatio–spectral signatures of HHG driven by the CVVBs. (a)–(c) High harmonic spectra. (d)–(f) Spatial intensity distributions of logarithmic scales. (g)–(i) Spatial distributions of normalized SAM densities along the
In time domain, the generated field in HHG is spatially structured attosecond pulse trains. To identify the real-space topology of the harmonics, we confine one specific frequency of harmonics and extract its spatial intensity distributions [Figs. 2(d)–2(f)]. It can be noticed that when the driver is a single-color CVVB, the harmonics manifest as annulus structures. By contrast, the intensity of harmonics with the two-color driving field is different, which reveals a Gaussian-like structure for ()th and ()th harmonic or annulus structure for th harmonic. Due to the peculiar polarization texture of the driving two-color beam, the th harmonic has significantly lower yield than the ()th harmonics. In Figs. 2(g)–2(i), we show the corresponding SAM distributions. For each harmonic, the SAM density is spatially varied. Upon closer inspection, in the case of the two-color CVVB driver, the SAM distribution of produced harmonics reveals a sixfold symmetrical distribution [Fig. 2(i)]. This unique distribution can be attributed to the inherent sixfold symmetry of the driving field, which is illustrated by the red profiles depicted in Fig. 1(c).
To form a skyrmionic field, local SAM should orient oppositely at the center of the field compared to the surrounding. Our simulation indicates that one cannot obtain an explicit EUV skyrmionic field with single-color CVVB drivers in HHG. This is because the recollision of electrons is rigidly forbidden in a single-color circularly polarized local field,15 resulting in an extremely low photon yield at the center of the harmonic beam [as shown in Figs. 2(d) and 2(e)]. Consequently, the spatial spin reversal is not complete. By contrast, such obstacles can be overcome in the ()th harmonics with a two-color driver. In this case, the yield at the center is substantial, and the SAM is oriented oppositely at the center compared with the surrounding, leading to a non-trivial real-space topology.
To prove that the created ()th harmonics from two-color CVVB fields are skyrmions, we then analyze their spatial vector texture. Conventionally, skyrmions are characterized by skyrmion number , defined as34, where and denote the vector fields and the defined region for constructing skyrmions in the plane. The vector fields of a two-dimensional skyrmion can be written as , , where and are polar coordinates in the focal plane, and and are spherical coordinates in the Poincaré sphere, as shown in Fig. 3(a). In this way, reads
Figure 3.Characterization of the EUV skyrmions in HHG. (a) Unit Poincaré sphere defined by the normalized Stokes parameters,
The skyrmion number is divided into two integers, and , in which describes the polarity of skyrmions, and indicates the vorticity of skyrmions, respectively. In order to create a skyrmionic field, and should both be non-zero. That is, the field vector has to be reversed, comparing the center with its confinement , and, meanwhile, it should also cover the unit Poincaré sphere as increases from to .
In our configuration, the focal light field can be approximated as a paraxial beam, and thus we explore the possibility of generating EUV fields with a skyrmionic structure in Stokes vector fields, . Each parameter of corresponds to an expectation value of a Pauli matrix for a photonic spin state. The skyrmionic boundary is selected to be , in which the focal 800- and 400-nm driving CVVBs are circularly polarized. We then inspect the topological character of the skyrmions by the stereographic projection of Stokes vectors [Figs. 3(b)–3(d)]. One can notice that the Stokes vectors point to the opposite direction at and , and their orientations change in a vortex type [Figs. 3(c) and 3(d)]. To obtain the specific topological charge, , we inspect the azimuthal and radial variation of the Stokes vector distribution [Fig. 3(b)]. It can be seen that the direction of the Stokes vector is changed from downward (upward) to upward (downward) for the 13th (14th) harmonic. This indicates that the polarities of the 13th and 14th harmonics are given by and , respectively. Besides, their transverse components rotate twice counterclockwise from to , meaning the vorticity of the skyrmion is given by . Taking both the polarity and vorticity of skyrmions into account, the skyrmion numbers can be determined to be and for these two high harmonics, respectively. Such EUV skyrmionic fields can be categorized as two-order skyrmions.3
It is interesting to note that the skyrmion number is dependent on harmonic order. In the picture of absorbing photons, the allowed photon absorption channel for specific harmonics can be written by (), in which and denote the number of photons absorbed from 800- and 400-nm CVVBs. Since the driving field is spatially structured, the allowed photon channel is spatially varied. However, the allowed photon absorption channels cannot be solely determined at the positions in , but the channel is well known at and . There, the channel is given by (5, 4) for the 13th harmonic. That is, the number of absorbed 800-nm photons is one more than the number of 400-nm photons. Therefore, the local spin states at and are the same as that of the 800-nm photons for the 13th harmonic [Fig. 3(e)]. Similarly, the corresponding photon channel of the 14th harmonic is (4, 5), and thus the local spin is the same as the 400-nm photon at and [Fig. 3(f)]. As shown in Fig. 1(b), the local spin of 800- and 400-nm CVVBs is inverse, leading to opposite polarities of the 13th and 14th EUV skyrmions. On the other hand, the vorticity of EUV skyrmions is influenced by the focal phase distribution of the 800- and 400-nm CVVBs, which is independent of the harmonic order. Hence, the vorticities are the same for the 13th and 14th skyrmions. As a result, the skyrmion numbers of the 13th and 14th skyrmions are opposite.
2.3 Complete Control Over the Texture of EUV Skyrmions
The behaviors and stabilization of skyrmions highly depend on their topological texture, and thus we then demonstrate how to control the topology of EUV skyrmions. We extend the driving fields from CVVBs to general full Poincaré beams, and we also resort to the synthetization of two-color (800 and 400 nm) beams. To this end, we start by giving an analytical expression for the focal field of full Poincaré beams as
where ( is the wavenumber, is the focal distance) and are introduced by optical focusing. modulates the complex amplitude of focal field, and it is expressed as
With Eq. (2), we then propose two empirical restriction rules for the geometrical parameters of 800- and 400-nm incident full Poincaré beams, which ensure the generated EUV fields in HHG exhibit skyrmionic characteristics. (i) Considering the spin constraints in HHG, the local SAMs of 800- and 400-nm full Poincaré beams should have opposite orientations in the focal plane. It requires that the third components of their Stokes vectors ( and ) exhibit opposite signs, leading to
(ii) The transverse components of Stokes vectors (i.e., and ) determine the direction of the long-axis of the focal electric field ellipse. To create regular spatial polarization distributions for harmonics, and need to be the same as and , respectively. According to Eq. (2), and are determined by the relative phase between right- and left-handed components, i.e., . Thus, the two-color full Poincaré beams are further confined by
Confined by the above rules, for each 800-nm full Poincaré beam, there is only one qualified 400-nm full Poincaré beam. Taking an 800-nm full Poincaré beam with and as an example, the topological charges of a qualified 400-nm beam should be and . Then, to visualize our controlling strategy of EUV skyrmions, we employ the concept of HyOPSP, which is comprised of two hybrid-order Poincaré spheres35 at different wavelengths. Each point on the sphere represents one specific full Poincaré beam. In Figs. 4(a) and 4(b), we display the corresponding HyOPSP, in which , , and are the redefined Stokes parameters serving as the sphere’s Cartesian coordinates. and are the azimuthal angle and pitching angle, respectively. We colorize some points on these two spheres (, for 800 nm and , for 400 nm). The points in the same colors construct the qualified two-color full Poincaré beams which can generate EUV skyrmions in HHG. For each point on the 800-nm hybrid-order Poincaré sphere, one can find a unique point on the 400-nm sphere that meets the above restriction rules [Eqs. (4) and (5)].
Figure 4.Control over the topological texture of EUV skyrmions. (a), (b) HyOPSP formed by an 800-nm Poincaré beam (
Then, we select some typical points on the HyOPSP to demonstrate control over topological textures of skyrmions. At point , the geometrical parameters of 800- and 400-nm full Poincaré beams are and , respectively. As shown in Figs. 4(c) and 4(d), the HHG simulations demonstrate that the generated EUV fields manifest as Bloch-type skyrmions, whose skyrmion numbers are (, ) for the 13th harmonic and (, ) for the 14th harmonic. Through changing the positions on the HyOPSP, the topological texture can be adjusted. At the points , the geometrical parameters are given by and , and the resulted harmonic field is revealed as the Néel-type skyrmions of (, ) [Fig. 4(e)]. Besides the Bloch-type and Neel-type skyrmions, one can also shape the EUV skyrmions with arbitrary topological textures, such as anti-skyrmion of (, ) [ and ] [Fig. 4(g)] and third-order skyrmion of (, ) [ and ] [Fig. 4(h)]. One can notice that the topology of EUV skyrmions in HHG is closely related with the geometrical parameters of the incident laser. The first two parameters of control the polarity and vorticity of EUV skyrmions, while the last two parameters of influence the skyrmionic texture (such as Néel-type and Bloch-type).
For comparison, if we select two-color full Poincaré beams that do not satisfy the restriction rules of , e.g., point on the 800-nm sphere [] and point on the 400-nm sphere [], one can see that the formed Stokes vectors of the EUV field distribute desultorily, and the HHG field is not a skyrmion [Fig. 4(f)]. Overall, the texture of EUV skyrmions can be controlled using the driving two-color full Poincaré beams with qualified geometrical parameters. Such a robust methodology shows its uniqueness in free-space optical skyrmions, which has never been unveiled for the optical skyrmions in guided modes or surface plasmon.
3 Discussion and Conclusion
This work demonstrates the generation and control of EUV skyrmions. On the one hand, the generation of EUV skyrmions is rooted in the transfer of photon angular momentum between driving fields and high harmonics. From the view of electron motion, the transfer of photon orbital angular momentum is realized by the spatially resolved recombination time of electrons in the interaction region. As for the transfer of photon SAM, it is first recorded by the ionized electron trajectory (or electron orbital angular momentum) with respect to nucleus and then taken by high harmonics at the instant of recombination. It should be noted that the spin-orbit interaction of light is critical in modifying the distribution of photon angular momentum in the focal plane, giving rise to the generation of EUV skyrmions in HHG. On the other hand, the control of EUV skyrmions is realized by adjusting the spatial mode of two-color driving fields. In principle, the selection of the two-color laser wavelength affects the spectrum of HHG, but it has little influence on the formation and control of EUV skyrmions. Compared with the single-color scenario, the two-color scenario has significant advantages. It introduces more degrees of freedom that allows one to finely control the spatial structure of EUV skyrmionic fields. Moreover, two-color fields enable the creation of bright circularly polarized HHG, where the spin states of driving fields can be imprinted onto different orders of harmonics. Hence, one can sculpt the spatial mode of three-dimensional Stokes parameters for emitted EUV fields [Figs. 4(c)–4(h)]. Here, the employment of HyOPSP presents a clear picture for the configuration of two-color full Poincaré, and it can be easily expanded to many fields, such as four-wave mixing.
To summarize, we have shown theoretically that HHG makes it possible to generate and manipulate the EUV free-space optical skyrmions. The robust control methodology shows its uniqueness in EUV optical skyrmions, which has never been unveiled for optical skyrmions in surface plasmon polaritons or guided modes. This work presents an important interface between strong-field physics and topology.36,37 The generation and control of EUV skyrmions open a door to access skyrmionic dynamics in the broad ultraviolet region. Magnetic skyrmion Hall effect has been experimentally discovered recently;38 however, to our knowledge, its optical analog has only been studied in theory.6 Since the EUV optical skyrmions reveals a regular spatial distribution of the photon spin-orbital state, its highly nonlinear process and ability in ultrahigh spatial resolution may enable the observation of the optical skyrmion Hall effect. Looking broadly, the EUV skyrmions endow one laser beam with totally opposite chiral responses when interacting with chiral matters, and, thus, it could offer new opportunities in the studies of laser-based spatial separation of enantiomers via photoionization. Furthermore, due to the elaborate spatial polarization structure of the EUV skyrmion, it also has implications in HHG-induced EUV nanoscale imaging.39
Biographies of the authors are not available.
References
[1] T. H. R. Skyrme. A unified field theory of mesons and baryons. Nucl. Phys., 31, 556-569(1962).
[2] Y. Tokura, N. Kanazawa. Magnetic skyrmion materials. Chem. Rev., 121, 2857-2897(2021).
[9] S. Gao et al. Paraxial skyrmionic beams. Phys. Rev. A, 102, 053513(2020).
[17] M. Zurch et al. Strong-field physics with singular light beams. Nat. Phys., 8, 743-746(2012).
[31] D. A. Kessler, I. Freund. Lissajous singularities. Opt. Lett., 28, 111-113(2003).
[32] W. Miao, G. Gbur. Design of Lissajous beams. Opt. Lett., 47, 297-300(2022).
[35] X. Yi et al. Hybrid-order Poincaré sphere. Phys. Rev. A, 91, 023801(2015).
[36] Y. Shen et al. Roadmap on spatiotemporal light fields(2022).
Set citation alerts for the article
Please enter your email address