modificatons in Sec III

This commit is contained in:
Pierre-Francois Loos 2023-01-31 16:11:36 +01:00
parent df448725d3
commit 7d09735a52

View File

@ -310,11 +310,10 @@ The similarity renormalization group method aims at continuously transforming a
Therefore, the transformed Hamiltonian Therefore, the transformed Hamiltonian
\begin{equation} \begin{equation}
\label{eq:SRG_Ham} \label{eq:SRG_Ham}
\bH(s) = \bU(s) \, \bH \, \bU^\dag(s), \bH(s) = \bU(s) \, \bH \, \bU^\dag(s)
\end{equation} \end{equation}
depends on a flow parameter $s$, such that $\bH(s=0)$ is the initial untransformed Hamiltonian and $\bH(s=\infty)$ is the (block)-diagonal Hamiltonian. depends on a flow parameter $s$, such that $\bH(s=0)$ is the initial untransformed Hamiltonian and $\bH(s=\infty)$ is the (block-)diagonal Hamiltonian.
An evolution equation for $\bH(s)$ can be easily obtained by deriving Eq.~\eqref{eq:SRG_Ham} with respect to $s$. An evolution equation for $\bH(s)$ can be easily obtained by differentiating Eq.~\eqref{eq:SRG_Ham} with respect to $s$, yielding the flow equation
This gives the flow equation
\begin{equation} \begin{equation}
\label{eq:flowEquation} \label{eq:flowEquation}
\dv{\bH(s)}{s} = \comm{\boldsymbol{\eta}(s)}{\bH(s)}, \dv{\bH(s)}{s} = \comm{\boldsymbol{\eta}(s)}{\bH(s)},
@ -323,8 +322,8 @@ where $\boldsymbol{\eta}(s)$, the flow generator, is defined as
\begin{equation} \begin{equation}
\boldsymbol{\eta}(s) = \dv{\bU(s)}{s} \bU^\dag(s) = - \boldsymbol{\eta}^\dag(s). \boldsymbol{\eta}(s) = \dv{\bU(s)}{s} \bU^\dag(s) = - \boldsymbol{\eta}^\dag(s).
\end{equation} \end{equation}
To solve this equation at a cost inferior to the one of diagonalizing the initial Hamiltonian, one needs to introduce an approximation for $\boldsymbol{\eta}(s)$. To solve this equation at a lower cost than the one of diagonalizing the initial Hamiltonian, one must introduce an approximate form for $\boldsymbol{\eta}(s)$.
Before defining such an approximation, we need to define what are the blocks to suppress to obtain a block-diagonal Hamiltonian. \titou{Before defining such an approximation, we need to define what are the blocks to suppress to obtain a block-diagonal Hamiltonian.
The Hamiltonian is separated into two parts as The Hamiltonian is separated into two parts as
\begin{equation} \begin{equation}
\bH(s) = \underbrace{\bH^\text{d}(s)}_{\text{diagonal}} + \underbrace{\bH^\text{od}(s)}_{\text{off-diagonal}}, \bH(s) = \underbrace{\bH^\text{d}(s)}_{\text{diagonal}} + \underbrace{\bH^\text{od}(s)}_{\text{off-diagonal}},
@ -332,9 +331,10 @@ The Hamiltonian is separated into two parts as
and, by definition, we have the following condition on $\bH^\text{od}$ and, by definition, we have the following condition on $\bH^\text{od}$
\begin{equation} \begin{equation}
\bH^\text{od}(s=\infty) = \boldsymbol{0}. \bH^\text{od}(s=\infty) = \boldsymbol{0}.
\end{equation} \end{equation}}
\PFL{Move this part at the start of the section.}
In this work, we will use Wegner's canonical generator which is defined as \cite{Wegner_1994} In this work, we consider Wegner's canonical generator which is defined as \cite{Wegner_1994}
\begin{equation} \begin{equation}
\boldsymbol{\eta}^\text{W}(s) = \comm{\bH^\text{d}(s)}{\bH(s)} = \comm{\bH^\text{d}(s)}{\bH^\text{od}(s)}. \boldsymbol{\eta}^\text{W}(s) = \comm{\bH^\text{d}(s)}{\bH(s)} = \comm{\bH^\text{d}(s)}{\bH^\text{od}(s)}.
\end{equation} \end{equation}
@ -343,21 +343,21 @@ If this generator is used, the following condition is verified \cite{Kehrein_20
\label{eq:derivative_trace} \label{eq:derivative_trace}
\dv{s}\text{Tr}\left[ \bH^\text{od}(s)^2 \right] \leq 0 \dv{s}\text{Tr}\left[ \bH^\text{od}(s)^2 \right] \leq 0
\end{equation} \end{equation}
which implies that the matrix elements of the off-diagonal part will decrease in a monotonic way. which implies that the matrix elements of the off-diagonal part decrease in a monotonic way throughout the transformation.
Even more, the coupling coefficients associated with the highest energy determinants are removed first as will be evidenced by the perturbative analysis after. Even more, the coupling coefficients associated with the highest-energy determinants are removed first as we shall evidence in the perturbative analysis below.
The main flaw of this generator is that it generates a stiff set of ODE which is therefore difficult to solve numerically. \ant{ref} The main drawback of this generator is that it generates a stiff set of ODE which is therefore difficult to solve numerically. \ant{ref}
However, here we will not tackle the full SRG problem but only consider analytical low-order perturbative expressions so we will not be affected by this problem. \cite{Evangelista_2014,Hergert_2016} However, here we will not tackle the full SRG problem but only consider analytical low-order perturbative expressions so we will not be affected by this problem. \cite{Evangelista_2014,Hergert_2016}
Let us now perform the perturbative analysis of the SRG equations. Let us now perform the perturbative analysis of the SRG equations.
For $s=0$, the initial problem is For $s=0$, the initial problem is
\begin{equation} \begin{equation}
\bH(0) = \bH^\text{d}(0) + \lambda ~ \bH^\text{od}(0) \bH(0) = \bH^\text{d}(0) + \lambda \bH^\text{od}(0)
\end{equation} \end{equation}
where $\lambda$ is the usual perturbation parameter and the off-diagonal part of the Hamiltonian has been defined as the perturber. where $\lambda$ is the usual perturbation parameter and the off-diagonal part of the Hamiltonian has been defined as the perturber.
For finite values of $s$, we have the following perturbation expansion of the Hamiltonian For finite values of $s$, we have the following perturbation expansion of the Hamiltonian
\begin{equation} \begin{equation}
\label{eq:perturbation_expansionH} \label{eq:perturbation_expansionH}
\bH(s) = \bH^{(0)}(s) + \lambda ~ \bH^{(1)}(s) + \lambda^2 \bH^{(2)}(s) + \dots \bH(s) = \bH^{(0)}(s) + \lambda ~ \bH^{(1)}(s) + \lambda^2 \bH^{(2)}(s) + \cdots
\end{equation} \end{equation}
Hence, the generator $\boldsymbol{\eta}(s)$ admits a perturbation expansion as well. Hence, the generator $\boldsymbol{\eta}(s)$ admits a perturbation expansion as well.
Then, one can collect order by order the terms in Eq.~\eqref{eq:flowEquation} and solve analytically the low-order differential equations. Then, one can collect order by order the terms in Eq.~\eqref{eq:flowEquation} and solve analytically the low-order differential equations.
@ -369,7 +369,7 @@ Then, one can collect order by order the terms in Eq.~\eqref{eq:flowEquation} an
Finally, the SRG formalism exposed above will be applied to $GW$. Finally, the SRG formalism exposed above will be applied to $GW$.
The first step is to define the diagonal and off-diagonal parts of the $GW$ effective Hamiltonian. The first step is to define the diagonal and off-diagonal parts of the $GW$ effective Hamiltonian.
As hinted at the end of section~\ref{sec:gw}, the diagonal and off-diagonal parts will be defined as As hinted at the end of Sec.~\ref{sec:gw}, the diagonal and off-diagonal parts are defined as
\begin{align} \begin{align}
\label{eq:diag_and_offdiag} \label{eq:diag_and_offdiag}
\bH^\text{d}(s) &= \bH^\text{d}(s) &=