Resilience of the topological phases to frustration

Recently it was highlighted that one-dimensional antiferromagnetic spin models with frustrated boundary conditions, i.e. periodic boundary conditions in a ring with an odd number of elements, may show very peculiar behavior. Indeed the presence of frustrated boundary conditions can destroy the local magnetic orders presented by the models when different boundary conditions are taken into account and induce novel phase transitions. Motivated by these results, we analyze the effects of the introduction of frustrated boundary conditions on several models supporting (symmetry protected) topological orders, and compare our results with the ones obtained with different boundary conditions. None of the topological order phases analyzed are altered by this change. This observation leads naturally to the conjecture that topological phases of one-dimensional systems are in general not affected by topological frustration.

Complex systems and their different ordered phases have always attracted a large interest, not only from a purely speculative point of view but also for the innumerable technological applications that exploit their characteristics [1][2][3] . In the middle of the last century, all the different phases of many-body systems obeying the laws of classical mechanics were classified, using the Landau theory 4 . According with this theory, different phases are characterized by different order parameters. Each one of them is uniquely associated to a particular kind of order that is related to the violation of a specific symmetry (spontaneous symmetry breaking 5 ). It is therefore evident that, in the framework of Landau theory, the key role is played by the symmetries of the system, while other aspects, such as the boundary conditions, fall into the background and, generally, are considered not to be relevant for the presence or the absence of a particular kind of order.
When physicists began to extend their interest to the study of many-body models obeying quantum mechanics, Landau theory was immediately borrowed to study also the quantum regime 6 . Notwithstanding its undoubted success, after a few years, it became clear that it was unable to catch all the aspects of the richer landscape that comes out from the different quantum complex systems. Indeed, nematic ordered phases 7 show a situation in which the violation of the same symmetry is associated to different and non-equivalent order parameters, depending on the model under analysis 8 , while topologically ordered phases 9,10 do not admit local order parameters, but rather non-local ones. This implies that the classification of different ordered phases in terms of the local symmetry violated by the order parameter, which is the cornerstone of Landau's theory, is not able to provide a complete classification of the quantum phases of matter.
Moreover, in recent years, even the assumption that some aspects of complex systems, such as boundary conditions, do not affect the onset of an ordered phase has been questioned. Indeed, it was observed that particular boundary conditions, named frustrated boundary conditions (FBCs), i.e. periodic boundary conditions with an odd number of lattice sites, when paired with antiferromagnetic interactions, induce a topological frustration that can change several properties of the system such as the energy gap and the entanglement entropy 11,12 . Moreover the presence of the topological frustration can affect the exponential/algebraic dependence of the two-points correlation functions on the distance 11,13 , and prevent the formation of the AFM ordered phase that characterizes the same model without frustration 13,14 . Depending on the symmetry of the model, the frustration can give birth to a first-order phase transition that separates a mesoscopic ferromagnetic phase from a peculiar antiferromagnetic order, with a magnetization slowly varying in space and incommensurate with the lattice length 15,16 . In other models, the effects of the topological frustration can be even stronger coming to the point of preventing the creation of any local order 17 and thus inducing a change in the nature of the macroscopic phases, and related phase transitions, in the system 18 .
Behind this violation of one of the main prescriptions of Landau theory is that FBCs induce a geometrical frustration in the system 19,20 . Frustration is an extremely broad phenomenon, that can have a quantum and/or a geometrical/topological origin with a huge spectrum of effects 8,21,22 . However, at its deep core, frustration always comes out from the impossibility to satisfy simultaneously all local constraints in a many-body system 23- 25  www.nature.com/scientificreports/ usual way to understand frustration is to think about a classical system made of dichotomous objects, i.e. a classical equivalent of quantum spins that can assume only two values, namely ±1 , whose behavior is governed by a Hamiltonian with a two-body short-range Ising-like antiferromagnetic interactions. While, locally, there is no problem to minimize an antiferromagnetic Ising bond, when the dichotomous objects are arranged in loops made by odd numbers of sites, at least one of such bonds needs to display a ferromagnetic alignment. Despite its simplicity, this example captures all the main aspects of frustration. Indeed, even in more complex lattices, the presence of geometrical frustration can be traced back to the existence of frustrated loops made by an odd number of antiferromagnetic bonds 19 . While some phenomenology on the effects of geometrical frustration on antiferromagnetic phases in quantum many-body has been accumulated, nothing has still been done for phases exhibiting topological orders. In the present paper, our goal is to fill this gap. Therefore, we will focus on several exactly solvable models known to exhibit a weak form of topological order, namely a Symmetry Protected Topological Order (SPT), which can be broken by a perturbation not respecting the symmetry of the Hamiltonian without closing the gap [26][27][28][29][30] , and we analyze the effect of imposing frustrated boundary conditions on them.
At first, in "The Cluster-Ising model", we study the so-called Cluster-Ising model 31,32 that is a one-dimensional model, in which a three-body cluster interaction competes with an antiferromagnetic Ising one. Despite its apparent complexity, the model admits an analytical solution obtained by mapping the spins degrees of freedom into spinless fermions using the Jordan-Wigner transformation 33 . Since there is an interplay between two different types of interactions, the model presents a transition between a phase dominated by the cluster interaction and one in which the leading term is the Ising one. While in the second region the system (usually) admits a magnetically ordered phase, when the cluster interactions dominate over the antiferromagnetic ones, the model is known to exhibit a symmetry protected topological order. We will show that, while the magnetic phase results to be deeply affected by the rising of topological frustration, the SPT ordered phase is insensitive to the particular choice of boundary conditions, a fact that we explain considering simple properties of the two and three body interactions. The fact that the sensitivity of the model to the change of boundary conditions is a function of the macroscopic phase and not only of the model considered, strongly suggests that such property must extend to all other models that present SPT ordered phases.
To test such a suggestion we turn, in "AKLT model", to the AKLT chain 34,35 , which is a one-dimensional model, with a SO(3) symmetric Hamiltonian, describing spin-1 degrees of freedom interacting antiferromagnetically. The model supports a unique ground state (up to boundary state degeneracies), separated from the rest of the eigenstates by a finite gap, and it exhibits exponentially decaying correlation functions. It is characterized by spins paired into valence-bonds and symmetry protected topological order 28,30 . The model with periodic boundary conditions, including FBCs, has been already studied in details 35,36 so we simply summarize and discuss the results relevant for this work showing that also in this model frustrated boundary conditions do not affect the model and its order.
Finally, in "Kitaev chain" we consider last of the 1D topological models in our survey, namely the Kitaev chain 37 . With open boundary conditions, the Kitaev chain is known to be exactly mappable to the quantum XY chain in a transverse magnetic field 38,39 , but moving to periodic BC there are some subtleties, that play an important role when we wish to analyze the effect of frustration. In all cases, the Kitaev chain can be reduced to a free fermionic problem, and therefore diagonalized analytically 39,40 . Exploiting this approach we highlight that, regardless of the macroscopic phase, the Kitaev chain is completely unaltered by moving from an even to an odd number of sites (or viceversa).
Thus in all models analyzed we find that their symmetry protected topological orders are either resilient to the topological frustration imposed by frustrated boundary conditions, or that the system is such that frustrated boundary conditions do not introduce any incompatibility in the simultaneous minimization of the local and global structure of the Hamiltonian. These observations naturally drives us to conjecture that SPT ordered phases (and, a fortiori, topological ordered phases) in one dimension are not affected by frustrated boundary conditions. This result could be, to some extent, counterintuitive. In fact, it shows that topological phases are not affected by the (real-space) topology of the system. Most of all, our conjecture places a further clear distinction between phases characterized by non-global order parameters, as the magnetic and the nematic ones, and topologically ordered ones, at least in one dimension.

