Fermionic CFTs and classifying algebras

We study fermionic conformal field theories on surfaces with spin structure in the presence of boundaries, defects, and interfaces. We obtain the relevant crossing relations, taking particular care with parity signs and signs arising from the change of spin structure in different limits. We define fermionic classifying algebras for boundaries, defects, and interfaces, which allow one to read off the elementary boundary conditions, etc. As examples, we define fermionic extensions of Virasoro minimal models and give explicit solutions for the spectrum and bulk structure constants. We show how the A- and D-type fermionic Virasoro minimal models are related by a parity-shift operation which we define in general. We study the boundaries, defects, and interfaces in several examples, in particular in the fermionic Ising model, i.e. the free fermion, in the fermionic tri-critical Ising model, i.e. the first unitary N = 1 superconformal minimal model, and in the supersymmetric Lee-Yang model, of which there are two distinct versions that are related by parity-shift.


JHEP06(2020)025
for example they have been used in the coupling of the Ising model to boundary magnetic fields [5,6], see also [7,8] for more recent applications to boundary renormalisation group flows. Our analysis shows they are a necessary and integral part of the description of fermionic CFTs with boundaries, defects and interfaces. We illustrate these ideas in the concrete cases of fermionic extensions of Virasoro minimal models which we define in section 4, and for which we give an explicit solution for the bulk structure constants. As is well known, for some values of the central charge there are two or more different Virasoro minimal models labelled by pairs of Lie algebras [9]. We show that, remarkably, the fermionic extensions of the (A m , A 4n−1 ) and (A m , D 2n+1 ) bosonic minimal models -and in particular their bulk structure constants [10,11] -are related by the parity-shift operation.
The fermionic extensions of the Virasoro minimal models include such important examples as the fermionic Ising model, i.e. the free fermion, and the fermionic tri-critical Ising model, i.e. the first unitary N = 1 superconformal minimal model, as well as many other theories with extended symmetries. Sections 5 and 6 of the paper are taken up with exploring our results in these situations and comparing our findings with the discussions already in the literature, in particular [12][13][14].

Bulk fields in fermionic CFT
In this section we explain how to describe fields of a fermionic CFT. In order to include the effects of the spin structure and the bosonic/fermionic nature of fields, we introduce a special type of topological line defect; all bulk fields are then connected to one of these defects, so that we think of them as disorder fields which sit at the starting point of the specific topological defect.
We start by describing the relevant properties of the topological defect and then use these to define OPE coefficients and to obtain the crossing constraint they have to satisfy. We show that, given one solution to the crossing constraint, one can obtain another solution by shifting the Ramond sector parity and modifying the given solution by signs.

The topological defect F
The spin structure on the worldsheet of a fermionic CFT is encoded by a topological defect which we call F . The technical details of this procedure are given in [1,15]. Here we do not need the full formalism and just state the properties we will use below.
We denote by H F the space of disorder fields that sit at the start of the topological defect F (figure 1a). Since F is topological, H F carries a representation of the holomorphic and anti-holomorphic copy of the Virasoro algebra.
The space H F is a super-vector space, that is, it is Z 2 -graded into an even and an odd component, We refer to this grading as parity, and for a homogeneous element φ ∈ H F we write |φ| ∈ {0, 1} -or just φ if no confusion can arise -for its parity. There is a second Z 2 - Figure 1. Properties of the topological line defect F . a) A field φ ∈ H F sits at the start of F . b) The OPE of the weight zero defect field π and a field φ multiplies φ by ±1, depending on its parity. c) The weight zero junction joining two T defects into one. d) Pushing the defect field π through the defect junction. e) Associativity relation for the junction field. f) Dragging the F defect through a field φ of spin grade ν φ inserts π ν φ . g) Rotating φ ∈ H ν φ F by 2π can be traded for an insertion of π ν φ +1 . In particular, one cannot just unwind an F -defect around φ, instead the tangent at the insertion point of φ has to remain fixed.
grading on H F whose components are called the Neveu-Schwarz and the Ramond sector, We will refer to this as the spin grading. , so that altogether, H F is Z 2 × Z 2 -graded. For a field φ ∈ H F that is homogeneous with respect to the spin grading we write ν φ ∈ {0, 1} for its degree.

JHEP06(2020)025
We will refer to fields in H F as bulk fields of the fermionic CFT. The defect F has the following properties: 1. There is a parity even weight zero 1 defect field π on F which implements parity on F in the sense that the OPE of π with a bulk field φ is (−1) φ φ (figure 1b). It will be convenient to write π 0 = 1 for the identity defect field and π 1 = π.
2. There is a parity even weight zero defect junction joining two in-coming F -defects into an out-going F -defect (figure 1c). This junction commutes with π in the sense shown in figure 1d, and it is associative as shown in figure 1e.
3. Taking a bulk field φ of spin grade ν φ past an F -defect results in the insertion of π ν φ as shown in figure 1f.
One can verify that the effect of a 2π-rotation of a bulk field φ can be replaced by an insertion of π ν φ +1 on the F -defect (figure 1g), see [1, lemma 4.7]. Thus we get Denote the conformal spin of a field by S φ := h φ −h φ . The above relation implies that φ is a parity-odd NS-field , Z; otherwise . (2.4)

Bulk structure constants
To define the OPE of bulk fields, we need to fix a convention for the spin structure in a neighbourhood of the fields. In the present formalism, this is done by requiring a particular pattern of defect lines. The convention we will use is that for x > y real, (2.5) where the φ i are a basis of primary fields in H F , and ∆ i = h i +h i denotes the scaling dimension of φ i . One could alternatively have chosen the convention that the F -defect starting at φ i passes below φ j . According to figure 1f, the two choices are related by (2.6) If we denote the structure constants computed in the F -defect-passes-below convention by C k ij , then the resulting relation isČ In the following we will stick to the convention in (2.5). 1 We use "weight zero" to mean that it behaves like a vacuum field, i.e. that it is annihilated by the translation operators L−1 andL−1. Figure 2. a) The arrangement of F -defects needed to substitute the OPE in the a → 0 limit. b) The corresponding configuration required in the b → 0 limit.

Crossing symmetry constraint
Consider a correlator f (a, b) of four primary bulk fields φ i , φ j , φ k , φ , inserted at positions a + ib, ib, a, 0 for some a, b > 0, Here it is understood that the F -defects are placed as in figure 2a, and that the ordering relevant for the parity signs is radial ordering. For simplicity, in the derivation of the crossing constraint (and in that of similar constraints below) we restrict ourselves to theories which are rational with respect to the Virasoro symmetry. Extended chiral algebras can be treated in the same way, but one has to account for two complications. Firstly, the fusing matrices will in general carry multiplicity labels. Secondly, the leading contribution in the OPE of two primary fields (with respect to the extended symmetry) may be a descendent field, which makes the definition of the OPE coefficients more involved.
The defect arrangement in figure 2a is such that in the a → 0 limit, we can substitute the bulk OPE right away. Hence in this limit we simply have The phase e −πiSp = e −πi(hp−hp) in (2.9) is determined by our convention on how to continue the OPE (2.5) to the configuration in figure 2 where the result of the OPE of fields φ i with φ j and φ l with φ k results in fields at positions 0 and ib. We choose to continue the OPE so that the defect lines stay in the same topological arrangement and do not cross through the field insertion points which determines the phase uniquely. In the b → 0 limit we have to rearrange the F -defect before we are allowed to substitute the OPE. The relevant configuration is shown in figure 2b. We have (2.10)

