Optical vortex knots and links via holographic metasurfaces

ABSTRACT Vortices arise in many natural phenomena as dark points of total destructive interference. Sometimes they form continuous lines and even enclosed loops with knotted or linked topologies in three spatial dimensions. Since the mathematical topology was introduced into physics, from hydrodynamics, condensed matter physics to photonics, and other modern physical fields, scientists have been exploring the related topological essences of vortex knots; hence, the topology is a forefront topic in different physical systems. Owing to the reliability and observability of light in free space, optical vortex knots and links are the most studied physical topologies. Here, we review some of these developments with a focus on optical vortex knots and links. We first introduce the brief historical perspective and structural properties of optical vortices. Then, we trace the progress on the theoretically constructing, experimentally generating, and characterizing methods of topological light fields. Wherein, we review recent developments of holographic metasurfaces and their applications in generating ultrasmall optical vortex knots. At last, we envision the possible challenges and prospects of topological light fields.


Introduction
Vortices as phase singularities are a common phenomenon in nature, ranging from water and atmosphere circulation to condensed state in electronic systems, and even spiral galaxies [1][2][3]. In a complex scalar wavefield, vortex singularities are as locations where phases are undefined, and the intensities of the field are zero. They were first postulated by Whewell in the 1830s based on the observation of ocean wave [4], which may be considered as a complex scalar field, that amphidromic points exist in the tides, as the crosses of tides. Beyond the wave phase, in quantum physics, vortices represent quantized flux circulations of some physical quantities as topological defects, determining the properties of quantized systems [5][6][7][8][9]. The nontrivial dynamics of vortices in classical and quantum systems have evoked extensive interests [10][11][12][13][14][15] and brought significant opportunities for the exploration of new wave transmission and the realization of accurate wavefield manipulation, such as directional energy transport [16], topological quantum computation [17], topological information [18] and so on.
The serious study of vortices began in 1974, when Nye and Berry studied the position-dependent amplitude of ultrasound reflected from the ice covering Antarctica [19]. In the recorded intensity and phase distribution of scattered ultrasound reflected from the experimental model, they observed intensity nulls surrounded with changed phase from 0 to 2π and defined such special singularities as 'wave dislocation', analogous to the planes of atoms around dislocations in crystals. In 1989, Coullet firstly used the term 'optical vortex' to describe a possible laser mode having phase singularity [20]. Subsequently, as the recognition of orbital angular momentum (OAM) [21], which manifests itself as an intrinsic property of photon, optical vortices have been intensively studied due to their enormous application potentials [22][23][24][25][26].
Simultaneously, particular interests were traced into three-dimensional (3D) structures of vortex. Since the topology (in mathematics, topology describes the nature of geometric figures or spaces that can remain unchanged after continuously changing shapes) was introduced into the physical domain, various topological configurations, and associated dynamics in disparate fields have been discovered [27][28][29][30][31][32][33][34][35], boosting the emergence of a large number of striking wavefields, and novel topological materials such as topological insulators, topological semimetals, and topological superconductors [36][37][38]. Novel physical phenomena related to topological states are also radiated from condensed electronic systems to fluid, phonons, and other systems [39][40][41][42].
In light field, optical vortices form curved lines, namely, nodal lines, and tangle three-dimensionally to configure topological structures, such as knots and links [3]. Such intriguing spatially structured light fields with phase singularities in real space, exhibit broad application prospects in optical manipulation, optical communication, and quantum information transmission [26,43,44]. Moreover, topology can be considered as the exploration of itself, i.e. 3D structural properties, as a new degree of freedom in photonics that is different from the traditional modulation parameters such as amplitude, phase, and polarization [45]. Knotted and linked light fields, as a new family of null solutions to Helmholtz equation in 3D free space, have attracted increasing interests because of their prominent spatial feature and the unique opportunity that combines the experimental and analytical studies of common topological properties in disparate matter systems [45][46][47][48].
Here, we report on the recent progress of 3D vortex, more specifically, optical vortex knots and links. The contents are organized as follows. Firstly, we begin with a brief historical perspective and structural features of vortex knots and links. Then we trace the progress on the generating and characterizing methods referring to topological light fields. We discuss recent developments in holographic metasurfaces: high efficiency, multi-parameter and wideband response, and their applications in generating ultrasmall optical vortex knots. At last, we give an overview on the outlook of anticipated future expansion of the topological light fields.

