Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Topological laser with higher-order corner states in the 2-dimensional Su-Schrieffer-Heeger model

Open Access Open Access

Abstract

A nonlinear non-Hermitian topological laser system based on the higher-order corner states of the 2-dimensional (2D) Su-Schrieffer-Heeger (SSH) model is investigated. The topological property of this nonlinear non-Hermitian system described by the quench dynamics is in accordance with that of a normal 2D SSH model. In the topological phase, all sites belonging to the topological corner states begin to emit stable laser light when a pulse is given to any one site of the lattice, while no laser light is emitted when the lattice is in the trivial phase. Furthermore, the next-nearest-neighbor (NNN) couplings are introduced into the strong-coupling unit cells of the 2D SSH model, which open a band gap in the continuous band structure. In the topological phase, similar to the case of 2D SSH model without NNN couplings, the corner sites can emit stable laser light due to the robustness of the higher-order corner states when the NNN couplings are regarded as the perturbation. However, amplitude of the stimulated site does not decay to zero in the trivial phase, because the existence of the NNN couplings in the strong-coupling unit cells make the lattice like one in the tetramer limit, and a weaker laser light is emitted by each corner.

© 2023 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement

1. Introduction

Topological insulators (TIs) [1,2], including Chern insulators [3,4], time-reversal-invariant TIs [5,6] and topological crystalline insulators [7], have attracted much interest in recent years due to their robust topological edge states, forming a hot topic in condensed-matter physics. Then, the concept to higher-order TIs (HOTIs) has been generalized [818]. For example, a second-order TI is a $d$-dimensional ($d$D) insulator which has $(d-1)$D topological boundary states but also $(d-2)$D topological boundary states. These topological edge states have been investigated or observed in electric circuits [1923], phonics [2426] and photonics models [2729].

Recently, laser systems based on topological edge states have been investigated through photonics models [3047]. A typical photonics model is a coupled-ring array made of active resonators which has been widely used in topological photonics [4850]. 2D TIs based on 2D lattices [32,33,3847], including Chern insulators, valley topological insulators, and topological crystal insulators, provide a well scheme to realizing large-scale topological laser arrays with a higher power. These models generate laser light along the edge belonging to edge modes of a 2D lattice. Furthermore, topological laser based on the higher-order corner states of kagome lattice [46] as well as corner states of 2D lattice with non-Hermitian skin effect [47] has been investigated, where stable laser light is emitted by the sites belonging to the topological corner states.

In this paper, we consider a laser system based on the 2D Su-Schrieffer-Heeger (SSH) model [51], which is composed of coupled vertical and horizontal stacks of SSH chains [52]. Compared with the kagome lattice [17,46], the 2D SSH model has a simpler geometry but a more special band structure, where the topological phase transition happens with the band inversion but without a band gap opened on either side of the phase transition point. Due to the continuousness of the band structure near the zero energy, for a finite structure of 2D SSH model in the topological phase, the zero-energy eigenvalues are buried within continuous bulk bands, with the second-order topological corner states called the bound states in the continuum (BICs) [53]. Through the dynamics of the laser system in the 2D SSH model, stable laser light is emitted by each corner site of the lattice in the topological phase, while no laser light emitted in the trivial phase. Besides, the next-nearest-neighbor (NNN) couplings in a coupled-ring array made of active resonators may not be ignored sometimes. The NNN couplings are introduced into the strong-coupling unit cells of the 2D SSH model [54,55], and the dynamics of the topological laser system in this situation is discussed. When the NNN couplings exist in the strong-coupling unit cells, a band gap in the band structure of the 2D SSH model in the topological phase will be opened, and the topological zero-energy modes will be located in this band gap. This characteristic is also absent in the kagome lattice for the band gap is opened beside the phase transition point, and NNN couplings cannot be introduced into the unit cell with three sites in the kagome lattice. Similar to the result of the 2D SSH model without NNN couplings in the topological phase, the corner sites can emit stable laser light. While the results of the two cases are quite different in the trivial phase. The NNN couplings make the 2D SSH model more like a set of tetramers, and laser lights can also be emitted by the corner sites when any one site is stimulated.

The remainder of this paper is organized as follows. In Section 2, we describe the 2D SSH model and investigate the dynamics of a topological laser system with nonlinear non-Hermitian gain and loss on it. In Section 3, the NNN couplings are introduced into the strong-coupling unit cells of the 2D SSH model, and dynamics of the laser system on it is also discussed. Conclusions are presented in Section 4.

2. Topological laser in 2D SSH model

The lattice of the 2D SSH model consisting of four sites per unit cell with dimerized nearest-neighbor couplings is shown in Fig. 1, where $\kappa _A$ is the intracell hopping and $\kappa _B$ is the intercell hopping. The Hamiltonian in the real space is written as

$$\begin{aligned}H=&\sum_{x,y}\Big[\kappa(1-({-}1)^{x}\lambda)a_{x,y}^{\dagger}a_{x+1,y}^{ }\\ &+\kappa(1-({-}1)^{y}\lambda)a_{x,y}^{\dagger}a_{x,y+1}^{ }\Big]+\mathrm{H.c.}, \end{aligned}$$
where $a_{x,y}^{\dagger}$ ($a_{x,y}^{ }$) is the creation (annihilation) operator of site $(x,y)$. Then the hoppings are $\kappa _A=\kappa (1+\lambda )$ and $\kappa _B=\kappa (1-\lambda )$. We consider a square lattice of the size $L\times L$, indicating that the square lattice has $L^2$ unit cells with the number of sites $N=4L^2$. The site numbers are marked as Fig. 1 shows with $A$, $B$, $C$ and $D$ the four corner sites.

 figure: Fig. 1.

Fig. 1. Schematics of a finite lattice of 2D SSH model. Each unit cell contains four sites, and $\kappa _A$ ($\kappa _B$) is the intracell (intercell) hopping.

Download Full Size | PDF

It is known the 2D SSH model has chiral symmetry as well as $C_{4v}$ symmetry, and band structure in the lattice is in a topological phase for $-1<\lambda <0$ and in a trivial phase for $0<\lambda <1$. In the following, we take a square lattice with the size $L=6$. The energy eigenvalues $E_{n}$ ($n=1,2,\ldots,N$) of such a finite structure, in units of $\kappa$, are shown in Fig. 2(a). The second-order topological corner states corresponding to the zero-energy modes (red lines) emerge for $-1<\lambda <0$, and the density of a corner state is shown in Fig. 2(b), where the density of the state becomes arbitrarily small except for the four corners. These corner states are BICs due to these zero-energy modes existing in the continuous band structure. Figures 2(c) and 2(d) are two bulk states of the 2D SSH model in the trivial phase with $\lambda =0.5$, and the densities of them are distributed on both corner and bulk sites.

 figure: Fig. 2.

