Electronic correlations in magnetized helical edge states coupled to s-wave superconductors

TL;DR

This paper investigates how electron-electron interactions influence correlation functions and spin transport in magnetized helical edge states coupled to s-wave superconductors, revealing regimes where pairing promotes or suppresses various correlations.

Contribution

It provides a detailed theoretical analysis of the effects of interactions on correlation functions and spin transport in magnetized helical states with superconductivity, using bosonization and RG techniques.

Findings

Superconducting pairing promotes spin-density-wave correlations in the partially mixed regime.

Attractive interactions suppress singlet and triplet pairing correlations.

Zeeman field enhances pairing correlations under repulsive interactions.

Abstract

We theoretically study the role of electron-electron interactions in one-dimensional magnetized helical states coupled to an s-wave superconductor. We consider a partially mixed helical (superhelical) regime, where the magnetic field (superconductivity) has a dominant contribution. Using bosonization and renormalization group techniques, it is shown how the interactions affect the correlation functions in the system. In the partially mixed helical regime, we find that superconducting pairing promotes spin-density-wave correlations, while singlet and triplet pairing correlations suppress, especially under attractive interactions. In contrast, in the superhelical regime, a perturbative Zeeman field enhances both spin singlet and spin triplet-x pairings under repulsive interactions. We calculate both logarithmic and residual corrections to charge-density-wave, spin-density-wave, and pairing correlations, revealing short- and long-range behaviors. We further investigate spin transport properties supplemented by a renormalization group analysis of the temperature and frequency dependence of spin conductivity. We also analyze the transport of momentum-spin-locked carriers in the presence of Zeeman-induced and pairing-induced gaps, uncovering the effect of interactions and the interplay between the two gaps in helical systems.

Equations

$${\mathcal{H}}_{hel}$$

$${\mathcal{H}}_z$$

$${\mathcal{H}}_s$$

$${\psi }^{\prime }(x)=\sum _{r=R,L}{\chi }_r(r{k}_F){\psi }_r(x)e^{ir{k}_{F}x},$$

$${\chi }_r(k)=\frac{1}{\sqrt{2}}\left(\begin{array}{c}r{e}^{i{\vartheta }_k}\\ e^{-i{\vartheta }_k}\end{array}\right),$$

$${\vartheta }_k=\frac{1}{2}\arctan \left(-\frac{{\Delta }_z}{v_{F}k}\right).$$

$${ϵ}_{\pm }={\Delta }_s\pm \sqrt{v_F^2{k}_x^2+{\Delta }_z^2}.$$

$${\mathcal{H}}_{int}=$$

$$+$$

$${\psi }_r^{(†)}(x)=\frac{1}{\sqrt{2\pi a_0}}{\eta }_r{e}^{\pm ir\sqrt{4\pi }{\varphi }_r(x)},$$

$$\begin{array}{rl}\mathcal{H}=\frac{v}{2}[\frac{1}{K}({\partial }_x\Phi {)}^2+K({\partial }_x\Theta {)}^2]& +\frac{{\tilde{\Delta }}_z}{2\pi a_0}{e}^{i\sqrt{4\pi }\Phi }+h.c.\\ & +\frac{{\Delta }_s}{\pi a_0}\sin (\sqrt{4\pi }\Theta ),\end{array}$$

$${\tilde{\Delta }}_z={\Delta }_z{e}^{-i\alpha x},\text{with}\alpha =|k_F|K^2-2k_F,$$

$$K=\sqrt{\frac{1-\frac{y_{fb}}{2}+\frac{y_4}{2}}{1+\frac{y_{fb}}{2}+\frac{y_4}{2}}},$$

$$v=v_F\sqrt{{\left(1+\frac{y_4}{2}\right)}^2-{\left(\frac{y_{fb}}{2}\right)}^2}.$$

$$\frac{dK(l)}{dl}$$

$$\frac{d{\mathcal{Y}}_z(l)}{dl}$$

$$\frac{d{\mathcal{Y}}_s(l)}{dl}$$

$${\mathcal{O}}_{CDW}(x)$$

$${\mathcal{O}}_{SDW}^i(x)$$

$${\mathcal{O}}_{SS}(x)$$

$${\mathcal{O}}_{TS}^i(x)$$

$${\mathcal{O}}_{CDW}(x,\tau )$$

$${\mathcal{O}}_{SDW}^x(x,\tau )$$

$${\mathcal{O}}_{SDW}^y(x,\tau )$$

$${\mathcal{O}}_{SDW}^z(x,\tau )$$

$${\mathcal{O}}_{SS}(x,\tau )$$

$${\mathcal{O}}_{TS}^y(x,\tau )$$

$${\mathcal{R}}_{CDW}(r_1,r_2)$$

$${\mathcal{R}}_{SDW}^{ij}(r_1,r_2)$$

$${\mathcal{R}}_{SS}(r_1,r_2)$$

$${\mathcal{R}}_{TS}^{ij}(r_1,r_2)$$

$${\mathcal{R}}_{CDW}(r)$$

$${\mathcal{R}}_{SDW}^{xx}(r)$$

$${\mathcal{R}}_{SDW}^{yy}(r)$$

$${\mathcal{R}}_{SDW}^{zz}(r)$$

$${\mathcal{R}}_{SDW}^{yx}(r)$$

$${\mathcal{R}}_{SS}(r)$$

$${\mathcal{R}}_{TS}^x(r)$$

$${\mathcal{R}}_{TS}^y(r)$$

$${\mathcal{R}}_{TS}^z(r)$$

$$F_2(r)$$

$$F_3(r)$$

Full text

I Introduction

Topological quantum computation has been one of the most exciting topics in condensed matter physics in recent decades due to its immuneness to local perturbations Kitaev2003 ; Nayak2008 . Majorana zero modes are promising candidates for realizing the required platforms Kitaev2001 ; Alicea2012 , which are exotic quasiparticles that obey non-Abelian statistics and can be emerged in topological superconductors Beenakker2013 ; sato2016 . The nonlocal feature of Majorana modes makes them robust against local sources of perturbations, allowing operations that are topologically protected Microsoft2025 . The first quantum processing unit based on topological superconducting qubits has been demonstrated by Microsoft’s Station Q group Microsoft2025 ; Microsoft ; Garisto .

