Weyl Fermions and Broken Symmetry Phases of Laterally Confined $^3$He Films

Broken symmetries in topological condensed matter systems have implications for the spectrum of Fermionic excitations confined on surfaces or topological defects. The Fermionic spectrum of confined (quasi-2D) $^3$He-A consists of branches of chiral edge states. The negative energy states are related to the ground-state angular momentum, $L_z = (N/2) \hbar$, for $N/2$ Cooper pairs. The power law suppression of the angular momentum, $L_z(T) \simeq (N/2)\,\hbar\,[1 - \frac{2}{3}(\pi T/\Delta)^2 ]$ for $0 \le T \ll T_c$, in the fully gapped 2D chiral A-phase reflects the thermal excitation of the chiral edge Fermions. We discuss the effects of wave function overlap, and hybridization between edge states confined near opposing edge boundaries on the edge currents, ground-state angular momentum and ground-state order parameter of superfluid $^3$He thin films. Under strong lateral confinement, the chiral A phase undergoes a sequence of phase transitions, first to a pair density wave (PDW) phase with broken translational symmetry at $D_{c2} \sim 16 \xi_0$. The PDW phase is described by a periodic array of chiral domains with alternating chirality, separated by domain walls. The period of PDW phase diverges as the confinement length $D\rightarrow D_{c_2}$. The PDW phase breaks time-reversal symmetry, translation invariance, but is invariant under the combination of time-reversal and translation by a one-half period of the PDW. The mass current distribution of the PDW phase reflects this combined symmetry, and originates from the spectra of edge Fermions and the chiral branches bound to the domain walls. Under sufficiently strong confinement a second-order transition occurs to the non-chiral ``polar phase'' at $D_{c1} \sim 9\xi_0$, in which a single p-wave orbital state of Cooper pairs is aligned along the channel.


