Vibrational quantisation of the B=7 Skyrmion

We consider the inclusion of the most important vibrational modes in the quantisation of the dodecahedral $B=7$ Skyrmion. In contrast to a rigid body quantisation, this formalism allows a spin $\frac{3}{2}$ state to lie below the spin $\frac{7}{2}$ state, in agreement with experimental data. There is also a low lying spin $\frac{1}{2}$ state and two spin $\frac{5}{2} $ states. We find that the excited spin $\frac{7}{2}$ state has a smaller root mean square charge radius than the other states. This prediction is an important signature of the Skyrme model, in conflict with more conventional nuclear models.


Introduction
The Skyrme model is a non-linear field theory of pions which admits soliton solutions called Skyrmions [1]. These are classically stable due to the topology of the system and each Skyrmion has a conserved topological charge, B. After quantisation Skyrmions are identified as nuclei with topological charge equal to baryon number.
The theory is non-renormalisable and so a first principles quantisation is beyond current methods. Instead, one must reduce the degrees of freedom in the problem to a finite number and quantise these. Each charge B Skyrmion may be separated into B charge one Skyrmions. These have six zero modes, three rotations and three translations. Thus to calculate quantities such as the binding energy of a nucleus one should take account of at least 6B degrees of freedom. Unfortunately this means quantising on a 6B dimensional space and little progress has been made, even for B = 2 [2]. Instead, one must select a subset of modes.
The simplest idea is to only include the zero modes of the Skyrmion, those transformations which leave the static energy unchanged. These are the rotations and isorotations (we stay in the centre of mass frame, allowing us to ignore translations). This procedure ignores vibrational modes, dynamical oscillations around the Skyrmion. Zero mode quantisation has had some success, such as reproducing the energy spectra of some light nuclei [3] and a natural description of the Hoyle state [4]. However, there are also some failures. For example, the binding energies are all much too large. This is to be expected when we truncate the degrees of freedom from 6B to 6.
Another failure of zero mode quantisation is the prediction of a spin 7 2 ground state for the 7 Be/ 7 Li isodoublet. The dodecahedral symmetry of the B = 7 Skyrmion rules out low energy states with spin 1 2 , 3 2 and 5 2 . In reality, experimental data show that all these states exist and the ground state has spin 3 2 . The first excited state of 7 Li has spin 1 2 and lies 0.5 MeV above the ground state while the spin 7 2 state is the second excited state lying 4.6 MeV above. In this paper we shall see that the inclusion of vibrational modes in the quantisation procedure resolves this problem.
The 7 Li and 7 Be nuclei are special. Among all nuclei with B < 30 they are the only ones that have an observed spin 7 2 state lying below the lowest spin 5 2 state. The B = 7 Skyrmion is also special. It has the largest finite symmetry group of any known Skyrmion with non-zero pion mass. We shall see that this large symmetry group is the reason why the spin 7 2 state has abnormally low energy. The 7 Li nucleus is usually described using a cluster model [5] which asserts that the nucleus is made of two interacting clusters. These are an alpha particle and a tritium nucleus. This model successfully reproduces the energy spectrum and some electrostatic properties of the nucleus. We shall see that the inclusion of vibrational modes in Skyrmion quantisation highlights a connection between the Skyrme model and the ideas of clustering.
This paper is organised as follows. In section 2 we review the Skyrme model and the structure of the B = 7 vibrational space. We discuss how one should include vibrations in the quantisation procedure and the effects of the Finkelstein-Rubinstein constraints in section 3. Details of the quantisation are laid out in section 4, alongside the results of our calculations and a comparison with the experimental data.
2 The B = 7 Skyrmion and its vibrational space