Meanwhile, conventional superconductors that have weak spin-orbit interactions, as well as intrinsic topological superconductors such as odd-parity pairing Ruhman2017 ; Maneo1994 ; Hor2010 ; Hsieh2012 , or superconductors without a center of symmetry that do not have a specific parity Sato2009 ; Tanaka2009 ; Bauer2004 ; Badica2005 , will not contribute to hosting Majorana modes. The idea of inducing topological superconductivity in hetero structures was formed by Fu and Kane Fu2008 . In this regard, an s-wave superconductor is placed near the topological surface states of a 3D topological insulator, during which the helical spin polarization Sato2009 of the surface states manifests a 2D effective p-wave superconductor Fu2008 ; Sato2010 ; Sau2010 ; Alicea2010 hosting Majorana states Lutchyn2010 . This issue has been extensively studied in many experimental works Williams2012 ; Wang2012 ; Yang2012 ; Wang2013 ; Cho2013 ; Oostinga2013 ; Finck2014 ; Snelder2014 . Moreover, recent experiments Mourik2012 ; Deng2016 ; Albrecht2016 ; Zhang2018 ; Nichele2017 ; Vaitiekenas2020 ; Aguado2017 have reported the evidence for Majorana modes in fabricated metal-superconductor or semiconductor-superconductor hybrid systems that host strong spin-orbit interactions Manchon2015 .

One-dimensional (1D) topological superconductivity can be engineered in systems hosting helical states that are proximitized by an s-wave superconductor in the presence Fu2008 ; Sato2009 ; Tanaka2009 ; Lutchyn2010 or absence Orth2015 ; Pedder2017 of a Zeeman field. The 1D helical states may arise as edge states in two-dimensional quantum spin Hall insulators Kane2005 ; Koing2007a ; Wu2006 ; Berneving2006 ; Qi2010 ; SpinHall . Furthermore, they can be found in quasi-1D systems with strong spin-orbit coupling, including semiconductor nanowires Lutchyn2010 , carbon nanotubes Klinovaja2011 ; Klinovaja2012a ; Izumida2017 ; Marganska2018 , bilayer graphene Klinovaja2012b , and graphene nanoribbons Klinovaja2013a . In these systems, the interplay of spin-orbit coupling, magnetic field, and proximity-induced superconductivity can establish topological superconductivity. Various setups have been proposed to probe topological superconductivity, including Josephson junctions Deacon2017 ; Pedder2017 ; Udupa2021 , superconductor–helical boundary lines Hart2014 ; Klinovaja2013b , and hybrid structures connecting multiple edge states Qingyun2017 ; Fleckenstein2021 ; Eriksson2015 .

In one dimension, electron-electron interactions give rise to the Luttinger liquid behavior 1Dreview ; Imambekov2012 , where correlation functions exhibit power-law scaling with exponents determined by the interaction strength Giamarchi2004 . The study of quantum correlation functions has been extensively investigated in many research areas LLRange ranging from spin chains Mahdavifar2015 to topological systems in the presence of superconductivity Ruhman2017 . The role of strong interactions in quasi-helical Majorana edge states Gangadharaiah2011 and in controlling the destructive effect of magnetic fields on superconductors has been pursued Klinovaja2014 ; Thakurathi2018 . Also, more exotic quasi-particles have been discovered Fendley2012 ; Orth2015 , in systems with strong electron interactions and a pairing mode such as Cooper-pair tunneling into a quantum Hall system Michelsen2020 , due to the proximity effect of a superconductor with helical quantum wires Michelsen2023 . The analysis of correlations is a powerful tool for understanding the behavior of dominant excitations in these systems, paving the way for future advances in the fields of spintronics Data1990 ; Fabian2007 , spin-orbitronics Fabian2007 ; Chappert2007 , and quantum computing Hutter2016 .

The relatively long coherence length and considerable relaxation time in superconducting systems make them ideal platforms for spin-based quantum manipulation. Additionally, the study of quantum spin phases Mahdavifar2024 opens new routes toward quantum information processing through the control of spin-spin correlations in partially mixed helical (PMH) states Hosseini2020 . The PMH regime, characterized by incomplete spin-momentum locking Soori2012 due to magnetizing helical states Wozny2018 , provides an interesting playground to explore the interplay of superconductivity, magnetism, and interactions. Motivated by these perspectives, it is interesting to investigate the behavior of spin and charge correlations in a 1D topological superconductivity formed at the interface between a conventional s-wave superconductor and a helical edge state subject to a Zeeman field Qi2011 .

In this work, we employ bosonization and renormalization group (RG) techniques to study the interaction effects in a 1D magnetized helical state coupled to an s-wave superconductor. We consider two complementary regimes: (i) a PMH state with a dominant Zeeman gap perturbed by proximity-induced superconductivity and (ii) a superhelical state with dominant superconducting pairing perturbed by a transverse Zeeman field. In both cases, we analyze many body quantum correlation functions such as the charge-density-wave (CDW), spin-density-wave (SDW), and superconducting pairing correlations. We also identify the leading instabilities in different interaction regimes for short- and long-range contributions.

Our results reveal that, under repulsive interactions, the perturbative Zeeman field can increase the superconducting correlations, particularly the triplet-x component. However, singlet pairing has negligible values. In contrast, attractive interactions stabilize singlet superconductivity while also amplifying triplet-x correlation as the leading phase. The competition between magnetic and superconducting gaps leads to non-trivial corrections to correlation functions, which we calculate both perturbatively and non-perturbatively. Furthermore, combining both the Kubo formalism and the Memory function method, we evaluate the frequency and temperature dependence of the spin conductivity and show how transport properties reflect the underlying correlated phases.

The paper is organized as follows. In Sec. II, we introduce the model system. Using a bosonization approach, gapped forms corresponding to both the Zeeman and superconducting proximity terms are derived. In Sec. III, we perform a RG analysis to examine the interplay between electron-electron interactions and the gap parameters. Section IV is devoted to the analysis of charge-spin density wave (CDW/SDW) and superconducting pairing correlation functions in various regimes of interaction and gap strengths. In Sec. V, we investigate corrections to these correlation functions in the strong, intermediate, and weak gap regimes, including logarithmic or perturbative approaches. In Sec. VI, the spin transport properties of the system is investigated, focusing on the frequency and temperature dependence of the spin conductivity using the Memory function approach and RG-improved results. Finally, we summarize our main findings in Sec. VII.

II Model

We consider a continuum model of a quantum helical state, described by Hamiltonian $\mathcal{H}_{hel}$ Koing2007a ; Wu2006 , subject to a transverse Zeeman field, with Hamiltonian $\mathcal{H}_{z}$ Hosseini2020 ; Soori2012 ; Wozny2018 . Additionally, electrons experience a proximity-induced pairing potential generated by a conventional s-wave superconductor, described by $\mathcal{H}_{s}$ Alicea2012 . So, the system, shown in Fig. (1), has the total Hamiltonian, $\mathcal{H}=\mathcal{H}_{hel}+\mathcal{H}_{z}+\mathcal{H}_{s}$, where

[TABLE]

Here, $v_{F}$ is the Fermi velocity, $k_{x}$ is the wave vector along the wire, and $\sigma_{x,y}$ are the Pauli matrices acting in spin space. $\Delta_{z}$ and $\Delta_{s}$ denote the Zeeman and superconducting pairing gaps, respectively.

The single-electron field operator $\psi^{\prime}(x)$ can be expanded in terms of right- ($r=R$) and left-moving ($r=L$) components as

[TABLE]

