TY - CHAP A1 - Noack, Matthias A1 - Reinefeld, Alexander A1 - Kramer, Tobias A1 - Steinke, Thomas T1 - DM-HEOM: A Portable and Scalable Solver-Framework for the Hierarchical Equations of Motion T2 - 2018 IEEE International Parallel and Distributed Processing Symposium Workshops (IPDPSW), 19th IEEE Int. Workshop on Parallel and Distributed Scientific and Engineering Computing (PDSEC 2018) N2 - Computing the Hierarchical Equations of Motion (HEOM) is by itself a challenging problem, and so is writing portable production code that runs efficiently on a variety of architectures while scaling from PCs to supercomputers. We combined both challenges to push the boundaries of simulating quantum systems, and to evaluate and improve methodologies for scientific software engineering. Our contributions are threefold: We present the first distributed memory implementation of the HEOM method (DM-HEOM), we describe an interdisciplinary development workflow, and we provide guidelines and experiences for designing distributed, performance-portable HPC applications with MPI-3, OpenCL and other state-of-the-art programming models. We evaluated the resulting code on multi- and many-core CPUs as well as GPUs, and demonstrate scalability on a Cray XC40 supercomputer for the PS I molecular light harvesting complex. Y1 - 2018 SN - 978-1-5386-5555-9 U6 - https://doi.org/10.1109/IPDPSW.2018.00149 N1 - Best Paper Award SP - 947 EP - 956 ER - TY - JOUR A1 - Kramer, Tobias A1 - Noack, Matthias A1 - Reinefeld, Alexander A1 - Rodríguez, Mirta A1 - Zelinskyi, Yaroslav T1 - Efficient calculation of open quantum system dynamics and time-resolved spectroscopy with Distributed Memory HEOM (DM-HEOM) JF - Journal of Computational Chemistry N2 - Time- and frequency resolved optical signals provide insights into the properties of light harvesting molecular complexes, including excitation energies, dipole strengths and orientations, as well as in the exciton energy flow through the complex. The hierarchical equations of motion (HEOM) provide a unifying theory, which allows one to study the combined effects of system-environment dissipation and non-Markovian memory without making restrictive assumptions about weak or strong couplings or separability of vibrational and electronic degrees of freedom. With increasing system size the exact solution of the open quantum system dynamics requires memory and compute resources beyond a single compute node. To overcome this barrier, we developed a scalable variant of HEOM. Our distributed memory HEOM, DM-HEOM, is a universal tool for open quantum system dynamics. It is used to accurately compute all experimentally accessible time- and frequency resolved processes in light harvesting molecular complexes with arbitrary system-environment couplings for a wide range of temperatures and complex sizes. Y1 - 2018 U6 - https://doi.org/doi:10.1002/jcc.25354 VL - 39 IS - 22 SP - 1779 EP - 1794 PB - Wiley Periodicals, Inc. ER - TY - JOUR A1 - Kramer, Tobias A1 - Noack, Matthias A1 - Reimers, Jeffrey R. A1 - Reinefeld, Alexander A1 - Rodríguez, Mirta A1 - Yin, Shiwei T1 - Energy flow in the Photosystem I supercomplex: comparison of approximative theories with DM-HEOM JF - Chemical Physics N2 - We analyze the exciton dynamics in PhotosystemI from Thermosynechococcus elongatus using the distributed memory implementation of the hierarchical equation of motion (DM-HEOM) for the 96 Chlorophylls in the monomeric unit. The exciton-system parameters are taken from a first principles calculation. A comparison of the exact results with Foerster rates and Markovian approximations allows one to validate the exciton transfer times within the complex and to identify deviations from approximative theories. We show the optical absorption, linear, and circular dichroism spectra obtained with DM-HEOM and compare them to experimental results. Y1 - 2018 U6 - https://doi.org/10.1016/j.chemphys.2018.05.028 VL - 515 SP - 262 EP - 271 PB - Elsevier B.V. ER - TY - JOUR A1 - Kramer, Tobias A1 - Noack, Matthias A1 - Baum, Daniel A1 - Hege, Hans-Christian A1 - Heller, Eric J. T1 - Dust and gas emission from cometary nuclei: the case of comet 67P/Churyumov-Gerasimenko JF - Advances in Physics: X N2 - Comets display with decreasing solar distance an increased emission of gas and dust particles, leading to the formation of the coma and tail. Spacecraft missions provide insight in the temporal and spatial variations of the dust and gas sources located on the cometary nucleus. For the case of comet 67P/Churyumov-Gerasimenko (67P/C-G), the long-term obser- vations from the Rosetta mission point to a homogeneous dust emission across the entire illuminated surface. Despite the homogeneous initial dis- tribution, a collimation in jet-like structures becomes visible. We propose that this observation is linked directly to the complex shape of the nucleus and projects concave topographical features into the dust coma. To test this hypothesis, we put forward a gas-dust description of 67P/C-G, where gravitational and gas forces are accurately determined from the surface mesh and the rotation of the nucleus is fully incorporated. The emerging jet-like structures persist for a wide range of gas-dust interactions and show a dust velocity dependent bending. Y1 - 2018 U6 - https://doi.org/10.1080/23746149.2017.1404436 VL - 3 IS - 1 SP - 1404436 ER - TY - CHAP A1 - Läuter, Matthias A1 - Kramer, Tobias A1 - Rubin, Martin A1 - Altwegg, Kathrin T1 - Gas production of comet 67P/Churyumov-Gerasimenko reconstructed from DFMS/COPS data T2 - Europlanet Science Congress N2 - We reconstruct the temporal evolution of surface emissions for the four major gas species H2O, CO2, CO, and O2 emitted during the 2015 apparition of comet 67P/Churyumov-Gerasimenko (67P/C-G). Measurements from the Double Focusing Mass Spectrometer (DFMS) of the Rosetta Orbiter Spectrometer for Ion and Neutral Analysis (ROSINA) and the COmet Pressure Sensor (COPS) are used to determine the gas sources on the surface with an inverse gas model for the entire coma. For all species, peak production rates and integrated production rates per orbit are evaluated separately for the northern and the southern hemisphere. Complemented with the total mass production, this allows us to estimate the dust-to-gas ratio of the emitted material. Y1 - 2018 VL - 12 SP - EPSC2018-515-1 PB - Europlanet Science Congress ER - TY - CHAP A1 - Kramer, Tobias A1 - Läuter, Matthias T1 - Modelling the inner coma of comet 67P/Churyumov-Gerasimenko T2 - Europlanet Science Congress N2 - Based on about 1 million of pressure measurements around comet 67P/Churyumov-Gerasimenko we reconstruct the gas emission across the entire nucleus. Dust particles are seeded in the gas model and the resulting dust distribution follows a daily pattern which agrees with observations if a uniform dust release across the entire sunlit surface is assumed. Y1 - 2018 VL - 12 SP - EPSC2018-115-1 PB - Europlanet Science Congress ER -