1,999
Views
0
CrossRef citations to date
0
Altmetric
Reviews

Optical simulation of various phenomena in curved space on photonic chips

, &
Article: 2153626 | Received 29 Mar 2022, Accepted 16 Nov 2022, Published online: 07 Dec 2022

ABSTRACT

Transformation optics have been an essential paradigm to manipulate electromagnetic waves on the subwavelength scale and have brought various functional photonic architectures into integrated photonic chips. On the other hand, in the spirit of analogical thinking, classical and quantum simulations of general relativity have been extensively studied in diverse physical systems. In this review, we summarize recent advances in analogical gravitation based on integrated photonic chips with the aid of transformation optics. Meanwhile, different types of transformation optical structures, such as gradient waveguides, metasurface waveguides, waveguides on curved space and gradient waveguide arrays, emulating a variety of phenomena in curved space are reviewed, including the gravitational lensing of black holes, Einstein rings, cosmic strings, the particle pair evolution near the event horizon and so on. Furthermore, perspectives for the study of analogical gravitation based on integrated photonic chips are discussed.

Graphical Abstract

1. Introduction

Analogy has been widely used in scientific research and dates back to the ancient Greek philosopher Aristotle, who was the first to lay the foundations of the theory of analogical reasoning. Analogical study, which allows scientists to go beyond their own professional scope to conduct interdisciplinary research, has been assumed to be one of the most vital and efficient tools to explore and understand nature. For instance, De Broglie proposed the genuine hypothesis on wave properties of material particles inspired by the analogy between geometric optics and Newtonian particle dynamics, expressed by the formal similarity between variational Fermat’s principle of ray optics and the least action principle of Newtonian mechanics. Another celebrated example is that Ernest Rutherford envisioned the atom in the spirit of the solar system, with electrons orbiting around a massive nucleus that occupy only a very small part of the atom. Moreover, it is worth mentioning that another notable analog case is topological physics [Citation1,Citation2], a current hot topic research, which originates from the concept of the topology in mathematics that deals with highly robust conserved quantities that do not change when physical objects are continuously deformed. Furthermore, topological phases, originally discovered in solid-state electron systems, have been expanded to other subfields of physics, such as optics and photonics [Citation3,Citation4], acoustic and mechanical systems [Citation5,Citation6].

In the spirit of analogical thinking, there has also been a fruitful intertwining between different branches of classical physical systems. Taking the early development of electrodynamics as an example, electric fields, magnetic fields and their laws are analogous to flow velocity fields and their laws, such as the introduction of concepts of vortex, source, sink, etc. Such analogues are essential to advance the development of electromagnetism. Additionally, another intriguing case is laboratory analogues of general relativity (GR). In the early 1920s, Gordon noticed that moving isotropic media appear to electromagnetic fields as certain effective space-time geometries [Citation7]. In 1960, Plebanski formulated the analogy [Citation8] between the macroscopic Maxwell’s equations in complex inhomogeneous media and free-space Maxwell’s equations for the background of an arbitrary spacetime metric, and especially proposed an original viewpoint that smooth deformation of space can be physically equivalent to inhomogeneous and anisotropic electromagnetic medium. Unfortunately, there are physical limitations to obtain such specially designed electromagnetic parameters for conventional materials made of atoms. In 2006, Pendry and Leondart independently and originally put forward transformation optics [Citation9,Citation10], which offers unparalleled opportunities for controlling light propagation through careful refractive index engineering using metamaterials. In particular, metamaterials [Citation11], artificially engineered materials on the subwavelength scale, can be designed to induce customized properties with arbitrarily effective parameters that originally did not exist, such as the negative refractive material [Citation12] and perfect imaging [Citation13–15]. One of the amazing applications of transformation optics is invisibility cloaks [Citation16–22], in which light is regarded as linear parallel geodesic rays in deformed spaces. Moreover, with the aid of transformation optics and metamaterials, various well-known phenomena related to GR have been emulated, ranging from artificial black holes [Citation23–26], Minkowski spacetimes [Citation27], electromagnetic wormholes [Citation28], and De Sitter space [Citation29] to cosmological inflation and redshift [Citation30,Citation31].

On the other hand, there has also been increasing interest in the analogy between classical physics and quantum physics. For example, based on the similarity between the paraxial optical wave equation describing the spatial propagation of a monochromatic light beam and the temporal Schrodinger equation for a quantum particle, a wide variety of coherent quantum effects encountered in atomic, molecular or condensed-matter physics have been mapped into optical waveguides. Given that understanding and accessing coherent wave phenomena in microscopic quantum systems is often a challenging task because of complications arising from many-body effects, decoherence, and the presence of time-dependent or nonlinear terms in the Schrodinger equation, such analogies in optical waveguides have the possibility of exploring coherent dynamical regimes not yet accessible in original quantum systems. For example, photonic lattices made of single-mode coupled waveguides have become a popular tool for simulating a plethora of quantum phenomena in condensed-matter systems, ranging from photonic Floquet topological insulators [Citation32] and Weyl physics [Citation33] using spiral waveguides, Landau levels using strain photonic lattices [Citation34], Aharonov-Bohm cages [Citation35] to the demonstration of flat-band states and dispersion-less light propagation [Citation36]. At the same time, such quantum-classical analogies in waveguides in turn promote the acquisition of novel integrated optical functions and devices, such as fractional discrete Fourier transforms [Citation37] and quantum logic gates [Citation38] on such photonic chips. It is emphasized that such quantum-optical analogies using optical waveguides have been well summarized, and several pertinent reviews have been given [Citation39,Citation40]. In addition to quantum-optical analogies, other quantum-classical analog systems have also been demonstrated, for example, exploiting water surface waves to emulate Hawking radiation [Citation41] and Penrose superradiance [Citation42], employing electric circuit networks [Citation43] to study higher-order skin-topological effects in condensed matter.

At the same time, the over past decades have witnessed remarkable progress in quantum simulations, Feynman’s innovative idea [Citation44] that one could employ a controllable quantum system to study another less controllable or inaccessible quantum system. A variety of quantum architectures [Citation45], such as atoms and ions, nuclear and electronic spins, superconducting circuits, and photons, with salient advantages of a high degree of controllability, novel detection possibilities and extreme physical parameter regimes, have been demonstrated as flexible quantum simulators for disposing of intractable problems in, e.g. condensed-matter physics [Citation46], high-energy physics [Citation47], quantum chemistry [Citation48], and cosmology [Citation49]. For instance, Hawking radiation [Citation50], an undetectable quantum effect near the event horizon of black holes owing to the extremely low radiation temperature and lack of astronomical observation tools, has been successfully emulated in various quantum simulators, ranging from ultracold atoms [Citation51], Bose–Einstein condensates [Citation52], trapped ions [Citation53], and Fermi-degenerate liquids [Citation54] to superconducting circuits [Citation55] and nonlinear optical media [Citation56,Citation57]. Additionally, the complex interacting many-body problems in natural condensed-matter physics, such as the Hubbard model [Citation58,Citation59], spin glasses [Citation60], quantum phase transitions [Citation61] and so on, have also been studied using quantum simulations with high efficiency. Encouragingly, what is more amazing is the recent advances of special quantum simulators for addressing certain mathematical tasks, such as the Sycamore quantum processor [Citation62], Gaussian boson sampling machine Jiuzhang [Citation63] and Zuchongzhi quantum processor [Citation64], which have fully exhibited quantum advantages to perform well-defined intractable problems compared to classical computers.

