Real fermion modes, impurity entropy, and nontrivial fixed points in the phase diagram of junctions of interacting quantum wires and topological superconductors

We discuss how to extend the impurity entropy to systems with boundary interactions depending on zero-mode real fermion operators (Majorana modes as well as Klein factors). As specific applications of our method, we consider a junction between N interacting quantum wires and a topological superconductor, as well as a Y-junction of three spinless interacting quantum wires. In addition we find a remarkable correspondence between the N=2 topological superconductor junction and the Y-junction. On one hand, this allows us to determine the range of the system parameters in which a stable phase of the N=2 junction is realized as a nontrivial, finite-coupling fixed point corresponding to the M-fixed point in the phase diagram of the Y-junction. On the other hand, it enables us to show the occurence of a novel ``planar'' finite-coupling fixed point in the phase diagram of the Y-junction. Eventually, we discuss how to set the system's parameters to realize the correspondence.


Introduction
Recently, considerable interest has arisen in junctions involving interacting quantum wires (QW's), both in the case of spinless [1,2,3,4,5,6], and of spinful wires [7,8,9,10]. This is mainly due to the fact that Landau's Fermi liquid paradigm typically breaks down in onedimensional interacting electronic systems, whose low-energy, long-wavelength properties are rather described by means of the Tomonaga-Luttinger liquid (TLL) framework [11,12]. Within the TLL-approach, tunneling processes at a junction are described in terms of nonlinear vertex operators of the bosonic fields, with nonuniversal scaling dimensions continuously depending on the "bulk" interaction parameters [13,14]. This opens the way to a plethora of nonperturbative features in the phase diagram of those systems, including the remarkable emergence of intermediate, finite-coupling fixed points (FCFP's), either describing phase transitions between different phases (repulsive fixed points), or novel, nontrivial phases of the junction (attractive fixed points), thus generalizing to multi-wire junctions the Kane-Fisher FCFP emerging at a junction between two spinful QW's [15].
In this context, the prediction that localized Majorana modes (MM's) can appear at a junction between a normal QW and a topological superconductor (TS) [16] has opened additional brandnew scenarios, as the direct coupling between a quantum wire and a localized MM can potentially give rise to relevant boundary interactions typically, not allowed at junctions between normal wires [17]. As a result, it has been possible to predict the emergence of novel FCFP's in the phase diagram of junctions between more-than-one interacting QW and TS's [18,19]. Moreover, due to ubiquity of the TLL-formalism, which successfully describes (junctions of) quantum spin chains [20,21,22,23], Josephson junction networks [24,25,26,27], as well as topological, Kondo-like systems [28,29,30,31], novel FCFP's have been predicted to emerge in the phase diagram of those systems, as well. Besides their theoretical interest, FCFP's have been argued to correspond to "decoherence-frustrated" phases, in which competing frustration effects can operate to reduce the unavoidable decoherence in the boundary quantum degrees of freedom coupled to the "bath" of bulk modes [32,33], thus making the junction, regarded as a localized quantum impurity, a good candidate to work as a frustration-protected quantum bit [24]. For this reason, it becomes of importance to search for FCFP's in the phase diagram of pertinently designed junctions of quantum wires.
An effective means to study the phase diagram of junctions of QW's is given by a cooperative combination of perturbative renormalization group (RG) approach [34,15,35] and of the delayed evaluation of boundary conditions (DEBC) technique [2,7]. In particular, as extensively discussed in e.g. Ref. [2], DEBC method is based on constructing the boundary operators allowed by symmetries at a certain fixed point: the emergence of (at least one) relevant operator is, therefore, evidence of the instability of that fixed point against some other one. In addition, as pointed out, for instance, in Ref. [2], for a junction of normal QW's and in Refs. [18,19] for junctions between interacting QW's and TS's, those methods are efficiently complemented by using the zero-temperature impurity entropy (IE) to characterize the fixed points of the junction. The IE was originally introduced as a mean to characterize and classify, at zero temperature, quantum impurity systems that are critical in the bulk at the fixed points of their boundary phase diagram [35,36]. Later on, it has been shown to correspond to the impurity contribution to the groundstate entanglement entropy, which is particularly suitable for characterizing the phases of 1+1-dimensional models via the Density Matrix Renormalization Group (DMRG) technique [37]. The exponentiated IE g yields the groundstate degeneracy of the system at a certain fixed point. In general, g corresponds to the specific value of the "g-function", which always decreases along renormalization group trajectories ("g-theorem"); therefore, if two boundary fixed points are characterized by groundstate degeneracies g 1 and g 2 , with g 1 > g 2 , provided there is the possibility of connecting the two of them with a renormalization group trajectory, this will always flow from g 1 to g 2 . [35,36,38,39,40].
Technically, to compute g one first constructs the partition function over a strip of length ℓ for the model at a boundary critical point. Letting A label the conformally invariant boundary conditions (CIBC's) characterizing a specific boundary critical point, one computes the partition function on the strips, Z AA , by assuming type A-CIBC's at both boundaries. Specifically, one gets Z AA = n exp − x n AA βu ℓ , where the sum is taken over energy eigenstates of the whole system, u is a velocity scale determined by the (critical) bulk of the system, β is the inverse temperature, and x n AA are dimensionless numbers typical of the system. The g-function is derived by sending ℓ → ∞ at finite β. In this limit, one obtains Z AA → g 2 A e πℓc 6βu , with the dimensionless number c being the conformal anomaly of the bulk critical theory. From the last result, g A can be readily extracted [41].
In this paper we employ a combined use of perturbative RG approach, DEBC-method and calculation of the g-function to study junctions of interacting QW's and TS's and to spell out the correspondence between this sort of junctions and the Y-junction of three interacting quantum wires (Y3J) studied in Ref. [2] in the bulk Z 3 -symmetric version and later on discussed in Ref. [42] in the general case in which the Z 3 symmetry is broken. In doing so, we have necessarily to take into account the emergence of real fermionic modes at our junctions. These are primarily provided by the localized MM's {γ j } emerging at the interface between a topological superconductor and a normal system [16,43,44]. In addition, real fermionic operators also appear as Klein factors (KF's), {Γ j }, which have to be introduced when employing the bosonization approach to interacting one dimensional fermionic systems, to recover the correct (anti)commutation relations between fermionic fields for different wires, as well as between the fermionic field for each wire and the γ j 's. In many cases, KF's play no role, as the relevant multi-point correlation functions of the fermionic fields either contain an equal number of creation and annihilation operators of fermions of the same kind (see, for instance, Ref. [45] and reference therein), or, at most, the net effect of KF's can be an extra minus sign which can be equally well accounted for by, for instance, redefining the zero-mode operators of the bosonic fields [46]. On the other hand, they are definitely essential to recover the correct phase diagram of e.g. a junction of three interacting QW's [2], as well as to correctly account for the hybridization between MM's and normal electronic modes in a conductor [47,48], at junctions between normal wires and TS's [17,18,49], and in the remarkable "topological" realization of the Kondo effect, in which MM's determine an effective impurity spin coupled to electronic modes by the normal contacts [50]. Thus, it is by now evident that they must be properly accounted for in applying the g-theorem to junctions of interacting quantum wires.
Besides generalizing the results of Refs. [18,2] respectively to a multiwire junction of interacting QW's and a TS's and to a non Z 3 -symmetric Y3J, we unveil the remarkable correspondence between the two models. In doing so, we prove how, extending the range of system's parameters with respect to the case discussed in Ref. [18], it is possible for a FCFP to correspond to the stable phase of the system. Remarkably, this completely reverses the scenario found within the range of parameters discussed in Ref. [18], where we proved that it is necessary to fine-tune the boundary couplings to the MM, to drive the RG flow towards the FCFP, which was unstable against more "trivial" fixed points. Apart being interesting per se, this result appears also of relevance for engineering stable phases with frustrated decoherence, potentially amenable for applications 3 to quantum computation. An important point to stress about the correspondence is that it holds despite the obvious observation that the Y3J has three KF's, while the N = 2 junction has just two KF's and one MM. In fact, as we discuss in the following, while one of the three fields of the Y3J decouples from the boundary interaction, its KF gets "left behind" and it plays the role of the MM. Therefore, the correspondence works perfectly well when considering a single junction in both systems. Remarkably, it also yields the right result when computing g at the fixed points of the two model. In this case, as we discuss in detail below, one has to resort to a two-boundary version of the corresponding model Hamiltonian. Thus, one cannot ignore the intrinsic difference between MM's, which can be assumed to be basically local in real space, and KF's, which on the contrary are global, as one naturally associates them with the whole extent of a QW in real space. Despite this, the correspondence works fine and allows for recovering a number of nontrivial results about the phase diagram of one model from what is known about the phase diagram of the other.
The paper is organized as follows: • In section 2 we introduce our procedure for computing the g-function in boundary models with real fermionic modes emerging in the boundary interaction Hamiltonian. To illustrate our procedure, here we apply it to a single interacting spinless quantum wire connected to two p-wave superconductors, by discussing in detail the subtleties in counting the degrees of freedom associated with the real fermionic modes and how to deal with them; • In section 3 we discuss the main features of the phase diagram and compute the impurity entropy at the fixed points of a junction between two quantum wires and a topological superconductor and of the asymmetric Y junction of three spinless interacting quantum wires.
In both cases, we mostly review known results [18,2,42] which, nevertheless, are important for the sake of the presentation of the following results; • In section 4 we discuss in detail the correspondence between a junction with two quantum wires and a topological superconductor and the asymmetric Y3J. In particular, we show how the results derived in Ref. [18] for the former system can shed light on the phase diagram of the Y3J in the case of asymmetric bulk, as well as boundary interaction and, conversely, how the results of Refs. [2,42] for the Y3J allow for extending the analysis of the phase diagram of the junction with two quantum wires and a topological superconductor to windows of values of the system's parameters which were not encompassed in the derivation of Ref. [18]; • In section 5, as a further application of our method for computing the g-function in boundary models with real fermion modes emerging at the boundary interaction, we generalize the results of Ref. [18] by discussing the fixed points in the phase diagram, and the corresponding calculation of the g-function, in a junction between N quantum wires and a topological superconductor.
• In section 6 we provide our conclusions and discuss possible further developments of our work.
• In the various appendices, we provide mathematical details of our derivation. To help following the various abbreviations, we list in table 1 the meaning of the ones we use most commonly throughout the paper.

Impurity entropy in a boundary model with real fermionic modes in the boundary interaction
When bosonizing more than one species of fermion operators in one dimension, real fermionic Klein factors must be introduced, to properly account for the anticommutation relations between operators corresponding to different species of fermions. Typically, KF's appear in boundary Hamiltonians describing junctions of one-dimensional quantum wires (which is appropriate, at points where different wires contact each other) and, in many cases, they strongly affect the boundary dynamics of the junction [45]. For instance, only by properly accounting for KF's in the boundary Hamiltonian, does one prove the emergence of a FCFP in the Y3J discussed in Refs. [1,2], or in its spinful version [7]. In addition to KF's, real fermion operators emerge as MM's at junctions between QW's and TS's [16]. The combined effect of KF's and MM's can eventually lead to the "Majorana-Klein" hybridization and, eventually, to a remarkable topological version of the Kondo effect [50,28,30], also discussed in its multi-channel realization [51,31,29], as well as to novel phases corresponding to FCFP's in the phase diagram of junctions between QW's and TS's [18,19]. Thus, despite their definition as a mathematical means for properly doing bosonization, KF's affect the boundary dynamics of a junction exactly as "physical" MM's do and, accordingly, they must be properly accounted for, when computing the IE of the junction. To demonstrate this point, in this section we compute the IE in a paradigmatic system given by a single spinless interacting wire connected to two TS's at its endpoints. This enables us to show how, in order to find results for the g-function consistent with the expected phase diagram of this systems, one has to count the degrees of freedom associated with real fermions at the system boundaries. In doing so, we face an additional subtelty, which was originally put forward in Ref. [52], which is strictly related to how to count real, zero-mode fermionic degrees of freedom.
In general, two real fermionic modes, say γ a and γ b , can be combined together into a complex (Dirac) fermionic mode a = 1 2 (γ a + iγ b ), which leads to a single fermion energy level, which can be either empty, or full, eventually resulting in an additional degeneracy factor of 2 in the 5 partition function. In a boundary theory, the procedure for computing the partition function (and, eventually extracting the g-function from the result) consists in making up a two-boundary version of the model Hamiltonian by mirroring the boundary interaction describing the junction at the other boundary of a finite-size (ℓ) version of the system (see Ref. [53] for details of the procedure). While this procedure unavoidably leads to a doubling of the MM's emerging at the interfaces (so they always contribute an even number of real fermionic modes) [17,18,49,19], when introducing KF's through bosonization of the normal wires, the final total number of real fermionic modes can either be even, or odd. When it is odd, one has to face an ambiguity about how to count the leftover real fermion, which is strictly related to the need to account for fermion parity conservation in the presence of real fermion operators [52]. To overcome such a difficulty, we introduce an additional decoupled "auxiliary" fermionic wire which, in bosonization language, is characterized by its own Luttinger parameterK and by its Klein factorΓ.Γ enters the total counting of degrees of freedom related to real fermions, by pairing with the real mode that is left over after all the other "physical" ones have paired into Dirac complex modes. Of course, we expect this to affect the actual value of the g function. Yet, as for any definition of entropy, what matters in the g-theorem is the entropy difference between two different fixed points or, which is the same, the ratio between the corresponding values of g. In fact, we expect our procedure to provide the correct result for the ratio and, to ground our speculation, in the following we provide a number of different examples of physical interest where we show that this is, in fact, the case. To illustrate our procedure, here we apply it to compute the g-function in a single interacting spinless quantum wire connected to two p-wave superconductors in their topological phase [16] at its endpoints. At low energies, the superconducting leads can be traded for two MM's γ L , γ R , respectively residing at the left-hand side and at the right-hand side of the QW. In addition, the interacting QW is effectively described by resorting to the bosonization approach, which we review in Appendix A, in terms of the bosonic fields φ(x), θ(x), whose dynamics is encoded in the Luttinger liquid Hamiltonian in Eq. (A.10). In particular, when expressing the chiral fermion operators ψ R (x), ψ L (x) in terms of φ(x) and θ(x), one sets with the KF Γ such that {Γ, γ L } = {Γ, γ R } = 0. Together with the MM's, the KF forms a set of three real fermionic zero-mode operators. This is the case in which, as we discuss above, to consistently count for the corresponding degrees of freedom, we introduce a second wire, described again by a Luttinger liquid Hamiltonian such as the one in Eq. (A.10), with bosonic fieldsφ(x),θ(x), with parametersū,K and, more importantly, requiring the introduction of a second KF,Γ. The second wire is fully decoupled from the rest of the system. Therefore, for any values of the boundary couplings,θ(x) is pinned at both x = 0 and x = ℓ and, accordingly, . Therefore,φ(x),θ(x) are decomposed in normal modes according to Eqs. (A.11), with velocity and Luttinger parametersū,K. The calculation of the factor that the auxiliary QW contribute the total partition function,Z, is, therefore, a straightforward exercise in elementary algebra: the result isZ 6 withq = e − βπū ℓ (see Appendix C for the definition of the Dedekind function η(q)). At the disconnected fixed point of the phase diagram the QW is fully decoupled from the TS's. Therefore, φ(x), θ(x) take the expansion in normal modes in Eqs. (A.11), as well, and they will accordingly contribute the total partition function by a factor analogous toZ in Eq. (2). In addition, we have to account for the (4) real fermionic zero mode operators. As at the disconnected fixed point they are fully decoupled from each other, as well as from the dynamical degrees of freedom of the bulk, according to the above argument, we expect them to contribute to the partition function a factor equal to 2 elevated to the total number of real fermions divided by 2 (that is, 4). As a result, the total partition function at the disconnected fixed point is given bŷ Using the standard asymptotic expansions of the Dedekind function and of the elliptic functions in the ℓ → ∞ (q → 1)-limit, one eventually finds that which implies, for the g-function at the disconnected fixed point Turning on the (two)-boundary coupling to the MM, taking into account the type N CIBC's at the disconnected fixed point, the two-boundary Hamiltonian H b takes the form When K > 1/2, H b is a relevant operator, which drives the system toward a fixed point in which φ(x) is pinned both at x = 0 and at x = ℓ, as φ(0) = √ πν 0 , φ(ℓ) = √ πν ℓ , with integer ν 0 , ν ℓ . Accordingly, φ(x) and θ(x) take the mode expansion in Eqs. (A.12) and, in addition, taking into account the boundary conditions, one gets H b = −2it {(−1) ν 0 γ L + (−1) ν ℓ γ R } Γ. A mode expansion such as the one in Eqs. (A.12) implies a corresponding contribution to the total partition function such as the one in Eq. (2), except for switching K with K −1 . In addition, minimizing H b locks together the real fermions Γ andγ = 1 by the condition that physical states are annihilated by the complex fermionic operator a = 1 2 (γ + iΓ). This condition leaves unpaired the real fermion operator η = − 1 which, together withΓ, determines an additional degeneracy factor of 2. Therefore, the partition function and the g-function at the corresponding fixed point are given by The derivation of Eq. (7) is a first example of implementation of the Majorana-Klein hybridization [50] in computing the g-function at the A fixed point. Indeed, the real-fermion zero-mode 7 operator we combined together withΓ to obtain the degeneracy factor of 2 originates from a linear combination of γ L and γ R , with coefficients determined by the boundary conditions. In addition, we note that, introducing the auxiliary wire on one hand allowed us to perform the calculations in an unambiguous way, on the other hand gave us back g Disc and g A up to an over-all arbitrary, multiplicative factor. To get rid of the factor, we normalize g A to g Disc by considering the ratio which, on one hand shows that the ratio is independent of the (arbitrary) parameterK, as it must be, on the other hand that, consistently with the g-theorem, the renormalization group trajectories are expected to flow from the disconnected to the A fixed point when K > 1/2. Later on in the paper, we generalize the results of this section to a junction between a generic number N of QW's and a topological superconductor. In this case, the relevant Luttinger param- πu , with u, K being the Luttinger parameter of each QW and U being the inter-wire interaction strength. (See Appendix A for a detailed derivation and discussion of the Luttinger parameters for the N-wire junction.) 3. Review of the phase diagram of the 2-wire junction with a topological superconductor for 1 2 < K σ < 1 and of the Y-junction of interacting quantum wires To complement the calculation of the g-function and to pave the way to the derivation of the correspondence between the N=2-wire junction with a TS and the Y3J, here we review known results about the phase diagram of the two systems. In doing so, we also compute the g-function at the various fixed points of the (generically asymmetric) Y3J, and demonstrate the consistency between our results and what is known about the phase diagram of that system [42]. We now begin by briefly reviewing the phase diagram of the N = 2-wire junction with a TS, which we extensively analyzed and discussed in Ref. [18] in the regime K σ < 1.

3.1.
Phase diagram of the 2-wire junction with a topological superconductor for K σ < 1 We refer to section 5 for an extensive derivation, from the microscopic Hamiltonian, of the Luttinger liquid description of a generic junction between N-QW's and a TS and for the analysis of the corresponding phase diagram. Here, we just review the main features of the phase diagram in the case N = 2 starting from the effective Luttinger liquid Hamiltonian, H 2,B = H 2,Bulk + H (1) b,B,2 , with the bulk Hamiltonian with the plasmon velocities and the Luttinger parameters defined in terms of the microscopic system parameters as per Eqs. (A.9) for N = 2, and φ ρ (x) = 1 being the fields in the two wires, and analogous formulas for θ ρ (x), θ σ (x) in terms of θ 1 (x), θ 2 (x). By construction, H 2,Bulk in Eq. (9) corresponds to the symmetric version 8 of the model Hamiltonian Ref. [18] taken for K 1 = K 2 = K, u 1 = u 2 = u, r = 1 and α = π/4. Accordingly, the chiral fermionic fields in the QW's are bosonized as (see Appendix A for details) with a = 1, 2 and Γ 1 , Γ 2 fermionic KF's. Having stated this, the phase diagram for K σ < 1 is readily recovered from the analysis of Ref. [18]. The simplest fixed point corresponds to setting all the boundary couplings t a = 0. This is the disconnected fixed point, which, in a sample defined over a segment of length ℓ, corresponds to pinning θ 1 (x), θ 2 (x) at both boundaries x = 0, ℓ. The corresponding value of the g-function can be computed as done in Ref. [18], yielding, as a special case of the general formula we derive in section 5, the result with the factor 2 stemming from the fourfold degeneracy due to the two MM's in the two-boundary version of the model and to the KF's associated to the two QW's. (Note that, differently from what was done in Ref. [18], here we do not count the additional degeneracy associated with the MM's at the outer boundaries of the superconducting lead, as they contribute to the g-function just an overall factor of 2, which drops from the physically meaningful ratios between g computed at two different fixed points). At the disconnected fixed point, the boundary Hamiltonian describing the coupling of the QW's to the MM is presented in bosonic coordinates as The scaling dimension of H (1) b,B,2 , d b , can be readily derived using the transformation in Eqs. (A.6). The result is d b = 1 4K ρ + 1 4K σ . For 1/2 < K < 1 and for U > 0, we find 1 2 < d b < 1, which implies that H (1) b,B,2 always corresponds to a relevant boundary interaction. In addition, other boundary interactions, though not present in the original ("bare") Hamiltonian, can be generated from the RG. The first one corresponds to boundary normal intra-wire backscattering, described by operators of the form ψ † R,a (0)ψ L,a (0)+h.c., which are not effective, due to Dirichlet boundary conditions on θ a (x) at x = 0. Then, one has inter-wire normal backscattering, corresponding to operators of the form V Normal, (1,2) σ . Finally, one has inter-channel pairing, corresponding to operators of the form V Pair, (1,2) , with scaling dimension d Pair,(1,2) = K −1 ρ . Thus, as long as 1/2 < K ρ , K σ < 1, no relevant operators are allowed at the disconnected fixed point but the boundary coupling to the MM, H (1) b,B,2 . The relevance of the operators in H (1) b,B,2 implies that, as soon as (at least one of) the t a 's are 0, the corresponding operator(s) trigger an RG flow away from the disconnected fixed point. In Ref. [18] it is shown that, for K σ < 1, the junction either flow towards a fixed point with type A(N) boundary conditions in channel-1(2) (the A ⊗ N fixed point), or towards the complementary, 9 N ⊗ A fixed point. In both cases, a straightforward implementation of our method yields, for the g function The result in Eq. (13) corresponds to the symmetric limit of the junction discussed in Ref. [18]. Note that it also fixes an error in that reference, though without affecting the final result. The important point is the over-all factor √ 2, which ensures that ρ 1,1 < 1 as long as 1/2 < K. It actually comes from the correct counting of the degrees of freedom associated with zero-mode real fermion operators. In particular, for the RG flow towards the A ⊗ N (N ⊗ A) fixed point to occur, the bare couplings must be such that t 1 > t 2 (t 1 < t 2 ). In this case, the RG makes the running coupling corresponding to the larger bare coupling constant flow all the way to ∞. Let us assume this is t 1 . Accordingly, to recover the A ⊗ N fixed point, one considers the two-boundary version of H (1) b,B,2 , in which t 1 → ∞. In fact, this implies "locking" two of the four real fermions (2 KF's plus 2 MM's) into a linear combination annihilating the physical states, leaving the other two decoupled from the boundary interaction, with a total degeneracy factor of 2 in the total partition function and, therefore, a factor √ 2 in g A⊗N [18]. The set of allowed boundary operators at the A ⊗ N fixed point includes the same operators we listed at the disconnected fixed point, though realized differently, and with different scaling dimensions, due to the change in the boundary conditions for the bosonic fields [18]. In particular, normal intra-channel 1 backscattering corresponds to an operator V Intra,1 ∝ cos[2 √ πθ 1 (0)], with scaling dimension d Intra,1 = d −1 b . At the same time, interchannel normal backscattering and inter-channel pairing are described by boundary operators that are linear combinations of V a, (1,2) and of V b, (1,2) , respectively given by V a, (1,2) (1,2) together with their Hermitean conjugates. The operators in Eqs. (14) have the same scaling dimension, d a, (1,2) = d b, (1,2) , which corresponds to the result in Eq. (D.26) of Ref. [18] in the symmetric limit r = 1, α = π/4. As a result, one obtains d a, (1,2) πu , which is always < 1 for 1/2 < K < 1 and U > 0. Thus, one obtains that 1+K ρ K σ K ρ +K σ > 1, as long as K σ < 1. Therefore, we conclude that both inter-channel normal backscattering and inter-channel pairing are described by irrelevant operators. Finally, an additional boundary operator arises from the residual boundary coupling of channel 2 to the MM. This corresponds to the term ∝ t 2 in Eq. (12). Despite the fact that it appears to correspond to a relevant operator, due to the hybridization between γ L and Γ 1 in the state that sets in at the A ⊗ N fixed point, it becomes effective only to order t 2 2 , corresponding to an operator V 2,Res ∝ cos[2 √ πφ 2 (0)], with scaling dimension d 2,Res = 4 K ρ +K σ > 1. This eventually proves that, for K σ < 1, the stable phase of the N = 2 junction either corresponds to the A ⊗ N, or to the N ⊗ A fixed point.
An additional possibility is provided by the A⊗A fixed point of Ref. [18], with type A boundary conditions in both channels. The g-function at the A ⊗ A fixed point can be readily derived either 10 from the analysis of Ref. [18] (up to an over-all √ 2, as discussed above), or from the general result of section 5, taken for N = 2 and for N a = 2, N n = 0. As a result, one obtains From Eqs. (15) one readily checks that, as long as K σ < 1, one obtains ρ A⊗A < 1, as well asρ A⊗A < 1, which implies that the A ⊗ A fixed point is unstable to both the N ⊗ N, as well as to the A ⊗ N fixed point. This is consistent with the DEBC results about the set of allowed boundary operators at the A ⊗ A fixed point. Indeed, implementing type A boundary conditions at x = 0 for both φ 1 (x) and φ 2 (x), we see that the boundary operators describing inter-channel normal backscattering, as well as inter-channel pairing, are in general expressed as linear combinations of the operators V a, (1,2) , V b, (1,2) , respectively given by together with their Hermitean conjugates. Their scaling dimensions are accordingly given by d a, (1,2) (1,2) = K ρ . Thus, we see that, for K σ < 1, they both correspond to relevant boundary interactions. Other boundary interaction terms are determined by the operatorsṼ Res,1 ,Ṽ Res,2 describing the residual coupling to the MM which, as discussed in detail in Appendix B, in this case can be effective to first-order in the boundary interaction strengths, different from what happens at the A⊗ N fixed point. In particular, on applying the bosonization procedure to the operators derived in Appendix B, one obtainsṼ The corresponding scaling dimensions are readily derived to be equal to each other and given by d Res,1 = d Res,2 = K ρ +K σ 4 . Given the definition of K ρ and K σ in section 2, we see that they are both relevant, as long as K σ < 1. Incidentally, we note that the other allowed boundary operators, corresponding to intra-channel boundary backscattering processes, are realized as . Accordingly, they have the same scaling dimension, d Intra,1 = d Intra,2 = K ρ + K σ and, therefore, they are both irrelevant, for 1/2 < K and U > 0.
The conclusion that, for K σ < 1, there are two equivalent stable fixed points in the phase diagram of the N = 2 junction (the A ⊗ N and the N ⊗ A fixed points discussed above) implies that there must be a phase transition between the two of them. In Ref. [18], the phase transition has been identified at a FCFP in the phase diagram of the junction, which is attractive along the line in parameter space corresponding, in the symmetric case, to t 1 = t 2 , and otherwise repulsive. To show this, an effective means is to resort to the perturbative RG approach within the ǫ-expansion method. Basically, one assumes that the junction parameters are such that d b = 1 − ǫ, with 0 < ǫ ≪ 1, and accordingly derives the RG equations to the first nonlinear order in the boundary couplings, 11 so as to recover nontrivial zeroes for the β-functions corresponding to the FCFP. For the details of the systematic derivation of the corresponding RG equations we refer to Ref. [18] in the specific case N = 2, as well as to Appendix D.1 for the generalization of the procedure to a generic N, while here we just quote the final result. Specifically, as extensively discussed in Appendix D.1, one introduces the dimensionless running couplingst a = t a τ ǫ 0 , with the cutoff τ 0 ∝ D −1 0 . D 0 being a high-energy (band) cut-off for the system. Therefore, letting the scale run from D 0 down to the scale parameter D < D 0 , one obtains that the corresponding RG trajectories of the running couplings are determined by the differential equations with l = ln(D 0 /D), D being the running energy scale, and the function F (ν) defined in Eq. (D.7). In general, for small initial values of thet i 's, Eq. (18) implies a growth of thet a along the RG trajectories. Along the symmetric line t 1 = t 2 in parameter space, this takes the system to the FCFP discussed in Ref. [18], which corresponds to the nontrivial zeroes of the right-hand sides Alternatively, if the initial condition lies off the symmetric line, the RG trajectories flow towards either the A ⊗ N, or the N ⊗ A, fixed point, according to whether, at D = D 0 , one has t 1 > t 2 , or t 1 < t 2 . While an exact description of the FCFP is still missing, within the ǫ-expansion method it is possible to estimate the corresponding value of the g function to leading order in the ǫ, obtaining [18] which implies g FCFP g N⊗N < 1, consistently with the RG flow from the disconnected fixed point to the FCFP, for ǫ > 0.

Phase diagram of the Y-junction of three spinless interacting normal wires
The Y3J has been introduced and extensively discussed in Ref. [2] in the fully Z 3 -symmetric case (bulk and boundary interaction). Later on, in Ref. [42], the effects of relaxing the bulk Z 3symmetry have been considered. Here, we consider the most general situation in which the Z 3 symmetry between the QW's can be broken by the boundary interaction, or by the bulk Hamiltonian [42], or both. Accordingly, we use as bulk Hamiltonian of the (asymmetric) Y3J, H Bulk , given by (20) with K, K 3 Luttinger parameters of the QW's and the velocity u set equal in all three channels to avoid unecessary formal complications.
In the absence of a boundary interaction, the g-function for the Y3J can be readily computed following the recipe of section 2 for N = 3. In particular, since N is odd, we add an auxiliary disconnected wire, with parametersū andK and Klein factorΓ, so to recover a total even number of KF's. Taking into account that, for a generic K 3 , the g-function at the disconnected fixed point, g Disc (K, K 3 ), is given by with the factor 2 due to the four real fermions that determine a total degeneracy of 4. Keeping all the θ a (0) pinned and turning on a (non Z 3 -symmetric) boundary interaction, one may readily present the corresponding boundary Hamiltonian H b by implementing the transformation matrix M N in Eq. (A.6) with N = 3 to resort to the center of mass-and to the relative-field basis, that is, by setting As a result, one obtains with t a+1,a , a = 1, 2, 3 being the boundary interaction strengths (assuming the convention a + 1 ≡ 1 for a = 3) and the {Γ a }'s being the three KF's required to bosonize the fermionic fields of the three wires, which shows that H b only depends on the relative fields ϕ 1 (x), ϕ 2 (x). Therefore, in constructing other boundary fixed points, we only act on the boundary conditions on ϕ 1 (x), ϕ 2 (x), which are type N at the disconnected fixed point but, in general, can change at other fixed points [2,42]. Incidentally, we note that the right-hand side of Eq. (23) is the leading boundary operator allowed at the disconnected fixed point. It is a linear combination of the operators V a,a+1 (plus their Hermitean conjugates), with scaling dimensions respectively given by So, a necessary condition for the disconnected fixed point not to be stable is that either one of the scaling dimensions above (or both) become < 1.
A first, alternative, fixed point is recovered by assuming type A boundary conditions for both ϕ 1 (x) and ϕ 2 (x). This corresponds to the D P fixed point of Refs. [2,42]. To stabilize D P , we introduce a two-boundary pairing potential V P , given by and eventually send ∆ → ∞. Sending ∆ → ∞, one pins ϕ 1 (x) and ϕ 2 (x) at both boundaries. Taking this into account, one determines the corresponding spectrum of the zero-mode operators and, repeating the calculation of the g-function at the D P fixed point, one obtains which yields the ratio An important comment about Eq. (26) is that, despite the fact that, for K = K 3 , neither g Disc (K, K 3 ), nor g D P (K, K 3 ), are equal to the values derived in Ref. [2], the ratio between the two of them is the same as one would get by using the results obtained there. This is due to the fact that, in our derivation, we count the degrees of freedom associated with the real KF's, including the auxiliary one and, in addition, do not restrict our derivation to the sector involving the relative fields only. Nevertheless, the ratio between the two of them is consistent with Ref. [2]. Clearly, this further enforces out intuition that, despite the arbitrary aspects of our procedure, the ratios between values of the g-function at different fixed points do always give back the right, physical result.

The leading dimension boundary interaction at the disconnected fixed point is H b in Eq. (23). It is a linear combination of operators with scaling dimensions
At variance, at the D P fixed point, the leading dimension boundary operators are given by [2,42] together with their Hermitean conjugates. To derive the corresponding scaling dimensions, one has to resort to the center of mass-and relative field basis by using Eqs. (22) and to take into account that Θ(0) is always pinned, as Φ(0) never appears in the boundary interaction [2]. As a result, one obtains [42] . From the scaling dimensions of the boundary operators we see that there is a K 3 /K-dependent window of values of K in which neither the disconnected, nor the D P fixed point, is stable. Specifically, this happens for 1 < K < 3 for K = K 3 [2] and, more generically, for for a generic K 3 [42]. The absence of time-reversal breaking in H b in Eq. (23) rules out the possibility of stabilizing the "chiral" χ ± fixed points: thus, one concludes that, for K min < K < K max , the stable phase of the system either corresponds to one of the asymmetric A a -points emerging in the Y3J when the Z 3 -symmetry between the channels is broken, or to a generically asymmetric version of the M-FCFP found in the Z 3 -symmetric Y3J in the time-reversal symmetric case [2].
The symmetries of the bulk Hamiltonian in Eq. (20) naturally lead to two different types of asymmetric fixed points: the A 3 fixed point corresponds to QW-3 disconnected from the junction, while the two-wire junction between QW's -1 and -2 is "healed" (which is a natural consequence of having K > 1, once QW-3 is disconnected from the junction [15]), and the (equivalent, up to swapping QW-1 and QW-2 with each other) A 1 and A 2 fixed points, in which respectively QW-1 and QW-2 are disconnected from the junction. Mathematically, disconnecting QW-a corresponds to imposing type N (type A) boundary conditions on φ a (x) (θ a (x)), as well as type ). Accordingly, the calculation of the corresponding value of the g-function can be readily carried out, providing the result yielding the ratios It is useful to also compute the ratios with g D P (K, K 3 ). The result is An effective mean to infer the stability of the A a fixed points against the disconnected and the D P fixed point consists in using Eqs. (29,30) in combination of the DEBC analysis of the corresponding allowed boundary operators. To construct the leading boundary perturbation at the A a fixed point, one considers the operators T a,a+1 ,T a,a+1 , respectively given by [42] together with their Hermitean conjugates. T a,a+1 andT a,a+1 respectively correspond to the boundary operators bilinear in the {ψ R,a , ψ L,a }'s given by ψ † R,a (0)ψ L,a+1 (0) and ψ † L,a (0)ψ R,a+1 (0). Once the appropriate CIBC's are implemented, they only depend on the linear combinations of the φ a (0)'s and of the θ a (0)'s that are not pinned at the corresponding fixed point. In particular, the CIBC's corresponding to the A a fixed point are recovered by pinning the arguments of both T a+1,a+2 and T a+1,a+2 . The T a,a+1 ,T a,a+1 -operators are the only operators that may become relevant at the A a fixed point, with scaling dimension [42]. As a result, from Eqs. (29,30) one finds that, in order for the A 3 fixed point to be stable with respect to both the disconnected and the D P fixed point, the condition 1 < K < 1 + 2K K 3 has to be satisfied. In 15 addition, there must be no relevant boundary operators allowed at A 3 in order for it to correspond to the actual stable fixed point of the Y3J. This leads to the additional condition d A 3 is stable against both A 1 and A 2 , as well. This is ultimately consistent with the results plotted in Fig.5 of Ref. [42], as well as with the observation that, a small enough K 3 /K eventually makes the interaction in wire-3 to be effectively repulsive, thus triggering the disconnection of this wire from the junction, in agreement with the known results about junctions of Luttinger liquids [34,15,2]. Conversely, in order for either A 1 , or A 2 , to be stable against the disconnected, as well as the D P fixed point, the condition K+K 3 2K 3 < K < 2(K+K 3 ) 2K+K 3 has to be met. In addition, the condition In particular, the above conditions yield g A 3 (K, Finally, we note that there are regions in parameter space in which, though one has that one of the conditions ρ A 3 (K, In this case, based on the well-grounded results of Ref. [2] about the Z 3 -symmetric Y3J, we expect that the stable phase of the system corresponds to a (possibly non-Z 3 -symmetric) FCFP, which generalize the M-FCFP of Ref. [2]. In the Z 3 -symmetric case K = K 3 , the emergence of the M-FCFP can be inferred from the perturbative RG equations in Eqs. (D.22) of Appendix D.2, given by with 0 < ǫ(= 1 − K −1 ) ≪ 1, and the parameters b and c estimated in Appendix D.2 to be b ≈ 26.32, c ≈ 16.45. An important point about the β-functions in Eqs. (32) is that they are over-all ∝ ǫ. As a result, the M-FCFP is found to lie att 2,1 =t 3,2 =t 1,3 = t * = 1/ √ b + 2c, independent of ǫ. On one hand, this result points in the right direction. Indeed, analytical [2], as well as numerical [53], results for the conductance tensor at the M-FCFP ultimately show that it has to be finite, as ǫ → 0. Had we found an M-FCFP at t * going to zero as ǫ → 0, we would unavoidably get a conductance tensor going to zero as ǫ → 0, as well, which would be incorrect [2,53]. On the other hand, since there is no "small parameter", such as ǫ, that can be used to control the coupling strengths at the FCFP's, one cannot really expect Eqs. (32) to be reliable to make quantitative predictions on e.g. the conductance tensor at the FCFP, or on the g-function (at variance with what happens for the junction between N QW's and a TS). Yet, besides the emergence of the M-FCFP itself, other remarkable conclusions can be derived from Eqs. (32), such as that the RG-trajectories always point towards the Z 3 -symmetric M-FCFP, that is, any asymmetry in the boundary couplings is an irrelevant perturbation of the RG flow trajectories. In fact, this is a remarkable feature that, in the Z 3 -symmetric case, the Y3J shares with the topological Kondo effect, in which the magnetic 16 impurity is realized in terms of localized MM's [28]. In the general case K 3 K, we rather refer to the corresponding generalization of Eqs. (32) that we provide in Eqs. (D.23) of Appendix D.2.
In particular, looking for nontrivial zeroes of theβ a -functions of Eqs. (D.23), we see that the predicted values of the running couplings corresponding to the M-FCFP are either characterized by an "easy plane" asymmetry for K 3 > K (which impliest 2,1, * <t 3,2, * =t 1,3, * ), or by an "easy axis" asymmetry in the complementary case, K 3 < K (which impliest 2,1, * >t 3,2, * =t 1,3, * ). In both cases, the flow towards the M-FCFP always requires the relevance of all the V a,a+1 -operators entering H b , as we discuss in detail in section 4, when spelling out the correspondence between the N = 2 junction and the Y3J.
4. Correspondence between an N = 2 junction with a topological superconductor and a Y-junction of three spinless quantum wires In this section we discuss in detail the various aspects of the correspondence between a junction with 2 QW's and a topological superconductor and the (generically asymmetric) Y3J. For the purpose of this work, the correspondence is of the utmost importance for several reason. First of all, it works as a sort of "model duality", allowing for recovering results about the phase diagram of one of the two systems from the known (and controlled) features of the phase diagram of the other, in the various regions of the system parameters. Moreover, the correspondence is useful in computing the g-function of one model from known results on the other one. About this point, it is worth stressing that, as in our work we attribute physical meaning only to the ratio between the g-function at different fixed points of the phase diagram, contributions from modes not entering the correspondence factorize and cancel, when computing the ratios, which enforces the reliability of the correspondence to computing the IE. Finally, as the correspondence requires defining MM's in the N = 2 junction in terms of KF's in the Y3J, and vice versa, it provides also strong evidence for the fact that both real fermionic modes have to be taken into account, and considered on the same footing, when computing the g-function, which is one of the main points we make here.
For clarity, in the following we split the presentation of the correspondence in two sub-sections. In sub-section 4.1, we explicitly construct the mapping from the N = 2 junction to the Y3J junction. This allows us to use known results about the phase diagram of the N = 2 junction [18] to unveil specific features in the phase diagram of the Y3J, such as emergence of "planar" FCFP's (that is, with one of the boundary coupling strengths set to 0). In sub-section 4.2, we derive the mapping from the Y3J to the N = 2 junction. Reversing the direction of the correspondence allows us to employ the known results about the phase diagram of the non-Z 3 -symmetric Y3J to derive the phase diagram of the N = 2 junction for K σ > 1, a range of values of the system's parameters which was not discussed in Ref. [18].

From the N = 2 junction with a topological superconductor to the Y-junction
We now consider the asymmetric Y3J with K and K 3 set so that V 1,2 , defined in Eq.
. For the asymmetric Y3J, the boundary coupling flow is determined by the perturbative RG Eqs. (D.23). Since V 1,2 is irrelevant,t 2,1 expected to renormalize to 0 for D 0 /D → ∞ and, as a result, one may recover the phase diagram of the Y3J in this regime by restricting the analysis to the planet 2,1 = 0 in parameter space. Settingt 2,1 = 0 in the second-and in the third lines of Eqs. (D.23), we obtain the system of two coupled RG equations given by Remarkably, Eqs. (33) can now be consistently dealt with within the ǫ-expansion method, by set- Expanding to linear order in ǫ and neglecting subleading contributions (in ǫ) to nonlinear terms in theβ-functions, according to Eqs. (33) and to the definition of the function B in Eq. (D.21), which implies that terms ∝ B at the right-hand side of Eqs. (33) are all ∝ ǫ, we trade Eqs. (33) for the system where F is defined in Eq. (D.7). Apparently, Eqs. (34) correspond to the perturbative RG equations of an N = 2 junction with 1 The correspondence is clearly not accidental. Indeed, on performing the canonical transformations followed by the rotation and analogous rotation from

one obtains, for the bulk Hamiltonian
Oncet 2,1 is set to 0 in the boundary Hamiltonian (which corresponds to dropping the term ∝ V 2,1 in Eq.(23)), H b becomes Apparently,φ χ (x),θ χ (x) fully decouple from H b in Eq. (38). Moreover, shiftingφ ρ (x) by a constant, so that π Finally, performing the reverse canonical rescaling given by and setting Γ 3 → γ L , we recover the Hamiltonian for the N = 2 junction with a topological superconductor, with Luttinger parameters given by Besides the mapping procedure involving H b , to further ground the correspondence we now extend it to all the allowed boundary operators in the Y3J and in the N = 2 junction, at each fixed point in the boundary phase diagram of the two systems that we discuss in section 3. Starting with the disconnected fixed point, due to the condition K < 1, the only relevant allowed boundary operators are V 2,3 and and V 3,1 entering H b in Eq. (38). Consistently with their scaling dimensions, these are identified with the operators at the second and third line of the table in appendix A.a of Ref. [42]. Another operator, which is irrelevant due to our choice of the system parameter, is the boundary operator V 1,2 of the asymmetric Y3J, with scaling dimension 1/K. Referring to the table in appendix A.a of Ref. [42], it apparently corresponds to any of the operators listed at the first line, taken at the disconnected fixed point. According to the analysis of section 3.1, its counterpart in the N = 2 junction is the normal boundary backscattering operator V Normal, (1,2) , with scaling dimension 1/K σ = 1/K. Additional boundary operators can potentially appear in the Y3J, which are quartic in the fermionic fields of the Y3J such as, for instance (in the notation of Ref. [42]) T = T LR 31 T RL 32 . Quartic operators do not appear in the table in appendix A.a 19 of Ref. [42], which only contains quadratic operators: to make them relevant a strong, bulk interchannel attractive interaction is required, which we exclude here, as we only focus on repulsive, inter-channel bulk interactions. Yet, to complete the correspondence with the N = 2 junction, we see that T is the second boundary operator that must be identified with the boundary pairing operator of the N = 2 junction, V Pair, (1,2) , with scaling dimension 1/K ρ . Moving to the A ⊗ N and to the N ⊗ A fixed points of the N = 2 junction [18], based on the analysis of section 3, we naturally identify them with respectively the A 1 -and the A 2 -asymmetric fixed point of the Y3J. To further corroborate our identification, we now show that it is realized as a one-to-one correspondence between boundary operators in the two systems. Here, we only discuss the correspondence between the A ⊗ N and the A 1 fixed point. The complementary one can be readily recoverd by symmetry. At the A ⊗ N fixed point of the N = 2 junction, the first pair of allowed boundary operators corresponds to boundary inter-channel backscattering/pairing between channels 1 and 2. The corresponding operators are realized as a linear combination of V a, (1,2) , V b, (1,2) , respectively given by V a, (1,2) plus their Hermitean conjugates. (Note that, due to the boundary conditions at the A⊗N fixed point, V a, (1,2) and V b, (1,2) do no more correspond respectively to normal boundary scattering and to boundary pairing, as they instead do at the disconnected fixed point -see the discussion after Eq. (12) in sub-section 3.1. Instead, as we state above, normal boundary scattering and boundary pairing operators are realized as linear combinations of the two of them.) Inverting the transformations above to get back to the fields of the Y3J, it is not difficult to check that the operators in Eqs. (42) respectively correspond to the T RL 21 -and to the T RR 21 -operators at the A 1 fixed point of the asymmetric Y3J, plus their Hermitean conjugates (see appendix A.d of Ref. [42]). This is further confirmed by the observation that, using the results of section 3 for the N = 2 junction, the scaling dimension of the operators in Eqs. (42) are the same and are given by d a, (1,2) 3 )K , which is the correct result for the T RL 21 -and for the T RR 21 -operators at the A 1 fixed point of the Y3J [42].
A second class of boundary operators at the A ⊗ N fixed point corresponds to intra-channel 1 normal backscattering processes, that is, to the operator Considering thatφ χ (x),θ χ (x) are fully decoupled from H b , it is natural to assume that, throughout the whole phase diagram of the system,θ χ (x) is pinned at x = 0. As a result, going again backwards along the sequence of transformations discussed above, we express V Intra,1 in terms of the fields of the Y3J at x = 0 as The right hand side of Eq. (44) corresponds to the T RL 22 -operator at the second line of the table at appendix A.d of Ref. [42] (plus its Hermitean conjugate), as witnessed by the perfect agreement 20 between the scaling dimension of that operator and the result of section 3 for the dimension of V Intra,1 , d Intra,1 = 2KK 3 K+K 3 . Finally, in the N = 2 junction one has the residual coupling of channel-2 to the MM. When properly accounting for the "Schrödinger cat" nature of the state formed out of the hybridization between the MM and the KF Γ 1 , one obtains, as corresponding boundary operator, V 2,Res given by [18] . The A ⊗ A fixed point of the N = 2 junction corresponds to pinning both φ ρ (x) and φ σ (x) at x = 0. Accordingly, one naturally identifies it with the D P fixed point of the Y3J. To double-check the identification between the two fixed point, we note that the leading boundary operators at the A ⊗ A fixed point of the N = 2 junction corresponds to the V b, (1,2) operator at the second line of Eq. (16), as well asṼ Res,1 andṼ Res,2 in Eqs. (17). Following the correspondence between the parameters of the N = 2 junction and the ones of the Y3J, one finds that the corresponding scaling dimensions are given by .
A comparison of the results in Eqs. (46) with the [18] and of section 3.1, which were reached under the assumption that K σ < 1).
In view of the perfect correspondence of the N = 2 junction and the asymmetric Y3J with K < 1, 2K > 1 + K/K 3 , one naturally concludes that an analog of the FCFP found in the N = 2 junction along the symmetric linet 2 =t 3 exists in the phase diagram of the Y3J, as well. In particular, we infer that, for K < 1, 2K > 1 + K/K 3 the stable fixed point of the Y3J is either the A 1 or the A 2 asymmetric fixed point of Ref. [42], depending on whether, at the reference scale, t 3,2 >t 1,3 , ort 3,2 <t 1,3 , or the FCFP of the N = 2 junction located, according to the analysis of Ref. [18], att 3,2, * =t 1,3, * = ǫ/F 2 − K −1 .
Before concluding this sub-section, two remarks are in order. First of all, we would like to stress that the condition K < 1, which we have assumed at the start of the discussion here, has the mere effect of making V 1,2 irrelevant, thus allowing for dropping terms ∝t 1 from the following discussion. In fact, this condition can be relaxed and one can extend all the conclusions we derive here to the case K > 1, as well, but only providedt 1 is fine-tuned to 0 from the very beginning, and remains =0 along the RG trajectories. Secondly, we would like to emphasize that, in order to 21 make the mapping effective, we identified one of the three KF's, specifically Γ 3 , with the MM γ L emerging at the N = 2 junction with a topological superconductor. This points out, once more, that, to make a comprehensive discussion of the physics of real fermions at junctions of onedimensional interacting electronic systems and/or topological superconductors, KF's and MM's have to be considered altogether as actual degrees of freedom, despite the apparent conventional definition of the former ones as a mere mathematical means to properly follow the bosonization procedure.

From the Y-junction to the N = 2 junction with a topological superconductor
In the previous section we used the known results about the N = 2 junction [18] to infer the emergence of a planar FCFP in the phase diagram of the asymmetric Y3J for K < 1, 1 2K + 1 2K 3 < 1. Nevertheless, the identification of K with the Luttinger parameter K σ of the N = 2 junction makes it impossible to directly extend the correspondence to the case K > 1. In fact, the analysis of Ref. [18] is limited to the regime K σ < 1. For K σ > 1, two key things happen. First, the identification of the argument ν of the F -function in the perturbative RG Eqs. (34) with 2 − K −1 implies that ν > 1 for K > 1. For ν > 1, F (ν) < 0, with the corresponding disappearance of the FCFP along the diagonal in thet 3,2 −t 1,3 -plane. Second, DEBC method shows the emergence of a relevant operator at the A ⊗ N, as well as at the N ⊗ A, fixed point, given by the inter-channel normal backscattering operator, V Normal, (1,2) , which implies that, unless one fine-tunes to 0 the coupling strength in front of V Normal,(1,2) , neither A ⊗ N, nor N ⊗ A, are stable fixed points anymore. To figure out what the phase diagram of the N = 2 junction looks like for K σ > 1, in this section we reformulate the correspondence with the Y3J, but this time to retrieve informations about the N = 2 junction from what is known about the phase diagram of the asymmetric Y3J.
To begin with, we refer to the disconnected fixed point. There, as stated above, for K σ > 1, the leading boundary operators for the N = 2 junction are the couplings of the two wires to the MM, V b,1 (2) , and the inter-channel normal backscattering operator, V Normal, (1,2) , given by of scaling dimensions respectively given by At a given U, the condition K σ > 1 is recovered by setting 1 2 Thus, we conclude that, for 1 2 < K < K * (U), the disconnected fixed point is unstable, with three allowed independent relevant boundary operators.
Moving to the A ⊗ N fixed point (and/or to the complementary N ⊗ A fixed point), we see that, referring to the operators V a, (1,2) , V b, (1,2) , their scaling dimension can be rewritten as d a, (1,2) which is clearly < 1 for K σ > 1, K ρ < 1. Therefore, we conclude that both inter-channel normal boundary backscattering, as well as boundary pairing, provide relevant perturbations at the A ⊗ N fixed point, as soon as K σ > 1. Finally, we readily see that the A ⊗ A fixed point is not stable either, as a relevant boundary perturbation is provided by the operator V b, (1,2) of Eq. (16), with scaling dimension K ρ < 1. To get some insight on the phase diagram of the N = 2 junction for K σ > 1, we now employ the correspondence with the Y3J. In order to do so, we start by assuming, for the time being, that V b,1 , V b,2 and V Normal, (1,2) all have the same scaling dimension, that is, d b = d Normal, (1,2) . At N = 2, this condition requires LetK(U) be the value of K that satisfies Eq. (50) at given U. We obtainK (U)U πu = 4 5 and, in addition, due to the assumption of a repulsive intra-wire interaction, which implies 1 2 <K 2 (U) < 1, we also get 1 and K σ =K(U) √ 5, that is, Accordingly, H 2,Bulk now takes the form In deriving the boundary phase diagram one has to work on the semi-infinite system. Thus, the upper integration bound (ℓ) in the integrals at the right-hand side of Eq. (51) must be sent to ∞. Taking this into account, we now perform the canonical transformation followed by a rescaling by √ 3 of the x coordinate in the integrals involving the center-of-mass fields. As a result, Eq. (51) becomes while, upon also redefining the field Φ(x) according to √ 3Φ(x) → √ 3Φ(x)+ π 2 , the most general boundary interaction at the disconnected fixed point,Ĥ b , can be written as a linear combination of V b,1 , V b,2 and V Normal, (1,2) aŝ The right-hand side of Eq. (54) corresponds to the boundary Hamiltonian of a Y3J with relative fields ϕ(x), Φ(x), boundary couplings t 1 , t 2 , v 1,2 and Luttinger parameters K = K 3 = K σ . In section 23 3.2, we argue that, as soon as all three of the boundary couplings are 0, for 1 < K σ < 3 (as is the case here), the system flows towards the Z 3 -symmetric M-FCFP of the symmetric Y3J [2]. Thus, we conclude that, for K σ > 1, as soon as v 1,2 0, the stable fixed point of the N = 2 junction is realized outside of the t 1 − t 2 -plane. It sets in at a finite value of the three boundary couplings t 1 =t 2 = v 1,2 = t * , with t * corresponding to the M-FCFP of the Z 3 -symmetric Y3J with Luttinger parameter K = K σ . The M-FCFP is the endpoint of RG flow lines that, were v 1,2 = 0, would instead end up at the A ⊗ A fixed point of the N = 2 junction. To complete our derivation, we now discuss the phase diagram of the N = 2 junction when K ρ = λK σ /3, with λ 1, and K σ > 1. In this case, going backwards along the mapping we derived in section 4.1, we see that the N = 2 junction maps onto the (bulk) asymmetric Y3J we review in section 3.2, with K = K σ and K 3 = 2λ 3−λ K. Accordingly, one obtains for the scaling dimensions (1,2) depending on whether λ < 1 (λ > 1). Apparently, the difference between d b,1 (2) and d Normal, (1,2) could, in principle, trigger RG trajectories either towards an A a -asymmetric fixed point, or towards an asymmetric version of the M-FCFP. As a function of λ and of K σ , one obtains for the scaling dimensions of the relevant boundary operators at the A 3 and at the A 1 , A 2 fixed points of the corresponding Y3J the expressions Given the assumption K σ > 1, based on the discussion of sections 3.1,3.2, in the following, referring to the fixed points of the Y3J, we assume that the stable fixed point has to be identified with Eventually, from the correspondence rules of section 4.1, we make the appropriate identifications with the fixed points of the N = 2 junction. A straightforward algebraic derivation leads us to conclude that, depending on whether λ > 1, or λ < 1, there is the possibility of stabilizing either the A 3 , or the A 1 , A 2 fixed points. In particular, without considering additional constraints on the various parameters, one would obtain • λ > 1: In this case, for K σ < 3 4 + 9 4λ 2 the stable fixed point corresponds to the M-FCFP of the Y3J. For 3 4 + 9 4λ 2 < K σ < 6 − 3 λ the stable fixed point would correspond to either the A 1 or the A 2 fixed point and, for K σ > 6 − 3 λ , to the D P fixed point. However, these last two possibilities are ruled out by the observation that, by definition, one has K σ = 3 λ K ρ < 3 λ . Therefore, one obtains that 3 4 + 9 4λ 2 < K σ ⇒ 1 3 < λ < 1, against the initial assumption. As a result, for λ > 1 only the M-FCFP of the Y3J corresponds to a stable fixed point of the N = 2 junction.
• λ < 1: In this case, the M-FCFP of the Y3J corresponds to the stable phase of the system for K σ < 6 − 3 λ . For 6 − 3 λ < K σ < 3 4 + 9 4λ 2 the stable fixed point corresponds to the A 3 fixed 24 point. For K σ > 3 4 + 9 4λ 2 , D P would become the stable fixed point. Again, we rule out this last possibility, due to the observation that, by definition, K −2 σ + K −2 ρ = 2(K j ) −2 where, to avoid confusion, we here use K j to mean the Luttinger parameter of each QW in the N = 2 junction in the absence of bulk, inter-wire interaction. At this stage, we are assuming, just as in Ref. [18], 1 2 < K j < 1 (later on we discuss an extension of our analysis to K j > 1 in the absence of inter-wire interaction). As a result, we find the condition 9 Therefore, the only allowed phase transition happens for 6 − 3 λ = K σ , where the stable fixed point of the system changes from the M-FCFP of the Y3J (for K σ < 6 − 3 λ ) to the A 3 fixed point (for K σ > 6 − 3 λ ).
In terms of the parameters of the N = 2 junction, such a fixed point corresponding to settinḡ v 1,2 → ∞, andt 1 =t 2 = 0. Accordingly, we see that it corresponds to the perfect healing of the junction between wires-1 and -2, with the MM decoupled from the two wires. To double-check its stability, we resort to DEBC-approach, by imposing type N boundary conditions on Φ(x), ϑ(x) at x = 0 and, accordingly, type A boundary conditions on ϕ(x), Θ(x). In this case, the leading boundary perturbation is indeed realized as a linear combination of the V b,1 and V b,2 -operators in Eqs. (47), that is, by the hybridization between the normal wires and the MM, which now take the form both with scaling dimension d b = 3+λK 2 σ 4λK σ > 1. Therefore, we conclude that, as soon as K σ > 6 − 3 λ , the system is attracted towards the A 3 -like fixed point, in which the MM is "pushed out" of the quantum wires, which hybridize with each other to an effectively uniform wire, out of which lies the decoupled MM.
In conclusion, we have shown that the condition K σ > 1 is enough to reverse the phase diagram of the junction between two quantum wires and a topological superconductor, with respect to the result derived in Ref. [18] for K σ < 1. Specifically, at K σ > 1, the FCFP corresponds to the true stable phase of the system and is eventually identified with the (in general non-Z 3 -symmetric) M-FCFP of the Y3J. Further increasing K σ with respect to K ρ may eventually trigger an additional phase transition towards a phase corresponding to the perfect healing of the junction between wires-1 and -2, with the MM decoupled from the two wires. To evidence the new phases we find in the N = 2 junction by means of the correspondence with the Y3J, in Fig. 1 we draw the phase diagram of the N = 2 junction for K ρ < 1 by including, in addition to what we found in Ref. [18] for K σ < 1, the phases emerging when the parameter K σ is > 1.
Before concluding this sub-section, it is worth remarking how the correspondence between the N = 2 junction and the Y3J also allows for inferring the phase diagram of the former system in the case of zero inter-wire bulk interaction (U = 0), and attractive intra-wire interaction (K > 1). To do so, we first of all note that U = 0 implies K ρ = K σ = K. From Eqs. (41), one sees that this condition is recovered in the K 3 K → ∞ limit, which yields K ρ = K σ = K. Accordingly, 25 Figure 1: Sketch of the phase diagram of the N = 2 junction for K ρ < 1. The red curve corresponds to 1 4K ρ + 1 4K σ = 1: it separates the phase corresponding to the disconnected fixed point ( 1 4K ρ + 1 4K σ < 1) from the other, "nontrivial" phases. For K σ < 1, the phase of the junction corresponds to either the A ⊗ N, or to the N ⊗ A, fixed point, according to whether the initial boundary coupling t 1 is larger, or smaller, than t 2 . In the symmetric case t 1 = t 2 , the system's phase corresponds to the FCFP of Ref. [18]. When K σ becomes > 1, a phase that maps onto a (generically non-Z 3 symmetric) deformation of the M-FCFP of the Y3J opens till, after crossing the green line, corresponding to the curve K σ 1 + K −1 ρ = 6, an additional phase opens (DMM), corresponding to the MM fully decoupled from the junction with the other two wires. one sees that, at the disconnected fixed point, the leading boundary operator is again provided by a generic linear combination of V 2,3 and of V 3,1 in Eq. (23) (plus their Hermitean conjugates), all with scaling dimensions d V 2,3 = d V 3,1 = 1 2K . Thus, we recover the expected result that the disconnected fixed point is unstable for K > 1/2 [18]. At the A ⊗ N fixed point, the most general allowed boundary interaction contains V a, (1,2) and V b, (1,2) in Eqs. (42), both with scaling dimension d a, (1,2) (1,2) = 1+K 2 2K , of the intra-channel 1 normal backscattering operator in Eq. (44), with scaling dimension d Intra,1 = 2K, and the operator V 2,Res in Eq. (45) describing the residual coupling to channel 2, with scaling dimension d 2,Res = 2 K . Among all those operators, the only one that can become relevant for K > 1 is V 2,Res , whose scaling dimension becomes < 1 as K > 2. Finally, at the A ⊗ A fixed point, the leading boundary interaction is provided by a linear combination of V Res,1 andṼ Res,2 in Eqs. (17), with corresponding scaling dimension d Res = K 2 . Putting the above results all together, one therefore infers that, for U = 0, K > 1, the stable fixed point of the N = 2 junction either corresponds to A ⊗ N, or to N ⊗ A, for K < 2, depending on the initial values of the boundary coupling strengths. For K > 2, A ⊗ A becomes the stable fixed point of the junction. Moreover, consistently with the discussion of sub-section 3.1, as well as with the results of Ref. [18], one expects A ⊗ N and N ⊗ A to be separated by some intermediate phase(s). Whether this corresponds to just a FCFP, as it happens for K < 1 [18], or to more than one FCFP's, or even to a continuous line of fixed points, cannot be firmly stated with our method in this range of values of system's parameters and, very likely, to discriminate among the various possible options will require resorting to a nonperturbative, numerical approach to the problem.  [18] for the N = 2 junction, Refs. [2,7] for the Y3J).
To summarize the correspondence between phases of the N = 2 junction and of the Y3J, in table 2 we provide a synoptic view of corresponding fixed points in the two models, using a different color (red, instead than blue) to highlight phases that in either model are predicted by means of the correspondence with the other.
To complement the results of the previous sub-sections, we now briefly discuss how the correspondence between the N = 2 junction and the Y3J has to be implemented in computing the g-function at corresponding fixed points of the two models.

Calculation of the g-function at corresponding fixed points
The simplest fixed point in both the N = 2 junction and in the Y3J is the disconnected one, in which all the boundary interaction strengths are set to 0. As from Eqs. (11,21), at the disconnected fixed point, one obtains g Disc = 2[K ρ K σ ] 1 4 in the N = 2 junction, and g Disc = 2[K 3 K 2K ] 1 4 in the Y3J. On comparing the two results, the first observation is that, to recover the over-all factor of 2 in the case of the Y3J, one has to include in the calculation the auxiliary KFΓ, as well. Besides that, the two results are apparently not related to each other via the correspondence between the Luttinger parameters in the two models in Eqs. (41). This is due to the fact that the Y3J g-function receives contributions from overall degrees of freedom not entering the correspondence with the N = 2 junction, that is, the auxiliary field and the center of mass field φ χ (x). On recomputing g Disc in the Y3J by dropping those contributions, one eventually obtains the asymmetric version of the result of Ref. [2], that is that is, exactly the result one obtains when inserting Eqs. (41) into the formula for g Disc in the N = 2 junction. Having stated the correspondence between the g-function at the disconnected fixed points, in the following we consider the g-function at alternative fixed points always normalized to g Disc , in both models.
The A ⊗ A fixed point in the N = 2 junction corresponds to the D P fixed point of the Y3J. Here, despite the counting of the real fermionic degrees of freedom working differently in the two models, the results for the g-function are again consistent with each other. While in Appendix B we discuss in detail the derivation of the corresponding degeneracy factor in the N = 2 junction, it is worth recalling how one recovers it in the Y3J. Setting for simplicity t 2,1 = 0, when mirroring H b in Eq. (39), one obtains its two-boundary version, H (2) b , given by This is the two-boundary Hamiltonian for the N = 2 junction, except that now one has γ L = γ R = Γ 3 . When discussing the degeneracy factor due to the zero-mode real fermion operators, in Appendix B we separately consider this case, concluding that, when both φ ρ (0) and φ σ (0) are properly pinned, one recovers a total number of three real-fermion zero modes which, put together withΓ, provide the degeneracy factor of 4 (2) to the total partition function (to the g-function).
Once the correct degeneracy factors have been taken into account, in the N = 2 junction one obtains that ρ A⊗A = g A⊗A /g N⊗N = [K ρ K σ ] − 1 4 , which, using Eqs. (41), one readily shows to be the same as the result of Eq. (26) for the Y3J.
When computing the g-function at the A ⊗ N fixed point of the N = 2 junction, we have to pertinently modify the result in Eq. (13), due to the identity u σ = u ρ = u, which is a direct consequence of the mapping from the Y3J. In this case, implementing the approach of Ref. [18], ore readily derives, using Eqs. (29), the identity which shows that, once normalized to the g-function at the disconnected fixed point, g A⊗N g Disc in the N = 2 junction is equal to g A 1 g Disc in the Y3J, as expected from the correspondence between the two models (note the apparent difference between the right-hand side of Eq. (59) for g A⊗N and the result in Eq. (13). This is due to the condition u ρ = u σ = u which naturally arises from the mapping and, in this case, takes the place of the formulas one generally derives from Eqs. (A.9) of Appendix A.) Before concluding this section, a comment is in order about the correspondence between real fermion operators in the N = 2 junction and in the Y3J. On one hand, we see that it is rather straightforward in the single-boundary version of the models, as, in that case, one simply uses the observation that the center of mass field of the Y3J decouples from the junction dynamics and, therefore, the "left-over" KF can be formally mapped onto the MM in the N = 2 junction. On the other hand, the correspondence is not anymore straight when resorting to the two-boundary version of the model Hamiltonian to compute g. In this case, the different nature of the MM's, which are local in real space, and of the KF's, which are global along the full extent of a QW, results, for instance, in that, while in the N = 2 junction one has a TSS at each end of the system 28 and, therefore, two MM's, in the Y3J one still has only the KF associated to the center of mass field. Remarkably, as we discuss above, this mismatch can be fixed by adding the auxiliary KF Γ to the count of the degrees of freedom of the Y3J, though without the possibility of rigorously extending the Hamiltonian mapping to the two-boundary systems, as well. Yet, the very fact that the results are the same in the two models, apparently further supports the extension of the actual correspondence between the N = 2 junction and the Y3J also to the g-function at corresponding fixed points, which can be eventually regarded as a double check of the results about the mapping of the phase diagram of one system to the other.

Phase diagram and impurity entropy of a junction between N quantum wires and a topological superconductor
As a further example of application of our method for computing the g-function, in this section we discuss the fixed points in the phase diagram, and the corresponding calculation of the IE, in a junction between N QW's and a topological superconductor. For the sake of simplicity, in the following we make the symmetric assumption that the Luttinger parameters u, K are the same for each QW. In this respect, this is a symmetric multiwire generalization of the junction discussed, for K σ < 1, N = 2 and (partially) for N = 3, in Ref. [18]. (An error occurred there in the final numerical estimate of the g-function, which we amend here; it did not affect the final conclusions.) Referring to the TLL-model Hamiltonian for the junction in Eq. (A.5) of Appendix A, in the following, we assume that both intra-wire and inter-wire bulk interactions are repulsive, which implies K < 1, U > 0. In addition, we assume that 1/2 < K, which is a necessary condition to assure the relevance of the boundary coupling to the MM [17,18]. In Fig.2 we provide a sketch of the junction between N interacting QW's and a TS in the single-boundary version (which we use to discuss the phase diagram) and in the two-boundary version (which we use to compute the g-function). We now provide a discussion of the phase diagram, which basically generalizes the analysis of Ref. [18] to a generic N.

Phase diagram of the N-wire junction with a topological superconductor
The simplest fixed point in the phase diagram corresponds to having all the wires disconnected from the TS. This implies type N boundary conditions at x = 0 for all the channels. When turning on nonzero couplings to the TS, {t a }, taking into account the boundary conditions, one may write the bosonized boundary Hamiltonian at the disconnected fixed point, H (1) b,B,N , in the form The scaling dimension of H (1)  b,B,N , d b , can be readily derived using the transformation in Eqs. (A.6). The result is d b = 1 2NK ρ + N−1 2NK σ . For 1/2 < K < 1 and for U > 0, we find 1 2 < d b < 1 ∀N, which implies that H (1) b,B,N always corresponds to a relevant boundary interaction for the range of parameters considered. Thus, having nonzero t a 's paves the way to the opening of new phases, corresponding to additional fixed points in the phase diagram of the junction. To discuss them and especially their stability under RG flow, one has to first identify the corresponding CIBC's, and then to employ them to construct, within the spirit of DEBC approach, all the allowed boundary operators at a given fixed point and eventually to check whether any of them corresponds to a relevant perturbation. This can be readily done within the imaginary-time framework, in which we describe the wires using the Euclidean action S Eff only depending on the fields at x = 0. To derive S Eff , one has to integrate over the bulk fields everywhere in real space, except at x = 0. Doing so, due to the duality between the φand the θ-fields on the semi-infinite line, S Eff can be either expressed in terms of the fields φ a (τ) = φ a (x = 0, τ), or of the fields θ a (τ) = θ a (x = 0, τ) [2], as with φ(Ω) = dτ e iΩτ [φ 1 (τ), φ 2 (τ), . . . , φ N (τ)] T , θ(Ω) = dτ e iΩτ [θ 1 (τ), θ 2 (τ), . . . , θ N (τ)] T , the matrix M N defined in Eq.(A.7) of Appendix A, and the matrix K N given by Eqs. (61) are true in general. The specific choice of either one of the actions in Eqs. (61) depends on the boundary conditions on the various fields. At the disconnected fixed point, all the θ-fields obey Dirichlet boundary conditions at x = 0. Therefore, we must use S Eff [{φ a }] in Eq. (61), with all the φ's different from 0. A straightforward calculation allows us to derive the scaling dimension of H (1) b,B,N , which is simply given by In addition, while normal intra-wire backscattering plays no role, due to Dirichlet boundary conditions on θ a (0), inter-wire backscattering and inter-wire pairing between channels a and b respectively correspond to the operators V Normal,(a,b) and V Pair,(a,b) , given by 30 plus their Hermitean conjugates. Their scaling dimensions are respectively given by (assuming Given the assumption that K σ < 1 (which we relaxed, when discussing the correspondence with the Y3J in section 4), we conclude that both d Normal,(a,b) and d Pair,(a,b) are > 1 and, accordingly, ,N provides the only relevant perturbation at the disconnected fixed point. To move away from the disconnected fixed point, in analogy to our derivation of Ref. [18], we employ the perturbative RG approach within ǫ-expansion method, which we briefly review in in Appendix D.1. The corresponding RG equations for the running couplings are given by with l = − ln(D/D 0 ), D being the running energy scale, D 0 the high-energy (band) cutoff ∼ τ −1 0 , and the function F (ν) defined in Eq. (D.7), with ν = 1 NK ρ − 1 NK σ . In general, Eq. (65) implies a growth of thet a along the RG trajectories. This may either take the system to some FCFP, which generalizes the one discussed in Ref. [18] for N = 2, 3, or to pinning N a φ a (0)'s, leaving the corresponding θ a (0) unpinned. Because of the symmetry between the channels, in the following we assume always that such a fixed point corresponds to pinning the first N a φ a (x)'s at x = 0, leaving the remaining N n unpinned (so that, for instance, the disconnected fixed point corresponds to N a = 0, N n = N).
Let us consider the N a = 1, N n = N − 1 fixed point. From Eq. (60), we see that this is recovered by sendingt 1 → ∞ and minimizing the corresponding contribution to H (1) b,B,N . Accordingly, besides pinning φ 1 (0), this also requires "locking" the system into a state either annihilated by the Dirac fermion a 1 = 1 2 [γ L + iΓ 1 ], or by a † 1 , depending on the value at which φ 1 (0) is pinned. Taking this into account, we may list the various allowed boundary operators at that fixed point. First of all, the intra-channel normal backscattering operator in channel-1 is realized as To derive the corresponding scaling dimension, d Intra,1 , in Eq. (61) we set to 0 all the θ a 's but θ 1 . As a result, we eventually find plus their Hermitean conjugates. To compute the corresponding scaling dimensions, we have to account for the Dirichlet boundary conditions on φ 1 (0). To do so, we get rid of the corresponding field in the Euclidean action, by emplying the "reduced" action S Eff, (1) [{φ b }], given by NĨ N−1 , and I N being the N-dimensional identity matrix andĨ N being the N-dimensional square matrix with all the entries equal to 1 but the ones at the diagonal, which are equal to 0. Accordingly, we find d a, (1,b) By mathematical recursion, one may show that which eventually leads to the final result For 1/2 < K ρ , K σ < 1 and for U > 0 one obtains d Inter,(1,b) > 1, thus showing the irrelevance of the corresponding operators. Finally, we consider the residual coupling to the MM. In this case, as discussed at length in Refs. [17,18], though the residual coupling between, say, channel-2 and the Majorana mode seems to provide a relevant perturbation, in fact, it does not, due to the condition that the physical states must either be annihilated by a 1 , or by a † 1 defined above. This makes an operator such as 2it 2 γΓ 2 cos[ √ πφ 2 (0)] become effective only to second order in t 2 , where it effectively behaves like an operator V 2,Res ∝ cos[2 √ πφ 2 (0)], with scaling dimension > 1 for 1/2 < K < 1 and for U > 0. Accordingly, this is an irrelevant operator, which leads us to conclude that, as long as K σ < 1, the stable phase of the N-wire junction with a topological superconductor always corresponds to a N a = 1, N n = N − 1 fixed point.
While our above analysis can in principle be readily extended to any N a ≥ 2, in the following we limit ourselves to the case N a = 2 to show how, in this case, at least two relevant boundary operators emerge at the corresponding fixed point. Eventually, this leads to the conclusion that the corresponding fixed point is not stable, consistently with the result of Ref. [18] for N = 2. The instability of fixed points with an N a ≥ 3 can eventually be inferred by means of similar arguments. 32 Assuming N a = 2, the key operators correspond to normal boundary backscattering/pairing involving channels-1 and -2. Within DEBC approach, they are readily recovered as a linear combination of the operators V a, (1,2) and V b, (1,2) , given by plus their Hermitean conjugates. The corresponding scaling dimensions are accordingly given by Both d a, (1,2) and d b, (1,2) are < 1, implying that boundary operators encoding normal inter-channel backscattering and pairing both correspond to relevant boundary interactions. Accordingly, we conclude that the N a = 2 fixed point is unstable and, by means of an obvious extension of the argument, that any fixed point with N a ≥ 3 is unstable, as well. In conclusion, we see that, also for N > 2, the only stable fixed points in the phase diagram of the N-wire junction are the N-ones with N a = 1. As those are all equivalent to each other, there must be intermediate FCFP's separating the corresponding phases. FCFP's have been argued to potentially host "decoherence-frustrated" phases with reduced decoherence effects in the boundary quantum degrees of freedom [32,33,24]. In our case, FCFP's are expected to emerge at the bifurcations between RG trajectories leading to any one of the stable N a = 1 fixed points. While we are not able to provide an exact conformal boundary field theory description of the FCFP's, we can still access them in the ǫ-expansion framework. Indeed, they emerge as nontrivial zeroes of the β-functions at the right-hand side of Eqs. (65), with the corresponding boundary couplings satisfying the equations By inspection, we see that there is only a solution with all thet * a 0 (FCFP N ), corresponding tō .
Next (assuming N ≥ 3), it is possible to have nontrivial solutions in which onet * a = 0, all the others being 0 (FCFP N−1 ). These are given bȳ Going ahead (assuming N ≥ 4), we find N(N − 1)/2 FCFP's in which twot * a = 0, with all the others being 0, etc. Remarkably, Eqs. (74) do not exhibit solutions with just one t * a 0 and all 33 the others being = 0, which gives us one more insight about the possible topology of the boundary phase diagram of the junction. To do so, we first of all note that, if M (≤ N) bare couplings are equal to each other, then the symmetry among them is not broken along the RG flow generated by Eqs. (65). In this case, we therefore expect the junction to flow towards the FCFP N . At variance, a slight breaking of the symmetry between the couplings does, in fact, take the system out of the FCFP N . To show this, let us assume that, in the vicinities of the FCFP N , the couplings are set so thatt 1,L = t * (N) − ρ, whilet a,L = t * (N) + σ for a = 2, . . . , N, with 0 < ρ ≪ 1 and 0 < σ ≪ 1. On linearizing Eqs. (65), one obtains Once integrated, setting From Eqs. (78) one infers that the RG trajectories flow towards the FCFP N−1 fixed point, obtained by decreasingt 1 and symmetrically increasing all the other couplings. Conversely, if both ρ(0) and σ(0) are < 0, then a direct flow to an N a = 1 fixed point is recovered. Now, the above analysis can be straightforwardly iterated, eventually generalizing to the N-wire junction the RG flow diagram derived in Ref. [18]. In Fig.3, we draw a sketch of the minimal flow diagram for the junction. We see that, for 1/2 < K < 1, the RG trajectories flow away from the disconnected fixed point, either towards one of the N a = 1 fixed points, or towards some FCFP, depending on the symmetry between the initial values of the boundary couplings. Eventually, reducing the symmetry between the boundary couplings implies a flow between different FCFP's, till, when all the symmetries are removed, the system flows towards one of the maximally stable N a = 1 fixed points. We now compute the g-function at the various fixed points of the N-wire junction, eventually arguing that the corresponding results are consistent with the expected topology of the phase diagram only provided one properly accounts for the real fermionic modes, which is at the heart of our approach.

Impurity entropy at the fixed points of a junction between N quantum wires and a topological
superconductor To begin with, let us consider the disconnected fixed point, corresponding to type N CIBC's at both boundaries in each channel. As a general remark we note that, in the junction we consider 2 N 1 t t t Figure 3: Sketch of the typical renormalization group flow diagram of a junction of N interacting quantum wires with a topological superconductor drawn for U > 0 and K σ < 1 (see text). In this range of parameters, we see that the disconnected fixed point, as well as the fixed points with N a > 1, are unstable against turning on a nonzero boundary interaction. In the presence of a symmetry between two, or more than two, boundary couplings, the renormalization group trajectories flow towards FCFP's (drawn in green and blue). Eventually, when one boundary coupling takes over all the others, the system flows towards one of the N a = 1 fixed points.
here, we have 2MM's (one at each boundary), as well as N KF's, for a total of 2 + N real fermionic modes. Consistent with the discussion of section 2, if N is odd, we introduce the auxiliary wire, with Luttinger parametersū,K, providing an extra Klein factorΓ, so to make the total number of real fermionic zero mode operators even. From Eqs. (A.8,A.13,A.14) and taking into account the degeneracy factor associated with zero-mode real fermionic operators, for N even, we eventually find for the partition function at the disconnected fixed point the result with q c,r = e − uc,r βπ ℓ , and m N,φ ), and the matrix B N,e given by At variance, for N odd, one obtains  To extract the g-function, we have to consider the partition function in the ℓ → ∞ limit at fixed β, which can be readily done by employing the Poisson summation formula in the form presented in Eq. (C.10). As a result, we obtain from which we eventually obtain for g Disc the result To generalize Eqs. (85) to a fixed point with type A CIBC's in the first N a channels and type N in the remaining N n ones, we refer to Eq. (A.15) of Appendix A for the spectrum of the zero-mode operators. Accordingly, for N even, we obtain that the zero-mode contribution to the total partition function is given by , (N odd) .

36
In Appendix B we discuss in detail the calculation of δ e,o [N a , N n ]. Here, we just quote the final result for the g-function, which is On normalizing g N a ,N n to g Disc , we obtain the ratio which gives back the result of Eq. (13) for N a = N n = 1 and the results of Eq. (15) for N a = 2, N n = 0. Besides the consistency check, a first important result is that one obtains as well as which is < 1 for K σ < 1. Thus, the systematic calculation of the g-function at fixed points with given CIBC's ultimately confirms the phase diagram emerging from the perturbative RG approach combined with DEBC method. There are N equivalent stable fixed points, corresponding to N a = 1, N n = N − 1. These are separated by FCFP's that are expected to lie along specific symmetry line in the boundary parameter space [18]. To compute the g-function at the FCFP's, we employ the ǫ-expansion method discussed in Ref. [18]. Specifically, we assume that d b = 1 − ǫ, with 0 < ǫ ≪ 1 and eventually find that the FCFP's are all located at values of the boundary parameters t * ∝ (F (2 − K −1 σ )) − 1 2 , with the function F defined in Eq. (B.41) of Ref. [18] and reviewed here, in Appendix D. In general, letting M(≤ N) be the number of finite couplings t * ,M at a FCFP, we find t * ,M = ǫ (M−1)F (ν) . To proceed with the calculation of the corresponding value of the g-function, g FCFP,M , we go through exactly the same derivation of appendix G of Ref. [18]. As a result, to leading order in the t * ,M , we find A remarkable consequence of Eq. (93) is that, since, given two integers M, M ′ both ≤ N, we find , that is, if an RG trajectory takes place between two FCFP's, it must take the system towards the fixed point with the lower value of M, consistently with the result of Ref. [18] in the case N = 3, U = 0. 37

Conclusions
We discuss the method to consistently compute the g-function at the boundary fixed points of the phase diagram of junctions between interacting quantum wires and/or topological superconductors, involving real fermionic modes in the corresponding boundary Hamiltonian (localized Majorana modes and/or Klein factors). We show that, in doing the calculation, one has to treat all of the real fermionic degrees of freedom on the same footing, which is apparently a version of the Majorana-Klein hybridization phenomenon, which requires KF's to be considered as actual "physical" degrees of freedom in exactly the same way as MM's, when describing junctions between interacting quantum wires and topological superconductors [28]. Incidentally, in our procedure for computing g, we also introduced a means to avoid ambiguities in counting the degrees of freedom associated with an odd total number of real fermions, by introducing an auxiliary wire, fully disconnected from the junction. The additional wire has the effect of providing an additional KF, which makes the total number of real fermionic modes always even. While affecting the value of g at a specific fixed points, our procedure eventually gives back the right value of the ratio between g computed at two different points.
By comparing the results of perturbative RG approach and DEBC-method with the explicit calculation of the g-function and the implications of the g-theorem, we have mapped out a remarkable correspondence between the N = 2 junction and the non-Z 3 -symmetric Y3J, for suitably chosen values of the system parameters. In particular, we have employed the correspondence to recover informations about the phase diagram of the former system from known results about the phase diagram of the latter, and vice versa. Concerning the N = 2 junction, we have shown that the condition K σ > 1 is enough to reverse its phase diagram with respect to the result of Ref. [18] for K σ < 1. Increasing K σ to values > 1, we proved that the FCFP corresponds to the true stable phase of the system and is identified with the M-FCFP of the Y3J and, eventually, that, for large enough values of K σ , the N = 2 junction undergoes a phase transition to a phase with perfect healing of the junction between wires-1 and -2, with the MM decoupled from the two wires. Conversely, for the Y3J, we demonstrated the emergence of a "planar" FCFP's (that is, with one of the boundary coupling strengths set to 0), which is a novel feature, so far not discussed for such systems. In addition, we were able to infer the phase diagram of the N = 2-junction at zero inter-wire interaction and for K > 1 in each wire, a regime which was not discussed in Ref. [18].
Despite being effective in deriving a number of results on the phase diagram of both systems, the correspondence between the N = 2 junction and the non-Z 3 -symmetric Y3J still presents a number of "critical" issues, which will have to be further analyzed, possibly with the help of a numerical approach to the problem, such as the one employed in Ref. [53]. In particular, issues related to our work that deserve to be further analyzed are: -The failure of the ǫ-expansion method to provide quantitative results about the FCFP, when applied to the Y3J in the Z 3 -symmetric limit. As we show in Appendix D.2, the perturbative βfunction for the boundary running coupling in the Y3J contain terms that are all ∝ ǫ, which makes the perturbative RG approach not reliable for extracting informations about the FCFP, at variance to what happens in the N = 2-junction [18]. (Yet, it must be stressed that, while not applicable in general, the ǫ-expansion method works fine for the Y3J, as well, in some range of values of the system's parameters, such as the one considered in section 4.1, leading to Eqs. (33,34) .) -The nature of the FCFP that our correspondence predicts in the N = 2 junction with no interwire interaction, 1 < K < 2 in each wire, and for symmetric boundary couplings. In particular, it would be extremely interesting to figure out whether there is still just one FCFP and whether it is continuously connected to the one we found in Ref. [18] for 1/2 < K < 1, or if there is more than one FCFP's, possibly of some intrinsically different nature; -Whether the the correspondence extends to the g-function, despite the fact that the physical nature of the real modes in the Y3J and in the N = 2 junction are fundamentally different, is just an accident, or can apply, possibly in some different form, in similar systems.
Apart from the ones listed above, from our results, there are a number of issues that we left over and should be properly addressed, such as the relation between the validity of the g-theorem in the presence of real fermionic modes and the conservation of the total fermion parity, the explicit calculation of the g-functions at FCFP's where the ǫ-expansion method fails, or the extension of our derivation to systems such as the "Majorana-Kondo devices", which at the same time encompass Majorana and (topological) Kondo physics [29]. These topics are outside of the range of this work and we plan to address them in forthcoming publications.
long-wavelength fermionic modes allows for expanding the lattice fermion operators as c j,a ∼ {e ik f j ψ R,a (x) + e −ik f j ψ L,a (x)}, with the Fermi momentum ±k f = arccos − µ 2J . Bosonizing the chiral fermionic fields requires introducing N pairs of canonically conjugate fields {φ R,a (x), φ L,a (x)} (a = 1, . . . , N). In the noninteracting limit, they can be expressed in terms of N pairs of chiral bosonic fields as with the chiral bosonic fields satisfying the algebra with the (canonically conjugate) fields {φ a (x), θ a (x)} described by the bulk Hamiltonian H N with the Luttinger parameter K and the plasmon velocity u determined by the intra-wire interaction V and by the Fermi velocity in the wires, v f . By means of an appropriate orthogonal transformation, H N,B can be separated into independent terms by rotating to the basis of the center-of-mass fields Φ(x), Θ(x) and the relative fields ϕ 1 (x), . . . , ϕ N−1 (x) and ϑ 1 (x), . . . , ϑ N−1 (x), defined as with the matrix M N only depending on N and given by In terms of the rotated fields, one obtains Eqs. (A.9) yield (u ρ , K ρ ) = (u N , K N ) and (u σ , K σ ) = (u 0 , K 0 ), with K n = K/ 1 + (n−1)U K πu and u n = u 1 + (n−1)U K πu . Note that, in particular, Eqs. (A.9) are consistently defined only as long as KU πu < 1, which is our over-all assumption in this work. Also, for a repulsive inter-wire interaction (U > 0), by definition one always has K ρ < 1. As stated in section 3.1, H N,B , for N = 2, corresponds to the bulk Hamiltonian of Ref. [18] taken for K 1 = K 2 = K, u 1 = u 2 = u, r = 1 and α = π/4. Following the notation of Ref. [18], we use N to denote open boundary conditions at both boundaries in a single bosonic channel corresponding to a disconnected wire, which implies pure normal reflection at both boundaries. This implies open boundary conditions for the lattice fermions [20] or, in terms of the bosonic fields, pinning of θ a (x) at both boundaries, as θ a (0) = √ πn 0,a , θ a (ℓ) = √ πn ℓ,a , n 0,a , n ℓ,a ∈ Z. As a simple, paradigmatic, example we consider a single field φ(x) which, together with its dual field θ(x), is described by the TLL Hamiltonian Imposing Neumann boundary conditions on φ(x) at both boundaries implies the mode expansions with the oscillator modes satisfying the algebra [α n , α † n ′ ] = nδ n,n ′ and the spectrum of the zeromode operators given by p φ = √ πm φ , with m φ relative integer. A mode expansion complementary to the one in Eqs. (A.11) is recovered when imposing type A boundary conditions on φ(x) at both boundaries, and, accordingly, type N boundary conditions on θ(x). In this case, one obtains with the eigenvalues of p θ equal to √ πm θ , and m θ relative integer. For the N-wire junction, when all the wires satisfy type N boundary conditions at both boundaries, one obtains the mode expansion with the oscillator modes satisfying the algebra [α b,n , α b ′ ,n ′ ] = nδ b,b ′ δ n+n ′ ,0 , and the spectrum of the zero-mode operators given by with m 1,φ , . . . , m N,φ relative integers. In the general case in which one has type A CIBC's at both boundaries in the first N a channels and type N CIBC's in the remaining N n ones, Eqs. (A.14) generalize to for the eigenvalues of the zero-mode operators of the (Φ, ϕ 1 , . . . , ϕ N−1 ) and of the (Θ, ϑ 1 , . . . , ϑ N−1 )fields respectively, with m 1,θ , . . . , m N a ,θ and m N a +1,φ , . . . , m N,φ relative integers.
Moving to PSF, for a single-variable function f (x), it is defined starting from the quantity F, given by Assuming that, ∀ R ∈ Λ, ∃ (n 1 , . . . , n d ) such that R = d i=1 n i R i , we therefore obtain withf (p 1 , . . . , p d ) being the multidimensional Fourier transform of f (x 1 , . . . , x d ). On explicitly performing the integral, we eventually get with Λ * being the dual lattice of Λ. Now, in the specific problem we consider, we typically obtain Eq. (C.12) is the key equation we use to compute the g-function at the various fixed points of the systems we study in our paper.
Appendix D. Review of the ǫ-expansion approach to junction of quantum wires In this appendix, we review the ǫ-expansion approach to the RG equations and to the calculation of the g-function at the FCFP's in a junction between N interacting quantum wires and a topological superconductor and in the Y3J. In the former case, we revisit and generalize the derivation discussed in Ref. [18]. In the latter case, we highlight the peculiarites of the ǫ-expansion applied to the Y3J, estimate the g-function at the FCFP of such a system [2] and eventually discuss the analogies with the N-wire junction.
with a b, et cetera. Following the derivation of Ref. [18], we also make the assumption that the coupling to the Majorana mode is a slightly relevant operator, that is, that we have 1 − (2d b ) −1 = ǫ, with 0 < ǫ ≪ 1. Now, to leading order in ǫ, we set d −1 b = 2 in the integrals involving the OPE's at the right hand side of Eq. (D.3). Accordingly, due to the remarkable identity after subtracting the free energy correction times H b , as discussed above, we see that, to leading order in ǫ, no renormalization of the t a 's arises that is ∝ t 3 a . Instead, a nonzero renormalization of the boundary coupling arises from corrections to the boundary action that can be derived by a straightforward generalization of the analysis done in Ref. [18] for N = 2. In particular, the relevant correction turns out to be given by with that is, Eq. (65) of the main text.