Next Article in Journal
Experimental Investigation on the Structural Performance of Single Span Hollow Core Slab under Successive Impact Loading
Previous Article in Journal
A Kinetic Model Considering Catalyst Deactivation for Methanol-to-Dimethyl Ether on a Biomass-Derived Zr/P-Carbon Catalyst
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Localization Properties of a Quasiperiodic Ladder under Physical Gain and Loss: Tuning of Critical Points, Mixed-Phase Zone and Mobility Edge

1
Physics and Applied Mathematics Unit, Indian Statistical Institute, 203 Barrackpore Trunk Road, Kolkata 700 108, India
2
Departamento de Física, FACI, Universidad de Tarapacá, Casilla 7D, Arica 1000000, Chile
3
Grupo de Física de Materiales, Universidad Pedagógica y Tecnológica de Colombia, Tunja 150003, Colombia
4
Laboratorio de Química Teórica y Computacional, Grupo de Investigación Química-Física Molecular y Modelamiento Computacional (QUIMOL), Facultad de Ciencias, Universidad Pedagógica y Tecnológica de Colombia, Tunja 150003, Colombia
5
Instituto de Alta Investigación, CEDENNA, Universidad de Tarapacá, Casilla 7D, Arica 1000000, Chile
*
Author to whom correspondence should be addressed.
Materials 2022, 15(2), 597; https://doi.org/10.3390/ma15020597
Submission received: 23 October 2021 / Revised: 4 January 2022 / Accepted: 7 January 2022 / Published: 13 January 2022
(This article belongs to the Section Electronic Materials)

Abstract

:
We explore the localization properties of a double-stranded ladder within a tight-binding framework where the site energies of different lattice sites are distributed in the cosine form following the Aubry–André–Harper (AAH) model. An imaginary site energy, which can be positive or negative, referred to as physical gain or loss, is included in each of these lattice sites which makes the system a non-Hermitian (NH) one. Depending on the distribution of imaginary site energies, we obtain balanced and imbalanced NH ladders of different types, and for all these cases, we critically investigate localization phenomena. Each ladder can be decoupled into two effective one-dimensional (1D) chains which exhibit two distinct critical points of transition from metallic to insulating (MI) phase. Because of the existence of two distinct critical points, a mixed-phase (MP) zone emerges which yields the possibility of getting a mobility edge (ME). The conducting behaviors of different energy eigenstates are investigated in terms of inverse participation ratio (IPR). The critical points and thus the MP window can be selectively controlled by tuning the strength of the imaginary site energies which brings a new insight into the localization aspect. A brief discussion on phase transition considering a multi-stranded ladder was also given as a general case, to make the present communication a self-contained one. Our theoretical analysis can be utilized to investigate the localization phenomena in different kinds of simple and complex quasicrystals in the presence of physical gain and/or loss.

1. Introduction

The phenomenon of Anderson localization [1,2] was explored almost seven decades ago and it remains a highly active area of research in the discipline of condensed matter physics [3,4] since the localization behavior helps us understand how the mobility of carriers is affected by disorder. After many years of theoretical prediction, experimental verification was performed using an optical lattice setup [5,6] and is nowadays being observed in many cases.
It is well known that for a 1D system with random (uncorrelated) site energies (diagonal disordered system), all energy eigenstates are exponentially localized regardless of the strength of the disorder [1,2]. This indicates that the critical disorder strength W c = 0 (W measures the disorder strength), and thus, such a system is quite trivial as we do not have any option to tune the localization behavior. The situation becomes very interesting when a correlation is imposed among the site energies. There exist many examples of correlated disordered systems that are widely used to study localization phenomena, and among them, the most notable system is probably the Aubry–André–Harper (AAH) [7,8,9,10,11,12,13,14,15,16,17] where site energies are correlated in the cosine form. Unlike in the 1D random disordered system, a 1D AAH chain provides a finite critical point of phase transition, described by the relation W c = 2 t (t being the nearest-neighbor hopping (NNH) strength) [8,18,19,20]. For W < W c , all energy eigenstates are perfectly conducting, whereas beyond this critical limit these states are localized. Thus, a metal-to-insulator (MI) transition occurs at W = W c . Though a finite critical point is obtained for a 1D AAH chain, one never encounters any mixture of both conducting and localized states in the energy band spectrum, i.e., a mixed-phase (MP) [18,20] does not appear. The existence of an MP zone is of course very crucial in many contexts, and it is the primary requirement to have a mobility edge (ME) which separates the conducting zone from the insulating one in the energy eigenspectrum [18,20,21,22,23,24].
To have the mixed-phase energy window and mobility edge, two basic requirements are (i) coupling at least two such 1D AAH chains to form a ladder-like geometry and (ii) diagonal coupling between the strands (see Figure 1). In this context, it is relevant to note that for a strictly 1D AAH chain, one can obtain the mixture of both conducting and insulating states provided the hopping of electrons is not limited to nearest-neighbor sites [25]. In that case, the critical disorder strength is no longer a constant; rather, it depends on the discrete energy eigenvalues. In our present work, we do not consider that situation, and it will maybe be taken into account in any of our forthcoming works, as our present focus is something different. For a two-stranded AAH ladder, a finite energy window is obtained for the mixed phase resulting in a mobility edge [18], and the primary aim of our work is to establish a prescription for tuning the MP window together with the ME.
To substantiate this fact, we included the effect of environmental interaction [26,27,28,29,30,31,32,33,34] in the system, which is actually unavoidable in many realistic situations. Phenomenologically, such an interaction is introduced by adding an imaginary site energy in a lattice site, which makes the system a non-Hermitian (NH) one [35,36,37,38,39,40]. Depending on the sign—which is positive or negative—of the imaginary site energy, we have physical gain or loss in the system. This physical gain/loss might be associated with energy or electrons or even bosonic particles as well, and has already been established in different contemporary works [41,42,43,44,45]. In our analysis, we add imaginary site energies at all distinct lattice sites, along with the cosine modulation following the AAH form, and studied the interplay between the environmental interaction (viz, the NH quantity) [46,47,48,49,50,51] and the correlated disorder on the localization phenomena. Both balanced and imbalanced NH distributions are taken into account [52,53] for a comprehensive analysis.
We simulate the two-stranded NH AAH ladder within a tight-binding (TB) framework, where both the strands are assumed to be identical in nature. Such a ladder can be decoupled into two effective 1D chains, and the localization phenomena of the parent ladder is investigated by critically analyzing the localization behavior of these individual 1D chains [18,20]. The interesting thing is that the two decoupled chains provide two distinct critical points of MI transition, and thus an MP zone appears. The existence of the MP window suggests the availability of an ME in the energy band spectrum. The conducting nature of each energy eigenstate is described by determining the inverse participation ration (IPR) which gives a very good measure as it is directly involved with the participation of particles at different lattice sites. For a conducting state, IPR becomes vanishingly small, whereas it goes to unity for the absolute localized state. Using IPRs of all the distinct energy states, we eventually computed average inverse participation ratio (AIPR) [20] which measures the conducting behavior of the system. AIPR 0 denotes the metallic phase [54,55], while for the insulating phase [56,57] AIPR 1 .
The key features that we want to address in this work are as follows. (i) The specific role of environmental interaction [58,59,60,61,62,63,64,65,66] (NH factor), which is introduced through imaginary site energies, on critical point of phase transition. It essentially gives us a possible tuning mechanism of the critical disorder strength. (ii) The modification of mixed-phase [67,68] energy region and ME, by means of the NH quantity. (iii) Decoupling of the parent ladder system into two effective 1D NH AAH chains. A general case considering a multi-stranded system is briefly discussed at the end, as an extension, for the sake of completeness. These features have not been well explored in the literature, to the best of our knowledge. We believe that the present analysis may provide some insights into the localization phenomena in different kinds of correlated disordered systems subjected to physical gain/loss, and certainly bring a new route for designing more controlled devices.
Following the above brief introduction (Section 1), the remaining parts of the paper are arranged as follows. In Section 2, we describe the NH AAH ladder and the distribution of physical gain and loss. Subsequently, we present the TB Hamiltonian, decoupling the prescription of the ladder and the methodology for studying the localization behavior. In Section 3, we present all the essential results and thoroughly discuss them one by one. The specific role of the NH factor is critically discussed. Finally, in Section 4, we summarize our key findings.

