ASTROPHYSICS

Astrophysics

Principal Investigator: Prof. Dr. Marcus Brüggen , University of Hamburg, Hamburg, Germany

HPC Platform used: JUWELS CPU and JUWELS BOOSTER at JSC

Local Project ID: nonequioutflows and SHOCKCLOUD

Turbulence is a ubiquitous phenomenon that affects everything ranging from blood flow in our arteries, via aircraft to processes that form stars such as our Sun. In particular, turbulence that moves faster than the speed of sound, so-called supersonic turbulence, is important in many astrophysical settings, for example in giant molecular clouds that are the birth places of stars and that are scattered throughout galaxies. However, many properties of supersonic turbulence are poorly understood.

Astrophysics

Principal Investigator: Dr. Thales Gutcke , University of Hawaii, Hilo, United States

HPC Platform used: SuperMUC-NG at LRZ

Local Project ID: pn73we

This project targets various central questions in modern astrophysics, including "What is the nature of dark matter?", "How do galaxies form and evolve?" and "Do we understand the extremes of the universe?". Dwarf galaxies provide a natural laboratory for confronting these questions as we will explain below. The formation of dwarf galaxies tracks an extreme situation in various ways. Dwarf galaxies are assumed to be the very first type of galaxy to form in the earliest Universe.

The aim of this project is to create a highly accurate cosmological, hydrodynamical simulation model that can produce extremely realistic representations of dwarf galaxies right into the centers, where dark matter models can be tested.

Astrophysics

Principal Investigator: Dr. Salvatore Cielo , Bayrische Akademie der Wissenschaften, Garching, Germany

HPC Platform used: SuperMUC-NG of LRZ

Local Project ID: pn49xu

Numerical sciences are experiencing a renaissance thanks to the spread of heterogeneous computing. The SYCL open standard unlocks GPGPUs, accelerators, multicore and vector CPUs, and advanced compiler features and technologies (LLVM, JIT), while offering intuitive C++ APIs for work-sharing and scheduling. The project allowed for the kick-off of DPEcho (short for Data-Parallel ECHO), a SYCL+MPI porting of the General-Relativity-Magneto-Hydrodynamic (GRMHD) OpenMP+MPI code ECHO, used to model instabilities, turbulence, propagation of waves, stellar winds and magnetospheres, and astrophysical processes around Black Holes, in Cartesian or any coded GR metric.

 

Astrophysics

Principal Investigator: Prof. Dr. Rainer Grauer , Institute for theoretical physics I, Ruhr-University Bochum, Bochum, Germany

HPC Platform used: JUWELS CPU of JSC

Local Project ID: SpecDyn

Currently, the DRESDYN experiment is under construction. It aims to find and understand the mechanisms behind the excitation of a magnetic field due to the precessing motion of a conducting fluid. Such a precessing fluid is, for example, inside the core of the Earth, which is believed to sustain the Earth's magnetic field. To study the flow fields as well as the magnetic fields and to predict optimal parameter regimes for the DRESDYN experiment, numerical simulations are performed on JUWELS CPU using the code SpecDyn.

Astrophysics

Principal Investigator: Prof. Dr. Stefanie Walch-Gassner , Physikalisches Institut, Universität zu Köln, Cologne, Germany

HPC Platform used: JUWELS CPU of JSC

Local Project ID: COSMOS

High-mass star formation is a highly complex and dynamic process involving a large number of physical mechanisms. To better interpret real-world star-forming regions, simulations of collapsing clouds are used. These simulations produce the distribution of matter within star-forming regions, taking into account the effects of gravitational contraction as well as the feedback from massive stars. These simulations are then used to compute synthetic telescope images which may be compared to observations made with instruments like the Atacama Large Millimeter Array (ALMA).

Astrophysics

Principal Investigator: Dr. Stefan Gottlöber , Leibniz-Institut für Astrophysik Potsdam, Germany

HPC Platform used: JUWELS of JSC

Local Project ID: chpo22

Near field cosmology is the theoretical and observational study of our neighbourhood in the universe. This is the main topic of focus for the CLUES collaboration (www.clues-project.org). Our neighbourhood is the best observed part of the universe where also tiny dwarf galaxies can be studied. The properties of these dwarfs reflect their early formation history and state of the universe at the Cosmic Dawn, when the first stars and galaxies formed. Studying them sheds new light on these times, known as Cosmic Dawn and Epoch of Reionisation.

Astrophysics

Principal Investigator: Prof. Dr. Wolfgang Hillebrandt , Max Planck Gesellschaft, Max-Planck-Institut für Astrophysik, Garching, Germany

HPC Platform used: JUQUEEN, JURECAVIS, JUROPA and JUWELS CPU of JSC

Local Project ID: CHMU14

In project CHMU14, challenging three-dimensional simulations of thermonuclear explosions of white dwarf stars near the Chandrasekhar-mass limit were conducted. These were followed by radiative transfer simulations that allow to predict observables. A comparison with astronomical data shows that such models can explain the subclass of Type Iax supernovae.

Astrophysics

Principal Investigator: James R. Beattie and Prof. Christoph Federrath , Research School of Astronomy and Astrophysics, Australian National University

HPC Platform used: SuperMUC at LRZ

Local Project ID: pn73fi

Supersonic, magnetised turbulence is ubiquitous in the interstellar medium of galaxies. Unlike incompressible turbulence, supersonic turbulence is not scale-free. The scale that marks the transition from supersonic to subsonic turbulence is the so-called sonic scale, which in the context of star formation may define the critical value for which regions inside of molecular gas clouds collapse under their own gravity to form stars.

Astrophysics

Principal Investigator: Dr. Denis Wittor , University of Hamburg

HPC Platform used: JUWELS at JSC

Local Project ID: hhh44

A team of researchers led by Dr. Denis Wittor at the University of Hamburg has been leveraging the high-performance power of the JUWELS supercomputer at the Jülich Supercomputing Centre (JSC) for deeper insight into radio relics, cloud of diffuse radio wave emission, often found in galaxy clusters.

Astrophysics

Principal Investigator: Dr. Jenny Sorce , Ludwigs-Maximillians-Universität München

HPC Platform used: SuperMUC-NG at LRZ

Local Project ID: Pn29me

Researchers at the Ludwigs-Maximillians Universität München are focused on studying the concentration of galaxies called the Coma, which is made up of more than 1,000 galaxies in our interstellar neighborhood. Using experimental facilities and world-class computng, the team was able to simulate the Coma cluster in unprecedented detail.

Astrophysics

Principal Investigator: Dr. Klaus Dolag , LMU München

HPC Platform used: SuperMUC-NG at LRZ

Local Project ID: pn68na

Our Cosmic Home, which is the local volume of the Universe centered on us, contains very prominently visible structures, extending over almost one billion light-years. Such structures, ranging from the Local Group over the Local Void and the most prominent galaxy clusters like Virgo, Perseus, Coma and many more, represent a formidable site where extremely detailed observations exist. Therefore, cosmological simulations of the formation of galaxies and galaxy clusters within the Local Universe, rather than any other, randomly selected part of the cosmic web, are perfect tools to test our formation and evolution theories of galaxies and galaxy clusters down to the details. However, at these detailed levels, such simulations are facing various…

Astrophysics

Principal Investigator: Prof. Luciano Rezzolla , Goethe Universität Frankfurt

HPC Platform used: SuperMUC-NG at LRZ

Local Project ID: pn56bi

Gravitational wave detectors such as LIGO, VIRGO, and KAGRA, have brought about an era of multi-messenger astronomy that has given new insights into the merger of binary compact objects. In all cases, the ability to constrain the characteristics of the compact objects is very limited, especially in the absence of an electro-magnetic (EM) counterpart. Gravitational wave events such as GW190425, however, present a very unique opportunity to study the mass gap regime where the binary could consist of either a black hole and a neutron star (BHNS) or a highly asymmetric neutron star binary (BNS).

Astrophysics

Principal Investigator: Prof. Dr. Volker Springel , Max Planck Institute for Astrophysics, Garching

HPC Platform used: SuperMUC-NG at LRZ

Local Project ID: pn34mo

The amazing progress in observational cosmology over the last decades has brought many surprises. Perhaps the most stunning is that we live in a Universe where most of the matter (~85%) is comprised of yet unidentified collisionless dark matter particles, while ordinary baryons produced in the Big Bang make up only a subdominant part (~15%). The real physical nature of dark energy, as well as the mass of the neutrinos which contribute a tiny admixture of “hot” dark matter, are profound and fundamental open questions in physics. To make further progress, this firmly established standard cosmological model will be subjected to precision tests in the coming years that are far more sensitive than anything done thus far.

Astrophysics

Principal Investigator: Minna Palmroth , Finnish Centre of Excellence in Research of Sustainable Space, Helsinki (Finland)

HPC Platform used: Hawk at HLRS

Local Project ID: SIMPLE

Space is the finest plasma laboratory one can reach, hence many of the fundamental and universal physics discoveries of to the fourth state of matter – plasma – root to space physics. The near-Earth space is the only place one can send spacecraft to study the variability of plasma ranging from meters to millions of kilometres and from milliseconds to hundreds of years. However, one can send only a few satellites on a few orbits, making near-Earth space environment modelling crucial. To model the near-Earth space accurately, one requires a good resolution for the 3D position space, and additional 3D space for particle distributions— demanding computing performance that easily can reach the limits of any available supercomputer. 

Astrophysics

Principal Investigator: Luciano Rezzolla , Institute for Theoretical Physics, Goethe University Frankfurt (Germany)

HPC Platform used: SuperMUC-NG of LRZ

Local Project ID: pn56bi

This ongoing project aims at investigating the long-term evolution of a merging binary system of two neutron stars. The investigation conducted within this project is well aligned with the past research conducted by the Relastro group in Frankfurt and is motivated by the gravitational-wave detection GW170817 and its electromagnetic counterpart, the
so-called kilonova. This kilonova signal is produced by the nuclear processes within the dense and neutron rich mass that is ejected during the merger. Since a lot of mass is ejected during the longterm postmerger evolution, it is crucial to investigate this part via state-of-the-art simulations in order to fully understand the observation.

Astrophysics

Principal Investigator: Daniel Seifried , I. Physikalisches Institut, Universität zu Köln

HPC Platform used: SuperMUC and SuperMUC-NG of LRZ

Local Project ID: pr94du

Researchers investigated the formation and evolution of molecular clouds, i.e. the nurseries of star formation, by means of 3D magneto-hydrodynamical simulations. These molecular clouds, which are embedded in a galactic disk like our Milky Way, were modelled with a high spatial resolution using a smart zoom-in approach relying on the adaptive mesh refinement technique. As the modelled molecular clouds were embedded in a realistic astrophysical environment, it was possible to study their detailed evolution, e.g. the impact of supernova explosions and radiation from nearby massive stars. Moreover, the research team modelled the chemical evolution of these clouds as well as their dynamics and complex internal structure.

Astrophysics

Principal Investigator: Anna Therese Phoebe Schauer , Zentrum für Astronomie, Universität Heidelberg

HPC Platform used: SuperMUC and SuperMUC-NG of LRZ

Local Project ID: pr53ka

Before the first stars formed more than 13 billion years ago, the gas of the Universe consisted of hydrogen, helium, and lithium only. Elements necessary for life, eg carbon or oxygen, are produced by stars, and it is of fundamental importance to understand how the first stars formed. With a large allocation on SuperMUC and SuperMUC-NG, state-of-the-art numerical simulations were performed to mimic these first star formation regions. In these high-resolution simulations, two effects – a so-called Lyman-Werner background and streaming velocities – that delay star formation globally were included. It could be demonstrated for the first time that the combination of both effects results in an even more delayed formation of the first stars.

Astrophysics

Principal Investigator: Maarit Käpylä , Aalto University, Department of Computer Science, Astroinformatics Group, Finland, and Max Planck Institute for Solar System Research, SOLSTAR group, Germany

HPC Platform used: SuperMUC-NG of LRZ

Local Project ID: pn98qu

To understand solar and stellar magnetic field evolution combining local and global numerical modelling with long-term observations is a challenging task: even with state-of-the-art computational methods and resources, the stellar parameter regime remains unattainable. Our goal is to relax some approximations, in order to simulate more realistic systems, and try to connect the results with theoretical predictions and state-of-the-art observations. Higher resolution runs undertaken in this project will bring us into an even more turbulent regime, in which we will be able to study, for the first time, the interaction of small- and large-scale dynamos in a quantitative way.

Astrophysics

Principal Investigator: Sebastiano Bernuzzi, Bernd Brügmann , Friedrich-Schiller-Universität Jena

HPC Platform used: SuperMUC and SuperMUC-NG of LRZ

Local Project ID: pn56zo

The project developed multiscale 3+1D simulations of binary neutron mergers in numerical general relativity for applications to multi-messenger astrophysics. It focused on two aspects: (i) the production of high-quality gravitational waveforms suitable for template design and data analysis, and (ii) the investigation of merger remnants and ejecta with sophisticated microphysics, magnetic-fields induced turbulent viscosity and neutrino transport schemes for the interpretation of kilonova signals. The simulations led to several breakthroughs in the first-principles modeling of gravitational-wave and electromagnetic signal, with direct application to LIGO-Virgo's GW170817 and counterparts observations. All data products are publicly released.

Astrophysics

Principal Investigator: Hubert Klahr , Max-Planck-Institut für Astronomie, Heidelberg (Germany)

HPC Platform used: JUQUEEN and JUWELS of JSC

Local Project ID: chhd19

MPIA scientists have developed a planetesimal formation model based on high-resolution hydro-dynamical simulations performed on JSC HPC systems. The simulations were used to model disk turbulence and its two effects on the dust, the mixing and diffusion of the dust on large scales but also the concentration of dust on small scales. This research helped to better understand the efficiency of these processes and to derive initial mass functions for planetesimals and gas giant planets to predict when and where planetesimals and Jupiter-like planets should form and of which size they will be. This is a fundamental step forward in understanding the formation of our own solar system as well as of the many planetary systems around other stars.

Astrophysics

Principal Investigator: Petri Käpylä , Institut für Astrophysik, Georg-August-Universität Göttingen

HPC Platform used: SuperMUC and SuperMUC-NG of LRZ

Local Project ID: pr27li

The outer layers of the Sun are convectively unstable such that heat and momentum are transported by material motions. These motions are thought to be responsible for the large-scale magnetism and differential rotation of the Sun. Employing a more realistic description of the heat conductivity in our simulations than in previous studies, we demonstrate that stellar convection is highly non-local. Furthermore, we found substantial formally stably stratified but fully mixed layers that can cover up to 40 per cent of the solar convection zone. These results are reshaping our picture of stellar convection.

Astrophysics

Principal Investigator: Jörg Büchner , Max-Planck-Institute for Solar System Research, Göttingen

HPC Platform used: SuperMUC and SuperMUC-NG

Local Project ID: pr74vi

The hot and dilute astrophysical plasmas - from Solar to galactic scales - are inherently turbulent. The turbulence determines transport and structure formation in accretion disks, in the interstellar medium, in clusters of galaxies as well as their observable radiation. Due to its routing in microscopic kinetic processes the turbulence of astrophysical plasmas is, however, not well understood, yet. Utilizing state-of-the-art microphysics particle-in-cell codes in this project self-consistent 3D electromagnetic kinetic simulations were performed to simulate the kinetic turbulence inherently linked with two fundamental processes of energy conversion in the Universe – collisionless shock waves and magnetic reconnection.

Astrophysics

Principal Investigator: Hans-Thomas Janka , Max Planck Institute for Astrophysics, Garching

HPC Platform used: SuperMUC and SuperMUC-NG of LRZ

Local Project ID: pn69ho, pr53yi

Core-collapse supernovae are among the most energetic events in the Universe and can be as bright as a galaxy. They mark the violent, explosive death of massive stars, whose iron cores collapse to the most exotic compact objects known as neutron stars and black holes. In this project self-consistent 3D simulations with state-of-the-art microphysics were performed for the explosion of a ~19 solar-mass star, whose final 7 minutes of convective oxygen-shell burning had been computed, too. It could be demonstrated that explosions by the neutrino-driven mechanism can produce powerful supernovae with energies, radioactive nickel ejecta, and neutron-star masses and kick velocities in agreement with observations, in particular Supernova 1987A.

Astrophysics

Principal Investigator: Hans-Thomas Janka , Max Planck Institute for Astrophysics, Garching

HPC Platform used: SuperMUC-NG of LRZ

Local Project ID: pr53yi

Core-collapse supernovae are among the most energetic events in the Universe and can be as bright as a galaxy. They mark the violent, explosive death of massive stars, whose iron cores collapse to the most exotic compact objects known as neutron stars and black holes. In this project self-consistent 3D simulations with state-of-the-art microphysics were performed for the explosion of a ~19 solar-mass star. It could be demonstrated that muon formation in the hot neutron star, which had been ignored in supernova models so far, leads to a faster onset of the explosion. The effects of muons thus over-compensate the delay of the explosion caused by low resolution, where numerical viscosity impedes the growth of hydrodynamic instabilities.

Astrophysics

Principal Investigator: Daniel Ceverino , Zentrum für Astronomie, Institut für Theoretische Astrophysik, Universität Heidelberg

HPC Platform used: SuperMUC and SuperMUC-NG of LRZ

Local Project ID: pr92za

The FirstLight project at LRZ is a large database of numerical models of galaxy formation that mimic a galaxy survey of the high-redshift Universe, before and after the Reionization Epoch. This is the largest sample of zoom simulations of galaxy formation with a spatial resolution better than 10 pc. This database improves our understanding of cosmic dawn. It sheds light on the distribution of gas, stars, metals and dust in the first galaxies. This mock survey makes predictions about the galaxy population that will be first observed with future facilities, such as the James Webb Space Telescope and the next generation of large telescopes.

Astrophysics

Principal Investigator: Luciano Rezzolla , Institute for Theoretical Physics, Goethe University FIAS – Frankfurt Institute for Advanced Studies

HPC Platform used: SuperMUC and SuperMUC-NG

Local Project ID: pr27ju

Two major events are responsible for what is considered the “golden age” of relativistic astrophysics. One is the detection of gravitational waves from merging neutron stars heralding the beginning of the multimessenger age. The other is the effort of the Event Horizon Telescope collaboration culminating in the first image of a black hole. Both events have been aided by simulations that require HPC. With this project, several studies could be conducted well alligned with these type of simulations expanding our knowledge about these important astrophysical events.

Astrophysics

Principal Investigator: Tim Dietrich , University of Potsdam, Dutch National Institut for Subatomic Physics Amsterdam

HPC Platform used: SuperMUC and SuperMUC-NG of LRZ

Local Project ID: pr48pu

Neutron stars are ultracompact stars in which densities above the nuclear saturation densities are reached and that provide one of the best laboratories to test nuclear physics principles. Within this project, researchers perform 3+1-dimensional numerical-relativity simulations studying the last few orbits before the merger of two of these stars. In fact, a binary neutron star merger is one of the most energetic phenomena in our Universe and is accompanied by a variety of electromagnetic signatures and with characteristic gravitational-wave signatures. With the help of these simulations existing theoretical models can be developed and verified and the growing field of multi-messenger astronomy is supported. 

Astrophysics

Principal Investigator: Jose Oñorbe , Universidad de Sevilla, Spain

HPC Platform used: JUWELS of JSC

Local Project ID: CPRA100

The intergalactic medium, the low density gas that lies between galaxies, contains vast information about how the universe evolved and when the first stars formed. In order to provide a solid theoretical platform for current and upcoming observations of the properties of this gas, a series of hydrodynamical cosmological simulations using the Nyx code were run on JUWELS on JSC. These simulations had an unprecedented dynamical range and used a novel approach to account for the challenging inhomogeneous radiation that must be included in these type of calculations.

Astrophysics

Principal Investigator: Wolfgang Hillebrandt , Max-Planck-Institut für Astrophysik, Garching b. München

HPC Platform used: JUWELS of JSC

Local Project ID: hmu14

Supernovae of Type Ia are modeled as thermonuclear explosions of a carbon-oxygen white dwarf stars. The way these trigger the explosive burning, however, is still unclear. This project performs hydrodynamic simulations that give insights into possible explosion mechanisms. With its pipeline extending from explosion simulation to the derivation of synthetic observables, the project allows for a direct comparison with astronomical observations thus scrutinizing the modeled scenarios.

