From d6b702581d4e50b4a16889d1448f3c431cbe38a5 Mon Sep 17 00:00:00 2001
From: ysyapa <ysyapa@ncsu.edu>
Date: Tue, 23 Jan 2024 18:58:03 +0000
Subject: [PATCH] Define potentials in p_space.jl

---
 EC_test.ipynb         | 5 +----
 p_space.jl            | 4 ++++
 p_space_test.ipynb    | 1 -
 simple_berggren.ipynb | 5 +----
 4 files changed, 6 insertions(+), 9 deletions(-)

diff --git a/EC_test.ipynb b/EC_test.ipynb
index 061d003..e442202 100644
--- a/EC_test.ipynb
+++ b/EC_test.ipynb
@@ -7,10 +7,7 @@
    "outputs": [],
    "source": [
     "using Plots, LinearAlgebra\n",
-    "include(\"p_space.jl\")\n",
-    "\n",
-    "g0(R, p, q) = (exp(-(1/4)*(p + q)^2*R^2)*(-1 + exp(p*q*R^2))*R)/(2*sqrt(π))\n",
-    "g1(R, p, q) = (exp(-(1/4)*(p + q)^2*R^2)*(2 + p*q*R^2 + exp(p*q*R^2)*(-2 + p*q*R^2)))/(2*p*sqrt(π)*q*R)"
+    "include(\"p_space.jl\")"
    ]
   },
   {
diff --git a/p_space.jl b/p_space.jl
index 3915930..2e5e2de 100644
--- a/p_space.jl
+++ b/p_space.jl
@@ -1,5 +1,9 @@
 using FastGaussQuadrature
 
+# Gaussian potentials in momentum space
+g0(R, p, q) = (exp(-(1/4)*(p + q)^2*R^2)*(-1 + exp(p*q*R^2))*R)/(2*sqrt(π))
+g1(R, p, q) = (exp(-(1/4)*(p + q)^2*R^2)*(2 + p*q*R^2 + exp(p*q*R^2)*(-2 + p*q*R^2)))/(2*p*sqrt(π)*q*R)
+
 function gausslegendre_shifted(a, b, n)
     scale = (b - a) / 2
     shift = (a + b) / 2
diff --git a/p_space_test.ipynb b/p_space_test.ipynb
index 6bdbcd0..b30664e 100644
--- a/p_space_test.ipynb
+++ b/p_space_test.ipynb
@@ -30,7 +30,6 @@
    "outputs": [],
    "source": [
     "# ResonanceEC: Eq. (21)\n",
-    "g1(α, p, q) = sqrt(α/(4*π))*((1/α+2/(p*q))*exp(-(p+q)^2/(4*α))+(1/α-2/(p*q))*exp(-(p-q)^2/(4*α)))\n",
     "V_pq(p, q) = -10 * g1(1, p, q)\n",
     "\n",
     "H = get_H_matrix(V_pq, p, w)\n",
diff --git a/simple_berggren.ipynb b/simple_berggren.ipynb
index 4e4e1d6..0608ce5 100644
--- a/simple_berggren.ipynb
+++ b/simple_berggren.ipynb
@@ -7,10 +7,7 @@
    "outputs": [],
    "source": [
     "using Plots, LinearAlgebra\n",
-    "include(\"p_space.jl\")\n",
-    "\n",
-    "g0(R, p, q) = (exp(-(1/4)*(p + q)^2*R^2)*(-1 + exp(p*q*R^2))*R)/(2*sqrt(π))\n",
-    "g1(R, p, q) = (exp(-(1/4)*(p + q)^2*R^2)*(2 + p*q*R^2 + exp(p*q*R^2)*(-2 + p*q*R^2)))/(2*p*sqrt(π)*q*R)"
+    "include(\"p_space.jl\")"
    ]
   },
   {