![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00126v1-figure-040843c0-6ac8-4d9e-8b7a-310b1b39a485.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00126v1-figure-d011d5c0-7230-4c11-a13e-f4b3eb71374d.jpg)

## Tensor network approach to real-time path integral

## Shinji Takeda ∗

Institute of Physics, Kanazawa University, Kanazawa 920-1192, Japan E-mail: takeda@hep.s.kanazawa-u.ac.jp

We present a tensor network representation of the path integral for the one-component real scalar field theory in 1+1 dimensional Minkowski space-time. It is numerically verified by comparing with the exact result in the non-interacting case.

37th International Symposium on Lattice Field Theory - Lattice2019 16-22 June 2019

Wuhan, China

∗ Speaker.

https://pos.sissa.it/

## 1. Introduction

The real-time dynamics is one of the most important topics in modern physics, however, the field theoretical approaches based on the path integral representation suffer from the severe sign problem. So far several approaches have been proposed to overcome the sign problem. For example, the complex Langevin method is used to study φ 4 scalar field theory [1] and the SU(2) gauge theory [2] in 3+1 dimensional Minkowski space-time. In the course of studies, however, it turns out that the convergence issue gets harder when the time extent is taken relatively larger. Recently, the generalized Lefschetz thimble algorithm is used in the study of the 1+1 dimensional φ 4 theory [5] with the Schwinger-Keldysh setup [3, 4]. This approach, however, also has a limitation on the lattice size and it seems hard to perform larger volume simulations.

Tensor network method can also avoid the sign problem and is another promising approach. The idea of tensor networks was invented in the condensate matter physics and developed with a help of quantum information [6, 7, 8, 9]. Recently, the idea and technique have been used in elementary particle physics. It is known that within a framework of the tensor networks there are two types of approaches: the Hamiltonian approach and the Lagrangian one. The former one uses a tensor network in an ansatz form of a trial wave function for the variational method, and the real-time dynamics is studied in various two dimensional models [10, 11, 12]. The letter approach is nothing but the path integral and we will use it in the following. The aim of the study here is to demonstrate how to use the tensor network technique when evaluating the real-time path integral. As a first trial, we work on 1+1 dimensional real scalar field theory with the Minkowskian metric. In the following, we assume the lattice units a = 1.

## 2. Tensor network representation

Weconsider 1+1 dimensional square lattice system whose coordinate is denoted by x =( x 0 , x 1 ) with x 0 = 0 1 2 , , , ..., T -1 and x 1 = 0 1 2 , , , ..., L -1 for time and space respectively. In the following T = L is assumed. The boundary condition for both directions is taken to be periodic just for a simplicity, though this is not a physically meaningful setting. The one-component real scalar fields on the lattice are denoted by φ x ∈ R .

The lattice action for the 1+1 dimensional real scalar φ 4 theory with the Minkowskian metric is given by

$$S [ \phi ] = \sum _ { x } \left [ \frac { 1 } { 2 } ( \phi _ { x + \hat { 0 } } - \phi _ { x } ) ^ { 2 } - \frac { 1 } { 2 } ( \phi _ { x + \hat { 1 } } - \phi _ { x } ) ^ { 2 } - V ( \phi _ { x } ) \right ]$$

with the potential term

$$V ( \phi ) = \frac { 1 } { 2 } ( m _ { 0 } ^ { 2 } - i \varepsilon ) \phi _ { x } ^ { 2 } + \frac { \lambda } { 4! } \phi _ { x } ^ { 4 },$$

where m 0 , λ ∈ R are the bare parameters and the Feynman prescription parameter ε ( ε ∈ R , ε &gt; 0) is introduced and this plays an important role when forming the tensor network representation as we will see later. In an actual computation, note that we will have to take ε → 0 limit for a bare physical quantity.

In order to use the tensor network method, first of all, one has to rewrite the path integral in terms of the tensor network representation. For that purpose, let us rewrite it as a product of local factors,

$$\mathcal { Z } \equiv \int [ d \phi ] e ^ { i S [ \phi ] } = \int [ d \phi ] \prod _ { x } H ^ { ( 0 ) } ( \phi _ { x }, \phi _ { x + \dagger } ) H ^ { ( 1 ) } ( \phi _ { x }, \phi _ { x + \dagger } ),$$