where $\psi_{r}(x)$ annihilates an electron moving in direction $r$ at position $x$ with Fermi momentum $k_{F}$. The spinor $\chi_{r}(k)$, containing the spin structure, is the eigenstate of the magnetized helical state, which is given by

[TABLE]

where the mixing angle $\vartheta_{k}$ encodes the Zeeman-induced spin rotation,

[TABLE]

In this basis, the energy spectrum of the system is

[TABLE]

The Zeeman gap partially disrupts the pure longitudinal spin channel of the helical states, causing slight spin mixing of the right- and left-moving quasi-particles. Figure 2 shows the expectation values of the x-component of helical electron spin, $S_{x}$, as functions of $\Delta_{z}$ and $k_{x}$ for the right movers (panel a) and the left movers (panel b). When the Zeeman gap $\Delta_{z}$ is zero, the spin expectation values of the right- and left-moving carriers are at the highest value in magnitude. This indicates that the spin quantization axis is parallel to the 1D conducting helical channel, i.e., the x-direction. For nonzero values of the Zeeman field, being in the y direction, the angle $\vartheta_{k}$ causes mixing between the different spin states near the helical edge, so that at $k_{x}=0$, the spin quantization axis of the helical wire rotates towards the direction of the Zeeman field. As $\Delta_{z}$ increases, the loss of spin alignment along the x-direction extends to larger momentum states resulting in the PMH state.

We incorporate the effects of electron-electron interactions through the interaction Hamiltonian Voit1995 ,

[TABLE]

where $g_{1}$, $g_{2}$, and $g_{4}$ denote the backscattering, dispersive, and forward scattering interaction strengths, respectively.

Using the standard bosonization method, we express fermionic operators in terms of bosonic fields. So, the right- and left-moving fermionic operators can be represented as Tsvelik

[TABLE]

where $a_{0}$ is a short-distance cutoff, $\phi_{r}$ are the chiral bosonic fields, and $\eta_{r}$ are Klein factors ensuring proper anticommutation relations. In this formalism, the full Hamiltonian $\mathcal{H}=\mathcal{H}_{hel}+\mathcal{H}_{z}+\mathcal{H}_{s}$ can be written as

[TABLE]

where $\Phi=\phi_{R}+\phi_{L}$, $\Theta=\phi_{R}-\phi_{L}$, and

[TABLE]

and the velocity and Luttinger parameter satisfy $vK=v_{F}$. The Luttinger liquid parameter $K$ and the velocity $v$ are given by

[TABLE]

[TABLE]

Here, we have defined $y_{fb,4}=\frac{g_{fb,4}}{4\pi v_{F}}$ and $y_{fb}=y_{2}-y_{1}$.

The Zeeman gap $\Delta_{z}$ and the superconducting gap $\Delta_{s}$ introduce sine-Gordon terms in the Hamiltonian (Eq. 10), which directly gap out the charge and spin sectors, respectively. Partial spin mixing affects the coupling between spin and charge degrees of freedom in collective excitations, resulting in competition controlled by both the Zeeman field and electron-electron interactions. For now, we focus on the simpler case $\vartheta_{k}=0$ to analyze the main form of the gap terms. Its effect will be included later. Furthermore, it should be noted that the proximity of helical charge carriers to a superconductor modifies the spin and consequently the charge sector, which is in sharp contrast to the usual cases. This effect arises from the helicity of the states, which influences the dominant phase of the system. As a result, the system may undergo a phase transition controlled by electron interactions, the strength of the Zeeman field, and superconducting pairing, as will be discussed in the following section.

III Renormalization Group analysis

In the presence of electron-electron interactions, the system parameters must be renormalized. This renormalization changes the relevance of the gap-opening terms from the Zeeman and superconducting couplings, thereby modifying the phase transition boundaries. The RG flow equations for the dimensionless helical Luttinger parameter $K(l)$ and the dimensionless coupling constants $\mathcal{Y}_{z}(l)$ and $\mathcal{Y}_{s}(l)$ are given by

[TABLE]

where $\mathcal{Y}_{z}$ and $\mathcal{Y}_{s}$ are the dimensionless Zeeman and superconducting gap parameters, respectively, defined as $\mathcal{Y}_{z,s}=\frac{\sqrt{2}\Delta_{z,s}a_{0}}{v}$.

The RG equations (14)-(16) describe how the Luttinger parameter $K$ and the gap amplitudes $\mathcal{Y}_{z}$ (Zeeman) and $\mathcal{Y}_{s}$ (superconductivity) evolve with the change of length scale $l$. From Eq. (14), one observes that the flow of $K$ is suppressed by the Zeeman gap, which enhances repulsive interactions (thus lowering $K$), whereas superconducting pairing tends to increase $K$, corresponding to effectively stronger attractive interactions. According to Eq. (15), the Zeeman gap $\mathcal{Y}_{z}$ increases under RG flow when $K<2$, indicating a relevant perturbation that opens a magnetic gap in the spectrum. In contrast, $\mathcal{Y}_{z}$ becomes irrelevant under RG flow if $K>2$. From Eq. (16), one can see that the superconducting gap $\mathcal{Y}_{s}$ is RG relevant when $K>\frac{1}{2}$, promoting superconductivity in regimes of strongly attractive interactions. Thus, the competition between these two gaps and their mutual feedback with the interaction parameter $K$ allow us to map the phase diagram of the system. Furthermore, we note that, unlike many previous studies Braunecker2009a ; Braunecker2009b ; Visuri2020 ; Citro2020 ; Giamarchi1991 , in our model both gap amplitudes remain finite even in the absence of interactions. This means that the gaps cannot act directly on each other. Instead, the presence of each gap alters $K$, and $K$ also controls the flow of both gaps.

Figure 3 illustrates the RG flow trajectories, determined by numerical evaluation of Eqs. (14)-(16), for the Zeeman and superconducting gaps as a function of the interaction parameter $K$. The Zeeman gap (orange curves) becomes relevant and flows to the strong gap regime at $K<2$, signaling a magnetically dominated phase. In contrast, the superconducting gap (blue curves) increases when $K>\frac{1}{2}$, corresponding to a superconducting phase stabilized by attractive interactions. As a result, below the separatrices, where the Zeeman (superconducting) gap is strong and relevant, the superconducting (Zeeman) gap is weak and irrelevant. Moreover, there are intermediate regimes above the separatrices of magnetism and superconductivity where their gaps can compete and coexist, particularly for $1/2<K<2$. Thus, when one of the gaps (Zeeman or superconducting) is dominant or even moderate, the other can be treated as a perturbation. Consequently, the mutual effects of interactions with both magnetism and superconductivity in helical states can give rise to the rich phase diagram.

IV Density and pairing correlation functions

In this section, we calculate the CDW and SDW correlation functions as well as the singlet (SS) and triplet (TS) superconducting pairing correlations for the helical states. Furthermore, we investigate how these correlations change in the presence of the combined pairing gap and Zeeman gap which enter the Hamiltonian through sine-Gordon terms arising from the superconductivity and Zeeman field to the helical state.

