TWQCD's dynamical DWF project
We present an overview of our project of simulation of unquenched lattice QCD with optimal domain-wall quarks, using a GPU cluster currently constituting of 16 units of Nvidia Tesla S1070 plus 64 graphic cards with Nvidia GTX285 (total 128 GPUs with 128 Teraflops peak), attaining sustained computing power of 15.36 Teraflops. The first production run in two-flavor QCD is on-going, using the Iwasaki gauge action on a set of lattices with sizes $ 16^3 \times (32,10,8,6,4) \times (16,32) $ at the lattice spacing $ a \sim 0.1$ fm, with eight sea quark masses down to $ m_\pi \simeq 200 $ MeV. We outline our simulation algorithm, and describe the present status of the production run. Preliminary results of pseudoscalar mass and decay constant are also presented.

title: The Dependence of Subhalo Abundance on Halo Concentration
abstract: Hierarchical structure formation implies that the number of subhalos within a dark matter halo depends not only on halo mass, but also on the formation history of the halo. This dependence on the formation history, which is highly correlated with halo concentration, can account for the super-Poissonian scatter in subhalo occupation at a fixed halo mass that has been previously measured in simulations. Here we propose a model to predict the subhalo abundance function for individual host halos that incorporates both halo mass and concentration. We combine results of cosmological simulations with a new suite of zoom-in simulations of Milky Way-mass halos to calibrate our model. We show that the model can successfully reproduce the mean and the scatter of subhalo occupation in these simulations. The implications of this correlation between subhalo abundance and halo concentration are further investigated. We also discuss cases in which inferences about halo properties can be affected if this correlation between subhalo abundance and halo concentration is ignored; in these cases, our model would give a more accurate inference. We propose that with future deep surveys, satellite occupation in the low-mass regime can be used to verify the existence of halo assembly bias.

title: GPU-Based Conjugate Gradient Solver for Lattice QCD with Domain-Wall Fermions
abstract: We present the first GPU-based conjugate gradient (CG) solver for lattice QCD with domain-wall fermions (DWF). It is well-known that CG is the most time-consuming part in the Hybrid Monte Carlo simulation of unquenched lattice QCD, which becomes even more computational demanding for lattice QCD with exact chiral symmetry. We have designed a CG solver for the general 5-dimensional DWF operator on NVIDIA CUDA architecture with mixed-precision, using the defect correction as well as the reliable updates algorithms. We optimize our computation by even-odd preconditioning in the 4D space-time lattice, plus several innovative techniques for CUDA kernels. For NVIDIA GeForce GTX 285/480, our CG solver attains 180/233 Gflops (sustained).

title: Topological susceptibility to the one-loop order in chiral perturbation theory
abstract: We derive the topological susceptibility to the one-loop order in the chiral effective theory of QCD, for an arbitrary number of flavors.

title: Topological susceptibility in two flavors lattice QCD with the optimal domain-wall fermion
abstract: We determine the topological susceptibility of the gauge configurations generated by lattice simulations using two flavors of optimal domain-wall fermion on the 16 3 × 32 lattice with length 16 in the fifth dimension, at the lattice spacing a ≃ 0.1 fm . Using the adaptive thick-restart Lanczos algorithm, we project the low-lying eigenmodes of the overlap Dirac operator, and obtain the topological charge of each configuration, for eight ensembles with pion masses in the range 220–550 MeV. From the topological charge, we compute the topological susceptibility and the second normalized cumulant. Our result of the topological susceptibility agrees with the sea-quark mass dependence predicted by the chiral perturbation theory and provides a determination of the chiral condensate, Σ MS ¯ ( 2 GeV ) = [ 259 ( 6 ) ( 7 ) MeV ] 3 , and the pion decay constant F π = 92 ( 12 ) ( 2 ) MeV .

