truncation

  • full name: tenpy.linalg.truncation

  • parent module: tenpy.linalg

  • type: module

Classes

Inheritance diagram of tenpy.linalg.truncation

TruncationError([eps, ov])

Class representing a truncation error.

Functions

decompose_theta_qr_based(old_qtotal_L, ...)

Performs a QR based decomposition of a matrix theta (= the wavefunction) and truncates it.

svd_theta(theta, trunc_par[, qtotal_LR, ...])

Performs SVD of a matrix theta (= the wavefunction) and truncates it.

truncate(S, options)

Given a Schmidt spectrum S, determine which values to keep.

Module description

Truncation of Schmidt values.

Often, it is necessary to truncate the number of states on a virtual bond of an MPS, keeping only the state with the largest Schmidt values. The function truncate() picks exactly those from a given Schmidt spectrum \(\lambda_a\), depending on some parameters explained in the doc-string of the function.

Further, we provide TruncationError for a simple way to keep track of the total truncation error.

The SVD on a virtual bond of an MPS actually gives a Schmidt decomposition \(|\psi\rangle = \sum_{a} \lambda_a |L_a\rangle |R_a\rangle\) where \(|L_a\rangle\) and \(|R_a\rangle\) form orthonormal bases of the parts left and right of the virtual bond. Let us assume that the state is properly normalized, \(\langle\psi | \psi\rangle = \sum_{a} \lambda^2_a = 1\). Assume that the singular values are ordered descending, and that we keep the first \(\chi_c\) of the initially \(\chi\) Schmidt values.

Then we decompose the untruncated state as \(|\psi\rangle = \sqrt{1-\epsilon}|\psi_{tr}\rangle + \sqrt{\epsilon}|\psi_{tr}^\perp\rangle\) where \(|\psi_{tr}\rangle = \frac{1}{\sqrt{1-\epsilon}} \sum_{a < \chi_c} \lambda_a|L_a\rangle|R_a\rangle\) is the truncated state kept (normalized to 1), \(|\psi_{tr}^\perp\rangle = \frac{1}{\sqrt{\epsilon}} \sum_{a >= \chi_c} \lambda_a |L_a\rangle|R_a\rangle\) is the discarded part (orthogonal to the kept part) and the truncation error of a single truncation is defined as \(\epsilon = 1 - |\langle \psi | \psi_{tr}\rangle |^2 = \sum_{a >= \chi_c} \lambda_a^2\).

Warning

For imaginary time evolution (e.g. with TEBD), you try to project out the ground state. Then, looking at the truncation error defined in this module does not give you any information how good the found state coincides with the actual ground state! (Instead, the returned truncation error depends on the overlap with the initial state, which is arbitrary > 0)

Warning

This module takes only track of the errors coming from the truncation of Schmidt values. There might be other sources of error as well, for example TEBD has also an discretization error depending on the chosen time step.