frequency tail. We find that $|F(q,\omega)|$ is a rapidly decreasing function of the number of particles involved in the excitation. In particular, we find for all $\Delta$ that the peak is completely dominated by two-particle (single particle-hole) and the tail by four-particle states (denoted by 2$p$ and 4$p$ states, respectively). Including up to eight-particle as well as bound states we verify using Density Matrix Renormalization Group (DMRG) that the sum rules are fulfilled with high accuracy corroborating our numerical results. By solving the BA equations for small $\Delta$ and infinite system size analytically we show that the width of the peak scales like $q^2$ for $h\neq 0$. Furthermore, we calculate the high-frequency tail analytically based on a parameter-free effective bosonic Hamiltonian. We demonstrate that our analytical results for the linewidth and the tail are in excellent agreement with our numerical data.

For a chain of length $N$ the longitudinal dynamical structure factor is defined by $$\begin{aligned}
S^{zz}\left(q,\omega\right)&=&\frac{1}{N}\sum_{j,j^{\prime}=1}^{N}e^{-iq\left(j-j^{\prime}\right)}\int_{-\infty}^{+\infty}\!\!\!\!\!\!\!\!
dt\,
e^{i\omega t}\left\langle
  S_{j}^{z}\left(t\right)S_{j^{\prime}}^{z}\left(0\right)\right\rangle
\nonumber \\
&=& \frac{2\pi}{N}\sum_{\alpha}\left|\left\langle
    0\left|S_{q}^{z}\right|\alpha\right\rangle
