We present release 2.0 of the ALPS (Algorithms and Libraries for Physics Simulations) project, an open source software project to develop libraries and application programs for the simulation of strongly correlated quantum lattice models such as quantum magnets, lattice bosons, and strongly correlated fermion systems. The code development is centered on common XML and HDF5 data formats, libraries to simplify and speed up code development, common evaluation and plotting tools, and simulation programs. The programs enable non-experts to start carrying out serial or parallel numerical simulations by providing basic implementations of the important algorithms for quantum lattice models: classical and quantum Monte Carlo (QMC) using non-local updates, extended ensemble simulations, exact and full diagonalization (ED), the density matrix renormalization group (DMRG) both in a static version and a dynamic time-evolving block decimation (TEBD) code, and quantum Monte Carlo solvers for dynamical mean field theory (DMFT). The ALPS libraries provide a powerful framework for programers to develop their own applications, which, for instance, greatly simplify the steps of porting a serial code onto a parallel, distributed memory machine. Major changes in release 2.0 include the use of HDF5 for binary data, evaluation tools in Python, support for the Windows operating system, the use of CMake as build system and binary installation packages for Mac OS X and Windows, and integration with the VisTrails workflow provenance tool. The software is available from our web server at http://alps.comp-phys.org/.

In the Mott insulating phase of the transition metal oxides, the effective orbital-orbital interaction is directional both in the orbital space and in the real space. We discuss a classical realization of directional coupling in two dimensions. Despite extensive degeneracy of the ground state, the model exhibits partial orbital ordering in the form of directional ordering of fluctuations at low temperatures stabilized by an entropy gap. Transition to the disordered phase is shown to be in the Ising universality class through exact mapping and multicanonical Monte Carlo simulations.PACS numbers: 64.60. Cn,05.10.Ln,75.30.Et,75.40.Cx Recently, there has been growing interest in the effects of orbital degeneracy in the physics of transition metal oxide insulators [1,2]. In these systems, the dominating energy scales for d-electrons on transition metal (TM) ions are the on-site Coulomb repulsion (which freezes out the charge degrees of freedom), the Hund's rule coupling, and the crystal field due to the surrounding oxygen ions. The latter two together determine the degeneracies and degrees of freedom of spin and orbital on each transition metal ions. Spins and orbitals on neighboring TM ions can then be coupled through the superexchange mechanism. In the case of orbitals, they are also coupled through the phonon-mediated cooperative JahnTeller mechanism [1]. These couplings determine the low temperature properties of these systems.Because orbital coupling is intrinsically directional[3], orbital ordering brings up some unusual questions. Especially interesting is when the coupling along a given bond direction is Ising like, but with different Ising axes along different bond directions [4]. The Hamiltonian for orbitals is then given bywhere τ is an isospin operator representing the orbital degree of freedom, and n ij is a unit vector giving the Ising axis for the bond ij . For example, for e g orbitals in Perovskite structures, n ij for the three different bond directions are coplanar and oriented relative to each other by 120• , giving rise to the so-called 120• Model [4,5,6]. This model is also applicable to t 2g orbitals on the three bonds of a honeycomb lattice as in planes of V 2 O 3 [7,8]. On the other hand, for t 2g orbitals on Perovskite structures, the relevant model is the Compass Model [1,4,10] with n ij = x, y, z for the three bond directions. The Compass Model may also feature as part of the spinspin coupling of t 2g orbitals when spin-orbit interaction is taken into account [11]. A common feature of both the Compass Model and the 120• Model is the competition between bonds in different directions, with the resulting frustration leading to macroscopic degeneracy of the classical ground state.In this Letter, we report analytical and numerical results on a classical version of (1) in 2D. The highly anisotropic coupling gives rise to interesting interplay between one-and two-dimensional ordering, between continuous and discrete spin physics, and between slow and fast modes. Our main result is that at low ...

We study doping on the Kagome lattice by exploring the t − J-model with variational Monte-Carlo. We use a number of Gutzwiller projected spin-liquid and valence bond-crystal states and compare their energies at several system-sizes. We find that introducing mobile holes drives the system away from the Spin-Liquid state proposed by Ran et al for the undoped system, towards a uniform state with zero-flux. On top of the uniform-state a VBC of the Hastings-type is formed for low doping. The results are compared to exact diagonalization on small clusters. This agrees well.

We investigate the ground state properties of the Kagome lattice t − J model at low doping. We propose a new state of matter. Our results suggest an interesting phase with a balance of spinexchange and kinetic exchange through the building blocks of stars and their internal hexagons, being linked by triangles. The particles favor to hop in the hexagons, while the spin-exchange is taking place mainly on the stars. The combined effect leads to the formation of the "Star of David" valence bond crystal. We discuss our result in connection with static impurities, and show the likely relevance to the diluted Kagome lattice Heisenberg model, describing actual compounds.