Fig. 2. (a) The energy spectrum of a finite lattice of 2D SSH model with the size $L=6$. (b) Density of a corner state corresponding to a zero-energy mode in the topological phase with $\lambda =-0.5$. (c),(d) Density of two bulk states in the trivial phase with $\lambda =0.5$.

Download Full Size | PDF

The finite lattice of the 2D SSH model can be realized through a coupled-ring array made of active resonators. The Hamiltonian $H$ can be written as a $4L^{2}\times 4L^{2}$ matrix $\{H_{mn}\}_{4L^{2}\times 4L^{2}}$, then the dynamics of a laser system on the coupled-ring array is expressed by

$$\mathrm{i}\frac{\mathrm{d}\psi_{n}}{\mathrm{d}t}=\sum_{m=1}^{4L^{2}}H_{nm}\psi_{m}-\mathrm{i}\gamma\left(1-\frac{\xi P_{n}}{1+|\psi_{n}|^{2}/\eta}\right)\psi_{n},$$
where $\psi _{n}$ ($n=1,2,3,\ldots,N$) is the amplitude of the site $n$, $H_{nm}$ is the matrix element of $H$, $\gamma$ is the loss in each resonator, $\gamma \xi$ can be treated as the optical gain via stimulated emission, $\eta$ is a nonlinear coefficient, and $P_{n}$ stands for the spatial profile of the pump. We consider the optical gains are introduced only at the corner sites, then we have
$$P_{n}=\sum_{n_{c}}\delta_{n,n_{c}},$$
where $n_{c}=n_{A},n_{B},n_{C},n_{D}$.

The property of such a nonlinear non-Hermitian system can be expressed by the quench dynamics [46,5659], where a pulse is given to one site, and its time evolution under Eq. (2) is explored. We numerically solve Eq. (2) under the initial condition $\psi _{1}(0)=1$, i.e. a pulse is given to the corner $A$ at $t=0$. A good signal to detect whether the system is topological or trivial is the saturated amplitude $|\psi _{1}|$ after enough time. With the loss or gain in each resonator, the system will reach a stationary state as the system evolves for enough time, and the time is chosen as $t=60$ here after our numerical analysis. The saturated amplitude $|\psi _{1}|$ for different losses $\gamma$ is plotted in Fig. 3. It is shown that the amplitude $|\psi _{1}|$ is a finite value when the 2D SSH model is in the topological phase, while it is almost zero in the trivial phase.

 figure: Fig. 3.

Fig. 3. The saturated amplitude $|\psi _{1}|$ at the corner $A$ after enough time. The parameters are set as $\kappa =1$, $\xi =2$ and $\eta =1$.

Download Full Size | PDF

Time evolution of the amplitude $|\psi _{n}|$ of the 2D SSH model with the initial condition $\psi _{1}(0)=1$ is plotted in Fig. 4. Due to the existence of the loss or gain in each resonator, the amplitude $|\psi _{n}|$ of each site reaches a stationary solution. In the topological phase, the density of state is mainly localized at the four corner sites, and the four saturated amplitudes are almost identical. While in the trivial phase, the four saturated amplitudes are almost zero. In the trivial phase, the corner states do not exist, and the bulk states are distributed on both corner and bulk sites as shown in Figs. 2(c) and 2(d). The system evolves under the bulk state causing the energy flowing into the bulk site, and the loss in each resonator quickly dissipates the energy.

 figure: Fig. 4.

Fig. 4. Time evolution of the amplitude $|\psi _{n}|$ of the 2D SSH model. The parameters are set as $\kappa =1$, $\gamma =1/2$, $\xi =2$, $\eta =1$, (a1)-(a5) $\lambda =-0.5$ for the topological phase and (b1)-(b5) $\lambda =0.5$ for the trivial phase.

Download Full Size | PDF

Figure 5 is density plot of time evolution of the amplitude $|\psi _{n}|$ of the side $AB$. For convenience, we re-mark the numbers of the sites along the side $AB$, which are marked from $1$ to $12$. The amplitudes $|\psi _{1}|$ and $|\psi _{12}|$ both reach to a finite value after enough time for the initial condition $\psi _{1}(0)=1$ and $\psi _{6}(0)=1$ in the topological phase in Figs. 5(a) and 5(c), while quickly decay to zero in the trivial phase in Figs. 5(b) and 5(d).

 figure: Fig. 5.

Fig. 5. Density plot of time evolution of the amplitude $|\psi _{n}|$ of the side $AB$. The sites on the side are marked from $1$ to $12$. (a),(b) The initial condition $\psi _{1}(0)=1$, and (c),(d) the initial condition is $\psi _{6}(0)=1$. The parameters are set as $\kappa =1$, $\gamma =1/2$, $\xi =2$, $\eta =1$, (a),(c) $\lambda =-0.5$ for the topological phase and (b),(d) $\lambda =0.5$ for the trivial phase.

Download Full Size | PDF

Time evolution of several amplitudes $|\psi _{n}|$ is plotted in Fig. 6. In the topological phase, the amplitude $|\psi _{1}|$ at the corner site $A$ is rapidly saturated as shown in Fig. 6(a), and other amplitudes is saturated with a delay, which implies the propagation of a wave with a tiny amplitude through the square lattice although it is hard to see the propagation of the wave in Fig. 5(a). The propagation of this wave induces a stimulated emission at the other three corner sites. The delay time of amplitude $|\psi _{144}|$ is longer than that of amplitude $|\psi _{22}|$, which means that the wave propagates through a longer path from corner site $A$ to $B$ than from corner site $A$ to $D$. Delay time of the saturated amplitudes (a) $|\psi _{B}|$ of the corner site $B$ and (b) $|\psi _{D}|$ of the corner site $D$ with different values of the hopping amplitude $\kappa$ is plotted in Fig. 7, and it is shown that increasing the value of $\kappa$ can shorten the delay time. Besides, in the trivial phase, the amplitude $|\psi _{1}|$ quickly decays to zero in Fig. 6(b), and there is no excitation at the other three corner sites.

 figure: Fig. 6.

Fig. 6. Time evolution of several amplitudes $|\psi _{n}|$. The parameters are set as $\kappa =1$, $\gamma =1/2$, $\xi =2$, $\eta =1$, (a) $\lambda =-0.5$ for the topological phase and (b) $\lambda =0.5$ for the trivial phase.

