TY - JOUR A1 - Schneck, Jakob A1 - Weiser, Martin A1 - Wende, Florian T1 - Impact of mixed precision and storage layout on additive Schwarz smoothers JF - Numerical Linear Algebra with Applications N2 - The growing discrepancy between CPU computing power and memory bandwidth drives more and more numerical algorithms into a bandwidth-bound regime. One example is the overlapping Schwarz smoother, a highly effective building block for iterative multigrid solution of elliptic equations with higher order finite elements. Two options of reducing the required memory bandwidth are sparsity exploiting storage layouts and representing matrix entries with reduced precision in floating point or fixed point format. We investigate the impact of several options on storage demand and contraction rate, both analytically in the context of subspace correction methods and numerically at an example of solid mechanics. Both perspectives agree on the favourite scheme: fixed point representation of Cholesky factors in nested dissection storage. Y1 - 2021 U6 - https://doi.org/10.1002/nla.2366 VL - 28 IS - 4 ER - TY - JOUR A1 - Schimunek, Johannes A1 - Seidl, Philipp A1 - Elez, Katarina A1 - Hempel, Tim A1 - Le, Tuan A1 - Noé, Frank A1 - Olsson, Simon A1 - Raich, Lluís A1 - Winter, Robin A1 - Gokcan, Hatice A1 - Gusev, Filipp A1 - Gutkin, Evgeny M. A1 - Isayev, Olexandr A1 - Kurnikova, Maria G. A1 - Narangoda, Chamali H. A1 - Zubatyuk, Roman A1 - Bosko, Ivan P. A1 - Furs, Konstantin V. A1 - Karpenko, Anna D. A1 - Kornoushenko, Yury V. A1 - Shuldau, Mikita A1 - Yushkevich, Artsemi A1 - Benabderrahmane, Mohammed B. A1 - Bousquet-Melou, Patrick A1 - Bureau, Ronan A1 - Charton, Beatrice A1 - Cirou, Bertrand C. A1 - Gil, Gérard A1 - Allen, William J. A1 - Sirimulla, Suman A1 - Watowich, Stanley A1 - Antonopoulos, Nick A1 - Epitropakis, Nikolaos A1 - Krasoulis, Agamemnon A1 - Itsikalis, Vassilis A1 - Theodorakis, Stavros A1 - Kozlovskii, Igor A1 - Maliutin, Anton A1 - Medvedev, Alexander A1 - Popov, Petr A1 - Zaretckii, Mark A1 - Eghbal-Zadeh, Hamid A1 - Halmich, Christina A1 - Hochreiter, Sepp A1 - Mayr, Andreas A1 - Ruch, Peter A1 - Widrich, Michael A1 - Berenger, Francois A1 - Kumar, Ashutosh A1 - Yamanishi, Yoshihiro A1 - Zhang, Kam Y. J. A1 - Bengio, Emmanuel A1 - Bengio, Yoshua A1 - Jain, Moksh J. A1 - Korablyov, Maksym A1 - Liu, Cheng-Hao A1 - Marcou, Gilles A1 - Glaab, Enrico A1 - Barnsley, Kelly A1 - Iyengar, Suhasini M. A1 - Ondrechen, Mary Jo A1 - Haupt, V. Joachim A1 - Kaiser, Florian A1 - Schroeder, Michael A1 - Pugliese, Luisa A1 - Albani, Simone A1 - Athanasiou, Christina A1 - Beccari, Andrea A1 - Carloni, Paolo A1 - D’Arrigo, Giulia A1 - Gianquinto, Eleonora A1 - Goßen, Jonas A1 - Hanke, Anton A1 - Joseph, Benjamin P. A1 - Kokh, Daria B. A1 - Kovachka, Sandra A1 - Manelfi, Candida A1 - Mukherjee, Goutam A1 - Muñiz-Chicharro, Abraham A1 - Musiani, Francesco A1 - Nunes-Alves, Ariane A1 - Paiardi, Giulia A1 - Rossetti, Giulia A1 - Sadiq, S. Kashif A1 - Spyrakis, Francesca A1 - Talarico, Carmine A1 - Tsengenes, Alexandros A1 - Wade, Rebecca C. A1 - Copeland, Conner A1 - Gaiser, Jeremiah A1 - Olson, Daniel R. A1 - Roy, Amitava A1 - Venkatraman, Vishwesh A1 - Wheeler, Travis J. A1 - Arthanari, Haribabu A1 - Blaschitz, Klara A1 - Cespugli, Marco A1 - Durmaz, Vedat A1 - Fackeldey, Konstantin A1 - Fischer, Patrick D. A1 - Gorgulla, Christoph A1 - Gruber, Christian A1 - Gruber, Karl A1 - Hetmann, Michael A1 - Kinney, Jamie E. A1 - Padmanabha Das, Krishna M. A1 - Pandita, Shreya A1 - Singh, Amit A1 - Steinkellner, Georg A1 - Tesseyre, Guilhem A1 - Wagner, Gerhard A1 - Wang, Zi-Fu A1 - Yust, Ryan J. A1 - Druzhilovskiy, Dmitry S. A1 - Filimonov, Dmitry A. A1 - Pogodin, Pavel V. A1 - Poroikov, Vladimir A1 - Rudik, Anastassia V. A1 - Stolbov, Leonid A. A1 - Veselovsky, Alexander V. A1 - De Rosa, Maria A1 - De Simone, Giada A1 - Gulotta, Maria R. A1 - Lombino, Jessica A1 - Mekni, Nedra A1 - Perricone, Ugo A1 - Casini, Arturo A1 - Embree, Amanda A1 - Gordon, D. Benjamin A1 - Lei, David A1 - Pratt, Katelin A1 - Voigt, Christopher A. A1 - Chen, Kuang-Yu A1 - Jacob, Yves A1 - Krischuns, Tim A1 - Lafaye, Pierre A1 - Zettor, Agnès A1 - Rodríguez, M. Luis A1 - White, Kris M. A1 - Fearon, Daren A1 - Von Delft, Frank A1 - Walsh, Martin A. A1 - Horvath, Dragos A1 - Brooks III, Charles L. A1 - Falsafi, Babak A1 - Ford, Bryan A1 - García-Sastre, Adolfo A1 - Yup Lee, Sang A1 - Naffakh, Nadia A1 - Varnek, Alexandre A1 - Klambauer, Günter A1 - Hermans, Thomas M. T1 - A community effort in SARS-CoV-2 drug discovery JF - Molecular Informatics KW - COVID-19 KW - drug discovery KW - machine learning KW - SARS-CoV-2 Y1 - 2023 U6 - https://doi.org/https://doi.org/10.1002/minf.202300262 VL - 43 IS - 1 SP - e202300262 ER - TY - JOUR A1 - Gorgulla, Christoph A1 - Nigam, AkshatKumar A1 - Koop, Matt A1 - Selim Çınaroğlu, Süleyman A1 - Secker, Christopher A1 - Haddadnia, Mohammad A1 - Kumar, Abhishek A1 - Malets, Yehor A1 - Hasson, Alexander A1 - Li, Minkai A1 - Tang, Ming A1 - Levin-Konigsberg, Roni A1 - Radchenko, Dmitry A1 - Kumar, Aditya A1 - Gehev, Minko A1 - Aquilanti, Pierre-Yves A1 - Gabb, Henry A1 - Alhossary, Amr A1 - Wagner, Gerhard A1 - Aspuru-Guzik, Alán A1 - Moroz, Yurii S. A1 - Fackeldey, Konstantin A1 - Arthanari, Haribabu T1 - VirtualFlow 2.0 - The Next Generation Drug Discovery Platform Enabling Adaptive Screens of 69 Billion Molecules JF - bioRxiv KW - preprint Y1 - 2023 U6 - https://doi.org/10.1101/2023.04.25.537981 ER - TY - JOUR A1 - Rodrigues Pela, Ronaldo A1 - Hsiao, Ching-Lien A1 - Hultman, Lars A1 - Birch, Jens A1 - Gueorguiev, Gueorgui Kostov T1 - Electronic and optical properties of core-shell InAlN nanorods: a comparative study via LDA, LDA-1/2, mBJ, HSE06, G0W0 and BSE methods JF - Physical Chemistry Chemical Physics Y1 - 2024 U6 - https://doi.org/10.1039/D3CP05295H VL - 26 SP - 7504 EP - 7514 PB - Royal Society of Chemistry ER - TY - JOUR A1 - Rodrigues Pela, Ronaldo A1 - Vona, Cecilia A1 - Lubeck, Sven A1 - Alex, Ben A1 - Gonzalez Oliva, Ignacio A1 - Draxl, Claudia T1 - Critical assessment of G0W0 calculations for 2D materials: the example of monolayer MoS2 JF - npj Comput. Mater. Y1 - 2024 U6 - https://doi.org/10.1038/s41524-024-01253-2 VL - 10 SP - 77 ER - TY - JOUR A1 - Trepte, Philipp A1 - Secker, Christopher A1 - Olivet, Julien A1 - Blavier, Jeremy A1 - Kostova, Simona A1 - Maseko, Sibusiso B A1 - Minia, Igor A1 - Silva Ramos, Eduardo A1 - Cassonnet, Patricia A1 - Golusik, Sabrina A1 - Zenkner, Martina A1 - Beetz, Stephanie A1 - Liebich, Mara J A1 - Scharek, Nadine A1 - Schütz, Anja A1 - Sperling, Marcel A1 - Lisurek, Michael A1 - Wang, Yang A1 - Spirohn, Kerstin A1 - Hao, Tong A1 - Calderwood, Michael A A1 - Hill, David E A1 - Landthaler, Markus A1 - Choi, Soon Gang A1 - Twizere, Jean-Claude A1 - Vidal, Marc A1 - Wanker, Erich E T1 - AI-guided pipeline for protein–protein interaction drug discovery identifies a SARS-CoV-2 inhibitor JF - Molecular Systems Biology N2 - Protein–protein interactions (PPIs) offer great opportunities to expand the druggable proteome and therapeutically tackle various diseases, but remain challenging targets for drug discovery. Here, we provide a comprehensive pipeline that combines experimental and computational tools to identify and validate PPI targets and perform early-stage drug discovery. We have developed a machine learning approach that prioritizes interactions by analyzing quantitative data from binary PPI assays or AlphaFold-Multimer predictions. Using the quantitative assay LuTHy together with our machine learning algorithm, we identified high-confidence interactions among SARS-CoV-2 proteins for which we predicted three-dimensional structures using AlphaFold-Multimer. We employed VirtualFlow to target the contact interface of the NSP10-NSP16 SARS-CoV-2 methyltransferase complex by ultra-large virtual drug screening. Thereby, we identified a compound that binds to NSP10 and inhibits its interaction with NSP16, while also disrupting the methyltransferase activity of the complex, and SARS-CoV-2 replication. Overall, this pipeline will help to prioritize PPI targets to accelerate the discovery of early-stage drug candidates targeting protein complexes and pathways. KW - Applied Mathematics KW - Computational Theory and Mathematics KW - General Agricultural and Biological Sciences KW - General Immunology and Microbiology KW - General Biochemistry, Genetics and Molecular Biology KW - Information Systems Y1 - 2024 U6 - https://doi.org/https://doi.org/10.1038/s44320-024-00019-8 SN - 1744-4292 VL - 20 IS - 4 SP - 428 EP - 457 PB - Springer Science and Business Media LLC ER - TY - JOUR A1 - Rodrigues Pelá, Ronaldo A1 - Draxl, Claudia T1 - Speeding up all-electron real-time TDDFT demonstrated by the exciting package JF - Comput. Phys. Commun. Y1 - 2024 U6 - https://doi.org/10.1016/j.cpc.2024.109292 VL - 304 SP - 109292 ER - TY - JOUR A1 - Fritsch, Daniel A1 - Schorr, Susan T1 - Disorder induced band gap lowering in kesterite type Cu2ZnSnSe4 and Ag2ZnSnSe4: a first-principles and special quasirandom structures investigation JF - Journal of Physics: Condensed Matter N2 - Abstract Quaternary chalcogenides, i.e. Cu2ZnSnS4, crystallising in the kesterite crystal structure have already been demonstrated as potential building blocks of thin film solar cells, containing only abundant elements and exhibiting power conversion efficiencies of about 14.9% so far. However, due to the potential presence of several structurally similar polymorphs, the unequivocal identification of their ground state crystal structures required the application of more elaborate neutron diffraction experiments. One particular complication arose from the later identified Cu–Zn disorder, present in virtually all thin film samples. Subsequently, it has been shown experimentally that this unavoidable Cu–Zn disorder leads to a band gap lowering in the respective samples. Additional theoretical investigations, mostly based on Monte-Carlo methods, tried to understand the atomistic origin of this disorder induced band gap lowering. Here, we present theoretical results from first-principles calculations based on density functional theory for the disorder induced band gap lowering in kesterite Cu2ZnSnSe4 and Ag2ZnSnSe4, where the Cu–Zn and Ag–Zn disorder is modelled via a supercell approach and special quasirandom structures. Results of subsequent analyses of structural, electronic, and optical properties are discussed with respect to available experimental results, and will provide additional insight and knowledge towards the atomistic origin of the observed disorder induced band gap lowering in kesterite type materials. Y1 - 2024 U6 - https://doi.org/10.1088/1361-648X/ad52de SN - 0953-8984 VL - 36 IS - 37 PB - IOP Publishing ER - TY - JOUR A1 - Groussin, Olivier A1 - Jorda, L. A1 - Attree, N. A1 - Birch, S.P.D. A1 - Bürger, Johanna A1 - Gutiérrez, Pedro A1 - Jindal, Abhinav S. A1 - Keller, H.U. A1 - Kramer, Tobias A1 - Lasagni Manghi, Riccardo A1 - Läuter, Matthias A1 - Markkanen, J. A1 - Marschall, R. A1 - Schuckart, C. T1 - Thermal environment and erosion of comet 67P/Churyumov-Gerasimenko JF - Astronomy and Astrophysics N2 - Aims. This paper focuses on how insolation affects the nucleus of comet 67P/Churyumov-Gerasimenko over its current orbit. We aim to better understand the thermal environment of the nucleus, in particular its surface temperature variations, erosion, relationship with topography, and how insolation affects the interior temperature for the location of volatile species (H2O and CO2). Methods. We have developed two thermal models to calculate the surface and subsurface temperatures of 67P over its 6.45-year orbit. The first model, with high resolution (300 000 facets), calculates surface temperatures, taking shadows and self-heating into account but ignoring thermal conductivity. The second model, with lower resolution (10 000 facets), includes thermal conductivity to estimate temperatures down to ∼3 m below the surface. Results. The thermal environment of 67P is strongly influenced by its large obliquity (52◦), which causes significant seasonal effects and polar nights. The northern hemisphere is the coldest region, with temperatures of 210–300 K. H2O is found in the first few centimetres, while CO2 is found deeper (∼2 m) except during polar night around perihelion, when CO2 accumulates near the surface. Cliffs erode 3–5 times faster than plains, forming terraces. The equatorial region receives maximum solar energy (8.5×109 J m−2 per orbit), with maximum surface temperatures of 300–350 K. On the plains, H2O is found in the first few centimetres, while CO2 is found deeper (∼2 m) and never accumulates near the surface. In the southern hemisphere, a brief intense perihelion heating raises temperatures to 350–400 K, which is followed by a 5-year polar night when surface temperatures drop to 55 K. Here H2O remains in the first few centimetres, while CO2 accumulates shallowly during polar night, enriching the region. Erosion is maximal in the southern hemisphere and concentrated on the plains, which explains the observed overall flatness of this hemisphere compared to the northern one. Over one orbit, the total energy from self-heating is 17% of the total energy budget, and 34% for thermal conduction. Our study contributes to a better understanding of the surface changes observed on 67P. Y1 - 2025 U6 - https://doi.org/10.1051/0004-6361/202452260 VL - 694 SP - A21 ER - TY - JOUR A1 - Straube, Arthur A1 - Olicón Méndez, Guillermo A1 - Winkelmann, Stefanie A1 - Höfling, Felix A1 - Engel, Maximilian T1 - Unfolding the geometric structure and multiple timescales of the urea-urease pH oscillator N2 - We study a two-variable dynamical system modeling pH oscillations in the urea–urease reaction within giant lipid vesicles – a problem that intrinsically contains multiple, well-separated timescales. Building on an existing, deterministic formulation via ordinary differential equations, we resolve different orders of magnitude within a small parameter and analyze the system's limit cycle behavior using geometric singular perturbation theory (GSPT). By introducing two different coordinate scalings – each valid in a distinct region of the phase space – we resolve the local dynamics near critical fold points, using the extension of GSPT through such singular points due to Krupa and Szmolyan. This framework enables a geometric decomposition of the periodic orbits into slow and fast segments and yields closed-form estimates for the period of oscillation. In particular, we link the existence of such oscillations to an underlying biochemical asymmetry, namely, the differential transport across the vesicle membrane. Y1 - 2025 ER - TY - JOUR A1 - Secker, Christopher A1 - Fackeldey, Konstantin A1 - Weber, Marcus A1 - Ray, Sourav A1 - Gorgulla, Christoph A1 - Schütte, Christof T1 - Novel multi-objective affinity approach allows to identify pH-specific μ-opioid receptor agonists JF - Journal of Cheminformatics N2 - Opioids are essential pharmaceuticals due to their analgesic properties, however, lethal side effects, addiction, and opioid tolerance are extremely challenging. The development of novel molecules targeting the μ-opioid receptor (MOR) in inflamed, but not in healthy tissue, could significantly reduce these unwanted effects. Finding such novel molecules can be achieved by maximizing the binding affinity to the MOR at acidic pH while minimizing it at neutral pH, thus combining two conflicting objectives. Here, this multi-objective optimal affinity approach is presented, together with a virtual drug discovery pipeline for its practical implementation. When applied to finding pH-specific drug candidates, it combines protonation state-dependent structure and ligand preparation with high-throughput virtual screening. We employ this pipeline to characterize a set of MOR agonists identifying a morphine-like opioid derivative with higher predicted binding affinities to the MOR at low pH compared to neutral pH. Our results also confirm existing experimental evidence that NFEPP, a previously described fentanyl derivative with reduced side effects, and recently reported β-fluorofentanyls and -morphines show an increased specificity for the MOR at acidic pH when compared to fentanyl and morphine. We further applied our approach to screen a >50K ligand library identifying novel molecules with pH-specific predicted binding affinities to the MOR. The presented differential docking pipeline can be applied to perform multi-objective affinity optimization to identify safer and more specific drug candidates at large scale. Y1 - 2023 U6 - https://doi.org/10.1186/s13321-023-00746-4 VL - 15 ER - TY - JOUR A1 - Höfling, Felix A1 - Dietrich, Siegfried T1 - Finite-size corrections for the static structure factor of a liquid slab with open boundaries JF - The Journal of Chemical Physics N2 - The presence of a confining boundary can modify the local structure of a liquid markedly. In addition, small samples of finite size are known to exhibit systematic deviations of thermodynamic quantities relative to their bulk values. Here, we consider the static structure factor of a liquid sample in slab geometry with open boundaries at the surfaces, which can be thought of as virtually cutting out the sample from a macroscopically large, homogeneous fluid. This situation is a relevant limit for the interpretation of grazing-incidence diffraction experiments at liquid interfaces and films. We derive an exact, closed expression for the slab structure factor, with the bulk structure factor as the only input. This shows that such free boundary conditions cause significant differences between the two structure factors, in particular, at small wavenumbers. An asymptotic analysis of this result yields the scaling exponent and an accurate, useful approximation of these finite-size corrections. Furthermore, the open boundaries permit the interpretation of the slab as an open system, supporting particle exchange with a reservoir. We relate the slab structure factor to the particle number fluctuations and discuss conditions under which the subvolume of the slab represents a grand canonical ensemble with chemical potential μ and temperature T. Thus, the open slab serves as a test-bed for the small-system thermodynamics in a μT reservoir. We provide a microscopically justified and exact result for the size dependence of the isothermal compressibility. Our findings are corroborated by simulation data for Lennard-Jones liquids at two representative temperatures. Y1 - 2020 U6 - https://doi.org/10.1063/5.0017923 VL - 153 SP - 054119 ER - TY - JOUR A1 - Ebrahimi Viand, Roya A1 - Höfling, Felix A1 - Klein, Rupert A1 - Delle Site, Luigi T1 - Theory and simulation of open systems out of equilibrium JF - The Journal of Chemical Physics N2 - We consider the theoretical model of Bergmann and Lebowitz for open systems out of equilibrium and translate its principles in the adaptive resolution simulation molecular dynamics technique. We simulate Lennard-Jones fluids with open boundaries in a thermal gradient and find excellent agreement of the stationary responses with the results obtained from the simulation of a larger locally forced closed system. The encouraging results pave the way for a computational treatment of open systems far from equilibrium framed in a well-established theoretical model that avoids possible numerical artifacts and physical misinterpretations. Y1 - 2020 U6 - https://doi.org/10.1063/5.0014065 VL - 153 SP - 101102 ER - TY - JOUR A1 - Stein, Lewin A1 - Straube, Florian A1 - Weinzierl, Stefan A1 - Lemke, Mathias T1 - Directional sound source modeling using the adjoint Euler equations in a finite-difference time-domain approach JF - Acoustical Society of America N2 - An adjoint-based approach for synthesizing complex sound sources by discrete, grid-based monopoles in finite-difference time-domain simulations is presented. Previously [Stein et al., 2019a, J. Acoust. Soc. Am. 146(3), 1774–1785] demonstrated that the approach allows to consider unsteady and non-uniform ambient conditions such as wind flow and thermal gradient in contrast to standard methods of numerical sound field simulation. In this work, it is proven that not only ideal monopoles but also realistic sound sources with complex directivity characteristics can be synthesized. In detail, an oscillating circular piston and a real 2-way near-field monitor are modeled. The required number of monopoles in terms of the SPL deviation between the directivity of the original and the synthesized source is analyzed. Since the computational effort is independent of the number of monopoles used for the synthesis, also more complex sources can be reproduced by increasing the number of monopoles utilized. In contrast to classical least-square problem solvers, this does not increase the computational effort, which makes the method attractive for predicting the effect of sound reinforcement systems with highly directional sources under difficult acoustic boundary conditions. Y1 - 2020 U6 - https://doi.org/https://doi.org/10.1121/10.0002425 VL - 148 IS - 5 SP - 3075 EP - 3085 PB - The Journal of the Acoustical Society of America ER -