## Nuclear deformation at finite temperature

Y. Alhassid, 1 C.N. Gilbreth, 1 and G.F. Bertsch 2 1 Center for Theoretical Physics, Sloane Physics Laboratory, Yale University, New Haven, CT 06520 2 Department of Physics and Institute of Nuclear Theory, Box 351560 University of Washington, Seattle, WA 98915

(Dated: August 4, 2014)

Deformation, a key concept in our understanding of heavy nuclei, is based on a mean-field description that breaks the rotational invariance of the nuclear many-body Hamiltonian. We present a method to analyze nuclear deformations at finite temperature in a framework that preserves rotational invariance. The auxiliary-field Monte-Carlo method is used to generate the statistical ensemble and calculate the probability distribution associated with the quadrupole operator. Applying the technique to nuclei in the rare-earth region, we identify model-independent signatures of deformation and find that deformation effects persist to higher temperatures than the spherical-todeformed shape phase-transition temperature of mean-field theory.

PACS numbers: 21.60.Cs, 21.60.Ka, 21.10.Ma, 02.70.Ss

Motivation .- Mean-field theory is a useful method for studying correlated many-body systems. However, it often breaks symmetries, making it difficult to compare its results with physical spectra that preserve these symmetries. In addition, although mean-field theory often predicts sharp phase transitions at finite temperature, they are washed out in finite-size systems. The challenge is to find tools to study the properties of finite-size systems within a framework that preserves the underlying symmetries while also allowing calculation of the quantities that describe symmetry breaking in mean-field theory.

In nuclear physics, this issue is especially important in the understanding of heavy deformed nuclei, which are of wide experimental and theoretical interest. The current theory of these nuclei is based on self-consistent mean-field (SCMF) theory, which predicts both spherical and deformed ground states [1] depending on the nucleus. SCMF is a convenient tool to study their intrinsic structure but it breaks rotational invariance, a prominent symmetry in nuclear spectroscopy. The occurrence of large deformations in the ground state and at low excitations gives rise to rotational bands and large electric quadrupole transition intensities between states within the bands. At higher excitations, much less is known experimentally. Characterization of this part of the spectrum is needed for accurate calculation of the nuclear level density, which is very sensitive to deformation and other structure effects; observed level densities in rare-earth nuclei at the neutron evaporation threshold vary by more than an order of magnitude [3]. In addition, nuclear fission is a phenomenon of shape dynamics, and calculation of fission rates for excited nuclei requires their level densities as a function of deformation [4].

Here we investigate nuclear deformation at finite temperature using the auxilliary-field Monte Carlo (AFMC) method, which is well suited to the study of the evolution of nuclear properties with excitation energy while preserving rotational invariance. In particular, we cal- culate the distribution of the quadrupole operator in the lab frame and demonstrate that it exhibits modelindependent signatures of deformation. We use moments of this distribution to calculate rotationally invariant observables, which allow us to extract effective values of the intrinsic deformation and its fluctuations. Deformations have been studied previously by the AFMC method, but with an ad hoc prescription to extract the intrinsic-frame properties [5]. The methods presented here should be applicable to other finite-size systems in which correlations beyond the mean field are important.

Methodology .- Formally, we can examine the statistical characteristics of nuclei at finite excitations by calculating the thermal expectation values of observables ˆ O associated with the property of interest, 〈 ˆ O 〉 = Tr( ˆ Oe -βH ˆ ) / Tr e -βH ˆ . Here β -1 is the temperature and ˆ H is the Hamiltonian, which we assume to be rotationally invariant. We denote operators in the many-particle space with a circumflex, to be distinguished from operators in the single-particle space, which are ordinary matrices, denoted by bold-face symbols. Also, we denote the trace over the full many-particle Fock space as Tr and the trace of matrices in the single-particle space by tr. The probability distribution of an operator ˆ , O P β ( o ) = Tr( δ O [ ˆ -o e ) -βH ˆ ] / Tr e -βH ˆ can be calculated using the Fourier representation of the δ function:

$$\underset { \text{old} } { \text{is} } \quad P _ { \beta } ( o ) = \frac { 1 } { \text{Tr} \, e ^ { - \beta \hat { H } } } \int _ { - \infty } ^ { \infty } \frac { d \varphi } { 2 \pi } e ^ { - i \varphi o } \, \text{Tr} \, \left ( e ^ { i \varphi \hat { O } } e ^ { - \beta \hat { H } } \right ). \quad ( 1 ) \\ \cdot \quad \quad - \quad \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \,.$$

