Search

We apply a large-scale summation of Feynman diagrams, including the class of parquet-diagrams {\em plus} important contributions outside the parquet class, for calculating effective pairing interactions and subsequently the superfluid gap in P-wave pairing in neutron matter. We employ realistic nucleon-nucleon interactions of the $v_8$ type and perform calculations up to a Fermi momentum of $1.8\,$fm$^{-1}$. We find that many-body correlations lead to a strong reduction of the spin-orbit interaction, and, therefore, to an almost complete suppression of the $^3$P$_2$ and $^3$P$_2$-$^3$F$_2$ gaps. We also find pairing in $^3$P$_0$ states; the strength of the pairing gap depends sensitively on the potential model employed. Our results for triplet pairing are relevant for assessing superfluidity in neutron star interiors, whose presence can affect the cooling of neutron stars.

We examine in detail two alternative descriptions of a system of $\alpha$ particles interacting via local interactions of different character, highlighting the fact that a faithful microscopic description of such systems demands a consistent treatment of both short- and long-range correlations. In preparation, we examine four different versions of modern microscopic many-body theory and conclude by emphasizing that these approaches, although {\it a priori} very different, actually lead to the same equations for their efficient application. The only quantity that depends on the formulation of many-body theory chosen is an {\it irreducible} interaction correction. In the language of Green's functions and Feynman diagrams, it is the set of both particle-particle and particle-hole irreducible diagrams, and in variational Jastrow-Feenberg theory it is determined by {\it multipartite correlations} and {\it elementary diagrams}. We apply these theoretical methods to the calculation of the energetics, structure, thermodynamics, and dynamics of $\alpha$ matter, as well as its condensate fraction. In dimensionless units, $\alpha$ matter appears to be remarkably similar to the much-studied $^4$He quantum fluid, its low-temperature properties now basically solved in the Jastrow-Feenberg framework. Accordingly, one can have confidence in the results of application of the same procedure to $\alpha$ matter. Even so, closer examination reveals significant differences between the physics of the two systems. Within an infinite nuclear medium, alpha matter is subject to a spinoidal instability. Extended mixtures of nucleons and alpha particles are yet to be given rigorous consideration in a corresponding theoretical framework.

We develop the parquet-diagram summation method for neutron matter interacting via potentials that include spin, tensor, and spin-orbit components. For that purpose, we derive an exact expression for the sum of all ring-diagrams in terms effective local particle-hole interactions involving the above four operators. The parquet equations are closed by deriving the spin-orbit contribution to that particle-hole interaction. We show that many-body correlations screen the bare spin-orbit potential considerably, and the corrections of that screened spin-orbit potential to the other three interaction channels are quite small. We apply our method to the calculation of the response of neutron matter to density and both longitudinal and transverse spin-dependent external fields.

The dispersion relation $\epsilon(k)$ of the elementary excitations of superfluid $^4$He has been measured at very low temperatures, from saturated vapor pressure up to solidification, using a high flux time-of-flight neutron scattering spectrometer equipped with a high spatial resolution detector (10$^5$ 'pixels'). A complete determination of $\epsilon(k)$ is achieved, from very low wave-vectors up to the end of Pitaeskii's plateau. The results compare favorably in the whole the wave-vector range with the predictions of the dynamic many-body theory (DMBT). At low wave-vectors, bridging the gap between ultrasonic data and former neutron measurements, the evolution with the pressure from anomalous to normal dispersion, as well as the peculiar wave-vector dependence of the phase and group velocities, are accurately characterized. The thermodynamic properties have been calculated analytically, developing Landau's model, using the measured dispersion curve. A good agreement is found below 0.85 K between direct heat capacity measurements and the calculated specific heat, if thermodynamically consistent power series expansions are used. The thermodynamic properties have also been calculated numerically; in this case, the results are applicable with excellent accuracy up to 1.3 K, a temperature above which the dispersion relation itself becomes temperature dependent., Comment: 35 pages, 59 figures; Supplemental Material: Text, Tables, data files

We apply parquet-diagram summation methods for the calculation of the superfluid gap in $S$-wave pairing in neutron matter for realistic nucleon-nucleon interactions such as the Argonne $v_6$ and the Reid $v_6$ potentials. It is shown that diagrammatic contributions that are outside the parquet class play an important role. These are, in variational theories, identified as so-called "commutator contributions". Moreover, using a particle-hole propagator appropriate for a superfluid system results in the suppression of the spin-channel contribution to the induced interaction. Applying these corrections to the pairing interaction, our results agree quite well with Quantum Monte Carlo data.