Download Full Size | PDF

 figure: Fig. 7.

Fig. 7. Delay time of the saturated amplitudes (a) $|\psi _{B}|$ and (b) $|\psi _{D}|$ with $\kappa =1/4, 1/2, 1,\ldots,2$. The parameters are set as $\kappa =1$, $\gamma =1/2$, $\xi =2$, $\eta =1$, and $\lambda =-0.5$

Download Full Size | PDF

Moreover, we study the tetramer limit ($\kappa _{A}=1$ and $\kappa _{B}=0$) in the trivial phase. Then the system is decomposed into a set of tetramers with only the intracell hopping $\kappa _A$. Taking the initial condition $\psi _{1}(0)=1$, we obtain the equations

$$\mathrm{i}\frac{\mathrm{d}\psi_{1}}{\mathrm{d}t}=\kappa_{A}(\psi_{2}+\psi_{3})-\mathrm{i}\gamma\left(1-\frac{\xi}{1+|\psi_{1}|^{2}/\eta}\right)\psi_{1},$$
$$\mathrm{i}\frac{\mathrm{d}\psi_{2}}{\mathrm{d}t}=\kappa_{A}(\psi_{1}+\psi_{4})-\mathrm{i}\gamma\psi_{2},$$
$$\mathrm{i}\frac{\mathrm{d}\psi_{3}}{\mathrm{d}t}=\kappa_{A}(\psi_{1}+\psi_{4})-\mathrm{i}\gamma\psi_{3},$$
$$\mathrm{i}\frac{\mathrm{d}\psi_{4}}{\mathrm{d}t}=\kappa_{A}(\psi_{2}+\psi_{3})-\mathrm{i}\gamma\psi_{4}.$$

We numerically solve Eqs. (4)–(7) and plot the time evolution of the amplitudes $|\psi _{n}|$ ($n=1,2,3,4$) in Figs. 8(a) and 8(b). Amplitude $|\psi _{1}|$ approaches a saturated value, which is consistent with Fig. 3 when $\lambda$ is close to $1$. Analytically stationary solutions of Eqs. (4)–(7) are too complex, so we remove the losses on sites $2$, $3$ and $4$ to obtain a simple analytical solution. Removing the losses on sites $2$, $3$ and $4$, Eqs. (4)–(7) turn into

$$\mathrm{i}\frac{\mathrm{d}\psi_{1}}{\mathrm{d}t}=\kappa_{A}(\psi_{2}+\psi_{3})-\mathrm{i}\gamma\left(1-\frac{\xi}{1+|\psi_{1}|^{2}/\eta}\right)\psi_{1},$$
$$\mathrm{i}\frac{\mathrm{d}\psi_{2}}{\mathrm{d}t}=\kappa_{A}(\psi_{1}+\psi_{4}),$$
$$\mathrm{i}\frac{\mathrm{d}\psi_{3}}{\mathrm{d}t}=\kappa_{A}(\psi_{1}+\psi_{4}),$$
$$\mathrm{i}\frac{\mathrm{d}\psi_{4}}{\mathrm{d}t}=\kappa_{A}(\psi_{2}+\psi_{3}).$$

Then, the stationary solutions of Eqs. (8)–(11) have a simple form and are given by

$$\begin{aligned} \psi_{1}&={-}\psi_{4}=\sqrt{\eta(\xi-1)},\\ \psi_{2}&={-}\psi_{3}=\mathrm{const}. \end{aligned}$$

The numerical solutions are shown in Figs. 8(c) and 8(d), amplitudes $|\psi _{1}|$ and $|\psi _{4}|$ approach a saturated value (12).

 figure: Fig. 8.

Fig. 8. Time evolution of the amplitudes $|\psi _{n}|$ of the tetramer in 2D SSH model (a),(b) with and (c),(d) without losses on sites $2$, $3$ and $4$. The parameters are set as $\kappa _{A}=1$, $\kappa _{B}=0$, $\gamma =1/2$, $\xi =2$ and $\eta =1$.

Download Full Size | PDF

3. 2D SSH model with NNN couplings

The NNN couplings in a coupled-ring array made of active resonators may not always be ignored as a matter of fact. Photonics systems in such a 2D form can be realized by a 2D coupled waveguide array [22] or a coupled-ring array [32,33,4850], and the sites are usually coupled by evanescent waves or link rings. In these cases, the coupling between two sites always depends on the distance between them. Hence, in the strong-coupling unit cells, the sites are close to each other, and the NNN couplings in a coupled-ring array may be enhanced. In this section, we consider the 2D SSH model with NNN couplings in the strong-coupling unit cells. As is shown in Fig. 9, the next-nearest-neighbor couplings $\chi$ exist in the unit cells with the hopping $\kappa _{B}$ for $\kappa _{A}<\kappa _{B}$ (topological), and in the unit cells with the hopping $\kappa _{A}$ for $\kappa _{A}>\kappa _{B}$ (trivial). In the following, we take $\chi =\kappa _{B}$ for $\kappa _{A}<\kappa _{B}$ and $\chi =\kappa _{A}$ for $\kappa _{A}>\kappa _{B}$. The size of the 2D SSH model NNN couplings is chosen to be $L=6$, which is the same as that of the finite lattice of 2D SSH model in Fig. 1. The sites numbers are also marked as in Fig. 1.

 figure: Fig. 9.

Fig. 9. 2D SSH model with NNN couplings $\chi$ in the strong-coupling unit cells. (a) The strong-coupling unit cells are cells with hopping $\kappa _{B}$, and (b) the strong-coupling unit cells are cells with hopping $\kappa _{A}$.

Download Full Size | PDF

Figure 10(a) is the energy spectrum of a finite lattice of 2D SSH model with NNN couplings in cells with coupling $\kappa _{A}$, and Fig. 10(b) is the energy spectrum of the model with NNN couplings in cells with coupling $\kappa _{B}$. The corner states corresponding to the zero-energy modes also emerge for $-1<\lambda <0$ and disappear for $0<\lambda <1$. Figures 10(c) and 10(d) give the energy eigenvalues of 2D SSH model without NNN couplings and with NNN couplings in the topological phase with $\lambda =-0.5$. It is known the corner states of 2D SSH model in the topological phase are BICs, and it can be found the energy spectrum in Fig. 10(c) is almost continuous and inset with the topological zero-energy modes (red points). It is not easy to distinguish the topological zero-energy eigenvalues and the zero-energy eigenvalues of the bulk states. Differently, the topological zero-energy modes are located in a band gap in Fig. 10(d). Moreover, introducing NNN couplings is also a method used to open a band gap in the continuous band structure of 2D SSH model. The density of a corner states corresponding to one of the zero-energy eigenvalues in Fig. 10(d) are shown in Fig. 11(a), and Fig. 11(b) is a bulk state of the trivial 2D SSH model with NNN couplings.

 figure: Fig. 10.