Optical vortex
The study of optical vortices can be dated back to the research of Goos-Hänchen shift, that Wolter discovered wave circulation induced by a street of vortices with alternating signs, which was in the phase pattern of interference field composed by the incident and reflected plane waves, as well as evanescent waves [3]. Subsequently, this circulation phenomenon associated with optical vortex was demonstrated to be quite general in light fields as the result of multiple wave interference. For instance, Braunbek recognized that three plane waves were sufficient to produce energy vortices, and found vortices in Sommerfeld's exact solution of a plane wave diffracted by an infinite halfplane. Meanwhile, the exploitation of beams bearing phase singularities was also underway. Remarkably, Vaughan and Willettes investigated the properties of light beams with phase singularities and their interference along the beam axis [49].
Optical vortices, as a ubiquitous phenomenon in interference, naturally occur in the superposition field [3,50]. In three spatial dimensions, phase singularities continuously distribute in the form of tangled vortex lines, that is, threads of darkness in light fields [51]. Take a random speckle as an example, which is the interference pattern generated by the random scattering of a coherent source and can be described as the superposition of a large number of random plane waves with different wavevectors, in which these 3D configurations of the vortex line are mainly divided into the following two categories [52].
The first type, passing through the entire volume as unbounded line, as the red lines shown in Figure 1(a), which exists typically in the superposition field of three plane waves with different wave vectors [3]. Since such kind of vortex lines extends on axis, or parallel to the axis, the topologies are structurally stable during beam propagating. Therefore, researches are mostly focused on the transverse modes with helical phase patterns and the associated energy flow distributions of light beams [53]. In this case, the vortex is characterized by a topological index, namely, topological charge, which depicts the change of phase φ along a circuit C enclosing the phase singularity in the lateral plane, quantized in units of 2π: 111 (1) where l is the so-called topological charge of optical vortex. Inspired by the exploitation of Vaughan and Willets [20], cylindrically symmetric higher-order beams, as transverse eigenmodes of laser cavity, were discovered carrying inherent optical vortices on axis. Thus, the topological charge is an essential parameter to characterize optical vortex and the carrier beam. Laguerre-Gauss (LG) beams and Bessel-Gauss (BG) beams [21,54], the socalled vortex beams, are examples. Figure 2 displays the transverse amplitude and phase patterns, as well as helix wavefronts of the LG (l,p) beams with different topological charges and radial indices (denoted as p), respectively [26].
In the past two decades, the properties of vortex beams have been extensively studied, since the breakthrough made by Allen [21]. He pointed  out that the LG beams carry intrinsic OAMs quantified as lℏ along their axis, analogous to the spin angular momentum σℏ with σ = ±1/2 denoting the spin quantum number. Due to the transfer possibility of OAM from photons to particles, the vortex beams were often used as optical spanner to manipulate particles [24]. On the other hand, because these photon OAMs can form qubits in high-dimensional Hilbert space, vortex beam can be further used for information transmission [22]. Particular interests have been concentrated on the interaction of angular momenta, i.e. spinorbit coupling, which produced marvelous evolutions of light fields such as spin Hall effect of light [55][56][57][58]. These related performances of optical vortices have been introduced into electric beams, atom beams, and even matter wave systems [59][60][61]. The other type, which approximates 27% of total vortex lines, connects themself to form closed loops [52], as shown in Figure 1(b), with complex topologies as links or knots, as shown in Figure 3. These configurations occur in the superposition of waves more than three and represent the global topology of light fields [62]. Since the vortex line always occurs on the intersection of the wavefront surfaces with all phase labels, these topologically structured vortex lines indicate that all wavefront surfaces have a knotted boundary curve, and are thus multiply connected. Although the vortex topology occupies a finite volume of space, it affects the entire wavefield. Therefore, such nontrivial topologies of optical vortex widely arising elsewhere in physics, such as flows in fluid dynamics [42], quantum condensates [63], and field theory [34], have received considerable concern.

