One-point Functions in Defect CFT and Integrability

We calculate planar tree level one-point functions of non-protected operators in the defect conformal field theory dual to the D3-D5 brane system with k units of the world volume flux. Working in the operator basis of Bethe eigenstates of the Heisenberg XXX_{1/2} spin chain we express the one-point functions as overlaps of these eigenstates with a matrix product state. For k=2 we obtain a closed expression of determinant form for any number of excitations, and in the case of half-filling we find a relation with the N\'eel state. In addition, we present a number of results for the limiting case of infinite k.


Introduction
The simplest probes of external heavy objects in a conformal field theory, such as Wilson or 't Hooft lines, surface operators or interfaces, are onepoint functions of local operators in the presence of the defect.By conformal symmetry, where z is the distance from x to the defect and ∆ is the scaling dimension of the operator O.The constant C in principle depends on the normalization of the operator at hand, but if the two-point function of O is unit-normalized, C is defined unambiguously.
Here we focus on a domain wall in N = 4 Super-Yang-Mills (SYM) theory which separates vacua with SU (N ) and SU (N − k) gauge groups [1].This defect originates from the D3-D5 brane intersection and is dual to a probe D5 brane in AdS 5 × S 5 with k units of electric flux on its world-volume [2].One-point functions of chiral operators in this [3] and in the closely related D3-D7 defect CFT [4], when continued to strong coupling perfectly agree with the predictions of the AdS/CFT duality.
We would like to make a connection with integrability and will thus consider expectation values of non-protected operators.It has proven useful in this context to study operators of large bare dimension, which correspond to long quantum spin chains.Conformal operators of this type, due to operator mixing, are linear combinations of a large number of field monomials.Efficient calculation of the classical expectation values for such operators becomes a non-trivial problem, which can only be solved by employing the full machinery of the Bethe ansatz.The one-point correlators are probably the simplest objects sensitive to the structure of the Bethe wavefunctions, and are thus ideally suited to probe integrability beyond the spectral data.

