Fermionic Molecular Dynamics (FMD) models a system of fermions by means of many-body states which are composed of antisymmetrized products of single-particle states. These consist of one or several Gaussians localized in coordinate and momentum space. The parameters specifying them are the dynamical variables of the model. As the repulsive core of the nucleon-nucleon interaction induces short range correlations which cannot be accommodated by a Slater determinant, a novel approach, the unitary correlation operator method (UCOM), is applied. The unitary correlator moves two particles away from each other whenever their relative distance is within the repulsive core. The time-dependent variational principle yields the equations of motion for the variables. Energies of the stationary ground states are calculated and compared to exact many-body results for nuclei up to Ca 48. Time-dependent solutions are shown for collisions between nuclei.
In this work a unitary correlation operator is presented that explicitly describes the short-ranged central and tensor correlations in the nuclear many-body system. These short-ranged correlations are induced by the repulsive core and the pronounced tensor force of realistic nucleon-nucleon interactions and cannot be described by the simple many-body states of a mean-field or shell model approach. The unitary correlation operator is discussed for the Argonne V18 and the Bonn-A interactions. Applying the correlation operator onto the Hamiltonian a common effective interaction for low momenta is obtained. Calculations for 4He using the one- and two-body part of the correlated Hamiltonian compare favorably with exact many-body methods. Calculations for 16O and 40Ca which are not possible with exact methods are performed using harmonic oscillator shell model states. The observed deviation from the experimental binding energies and radii is attributed to the missing three-body forces. The correlated interaction in a basis-free representation is used in the Fermionic Molecular Dynamics model to calculate the nuclei of the p- and the sd-shell. In this model the uncorrelated many-particle state is given as a Slaterdeterminant of Gaussian wave-packets with spin and isospin degrees of freedom which allows a consistent description of spherically symmetric, intrinsically deformed and clustered nuclei.
Within the framework of fermionic molecular dynamics a method is developed to better account for long-range tensor correlations in nuclei when working with a single Slater determinant. Single-particle states with mixed isospin and broken parity build up an intrinsic Slater determinant, which is then charge and parity projected. By minimizing the energy of this many-body state with respect to the parameters of the single-particle states and projecting afterward on angular momentum, ground-state energies are obtained that are systematically lower than corresponding Hartree-Fock results. The realistic Argonne V18 potential is used and short-range correlations are treated with the unitary correlation operator method. Comparison with exact few-body calculations shows that in
4He about one-fifth of the correlation energy from long-range correlations are accounted for. These correlations, which extend over the whole nucleus, are visualized with the isospin and spin-isospin density of the intrinsic state. The divergence of the spin-isospin density, the source for pion fields, turns out to be of dipole nature.
A theory is formulated to derive an effective interaction in nuclei in the framework of the unitarymodel-operator approach.A unitary operator is introduced to describe two-body correlations, and it is determined from the condition that the effective interaction should be decoupled between the lowmomentum model space and the high-momentum space.The unitary-transformed Hamiltonian is represented in a cluster expansion form.The one-and two-body cluster terms are mainly discussed.The theory is applied to the calculation of the energies of the ground-and one-body states with the Reid soft core potential.A large gain of the energy has been obtained for the ground state in comparison with the result of the usual G-matrix theory.As for the one-body energies obtained, agreement with the experimental values is fairly good.Although the discrepancies with the experiments still remain at the 1O~ 15 % level, the present approach shows to have some desirable features compared with the nuclear many-body theories proposed so far. § 1. IntroductionA fundamental problem in nuclear physics is to clarify the properties of nuclei in terms of the interaction among their constituent particles.Although there have been many theories or approaches, we have not yet reached the goal of reproducing successfully the experimental data.The investigation of nuclear properties, in the usual nonrelativistic approach, proceeds in two steps: The first one is to determine the interaction among nucleons, and the second is to solve the many-body SchrOdinger equation.Concerning the nucleon-nucleon force, it seems that the phenomenological potential, such as the Reid potential, would be incomplete and one should wait until quark dynamics, probably quantum chromodynamics (QCD), solves finally the problem o~ nuclear force.It is of no doubt that in order to fit the experiments we need more accurate knowledge on the nuclear force including the threebody force.However, if we want to reach the goal, we need also a good many-body theory.Regardless of the accuracy of the nuclear force, we shall not be able to solve the problem unless we have a way of constructing an accurate solution for a many-body system.We believe that many-body theory should share the responsibility for the discrepancy between theory and experiment.It will be sure that nuclear many-body theory is still to be improved.Among many theories for nuclear many-body problem, the most widely applied method is, needless to say, the linked-diagram perturbation theory.This approach may also be referred to as the folded-diagram theory based on the G-matrix interaction especially due 'to Brandow 1l and Kuo et al. 2l An outstanding problem in this approach is the convergence of the perturbation expansion.Since Schucan and Weidenmliller 3 ) pointed out the divergence of the perturbation expansion due to the presence of intruder states, this difficulty has long been a serious problem.Present status of this problem is
The short range repulsion between nucleons is treated by a unitary correlation operator which shifts the nucleons away from each other whenever their uncorrelated positions are within the replusive core. By formulating the correlation as a transformation of the relative distance between particle pairs, general analytic expressions for the correlated wave functions and correlated operators are given. The decomposition of correlated operators into irreducible n-body operators is discussed. The one- and two-body-irreducible parts are worked out explicitly and the contribution of three-body correlations is estimated to check convergence. Ground state energies of nuclei up to mass number A=48 are calculated with a spin-isospin-dependent potential and single Slater determinants as uncorrelated states. They show that the deduced energy- and mass-number-independent correlated two-body Hamiltonian reproduces all "exact" many-body calculations surprisingly well.
There is no unique and widely accepted definition of the complexity measure (CM) of a many-fermion wave function in the presence of interactions. The simplest many-fermion wave function is a Slater determinant. In shell-model or configuration interaction (CI) and other related methods, the state is represented as a superposition of a large number of Slater determinants, which in the case of CI calculations reaches about 20 billion terms [Johnson, arXiv:1809.07869]. Although in practice this number has been used as a CM for decades, it is ill defined: it is not unique, and it depends on the particular type and the number of single-particle wave functions used to construct the Slater determinants. The canonical wave functions and/or natural orbitals [L\"owdin, Adv. Phys. 5, 1 (1956); L\"owdin and Shull, Phys. Rev. 101, 1730 (1956); Bardeen et al., Phys. Rev. 108, 1175 (1957); N. N. Bogoljubov, Il Nuovo Cimento 7, 794 (1958); Valatin, Il Nuovo Cimento 7, 843 (1958); de Gennes, Superconductivity of Metals and Alloys (CRC Press, Boca Raton, FL, 1999); Ring and Schuck, The Nuclear Many-Body Problem, 1st ed. (Springer-Verlag, Berlin, 2004)] and their corresponding occupation probabilities are intrinsic properties of any many-body wave function, irrespective of the representation, and they provide a unique solution to characterize the CM. The non-negative orbital entanglement entropy, which vanishes for a Slater determinant, provides the simplest CM, while a more complete measure of complexity is the entanglement spectrum. We illustrate these aspects in the case of a complex nonequilibrium time-dependent process, induced nuclear fission described within a real-time density functional theory framework extended to superfluid systems, which can describe simultaneously the long-range and the short-range correlations between fermions. The orbital entanglement entropy of the fissioning nucleus illustrates the localization mechanism of the many-body wave function in Fock and/or Hilbert space. The (minimal) number of Slater determinants required to represent such a complex many-body wave function with a well-defined number of particles in the case presented here is about
10500. The realistic case of the highly nonequilibrium nuclear fission process illustrated here is equivalent to a system of $23.328\ifmmode\times\else\texttimes\fi{}{10}^{9}$ interacting quantum spin-1/2 particles, a very large system for the study of quantum entanglement.