I. INTRODUCTION
The superfluid phases of liquid 3 He are BCS condensates with topologically non-trivial ground states. The topology of the ground state of 3 He has sparked research into the properties of liquid 3 He in confined geometries. [1][2][3][4][5][6][7][8][9] In parallel with new theoretical investigations, the physical properties of superfluid 3 He confined in one-or two dimensions have become accessible through a combination of nano-scale device fabrication, precision NMR and transport measurements at low temperatures. [10][11][12] For thin films, or 3 He confined in thin rectangular cavities, the ground state of 3 He is the chiral ABM state, 11,13,14 defined by the spin-triplet, p-wave order parameter, d(p) = ∆ẑ (p x ± ip y ), for a condensate of Cooper pairs with orbital angular momentum L z = ±h quantized along the normal to the film surface, and spin projection S z = 0. 15,16,1 This phase breaks timereversal symmetry as well as parity, and is realized at all pressures below melting pressure in films with thickness w ≲ 300 nm. Experimental confirmation of broken time-reversal and mirror symmetry of 3 He-A was recently made by the observation of an anomalous Hall effect for electron bubbles moving in 3 He-A films. 17,18 Thus, thin films and confined geometries are ideal for investigating the topological excitations of the chiral phase of superfluid 3 He.
The ABM state, and its generalization to higher angular momenta, is also a model for topological superconductors with a chiral ground state, including the model for the superconducting state of Sr 2 RuO 4 originally proposed as an electronic analog of superfluid 3 He-A, 19,2 as well as * sauls@lsu.edu 1 The spin-quantization axis is aligned with the orbital angular momentum in order to minimize the nuclear dipolar energy. 2 The identification of the pairing symmetry and ground state of Sr 2 RuO 4 , including whether or not it is chiral, as well as whether it is spin-triplet or spin-singlet, has been a focus of several research groups. [20][21][22][23][24][25][26] the multi-component, heavy fermion superconductor, UPt 3 , which exhibits broken time-reversal symmetry in the lowtemperature B-phase; a spin-triplet, chiral order parameter of the form, 27 d(p) = ∆ẑ p z (p x +ip y ) 2 , is consistent with Josephson interferometry 28 and polar Kerr effect measurements, 29 and is predicted to exhibit chiral edge states, an intrinsic thermal Hall effect, 30 and other novel transport properties. 24 There is also broad interest in inhomogeneous superconducting ground states, often associated with competing orders or pairing breaking mechanisms, 8,[31][32][33] as well as novel pair density wave phases in chiral superconductors 34 and paired fractional quantum Hall fluids. 35 Here we predict novel inhomogeneous ground states to emerge, via a sequence of phase transitions, when the chiral phase of superfluid 3 He is confined along one direction of its orbital motion. Confinment may also provide a technique for creating novel ground states in unconventional and topological superconductors.
Much is now understood about topological superfluids and superconductors, particularly superfluid 3 He. 36,37 Many new phases of confined 3 He with novel broken symmetries have been predicted based on Ginzburg-Landau theory that incorporates strong-coupling corrections to weak-coupling BCS theory for spin-triplet, p-wave pairing and boundary scattering and confinement. [38][39][40][41] For recent reviews of broken symmetry phases of superfluid 3 He in confined geometries and emergent topology in these phases see Refs. 8,9, for recent experimental studies of new phases of confined superfluid 3 He. In the 2D limit the chiral ABM state is fully gapped and belongs to a topological class related to that of integer quantum Hall systems, but with a topological index defined in terms of the Nambu-Bogoliubov Hamiltonian. [45][46][47] The topology of the 2D chiral ABM state requires the presence of gapless Weyl fermions confined on the edge of a thin film of superfluid 3 He-A, or a domain wall separating degenerate topologically distinct phases with opposite chirality. 48 For an isolated boundary a single branch of Weyl fermioms disperses linearly with momentum p || along the boundary, i.e. ε(p || ) = c p || , where the velocity of the Weyl Fermions, c = v f (∆/2E f ), is much smaller than the group velocity of Fermions in the normal phase of the Fermi liquid.
Broken time-reversal symmetry is reflected in the asymmetry of the Weyl branch under time-reversal, ε(−p || ) = −ε(p || ), which implies the existence of nonvanishing ground-state edge current derived from the occupation of the negative energy states. For superfluid 3 He-A confined in a thin cylindrically symmetric cavity, the edge sheet current, J = 1 4 nh, is the origin of the ground-state angular momentum, L z = (N/2)h, predicted by McClure and Takagi based on symmetry properties of an N-particle BCS ground-state of chiral p-wave pairs confined by a cylindrically symmetric potential. 49 Many unique features of 3 He-A, particularly anomalous transport properties, are directly related to Weyl Fermions confined on edges, domain walls, vortices or bound to electron bubbles. 17,18,[50][51][52] We consider the effects of hybridization between distinct branches of Weyl Fermions confined on neighboring boundaries and domain walls. We focus on thin (thickness w ≪ ξ 0 ) films of 3 He-A that are laterally confined, as shown for example in Fig. 1; D is the lateral confinement dimension. In Sec. II we relate bound-state formation to the multiple scattering problem of Nambu-Bogoliubov Fermions in a laterally confined chiral superfluid film. We develop a quasiclassical analysis for Bogoliubov quasiparticles in a laterally confined film, map the multiple-reflection of wave packets onto ballistic porpagation of wavepackets through a periodic array of domain walls separating degenerate, time-reversed domains of chiral order parameter. The sub-gap and continuum spectra exhibit a band-structure resulting from hybridization of counterpropagating Weyl Fermions on opposite edges. In Sec. III we employ a quasiclassical Green's function formalism to calculate the local density of states (DOS), resolved in both energy and momentum, as well as the edge currents and order parameter. Self-consistent solutions for the order parameter and spectral function (Sec. IV) are essential for obtaining and interpreting the ground-state of 3 He films under strong confinement. In section V, we show that hybridization of the Weyl branches leads to suppression of the chiral ABM state, and a reduction of the edge current. Even weak lateral confinement (D ≫ ξ 0 where ξ 0 =hv f /2π k B T c ≈ 200 − 800Å is the superfluid correlation length) suppresses the chiral ABM state and stabilizes a non-chiral polar state, d(p) = ∆ẑ p y for T ≲ T c . In Sec. VI we develop linear instability analysis to identify phase transitions and phase boundaries as a function of temperature T and lateral confinement D.
The phase diagram for 3 He films as a function of temperature and lateral confinement is obtained based on numerical calculations of the thermodynamic potential (Sec. IX) from selfconsistent solutions for the order parameter and spectral function. A central result of this report is the prediction of a pair density wave (PDW) phase of 3 He that spontaneously breaks translational symmetry, in addition to parity and time-reversal symmetries, and is stable over a wide range of temperature and confinement (Sec. VII). As a function of confinement, the PDW phase evolves from the polar phase as a single-Q mode chiral instability at D c 1 ≃ 9ξ 0 into to a periodic array of domain walls separating pairs of time-reversed chiral domains. The onset of the PDW phase from the translationally invariant chiral ABM state defines the upper critical confinement length, D c 2 ≈ 16ξ 0 . The structure, evolution and signatures of the PDW phase are discussed in Secs. VII-VIII. We begin with a brief introduction to the Fermionic spectrum and boundary current on an isolated edge of a thin film of 3 He-A.
for the two-component particle (u(r)) and hole (v(r)) spinor wavefunctions. The kinetic energy, ξ (p) = |p| 2 /2m * − µ, is measured relative to the chemical potential, µ, while∆(p) is the 2 × 2 spin-matrix order parameter, and p →¯h i ∇ in the coordinate basis. For the chiral ABM state the orbital order parameter is where {m ,n ,l} is an orthogonal triad defining the orbital orientation of the Cooper pairs with chirality ±1. The real vector d is the direction in spin space along which Cooper pairs have zero spin projection. For a dipole-locked 3 He-A filmd ∥l =ẑ. Thus, the spin state of the pairs, i⃗ σ σ y ·d = σ x , is the ESP state with equal amplitudes for spins polarized along +x and −x.
The 3D bulk A-phase has gapless Fermions for momenta along the nodal directions, p || ±l. 54 These zero-energy Fermions are protected by the integer phase winding of ∆ ± ∝ e ±iϕp about the chiral axisl. For thin films of 3 He-A, dimensional confinement discretizes the 3D Fermi sphere into the set of 2D Fermi disks; the Fermi surface is the set of boundaries of the Fermi disks. The resulting chiral phase is fully gapped on each disk. In the 2D limit the topology of 3 He-A can be expressed via the Bogoliubov Hamiltonian written in the momentum and iso-spin representation using the Nambu-Pauli matrices, ⃗ τ, for both ESP components, with c = ∆/p f , which emerge as the velocity of Weyl Fermions confined on a boundary of 3 He-A. The dispersion relation for excitations in the bulk of the chiral condensate is obtained from H 2 = |⃗ m(p)| 2 = ξ (p) 2 + c 2 |p| 2 . Note that |⃗ m(p)| ̸ = 0 for all p provided µ ̸ = 0. The point µ = 0 defines a quantum critical point separating a chiral BEC (µ < 0) and a chiral BCS (µ > 0) condensate. For 3 He-A the chemical potential is large and positive, µ = E f > 0, thus the low-energy excitations are confined near the Fermi surface with |p| ≈ p f , with bulk excitation energies ε (±) (p) = ± ξ (p) 2 + ∆ 2 .

Andreev's Equation
At a domain wall or an interface the Chern number changes abruptly. A discontinous change in the topological index forces the excitation gap to close and a spectrum of gapless Fermions to be confined on the boundary. In 3 He-A such a change in topology is also accompanied by the breaking of chiral Cooper pairs, with one or more components of the order parameter suppressed at the boundary. Thus, we consider a more general parametrization, ∆(r, p) = (∆ 1 (r) p x + i ∆ 2 (r) p y ) /p f , where the two order parameter components heal to their bulk values far from the boundary, ∆ 1,2 (r) → ∆. The Hamiltonian generalizes to for both ESP components defined byd =ẑ, with the pair potentials interpreted as symmetrized differential operators, A quasiclassical approximation to the Bogoliubov equations is obtained by separating the Bogoliubov spinor, | ϕ ⟩ ≡ (u , v) T , into fast-and slow spatial variations: | ϕ ⟩ = e ip fp ·r/h | ψ ⟩, where the quasi-classical spinor amplitudes, | ψ ⟩ ≡ (u p , v p ) T , vary in space on a length scale L ≳ ξ . By retaining the leading order terms inh/p f L ≪ 1 we obtain Andreev's equation, 55,56 ih where ∆(p, r) = ∆ 1 (r)p x + i ∆ 2 (r)p y ,p x,y are the direction cosines of the Fermi momentum, p = p fp , and Fermi velocity, v p = v fp . The latter define straight-line trajectories in classical phase space for the propagation of wavepackets of Bogoliubov Fermions, i.e. coherent superpositions of particles and holes with amplitudes given by | ψ ⟩.