We develop a manifestly microscopic method to deal with strongly interacting nuclear systems that have different interactions in spin-singlet and spin-triplet states. In a first step we analyze variational wave functions that have been suggested to describe such systems, and demonstrate that the so-called commutator contributions can have important effects whenever the interactions in the spin-singlet and the spin-triplet states are very different. We then identify these contributions as terms that correspond, in the language of perturbation theory, to non-parquet diagrams. We include these diagrams in a way that is suggested by the Jastrow-Feenberg approach and show that the corrections from non-parquet contributions are, at short distances, larger than all other many-body effects.

We apply ideas of the parquet-diagram and optimized Fermi-hypernetted chain methods to determine the short-range structure of the pair wave function in neutron matter and compare these with Bethe-Goldstone results and those of low-order variational calculations. It is shown that the induced interaction, describing the exchange of density, spin, and tensor fluctuations, has a profound influence on the short-ranged structure of the pair wave function and, hence, on effective interactions in neutron matter.

We develop the variational/parquet diagram approach to the structure of nuclear systems with strongly state-dependent interactions. For that purpose, we combine ideas of the general Jastrow-Feenberg variational method and the local parquet-diagram theory for bosons with state-dependent interactions (R. A. Smith and A. D. Jackson, Nucl. Phys. {\bf 476}, 448 (1988)). The most tedious aspect of variational approaches, namely the symmetrization of an operator dependent variational wave function, is thereby avoided. We carry out calculations for neutron matter interacting via the Reid and Argonne $v_6$ models of the nucleon-nucleon interaction. While the equation of state is a rather robust quantity that comes out reasonably well even in very simplistic approaches, we show that effective interactions, which are the essential input for calculating dynamic properties, depend sensitively on the quality of the treatment of the many-body problem.

The dynamic structure function $S(k,\omega)$ informs about the dispersion and damping of excitations. We have recently (Phys. Rev. B {\bf 97}, 184520 (2018)) compared experimental results for $S(k,\omega)$ from high-precision neutron scattering experiment and theoretical results using the ``dynamic many-body theory'' (DMBT), showing excellent agreement over the whole experimentally accessible pressure regime. This paper focuses on the specific aspect of the propagation of low-energy phonons. We report calculations of the phonon mean-free path and phonon life time in liquid \he4 as a function of wave length and pressure. Historically, the question was of interest for experiments of quantum evaporation. More recently, there is interest in the potential use of $^4$He as a detector for low-energy dark matter (K. Schulz and Kathryn M. Zurek, Phys. Rev. Lett. {\bf 117}, 121302 (2016)). While the mean free path of long wave length phonons is large, phonons of intermediate energy can have a short mean free path of the order of $\mu$m. Comparison of different levels of theory indicate that reliable predictions of the phonon mean free path can be made only by using the most advanced many--body method available, namely, DMBT.

We report calculations of equation of state of a model system, representative for a cold Fermi gas, of particles interacting via the P\"oschl-Teller interaction. In successively more sophisticated calculations, we determine the importance of correlations and non-localities. We show that the gas displays, at relatively low density, an instability indicated by a divergence of the in-medium scattering length which occurs well before the divergence of the vacuum scattering length and the spinodal density. We also calculate the pairing gap and show that non-local correlations can change the pairing gap by almost an order of magnitude., Comment: 9 pages, 2 figures, QFS 2018 in Tokyo

The dynamic structure factor of superfluid $^4$He has been investigated at very low temperatures by inelastic neutron scattering. The measurements combine different incoming energies resulting in an unprecedentedly large dynamic range with excellent energy resolution, covering wave vectors $Q$ up to 5 \AA$^{-1}$ and energies $\omega$ up to 15 meV. A detailed description of the dynamics of superfluid $^4$He is obtained from saturated vapor pressure up to solidification. The single-excitation spectrum is substantially modified at high pressures, as the maxon energy exceeds the roton-roton decay threshold. A highly structured multi-excitation spectrum is observed at low energies, where clear thresholds and branches have been identified. Strong phonon emission branches are observed when the phonon or roton group velocities exceed the sound velocity. The spectrum is found to display strong multi-excitations whenever the single-excitations face disintegration following Pitaevskii's type a or b criteria. At intermediate energies, an interesting pattern in the dynamic structure factor is observed in the vicinity of the recoil energy. All these features, which evolve significantly with pressure, are in very good agreement with the Dynamic Many-body calculations, even at the highest densities, where the correlations are strongest.