2. Non-Hermitian AAH Ladder, Decoupling Procedure and the Theoretical Framework

2.1. Two-Stranded NH AAH Ladder and the Tight-Binding Hamiltonian

Let us start with Figure 1, where the schematic diagrams of two-stranded ladders are shown. In each case, two identical strands are vertically stacked to form a ladder, where the site energies are modulated by the cosine form obeying the well-known AAH model. The two colored balls are used to refer two kinds of imaginary site energies that are added into the parent AAH lattices, i.e., the base site energies are described in the AAH form where the imaginary part is added, to make the system a non-Hermitian one. In one type of site, a positive imaginary quantity is added (orange site), described as ‘ i δ ’ where δ measures the strength of the NH factor. This is referred as ‘physical gain’ in the system, which can be the gain of energy or we can also call it as the gain of particles. In the other type of site, represented by green ball, a negative energy is included in the form ‘ i δ ’. This is known as ‘physical loss’ in the system. Both the gain and loss in the system are involved with the interaction of the physical system with the environment.
We can arrange the green and orange balls, or more precisely, we can say that the gain and loss terms can be added into the system, in many ways, but here, we consider only four distinct configurations and they are sequentially represented in Figure 1. When the green and orange balls are alternately arranged, we obtain the most stable configuration, and it is called the balanced type-I configuration (see Figure 1a). The gain/loss configurations used in Figure 1b,c are referred as imbalanced type-I and imbalanced type-II, respectively. In one case, all sites are subjected to a constant gain term, while for the other ladder a loss factor is added into the lattice sites. Finally, the ladder shown in Figure 1d is called balanced type-II, where in one half, the gain term is added and the other half loss factor is included. For all these four types of NH AAH ladders, we critically investigate the localization properties and specific role of the NH factor.
Any such ladder can be simulated within a tight-binding framework, and the general form of the TB Hamiltonian looks like:
H = n , m ϵ n , m e f f c n , m c n , m + n , m t l ( c n + 1 , m c n , m + h . c . ) + n , m t p ( c n , m + 1 c n , m + h . c . ) + n , m t d [ ( c n + 1 , m + 1 c n , m + c n + 1 , m c n , m + 1 ) + h . c . ] .
The integer n is used to describe the lattice sites in each strand which runs from 1 to N (N being the total number of lattice sites in each single strand), whereas the strand index is described by the integer m and it runs from 1 to M (M denotes the total number of strands). c n , m , c n , m are the conventional Fermionic operators at the lattice site ( n , m ) (nth site of mth strand). The parameters t l , t p and t d describe the intra-strand, inter-strand and diagonal hopping integrals, respectively, as clearly depicted in Figure 1. ϵ n , m e f f is the effective site energy which is the sum of two terms: one is associated with the cosine modulation, i.e., the AAH term and the other is the imaginary site energy in the form of + i δ or i δ . Since the AAH site energies are taken to be identical for both strands, we can express ϵ n , m e f f as
ϵ n , m e f f = ϵ n e f f = ϵ n + δ n
where:
ϵ n = W cos ( 2 π b n ) .
The index m in effective site energy is removed without any loss of generality, as for all the strands, ϵ n s have an identical form. W measures the AAH modulation strength, which is usually referred to as the correlated disorder strength, and b is an incommensurate quantity. In our analysis, we chose b = ( 1 + 5 ) / 2 (golden mean) [13,69,70]. This typical value is widely used in the literature, though any other irrational number can also be considered, and the physics will remain unchanged. δ n represents the gain or loss term, viz, δ n = i δ or δ n = i δ depending on the arrangement of the orange and/or green colored balls in the AAH ladders (see Figure 1).

2.2. Decoupling of Two-Stranded AAH Ladder into Two Effective 1D Chains

