62 lines
1.3 KiB
Plaintext
62 lines
1.3 KiB
Plaintext
{
|
||
"cells": [
|
||
{
|
||
"cell_type": "code",
|
||
"execution_count": null,
|
||
"metadata": {},
|
||
"outputs": [],
|
||
"source": [
|
||
"using Plots, LinearAlgebra\n",
|
||
"include(\"p_space.jl\")"
|
||
]
|
||
},
|
||
{
|
||
"cell_type": "code",
|
||
"execution_count": null,
|
||
"metadata": {},
|
||
"outputs": [],
|
||
"source": [
|
||
"vertices = (0, 1 - 0.5im, 2, 6)\n",
|
||
"subdivisions = 64\n",
|
||
"p, w = get_mesh(vertices, subdivisions)\n",
|
||
"\n",
|
||
"scatter(real.(p), imag.(p))"
|
||
]
|
||
},
|
||
{
|
||
"cell_type": "code",
|
||
"execution_count": null,
|
||
"metadata": {},
|
||
"outputs": [],
|
||
"source": [
|
||
"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",
|
||
"evals = eigen(H).values\n",
|
||
"\n",
|
||
"E_target = 0.258 - 0.164im\n",
|
||
"E = evals[argmin(norm.(evals .- E_target))]\n",
|
||
"\n",
|
||
"print(\"E = $E\")\n",
|
||
"scatter(real.(evals), imag.(evals), xlim = (0,1))"
|
||
]
|
||
}
|
||
],
|
||
"metadata": {
|
||
"kernelspec": {
|
||
"display_name": "Julia 1.9.0",
|
||
"language": "julia",
|
||
"name": "julia-1.9"
|
||
},
|
||
"language_info": {
|
||
"file_extension": ".jl",
|
||
"mimetype": "application/julia",
|
||
"name": "julia",
|
||
"version": "1.9.0"
|
||
}
|
||
},
|
||
"nbformat": 4,
|
||
"nbformat_minor": 2
|
||
}
|