The Cluster-Ising model
The model. Let us start by introducing the Cluster-Ising model. We consider a system made of spins-1 2 , in which a two-body antiferromagnetic Ising pairing competes with a three-body cluster interaction. The Hamiltonian of such model reads where σ α j , for α = x, y, z , are the Pauli operators on the j-th spin and the parameter φ allows changing the relative weight between the two terms. Since our goal is to study the effect of topological frustration, we assume periodic boundary conditions σ α j+N = σ α j and that N is an odd number ( N = 2M + 1 ). Due to the existence of an analytical solution, the family of spin-1/2 cluster models was intensively studied in the past years 31,32,[41][42][43][44] . For the Cluster-Ising model in Eq. (1) it is well-known that for φ ∈ (−π/4, π/4) the model is in an antiferromagnetic Ising phase, for φ ∈ (3π/4, 5π/4) in a ferromagnetic phase, while for φ ∈ (π/4, 3π/4) and φ ∈ (−3π/4, −π/4) , where the many-body interactions dominate over the Ising ones, the model exhibits a symmetry protected www.nature.com/scientificreports/ topological order. Such order can be characterized by a non-zero expectation value for the non-local string operator, defined as Before we start the detailed solution of the model based on the Jordan-Wigner transformation, let us make some general considerations. Both from the seminal Toulouse's works 19,20 for classical models and from their generalization to the quantum regime 24,25 , we have that, to determine whether or not a model is geometrically or topologically frustrated, we need two elements. The first of these elements is a prototype model, i.e. a model in which frustration is absent. The second element is a set of local unitary operators, by which we try to reduce the model under analysis to the prototype one. If it is possible to find such a set of local unitary operations that map our Ising-like model into the prototype one, then the system is free from geometrical frustration. Otherwise, we have that the system is geometrically or topologically frustrated.
In the case of an antiferromagnetic Ising model the prototype model is the model in which each bond is turned into a ferromagnetic one. Let us now focus on the case of a one-dimensional lattice with periodic boundary conditions in which all bonds between neighboring spins are antiferromagnetic. Such model is exactly the Ising term of the Hamiltonian in Eq. (1) and can be written If N is even, the sign of all the terms σ x j σ x j+1 can be inverted simply by inverting every second spin, starting from j = 1 , hence reducing it to a purely ferromagnetic Ising model and thus proving that the model is not frustrated. On the contrary, with odd N, it is impossible to bipartite the system in two sets made by the same number of elements making it impossible to find a suitable transformation, acting on only one of the two subsets, that transforms it into a purely ferromagnetic model, hence proving the presence of frustration.
Let us now consider the cluster term of the Eq. (1), i.e. j σ y j−1 σ z j σ y j+1 . Differently from the Ising case, inverting every single spin of the lattice through a unitary operation generated by the spin operators σ y j it is always possible to change the global sign of the model. As a consequence, the three-body cluster interaction is not expected to show any geometrical frustration. This should be contrasted with the case of cluster interaction with an even number of sites, such as the 2-Cluster-Ising chain 43,44 , which do not have topological order. According to the argument above, in these cases we expect that frustration will be effective and indeed we have proven that to be the case 17,18 .
When we have the simultaneous presence of both the interactions, arguments like those just made become more complex. However, in the parameter region in which the Ising type interaction dominates over the cluster one, taking inspiration from arguments such as adiabatic deformation 45 , we expect the system to present signatures of the presence of geometrical frustration. Viceversa, in the region dominated by cluster-type interaction, these signatures can be expected to be absent. Such region is expected to show a symmetry protected topological order 31 , and therefore it is expected to be unaffected by the presence of frustration.
Let us now turn to the exact solution of the Cluster-Ising model in presence of frustrated boundary conditions. It is well known 31,43 that it can be diagonalized exactly, using the Jordan-Wigner transformation that maps spins into spinless fermions. In the process of diagonalization, details can be found in the Supplementary Information 1, the Hamiltonian is divided in the two parity sectors of � z = N j=1 σ z j , and in each sector the Hamiltonian can be written in terms of free fermionic operators where a q are Bogoliubov fermions. The fermionic momenta q in Eq. (5) belong to two different sets, respectively q ∈ Ŵ + = { 2π N (k + 1 2 )} for the even parity sector ( z = 1 ) and q ∈ Ŵ − = { 2π N k} for the odd one ( z = −1 ), where, in both cases, k runs over all integers between 0 and N − 1.
To each fermionic momentum is associated an energy, given by It is worth noting that the momenta 0 ∈ Ŵ − and π ∈ Ŵ + (since we study the case in which N is odd), are different from the others because: (a) they do not have a corresponding opposite momentum; b) their energies can be negative.
ε q = 2 1 + sin 2φ cos 3q ∀ q � = 0, π , www.nature.com/scientificreports/ From Eq. (6) it is easy to determine the ground states of the system starting from the vacuum of Bogoliubov fermions in the two sectors ( |0 ± � ), which, by construction, have positive parity z = 1 , and taking into account the modes with negative energy and the parity requirements. We are going to examine the antiferromagnetic phase and the cluster phase separately. To compute the ground state expectation of observables it is going to be convenient to use the Majorana fermions, defined as which are related to the spin operators as On the basis of Wick's theorem the expectation values of observables are determined by the two-point correlators of Majorana fermions.
The Majorana correlators in the ground state (9) are found to be while, for N not divisible by 3, the correlators can be reproduced from the previous expressions by taking formally u 2 = u 3 = 0 . The exact finite-N result for the last correlator can be found in the Supplementary Information 1.
In the correlators we can see corrections of order 1/N, that would not be present in a model without topological frustration. Although they vanish in the thermodynamic limit, they cannot be neglected, as they can influence the spin-correlation functions at large distances, i.e. for a distance that scales with the dimension of the system. Indeed, as shown in the Supplementary Information 1, we find the spin-correlation functions The correlations at large distances r decay linearly and for the most distant spins on the ring, separated by r ≈ N/2 , they vanish in the thermodynamic limit. Coherently also the magnetization order parameter vanishes in the presence of FBCs and antiferromagnetic interactions. Thus, topological frustration affects the antiferromagnetic phase of the Cluster-Ising chain, closing the energy gap and destroying the magnetization and spincorrelations at large distances, similarly to several other frustrated chains [11][12][13] .
Cluster phase. We find that for φ ∈ (π/4, 3π/4) the ground state is |g� = |0 + � and it is separated by a finite energy gap from the excited states above it. Similarly, for φ ∈ (−3π/4, π/4) the ground state is |g� = a † 0 |0 − � , also with a finite energy gap above it. These two regions of the model's phase space are those known to display symmetry protected topological order 31,43 .
As shown in the Supplementary Information 1, the correlators of Majorana fermions are �A j A l � g = �B j B l � g = δ jl and www.nature.com/scientificreports/ which is the same result as without frustration 31,43 , i.e. the same result would be obtained for even N. There is no corrections of order 1/N as in the antiferromagnetic phase. The consequence of this result is that the expectation value of any bulk observable remains the same with frustrated boundary conditions, generic periodic BC or, arguably, open BC. Indeed, in the antiferromagnetic phase the effect of frustration arises as a correction to the Majorana correlation function. Since, any observable can be expressed in terms of Pauli spin-operators, while Pauli spin-operators can be expressed as a product of Majorana fermions, the expectation value of any observable can be expressed as an expectation of a product of Majorana fermions, which is by Wick theorem determined by two-point correlators of Majorana fermions. Therefore, since the two-point correlators of Majorana fermions are the same as without frustration for large N, and since the same applies for the Jordan-Wigner transformation (7), so is the expectation value of any bulk observable.
In particular, in the Supplementary Information 1, we compute the expectation value of the string operator. We obtain as without frustration 43 . On the contrary, in the topologically ordered phases, the expectation values of the operators σ α j for α = x, y, z are zero. Namely, since the ground state does not break the parity symmetry z of the model we have immediately �σ x j � = 0, �σ y j � = 0 , while the relation �σ z j � = 0 follows from the equality σ z j = −ıA j B j and the property that the corresponding integral in (12) vanishes. The non-zero expectation value of the non-local string operator, and zero expectation value of local observables, such as spin operators, characterizes the symmetry protected topological ordered phase in the Cluster-Ising model. It is, hence, proved that such a phase is not affected by topological frustration.