As mentioned above, the classic and quantum simulations have cross-penetrated in different branches of physics and even various subjects of science, and have gained numerous and fruitful achievements with well-demonstrated vigorous vitality. For the sake of brevity, the aim of this view, hence, is not intended to cover all the related advances but to focus on the current progress on optical simulations of various classical and quantum phenomena about GR emulated on photonic chips using photons. Stirred by metamaterials and transformation optics, we introduce several concrete optical architectures on photonic chips to emulate some interesting classic and quantum phenomena related to GR, aim for a comprehensive treatment that incorporates state-of-the-art experiments, and especially point the way toward potential applications. This review is organized as follows: In Section 2 we briefly review the underlying principle of optical simulation of GR. In Section 3, we briefly classify optical waveguides on chips, which are judiciously designed and leveraged to carry on the research of optical simulation of GR. In Section 4, the optical simulation of classical GR phenomena using various types of photonic chips are reviewed, such as gradient waveguides, metasurface waveguides and waveguides on curved space. In Section 5, optical simulations of phenomena about quantum fields in curved space used gradient waveguide arrays are highlighted. In Section 6, we summarize the outlook and future of optical simulation and quantum simulation on photonic chips.

2. The underlying principle of the optical simulation of GR

First, we briefly review the principle underlying the GR-optical analogies, which is based on the form invariance of Maxwell’s equations between the complex inhomogeneous media and the background of an arbitrary spacetime metric. It is emphasized that the specific description of the invariance has been well summarized in [Citation8,Citation25,Citation65]. In general, a geometry of the space-time is described by the metric , where we take with

Greek indices running from 0 to 3 to denote the space-time coordinates. Latin indices indicate the spatial coordinates and run from 1 to 3, whereas describes time measured in spatial units with being the speed of light in vacuum. The metric tensor may vary as a function of the coordinates, because space-time may be curved or because curved coordinates are used in flat space. The determinant of measures the four-dimensional volume of an infinitesimal space-time element. According to customary in general relativity, we denote the matrix inverse of by , where the position of the indices indicates that is co-variant and is contra-variant under coordinate transformations. We give the free-space Maxwell equations in Cartesian components:

(1)
(2)

where is the completely antisymmetric Levi-Civita tensor. The spatial indices indicate that in this representation and form the components of vectors that are co-variant under purely spatial transformations, whereas and constitute contra-variant spatial vectors. The charge density and the current density are given by and of the four-current . In empty but possibly curved space, the electromagnetic fields are connected by the constitutive equations in SI units.

(3)

with the symmetric matrices and and the vector given as

(4)

We can see that empty but curved space-time manifolds or curved coordinates appear as effective dielectric media with anisotropic parameters and , especially with a magneto-electric coupling vector . Since the material representation Equationequation (4) requires both electric and magnetic response and the magneto-electric coupling parameter, it is a huge challenge to realize these designed material parameters when the corresponding relation was firstly proposed. Fortunately, with the advent of metamaterials, the material parameters corresponding to some metrics of space-time have been realized.

Among these complex metrics of space-time, an isotropic central symmetric diagonal metrics has attracted attention thanking to easy experimental realization. For a class of centrally symmetric metrics that can be written under coordinate transform in an isotropic form where , and . Note that both coordinate systems represent the same space-time. This equivalence allows the rendering of the material response to be reduced from anisotropic Equationequation (4) to isotropic with an effective refraction index:

(5)

that could be achieved with pure dielectric non-dissipative and non-dispersive materials.

3. The classification of optical waveguide on chips for optical simulations

As is well known, photons, as information carriers, play an important role in the current information society. High efficiency manipulation of electromagnetic waves or photons on chips has always been a pursued goal in the optical community. In particular, the high development of the current information society requires increasingly advanced photonic integration technology [Citation66]. Fortunately, a range of integrated photonics platforms, e.g. silica-on-insulator [Citation67], laser-writing silica [Citation68], silicon nitride [Citation69], lithium niobite-on-insulator [Citation70], gallium arsenide [Citation71], indium phosphide [Citation72] and many others, have exhibited a strong capacity for the manipulation of photons with high efficiency to obtain the desired effective electromagnetic parameters.

In general, optical waveguides are fundamental elements in integrated optics. Essentially, all integrated devices are built on waveguides. To some extent, the properties of optical waveguides determine the performances of integrated components on the chip. Structural design and refractive engineering for optical waveguides, thus, are important for the realization of novel and high-performance devices. On the other hand, recent years have seen a boom in metamaterial research. Hence, there is also surging motivation to bring metamaterials into optical waveguide technologies and on-chip architectures. This method has new degrees of freedom to control the flow of light in integrated photonic devices. In particular, with the aid of ever-improving nanofabrication techniques, such as electron beam lithography, focused ion beam, and deep ultraviolet lithography, subwavelength metamaterial structures can be fabricated in the same lithography step as conventional waveguides by using manufacturing processes that are well established in the semiconductor electronics industry, thus making their integration straightforward. A plethora of advanced integrated devices on chips with unprecedented performance have been well demonstrated, e.g. polarization beam splitters [Citation73] and rotators [Citation74,Citation75], waveguide crossings and bends [Citation76,Citation77], mode converts [Citation78] and multiplexers [Citation79,Citation80]. At the same time, the emerging concept of transformation optics, a design method for controlling photons at our discretion, opens up unique possibilities for advancing the integration of complex functionalities in photonic circuits, such as a flattened Luneburg lens [Citation81], Maxwell’s fish-eye lens [Citation82] and dual-function ‘Janus’ devices [Citation83].

Optical waveguides can be classified according to their geometry, mode structure, refractive index distribution, and material (see ). In general, optical waveguides comprise a longitudinally extended high-index medium called the Core, which are transversely surrounded by a low-index medium called the Cladding. In the following, invoked by the design concept of transformation optics and the mature development of fabrication at the micro-nanoscale, we can judiciously manipulate either the core layer or cladding layer or both in optical waveguides to realize distinct optical functions. Intriguingly, with the aid of nontrivial optical waveguides, we can investigate the large-scale gravitational astronomical phenomena predicted by GR, which are either inaccessible or cannot be studied repeatedly and carefully. On the other hand, it is noted that such devices inspired by GR can have unique advantages in some aspects compared to conventional optical devices.

Figure 1. Different type waveguides for optical simulations inspired by general relativity.

Figure 1. Different type waveguides for optical simulations inspired by general relativity.

4. Optical simulation of classical GR phenomena

Gradient waveguides

Recently, gradient waveguides whose core thickness varies at our discretion have attracted considerable attention in the manipulation of photons on integrated chips. For instance, Thomas Zentgraf et al [Citation81]. explored greyscale lithography to adiabatically tailor the topology of a dielectric layer adjacent to a metal surface to demonstrate a plasmonic Luneburg lens and Eaton lens that can, respectively, focus and bend surface plasmon polaritons (SPPs). Shuyi Li et al [Citation82]. also used greyscale electron-beam lithography to adiabatically change the profile of a Si layer on a commercial silicon-on-insulator wafer to achieve a Maxwell fisheye lens for multimode waveguide crossing, which has manifested an advantage in on-chip multimode routing and communications with a highly integrated and large capacity. Hyuck Choo et al [Citation84]. experimentally demonstrated the achievement of highly efficient nanofocusing of electromagnetic waves into the deep subwavelength regime by leveraging a three-dimensional Au–SiO2–Au taper with a carefully engineered SiO2 layer. Xiangyang Wang et al [Citation85]. exploited a gradient polymer waveguide to achieve the Mikaelian lens to simultaneously obtain a self-focusing effect in the geometric-optics limit and a Talbot effect in the wave-optics limit. In particular, the Talbot effect constructed in such a lens can be used to transfer a digital field pattern without diffraction and has potential applications in digital coding communications. Additionally, optical cloaking [Citation86] and optical super-imaging [Citation87] have also been achieved with the aid of gradient waveguides.

