Superconductivity in hole-doped C60 from electronic correlations

We derive a model for the highest occupied molecular orbital band of a C60 crystal which includes on-site electron-electron interactions. The form of the interactions are based on the icosahedral symmetry of the C60 molecule together with a perturbative treatment of an isolated C60 molecule. Using this model we do a mean-field calculation in two dimensions on the [100] surface of the crystal. Due to the multi-band nature we find that electron-electron interactions can have a profound effect on the density of states as a function of doping. The doping dependence of the transition temperature can then be qualitatively different from that expected from simple BCS theory based on the density of states from band structure calculations

Unmixing in Random Flows

We consider particles suspended in a randomly stirred or turbulent fluid. When effects of the inertia of the particles are significant, an initially uniform scatter of particles can cluster together. We analyse this 'unmixing' effect by calculating the Lyapunov exponents for dense particles suspended in such a random three-dimensional flow, concentrating on the limit where the viscous damping rate is small compared to the inverse correlation time of the random flow (that is, the regime of large Stokes number). In this limit Lyapunov exponents are obtained as a power series in a parameter which is a dimensionless measure of the inertia. We report results for the first seven orders. The perturbation series is divergent, but we obtain accurate results from a Pade-Borel summation. We deduce that particles can cluster onto a fractal set and show that its dimension is in satisfactory agreement with previously reported in simulations of turbulent Navier-Stokes flows. We also investigate the rate of formation of caustics in the particle flow.Comment: 39 pages, 8 figure

Recurrent Variational Approach to the Two-Leg Hubbard Ladder

We applied the Recurrent Variational Approach to the two-leg Hubbard ladder. At half-filling, our variational Ansatz was a generalization of the resonating valence bond state. At finite doping, hole pairs were allowed to move in the resonating valence bond background. The results obtained by the Recurrent Variational Approach were compared with results from Density Matrix Renormalization Group.Comment: 10 pages, 14 Postscript figure

A class of ansatz wave functions for 1D spin systems and their relation to DMRG

We investigate the density matrix renormalization group (DMRG) discovered by White and show that in the case where the renormalization eventually converges to a fixed point the DMRG ground state can be simply written as a ``matrix product'' form. This ground state can also be rederived through a simple variational ansatz making no reference to the DMRG construction. We also show how to construct the ``matrix product'' states and how to calculate their properties, including the excitation spectrum. This paper provides details of many results announced in an earlier letter.Comment: RevTeX, 49 pages including 4 figures (macro included). Uuencoded with uufiles. A complete postscript file is available at http://fy.chalmers.se/~tfksr/prb.dmrg.p

Fixed Point of the Finite System DMRG

The density matrix renormalization group (DMRG) is a numerical method that optimizes a variational state expressed by a tensor product. We show that the ground state is not fully optimized as far as we use the standard finite system algorithm, that uses the block structure B**B. This is because the tensors are not improved directly. We overcome this problem by using the simpler block structure B*B for the final several sweeps in the finite iteration process. It is possible to increase the numerical precision of the finite system algorithm without increasing the computational effort.Comment: 6 pages, 4 figure

The critical behaviour of the 2D Ising model in Transverse Field; a Density Matrix Renormalization calculation

We have adjusted the Density Matrix Renormalization method to handle two dimensional systems of limited width. The key ingredient for this extension is the incorporation of symmetries in the method. The advantage of our approach is that we can force certain symmetry properties to the resulting ground state wave function. Combining the results obtained for system sizes up-to $30 \times 6$ and finite size scaling, we derive the phase transition point and the critical exponent for the gap in the Ising model in a Transverse Field on a two dimensional square lattice.Comment: 9 pages, 8 figure

Conformations of Linear DNA

We examine the conformations of a model for under- and overwound DNA. The molecule is represented as a cylindrically symmetric elastic string subjected to a stretching force and to constraints corresponding to a specification of the link number. We derive a fundamental relation between the Euler angles that describe the curve and the topological linking number. Analytical expressions for the spatial configurations of the molecule in the infinite- length limit were obtained. A unique configuraion minimizes the energy for a given set of physical conditions. An elastic model incorporating thermal fluctuations provides excellent agreement with experimental results on the plectonemic transition.Comment: 5 pages, RevTeX; 6 postscript figure

A Density Matrix Algorithm for 3D Classical Models

We generalize the corner transfer matrix renormalization group, which consists of White's density matrix algorithm and Baxter's method of the corner transfer matrix, to three dimensional (3D) classical models. The renormalization group transformation is obtained through the diagonalization of density matrices for a cubic cluster. A trial application for 3D Ising model with m=2 is shown as the simplest case.Comment: 15 pages, Latex(JPSJ style files are included), 8 ps figures, submitted to J. Phys. Soc. Jpn., some references are correcte

Magnetic Field Effect in Josephson tunneling between d-Wave Superconductors

The magnetic field effect in the Josephson tunneling between two d-wave superconductors are investigated. When the crystal orientation of one (or each) superconductor relative to the interface normal is such that midgap states exist at the interface, there is a component of the tunneling current due to the midgap states. For a junction with a flat {001}|{110} or {100}|{110} interface, this component is the predominant contribution to the current. The predicted current-field dependence differs entirely from the conventional Fraunhofer pattern, in agreement with a published measurement. This is because, apart from the Fraunhofer factor, the critical current depends on the magnetic field B through the current density also which is a linear function of B for weak B.Comment: 5 pages, 2 figure