Showing 1–50 of 57 results for author: Hartmann, C

Quantum Annealing based Power Grid Partitioning for Parallel Simulation

Authors: Carsten Hartmann, Junjie Zhang, Carlos D. Gonzalez Calaza, Thiemo Pesch, Kristel Michielsen, Andrea Benigni

Abstract: Graph partitioning has many applications in powersystems from decentralized state estimation to parallel simulation. Focusing on parallel simulation, optimal grid partitioning minimizes the idle time caused by different simulation times for the sub-networks and their components and reduces the overhead required to simulate the cuts. Partitioning a graph into two parts such that, for example, the c… ▽ More Graph partitioning has many applications in powersystems from decentralized state estimation to parallel simulation. Focusing on parallel simulation, optimal grid partitioning minimizes the idle time caused by different simulation times for the sub-networks and their components and reduces the overhead required to simulate the cuts. Partitioning a graph into two parts such that, for example, the cut is minimal and the subgraphs have equal size is an NP-hard problem. In this paper we show how optimal partitioning of a graph can be obtained using quantum annealing (QA). We show how to map the requirements for optimal splitting to a quadratic unconstrained binary optimization (QUBO) formulation and test the proposed formulation using a current D-Wave QPU. We show that the necessity to find an embedding of the QUBO on current D-Wave QPUs limits the problem size to under 200 buses and notably affects the time-to-solution. We finally discuss the implications on near-term implementation of QA in combination to traditional CPU or GPU based simulation. △ Less

Submitted 7 August, 2024; originally announced August 2024.

Comments: 10 pages, 5 figures

Computationally feasible bounds for the free energy of nonequilibrium steady states, applied to simple models of heat conduction

Authors: Luigi Delle Site, Carsten Hartmann

Abstract: In this paper we study computationally feasible bounds for relative free energies between two many-particle systems. Specifically, we consider systems out of equilibrium that do not necessarily satisfy a fluctuation-dissipation relation, but that nevertheless admit a nonequilibrium steady state that is reached asymptotically in the long-time limit. The bounds that we suggest are based on the well-… ▽ More In this paper we study computationally feasible bounds for relative free energies between two many-particle systems. Specifically, we consider systems out of equilibrium that do not necessarily satisfy a fluctuation-dissipation relation, but that nevertheless admit a nonequilibrium steady state that is reached asymptotically in the long-time limit. The bounds that we suggest are based on the well-known Bogoliubov inequality and variants of Gibbs' and Donsker-Varadhan variational principles. As a general paradigm, we consider systems of oscillators coupled to heat baths at different temperatures. For such systems, we define the free energy of the system relative to any given reference system (that may or may not be in thermal equilibrium) in terms of the Kullback-Leibler divergence between steady states. By employing a two-sided Bogoliubov inequality and a mean-variance approximation of the free energy (or cumulant generating function, we can efficiently estimate the free energy cost needed in passing from the reference system to the system out of equilibrium (characterised by a temperature gradient). A specific test case to validate our bounds are harmonic oscillator chains with ends that are coupled to Langevin thermostats at different temperatures; such a system is simple enough to allow for analytic calculations and general enough to be used as a prototype to estimate, e.g., heat fluxes or interface effects in a larger class of nonequilibrium particle systems. △ Less

Submitted 16 March, 2024; originally announced March 2024.

Comments: 14 pages, 7 figures

MSC Class: 82C31 (Primary); 80A19; 37M25 (Secondary)

Synchronized states of power grids and oscillator networks by convex optimization

Authors: Carsten Hartmann, Philipp C. Böttcher, David Gross, Dirk Witthaut

Abstract: Synchronization is essential for the operation of AC power systems: All generators in the power grid must rotate with fixed relative phases to enable a steady flow of electric power. Understanding the conditions for and the limitations of synchronization is of utmost practical importance. In this article, we propose a novel approach to compute and analyze the stable stationary states of a power gr… ▽ More Synchronization is essential for the operation of AC power systems: All generators in the power grid must rotate with fixed relative phases to enable a steady flow of electric power. Understanding the conditions for and the limitations of synchronization is of utmost practical importance. In this article, we propose a novel approach to compute and analyze the stable stationary states of a power grid or an oscillator network in terms of a convex optimization problem. This approach allows to systematically compute \emph{all} stable states where the phase difference across an edge does not exceed $π/2$.Furthermore, the optimization formulation allows to rigorously establish certain properties of synchronized states and to bound the error in the widely used linear power flow approximation. △ Less

Submitted 14 March, 2024; originally announced March 2024.

Comments: 19 pages, 7 figures

Risk-neutral limit of adaptive importance sampling of random stopping times

Authors: Carsten Hartmann, Annika Jöster

Abstract: We discuss importance sampling of exit problems that involve unbounded stopping times; examples are mean first passage times, transition rates or committor probabilities in molecular dynamics. The naive application of variance minimization techniques can lead to pathologies here, including proposal measures that are not absolutely continuous to the reference measure or importance sampling estimato… ▽ More We discuss importance sampling of exit problems that involve unbounded stopping times; examples are mean first passage times, transition rates or committor probabilities in molecular dynamics. The naive application of variance minimization techniques can lead to pathologies here, including proposal measures that are not absolutely continuous to the reference measure or importance sampling estimators that formally have zero variance, but that produce infinitely long trajectories. We illustrate these issues with simple examples and discuss a possible solution that is based on a risk-sensitive optimal control framework of importance sampling. △ Less

Submitted 13 February, 2024; originally announced February 2024.

Comments: 7 pages, 2 figures

MSC Class: 93E20 (primary); 65C05; 65C20 (secondary)

Taking the Next Step with Generative Artificial Intelligence: The Transformative Role of Multimodal Large Language Models in Science Education

Authors: Arne Bewersdorff, Christian Hartmann, Marie Hornberger, Kathrin Seßler, Maria Bannert, Enkelejda Kasneci, Gjergji Kasneci, Xiaoming Zhai, Claudia Nerdel

Abstract: The integration of Artificial Intelligence (AI), particularly Large Language Model (LLM)-based systems, in education has shown promise in enhancing teaching and learning experiences. However, the advent of Multimodal Large Language Models (MLLMs) like GPT-4 with vision (GPT-4V), capable of processing multimodal data including text, sound, and visual inputs, opens a new era of enriched, personalize… ▽ More The integration of Artificial Intelligence (AI), particularly Large Language Model (LLM)-based systems, in education has shown promise in enhancing teaching and learning experiences. However, the advent of Multimodal Large Language Models (MLLMs) like GPT-4 with vision (GPT-4V), capable of processing multimodal data including text, sound, and visual inputs, opens a new era of enriched, personalized, and interactive learning landscapes in education. Grounded in theory of multimedia learning, this paper explores the transformative role of MLLMs in central aspects of science education by presenting exemplary innovative learning scenarios. Possible applications for MLLMs could range from content creation to tailored support for learning, fostering competencies in scientific practices, and providing assessment and feedback. These scenarios are not limited to text-based and uni-modal formats but can be multimodal, increasing thus personalization, accessibility, and potential learning effectiveness. Besides many opportunities, challenges such as data protection and ethical considerations become more salient, calling for robust frameworks to ensure responsible integration. This paper underscores the necessity for a balanced approach in implementing MLLMs, where the technology complements rather than supplants the educator's role, ensuring thus an effective and ethical use of AI in science education. It calls for further research to explore the nuanced implications of MLLMs on the evolving role of educators and to extend the discourse beyond science education to other disciplines. Through the exploration of potentials, challenges, and future implications, we aim to contribute to a preliminary understanding of the transformative trajectory of MLLMs in science education and beyond. △ Less

Submitted 3 September, 2024; v1 submitted 1 January, 2024; originally announced January 2024.

Comments: revised version 2. September 2024

Scaling law for the size dependence of a finite-range quantum gas

Authors: Luigi Delle Site, Carsten Hartmann

Abstract: In a recent work [Reible et al., Phys. Rev. Res. 5, 023156, 2023], it has been shown that the mean particle-particle interaction across an ideal surface that divides a system into two parts, can be employed to estimate the size dependence for the thermodynamic accuracy of the system. In this work we propose its application to systems with finite range interactions that models a dense quantum gases… ▽ More In a recent work [Reible et al., Phys. Rev. Res. 5, 023156, 2023], it has been shown that the mean particle-particle interaction across an ideal surface that divides a system into two parts, can be employed to estimate the size dependence for the thermodynamic accuracy of the system. In this work we propose its application to systems with finite range interactions that models a dense quantum gases and derive an approximate size-dependence scaling law. In addition, we show that the application of the criterion is equivalent to the determination of a free energy response to a perturbation. The latter result confirms the complementarity of the criterion to other estimates of finite-size effects based on direct simulations and empirical structure or energy convergence criteria. △ Less

Submitted 24 January, 2024; v1 submitted 30 August, 2023; originally announced August 2023.

Comments: 14 pages, 1 figure

MSC Class: 82B10; 82D03

Transgressing the boundaries: towards a rigorous understanding of deep learning and its (non-)robustness

Authors: Carsten Hartmann, Lorenz Richter