Interestingly, gradient waveguides with high efficiency to model light can be utilized to mimic large-scale astrophysical phenomena in curved spacetime, such as celestial mechanics, Einstein’s field equations, gravitational lensing and the event horizons of black holes. Chong Sheng et al [Citation88]. judiciously embedded a microsphere into a planar polymer waveguide formed during a controlled spin-coating process (see (b)). Due to surface tension effects, the polymer waveguide layer around the microsphere is distorted, resulting in a continuous change in the waveguide effective refractive index that, under certain conditions, can mimic the curved spacetimes caused by strong gravitational fields. In such a gradient waveguide, the lensing and asymptotic capture of the incident light in an unstable circular orbit that corresponds to the photon sphere of a compact stellar object are observed (see (c, d)). Following the same idea, the gradient waveguide can also demonstrate Einstein’s Rings [Citation89], which is a famous phenomenon predicted by GR and observed in astronomy. The Einstein ring phenomenon occurs when light from a point source is deformed by a mass distribution through gravitation lensing that causes the appearance of a ring around the mass distribution. In such a setting, the spherical wavefront produced by the arc grating emulates the wave radiated outwards from a point source. When this 1D spherical wavefront passes by the star, it is focused, and the beam width changes as a function of the propagation distance, as extracted from the experimental data. Typical results for two different ‘stars’ (microdroplets with two different radii) are fabricated, as shown in (g, h). As the radius of the ‘star’ is larger, the convergence of the beam is more extreme, but the final beam is wider. At this point, it is intriguing to compare emulation results with Einstein’s prediction. The Einstein formula for the angular diameter of the virtual ring depends on the convergence angle , the radius of the mass distribution and the distance between the center of the mass distribution to the observation point z. To conform with the Einstein formula, the relative angular radius between the two mass distributions (two samples) was calculated and agreed well with theory (see (i)). Aside from gravitational lensing, Igor I. Smolyaninov theoretically proposed to utilize a gradient waveguide for emulating the Unruh effect [Citation90,Citation91], that is, an accelerating observer will treat its surrounding environment as a thermal radiation bath (see (j)). In this theory, when photons launched into a specially designed optical gradient waveguide will act as massive quasi-particles and experience strong acceleration, such a waveguide can be used as a thermometer that would measure the Unruh temperature. Subsequently, Hui Ge et al. exploited a gradient optical fiber [Citation92] to validate this Unruh effect using a leaky phenomenon within the optical fibers akin to bremsstrahlung (see (l)). Additionally, the emulated curved space in gradient waveguides can also be utilized as a tool of wavefront shaping in waveguide settings. For instance, the transformation of a broad Gaussian beam (FWHM ) into a very narrow nondiffracting beam (FWHM ) with the propagation of was achieved in a predesigned gradient waveguide by drawing on intuition from GR [Citation89]. A gradient waveguide whose effective index satisfies a gravitational field profile in the form of a linear potential could result in the oscillation between an incident broad Gaussian beam and an accelerating Airy beam [Citation93].

Figure 2. Gradient waveguides for emulating various general relativity phenomena. (a)-(d) Gravitational lensing of artificial black holes [Citation88]. (a) Depiction of light deflection by the gravitational field of a massive stellar object. (b) Schematic view of the microstructured optical waveguide formed around a microsphere and used to emulate the deflection of light by a gravitational field. (c, d) Scattered field intensity observed in the experiment and calculated using a full-wave finite-difference frequency-domain (FDFD) electromagnetic code. (e)-(i) Einstein rings [Citation89]. (e) Einstein’s vision: light from a point source is focused by a gravitational lens, and is subsequently observed as a virtual ring around the mass distribution. (f) Schematic view of the fabricated inhomogeneous waveguide. (g, h) Experimental results showing a spherical wave passing through two domes with different radii. (i) Fit to Einstein’s formula for the angular diameter of the Einstein rings. (j) Unruh effect in tapered waveguides [Citation91]. (j) Massive body thrown up in a gravitational field that can be emulated by waveguide tapering. (k)-(l) The bremsstrahlung in a tapered optical fiber [Citation92]. (k) A tapered fiber. (l) The radiation phenomenon generated by a tapered fiber.

Figure 2. Gradient waveguides for emulating various general relativity phenomena. (a)-(d) Gravitational lensing of artificial black holes [Citation88]. (a) Depiction of light deflection by the gravitational field of a massive stellar object. (b) Schematic view of the microstructured optical waveguide formed around a microsphere and used to emulate the deflection of light by a gravitational field. (c, d) Scattered field intensity observed in the experiment and calculated using a full-wave finite-difference frequency-domain (FDFD) electromagnetic code. (e)-(i) Einstein rings [Citation89]. (e) Einstein’s vision: light from a point source is focused by a gravitational lens, and is subsequently observed as a virtual ring around the mass distribution. (f) Schematic view of the fabricated inhomogeneous waveguide. (g, h) Experimental results showing a spherical wave passing through two domes with different radii. (i) Fit to Einstein’s formula for the angular diameter of the Einstein rings. (j) Unruh effect in tapered waveguides [Citation91]. (j) Massive body thrown up in a gravitational field that can be emulated by waveguide tapering. (k)-(l) The bremsstrahlung in a tapered optical fiber [Citation92]. (k) A tapered fiber. (l) The radiation phenomenon generated by a tapered fiber.

Metasurfae waveguides

Metasurfaces, two-dimensional metamaterials, have been important tools to control electromagnetic waves in both far-field and near-field regimes. Compared with a three-dimensional bulky metamaterial, this artificial plane material with subwavelength thickness greatly reduces fabrication time and mitigates fabrication complexity. Additionally, traditional metamaterials usually control the wavefront of an electromagnetic wave by accumulating the phase by propagating at a distance far larger than the wavelength. However, a metasurface can efficiently manipulate the wavefront of an incident electromagnetic wave through just the subwavelength propagation distance. Therefore, this can largely alleviate the propagation loss. A more detailed discussion of remarkable progress in metasurfaces for controlling electromagnetic waves on integrated chips have been well summarized in reference [Citation94,Citation95].

Interestingly, based on the high efficiency of controlling light by using metasurfaces, a metasurface waveguide can be utilized to conduct research on analogical gravitation. Theorists predicted that some topological defects may have formed during a symmetry breaking of the Higgs vacuum field in the early universe when the topology of the vacuum manifold associated with this symmetry breaking was not simply connected. Topological defects of spacetime include 0-dimensional monopoles, 1-dimensional cosmic strings and 2-dimensional domain walls. Fortunately, stimulated by the discovery of gravitational waves [Citation96,Citation97], a stochastic gravitational wave background possibly emitted by a cosmic-string network [Citation98,Citation99] and cosmic domain walls [Citation100] generated in the early Universe has been reported. Aside from gravitational waves, many other interesting nontrivial phenomena about these cosmic topological defects [Citation101–108] also have been predicted. Interestingly, C. Sheng et al [Citation109]. experimentally demonstrated definite photon deflections caused by emulating cosmic strings with the aid of metasurface waveguides, whose structure combines a slab waveguide with a rotational metasurface (see )). Because cosmic strings can change the topological structure of surrounding space-time to result in the defect or surplus of space-time angle, when light propagates in such a topological space-time, it will produce a definite deflection angle ( is the gravitational constant, is the linear mass density of the cosmic string), which is only determined by mass density parameter that originates from the intrinsic topological properties of the cosmic string, and which is regardless of the incident position, propagation direction, wavelength and polarization direction of photons. If , the string has a positive mass density and carries positive curvature at the origin, and the photon undergoes an attractive force and moves toward the origin of the string. When , the string has a negative mass density and carries negative curvature at the origin, and the photon undergoes a repulsive force and moves away from the origin of the string. By judiciously controlling the parameters of rotational metasurfaces, robust definite photon deflection in both positive and negative topological space has been observed as shown in ). This remarkably robust photon deflection can be used as a new kind of omnidirectional lens, which can bend light without changing the beam profiles. Such a unique property can be used in optical focusing, imaging, and information transfer.