JHEP06(2020)025
Here we used the associativity of the junction field and the effect of dragging an F -defect through a field, see figure 1e, 1f. The π-insertion contributes the sign factor (−1) ν j φ k . The resulting graph of F -defects can be brought to the form in figure 2b by using the associativity relation once more. We see that bringing the spin structure to the form required for the OPE contributes a sign factor (−1) ν j φ k . Another sign arises from parity as the order of φ j and φ k changes.
Altogether we get (2.11) The rest of the computation is a standard manipulation of conformal blocks. The overall result iŝ (2.12) Here, the entries "i" and "ī" in the two F-matrices refer to the holomorphic and antiholomorphic conformal weight (h i ,h i ) of φ i , etc.

Symmetry properties of structure constants
The four-point crossing relation determines how the structure constants behave under permutation of indices. DefineĈ We obtain the following two relations by setting φ j = 1 and φ = 1 in (2.12), respectively: (2.14) In particular, theĈ ijk are cyclically symmetric. Using this we can rewrite the second equality above as a relation betweenĈ jki andĈ kji . After relabelling and dividing byĈ 1 kk this implieŝ In step ( * ) we used that e 2πiS i = (−1) φ i (ν i +1) , etc., from (2.3) and the fact that the OPE preserves the parity and spin gradings. For the description of the classifying algebra below it will be important that structure constants involving only spinless fields are symmetric: To see this, note that the only situation in which both expressions in (2.15) produce a sign is φ i = φ j = 1 and ν i = ν j = 0. But then (2.4) imposes S i , S j ∈ Z + 1 2 , which we excluded. JHEP06(2020)025

Parity shift of the Ramond sector
As we have seen in section 2.1, the state space H F of a fermionic CFT carries two Z 2gradings, namely even/odd and NS/R. Accordingly it splits into four direct summands, As before, for a primary bulk field φ i , its parity is |φ i | ∈ {0, 1} (with 0 being even) and its spin grade is ν i ∈ {0, 1} (with 0 being NS). LetĈ k ij be a solution of the bulk crossing relation (2.12).
Given such a solution, we can construct a new fermionic CFT as follows. The new state space H F agrees with the old one, except for a parity shift in the Ramond sector: That is, the old and new gradings are related by The new structure constants are related to the old ones by signs, It is a straightforward computation to see that C k ij again solves the crossing relation (2.12) (and hence also has the symmetry properties stated in section 2.4).
Applying the shift operation twice produces the structure constants (−1) ν i ν jĈ k ij . The extra sign can be absorbed into the normalisation of the fields (e.g. multiply all Ramond fields by i), so that this reproduces the theory one started from. 2 The parity shifted theory may or may not be isomorphic to the unshifted theory. For example, parity shifting the fermionic Ising model produces an equivalent theory, while parity shifting the fermionic minimal model FM (3,8) produces an inequivalent theory, see sections 5.1 and 6.2.

Classifying algebras
In this section we use the bulk structure constants of a fermionic CFT to define several types of fermionic classifying algebras. These are semisimple super-algebras, graded by parity, whose direct summands are in 1-1 correspondence to boundary conditions (section 3.1), defects, or interfaces (section 3.2), depending on the algebra under consideration. In section 3.3 we observe that either one of these classifying algebras in a fermionic CFT is isomorphic to the corresponding algebra in the parity shifted CFT as an ungraded algebra, but typically not as a super-algebra. 2 It is possible to modify (2.20) so that applying the parity shift twice gives back precisely the structure constants one started from. For example, with γa = e −πi(Sa+φa/2) one can set C k ij = (−1) where in the last step we used (2.15).

Conventions for bulk-boundary OPE and boundary OPE
We now consider the theory on the upper half plane with some conformally invariant boundary condition α placed on the real line. As for bulk fields, boundary fields will also serve as starting point for an F -defect. We denote the space of boundary fields by H (3.1) Our convention for the defect arrangement near an insertion of a boundary field ψ ∈ H (α) The bulk-boundary structure constants B φ ψ are defined as, for a primary bulk field where the sum runs over a basis of primary fields ψ in H (α) F . In pictures the above relation looks as follows, (3.4) We will also use the OPE of boundary fields, for which our convention is, for x > y, where the sum runs over a basis of primary fields ψ r in H (α) F . We will only consider boundary conditions for which the bulk-boundary OPE and the boundary OPE preserve parity.

Bulk-boundary crossing relation
Let g(a, b) stand for the correlator of two bulk fields and one boundary field shown in figure 3. The crossing relation is obtained by comparing the a → 0 and b → 0 limit. Figure 3. Two bulk fields and one boundary field together with their defect lines as used in computing the crossing constraint linking bulk-boundary couplings and bulk structure constants. In the correlator g(a, b), the boundary field ψ z inserted at L is moved off to ∞.
It turns out that this relation does not involve parity signs or spin structure signs. The computation is thus standard and we will be very brief. The asymptotic behaviour of g(a, b) in the two limits is The crossing relation now follows from a computation with five-point Virasoro conformal blocks similar to those in [16], and the result is If the boundary field ψ z has weight zero (but could still be of either parity, see next subsection), this relation simplifies to where we used that B k z = 0 requires h k = h k . If the boundary fields ψ x , ψ y , ψ z all have weight zero (but could again be of either parity), this relation simplifies further to The corresponding relation for bosonic theories was derived in [17].

Boundary classifying algebra: bosonic case
Let us first recall the situation when there are just bosonic (i.e. parity even) fields and consider the space of weight zero boundary fields which can arise in the bulk-boundary OPE. Looking at the right hand side of the crossing constraint (3.9), this is symmetric under i ↔ j and as a consequence c z xy must be symmetric under x ↔ y. This means that JHEP06(2020)025 the weight zero boundary fields which arise in the bulk-boundary OPE form a commutative, associative algebra. If one restricts one's attention to situations where this algebra is semisimple, by the Wedderburn theorem it must be equal to a direct sum of copies of C, each of which can be considered as the identity operator on an elementary boundary condition. Due to this observation, it is enough to consider elementary boundary conditions, i.e. boundary conditions on which there is a single boundary field of weight 0 (the identity field). The bosonic version of (3.9) reduces to This equation can be simply seen as a set of equations for the constants B i , but as shown in [3,18] this is also the defining relation for a commutative associative algebra with generators B k , the (bosonic) boundary classifying algebra. From this point of view, the elementary boundary conditions correspond to the one-dimensional representations of this algebra, or, equivalently, to the number of summands C in the direct sum decomposition of (3.10).