Weyl Fermions on the edge
At an edge boundary, or a domain wall, the topology of 3 He-A changes abruptly leading to a spectrum of gapless Weyl Fermions propagating along the edge with dispersion relation where the ∓ sign corresponds to the chirality, ±l, of the ground state. A T = 0 the negative energy states are occupied, and generate a ground-state edge current, j(x ⊥ ), parallel to the edge, where x ⊥ is the distance into 3 He-A normal to the surface. The total sheet current, including the back action from the negative energy continuum to the formation of the gapless edge states, is given by 4 where n is the density of 3 He atoms, andê ∥,⊥ are unit vectors tangent (∥) and outward normal (⊥) to the boundary. The sign of the current is defined by the chirality of the ground state. For a thin film (w ≪ ξ 0 ) of 3 He-A with chirality ±ẑ confined in a cylindrical cavity with a macroscopic radius R ≫ ξ 0 , a sheet current circulates on the boundary given by J = ± 1 4 nhê ϕ . This current gives rise to the ground-state angular momentum, where N is the total number of 3 He atoms. 4,57 Symmetry Protection of the Edge Current McClure and Takagi (MT) obtained the above prediction for the ground-state angular momentum of 3 He-A by considering an N-particle BCS ground state (condensate) of chiral p-wave pairs. 49 For a cylindrically symmetric confining potential they showed that a condensate wave function of N/2 chiral pairs is an eigenfunction of the angular momentum operator with eigenvalue L z = (N/2)h, and that this result is independent of the radial size of the pairs, i.e. the result is valid in both the BEC and BCS limits. Since a BEC of chiral p-wave molecules is topologically trivial (N 2D = 0), while the chiral BCS condensate is characterized by a non-trivial topological winding number N 2D = ±1, the MT result, L z = (N/2)h, is not simply a consequence of the topology of the bulk Hamiltonian in Eq. 4, nor is the edge current in the BCS condensate protected by topology. Rather, the MT result implies that the ground-state angular momentum of (N/2)h is protected by symmetry. Indeed, for a BCS condensate of chiral pairs confined by a nonspecular boundary that violates axial symmetry on any length scale L ≪ ξ 0 , the ground-state edge current, and thus the angular momentum, is reduced compared to the MT result. 4 By contrast, for a chiral BEC with a non-specular boundary the ground-state angular momentum is still (N/2)h because the angular momentum is confined within each molecule. Thus, the total angular momentum of the BEC is thenh times the volume integral of the condensate density which is insensitive to the boundary conditions. 58 Nevertheless, the topology of the chiral BCS condensate plays a key role in determining the origin of the ground-state angular momentum. The point µ = 0 separating BEC and BCS condensation is a phase transition in which the topology of the Bogoliubov Hamiltonian changes discontinuously. The change in topology is accompanied by "spectral flow" in which an odd number of branches of Fermionic excitations appear as µ is "adiabatically tuned" through zero from µ < 0 to µ > 0. 59,60 In the 2D limit the spectral flow yields a single branch of chiral Fermions, which disperses through the Fermi level. The zero energy mode is protected by the non-trivial topology of the bulk Hamiltonian, while the occupation of the negative energy states leads to a ground-state edge current. However, the magnitude of that current depends on the surface boundary conditions, and thus the symmetry of the confining potential, and is equal to the MT result only for specular reflection by a cylindrically symmetric potential.

II. LATERALLY CONFINED 3 HE FILMS
A key feature of the edge current is that the direction of the edge current is defined by the chirality and the outward normal to the surface, i.e. J ∼l ×ê ⊥ . Thus, for a toroidal geometry the currents counter propagate on the inner and outer edges as shown in Fig. 1. If the inner and outer edges are both specularly reflecting, then J 2 = −J 1 = 1 4 nhê ϕ , and the ground-state angular momentum is then L = e r × 2π R 2 1 J 1 + 2π R 2 2 J 2 w = (N/2)hẑ. Thus, we recover the MT value as a result of compensation of the angular momentum from the outer edge current by that from counter-propagating inner edge current. However, this result can only be valid when there is no overlap between the states on the inner and outer edges. So, "what is the effect of lateral confinement (D = R 2 − R 1 ) on the ground state of 3 He films?", and at "what length scales D are the ground state properties of the film strongly modified?" In what follows we neglect the effects of boundary curvature and consider a thin film, laterally confined in a narrow channel with parallel edges separated by distance D as shown Fig. 2. Our analysis also assumes zero magnetic field. However, all of the ground states of laterally confined, thin 3 He films that we report are equal spin pairing states. Thus, the magnetic field does not play a role in the relative stability of the different orbital ground states of thin 3 He films under lateral confinement. The nuclear Zeeman and dipole energies are however important for determining the NMR signatures of these phases, but this is outside the scope of this report.
We first investigate the effects of lateral confinement by solving Andreev's equation (Eq. 8) subject to boundary conditions connecting the Nambu spinor amplitudes defined on classical trajectories. Wave packets of normal-state particles and holes propagate at the Fermi velocity along classical, straight-line trajectories defined by the direction of the group velocityp. The wave packets undergo specular reflection at a translationally invariant boundary. 57,61,62 Thus, for specular reflection, p → p = p − 2n (n · p), where n is the normal to the boundary directed into 3 He, the Nambu spinors obey the conditions, for parallel boundaries as shown in Fig. 2 with x being the coordinate normal to the boundaries, and y is the coordinate along the channel, in which casep x = −p x andp y =p y . Note that we have assumed that translational symmetry is preserved along the laterally confined film by restricting the order parameter to depend only on x. Andreev's equations reduce to coupled scalar equations, where ∂ ≡hv fp · ∇ andṽ ≡ σ x v. where . Periodicity allows us to express the Andreev amplitudes as a Bloch wave solution of the form, where k is the Bloch wave number, U(x) and V (x) are periodic functions that satisfy and boundary conditions: The solution of Eqs. 16 and 17, with periodic boundary conditions, can be carried out (see Appendix A) to obtain the eigenfunctions, U(x) and V (x), and the eigenvalues, ε(p || , k), for both positive and negative energy branches of the spectrum, which depend on the conserved momentum along the channel, p || ∈ [−p f , +p f ], and the Bloch wave number, k ∈ [−π/2D, +π/2D], defined on the first Brilloun zone.
The branches of the sub-gap spectrum, |ε| < ∆, are the solutions of whereh v x λ (ε) = √ ∆ 2 − ε 2 and v x = v fpx . For k = 0 we obtain two linearly dispersing branches, corresponding to Weyl fermions with opposite chirality on the opposing boundaries (see Fig. 4). Hybridization of the pair of Weyl branches is evident for k ̸ = 0 with band formation of the levels at fixed p || ; the band edges are defined by the solutions of Eq. 18 at k = 0 and the Brilloun zone boundary, k = π/2D. The asymptote, ε + = + c p || , is the k = 0 dispersion for states confined on the upper edge x = D, while ε − = −c p || corresponds to states confined near x = 0. In the limit D ≫hv f /2∆, the maximum bandwidth is exponentially small, w ≃ ∆ e −2D ∆/hv f , which reflects the exponential decay of edge states into the channel. The bandwidth increases with confinement and is largest for states with p x = p f , i.e. trajectories normal to the confining boundaries. The bandwidths are exponentially small for trajectories nearly grazing incidence, i.e. p || ≈ p f , which follows from the divergence of the effective period, 2D eff = 2D/p x , for near grazing incidence trajectories.

