Formation probabilities and Shannon information and their time evolution after quantum quench in transversefield XYchain
Abstract
We first provide a formula to calculate the probability of occurrence of different configurations (formation probabilities) in a generic free fermion system. We then study the scaling of these probabilities with respect to the size in the case of critical transversefield XYchain in the bases. In the case of the transverse field Ising model, we show that all the ”crystal” configurations follow the formulas expected from conformal field theory (CFT). In the case of critical chain, we show that the only configurations that follow the formulas of the CFT are the ones which respect the filling factor of the system. By repeating all the calculations in the presence of open and periodic boundary conditions we find further support to our classification of different configurations. Using the developed technique, we also study Shannon information of a subregion in our system. In this respect we distinguish particular configurations that are more important in the study of the scaling limit of the Shannon information of the subsystem. Finally, we study the evolution of formation probabilities, Shannon information and Shannon mutual information after a quantum quench in free fermion system. In particular, for the intial state considered in this paper, we demonstrate that the Shannon information after quantum quench first increases with the time and then saturates at time , where is the size of the subsystem.
pacs:
03.67.Mn,11.25.Hf, 05.70.JkI Introduction
Studying correlation functions in manybody systems has been considered one of the main topics in statistical mechanics and condensed matter physics for many years. Although, for a longtime the main quantities of interest were the correlation functions of local observables the recent interest in calculating nonlocal quantities, especially the entanglement entropy, has made significant changes. One of the main reasons for this interest (at least in dimensional critical systems) is that by calculating some of the nonlocal observables one can derive the central charge of the system without referring to the velocity of sound, for the case of entanglement entropy see CC2009 (). Another nonlocal quantity which has been studied for many years with Bethe ansatz techniques and some other methods is emptiness formation probability Essler (); U1 (); Shiroishi (); AbanovKorepin; Franchini (); Stephan2013 (). In the case of spin chains, it is the probability of finding a sequence of up spins in the system (note that almost all of the studies in this regard concentrate on the bases). These studies show that this probability, with respect to the sequence size, decreases like a Gaussian in the case of systems with symmetry U1 () and exponentially in other cases Shiroishi (); Franchini (). In the critical cases, the Gaussian and exponential are accompanied with a powerlaw decay with a universal exponent. In a recent development Stephan2013 () it was shown that for those critical systems without symmetry this universal exponent is dependent on the central charge of the system. The argument is based on connecting the configuration of all spins up to some sort of boundary conformal field theory. One should notice that the argument works just for those bases that can be connected to boundary conformal field theory.
In an apparently connected studies recently many authors investigated the Shannon information of quantum systems in different systems Stephan2009 (); Stephan2010 (); Oshikawa2010 (); Stephan2011 (); Alet2013 (); Alet2014a (); Alet2014b (); Alet2014c (); Trebst2014 (); Trebst2015 (); Jose (). The Shannon information is defined as follows: Consider the normalized ground state eigenfunction of a quantum spin chain Hamiltonian , expressed in a particular local bases , where are the eigenvalues of some local operators defined on the lattice sites. The Shannon entropy is defined as
(1) 
where is the probability of finding the system in the particular configuration given by . As it is quite clear this quantity is bases dependent and to calculate it one needs in principle to know the probability of occurrence of all the configurations. The number of all configurations increases exponentially with the size of the system and that makes analytical and numerical calculations of this quantity quite difficult. Note that the emptiness formation probability is just one of the whole possible configurations. The Shannon information of the system changes like a volume law with respect to the size of the system so in principle one does not expect to extract any universal information by studying the leading term. The universal quantities should come from the subleading terms. To study subleading terms, it is useful to define yet another quantity called Shannon mutual information. By considering local bases it is always possible to decompose the configurations as a combination of the configurations inside and outside of a subregion as . Then one can define the marginal probabilities as and for the subregion and its complement . Then the Shannon mutual information is
(2) 
where and are the Shannon informations of the subregions and . From now on instead of using we will use just . Since in the above quantity the volume part of the Shannon entropy disappears Shannon mutual information provides a useful technique to study the subleading terms. Note that one can similarly define the above quantity also for two regions and , i.e. that are not necessarily complement of each other. The Shannon mutual information has been studied in many classical Wolf2008 (); wilms (); Grassberger (); Rittenberg2014 () and quantum systems Um (); AR2013 (); Stephan2014 (); AR2014 (); AR2015 (). Numerical studies on variety of different periodic quantum critical spin chains show that for particular bases (so called conformal bases) we have AR2013 (); AR2014 (); AR2015 ()
(3) 
where and are the total size and subsystem size respectively and is very close to , with the central charge of the system. Note that if one takes an arbitrary base (nonconformal bases) the coefficient is nothing to do with the central charge. It is worth mentioning that based on AR2014 () the conformal bases are those bases that can be connected to some sort of boundary conformal field theory in the sense of Stephan2013 (). For example in the transverse field Ising model the and bases are the conformal bases. Note that if one considers the Shannon information of the subsystem we will have . Since to extract the Shannon information one needs to use all the probabilities the only way to consider all of them in the numerical calculations is exact diagonalization. This makes the numerical calculation for large sizes very difficult. The results of the articles AR2013 (); AR2014 (); AR2015 () are all for periodic systems with whenever there are spin onehalf system and smaller sizes for systems with bigger spins. In a recent work Stephan2014 () the author was able to consider an infinite system and study the Shannon information of the subsystem up to the size . It was concluded that for the chain is with but for the Ising model although it is very close to with the results are suggesting that probably is not exactly connected to the central charge. Notice that all of the above calculations are done by considering periodic boundary conditions for the connected regions and . It is not clear how the equation (3) might change if one considers open boundary conditions. Finally it is worth mentioning that some of the the above results have recently been extended to disconnected regions in Rajabpour2015d ().
Motivated by the studies of emptiness formation probability and Shannon information of the subsystem we study here some related quantities. First of all, as it is natural one might be interested in studying the scaling limit of some other configurations with respect to the size of the subsystem. For example, consider an antiferromagnetic configuration in the Ising model or any other configuration with pattern. It is very important to know that these configurations are also flowing to some sort of boundary conformal field theory or not. This study will clearly also help to understand the nature of the Shannon mutual information. In addition, this kind of studies also are very useful in the calculations of postmeasurement entanglement entropy and localizable entanglement entropy Rajabpour2015b (); Rajabpour2015c (). Having the above motivations in mind, we study formation probabilities , Shannon information and their evolution after a quantum quench in the quantum chain in the bases.
The outline of the paper is as follows: In the next section we will first define our system of interest, i.e. chain and then we will provide a method to calculate the probability of any configuration in the free fermionic systems. In section three we will list all the known analytical results regarding emptiness formation probability for infinite systems and also finite systems with periodic and open boundary conditions. Explicit distinguishment is made between critical Ising model and chain with symmetry. In section IV, we will define many different configurations with specially defined pattern and calculate their corresponding probabilities numerically. Here again, we discuss Ising and universalities separately for infinite systems and for the systems with boundary. We also discuss configurations that do not have any pattern. In section V, we study the Shannon information in the transverse field critical Ising model and also critical chain. We will classify the configurations based on their magnetization and show that in principle just a small part of the configurations can have a finite contribution to Shannon information in the scaling limit. Section VI is devoted to the evolution of formation probabilities and Shannon information after a quantum quench. We prepare the system in a particular state and then we let it evolve with another Hamiltonian and study the time evolution of the formation probabilities and especially Shannon information. Finally, the last section is about our conclusions and possible future works.
Ii Formation probabilities from reduced density matrix
In this section, we first define the system of interest and after that using the reduced density matrix of this system we will find a very efficient method to calculate formation probabilities for systems that can be mapped to free fermions. The Hamiltonian of the XYchain is as follows:
(4) 
After using JordanWigner transformation, i.e. and with and for open and periodic boundary conditions respectively the Hamiltonian will have the following form:
(5) 
The above Hamiltonian has a very rich phase space with different critical regions Mccoy1971 (). In figure 1 we show different critical regions of the system. To calculate probability of formations for different patterns we first write the reduced density matrix of a block of spins by using block Green matrices. Following Peschela (); Barthel2008 () we first define the operators
(6) 
The block Green matrix is defined as
(7) 
The elements of the Green matrix can be calculated following Lieb () and we will mention their explicit form for open and periodic boundary conditions later.
To calculate the reduced density matrix after partial measurement we need to first define fermionic coherent states. They can be defined as follows
(8) 
where ’s are Grassmann numbers following the properties: and . Then it is easy to show that
(9) 
Using the coherent states (8) the reduced density matrix has the following form Barthel2008 ()
(10)  
where . One can use the above formula to extract the formation probability for arbitrary configuration in the bases as follows: first of all to extract the probability of particular configuration we need to look to the diagonal elements of . When the spin in the direction is up in the fermionic representation it can be understood as the lack of a fermion in that site which in the language of coherent states means that the corresponding is zero in the equation (10). After putting some of the ’s equal to zero one will have a new reduced density matrix in the coherent state basis with this constraint that some of the spins are fixed to be up. In other words in the equation (10) instead of we will have which is a submatrix of the matrix . The elements of the new reduced density matrix will be , where we put the ’s corresponding to the sites filled with fermions equal to . To extract the probability of formation one just needs to integrate over all the ’s that correspond to the down spins. In other words after using formulas of the Grassmann Gaussian integrals the formation probability will have the following formula
(11) 
where is the minor of the matrix corresponding to the configuration . Notice that we just need principal minors of the matrix . Since the sum of all the principal minors of the matrix is equal to the normalization is ensured. We have number of rankk minors for matrix with size . Summing over the number of all principal minors, one can obtain which is the number of all possible configurations. Configurations with the same minor rank have the same number of up spins and by increasing we actually increase the number of down spins in the corresponding configuration. For example is the case with all spins up (down) and sometimes it is called emptiness formation. The above formula gives a very efficient way to calculate the formation probabilities in numerical approach. Since, as we will show, all of these probabilities are exponentially small with respect to the size of the subsystem it is much easier to work with the logarithm of them and define logarithmic formation probabilities
(12) 
All of the calculations done in this paper are based on the formula (11) and are performed using Mathematica. In the next sections, we will study the formation probabilities for different configurations with the crystal order (pattern formation probabilities) with respect to the size of the subsystem for some particular critical regions of the system.
Iii Emptiness formation probability: known results
Before presenting our results, we first review here the wellknown facts regarding emptiness formation probability ( and ). For reasons that will be clear in the next section, we will call these two configurations and respectively. Using FisherHartwig theorem the emptiness formation probability was already exhaustively studied in Franchini (). The results were generalized to arbitrary conformal critical systems in Stephan2013 (). Due to the symmetry there is an important difference between the XX critical line and the Ising critical point. For this reason, we report the results regarding these two possibilities separately.
iii.1 Critical Ising point
We list here the results regarding the logarithmic emptiness formation probability at the critical Ising point ().
Configuration , i.e. :
This configuration corresponds to and has the highest probability and using the equation (11) one can easily show that
(13) 
The above formula is valid independent of the boundary conditions and the size of the system. For the infinite system at the critical Ising point the matrix has the following form:
(14) 
Since the above matrix is a Toeplitz matrix using FisherHartwig conjecture in Franchini () it was shown that the logarithmic probability for this configuration changes with the subsystem size as follows:
(15) 
where , and and are some nonuniversal numbers. At the critical point of the Ising model these numbers are known , with the Catalan constant and . The term is the result of the paper Stephan2013 (). Using conformal field theory techniques in Stephan2013 () it was argued that for generic critical systems the coefficient of the logarithm should be , where is the central charge of the critical system. For the Ising universality class .
When the size of the total system is finite depending on the form of the boundary conditions, periodic or open; has the following two forms
(16)  
Notice that for the first equation reduces to (14) and the second equation will give the result for a semiinfinite chain. The results for emptiness formation probabilities for the above two cases are Stephan2013 ()
(18)  
(19) 
where and . It is also conjectured that and . The above two equations are derived using boundary conformal field theory techniques and in principle they are valid in the Ising case because configuration is related to the free conformal boundary condition in the conformal Ising model.
Configuration , i.e. :
This case, which is also studied in Franchini () and Stephan2013 (), corresponds to and has the lowest probability. One can easily show that
(20) 
For the infinite system it follows similar formula as (15) with also an extra term, in other words,
(21) 
where . At the critical Ising point and and is unknown. The oscillating term is mathematically explained by using generalized FisherHartwig conjecture in Franchini (). To the best of our knowledge its presence at the critical point has not been understood by physical arguments footnote1 (). Our numerical results in the next section will show that the term is present whenever the parity of the number of down spins changes with the subsystem size. The term is very important to be considered in numerical calculations to get reliable results for which is the universal and the most interesting term.
When the size of the system is finite depending on the type of the boundary conditions boundary changing operators can play an important role. The following formulas are presented in Stephan2013 ():
(22)  
where , and with being the conformal weight of the boundary changing operator. The dots are the subleading terms.
iii.2 XX critical line
The critical XX chain has symmetry which as it is already discussed extensively in the literature is the main reason for having Gaussian decaying emptiness formation probability U1 (); Stephan2013 (). Since in this model the equation (10) has simpler form
(24) 
where and . Finally we have
(25) 
The form of the matrices in the periodic and open cases are Fagotti ():
(26)  
(27) 
where is the Fermi momentum and with as the closest integer larger (smaller) than .
The all spins up and down configurations do not lead to conformal boundary conditions and so none of the equations that we mentioned in the last subsection are valid. However, using Widom conjecture it is already known that, see for example Shiroishi (), the probabilities for both and show Gaussian behavior. For systems with symmetry one expects the following behavior for logarithmic emptiness formation probability U1 ():
(28) 
where for critical XX chain.
Iv Logarithmic pattern formation probabilities
In this section, we study the logarithmic pattern formation probability defined as with respect to the size of the subsystem. The easiest configurations to study are those that have some kind of crystal structure. Although everything is already known and checked numerically for the emptiness formation probabilities we will also report the results concerning these cases as benchmarks. Here we introduce the configurations that we studied numerically. None of these configurations have been considered before in the literature.
We study here the configurations with with crystal pattern and we call them configurations and for some ranks we will study the two most basic cases. For example we will study