title: SUSSING MERGER TREES: the influence of the halo finder
abstract: Merger tree codes are routinely used to follow the growth and merger of dark matter haloes in simulations of cosmic structure formation. Whereas in Srisawat et. al. we compared the trees built using a wide variety of such codes, here we study the influence of the underlying halo catalogue upon the resulting trees. We observe that the specifics of halo finding itself greatly influences the constructed merger trees. We find that the choices made to define the halo mass are of prime importance. For instance, amongst many potential options different finders select self-bound objects or spherical regions of defined overdensity, decide whether or not to include substructures within the mass returned and vary in their initial particle selection. The impact of these decisions is seen in tree length (the period of time a particularly halo can be traced back through the simulation), branching ratio (essentially the merger rate of subhaloes) and mass evolution. We therefore conclude that the choice of the underlying halo finder is more relevant to the process of building merger trees than the tree builder itself. We also report on some built-in features of specific merger tree codes that (sometimes) help to improve the quality of the merger trees produced.

title: Measurements of jet shape in 2.76 TeV PbPb collisions with CMS
abstract: The fragmentation properties of inclusive jets produced in PbPb collisions at a center of mass energy of s N N = 2.76 TeV are characterized by measuring differential jet shapes. The jet shape is defined as the fractional jet transverse momentum carried by charged particles emitted at a distance r from the jet axis, with jets reconstructed using the anti- k T clustering algorithm. The jet shapes are measured for reconstructed jet p T , jet > 100 GeV / c and charged hardons with p T > 1 GeV / c . Different background subtraction methods are employed to remove the contribution from the underlying events. The results obtained for PbPb collisions as a function of collision centrality are compared to reference distributions from pp data collected at the same collision energy. For the most central collisions, the jet shapes are found to be modified at large distance from the jet axis.

title: Inclusive $J/\psi$ production in $pp$ collisions at $\sqrt{s} = 2.76$ TeV
abstract: The ALICE Collaboration has measured inclusive J/psi production in pp collisions at a center of mass energy sqrt(s)=2.76 TeV at the LHC. The results presented in this Letter refer to the rapidity ranges |y|<0.9 and 2.5

title: CHARM QUARK ENERGY LOSS IN DENSE MATTER WITHIN THE LIGHT-CONE PATH INTEGRAL APPROACH
abstract: The energy loss of heavy quarks traversing color dense matter is calculated with an analytical formula derived within the light-cone path integral (LCPI) approach. We find that the energy loss mechanism is dominated by gluon radiation and induces a suppression pattern of charm quark different from the suppression of light quarks. We find also that this radiative energy-loss is proportional to L2 for energetic quarks, but changes to a L dependance with decreasing quark energy, where L is the length of the traversed medium.

title: High-pT particle correlations in pp collisions at LHC/ALICE*
abstract: Two-particle correlation triggered by high-pT particles allows us to study hard scattering phenomena when full jet reconstruction is challenging. An analysis of the first ALICE pp data where charged and neutral particles isolated or not are used as trigger particles is presented. The two-particle correlation between the trigger (t) and the associate (a) particles is studied as a function of the imbalance parameter and interpreted in terms of the jet fragmentation function.

title: Two particle correlations: A probe of the LHC QCD medium
abstract: The properties of γ-jet pairs emitted in heavy-ion collisions provide an accurate mean to perform a tomographic measurement of the medium created in the collision through the study of the medium modified jet properties. The idea is to measure the distribution of hadrons emitted on the opposite side of the direct photon. The feasibility of such measurements is studied by applying the approach on the simulation data, we have demonstrated that this method allows us to measure, with a good approximation, both the jet fragmentation and the back-to-back azimuthal alignment of the direct photon and the jet. Comparing these two observables measured in pp collisions with the ones measured in AA collisions reveals the modifications induced by the medium on the jet structure and consequently allows us to infer the medium properties. In this contribution, we discuss a first attempt of such measurements applied to real proton-proton data from the ALICE experiment.