Eq. (1) is well-known for one-body observables ˆ O that commute with the Hamiltonian, e.g., the number operator and the z -component of the angular momentum [6].

Nuclear shape is different in that the relevant operators, e.g., the quadrupole operators, do not commute with the Hamiltonian. Nevertheless, it is possible with Eq. (1) to define the distribution of quantum-mechanical observables that carry information about deformation as

well as energy. The distribution (1) can be expressed in terms of the many-particle eigenstates of ˆ and O ˆ H as

$$P _ { \beta } ( o ) = \sum _ { n } \delta ( o - o _ { n } ) \sum _ { m } \langle o, n | e, m \rangle ^ { 2 } e ^ { - \beta e _ { m } }. \quad ( 2 ) \quad \ e ^ { - \beta \tilde { h } } \\ \quad \text{$\text{$\text{$\text{$\text{$nron}$} } }. \quad. \quad \text{$\text{$\text{$\text{$nron}$} } }.$$

Here | o, n 〉 are eigenstates of ˆ O satisfying ˆ O o, n | 〉 = o n | o, n 〉 and similarly for | e, m 〉 . Eq. (2) is valid whether or not the operators ˆ and O ˆ H commute. When they do commute, they share a common basis of eigenstates such that 〈 o, n e, m | 〉 = δ m,n and the distribution (2) reduces to its more familiar form P β ( o ) = ∑ n δ o ( -o n ) e -βe n . Note that in a finite model space the eigenvalues o n form a discrete set and P β ( o ) is a finite sum of δ functions.

In this work we consider the observable ˆ O to be the spectroscopic mass quadrupole operator ˆ Q 20 = ∑ ( i 2 z 2 i -x 2 i -y 2 i ) where the sum is taken over all nucleons. The probability distribution P β ( q ) of ˆ Q 20 is defined as in Eq. (1) with ˆ = O ˆ Q 20 and o = . q

As we will show, this distribution can be accurately computed by the AFMC method. However, the intrinsicframe properties are not directly accessed by the operator ˆ Q 20 , which is a laboratory-frame observable. We shall demonstrate in this work that nevertheless the distribution P β ( q ) is sensitive to deformation effects and that the main properties of the deformation in the intrinsic frame can be recovered from moments of this distribution.

Intrinsic frame quantities may be defined in terms of the expectation values of rotationally invariant combinations of the quadrupole tensor operator ˆ Q 2 µ ( µ = -2 , . . . , 2) [7, 8]. The lowest-order invariant is quadratic, ˆ Q · ˆ Q = ∑ µ ( -) µ ˆ Q 2 µ ˆ Q 2 -µ . There is one third-order invariant defined by coupling three quadrupole operators to angular momentum zero, ( ˆ Q × ˆ ) Q · ˆ Q = √ 5 ∑ µ ,µ ,µ 1 2 3 ( 2 2 2 µ 1 µ 2 µ 3 ) ˆ Q 2 µ 1 ˆ Q 2 µ 2 ˆ Q 2 µ 3 . The fourthand fifth-order invariants are also unique [9] and we define them as ( ˆ Q Q · ˆ ) 2 and ( ˆ Q Q · ˆ )(( ˆ Q × ˆ ) Q Q · ˆ ), respectively. When the invariant is unique at a given order, its expectation value can be computed directly from the lab-frame moments of ˆ Q 20 , defined by 〈 ˆ Q n 20 〉 β = ∫ q n P β ( q dq ) . The conversion factors are given in Table I.

TABLE I: First line: the ratio of the expectation value of the invariant of order n (see text) to the n -th moment of ˆ Q 20 . Second line: the n -th moment of ˆ Q 20 for the rigid rotor in units of q n 0 ( q 0 is the rotor's intrinsic quadrupole moment).

| n         | 2   | 3                | 4    | 5                       |
|-----------|-----|------------------|------|-------------------------|
| invariant | 5   | - 5(7 / 2) 1 / 2 | 35/3 | - (11 / 2)(7 / 2) 1 / 2 |
| rotor     | 1/5 | 2/35             | 3/35 | 4/77                    |

AFMC .- We shall use the AFMC to evaluate the distribution in Eq. (1) for ˆ = O ˆ Q 20 . AFMC is arguably the most powerful computational tool for finding the ground states and thermal properties in large-dimension manyparticle spaces. It is based on the Hubbard-Stratonovich representation [10] of the imaginary-time propagator, e -βH ˆ = ∫ D σ G U [ ] σ ˆ σ , where D σ [ ] is the integration measure, G σ ( ) is a Gaussian weight, and ˆ U σ is a one-body propagator of non-interacting nucleons moving in auxiliary fields σ . Practical implementations require that the Hamiltonian be restricted to oneand two-body terms, and that the two-body terms have the so-called good sign [11]. The method has been applied to nuclei in the framework of the configuration-interaction shell model [12-14], where it is called the shell-model Monte Carlo (SMMC). It has been particularly successful in calculating statistical properties of nuclei such as level densities [15]. The distribution of ˆ Q 20 is obtained from the Monte Carlo sampling of fields σ as a ratio of averages

$$\overset { \nearrow } { 1 } \\ y$$

Here 〈 X 〉 W = ∫ D σ W X / [ ] σ σ ∫ D σ W [ ] σ , where W σ = G σ | Tr ˆ U σ | is used for the Monte Carlo sampling and Φ σ = Tr ˆ U / σ | Tr ˆ U σ | is the Monte Carlo sign function.

For a given ˆ U σ , we carry out the ˆ Q 20 projection using a discretized version of the Fourier decomposition in Eq. (1). We take an interval [ -q max , q max ] and divide it into 2 M +1 equal intervals of length ∆ q = 2 q max / (2 M + 1). We define q m = m q ∆ , where m = -M,...,M , and approximate the quadrupole-projected trace in (3) by

$$\overset { \mathbf V ^ { \mu } } { \text{lract} }, & \overset { \mathbf V } { \text{lract} }, \\ \overset { \mathbf V } { \text{order} } & \overset { \mathbf V } { \text{Tr} } \left ( \delta ( \hat { Q } _ { 2 0 } - q _ { m } ) U _ { \sigma } \right ) \approx \frac { 1 } { 2 q _ { \max } } \sum _ { k = - M } ^ { M } e ^ { - i \varphi _ { k } q _ { m } } \text{Tr} ( e ^ { i \varphi _ { k } \hat { Q } _ { 2 0 } } \hat { U } _ { \sigma } ) \, \\ \overset { \mathbf V } { \text{oper} } & \overset { \mathbf V } { \text{where} } \ \overset { \mathbf V } { \text{$\mathbf V} }, \ - \ \pi k \cdot / \alpha \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

We found the thermalization of ˆ Q n 20 to be slow with the pure Metropolis sampling. This can be overcome by augmenting the Metropolis-generated configurations by rotating them through a properly chosen set of N Ω rotation angles Ω. In practice, it is easier to rotate the observables, i.e., we replace 〈 e iϕQ ˆ 20 〉 σ by 1 N Ω ∑ j 〈 e iϕQ ˆ 20 (Ω ) j 〉 σ . ˆ ˆ ˆ ˆ 1 ˆ being the rotation op- where ϕ k = πk/q max ( k = -M,...,M ). Since ˆ Q 20 is a one-body operator and ˆ U σ is a one-body propagator, the Fock space many-particle traces on the r.h.s. of Eq. (4) reduce to determinants in the single-particle space Tr ( e iϕ k ˆ Q 20 ˆ U σ ) = det ( 1 + e iϕ k Q 20 U σ ) . Here Q 20 and U σ are the matrices representing, respectively, ˆ Q 20 and ˆ U σ , in the single-particle space. In practice, projections are carried on the neutron and proton number operators as well to fix the Z and N of the ensemble [14].

Here Q 20 (Ω) = RQ 20 R -with R erator for angle Ω. Details will be given elsewhere.

We next discuss a few simple examples that can be treated analytically or nearly so.

Rigid rotor .- As a first simple example, we consider an axially symmetric rigid rotor with an intrinsic quadrupole

FIG. 1: The ground-state distribution P g s . . ( q ) vs q/q 0 for a prolate rotor with intrinsic quadrupole moment q 0 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0081v1-figure-65372637-b99f-4a4e-94af-b537d6f5e338.jpg)

moment q 0 in its ground state. The distribution of its spectroscopic quadrupole operator in the laboratory frame Q 20 = q 0 (3 cos 2 θ -1) / 2 can be calculated in closed form. For a prolate rotor ( q 0 &gt; 0)

$$P _ { g. s. } ( q ) = \begin{cases} \begin{array} { c c } \left ( \sqrt { 3 } q _ { 0 } \sqrt { 1 + 2 \frac { q } { q _ { 0 } } } \right ) & \text{for } - \frac { q _ { 0 } } { 2 } \leq q \leq q _ { 0 } \\ 0 & \text{otherwise } \end{array} \cdot ( 5 ) & \text{groun } \\ \text{for } \mu \\ \text{and} \end{cases}$$

This distribution is shown in Fig. 1. The oblate rotor ( q 0 &lt; 0) distribution is obtained from (5) by replacing q with -q and q 0 with | q 0 | . The moments of the distribution (5) can be calculated from a simple recursion relation; their values for 2 ≤ n ≤ 5 are given in Table I.

20 Ne .-As a simple illustration in nuclear spectroscopy, we consider the light deformed nucleus 20 Ne. The orbital part of the single-particle wave functions are taken to be the states of the N = 2 harmonic oscillator shell, i.e., the sd -shell. The single-particle eigenvalues of Q 20 are -2, 1, and 4 (in units of b 2 [16]) with degeneracies of 6, 4 and 2, respectively. The many-particle eigenvalues of ˆ Q 20 for 20 Ne in the valence sd -shell thus range from -8 to 16 with a uniform spacing of 3. The distribution P β ( q ) at β = 0 is just the distribution of these eigenvalues.

We have used this nucleus as a simple test of the AFMC. Here we take the single-particle energies according to the USD interaction [17] and consider an attractive quadrupole-quadrupole interaction -χQ ˜ · ˜ , Q with ˜ Q 2 µ = ∑ i r 2 i Y 2 µ (ˆ ) r i and χ = 8 π 5 38 5 . A 5 / 3 MeV /b 4 [18]. In Fig. 2 we show the quadrupole distribution of the 20 Ne ground state. The discrete nature of the many-particle eigenvalues of ˆ Q 20 is evident; the distribution is a set δ functions at integers -8 , -5 , . . . , 13 16. The envelope of , the strengths has the skewed shape that looks qualitatively similar to the prolate rigid-rotor distribution.

SCMF .- It is instructive to compare our results with those of the thermal SCMF, e.g., the finite-temperature Hartree-Fock-Bogoliubov (HFB) approximation. The HFB solution is characterized by temperature-dependent one-body density matrix ρ β and pairing tensor κ β . In general, two types of phase transitions can occur vs

FIG. 2: The AFMC ground-state quadrupole distribution P g s . . ( q ) for 20 Ne. The sharp δ -like peaks demonstrate the discrete nature of the spectrum of ˆ Q 20 and their envelope resembles the prolate rigid-rotor distribution in Fig. 1.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0081v1-figure-7ab77bc5-c496-4fd1-8aaf-38293e563934.jpg)