Figure 3. Metasurface waveguides for emulating various general relativity phenomena. (a)-(h) Definite photon deflections by cosmic stings [Citation109]. (a, e) The schematic of cosmic string with negative and positive mass density using metasurface waveguides. (b)-(d) Photon deflection observed in the experiment under different impact parameters with negative mass density. (f)-(h) Photon deflection observed in the experiment under different impact parameters with positive mass density. (i, j) Cherenkov surface plasmon wakes [Citation110]. (i) Scanning electron micrograph of the nanostructures used to generate the SPP wakes. (j) Near-field scanning optical microscopy images and the calculated results. (k-m) The Bremsstrahlung radiation of moving particles using metsurfaces [Citation111]. (k) Geometric picture describing the same event in different coordinates through mimicking Bremsstrahlung radiation of moving particles (blue ball and purple ball) in flat space-time and curvilinear space-time respectively. (l1, m1) SEM images of different samples. (l2, l3) The simulations and experiments for the case of l1. (m2, m3) The simulations and experiments for the case of m1.

Figure 3. Metasurface waveguides for emulating various general relativity phenomena. (a)-(h) Definite photon deflections by cosmic stings [Citation109]. (a, e) The schematic of cosmic string with negative and positive mass density using metasurface waveguides. (b)-(d) Photon deflection observed in the experiment under different impact parameters with negative mass density. (f)-(h) Photon deflection observed in the experiment under different impact parameters with positive mass density. (i, j) Cherenkov surface plasmon wakes [Citation110]. (i) Scanning electron micrograph of the nanostructures used to generate the SPP wakes. (j) Near-field scanning optical microscopy images and the calculated results. (k-m) The Bremsstrahlung radiation of moving particles using metsurfaces [Citation111]. (k) Geometric picture describing the same event in different coordinates through mimicking Bremsstrahlung radiation of moving particles (blue ball and purple ball) in flat space-time and curvilinear space-time respectively. (l1, m1) SEM images of different samples. (l2, l3) The simulations and experiments for the case of l1. (m2, m3) The simulations and experiments for the case of m1.

In addition to cosmic strings, the accelerating particles in general relativity can be well-emulated by leveraging a one-dimensional metasurface that generates a predefined caustic of SPPs on a metal surface. Genevet et al [Citation110]. utilized SPPs generated on a metasurface through spin-orbit interaction as an analogous system to study Cherenkov radiation from a charged particle moving with a constant high speed. F. Zhong et al [Citation111]. also explored the geometric phase of SPPs caused by the spin-orbit interaction of photons to emulate accelerating frames in general relativity and obtained a new kind of accelerating beam, the Rindler beam, based on the Rindler metric in gravity. Moreover, compared to metamaterials with the limited scope of the effective index, the geometric phase of SPPs with the aid of metasurfaces provides a freedom to control the effective index that satisfies the requirement of the large effective index for emulating curved space.

Waveguides on curved space

In contrast to engineering space-dependent electric and magnetic parameters using metamaterials or metasurface to observe curved-space phenomena in the laboratory, another route for such studies is to create curved space by engineering the geometry of the space itself. This idea [Citation112], suggested in 1981, started by exploring the dynamics of a free quantum particle rigidly bounded to a curved surface, resulting in the Schrodinger equation containing a linear potential that is a function of the curvature.

Likewise, these ideas were brought to electromagnetic waves and first demonstrated by Batz and Peschel [Citation113], who explored the dynamics of light propagating within a thin-film waveguide covering the curved surface area of a three-dimensional body, effectively creating 2D curved space for the light. These works [Citation113,Citation114] clearly show that when a Gaussian beam propagates on a surface with constant positive curvature, periodic refocusing, self-imaging, and diffractionless propagation can be observed as shown in ). In contrast, light spreads exponentially on surfaces with constant negative Gaussian curvature (see )). Furthermore, the evolution of speckle patterns propagating along two-dimensional surfaces of constant positive and negative Gaussian curvature are utilized to study the influence of intrinsic curvature on HBT measurements [Citation115]. In addition to Gaussian beams, shape-preserving accelerating beams propagating on a curved space [Citation116,Citation117] with a constant positive curvature have been studies, which found that unlike accelerating beams in flat space, these wave packets change their acceleration trajectory due to the interplay between interference effects and the space curvature, and focus and defocus periodically due to the spatial curvature of the medium in which they propagate (see )). Additionally, C.N. Xu et al [Citation118]. theoretically studied the Gouy phase of light beams in both the paraxial and nonparaxial regimes on a two-dimensional curved surface by generalizing the angular spectrum method and found that there exists an extra phase shift introduced by the curvature of the surface. G.H. Liang et.al devised an air helicoid waveguide [Citation119] which was obtained with an air layer sandwiched by two metal helicoid surfaces to mimic an expanding universe. W.F. Ding et al [Citation120]. utilized a curved space to study laser propagation in a Rindler accelerated reference frame based on matrix optics.

Figure 4. Waveguides on curved space for emulating various general relativity phenomena. (a, b) An initial Gaussian beam on a positively and negatively curved surface [Citation114]. (a) Experimental realizations of light propagation on a sphere with a self-focusing. (b) Experimental realizations of light propagation on a hyperbolic with exponentially spreading. (c) Accelerating Wave Packets in Curved Space [Citation117]. (d) The speckle evolution along the hourglass [Citation115]. (e, f) Control of light by curved space in nanophotonic structures. (e) Schematic of the coupling scheme of the light to the paraboloid waveguide [Citation121]. (f) Curvature effects on diffraction. (g-n) Simulations of wormhole using curved optical spaces [Citation122]. (g) The embedding diagram of the 2D wormhole. (h) The image of the sample fabricated with 3D printing technique. (k-n) The calculated light geodesic lines with four impact parameters for the case that the light enter the wormhole through the on plane and exit through the other plane.

Figure 4. Waveguides on curved space for emulating various general relativity phenomena. (a, b) An initial Gaussian beam on a positively and negatively curved surface [Citation114]. (a) Experimental realizations of light propagation on a sphere with a self-focusing. (b) Experimental realizations of light propagation on a hyperbolic with exponentially spreading. (c) Accelerating Wave Packets in Curved Space [Citation117]. (d) The speckle evolution along the hourglass [Citation115]. (e, f) Control of light by curved space in nanophotonic structures. (e) Schematic of the coupling scheme of the light to the paraboloid waveguide [Citation121]. (f) Curvature effects on diffraction. (g-n) Simulations of wormhole using curved optical spaces [Citation122]. (g) The embedding diagram of the 2D wormhole. (h) The image of the sample fabricated with 3D printing technique. (k-n) The calculated light geodesic lines with four impact parameters for the case that the light enter the wormhole through the on plane and exit through the other plane.