The decoupling of the ladder into individual 1D chains is an important perspective for analyzing the mixed-phase window, mobility edge phenomenon, and the overall localization behavior of a full system. The decoupling can be made in several ways, as already noted by different groups including us [18,19,20]. Here, we use a simple prescription that is described as follows.
Starting from the Schrödinger equation H | Ψ = E | Ψ , we can write the tight-binding difference equation as
( E ϵ n e f f ) Ψ n , m = t l ( Ψ n + 1 , m + Ψ n 1 , m ) + t p ( Ψ n , m + 1 + Ψ n , m 1 ) + t d ( Ψ n + 1 , m + 1 + Ψ n 1 , m + 1 + Ψ n + 1 , m 1 + Ψ n 1 , m 1 )
where Ψ n , m ’s are the wave amplitudes. As the site energies are identical along each rung (referred as the y direction), we can write a plane wave solution along the y direction, and hence:
Ψ n , m = B ψ n q e i m k y a
where B is a constant factor, ψ n q is the wave amplitude for the x direction, k y is the wave vector, q is the strand index and a is the lattice spacing. Substituting such wave amplitude in the above equation, we have:
( E ϵ n e f f ) B ψ n q e i m k y a = t l ( B ψ n + 1 q e i m k y a + B ψ n 1 q e i m k y a ) + t p ( B ψ n q e i ( m + 1 ) k y a + B ψ n q e i ( m 1 ) k y a ) + t d ( B ψ n + 1 q e i ( m + 1 ) k y a + B ψ n 1 q e i ( m + 1 ) k y a + B ψ n + 1 q e i ( m 1 ) k y a + B ψ n 1 q e i ( m 1 ) k y a ) .
Now, eliminating the term B e i m k y a from both sides of the above equation, we obtain:
( E ϵ n e f f ) ψ n q = t l ( ψ n + 1 q + ψ n 1 q ) + t p ( ψ n q e i k y a + ψ n q e i k y a ) + t d ( ψ n + 1 q e i k y a + ψ n 1 q e i k y a + ψ n + 1 q e i k y a + ψ n 1 q e i k y a ) .
Doing some simple algebraic steps, we reach the expression:
E ( ϵ n e f f + 2 t p cos ( k y a ) ψ n q = t l + 2 t d cos ( k y a ) ( ψ n + 1 q + ψ n 1 q ) .
The wave vector k y is governed by the relation k y a = q π / ( M + 1 ) . Now, considering M = 2 ( q = 1 , 2), we eventually reach the following equations:
[ E { W cos ( 2 π b n ) + δ n + t p } ] ψ n 1 =   t l + t d [ ψ n + 1 1 + ψ n 1 1 ] ,
[ E { W cos ( 2 π b n ) + δ n t p } ] ψ n 2 =   t l t d [ ψ n + 1 2 + ψ n 1 2 ] .
These equations can be further expressed in compact forms as
[ E ϵ n 1 , e f f ] ψ n 1 = t 1 , e f f [ ψ n + 1 1 + ψ n 1 1 ]
[ E ϵ n 2 , e f f ] ψ n 2 = t 2 , e f f [ ψ n + 1 2 + ψ n 1 2 ]
The effective site energies are:
ϵ n 1 , e f f = W cos ( 2 π b n ) + δ n + t p
and:
ϵ n 2 , e f f = W cos ( 2 π b n ) + δ n t p ,
and the effective hopping integrals are: t 1 , e f f = ( t l + t d ) and t 2 , e f f = ( t l t d ) . Equations (10a) and (10b) represent two effective 1D chains which are obtained by decoupling the parent NH AAH ladder. Now, depending on the distribution of the gain/loss term ( δ n ), we obtain different effective site energies for the decoupled chains.

2.3. Theoretical Prescription

We compute the localization behavior of each energy eigenstate by studying the inverse participation ratio. For any normalized eigenstate (say) | ϕ q = n = 1 N | Ψ n q , the IPR is defined as
I P R q = n = 1 N | Ψ n q | 4 .
This quantity clearly describes whether a state is conducting or non-conducting in behavior. For the conducting state, we obtain IPR q 0 , while it becomes to unity for the non-conducting sate. It is usually quite difficult to have IPR q 1 when N is finite, since only in the asymptotic limit ( N ) can we attain the ‘absolute localized phase’.
Determining all the IPR q s for the individual energy eigenstates, we finally compute the average inverse participation ratio that describes the conducting properties of the system. The AIPR [20,54,55,56,57] is defined as
A I P R = I P R = 1 N q = 1 N I P R q .
For the conducting phase, AIPR 0 , while AIPR 1 denotes the localized phase of the system. For the intermediate values of AIPR, we obtain the mixed phase.

3. Numerical Results and Discussion

In what follows, we present and discuss our results which include the energy band spectra, the control of the critical point of transition from metallic to insulating phase, the regulation of mixed-phase window and related issues. From now on, we abbreviate ‘balanced type’ as BT and ‘imbalanced type’ as IT, for simplification, and here we thoroughly investigate the results of four different types of NH AAH ladders (BT-I, IB-I, IB-II and BT-II), those are schematically shown in Figure 1.
Before presenting the results, let us mention the parameter values which are kept constant throughout the computation. We set t l = 1.5 eV, t p = 1 eV and t d = 0.5 eV. The number of lattice sites in each strand of the ladder is fixed at N = 250 , unless specified. The other parameters that are not so common are mentioned in the appropriate places. All the energies are measured in unit of electron-volt (eV).

3.1. Two-Stranded AAH Ladders in Presence of Physical Gain and/or Loss

The central focus of our work was to investigate the critical role of physical gain and/or loss on electronic localization. Before coming to that part, it is indeed necessary to analyze the energy level diagram, as this is the most fundamental thing to understand a physical system. From our mathematical analysis, we find that a two-stranded ladder is decoupled into two effective 1D NH AAH chains, described by the relations given in Equations (10a) and (10b). Both these two decoupled chains obtain effective site energies and NNH integrals. The chain possessing the higher effective hopping integral (viz, t l + t d ) is called chain-1, while the other chain with the lower hopping strength (i.e., t l t d ) is called chain-2.

3.1.1. Eigenvalue Spectrum

In the presence of the NH factor, the eigenvalues become complex, and in order to have a clearer picture, in Figure 2, we show the variations of both the real and imaginary parts of the eigenvalues as a function of δ , and for all the four different types of ladders. The orange and blue colors are involved with the chain-1 and chain-2, respectively. Several interesting features emerge from the band spectra. For the BT-I case, a significant change in real eigenvalues with δ was obtained when δ is quite small (Figure 2a), while the effect becomes relatively weaker for higher δ . This is primarily due to the interplay between AAH potential and the NH quantity. For low δ , both these factors contribute and the effect of AAH modulation becomes gradually diminished with increasing δ . Because of this fact, we initially obtain three sub-bands in the real eigenvalue spectrum which is the generic feature of AAH modulation, while this behavior almost vanishes when δ is reasonably large. Instead of three sub-bands, we almost obtain a continuous band. As alternate sites contain gain and loss terms for this BT-I case, the real eigenvalues are practically constant beyond some critical limit of δ . The finite separation between the orange and blue bands are associated with the effective tight-binding parameters of the decoupled chains.
For the other two ladders denoted by IT-I and IT-II, the real eigenvalues are insensitive to δ (as can be seen in Figure 2c,e), which is easy to understand since all the sites contain a constant gain or loss term. Finally, for the ladder where the gain term is included in one half and the loss term is added in the other half (BT-II ladder), a slight variation in the real spectrum is obtained, and that is also for the smaller δ (Figure 2g).
The variations of imaginary eigenvalues with δ are also quite interesting. For the IT-I and IT-II ladders, the eigenvalues of the decoupled chains (chain-1 and chain-2) completely merge with each other (Figure 2d,f), and the increment or reduction in eigenvalues with δ is very simple to understand since all sites contain + i δ and i δ , respectively, for the IT-I and IT-II ladders, respectively. Thus, the variation follows a straight line with δ . Similarly to the real eigenvalue spectrum, the imaginary eigenvalues are also highly sensitive to the NH factor for the BT-I ladder, and this effect is more prominent for a smaller δ (Figure 2b). Furthermore, finally, for the BT-II ladder, the imaginary eigenvalues follow almost linear dependence both along the positive and negative directions (Figure 2h). This is solely associated with the two halves of the ladder that contain positive and negative imaginary site energies.
Thus, what we find from the real and imaginary eigenvalue spectra shown in Figure 2 is that for the AAH ladder with balanced gain and loss distribution, the eigenvalues are most sensitive and a more complex pattern is obtained when δ is relatively small, i.e., comparable to the AAH modulation strength. This is due to the combined effects of these two factors (viz, W and δ ).

3.1.2. Tuning of Critical Points, Mixed-Phase Zone and Mobility Edge

Now, we focus on the central feature of our analysis, i.e., the tuning of critical points and the localization phenomena by means of the NH parameter δ .
In Figure 3, we show the dependence of AIPR as a function of the correlated disorder strength W for all four AAH ladders, at three distinct values of δ . In each spectrum, the AIPR- δ curves associated with chain-1 (black line) and chain-2 (red line) are superimposed. It is observed that the two chains exhibit two separate critical points of transition from the metallic phase to the insulating one. Here, we are not expecting a sharp jump of AIPR from its almost zero value to unity, since the system size is finite, as mentioned previously. However, the signature of the phase transition is clearly reflected. Because of the two different critical points, we obtain a window of W for which all the eigenstates of one chain (chain-1) are completely extended in behavior, while for the other chain (chain-2), all the states are localized. As a result of this, a mixed-phase window appears (lightly shaded region), and at the same time, we obtain a mobility edge that creates a partition among the conducting and insulating states. When δ = 0 , the critical disorder strengths of the transitions of the two decoupled chains are: W c 1 = 2 ( t l + t d ) and W c 2 = 2 ( t l t d ) (two times the effective NNH strength for the diagonal AAH chain). This is well known in the literature, and has been discussed by many groups including ourselves [18,19,20]. The main goal of this communication was to verify whether the transition points can be tuned by means of environmental interaction δ , in the presence of the cosine modulation. If yes, then it will definitely be an interesting result and might be utilized in several ways, especially in the context of localization phenomena.
Carefully looking into the spectra given in the second and third rows of Figure 3, we can see that, for the IT-I and IT-II ladders, there is absolutely no change in the transition points, and thus the MP window, due to δ . A very slight change can be noticed for the BT-II ladder (last row), and an appreciable modification is only obtained for the BT-I AAH ladder (top row). These features can be explained as follows. In the presence of the physical gain and/or loss, a phase factor is introduced at each and every lattice site. Depending on the gain or loss term, the probability of finding a particle at different lattice sites becomes increased or decreased. When all sites possess either a gain or loss factor, for instance, in the IT-I and IT-II ladders, the phases are the same at all the sites which essentially do not produce any new scattering and hence there is no effect on the localization. In the BT-II ladder, a phase mismatch only occurs at one particular place (the center of the ladder), where the sign reversal of δ takes place. Since at all other regions the phases are identical, the scattering of particles only at the center does not produce any appreciable effect, and it becomes much weaker with the increasing ladder length. The noticeable effect of δ is only observed in BT-I ladder as in this case, a regular phase change takes place in every alternate lattice sites, and thus the scattering becomes stronger which yields the localization. Naturally, the critical points of phase transition for both the two decoupled chains, i.e., W c 1 and W c 2 become reduced with the increase in the interaction parameter δ , as it enhances more scattering at the lattice sites.
The results presented in Figure 3 clearly indicate that the BT-I ladder is most sensitive on the interaction parameter δ . Therefore, in the rest of our analysis, we only concentrate on the BT-I ladder.
In order to see the direct effect of δ on different eigenstates and for a better viewing of the mobility edge, in Figure 4, we show the density plot of IPRs of all the eigenstates as functions of the eigenvalues and the AAH modulation strength. The two columns are associated with the two decoupled chains of the BT-1 ladder. The results for δ = 0 are also included (where all the eigenvalues are real in nature), to obtain a comparative analysis on IPRs when δ is finite. For δ = 0 , we see that in a wide range of W, all the eigenstates of chain-1 exhibit conducting behavior (see Figure 4b), while for the other chain (chain-2), all the states become localized at a low W (Figure 4a). This reveals the co-existence of the conducting and localized eigenstate, and thus the mobility edge. The appearance of the mixed phase is fully consistent with the results discussed in Figure 3 for the BT-I ladder. For non-zero δ , the density plots are quite anomalous, especially for the imaginary case eigenvalues, though the band splitting and gapped-like spectrum still persists. These are the generic features of AAH systems. What we clearly find is that states become localized at a much lesser W when we include the effect of δ which clearly corroborates our findings in Figure 3 for the BT-I ladder. The physical explanation behind the reduction in critical points of both the two decoupled chains with δ was already given earlier.
Based on the characteristic features of Figure 3 and Figure 4, it is indeed necessary to check how the critical points are affected by a continuous variation of δ . The results are presented in Figure 5, where the variations of W c 1 and W c 2 are shown as a function of δ for the BT-I ladder. We consider a large range of δ , and compute the critical points of transitions for the two decoupled chains. For both these chains, the critical disorder strengths varies continuously, providing a smooth variation with respect to δ . The green circles and the blue squares represent the numerically obtained data points, and these data are fitted well with closed mathematical forms in terms of δ . Here, it is relevant to note that there is no physical reason for getting this dependence of the critical point with δ . It is always easier to understand any property if we have a closed form or mathematical expression in terms of the variable parameter(s) associated with the system. We can therefore even estimate the value of the critical point even for any arbitrary values of δ . For our case, we find that the data points W c 1 and W c 2 are fitted well with these functional forms. One can of course fit the data with other functional form(s), but we need to focus on the best fit, and here we confirm these relations through numerous attempts. The results of Figure 5 reveal that the difference between the two curves is initially large, but the difference gradually decreases with the increasing δ . The gap practically disappears for a sufficiently large δ (not shown here in the figure), which we confirm through our detailed numerics.
This is quite easily understandable, as for a much higher δ , both the two chains become localized for a vanishingly small W. Due to the reduction in critical disorder strengths with δ , the width of the mixed-phase window becomes narrowed providing a shifting of the mobility edge. Eventually, at the stage when the width becomes small to the point of vanishing, it is difficult to achieve a mobility edge since the coexistence of both the localized and extended states is the primary requirement to have the mobility edge.
To make the present analysis a self-contained one, here we focus on the scaling behavior of AIPR on the system size N. The results are presented in Figure 6 for the BT-I ladder under different input conditions of δ and W. Very nicely, we find that for each and every set of input parameters, AIPR follows a nice pattern with N, and therefore, we can fit the data using a simple closed form (given in the inset of each sub-figure).

3.2. M-Stranded Ladder: A General Case

To make the present analysis a self-contained one, we finally focus on an M-stranded ladder. The value of M can be anything. Any such ladder is constructed by vertically stacking M number of identical AAH NH chains, and the neighboring chains are coupled in the same way as was discussed above for the two-stranded ladder (viz, vertically as well as diagonally). As the site energies in each rung, possessing M number of lattice sites, are identical, we can write a plane wave solution along the y direction, and following the same mathematical prescription given in Section 2.2, we can decouple M-stranded ladder into M-distinct effective 1D chains. The different equation of M such chains involving the wave amplitudes are written in a compact form as
E W cos ( 2 π b n ) + δ n + 2 t p cos q π M + 1 ψ n q = t l + 2 t d cos q π M + 1 ψ n + 1 q + ψ n 1 q
where the integer q runs from 1 to M. Once we obtain the decoupled chains, we can easily determine the critical disorder strengths W c m for the M different effective chains by computing AIPR as a function of the AAH modulation strength W.
In Figure 7, we plot the critical points of transition for different stranded AAH ladders at four typical values of the NH parameter δ . The results are computed by varying M from 1 to 10, where we only choose the BT-I case since this type of ladder is most sensitive to δ . The values of the critical disorder strengths are represented by the colored dot points in each spectrum. At a first glance, we see that the dot points follow an almost regular pattern, but a careful inspection reveals that the uniformity starts to become distorted with the increasing strength of δ . Most importantly, we find that the values of W c m gradually decrease with δ . The reduction in critical points is solely associated with more scattering of particles in the regime of higher δ . The key feature that emerges from the spectra is that, with the increasing number of strands, more and more distinct critical points are generated which yield multiple mixed-phase zones in the energy spectrum. Because of this, we have a finite possibility of obtaining multiple mobility edges. Moreover, we can tune the mixed-phase windows as well as the mobility edges by means of the interaction parameter δ . These aspects are highly significant to achieve high to low conducting switching action and vice versa, and to selectively tune the localization properties with the help of δ .

4. Closing Remarks

In this work, we investigated the localization properties of Aubry–André–Harper ladder networks in the presence of environmental interaction. The interaction of the system with the environment was phenomenologically incorporated by adding imaginary site energies, and depending on the sign, the interaction is referred as physical gain (positive imaginary) or physical loss (negative imaginary). We have considered four distinct configurations associated with gain/loss terms which are BT-I, IT-I, IT-II and BT-II, and for all these cases, we critically analyzed the localization phenomena. We elaborately discussed different cases of two-stranded Aubry–André–Harper ladders, and at the end, we briefly illustrated the behavior of a multi-stranded ladder as a general case. All these ladders were simulated within a tight-binding framework, which is most convenient for describing a physical system, and the localization aspects were analyzed by calculating the average inverse participation ratio. The key features and new aspects of our work are as follows:
  • For the first time, to the best of our knowledge, localization phenomena were investigated considering a multi-stranded ladder network in the presence of environmental interaction;
  • A simple prescription for decoupling any arbitrary stranded ladder system was provided. This helped us detect the mixed-phase window(s) and the mobility edge(s);
  • From a comparative analysis, we established that when the gain and loss terms are added in alternate sites, maximum scattering occurs which leads to a significant impact on the localization behavior. For the ladders with either gain or loss terms, there is no effect on localization. A minor effect was noticed for the balanced type-II ladder which is very easy to understand;
  • The mixed-phase window(s), and thus, the mobility edge(s) can be monitored by means of the interaction factor δ . This is an interesting observation which might be helpful in deriving controlled transport properties;
  • Our analysis is not specific to any particular system—rather, we can easily extend it to any other quasicrystals in the presence of a non-Hermitian factor that exhibits a fragmented and gapped spectrum.

Author Contributions

This paper was written based as a collaborative work by the authors. S.R. and S.K.M. conceived the project. S.R., L.M.P. and J.H.O.S. performed the numerical calculations. S.R., D.L. and S.K.M. analyzed the data. All authors have read and agreed to the published version of the manuscript.

Funding

S.R. would like to thank CSIR, India (09/093(0178)/2017-EMR-I), for providing their research fellowship. The research of S.K.M. is supported by DST-SERB, Government of India (Project File Number: EMR/2017/000504). L.M.P. acknowledges financial support from ANID through Convocatoria Nacional Subvención a Instalación en la Academia Convocatoria A n ¯ o 2021, Grant SA77210040. D.L. acknowledges partial financial support from FONDECYT 1180905. D.L. also acknowledges partial financial support from Centers of Excellence with BASAL/ANID financing Grant AFB180001, CEDENNA. J.H.O.S. acknowledges to Universidad Pedagógica y Tecnológica de Colombia.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Anderson, P.W. Absence of Diffusion in Certain Random Lattices. Phys. Rev. 1958, 109, 1492. [Google Scholar] [CrossRef]
  2. Lee, P.A.; Ramakrishnan, T.V. Disordered electronic systems. Rev. Mod. Phys. 1985, 57, 287. [Google Scholar] [CrossRef]
  3. Nandkishore, R.; Huse, D.A. Many-Body Localization and Thermalization in Quantum Statistical Mechanics. Annu. Rev. Comdens. Matter Phys. 2015, 6, 15. [Google Scholar] [CrossRef] [Green Version]
  4. Vosk, R.; Huse, D.A.; Altman, E. Theory of the Many-Body Localization Transition in One-Dimensional Systems. Phys. Rev. X 2015, 5, 031032. [Google Scholar] [CrossRef] [Green Version]
  5. Billy, J.; Josse, V.; Zuo, Z.; Bernard, A.; Hambrecht, B.; Lugan, P.; Clément, D.; Sanchez-Palencia, L.; Bouyer, P.; Aspect, A. Direct observation of Anderson localization of matter waves in a controlled disorder. Nature 2008, 453, 891. [Google Scholar] [CrossRef]
  6. Roati, G.; Errico, C.D.; Fallani, L.; Fattori, M.; Fort, C.; Zaccanti, M.; Modugno, G.; Modugno, M.; Inguscio, M. Anderson localization of a non-interacting Bose–Einstein condensate. Nature 2008, 453, 895. [Google Scholar] [CrossRef] [Green Version]
  7. Harper, P.G. The General Motion of Conduction Electrons in a Uniform Magnetic Field, with Application to the Diamagnetism of Metals. Proc. Phys. Soc. Lond. Sect. A 1955, 68, 874. [Google Scholar] [CrossRef]
  8. Aubry, S.; Andŕe, G. Analyticity breaking and Anderson localization in incommensurate lattices. Ann. Isr. Phys. Soc. 1980, 3, 133. [Google Scholar]
  9. Ganeshan, S.; Sun, K.; Sarma, S.D. Topological Zero-Energy Modes in Gapless Commensurate Aubry-André-Harper Models. Phys. Rev. Lett. 2013, 110, 180403. [Google Scholar] [CrossRef] [Green Version]
  10. An, F.A.; Padavić, K.; Meier, E.J.; Hegde, S.; Ganeshan, S.; Pixley, J.H.; Vishveshwara, S.; Gadway, B. Interactions and Mobility Edges: Observing the Generalized Aubry-André Model. Phys. Rev. Lett. 2021, 126, 040603. [Google Scholar] [CrossRef]
  11. Dey, S.; Daw, D.; Maiti, S.K. Flux-driven circular current and near-zero field magnetic response in an Aubry ring: High-to-low conducting switching action. Europhys. Lett. 2020, 129, 47002. [Google Scholar] [CrossRef]
  12. Patra, M.; Maiti, S.K.; Sil, S. Engineering magnetoresistance: A new perspective. J. Phys. Condens. Matter 2019, 31, 355303. [Google Scholar] [CrossRef] [Green Version]
  13. Maiti, S.K.; Sil, S.; Chakrabarti, A. Phase controlled metal–insulator transition in multi-leg quasiperiodic optical lattices. Ann. Phys. 2017, 382, 150. [Google Scholar] [CrossRef] [Green Version]
  14. Saha, S.; Maiti, S.K.; Karmakar, S.N. Multiple mobility edges in a 1D Aubry chain with Hubbard interaction in presence of electric field: Controlled electron transport. Physica E 2016, 83, 358. [Google Scholar] [CrossRef] [Green Version]
  15. Dey, M.; Chakraborty, S.; Maiti, S.K. New Route to Enhanced Figure of Merit at Nano Scale: Effect of AAH Modulation. J. Phys. D Appl. Phys. 2022, 55, 085302. [Google Scholar] [CrossRef]
  16. Koley, A.; Maiti, S.K.; Pérez, L.M.; Silva, J.H.O.; Laroze, D. Possible Routes to Obtain Enhanced Magnetoresistance in a Driven Quantum Heterostructure with a Quasi-Periodic Spacer. Micromachines 2021, 12, 1021. [Google Scholar] [CrossRef]
  17. Koley, A.; Maiti, S.K.; Ojeda Silva, J.H.; Laroze, D. Spin Dependent Transport through Driven Magnetic System with Aubry-André-Harper Modulation. Appl. Sci. 2021, 11, 2309. [Google Scholar] [CrossRef]
  18. Sil, S.; Maiti, S.K.; Chakrabarti, A. Metal-Insulator Transition in an Aperiodic Ladder Network: An Exact Result. Phys. Rev. Lett. 2008, 101, 076803. [Google Scholar] [CrossRef] [Green Version]
  19. Roy, S.; Maiti, S.K. Tight-binding quantum network with cosine modulations: Electronic localization and delocalization. Eur. Phys. J. B 2019, 92, 267. [Google Scholar] [CrossRef]
  20. Rossigonolo, M.; Dell’Anna, L. Localization transitions and mobility edges in coupled Aubry-André chains. Phys. Rev. B 2019, 99, 054211. [Google Scholar] [CrossRef] [Green Version]
  21. Maiti, S.K.; Nitzan, A. Mobility edge phenomenon in a Hubbard chain: A mean field study. Phys. Lett. A 2013, 377, 1205. [Google Scholar] [CrossRef] [Green Version]
  22. Sarkar, S.; Maiti, S.K. Localization to delocalization transition in a double stranded helical geometry: Effects of conformation, transverse electric field and dynamics. J. Phys. Condens. Matter 2020, 32, 505301. [Google Scholar] [CrossRef] [PubMed]
  23. Mal, B.; Banerjee, M.; Maiti, S.K. Magnetotransport in fractal network with loop sub-structures: Anisotropic effect and delocalization. Phys. Lett. A 2020, 384, 126378. [Google Scholar] [CrossRef]
  24. Banerjee, M.; Mal, B.; Maiti, S.K. Unconventional localization phenomena in a spatially non-uniform disordered material. Physica E 2019, 106, 312. [Google Scholar] [CrossRef]
  25. Biddle, J.; Sarma, S.D. Predicted mobility edges in one-dimensional incommensurate optical lattices: An exactly solvable model of Anderson localization. Phys. Rev. Lett. 2010, 104, 070601. [Google Scholar] [CrossRef] [Green Version]
  26. Dities, F.M. The decay of quantum systems with a small number of open channels. Phys. Rep. 2000, 339, 215. [Google Scholar] [CrossRef]
  27. Amir, A.; Oreg, Y.; Imry, Y. Decays in quantum hierarchical models. Phys. Rev. A 2008, 77, 050101. [Google Scholar] [CrossRef] [Green Version]
  28. Sokolov, V.V.; Zelevinsky, V.G. Collective dynamics of unstable quantum states. Ann. Phys. 1992, 216, 323. [Google Scholar] [CrossRef]
  29. Rotter, I. A continuum shell model for the open quantum mechanical nuclear system. Rep. Prog. Phys. 1991, 54, 635. [Google Scholar] [CrossRef]
  30. Sokolov, V.V.; Zelevinsky, V.G. Dynamics and statistics of unstable quantum states. Nucl. Phys. A 1989, 504, 562. [Google Scholar] [CrossRef]
  31. Sokolov, V.V.; Zelevinsky, V.G. On a statistical theory of overlapping resonances. Phys. Lett. B 1988, 202, 10. [Google Scholar] [CrossRef]
  32. Hatano, N.; Nelson, D.R. Localization Transitions in Non-Hermitian Quantum Mechanics. Phys. Rev. Lett. 1996, 77, 570. [Google Scholar] [CrossRef] [Green Version]
  33. Hatano, N.; Nelson, D.R. Vortex pinning and non-Hermitian quantum mechanics. Phys. Rev. B 1997, 56, 8651. [Google Scholar] [CrossRef] [Green Version]
  34. Carmele, A.; Heyl, M.; Kraus, C.; Dalmonte, M. Stretched exponential decay of Majorana edge modes in many-body localized Kitaev chains under dissipation. Phys. Rev. B 2015, 92, 195107. [Google Scholar] [CrossRef] [Green Version]
  35. Bender, C.M. Making sense of non-Hermitian Hamiltonians. Rep. Prog. Phys. 2007, 70, 947. [Google Scholar] [CrossRef] [Green Version]
  36. Bender, C.M.; Boettcher, S. Real Spectra in Non-Hermitian Hamiltonians Having PT Symmetry. Phys. Rev. Lett. 1998, 80, 5243. [Google Scholar] [CrossRef] [Green Version]
  37. Bender, C.M.; Boettcher, S.; Meisinger, P.N. PT-symmetric quantum mechanics. J. Math. Phys. 1999, 40, 2201. [Google Scholar] [CrossRef] [Green Version]
  38. Marinello, G.; Pato, M.P. Random non-Hermitian tight-binding models. J. Phys. Conf. Ser. 2016, 738, 012040. [Google Scholar] [CrossRef]
  39. Burke, P.C.; Wiersig, J.; Haque, M. Non-Hermitian scattering on a tight-binding lattice. Phys. Rev. A 2020, 102, 012212. [Google Scholar] [CrossRef]
  40. Liu, T.; Cheng, S.; Guo, H.; Gao, X. Fate of Majorana zero modes, exact location of critical states, and unconventional real-complex transition in non-Hermitian quasiperiodic lattices. Phys. Rev. B 2021, 103, 104203. [Google Scholar] [CrossRef]
  41. Ramezani, H.; Kottos, T.; El-Ganainy, R.; Christodoulides, D.N. Unidirectional nonlinear PT-symmetric optical structures. Phys. Rev. A 2010, 82, 043803. [Google Scholar] [CrossRef] [Green Version]
  42. Rüter, C.E.; Makris, K.G.; El-Ganainy, R.; Christodoulides, D.N.; Segev, M.; Kip, D. Observation of parity–time symmetry in optics. Nat. Phys. 2010, 6, 192. [Google Scholar] [CrossRef] [Green Version]
  43. Bendix, O.; Fleischmann, R.; Kottos, T.; Shapiro, B. Exponentially Fragile PT-Symmetry in Lattices with Localized Eigenmodes. Phys. Rev. Lett. 2009, 103, 030402. [Google Scholar] [CrossRef] [Green Version]
  44. Graefe, E.M.; Korsch, H.J.; Niederle, A.E. Mean-Field Dynamics of a Non-Hermitian Bose-Hubbard Dimer. Phys. Rev. Lett. 2008, 101, 150408. [Google Scholar] [CrossRef] [Green Version]
  45. West, C.T.; Kottos, T.; Prosen, T. PT-Symmetric Wave Chaos. Phys. Rev. Lett. 2010, 104, 054102. [Google Scholar] [CrossRef] [Green Version]
  46. Bendix, O.; Fleischmann, R.; Kottos, T.; Shapiro, B. Optical structures with local PT-symmetry. J. Phys. A Math. Theor. 2010, 43, 265305. [Google Scholar] [CrossRef] [Green Version]
  47. Oberthaler, M.K.; Abfalterer, R.; Bernet, S.; Schmiedmayer, J.; Zeilinger, A. Atom Waves in Crystals of Light. Phys. Rev. Lett. 1996, 77, 4980. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Luo, X.B.; Huang, J.H.; Zhong, H.H.; Qin, X.Z.; Xie, Q.T.; Kivshar, Y.S.; Lee, C.H. Pseudo-parity-time symmetry in optical systems. Phys. Rev. Lett. 2013, 110, 243902. [Google Scholar] [CrossRef] [PubMed]
  49. Chong, Y.D.; Ge, L.; Stone, A.D. PT-Symmetry Breaking and Laser-Absorber Modes in Optical Scattering Systems. Phys. Rev. Lett. 2011, 106, 093902. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Sukhorukov, A.A.; Xu, Z.; Kivshar, Y.S. Nonlinear suppression of time reversals in PT-symmetric optical couplers. Phys. Rev. A 2010, 82, 043818. [Google Scholar] [CrossRef] [Green Version]
  51. Klaiman, S.; Günther, U.; Moiseyev, N. Visualization of Branch Points in PT-Symmetric Waveguides. Phys. Rev. Lett. 2008, 101, 080402. [Google Scholar] [CrossRef] [Green Version]
  52. Levi, E.; Heyl, M.; Lesanovsky, I.; Garrahan, J.P. Robustness of Many-Body Localization in the Presence of Dissipation. Phys. Rev. Lett. 2016, 116, 237203. [Google Scholar] [CrossRef]
  53. Zeng, Q.B.; Chen, S.; Lü, R. Anderson localization in the non-Hermitian Aubry-André-Harper model with physical gain and loss. Phys. Rev. A 2017, 95, 062118. [Google Scholar] [CrossRef] [Green Version]
  54. Thouless, D.J. Electrons in disordered systems and the theory of localization. Phys. Rep. 1974, 13, 93. [Google Scholar] [CrossRef]
  55. Wegner, F.J. Inverse participation ratio in 2+ϵ dimensions. Z. Phys. B 1980, 36, 209. [Google Scholar] [CrossRef]
  56. Ganeshan, S.; Pixley, J.H.; Sarma, S.D. Nearest Neighbor Tight Binding Models with an Exact Mobility Edge in One Dimension. Phys. Rev. Lett. 2015, 114, 146601. [Google Scholar] [CrossRef] [Green Version]
  57. Biddle, J.; Wang, B.; Priour, D.J., Jr.; Sarma, S.D. Localization in one-dimensional incommensurate lattices beyond the Aubry-André model. Phys. Rev. A 2009, 80, 021603. [Google Scholar] [CrossRef] [Green Version]
  58. Guo, A.; Salamo, G.J.; Duchesne, D.; Morandotti, R.; Volatier-Ravat, M.; Aimez, V.; Siviloglou, G.A.; Christodoulides, D.N. Observation of PT-Symmetry Breaking in Complex Optical Potentials. Phys. Rev. Lett. 2009, 103, 093902. [Google Scholar] [CrossRef] [Green Version]
  59. Amir, A.; Hatano, N.; Nelson, D.R. Non-Hermitian localization in biological networks. Phys. Rev. E 2016, 93, 042310. [Google Scholar] [CrossRef] [Green Version]
  60. Feng, L.; Ayache, M.; Huang, J.; Xu, Y.-L.; Lu, M.-H.; Chen, Y.-F.; Fainman, Y.; Scherer, A. Nonreciprocal light propagation in a silicon photonic circuit. Science 2011, 333, 729. [Google Scholar] [CrossRef] [Green Version]
  61. Regensburger, A.; Bersch, C.; Miri, M.-A.; Onishchukov, G.; Christodoulides, D.N.; Peschel, U. Parity-time synthetic photonic lattices. Nature 2012, 488, 167. [Google Scholar] [CrossRef]
  62. Bender, C.M.; Berntson, B.K.; Parker, D.; Samuel, E. Observation of PT-phase transition in a simple mechanical system. Am. J. Phys. 2013, 81, 173. [Google Scholar] [CrossRef] [Green Version]
  63. Mejia-Cortes, C.; Molina, M.I. Interplay of disorder and PT symmetry in one-dimensional optical lattices. Phys. Rev. A 2015, 91, 033815. [Google Scholar] [CrossRef]
  64. Jovic, D.M.; Denz, C.; Belic, M.R. Anderson localization of light in PT-symmetric optical lattices. Opt. Lett. 2012, 37, 4455. [Google Scholar] [CrossRef]
  65. Schindler, J.; Li, A.; Zheng, M.C.; Ellis, F.M.; Kottos, T. Experimental study of active LRC circuits with PT symmetries. Phys. Rev. A 2011, 84, 040101. [Google Scholar] [CrossRef] [Green Version]
  66. Hofmann, T.; Helbig, T.; Schindler, F.; Salgo, N.; Brzezinska, M.; Greiter, M.; Kiessling, T.; Wolf, D.; Vollhardt, A.; Kabasi, A.; et al. Reciprocal skin effect and its realization in a topolectrical circuit. Phys. Rev. Res. 2020, 2, 023265. [Google Scholar] [CrossRef]
  67. Zeng, Q.; Yang, Y.; Xu, Y. Topological phases in non-Hermitian Aubry-André-Harper models. Phys. Rev. B 2020, 101, 020201(R). [Google Scholar] [CrossRef] [Green Version]
  68. Yahyavi, M.; Hetényi, B.; Tanatar, B. Generalized Aubry-André-Harper model with modulated hopping and p-wave pairing. Phys. Rev. B 2019, 100, 064202. [Google Scholar] [CrossRef] [Green Version]
  69. Patra, M.; Maiti, S.K. Controlled charge and spin current rectifications in a spin polarized device. J. Magn. Magn. Mater. 2019, 484, 408. [Google Scholar] [CrossRef] [Green Version]
  70. Patra, M.; Maiti, S.K. Externally controlled high degree of spin polarization and spin inversion in a conducting junction: Two new approaches. Sci. Rep. 2017, 7, 14313. [Google Scholar] [CrossRef]
Figure 1. (Color online). Schematic diagrams of non-Hermitian Aubry–André–Harper ladders following the four different arrangements of the imaginary site energies those are respectively presented in (a), (b), (c) and (d). The green ball corresponds to the lattice site where negative imaginary energy (loss term) is added, while the orange ball indicates the site where positive imaginary energy is included (gain term). The parameters t l and t p refer to the intra- and inter-strand hopping integrals, respectively, and t d describes the diagonal hopping between the two strands.
Figure 1. (Color online). Schematic diagrams of non-Hermitian Aubry–André–Harper ladders following the four different arrangements of the imaginary site energies those are respectively presented in (a), (b), (c) and (d). The green ball corresponds to the lattice site where negative imaginary energy (loss term) is added, while the orange ball indicates the site where positive imaginary energy is included (gain term). The parameters t l and t p refer to the intra- and inter-strand hopping integrals, respectively, and t d describes the diagonal hopping between the two strands.
Materials 15 00597 g001
Figure 2. (Color online). Variation of the real (first column (a,c,e,g)) and imaginary (second column (b,d,f,h)) eigenvalues with the non-Hermitian parameter δ for the four different ladders, where the 1st, 2nd, 3rd and 4th rows correspond to the BT-I, IT-I, IT-II and BT-II cases, respectively. Here, we fix the strength of the Aubry–André–Harper site potential W = 1 in each case. The orange and blue colors are associated with the decoupled chain-1 and chain-2, governed by the relations Equations (10a) and (10b), respectively.
Figure 2. (Color online). Variation of the real (first column (a,c,e,g)) and imaginary (second column (b,d,f,h)) eigenvalues with the non-Hermitian parameter δ for the four different ladders, where the 1st, 2nd, 3rd and 4th rows correspond to the BT-I, IT-I, IT-II and BT-II cases, respectively. Here, we fix the strength of the Aubry–André–Harper site potential W = 1 in each case. The orange and blue colors are associated with the decoupled chain-1 and chain-2, governed by the relations Equations (10a) and (10b), respectively.
Materials 15 00597 g002
Figure 3. (Color online). Dependence of the average inverse participation ratio with the Aubry–André–Harper modulation strength W, where in each spectrum, the black and red lines are associated with chain-1 and chain-2, respectively. The results are shown (al) for three typical δ values (placed in the three different rows), where the four different rows (top to bottom) correspond to the BT-I, IT-I, IT-II and BT-II ladders, respectively. The light-shaded region denotes the mixed-phase window, that is tuned by means of δ .
Figure 3. (Color online). Dependence of the average inverse participation ratio with the Aubry–André–Harper modulation strength W, where in each spectrum, the black and red lines are associated with chain-1 and chain-2, respectively. The results are shown (al) for three typical δ values (placed in the three different rows), where the four different rows (top to bottom) correspond to the BT-I, IT-I, IT-II and BT-II ladders, respectively. The light-shaded region denotes the mixed-phase window, that is tuned by means of δ .
Materials 15 00597 g003
Figure 4. (Color online). Energy resolved inverse participation ratio. Density plot of inverse participation ratios of different eigenstates as functions of the Aubry–André–Harper modulation strength W and the eigenvalues for the BT-I ladder, where the first (a,c,e) and second columns (b,d,f) are associated with the chain-2 and chain-1, respectively. For δ = 0 , we only have real eigenvalues, whereas for non-zero δ , we have both real and imaginary parts, and we plot them into two separate rows (second and third).
Figure 4. (Color online). Energy resolved inverse participation ratio. Density plot of inverse participation ratios of different eigenstates as functions of the Aubry–André–Harper modulation strength W and the eigenvalues for the BT-I ladder, where the first (a,c,e) and second columns (b,d,f) are associated with the chain-2 and chain-1, respectively. For δ = 0 , we only have real eigenvalues, whereas for non-zero δ , we have both real and imaginary parts, and we plot them into two separate rows (second and third).
Materials 15 00597 g004
Figure 5. (Color online). Dependence of two critical points W c 1 and W c 2 on the interaction parameter δ for the BT-I ladder. The green circles and blue squares, associated with the chain-1 and chain-2, are numerically computed, and these data points are fitted with appropriate mathematical expressions which are given in the inset. The two expressions give a very good fit of the numerical data.
Figure 5. (Color online). Dependence of two critical points W c 1 and W c 2 on the interaction parameter δ for the BT-I ladder. The green circles and blue squares, associated with the chain-1 and chain-2, are numerically computed, and these data points are fitted with appropriate mathematical expressions which are given in the inset. The two expressions give a very good fit of the numerical data.
Materials 15 00597 g005
Figure 6. (Color online). Scaling behavior. Variation of the average inverse participation ratio on the number of lattice sites N, both in the absence (left column (ac)) and presence (right column (df)) of the interaction parameter δ for the BT-I ladder. The red color is associated with chain-2, while for the other chain, the black color is used, similarly to in Figure 3. We selectively choose W such that we can only obtain the conducting phase, or the mixed phase or the fully localized one. These three different cases are presented in the top, middle and bottom rows respectively. In each case, we provide a closed-form analytical relation that fits very well with the numerically extracted data.
Figure 6. (Color online). Scaling behavior. Variation of the average inverse participation ratio on the number of lattice sites N, both in the absence (left column (ac)) and presence (right column (df)) of the interaction parameter δ for the BT-I ladder. The red color is associated with chain-2, while for the other chain, the black color is used, similarly to in Figure 3. We selectively choose W such that we can only obtain the conducting phase, or the mixed phase or the fully localized one. These three different cases are presented in the top, middle and bottom rows respectively. In each case, we provide a closed-form analytical relation that fits very well with the numerically extracted data.
Materials 15 00597 g006
Figure 7. (Color online). Dependence of critical points W c m of phase transition on the non-Hermitian parameter δ for different stranded Aubry–André–Harper ladders. We compute the results for the BT-I case only by varying the number of strands M = 1 to M = 10 . The colored dot points represent the critical disorder strengths, which follow quite a regular pattern. Four different δ values are taken into account, as typical examples, and the results are presented in (a), (b), (c) and (d).
Figure 7. (Color online). Dependence of critical points W c m of phase transition on the non-Hermitian parameter δ for different stranded Aubry–André–Harper ladders. We compute the results for the BT-I case only by varying the number of strands M = 1 to M = 10 . The colored dot points represent the critical disorder strengths, which follow quite a regular pattern. Four different δ values are taken into account, as typical examples, and the results are presented in (a), (b), (c) and (d).
Materials 15 00597 g007
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Roy, S.; Maiti, S.K.; Pérez, L.M.; Silva, J.H.O.; Laroze, D. Localization Properties of a Quasiperiodic Ladder under Physical Gain and Loss: Tuning of Critical Points, Mixed-Phase Zone and Mobility Edge. Materials 2022, 15, 597. https://doi.org/10.3390/ma15020597

AMA Style

Roy S, Maiti SK, Pérez LM, Silva JHO, Laroze D. Localization Properties of a Quasiperiodic Ladder under Physical Gain and Loss: Tuning of Critical Points, Mixed-Phase Zone and Mobility Edge. Materials. 2022; 15(2):597. https://doi.org/10.3390/ma15020597

Chicago/Turabian Style

Roy, Souvik, Santanu K. Maiti, Laura M. Pérez, Judith Helena Ojeda Silva, and David Laroze. 2022. "Localization Properties of a Quasiperiodic Ladder under Physical Gain and Loss: Tuning of Critical Points, Mixed-Phase Zone and Mobility Edge" Materials 15, no. 2: 597. https://doi.org/10.3390/ma15020597

APA Style

Roy, S., Maiti, S. K., Pérez, L. M., Silva, J. H. O., & Laroze, D. (2022). Localization Properties of a Quasiperiodic Ladder under Physical Gain and Loss: Tuning of Critical Points, Mixed-Phase Zone and Mobility Edge. Materials, 15(2), 597. https://doi.org/10.3390/ma15020597

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop