Provide numerical result of s-wave Caroli-de Gennes-Matricon mode (CdGM mode) at $T/T_c=0.3, 0.5, 0.8$. You can get eigenenergy, eigenfunction and pair potential at each temperature.

It is known that there is low-energy excitation levels inside the s-wave vortex core in the superconductor. These states are got by solving following Bogoliubov-de Gennes equation(BdG eq) self-consistently.

$$ \left[-\frac{1}{2k_{F}\xi_{0}}\nabla^{2}-\mu\right]\mathcal{U}{q}(\boldsymbol{r})+\Delta(\boldsymbol{r})\mathcal{V}{q}(\boldsymbol{r}) = E_{q}\mathcal{U}_{q}(\boldsymbol{r}) $$

$$ \left[\frac{1}{2k_{F}\xi_{0}}\nabla^{2}+\mu\right]\mathcal{V}{q}(\boldsymbol{r})+\Delta^{*}(\boldsymbol{r})\mathcal{U}{q}(\boldsymbol{r}) = E_{q}\mathcal{V}_{q}(\boldsymbol{r})
$$

$$ \Delta(\boldsymbol{r})=g\sum_{E_{q}\leq E_{\mathrm{c}}} \mathcal{U}{q}(r)\mathcal{V}{q}^{*}(r)[1-2f(E_{q})] $$

Here, BdG eq is rewritten in dimensionless form using Pippard length $\xi_{0} = \hbar v_{F}/\Delta_{0}$ and zero-temperature bulk gap $\Delta_{0}$. $f(E_{q})$ is Fermi distribution function. Solutions in an isolated vortex, especially CdGM mode is given by following form.

Here, $n$ corresponds to angular momentum quantum number, i.e. CdGM mode is characterized by this number. In this library, the range of $n$ is integers of in $[-70, 69]$. Note that the part of $u_{n}(r), v_{n}(r)$ in the right side of above formula is one of the target of this library, not the left side of it.

https://webpark1378.sakura.ne.jp/