where the local factors for the time-direction H ( 0 ) and for the space H ( 1 ) are explicitly given by 1

$$H ^ { ( 0 ) } ( \phi, \phi ^ { \prime } ) & = \exp \left [ - i \phi \phi ^ { \prime } - \frac { i } { 4 } V ( \phi ) - \frac { i } { 4 } V ( \phi ^ { \prime } ) \right ], \\ _ { \nu \nu }. \. & \quad. \coprod \. & i \quad. \. \.$$

$$H ^ { ( 1 ) } ( \phi, \phi ^ { \prime } ) = \exp \overset { \lambda } { \left [ + i \phi \phi ^ { \prime } - \frac { i } { 4 } V ( \phi ) - \frac { i } { 4 } V ( \phi ^ { \prime } ) \overset { \lambda } { \right ] } } \,.$$

They are considered as a two-variable function. An important property of the local factors is that thanks to the presence of the Feynman prescription parameter they are Hilbert-Schmidt operator

$$\int _ { - \infty } ^ { \infty } d \phi d \phi ^ { \prime } | H ^ { ( \mu ) } ( \phi, \phi ^ { \prime } ) | ^ { 2 } < \infty \quad \text{ for } \mu = 0, 1.$$

This fact immediately means that they are compact operator and can be expanded as 2

$$H ^ { ( \mu ) } ( \phi, \phi ^ { \prime } ) = \sum _ { k = 0 } ^ { \infty } \Psi _ { k } ^ { ( \mu ) } ( \phi ) \, \sigma _ { k } \Phi _ { k } ^ { ( \mu ) * } ( \phi ^ { \prime } ),$$

where Ψ and Φ are orthonormal basis function satisfying

$$\int _ { - \infty } ^ { \infty } d \phi \Psi _ { m } ^ { * } ( \phi ) \Psi _ { n } ( \phi ) & = \int _ { - \infty } ^ { \infty } d \phi \, \Phi _ { m } ^ { * } ( \phi ) \Phi _ { n } ( \phi ) = \delta _ { m n }, & ( 2. 8 ) \\ \sum _ { k = 0 } ^ { \infty } \Psi _ { k } ( \phi ) \Psi _ { k } ^ { * } ( \phi ^ { \prime } ) & = \sum _ { k = 0 } ^ { \infty } \Phi _ { k } ( \phi ) \Phi _ { k } ^ { * } ( \phi ^ { \prime } ) = \delta ( \phi - \phi ^ { \prime } ),$$

and σ k represents singular value and is non-negative [13, 14].

As a next step let us see how to obtain the basis functions and the singular values. First, we expand H ( µ )

$$H ^ { ( \mu ) } ( \phi, \phi ^ { \prime } ) = \sqrt { 2 \pi } \sum _ { m, n = 0 } ^ { \infty } \psi _ { m } ( \phi ) X _ { m n } ^ { ( \mu ) } \psi _ { n } ( \phi ^ { \prime } ) & & ( 2. 1 0 )$$

in terms of the Hermite function ψ m which is the eigen-function of the harmonic oscillator

$$\psi _ { n } ( x ) = \frac { 1 } { \sqrt { \pi ^ { 1 / 2 } n! 2 ^ { n } } } H _ { n } ( x ) e ^ { - x ^ { 2 } / 2 }, \quad \ \ H _ { n } ( x ) \, \colon \text{Hermite polynomials}, \quad \ \ ( 2. 1 1 )$$

and satisfies the orthonormal and the completeness relations as in eq.(2.8) and (2.9). Basically H ( µ ) are an oscillating function thus it is reasonable to use ψ m as a basis since the latter is also such a

1 Note that the phase factors exp ( i φ 2 x / 2 ) associated with the kinetic terms has canceled in eq.(2.4) and (2.5).

2 One can show that the singular value is independent of µ = 0 1. Moreover, it can be shown that it is independent , of the bare parameters m 0 and λ since they are included in the phase factor in eq.(2.4) and (2.5). Therefore the singular value is a function of only ε , while the bare parameter dependence are captured in the basis functions.

function. By using the orthonormal property of ψ m , the coefficient matrix Xmn ( µ ) , which is a complex valued symmetric matrix, is obtained as