The Skyrme Model
The Skyrme model can be defined in terms of the three pion fields, π(t, x). These are combined into an SU (2)-valued field where τ are the Pauli matrices and σ is an auxiliary field which satisfies σ 2 + π · π = 1. This ensures that U ∈ SU (2). Many quantities are most easily expressed in terms of the right current R µ = (∂ µ U )U † . The Lagrange density is given by where F π is the pion decay constant, e is a dimensionless parameter and m π is the pion mass. It is more natural to work in Skyrme units. In these, the energy and length units are F π /4e and 2/eF π respectively. The Lagrangian becomes where m = 2m π /eF π is the dimensionless pion mass.
A Skyrmion is a solution of the field equations which minimises the static energy. This is interpreted as the classical mass of the Skyrmion and is given by For this to be finite the Skyrme field must take a constant value, U = 1 2 , at spatial infinity. This one point compactification of space means that U is a map from R 3 ∪ {∞} ∼ = S 3 to SU (2), which is topologically equivalent to S 3 . These maps are labelled by an integer as π 3 (S 3 ) = Z. The integer is identified with the baryon number, B, and can be calculated explicitly from the Skyrme field, where B is the baryon density.
To visualise a Skyrmion we plot a surface of constant baryon density. This is then coloured to express the direction of the pion field,π, as it varies over the surface. We use the same colouring scheme as in [6]. The Skyrmion is coloured white/black whenπ 3 equals ±1 and red, green and blue whenπ 1 + iπ 2 is equal to 1, exp(2πi/3) and exp(4πi/3) respectively.

