# ALPS 2 Tutorials:DWA-01 SpecificHeat

From ALPS

# Specific Heat

In this tutorial, we will illustrate how to calculate the specific heat of a **2D boson Hubbard model** of square lattice geometry using the directed worm algorithm.

## with fixed chemical potential

### Using the command line

We consider, for instance, the parameter file parm1a that sets up the directed worm algorithm for simulating the boson Hubbard model of a 20^{2} lattice with square geometry.

LATTICE="square lattice" MODEL="boson Hubbard" L=20 Nmax=20 t=1. U=16. mu=32. THERMALIZATION=10000 SWEEPS=100000 SKIP=400 {T=0.7} {T=0.8} {T=0.9} {T=1.0} {T=1.1} {T=1.2} {T=1.3} {T=1.4} {T=1.5} {T=1.6} {T=1.7} {T=1.8} {T=1.9} {T=2.0}

We want to