@article{KramerNoackBaumetal., author = {Kramer, Tobias and Noack, Matthias and Baum, Daniel and Hege, Hans-Christian and Heller, Eric J.}, title = {Homogeneous dust emission and jet structure near active cometary nuclei: the case of 67P/Churyumov-Gerasimenko}, abstract = {We compute trajectories of dust grains starting from a homogeneous surface activity-profile on a irregularly shaped cometary nucleus. Despite the initially homogeneous dust distribution a collimation in jet-like structures becomes visible. The fine structure is caused by concave topographical features with similar bundles of normal vectors. The model incorporates accurately determined gravitational forces, rotation of the nucleus, and gas-dust interaction. Jet-like dust structures are obtained for a wide range of gas-dust interactions. For the comet 67P/Churyumov-Gerasimenko, we derive the global dust distribution around the nucleus and find several areas of agreement between the homogeneous dust emission model and the Rosetta observation of dust jets, including velocity-dependent bending of trajectories.}, language = {en} } @article{KramerNoack, author = {Kramer, Tobias and Noack, Matthias}, title = {Prevailing Dust-transport Directions on Comet 67P/Churyumov-Gerasimenko}, series = {The Astrophysical Journal Letters}, volume = {813}, journal = {The Astrophysical Journal Letters}, number = {2}, doi = {10.1088/2041-8205/813/2/L33}, pages = {L33}, abstract = {Dust transport and deposition behind larger boulders on the comet 67P/Churyumov-Gerasimenko (67P/C-G) have been observed by the Rosetta mission. We present a mechanism for dust-transport vectors based on a homogeneous surface activity model incorporating in detail the topography of 67P/C-G. The combination of gravitation, gas drag, and Coriolis force leads to specific dust transfer pathways, which for higher dust velocities fuel the near-nucleus coma. By distributing dust sources homogeneously across the whole cometary surface, we derive a global dust-transport map of 67P/C-G. The transport vectors are in agreement with the reported wind-tail directions in the Philae descent area.}, language = {en} } @misc{KramerNoack, author = {Kramer, Tobias and Noack, Matthias}, title = {On the origin of inner coma structures observed by Rosetta during a diurnal rotation of comet 67P/Churyumov-Gerasimenko.}, issn = {1438-0064}, doi = {10.3847/2041-8205/823/1/L11}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-59225}, abstract = {The Rosetta probe around comet 67P/Churyumov-Gerasimenko (67P) reveals an anisotropic dust distribution of the inner coma with jet-like structures. The physical processes leading to jet formation are under debate, with most models for cometary activity focusing on localised emission sources, such as cliffs or terraced regions. Here we suggest, by correlating high-resolution simulations of the dust environment around 67P with observations, that the anisotropy and the background dust density of 67P originate from dust released across the entire sunlit surface of the nucleus rather than from few isolated sources. We trace back trajectories from coma regions with high local dust density in space to the non-spherical nucleus and identify two mechanisms of jet formation: areas with local concavity in either two dimensions or only one. Pits and craters are examples of the first case, the neck region of the bilobed nucleus of 67P for the latter one. The conjunction of multiple sources in addition to dust released from all other sunlit areas results in a high correlation coefficient (∼0.8) of the predictions with observations during a complete diurnal rotation period of 67P.}, language = {en} } @article{KramerNoack, author = {Kramer, Tobias and Noack, Matthias}, title = {On the origin of inner coma structures observed by Rosetta during a diurnal rotation of comet 67P/Churyumov-Gerasimenko}, series = {The Astrophysical Journal Letters}, volume = {823}, journal = {The Astrophysical Journal Letters}, number = {1}, doi = {10.3847/2041-8205/823/1/L11}, pages = {L11}, abstract = {The Rosetta probe around comet 67P/Churyumov-Gerasimenko (67P) reveals an anisotropic dust distribution of the inner coma with jet-like structures. The physical processes leading to jet formation are under debate, with most models for cometary activity focusing on localized emission sources, such as cliffs or terraced regions. Here we suggest, by correlating high-resolution simulations of the dust environment around 67P with observations, that the anisotropy and the background dust density of 67P originate from dust released across the entire sunlit surface of the nucleus rather than from few isolated sources. We trace back trajectories from coma regions with high local dust density in space to the non-spherical nucleus and identify two mechanisms of jet formation: areas with local concavity in either two dimensions or only one. Pits and craters are examples of the first case; the neck region of the bi-lobed nucleus of 67P is an example of the latter case. The conjunction of multiple sources, in addition to dust released from all other sunlit areas, results in a high correlation coefficient (~0.8) of the predictions with observations during a complete diurnal rotation period of 67P.}, language = {en} } @misc{KramerNoackBaumetal.2017, author = {Kramer, Tobias and Noack, Matthias and Baum, Daniel and Hege, Hans-Christian and Heller, Eric J.}, title = {Dust and gas emission from cometary nuclei: the case of comet 67P/Churyumov-Gerasimenko}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-66338}, year = {2017}, abstract = {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.}, language = {en} } @article{KramerNoackBaumetal., author = {Kramer, Tobias and Noack, Matthias and Baum, Daniel and Hege, Hans-Christian and Heller, Eric J.}, title = {Dust and gas emission from cometary nuclei: the case of comet 67P/Churyumov-Gerasimenko}, series = {Advances in Physics: X}, volume = {3}, journal = {Advances in Physics: X}, number = {1}, doi = {10.1080/23746149.2017.1404436}, pages = {1404436}, abstract = {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.}, language = {en} } @inproceedings{NoackReinefeldKrameretal., author = {Noack, Matthias and Reinefeld, Alexander and Kramer, Tobias and Steinke, Thomas}, title = {DM-HEOM: A Portable and Scalable Solver-Framework for the Hierarchical Equations of Motion}, series = {2018 IEEE International Parallel and Distributed Processing Symposium Workshops (IPDPSW), 19th IEEE Int. Workshop on Parallel and Distributed Scientific and Engineering Computing (PDSEC 2018)}, booktitle = {2018 IEEE International Parallel and Distributed Processing Symposium Workshops (IPDPSW), 19th IEEE Int. Workshop on Parallel and Distributed Scientific and Engineering Computing (PDSEC 2018)}, isbn = {978-1-5386-5555-9}, doi = {10.1109/IPDPSW.2018.00149}, pages = {947 -- 956}, abstract = {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.}, language = {en} } @article{KramerNoackReinefeldetal., author = {Kramer, Tobias and Noack, Matthias and Reinefeld, Alexander and Rodr{\´i}guez, Mirta and Zelinskyi, Yaroslav}, title = {Efficient calculation of open quantum system dynamics and time-resolved spectroscopy with Distributed Memory HEOM (DM-HEOM)}, series = {Journal of Computational Chemistry}, volume = {39}, journal = {Journal of Computational Chemistry}, number = {22}, publisher = {Wiley Periodicals, Inc.}, doi = {doi:10.1002/jcc.25354}, pages = {1779 -- 1794}, abstract = {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.}, language = {en} } @article{KramerNoackReimersetal., author = {Kramer, Tobias and Noack, Matthias and Reimers, Jeffrey R. and Reinefeld, Alexander and Rodr{\´i}guez, Mirta and Yin, Shiwei}, title = {Energy flow in the Photosystem I supercomplex: comparison of approximative theories with DM-HEOM}, series = {Chemical Physics}, volume = {515}, journal = {Chemical Physics}, publisher = {Elsevier B.V.}, doi = {10.1016/j.chemphys.2018.05.028}, pages = {262 -- 271}, abstract = {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.}, language = {en} }