We report a theoretical analysis of variational wave functions for the BCS pairing problem. Starting with a Jastrow-Feenberg (or, in a more recent language "fixed-node") wave function for the superfluid state, we develop the full optimized Fermi-Hypernetted Chain (FHNC-EL) equations which sum a local approximation of the parquet-diagrams. Close examination of the procedure reveals that it is essential to go beyond the usual Jastrow-Feenberg approximation to guarantee the correct stability range., Comment: 12 pages, 3 figures

We report calculations of the superfluid pairing gap in neutron matter for the $^1S_0$ components of the Reid soft-core $V_6$ and the Argonne $V_{4}'$ two-nucleon interactions. Ground-state calculations have been carried out using the central part of the operator-basis representation of these interactions to determine optimal Jastrow-Feenberg correlations and corresponding effective pairing interactions within the correlated-basis formalism (CBF), the required matrix elements in the correlated basis being evaluated by Fermi hypernetted-chain techniques. Different implementations of the Fermi-Hypernetted Chain Euler-Lagrange method (FHNC-EL) agree at the percent level up to nuclear matter saturation density. For the assumed interactions, which are realistic within the low density range involved in $^1S_0$ neutron pairing, we did not find a dimerization instability arising from divergence of the in-medium scattering length, as was reported recently for simple square-well and Lennard-Jones potential models (Phys. Rev. A {\bf 92}, 023640 (2015))., Comment: 25 pages, 6 figures

The dynamics of superfluid 4He at and above the Landau quasiparticle regime is investigated by high precision inelastic neutron scattering measurements of the dynamic structure factor. A highly structured response is observed above the familiar phonon-maxon-roton spectrum, characterized by sharp thresholds for phonon-phonon, maxon-roton and roton-roton coupling processes. The experimental dynamic structure factor is compared to the calculation of the same physical quantity by a Dynamic Many-body theory including three-phonon processes self-consistently. The theory is found to provide a quantitative description of the dynamics of the correlated bosons for energies up to about three times that of the Landau quasiparticles., Comment: 5 pages, 3 figures

We present ground state calculations for low-density Fermi gases described by two model interactions, an attractive square-well potential and a Lennard-Jones potential, of varying strength. We use the optimized Fermi-Hypernetted Chain integral equation method which has been proved to provide, in the density regimes of interest here, an accuracy better than one percent. We first examine the low-density expansion of the energy and compare with the exact answer by Huang and Yang (H. Huang and C. N. Yang, {\em Phys. Rev.\/} {\bf 105}, 767 (1957)). It is shown that a locally correlated wave function of the Jastrow-Feenberg type does not recover the quadratic term in the expansion of the energy in powers of $\a0\KF$, where $\a0$ is the vacuum $s$-wave scattering length and $\KF$ the Fermi wave number. The problem is cured by adding second-order perturbation corrections in a correlated basis. Going to higher densities and/or more strongly coupled systems, we encounter an instability of the normal state of the system which is characterized by a divergence of the {\em in-medium\/} scattering length. We interpret this divergence as a phonon-exchange driven dimerization of the system, similar to what one has at zero density when the vacuum scattering length $\a0$ diverges. We then study, in the stable regime, the superfluid gap and its dependence on the density and the interaction strength. We identify two different corrections to low-density expansions: One is medium corrections to the pairing interaction, and the other one finite-range corrections. We show that the most important finite-range corrections are a direct manifestation of the many-body nature of the system., Comment: 24 pages, 18 figures

We present a theoretical study of the dynamic structure function of a resonantly interacting two-component Fermi gas at zero temperature. Our approach is based on dynamic many-body theory able to describe excitations in strongly correlated Fermi systems. The fixed-node diffusion Monte Carlo method is used to produce the ground-state correlation functions which are used as an input for the excitation theory. Our approach reproduces recent Bragg scattering data in both the density and the spin channel. In the BCS regime, the response is close to that of the ideal Fermi gas. On the BEC side, the Bose peak associated with the formation of dimers dominates the density channel of the dynamic response. When the fraction of dimers is large our theory departs from the experimental data, mainly in the spin channel., Comment: 6 pages

We develop a systematic theory of multi-particle excitations in strongly interacting Fermi systems. Our work is the generalization of the time-honored work by Jackson, Feenberg, and Campbell for bosons, that provides, in its most advanced implementation, quantitative predictions for the dynamic structure function in the whole experimentally accessible energy/momentum regime. Our view is that the same physical effects -- namely fluctuations of the wave function at an atomic length scale -- are responsible for the correct energetics of the excitations in both Bose and Fermi fluids. Besides a comprehensive derivation of the fermion version of the theory and discussion of the approximations made, we present results for homogeneous He-3 and electrons in three dimensions. We find indeed a significant lowering of the zero sound mode in He-3 and a broadening of the collective mode due to the coupling to particle-hole excitations in good agreement with experiments. The most visible effect in electronic systems is the appearance of a ``double-plasmon'' excitation., Comment: submitted to Phys. Rev. B

The eigenvalue-function pair of the 3D Schr\"odinger equation can be efficiently computed by use of high order, imaginary time propagators. Due to the diffusion character of the kinetic energy operator in imaginary time, algorithms developed so far are at most fourth-order. In this work, we show that for a grid based algorithm, imaginary time propagation of any even order can be devised on the basis of multi-product splitting. The effectiveness of these algorithms, up to the 12$^{\rm th}$ order, is demonstrated by computing all 120 eigenstates of a model C$_{60}$ molecule to very high precisions. The algorithms are particularly useful when implemented on parallel computer architectures., Comment: 8 pages, 3 figures

We investigate the effects of impurities and changing ring geometry on the energetics of quantum rings under different magnetic field strengths. We show that as the magnetic field and/or the electron number are/is increased, both the quasiperiodic Aharonov-Bohm oscillations and various magnetic phases become insensitive to whether the ring is circular or square in shape. This is in qualitative agreement with experiments. However, we also find that the Aharonov-Bohm oscillation can be greatly phase-shifted by only a few impurities and can be completely obliterated by a high level of impurity density. In the many-electron calculations we use a recently developed fourth-order imaginary time projection algorithm that can exactly compute the density matrix of a free-electron in a uniform magnetic field., Comment: 8 pages, 7 figures, to appear in PRB

We investigate positron annihilation in electron liquid as a case study for many-body theory, in particular the optimized Fermi Hypernetted Chain (FHNC-EL) method. We examine several approximation schemes and show that one has to go up to the most sophisticated implementation of the theory available at the moment in order to get annihilation rates that agree reasonably well with experimental data. Even though there is basically just one number to look at, the electron-positron pair distribution function at zero distance, it is exactly this number that dictates how the full pair distribution behaves: In most cases, it falls off monotonously towards unity as the distance increases. Cases where the electron-positron pair distribution exhibits a dip are precursors to the formation of bound electron--positron pairs. The formation of electron-positron pairs is indicated by a divergence of the FHNC-EL equations, from this we can estimate the density regime where positrons must be localized. This occurs in our calculations in the range 9.4 <= r_s <=10, where r_s is the dimensionless density parameter of the electron liquid., Comment: To appear in Phys. Rev. B (2003)

We use structural information from diffusion Monte Carlo calculations for two--dimensional 3He to calculate the effective mass. Static effective interactions are constructed from the density-- and spin structure functions using sumrules. We find that both spin-- and density-- fluctuations contribute about equally to the effective mass. Our results show, in agreement with recent experiments, a flattening of the single--particle self--energy with increasing density, which eventually leads to a divergent effective mass., Comment: 4 pages, accepted in PRL

We examine the energetics of bulk and layer-roton excitations of 4He in various porous medial such as aerogel, Geltech, or Vycor, in order to find out what conclusions can be drawn from experiments on the energetics about the physisorption mechanism. The energy of the layer-roton minimum depends sensitively on the substrate strength, thus providing a mechanism for a direct measurement of this quantity. On the other hand, bulk-like roton excitations are largely independent of the interaction between the medium and the helium atoms, but the dependence of their energy on the degree of filling reflects the internal structure of the matrix and can reveal features of 4He at negative pressures. While bulk-like rotons are very similar to their true bulk counterparts, the layer modes are not in close relation to two-dimensional rotons and should be regarded as a third, completely independent kind of excitation.

We design models for helium in matrices like aerogel, Vycor or Geltech from a manifestly microscopic point of view. For that purpose, we calculate the dynamic structure function of 4He on Si substrates and between two Si walls as a function of energy, momentum transfer, and the scattering angle. The angle--averaged results are in good agreement with the neutron scattering data; the remaining differences can be attributed to the simplified model used here for the complex pore structure of the materials. A focus of the present work is the detailed identification of coexisting layer modes and bulk--like excitations, and, in the case of thick films, ripplon excitations. Involving essentially two--dimensional motion of atoms, the layer modes are sensitive to the scattering angle., Comment: Phys. Rev. B (2003, in press)