Abstract: The recent advances in machine learning in various fields of applications can be largely attributed to the rise of deep learning (DL) methods and architectures. Despite being a key technology behind autonomous cars, image processing, speech recognition, etc., a notorious problem remains the lack of theoretical understanding of DL and related interpretability and (adversarial) robustness issues. Un… ▽ More The recent advances in machine learning in various fields of applications can be largely attributed to the rise of deep learning (DL) methods and architectures. Despite being a key technology behind autonomous cars, image processing, speech recognition, etc., a notorious problem remains the lack of theoretical understanding of DL and related interpretability and (adversarial) robustness issues. Understanding the specifics of DL, as compared to, say, other forms of nonlinear regression methods or statistical learning, is interesting from a mathematical perspective, but at the same time it is of crucial importance in practice: treating neural networks as mere black boxes might be sufficient in certain cases, but many applications require waterproof performance guarantees and a deeper understanding of what could go wrong and why it could go wrong. It is probably fair to say that, despite being mathematically well founded as a method to approximate complicated functions, DL is mostly still more like modern alchemy that is firmly in the hands of engineers and computer scientists. Nevertheless, it is evident that certain specifics of DL that could explain its success in applications demands systematic mathematical approaches. In this work, we review robustness issues of DL and particularly bridge concerns and attempts from approximation theory to statistical learning theory. Further, we review Bayesian Deep Learning as a means for uncertainty quantification and rigorous explainability. △ Less

Submitted 5 July, 2023; originally announced July 2023.

Magneto-active composites with locally tailored stiffness produced by laser powder bed fusion

Authors: Kilian Schäfer, Matthias Lutzi, Muhammad Bilal Khan, Lukas Schäfer, Konstantin Skokov, Imants Dirba, Sebastian Bruns, Iman Valizadeh, Oliver Weeger, Claas Hartmann, Mario Kupnik, Esmaeil Adabifiroozjaei, Leopoldo Molina-Luna, Oliver Gutfleisch

Abstract: Additive manufacturing technologies enable the production of complex and bioinspired shapes using magneto-responsive materials, which find diverse applications in soft robotics. Particularly, the development of composites with controlled gradients in mechanical properties offers new prospects for advancements in magneto-active materials. However, achieving such composites with gradients typically… ▽ More Additive manufacturing technologies enable the production of complex and bioinspired shapes using magneto-responsive materials, which find diverse applications in soft robotics. Particularly, the development of composites with controlled gradients in mechanical properties offers new prospects for advancements in magneto-active materials. However, achieving such composites with gradients typically involves complex multi-material printing procedures. In this study, a single-step laser powder bed fusion (LPBF) process is proposed that enables precise local adjustments of the mechanical stiffness within magneto-active composites. By utilizing distinct laser parameters in specific regions of a composite containing thermoplastic polyurethane and atomized magnetic powder derived from hard magnetic Nd-Fe-B, the stiffness of the composite can be modified within the range of 2 to 22 MPa. Various magneto-responsive actuators with locally tailored stiffness are fabricated and their magnetic performance is investigated. The enhanced response exhibited by actuators with locally adjusted mechanical properties in comparison to their homogeneous counterparts with identical geometries is shown. As a demonstration of a biomedical application, a magnetically responsive stent with localized adjustment is presented with the ability to meet specific requirements in terms of geometry and local stiffness based on an individual's anatomy and disease condition. The proposed method presents an approach for creating functionally graded materials using LPBF, not only for magneto-active materials but also for several other structural and functional materials. △ Less

Submitted 20 December, 2023; v1 submitted 4 May, 2023; originally announced May 2023.

Journal ref: Additive Manufacturing 79 (2024) 103905

Acousto-Optic Modulation in Ambient Air

Authors: Yannick Schrödel, Claas Hartmann, Tino Lang, Jiaan Zheng, Max Steudel, Matthias Rutsch, Sarper H. Salman, Martin Kellert, Mikhail Pergament, Thomas Hahn-Jose, Sven Suppelt, Jan Helge Dörsam, Anne Harth, Wim P. Leemans, Franz X. Kärtner, Ingmar Hartl, Mario Kupnik, Christoph M. Heyl

Abstract: Control over intensity, shape, direction, and phase of coherent light is essential in numerous fields, reaching from gravitational wave astronomy over quantum metrology and ultrafast sciences to semi-conductor fabrication. Modern laser optics, however, frequently demands parameter regimes where either the wavelength or the optical power restricts control due to linear absorption, light-induced dam… ▽ More Control over intensity, shape, direction, and phase of coherent light is essential in numerous fields, reaching from gravitational wave astronomy over quantum metrology and ultrafast sciences to semi-conductor fabrication. Modern laser optics, however, frequently demands parameter regimes where either the wavelength or the optical power restricts control due to linear absorption, light-induced damage or optical nonlinearity. The properties of solid media, upon which most photonic control schemes rely, impose these limitations. We propose to circumvent these constraints using gaseous media tailored by high-intensity ultrasound waves. We demonstrate a first implementation of this approach by deflecting ultrashort laser pulses using ultrasound waves in ambient air, entirely omitting transmissive solid media. At optical peak powers of 20 GW exceeding previous limits of solid-based acousto-optic modulation by about three orders of magnitude, we reach a deflection efficiency greater than 50% while preserving excellent beam quality. Our approach is not limited to laser pulse deflection via acousto-optic modulation: gas-phase photonic schemes controlled by sonic waves can prospectively be translated to various optical methods, e.g., lenses or waveguides, rendering them effectively invulnerable against damage and opening up new spectral regions. △ Less

Submitted 28 April, 2023; v1 submitted 13 April, 2023; originally announced April 2023.

Comments: 23 pages, including 9 pages of main text and 9 pages of Methods and Extended Data, 3 figures, 3 Extended Data figures, 1 Extended Data table

Report number: PUBDB-2023-01719

Long Short-Term Memory Neural Network for Temperature Prediction in Laser Powder Bed Additive Manufacturing

Authors: Ashkan Mansouri Yarahmadi, Michael Breuß, Carsten Hartmann

Abstract: In context of laser powder bed fusion (L-PBF), it is known that the properties of the final fabricated product highly depend on the temperature distribution and its gradient over the manufacturing plate. In this paper, we propose a novel means to predict the temperature gradient distributions during the printing process by making use of neural networks. This is realized by employing heat maps prod… ▽ More In context of laser powder bed fusion (L-PBF), it is known that the properties of the final fabricated product highly depend on the temperature distribution and its gradient over the manufacturing plate. In this paper, we propose a novel means to predict the temperature gradient distributions during the printing process by making use of neural networks. This is realized by employing heat maps produced by an optimized printing protocol simulation and used for training a specifically tailored recurrent neural network in terms of a long short-term memory architecture. The aim of this is to avoid extreme and inhomogeneous temperature distribution that may occur across the plate in the course of the printing process. In order to train the neural network, we adopt a well-engineered simulation and unsupervised learning framework. To maintain a minimized average thermal gradient across the plate, a cost function is introduced as the core criteria, which is inspired and optimized by considering the well-known traveling salesman problem (TSP). As time evolves the unsupervised printing process governed by TSP produces a history of temperature heat maps that maintain minimized average thermal gradient. All in one, we propose an intelligent printing tool that provides control over the substantial printing process components for L-PBF, i.e.\ optimal nozzle trajectory deployment as well as online temperature prediction for controlling printing quality. △ Less

Submitted 30 January, 2023; originally announced January 2023.

Combining advanced photoelectron spectroscopy approaches to analyse deeply buried GaP(As)/Si(100) interfaces: Interfacial chemical states and complete band energy diagrams

Authors: O. Romanyuk, A. Paszuk, I. Gordeev, R. G. Wilks, S. Ueda, C. Hartmann, R. Félix, M. Bär, C. Schlueter, A. Gloskovskii, I. Bartoš, M. Nandy, J. Houdková, P. Jiříček, W. Jaegermann, J. P. Hofmann, T. Hannappel

Abstract: The epitaxial growth of the polar GaP(100) on the nonpolar Si(100) substrate suffers from inevitable defects at the antiphase domain boundaries, resulting from mono-atomic steps on the Si(100) surface. Stabilization of Si(100) substrate surfaces with arsenic is a promising technological step enabling the preparation of Si substrates with double atomic steps and reduced density of the APDs. In this… ▽ More The epitaxial growth of the polar GaP(100) on the nonpolar Si(100) substrate suffers from inevitable defects at the antiphase domain boundaries, resulting from mono-atomic steps on the Si(100) surface. Stabilization of Si(100) substrate surfaces with arsenic is a promising technological step enabling the preparation of Si substrates with double atomic steps and reduced density of the APDs. In this paper, 4-50 nm thick GaP epitaxial films were grown on As-terminated Si(100) substrates with different types of doping, miscuts, and As-surface termination by metalorganic vapor phase epitaxy. The GaP(As)/Si(100) heterostructures were investigated by X-ray photoelectron spectroscopy (XPS) combined with gas cluster ion beam (GCIB) sputtering and by hard X-ray photoelectron spectroscopy (HAXPES). We found residuals of arsenic atoms in the GaP lattice (0.2-0.3 at.%) and a localization of As atoms at the GaP(As)/Si(100) interface (1 at.%). Deconvolution of core level peaks revealed interface core level shifts. In As core levels, chemical shifts between 0.5-0.8 eV were measured and identified by angle-resolved XPS measurements. Similar valence band offset (VBO) values of 0.6 eV were obtained, regardless of the doping type of Si substrate, Si substrate miscut or type of As-terminated Si substrate surface. The band alignment diagram of the heterostructure was deduced. △ Less

Submitted 17 June, 2022; originally announced June 2022.

Two-sided Bogoliubov inequality to estimate finite size effects in quantum molecular simulations

Authors: Benedikt Reible, Carsten Hartmann, Luigi Delle Site

