EDN: BVFFHK УДК 538.945
Fermion Parity of Phases Supporting Multiple Majorana Modes in a Superconducting Nanowire
Alexander Gamov*
Siberian Federal University Krasnoyarsk, Russian Federation
Anton O. ZlotnikoV
Kirensky Institute of Physics Federal Research Center KSC SB RAS Krasnoyarsk, Russian Federation Siberian Federal University Krasnoyarsk, Russian Federation
Received 10.09.2023, received in revised form 25.10.2023, accepted 28.11.2023 Abstract. The fermion parity of the ground state is determined in various topological phases of the semiconducting nanowire under external magnetic field with proximity-induced superconductivity and strong spin-orbit interaction. Electron hopping as well as spin-flip hopping due to spin-orbit coupling and superconducting pairings in the second coordination sphere are taken into account. The connection between the fermion parity and the parity of the BDI topological invariant is shown. The formation of topological phases with three and four pairs of Majorana modes has been demonstrated. Keywords: Fermion parity, superconducting nanowire, Majorana modes, topological invariant, topo-logical phase diagram.
Citation: A. Gamov, A.O. Zlotnikov, Fermion Parity of Phases Supporting Multiple Majorana Modes in a Superconducting Nanowire, J. Sib. Fed. Univ. Math. Phys., 2023, 16(6), 820-829. EDN: BVFFHK.
Majorana modes (MMs) that emerge in topological superconductors are of both fundamental and practical interest. This is due to their spatial nonlocality and the potential for implementing non-Abelian statistics through permutations. These characteristics suggest the possibility of utilizing MMs for quantum computing, providing protection against decoherence processes [1-4]. Permutations of non-Abelian anyons result in alterations to the system's ground state wave function. Unlike Abelian anyons (such as fermions and bosons), this change is not solely confined to the introduction of a phase factor.
One of the most renowned systems where the realization of MMs is feasible involves an InSb (or InAs) semiconductor nanowire brought into contact with an s-wave superconductor. For instance, this could be achieved by depositing a 3-5 mm thick layer of aluminum on the nanowire's surface and subjecting it to an external magnetic field [5,6]. The formation of MMs in such heterostructures was initially predicted in several studies [7-9]. The joint implementation of strong spin-orbit coupling, spin-singlet s-wave superconductivity, and a magnetic field allows the realization of a Kitaev chain regime with effective triplet p-wave superconducting pairings [10]
* agamow96@gmail.com https://orcid.org/0009-0000-3987-4248 tzlotn@iph.krasn.ru https://orcid.org/0000-0003-3195-1397 © Siberian Federal University. All rights reserved
in the system under consideration. It's essential to note that the experiments have yet to conclusively confirm that the observed signals, such as the conductance peak at zero bias voltage, are precisely associated with MMs and not with trivial Andreev bound states [11,12].
Typically, in superconducting nanowires, a Majorana bound state is considered, consisting of two MMs localized at opposite ends of the wire. However, the one-dimensional (1D) Hamiltonian used to describe superconducting wires belongs to the BDI symmetry class [13-15], enabling not only non-trivial topology in 1D systems but also the formation of multiple MMs on a single end of a finite-dimensional chain. In this scenario, different MMs on one edge of the chain do not hybridize with each other due to the effective symmetry similar to time-reversal symmetry [16]. Indeed, the realization of a phase with two pairs of MMs (two modes at each edge of the chain) was shown in Refs. [16,17] when considering extended s-wave superconducting pairings. This implies the presence of pairing interaction between electrons at one site of the chain, as well as between nearest sites. The transition to such a phase can be induced by on-site Hubbard repulsion of electrons [17].
Previously, the formation of topological phases with multiple edge MMs was also shown in chains containing magnetic atoms [18] and under the influence of periodically time-varying external fields [19]. The interest in topological phases with several MMs stems from their realization with arbitrarily small Zeeman splitting of spin subbands (in zero field, they transform into Majorana-Kramers doublets [16]). Additionally, it offers the potential to permute multiple MMs simultaneously in quantum wires of T-, X-, and Y-structures, as was shown for single MMs in [2].
It is natural to assume that in the mentioned systems there also exist non-zero amplitudes of electron hopping, including spin-flip hopping due to spin-orbit coupling, and superconducting pairings between next-to-nearest neighbors in the chain. In this work, considering the processes in the second coordination sphere (CS), we demonstrate the formation of phases with three and four pairs of MMs in the model of a superconductor-semiconductor nanowire with spin-orbit coupling, placed in a magnetic field. Notably, for the classification of regions in the topological phase diagrams, a clear method is employed in addition to the BDI topological invariant [14,16]. This method is associated with determining the fermion parity (FP) of the system ground state and was previously proposed in Refs. [20-22].
1. Model of a superconducting nanowire in the presence of processes in the second coordination sphere
The Hamiltonian of the model of a superconducting nanowire with open boundary conditions in the tight binding approximation has the form [21]:
N N-1
HW = £ — - h)CL Cna + £ t(C™ Cn+1a + Cl+ia Cna ) + a,n=1 a,n=1
+ £ {(-*)(CUCn+n - C^Cn+it) + H.c] + £ [AcCitCl + H.c] + (1)
n=1 n=1
N-1
+ £ [A(CntCn +14 + Cn+1tC^)+ H.c].
The first term in (1) characterizes the electron energy measured from the chemical potential M taking into account the Zeeman splitting h in an external magnetic field. The second term describes hoppings between nearest sites with amplitude t. The third and fourth terms define the Rashba spin-orbit interaction with amplitude a and on-site superconducting pairings with amplitude Ao, respectively. The last term is due to superconducting pairings between electrons at nearby sites with the amplitude A. It was shown in [16,17] that the formation of phases with two pairs of MMs is possible at 2A > A0. To consider this regime in further calculations, it is convenient to set A0 = 0, keeping in mind that the results will be preserved at a qualitative level at A0 =0 and 2A > A0. The parameter no determines the sign and takes the following values depending on the projection of the spin momentum a: no = +1 (at a =t), no = —1 (at a =4-).
The Hamiltonian (1) describes only processes in the first CS. Therefore, to take into account long-range hopping, including spin-flip processes, as well as superconducting pairings in the second CS, we add the following terms:
N-2
N-2
H2 = J2 t2(CL Cn+2o + Cl+2o Cno ) + £ [(-a2)(CltCn+2l — Cn+2t) + Hc
o,n=l n=1
N-2
+ £ [A2(CltCl +24 + Cl+2tC4.)+ H.c. n=l
+
(2)
The first part of (2) determines the hopping between next nearest sites with amplitude t2, the second and third terms are responsible for the effective spin-orbit interaction with the parameter a2 and superconducting pairing with the parameter A2 for the second CS. As a result, the Hamiltonian of the model under consideration is written as:
H = Hw + H2
(3)
2. Fermion parity and topological phase diagram
The Hamiltonian of the model in the reciprocal space under periodic boundary conditions takes the form:
h = £ 6o CI Cko + £[/
Ak C%C-
k4 + Afc C-kiCl
where
k4Ckt
+ y^ iak k
CktCk4 - Ck4Ckt
(4)
£ko
tk — M — noh, tk = 2t cos(k) + 2t2 cos(2k),
ak = 2asin(k) + 2a2 sin(2k), Ak = 2Acos(k) + 2A2 cos(2k).
(5)
(6) (7)
It was previously shown that the fermion parity (FP) of the ground state of a superconducting system is related to the structure of the Hamiltonian (1) at special points of the first Brillouin zone [20-22]. Here, the special points are considered as those values of the quasimomentum at which either the superconducting order parameter Ak or the spin-orbit coupling integral ak vanishes while having non-zero amplitudes a, A, a2, A2. The identification of such special points is justified by the fact that, as it will be shown later, the Bogoliubov transformations have a different form for these points compared to other BZ points. It is worth noting that for models,
like Bardeen, Cooper, Schrieffer (BCS) model, special points are limited to the nodal points of the superconducting order parameter, where Ak = 0. When considering the spin-orbit interaction, the energy spectrum takes a form distinct from BCS theory, although the superconducting gap in the spectrum still closes at the nodal points upon tuning model parameters. However, it can be easily shown in the considered model, that the condition ak =0 also results in the closing the gap in the elementary excitation spectrum, which could be the cause of a topological transition with a possible change in the ground state FP. In this case, the gap is closed at points where ak =0 but the superconducting order parameter Ak is non-zero. It is important to emphasize that it applies specifically to the closure of the superconducting gap, as in the absence of superconductivity the spectrum becomes gapless on the Fermi contour (lines in BZ).
From (6) and (7) it is obtained that the parameter Ak vanishes at the following points:
K± = * - arccos (A±) , -K±; A± = 1A ± ^(A^ +8 G -1, ^ (8)
and the spin-orbit coupling integral ak vanishes at the points
a
Q = 0, n, P = n - arccos (a), -P; a =- e [-1,1]. (9)
2a2
Thus, the Hamiltonian (4) is modified to the following sum of quadratic forms:
H = h(0) + h(n) + r(P)+ r(K+) + r(K-)+ E r(k), (10)
k>0,k=Q, P,K±
where r(k) = h(k) + h(-k) (this definition is also valid for points P, K±), and
h(Q) = £ Q CQ a CQa + AQCQt CQl + A*QCQiCQt, (11)
a
h(P) = £ £p* Cp a C pa + Ap CpfC- Pi + AP C-PiCpt, (12)
a
h(K± ) = £K±a CK ± a CK±a + iaK± CK± tCK±i - CK±iCK±t , (13)
h(k) = E Cka + Ak CltC- + Ak + iak
Cl\CH - CiiCk\
(14)
In [22], it was demonstrated that the eigenfunction of the term r(k) (refer to (14)) can always be expressed as a superposition of states with an even number of fermions, akin to the wave function in BCS theory. Consequently, the structure and FP of the system's ground state are determined by the eigenvectors of the operators h(Q), r(P), and r(K±).
As delineated in [22], the eigenvectors of the term of h(Q) comprise two states: one with a single fermion, denoted as \FQa) = CQa|0), and two with an even number of fermions, represented as \$q±) = [aQ + B±CQtCQ\0). Of particular interest is the scenario where the state containing one fermion possesses the lowest energy. For h > 0 (which is the case always considered here), this reduces to the condition EQ^ < EQ _. Consequently, the following inequalities emerge:
h > J{2t2 + 2t - n}2 + \2A2 + 2A\2,
,- (15)
h>^ {2t2 - 2t - n}2 + \2A2 - 2A\2.
The term r(P) has the BCS form, and its eigenvectors also consist of an even number of fermions, leaving FP of the system's ground state wave function unaffected. In contrast, the term h(K±) possesses an eigenvector that can be expressed as a superposition of states with a single fermion: (aK±CK±o + bK±CK±o^j |0) (a represents the opposite spin projection to a).
Nevertheless, it's evident that the state (a-K±C— K±o + b-K±C-K±o^j |0), which is appropriate for the term h(—K±), has the same energy. Therefore, even in cases where these states have lower energies than states with an even number of fermions, FP of the ground state wave function remains unchanged because the one-fermion states at points K± and —K± are filled simultaneously. Consequently, the conditions for achieving the ground state of a system with odd FP are dictated by the inequalities (15).
It is important to note that using the quadratic forms highlighted in (10), it is feasible to describe transitions that change the ground state wave function even when FP remains the same. By applying the Bogoliubov transformation, we can determine the excitation energies at the selected k points that may change signs with parameter variations.
Elk = J(tk — m)2 + |A2k | — h, k = Q,P, —P,
V ,- (16)
E2,3k = (—M) ±Va2k + h2, k = K±, —K±.
From the considerations given above, it follows that FP of the ground state for a given set of parameters coincides with the parity of the number of negative energies (16). The energies E1q at points 0 and n assume negative values under the conditions of (15), respectively. We will now outline the inequalities that lead to the remaining energies becoming negative:
Eip, —p < 0, h>
E2K±, — K± < 0, h<
E3K± — K± < 0, h>
E3K±, — K± < 0, M>
2 2 ^ 1/2 a2 — 1/2) — 2at — m] + |4A2(a2 — 1/2) — 2aA| j , (17)
■*2}1/2
2 } 1/2 ^ j
M > Mc, (18)
M < —Mc, (19)
(20)
where Mc = 2^(1 — A±)(a — 2a2^±).
Whenever the sign of any of the energies changes, it leads to a modification in the ground state wave function. It is easy to verify that the same conditions correspond to the closure of the gap in the elementary excitations spectrum at special points such as Q, ±P, K±, —K±, while at other points the spectrum remains gapped until the amplitudes of superconducting pairings are non-zero.
In further calculations, we will assume that all parameters within the same CS are equal to each other. In this context, we will consider A1 = a = A = t and A2 = a2 = A2 = t2 and explore the case where A2 = 1.5A1, as it allows for the maximum number of MMs to be realized within the given model. Moreover, to establish the relationship between FP and the parity of the topological invariant, we compute the value of the topological index BDI [14,16]:
—i fk= dz(k) det(Q(k)) , , ..
Nbdi =— —7jt-, z(k) =.,.,!, |, Q(k) = a0 — haz — akay — lAkay, (21) n Jk=0 z(k) |det(Q(k))|
where a0 is the identity matrix, ay, az are Pauli matrices.
0 m-1--1---1-™
-5 0 5
P
Fig. 1. Topological phase diagram in the variables chemical potential n — Zeeman splitting h (in units of Ai) for the amplitude A2 = 1.5Ai. The lines are obtained from the conditions of changing the sign of the energies of one-fermion states at special points of the Brillouin zone. The value of the topological index is marked in each phase. Different colors indicate different fermion parities of the ground state wave function
Fig. 1 illustrates parameter lines h(^) where the energy sign changes (16). As indicated in Fig. 1, these lines delineate the regions with distinct values of the topological index. Additionally, these regions are characterized by distinct numbers of negative energies (16). The solid line depicts the dependence (15). Traversing this line results in an energy sign change and a parity shifts in the count of negative energies. The other dependencies, determined from the inequalities (17)-(19), are shown by the dashed lines in Fig. 1. The dashed line obtained from (17) is highlighted in purple, the dependencies (18) and (19) are shown by green and red lines, respectively (see color version). Crossing these lines does not change FP but does modify the form of the ground state wave function. Similarly, when passing through the dashed lines, the topological index changes without affecting its parity. Consequently, the solid line serves as the boundary of the regions in the topological phase diagram with different FP of the wave function of the ground state. The regions highlighted in grey (turquoise in color version) have even FP, and the light (yellow) ones have odd FP.
It can be concluded that the parity of the topological invariant coincides with the parity of the number of one-fermion states that have negative energy in a given region. In this work, we do not specify the precise forms of the ground state wave functions. Their construction is similar to the method described in [22].
It is seen in Fig. 1, that the phases with NBDI = 3 and NBDI = 4 are formed in the presence of processes in the second CS. It should be noted that the NBDI = 3 phase appear even when A2 < A1 (approximately for A2 > 0.5A1). And this phase can be realized at low values of Zeeman splitting 0.1t). On the other hand, the NBDI = 4 phase is implemented for A2 > A1.
3. Majorana modes in phases with NBDI = 3 and NBDI = 4
To describe features of MMs in the new phases, we consider boundary conditions and assume that the number of sites is N original Fermi operators to Bogoliubov operators is carried out:
N , ,
aj = E UjufCnt + UjniCnl + VjniCni + vjnfCnt
n=1 N
a] = E ujn-iCU + uUCU + vjntCnt + vUCm
n=1
the Hamiltonian (3) with open = 150. The transition from the
j = [-2N, 2N] .
The coefficients ujna and vjna in this case are components of the eigenfunctions of the Hamil-tonian matrix. They can be used to obtain the distribution of quasiparticle operators over the sites of the chain. However in order to spatially separate MMs in the chain, it is convenient to use Majorana operators:
b = aj + a\ = Y, Una Y Ana + Zna YBna, n,a
b" = i (a\ - aA =Y1 Una Y Ana + Zna YB
^ ' n a
(23)
i^na 2Re(ujnCT + vjna )? zna 2^-m(ujna vjna )?
OJna 2-m(ujna + vjna)? zna 2Re(ujna vjna)?
where YAna = Cna + C^a, YBna = ¿(CL — Cna) are Majorana self-conjugate operators.
Fig. 2. Distributions of Majorana modes depending on the chain site for the first (a), second (b), third (c) and fourth (d) excitation energies in the phase with NBDI = 3. The ratio between amplitudes in the second and first coordination spheres is A2/A1 = 2
In Fig. 2, the distributions of MMs on the chain sites (coefficients uno and ujno in (23)) are presented for the phase with NBDI = 3. The dependencies for coefficients zno and zno are qualitatively the same. Figures (a-c) reveal that the formation of MMs with edge localization occurs only for the first three excitation energies. These three excitations have near-zero energies (counted from the Fermi level). In the case of the fourth excitation with nonzero energy (d), the probability density of quasiparticles is minimal at the edges of the chain and increases towards to its middle. The probability density of edge states in NBDI = 3 phase exhibits an exponential decay with oscillations.
Fig. 3 shows the distributions of MMs for the phase with NBDI = 4. As evident from figures (a-d), MMs are formed for all four excitation energies in this phase, resulting in the formation of four pairs of MMs with near-zero energies.
Fig. 3. Site distributions of Majorana modes for the first (a), second (b), third (c) and fourth (d) excitation energies in the phase with NBDI = 4 (\2/\\ = 2)
СС^^ЮП
The model of the superconducting nanowire under external magnetic field is considered with the periodic as well as open boundary conditions taking into account electron hoppings, including spin-flip hoppings due to Rashba spin-orbit coupling, as well as superconducting pairings between nearest and next-nearest neighbors. It is shown for the periodic boundary conditions that different phases on the topological phase diagram can be characterized by the number of occupied one-fermion states that have negative energy at special points of the Brillouin zone. As a result, the region of formation of odd fermion parity of the ground state is determined. It is obtained that the parity of the ground states coincides with the parity of the topological invariant, which in turn indicates the number of Majorana bound states, each of which contains one pair of Majorana modes localized at opposite edges of the chain. We show that taking into account the processes between next-nearest neighbors in a chain with open boundary conditions leads to the formation of phases with three and four pairs of Majorana modes.
The work was partially carried out within the state assignment of Kirensky Institute of Physics. Gamov A. acknowledges the support of the Theoretical Physics and Mathematics Advancement Foundation "BASIS".
References
[1] C.Nayak, S.H.Simon, A.Stern, M.Freedman, S.Das Sarma, Non-Abelian Anyons and Topological Quantum Computation, Rev. Mod. Phys, 80(2008), 1083-1159.
[2] J.Alicea, Y.Oreg, G.Refael, F. von Oppen, M.P.A.Fishser, Non-Abelian statistics and topological quantum information processing in 1D wire networks, Nature Physics, 7(2011), 412-417. DOI: 10.1038/nphys1915
[3] T.Karzig, F.Pientka, G.Rafael, F.von Oppen, Shortcuts of non-Abelian braidings, Phys. Rev. B, 91(2015), 201102. DOI: 10.1103/PhysRevB.91.201102
[4] F.Harper, A Pushp, R.Roy, Majorana braiding in realistic nanowire Y-junctions and tuning forks, Phys. Rev. Res., 1(2019), 033207. DOI: 10.1103/PhysRevResearch.1.033207
[5] M.T.Deng, S.Vaitiekenas, E.B.Hansen, J.Danon, M.Leijnse, K.Flensberg, J.Nygard, P.Krogstrup, C.M. Marcus, Majorana bound state in a coupled quantum-dot hybrid-nanowire system, Science, 354(2016), 1557-1562. DOI: 10.1126/science.aaf3961
[6] Microsoft Quantum et al., InAs-Al Hybrid Devices Passing the Topological Gap Protocol, Phys. Rev. B, 107(2023), 245423. DOI 10.1103/physrevb.107.245423
[7] J.D.Sau, R.M.Lutchyn, S.Tewari, S.D.Sarma, Generic New Platform for Topological Quantum Computation Using Semiconductor Heterostructures, Phys. Rev. Lett., 104(2010), 040502. DOI: 10.1103/PhysRevLett.104.040502
[8] R.M.Lutchyn, J.D.Sau, S.Das Sarma, Majorana Fermions and a Topological Phase Transition in Semiconductor-Superconductor Heterostructures, Phys. Rev. Lett., 105(2010), 077001. DOI: 10.1103/PhysRevLett.105.077001
[9] Y.Oreg, G.Refael, F. von Oppen, Helical Liquids and Majorana Bound States in Quantum Wires, Phys. Rev. Lett., 105(2010), 177002. DOI:10.1103/PhysRevLett.105.177002
[10] A.Yu.Kitaev, Unpaired Majorana fermions in quantum wires, Phys. Usp., 44(2001), 131-136.
[11] C.-X.Liu, J.D.Sau, T.D.Stanescu, S.D.Sarma, Andreev bound states versus Majorana bound states in quantum dot-nanowire-superconductor hybrid structures: Trivial versus topological zero-bias conductance peaks, Phys. Rev. B, 96(2017), 075161.
DOI: 10.1103/PhysRevB.96.075161
[12] S.Das Sarma, In Search of Majorana, Nat. Phys., 19(2023), 165-170. DOI: 10.1038/s41567-022-01900-9
[13] S.Ryu, A.P.Schnyder, A.Furusaki, A.W.W.Ludwig, Topological Insulators and Superconductors: Tenfold Way and Dimensional Hierarchy, New J. Phys., 12(2010), 065010.
DOI: 10.1088/1367-2630/12/6/065010
[14] S.Tewari, J.D.Sau, Topological Invariants for Spin-Orbit Coupled Superconductor Nanowires, Phys. Rev. Lett., 109(2012), 150408. DOI: 10.1103/PhysRevLett.109.150408
[15] Y. Niu, S.B.Chung, C.-H.Hsu, I.Mandal, S.Raghu, S.Chakravarty, Majorana zero modes in a quantum Ising chain with longer ranged interactions, Phys. Rev. B, 85(2012), 035110. DOI: 10.1103/PhysRevB.85.035110
[16] C.L.M. Wong, K.T.Law, Majorana Kramers Doublets in dx2-y2-Wave Superconductors with Rashba Spin-Orbit Coupling, Phys. Rev. B, 86(2012), 184516.
DOI: 10.1103/PhysRevB.86.184516
[17] S.V.Aksenov, A.O.Zlotnikov, M.S.Shustin, Strong Coulomb interactions in the problem of Majorana modes in a wire of the non-trivial topological class BDI, Phys. Rev. B, 101(2020), 125431. DOI: 10.1103/PhysRevB.101.125431
[18] A.A.Bespalov, Tuning the Topological state of a helical atom chain via a Josepson phase, Phys. Rev. B, 106(2022), 134503. DOI: 10.1103/PhysRevB.106.134503
[19] H.Wu, S.Wu, L.Zhou, Floquet Topological Superconductors with Many Majorana Edge Modes: Topological Invariants, Entanglement Spectrum and Bulk-Edge Correspondence, New J. Phys., 25(2023), 083042. DOI: 10.1088/1367-2630/acf0e3
[20] N.Read, D.Green, Paired States of Fermions in Two Dimensions with Breaking of Parity and Time-Reversal Symmetries and the Fractional Quantum Hall Effect, Phys. Rev. B, 61(2000), 10267. DOI: 10.1103/PhysRevB.61.10267
[21] V.V.Val'kov, M.S.Shustin, S.V.Aksenov, A.O.Zlotnikov, A.D.Fedoseev, V.A.Mitskan, M.Yu.Kagan, Topological superconductivity and Majorana states in low-dimensional systems, Phys. Usp., 65(2022), 2-39. DOI 10.3367/UFNe.2021.03.038950
[22] V.V.Val'kov, V.A.Mitskan, M.S.Shustin, Ground-State Fermion Parity and Caloric Properties of a Superconducting Nanowire, J. Exp. Theor. Phys., 129(2019), 426-437.
DOI: 10.1134/S1063776119080144
Фермионная четность в фазах с множеством майорановских мод в сверхпроводящей нанопроволоке
Александр Гамов
Сибирский федеральный университет Красноярск, Российская Федерация
Антон О. Злотников
Kirensky Institute of Physics Федеральный исследовательский центр СО РАН Красноярск, Российская Федерация Сибирский федеральный университет Красноярск, Российская Федерация
Аннотация. Для полупроводниковой нанопроволоки, помещенной во внешнее магнитное поле, с наведенной сверхпроводимостью и сильным спин-орбитальным взаимодействием определена фермионная четность основного состояния в различных топологических фазах при учете перескоков электронов, включая перескоки с переворотом спина за счет спин-орбитальной связи, и сверхпроводящих спариваний во второй координационной сфере. Показана связь фермионной четности и четности BDI топологического инварианта. Продемонстрировано формирование топологических фаз с тремя и четырьмя парами майорановских мод.
Ключевые слова: фермионная четность, сверхпроводящая нанопроволока, майорановские моды, топологический инвариант, топологическая фазовая диаграмма.