Optical vortex knots and links
Knots are defined mathematically as closed curves in the 3D free space. The first description of vortex knot can be tracked back to the hypothesis of Lord Kelvin that atoms were naturally knotted vortex loops in the ether. Although this hypothesis has been falsified, the remarkable topology theory has inspired the development in mathematical knot and many physical fields, where knotted vortices have long been thought to play an important role. In recent years, it has become experimentally accessible in disparate physical systems [31,42,48]. Therein, the topological skeleton unites those many seemingly unrelated physical phenomena. The construction of vortex lines and related topology modulation in these systems is an area of intense theoretical and experimental investigation.
As the initial conceive of introducing topology concept, knot has been studied as a fundamental conserved quantity in fluid mechanics for a long time, but this conservation characterization in the presence of dissipation has proven difficult to resolve [42]. Although the topology conservation in conventional fluids is unprocurable, due to the inherent nonlinearity of most dynamic fields, this challenging performance has been achieved in superfluid [64]. In contrast, universality and the stability of linear homogeneous media in optical system have evoked tremendous interests to exploit the optical counterparts.
As originally emphasized by Nye and Berry, optical vortices are naturally curved lines embedded in 3D wavefields [19]. Optical vortex lines are naturally found throughout optical speckle [52], occurring as dark points where the singularity lines intersect the imaged plane. In 3D optical speckle, these vortex lines form fractal tangles, percolating through space with many closed loops that are linked together. However, how to construct stable and isolated vortex knots in a limited space, especially in a propagating light beam, is a prerequisite to reveal the topology property and explore practical applications.

Design methods
Tying isolated optical vortex knots and links, the theoretical challenge lies in finding the topologically nontrivial solutions to complex scalar Helmholtz equation in all three dimensions. In 2001, Berry and Dennis first answered this question by theoretically constructing explicit and exact solutions of Helmholtz equation [51], in which the optical vortex lines can be knotted or linked, via the unfolding operation of a perturbation on unstable vortex loops. The wave function of this process can be described as Ψ(r) = A(r)exp (ilϕ) + Ψ perturb (r), where r and ϕ were the cylindrical coordinates. For the explicit construction, the superposition of Bessel beams, i.e. A(r)exp(ilϕ), which were exact solutions themselves, was chosen to implement a highorder dislocation loop, then a perturbated Gaussian background, i.e. Ψ perturb (r), was introduced to unfold this single vortex loop into knotted or linked vortex line. Figure 4 schematically shows the unfolding operation [51]. Subsequently, the LG set of optical modes that contain on-axial vortex lines with different waist widths, as a convenient basis for specifying modal superpositions, was utilized by Leach to experimentally realize the same topological configuration [31]. The other vortex topologies have also been generated in a specific superposition of beams, e.g. linked and knotted synthetic magnetic fields [65].
As shown in Figure 4(b), these previous vortex knots and links must be threaded by infinite vortex lines, which destroy the topology configuration as a single unit and underlying wavefront topology of interest. To construct isolated vortex knots and links, the Milnor map that describes complex scalar fields with knotted zero lines, as polynomial expansions around singularities in high-dimensional spaces, was introduced to mathematically realize these vortex knots and links [48,66]. As such, a periodic braid enclosed in a cylinder was smoothly mapped onto the solid torus under the Milnor mapping that the top and bottom of the cylinder bent around to meet each other, as schematically shown in Figures 5(a) and 5(d), so that the braided vortex line was spliced into an isolated loop residing on a torus, e.g. torus knot, with the height coordinate h replaced by the torus azimuth. This typical configuration is defined by two coprime integer numbers (n,m) in such a way the vortex line winds n times around a circle inside the torus and m times around a line through the hole in the torus [34]. Figures 5(b) and 5(e) plot the corresponding trefoil and figure-8 knots, whose phase distributions at the z = 0 plane are shown in Figures 5(c) and 5(f) [48].
In this kind of modal superpositions, the topological behavior of the vortex lines can be predicted with polynomial beams. Yet the linked vortex lines intertwine within regions of very low optical intensity; thus, the reconnection is extremely sensitive to perturbations in the constituent beams. In order to realize the polynomial beams physically, in practice, they are replaced by the realizable beams that satisfy the optical wave equation, e.g. such as LG beams, so that the vortex knot within superposition field has a finite scale. Moreover, these LG components should have a same Gaussian background [48]. Consequently, the z-dependent Gouy phase shift of these LG components induces variant vortex to construct knots in a confined space [67]. Therefore, the practical generation of these topological features with robust topology relies on the numerical optimization of the complex mode coefficients and the size of Gaussian background. Figure 6 shows some vortex knots constructed by polynomial beams with optimized complex coefficients [45].
In the linear homogeneous media, the diffraction and dispersion may diffuse the field structure and reduce the interference contrast. For instance, under the effect of a diffractive background, the vortex and knot corresponding to the destructive dark line are non-separable after a certain diffraction distance. To overcome this issue, Veretenov studied a new class of 3D topological dissipative solitons in homogeneous laser media with fast saturable absorption, namely, hula-hoop solitons [68], which combined the topological structures like vortex knots and the stability of solitons, supplying the potentials of optical information processing, including coding of topologically stable 3D symbols.
Besides, there also are exotic time-dependent solutions of Maxwell's equations, whose electric and magnetic field lines preserve topological configurations [34,41,69,70]. For example, Irvine and Bouwmeester studied the exact knotted solutions in free space based on a topological construction known as a Hopf fibration, in which all electric and magnetic field lines were closed loops and any two electric or magnetic field lines encode into knots and links [47]. Importantly, although these electric and magnetic fields expanded and deformed along with time, their topologies were preserved, as the helicity of the field configuration kept a conserved quantity.

SLM realizations
Owing to the requirement of exact control on the complex coefficients of vortex carrying and perturbed components, holography is the most convenient and accurate approach [71]. To implement the hologram design, a programmable spatial light modulator (SLM) is the common device [72]. Such device is composed of liquid crystal arrays, which introduce spatially dependent phase retardation or amplitude modulation across the beam profile, by uploading calculated gray patterns. Compared with a static hologram, the SLM allows controlling the beam parameters dynamically. However, limited by the response freedom of degree of wavefront parameter, e.g. most of the commercial SLM products are of phase-only type, it is necessary to design a phase-only hologram, which can control not only the phase structure of the diffracted beam but also its intensity [46].
In the experiment, the holograms uploaded into the SLM are calculated from the complex amplitude distributions of the knotted or linked fields in their waist plane. The SLM regenerates appropriate interference patterns, pixel by pixel, and the desired topologies will be reconstructed within light fields nearby the output plane of a 4 f system [46]. For example, Figures 7(a) and 7(b) illustrate the construction of a Hopf link and trefoil knot based on the LG polynomial beams with a finite width [73]. Suppose that these LG components are modulated by a perturbation term and Gaussian background, the knotted and linked fields at the z = 0 plane thus can be described as ψ link, knot = Σa i LG(l,p,w 0 ), with a i the coefficient of LG field, and w 0 the waist width that determines the 3D scale of isolated topology, and the corresponding amplitude and phase distributions are shown in  It is noteworthy that the size of the Gaussian background in polynomial beams significantly affects the topological configuration [48]. For the holographic method, the incident beam size can likewise affect the topological structure. Figure 8 displays the numerically configured Hopf links under the illuminations of Gaussian beams with different waists [73]. As shown in Figure 8(a), for a plane wave incidence, an isolated vortex link occurs because of the pre-modulated waist in the hologram. As the decrease of the waist radius, extra vortex lines appear and go close to the isolated vortex link. Once the incident beam waist is small enough, extra vortex lines reconnect with the isolated vortex lines, and then the isolated topological link is broken. This also proves that the holographic method guarantees sufficient generation robustness.
For the common generation approach based on holography, the reconstruction quality is closely dependent on the response of device pixel. For the SLM, relatively large pixel size induces the disadvantages such as misalignment and pixelation, which limit the construction perfection. Therefore, efficient optical device with smaller pixel is necessary to overcome such issues.

