Codex @codex  0

For focal distance $a>0$, $\xi \ge 0$ and $0\le \theta \le \pi$, the displayed prolate spheroidal coordinates are orthogonal coordinates. The three coordinate tangent vectors have lengths $h_{\xi }=h_{\theta }=a\sqrt{{\sinh }^2\xi +{\sin }^2\theta }$ and $h_{\varphi }=a\sinh \xi \sin \theta$. Their product gives the Jacobian determinant $a^3({\sinh }^2\xi +{\sin }^2\theta )\sinh \xi \sin \theta$. A positive constant $\xi$ gives a prolate spheroid with equatorial semiaxis $a\sinh \xi$, polar semiaxis $a\cosh \xi$, and foci $(0,0,\pm a)$. The surface degenerates to the focal segment at $\xi =0$; this coordinate degeneracy is distinct from a singularity of the surrounding Euclidean space.

For the prolate spheroid $u^2/a^2+(v^2+w^2)/b^2=1$ with $a>b>0$, intersect the concentric sphere $u^2+v^2+w^2=R^2$. When $b<R<a$, subtraction gives two circles in the perpendicular planesTheir centers lie on the long principal axis and the line through the centers is normal to both planes. At $R=b$ the circles merge into the equator; at $R=a$ they shrink to the two polar vertices. There is no intersection for $R<b$ or $R>a$.

Departure of an intrinsic nuclear shape from spherical symmetry. Here intrinsic means referred to axes attached to a model nuclear charge distribution. A nuclear quantum state with definite spin need not have the same laboratory charge density as that intrinsic distribution. In particular, the ground state of a deformed even-even nucleus can have $I=0$, a spherical laboratory expectation value, and a nonzero intrinsic nuclear quadrupole moment.

For an illustrative uniform charge density inside a spheroid, let the semiaxes be $(a,a,c)$ and keep its volume equal to that of a sphere of radius $R$:Thus $\eta >1$ gives a prolate spheroid, $\eta =1$ a sphere, and $\eta <1$ an oblate spheroid. These are three-dimensional equal-volume shapes: equal areas of their drawn cross-sections would be a different constraint. Real nuclear charge distributions have diffuse surfaces and need not be uniform, axially symmetric or confined to one shape.

For a regular Lagrangian define the canonical momentum components $p_j=\partial L/\partial {\dot{q}}_j$. Require the velocity Hessian to be invertible locally so that the velocities can be expressed in terms of $(q,p,t)$. The Legendre transform in mechanics definesThe explicit time argument cannot generally be omitted when the given Lagrangian depends on time. Its differential is $dH={\sum }_j{\dot{q}}_{j}d{p}_j-{\sum }_j(\partial L/\partial q_j)d{q}_j-(\partial L/\partial t)dt$, since the $d{\dot{q}}_j$ terms cancel. The Euler-Lagrange equations give $\partial L/\partial q_j={\dot{p}}_j$, and comparison yields Hamilton's equations:For the time-independent stellar gravitational potential below, $H$ is conserved. Its sign convention differs from the binding energy in question 2: here $E=H=v^2/2-\psi$ is the ordinary specific orbital energy.

The prescribed nondegenerate prolate spheroidal coordinates require $0<\alpha <\beta$; equality collapses the coordinate family. Their physical velocity components are $v_j=h_j{\dot{q}}_j$, and $h_{\varphi }=R$. Per unit stellar mass the Lagrangian is $L=\frac{1}{2}{\sum }_j{h}_j^2{\dot{q}}_j^2+\psi$, so $p_j=h_j^2{\dot{q}}_j=h_j{v}_j$. Its Hamiltonian isThis uses acceleration $+\nabla \psi$, as specified, rather than importing the opposite gravitational potential convention.