Abstract: We generalise the two-sided Bogoliubov inequality for classical particles from [L. Delle Site et al., J.Stat.Mech.Th.Exp. 083201 (2017)] to systems of quantum particles. As in the classical set-up, the inequality leads to upper and lower bounds for the free energy difference associated with the partitioning of a large system into smaller, independent subsystems. From a thermodynamic modelling poin… ▽ More We generalise the two-sided Bogoliubov inequality for classical particles from [L. Delle Site et al., J.Stat.Mech.Th.Exp. 083201 (2017)] to systems of quantum particles. As in the classical set-up, the inequality leads to upper and lower bounds for the free energy difference associated with the partitioning of a large system into smaller, independent subsystems. From a thermodynamic modelling point of view, the free energy difference determines the finite size correction needed to consistently treat a small system as a representation of a large system. Applications of the bounds to quantify finite size effects are ubiquitous in physics, chemistry, material science, or biology, to name just a few; in particular it is relevant for molecular dynamics simulations in which a small portion of a system is usually taken as representative of the idealized large system. △ Less

Submitted 5 September, 2022; v1 submitted 19 November, 2021; originally announced November 2021.

Journal ref: Letters in Mathematical Physics, 112, 97 (2022)

Origin of interface limitation in CuInS$_2$ based solar cells

Authors: Mohit Sood, Jakob Bombsch, Alberto Lomuscio, Sudhanshu Shukla, Alberto Lomuscio, Claudia Hartmann, Johannes Frisch, Wolfgang Bremsteller, Shigenori Ueda, Regan G. Wilks, Marcus Bär, Susanne Siebentritt

Abstract: Copper indium disulfide (CuInS$_2$) grown under Cu-rich conditions exhibits high optical quality but suffers predominantly from charge carrier interface recombination resulting in poor solar cell performance. An unfavorable cliff like conduction band alignment at the buffer/CuInS$_2$ interface could be a possible cause of enhanced interface recombination in the device. In this work, we exploit dir… ▽ More Copper indium disulfide (CuInS$_2$) grown under Cu-rich conditions exhibits high optical quality but suffers predominantly from charge carrier interface recombination resulting in poor solar cell performance. An unfavorable cliff like conduction band alignment at the buffer/CuInS$_2$ interface could be a possible cause of enhanced interface recombination in the device. In this work, we exploit direct and inverse photoelectron spectroscopy together with electrical characterization to investigate the cause of interface recombination in Zn(O,S)/CuInS$_2$ devices. Temperature-dependent current-voltage analysis indeed reveal an activation energy of the dominant charge carrier recombination path, considerably smaller than the absorber bandgap, confirming the dominant recombination channel to be present at the Zn(O,S)/CuInS2 interface. However, photoelectron spectroscopy measurements indicate a small spike like conduction band offset of 0.1 eV at the Zn(O,S) CuInS$_2$ interface, excluding an unfavorable energy level alignment to be the prominent cause for strong interface recombination. The observed band bending upon interface formation also rules out Fermi level pinning as the main reason, leaving near-interface defects (as recently observed in Cu-rich CuInSe2)1 as the likely reason for the performance limiting interface recombination. △ Less

Submitted 13 October, 2021; originally announced October 2021.

$α$-Hypergeometric Uncertain Volatility Models and their Connection to 2BSDEs

Authors: Zaineb Mezdoud, Carsten Hartmann, Mohamed Riad Remita, Omar Kebiri

Abstract: In this article we propose a $α$-hypergeometric model with uncertain volatility (UV) where we derive a worst-case scenario for option pricing. The approach is based on the connexion between a certain class of nonlinear partial differential equations of HJB-type (G-HJB equations), that govern the nonlinear expectation of the UV model and that provide an alternative to the difficult model calibratio… ▽ More In this article we propose a $α$-hypergeometric model with uncertain volatility (UV) where we derive a worst-case scenario for option pricing. The approach is based on the connexion between a certain class of nonlinear partial differential equations of HJB-type (G-HJB equations), that govern the nonlinear expectation of the UV model and that provide an alternative to the difficult model calibration problem of UV models, and second-order backward stochastic differential equations (2BSDEs). Using asymptotic analysis for the G-HJB equation and the equivalent 2BSDE representation, we derive a limit model that provides an accurate description of the worst-case price scenario in cases when the bounds of the UV model are slowly varying. The analytical results are tested by numerical simulations using a deep learning based approximation of the underlying 2BSDE. △ Less

Submitted 16 August, 2021; originally announced August 2021.

Comments: 15 pages, 1 figure

MSC Class: 91G30; 35Q93

Choice of Damping Coefficient in Langevin Dynamics

Authors: Robert D. Skeel, Carsten Hartmann

Abstract: This article considers the application of Langevin dynamics to sampling and investigates how to choose the damping parameter in Langevin dynamics for the purpose of maximizing thoroughness of sampling. Also, it considers the computation of measures of sampling thoroughness. This article considers the application of Langevin dynamics to sampling and investigates how to choose the damping parameter in Langevin dynamics for the purpose of maximizing thoroughness of sampling. Also, it considers the computation of measures of sampling thoroughness. △ Less

Submitted 22 June, 2021; originally announced June 2021.

Comments: 21 pages, 6 figures

MSC Class: 82M37; 65C05

Accurate and Efficient Time Series Matching by Season- and Trend-aware Symbolic Approximation -- Extended Version Including Additional Evaluation and Proofs

Authors: Lars Kegel, Claudio Hartmann, Maik Thiele, Wolfgang Lehner

Abstract: Processing and analyzing time series data\-sets have become a central issue in many domains requiring data management systems to support time series as a native data type. A crucial prerequisite of these systems is time series matching, which still is a challenging problem. A time series is a high-dimensional data type, its representation is storage-, and its comparison is time-consuming. Among th… ▽ More Processing and analyzing time series data\-sets have become a central issue in many domains requiring data management systems to support time series as a native data type. A crucial prerequisite of these systems is time series matching, which still is a challenging problem. A time series is a high-dimensional data type, its representation is storage-, and its comparison is time-consuming. Among the representation techniques that tackle these challenges, the symbolic aggregate approximation (SAX) is the current state of the art. This technique reduces a time series to a low-dimensional space by segmenting it and discretizing each segment into a small symbolic alphabet. However, SAX ignores the deterministic behavior of time series such as cyclical repeating patterns or trend component affecting all segments and leading to a distortion of the symbolic distribution. In this paper, we present a season- and a trend-aware symbolic approximation. We show that this improves the symbolic distribution and increase the representation accuracy without increasing its memory footprint. Most importantly, this enables a more efficient time series matching by providing a match up to three orders of magnitude faster than SAX. △ Less

Submitted 11 October, 2021; v1 submitted 31 May, 2021; originally announced May 2021.

Reachability Analysis of Randomly Perturbed Hamiltonian Systems

Authors: Carsten Hartmann, Lara Neureither, Markus Strehlau

Abstract: In this paper, we revisit energy-based concepts of controllability and reformulate them for control-affine nonlinear systems perturbed by white noise. Specifically, we discuss the relation between controllability of deterministic systems and the corresponding stochastic control systems in the limit of small noise and in the case in which the target state is a measurable subset of the state space.… ▽ More In this paper, we revisit energy-based concepts of controllability and reformulate them for control-affine nonlinear systems perturbed by white noise. Specifically, we discuss the relation between controllability of deterministic systems and the corresponding stochastic control systems in the limit of small noise and in the case in which the target state is a measurable subset of the state space. We derive computable expression for hitting probabilities and mean first hitting times in terms of empirical Gramians, when the dynamics is given by a Hamiltonian system perturbed by dissipation and noise, and provide an easily computable expression for the corresponding controllability function as a function of a subset of the state variables. △ Less

Submitted 1 August, 2021; v1 submitted 5 May, 2021; originally announced May 2021.

Comments: 8 pages, 3 figures

MSC Class: 93E20; 93B03; 82C31

Stochastic gradient descent and fast relaxation to thermodynamic equilibrium: a stochastic control approach

Authors: Tobias Breiten, Carsten Hartmann, Lara Neureither, Upanshu Sharma

