An investigation of the 2D attractive Hubbard model

We present an investigation of the 2D attractive Hubbard model, considered as an effective model relevant to superconductivity in strongly interacting electron systems. We use both hybrid Monte-Carlo simulations and existing hopping parameter expansions to explore the low temperature domain. The increase of the static S-wave pair correlation with decreasing temperature, which depends weakly on the band filling in the explored temperature range, is analyzed in terms of an expected Kosterlitz-Thouless superconducting transition. Using both our data and previously published results, we show that the evidence for this transition is weak: If it exists, its temperature is very low. The number of unpaired electrons remains nearly constant with temperature at fixed attractive potential strength. In contrast, the static magnetic susceptibility decreases fast with temperature, and cannot be related only to pair formation. We introduce a method by which the Padé approximants of the existing series for the susceptibility give sensible results down to rather low temperature region, as shown by comparison with our numerical data. Received: 30 October 1996 / Revised: 23 October 1997 / Accepted: 29 January 1998     

Combinational simulations of free energy of polymer solution

It is conceptually proposed that the total entropy of polymer solution is contributed from two distinct parts: the positional and the oomformational. The former can be represented analytically, while the latter can be simulated with the random self-avoiding walk model on the simple cubic lattice for multichain systems. The obtained results indicated that both the conformational entropy and the mixing heat are consistent with the scaling laws wry well.     

Recognition of ring pairs in the data analysis of Cherenkov detectors

We present a new algorithm for the identification of ring pairs in Cherenkov detectors. The approach combines linear regression methods for the construction of an appropriate reference quantity, methods from statistics for the evaluation of hypotheses and Monte Carlo simulation runs for the tuning of threshold levels.     

High-precision Monte Carlo test of the conformai-invariance predictions for two-dimensional mutually avoiding walks

Let _l be the critical exponent associated with the probability thatl independentN-step ordinary random walks, starting at nearby points, are mutually avoiding. Using Monte Carlo methods combined with a maximum-likelihood data analysis, we find that in two dimensions ₂=0.6240±0.0005±0.0011 and ₃=1.4575±0.0030±0.0052, where the first error bar represents systematic error due to corrections to scaling (subjective 95% confidence limits) and the second error bar represents statistical error (classical 95% confidence limits). These results are in good agreement with the conformal-invariance predictions ₂=5/8 and ₃=35/24.     

Branched polymers in a wedge geometry in three dimensions

We investigate the statistical and dimensional properties of uniform star polymers attached by the branching vertex of degreef in a wedge geometry in three dimensions and described by the wedge angles and. We show that the growth constant is equal to ^f , where is the self-avoiding walk limit. Thef and (, ) dependences of the corresponding critical exponent _f (, ) are studied using Monte Carlo techniques. In the casef=1, our results are compared with existing predictions obtained from series expansion and renormalization group methods. We have also estimated the amplitudes for the mean square radius of gyration and the mean square end-to-end branch length. Our results for the ratio of the mean square radius of gyration of anf-star to that of a linear polymer of the same degree of polymerization attached in a similar wedge, and the analogous ratio for the mean square end-to-end branch length, are consistent with these ratios being lattice-independent quantities.     

Superdegenerate point in FCC phase diagram: CVM and Monte Carlo investigations

We investigate the topology of the phase diagram of binary alloys on the fee lattice with first-neighbor antiferromagnetic interactions around the superdegenerate point, where the L1₀ and L1₂ phases meet. We treat the system as a hard-constraint lattice gas, following a procedure previously described by Lebowitzet al. We perform cluster variation method calculations in theT0 limit and Monte Carlo simulations directly atT=0 K on the ground states of the superdegenerate point. We find that: (i) there is no disordered phase in the neighborhood of this point; (ii) a phase L for which two of the four cubic sublattices have the same average occupancy and each of the two others are different appears between L1₀ and L1₂; (iii) the transition L/L1₂ is of first order.     