Domain wall and spin chains
The D3-D5 intersection defect in N = 4 SYM has the following semiclassical description at weak coupling.On the one side of the domain wall, the gauge symmetry is broken from SU (N ) to SU (N − k) by an infinite scalar vev.On the other side the scalar fields relax to zero according to their classical equations of motion: For a supersymmetric defect, the solution also satisfies the first-order Nahm equations [5]: which automatically imply (2.1).The solution describing the D3-D5 intersection is [6]: where the three k × k matrices t i satisfy and consequently realize the unitary k-dimensional representation of su (2).
The one-point functions, to the first approximation, are obtained by simply replacing quantum fields in the operator with their classical expectation values [3,4].To get a non-zero answer the operators must be built from scalar fields, and we will consider the most general such operators that do not contain derivatives: (2.5) The SO (6) tensor Ψ is cyclically symmetric because of the trace condition.These operators form a closed sector at one loop, and their mixing is described by an integrable SO(6) spin-chain Hamiltonian, wherein the tensor Ψ plays the role of the wave function in the spin-chain Hilbert space.The anomalous part of the dilatation generator (the mixing matrix) at one loop contains only nearest-neighbor interactions [7]: where λ = g 2 N is the 't Hooft coupling of the SYM theory, and P lm and K lm are permutation and trace operators acting on sites l and m of the spin chain: This result is not modified by the presence of the defect [8].Notice, however, that the latter reference deals with a probe brane set-up without fluxes (corresponding to k = 0) but the ultraviolet divergencies of the theory should be the same when the classical fields are turned on.The Hamiltonian (2.6) is a member of an infinite hierarchy of commuting charges responsible for the integrability of the model.The third charge1 of the hierarchy acts on three neighboring spins: Unlike the Hamiltonian, the third charge is parity-odd, and changes sign under the inversion of the spin chain orientation2 .The spectrum of the spin chain therefore contains parity pairs with degenerate energy and opposite values of Q 3 , as well as unpaired states with vanishing Q 3 .
The defect CFT contains also operators localized on the domain wall.These operators are described by an integrable open spin chain [8] and are dual to open strings with ends attached to the D5 brane.By considering one-point functions of the bulk operators we are, in a sense, dealing with the same string diagram but viewed as an absorption of a closed string by the D5 brane.In string theory the two descriptions should be related by t−s channel duality, and it would be interesting to understand how they are related at weak coupling.
By substituting (2.3) into (2.5)we find that the one-point function is proportional to the inner product of the spin-chain state Ψ that characterizes the operator and the state with the wave function (2.10) MPS here stands for the 'Matrix Product State', the term that will be explained below.The defect thus maps to a particular state in the spin-chain Hilbert space.We may interpret this state as a weak-coupling counterpart of the boundary state that describes the D5 brane in closed string theory.
Recovering the normalization factor that makes the bulk two-point function of O unit-normalized, we get for the structure constant: . (2.11) What can be said about the state associated with the defect?It is not an eigenstate of the spin-chain Hamiltonian.We do not get anything nice when apply (2.6) to (2.10).However, the third charge of the integrable hierarchy acts in a simple way and actually annihilates the defect state: (2.12) The proof is given in appendix A. This property leads to a selection rule for the one-point functions, since the overlap with MPS vanishes for all states that carry Q 3 ≠ 0.
To further simplify the problem we consider the SU (2) subsector composed of operators which are built from two complex scalars (2.13) The SU (2) sector is closed to all loop orders, and at the leading order is described by the Heisenberg spin chain.When restricted to the SU (2) sector, the spin-chain state associated with the defect becomes The index a is introduced here to distinguish the "auxiliary" space of color indices of t i from the quantum space spanned by ↑⟩, ↓⟩ on each site of the spin chain.The defect state (2.14) can be obtained by applying an operator, which we can call the defect operator, to the ferromagnetic ground state of the spin chain: The defect operator is not uniquely defined, because there are many operators that annihilate the ground state.We can choose it in the form where σ i l are the Pauli matrices acting on the l-th site of the spin chain, s is an arbitrary complex number, and t can be any k × k matrix.For instance, taking s = 0 and t = t 2 , we find: which takes particularly simple form for k = 2, with t 1 = σ 3 2 and t 2 = σ 1 2: States of the form (2.14) are known as the Matrix Product States, and were extensively studied in the condensed-matter literature [9], in particular to model quantum entanglement in one-dimensional systems.The operators (2.16), (2.17) and (2.18) are usually called the Matrix Product Operators.
In analogy to the algebraic Bethe ansatz (ABA) [10] the construction of the MPS uses the auxiliary space which threads through all sites of the spin chain.Interestingly, here the auxiliary space has a direct physical meaning of the color SU (N ) representation in the underlying gauge theory.
The conformal operators in the SU (2) sector are labelled by zero-momentum eigenstates of the Heisenberg Hamiltonian.In the ABA framework, the eigenfunctions are constructed by applying creation operators B(u) to the ferromagnetic vacuum of the spin chain: Of course this way to impose the zero-momentum constraint is too restrictive and there are zero-momentum Bethe states in which rapidities are not balanced pairwise.These states form degenerate parity pairs related by reflection of all rapidities.Such paired states, however, carry a non-zero Q 3 and have zero overlap with the defect state as a consequence of (2.12).We can thus concentrate on the fully balanced, unpaired states of the form (2.20).
Our goal is to calculate . (

2.21)
There is a considerable literature on overlaps of Bethe states in integrable systems (see [11] for reviews), which in many cases admit compact determinant representation.The most famous examples are the Gaudin norm of an ABA state [12], which is a part of the expression we need to evaluate, and the overlap of the on-shell and off-shell Bethe states [13].Overlaps of Bethe states with MPS have not been studied so far, to the best of our knowledge.
From known results the one that comes closest to our setup is the overlap of an arbitrary Bethe state with the Néel state, which was calculated in [14] and transformed into a convenient determinant form in [15].Bethe-state overlaps are playing an important rôle in the gauge/string integrability.The three-point functions in the N = 4 SYM at weak coupling can be expressed as generalized overlaps of Bethe states [16,17,18,19] and can be rendered into a compact determinant form [20], which is particularly useful in the semiclassical thermodynamic limit [18,21].An interesting question is whether the one-point overlap (2.21) also admits a determinant representation.
In this paper we investigate this question in the simplest case when the auxiliary space has dimension two (k = 2).We have found that the answer is affirmative, and moreover the result is given by exactly the same determinant formula as the overlap with the Néel state [15], upon relaxing the half-filling condition M = L 2 necessary to make the Néel overlap non-zero.The final result is written in terms of the matrices of size M 2 × M 2: and The structure constant (2.21) is given by the ratio of two determinants: . (2.24) When M = L 2, this formula coincides exactly with the expression for overlap between a half-filled Bethe eigenstate and the Néel state given in [15].Although the MPS is different from the Néel state, even if restricted to equal number of up and down spins, this is not a coincidence.We were able to show that the MPS is cohomologically equivalent to the Néel state at half filling and consequently has the same overlaps with all half-filled Bethe eigenstates.The result above then follows from the derivation in [14,15], for M = L 2. When M < L 2, this formula is a conjecture which we have extensively checked.We have also identified a natural generalization of the Néel state away from half-filling, which lies in the same cohomology class as the definite-spin projection of MPS.
In section 3 we introduce the tools necessary for our computation, namely the Bethe ansatz and an explicit realization of a set of k × k matrices which constitute a unitary k-dimensional representation of SU (2).Subsequently, in section 4 we sketch our computations and present the results.In section 5 we discuss the relationship between the MPS and the Néel state and introduce generalized Néel states with unequal number of up and down spins.Section 6 contains a discussion of the thermodynamical limit and section 7 some comments on the string theory observables dual to the one-point functions of the defect CFT.Finally section 8 contains our conclusion.