/negationslash temperature, a pairing transition and a deformed-tospherical shape transition [19-21]. A shape phase transition is also the generic result of a Landau theory in which the order parameter is the quadrupole deformation tensor [22]. The vast majority of deformed HFB ground states are axially symmetric [23], i.e., 〈 ˆ Q 2 µ 〉 = 0 for µ = 0. The second-order invariant 〈 ˆ Q · ˆ Q 〉 may be calculated in HFB by using Wick's theorem

$$\intertext { t o r } \underset { \mu } { \mathrm i n g } \quad \langle \hat { Q } \cdot \hat { Q } \rangle = Q _ { 0 } ^ { 2 } + \sum _ { \mu } ^ { \kappa } ( - ) ^ { \mu } \text{tr} \left [ Q _ { 2 \mu } \left ( 1 - \rho _ { \beta } \right ) \mathbf Q _ { 2 - \mu } \rho _ { \beta } \right ] \\ \intertext { s i o n } \mathrm i. \quad \ \ \ \ + \sum _ { \mu } ( - ) ^ { \mu } \text{tr} \left [ Q _ { 2 \mu } \kappa _ { \beta } \, \mathbf Q _ { 2 - \mu } ^ { T } \kappa _ { \beta } ^ { * } \right ] \, \ \ ( 6 ) \\ \intertext { e c - } \mathrm i \text{Ne. } \quad \text{ where } Q _ { 0 } \equiv \text{tr} ( Q _ { 2 0 } \rho _ { \beta } ) \text{ is the intrinsic axial quadrupole}$$