Metasurface realizations
Anisotropic metasurfaces, which work as nanoscale birefringent phaseshifters, in recent years have been substantially studied [74], showing application prospects in achromatic metalens [75], polarization-dependent holographic display [76], broadband response [77], and efficient optical devices [78] et al. This artificially planar material enables the simultaneous modulation of local phase, amplitude, and polarization of light field in subwavelength scale, molds optical wavefront into arbitrary mode, providing a new possibility for the miniaturization and integration of modern optical devices [79][80][81][82][83]. Typically, various metasurfaces have been proposed with advantages in applications associated with optical vortex, such as OAM holography [84], multiplexing, and demultiplexing [85] et al.
Compared with conventional technologies referring to SLMs, holography based on metasurfaces exhibits significant advantages in the control of polarization and large view angle [86], because of their inherent properties of ultra-small pixel, polarization response, and flexible phase gradient. More recently, some subtle metasurfaces have been proposed to the construction of topological optical fields with ultrasmall scale [87][88][89][90]. Typically, a metal type metasurface that consists of rectangular gold nanorods was proposed to realize the holographic reconstruction of vortex knots and links [87]. This  [88]., copyright 2020, WILEY-VCH.
reflective metal metasurface was designed based on the Pancharatnam-Berry (PB) phase, according to the pure phase-encoded holograms calculated from the complex amplitude of the Hopf linked or trefoil knotted light field at the z = 0 plane. The schematic of tying ultrasmall vortex knot and link is shown in Figure 9(a).
When a circularly polarized light was incident vertically, it produced topological knot or link that was six orders of magnitude smaller than the SLM method. Remarkably, due to the spin-orbit interaction, the angular position of the reconstructed topology was dependent on the incident circular polarization state [see Figure 9(c)]. Moreover, because of the PB phase modulation mechanism, this kind of metasurface had the advantage of broadband response at the NIR frequency with high efficiency, as shown in Figure 9(d).
A holographic all-dielectric metasurface based on the combination of PB phase and propagation phase was proposed, to compact and reconstruct multichannel topological configurations [88]. The knotted and linked fields at the z = 0 plane were decomposed into LG polynomial components, as a slice of a Milnor polynomial. Such two typical fields were encoded into respective holograms via the inverse sinc functional phaseonly encoding technique, and then these two holograms were combined onto the dielectric metasurface composed of rectangle pillars with variant  [100]., copyright 2019, The Royal Society of Chemistry; (e)-(f) ref [89]., copyright 2020, WILEY-VCH. geometry and orientation angle, based on the independent spin modulation mechanism.
As Figure 10(a) shows, for the incidence of a linearly polarized beam, polarization-switchable topological configuration with ultras-mall size, e.g. trefoil knot and Hopf link, were configured. Likewise, by using the phase singularity locating method, as to be described in Section 7, such small  [72]., Copyright 2016, nature publishing group; (c) -(h) ref [103]., copyright 2020, optical society of America. topologies were also accurately mapped and analyzed. An experimental map trefoil knot is shown in Figure 10(c). Importantly, although the topological configuration occupied a finite volume of space, as the isolated link shown in Figure 10(d), this topology affected the entire wavefield, even the far-field structure. Figures 10(e) and 10(f) compare the intensity and phase distributions of field beyond the isolated link and knot, which are obtained from the far-field detection and the Fourier transform of field at the z = 0 plane, respectively. These ultra-small topological fields and the vortex structures show robust configuration conservation property.
For the optical vortex knots and links generated from holographic metasurface, the configuration scale is determined by the metasurface dimension. In other words, the advantage of metasurface in dimension enables the possibility of further decreasing the magnitude of topological configuration. As the theory predicts, the topological structure can be modulated further to the scale of several wavelengths by adjusting the fabrication device, on the premise that the integrity of the phase profile can be reconstructed. For instance, the acoustic vortex knots were recently realized via holographic metamaterials with only 24 × 24 pixels [90]. The measured Hopf link and trefoil knot configurations have been reduced to about 4λ transversely, and about 8λ longitudinally, at wavelength scale.
The research on the polarization singularity has long history likewise the phase singularity, because of the crucial role of polarization response in lightmatter interaction. In the past decade, the vector beams embedded with onaxis polarization singularities [91], namely, polarization vortices, whose polarization azimuths are undefined, akin to the first type vortex beam, have attracted tremendous interests, because of their unique advantages for the construction of light fields that break the diffraction limitation and present delicate spatial and polarization structures [92][93][94][95][96]. Furthermore, some typical topological structures, e.g. polarization Möbius strip [97], have been constructed by the tightly focusing of vector fields. Recently, 3D polarization knots composed by vortex singularities of complex Stokes field s 1 + is 2 , such as the torus knots and links as well as the complicated figure-8 knot, have been generated based on the superposition of a traditional knotted vortex beam and a Gaussian beam with opposite spin states [98,99].
As the advantages of more degrees of freedom and miniaturization, the knot and link with controllable topological configurations formed by the polarization singularity have been achieved with the help of metasurfaces. Figure 11(a) shows an ultrathin all-dielectric metasurface that can achieve efficient generation and transformation of two arbitrary 3D polarization topologies [100]. This metasurface consists of two distinct elements, which are designed for the generation of tightly focused Poincaré beams composed of Gauss and LG components with opposite spin states. For the illuminations of the LCP or RCP light beams, it produces a Möbius strip with three or four half-twists. Similarly, a composite metasurface that enables complete control of arbitrary spatially varying polarization fields has been proposed to generate polarization knots [89]. This metasurface was designed according to a coherent pixel polarization metahologram, which can be used to generate angle-multiplexed output fields with orthogonal polarization pairs to construct polarization knots. As shown in Figure 11(f), three azimuth-controllable polarization knot profiles where the C-points of the polarization singularities are unchanged were experimentally realized.

