Density Matrix Renormalization Group and the Nuclear Shell Model

We describe the use of the Density Matrix Renormalization Group method as a means of approximately solving large-scale nuclear shell-model problems. We focus on an angular-momentum-conserving variant of the method and report test results for the nucleus $^{48}Cr$. The calculation is able to reproduce both the ground state energy and the energy of the first excited state, by diagonalizing matrices much smaller than those of the full shell model.Comment: 7 pages, 3 figures; To appears in Phys. Rev.

"Almost stable" matchings in the Roommates problem

An instance of the classical Stable Roommates problem (SR) need not admit a stable matching. This motivates the problem of finding a matching that is “as stable as possible”, i.e. admits the fewest number of blocking pairs. In this paper we prove that, given an SR instance with n agents, in which all preference lists are complete, the problem of finding a matching with the fewest number of blocking pairs is NP-hard and not approximable within n^{\frac{1}{2}-\varepsilon}, for any \varepsilon>0, unless P=NP. If the preference lists contain ties, we improve this result to n^{1-\varepsilon}. Also, we show that, given an integer K and an SR instance I in which all preference lists are complete, the problem of deciding whether I admits a matching with exactly K blocking pairs is NP-complete. By contrast, if K is constant, we give a polynomial-time algorithm that finds a matching with at most (or exactly) K blocking pairs, or reports that no such matching exists. Finally, we give upper and lower bounds for the minimum number of blocking pairs over all matchings in terms of some properties of a stable partition, given an SR instance I

A New Approach to Large-Scale Nuclear Structure Calculations

A new approach to large-scale nuclear structure calculations, based on the Density Matrix Renormalization Group (DMRG), is described. The method is tested in the context of a problem involving many identical nucleons constrained to move in a single large-j shell and interacting via a pairing plus quadrupole interaction. In cases in which exact diagonalization of the hamiltonian is possible, the method is able to reproduce the exact results for the ground state energy and the energies of low-lying excited states with extreme precision. Results are also presented for a model problem in which exact solution is not feasible.Comment: 6 pages + 1 eps figur

Dobinski-type relations and the Log-normal distribution

We consider sequences of generalized Bell numbers B(n), n=0,1,... for which there exist Dobinski-type summation formulas; that is, where B(n) is represented as an infinite sum over k of terms P(k)^n/D(k). These include the standard Bell numbers and their generalizations appearing in the normal ordering of powers of boson monomials, as well as variants of the "ordered" Bell numbers. For any such B we demonstrate that every positive integral power of B(m(n)), where m(n) is a quadratic function of n with positive integral coefficients, is the n-th moment of a positive function on the positive real axis, given by a weighted infinite sum of log-normal distributions.Comment: 7 pages, 2 Figure

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

New Discrete Basis for Nuclear Structure Studies

A complete discrete set of spherical single-particle wave functions for studies of weakly-bound many-body systems is proposed. The new basis is obtained by means of a local-scale point transformation of the spherical harmonic oscillator wave functions. Unlike the harmonic oscillator states, the new wave functions decay exponentially at large distances. Using the new basis, characteristics of weakly-bound orbitals are analyzed and the ground state properties of some spherical doubly-magic nuclei are studied. The basis of the transformed harmonic oscillator is a significant improvement over the harmonic oscillator basis, especially in studies of exotic nuclei where the coupling to the particle continuum is important.Comment: 13 pages, RevTex, 6 p.s. figures, submitted to Phys. Rev.

Polynomial iterative algorithms for coloring and analyzing random graphs

We study the graph coloring problem over random graphs of finite average connectivity $c$. Given a number $q$ of available colors, we find that graphs with low connectivity admit almost always a proper coloring whereas graphs with high connectivity are uncolorable. Depending on $q$, we find the precise value of the critical average connectivity $c_q$. Moreover, we show that below $c_q$ there exist a clustering phase $c\in [c_d,c_q]$ in which ground states spontaneously divide into an exponential number of clusters. Furthermore, we extended our considerations to the case of single instances showing consistent results. This lead us to propose a new algorithm able to color in polynomial time random graphs in the hard but colorable region, i.e when $c\in [c_d,c_q]$.Comment: 23 pages, 10 eps figure

The Density Matrix Renormalization Group for finite Fermi systems

The Density Matrix Renormalization Group (DMRG) was introduced by Steven White in 1992 as a method for accurately describing the properties of one-dimensional quantum lattices. The method, as originally introduced, was based on the iterative inclusion of sites on a real-space lattice. Based on its enormous success in that domain, it was subsequently proposed that the DMRG could be modified for use on finite Fermi systems, through the replacement of real-space lattice sites by an appropriately ordered set of single-particle levels. Since then, there has been an enormous amount of work on the subject, ranging from efforts to clarify the optimal means of implementing the algorithm to extensive applications in a variety of fields. In this article, we review these recent developments. Following a description of the real-space DMRG method, we discuss the key steps that were undertaken to modify it for use on finite Fermi systems and then describe its applications to Quantum Chemistry, ultrasmall superconducting grains, finite nuclei and two-dimensional electron systems. We also describe a recent development which permits symmetries to be taken into account consistently throughout the DMRG algorithm. We close with an outlook for future applications of the method.Comment: 48 pages, 17 figures Corrections made to equation 19 and table

Boson mappings and four-particle correlations in algebraic neutron-proton pairing models

Neutron-proton pairing correlations are studied within the context of two solvable models, one based on the algebra SO(5) and the other on the algebra SO(8). Boson-mapping techniques are applied to these models and shown to provide a convenient methodological tool both for solving such problems and for gaining useful insight into general features of pairing. We first focus on the SO(5) model, which involves generalized T=1 pairing. Neither boson mean-field methods nor fermion-pair approximations are able to describe in detail neutron-proton pairing in this model. The analysis suggests, however, that the boson Hamiltonian obtained from a mapping of the fermion Hamiltonian contains a pairing force between bosons, pointing to the importance of boson-boson (or equivalently four-fermion) correlations with isospin T=0 and spin S=0. These correlations are investigated by carrying out a second boson mapping. Closed forms for the fermion wave functions are given in terms of the fermion-pair operators. Similar techniques are applied -- albeit in less detail -- to the SO(8) model, involving a competition between T=1 and T=0 pairing. Conclusions similar to those of the SO(5) analysis are reached regarding the importance of four-particle correlations in systems involving neutron-proton pairing.Comment: 31 pages, Latex, 3 Postscript figures, uses epsf.sty, submitted to Physical Review