Optimisation of quantum Monte Carlo wave function: steepest descent method

We have employed the steepest descent method to optimise the variational ground state quantum Monte Carlo wave function for He, Li, Be, B and C atoms. We have used both the direct energy minimisation and the variance minimisation approaches. Our calculations show that in spite of receiving insufficient attention, the steepest descent method can successfully minimise the wave function. All the derivatives of the trial wave function respect to spatial coordinates and variational parameters have been computed analytically. Our ground state energies are in a very good agreement with those obtained with diffusion quantum Monte Carlo method (DMC) and the exact results.Comment: 13 pages, 3 eps figure

Slow Switching in Globally Coupled Oscillators: Robustness and Occurrence through Delayed Coupling

The phenomenon of slow switching in populations of globally coupled oscillators is discussed. This characteristic collective dynamics, which was first discovered in a particular class of the phase oscillator model, is a result of the formation of a heteroclinic loop connecting a pair of clustered states of the population. We argue that the same behavior can arise in a wider class of oscillator models with the amplitude degree of freedom. We also argue how such heteroclinic loops arise inevitably and persist robustly in a homogeneous population of globally coupled oscillators. Although the heteroclinic loop might seem to arise only exceptionally, we find that it appears rather easily by introducing the time-delay in the population which would otherwise exhibit perfect phase synchrony. We argue that the appearance of the heteroclinic loop induced by the delayed coupling is then characterized by transcritical and saddle-node bifurcations. Slow switching arises when the system with a heteroclinic loop is weakly perturbed. This will be demonstrated with a vector model by applying weak noises. Other types of weak symmetry-breaking perturbations can also cause slow switching.Comment: 10 pages, 14 figures, RevTex, twocolumn, to appear in Phys. Rev.

Superconductivity in the fullerenes

Wetensch. publicatieFaculteit der Wiskunde en Natuurwetenschappe

Effect of randomness and anisotropy on Turing patterns in reaction-diffusion systems

We study the effect of randomness and anisotropy on Turing patterns in reaction-diffusion systems. For this purpose, the Gierer-Meinhardt model of pattern formation is considered. The cases we study are: (i)randomness in the underlying lattice structure, (ii)the case in which there is a probablity p that at a lattice site both reaction and diffusion occur, otherwise there is only diffusion and lastly, the effect of (iii) anisotropic and (iv) random diffusion coefficients on the formation of Turing patterns. The general conclusion is that the Turing mechanism of pattern formation is fairly robust in the presence of randomness and anisotropy.Comment: 11 pages LaTeX, 14 postscript figures, accepted in Phys. Rev.

Balancing Minimum Spanning and Shortest Path Trees

This paper give a simple linear-time algorithm that, given a weighted digraph, finds a spanning tree that simultaneously approximates a shortest-path tree and a minimum spanning tree. The algorithm provides a continuous trade-off: given the two trees and epsilon > 0, the algorithm returns a spanning tree in which the distance between any vertex and the root of the shortest-path tree is at most 1+epsilon times the shortest-path distance, and yet the total weight of the tree is at most 1+2/epsilon times the weight of a minimum spanning tree. This is the best tradeoff possible. The paper also describes a fast parallel implementation.Comment: conference version: ACM-SIAM Symposium on Discrete Algorithms (1993

Magnetic Properties of Undoped C60C_{60}

The Heisenberg antiferromagnet, which arises from the large $U$ Hubbard model, is investigated on the $C_{60}$ molecule and other fullerenes. The connectivity of $C_{60}$ leads to an exotic classical ground state with nontrivial topology. We argue that there is no phase transition in the Hubbard model as a function of $U/t$, and thus the large $U$ solution is relevant for the physical case of intermediate coupling. The system undergoes a first order metamagnetic phase transition. We also consider the S=1/2 case using perturbation theory. Experimental tests are suggested.Comment: 12 pages, 3 figures (included

Analytical results for coupled map lattices with long-range interactions

We obtain exact analytical results for lattices of maps with couplings that decay with distance as $r^{-\alpha}$. We analyze the effect of the coupling range on the system dynamics through the Lyapunov spectrum. For lattices whose elements are piecewise linear maps, we get an algebraic expression for the Lyapunov spectrum. When the local dynamics is given by a nonlinear map, the Lyapunov spectrum for a completely synchronized state is analytically obtained. The critical lines characterizing the synchronization transition are determined from the expression for the largest transversal Lyapunov exponent. In particular, it is shown that in the thermodynamical limit, such transition is only possible for sufficiently long-range interactions, namely, for $\alpha\le alpha_c<d$, where $d$ is the lattice dimension.Comment: 4 pages, 2 figures, corrections included. Phys. Rev. E 68, 045202(R) (2003); correction in pres

Vibrational signatures for low-energy intermediate-sized Si clusters

We report low-energy locally stable structures for the clusters Si20 and Si21. The structures were obtained by performing geometry optimizations within the local density approximation. Our calculated binding energies for these clusters are larger than any previously reported for this size regime. To aid in the experimental identification of the structures, we have computed the full vibrational spectra of the clusters, along with the Raman and IR activities of the various modes using a recently developed first-principles technique. These represent, to our knowledge, the first calculations of Raman and IR spectra for Si clusters of this size

High pressure diamond-like liquid carbon

We report density-functional based molecular dynamics simulations, that show that, with increasing pressure, liquid carbon undergoes a gradual transformation from a liquid with local three-fold coordination to a 'diamond-like' liquid. We demonstrate that this unusual structural change is well reproduced by an empirical bond order potential with isotropic long range interactions, supplemented by torsional terms. In contrast, state-of-the-art short-range bond-order potentials do not reproduce this diamond structure. This suggests that a correct description of long-range interactions is crucial for a unified description of the solid and liquid phases of carbon.Comment: 4 pages, 5 figure

Surface reconstruction induced geometries of Si clusters

We discuss a generalization of the surface reconstruction arguments for the structure of intermediate size Si clusters, which leads to model geometries for the sizes 33, 39 (two isomers), 45 (two isomers), 49 (two isomers), 57 and 61 (two isomers). The common feature in all these models is a structure that closely resembles the most stable reconstruction of Si surfaces, surrounding a core of bulk-like tetrahedrally bonded atoms. We investigate the energetics and the electronic structure of these models through first-principles density functional theory calculations. These models may be useful in understanding experimental results on the reactivity of Si clusters and their shape as inferred from mobility measurements.Comment: 9 figures (available from the author upon request) Submitted to Phys. Rev.