More interestingly, such waveguides on curved space have been extensively utilized to emulate some general relativity phenomena about black holes and wormholes. R. Bekenstein et al [Citation121]. investigated the evolution of light in a paraboloid structure inspired by the Schwarzschild metric describing the space surrounding a massive black hole. In contrast to 2D fabrication techniques imported from microelectronics, this work employs the multiphoton polymerization technique to construct a special 3D waveguide whose size is comparable to the wavelength scale ()). Such a curved nanophotonic structure can not only control the diffraction properties and the phase and group velocities of wavepackets propagating within it but also exhibit tunneling through an electromagnetic bottleneck (see )). In addition to the Schwarzschild metric of black holes, R. Q. He et al [Citation122]. also employed a waveguide on curved space inspired by the Morris-Thorne traversable wormhole metric to show that a giant tidal force induces the divergence of optical deflection at the throat of the wormhole (see )). In this work, they gradually changed the impact parameter compared to the radius of the wormhole’s radius and compared three different situations. (1) When , the calculated light rays do not pass through the wormhole. In the process, the deflection angle increases with the decrease in . (2) When , all three rays pass through the wormhole. They enter the wormhole through the one plane and exit through the other plane (see (i-m)). (3) When , the rays are trapped within the stable orbit at the throat of the wormhole (see (n)). Furthermore, considering the incident Gaussian beam with a certain width, the spread width of the wavepacket caused by the tidal force is obtained by comparing the divergence deflection angle for both edges of the incident beam. Additionally, similar wormholes have also been emulated in waveguides on curved space for elastic waves [Citation123].

At the same time, waveguides on curved space have also inspired optical functional devices. One of the limitations of metamaterials that restricts us from exploring more complicated optical devices is the demand for either tremendous change or a singularity in the material parameter, which requires a high level of material fabrication. L. Xu et al [Citation124]. utilized waveguides on curved space to emulate topological defects with a singularity and achieve bending optical devices with a constant angle. Moreover, geodesic lenses [Citation125], which are well-known absolute instruments, were obtained by leveraging such waveguides on curved space.

5. Optical simulation of a quantum field in a curved space

Gradient waveguide arrays

Over the last decades, arrays of evanescently coupled waveguides have been brought into focus as a particular platform for manipulating electromagnetic waves, in which the dispersion and diffraction of propagating light can be specifically modified. Moreover, it turns out that the light evolution in these systems shares fundamental similarities to the quantum evolution of particle functions in the context of relativity, so that waveguide arrays can behave as a toy model system for emulating quantum fields related to GR. For instance, F. Dreisow et al [Citation126]. exploited a one-dimensional binary waveguide array to simulate the Zitterbewegung (trembling motion) of a free Dirac electron, whose dispersion curves in the vicinity of the edges of the first Brillouin zone are mathematically equivalent to the typical hyperbolic energy-momentum dispersion relation for positive-energy and negative-energy branches of a freely moving relativistic massive particle. Likewise, massless Dirac particles [Citation127], the vacuum instability and pair production [Citation128,Citation129], neutrino oscillations [Citation130] and Majorana fermions [Citation131] have also been well simulated in waveguide arrays in the spirit of the similar principle. However, these simulations belong to the realm of classic photonic simulation. It is appealing to consider the quantum aspect of analogical gravitation with the aid of quantum optics.

Fortunately, recent years have witnessed remarkable advances in quantum optics technology. Y. Wang et al [Citation132]. proposed and experimentally demonstrated a quantum evolution of fermions in close proximity to an artificial black hole, which was emulated using bilayer waveguide arrays fabricated by a femtosecond laser direct writing technique (see )). Due to the high controllability over the coupling and the on-site energy among waveguides, one-dimensional artificial black holes with the Schwarzschild metric can be constructed using evanescently coupled waveguide arrays with inhomogeneous hopping coefficients inspired by the concept of transformation optics. Compared to the linear time evolution in flat space, the dynamic behavior of a single-photon wave packet near the event horizon of the black hole has an exponential form over time, whose exponential index depends on the curvature of the black hole. More interestingly, the quantum evolution process of a single-photon wave packet in such predesigned waveguide arrays that emulate the event horizon of black hole exhibits a good analogue of Hawking Radiation: a single-photon wave packet with positive energy escapes from the black hole while negative energy is captured as shown in ). Moreover, it is an attractive question about the quantum evolution of entangled photons near the event horizon of black holes. C. Sheng et al [Citation133]. theoretically studied the quantum walks of single photons, two indistinguishable photons and entangled photons in gradient waveguide arrays that emulate the event horizon of black holes. The optical trapping of single photons and two indistinguishable photons in such nonuniform lattices conforms to a well-known classical physical recognition due to the strong gravitational force of black holes (see )). Counterintuitively, there is an optical escape for path-entangled photons for which one photon is captured, while the other photon escapes ()). Intriguingly, the counterintuitive phenomenon has a distinct escape mechanism compared to Hawking radiation, which is wholly due to quantum interference. Moreover, entanglement entropy is also studied in such a toy model system, which clearly shows that the maximally entangled state becomes less entangled and monotonically decreases with the increasing curvature of black holes.

Figure 5. Quantum simulations of general relativity phenomena with waveguide arrays. (a, b) Quantum simulation of particle pair creation near the event horizon [Citation132]. (a) The schematic of mapping the behavior of fermion pair into a waveguide array. (b) The imagined output probability distribution of single-photon wave packet and the output probability distributions with different excited positions in the same lattice. (c, d) Quantum walks of single photons and entangled photons in nonuniform waveguide arrays corresponding to a noninertial frame with the Rindler metric [Citation133]. (c) Quantum walks of single photons exhibits that photons are captured by the event horizon. (d) Quantum walks of path-entangled photons exhibits that one photon is captured, while the other photon escapes.

Figure 5. Quantum simulations of general relativity phenomena with waveguide arrays. (a, b) Quantum simulation of particle pair creation near the event horizon [Citation132]. (a) The schematic of mapping the behavior of fermion pair into a waveguide array. (b) The imagined output probability distribution of single-photon wave packet and the output probability distributions with different excited positions in the same lattice. (c, d) Quantum walks of single photons and entangled photons in nonuniform waveguide arrays corresponding to a noninertial frame with the Rindler metric [Citation133]. (c) Quantum walks of single photons exhibits that photons are captured by the event horizon. (d) Quantum walks of path-entangled photons exhibits that one photon is captured, while the other photon escapes.

6. Perspectives and conclusion

Over the years, integrated photonic chips have emerged as a well-controlled toy model system to study various large-scale gravitational astronomical phenomena in the tabletop environment. Although such analogies are just simulations, the underlying physical principle of these phenomena in both analogical systems and the real universe is supported by similar mathematical equations. Many phenomena predicted by GR are difficult to be observed and cannot be repeatedly studied. Thanks to analogical studies in the laboratory environment, one can accurately and repeatedly investigate these phenomena to verify the backing physics. Furthermore, the study of analog gravity based on photonic chips may address some more challenging works. For example, quantum gravity is a great challenge field in modern physics that seeks to describe gravity based on the principles of quantum physics and finally needs unification of quantum mechanics and general relativity. Hence, with the aid of the study of analogical gravitation based on integrated photonic chips, possible experimental works on quantum gravity will be carried out. It is noticed that a photon itself is a good quantum system and can be constructed into various types of quantum entangled states. Meanwhile, various types of curved space can be well-emulated on integrated photonic chips with the aid of metamaterials and transformation optics, such as De Sitter space and Anti-De Sitter space. The combination between the quantum entangle state and the emulated curved space on integrated photonic chips may bring much meaningful thinking to explore and understand the essence of gravity.