To proceed, we use the operators corresponding to the above-mentioned physical quantities Giamarchi1986 ; Giamarchi1988 ; Giamarchi1989 ,

[TABLE]

where $\psi^{\prime}(x)$ is the single-electron field operator defined in Eq. (4) and $\sigma_{i}$ are Pauli matrices acting in the spin space. We note that the spin modes of the right-moving and left-moving electrons correspond to spin-up and spin-down states, respectively. In the bosonization framework, the CDW, SDW, and superconducting pairing operators can be expressed in terms of the bosonic fields as,

[TABLE]

We use the definition of the density and pairing correlations as Giamarchi2004 :

[TABLE]

where $T_{\tau}$ is the imaginary-time ordering operator, $r_{1}=(x_{1},v\tau_{1})$, $r_{2}=(x_{2},v\tau_{2})$, and $i,j=x,y,z$. In the following, we set $r_{1}=r$ and $r_{2}=0$. Using the bosonized form of operators (21)–(26), the correlation functions with respect to $\mathcal{H}_{hel}$ can be expressed as

[TABLE]

where

[TABLE]

In the following section, we incorporate the effects of superconducting pairing and Zeeman perturbations and calculate their impact on the correlation functions. In the weak-gap regime, these corrections are treated with two complementary approaches based on whether the operators contain exponential or non-exponential terms. For the exponential (long-range) components, we apply logarithmic corrections derived from the RG equations Giamarchi1989 . For the non-exponential terms (short-range), we use a perturbative expansion method up to second order in the gap parameters Citro2020 .

V Corrections to correlation functions in the strong, intermediate, and weak regimes

As already discussed above, a weak gap in one channel corresponds to a strong gap in the other. We subsequently can evaluate corrections introduced by a weak pairing (Zeeman) gap to a helical wire with an intermediate Zeeman (pairing) gap in the following Sections (V.1 and (V.2)). This approach allows us to identify marginal boundaries for each gap in the weak regime while accounting for the presence of the relevant gap.

V.1 Zeeman-induced PMH egde perturbed by superconductivity

The presence of a superconductor near a helical state experiencing a Zeeman field can significantly modify the dominant fluctuations in the system. Here, we consider two cases: Strong Zeeman and intermediate Zeeman, and examine corrections introduced by superconducting perturbations.

V.1.1 Strong Zeeman with Weak Pairing

In the strong Zeeman regime, the bosonic field $\Phi$ becomes ordered, causing its dual field $\Theta$ to fluctuate rapidly. As a consequence, long-range correlations involving $\Theta$ vanish, suppressing pairing correlations. Following the calculations in Ref. Bakhshipour2024 , only the $x$- and $y$-components of the SDW survive, modified by the strong Zeeman gap as

[TABLE]

where

[TABLE]

The $x$-component represents short-range correlations similar to the helical case but with renormalized Luttinger parameters $K^{*}$ and velocity $v^{*}$ affected by the Zeeman gap. The averaging over the pinned field $\Phi$ yields a constant value for the $y$-component, indicating its stability.

V.1.2 Intermediate Zeeman with Weak Pairing: Zeeman-Affected Hamiltonian

We now consider the PMH state formed by an intermediate Zeeman field and study the effect of weak superconductivity. In this case, the sine-Gordon-type Zeeman term can be approximated by a quadratic form, assuming the Zeeman gap lies in a moderately relevant regime. Including the spin mixing angle $\vartheta_{k_{F}}$, the effective Hamiltonian reads

[TABLE]

where the magnetized Luttinger parameter $\tilde{K}$ and velocity $\tilde{v}$ are given by

[TABLE]

with $\tilde{v}_{F}=v_{F}\sqrt{1-\left(\frac{\Delta_{z}}{\epsilon_{F}}\right)^{2}}$ the Fermi velocity of the magnetized helical modes. The superconducting gap is effectively $\Delta=\Delta_{s}\sin(2\vartheta_{k_{F}})$.

The corresponding RG flow equations are

[TABLE]

where $\mathcal{Y}=\frac{4\Delta a_{0}}{\tilde{v}}$ is the dimensionless superconducting gap parameter. Equation (47) shows that $\mathcal{Y}$ becomes relevant for $\tilde{K}>\frac{1}{2}$, where $\tilde{K}$ is the Zeeman-renormalized Luttinger parameter. We further relate the magnetized Luttinger parameter $\tilde{K}$ to the interactions by introducing the magnetized-interaction parameter $\mathcal{Y}_{\vartheta_{k_{F}}}$:

[TABLE]

where $\mathcal{Y}_{\vartheta_{k_{F}}}=y_{bf}-1$ serves as an effective interaction parameter connecting the electron-electron interactions to $\tilde{K}$. Here, $y_{bf}$ is taken to equal $y_{4}$. Strong repulsive interactions drive the superconducting gap into the irrelevant regime. When $y_{bf}$ crosses the threshold of 1, and repulsive interactions weaken, the gap exits the marginal regime and becomes relevant. Without interactions, the system experiences a relevant gap. However, attractive interactions cause the superconducting gap to become more pronounced, pushing it deeper into the relevant regime as the attraction becomes stronger.

The behavior of the inverse Luttinger parameter $\tilde{K}^{-1}$ as a function of both the Zeeman gap $\Delta_{z}$ and the electron-electron interaction parameter $y_{bf}$ is illustrated in Fig. 4. In panel (a), which is for the repulsive interaction regime ($y_{bf}>0$), $\tilde{K}^{-1}$ increases monotonically with $y_{bf}$. It also shows a modest increase with an increase in the Zeeman gap $\Delta_{z}$. So, increasing the Zeeman gap shifts the boundary for the relevance of the superconducting gap to smaller values of the repulsive interaction. Panel (b), which is for the attractive regime ($y_{bf}<0$), shows that $\tilde{K}^{-1}$ is suppressed as the attraction (more negative $y_{bf}$) and $\Delta_{z}$ increase. So, in contrast to panel (a), an increase in the Zeeman gap enhances the strength of the superconducting gap, supporting the emergence of a robust pairing phase even for weaker attractive interactions.

V.1.3 Zeeman-Affected Density and Pairing Correlation Functions

In this subsection, we calculate the charge and spin correlation functions as well as the pairing correlations for the Zeeman-induced PMH states. We further examine the fate of these correlations in the presence of the combined Zeeman-pairing gap, which enters as sine-Gordon term into the Hamiltonian, arising from the superconducting coupling with the magnetized helical wire. In this scenario, the spins experience an additional gap along with the Zeeman gap.

Equations (17)–(20) yield the correlation operators for the PMH edge as

[TABLE]

Inserting Eqs. (50)–(56) into the correlation function definitions (27)–(30), the magnetized correlations with respect to $\mathcal{H}_{PMH}$ are obtained as

[TABLE]