The vibrational space of the B = 7 Skyrmion
The B = 7 Skyrmion has dodecahedral symmetry as seen in figure 1. There is D 5 symmetry around each face of the Skyrmion and D 3 symmetry around each vertex. These, alongside the additional reflection symmetry, generate the full symmetry group of the Skyrmion Y h .
The vibrational space was numerically generated and studied in [7] by considering small perturbations around the B = 7 Skyrmion. Two low frequency modes were found, one of which had a clear physical interpretation and a clean peak in the power spectrum. We will assume that this is the lowest energy vibrational mode and exclude all others from our analysis. Each point in the vibrational space corresponds to a deformed Skyrme configuration. Our aim is to understand the structure of the space and to find subspaces where the Skyrme configurations have enhanced symmetry.
The vibrational mode we consider has five fold degeneracy and so spans a 5-dimensional vibrational space which we denote V 5 . Each point v ∈ V 5 corresponds to a quadrupole deformation tensor of the Skyrmion, Q(v). There is a natural mapping from a hyperplane in R 6 (isomorphic to V 5 ) to the space of quadrupole tensors. It is where v satisfies (1, 1, 1, 0, 0, 0)·v = 0 to ensure the quadrupole tensor is traceless. We can add vectors on the hyperplane; this is equivalent to adding the quadrupole tensors in R 3 . We choose the normalisation so that a unit vectorv maps to a quadrupole which satisfies Q ij Q ij = 1. Each quadrupole tensor, Q(v), has an associated symmetry group which acts on R 3 . Any symmetry shared by the quadrupole tensor and the B = 7 Skyrmion is a symmetry of the Skyrme configuration at the point v.
In [7] it was found that the vibration we consider preserves the Skyrmion's D 5 symmetry along certain lines in V 5 . Physically, this vibration pulls on two opposite faces of the dodecahedron and breaks the Skyrmion into three clusters: a B = 3 torus sandwiched between two B = 2 tori. This can happen in six ways as there are six pairs of faces on the Skyrmion. Hence there are six special lines in V 5 which preserve D 5 symmetry. They are evenly spaced and are aligned with the vertices of a regular 5-simplex. We must position the 5-simplex in V 5 so that each vertex, v a , maps to a quadrupole tensor which is circle invariant around the axis passing through the Skyrmion faces that are being pulled upon. This ensures that the Skyrme configuration at v a preserves D 5 symmetry. We use the Veronese mapping to help us. This is a map from RP 2 to a 2-dimensional subspace of V 5 . Explicitly it takes This then maps to a quadrupole via (2.6) which is circle invariant around (x 1 , x 2 , x 3 ). For example, the Skyrmion has D 5 symmetry around the axis x 1 = (0, 0, 1). This goes, via the Veronese mapping, to the 6-vector which maps to the quadrupole (2.9) Figure 2: A vibration in V 5 which preserves D 5 symmetry. The parameter λ measures the amplitude of the vibration. This figure was generated using the gradient flow approximation to dynamics. The minimum energy Skyrmion is at λ = 0. This deforms into a ring-like configuration for λ < 0 and three clusters for λ > 0.
This is circle invariant around x 1 as desired. Repeating this process, we may generate the vertices of the 5-simplex in V 5 from the lines which pass through the faces of dodecahedron. This procedure has the corollary that all six vertices of the 5-simplex lie on the 2-dimensional Veronese surface. We denote the 5-simplex vertices as v a ∈ V 5 and the corresponding quadrupole tensors Q a ; these are circle invariant around x a . Any configuration which lies on the line λv a ∈ V 5 , λ ∈ R has D 5 symmetry. The parameter λ is the amplitude of the vibration. For λ > 0 the Skyrmion deforms as described above: a pair of opposite faces are pulled upon. When λ < 0 the faces are pushed together and the Skyrmion flattens out. The full vibration is displayed in figure 2.
We may use the geometry of the 5-simplex to find additional symmetric subspaces in V 5 . The planes passing through an edge of the simplex can be written as where a = b and µ, ν ∈ R. The corresponding quadrupole has a C 2 symmetry, shared with the B = 7 Skyrmion, about the axis x a × x b . This is enhanced to a D 2 symmetry when µ = ν.
The 5-simplex has 20 triangular faces. A line passing through the centre of a face takes the form where a = b = c. In fact, this line passes through two triangular faces which are dual to each other. Thus there are only ten distinct lines. The quadrupole tensor derived from (2.11) has only two distinct eigenvalues. Thus it is circle invariant around the eigenvector of the non-degenerate eigenvalue. This eigenvector passes through a vertex of the B = 7 Skyrmion which has D 3 symmetry. Thus the Skyrme configurations on these 10 lines in V 5 retain D 3 symmetry. Note that, since these quadrupoles are circle invariant, these The analysis so far is based on small perturbations around the Skyrmion. We believe Figure 4: A vibration in V 5 which preserves C 3 symmetry. The parameter λ measures the amplitude of the vibration. The minimum energy Skyrmion is at λ = 0. This deforms into seven individual distorted Skyrmions for λ < 0 and two clusters for λ > 0.
that there will be a bifurcation where the exact symmetries discussed above will break. For example, the D 3 symmetry will break to a C 3 symmetry. This allows the asymptotic configuration in V 5 to be a two-cluster system consisting of a B = 3 and B = 4 Skyrmion. This has lower energy than the three-cluster system described in the previous paragraph. These are hard to distinguish near the origin of V 5 and so the difference will not be apparent in our analysis. Thus we shall assume that the Skyrme configuration along this vibration, at large amplitudes, will be the C 3 symmetric 3+4 cluster configuration instead of the D 3 symmetric 3 + 1 + 3 configuration. The entire vibration is displayed in figure 4.