Fig. 10. Energy spectrum of a finite lattice of 2D SSH model (a) with NNN couplings $\chi$ in cells with coupling $\kappa _{A}$ and (b) with NNN couplings $\chi$ in cells with coupling $\kappa _{B}$. Energy eigenvalues of 2D SSH model of size $L=6$ (c) without NNN couplings and (d) with NNN couplings in the topological phase. Red points are the zero-energy eigenvalues of the topological corner states. The parameters in (c),(d) are set as $\kappa =1$ and $\lambda =-0.5$.

Download Full Size | PDF

 figure: Fig. 11.

Fig. 11. (a) The density of a corner state corresponding to one of the zero-energy eigenvalues in Fig. 10(d). (b) The density of a bulk state of the 2D SSH model with NNN couplings in the trivial phase ($\lambda =0.5$).

Download Full Size | PDF

With the initial condition $\psi _{1}(0)=1$, time evolution of the amplitude $|\psi _{n}|$ of the 2D SSH model with NNN couplings is plotted in Fig. 12. Similar to the results of the topological phase in Fig. 4, the amplitude $|\psi _{n}|$ of each site also reaches a stationary solution due to the robustness of the higher-order corner states. While in the trivial phase, the four saturated amplitudes of the corner sites do not decay to zero, and also reach a nonzero stationary solution, which is different from that in Fig. 4 for trivial phase. The results of trivial phase are caused by the corner-like state shown in Fig. 11(b). As is shown in Fig. 11(b), although this eigenstate is a bulk state, its density is mainly distribute on the sites of the cells at four corners. And the system in the trivial phase will evolve into a stationary state like this corner-like state when the first site is stimulated.

 figure: Fig. 12.

Fig. 12. Time evolution of the amplitude $|\psi _{n}|$ of the 2D SSH model with NNN couplings. The parameters are set as $\kappa =1$, $\gamma =1/2$, $\xi =2$, $\eta =1$, (a1)-(a5) $\lambda =-0.5$ for the topological phase and (b1)-(b5) $\lambda =0.5$ for the trivial phase.

Download Full Size | PDF

The density plot of time evolution of the amplitude $|\psi _{n}|$ of the side $AB$ of the 2D SSH model with NNN couplings is shown in Fig. 13. The amplitudes $|\psi _{1}|$ and $|\psi _{12}|$ also reach to a finite value after enough time in the topological phase in Figs. 13(a) and 13(c), like the results in Figs. 5(a) and 5(c). However, this phenomenon also occurs in the trivial phase when NNN couplings are introduced into the model, which is different from the results of model without NNN couplings.

 figure: Fig. 13.

Fig. 13. Density plot of time evolution of the amplitude $|\psi _{n}|$ of the side $AB$ of the 2D SSH model with NNN couplings. The sites on the side are marked from $1$ to $12$. (a),(b) The initial condition $\psi _{1}(0)=1$, and (c),(d) the initial condition is $\psi _{6}(0)=1$. The parameters are set as $\kappa =1$, $\gamma =1/2$, $\xi =2$, $\eta =1$, (a),(c) $\lambda =-0.5$ for the topological phase and (b),(d) $\lambda =0.5$ for the trivial phase.

Download Full Size | PDF

Time evolution of several amplitudes $|\psi _{n}|$ is plotted in Fig. 14. The results of the topological phase with NNN couplings in Fig. 14(a) are almost the same as the results without NNN couplings in Fig. 6(a), which implies the robustness of the higher-order corner states by regarding the NNN couplings as the perturbation. However, different from the trivial case in Fig. 6(b), the amplitudes $|\psi _{n}|$ do not decay to zero. Moreover, the amplitude $|\psi _{2}|$ of site near the corner site $A$ is increased by the NNN couplings. A possible reason is the NNN couplings enhance the interaction of the sites in a same unit cell, making the model is more like a 2D SSH model in the tetramer limit.

 figure: Fig. 14.

Fig. 14. Time evolution of several amplitudes $|\psi _{n}|$ of the 2D SSH model with NNN couplings. The parameters are set as $\kappa =1$, $\gamma =1/2$, $\xi =2$, $\eta =1$, (a) $\lambda =-0.5$ for the topological phase and (b) $\lambda =0.5$ for the trivial phase.

Download Full Size | PDF

Next, we compare the result of Fig. 14(b) with that of the tetramer with NNN couplings. In the tetramer limit, the equations of the tetramer are

$$\mathrm{i}\frac{\mathrm{d}\psi_{1}}{\mathrm{d}t}=\kappa_{A}(\psi_{2}+\psi_{3}+\psi_{4})-\mathrm{i}\gamma\left(1-\frac{\xi}{1+|\psi_{1}|^{2}/\eta}\right)\psi_{1},$$
$$\mathrm{i}\frac{\mathrm{d}\psi_{2}}{\mathrm{d}t}=\kappa_{A}(\psi_{1}+\psi_{3}+\psi_{4})-\mathrm{i}\gamma\psi_{2},$$
$$\mathrm{i}\frac{\mathrm{d}\psi_{3}}{\mathrm{d}t}=\kappa_{A}(\psi_{1}+\psi_{2}+\psi_{4})-\mathrm{i}\gamma\psi_{3},$$
$$\mathrm{i}\frac{\mathrm{d}\psi_{4}}{\mathrm{d}t}=\kappa_{A}(\psi_{1}+\psi_{2}+\psi_{3})-\mathrm{i}\gamma\psi_{4}.$$

The numerically solutions of Eqs. (13)–(16) with the initial condition $\psi _{1}(0)=1$ are plotted in Fig. 15(a). Amplitude $|\psi _{1}|$ approaches a saturated value, as well as the amplitude $|\psi _{2}|$. Due to the symmetry of the tetramer with NNN couplings, we have $|\psi _{2}|=|\psi _{3}|=|\psi _{4}|$. Comparing the solutions of amplitude $|\psi _{1}|$ and $|\psi _{2}|$ in Fig. 15(a) with that in Fig. 14(b) for the trivial phase, we find the results of the two cases are very close, indicating the NNN couplings make the tetramerization play a leading role in the dynamics of 2D SSH model in the trivial phase when taking $\chi =\kappa _{A}$. Equations of the tetramer by removing the losses on sites $2$, $3$ and $4$ are written as

