Hardware and software status of QCDOC

QCDOC is a massively parallel supercomputer whose processing nodes are based on an application-specific integrated circuit (ASIC). This ASIC was custom-designed so that crucial lattice QCD kernels achieve an overall sustained performance of 50% on machines with several 10,000 nodes. This strong scalability, together with low power consumption and a price/performance ratio of $1 per sustained MFlops, enable QCDOC to attack the most demanding lattice QCD problems. The first ASICs became available in June of 2003, and the testing performed so far has shown all systems functioning according to specification. We review the hardware and software status of QCDOC and present performance figures obtained in real hardware as well as in simulation.Comment: Lattice2003(machine), 6 pages, 5 figure

The Fall of the Quasar Population

We derive quantitative predictions of the optical and X-ray luminosity functions (LF) for QSs in the redshift range $z<3$. Based on BH paradigm, we investigate how the accretion is controlled by the surrounding structures, as these grow hierarchically. We argue that for $z < 3$ efficient black hole fueling is triggered by the encounters of a gas-rich host with its companions in a group. The dispersion of the dynamical parameters in the encounters produces a double power-law LF. Strong luminosity evolution (LE) is produced as these encounters deplete the gas supply in the host; an additional, milder density evolution obtains since the interactions become progressively rarer as the groups grow richer but less dense. From the agreement with the optical and the X-ray data, we conclude that the evolution of the bright quasars is articulated in two ways. Earlier than $z~3$ the gas-rich protogalaxies grow by merging, which also induces parallel growth of central holes accreting at Eddington rates. In the later era of group assemblage the host encounters with companions drive onto already existing holes further but meager accretion; these consume the gas in the hosts, and cause supply-limited emissions which are intermittent, go progressively sub-Eddington and peter out. Then other fueling processes come to the foreground; we discuss the faint emissions, especially noticeable in X-rays, which are expected when hosts in the field cannibalize satellite galaxies with their meager gas contents.Comment: 12 pages Latex + 3 EPS figures, ApJ in press, we have corrected the previous printing problems with the style \ca

Lattice determination of the K→(ππ)I=2K \to (\pi\pi)_{I=2} Decay Amplitude A2A_2

We describe the computation of the amplitude A_2 for a kaon to decay into two pions with isospin I=2. The results presented in the letter Phys.Rev.Lett. 108 (2012) 141601 from an analysis of 63 gluon configurations are updated to 146 configurations giving Re$A_2=1.381(46)_{\textrm{stat}}(258)_{\textrm{syst}} 10^{-8}$ GeV and Im$A_2=-6.54(46)_{\textrm{stat}}(120)_{\textrm{syst}}10^{-13}$ GeV. Re$A_2$ is in good agreement with the experimental result, whereas the value of Im$A_2$ was hitherto unknown. We are also working towards a direct computation of the $K\to(\pi\pi)_{I=0}$ amplitude $A_0$ but, within the standard model, our result for Im$A_2$ can be combined with the experimental results for Re$A_0$, Re$A_2$ and $\epsilon^\prime/\epsilon$ to give Im$A_0/$Re$A_0= -1.61(28)\times 10^{-4}$ . Our result for Im\,$A_2$ implies that the electroweak penguin (EWP) contribution to $\epsilon^\prime/\epsilon$ is Re$(\epsilon^\prime/\epsilon)_{\mathrm{EWP}} = -(6.25 \pm 0.44_{\textrm{stat}} \pm 1.19_{\textrm{syst}}) \times 10^{-4}$.Comment: 59 pages, 11 figure

Opening the Rome-Southampton window for operator mixing matrices

We show that the running of operators which mix under renormalization can be computed fully non-perturbatively as a product of continuum step scaling matrices. These step scaling matrices are obtained by taking the "ratio" of Z matrices computed at different energies in an RI-MOM type scheme for which twisted boundary conditions are an essential ingredient. Our method allows us to relax the bounds of the Rome-Southampton window. We also explain why such a method is important in view of the light quark physics program of the RBC-UKQCD collaborations. To illustrate our method, using n_f=2+1 domain-wall fermions, we compute the non-perturbative running matrix of four-quark operators needed in K->pipi decay and neutral kaon mixing. Our results are then compared to perturbation theory.Comment: 8 pages, 7 figures. v2: PRD version, minor changes and few references adde

QSO clustering and the AAT 2dF redshift survey