Characterizing topological configurations
Several methods have been proposed to observe the vortex topologies, of which the key point is how to locate the phase singularity accurately. The most widely used method is detecting the locations of dark cores caused by singularities [31,46], that is, intensity locating method. For this kind of method, consider that the fine details of the vortex lines lie in regions of near darkness and in order to measure the vortex positions within the beam cross-section directly, it is necessary to over-saturate the record CCD. Figure 12(a) illustrates the experimental results of the intensity distributions of a link threaded by two vortex lines at the z = 0 plane, where the inset shows the overexposure intensity pattern. As shown, the resulted image is mainly white, but contains points of darkness corresponding to the vortices as they intersected the plane of the CCD. After the z-scanning detection, the combination of transverse vortex coordinates measured in neighbouring planes enables the configuration of 3D topological structures. This overexposure intensity method has a rapid speed to locate the singularities. However, the locating accuracy is unsatisfactory, in an ultra-small topological light field, because the dark spots are not so small enough when the background field is weak, and then several singularities that close together overlap, as the results shown in Figure 12(h).
One improved method that locates the positions of phase singularities from the forked stripe within interferogram produced by the vortex field and a reference planar wave was proposed [72], as shown in Figure 12(b). The locating precision of this method is mainly limited by the period of interference fringes. Further, a method based on phase pattern has been proposed, where the singularities are located from considering the phase change around a small loop of pixels [101]. Likewise, this method is inaccurate for the case of multiple singularities coexisting in an ultra-small topological light field. To measure the phase distribution of the entire optical field, the phase shift method is necessary [102], which requires a multi-step measurement of each plane, greatly decreasing the measuring speed and increasing the reconstruction error.
Recently, a more accurate and rapid method to experimentally map the topology configuration of optical vortex has been reported [103], in which the phase distribution in any transverse plane was firstly measured by the digital holography, and the location of phase singularity was accurately obtained by utilizing a numerical search algorithm based on local spatial frequency. Meanwhile, self-calibration property between multi-plane singularities was adopted for 3D reconstruction of vortex lines. Figures 12(d-f) shows the locating process of phase singularities based on the digital holography. In contrast to the overexposure intensity method, this method has higher precision in locating the inflexion of vortex lines, as the experiment results shown in Figures 12(g-i).

