BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
G. F. Bertsch appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher 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.
The codes HFgradZ.py and HFgradT.py find axially symmetric minima of a Hartree-Fock energy functional for a Hamiltonian supplied in a shell-model basis. The functional to be minimized is the Hartree-Fock energy for zero-temperature properties or the Hartree-Fock grand potential for finite-temperature properties. Various constraints may be imposed on the minima.
The configuration-interaction shell model approach provides an attractive framework for the calculation of nuclear level densities in the presence of correlations, but the large dimensionality of the model space has hindered its application in mid-mass and heavy nuclei. The shell model Monte Carlo (SMMC) method permits calculations in model spaces that are many orders of magnitude larger than spaces that can be treated by conventional diagonalization methods. We discuss recent progress in the SMMC approach to level densities, and in particular the calculation of level densities in heavy nuclei. We calculate the distribution of the axial quadrupole operator in the laboratory frame at finite temperature and demonstrate that it is a model-independent signature of deformation in the rotational invariant framework of the shell model. We propose a method to use these distributions for calculating level densities as a function of intrinsic deformation.
We proposed a calculational framework for describing induced fission that avoids the Bohr-Wheeler assumption of well-defined fission channels. The building blocks of our approach are configurations that form a discrete, orthogonal basis and can be characterized by both energy and shape. The dynamics is to be determined by interaction matrix elements between the states rather than by a Hill-Wheeler construction of a collective coordinate. Within our approach, several simple limits can be seen: diffusion; quantized conductance; and ordinary decay through channels. The specific proposal for the discrete basis is to use the $K^π$ quantum numbers of the axially symmetric Hartree-Fock approximation to generate the configurations. Fission paths would be determined by hopping from configuration to configuration via the residual interaction. We show as an example the configurations needed to describe a fictitious fission decay $^{32}{\rm S} \rightarrow ^{16}{\rm O} + ^{16}{\rm O}$. We also examine the geometry of the path for fission of $^{236}$U, measuring distances by the number of jumps needed to go to a new $K^π$ partition.
We assess the accuracy of finite-temperature mean-field theory using as a standard the Hamiltonian and model space of the shell model Monte Carlo calculations. Two examples are considered: the nucleus $^{162}$Dy, representing a heavy deformed nucleus, and $^{148}$Sm, representing a nearby heavy spherical nucleus with strong pairing correlations. The errors inherent in the finite-temperature Hartree-Fock and Hartree-Fock-Bogoliubov approximations are analyzed by comparing the entropies of the grand canonical and canonical ensembles, as well as the level density at the neutron resonance threshold, with shell model Monte Carlo (SMMC) calculations, which are accurate up to well-controlled statistical errors. The main weak points in the mean-field treatments are seen to be: (i) the extraction of number-projected densities from the grand canonical ensembles, and (ii) the symmetry breaking by deformation or by the pairing condensate. In the absence of a pairing condensate, we confirm that the usual saddle-point approximation to extract the number-projected densities is not a significant source of error compared to other errors inherent to the mean-field theory. We also present an alternative formulation of the saddle-point approximation that makes direct use of an approximate particle-number projection and avoids computing the usual three-dimensional Jacobian of the saddle-point integration. We find that the pairing condensate is less amenable to approximate particle-number projection methods due to the explicit violation of particle-number conservation in the pairing condensate. Nevertheless, the Hartree-Fock-Bogoliubov theory is accurate to less than one unit of entropy for $^{148}$Sm at the neutron threshold energy, which is above the pairing phase transition.
We suggest a small set of fission observables to be used as test cases for validation of theoretical calculations. The purpose is to provide common data to facilitate the comparison of different fission theories and models. The proposed observables are chosen from fission barriers, spontaneous fission lifetimes, fission yield characteristics, and fission isomer excitation energies.
Levels densities of independent-particle Hamiltonians can be calculated easily by using the real-time representation of the evolution operator together with the fast Fourier transform. We describe the method and implement it with a set of Python programs. Examples are provided for the total and partial levels densities of a heavy deformed nucleus (Dy-164). The partial level densities that may be calculated are the projected ones on neutron number, proton number, azimuthal angular momentum, and parity.
Deformation, a key concept in our understanding of heavy nuclei, is based on a mean-field description that breaks the rotational invariance of the nuclear many-body Hamiltonian. We present a method to analyze nuclear deformations at finite temperature in a framework that preserves rotational invariance. The auxiliary-field Monte-Carlo method is used to generate the statistical ensemble and calculate the probability distribution associated with the quadrupole operator. Applying the technique to nuclei in the rare-earth region, we identify model-independent signatures of deformation and find that deformation effects persist to higher temperatures than the spherical-to-deformed shape phase-transition temperature of mean-field theory.
For 75 years the theory of nuclear fission has been based on the existence of a collective coordinate associated with the nuclear shape, an assumption required by the Bohr-Wheeler formula as well as by the R-matrix theory of fission. We show that it is also possible to formulate the theory without the help of collective coordinates. In the new formulation, fission is facilitated by individual states in the barrier region rather than channels over the barrier. In a certain limit the theory reduces to a formula closely related to the formula for electronic conductance through resonant tunneling states. In contrast, conduction through channels gives rise to a staircase excitation function that is well-known in nanoscale electronics but has never been seen in nuclear fission.
We calculate the energy deposition by very short laser pulses in SiO_2 (alpha-quartz) with a view to establishing systematics for predicting damage and nanoparticle production. The theoretical framework is time-dependent density functional theory, implemented by the real-time method in a multiscale representation. For the most realistic simulations we employ a meta-GGA Kohn-Sham potential similar to that of Becke and Johnson. We find that the deposited energy in the medium can be accurately modeled as a function of the local electromagnetic pulse fluence. The energy-deposition function can in turn be quite well fitted to the strong-field Keldysh formula for a range of intensities from below the melting threshold to well beyond the ablation threshold. We find reasonable agreement between the damage threshold and the energy required to melt the substrate. The ablation threshold estimated by the energy to convert the substrate to an atomic fluid is higher than the measurement, indicating significance of nonthermal nature of the process. A fair agreement is found for the depth of the ablation.
We calculate the dielectric response of crystalline silicon following irradiation by a high-intensity laser pulse, modeling the dynamics by time-dependent density functional theory (TDDFT). The pump-probe measurements are numerically simulated by solving the time-dependent Kohn-Sham equation with the pump and probe fields included as external fields. As expected, the excited silicon shows features of a particle-hole plasma in its response. We compare the calculated response with a thermal model and with a simple Drude model. The thermal model requires only a static DFT calculation to prepare electronically excited matter and agrees rather well with the TDDFT for the same particle-hole density. The Drude model with two fitted parameters (electron effective mass and collision time) also shows fair agreement at the lower excitation energies; the fitted effective masses are consistent with carrier-band dispersions. The extracted Drude lifetimes range from 6 fs at weak pumping fields to much lower values at high fields. However, we find that the Drude model does not give a good fit to the imaginary dielectric function at the highest fields. Comparing the thermal model with the Drude, we find that the extracted lifetimes are in the same range, 1-13 fs depending on the temperature. These short Drude lifetimes show that strong damping is possible in the TDDFT, despite the absence of electron scattering. One significant difference between the TDDFT response and the other models is that the response to the probe pulse depends on the polarization of the pump pulse. We also find that the imaginary part of the dielectric function can be negative, particularly for the parallel polarization of pump and probe fields.
The Gallagher-Moszkowski rule in the spectroscopy of odd-odd nuclei imposes a new spin constraint on the energy functionals for self-consistent mean field theory. The commonly used parameterization of the effective three-body interaction in the Gogny and Skyrme families of energy functionals is ill-suited to satisfy the spin constraint. In particular, the Gogny parameterization of the three-body interaction has the opposite spin dependence to that required by the observed spectra. The two-body part has a correct sign, but in combination the rule is violated as often as not. We conclude that a new functional form is needed for the effective three-body interaction that can take into better account the different spin-isospin channels of the interaction.
We derive a time-dependent density functional theory appropriate for calculating the near-edge X-ray absorption spectrum in molecules and condensed matter. The basic assumption is to increase the space of many-body wave functions from one Slater determinant to two. The equations of motion derived from Dirac's variational principle provide an exact solution for the linear response when the interaction Hamiltonian has only a core-electron field. The equations can be solved numerically nearly as easily as the ordinary real-time time-dependent Kohn-Sham equations. We carry out the solution under conditions that permit comparison with the expected power-law behavior. Our extracted power-law exponents are similar to those derived by Nozieres and DeDominicis, but are not in quantitative agreement. We argue that our calculational method can be readily generalized to density functionals that take into account the more general electron-electron interactions that are needed for treating dynamic effects such as plasmon excitations.
Spectroscopic observables such as electromagnetic transitions strengths can be related to the properties of the intrinsic mean-field wave function when the latter are strongly deformed, but the standard rotational formulas break down when the deformation decreases. Nevertheless there is a well-defined, non-zero, spherical limit that can be evaluated in terms of overlaps of mean-field intrinsic deformed wave functions. We examine the transition between the spherical limit and strongly deformed one for a range of nuclei comparing the two limiting formulas with exact projection results. We find a simple criterion for the validity of the rotational formula depending on $<Δ\vec{J}^2>$, the mean square fluctuation in the angular momentum of the intrinsic state. We also propose an interpolation formula which describes the transition strengths over the entire range of deformations, reducing to the two simple expressions in the appropriate limits.
The time-dependent density functional theory (TDDFT) is the leading computationally feasible theory to treat excitations by strong electromagnetic fields. Here the theory is applied to coherent optical phonon generation produced by intense laser pulses. We examine the process in the crystalline semimetal antimony (Sb), where nonadiabatic coupling is very important. This material is of particular interest because it exhibits strong phonon coupling and optical phonons of different symmetries can be observed. The TDDFT is able to account for a number of qualitative features of the observed coherent phonons, despite its unsatisfactory performance on reproducing the observed dielectric functions of Sb. A simple dielectric model for nonadiabatic coherent phonon generation is also examined and compared with the TDDFT calculations.
We analyze the effect of pairing on particle transport in time-dependent theories based on the Hartree-Fock-Bogoliubov (HFB) or BCS approximations. The equations of motion for the HFB density matrices are unique and the theory respects the usual conservation laws defined by commutators of the conserved quantity with the Hamiltonian. In contrast, the theories based on the BCS approximation are more problematic. In the usual formulation of TDHF+BCS, the equation of continuity is violated and one sees unphysical oscillations in particle densities. This can be ameliorated by freezing the occupation numbers during the evolution in TDHF+BCS, but there are other problems with the BCS that make it doubtful for reaction dynamics. We also compare different numerical implementations of the time-dependent HFB equations. The equations of motion for the $U$ and $V$ Bogoliubov transformations are not unique, but it appears that the usual formulation is also the most efficient. Finally, we compare the time-dependent HFB solutions with numerically exact solutions of the two-particle Schrodinger equation. Depending on the treatment of the initial state, the HFB dynamics produces a particle emission rate at short times similar to that of the Schrodinger equation. At long times, the total particle emission can be quite different, due to inherent mean-field approximation of the HFB theory.
As part of a program to study odd-A nuclei in the Hartree-Fock-Bogoliubov (HFB) theory, we have developed a new calculational tool to find the HFB minima of odd-A nuclei based on the gradient method and using interactions of Gogny's form. The HFB minimization includes both time-even and time-odd fields in the energy functional, avoiding the commonly used "filling approximation". Here we apply the method to calculate neutron pairing gaps in some representative isotope chains of spherical and deformed nuclei, namely the Z=8,50 and 82 spherical chains and the Z=62 and 92 deformed chains. We find that the gradient method is quite robust, permitting us to carry out systematic surveys involving many nuclei. We find that the time-odd field does not have large effect on the pairing gaps calculated with the Gogny D1S interaction. Typically, adding the T-odd field as a perturbation increases the pairing gap by ~100 keV, but the re-minimization brings the gap back down. This outcome is very similar to results reported for the Skyrme family of nuclear energy density functionals. Comparing the calculated gaps with the experimental ones, we find that the theoretical errors have both signs implying that the D1S interaction has a reasonable overall strength. However, we find some systematic deficiencies comparing spherical and deformed chains and comparing the lighter chains with the heavier ones. The gaps for heavy spherical nuclei are too high, while those for deformed nuclei tend to be too low. The calculated gaps of spherical nuclei show hardly any A-dependence, contrary to the data. Inclusion of the T-odd component of the interaction does not change these qualitative findings.
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 present a computational methodology for a theory of the lowest octupole excitations applicable to all even-even nuclei beyond the lightest. The theory is the well-known generator-coordinate extension (GCM) of the Hartree-Fock-Bogoliubov self-consistent mean field theory (HFB). We use the discrete-basis Hill-Wheel method (HW) to compute the wave functions with an interaction from the Gogny family of Hamiltonians. Comparing to the compiled experimental data on octupole excitations, we find that the performance of the theory depends on the deformation characteristics of the nucleus. For nondeformed nuclei, the theory reproduces the energies to about 20 % apart from an overall scale factor of about 1.6. The performance is somewhat poorer for (quadrupole) deformed nuclei, and for both together the dispersion of the scaled energies about the experimental values is about 25 %. This compares favorably with the performance of similar theories of the quadrupole excitations. Nuclei having static octupole deformations in HFB form a special category. These nuclei have the smallest measured octupole excitation energies as well as the smallest predicted energies. However, in these cases the energies are seriously underpredicted by the theory. We find that a simple two-configuration approximation, the Minimization After Projection method, (MAP) is almost as accurate as the full HW treatment, provided that the octupole-deformed nuclei are omitted from the comparison. This article is accompanied by a tabulation of the predicted octupole excitations for 818 nuclei extending from dripline to dripline, computed with several variants of the Gogny interaction.
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 long-range correlations for nuclear masses have a significant effect on the synthesis of heavy elements by the r-process. As calculated by Delaroche et al. [1], these correlations suppress magic number effects associated with minor shells. This impacts the calculated abundances before the third r-process peak (at mass number A~195), where the abundances are low and form a trough. This trough and the position of the third abundance peak are strongly affected by the masses of nuclei in the transition region between deformed and spherical. Based on different astrophysical environments, our results demonstrate that a microscopic theory of nuclear masses including correlations naturally smoothens the separation energies, thus reducing the trough and improving the agreement with observed solar system abundances.
We present a pfaffian formula for projection and symmetry restoration for wave functions of the general Bogoliubov form, including quasiparticle excited states and linear combinations of them. This solves a long-standing problem in calculating states of good symmetry, arising from the sign ambiguity of the commonly used determinant formula. A simple example is given of projecting good particle number and angular momentum from a Bogoliubov wave function in the Fock space of a single j-shell.
We apply the coupled dynamics of time-dependent density functional theory and Maxwell equations to the interaction of intense laser pulses with crystalline silicon. As a function of electromagnetic field intensity, we see several regions in the response. At the lowest intensities, the pulse is reflected and transmitted in accord with the dielectric response, and the characteristics of the energy deposition is consistent with two-photon absorption. The absorption process begins to deviate from that at laser intensities ~ 10^13 W/cm^2, where the energy deposited is of the order of 1 eV per atom. Changes in the reflectivity are seen as a function of intensity. When it passes a threshold of about 3 \times 1012 W/cm2, there is a small decrease. At higher intensities, above 2 \times 10^13 W/cm^2, the reflectivity increases strongly. This behavior can be understood qualitatively in a model treating the excited electron-hole pairs as a plasma.
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 apply the adiabatic time-dependent density functional theory to magnetic ci the real-space, real-time computational method. The standard formulas for the MCD response and its A and B terms are derived from the observables in the time-dependent wave function. We find the real time method is well suited for calculating the overall spectrum, particularly at higher excitation energies where individual excited states are numerous and overlapping. The MCD sum rules are derived and interpreted in the real-time formalism; we find that they are very useful for normalization purposes and assessing the accuracy of the theory. The method is applied to MCD spectrum of C-60 using the adiabatic energy functional from the local density approximation. The theory correctly predicts the signs of the A and B terms for the lowest allowed excitations. However, the magnitudes of the terms only show qualitative agreement with experiment.
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.
Evaluation of pfaffians arises in a number of physics applications, and for some of them a direct method is preferable to using the determinantal formula. We discuss two methods for the numerical evaluation of pfaffians. The first is tridiagonalization based on Householder transformations. The main advantage of this method is its numerical stability that makes unnecessary the implementation of a pivoting strategy. The second method considered is based on Aitken's block diagonalization formula. It yields to a kind of LU (similar to Cholesky's factorization) decomposition (under congruence) of arbitrary skew-symmetric matrices that is well suited both for the numeric and symbolic evaluations of the pfaffian. Fortran subroutines (FORTRAN 77 and 90) implementing both methods are given. We also provide simple implementations in Python and Mathematica for purpose of testing, or for exploratory studies of methods that make use of pfaffians.
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.
Fermion systems with more than two components can exhibit pairing condensates of much more complex structure than the well-known single BCS condensate of spin-up and spin-down fermions. In the framework of the exactly solvable SO(8) Richardson-Gaudin model with SU(4)-symmetric Hamiltonians, we show that the BCS approximation remains valid in the thermodynamic limit of large systems for describing the ground state energy and the canonical and quasiparticle excitation gaps. Correlations beyond BCS pairing give rise to a spectrum of collective excitations, but these do not affect the bulk energy and quasiparticle gaps.
The nuclear pairing condensate is expected to change character from spin-singlet to spin-triplet when the nucleus is very large and the neutron and proton numbers $Z,N$ are equal. We investigate the transition between these two phases within the framework of the Hartree-Fock-Bogoliubov equations, using a zero-range interaction to generate the pairing. We confirm that extremely large nucleus would indeed favor triplet pairing condensates, with the Hamiltonian parameters taken from known systematics. The favored phase is found to depend on the specific orbitals at the Fermi energy. The smallest nuclei with a well-developed spin-triplet condensate are in the mass region A ~ 130-140.
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.
We apply the configuration-interaction method to calculate the spectra of two-component Fermi systems in a harmonic trap, studying the convergence of the method at the unitary interaction limit. We find that for a fixed regularization of the two-body interaction the convergence is exponential or better in the truncation parameter of the many-body space. However, the conventional regularization is found to have poor convergence in the regularization parameter, with an error that scales as a low negative power of this parameter. We propose a new regularization of the two-body interaction that produces exponential convergence for systems of three and four particles. From the systematics, we estimate the ground-state energy of the four-particle system to be (5.05 +- 0.024)hbar omega.
We investigate the accuracy of several schemes to calculate ground-state correlation energies using the generator coordinate technique. Our test-bed for the study is the $sd$ interacting boson model, equivalent to a 6-level Lipkin-type model. We find that the simplified projection of a triaxial generator coordinate state using the $S_3$ subgroup of the rotation group is not very accurate in the parameter space of the Hamiltonian of interest. On the other hand, a full rotational projection of an axial generator coordinate state gives remarkable accuracy. We also discuss the validity of the simplified treatment using the extended Gaussian overlap approximation (top-GOA), and show that it works reasonably well when the number of boson is four or larger.