1 Introduction
In recent years, the study of Hermitian systems has gained increasing attention, as it holds significant implications for topology research. This includes categorizing and defining topological phases [
1−
5], as well as its application in various fields such as superconductors [
6−
8], cold atoms [
9−
14], and solid states [
15,
16]. However, the non-Hermitian system is more closely related to real systems and widely exists because losses or gains may be inevitable, which brings challenges to basic topological classifications [
16−
18] and corresponding boundary correspondence [
20−
29]. Nevertheless, interesting non-Hermitian topological phenomena have also been observed, such as the occurrence of amplified or attenuated modes in the valley states [
30], the degenerated Weyl evolves into the Weyl exceptional ring [
31], topological interface generated by gain and loss [
32,
33], or non-Hermitian skin effect protected by spectrum winding numbers [
34−
38]. By allowing the Hamiltonian to be non-Hermitian, we have expanded the different symmetries protected phases [
17,
18]. It is commonly believed that dissipation in topological systems is detrimental, but in reality, non-Hermitianity can play a crucial role, not just as a source of disturbance, such as for robust target boundary modes [
39].
Topological modes are usually measured by spectral measurements, but topological phases may be embedded in boundary states in the bulk spectrum, which makes it impossible to accurately distinguish between topological phases, especially in higher-order of topology [
40,
41]. Recently, research on topological insulator crystals has revealed that crystal defects can serve as effective probes for exploring the topology of crystals beyond the edge, a finding that has been experimentally confirmed by [
42−
48]. The local density of states (LDOSs) of acoustic systems can be obtained through direct measurement of acoustic metamaterials [
47,
48]. This is advantageous for the classification of topological acoustic crystals, but the exploration of the possibility of non-Hermitian LDOSs remains an open question.
In this work, we present a simulation demonstration of the topological monomodes corresponding to non-Hermitian acoustic system, which is determined by observing the fractional mode charge of the system. We propose a novel coupled acoustic cavity structure that introduces non-Hermitian losses in the coupled cavity, thereby enabling the realization of topological monomode, even in high-dimensional structures. We have also observed that the topological properties of this structure can be analyzed through the fractional mode charge. We have verified the existence of this topological effect both theoretically and through simulations, and we have been able to observe the topological bandgap and its associated edge and corner states. Moreover, we compared the bandgap characteristics under different loss configurations and found that even disturbance could lead to the appearance of topological states under specific conditions. This work has significant implications for understanding the topological properties of non-Hermitian acoustic crystals and exploring applications of high-order non-Hermitian topological states.
2 Model and simulation
2.1 1D SSH model
We first consider a typical example is the 1D SSH model, where units are composed of two closely coupled cavities and waveguides. As illustrated in Fig.1(a), the system can be tuned to modify its energy band structure and band gaps by controlling the sizes of internal coupling
t1 and external coupling
t2, which is accomplished through dimerization coupling [
49,
50]. Then we add an additional loss is artificially introduced in one of the resonators, topological monomode emerges in Fig.1(b), which depends on the configuration of the losses. In tight-binding model, the Hamiltonian of one-dimensional finite chain is written as:
Here, represents the on-site potential, the interunit-cell coupling strength is stronger than the intraunit-cell coupling, which can be represented by t2 = 1.33∗t1. The loss can be simulated by the imaginary part of the dispersion.
2.2 Simulation calculation of local density of states
In simulations, we employed the Pressure-Sound module of COMSOL Multiphysics. For all the simulations, the density of the background medium air was set to 1.22 kg/m3, and the real part of the sound speed in air was set to 340 m/s. By adjusting the sound speed, we fitted the simulated spectrum of the resonator with only background loss (c = 340 + 2.2i m/s) and the simulated spectrum of the resonator with additional loss (c = 340 + 85.0i m/s). For the calculation of the band structure, we adopted periodic boundary conditions at the outer boundary, while the other boundaries were considered as hard boundaries. In the calculation of the finite-element eigenfrequency [Fig.1(a) and (b)], all boundaries were assumed to be hard boundaries.
We calculated the eigenstates of the 1D chain, and the 12 eigenstates correspond to the number of cavities. The edge states are located in the bandgap, and the frequencies of the edge states on the left and right boundaries are approximately the same in Fig.1(c). After adding loss, the frequency of the right edge state is shifted in Fig.1(d), which is because the loss breaks the symmetry of the one-dimensional chain. When the loss tends to be finite, one of these two zero modes acquires an energy with a negative imaginary part, while the other remains at zero energy. This phenomenon arises from the fact that we have introduced loss to only one of the two sublattices. Each edge state exists on only one sublattice, so they exhibit different energy properties. In previous studies, we have known that local state density can characterize topological states in Hermitian systems, especially in defect structures [
42−
44]. Then, we calculate the local state density, with calculation details referring to Ref. [
47]. The acoustic LDOS is obtained from
where is the mass density of air, is the acoustic particle velocity on the surface of monopole source, and is the surface area of the source. The velocity stands for the oscillating velocity of the air flows. The fractional mode charge can be obtained by integrating the LDOS above the edge state frequency in the band gap.
3 Results and discussion
3.1 The fractional mode charge of 1D topological monomode
Theoretical and simulation results show that although the bandgap is small, and the boundary states and body modes have hybridization, the fractional charge modes on the left and right edges are approximately 1, while the fractional charge modes in the middle are fractional in Fig.1(e). This feature allows us to distinguish boundary states through local state density in Fig.1(g). We also simulated topological monomode in non-Hermitian systems in Fig.1(f). By adding loss in the cavity, a monomode can be formed, and this phenomenon can be observed through local state density as well, which is similar to Hermitian systems. Specifically, the local state density of the left boundary state changes from outward diffusion in Fig.1(h), which is due to the damage of the edge modes by loss. We can also find that the fractional charge mode on the right side is significantly greater than the left side, indicating that we detect the topological monomodes. This may bring new applications to robust boundary states.
3.2 The fractional mode charge of 2D topological monomode
To further investigate the impact of loss on topology, we extend the model to two dimensions. We first solve the conventional two-dimensional (2D) SSH model and find that the eigenstate is localized at the four corners in Fig.2(a). Then, we try different operations, introducing loss along the x and y directions of each cavity, to reduce the symmetry and form diagonal eigenmodes in Fig.2(b). By further increasing the amount of loss, we break the symmetry even more, resulting in a eigenmode in Fig.2(c), which is referred to as a monomode. The existence of monomode in non-Hermitian systems is theoretically and simulatively proven through the introduction of loss in the selected crystal structure and near the corresponding edges, where one topological edge mode will decay over time. We have already simulatively realized these monomode in acoustic coupled systems. This result shows that under appropriate conditions, we can effectively control and adjust the symmetry and the number of corner modes of the acoustic coupled systems.
We are more interested in characterizing the local state density as a means of addressing this issue. When the corner and body modes are hybridized, the boundary state can be entirely located within the bandgap, making it difficult to distinguish it from the bulk state using the bulk-boundary correspondence principle. Therefore, we also calculate the local state density of each cavity in Fig.2(d). The local density of state in the corner state is significantly distinct from that in the bulk state. As shown in Fig.2(e), when loss is introduced into the system, the local state density of the double-mode corner state changes, with the peak transforming from a single to a double-peaked state. This result is attributed to the disappearance of symmetry, as evidenced by both theoretical and simulation findings in Fig.2(g) and (f). The topological monomodes induced by the introduced losses through acoustic local density of states simulated are directly observed, leading to the appearance of edge states and intermediate gap corner states, which are typical features of higher-order topology [
33].
3.3 The fractional mode charge of 2D robust topological monomode
Previous research indicates that disorder typically has a significant impact on altering the energy band structure within quantum systems [
51,
52]. We hope that topological monomode can also maintain their integrity without being perturbed. As a comparison, we conduct simulations to investigate the effects of random perturbations on loss configurations. To evaluate the robustness of this structure, we have introduced random losses to each cavity, which is consistent with the real-life situation, as the experimental process inevitably introduces losses. The disturbance intensity is
E00.25
rand(−1,1)
i, and the simulation results demonstrate that topological monomode exhibit strong robustness to perturbations in Fig.3(a) and (b). Subsequently, fractional mode charge calculations were carried out at all positions on the topological crystal, and it was found that while the fractional mode charge changed in Fig.3(c), the corner state characteristics remained evident, further demonstrating the robustness of the structure to perturbations. The simulated fractional mode charge is presented in Fig.3(d) and (e). Since each disturbance strength is random, the simulation and simulation may not be consistent. We found that although there are some fluctuations in the peaks, the overall intensity remains basically unchanged. This indicates that although there are disturbances, the overall system remains stable, which is consistent with the inevitable introduction of losses or gains in real life.
3.4 The fractional mode charge of 3D topological monomode
Most crystals have a natural three-dimensional structure. To further demonstrate the practicality of monomode, we have extended it to a three-dimensional structure, as shown in Fig.4(a). Similar to the previous two-dimensional structure, we display the topological state of the loss structure within the three-dimensional structure. The three-dimensional structure exhibits a surface state that complicates the band structure compared to the two-dimensional structure. In addition to incorporating loss into the cavity, we have observed diagonal-symmetric corner modes in Fig.4(a) and (b), making it easier to design practical acoustic structures. We have also simulated the local density of states in the three-dimensional structure in Fig.4(c), similar to the two-dimensional structure, where the corner state displays a single peak; however, due to a symmetry break, the corner state becomes a double peak in Fig.4(d) and (e). By integrating the local density of states in the corresponding frequency domain, we obtained fractional mode charge, which indicates that the corner charge is much higher than that of the boundary state and bulk state. Therefore, we can verify the topological monomode under non-Hermitian systems through the local density of state.
4 Conclusions
In summary, we have observed topological monomodes arising from the addition of losses in non-Hermitian acoustic coupled systems. We have verified the acoustic topological insulators under different loss configurations, starting from low-dimensional structures by introducing losses in selected cavities with adequate proximity to the corresponding edges. We first achieved topological monomode in one-dimensional structures and then extended it to high-dimensional structures. More importantly, we found that topological monomodes of non-Hermitian acoustic coupled systems can be characterized by fractional mode charges, even in higher-order structures. Lastly, we evaluated the robustness of topological monomodes to disturbances, and the experimental results showed that these topological monomodes exhibited robustness and did not hybridize, pointing to huge potential applications, as non-Hermiticity can be externally modulated, reconfigurable, and active. Additionally, the tunability of losses opens up great potential for future reconfigurable devices [
53,
54]. These results may be applied to explore higher-order non-Hermitian models and provide a flexible and useful scheme for characterizing topological states under non-Hermiticity.