III. PROPAGATOR AND SPECTRAL FUNCTION
By solving the Andreev equations, we obtain the sub-gap eigenstates for the confined film. The states are coherent superpositions of bound states confined at both edges with opposite group velocities, which decay into the channel region on the length scale ∼ ξ 0 . The overlap of the edgestate wave functions does not open a gap in the subgap spectrum; rather the isolated bound states with spectral function 2π∆δ (ε − ε ± ) develops into a band with a bandwidth dependent on both confinement length D and trajectory p. In what follows we discuss the origin of the additional states that form the sub-gap band structure, and the spectral weight redistribution under confinement. The edge-state spectral function is obtained from the local, retarded single-particle Green's function, which in the quasiclassical limit satisfies Eilenberger's transport equation, 63 and normalization condition, G 2 = −π 2 1, where ∆ is the pairing self-energy and G is the quasiclassical Green's function in Nambu space, and A , B = A B − B A. In this section we present results for the quasiclassical propagator in laterally confined 3 He-A thin films for a spatially uniform order parameter ∆(p). These results provide insight into the spectrum of the laterally confined chiral phase. Self-consistent numerical solutions of the quasiclassical transport equation are discussed in Secs. IV-IX to follow.
For a laterally confined film, the solution for G (see Appendix B) is where G 0 is the bulk propagator and G ± are the Nambu matrices representing the states generated by the edges. Coefficients C 1 and C 2 are obtained by matching boundary conditions at both edges for the trajectory p and its specularly reflected counterpart p. The factor M is given by The physical Green's function satisfies Eilenberger's normalization condition, which is guaranteed by the commuatation relations, The propagator G can also be expanded in Pauli matrices in The off-diagonal propagators, F 1,2 (r, p), contain information about the Cooper pair spectrum, while the diagonal component, determines the local fermionic density of states, where Λ ≡ √ ∆ 2 − ε 2 , ∆ 1 = ∆p x is thep x component of the order parameter for a constant gap amplitude, and ∆ 2 = ∓∆p y is thep y orbital component. Here and thoughout the manuscript the notation, ∆ 1,2 refers to the τ 1 ,τ 2 component of ∆, respectively. For the undistorted chiral A phase, ∆ 1 = ∆p x and ∆ 2 = ∓p y for the ground state with L z = ±h. For laterally confined 3 He the amplitudes of the two components differ. Thus, for the laterally confined A-phase with right-handed chirality, and withx defining the normal to the confining edge, is the amplitude of the orbital state that is normal (parallel) to the pair-breaking edge. We address the effects of anisotropy induced by pairbreaking in Secs. IV-IX.
In the limit, D/ξ 0 → ∞, the factor M → 1, and the τ 3 component reduces to For each p, the imaginary part of G ∞ for ε < ∆ only comes from two simple poles: One at ε = ∆ 2 with spectral weight proportional to e −2∆x/v f , corresponds to the edge states confined near x = 0, and the other pole The spectral function is shown as a function of parallel momentum p ∥ and energy ε. The energy ranges from −5∆ to 5∆, which contains both the bound-state spectrum shown in detail in Fig. 4, and the continuum |ε| > ∆. For each p ∥ , the continuum forms a band structure with band gaps generated by hybridization and multiple reflection.
band structure of sub-gap spectrum is encoded in the normalization factor M, which becomes imaginary when the righthand side of Eq. 22 is negative. For sub-gap states, The range in energy over which the sub-gap spectral density is nonzero is identical to the bandwidth of bound-state energies obtained from the solution of Andreev's equations discussed in Sec. II. Note that M not only modifies the edge contributions to the propagator G ± , but also the bulk term proportional to G 0 . Multiple reflections not only reorganize the bound states by re-distributing spectral weight among themselves, but the formation of sub-gap band is also connected with changes in continuum spectrum with |ε| > ∆. The spectral density at x = D for D = 8ξ 0 is shown in Fig 4 for the sub-gap spectrum, and in Fig 5 over a range of energies into the continuum. Figure 4 shows the sub-gap spectrum |ε| < ∆ at x = D. For each p ∥ , the chiral bound states form a band. The density of states is maximum at the upper and lower band edges, resulting from Van-Hove singularities of the one dimensional periodic pair potential. The bandwidth is largest for the normal incidence where the period is smallest. The spectral density is maximum for branch ε = cp ∥ , and exponentially small for the opposite chirality ε = −cp ∥ corresponding to the hybridization with Weyl fermions on the opposing edge x = 0. The continuum spectrum for |ε| > ∆ shown in Fig. 5, also exhibits energy bands. Finite band gaps develop for each p ∥ . The band structure in both sub-gap and continuum states results from the overlap of edges states, and the periodicity due to multiple reflection of quasiclassical trajectories.