Boundary classifying algebra: fermionic case
In the fermionic case, we equally take the coefficientsĈ k ij F k0 j j i i in (3.9) as the structure constants of an algebra with generators B i , where the indices i, j, k now run over all spinless fields of either parity. We denote this algebra byB and call it the fermionic boundary classifying algebra.
In the examples we consider in sections 5 and 6, this is indeed an associative algebra which is in addition semisimple. In [2] we will show that (3.11) defines a semisimple associative algebra for fermionic rational CFTs in general. While theĈ ij k are not necessarily symmetric under i ↔ j, we saw in (2.16) that they are symmetric for all i, j and k which have spin 0 and which can therefore couple to boundary fields of weight 0. It follows that the algebraB is commutative, Since the bulk fields are graded by their parity, we can also viewB as a super-algebra with the parity of the generator B i being the same as that of the field φ i . As the above commutativity relation does not involve parity signs,B is in general not super-commutative. Note thatB is in fact bi-graded, by parity and by spin, since the bulk structure constantŝ C ij k preserve both these gradings.
Let us assume thatB is semisimple. SinceB is commutative (rather than supercommutative), the super-Wedderburn theorem [19, corollary 2.12] states that the algebra splits into a sum of copies of C 1|0 = C and the Clifford algebra Q(2) ≡ C 1 which is twodimensional with one odd generator, a, satisfying a 2 = 1. In other words,B is equal to a sum of m copies of C with generators e α (where α labels the copy) and n copies of C 1 JHEP06(2020)025 with even generator f β and odd generator a β , (3.12) These generators satisfy the relations After analysingB in some detail, let us return to the bulk-boundary crossing constraint (3.9). Analogous to the bosonic case, we conclude that the weight zero boundary fields which arise in the bulk-boundary OPE form an associative super-algebra which is commutative but not necessarily super-commutative, and which is in addition semisimple. 4 Again by the super-Wedderburn theorem, this algebra of boundary fields decomposes into a direct sum of C's and C 1 's. We will call a boundary condition of the fermionic CFT elementary if its weight zero boundary fields consist of either exactly C or exactly C 1 . Note that in the standard basis of C and C 1 , the boundary structure constants c z xy in (3.9) are either 0 or 1.
In the case of C 1 there is a weight zero boundary fermion a which satisfies a 2 = 1. Such weight-zero fermions are well-known and appear in treatments of the boundary Ising model, for example in [5][6][7][8].
We note that as an ungraded algebra, C 1 is isomorphic to C⊕C, but the corresponding generators would not have a fixed parity. Consequently, it is not possible to split a boundary condition with weight zero field content C 1 into two more elementary boundaries without breaking parity-preservation of the bulk-boundary OPE.
We conclude that the indices α in (3.12) label the different elementary boundary conditions: Each pair {f β , a β } corresponds to a boundary condition which supports an odd weight 0 field, while each generator e α corresponds to a boundary condition which does not. The bulk-boundary structure constants are given by the action of the algebra elements B i on the generators, Here, we added the superscripts (α), (β) to distinguish different solutions to (3.9).
Analogous to the bosonic case, there is a relation between elementary boundary conditions and representations of the fermionic boundary classifying algebra. Namely, Both equivalences follow from the description of elementary boundary conditions in terms of direct summands in (3.12). Note, however, that elementary boundary conditions do not JHEP06(2020)025 correspond to irreducible Z 2 -gradedB-modules up to only even isomorphism (this would count summands C with a factor of two and summands C 1 with a factor of one). 5 As can be seen from the structure of the boundary classifying algebra of fermionic theories, restricting to the bosonic generators and removing the fermionic generators reduces each copy of C 1 = C 1|1 to just C. In this way we see that the elementary boundary conditions of a fermionic theory are in 1-1 correspondence with the elementary boundary conditions of its bosonic projection (the related "spin theory"); consideration of the fermionic generators allows one to see which boundary conditions support a fermionic weight zero field and are hence "supersymmetric" as we discuss below.

Bulk-defect OPE and defect classifying algebra
Topological defects in a conformal field theory C have been long known to be equivalent to a particular class of conformal boundary conditions (a special case of so-called permutation boundary conditions [20]) on the doubled model C × C, and so can be studied from that perspective. It can, however, also be useful to consider them in their own right and that is what we do here.
A topological defect will always have a non-trivial set of defect fields, as it will always support the full set of bulk fields; in the case of the trivial identity defect this is exactly the space of defect fields, but in general it will be larger.
It does not add any difficulty to generalise the situation slightly to topological interfaces with a CFT C above the interface and a CFT C below the interface. Below we will sometimes use the term defect to include both situations, with C and C being equal or different.
The sewing constraint for topological defects is almost identical to the bulk crossing symmetry constraint (2.12), with the replacement of some bulk structure constants by bulk-defect structure constants and defect-defect structure constants. We shall here focus only on weight zero fields on the defect as this is sufficient to derive the defect classifying algebra. The bulk-defect structure constants in the OPE of two bulk fields φ i (above the defect) and φ k (below the defect) to a weight zero defect field ϑ z will be denoted by D ik z and the OPE of weight-zero defect fields will be taken to have structure constants c xy z : The result is the defect sewing constraint

JHEP06(2020)025
As in the boundary case, we can consider the couplings to the weight zero fields as generators D ik of an algebraD -the fermionic defect classifying algebra -defined in terms of the constants on the right hand side of (3.17), Here, the pairs ik, j and pq have to satisfy the condition that for each pair xy the bulk fields φ x and φ y have the same left conformal weight and the same right conformal weight. This implies, for example, that F p0 The algebraD is commutative and associative (at least in fermionic rational CFTs [2]). We take the generators D ij to inherit the parity of the product of fields φ i φ j , so thatD is also a super-algebra (not necessarily super-commutative). As in the boundary situation, we assume thatD is semisimple, in which case it decomposes into a direct sum of copies of C and C 1 which are in 1-1 correspondence with the solutions of the defect sewing constraints and hence in correspondence with the defect operators (up to an overall sign) and elementary defect conditions.
In the case C = C , one generic solution for the D ik z is provided by the trivial defect. The only weight zero field on the trivial defect is the identity bulk field 1, and we can set The bulk-defect crossing relation (3.17) then turns into a special case of the bulk crossing relation (2.12).
Interfaces between fermionic and bosonic theories can also be treated by the above classifying algebra by simply choosing one of the two CFTs to be purely even.
In the purely bosonic case, classifying algebras for topological defects were studied in [4].
As opposed to case of boundary conditions, the number of defects of a fermionic theory is larger than the number of defects of its bosonic projection since there is a bosonic generator D φφ of D for each field φ, whether it is bosonic or fermionic. Each defect of the bosonic theory can be associated with one or more defects in the fermionic theory (as it happens, in the examples we have looked at, each bosonic defect is associated to two defects in the fermionic theory).

Classifying algebras in the parity shifted theory
When comparing the various classifying algebras between a fermionic theory and its parityshifted version, one finds that they are isomorphic as ungraded algebras, and that one can give an isomorphism by a simple rescaling of the generators.
However, we stress that the fermionic classifying algebras will in general not be isomorphic as super-algebras (with Z 2 -grading given by parity). We will see this explicitly in the example of FM (3,8) in section 6.2.

Fermionic boundary classifying algebra
Let us write B i for the generators in the unshifted theory as in section 3.1.4 and denote the structure constants of the fermionic boundary classifying algebra by β k ij , such that (3.11) becomes B i B j = k β k ij B k . The bulk structure constants in the parity-shifted theory JHEP06(2020)025 are given in (2.20). If we denote the generators and structure constants of the classifying algebra of the parity shifted theory by B i and β k ij , we get In general, the factor (−1) ν i φ j cannot be absorbed into a rescaling, 6 but since certain of the β k ij are zero, here this will be possible. Explicitly, we may identify The scalar coefficients λ are not unique, and we just exhibit one possible solution. To verify the above equality, first note that λ i λ j /λ k = (−1) φ i φ j . One thus needs to check that Since all NS-sector generators in the boundary classifying algebra are parity-even, we only need to consider the case that B i and B j are from the R-sector, and that B i is even and B j is odd. But then β k ij = 0, as B k would then need to be odd and in the NS-sector.