where the functions $\tilde{F}_{2}(r)$ and $\tilde{F}_{3}(r)$ are defined as

[TABLE]

We note that the effect of the Zeeman field is encoded in the Luttinger parameters $\tilde{K}$ and $\tilde{v}$ as well as in the spin mixing angle $\vartheta_{k_{F}}$.

SDW states naturally emerge in low-dimensional electron systems due to Fermi surface nesting and strong electron-electron correlations. The x-phase, in particular, is characterized by a helical spin modulation resulting from spontaneous breaking of spin-rotational symmetry. Application of an external Zeeman field destabilizes this helical spin configuration by explicitly breaking helicity conservation. This symmetry breaking manifests as deformation or suppression of the spin texture, depending sensitively on the magnitude and orientation of the field.

Figure 5 shows the short-range behavior of the corrected correlation functions for CDW and SDW orders in the PMH regime versus the interaction strength $g$, distance $x$, and Zeeman gap $\Delta_{z}$. The first row, panels (a)-(c), corresponds to attractive interactions, $g<0$; the second row, panels (d)-(f), to repulsive interactions, $g>0$.

Panels (a) and (d) present the short-range CDW amplitudes for attractive and repulsive interactions, respectively. The CDW amplitude exhibits a pronounced peak at large negative $g$ and small separations, signifying robust local charge order when interactions are strongly attractive and the Zeeman gap is weak. The CDW correlation function decreases with increasing Zeeman field, and this suppression is further enhanced when the Zeeman term becomes more relevant due to electron-electron interactions. Under repulsive interactions, the maximum amplitude of the CDW shifts toward small positive values of $g$. In this regime as well, increasing the Zeeman field leads to a reduction in the CDW correlation. Stronger repulsion and larger Zeeman gap further reinforce this suppression. In both cases, increasing the separation suppresses the CDW rapidly, establishing that the enhancement of short-range CDW correlations is confined to a strong coupling and a minimal Zeeman field.

Panels (b) and (e) show the corresponding short-range amplitudes for the SDWxx in attractive and repulsive regimes, respectively. Throughout the parameter space, attractive and repulsive interactions result in SDWxx amplitudes that are smaller, by more than three orders of magnitude, than those of the CDW channel. Interestingly, unlike in the CDW case, a strong peak takes place at a small negative interaction corresponding to a deeper Zeeman gap. This peak diminishes as the system goes toward stronger attractive interactions with the weaker Zeeman gap. Moreover, SDWxx suppresses with increasing Zeeman strength in both interacting regimes. This non-monotonic behavior implies the subtle competition between external magnetic field and the electron-electron interaction. However, similar to the CDW case, in both attractive and repulsive regimes the amplitudes remain negligible and decrease rapidly with increasing distance.

Panels (c) and (f) display the behavior of the SDWyy correlation in attractive and repulsive regimes, respectively. The SDWyy channel exhibits negligible amplitudes seen in the longitudinal spin sector, regardless of whether interactions are attractive or repulsive. Here, a large enhancement is observed at both large attractive interaction and large Zeeman field. This is in contrast to usual cases where a strong Zeeman field tends to break spin rotational symmetry and suppress transverse spin order. However, in the present case, the renormalization of the Zeeman coupling due to interactions can reinforce this component of the SDW. Consequently, even at fixed Zeeman strength, stronger attractive interactions enhance effective magnetic fields felt by spins, potentially stabilizing or restoring long-range correlations weakened by external perturbations. Furthermore, in all cases, the amplitude is quickly diminished as the distance is increased. Thus, local transverse spin-density fluctuations remain well below the leading CDW response in both interaction regimes and cannot develop into dominant short-range order.

Figure 6 provides an analysis of the long-range behavior of the density wave and superconducting pairing correlations within the PMH state. The figure shows the correlation amplitudes as a function of spatial separation $x$, interaction strength $g$, and Zeeman gap $\Delta_{z}$. The results are organized into two main regimes: attractive interactions (panels (a)-(e)) and repulsive interactions (panels (f)-(j)).

Panels (a) and (f) show the behavior of the long-range CDW correlations for attractive and repulsive interactions, respectively. For repulsive interactions, the CDW amplitude is most prominent, reaching a maximum at strong coupling and highest $\Delta_{z}$. It decreases rapidly as $\Delta_{z}$ decreases. In the attractive case, the CDW remains vanishingly small over a broad range of parameters. But, its enhancement for strong Zeeman gaps at weak coupling is less pronounced.

Panels (b) and (g) show the SDWyy component. For attractive interactions (panel b), the SDWyy is highly suppressed, maintaining only negligible amplitudes at weak coupling. Such amplitudes are also less sensitive to $\Delta_{z}$, indicating that Zeeman fields cannot induce substantial spin ordering in this regime. In contrast, under repulsion (panel g), the SDWyy correlation exhibits a large amplitude at strong couplings. Moreover, increasing the Zeeman field strengthens the long-range SDWyy correlations and promotes a dominant ordering tendency, possibly leading to a leading spin phase.

Panels (c) and (h) present the SDWzz correlations. Under attractive interactions (panel c), the SDWzz remains almost suppressed, revealing substantial suppression by magnetic effect and attractive correlation. It generally remains subdominant throughout the parameter space. For repulsive interactions (panel h), interestingly, this component shows a similar profile compared with the transverse SDWyy, with a maximum appearing at larger $\Delta_{z}$ and strong repulsion. Although its amplitude remains smaller compared to the SDW in the $y$-direction.

Panels (d) and (i), show the singlet superconducting (SS) correlation amplitude. With attractive interactions (panel d), the SS channel is the strongest in the regime of strong attraction and a high Zeeman field. This implies the stability of conventional superconducting pairing. The amplitude reduces sharply as $\Delta_{z}$ decrease. Under repulsive interactions (panel i), the singlet correlation is completely suppressed, remaining essentially near zero across the explored regime. But still the Zeeman field promotes this correlation.

Panels (e) and (j) illustrate the triplet-x superconducting component. Notably, in the attractive regime (panel e), the triplet correlation is strongly enhanced by increasing the Zeeman field, yielding a distinct ridge in parameter space. In this regime, triplet-x component is a dominant phase. This is consistent with the theoretical expectation that attractive interactions and finite Zeeman splitting favor the stabilization of triplet pairing. However, under repulsion (panel j), the triplet amplitude is consistently small across all parameters, indicating that triplet pairing does not emerge in this regime. It also is insensitive to the Zeeman field in this regime. Remarkably, the predominance of triplet-$x$ pairing also suggests that the Zeeman field may induce a topologically nontrivial superconducting phase characterized by unconventional triplet pairing.