Abstract: We study the convergence to equilibrium of an underdamped Langevin equation that is controlled by a linear feedback force. Specifically, we are interested in sampling the possibly multimodal invariant probability distribution of a Langevin system at small noise (or low temperature), for which the dynamics can easily get trapped inside metastable subsets of the phase space. We follow [Chen et al.,… ▽ More We study the convergence to equilibrium of an underdamped Langevin equation that is controlled by a linear feedback force. Specifically, we are interested in sampling the possibly multimodal invariant probability distribution of a Langevin system at small noise (or low temperature), for which the dynamics can easily get trapped inside metastable subsets of the phase space. We follow [Chen et al., J. Math. Phys. 56, 113302, 2015] and consider a Langevin equation that is simulated at a high temperature, with the control playing the role of a friction that balances the additional noise so as to restore the original invariant measure at a lower temperature. We discuss different limits as the temperature ratio goes to infinity and prove convergence to a limit dynamics. It turns out that, depending on whether the lower ("target") or the higher ("simulation") temperature is fixed, the controlled dynamics converges either to the overdamped Langevin equation or to a deterministic gradient flow. This implies that (a) the ergodic limit and the large temperature separation limit do not commute in general, and that (b) it is not possible to accelerate the speed of convergence to the ergodic limit by making the temperature separation larger and larger. We discuss the implications of these observation from the perspective of stochastic optimisation algorithms and enhanced sampling schemes in molecular dynamics. △ Less

Submitted 8 March, 2021; originally announced March 2021.

Comments: 23 pages, 5 figures

MSC Class: 82C31 (Primary); 82M31; 37M25 (Secondary)

Nonasymptotic bounds for suboptimal importance sampling

Authors: Carsten Hartmann, Lorenz Richter

Abstract: Importance sampling is a popular variance reduction method for Monte Carlo estimation, where a notorious question is how to design good proposal distributions. While in most cases optimal (zero-variance) estimators are theoretically possible, in practice only suboptimal proposal distributions are available and it can often be observed numerically that those can reduce statistical performance signi… ▽ More Importance sampling is a popular variance reduction method for Monte Carlo estimation, where a notorious question is how to design good proposal distributions. While in most cases optimal (zero-variance) estimators are theoretically possible, in practice only suboptimal proposal distributions are available and it can often be observed numerically that those can reduce statistical performance significantly, leading to large relative errors and therefore counteracting the original intention. In this article, we provide nonasymptotic lower and upper bounds on the relative error in importance sampling that depend on the deviation of the actual proposal from optimality, and we thus identify potential robustness issues that importance sampling may have, especially in high dimensions. We focus on path sampling problems for diffusion processes, for which generating good proposals comes with additional technical challenges, and we provide numerous numerical examples that support our findings. △ Less

Submitted 18 February, 2021; originally announced February 2021.

Model reduction and uncertainty quantification of multiscale diffusions with parameter uncertainties using nonlinear expectations

Authors: Hafida Bouanani, Carsten Hartmann, Omar Kebiri

Abstract: In this paper we study model reduction of linear and bilinear quadratic stochastic control problems with parameter uncertainties. Specifically, we consider slow-fast systems with unknown diffusion coefficient and study the convergence of the slow process in the limit of infinite scale separation. The aim of our work is two-fold: Firstly, we want to propose a general framework for averaging and hom… ▽ More In this paper we study model reduction of linear and bilinear quadratic stochastic control problems with parameter uncertainties. Specifically, we consider slow-fast systems with unknown diffusion coefficient and study the convergence of the slow process in the limit of infinite scale separation. The aim of our work is two-fold: Firstly, we want to propose a general framework for averaging and homogenisation of multiscale systems with parametric uncertainties in the drift or in the diffusion coefficient. Secondly, we want to use this framework to quantify the uncertainty in the reduced system by deriving a limit equation that represents a worst-case scenario for any given (possibly path-dependent) quantity of interest. We do so by reformulating the slow-fast system as an optimal control problem in which the unknown parameter plays the role of a control variable that can take values in a closed bounded set. For systems with unknown diffusion coefficient, the underlying stochastic control problem admits an interpretation in terms of a stochastic differential equation driven by a G-Brownian motion. We prove convergence of the slow process with respect to the nonlinear expectation on the probability space induced by the G-Brownian motion. The idea here is to formulate the nonlinear dynamic programming equation of the underlying control problem as a forward-backward stochastic differential equation in the G-Brownian motion framework (in brief: G-FBSDE), for which convergence can be proved by standard means. We illustrate the theoretical findings with two simple numerical examples, exploiting the connection between fully nonlinear dynamic programming equations and second-order BSDE (2BSDE): a linear quadratic Gaussian regulator problem and a bilinear multiplicative triad that is a standard benchmark system in turbulence and climate modelling. △ Less

Submitted 9 February, 2021; originally announced February 2021.

Comments: 22 pages, 4 figures

MSC Class: 60G65 (Primary); 70K70; 93E20 (Secondary)

Photochromism and influence of point defect charge states on optical absorption in aluminum nitride (AlN)

Authors: Ivan Gamov, Carsten Hartmann, Thomas Straubinger, Matthias Bickermann

Abstract: In this work, we study the absorption properties of AlN in the range of 1.5 - 5.5 eV, as well as the meta-stable change in absorption induced by ultraviolet (UV) irradiation (photochromism), and the restoration of the initial state under the action of irradiation of 2 - 4 eV or elevated temperatures. UV irradiation results in a decrease of the absorption coefficient from 110 cm$^{-1}$ to 55 cm… ▽ More In this work, we study the absorption properties of AlN in the range of 1.5 - 5.5 eV, as well as the meta-stable change in absorption induced by ultraviolet (UV) irradiation (photochromism), and the restoration of the initial state under the action of irradiation of 2 - 4 eV or elevated temperatures. UV irradiation results in a decrease of the absorption coefficient from 110 cm$^{-1}$ to 55 cm$^{-1}$ at 4.7 eV while in the visible range the absorption coefficient increases from values below 5 cm$^{-1}$ to ~35 cm$^{-1}$. Measurements with two linear polarizations, parallel and perpendicular to the optical axis, provide determination of several different absorption bands at 2.6, 2.8, 3.4, 4.0, 4.5, and 4.8 eV. The bands at 2.6 eV and 3.4 eV identify the defect levels near to the valence band, while the band peaking at 2.8 eV is related to the conduction band. The photochromism allows controlling the absorption of light in two related spectral ranges because the decrease of UV absorption and increase of the visible absorption are related to switching the charge state of the same defects. △ Less

Submitted 17 January, 2021; originally announced January 2021.

Comments: 37 pages, 8 figures, 2 tables

Machine Learning-based Cardinality Estimation in DBMS on Pre-Aggregated Data

Authors: Lucas Woltmann, Claudio Hartmann, Dirk Habich, Wolfgang Lehner

Abstract: Cardinality estimation is a fundamental task in database query processing and optimization. As shown in recent papers, machine learning (ML)-based approaches can deliver more accurate cardinality estimations than traditional approaches. However, a lot of example queries have to be executed during the model training phase to learn a data-dependent ML model leading to a very time-consuming training… ▽ More Cardinality estimation is a fundamental task in database query processing and optimization. As shown in recent papers, machine learning (ML)-based approaches can deliver more accurate cardinality estimations than traditional approaches. However, a lot of example queries have to be executed during the model training phase to learn a data-dependent ML model leading to a very time-consuming training phase. Many of those example queries use the same base data, have the same query structure, and only differ in their predicates. Thus, index structures appear to be an ideal optimization technique at first glance. However, their benefit is limited. To speed up this model training phase, our core idea is to determine a predicate-independent pre-aggregation of the base data and to execute the example queries over this pre-aggregated data. Based on this idea, we present a specific aggregate-enabled training phase for ML-based cardinality estimation approaches in this paper. As we are going to show with different workloads in our evaluation, we are able to achieve an average speedup of 63 with our aggregate-enabled training phase. △ Less

Submitted 19 May, 2020; originally announced May 2020.

Comments: 10 pages, technical report

Journal ref: Datenbank-Spektrum 22 (2022) 1-13

Error bounds for model reduction of feedback-controlled linear stochastic dynamics on Hilbert spaces

Authors: Simon Becker, Carsten Hartmann, Martin Redmann, Lorenz Richter

Abstract: We analyze structure-preserving model order reduction methods for Ornstein-Uhlenbeck processes and linear S(P)DEs with multiplicative noise based on balanced truncation. For the first time, we include in this study the analysis of non-zero initial conditions. We moreover allow for feedback-controlled dynamics for solving stochastic optimal control problems with reduced-order models and prove novel… ▽ More We analyze structure-preserving model order reduction methods for Ornstein-Uhlenbeck processes and linear S(P)DEs with multiplicative noise based on balanced truncation. For the first time, we include in this study the analysis of non-zero initial conditions. We moreover allow for feedback-controlled dynamics for solving stochastic optimal control problems with reduced-order models and prove novel error bounds for a class of linear quadratic regulator problems. We provide numerical evidence for the bounds and discuss the application of our approach to enhanced sampling methods from non-equilibrium statistical mechanics. △ Less

Submitted 17 March, 2022; v1 submitted 12 December, 2019; originally announced December 2019.

Comments: comments welcome

Coarse-graining of non-reversible stochastic differential equations: quantitative results and connections to averaging

Authors: Carsten Hartmann, Lara Neureither, Upanshu Sharma

Abstract: This work is concerned with model reduction of stochastic differential equations and builds on the idea of replacing drift and noise coefficients of preselected relevant, e.g. slow variables by their conditional expectations. We extend recent results by Legoll & Lelièvre [Nonlinearity 23, 2131, 2010] and Duong et al. [Nonlinearity 31, 4517, 2018] on effective reversible dynamics by conditional exp… ▽ More This work is concerned with model reduction of stochastic differential equations and builds on the idea of replacing drift and noise coefficients of preselected relevant, e.g. slow variables by their conditional expectations. We extend recent results by Legoll & Lelièvre [Nonlinearity 23, 2131, 2010] and Duong et al. [Nonlinearity 31, 4517, 2018] on effective reversible dynamics by conditional expectations to the setting of general non-reversible processes with non-constant diffusion coefficient. We prove relative entropy and Wasserstein error estimates for the difference between the time marginals of the effective and original dynamics as well as an entropy error bound for the corresponding path space measures. A comparison with the averaging principle for systems with time-scale separation reveals that, unlike in the reversible setting, the effective dynamics for a non-reversible system need not agree with the averaged equations. We present a thorough comparison for the Ornstein-Uhlenbeck process and make a conjecture about necessary and sufficient conditions for when averaged and effective dynamics agree for nonlinear non-reversible processes. The theoretical results are illustrated with suitable numerical examples. △ Less

Submitted 4 March, 2020; v1 submitted 14 November, 2019; originally announced November 2019.

Variational approach to rare event simulation using least-squares regression

Authors: Carsten Hartmann, Omar Kebiri, Lara Neureither, Lorenz Richter

Abstract: We propose an adaptive importance sampling scheme for the simulation of rare events when the underlying dynamics is given by a diffusion. The scheme is based on a Gibbs variational principle that is used to determine the optimal (i.e. zero-variance) change of measure and exploits the fact that the latter can be rephrased as a stochastic optimal control problem. The control problem can be solved by… ▽ More We propose an adaptive importance sampling scheme for the simulation of rare events when the underlying dynamics is given by a diffusion. The scheme is based on a Gibbs variational principle that is used to determine the optimal (i.e. zero-variance) change of measure and exploits the fact that the latter can be rephrased as a stochastic optimal control problem. The control problem can be solved by a stochastic approximation algorithm, using the Feynman-Kac representation of the associated dynamic programming equations, and we discuss numerical aspects for high-dimensional problems along with simple toy examples. △ Less

Submitted 16 April, 2019; v1 submitted 26 January, 2019; originally announced January 2019.

Comments: 28 pages, 7 figures

MSC Class: 65C05 (primary); 65C30; 92C40 (secondary)

Infinite-dimensional bilinear and stochastic balanced truncation with error bounds

Authors: Simon Becker, Carsten Hartmann

Abstract: Along the ideas of Curtain and Glover, we extend the balanced truncation method for infinite-dimensional linear systems to bilinear and stochastic systems. Specifically , we apply Hilbert space techniques used in many-body quantum mechanics to establish error bounds for the truncated system and prove convergence results. The functional analytic setting allows us to obtain mixed Hardy space error b… ▽ More Along the ideas of Curtain and Glover, we extend the balanced truncation method for infinite-dimensional linear systems to bilinear and stochastic systems. Specifically , we apply Hilbert space techniques used in many-body quantum mechanics to establish error bounds for the truncated system and prove convergence results. The functional analytic setting allows us to obtain mixed Hardy space error bounds for both finite-and infinite-dimensional systems, and it is then applied to the model reduction of stochastic evolution equations driven by Wiener noise. △ Less

Submitted 26 November, 2018; v1 submitted 13 June, 2018; originally announced June 2018.

Comments: 37 pages

Exciton-Dominated Core-Level Absorption Spectra of Hybrid Organic-Inorganic Lead Halide Perovskites

Authors: Christian Vorwerk, Claudia Hartmann, Caterina Cocchi, Golnaz Sadoughi, Severin N. Habisreutinger, Roberto Félix, Regan G. Wilks, Henry J. Snaith, Marcus Bär, Claudia Draxl

Abstract: In a combined theoretical and experimental work, we investigate X-ray Absorption Near-Edge Structure (XANES) spectroscopy of the I $L_3$ and the Pb $M_5$ edges of the methylammonium lead iodide ($\textrm{MAPbI}_3$) hybrid inorganic-organic perovskite and its binary phase $\textrm{PbI}_2$. The absorption onsets are dominated by bound excitons with sizable binding energies of a few hundred meV and p… ▽ More In a combined theoretical and experimental work, we investigate X-ray Absorption Near-Edge Structure (XANES) spectroscopy of the I $L_3$ and the Pb $M_5$ edges of the methylammonium lead iodide ($\textrm{MAPbI}_3$) hybrid inorganic-organic perovskite and its binary phase $\textrm{PbI}_2$. The absorption onsets are dominated by bound excitons with sizable binding energies of a few hundred meV and pronounced anisotropy. The spectra of both materials exhibit remarkable similarities, suggesting that the fingerprints of core excitations in $\textrm{MAPbI}_3$ are essentially given by its inorganic component, with negligible influence from the organic groups. The theoretical analysis complementing experimental observations provides the conceptual insights required for a full characterization of this complex material. △ Less

Submitted 18 April, 2018; originally announced April 2018.

Journal ref: J. Phys. Chem. Lett., 2018, 9, pp 1852-1858

Jarzynski's equality, fluctuation theorems, and variance reduction: Mathematical analysis and numerical algorithms

Authors: Carsten Hartmann, Christof Schuette, Wei Zhang

Abstract: In this paper, we study Jarzynski's equality and fluctuation theorems for diffusion processes. While some of the results considered in the current work are known in the (mainly physics) literature, we review and generalize these nonequilibrium theorems using mathematical arguments, therefore enabling further investigations in the mathematical community. On the numerical side, variance reduction ap… ▽ More In this paper, we study Jarzynski's equality and fluctuation theorems for diffusion processes. While some of the results considered in the current work are known in the (mainly physics) literature, we review and generalize these nonequilibrium theorems using mathematical arguments, therefore enabling further investigations in the mathematical community. On the numerical side, variance reduction approaches such as importance sampling method are studied in order to compute free energy differences based on Jarzynski's equality. △ Less

Submitted 8 April, 2019; v1 submitted 25 March, 2018; originally announced March 2018.

Comments: journal version

MSC Class: 82C05; 60J60

Journal ref: J. Stat. Phys., 2019

Cage occupancies of methane hydrates: Results from synchrotron X-ray diffraction and Raman spectroscopy

Authors: Junfeng Qin, Christiane D. Hartmann, Werner F. Kuhs

Abstract: An accurate knowledge of cage occupancy of methane is central for understanding the physical-chemical properties of gas hydrates, the actual inventory of natural gas in hydrate deposits and the description of gas exchange processes. Here we report the absolute cage occupancies, the cage occupancy ratios and hydration numbers of the synthetic CH4-H2O and CH4-D2O hydrates formed from the ice-gas sys… ▽ More An accurate knowledge of cage occupancy of methane is central for understanding the physical-chemical properties of gas hydrates, the actual inventory of natural gas in hydrate deposits and the description of gas exchange processes. Here we report the absolute cage occupancies, the cage occupancy ratios and hydration numbers of the synthetic CH4-H2O and CH4-D2O hydrates formed from the ice-gas system under different pressures and temperatures. The results were obtained from Rietveld refinement using high-resolution synchrotron X-ray powder diffraction patterns and from Raman spectroscopic measurements. The small-cage occupancies of methane in the deuterated hydrates are found to be slightly higher than in the hydrogenated form, likely due to their different lattice constants. The CH4 occupancy in the small cages agrees fairly well with the predictions of CSMGem at the formation pressure of 3.5 MPa, but with the increasing formation pressure the disagreement grows up to 11 percent. While some deficiency of the prediction model cannot be excluded, the observed discrepancy may well be due to experimental difficulties of reaching true equilibrium at higher pressures. The experimentally determined large-to-small cage occupancy ratios of the synthetic and natural CH4 hydrates formed from the water-gas system are consistently higher than the results of CSMGem calculations. Possible reasons for these discrepancies will be discussed. △ Less

Submitted 6 March, 2018; originally announced March 2018.

Comments: Unpublished contribution to the 8th International Conference on Gas Hydrates (ICGH-8), Beijing, China, 28 July - 1 August 2015 (was only available to the conference participants)

Adaptive importance sampling with forward-backward stochastic differential equations

Authors: Omar Kebiri, Lara Neureither, Carsten Hartmann

Abstract: We describe an adaptive importance sampling algorithm for rare events that is based on a dual stochastic control formulation of a path sampling problem. Specifically, we focus on path functionals that have the form of cumulate generating functions, which appear relevant in the context of, e.g.~molecular dynamics, and we discuss the construction of an optimal (i.e. minimum variance) change of measu… ▽ More We describe an adaptive importance sampling algorithm for rare events that is based on a dual stochastic control formulation of a path sampling problem. Specifically, we focus on path functionals that have the form of cumulate generating functions, which appear relevant in the context of, e.g.~molecular dynamics, and we discuss the construction of an optimal (i.e. minimum variance) change of measure by solving a stochastic control problem. We show that the associated semi-linear dynamic programming equations admit an equivalent formulation as a system of uncoupled forward-backward stochastic differential equations that can be solved efficiently by a least squares Monte Carlo algorithm. We illustrate the approach with a suitable numerical example and discuss the extension of the algorithm to high-dimensional systems. △ Less

Submitted 29 January, 2019; v1 submitted 14 February, 2018; originally announced February 2018.

MSC Class: 65C05 (primary); 65C30; 92C40 (secondary)

Singularly perturbed forward-backward stochastic differential equations: application to the optimal control of bilinear systems

Authors: Omar Kebiri, Lara Neureither, Carsten Hartmann

Abstract: We study linear-quadratic stochastic optimal control problems with bilinear state dependence for which the underlying stochastic differential equation (SDE) consists of slow and fast degrees of freedom. We show that, in the same way in which the underlying dynamics can be well approximated by a reduced order effective dynamics in the time scale limit (using classical homogenziation results), the a… ▽ More We study linear-quadratic stochastic optimal control problems with bilinear state dependence for which the underlying stochastic differential equation (SDE) consists of slow and fast degrees of freedom. We show that, in the same way in which the underlying dynamics can be well approximated by a reduced order effective dynamics in the time scale limit (using classical homogenziation results), the associated optimal expected cost converges in the time scale limit to an effective optimal cost. This entails that we can well approximate the stochastic optimal control for the whole system by the reduced order stochastic optimal control, which is clearly easier to solve because of lower dimensionality. The approach uses an equivalent formulation of the Hamilton-Jacobi-Bellman (HJB) equation, in terms of forward-backward SDEs (FBSDEs). We exploit the efficient solvability of FBSDEs via a least squares Monte Carlo algorithm and show its applicability by a suitable numerical example. △ Less

Submitted 20 March, 2018; v1 submitted 14 February, 2018; originally announced February 2018.

MSC Class: 49M29 (primary); 93C70 (secondary)

Learning Inverse Statics Models Efficiently

Authors: Rania Rayyes, Daniel Kubus, Carsten Hartmann, Jochen Steil

Abstract: Online Goal Babbling and Direction Sampling are recently proposed methods for direct learning of inverse kinematics mappings from scratch even in high-dimensional sensorimotor spaces following the paradigm of "learning while behaving". To learn inverse statics mappings - primarily for gravity compensation - from scratch and without using any closed-loop controller, we modify and enhance the Online… ▽ More Online Goal Babbling and Direction Sampling are recently proposed methods for direct learning of inverse kinematics mappings from scratch even in high-dimensional sensorimotor spaces following the paradigm of "learning while behaving". To learn inverse statics mappings - primarily for gravity compensation - from scratch and without using any closed-loop controller, we modify and enhance the Online Goal Babbling and Direction Sampling schemes. Moreover, we exploit symmetries in the inverse statics mappings to drastically reduce the number of samples required for learning inverse statics models. Results for a 2R planar robot, a 3R simplified human arm, and a 4R humanoid robot arm clearly demonstrate that their inverse statics mappings can be learned successfully with our modified online Goal Babbling scheme. Furthermore, we show that the number of samples required for the 2R and 3R arms can be reduced by a factor of at least 8 and 16 resp. -depending on the number of discovered symmetries. △ Less

Submitted 17 October, 2017; originally announced October 2017.

WavePacket: A Matlab package for numerical quantum dynamics. II: Open quantum systems, optimal control, and model reduction

Authors: Burkhard Schmidt, Carsten Hartmann

Abstract: WavePacket is an open-source program package for numeric simulations in quantum dynamics. It can solve time-independent or time-dependent linear Schrödinger and Liouville-von Neumann-equations in one or more dimensions. Also coupled equations can be treated, which allows, e.g., to simulate molecular quantum dynamics beyond the Born-Oppenheimer approximation. Optionally accounting for the interacti… ▽ More WavePacket is an open-source program package for numeric simulations in quantum dynamics. It can solve time-independent or time-dependent linear Schrödinger and Liouville-von Neumann-equations in one or more dimensions. Also coupled equations can be treated, which allows, e.g., to simulate molecular quantum dynamics beyond the Born-Oppenheimer approximation. Optionally accounting for the interaction with external electric fields within the semi-classical dipole approximation, WavePacket can be used to simulate experiments involving tailored light pulses in photo-induced physics or chemistry. Being highly versatile and offering visualization of quantum dynamics 'on the fly', WavePacket is well suited for teaching or research projects in atomic, molecular and optical physics as well as in physical or theoretical chemistry. Building on the previous Part I which dealt with closed quantum systems and discrete variable representations, the present Part II focuses on the dynamics of open quantum systems, with Lindblad operators modeling dissipation and dephasing. This part also describes the WavePacket function for optimal control of quantum dynamics, building on rapid monotonically convergent iteration methods. Furthermore, two different approaches to dimension reduction implemented in WavePacket are documented here. In the first one, a balancing transformation based on the concepts of controllability and observability Gramians is used to identify states that are neither well controllable nor well observable. Those states are either truncated or averaged out. In the other approach, the H2-error for a given reduced dimensionality is minimized by H2 optimal model reduction techniques, utilizing a bilinear iterative rational Krylov algorithm. △ Less

Submitted 15 January, 2018; v1 submitted 19 July, 2017; originally announced July 2017.

Journal ref: Comp. Phys. Comm. 228, 229-244 (2018)

Model reduction of controlled Fokker--Planck and Liouville-von Neumann equations

Authors: Peter Benner, Tobias Breiten, Carsten Hartmann, Burkhard Schmidt

Abstract: Model reduction methods for bilinear control systems are compared by means of practical examples of Liouville-von Neumann and Fokker--Planck type. Methods based on balancing generalized system Gramians and on minimizing an H2-type cost functional are considered. The focus is on the numerical implementation and a thorough comparison of the methods. Structure and stability preservation are investiga… ▽ More Model reduction methods for bilinear control systems are compared by means of practical examples of Liouville-von Neumann and Fokker--Planck type. Methods based on balancing generalized system Gramians and on minimizing an H2-type cost functional are considered. The focus is on the numerical implementation and a thorough comparison of the methods. Structure and stability preservation are investigated, and the competitiveness of the approaches is shown for practically relevant, large-scale examples. △ Less

Submitted 26 June, 2017; originally announced June 2017.

Journal ref: J. Comp. Dyn. 7 (1), 1-33 (2020)

Partitioning a macroscopic system into independent subsystems

Authors: Luigi Delle Site, Giovanni Ciccotti, Carsten Hartmann

Abstract: We discuss the problem of partitioning a macroscopic system into a collection of independent subsystems. The partitioning of a system into replica-like subsystems is nowadays a subject of major interest in several field of theoretical and applied physics, and the thermodynamic approach currently favoured by practitioners is based on a phenomenological definition of an interface energy associated w… ▽ More We discuss the problem of partitioning a macroscopic system into a collection of independent subsystems. The partitioning of a system into replica-like subsystems is nowadays a subject of major interest in several field of theoretical and applied physics, and the thermodynamic approach currently favoured by practitioners is based on a phenomenological definition of an interface energy associated with the partition, due to a lack of easily computable expressions for a microscopic (i.e.~particle-based) interface energy. In this article, we outline a general approach to derive sharp and computable bounds for the interface free energy in terms of microscopic statistical quantities. We discuss potential applications in nanothermodynamics and outline possible future directions. △ Less

Submitted 6 June, 2017; v1 submitted 31 March, 2017; originally announced March 2017.

Comments: This is an author-created, un-copyedited version of an article accepted for publication in JSTAT

MSC Class: 37M25; 82B05; 82B24; 94A17

Girsanov reweighting for path ensembles and Markov state models

Authors: Luca Donati, Carsten Hartmann, Bettina G. Keller

Abstract: The sensitivity of molecular dynamics on changes in the potential energy function plays an important role in understanding the dynamics and function of complex molecules.We present a method to obtain path ensemble averages of a perturbed dynamics from a set of paths generated by a reference dynamics. It is based on the concept of path probability measure and the Girsanov theorem, a result from sto… ▽ More The sensitivity of molecular dynamics on changes in the potential energy function plays an important role in understanding the dynamics and function of complex molecules.We present a method to obtain path ensemble averages of a perturbed dynamics from a set of paths generated by a reference dynamics. It is based on the concept of path probability measure and the Girsanov theorem, a result from stochastic analysis to estimate a change of measure of a path ensemble. Since Markov state models (MSM) of the molecular dynamics can be formulated as a combined phase-space and path ensemble average, the method can be extended toreweight MSMs by combining it with a reweighting of the Boltzmann distribution. We demonstrate how to efficiently implement the Girsanov reweighting in a molecular dynamics simulation program by calculating parts of the reweighting factor "on the fly" during the simulation, and we benchmark the method on test systems ranging from a two-dimensional diffusion process to an artificial many-body system and alanine dipeptide and valine dipeptide in implicit and explicit water. The method can be used to study the sensitivity of molecular dynamics on external perturbations as well as to reweight trajectories generated by enhanced sampling schemes to the original dynamics. △ Less

Submitted 16 March, 2017; originally announced March 2017.

Magnetic anisotropy in surface-supported single-ion lanthanide complexes

Authors: Paul Stoll, Matthias Bernien, Daniela Rolf, Fabian Nickel, Qingyu Xu, Claudia Hartmann, Tobias R. Umbach, Jens Kopprasch, Janina N. Ladenthin, Enrico Schierle, Eugen Weschke, Constantin Czekelius, Wolfgang Kuch, Katharina J. Franke

Abstract: Single-ion lanthanide-organic complexes can provide stable magnetic moments with well-defined orientation for spintronic applications on the atomic level. Here, we show by a combined experimental approach of scanning tunneling microscopy and X-ray absorption spectroscopy that dysprosium-tris(1,1,1-trifluoro-4-(2-thienyl)-2,4butanedionate) (Dy(tta)$_3$) complexes deposited on a Au(111) surface unde… ▽ More Single-ion lanthanide-organic complexes can provide stable magnetic moments with well-defined orientation for spintronic applications on the atomic level. Here, we show by a combined experimental approach of scanning tunneling microscopy and X-ray absorption spectroscopy that dysprosium-tris(1,1,1-trifluoro-4-(2-thienyl)-2,4butanedionate) (Dy(tta)$_3$) complexes deposited on a Au(111) surface undergo a molecular distortion, resulting in distinct crystal field symmetry imposed on the Dy ion. This leads to an easy-axis magnetization direction in the ligand plane. Furthermore, we show that tunneling electrons hardly couple to the spin excitations, which we ascribe to the shielded nature of the $4f$ electrons. △ Less

Submitted 30 November, 2016; originally announced November 2016.

Journal ref: Phys. Rev. B 94, 224426 (2016)

The $Z$ decay width in the SMEFT: $y_t$ and $λ$ corrections at one loop

Authors: Christine Hartmann, William Shepherd, Michael Trott

Abstract: We calculate one loop $y_t$ and $λ$ dependent corrections to $\barΓ_Z,\bar{R}_f^0$ and the partial $Z$ widths due to dimension six operators in the Standard Model Effective Field Theory (SMEFT), including finite terms. We assume $\rm CP$ symmetry and a $\rm U(3)^5$ symmetry in the UV matching onto the dimension six operators, dominantly broken by the Standard Model Yukawa matrices. Corrections to… ▽ More We calculate one loop $y_t$ and $λ$ dependent corrections to $\barΓ_Z,\bar{R}_f^0$ and the partial $Z$ widths due to dimension six operators in the Standard Model Effective Field Theory (SMEFT), including finite terms. We assume $\rm CP$ symmetry and a $\rm U(3)^5$ symmetry in the UV matching onto the dimension six operators, dominantly broken by the Standard Model Yukawa matrices. Corrections to these observables are predicted using the input parameters $\{\hatα_{ew}, \hat{M}_Z, \hat{G}_F, \hat{m}_t, \hat{m}_h\}$ extracted with one loop corrections in the same limit. We show that at one loop the number of SMEFT parameters contributing to the precise LEPI pseudo-observables exceeds the number of measurements. As a result the SMEFT parameters contributing to LEP data are formally unbounded when the size of loop corrections are reached until other data is considered in a global analysis. The size of these loop effects is generically a correction of order $\sim\%$ to leading effects in the SMEFT, but we find multiple large numerical coefficients in our calculation at this order. We use a $\rm \overline{MS}$ scheme, modified for the SMEFT, for renormalization. Some subtleties involving novel evanescent scheme dependence present in this result are explained. △ Less

Submitted 11 July, 2017; v1 submitted 29 November, 2016; originally announced November 2016.

Comments: 38 pages + refs, 11 figures, v3: Acknowledgement updated with grant numbers

Optimal control of Markov jump processes : Asymptotic analysis, algorithms and applications to the modelling of chemical reaction systems

Authors: Wei Zhang, Carsten Hartmann, Max von Kleist

Abstract: Markov jump processes are widely used to model natural and engineered processes. In the context of biological or chemical applications one typically refers to the chemical master equation (CME), which models the evolution of the probability mass of any copy-number combination of the interacting particles. When many interacting particles ("species") are considered, the complexity of the CME quickly… ▽ More Markov jump processes are widely used to model natural and engineered processes. In the context of biological or chemical applications one typically refers to the chemical master equation (CME), which models the evolution of the probability mass of any copy-number combination of the interacting particles. When many interacting particles ("species") are considered, the complexity of the CME quickly increases, making direct numerical simulations impossible. This is even more problematic when one aims at controlling the Markov jump processes defined by the CME. In this work, we study both open loop and feedback optimal control problems of the Markov jump processes in the case that the controls can only be switched at fixed control stages. Based on Kurtz's limit theorems, we prove the convergence of the respective control value functions of the underlying Markov decision problem as the copy numbers of the species go to infinity. In the case of the optimal control problem on a finite time-horizon, we propose a hybrid control policy algorithm to overcome the difficulties due to the curse of dimensionality when the copy number of the involved species is large. Two numerical examples demonstrate the suitability of both the analysis and the proposed algorithms. △ Less

Submitted 4 July, 2017; v1 submitted 1 December, 2015; originally announced December 2015.

Comments: 39 pages

MSC Class: 60J27 (Primary); 93E20 (Secondary)

The structure and cage filling of gas hydrates as established by synchrotron powder diffraction data

Authors: Christiane D. Hartmann, Susanne Hemes, Andrzej Falenty, Werner F. Kuhs

Abstract: Synchrotron powder diffraction data were measured on various gas hydrates prepared at different fugacities and recovered to low temperatures. The data extend to sin θ/λ values of at least 1.01 Å-1, in some cases up to 1.48 Å-1, i.e. at least 3 times further in reciprocal space than the best single crystal data obtained so far. Here the results of CO2-hydrate are presented in detail. Structural inf… ▽ More Synchrotron powder diffraction data were measured on various gas hydrates prepared at different fugacities and recovered to low temperatures. The data extend to sin θ/λ values of at least 1.01 Å-1, in some cases up to 1.48 Å-1, i.e. at least 3 times further in reciprocal space than the best single crystal data obtained so far. Here the results of CO2-hydrate are presented in detail. Structural information is obtained from Rietveld refinements to unprecedented precision. All H-bonded O-O distances are found to be close to 2.75-2.76 Å, i.e. very close to the values in ice Ih; differences between a hydrogenated and deuterated host lattice are insignificant. Cage occupancies were obtained for the first time in free refinements together with positional parameters of the guest molecule and their atomic displacement parameters. There is good agreement between the experimental cage occupancies and prediction results from CSMGem. △ Less

Submitted 12 November, 2015; originally announced November 2015.

Comments: Unpublished contribution to the 7th International Conference on Gas Hydrates (ICGH-7), Edinburgh, UK, 17-21 July 2011 (was only available to the conference participants)

Higgs decay to two photons at one-loop in the SMEFT

Authors: Christine Hartmann, Michael Trott

Abstract: We present the calculation of the $\rm CP$ conserving contributions to $Γ(h \rightarrow γγ)$, from dimension six operators at one-loop order, in the linear Standard Model Effective Field Theory. We discuss the impact of these corrections on interpreting current and future experimental bounds on this decay. We present the calculation of the $\rm CP$ conserving contributions to $Γ(h \rightarrow γγ)$, from dimension six operators at one-loop order, in the linear Standard Model Effective Field Theory. We discuss the impact of these corrections on interpreting current and future experimental bounds on this decay. △ Less

Submitted 7 May, 2018; v1 submitted 13 July, 2015; originally announced July 2015.

Comments: 6pp, 1 fig. V4: Typos corrected and comment added

Journal ref: Phys. Rev. Lett. 115, 191801 (2015)

On one-loop corrections in the standard model effective field theory; the $Γ(h \rightarrow γ\, γ)$ case

Authors: Christine Hartmann, Michael Trott

Abstract: We calculate one loop contributions to $Γ(h \rightarrow γ\, γ)$ from higher dimensional operators, in the Standard Model Effective Field Theory (SMEFT). Some technical challenges related to determining Electroweak one loop "finite terms" are discussed and overcome. Although we restrict our attention to $Γ(h \rightarrow γ\, γ)$, several developments we report have broad implications. Firstly, the r… ▽ More We calculate one loop contributions to $Γ(h \rightarrow γ\, γ)$ from higher dimensional operators, in the Standard Model Effective Field Theory (SMEFT). Some technical challenges related to determining Electroweak one loop "finite terms" are discussed and overcome. Although we restrict our attention to $Γ(h \rightarrow γ\, γ)$, several developments we report have broad implications. Firstly, the running of the vacuum expectation value modifies the $\log(μ)$ dependence of processes in a manner that is not captured in some past SMEFT Renormalization Group (RG) calculations. Secondly, higher dimensional operators can source ghost interactions in $R_ξ$ gauges due to a modified gauge fixing procedure. Lastly, higher dimensional operators can contribute with pure finite terms at one loop in a manner that is not anticipated in a RG analysis. These results cast recent speculation on the nature of one loop corrections in the SMEFT in an entirely new light. △ Less

Submitted 11 July, 2017; v1 submitted 11 May, 2015; originally announced May 2015.

Comments: 25 pages, 5 figures, v3: Acknowledgement updated with grant numbers

Importance sampling in path space for diffusion processes with slow-fast variables

Authors: Carsten Hartmann, Christof Schütte, Marcus Weber, Wei Zhang

Abstract: Importance sampling is a widely used technique to reduce the variance of a Monte Carlo estimator by an appropriate change of measure. In this work, we study importance sam- pling in the framework of diffusion process and consider the change of measure which is realized by adding a control force to the original dynamics. For certain exponential type expectation, the corresponding control force of t… ▽ More Importance sampling is a widely used technique to reduce the variance of a Monte Carlo estimator by an appropriate change of measure. In this work, we study importance sam- pling in the framework of diffusion process and consider the change of measure which is realized by adding a control force to the original dynamics. For certain exponential type expectation, the corresponding control force of the optimal change of measure leads to a zero-variance estimator and is related to the solution of a Hamilton-Jacobi-Bellmann equation. We focus on certain diffu- sions with both slow and fast variables, and the main result is that we obtain an upper bound of the relative error for the importance sampling estimators with control obtained from the limiting dynamics. We demonstrate our approximation strategy with an illustrative numerical example. △ Less

Submitted 12 April, 2016; v1 submitted 27 February, 2015; originally announced February 2015.

Comments: 41 pages

MSC Class: 65C05

Journal ref: Probab. Theory Related Fields, 170(1), pp.177-228, 2018

A Holistic Approach for Modeling and Synthesis of Image Processing Applications for Heterogeneous Computing Architectures

Authors: Christian Hartmann, Anna Yupatova, Marc Reichenbach, Dietmar Fey, Reinhard German

Abstract: Image processing applications are common in every field of our daily life. However, most of them are very complex and contain several tasks with different complexities which result in varying requirements for computing architectures. Nevertheless, a general processing scheme in every image processing application has a similar structure, called image processing pipeline: (1) capturing an image, (2)… ▽ More Image processing applications are common in every field of our daily life. However, most of them are very complex and contain several tasks with different complexities which result in varying requirements for computing architectures. Nevertheless, a general processing scheme in every image processing application has a similar structure, called image processing pipeline: (1) capturing an image, (2) pre-processing using local operators, (3) processing with global operators and (4) post-processing using complex operations. Therefore, application-specialized hardware solutions based on heterogeneous architectures are used for image processing. Unfortunately the development of applications for heterogeneous hardware architectures is challenging due to the distribution of computational tasks among processors and programmable logic units. Nowadays, image processing systems are started from scratch which is time-consuming, error-prone and inflexible. A new methodology for modeling and implementing is needed in order to reduce the development time of heterogenous image processing systems. This paper introduces a new holistic top down approach for image processing systems. Two challenges have to be investigated. First, designers ought to be able to model their complete image processing pipeline on an abstract layer using UML. Second, we want to close the gap between the abstract system and the system architecture. △ Less

Submitted 26 February, 2015; originally announced February 2015.

Comments: Presented at DATE Friday Workshop on Heterogeneous Architectures and Design Methods for Embedded Image Systems (HIS 2015) (arXiv:1502.07241)

Report number: DATEHIS/2015/06

Pseudo generators for under-resolved molecular dynamics

Authors: Andreas Bittracher, Carsten Hartmann, Oliver Junge, Péter Koltai

Abstract: Many features of a molecule which are of physical interest (e.g. molecular conformations, reaction rates) are described in terms of its dynamics in configuration space. This article deals with the projection of molecular dynamics in phase space onto configuration space. Specifically, we study the situation that the phase space dynamics is governed by a stochastic Langevin equation and study its re… ▽ More Many features of a molecule which are of physical interest (e.g. molecular conformations, reaction rates) are described in terms of its dynamics in configuration space. This article deals with the projection of molecular dynamics in phase space onto configuration space. Specifically, we study the situation that the phase space dynamics is governed by a stochastic Langevin equation and study its relation with the configurational Smoluchowski equation in the three different scaling regimes: Firstly, the Smoluchowski equations in non-Cartesian geometries are derived from the overdamped limit of the Langevin equation. Secondly, transfer operator methods are used to describe the metastable behaviour of the system at hand, and an explicit small-time asymptotics is derived on which the Smoluchowski equation turns out to govern the dynamics of the position coordinate (without any assumptions on the damping). By using an adequate reduction technique, these considerations are then extended to one-dimensional reaction coordinates. Thirdly, we sketch three different approaches to approximate the metastable dynamics based on time-local information only. △ Less

Submitted 19 May, 2015; v1 submitted 4 February, 2015; originally announced February 2015.

MSC Class: 37A60; 82C31; 70G10 (Primary); 47D07; 47A55; 58J65 (Secondary)

Molecular Dynamics in a Grand Ensemble: Bergmann-Lebowitz model and Adaptive Resolution Simulation

Authors: Animesh Agarwal, Jinglong Zhu, Carsten Hartmann, Han Wang, Luigi Delle Site

Abstract: This article deals with the molecular dynamics simulation of open systems that can exchange energy and matter with a reservoir; the physics of the reservoir and its interactions with the system are described by the model introduced by Bergmann and Lebowitz.Despite its conceptual appeal, the model did not gain popularity in the field of molecular simulation and, as a consequence, did not play a rol… ▽ More This article deals with the molecular dynamics simulation of open systems that can exchange energy and matter with a reservoir; the physics of the reservoir and its interactions with the system are described by the model introduced by Bergmann and Lebowitz.Despite its conceptual appeal, the model did not gain popularity in the field of molecular simulation and, as a consequence, did not play a role in the development of open system molecular simulation techniques, even though it can provide the conceptual legitimation of simulation techniques that mimic open systems. We shall demonstrate that the model can serve as a tool to devise both numerical procedures and conceptual definitions of physical quantities that cannot be defined in a straightforward way by systems with a fixed number of molecules. In particular, we discuss the utility of the Bergmann-Lebowitz (BL) model for the calculation of equilibrium time correlation functions within the Grand Canonical Adaptive Resolution method (GC-AdResS) and report numerical results for the case of liquid water. △ Less

Submitted 27 May, 2015; v1 submitted 2 December, 2014; originally announced December 2014.

Comments: 31 pages, 6 figures

Journal ref: New J. Phys. 17, 083042 (2015)

Besov regularity of solutions to the p-Poisson equation

Authors: Stephan Dahlke, Lars Diening, Christoph Hartmann, Benjamin Scharf, Markus Weimar

Abstract: In this paper, we study the regularity of solutions to the $p$-Poisson equation for all $1<p<\infty$. In particular, we are interested in smoothness estimates in the adaptivity scale $ B^σ_τ(L_τ(Ω))$, $1/τ= σ/d+1/p$, of Besov spaces. The regularity in this scale determines the order of approximation that can be achieved by adaptive and other nonlinear approximation methods. It turns out that, espe… ▽ More In this paper, we study the regularity of solutions to the $p$-Poisson equation for all $1<p<\infty$. In particular, we are interested in smoothness estimates in the adaptivity scale $ B^σ_τ(L_τ(Ω))$, $1/τ= σ/d+1/p$, of Besov spaces. The regularity in this scale determines the order of approximation that can be achieved by adaptive and other nonlinear approximation methods. It turns out that, especially for solutions to $p$-Poisson equations with homogeneous Dirichlet boundary conditions on bounded polygonal domains, the Besov regularity is significantly higher than the Sobolev regularity which justifies the use of adaptive algorithms. This type of results is obtained by combining local Hölder with global Sobolev estimates. In particular, we prove that intersections of locally weighted Hölder spaces and Sobolev spaces can be continuously embedded into the specific scale of Besov spaces we are interested in. The proof of this embedding result is based on wavelet characterizations of Besov spaces. △ Less

Submitted 19 August, 2014; originally announced August 2014.

Comments: 45 pages, 2 figures

Report number: Preprint Nr. 2014-07, Philipps-Universit\"at Marburg, 2014 MSC Class: 35B35; 35J92; 41A25; 41A46; 46E35; 65M99; 65T60

Optimal control of multiscale systems using reduced-order models

Authors: Wei Zhang, Juan C. Latorre, Grigorios A. Pavliotis, Carsten Hartmann

Abstract: We study optimal control of diffusions with slow and fast variables and address a question raised by practitioners: is it possible to first eliminate the fast variables before solving the optimal control problem and then use the optimal control computed from the reduced-order model to control the original, high-dimensional system? The strategy "first reduce, then optimize"--rather than "first opti… ▽ More We study optimal control of diffusions with slow and fast variables and address a question raised by practitioners: is it possible to first eliminate the fast variables before solving the optimal control problem and then use the optimal control computed from the reduced-order model to control the original, high-dimensional system? The strategy "first reduce, then optimize"--rather than "first optimize, then reduce"--is motivated by the fact that solving optimal control problems for high-dimensional multiscale systems is numerically challenging and often computationally prohibitive. We state sufficient and necessary conditions, under which the "first reduce, then control" strategy can be employed and discuss when it should be avoided. We further give numerical examples that illustrate the "first reduce, then optmize" approach and discuss possible pitfalls. △ Less

Submitted 13 June, 2014; originally announced June 2014.

Meshless discretization of LQ-type stochastic control problems

Authors: Ralf Banisch, Carsten Hartmann

Abstract: We propose a novel Galerkin discretization scheme for stochastic optimal control problems on an indefinite time horizon. The control problems are linear-quadratic in the controls, but possibly nonlinear in the state variables, and the discretization is based on the fact that problems of this kind can be transformed into linear boundary value problems by a logarithmic transformation. We show that t… ▽ More We propose a novel Galerkin discretization scheme for stochastic optimal control problems on an indefinite time horizon. The control problems are linear-quadratic in the controls, but possibly nonlinear in the state variables, and the discretization is based on the fact that problems of this kind can be transformed into linear boundary value problems by a logarithmic transformation. We show that the discretized linear problem is dual to a Markov decision problem, the precise form of which depends on the chosen Galerkin basis. We prove a strong error bound in $L^{2}$ for the general scheme and discuss two special cases: a variant of the known Markov chain approximation obtained from a basis of characteristic functions of a box discretization, and a sparse approximation that uses the basis of committor functions of metastable sets of the dynamics; the latter is particularly suited for high-dimensional systems, e.g., control problems in molecular dynamics. We illustrate the method with several numerical examples, one being the optimal control of Alanine dipeptide to its helical conformation. △ Less

Submitted 28 September, 2013; originally announced September 2013.

Comments: 38 pages, 6 figures

MSC Class: 82C31; 65C05; 65C40; 65N30

Grand-canonical-like molecular-dynamics simulations by using an adaptive-resolution technique

Authors: Han Wang, Carsten Hartmann, Christof Schütte, Luigi Delle Site

Abstract: In this work, we provide a detailed theoretical analysis, supported by numerical tests, of the reliability of the adaptive resolution simulation (AdResS) technique in sampling the Grand Canonical ensemble. We demonstrate that the correct density and radial distribution functions in the hybrid region, where molecules change resolution, are two necessary conditions for considering the atomistic and… ▽ More In this work, we provide a detailed theoretical analysis, supported by numerical tests, of the reliability of the adaptive resolution simulation (AdResS) technique in sampling the Grand Canonical ensemble. We demonstrate that the correct density and radial distribution functions in the hybrid region, where molecules change resolution, are two necessary conditions for considering the atomistic and coarse-grained regions in AdResS equivalent to subsystems of a full atomistic system with an accuracy up to the second order with respect to the probability distribution of the system. Moreover, we show that the work done by the thermostat and a thermodynamic force in the transition region is formally equivalent to balance the chemical potential difference between the different resolutions. From these results follows the main conclusion that the atomistic region exchanges molecules with the coarse-grained region in a Grand Canonical fashion with an accuracy up to (at least) second order. Numerical tests, for the relevant case of liquid water at ambient conditions, are carried out to strengthen the conclusions of the theoretical analysis. Finally, in order to show the computational convenience of AdResS as a Grand Canonical set up, we compare our method to the Insertion Particle Method (IMP) in its most efficient computational implementation. This fruitful combination of theoretical principles and numerical evidence candidates the adaptive resolution technique as a natural, general and efficient protocol for Grand Canonical Molecular Dynamics for the case of large systems. △ Less

Submitted 21 January, 2013; originally announced January 2013.

Journal ref: Phys. Rev. X 3, 011018 (2013)