We consider the Rational Hybrid Monte Carlo algorithm for performing exact 2+1 flavour fermion simulations. The specific cases of ASQTAD and domain wall fermions are considered. We find that in both cases the naive performance is similar to conventional hybrid algorithms.
Filtering algorithms for two degenerate quark flavours have advanced to the point that, in 2+1 flavour simulations, the cost of the strange quark is significant compared with the light quarks. This makes efficient filtering algorithms for single flav
our actions highly desirable, in particular when considering 1+1+1 flavour simulations for QED+QCD. Here we discuss methods for filtering the RHMC algorithm that are implemented within BQCD, an open-source Fortran program for Hybrid Monte Carlo simulations.
We conclude our investigation on the QCD equation of state (EoS) with 2+1 staggered flavors and one-link stout improvement. We extend our previous study [JHEP 0601:089 (2006)] by choosing even finer lattices. These new results [for details see arXiv:
1007.2580] support our earlier findings. Lattices with N_t=6,8 and 10 are used, and the continuum limit is approached by checking the results at N_t=12. A Symanzik improved gauge and a stout-link improved staggered fermion action is taken; the light and strange quark masses are set to their physical values. Various observables are calculated in the temperature (T) interval of 100 to 1000~MeV. We compare our data to the equation of state obtained by the hotQCD collaboration.
Finite temperature charmonium spectral functions in the pseudoscalar and vector channels are studied in lattice QCD with 2+1 flavours of dynamical Wilson quarks, on fine isotropic lattices (with a lattice spacing of 0.057 fm), with a non-physical pio
n mass of $m_{pi} approx$ 545 MeV. The highest temperature studied is approximately $1.4 T_c$. Up to this temperature no significant variation of the spectral function is seen in the pseudoscalar channel. The vector channel shows some temperature dependence, which seems to be consistent with a temperature dependent low frequency peak related to heavy quark transport, plus a temperature independent term at omega>0. These results are in accord with previous calculations using the quenched approximation.
We describe an implementation of the Rational Hybrid Monte Carlo (RHMC) algorithm for dynamical computations with two flavours of staggered quarks. We discuss several variants of the method, the performance and possible sources of error for each of t
hem, and we compare the performance and results to the inexact R algorithm.
We present the first results for the Kl3 form factor from simulations with 2+1 flavours of dynamical domain wall quarks. Combining our result, namely f_+(0)=0.964(5), with the latest experimental results for Kl3 decays leads to |V_{us}|=0.2249(14), r
educing the uncertaintity in this important parameter. For the O(p^6) term in the chiral expansion we obtain Delta f=-0.013(5).