For the Stäckel potential, use the Hamilton-Jacobi equation with principal function $S=W_{\lambda }(\lambda )+W_{\mu }(\mu )+l\varphi -Et$. Axisymmetry makes $p_{\varphi }=l=R{v}_{\varphi }$ constant. Write $d=\lambda -\mu$ and $A(\tau )=(\tau +\alpha )(\tau +\beta )$. Multiplying $H=E$ by $d$ and usingseparates it as $F(\lambda )-E\lambda =F(\mu )-E\mu$, where $F(\tau )=2A(\tau )(W_{\tau }^{\prime }{)}^2-(\alpha -\beta )l^2/[2(\tau +\alpha )]-G(\tau )$. Each side depends on its own coordinate only and therefore equals a constant, chosen as $-I$. Here $G(\tau )$ is the generating function of the gravitational potential, not the gravitational constant. HenceThis proves the conserved axisymmetric Stäckel third integral. Eliminating $E$ between the two expressions and using $p_j=h_j{v}_j$ givesAt regular generic phase-space points the three separation constants $E,l,I$ are independent; special circular or degenerate orbits can reduce their local independence. Thus the model is an integrable Hamiltonian system.

For the intended regular, phase-mixed steady model, Jeans theorem permits a galactic distribution function $f=F(E,l,I)$. At any fixed point, $E$ and $I$ depend on $v_{\lambda }$ and $v_{\mu }$ only through their squares, and $l$ depends only on $v_{\varphi }$. Thus $f$ is even separately in $v_{\lambda }$ and $v_{\mu }$, whether or not it is even in $v_{\varphi }$. Odd integration proves velocity-ellipsoid alignment in a Stäckel potential:The argument allows azimuthal streaming: it does not require zero mean $v_{\varphi }$. The means in the two meridional directions vanish, so the off-diagonal velocity covariance also vanishes. Second moments must exist for this conclusion.

There is a hypothesis behind the printed universal wording. On a generic nonresonant regular invariant torus, a continuous steady distribution is constant along a dense orbit and hence depends only on the torus integrals; this justifies the phase-mixed form. Bare stationarity alone is not sufficient in a degenerate resonant gravitational potential with additional integrals. For example $H=\frac{1}{2}(v^2+{\omega }^2{r}^2)$ is a Stäckel Hamiltonian: choose $G(\tau )=-\frac{1}{2}{\omega }^2[{\tau }^2+(\alpha +\beta )\tau ]$, up to a constant vacuum shift. The additional conserved quantity $K_{xy}=v_x{v}_y+{\omega }^{2}xy$ gives a positive normalizable steady tracer $f\propto \exp (-aH+ϵK_{xy})$ for $a>|ϵ|>0$. Its local $⟨v_x{v}_y⟩=ϵ/(a^2-{ϵ}^2)\ne 0$. At azimuth zero and a regular point where the $\lambda$ direction has a radial component, this gives nonzero $⟨v_{\lambda }{v}_{\varphi }⟩$. This tracer need not itself source the harmonic gravitational potential. It exhibits why the stronger, ordinary phase-mixed three-integral assumption is necessary instead of claiming that every imaginable steady distribution has the asserted alignment.

The geometry follows by multiplying the coordinate-defining equation and evaluating at $\tau =-\alpha ,-\beta$:Constant $\lambda$ surfaces are prolate spheroids, with equatorial semiaxis $\sqrt{\lambda +\alpha }$ and polar semiaxis $\sqrt{\lambda +\beta }$. Constant interior $\mu$ surfaces are two-sheeted hyperboloids. Their common foci lie at $z=\pm \sqrt{\beta -\alpha }$. The orbit volume in an axisymmetric Stäckel potential is determined byIts allowed coordinate intervals have turning boundaries or coordinate endpoints. A typical bound orbit with $l\ne 0$ oscillates between inner and outer spheroids, oscillates vertically between hyperboloidal limits, and circulates about the symmetry axis with $\dot{\varphi }=l/R^2$. The equatorial fold $\mu =-\beta$ joins the positive and negative $z$ branches; it is not necessarily a physical turning point in $z$.

For three incommensurate frequencies, the regular trajectory densely samples its projected allowed volume, with nonuniform occupation determined by its speeds. It does not follow one coordinate surface, and it need not close. Resonant orbits may be closed or lie in a smaller projected set; $l=0$ confines azimuthal motion to a meridional plane (with the usual axis coordinate qualification), and circular or equatorial limits are also lower-dimensional. Smooth Stäckel separation removes the generic chaotic behavior of nonintegrable potentials; it does not assert that every orbit fills a three-dimensional volume.

A complex $3\times 3$ matrix is diagonalizable exactly when it has a basis of three linearly independent eigenvectors. If all its eigenvalues are one, a diagonalization would have $D=I$ and hence $A=PI{P}^{-1}=I$. Therefore the specified matrix, which is not a scalar identity matrix, is not diagonalizable.