We review previous results on the clustering and environments of QSOs. We show that the correlation length for QSOs derived from existing surveys is r~5/h Mpc, similar to the observed correlation length for field galaxies at the present epoch. The galaxy environment for z<1 radio-quiet QSOs is also consistent with field galaxies. The evolution of the QSO correlation length with redshift is currently uncertain, largely due to the small numbers of QSOs (~2000) in surveys suitable for clustering analysis. We report on intial progress with the AAT 2dF QSO redshift survey, which, once completed will comprise almost 30000 QSOs. With over 1000 QSOs already observed, it is already the largest single homogeneous QSO survey. We discuss prospects for deriving limits on cosmological parameters from this survey, and on the evolution of large-scale structure in the Universe.Comment: Invited talk at RS meeting on 'Large Scale Structure in the Universe' held at the Royal Society on 25-26 March 1998 14 pages, 11 figre

The NINJA-2 catalog of hybrid post-Newtonian/numerical-relativity waveforms for non-precessing black-hole binaries

The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational wave data analysis communities. The purpose of NINJA is to study the sensitivity of existing gravitational-wave search and parameter-estimation algorithms using numerically generated waveforms, and to foster closer collaboration between the numerical relativity and data analysis communities. The first NINJA project used only a small number of injections of short numerical-relativity waveforms, which limited its ability to draw quantitative conclusions. The goal of the NINJA-2 project is to overcome these limitations with long post-Newtonian - numerical relativity hybrid waveforms, large numbers of injections, and the use of real detector data. We report on the submission requirements for the NINJA-2 project and the construction of the waveform catalog. Eight numerical relativity groups have contributed 63 hybrid waveforms consisting of a numerical portion modelling the late inspiral, merger, and ringdown stitched to a post-Newtonian portion modelling the early inspiral. We summarize the techniques used by each group in constructing their submissions. We also report on the procedures used to validate these submissions, including examination in the time and frequency domains and comparisons of waveforms from different groups against each other. These procedures have so far considered only the $(ell,m)=(2,2)$ mode. Based on these studies we judge that the hybrid waveforms are suitable for NINJA-2 studies. We note some of the plans for these investigations

Standard-model prediction for direct CP violation in K→ππK\to\pi\pi decay

We report the first lattice QCD calculation of the complex kaon decay amplitude $A_0$ with physical kinematics, using a $32^3\times 64$ lattice volume and a single lattice spacing $a$, with $1/a= 1.3784(68)$ GeV. We find Re$(A_0) = 4.66(1.00)(1.26) \times 10^{-7}$ GeV and Im$(A_0) = -1.90(1.23)(1.08) \times 10^{-11}$ GeV, where the first error is statistical and the second systematic. The first value is in approximate agreement with the experimental result: Re$(A_0) = 3.3201(18) \times 10^{-7}$ GeV while the second can be used to compute the direct CP violating ratio Re$(\varepsilon'/\varepsilon)=1.38(5.15)(4.59)\times 10^{-4}$, which is $2.1\sigma$ below the experimental value $16.6(2.3)\times 10^{-4}$. The real part of $A_0$ is CP conserving and serves as a test of our method while the result for Re$(\varepsilon'/\varepsilon)$ provides a new test of the standard-model theory of CP violation, one which can be made more accurate with increasing computer capability.Comment: 9 pages, 3 figures. Updated to match published versio

Circular economy: Questions for responsible minerals, additive manufacturing and recycling of metals

© 2014 by the authors; licensee MDPI, Basel, Switzerland. The concept of the circular economy proposes new patterns of production, consumption and use, based on circular flows of resources. Under a scenario where there is a global shift towards the circular economy, this paper discusses the advent of two parallel and yet-to-be-connected trends for Australia, namely: (i) responsible minerals supply chains and (ii) additive manufacturing, also known as 3D production systems. Acknowledging the current context for waste management, the paper explores future interlinked questions which arise in the circular economy for responsible supply chains, additive manufacturing, and metals recycling. For example, where do mined and recycled resources fit in responsible supply chains as inputs to responsible production? What is required to ensure 3D production systems are resource efficient? How could more distributed models of production, enabled by additive manufacturing, change the geographical scale at which it is economic or desirable to close the loop? Examples are given to highlight the need for an integrated research agenda to address these questions and to foster Australian opportunities in the circular economy