Likelihood estimation of soft-core interaction potentials for Gibbsian point patterns

The likelihood method is developed for the analysis of socalled regular point patterns. Approximating the normalizing factor of Gibbs canonical distribution, we simultaneously estimate two parameters, one for the scale and the other which measures the softness (or hardness), of repulsive interactions between points. The approximations are useful up to a considerably high density. Some real data are analyzed to illustrate the utility of the parameters for characterizing the regular point pattern.     

Computer simulation of some dynamical properties of the Lorentz gas

We carried out molecular dynamics simulations of a Lorentz gas, consisting of a lone hydrogen molecule moving in a sea of stationary argon atoms. A Lennard-Jones form was assumed for the H₂-Ar potential. The calculations were performed at a reduced temperatureK ^* =kT/_{H 2–Ar} = 4.64 and at reduced densities ^*= _{Ar Ar} ³ in the range 0.074–0.414. The placement of Ar atoms was assumed to be random rather than dictated by equilibrium considerations. We followed the trajectories of many H₂ molecules, each of which is assigned in turn a velocity given by the Maxwell-Boltzmann distribution at the temperature of the simulation. Solving the equations of motion classically, we obtained the translational part of the incoherent dynamic structure factor for the H₂ molecule,S _tr(q, ). This was convoluted with the rotational structure factorS _rot(q, ) calculated assuming unhindered rotation to obtain the total structure factorS(q, ). Our results agree well with experimental data on this function obtained by Egelstaffet al. At the highest density ( ^*=0.414) we studied the dependence ofS(q, ) on system size (number of Ar atoms), number of H₂ molecules for which trajectories are generated, and the length of time over which these trajectories are followed.     

Surface properties of finite classical Coulomb systems: Debye-Hückel approximation and computer simulations

We report analytical and numerical studies of surface correlations in finite, homogeneously polarizable, classical Coulomb systems placed in an insulating or conducting environment. Their purpose is to understand the phenomenological, shape-dependent laws of electrostatics, from the point of view of statistical mechanics; we focus on the knowledge of the dielectric susceptibility of the system, a quantity proportional to the equilibrium fluctuation of the system's instantaneous polarization per unit volume. This goal has been achieved for a system in a conducting state. The picture is that the shape-dependent part of the susceptibilities results from the action of unbounded observables (the second moments of the instantaneous polarization of the system) on long-range surface correlations and that the relations of electrostatics are verified by means of shape-dependent thermodynamic limits. This picture is supported (i) by exact solutions and asymptotic analysis of the Debye-Hückel approximation of multicomponent plasmas in disks and spheres with insulating and conducting environment and also in ellipses in a vacuum, and (ii) by computer simulations of a one-component plasma in a disk with different environments, notably a conducting environment with permeable and impermeable wall. These observations have revealed for the first time the reason why the susceptibility of a conducting disk in a conductor with impermeable walls diverges linearly with the radius of the disk: this is due to the occurrence of long-range radial correlations in the conductor. These findings are quantitatively interpreted in terms of a novel canonical Debye-Huckel approximation as contrasted to the ordinary grand canonical version. Lastly a fresh look at the problem of the surface correlations of a conductor in a vacuum, which places the observer close to the surface of the conductor but in the vacuum, is presented and applied to the disk, the ellipse, the cylinder, the sphere, and the wedge.     

On the spectrum of the Heisenberg Hamiltonian

The quantum, antiferromagnetic, spin-1/2 Heisenberg Hamiltonian on thed-dimensional cubic lattice ^d is considered for any dimensiond. First the anisotropic case is considered for small transversal coupling and a convergent expansion is given for a family of eigenprojections which is complete in all finite-volume truncations. Then the general case is considered, for which an upper bound to the ground-state energy is given which is optimal for strong enough anisotropy. This bound is expressed through a functional involving the statistical expectation value at finite temperature of a certain correlation function of an Ising model defined on the lattice ^d itself.