It is important to note that in the strong-coupling regime of Zeeman (see Section V.1.1), these pairing excitations are absent. So, only the $x$- and $y$-phases survive, with the $y$-phase dominating. In contrast, at finite Zeeman fields, we observe an enhancement in not only the CDW but also the $z$ and mixed components of the SDW correlations, as well as the singlet and triplet pairing correlations. Even more interestingly, at the relevancy boundary of the Zeeman gap, corresponding to the weak Zeeman case discussed in Section V.2.2, the Zeeman-induced logarithmic correction to the correlations notably suppresses the $y$-phase while enhancing the pairing correlations.

V.1.4 Logarithmic Corrections Arising from Pairing

In this section, we calculate the logarithmic corrections to the main correlation functions arising from the marginal pairing gap. We rewrite the RG flow equations (47) and (48) in terms of the magnetized-interaction parameter $\mathcal{Y}_{\vartheta_{k_{F}}}$ as

[TABLE]

The above equations are obtained by expansion up to second order in the coupling constants $\mathcal{Y}$ and $\mathcal{Y}_{\vartheta_{k_{F}}}$. To find the RG flow trajectories, we multiply Eq. (68) by $\mathcal{Y}(l)$ and Eq. (69) by $\mathcal{Y}_{\vartheta_{k_{F}}}(l)$, leading to the conserved quantity

[TABLE]

Here, $A_{z}=0$ corresponds to the boundary between relevant and irrelevant gap regimes, while $A_{z}\neq 0$ characterizes exact regimes. Inserting this constant of motion into Eqs. (47) and (48) and integrating yields

[TABLE]

Using these solutions, we evaluate the pairing-induced corrections to the density correlation functions as

[TABLE]

where

[TABLE]

The corresponding correction for pairing correlations reads

[TABLE]

with

[TABLE]

At the marginal boundary $\mathcal{Y}_{\vartheta_{k_{F}}}=\mathcal{Y}$, the solutions simplify to

[TABLE]

and in the limit $r\to\infty$,

[TABLE]

Since the marginal boundary of superconductivity shifts towards weaker repulsive interactions with increasing Zeeman strength, these logarithmic corrections appear in the correlations for weaker repulsion. According to Eqs. (73)–(77), all density wave correlations exhibit the same power-law decay characteristic of the topological insulator edge, with superconductivity introducing an incremental logarithmic correction. In contrast, pairing correlations, which decay faster, receive decreasing logarithmic corrections from superconductivity. This suggests that in a 1D topological superconductor with relevant Zeeman and marginal superconductivity, density wave excitations dominate, and superconductivity enhances their strength.

In the irrelevant regime ($A_{z}l\gg 1$), using the general solutions (71) and (72) to integrate $L_{1}$ and $L_{2}$, we find

[TABLE]

No logarithmic corrections appear, and the results resemble the non-perturbative case but with additional coefficients depending on the magnetized interaction parameter and the renormalized Luttinger parameter:

[TABLE]

Here, $l=\log\frac{r}{a_{0}}$ defines the characteristic length scale separating regimes dominated by logarithmic corrections and non-perturbative renormalization of parameters.

V.1.5 Pairing-Induced Correction of Short-Range Correlations in the PMH Edge

In this subsection, we calculate corrections to the short-range components of the density correlation functions arising from perturbation theory up to second order in the superconducting gap parameter. This approach allows us to study the effects of an irrelevant gap (in the RG sense) on the residual parts of the correlation functions. Ultimately, the full correction to each correlation function is expressed as a superposition of logarithmic corrections (treated in the previous section) and these residual corrections.

Since the logarithmic corrections were derived in real space and time, we combine those results here with Fourier space calculations. Corrections to the first term of the CDW operator (Eq. (50)) and the SDW operator in the $x$ direction (Eq. (51)) are defined as

[TABLE]

and

[TABLE]

The expected values in these expressions can be evaluated by known correlators of the form

[TABLE]

and

[TABLE]

As shown in Refs. Bakhshipour2024 ; Citro2020 , after evaluation the correction terms assume the well-known form involving Hartree and Fock averages:

[TABLE]

[TABLE]

It can be shown that the Hartree term does not contribute to the correction when $\tilde{K}^{-1}>2$, whereas the Fock term exhibits physically relevant behavior in this regime. Specifically, the short-range correlation corrections in the frequency space behave as

[TABLE]

Notably, the charge density phase at $\tilde{K}^{-1}>2$ is unaffected by superconducting corrections, so the correction term becomes negligible as $\omega\to 0$. For $\tilde{K}^{-1}<2$, corrections become significant. The SDW phase in the $x$ direction undergoes superconducting corrections for $\tilde{K}^{-1}<5/2$, indicating competition between helicity-aligned SDW short-range order and long-range density wave excitations under pairing. The $y$-component of SDW behaves similarly to CDW but is weighted by $\sin^{2}(2\vartheta_{k_{F}})$.

In topological superconducting PMH states with relevant Zeeman and marginal superconductivity, long-range density wave excitations under repulsive interactions are enhanced by superconductivity. For the short-range $x$-component SDW, the relevancy boundary shifts toward repulsive interactions beyond $1$, implying that stronger Zeeman fields are needed to suppress superconductivity in this channel, consistent with discussions in Section V.1.2.

V.2 Pairing-induced superhelical state perturbed by Zeeman effect

In this part, we consider the pairing gap to be relevant. This means that the Zeeman gap dwells near the boundary of its weak regime. We then calculate the logarithmic correction of this marginal Zeeman gap to the superhelical wire. Here, the field operators are taken with $\vartheta_{k_{F}}=0$. The effective superconducting helical state Hamiltonian $\mathcal{H}=\mathcal{H}_{s-hel}+\mathcal{H}_{z}$ reads

[TABLE]

The superconducting Luttinger liquid parameter $\bar{K}$ and velocity $\bar{v}$ are given by

[TABLE]

[TABLE]

Here, $\bar{v}_{F}=\frac{\epsilon_{F}\pm\Delta_{s}}{\pm k_{F}}$ is the effective Fermi velocity.

The RG flow equations for the Zeeman coupling and the superconducting Luttinger parameter can be obtained as

[TABLE]

where $\mathcal{Y}_{z}=\frac{\sqrt{2}\Delta_{z}a_{0}}{\bar{v}}$ is the dimensionless Zeeman gap coupling. According to Eq. (102), the Zeeman coupling $\mathcal{Y}_{z}$ becomes relevant when $\bar{K}<2$. It is observed that, unlike Eq. (15), the scaling of the Zeeman coupling coefficient here depends on the pairing-renormalized Luttinger parameter, which itself is influenced not only by the electron-electron interactions but also by the superconducting gap.

Moreover, we find a relation between the Luttinger parameter $\bar{K}$ and the interactions by defining the pairing-gap-tuning interaction parameter $\bar{\mathcal{Y}}_{int}$ as

[TABLE]

where $\bar{\mathcal{Y}}_{int}=1+y_{bf}$. Here, in the absence of a pairing gap, as discussed in Section III, the Zeeman term becomes relevant for both weak attractive and repulsive interactions with strengths less than 1. However, once superconducting pairing sets in, these boundaries shift.