where Q 0 ≡ tr( Q 20 ρ β ) is the intrinsic axial quadrupole moment. The remaining terms on the r.h.s. of (6) represent the contributions due to quantal and thermal fluctuations. We shall compare our AFMC results for rareearth nuclei with the HFB theory in the next section.

Rare-earth nuclei .- Here we present results for rareearth nuclei. The single-particle orbitals are taken from a Woods-Saxon potential plus spin-orbit interaction; they span the 50 -82 shell plus 1 f 7 / 2 orbital for protons and the 82 -126 shell plus 0 h 11 / 2 , 1 g 9 / 2 orbitals for neutrons. We use the same interaction as in Refs. [24, 25]. The quadrupole moments are scaled by a factor of 2 to account for the model space truncation.

We first examine 154 Sm, a strongly deformed nucleus with an intrinsic quadrupole moment of Q 0 ∼ 1600 fm , as determined experimentally from in-band electric 2 quadrupole transitions [26]. AFMC P β ( q ) distributions are shown in Fig. 3 at three temperatures. The distributions appear continuous because the many-particle eigenvalues of ˆ Q 20 are closely spaced. At the lowest temperature of T = 0 1 MeV (bottom panel), . e -βH ˆ effectively projects out the ground-state band. We observe the characteristic skewed distribution of the prolate rotor. The dashed line is the rotor distribution (5) with q 0 taken at the HFB value of Q 0 . The middle panel is the distribution at the HFB shape transition temperature, T = 1 14 .

