SciPost Submission Page
Counterdiabatic Hamiltonians for multistate LandauZener problem
by Kohji Nishimura, Kazutaka Takahashi
This Submission thread is now published as
Submission summary
Authors (as registered SciPost users):  Kohji Nishimura 
Submission information  

Preprint Link:  https://arxiv.org/abs/1805.06662v2 (pdf) 
Date accepted:  20180919 
Date submitted:  20180912 02:00 
Submitted by:  Nishimura, Kohji 
Submitted to:  SciPost Physics 
Ontological classification  

Academic field:  Physics 
Specialties: 

Approach:  Theoretical 
Abstract
We study the LandauZener transitions generalized to multistate systems. Based on the work by Sinitsyn et al. [Phys. Rev. Lett. 120, 190402 (2018)], we introduce the auxiliary Hamiltonians that are interpreted as the counterdiabatic terms. We find that the counterdiabatic Hamiltonians satisfy the zero curvature condition. The general structures of the auxiliary Hamiltonians are studied in detail and the timeevolution operator is evaluated by using a deformation of the integration contour and asymptotic forms of the auxiliary Hamiltonians. For several spin models with transverse field, we calculate the transition probability between the initial and final ground states and find that the method is useful to study nonadiabatic regime.
Published as SciPost Phys. 5, 029 (2018)
Author comments upon resubmission
Reply to the referees:
We first would like to thank the referees for reviewing our manuscript. We appreciate their useful comments and recommendations. We have modified the manuscript following the comments. We do hope this revised version deserves publication in SciPost. Below is our response for the comments.
Report 1:
Strengths: 1. Technical correctness
We have found a slightly incorrect description in Section 2.2. $\xi$ introduced in this section is not necessarily equal to the counterdiabatic term. It possibly includes additional contributions. However, the equation to be solved in the following analysis, the zero curvature condition, is unchanged. Therefore, this correction does not affect our conclusion. We also mention that the present analysis is not necessarily restricted to the real symmetric Hamiltonian.
Since the novelty here is not very clear, I'd suggest the authors to denote it explicitly and possibly think of other means to enhance the main message in comparison to previous works in this area.
The main claim of this paper is that the auxiliary Hamiltonians introduced by Sinitsyn et al. (reference [22] in the new manuscript) can be interpreted as the counterdiabatic terms of Shortcuts to Adiabaticity. While the introduction of the auxiliary Hamiltonians is artificial and the physical meaning of these Hamiltonians is unclear in the previous study, we have succeeded in giving a natural interpretation of these Hamiltonians by identifying them as the counterdiabatic terms. Our method can be applied to general LandauZener Hamiltonians ($\hat{H} = J\hat{Z}+\gamma t \hat{X}$, where $Z$ and $X$ are arbitrary Hermitian operators), hence it enables us to calculate transition probabilities approximately even if the Hamiltonian is not solvable. Moreover, we may possibly generalize our discussion to Hamiltonians which are not limited to the present LandauZener type.
Following the referee's comment, we have modified the introduction to stress the novelty of the present study.
In the introduction the authors write: "The same type of the Hamiltonian is used in the method of quantum annealing [13–17]". This sentence is very general and not entirely precise. I'd encourage the authors to elaborate more on this topic.
With respect to the latter point, but also in general, I'd like to recommend the authors to use "Quantum Spin Glasses, Annealing and Computation" by Tanaka, Tamura and Chakrabarti (e.g. Chs. 57). Additional references to the literature would also be appreciated.
Thanks for the comment. When performing quantum annealing, one uses the Hamiltonian $H(t) = J\hat{Z}+\Gamma(t) \hat{X}$ where the operator $Z$ represents the problem of the combinatorial optimization and the operator $X$ is a transverse magnetic field. The parameter $\Gamma(t)$ is decreased slowly in the course of the time evolution to obtain the optimal solution. If we consider the schedule which is linear in $t$, $\Gamma(t) = \gamma t$, and consider the time evolution from $t=\infty$ to $t=\infty$, the protocol is the same as LandauZener protocol. Due to the similarity of these two protocols, we expect that LandauZener picture is applicable to analysing the performance of quantum annealing, which was actually discussed in literatures such as Refs.[20, 21] in the new manuscript.
Following the referee's suggestion, we have added a brief introduction to quantum annealing and several references ([17], [19] and [20] in the new manuscript).
Report 2:
Weaknesses: The content lacks novelty. There is no result that adds to the general scenario known already.
Thanks for the recommendation and the comment. We think our work is of importance due to the following three reasons.
First, no previous studies exist that consider the relation between the auxiliary Hamiltonians introduced by Sinitsyn et al. (reference [22] in the new manuscript) and the counterdiabatic terms of Shortcuts to Adiabaticity. We have suggested that the auxiliary Hamiltonians can be interpreted as the counterdiabatic terms. Accordingly, we have succeeded in giving a natural interpretation to the auxiliary Hamiltonians. In the previous study, the auxiliary Hamiltonians were introduced artificially, whose physical meaning was unclear.
Second, our theory can be applied to the general LandauZener Hamiltonians which are not necessarily solvable. The previous studies ([2224] in the new manuscript) treated solvable Hamiltonians. We can develop a systematic approximation method to treat the transition probabilities of the nonsolvable Hamiltonian, as discussed in Section 4.
Finally, our formulation may possibly be applicable to the Hamiltonian which is not limited to LandauZener type (e.g. the Hamiltonian $\hat{H} = J\hat{Z}+\gamma f(t) \hat{X}$ with arbitrary time dependence $f(t)$) since we can discuss the correspondence of the auxiliary Hamiltonians to the counterdiabatic terms for arbitrary Hamiltonians, in principle.
We have added some explanations to the manuscript so that the novelty of the present study becomes clear.
Strengths: 2. The presentation is rigorous and clear.
As stated in the reply to the report 1, we have modified an incorrect description in our manuscript, which does not affect our conclusion.
Requested changes: 1. The language has to be revised thoroughly.
Thanks for the advice. We have checked the manuscript thoroughly to improve the presentation.
Report 3:
Strengths: 1. Analyses are clear.
As stated in the reply to the report 1 and report 2, we have modified an incorrect description in our manuscript, which does not affect our conclusion.
Comment. The claim of this manuscript sounds vague. To my understanding, the main point of this manuscript should be the invention of an approximate computation method for the transition probability between the initial and final ground states in a generic multilevel LandauZenertype model. However, the manuscript stresses auxiliary Hamiltonians rather than transition probabilities. I suggest that Introduction should be revised so that the letter is stressed more.
Thanks for the useful comment. The main claim we would like to stress in this paper is the correspondence between the auxiliary Hamiltonians and the counterdiabatic terms appeared in Shortcuts to Adiabaticity. However, as the referee pointed out, we consider that a systematic calculation of the transition probabilities is as important as the above claim. Our theory is not only for a formal discussion but also for a practical use.
According to the request, we have modified the introduction to stress the importance of the second point.
Question. It is stated that Eq. (45) is exactly the same as the BrundoblerElser (BE) formula in Sec. 4. This statement is surprising and mysterious to me. This is because, while the BE formula is exact, Eq. (45) is obtained after the zeroth order approximation in g of $Z_1(\tau)$ and the secondorder cumulant expansion. The statement implies that, although the approximation used here should be valid for small $g$ and $\delta$, the resultant formula is exact for any $g$ and $\delta$. Why does this happen? I would like to ask the authors to solve this mystery.
This is a highly nontrivial result and we have not been able to explain why this happens. We suppose that this behavior comes from the fact that the exponent of the BE formula depends only on the factor $\delta^2 (=J^2/\gamma)$ and does not depend on other degrees of $\delta$. The exact mechanism cannot be explained in our series expansion method and a nonperturbative approach would be needed to answer this question completely.
We have added a comment below Eq. (47) in the new manuscript.
List of changes
1. Modified fifth (paragraph of Eq.(4)) and seventh paragraphs in the introduction to stress the novelty of the present study.
2. Added a brief introduction to quantum annealing in the third paragraph of the introduction.
3. Added references [17], [19], [20] on quantum annealing.
4. Added a sentence below Eq. (45) on the nontriviality of the equation.
5. Modified the second and third paragraphs in section 2.2.
6. Presentation improved, and grammatical errors and typos corrected.