$$\mathrm{i}\frac{\mathrm{d}\psi_{1}}{\mathrm{d}t}=\kappa_{A}(\psi_{2}+\psi_{3}+\psi_{4})-\mathrm{i}\gamma\left(1-\frac{\xi}{1+|\psi_{1}|^{2}/\eta}\right)\psi_{1},$$
$$\mathrm{i}\frac{\mathrm{d}\psi_{2}}{\mathrm{d}t}=\kappa_{A}(\psi_{1}+\psi_{3}+\psi_{4}),$$
$$\mathrm{i}\frac{\mathrm{d}\psi_{3}}{\mathrm{d}t}=\kappa_{A}(\psi_{1}+\psi_{2}+\psi_{4}),$$
$$\mathrm{i}\frac{\mathrm{d}\psi_{4}}{\mathrm{d}t}=\kappa_{A}(\psi_{1}+\psi_{2}+\psi_{3}).$$

Then, the stationary solutions of Eqs. (17)–(20) are given by

$$\begin{aligned} \psi_{1}&=\sqrt{\frac{\eta(3i\kappa_{A}+2\gamma)[3i\kappa_{A}+2\gamma(\xi-1)]}{4\gamma^{2}+9\kappa_{A}^{2}}},\\ \psi_{2}&=\psi_{3}=\psi_{4}=\mathrm{const}, \end{aligned}$$
and the numerical solutions are shown in Fig. 15(b).

 figure: Fig. 15.

Fig. 15. Time evolution the amplitudes $|\psi _{n}|$ of the tetramer with NNN couplings. Tetramer (a) with and (b) without losses on sites $2$, $3$ and $4$. The parameters are set as $\kappa _{A}=1$, $\kappa _{B}=0$, $\gamma =1/2$, $\xi =2$ and $\eta =1$.

Download Full Size | PDF

4. Conclusions

In conclusion, we have investigated a topological laser system in coupled-ring array made of active resonators corresponding to the 2D SSH model. Through the quench dynamics, we find the nonlinear non-Hermitian 2D SSH model is in the topological phase for $-1<\lambda <0$ and in the trivial phase for $0<\lambda <1$, as same as the normal 2D SSH model. When a pulse is given to any one site of the lattice in the topological phase, the amplitudes of all sites belonging to the topological corner states will reach a finite value after a delay time, meaning that these sites begin to emit stable laser light. The time delay can be controlled by tuning the hopping amplitude $\kappa$. In the trivial phase, the amplitude of the stimulated site decays rapidly, and no laser light is emitted by other sites. However, the stationary solutions of the amplitudes of the four sites in the tetramer limit are nonzero.

Furthermore, for a 2D coupled waveguide array [22] or a coupled-ring array [32,33,4850], the coupling between two sites always depends on the distance between them. So, in the strong-coupling unit cells, where sites are close to each other, the NNN couplings in a coupled-ring array may not be ignored. We introduce NNN couplings into the strong-coupling unit cells of 2D SSH model and discuss the dynamics of the topological laser system based on it. In the topological phase, it is found that the corner sites begin to emit stable laser light when any one site is stimulated, which is similar to the results of the 2D SSH model without NNN couplings. While in the trivial phase, the corner sites can also emit stable laser light although weaker than that in the topological phase. The reason is that tetramerization plays a leading role in the dynamics of the laser system due to the existence of the NNN couplings in the strong-coupling unit cells.

Funding

National Key Research and Development Program of China (2021YFA1400602); National Natural Science Foundation of China (12174288, 12204352, 12274326).

Acknowledgments

This work is supported by the National Natural Science Foundation of China (Grants No. 12174288, No. 12204352 and No. 12274326) and the National Key R & D Program of China (Grant No. 2021YFA1400602).

Disclosures

The authors declare no conflicts of interest.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

References

1. M. Z. Hasan and C. L. Kane, “Colloquium: Topological insulators,” Rev. Mod. Phys. 82(4), 3045–3067 (2010). [CrossRef]  

2. X.-L. Qi and S.-C. Zhang, “Topological insulators and superconductors,” Rev. Mod. Phys. 83(4), 1057–1110 (2011). [CrossRef]  

3. K. v. Klitzing, G. Dorda, and M. Pepper, “New method for high-accuracy determination of the fine-structure monstant based on quantized Hall resistance,” Phys. Rev. Lett. 45(6), 494–497 (1980). [CrossRef]  

4. F. D. M. Haldane, “Model for a quantum Hall effect without Landau Levels: Condensed-matter realization of the ‘parity anomaly’,” Phys. Rev. Lett. 61(18), 2015–2018 (1988). [CrossRef]  

5. C. L. Kane and E. J. Mele, “Quantum spin Hall effect in graphene,” Phys. Rev. Lett. 95(22), 226801 (2005). [CrossRef]  

6. B. A. Bernevig, T. L. Hughes, and S.-C. Zhang, “Quantum spin Hall effect and topological phase transition in HgTe quantum wells,” Science 314(5806), 1757–1761 (2006). [CrossRef]  

7. L. Fu, “Topological crystalline insulators,” Phys. Rev. Lett. 106(10), 106802 (2011). [CrossRef]  

8. F. Zhang, C. L. Kane, and E. J. Mele, “Surface state magnetization and chiral edge states on topological insulators,” Phys. Rev. Lett. 110(4), 046404 (2013). [CrossRef]  

9. W. A. Benalcazar, B. A. Bernevig, and T. L. Hughes, “Quantized electric multipole insulators,” Science 357(6346), 61–66 (2017). [CrossRef]  

10. Y. Peng, Y. Bao, and F. von Oppen, “Boundary Green functions of topological insulators and superconductors,” Phys. Rev. B 95(23), 235143 (2017). [CrossRef]  

11. J. Langbehn, Y. Peng, L. Trifunovic, F. von Oppen, and P. W. Brouwer, “Reflection-symmetric second-order topological insulators and superconductors,” Phys. Rev. Lett. 119(24), 246401 (2017). [CrossRef]  

12. Z. Song, Z. Fang, and C. Fang, “(d-2)-dimensional edge states of rotation symmetry protected topological states,” Phys. Rev. Lett. 119(24), 246402 (2017). [CrossRef]  

13. W. A. Benalcazar, B. A. Bernevig, and T. L. Hughes, “Electric multipole moments, topological multipole moment pumping, and chiral hinge states in crystalline insulators,” Phys. Rev. B 96(24), 245115 (2017). [CrossRef]  