AKLT model
The AKLT model 34,35 is a one dimensional model, describing spin-1 degrees of freedom interacting antiferromagnetically. It is defined by the SO(3) symmetric Hamiltonian The FBCs are achieved by imposing an odd number of lattice sites N = 2M + 1 and periodic BC � S j+N = � S j . Naively, we can expect that FBCs induce frustration in the AKLT model, by considering its classical limit, and thus neglecting the second term in Eq. (14). With periodic boundary conditions and even N the energy is clearly minimized by the perfectly staggered configuration S j = (−1) j or its flipped version. However, for odd N, corresponding to FBCs, such a configuration is impossible because the lattice is not bipartite. Therefore, not all local terms in the Hamiltonian can be minimized simultaneously, and the Hamiltonian is frustrated.
Turning back to the quantum case, the AKLT Hamiltonian in Eq. (14) acquires a (symmetry protected) topological order and it is thus an ideal candidate to test our conjecture that this property protects it from the effect of FBCs. It is easy to check that Eq. (14) can be written as a sum of projectors where P (2) ( S j , S l ) projects the state of two spin-1 at sites j and l into their spin-2 representation.
The AKLT model with both open and periodic BC has been studied in detail in 35 (for a more pedagogical approach see the book 36 ). It is known that the ground state is unique with PBC, and four-fold degenerate with open BC, with this degeneracy related to the existence of edge states and thus not influencing the expectation values of bulk observables, which are the same in the different ground states (similarly to what happens in the Cluster-Ising model and the Kitaev chain that we will analyze later). The AKLT ground state, with periodic boundary conditions, is a valence-bond state. To show this, one represents each spin-1 degree of freedom through two spin-1/2 in their triplet representation. Then, these spin-1/2 on neighboring sites are paired as singlet to prevent the ferromagnetic alignment penalized by the Hamiltonian. Denoting by |α� j , |β� j the two spin-1/2 on the j-th site, the valence bond state is . We then construct the ground state of (14) as www.nature.com/scientificreports/ where projects the two spin-1/2 into their spin-1 representation, with the Clebsch-Gordan coefficients P +1 ≡ 1 0 0 0 , P 0 ≡ 1 √ 2 0 1 1 0 , and P −1 ≡ 0 0 0 1 , and where the N + 1 site is identified with the first, since periodic boundary conditions are assumed. As already discussed in 35 , the only difference between having even or odd N is the need of a different index contraction in the spinorial representation of the valence bond state. However, the even-odd choice does not change the gapped nature of the system and the bulk behavior of the correlation functions 35,36 .
In particular, the string order parameters that encode the topological order of the AKLT model for α = x, y, z , remain non-zero and unchanged. We thus conclude that the (symmetry protected) topological order of the AKLT model is not affected by frustrated boundary conditions.