We describe a rapidly converging algorithm for solving the Kohn--Sham equations and equations of similar structure that appear frequently in calculations of the structure of inhomogeneous electronic many--body systems. The algorithm has its roots the Hohenberg-Kohn theorem and solves directly for the electron density; single--particle wave functions are only used as auxiliary quantities. The method has been implemented for symmetric ``slabs'' of jellium as well as for spherical jellium clusters. Starting from very rough guesses for the initial electron density, convergence is reached within a few iterations. The iterations are driven by the static electric susceptibility., Comment: 13 pages, 2 figures

We calculate microscopically the properties of the dilute He-3 component in a He-3/--He-4 mixture. These depend on both, the dominant interaction between the impurity atom and the background, and the Fermi liquid contribution due to the interaction between the constituents of the He-3 component. We first calculate the dynamic structure function of a He-3 impurity atom moving in He-3. From that we obtain the excitation spectrum and the momentum dependent effective mass. The pole strength of this excitation mode is strongly reduced from the free particle value in agreement with experiments; part of the strength is distributed over high frequency excitations. Above k > 1.7$\AA$^{-1}$ the motion of the impurity is damped due to the decay into a roton and a low energy impurity mode. Next we determine the Fermi--Liquid interaction between He-4 atoms and calculate the pressure-- and concentration dependence of the effective mass, magnetic susceptibility, and the He-3--He-3 scattering phase shifts. The calculations are based on a dynamic theory that uses, as input, effective interactions provided by the Fermi hypernetted--chain theory. The relationship between both theories is discussed. Our theoretical effective masses agree well with recent measurements by Yorozu et al. (Phys. Rev. B 48, 9660 (1993)) as well as those by R. Simons and R. M. Mueller (Czekoslowak Journal of Physics Suppl. 46, 201 (1996)), but our analysis suggests a new extrapolation to the zero-concentration limit. With that effective mass we also find a good agreement with the measured Landau parameter F_0^a., Comment: 47 pages, 15 figures

We develop a first principles, microscopic theory of impurity atom scattering from inhomogeneous quantum liquids such as adsorbed films, slabs, or clusters of He-4. The theory is built upon a quantitative, microscopic description of the ground state of both the host liquid as well as the impurity atom. Dynamic effects are treated by allowing all ground-state correlation functions to be time-dependent. Our description includes both the elastic and inelastic coupling of impurity motion to the excitations of the host liquid. As a specific example, we study the scattering of He-3 atoms from adsorbed He-4 films. We examine the dependence of ``quantum reflection'' on the substrate, and the consequences of impurity bound states, resonances, and background excitations for scattering properties. A thorough analysis of the theoretical approach and the physical circumstances point towards the essential role played by inelastic processes which determine almost exclusively the reflection probabilities. The coupling to impurity resonances within the film leads to a visible dependence of the reflection coefficient on the direction of the impinging particle., Comment: 36 pages, 16 figures

Recent measurements by Yorozu et al. (S. Yorozu, H. Fukuyama, and H. Ishimoto, Phys. Rev. B 48, 9660 (1993)) as well as by Simons and Mueller (R. Simons and R. M. Mueller, Czhechoslowak Journal of Physics Suppl. 46, 201 (1976)) have determined the effective mass of He-3 atoms in a He-3/He-4 mixture with great accuracy. We here report theoretical calculations for the dependence of that effective mass on the He-3 concentration. Using correlated basis functions perturbation theory to infinite order to compute effective interactions in the appropriate channels, we obtain good agreement between theory and experiment., Comment: 4 pages, 1 figure

We develop a microscopic theory of the scattering, transmission, and sticking of 4He atoms impinging on a superfluid 4He slab at near normal incidence, and inelastic neutron scattering from the slab. The theory includes coupling between different modes and allows for inelastic processes. We find a number of essential aspects that must be observed in a physically meaningful and reliable theory of atom transmission and scattering; all are connected with multiparticle scattering, particularly the possibility of energy loss. These processes are (a) the coupling to low-lying (surface) excitations (ripplons/third sound) which is manifested in a finite imaginary part of the self energy, and (b) the reduction of the strength of the excitation in the maxon/roton region.