179 research outputs found
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 . 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 . Given a number 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 , we find the precise value
of the critical average connectivity . Moreover, we show that below
there exist a clustering phase 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 .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
- …