Quantisation on the vibrational space
Our aim is to quantise the B = 7 Skyrmion taking the lowest energy vibrational mode into account. The manifold we quantise must contain all configurations in V 5 including those generated by rotations and isorotations. Explicitly the manifold is where D is a finite group encoding the dodecahedral symmetry of the Skyrmion. We can think of this manifold as a family of Skyrme configurations parametrised by vibrational (s), rotational (φ, θ, ψ) and isorotational (α, β, γ) coordinates. The angular coordinates are two sets of Euler angles. To quantise we promote all these parameters to dynamical degrees of freedom by allowing them to depend on time. This ansatz allows us to define angular velocities b, isoangular velocities a and find the kinetic energy where g(s) is the metric on N , which depends on the Skyrme configuration at s.
With the kinetic energy written in this way, the quantum kinetic operator is well known [8]. It is proportional to the Laplace-Beltrami operator, ∆. Explicitly where |g| is the determinant of the metric. This preserves the classical symmetries of the kinetic energy after quantisation.
The potential energy, V (s), is the mass (2.4) of the configuration at s. The Hamiltonian on N is To find bound states of definite energy we solve the stationary Schrödinger equation arising from this Hamiltonian, Formally, the wavefunction Ψ is a section of a complex line bundle over N .
There are constraints on Ψ which encode the fact that nucleons are fermions. These are the Finkelstein-Rubinstein (FR) constraints [9]. They can be written in terms of the classical symmetries of the Skyrmion. For example, the B = 7 Skyrmion is invariant under a 2π/5 rotation around the 3-axis followed by a −4π/5 isorotation around the 3axis in isospace. In operator form, this C 5 symmetry puts the following constraint on the where is a vector which passes through one of the dodecahedron's vertices while ξ and n 2 define the isorotation required to return the Skyrmion to its original colouring. The FR signs can be calculated using the rational map ansatz [10]. The procedure is set out in [11].
The constraints (3.6) and (3.7) both apply when the Skyrmion has dodecahedral symmetry. This occurs at the origin of V 5 . For a generic point s there is no symmetry and thus no constraints. In the previous Section we found lines in V 5 which had enhanced symmetry. One set of these preserved D 5 symmetry. Thus, on these lines, only a constraint such as (3.6) applies, as well as an additional constraint which enhances the C 5 symmetry to D 5 . Another set of lines preserved C 3 symmetry, meaning the wavefunction must satisfy a constraint such as (3.7) on these. Now the problem is formulated. To include the lowest vibrational mode when studying the states of 7 Li/ 7 Be we must solve (3.5), an 11-dimensional Schrödinger equation, subject to (3.6) and (3.7) at s = 0, just a constraint such as (3.6) on six lines in the vibrational space (representing the D 5 preserving directions in V 5 ) and a constraint such as (3.7) on ten lines. There are further constraints on the edges of the 5-simplex. To set up and solve this problem rigorously is too hard and so we will make some simplifying assumptions below.
In certain cases the metric, g, will simplify due to the symmetries of the system. In particular the kinetic operator can separate into a part which only acts via the rotational coordinates and a part which only acts via the vibrational ones. We denote this as (3.8) We can then solve the Schrödinger equation using separation of variables where we call u the vibrational wavefunction and Θ the rotational wavefunction. The latter is a solution of the rigid body Schrödinger equation. This problem has been studied extensively, see [3] for details. The solutions are tensor products of Wigner Dfunctions and each solution has six conserved quantities: total spin (J 2 ), total isospin (I 2 ), body fixed spin (L 3 ) and isospin (K 3 ) projections and space fixed spin (J 3 ) and isospin (I 3 ) projections. The space fixed projections do not affect the energy spectrum and as such they are often suppressed in the bra-ket notation where we denote the state One may satisfy the FR constraints using the rotational wavefunction by taking appropriate linear combinations of the Wigner functions. Often the constraints rule out certain spin states. The dodecahedral symmetry of the B = 7 Skyrmion rules out states with J = 1 2 , 3 2 and 5 2 and I = 1 2 . However, the constraints apply to the entire wavefunction Ψ, not just the rotational part Θ. We may alternatively satisfy the FR constraints at the origin of V 5 by insisting that Ψ vanishes there. So there are two ways to satisfy the FR constraints at s = 0: (a) The rotational wavefunction, Θ, is permitted by both the FR constraints. There are no restrictions on the vibrational wavefunction, u(s).
(b) The vibrational wavefunction is zero at the origin, u(0) = 0. There are fewer restrictions on Θ.
These two options also apply on any subspace of vibrational space with enhanced symmetry.
We now look for the low energy states. A rotational wavefunction with spin 7 2 is allowed everywhere in vibrational space. Thus the corresponding wavefunction, Ψ J= 7 2 , can be of type (a). There are no spin 3 2 states allowed at the origin of V 5 . Thus the spin 3 2 wavefunction is of type (b) and u J= 3 2 must vanish at s = 0. We can schematically calculate the energy difference of these states using a harmonic approximation. The spin 7 2 vibrational wavefunction is in the ground state of V 5 while the spin 3 2 wavefunction must be excited in one direction so that it vanishes at the origin. So the spin 3 2 state has one unit of vibrational energy more than the spin 7 2 state. However it will have less rotational energy as the spin is smaller. The ordering of these states depends on the relative energy contributions from vibrations and rotations. We make the approximation that all other vibrations contribute equally to the states. Thus to compare these low lying states we only need to account for the vibration in one direction in V 5 , the smallest energy direction.
In the harmonic approximation the direction of the vibration does not matter as the potential on V 5 is isotropic. However in the full model the direction will be important. A generic direction in V 5 will break the B = 7 Skyrmion into seven B = 1 Skyrmions. This has high potential energy compared to the break up into clusters we saw in Section 2 where the Skyrmion could break into fewer, higher-charge Skyrmions. Thus we believe that the smallest vibration will not be in a generic direction. Instead it will be along one of the high symmetry directions previously discussed.