Fermionic defect classifying algebra
Let us use the letter "U " for the unshifted theory and "P " for the parity shifted theory. Here we will give the relation between the defect classifying algebras of types U -U , U -P , P -U and P -P , where the first letter refers to the theory in the upper half plane. The U -U -case is as given in section 3.2 and will be our reference case: with β (U U ) pq ik,j determined by (3.18). A computation similar to the boundary classifying algebra gives (all parities are stated with respect to the unshifted theory) To verify these equations one needs to use that for a generator D xy one always has ν x = ν y (but not necessarily φ x = φ y ) as otherwise the bulk-defect OPE does not contain 1 or a.

Fermionic Virasoro minimal models
The examples we shall consider are all Virasoro minimal, i.e. the Hilbert space is formed from a finite set of Virasoro representations. This might seem odd, as it might be more natural to start with a super-algebra for which there is naturally a fermionic interpretation, but as we shall see it provides ample examples, not only infinite series which extend the bosonic Virasoro minimal models but also the simplest example, a free fermion, as well as examples with super-Virasoro symmetry and extended (W-algebra type) fermionic algebras.
In this section we give the bulk structure constants for fermionic A-and D-type minimal models, in section 5 we treat the free fermion in detail, and in section 6 we give further Virasoro examples.
Consider the Virasoro algebra at the minimal model central charge c(p, q) = 1 − 6(p − q) 2 /(pq) with p, q coprime integers greater than 1. We recall that the possible Virasoro representations are labelled by two integers (r, s) where 1 ≤ r < p and 1 ≤ s < q with the identification (r, s) (p − r, q − s); for more details see [21].
The fermionic generator G we consider has Kac-labels (1, q − 1) ∼ (p − 1, 1) and is the unique non-trivial simple current at that central charge. Note that for G to indeed be different from the vacuum representation we actually need p, q > 2. The conformal weight of G is h G = 1 4 (p − 2)(q − 2). The condition that h G ∈ Z + 1 2 amounts to p = 2n + 1 , q = 4k , or equivalently p = 4k , q = 2n + 1 , . The first few values of h G are realised in the models in table 1. The next value of h G realised in a unitary model is h G = 15 2 at central charge c(7, 8) = 25 28 . Recall that the possible different bosonic field theories with this central charge are labelled by a pair of simply-laced Lie algebras with Coxeter numbers p, q [9]. We will give two explicit solutions to the fermionic crossing relation (2.12). One is an extension of the minimal model of type M (A p−1 , A q−1 ) which we will call the fermionic Virasoro minimal model FM (A p−1 , A q−1 ) or FM (p, q) for short. The other is an extension of the minimal model M (A p−1 , D q/2+1 ) which we will call FM (A p−1 , D q/2+1 ) or FM (p, q) for short. By an extension we mean that the fermionic model contains the full bosonic field theory as a sub-theory. The two models FM (p, q) and FM (p, q) are obtained from each other by shifting the Ramond sector parity as in section 2.5.

A-type fermionic models
Let I be an indexing set for the Kac-table modulo its Z 2 -identification (r, s) ∼ (p−r, q −s), and let M a , a ∈ I be the corresponding irreducible Virasoro representation. In the setting of [1, section 5], the defect F relevant for A-type fermionic models is the superposition of the identity defect and a parity shifted version of the topological defect labelled by the representation G = M (1,q−1) . (What we call F here is called A in [1].) The splitting into NS-and R-sector depends on the sign of the ratio of S-matrices S G,a /S G,0 ∈ {±1}, with +1 being the NS-sector and −1 the R-sector. Explicitly, for a = (r, s) we have S G,a /S G,0 = (−1) qr+ps+1 so that the set I splits as The state space H F of FM (A p−1 , A q−1 ) splits into the four sectors in (2.17) as follows: where Ga is the result of the fusion product of G and a. That is, for a = (r, s) we have contains a holomorphic field of weight (h G , 0) and an anti-holomorphic one of weight (0, h G ). We will use the notation φ e a for the primary field in H ev F of conformal weights (h a , h a ), and φ o a for the primary field in H odd F of conformal weights (h a , h Ga ). This notation is slightly asymmetric in that φ o G is the field of weight (h G , 0) while φ o 1 is that of weight (0, h G ). Suppose we take q = 4k in (4.1). Then Ga = a if and only if a = (r, q 2 ) = (r, 2k), so that for Kac labels of this form we have a ∈ I R and fields with conformal weights (h a , h a ) occur with multiplicity two, once as φ e a ∈ H R,ev F and once as φ o a ∈ H R,odd F . For the structure constants we use the notation where α, β ∈ {e, o} and α + β stands for the parity of the product. One solution to the bulk crossing relation (2.12) is given bŷ The λ e/o i ∈ C × are normalisation constants which can be chosen at will. The proof that these constants indeed solve the crossing constraint will be given in [2]. There it will also be shown that one can find a topological defect F with the required properties.
One standard normalisation is to make a choice of square roots With this choice one hasĈ
Bosonic subtheory. The bosonic subtheory of FM (p, q) is the restriction to the parity even subspace of H F . This can alternatively be understood as the result of summing over spin structures. Explicitly we have which agrees with the state space of the A-type (bosonic) minimal model M (A p−1 , A q−1 ).
Restricting the structure constants to the even subsector and choosing the normalisation λ e a = S 0a /S 00 , the solution (4.5) precisely recovers the bulk structure constants for the A-type models in terms of F -matrices as given in [10], Note that with these conventions, the coupling to the identity field is not normalised to 1.

D-type fermionic models
We now consider the theory FM (p, q) obtained from FM (p, q) by shifting parity in the Ramond sector as in section 2.5. We will assume that p is odd and q is even, which according to (4.1) implies that q ∈ 4Z. We will also denote FM (p, q) as FM (A p−1 , D q/2+1 ) which will be justified later by the restriction to the bosonic subtheory. As opposed to F M (p, q), shifting the antiholomorphic label by G or not no longer corresponds to the field being odd. We therefore label the fields as, for a ∈ I, where u stands for "unshifted" and s for "shifted" (referring to the shift by G). The spin grading ofφ u/s a is ν a as before, but the parity has changed, |φ u a | = ν a , |φ s a | = ν a + 1 . As in the A-type case, for a = (r, q 2 ) we have a ∈ I R , Ga = a, and bothφ u a ∈ H R,odd F and φ s a ∈ H R,ev F have conformal weights (h a , h a ), so that these weights occur with multiplicity two.

