Differences

This shows you the differences between two versions of the page.

--- 数値計算ライブラリ:itensor:1サイト [2017/06/13 18:13] – [プログラムの解説] koudai
+++ 数値計算ライブラリ:itensor:1サイト [2021/06/27 22:04] (current) – external edit 127.0.0.1
@@ Line 2: / Line 2: @@
 もっとも簡単な例として、シングルサイトにおけるスピン演算子と波動関数を定義する。
 ====== 理論 ======
@@ Line 54: / Line 52: @@
 となる。
+状態$| \psi \rangle$におけるスピン演算子の期待値は
+\begin{equation}
+ \langle \psi | S^x | \psi \rangle = \sum_{s, s'} \psi^*_s S^x_{s, s'} \psi_{s'}
+\end{equation}
+となる。
 ====== ソースコード ======
@@ Line 138: / Line 141: @@
     }
 </sxh>
-Makefileの APP=one を APP=01_tutorial_code に書き換える。
-<sxh; first-line: 1; highlight: [6]; title: Makefile>
-LIBRARY_DIR=../..
-ifdef app
-APP=$(app)
-else
-APP=01_tutorial_code
-endif
-CCFILES=$(APP).cc
-#################################################################
-#################################################################
-#################################################################
-#################################################################
-include $(LIBRARY_DIR)/this_dir.mk
-include $(LIBRARY_DIR)/options.mk
-TENSOR_HEADERS=$(LIBRARY_DIR)/itensor/all.h
-#Mappings --------------
-OBJECTS=$(patsubst %.cc,%.o, $(CCFILES))
-GOBJECTS=$(patsubst %,.debug_objs/%, $(OBJECTS))
-#Rules ------------------
-%.o: %.cc $(HEADERS) $(TENSOR_HEADERS)
-	$(CCCOM) -c $(CCFLAGS) -o $@ $<
-.debug_objs/%.o: %.cc $(HEADERS) $(TENSOR_HEADERS)
-	$(CCCOM) -c $(CCGFLAGS) -o $@ $<
-#Targets -----------------
-build: $(APP)
-debug: $(APP)-g
-$(APP): $(OBJECTS) $(ITENSOR_LIBS)
-	$(CCCOM) $(CCFLAGS) $(OBJECTS) -o $(APP) $(LIBFLAGS)
-$(APP)-g: mkdebugdir $(GOBJECTS) $(ITENSOR_GLIBS)
-	$(CCCOM) $(CCGFLAGS) $(GOBJECTS) -o $(APP)-g $(LIBGFLAGS)
-clean:
-	rm -fr .debug_objs *.o $(APP) $(APP)-g
-mkdebugdir:
-	mkdir -p .debug_objs
-</sxh>
-プログラムのコンパイルと実行は、プログラムのあるディレクトリに入って次のようにする。
-> \$ make
-> \$ ./01_tutorial_code
 ====== プログラムの解説 ======
@@ Line 204: / Line 147: @@
 求めたい演算子や状態に対応するテンソルを定義していくことで計算を進める。
 テンソルの定義やテンソル同士の演算に便利なクラスが、ITensorではライブラリとして提供されている。
+===== ライブラリのインクルード =====
 <sxh C++; first-line: 1>
@@ Line 212: / Line 157: @@
 all.hはITensorのすべてのクラスがインクルードされるが、プログラムが重くなったりコンパイルに時間がかかったりするので、慣れてくれば必要なヘッダファイルだけインクルードすればよい。
-テンソルの添え字$s$を定義する。
+===== テンソルの添字の定義 =====
+テンソルの添字$s$を定義する。
 <sxh C++; first-line: 12>
     auto s = Index("s",2);
 </sxh>
-クラスautoにより、添字オブジェクトsを定義する。
+クラスIndexにより、添字オブジェクトsを定義する。
 書き方は Index("[添字に使う記号]", [添字の次元]) である。
 今はアップとダウンからなるので、添え字の次元は2である。
-さて、添え字の定義ができたので、状態と演算子を構成する。
+===== テンソルの定義 ======
+添え字の定義ができたので、状態と演算子を構成する。
 状態$\psi_s$を次のように作成する。
 <sxh C++; first-line: 15>
@@ Line 238: / Line 187: @@
 </sxh>
-次にスピン演算子$S^x_{s, s'}$, $S^z_{s,s'}$を作成する。
+次にスピン演算子$S^z_{s,s'}$および$S^x_{s, s'}$を作成する。
 階のテンソルを作りたい場合は次のようにする。
 <sxh C++; first-line: 30>
@@ Line 244: / Line 193: @@
     auto Sx = ITensor(s,prime(s));
 </sxh>
-メソッド .prime(s) は、次元が添字$s$と同じ添字$s'$を生成する。
+prime(s) によって次元が添字$s$と同じである添字$s'$を生成する。
 このままだと演算子のすべての成分がゼロのままなので、値を設定する。
 <sxh C++; first-line: 33>
@@ Line 253: / Line 202: @@
     Sx.set(s(2),prime(s)(1),+0.5);
 </sxh>
+===== テンソル同士の積 =====
+テンソル同士の積は * によって計算できる。
+<sxh C++; first-line: 46>
+    ITensor phi = Sx * psi;
+</sxh>
+ここで注意したいのが、ここで計算するのは
+\begin{equation}
+ \phi_{s'} = \sum_{s} S^x_{s,s'} \psi_s
+\end{equation}
+だということである。
+積を定義したとき、同じ添字に関して和をとることになる。$\phi_{s'}$を$\phi_s$とするためには
+<sxh C++; first-line: 48>
+    phi.noprime();
+</sxh>
+とすればよい。
+===== 期待値 =====
+最後に、演算子の期待値の求め方を紹介する。
+状態$\psi_s$とエルミート共役$\psi^*_{s'}$は次のように作る。
+<sxh C++; first-line: 68>
+    auto cpsi = dag(prime(psi));
+</sxh>
+スピン演算子の期待値は積 cpsi * Sz * psi あるいは cpsi * Sx * psi により得られるが、このままだとオブジェクトのままなので、実数として扱うためにメソッド.realを使う。
+<sxh C++; first-line: 68>
+    Real zz = (cpsi * Sz * psi).real();
+    Real xx = (cpsi * Sx * psi).real();
+</sxh>
+ここで計算したのは
+\begin{equation}
+ \sum_{s,s'} \psi^*_{s'} S^x_{s, s'} \psi_{s}
+\end{equation}
+である（添字の順番に注意）。