Quantising the B = 7 Skyrmion
We will now quantise the B = 7 Skyrmion taking a single vibrational mode into account. To do this we must decide on the direction of the vibration in V 5 . We shall assume the lowest energy direction is along either a C 3 preserving line or a D 5 preserving line since these have low energy configurations asymptotically.
The symmetry present on these lines restricts the form of the metric which is 7dimensional. It is standard convention to split the metric into submatrices. We follow [12] and write where U , W and V are 3 × 3 matrices and Λ is a scalar. The kinetic energy is invariant under the action of the symmetry group of the vibration. This restriction means that, along the symmetric lines in V 5 , the cross terms in the metric vanish and the kinetic energy is separable in the sense described in the previous Section. As such the wavefunction takes the form (3.9) on these lines, with the vibrational parameter s now 1-dimensional.
Consider a rotational state with spin J and denote the rotational energy contribution E J (s) so that Note that the rotational energy contribution is a function of s through its dependence on the moments of inertia which vary as the Skyrmion deforms. Then the Schrödinger equation (3.5) reduces to the 1-dimensional equation To solve this we must first generate g(s), V (s) and E J (s). We will do this using gradient flow.
Gradient flow generates a path of steepest descent in field space. We use the separated Skyrmion clusters as initial configurations which are then evolved in gradient flow time τ according to where π are the pion fields and M 7 is the potential energy (2.4). This flow reduces the potential energy of the system and ends at a stationary point of field space. The fields π(τ ) approximate the Skyrme configurations along a half line in V 5 . The solution of (4.4) is beyond analytic calculation and so we must use a numerical code to calculate the flow. The energy V (τ ) and the metric g(τ ) are calculated at numerous points during the process.
The metric at time τ can be expressed in terms of the currents R i = (∂ i U )U −1 and The moments of inertia and Λ are given by Gradient flow time is an unnatural parameter when the Skyrmion clusters are widely separated and near the dodecahedral configuration. Thus, once we have found our quantities numerically we change variables to the geodesic distance, s [13]. This can be defined in terms of the vibrational kinetic energy by demanding which means that There are several advantages to this new coordinate. First, the geodesic distance is related to the cluster separation, r, asymptotically. We can calculate how the moments of inertia vary with r and this gives an asymptotic check of the numerics. Additionally we are able to add an analytic tail to the numerically derived potential and moments of inertia. Further, we may now calculate the harmonic frequency near the origin of V 5 and compare it to what was calculated in [7]. We find the frequency to be 0.34 compared with 0.302 as found in [7]. These are approximately the same, showing the methods are consistent. The small difference is likely due to the different pion masses used. Finally, the new coordinate simplifies the Schrödinger equation (4.3). It now reads From now on, s will refer exclusively to the geodesic distance.