JHEP06(2020)025
The structure constants of FM (p, q) are obtained from (4.5) by the transformation (2.20): (4.12) Bosonic subtheory. The even subspace of H F is which is the state space of the D-type minimal model M (A p−1 , D q/2+1 ). Note that since q ∈ 4Z, this model is always of D odd -type (a permutation modular invariant). When restricting the structure constants to the even subsector, i.e. to φ u a , a ∈ I NS and φ s a , a ∈ I R , the expression (4.12) reproduces the structure constants of the D-series bosonic Virasoro minimal model found in [11].
Since the comparison requires a bit of calculation, we give some details. The expression in [11] is (4.14) The notation is as follows. The indices α, β, γ take values in {e, o}. The index i e corresponds to the field φ u i and i o to φ s i . Furthermore, ω e = (1, 1) and ω o = G. For i e we haveh i = h i and for i o we haveh i = h Gi . The label r is given by r = j if α = e and r = Gj if α=o. The constants (ω) B ωα iα , etc., are certain bulk-boundary structure constants which will be absorbed into the normalisation coefficients (together with a phase for the G-shifted fields), The label µ u is the Kac-label (1, q 2 ), which is a fixed point for G. The F -matrix entry F µu ωγ µu µu ωα ω β is either 0 or 1, depending on whether the Z 2 fusion rules in the OPE are obeyed. This grading rule holds by construction for (4.12), so that this F -matrix coefficient can be dropped. One can now check sector by sector that the restriction of (4.12) to the even subsector agrees with (4.14). In sector "su", i.e. for α = o, β = e, this requires an F -matrix identity: Finally, note that since we restrict to the even subsector, for φ u a we have ν a = 0 while for φ s a we have ν a = 1. Accordingly, the factor (−1) νa is equal to 1 in C (su) c ab and equal to −1 in C (ss) c ab . In the latter case, the minus sign cancels against e −2πih G . Finally, let us note that fermionic minimal model FM (p, q) and its parity shifted cousin FM (p, q) are non-isomorphic (as graded theories) whenever the corresponding D odddiagram differs from an A-diagram. Since q ∈ 4Z, this happens for q ≥ 8. On the other hand, the fermionic Ising model FM (3,4), which we will treat in detail in the next section, is isomorphic to FM (3, 4).

The Ising model and the free fermion
The Ising model is the Virasoro minimal model M (3, 4) and the relevant data is given in appendix A.1. The fermionic Ising model FM (3,4) is really the theory of the free fermion -the field φ o = ψ is a free holomorphic fermion and φ o 1 =ψ is a free anti-holomorphic fermion. Altogether, the Virasoro primary fields in the fermionic model are given in table 2.

Bulk structure constants
The structure constants of the free fermion can be read off from (4.5), but in this simple model it is easy to compute them directly from the crossing constraint (2.12), and this is what we will do.
For all fields we set C 1 φφ = 1, which fixes all normalisations up to signs. The structure constants of the even fields are those of the Ising model and are of course well-known: up to the symmetry properties (2.14) the only remaining structure constant is which also fixes the sign-freedom in the normalisation of . The structure constants involving even and odd fields in the fermionic model have been considered before and some are given in [6] and [21, section 12.3.3], but these are only partial results and it is not clear how the various signs were chosen nor how the full consistency could be checked. Using our formalism, we state and solve the sewing constraints, and taking the same normalisation of the primaries as [21] we agree with the partial results stated there. The odd fields in the fermionic model are ψ,ψ and µ. Setting i = j = ψ, k = l =ψ in bulk crossing relation (2.12) results in (Ĉ ψψ ) 2 = −Ĉ 1 ψψĈ 1 ψψ = −1. We link the remaining JHEP06(2020)025 free signs in the normalisation of ψ andψ by setting = iψψ, or, in other words, The symmetry properties (2.14) relate any two permutations of the three primaries in a structure constants. Up to such permutations, the remaining structure constants involving two odd fields (and no identity field) arê To fix their value, we consider the crossing constraint (2.12) for the following four choices of parameters: a) i = ψ, k =ψ, j = = σ and q = ; b) the same with j = = µ; c) i = j = σ, k = l = µ, q = ψ; and finally d) the same with q =ψ, Using (2.14) we get the following relations between the structure constants appearing in the above constraints and those listed in (5.3), Consider for a moment the parity shifted theory FM (3, 4). Its structure constants are related to the ones given above as in (2.20). This turns out to be an equivalent theory, and a choice of parity-grading preserving isomorphism from FM (3, 4) (whose fields we denote by σ, ψ, etc.) to FM (3, 4) is NS-sector: R-sector: This can be interpreted as the action of an invertible interface and agrees with the solution to the (ungraded) fermionic defect classifying algebra called "duality 2 " in table 8 below.
JHEP06(2020)025 Table 3. Solutions to the bulk-boundary sewing constraints for the fermionic Ising model.

Boundary classifying algebra and boundary conditions
In the fermionic Ising model case, only the three spinless fields ( , σ, µ) can couple to a weight zero field on the boundary and of these, the two even fields ( , σ) can only couple to a bosonic field of weight zero while the odd field (µ) can only couple to a fermionic field of weight zero. We shall denote a possible fermionic weight zero field by a and normalise its two point function to aa = 1, so that c 1 aa = 1. As a consequence, we would like to find the sewing constraints on the following set of bulk-boundary structure constants: and we have six sewing constraints of the form (3.9), namely These equations have four solutions, given in table 3, but the last two can be identified by the change of normalisation a → −a, and so there are three inequivalent solutions to the boundary sewing constraints, two fixed boundary conditions which do not require a fermionic weight zero field, and one free boundary condition which does require a fermionic weight zero field. From (5.10) we can also read off the fermionic boundary classifying algebraB in (3.11) for the free fermion:

JHEP06(2020)025
This is a super-algebra and the parity of each generator agrees with that of the corresponding bulk field: B , B σ are even and B µ is odd. Now we can more easily identify the elementary boundary conditions by expressingB as a direct sum of copies of C and the Clifford algebra C 1 which correspond directly to the elementary boundary conditions without overcounting the free boundary condition (cf. the discussion in section 3. 1.4), where the indices +, −, f correspond to fixed up, fixed down and free. Explicitly, One can check that these satisfy the algebra (3.13) and from the action of the generators B , B σ and B µ , we can read off the representations: where the 2 × 2 matrices 1 and A are It is easy to see that the matrix representation of the free boundary condition is reducible and (by diagonalising A) splits into exactly the two solution in table 3 but this does not respect the even-odd nature of the generators and leads to the over-counting of the free boundary condition in table 3; only when ensuring that the representation ofB respects the even-odd grading do we get the correct counting (cf. (3.15) for the precise statement).

Boundary field content and boundary states
It is also instructive to consider the contribution of non-zero weight boundary fields as this will determine the full boundary field content of each boundary condition as well as give the gluing conditions for the fermion fields on the boundaries.