Setting up the computation
Although the construction of the defect state has a strong resemblance with certain elements of the algebraic Bethe ansatz we have found it most convenient to evaluate the overlaps by using the Bethe ansatz in its coordinate space version which we will summarize below, see for instance [22].Hereafter we will present the explicit representations of SU (2) that we will make use of in our computations.

The coordinate Bethe Ansatz
The eigenstates of the dilatation operator restricted to the SU (2) sector are in one-to-one correspondence with eigenstates of the Heisenberg XXX spin chain.In this section we introduce this model and discuss its solution via the coordinate Bethe ansatz.

Model
The XXX spin chain is a one-dimensional lattice model consisting of L spin- 1  2 particles.Therefore, the Hilbert space is ⊗ L C 2 , where each C 2 is spanned by ↑⟩, ↓⟩.The Hamiltonian describes a standard nearest neighbor spin-spin interaction with periodic boundary conditions L + 1 ≡ 1.For simplicity let us also introduce the the usual raising and lowering operators S ± such that Expressing the permutation operator in terms of spin operators one can see that (2.6) reduces to (3.1) in the SU (2) subsector, up to normalization.The (coordinate) Bethe ansatz gives us a method to diagonalize this Hamiltonian and to compute its spectrum.

Bethe eigenstates
The first step of the Bethe ansatz is to introduce a vacuum state This vacuum state is trivially an eigenstate of the Hamiltonian.The other eigenstates will also have down-spins on various sites.The Bethe ansatz postulates that these eigenstates are of a plane wave type.More precisely, each flipped spin behaves like a quasi-particle referred to as a magnon.These magnons propagate along the spin chain with some definite momentum p.
The Bethe eigenstate for a chain of length L describing M magnons, is of the form where N is an overall normalization.The sum over σ runs over all permutations of M elements.Furthermore, the factors θ parameterize the twomagnon S-matrix via It is worthwhile to note that, up to an overall normalization, the Bethe vector (3.4)only depends on the S-matrix S rather than the phase θ.In the remainder we will choose the normalization N such that the term e ip i n i (i.e. the term with σ = 1) in (3.4) appears with unit coefficient.In other words, we will set N = e − ∑ j<k θ jk 2 .
Bethe equations Finally, the state (3.4) should respect the correct boundary conditions, i.e. it should be periodic.Imposing periodicity results in a set of equations on the momenta of the magnons, called the Bethe equations When the momenta satisfy these Bethe ansatz equations, it is easy to check that the state (3.4) is an eigenstate of the Hamiltonian with eigenvalue where u = 1 2 cot(p 2) is the rapidity.In order for a Bethe eigenstate to represent a single trace gauge theory operator it is furthermore necessary that the momenta of its excitations add up to an integer multiple of 2π.This is required to account for the cyclicity properties of the trace, i.e.
Finally, notice that our Bethe states (3.4) (with N = e − ∑ j<k θ jk 2 ) are not normalized to unity.These coordinate space Bethe eigenstates can be related to the eigenstates of the algebraic Bethe ansatz approach in the following way (see, for example, [16]) This, in conjunction with the Gaudin formula [12] for the norm of {u i }⟩, fixes the normalization of coordinate Bethe ansatz eigenstates.
Overlap Let us now continue by computing the overlap between the Bethe states and the defect state ⟨MPS ⃗ p⟩.Inserting the M -magnon state (3.4) into (2.9)yields where the t i form the standard k-dimensional irreducible representation of su(2).However, for practical computations it is more convenient to take which will clearly yield the same results.

