BZPEERReading, trust and collaboration for researchDiscover
Workspaces
Network
Opportunities
Account
Researcher profile
G. F. Bertsch contributes to research discovery and scholarly infrastructure.
Trust snapshot
Actions
Identity and collaboration
Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.
Log in to claimDirect collaboration
Claim this author entity first to unlock direct invitations.
Research graph
Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
With a view to applying the Generator Coordinate Method to large configuration spaces, we propose a simple approximate formula to compute diabatic many-body matrix elements without having to evaluate two-body interaction matrix elements. The method is illustrated with two analytically solvable Hamiltonians based on the harmonic oscillator.
Since its beginnings, fission theory has asumed that low-energy induced fission takes place through transition-state channels at the barrier tops. Neverthess, up to now there is no microscopic theory applicable to those conditions. We suggest that modern reaction theory is suitable for this purpose, and propose a methodology based on a configuration-interaction framework using the Generator Coordinate Method (GCM). Simple reaction-theoretic models are constructed with the Gaussian Overlap Approximation (GOA) to parameterize both the dynamics within the channels and their incoherent couplings to states outside the barrier. The physical characteristics of the channels examined here are their effective bandwidths and the quality of the coupling to compound-nucleus states as measured by the transmission factor $T$. We also investigate the spacing of GCM states with respect to their degree of overlap. We find that a rather coarse mesh provides an acceptable accuracy for estimating the bandwidths and transmission factors. The common numerical stability problem in using the GCM is avoided due to the choice of meshes and the finite bandwidths of the channels. The bandwidths of the channels are largely controlled by the zero-point energy with respect to the collective coordinate in the GCM configurations.
We propose a configuration-interaction (CI) representation to calculate induced nuclear fission with explicit inclusion of nucleon-nucleon interactions in the Hamiltonian. The framework is designed for easy modeling of schematic interactions but still permits a straightforward extension to realistic ones. As a first application, the model is applied to branching ratios between fission and capture in the decay modes of excited fissile nuclei. The ratios are compared with the Bohr-Wheeler transition-state theory to explore its domain of validity. The Bohr-Wheeler theory assumes that the rates are insensitive to the final-state scission dynamics; the insensitivity is rather easily achieved in the CI parameterizations. The CI modeling is also capable of reproducing the branching ratios of the transition-state hypothesis which is one of the key ingredients in the present-day theory of induced fission.
We investigate microscopically the tunneling dynamics in spontaneous fission of atomic nuclei. To this end, we employ a schematic solvable model with a pairing-plus-quadrupole interaction. The spontaneous decay of a system is simulated by introducing a small imaginary part to the energy of a fission doorway state. We show that the many-body Hamiltonian can be reduced to an effective 2$\times$2 Hamiltonian, from which one can derive a simple approximate formula for the decay width. We particularly investigate the applicability of the adiabatic approximation, which has often been used in the literature. With typical value of the parameters, we find that the adiabatic approximation works within a factor of around 5 for the decay width.
The usual derivations of the S and K matrices for two-particle reactions proceed through the Lippmann-Schwinger equation with formal definitions of the incoming and outgoing scattering states. Here we present an alternative derivation that is carried out completely in the Hamiltonian representation, using a discrete basis of configurations for the scattering channels as well as the quasi-bound configurations of the combined fragments. We use matrix algebra to derive an explicit expression for the K matrix in terms of the Hamiltonian of the internal states of the compound system and the coupling between the channels and the internal states. The formula for the K matrix includes explicitly a real dispersive shift matrix to the internal Hamiltonian that is easily computed in the formalism. That expression is applied to derive the usual form of the S matrix as a sum over poles in the complex energy plane. Some extensions and limitations of the discrete-basis Hamiltonian formalism are discussed in the concluding remarks and in the Appendix.
The Newtonian concept of force may be useful in some aspects of the dynamics of many-particle quantum systems such as fissioning nuclei. Following Ehrenfest's method, we show that the quantum kinetic force between parts of an extended quantum system can be described by an operator acting on the boundary between the two subsystems. The contribution to the force due to a short-ranged particle interaction can also be treated in the same way. This includes interaction effects treated in density functional theory by local functionals. The force operators are applied to several simple models to demonstrate the method.
We investigate the dynamics of spontaneous fission in a configuration-interaction (CI) approach. In that formalism the decay rate is governed by an effective interaction coupling the ground-state configuration and a fission doorway configuration, with the interaction strength determined by inverting a high-dimensioned CI Hamiltonian matrix that may have a block-tridiagonal structure. It is shown that the decay rate decreases exponentially with the number of blocks at a rate determined by the largest eigenvalue of a matrix in the block space for Hamiltonians with identical off-diagonal blocks. The theory is greatly simplified by approximations similar in spirit to the adiabatic and the least-action approximations in continuum representations. Here each block is replaced by a single matrix element. While the adiabatic reduction underestimates the coupling, a reduction based on a maximum-coupling approximation works well in a schematic CI model.
The K-matrix formalism is applied to a schematic model for nuclear fission. The purpose is to explore the dependence of observables on the assumptions made about the configuration space and nucleon interaction in the Hamiltonian of the fissile nucleus. As expected, branching ratios in induced fission are found to depend sensitively on the character of the residual interaction, whether it is pairing in form or taken from a random ensemble. On the other hand, the branching ratio is not much affected by the presence of additional configurations that do not introduce new fission paths.
A computer code is presented for solving the equations of Hartree-Fock-Bogoliubov (HFB) theory by the gradient method, motivated by the need for efficient and robust codes to calculate the configurations required by extensions of HFB such as the generator coordinate method. The code is organized with a separation between the parts that are specific to the details of the Hamiltonian and the parts that are generic to the gradient method. This permits total flexibility in choosing the symmetries to be imposed on the HFB solutions. The code solves for both even and odd particle number ground states, the choice determined by the input data stream. Application is made to the nuclei in the $sd$-shell using the USDB shell-model Hamiltonian.
We recently proposed a nonlocal form for the 3-body induced interaction that is consistent with the Fock space representation of interaction operators but leads to a fractional power dependence on the density. Here we examine the implications of the nonlocality for the excitation spectrum. In the two-component weakly interacting Fermi gas, we find that it gives an effective mass that is comparable to the one in many-body perturbation theory. Applying the interaction to nuclear matter, it predicts a large enhancement to the effective mass. Since the saturation of nuclear matter is partly due to the induced 3-body interaction, fitted functionals should treat the effective mass as a free parameter, unless the two- and three-body contributions are determined from basic theory.
We show that the Bogoliubov-de Gennes equations for nuclear ground-state wave functions support solutions in which the condensate has a mixture of spin-singlet and spin-triplet pairing. We find that such mixed-spin condensates do not occur when there are equal numbers of neutrons and protons, but only when there is an isospin imbalance. Using a phenomenological Hamiltonian, we predict that such nuclei may occur in the physical region within the proton dripline. We also solve the Bogoliubov-de Gennes equations with variable constraints on the spin-singlet and spin-triplet pairing amplitudes. For nuclei that exhibit this new pairing behavior, the resulting energy surface can be rather soft, suggesting that there may be low-lying excitations associated with the spin mixing.
We show that time-dependent density functional theory (TDDFT) is applicable to coherent optical phonon generation by intense laser pulses in solids. The two mechanisms invoked in phenomenological theories, namely impulsively stimulated Raman scattering and displacive excitation, are present in the TDDFT. Taking the example of crystalline Si, we find that the theory reproduces the phenomena observed experimentally: dependence on polarization, strong growth at the direct band gap, and the change of phase from below to above the band gap. We conclude that the TDDFT offers a predictive ab initio framework to treat coherent optical phonon generation.
Density functionals for nuclei usually include an effective 3-body interaction that depends on a fractional power of the density. Using insights from the many-body theory of the low-density two-component Fermi gas, we consider a new, nonlocal, form for the energy functional that is consistent with the Fock space representation of interaction operators. In particular, there is a unique spatially nonlocal generalization of the contact form of the interaction that preserves the density-to-the-seven-thirds dependence required by the many-body theory. We calculate the ground state energies for particles in a harmonic trap using the nonlocal induced 3-body interaction, and compare them to numerically accurate Green's Function Monte Carlo calculations. Using no free parameters, we find that a nonlocality in the space domain provides a better description of the weak-coupling regime than the local-density approximation.
We study odd-even mass staggering of nuclei, also called pairing gaps, using a Skyrme self-consistent mean-field theory and a numerically exact treatment of the pairing Hamiltonian. We find that the configuration-space Monte Carlo method proposed by Cerf and Martin offers a practical computational procedure to carry out the numerical solutions in large-dimensional model spaces. Refitting the global strength of the pairing interaction for 443 neutron pairing gaps in our number-conserving treatment, we find the correction to the pairing correlation energies and pairing gaps to have rms values of 0.6 MeV and 0.12 MeV, respectively. The exact treatment provides a significant improvement in the fit to experimental gaps, although it is partially masked by a larger rms error due to deficiencies in other aspects of the theory such as the mean-field energy functional.
We study the conductance of an interconnect between two graphene leads formed by a single-atom carbon chain. Its dependence on the chemical potential and the number of atoms in the chain is qualitatively different from that in the case of normal metal leads. Electron transport proceeds via narrow resonant states in the wire. The latter arise due to strong reflection at the junctions between the chain and the leads, which is caused by the small density of states in the leads at low energy. The energy dependence of the transmission coefficient near resonance is asymmetric and acquires a universal form at small energies. We find that in the case of leads with the zigzag edges the dispersion of the edge states has a significant effect on the device conductance.
A systematic study of low energy nuclear structure at normal deformation is carried out using the Hartree-Fock-Bogoliubov theory extended by the Generator Coordinate Method and mapped onto a 5-dimensional collective quadrupole Hamiltonian. Results obtained with the Gogny D1S interaction are presented from dripline to dripline for even-even nuclei with proton numbers Z=10 to Z=110 and neutron numbers N less than 200. The properties calculated for the ground states are their charge radii, 2-particle separation energies, correlation energies, and the intrinsic quadrupole shape parameters. For the excited spectroscopy, the observables calculated are the excitation energies and quadrupole as well as monopole transition matrix elements. We examine in this work the yrast levels up to J=6, the lowest excited 0^+ states, and the two next yrare 2^+ states. The theory is applicable to more than 90% of the nuclei which have tabulated measurements. The data set of the calculated properties of 1712 even-even nuclei, including spectroscopic properties for 1693 of them, are provided in CEA website and EPAPS repository with this article \cite{epaps}.
We survey odd-even nuclear binding energy staggering using density functional theory with several treatments of the pairing interaction including the BCS, Hartree-Fock-Bogoliubov, and the Hartree-Fock-Bogoliubov with the Lipkin-Nogami approximation. We calculate the second difference of binding energies and compare with 443 measured neutron energy differences in isotope chains and 418 measured proton energy differences in isotone chains. The particle-hole part of the energy functional is taken as the SLy4 Skyrme parametrization and the pairing part of the functional is based on a contact interaction with possible density dependence. An important feature of the data, reproduced by the theory, is the sharp gap quenching at magic numbers. With the strength of the interaction as a free parameter, the theory can reproduce the data to an rms accuracy of about 0.25 MeV. This is slightly better than a single-parameter phenomenological description but slightly poorer than the usual two-parameter phenomenological form C/A^alpha . The following conclusions can be made about the performance of common parametrization of the pairing interaction: (i) there is a weak preference for a surface-peaked neutron-neutron pairing, which might be attributable to many-body effects; (ii) a larger strength is required in the proton pairing channel than in the neutron pairing channel; (iii) pairing strengths adjusted to the well-known spherical isotope chains are too weak to give a good overall fit to the mass differences.