IV. SELF-CONSISTENT SOLUTION
The existence of boundaries and the formation of edge states modifies the pair propagators from their bulk values, which then leads to a modification of the off-diagonal mean-field paring energy, or order parameter. Specifically, the order parameter component associated with the p-wave orbital that is normal to the boundaries is suppressed to zero at the edge; F 1 = 0 at both edges as can be seen from Eq. 34. The other orbital component will either be enhanced or suppressed depending on the conditions imposed by boundary scattering. For a specularly reflecting edge, the off-diagonal self-energies satisfy the following "gap equations", where the integration is over the Fermi surface, which in the 2D limit is angular integration around a circle of radius p f . The p-wave pairing interaction in the 2D limit, g(p, p ′ ) = 2g 1 p · p ′ , and the pairing bandwidth, ω c , are eliminated in favor of the bulk transition temperature, T c , using the linearized bulk gap equation for the pairing instability, 1/g 1 = K(T c ), where is a digamma function. More generally, Eq. 32 is used to regulate the log-divergent Matsubara sums in Eqs. 30, 31 in order to eliminate the cutoff and to express these functions in terms of the scaled temperature, T /T c , for i = 1, 2. The off-diagonal Matsubara Green's functions are obtained from the retarded propagators by analytic continuation (ε → iε n ). For constant gap amplitudes, ∆ 1,2 , analytical results for F 1,2 are given by, order parameters, provide considerable insight into the subgap and continuum spectra, as well as semi-quantitative results for the spectral functions, order parameter and edge currents.
The gap equations (Eqs. 30, 31) are solved self-consistently with Eilenberger's transport equations. 64,65 The resulting gap profiles are shown in Fig. 6. There is strong pair breaking of the amplitude for the orbital component normal to the edges; ∆ ⊥ (x) vanishes at both edges and is suppressed throughout the confined channel. Pair breaking is also reflected in the spectral density. In Fig. 7 we show the local density of states for the self-consistently determined order parameter, as well as results for a suppressed, but spatially uniform amplitude, ∆ ⊥ = 0.5∆. As ∆ ⊥ (x) decreases the bandwidth centered around zero energy increases, and band gap to the continuum reduces. The basic structure of the spectrum, including the Van Hove singularities is already captured by the non-selfconsistent order parameter.