Applications
The light embedded with optical vortex has OAM and associated transverse energy flow and thus can exert a torque, which has been widely exploited in particle manipulations as light spanners [104,105]. Knotted optical fields supply the possibility of extended and knotted optical manipulations in three dimensions. This hypothesis was recently realized by light fields with knotted trajectories of bright spots [106], whose intensities are complementary to the knotted vortex lines, with multi-particle transport in knot circuits and networks under programmable paths [107].
Knot theory has already influenced particular areas of physics, e.g. topological quantum computing, where knot can be used as 'quantum core'. Entanglement of linked optical vortex loops in the light, within macroscopic and finite volumes, was realized by spontaneous parametric downconversion [108]. This entanglement of photons in complex 3D topological states suggests the possibility of entanglement of similar features in other systems. In addition, the knotted fields can act as synthetic gauge fields on neutral atomic gas, and transfer the topological properties to the ground state of a 3D Bose-Einstein condensate [65].

Conclusion and outlook
With the deepening of research on light field modulation, novel light fields with distinct configurations have been continuously developed, which sufficiently support for expanding the application of structured light fields. As above mentioned, in order to make the polynomial beams physically realizable, the polynomial is multiplied by a Gaussian amplitude function with width. Therefore, the size of the topology is order of the waist radius transversely (10 3 -10 4 λ), but of the Rayleigh length longitudinally (10 6 -10 8 λ) [109]. This extreme aspect ratio is impractical to imprint into material systems, such as ferroelectrics [110] and chiral nematic colloids [111]. Fortunately, with the development of sub-wavelength and near-field wavefront control technologies such as metasurfaces, and nonparaxial modulation methods such as tightly focusing, the configuration dimension can be further decreased.
Once the topological configuration is compressed to the wavelength scale and even smaller, the imprint possibility between it and the material structure will greatly promote the development of light-matter interaction. Remarkably, this ultra-small-structured light field provides the possibility to construct knotted structures in photosensitive materials, and even to control the velocity of a vortex knot excited in the superfluid wave function of a Bose-Einstein condensate [63].
As topologically stable objects in field theories, knots have been put forward to explain various persistent phenomena in systems ranging from atoms and molecules to cosmic textures in the universe. The topological manipulation observed in other wave system [112], such as the tying and untying knots that occur in hydrodynamics [113], can be realized in optical system, utilizing the holographic reconstruction method based on polynomial beams. Topology switch and even the continuous transformation will bring new opportunities for topological information processing.

Disclosure statement
No potential conflict of interest was reported by the authors.