Astrophysics

Principal Investigator: Friedrich Röpke , Zentrum für Astronomie der Universität Heidelberg, Institut für Theoretische Astrophysik und Heidelberger Institut für Theoretische Studien

HPC Platform used: JUQUEEN and JUWELS of JSC

Local Project ID: chwb07

Classical stellar models are formulated in one spatial dimension and parameterize dynamical multidimensional effects. While successful in a qualitative description of how stars evolve, such models lack predictive power. Multidimensional hydrodynamic simulations of critical phases and processes are still extremely challenging but have become feasible due to improved numerical techniques and increasing computational power. This project performs such simulations aiming at an improved understanding of the physics ruling stellar structure and evolution. As an example, a simulation of convective helium-shell burning in a massive star is discussed.

Astrophysics

Principal Investigator: Jenny G. Sorce(1), Klaus Dolag(2) , (1) Leibniz-Institut für Astrophysik Potsdam/AIP (Germany) and Centre de Recherche Astrophysique de Lyon (France), (2) Universitäts-Sternwarte, Ludwig-Maximilians-Universität München (Germany)

HPC Platform used: SuperMUC of LRZ

Local Project ID: pr74do

The neighbourhood in the immediate vicinity of the Milky Way is known as the “Local Group”. It is a binary system composed of two averaged sized galaxies (the Milky Way and Andromeda) dominating a volume that is roughly ~7 Mpc in diameter. At a distance of around 15Mpc, the Virgo cluster comes into view as the main defining feature of our neighbourhood on these scales. Beyond Virgo, a number of well known and well observed clusters like Centaurus, Fornax, Hydra, Norma, Perseus and Coma dominate the night sky. This is our cosmic neighbourhood. The goal of this project is, for the first time, to perform targeted, state of the art hydro-dynamical simulations covering this special region of the universe and to compare the results with various…

Astrophysics

Principal Investigator: Felix Spanier(1), Anne Stockem-Novo(2) , (1) Karlsruhe Institut für Technologie, Eggenstein-Leopoldshafen (Germany), (2) Ruhr-Universität Bochum (Germany)

HPC Platform used: SuperMUC of LRZ

Local Project ID: pr74se

Active galactic nuclei (AGN) are powerful emitters of photons in energy ranges from few millielectron volts (meV) to several teraelectron volts (TeV). These sources show variabilities as fast as a few minutes. It is believed that the emission originates from particles accelerated in shock waves in the jet of AGN. Observational data, however, is too sparse to constrain radiation models. Therefore, light curves (i.e. temporal data) are used to constrain models further. Using the Particle-in-Cell method to investigate shock collisions, this project aims at gaining more detailed insight into a special case of variability.

Astrophysics

Principal Investigator: Klaus Dolag , Universitäts-Sternwarte, Ludwig-Maximilians-Universität München (Germany)

HPC Platform used: SuperMUC of LRZ

Local Project ID: pr83li, pr86re

