-
Gauge Invariance of Equilibrium Statistical Mechanics
Authors:
Johanna Müller,
Sophie Hermann,
Florian Sammüller,
Matthias Schmidt
Abstract:
We identify a recently proposed shifting operation on classical phase space as a gauge transformation for statistical mechanical microstates. The infinitesimal generators of the continuous gauge group form a non-commutative Lie algebra, which induces exact sum rules when thermally averaged. Gauge invariance with respect to finite shifting is demonstrated via Monte Carlo simulation in the transform…
▽ More
We identify a recently proposed shifting operation on classical phase space as a gauge transformation for statistical mechanical microstates. The infinitesimal generators of the continuous gauge group form a non-commutative Lie algebra, which induces exact sum rules when thermally averaged. Gauge invariance with respect to finite shifting is demonstrated via Monte Carlo simulation in the transformed phase space which generates identical equilibrium averages. Our results point towards a deeper basis of statistical mechanics than previously known and they offer avenues for systematic construction of exact identities and of sampling algorithms.
△ Less
Submitted 27 June, 2024;
originally announced June 2024.
-
Neural force functional for non-equilibrium many-body colloidal systems
Authors:
Toni Zimmerman,
Florian Sammüller,
Sophie Hermann,
Matthias Schmidt,
Daniel de las Heras
Abstract:
We combine power functional theory and machine learning to study non-equilibrium overdamped many-body systems of colloidal particles at the level of one-body fields. We first sample in steady state the one-body fields relevant for the dynamics from computer simulations of Brownian particles under the influence of randomly generated external fields. A neural network is then trained with this data t…
▽ More
We combine power functional theory and machine learning to study non-equilibrium overdamped many-body systems of colloidal particles at the level of one-body fields. We first sample in steady state the one-body fields relevant for the dynamics from computer simulations of Brownian particles under the influence of randomly generated external fields. A neural network is then trained with this data to represent locally in space the formally exact functional mapping from the one-body density and velocity profiles to the one-body internal force field. The trained network is used to analyse the non-equilibrium superadiabatic force field and the transport coefficients such as shear and bulk viscosities. Due to the local learning approach, the network can be applied to systems much larger than the original simulation box in which the one-body fields are sampled. Complemented with the exact non-equilibrium one-body force balance equation and a continuity equation, the network yields viable predictions of the dynamics in time-dependent situations. Even though training is based on steady states only, the predicted dynamics is in good agreement with simulation results. A neural dynamical density functional theory can be straightforwardly implemented as a limiting case in which the internal force field is that of an equilibrium system. The framework is general and directly applicable to other many-body systems of interacting particles following Brownian dynamics.
△ Less
Submitted 5 June, 2024;
originally announced June 2024.
-
PHOBIC: Perfect Hashing with Optimized Bucket Sizes and Interleaved Coding
Authors:
Stefan Hermann,
Hans-Peter Lehmann,
Giulio Ermanno Pibiri,
Peter Sanders,
Stefan Walzer
Abstract:
A minimal perfect hash function (MPHF) maps a set of n keys to {1, ..., n} without collisions. Such functions find widespread application e.g. in bioinformatics and databases. In this paper we revisit PTHash - a construction technique particularly designed for fast queries. PTHash distributes the input keys into small buckets and, for each bucket, it searches for a hash function seed that places i…
▽ More
A minimal perfect hash function (MPHF) maps a set of n keys to {1, ..., n} without collisions. Such functions find widespread application e.g. in bioinformatics and databases. In this paper we revisit PTHash - a construction technique particularly designed for fast queries. PTHash distributes the input keys into small buckets and, for each bucket, it searches for a hash function seed that places its keys in the output domain without collisions. The collection of all seeds is then stored in a compressed way. Since the first buckets are easier to place, buckets are considered in non-increasing order of size. Additionally, PTHash heuristically produces an imbalanced distribution of bucket sizes by distributing 60% of the keys into 30% of the buckets. Our main contribution is to characterize, up to lower order terms, an optimal distribution of expected bucket sizes. We arrive at a simple, closed form solution which improves construction throughput for space efficient configurations in practice. Our second contribution is a novel encoding scheme for the seeds. We split the keys into partitions. Within each partition, we run the bucket distribution and search step. We then store the seeds in an interleaved way by consecutively placing the seeds for the i-th buckets from all partitions. The seeds for the i-th bucket of each partition follow the same statistical distribution. This allows us to tune a compressor for each bucket. Hence, we call our technique PHOBIC - Perfect Hashing with Optimized Bucket sizes and Interleaved Coding. Compared to PTHash, PHOBIC is 0.17 bits/key more space efficient for same query time and construction throughput. We also contribute a GPU implementation to further accelerate MPHF construction. For a configuration with fast queries, PHOBIC-GPU can construct a perfect hash function at 2.17 bits/key in 28 ns per key, which can be queried in 37 ns on the CPU.
△ Less
Submitted 29 April, 2024;
originally announced April 2024.
-
Design and experimental validation of a finite-size labyrinthine metamaterial for vibro-acoustics: enabling upscaling towards large-scale structures
Authors:
Svenja Hermann,
Kévin Billon,
Alina-Michelle Parlak,
Jeanette Orlowsky,
Manuel Collet,
Angela Madeo
Abstract:
In this paper, we present the design and experimental validation of a labyrinthine metamaterial for vibro-acoustic applications. Based on a 2D unit cell, different designs of finite-size metamaterial specimens in a sandwich configuration including two plates are proposed. The design phase includes an optimisation based on Bloch-Floquet analysis with the aims of maximising the band gap and extrudin…
▽ More
In this paper, we present the design and experimental validation of a labyrinthine metamaterial for vibro-acoustic applications. Based on a 2D unit cell, different designs of finite-size metamaterial specimens in a sandwich configuration including two plates are proposed. The design phase includes an optimisation based on Bloch-Floquet analysis with the aims of maximising the band gap and extruding the specimens in the third dimension while keeping the absorption properties almost unaffected. By manufacturing and experimentally testing finite-sized specimens, we assess their capacity to mitigate vibrations in vibro-impact tests. The experiments confirm a band gap in the low- to mid-frequency range. Numerical models are employed to validate the experiments and to examine additional vibro-acoustic load cases. The metamaterial's performances are compared to benchmark solutions, usually employed for noise and vibration mitigation, showing a comparable efficacy in the band gap region. To eventually improve the metamaterial's performance, we optimise its interaction with the air and test different types of connections between the metamaterial and the homogeneous plates. This finally leads to metamaterial samples largely exceeding the benchmark performances in the band gap region and reveal the potential of interfaces for performance optimisation of composed structures.
△ Less
Submitted 16 April, 2024;
originally announced April 2024.
-
Hyper-density functional theory of soft matter
Authors:
Florian Sammüller,
Silas Robitschko,
Sophie Hermann,
Matthias Schmidt
Abstract:
We present a scheme for investigating arbitrary thermal observables in spatially inhomogeneous equilibrium many-body systems. Extending the grand canonical ensemble yields any given observable as an explicit hyper-density functional. Associated local fluctuation profiles follow from an exact hyper-Ornstein-Zernike equation. While the local compressibility and simple observables permit analytic tre…
▽ More
We present a scheme for investigating arbitrary thermal observables in spatially inhomogeneous equilibrium many-body systems. Extending the grand canonical ensemble yields any given observable as an explicit hyper-density functional. Associated local fluctuation profiles follow from an exact hyper-Ornstein-Zernike equation. While the local compressibility and simple observables permit analytic treatment, complex order parameters are accessible via simulation-based supervised machine learning of neural hyper-direct correlation functionals. We exemplify efficient and accurate neural predictions for the cluster statistics of hard rods, square well rods, and hard spheres. The theory allows to treat complex observables, as is impossible in standard density functional theory.
△ Less
Submitted 25 July, 2024; v1 submitted 12 March, 2024;
originally announced March 2024.
-
Experimental analysis of the TRC benchmark system
Authors:
Arati Bhattu,
Svenja Hermann,
Nidhal Jamia,
Florian Müller,
Maren Scheel,
Christoph Schwingshackl,
H. Nevzat Özgüven,
Malte Krack
Abstract:
The Tribomechadynamics Research Challenge (TRC) was a blind prediction of the vibration behavior of a thin plate clamped on two sides using bolted joints. The first bending mode's natural frequency and damping ratio were requested as function of the amplitude, starting from the linear regime until high levels, where both frictional contact and nonlinear bending-stretching coupling become relevant.…
▽ More
The Tribomechadynamics Research Challenge (TRC) was a blind prediction of the vibration behavior of a thin plate clamped on two sides using bolted joints. The first bending mode's natural frequency and damping ratio were requested as function of the amplitude, starting from the linear regime until high levels, where both frictional contact and nonlinear bending-stretching coupling become relevant. The predictions were confronted with experimental results in a companion paper; the present article addresses the experimental analysis of this benchmark system. Amplitude-dependent modal data was obtained from phase resonance and response controlled tests. An original variant of response controlled testing is proposed: Instead of a fixed frequency interval, a fixed phase interval is analyzed. This way, the high excitation levels required outside resonance, which could activate unwanted exciter nonlinearity, are avoided. Consistency of testing methods is carefully analyzed. Overall, these measures have permitted to gain high confidence in the acquired modal data. The different sources of the remaining uncertainty were further analyzed. A low reassembly-variability but a moderate time-variability were identified, where the latter is attributed to some thermal sensitivity of the system. Two nominally identical plates were analyzed, which both have an appreciable initial curvature, and a significant effect on the vibration behavior was found depending on whether the plate is aligned/misaligned with the support structure. Further, a 1:2 nonlinear modal interaction with the first torsion mode was observed, which only occurs in the aligned configurations.
△ Less
Submitted 12 March, 2024;
originally announced March 2024.
-
Noether invariance theory for the equilibrium force structure of soft matter
Authors:
Sophie Hermann,
Florian Sammüller,
Matthias Schmidt
Abstract:
We give details and derivations for the Noether invariance theory that characterizes the spatial equilibrium structure of inhomogeneous classical many-body systems, as recently proposed and investigated for bulk systems [F. Sammüller $\textit{et al.}$, Phys. Rev. Lett. $\textbf{130}$, 268203 (2023)]. Thereby an intrinsic thermal symmetry against a local shifting transformation on phase space is ex…
▽ More
We give details and derivations for the Noether invariance theory that characterizes the spatial equilibrium structure of inhomogeneous classical many-body systems, as recently proposed and investigated for bulk systems [F. Sammüller $\textit{et al.}$, Phys. Rev. Lett. $\textbf{130}$, 268203 (2023)]. Thereby an intrinsic thermal symmetry against a local shifting transformation on phase space is exploited on the basis of the Noether theorem for invariant variations. We consider the consequences of the shifting that emerge at second order in the displacement field that parameterizes the transformation. In a natural way the standard two-body density distribution is generated. Its second spatial derivative (Hessian) is thereby balanced by two further and different two-body correlation functions, which respectively introduce thermally averaged force correlations and force gradients in a systematic and microscopically sharp way into the framework. Separate exact self and distinct sum rules hold expressing this balance. We exemplify the validity of the theory on the basis of computer simulations for the Lennard-Jones gas, liquid, and crystal, the Weeks-Chandler-Andersen fluid, monatomic Molinero-Moore water at ambient conditions, a three-body-interacting colloidal gel former, the Yukawa and soft-sphere dipolar fluids, and for isotropic and nematic phases of Gay-Berne particles. We describe explicitly the derivation of the sum rules based on Noether's theorem and also give more elementary proofs based on partial phase space integration following Yvon's theorem.
△ Less
Submitted 23 March, 2024; v1 submitted 26 January, 2024;
originally announced January 2024.
-
Why neural functionals suit statistical mechanics
Authors:
Florian Sammüller,
Sophie Hermann,
Matthias Schmidt
Abstract:
We describe recent progress in the statistical mechanical description of many-body systems via machine learning combined with concepts from density functional theory and many-body simulations. We argue that the neural functional theory by Sammüller et al. [Proc. Nat. Acad. Sci. 120, e2312484120 (2023)] gives a functional representation of direct correlations and of thermodynamics that allows for t…
▽ More
We describe recent progress in the statistical mechanical description of many-body systems via machine learning combined with concepts from density functional theory and many-body simulations. We argue that the neural functional theory by Sammüller et al. [Proc. Nat. Acad. Sci. 120, e2312484120 (2023)] gives a functional representation of direct correlations and of thermodynamics that allows for thorough quality control and consistency checking of the involved methods of artificial intelligence. Addressing a prototypical system we here present a pedagogical application to hard core particle in one spatial dimension, where Percus' exact solution for the free energy functional provides an unambiguous reference. A corresponding standalone numerical tutorial that demonstrates the neural functional concepts together with the underlying fundamentals of Monte Carlo simulations, classical density functional theory, machine learning, and differential programming is available online at https://github.com/sfalmo/NeuralDFT-Tutorial.
△ Less
Submitted 14 February, 2024; v1 submitted 7 December, 2023;
originally announced December 2023.
-
Hyperforce balance via thermal Noether invariance of any observable
Authors:
Silas Robitschko,
Florian Sammüller,
Matthias Schmidt,
Sophie Hermann
Abstract:
Noether invariance in statistical mechanics provides fundamental connections between the symmetries of a physical system and its conservation laws and sum rules. The latter are exact identities that involve statistically averaged forces and force correlations and they are derived from statistical mechanical functionals. However, the implications for more general observables and order parameters ar…
▽ More
Noether invariance in statistical mechanics provides fundamental connections between the symmetries of a physical system and its conservation laws and sum rules. The latter are exact identities that involve statistically averaged forces and force correlations and they are derived from statistical mechanical functionals. However, the implications for more general observables and order parameters are unclear. Here, we demonstrate that thermally averaged classical phase space functions are associated with exact hyperforce sum rules that follow from translational Noether invariance. Both global and locally resolved identities hold and they relate the mean gradient of a phase-space function to its negative mean product with the total force. Similar to Hirschfelder's hypervirial theorem, the hyperforce sum rules apply to arbitrary observables in equilibrium. Exact hierarchies of higher-order sum rules follow iteratively. As applications we investigate via computer simulations the emerging one-body force fluctuation profiles in confined liquids. These local correlators quantify spatially inhomogeneous self-organization and their measurement allows for the development of stringent convergence tests and enhanced sampling schemes in complex systems.
△ Less
Submitted 7 February, 2024; v1 submitted 23 August, 2023;
originally announced August 2023.
-
Active crystallization from power functional theory
Authors:
Sophie Hermann,
Matthias Schmidt
Abstract:
We address the gas, liquid, and crystal phase behaviour of active Brownian particles in three dimensions. The nonequilibrium force balance at coexistence leads to equality of state functions for which we use power functional approximations. Motility-induced phase separation starts at a critical point and quickly becomes metastable against active freezing for Péclet numbers above a nonequilibrium t…
▽ More
We address the gas, liquid, and crystal phase behaviour of active Brownian particles in three dimensions. The nonequilibrium force balance at coexistence leads to equality of state functions for which we use power functional approximations. Motility-induced phase separation starts at a critical point and quickly becomes metastable against active freezing for Péclet numbers above a nonequilibrium triple point. The mean swim speed acts as a state variable, similar to the density of depletion agents in colloidal demixing. We obtain agreement with recent simulation results and correctly predict the strength of particle number fluctuations in active fluids.
△ Less
Submitted 5 February, 2024; v1 submitted 21 August, 2023;
originally announced August 2023.
-
Neural functional theory for inhomogeneous fluids: Fundamentals and applications
Authors:
Florian Sammüller,
Sophie Hermann,
Daniel de las Heras,
Matthias Schmidt
Abstract:
We present a hybrid scheme based on classical density functional theory and machine learning for determining the equilibrium structure and thermodynamics of inhomogeneous fluids. The exact functional map from the density profile to the one-body direct correlation function is represented locally by a deep neural network. We substantiate the general framework for the hard sphere fluid and use grand…
▽ More
We present a hybrid scheme based on classical density functional theory and machine learning for determining the equilibrium structure and thermodynamics of inhomogeneous fluids. The exact functional map from the density profile to the one-body direct correlation function is represented locally by a deep neural network. We substantiate the general framework for the hard sphere fluid and use grand canonical Monte Carlo simulation data of systems in randomized external environments during training and as reference. Functional calculus is implemented on the basis of the neural network to access higher-order correlation functions via automatic differentiation and the free energy via functional line integration. Thermal Noether sum rules are validated explicitly. We demonstrate the use of the neural functional in the self-consistent calculation of density profiles. The results outperform those from state-of-the-art fundamental measure density functional theory. The low cost of solving an associated Euler-Lagrange equation allows to bridge the gap from the system size of the original training data to macroscopic predictions upon maintaining near-simulation microscopic precision. These results establish the machine learning of functionals as an effective tool in the multiscale description of soft matter.
△ Less
Submitted 27 September, 2023; v1 submitted 10 July, 2023;
originally announced July 2023.
-
Practical Non-Invasive Probing Attacks Against Novel Carbon-Nanotube-Based Physical Unclonable Functions
Authors:
Nikolaos Athanasios Anagnostopoulos,
Alexander Braml,
Nico Mexis,
Florian Frank,
Simon Böttger,
Martin Hartmann,
Sascha Hermann,
Elif Bilge Kavun,
Stefan Katzenbeisser,
Tolga Arul
Abstract:
As the number of devices being interconnected increases, so does also the demand for (lightweight) security. To this end, Physical Unclonable Functions (PUFs) have been proposed as hardware primitives that can act as roots of trust and security. Recently, a new type of PUF based on Carbon NanoTubes (CNTs) has been proposed. At the same time, attacks and testing based on direct electrical probing a…
▽ More
As the number of devices being interconnected increases, so does also the demand for (lightweight) security. To this end, Physical Unclonable Functions (PUFs) have been proposed as hardware primitives that can act as roots of trust and security. Recently, a new type of PUF based on Carbon NanoTubes (CNTs) has been proposed. At the same time, attacks and testing based on direct electrical probing appear to be moving towards non-invasive techniques. In this context, this work attempts to examine the potential for practical non-invasive probing attacks against the CNT-PUF, a novel PUF based on CNTs. Our results indicate that direct probing might potentially compromise the security of this PUF. Nevertheless, we note that this holds true only in the case that the attacker can directly probe the wire corresponding to the secret value of each CNT-PUF cell. Thus, we can conclude that the examined CNT-PUFs are rather resilient to direct probing attacks, that non-invasive probing methods appear to be promising for testing such PUFs, and that, in order for the attacker to gain the full-length value of the secret, all the relevant channels would need to be probed. Nevertheless, as our work proves, practical non-invasive attacks against the CNT-PUF are feasible and adequate countermeasures need to be employed in order to address this issue.
△ Less
Submitted 3 July, 2023;
originally announced July 2023.
-
A sustainable infrastructure concept for improved accessibility, reusability, and archival of research software
Authors:
Timo Koch,
Dennis Gläser,
Anett Seeland,
Sarbani Roy,
Katharina Schulze,
Kilian Weishaupt,
David Boehringer,
Sibylle Hermann,
Bernd Flemisch
Abstract:
Research software is an integral part of most research today and it is widely accepted that research software artifacts should be accessible and reproducible. However, the sustainable archival of research software artifacts is an ongoing effort. We identify research software artifacts as snapshots of the current state of research and an integral part of a sustainable cycle of software development,…
▽ More
Research software is an integral part of most research today and it is widely accepted that research software artifacts should be accessible and reproducible. However, the sustainable archival of research software artifacts is an ongoing effort. We identify research software artifacts as snapshots of the current state of research and an integral part of a sustainable cycle of software development, research, and publication. We develop requirements and recommendations to improve the archival, access, and reuse of research software artifacts based on installable, configurable, extensible research software, and sustainable public open-access infrastructure. The described goal is to enable the reuse and exploration of research software beyond published research results, in parallel with reproducibility efforts, and in line with the FAIR principles for data and software. Research software artifacts can be reused in varying scenarios. To this end, we design a multi-modal representation concept supporting multiple reuse scenarios. We identify types of research software artifacts that can be viewed as different modes of the same software-based research result, for example, installation-free configurable browser-based apps to containerized environments, descriptions in journal publications and software documentation, or source code with installation instructions. We discuss how the sustainability and reuse of research software are enhanced or enabled by a suitable archive infrastructure. Finally, at the example of a pilot project at the University of Stuttgart, Germany -- a collaborative effort between research software developers and infrastructure providers -- we outline practical challenges and experiences
△ Less
Submitted 26 January, 2023;
originally announced January 2023.
-
Perspective: How to overcome dynamical density functional theory
Authors:
Daniel de las Heras,
Toni Zimmermann,
Florian Sammüller,
Sophie Hermann,
Matthias Schmidt
Abstract:
We argue in favour of developing a comprehensive dynamical theory for rationalizing, predicting, designing, and machine learning nonequilibrium phenomena that occur in soft matter. To give guidance for navigating the theoretical and practical challenges that lie ahead, we discuss and exemplify the limitations of dynamical density functional theory. Instead of the implied adiabatic sequence of equi…
▽ More
We argue in favour of developing a comprehensive dynamical theory for rationalizing, predicting, designing, and machine learning nonequilibrium phenomena that occur in soft matter. To give guidance for navigating the theoretical and practical challenges that lie ahead, we discuss and exemplify the limitations of dynamical density functional theory. Instead of the implied adiabatic sequence of equilibrium states that this approach provides as a makeshift for the true time evolution, we posit that the pending theoretical tasks lie in developing a systematic understanding of the dynamical functional relationships that govern the genuine nonequilibrium physics. While static density functional theory gives a comprehensive account of the equilibrium properties of many-body systems, we argue that power functional theory is the only present contender to shed similar insights into nonequilibrium dynamics, including the recognition and implementation of exact sum rules that result from the Noether theorem. As~a~demonstration of the power functional point of view, we consider an idealized steady sedimentation flow of the three-dimensional Lennard-Jones fluid and machine-learn the kinematic map from the mean motion to the internal force field. The trained model is capable of both predicting and designing the steady state dynamics universally for various target density modulations. This demonstrates the significant potential of using such techniques in nonequilibrium many-body physics and overcomes both the conceptual constraints of dynamical density functional theory as well as the limited availability of its analytical functional approximations.
△ Less
Submitted 23 March, 2023; v1 submitted 28 January, 2023;
originally announced January 2023.
-
Noether-Constrained Correlations in Equilibrium Liquids
Authors:
Florian Sammüller,
Sophie Hermann,
Daniel de las Heras,
Matthias Schmidt
Abstract:
Liquid structure carries deep imprints of an inherent thermal invariance against a spatial transformation of the underlying classical many-body Hamiltonian. At first order in the transformation field the Noether theorem yields the local force balance. Three distinct two-body correlation functions emerge at second order, namely the standard two-body density, the localized force-force correlation fu…
▽ More
Liquid structure carries deep imprints of an inherent thermal invariance against a spatial transformation of the underlying classical many-body Hamiltonian. At first order in the transformation field the Noether theorem yields the local force balance. Three distinct two-body correlation functions emerge at second order, namely the standard two-body density, the localized force-force correlation function, and the localized force gradient. An exact Noether sum rule interrelates these correlators. Simulations of Lennard-Jones, Yukawa, soft-sphere dipolar, Stockmayer, Gay-Berne and Weeks-Chandler-Andersen liquids, of monatomic water and of a colloidal gel former demonstrate their fundamental role in the characterization of spatial structure.
△ Less
Submitted 7 June, 2023; v1 submitted 26 January, 2023;
originally announced January 2023.
-
LuSEE 'Night': The Lunar Surface Electromagnetics Experiment
Authors:
Stuart D. Bale,
Neil Bassett,
Jack O. Burns,
Johnny Dorigo Jones,
Keith Goetz,
Christian Hellum-Bye,
Sven Hermann,
Joshua Hibbard,
Milan Maksimovic,
Ryan McLean,
Raul Monsalve,
Paul O'Connor,
Aaron Parsons,
Marc Pulupa,
Rugved Pund,
David Rapetti,
Kaja M. Rotermund,
Ben Saliwanchik,
Anze Slosar,
David Sundkvist,
Aritoki Suzuki
Abstract:
The Lunar Surface Electromagnetics Explorer 'LuSEE Night' is a low frequency radio astronomy experiment that will be delivered to the farside of the Moon by the NASA Commercial Lunar Payload Services (CLPS) program in late 2025 or early 2026. The payload system is being developed jointly by NASA and the US Department of Energy (DOE) and consists of a 4 channel, 50 MHz Nyquist baseband receiver sys…
▽ More
The Lunar Surface Electromagnetics Explorer 'LuSEE Night' is a low frequency radio astronomy experiment that will be delivered to the farside of the Moon by the NASA Commercial Lunar Payload Services (CLPS) program in late 2025 or early 2026. The payload system is being developed jointly by NASA and the US Department of Energy (DOE) and consists of a 4 channel, 50 MHz Nyquist baseband receiver system and 2 orthogonal $\sim$6m tip-to-tip electric dipole antennas. LuSEE Night will enjoy standalone operations through the lunar night, without the electromagnetic interference (EMI) of an operating lander system and antipodal to our noisy home planet.
△ Less
Submitted 24 January, 2023;
originally announced January 2023.
-
Comparative study of force-based classical density functional theory
Authors:
Florian Sammüller,
Sophie Hermann,
Matthias Schmidt
Abstract:
We reexamine results obtained with the recently proposed density functional theory framework based on forces (force-DFT) [Tschopp et al., Phys. Rev. E 106, 014115 (2022)]. We compare inhomogeneous density profiles for hard sphere fluids to results from both standard density functional theory and from computer simulations. Test situations include the equilibrium hard sphere fluid adsorbed against a…
▽ More
We reexamine results obtained with the recently proposed density functional theory framework based on forces (force-DFT) [Tschopp et al., Phys. Rev. E 106, 014115 (2022)]. We compare inhomogeneous density profiles for hard sphere fluids to results from both standard density functional theory and from computer simulations. Test situations include the equilibrium hard sphere fluid adsorbed against a planar hard wall and the dynamical relaxation of hard spheres in a switched harmonic potential. The comparison to grand canonical Monte Carlo simulation profiles shows that equilibrium force-DFT alone does not improve upon results obtained with the standard Rosenfeld functional. Similar behavior holds for the relaxation dynamics, where we use our event-driven Brownian dynamics data as benchmark. Based on an appropriate linear combination of standard and force-DFT results, we investigate a simple hybrid scheme which rectifies these deficiencies in both the equilibrium and the dynamical case. We explicitly demonstrate that although the hybrid method is based on the original Rosenfeld fundamental measure functional, its performance is comparable to that of the more advanced White Bear theory.
△ Less
Submitted 3 March, 2023; v1 submitted 4 December, 2022;
originally announced December 2022.
-
Reducing the background in X-ray imaging detectors via machine learning
Authors:
D. R. Wilkins,
S. W. Allen,
E. D. Miller,
M. Bautz,
T. Chattopadhyay,
R. Foster,
C. E. Grant,
S. Hermann,
R. Kraft,
R. G. Morris,
P. Nulsen,
G. Schellenberger
Abstract:
The sensitivity of astronomical X-ray detectors is limited by the instrumental background. The background is especially important when observing low surface brightness sources that are critical for many of the science cases targeted by future X-ray observatories, including Athena and future US-led flagship or probe-class X-ray missions. Above 2keV, the background is dominated by signals induced by…
▽ More
The sensitivity of astronomical X-ray detectors is limited by the instrumental background. The background is especially important when observing low surface brightness sources that are critical for many of the science cases targeted by future X-ray observatories, including Athena and future US-led flagship or probe-class X-ray missions. Above 2keV, the background is dominated by signals induced by cosmic rays interacting with the spacecraft and detector. We develop novel machine learning algorithms to identify events in next-generation X-ray imaging detectors and to predict the probability that an event is induced by a cosmic ray vs. an astrophysical X-ray photon, enabling enhanced filtering of the cosmic ray-induced background. We find that by learning the typical correlations between the secondary events that arise from a single primary, machine learning algorithms are able to successfully identify cosmic ray-induced background events that are missed by traditional filtering methods employed on current-generation X-ray missions, reducing the unrejected background by as much as 30 per cent.
△ Less
Submitted 16 August, 2022;
originally announced August 2022.
-
Force balance in thermal quantum many-body systems from Noether's theorem
Authors:
Sophie Hermann,
Matthias Schmidt
Abstract:
We address the consequences of invariance properties of the free energy of spatially inhomogeneous quantum many-body systems. We consider a specific position-dependent transformation of the system that consists of a spatial deformation and a corresponding locally resolved change of momenta. This operator transformation is canonical and hence equivalent to a unitary transformation on the underlying…
▽ More
We address the consequences of invariance properties of the free energy of spatially inhomogeneous quantum many-body systems. We consider a specific position-dependent transformation of the system that consists of a spatial deformation and a corresponding locally resolved change of momenta. This operator transformation is canonical and hence equivalent to a unitary transformation on the underlying Hilbert space of the system. As a consequence, the free energy is an invariant under the transformation. Noether's theorem for invariant variations then allows to derive an exact sum rule, which we show to be the locally resolved equilibrium one-body force balance. For the special case of homogeneous shifting, the sum rule states that the average global external force vanishes in thermal equilibrium.
△ Less
Submitted 17 October, 2022; v1 submitted 31 July, 2022;
originally announced August 2022.
-
Variance of fluctuations from Noether invariance
Authors:
Sophie Hermann,
Matthias Schmidt
Abstract:
The strength of fluctuations, as measured by their variance, is paramount in the quantitative description of a large class of physical systems, ranging from simple and complex liquids to active fluids and solids. Fluctuations originate from the irregular motion of thermal degrees of freedom and statistical mechanics facilitates their description. Here we demonstrate that fluctuations are constrain…
▽ More
The strength of fluctuations, as measured by their variance, is paramount in the quantitative description of a large class of physical systems, ranging from simple and complex liquids to active fluids and solids. Fluctuations originate from the irregular motion of thermal degrees of freedom and statistical mechanics facilitates their description. Here we demonstrate that fluctuations are constrained by the inherent symmetries of the given system. For particle-based classical many-body systems, Noether invariance at second order in the symmetry parameter leads to exact sum rules. These identities interrelate the global force variance with the mean potential energy curvature. Noether invariance is restored by an exact balance between these distinct mechanisms. The sum rules provide a practical guide for assessing and constructing theories, for ensuring self-consistency in simulation work, and for providing a systematic pathway to the theoretical quantification of fluctuations.
△ Less
Submitted 14 September, 2022; v1 submitted 29 March, 2022;
originally announced March 2022.
-
Force density functional theory in- and out-of-equilibrium
Authors:
S. M. Tschopp,
F. Sammüller,
S. Hermann,
M. Schmidt,
J. M. Brader
Abstract:
When a fluid is subject to an external field, as is the case near an interface or under spatial confinement, then the density becomes spatially inhomogeneous. Although the one-body density provides much useful information, a higher level of resolution is provided by the two-body correlations. These give a statistical description of the internal microstructure of the fluid and enable calculation of…
▽ More
When a fluid is subject to an external field, as is the case near an interface or under spatial confinement, then the density becomes spatially inhomogeneous. Although the one-body density provides much useful information, a higher level of resolution is provided by the two-body correlations. These give a statistical description of the internal microstructure of the fluid and enable calculation of the average interparticle force, which plays an essential role in determining both the equilibrium and dynamic properties of interacting fluids. We present a theoretical framework for the description of inhomogeneous (classical) many-body systems, based explicitly on the two-body correlation functions. By consideration of local Noether-invariance against spatial distortion of the system we demonstrate the fundamental status of the Yvon-Born-Green (YBG) equation as a local force-balance within the fluid. Using the inhomogeneous Ornstein-Zernike equation we show that the two-body correlations are density functionals and, thus, that the average interparticle force entering the YBG equation is also a functional of the one-body density. The force-based theory we develop provides an alternative to standard density functional theory for the study of inhomogeneous systems both in- and out-of-equilibrium. We compare force-based density profiles to the results of the standard potential-based (dynamical) density functional theory. In equilibrium, we confirm both analytically and numerically that the standard approach yields profiles that are consistent with the compressibility pressure, whereas the force-density functional gives profiles consistent with the virial pressure. The structure of the theory offers deep insights into the nature of correlation in dense and inhomogeneous systems.
△ Less
Submitted 13 July, 2022; v1 submitted 3 March, 2022;
originally announced March 2022.
-
Why Noether's Theorem applies to Statistical Mechanics
Authors:
Sophie Hermann,
Matthias Schmidt
Abstract:
Noether's Theorem is familiar to most physicists due its fundamental role in linking the existence of conservation laws to the underlying symmetries of a physical system. Typically the systems are described in the particle-based context of classical mechanics or on the basis of field theory. We have recently shown [Commun. Phys. $\textbf{4}$, 176 (2021)] that Noether's reasoning also applies to th…
▽ More
Noether's Theorem is familiar to most physicists due its fundamental role in linking the existence of conservation laws to the underlying symmetries of a physical system. Typically the systems are described in the particle-based context of classical mechanics or on the basis of field theory. We have recently shown [Commun. Phys. $\textbf{4}$, 176 (2021)] that Noether's reasoning also applies to thermal systems, where fluctuations are paramount and one aims for a statistical mechanical description. Here we give a pedagogical introduction based on the canonical ensemble and apply it explicitly to ideal sedimentation. The relevant mathematical objects, such as the free energy, are viewed as functionals. This vantage point allows for systematic functional differentiation and the resulting identities express properties of both macroscopic average forces and molecularly resolved correlations in many-body systems, both in and out-of-equilibrium, and for active Brownian particles. To provide further background, we briefly describe the variational principles of classical density functional theory, of power functional theory, and of classical mechanics.
△ Less
Submitted 2 March, 2022; v1 submitted 21 September, 2021;
originally announced September 2021.
-
Noether's Theorem in Statistical Mechanics
Authors:
Sophie Hermann,
Matthias Schmidt
Abstract:
Noether's calculus of invariant variations yields exact identities from functional symmetries. The standard application to an action integral allows to identify conservation laws. Here we rather consider generating functionals, such as the free energy and the power functional, for equilibrium and driven many-body systems. Translational and rotational symmetry operations yield mechanical laws. Thes…
▽ More
Noether's calculus of invariant variations yields exact identities from functional symmetries. The standard application to an action integral allows to identify conservation laws. Here we rather consider generating functionals, such as the free energy and the power functional, for equilibrium and driven many-body systems. Translational and rotational symmetry operations yield mechanical laws. These global identities express vanishing of total internal and total external forces and torques. We show that functional differentiation then leads to hierarchies of local sum rules that interrelate density correlators as well as static and time direct correlation functions, including memory. For anisotropic particles, orbital and spin motion become systematically coupled. The theory allows us to shed new light on the spatio-temporal coupling of correlations in complex systems. As applications we consider active Brownian particles, where the theory clarifies the role of interfacial forces in motility-induced phase separation. For active sedimentation, the center-of-mass motion is constrained by an internal Noether sum rule.
△ Less
Submitted 13 August, 2021; v1 submitted 27 May, 2021;
originally announced May 2021.
-
Phase separation of active Brownian particles in two dimensions: Anything for a quiet life
Authors:
Sophie Hermann,
Daniel de las Heras,
Matthias Schmidt
Abstract:
Active Brownian particles display self-propelled movement, which can be modelled as arising from a one-body force. Although their interparticle interactions are purely repulsive, for strong self propulsion the swimmers phase separate into dilute and dense phases. We describe in detail a recent theory (Phys. Rev. E 100, 052604 (2019); Phys. Rev. Lett. 128, 26802 (2019)) for such motility induced ph…
▽ More
Active Brownian particles display self-propelled movement, which can be modelled as arising from a one-body force. Although their interparticle interactions are purely repulsive, for strong self propulsion the swimmers phase separate into dilute and dense phases. We describe in detail a recent theory (Phys. Rev. E 100, 052604 (2019); Phys. Rev. Lett. 128, 26802 (2019)) for such motility induced phase-separation. Starting from the continuity equation and the force density balance, the description is based on four superadiabatic contributions to the internal force density. Here the superadiabatic forces are due to the flow in the system and they act on top of the adiabatic forces that arise from the equilibrium free energy. Phase coexistence is described by bulk state functions and agrees quantitatively with Brownian dynamics simulation results from the literature. We describe in detail all analytical steps to fully resolve the spatial and orientational dependence of the one-body density and current. The decomposition into angular Fourier series leads to coupling of total density, polarization and all higher modes. We describe the power functional approach, including the kinematic dependence of the superadiabatic force fields and the quiet life effect that pushes particles from fast to slow regions, and hence induces the phase separation.
△ Less
Submitted 24 March, 2021; v1 submitted 5 March, 2021;
originally announced March 2021.
-
An Environment for Sustainable Research Software in Germany and Beyond: Current State, Open Challenges, and Call for Action
Authors:
Hartwig Anzt,
Felix Bach,
Stephan Druskat,
Frank Löffler,
Axel Loewe,
Bernhard Y. Renard,
Gunnar Seemann,
Alexander Struck,
Elke Achhammer,
Piush Aggarwal,
Franziska Appel,
Michael Bader,
Lutz Brusch,
Christian Busse,
Gerasimos Chourdakis,
Piotr W. Dabrowski,
Peter Ebert,
Bernd Flemisch,
Sven Friedl,
Bernadette Fritzsch,
Maximilian D. Funk,
Volker Gast,
Florian Goth,
Jean-Noël Grad,
Sibylle Hermann
, et al. (18 additional authors not shown)
Abstract:
Research software has become a central asset in academic research. It optimizes existing and enables new research methods, implements and embeds research knowledge, and constitutes an essential research product in itself. Research software must be sustainable in order to understand, replicate, reproduce, and build upon existing research or conduct new research effectively. In other words, software…
▽ More
Research software has become a central asset in academic research. It optimizes existing and enables new research methods, implements and embeds research knowledge, and constitutes an essential research product in itself. Research software must be sustainable in order to understand, replicate, reproduce, and build upon existing research or conduct new research effectively. In other words, software must be available, discoverable, usable, and adaptable to new needs, both now and in the future. Research software therefore requires an environment that supports sustainability. Hence, a change is needed in the way research software development and maintenance are currently motivated, incentivized, funded, structurally and infrastructurally supported, and legally treated. Failing to do so will threaten the quality and validity of research. In this paper, we identify challenges for research software sustainability in Germany and beyond, in terms of motivation, selection, research software engineering personnel, funding, infrastructure, and legal aspects. Besides researchers, we specifically address political and academic decision-makers to increase awareness of the importance and needs of sustainable research software practices. In particular, we recommend strategies and measures to create an environment for sustainable research software, with the ultimate goal to ensure that software-driven research is valid, reproducible and sustainable, and that software is recognized as a first class citizen in research. This paper is the outcome of two workshops run in Germany in 2019, at deRSE19 - the first International Conference of Research Software Engineers in Germany - and a dedicated DFG-supported follow-up workshop in Berlin.
△ Less
Submitted 5 May, 2020; v1 submitted 27 April, 2020;
originally announced May 2020.
-
Active interface polarization is a state function
Authors:
Sophie Hermann,
Matthias Schmidt
Abstract:
We prove three exact sum rules that relate the polarization of active Brownian particles to their one-body current: (i) The total polarization vanishes, provided that there is no net flux through the boundaries, (ii) at any planar wall the polarization is determined by the magnitude of the bulk current, and (iii) the total interface polarization between phase-separated fluid states is rigorously d…
▽ More
We prove three exact sum rules that relate the polarization of active Brownian particles to their one-body current: (i) The total polarization vanishes, provided that there is no net flux through the boundaries, (ii) at any planar wall the polarization is determined by the magnitude of the bulk current, and (iii) the total interface polarization between phase-separated fluid states is rigorously determined by the gas-liquid current difference. This result precludes the influence of the total interface polarization on active bulk coexistence and questions the proposed coupling of interface to bulk.
△ Less
Submitted 17 March, 2020;
originally announced March 2020.
-
Non-negative Interfacial Tension in Phase-Separated Active Brownian Particles
Authors:
Sophie Hermann,
Daniel de las Heras,
Matthias Schmidt
Abstract:
We present a microscopic theory for the nonequilibrium interfacial tension $γ_{\rm gl}$ of the free interface between gas and liquid phases of active Brownian particles. The underlying square gradient treatment and the splitting of the force balance in flow and structural contributions is general and applies to inhomogeneous nonequilibrium steady states. We find $γ_{\rm gl}\geq 0$, which opposes c…
▽ More
We present a microscopic theory for the nonequilibrium interfacial tension $γ_{\rm gl}$ of the free interface between gas and liquid phases of active Brownian particles. The underlying square gradient treatment and the splitting of the force balance in flow and structural contributions is general and applies to inhomogeneous nonequilibrium steady states. We find $γ_{\rm gl}\geq 0$, which opposes claims by Bialké et al. [Phys. Rev. Lett. $\bf{115}$, 098301 (2015)] and delivers the theoretical justification for the widely observed interfacial stability in active Brownian dynamics many-body simulations.
△ Less
Submitted 28 January, 2020; v1 submitted 9 December, 2019;
originally announced December 2019.
-
Non-reciprocal hidden degrees of freedom: A unifying perspective on memory, feedback, and activity
Authors:
Sarah A. M. Loos,
Simon M. Hermann,
Sabine H. L. Klapp
Abstract:
We show that memory, feedback, and activity are all describable by the same unifying concept, that is non-reciprocal (NR) coupling. We demonstrate that characteristic thermodynamic features of these intrinsically nonequilibrium systems are reproduced by low-dimensional Markovian networks with NR coupling, which we establish as minimal models for such complex systems. NR coupling alone implies a vi…
▽ More
We show that memory, feedback, and activity are all describable by the same unifying concept, that is non-reciprocal (NR) coupling. We demonstrate that characteristic thermodynamic features of these intrinsically nonequilibrium systems are reproduced by low-dimensional Markovian networks with NR coupling, which we establish as minimal models for such complex systems. NR coupling alone implies a violation of the fluctuation-dissipation relation, which is inevitably connected to entropy production, i.e., irreversibility. Hiding the NR coupled degrees of freedom renders non-Markovian, one-variable Langevin descriptions with complex types of memory, for which we find a generalized second law involving information flow. We demonstrate that non-monotonic memory is inextricably linked to NR coupling. Furthermore, we discuss discrete time delay as the infinite-dimensional limit, and find a divergent entropy production, corresponding to unbounded cost for precisely storing a Brownian trajectory.
△ Less
Submitted 18 October, 2019;
originally announced October 2019.
-
Phase coexistence of active Brownian particles
Authors:
Sophie Hermann,
Philip Krinninger,
Daniel de las Heras,
Matthias Schmidt
Abstract:
We investigate motility-induced phase separation of active Brownian particles, which are modeled as purely repulsive spheres that move due to a constant swim force with freely diffusing orientation. We develop on the basis of power functional concepts an analytical theory for nonequilibrium phase coexistence and interfacial structure. Theoretical predictions are validated against Brownian dynamics…
▽ More
We investigate motility-induced phase separation of active Brownian particles, which are modeled as purely repulsive spheres that move due to a constant swim force with freely diffusing orientation. We develop on the basis of power functional concepts an analytical theory for nonequilibrium phase coexistence and interfacial structure. Theoretical predictions are validated against Brownian dynamics computer simulations. We show that the internal one-body force field has four nonequilibrium contributions: (i) isotropic drag and (ii) interfacial drag forces against the forward motion, (iii) a superadiabatic spherical pressure gradient and (iv) the quiet life gradient force. The intrinsic spherical pressure is balanced by the swim pressure, which arises from the polarization of the free interface. The quiet life force opposes the adiabatic force, which is due to the inhomogeneous density distribution. The balance of quiet life and adiabatic forces determines bulk coexistence via equality of two bulk state functions, which are independent of interfacial contributions. The internal force fields are kinematic functionals which depend on density and current, but are independent of external and swim forces, consistent with power functional theory. The phase transition originates from nonequilibrium repulsion, with the agile gas being more repulsive than the quiet liquid.
△ Less
Submitted 11 November, 2019; v1 submitted 14 October, 2019;
originally announced October 2019.
-
Active ideal sedimentation: Exact two-dimensional steady states
Authors:
Sophie Hermann,
Matthias Schmidt
Abstract:
We consider an ideal gas of active Brownian particles that undergo self-propelled motion and both translational and rotational diffusion under the influence of gravity. We solve analytically the corresponding Smoluchowski equation in two space dimensions for steady states. The resulting one-body density is given as a series, where each term is a product of an orientation-dependent Mathieu function…
▽ More
We consider an ideal gas of active Brownian particles that undergo self-propelled motion and both translational and rotational diffusion under the influence of gravity. We solve analytically the corresponding Smoluchowski equation in two space dimensions for steady states. The resulting one-body density is given as a series, where each term is a product of an orientation-dependent Mathieu function and a height-dependent exponential. A lower hard wall is implemented as a no-flux boundary condition. Numerical evaluation of the suitably truncated analytical solution shows the formation of two different spatial regimes upon increasing Peclet number. These regimes differ in their mean particle orientation and in their variation of the orientation-averaged density with height.
△ Less
Submitted 22 December, 2017;
originally announced December 2017.
-
The large core limit of spiral waves in excitable media: A numerical approach
Authors:
Sebastian Hermann,
Georg A. Gottwald
Abstract:
We modify the freezing method introduced by Beyn & Thuemmler, 2004, for analyzing rigidly rotating spiral waves in excitable media. The proposed method is designed to stably determine the rotation frequency and the core radius of rotating spirals, as well as the approximate shape of spiral waves in unbounded domains. In particular, we introduce spiral wave boundary conditions based on geometric ap…
▽ More
We modify the freezing method introduced by Beyn & Thuemmler, 2004, for analyzing rigidly rotating spiral waves in excitable media. The proposed method is designed to stably determine the rotation frequency and the core radius of rotating spirals, as well as the approximate shape of spiral waves in unbounded domains. In particular, we introduce spiral wave boundary conditions based on geometric approximations of spiral wave solutions by Archimedean spirals and by involutes of circles. We further propose a simple implementation of boundary conditions for the case when the inhibitor is non-diffusive, a case which had previously caused spurious oscillations.
We then utilize the method to numerically analyze the large core limit. The proposed method allows us to investigate the case close to criticality where spiral waves acquire infinite core radius and zero rotation frequency, before they begin to develop into retracting fingers. We confirm the linear scaling regime of a drift bifurcation for the rotation frequency and the core radius of spiral wave solutions close to criticality. This regime is unattainable with conventional numerical methods.
△ Less
Submitted 30 March, 2010;
originally announced March 2010.