Representations of su(2)
Let us spell out the explicit representation for the su(2) generators t i that we will use and derive some useful relations.
Definition Consider the k-dimensional complex vector space generated by the basis vectors E i .Define the standard matrix unities E i j that are zero everywhere except for a 1 at position (i, j), such that they satisfy If we introduce the following constants and consider the matrices then we obtain the standard k-dimensional su(2) representation by defining It is easy to check that these matrices satisfy the su(2) commutation relations (2.4).Note that all these matrices are traceless.
Automorphisms Let us introduce two similarity transformations It is easy to show that under these transformations Hence, they provide a trivial automorphism of the algebra.

Results
In this section we present a number of explicit results for the overlap (3.11).

L or M odd
If L or M is odd, the overlap vanishes.This follows directly from the automorphisms (3.16).Indeed, for any state of the form tr[t n 1 3 t 1 t n 2 3 . ..], containing M t 1 's and L t's we have by cyclicity of the trace and This implies that the expression tr[t n 1 3 t 1 t n 2 3 . ..], and hence the overlap (3.11), is only non-vanishing if L and M are both even.

Vacuum, M = 0
From (3.14) we see that t 3 is a diagonal matrix with entries 1 2 (k − 2i + 1) for i = 1, . . ., k. From this, it immediately follows that for the vacuum state (3.3) the overlap (3.11) reduces to The resulting sum can be evaluated to a combination of ζ-functions Taking the k → ∞ limit of the explicit expression for ⟨MPS 0⟩ yields This agrees with the large k behavior which was found previously in [3,4].