$$X _ { m n } ^ { ( \mu ) } = \frac { 1 } { \sqrt { 2 \pi } } \int _ { - \infty } ^ { \infty } d \phi d \phi ^ { \prime } \psi _ { m } ( \phi ) H ^ { ( \mu ) } ( \phi, \phi ^ { \prime } ) \psi _ { n } ( \phi ^ { \prime } ). & & ( 2. 1 2 ) \\. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &.. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &.$$

An evaluation of the double integral for such an oscillating function is rather technically demanding, thus here we use a trick in the following. In order to separate the integration into two parts, we use a formula 3

$$e ^ { \mp i \phi \phi ^ { \prime } } = \sqrt { 2 \pi } \sum _ { n = 0 } ^ { \infty } ( \mp i ) ^ { n } \psi _ { n } ( \phi ) \psi _ { n } ( \phi ^ { \prime } ).$$

By using the formula together with eq.(2.4) and (2.5), Xmn ( µ ) is given by

$$X _ { m n } ^ { ( \mu ) } = \sum _ { k = 0 } ^ { \infty } ( - 1 ) ^ { k \delta _ { \mu 0 } k } G _ { m k } G _ { n k } \quad \text{with} \quad G _ { m n } = \int _ { - \infty } ^ { \infty } d \phi \, \psi _ { m } ( \phi ) \, \psi _ { n } ( \phi ) \exp \left [ - \frac { i } { 4 } V ( \phi ) \right ]. \quad ( 2. 1 4 )$$

The resulting single integration for Gmn , which is again a complex valued symmetric matrix, can be exactly carried out for the free case and then Gmn is given by the hypergeometric function, while for the interacting case it may be estimated by some numerical integration scheme with a suitable deformation of integration path. After obtaining Xmn ( µ ) numerically, we apply SVD to it,

$$X _ { m n } ^ { ( \mu ) } = \sum _ { k = 0 } ^ { \infty } U _ { m k } ^ { ( \mu ) } \sigma _ { k } ( V ^ { ( \mu ) \dagger } ) _ { k n }, \quad \text{with} \quad U ^ { ( \mu ) }, V ^ { ( \mu ) } \, \colon \, \text{unitary matrix.}$$

By inserting the above equation into eq.(2.10), H ( µ ) are given as

$$\text{By inserting the above equation into eq.} ( 2. 1 0 ), H ^ { \nu ^ { \prime } } \text{ are given as} \\ H ^ { ( \mu ) } ( \phi, \phi ^ { \prime } ) & = \sqrt { 2 \pi } \sum _ { m, n = 0 } ^ { \infty } \psi _ { m } ( \phi ) \left ( \sum _ { k = 0 } ^ { \infty } U _ { m k } ^ { ( \mu ) } \sigma _ { k } ( V ^ { ( \mu ) \dagger } ) _ { k n } \right ) \psi _ { n } ( \phi ^ { \prime } ) \\ & = \sqrt { 2 \pi } \sum _ { k = 0 } ^ { \infty } \left ( \sum _ { m = 0 } ^ { \infty } \psi _ { m } ( \phi ) U _ { m k } ^ { ( \mu ) } \right ) \sigma _ { k } \left ( \sum _ { n = 0 } ^ { \infty } ( V ^ { ( \mu ) \dagger } ) _ { k n } \psi _ { n } ( \phi ^ { \prime } ) \right ) \\ & = \sqrt { 2 \pi } \sum _ { k = 0 } ^ { \infty } \Psi _ { k } ^ { ( \mu ) } ( \phi ) \sigma _ { k } \Phi _ { k } ^ { ( \mu ) * } ( \phi ^ { \prime } ), & & ( 2. 1 6 ) \\ \text{where in the last step, we have defined new basis functions}$$

where in the last step, we have defined new basis functions

$$\Psi _ { k } ^ { ( \mu ) } ( \phi ) = \sum _ { m = 0 } ^ { \infty } \psi _ { m } ( \phi ) U _ { m k } ^ { ( \mu ) }, \quad \Phi _ { k } ^ { ( \mu ) * } ( \phi ) = \sum _ { n = 0 } ^ { \infty } ( V ^ { ( \mu ) \dagger } ) _ { k n } \psi _ { n } ( \phi ). \quad \ ( 2. 1 7 )$$