On the other hand, the study of analogical gravitation based on integrated photonic chips in turn has inspired various optical functional devices on chips, ranging from optical omnidirectional absorption cavities akin to black holes, nondiffractive optical beams inspired by the Einstein ring to definite nontrivial electromagnetic scatterings inspired by cosmic strings. It is noted that these gravitational inspired optical devices are achieved with the aid of the manipulation of the spatial dimension of gravitational field. The temporal dimension is another important and essential variable to describe curved spacetime. At the same time, spatiotemporal manipulation very recently has emerged as a new paradigm for wave manipulation using time-varying media [Citation134]. Therefore, if we account for the temporal dimension in the design of optical devices on integrated chips, many more fascinating optical devices and applications, inspired by the study of analogical gravitation, may be achieved, such as the nonlinear frequency conversion and a reconfigurable quantum simulator.

In conclusion, on the one hand, the research of analogical gravitation on integrated photonic chips may promote the understanding of the basic physics of gravity. On the other hand, such studies inspire us to go beyond traditional photonic devices to bring new applications. We gave a brief overview, as well as perspective, of various GR phenomena based on the integrated photonic chip. With the aid of analogical thinking about gravity, we believe that tremendous new opportunities both in basic physics of gravity and optical devices inspired by gravity still wait to be explored.

Disclosure statement

No potential conflict of interest was reported by the author(s).

Additional information

Funding

This work was financially supported by the National Natural Science Foundation of China (Grant Nos. 11690033, 12174187, 92163216, 92150302, 12147103,61425018, 11621091,11974178), and the Fundamental Research Fund for the Central Universities, China (Grant No. 14380139, 14380191).

