Source code for tdg.references.PRA93023602
import torch
import tdg.conventions
from tdg.references.citation import Citation
citation = Citation(
'Galea, Dawkins, Gandolfi, and Gezerlis, Phys. Rev. A 93, 023602 (2016)',
'Galea:2015vdy')
[docs]def figure_6():
r'''
Figure 6 shows :math:`E/N - \epsilon_b/2` in units of :math:`E_{FG}` as a function of :math:`\log k_F a_{2D}`.
This returns two tensors, the first for the blue circles and the latter for the purple triangles. Each tensor is two rows; the first row being the x-coordinate and the second row the y.
'''
citation('Data lifted from Figure 6')
blue_circles = torch.tensor([
[ +3.0026, +0.7146 ],
[ +2.1522, +0.6029 ],
[ +1.9998, +0.5794 ],
[ +1.4963, +0.4821 ],
[ +1.4442, +0.4678 ],
[ +1.0005, +0.3713 ],
[ +0.7489, +0.3166 ],
[ +0.4954, +0.2640 ],
[ +0.0025, +0.1821 ],
[ -0.5021, +0.1343 ],
[ -1.0024, +0.1092 ],
]).T
purple_triangles = torch.tensor([
[ +3.0061, +0.7150 ],
[ +2.1542, +0.6202 ],
[ +1.4449, +0.5252 ],
[ +0.7498, +0.4580 ],
[ +0.0027, +0.3678 ],
[ -0.5003, +0.3409 ],
[ -0.9991, +0.3333 ],
]).T
return blue_circles, purple_triangles
[docs]def conventional_figure_6():
r'''
The authors only imply the meaning of :math:`\epsilon_b` at the beginning of section VI, where they state the mean-field BCS result :math:`E_{BCS} = E_{FG} + \epsilon_b/2`.
If we use (99) of Ref. :cite:`Beane:2022wcn`, presumably :math:`\epsilon_b/2E_{FG} = -\alpha`.
'''
blue_circles, purple_triangles = figure_6()
blue_circles[0] = -1. / tdg.conventions.from_geometric.log_ka(blue_circles[0])
purple_triangles[0] = -1. / tdg.conventions.from_geometric.log_ka(purple_triangles[0])
return blue_circles, purple_triangles
[docs]def energy_comparison(ax, **kwargs):
blue_circles, purple_triangles = conventional_figure_6()
ax.plot(
blue_circles[0],
blue_circles[1],
marker='o', color='blue', linestyle='none',
label='Optimized Jastrow-BCS [Galea (2016)]'
)
ax.plot(
purple_triangles[0],
purple_triangles[1],
marker='^', color='purple', linestyle='none',
label='Jastrow-Slater [Galea (2016)]'
)