Linear response strength functions with iterative Arnoldi diagonalization

We report on an implementation of a new method to calculate RPA strength functions with iterative non-hermitian Arnoldi diagonalization method, which does not explicitly calculate and store the RPA matrix. We discuss the treatment of spurious modes, numerical stability, and how the method scales as the used model space is enlarged. We perform the particle-hole RPA benchmark calculations for double magic nucleus 132Sn and compare the resulting electromagnetic strength functions against those obtained within the standard RPA.Comment: 9 RevTeX pages, 11 figures, submitted to Physical Review

Continuity equation and local gauge invariance for the N3LO nuclear Energy Density Functionals

Background: The next-to-next-to-next-to-leading order (N3LO) nuclear energy density functional extends the standard Skyrme functional with new terms depending on higher-order derivatives of densities, introduced to gain better precision in the nuclear many-body calculations. A thorough study of the transformation properties of the functional with respect to different symmetries is required, as a step preliminary to the adjustment of the coupling constants. Purpose: Determine to which extent the presence of higher-order derivatives in the functional can be compatible with the continuity equation. In particular, to study the relations between the validity of the continuity equation and invariance of the functional under gauge transformations. Methods: Derive conditions for the validity of the continuity equation in the framework of time-dependent density functional theory. The conditions apply separately to the four spin-isospin channels of the one-body density matrix. Results: We obtained four sets of constraints on the coupling constants of the N3LO energy density functional that guarantee the validity of the continuity equation in all spin-isospin channels. In particular, for the scalar-isoscalar channel, the constraints are the same as those resulting from imposing the standard U(1) local-gauge-invariance conditions. Conclusions: Validity of the continuity equation in the four spin-isospin channels is equivalent to the local-gauge invariance of the energy density functional. For vector and isovector channels, such validity requires the invariance of the functional under local rotations in the spin and isospin spaces.Comment: 12 Latex pages, submitted to Physical Review

A Substorm Onset Observed by the POLAR Spacecraft in Conjunction with the IMAGE Chain