In this way, the basis functions and the singular values are obtained. In an actual calculation, one cannot deal with the infinite range of index for the matrix, thus one has to truncate it and should monitor the truncation errors of some physical quantities. Note that the singular values of H ( µ ) are the same as those of Xmn ( µ ) .

Final step is to obtain the tensor. By collecting all building blocks, the tensor is formed as

$$T _ { i j k l } = 2 \pi \sqrt { \sigma _ { i } \sigma _ { j } \sigma _ { k } \sigma _ { l } } \int _ { - \infty } ^ { \infty } d \phi \Psi _ { i } ^ { ( 0 ) } ( \phi ) \Psi _ { j } ^ { ( 1 ) } ( \phi ) \Phi _ { k } ^ { ( 0 ) * } ( \phi ) \Phi _ { l } ^ { ( 1 ) * } ( \phi ). \quad ( 2. 1 8 ) \\ \frac { 3 \tau \tau \colon \max \max \max \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \max \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \max \max \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \max \min \max \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \max \max \max \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \max \max \min \max \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \max \max \max \max \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \min \sum \text{$\ }$}$$

3 This may be derived from e ixp = √ 2 π 〈 x p | 〉 = √ 2 π ∑ ∞ n = 0 〈 x n | 〉〈 n p | 〉 and the Fourier transformation of the Hermite function.

Figure 1: Left panel: Singular values σ k of H ( µ ) ( φ φ , ′ ) (or equivalently Xmn ( µ ) ) with ε = 0 1, 0 5 and 1. Right . . panel: The absolute value of Re [ Tiiii ] and the imaginary part at λ = 0 and ε = 0 5 with various . m 0 values.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00126v1-figure-fe01260e-f3b5-43e2-8433-72779c0696bc.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00126v1-figure-c66432f5-10c8-4f06-be29-e3fb8d62acb1.jpg)

Note that the numbering of the basis function turns out to be an index of tensor i , j , k and l . Basically, the range of the indices is infinity but to put it on a computer one has to truncate it. The truncation order is denoted as N and the index range is restricted as 0 ≤ i , j , k l , ≤ N in the following. The integration in eq.(2.18) can be carried out exactly and it is given by the hypergeometric function again.

## 3. Numerical results

The resulting singular values σ k are shown in Fig. 1 (left). For large ε where a damping in the integrand gets stronger, the hierarchy is sharper as expected, while for small ε , the hierarchy is not clear and the information compression is not well realized.

To see the hierarchy structure of the tensor elements, Tiiii are shown in Fig. 1 (right). Basically the size of elements tends to be exponentially small for larger i . This behavior allows us to truncate the tensor index.

In order to check whether the initial tensor in eq.(2.18) is correctly made, we compute the path integral on 2 × 2 lattice with the periodic boundary condition. We show the real part of ln Z / 4 in Fig. 2 (left). Here note that we did not use a coarse-graining, thus this is purely a check of the initial tensor itself. The exact result and the numerical results with some truncation orders are consistent with each other at this scale. Right panel of Fig. 2 shows the relative deviation from the exact results

∣ ∣ By increasing N , the deviation δ tends to be smaller. Note that for smaller ε , the deviation is larger since the hierarchy of the singular value is not so good as seen before.

$$\delta = \left | \frac { \text{Re} [ \ln \mathcal { Z } _ { \text{exact} } ] - \text{Re} [ \ln \mathcal { Z } _ { N } ] } { \text{Re} [ \ln \mathcal { Z } _ { \text{exact} } ] } \right |. \quad \quad \quad ( 3. 1 ) \\ \text{deviation} \, \delta \, \text{ tends to be smaller. Note that for smaller $\varepsilon$, the deviation is larger} \\ \text{$a$ the median is not $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ and $a$ respectively}$$

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00126v1-figure-8c9b0822-d3da-4aa5-b38d-86d9e33df457.jpg)

ε

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00126v1-figure-92ca97e8-1c5a-44e1-93e7-2a0ad3d142a3.jpg)

