diff options
40 files changed, 2555 insertions, 0 deletions
diff --git a/Jauslin_LMU_2021.tex b/Jauslin_LMU_2021.tex new file mode 100644 index 0000000..ab433b8 --- /dev/null +++ b/Jauslin_LMU_2021.tex @@ -0,0 +1,398 @@ +\documentclass{ian-presentation} + +\usepackage[hidelinks]{hyperref} +\usepackage{graphicx} +\usepackage{array} +\usepackage{xcolor} + + +\definecolor{ipurple}{HTML}{4B0082} +\definecolor{iyellow}{HTML}{DAA520} +\definecolor{igreen}{HTML}{32CD32} +\definecolor{iblue}{HTML}{4169E1} +\definecolor{ired}{HTML}{DC143C} + +\definecolor{highlight}{HTML}{328932} +\definecolor{highlight}{HTML}{981414} + +\begin{document} +\pagestyle{empty} +\hbox{}\vfil +\bf\Large +\hfil An effective equation to study Bose gases\par +\smallskip +\hfil at all densities\par +\vfil +\large +\hfil Ian Jauslin +\normalsize +\vfil +\hfil\rm joint with {\bf Eric A. Carlen}, {\bf Markus Holzmann}, {\bf Elliott H. Lieb}\par +\vfil +arXiv:{\tt\ \parbox[b]{3cm}{ + \href{https://arxiv.org/abs/1912.04987}{1912.04987}\par + \href{https://arxiv.org/abs/2010.13882}{2010.13882}\par + \href{https://arxiv.org/abs/2011.10869}{2011.10869} +}} +\hfill{\tt \href{http://ian.jauslin.org}{http://ian.jauslin.org}} +\eject + +\setcounter{page}1 +\pagestyle{plain} + +\title{Bose-Einstein condensation} +\begin{itemize} +\item System of many Bosons, e.g. {\color{highlight}Helium} atoms, {\color{highlight}Rubidium} atoms, etc... +\item {\color{highlight}Bose-Einstein condensate}: most particles are in the same quantum state. +\item Related to the phenomena of {\color{highlight}superfluidity} (flow with zero viscocity) and {\color{highlight}superconductivity} (currents with zero resistance). +\item Predicted theoretically in {\color{highlight}1924-1925}, experimentally observed in {\color{highlight}1995}. +\item Mathematical understanding: still {\color{highlight}no proof} of the existence of a condensate (at finite density, in the presence of interactions and in the continuum). +\end{itemize} +\vfill +\eject + +\title{Repulsive Bose gas} +\begin{itemize} + \item Potential: {\color{highlight}$v(r)\geqslant 0$} and {\color{highlight}$v\in L_1(\mathbb R^3)$}, Hamiltonian: + $$ + H_N:= + -\frac12\sum_{i=1}^N\Delta_i + +\sum_{1\leqslant i<j\leqslant N}v(|x_i-x_j|) + $$ + \vskip-5pt + \item Ground state: $\psi_0$, energy $E_0$. + \item Observables in the {\color{highlight}thermodynamic limit}: ground state energy per particle and condensate fraction: $P_i$: projection onto condensate state + $$ + e_0:=\lim_{\displaystyle\mathop{\scriptstyle V,N\to\infty}_{{\color{highlight}\frac NV=\rho}}}\frac{E_0}N + ,\quad + \eta_0:=\lim_{\displaystyle\mathop{\scriptstyle V,N\to\infty}_{{\color{highlight}\frac NV=\rho}}}\frac1N\sum_{i=1}^N\left<\psi_0\right|P_i\left|\psi_0\right> + . + $$ +\end{itemize} +\vfill +\eject + +\title{Low density} +\begin{itemize} + \item Bogolyubov theory: {\color{highlight}approximation scheme} that reduces the problem to an effective {\color{highlight}1-particle problem}. + \item Predictions \href{https://doi.org/10.1103/PhysRev.106.1135}{[Lee, Huang, Yang, 1957]}: + \begin{itemize} + \item Energy: + $$ + {\color{highlight}e_0=2\pi\rho a\left(1+\frac{128}{15\sqrt\pi}\sqrt{\rho a^3}+o(\sqrt\rho)\right)} + $$ + \vskip-10pt + \item Condensate fraction: + $$ + {\color{highlight}1-\eta_0\sim\frac{8\sqrt{\rho a^3}}{3\sqrt\pi}} + $$ + \end{itemize} +\end{itemize} +\vfill +\eject + +\title{Low density} +\begin{itemize} + \item Energy asymptotics: {\color{highlight} proved}: + \href{https://doi.org/10.1103/PhysRevLett.80.2504}{[Lieb, Yngvason, 1998]}, + \href{https://doi.org/10.1007/s10955-009-9792-3}{[Yau, Yin, 2009]}, + \href{https://doi.org/10.4007/annals.2020.192.3.5}{[Fournais, Solovej, 2020]}. + + \item Condensate fraction: {\color{highlight}still open} in the theormodynamic limit, but there are proofs of condensation in the Gross-Pitaevskii regime (ultra-dilute): + \href{https://doi.org/10.1103/PhysRevLett.88.170409}{[Lieb, Seiringer, 2002]}, + \href{https://doi.org/10.1007/s00220-017-3016-5}{[Boccato, Brennecke, Cenatiempo, Schlein, 2018]}. + +\end{itemize} +\vfill +\eject + +\title{High density} +\begin{itemize} + \item [Bogolyubov, 1947]: if $\hat v\geqslant 0$. + $$ + {\color{highlight}e_0\sim\frac\rho2\int v} + $$ + {\color{highlight}Hartree} (mean field) energy. + \item {\color{highlight}Proved} in \href{https://doi.org/10.1103/PhysRev.130.2518}{[Lieb, 1963]}. + \item Condensate fraction + $$ + \eta\to1 + $$ + {\color{highlight}open}. +\end{itemize} +\vfill +\eject + +\title{Energy as a function of density} +For $v(x)=e^{-|x|}$: + +\hfil\includegraphics[height=5.5cm]{erho_lowhigh.pdf} +\vfill +\eject + +\addtocounter{page}{-1} +\title{Energy as a function of density} +For $v(x)=e^{-|x|}$: + +\hfil\includegraphics[height=5.5cm]{erho_effective.pdf} +\vfill +\eject + +\title{Derivation of the equation} +\begin{itemize} + \item \href{https://doi.org/10.1103/PhysRev.130.2518}{[Lieb, 1963]}. + \item Integrate $H_N\psi_0=E_0\psi_0$: + $$ + \int dx_1\cdots dx_N\ + \left( + -\frac12\sum_{i=1}^N\Delta_i\psi_0 + +\sum_{1\leqslant i<j\leqslant N} v(x_i-x_j)\psi_0 + \right) + =E_0\int dx_1\cdots dx_N\ \psi_0 + $$ + \item Therefore, + $$ + \frac{N(N-1)}2\int dx_1dx_2\ v(x_1-x_2)\frac{\int dx_3\cdots dx_N\ \psi_0}{\int dx_1\cdots dx_N\ \psi_0} + =E_0 + $$ +\end{itemize} +\vfill +\eject + +\title{Derivation of the equation} +\begin{itemize} + \item Thus, + $$ + \frac{E_0}N=\frac{N-1}{2V}\int dx\ v(x)g_2(0,x) + $$ + \item {\color{highlight}$\psi_0\geqslant 0$}, so it can be thought of as a probability distribution. + \item $g_n$: {\color{highlight}correlation functions} of $V^{-N}\psi_0$ + $$ + g_n(x_1,\cdots,x_n):=\frac{V^n\int dx_{n+1}\cdots dx_N\ \psi_0(x_1,\cdots,x_N)}{\int dx_1\cdots dx_N\ \psi_0(x_1,\cdots,x_N)} + $$ +\end{itemize} +\vfill +\eject + +\title{Hierarchy} +\vskip-10pt +\begin{itemize} + \item Equation for $g_2$: integrate $H_N\psi_0=E_0\psi_0$ with respect to $x_3,\cdots,x_N$: + $$ + \begin{array}{>\displaystyle l} + -\frac12(\Delta_x+\Delta_y) g_2(x,y) + +\frac{N-2}V\int dz\ (v(x-z)+v(y-z))g_3(x,y,z) + \\[0.5cm]\hfill + +v(x-y)g_2(x,y) + +\frac{(N-2)(N-3)}{2V^2}\int dzdt\ v(z-t)g_4(x,y,z,t) + =E_0g_2(x,y) + \end{array} + $$ + \item Factorization {\color{highlight}assumption}: + $$ + g_3(x_1,x_2,x_3)=g_2(x_1,x_2)g_2(x_1,x_3)g_2(x_2,x_3) + $$ + $$ + g_4(x_1,x_2,x_3,x_4)=\prod_{i<j}(g_2(x_i,x_j)+O(V^{-1})) + $$ +\end{itemize} +\vfill +\eject + +\title{{\color{ipurple}Big equation}} +\begin{itemize} +\item In the thermodynamic limit, if ${\color{highlight}u(x):=1-g_2(0,x)}$, + $$ + -\Delta u(x) + = + (1-u(x))\left(v(x)-2\rho K(x)+\rho^2 L(x)\right) + $$ + $$ + K:= + u\ast S + ,\quad + S(y):=(1-u(y))v(y) + $$ + $$ + L:= + u\ast u\ast S + -2u\ast(u(u\ast S)) + +\frac12 + \int dydz\ u(y)u(z-x)u(z)u(y-x)S(z-y) + . + $$ + + \item {\color{ipurple}``Big'' equation}: + $$ + L\approx + u\ast u\ast S + -2u\ast (u(u\ast S)) + $$ +\end{itemize} +\vfill +\eject + +\title{{\color{iblue}Simple equation}} +\vskip-10pt +\begin{itemize} + \item Further approximate $S(x)\approx\frac{2e}\rho\delta(x)$ and $u\ll 1$. + \item Simple equation + $$ + {\color{iblue}-\Delta u(x)=(1-u(x))v(x)- 4eu(x)+2e\rho\ u\ast u(x)} + $$ + $$ + {\color{iblue}e=\frac\rho2\int dx\ (1-u(x))v(x)} + $$ + \item {\bf Theorem 1}: + If $v(x)\geqslant 0$ and $v\in L_1\cap L_2(\mathbb R^3)$, then the {\color{iblue}simple equation} has an {\color{highlight}integrable solution} (proved constructively), with $0\leqslant u\leqslant 1$. +\end{itemize} +\vfill +\eject + +\title{Energy for the {\color{iblue}simple equation}} +\vskip-10pt +\begin{itemize} + \item {\bf Theorem 2}: + $$ + \frac{e}{\rho}\mathop{\longrightarrow}_{\rho\to\infty}\frac12\int dx\ v(x) + $$ + (note that there is no condition that $\hat v\geqslant 0$). + This coincides with the {\color{highlight}Hartree energy}. + \item {\bf Theorem 3}: + $$ + e=2\pi\rho a\left(1+\frac{128}{15\sqrt\pi}\sqrt{\rho a^3}+o(\sqrt\rho)\right) + $$ + This coincides with the {\color{highlight}Lee-Huang-Yang prediction}. +\end{itemize} +\vfill +\eject + +\title{Energy} +$v(x)=e^{-|x|}$, Blue: {\color{iblue}simple equation}; purple: {\color{ipurple}big equation}; red: {\color{ired}Monte Carlo} + +\hfil\includegraphics[height=5.5cm]{erho_fulleq.pdf} +\vfill +\eject + +\title{Energy} +$v(x)=e^{-|x|}$, Blue: {\color{iblue}simple equation}; red: {\color{ired}Jastrow}; purple: {\color{ipurple}big equation} + +\hfil\includegraphics[height=5.5cm]{cmp_energy1.pdf} +\vfill +\eject + +\title{Condensate fraction} +\begin{itemize} + \item Add a parameter $\mu$ to the Hamiltonian: + $$ + H_N(\mu):= + -\frac12\sum_{i=1}^N\Delta_i + +\sum_{1\leqslant i<j\leqslant N}v(x_i-x_j) + {\color{highlight}-\mu\sum_{i=1}^NP_i} + $$ + \item Projection onto condensate wavefunction: $P_i$. + \item Condensate fraction: + $$ + {\color{highlight}\eta_0}:=\frac1N\left<\psi_0\right|\sum_{i=1}^NP_i\left|\psi_0\right> + =-\frac1N\partial_\mu \left<\psi_0\right|H_N(\mu)\left|\psi_0\right>|_{\mu_0} + \equiv + {\color{highlight}-\partial_\mu e_0(\mu)|_{\mu=0}} + $$ +\end{itemize} +\vfill +\eject + +\title{Condensate fraction} +\begin{itemize} + \item {\bf Theorem 4}: + For the {\color{iblue}simple equation}, as $\rho\to0$ + $$ + 1-\eta\sim\frac{8\sqrt{\rho a^3}}{3\sqrt\pi} + $$ + which coincides with {\color{highlight}Bogolyubov's prediction}. +\end{itemize} +\vfill +\eject + +\title{Condensate fraction} +$v(x)=e^{-|x|}$, Blue: {\color{iblue}simple equation}; purple: {\color{ipurple}big equation}; red: {\color{ired}Monte Carlo} + +\hfil\includegraphics[height=5.5cm]{condensate.pdf} +\vfill +\eject + +\title{Two point correlation function} +$v(x)=16e^{-|x|}$, Blue: {\color{iblue}simple equation}; purple: {\color{ipurple}big equation}; red: {\color{ired}Monte Carlo} + +\hfil\includegraphics[height=5.5cm]{correlation.pdf} +\vfill +\eject + +\title{The uniqueness problem} + $$ + -\Delta u(x)=(1-u(x))v(x)- 4eu(x)+2e\rho\ u\ast u(x) + ,\quad + e=\frac\rho2\int dx\ (1-u(x))v(x) + $$ +\begin{itemize} + \item Change the point of view: {\color{highlight}fix $e>0$}, and compute $\rho$ and $u$. + + \item {\color{highlight}Iteration}: $u_0=0$, + $$ + (-\Delta+4e+v)u_n=v+2e\rho_{n-1}u_{n-1}\ast u_{n-1} + ,\quad + \rho_n:=\frac{2e}{\int dx\ (1-u_n(x))v(x)} + . + $$ +\end{itemize} +\vfill +\eject + +\title{The uniqueness problem} +\begin{itemize} + \item {\bf Lemma}: $u_n(x)$ is an {\color{highlight}increasing} sequence, and is {\color{highlight}bounded} $u_n(x)\leqslant 1$. + It converges to a function $u$, which is the {\color{highlight}unique} integrable solution of the equation {\color{highlight}with $e$ fixed}. + + \item {\bf Lemma}: $e\mapsto\rho(e)$ is {\color{highlight}continuous}, and $\rho(0)=0$ and $\rho(\infty)=\infty$, which allows us to compute solutions for the problem at fixed $\rho$. + + \item We thus have a {\color{highlight}restricted} notion of uniqueness. + The full uniqueness would follow from a proof that $e\mapsto\rho(r)$ is {\color{highlight}monotone increasing} (which must be true for the physics to make sense). +\end{itemize} + +\vfill +\eject + +\title{Limitations of the simple and big equations} +\begin{itemize} + \item Only works at high densities for {\color{highlight}$\hat v\geqslant 0$}. + \item Less accurate for {\color{highlight}large potentials}: for $v(x)=16e^{-|x|}$, + + \hfil\includegraphics[width=5.5cm]{energy16.pdf} + \hfil\includegraphics[width=5.5cm]{condensate16.pdf} +\end{itemize} +\vfill +\eject + +\title{Conclusions and outlook} +\begin{itemize} + \item Two {\color{highlight}effective equations}: the {\color{ipurple}big equation} and the {\color{iblue}simple equation}, which are {\color{highlight}non-linear 1-particle equations}. + \item Reproduce the known results for both {\color{highlight}small and large densities}. + \item Their derivation is {\color{highlight}different from Bogolyubov theory}, so they may give new insights onto studying the Bose gas in these asymptotic regimes. + \item The {\color{ipurple}big equation} is {\color{highlight}quantitatively accurate} at intermediate densities. + \item This opens up the possibility of studying the physics of the {\color{highlight}Bose gas at intermediate densities}. +\end{itemize} +\vfill +\eject + +\title{Open problems} +\begin{itemize} + \item Analysis of the {\color{iblue}simple equation}: {\color{highlight}Monotonicity} of $e(\rho)$, and {\color{highlight}convexity} of $\rho e(\rho)$. (So far, we have proofs for small and large $\rho$.). Similarly, prove that $0\leqslant\eta\leqslant 1$. (We have a proof for small $\rho$.) + + \item Analysis of the {\color{ipurple}big equation}: everything is still open. + + \item Relate these equations back to the {\color{highlight}many-body Bose gas}. + + \item Other setups: {\color{highlight}trapping potential}. +\end{itemize} + +\end{document} diff --git a/Makefile b/Makefile new file mode 100644 index 0000000..c885343 --- /dev/null +++ b/Makefile @@ -0,0 +1,45 @@ +PROJECTNAME=$(basename $(wildcard *.tex)) +LIBS=$(notdir $(wildcard libs/*)) +FIGS=$(notdir $(wildcard figs/*.fig)) + +PDFS=$(addsuffix .pdf, $(PROJECTNAME)) +SYNCTEXS=$(addsuffix .synctex.gz, $(PROJECTNAME)) + +all: $(PROJECTNAME) + +$(PROJECTNAME): $(LIBS) $(FIGS) + pdflatex -file-line-error $@.tex + pdflatex -synctex=1 $@.tex + +$(SYNCTEXS): $(LIBS) $(FIGS) + pdflatex -synctex=1 $(patsubst %.synctex.gz, %.tex, $@) + +$(LIBS): + ln -fs libs/$@ ./ + + +$(FIGS): + make -C figs/$@ + for pdf in $$(find figs/$@/ -name '*.pdf'); do ln -fs "$$pdf" ./ ; done + + +clean-aux: clean-figs-aux + rm -f $(addsuffix .aux, $(PROJECTNAME)) + rm -f $(addsuffix .log, $(PROJECTNAME)) + rm -f $(addsuffix .out, $(PROJECTNAME)) + +clean-libs: + rm -f $(LIBS) + +clean-figs: + $(foreach fig,$(addprefix figs/, $(FIGS)), make -C $(fig) clean; ) + rm -f $(notdir $(wildcard figs/*.fig/*.pdf)) + +clean-figs-aux: + $(foreach fig,$(addprefix figs/, $(FIGS)), make -C $(fig) clean-aux; ) + + +clean-tex: + rm -f $(PDFS) $(SYNCTEXS) + +clean: clean-aux clean-tex clean-libs clean-figs @@ -0,0 +1,33 @@ +This directory contains the source files to typeset the presentation, and +generate the figures. This can be accomplished by running + make + +This document uses a custom class file, located in the 'libs' directory, which +defines a number of commands. + + +* Dependencies: + + pdflatex + TeXlive packages: + amsfonts + array + graphics + hyperref + latex + pgf + standalone + GNU make + gnuplot + +* Files: + + Jauslin_LMU_2021.tex: + main LaTeX file + + libs: + custom LaTeX class file + + figs: + source code for the figures + diff --git a/figs/condensate.fig/Makefile b/figs/condensate.fig/Makefile new file mode 100644 index 0000000..c45de65 --- /dev/null +++ b/figs/condensate.fig/Makefile @@ -0,0 +1,34 @@ +PROJECTNAME=condensate + +PDFS=$(addsuffix .pdf, $(PROJECTNAME)) +DATS= simpleq.dat bigeq.dat + +all: $(PDFS) + +$(PDFS): $(DATS) + gnuplot $(patsubst %.pdf, %.gnuplot, $@) > $(patsubst %.pdf, %.tikz.tex, $@) + pdflatex -jobname $(basename $@) -file-line-error $(patsubst %.pdf, %.tikz.tex, $@) + +simpleq.dat: + julia $(SIMPLEQ)/main.jl -p "tolerance=1e-11;order=100;maxiter=100;v_a=1" -M simpleq condensate_fraction > $@ +bigeq.dat: + julia -p 8 $(SIMPLEQ)/main.jl -p "eq=bigeq;N=12;P=8;v_a=1" -M medeq condensate_fraction > $@ + +install: $(PDFS) + cp $^ $(INSTALLDIR)/ + +$(LIBS): + ln -fs libs/$@ ./ + +clean-libs: + rm -f $(LIBS) + +clean-aux: + rm -f $(addsuffix .aux, $(PROJECTNAME)) + rm -f $(addsuffix .log, $(PROJECTNAME)) + rm -f $(addsuffix .tex, $(PROJECTNAME)) + +clean-tex: + rm -f $(PDFS) + +clean: clean-libs clean-aux clean-tex diff --git a/figs/condensate.fig/bigeq.dat b/figs/condensate.fig/bigeq.dat new file mode 100644 index 0000000..e48cbc3 --- /dev/null +++ b/figs/condensate.fig/bigeq.dat @@ -0,0 +1,100 @@ + 1.000000000000000e-06 1.942109602248065e-03 1.702970665674168e-12 + 1.202264434617413e-06 2.144142093189251e-03 1.178808810520995e-13 + 1.445439770745928e-06 2.370859341242408e-03 1.192336354552407e-12 + 1.737800828749376e-06 2.623985546948715e-03 4.076932518598860e-12 + 2.089296130854041e-06 2.904135074720890e-03 3.065737407737518e-12 + 2.511886431509582e-06 3.210803051150611e-03 6.684762465649123e-13 + 3.019951720402019e-06 3.543752141002869e-03 4.093999456324028e-14 + 3.630780547701017e-06 3.904216206875963e-03 4.822691376657683e-14 + 4.365158322401657e-06 4.294899506208919e-03 3.259674352110494e-13 + 5.248074602497723e-06 4.719386978665083e-03 1.028819821235890e-12 + 6.309573444801930e-06 5.181659336263778e-03 1.861524874355775e-12 + 7.585775750291836e-06 5.685862176450733e-03 2.488338274080411e-12 + 9.120108393559096e-06 6.236245567146821e-03 2.709450451184446e-12 + 1.096478196143185e-05 6.837185816392036e-03 2.543679948619913e-12 + 1.318256738556407e-05 7.493237426386802e-03 2.102973889594259e-12 + 1.584893192461114e-05 8.209187621063666e-03 1.527249098895665e-12 + 1.905460717963248e-05 8.990099289731322e-03 9.540538061535407e-13 + 2.290867652767775e-05 9.841336884630313e-03 4.976651096293187e-13 + 2.754228703338163e-05 1.076857581585209e-02 2.551631314103556e-13 + 3.311311214825908e-05 1.177779891236849e-02 2.812676798242623e-13 + 3.981071705534969e-05 1.287528340380401e-02 6.307586000801811e-13 + 4.786300923226380e-05 1.406757951946043e-02 1.319160137450734e-12 + 5.754399373371566e-05 1.536147868007941e-02 2.314118180410572e-12 + 6.918309709189363e-05 1.676396654008137e-02 3.536707073357365e-12 + 8.317637711026709e-05 1.828215421336748e-02 4.866668455989243e-12 + 1.000000000000000e-04 1.992317980401387e-02 6.169192714323269e-12 + 1.202264434617413e-04 2.169407203311732e-02 7.386129325640002e-12 + 1.445439770745928e-04 2.360156902423473e-02 8.592071601482435e-12 + 1.737800828749376e-04 2.565188877539273e-02 4.878829234412997e-15 + 2.089296130854039e-04 2.785045276854119e-02 3.942521128486340e-15 + 2.511886431509580e-04 3.020156819707873e-02 6.234437426160288e-15 + 3.019951720402016e-04 3.270807501963133e-02 1.399974577150145e-14 + 3.630780547701014e-04 3.537096132498831e-02 5.005525484699284e-15 + 4.365158322401661e-04 3.818894722948108e-02 5.456358250864572e-15 + 5.248074602497728e-04 4.115803806089322e-02 4.917806683980876e-15 + 6.309573444801930e-04 4.427105486160373e-02 6.680870445775179e-15 + 7.585775750291835e-04 4.751716358764096e-02 1.779168320111825e-14 + 9.120108393559096e-04 5.088143908449703e-02 9.760013767224759e-15 + 1.096478196143185e-03 5.434450964117910e-02 9.521419345061477e-15 + 1.318256738556408e-03 5.788232915180864e-02 7.574188114522546e-15 + 1.584893192461114e-03 6.146611885354680e-02 1.039377505239065e-14 + 1.905460717963246e-03 6.506251478122342e-02 8.176786147785561e-15 + 2.290867652767772e-03 6.863395318622874e-02 8.338301801453393e-15 + 2.754228703338166e-03 7.213932061012299e-02 7.581270910814022e-15 + 3.311311214825911e-03 7.553488145864165e-02 7.156938338864503e-15 + 3.981071705534973e-03 7.877547091592903e-02 6.929415027760092e-15 + 4.786300923226385e-03 8.181590892502735e-02 8.674664895759950e-15 + 5.754399373371567e-03 8.461255959422054e-02 7.526793721152772e-15 + 6.918309709189363e-03 8.712493630911380e-02 4.144857506288649e-15 + 8.317637711026709e-03 8.931723938183408e-02 3.530020645632524e-15 + 1.000000000000000e-02 9.115971178417075e-02 6.924046664381005e-15 + 1.202264434617413e-02 9.262971027971965e-02 6.034088541831901e-15 + 1.445439770745928e-02 9.371241367329312e-02 4.885148310946663e-15 + 1.737800828749376e-02 9.440112413816591e-02 2.917320379507791e-15 + 2.089296130854041e-02 9.469715660989715e-02 2.610686638602541e-15 + 2.511886431509582e-02 9.460934902764143e-02 7.354666705237897e-15 + 3.019951720402019e-02 9.415325731511556e-02 5.796225767460134e-15 + 3.630780547701010e-02 9.335011969104240e-02 6.017084081620531e-15 + 4.365158322401656e-02 9.222568368303644e-02 5.302998003473651e-15 + 5.248074602497723e-02 9.080898675648411e-02 4.039838886119182e-15 + 6.309573444801930e-02 8.913117007641815e-02 4.125816885656187e-15 + 7.585775750291836e-02 8.722438779915097e-02 4.273284938741890e-15 + 9.120108393559097e-02 8.512085476424038e-02 2.623094380645173e-15 + 1.096478196143185e-01 8.285205637207869e-02 4.927649012520671e-15 + 1.318256738556407e-01 8.044812781958981e-02 2.420044464772237e-15 + 1.584893192461114e-01 7.793739685903543e-02 2.149025610829404e-15 + 1.905460717963248e-01 7.534607518083469e-02 3.388825035919124e-15 + 2.290867652767775e-01 7.269807815050612e-02 2.822970856541473e-15 + 2.754228703338169e-01 7.001495034676258e-02 4.182630891668936e-15 + 3.311311214825908e-01 6.731587439715549e-02 2.904212967331170e-15 + 3.981071705534969e-01 6.461774223707760e-02 5.477274398838064e-15 + 4.786300923226380e-01 6.193527047174405e-02 4.783552039455272e-15 + 5.754399373371566e-01 5.928114448311252e-02 4.439456789833779e-15 + 6.918309709189363e-01 5.666617892504009e-02 3.191356061395679e-15 + 8.317637711026709e-01 5.409948505148902e-02 4.658959396388796e-15 + 1.000000000000000e+00 5.158863779039475e-02 2.496311679602705e-15 + 1.202264434617413e+00 4.913983755302447e-02 2.503975647353570e-15 + 1.445439770745928e+00 4.675806345070745e-02 1.980624000103322e-15 + 1.737800828749376e+00 4.444721590565470e-02 3.513326611439523e-15 + 2.089296130854041e+00 4.221024763556332e-02 2.128775949903540e-15 + 2.511886431509582e+00 4.004928271457441e-02 2.084732824400213e-15 + 3.019951720402019e+00 3.796572391712868e-02 2.057863428784457e-15 + 3.630780547701010e+00 3.596034888344765e-02 2.149696722374579e-15 + 4.365158322401657e+00 3.403339584700401e-02 2.506590750983091e-15 + 5.248074602497723e+00 3.218463976982702e-02 2.666959749115830e-15 + 6.309573444801930e+00 3.041345976812730e-02 1.843390243429058e-15 + 7.585775750291836e+00 2.871889869982439e-02 2.075866723337706e-15 + 9.120108393559097e+00 2.709971574419922e-02 2.614004879642500e-15 + 1.096478196143185e+01 2.555443274573189e-02 2.784239762991527e-15 + 1.318256738556407e+01 2.408137502873176e-02 4.417788851818592e-15 + 1.584893192461114e+01 2.267870732135474e-02 2.346418536174242e-15 + 1.905460717963248e+01 2.134446535879641e-02 3.159630166116848e-15 + 2.290867652767775e+01 2.007658366786368e-02 3.067450163255898e-15 + 2.754228703338169e+01 1.887291997221849e-02 4.085090221270331e-15 + 3.311311214825908e+01 1.773127660229970e-02 4.662766446257259e-15 + 3.981071705534969e+01 1.664941924664940e-02 5.219648685301088e-15 + 4.786300923226381e+01 1.562509334029200e-02 5.188731853763820e-15 + 5.754399373371567e+01 1.465603834912039e-02 5.298056449394350e-15 + 6.918309709189363e+01 1.374000017615125e-02 6.088475789849585e-15 + 8.317637711026708e+01 1.287474188620420e-02 5.607884807044724e-15 diff --git a/figs/condensate.fig/condensate.gnuplot b/figs/condensate.fig/condensate.gnuplot new file mode 100644 index 0000000..f1c689d --- /dev/null +++ b/figs/condensate.fig/condensate.gnuplot @@ -0,0 +1,36 @@ +set ylabel "$\\eta$" norotate +set xlabel "$\\rho$" + +set yrange [0.9:1] +set ytics 0.91, 0.03 +set mytics 3 + +set xrange [1e-6:100] +set xtics 1e-6, 100, 100 +set xtics add ("$10^{-6}$" 0.000001, "$10^{-4}$" 0.0001, "$10^{-2}$" 0.01, "$1$" 1.0, "$10^2$" 100) +unset mxtics + +# default output canvas size: 12.5cm x 8.75cm +set term lua tikz size 8,6 standalone + +set key off + +# set linestyle +set style line 1 linetype rgbcolor "#4169E1" linewidth 2 +set style line 2 linetype rgbcolor "#DC143C" linewidth 2 +set style line 3 linetype rgbcolor "#32CD32" linewidth 2 +set style line 4 linetype rgbcolor "#4B0082" linewidth 2 +set style line 5 linetype rgbcolor "#DAA520" linewidth 2 + +set pointsize 0.6 + +set label at 0.000005,0.91 "Bogolyubov" textcolor "#32CD32" + +set logscale x + +a=1.25435641059 + +plot "simpleq.dat" u 1:(1-$2) w l ls 1, \ + "bigeq.dat" u 1:(1-$2) w l ls 4 , \ + 1-8/(3*sqrt(pi))*sqrt(x*a**3) linestyle 3 dashtype "." ,\ + "holzmann_2019-12-25.dat" u 1:(1-$3) ls 2 diff --git a/figs/condensate.fig/holzmann_2019-12-25.dat b/figs/condensate.fig/holzmann_2019-12-25.dat new file mode 100644 index 0000000..96689cb --- /dev/null +++ b/figs/condensate.fig/holzmann_2019-12-25.dat @@ -0,0 +1,11 @@ +## data from M. Holzmann, 2019-09-22 +# rho E0 n0 +1e-6 7.902e-6 +1e-4 8.3441e-4 0.0171 +1e-3 9.13384e-3 0.0481 +1e-2 1.061073e-1 0.0871 +1e-1 1.19182e+0 0.0851 +1e-0 1.24302e+1 0.0587 +1e+1 1.25442e+2 0.0319 +5e+1 6.28032e+2 0.021 + diff --git a/figs/condensate.fig/simpleq.dat b/figs/condensate.fig/simpleq.dat new file mode 100644 index 0000000..348121c --- /dev/null +++ b/figs/condensate.fig/simpleq.dat @@ -0,0 +1,100 @@ + 1.000000000000000e-06 2.095081711925015e-03 + 1.202264434617413e-06 2.295235280091131e-03 + 1.445439770745928e-06 2.514296186738034e-03 + 1.737800828749376e-06 2.754005992218143e-03 + 2.089296130854041e-06 3.016257261240048e-03 + 2.511886431509582e-06 3.303104592145314e-03 + 3.019951720402019e-06 3.616775975726918e-03 + 3.630780547701017e-06 3.959684344305409e-03 + 4.365158322401657e-06 4.334439155646415e-03 + 5.248074602497723e-06 4.743857767750661e-03 + 6.309573444801930e-06 5.190976310677877e-03 + 7.585775750291836e-06 5.679059686385030e-03 + 9.120108393559096e-06 6.211610171947045e-03 + 1.096478196143185e-05 6.792374027485015e-03 + 1.318256738556407e-05 7.425345312781278e-03 + 1.584893192461114e-05 8.114765932045185e-03 + 1.905460717963248e-05 8.865120699464721e-03 + 2.290867652767775e-05 9.681125931789915e-03 + 2.754228703338163e-05 1.056770975064344e-02 + 3.311311214825908e-05 1.152998189333092e-02 + 3.981071705534969e-05 1.257319039145271e-02 + 4.786300923226380e-05 1.370266197723624e-02 + 5.754399373371566e-05 1.492372253381352e-02 + 6.918309709189363e-05 1.624159332185176e-02 + 8.317637711026709e-05 1.766125812610249e-02 + 1.000000000000000e-04 1.918729591875359e-02 + 1.202264434617413e-04 2.082367319967555e-02 + 1.445439770745928e-04 2.257348995912804e-02 + 1.737800828749376e-04 2.443867336254203e-02 + 2.089296130854039e-04 2.641961398528058e-02 + 2.511886431509580e-04 2.851474099050601e-02 + 3.019951720402016e-04 3.072003537104471e-02 + 3.630780547701014e-04 3.302848464372067e-02 + 4.365158322401661e-04 3.542948859291789e-02 + 5.248074602497728e-04 3.790823417663501e-02 + 6.309573444801930e-04 4.044506877087432e-02 + 7.585775750291835e-04 4.301491450154334e-02 + 9.120108393559096e-04 4.558678199565728e-02 + 1.096478196143185e-03 4.812345827946150e-02 + 1.318256738556408e-03 5.058145864328128e-02 + 1.584893192461114e-03 5.291134289200051e-02 + 1.905460717963246e-03 5.505849827399160e-02 + 2.290867652767772e-03 5.696447960043553e-02 + 2.754228703338166e-03 5.856896678930269e-02 + 3.311311214825911e-03 5.981234782238041e-02 + 3.981071705534973e-03 6.063886048364867e-02 + 4.786300923226385e-03 6.100013371714988e-02 + 5.754399373371567e-03 6.085886970213172e-02 + 6.918309709189363e-03 6.019231791545348e-02 + 8.317637711026709e-03 5.899513433487237e-02 + 1.000000000000000e-02 5.728121486315967e-02 + 1.202264434617413e-02 5.508415858187780e-02 + 1.445439770745928e-02 5.245615721904299e-02 + 1.737800828749376e-02 4.946530746709376e-02 + 2.089296130854041e-02 4.619156855992348e-02 + 2.511886431509582e-02 4.272179185226004e-02 + 3.019951720402019e-02 3.914438480508303e-02 + 3.630780547701010e-02 3.554420649152192e-02 + 4.365158322401656e-02 3.199821953636445e-02 + 5.248074602497723e-02 2.857226663349834e-02 + 6.309573444801930e-02 2.531914059520501e-02 + 7.585775750291836e-02 2.227792229852985e-02 + 9.120108393559097e-02 1.947440831436047e-02 + 1.096478196143185e-01 1.692235880228233e-02 + 1.318256738556407e-01 1.462526717442258e-02 + 1.584893192461114e-01 1.257837345877809e-02 + 1.905460717963248e-01 1.077069465848733e-02 + 2.290867652767775e-01 9.186909564013418e-03 + 2.754228703338169e-01 7.808998508205012e-03 + 3.311311214825908e-01 6.617591964706389e-03 + 3.981071705534969e-01 5.593022042416338e-03 + 4.786300923226380e-01 4.716097668760250e-03 + 5.754399373371566e-01 3.968639441621027e-03 + 6.918309709189363e-01 3.333816454742892e-03 + 8.317637711026709e-01 2.796327604398189e-03 + 1.000000000000000e+00 2.342466331442259e-03 + 1.202264434617413e+00 1.960102269258443e-03 + 1.445439770745928e+00 1.638607138480099e-03 + 1.737800828749376e+00 1.368746272215081e-03 + 2.089296130854041e+00 1.142551829784601e-03 + 2.511886431509582e+00 9.531892670181336e-04 + 3.019951720402019e+00 7.948250154423670e-04 + 3.630780547701010e+00 6.625005224425910e-04 + 4.365158322401657e+00 5.520157121243856e-04 + 5.248074602497723e+00 4.598234179345307e-04 + 6.309573444801930e+00 3.829352915009887e-04 + 7.585775750291836e+00 3.188389988971121e-04 + 9.120108393559097e+00 2.654260840747906e-04 + 1.096478196143185e+01 2.209296357978901e-04 + 1.318256738556407e+01 1.838707817059403e-04 + 1.584893192461114e+01 1.530130079086820e-04 + 1.905460717963248e+01 1.273233330370514e-04 + 2.290867652767775e+01 1.059394293104648e-04 + 2.754228703338169e+01 8.814186420712528e-05 + 3.311311214825908e+01 7.333072418003886e-05 + 3.981071705534969e+01 6.100596972384632e-05 + 4.786300923226381e+01 5.075095479639026e-05 + 5.754399373371567e+01 4.221862079919067e-05 + 6.918309709189363e+01 3.511994492291246e-05 + 8.317637711026708e+01 2.921428437169262e-05 diff --git a/figs/condensate16.fig/.energy16.gnuplot.swp b/figs/condensate16.fig/.energy16.gnuplot.swp Binary files differnew file mode 100644 index 0000000..d23db08 --- /dev/null +++ b/figs/condensate16.fig/.energy16.gnuplot.swp diff --git a/figs/condensate16.fig/2020-05-22+16.dat b/figs/condensate16.fig/2020-05-22+16.dat new file mode 100644 index 0000000..c28285d --- /dev/null +++ b/figs/condensate16.fig/2020-05-22+16.dat @@ -0,0 +1,7 @@ +# rho E0 n0 + 1e-3 0.05195 0.61 + 5e-3 0.449 0.40 + 1e-2 1.1335 0.34 + 2e-2 2.748 0.33 + 1e-1 17.797 0.43 + 1e+0 197.35 0.70 diff --git a/figs/condensate16.fig/Makefile b/figs/condensate16.fig/Makefile new file mode 100644 index 0000000..d44f73c --- /dev/null +++ b/figs/condensate16.fig/Makefile @@ -0,0 +1,35 @@ +PROJECTNAME=condensate16 + +SIMPLEQ=~/Work/Research/2018+bose_gas/cmp/simpleq + +DATS=simpleq16.dat medeq16.dat +PDFS=$(addsuffix .pdf, $(PROJECTNAME)) +TEXS=$(addsuffix .tikz.tex, $(PROJECTNAME)) + +all: $(PDFS) + +$(PDFS): $(DATS) + gnuplot $(patsubst %.pdf, %.gnuplot, $@) > $(patsubst %.pdf, %.tikz.tex, $@) + pdflatex -jobname $(basename $@) -file-line-error $(patsubst %.pdf, %.tikz.tex, $@) + +simpleq16.dat: + julia $(SIMPLEQ)/main.jl -p "tolerance=1e-11;order=100;maxiter=100;v_a=16" -M simpleq condensate_fraction > $@ +medeq16.dat: + julia -p 8 $(SIMPLEQ)/main.jl -p "aK=1;aL1=1;aL2=1;aL3=0;bK=1;bL1=1;bL2=1;gK=1;gL2=1;gL3=0;N=12;P=8;v_a=16" -M medeq condensate_fraction > $@ + + +install: $(PDFS) + cp $^ $(INSTALLDIR)/ + +clean-aux: + rm -f $(addsuffix .tikz.tex, $(PROJECTNAME)) + rm -f $(addsuffix .aux, $(PROJECTNAME)) + rm -f $(addsuffix .log, $(PROJECTNAME)) + +clean-dat: + rm -f $(DATS) + +clean-tex: + rm -f $(PDFS) + +clean: clean-aux clean-tex diff --git a/figs/condensate16.fig/condensate16.gnuplot b/figs/condensate16.fig/condensate16.gnuplot new file mode 100644 index 0000000..e94eb1a --- /dev/null +++ b/figs/condensate16.fig/condensate16.gnuplot @@ -0,0 +1,34 @@ +set ylabel "$\\eta$" norotate offset -1,0 +set xlabel "$\\rho$" + +set xtics 1e-6, 100, 100 +set xtics add ("$10^{-6}$" 0.000001, "$10^{-4}$" 0.0001, "$10^{-2}$" 0.01, "$1$" 1.0, "$10^2$" 100) +unset mxtics + +#set ytics 0.6, 0.1 +set mytics 2 + +set xrange [0.000001:100] + +# default output canvas size: 12.5cm x 8.75cm +set term lua tikz size 8,6 standalone + +set key off + + +# set linestyle +set style line 1 linetype rgbcolor "#4169E1" linewidth 2 +set style line 2 linetype rgbcolor "#DC143C" linewidth 2 +set style line 3 linetype rgbcolor "#32CD32" linewidth 2 +set style line 4 linetype rgbcolor "#4B0082" linewidth 2 +set style line 5 linetype rgbcolor "#DAA520" linewidth 2 + +set pointsize 1 + +set logscale x + +plot \ + "simpleq16.dat" using 1:($3<1e-5 ? 1-$2 : 1/0) with lines ls 1 ,\ + "medeq16.dat" using 1:($3<1e-5 ? 1-$2 : 1/0) with lines ls 4 ,\ + "2020-05-22+16.dat" using 1:3 with points ls 2 + diff --git a/figs/condensate16.fig/medeq16.dat b/figs/condensate16.fig/medeq16.dat new file mode 100644 index 0000000..16ef72d --- /dev/null +++ b/figs/condensate16.fig/medeq16.dat @@ -0,0 +1,100 @@ + 1.000000000000000e-06 1.166000749536307e-02 1.397500557749096e-14 + 1.202264434617413e-06 1.294150659557208e-02 1.638378932477403e-14 + 1.445439770745928e-06 1.433395863576837e-02 2.068669778175053e-14 + 1.737800828749376e-06 1.585523520080719e-02 9.413457244472722e-15 + 2.089296130854041e-06 1.752468227078358e-02 2.147730951610164e-14 + 2.511886431509582e-06 1.936219857512355e-02 5.079185685242370e-14 + 3.019951720402019e-06 2.138821025238190e-02 4.890613903920614e-14 + 3.630780547701017e-06 2.362406333962305e-02 3.239855276438957e-14 + 4.365158322401657e-06 2.609258293662487e-02 6.114316674680731e-14 + 5.248074602497723e-06 2.881866901219256e-02 2.406593746021106e-14 + 6.309573444801930e-06 3.182986273284136e-02 5.235506117228457e-14 + 7.585775750291836e-06 3.515686411604940e-02 2.215175470514193e-14 + 9.120108393559096e-06 3.883401898804675e-02 1.629554526288622e-14 + 1.096478196143185e-05 4.289981143876489e-02 4.459384031826223e-15 + 1.318256738556407e-05 4.739739018972533e-02 7.959533672028488e-14 + 1.584893192461114e-05 5.237512938256743e-02 1.856627535995191e-14 + 1.905460717963248e-05 5.788719492270885e-02 1.690345207320381e-14 + 2.290867652767775e-05 6.399407839575801e-02 8.321647126329251e-15 + 2.754228703338163e-05 7.076308123314866e-02 1.217647346338456e-14 + 3.311311214825908e-05 7.826876839555350e-02 8.937613010647126e-15 + 3.981071705534969e-05 8.659343186556372e-02 1.053125999750607e-14 + 4.786300923226380e-05 9.582758472268986e-02 2.683673474487137e-14 + 5.754399373371566e-05 1.060704517836634e-01 1.137789331779729e-14 + 6.918309709189363e-05 1.174303688455619e-01 2.458170758492918e-14 + 8.317637711026709e-05 1.300249907645473e-01 3.251467169362365e-14 + 1.000000000000000e-04 1.439812489667652e-01 1.604790218801148e-14 + 1.202264434617413e-04 1.594350485585813e-01 1.208705427621754e-14 + 1.445439770745928e-04 1.765306819395882e-01 9.949882875987985e-15 + 1.737800828749376e-04 1.954198298344809e-01 9.175378143390922e-15 + 2.089296130854039e-04 2.162598861818270e-01 2.214874949422118e-14 + 2.511886431509580e-04 2.392112901954234e-01 8.309109937206045e-14 + 3.019951720402016e-04 2.644336112806493e-01 2.396263674278965e-13 + 3.630780547701014e-04 2.920802203955168e-01 7.677444642221211e-13 + 4.365158322401661e-04 3.222913892466672e-01 2.316125381119014e-12 + 5.248074602497728e-04 3.551856019023070e-01 6.550410386633497e-12 + 6.309573444801930e-04 3.908488535122212e-01 1.868355419086652e-15 + 7.585775750291835e-04 4.293218114670888e-01 2.423897582964416e-15 + 9.120108393559096e-04 4.705848954147140e-01 6.448772157311636e-15 + 1.096478196143185e-03 5.145415573312192e-01 5.163654832442502e-15 + 1.318256738556408e-03 5.610003206868235e-01 9.938237834818324e-15 + 1.584893192461114e-03 6.096564869738526e-01 1.561488943107384e-14 + 1.905460717963246e-03 6.600748309932439e-01 1.480152300205420e-14 + 2.290867652767772e-03 7.116750527032143e-01 2.186407540038661e-14 + 2.754228703338166e-03 7.637221822124044e-01 8.674104990124533e-15 + 3.311311214825911e-03 8.153244713068495e-01 3.829444151273552e-14 + 3.981071705534973e-03 8.654414602027135e-01 1.538039001478893e-14 + 4.786300923226385e-03 9.129047921304377e-01 6.821521006308993e-15 + 5.754399373371567e-03 9.564538682725248e-01 2.800593843254869e-15 + 6.918309709189363e-03 9.947874470918988e-01 4.113794640712594e-15 + 8.317637711026709e-03 1.026630567329738e+00 3.133249774936964e-15 + 1.000000000000000e-02 1.050813611431618e+00 4.457911153045777e-15 + 1.202264434617413e-02 1.066357266225810e+00 6.422681440463997e-15 + 1.445439770745928e-02 1.072554343989952e+00 4.370320881832804e-15 + 1.737800828749376e-02 1.069037933552943e+00 3.253121018609543e-15 + 2.089296130854041e-02 1.055825877502355e+00 3.370516638785516e-15 + 2.511886431509582e-02 1.033334096440169e+00 5.365630506272522e-15 + 3.019951720402019e-02 1.002355180715882e+00 3.715496417945286e-15 + 3.630780547701010e-02 9.640034053554465e-01 3.496056342568733e-15 + 4.365158322401656e-02 9.196325847132431e-01 1.506892881848512e-15 + 5.248074602497723e-02 8.707370299052170e-01 1.746468145174605e-15 + 6.309573444801930e-02 8.188475009356775e-01 1.661095648248939e-15 + 7.585775750291836e-02 7.654340485028318e-01 2.123676640240517e-15 + 9.120108393559097e-02 7.118255460316097e-01 4.350828619157266e-12 + 1.096478196143185e-01 6.591523262981006e-01 1.196771616697282e-12 + 1.318256738556407e-01 6.083146270922432e-01 3.018342123439392e-13 + 1.584893192461114e-01 5.599756911772029e-01 6.985575344309217e-14 + 1.905460717963248e-01 5.145752959414555e-01 2.023128031947868e-14 + 2.290867652767775e-01 4.723576806643551e-01 7.636489971812963e-14 + 2.754228703338169e-01 4.334072961327285e-01 3.450372875056261e-14 + 3.311311214825908e-01 3.976865103740502e-01 1.050825065430309e-14 + 3.981071705534969e-01 3.650709469150959e-01 1.016139346429647e-14 + 4.786300923226380e-01 3.353798855663381e-01 6.058578163892193e-15 + 5.754399373371566e-01 3.084006551579872e-01 3.382728454855092e-15 + 6.918309709189363e-01 2.839070252124540e-01 2.709608246018857e-15 + 8.317637711026709e-01 2.616722559700707e-01 1.878556843412201e-15 + 1.000000000000000e+00 2.414777667518989e-01 1.223063234010325e-15 + 1.202264434617413e+00 2.231184278236290e-01 1.993050434866885e-15 + 1.445439770745928e+00 2.064053790451102e-01 3.043149243838522e-15 + 1.737800828749376e+00 1.911671203240175e-01 2.998782498464775e-15 + 2.089296130854041e+00 1.772494557593911e-01 2.611515016643540e-15 + 2.511886431509582e+00 1.645147256322969e-01 1.646713357943204e-15 + 3.019951720402019e+00 1.528406355477221e-01 2.632388911909736e-15 + 3.630780547701010e+00 1.421188932690029e-01 3.725804901562014e-15 + 4.365158322401657e+00 1.322537904924787e-01 3.933266540669057e-15 + 5.248074602497723e+00 1.231608149876170e-01 3.818189041855540e-15 + 6.309573444801930e+00 1.147653429926582e-01 3.557698468765215e-15 + 7.585775750291836e+00 1.070014379973335e-01 3.200776178020959e-15 + 9.120108393559097e+00 9.981076672752043e-02 2.768978034961086e-15 + 1.096478196143185e+01 9.314163373013455e-02 2.784644153002823e-15 + 1.318256738556407e+01 8.694813025242200e-02 2.561848808069188e-15 + 1.584893192461114e+01 8.118938938383906e-02 3.342103900653102e-15 + 1.905460717963248e+01 7.582893689428576e-02 3.061126283356729e-15 + 2.290867652767775e+01 7.083412612140300e-02 3.012741771786859e-15 + 2.754228703338169e+01 6.617564611549075e-02 3.540875919778527e-15 + 3.311311214825908e+01 6.182709462529966e-02 3.384458193003484e-15 + 3.981071705534969e+01 5.776460991462094e-02 3.313778691885289e-15 + 4.786300923226381e+01 5.396655640631563e-02 4.116261947533335e-15 + 5.754399373371567e+01 5.041325862799244e-02 3.799276533234762e-15 + 6.918309709189363e+01 4.708677706941806e-02 4.136504038005983e-15 + 8.317637711026708e+01 4.397071963615008e-02 4.852059592739983e-15 diff --git a/figs/condensate16.fig/simpleq16.dat b/figs/condensate16.fig/simpleq16.dat new file mode 100644 index 0000000..c19b5fd --- /dev/null +++ b/figs/condensate16.fig/simpleq16.dat @@ -0,0 +1,100 @@ + 1.000000000000000e-06 1.193252238402537e-02 1.001128563083394e-13 + 1.202264434617413e-06 1.310706858781873e-02 1.166419917451259e-13 + 1.445439770745928e-06 1.439955124980759e-02 1.641671221071265e-13 + 1.737800828749376e-06 1.582224960123280e-02 1.318046204468745e-13 + 2.089296130854041e-06 1.738879830343562e-02 6.465028107789633e-14 + 2.511886431509582e-06 1.911434871723211e-02 1.282631449030218e-13 + 3.019951720402019e-06 2.101575037740115e-02 1.828007544593001e-13 + 3.630780547701017e-06 2.311175508419431e-02 8.716203682159707e-14 + 4.365158322401657e-06 2.542324601386885e-02 2.567118665713670e-13 + 5.248074602497723e-06 2.797349426852578e-02 1.909114183041957e-13 + 6.309573444801930e-06 3.078844511607499e-02 4.832552648271436e-13 + 7.585775750291836e-06 3.389703567082637e-02 9.663351363715864e-13 + 9.120108393559096e-06 3.733154490632996e-02 1.532273658747508e-12 + 1.096478196143185e-05 4.112797535363494e-02 2.557647979444216e-12 + 1.318256738556407e-05 4.532646335311902e-02 4.227854029289540e-12 + 1.584893192461114e-05 4.997171084390572e-02 7.262584516187726e-12 + 1.905460717963248e-05 5.511342573964625e-02 1.581480833654367e-13 + 2.290867652767775e-05 6.080674905753369e-02 1.081736075443477e-13 + 2.754228703338163e-05 6.711263388139645e-02 9.892826615755524e-14 + 3.311311214825908e-05 7.409812223095856e-02 1.000146289881600e-13 + 3.981071705534969e-05 8.183643867649469e-02 1.728836305265867e-13 + 4.786300923226380e-05 9.040678111314078e-02 4.874749473980382e-14 + 5.754399373371566e-05 9.989363584872350e-02 7.675596130065904e-14 + 6.918309709189363e-05 1.103853719968644e-01 3.488065867913569e-14 + 8.317637711026709e-05 1.219717753838641e-01 4.074847689970330e-13 + 1.000000000000000e-04 1.347400631725403e-01 8.837677779017190e-14 + 1.202264434617413e-04 1.487687813292578e-01 1.141726470797567e-13 + 1.445439770745928e-04 1.641188444089739e-01 2.367473489955975e-13 + 1.737800828749376e-04 1.808208708129144e-01 9.521496171923438e-14 + 2.089296130854039e-04 1.988579762291830e-01 2.852315979737347e-13 + 2.511886431509580e-04 2.181434521377151e-01 8.027724386154160e-13 + 3.019951720402016e-04 2.384934863352307e-01 1.058787720274644e-13 + 3.630780547701014e-04 2.595965448975509e-01 1.539765713616111e-13 + 4.365158322401661e-04 2.809834568202456e-01 6.040763684158786e-14 + 5.248074602497728e-04 3.020055785198434e-01 8.733135409851157e-13 + 6.309573444801930e-04 3.218319564612944e-01 3.684231597353426e-14 + 7.585775750291835e-04 3.394783780636358e-01 1.647277211646135e-13 + 9.120108393559096e-04 3.538788243847390e-01 2.032722085217800e-13 + 1.096478196143185e-03 3.640005594277430e-01 2.838665806426261e-13 + 1.318256738556408e-03 3.689881894227810e-01 9.260531251903156e-14 + 1.584893192461114e-03 3.683050651268248e-01 6.326291032696907e-14 + 1.905460717963246e-03 3.618322423397615e-01 4.633167297377464e-14 + 2.290867652767772e-03 3.498937848449674e-01 1.658072106704119e-13 + 2.754228703338166e-03 3.332009525742046e-01 1.122897591215287e-13 + 3.311311214825911e-03 3.127348598651754e-01 5.692534927020893e-14 + 3.981071705534973e-03 2.896035594654576e-01 5.657153215877982e-14 + 4.786300923226385e-03 2.649092240010952e-01 9.048650405950154e-14 + 5.754399373371567e-03 2.396485623749096e-01 8.607211052131642e-14 + 6.918309709189363e-03 2.146539478886141e-01 1.595875861854557e-13 + 8.317637711026709e-03 1.905709675646656e-01 1.531883433865405e-13 + 1.000000000000000e-02 1.678624919366197e-01 1.407687268306203e-13 + 1.202264434617413e-02 1.468287025622779e-01 6.748557036716531e-14 + 1.445439770745928e-02 1.276344869940322e-01 6.869589239332279e-14 + 1.737800828749376e-02 1.103383284494799e-01 2.437733106567872e-14 + 2.089296130854041e-02 9.491924304319868e-02 1.064344070654439e-13 + 2.511886431509582e-02 8.130009294494668e-02 9.238668475341097e-15 + 3.019951720402019e-02 6.936674066697997e-02 1.020232616470553e-13 + 3.630780547701010e-02 5.898315634460066e-02 5.074743271963278e-14 + 4.365158322401656e-02 5.000290867811857e-02 1.028088279693266e-13 + 5.248074602497723e-02 4.227759058635282e-02 4.624881313939891e-12 + 6.309573444801930e-02 3.566274140971409e-02 1.852661419803306e-12 + 7.585775750291836e-02 3.002178399074026e-02 7.258972029472108e-13 + 9.120108393559097e-02 2.522842899157464e-02 2.891235204847115e-13 + 1.096478196143185e-01 2.116792739912998e-02 1.401726396203463e-13 + 1.318256738556407e-01 1.773748368422931e-02 5.743044386657572e-14 + 1.584893192461114e-01 1.484608018800232e-02 3.606825259935793e-14 + 1.905460717963248e-01 1.241390929900180e-02 1.892505250290776e-13 + 2.290867652767775e-01 1.037156402191131e-02 2.083768656818398e-14 + 2.754228703338169e-01 8.659099268369731e-03 4.075384057289656e-14 + 3.311311214825908e-01 7.225045018893100e-03 2.600299498350245e-14 + 3.981071705534969e-01 6.025427676628160e-03 1.413936909490012e-14 + 4.786300923226380e-01 5.022836651218588e-03 2.463069066912327e-14 + 5.754399373371566e-01 4.185558640329008e-03 8.558078927823013e-15 + 6.918309709189363e-01 3.486791398547768e-03 6.226596123328673e-14 + 8.317637711026709e-01 2.903941231756474e-03 7.123438110270025e-14 + 1.000000000000000e+00 2.418003342628212e-03 8.748669232588792e-14 + 1.202264434617413e+00 2.013020888682947e-03 6.564186088832174e-14 + 1.445439770745928e+00 1.675616708432474e-03 6.210366709884162e-14 + 1.737800828749376e+00 1.394590731209103e-03 5.060423878974540e-15 + 2.089296130854041e+00 1.160575807795245e-03 4.504831194472809e-15 + 2.511886431509582e+00 9.657448526444926e-04 5.178886653500039e-14 + 3.019951720402019e+00 8.035626066249569e-04 5.546884105120939e-14 + 3.630780547701010e+00 6.685758928573441e-04 4.099740962135071e-15 + 4.365158322401657e+00 5.562368657074897e-04 5.461022666145280e-14 + 5.248074602497723e+00 4.627543904435245e-04 2.773788643238822e-15 + 6.309573444801930e+00 3.849693047135169e-04 5.605531588703275e-14 + 7.585775750291836e+00 3.202498832578166e-04 1.548482473273102e-15 + 9.120108393559097e+00 2.664043442730844e-04 5.542616472939282e-15 + 1.096478196143185e+01 2.216076961608340e-04 5.416841008960720e-15 + 1.318256738556407e+01 1.843406276770280e-04 1.051948766670258e-13 + 1.584893192461114e+01 1.533384958473269e-04 3.020416191025176e-15 + 1.905460717963248e+01 1.275487688416503e-04 1.054514720940622e-13 + 2.290867652767775e+01 1.060955403032162e-04 1.054818735884649e-13 + 2.754228703338169e+01 8.824995249206293e-05 1.068333765010024e-13 + 3.311311214825908e+01 7.340555293075552e-05 5.311424682002792e-14 + 3.981071705534969e+01 6.105776757545735e-05 8.573716686070972e-15 + 4.786300923226381e+01 5.078680698963938e-05 5.255692893479381e-14 + 5.754399373371567e+01 4.224343422123316e-05 5.412915990778012e-16 + 6.918309709189363e+01 3.513711727611226e-05 5.629397185139020e-14 + 8.317637711026708e+01 2.922616801548849e-05 8.852629630882290e-12 diff --git a/figs/correlation.fig/2020-08-27+correlation.dat b/figs/correlation.fig/2020-08-27+correlation.dat new file mode 100644 index 0000000..c436b89 --- /dev/null +++ b/figs/correlation.fig/2020-08-27+correlation.dat @@ -0,0 +1,102 @@ +# Potential: 16 e^{-|x|}, rho=0.02 +# |x| 2-pt correlation + 0.04641589 0.00532330 0.00040199 0.00532330 0.00040199 0.00000000 0.00000000 0.00000000 0.00000000 + 0.13924767 0.00537264 0.00024875 0.00537264 0.00024875 0.00000000 0.00000000 0.00000000 0.00000000 + 0.23207944 0.00605816 0.00017567 0.00605816 0.00017567 0.00000000 0.00000000 0.00000000 0.00000000 + 0.32491122 0.00727112 0.00014500 0.00727112 0.00014500 0.00000000 0.00000000 0.00000000 0.00000000 + 0.41774300 0.00896049 0.00013083 0.00896049 0.00013083 0.00000000 0.00000000 0.00000000 0.00000000 + 0.51057477 0.01116990 0.00012069 0.01116990 0.00012069 0.00000000 0.00000000 0.00000000 0.00000000 + 0.60340655 0.01418247 0.00011223 0.01418247 0.00011223 0.00000000 0.00000000 0.00000000 0.00000000 + 0.69623833 0.01828721 0.00011014 0.01828721 0.00011014 0.00000000 0.00000000 0.00000000 0.00000000 + 0.78907010 0.02370094 0.00011191 0.02370094 0.00011191 0.00000000 0.00000000 0.00000000 0.00000000 + 0.88190188 0.03049612 0.00012453 0.03049612 0.00012453 0.00000000 0.00000000 0.00000000 0.00000000 + 0.97473366 0.03907954 0.00012112 0.03907954 0.00012112 0.00000000 0.00000000 0.00000000 0.00000000 + 1.06756543 0.04977163 0.00012676 0.04977163 0.00012676 0.00000000 0.00000000 0.00000000 0.00000000 + 1.16039721 0.06291596 0.00012880 0.06291596 0.00012880 0.00000000 0.00000000 0.00000000 0.00000000 + 1.25322899 0.07860896 0.00013522 0.07860896 0.00013522 0.00000000 0.00000000 0.00000000 0.00000000 + 1.34606076 0.09733963 0.00014306 0.09733963 0.00014306 0.00000000 0.00000000 0.00000000 0.00000000 + 1.43889254 0.11959686 0.00015342 0.11959686 0.00015342 0.00000000 0.00000000 0.00000000 0.00000000 + 1.53172432 0.14519143 0.00016228 0.14519143 0.00016228 0.00000000 0.00000000 0.00000000 0.00000000 + 1.62455609 0.17429949 0.00016887 0.17429949 0.00016887 0.00000000 0.00000000 0.00000000 0.00000000 + 1.71738787 0.20718710 0.00017469 0.20718710 0.00017469 0.00000000 0.00000000 0.00000000 0.00000000 + 1.81021965 0.24369539 0.00017395 0.24369539 0.00017395 0.00000000 0.00000000 0.00000000 0.00000000 + 1.90305142 0.28352716 0.00017471 0.28352716 0.00017471 0.00000000 0.00000000 0.00000000 0.00000000 + 1.99588320 0.32631191 0.00017370 0.32631191 0.00017370 0.00000000 0.00000000 0.00000000 0.00000000 + 2.08871498 0.37151918 0.00017082 0.37151918 0.00017082 0.00000000 0.00000000 0.00000000 0.00000000 + 2.18154675 0.41909276 0.00017134 0.41909276 0.00017134 0.00000000 0.00000000 0.00000000 0.00000000 + 2.27437853 0.46853305 0.00017701 0.46853305 0.00017701 0.00000000 0.00000000 0.00000000 0.00000000 + 2.36721031 0.51895440 0.00017445 0.51895440 0.00017445 0.00000000 0.00000000 0.00000000 0.00000000 + 2.46004208 0.56983458 0.00017914 0.56983458 0.00017914 0.00000000 0.00000000 0.00000000 0.00000000 + 2.55287386 0.62086170 0.00018145 0.62086170 0.00018145 0.00000000 0.00000000 0.00000000 0.00000000 + 2.64570564 0.67087338 0.00017118 0.67087338 0.00017118 0.00000000 0.00000000 0.00000000 0.00000000 + 2.73853741 0.71925531 0.00017770 0.71925531 0.00017770 0.00000000 0.00000000 0.00000000 0.00000000 + 2.83136919 0.76584767 0.00016606 0.76584767 0.00016606 0.00000000 0.00000000 0.00000000 0.00000000 + 2.92420097 0.80980633 0.00016880 0.80980633 0.00016880 0.00000000 0.00000000 0.00000000 0.00000000 + 3.01703274 0.85064291 0.00017081 0.85064291 0.00017081 0.00000000 0.00000000 0.00000000 0.00000000 + 3.10986452 0.88811443 0.00016890 0.88811443 0.00016890 0.00000000 0.00000000 0.00000000 0.00000000 + 3.20269630 0.92186331 0.00017448 0.92186331 0.00017448 0.00000000 0.00000000 0.00000000 0.00000000 + 3.29552807 0.95179423 0.00017224 0.95179423 0.00017224 0.00000000 0.00000000 0.00000000 0.00000000 + 3.38835985 0.97798233 0.00016981 0.97798233 0.00016981 0.00000000 0.00000000 0.00000000 0.00000000 + 3.48119163 1.00058535 0.00016426 1.00058535 0.00016426 0.00000000 0.00000000 0.00000000 0.00000000 + 3.57402340 1.01951391 0.00016349 1.01951391 0.00016349 0.00000000 0.00000000 0.00000000 0.00000000 + 3.66685518 1.03469645 0.00015961 1.03469645 0.00015961 0.00000000 0.00000000 0.00000000 0.00000000 + 3.75968696 1.04648197 0.00015142 1.04648197 0.00015142 0.00000000 0.00000000 0.00000000 0.00000000 + 3.85251873 1.05521888 0.00014825 1.05521888 0.00014825 0.00000000 0.00000000 0.00000000 0.00000000 + 3.94535051 1.06104242 0.00015316 1.06104242 0.00015316 0.00000000 0.00000000 0.00000000 0.00000000 + 4.03818229 1.06451223 0.00015174 1.06451223 0.00015174 0.00000000 0.00000000 0.00000000 0.00000000 + 4.13101406 1.06585125 0.00014890 1.06585125 0.00014890 0.00000000 0.00000000 0.00000000 0.00000000 + 4.22384584 1.06532590 0.00014168 1.06532590 0.00014168 0.00000000 0.00000000 0.00000000 0.00000000 + 4.31667762 1.06332074 0.00014668 1.06332074 0.00014668 0.00000000 0.00000000 0.00000000 0.00000000 + 4.40950939 1.05993711 0.00014682 1.05993711 0.00014682 0.00000000 0.00000000 0.00000000 0.00000000 + 4.50234117 1.05561851 0.00014625 1.05561851 0.00014625 0.00000000 0.00000000 0.00000000 0.00000000 + 4.59517295 1.05060642 0.00013572 1.05060642 0.00013572 0.00000000 0.00000000 0.00000000 0.00000000 + 4.68800472 1.04534502 0.00013193 1.04534502 0.00013193 0.00000000 0.00000000 0.00000000 0.00000000 + 4.78083650 1.03971938 0.00012073 1.03971938 0.00012073 0.00000000 0.00000000 0.00000000 0.00000000 + 4.87366828 1.03396128 0.00012257 1.03396128 0.00012257 0.00000000 0.00000000 0.00000000 0.00000000 + 4.96650005 1.02834322 0.00011721 1.02834322 0.00011721 0.00000000 0.00000000 0.00000000 0.00000000 + 5.05933183 1.02289857 0.00012010 1.02289857 0.00012010 0.00000000 0.00000000 0.00000000 0.00000000 + 5.15216361 1.01760580 0.00012241 1.01760580 0.00012241 0.00000000 0.00000000 0.00000000 0.00000000 + 5.24499538 1.01265617 0.00011612 1.01265617 0.00011612 0.00000000 0.00000000 0.00000000 0.00000000 + 5.33782716 1.00816416 0.00010902 1.00816416 0.00010902 0.00000000 0.00000000 0.00000000 0.00000000 + 5.43065894 1.00407768 0.00010997 1.00407768 0.00010997 0.00000000 0.00000000 0.00000000 0.00000000 + 5.52349071 1.00065299 0.00010948 1.00065299 0.00010948 0.00000000 0.00000000 0.00000000 0.00000000 + 5.61632249 0.99773671 0.00011169 0.99773671 0.00011169 0.00000000 0.00000000 0.00000000 0.00000000 + 5.70915427 0.99521455 0.00011022 0.99521455 0.00011022 0.00000000 0.00000000 0.00000000 0.00000000 + 5.80198604 0.99313289 0.00010491 0.99313289 0.00010491 0.00000000 0.00000000 0.00000000 0.00000000 + 5.89481782 0.99149782 0.00010175 0.99149782 0.00010175 0.00000000 0.00000000 0.00000000 0.00000000 + 5.98764960 0.99032198 0.00009790 0.99032198 0.00009790 0.00000000 0.00000000 0.00000000 0.00000000 + 6.08048137 0.98958110 0.00009610 0.98958110 0.00009610 0.00000000 0.00000000 0.00000000 0.00000000 + 6.17331315 0.98910112 0.00009612 0.98910112 0.00009612 0.00000000 0.00000000 0.00000000 0.00000000 + 6.26614493 0.98893555 0.00009629 0.98893555 0.00009629 0.00000000 0.00000000 0.00000000 0.00000000 + 6.35897670 0.98909456 0.00009693 0.98909456 0.00009693 0.00000000 0.00000000 0.00000000 0.00000000 + 6.45180848 0.98948384 0.00009719 0.98948384 0.00009719 0.00000000 0.00000000 0.00000000 0.00000000 + 6.54464026 0.99007466 0.00009778 0.99007466 0.00009778 0.00000000 0.00000000 0.00000000 0.00000000 + 6.63747203 0.99069454 0.00009175 0.99069454 0.00009175 0.00000000 0.00000000 0.00000000 0.00000000 + 6.73030381 0.99152323 0.00008887 0.99152323 0.00008887 0.00000000 0.00000000 0.00000000 0.00000000 + 6.82313559 0.99218652 0.00008729 0.99218652 0.00008729 0.00000000 0.00000000 0.00000000 0.00000000 + 6.91596736 0.99305844 0.00008686 0.99305844 0.00008686 0.00000000 0.00000000 0.00000000 0.00000000 + 7.00879914 0.99402451 0.00008778 0.99402451 0.00008778 0.00000000 0.00000000 0.00000000 0.00000000 + 7.10163092 0.99501487 0.00008386 0.99501487 0.00008386 0.00000000 0.00000000 0.00000000 0.00000000 + 7.19446269 0.99598405 0.00007919 0.99598405 0.00007919 0.00000000 0.00000000 0.00000000 0.00000000 + 7.28729447 0.99683782 0.00008214 0.99683782 0.00008214 0.00000000 0.00000000 0.00000000 0.00000000 + 7.38012625 0.99779938 0.00007994 0.99779938 0.00007994 0.00000000 0.00000000 0.00000000 0.00000000 + 7.47295802 0.99856267 0.00008169 0.99856267 0.00008169 0.00000000 0.00000000 0.00000000 0.00000000 + 7.56578980 0.99930080 0.00008203 0.99930080 0.00008203 0.00000000 0.00000000 0.00000000 0.00000000 + 7.65862158 1.00001485 0.00008598 1.00001485 0.00008598 0.00000000 0.00000000 0.00000000 0.00000000 + 7.75145335 1.00049762 0.00008693 1.00049762 0.00008693 0.00000000 0.00000000 0.00000000 0.00000000 + 7.84428513 1.00095153 0.00008531 1.00095153 0.00008531 0.00000000 0.00000000 0.00000000 0.00000000 + 7.93711691 1.00133388 0.00008126 1.00133388 0.00008126 0.00000000 0.00000000 0.00000000 0.00000000 + 8.02994868 1.00152882 0.00007855 1.00152882 0.00007855 0.00000000 0.00000000 0.00000000 0.00000000 + 8.12278046 1.00163599 0.00007736 1.00163599 0.00007736 0.00000000 0.00000000 0.00000000 0.00000000 + 8.21561224 1.00175148 0.00007549 1.00175148 0.00007549 0.00000000 0.00000000 0.00000000 0.00000000 + 8.30844401 1.00172584 0.00007298 1.00172584 0.00007298 0.00000000 0.00000000 0.00000000 0.00000000 + 8.40127579 1.00175558 0.00007456 1.00175558 0.00007456 0.00000000 0.00000000 0.00000000 0.00000000 + 8.49410757 1.00175353 0.00007550 1.00175353 0.00007550 0.00000000 0.00000000 0.00000000 0.00000000 + 8.58693934 1.00167392 0.00007070 1.00167392 0.00007070 0.00000000 0.00000000 0.00000000 0.00000000 + 8.67977112 1.00158009 0.00006882 1.00158009 0.00006882 0.00000000 0.00000000 0.00000000 0.00000000 + 8.77260290 1.00144512 0.00006669 1.00144512 0.00006669 0.00000000 0.00000000 0.00000000 0.00000000 + 8.86543467 1.00130747 0.00006420 1.00130747 0.00006420 0.00000000 0.00000000 0.00000000 0.00000000 + 8.95826645 1.00117137 0.00006419 1.00117137 0.00006419 0.00000000 0.00000000 0.00000000 0.00000000 + 9.05109823 1.00109365 0.00006530 1.00109365 0.00006530 0.00000000 0.00000000 0.00000000 0.00000000 + 9.14393000 1.00094515 0.00006273 1.00094515 0.00006273 0.00000000 0.00000000 0.00000000 0.00000000 + 9.23676178 1.00076430 0.00006358 1.00076430 0.00006358 0.00000000 0.00000000 0.00000000 0.00000000 diff --git a/figs/correlation.fig/Makefile b/figs/correlation.fig/Makefile new file mode 100644 index 0000000..fd4b6c3 --- /dev/null +++ b/figs/correlation.fig/Makefile @@ -0,0 +1,36 @@ +PROJECTNAME=correlation + +SIMPLEQ=~/Work/Research/2018+bose_gas/cmp/simpleq + +DATS=bigeq.dat mueq.dat simpleq.dat +PDFS=$(addsuffix .pdf, $(PROJECTNAME)) +TEXS=$(addsuffix .tikz.tex, $(PROJECTNAME)) + +all: $(PDFS) + +$(PDFS): $(DATS) + gnuplot $(patsubst %.pdf, %.gnuplot, $@) > $(patsubst %.pdf, %.tikz.tex, $@) + pdflatex -jobname $(basename $@) -file-line-error $(patsubst %.pdf, %.tikz.tex, $@) + +simpleq.dat: + julia $(SIMPLEQ)/main.jl -p "rho=2e-2;minlrho=-6;nlrho=60;eq=simpleq;N=12;P=8;xmax=10" -M medeq 2pt-simpleq > $@ +mueq.dat: + julia $(SIMPLEQ)/main.jl -p "rho=2e-2;minlrho=-6;nlrho=60;eq=mueq;N=12;P=8;xmax=10" -M medeq 2pt > $@ +bigeq.dat: + julia $(SIMPLEQ)/main.jl -p "rho=2e-2;minlrho=-6;nlrho=60;v_a=16;eq=bigeq;N=12;P=8;xmax=10" -M medeq 2pt > $@ + +install: $(PDFS) + cp $^ $(INSTALLDIR)/ + +clean-aux: + rm -f $(addsuffix .tikz.tex, $(PROJECTNAME)) + rm -f $(addsuffix .aux, $(PROJECTNAME)) + rm -f $(addsuffix .log, $(PROJECTNAME)) + +clean-dat: + rm -f $(DATS) + +clean-tex: + rm -f $(PDFS) + +clean: clean-aux clean-tex diff --git a/figs/correlation.fig/bigeq.dat b/figs/correlation.fig/bigeq.dat new file mode 100644 index 0000000..07b50b3 --- /dev/null +++ b/figs/correlation.fig/bigeq.dat @@ -0,0 +1,100 @@ + 1.000000000000000e-01 -7.492541163311106e-06 + 2.000000000000000e-01 -8.016902978356378e-06 + 3.000000000000000e-01 -8.679554787819034e-06 + 4.000000000000000e-01 -9.411241330792030e-06 + 5.000000000000000e-01 -1.015424596606049e-05 + 6.000000000000000e-01 -1.087711166991134e-05 + 7.000000000000000e-01 -1.151894832197664e-05 + 8.000000000000000e-01 -1.194142990150118e-05 + 9.000000000000000e-01 -1.201566606206566e-05 + 1.000000000000000e+00 -1.155831226525352e-05 + 1.100000000000000e+00 -1.032473715914746e-05 + 1.200000000000000e+00 -8.112234090907045e-06 + 1.300000000000000e+00 -4.727643509494950e-06 + 1.400000000000000e+00 2.139117276462876e-08 + 1.500000000000000e+00 6.344851444408795e-06 + 1.600000000000000e+00 1.444909419071485e-05 + 1.700000000000000e+00 2.442974443258903e-05 + 1.800000000000000e+00 3.635217589273449e-05 + 1.900000000000000e+00 5.028207520646167e-05 + 2.000000000000000e+00 6.618981856920616e-05 + 2.100000000000000e+00 8.394946793146488e-05 + 2.200000000000000e+00 1.033337991622068e-04 + 2.300000000000000e+00 1.241127999309069e-04 + 2.400000000000000e+00 1.460739157581073e-04 + 2.500000000000000e+00 1.688969605587932e-04 + 2.600000000000000e+00 1.922338554815170e-04 + 2.700000000000000e+00 2.157320193004696e-04 + 2.800000000000000e+00 2.390057792097371e-04 + 2.900000000000000e+00 2.617182666847289e-04 + 3.000000000000000e+00 2.835764720628177e-04 + 3.100000000000000e+00 3.043594085466944e-04 + 3.200000000000000e+00 3.238609190134170e-04 + 3.300000000000000e+00 3.418163389612801e-04 + 3.400000000000000e+00 3.580645020981861e-04 + 3.500000000000000e+00 3.725706343056383e-04 + 3.600000000000000e+00 3.853364786559161e-04 + 3.700000000000000e+00 3.963643140045606e-04 + 3.800000000000000e+00 4.056288039025831e-04 + 3.900000000000000e+00 4.131965917181045e-04 + 4.000000000000000e+00 4.192087182163521e-04 + 4.100000000000000e+00 4.237754997457688e-04 + 4.200000000000000e+00 4.270396577532589e-04 + 4.300000000000000e+00 4.291417424125078e-04 + 4.400000000000000e+00 4.301990718613866e-04 + 4.500000000000000e+00 4.303462807997311e-04 + 4.600000000000000e+00 4.297231551355101e-04 + 4.700000000000000e+00 4.285116447871562e-04 + 4.800000000000000e+00 4.268572864513014e-04 + 4.900000000000000e+00 4.248196795555219e-04 + 5.000000000000000e+00 4.224865659633607e-04 + 5.100000000000000e+00 4.199715983241658e-04 + 5.200000000000000e+00 4.173825684782714e-04 + 5.300000000000000e+00 4.147934916281848e-04 + 5.400000000000000e+00 4.122311381082057e-04 + 5.500000000000000e+00 4.097622044961587e-04 + 5.600000000000000e+00 4.074459723215993e-04 + 5.700000000000000e+00 4.052871021874686e-04 + 5.800000000000000e+00 4.033031600836669e-04 + 5.900000000000000e+00 4.015105996408118e-04 + 6.000000000000000e+00 3.999275938626445e-04 + 6.100000000000000e+00 3.985668439209783e-04 + 6.200000000000000e+00 3.974177250112088e-04 + 6.300000000000000e+00 3.964879713397231e-04 + 6.400000000000000e+00 3.957576409649167e-04 + 6.500000000000000e+00 3.951823102306802e-04 + 6.600000000000000e+00 3.947646482049149e-04 + 6.700000000000000e+00 3.945221470410044e-04 + 6.800000000000000e+00 3.944573880683813e-04 + 6.900000000000000e+00 3.945250250489386e-04 + 7.000000000000000e+00 3.946577269589084e-04 + 7.100000000000000e+00 3.948557469149244e-04 + 7.200000000000000e+00 3.951364830948637e-04 + 7.300000000000000e+00 3.954938002278550e-04 + 7.400000000000000e+00 3.959164433740596e-04 + 7.500000000000000e+00 3.963714179081655e-04 + 7.600000000000000e+00 3.968325233124400e-04 + 7.700000000000000e+00 3.972813752787329e-04 + 7.800000000000000e+00 3.977122476507974e-04 + 7.900000000000000e+00 3.981522766684607e-04 + 8.000000000000000e+00 3.985956172957380e-04 + 8.100000000000000e+00 3.990075357529208e-04 + 8.199999999999999e+00 3.993746815454193e-04 + 8.300000000000001e+00 3.996997792152681e-04 + 8.400000000000000e+00 4.000003941078366e-04 + 8.500000000000000e+00 4.002682692032775e-04 + 8.600000000000000e+00 4.004855271470240e-04 + 8.699999999999999e+00 4.006719309576352e-04 + 8.800000000000001e+00 4.008334081453375e-04 + 8.900000000000000e+00 4.009580309304248e-04 + 9.000000000000000e+00 4.010411794697340e-04 + 9.100000000000000e+00 4.010836381495669e-04 + 9.199999999999999e+00 4.011041322566935e-04 + 9.300000000000001e+00 4.011069223510944e-04 + 9.400000000000000e+00 4.010881050199837e-04 + 9.500000000000000e+00 4.010582610257124e-04 + 9.600000000000000e+00 4.010069812998826e-04 + 9.699999999999999e+00 4.009266657300464e-04 + 9.800000000000001e+00 4.008319709508709e-04 + 9.900000000000000e+00 4.007423209355591e-04 + 1.000000000000000e+01 4.006713168605734e-04 diff --git a/figs/correlation.fig/correlation.gnuplot b/figs/correlation.fig/correlation.gnuplot new file mode 100644 index 0000000..9be226f --- /dev/null +++ b/figs/correlation.fig/correlation.gnuplot @@ -0,0 +1,26 @@ +set ylabel "$\\frac{C_2}{\\rho^2}$" norotate offset -1,0 +set xlabel "$|\\mathbf x|$" + +# default output canvas size: 12.5cm x 8.75cm +set term lua tikz size 8,6 standalone + +unset key + + +# set linestyle +set style line 1 linetype rgbcolor "#4169E1" linewidth 2 +set style line 2 linetype rgbcolor "#DC143C" linewidth 1 +set style line 3 linetype rgbcolor "#32CD32" linewidth 2 +set style line 4 linetype rgbcolor "#4B0082" linewidth 2 +set style line 5 linetype rgbcolor "#DAA520" linewidth 2 +set style line 6 linetype rgbcolor "#555500" linewidth 2 + +set pointsize 0.5 + +rho=2e-2 + +plot \ + "simpleq.dat" using 1:($2/rho**2) with lines ls 1,\ + "bigeq.dat" using 1:($2/rho**2) with lines ls 4,\ + "2020-08-27+correlation.dat" using 1:2 ls 2 + diff --git a/figs/correlation.fig/mueq.dat b/figs/correlation.fig/mueq.dat new file mode 100644 index 0000000..4662058 --- /dev/null +++ b/figs/correlation.fig/mueq.dat @@ -0,0 +1,100 @@ + 1.000000000000000e-01 2.130418403340800e-04 + 2.000000000000000e-01 2.145091238125419e-04 + 3.000000000000000e-01 2.168961803510115e-04 + 4.000000000000000e-01 2.200881047098373e-04 + 5.000000000000000e-01 2.239499052926825e-04 + 6.000000000000000e-01 2.283477782850422e-04 + 7.000000000000000e-01 2.331717114746654e-04 + 8.000000000000000e-01 2.383427545804145e-04 + 9.000000000000000e-01 2.437959663709728e-04 + 1.000000000000000e+00 2.494792780364908e-04 + 1.100000000000000e+00 2.553458968242682e-04 + 1.200000000000000e+00 2.613430666435472e-04 + 1.300000000000000e+00 2.674185831682237e-04 + 1.400000000000000e+00 2.735256018683541e-04 + 1.500000000000000e+00 2.796262204436459e-04 + 1.600000000000000e+00 2.856899925541691e-04 + 1.700000000000000e+00 2.916874047060339e-04 + 1.800000000000000e+00 2.975929604366267e-04 + 1.900000000000000e+00 3.033849364510787e-04 + 2.000000000000000e+00 3.090415725152513e-04 + 2.100000000000000e+00 3.145428049252462e-04 + 2.200000000000000e+00 3.198715027679640e-04 + 2.300000000000000e+00 3.250164750188606e-04 + 2.400000000000000e+00 3.299712599259801e-04 + 2.500000000000000e+00 3.347290403741292e-04 + 2.600000000000000e+00 3.392840318356562e-04 + 2.700000000000000e+00 3.436317387754585e-04 + 2.800000000000000e+00 3.477685463904850e-04 + 2.900000000000000e+00 3.516938041693960e-04 + 3.000000000000000e+00 3.554093745145367e-04 + 3.100000000000000e+00 3.589193956234003e-04 + 3.200000000000000e+00 3.622282098184321e-04 + 3.300000000000000e+00 3.653385452150326e-04 + 3.400000000000000e+00 3.682546090447632e-04 + 3.500000000000000e+00 3.709826773383462e-04 + 3.600000000000000e+00 3.735297109279578e-04 + 3.700000000000000e+00 3.759027561267685e-04 + 3.800000000000000e+00 3.781082210370979e-04 + 3.900000000000000e+00 3.801535126992011e-04 + 4.000000000000000e+00 3.820465394300752e-04 + 4.100000000000000e+00 3.837942469843196e-04 + 4.200000000000000e+00 3.854041870302811e-04 + 4.300000000000000e+00 3.868843482990351e-04 + 4.400000000000000e+00 3.882427613658928e-04 + 4.500000000000000e+00 3.894876304891240e-04 + 4.600000000000000e+00 3.906264757819791e-04 + 4.700000000000000e+00 3.916663187842789e-04 + 4.800000000000000e+00 3.926125919231407e-04 + 4.900000000000000e+00 3.934690910984855e-04 + 5.000000000000000e+00 3.942407679577102e-04 + 5.100000000000000e+00 3.949344046916439e-04 + 5.200000000000000e+00 3.955581225640474e-04 + 5.300000000000000e+00 3.961201895403179e-04 + 5.400000000000000e+00 3.966275873395823e-04 + 5.500000000000000e+00 3.970862020080399e-04 + 5.600000000000000e+00 3.974994348438934e-04 + 5.700000000000000e+00 3.978677396766229e-04 + 5.800000000000000e+00 3.981909364370117e-04 + 5.900000000000000e+00 3.984697205807291e-04 + 6.000000000000000e+00 3.987072113448590e-04 + 6.100000000000000e+00 3.989094136349676e-04 + 6.200000000000000e+00 3.990841977712767e-04 + 6.300000000000000e+00 3.992400717207986e-04 + 6.400000000000000e+00 3.993830662220494e-04 + 6.500000000000000e+00 3.995150034938932e-04 + 6.600000000000000e+00 3.996343573920134e-04 + 6.700000000000000e+00 3.997371482013122e-04 + 6.800000000000000e+00 3.998190698185217e-04 + 6.900000000000000e+00 3.998776637317933e-04 + 7.000000000000000e+00 3.999143485918603e-04 + 7.100000000000000e+00 3.999356186244626e-04 + 7.200000000000000e+00 3.999501739938431e-04 + 7.300000000000000e+00 3.999653908411677e-04 + 7.400000000000000e+00 3.999852592754089e-04 + 7.500000000000000e+00 4.000092335895207e-04 + 7.600000000000000e+00 4.000335522452741e-04 + 7.700000000000000e+00 4.000532788642888e-04 + 7.800000000000000e+00 4.000645632283129e-04 + 7.900000000000000e+00 4.000663594103156e-04 + 8.000000000000000e+00 4.000595166185636e-04 + 8.100000000000000e+00 4.000460784804748e-04 + 8.199999999999999e+00 4.000290312335042e-04 + 8.300000000000001e+00 4.000113855708080e-04 + 8.400000000000000e+00 3.999957999911112e-04 + 8.500000000000000e+00 3.999838939317796e-04 + 8.600000000000000e+00 3.999765161553756e-04 + 8.699999999999999e+00 3.999743111385341e-04 + 8.800000000000001e+00 3.999765150548835e-04 + 8.900000000000000e+00 3.999807417161546e-04 + 9.000000000000000e+00 3.999837823577945e-04 + 9.100000000000000e+00 3.999825068985103e-04 + 9.199999999999999e+00 3.999752159160609e-04 + 9.300000000000001e+00 3.999619448716887e-04 + 9.400000000000000e+00 3.999446827130353e-04 + 9.500000000000000e+00 3.999270320511719e-04 + 9.600000000000000e+00 3.999123846519174e-04 + 9.699999999999999e+00 3.999031133009585e-04 + 9.800000000000001e+00 3.999002163065477e-04 + 9.900000000000000e+00 3.999030479451310e-04 + 1.000000000000000e+01 3.999097057656700e-04 diff --git a/figs/correlation.fig/simpleq.dat b/figs/correlation.fig/simpleq.dat new file mode 100644 index 0000000..21074ac --- /dev/null +++ b/figs/correlation.fig/simpleq.dat @@ -0,0 +1,100 @@ + 1.000000000000000e-01 2.576832012945544e-04 + 2.000000000000000e-01 2.591801994046203e-04 + 3.000000000000000e-01 2.614522778615415e-04 + 4.000000000000000e-01 2.643583589076412e-04 + 5.000000000000000e-01 2.677681222327044e-04 + 6.000000000000000e-01 2.715846103272726e-04 + 7.000000000000000e-01 2.757225960482556e-04 + 8.000000000000000e-01 2.800929817186738e-04 + 9.000000000000000e-01 2.846434285234684e-04 + 1.000000000000000e+00 2.893199212355926e-04 + 1.100000000000000e+00 2.940603519412763e-04 + 1.200000000000000e+00 2.988277827132428e-04 + 1.300000000000000e+00 3.035909908501411e-04 + 1.400000000000000e+00 3.083202564269222e-04 + 1.500000000000000e+00 3.129828668185841e-04 + 1.600000000000000e+00 3.175489399236783e-04 + 1.700000000000000e+00 3.220119989331659e-04 + 1.800000000000000e+00 3.263613761073910e-04 + 1.900000000000000e+00 3.305748410709799e-04 + 2.000000000000000e+00 3.346389331284919e-04 + 2.100000000000000e+00 3.385480886964766e-04 + 2.200000000000000e+00 3.423054676538781e-04 + 2.300000000000000e+00 3.459073934204338e-04 + 2.400000000000000e+00 3.493431202324514e-04 + 2.500000000000000e+00 3.526146811685762e-04 + 2.600000000000000e+00 3.557242030505029e-04 + 2.700000000000000e+00 3.586726118760635e-04 + 2.800000000000000e+00 3.614662114908203e-04 + 2.900000000000000e+00 3.641083119885910e-04 + 3.000000000000000e+00 3.666018401752819e-04 + 3.100000000000000e+00 3.689464671595334e-04 + 3.200000000000000e+00 3.711454942522368e-04 + 3.300000000000000e+00 3.732138521709597e-04 + 3.400000000000000e+00 3.751593146275350e-04 + 3.500000000000000e+00 3.769810843056774e-04 + 3.600000000000000e+00 3.786802971063406e-04 + 3.700000000000000e+00 3.802632360936313e-04 + 3.800000000000000e+00 3.817430884369813e-04 + 3.900000000000000e+00 3.831259420996993e-04 + 4.000000000000000e+00 3.844122336094200e-04 + 4.100000000000000e+00 3.856075868184363e-04 + 4.200000000000000e+00 3.867157857196667e-04 + 4.300000000000000e+00 3.877421883466847e-04 + 4.400000000000000e+00 3.886954646057492e-04 + 4.500000000000000e+00 3.895824529571015e-04 + 4.600000000000000e+00 3.904083035655177e-04 + 4.700000000000000e+00 3.911716727892475e-04 + 4.800000000000000e+00 3.918728311874257e-04 + 4.900000000000000e+00 3.925195508586967e-04 + 5.000000000000000e+00 3.931161391707307e-04 + 5.100000000000000e+00 3.936644453388421e-04 + 5.200000000000000e+00 3.941672724782094e-04 + 5.300000000000000e+00 3.946305386436556e-04 + 5.400000000000000e+00 3.950633548128952e-04 + 5.500000000000000e+00 3.954674761268850e-04 + 5.600000000000000e+00 3.958417458083846e-04 + 5.700000000000000e+00 3.961876283635188e-04 + 5.800000000000000e+00 3.965034237988770e-04 + 5.900000000000000e+00 3.967877244043233e-04 + 6.000000000000000e+00 3.970408779450106e-04 + 6.100000000000000e+00 3.972664873238014e-04 + 6.200000000000000e+00 3.974725875030283e-04 + 6.300000000000000e+00 3.976651525934551e-04 + 6.400000000000000e+00 3.978507334953781e-04 + 6.500000000000000e+00 3.980343736709574e-04 + 6.600000000000000e+00 3.982114768954074e-04 + 6.700000000000000e+00 3.983727245071801e-04 + 6.800000000000000e+00 3.985101207624151e-04 + 6.900000000000000e+00 3.986236489855360e-04 + 7.000000000000000e+00 3.987203589794819e-04 + 7.100000000000000e+00 3.988044480206135e-04 + 7.200000000000000e+00 3.988806206101930e-04 + 7.300000000000000e+00 3.989553474459426e-04 + 7.400000000000000e+00 3.990321098828426e-04 + 7.500000000000000e+00 3.991122480017219e-04 + 7.600000000000000e+00 3.991929365132253e-04 + 7.700000000000000e+00 3.992694614994354e-04 + 7.800000000000000e+00 3.993370599658504e-04 + 7.900000000000000e+00 3.993900657750148e-04 + 8.000000000000000e+00 3.994292103851801e-04 + 8.100000000000000e+00 3.994602767744539e-04 + 8.199999999999999e+00 3.994873962055698e-04 + 8.300000000000001e+00 3.995127852700789e-04 + 8.400000000000000e+00 3.995366010201892e-04 + 8.500000000000000e+00 3.995614288866655e-04 + 8.600000000000000e+00 3.995904535022684e-04 + 8.699999999999999e+00 3.996219698319997e-04 + 8.800000000000001e+00 3.996546231782575e-04 + 8.900000000000000e+00 3.996876690865400e-04 + 9.000000000000000e+00 3.997185817907380e-04 + 9.100000000000000e+00 3.997442613215054e-04 + 9.199999999999999e+00 3.997609142987354e-04 + 9.300000000000001e+00 3.997685087804397e-04 + 9.400000000000000e+00 3.997701313603512e-04 + 9.500000000000000e+00 3.997687010738713e-04 + 9.600000000000000e+00 3.997695942562074e-04 + 9.699999999999999e+00 3.997767887499775e-04 + 9.800000000000001e+00 3.997899597617045e-04 + 9.900000000000000e+00 3.998063543389803e-04 + 1.000000000000000e+01 3.998225590499802e-04 diff --git a/figs/energy.fig/Makefile b/figs/energy.fig/Makefile new file mode 100644 index 0000000..e8d3f12 --- /dev/null +++ b/figs/energy.fig/Makefile @@ -0,0 +1,28 @@ +PROJECTNAME=erho_effective erho_fulleq erho_lowhigh + +PDFS=$(addsuffix .pdf, $(PROJECTNAME)) + +all: $(PDFS) + +$(PDFS): + gnuplot $(patsubst %.pdf, %.gnuplot, $@) > $(patsubst %.pdf, %.tex, $@) + pdflatex -file-line-error $(patsubst %.pdf, %.tex, $@) + +install: $(PDFS) + cp $^ $(INSTALLDIR)/ + +$(LIBS): + ln -fs libs/$@ ./ + +clean-libs: + rm -f $(LIBS) + +clean-aux: + rm -f $(addsuffix .aux, $(PROJECTNAME)) + rm -f $(addsuffix .log, $(PROJECTNAME)) + rm -f $(addsuffix .tex, $(PROJECTNAME)) + +clean-tex: + rm -f $(PDFS) + +clean: clean-libs clean-aux clean-tex diff --git a/figs/energy.fig/erho.dat b/figs/energy.fig/erho.dat new file mode 100644 index 0000000..837439f --- /dev/null +++ b/figs/energy.fig/erho.dat @@ -0,0 +1,104 @@ + # v(x)=exp(-|x|), order=100, tolerance=1e-14 + # + # rho energy + # + 1.00000000e-06 7.93445832e-06 + 1.20226443e-06 9.54545091e-06 + 1.44543977e-06 1.14842358e-05 + 1.73780083e-06 1.38177301e-05 + 2.08929613e-06 1.66265809e-05 + 2.51188643e-06 2.00080034e-05 + 3.01995172e-06 2.40792133e-05 + 3.63078055e-06 2.89815779e-05 + 4.36515832e-06 3.48856405e-05 + 5.24807460e-06 4.19972045e-05 + 6.30957344e-06 5.05647046e-05 + 7.58577575e-06 6.08881446e-05 + 9.12010839e-06 7.33299394e-05 + 1.09647820e-05 8.83280780e-05 + 1.31825674e-05 1.06412115e-04 + 1.58489319e-05 1.28222614e-04 + 1.90546072e-05 1.54534805e-04 + 2.29086765e-05 1.86287394e-04 + 2.75422870e-05 2.24617677e-04 + 3.31131121e-05 2.70904370e-04 + 3.98107171e-05 3.26819905e-04 + 4.78630092e-05 3.94394317e-04 + 5.75439937e-05 4.76093386e-04 + 6.91830971e-05 5.74914264e-04 + 8.31763771e-05 6.94502595e-04 + 1.00000000e-04 8.39296065e-04 + 1.20226443e-04 1.01470042e-03 + 1.44543977e-04 1.22730542e-03 + 1.73780083e-04 1.48514988e-03 + 2.08929613e-04 1.79804690e-03 + 2.51188643e-04 2.17798310e-03 + 3.01995172e-04 2.63960836e-03 + 3.63078055e-04 3.20083623e-03 + 4.36515832e-04 3.88357954e-03 + 5.24807460e-04 4.71465030e-03 + 6.30957344e-04 5.72685888e-03 + 7.58577575e-04 6.96035420e-03 + 9.12010839e-04 8.46425370e-03 + 1.09647820e-03 1.02986208e-02 + 1.31825674e-03 1.25368572e-02 + 1.58489319e-03 1.52685874e-02 + 1.90546072e-03 1.86031268e-02 + 2.29086765e-03 2.26736380e-02 + 2.75422870e-03 2.76420977e-02 + 3.31131121e-03 3.37052167e-02 + 3.98107171e-03 4.11014850e-02 + 4.78630092e-03 5.01195435e-02 + 5.75439937e-03 6.11081308e-02 + 6.91830971e-03 7.44879092e-02 + 8.31763771e-03 9.07655413e-02 + 1.00000000e-02 1.10550483e-01 + 1.20226443e-02 1.34575061e-01 + 1.44543977e-02 1.63718539e-01 + 1.73780083e-02 1.99036030e-01 + 2.08929613e-02 2.41793284e-01 + 2.51188643e-02 2.93508614e-01 + 3.01995172e-02 3.56003446e-01 + 3.63078055e-02 4.31463277e-01 + 4.36515832e-02 5.22511186e-01 + 5.24807460e-02 6.32296416e-01 + 6.30957344e-02 7.64601077e-01 + 7.58577575e-02 9.23968602e-01 + 9.12010839e-02 1.11585830e+00 + 1.09647820e-01 1.34683126e+00 + 1.31825674e-01 1.62477387e+00 + 1.58489319e-01 1.95916654e+00 + 1.90546072e-01 2.36140673e+00 + 2.29086765e-01 2.84519724e+00 + 2.75422870e-01 3.42701297e+00 + 3.31131121e-01 4.12666196e+00 + 3.98107171e-01 4.96795983e+00 + 4.78630092e-01 5.97954056e+00 + 5.75439937e-01 7.19583118e+00 + 6.91830971e-01 8.65822353e+00 + 8.31763771e-01 1.04164830e+01 + 1.00000000e+00 1.25304420e+01 + 1.20226443e+00 1.50720365e+01 + 1.44543977e+00 1.81277534e+01 + 1.73780083e+00 2.18015742e+01 + 2.08929613e+00 2.62185132e+01 + 2.51188643e+00 3.15288712e+01 + 3.01995172e+00 3.79133506e+01 + 3.63078055e+00 4.55892039e+01 + 4.36515832e+00 5.48176270e+01 + 5.24807460e+00 6.59126465e+01 + 6.30957344e+00 7.92518063e+01 + 7.58577575e+00 9.52890139e+01 + 9.12010839e+00 1.14569987e+02 + 1.09647820e+01 1.37750822e+02 + 1.31825674e+01 1.65620322e+02 + 1.58489319e+01 1.99126836e+02 + 1.90546072e+01 2.39410530e+02 + 2.29086765e+01 2.87842186e+02 + 2.75422870e+01 3.46069846e+02 + 3.31131121e+01 4.16074894e+02 + 3.98107171e+01 5.00239474e+02 + 4.78630092e+01 6.01427558e+02 + 5.75439937e+01 7.23082394e+02 + 6.91830971e+01 8.69343678e+02 + 8.31763771e+01 1.04518842e+03 diff --git a/figs/energy.fig/erho_effective.gnuplot b/figs/energy.fig/erho_effective.gnuplot new file mode 100644 index 0000000..f26252e --- /dev/null +++ b/figs/energy.fig/erho_effective.gnuplot @@ -0,0 +1,38 @@ +set ylabel "$\\displaystyle\\frac{e}{\\rho}$" norotate offset -1,0 +set xlabel "$\\rho$" + +set xtics 1e-6, 100, 100 +set xtics add ("$10^{-6}$" 0.000001, "$10^{-4}$" 0.0001, "$10^{-2}$" 0.01, "$1$" 1.0, "$10^2$" 100) +unset mxtics + +set ytics 8,1 +set mytics 2 + +set xrange [0.000001:100] +set yrange [7.5:13.5] + +# default output canvas size: 12.5cm x 8.75cm +set term lua tikz size 8,6 standalone + +set key off + + +# set linestyle +set style line 1 linetype rgbcolor "#4169E1" linewidth 3 +set style line 2 linetype rgbcolor "#DC143C" linewidth 3 +set style line 3 linetype rgbcolor "#32CD32" linewidth 3 +set style line 4 linetype rgbcolor "#4B0082" linewidth 3 +set style line 5 linetype rgbcolor "#DAA520" linewidth 3 + +set pointsize 1 + +set logscale x + +a=1.25435641059 + +set label at 0.0001,2*pi*a*1.03 "LHY" textcolor "#DAA520" +set label at 0.1,4.1*pi "Hartree" textcolor "#32CD32" + +plot "erho.dat" using 1:($2/$1) with lines linestyle 1 ,\ + 2*pi*a*(1+128/15/sqrt(pi)*sqrt(x*a**3)) linestyle 5 dashtype "." ,\ + 4*pi linestyle 3 dashtype "." diff --git a/figs/energy.fig/erho_fulleq.dat b/figs/energy.fig/erho_fulleq.dat new file mode 100644 index 0000000..351d1b9 --- /dev/null +++ b/figs/energy.fig/erho_fulleq.dat @@ -0,0 +1,100 @@ + 1.000000000000000e-06 7.943042518656766e-06 0 + 1.202264434617413e-06 9.550219866068997e-06 0 + 1.445439770745928e-06 1.148608224097378e-05 1 + 1.737800828749376e-06 1.381647407426905e-05 1 + 2.089296130854041e-06 1.662160583174429e-05 0 + 2.511886431509582e-06 2.000065135391551e-05 0 + 3.019951720402019e-06 2.406839054942614e-05 0 + 3.630780547701017e-06 2.896662133595565e-05 1 + 4.365158322401657e-06 3.486570260772683e-05 1 + 5.248074602497723e-06 4.197087478032896e-05 1 + 6.309573444801930e-06 5.052956242205629e-05 1 + 7.585775750291836e-06 6.084048678334965e-05 1 + 9.120108393559096e-06 7.326450648917462e-05 1 + 1.096478196143185e-05 8.823760552067889e-05 1 + 1.318256738556407e-05 1.062866339996179e-04 1 + 1.584893192461114e-05 1.280484412743087e-04 1 + 1.905460717963248e-05 1.542931396948374e-04 1 + 2.290867652767775e-05 1.859523935537466e-04 1 + 2.754228703338163e-05 2.241538298379074e-04 1 + 3.311311214825908e-05 2.702629184292266e-04 1 + 3.981071705534969e-05 3.259339774051344e-04 1 + 4.786300923226380e-05 3.931723354869787e-04 1 + 5.754399373371566e-05 4.744101435386617e-04 1 + 6.918309709189363e-05 5.725988899895806e-04 1 + 8.317637711026709e-05 6.913223659756494e-04 1 + 1.000000000000000e-04 8.349346745163116e-04 1 + 1.202264434617413e-04 1.008728918608530e-03 1 + 1.445439770745928e-04 1.219143475424315e-03 1 + 1.737800828749376e-04 1.474014314386416e-03 1 + 2.089296130854039e-04 1.782883700975914e-03 1 + 2.511886431509580e-04 2.157377911005317e-03 1 + 3.019951720402016e-04 2.611669336976909e-03 1 + 3.630780547701014e-04 3.163041683927843e-03 1 + 4.365158322401661e-04 3.832580921405091e-03 1 + 5.248074602497728e-04 4.646019387820095e-03 1 + 6.309573444801930e-04 5.634766055532872e-03 1 + 7.585775750291835e-04 6.837162598523666e-03 1 + 9.120108393559096e-04 8.300012713479526e-03 1 + 1.096478196143185e-03 1.008044129788173e-02 1 + 1.318256738556408e-03 1.224815077758103e-02 1 + 1.584893192461114e-03 1.488815434797412e-02 1 + 1.905460717963246e-03 1.810408049926405e-02 1 + 2.290867652767772e-03 2.202216044849479e-02 1 + 2.754228703338166e-03 2.679603069994577e-02 1 + 3.311311214825911e-03 3.261250778551406e-02 1 + 3.981071705534973e-03 3.969852233927642e-02 1 + 4.786300923226385e-03 4.832943621884326e-02 1 + 5.754399373371567e-03 5.883901075761193e-02 1 + 6.918309709189363e-03 7.163134806415740e-02 1 + 8.317637711026709e-03 8.719519272244317e-02 1 + 1.000000000000000e-02 1.061210611088454e-01 1 + 1.202264434617413e-02 1.291217632376783e-01 1 + 1.445439770745928e-02 1.570570011615084e-01 1 + 1.737800828749376e-02 1.909628719282619e-01 1 + 2.089296130854041e-02 2.320872756251715e-01 1 + 2.511886431509582e-02 2.819324351385837e-01 1 + 3.019951720402019e-02 3.423059810454135e-01 1 + 3.630780547701010e-02 4.153823512868898e-01 1 + 4.365158322401656e-02 5.037766102499891e-01 1 + 5.248074602497723e-02 6.106332154610473e-01 1 + 6.309573444801930e-02 7.397327648968161e-01 1 + 7.585775750291836e-02 8.956203622431048e-01 1 + 9.120108393559097e-02 1.083759966304908e+00 1 + 1.096478196143185e-01 1.310719973979596e+00 1 + 1.318256738556407e-01 1.584396355238566e+00 1 + 1.584893192461114e-01 1.914280945866737e+00 1 + 1.905460717963248e-01 2.311784046064119e+00 1 + 2.290867652767775e-01 2.790622317703228e+00 1 + 2.754228703338169e-01 3.367285183378040e+00 1 + 3.311311214825908e-01 4.061595587917465e+00 1 + 3.981071705534969e-01 4.897384186423424e+00 1 + 4.786300923226380e-01 5.903299886083452e+00 1 + 5.754399373371566e-01 7.113784322007290e+00 1 + 6.918309709189363e-01 8.570243441380599e+00 1 + 8.317637711026709e-01 1.032245608693097e+01 1 + 1.000000000000000e+00 1.243026753303222e+01 1 + 1.202264434617413e+00 1.496562561057760e+01 1 + 1.445439770745928e+00 1.801502869576436e+01 1 + 1.737800828749376e+00 2.168246883802093e+01 1 + 2.089296130854041e+00 2.609297014699965e+01 1 + 2.511886431509582e+00 3.139684282320142e+01 1 + 3.019951720402019e+00 3.777479758723231e+01 1 + 3.630780547701010e+00 4.544409456028400e+01 1 + 4.365158322401657e+00 5.466593586138625e+01 1 + 5.248074602497723e+00 6.575435350738108e+01 1 + 6.309573444801930e+00 7.908689506783082e+01 1 + 7.585775750291836e+00 9.511747067219883e+01 1 + 9.120108393559097e+00 1.143917984731058e+02 1 + 1.096478196143185e+01 1.375659740282897e+02 1 + 1.318256738556407e+01 1.654287952107530e+02 1 + 1.584893192461114e+01 1.989286012247062e+02 1 + 1.905460717963248e+01 2.392055363411829e+02 0 + 2.290867652767775e+01 2.876303621855438e+02 0 + 2.754228703338169e+01 3.458511552270688e+02 0 + 3.311311214825908e+01 4.158492765508882e+02 0 + 3.981071705534969e+01 5.000067330413117e+02 0 + 4.786300923226381e+01 6.011878314695363e+02 0 + 5.754399373371567e+01 7.228352310513908e+02 0 + 6.918309709189363e+01 8.690894472935229e+02 0 + 8.317637711026708e+01 1.044926046890143e+03 0 diff --git a/figs/energy.fig/erho_fulleq.gnuplot b/figs/energy.fig/erho_fulleq.gnuplot new file mode 100644 index 0000000..2927ead --- /dev/null +++ b/figs/energy.fig/erho_fulleq.gnuplot @@ -0,0 +1,37 @@ +set ylabel "$\\displaystyle\\frac{e}{\\rho}$" norotate offset -1,0 +set xlabel "$\\rho$" + +set xtics 1e-6, 100, 100 +set xtics add ("$10^{-6}$" 0.000001, "$10^{-4}$" 0.0001, "$10^{-2}$" 0.01, "$1$" 1.0, "$10^2$" 100) +unset mxtics + +set ytics 8,1 +set mytics 2 + +set xrange [0.000001:100] +set yrange [7.5:13.5] + +# default output canvas size: 12.5cm x 8.75cm +set term lua tikz size 8,6 standalone + +set key off + + +# set linestyle +set style line 1 linetype rgbcolor "#4169E1" linewidth 3 +set style line 2 linetype rgbcolor "#DC143C" linewidth 3 +set style line 3 linetype rgbcolor "#32CD32" linewidth 3 +set style line 4 linetype rgbcolor "#4B0082" linewidth 3 +set style line 5 linetype rgbcolor "#DAA520" linewidth 3 + +set pointsize 1 + +set logscale x + +a=1.25435641059 + +plot "erho.dat" using 1:($2/$1) with lines linestyle 1 ,\ + "erho_fulleq.dat" using 1:($2/$1) with lines linestyle 4 ,\ + 2*pi*a*(1+128/15/sqrt(pi)*sqrt(x*a**3)) linestyle 5 dashtype "." ,\ + 4*pi linestyle 3 dashtype "." ,\ + "holzmann_2019-12-25.dat" using 1:($2/$1) with points linestyle 2 diff --git a/figs/energy.fig/erho_lowhigh.gnuplot b/figs/energy.fig/erho_lowhigh.gnuplot new file mode 100644 index 0000000..9d2a590 --- /dev/null +++ b/figs/energy.fig/erho_lowhigh.gnuplot @@ -0,0 +1,37 @@ +set ylabel "$\\displaystyle\\frac{e}{\\rho}$" norotate offset -1,0 +set xlabel "$\\rho$" + +set xtics 1e-6, 100, 100 +set xtics add ("$10^{-6}$" 0.000001, "$10^{-4}$" 0.0001, "$10^{-2}$" 0.01, "$1$" 1.0, "$10^2$" 100) +unset mxtics + +set ytics 8,1 +set mytics 2 + +set xrange [0.000001:100] +set yrange [7.5:13.5] + +# default output canvas size: 12.5cm x 8.75cm +set term lua tikz size 8,6 standalone + +set key off + + +# set linestyle +set style line 1 linetype rgbcolor "#4169E1" linewidth 3 +set style line 2 linetype rgbcolor "#DC143C" linewidth 3 +set style line 3 linetype rgbcolor "#32CD32" linewidth 3 +set style line 4 linetype rgbcolor "#4B0082" linewidth 3 +set style line 5 linetype rgbcolor "#DAA520" linewidth 3 + +set pointsize 1 + +set logscale x + +a=1.25435641059 + +set label at 0.0001,2*pi*a*1.03 "LHY" textcolor "#DAA520" +set label at 0.1,4.1*pi "Hartree" textcolor "#32CD32" + +plot 2*pi*a*(1+128/15/sqrt(pi)*sqrt(x*a**3)) linestyle 5 dashtype "." ,\ + 4*pi linestyle 3 dashtype "." diff --git a/figs/energy.fig/holzmann_2019-12-25.dat b/figs/energy.fig/holzmann_2019-12-25.dat new file mode 100644 index 0000000..96689cb --- /dev/null +++ b/figs/energy.fig/holzmann_2019-12-25.dat @@ -0,0 +1,11 @@ +## data from M. Holzmann, 2019-09-22 +# rho E0 n0 +1e-6 7.902e-6 +1e-4 8.3441e-4 0.0171 +1e-3 9.13384e-3 0.0481 +1e-2 1.061073e-1 0.0871 +1e-1 1.19182e+0 0.0851 +1e-0 1.24302e+1 0.0587 +1e+1 1.25442e+2 0.0319 +5e+1 6.28032e+2 0.021 + diff --git a/figs/energy16.fig/2020-05-22+16.dat b/figs/energy16.fig/2020-05-22+16.dat new file mode 100644 index 0000000..c28285d --- /dev/null +++ b/figs/energy16.fig/2020-05-22+16.dat @@ -0,0 +1,7 @@ +# rho E0 n0 + 1e-3 0.05195 0.61 + 5e-3 0.449 0.40 + 1e-2 1.1335 0.34 + 2e-2 2.748 0.33 + 1e-1 17.797 0.43 + 1e+0 197.35 0.70 diff --git a/figs/energy16.fig/Makefile b/figs/energy16.fig/Makefile new file mode 100644 index 0000000..964a8c5 --- /dev/null +++ b/figs/energy16.fig/Makefile @@ -0,0 +1,35 @@ +PROJECTNAME=energy16 + +SIMPLEQ=~/Work/Research/2018+bose_gas/cmp/simpleq + +DATS=simpleq16.dat medeq16.dat +PDFS=$(addsuffix .pdf, $(PROJECTNAME)) +TEXS=$(addsuffix .tikz.tex, $(PROJECTNAME)) + +all: $(PDFS) + +$(PDFS): $(DATS) + gnuplot $(patsubst %.pdf, %.gnuplot, $@) > $(patsubst %.pdf, %.tikz.tex, $@) + pdflatex -jobname $(basename $@) -file-line-error $(patsubst %.pdf, %.tikz.tex, $@) + +simpleq16.dat: + julia $(SIMPLEQ)/main.jl -p "tolerance=1e-11;order=100;maxiter=100;v_a=16" -M simpleq energy_rho > $@ +medeq16.dat: + julia -p 8 $(SIMPLEQ)/main.jl -p "aK=1;aL1=1;aL2=1;aL3=0;bK=1;bL1=1;bL2=1;gK=1;gL2=1;gL3=0;N=12;P=8;v_a=16" -M medeq energy_rho > $@ + + +install: $(PDFS) + cp $^ $(INSTALLDIR)/ + +clean-aux: + rm -f $(addsuffix .tikz.tex, $(PROJECTNAME)) + rm -f $(addsuffix .aux, $(PROJECTNAME)) + rm -f $(addsuffix .log, $(PROJECTNAME)) + +clean-dat: + rm -f $(DATS) + +clean-tex: + rm -f $(PDFS) + +clean: clean-aux clean-tex diff --git a/figs/energy16.fig/energy16.gnuplot b/figs/energy16.fig/energy16.gnuplot new file mode 100644 index 0000000..dc996c5 --- /dev/null +++ b/figs/energy16.fig/energy16.gnuplot @@ -0,0 +1,34 @@ +set ylabel "$\\displaystyle\\frac{e}{\\rho}$" norotate offset -1,0 +set xlabel "$\\rho$" + +set xtics 1e-6, 100, 100 +set xtics add ("$10^{-6}$" 0.000001, "$10^{-4}$" 0.0001, "$10^{-2}$" 0.01, "$1$" 1.0, "$10^2$" 100) +unset mxtics + +#set ytics 0.6, 0.1 +set mytics 2 + +set xrange [0.000001:100] + +# default output canvas size: 12.5cm x 8.75cm +set term lua tikz size 8,6 standalone + +set key off + + +# set linestyle +set style line 1 linetype rgbcolor "#4169E1" linewidth 2 +set style line 2 linetype rgbcolor "#DC143C" linewidth 2 +set style line 3 linetype rgbcolor "#32CD32" linewidth 2 +set style line 4 linetype rgbcolor "#4B0082" linewidth 2 +set style line 5 linetype rgbcolor "#DAA520" linewidth 2 + +set pointsize 1 + +set logscale x + +plot \ + "simpleq16.dat" using 1:($3<1e-5 ? ($2/$1) : 1/0) with lines ls 1 ,\ + "medeq16.dat" using 1:($3<1e-5 ? ($2/$1) : 1/0) with lines ls 4 ,\ + "2020-05-22+16.dat" using 1:($2/$1) with points ls 2 + diff --git a/figs/energy16.fig/medeq16.dat b/figs/energy16.fig/medeq16.dat new file mode 100644 index 0000000..b670ef5 --- /dev/null +++ b/figs/energy16.fig/medeq16.dat @@ -0,0 +1,100 @@ + 1.000000000000000e-06 2.554068644475559e-05 1.397500557749096e-14 + 1.202264434617413e-06 3.080347047701689e-05 1.638378932477403e-14 + 1.445439770745928e-06 3.716098929256974e-05 2.068669778175053e-14 + 1.737800828749376e-06 4.484410409494435e-05 9.413457244472722e-15 + 2.089296130854041e-06 5.413333030032436e-05 2.147730951610164e-14 + 2.511886431509582e-06 6.536981275090576e-05 5.079185685242370e-14 + 3.019951720402019e-06 7.896881202049198e-05 4.890613903920614e-14 + 3.630780547701017e-06 9.543630295201912e-05 3.239855276438957e-14 + 4.365158322401657e-06 1.153894399986016e-04 6.114316674680731e-14 + 5.248074602497723e-06 1.395818383298505e-04 2.406593746021106e-14 + 6.309573444801930e-06 1.689348660488532e-04 5.235506117228457e-14 + 7.585775750291836e-06 2.045764562134621e-04 2.215175470514193e-14 + 9.120108393559096e-06 2.478893480268333e-04 1.629554526288622e-14 + 1.096478196143185e-05 3.005711817436355e-04 4.459384031826223e-15 + 1.318256738556407e-05 3.647095370618017e-04 7.959533672028488e-14 + 1.584893192461114e-05 4.428758660698791e-04 1.856627535995191e-14 + 1.905460717963248e-05 5.382433893959377e-04 1.690345207320381e-14 + 2.290867652767775e-05 6.547354778876219e-04 8.321647126329251e-15 + 2.754228703338163e-05 7.972129388261395e-04 1.217647346338456e-14 + 3.311311214825908e-05 9.717111086044649e-04 8.937613010647126e-15 + 3.981071705534969e-05 1.185740914212854e-03 1.053125999750607e-14 + 4.786300923226380e-05 1.448672358365188e-03 2.683673474487137e-14 + 5.754399373371566e-05 1.772224545950021e-03 1.137789331779729e-14 + 6.918309709189363e-05 2.171093851718767e-03 2.458170758492918e-14 + 8.317637711026709e-05 2.663761724994058e-03 3.251467169362365e-14 + 1.000000000000000e-04 3.273536719674966e-03 1.604790218801148e-14 + 1.202264434617413e-04 4.029902651256568e-03 1.208705427621754e-14 + 1.445439770745928e-04 4.970267643766135e-03 9.949882875987985e-15 + 1.737800828749376e-04 6.142238932807714e-03 9.175378143390922e-15 + 2.089296130854039e-04 7.606587760365588e-03 2.214874949422118e-14 + 2.511886431509580e-04 9.441120133796742e-03 8.309109937206045e-14 + 3.019951720402016e-04 1.174573571624843e-02 2.396263674278965e-13 + 3.630780547701014e-04 1.464904218848119e-02 7.677444642221211e-13 + 4.365158322401661e-04 1.831699983067385e-02 2.316125381119014e-12 + 5.248074602497728e-04 2.296420442581745e-02 6.550410386633497e-12 + 6.309573444801930e-04 2.886857870385148e-02 1.868355419086652e-15 + 7.585775750291835e-04 3.639043452656086e-02 2.423897582964416e-15 + 9.120108393559096e-04 4.599708799630445e-02 6.448772157311636e-15 + 1.096478196143185e-03 5.829445137902425e-02 5.163654832442502e-15 + 1.318256738556408e-03 7.406727707716829e-02 9.938237834818324e-15 + 1.584893192461114e-03 9.432997118636552e-02 1.561488943107384e-14 + 1.905460717963246e-03 1.203901027592137e-01 1.480152300205420e-14 + 2.290867652767772e-03 1.539268821910583e-01 2.186407540038661e-14 + 2.754228703338166e-03 1.970869439988538e-01 8.674104990124533e-15 + 3.311311214825911e-03 2.525997349448120e-01 3.829444151273552e-14 + 3.981071705534973e-03 3.239146932311696e-01 1.538039001478893e-14 + 4.786300923226385e-03 4.153622628125486e-01 6.821521006308993e-15 + 5.754399373371567e-03 5.323407233021807e-01 2.800593843254869e-15 + 6.918309709189363e-03 6.815309908475140e-01 4.113794640712594e-15 + 8.317637711026709e-03 8.711421688804751e-01 3.133249774936964e-15 + 1.000000000000000e-02 1.111191929557474e+00 4.457911153045777e-15 + 1.202264434617413e-02 1.413827980332061e+00 6.422681440463997e-15 + 1.445439770745928e-02 1.793700009721424e+00 4.370320881832804e-15 + 1.737800828749376e-02 2.268395559532093e+00 3.253121018609543e-15 + 2.089296130854041e-02 2.858958045662395e+00 3.370516638785516e-15 + 2.511886431509582e-02 3.590510376831549e+00 5.365630506272522e-15 + 3.019951720402019e-02 4.493013035241590e+00 3.715496417945286e-15 + 3.630780547701010e-02 5.602190865095816e+00 3.496056342568733e-15 + 4.365158322401656e-02 6.960668080208332e+00 1.506892881848512e-15 + 5.248074602497723e-02 8.619356424678557e+00 1.746468145174605e-15 + 6.309573444801930e-02 1.063914751302661e+01 1.661095648248939e-15 + 7.585775750291836e-02 1.309296777581257e+01 2.123676640240517e-15 + 9.120108393559097e-02 1.606826398387775e+01 4.350828619157266e-12 + 1.096478196143185e-01 1.966999980308274e+01 1.196771616697282e-12 + 1.318256738556407e-01 2.402425994094341e+01 3.018342123439392e-13 + 1.584893192461114e-01 2.928257887420183e+01 6.985575344309217e-14 + 1.905460717963248e-01 3.562713643064021e+01 2.023128031947868e-14 + 2.290867652767775e-01 4.327699319810110e+01 7.636489971812963e-14 + 2.754228703338169e-01 5.249557556550972e+01 3.450372875056261e-14 + 3.311311214825908e-01 6.359966407666808e+01 1.050825065430309e-14 + 3.981071705534969e-01 7.697019095971760e+01 1.016139346429647e-14 + 4.786300923226380e-01 9.306521489952080e+01 6.058578163892193e-15 + 5.754399373371566e-01 1.124355155175409e+02 3.382728454855092e-15 + 6.918309709189363e-01 1.357433392373146e+02 2.709608246018857e-15 + 8.317637711026709e-01 1.637849353682065e+02 1.878556843412201e-15 + 1.000000000000000e+00 1.975176500593190e+02 1.223063234010325e-15 + 1.202264434617413e+00 2.380925007021913e+02 1.993050434866885e-15 + 1.445439770745928e+00 2.868933397000474e+02 3.043149243838522e-15 + 1.737800828749376e+00 3.455839406222755e+02 2.998782498464775e-15 + 2.089296130854041e+00 4.161646092494719e+02 2.611515016643540e-15 + 2.511886431509582e+00 5.010402460546079e+02 1.646713357943204e-15 + 3.019951720402019e+00 6.031021762840705e+02 2.632388911909736e-15 + 3.630780547701010e+00 7.258265322526299e+02 3.725804901562014e-15 + 4.365158322401657e+00 8.733925356842659e+02 3.933266540669057e-15 + 5.248074602497723e+00 1.050824705077102e+03 3.818189041855540e-15 + 6.309573444801930e+00 1.264163827184745e+03 3.557698468765215e-15 + 7.585775750291836e+00 1.520672510488972e+03 3.200776178020959e-15 + 9.120108393559097e+00 1.829082315295077e+03 2.768978034961086e-15 + 1.096478196143185e+01 2.199890869853228e+03 2.784644153002823e-15 + 1.318256738556407e+01 2.645719082839147e+03 2.561848808069188e-15 + 1.584893192461114e+01 3.181740607495819e+03 3.342103900653102e-15 + 1.905460717963248e+01 3.826198171331817e+03 3.061126283356729e-15 + 2.290867652767775e+01 4.601024341151761e+03 3.012741771786859e-15 + 2.754228703338169e+01 5.532587846940884e+03 3.540875919778527e-15 + 3.311311214825908e+01 6.652590860664027e+03 3.384458193003484e-15 + 3.981071705534969e+01 7.999147762757379e+03 3.313778691885289e-15 + 4.786300923226381e+01 9.618082104785730e+03 4.116261947533335e-15 + 5.754399373371567e+01 1.156448590156131e+04 3.799276533234762e-15 + 6.918309709189363e+01 1.390459431262479e+04 4.136504038005983e-15 + 8.317637711026708e+01 1.671803950512843e+04 4.852059592739983e-15 diff --git a/figs/energy16.fig/medeq4.dat b/figs/energy16.fig/medeq4.dat new file mode 100644 index 0000000..ffaa5eb --- /dev/null +++ b/figs/energy16.fig/medeq4.dat @@ -0,0 +1,100 @@ + 1.000000000000000e-06 1.627424240030753e-05 9.007835427072701e-12 + 1.202264434617413e-06 1.959943745344100e-05 1.671351637616213e-12 + 1.445439770745928e-06 2.360770720466136e-05 4.387069667322242e-12 + 1.737800828749376e-06 2.844048299127554e-05 3.371679780930921e-14 + 2.089296130854041e-06 3.426881796805706e-05 3.499980814968632e-14 + 2.511886431509582e-06 4.129969708194542e-05 7.218204810983588e-15 + 3.019951720402019e-06 4.978372244366768e-05 8.253688999006713e-15 + 3.630780547701017e-06 6.002448197752047e-05 1.233139066041191e-14 + 4.365158322401657e-06 7.238998092442059e-05 6.922441203895261e-15 + 5.248074602497723e-06 8.732660537970414e-05 2.847530345897527e-14 + 6.309573444801930e-06 1.053761990232000e-04 1.507914112842546e-14 + 7.585775750291836e-06 1.271969739658683e-04 7.749491820001025e-12 + 9.120108393559096e-06 1.535891510001255e-04 3.420455433776017e-12 + 1.096478196143185e-05 1.855264423691571e-04 1.192194927362554e-12 + 1.318256738556407e-05 2.241947631453101e-04 4.676335170879237e-13 + 1.584893192461114e-05 2.710399006406206e-04 6.056007240771575e-13 + 1.905460717963248e-05 3.278263043888957e-04 1.171078624784598e-12 + 2.290867652767775e-05 3.967097065782638e-04 1.794813944769702e-12 + 2.754228703338163e-05 4.803269751408225e-04 2.251917988222614e-12 + 3.311311214825908e-05 5.819074780693888e-04 2.415044085011130e-12 + 3.981071705534969e-05 7.054113479098947e-04 2.308357874960898e-12 + 4.786300923226380e-05 8.557014437731245e-04 2.127244943807619e-12 + 5.754399373371566e-05 1.038757596087362e-03 2.664244293528506e-12 + 6.918309709189363e-05 1.261943989591930e-03 5.395675890368970e-12 + 8.317637711026709e-05 1.534343421852030e-03 1.539374470698132e-14 + 1.000000000000000e-04 1.867175828671813e-03 1.597016530606916e-14 + 1.202264434617413e-04 2.274323092769955e-03 9.714938856241721e-15 + 1.445439770745928e-04 2.772987990027553e-03 1.269495496436939e-14 + 1.737800828749376e-04 3.384522468901346e-03 2.589592076658382e-14 + 2.089296130854039e-04 4.135469643992774e-03 6.402903743410984e-15 + 2.511886431509580e-04 5.058875303255412e-03 1.841091338474481e-14 + 3.019951720402016e-04 6.195938807158513e-03 7.894587486448972e-15 + 3.630780547701014e-04 7.598090441518367e-03 1.332917339601970e-14 + 4.365158322401661e-04 9.329602997799082e-03 1.274924502157576e-14 + 5.248074602497728e-04 1.147086994497691e-02 9.168558395794375e-15 + 6.309573444801930e-04 1.412251122865532e-02 1.086521688806541e-14 + 7.585775750291835e-04 1.741050046104560e-02 1.304008186909502e-14 + 9.120108393559096e-04 2.149254371940543e-02 7.577632567649028e-15 + 1.096478196143185e-03 2.656597967339182e-02 1.354559449592211e-14 + 1.318256738556408e-03 3.287751233543990e-02 7.370000623409608e-15 + 1.584893192461114e-03 4.073513029684535e-02 6.378750030476871e-15 + 1.905460717963246e-03 5.052260913087291e-02 3.698643037579531e-15 + 2.290867652767772e-03 6.271703699425088e-02 7.017835372202290e-15 + 2.754228703338166e-03 7.790984958491401e-02 4.177399466817889e-15 + 3.311311214825911e-03 9.683191476906949e-02 4.287585917545763e-15 + 3.981071705534973e-03 1.203832784973463e-01 7.313787959453048e-15 + 4.786300923226385e-03 1.496682853487706e-01 4.816305546585549e-15 + 5.754399373371567e-03 1.860369358262502e-01 5.546470847879690e-15 + 6.918309709189363e-03 2.311335564559619e-01 3.744556832996796e-15 + 8.317637711026709e-03 2.869541543528080e-01 2.252401346512346e-15 + 1.000000000000000e-02 3.559142170277424e-01 5.621737354494605e-15 + 1.202264434617413e-02 4.409292057883803e-01 4.681839933973202e-15 + 1.445439770745928e-02 5.455105753566664e-01 1.546699707615607e-15 + 1.737800828749376e-02 6.738808279465798e-01 5.874098214125763e-15 + 2.089296130854041e-02 8.311118738264928e-01 5.149027037818788e-15 + 2.511886431509582e-02 1.023291828096413e+00 5.696277355342162e-15 + 3.019951720402019e-02 1.257726345338623e+00 4.432800724331350e-15 + 3.630780547701010e-02 1.543181717371867e+00 4.005749823244296e-15 + 4.365158322401656e-02 1.890178288183423e+00 2.239736734668354e-15 + 5.248074602497723e-02 2.311344340711753e+00 2.245605809580106e-15 + 6.309573444801930e-02 2.821842558067281e+00 1.645883940583995e-15 + 7.585775750291836e-02 3.439883539443798e+00 3.304227952844349e-15 + 9.120108393559097e-02 4.187343748169733e+00 2.467510197024562e-15 + 1.096478196143185e-01 5.090508786399249e+00 8.482875770705394e-16 + 1.318256738556407e-01 6.180967143946303e+00 2.119874330239340e-15 + 1.584893192461114e-01 7.496684696001272e+00 5.665629796095111e-15 + 1.905460717963248e-01 9.083296393338594e+00 3.907507356933792e-15 + 2.290867652767775e-01 1.099565900285665e+01 4.575037511918055e-15 + 2.754228703338169e-01 1.329971766360591e+01 2.883241152474376e-15 + 3.311311214825908e-01 1.607474972409581e+01 5.147537259204491e-15 + 3.981071705534969e-01 1.941606218332745e+01 2.362757842000826e-15 + 4.786300923226380e-01 2.343823450777575e+01 2.050801375388800e-15 + 5.754399373371566e-01 2.827901716541566e+01 2.392795176522673e-15 + 6.918309709189363e-01 3.410401853827422e+01 2.477737430323620e-15 + 8.317637711026709e-01 4.111233970707951e+01 1.553492938204756e-15 + 1.000000000000000e+00 4.954334885903700e+01 1.366441830606871e-15 + 1.202264434617413e+00 5.968482585078009e+01 3.691739631003516e-15 + 1.445439770745928e+00 7.188275408317689e+01 1.597349795336168e-15 + 1.737800828749376e+00 8.655309289979826e+01 1.415290168167683e-15 + 2.089296130854041e+00 1.041959311141499e+02 2.406080183335882e-15 + 2.511886431509582e+00 1.254125032956741e+02 1.889598048843391e-15 + 3.019951720402019e+00 1.509256478582651e+02 1.758967767503912e-15 + 3.630780547701010e+00 1.816044031124912e+02 2.441779747875737e-15 + 4.365158322401657e+00 2.184935782492489e+02 1.921466737271511e-15 + 5.248074602497723e+00 2.628493055096012e+02 1.991374998750417e-15 + 6.309573444801930e+00 3.161817833236910e+02 2.268079781904007e-15 + 7.585775750291836e+00 3.803066648964960e+02 2.072951425856464e-15 + 9.120108393559097e+00 4.574068409063933e+02 3.147983516813687e-15 + 1.096478196143185e+01 5.501067186746773e+02 4.623214476237747e-15 + 1.318256738556407e+01 6.615615253958870e+02 5.260378128143234e-15 + 1.584893192461114e+01 7.955646742599187e+02 3.587644312061816e-15 + 1.905460717963248e+01 9.566768469440615e+02 3.861851068125671e-15 + 2.290867652767775e+01 1.150381184921633e+03 8.683260209845366e-12 + 2.754228703338169e+01 1.383269870469561e+03 6.646886653590585e-12 + 3.311311214825908e+01 1.663268446391906e+03 5.054756332839416e-12 + 3.981071705534969e+01 1.999905507656653e+03 3.819832671531302e-12 + 4.786300923226381e+01 2.404636942067389e+03 2.870439604684976e-12 + 5.754399373371567e+01 2.891235753296954e+03 2.143827953509087e-12 + 6.918309709189363e+01 3.476260731260002e+03 1.592331716665501e-12 + 8.317637711026708e+01 4.179619917834820e+03 1.176753139692823e-12 diff --git a/figs/energy16.fig/simpleq16.dat b/figs/energy16.fig/simpleq16.dat new file mode 100644 index 0000000..2bc4ebc --- /dev/null +++ b/figs/energy16.fig/simpleq16.dat @@ -0,0 +1,100 @@ + 1.000000000000000e-06 2.561870893688601e-05 1.001128563083394e-13 + 1.202264434617413e-06 3.091259293978515e-05 1.166419917451259e-13 + 1.445439770745928e-06 3.731355959329154e-05 1.641671221071265e-13 + 1.737800828749376e-06 4.505737456991875e-05 1.318046204468745e-13 + 2.089296130854041e-06 5.443138706366596e-05 6.465028107789633e-14 + 2.511886431509582e-06 6.578626096553237e-05 1.282631449030218e-13 + 3.019951720402019e-06 7.955051787513295e-05 1.828007544593001e-13 + 3.630780547701017e-06 9.624861222940964e-05 8.716203682159707e-14 + 4.365158322401657e-06 1.165234581314891e-04 2.567118665713670e-13 + 5.248074602497723e-06 1.411645865934592e-04 1.909114183041957e-13 + 6.309573444801930e-06 1.711434503024673e-04 4.832552648271436e-13 + 7.585775750291836e-06 2.076578368881950e-04 9.663351363715864e-13 + 9.120108393559096e-06 2.521879364891216e-04 1.532273658747508e-12 + 1.096478196143185e-05 3.065673834006124e-04 2.557647979444216e-12 + 1.318256738556407e-05 3.730736204970770e-04 4.227854029289540e-12 + 1.584893192461114e-05 4.545433088768755e-04 7.262584516187726e-12 + 1.905460717963248e-05 5.545203471463022e-04 1.581480833654367e-13 + 2.290867652767775e-05 6.774465443911747e-04 1.081736075443477e-13 + 2.754228703338163e-05 8.289083412906583e-04 9.892826615755524e-14 + 3.311311214825908e-05 1.015957515485774e-03 1.000146289881600e-13 + 3.981071705534969e-05 1.247529980092392e-03 1.728836305265867e-13 + 4.786300923226380e-05 1.534995189341161e-03 4.874749473980382e-14 + 5.754399373371566e-05 1.892880118596425e-03 7.675596130065904e-14 + 6.918309709189363e-05 2.339827381972870e-03 3.488065867913569e-14 + 8.317637711026709e-05 2.899868231941961e-03 4.074847689970330e-13 + 1.000000000000000e-04 3.604119801371417e-03 8.837677779017190e-14 + 1.202264434617413e-04 4.493054270740071e-03 1.141726470797567e-13 + 1.445439770745928e-04 5.619538266005139e-03 2.367473489955975e-13 + 1.737800828749376e-04 7.052906297447757e-03 9.521496171923438e-14 + 2.089296130854039e-04 8.884414234736409e-03 2.852315979737347e-13 + 2.511886431509580e-04 1.123451720728305e-02 8.027724386154160e-13 + 3.019951720402016e-04 1.426252568892359e-02 1.058787720274644e-13 + 3.630780547701014e-04 1.817930030054440e-02 1.539765713616111e-13 + 4.365158322401661e-04 2.326372414304216e-02 6.040763684158786e-14 + 5.248074602497728e-04 2.988370094092252e-02 8.733135409851157e-13 + 6.309573444801930e-04 3.852231787659451e-02 3.684231597353426e-14 + 7.585775750291835e-04 4.980954220480775e-02 1.647277211646135e-13 + 9.120108393559096e-04 6.455939366413402e-02 2.032722085217800e-13 + 1.096478196143185e-03 8.381204323959957e-02 2.838665806426261e-13 + 1.318256738556408e-03 1.088799379503562e-01 9.260531251903156e-14 + 1.584893192461114e-03 1.413971109762464e-01 6.326291032696907e-14 + 1.905460717963246e-03 1.833715077279505e-01 4.633167297377464e-14 + 2.290867652767772e-03 2.372414190490330e-01 1.658072106704119e-13 + 2.754228703338166e-03 3.059386735352517e-01 1.122897591215287e-13 + 3.311311214825911e-03 3.929626900151192e-01 5.692534927020893e-14 + 3.981071705534973e-03 5.024705334214158e-01 5.657153215877982e-14 + 4.786300923226385e-03 6.393887201581897e-01 9.048650405950154e-14 + 5.754399373371567e-03 8.095528700733856e-01 8.607211052131642e-14 + 6.918309709189363e-03 1.019881656036838e+00 1.595875861854557e-13 + 8.317637711026709e-03 1.278592057372205e+00 1.531883433865405e-13 + 1.000000000000000e-02 1.595463790063535e+00 1.407687268306203e-13 + 1.202264434617413e-02 1.982162017634850e+00 6.748557036716531e-14 + 1.445439770745928e-02 2.452629079913038e+00 6.869589239332279e-14 + 1.737800828749376e-02 3.023558061280176e+00 2.437733106567872e-14 + 2.089296130854041e-02 3.714963622955766e+00 1.064344070654439e-13 + 2.511886431509582e-02 4.550868731266756e+00 9.238668475341097e-15 + 3.019951720402019e-02 5.560129829441683e+00 1.020232616470553e-13 + 3.630780547701010e-02 6.777427747431538e+00 5.074743271963278e-14 + 4.365158322401656e-02 8.244457369990231e+00 1.028088279693266e-13 + 5.248074602497723e-02 1.001135596595401e+01 4.624881313939891e-12 + 6.309573444801930e-02 1.213841833855832e+01 1.852661419803306e-12 + 7.585775750291836e-02 1.469815685330527e+01 7.258972029472108e-13 + 9.120108393559097e-02 1.777777626084305e+01 2.891235204847115e-13 + 1.096478196143185e-01 2.148214745283791e+01 1.401726396203463e-13 + 1.318256738556407e-01 2.593738134860670e+01 5.743044386657572e-14 + 1.584893192461114e-01 3.129512458927278e+01 3.606825259935793e-14 + 1.905460717963248e-01 3.773772331343910e+01 1.892505250290776e-13 + 2.290867652767775e-01 4.548443084468321e+01 2.083768656818398e-14 + 2.754228703338169e-01 5.479887064684699e+01 4.075384057289656e-14 + 3.311311214825908e-01 6.599800861142323e+01 2.600299498350245e-14 + 3.981071705534969e-01 7.946294008944992e+01 1.413936909490012e-14 + 4.786300923226380e-01 9.565185881847951e+01 2.463069066912327e-14 + 5.754399373371566e-01 1.151156491263624e+02 8.558078927823013e-15 + 6.918309709189363e-01 1.385166320450518e+02 6.226596123328673e-14 + 8.317637711026709e-01 1.666511032772171e+02 7.123438110270025e-14 + 1.000000000000000e+00 2.004764299779804e+02 8.748669232588792e-14 + 1.202264434617413e+00 2.411436284338656e+02 6.564186088832174e-14 + 1.445439770745928e+00 2.900365312193817e+02 6.210366709884162e-14 + 1.737800828749376e+00 3.488188766344850e+02 5.060423878974540e-15 + 2.089296130854041e+00 4.194909228031429e+02 4.504831194472809e-15 + 2.511886431509582e+00 5.044575129173239e+02 5.178886653500039e-14 + 3.019951720402019e+00 6.066099077704147e+02 5.546884105120939e-14 + 3.630780547701010e+00 7.294241701985674e+02 4.099740962135071e-15 + 4.365158322401657e+00 8.770794492790153e+02 5.461022666145280e-14 + 5.248074602497723e+00 1.054600189286513e+03 2.773788643238822e-15 + 6.309573444801930e+00 1.268027102525078e+03 5.605531588703275e-14 + 7.585775750291836e+00 1.524622723934775e+03 1.548482473273102e-15 + 9.120108393559097e+00 1.833118542128751e+03 5.542616472939282e-15 + 1.096478196143185e+01 2.204012116286832e+03 5.416841008960720e-15 + 1.318256738556407e+01 2.649924289260801e+03 1.051948766670258e-13 + 1.584893192461114e+01 3.186028652212697e+03 3.020416191025176e-15 + 1.905460717963248e+01 3.830567874697323e+03 1.054514720940622e-13 + 2.290867652767775e+01 4.605474469975405e+03 1.054818735884649e-13 + 2.754228703338169e+01 5.537117119087241e+03 1.068333765010024e-13 + 3.311311214825908e+01 6.657197949756010e+03 5.311424682002792e-14 + 3.981071705534969e+01 8.003831302911970e+03 8.573716686070972e-15 + 4.786300923226381e+01 9.622840695326604e+03 5.255692893479381e-14 + 5.754399373371567e+01 1.156931811166827e+04 5.412915990778012e-16 + 6.918309709189363e+01 1.390949868589013e+04 5.629397185139020e-14 + 8.317637711026708e+01 1.672301456399418e+04 8.852629630882290e-12 diff --git a/figs/energy_relative_error.fig/2020-10-15+energy.dat b/figs/energy_relative_error.fig/2020-10-15+energy.dat new file mode 100644 index 0000000..8a6566c --- /dev/null +++ b/figs/energy_relative_error.fig/2020-10-15+energy.dat @@ -0,0 +1,13 @@ +1e-6 7.90e-6 7.91e-6 +1e-4 8.344e-4 8.345e-4 0.017 0.018 8.355e-4 8.356e-4 0.015 0.016 +1e-3 9.1338e-3 9.1339e-3 0.048 0.049 9.142e-3 9.143e-3 0.045 0.046 +5e-3 5.0613e-2 5.0614e-2 0.080 0.081 5.0648e-2 5.0649e-2 0.070 0.071 +1e-2 1.06107e-1 1.06108e-1 0.087 0.088 1.0627e-1 1.0628e-1 0.073 0.074 +1.5e-2 1.6337e-1 1.6338e-1 0.090 0.091 1.6342e-1 1.6343e-1 0.088 0.089 +2e-2 2.2162e-1 2.2163e-1 0.093 0.094 2.2169e-1 2.2170e-1 0.090 0.091 +5e-2 5.8044e-1 5.8045e-1 0.091 0.092 5.8046e-1 5.8047e-1 0.091 0.092 +1e-1 1.1918e+0 1.1919e+0 0.088 0.089 1.1919e+0 1.1920e+0 0.085 0.086 +1e-0 1.2430e+1 1.2431e+1 0.058 0.059 1.245e+1 1.246e+1 0.05 0.06 +1e+1 1.2544e+2 1.2545e+2 0.031 0.032 +5e+1 6.2803e+2 6.2804e+2 0.02 0.03 + diff --git a/figs/energy_relative_error.fig/Makefile b/figs/energy_relative_error.fig/Makefile new file mode 100644 index 0000000..93e03f5 --- /dev/null +++ b/figs/energy_relative_error.fig/Makefile @@ -0,0 +1,51 @@ +PROJECTNAME=cmp_energy1 + +SIMPLEQ=~/Work/Research/2018+bose_gas/cmp/simpleq + +DATS= simpleq1.dat mueq1.dat bigeq1.dat +PDFS=$(addsuffix .pdf, $(PROJECTNAME)) +TEXS=$(addsuffix .tikz.tex, $(PROJECTNAME)) + +all: $(PDFS) + +$(PDFS): $(DATS) + gnuplot $(patsubst %.pdf, %.gnuplot, $@) > $(patsubst %.pdf, %.tikz.tex, $@) + pdflatex -jobname $(basename $@) -file-line-error $(patsubst %.pdf, %.tikz.tex, $@) + +bigeq1.dat: + for rho in 1e-6 1e-4 1e-3 5e-3 1e-2 1.5e-2 2e-2 5e-2 1e-1 1.0 1e1 5e1; do \ + echo -n $$rho " " >> $@-tmp ;\ + julia $(SIMPLEQ)/main.jl -p "N=12;P=8;J=10;maxiter=100;v_a=1;rho=$$rho;minlrho=-6;nlrho=50;eq=bigeq" -M medeq energy >> $@-tmp ;\ + done + paste $@-tmp 2020-10-15+energy.dat > $@ + rm $@-tmp +simpleq1.dat: + for lrho in -6 -4 -3 -2.3010299956639813 -2 -1.8239087409443189 -1.6989700043360187 -1.3010299956639813 -1 0 1 1.6989700043360187; do\ + julia $(SIMPLEQ)/main.jl -p "tolerance=1e-11;order=100;maxiter=100;v_a=1;minlrho=$$lrho;nlrho=1" -M simpleq energy_rho >> $@-tmp ;\ + done + paste $@-tmp 2020-10-15+energy.dat > $@ + rm $@-tmp +mueq1.dat: + for lrho in -6 -4 -3 -2.3010299956639813 -2 -1.8239087409443189 -1.6989700043360187 -1.3010299956639813 -1 0 1 1.6989700043360187; do\ + julia $(SIMPLEQ)/main.jl -p "tolerance=1e-11;order=100;maxiter=100;v_a=1;minlrho=$$lrho;nlrho=1" -M mueq energy_rho >> $@-tmp ;\ + done + paste $@-tmp 2020-10-15+energy.dat > $@ + rm $@-tmp + + +install: $(PDFS) + cp $^ $(INSTALLDIR)/ + +clean-aux: + rm -f $(addsuffix .tikz.tex, $(PROJECTNAME)) + rm -f $(addsuffix .aux, $(PROJECTNAME)) + rm -f $(addsuffix .log, $(PROJECTNAME)) + +clean-dat: + rm -f $(DATS) + rm -f Abar.dat + +clean-tex: + rm -f $(PDFS) + +clean: clean-aux clean-tex diff --git a/figs/energy_relative_error.fig/bigeq1.dat b/figs/energy_relative_error.fig/bigeq1.dat new file mode 100644 index 0000000..ccab7cc --- /dev/null +++ b/figs/energy_relative_error.fig/bigeq1.dat @@ -0,0 +1,13 @@ +1e-6 7.933345669932108e-06 1.702110068633820e-12 1e-6 7.90e-6 7.91e-6 +1e-4 8.349333596092273e-04 1.112527655484278e-13 1e-4 8.344e-4 8.345e-4 0.017 0.018 8.355e-4 8.356e-4 0.015 0.016 +1e-3 9.146316224214820e-03 2.077370190827862e-12 1e-3 9.1338e-3 9.1339e-3 0.048 0.049 9.142e-3 9.143e-3 0.045 0.046 +5e-3 5.063551121414726e-02 8.398015568932763e-12 5e-3 5.0613e-2 5.0614e-2 0.080 0.081 5.0648e-2 5.0649e-2 0.070 0.071 +1e-2 1.061166635113156e-01 5.017547334745932e-15 1e-2 1.06107e-1 1.06108e-1 0.087 0.088 1.0627e-1 1.0628e-1 0.073 0.074 +1.5e-2 1.633545333357252e-01 2.212732350724188e-15 1.5e-2 1.6337e-1 1.6338e-1 0.090 0.091 1.6342e-1 1.6343e-1 0.088 0.089 +2e-2 2.216006926766234e-01 4.554921896809050e-15 2e-2 2.2162e-1 2.2163e-1 0.093 0.094 2.2169e-1 2.2170e-1 0.090 0.091 +5e-2 5.805467427278628e-01 3.333614946727334e-15 5e-2 5.8044e-1 5.8045e-1 0.091 0.092 5.8046e-1 5.8047e-1 0.091 0.092 +1e-1 1.191928562890024e+00 2.450444415773880e-15 1e-1 1.1918e+0 1.1919e+0 0.088 0.089 1.1919e+0 1.1920e+0 0.085 0.086 +1.0 1.243051761633436e+01 3.791416729017179e-12 1e-0 1.2430e+1 1.2431e+1 0.058 0.059 1.245e+1 1.246e+1 0.05 0.06 +1e1 1.254476421597336e+02 2.677873467933947e-13 1e+1 1.2544e+2 1.2545e+2 0.031 0.032 +5e1 nan nan 5e+1 6.2803e+2 6.2804e+2 0.02 0.03 + diff --git a/figs/energy_relative_error.fig/cmp_energy1.gnuplot b/figs/energy_relative_error.fig/cmp_energy1.gnuplot new file mode 100644 index 0000000..c2b1afd --- /dev/null +++ b/figs/energy_relative_error.fig/cmp_energy1.gnuplot @@ -0,0 +1,37 @@ +set title "$v(\\mathbf x)=e^{-|\\mathbf x|}$" +set ylabel "energy relative error" offset +3 +set xlabel "$\\rho$" + +set xtics 1e-6, 100, 100 +set xtics add ("$10^{-6}$" 0.000001, "$10^{-4}$" 0.0001, "$10^{-2}$" 0.01, "$1$" 1.0, "$10^2$" 100) +unset mxtics +set xrange [1e-6:1e2] + +set ytics 1e-5, 10, 1e-1 +set ytics add ("$10^{-5}$" 1e-5, "$10^{-4}$" 1e-4, "$10^{-3}$" 1e-3, "$10^{-2}$" 1e-2, "$10^{-1}$" 1e-1) +set yrange [1e-5:1e-1] + +# default output canvas size: 12.5cm x 8.75cm +set term lua tikz size 8,6 standalone + +unset key + + +# set linestyle +set style line 1 linetype rgbcolor "#4169E1" linewidth 2 +set style line 2 linetype rgbcolor "#DC143C" linewidth 2 +set style line 3 linetype rgbcolor "#32CD32" linewidth 2 +set style line 4 linetype rgbcolor "#4B0082" linewidth 2 +set style line 5 linetype rgbcolor "#DAA520" linewidth 2 + +set pointsize 1 + +max(x,y)=(x>y ? x : y) + +set logscale xy + +plot \ + "simpleq1.dat" using 1:($3<1e-5 ? max(abs($2-$5)/$5,abs($2-$6)/$6) : 1/0) every ::1 with lp ls 1 ,\ + "bigeq1.dat" using 1:($3<1e-5 ? max(abs($2-$5)/$5,abs($2-$6)/$6) : 1/0) every ::1 with lp ls 4 ,\ + "2020-10-15+energy.dat" using 1:(max(max(abs($6-$2)/$2,abs($6-$3)/$3),max(abs($7-$2)/$2,abs($7-$3)/$3))) every ::1 with lp ls 2 + diff --git a/figs/energy_relative_error.fig/mueq1.dat b/figs/energy_relative_error.fig/mueq1.dat new file mode 100644 index 0000000..6680ed3 --- /dev/null +++ b/figs/energy_relative_error.fig/mueq1.dat @@ -0,0 +1,13 @@ + 1.000000000000000e-06 7.934286257047026e-06 3.235575358862730e-14 1e-6 7.90e-6 7.91e-6 + 1.000000000000000e-04 8.377752370367244e-04 3.201262322514326e-14 1e-4 8.344e-4 8.345e-4 0.017 0.018 8.355e-4 8.356e-4 0.015 0.016 + 1.000000000000000e-03 9.227516472314044e-03 5.926882960328585e-12 1e-3 9.1338e-3 9.1339e-3 0.048 0.049 9.142e-3 9.143e-3 0.045 0.046 + 4.999999999999999e-03 5.110614907540598e-02 1.069784208067329e-14 5e-3 5.0613e-2 5.0614e-2 0.080 0.081 5.0648e-2 5.0649e-2 0.070 0.071 + 1.000000000000000e-02 1.069255387977072e-01 1.572627931557844e-14 1e-2 1.06107e-1 1.06108e-1 0.087 0.088 1.0627e-1 1.0628e-1 0.073 0.074 + 1.500000000000000e-02 1.643925118385577e-01 2.494863657028789e-14 1.5e-2 1.6337e-1 1.6338e-1 0.090 0.091 1.6342e-1 1.6343e-1 0.088 0.089 + 2.000000000000000e-02 2.228019992986293e-01 1.448625206361080e-14 2e-2 2.2162e-1 2.2163e-1 0.093 0.094 2.2169e-1 2.2170e-1 0.090 0.091 + 5.000000000000000e-02 5.821600306289594e-01 2.429130993459621e-14 5e-2 5.8044e-1 5.8045e-1 0.091 0.092 5.8046e-1 5.8047e-1 0.091 0.092 + 1.000000000000000e-01 1.193625606724180e+00 9.289627651289616e-14 1e-1 1.1918e+0 1.1919e+0 0.088 0.089 1.1919e+0 1.1920e+0 0.085 0.086 + 1.000000000000000e+00 1.243092200540192e+01 3.956891809887943e-12 1e-0 1.2430e+1 1.2431e+1 0.058 0.059 1.245e+1 1.246e+1 0.05 0.06 + 1.000000000000000e+01 1.254464788052155e+02 1.136222076562384e-14 1e+1 1.2544e+2 1.2545e+2 0.031 0.032 + 4.999999999999999e+01 6.280479158155819e+02 4.085062088184932e-15 5e+1 6.2803e+2 6.2804e+2 0.02 0.03 + diff --git a/figs/energy_relative_error.fig/simpleq1.dat b/figs/energy_relative_error.fig/simpleq1.dat new file mode 100644 index 0000000..0a9946c --- /dev/null +++ b/figs/energy_relative_error.fig/simpleq1.dat @@ -0,0 +1,13 @@ + 1.000000000000000e-06 7.934458317554763e-06 1.386474846476237e-14 1e-6 7.90e-6 7.91e-6 + 1.000000000000000e-04 8.392960650211980e-04 4.992772141681879e-14 1e-4 8.344e-4 8.345e-4 0.017 0.018 8.355e-4 8.356e-4 0.015 0.016 + 1.000000000000000e-03 9.335885024335848e-03 9.476600935679491e-12 1e-3 9.1338e-3 9.1339e-3 0.048 0.049 9.142e-3 9.143e-3 0.045 0.046 + 4.999999999999999e-03 5.253256909806337e-02 2.904203225041083e-14 5e-3 5.0613e-2 5.0614e-2 0.080 0.081 5.0648e-2 5.0649e-2 0.070 0.071 + 1.000000000000000e-02 1.105504827786785e-01 3.695967730970493e-14 1e-2 1.06107e-1 1.06108e-1 0.087 0.088 1.0627e-1 1.0628e-1 0.073 0.074 + 1.500000000000000e-02 1.702888708763101e-01 2.044316796426063e-14 1.5e-2 1.6337e-1 1.6338e-1 0.090 0.091 1.6342e-1 1.6343e-1 0.088 0.089 + 2.000000000000000e-02 2.309054339970606e-01 1.984104771646051e-14 2e-2 2.2162e-1 2.2163e-1 0.093 0.094 2.2169e-1 2.2170e-1 0.090 0.091 + 5.000000000000000e-02 6.014214358150070e-01 1.045253287889987e-14 5e-2 5.8044e-1 5.8045e-1 0.091 0.092 5.8046e-1 5.8047e-1 0.091 0.092 + 1.000000000000000e-01 1.226000446492040e+00 3.630493925691041e-15 1e-1 1.1918e+0 1.1919e+0 0.088 0.089 1.1919e+0 1.1920e+0 0.085 0.086 + 1.000000000000000e+00 1.253044200066420e+01 1.230015282154983e-13 1e-0 1.2430e+1 1.2431e+1 0.058 0.059 1.245e+1 1.246e+1 0.05 0.06 + 1.000000000000000e+01 1.256270224354739e+02 6.086298884133178e-15 1e+1 1.2544e+2 1.2545e+2 0.031 0.032 + 4.999999999999999e+01 6.282817756461592e+02 1.506690451671515e-15 5e+1 6.2803e+2 6.2804e+2 0.02 0.03 + diff --git a/libs/ian-presentation.cls b/libs/ian-presentation.cls new file mode 100644 index 0000000..91bd487 --- /dev/null +++ b/libs/ian-presentation.cls @@ -0,0 +1,187 @@ +%% +%% Ian's presentation class +%% + +%% TeX format +\NeedsTeXFormat{LaTeX2e}[1995/12/01] + +%% class name +\ProvidesClass{ian-presentation}[2017/09/29] + +\def\ian@defaultoptions{ + \pagestyle{plain} + \RequirePackage{color} + \RequirePackage{amssymb} +} + +%% paper dimensions +\setlength\paperheight{240pt} +\setlength\paperwidth{427pt} + +%% fonts +\input{size11.clo} +\DeclareOldFontCommand{\rm}{\normalfont\rmfamily}{\mathrm} +\DeclareOldFontCommand{\sf}{\normalfont\sffamily}{\mathsf} +\DeclareOldFontCommand{\tt}{\normalfont\ttfamily}{\mathtt} +\DeclareOldFontCommand{\bf}{\normalfont\bfseries}{\mathbf} +\DeclareOldFontCommand{\it}{\normalfont\itshape}{\mathit} + +%% text dimensions +\textheight=208pt +\textwidth=370pt +\hoffset=-1in +\voffset=-1in +\oddsidemargin=24pt +\evensidemargin=24pt +\topmargin=8pt +\headheight=0pt +\headsep=0pt +\marginparsep=0pt +\marginparwidth=0pt +\footskip=16pt + + +%% remove default skips +\parindent=0pt +\parskip=0pt +\baselineskip=0pt + +%% something is wrong with \thepage, redefine it +\gdef\thepage{\the\c@page} + +%% correct vertical alignment at the end of a document +\AtEndDocument{ + % save total slide count + \immediate\write\@auxout{\noexpand\gdef\noexpand\slidecount{\thepage}} + \vfill + \eject +} + + +%% footer +\def\ps@plain{ + \def\@oddhead{} + \def\@evenhead{\@oddhead} + \def\@oddfoot{\tiny\hfill\thepage/\safe\slidecount\hfill} + \def\@evenfoot{\@oddfoot} +} +\def\ps@empty{ + \def\@oddhead{} + \def\@evenhead{\@oddhead} + \def\@oddfoot{} + \def\@evenfoot{\@oddfoot} +} + + +%% title of slide +\def\title#1{ + \hfil{\bf\large #1}\par + \hfil\vrule width0.75\textwidth height0.3pt\par + \vskip5pt +} + + +%% hyperlinks +% hyperlinkcounter +\newcounter{lncount} +% hyperref anchor +\def\hrefanchor{% +\stepcounter{lncount}% +\hypertarget{ln.\thelncount}{}% +} + +%% define a command and write it to aux file +\def\outdef#1#2{% + % define command% + \expandafter\xdef\csname #1\endcsname{#2}% + % hyperlink number% + \expandafter\xdef\csname #1@hl\endcsname{\thelncount}% + % write command to aux% + \immediate\write\@auxout{\noexpand\expandafter\noexpand\gdef\noexpand\csname #1\endcsname{\csname #1\endcsname}}% + \immediate\write\@auxout{\noexpand\expandafter\noexpand\gdef\noexpand\csname #1@hl\endcsname{\thelncount}}% +} + +%% can call commands even when they are not defined +\def\safe#1{% + \ifdefined#1% + #1% + \else% + {\color{red}\bf?}% + \fi% +} + + +%% itemize +\newlength\itemizeskip +% left margin for items +\setlength\itemizeskip{20pt} +\newlength\itemizeseparator +% space between the item symbol and the text +\setlength\itemizeseparator{5pt} +% penalty preceding an itemize +\def\itemizepenalty{0} +% counter counting the itemize level +\newcounter{itemizecount} + +% item symbol +\def\itemizept#1{ + \ifnum#1=1 + \textbullet + \else + $\scriptstyle\blacktriangleright$ + \fi +} + +\newlength\current@itemizeskip +\setlength\current@itemizeskip{0pt} +\def\itemize{ + \par\penalty\itemizepenalty\medskip\penalty\itemizepenalty + \addtocounter{itemizecount}{1} + \addtolength\current@itemizeskip{\itemizeskip} + \leftskip\current@itemizeskip +} +\def\enditemize{ + \addtocounter{itemizecount}{-1} + \addtolength\current@itemizeskip{-\itemizeskip} + \par\leftskip\current@itemizeskip + \medskip +} +\newlength\itempt@total +\def\item{ + \settowidth\itempt@total{\itemizept\theitemizecount} + \addtolength\itempt@total{\itemizeseparator} + \par + \medskip + \hskip-\itempt@total\itemizept\theitemizecount\hskip\itemizeseparator +} + +%% enumerate +\newcounter{enumerate@count} +\def\enumerate{ + \setcounter{enumerate@count}0 + \let\olditem\item + \let\olditemizept\itemizept + \def\item{ + % counter + \stepcounter{enumerate@count} + % set header + \def\itemizept{\theenumerate@count.} + % hyperref anchor + \hrefanchor + % define tag (for \label) + \xdef\tag{\theenumerate@count} + \olditem + } + \itemize +} +\def\endenumerate{ + \enditemize + \let\item\olditem + \let\itemizept\olditemizept +} + + +%% end +\ian@defaultoptions + +\endinput |