V. EDGE CURRENT SUPPRESSION
The propagator, G(p, x; ε), determines the mass current density, which in the Matsubara representation is given by The solution for G includes the current resulting from the subgap edge spectrum as well as the back action of continuum to their formation. The resulting integrated sheet current for a single edge is then J = 1 4 nh, at zero temperature. 4 Due to the chirality of the edge states the mass current flows in the positive +y direction near the edge at x = D, and in −y direction along the edge at x = 0. Thus, in the laterally confined geometry with two edges the current density vanishes at the mid point of the channel. Thus, we define the sheet current for a single edge as the integrated current density in half of the channel, The reduction of the edge current at any temperature can be expressed as |J| = 1 4 nh × Y (T, D). In the limit D → ∞, the reduction is due to thermal excitations and is given by, where the sum runs over the Matsubara energies, ε n = (n +  66,4 The edge current calculated self-consistently as a function of D and T is shown in Fig. 8. For weak confinement, D = 15ξ 0 , the current decreases monotonically with temperature similar to that for D → ∞; however, J(0, D) is suppressed relative to 1 4 nh, and J(T, D) vanishes at T c 2 = 0.89T c , i.e. below the superfluid transition at T c . The vanishing of the edge current for T c 2 < T ≤ T c indicates that the high-temperature superfluid phase under confinement is a non-chiral phase. For 10.4ξ 0 ≲ D < 12ξ 0 , the mass current onsets at a relatively low T c 2 , increases to maximum and then decreases. This nonmonotonic behavior implies that the chiral phase is also suppressed at low temperatures for sufficiently strong confinement. For D ≈ 10ξ 0 the low-temperature phase is also a nonchiral phase, which undergoes a transition to a chiral phase at finite temperature (T c 3 = 0.23T c for D = 10ξ 0 ), then re-enters a non-chiral phase at a higher temperature (T c 2 = 0.62T c for D = 10ξ 0 ). The temperature window, T c 3 < T < T c 2 , for a stable chiral phase decreases with confinement. Below a critical thickness, the mass current is zero and the superfluid phase is non-chiral at all temperatures.

VI. PHASE TRANSITIONS UNDER CONFINEMENT
Based on the supression of the edge currents, we conclude that for the confined channel, the chiral A phase is suppressed near T c in favor of a non-chiral phase, but for weak to moderate confinement undergoes a second order transition at T c 2 to a chiral phase. The only possible non-chiral superfluid phase that is not suppressed by lateral confinement is the polar phase with only one p-wave orbital component, i∆ ∥ p y , aligned along the direction of the channel. The component ∝ p x is completely suppressed. To investigate the polar to chiral phase transition, we carry out linear stability analysis starting from the polar state, i∆ ∥ p y . In particular, we allow for nucleation of additional order parameter components at a second-order transition, where |A 1 |, |A 2 | ≪ ∆ ∥ are perturbations. For specular boundary scattering, the amplitudes satisfy the boundary conditions: Eilenberger's equations are solved to linear order in A 1 and A 2 . In zeroth order, the polar phase order parameter and the Green's function are given by where ∆ ∥ is computed self-consistently by solving the nonlinear gap equation, It is convenient to separate the first-order corrections to the order parameter into real and imaginary terms The corresponding corrections to the propagator obey the linearized transport equation and normalization condition, The resulting Cooper pair propagators for the first-order per- turbations satisfy second-order differential equations,

Polar Chiral
I . (46) where ω n = ε 2 n + ∆ 2 ∥ p 2 y and ∂ = v p · ∇ is the directional derivative along trajectory p. The corresponding gap equation for ∆ (1) is given by

A. Polar to A phase transition
The linear instability equations have a solution corresponding to a second-order transition from the polar phase to the chiral A phase. The unstable amplitude is determined by Eq. 45 and ∆ (1) R , and is translationally invariant along the channel direction (Q y = 0), i.e. the only spatial dependence is in x direction imposed by the edge boundary conditions. The Fourier expansion of the unstable order parameter is and the solution of Eq. (45) becomes, The linearized gap equation (Eq. 33) can also be Fourier transformed, and reduces to the eigenvalue equation,

B. Polar to PDW phase transition
The reentrant behavior is a hint that an inhomogeneous phase, one that accomondates both polar and chiral order but breaks translational symmetry along the channel, may be favored over a range of lateral confinement dimensions D. 65 In order to investigate this possibility we consider a more general form of the order parameter for the linear instability analysis that breaks translational symmetry along the channel, and which includes both real and imaginary parts of ∆ (1) , i.e. spatially modulated currents. In particular, we start from Eq. 39, then represent the complex amplitudes, A 1,2 (r), in Fourier space as where A R,I 1,2 are the real (R) and imaginary (I) parts of A 1,2 . Equation (45) becomes, and similarly for Eq. (45), The resulting homogeneous equations for a R,I 1,2 decouple for each mode Q, but the amplitudes for the p x and p y Cooper pairs are in general coupled, where t = T /T c and the coefficients, I R,I i j (Q), are given by The first non-trivial solution to Eq. 54 as a function of T or D signals the nucleation of a new superfluid phase from the translationally invariant polar phase. The spectrum of eigenvalues can be expressed in terms of a spectrum of wavevectors, Q, that satisfy the eigenvalue equation obtained from the vanishing determinants for the real and imaginary components of the linearized gap equations, We find that for all temperatures, the real sector of Eqs. 57, with both Q x,y ̸ = 0, and Q y ∼ π 15ξ 0 , gives the smallest critical dimension, D c 1 , for the transition from the polar phase. The resulting phase spontaneously breaks the y-translation symmetry along the channel. At the instability point a single mode becomes unstable leading to spatial modulation of the order parameter defined by the wavevector, Q = (Q x , Q y ), of the unstable mode. In particular, the real part of the Cooper pair amplitude with orbital wave function p x develops a standing wave along the y direction, ∆ 1 ∝ sin(Q y y), i.e. a pair density wave (PDW). 65 The critical curve for the polar to PDW transition, D c 1 (T ) (blue curve in Fig. 10) pre-empts a polar to A-phase transition (red line), with Q y = 0 enforced, at all temperatures. Close to the transition at D c 1 (T ) the real part of the amplitude, ∆ 1 , of the p x orbital changes sign along the channel with period 2π/Q y (T ), whereas the p y orbital has a large imaginary polar component, i∆ ∥ p y , modulated by a small real component with the same period.

VII. PDW ORDER PARAMETER
For D ≳ D c 1 (T ), there is a proliferation of modes with Q y ̸ = 0 that become unstable. The resulting order parameter evolves continuously into a multi-Q PDW phase. The structure of the order parameter in the PDW phase requires a self-consistent analysis that incorporates the mode-mode coupling. As we show here the PDW order parameter evolves for increasing D > D c 1 (T ) towards a periodic array of chiral domains separated by domain walls. The latter are also chiral in that the domain walls support currents flowing along the domain walls.
The analysis is an extension of that described in Sec. IV for the order parameter of the laterally confined chiral Aphase. The long period of the PDW phase, combined with the fact that there is no bulk order within the confined geometry requires a large computational cell, L ≈ 100ξ 0 , with integration along multiply-reflected classical trajectories. We start with D ≳ D c 1 (T ) from a "seed" order parameter which is the polar phase with random noise added at the level of 10 −4 ∆ ∥ , which is unstable to the formation of the PDW phase. The resulting self-consistent order parameter for T = 0.7T c and D = 10.5ξ 0 is shown in Fig 11. It is parametrized as ∆(r, p) = ∆ ⊥ (r)p x + ∆ || (r)p y , where ∆ ⊥ (r) and ∆ || (r) are both complex amplitudes for the orbital components, p x and p y , respectively. The resulting order parameter is a periodic structure with a period L ≈ 30ξ 0 , which corresponds to the unstable wave vector, Q y ∼ π/15ξ 0 , predicted in single-mode instability analysis. Note that the maximum amplitude of ∆ ⊥ is at the scale of 0.08[2πT c ], around one-half of the bulk amplitude at that temperature. The phase changes by π across each domain wall centered at the half period. The resulting ordered phase is both chiral and periodic along the channel with the sign of chirality alternating each half period.
The structure and evolution of the PDW phase is most clearly shown in terms of slices of the order parameter amplitudes plotted as a function of the channel coordinate, y, for several lateral positions, x. Figure 12 shows the real and imaginary parts for both obital components: ∆ ⊥ (r) = ∆ R ⊥ (r) + i∆ I ⊥ (r) and ∆ || (r) = ∆ R || (r) + i∆ I || (r). All components exhibit periodic structure, reflecting the underlying single-mode instability with Q y determining the period near the Polar to PDW phase transition. The real part of the p x orbital, ∆ R ⊥ (r), develops large amplitude modulations and changes sign within one period. By contrast the imaginary part of the p y orbital amplitude, ∆ I || (r), is nearly constant, exhibiting small oscillations about what would otherwise be the bulk amplitude of the polar phase. The resulting structure is a periodic array of chiral domains, with the sign of the chirality switching every half period, i.e. an "anti-ferromagnetic" chiral phase topologically equivalent to . . .
Both orbital components ∆ ⊥ (r) and ∆ || (r) have much larger , which agrees with the decoupling of real and imaginary components in the linear stability analysis discussed in Sec. VI B and the fact that the first unstable mode occurs in the real sector of the order parameter fluctuations. The much smaller flucutaions of the imaginary components are driven by mode-mode coupling once the PDW phase is established for D ≳ D c 1 (T ). In particular, for the amplitude of the p y orbital, ∆ I || (r), exhibits only small oscillations about the amplitude of the parent polar phase. The amplitude ∆ R || is largely responsible for the periodic amplitude and phase modulations shown in the right panels of Fig 11. Note also that the amplitude for the p x orbital, ∆ ⊥ (r), is maximum in the center of a channel, and vanishes at the edges, while in the case of the real part of the amplitude of the p y orbital, ∆ R || (r), the converse is the case. These features are a result of specularly reflecting boundary conditions at both edges.
For weaker confinement D = 11.5ξ 0 at the same temperature, T = 0.7T c , self-consistent calculations of the order parameter lead to increased periods of the PDW phase. As we move well away from the Polar to PDW transition at D c 1 (T ), the structure of the order parameter develops soliton-like structures, with nearly homogeneous ordered regions confined between domain walls, the latter roughly 10ξ 0 in width. In this limit the structure of the order parameter is dominated by the coupling of many unstable modes with different wavevectors, Q.
The region of the channel shown in Fig. 13 spans 60ξ 0 ≈ 6µm and represents a half period of the PDW phase. The spatial profiles of the order parameter away from the lateral domain walls centered at y = 0 are identical with those of laterally confined chiral A phase with the sign of the chirality switching as one crosses the domain wall. The period increases dramatically with weaker confinement, and as D → D c 2 (T ) diverges towards the system size. Thus, the weak-confinement transition is the locus of points D c 2 (T ) at which it becomes energetically possible for a single domain wall to exist in the channel separating degenerate chiral ground states.
Finally, we note that there are similiarities in the energetics FIG. 13. Self-consistent order parameter amplitudes for T = 0.7T c and D = 11.5ξ 0 . The red curves correspond to the center at x = D/2. Blue curves correspond the edge at x = 0. The order parameters is shown for one-half of a period, and exhibits soliton-like domain walls as well as a "breather mode" confined on the domain wall. that lead to stabilization of the PDW phase in laterally confined very thin (quasi-2D) 3 He films considered in this report, and the "stripe phase" that develops from confinement in one dimension from bulk 3 He-B. 40,65 . In particular, the transition from the chiral A phase to the chiral PDW phase under lateral confinement of very thin 3 He-A films, and the transition from the B-phase of 3 He to the stripe phase, both result from the competition between the cost in energy of a domain wall separating degenerate states -2D 3 He-A with opposite chirality in laterally confined films and degenerate 3D B states of 3 He-B confined in one dimension -and the gain in condensation energy at the intersection of a domain wall and the boundariesconfining edges for the PDW phase or confining surfaces for the stripe phase. For both cases it is favorable for a domain wall to penetrate for sufficiently strong confinement in order to recover lost condensation energy on the boundaries. However, the structure of the PDW phase, and its broken symmetries, are distinctly different from that of the stripe phase of thick 3 He films. The PDW phase breaks time-reversal symmetry, while the stripe phase is time-reversal symmetric. Also, the transition from the PDW to the polar phase is second order, while the stripe to A phase is first order. There appears to be an additional similarity if we consider the currents that flow on the domain walls and near the confining boundaries -mass currents for the PDW phase and spin currents for the stripe phase. In the section that follows we discuss the structure of the chiral domain walls that define the PDW phase near D c 2 (T ), including the mass current near the intersection of the domain wall and the edges of the laterally confined film, and how current conservation, ∇ · j(r) = 0, is maintained. Silveri et al. discuss a similar situation that arises for spin-current conservation in the case of domain walls separating degenerate B-phases that define the onset of the stripe phase. 67 orbital across the domain (located at y = 0), i.e. ∆(y < 0) ∝ (−p x +ip y ) and ∆(y > 0) ∝ (p x +ip y ), and (ii) the p y -type DW where the sign of the amplitude of the p y orbital changes sign, ∆(y < 0) ∝ (p x − ip y ) and ∆(y > 0) ∝ (p x + ip y ). Even though the left (or right) chiral domains have the same Chern numbers for both types of domain walls, the DWs are physically distinct structures. For both DWs the Chern number changes from N = −1 to N = +1. Thus, both DWs support a branch of chiral Fermions and a ground-state current confined on the DW. However, the ground-state currents that flow along the p x -type DW and the p y -type DW at y = 0 are oppositely directed; the p y -type DW has current flowing along the domain wall in the direction consistent with chiralities of both chiral domains, while for the p x -type DW the current flows opposite the chirality. 50 Analysis of the energy per unit length of these two DWs shows that the p x -type DW has a lower free energy per unit length, and is thus the energetically favored DW, at least for weakly confined chiral domains. 50 Our analysis confirms that result. For a single p x -type DW in laterally confined 3 He-A, e.g. the PDW phase near the upper critical confinement length, D c 2 , the currents along the edges are determined by the chirality of each chiral domain. Thus, the energetically favored p x -type DW there appears to be a violation of current conservation at the junctions of the DW and the edges.

VIII. CHIRAL DOMAIN WALLS
The apparent violation of current conservation is resolved by the self-consistent solution for the order parameter and excitation spectrum for the p x -type single domain wall crossing the channel. A self-consistent calculation of the structure of a single, laterally confined p x -type DW is carried out for T = 0.5T c and D = 20ξ 0 . The results include the current density at every point within the channel. We initialize the order parameter to be −∆ ⊥ (x)p x + i∆ ∥ (x)p y in the region y < 0, and +∆ ⊥ (x)p x + i∆ ∥ (x)p y in the region y > 0, with the domain wall centered at y = 0. Note that ∆ ⊥,|| (x) are real for the solutions far from the DW. The order parameter is calculated by self-consistent solution of the Eilenberger and gap equations (Eqs. 20, 30, 31), subject to the asymptotic boundary conditions for y → ±∞ determined by the two possible degenerate  Figure 14 shows the current density in a region containing a single p x -type DW. As expected, the current confined near y = 0 flows in the +x direction with a magnitude that vanishes on approaching the edges at x = 0 and x = 20ξ 0 . Current conservation is maintained by the appearance of two branches of current flowing in the −x direction at distances of y ≈ ±5ξ 0 from the DW, smoothly connecting to the currents on both the top and bottom edges as shown in Fig. 14. Finally, we note that the complex current pattern associated with the laterally confined chiral DW includes supercurrents associated with the gradient of the global phase of the PDW order parameter that are induced in the central region of the channel. The resulting current density satisfies current conservation at every point within the computational cell, i.e. ∇ · J(r) = 0 to the accuracy of our numerics which was set by an order parameter residual no greater than 10 −4 at each point in the computational grid.
Finer detail of the current density of the chiral DW confined in the channel is shown in Fig. 15. The component j x exhibits the same sequence of current reversal on traversing the DW for all of the cuts defined by x. Indeed, the p x -type DW appears to support three branches of chiral Fermions, two branches with chirality ν = −1 centered at y = ±5ξ 0 , and a third branch with chirality ν = +1 centered at y = 0, for a net chirality of ∑ ν = −1. This is then consistent with the index we expect for the p y -type DW with a single chiral branch bound to the DW.

IX. STABILITY OF THE PDW PHASE NEAR D c 2
The self-consistent solutions for the order parameter, in a region where the chiral A phase has lower energy than the polar phase, yields the PDW phase with broken translational symmetry. The PDW phase spans a wide range of channel widths, D c 1 ≤ D ≤ D c 2 . The period of the PDW increases with increasing D (weaker confinement) and diverges as D → D c 2 (T ); i.e. the period becomes of order the system size. The upper critical channel dimension, D c2 (T ), can be defined as the value of D at which it costs zero energy to introduce a DW wall separating two otherwise degenerate domains of laterally confined 3 He-A. In the limit D ≲ D c 2 , the question is "how is it possible for a domain wall to lower the total energy of 3 He-A in a laterally confined geometry?" The answer is that in a laterally confined geometry there is competition between the DW energy and the pair-breaking energy -the occupied branch of chiral Fermions -at edges. Each edge can be shown to be mathematically equivalent to one-half of a p y -type DW. Furthermore, the energy cost of pair-breaking on the edge, equivalent to that of the p y -type DW, is reduced within several coherence lengths of the points of intersection of the lateral DW (p x -type) and the confining edges. Since the energy per unit length of the p y -type DW is greater than that for the p x -type DW, there is a critical channel dimension at which the energy is lowered by entry of a lateral p x -type DW. Once one DW can lower the energy, DWs proliferate the channel. Interactions between DWs regulate the DW density, and thereby determine the period of the PDW phase for D ≲ D c 2 .
The stability of the PDW phase very near D c 2 (T ) is analyzed on the basis of a quasiclassical reduction of the Luttinger-Ward (LW) free-energy density for confined phases of 3 He. The resulting LW free energy is then a functional of the quasiclassical Green's function, G, and mean-field pairing self energy, ∆, 64 where the functional ∆Φ[ G] is defined by an digrammatic expansion of the free energy functional in terms of effective interactions of low-energy quasiparticles and Cooper pairs. For our analysis we retain the leading-order weak-coupling BCS corrections to the LW free-energy for the normal Fermi-liquid, c.f. Ref. 68.
The operation, Sp, is short-hand for Tr 4 X(p, r; ε n ) .
The stationary conditions for the LW functional generate both Eilenberger's transport equation for the propagator, G, as well as the self-consistency equation for the off-diagonal selfenergy, i.e. the "gap equation". The auxiliary propagator, G λ , is a function of the variable coupling constant, λ , and is defined as the solution of the Eilenberger equation with a rescaled interaction, or equivalently ∆ λ ≡ λ ∆(p, r). For λ = 1, corresponding to the full pairing interaction, we obtain the self-consistent propagator, G, and self-energy, ∆, whereas λ = 0 corresponds to the normal state with ∆ = 0 and G 0 = −iπsgn(ε n ) τ 3 .
Using the self-consistency condition from variations of the LW functional with respect to the propagator, and using the linearized gap equation to regulate the log-divergent Matsubara sums, the LW functional can be reduced to a functional of the off-diagonal pairing self-energy, where the free-energy density is dλ dS p T ∑ n ∆(p, r)F * λ (p, r; ε n ) + ∆ * (p, r)F λ (p, r; ε n ) − π |ε n | |∆(p, r)| 2 , (60) where f λ (p, r; ε n ) is the solution of the Eilenberger transport equation defined with the re-scaled self-energy, ∆ λ . This functional generalizes the GL free energy functional to all temperatures and to length scalesh/p f ≪ D ≲ ξ 0 .
The difference in free-energy density between the PDW phase and the laterally confined A phase is shown in Fig. 16. The free energy density of the PDW phase is high at the center of domain wall, but the PDW phase gains in condensation energy at the junctions of the DW and the edges. The gain in energy at the junctions dominates as the width of the channel is reduced. Below the critical channel dimension, D c2 , the PDW phase has lower free energy, and is thus the stable ground state. The numerical results for the upper critical channel width, D c2 (T ), are shown in Fig. 10 as the blue squares. The resulting PDW phase is the thermodynamically stable phase over a wide range of channel widths (∆D ≈ 7ξ 0 ≈ 5.6 µm) and temperatures for laterally confined 3 He.

X. SUMMARY AND OUTLOOK
The ground state of thin films of 3 He is the chiral Anderson-Morel state ( 3 He-A), which breaks both time-reversal symmetry and mirror symmetry. The 2D A phase is a topological phase with Chern number N = ±1 for the condensate of each Fermi disk. A consequence of the topology of 3 He-A films is a spectrum of chiral edge states, their signature being a ground-state edge current. Under lateral confinement, edge states on opposing boundaries hybridize and develop a an intricate band structure for both the sub-gap and continuum spectrum. Hybridization of the edge spectrum leads to a reduction of the maxiumn edge current at any temperature. The non-monotomic temperature depdence of the edge current was the first signature of new ground states of superfluid 3 He films under strong confinement. Indeed thin films of 3 He undergo a sequence of phase transitions. At T = 0, starting from the chiral A phase in the limit of laterally unconfined films (D → ∞), superfluid 3 He undergoes a second-order transition to a PDW phase with broken translational symmetry at D c2 ∼ 16ξ 0 . At this scale the PDW phase is a periodic array of chiral domains with alternating chirality, separated by domain walls, the latter also are chiral. The PDW phase breaks time-reversal symmetry, translation invariance, but is invariant under the combination of time-reversal and translation by a one-half period of the PDW. Under stronger confinement a second-order phase transition occurs to a non-chiral "polar phase" at D c1 ∼ 9ξ 0 , in which a single p-wave orbital state of Cooper pairs is aligned along the channel. This remarkable sequence of phase transitions is also mapped out in the full temperature-confinement phase diagram.
The role of disorder on the confining boundaries has not been addressed in terms of the phase diagram for laterally confined 3 He, but in terms of testing the theoretical predictions of new phases of laterally confined 3 He thin films, fortunately surface disorder can be masked by pre-plating the confining surfaces with a few monolayers of superfluid 4 He, which has been demonstrated to lead to specular surface scattering. 44 In this limit thin films of superfluid 3 He-A are expected to persist well into the quasi-2D regimeh/p f ≪ w ≪ ξ 0 , at least until the effects of dimensional confinement modify the pairing interaction.

(A14)
Using cos(2φ ε ) = (2ε 2 /∆ 2 ) − 1 and cos(2φp) =p 2 x −p 2 y , Eq. A14 can be expressed as which is Eq. 18. In the limit D → ∞ the sub-gap bands collapse to the two chiral branches, ε ± = ±∆p y . The branch with ε + = +∆p y is confined on the upper edge x = D, while the branch with dispersion ε − = −∆p y corresponds to the edge state on the lower edge x = 0. For D finite, the sub-gap bandwidth increases for stronger confinement; the exponential dependence of bandwidth originates from the exponential decay of edge states. The momentum dependence of the sub-gap band shows that the bandwidth is largest for normal incidence p x = p f , and is smallest for grazing incidence p x = 0. This agrees with the mapping to a periodic array of domain walls; a trajectory normal to the edges has the smallest period, 2D, while trajectories near grazing incidence have a diverging period.

Appendix B: Solutions for the Quasiclassical Propagator
Here we provide the analysis for the solution of the Eilenberger equation and normalization condition for a piece-wise constant chiral order parameter. We convert the Nambu matrix representation of the Eilenberger equation to coupled differen-tial equations for a vector representation of propagator, 4 For a constant order parameter amplitude, we have ∆ 1 (p) = ∆p x , ∆ 2 (p) = ∆p y . We can express | G ⟩ in terms of the eigenvectors of M, M| µ ⟩ = µ| µ ⟩. The eigenvector with µ = 0 is This solution corresponds to the bulk equilibrium propagator, The eigenvectors corresponding to eigenvalues µ = ±λ are where λ 1 ≡ ∆ 1 (p) 2 − ε 2 . These eigenvectors correspond to Nambu matrices, These are the base matrices defining the propagator in Eq. 21.
In the vector represenation, the solutions for the incident and reflected trajectories are which yields C 1 out = C 1 in ≡ C 1 and C 2 out = C 2 in ≡ C 2 , with The coefficients defining the propagator in Eq. 21 are