14. F. Schindler, A. M. Cook, M. G. Vergniory, Z. Wang, S. S. P. Parkin, B. A. Bernevig, and T. Neupert, “Higher-order topological insulators,” Sci. Adv. 4(6), eaat0346 (2018). [CrossRef]  

15. M. Lin and T. L. Hughes, “Topological quadrupolar semimetals,” Phys. Rev. B 98(24), 241103 (2018). [CrossRef]  

16. M. Ezawa, “Magnetic second-order topological insulators and semimetals,” Phys. Rev. B 97(15), 155305 (2018). [CrossRef]  

17. M. Ezawa, “Higher-order topological insulators and semimetals on the breathing kagome and pyrochlore lattices,” Phys. Rev. Lett. 120(2), 026801 (2018). [CrossRef]  

18. M. Ezawa, “Topological switch between second-order topological insulators and topological crystalline insulators,” Phys. Rev. Lett. 121(11), 116801 (2018). [CrossRef]  

19. C. H. Lee, S. Imhof, C. Berger, F. Bayer, J. Brehm, L. W. Molenkamp, T. Kiessling, and R. Thomale, “Topolectrical circuits,” Commun. Phys. 1(1), 39 (2018). [CrossRef]  

20. S. Imhof, C. Berger, F. Bayer, J. Brehm, L. Molenkamp, T. Kiessling, F. Schindler, C. H. Lee, M. Greiter, T. Neupert, and R. Thomale, “Topolectrical-circuit realization of topological corner modes,” Nat. Phys. 14(9), 925–929 (2018). [CrossRef]  

21. M. Ezawa, “Higher-order topological electric circuits and topological corner resonance on the breathing kagome and pyrochlore lattices,” Phys. Rev. B 98(20), 201402 (2018). [CrossRef]  

22. H. Yang, Z.-X. Li, Y. Liu, Y. Cao, and P. Yan, “Observation of symmetry-protected zero modes in topolectrical circuits,” Phys. Rev. Res. 2(2), 022028 (2020). [CrossRef]  

23. J. Wu, X. Huang, J. Lu, Y. Wu, W. Deng, F. Li, and Z. Liu, “Observation of corner states in second-order topological electric circuits,” Phys. Rev. B 102(10), 104109 (2020). [CrossRef]  

24. M. S.-Garcia, V. Peri, R. Susstrunk, O. R. Bilal, T. Larsen, L. G. Villanueva, and S. D. Huber, “Observation of a phononic quadrupole topological insulator,” Nature (London) 555(7696), 342–345 (2018). [CrossRef]  

25. Z. Liu, “Trapping sound at corners,” Nat. Mater. 18(2), 98–99 (2019). [CrossRef]  

26. H. Xue, Y. Yang, F. Gao, Y. Chong, and B. Zhang, “Acoustic higher-order topological insulator on a kagome lattice,” Nat. Mater. 18(2), 108–112 (2019). [CrossRef]  

27. B. Y. Xie, H. F. Wang, H.-X. Wang, X. Y. Zhu, J.-H. Jiang, M. H. Lu, and Y. F. Chen, “Second-order photonic topological insulator with corner states,” Phys. Rev. B 98(20), 205147 (2018). [CrossRef]  

28. A. Cerjan, M. Jürgensen, W. A. Benalcazar, S. Mukherjee, and M. C. Rechtsman, “Observation of a higher-order topological bound state in the continuum,” Phys. Rev. Lett. 125(21), 213901 (2020). [CrossRef]  

29. M. Li, Y. Wang, T. Sang, H. Chu, Y. Lai, and G. Yang, “Experimental observation of multiple edge and corner states in photonic slabs heterostructures,” Photonics Res. 10(1), 197–204 (2022). [CrossRef]  

30. P. St-Jean, V. Goblot, E. Galopin, A. Lemaitre, T. Ozawa, L. Le Gratiet, I. Sagnes, J. Bloch, and A. Amo, “Lasing in topological edge states of a one-dimensional lattice,” Nat. Photonics 11(10), 651–656 (2017). [CrossRef]  

31. Y. Ota, R. Katsumi, K. Watanabe, S. Iwamoto, and Y. Arakawa, “Topological photonic crystal nanocavity laser,” Commun. Phys. 1(1), 86 (2018). [CrossRef]  

32. G. Harari, M. A. Bandres, Y. Lumer, M. C. Rechtsman, Y. D. Chong, M. Khajavikhan, D. N. Christodoulides, and M. Segev, “Topological insulator laser: Theory,” Science 359(6381), eaar4003 (2018). [CrossRef]  

33. M. A. Bandres, S. Wittek, G. Harari, M. Parto, J. Ren, M. Segev, D. N. Christodoulides, and M. Khajavikhan, “Topological insulator laser: Experiments,” Science 359(6381), eaar4005 (2018). [CrossRef]  

34. M. Parto, S. Wittek, H. Hodaei, G. Harari, M. A. Bandres, J. Ren, M. C. Rechtsman, M. Segev, D. N. Christodoulides, and M. Khajavikhan, “Edge-Mode Lasing in 1D Topological Active Arrays,” Phys. Rev. Lett. 120(11), 113901 (2018). [CrossRef]  

35. H. Zhao, P. Miao, M. H. Teimourpour, S. Malzard, R. El-Ganainy, H. Schomerus, and L. Feng, “Topological hybrid silicon microlasers,” Nat. Commun. 9(1), 981 (2018). [CrossRef]  

36. S. Malzard and H. Schomerus, “Nonlinear mode competition and symmetry-protected power oscillations in topological lasers,” New J. Phys. 20(6), 063044 (2018). [CrossRef]  

37. S. Malzard, E. Cancellieri, and H. Schomerus, “Topological dynamics and excitations in lasers and condensates with saturable gain or loss,” Opt. Express 26(17), 22506–22518 (2018). [CrossRef]  

38. Y. V. Kartashov and D. V. Skryabin, “Two-dimensional topological polariton laser,” Phys. Rev. Lett. 122(8), 083902 (2019). [CrossRef]  

39. M. Seclì, M. Capone, and I. Carusotto, “Theory of chiral edge state lasing in a two-dimensional topological system,” Phys. Rev. Res. 1(3), 033148 (2019). [CrossRef]  

40. S. Gündoğdu, J. Thingna, and D. Leykam, “Edge mode bifurcations of two-dimensional topological lasers,” Opt. Lett. 45(13), 3673–3676 (2020). [CrossRef]  