The C 3 direction
The initial configuration for the C 3 direction is constructed using a symmetrised product ansatz of a B = 3 Skyrmion with a B = 4 Skyrmion. These are orientated as in figure 4. The C 3 symmetry constrains the form of the metric. We find that U , V and W are all diagonal. Further We have set Λ = 1 by choosing our parameter to be the geodesic distance.
We now look at specific rotational wavefunctions. Although this direction in V 5 only has C 3 symmetry, it has approximate D 3 symmetry near the origin. This means that a wavefunction disallowed by D 3 symmetry would have extra constraints imposed on it nearby in the full vibrational space. This would increase the energy of the state. Thus we focus on states which are allowed by D 3 symmetry. The rotational wavefunctions we consider are presented in table 1.
Consider the spin 3 2 state. The full wavefunction is of the form The vibrational wavefunction must satisfy u3 2 (0) = 0 so that the full wavefunction Ψ3 2 is consistent with the additional FR constraint (3.6) at s = 0. The derivative must be non-zero here or the vibrational wavefunction will be trivial everywhere. Inserting (4.13) into the Schrödinger equation (3.5) we find that u3 (4.14) We would like to understand the contributions from rotations and vibrations separately.
There is no unique way to split the energy; we choose to define the rotational energy contribution as the rigid body energy of the undeformed Skyrmion, E J (0). We may then split the energy E into three parts: the classical mass of the Skyrmion with each vibrational wavefunction being zero at s = 0, just like the spin 3 2 case. There are three spin 7 2 states. We focus on the lowest energy state. This is a linear combination of the three states,        to ensure that Ψ7 2 satisfies the additional FR constraint (3.6) at s = 0. This gives three uncoupled Schrödinger equations for u, v and w. Generally these three independent equations will not produce a shared energy eigenstate as the effective potential is different for each component rotational wavefunction. However, we can obtain a shared eigenvalue by enforcing an additional boundary condition that the probability distribution is maximal at the origin. This gives This condition, alongside (4.18), produces a discrete eigenvalue spectrum.
The spin 1 2 state is similar to the spin 3 2 state and takes a form analogous to (4.13) with the same conditions on the vibrational wavefunction. However, it has additional constraints in the full vibrational space. It must vanish on the D 5 preserving lines in V 5 due to the FR constraints. The wavefunction we construct is concentrated along a C 3 direction. This direction is maximally far away from the D 5 lines. This can be seen geometrically: the C 3 lines go through the centre of the 5-simplex faces while the D 5 lines pass through the vertices. Thus the wavefunction we construct should already be small on the D 5 lines. A modification is required to make the wavefunction vanish on the D 5 lines which will cost energy. Thus, we expect the true spin 1 2 state to have higher energy than what is calculated here.

Calibration of the model
Before comparing our results to experimental data we must calibrate the model. All previous calibrations are based on zero mode quantisation and as such we don't necessarily expect our choice of parameters to match previous studies. The vibrational energy contribution is of order while the rotational energy contribution is of order 2 . Thus the relative energies of the states will be sensitive to the value of .
In figure 5 the quantum energy of each state is plotted (in Skyrme units) for various values of . The most important feature of the plot is that the spin 7 2 state increases in energy, relative to the other states, as increases. This is because the spin 7 2 state has the largest rotational energy and the smallest vibrational energy; rotational effects dominate for large while vibrational effects dominate for small . To match experimental data the spin 7 2 state must lie between the spin 3 2 state and the first spin 5 2 state. This occurs when 55 < < 65 , (4.20) and as such we demand that lies in this interval. For illustrative purposes we fix = 60. We are left to choose the value of F π , as fixes the dimensionless constant e through the identity = 2e 2 . We will consider two alternative calibrations: Parameter choice (i) gives a good fit to the size of the gaps in the energy spectrum but underestimates the total mass of the Skyrmion compared to the total mass of 7 Li. Choice (ii) gives a reasonable value of the total mass but overestimates the gaps in the spectrum. We have fixed the dimensionless pion mass m to 1 throughout.