Configurations with :
 a

 b


Configurations with :
 a

 b

All the configurations with the same belong to the cases with an equal rank of the minor in the equation (11). Note that in all of the upcoming numerical calculations in every step we increase the size of the subsystem with a number which is devidable to the length of the base of the corresponding configuration. For some particular ’s the a and b configurations differ by the parity effect. For example, in depending on or with the subsystem has even or odd number of down spins. This means that the parity of the number of down spins changes with the subsystem size for this configuration. However, for this is not the case because in order to have ”perfect crystal” in the subsystem we need to consider a subsystem with with which has always even number of down spins inside. Because of this difference in parity effect for and we expect different subleading behavior for these two cases. Finally notice that one can simply define configurations like and by simply taking bigger bases for the crystals. For example, can be understood as a configuration with the base: three down spins and then three up spins.
iv.1 Transversefield Ising chain
Using the equation (11) and (14) we first studied the crystal configurations introduced in the previous subsection for the case of infinite chain. To calculate the formation probability for every configuration we first use the matrix introduced in (14) to find the matrix . Then for every configuration we use an appropriate minor to calculate the corresponding probability in (11). For example, in the case of this can be done by just finding a minor of which can be derived by calculating the determinant of a submatrix obtained from by removing every other row and column. The results for and (the most interesting quantities in this study) are shown in the Table I. Based on the numerical results one can derive the following conclusions regarding crystal configurations:

Whenever the parity of the number of down spins in a configuration changes with respect to the size of the subsystem we have the oscillating term . For example, in the case of we have the subleading term but correction is absent in . It appears again for in the form of . Generalization to other configurations is strightforward. .

Although in some cases for bigger is smaller than with smaller in average increases with .
Configuration  

We then studied the same configurations for the periodic boundary condition. In the figure 2 it is shown that all of the configurations follow the formula (18) . The case of the open boundary condition is more tricky and depending on the configuration we have two possibilities: When the parity of the number of down spins is independent of the size of the subsystem (for example in the configurations and ) we have the formula (19) but when we have the possibility of having odd or even number of down spins in a configurations (for example ) we have the formula (III.1). The results are shown in the Figure 3. This behavior could be anticipated based on the difference between configurations and that we discussed before. It looks like that the boundary changing operator plays a role whenever there is the parity effect in the configuration. Looking to the problem in the language of Euclidean two dimensional classical system one can argue that in the case of open boundary condition we have a strip with a slit on it Stephan2013 (). However, the boundary conditions on the boundary of strip can be different from the boundary condition on the slit, consequently, one needs to consider boundary changing operator on the point where the boundary condition changes. However, in general it is not clear which configurations lead to different boundary conditions on the slit and on the boundary of the strip. Our numerical results indeed give a hint that depending on the bahavior of the parity of the number of down spins in a configuration the conformal boundary condition on the slit can be different. In the next two subsections, we will first comment on the validity of the above results in other cases such as noncrystal configurations. Then we will also indicate the possible universal behavior of our results.
iv.1.1 Logarithmic formation probability of noncrystal configurations
The number of crystal configurations is much smaller than the number of the whole configurations. In fact, the number of crystal configurations grows polynomially with the subsystem size but the number of whole configurations grows exponentially. However, numerically it is very simple to check the formula for many configurations that have a small deviation from the crystal states. For example, one can consider the case with all spins up except one and calculate the logarithmic formation probability using the equation (11). It is clear that one does not expect the result be any different from the equation (15) and indeed numerical results confirm this expectation. The important conclusion of this numerical exercise is that there are many configurations ”close” to crystal configurations that indeed follow either the equation (15) or equation (21) with all having the same ’s but different ’s.
The above results strongly suggest that all of the crystal and noncrystal configurations discussed in this section are flowing to some sort of conformal boundary conditions in the scaling limit.
iv.1.2 Universality
To check that the above results are the properties of the Ising universality class we also studied the critical XYchain which has also central charge . The Green matrix, in this case, is given by
Our numerical results depicted in the Figure 4 show that the coefficient of the logarithmic term is a universal quantity which means that it has a fixed value on the critical XYline. The coefficient of the linear term changes by varying which indicates its nonuniversal nature.
iv.2 XX chain
We repeated the calculations of the last section for also critical XX chain. The central charge of the system is . The results of logarithmic formation probabilities for different magnetic field are shown in the Table II and Table III. Based on the numerical calculations we conclude the followings:

The configurations with follow the equation (15) with . This means that in the scaling limit most probably all of these configurations flow to some sort of bounday conformal conditions. Note that as far as there is no boundary changing operator in the system the equation (15) is valid for any CFT independent of its structure.

All the other configurations follow the equation (28) with which is different for different configurations.
As we mentioned earlier XX chain has a symmetry which means that the number of particles is conserved. The only configurations that respect this symmetry in the subsystem level are the configurations with . Any injection of the particles into the subsystem changes drastically the formation probability. In the case of this phenomena is already explained in Stephan2013 () based on arctic phenomena in the dimer model. It is quite natural to expect that similar structure is valid for all the configurations with . Note that based on our results the coefficient of the is dependent and strictly speaking is not a universal quantity.
Configuration  

Configuration  

We also studied the finite size effect in this model. The results of the numerical calculations of periodic boundary condition for are shown in the Figure 5. It is shown that all of the configurations with follow the formula (18) with .
We also repeated the same calculations for open boundary conditions. We first performed the calculations for semiinfinite open chain and fitted the results to (15) and extracted the . The coefficient of the logarithm not only depends on the rank of the configuration but also to the configuration itself. It also changes with . We were not able to find any universal feature in this case.
The above results suggest that most probably all of the crystal configurations with in the periodic boundary condition flow to a boundary conformal field theory. In the language of Luttinger liquid, the corresponding boundary condition should be the Dirichlet boundary condition Stephan2013 (). The case of the open boundary condition is intriguing and we leave it as an open problem.
V Shannon information of a subsystem
In this section, we study Shannon information of a subsystem in transversefield Ising model and chain. For both models, the Shannon information is already calculated in Stephan2014 () up to the size which it seems to be the current limit for classical computers. The reason that we are interested in revisiting this quantity is to have a more detailed study of the contribution of different configurations. This will give an interesting insight regarding the possible scaling limit for this quantity.
v.1 Critical Ising
In the last section, we studied many different configurations in the critical Ising model and we found that all of them follow with . The natural expectation is that if we plug this formula in the definition of the Shannon information we get
(29) 
where the dots are the subleading terms. The above formula is consistent with AR2013 (). However, one should be careful that although there are a lot of crystal configurations (polynomial number of them) and ”close” to crystal configurations that are connected to the central charge it is absolutely not clear what is going to happen in the scaling limit. For example we repeated the calculations of Stephan2014 () and realized that extraction of the coefficient of the in the above equation is indeed very difficult, see appendix. Here we show where one should look for the most important configurations. After a bit of inspection and numerical check, one can see that the configuration with the highest probability is the . Although the proof of the above statment doesn’t look straightforward one can understand it qualitatively by starting from the ground state of the Ising model with and approaching to the critical point . The ground state of the Ising model with is made of a configuration with all spins up. When we decrease the transverse magnetic field the other configurations start to appear in the ground state. Although the amplitude of the configuration with all spins up decreases by decreasing it still remains always bigger than the other configurations. Another way to look at this phenomena is by looking at the variation of the expectation value of the Hamiltonian for different configurations . It is easy to see that is minimum for the configuration with all spins up. This simply means that most probably when one construct the ground state of the Ising model using the variational techniques this configuration playes the most important role. The least important configuration is with the lowest probability. This can also be understood with the same heuristic argument as above. For every rank of the minors, as we discussed, we have number of configurations which means that for every for large we have with number of configurations. It is obvious that the number of configurations in every rank should be high enough to compensate the exponential decrease of probabilities. We realized that in every rank the configurations a has the lowest probabilities. One can again understand this fact using the variational argument. In this case it is much better to make first the canonical transformation: and in the Hamiltonian of the crtical Ising model. Then one can simply argue that is big if there are a lot of domain walls, i.e. in the system which is the case for the configurations a. Other important configurations are the configurations which divide the subsystem to two connected regions with in one part all the spins are up and in the other part all the spins are down. These configurations are interesting because they have the biggest probabilities among all the configurations corresponding to their minor rank. Note that in this case we have just one domain wall. It is not difficult to see that the probability of all of these configurations decay exponentially with the following coefficient :
(30) 
where is the Catalan constant. In the two extreme points, we recover the previous results. We also checked the validity of the above formula numerically. Having the biggest and smallest probabilities for every rank, we can now easily read the most important ranks. In Figure 6 we depicted the and for different configurations. We also depicted the graph of the number of configurations in every rank. The Figure clearly show that the configurations with can not have any significant contribution in the scaling limit because the number of configurations is not enough to compensate the exponential decay of the probabilities. A similar story seems to be valid also for the values of close to zero. The reason is that the number of configurations with small is such low that can not compensate exponential decay of the probabilities in this region to have a significant contribution in the Shannon information. Just the region between the points that the two lines cut each other will most likely survive in the scaling limit.
The numerical results indeed prove our expectation. In Figure 7 we depicted the contribution of every rank in Shannon information for two different sizes. As it is quite clear the most important contributions come from . The contribution of the configurations with is exponentially small. This means that ignoring a lot of configurations will produce a very small amount of error in the final result of Shannon information. To quantify this argument we calculated the amount of error in the evaluation of the Shannon information if we just keep the configurations with ranks up to . Suppose is the contribution of the configurations with all ranks equal or smaller than . Then the error of truncation can be calculated by . Interestingly we found that the logarithm of the error function is a linear function of , see Figure 8. In other words
(31) 
where is equal to zero and infinity for and respectively. for the other values are shown in the inset of the Figure 8. The above formula shows that one can calculate Shannon information with a controllable accuracy by ignoring nonimportant configurations. Although the above truncation method help to calculate the Shannon information with good accuracy (especially the coefficient of the linear term ) it is still not good enough to calculate the coefficient of the logarithm with controllable precision.
v.2 XX chain
The Shannon information of the subsystem in the XXchain is already discussed in Stephan2014 () and based on numerical results it is concluded that the equation (30) is valid with which is consistent with the conjecture in AR2013 (). Here we just comment on the contribution of different ranks which shows very different behavior from the transverse field Ising chain case. First of all, as we discussed in the previous section when the external field is zero the only configurations that decay exponentially are those that respect the half filling structure of the total system. The rest of the configurations scale like a Gaussian which simply indicates that their contribution is very small in the Shannon information. This is simply because the number of these configurations scale just exponentially. Based on this simple fact one can anticipate that the only configurations that can survive in the scaling limit are those with . Numerical results depicted in the Figure 9 indeed support this idea. Although the is only one among possible minor ranks the number of configurations with this rank is highest with respect to the others which can be one of the reasons that one can obtain a good estimate for the coefficient of the logarithm in (30) with relatively modest sizes.
Vi Evolution of Shannon and mutual information after global quantum quench
Inspired by experimental motivations, the field of quantum nonequilibrium systems has enjoyed a huge boost in the recent decadeexperimental (). One of the interesting directions in this field is the study of information propagation after quantum quench, see for example CCevolution (); Plenio (); LRbound (); fidelity (); Liu (); MM (). Based on semiclassical arguments and also using LiebRobinson bound it is shown CCevolution (); LRbound () that in one dimensional integrable system one can understand the evolution of entanglement entropy of a subsystem based on quasiparticle picture CCevolution (). The argument is as follows: after the quench, there is an extensive excess in energy which appears as quasiparticles that propagate in time. The quasiparticles emitted from nearby points are entangled and they are responsible for the linear increase of the entanglement entropy of a subsystem with respect to the rest. In this section, we first study the time evolution of formation probabilities and subsequently Shannon and mutual information after a quantum quench. One can consider this section as a complement to the other studies of information propagation after quantum quench. To keep the discussion as simple as possible, we will concentrate on the most simple case of XXchain or free fermions. Following PeschelEisler () consider the Hamiltonian
(32) 
The time evolution of the correlation functions in the half filling are given as
(33) 
where, is the Bessel function of the first kind. Here we consider the dimerized initial conditions with and . The dimerized nature of the initial state will help later to consider different possibilities for the initial Shannon mutual information. Then at time zero we change the Hamiltonian to and let it evolve. The time evolution of the correlation matrix is given by PeschelEisler ()
(34) 
To calculate the time evolution of the probability of different configurations one just needs to use the above formula in (25). The results for few configurations are shown in the Figure 10. Of course since the sum of all the probabilities should be equal to one some of the probabilities increase with time and some decrease. All the probabilities change rapidly up to time and after that saturate. One can also simply calculate the evolution of the Shannon information with the tools of previous sections. In Figure 11 we depicted the evolution of Shannon information of a subsystem with respect to the time . The numerical results show an increase in the Shannon information up to time and then saturation. This is similar to what we usually have in the study of the time evolution of von Neumann entanglement entropy after quantum quench CCevolution (). However, one should be careful that in contrast to the von Neumann entropy the Shannon information of the subsystem is not a measure of correlation between the two subsystems. In addition, the increase in the Shannon information of the subsystem is not linear as the evolution of the von Neumann entanglement entropy. Our numerical results indicate that apart from a small regime at the beginning the Shannon information increases as
(35) 
where and and are positive independent quantities.
To study the time evolution of correlations, it is much better to study another quantity, Shannon mutual information of two subsystems. To investigate this quantity we first studied the time evolution of Shannon mutual information of a couple of dimers located far from each other. The results depicted in the Figure 12 show that the Shannon mutual information of the dimers are zero up to time and after that increases rapidly and then again decays slowly. This picture is consistent with the quasiparticle picture. The two regions are not correlated up to time that the quasiparticles emitted from the middle point reach each dimerfootnote3 (). However, the similarity between the evolution of the von Neumann entropy and mutual Shannon information ends here. To elaborate on that we consider the mutual information of two adjacent regions with sizes . Because of the dimerized nature of the initial state there are two possibilities for choosing the subsystems: at time zero at the boundary between the two subsystems there can be a dimer or not. In the first case at time zero the Shannon mutual information between the two subsystems is not zero but in the second case it is zero. In the second case naturally one expects an overall increase in the mutual information but in the first case a priory it is not clear that the mutual information should increase or decrease. In the Figure 13 we have depicted the results of the numerics for the two adjacent subsystems for different sizes. The numerical results show that for the uncorrelated initial conditions the Shannon mutual information first increases rapidly and then it decays and finally saturates at time . In the correlated case, we have overall decay in the mutual information and finally the saturation again at time . This behavior is very different from the quantum mutual information of the same regions which for the considered initial states first increases linearly and then saturates at time . The interesting phenomena is that after a short initial regime that the Shannon mutual information is initial state dependent the system enters to a regime that this quantity is completely independent of initial state and it decreases ”almost” linearly and then saturates. This can be also easily seen from the equation (35), where we can simply drive
(36) 
The saturation regime is independent of the size of the subsystem, this is simply because in the equlibrium regime the Shannon mutual information follows the arealaw Wolf2008 () and so it is independent of the volume of the subsystems.
Vii Conclusions
In this paper, we employed Grassmann numbers to write the probability of occurrence of different configurations in free fermion systems with respect to the minors of a particular matrix. The formula gives a very efficient method to study the scaling properties of logarithmic formation probabilities in the critical XYchain. In particular, we showed that the logarithmic formation probabilities of crystal configurations are given by the CFT formulas for the critical transverse field Ising model. This is checked by studying the probabilities in the infinite and finite (periodic and open boundary conditions) chain. In the case of critical XXchain which has a symmetry just the configurations with follow the CFT formulas. The rest of the configurations decay like a Gaussian and do not show much universal behavior. We also studied the Shannon information of a subsystem in the transverse field Ising model and XXchain. In particular, for the Ising model, we showed that in the scaling limit just the configurations with a high number of up spins contribute to the scaling of the Shannon information. In principle, if one considers all the configurations, with our method one can not calculate the Shannon information with classical computers for sizes bigger than in a reasonable time. However, if one admits a controllable error in the calculation of Shannon information it is possible to hire the results of section V to go to higher sizes. It would be very nice to extend this aspect of our calculations further to calculate the universal quantities in the Shannon information with higher accuracy. For example, one interesting direction can be finding an explicit formula for the sum of different powers of principal minors of a matrix. This kind of formulas can be very useful to calculate analytically or numerically the Rényi entropy of the subsystem.
Finally, we also studied the evolution of formation probabilities after quantum quench in free fermion system. In this case, we prepared the system in the dimer configuration and then we let it evolve with homogeneous Hamiltonian. The evolution of Shannon information of a susbsytem shows a very similar behavior as the evolution of entanglement entropy after a quantum quench. Especially our calculations show that the saturation of the Shannon information of the subsystem occurs at the same time as the entanglement entropy. This is probably not surprising because the is also the time that the reduced density operator saturates.
It will be very nice to extend our calculations in few other directions. One direction can be investigating the evolution of mutual information after local quantum quenches as it is done extensively in the studies of the entanglement entropy CClocal (); MP (); Karevski (). The other interesting direction can be calculating the same quantities in other bases, especially those bases that do not have any direct connection to CFT.
Acknowledgments.
We thank M Fagotti and F Alcaraz for many useful discussions. The work of MAR was supported in part by CNPq. KN acknowledges support from the National Science Foundation under grant number PHY1314295. MAR also thanks ICTP for hospitality during a period which part of this work was completed.
Appendix A: Shannon information for critical transversefield Ising model
In this Appendix, we will provide more details regarding the Shannon information of transverse critical Ising chain and XX chain. The data regarding Shannon information for a subsystem with length up to is listed in the Table A1. Having the data, we checked many different functions with different parameters to study the coefficient of the logarithm. Needless to say increasing the possible parameters can make a difference in the final result. In Stephan2014 () the results fitted to
(37) 
show that the best value is . This can be also checked using the data provided in the Table A1. It is worth mentioning that one can also get reasonable results using the data up to for some formation probabilities (not all) if we consider extra terms in the fitting procedure. Although we found that the equation (37) is the most stable fit with the least standard deviation based on our results in the main text we found it is hard to exclude the term because it is present in all the configurations studied there. If one includes this term and does not add the terms the coefficient will be . If one keeps all the terms the result will be . The final conclusion is that as far as one justifies the presence of the terms in the Shannon information formula the best value for with the current available data is .
l  Shannon  l  Shannon 