FIG. 3: Probability distributions P β ( q ) for 154 Sm at T = 0 1 MeV, . T = 1 14 MeV (shape transition temperature) and . T = 4 MeV. The low-temperature distribution is compared with the rigid-rotor distribution (dashed line) and reflects the strongly deformed character of this nucleus.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0081v1-figure-e2db5809-89d7-4d61-8988-11e36f4437ed.jpg)

MeV. The distribution is less skewed, but nevertheless it retains some trace of a prolate character. The HFB excitation energy at this temperature is about 25 MeV, much higher than energies of interest for spectroscopy and for the neutron-capture reaction. The top panel shows the distribution at T = 4 MeV. At this high excitation the distribution is featureless and close to a Gaussian.

We have also calculated P β ( q ) for 148 Sm, which is spherical in its HFB ground state. They are more symmetric and change less with temperature, consistent with the absence of a coherent quadrupole moment.

Invariants .- Fig. 4 shows the second-order invariant 〈 ˆ Q · ˆ Q 〉 vs temperature T for 148 Sm and 154 Sm. The AFMC results (circles) are compared with the HFB results (dashed lines) of Eq. (6). In HFB, 〈 ˆ Q Q · ˆ 〉 for 148 Sm can be entirely attributed to the fluctuation terms in (6). There is a small kink at T = 0 4 MeV associated with . the pairing transition, but by and large the curve is flat. The same is true of the AFMC curve. In contrast, 〈 ˆ Q Q · ˆ 〉 in 154 Sm is very different at low temperatures. In HFB, the intrinsic quadrupole moment Q 0 is large, and it persists up to a temperature of the order of 1 MeV, close to the spherical-to-deformed phase-transition temperature. The AFMC results are in semiquantitative agreement at the lowest temperatures showing that the coherent intrinsic quadrupole moment is not an artifact of the HFB. The sharp kink characterizing the HFB shape transition [19, 20] is washed out, as is expected in a finite-size system. Nevertheless a signature of this phase transition remains in the rapid decrease of 〈 Q · Q 〉 with temperature. In AFMC deformation effects survive well above

FIG. 4: 〈 Q Q · 〉 vs temperature T for the spherical 148 Sm (left) and the deformed 154 Sm (right). The AFMC results (solid circles) are compared with the HFB results (dashed lines).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0081v1-figure-9c9f55e2-aeb4-4313-a241-f6a0c34f28c4.jpg)