Boundary field content
Given the bulk field content in the fermionic Ising model, one only has to consider h = 1 2 boundary fields as no other couplings are possible. The first result is that there must always exist a fermionic weight 1 2 field, which we denote µ B and which could be identified with the boundary disorder field. The reason is that the bulk fermionic fields can only couple to a weight 1 2 field on the boundary. If we JHEP06(2020)025 denote the bulk-boundary coupling of ψ as B ψ µ B and if the boundary field two point function normalisation is c µ B µ B 1 , then the sewing constraint (3.7) requires (B ψ µ B ) 2 c µ B µ B 1 = 1 and hence we see that not only must such a field exist but that the bulk fermions must couple to it. We shall not consider the possibility of two distinct boundary fermions but instead suppose that both bulk fields ψ andψ couple to the same boundary field. The next question is whether there is also an even (bosonic) weight 1 2 boundary field σ B (which could be identified with the boundary spin). If there is such a field then its OPE with µ B must necessarily be to an odd weight 0 field and conversely the OPE of an odd weight zero field with µ B must be an even weight 1 2 field σ B . If we normalise a 2 = c σ B σ B 1 = c µ B µ B 1 = 1 then the OPE algebra is forced by the boundary sewing constraints to have the form where λ 4 = 1. The values λ = ±1 turn out to be inconsistent with the bulk-boundary sewing constraints and λ = ±i are equivalent under a field redefinition, hence we make the choice λ = i from here on. With this boundary field algebra and the seven structure constants   solutions, there are no couplings to the fields a and σ B and so these fields can be consistently excluded from the set of boundary fields, as expected.
Only three of the eight solutions are physically distinct. Namely, for each of the sign choices ζ, ξ ∈ {±1} we can redefine the boundary fields as a → ζ a, σ B → ξ σ B , µ B → ζξ µ B . This agrees with the three physically distinct solutions found in table 3.
What is perhaps surprising is that the full set of sewing constraints in table 4 also defines a commutative algebra [with identity] with generators {1, with the boundary classifying algebra as a sub-algebra, and the solutions in table 5 are the eight one-dimensional representations of this algebra.
As with the classifying algebra, this commutative algebra can also be viewed as a superalgebra with even generators {1, }, and odd generators {B µ a , B µ µ B , B ψ µ B , Bψ µ B }, and the eight one-dimensional representations in table 5 combine into four representations of this super-algebra on C 1|1 , of which the final two are again related by σ B → −σ B , µ B → −µ B and so are equivalent physically; these are given in table 6. This table also shows that the fermions have opposite gluing conditions on fixed and free boundary conditions. On the free boundary condition with B = −1, the bulk fermions obey ψ = −ψ on the boundary; on the fixed boundary conditions they obey ψ =ψ.

Boundary states
Since we are working with theories on a fixed spin structure, each boundary condition b will be associated to two boundary states, one state |b NS entirely in the NS sector and one state |b R entirely in the R sector. In our approach, correlators are overall even linear forms, so that a one-point correlator on a disc is an even linear map H F → C. This implies that boundary states are purely even.
The bulk-boundary couplings to the identity field (cf. table 3) determine the overlap of the boundary states with each primary bulk field, normalised by the overlap with the vacuum. Namely, for primary bulk fields φ i in the NS sector and φ j in the R sector we have This fixes the expansion of |b NS,R in terms of Virasoro-Ishibashi states |1 , |σ , | up to the overall constant 0|b N S . To determine 0|b NS (up to a sign) we will make use of the fact that we know the boundary field content for each boundary condition: This means that the NS boundary states have to satisfy where as usual q = exp(2πiτ ),q = exp(−2πi/τ ) and χ h are the characters of the Virasoro highest weight representations of weight h. Putting all this together, we arrive at (up to an overall undetermined sign for each boundary condition) The overlaps between the R-sector boundary states must give the supertrace over the field contents, and the above states correctly give these supertraces, As briefly mentioned in the end of section 3.1.4, one can think of boundary states as coming in pairs which differ by an overall sign; the three boundary states presented above are then characterised by the condition 0|b NS > 0.
Note that the NS components in (5.21) agree with the conjectures in [12] for the NS-sector of the boundary states, although that paper does not correctly account for the fermionic fields in this model nor the Ramond sectors. Boundary states in the NS-sector of free fermions are also considered in [22], and the boundary states for a single Majorana fermion given in appendix D there agree with the NS-component in (5.21) up to normalisation.

JHEP06(2020)025
Finally, we note that our boundary states are related to the boundary states for the purely bosonic Ising model by taking a superposition of the NS-and R-sector boundary states,

Classifying algebra and defect conditions
We consider now the defect classifying algebra of the standard Ising model and the fermionic version, assuming that there is a single bosonic field of weight zero (denoted 1) on the defect and at most one fermionic defect field of weight zero denoted a which satisfies a 2 = 1. The bulk-defect structure constants for the non-identity bulk fields in the Ising case are {D 1 , D σσ 1 }, and for the fermionic case they are {D 1 , D σσ 1 , D µµ 1 , D σµ a , D µσ a , D ψψ 1 , Dψψ 1 }. The sewing constraints are given by the multiplication table 7.
Viewed as equations, there are eight solutions to the defect sewing constraints, given in table 8.
We can view the bulk-defect structure constants as algebra generators of a super-algebraD, the fermionic defect classifying algebra, with odd generators {D σµ a , D µσ a } and the rest being even. As an ungraded algebraD is commutative, and the solutions in table 8 are its eight one-dimensional representations. When viewed as representations of the superalgebra, the first four are one-dimensional representations on C 1|0 while the second four combine to form two two-dimensional representations on C 1|1 . Accordingly,D decomposes into four copies of C (with generators e α ) and two copies of C 1 (with generators {f α , a α }) JHEP06(2020)025 Table 8. Solutions for the bulk-defect structure constants in the fermionic Ising model.
satisfying the same algebra (3.13). The explicit forms are:

Defect field content and defect operators
Topological defects can also be described by operators on the Hilbert space of the theory which commute with the Virasoro algebra, and hence are sums of intertwiners between equivalent Virasoro representations. The bulk-defect structure constants are proportional to the coefficients of these intertwiners and hence the defect operator is defined, up to an overall scalar multiple, by the bulk-defect structure constants. This overall multiple can itself be fixed (up to a sign) by the requirement that it correctly determines the field content on the defect, and in particular that it correctly determines the dimension of the space of zero-weight fields on the defect. We illustrate this in the case of the bosonic and fermionic Ising defects. The bosonic Hilbert space is

JHEP06(2020)025
where M a denotes the corresponding irreducible Virasoro representation. This means a topological defect operator takes the form D = αP 00 + βP + γP σσ . (5.26) The constants D ij 1 are given by 27) and hence the defect is fixed up to a scalar, The field content on the defect is then given by the modular transform of the trace on the cylinder, 29) The upshot is that λ = ±1 for the identity and spin defects and ± √ 2 for the duality defect, cf. [23,24].
The same ideas can be applied to the defects in the fermionic model, with the observation that the defects come in two versions, depending on the spin structure, so that there are separate defect operators acting on the Neveu-Schwarz and Ramond sectors of the Hilbert space. In the fermionic model the Hilbert space is The defect operators in the two sectors are then given by the bulk-defect structure constants up to an overall constant The main difference is that the space of weight zero fields on the defects can now be either one or two, depending on whether the defect supports the fermionic weight zero field a or not. When we perform the calculations we find thatD NS are in fact identical with the defect operators proposed in [12]. We give the explicit forms in table 9.
We note that the defects of the bosonic Ising model are given by the combinations Table 9. Explicit forms of the defect operators in the fermionic Ising model. Table 10. Interface sewing constraints for an interface between the standard and fermionic Ising models.

Interfaces between Ising and fermionic Ising
As an example, consider interfaces between the bosonic Ising model and the fermionic version. Let us assume that the Ising model is in the upper half plane and the fermionic model in the lower half plane. This means that the possible bulk-defect structure constants are D ij x where i takes values in { , σ}, j takes values in { , σ, µ} and x ∈ {1, a}. The interface sewing constraints are exactly the appropriate subset of the fermionic defect sewing constraints in table 7, as given in table 10. There are four solutions to the sewing constraints, which constitute the four one-dimensional representations of the commutative interface classifying algebra and which also form two one-dimensional representations on C 1|0 and one two-dimensional representation on C 1|1 of the interface classifying super-algebra, as in table 11.

Further Virasoro examples
In this section we will consider a few examples of fermionic minimal models by increasing weight of the generator G as listed in table 1. The smallest value is h G = 1/2 in the single model FM (4, 3) = FM (4, 3) which is the free fermion and was already treated in section 5.
The next value is h G = 3/2, so that the fermionic theory has super-Virasoro symmetry.  Table 11. Solutions to the interface sewing constraints for an interface between the standard and fermionic Ising models.