1  0.473946633733778  21  9.094267377324401 
2  0.925441055292197  22  9.520258511384927 
3  1.367970612016317  23  9.946131954351737 
4  1.805854593071358  24  10.37189747498959 
5  2.240889870728481  25  10.79756367448224 
6  2.674003797245196  26  11.22313816533366 
7  3.105734740754158  27  11.64862771729621 
8  3.536422963908594  28  12.07403837729498 
9  3.966297046625437  29  12.49937556879184 
10  4.395517906953372  30  12.92464417475445 
11  4.824203084194648  31  13.34984860684562 
12  5.252441034545332  32  13.77499286454422 
13  5.473946633733777  33  14.21026702317442 
14  6.107833679024358  34  14.62511508430369 
15  6.535085171703405  35  15.05009939651494 
16  6.962089515106671  36  15.47503630258492 
17  7.388875612253789  37  15.89992835887542 
18  7.815467577831834  38  16.32477792018708 
19  8.241885740227190  39  16.74960160654153 
20  8.668147394540807 
References
 (1) P. Calabrese, J. Cardy, J. Phys. A 42:504005 (2009)
 (2) V. E. Korepin , A. G. Izergin , F. H. L. Essler and D. B. Uglov , Phys. Lett. A 190 182 (1994) and F. H. L. Essler, H. Frahm, A. G. Izergin and V. E. Korepin, Comm. Math. Phys. 174. 191 (1995)
 (3) N. Kitanine, JM. Maillet , N. Slavnov and V. Terras, J. Phys. A: Math. Gen 35, L385 (2002); N. Kitanine, JM. Maillet , N. Slavnov and V. Terras, J. Phys. A: Math. Gen. 35 L753 (2002); A. G. Abanov and V. E. Korepin, Nucl. Phys. B. 647 565 (2002); V. E. Korepin , S. Lukyanov, Y. Nishiyama and M. Shiroishi, Phys. Lett. A 312 21 (2003); K. K. Kozlowski, J. Stat. Mech. P02006 (2008); L. Cantini, J. Phys. A: Math. Theor. 45 135207 (2012)
 (4) M. Shiroishi, M. Takahashi and Y. Nishiyama, J. Phys. Soc. Jpn. 70 3535 (2001)
 (5) A. G. Abanov, F. Franchini, Phys. Lett. A, 316, 342 (2003) and J. Phys. A 38, 5069 (2005)
 (6) JM Stéphan, J. Stat. Mech. P05010 (2014)
 (7) JM Stéphan, S. Furukawa, G. Misguich, and V. Pasquier, Phys. Rev. B, 80, 184421 (2009).
 (8) JM Stéphan, G. Misguich, and V. Pasquier, Phys. Rev. B, 82, 125455 (2010);
 (9) M. Oshikawa [arXiv:1007.3739]
 (10) JM Stéphan, G. Misguich, and V. Pasquier, Phys. Rev. B 84, 195128 (2011).
 (11) D. J. Luitz, F. Alet, N. Laflorencie, Phys. Rev. Lett. 112, 057203 (2014)
 (12) D. J. Luitz, F. Alet, N. Laflorencie, Phys. Rev. B 89, 165106 (2014)
 (13) D. J. Luitz, N. Laflorencie, F. Alet, J. Stat. Mech. P08007 (2014)
 (14) D. J. Luitz, X. Plat, N. Laflorencie, F. Alet, Phys. Rev. B 90, 125105 (2014)
 (15) M. Hermanns, S. Trebst, Phys. Rev. B 89, 205107 (2014)
 (16) J. Helmes, JM. Stéphan, S. Trebst, Phys. Rev. B 92, 125144 (2015)
 (17) J. C. Getelina, F. C. Alcaraz, J. A. Hoyos, Phys. Rev. B 93, 045136 (2016)
 (18) M. M. Wolf, F. Verstraete, M. B. Hastings, J. I. Cirac, Phys. Rev. Lett. 100, 070502 (2008).
 (19) J. Wilms, M. Troyer, F. Verstraete, J. Stat. Mech. P10011 (2011)and J. Wilms, J. Vidal, F. Verstraete, S. Dusuel, J. Stat. Mech. P01023 (2012)
 (20) H. W. Lau, P. Grassberger, Phys. Rev. E 87, 022128 (2013)
 (21) O. Cohen, V. Rittenberg, T. Sadhu, J. Phys. A: Math. Theor. 48, 055002 (2015)
 (22) J. Um, H. Park, H. Hinrichsen, J. Stat. Mech. P10026 (2012)
 (23) F. C. Alcaraz, M. A. Rajabpour, Phys. Rev. Lett. 111, 017201 (2013)
 (24) JM Stéphan, Phys. Rev. B 90, 045424 (2014)
 (25) F. C. Alcaraz, M. A. Rajabpour, Phys. Rev. B, 90, 075132 (2014)
 (26) F. C. Alcaraz, M. A. Rajabpour, Phys. Rev. B, 91, 155122 (2015)
 (27) M. A. Rajabpour, EPL, 112, 66001 (2015)
 (28) M. A. Rajabpour, Phys. Rev. B 92, 075108 (2015)
 (29) M. A. Rajabpour, [arXiv:1512.03940]
 (30) E. Barouch and B. McCoy, Phys. Rev. A 3, 786 (1971)
 (31) M. C. Chung and I. Peschel, Phys. Rev. B, 62, 4191 (2000) and I. Peschel, J. Phys. A: Math. Gen. 36, L205 (2003)
 (32) T. Barthel, and U. Schollwock, Phys. Rev. Lett. 100, 100601 (2008)
 (33) E. Lieb, T. Schultz, and D. Mattis, Ann. Phys. (NewYork) 16, 407 (1961)
 (34) Similar oscillating behavior exist also in noncritical regime with high transverse magnetic field. In this case the ffect is already understood based on superconducting correlations of real fermions in Franchini ().
 (35) M. Fagotti, ”Entanglement and Correlations in exactly solvable models”, unpublished PhD thesis (2011), ETD UniversitÃ di Pisa.
 (36) Although not shown in the table it is worth mentioning that for those configurations that depending on the parity of the number of down spins we have oscillating behavior the extracted value of is not exactly the same for those cases that is odd or even. However we conjecture that this is just a finite size effect and the formula (21) is correct.
 (37) M. Greiner, O. Mandel, T. W. HÃ¤nsch, and I. Bloch, Nature 419, 51 (2002); T. Kinoshita, T. Wenger, D. S. Weiss, Nature 440, 900 (2006); S. Hofferberth, I. Lesanovsky, B. Fischer, T. Schumm, and J. Schmiedmayer, Nature 449, 324 (2007) and S. Trotzky Y.A. Chen, A. Flesch, I. P. McCulloch, U. SchollwÃ¶ck, J. Eisert, and I. Bloch, Nature Phys. 8, 325 (2012)
 (38) P. Calabrese, J. Cardy, J. Stat. Mech. P04010 (2005)
 (39) M. B. Plenio, J. Hartley and J. Eisert, New J. Phys. 6, 36 (2004); A. Perales and M. B. Plenio, J. Opt. B: Quantum Semi class. Opt. 7, S601S609 (2005); R. G. Unanyan and M. Fleischhauer, Phys. Rev. A 90, 062330 (2014); G. De Chiara, S. Montangero, P. Calabrese and R. Fazio, J. Stat. Mech. 0603, P001 (2006) and M. Fagotti and P. Calabrese 2008 Phys. Rev. A 78 010306 (2008)
 (40) S. Bravyi, M. B. Hastings, and F. Verstraete, Phys. Rev.Lett. 97, 050401 (2006) and J. Eisert, T. J. Osborne, Phys. Rev. Lett. 97, 150404 (2006)
 (41) J. Dubail, JM Stéphan, J. Stat. Mech. L03002 (2011)
 (42) H. Liu and S. J. Suh, Phys. Rev. Lett. 112, 011601 (2014)
 (43) M. Fagotti, M. Collura, [arXiv:1507.02678]
 (44) V. Eisler and I. Peschel, J. Stat. Mech. P06005 (2007) and J. Phys. A: Math. Theor. 42, 504003 (2009)
 (45) Recently in MM () it has been pointed out that because of the presence of integrals of motions instead of two quasiparticles one needs to consider four of them. It is not clear for us how this can change the results of our investigation.
 (46) P. Calabrese, J. Cardy, J. Stat. Mech. P10004 (2007)
 (47) M. Collura, P. Calabrese, J. Phys. A: Math. Theor. 46, 175001 (2013)
 (48) V. Eisler, D. Karevski, T. Platini, I. Peschel, J. Stat. Mech. P01023 (2008)