We present observations of the Polar spacecraft of magnetospheric substorm signatures in the plasma sheet midway along auroral field lines between the ionosphere and the equatorial plasma sheet. An isolated substorm on October 17, 1997, is studied in detail. The onset time at 2040:42 UT was defined to be almost simultaneous (within 20 sec) on the ground and at Polar altitude (~ 4RE. At Polar, the onset was manifested as strong magnetic field variations and plasma flow bursts with amplitudes exceeding 20 nT and 100 km/s, respectively. Bursts of parallel Poynting flux of ~ 0.5 ergs/cm2 were related to these variations, and they were predominantly directed toward the ionosphere. In addition, a sequence of weak magnetic field variations and flow bursts were observes at Polar ~ 10 min before the onset. Associated with this, a weak westward electrojet was observed on the ground. We discuss the importance of these observations in the physics of substorm onset

Neutron-Proton Correlations in an Exactly Solvable Model

We examine isovector and isoscalar neutron-proton correlations in an exactly solvable model based on the algebra SO(8). We look particularly closely at Gamow-Teller strength and double beta decay, both to isolate the effects of the two kinds of pairing and to test two approximation schemes: the renormalized neutron-proton QRPA (RQRPA) and generalized BCS theory. When isoscalar pairing correlations become strong enough a phase transition occurs and the dependence of the Gamow-Teller beta+ strength on isospin changes in a dramatic and unfamiliar way, actually increasing as neutrons are added to an N=Z core. Renormalization eliminates the well-known instabilities that plague the QRPA as the phase transition is approached, but only by unnaturally suppressing the isoscalar correlations. Generalized BCS theory, on the other hand, reproduces the Gamow-Teller strength more accurately in the isoscalar phase than in the usual isovector phase, even though its predictions for energies are equally good everywhere. It also mixes T=0 and T=1 pairing, but only on the isoscalar side of the phase transition.Comment: 13 pages + 11 postscript figures, in RevTe

Fully-Renormalized QRPA fulfills Ikeda sum rule exactly

The renormalized quasiparticle-RPA is reformulated for even-even nuclei using restrictions imposed by the commutativity of the phonon creation operator with the total particle number operator. This new version, Fully-Renormalized QRPA (FR-QRPA), is free from the spurious low-energy solutions. Analytical proof is given that the Ikeda sum rule is fullfiled within the FR-QRPA.Comment: 9 page

Extension of random-phase approximation preserving energy weighted sum rules: an application to a 3-level Lipkin model

A limitation common to all extensions of random-phase approximation including only particle-hole configurations is that they violate to some extent the energy weighted sum rules. Considering one such extension, the improved RPA (IRPA), already used to study the electronic properties of metallic clusters, we show how it can be generalized in order to eliminate this drawback. This is achieved by enlarging the configuration space, including also elementary excitations corresponding to the annihilation of a particle (hole) and the creation of another particle (hole) on the correlated ground state. The approach is tested within a solvable 3-level model.Comment: 2 figure

Reconciliation of the Substorm Onset Determined on the Ground and at the Polar spacecraft

An isolated substorm on Oct. 17, 1997 during a close conjunction of the Polar spacecraft and the ground-based MIRACLE network is studied in detail. We identify signatures of substorm onset in the plasma sheet midway between the ionosphere and the equatorial plasma sheet, determine their timing relative to the ground signatures, and discuss their counterparts on the ground and in the equatorial plasma sheet. The substorm onset is determined as the negative bay onset at 2040:42(≠ 5 sec) UT coinciding with the onset of auroral precipitation, energization of plasma sheet electrons at Polar, and strong magnetic field variations perpendicular to the ambient field. Such accurate timing coincidence is consistent with the Alfvén transit time between Polar and the ionosphere. Furthermore, the timing of other field and particle signatures at Polar showed clear deviations from the onset time (≠ 2 min). This suggests that the sequence of these signatures around the onset time can be used to validate the signatures predicted by various substorm onset models

A large Hilbert space QRPA and RQRPA calculation of neutrinoless double beta decay

A large Hilbert space is used for the calculation of the nuclear matrix elements governing the light neutrino mass mediated mode of neutrinoless double beta decay of Ge76, Mo100, Cd116, Te128 and Xe136 within the proton-neutron quasiparticle random phase approximation (pn-QRPA) and the renormalized QRPA with proton-neutron pairing (full-RQRPA) methods. We have found that the nuclear matrix elements obtained with the standard pn-QRPA for several nuclear transitions are extremely sensitive to the renormalization of the particle-particle component of the residual interaction of the nuclear hamiltonian. Therefore the standard pn-QRPA does not guarantee the necessary accuracy to allow us to extract a reliable limit on the effective neutrino mass. This behaviour, already known from the calculation of the two-neutrino double beta decay matrix elements, manifests itself in the neutrinoless double-beta decay but only if a large model space is used. The full-RQRPA, which takes into account proton-neutron pairing and considers the Pauli principle in an approximate way, offers a stable solution in the physically acceptable region of the particle-particle strength. In this way more accurate values on the effective neutrino mass have been deduced from the experimental lower limits of the half-lifes of neutrinoless double beta decay.Comment: 19 pages, RevTex, 1 Postscript figur

Two neutrino positron double beta decay of 106^{106}Cd for 0+→0+0^+ \to 0^+ transition

The two neutrino positron double beta decay of $^{106}$Cd for 0$^{+} \to$ 0$^{+}$ transition has been studied in the Hartree-Fock-Bogoliubov model in conjunction with the summation method. In the first step, the reliability of the intrinsic wave functions of $^{106}$Cd and $^{106}$Pd nuclei has been tested by comparing the theoretically calculated results for yrast spectra, reduced $B(E2$:$0^{+}\to 2^{+})$ transition probabilities, quadrupole moments $Q(2^{+})$ and gyromagnetic factors $g(2^{+})$ with the available experimental data. In the second step, the nuclear transition matrix element $M_{2\nu}$ and the half-life $T_{1/2}^{2\nu}$ for 0$^{+}\to$ 0$^{+}$ transition have been calculated with these wave functions. Moreover, we have studied the effect of deformation on nuclear transition matrix element M$_{2\nu}.$Comment: 20 page

Single- and double-beta decay Fermi-transitions in an exactly solvable model

An exactly solvable model suitable for the description of single and double-beta decay processes of the Fermi-type is introduced. The model is equivalent to the exact shell-model treatment of protons and neutrons in a single j-shell. Exact eigenvalues and eigenvectors are compared to those corresponding to the hamiltonian in the quasiparticle basis (qp) and with the results of both the standard quasiparticle random phase approximation (QRPA) and the renormalized one (RQRPA). The role of the scattering term of the quasiparticle hamiltonian is analyzed. The presence of an exact eigenstate with zero energy is shown to be related to the collapse of the QRPA. The RQRPA and the qp solutions do not include this zero-energy eigenvalue in their spectra, probably due to spurious correlations. The meaning of this result in terms of symmetries is presented.Comment: 29 pages, 9 figures included in a Postsript file. Submitted to Physcal Review