JHEP06(2020)025
The N = 1 superconformal minimal model values of c are c(p, q) = 15/2 − 3p/q − 3q/p, parametrised by two integers p, q with p − q even and p, q ≥ 2 [25,26]; we shall denote them generically by SM (p, q). The possible modular invariant partition functions for the unitary models |p − q| = 2 have been classified in [27] but, as noted there, modular invariance does not fix the partition function uniquely and it is determined only up to a constant which is fixed by the parities of the states with h = c/24.
The final value we consider in any detail is h G = 5/2 in FM (4, 7) in section 6.3. This final case was first noted in [28] and a longer list is given in [29] where chiral algebras which extend the Virasoro algebra by a single fermionic field are considered. FM (4, 7) is a reduction the WB (0, 2) algebra (a.k.a. the fermionic WB 2 algebra) at a value of c at which the spin 4 field decouples.
These examples in fact all fit into two infinite series of fermionic models with a current of spin (2k − 1)/2, and these are FM (4, 2k + 1) and FM (3, 4k). For k ≥ 3, these two series are all special cases of the fermionic W-algebra WB (0, k − 1) (a.k.a. the fermionic WB k−1 algebra first introduced in [30]) in which all but the fermionic W-algebra field decouple.

Fermionic tri-critical Ising model
The fermionic tri-critical Ising model is the second in the series of fermionic extensions of the minimal models and some essential data is given in appendix A.2.