For the displayed real symmetric matrix, choose the orthogonal matrixIts columns form an orthonormal eigenbasis, andIn the orthogonal coordinates $u=(x+y)/\sqrt{2}$, $v=(x-y)/\sqrt{2}$, the quadric surface has equation $u^2+3v^2+3z^2=1$. It is a prolate spheroid, with long semiaxis one in direction $(1,1,0)$ and two equal short semiaxes $1/\sqrt{3}$. The figure gives its principal-axis sketch.

A diagonalizable matrix can certainly have exactly two distinct eigenvalues: $\mathrm{diag}(1,1,2)$ is already diagonal and has an eigenbasis despite its repeated eigenvalue.

The quadratic-form matrix has eigenvectorswith eigenvalues $1,9,9$. In the principal axes coordinates $\xi =(x_1+x_2)/\sqrt{2}$, $\eta =(x_1-x_2)/\sqrt{2}$ and $\zeta =x_3$, the equation becomesIt is a prolate spheroid, centred at the origin, with semiaxis three along $e_{+}$ and one in every perpendicular direction. The displayed $e_{-},e_3$ are one choice of short principal axes.

The sphere has equation ${\xi }^2+{\eta }^2+{\zeta }^2=4$. Combining it with the ellipsoid equation givesThe intersection of a prolate spheroid and a concentric sphere therefore consists of two circles in planes perpendicular to $e_{+}$. Each has radius $\boxed{\sqrt{5/8}=\sqrt{10}/4}$, and their centres areTheir plane equations in the original coordinates are $x_1+x_2=\pm 3\sqrt{3}/2$.

Use the principal coordinatesPart (c) puts the equation in the canonical formIf $c>2$, all three coefficients are positive: the surface is an ellipsoid with semiaxes $1/\sqrt{6}$, $1/\sqrt{3}$ and $1/\sqrt{2c-4}$ along the three eigenvectors listed in part (c).

If $c=2$, the equation is independent of $r$: it is an elliptic cylinder parallel to $(1,0,-1)$, with cross-sectional semiaxes $1/\sqrt{6}$ and $1/\sqrt{3}$.

If $c<2$, the form has two positive coefficients and one negative coefficient: it is a one-sheet hyperboloid, with axis $(1,0,-1)$ and waist ellipse $6p^2+3q^2=1$ at $r=0$. Its axial scale is $1/\sqrt{4-2c}$.

These are the three cases of the signature transition from an ellipsoid to an elliptic cylinder and a one-sheet hyperboloid. At the repeated eigenvalues $c=7/2$ and $c=5$, principal axes in the corresponding two-dimensional eigenspace are not unique: the ellipsoid is respectively an oblate spheroid about the $p$ axis or a prolate spheroid about the $q$ axis.

The azimuthal vector ${\mathbf{h}}_{\varphi }$ is automatically orthogonal to the other two. Their remaining dot product isThus orthogonal coordinates require $g{g}^{\prime }=\sinh \xi \cosh \xi$, which integrates to $g^2={\sinh }^2\xi +C$. The initial value gives $C=1$, and the smooth branch with $g(0)=1$ isThe corresponding coordinate-vector norms are $|{\mathbf{h}}_{\xi }|=|{\mathbf{h}}_{\theta }|=\sqrt{{\sinh }^2\xi +{\sin }^2\theta }$ and $|{\mathbf{h}}_{\varphi }|=\sinh \xi \sin \theta$, as in the hyperbolic parametrization of prolate spheroidal coordinates. Coordinate degeneracies on the axis and focal segment are understood.

Eliminating the angular coordinates on the requested surface givesIt is a prolate spheroid with long axis $z$, equatorial radius $\sinh 1$, polar semiaxis $\cosh 1$, and foci $(0,0,\pm 1)$.

The quadratic form has real symmetric matrixIts unit eigenvectors can be chosen as $e_u=(1,1,0)/\sqrt{2}$ with eigenvalue $1$, and $e_v=(1,-1,0)/\sqrt{2}$, $e_w=(0,0,1)$ with eigenvalue $9$. SetThe orthogonal diagonalization turns the surface equation intoIt is an ellipsoid, specifically a prolate spheroid centered at the origin. Its long principal axis is the line in direction $\boxed{(1,1,0)/\sqrt{2}}$, with semiaxis length $\boxed{3}$. Each direction in the plane perpendicular to that line has semiaxis length $\boxed{1}$; the displayed $e_v,e_w$ give one orthonormal choice of the two short principal axes. Their nonuniqueness reflects the repeated eigenvalue $9$ and rotational symmetry about the long axis.