References

  • Hasan MZ, Kane CL. Colloquium: topological insulators. Rev. Mod. Phys. 2010;82:3045.
  • Vonklitzing K, Dorda G, Pepper M. New method for high-accuracy determination of the fine-structure constant based on quantized hall resistance. Phys. Rev. Lett. 1980;45:494.
  • Lu L, Joannopoulos JD, Soljacic M. Topological photonics. Nat. Photonics. 2014;8:821.
  • Ozawa T, Price HM, Amo A, et al. Topological photonics. Rev. Mod. Phys. 2019;91:015006.
  • Huber SD. Topological mechanics. Nat. Phys. 2016;12:621.
  • Ma G, Xiao M, Chan CT. Topological phases in acoustic and mechanical systems. Nat. Rev. Phys. 2019;1:281.
  • Gordon W. Light propagation according to the relativity theory. Ann. Phys. 1923;72:421.
  • Plebanski J. Electromagnetic waves in gravitational fields. Phys. Rev. 1960;118:1396.
  • Pendry JB, Schurig D, Smith DR. Controlling electromagnetic fields. Science. 2006;312:1780.
  • Leonhardt U. Optical conformal mapping. Science. 2006;312:1777.
  • Pendry JB, Holden AJ, Robbins DJ, et al. Magnetism from conductors and enhanced nonlinear phenomena. IEEE Trans. Microwave Theory Tech. 1999;47:2075.
  • Shelby RA, Smith DR, Schultz S. Experimental verification of a negative index of refraction. Science. 2001;292:77.
  • Pendry JB. Negative refraction makes a perfect lens. Phys. Rev. Lett. 2000;85:3966.
  • Fang N, Lee H, Sun C, et al. Sub-diffraction-limited optical imaging with a silver superlens. Science. 2005;308:534.
  • Liu Z, Lee H, Xiong Y, et al. Far-field optical hyperlens magnifying sub-diffraction-limited objects. Science. 2007;315:1686.
  • Schurig D, Mock JJ, Justice BJ, et al. Metamaterial electromagnetic cloak at microwave frequencies. Science. 2006;314:977.
  • Li J, Pendry JB. Hiding under the carpet: a new strategy for cloaking. Phys. Rev. Lett. 2008;101:203901.
  • Liu R, Ji C, Mock JJ, et al. Broadband ground-plane cloak. Science. 2009;323:366.
  • Ma HF, Cui TJ. Three-dimensional broadband ground-plane cloak made of metamaterials. Nat. Commun. 2010;1:21.
  • Valentine J, Li J, Zentgraf T, et al. An optical cloak made of dielectrics. Nat. Matter. 2009;8:568.
  • Gabrielli LH, Cardenas J, Poitras CB, et al. Silicon nanostructure cloak operating at optical frequencies. Nat. Photonics. 2009;3:461.
  • Chen H, Zheng B, Shen L, et al. Ray-optics cloaking devices for large objects in incoherent natural light. Nat. Commun 2013;4:2652.
  • Narimanov EE, Kildishev AV. Optical black hole: broadband omnidirectional light absorber. Appl. Phys. Lett. 2009;95:041106.
  • Cheng Q, Cui TJ, Jiang WX, et al. An omnidirectional electromagnetic absorber made of metamaterials. New J. Phys. 2010;12:063006.
  • Genov DA, Zhang S, Zhang X. Mimicking celestial mechanics in metamaterials. Nat. Phys. 2009;5:687.
  • Chen H, Miao R-X, Li M. Transformation optics that mimics the system outside a Schwarzschild black hole. Opt. Express. 2010;18:15183.
  • Smolyaninov II, Narimanov EE. Metric signature transitions in optical metamaterials. Phys. Rev. Lett. 2010;105:067402.
  • Greenleaf A, Kurylev Y, Lassas M, et al. Electromagnetic wormholes and virtual magnetic monopoles from metamaterials. Phys. Rev. Lett. 2007;99:183901.
  • Li M, Miao R-X, Pang Y. Casimir energy, holographic dark energy and electromagnetic metamaterial mimicking de Sitter. Phys. Lett. B. 2010;689:55.
  • Smolyaninov II, Hung Y-J. Modeling of time with metamaterials. J. Opt. Soc. Am. B. 2011;28:1591.
  • Ginis V, Tassin P, Craps B, et al. Frequency converter implementing an optical analogue of the cosmological redshift. Opt. Express. 2010;18:5350.
  • Rechtsman MC, Zeuner JM, Plotnik Y, et al. Photonic Floquet topological insulators. Nature. 2013;496:196.
  • Cerjan A, Huang S, Wang M, et al. Experimental realization of a Weyl exceptional ring. Nat. Photonics. 2019;13:623.
  • Rechtsman MC, Zeuner JM, Tuennermann A, et al. Strain-induced pseudomagnetic field and photonic Landau levels in dielectric structures. Nat. Photonics. 2013;7:153.
  • Mukherjee S, Di Liberto M, Oehberg P, et al. Experimental observation of aharonov-Bohm cages in photonic lattices. Phys. Rev. Lett. 2018;121:075502.
  • Mukherjee S, Spracklen A, Choudhury D, et al. Observation of a localized flat-band state in a photonic lieb lattice. Phys. Rev. Lett. 2015;114:245504.
  • Weimann S, Perez-Leija A, Lebugle M, et al. Implementation of quantum and classical discrete fractional Fourier transforms. Nat. Commun. 2016;7:11027.
  • Zeuner J, Sharma AN, Tillmann M, et al. Integrated-optics heralded controlled-NOT gate for polarization-encoded qubits. NPJ Quantum Inf. 2018;4:13.
  • Longhi S. Quantum-optical analogies using photonic structures. Laser Photonics Rev. 2009;3:243.
  • Szameit A, Nolte S. Discrete optics in femtosecond-laser-written photonic structures. J. Phys. B. 2010;43:163001.
  • Weinfurtner S, Tedford EW, Penrice MCJ, et al. Measurement of stimulated hawking emission in an analogue system. Phys. Rev. Lett. 2011;106:021302.
  • Torres T, Patrick S, Coutant A, et al. Rotational superradiant scattering in a vortex flow. Nat. Phys. 2017;13:833.
  • Zou D, Chen T, He W, et al. Observation of hybrid higher-order skin-topological effect in non-hermitian topolectrical circuits. Nat. Commun. 2021;12:7201.
  • Feynman RP. Simulating physics with computers. Int. J. Theor. Phys. 1982;21:467.
  • Georgescu IM, Ashhab S, Nori F. Quantum simulation. Rev. Mod. Phys. 2014;86:153.
  • Regal CA, Greiner M, Jin DS. Observation of resonance condensation of fermionic atom pairs. Phys. Rev. Lett. 2004;92:040403.
  • Buchler HP, Hermele M, Huber SD, et al. Atomic quantum simulator for lattice gauge theories and ring exchange models. Phys. Rev. Lett. 2005;95:040402.
  • Kassal I, Jordan SP, Love PJ, et al. Polynomial-time quantum algorithm for the simulation of chemical dynamics. Proc. Natl. Acad. Sci. USA. 2008;105:18681.
  • Fischer UR, Schutzhold R. Quantum simulation of cosmic inflation in two-component bose-Einstein condensates. Phys. Rev. A. 2004;70:063615.
  • Hawking SW. Black-hole explosions. Nature. 1974;248:30.
  • Hu J, Feng L, Zhang Z, et al. Quantum simulation of Unruh radiation. Nat. Phys. 2019;15:785.
  • Steinhauer J. Observation of quantum Hawking radiation and its entanglement in an analogue black hole. Nat. Phys. 2016;12:959.
  • Horstmann B, Reznik B, Fagnocchi S, et al. Hawking radiation from an acoustic black hole on an ion ring. Phys. Rev. Lett. 2010;104:250403.
  • Giovanazzi S. Hawking radiation in sonic black holes. Phys. Rev. Lett. 2005;94:061302.
  • Nation PD, Blencowe MP, Rimberg AJ, et al. Analogue hawking radiation in a dc-SQUID array transmission line. Phys. Rev. Lett. 2009;103:087004.
  • Drori J, Rosenberg Y, Bermudez D, et al. Observation of stimulated hawking radiation in an optical analogue. Phys. Rev. Lett. 2019;122:010404.
  • Philbin TG, Kuklewicz C, Robertson S, et al. Fiber-optical analog of the event horizon. Science. 2008;319:1367.
  • Jaksch D, Bruder C, Cirac JI, et al. Cold bosonic atoms in optical lattices. Phys. Rev. Lett. 1998;81:3108.
  • Yang B, Sun H, Ott R, et al. Observation of gauge invariance in a 71-site bose-Hubbard quantum simulator. Nature. 2020;587:392.
  • Lidar DA, Biham O. Simulating ising spin glasses on a quantum computer. Physl. Rev. E. 1997;56:3661.
  • Greiner M, Mandel O, Esslinger T, et al. Quantum phase transition from a superfluid to a mott insulator in a gas of ultracold atoms. Nature. 2002;415:39.
  • Arute F, Arya K, Babbush R, et al. Quantum supremacy using a programmable superconducting processor. Nature. 2019;574:505.
  • Zhong H-S, Wang H, Deng Y-H, et al. Quantum computational advantage using photons. Science. 2020;370:1460.
  • Wu Y, Bao W-S, Cao S, et al. Strong quantum computational advantage using a superconducting quantum processor. Phys. Rev. Lett. 2021;127:180501.
  • Leonhardt U, Philbin TG. General relativity in electrical engineering. New J. Phys. 2006;8:247.
  • Wang J, Sciarrino F, Laing A, et al. Integrated photonic quantum technologies. Nat. Photonics. 2020;14:273.
  • Politi A, Cryan MJ, Rarity JG, et al. Silica-on-silicon waveguide quantum circuits. Science. 2008;320:646.
  • Crespi A, Ramponi R, Osellame R, et al. Integrated photonic quantum gates for polarization qubits. Nat. Commun. 2011;2:566.
  • Lu X, Li Q, Westly DA, et al. Chip-integrated visible-telecom entangled photon pair source for quantum communication. Nat.Phys. 2019;15:373.
  • Wang C, Zhang M, Chen X, et al. Integrated lithium niobate electro-optic modulators operating at CMOS-compatible voltages. Nature. 2018;562:101.
  • Wang J, Santamato A, Jiang P, et al. Gallium arsenide (GaAs) quantum photonic waveguide circuits. Opt. Commun. 2014;327:49.
  • Abellan C, Amaya W, Domenech D, et al. Quantum entropy source on an InP photonic integrated circuit for random number generation. Optica. 2016;3:989.
  • Xu H, Dai D, Shi Y. Ultra-broadband and ultra-compact on-chip silicon polarization beam splitter by using hetero-anisotropic metamaterials. Laser Photonics Rev. 2019;13:1800349.
  • Wu S, Xiao J. Compact polarization rotator for silicon-based cross-slot waveguides using subwavelength gratings. Appl. Opt. 2017;56:4892.
  • Xu H, Shi Y. Subwavelength-grating-assisted silicon polarization rotator covering all optical communication bands. Opt. Express. 2019;27:5588.
  • Gabrielli LH, Liu D, Johnson SG, et al. On-chip transformation optics for multimode waveguide bends. Nat. Commun. 2012;3:1217.
  • Wang Y, Sheng C, Liu H, et al. Transformation bending device emulated by graded-index waveguide. Opt. Express. 2012;20:13006.
  • Li Z, Kim M-H, Wang C, et al. Controlling propagation and coupling of waveguide modes using phase-gradient metasurfaces. Nature Nanotech. 2017;12:675.
  • Liu Y, Xu K, Wang S, et al. Arbitrarily routed mode-division multiplexed photonic circuits for dense integration. Nat. Commun. 2019;10:3263.
  • Mohanty A, Zhang M, Dutt A, et al. Quantum interference between transverse spatial waveguide modes. Nat. Commun. 2017;8:14010.
  • Zentgraf T, Liu Y, Mikkelsen MH, et al. Plasmonic luneburg and Eaton lenses. Nature Nanotech. 2011;6:151.
  • Li S, Zhou Y, Dong J, et al. Universal multimode waveguide crossing based on transformation optics. Optica. 2018;5:1549.
  • Zentgraf T, Valentine J, Tapia N, et al. An optical “Janus” device for integrated photonics. Adv. Mater. 2010;22:2561.
  • Choo H, Kim M-K, Staffaroni M, et al. Nanofocusing in a metal-insulator-metal gap plasmon waveguide with a three-dimensional linear taper. Nat. Photonics. 2012;6;838.
  • Wang X, Chen H, Liu H, et al. Self-focusing and the Talbot effect in conformal transformation optics. Phys. Rev. Lett. 2017;119:033902.
  • Smolyaninov II, Smolyaninova VN, Kildishev AV, et al. Anisotropic metamaterials emulated by tapered waveguides: application to optical cloaking. Phys. Rev. Lett. 2009;102:213901.
  • Cang H, Salandrino A, Wang Y, et al. Adiabatic far-field sub-diffraction imaging. Nat. Commun. 2015;6:7942.
  • Sheng C, Liu H, Wang Y, et al. Trapping light by mimicking gravitational lensing. Nat. Photonics. 2013;7:902.
  • Sheng C, Bekenstein R, Liu H, et al. Wavefront shaping through emulated curved space in waveguide settings. Nat. Commun. 2016;7:10747.
  • Smolyaninov II. Unruh effect in a waveguide. Phys. Lett. A. 2008;372:5861.
  • Smolyaninov II. Giant Unruh effect in hyperbolic metamaterial waveguides. Opt. Lett. 2019;44:2224.
  • Ge H, Sheng C, Zhu S, et al. Observation of the acceleration of light in a tapered optical fiber. Opt. Express. 2021;29:27212.
  • Wang X, Liu H, Sheng C, et al. Reversible wavefront shaping between gaussian and airy beams by mimicking gravitational field. J. Opt. 2018;20:024015.
  • Meng Y, Chen Y, Lu L, et al. Optical meta-waveguides for integrated photonics and beyond. Light-Science & Applications. 2021;10:235.
  • Sheng C, Liu H, Zhu S. Transformation optics based on metasurfaces. Sci. Bull. 2019;64:793.
  • Abbott BP, Abbott R, Abbott TD, et al. Observation of gravitational waves from a binary black hole merger. Phys. Rev. Lett. 2016;116:061102.
  • Abbott BP, Abbott R, Abbott TD, et al. GW170817: observation of gravitational waves from a binary neutron star inspiral. Phys. Rev. Lett. 2017;119:161101.
  • Blasi S, Brdar V, Schmitz K. Has NANOGrav found first evidence for cosmic strings? Phys. Rev. Lett. 2021;126:041305.
  • Ellis J, Lewicki M. Cosmic string interpretation of NANOGrav pulsar timing data. Phys. Rev. Lett. 2021;126:041304.
  • Liu J, Cai R-G, Guo Z-K. Large anisotropies of the stochastic gravitational wave background from cosmic domain walls. Phys. Rev. Lett. 2021;126:141303.
  • Zhang Y-L, Pendry JB, Lei DY. Radial anisotropy from a geometric viewpoint: topological singularity and effective medium realization. Phys. Rev. B. 2017;96:035430.
  • Fernandez-Nunez I, Bulashenko O. Emergence of fresnel diffraction zones in gravitational lensing by a cosmic string. Phys. Lett. A. 2017;381:1764.
  • Zhang Y-L, Dong X-Z, Zheng M-L, et al. Steering electromagnetic beams with conical curvature singularities. Opt. Lett. 2015;40:4783.
  • Fernandez-Nunez I, Bulashenko O. Wave propagation in metamaterials mimicking the topology of a cosmic string. J. Opt. 2018;20:045603.
  • Alford MG, Wilczek F. Aharonov-Bohm interaction of cosmic strings with matter. Phys. Rev. Lett. 1989;62:1071.
  • Yang Y. Cosmic strings in a product abelian gauge field theory. Nucl. Phys. B. 2014;885:25.
  • Gerbert PD. Fermions in an aharonov-Bohm field and cosmic strings. Phys. Rev. D. 1989;40:1346.
  • Mackay TG, Lakhtakia A. Towards a metamaterial simulation of a spinning cosmic string. Phys. Lett. A. 2010;374:2305.
  • Sheng C, Liu H, Chen H, et al. Definite photon deflections of topological defects in metasurfaces and symmetry-breaking phase transitions with material loss. Nat. Commun. 2018;9:4271.
  • Genevet P, Wintz D, Ambrosio A, et al. Controlled steering of Cherenkov surface plasmon wakes with a one-dimensional metamaterial. Nature Nanotech. 2015;10:804.
  • Zhong F, Li J, Liu H, et al. Controlling surface plasmons through covariant transformation of the spin-dependent geometric phase between curved metamaterials. Phys. Rev. Lett. 2018;120:243901.
  • Dacosta RCT. Quantum-mechanics of a constrained particle. Phys. Rev. A. 1981;23:1982.
  • Batz S, Peschel U. Linear and nonlinear optics in curved space. Phys. Rev. A. 2008;78:043821.
  • Schultheiss VH, Batz S, Szameit A, et al. Optics in curved space. Phys. Rev. Lett. 2010;105:143901.
  • Schultheiss VH, Batz S, Peschel U. Hanbury brown and Twiss measurements in curved space. Nat. Photonics. 2016;10:106.
  • Bekenstein R, Nemirovsky J, Kaminer I, et al. Shape-preserving accelerating electromagnetic wave packets in curved space. Phys. Rev. X. 2014;4,:011038.
  • Patsyk A, Bandres MA, Bekenstein R, et al. Observation of accelerating wave packets in curved space. Phys. Rev. X. 2018;8:011001.
  • Xu C, Wang L-G. Gouy and spatial-curvature-induced phase shifts of light in two-dimensional curved space. New J. Phys. 2019;21:113013.
  • Liang G-H, Cai RG, Ma Y-Z, et al. Mimicking an expanding universe by optical interference in a helicoid waveguide. Opt. Express. 2020;28:11406.
  • Ding W, Wang Z. Laser propagation in a Rindler accelerated reference frame based on matrix optics. Opt. Express. 2021;29:28631.
  • Bekenstein R, Kabessa Y, Sharabi Y, et al. Control of light by curved space in nanophotonic structures. Nat. Photonics. 2017;11:664.
  • He GHLRQ, Zhu SN, Liu H. Simulation of giant tidal force of wormhole using curved optical spaces. Physcial Review Research. 2020;2:013237.
  • Zhu J, Liu Y, Liang Z, et al. Elastic waves in curved space: mimicking a wormhole. Phys. Rev. Lett. 2018;121:234301.
  • Xu L, He R, Yao K, et al. Conformal singularities and topological defects from inverse transformation optics. Phys. Rev. Appl. 2019;11:034072.
  • Xu L, Wang X, Tyc T, et al. Light rays and waves on geodesic lenses. Photon. Res. 2019;7:1266.
  • Dreisow F, Heinrich M, Keil R, et al. Classical simulation of relativistic zitterbewegung in photonic lattices. Phys. Rev. Lett. 2010;105:143902.
  • Zeuner JM, Efremidis NK, Keil R, et al. Optical analogues for massless dirac particles and conical diffraction in one dimension. Phys. Rev. Lett. 2012;109:023602.
  • Dreisow F, Longhi S, Nolte S, et al. Vacuum instability and pair production in an optical setting. Phys. Rev. Lett. 2012;109:110401.
  • Koke C, Noh C, Angelakis DG. Dirac equation in 2-dimensional curved spacetime, particle creation, and coupled waveguide arrays. Ann. Phys. 2016;374:162.
  • Marini A, Longhi S, Biancalana F. Optical simulation of neutrino oscillations in binary waveguide arrays. Phys. Rev. Lett. 2014;113:150401.
  • Keil R, Noh C, Rai A, et al. Optical simulation of charge conservation violation and Majorana dynamics. Optica. 2015;2:454.
  • Wang Y, Sheng C, Lu Y-H, et al. Quantum simulation of particle pair creation near the event horizon. Natl. Sci. Rev. 2020;7:1476.
  • Sheng C, Huang C, Yang R, et al. Simulating the escape of entangled photons from the event horizon of black holes in nonuniform optical lattices. Phys. Rev. A. 2021;103:033703.
  • Galiffi E, Tirole R, Yin S, et al. Photonics of time-varying media. Advanced Photonics. 2022;4:014002.