the transition temperature, in that 〈 Q Q · 〉 continues to be enhanced beyond its uncorrelated mean-field value.

The secondand third-order invariants can be used to define effective values of the intrinsic shape parameters β, γ [27] of the collective Bohr model [28, Sec. 6B-1a]. The model assumes an intrinsic frame in which the quadrupole deformation parameters α 2 µ = √ 5 π Q 〈 ˆ 2 µ 〉 / 3 r 2 0 A 5 / 3 are expressed as α 20 = β cos γ , α 22 = α 2 -2 = 1 √ 2 β sin γ , and α 2 ± 1 = 0. Effective β and γ can then be determined from the corresponding invariants

$$\text{-} \text{$\text{ can use common from corresponding invations}$} \\ \text{much} \\ \text{id for} \quad \beta = \frac { \sqrt { 5 \pi } } { 3 r _ { 0 } ^ { 2 } A ^ { 5 / 3 } } \langle \hat { Q } \cdot \hat { Q } \rangle ^ { 1 / 2 } \, ; \, \cos 3 \gamma = - \sqrt { \frac { 7 } { 2 } } \frac { \langle ( \hat { Q } \times \hat { Q } ) \cdot \hat { Q } \rangle } { \langle \hat { Q } \cdot \hat { Q } \rangle ^ { 3 / 2 } } \,. \\ \text{in the} \quad \text{$\text{ in addition we can extract a measure } \Lambda \, \hat { Q } \, of the flatten}.$$

In addition, we can extract a measure ∆ β of the fluctuations in β using the second- and fourth-order invariants

$$\mathbf n _ { \mathbf t } \quad & ( \Delta \beta / \beta ) ^ { 2 } = \left [ \langle ( \hat { Q } \cdot \hat { Q } ) ^ { 2 } \rangle - \langle \hat { Q } \cdot \hat { Q } \rangle ^ { 2 } \right ] ^ { 1 / 2 } / \langle \hat { Q } \cdot \hat { Q } \rangle \,. \quad ( 8 ) \\ \mathbf n t \quad & \text{The invariants themselves are calculated from the $\mathbf moc$}$$

The invariants themselves are calculated from the moments of P β ( q ) using the relations in Table I. As expected, the deformed 154 Sm has a larger deformation β than 148 Sm (0 232 vs 0 137), but a smaller deforma-. . tion angle γ (13 4 . ◦ vs 21 6 ) that is . ◦ closer to an axial shape. The deformed nucleus is more rigid in that it has a smaller ∆ β/β , 0 51 for . 154 Sm vs 0 72 for . 148 Sm.

Summary .- We have demonstrated that the distribution of the axial quadrupole operator can be computed in the AFMC method, and that it conveys important information about deformation and the intrinsic shapes of nuclei at finite temperature. In particular, the expectation values of β 2 , β 3 cos 3 γ and the fluctuation in β 2 can be extracted as a function of temperature. With these moments, it should be possible to construct models of the joint level density distribution ρ β, E ( x ) = ρ E ( x ) P E x ( β ), where ρ E ( x ) is the total level density and P E x ( β ) is the intrinsic shape distribution at excitation energy E x . This joint distribution is an important component in the theory of fission and will be discussed in a future publication.

Acknowledgments .-Y. Alhassid acknowledges the hospitality of the Institute of Nuclear Theory in Seattle, where part of this work was completed. We thank H. Nakada for the HFB code. This work was supported in part by the U.S. DOE grant Nos. DE-FG02-91ER40608 and DE-FG02-00ER411132. Computational cycles were provided by the NERSC high performance computing facility and by the facilities of the Yale University Faculty of Arts and Sciences High Performance Computing Center.