Excited states 4.3.1 General considerations
We first notice that the defect state MPS⟩ is a cyclically invariant state (due to the cyclic nature of its expansion coefficients).This implies that where U = e i P is the lattice translation operator and P the momentum operator.From this we conclude that the overlap vanishes unless ⃗ p⟩ is a zeromomentum state.
Secondly, we notice that for an even number of excitations MPS⟩ is invariant under an operation traditionally denoted as parity, see for instance [23].Its action on a spin state is defined by where t i ∈ {↓, ↑}.The invariance of MPS⟩under this transformation follows from the invariance of its expansion coefficients under a similar operation performed on the matrices inside the traces.By an argument similar to the one above it follows that the overlap vanishes unless the Bethe eigenstate has positive parity.It is well-known that the eigenstates of the Heisenberg spin chain can be chosen to be eigenstates of a definite parity.In particular, the so-called un-paired eigenstates for which the Bethe rapidities fulfill that {u i } = {−u i } are automatically eigenstates with parity equal to (−1) M (L+1) .Moreover, as discussed in section 2, we find that only these unpaired state can have a non-trivial overlap with the classical function.This follows from the fact that the unpaired states are exactly the states that are annihilated by the odd charges Q 2n+1 .

Two excitations, M = 2
By using the cyclicity of the trace, we can rewrite the overlap (3.11) as a sum of terms of the form We can evaluate this trace by implementing the explicit expressions for t i (3.14) The definition of the matrix unities then allows us to work out the trace Thus, for M = 2, the Bethe states are mapped to The sums over m, n can easily be done and we find the following formula for the overlap 12) Notice that the above expression has to be evaluated with care in case k is odd due to a superficial pole at i = 1 2 (k − 1).By using that ⟨MPS p 1 , p 2 ⟩ is invariant if we redefine the summation via i → k − i it is easy to check that upon substituting the Bethe equations (3.6) the overlap vanishes unless p 1 + p 2 = 0 where the above expression has a pole.Then, imposing the vanishing of the total momentum and setting p 1 = −p 2 = p from the beginning gives us the following one-point function For k = 2 this reduces to 2 1−L Lu −1 (u − i 2 ).

General M
In the following we will derive some results for a general even number of excitations M .In particular, for the case k = 2, we will give a closed formula of determinant form, valid for any even M .k = 2 For k = 2 computing the overlap simplifies due to the identities The anti-commutator identity means that we can order the generators in the trace (possibly at the cost of a sign) and the first identity implies that we can take all the powers in the trace mod 2. In particular, we can simplify (3.11) to The above sum can be evaluated as follows In agreement with our general discussion, cf.section 2, we find that the only Bethe eigenstates that give a non-zero overlap function are states with momentum configurations of the form For these states one can write the overlap function in a compact form as the determinant of a matrix.Define the following function and the following M 2 × M 2 matrix then the overlap function is given by