The outcome of a large set of cosmological, hydro-dynamical simulations from the project Magneticum now became made available to the general community through operating a cosmological simulation web portal. Users are able to access data products extracted from the simulations via a user-friendly web interface, browsing through visualizations of cosmological structures while guided by meta data queries helping to select galaxy clusters and galaxy groups of interest. Several services are available for the users: (I) ClusterInspect; (II) SimCut (raw data access); (III) Smac (2D maps); (IV) Phox (virtual X-ray observations, taking the specifications of various, existing and future X-ray telescopes into account.

Astrophysics

Principal Investigator: Jenny Sorce , Leibniz-Institut für Astrophysik Potsdam (Germany) and Centre de Recherche Astrophysique de Lyon (France)

HPC Platform used: SuperMUC of LRZ

Local Project ID: pr74je

Galaxy clusters are large reservoirs of galaxies. As such they are perfect objects of studies to unravel the mysteries of galaxy formation and evolution in dense environments. At ~50 million light-years away from Earth, the Virgo cluster, a gathering of more than a thousand galaxies, is our closest cluster-neighbour. Its proximity permits deep observations. Cosmological numerical simulations of the cluster constitute the numerical counterparts to be compared with observations to test our theoretical models. In such simulations, dark matter (nature of most of the matter in the Universe) and baryons (visible matter) follow physical laws to reproduce our closest cluster-neighbour and its galaxies in a simulated box across cosmic time.

Astrophysics

Principal Investigator: Hans-Thomas Janka , Max-Planck-Institut für Astrophysik, Garching (Germany)

HPC Platform used: SuperMUC of LRZ

Local Project ID: pr53yi

Traditionally, numerical simulations of core-collapse supernovae have been performed with spherically symmetric initial models for the progenitor stars, because stellar evolution is computed with this restriction. Recently, however, it has been demonstrated that pre-collapse asymmetries in the convectively burning oxygen shell can have an impact on the explosion by enhancing turbulence behind the supernova shock. In this project researchers simulated the final seven minutes of oxygen burning and the subsequent collapse of a 19 solar-mass star in order to investigate the consequences of pre-collapse asymmetries for the supernova explosion.

Astrophysics

Principal Investigator: Christoph Federrath (1), Ralf S. Klessen (2) , (1) Research School of Astronomy and Astrophysics, Australian National University (ANU), (2) Zentrum für Astronomie, Institut für Theoretische Astrophysik und Interdisziplinäres Zentrum für Wissenschaftliches Rechnen, Universität Heidelberg (Germany)

HPC Platform used: SuperMUC of LRZ

Local Project ID: pr32lo

Understanding turbulent gases and fluids is critical for a wide range of terrestrial and astrophysical applications. Here we present the world's largest turbulence simulation to date. This GCS Large-Scale Project on SuperMUC consumed 45 million core hours and produced 2 PB of data. It is the first and only simulation to bridge the scales from supersonic (Mach > 1) to subsonic (Mach < 1) flow and resolves the sonic scale (where the Mach number = 1). The sonic scale is a key ingredient for star formation models and may determine the size of filamentary structures in the interstellar medium.

Astrophysics

Principal Investigator: Joakim Rosdahl , Centre de Recherche Astrophysique de Lyon (France)

HPC Platform used: SuperMUC of LRZ

Local Project ID: pr53na

The formation of the first galaxies marked the end of the cosmological dark ages. Radiation from the first stars ionized and heated inter-galactic gas. As these ionized gas bubbles grew and percolated, the whole Universe was transformed from a dark, cold, neutral state into a hot ionized one, about a billion years after the Big Bang. The SPHINX cosmological radiation-hydrodynamics simulations of the first billion years are designed to understand the formation of the first galaxies and how they contributed to reionization via the interplay of star formation, stellar radiation, and powerful supernova explosions that disrupt galaxies and allow their radiation to escape into inter-galactic space.

Astrophysics

Principal Investigator: Hans-Thomas Janka , Max-Planck-Institut für Astrophysik, Garching (Germany)

HPC Platform used: SuperMUC of LRZ

Local Project ID: pr62za

The "ray-by-ray" approximation is a widely used simplification of the time-dependent, six-dimensional transport of all neutrino species in core-collapse supernovae. It reduces the dimensionality of the computationally challenging problem by assuming that non-radial flux components are negligible. This leads to the solution of three-dimensional (radius-, energy-, and angle-dependent) transport equations for all angular directions of the spatial polar grid. Such a task can be extremely efficiently parallelized also on huge numbers of computing cores. In this project 3D simulations were performed to test this approximation and could demonstrate its validity.

Astrophysics

Principal Investigator: Stefanie Walch-Gassner , I. Physikalisches Institut, Universität zu Köln

HPC Platform used: SuperMUC of LRZ

Local Project ID: pr62su

Molecular clouds form out of the diffuse interstellar medium (ISM) within galactic disks and continuously accrete gas and interact with their surroundings as they evolve. Hence the evolution of turbulent, filamentary molecular clouds has to be modeled at the same time as the surrounding multiphase ISM. In the SILCC-ZOOM project, we simulate molecular cloud formation, the star formation within them, and their subsequent dispersal by stellar feedback on sub-parsec scales in 3D, AMR, MHD simulations with the FLASH code including self-gravity, radiative transfer, and a chemical network.

Astrophysics

Principal Investigator: Luciano Rezzolla , Institute for Theoretical Physics, Goethe University Frankfurt (Germany)

HPC Platform used: SuperMUC of LRZ

Local Project ID: pr62do

The age of multi-messenger gravitational wave astronomy has arrived. The simultaneous detection of gravitational and electromagnetic waves from merging neutron stars has illustrated the importance of having high resolution numerical relativity simulations, performed on SuperMUC, available to disentangle the complex interplay of nuclear physics, neutrino physics, and strong field gravity. Using these simulations, it is possible to study matter at densities unreachable with terrestrial experiments and determine the origin of the heavy elements in the universe.

Astrophysics

Principal Investigator: Tim Dietrich , Max Planck Institute for Gravitational Physics, Potsdam (Germany)

HPC Platform used: SuperMUC of LRZ

Local Project ID: pr48pu

The collision of two neutron stars is one of the most violent events in the Universe. The extreme conditions, with densities of about one hundred million tons per cubic centimeter and gravity hundred billion times that of Earth gravity, cannot be tested on Earth, which makes these events a perfect laboratory to study matter at extreme limits. Using advanced numerical relativity simulations, scientists study the phenomena close to the merger of the two neutron stars to extract information about the emitted gravitational wave and electromagnetic signals.

Astrophysics

Principal Investigator: David Hilditch and Bernd Brügmann , Theoretisch-Physikalisches Institut, Friedrich-Schiller-Universität Jena (Germany)

HPC Platform used: SuperMUC of LRZ

Local Project ID: pr87nu

General relativity describes the gravitational interaction as the curvature of spacetime. This involves complicated partial differential equations, and consequently extreme scenarios can be treated only by numerical simulations. In this project spacetimes close to the critical threshold of black hole formation were evolved on SuperMuc. These computations were performed using bamps, a new massively parallel code for numerical relativity. The spacetimes constructed constitute the most extreme regime imaginable - that in which cosmic censorship itself may be violated and the black hole singularity could be seen by distant observers. 

Astrophysics

Principal Investigator: Hubert Klahr , Max-Planck-Institut für Astronomie, Heidelberg (Germany)

HPC Platform used: JUQUEEN of JSC

Local Project ID: hhd19

Planetesimals are kilometre-sized planetary building blocks in the early solar system. Scientists pioneered a scenario in which turbulent concentrations of the icy and dusty material leads to sufficiently large densities in which self-gravity dominates over gas shear and tidal forces of the star. As a consequence, the material collapses spontaneously under its own weight into planetesimals. Therefore, the motion of many millions of particles in magneto-hydro-dynamically and particle driven turbulence and include the gravity among gas and particles are all simulated in one huge simulation. The goal is to link the observations of dust around young stars in a quantified way to an initial mass distribution of planetesimals.

Astrophysics

Principal Investigator: Alfred Müller , Institut für Atom- und Molekülphysik, Universität Giessen (Germany)

HPC Platform used: Hazel Hen of HLRS

Local Project ID: PAMOP/PAMOP2

An international group of scientists leverages high-performance computing to support current and future measurements of atomic photoionization cross-sections at various synchrotron radiation facilities, ion-atom collision experiments, together with plasma, fusion and astrophysical applications. In their work they solve the Schrödinger or Dirac equation using the R-matrix or R-matrix with pseudo-states approach from first principles. Cross-sections and rates for radiative charge transfer, radiative association, and photodissociation collision processes between atoms and ions of interest for several astrophysical applications are presented.

Astrophysics

Principal Investigator: Dylan Nelson(1) and Annalisa Pillepich(2) , (1) MPA Garching (Germany), (2) MPIA Heidelberg (Germany)

HPC Platform used: Hazel Hen of HLRS

Local Project ID: GCS-dwar

Modern simulations of galaxy formation, which simultaneously follow the co-evolution of dark matter, cosmic gas, stars, and supermassive black holes, enable us to directly calculate the observable signatures that arise from the complex process of cosmic structure formation. TNG50 is an unprecedented ‘next generation’ cosmological, magneto-hydrodynamical simulation -- the third and final volume of the IllustrisTNG project. It captures spatial scales as small as ~100 parsecs, resolving the interior structure of galaxies, and incorporates a comprehensive model for galaxy formation physics.

Astrophysics

Principal Investigator: Volker Springel , Heidelberg Institute for Theoretical Studies, Heidelberg University, and Max-Planck Institute for Astrophysics (Germany)

HPC Platform used: Hazel Hen of HLRS

Local Project ID: GCS-ILLU

Hydrodynamical simulations of galaxy formation have now reached sufficient physical fidelity to allow detailed predictions for their formation and evolution over cosmic time. The aim of this project is to carry out a new generation of structure formation simulations, IllustrisTNG, that reach sufficient volume to make accurate predictions for clustering on cosmologically relevant scales, while at the same time being able to compute detailed galaxy morphologies, the enrichment of diffuse gas with metals, and the amplification of magnetic fields during structure growth.

Astrophysics

Principal Investigator: Ana-Catalina Plesa , Institute of Planetary Research, Planetary Physics, German Aerospace Center/DLR, Berlin (Germany)

HPC Platform used: Hornet of HLRS

Local Project ID: MATHECO

The large amount of data returned by several space missions to the terrestrial planets has greatly improved our understanding of the similarities and differences between the innermost planets of our Solar System. Nevertheless, their interior remains poorly known since most of the data is related to surface processes. In the absence of direct data of the interior evolution of terrestrial planets, numerical simulations of mantle convection are an important mean to reconstruct the thermal and chemical history of the interior of the Earth, Moon, Mercury, Venus and Mars. In this project, run on Hornet of HLRS, researchers used the mantle convection code Gaia to model the thermal evolution of terrestrial planets and in particular the early stage…

Astrophysics

Principal Investigator: Ralf Klessen(1), Christoph Federrath (2) , (1) Universität Heidelberg, Germany (2) Australian National University

HPC Platform used: SuperMUC of LRZ

Local Project ID: pr48pi

Interstellar turbulence shapes the structure of the multi-phase interstellar medium (ISM) and is a key process in the formation of molecular clouds as well as the build-up of star clusters in their interior. The key ingredient for our theoretical understanding of ISM dynamics and stellar birth is the sonic scale in the turbulent cascade, which marks the transition from supersonic to subsonic turbulence and produces a break in the turbulence power spectrum. To measure this scale and study the sonic transition region in detail, scientists, for the first time, ran a simulation with the unprecedented resolution of 10,0483 grid cells.

Astrophysics

Principal Investigator: Hans-Thomas Janka , Max Planck Institute for Astrophysics, Garching (Germany)

HPC Platform used: SuperMUC of LRZ

Local Project ID: pr48de

Recently the first three-dimensional simulations have confirmed the long-standing hypothesis that the neutrino-driven mechanism, supported by violent hydrodynamic instabilities and turbulent mass flows, can explain supernova explosions of stars with more than 8−10 solar masses. Further consolidation of this mechanism and a deeper theoretical understanding of its functioning require the exploration of a broader variety of progenitor stars and of dependences on the initial conditions prior to iron-core collapse. In this Gauss project the influence of stellar rotation, perturbations in the convective oxygen-burning layer, and of large mass-infall rates due to high core compactness in very massive progenitor stars were explored.

Astrophysics

Principal Investigator: Wolfram Schmidt , Institut für Astrophysik, Universität Göttingen, and Hamburger Sternwarte, Universität Hamburg (Germany)

HPC Platform used: SuperMUC of LRZ

Local Project ID: pr84vo

The modelling of star formation and feedback processes such as supernova explosions is a longstanding problem in numerical simulations of cosmological structure formation because the internal structure of galaxies cannot be resolved in sufficient detail even on very powerful supercomputers. For this reason, star formation and stellar feedback are treated as so-called subgrid physics. The aim of our project is to combine standard recipes for star formation in simulations on cosmological scales with a subgrid-scale model for numerically unresolved turbulence, which allows us to study the influence of turbulence on star formation and the mixing of metals expelled by supernova explosions in galaxies. It is believed that, in addition to…

Astrophysics

Principal Investigator: Hans-Thomas Janka , Max Planck Institute for Astrophysics, Garching (Germany)

HPC Platform used: SuperMUC of LRZ

Local Project ID: pr48ra

First self-consistent, first-principle simulations in three dimensions have provided support for the viability of the neutrino-driven mechanism as an explanation of supernova explosions of stars with more than 8−10 solar masses. While these results respresent fundamentally important progress in our understanding of how massive stars terminate their lives, the enormous complexity and computational demand of the involved neutrino physics set severe resolution limitations to current full-scale supernova models. In this project, the numerical convergence of the present simulations were investigated.

Astrophysics

Principal Investigator: Filippo Galeazzi, Luciano Rezzolla , Institute for Theoretical Physics, Goethe University Frankfurt (Germany)

HPC Platform used: SuperMUC of LRZ

Local Project ID: pr48fa

Leveraging the HPC infrastructure of LRZ, researchers at the Goethe University in Frankfurt/Main employ a series of in-house developed cutting-edge numerical methods to simulate in full general relativity the inspiral, merger, and collapse of neutron stars. The computationally intense, fully parallel simulations incorporate relativistic hydrodynamics, nuclear finite-temperature equations of state, and an approximate treatment of neutrino emission and absorption. The results, obtained by measuring gravitational waves, can provide important information on the properties of matter at nuclear densities.

Astrophysics

Principal Investigator: Stefan Gottlöber , Leibniz-Institut für Astrophysik Potsdam (Germany)

HPC Platform used: SuperMUC of LRZ

Local Project ID: h009z

An international team of scientists performed a series of Constrained Simulations to study Near Field Cosmology. These high-resolution simulations allowed the astrophysicists, for the first time, to study the formation of the Local Group in the right cosmic environment. 

Astrophysics

Principal Investigator: (1)Bernd Brügmann, (2)Tim Dietrich , (1)Friedrich-Schiller-University, Jena, (2)Max Planck Institute for Gravitational Physics (Albert Einstein Institute), Potsdam-Golm (Germany)

HPC Platform used: SuperMUC of LRZ

Local Project ID: pr48pu

The recent observations of gravitational waves (GWs) marked a breakthrough and inaugurated the field of GW astronomy. To extract information from a detection, the measured signal needs to be cross-correlated with a template family. However, due to the nonlinearity of Einstein’s equations, numerical simulations have to be used to study systems with gravitational fields strong enough to emit GWs. This project focused on the simulation of systems consisting of two neutron stars and investigated the effect of the mass ratio and the influence of the spin of the individual stars.

Astrophysics

Principal Investigator: Ulrich Hansen , Institut für Geophysik, Westfälische Wilhelms-Universität Münster (Germany)

HPC Platform used: JUQUEEN of JSC

Local Project ID: hms18

Researchers at the University of Münster investigated precession driven flows in planets by direct numerical simulations on the JUQUEEN cluster. Precession of the rotation axis is an often neglected driving mechanism for flows in planetary cores, a field of research were other scientists mainly focus on the influence of thermal or chemical effects. As an additional complication that moves the models closer to the physical reality, the project considered the spheroidal shape of the planet, whereas previous research has been focused on the idealized case of a perfect sphere.

Astrophysics

Principal Investigator: Bruno Giacomazzo , University of Trento and INFN-TIFPA, Trento, Italy

HPC Platform used: SuperMUC of LRZ

Local Project ID: GRSimStar

What can we learn from some of the most powerful explosions in the universe? Researchers in Italy, USA, and Japan joined forces to study, via computer simulations in general relativity, what happens when two neutron stars in a binary system finally merge. Besides black holes, neutron stars are the most compact objects ever observed. Their collisions can produce bright electromagnetic emission and strong gravitational waves. Understanding how to relate the different signals with the properties of neutron stars may allow us to understand how matter behaves in conditions so extreme that cannot be reproduced on Earth.

Astrophysics

Principal Investigator: Yannick Bahé , Max Planck Institut for Astrophysics, Garching (Germany)

HPC Platform used: Hazel Hen of HLRS

Local Project ID: gcs-hydra

Why do galaxies that live in the enormous structures known as galaxy clusters look different from normal, isolated galaxies, like our Milky Way? To answer this question, astrophysicists have created the Hydrangea simulations, a suite of 24 high-resolution cosmological hydrodynamical simulations of galaxy clusters. Containing over 20,000 cluster galaxies in unprecedented detail and accuracy, these simulations are giving astrophysicists a powerful tool to understand how galaxies have formed and evolved in one of the most extreme environments of our Universe.

Astrophysics

Principal Investigator: Peter Hauschildt , Hamburger Sternwarte, Universität Hamburg (Germany)

HPC Platform used: Hazel Hen of HLRS

Local Project ID: PP14112588

Understanding the light emitted by (magnetically) active cool stars (‘M dwarfs’) is a major challenge for astrophysics. In this project, scientists use their PHOENIX/3D code to simulate the light emitted by a ‘box’ inside the outer layers of an active M dwarf in detail. The temperatures and pressures inside the box are taken from an existing gas dynamics simulation (including magnetic field effects) by S. Wedemeyer (Oslo). The computational requirements of detailed non-equilibrium 3D radiative transfer simulations are staggering and require the largest supercomputers on Earth.

Astrophysics

Principal Investigator: Peter Hauschildt , Hamburger Sternwarte, Universität Hamburg (Germany)

HPC Platform used: JUQUEEN of JSC

Local Project ID: hhh15

Researchers finished the implementation and verification of a 3D non-local thermodynamic equilibrium (NLTE/3D) module for the PHOENIX/3D model atmosphere simulation code. The methods were extended to also allow NLTE modelling of molecular lines (here: CO) and then used to model the radiation from parameterized star-spots to investigate the effects of detailed 3D radiation transport on observables.

Astrophysics

Principal Investigator: Minna Palmroth , Earth Observation Finnish Meteorological Institute, Helsinki (Finland)

HPC Platform used: Hazel Hen of HLRS

Local Project ID: PP14112573

Space weather is an increasingly important aspect for our technology-dependent society. Modelling space weather is difficult, however, a Finnish team has succeeded in something that was said to be impossible: an accurate simulation of the large-scale near-Earth space environment. PRACE Tier-0 grant from Hazel Hen (HLRS, Stuttgart) both allowed the Vlasiator team to discover new space physics phenomena, and significantly helped in the acceptance of the second European Research Council grant awarded to the project PI in fall 2015.

Astrophysics

Principal Investigator: David Weir , Department of Mathematics and Natural Sciences, University of Stavanger (Norway)

HPC Platform used: Hazel Hen of HLRS

Local Project ID: PP14112721

Gravitational waves are ripples in spacetime, predicted by Einstein already a century ago. With the announcement earlier this year that gravitational waves had been successfully detected from two black holes merging, attention now turns to other potential sources of gravitational waves. Such sources include dramatic events that may have occurred very early in the history of the universe. Understanding these other sources also informs the design of future gravitational wave detectors, such as the European Space Agency (ESA) project eLISA.

Astrophysics

Principal Investigator: Giovanni Lapenta , KU Leuven (Belgium)

HPC Platform used: SuperMUC of LRZ

Local Project ID: pr87di

The process of magnetic reconnection — when magnetic fields in plasma reconfigure and explosively release thermal and kinetic energy — is only just beginning to be understood. Professor Giovanni Lapenta has been carrying out simulations on SuperMUC of how these events can cause chain reactions that very quickly fill vast volumes of space. This data is now being verified with the recent NASA Magnetospheric MultiScale Mission that is measuring magnetic reconnection events around the Earth.

Astrophysics

Principal Investigator: Mats Carlsson , Institute of Theoretical Astrophysics, University of Oslo (Norway)

HPC Platform used: SuperMUC of LRZ

Local Project ID: pr85wo

A simulation project, run on SuperMUC, targets the intrinsic physics of the chromosphere in order to understand its mass and energy budgets and transfer mechanisms. Elucidating these is a principal quest of solar physics, a necessary step towards better space-weather prediction, and of interest to general astrophysics using the Sun as a close-up Rosetta-Stone star and to plasma physics using the Sun and heliosphere as a nearby laboratory. The project aims at a breakthrough in our understanding of the solar chromosphere by developing sophisticated radiation-magnetohydrodynamic simulations.

Astrophysics

Principal Investigator: Felix Spanier , Institut für Theoretische Physik und Astrophysik, Universität Würzburg (Germany)

HPC Platform used: SuperMUC of LRZ

An international research collaboration led by the University of Würzburg delved into the subjects of turbulence and particle acceleration in the solar wind by performing highly complex numerical simulations leveraging the particle-in-cell (PiC) approach, a technique used to solve a certain class of partial differential equations thus capable of studying these phenomena. In order to model the complex system of different waves, particles and electromagnetic fields self-consistently, the use of massive computing power such as provided by high performance computing system SuperMUC is inevitable.

Astrophysics

Principal Investigator: Anne Stockem Novo , Ruhr-University Bochum (Germany)

HPC Platform used: JUQUEEN of JSC

Local Project ID: hbo41

The acceleration of charged particles is still one of the most important problems in astrophysics. Cosmic rays, which mainly consist of protons, show a broad spectrum with energies up to 1021 eV, which can be produced in collisionless shocks. However, many questions are still open regarding the acceleration process and the process of shock formation. To study this complex process with non-linear methods, researchers used JUQUEEN to investigate different aspects of the shock formation process and further applications.

Astrophysics

Principal Investigator: Hubertus Klahr , Max-Planck-Institut für Astronomie, Heidelberg (Germany)

HPC Platform used: JUQUEEN of JSC

Local Project ID: hhd19

Scientists of the Max Planck Institute for Astronomy in Heidelberg are using the HPC infrastructure of the Jülich Supercomputing Centre for extensive magneto-hydro-dynamical and million particle simulations of protoplanetary disks to study their evolution and properties. Findings are helping the researchers to understand the processes leading to the formation of planets, moons and asteroids. Their investigations will help to explain the observed diversity in planetary systems and in our own solar system.

Astrophysics

Principal Investigator: Ana-Catalina Plesa , German Aerospace Center/DLR, Berlin (Germany)

HPC Platform used: Hornet of HLRS

Local Project ID: MATHECO

Scientists of the German Aerospace Center Berlin (DLR) exploited the computing capacity of the petascale system Hornet of HLRS to study the convective dynamics and evolution of planetary interiors. The goal of the large-scale simulation project MATHECO (MAntle THErmo-chemical COnvection Simulations), which scaled to 54,000 compute cores of the supercomputer Hornet, was to gain further insights into the cooling history of planets and its influences on volcanic and tectonic surface processes.

Astrophysics

Principal Investigator: Patrick Kilian , Lehrstuhl für Astronomie, Universität Würzburg (Germany)

HPC Platform used: Hermit of HLRS

Local Project ID: EAINRS

Observations show that Earth is constantly bombarded by highly energetic particles that are called cosmic rays. A possible explanation for the origin of the cosmic rays as well as their energy distribution is particle acceleration at shock fronts. Several different physical processes take place there, but due to the large astrophysical distances it is, unfortunately, impossible to study these in-situ. One way out is large scale computer simulations.

Astrophysics

Principal Investigator: Dr. Andreas Pawlik , Max-Planck-Institut für Astrophysik, Garching (Germany)

HPC Platform used: SuperMUC of LRZ

Local Project ID: pr83le

A multi-million compute hours allocation by the Gauss Centre for Supercomputing on HPC system SuperMUC of the Leibniz Supercomputing Centre (LRZ) was used to carry out Aurora, a new set of radiation-hydrodynamical simulations of galaxy formation during reionization. Numerical simulations have emerged as the most powerful tools for the ab initio theoretical treatment of reionization. 

Astrophysics

Principal Investigator: Hans-Thomas Janka , Max-Planck-Institut für Astrophysik, Garching (Germany)

HPC Platform used: SuperMUC of LRZ

Local Project ID: pr86la

The Stellar Core-Collapse Group at the Max Planck Institute for Astrophysics (MPA) is able to conduct the presently most advanced 3D supernova simulations thanks to a suitably constructed description of the neutrino physics and a highly efficient, extremely well parallelized numerical implementation on petascale system SuperMUC. Because neither experiments nor direct observations can reveal the processes at the center of exploding stars, highly complex numerical simulations are indispensable to develop a deeper and quantitative understanding of this hypothetical “neutrino-driven explosion mechanism”, whose solid theoretical foundation is still missing.

Astrophysics

Principal Investigator: Luís O. Silva , GoLP/Instituto de Plasmas e Fusão Nuclear, Lisboa (Portugal)

HPC Platform used: JUQUEEN of JSC

Local Project ID: PRA_2015

A team of researchers has performed the largest simulations of unmagnetised shocks driven in astrophysical conditions to determine the parameters required to excite shocks in the laboratory, and studied the set of complex and nonlinear phenomena involved in these scenarios, such as magnetic field generation and particle acceleration. The simulations have been performed with the state-of-the-art particle-in-cell code OSIRIS.

Astrophysics

Principal Investigator: Stefanie Walch , Universität zu Köln (Germany)

HPC Platform used: SuperMUC of LRZ

Local Project ID: pr62su

A European team of scientists from Cologne, Garching, Heidelberg, Prague and Zurich used GCS HPC resources to model representative regions of disk galaxies using adaptive, three-dimensional simulations at unprecedented resolution and with the necessary physical complexity to follow the full life-cycle of molecular clouds. They aim to provide a self-consistent answer as to how stellar feedback regulates the star formation efficiency of a galaxy, how molecular clouds are formed and destroyed, and how galactic outflows are driven.

Astrophysics

Principal Investigator: Robi Banerjee , Hamburger Sternwarte (Germany)

HPC Platform used: JUQUEEN of JSC

Local Project ID: hhd14

(Using GCS HPC resources, a group of scientist from a number of international institutes were able to prove that very weak magnetic fields can be efficiently amplified during different stages of cosmic evolution. 

Astrophysics

Principal Investigator: Friedrich Röpke , Universität Würzburg (Germany)

HPC Platform used: JUQUEEN of JSC

Local Project ID: hmu13

(Type Ia supernovae, gigantic astrophysical explosions that completely disrupt one star and shine brighter than an entire galaxy consisting of 100 billion stars, have been successfully used to measure distances in the Universe. But what are the stars that give birth to Type Ia supernovae? The answer to this question remains elusive despite advances in modelling and observing these cosmic events over the past decades. From the perspective of theoretical modelling, only detailed multi-dimensional simulations of the explosion process on the most powerful supercomputers offer a way to tackle this long-standing problem. 

Astrophysics

Principal Investigator: Dr. Sascha Husa , Universitat de les Illes Balears (Spain)

HPC Platform used: SuperMUC of LRZ

Local Project ID: pr86va

A team of about 20 scientists working in Europe, India, South Africa and the USA have been involved in an Astrophysics simulation project calculated on LRZ system SuperMUC. The obtained results will allow the efficient detection and identification of gravitational wave events, e.g. to tell apart black holes from neutron stars.

Astrophysics

Principal Investigator: Volker Springel , HITS, Universität Heidelberg (Germany)

HPC Platform used: SuperMUC of LRZ

Local Project ID: pr85je

An international team of scientists at the Heidelberg Institute for Theoretical Studies (HITS), MIT, Harvard University and the University of Cambridge has carried out the “Illustris Simulation” on the SuperMUC and CURIE supercomputers, and created the largest and most sophisticated computational model of cosmic structure formation thus far.

Astrophysics

Principal Investigator: Klaus Dolag , University Observatory Munich (Germany)

HPC Platform used: SuperMUC of LRZ

Local Project ID: pr83li

In project Magneticum, scientists perform simulations of which the most computational intensive one covers a cosmic volume of 1 Gpc3. This allows the researchers, for the very first time, to self consistently study galaxy clusters and groups, galaxies, and active galaxy nuclei (AGNs) within an enormously large volume of the Universe.

Astrophysics

Principal Investigator: Luciano Rezzolla , Max Planck Institute for Gravitational Physics (Albert Einstein Institute), Potsdam-Golm (Germany)

HPC Platform used: SuperMUC of LRZ

Local Project ID: pr32pi

Scientists from the Albert Einstein Institute in Potsdam/Germany used HPC system SuperMUC to study the dynamics of compact-object binaries, i.e. neutron stars and black holes, and to improve our understanding of strong gravity. 

Astrophysics

Principal Investigator: Minna Palmroth , Finnish Meteorological Institute, Helsinki

HPC Platform used: Hermit of HLRS

Local Project ID: PP12061111

The HPC resources of HLRS Stuttgart enabled the world’s first global runs of the near-Earth space using a hybrid-Vlasov approach at highest resolutions. 

Astrophysics

Principal Investigator: Sven Bingert , Max-Planck-Institut für Sonnensystemforschung, Göttingen (Germany)

HPC Platform used: Hermit of HLRS

Local Project ID: BRUSH

Scientists at the Max Planck Institute for Solar System Research in Göttingen employed a three-dimensional numerical model on GCS supercomputer Hermit of HLRS Stuttgart to investigate the heating process of the highly structured and dynamic corona. 

Astrophysics

Principal Investigator: Paolo Padoan , Instituto de Ciencias del Cosmos (ICC), Universidad de Barcelona (Spain)

HPC Platform used: SuperMUC of LRZ

Local Project ID: pr86li

An international team of scientists used GCS supercomputing resources to study the evolution and fragmentation of clouds into stars. The degree of complexity, resulting from the mutual interaction of magnetic fields, gravity, and supersonic turbulence, is such that no complete theory of star formation is available to date. The best way to tackle this problem is to use powerful supercomputers such as petascale system SuperMUC of LRZ.

Astrophysics

Principal Investigator: Gianluigi Bodo , INAF Astrophysical Observatory of Torino (Italy)

HPC Platform used: JUQUEEN of JSC

Local Project ID: PRA057

Accretion discs can power some of the most energetic phenomena in the universe and understanding how they work is very important for the comprehension of different astrophysical problems like how stars are formed or what happens in the central cores of galaxies.

Astrophysics

Principal Investigator: Dr. Ilian Iliev , Astronomy Centre, Physics and Astronomy, University of Sussex (U.K.)

HPC Platform used: SuperMUC of LRZ

Local Project ID: pr86be

Scientists used GCS supercomputing resources to perform a series of state-of-the-art constrained simulations of the first galaxies and their radiative effects in our Local Universe, among the largest and most detailed of their kind, following tens of billions of particles, while also modelling the complex physics involved in this process.

Astrophysics

Principal Investigator: Steffen Heß , Leibniz-Institut für Astrophysik Potsdam/AIP (Germany)

HPC Platform used: JUQUEEN of JSC

Local Project ID: hpo08

Scientists improved and combined methods to simulate the formation of the actual distribution of galaxies and galaxy clusters which allowed them to simulate the density distribution in the local universe up to distances of 670 million light-years.

Astrophysics

Principal Investigator: Dr. Christoph Federrath , Monash University, Clayton (Australia)

HPC Platform used: SuperMUC of LRZ

Local Project ID: h1343

To advance the so far limited knowledge of density probability distribution function and the power spectrum of compressible, supersonic turbulence, a team of astrophysicists compared hydrodynamic models with numerical resolutions of 2563–40963 mesh points and with two distinct driving mechanisms, solenoidal (divergence-free) driving and compressive (curl-free) driving. By doing so, the scientists ran the world's largest simulation of supersonic turbulence on GCS supercomputers. 

Astrophysics

Principal Investigator: Christina Korntreff , Jülich Supercomputing Centre (Germany)

HPC Platform used: JUROPA of JSC

Local Project ID: hbn30

Using hydrodynamic simulations on GCS supercomputers, a team of scientists from the Max-Planck-Institute for Radio Astronomy, Bonn, and Jülich Supercomputing Centre investigated the influences of surrounding gas onto the period of binary systems. 

Astrophysics

Principal Investigator: Ewald Müller , Max-Planck-Institut für Astrophysik, Garching/München (Germany)

HPC Platform used: SuperMUC of LRZ

Local Project ID: pr86la

Massive stars end their lives as core-collapse supernovae when the stellar core implodes to a neutron star and the stellar envelope is expelled. Using computer models, we have simulated the mixing processes occurring during the explosion without assuming any symmetry.

Astrophysics

Principal Investigator: Frederic Bournaud , CEA Saclay (France)

HPC Platform used: SuperMUC of LRZ

Local Project ID: pr86di

Numerical simulations are a crucial tool to understand the physics of gas turbulence and star formation: there is no analytic theory. More than six decades of spatial scales need to be described, which is best done with "adaptive resolution" codes on supercomputers. 

Astrophysics

Principal Investigator: Wolfram Schmidt , Institut für Astrophysik, Universität Göttingen (Germany)

HPC Platform used: SuperMUC of LRZ

Local Project ID: pr74bi

One of the cutting-edge problems in current astrophysical research is the formation and evolution of galaxies similar to our Milky Way Galaxy. 

The JUPITER Research and Early Access Program (JUREAP) ensures the smooth launch of the first European exascale supercomputer, JUPITER. The GCS Exascale Pioneer Call provides sufficient computing time for selected projects to enable groundbreaking research for the German scientific community. The computing time resources are twofold: On the one hand, this call distribute resources after JUPITER is officially operational until the end of October 2025. On the other hand, it guarantees successful projects early access to JUPITER during build-up.

 

Here, we present the list of the approved GCS Exascale Pioneer Projects.

Understanding lubrication at extreme conditions is key to efficient, sustainable mechanical systems. In this context, this project deals with nanoscale lubrication, revealing how molecular dynamics simulations guide better models for friction and lubrication. Breakthroughs include a novel viscosity-pressure relationship for hydrocarbons, a lubrication model with improved boundary slip laws, and molecular insights into lubricant behavior, offering transformative tools for engineering high-performance machinery.

Fluids are key agents in many geological processes of the Earth's crust and upper mantle. Despite their importance for geological and technological processes, their thermodynamic and physical properties are not well constrained at many of the relevant conditions, especially in the supercritical state. In this project, we collaborate with experimentalists and thermodynamicists to study properties of hydrothermal fluids in a wide range of densities and temperatures. The main goals of the simulations are the development of molecular structure models including electronic and vibrational properties and prediction of thermodynamic properties such as solute dissociation constants and partial molar volumes.

The air stream in a gas turbine is firstly compressed and delivered to the combustion chamber, where fuel is mixing in and burnt, releasing a tremendous amount of heat. The hot turbulent bumt gases expand through the turbine placed downstream and the exhaust nozzle. Over the last decades, the turbine inlet temperature has increased because this leads to a higher efficiency of the gas turbine. The temperature of the hot gas of the combustion chamber (2,200 °C) and turbine section (1,700 °C) surpasses the material's maximum temperature limit (900 °C). In order to safeguard the metal walls from damage, they are covered by a ceramic thermal barrier coating (TBC) but this is not sufficient to protect the metal components from overheating.

Experimental advancements within the last two decades have enabled unprecedented control of quantum systems, posing outstanding challenges for their theoretical description. Our project is based on a novel computational strategy at the intersection of machine learning and quantum physics, utilizing artificial neural networks to efficiently represent quantum wave functions. By leveraging supercomputing resources from FZ Jülich and the Gauss Centre for Supercomputing, we have advanced the theoretical understanding of strongly interacting systems in two dimensions, including the first demonstration of the quantum Kibble-Zurek mechanism.

The IHPms21 project combined advanced computer simulations with experimental techniques to develop new materials for future microelectronics that are compatible with silicon technology. Using ab initio density functional calculations, the project helped interpret experimental results and guide further research. The team achieved three key outcomes: first, they explained why germanium surfaces behave differently depending on their orientation during graphene growth; second, they uncovered how multilayer hexagonal boron nitride (hBN) can grow from an inert gas, despite its chemical inactivity; and third, they analyzed photoemission spectra to reveal the presence of ultrathin β-Ga2O3 films on the surface of ZnGa4O4 crystals.

Liquid metals like Gallium (Ga) are a promising platform for catalytic devices such as SCALMS (Supported Catalytically Active Liquid Metal Solutions). Ga develops an oxidized surface layer (Ga₂O₃), which is known to have a major impact on droplet dynamics and technological performance.

We simulate droplets via a coupled Immersed Boundary Lattice-Boltzmann (IBLB) method, for which we introduce a generalized model for elastic properties of the membrane, to cover properties of oxidized droplets and beyond [1].

Membrane topology transformations – such as scission, fusion, and pore formation – are driven by membrane tension, curvature stress, and lipid dynamics, playing critical roles in exocytosis and organelle division. The final stage of cellular compartment division involves the scission of a highly constricted membrane neck. Using self-consistent field theory (SCFT), we explore the mechanisms of scission in single- and double-membrane neck structures.

 

The Earth is a dynamic system, various physical processes lead to deformations of the Earths surface or mass transport in its interior. Quantifying the changes with the help of measurements is a key task of geodesy, for instance to make signals attributed to climate change visible. For this, reference systems and reference surfaces are required. E.g. sea level rise refers to the so called Mean Sea Surface (MSS), or the geoid as an equipotential surface is required to show mass transport. The reference surfaces can be determined from hundreds of million measurements collected by satellites. Due to the characteristics of the collected data and the complexity of the surfaces, high performance computing is required for the numerical analysis

In July 2021, a devastating flood hit Central and Western Europe, causing severe damage, especially in the Ahr region in Germany. Researchers at the University of Bonn investigated the role of soil moisture in intensifying this extreme event. Using the JUWELS supercomputer at Forschungszentrum Jülich, they simulated varying soil moisture conditions to assess its impact on precipitation. The findings suggest that land surfaces contributed significantly to the heavy rainfall, with potential for even more precipitation under wetter soil moisture conditions. These insights can help to improve understanding of land-atmosphere interactions and disentangle drivers of extreme events.

Irradiation modeling is a crucial aspect of integrated photovoltaic (PV) system yield prediction and optimizations of the design and dimensions of PV systems. Shading models typically rely on high-resolution topography data that includes buildings and vegetation. However, the cost and limited availability of such data pose significant challenges. In this project, we took an innovative approach by considering an alternative source of topography data: maps. Specifically, we focused on the OpenStreetMap (OSM) due to its open-data availability. While maps cannot be directly used for irradiance modeling, we explored novel approaches to overcome this challenge.

Bulk metallic glasses (BMGs) are known to have remarkable mechanical properties, such as high tensile strength, elasticity, and yield strength, which surpass those of many crystalline and polycrystalline metals. These properties make BMGs highly promising candidates for applications requiring materials that can withstand high and complex mechanical stress. However, BMGs have drawbacks; they show strain softening, resulting in localized deformation in the form of shear transformation zones that later lead to the formation of shear bands. This strain-softening characteristic limits their broader application potential, as it can lead to surface defects and, ultimately, fracture.

Prof. Dr. Holger Gohlke and Jesko Kaiser investigated the binding of resensitizers in the nicotinic acetylcholine receptor as a potential treatment option for nerve agent poisoning. They identified a potential allosteric binding site, explaining the experimentally observed effect on the receptor. Based on these results, the researchers identified novel analogs with improved properties and new lead structures with improved affinity compared to MB327, potentially acting as new starting points to ultimately close the gap in nerve agent poisoning treatment.

Chiral amines, a group of small chemicals, are central building blocks to a variety of fine chemical products. These include agrochemicals and pharmaceuticals such as Sitagliptin, a potent drug used to treat type II diabetes. Accordingly, biotech and pharmaceutical companies are highly interested in the efficient and sustainable production of these compounds. A group of enzymes already in use to fill this need are Transaminases (TAs). In this project, Prof. Dr. Gohlke and Steffen Docter investigated the thermal unfolding behavior of two sets of TA variants of fold type I and IV families of PLP-dependant enzymes by simulating rigid cluster decompositions using Constraint Network Analysis (CNA).

The efficiency of any opto-electronic device, such as a solar cell, light emitting diode, or photodetector, is intrinsically linked to the nature of the electronic quantum states of the photoactive material. For a deeper understanding and targeted development of new devices, an improved theoretical description of bound electronic excitations, i.e., excitons and trions, is crucial.

The exponential scaling in computational cost when simulating quantum many-body systems poses a significant challenge to their understanding. Variational methods, that approximate the state of the quantum system using an ansatz function, promise to lower that computational cost while making no approximations about the interactions that occur in the system. A novel type of ansatz function, which we explore thoroughly in this project, uses neural networks to approximate the state of the system. As is usual in deep learning, we rely heavily on the use of GPUs during execution, thereby making the use of the JUWELS Booster Module a necessity.

Researchers from Heinrich Heine University Düsseldorf have investigated the interaction of high-intensity laser pulses with matter using particle-in-cell simulations. Their research has led to a novel mechanism for compact ion acceleration, a method to generate spin-polarized ion beams, and a potential path to probe quantum electrodynamics.

We are all made of atoms, different types of atoms, and different combinations of them, which, by their turn, are composed of a cloud of electrons and a nucleus. A nucleus contains at least one proton in its simplest form, the hydrogen atom. Comprehending the proton, the origin of its measured properties, like its mass and electric charge, and its structure is, thus, one of the most important endeavors of the physical sciences. How can we probe/see the proton and its structure?

In a breakthrough that could revolutionize particle accelerators, scientists have discovered how to better control high-energy electron beams using ultra-powerful lasers. This new understanding delves deep into the complex dance between intense laser pulses and the plasma they create, revealing the subtle mechanisms that influence electron beam stability.

The mitigation of aircraft noise is a major goal of the society to reduce the harmful effects on the human health and cognitive performance when exposed to a pervasive noise level. Although several events in the past have temporarily reduced the air traffic, a long-term constant growing rate if 4- 8%p.a. of passengers has been observed in the recent years. New concepts for on-demand Urban Air Mobility evolve and thus, additionally implying an extension of urban areas. Coping with this trend, the ACARE 1 defined ambitious goals of Europe’s vision for aviation for the year 2050 in the Flightpath 2050.

Reducing drag in engineering type flows is of paramount importance. Various approaches and configurations were tackled in the past, mostly, however, dealing with incompressible flow. In this project, researchers at University of Siegen were investigating a specific oscillation type control method for sub- and supersonic channel flow, a configuration where the fluid domain is restricted by cooled lower and upper walls.

Periodically driven ultracold atomic systems can be used to engineer topologically nontrivial phases, and can give rise to anomalous topological phases with chiral edge modes in the presence of a trivial bulk (AFTI). We have investigated the role of additional quenched disorder and two-particle interactions on this state. Within an exact diagonalization study we have found signatures of the anomalous Floquet Anderson insulator (AFAI) phase within an experimentally realized model by calculating several indicators [1]. It supports quantized charge pumping through chiral edge states, while the bulk states remain completely localized. Moreover, we have developed an efficient new algorithm for calculating higher-order Chern numbers [2].

The internal structure of the proton and neutron (collectively known as the nucleon), which form the building blocks of atomic nuclei, still poses many open questions. Not only is it not completely understood how the nucleon’s spin and momentum are composed of those of its constituent particles (the quarks and gluons), but even its size is subject to significant uncertainty arising from discrepancies between different determinations: there is a decade-old inconsistency between the electric charge radius of the proton as obtained from scattering experiments in good agreement with the value from hydrogen spectroscopy on the one hand, and the most accurate determination from the spectroscopy of muonic hydrogen on the other. This significant…

Life at the molecular level is driven by the interplay of many biomolecules. Much like man-made machines in everyday life, they need to move, rotate, react to signals or use and provide resources. Unlike man-made machines, however, they function at the atomic level so directly observing their workings is impossible as they are invisible both to the naked eye and regular optic microscopes. Specific highly specialised equipment can provide insight into the inner working of these atomic-sized machines, but such equipment is very expensive and the required wet-lab setups can be highly involved.

The helicopters & aeroacoustics group of the Insitute of Aerdynamics and Gasdynamics at the University of Stuttgart continues to develop their well-established and validated rotorcraft simulation framework. In addition to vibration prediction, noise reduction, and maneuver flight developments, new application areas like air taxis and distributed propulsion emerge out of industrial needs and fundamental research questions.

A direct numerical simulation (DNS) of a turbulent air/ethylene jet was conducted to further understand soot oxidation at relevant combustion regimes. The DNS computational domain comprised 1.5 billion points, which integrated a detailed soot model and a chemical kinetic mechanism that involved 41 chemical species. Diverging from previous works primarily focused on soot formation, this project investigates the later stages of soot evolution, particularly its oxidation in turbulent flames. The roles of OH radicals and molecular oxygen in oxidizing soot particles, along with their distribution across mixture fraction space, were analyzed. Leveraging the dataset, an assessment of existing subfilter models for soot-gas phase interaction,…

FitMultiCell is a computational pipeline developed by Prof. Dr.-Ing. Jan Hasenauer's team to tackle the complexity of simulating and fine-tuning biological tissues. This tool streamlines the creation, simulation, and calibration of biological models that imitate cellular interactions within tissues. The pipeline offers a user-friendly platform for researchers to conduct analyses on supercomputers like JUWELS. FitMultiCell's flexibility and power are demonstrated in studies on viral infections, tumor growth, and organ regeneration, proving its efficiency in refining models to match experimental data. Furthermore, enhancements for handling data outliers and scalability ensure FitMultiCell's robust application in diverse research fields.

Lattice Boltzmann method (LBM) with phase-field model has been performed to investigate the growth habit of a single ice crystal. Given the multitude of growth habits, pronounced sensitivity to ambient conditions and wide range of scales involved, snowflake crystals are particularly challenging. Only few models are able to reproduce the diversity observed regarding snowflake morphology. It is particularly difficult to perform reliable numerical simulations of snow crystals. Here, we present a modified phase-field model that describes vapor-ice phase transition through anisotropic surface tension, surface diffusion, condensation, and water molecule depletion rate.

Although MD is a widely used and accepted method, there is often the need of comparison to experiment to validate the reliability of the findings. The development of supercomputers as well as the optimization of the used simulation code led to enormous changes in the achievable dimensions. Nevertheless, the simulation of an aluminum cube of 1 µm side length would contain about 60 billion atoms. To simulate this for a duration of 1 µs simulated time, with time step 1 fs, it would take 47,248 core years on a Cray XT5 system, based on the 1.49e-6 sec/atom/time step determined at the benchmark on the LAMMPS webpage [2].

The regular ups and downs of tides are phenomena obvious to any observer of the sea. Less known is that ocean tides also undergo small changes over time, for reasons that are not yet fully understood. In this project, large simulations with a global three-dimensional ocean model were performed to understand the extent to which ocean warming and the resultant increase in the vertical density structure have contributed to changes in the largest tidal wave, M2, from 1993 to 2020. Evidence was found that upper ocean warming is the leading cause for present weakening of the size of M2 across entire ocean basins. In turn, more tidal energy is currently being transferred from M2 to three-dimensional waves in the ocean’s interior.

Collisions of protons and pions are usually observed and measured in particle accelerators. Thanks to today’s powerful supercomputers we can study these elementary particles also in theory, namely based on the core principles of Quantum Chromodynamics. By simulating the fundamental quark and gluon fields on a space-time lattice not only can we investigate why protons (and pions and many other particles) emerge at all from the strong force, but also their reaction with each other, for example in an elastic collision. And sometimes such collisions bring forth entirely new, short-lived particles, like the Δ resonance. Our project is dedicated to applying the Lattice QCD method to track from fundamental quarks and gluons to the Δ particle.

A direct numerical simulation (DNS) with finite rate chemistry has been performed to investigate the influence of flame-wall interaction (FWI) on carbon monoxide (CO) emissions in very lean turbulent premixed methane flames. CO emissions are affected by the mean strain rate of the turbulent flow, the FWI, and the interactions of the flame with the recirculation zones of the flow. The CO production and consumption in the turbulent flame differ strongly from the reaction rates in a freely propagating flame.

Geothermal energy is vital for renewable power and heating. To improve project safety and efficiency, scientists employ various methods to understand subsurface processes. Monitoring earthquakes is key and reveals how the subsurface reacts to different factors. However, interpreting seismic recordings is challenging due to complex interactions and background noise. Underground processes are complicated by fluids and fault systems. Our project uses computer simulations to analyze seismic waves, enhancing our ability to pinpoint small earthquakes accurately. This helps understanding seismic triggers and reducing hazards. Additionally, we utilize recorded background noise to directly investigate subsurface structures, maximizing data insights.

Gallium oxide (Ga2O3), a transparent semiconducting oxide with a wide bandgap of around 4.9 eV, has emerged as a promising candidate for future applications in electronics (Schottky barrier diodes, field-effect transistors), optoelectronics (solar- and visible-blind photodetectors, flame detectors, light emitting diodes, touch screens), and sensing systems (gas sensors, nuclear radiation detectors) [2]. The monoclinic β phase is its most stable and studied polymorph. Compared to the bulk properties, research on its surface properties is still sparse. However, these play a crucial role in many processes and applications, such as epitaxial growth and electrical contacts.

 

Ammonia (NH3) is a versatile compound that finds applications in fertilizer production and fiber manufacturing. The industrial synthesis of ammonia relies on high temperatures and pressures, and requires large amounts of energy while emitting greenhouse gases. A more promising alternative is the electrochemical nitrogen reduction reaction (NRR), which offers a more efficient and environmentally friendly way for ammonia synthesis. Researchers from Technical University of Munich (TUM) have taken this concept further by utilizing artificial-intelligence methods to theoretically design and screen appropriate catalysts for the NRR process. These catalysts, known as single-atom catalysts, play a crucial role in facilitating the reaction.

The project "High performance computational homogenization software for multi-scale problems in solid mechanics" focusses on simulating Advanced High-Strength Steels (AHSS) using computational methods that consider the microscale grain structure. The virtual laboratory relies on high-performance computing and robust numerical methods to predict steel behavior before experimental testing. Computational homogenization reduces the number of degrees of freedom drastically and introduces natural algorithmic parallelism. The scalability of new nonlinear solution methods, as well as the FE^2 software FE2TI was demonstrated under production conditions. A complete virtual Nakajima test was performed leveraging the power of modern supercomputers.

More food for a growing human demand needs more water to produce it. Nevertheless, global water resources are limited. Without appropriate action towards a most efficient and most sustainable water use, the world will run into a severe water crisis.

In the BMBF research project ViWA we show how High Performance Computing using SuperMUC-NG can open new ways to create the necessary knowledge for action towards more efficient and sustainable water use in agriculture. Complex global crop growth simulations based on climate and environmental data show how water could globally be saved through better farm management. Comparing simulations with actual global crop growth observations using Sentinel-2 satellites create a global monitoring system,…

Halide perovskites are booming as absorber materials for solar cells: they are cheap and easy to make in the lab while delivering devices with efficiencies of converting the energy of sunlight into electricity. An interesting, more fundamental aspect of these materials regards their atomic motions, which are unusual compared to other solar materials. Researchers at the Technical University of Munich used the power of the JUWELS cluster to investigate the impact of these atomic motions on the absorption of sunlight in halide perovskites. Their findings were intriguing: the strong atomic motions do not hinder but are in actuality a beneficial feature for the efficient collection of sunlight in halide perovskites.

Heinrich Heine University Researchers use JUWELS to study reactive metabolites in their pursuit of new biotechnological applications.

A research team led by Prof. Dr. Holger Gohlke at the Heinrich Heine University of Düsseldorf is a long-time user of the Jülich Supercomputing Centre’s (JSC’s) world-class high-performance computing infrastructure. The team has recently employed JSC’s JUWELS supercomputer to study a select class of enyzmes that play an outsized role in metabolizing chemical compounds coming from outside the body.

Controlled and reversible opening and forming of chemical bonds allows to switch material properties by light or mechanical load. The controlled change of material properties to enable different functionalities is a promising route to the design of so-called programmable materials that allow the tailored control of materials functions by well-defined external stimuli. A system that can reversibly form and break bonds is the molecule anthracene. Two anthracene molecules can bind together upon stimulation by UV-light. Heating in turn leads to the release of the formed bonds and thus to the regeneration of the initial state. Here we show that mechanical forces considerably accelerate this backreaction and do not lead to irreversible bond…

The MIQS project aims at uncovering demanding orbital-based mechanisms in quantum materials driven by strong electron correlation. First-principles many-body approaches are employed to tackle the challenging electronic states in systems such as superconducting nickelates and layered van der Waals magnets. Complex electronic phases are explored on a realistic level by combining density functional theory and dynamical mean-field theory methods on an equal footing. The high computational power of the JUWELS is needed to address the intriguing many-body physics subject to a large number of degrees of freedom at different temperature scales. Predictions and fathom design routes for novel materials and architecture is an essential part of the…

Due to a warming atmosphere and ocean, accelerated melting of the Greenland ice sheet and glaciers contribute increasingly to global sea level rise. We investigate this effect by reproducing observed sea level, temperature and salinity of the northern North Atlantic Ocean in a numerical ocean model. We compared different model simulations to situ observations and satellites data. Adding realistic Greenland melting results in a better model agreement with data, especially in Baffin Bay. Our study suggests that further work should focus on improving model resolution souch that small-scale processes can be well represented.

It has been a long-standing dream in nuclear physics to study nuclei like, for instance, carbon directly from Quantum Chromodynamics (QCD), the underlying fundamental theory of strong interactions. Such an endeavor is very challenging both, methodically and numerically. Towards this goal physicists from the European Twisted Mass Collaboration and in particular the University of Bonn have investigated two- and three-hadron systems using the approach of Lattice QCD.

Using a combination of computational simulations and experiments, researchers at the Ruhr University Bochum are investigating the complex dynamics that govern how cracks propagate in brittle and quasi-brittle materials, such as glass and hard rock. The work has implications both for mining and mineral extraction, as well as designing safer buildings.

A research team at the Heidelberg Institute for Theoretical Studies and Heidelberg University is using the power of high-performance computing (HPC) to better understand how collagen—the most common protein in our body—transports shock and other forces toward its weakest molecular links, giving researchers deeper insight into understanding how collagen in tendons absorbs stress and how this can prevent larger injuries

A team of researchers led by TU Ilmenau Professor Jörg Schumacher have been using the JUWELS supercomputer at the Jülich Supercomputing Centre to run highly detailed direct numerical simulations (DNS) of turbulent flows at the so-called mesoscale—the intermediate range where both small-scale turbulent fluid interactions and large-scale fluid dynamics converge.

A research team led by Prof. Szabolcs Borsányi, long-time users of Gauss Centre for Supercomputing (GCS) resources, have leveraged GCS’s world-class computing resources in pursuit of furthering our understanding of the most fundamental building blocks of matter and their respective roles in how the universe came to be.

A research team based at the University of Wuppertal has benefited from generous shares of Gauss Centre for Supercomputing (GCS) resources. Participating in many consortia involved in gaining a fundamental understanding of the universe’s most basic building blocks, the team combines numerical theory with experiment in pursuit of a richer understanding of how the universe and all that is in it came to be.

With the help of world-class supercomputing resources from the Gauss Centre for Supercomputing (GCS), a team of researchers led by Prof. Zoltan Fodor at the University of Wuppertal has continued to advance the state-of-the-art in elementary particle physics.

A research team led by Prof. Holger Gohlke at the Heinrich Heine University Düsseldorf is using supercomputing resources at the Jülich Supercomputing Centre (JSC) to better understand so-called hyperpolarization-activated cyclic nucleotide–gated (HCN) channels, which serve as crucial ion channels in the membrane for controlling electric pulses in the brain and heart, among other fundamental processes in the body.

Particle accelerators are among the world’s most effective methods for experiments in materials science and physics. High-intensity, laser-based accelerators are novel accelerator-concepts which are much more compact compared to conventional accelerator facilities. As next-generation facilities with even more powerful lasers begin to come online, researchers must reckon with how these devices can alter plasmas contained in these accelerators through so-called quantum electrodynamic (QED) effects. Researchers predicted how lasers in these facilities would behave, and researchers are now leveraging high-performance computing (HPC) to model these QED effects and compare with experimental data.

Using HPC resources at the High-Performance Computing Center Stuttgart, a team led by Prof. Britta Nestler is improving sintering—a common process in advanced manufacturing.

The Scalasca project brings together HPC experts in pursuit of new ways to measure and improve performance for increasingly large, heterogeneous architectures.

Using the JUWELS supercomputer at the Jülich Supercomputing Centre, researchers are simulating the so-called Brout-Englert-Higgs mechanism, or how elementary particles acquire mass.

For decades, researchers have turned to the twin power of state-of-the-art particle accelerator facilities and world-class supercomputing facilities to better understand the mysterious world of subatomic particles. These particles are very short lived and are hard to detect with even the most advanced technologies. In recent years, researchers have used the Large Hadron Collider at CERN, among other facilities, to discover new charmed baryons.

Using the JURECA supercomputer, a team of University of Cologne researchers led by Prof. Dr. Joseph Kambeitz is simulating biophysical processes in the brain in pursuit of better understanding what leads to schizophrenia in patients. 

A research team led by Prof. Dr. Holger Gohlke and Dr. Carlos Navarro-Retamal have been using the JUWELS supercomputer at the Jülich Supercomputing Centre (JSC) to better understand how plants respond to changes in their environment at a molecular level. Specifically, the team used JUWELs to simulate how the TPC protein—also prevalent in the human body—helps facilitate information sharing between different parts of a plant in responding to changes in temperature, light, or other conditions that can affect growth. In order to gain a fundamental understanding of the process, the researchers ran computationally intensive molecular dynamics simulations of up to 600,000 atoms.

A team of researchers led by Prof. Dr. Holger Gohlke and Till El Harrar have been using high-performance computing (HPC) resources at the Jülich Supercomputing Centre (JSC) to better understand how aqueous ionic liquids and seawater interact with enzymes relevant for a host of biotechnological applications. Recently, the team focused on how aqueous ionic liquids—reminiscent to molten salts, certain types of mineral-rich hydrothermal waters and the like—impact behavior of the enzymes Lipase A from Bacillus subtilis. The team published three papers on its research.

A team of researchers led by Prof. Xiaoxiang Zhu at the Technical University of Munich are using high-performance computing resources at the Leibniz Supercomputing Centre to create the first-ever 3D/4D dataset on urban morphology of settlements, joining traditional remote sensing data with social media content.

Polymers are a broad class of materials: From nylon and rubber to materials for advanced material design, polymers are long chains of repeating units. Diblock copolymers consist of two halves that repel each other and self-assemble into different phases, creating shapes such as cylinders at the molecular scale. These cylinders arrange as parallel within an individual grain but, on large scales, there are multiple grains that differ in the orientation of their cylinders.

A group of researchers from the Fritz Haber Institute and Aarhus University in Denmark have leveraged the power of the JUWELS supercomputer at the Jülich Supercomputing Centre (JSC) to develop a machine learning algorithm that helps predict how specific molecules bind to the surface of a catalyst. Catalysts play an essential role in many chemical processes, and how specific molecules interact with these materials can influence the efficiency, effectiveness, and safety of chemical reactions at an industrial scale.

A research team led by Prof. Frithjof Karsch at Bielefeld University has been using the JUWELS supercomputer at the Jülich Supercomputing Centre (JSC) as part of the international HOTQCD collaboration to better understand the conditions under which particles made of protons, neutrons, and pions go through phase transitions, and how those changes impact the system’s behavior and give rise to new forms of matter, such as quark-gluon plasma.

Using high-performance computing (HPC) resources at the Jülich Supercomputing Centre, a team of researchers led by Technical University of Dortmund Professor Frithjof Anders is gaining a better understanding of electrons’ behaviors in so-called quantum dots.

In this long term project, which lasted for 6 years and had two stages, we computed the leading order hadronic vacuum polarization contribution to the anomalous magnetic moment of the muon, aLO−HVP, using lattice quantum field theory.

Simulations at the atomic level employing supercomputers are a powerful instrument to probe and design new materials or understand better already known ones. Such findings and discoveries could potentially lead to new landscapes regarding the usage of our planet’s resources. Here we present results of our quantum-mechanical simulations investigating the band gap properties of silicon (Si) and germanium (Ge) alloys with cubic and tetragonal symmetries. These alloys’ investigations were inspired by their synthesis using high pressure techniques.

This project collates individual applications from the Fluid Dynamics group at Duisburg-Essen University. The subprojects include the investigation of phenomena from the fields of nanoparticle synthesis, supersonic flows and stratified burners using LES, DNS and direct chemistry.

Quarks are the constituents of the massive basic building blocks of visible matter. These building blocks are the hadrons, more precisely protons and neutrons, which are about 2,000 times heavier than electrons. It was a heroic effort to determine the nature of the phase transition for physical quark masses, which our group carried out in 2006 and published in Nature. This finding has fundamental consequences for the early universe and for the possible remnants we might detect even today. Note however, the result has not only relevance for the early universe (Big Bang) but also for heavy ion collisions (Little Bang), which are carried out at the RHIC (Brookhaven, USA) and LHC (Geneva, Switzerland) accelerators.

Zinc ions have shown antiviral properties, but a key issue for their use for antiviral therapy is its difficulty, as a divalent metal ion, to cross the cell membrane and thus reach its targets inside the cell. A variety of ligands, including the FDA approved drug chloroquine (CQ), form complexes with these ions have been proposed to assist zinc permeation, possibly promoting the combined beneficial action of both zinc ions and the drugs against the virus. Here, we studied the permeation of chloroquine and the interaction of the drug with zinc ions in aqueous solution. For the latter, we take advantage of highly scalable ab initio molecular dynamics simulations to explore the diverse coordination chemistry of zinc ions.

The defining properties of our numerical research in the domain of correlated electron systems are the notions of emergence and criticality. Emergence only occurs in the thermodynamic limit where the volume of the system is taken to infinity at constant particle number. To investigate this phenomena we use the Algorithms for Lattice Fermion implementation of the auxiliary field quantum Monte Carlo algorithm that allows to simulate a large variety of model systems on importance in the solid state.

 

A tremendous variety of physical phenomena involve turbulence, such as the dynamics of the atmosphere or the oceans, avian and airplane flight, fish and boats, sailing, heating and ventilation, and even galaxy formation. Turbulent flow is characterized by chaotic swirling movements that vary widely in size, from sub-millimeter, over the extent of storm clouds, to galactic scales. The interaction of the chaotic movements on different scales makes it challenging to simulate and understand turbulent flows. Turbulent thermal convection plays an important role in a wide range of natural and industrial settings, from astrophysical and geophysical flows to process engineering.

Turbulence has been a topic of research for many decades and finds its applications in many aspects of life. Still, turbulence is not fully understood up until today. The Navier-Stokes equations, which are used to describe the motion of viscous fluids, do not have a general analytical solution. Consequently, many researchers work with specific canonical cases to understand turbulence better. In the recent years as computers became increasingly powerful, more and more direct numerical simulations (DNS) have been conducted to solve turbulent flows. The goal of this project is to validate the scaling laws that have been derived for a turbulent round jet using Lie-symmetry analysis with numerical data.

Complex magnetic textures and localized particle-like structures on the nanometer scale such as chiral magnetic skyrmions with non-trivial topological properties are nowadays the most studied objects in the field of nanomagnetism. They offer the promise of new data storage and data processing technologies ranging from racetrack memories to memristive switches for neuromorphic computing. In this project we extend the realm of DFT calculations for magnetic systems to significantly larger setups from the basic description provided by our ab-initio code FLEUR

In high-precision low energy particle physics experiments, small but significant discrepancies have been found when compared to expectations from theory. This has substantially increased the interest in precision nucleon structure measurements. Theoretically, strong interaction phenomena are governed by Quantum Chromodynamics (QCD), and at energy scales relevant to the proton structure, the only known way of studying QCD from first principles is via large scale simulations using the lattice formulation.

The recent rise of renewable energy sources is promoting the use of hydrogen as a carbon-free energy carrier. One possibility to harness the energy stored in hydrogen is its usage in thermochemical energy conversion processes such as in gas turbines, industrial burners, or internal combustion engines. However, in contrast to conventional fuels, lean hydrogen/air flames are prone to thermodiffusive instabilities, which can substantially change flame dynamics, heat release rates, and flame speeds. To improve prediction capabilities of Large Eddy Simulations of hydrogen/air flames, detailed data of such flames are needed for model development and validation. However, only rare data of three-dimensional thermodiffusively unstable flames that…

Quantum Chromodynamics (QCD) is the sector of the standard model describing the strong nuclear force, which binds quarks and gluons inside hadrons. The theory confines these constituents, which are never observed directly in experiment. In this project the researchers study charmonium, a system containing a charm quark-anti-quark pair.

As a carbon-free fuel, hydrogen has the potential of emerging as the leading energy carrier for next-generation, zero-carbon power generation, and hence has received considerable attention. Hydrogen can offer significant benefits over hydrocarbon fuels, such as wide flammability range, low ignition energy, and high diffusivity. However, the use of hydrogen in gas turbines poses considerable challenges, such as the risk of flashback due to its high flame speed, which adversely affects the performance of hydrogen combustion. Flashback, a problem that occurs in premixed combustors, is the upstream propagation of the flame from the combustor into the premixing tube due to the change in mass flow rate, which could change the combustion process…

G-protein coupled receptors (GPCRs) are membrane proteins that transmit the effects of extracellular ligands to effect changes in the intracellular G-protein signaling system. Approximately 800 GPCRs are encoded in the human genome and approximately half of all marketed drugs target GPCRs. Crystal structures often deviate from the natural system: Proteins, especially membrane-bound ones, do not necessarily crystallize in their biologically active structures and the measures needed to obtain suitable GPCR crystals tend to increase the diversity between the natural environment and the crystal. It is within this context that molecular-dynamics simulations play a special role in GPCR research as a full-value complement to experimental studies.

Today, car manufacturers rely on CFD tools to design and optimise spark-ignition engines. However, current models of turbulent combustion—which are built based on the assumptions of the flamelet regime—lose their predictivity when used to simulate a highly diluted or ultra-lean combustion involving high turbulent intensities. Yet the combustion in a diluted boosted spark-ignition engine shifts from the flamelet to the thin reaction zone (TRZ) regime. This research project performed direct numerical simulations of premixed C8H18/air statistically flat flame interacting with a turbulent flow field. Results were analysed to develop a combustion model suitable for combustion in the TRZ regime based on the formalism of the coherent flame model.

Solution crystallization and dissolution are of fundamental importance for science and industry. In this project, molecular dynamics simulations were used to study these processes at the molecular scale. By following the motion of molecules towards and away from the crystal surface over short periods of time the intrinsic kinetic behavior that governs the growth and dissolution can be extracted. The obtained information is then used for parametrization of other methods such as kinetic Monte Carlo and continuum simulations to study the dynamics of the crystal surface from the nanoscale up to the microscale and beyond, where the theoretical results would be industrially relevant and easily comparable to experimental results.

Deceleration of a supersonic flow in a channel by shocks and interaction with the turbulent boundary layer leads to the formation of a complex array of shocks, subsonic and supersonic regions, and recirculation zones. In this project, high-fidelity and well-resolved large-eddy simulations (LES) of such a fully turbulent (Reδ≈105) pseudo-shock system were performed and compared with experimental data. Particular attention is paid to the occurrence of flow instabilities (such as shock motion, shock-boundary layer interaction, and symmetry breaking of the shock system), mixing behaviour in the transonic shear layer, and a comparison with sophisticated RANS turbulence models.

Leveraging the computing power of HPC systems SuperMUC and SuperMUC-NG hosted at LRZ, researchers of the Munich University of Applied Sciences investigated the piezoelectric properties of ferroelectric hafnia and zirconia, which represent a novel material class based on the fluorite crystal structure. If properly doped, such thin films show large strain effects in field induced phase transitions. A large number of doped supercells were investigated with density functional theory to find the most appropriate dopants.

Micro-scale directional grooves with spanwise heterogeneity can induce large-scale vortices across the boundary layer, which is of great importance to both theoretical research and industrial applications. The direct numerical simulation approach was adopted in this project to explore flow structure and control mechanism of convergent-divergent (C-D) riblets, as well as the impact of their spacing, wavelength and height. The results show that the C-D riblets produce a well-defined secondary flow motion characterised by a pair of weak large-scale counter-rotating vortices. This roll mode can play a key role in supressing separation when the flow undergoes adverse pressure gradients, but it may also lead to the increase of friction drag.

Clouds and precipitation are the major source of uncertainty in numerical predictions of weather and climate. A common analysis of polarimetric radar observations and synthetic radar data from numerical simulations provides new methods to evaluate models. Using the Terrestrial Systems Modeling Platform, researchers conducted ensemble simulations for multiple summertime storms over north-western Germany. The simulated cloud processes were compared in the radar space using a forward operator with the measurements from X-band polarimetric radars. In addition, sensitivity studies were conducted using different background aerosol states and land cover types in the model to better understand land-aerosol-cloud-precipitation interactions.

Geological processes are generally quite complex and occur under a wide range of thermodynamic conditions. The structure and the properties of crystalline and non-crystalline phases in the Earth’s interior are often not accessible directly and must be investigated by experiments and by numerical simulations. In this project, we use predictive molecular simulation approaches to establish relations between structural properties of relevant phases, in particular oxide and silicate glasses and melts and aqueous fluids, at high temperatures and high pressures and their respective thermodynamic and physical properties.

The open-source software framework waLBerla provides a common basis for stencil codes on structured grids with special focus on computational fluid dynamics with the lattice Boltzmann method. Other codes that build upon the waLBerla core are the particle dynamics module MESA-PD and the finite element framework HYTEG. Various contributors have used waLBerla to simulate a multitude of applications, such as multiphase fluid flows, electrokinetic flows, phase-field methods and fluid-particle interaction phenomena. The software design of waLBerla is specifically aimed to exploit massively parallel computing architectures with highest efficiency.

Most biological functions are mediated by conformational changes and specific association of protein molecules. Atomistic simulations are ideal to study the molecular details of such systems. However, often the associated timescales are beyond the maximum simulation times that can be reached even on supercomputers. In this project, researchers developed and tested advanced sampling simulations to accelerate protein domain motions and association of partner molecules. These techniques allow to study domain motions and association of protein molecules on currently accessible time scales. They were successfully applied to study the Hsp90 chaperone protein and to several protein-protein and protein-peptide systems of biological importance.

The conformations of ubiquitin chains are crucial for the so-called ubiquitin code, i.e. the selective signaling of ubiquitylated proteins for different fates in the eukaryotic cellular system. Extensive molecular dynamics simulations at two resolution levels were carried out for ubiquitin di-, tri- and tetramers of all possible linkage types. Analyzing the resulting, exceedingly large high-dimensional data sets was made possible by combining highly efficient neural network based dimensionality reduction with density based clustering and a metric to compare conformational spaces. The so obtained conformational characteristics of ubiquitin chains could be correlated with linkage-type and chain-length dependent experimental observations.

In the framework of the ASCETE (Advanced Simulation of Coupled Earthquake and Tsunami Events) project, the computational seismology group of LMU Munich and the high performance computing group of TUM jointly used the SuperMUC HPC infrastructures for running large-scale modeling of earthquake rupture dynamics and tsunami propagation and inundation, to gain insight into earthquake physics and to better understand the fundamental conditions of tsunami generation. The project merges a variety of methods and topics, of which we highlight selected results and impacts in the following sections.

Rotorcraft are regularly operating in ground effect over moving ship decks or on hillsides. However, only a very limited amount of research has been done to investigate the complex three-dimensional flow fields in these flight conditions and the resulting changes in rotor performance. Therefore, a hovering rotor in non-parallel ground effect was simulated in this project. URANS CFD simulations were made using various turbulence models to gain insight into the three-dimensional flow field, the rotor tip vortex evolution and the velocity distribution in the rotor plane. Best agreement with available experimental data was seen with a Reynolds stress model. Overall, the flow field was most affected close to the rotor hub and on the uphill side.

The ExaHyPE SuperMUC-NG project accompanied the corresponding Horizon 2020 project to develop the ExaHyPE engine, a software package to solve hyperbolic systems of partial differential equations (PDEs) using high-order discontinuous Galerkin (DG) discretisation on tree-structured adaptive Cartesian meshes. Hyperbolic conservation laws model a wide range of phenomena and processes in science and engineering – together with a suite of example models, an international multi-institutional research team developed two large demonstrator applications that tackle grand challenge scenarios from earthquake simulation and from relativistic astrophysics.

Quantum Chromodynamics (QCD) is the theory of strong interactions. It explains how quarks and gluons form the composite particles called hadrons which are observed in nature. Hadrons can be studied by means of computer simulations of QCD discretized on a Euclidean lattice. This project focuses on hadrons formed by heavy quarks. The question addressed is the relevance of including virtual charm-quark effects in lattice QCD simulations. This dynamics is challenging since it requires small values of the lattice spacing for reliable extrapolations to zero lattice spacing. It is found that its effects are at the sub-percent level even for quantities like the decay constants of charmonium at an energy scale of about half of the proton mass.

Today, large-scale computations of lattice QCD can easily reach a precision of 1% and below—a level at which it is necessary to factor in isospin breaking arising from 1. the presence of the electromagnetic interaction, and 2. the mass difference between up and down quarks. The most prominent consequence of these effects is the mass difference of the neutron and the proton as its numerical value influences the stability of matter: were this difference a bit different from what is measured in experiments, matter would become unstable so that no atoms, molecules and more complex structures could be formed. It was successfully demonstrated that this mass difference can be computed in a common lattice framework of a full QCD + QED calculation.

Investigations of different approaches to transition modelling on rotors were undertaken, including comparison to experimental data and results of other European CFD codes. For flows at Reynolds numbers below 500,000 the transition transport models predict unphysically large areas of laminar flow compared to the experimental data. A new boundary layer transition model was developed to improve the transition prediction for a wide range of parameters crucial to external aerodynamics. The new model was implemented into the DLR TAU code and works on either structured or unstructured grids. The agreement of the new model with the experimental data is significantly improved compared to the results of the basic transition transport model.

This project focuses on the modelling and physical understanding of 3D turbulent natural convection of non-Newtonian fluids in enclosures. This topic has wide relevance in engineering applications such as preservation of canned foods, polymer and chemical processing, bio-chemical synthesis, solar and nuclear energy, thermal energy storages. Different aspects of non-Newtonian fluids have been analysed in the course of this work: The behaviour of yield stress fluids in cubical enclosures, 2D and 3D Rayleigh-Bénard convection of power-law fluids in cylindrical and annular enclosures and finally the investigation of Prandtl number (Pr) effects near active walls on the velocity gradient and flow topologies.

Investigating hadron structure, how the quark and gluon constituents account for the properties of hadrons (which include neutrons and protons), is challenging due to the nature of the strong interaction. However, such information is crucial for exploiting experiments that are searching for evidence of the physics that lies beyond our current understanding of particle physics (that is encapsulated in the Standard Model) as these experiments often involve protons and neutrons in some way. Hadron structure observables can be computed via large-scale numerical calculations. This project determines key quantities on a fine lattice with physical quark masses, enabling reliable results to be extracted.

Understanding the response of silicon and diamond to shear deformation is crucial to improve the performance of nanodevices and low friction coatings. Atomic length scale simulations show that the two materials differ significantly in their amorphization-mediated wear behavior: Externally applied pressure favors the wear of silicon, while it reduces the wear of diamond. For silicon, a shear-induced recrystallization process opposes amorphization. By choosing suitable orientations of two silicon crystals in the sliding contact, the combination of both phase transformations can be exploited to grow silicon crystals with nanoscale precision.

In this project the flow in partially-filled pipes is investigated. This flow can be seen as a model flow for rivers and waste-water channels and represents a fundamental flow problem that is not yet fully understood. Nevertheless, there have neither been any high-resolution simulations nor well resolved experiments reported in literature to date for this flow configuration. In this project highly resolved 3D-simulations are performed which help further understanding narrow open-channel and partially-filled pipe flows. The analysis concentrates on the origin of the mean secondary flow and the role of coherent structures as well as on the time-averaged and instantaneous wall shear stress.

Researchers of Koc University, Istanbul, performed extensive large-scale direct numerical simulations of turbulent bubbly channel flows to examine combined effects of surfactant and viscoelasticity by using a fully parallelized 3D finite-difference/front-tracking method. The insights achieved shed light for the first time on the intricate interactions of soluble surfactant and viscoelasticity in complex turbulent bubbly flows and reveal their effects on friction drag in channels. The results are expected to guide practitioners in engineering designs such as heat exchangers and pipelines. 

Synthetic or biological amphiphiles self-assemble into spatially modulated structures on the nanoscale with applications ranging from etch masks in semiconductor fabrication, over porous membranes for separation or energy applications, to the compartmentalization of living cells. Often, such systems do not reach thermal equilibrium but, instead, the structures are dictated by processing and kinetic pathways. These molecular simulations provide insight into the correlation between molecular structure and collective dynamics that alter the self-assembly. Two results are being highlighted: (i) the kinetically accessible states in the course of directed self-assembly and (ii) the kinetic pathway of the fusion of two apposing lipid membranes.

The quality of surface water typically depends upon a complex interplay between physical, chemical and biological factors which are far from being completely understood. Most practical water quality predictions for rivers or streams rely on various simplifications esp. with regards to the turbulent flow conditions. This project aims at pushing the modeling boundary further by performing massively-parallel computer simulations which resolve all scales of hydrodynamic turbulence in river-like flows, the micro-scale flow around rigid, mobile particles, and the concentration field of suspended bacteria. The data obtained helps quantifying the shortcomings of simpler currently used prediction models and will contribute to their improvement.

Combustion instabilities in rocket thrust chambers pose a serious risk for the development of future launch vehicles as they can’t be predicted reliably by numerical simulations. To better understand the interaction between the flames and acoustic waves inside a combustion chamber, this project numerically investigates the flame response to forced transversal excitation by using Detached-Eddy simulations. In a first step, the eigenmodes of a model combustion chamber are determined from an impulse-response and they are compared to experimental results. We then investigate a specific mode coupling scenario in which the oxygen injector longitudinal eigenmode is adjusted to match the dominant transversal combustion chamber eigenmode.

Focused ion beams can be used to pattern 2D materials and ultimately to create arrays of nanoscale pores in atomically thin membranes for various technologies such as DNA sequencing, water purification and separation of chemical species. Among 2D materials, transition metal dichalcogenides, and specifically, MoS2, are of particular interest due to their spectacular physical properties, which make them intriguing candidates for various electronic, optical and energy conversion applications. Findings achieved by running large-scale molecular dynamics simulations to study the response of MoS2 monolayer to cluster ion irradiation suggest new opportunities for the creation of 2D nanoporous membranes with an atomically thin nature.

While for the design point operation of centrifugal pumps an essentially steady flow field is present, the flow field gets increasingly unsteady towards off-design operation. Particular pump types as e.g. single-blade or positive displace pumps show a high unsteadiness even in the design point operation. Simulation results for the highly unsteady and turbulent flow in a centrifugal pump are presented. For statistical turbulence models an a-priori averaged turbulence spectrum is assumed, and limitations of these state-of-the-art models are discussed. Since the computational effort of a scale-resolving Large-Eddy-Simulation is tremendous, the potential of scale-adaptive turbulence models is highlighted.

Metal hydrides have become of great scientific interest as high-temperature superconducting materials at high pressure, with hydrogen-hydrogen interactions suspected as critical in this behavior. Here, nuclear magnetic resonance experiments and electronic structure calculations are combined to explore the compression behavior of FeH and Cu2H, and results show that within the hydrides a connected hydrogen network forms at significantly larger H-H distances than previously assumed. The network leads to an increased contribution of hydrogen electrons to metallic conduction, and seems to induce a significantly enhanced diffusion of protons.

Flows over the curved surface of wings, cars, turbine blades in gas turbines and impeller blades in pumps have curved streamlines. The influence of streamline curvature on flows, drag and also heat transfer in flows is substantial to large. However, engineering models have difficulties in correctly predicting flows over curved surfaces and our knowledge on streamline curvature influences on flows is still limited. In this project, turbulent flows in moderately to strongly curved channels are studied by highly accurate, large-scale numerical simulations fully resolving the turbulent fluid motions. These give important insights into streamline curvature influences on flows, and produce data that form the basis for better engineering models.

Impressive progress has recently been made in machine learning where learning capabilities at (super-)human level can now be produced in non-spiking artificial neural networks. A critical challenge for machine learning is the large number of samples required for training. This project investigated new high-throughput methods across various domains for biologically based spiking neuronal networks. Sub-projects explored tools and learning algorithms to study and enhance learning performance in biological neural networks and to equip variants of data driven models with fast learning capabilities. Applications of these learning techniques in neuromorphic hardware and design for their future application in neurorobotics were also included.

The aerodynamics of generic space launch vehicles, in particular the flow field at the bottom of the vehicle, at transonic conditions are investigated  numerically using hybrid RANS-LES methods. The focus of the project is the investigation of the impact of hot plumes and hot walls on the flow field. It is found that both higher plume velocities and higher wall temperatures shift the reattachment location downstream, leading to a stronger interaction of shear layer and plume. An additional contribution in the pressure spectral content is observed that exhibits a symmetric pressure footprint. The increased wall temperature leads to reduced radial forces on the nozzle structure due to a slower development of turbulent structures.

At high temperatures the nuclear matter melts into a plasma state. This phase transition is expected to have a “critical point” for systems which have increasingly more protons than antiprotons. The search for this elusive critical point on the QCD phase diagram is one of the greatest challenges in today’s high energy physics, both in theory and in experiment. The calculations of the theory at non-zero densities in supercomputers are hampered by the sign-problem. In this project multiple research tracks were pursued and the methods that deal with the sign-problem and search for signals of the critical point on the phase diagram were developed.

The need to reduce the skin-friction drag of aerodynamic vehicles is of paramount importance. Nominally 50% of the total energy consumption of an aircraft or high-speed train is due to skin-friction drag. Reducing skin-friction drag reduces fuel consumption and transport emissions, leading to vast economic savings and wider health and environmental benefits. In this project, wall-normal blowing is combined with a Bayesian Optimisation framework in order to find the optimal parameters to generate net energy savings over a turbulent boundary layer. It is found that wall-normal blowing with amplitudes of less than 1% of the freestream velocity of the boundary layer can generate a drag reduction of up to 80% with up to 5% of energy saving.

Protons are composite particles: bound states of quarks and gluons, as described by the theory of quantum chromodynamics (QCD). Using lattice QCD, we know in principle how to use supercomputers to compute various properties of the proton such as its radius and magnetic moment, however this is very challenging in practice. A major part of this project was devoted to developing and studying methods for more reliable calculations, in particular for obtaining more accurate results in a finite box and for better isolation of proton states.

As most notorious greenhouse gas, CO2 emissions prevail as high as about 364 million tons carbon with the concentration reaching over 400 ppm in the atmosphere. A drastic reduction of CO2 is urgently necessary for sustainable growth and to fight climate change. The electrochemical reduction of CO2 (CO2RR) is a promising approach to utilize renewable electricity to convert CO2 into chemical energy carriers at ambient conditions and in small-scale decentralized operation. Researchers from Technical University of Munich have employed an active-site screening approach and proposed carbon-rich molybdenum carbides as a promising CO2RR catalyst to produce methanol.

A multi-institutional team of researchers is developing a data assimilation framework for coupled atmosphere-land-surface-groundwater models. These coupled models, which potentially allow a more accurate description of the coupled terrestrial water and energy fluxes, in particular fluxes across compartments, are affected by large uncertainties related to uncertain input parameters, initial conditions and boundary conditions. Data assimilation can alleviate these limitations and this project is focused in particular on the value of coupled data assimilation which means that observations in one compartment (e.g., subsurface) are used to update states, and possibly also parameters, in another compartment (e.g., land surface).

GPCRs sit in the cell membrane and transmit signals from the outside of the cell to its interior. Currently, drugs targeting these receptors only work by mimicking ligands, i.e. they activate or inhibit the receptors by changing their conformation. If the GPCR adopts an active conformation, it can bind proteins on the intracellular side of the cellular membrane, which then transmit the signal inside the cell. In this study, we investigated how a protein that stops the GPCR from signaling, interacts with a prototypical GPCR. We discovered that specific lipids can modify how signals are transmitted by modifying the way of interaction between the GPCR and arrestin. In the future this could enable the discovery of a new kind of drugs for GPCRs.

Helicopters and other rotorcraft like future air taxis generate substantial sound, placing a noise burden on the community. Advanced simulation capabilities developed at IAG over the last decades enable the prediction of aeroacoustics together with aerodynamics and performance, and thus allow an accurate and reliable assessment of different concepts long before first flight. Consequently, this technology serves to identify promising radical configurations initially as well as to further optimize designs decided on at later stages of the development process. Conventional helicopters may benefit from these tools as much as breakthrough layouts in the highly dynamic Urban Air Mobility sector.

Hydrogen-enriched fuels can reduce the CO2 emissions of gas turbines. However, the presence of hydrogen in fuel mixtures can also lead to undesirable phenomena like flashback. Swirling combustors can take advantage of an axial air injection to increase their resistance against flashback. Such an example is the swirl-stabilized presented in experiments at the TU Berlin. The axial momentum ratio between the fuel jets and the air was found to control flashback resistance. This experimental hypothesis motivates the present study where large-eddy simulations of the combustion system are carried out to study the physics behind flashback phenomena in hydrogen gas turbine combustors.

In this project the most inner structure of the proton has been deciphered through a large-scale numerical simulation of quantum chromodynamics. This could be achieved by novel algorithms developed by the project team. In particular, the project made a large leap forward to solve the spin puzzle of the proton. While theory predicted a dominant contribution to the spin of the proton from the quarks, in experiments it was found that this contribution is surprisingly small. The research team found out that it is actually the gluon which is contributing a large fraction of the spin. Although still a number of systematic uncertainties have to be fixed, this is a most remarkable result which will lead to eventually resolve the proton spin puzzle.

Data sent over the internet relies on public key cryptographical systems to remain secure. A project under leadership of Dr. Paul Zimmermann of the French National Institute for computer science and applied mathematics (INRIA), run on HPC system JUWELS of the Jülich Supercomputing Centre, has been carrying out record computations of integer factorisation and the discrete logarithm problem, the results of which are used as a benchmark for setting the length of the keys needed to keep such systems secure.

In this project, the biophysics of Photosynthesis are probed employing high-performance computing. Photosynthesis is based on the Sun light and fuels the metabolic pathways of numerous organisms in our biosphere. However, fluctuations in the light intensity or quality are expected due to the diurnal cycle, or the environmental conditions and could be detrimental to plants. Absorption of light and tunnelling of the associated energy towards the reaction centres of the photosynthetic apparatus are finely-tuned within a well-orchestrated photoprotective mechanism. The atomic-scale details of this mechanism is probed by computational biophysics, with applications on the increase of crop yields and artificial photosynthesis.

It is well-known that the catalytic properties of metals may extend beyond their melting point. Recently, this has been exploited to grow high-quality 2D materials such as graphene. To improve our understanding of the growth mechanism on liquid metal catalysts, researchers at the Technical University of Munich have employed a multi-scale modelling approach. Here, detailed simulations of various building blocks for the final graphene sheet such as simple hydrocarbons and smaller graphene flakes on solid and liquid Cu surfaces have been carried out. The insights from these simulations were then used to propose a mesoscopic model for the dynamics of graphene growth on molten Cu based on capillary and electrostatic interactions.

Computational fluid dynamics (CFD) simulations play an important role in today’s science and technology. Therefore, it is crucial to validate its underlying methods and models. This can be done by experiments or with molecular dynamics (MD) simulations, but in some cases only the latter are applicable. Since MD simulations follow the motion of each molecule individually, they are computationally very demanding, but they rest on an excellent physical basis. In this project, large systems of several hundred million atoms are considered to study the thermo- and hydrodynamic behavior of fluids during shock wave propagation, droplet coalescence and injection. The results are compared to that of macroscopic numerical methods.

Super-Yang-Mills theory is a central building block for supersymmetric extensions of the Standard Model. While the weakly coupled sector can be treated within perturbation theory, the strongly coupled sector must be dealt with a non-perturbative approach. Lattice regularizations provide such an approach but they break supersymmetry and hence the mass degeneracy within a supermultiplet. Researchers of Uni Jena study N=1 supersymmetric SU(3) Yang-Mills theory with a lattice Dirac operator with an additional parity mass. They show that a special 45° twist effectively removes the mass splitting at finite lattice spacing–thus improves the continuum extrapolation—and that the DDαAMG algorithm accelerates such lattice calculations considerably.

Porous media are everywhere. When Reynolds numbers in pores are large, the unsteady inertial effects become important giving rise to the onset of turbulence, for instance in packed bed catalysis, gas turbine cooling, and pebble-bed high-temperature nuclear reactors. Access to detailed flow measurements is very challenging due to the inherent space constraints of the porous media. Therefore, a research group of the Institute of Aerospace Thermodynamics (ITLR) at University of Stuttgart uses high-fidelity direct numerical simulation (DNS) to investigate the physics of fluids inside porous media which serves for industrial 3D-printed porous media design.

The Universe was just a few microseconds old when the gradually cooling matter organized itself into the massive particles that form much of the visible matter today, protons and neutrons. The fascinating world of the hot Universe is recreated in huge particle accelerators. These experiments go beyond the study of the primordial world, as they can probe a whole new dimension by tuning the balance of particles vs antiparticles. This imbalance, the net baryon density, could be tuned to the extreme, as that can be found in neutron stars. Researchers of Uni Wuppertal launched a large scale simulation project to calculate how baryon density impacts temperature where today's particles emerge from the primordial plasma.

Combustion remains the most important process for power generation and more research is needed to reduce future pollutant emissions. However, combustion is governed by thermo-chemical processes that interact over a wide range of length and time scales. Detailed simulations are of high interest to gain more information about flames. Two examples of large-scale simulations of challenging flame setups are given: The thermo-diffusive instabilities of hydrogen flames as well as the interplay between turbulent flow and flames. A special method for investigating the local dynamics of flames, called flame particle tracking, has been implemented specifically for large parallel clusters for high performance computing to further evaluate these cases.

A supramolecular polymer (SMP) has functional groups which interact with each other to form a physical bond. In contrast to chemical bonds, the bond formation in an SMP is reversible and the resulting aggregate morphology in a SMP melt thermally fluctuates. For functional groups allowing only a pairwise association, a ring aggregate is highly important as a ring topologically reduces the mobility of surrounding linear polymers by threading. Using molecular dynamics simulations of SMPs, the effect of ring aggregates on the system relaxation time governing rheological response was investigated. It was shown that the presence of ring aggregates slows down rheological response as measured by a reduction of the so-called entanglement length.

Heat shock protein 90 (Hsp90) is a molecular chaperone essential for the folding and stabilization of a wide variety of client proteins in eukaryotes. Many of these processes are associated with cancer and other diseases, making Hsp90 an attractive drug target. Hsp90 is a highly flexible protein that can adopt a wide range of distinct conformational states, which in turn are tightly coupled to the enzyme’s ATPase activity. In this project, atomistic molecular dynamics simulations, free energy calculations, and hybrid quantum mechanics/classical mechanics simulations were performed on both monomeric and full-length dimeric Hsp90 models to probe how long-range effects in the global Hsp90 structure regulate ATP-binding and hydrolysis.

The nucleon has an extremely complicated many-body wave function because QCD is very strongly coupled, very non-linear and characterized by massive quantum fluctuations. Its investigation started with collinear processes and has by now progressed to non-collinear ones. The latter are characterized by non-trivial parallel transport, leading to observable effects. Many of these are described by TMDs the properties of which are not yet well understood and are planned to be studied at the new Electron Ion Collider. We have calculated one of the most important of these properties on the lattice. Only in 2020, first lattice calculations, all using alternative approaches, of this quantity were published. All results agree within error.

The SuperMUC-NG is being used to simulate materials from first-principles, materials ranging from active materials important to technology to planetary materials that govern, for example, Earth’s magnetic field. Solid and liquid iron at conditions of Earth’s core have been simulated, and transport properties such as electrical and thermal conductivity were computed to constrain the properties that govern Earth’s dynamo. At much lower pressures, filled ices, which are believed to form in the interior of water planets such as Titan, and carbon solubility in silicates melts in the mantle of the Earth were studied. Three new class of materials were developed computationally: polar metallocenes, ferroelectric clathrates, and polar oxynitrides.

Using a canonical jet in cross flow (JICF) flame configuration, researchers of TU Darmstadt performed a high-resolution DNS study concerning differential diffusion effects and mixing characteristics during hydrogen combustion. The investigations in the hydrogen JICF configurations were twofold. First, a detailed analysis of the DNS data was to yield a fundamental understanding of mixing characteristics in the JICF configuration and differential diffusion effects. Second, commonly applied tabulated chemistry approaches and their capability of predicting differential diffusion were to be validated against the DNS data. The latter, which is of highly practical interest for a related project, was the final target of this project.

Generating high energy ions by irradiating an ultra-intense laser pulse on a foil-coated foam-like double-layer plasma target is investigated with the help of particle-in-cell simulations. The foil is ultra-thin so that the incident laser pulse can penetrate through it. The acceleration of ions happens in the foam when the density of the foam is in the laser-induced relativistic-transparent regime. Simulations show that a proton beam with peak energy beyond 150 MeV is generated by using a 16 Joule laser pulse. The laser pulse used in the simulation is already available, and the targets can be prepared with the current technology. This simulation work provides helpful information for the further experiments and related applications.

The African Continent will be severely hit by climate change. A necessary building brick for counteraction are reliable projections of the African climate of our century. The CORDEX CORE initiative is designed to provide such information for the CORDEX CORE regions, among them CORDEX CORE Africa. IMK-TRO contributed to this with an ensemble of presently ten regional climate simulations performed on the Hazel Hen at HLRS Stuttgart. Results indicate dramatic changes especially in precipitation. The simulations presented here will be part of the IPCC AR6 atlas of regional climate change and the CORDEX data repository. They will be freely available for impact, adaptation and mitigation studies.

Complex I is the largest and most intricate respiratory enzyme, which couples the free energy released from quinone reduction to transfer protons across a biological membrane. Recent X-ray structures of bacterial and eukaryotic complex I have advanced our understanding of the enzyme’s function, but the mechanism of its long-range energy conversion remains unsolved. Here, we use atomistic molecular dynamics simulations and free energy calculations to study how the protonation state, hydration dynamics, and conformational dynamics of complex I regulate its proton pumping activity. Our simulations mimic transient states in the enzyme’s pumping cycle to draw a molecular picture of the protonation signals along the membrane domain of complex I.

Quantifying the dynamics of basins across diverse time and space scales is one challenge faced by earth scientists. To understand their response to natural or man-made forcing is crucial to constrain the state and fate of georesources and hazards related to their exploitation. In this project, we developed and used a hybrid scalable modelling approach combining deterministic and probabilistic modules to improve our comprehension of the complex nonlinear dynamics of this specific terrestrial compartment interacting with the other geo-hydro-atmosphere systems making up the system Earth.

The helicopters & aeroacoustics group of the Institute of Aerodynamics and Gas Dynamics at the University of Stuttgart continues to develop their well-established and validated rotorcraft simulation framework. Vibration prediction and noise reduction are currently the focus of research, and progress into manoeuvre flight situations is on the way. For two decades, high-performance computing leverged within the HELISIM project has enabled improvements for conventional helicopters as much as for the upcoming eVTOLs, commonly known as air taxis, in terms of performance, comfort, and efficiency. Community acceptance will be fostered via noise reduction and safety enhancements, made possible by this research project.

Lattice QCD enables calculation of many details of quark-gluon bound states like the proton. One first parameterizes all properties of, e.g., the proton by certain parameters and functions. Next, one links experimental observables to these quantities and clarifies their meaning. In recent years, lattice calculations have become a valid alternative to performing experiments to determine these quantities. We have calculated the quantity d2, which characterizes certain spin-dependent effects and is linked to the color force exerted on quarks in a proton or neutron. Non-trivial renormalization properties make this an especially difficult quantity to calculate, but this project was successful in doing so with results that agree with experiment.

Turbulent convection is one essential process to transport heat in fluid flows. In many of the astrophysical or technological applications of convection the working fluid is characterized by a very low dimensionless Prandtl number which relates the kinematic viscosity of the fluid to its temperature diffusivity. Two important cases are turbulent convection in the Sun and turbulent heat transfer in the cooling blankets of nuclear fusion reactors. Massively parallel simulations of the simplest setting of a turbulent convection flow, Rayleigh-Bénard convection in a layer or a straight duct that is uniformly heated from below and cooled from above, help to understand the basic heat transfer mechanisms that these applications have in common.

Recent cosmological observations tell us that only a small part of the matter content of the Universe is coming from ordinary particles, e.g. protons and neutrons. We call the rest dark matter. But what constitutes this invisible ingredient of the Universe? A possible candidate is the so called axion, for which a mass limit was worked out in the prequel of this project. To learn more on the features of this hypothetical particle its dynamics was investigated through a link to the strong interactions.

A new active surface actuation technique to reduce the friction drag of turbulent boundary layers is applied to the flow around an aircraft wing section. Through the interaction of the transversal traveling surface wave with the turbulent flow structures, the skin-friction on the surface can be considerably reduced. Highly-resolved large-eddy simulations are conducted to investigate the influence of the surface actuation technique on the turbulent flow field around an airfoil at subsonic flow conditions. The active technique, which previously was only tested in generic scenarios, achieves a considerable decrease of the airfoil drag.

This project explores laminar-turbulent transition, turbulence, and flow control in boundary layers at various flow speeds from the subsonic to the hypersonic regime. The physical problems under investigation deal with prediction of laminar-turbulent transition on airfoils for aircraft, prediction of critical roughness heights in laminar boundary layers, turbulent drag reduction, the origins of turbulent superstructures in turbulent flows, the use of roughness patterns for flow control, effusion cooling in laminar and turbulent supersonic boundary-layer flow, DNS of disturbance receptivity on a swept wing at high Reynolds numbers, and plasma actuator design for active control of disturbances in a swept-wing flow.

A large amount of the energy needed to push fluids through pipes worldwide is dissipated by viscous turbulence in the vicinity of solid walls. Therefore the study of wall-bounded turbulent flows is not only of theoretical interest but also of practical importance for many engineering applications. In wall-bounded turbulence the energy of the turbulent fluctuations is distributed among different scales. The largest energetic scales are denoted as superstructures or very-large-scale motions (VLSMs). In our project we carry out direct numerical simulations (DNSs) of turbulent pipe flow aiming at the understanding of the energy exchange between VLSMs and the small-scale coherent.

Large scale simulations are particularly valuable and important for a better understanding of coupled multi physics problems describing a large class of physical phenomena. This research project focuses on the development of new numerical methods for efficiently solving coupled non-linear and time-dependent fluid flow problems on a large scale. In particular, two applications are considered. Namely, the Navier–Stokes equations coupled to a transport equation describing diluted polymers and geodynamical model problems which involve non-linearities in the viscosity. The goal is to develop new methods for solving these problems, evaluating their performance and scalability, and to perform simulations based on these new methods.

Using DLR’s finite-volume solver TAU, researchers of the Institute for Aerodynamics and Flow Technology at DLR Göttingen numerically investigated the vortex system of four rotating and pitching DSA-9A blades. The computations were validated against experimental data gathered using particle image velocimetry (PIV) carried out at the rotor test facility in Göttingen. Algorithms deriving the vortex position, swirl velocity, circulation and core radius were implemented. Hover-like conditions with a fixed blade pitch were analyzed giving a good picture of the static vortex system. These results are used to understand the vortex development for the unsteady pitching conditions, which can be described as a superposition of static vortex states.

The availabiltiy of ultra-short, high-power lasers has led to greater interest in their potential use for accelerators, as the charge separation in plasmas can induce enormous electromagnetic field strengths on a sub-micrometer scale. With the high performance and extreme scalability of the Plasma Simulation Code (PSC) for fully kinetic simulations, a wide field of applications was researched: From ions for medical purposes (Ion Wave Breaking Acceleration and Mass-Limited Targets) to breakthrough Lepton acceleration by proton-driven wakefields (AWAKE), all the way to radiation generation (attosecond X-ray pulses from Ultra-Thin Foils). Even QED based approaches were covered in this project.

Although Quantum Chromodynamics (QCD) has long been established as the correct theory of the subatomic strong interaction, obtaining quantitative predictions from it often represents a challenging computational task. In this project, large-scale lattice QCD simulations are used to determine structural properties of protons and neutrons. The lattice approach to QCD amounts to discretizing space-time and applying importance-sampling techniques to the path-integral representation of QCD. One specific observable under scrutiny in this project is the “scalar matrix element” of the proton, which provides a quantitative answer to the question of “How much would the proton mass change if the light quark masses changed by a small amount?”.

As a Lagrangian method, Smoothed Particle Hydrodynamics (SPH) has been explored and demonstrated for a wide range of applications. Several open-source frameworks exist for the large-scale parallel simulation of particle-based methods in which the resolution of simulation is fixed. Some preliminary work has also been published to tackle the difficulties encountered in extending codes with adaptive-resolution capability. However, the support for fully parallelized adaptive-resolution in distributed systems is generally still limited in the aforementioned codes. This research project focuses on an alternative approach by introducing a new multi-resolution parallel framework employing several algorithms from previous work.

To unravel the complexity of the solid state, researchers from the University of Würzburg have mastered very different and complementary methods. Density functional theory in the local density approximation with added dynamical local interactions using the dynamical mean-field approximation has the merit of being material dependent since one can include the chemical constituents of materials. Spacial and temporal fluctuations are crucial to understand e.g. the Iridates, a topic that is explored with the new pseudo-fermion functional renormalization group. Another aspect of this research are realistic quantum Monte Carlo simulations of free standing graphene aiming to elucidate the role of electronic correlations.

The heat transfer in the stagnation region of an impinging jet at given jet to distance ratio, Re-number and Temperature ratio also depend on the turbulent inflow characteristics. Using Direct Numerical Simulations, the Nusselt-number distribution as well as the turbulent statistics close to the heated wall have been investigated. At first a calculation has been done comparing the results with published DNS and experiments from Dairay et al. (2015). Since in their paper not all necessary turbulence values were given, the missing values (e.g. turbulent length scale) had to be adjusted in order to fit their results. A good agreement has been found of our calculations with their DNS and experiments.

In order to simulate compressible multi-phase flows at extreme ambient conditions, researchers from the Institute of Aerodynamics and Gas Dynamics have developed a multi-phase flow solver based on the discontinuous Galerkin spectral element method in conjunction with an efficient tabulation technique for highly accurate equations of state. The aim of this development is the simulation of phase transition, droplet dynamics and large-scale multi-component phenomena at pressures and temperatures near the critical point. Simulations of liquid fuel injections and shock-drop interactions have been performed on the HPC systems installed at the High-Performance Computing Center Stuttgart (HLRS).

The Standard Model of Particle Physics is a highly successful theoretical framework for the treatment of fundamental interactions, but fails to explain phenomena such as dark matter or the abundance of matter over antimatter. Precision observables, such as the anomalous magnetic moment of the muon, aμ, play a central role in the search for “New Physics”. A promising hint is provided by the persistent tension of 3.7 standard deviations between the theoretical estimate for aμ and its experimental determination. In our project we employ the methodology of lattice QCD to compute the hadronic contributions to aμ from first principles. In the long run, our results will supersede the estimates based on data-driven approaches and hadronic models.

The aerodynamic flow field around helicopters is challenging to simulate due to complex configurations in relative motion. In an effort to evolve computational fluid dynamics (CFD) technology to new levels of accuracy, reliability, and parallelization efficiency, the helicopter & aeroacoustics group at the IAG of University of Stuttgart employs advanced, high-order Discontinuous Galerkin (DG) methods to help solve difficult rotorcraft-based engineering applications. Complex geometries, curved surfaces, relative motion with elaborate kinematics, and fluid-structure coupling to blade dynamics call for sophisticated techniques within the simulation tool chain to account for all important physical phenomena relevant to the field of study.

To avoid dangerous climate change, we have to reduce the emission of greenhouse gases radically. This requires – among other measures – an increase of renewable sources of energy like solar and wind. In 2019, already a quarter of Germanys electricity demand has been met by wind power. In order to increase this share, one has to develop sites in hilly terrain. High resolution models are required to assess the suitability of candidate sites with respect to turbulence intensity, power production and variability. This project supports the development of the test-site WINSENT, which is located on the Swabian Alp near Stuttgart.

The water electrolysis in Proton Exchange Membrane (PEM) cells is fitting plenty of industrial requirements. The main drawback of PEM cells however is the overpotential of the oxygen evolution reaction. In its acidic environment iridium dioxide (IrO2) is currently the only stable catalyst. Yet the low abundance of iridium makes a reduction of its loading inevitable. One approach to decrease the catalyst loading is the use of nanoparticles. For catalyst optimization a general understanding of shape and surface structure of these nanoparticles is required. In this project a protocol has been developed to generate and simulate IrO2 nanoparticles based on energies of slab models and to provide insights regarding stability and structure.

This project looked into various strategies to couple domains of distinct physical and numerical properties to tackle direct aero-acoustic simulations. A turbulent flow around an airfoil and the emitted sound waves in a large area of interest was simulated. Different physical effects can be observed in spatially separated domains and the appropriate equation systems are solved in each one using the best fitting numerical discretization. The main focus of the project was the evaluation of different coupling methods to enable this partitioned simulation on massively parallel systems.

Carbon nitride materials have attracted vast interest in the field of photocatalytic water splitting. However, the underlying mechanism is not fully understood. Herein, results are being reported from large-scale first-principles simulations for the specific electron- and proton-transfer processes in the photochemical oxidation of liquid water with heptazine-based photocatalysts. The results reveal that heptazine possesses energy levels that are suitable for the water oxidation reaction. Moreover, the critical role of the solvent in the overall water-splitting cycle is shown. A simple model is developed to describe the water oxidation mechanism.

Cells communicate with each other through biochemical as well as mechanical signals. Essential biological processes such as cell division are critically steered by the tension across the cell-cell contacts. Using extensive molecular dynamics simulations, scientists analyzed the underlying molecular principles of mechano-sensing at cell-cell contacts. These simulations can give first insights into how proteins present at the cell-cell contact change their structure and localization and thereby help to sense mechanical stimuli. The findings can help understanding the mechanisms by which tissues, e.g. skin, grow along the direction of pulling forces which were applied by adding virtual springs into the simulation system.

The general interest of the researchers of this project is in disordered thin film superconductors within the Boguliubov-deGennes (BdG) theory of the attractive-U Hubbard model in the presence of on-site disorder; the sc-fields are the particle density n(r) and the gap function ∆(r). For this case, system sizes unprecedented in earlier work are being reached. They allow to study phenomena emerging at scales substantially larger than the lattice constant, such as the interplay of multifractality and interactions, or the formation of superconducting islands. For example, it is being observed that the coherence length exhibits a nonmonotonic behavior with increasing disorder strength already at moderate interaction strength.

Using the vast computing power of the HPC system JUWELS of JSC, an international team of physicists – the HotQCD Collaboration – simulates almost massless quarks and reveals another piece in the puzzle of how hot quarks and gluons behave under extreme thermal conditions.

High-energy ion-beam therapy of tumours has many advantages compared with conventional radiation therapy. Ion beams generated by synchrotron accelerators have been used in many medical institutions. However, a synchrotron accelerator has a large footprint (soccer field size) and is expensive. With the rapid development of high power laser technology, a laser-plasma ion accelerator is a more compact (table-size) and inexpensive alternative. Ion wave breaking acceleration which happens in laser-driven foam-like plasma targets is a promising regime for designing controllable high-energy high-quality ion accelerators. To gain a deeper knowledge on it, researchers carried out 3D simulations on SuperMUC using the Plasma-Simulation-Code (PSC).

The approximately 800 G-protein coupled receptors (GPCRs) in the human genome regulate communication across cell walls. They are targeted by approximately 40% of all marketed drugs. The project uses molecular-dynamics (MD) simulations to investigate ligand binding and receptor activation processes in GPCRs. An activation index for Class A GPCRs has been developed from a series of µsec molecular-dynamics simulations and tested for 275 published X-ray structures.  Binding of the α-domain of G proteins to GPCRs has also been characterized in detail. Metadynamics simulations in conjunction with unbiased MD simulations demonstrated the effect of mutations on the GPCR-ligand interaction in the histamine H1 receptor.

Roughness of many natural and engineered surfaces follows a scaling law called self-affine scaling. In project chka18, the origins of self-affine have been investigated using Molecular Dynamics simulations. It was shown that the self-affine roughness emerges naturally during deformation of initially flat surfaces in different materials.

Numerical models are excellent tools to improve our understanding of atmospheric processes across scales since they provide a consistent 4D representation of the atmosphere. Project WRFSCALE consists of different sub-projects, applying the Weather Research and Forecasting (WRF) model at resolutions between 3 km and 100 m, performing investigations in the fields of data assimilation, bio-geoengineering and boundary layer research. By increasing the resolution to 100 m, the model starts to explicitly resolve the representation of turbulence. With such simulations and comparisons to high-resolution observations, it is the aim to better understand the turbulent boundary layer and its interaction with the underlying land surface.

The Atlantic Meridional Overturning Circulation transports warm tropical surface water towards northern Europe and returns cold water at depth to the world’s ocean. At the same time it plays a significant role in the global carbon cycle through the ocean’s ability to dissolve carbon dioxide. This overturning is thus of great climatic importance, but a complete picture of its driving forces has not yet emerged due to several observational and theoretical challenges. Using realistic coarse and high resolution ocean models, scientists investigated the ocean response to changes in wind stress and the ability of meso-scale eddy parameterisations to simulate that response.

In order to quantify the uncertainty due to stochastic input in computer fluid dynamic simulations, researchers from the Institute of Aerodynamics and Gas Dynamics developed an Uncertainty Quantification framework and applied it to direct noise computations of aeroacoustic cavity flows. Simulations have been performed with the discontinuous Galerkin spectral element method on HPC system Hazel Hen at the High Performance Computing Center Stuttgart (HLRS). The aim of this investigation is to gain insight into the sensitivity of uncertain input with respect to the acoustic results and to get a reliable comparison between numerical and experimental results.

The simulation of turbulent, partially premixed flames constitutes a challenge due to the complex interplay of the mixing process of fuel and oxidizer, chemical reactions and turbulent flow. Therefore, a detailed numerical simulation of an experimentally investigated flame of laboratory scale has been performed, which allows to study these fundamental interactions in great detail. The results have been compiled into a database which aids the improvement of future combustion simulations. The simulation has been performed with an in-house solver based on OpenFOAM, which includes several performance optimizations to maximize the hardware utilization on supercomputers.

Combustion noise is an undesirable, but unavoidable by-product of turbulent combustion in, e.g., stationary gas turbines or aeronautical engines. This project combines Large Eddy Simulation (LES) of turbulent, reacting flow with advanced System Identification (SI) – a form of supervised machine learning –  to infer reduced-order models of combustion noise. Models for the source of noise on the one hand, and the flame dynamic response to acoustic perturbations on the other, are estimated to make possible the flexible and computationally efficient prediction of combustion noise across a wide variety of combustor configurations.

In nuclear fusion experiments, researchers routinely heat hot gases up to temperatures of 100 million degrees in order to create the conditions needed for energy-producing fusion reactions. Turbulence is one of the main obstacles on the way to sustaining these conditions reliably. A particular challenge is found in the plasma edge, where turbulence is suppressed by a self-organized transport barrier. Researchers from the Max-Planck Institute for Plasma Physics have made important progress to understanding the turbulence in this region, leveraging resources provided by the Gauss Centre for Supercomputing.

New wind harnessing generators that gather energy through a phenomenon known as vortex-induced vibrations could represent a new frontier for renewable energy. Researchers of the Barcelona Supercomputing Centre have been using high-performance computing system SuperMUC of the Leibniz Supercoputing Centre to help advance this technology.

Recent direct numerical simulations in closed slender Rayleigh-Bénard convection cells advanced to Rayleigh numbers of Ra = 1015 which were never obtained before and reveal a classical turbulent transport law for the heat transfer from the bottom to the top of the cell which is based on the concept of marginally stable boundary layers. Our simulations were able to resolve the complex dynamics inside the thin boundary layers at the top and bottom plates of the convection cell and to determine a steady increase of the turbulent fluctuations without an abrupt transition near the wall for a range of 8 orders of magnitude in Rayleigh number.

Understanding the internal structure of the nucleon is an active field of research with important phenomenological implications in high-energy, nuclear and astroparticle physics. Nucleon structure functions and their derivatives, parton distribution functions (PDFs) and generalized parton distribution functions (GPDs), teach us how the nucleon is built from quarks and gluons, and how QCD works. Beyond that, the cross section for hadron production at the LHC relies upon a precise knowledge of PDFs.

Within this project, the goal is to study and develop novel approaches to boost the performance of thin film solar cells. For this, 3D optical simulation of the photovoltaic devices is performed by discretizing Maxwell’s equations. A sophisticated light management is important to construct thin-film solar cells with optimal efficiency. The light management is based on suitable nano structures of the different layers and materials with optimized optical properties. The design, development and test of new solar cell prototypes with respect to an optimal light management are time consuming processes. For this reason, suitable models and simulation techniques are required for the analysis of optical properties within thin-film solar cells.

Recent developments in direct injection systems aim at increasing the rail pressures to more than 3000 bar for Diesel and 1000 bar for gasoline, to enhance liquid break-up and mixing which in turn improves combustion and reduces emissions. Higher flow accelerations, however, imply thermo-hydrodynamic effects, e.g. cavitation, which occurs when the pressure locally drops below saturation conditions and the liquid vaporizes. The subsequent collapse of such vapor structures causes the emission of strong shock-waves leading to material erosion. But cavitation can also be beneficial by promoting primary jet break-up, thus the ability to predict cavitation and cavitation erosion during the early stages of design of fuel injectors is desirable.

The proteasome is a large biomolecular complex responsible for protein degradation. Recent experimental data revealed that there is an allosteric communication between a core and regulatory parts of the proteasome. In the project, researchers have used atomistic simulations to study molecular details of the allosteric signal – in their study triggered by a covalent inhibitor. While the inhibitor causes only subtle structural changes, the proteasome-wide fluctuation changes may explain the self-regulation of the biomolecular machine.

Aircraft wake vortices pose a potential threat to following aircraft. Highly resolving numerical simulations provide valuable in-sights in the physics of wake vortex behaviour during different flight phases and under various environmental conditions. Hybrid simulation techniques introduce the flowfield around detailed aircraft geometries into an atmospheric environment that controls the vortical aircraft wake until its decay. The vision of virtual flight in a realistic environment is addressed by the two-way coupling of two separate flow solvers. To mitigate the risk of wake encounters and thereby to improve runway capacity, so-called plate lines have been developed and tested at Vienna airport.

Using HPC system resources available at the Jülich Supercomputing Centre, scientists of the Institute for Theoretical Physics of the Goethe-Universität in Frankfurt/Germany are performing extensive simulations to theoretically predict the properties of the phase transition from nuclear matter to a quark gluon plasma state.

In order to analyse the complex flow in rotating turbomachinery components, researchers from the Institute for Aerodynamics and Gas Dynamics performed high fidelity, large-scale turbulent flow computations of stator-rotor interactions using the discontinuous Galerkin spectral element method on the HPC system Hazel Hen at the High Performance Computing Center Stuttgart (HLRS). The aim of this investigation is to gain insight into the intricate time-dependent behaviour of these flows and to inform future design improvements.

By applying approaches based on computational chemistry, researchers at the University of Marburg are addressing the challenge of designing functional materials in a novel way. Using computing resources at the High-Performance Computing Center Stuttgart, the scientists under leadership of Dr. Ralf Tonner model phenomena that happen at the atomic and subatomic scale to understand how factors such as molecular structure, electronic properties, chemical bonding, and interactions among atoms affect a material's behaviour.

Project DEFTD is focused on large scale computer simulations of the atomic, electronic and magnetic properties of novel materials for energy applications, first of all, fuel cells transforming chemical energy into electricity, and batteries. Understanding of a role of dopants and defects is a key for prediction of improvement of device performance which is validated later on experimentally. Addressing realistic operational conditions is achieved via combination with ab initio thermodynamics. The state of the art first principles calculations of large and low symmetry are very time consuming and need use of supercomputer technologies as provided at HLRS in Stuttgart.

In the search of new physics, some proposed models fall into the category of nearly conformal Strongly Coupled Gauge Theories (SCGTs). Such theories are identified by the almost existence of non-trivial zero (pseudo infrared fixed point) in their beta functions. In this project, the Lattice Higgs Collaboration quantitatively investigates the beta function of nearly conformal SCGTs and observes how the beta function depends on the number of fermion flavors and representations. This provides insight of how SCGTs approach near conformality, which is crucial in the identification of models suitable for the development of new physics.

With constantly growing fuel prices and toughening of environmental legislation, the vehicle industry is struggling to reduce fuel consumption and decrease emission levels for the new and existing vehicles. One way to achieve this goal is to improve aerodynamic performance by decreasing aerodynamic resistance. Leveraging HPC resources, researchers of the Technical University of Munich conducted a wide range of studies with the aim to improve modeling techniques, develop a profound understanding for flow phenomena, and optimize vehicle shapes.

This group from the Helmholtz-Institute Erlangen-Nürnberg performed simulations, both on a coarse-grained and a molecular level of detail, elucidating how so-called antagonistic salts, consisting of a large anion and a small cation, trigger the spontaneous formation of highly regular, nanometer sized structures in water/oil mixtures. Due to their size difference the small cations accumulate in the water phase while the large anions go to the oil phase. The resulting electrostatic interactions between the phases can lead to long-range ordering.

Ziegler-Natta catalysts are important for industry, but determining exactly how they work is difficult due to their complex nature which involves a number of different active compounds on nano-sized structures. Researchers of the University of Turin led by Dr. Maddalena D’Amore have been using Density Functional Theory (DFT) to try to find out more about these types of systems.

Recently, European legislative bodies have imposed significant restrictions on the emission level of Diesel injections systems, thus challenging car manufacturers and suppliers to reduce pollution. Improvement of the combustion process and spray quality has become a main objective in fulfilling those policies, which is mainly achieved by increasing injection pressures. Therefore, understanding internal nozzle flows has become a key aspect in designing efficient and durable Diesel injection systems. Since to this date quantitative experimental investigations are challenging, researchers use HPC technologies and computational fluid dynamics to complement experimental findings by providing additional information about the flow topology.

Understanding the most inner structure of matter has been a driving force of science sine the idea of an "atom" by the antique Greeks. And, with todays supercomputer power we are now in the fascinating position to finally reveal what holds the world together. As a most important step in this direction, in this project, basic properties of the proton, e.g. the spin, the angular momentum and the quark and gluon content as well as their distribution within the proton have been calculated. This constitutes a pioneering step to understand the nature of matter, the very early universe and ultimatley to answer the question where we are coming from.

A series of highly resolved direct numerical simulations (DNSs) of temporally evolving turbulent non-premixed jet flames was conducted on the SuperMUC of LRZ. Two promising approaches were used to analyze the databases. The first approach, on-the-fly tracking flamelet structure, helps to understand the effects of neglecting tangential diffusion (TD) on the performance of classical flamelet models. The second approach - dissipation elements – helps to develop possible closure strategies for including flame-tangential effects in the flamelet models. Moreover, TD was used as an important performance indicator to assess tabulation strategies, differential diffusion effects, and Soret effects in turbulent non-premixed combustion.

A high fidelity, high Reynolds number direct numerical simulation (DNS) of a planar temporally evolving non-premixed jet flame was performed on the new supercomputer JUWELS of JSC. The DNS enabled the detailed investigation of combustion conditions with a high level of scale interaction between combustion chemistry and turbulence. Furthermore, the simulation was instrumental in understanding how the structure of scalar fields is affected by heat release in non-premixed flames. The insights gained from the DNS are instrumental in the development of new combustion models with the goal of improving the accuracy of simulations of real-world engineering applications.

An effective cooling of the gas turbine components subject to high thermal stresses is vital for the success of new engine and combustion concepts aiming at achieving further improvements in the energy conversion efficiency of the overall machine. The use of pulsating impinging jets - which enlarge vortex structures naturally occurring in the impinging jet flow when no pulsation is enforced - is a promising approach to develop a substantially more performant cooling system. To gain a deeper understanding of how the vortex system behaves under realistic conditions, researchers performed a DNS of a non-pulsating impinging jet flow with fully turbulent inflow conditions and compared its results with a reference case with a laminar inflow.

Intramembrane proteases control the activity of membrane proteins and occur in all organisms. A prime example is g-secretase, cleaving the amyloid precursor protein, whose misprocessing is related to onset and progression of Alzheimer's disease. Since a protease's biological function depends on its substrate spectrum, it is essential to study the repertoire of natural substrates as well as determinants and mechanisms of substrate recognition and cleavage—which is the aim of this collaborative research project. Conformational flexibility of substrate and enzyme plays an essential role for recognition, complex formation and subsequent relaxation steps leading to cleavage and product release.

The development of novel sustainable biocatalytic processes requires systematic studies of the molecular interactions between enzymes, substrates, and solvents. Based on the HLRS HPC infrastructure, comprehensive molecular simulations were performed to investigate substrate binding in enzymatic reaction systems.

Studying the mechanochemistry of disulfide systems upon nucleophilic attack is a very rich field where each system requires computing resources and CPU time that can only be provided by very powerful supercomputers such as provided by the Gauss Centre for Supercomputing. Simulations run on JUQUEEN of JSC in the course of this project offered a wealth of surprises and novel insights into mechanochemical reactions. While they resulted in discovering unexpected reaction mechanisms, they - amongst others - brought to light an unknown phenomenon with respect to splitting disulphide bonds in water.

Diabetes reaches epidemic proportions with a major and growing economic impact on the society. An effective treatment requires atomic-level understanding of how insulin acts on cells. Using molecular dynamics simulations, an international team of researchers studied the process of insulin binding to its receptor and the resulting structural changes at atomic scale with cryogenic election microscopy and atomistic MD simulation. The results of these studies were recently published in the Journal of Cell Biology.

Convection in the Earth’s mantle is the driving force behind large scale geologic activity such as plate tectonics and continental drift. As such it is related to phenomena like e.g. earthquakes, mountain building, and hot-spot volcanism. Laboratory experiments naturally fail to reproduce the pressures and temperatures in the mantle, thus simulation is a key ingredient in the research of mantle convection. However, since simulating convection in the Earth’s mantle is a very resource consuming HPC application as it requires extremely large grids and many time steps in order to allow models with realistic geological parameters, researchers turn towards GCS supercomputers to tackle this challenge.

Researchers of Forschungszentrum Jülich used the computing resources of high-performance computing system JUQUEEN of JSC to improve the understanding of the QCD transition.

Researchers of the Institute of Aerodynamics at RWTH Aachen University used large-eddy simulation and computational aeroacoustics methods to analyze noise sources in turbulent flames and the interaction of the resulting acoustic waves with the flame and the turbulent flow field. To achieve accurate results of the flow and the acoustic field highly resolved large-scale simulations with several hundred million mesh points are necessary. The simulation results have given new insights into fundamental sound-generation mechanisms and their phase-relationship that are important for the prediction and control of thermoacoustic instabilities, and ultimately, the development of more efficient and gas turbines with lower pollutant emissions.

It is a long lasting dream in nuclear physics to study nuclei like, for instance, carbon directly from Quantum Chromodynamics (QCD), the underlying fundamental theory of strong interactions. Such an endeavor is very challenging both, methodically and numerically. Towards this goal physicists from the European Twisted Mass Collaboration and in particular the University of Bonn have started to investigate two hadron systems using the approach of Lattice QCD.

Researchers carried out density functional theory defect calculations of materials relevant in energy applications. They calculated Raman spectra of LiCoO2 which allow to follow the structural evolution during charging and discharging of this important class of lithium-ion battery cathode materials and to understand what can lead to their failure. Furthermore, the effect of defects forming on a dissolving metastable surface on the (photo)electrocatalytic performance were calculated, and the team worked on novel computational methods applied to defects that will enable DFT calculations of defects with a similar accuracy than state-of-the-art methods, however at a much-reduced computational cost.

ATP synthase is an enzyme found in organisms ranging from primitive bacteria to some of the most complex lifeforms, such as humans. Its energetic efficiency is unrivalled, but not well understood. Researchers of Gdansk University of Technology have been using HPC to study this remarkable enzyme at a level of detail never seen before.

In recent years, hydroelectric power plants have received increased attention for the role they play in integrating volatile renewable energies that contribute to stabilizing the electrical grid. One major issue, though, is rooted in running turbines under conditions they were not originally designed for, leading to undesirable flow phenomena. With the standard modeling approaches that are typically used in industry simulations of hydroelectric turbines, simulation accuracy in scenarios where the turbine is used off-design is rather poor. The goal of this project is to increase simulation accuracy by the selection of suitable modeling approaches and the use of a fine mesh resolution, which is only possible by the use of supercomputers.

Rotating convection is ubiquitous in geophysical systems. In generates the Earth magnetic field, stirs the deep atmospheres of giant planets and possibly also drives their strong surface winds. A thorough understanding of these objects requires comprehensive insight into the physics of turbulent convective flows that are strongly constrained by Coriolis forces. Numerical simulations reveal the full three-dimensional structure of the flow, and can be used to guide theoretical modeling.

A GCS large-scale project under leadership of Dr.-Ing. Cifuentes of the University of Duisburg-Essen aims at understanding the physics of entrainment in turbulent premixed flames. This research characterizes the entrainment processes through the study and comparison of the flame front and the enstrophy interface. This is an essential issue in reactive turbulent flows, because a better understanding of the dynamics of the flame front and the enstrophy interface leads to better predictions of flame instabilities and scalar structures.

In this project, researchers use state of the art fermion quantum Monte Carlo methods to understand emergent collective phenomena in correlated electron system. The scientists define and study theoretical models where topology emerges and leads to novel particles at quantum critical points. The flexibility of their approach also makes it possible to study the physics of magnetic moments in a metallic environment. This could, for instance, enable theoretical experiments for understanding magnetic adatoms on metallic surfaces. In this report, a succinct account of the ALF (Algorithms, Lattice, Fermions) program package, which was developed by the scientists, as well as a summary of selected research projects is provided.

In eukaryotes, conversion of foodstuff into electrochemical energy takes place in mitochondria by enzymes of the respiratory chain. Cytochrome c oxidase (CcO) reduces oxygen to water and pumps protons across the membrane. In this project, we elucidated how reduction of metal co-factors in CcO control the proton transfer dynamics. By combining atomistic MD simulations with hybrid QM/MM free energy calculations, we elucidated the location of a transient proton loading site near the active site, and identified how proton channels are activated during the different steps of the catalytic cycle.

Computational mechanistic modelling using systems of ordinary differential equations (ODE) has become an integral tool in systems biology. Parameters of such models are often not known in advance and need to be inferred from experimental data, which is computationally very expensive. The SuperMUC supercomputer enabled researchers from the Helmholtz Zentrum Munich to evaluate state-of-the-art algorithms and to develop novel, more efficient algorithms for parameter estimation from large datasets and relative measurements.

Shock-related buffeting is a phenomenon that occurs when air passes over the wing of an aeroplane under extreme conditions and can have profound consequences for how wings are engineered and their durability. Leveraging the computing capacities of HPC system Hazel Hen, researchers at the University of Southampton have been investigating this phenomenon using direct numerical simulations.

Spray painting is the most common application technique in coating technology. Typical atomizers used in spray coating industries are such as High-speed rotary bell and spray guns with compressed air. High-speed rotary bell atomizers provide an excellent paint film quality as well as high transfer efficiencies (approx. 90%) due to electrostatic support. Small and medium-sized enterprises continue, however, to use compressed air atomizers, although they no longer meet today's requirements from an economic and environmental point of view. It is very important to understand the atomization mechanisms of these two kinds of atomizers, in order to improve the paint quality, to reduce the overspray and to optimize the coating process.

Many wall-bounded flows in nature and technology are affected by the surface roughness of the wall. In some cases, this has adverse effects, e.g. drag increase leading to higher fuel costs; in others, it is beneficial for mixing enhancement or transfer properties. Computationally, it is notoriously difficult to simulate these flows because of the vast separation of scales in highly turbulent flows and the challenges involved in handling complex geometries. The studies are carried out in two paradigmatic and complementary systems in turbulence research, Taylor-Couette and Rayleigh-Bénard flow.

Nucleons make up more than 99% of the mass of ordinary matter. Computing their properties from first principles, i.e. the theory of Quantum Chromodynamics, is complicated by the non-linear nature of the underlying equations. Only by using supercomputers can we attempt to compute these quantities with the necessary precision. Beyond shedding light on the nature of the nucleons, the results help to resolve some long-standing puzzles in nucleon structure physics and restrict possible models of physics beyond the Standard Model.

Direct bandgap silicon can be the key to integrate both electronic and optical functionalities on a silicon platform. Despite considerable effort, achieving light emission from group IV semiconductors has remained unattainable until now. Very recently, ab initio calculations combined with experiments could prove that Ge-rich hexagonal crystal phases of SixGe1-x feature a direct bandgap, tunable in a frequency range coinciding with the low loss window for optical fiber communications. Efficient light emission from direct band gap SiGe could also be shown. Further calculations explore how to engineer light emission by strain and alloying.

Molecular modelling and simulation is an established method for describing and predicting thermodynamic properties of fluids. This project examines interfacial properties of fluids, their contact with solid materials, interfacial fluctuations and finite-size effects, linear transport coefficients in the bulk and at interfaces and surfaces as well as transport processes near and far from equilibrium. These phenomena are investigated by massively-parallel molecular dynamics simulation based on quantitatively reliable classical-mechanical force fields.

Flow noise during takeoff and landing of commercial aircraft can be substantially reduced by the use of porous surface layers in suitable sections of the airfoil. However, porosity and roughness of surfaces tend to have an adverse effect on the boundary layer and thus on the lift of wings. This motivates the need to be able to predict the aerodynamic effects of porous segments of the wing surface by numerical methods. Due to the inherent requirements of resolving both the turbulence on the scale of an airfoil and the flow inside the pore-scale resolved porous medium, the simulations run on SuperMUC required more than a billion grid nodes on a locally refined three-dimensional mesh.

The interaction between fluids and structures (fluid structure interaction/FSI) is a topic of interest in many science fields. In addition to experimental investigations, numerical simulations have become a valuable tool to foresee complex flow phenomena such as vortex shedding, transition and separation or critical stresses in the structure exposed to the flow. In civil engineering, e.g., structures are exposed to strong variations of the wind, particularly wind gusts, and such high loads can ultimately lead to a complete destruction of the structure. Scientists are leveraging HPC technologies in order to model wind gusts and to comprehend their impact on the FSI phenomenon.

Turbulent thermal convection plays an essential role in a wide range of natural and industrial settings, from astrophysical and geophysical flows to process engineering. While heat transfer in industrial applications takes place in confined systems, the aspect ratio in many natural instances of convection is huge. Interestingly, flow organization on enormous scales is observed in, for example, oceanic and atmospheric convection. However, our physical understanding of the formation of turbulent superstructures is limited. In this project, we analyze the flow organization within turbulent superstructures and show that their size increases when the thermal driving is increased.

Nuclear matter changes at high temperatures from a gas of hadrons into a quark-gluon plasma. For sufficiently high temperatures this quark-gluon plasma can be described in terms of effective field theory calculations assuming weak coupling. In this project, scientists calculate the QCD Equation of State and the free energies of heavy quark systems using Lattice QCD, a Markov Chain Monte Carlo approach for solving the QCD path integral numerically in an imaginary time formalism. By comparing the continuum extrapolated results to weak-coupling calculations in different EFT frameworks, their applicability is being established.

The efficient mixing of fuel and oxidizer is essential in modern combustion engines. Especially in supersonic combustion the rapid mixing of fuel and oxidizer is of crucial importance as the detention time of the fuel-oxidizer mixture in the combustion chamber is only a few milliseconds. The shock-induced Richtmyer-Meshkov instability (RMI) promotes mixing and thus has the potential to increase the burning efficiency of supersonic combustion engines. To study the interaction between RMI and shock-induced reaction waves, which affects the flow field evolution und the mixing significantly, researchers leveraged HPC system SuperMUC to run 3D simulations of reacting shock-bubble interaction.

Turbulent convection flows in nature and technology often show prominent and nearly regular patterns on their largest scales which we term turbulent superstructures. Their appearance challenges the classical picture of turbulence in which a turbulent flow is considered as a tangle of chaotically moving vortices. Examples for superstructures in nature are cloud streets in the atmosphere or the granulation at the surface of the Sun. In several applications, this structure formation is additionally affected by magnetic fields. Our understanding of the origin of turbulent superstructures and their role for the turbulent transport is presently still incomplete and will be improved by direct numerical simulations of turbulent convection.

In this study the axion particles are investigated numerically. To guide experimental searches of the axion particle, its mass needs to be estimated theoretically. For this one needs to study the creation mechanisms of the axions in the early universe. The axion fields can form topological defects known as cosmological strings, which are highly energetic string-like excitations which decay into axion particles. The axions are created in a phase transition where a large part of the energy builds strings. In this study we follow the fate of the axion-string network to understand how the axion abundance in the universe is created.

Researchers of the Institute of Aerodynamics (AIA) at RWTH Aachen University conducted large-scale benchmark simulations on supercomputer Hazel Hen of the High-Performance Computing Center Stuttgart to analyze the interaction of non-spherical particles with turbulent flows. These simulations provide a unique data base for the development of simple models which can be applied to study complex engineering problems. Such models are required in a larger research framework to improve the efficiency of pulverized coal and biomass combustion to significantly reduce the CO2 emissions.

Designing new enzymes is a grand challenge for modern biochemistry, and there are few examples for artificial enzymes with significant catalytic rate accelerations. We have developed a new method for computational enzyme design where we mimic evolution in nature and randomly mutate amino acids using a Metropolis Monte Carlo (MC) procedure. The aim of the method is to identify substitutions that increase the catalytic activity of enzymes. We probe the catalytic activity by quantum mechanics/classical mechanics (QM/MM) calculations, which are important for accurately modeling chemical reactions.

Transient mixing and ignition play a significant role in many systems, where combustion efficiency and emissions are controlled by ignition and mixing dynamics. In the present work, high fidelity simulations of a pulsed fuel injection system are carried out using state of the art numerical tools and high-performance computing. The results contain all parameters that affect ignition dynamics and are mined and analyzed. The physics of transient reactive turbulent jets are thus identified and presented that partners in industry and academia can improve their understanding of the process and work on the design of better combustion devices.

Using the high-performance computing resources available at the Jülich Supercomputing Centre, scientists computed the mass difference between the up and down quarks. The result has been published in Physical Review Letters.

Regional climate simulations at the convection-permitting scale (< 4 km) have the potential to improve seasonal forecasts, especially where complex topography hinders global models. Due to high computational costs, tests using state-of-the-art ensemble forecasts have not been performed yet. In this one-year case study, a Weather Research and Forecasting (WRF) multi-physics ensemble was used to downscale the SEAS5 ensemble forecast over the Horn of Africa. Reliability of precipitation prediction is improved, although the global model’s biases in temperature and precipitation are not reduced. Measurable added value against the global model is provided for intense precipitation statistics over the Ethiopian highlands.

Granular matter is typically the result of random pattern formation in a solid, like breaking a glass or pulverizing a rock into pieces of variable sizes. Faraday waves are patterns that appear on a fluid that is perturbed by an external drive that oscillates in resonance. Faraday waves aren't random; in contrast to granular matter, these waves are regular, standing, periodic patterns, seen for instance in liquids in a vessel that is shaken. Surprisingly, granulation and Faraday waves can exist in quantum systems too and, even more surprisingly, they can be produced in the same quantum system: in a gas of trapped atoms cooled very close to absolute zero temperature. When the strength of interactions between atoms is modulated, a Faraday…

The electric dipole moment of the neutron, measuring the distance of positive and negative charge density in the neutron as shown in the image (left), provides a unique and sensitive probe to physics beyond the Standard Model. It has played an important part over many decades in shaping and constraining numerous models of CP violation. QCD allows for CP-violating effects that propagate into the hadronic sector via the so-called θ term Sθ in the action, S = S + Sθ, with Sθ = i θ Q, where Q is the topological charge. In this project the electric dipole mom