TY - JOUR A1 - Charron, Nicholas A1 - Musil, Félix A1 - Guljas, Andrea A1 - Chen, Yaoyi A1 - Bonneau, Klara A1 - Pasos-Trejo, Aldo A1 - Jacopo, Venturin A1 - Daria, Gusew A1 - Zaporozhets, Iryna A1 - Krämer, Andreas A1 - Templeton, Clark A1 - Atharva, Kelkar A1 - Durumeric, Aleksander A1 - Olsson, Simon A1 - Pérez, Adrià A1 - Majewski, Maciej A1 - Husic, Brooke A1 - Patel, Ankit A1 - De Fabritiis, Gianni A1 - Noé, Frank A1 - Clementi, Cecilia T1 - Navigating protein landscapes with a machine-learned transferable coarse-grained model JF - Nature Chemistry N2 - The most popular and universally predictive protein simulation models employ all-atom molecular dynamics (MD), but they come at extreme computational cost. The development of a universal, computationally efficient coarse-grained (CG) model with similar prediction performance has been a long-standing challenge. By combining recent deep learning methods with a large and diverse training set of all-atom protein simulations, we here develop a bottom-up CG force field with chemical transferability, which can be used for extrapolative molecular dynamics on new sequences not used during model parametrization. We demonstrate that the model successfully predicts folded structures, intermediates, metastable folded and unfolded basins, and the fluctuations of intrinsically disordered proteins while it is several orders of magnitude faster than an all-atom model. This showcases the feasibility of a universal and computationally efficient machine-learned CG model for proteins. Y1 - 2025 U6 - https://doi.org/10.1038/s41557-025-01874-0 VL - 17 SP - 1284 EP - 1292 ER - TY - JOUR A1 - Krishnadas, Anirudh A1 - Charron, Nicolas A1 - Fournier, Rene T1 - Melting Transitions in Small Aluminum Clusters Simulated with Energies Approaching DFT Accuracy. JF - ACS Publications- Journal of Chemical Theory and Computation N2 - We describe a computational framework for modelling melting-like transitions in atomic clusters that combines first-principles energy calculations, global optimization, and machine-learned interatomic potentials. A diverse set of configurations is generated by global optimization, and energies are calculated by Density Functional Theory. The energies are fitted to an accuracy of 10 meV/atom or better with an Allegro E(3)-equivariant neural network potential. The model allows efficient parallel tempering Monte Carlo simulations with near DFT-level accuracy. This methodology was validated by simulating Na_20 and comparing it to earlier experimental and computational results. We used it to study melting-like transitions in Al_n+ clusters (n=9 to 16), and Al_n and Al_n^- (n=12, 13,14). The simulated heat capacity of these clusters, in particular Al_16+, are in qualitative agreement with experiments. The melting point of Al_n+ clusters with n=11-16 are well above the bulk melting point (934 K). The closed-shell Al_13- species has an exceptionally high melting point, close to 2100 K. Y1 - 2025 ER - TY - JOUR A1 - Krishnadas, Anirudh A1 - Moshi, Maryam A1 - Quintana, Ramon Alain Miranda A1 - Fournier, Rene T1 - Geometric Features and a Neural Network Classifier for Detecting Melting-Like Transitions in Clusters. JF - AIP Publications - The Journal of Chemical Physics N2 - Melting-like transitions in clusters are normally identified by a peak in the heat capacity curve C(T ) at T = Tc. Computing C(T ) requires costly simulations with millions of steps. We discuss four easily calculated functions of temperature that help detect and characterize melting-like transitions. The first, f1 (or WU ), is the width of the potential energy distribution which shows an abrupt increase near Tc. The other three are statistics of the ordered set of N(N − 1)/2 interatomic distances rij: (i) f2 is a measure of dissimilarity to the lowest energy configuration, or global minimum; (ii) f3 is the effective number of rij’s found in a small interval centered around (r1 + r2)/2 where r1, r2 are the positions of the first two peaks in the pair distribution function; and (iii) f4 is a measure of non-uniformity in the distribution of the ri j’s. Numerical tests with empirical potentials that model three types of bonding (van der Waals, covalent, and metallic) show that f1, f2, f3, and f4 produce estimates for the middle of the melting region in general agreement with Tc. An Artificial Neural Network (ANN) classifier that takes, as inputs, f2, f3, and many variants of f4, is used to calculate the solid fraction FS(T) and find the solid-liquid coexistence region between freezing and melting temperatures, [Tf , Tm]. Inflection points in f3(T ) and FS(T ) are very sensitive indicators of phase transitions. Estimates of Tc obtained from them converge one to three orders of magnitude faster, in simulation time, than those obtained with C(T ). Y1 - 2025 ER - TY - JOUR A1 - Clausecker, Robert A1 - Lemire, Daniel A1 - Schintke, Florian T1 - Faster Positional‐Population Counts for AVX2, AVX‐512, and ASIMD JF - Concurrency and Computation: Practice and Experience Y1 - 2025 U6 - https://doi.org/10.1002/cpe.70435 SN - 1532-0626 VL - 37 IS - 27-28 PB - Wiley ER - TY - JOUR A1 - Hölter, Arne A1 - Lemke, Mathias A1 - Weinzierl, Stefan A1 - Stein, Lewin T1 - Non-Reflecting Characteristic Boundary Conditions for Adjoint Time-Domain Acoustic Simulations JF - Journal of Theoretical and Computational Acoustics N2 - Accurate acoustic simulations in the free field require non-reflective boundary conditions to suppress spurious reflections at the computational domain boundaries. Although several characteristic-based formulations for direct (forward) simulations have been proposed in recent decades, the adjoint formulations of such characteristic-based boundary conditions (CBCs) have received limited atten- tion in the literature and lack a comprehensive analysis. This paper presents the derivation and evaluation of adjoint CBCs complementing the existing direct CBCs. Both the forward and adjoint CBCs are applied to the (nonlinear) Euler equations and linear acoustic equations in time-domain simulations. In this manner, the CBCs are investigated and subsequently compared to assess their respective accuracy and consistency. The CBCs were implemented using both a single-point and a zonal approach, with the former optionally combined with a sponge layer. Both approaches yielded comparable results in direct and adjoint simulations, while the zonal CBC exhibited improved ac- curacy at lower frequencies. Across the evaluated frequency range of approximately 350 – 5600 Hz, spurious reflections were attenuated by up to −70dB in both forward and adjoint cases, demon- strating the effectiveness and consistency of the proposed boundary treatment. Y1 - 2025 ER - TY - JOUR A1 - Höfling, Felix A1 - Straube, Arthur T1 - Langevin equations and a geometric integration scheme for the overdamped limit of rotational Brownian motion of axisymmetric particles JF - Physical Review Research N2 - The translational motion of anisotropic or self-propelled colloidal particles is closely linked with the particle’s orientation and its rotational Brownian motion. In the overdamped limit, the stochastic evolution of the orientation vector follows a diffusion process on the unit sphere and is characterized by an orientation-dependent (“multiplicative”) noise. As a consequence, the corresponding Langevin equation attains different forms depending on whether Itō’s or Stratonovich’s stochastic calculus is used. We clarify that both forms are equivalent and derive them in a top-down approach from a geometric construction of Brownian motion on the unit sphere, based on infinitesimal random rotations. Our approach suggests further a geometric integration scheme for rotational Brownian motion, which preserves the normalization constraint of the orientation vector exactly. We show that a simple implementation of the scheme, using Gaussian random rotations, converges weakly at order 1 of the integration time step, and we outline an advanced variant of the scheme that is weakly exact for an arbitrarily large time step. Due to a favorable prefactor of the discretization error, already the Gaussian scheme allows for integration time steps that are one order of magnitude larger compared to a commonly used algorithm for rotational Brownian dynamics simulations based on projection on the constraining manifold. For torques originating from constant external fields, we prove by virtue of the Fokker-Planck equation that the constructed diffusion process satisfies detailed balance and converges to the correct equilibrium distribution. The analysis is restricted to time-homogeneous rotational Brownian motion (i.e., a single rotational diffusion constant), which is relevant for axisymmetric particles and also chemically anisotropic spheres, such as self-propelled Janus particles. Y1 - 2025 U6 - https://doi.org/10.1103/wzdn-29p4 VL - 7 SP - 043034 ER - TY - JOUR A1 - Marcelli, Giorgia A1 - Bottinelli Montandon, Tecla A1 - Ebrahimi Viand, Roya A1 - Höfling, Felix T1 - Fluid flow inside slit-shaped nanopores: the role of surface morphology at the molecular scale JF - The Journal of Chemical Physics N2 - Non-equilibrium molecular dynamics (NEMD) simulations of fluid flow have highlighted the peculiarities of nanoscale flows compared to classical fluid mechanics; in particular, boundary conditions can deviate from the no-slip behavior at macroscopic scales. For fluid flow in slit-shaped nanopores, we demonstrate that surface morphology provides an efficient control on the slip length, which approaches zero when matching the molecular structures of the pore wall and the fluid. Using boundary-driven, energy-conserving NEMD simulations with a pump-like driving mechanism, we examine two types of pore walls—mimicking a crystalline and an amorphous material—that exhibit markedly different surface resistances to flow. The resulting flow velocity profiles are consistent with Poiseuille theory for incompressible, Newtonian fluids when adjusted for surface slip. For the two pores, we observe partial slip and no-slip behavior, respectively. The hydrodynamic permeability corroborates that the simulated flows are in the Darcy regime. However, the confinement of the fluid gives rise to an effective viscosity below its bulk value; wide pores exhibit a crossover between boundary and bulk-like flows. In addition, the thermal isolation of the flow causes a linear increase in fluid temperature along the flow, which we relate to strong viscous dissipation and heat convection, utilizing conservation laws of fluid mechanics. Noting that the investigated fluid model does not form droplets, our findings challenge the universality of previously reported correlations between slippage, solvophobicity, and a depletion zone. Furthermore, they underscore the need for molecular-scale modeling to accurately capture the fluid dynamics near boundaries and in nanoporous materials, where macroscopic models may not be applicable. Y1 - 2025 U6 - https://doi.org/10.1063/5.0246573 VL - 162 SP - 104101 ER - TY - JOUR A1 - Zhong, Xin A1 - Höfling, Felix A1 - John, Timm T1 - Hydrogen diffusion in garnet: insights from atomistic simulations JF - Geochemistry Geophysics Geosystems N2 - Garnet has been widely used to decipher the pressure-temperature-time history of rocks, but its physical properties such as elasticity and diffusion are strongly affected by trace amounts of hydrogen. Experimental measurements of H diffusion in garnet are limited to room pressure. We use atomistic simulations to study H diffusion in perfect and defective garnet lattices, focusing on protonation defects at the Si and Mg sites, which are shown to be energetically favored. Transient trapping of H renders ab-initio simulations of H diffusion computationally challenging, which is overcome with machine learning techniques by training a deep neural network that encodes the interatomic potential. Our results from such deep potential molecular dynamics (DeePMD) simulations show high mobility of hydrogen in defect-free garnet lattices, whereas H diffusivity is significantly diminished in defective lattices. Tracer simulations focusing on H alone highlight the vital role of atomic vibrations of heavier atoms like Mg on the release of H atoms. Two regimes of H diffusion are identified: a diffuser-dominated regime at high hydrogen content with low activation energies due to saturation of vacancies by hydrogen, and a vacancy-dominated regime at low hydrogen content with high activation energies due to trapping of H atoms at vacancy sites. These regimes account for experimental observations, such as a H-concentration dependent diffusivity and the discrepancy in activation energy between deprotonation and D-H exchange experiments. This study underpins the crucial role of vacancies in H diffusion and demonstrates the utility of machine-learned interatomic potentials in studying kinetic processes in the Earth's interior. Y1 - 2025 U6 - https://doi.org/10.1029/2024GC011951 VL - 26 SP - e2024GC011951 ER - TY - CHAP A1 - Christgau, Steffen A1 - Dylan, Everingham A1 - Lübke, Max A1 - De Lucia, Marco A1 - Puhan, Danny A1 - Schelten, Niklas A1 - Schnor, Bettina A1 - Signer, Hannes A1 - Spazier, Johannes A1 - Stabernack, Benno A1 - Steinert, Fritjof A1 - Yahdzhyiev, Serhii T1 - On the Usability and Energy Efficiency of High-Level Synthesis for FPGA-based Network-Attached Accelerators T2 - 2025 IEEE International Parallel and Distributed Processing Symposium Workshops (IPDPSW) N2 - Heterogeneity in high performance computing systems is one of the most promising approaches towards more energy-efficient computing on one hand and satisfying the raising demand of global computation capacity on the other hand. Besides the well-known key components like CPUs and GPGPUs are domain-specific accelerators like TPUs, FPGAs well known for their energy efficiency. This is especially true for highly specialized use cases. Network-attached accelerators promise more scalability and flexibility for FPGA usage in HPC environments. Easy and efficient programming of those accelerators is, however, still an open issue. Based on a framework for such accelerators which enables decoupling of FPGAs from their host system, we present a workflow using High-Level Synthesis (HLS) to offload application kernels to them. We evaluate this approach against a conventional Hardware Description Language (HDL) based workflow. In addition, we introduce the energy measurement tool EMA and assess the energy efficiency of both HLS and HDL design. Y1 - 2025 U6 - https://doi.org/10.1109/IPDPSW66978.2025.00139 SP - 886 EP - 895 ER - TY - JOUR A1 - Läuter, Matthias A1 - Kramer, Tobias T1 - Rotation dynamics and torque efficiency of cometary nuclei JF - Astronomy & Astrophysics N2 - The dynamics of a rigid cometary nucleus is described by the evolutions of its center-of-mass and of its rotation state. Solar irradiation that reaches the surface of a cometary nucleus causes the sublimation of volatiles that form the coma around the nucleus. The sublimation process transfers linear momentum and rotational angular momentum from the nucleus to the surrounding space, and thus affects the dynamics via nongravitational forces and nongravitational torques. With the exception of close approaches to planets, these torques exert the dominant influence on the rotation states of cometary nuclei. The Rosetta mission 2014-2016 accompanying comet 67P/Churyumov-Gerasimenko provides the longest continuous observational data to track its rotation state. In particular, the data set encompasses the direction of the angular velocity, denoted by ω, and the angular frequency |ω|over a time period of approximately 700 days. The observed change of the rotation state is not explained by a low heat conductivity thermophysical model in combination with a homogeneous surface ice coverage of comet 67P. Spatially and/or temporally varying weights for effective active fraction with respect to a prescribed set of surface regions provide a potential solution to this problem. Here, we present a methodology for classifying the surface based on vectorial efficiency of the torque. On any cometary surface without geometric symmetry, the methodology highlights the decomposition into eight characteristic regions that encode the signs of torque efficiency with respect to all vector components. This decomposition is divided into two subsets of four regions each of which is located in one of both hemispheric regions. We analyze in detail rotation states close to lowest energy and different thermophysical models, and we discuss how the uncertainties of observations affect the model parameters. We study the occurrence of these regions for an oblate ellipsoid, a nearprolate ellipsoid, a bilobed shape, and a shape model analogous to that of comet 67P. The sensitivity analysis for comet 67P indicates that the observations constrain only one of the eight weights uniquely. The other directions are poorly constrained and show the limitation of the rotational data to determine the regional activity on comet 67P. Y1 - 2025 U6 - https://doi.org/10.1051/0004-6361/202553845 VL - 699 SP - A75 ER - TY - THES A1 - el Bahnasawy, Henry T1 - Konzeption und Entwicklung einer Plattform zur geographisch verteilten Bereitstellung von latenzkritischen Anwendungen Y1 - 2025 ER - TY - JOUR A1 - Rossi, Thomas A1 - Qiao, Lu A1 - Dykstra, Conner P. A1 - Rodrigues Pelá, Ronaldo A1 - Gnewkow, Richard A1 - Wallick, Rachel F. A1 - Burke, John H. A1 - Nicholas, Erin A1 - March, Anne-Marie A1 - Doumy, Gilles A1 - Buchholz, D. Bruce A1 - Deparis, Christiane A1 - Zuñiga-Pérez, Jesus A1 - Weise, Michael A1 - Ellmer, Klaus A1 - Fondell, Mattis A1 - Draxl, Claudia A1 - van der Veen, Renske T1 - Dynamic control of X-ray core-exciton resonances by Coulomb screening in photoexcited semiconductors JF - Communications Materials Y1 - 2025 U6 - https://doi.org/10.1038/s43246-025-00909-w VL - 6 SP - 191 ER - TY - JOUR A1 - Hawthorne, Felipe A1 - Raulino, Paulo R. E. A1 - Rodrigues Pelá, Ronaldo A1 - Woellner, Cristiano F. T1 - Efficient and Accurate Machine Learning Interatomic Potential for Graphene: Capturing Stress–Strain and Vibrational Properties JF - The Journal of Physical Chemistry C Y1 - 2025 U6 - https://doi.org/10.1021/acs.jpcc.5c03470 VL - 129 SP - 16319 EP - 16326 ER - TY - CHAP A1 - Lößer, Ansgar A1 - Witzke, Joel A1 - Schintke, Florian A1 - Scheuermann, Björn T1 - BottleMod: Modeling Data Flows and Tasks for Fast Bottleneck Analysis T2 - Proceedings of the 16th ACM/SPEC International Conference on Performance Engineering (ICPE 2025) Y1 - 2025 U6 - http://nbn-resolving.de/urn/resolver.pl?https://doi.org/10.1145/3676151.3719382 SP - 11 EP - 18 PB - Association for Computing Machinery ER -