In Fig. 7, the dependence of the superconductivity-renormalized Luttinger parameter $\bar{K}$ on the pairing gap $\Delta_{s}$ and the interaction parameter $y_{bf}$ is shown for the repulsive regime, $y_{bf}>0$, (panel a) and the attractive regime, $y_{bf}<0$, (panel b). As shown in Fig. 7(a), increasing the repulsive interaction strength reduces $\bar{K}$ so that it takes values below unity. This indicates that repulsion suppresses $\bar{K}$, resulting in stabilization of the Zeeman gap in the system. Increasing $\Delta_{s}$, slightly reduces $\bar{K}$. In contrast, Fig. 7(b) illustrates the influence of attractive interactions on the Zeeman gap. As the magnitude of the attractive interaction increases, the value of $\bar{K}$ increases, so that it reaches values above unity. Similarly, a larger $\Delta_{s}$ increases the value of $\bar{K}$. Therefore, in this regime, the system becomes unstable toward the formation of superconductivity over magnetism.

V.2.1 Strong Pairing with Weak Zeeman Field

A comprehensive view of the phase diagram requires exploring the strong-coupling regime of pairing. Notably, in many systems—such as the Kitaev model—this regime naturally leads the system into a trivial topological phase Kitaev2001 ; Shi2023 ; Cinnirella2025 . In proximity-induced systems, such as nanowires or two-dimensional heterostructures, a sufficiently strong parent superconductor may induce strong pairing in the wire. Investigating this regime becomes particularly important in the presence of magnetic fields or disorder.

According to our analysis, in the strong pairing limit, the bosonic field $\Theta$ becomes ordered. As a result, all power-law correlations involving the dual field $\Phi$ are suppressed. The CDW operator, which involves the gradient of the bosonic field $\Phi$, benefits from its strong fluctuations and maintains its Fermi-like behavior with parameters renormalized by superconductivity. The SDWxx vanishes, while the singlet and triplet-$x$ superconducting correlations still exhibit power-law decay, though slower than before, eventually saturating to a finite value.

Topologically, in contrast to the weak pairing regime, where the system enters a topological superconducting phase with Majorana edge states and dominant singlet correlations, in the strong pairing regime, the system typically flows into a trivial topological phase. Due to strong pairing, the fermions become locally bound and no edge states emerge, leading to dominant short-range and local correlations. Interestingly, in our study, despite the strong coupling, the singlet and triplet-$x$ superconducting correlations remain the dominant ones.

V.2.2 Logarithmic Corrections Arising from Zeeman

We use the pairing-gap-dependent interaction parameter, denoted as $\bar{\mathcal{Y}}_{int}$, and rewrite the RG flow equations for the Zeeman coupling $\mathcal{Y}_{z}$ and the interaction parameter by eliminating the Luttinger parameter $\bar{K}$. Keeping terms up to second order, the flow equations read:

[TABLE]

From these equations, we can obtain a constant motion defining the marginal Zeeman boundary in the presence of superconductivity:

[TABLE]

At this boundary (the separatrix), logarithmic corrections to the density and pairing correlation functions arise. The correction to the SDW correlation along the $y$ obtains

[TABLE]

The $zz$ component and mixed components $\mathcal{R}^{yz,\mathrm{log},\mathrm{Ze}}{SDW}(r)=-\mathcal{R}^{zy,\mathrm{log},\mathrm{Ze}}{SDW}(r)$ behave similarly but with different coefficients. The pairing correlations along the separatrix behave as

[TABLE]

At this superhelical-Zeeman boundary, the superconducting coupling tends toward strong coupling, and the dominant superconducting phase undergoes additional logarithmic corrections due to the Zeeman field. As a result, SDW correlations in the helical superconductor are suppressed by the Zeeman effect, while the Zeeman field incrementally enhances the pairing correlations, stabilizing a stronger helical topological superconducting phase.

V.2.3 Zeeman-Induced Corrections to Short-Range Correlations in the Superhelical Edge

The corrections to the short-range components of the correlation functions induced by the marginal Zeeman gap can be approximated in frequency space as

[TABLE]

At low frequencies, the Zeeman term contributes significantly to the leading behavior on the charge phase only when $\bar{K}<5/2$. In this situation, increasing the superconducting gap in the presence of repulsive interactions drives $\bar{K}$ to smaller values, which deepens the Zeeman gap. As a result, the enhancement of the superconducting gap makes the Zeeman-induced corrections to the charge sector more pronounced, making this phase more competitive. In the $x$-phase, a notable correction appears when $\bar{K}<2$. This analysis demonstrates that, under short-range interactions, the charge density wave phase can strongly compete with the long-range correlations discussed earlier.

VI Spin Transport in the Magnetized Superconducting Helical Edge

In this section, we investigate the spin transport properties of the helical edge state in the presence of both Zeeman and superconducting gaps. The Zeeman gap modifies the carrier velocity and, due to the helical feature of the states, mixes the spin state, in particular near $k=0$. In the presence of superconducting pairing, the additional spin gap opens through the sine-Gordon term which depends on the spin bosonic field $\Theta$ that may impact spin transport. In what follows, in order to analyze these effects, we calculate the spin conductivity employing the Memory function formalism Giamarchi1991 . Furthermore, to capture behavior over a broad temperature range, we present RG-improved Memory function results.

VI.1 The frequency dependence of spin conductivity

Based on the Kubo formula for conductivity Mahan , the conductivity can be written as Giamarchi1991

[TABLE]

where

[TABLE]

is the Memory function, and $\chi(\omega,T)$ is the current-current correlation function. Following Ref. Giamarchi2004 , we find the Memory function at $T=0$ as

[TABLE]

Substituting Eq. (113) into Eq. (111) shows that the frequency dependence of the conductivity takes the form

[TABLE]

By taking the real part of the conductivity, one arrives at

[TABLE]

where $\mathcal{D}$ is the stiffness parameter and the regular part is given by

[TABLE]

VI.2 The temperature dependence of spin conductivity

Using a similar approach to the previous section and following Refs. Giamarchi2004 ; Visuri2020 , we find the spin conductivity at $\omega=0$ as

[TABLE]

One can explicitly see that the temperature dependence of the conductivity follows a nonuniversal power law, with the exponent controlled by interactions and the Zeeman gap, in contrast to other works.

It is important to note that the unusual power-law behavior obtained for the spin conductivity is valid only in the high-temperature limit. To extend our results to lower temperature range, we invoke the RG analysis. Considering the RG flow of the gap parameter in Eq. (47), this parameter scales as

[TABLE]

Here, we assume that $\tilde{K}(l)$ varies slowly with $l$. In this scenario, temperature serves as the energy cutoff, so the renormalization should proceed down to a cutoff of the order of the thermal length, $L_{T}\sim\frac{v}{T}$, corresponding to $e^{l^{*}}\sim\frac{L_{T}}{a_{0}}$. By stopping the scaling at $L_{T}$, we use the renormalized values $\Delta(L_{T})$, $\tilde{K}(L_{T})$, and $\tilde{v}(L_{T})$ as the temperature-dependent quantities in the Memory function and in the conductivity, instead of the bare parameters $\Delta$, $\tilde{K}$, and $\tilde{v}$. The RG-improved expression for the conductivity then reads

