Source author record

Thorsten Kurth

Thorsten Kurth appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

14works
12topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2022arXiv

FourCastNet: A Global Data-driven High-resolution Weather Model using Adaptive Fourier Neural Operators

FourCastNet, short for Fourier Forecasting Neural Network, is a global data-driven weather forecasting model that provides accurate short to medium-range global predictions at $0.25^{\circ}$ resolution. FourCastNet accurately forecasts high-resolution, fast-timescale variables such as the surface wind speed, precipitation, and atmospheric water vapor. It has important implications for planning wind energy resources, predicting extreme weather events such as tropical cyclones, extra-tropical cyclones, and atmospheric rivers. FourCastNet matches the forecasting accuracy of the ECMWF Integrated Forecasting System (IFS), a state-of-the-art Numerical Weather Prediction (NWP) model, at short lead times for large-scale variables, while outperforming IFS for variables with complex fine-scale structure, including precipitation. FourCastNet generates a week-long forecast in less than 2 seconds, orders of magnitude faster than IFS. The speed of FourCastNet enables the creation of rapid and inexpensive large-ensemble forecasts with thousands of ensemble-members for improving probabilistic forecasting. We discuss how data-driven deep learning models such as FourCastNet are a valuable addition to the meteorology toolkit to aid and augment NWP models.

preprint2022arXiv

FourCastNet: Accelerating Global High-Resolution Weather Forecasting using Adaptive Fourier Neural Operators

Extreme weather amplified by climate change is causing increasingly devastating impacts across the globe. The current use of physics-based numerical weather prediction (NWP) limits accuracy due to high computational cost and strict time-to-solution limits. We report that a data-driven deep learning Earth system emulator, FourCastNet, can predict global weather and generate medium-range forecasts five orders-of-magnitude faster than NWP while approaching state-of-the-art accuracy. FourCast-Net is optimized and scales efficiently on three supercomputing systems: Selene, Perlmutter, and JUWELS Booster up to 3,808 NVIDIA A100 GPUs, attaining 140.8 petaFLOPS in mixed precision (11.9%of peak at that scale). The time-to-solution for training FourCastNet measured on JUWELS Booster on 3,072GPUs is 67.4minutes, resulting in an 80,000times faster time-to-solution relative to state-of-the-art NWP, in inference. FourCastNet produces accurate instantaneous weather predictions for a week in advance, enables enormous ensembles that better capture weather extremes, and supports higher global forecast resolutions.

preprint2020arXiv

$F_K / F_π$ from Möbius domain-wall fermions solved on gradient-flowed HISQ ensembles

We report the results of a lattice quantum chromodynamics calculation of $F_K/F_π$ using Möbius domain-wall fermions computed on gradient-flowed $N_f=2+1+1$ highly-improved staggered quark (HISQ) ensembles. The calculation is performed with five values of the pion mass ranging from $130 \lesssim m_π\lesssim 400$ MeV, four lattice spacings of $a\sim 0.15, 0.12, 0.09$ and $0.06$ fm and multiple values of the lattice volume. The interpolation/extrapolation to the physical pion and kaon mass point, the continuum, and infinite volume limits are performed with a variety of different extrapolation functions utilizing both the relevant mixed-action effective field theory expressions as well as discretization-enhanced continuum chiral perturbation theory formulas. We find that the $a\sim0.06$ fm ensemble is helpful, but not necessary to achieve a subpercent determination of $F_K/F_π$. We also include an estimate of the strong isospin breaking corrections and arrive at a final result of $F_{K^\pm}/F_{π^\pm} = 1.1942(45)$ with all sources of statistical and systematic uncertainty included. This is consistent with the Flavour Lattice Averaging Group average value, providing an important benchmark for our lattice action. Combining our result with experimental measurements of the pion and kaon leptonic decays leads to a determination of $|V_{us}|/|V_{ud}| = 0.2311(10)$.

preprint2019arXiv

Lattice QCD Determination of $g_A$

