TY - JOUR A1 - Delle Site, Luigi A1 - Krekeler, Christian A1 - Whittaker, John A1 - Agarwal, Animesh A1 - Klein, Rupert A1 - Höfling, Felix T1 - Molecular Dynamics of Open Systems: Construction of a Mean‐Field Particle Reservoir JF - Advanced Theory and Simulations N2 - The simulation of open molecular systems requires explicit or implicit reservoirs of energy and particles. Whereas full atomistic resolution is desired in the region of interest, there is some freedom in the implementation of the reservoirs. Here, a combined, explicit reservoir is constructed by interfacing the atomistic region with regions of point-like, non-interacting particles (tracers) embedded in a thermodynamic mean field. The tracer molecules acquire atomistic resolution upon entering the atomistic region and equilibrate with this environment, while atomistic molecules become tracers governed by an effective mean-field potential after crossing the atomistic boundary. The approach is extensively tested on thermodynamic, structural, and dynamic properties of liquid water. Conceptual and numerical advantages of the procedure as well as new perspectives are highlighted and discussed. Y1 - 2019 U6 - https://doi.org/10.1002/adts.201900014 VL - 2 SP - 1900014 ER - TY - JOUR A1 - Dibak, Manuel A1 - Fröhner, Christoph A1 - Noé, Frank A1 - Höfling, Felix T1 - Diffusion-influenced reaction rates in the presence of pair interactions JF - The Journal of Chemical Physics N2 - The kinetics of bimolecular reactions in solution depends, among other factors, on intermolecular forces such as steric repulsion or electrostatic interaction. Microscopically, a pair of molecules first has to meet by diffusion before the reaction can take place. In this work, we establish an extension of Doi’s volume reaction model to molecules interacting via pair potentials, which is a key ingredient for interacting-particle-based reaction–diffusion (iPRD) simulations. As a central result, we relate model parameters and macroscopic reaction rate constants in this situation. We solve the corresponding reaction–diffusion equation in the steady state and derive semi- analytical expressions for the reaction rate constant and the local concentration profiles. Our results apply to the full spectrum from well-mixed to diffusion-limited kinetics. For limiting cases, we give explicit formulas, and we provide a computationally inexpensive numerical scheme for the general case, including the intermediate, diffusion-influenced regime. The obtained rate constants decompose uniquely into encounter and formation rates, and we discuss the effect of the potential on both subprocesses, exemplified for a soft harmonic repulsion and a Lennard-Jones potential. The analysis is complemented by extensive stochastic iPRD simulations, and we find excellent agreement with the theoretical predictions. Y1 - 2019 U6 - https://doi.org/10.1063/1.5124728 VL - 151 SP - 164105 ER - TY - JOUR A1 - Straube, Arthur A1 - Kowalik, Bartosz G. A1 - Netz, Roland R. A1 - Höfling, Felix T1 - Rapid onset of molecular friction in liquids bridging between the atomistic and hydrodynamic pictures JF - Commun. Phys. N2 - Friction in liquids arises from conservative forces between molecules and atoms. Although the hydrodynamics at the nanoscale is subject of intense research and despite the enormous interest in the non-Markovian dynamics of single molecules and solutes, the onset of friction from the atomistic scale so far could not be demonstrated. Here, we fill this gap based on frequency-resolved friction data from high-precision simulations of three prototypical liquids, including water. Combining with theory, we show that friction in liquids emerges abruptly at a characteristic frequency, beyond which viscous liquids appear as non-dissipative, elastic solids. Concomitantly, the molecules experience Brownian forces that display persistent correlations. A critical test of the generalised Stokes–Einstein relation, mapping the friction of single molecules to the visco-elastic response of the macroscopic sample, disproves the relation for Newtonian fluids, but substantiates it exemplarily for water and a moderately supercooled liquid. The employed approach is suitable to yield insights into vitrification mechanisms and the intriguing mechanical properties of soft materials. Y1 - 2020 U6 - https://doi.org/10.1038/s42005-020-0389-0 VL - 3 SP - 126 PB - Nature ER - TY - GEN A1 - Mokhtari, Zahra A1 - Patterson, Robert I. A. A1 - Höfling, Felix T1 - Spontaneous trail formation in populations of auto-chemotactic walkers N2 - We study the formation of trails in populations of self-propelled agents that make oriented deposits of pheromones and also sense such deposits to which they then respond with gradual changes of their direction of motion. Based on extensive off-lattice computer simulations aiming at the scale of insects, e.g., ants, we identify a number of emerging stationary patterns and obtain qualitatively the non-equilibrium \add{state} diagram of the model, spanned by the strength of the agent--pheromone interaction and the number density of the population. In particular, we demonstrate the spontaneous formation of persistent, macroscopic trails, and highlight some behaviour that is consistent with a dynamic phase transition. This includes a characterisation of the mass of system-spanning trails as a potential order parameter. We also propose a dynamic model for a few macroscopic observables, including the sub-population size of trail-following agents, which captures the early phase of trail formation. T3 - ZIB-Report - 21-36 Y1 - 2021 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-84466 SN - 1438-0064 ER - TY - JOUR A1 - Pathania, Yogyata A1 - Chakraborty, Dipanjan A1 - Höfling, Felix T1 - Continuous demixing transition of binary liquids: finite-size scaling from the analysis of sub-systems JF - Advances Theory and Simulation Y1 - 2021 U6 - https://doi.org/10.1002/adts.202000235 VL - 4 SP - 2000235 ER - TY - JOUR A1 - Frömberg, Daniela A1 - Höfling, Felix T1 - Generalized master equation for first-passage problems in partitioned spaces JF - Journal of Physics A: Mathematical and Theoretical Y1 - 2021 U6 - https://doi.org/10.1088/1751-8121/abf2ec VL - 54 SP - 215601 ER - TY - JOUR A1 - Klein, Rupert A1 - Ebrahimi Viand, Roya A1 - Höfling, Felix A1 - Delle Site, Luigi T1 - Nonequilibrium induced by reservoirs: Physico-mathematical model and numerical tests JF - Advances Theory and Simulation Y1 - 2021 U6 - https://doi.org/10.1002/adts.202100071 VL - 4 SP - 2100071 ER - TY - JOUR A1 - Gholami, Abbas A1 - Höfling, Felix A1 - Klein, Rupert A1 - Delle Site, Luigi T1 - Thermodynamic relations at the coupling boundary in adaptive resolution simulations for open systems JF - Advances Theory and Simulation Y1 - 2021 U6 - https://doi.org/10.1002/adts.202000303 VL - 4 SP - 2000303 ER - TY - JOUR A1 - Mokhtari, Zahra A1 - Patterson, Robert I. A. A1 - Höfling, Felix T1 - Spontaneous trail formation in populations of auto-chemotactic walkers JF - New Journal of Physics Y1 - 2022 U6 - https://doi.org/10.1088/1367-2630/ac43ec VL - 24 SP - 013012 ER - TY - JOUR A1 - Höfling, Felix A1 - Dietrich, Siegfried T1 - Structure of liquid–vapor interfaces: perspectives from liquid state theory, large-scale simulations, and potential grazing-incidence X-ray diffraction JF - The Journal of Chemical Physics N2 - Grazing-incidence X-ray diffraction (GIXRD) is a scattering technique which allows one to characterize the structure of fluid interfaces down to the molecular scale, including the measurement of the surface tension and of the interface roughness. However, the corresponding standard data analysis at non-zero wave numbers has been criticized as to be inconclusive because the scattering intensity is polluted by the unavoidable scattering from the bulk. Here we overcome this ambiguity by proposing a physically consistent model of the bulk contribution which is based on a minimal set of assumptions of experimental relevance. To this end, we derive an explicit integral expression for the background scattering, which can be determined numerically from the static structure factors of the coexisting bulk phases as independent input. Concerning the interpretation of GIXRD data inferred from computer simulations, we account also for the finite sizes of the bulk phases, which are unavoidable in simulations. The corresponding leading-order correction beyond the dominant contribution to the scattered intensity is revealed by asymptotic analysis, which is characterized by the competition between the linear system size and the X-ray penetration depth in the case of simulations. Specifically, we have calculated the expected GIXRD intensity for scattering at the planar liquid--vapor interface of Lennard-Jones fluids with truncated pair interactions via extensive, high-precision simulations. The reported data cover interfacial and bulk properties of fluid states along the whole liquid--vapor coexistence line. A sensitivity analysis demonstrates the robustness of our findings concerning the detailed definition of the mean interface position. We conclude that previous claims of an enhanced surface tension at mesoscopic scales are amenable to unambiguous tests via scattering experiments. Y1 - 2024 U6 - https://doi.org/10.1063/5.0186955 VL - 160 SP - 104107 ER - TY - CHAP A1 - Skoblin, Viktor A1 - Höfling, Felix A1 - Christgau, Steffen T1 - Gaining Cross-Platform Parallelism for HAL’s Molecular Dynamics Package using SYCL T2 - 29. PARS-Workshop 2023 N2 - Molecular dynamics simulations are one of the methods in scientific computing that benefit from GPU acceleration. For those devices, SYCL is a promising API for writing portable codes. In this paper, we present the case study of HAL’s MD package that has been successfully migrated from CUDA to SYCL. We describe the different strategies that we followed in the process of porting the code. Following these strategies, we achieved code portability across major GPU vendors. Depending on the actual kernels, both significant performance improvements and regressions are observed. As a side effect of the migration process, we obtained impressing speedups also for execution on CPUs. Y1 - 2023 SN - 0177-0454 VL - 36 ER -