Kitaev chain
The Kitaev chain 37 is a model of a spinless fermions topological superconductor for which the Hamiltonian reads where µ is the chemical potential, w is the hopping amplitude and is the superconducting gap. As we have done so far, also with the Kitaev model we will focus on the case of FBCs, i.e. with periodic boundary conditions ( c j+N = c j ) and an odd number of lattice sites N = 2M + 1 . However, before we start the analysis of the case with periodic boundary conditions, let us summarize the main results that were obtained for the open ones.
It is well-known that the Kitaev chain with open boundary conditions, that can be obtained from Eq. (21) restricting the range of the second sum up to j = N − 1 , can be mapped, inverting the Jordan-Wigner transformation in Eq. (3), to the quantum XY chain in transverse field 39,46 Hence, exactly as the correspondent spin model, in the thermodynamic limit, it shows a phase transitions at µ = ±2w . Such quantum phase transition separates a topologically trivial phase for |µ| > 2|w| without edge modes in the open chain, from a topologically ordered phase |µ| < 2|w| characterized by the presence of Majorana edge modes. The latter, for the spin chain, was shown to be affected by FBCs when the Ising term promotes an AFM order, with a phenomenology similar to that of the Cluster-Ising chain we discussed in "Antiferromagnetic phase". Due to the relation between the XY and the Kitaev chain, it is natural to wonder if the latter can also be affected by FBCs, although we cannot compare the spinless model to a related classical one which is affected by frustration.
Moving from open to periodic boundary conditions, the exact equivalence between the fermionic and the spin model ceases to exist. The reason for such a quite surprising result is connected to the fact that the Jordan-Wigner transformations break the invariance under spatial translation, by selecting a reference site. This implies that the interaction terms between the first and the last spin of a chain are no more mapped in a standard two-body fermionic term, but in a string term that makes it impossible to map a short-range fermionic model into a short range spin model. To provide an example, we have that the term σ x N σ x 1 is mapped into the string where z is the parity operator along the z axis that has support on the whole lattice. A similar result stands also for σ y N σ y 1 . When either w ± � < 0 (and |µ| < 2|w| ), the XY chain, with an odd number of sites, becomes frustrated: the energy gap above the ground states closes (algebraically) in the thermodynamic limit 11 , the correlation functions acquire peculiar algebraic corrections and the entanglement entropy violates the area law 12 , and the AFM local order is replaced by either a ferromagnetic mesoscopic order or by a AFM incommensurate modulated one (with a quantum phase transition separating the two) 13,15 .
From a physical perspective, the main difference between the Kitaev and the XY chain is that only the former supports (symmetry protected) topological order. Given the close relationship between the two models and the (17)  www.nature.com/scientificreports/ fact that the effects of frustration have already been established for the spin chain, we may wonder whether there is an even-odd effect also in the Kitaev chain. As presented in the Supplementary Information 1, exploiting the approach illustrated in Ref. 39 , we can diagonalize the Kitaev chain with periodic boundary conditions and odd N obtaining where a q are Bogoliubov fermions, whose momenta belong to the set Ŵ = { 2π N k} , with k running over integers between 0 and N − 1 . It is worth noting that, assuming N to be odd, k = π is not an allowed momentum in this model. The dispersion relation is given by Similarly to the Cluster-Ising model case, the mode q = 0 is special because it is the only one in which the energies can be negative. The eigenstates of the model are constructed by populating the vacuum state |0� . Taking into account the dispersion relation, it is easy to see that the ground state of the Kitaev chain with periodic frustrated boundary conditions, is always non-degenerate, with a finite energy gap above it, except at the phase transition points µ = ±2w , where the spectrum becomes gapless and relativistic. Similarly to the Cluster-Ising model, in the topologically ordered phase |µ| < 2|w| , the ground state degeneracy, in the thermodynamic limit, is different from the one that characterizes the model with open boundary conditions. While the former is two-fold, with PBCs the ground state is unique. Nevertheless, the expectation values of bulk observables are the same, as we now show.
For this purpose, we define Majorana fermions using Eq. (7). Then, similarly to the Cluster-Ising model, all operators acting on the Fock space generated by c † j can be expressed in terms of Majorana fermions, and using the Wick theorem it follows that the ground state expectation value of any observable is determined by the twopoint correlators of Majorana fermions. As shown in the Supplementary Information 1, the two-point correlators of Majorana fermions in the ground state, in all parameter regions of the model, are �A j A l � g = �B j B l � g = δ jl and Hence, for large N, the expression of the Majorana correlation functions do not depend on the boundary conditions (the result for free boundary conditions follows from the equivalent spin chain 46 ). Therefore, the expectation values of all bulk observables in the Kitaev chain with FBCs are equal to those in other settings and no difference emerges when w ± � < 0 . We conclude, in particular, that topological order in the Kitaev chain is, as in the Cluster-Ising model, not affected by frustrated boundary conditions.