\right|^{2}\delta\left(\omega-E_{\alpha}\right) \;.
\label{strucFac}\end{aligned}$$ Here $S_{q}^{z}=\sum_{j}S_{j}^{z}e^{-iqj}$ and $\left|\alpha\right\rangle $ is an eigenstate with energy $E_{\alpha}$ above the ground state energy. For a finite system, $S^{zz}\left(q,\omega\right)$ at fixed $q$ is a sum of $\delta$-peaks at the energies of the eigenstates. In the thermodynamic limit $N\rightarrow\infty$, the spectrum is continuous and $S^{zz}\left(
frequency tail. We find that $ |F(q,\omega)|$ is a rapidly decreasing affair of the numeral of particles involved in the excitation. In finical, we find for all $ \Delta$ that the peak is wholly dominated by two - particle (single atom - fix) and the tail by four - atom states (denoted by 2$p$ and 4$p$ states, respectively). Including up to eight - atom as well as bound state we verify using Density Matrix Renormalization Group (DMRG) that the kernel rules are fulfilled with high accuracy corroborate our numerical results. By solving the BA equation for small $ \Delta$ and infinite system size analytically we show that the width of the peak scale like $ q^2 $ for $ h\neq 0$. Furthermore, we calculate the high - frequency tail analytically based on a parameter - free effective bosonic Hamiltonian. We demonstrate that our analytical results for the linewidth and the tail are in excellent agreement with our numerical data. 

 For a chain of length $ N$ the longitudinal dynamical structure factor is specify by $ $ \begin{aligned } 
 S^{zz}\left(q,\omega\right)&=&\frac{1}{N}\sum_{j, j^{\prime}=1}^{N}e^{-iq\left(j - j^{\prime}\right)}\int_{-\infty}^{+\infty}\!\!\!\!\!\!\!\! 
 dt\, 
 e^{i\omega t}\left\langle 
   S_{j}^{z}\left(t\right)S_{j^{\prime}}^{z}\left(0\right)\right\rangle 
 \nonumber \\ 
 & = & \frac{2\pi}{N}\sum_{\alpha}\left|\left\langle 
     0\left|S_{q}^{z}\right|\alpha\right\rangle 
 \right|^{2}\delta\left(\omega - E_{\alpha}\right) \; . 
 \label{strucFac}\end{aligned}$$ Here $ S_{q}^{z}=\sum_{j}S_{j}^{z}e^{-iqj}$ and $ \left|\alpha\right\rangle $ is an eigenstate with energy $ E_{\alpha}$ above the ground department of state energy. For a finite system, $ S^{zz}\left(q,\omega\right)$ at fixed $ q$ is a sum of $ \delta$-peaks at the department of energy of the eigenstates. In the thermodynamic limit $ N\rightarrow\infty$, the spectrum is continuous and $ S^{zz}\left (
freeuency tail. We find that $|F(q,\omega)|$ is a rcpidly vecreasjng funcgion of the number of particpew invilved in the excitatiov. In partpcular, we finv for all $\Delta$ vgat the peak jd cok'letely dominatgd by two-pardicle (single pdrgiele-hole) and the tail by four-particle states (dfnoted by 2$p$ anq 4$p$ xeatea, respectively). Including up to eiggt-partible as well as boind states we verify using Dendity Matrix Renormwlization Gtkup (EMRG) that thd sum rules are fulfilmed with high accuracy corroborxting our numerucql ggsults. By sooving the BA equabpons fos small $\Delta$ and inflnite sywtem size analyticallb we show that the wydth of tve peak scales like $q^2$ for $v\neq 0$. Fufrhefmode, ws calcklave the high-rrequency tqil analytically baxeq on a parametsr-free esfective bosonic Hamiltonian. We demonstgate that our analytical rewults for the linewidjh and the tail are in excellent agreement with our numericdl dava.

Wor a zyaln of length $N$ the longitudinal dynamical strtdtirv factor is defincd by $$\begin{aligned}
X^{zx}\lrst(q,\omega\right)&=&\wrac{1}{N}\sbj_{j,n^{\prime}=1}^{N}e^{-iq\left(j-j^{\prlme}\righj)}\int_{-\indty}^{+\infty}\!\!\!\!\!\!\!\!
du\,
e^{i\omrga t}\left\langle
  S_{j}^{z}\left(t\ritht)S_{j^{\prime}}^{z}\lvft(0\rught)\right\rangle
\nonbmber \\
&=& \frac{2\pn}{N}\sum_{\akpha}\lrft|\left\langle
    0\left|S_{q}^{v}\right|\zlpha\right\rwngle
\righf|^{2}\aelta\left(\omega-E_{\auphs}\rhght) \;.
\labtu{strucFac}\end{alignqd}$$ Here $S_{w}^{z}=\suk_{j}S_{j}^{z}e^{-idj}$ amd $\lefe|\alpha\righh\rangle $ is an eigenstatf witk enesgy $E_{\alpha}$ above the ground state energy. Hpr a finite xyvtek, $S^{zz}\lefc(q,\omegs\right)$ at fixqd $q$ is a sum pf $\deltc$-peaks at the envrgies of the eigenseates. In the djermodynamic limit $N\wighrarriw\infty$, ghe spectrum ix continuous and $S^{zz}\left(
frequency tail. We find that $|F(q,\omega)|$ is decreasing of the of particles involved we for all $\Delta$ the peak is dominated by two-particle (single particle-hole) and tail by four-particle states (denoted by 2$p$ and 4$p$ states, respectively). Including up eight-particle as well as bound states we verify using Density Matrix Renormalization Group that sum are with high accuracy corroborating our numerical results. By solving the BA equations for small $\Delta$ and system size analytically we show that the width the peak scales like for $h\neq 0$. Furthermore, we the tail analytically on parameter-free bosonic Hamiltonian. We that our analytical results for the linewidth and the tail are in excellent agreement with our numerical For a length $N$ longitudinal structure is defined by dt\, e^{i\omega t}\left\langle S_{j}^{z}\left(t\right)S_{j^{\prime}}^{z}\left(0\right)\right\rangle \nonumber \\ \right|^{2}\delta\left(\omega-E_{\alpha}\right) \;. \label{strucFac}\end{aligned}$$ Here $S_{q}^{z}=\sum_{j}S_{j}^{z}e^{-iqj}$ and $\left|\alpha\right\rangle $ an eigenstate energy $E_{\alpha}$ above the ground state For a finite system, $S^{zz}\left(q,\omega\right)$ at fixed $q$ a sum of $\delta$-peaks at the energies of the eigenstates. In the thermodynamic limit $N\rightarrow\infty$, is continuous and $S^{zz}\left(
frequency tail. We find that $|F(q,\Omega)|$ is a raPidly DecReaSiNg fuNctiOn of the number oF PartIcles involved in the exciTatioN. IN PartICuLar, we Find for ALl $\dELta$ ThAt The PeAK iS compLetEly domiNated by two-ParTiCle (single parTIcLe-hole) and tHe tAil by four-parTicLe statEs (DenOTed by 2$P$ anD 4$p$ staTes, resPEctiveLy). IncludiNg UP to eigHT-particLE As Well As bound states we veRIfY Using Density MaTrix ReNoRMaLIZatIon group (DMRG) tHaT the sUM rules aRE fULFIllED with high accuRacy corroboRAtiNg our nUmEriCAl resuLts. By SoLVinG the BA equatIons For small $\DElta$ anD InfinitE System sIze anaLytIcaLly wE ShOw ThaT tHE wiDTh Of tHE peAk scales LiKe $Q^2$ for $h\Neq 0$. FURTHErmoRe, wE calCulatE the high-frequEncY taiL AnaLyticAlly bAsed On A paraMeter-fRee efFeCtive bosonic HamIltoNian. We demOnsTrAte ThAt our ANalytiCal ResUlts for The lineWIdtH aND THe Tail are in excellent AgREEmEnt with oUr numeRIcAl DAta.

For a cHaIn oF lenGTH $N$ the LongITuDinal dynAmical STrUcTure facToR is defInEd bY $$\beGin{alIGned}
s^{zz}\lefT(q,\omega\rIght)&=&\fRAc{1}{N}\sum_{j,j^{\prime}=1}^{N}E^{-Iq\left(j-j^{\prime}\RIgHT)}\InT_{-\InftY}^{+\inFty}\!\!\!\!\!\!\!\!
dt\,
e^{i\omegA t}\leFT\lanGle
  S_{J}^{Z}\lEft(T\Right)s_{j^{\priMe}}^{Z}\LeFT(0\right)\right\rangle
\noNuMber \\
&=& \frAc{2\pi}{N}\Sum_{\alpha}\left|\lEft\langle
    0\lEFT|s_{q}^{z}\right|\AlphA\RiGHt\rangle
\right|^{2}\dElta\lEft(\omega-E_{\aLPha}\right) \;.
\Label{StrucFac}\End{aligneD}$$ hEre $S_{q}^{z}=\suM_{j}S_{J}^{z}e^{-Iqj}$ And $\LEFt|\Alpha\right\ranGLE $ is aN eIgenstaTe wIth enerGy $E_{\AlpHa}$ aBovE tHe ground sTate enerGy. foR a FiNitE systEM, $S^{zz}\left(Q,\oMegA\rIghT)$ at fiXEd $q$ is a Sum of $\DeltA$-pEaKS at The enerGIeS OF the EiGeNstaTes. in The thErmoDYnaMic limiT $N\rightarRow\INfty$, ThE sPectrum Is continuous aNd $s^{zz}\left(
frequency tail. We find th at $|F(q,\ omega )|$  is  a  rap idly  decreasing fu n ctio n of the number of par ticle si nvol v ed  in t he exci t at i o n.In  p art ic u la r, we  fi nd forall $\Delt a$th at the peaki scompletely  do minated by t wo- partic le  (s i nglepar ticle -hole) and th e tail by  f o ur-par t icle st a t es  (de noted by 2$p$ and 4$ p $ states, resp ective ly ) .I n clu din g up to ei gh t-par t icle as we l l asb ound states w e verify us i ngDensit yMat r ix Ren ormal iz a tio n Group (DM RG)that thesum ru l es aref ulfille d with  hi ghaccu r ac ycor ro b ora t in g o u r n umerical  r es ults.  Bys o l v ingthe  BAequat ions for smal l $ \Del t a$and i nfini te s ys tem s ize an alyti ca lly we show tha t th e width o f t he  pe ak  scal e s like  $q ^2$  for $h \neq 0$ .  Fu rt h e r mo re, we calculate t he h ig h-freque ncy ta i lan a lyticall ybas ed o n a par amet e r- free eff ective bo so nic Ham il tonian .Wedem onstr a te t hat ou r analyt icalr esults for the linewidth and th e ta i l ar e i n excellent  agr e emen t wi t hour numer icalda t a. 
For a chain of len gt h $N$the l ongitudinal d ynamical s t r u cture fa ctor is defined by $$\ begin {aligned}S ^{zz}\le ft(q, \omega\r ight)&=&\ f r ac{1}{N} \su m_{ j,j ^{\ p r im e}=1}^{N}e^{- i q \lef t( j-j^{\p rim e}\righ t)} \in t_{ -\i nf ty}^{+\in fty}\!\! \! \! \! \! \!\ !
dt\ , 
e^{i\om eg a t }\ lef t\lan g le
  S _{j}^ {z}\ le ft ( t\r ight)S_ { j^ { \ prim e} }^ {z}\ lef t( 0\rig ht)\ r igh t\rangl e
\nonumb er\ \
&= &\f rac{2\p i}{N}\sum_{\a lp ha}\left|\ le ft\ langle    0\lef t|S_{q}^{z}\right|\alph a \right\ ran gle
\ righ t|^{2}\de lta \left( \om e ga-E_{ \alpha }\rig ht ) \ ; . 
\lab e l {s tru cF ac}\end{al i g ned }$$ H er e $S _{q}^{z }=\sum_{j}S_{j}^{z } e^{ -iqj}$ and $\ lef t|\a l p ha \ri g ht \ ran gl e  $i s  an eigenstatewith energ y$ E_ {\alpha}$a bov ethe gro und sta te en e rgy. Fo r a finit e system,  $ S^{z z } \le ft(q,\omeg a\right) $ at fixe d  $q$i sa sum  of  $\del ta $-p eaksat the ene rgies  of th eeigens tates .In the t hermodynamic limit $N\r ightar row\i nft y$, the s pec t rum  is conti nuou s and $S^{ zz} \le ft(
frequency tail._We find_that $|F(q,\omega)|$ is a_rapidly decreasing_function_of the_number_of particles involved_in the excitation._In particular, we find_for all $\Delta$_that_the peak is completely dominated by two-particle (single particle-hole) and the tail by four-particle_states_(denoted by_2$p$_and_4$p$ states, respectively). Including up_to eight-particle as well as_bound states_we verify using Density Matrix Renormalization Group (DMRG)_that_the sum rules_are fulfilled with high accuracy corroborating our numerical results._By solving the BA equations for_small $\Delta$ and_infinite_system_size analytically we show_that the width of the peak_scales like $q^2$ for $h\neq 0$._Furthermore, we calculate the high-frequency tail analytically_based on a parameter-free effective bosonic_Hamiltonian. We demonstrate that our_analytical results_for the linewidth and the_tail are in_excellent agreement_with our numerical_data.

For a chain of length $N$_the longitudinal dynamical_structure factor is defined by $$\begin{aligned}
S^{zz}\left(q,\omega\right)&=&\frac{1}{N}\sum_{j,j^{\prime}=1}^{N}e^{-iq\left(j-j^{\prime}\right)}\int_{-\infty}^{+\infty}\!\!\!\!\!\!\!\!
dt\,
e^{i\omega_t}\left\langle
_ S_{j}^{z}\left(t\right)S_{j^{\prime}}^{z}\left(0\right)\right\rangle
\nonumber \\
&=&_\frac{2\pi}{N}\sum_{\alpha}\left|\left\langle
__ _0\left|S_{q}^{z}\right|\alpha\right\rangle
\right|^{2}\delta\left(\omega-E_{\alpha}\right) \;.
\label{strucFac}\end{aligned}$$ Here_$S_{q}^{z}=\sum_{j}S_{j}^{z}e^{-iqj}$_and $\left|\alpha\right\rangle_$_is an eigenstate with energy $E_{\alpha}$_above_the ground state energy. For a finite_system, $S^{zz}\left(q,\omega\right)$ at fixed_$q$_is a sum of_$\delta$-peaks at the energies of_the eigenstates. In the thermodynamic limit_$N\rightarrow\infty$, the_spectrum is_continuous and $S^{zz}\left(