Fermionic TCIM boundary classifying algebra
The boundary conditions of the tri-critical Ising model in both bosonic and fermionic models have been studied before. For the bosonic case see [31], and for the fermionic case see for example [13,14]. Here we will study boundary conditions of the fermionic model via the fermionic classifying algebra.
In the fermionic tri-critical Ising model, there are eight spinless bulk fields which can hence couple to a weight zero boundary field, as in table 12.
As before, we assume that the only possible weight zero boundary fields are the identity 1 and an odd field a, and so the bulk-boundary couplings fields are The sewing constraints can again be considered as the relations in an 8 dimensional commutative algebra which also has the form of a non-supercommutative super-algebra with even generators { 1, B , B    algebra with respect to the spin grading, so it is in fact a bi-graded algebra. The generators are shown along with their gradings φ andφ = φ + ν in table 13. There are 28 sewing constraints for these structure constants which have eight solutions corresponding to the eight one-dimensional representations of the commutative algebra with generators (6.1).
When viewed as super-algebra, the boundary classifying algebra decomposes into four copies of C and two copies of C 1 , giving six boundary conditions in all. This means the eight one-dimensional representations combine to form 4 one-dimensional representations on C 1|0 and two two-dimensional representations on C 1|1 , as shown in table 14.
The boundary states in the NS and R sectors are given by taking linear combinations determined by table 14 with an overall normalisation, The normalisation is fixed (up to a sign) by the requirement that the number of weight-zero fields is 1 on the boundaries corresponding to one-dimensional representations and 2 on the boundaries corresponding to two-dimensional representations.
6.1.2 Comparison to [13,14] There has been a considerable amount of work on boundary conditions of superconformal field theories. Here we compare our results with those of Nepomechie in [14] which discusses boundary states in the tri-critical Ising model M (4, 5) and which of these are "supersymmetric". The boundary states discussed in that paper include both NS-and R-sectors and a boundary condition is said to be supersymmetric if the partition function JHEP06(2020)025 on a cylinder is a sum of characters of the super Virasoro algebra, rather than simply a sum of characters of the Virasoro algebra.
Since the fermionic theory FM (4, 5) we consider includes the generators of the super-Virasoro algebra, the partition function calculated using the overlaps of the boundary states in the NS-sector will always be sums of characters of the super-Virasoro algebra. However, these are not the partition functions discussed in [13] which are instead the average over the two spin structures, that is the average of the overlaps between the NS-sectors and the R-sectors. Since the R-sectors will contribute the supertrace over a super-Virasoro representation, not a trace, any contribution from the R-sector will stop the partition function being a sum of super-Virasoro characters. Hence we see that "supersymmetric" boundary conditions in the sense of [14] correspond to boundary conditions with zero Rsector boundary state, that is boundary conditions for which the algebra of weight zero boundary fields is C 1 ; if the weight zero boundary fields are simply C then the boundary condition is "non-supersymmetric". The latter boundary conditions fall into pairs related by the spin symmetry and which give "supersymmetric" boundary conditions when taken as a superposition.
In this sense we find, as in [14], two "supersymmetric" boundary conditions and four "non-supersymmetric" ones.

Fermionic TCIM defect classifying algebra
The defect classifying algebra of FM (4, 5) has one generator D ij for each pair of fields φ i , φ j whose operator product on the defect includes a field of weight zero. There are 6 such pairs where i and j are both even which generate the defect classifying algebra of the bosonic M (4, 5). There are 6 more pairs for which i and j are both odd, but for which D ij is again therefore even. Together these generate the 12 dimensional even subalgebra ofD. There are further 4 pairs where one of i and j is even and the other odd and these span the four dimensional odd part ofD, see table 15. As a vector space,D is equivalent to C 12|4 . As a super-algebra,D splits into 8 copies of C and four copies of C 1 giving 12 defects in total, twice as many as the purely bosonic defects of M (4, 5), with the differences occurring in sectors that are not in M (4, 5).

JHEP06(2020)025
Even generators D 11 Odd generators D σµ a , D µσ a , D σ µ a , D µ σ a   If we restrict i and j to the NS sector alone, the algebra is instead equivalent to C 8|0 , i.e. it is purely even and the corresponding eight solutions are formally the same as those found in [12].
Up to now, the question has been entirely one of calculating the bulk-defect structure constants. The final step is to fix the normalisation of the defect operators (up to a sign) to reproduce the correct counting of weight zero fields from the torus expectation value of the defect operator which is easily done.

Fermionic TCIM interface classifying algebra
The only interface we can consider is that between M (4, 5) and F M (4,5). Since the fields in M (4, 5) form a subalgebra of the fields in F M (4, 5), the interface classifying algebra is a subalgebra of the defect classifying algebra of F M (4,5). From table 15, there are 6 surviving even generators and 2 surviving odd generators. This means the interface classifying algebra is C 6|2 as a super-vector space, splitting into four copies of C and two copies of C 1 as a graded algebra, giving 6 interfaces in total. The representations and their weights are given in table 16. As usual the Kac-table includes two copies of each representation and in this case we remove the degeneracy by considering only the representations of type (1, s).
The value c = −21/4 corresponds to SM (2, 8) and has been looked at before, see e.g. [32][33][34]. It has been identified as the supersymmetric Lee-Yang model. The issue of different models at the same central charge does not seem to have been considered in these works -they assume that the bosonic projection is the diagonal invariant M (3,8), and hence in our language identify the superconformal theory as FM (3,8).    The representation content in the even and odd Ramond sector of FM (3, 8) and FM (3,8) now differs, and so these models cannot be graded-isomorphic. From the point of view of the super-Virasoro algebra, the parity of the Ramond sector ground state has changed between these two models, cf. table 17, where the R-ground state has label (1, 2).
The structure of the boundary classifying algebras also differs in the two models. According to the general theory in section 3.3, there are equal numbers of bulk-boundary structure constants in the two models, but the parities are different. In FM (3, 8), where the Ramond ground state is even, there are 7 even and 1 odd generators. This leads to a classifying algebra which is C 7|1 as a super-vector space and which decomposes into 6 copies JHEP06(2020)025 of C and 1 copy of C 1 giving seven fundamental boundary conditions. From the point of view of the superconformal algebra, this means there is 1 "supersymmetric" boundary condition in the sense of [13] and 6 "non-supersymmetric" ones. In FM (3,8), where the Ramond ground state is odd, there are 5 even and 3 odd generators, leading to a classifying algebra which is C 5|3 as a super-vector space and which decomposes into 2 copies of C and 3 copies of C 1 giving five fundamental boundary conditions, 3 of which are "supersymmetric" and 2 are not.
As expected, the boundary conditions of FM (3,8) are in 1-1 correspondence with the 7 boundary conditions of M (A 2 , A 7 ), and those of FM (3,8) are in 1-1 correspondence with the 5 boundary conditions of M (A 2 , D 5 ) (see [35] for the boundary conditions of minimal models).
The two boundary classifying algebras can be made identical [as ungraded algebras] by a suitable rescaling of the fields in FM (3, 8) (or equivalently of the generators of the algebra). The generators and their gradings are shown in table 19 together with the rescaling of the generators that makes the algebras identical.
There are now three different fermionic defect classifying algebras, classifying the defects in FM (3,8), defects in FM (3,8) and interfaces between FM (3, 8) and FM (3,8). Again, according to the general theory in section 3.3, these can be made equal as ungraded algebras by a simple rescaling of the generators by phases, but they are not all equivalent as graded algebras. The two defect classifying algebras have the same gradings, but the interface algebra has a different grading. These are shown in table 20. The two defect classifying algebras are C 14|2 as vector spaces, decomposing into 12 copies of C and two copies of C 1 as graded algebras, giving 14 defects in each of these two theories. The interface algebra is C 10|6 as a vector space and 4 copies C and 6 copies of C 1 as an algebra giving 10 interfaces between the two theories. When we descend to the bosonic A and D theories, the defects/interfaces in the fermionic theories are identified in pairs (as in subsection 6.1.3, the pairs only differ in sectors that are not in the bosonic theories) and give 7 defects in the each of the A and D theories and 5 interfaces between the A and D theories, in agreement with the computation in terms of traces of products of bulk modular invariant matrices, see [24] and [36, remark 5.19].

6.3
The FM (4, 7) = WB(0, 2) 5,7 example Finally, the only fermionic minimal model to have a weight 5/2 current is FM (4, 7) with central charge c = −13/14. This turns out to be a restriction of the "fermionic" Walgebra WB (0, 2). This algebra, also known as the "fermionic" WB 2 algebra, extends JHEP06(2020)025   (3,8). The parity of a given generator in either of the two defect algebras is called φ, and in the interface algebra it is calledφ, with the identifications used in (3.22). the Virasoro algebra by primary fields of weights 4 and 5/2. Its structure constants were worked out explicitly in [37] and it can be seen that the field of weight 4 decouples from the algebra at this value of c. The minimal model of this algebra has 6 representations labelled by [rs; r s ] where r, s, r , s ≥ 1, 2r + s ≤ 3 and 2r + s ≤ 5. Three of these are irreducible as Virasoro algebra representations and in the Ramond sector and give the three R sector representations of FM (4, 7); three are reducible, splitting each into two Virasoro algebra representations and these together are the six NS sector representations of FM (4, 7). The rest of the analysis is straightforward -there are 9 boundary conditions of FM (4, 7) of which three are invariant under the W-algebra automorphism W → −W and three are related in pairs. This is exactly as in the fermionic TCIM, where there are two "supersymmetric" boundary conditions and two pairs which are related by G → −G.

Conclusions
We have defined fermionic conformal field theories and their classifying algebras, defined fermionic extensions of the Virasoro minimal models, found explicit expressions for all the bulk structure constants (of both Neveu-Schwarz and Ramond fields) of these models, and have given numerous examples in these cases.
We have found that it is natural for certain boundary conditions and defects to support a weight zero fermionic field which has up to now been introduced in an ad hoc manner.
We have also found that there is a natural parity-shift operation which can relate different theories. On the one hand, this relates bosonic theories as the projections of parity-shifted fermionic theories, and on the other hand this means that there are hitherto unconsidered fermionic theories to be looked at.
There are quite a few questions that are unresolved and which suggest new lines of enquiry.
Firstly, we showed that the full set of bulk-boundary structure constants in the Ising model defined an algebra, not just the couplings to boundary fields of weight zero. This seems to merit further investigation, even in the purely bosonic case.

JHEP06(2020)025
Secondly, we showed that the fermionic extensions of the A-and D-invariant minimal models are related by a parity-shift operation. This leads one to wonder if one can define extensions in other cases in which the extended algebra would be bosonic, such as the A-and D-invariants of M (5, 6), namely the tetra-critical Ising model and the 3-state Potts model respectively. It would also be interesting to investigate possible fermionic extensions of the exceptional invariants; one of these, M (A 4 , E 6 ), is a product theory M (A 1 , A 4 ) × M (A 2 , A 3 ) = M (2, 5) × M (3, 4) [38], with the obvious proposal that F M (A 4 , E 6 ) = M (2, 5) × F M (3,4). It would be good to have a general understanding of these models.
It would also be good to give the structure constants for the field theories on the boundaries, defects and interfaces that we have found which would be a necessary first step to discuss their perturbations and the resulting renormalisation group flows.
In [39], the moduli space of c = 1 CFTs with fermions was investigated (see figure 2 there), it would be interesting to look at this problem from our perspective.
Finally, one consequence of our construction relating F M (p, q) to FM (p, q) is that when these are also superconformal field theories, their partition functions differ in the sign of the "bottom component" of the super-partition function, that is Z RR = TR R ( (−1) F ) which is the trace of (−1) F on the highest weight space of the Ramond fields [40]. This cannot be determined on the grounds of modular invariance. We have found that in our construction, Z RR = +1 for FM (3,8) and Z RR = −1 for FM (3,8) which are hence two inequivalent superconformal field theories at c = −21/4. It is a curious fact that this implies that in this model the Virasoro characters satisfy χ 1,2 = χ 1,6 + 1, as is easily checked from the character formulae in [21]. The value Z RR = 0 for the TCIM was already observed in [40]. This leaves open now the question of investigating the superconformal models in which this parity-shift relates two inequivalent field theories, and how this might affect previous results on boundary conditions, boundary perturbations, etc.

JHEP06(2020)025
The bulk field content is as follows, giving both the conventional name from the freefermionic extension of the Ising model and the names following the conventions of sec- The F-matrix entries we need are: In the discussion of boundary states and defects we make use of the modular transformations of Virasoro characters, with q = exp(2πiτ ) andq = exp(−2πi/τ ), χ 0 (q) = 1 2 χ 0 (q) + 1 √ 2 χ σ (q) + 1 2 χ (q) ,

A.2 Tri-critical Ising data
The tri-critical Ising model is the Virasoro minimal model M (4, 5). The bulk model is described in [41] and the conformal boundary conditions in [31]. In table 21 we list the even and odd sectors in the same manner as for the Ising model, using mostly the naming conventions of [41] and [42] ( [43] use t instead of ). The even sector is the field content of the bosonic tri-critical Ising model. Note that the field G is conventionally normalised to G|G = 3c/2 = 21/20 (as opposed to our G|G = 1) and so we only have G ∝ φ o 1 . Likewise, ∝ iGḠ. As pointed out in [43], the tri-critical Ising model is related to the first non-trivial unitary superconformal minimal model SM (3,5) which is the associated fermionic model.