41. W. Noh, H. Nasari, H.-M. Kim, Q. Le-Van, Z. Jia, C.-H. Huang, and B. Kanté, “Experimental demonstration of single-mode topological valley-Hall lasing at telecommunication wavelength controlled by the degree of asymmetry,” Opt. Lett. 45(15), 4108–4111 (2020). [CrossRef]  

42. Y. Zeng, U. Chattopadhyay, B. Zhu, B. Qiang, J. Li, Y. Jin, L. Li, A. G. Davies, E. H. Linfield, B. Zhang, Y. D. Chong, and Q. J.Wang, “Electrically pumped topological laser with valley edge modes,” Nature (London) 578(7794), 246–250 (2020). [CrossRef]  

43. J. H. Choi, W. E. Hayenga, Y. G. Liu, M. Parto, B. Bahari, D. N. Christodoulides, and M. Khajavikhan, “Room temperature electrically pumped topological insulator lasers,” Nat. Commun. 12(1), 3434 (2021). [CrossRef]  

44. Z. Yang, E. Lustig, G. Harari, Y. Plotnik, Y. Lumer, M. A. Bandres, and M. Segev, “Mode-locked Topological Insulator Laser Utilizing Synthetic Dimensions,” Phys. Rev. X 10(1), 011059 (2020). [CrossRef]  

45. N. Ishida, Y. Ota, W. Lin, T. Byrnes, Y. Arakawa, and S. Iwamoto, “A large-scale single-mode array laser based on a topological edge mode,” Nanophotonics 11(9), 2169–2181 (2022). [CrossRef]  

46. M. Ezawa, “Nonlinear non-Hermitian higher-order topological laser,” Phys. Rev. Res. 4(1), 013195 (2022). [CrossRef]  

47. W. X. Teo, W. Zhu, and J. Gong, “Tunable two-dimensional laser arrays with zero-phase locking,” Phys. Rev. B 105(20), L201402 (2022). [CrossRef]  

48. M. Hafezi, E. A. Demler, M. D. Lukin, and J. M. Taylor, “Robust optical delay lines with topological protection,” Nat. Phys. 7(11), 907–912 (2011). [CrossRef]  

49. M. Hafezi, S. Mittal, J. Fan, A. Migdall, and J. M. Taylor, “Imaging topological edge states in silicon photonics,” Nat. Photonics 7(12), 1001–1005 (2013). [CrossRef]  

50. S. Mittal, V. V. Orre, D. Leykam, Y. D. Chong, and M. Hafezi, “Photonic anomalous quantum Hall effect,” Phys. Rev. Lett. 123(4), 043201 (2019). [CrossRef]  

51. F. Liu and K. Wakabayashi, “Novel topological phase with a zero Berry curvature,” Phys. Rev. Lett. 118(7), 076803 (2017). [CrossRef]  

52. W. P. Su, J. R. Schrieffer, and A. J. Heeger, “Solitons in polyacetylene,” Phys. Rev. Lett. 42(25), 1698–1701 (1979). [CrossRef]  

53. W. A. Benalcazar and A. Cerjan, “Bound states in the continuum of higher-order topological insulators,” Phys. Rev. B 101(16), 161116 (2020). [CrossRef]  

54. P. A. Ivanova, N. A. Olekhno, V. I. Kachin, D. V. Zhirihin, P. S. Seregin, and M. A. Gorlach, “Realizing topological corner states in two-dimensional Su-Schrieffer-Heeger model with next-nearest neighbor couplings,” J. Phys.: Conf. Ser. 1695(1), 012142 (2020). [CrossRef]  

55. N. A. Olekhno, A. D. Rozenblit, V. I. Kachin, A. A. Dmitriev, O. I. Burmistrov, P. S. Seregin, D. V. Zhirihin, and M. A. Gorlach, “Experimental realization of topological corner states in long-range-coupled electrical circuits,” Phys. Rev. B 105(8), L081107 (2022). [CrossRef]  

56. M. Ezawa, “Nonlinearity-induced transition in nonlinear Su-Schrieffer-Heeger model and a nonlinear higher-order topological system,” Phys. Rev. B 104(23), 235420 (2021). [CrossRef]  

57. M. Ezawa, “Quench dynamics and bulk-edge correspondence in nonlinear mechanical systems,” J. Phys. Soc. Jpn. 90(11), 114605 (2021). [CrossRef]  

58. M. Ezawa, “Topological edge states and bulk-edge correspondence in dimerized Toda lattice,” J. Phys. Soc. Jpn. 91(2), 024703 (2022). [CrossRef]  

59. M. Ezawa, “Nonlinear topological phase transitions in the dimerized sine-Gordon model,” Phys. Rev. B 105(16), 165418 (2022). [CrossRef]  

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (15)