- [1] The SU (3) representation provides an alternate way to study deformation effects [2], but it has not been widely applied to heavy nuclei.
- [2] C.E. Vargas, J.G. Hirsch, and J.P. Draayer, Nucl. Phys. A 697 , 655 (2002).
- [3] The RIPL3 database containing tables of experimental level densities can be accessed at https://www-nds.iaea.org/RIPL-3/resonances/ .
- [4] J.C. Pei, W. Nazarewicz, J.A. Sheikh, and A.K. Kerman, Phys. Rev. Lett. 102 , 192501 (2009).
- [5] Y. Alhassid, G.F. Bertsch, D.J. Dean, and S.E. Koonin, Phys. Rev. Lett. 77 , 1444 (1996).
- [6] Y. Alhassid, S. Liu, and H. Nakada, Phys. Rev. Lett. 99 , 162504 (2007).
- [7] K. Kumar, Phys. Rev. Lett. 28 , 249 (1972).
- [8] D. Cline, Ann. Rev. Nucl. Part. Sci. 36 , 683 (1986).
- [9] Strictly speaking the uniqueness of the fourth- and fifthorder invariants holds when the operators ˆ Q 2 µ commute among themselves. The quadrupole operators commute in coordinate space but generally not in a truncated shell model space. However, the effects due to their noncommutation is small.
- [10] J. Hubbard, Phys. Rev. Lett. 3 , 77 (1959); R. L. Stratonovich, Dokl. Akad. Nauk SSSR [Sov. Phys. -
- Dokl.] 115 , 1097 (1957).
- [11] Small bad-sign interaction terms can be treated using the extrapolation method of Ref. [12].
- [12] Y. Alhassid, D.J. Dean, S.E. Koonin, G. Lang, and W.E. Ormand, Phys. Rev. Lett. 72 , 613 (1994).
- [13] G.H. Lang, C.W. Johnson, S.E. Koonin, and W.E. Ormand, Phys. Rev. C 48 , 1518 (1993).
- [14] S.E. Koonin, D.J. Dean, and K. Langanke, Phys. Rep. 278 , 2 (1997); Y. Alhassid, Int. J. Mod. Phys. B 15 , 1447 (2001).
- [15] H. Nakada and Y. Alhassid, Phys. Rev. Lett. 79 , 2939 (1997); Y. Alhassid, S. Liu, and H. Nakada, Phys. Rev. Lett. 83 , 4265 (1999).
- [16] The oscillator length parameter b is defined to give a ground-state wave function of the form ∼ exp( -r 2 / b 2 2 ).
- [17] B.A. Brown and W.A. Richter, Phys. Rev. C 74 , 034315 (2006).
- [18] B. Lauritzen and G.F. Bertsch, Phys. Rev. C 39 , 2412 (1989).
- [19] A.L. Goodman, Phys. Rev. C 33 , 2212 (1986).
- [20] V. Martin, J. L. Egido and L. M. Robledo, Phys. Rev. C 68 , 034327 (2003).
- [21] B. K. Agrawal, T. Sil, J. N. De, and S. K. Samaddar, Phys. Rev. C 62 , 044307 (2000).
- [22] Y. Alhassid, S. Levit, and J. Zingman, Phys. Rev. Lett. 57 , 539 (1986).
- [23] J.-P. Delaroche, et al., Phys. Rev. C 81 , 014303 (2010).
- [24] Y. Alhassid, L. Fang and H. Nakada, Phys. Rev. Lett. 101 , 082501 (2008).
- [25] C. Ozen, Y. Alhassid, and H. Nakada, Phys. Rev. Lett. ¨ 110 , 042502 (2013).
- [26] S. Raman, C.W. Nestor, and P. Tikkanen, At. Data Nucl. Data Tables 78 , 1 (2001).
- [27] We use the symbol β for both the inverse temperature and the dimensionless deformation parameter; the meaning should be clear from the context.
- [28] A. Bohr and B. Mottelson, Nuclear Structure , Vol. II (Benjamin, 1975).
Y. Alhassid C. N. Gilbreth G. F. Bertsch
Nuclear deformation at finite temperature
Deformation, a key concept in our understanding of heavy nuclei, is based on
a mean-field description that breaks the rotational invariance of the nuclear
many-body Hamiltonian. We present a method to analyze nuclear deformations at
finite temperature in a framework that preserves rotational invariance. The
auxiliary-field Monte-Carlo method is used to generate the statistical ensemble
and calculate the probability distribution associated with the quadrupole
operator. Applying the technique to nuclei in the rare-earth region, we
identify model-independent signatures of deformation and find that deformation
effects persist to higher temperatures than the spherical-to-deformed shape
phase-transition temperature of mean-field theory.