Figure 2: Left panel: Re ln [ Z ] / V on 2 × 2 lattice as a function of ε with various truncation orders N = 16 -40 at m 0 = 0 1 and . λ = 0. The black line is the exact result. Right panel: The relative deviation for Re ln [ Z ] with the same parameter set as the left one.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00126v1-figure-41acd23f-097e-4ba4-a920-b0af62774f5e.jpg)

ε

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00126v1-figure-26925286-9753-4f5f-8102-bc9322154091.jpg)

ε

Figure 3: Left panel: Re ln [ Z ] / V on ( 1024 ) 2 lattice as a function of ε . Right panel: The associated relative deviation. The parameter set is the same as that of Fig. 2.

The path integral at larger volumes can be obtained by using the tensor renormalization group algorithm [15]. The results in Fig. 3 show that the coarse-graining step introduces visible errors. This part should be improved by using some sophisticated algorithms in future .

## 4. Summary and outlook

Wehave derived the tensor network representation for the scalar field theory with the Minkowskian metric and numerically checked the validity of the formulation for the free case. One of the important point in the procedure is to expand the local factor H ( µ ) in terms of the orthonormal basis functions. The another key point is to introduce the Feynman prescription ε that not only provides a damping factor in the the path integral but also dictates the hierarchy of the singular values of H ( µ ) . Note that although there is no sign problem in the formulation here, another problem emerges as an "information incompressibility problem" which means that the hierarchy of the singular value tends to be worse for smaller ε . Furthermore we have to eliminate the regulator ( ε → 0) in the end of the calculation. This is what we have to pay a price in our formulation.

As future outlooks, there are many things to do. For example, when the ε is small, the initial tensor has large truncation errors thus one may have to improve the initial tensor itself. Secondly, when making the initial tensor for the interacting case, one needs an efficient numerical method to evaluate the oscillating integral, say the steepest descent method and so on. Thirdly, instead of the Feynman prescription, one may use the tilted time axis as a regulator. Finally, one may straightforwardly extend to the Euclidean space as well as the Schwinger-Keldysh formulation in a similar way. Once such a setup is formulated, the transportation coefficients will be accessible.

This work is supported in part by JSPS KAKENHI Grant Numbers JP17K05411 and MEXT as 'Exploratory Challenge on Post-K computer' (Frontiers of Basic Science: Challenging the Limits).

## References

- [1] J. Berges and I.-O. Stamatescu, Phys. Rev. Lett. 95 , 202003 (2005)
- [2] J. Berges, S. Borsanyi, D. Sexty and I.-O. Stamatescu, Phys. Rev. D 75 , 045007 (2007)
- [3] J. S. Schwinger, J. Math. Phys. 2 , 407 (1961).
- [4] L. V. Keldysh, Zh. Eksp. Teor. Fiz. 47 , 1515 (1964) [Sov. Phys. JETP 20 , 1018 (1965)].
- [5] A. Alexandru, G. Basar, P. F. Bedaque and G. W. Ridgway, Phys. Rev. D 95 , no. 11, 114501 (2017)
- [6] S. R. White, Phys. Rev. Lett. 69 , 2863 (1992).
- [7] S. Ostlund and S. Rommer, Phys. Rev. Lett. 75 , 3537 (1995)
- [8] F. Verstraete and J. I. Cirac, cond-mat/0407066.
- [9] U. Schollwock, Rev. Mod. Phys. 77 , 259 (2005)
- [10] T. Pichler, M. Dalmonte, E. Rico, P. Zoller and S. Montangero, Phys. Rev. X 6 , no. 1, 011023 (2016)
- [11] S. Kuehn, J. I. Cirac and M. C. Banuls, JHEP 1507 , 130 (2015)
- [12] D. Lin et al, PoS(LATTICE2019)022.
- [13] W. Lay and J. Rundnick, Phys. Rev. Lett. 88 , 057203 (2002)
- [14] Y. Shimizu, Mod. Phys. Lett. A 27 , 1250035 (2012).
- [15] M. Levin and C. P. Nave, Phys. Rev. Lett. 99 , no. 12, 120601 (2007)
Shinji Takeda
Tensor network approach to real-time path integral
We present a tensor network representation of the path integral for the
one-component real scalar field theory in 1+1 dimensional Minkowski space-time.
It is numerically verified by comparing with the exact result in the
non-interacting case.