Results
We solved the Schrödinger equation (4.15) for all states in discussed in section 4.1. The numerically generated vibrational wavefunction u(s), potential V (s) − M 7 and effective potential V eff (s) for each state is plotted in table 2. We also note the classical mass of the Skyrmion M 7 , the energy contribution from rotations E J (0) and the contribution from vibrations vib , as well as the total energy of each state E. Our results are then compared to experimental data in table 3 for each calibration (i) and (ii).
The results are promising. All of the states considered are seen experimentally. The ordering is correct apart from the spin 1 2 and 3 2 states. We argued earlier that the spin 1 2 state has higher energy than our calculation suggests as it must vanish in a subspace of  the spin 3 2 and 5 2 states. The size of the gaps in the energy spectrum are reasonable for calibration (i) and much too large for calibration (ii). The ratios of the energy gaps between states are independent of this choice, though do depend on . We find that which is reasonably close to the experimental result, 1.44.
The second spin 5 2 state has very high energy. This can be understood by considering the body-fixed spin classically. The highly excited state has |L 3 | = 5 2 . This means that the spin is around the 3-axis. This gives a large energy contribution since the Skyrmion is prolate in this direction. The lower energy spin 5 Table 3: A comparison of the experimentally obtained energy spectrum of 7 Li (column 1) with the results from our calculation using Calibration (i) (column 2) and Calibration (ii) (column 3). The experimental data is from [14].
to rotate about an axis orthogonal to the prolate one. The state we have found probably does not correspond to the experimental state we have compared it to in table 3. In the cluster model [5] this state has a different structure than the others. It is described by a neutron orbiting a 6 Li nucleus. Thus, it could be that we only see this spin state if we include a vibration which can split the B = 7 Skyrmion into these clusters.
The next three experimental states of 7 Li have spin 7 2 , 3 2 and 3 2 . These have natural descriptions in this model. The 7 2 state is orthogonal to the one allowed by the dodecahedron and has a single excited vibration in V 5 . The excited spin 3 2 states will have isospin 3 2 , a possibility we neglected for simplicity. This would also describe the ground states of 7 B and 7 He which have spin 3 2 . The mass of the 7 Li nucleus is 6535 MeV. Calibration (ii) gives the total mass of the ground state to be 6404 MeV which is very close to the experimental value. Calibration (i) gives a much smaller value, only 2764 MeV. There are several ways this could be remedied. First, we have only taken one of the Skyrmion's vibrational modes into account. There are approximately 6B modes, all of which will contribute to the energy. The Casimir energy contribution is also large, a 40% correction in the B = 1 sector [15]. Finally, the Lagrangian may be altered to include a 6th order term which can be chosen to contribute positively to the mass. When this term is the same order as the other terms in the Lagrangian, the Skyrmion solutions do not change significantly [16]. Thus the calculation in this paper would not vary greatly except for the total energy. These three factors could combine to give a reasonable value for the total mass. They also highlight the uncertainty in calculations of total mass in the Skyrme model.
Inclusion of the C 3 vibration has given us a good model of the spin 3 2 , 7 2 and 5 2 states of the 7 Li/ 7 Be isodoublet. Further, it brings us closer to the cluster model of nuclei. This is apparent when we plot the classical baryon density at the maximum of the vibrational wavefunctions which shows us a classical approximation of the quantum state, before rotational averaging. These are plotted in figure 6. We see that the spin 3 2 state exhibits clustering while the spin 7 2 state does not. This goes against conventional wisdom in the cluster model where the ground state is generally the most isotropic.