.20)
We have confirmed this formula by explicit computations up to and including the case of eight excitations.Upon translating to the algebraic Bethe ansatz framework (cf.(3.9)), using the Gaudin formula for the norm, and applying elementary determinant identities, we arrive at the aforementioned result (2.24).
Large k Let us have a closer look at the leading order large k expansion for any number of excitations.One can show that for M excitations tr(t n 1 −1 This can be seen as follows.First, in the large k limit tr(t ) can be rewritten as a Riemann sum and integration then leads to the following identity Second, from the defining commutation relations of su(2) it can be seen that any distribution of t 3 , t ± under the trace can be ordered as (4.22) at the cost of terms of lower order in k.Then (4.21) follows by expressing t 1 in terms of t ± as in (3.15).This means that the large k limit of the overlap function reduces to It is easy to check that for M = 0 it reduces to the large k behavior we found for the vacuum state (4.5).However, for M ≠ 0 something unusual happens.Notice that (4.23) can be expressed as the inner product of the Bethe state (3.4) with the fully symmetrized state that has M spins down.Such a state can be expressed as the lowering operator S − acting on vacuum M times, i.e. ∆ (L) (S − ) M 0⟩.Thus, we can re-express the overlap as where ∆ is the coproduct.However, due to the fact that Bethe states are highest weight states, the above vanishes.In other words, the inclusion of excitations lowers the order of the overlap for large k.
In order to gain a better understanding of this phenomenon, let us look at the large k behavior for M = 0, 2, 4. We study the large k behavior by explicitly evaluating the relevant overlap function for a large range of values of L, k.The overlap will be a polynomial in k of degree at most L + 1 with coefficients that are rational functions of L. Letting L run from 2 to 20 and k from 2 to 30 allowed us to fix the relevant coefficients.In general, we find that the large k behavior is of the form The coefficient β (0) is constant and can be read off from (4.23).For M = 0 the first few β (m) are constant and from (4.4) the large k behavior is easily found to be L,2 is constant it vanishes by the same arguments as the leading order.For the other terms, the factors of n i can be written as derivatives of momenta p i when calculating the explicit overlap function.This allows us to evaluate the overlap (4.25) to the relevant order.Again we find that upon using the Bethe equations that it vanishes unless we impose pairwise momentum conservation.Doing this, we find for the next two terms Large k behavior of the one-point functions for M = 0, 2, 4 excitations.The order at which the expansion starts is k L+1−M .
Notice that, in contradistinction to the vacuum, there is a contribution at an even order.Finally, the next non-trivial contribution is Starting from k L−1 terms appear at both even and odd orders.Next, we turn to four excitations M = 4.It can be shown that the first order for M = 4 particles that contributes is k L−3 .This seems to indicate that the order at which the large k expansion begins is k L−M +1 .The first non-trivial coefficient for four particles can be computed along the same lines as for M = 2 and we find The general structure of the contributions is indicated in Table 1.

Matrix product and Néel states
In this section we elucidate the relationship between the matrix product and the Néel states.This will allow us to prove equation (2.24) for M = L 2. The Néel state is the vacuum of the classical (Ising) anti-ferromagnet: The state has equal number of up and down spins (we assume that the length L of the spin chain is even).On the other hand, the matrix product state has components with any even number of up and down spins.Since the total spin in conserved, it is ( To facilitate the bookkeeping, it is convenient to introduce the generalized MPS: where z is a complex number.Then, We can also generalize the Néel state to the case of an arbitrary even number of down spins: This looks like a descendant of the ground state, and would have been such, if not for the constraint that spin-flips hop by an even number of sites.Obviously, Another state that we shall deal with is a hybrid between the generalized Néel and MPS3 : where the product is understood in the cyclic sense, such that l m+1 ≡ l 1 and l = L + k is identified with l = k.Here π ± are chiral projectors in the auxiliary space: For instance, in the representation where t i = σ i 2, these are the ordinary spin-up/spin-down projectors: The generalized MPS can be pictured as a collection of m domains, separated by domain walls.Each domain wall carries a down spin in the quantum space and π ± projector in the auxiliary space.The sign of z flips across each domain wall (fig.1).Since π ± are projectors, the trace over the auxiliary space decomposes onto the product of matrix elements for each of the domains.The chirality of projectors enforces the domains to contain odd number of sites each.
The definite-spin projections of the generalized MPS, interpolate between the definite-spin components of the MPS and the generalized Néel states (5.5).Indeed, All these different states are related to each other, and in fact can be all expressed through the basic MPS (2.14) by simple projection and spinlowering operations.In particular, we will find that definite-spin components of the MPS are cohomologically equivalent to the generalized Néel states: where S i is the total spin operator, and S − is its lowering component that flips in turn all the spins in the chain with weight one.Since Bethe states are highest-weight: their overlaps with the MPS and the Néel states coincide: The determinant representation (2.24) in the case of M = L 2 then follows from the known overlap between the Bethe states and the ordinary Néel state [14,15].For other M , the overlap is given by the same equation, which we believe is a new result, that would be interesting to prove, either directly in the MPS representation or using its cohomological equivalence to the generalized Néel states (5.5).Now we proceed to prove (5.12).The proof rests on the following identity: Though not entirely obvious, this equation can be derived in a rather straightforward way.Both S − and d dz, when acting on MPS(z)⟩, produce l terms, where the l-th spin is flipped, in the former case with the coefficient t 1 and the latter case with the coefficient t 2 .Altogether, the action of id dz + S − creates a defect, a down spin accompanied by t + , where (5.16) Now, taking into account that we find that (5.15) immediately follows.
Applying the spin projection (5.10) to both sides of (5.15) we can express the generalized MPS through the ordinary one: The cohomological equivalence of the Néel states and the MPS state (5.12) is just a particular case of this relationship.

Classical limit
If the thermodynamic limit L → ∞ is accompanied by populating the spin chain with a large number of low-energy magnons, such that M L and u j L are kept fixed as L → ∞, the spin-chain states become semiclassical [24,25].
Oftentimes one can directly compare spin-chain results in this regime to classical string theory in AdS 5 × S 5 [25,26], even though the two approximations are supposed to work in the opposite range of the 't Hooft coupling.
In the scaling limit the Bethe roots concentrate on a number of cuts in the complex plane and are characterized by the density We are interested in symmetric configurations, due to the selection rules for the one-point function, and define the density by summing only over the right movers which constitute one half of all Bethe roots.The density satisfies an integral equation 2 ⨏ dy ρ(y) where n l are (positive) integer mode numbers, one integer for each arc of the Bethe root distribution.The general solution to these equations can be written in terms of Abelian integrals on an algebraic curve that characterizes a particular semiclassical state of the spin chain [27].We may ask how the overlap (2.21) behaves in this scaling limit.The non-trivial dependence of the overlap on the Bethe roots enters through the determinants of G ± , which are structurally similar to the Gaudin determinant.The thermodynamic limit of the latter was analyzed in [18] with the result that the leading contribution comes from the near-diagonal matrix elements, with i − j ≪ L. But in the ratio of determinants that enters the overlap formula (2.24) this contribution simply cancels, because the difference between the near-diagonal matrix elements of G + and G − is of order 1 (u j + u k ) 2 ∼ 1 L 2 and vanishes in the thermodynamic limit.One may then expect that the ratio approaches 1, with corrections of order 1 L.However, the situation is more subtle, and the ratio in fact approaches a finite constant value different from one: The coefficient K is given by the ratio of functional determinants: where are operators that act in the space of functions defined on the same set of arcs in the complex plane as the Bethe root density ρ(y).
The residual dependence on the density of Bethe roots arises because the original discrete determinants in (2.24) have a set of nearly zero modes, as already noticed in [18].These modes corresponds to vectors f j that are approximately constant on the scale j − k ≪ L. For such vectors the summation can be simply replaced by integration, and the matrices K ± and G ± in (2.22), (2.23) become integral operators: where Using the classical Bethe equations (6.2), the G ± operators can be further simplified to (6.5).
Apart from a trivial kinematic factor, the structure constant C u does not exponentiate in the thermodynamic limit.This is perhaps an indication that the limit of large M and L, at k = 2, is not really classical on the string side.Indeed, the natural classical limit in string theory would also involve taking k large (natural scaling is k ∼ √ λ at strong coupling [3]).We postpone a detailed study of this limit for future work, and just make a few general comments on possible comparison to string theory in the next section.

Comparison to string theory
Earlier studies of chiral primary operators have shown that one can expect an agreement between one-point functions calculated in gauge theory and onepoint functions calculated in string theory to leading order in the parameter λ k 2 in a double scaling limit where both λ and k are sent to infinity but the ratio λ k 2 is kept fixed and small.Hence, for this purpose one would mainly be interested in large representations.
The calculation of one-point functions on the string theory side was previously carried out in the case of chiral primary operators and involves computing the fluctuation of the probe D5 brane action due to fluctuations in the background supergravity fields when a source corresponding to the operator in question is inserted on the AdS boundary [3,4].The computation involved is completely analogous to the computation of a three-point function involving a chiral primary operator and two giant gravitons [28], and follows a general scheme of computing one point functions in the presence of a heavy probe, such as Wilson loops [29] or the three-point function of two heavy and one light operators [30].Performing the calculation of one-point functions involving other types of operators would require other techniques.One type of operators one could dream of considering could be BMN operators (i.e.two-excitation operators, considered in subsection 4.3.2).The string theory dual of these were given in [31].Another example could be the operator dual to a folded spinning string with two angular momenta on S 5 .This operator is characterized by its M ∼ O(L) Bethe roots being distributed on two arches placed symmetrically around zero [25,26] and belongs to the class of operators which have a non-vanishing overlap with the defect operators, cf.section 2. Both for BMN-and spinning string types of operators, however, it appears that the string theoretical calculation of the one-point function would be of a similar complexity as the computation of a three-point function involving three heavy operators.

Conclusion
We have seen a strong indication that the integrable structures underlying the duality between N = 4 SYM and type IIB string theory on AdS 5 ×S 5 leave an imprint on the correlation functions of the defect CFT derived from the D3-D5 probe-brane set-up with internal gauge field flux, k.We have concentrated our efforts on the calculation of one-point functions of non-protected operators and we have proposed, for k = 2, a closed expression of determinant form for the one-point function of Bethe eigenstates, based on explicit computations involving states with up to eight excitations.Furthermore, for half filling we have proved the formula by relating the matrix product state to the Néel state.Needless to say that it would be very interesting to construct a proof of the formula in the general case.
The formulation of the one-point function as an overlap involving a matrix product state could indicate interesting connections to problems in condensed matter physics.In addition, there are numerous other directions of investigation which could lead to further insights on the theme touched upon here.One-and multi-point correlation functions of defect CFT's with dual gauge field flux could be studied for higher values of k, to higher loop orders and for other probe brane set-ups, such as the D3-D7 case.Finally, it would obviously be very interesting if one could match any of these quantities with quantities derived in the dual string theory picture.
There are other cases in which heavy probes create a coherent field configuration in the CFT vacuum, which at weak coupling can be studied by semiclassical methods.This is the case for the 't Hooft loops [32], surface operators [33], and domain-wall defects [34].It would be interesting to investigate the spin-chain representation of one-point functions in these cases as well.

A Action of third charge on defect state
In this appendix we prove eq. (2.12).This is most easily done graphically.The Hamiltonian density H lm and the third charge Q lmn are shown in fig. 2  depicted in fig. 3. The total charge vanishes upon summation over l, which should be clear from the figure.

(4. 26 )
Notice that the even orders vanish.However, starting from M = 2 the coefficients become non-trivial.Let us list the first few β (m) L,2 and describe their contribution.If we denote n ij = n i − n j , then β
00082.The work of K.Z. was supported by the Marie Curie network GATIS of the European Union's FP7 Programme under REA Grant Agreement No 317089, by the ERC advanced grant No 341222, by the Swedish Research Council (VR) grant 2013-4329, and by RFBR grant 15-01-99504.

Figure 3 :
Figure 3: The result of application of the third charge to the defect state.The horizontal bar denotes the trace over the auxiliary space.The active sites are shown in thick black lines, while the spectator sites, unaffected by Q, are shown in blue.
. Applying Q l−1,l,l+1 to the defect state tr t i 1 ...t i L , we get(Q ⋅ MPS) ijk = δ kj t s t s t i − δ ij t k t s t s + 2δ ik t j t s t s − 2δ ik t s t s t j +2δ ij t s t k t s − 2δ jk t s t i t s + 4t k t i t j − 4t j t k t i ,where i, j, k are indices on sites l − 1, l, l + 1 and we have suppressed the rest of the wavefunction unaffected by the operator.Using the commutation relations(2.4)thiscanbe brought to the form(Q ⋅ MPS) ijk = δ ij t s t s t k + 2δ ij t k + 4iε ijs t s t k − δ jk t i t s t s − 2δ jk t i − 4it i ε jks t s ,