The nucleon axial coupling, $g_A$, is a fundamental property of protons and neutrons, dictating the strength with which the weak axial current of the Standard Model couples to nucleons, and hence, the lifetime of a free neutron. The prominence of $g_A$ in nuclear physics has made it a benchmark quantity with which to calibrate lattice QCD calculations of nucleon structure and more complex calculations of electroweak matrix elements in one and few nucleon systems. There were a number of significant challenges in determining $g_A$, notably the notorious exponentially-bad signal-to-noise problem and the requirement for hundreds of thousands of stochastic samples, that rendered this goal more difficult to obtain than originally thought. I will describe the use of an unconventional computation method, coupled with "ludicrously'" fast GPU code, access to publicly available lattice QCD configurations from MILC and access to leadership computing that have allowed these challenges to be overcome resulting in a determination of $g_A$ with 1% precision and all sources of systematic uncertainty controlled. I will discuss the implications of these results for the convergence of $SU(2)$ Chiral Perturbation theory for nucleons, as well as prospects for further improvements to $g_A$ (sub-percent precision, for which we have preliminary results) which is part of a more comprehensive application of lattice QCD to nuclear physics. This is particularly exciting in light of the new CORAL supercomputers coming online, Sierra and Summit, for which our lattice QCD codes achieve a machine-to-machine speed up over Titan of an order of magnitude.

preprint2018arXiv

Simulating the weak death of the neutron in a femtoscale universe with near-Exascale computing

The fundamental particle theory called Quantum Chromodynamics (QCD) dictates everything about protons and neutrons, from their intrinsic properties to interactions that bind them into atomic nuclei. Quantities that cannot be fully resolved through experiment, such as the neutron lifetime (whose precise value is important for the existence of light-atomic elements that make the sun shine and life possible), may be understood through numerical solutions to QCD. We directly solve QCD using Lattice Gauge Theory and calculate nuclear observables such as neutron lifetime. We have developed an improved algorithm that exponentially decreases the time-to solution and applied it on the new CORAL supercomputers, Sierra and Summit. We use run-time autotuning to distribute GPU resources, achieving 20% performance at low node count. We also developed optimal application mapping through a job manager, which allows CPU and GPU jobs to be interleaved, yielding 15% of peak performance when deployed across large fractions of CORAL.

preprint2016arXiv

Neutrinoless double beta decay from lattice QCD

While the discovery of non-zero neutrino masses is one of the most important accomplishments by physicists in the past century, it is still unknown how and in what form these masses arise. Lepton number-violating neutrinoless double beta decay is a natural consequence of Majorana neutrinos and many BSM theories, and many experimental efforts are involved in the search for these processes. Understanding how neutrinoless double beta decay would manifest in nuclear environments is key for understanding any observed signals. In these proceedings we present an overview of a set of one- and two-body matrix elements relevant for experimental searches for neutrinoless double beta decay, describe the role of lattice QCD calculations, and present preliminary lattice QCD results.

preprint2016arXiv

Two-Nucleon Higher Partial-Wave Scattering from Lattice QCD

We present a determination of nucleon-nucleon scattering phase shifts for l >= 0. The S, P, D and F phase shifts for both the spin-triplet and spin-singlet channels are computed with lattice Quantum ChromoDynamics. For l > 0, this is the first lattice QCD calculation using the Luscher finite-volume formalism. This required the design and implementation of novel lattice methods involving displaced sources and momentum-space cubic sinks. To demonstrate the utility of our approach, the calculations were performed in the SU(3)-flavor limit where the light quark masses have been tuned to the physical strange quark mass, corresponding to m_pi = m_K ~ 800 MeV. In this work, we have assumed that only the lowest partial waves contribute to each channel, ignoring the unphysical partial wave mixing that arises within the finite-volume formalism. This assumption is only valid for sufficiently low energies; we present evidence that it holds for our study using two different channels. Two spatial volumes of V ~ (3.5 fm)^3 and V ~ (4.6 fm)^3 were used. The finite-volume spectrum is extracted from the exponential falloff of the correlation functions. Said spectrum is mapped onto the infinite volume phase shifts using the generalization of the Luscher formalism for two-nucleon systems.

preprint2015arXiv

High-Performance I/O: HDF5 for Lattice QCD

Practitioners of lattice QCD/QFT have been some of the primary pioneer users of the state-of-the-art high-performance-computing systems, and contribute towards the stress tests of such new machines as soon as they become available. As with all aspects of high-performance-computing, I/O is becoming an increasingly specialized component of these systems. In order to take advantage of the latest available high-performance I/O infrastructure, to ensure reliability and backwards compatibility of data files, and to help unify the data structures used in lattice codes, we have incorporated parallel HDF5 I/O into the SciDAC supported USQCD software stack. Here we present the design and implementation of this I/O framework. Our HDF5 implementation outperforms optimized QIO at the 10-20% level and leaves room for further improvement by utilizing appropriate dataset chunking.

preprint2015arXiv

Nuclear Parity Violation from Lattice QCD