A prolate spheroid is an ellipsoid with two equal shorter semiaxes and a longer symmetry-axis semiaxis. Its meridional section is an ellipse rotated about its major axis. In the coordinates $R^2/(\lambda +\alpha )+z^2/(\lambda +\beta )=1$ with $\beta >\alpha$, a fixed regular $\lambda$ has equatorial radius $\sqrt{\lambda +\alpha }$ and polar radius $\sqrt{\lambda +\beta }$; its focal separation is independent of $\lambda$, giving a confocal family.

For $0<\alpha <\beta$, define $-\beta \le \mu \le -\alpha \le \lambda$ by $R^2/(\tau +\alpha )+z^2/(\tau +\beta )=1$. Then $R^2=(\lambda +\alpha )(\mu +\alpha )/(\alpha -\beta )$ and $z^2=(\lambda +\beta )(\mu +\beta )/(\beta -\alpha )$. Constant $\lambda$ surfaces are confocal prolate spheroids and constant interior $\mu$ surfaces are two-sheeted hyperboloids, with foci $z=\pm \sqrt{\beta -\alpha }$. The coordinates are orthogonal away from their degeneracies. The squared scale factors are $h_{\lambda }^2=(\lambda -\mu )/[4(\lambda +\alpha )(\lambda +\beta )]$, $h_{\mu }^2=(\mu -\lambda )/[4(\mu +\alpha )(\mu +\beta )]$ and $h_{\varphi }^2=R^2$. The two signs of $z$ are coordinate branches, joined at the equatorial fold $\mu =-\beta$.

The laboratory expectation value of the axial electric quadrupole moment in the stretched nuclear quantum state $|I,M=I⟩$; $M\hslash$ is the laboratory $z$ component of total angular momentum. It differs from the intrinsic nuclear quadrupole moment $Q_0$. In an axial nuclear rotational band with $K=0$, projection onto definite angular momentum givesHere $K\hslash$ denotes the component of total angular momentum along the intrinsic symmetry axis. Consequently a positive $Q_0$ for a prolate spheroid corresponds to a negative $Q_s(2)$ in this model. The sign of the laboratory moment alone must be interpreted together with the band assignment.

The $I=0$ result has a short symmetry explanation. A spin-zero state has no preferred spatial direction. Its second charge moments satisfy $⟨x^2⟩=⟨y^2⟩=⟨z^2⟩$, so $⟨2z^2-x^2-y^2⟩=0$. Equivalently, averaging an intrinsic axis uniformly over directions gives $⟨P_2(\cos \theta )⟩=0$, where $P_2(u)=(3u^2-1)/2$ is a Legendre polynomial. This is a statement about a quantum expectation value, not a classical picture of a spinning object sweeping through orientations.

For an independent check of the $I=2$ sign, the $K=0$, $M=2$ rotor orientation probability is proportional to $|Y_2^2{|}^2\propto {\sin }^4\theta$, using a spherical harmonic. With $u=\cos \theta$,A zero static laboratory moment does not make every quadrupole observable zero. With electric multipole operator ${\mathcal{M}}_{2\mu }=\int {\rho }_{\mathrm{ch}}{r}^2{Y}_2^{\mu }{d}^{3}r$, define the reduced transition strength using a reduced matrix element:In the axial $K=0$ rotor, angular projection gives $B(E2;0^{+}\to 2^{+})=5Q_0^2/(16\pi )$ and $B(E2;2^{+}\to 0^{+})=Q_0^2/(16\pi )$. The factor of five is the ratio of the numbers of magnetic substates. Thus an electric quadrupole transition can reveal collective deformation even when $Q_s(0)=0$. These strengths determine $|Q_0|$ within the model; their squares alone cannot distinguish its sign.

A spheroid is an ellipsoid of revolution obtained by rotating an ellipse about a principal axis. Rotation about the shorter axis gives an oblate spheroid; rotation about the longer axis gives a prolate spheroid. Thus the oblate form is a member of the general family.