[TABLE]

where $\mathbf{B}(x,y)$ is the Beta function. Note that, at high temperatures, the gap coefficient acts as a constant value and does not undergo rescaling.

Numerically evaluated spin conductivity versus temperature is shown in Fig. 8(a) and 8(b) for repulsive and attractive interactions, respectively. Each diagram displays the spin conductivity of the super-PMH edge across three temperature regimes: the high-temperature regime ($T\gg\Delta$, dashed lines), the intermediate regime just above the gap ($T\gtrsim\Delta$, dot-dashed lines), and the low-temperature regime ($T<\Delta$, solid lines). Conductivities in these regimes are obtained using the standard memory function method, RG-improved memory function, and an exponential decay model for electrons, respectively. All quantities are shown in dimensionless form, with the Boltzmann constant $k_{B}$ set to unity.

In both cases, the spin conductivity increases with decreasing temperature when $T\gg\Delta$. As $T$ approaches $\Delta$, the gap coefficient $\Delta=\Delta_{s}\sin[2\vartheta_{k_{F}}]$ becomes significant, leading to a marked reduction in conductivity. For $T<\Delta$, the memory function approach is no longer valid and the conductivity of free quasiparticles is exponentially suppressed as $\sigma_{s}(T)\propto e^{-\Delta/T}$ Visuri2020 ; Hsu2018 ; Rice1976 ; Nattermann2003 .

Two further trends are apparent: first, the spin conductivity decreases systematically with increasing Zeeman field in both regimes. Second, for attractive interactions [Fig. 8(b)], the Zeeman gap substantially enhances the superconducting gap, resulting in a much sharper drop in conductivity—approximately three orders of magnitude lower than in the repulsive case [Fig. 8(a)] at temperatures of the order of the gap.

VI.3 Spin Conductance

By considering the length parametric regime $L_{T}\ll L_{\Delta},L$, we calculate the spin conductance of a finite edge connected to non-interacting leads in the vicinity of conventional superconductivity. Since the leads introduce resistance at the contact with the 1D edge, this value must be included in the calculation of the total conductance, in addition to the resistance of the wire itself Visuri2020 :

[TABLE]

The spin conductance is then given by

[TABLE]

where we have set $R_{\text{contact}}=1$ and used the approximation that $R_{\text{edge}}$ is very small. The resistance is related to the conductivity by the specific resistivity $R=\rho L$, where $\rho=1/\sigma$. This yields the following expression for the conductance:

[TABLE]

It can be seen that when $\tilde{K}^{-1}>3/2$, decreasing the temperature reduces the correction of the conductance. In contrast, when $\tilde{K}^{-1}<3/2$, the gap term becomes relevant in the RG sense, and lowering the temperature increases the correction, resulting in a decrease in the conductance of the super-PMH system. This shows that strong attractive interactions significantly enhance the effectiveness of the superconducting gap. In this regime, increasing the Zeeman gap reduces the threshold of interaction strength needed for the superconducting gap to affect the conductance correction. This means that since a stronger Zeeman field in the presence of attractive interactions further enhances the superconducting gap, even relatively weak attraction can produce the same impact on conductance as would otherwise require much stronger interactions. As discussed in Section V.1.2, for $3/2<\tilde{K}^{-1}<2$, the system experiences weak repulsion, and the edge state is in a relevant gap regime. However, in this case, the superconducting gap does not reduce the conductance (Eq. (122)). Moreover, increasing the Zeeman gap further broadens the regime where the superconducting gap is ineffective in reducing the conductance, making successful conductance correction possible just at much weaker electron repulsion.

Figure 9 illustrates the temperature dependence of the spin conductance for the super-PMH edge in various interaction regimes. In panels (a) and (b), we show the normalized spin conductance $G_{s}$ as a function of temperature for repulsive ($\widetilde{K}^{-1}>1$) and attractive ($\widetilde{K}^{-1}<1$) interactions, respectively, for various values of the Zeeman gap. For strong repulsive interactions $\widetilde{K}^{-1}>3/2$ (see Fig. 9(a)) the conductance remains close to its quantized value as temperature decreases, reflecting the irrelevance of the superconducting gap and nearly ideal spin transport. In contrast, as the interaction becomes more attractive $\widetilde{K}^{-1}<3/2$ (see Fig. 9(b)), the superconducting gap term becomes relevant in the renormalization group sense. In this regime, lowering the temperature leads to a pronounced suppression of spin conductance, consistent with the opening of a spin gap in the excitation spectrum.

Importantly, increasing the Zeeman gap in the presence of attractive interactions lowers the critical interaction strength necessary for the superconducting gap to dominate. As a result, even relatively weak attractive interactions can substantially suppress the spin conductance, similarly to the effect produced by stronger attraction at lower Zeeman fields. This interplay is manifested in the shifting conductance curves in Figure 9 as functions of both the interaction and Zeeman field.

VII Summary and Conclusions

In this work, we have investigated the interplay between electronic correlations, Zeeman coupling, and s-wave superconducting proximity in 1D helical edge states. Employing bosonization and RG methods, we have analyzed how the competition between magnetic and superconducting gaps, together with electron-electron interactions, determines the dominant ordering tendencies and correlation functions in both the PMH and superhelical regimes.

Our analysis shows that repulsive interactions generally favor magnetic order, while attractive interactions enhance superconducting pairing. The relevance or irrelevance of the Zeeman and pairing gaps, imposed by the Luttinger parameter, leads to a rich phase diagram in which the dominant phase can be tuned by controlling the interaction strength, the Zeeman field, or the superconducting proximity.

We have derived analytical expressions for the density and spin correlation functions, identifying regimes where CDW, SDW, singlet, or triplet pairing correlations dominate. In particular, we have shown that logarithmic corrections—arising from the competition between marginal or nearly-marginal gap terms—can substantially modify the scaling of correlation functions and influence the stability of different phases. We further analyzed the corrections to short-range correlations using perturbation theory, and clarified the conditions under which the charge sector or triplet superconducting order can become dominant. Our results indicate that the coexistence of CDW and triplet pairing correlations can be stabilized by appropriately tuning the Zeeman gap and electron-electron interactions.

The spin transport properties of the helical edge is studied via the Kubo and Memory function approaches. We found that the frequency and temperature dependence of the spin conductivity is influenced by the interaction-dependent exponents. It is shown also that the RG flow of the gap parameters leads to significant renormalization of the low-energy transport. Our results for the spin conductance of a finite edge connected to non-interacting leads also indicate the role of contact resistance and show how superconductivity and magnetism along with mediating interactions can combine to modify the observable response.