The electroweak interaction at the level of quarks and gluons are well understood from precision measurements in high energy collider experiments. Relating these fundamental parameters to Hadronic Parity Violation in nuclei however remains an outstanding theoretical challenge. One of the most interesting observables in this respect is the parity violating hadronic neutral current: it is hard to measure in collider experiments and is thus the least constrained observable of the Standard Model. Precision measurements of parity violating transitions in nuclei can help to improve these constraints. In these systems however, the weak interaction is masked by effects of the seven orders of magnitude stronger non-perturbative strong interaction. Therefore, in order to relate experimental measurements of the parity violating pion-nucleon couplings to the fundamental Lagrangian of the SM, these non-perturbative effects have to be well understood. In this paper, we are going to present a Lattice QCD approach for computing the $ΔI{=}2$ parity violating matrix element in proton proton scattering. This process does not involve disconnected diagrams in the isospin symmetric limit and is thus a perfect testbed for studying the feasibility of the more involved calculation of the parity violating pion-nucleon coupling.

preprint2015arXiv

Two-nucleon scattering in multiple partial waves

We determine scattering phase shifts for S,P,D, and F partial wave channels in two-nucleon systems using lattice QCD methods. We use a generalization of Luscher's finite volume method to determine infinite volume phase shifts from a set of finite volume ground- and excited-state energy levels on two volumes, V=(3.4 fm)^3 and V=(4.5 fm)^3. The calculations are performed in the SU(3)-flavor limit, corresponding to a pion mass of approximately 800 MeV. From the energy dependence of the phase shifts we are able to extract scattering parameters corresponding to an effective range expansion.

preprint2013arXiv

Lattice QCD at the physical point meets SU(2) chiral perturbation theory

We perform a detailed, fully-correlated study of the chiral behavior of the pion mass and decay constant, based on 2+1 flavor lattice QCD simulations. These calculations are implemented using tree-level, O(a)-improved Wilson fermions, at four values of the lattice spacing down to 0.054 fm and all the way down to below the physical value of the pion mass. They allow a sharp comparison with the predictions of SU(2) chiral perturbation theory (χPT) and a determination of some of its low energy constants. In particular, we systematically explore the range of applicability of NLO SU(2) χPT in two different expansions: the first in quark mass (x-expansion), and the second in pion mass (ξ-expansion). We find that these expansions begin showing signs of failure around M_π=300 MeV for the typical percent-level precision of our N_f=2+1 lattice results. We further determine the LO low energy constants (LECs), F=88.0 \pm 1.3\pm 0.3 and B^\msbar(2 GeV)=2.58 \pm 0.07 \pm 0.02 GeV, and the related quark condensate, Σ^\msbar(2 GeV)=(271\pm 4\pm 1 MeV)^3, as well as the NLO ones, l_3=2.5 \pm 0.5 \pm 0.4 and l_4=3.8 \pm 0.4 \pm 0.2, with fully controlled uncertainties. We also explore the NNLO expansions and the values of NNLO LECs. In addition, we show that the lattice results favor the presence of chiral logarithms. We further demonstrate how the absence of lattice results with pion masses below 200 MeV can lead to misleading results and conclusions. Our calculations allow a fully controlled, ab initio determination of the pion decay constant with a total 1% error, which is in excellent agreement with experiment.

preprint2012arXiv

Systematic errors in partially-quenched QCD plus QED lattice simulations

At the precision reached in current lattice QCD calculations, electromagnetic effects are becoming numerically relevant. Here, electromagnetic effects are included by superimposing $\mathrm{U}(1)$ degrees of freedom on $N_f = 2+1$ QCD configurations from the Budapest-Marseille-Wuppertal Collaboration. We present preliminary results for the electromagnetic corrections to light pseudoscalars mesons masses and discuss some of the associated systematic errors.

preprint2010arXiv

Scaling study for 2 HEX smeared fermions: hadron and quark masses

The goal of this study is to investigate the scaling behaviour of our 2 HEX action. For this purpose, we compute the $N_f=3$ spectrum and compare the results to our 6 EXP action. We find a large scaling window up to $\sim 0.15\,\mathrm{fm}$ along with small scaling corrections at the 2%-level and full compatibility with our previous study. As a second important observable to be tested for scaling, we chose the non-perturbatively renormalized quenched strange quark mass. Here we find a fairly flat scaling with a broad scaling range up to $\simeq 0.15\,\mathrm{fm}$ and perfect agreement with the literature.