Fig. 1.
Fig. 1. Schematics of a finite lattice of 2D SSH model. Each unit cell contains four sites, and $\kappa _A$ ($\kappa _B$) is the intracell (intercell) hopping.
Fig. 2.
Fig. 2. (a) The energy spectrum of a finite lattice of 2D SSH model with the size $L=6$. (b) Density of a corner state corresponding to a zero-energy mode in the topological phase with $\lambda =-0.5$. (c),(d) Density of two bulk states in the trivial phase with $\lambda =0.5$.
Fig. 3.
Fig. 3. The saturated amplitude $|\psi _{1}|$ at the corner $A$ after enough time. The parameters are set as $\kappa =1$, $\xi =2$ and $\eta =1$.
Fig. 4.
Fig. 4. Time evolution of the amplitude $|\psi _{n}|$ of the 2D SSH model. The parameters are set as $\kappa =1$, $\gamma =1/2$, $\xi =2$, $\eta =1$, (a1)-(a5) $\lambda =-0.5$ for the topological phase and (b1)-(b5) $\lambda =0.5$ for the trivial phase.
Fig. 5.
Fig. 5. Density plot of time evolution of the amplitude $|\psi _{n}|$ of the side $AB$. The sites on the side are marked from $1$ to $12$. (a),(b) The initial condition $\psi _{1}(0)=1$, and (c),(d) the initial condition is $\psi _{6}(0)=1$. The parameters are set as $\kappa =1$, $\gamma =1/2$, $\xi =2$, $\eta =1$, (a),(c) $\lambda =-0.5$ for the topological phase and (b),(d) $\lambda =0.5$ for the trivial phase.
Fig. 6.
Fig. 6. Time evolution of several amplitudes $|\psi _{n}|$. The parameters are set as $\kappa =1$, $\gamma =1/2$, $\xi =2$, $\eta =1$, (a) $\lambda =-0.5$ for the topological phase and (b) $\lambda =0.5$ for the trivial phase.
Fig. 7.
Fig. 7. Delay time of the saturated amplitudes (a) $|\psi _{B}|$ and (b) $|\psi _{D}|$ with $\kappa =1/4, 1/2, 1,\ldots,2$. The parameters are set as $\kappa =1$, $\gamma =1/2$, $\xi =2$, $\eta =1$, and $\lambda =-0.5$
Fig. 8.
Fig. 8. Time evolution of the amplitudes $|\psi _{n}|$ of the tetramer in 2D SSH model (a),(b) with and (c),(d) without losses on sites $2$, $3$ and $4$. The parameters are set as $\kappa _{A}=1$, $\kappa _{B}=0$, $\gamma =1/2$, $\xi =2$ and $\eta =1$.
Fig. 9.
Fig. 9. 2D SSH model with NNN couplings $\chi$ in the strong-coupling unit cells. (a) The strong-coupling unit cells are cells with hopping $\kappa _{B}$, and (b) the strong-coupling unit cells are cells with hopping $\kappa _{A}$.
Fig. 10.
Fig. 10. Energy spectrum of a finite lattice of 2D SSH model (a) with NNN couplings $\chi$ in cells with coupling $\kappa _{A}$ and (b) with NNN couplings $\chi$ in cells with coupling $\kappa _{B}$. Energy eigenvalues of 2D SSH model of size $L=6$ (c) without NNN couplings and (d) with NNN couplings in the topological phase. Red points are the zero-energy eigenvalues of the topological corner states. The parameters in (c),(d) are set as $\kappa =1$ and $\lambda =-0.5$.
Fig. 11.
Fig. 11. (a) The density of a corner state corresponding to one of the zero-energy eigenvalues in Fig. 10(d). (b) The density of a bulk state of the 2D SSH model with NNN couplings in the trivial phase ($\lambda =0.5$).
Fig. 12.
Fig. 12. Time evolution of the amplitude $|\psi _{n}|$ of the 2D SSH model with NNN couplings. The parameters are set as $\kappa =1$, $\gamma =1/2$, $\xi =2$, $\eta =1$, (a1)-(a5) $\lambda =-0.5$ for the topological phase and (b1)-(b5) $\lambda =0.5$ for the trivial phase.
Fig. 13.
Fig. 13. Density plot of time evolution of the amplitude $|\psi _{n}|$ of the side $AB$ of the 2D SSH model with NNN couplings. The sites on the side are marked from $1$ to $12$. (a),(b) The initial condition $\psi _{1}(0)=1$, and (c),(d) the initial condition is $\psi _{6}(0)=1$. The parameters are set as $\kappa =1$, $\gamma =1/2$, $\xi =2$, $\eta =1$, (a),(c) $\lambda =-0.5$ for the topological phase and (b),(d) $\lambda =0.5$ for the trivial phase.
Fig. 14.
Fig. 14. Time evolution of several amplitudes $|\psi _{n}|$ of the 2D SSH model with NNN couplings. The parameters are set as $\kappa =1$, $\gamma =1/2$, $\xi =2$, $\eta =1$, (a) $\lambda =-0.5$ for the topological phase and (b) $\lambda =0.5$ for the trivial phase.
Fig. 15.
Fig. 15. Time evolution the amplitudes $|\psi _{n}|$ of the tetramer with NNN couplings. Tetramer (a) with and (b) without losses on sites $2$, $3$ and $4$. The parameters are set as $\kappa _{A}=1$, $\kappa _{B}=0$, $\gamma =1/2$, $\xi =2$ and $\eta =1$.

Equations (21)

Equations on this page are rendered with MathJax. Learn more.

H = x , y [ κ ( 1 ( 1 ) x λ ) a x , y a x + 1 , y + κ ( 1 ( 1 ) y λ ) a x , y a x , y + 1 ] + H . c . ,
i d ψ n d t = m = 1 4 L 2 H n m ψ m i γ ( 1 ξ P n 1 + | ψ n | 2 / η ) ψ n ,
P n = n c δ n , n c ,
i d ψ 1 d t = κ A ( ψ 2 + ψ 3 ) i γ ( 1 ξ 1 + | ψ 1 | 2 / η ) ψ 1 ,
i d ψ 2 d t = κ A ( ψ 1 + ψ 4 ) i γ ψ 2 ,
i d ψ 3 d t = κ A ( ψ 1 + ψ 4 ) i γ ψ 3 ,
i d ψ 4 d t = κ A ( ψ 2 + ψ 3 ) i γ ψ 4 .
i d ψ 1 d t = κ A ( ψ 2 + ψ 3 ) i γ ( 1 ξ 1 + | ψ 1 | 2 / η ) ψ 1 ,
i d ψ 2 d t = κ A ( ψ 1 + ψ 4 ) ,
i d ψ 3 d t = κ A ( ψ 1 + ψ 4 ) ,
i d ψ 4 d t = κ A ( ψ 2 + ψ 3 ) .
ψ 1 = ψ 4 = η ( ξ 1 ) , ψ 2 = ψ 3 = c o n s t .
i d ψ 1 d t = κ A ( ψ 2 + ψ 3 + ψ 4 ) i γ ( 1 ξ 1 + | ψ 1 | 2 / η ) ψ 1 ,
i d ψ 2 d t = κ A ( ψ 1 + ψ 3 + ψ 4 ) i γ ψ 2 ,
i d ψ 3 d t = κ A ( ψ 1 + ψ 2 + ψ 4 ) i γ ψ 3 ,
i d ψ 4 d t = κ A ( ψ 1 + ψ 2 + ψ 3 ) i γ ψ 4 .
i d ψ 1 d t = κ A ( ψ 2 + ψ 3 + ψ 4 ) i γ ( 1 ξ 1 + | ψ 1 | 2 / η ) ψ 1 ,
i d ψ 2 d t = κ A ( ψ 1 + ψ 3 + ψ 4 ) ,
i d ψ 3 d t = κ A ( ψ 1 + ψ 2 + ψ 4 ) ,
i d ψ 4 d t = κ A ( ψ 1 + ψ 2 + ψ 3 ) .
ψ 1 = η ( 3 i κ A + 2 γ ) [ 3 i κ A + 2 γ ( ξ 1 ) ] 4 γ 2 + 9 κ A 2 , ψ 2 = ψ 3 = ψ 4 = c o n s t ,
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.