Variational Methods for Nuclear Systems with Dynamical Mesons

We derive a model Hamiltonian whose ground state expectation value of any two-body operator coincides with that obtained with the Jastrow correlated wave function of the many-body Fermi system. Using this Hamiltonian we show that the variational principle can be extended to treat systems with dynamical mesons, even if in this case the concept of wave function looses its meaning

Computer simulation study of the closure relations in hard sphere fluids

We study, using Monte Carlo simulations, the cavity and the bridge functions of various hard sphere fluids: one component system, equimolar additive and non additive binary mixtures. In particular, we numerically check the assumption of local dependency of the bridge functions from the indirect correlation functions, on which most of the existing integral equation theories hinge. We find that this condition can be violated either in the region around the first and second neighbors shell, or inside the hard core, for the systems here considered. The violations manifest themselves clearly in the so called Duh-Haymet plots of the bridge functions versus the indirect correlation functions and become amplified as the coupling of the system increases.Comment: 24 pages, 13 figure

Quantum Monte Carlo Algorithm Based on Two-Body Density Functional Theory for Fermionic Many-Body Systems: Application to 3He

We construct a quantum Monte Carlo algorithm for interacting fermions using the two-body density as the fundamental quantity. The central idea is mapping the interacting fermionic system onto an auxiliary system of interacting bosons. The correction term is approximated using correlated wave functions for the interacting system, resulting in an effective potential that represents the nodal surface. We calculate the properties of 3He and find good agreement with experiment and with other theoretical work. In particular, our results for the total energy agree well with other calculations where the same approximations were implemented but the standard quantum Monte Carlo algorithm was usedComment: 4 pages, 3 figures, 1 tabl

Phase diagram of the penetrable square well-model

We study a system formed by soft colloidal spheres attracting each other via a square-well potential, using extensive Monte Carlo simulations of various nature. The softness is implemented through a reduction of the infinite part of the repulsive potential to a finite one. For sufficiently low values of the penetrability parameter we find the system to be Ruelle stable with square-well like behavior. For high values of the penetrability the system is thermodynamically unstable and collapses into an isolated blob formed by a few clusters each containing many overlapping particles. For intermediate values of the penetrability the system has a rich phase diagram with a partial lack of thermodynamic consistency.Comment: 6 pages and 5 figure

A Numerical Test of a High-Penetrability Approximation for the One-Dimensional Penetrable-Square-Well Model

The one-dimensional penetrable-square-well fluid is studied using both analytical tools and specialized Monte Carlo simulations. The model consists of a penetrable core characterized by a finite repulsive energy combined with a short-range attractive well. This is a many-body one-dimensional problem, lacking an exact analytical solution, for which the usual van Hove theorem on the absence of phase transition does not apply. We determine a high-penetrability approximation complementing a similar low-penetrability approximation presented in previous work. This is shown to be equivalent to the usual Debye-H\"{u}ckel theory for simple charged fluids for which the virial and energy routes are identical. The internal thermodynamic consistency with the compressibility route and the validity of the approximation in describing the radial distribution function is assessed by a comparison against numerical simulations. The Fisher-Widom line separating the oscillatory and monotonic large-distance behavior of the radial distribution function is computed within the high-penetrability approximation and compared with the opposite regime, thus providing a strong indication of the location of the line in all possible regimes. The high-penetrability approximation predicts the existence of a critical point and a spinodal line, but this occurs outside the applicability domain of the theory. We investigate the possibility of a fluid-fluid transition by Gibbs ensemble Monte Carlo techniques, not finding any evidence of such a transition. Additional analytical arguments are given to support this claim. Finally, we find a clustering transition when Ruelle's stability criterion is not fulfilled. The consequences of these findings on the three-dimensional phase diagrams are also discussed.Comment: 17 pages, 12 figures; to be published in JC

Generating functionals, consistency, and uniqueness in the integral equation theory of liquids

We discuss and illustrate through numerical examples the relations between generating functionals, thermodynamic consistency (in particular the virial-free energy one), and uniqueness of the solution, in the integral equation theory of liquids. We propose a new approach for deriving closures automatically satisfying such characteristics. Results from a first exploration of this program are presented and discussed.Comment: 27 pages, 5 figure

The spectra of mixed 3^3He-4^4He droplets

The diffusion Monte Carlo technique is used to calculate and analyze the excitation spectrum of $^3$He atoms bound to a cluster of $^4$He atoms, by using a previously determined optimum filling of single-fermion orbits with well defined orbital angular momentum $L$, spin $S$ and parity quantum numbers. The study concentrates on the energies and shapes of the three kinds of states for which the fermionic part of the wave function is a single Slater determinant: maximum $L$ or maximum $S$ states within a given orbit, and fully polarized clusters. The picture that emerges is that of systems with strong shell effects whose binding and excitation energies are essentially determined over configuration at fixed number of particles and spin, i.e., by the monopole properties of an effective Hamiltonian.Comment: 14 pages, 15 figure

Phase oscillations in superfluid 3He-B weak links

Oscillations in quantum phase about a mean value of $\pi$, observed across micropores connecting two \helium baths, are explained in a Ginzburg-Landau phenomenology. The dynamics arises from the Josephson phase relation,the interbath continuity equation, and helium boundary conditions. The pores are shown to act as Josephson tunnel junctions, and the dynamic variables are the inter bath phase difference and fractional difference in superfluid density at micropores. The system maps onto a non-rigid, momentum-shortened pendulum, with inverted-orientation oscillations about a vertical tilt angle $\phi = \pi$, and other modes are predicted