The D 5 direction
To simplify this calculation we orient the Skyrmion as it is in figure 1. The D 5 symmetry constrains U , V and W to be diagonal with The spin states allowed by D 5 symmetry are presented in table 4. The spin 7 2 state which is allowed by the dodecahedron is given bỹ (4.23) We solve the Schrödinger equation and present the numerically generated vibrational wavefunction u(s), potential V (s) − M 7 and effective potential V eff (s) for each state in table 5. We also note the classical mass of the Skyrmion M 7 , the energy contribution from rotations E J (0) and the contribution from vibrations vib , as well as the total energy of each state E. We see that the C 3 direction produces lower energy states and should be considered the lower energy direction in V 5 . The states arising from the D 5 direction have higher energy than any experimentally discovered state and so are not relevant to the known energy spectrum of the 7 Li/ 7 Be isodoublet. Earlier it was noted that the direction of the vibration in V 5 is important. These results highlight this fact.

The root mean square matter radius
We saw in figure 6 that the different states appear to have different sizes. The simplest quantitative measure of the size of a nucleus is the root mean square (rms) matter radius, r m . We can calculate this for each value of s by taking the square root of where ρ(x, s) is the energy density of the Skyrme configuration at s. For a given state, the rms matter radius is then where we have taken the vibrational wavefunctions to be normalised. We find that the matter radius of the spin 3 2 state, in Skyrme units, is r m 3 2 = 1.85. (4.26) Experiments are unable to measure the matter radius directly. However in most nuclei the matter and charge radii are very similar. Thus we compare (4.26) to the experimentally determined charge radius, 2.444 fm. Our result depends on our choice of F π . Calibration (i) gives a matter radius of 2.22 fm, close to the experimental value. However calibration (ii) gives a very small radius, 0.96 fm. Earlier we found that calibration (i) gave a better match to the energy spectrum. This result adds weight to the idea that it is the better choice. Regardless, ratios of lengths are independent of F π . As such we can compare the matter radii for the spin 7 2 and spin 3 2 states and have more trust in the result. We find that r m 3 2 r m 7 2 = 1.07 . Thus we predict that the ground state of 7 Li is 7% larger than the second excited state, which has spin 7 2 . The rms charge radius of an excited state is difficult to measure experimentally. As such there is no data to confirm our prediction. This is an important signature for the Skyrme model as this prediction is in conflict with the standard cluster model and shell model predictions.

Conclusion and Outlook
In this paper we have considered the inclusion of vibrational modes in the quantisation of the B = 7 Skyrmion. We argued that to understand the low lying states of the 7 Li/ 7 Be isodoublet one can truncate to quantisation along a 1-dimensional line in vibrational space, V 5 . The space has a rich structure best understood using the geometry of a 5simplex. Using this, we picked special directions in the space to quantise along. The calculation gives a reasonable energy spectrum, much closer to the experimental data than had previously been found using zero mode quantisation. Most importantly, the spectrum includes all experimentally seen states and has the spin 7 2 state lying above the spin 3 2 state. During the quantisation procedure some cluster structure emerged. We found that the Skyrmion picks out the C 3 direction as the lowest energy direction. This is remarkable since this is the vibration used in the basic 4 + 3 cluster model. This brings the Skyrme model closer to the cluster models which are used widely in nuclear physics. The advantage of the Skyrme model is that the dynamics of the clusters are fully determined by the Lagrangian. They can merge smoothly into the B = 7 Skyrmion or be infinitely separated; our formalism takes account of all configurations in between.
We predict that the excited spin 7 2 state of 7 Li is smaller than the spin 3 2 ground state. The result depends crucially on the dodecahedral symmetry of the B = 7 Skyrmion. This symmetry appears to persist in modified Skyrme models except in extreme BPS models [17] [18]. Thus this prediction is an important signature for soliton models of finite nuclei.
Vibrational modes have the capacity to fix many issues in the Skyrme model including the high binding energies and small radii found using zero mode quantisation. They also have a fascinating and rich geometric structure. For these reasons alone, more work should be done to understand the vibrational spaces of Skyrmions. It is somewhat surprising that their inclusion leads to a resolution of problems in the B = 7 sector. Hopefully a similar analysis in other sectors can produce more surprises.