Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: the quantum ising chain is a useful theoretical laboratory in which fundamental issues concerning quantum many - body systems can be thoroughly investigated , exploiting the exact knowledge of several features of its phase diagram and quantum correlations . many results for its low - energy properties have been derived in the quantum ordered and disordered phases , and , in particular , at the quantum critical point separating the two phases , both in the thermodynamic limit and in the finite - size scaling ( fss ) limit with several types of boundary conditions , see , e.g. , refs . @xcite and references therein . in this paper we extend these analytic results to the case in which boundary fields are present . we present a detailed study of the quantum ising chain in a transverse magnetic field @xcite in the presence of magnetic fields coupled with the order - parameter spin operator , located at the boundaries of the chain . we assume the two magnetic fields to have the same strength and consider two cases : ( i ) the two fields are parallel ; ( ii ) the two fields are oppositely aligned . as expected , in the quantum paramagnetic phase the boundary fields do not change the large - size behavior of low - energies quantities . at the critical point , bulk behavior is independent of boundary conditions . however , the magnetic fields induce a surface phase transition with a corresponding scaling behavior . the quantum ferromagnetic phase is more interesting . if the boundary fields have opposite direction , one observes two different bulk phases . for small magnetic fields , the ground state is ferromagnetic as it occurs in the absence of boundary interactions . on the other hand , kink states @xcite are the relevant low - energy excitations for large boundary fields . the two different phases are separated by a continuous transition that is only characterized by the nature of the two coexisting phases and , indeed , the same transition occurs in ising rings in the presence of a localized link defect @xcite .    because of the quantum - to - classical mapping , our results can also be applied to the two - dimensional ising model , and , more generally , to any model in the ising universality class , in a strip geometry . in the two - dimensional case the quantum transition between the ferromagnetic and the kink phase corresponds to the wetting transition , which , in the context of a strip geometry , is sometimes identified as an interface localization / delocalization transition @xcite . our results , therefore , provide new analytic expressions for two - dimensional ising systems at the wetting transition in the presence of boundary magnetic fields . the paper is organized as follows . in sec . [ model ] we introduce the one - dimensional quantum ising chain with boundary fields . in sec . [ fermmodel ] we compute the low - energy spectrum by exploiting the equivalent quadratic fermionic formulation of the hamiltonian @xcite . explicit calculations are reported in the following sections . [ paramag ] , [ critpoint ] and [ magnphase ] report results for the quantum paramagnetic phase , at the quantum critical point and for the ordered magnetized phase , respectively . [ fktra ] is devoted to the study of the magnet - to - kink transition driven by the boundary fields in the ordered phase . we obtain exact results for the gap , the magnetization profile , the two - point spin - spin correlation function , and the entanglement entropy . in sec . [ comp-2d ] these results are compared with the existing ones for two - dimensional classical models . finally , in sec .  [ conclu ] we summarize the main results and draw our conclusions . a number of appendices report the derivations of some of the results . the hamiltonian of the quantum ising chain in a transverse magnetic field is given by @xmath1 where @xmath2 are the pauli matrices . in the following we assume ferromagnetic nearest - neighbor interactions with @xmath3 , and @xmath4 . the ising chain undergoes a continuous transition at @xmath5 @xcite , separating a quantum ferromagnetic phase ( @xmath6 ) from a quantum paramagnetic phase ( @xmath7 ) . in this paper we investigate the effects of boundary magnetic fields aligned along the @xmath8 axis . they give rise to an additional energy term @xmath9 that is added to hamiltonian ( [ isc ] ) . in the following we shall consider the model in the specific case of equal parallel boundary fields ( pbf ) @xmath10 and of equal oppositely - aligned boundary fields ( obf ) @xmath11 it is not restrictive to assume @xmath12 in both cases . we will often use a basis in which @xmath13 is diagonal . states will be labelled as @xmath14 , where @xmath15 is the eigenvalues of @xmath16 . signs will be fixed so that @xmath17 has the form @xmath18 in this basis . in the absence of boundary fields , the hamiltonian commutes with the generator @xmath19 of the @xmath20 transformations @xmath21 and with the operator @xmath22 of the reflection transformations defined by @xmath23 . these operators do not generally commute with @xmath24 . however , note that @xmath25=0 $ ] for obf .    in this paper we analyze the low - energy spectrum of the model . in particular , we compute the energy differences between the lowest states and the ground state @xmath26 ( here @xmath27 are the energy eigenvalues ordered so that @xmath28 , and @xmath29 is the gap ) , the local magnetization and the two - point correlation function , @xmath30 for obf , because of the symmetry under @xmath31 , we have @xmath32 and , as a consequence , the average magnetization @xmath33 always vanishes . we also define the integrated correlation @xmath34 and the correlation length @xmath35 with respect to the center of the chain : @xmath36 to determine the spectrum of hamiltonian ( [ isc ] ) , we extend the model , considering two additional spins located in 0 and @xmath37 and the hamiltonian @xmath38 this is the ising hamiltonian with two different couplings on the boundary links and zero transverse field on the boundaries . let us note that the hamiltonian @xmath39 commutes with both @xmath40 and @xmath41 , which can therefore be simultaneously diagonalized . the hilbert space can be divided into four sectors , which we label as @xmath42 , @xmath43 , @xmath44 and @xmath45 , where @xmath46 are the eigenvalues of @xmath40 and @xmath41 . the restriction of @xmath39 to each sector gives rise to the hamiltonian @xmath47 , defined in eq . ( [ isc ] ) , with a boundary term of the form ( [ hb ] ) . hence , the spectrum of @xmath39 also provides the spectrum of @xmath48 . let us also note that @xmath39 is @xmath49 symmetric . indeed , if @xmath50 , then @xmath51 = 0 $ ] . since @xmath52 does not commute with @xmath40 and @xmath41 , the spectrum is necessarily degenerate . moreover , since @xmath52 anticommutes with these two boundary operators , @xmath52 maps sector @xmath42 to @xmath45 and @xmath44 to @xmath43 , so that the restriction of @xmath39 to @xmath53 allows us to compute the full spectrum of @xmath39 .    to compute the spectrum of hamiltonian ( [ def - hextended ] ) , we follow ref .  @xcite . we first perform a jordan - wigner transformation , defining fermionic operators @xmath54 and @xmath55 @xmath56 where @xmath57 . these relations can be inverted , obtaining @xmath58 thus , hamiltonian ( [ def - hextended ] ) becomes @xmath59 with @xmath60 for @xmath61 . in this formalism @xmath62 and @xmath63 . to go further , let us rewrite @xmath39 as @xmath64,\ ] ] with @xmath65 and @xmath66 symmetric and antisymmetric matrices , respectively . then , we perform a bogoliubov transformation . we introduce new canonical fermionic variables @xmath67 where @xmath68 and @xmath69 are fixed by the requirement that @xmath47 takes the form @xmath70 with @xmath71 following ref . @xcite , we define the vectors @xmath72 the variables @xmath73 satisfy canonical anticommutation relations if the vectors @xmath74 form an orthonormal basis , and so does the set @xmath75 . the vectors @xmath75 satisfy @xmath76 which turns the determination of the energies @xmath77 into an eigenvalue problem . if @xmath77 does not vanish , @xmath74 is given by @xmath78 if @xmath77 is zero , @xmath74 is the null eigenvector of @xmath79 . it is also possible to evaluate the constant @xmath80 in eq . ( [ h1-diag ] ) , which provides the energy of the ground state : @xmath81 the matrix @xmath82 is given by ( we write the @xmath83 matrix @xmath84 for @xmath85 , the generalization to any @xmath0 being obvious ) @xmath86 this matrix @xmath84 has a zero eigenvalue , @xmath87 , with eigenvector @xmath88 . correspondingly , we obtain @xmath89 and @xmath90 which depends only on the boundary fermion operators . note that it can be rewritten as a combination of @xmath40 , @xmath41 , and @xmath52 , hence it does not represent an additional symmetry of the hamiltonian . the presence of a zero eigenvalue was expected , as the spectrum is degenerate . if @xmath91 and @xmath92 do not vanish , since there is only one fermion operator with zero energy , the spectrum is doubly degenerate . let us now consider the nonzero eigenvalues . since the vectors @xmath75 are an orthonormal set , orthogonality with @xmath93 implies @xmath94 . it follows analogously @xmath95 . therefore , we have for @xmath96 @xmath97 where the dots represent a polynomial in @xmath98 and @xmath99 with @xmath100 . it is then easy to verify that @xmath101 = 0.\ ] ] it follows that , if a state @xmath102 satisfies @xmath103 , then we have @xmath104 . analogously , if @xmath105 , we have latexmath:[$\sigma_{l+1}^{(1)}\eta_k |\psi\rangle = s_1 \eta_k    sector @xmath46 , then @xmath107 belongs to sector @xmath108 .    to conclude , we should determine to which sectors the two degenerate ground states belong . if @xmath91 and @xmath92 are both positive , we expect the ground - state configurations to be ordered . therefore , if on one boundary the spin is directed in the positive @xmath8 direction , we expect the same to occur at the other end of the chain . therefore , we conclude that the ground states belong to sectors @xmath109 and @xmath45 . this identification is supported by exact diagonalization . we can now classify the states . we only consider the states with @xmath110 , to avoid the double degeneracy . then , if @xmath111 we have :    * in the sector @xmath112 , @xmath113 , the lowest energy state is the ground state of the hamiltonian and all states are obtained as @xmath114 , @xmath115 , with @xmath116 even . the first excited state is @xmath117 and the energy gap is @xmath118 . * in the sector @xmath112 , @xmath119 , the lowest energy state is the first excited state @xmath120 of the hamiltonian @xmath39 and all states are obtained as @xmath121 , @xmath115 , with @xmath116 odd . in particular , the first excited state in the sector is @xmath122 , so that the energy gap is @xmath123 . the first case is relevant when considering @xmath124 in the boundary hamiltonian @xmath24 , which are identified with @xmath91 and @xmath92 , respectively . the second case is relevant when @xmath125 and @xmath126 have opposite signs : the corresponding gap is obtained by identifying @xmath127 and @xmath128 . the energy spectrum of hamiltonian @xmath39 can be computed exactly , as discussed in [ app.a ] . here we will focus on the case @xmath129 , which allows us to derive the spectrum for parallel or opposite equal magnetic fields at the boundary , cf . ( [ pbfdef ] ) and ( [ obfdef ] ) , respectively . neglecting the irrelevant zero eigenvalue , we obtain @xmath37 elementary excitations whose energy @xmath130 is given by @xmath131 where @xmath132 are the @xmath37 solutions of the equation @xmath133 \sin k l = \label{secular - equation}\\ & & \quad g ( 1 - 2 g \cos k + g^2 - \zeta^4 ) \sin [ k ( l+1 ) ] . \nonumber\end{aligned}\ ] ] we should consider real solutions in @xmath134 ( correspondingly @xmath135 ) , purely imaginary solutions @xmath136 with @xmath137 ( correspondingly @xmath138 ) , and solutions @xmath139 ( correspondingly @xmath140 ) . equation  ( [ secular - equation ] ) has also a spurious solution for @xmath141 for all values of the parameters , which should be discarded . for @xmath142 , there is also a second solution with @xmath143 , which corresponds to a true excitation of the model . for @xmath144 we obtain the equation appropriate for open boundary conditions ( obc ) @xcite , @xmath145 and @xmath146 , which implies an additional zero mode for @xmath39 . this is not surprising as the spectrum has a fourfold degeneracy for @xmath147 : the four sectors are equivalent in the absence of boundary fields . the real solutions of eq . ( [ secular - equation ] ) are obtained by solving @xmath148\cos k + 2 g^2 \cos^2 k}. \label{realtankl}\ ] ] the imaginary solutions are obtained by setting @xmath149 . ( [ secular - equation ] ) then becomes @xmath150 ^ 2 \over               ( g - e^h ) ( e^h g - 1 + \zeta^2)^2 } . \label{equazione - j1eqj2-small}\ ] ] the corresponding energies @xmath151 are always smaller than those corresponding to real momenta , hence the localized states are the most important ones in the determination of the low - energy spectrum . for @xmath152 ( paramagnetic phase ) the lowest energy state is a localized state , i.e. , a solution of eq . ( [ equazione - j1eqj2-small ] ) with @xmath137 . for large values of @xmath0 we must consider the poles in the right - hand side of the equation . they are @xmath153 for @xmath152 , the second equation does not have any solution . the relevant pole is the first one , so that the relevant solution can be written as @xmath154 . for large values of @xmath0 , substituting this expression in eq . ( [ equazione - j1eqj2-small ] ) , we obtain @xmath155 correspondingly , we find @xmath156 the energy @xmath157 vanishes with exponentially small corrections . the energy of the second excited state corresponds to the lowest value @xmath158 of @xmath159 satisfying eq . ( [ realtankl ] ) . we obtain @xmath160 and @xmath161 the system is therefore gapped , and the same gap is obtained for both pbf and obf modulo exponentially small corrections , since @xmath162 : @xmath163 note that the same result would have been obtained for obc , confirming that boundary conditions and/or boundary fields are irrelevant in the paramagnetic quantum phase . at the critical point @xmath165 , there are no localized solutions , hence the two lowest energies @xmath166 are obtained by considering eq . ( [ realtankl ] ) . the equation becomes @xmath167 the lowest - energy solutions correspond to momenta that scale as @xmath168 . therefore , we can expand the right - hand side in powers of @xmath159 , obtaining @xmath169 which shows that @xmath170 must diverge for @xmath171 . this fixes @xmath172 with @xmath173 including the corrections , we obtain for the two lowest momenta @xmath174 so that @xmath175 these results imply that the pbf and obf gaps are @xmath176 respectively . note that the amplitude of the @xmath177 term differs from that obtained in the obc case ( corresponding to @xmath178 ) @xcite , @xmath179 indicating that @xmath144 is a surface critical point .    as we have discussed in ref . @xcite in the obc case , the corrections of order @xmath177 to the leading behavior can be interpreted as the effect of a nonlinear scaling field associated with @xmath0 : the leading irrelevant boundary operator gives corrections that scale as @xmath180 . the same holds in the presence of boundary fields . indeed , we can define a new length scale @xmath181 so that @xmath182 using the previously reported expressions we obtain @xmath183 note that the same rescaling applies both to pbf and obf , a result which is not obvious in the original formulation , but which has a natural explanation in the present setting in which we consider @xmath39 . indeed , the rescaling should cancel the leading correction in the whole low - energy spectrum of @xmath39 , hence in all low - energy excitation energies @xmath77 . therefore , it should apply to both pbf and obf . since @xmath144 is a surface critical point and the boundary term is a relevant perturbation , we can study the crossover behavior close to @xmath144 . the relevant scaling variable is @xmath184 . therefore , in the limit @xmath185 , @xmath171 at fixed @xmath186 we have @xmath187 where the function @xmath188 depends on the type , obf or pbf , of boundary fields . in the scaling limit defined above , momenta scale as @xmath189 . using eq . ( [ seceq - g1 ] ) , we obtain an implicit equation for @xmath190 : @xmath191 therefore , to compute the scaling function , one determines the two lowest values of @xmath190 satisfying eq . ( [ scaling - b - boundary ] ) and uses them to determine the gap as before . the scaling curves are reported in fig . [ fss - crossover - etab ] . let us now consider the scaling behavior in a neighborhood of the critical point @xmath165 . the corresponding scaling variable is @xmath192 the relevant equations are obtained from eqs . ( [ realtankl ] ) and ( [ equazione - j1eqj2-small ] ) by expanding the right - hand sides in the limit @xmath193 , @xmath171 , @xmath194 at fixed @xmath195 and @xmath196 . as expected , as long as @xmath197 , we obtain a result that is independent of @xmath198 : the fss behavior only depends on the nature of the boundary conditions , but not on the specific values of the boundary fields .    a localized state exists only for @xmath199 . if we set @xmath200 , the parameter @xmath201 is a solution of the equation @xmath202 the corresponding excitation energy is @xmath203 note that @xmath204 for @xmath205 , so that @xmath206 in the limit . on the other hand , for @xmath207 we obtain @xmath208 so that @xmath209 , consistently with eq . ( [ cale1_para ] ) . the propagating modes can be written as @xmath210 , where @xmath211 are the positive solutions of the equation @xmath212 with @xmath213 ( the solution with @xmath214 only exists for @xmath215 ) . the corresponding excitation energies are then @xmath216 . using these results we obtain for the energy gap in the fss limit : @xmath217 the two signs refer to the case of parallel ( plus sign ) and opposite ( minus sign ) boundary fields . the curves are reported in fig . [ fss - delta1 ] together with the corresponding obc scaling function . in the paramagnetic phase all curve have a similar shape  they coincide asymptotically  while in the phase @xmath218 ( @xmath219 ) , the behavior is clearly different . for @xmath220 , we have @xmath221 , @xmath222 , and @xmath223 . it is interesting to discuss the scaling corrections , that can be easily shown to decay as @xmath177 , as it occurs for @xmath144 . in ref . @xcite we argued that such corrections are not associated with irrelevant boundary operators , which give rise to corrections of order @xmath180 , but that they should be interpreted as due to nonlinear scaling fields . the same holds for @xmath197 . first , one should consider the effective length @xmath181 defined in eq . ( [ leff - def ] ) . then , one should consider the nonlinear scaling field @xmath224 associated with @xmath225 and , finally , rescale the energies with the _ sound velocity _ @xmath226 . the nonlinear scaling fields do not depend on the boundary conditions , hence we can use the expressions reported in ref . @xcite : @xmath227 then , we define the new scaling variable @xmath228 and express @xmath211 and @xmath201 as functions of @xmath229 . a somewhat lengthy calculation shows that @xmath230 can be written as in eq . ( [ gap - gcritico ] ) , with @xmath229 replacing @xmath195 , with corrections that scale as @xmath231 . as expected , the leading corrections cancel out . [ magnphase ]    let us now consider the case @xmath218 . in the infinite - volume limit , there exist two degenerate ground states @xmath232 and @xmath233 , which differ by the value of the magnetization  @xcite @xmath234 where @xmath235 is a bulk magnetic field applied along the @xmath8 direction . for a chain of finite size @xmath0 , the degeneracy is lifted . in the obc case , the energy difference between the two lowest - energy states vanishes exponentially as @xmath0 increases @xcite : @xmath236 . \label{deltaobc}\end{aligned}\ ] ] moreover , @xmath237 by symmetry , and @xmath238 these are the standard features of a quantum ferromagnetic phase with a spontaneously broken @xmath49 symmetry . we now extend the discussion to the case in which boundary fields are present , considering equal pbf and obf . to identify localized states , we consider eq . ( [ equazione - j1eqj2-small ] ) . the relevant solutions are associated with the poles in the right - hand side of the equation , see eq . ( [ equazione - poli ] ) . it is immediate to verify that there are no poles for @xmath239 , while a pole exists in the opposite case . therefore , we distinguish three different cases , depending on whether @xmath198 is smaller , larger , or equal to @xmath240        for @xmath241 , the two lowest - energy excitations corresponding to @xmath157 and @xmath242 are localized . we have @xmath243 ( @xmath137 ) with @xmath244 where @xmath245 and @xmath201 is a correction term . let us first obtain a large-@xmath0 expression for @xmath201 that is uniform as @xmath246 . this requires particular care as the right - hand side of eq . ( [ equazione - j1eqj2-small ] ) has a different singular behavior for @xmath144 and @xmath247 . for @xmath144 , it has a simple pole as @xmath248 . on the other hand , for @xmath249 , it has a double pole for @xmath250 .    to compute @xmath201 we set @xmath251 in all terms of eq . ( [ equazione - j1eqj2-small ] ) that are not relevant for the singular behavior as @xmath171 , i.e. , we write @xmath252 with @xmath253 equation ( [ eq - magn-1 ] ) is a second - order algebraic equation of @xmath254 , whose solutions are @xmath255 the energies @xmath256 are correspondingly @xmath257 and @xmath258 , @xmath259 . there are now two interesting cases . first , we take the limit @xmath260 at fixed @xmath261 . we can then simplify eq . ( [ epm - general ] ) , obtaining @xmath262 these results imply @xmath263 in the pbf case the gap is finite and proportional to @xmath198 for small fields . instead , @xmath264 vanishes exponentially , as in the obc case ( which corresponds to vanishing boundary fields ) . a second possibility consists in considering the fss limit around @xmath144 . as discussed in ref . @xcite , the ratio @xmath265 is expected to become a function of @xmath266 ( to compare with the formulae of ref . @xcite note that the energy associated with the perturbation is here proportional to @xmath198 , without additional factors of @xmath0 ) . in this limit @xmath267 is proportional to @xmath268 . using the expression for @xmath269 reported in eq . ( [ deltaobc ] ) , we obtain @xmath270 so that @xmath271 note that the scaling function agrees with that predicted in ref . we can also consider the ratio @xmath272 in the same limit , but in this case we obtain that the ratio is simply equal to 1 , i.e. , independent of @xmath268 . the reason can be easily understood . for obf the energy associated with the perturbation is zero  the contributions of the two boundaries cancel  hence @xmath273 identically .      for @xmath274 there are no localized states . the secular equation ( [ realtankl ] ) can be rewritten as @xmath275 for @xmath171 , the solutions can be parametrized as @xmath276 it follows that @xmath277 for large values of @xmath0 we obtain finally @xmath278 note that the correction term @xmath279 diverges as @xmath280 and has a finite limit as @xmath281 . moreover , the explicit form of the hamiltonian is only relevant for the corrections . it follows @xmath282 we obtain finally @xmath283 for obf , the gap vanishes for @xmath171 . however , while the approach is exponential for @xmath241 , we have @xmath284 for large fields . for pbf , the system is gapped as for @xmath285 . however , in this case corrections are of order @xmath180 . the gap equation is particularly simple . we obtain a solution @xmath287 and @xmath288 hence @xmath289 with no @xmath0 dependence . as for the other states we obtain @xmath290 with @xmath291 . it follows @xmath292 , as for @xmath293 . however , corrections , of order @xmath180 , are a factor of five smaller compared to the large - field case . for obf we obtain the asymptotic behavior @xmath294 as in the large - field case , @xmath264 scales as @xmath180 . however , the prefactor differs by a factor of 3 from that obtained for @xmath295 , cf . ( [ deltalk ] ) . we wish now to compute the ground - state magnetization and the corresponding two - point spin - spin correlation function in the magnetized phase , i.e. , for @xmath296 . in the pbf case , the system is magnetized , hence we expect @xmath297 and @xmath298 to be both independent of @xmath299 and @xmath300 , except possibly close to the boundaries , and to be equal to @xmath301 and @xmath302 , respectively , where @xmath301 is given in eq . ( [ modef ] ) . thus , the boundary magnetic fields have the only role of lifting the degeneracy . the obf case is more interesting . in [ app.b ] we have fully characterized the ground state of the model for small values of @xmath303 . we find that for it is ferromagnetic for @xmath285 ( this is consistent with the results for the gap discussed above ) . for @xmath293 kink states @xcite are the relevant low - energy excitations , while for @xmath304 the ground state is a superposition of kink and ferromagnetic states . the exact knowledge of the ground state allows us to compute the ground - state magnetization and the correlation function perturbatively in @xmath303 . for obf we find that @xmath305 vanishes in the low - field case except at the boundaries , depends linearly on @xmath299 for @xmath286 , and varies in a simple fashion for @xmath274 . we expect the same space dependence for all finite values of @xmath306 in the large-@xmath0 universal limit , modulo a multiplicative renormalization . if we set @xmath307 and @xmath308 , so that the center of the chain corresponds to @xmath309 , we predict @xmath310   \displaystyle{m(x)\over m_0 } = \displaystyle{x\over \ell }      & \qquad { \zeta = \zeta_c } , \\[2 mm ]   \displaystyle{m(x)\over m_0 } = 0 & \qquad { \zeta < \zeta_c}. \end{array } \label{magnetization - scaling}\ ] ] these relation should hold in the limit @xmath311 , @xmath312 at fixed @xmath313 , except possibly at the boundary , i.e. , for @xmath314 . here @xmath301 is the bulk magnetization of the system in the infinite - volume limit defined in eq . ( [ modef ] ) . it is interesting to note that the linear behavior of the magnetization for @xmath315 is due to the fact that the ground state is translationally invariant : the space probability of the kinks is independent of @xmath8 . due to the presence of the boundary fields , the translation invariance of the kink space distribution is not trivial and , indeed , it does not hold for @xmath295 . to check these predictions and study the approach to the asymptotic behavior , we perform numerical simulations using the density matrix renormalization - group ( dmrg ) method  @xcite . for convenience , we consider chains with odd @xmath0 , setting @xmath316 ( since we are considering ferromagnetic interactions , the low - energy properties do not depend on this feature , unlike the antiferromagnetic case , see , e.g. , ref . the results for the magnetization are shown in fig . [ magnetization - fasi ] . for @xmath285 , the magnetization vanishes in a large interval around @xmath317 , assuming positive and negative values only close to the boundaries . note that the region in which @xmath318 shrinks as @xmath0 increases , conferming the validity of eq . ( [ magnetization - scaling ] ) in the large-@xmath0 limit at fixed @xmath313 . for @xmath293 we show results for two values of @xmath303 . once we renormalize @xmath319 by using @xmath301 , i.e. , we consider the ratio @xmath320 , results fall one on top of the other , confirming universality . the results are in full agreement with the expressions ( [ magnetization - scaling ] ) . the same arguments apply to the correlation function @xmath298 . extending the small-@xmath303 results of [ app.b ] to the whole low-@xmath303 phase , we predict @xmath321 in the large-@xmath0 limit . therefore , the correlation function should be constant in the magnetized phase , @xmath322 , it should depend linearly on @xmath323 for @xmath304 , while it should have a sinusoidal dependence in the kink phase . again , we check these predictions by comparing them with dmrg data . they are reported in fig . [ gx - fasi ] . the dmrg data approach the theoretical predictions with increasing @xmath0 , confirming the theoretical predictions . as already discussed in ref . @xcite , in the case of obf and for @xmath218 , an interesting universal scaling behavior is observed for @xmath280 . such a transition is uniquely characterized by the nature of two phases occurring for @xmath285 and @xmath324 : on one side the ground state is ferromagnetic , while on the other side kink states are the relevant low - energy excitations . this characterization is confirmed by the analysis of the ising chain on a ring with a bond defect : at the boundary of the two phases one observes @xcite a behavior analogous to the one discussed here . here , we will report the computation of the energy gap the result was already presented in ref . @xcite  and additional results for the local magnetization , correlation function , and entanglement entropy . moreover , we will discuss the density of the fermionic modes in the equivalent fermionic picture of the model @xcite , showing that the transition can be interpreted as a localization transition of the fermionic states at the boundaries . the crossover behavior around the transition is parametrized in terms of the scaling variable @xcite @xmath325 then , in the large-@xmath0 limit we have @xmath326 in general , the scaling functions @xmath327 and @xmath328 are expected to be universal modulo a @xmath303-dependent normalization ( which is independent of the observable considered ) of the argument @xmath329 . the @xmath303 dependent prefactor appearing in eq . ( [ etasdef ] ) is the required nonuniversal factor . it is determined in sec . [ sec7.1 ] , analyzing the scaling behavior of the gap @xmath264 . according to renormalization - group theory , this normalization factor is independent of the observable considered , hence also the scaling functions of other observables should be independent of @xmath303 apart from a multiplicative factor ( not present if renormalization - group invariant ratios are considered ) , once @xmath329 is used as scaling variable . the results we will present for the magnetization and the two - point correlation function confirm this prediction .    beside the gap @xmath330 , we also consider the integrated correlation @xmath34 defined in eq . ( [ chixidef ] ) . it satisfies the scaling relation @xmath331 where @xmath332 is universal , i.e. , @xmath303 independent . these results allow us to define the rg dimensions of the different operators . in general , we expect @xmath333,\qquad \chi = a_\chi l^{d - y_\sigma } f_\chi [ l^{y_\zeta } ( \zeta-\zeta_c ) ] . \label{scalzetas}\ ] ] comparison with our results allows us to identify @xmath334 , @xmath335 , and @xmath336 . the latter result is consistent with the fact that a bulk magnetic field coupled to @xmath337 drives a first - order quantum transition , see , e.g. , ref . @xcite .    by defining @xmath329 as in eq . ( [ etasdef ] ) we have implicitly assumed that we are fixing @xmath303 and varying @xmath198 close to the critical value @xmath338 . in some contexts , however , it is more natural to consider a different point of view . the boundary field strength @xmath198 is fixed ( it should satisfy @xmath339 ) and @xmath303 is varied in a neighborhood of the critical value @xmath340 . in this case we have @xmath341 where the @xmath198 dependent prefactor guarantees that scaling functions are the same for all values of @xmath198 .      to determine the scaling behavior of the energy gap @xmath342 , we first simplify the secular equation ( [ secular - equation ] ) , by taking the limit @xmath343 at fixed @xmath329 . in this limit @xmath159 scales as @xmath177 , so that we set @xmath344 . for @xmath171 the secular equation becomes @xmath345 expressing @xmath346 as a function of @xmath347 , we find that the solutions of eq . ( [ sec - etac-1 ] ) satisfy one of the two equations @xmath348 if @xmath349 and @xmath350 are the positive solutions of eqs . ( [ eq - tanz2 - 1 ] ) and ( [ eq - tanz2 - 2 ] ) , respectively , it is easy to verify that @xmath351 for @xmath352 and @xmath353 . for @xmath354 we have instead @xmath355 . for @xmath356 there are also localized solutions with @xmath357 . the parameter @xmath358 satisfies one of the two equations @xmath359 eq . ( [ eq - tanw2 - 2 ] ) has a positive solution @xmath360 for @xmath356 , while eq . ( [ eq - tanw2 - 1 ] ) has a positive solution @xmath361 satisfying @xmath362 only for @xmath363 . these results allow us to determine the scaling function defined in eq . ( [ sca ] ) . we obtain @xmath364 for @xmath365 , @xmath366 for @xmath367 and @xmath368 for @xmath363 . we can easily determine the asymptotic behaviors . for small @xmath329 we have @xmath369 while for @xmath370 we have @xmath371 finally , for @xmath372 we have @xmath373 the theoretical curve is shown in fig . [ deetac ] . the dmrg data confirm the theoretical calculation : the numerical estimates of the ratio @xmath374 clearly approach the scaling curve @xmath375 with increasing @xmath0 . let us now consider the behavior of the correlations close to @xmath376 . the local magnetization @xmath297 and the correlation function @xmath298 are computed exactly in the limit @xmath377 in [ app.b4 ] . as we have already done in sec . [ sec6.2 ] , we extend these expressions to all values of @xmath303 satisfying @xmath6 , by replacing @xmath297 with @xmath378 and @xmath298 with @xmath379 . therefore , we predict for @xmath365 @xmath380,\ ] ] where @xmath381 is the smallest positive solution ( @xmath382 ) of eq . ( [ eq - tanz2 - 2 ] ) , @xmath383 , and @xmath308 , as before . for @xmath356 we have instead @xmath384,\ ] ] where @xmath358 is the solution ( @xmath385 ) of eq . ( [ eq - tanw2 - 2 ] ) . as for the correlation function @xmath298 it satisfies eq . ( [ gij - m ] ) . it is important to stress that the previous expressions , although derived for small values of @xmath303 , are expected to be _ exact _ for any @xmath218 . for instance , since @xmath320 is a dimensionless renormalization - group invariant ratio , it should scale as @xmath386,\ ] ] where @xmath387 is a universal , hence @xmath303 independent , function . the specific features of the model enter in this expression only through the normalization nonuniversal factor @xmath388 and the critical value @xmath338 . the factor @xmath388 does not depend on the observable , hence it can be determined by using any quantity . we use the gap @xmath264 , for which we are able to obtain exact results for any @xmath218 , defining @xmath389 . therefore , @xmath387 as a function of @xmath329 and @xmath390 is @xmath303 independent . therefore , our results , obtained in the limit @xmath391 , hold for any value of @xmath303 satisfying @xmath218 . the previous expressions for @xmath319 hold in the limit @xmath311 , @xmath171 at fixed @xmath313 and smoothly interpolate between the expressions valid in the different phases . note that we have @xmath392 for @xmath393 in all cases . this result does not however necessarily apply to the boundary points , i.e. , for @xmath394 , as the magnetization is not continuous at the boundary . for instance , in the magnetized phase , @xmath395 in all interior points , while @xmath396 and @xmath397 are nonuniversal functions of the field strength @xmath198 . a graph of @xmath320 for several values of @xmath329 as a function of @xmath313 is reported in fig . [ magnxcr ] . as expected , for @xmath398 the magnetization approaches zero exponentially in an interval that is centered in the middle of the chain and that widens as @xmath399 increases . more precisely , for @xmath372 we have @xmath400 and @xmath401 which shows that @xmath319 vanishes except in two tiny intervals close to the boundaries of width @xmath402 .    to verify the prediction ( [ gij - m ] ) for @xmath298 we compute the scaling functions associated with @xmath34 and @xmath35 . they are reported in [ app.b4 ] and compared with dmrg data in figs . [ chietacr ] and [ xietacr ] . results with different values of @xmath303 fall on top of the theoretical curve when plotted versus @xmath329 . this confirms our predictions for @xmath298 and the magnetization , as well as the correctness of the nonuniversal prefactor appearing in the definition of @xmath329 . it is interesting to reinterpret our results in the equivalent fermionic picture of the model . in the ferromagnetic phase , i.e. , for @xmath403 , the lowest eigenstates are superpositions of majorana fermionic states localized at the boundaries @xcite . in finite systems , their overlap does not vanish , giving rise to the splitting @xmath404 . the coherence length @xmath405 diverges at the kink - to - magnet transitions as @xmath406 , a behavior analogous to that observed at the order - disorder transition @xmath407 where @xmath408 .    within this picture , it is natural to consider @xmath409 which can interpreted as the particle density in the fermionic representation of the ising chain , see eq . ( [ sigma - c - relation ] ) . its large-@xmath0 limit is given by  @xcite @xmath410 we consider the difference @xmath411 where @xmath383 , as before . such a quantity has been computed for small values of @xmath303 in [ app.b4 ] . the universality argument we have used for @xmath297 and @xmath298 can be applied to @xmath412 . we therefore predict @xmath413 for @xmath365 and @xmath414 for @xmath356 , where @xmath381 and @xmath358 are determined as in the case of the local magnetization . to verify this prediction we compute @xmath415 numerically for @xmath416 and @xmath417 . the numerical results are compared with theory in fig . [ fig : deltan0 ] . the agreement confirms the analytic prediction .    the subtracted fermionic density @xmath418 is reported in fig .  [ deltanxcr ] as a function of @xmath313 . for positive values of @xmath329 it has a maximum at the origin and shows that fermions are delocalized . for @xmath419 , @xmath420 is constant , while for negative values of @xmath329 , there is a significant enhancement at the boundary . for @xmath372 , as in the case of the magnetization , @xmath418 is different from 0 only close to the boundary , in an interval of width @xmath421 . it is interesting to consider the behavior of @xmath422 at the boundaries , i.e. , for @xmath423 , which is shown in fig . [ deltanb ] . in the kink phase , such a quantity scales as @xmath424 for @xmath370 . on the other hand it diverges as @xmath425 for @xmath372 , that is when the magnetized phase is approached . as such , it represents a physically relevant order parameter for the transition , which distinguishes the two different phases .      finally , we discuss the behavior of the half - chain entanglement entropy  @xcite . we divide the chain into two connected parts @xmath426 $ ] and @xmath427 $ ] and consider the von neumann entropy @xmath428 } \ln \rho_{[-\ell,0 ] } , \label{vnen}\ ] ] where @xmath429}$ ] is the reduced density matrix of the subsystem @xmath426 $ ] , i.e. , @xmath430}={\rm tr}_{[1,\ell ] } \;\rho , \ ] ] and @xmath431 is the density matrix of the ground state . the entanglement entropy @xmath432 approaches a constant in the large-@xmath0 limit , except for the critical case @xmath5 where it increases logarithmically  @xcite . its large-@xmath0 limit for obc is given by  @xcite @xmath433 , \nonumber\\ & & i(z ) = \int_0 ^ 1 { dx\over \sqrt{(1-x^2)(1-z^2 x^2 ) } } . \label{sobc}\end{aligned}\ ] ] in fig . [ entaen ] we plot the difference @xmath434 which shows the scaling behavior @xmath435 where @xmath436 is a universal function apart from a multiplicative nonuniversal constant . these results imply that the large-@xmath0 limit of the half - chain entanglement entropy does not depend on the boundary conditions , being the same for obc and obf . however , the @xmath437 corrections are affected by the boundary term ( [ hb ] ) , giving rise to a nontrivial universal behavior around @xmath376 . the results we have obtained allow us to derive exact scaling functions for the classical two - dimensional ising model in a strip @xmath438 in the presence of two equal magnetic fields on the boundaries . this model has been extensively studied , especially in the case of opposite magnetic fields . in this case , the wetting ( or interface localization - delocalization ) transition , which is the analog of the magnet - to - kink transition we are considering here , has been extensively studied @xcite . we can therefore check some of our results , comparing them with the existing ones in the literature . the results for the gap @xmath264 obtained in sec . [ sec7.1 ] can be used to derive the scaling behavior of the longitudinal correlation length @xmath439 , defined by @xmath440 . such a quantity was computed in the restrictive solid - on - solid model in ref . their result is completely consistent with ours . if our scaling quantities @xmath381 and @xmath358 defined in sec . [ sec7.1 ] are identified with the quantities @xmath441 , @xmath442 appearing in ref . @xcite , we obtain exactly the same expression for @xmath443 . several papers study the magnetization profile , which gives information on the behavior of the interface between the two coexisting phases . the behavior in the large - field phase was determined in ref . @xcite by using the restrictive solid - on - solid model . the result is in agreement with our prediction ( [ magnetization - scaling ] ) for @xmath293 . the linear behavior observed at @xmath304 was also derived at the wetting transition by using the solid - on - solid model @xcite . no exact results are instead available for the magnetization profiles in the whole crossover region , which , in the classical case , is usually parametrized in terms of @xmath444l$ ] , where @xmath235 is the boundary field and @xmath445 the corresponding wetting temperature ( in the quantum - to - classical mapping @xmath303 corresponds to @xmath22 , hence this scaling corresponds to varying @xmath303 at fixed @xmath198 ) . however , the numerical results for the magnetization reported in refs .  @xcite ( and , in particular , fig .  2 of ref . @xcite ) appear to be in full agreement with our prediction shown in fig . [ magnxcr ] . we consider the one - dimensional quantum ising chain in the presence of a transverse magnetic field @xmath303 @xcite and of boundary magnetic fields aligned with the order - parameter spin operator . in particular , we consider fields with equal strength @xmath198 , which have either the same ( pbf ) or the opposite ( obf ) direction , cf . ( [ pbfdef ] ) and ( [ obfdef ] ) . we assume the chain to have length @xmath0 and derive the finite - size behavior of several quantities in the limit @xmath171 as a function of @xmath198 . we derive analytic predictions for the gap in all phases as a function of the boundary field strength . in the paramagnetic quantum phase , the leading behavior is independent of the boundary conditions . at the critical transition @xmath165 , finite - size scaling depends on the boundary fields . the point @xmath144 , corresponding to open boundary conditions , is a surface critical point . we derive analytic expressions for the scaling functions that parametrize the crossover behavior for @xmath193 and any @xmath198 ( bulk critical behavior ) , and for @xmath165 and @xmath246 ( surface critical behavior ) . in the quantum ferromagnetic phase @xmath218 , if the boundary fields are oppositely aligned , the finite - size behavior drastically changes as @xmath198 varies @xcite . for small @xmath198 , the system is in a ferromagnetic phase in which the gap decreases exponentially with the system size . on the other hand , for large fields , kink propagating modes of momenta @xmath177 are the relevant low - energy excitations , so that the gap decreases as @xmath446 . the two phases are separated by a critical transition at @xmath286 . close to this transition , low - energy properties show a universal scaling behavior in terms of the scaling variable @xmath447 . the transition is only characterized by the nature of the coexisting phases . indeed , the same transition occurs in ising rings in the presence of a localized bond defect @xcite . it is interesting to interpret our results in the context of two - dimensional classical ising models , using the quantum - to - classical mapping . as is well known @xcite , the quantum ising chain of length @xmath0 corresponds to a classical ising model on a strip of width @xmath0 , the transverse field playing the role of the temperature . therefore , our results provide exact predictions for the ising model in an @xmath448 geometry with boundary fields . our results for the gap at the critical point can be used to predict the fss behavior of the longitudinal exponential correlation function @xmath449 as a function of @xmath450 and of @xmath451 , where @xmath65 and @xmath66 are nonuniversal model - dependent constants . moreover , our extensive results at the magnet - to - kink transition for @xmath218 for obf provide predictions for the finite - size scaling functions at the classical wetting transition ( sometimes also called interface localization / delocalization transition @xcite ) in two dimensions . in particular , we have exact predictions for the fss behavior of the magnetization profile and of the correlation function of the layer magnetization in the whole crossover region close to the transition .    in the context of the wetting transition , one may also consider a second important surface parameter , which is usually referred to as surface enhancement @xcite . while a boundary magnetic field is the most relevant magnetic perturbation that breaks the @xmath20 symmetry at the boundary , surface enhancement is the most relevant ( energy - like ) perturbation that is @xmath20 invariant . in the context of the quantum chain , such an interaction can be mimicked by adding a term ( @xmath452 , localized at the boundary , to the general ising - chain hamiltonian . it is not difficult to generalize the results presented here to this more general case ( results for the spectrum in the absence of boundary magnetic fields appear in refs . @xcite ) . at the critical point , this would allow us to study the behavior at the ordinary and extraordinary surface transitions ( in the absence of surface enhancement the point @xmath165 , @xmath144 corresponds to the so - called special transition @xcite ) . moreover , a complete understanding of the phase behavior for @xmath218 in the absence of a bulk magnetic field should be possible . in the quantum case , these results would also be relevant for ising quantum rings with localized site defects , a distinct case with respect to that considered in ref . in this appendix we collect results on the spectrum of some matrices that are relevant for the discussion . we begin by considering the @xmath453 matrix ( we write it for @xmath454 ) @xmath455 to determine its eigenvalues , we compute @xmath456 if we set @xmath457 , the determinant @xmath458 satisfies the recursion relation @xmath459 the solution is @xmath460 where @xmath461 and @xmath462 are arbitrary constant . if we require @xmath463 and @xmath464 , we obtain @xmath465 it follows @xmath466 the eigenvalues are solutions of the equation @xmath467 . we obtain @xmath468 eigenvalues given by @xmath469 with @xmath470 . let us now compute the eigenvectors . let @xmath471 be the eigenvector corresponding to eigenvalue @xmath472 . it satisfies the relations @xmath473 the solution of the recursion relation ( [ eigenvectors - relation-2 ] ) is @xmath474 where @xmath475 and @xmath476 are two arbitrary constants . requiring @xmath477 , we obtain @xmath478 , so that @xmath479 where we have rescaled @xmath475 for convenience . requiring the eigenvector to have unit norm , we obtain @xmath480 we also determine the eigenvalues of the matrix @xmath481 which differs from @xmath65 only for two matrix elements : @xmath482 and @xmath483 are equal to @xmath484 . setting again @xmath485 and @xmath457 , we obtain @xmath486 using the recursion relation ( [ recursion - an ] ) , we obtain finally @xmath487 where @xmath488 .      to determine the spectrum of hamiltonian ( [ def - hextended ] ) , we compute the eigenvalues of the matrix @xmath84 defined in eq . ( [ matrixc - def ] ) . we write them as @xmath489 note that nothing forbids @xmath490 to be larger than @xmath491 or smaller that @xmath492 , hence @xmath493 is not necessarily bounded between @xmath494 and @xmath495 . hence , we must consider real values of @xmath159 with @xmath496 and complex values of the form @xmath497 and @xmath498 with @xmath137 . we are not interested in the zero eigenvalue of @xmath84 , hence we will only consider the @xmath499 matrix formed by the first @xmath500 rows and columns , which we label with the symbol @xmath501 . then , we introduce a matrix @xmath502 defined as ( here @xmath503 ) @xmath504 the corresponding secular equations are @xmath505 the determinant @xmath506 has already been computed in [ app.a1 ] : @xmath507 to determine @xmath508 , we expand the determinant with respect to the first and the last rows . we obtain the following recursion relation : @xmath509 a simpler expression is obtain by eliminating @xmath510 and @xmath511 , and expressing the result in terms of @xmath512 and @xmath513 . this can be done by using the recursion relation satisfied by the determinants @xmath506 , see eq . ( [ recursion - an ] ) with @xmath514 . we obtain finally @xmath515 d_{n-2}. \nonumber \end{aligned}\ ] ] note that the recursion relation ( [ rec - cn - eta1-eta2 ] ) is symmetric with respect to the exchange of @xmath91 and @xmath92 . the eigenvalues of @xmath516 are obtained by requiring @xmath517 . setting @xmath518 and @xmath519 , we obtain eq . ( [ secular - equation ] ) . note that we have simplified @xmath520 in the denominators , hence eq . ( [ secular - equation ] ) has a spurious solution for @xmath143 . a true solution with @xmath141 occurs only when the coefficients of @xmath512 and @xmath513 in eq . ( [ rec - cn - eta1-eta2 ] ) both vanish for @xmath141 . this only occurs when @xmath521 . the low - energy behavior of the model with hamiltonian ( [ isc ] ) can be understood analytically in the limit @xmath391 . for @xmath522 , it is trivial to compute the spectrum of hamiltonian ( [ isc ] ) in the presence of the magnetic boundary term ( [ hb ] ) , since @xmath47 is diagonal in the basis in which @xmath16 is diagonal . assuming @xmath523 and @xmath524 , there are two family of states that control the low - energy spectrum . first , we should consider the magnetized states @xmath525 then , we should also consider the @xmath526 states ( we call them kink states @xcite ) @xmath527 if @xmath528 is the hamiltonian for @xmath529 we have ( @xmath3 ) @xmath530 the value of @xmath198 determines which of the states is the ground state of the system . for @xmath339 , the ground state is doubly degenerate and spanned by @xmath531 and @xmath532 , while for @xmath533 the @xmath526 kink states are the lowest - energy ones . for @xmath534 , the magnetized and the kink states are degenerate . we wish now to determine how this picture changes when the perturbation @xmath535 is added . let us first understand the behavior for @xmath339 . for @xmath522 the ground state is doubly degenerate . such a degeneracy is lifted by perturbation ( [ hg - def ] ) . note , however , that @xmath536 vanishes for any @xmath537 . therefore , the energy gap is proportional to @xmath538 , in agreement with the exact results obtained in sec . [ sec5 ] . let us now consider the corrections to the ground state which are proportional to @xmath303 . at order @xmath303 , @xmath531 mixes with the states generated by @xmath539 , so that we can write @xmath540 where @xmath541 , @xmath542 , @xmath201 are of order @xmath303 and state @xmath543 is defined so that @xmath544 for @xmath545 and @xmath546 . requiring @xmath547 to be normalized , we immediately obtain @xmath548 . the coefficients @xmath541 , @xmath542 , and @xmath201 are fixed by the eigenvalue condition @xmath549 we obtain @xmath550 the analysis for @xmath532 is analogous and leads to @xmath551 . once the degeneracy is lifted , the ground state should be an eigenstate of @xmath552 , hence it should be of the form @xmath553 we are now in the position to compute correlation functions on the ground state . for the local magnetization , we obtain @xmath554 this result shows the the ground state does not show a local magnetization , except at the boundaries . for @xmath555 , the local magnetization @xmath297 vanishes as in the zero - field case . the applied field @xmath198 is too small to destabilize the low - energy state . we can analogously compute the correlation function @xmath298 : @xmath556 where @xmath557 . except on the boundaries , the correlation function is equal to the square of the bulk magnetization @xmath302 , as in the zero - field case . let us now consider the large - field region in which the ground state is @xmath526 degenerate and the relevant states are the kink states ( [ kink - states ] ) . the degeneracy is lifted when @xmath558 is included . the full hamiltonian restricted to the subspace spanned by the kink states has the form ( we write it for @xmath559 ) @xmath560 where @xmath561 . the eigenvalues and eigenvectors of this matrix have been computed in [ app.a1 ] . there are @xmath562 energy states @xmath563 with @xmath564 . the ground state corresponds to @xmath565 . the corresponding eigenfunctions are @xmath566 with @xmath567 . using this expression , we can compute the average of @xmath16 . for the average on state @xmath116 , we obtain @xmath568 where @xmath569 . if we now consider the ground state , in the large-@xmath0 limit , we can rewrite eq . ( [ magnetization ] ) as @xmath570 if we define @xmath383 and @xmath308 we obtain the more symmetric form @xmath571 analogously , we can compute the correlation function . we obtain @xmath572 where @xmath573 . focusing again on the ground state and taking the limit @xmath260 , if @xmath574 , @xmath575 , we obtain @xmath576      we should now discuss the behavior of the system in the intermediate case in which kink states and ferromagnetic states are degenerate . for @xmath522 this case corresponds to @xmath577 . as discussed in sec . [ sec5 ] , when the @xmath303 dependent term is added , the intermediate case corresponds to @xmath578 . therefore , to discuss the intermediate case , we set @xmath579 and consider the model restricted to the @xmath500 subspace spanned by @xmath580 , @xmath581 ( @xmath582 ) , @xmath583 . if @xmath584 is the energy of the kink states , the hamiltonian can be written as ( here @xmath503 ) @xmath585 the eigenvalues can be computed using the results of [ app.a1 ] . indeed , @xmath47 has the same form as matrix @xmath586 defined in eq . ( [ matrix - atilde - def ] ) . since @xmath587 , the secular equation ( [ recursion - tildean ] ) becomes @xmath588 which implies @xmath589 as the energy of each mode is given by @xmath590 , the ground state is obtained by taking @xmath591 . the ground - state energy is @xmath592 and the corresponding eigenvector is simply @xmath593 it is easy to compute the local magnetization and the correlation function . we find @xmath594 which shows that @xmath297 varies linearly between @xmath595 and @xmath596 . if we set , as usual , @xmath383 and @xmath308 we obtain in the large-@xmath0 limit @xmath597 the two - point correlation function is also easily computed @xmath598      let us now study the crossover behavior . we consider again the same basis as in [ app.b3 ] . the hamiltonian becomes ( here @xmath503 ) @xmath599 the secular equation is obtained by using the results of [ app.a1 ] , eq . ( [ recursion - tildean ] ) . we set @xmath600 , @xmath601 , @xmath602 , and @xmath590 . since @xmath603 , we obtain @xmath604   \sin k ( l+1 ) = 0.\ ] ] in the kink phase in which @xmath365 , @xmath159 varies between @xmath605 for @xmath419 and @xmath606 for @xmath370 . hence , we write @xmath607 . expanding the secular equation to order @xmath446 we obtain @xmath608 which coincides with eq . ( [ sec - etac-1 ] ) . the solutions of this equation either solve @xmath609 or @xmath610 . the ground state corresponds to the lowest value of @xmath381 that is a solution of eq . ( [ gap - appb ] ) , hence it satisfies @xmath609 . note the limiting values : @xmath611 we wish now to compute the ground - state eigenfunction , which we express as @xmath612 in the basis @xmath613 , @xmath614 , defined in [ app.b3 ] . if @xmath615 , the coefficients satisfy the relations @xmath616 the solution of eq .  ( [ rec - alpha - scaling-2 ] ) is @xmath617 the constants @xmath475 and @xmath476 are fixed by condition ( [ rec - alpha - scaling ] ) and by the normalization condition @xmath618 . in the large-@xmath0 limit these two conditions give @xmath619 we can now compute @xmath297 and @xmath298 . for the local magnetization we have @xmath620 in the large-@xmath0 limit , using @xmath621 , we obtain @xmath622 where @xmath383 , @xmath308 . note that the prefactor guarantees that @xmath297 is equal to @xmath623 at the two boundaries . we can also compute the correlation function , which satisfies the relation @xmath624 correspondingly , we obtain @xmath625 and @xmath626 let us finally consider @xmath627 defined in eq . ( [ defden ] ) . for @xmath377 , @xmath628 , so that @xmath629 . using the definition we obtain @xmath630 substituting the explicit expression of @xmath631 and using @xmath609 , we find @xmath632 the above - reported results apply to the kink phase @xmath352 . for @xmath356 the ground state is a localized state with @xmath633 , where @xmath358 satisfies the equation @xmath634 the solutions of this equation satisfy either @xmath635 or @xmath636 . since @xmath637 , the ground state is obtained by considering the largest positive solution of the equation . it is then easy to prove that it satisfies @xmath635 . the solution for @xmath356 is then obtained by analytic continuation of the results obtained for @xmath365 . it is enough to replace @xmath381 with @xmath638 . we therefore obtain @xmath639 note that @xmath640 for @xmath372 , so that @xmath641 in this limit . 99                    g. g. cabrera and r. jullien , universality of finite - size scaling : the role of the boundary conditions , phys . rev . lett . * 57 * , 393 ( 1986 ) ; role of boundary conditions in the finite - size ising model , phys . b * 35 * , 7062 ( 1987 ) .                        a. dutta , u. divakaran , d. sen , b. k. chakrabarti , t. f. rosenbaum , and g. aeppli , quantum phase transitions in transverse field spin models : from statistical physics to quantum information , arxiv:1012.0653v3 ; a. dutta , g. aeppli , b. k. chakrabarti , u. divakaran , t. f. rosenbaum and d. sen , _ quantum phase transitions in transverse field spin models : from statistical physics to quantum information _ ( cambridge university press , cambridge , 2015 ) .
Abstract:

Answer: we present a detailed study of the finite one - dimensional quantum ising chain in a transverse field in the presence of boundary magnetic fields coupled with the order - parameter spin operator . we consider two magnetic fields located at the boundaries of the chain that have the same strength and that are aligned in the same or in the opposite direction . we derive analytic expressions for the gap in all phases for large values of the chain length @xmath0 , as a function of the boundary field strength . we also investigate the behavior of the chain in the quantum ferromagnetic phase for oppositely aligned fields , focusing on the magnet - to - kink transition that occurs at a finite value of the magnetic field strength . at this transition we compute analytically the finite - size crossover functions for the gap , the magnetization profile , the two - point correlation function , and the density of fermionic modes . as the magnet - to - kink transition is equivalent to the wetting transition in two - dimensional classical ising models , our results provide new analytic predictions for the finite - size behavior of ising systems in a strip geometry at this transition .