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

High Energy Physics from High Performance Computing

We discuss Quantum Chromodynamics calculations using the lattice regulator. The theory of the strong force is a cornerstone of the Standard Model of particle physics. We present USQCD collaboration results obtained on Argonne National Lab's Intrepid supercomputer that deepen our understanding of these fundamental theories of Nature and provide critical support to frontier particle physics experiments and phenomenology.Comment: Proceedings of invited plenary talk given at SciDAC 2009, San Diego, June 14-18, 2009, on behalf of the USQCD collaboratio

The K→(ππ)I=2K\to(\pi\pi)_{I=2} Decay Amplitude from Lattice QCD

We report on the first realistic \emph{ab initio} calculation of a hadronic weak decay, that of the amplitude $A_2$ for a kaon to decay into two \pi-mesons with isospin 2. We find Re$A_2=(1.436\pm 0.063_{\textrm{stat}}\pm 0.258_{\textrm{syst}})\,10^{-8}\,\textrm{GeV}$ in good agreement with the experimental result and for the hitherto unknown imaginary part we find {Im}$\,A_2=-(6.83 \pm 0.51_{\textrm{stat}} \pm 1.30_{\textrm{syst}})\,10^{-13}\,{\rm GeV}$. Moreover combining our result for Im\,$A_2$ with experimental values of Re\,$A_2$, Re\,$A_0$ and $\epsilon^\prime/\epsilon$, we obtain the following value for the unknown ratio Im\,$A_0$/Re\,$A_0$ within the Standard Model: $\mathrm{Im}\,A_0/\mathrm{Re}\,A_0=-1.63(19)_{\mathrm{stat}}(20)_{\mathrm{syst}}\times10^{-4}$. One consequence of these results is that the contribution from Im\,$A_2$ to the direct CP violation parameter $\epsilon^{\prime}$ (the so-called Electroweak Penguin, EWP, contribution) is Re$(\epsilon^\prime/\epsilon)_{\mathrm{EWP}} = -(6.52 \pm 0.49_{\textrm{stat}} \pm 1.24_{\textrm{syst}}) \times 10^{-4}$. We explain why this calculation of $A_2$ represents a major milestone for lattice QCD and discuss the exciting prospects for a full quantitative understanding of CP-violation in kaon decays.Comment: 5 pages, 1 figur

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

Atmospheric drivers of melt on Larsen C Ice Shelf: Surface energy budget regimes and the impact of foehn

Recent ice shelf retreat on the east coast of the Antarctic Peninsula has been principally attributed to atmospherically driven melt. However, previous studies on the largest of these ice shelves—Larsen C—have struggled to reconcile atmospheric forcing with observed melt. This study provides the first comprehensive quantification and explanation of the atmospheric drivers of melt across Larsen C, using 31-months' worth of observations from Cabinet Inlet, a 6-month, high-resolution atmospheric model simulation and a novel approach to ascertain the surface energy budget (SEB) regime. The dominant meteorological controls on melt are shown to be the occurrence, strength, and warmth of mountain winds called foehn. At Cabinet Inlet, foehn occurs 15% of the time and causes 45% of melt. The primary effect of foehn on the SEB is elevated turbulent heat fluxes. Under typical, warm foehn conditions, this means elevated surface heating and melting, the intensity of which increases as foehn wind speed increases. Less commonly—due to cooler-than-normal foehn winds and/or radiatively warmed ice—the relationship between wind speed and net surface heat flux reverses. This explains the seemingly contradictory results of previous studies. In the model, spatial variability in cumulative melt across Larsen C is largely explained by foehn, with melt maxima in inlets reflecting maxima in foehn wind strength. However, most accumulated melt (58%) occurs due to solar radiation in the absence of foehn. A broad north-south gradient in melt is explained by the combined influence of foehn and non-foehn conditions

Decreasing intensity of open-ocean convection in the Greenland and Iceland seas

The air–sea transfer of heat and fresh water plays a critical role in the global climate system. This is particularly true for the Greenland and Iceland seas, where these fluxes drive ocean convection that contributes to Denmark Strait overflow water, the densest component of the lower limb of the Atlantic Meridional Overturning Circulation (AMOC). Here we show that the wintertime retreat of sea ice in the region, combined with different rates of warming for the atmosphere and sea surface of the Greenland and Iceland seas, has resulted in statistically significant reductions of approximately 20% in the magnitude of the winter air–sea heat fluxes since 1979. We also show that modes of climate variability other than the North Atlantic Oscillation (NAO) are required to fully characterize the regional air–sea interaction. Mixed-layer model simulations imply that further decreases in atmospheric forcing will exceed a threshold for the Greenland Sea whereby convection will become depth limited, reducing the ventilation of mid-depth waters in the Nordic seas. In the Iceland Sea, further reductions have the potential to decrease the supply of the densest overflow waters to the AMOC

Foehn jets over the Larsen C Ice Shelf, Antarctica

Previously unknown foehn jets have been identified to the east of the Antarctic Peninsula (AP) above the Larsen C Ice Shelf. These jets have major implications for the east coast of the AP, a region of rapid climatic warming and where two large sections of ice shelf have collapsed in recent years. During three foehn events across the AP, leeside warming and drying is seen in new aircraft observations and simulated well by the Met Office Unified Model (MetUM) at ∼1.5 km grid spacing. In case A, weak southwesterly flow and an elevated upwind inversion characterise a highly nonlinear flow regime with upwind flow blocking. In case C strong northwesterly winds characterise a relatively linear case with little upwind flow blocking. Case B resides somewhere between the two in flow regime linearity. The foehn jets – apparent in aircraft observations where available and MetUM simulations of all three cases – are mesoscale features (up to 60 km in width) originating from the mouths of leeside inlets. Through back trajectory analysis they are identified as a type of gap flow. In cases A and B the jets are distinct, being strongly accelerated relative to the background flow, and confined to low levels above the Larsen C Ice Shelf. They resemble the ‘shallow foehn’ of the Alps. Case C resembles a case of ‘deep foehn’, with the jets less distinct. The foehn jets are considerably cooler and moister relative to adjacent regions of calmer foehn air. This is due to a dampened foehn effect in the jet regions: in case A the jets have lower upwind source regions, and in the more linear case C there is less diabatic warming and precipitation along jet trajectories due to the reduced orographic uplift across the mountain passes