Conclusions
We have presented an analysis of the effects of frustrated boundary conditions on different (symmetry protected) topologically ordered phases characterizing one-dimensional models. At first, we have focused on the onedimensional Cluster-Ising model with an odd number of spins and periodic boundary conditions. We presented general arguments by which the symmetry protected topological order of the cluster phase is not expected to be affected by frustration, while the antiferromagnetic phase is. These speculative arguments have been confirmed by our analytic results. While in the antiferromagnetic phase FBCs close the energy gap and destroy both the spin-correlations at large distances and the magnetization, the string order parameter in the cluster phase is not affected by FBCs. The property that the effects of topological frustration are lost when the cluster interactions start dominating over the antiferromagnetic ones is, in some extent, similar to the situation in the frustrated Ising model 11,12 , where the effects of frustration are suppressed by increasing the magnetic field, resulting in the resilience of the paramagnetic phase to geometrical frustration.
Our results on the Cluster-Ising model are even more interesting if we observe that we can directly transfer our general arguments to the m-cluster Ising model, studied in Refs. 43,44 , which consists of m-body cluster interaction competing with the antiferromagnetic Ising pairing. It is known that for any odd m, the model is characterized by a symmetry protected topologically ordered phase, which is, by our general arguments, expected not to be affected by geometrical frustration. A question that arises naturally already at this point is whether FBCs do not affect only models with cluster interactions or it is a general property of systems supporting topological order. For this reason, we have extended our analysis to two additional one-dimensional models that also exhibit SPT, such as the AKLT model and the Kitaev chain. In both cases, the topological phase is unaffected by FBCs, as is also shown by the fact that a typical effect of frustration is the closing of the energy gap [11][12][13]15 , which remains open in the topological phases of all the models analyzed. Since we have found in various models that the topologically ordered phases analyzed so far are not affected by FBCs, we arrive naturally at the conjecture: frustrated boundary conditions do not affect symmetry protected topological phases of one dimensional systems. However, whether the conjecture is correct or not it remains to be established for certain. One could also think that only Landau local orders are thus sensitive (and fragile) to topological frustration: in 18 we show that this is not the case and that nematic order can be destroyed by FBCs.