Equilibrium Thermodynamics of Lattice QCD

Lattice QCD allows us to simulate QCD at non-zero temperature and/or densities. Such equilibrium thermodynamics calculations are relevant to the physics of relativistic heavy-ion collisions. I give a brief review of the field with emphasis on our work.Comment: 15 pages, 9 figures. Talk presented at SCGT06, Nagoya, Japan. Version 2 includes minor modifications to reference work not covered in version

The distribution of Mahler's measures of reciprocal polynomials

We study the distribution of Mahler's measures of reciprocal polynomials with complex coefficients and bounded even degree. We discover that the distribution function associated to Mahler's measure restricted to monic reciprocal polynomials is a reciprocal (or anti-reciprocal) Laurent polynomial on [1,\infty) and identically zero on [0,1). Moreover, the coefficients of this Laurent polynomial are rational numbers times a power of \pi. We are led to this discovery by the computation of the Mellin transform of the distribution function. This Mellin transform is an even (or odd) rational function with poles at small integers and residues that are rational numbers times a power of \pi. We also use this Mellin transform to show that the volume of the set of reciprocal polynomials with complex coefficients, bounded degree and Mahler's measure less than or equal to one is a rational number times a power of \pi.Comment: 13 pages. To be published in Int. J. Math. Math. Sc

Fatigue of friction stir welded 2024-T351 aluminium alloy

Fatigue failure characteristics of friction stir welds in 13mm gauge 2024-T351 plate have been assessed. Failure occurred from either the weld region (nugget/flow arm) or from the material immediately surrounding the weld. Fatigue failure from the surrounding material was essentially conventional, initiating from large S-phase intermetallic particles and growing in a macroscopic mode I manner. Corresponding fatigue lives were seen to be comparable to parent plate and results previously reported for similar welds in thinner plate. Failure over the weld region was identified with discontinuities in the macroscopic flow pattern of the weld flow arm. Subsequent crack growth showed pronounced macroscopic crack deflection around the ‘onion ring’ structure of the weld nugget. The bands making up the onion rings were identified with variations in local hardness levels, consistent with a mechanical contribution to the crack deflection process

The Ginibre ensemble of real random matrices and its scaling limits

We give a closed form for the correlation functions of ensembles of a class of asymmetric real matrices in terms of the Pfaffian of an antisymmetric matrix formed from a $2 \times 2$ matrix kernel associated to the ensemble. We apply this result to the real Ginibre ensemble and compute the bulk and edge scaling limits of its correlation functions as the size of the matrices becomes large.Comment: 47 pages, 8 figure

Extremal laws for the real Ginibre ensemble

The real Ginibre ensemble refers to the family of $n\times n$ matrices in which each entry is an independent Gaussian random variable of mean zero and variance one. Our main result is that the appropriately scaled spectral radius converges in law to a Gumbel distribution as $n\rightarrow\infty$. This fact has been known to hold in the complex and quaternion analogues of the ensemble for some time, with simpler proofs. Along the way we establish a new form for the limit law of the largest real eigenvalue.Comment: Published in at http://dx.doi.org/10.1214/13-AAP958 the Annals of Applied Probability (http://www.imstat.org/aap/) by the Institute of Mathematical Statistics (http://www.imstat.org

Complex Langevin Simulations of QCD at Finite Density -- Progress Report

We simulate lattice QCD at finite quark-number chemical potential to study nuclear matter, using the complex Langevin equation (CLE). The CLE is used because the fermion determinant is complex so that standard methods relying on importance sampling fail. Adaptive methods and gauge-cooling are used to prevent runaway solutions. Even then, the CLE is not guaranteed to give correct results. We are therefore performing extensive testing to determine under what, if any, conditions we can achieve reliable results. Our earlier simulations at $\beta=6/g^2=5.6$, $m=0.025$ on a $12^4$ lattice reproduced the expected phase structure but failed in the details. Our current simulations at $\beta=5.7$ on a $16^4$ lattice fail in similar ways while showing some improvement. We are therefore moving to even weaker couplings to see if the CLE might produce the correct results in the continuum (weak-coupling) limit, or, if it still fails, whether it might reproduce the results of the phase-quenched theory. We also discuss action (and other dynamics) modifications which might improve the performance of the CLE.Comment: Talk presented at Lattice 2017, Granada, Spain and submitted to proceedings. 8 pages, 4 figure